Estimation of primate speciation dates using local molecular clocks

Robust time estimation reconciles views of the antiquity of placental mammals

BACKGROUND

Molecular studies have reported divergence times of modern placental orders long before the Cretaceous-Tertiary boundary and far older than paleontological data. However, this discrepancy may not be real, but rather appear because of the violation of implicit assumptions in the estimation procedures, such as non-gradual change of evolutionary rate and failure to correct for convergent evolution.

METHODOLOGY/PRINCIPAL FINDINGS

New procedures for divergence-time estimation robust to abrupt changes in the rate of molecular evolution are described. We used a variant of the multidimensional vector space (MVS) procedure to take account of possible convergent evolution. Numerical simulations of abrupt rate change and convergent evolution showed good performance of the new procedures in contrast to current methods. Application to complete mitochondrial genomes identified marked rate accelerations and decelerations, which are not obtained with current methods. The root of placental mammals is estimated to be approximately 18 million years more recent than when assuming a log Brownian motion model. Correcting the pairwise distances for convergent evolution using MVS lowers the age of the root about another 20 million years compared to using standard maximum likelihood tree branch lengths. These two procedures combined revise the root time of placental mammals from around 122 million years ago to close to 84 million years ago. As a result, the estimated distribution of molecular divergence times is broadly consistent with quantitative analysis of the North American fossil record and traditional morphological views.

CONCLUSIONS/SIGNIFICANCE

By including the dual effects of abrupt rate change and directly accounting for convergent evolution at the molecular level, these estimates provide congruence between the molecular results, paleontological analyses and morphological expectations. The programs developed here are provided along with sample data that reproduce the results of this study and are especially applicable studies using genome-scale sequence lengths.

Chromosomal speciation of humans and chimpanzees revisited: studies of DNA divergence within inverted regions

The human and chimpanzee karyotypes are distinguishable in terms of nine pericentric inversions. According to the recombination suppression model of speciation, these inversions could have promoted the process of parapatric speciation between hominoid populations ancestral to chimpanzees and humans. Were recombination suppression to have occurred in inversion heterozygotes, gene flow would have been reduced, resulting in the accumulation of genetic incompatibilities leading to reproductive isolation and eventual speciation. In an attempt to detect the molecular signature of such events, the sequence divergence of non-coding DNA was compared between humans and chimpanzees. Precise knowledge of the locations of the inversion breakpoints permitted accurate discrimination between inverted and non-inverted regions. Contrary to the predictions of the recombination suppression model, sequence divergence was found to be lower in inverted chromosomal regions as compared to non-inverted regions, albeit with borderline statistical significance. Thus, no signature of recombination suppression resulting from inversion heterozygosity appears to be detectable by analysis of extant human and chimpanzee non-coding DNA. The precise delineation of the inversion breakpoints may nevertheless still prove helpful in identifying potential speciation-relevant genes within the inverted regions.

Understanding the recent evolution of the human genome: insights from human-chimpanzee genome comparisons

The sequencing of the chimpanzee genome and the comparison with its human counterpart have begun to reveal the spectrum of genetic changes that has accompanied human evolution. In addition to gross karyotypic rearrangements such as the fusion that formed human chromosome 2 and the human-specific pericentric inversions of chromosomes 1 and 18, there is considerable submicroscopic structural variation involving deletions, duplications, and inversions. Lineage-specific segmental duplications, detected by array comparative genomic hybridization and direct sequence comparison, have made a very significant contribution to this structural divergence, which is at least three-fold greater than that due to nucleotide substitutions. Since structural genomic changes may have given rise to irreversible functional differences between the diverging species, their detailed analysis could help to identify the biological processes that have accompanied speciation. To this end, interspecies comparisons have revealed numerous human-specific gains and losses of genes as well as changes in gene expression. The very considerable structural diversity (polymorphism) evident within both lineages has, however, hampered the analysis of the structural divergence between the human and chimpanzee genomes. The concomitant evaluation of genetic divergence and diversity at the nucleotide level has nevertheless served to identify many genes that have evolved under positive selection and may thus have been involved in the development of human lineage-specific traits. Genes that display signs of weak negative selection have also been identified and could represent candidate loci for complex genomic disorders. Here, we review recent progress in comparing the human and chimpanzee genomes and discuss how the differences detected have improved our understanding of the evolution of the human genome.

Phylogeography of the Western Lyresnake (Trimorphodon biscutatus): testing aridland biogeographical hypotheses across the Nearctic-Neotropical transition

The Western Lyresnake (Trimorphodon biscutatus) is a widespread, polytypic taxon inhabiting arid regions from the warm deserts of the southwestern United States southward along the Pacific versant of Mexico to the tropical deciduous forests of Mesoamerica. This broadly distributed species provides a unique opportunity to evaluate a priori biogeographical hypotheses spanning two major distinct biogeographical realms (the Nearctic and Neotropical) that are usually treated separately in phylogeographical analyses. I investigated the phylogeography of T. biscutatus using maximum likelihood and Bayesian phylogenetic analysis of mitochondrial DNA (mtDNA) from across this species' range. Phylogenetic analyses recovered five well-supported clades whose boundaries are concordant with existing geographical barriers, a pattern consistent with a model of vicariant allopatric divergence. Assuming a vicariance model, divergence times between mitochondrial lineages were estimated using Bayesian relaxed molecular clock methods calibrated using geological information from putative vicariant events. Divergence time point estimates were bounded by broad confidence intervals, and thus these highly conservative estimates should be considered tentative hypotheses at best. Comparison of mtDNA lineages and taxa traditionally recognized as subspecies based on morphology suggest this taxon is comprised of multiple independent lineages at various stages of divergence, ranging from putative secondary contact and hybridization to sympatry of 'subspecies'.

TOWARD A MOLECULAR PHYLOGENY FOR PEROMYSCUS: EVIDENCE FROM MITOCHONDRIAL CYTOCHROME-b SEQUENCES

One hundred DNA sequences from the mitochondrial cytochrome-b gene of 44 species of deer mice (Peromyscus (sensu stricto), 1 of Habromys, 1 of Isthmomys, 2 of Megadontomys, and the monotypic genera Neotomodon, Osgoodomys, and Podomys were used to develop a molecular phylogeny for Peromyscus. Phylogenetic analyses (maximum parsimony, maximum likelihood, and Bayesian inference) were conducted to evaluate alternative hypotheses concerning taxonomic arrangements (sensu stricto versus sensu lato) of the genus. In all analyses, monophyletic clades were obtained that corresponded to species groups proposed by previous authors; however, relationships among species groups generally were poorly resolved. The concept of the genus Peromyscus based on molecular data differed significantly from the most current taxonomic arrangement. Maximum-likelihood and Bayesian trees depicted strong support for a clade placing Habromys, Megadontomys, Neotomodon, Osgoodomys, and Podomys within Peromyscus. If Habromys, Megadontomys, Neotomodon, Osgoodomys, and Podomys are regarded as genera, then several species groups within Peromyscus (sensu stricto) should be elevated to generic rank. Isthmomys was associated with the genus Reithrodontomys; in turn this clade was sister to Baiomys, indicating a distant relationship of Isthmomys to Peromyscus. A formal taxonomic revision awaits synthesis of additional sequence data from nuclear markers together with inclusion of available allozymic and karyotypic data.

Dates from the molecular clock: how wrong can we be?

Large discrepancies have been found in dates of evolutionary events obtained using the molecular clock. Twofold differences have been reported between the dates estimated from molecular data and those from the fossil record; furthermore, different molecular methods can give dates that differ 20-fold. New software attempts to incorporate appropriate allowances for this uncertainty into the calculation of the accuracy of date estimates. Here, we propose that these innovations represent welcome progress towards obtaining reliable dates from the molecular clock, but warn that they are currently unproven, given that the causes and pattern of the discrepancies are the subject of ongoing research. This research implies that many previous studies, even some of those using recently developed methods, might have placed too much confidence in their date estimates, and their conclusions might need to be revised.

Structural divergence between the human and chimpanzee genomes

The structural microheterogeneity evident between the human and chimpanzee genomes is quite considerable and includes inversions and duplications as well as deletions, ranging in size from a few base-pairs up to several megabases (Mb). Insertions and deletions have together given rise to at least 150 Mb of genomic DNA sequence that is either present or absent in humans as compared to chimpanzees. Such regions often contain paralogous sequences and members of multigene families thereby ensuring that the human and chimpanzee genomes differ by a significant fraction of their gene content. There is as yet no evidence to suggest that the large chromosomal rearrangements which serve to distinguish the human and chimpanzee karyotypes have influenced either speciation or the evolution of lineage-specific traits. However, the myriad submicroscopic rearrangements in both genomes, particularly those involving copy number variation, are unlikely to represent exclusively neutral changes and hence promise to facilitate the identification of genes that have been important for human-specific evolution.

Ribosomal DNA in the grasshopper Podisma pedestris: escape from concerted evolution

Eukaryote nuclear ribosomal DNA (rDNA) typically exhibits strong concerted evolution: a pattern in which several hundred rDNA sequences within any one species show little or no genetic diversity, whereas the sequences of different species diverge. We report a markedly different pattern in the genome of the grasshopper Podisma pedestris. Single individuals contain several highly divergent ribosomal DNA groups. Analysis of the magnitude of divergence indicates that these groups have coexisted in the Podisma lineage for at least 11 million years. There are two putatively functional groups, each estimated to be at least 4 million years old, and several pseudogene groups, many of which are transcribed. Southern hybridization and real-time PCR experiments show that only one of the putatively functional types occurs at high copy number. However, this group is scarcely amplified under standard PCR conditions, which means that phylogenetic inference on the basis of standard PCR would be severely distorted. The analysis suggests that concerted evolution has been remarkably ineffective in P. pedestris. We propose that this outcome may be related to the species' exceptionally large genome and the associated low rate of deletion per base pair, which may allow pseudogenes to persist.

[Molecular dating in the genomic era]

The comparison of DNA and protein sequences of extant species might be informative for reconstructing the chronology of evolutionary events on Earth. A phylogenetic tree inferred from molecular data directly depicts the evolutionary affinities of species and indirectly allows estimating the age of their origin and diversification. Molecular dating is achieved by assuming the molecular clock hypothesis, i.e., that the rate of change of nucleotide and amino acid sequences is on average constant over geological time. If paleontological calibrations are available, then absolute divergence times of species can be estimated. However, three major difficulties potentially hamper molecular dating : (1) a limited sample of genes and organisms, (2) a limited number of fossil references, and (3) pervasive variations of molecular evolutionary rates among genomes and species. To circumvent these problems, different solutions have been recently proposed. Larger data sets are built with more genes and more species sampled through the mining of an increasing number of genomes. Moreover, independent key fossils are identified to calibrate molecular clocks, and the uncertainty on their age is integrated in subsequent analyses. Finally, models of molecular rate variations are constructed, and incorporated in the so-called relaxed molecular clock approaches. As an illustration of these improvements, we mention that the debated age of the animal (bilaterian metazoans) diversification may have occurred between 642-761 million years ago (Mya), roughly 100 Ma before the Cambrian explosion. Among mammals, the initial diversification of major placental groups may have taken place around 100 Mya, well before the Cretaceous/Tertiary boundary marking the extinction of dinosaurs.

Characterization of the human lineage-specific pericentric inversion that distinguishes human chromosome 1 from the homologous chromosomes of the great apes

The human and chimpanzee genomes are distinguishable in terms of ten gross karyotypic differences including nine pericentric inversions and a chromosomal fusion. Seven of these large pericentric inversions are chimpanzee-specific whereas two of them, involving human chromosomes 1 and 18, were fixed in the human lineage after the divergence of humans and chimpanzees. We have performed detailed molecular and computational characterization of the breakpoint regions of the human-specific inversion of chromosome 1. FISH analysis and sequence comparisons together revealed that the pericentromeric region of HSA 1 contains numerous segmental duplications that display a high degree of sequence similarity between both chromosomal arms. Detailed analysis of these regions has allowed us to refine the p-arm breakpoint region to a 154.2 kb interval at 1p11.2 and the q-arm breakpoint region to a 562.6 kb interval at 1q21.1. Both breakpoint regions contain human-specific segmental duplications arranged in inverted orientation. We therefore propose that the pericentric inversion of HSA 1 was mediated by intra-chromosomal non-homologous recombination between these highly homologous segmental duplications that had themselves arisen only recently in the human lineage by duplicative transposition.

Strikingly variable divergence times inferred across an Amazonian butterfly 'suture zone'

'Suture zones' are areas where hybrid and contact zones of multiple taxa are clustered. Such zones have been regarded as strong evidence for allopatric divergence by proponents of the Pleistocene forest refugia theory, a vicariance hypothesis frequently used to explain diversification in the Amazon basin. A central prediction of the refugia and other vicariance theories is that the taxa should have a common history so that divergence times should be coincident among taxa. A suture zone for Ithomiinae butterflies near Tarapoto, NE Peru, was therefore studied to examine divergence times of taxa in contact across the zone. We sequenced 1619bp of the mitochondrial COI/COII region in 172 individuals of 31 species from across the suture zone. Inferred divergence times differed remarkably, with divergence between some pairs of widespread species (each of which may have two or more subspecies interacting in the zone, as in the genus Melinaea) being considerably less than that between hybridizing subspecies in other genera (for instance in Oleria). Our data therefore strongly refute a simple hypothesis of simultaneous vicariance and suggest that ongoing parapatric or other modes of differentiation in continuous forest may be important in driving diversification in Amazonia.

Primate molecular divergence dates

With genomic data, alignments can be assembled that greatly increase the number of informative sites for analysis of molecular divergence dates. Here, we present an estimate of the molecular divergence dates for all of the major primate groups. These date estimates are based on a Bayesian analysis of approximately 59.8 kbp of genomic data from 13 primates and 6 mammalian outgroups, using a range of paleontologically supported calibration estimates. Results support a Cretaceous last common ancestor of extant primates (approximately 77 mya), an Eocene divergence between platyrrhine and catarrhine primates (approximately 43 mya), an Oligocene origin of apes and Old World monkeys (approximately 31 mya), and an early Miocene (approximately 18 mya) divergence of Asian and African great apes. These dates are examined in the context of other molecular clock studies.

Species diversification patterns in the Polynesian jumping spider genus Havaika Prószyński, 2001 (Araneae, Salticidae)

Hotspot archipelagoes provide exceptional models for the study of the evolutionary process, due to the effects of isolation and topographical diversity in inducing the formation of unique biotic assemblages. In this paper, we examine the evolutionary patterns exhibited by the jumping spider genus Havaika Prószyński, 2001 in the Polynesian islands of the Hawaiian and Marquesas chains. To date, systematic research on Havaika has been seriously limited by the poor taxonomic knowledge on the group, which was based on a handful of specimens that showed continuous variability and lacked clear-cut diagnostic characters. Here, we circumvent this problem by inferring a phylogeny based on DNA sequences of several fragments including both mitochondrial (protein coding cytochrome oxidase I, NAD1 dehydrogenase, ribosomal 16S, and tRNA leu) and nuclear (internal transcribed spacer 2) genes, and a statistical morphological analyses of a large sample of specimens. Results suggest that the Marquesan and Hawaiian Havaika may be the result of independent colonizations. Furthermore, data provide little support for the standard "progression rule" (evolution in the direction of older to younger islands) in Hawaiian Islands. This may be explained by a recent arrival of the group: age estimates of the different lineages suggest that Havaika colonized the Hawaiian Islands after most of the extant islands were already formed. The lack of clear-cut diagnostic characters among species may also be explained by the recent origin of the group since molecular data do not provide any evidence of hybridization among lineages. Quantitative morphological data coupled with the phylogenetic information allow us to reevaluate the current limitation of Havaika taxonomy. Molecular data support the existence of at least four different evolutionary lineages that are further morphologically diagnosable. However, genealogical relationships are better predicted by geographical affinity (i.e. island) than by morphological characters used in the original descriptions of the species. A pattern of size segregation linked to largely overlapping distributions of some of the species hints at a potential involvement of competition in generating morphological diversity. This study contributes to our understanding on the origin and shaping of the biodiversity of oceanic islands and sets the stage for more detailed studies on particular aspects of these previously overlooked spiders.

An empirical examination of the utility of codon-substitution models in phylogeny reconstruction

Models of codon substitution have been commonly used to compare protein-coding DNA sequences and are particularly effective in detecting signals of natural selection acting on the protein. Their utility in reconstructing molecular phylogenies and in dating species divergences has not been explored. Codon models naturally accommodate synonymous and nonsynonymous substitutions, which occur at very different rates and may be informative for recent and ancient divergences, respectively. Thus codon models may be expected to make an efficient use of phylogenetic information in protein-coding DNA sequences. Here we applied codon models to 106 protein-coding genes from eight yeast species to reconstruct phylogenies using the maximum likelihood method, in comparison with nucleotide- and amino acid-based analyses. The results appeared to confirm that expectation. Nucleotide-based analysis, under simplistic substitution models, were efficient in recovering recent divergences whereas amino acid-based analysis performed better at recovering deep divergences. Codon models appeared to combine the advantages of amino acid and nucleotide data and had good performance at recovering both recent and deep divergences. Estimation of relative species divergence times using amino acid and codon models suggested that translation of gene sequences into proteins led to information loss of from 30% for deep nodes to 66% for recent nodes. Although computational burden makes codon models unfeasible for tree search in large data sets, we suggest that they may be useful for comparing candidate trees. Nucleotide models that accommodate the differences in evolutionary dynamics at the three codon positions also performed well, at much less computational cost. We discuss the relationship between a model's fit to data and its utility in phylogeny reconstruction and caution against use of overly complex substitution models.

Relaxed phylogenetics and dating with confidence

In phylogenetics, the unrooted model of phylogeny and the strict molecular clock model are two extremes of a continuum. Despite their dominance in phylogenetic inference, it is evident that both are biologically unrealistic and that the real evolutionary process lies between these two extremes. Fortunately, intermediate models employing relaxed molecular clocks have been described. These models open the gate to a new field of "relaxed phylogenetics." Here we introduce a new approach to performing relaxed phylogenetic analysis. We describe how it can be used to estimate phylogenies and divergence times in the face of uncertainty in evolutionary rates and calibration times. Our approach also provides a means for measuring the clocklikeness of datasets and comparing this measure between different genes and phylogenies. We find no significant rate autocorrelation among branches in three large datasets, suggesting that autocorrelated models are not necessarily suitable for these data. In addition, we place these datasets on the continuum of clocklikeness between a strict molecular clock and the alternative unrooted extreme. Finally, we present analyses of 102 bacterial, 106 yeast, 61 plant, 99 metazoan, and 500 primate alignments. From these we conclude that our method is phylogenetically more accurate and precise than the traditional unrooted model while adding the ability to infer a timescale to evolution.

The pattern and timing of diversification of Philippine endemic rodents: evidence from mitochondrial and nuclear gene sequences

The 22 genera and 64 species of rodents (Muridae: Murinae) distributed in the Philippine Islands provide a unique opportunity to study patterns and processes of diversification in island systems. Over 90% of these rodent species are endemic to the archipelago, but the relative importance of dispersal from the mainland, dispersal within the archipelago, and in situ differentiation as explanations of this diversity remains unclear, as no phylogenetic hypothesis for these species and relevant mainland forms is currently available. Here we report the results of phylogenetic analyses of the endemic Philippine murines and a wide sampling of murine diversity from outside the archipelago, based on the mitochondrial cytochrome b gene and the nuclear-encoded IRBP exon 1. Analysis of our combined gene data set consistently identified five clades comprising endemic Philippine genera, suggesting multiple invasions of the archipelago. Molecular dating analyses using parametric and semiparametric methods suggest that colonization occurred in at least two stages, one ca. 15 Mya, and another 8 to 12 million years later, consistent with the previous recognition of "Old" and "New" endemic rodent faunas. Ancestral area analysis suggests that the Old Endemics invaded landmasses that are now part of the island of Luzon, whereas the three New Endemic clades may have colonized through either Mindanao, Luzon, or both. Further, our results suggest that most of the diversification of Philippine murines took place within the archipelago. Despite heterogeneity between nuclear and mitochondrial genes in most model parameters, combined analysis of the two data sets using both parsimony and likelihood increased phylogenetic resolution; however, the effect of data combination on support for resolved nodes was method dependent. In contrast, our results suggest that combination of mitochondrial and nuclear data to estimate relatively ancient divergence times can severely compromise those estimates, even when specific methods that account for rate heterogeneity among genes are employed. [Biogeography; divergence date estimation; mitochondrial DNA; molecular systematics; Murinae; nuclear exon; Philippines; phylogeny.].

Subfamilies of cpmA, a gene involved in circadian output, have different evolutionary histories in cyanobacteria

The cpmA gene mediates an output signal in the cyanobacterial circadian system. This gene and its homologues are evolutionarily old, and occur in some non-photosynthetic bacteria and archaea as well as in cyanobacteria. The gene has two functional domains that differ drastically in their level of polymorphism: the N-terminal domain is much more variable than the PurE homologous C-terminal domain. The phylogenetic tree of the cpmA homologues features four main clades (C1-C4), two of which (C1 and C3) belong to cyanobacteria. These cyanobacterial clades match respective ones in the previously reported phylogenetic trees of the other genes involved in the circadian system. The phylogenetic analysis suggested that the C3 subfamily, which comprises the genes from the cyanobacteria with the kaiBC-based circadian system, experienced a lateral transfer, probably from evolutionarily old proteobacteria about 1,000 million years ago. The genes of this subfamily have a significantly higher nonsynonymous substitution rate than those of C1 (2.13 x 10(-10) and 1.53 x 10(-10) substitutions per nonsynonymous site per year, respectively). It appears that the functional and selective constraints of the kaiABC-based system have slowed down the rate of sequence evolution compared to the cpmA homologues of the kaiBC-based system. On the other hand, the differences in the mutation rates between the two cyanobacterial clades point to the different functional constraints of the systems with or without kaiA.

The chimpanzee-specific pericentric inversions that distinguish humans and chimpanzees have identical breakpoints in Pan troglodytes and Pan paniscus

Seven of nine pericentric inversions that distinguish human (HSA) and chimpanzee karyotypes are chimpanzee-specific. In this study we investigated whether the two extant chimpanzee species, Pan troglodytes (common chimpanzee) and Pan paniscus (bonobo), share exactly the same pericentric inversions. The methods applied were FISH with breakpoint-spanning BAC/PAC clones and PCR analyses of the breakpoint junction sequences. Our findings for the homologues to HSA 4, 5, 9, 12, 16, and 17 confirm for the first time at the sequence level that these pericentric inversions have identical breakpoints in the common chimpanzee and the bonobo. Therefore, these inversions predate the separation of the two chimpanzee species 0.86-2 Mya. Further, the inversions distinguishing human and chimpanzee karyotypes may be regarded as early acquisitions, such that they are likely to have been present at the time of human/chimpanzee divergence. According to the chromosomal speciation theory the inversions themselves could have promoted human speciation.

Analysis of transitions at two-fold redundant sites in mammalian genomes. Transition redundant approach-to-equilibrium (TREx) distance metrics

Background

The exchange of nucleotides at synonymous sites in a gene encoding a protein is believed to have little impact on the fitness of a host organism. This should be especially true for synonymous transitions, where a pyrimidine nucleotide is replaced by another pyrimidine, or a purine is replaced by another purine. This suggests that transition redundant exchange (TREx) processes at the third position of conserved two-fold codon systems might offer the best approximation for a neutral molecular clock, serving to examine, within coding regions, theories that require neutrality, determine whether transition rate constants differ within genes in a single lineage, and correlate dates of events recorded in genomes with dates in the geological and paleontological records. To date, TREx analysis of the yeast genome has recognized correlated duplications that established a new metabolic strategies in fungi, and supported analyses of functional change in aromatases in pigs. TREx dating has limitations, however. Multiple transitions at synonymous sites may cause equilibration and loss of information. Further, to be useful to correlate events in the genomic record, different genes within a genome must suffer transitions at similar rates.

Results

A formalism to analyze divergence at two fold redundant codon systems is presented. This formalism exploits two-state approach-to-equilibrium kinetics from chemistry. This formalism captures, in a single equation, the possibility of multiple substitutions at individual sites, avoiding any need to "correct" for these. The formalism also connects specific rate constants for transitions to specific approximations in an underlying evolutionary model, including assumptions that transition rate constants are invariant at different sites, in different genes, in different lineages, and at different times. Therefore, the formalism supports analyses that evaluate these approximations.

Transitions at synonymous sites within two-fold redundant coding systems were examined in the mouse, rat, and human genomes. The key metric (f₂), the fraction of those sites that holds the same nucleotide, was measured for putative ortholog pairs. A transition redundant exchange (TREx) distance was calculated from f₂ for these pairs. Pyrimidine-pyrimidine transitions at these sites occur approximately 14% faster than purine-purine transitions in various lineages. Transition rate constants were similar in different genes within the same lineages; within a set of orthologs, the f₂ distribution is only modest overdispersed. No correlation between disparity and overdispersion is observed. In rodents, evidence was found for greater conservation of TREx sites in genes on the X chromosome, accounting for a small part of the overdispersion, however.

Conclusion

The TREx metric is useful to analyze the history of transition rate constants within these mammals over the past 100 million years. The TREx metric estimates the extent to which silent nucleotide substitutions accumulate in different genes, on different chromosomes, with different compositions, in different lineages, and at different times.

Isolation and characterization of a chimpanzee alphaherpesvirus

Although both beta- and gammaherpesviruses indigenous to great-ape species have been isolated, to date all alphaherpesviruses isolated from apes have proven to be human viruses [herpes simplex virus types 1 (HSV1) and 2 (HSV2) or varicella-zoster virus]. If the alphaherpesviruses have co-evolved with their host species, some if not all ape species should harbour their own alphaherpesviruses. Here, the isolation and characterization of an alphaherpesvirus from a chimpanzee (ChHV) are described. Sequencing of a number of genes throughout the ChHV genome indicates that it is collinear with that of HSV. Phylogenetic analyses place ChHV in a clade with HSV1 and HSV2, the alphaherpesviruses of Old World monkeys comprising a separate clade. Analysis of reactivity patterns of HSV2-immune human sera and ChHV-immune chimpanzee sera by competition ELISA support this relationship. Phylogenetic analyses also place ChHV rather than HSV1 as the closest relative of HSV2.

Mitochondrial molecular clocks and the origin of the major Otocephalan clades (Pisces: Teleostei): A new insight

The Otocephala, a clade including ostariophysan and clupeomorph teleosts, represents about a quarter of total fish species diversity, with about 1000 genera and more than 7000 species. A series of recent papers have defended that the origin of this clade and of its major groups may be significantly older than the oldest fossils of each of these groups suggest. Some of these recent papers explicitly defend a Pangean origin for some otocephalan groups such as the Siluriformes or Cypriniformes. To know whether or not the otocephalans as a whole, and particularly the mainly freshwater, cosmopolitan otophysans could have originated before the splitting of the Pangean supercontinent is of extreme importance, since otophysan fishes are among the most useful animal groups for the determination of historical continental relationships. In the present work we examined divergence times for each major otocephalan group by an analysis of complete mtDNA sequences, in order to investigate if these divergence times support the hypotheses advanced in recent studies. The complete mtDNA sequences of nine representative non-otocephalan fish species and of twenty-one representative otocephalan species was compared. The present study is thus, among the studies dealing with molecular divergence times of teleosts, the one in which a greater number of otocephalan species are included. The divergence times obtained support that the major otocephalan groups had a much older origin than the oldest fossil records available for these groups suggest. The origin of the Otocephala is estimated as having occurred about 282 Mya, with the origin of the Otophysi being estimated at about 251 Mya.

Breakpoint analysis of the pericentric inversion distinguishing human chromosome 4 from the homologous chromosome in the chimpanzee (Pan troglodytes)

The study of breakpoints that occurred during primate evolution promises to yield valuable insights into the mechanisms underlying chromosome rearrangements in both evolution and pathology. Karyotypic differences between humans and chimpanzees include nine pericentric inversions, which may have potentiated the parapatric speciation of hominids and chimpanzees 5-6 million years ago. Detailed analysis of the respective chromosomal breakpoints is a prerequisite for any assessment of the genetic consequences of these inversions. The breakpoints of the inversion that distinguishes human chromosome 4 (HSA4) from its chimpanzee counterpart were identified by fluorescence in situ hybridization (FISH) and comparative sequence analysis. These breakpoints, at HSA4p14 and 4q21.3, do not disrupt the protein coding region of a gene, although they occur in regions with an abundance of LINE and LTR-elements. At 30 kb proximal to the breakpoint in 4q21.3, we identified an as yet unannotated gene, C4orf12, that lacks an homologous counterpart in rodents and is expressed at a 33-fold higher level in human fibroblasts as compared to chimpanzee. Seven out of 11 genes that mapped to the breakpoint regions have been previously analyzed using oligonucleotide-microarrays. One of these genes, WDFY3, exhibits a three-fold difference in expression between human and chimpanzee. To investigate whether the genomic architecture might have facilitated the inversion, comparative sequence analysis was used to identify an approximately 5-kb inverted repeat in the breakpoint regions. This inverted repeat is inexact and comprises six subrepeats with 78 to 98% complementarity. (TA)-rich repeats were also noted at the breakpoints. These findings imply that genomic architecture, and specifically high-copy repetitive elements, may have made a significant contribution to hominoid karyotype evolution, predisposing specific genomic regions to rearrangements.

Genome-wide acceleration of protein evolution in flies (Diptera)

BACKGROUND

The rate of molecular evolution varies widely between proteins, both within and among lineages. To what extent is this variation influenced by genome-wide, lineage-specific effects? To answer this question, we assess the rate variation between insect lineages for a large number of orthologous genes.

RESULTS

When compared to the beetle Tribolium castaneum, we find that the stem lineage of flies and mosquitoes (Diptera) has experienced on average a 3-fold increase in the rate of evolution. Pairwise gene comparisons between Drosophila and Tribolium show a high correlation between evolutionary rates of orthologous proteins.

CONCLUSION

Gene specific divergence rates remain roughly constant over long evolutionary times, modulated by genome-wide, lineage-specific effects. Among the insects analysed so far, it appears that the Tribolium genes show the lowest rates of divergence. This has the practical consequence that homology searches for human genes yield significantly better matches in Tribolium than in Drosophila. We therefore suggest that Tribolium is better suited for comparisons between phyla than the widely employed dipterans.

Molecular dating of caprines using ancient DNA sequences of Myotragus balearicus, an extinct endemic Balearic mammal

BACKGROUND

Myotragus balearicus was an endemic bovid from the Balearic Islands (Western Mediterranean) that became extinct around 6,000-4,000 years ago. The Myotragus evolutionary lineage became isolated in the islands most probably at the end of the Messinian crisis, when the desiccation of the Mediterranean ended, in a geological date established at 5.35 Mya. Thus, the sequences of Myotragus could be very valuable for calibrating the mammalian mitochondrial DNA clock and, in particular, the tree of the Caprinae subfamily, to which Myotragus belongs.

RESULTS

We have retrieved the complete mitochondrial cytochrome b gene (1,143 base pairs), plus fragments of the mitochondrial 12S gene and the nuclear 28S rDNA multi-copy gene from a well preserved Myotragus subfossil bone. The best resolved phylogenetic trees, obtained with the cytochrome b gene, placed Myotragus in a position basal to the Ovis group. Using the calibration provided by the isolation of Balearic Islands, we calculated that the initial radiation of caprines can be dated at 6.2 +/- 0.4 Mya. In addition, alpine and southern chamois, considered until recently the same species, split around 1.6 +/- 0.3 Mya, indicating that the two chamois species have been separated much longer than previously thought.

CONCLUSION

Since there are almost no extant endemic mammals in Mediterranean islands, the sequence of the extinct Balearic endemic Myotragus has been crucial for allowing us to use the Messinian crisis calibration point for dating the caprines phylogenetic tree.

Population size and molecular evolution on islands

The nearly neutral theory predicts that the rate and pattern of molecular evolution will be influenced by effective population size (Ne), because in small populations more slightly deleterious mutations are expected to drift to fixation. This important prediction has not been widely empirically tested, largely because of the difficulty of comparing rates of molecular evolution in sufficient numbers of independent lineages which differ only in Ne. Island endemic species provide an ideal test of the effect of Ne on molecular evolution because species restricted to islands frequently have smaller Ne than closely related mainland species, and island endemics have arisen from mainland lineages many times in a wide range of taxa. We collated a dataset of 70 phylogenetically independent comparisons between island and mainland taxa, including vertebrates, invertebrates and plants, from 19 different island groups. The rate of molecular evolution in these lineages was estimated by maximum likelihood using two measures: overall substitution rate and the ratio of non-synonymous to synonymous substitution rates. We show that island lineages have significantly higher ratios of non-synonymous to synonymous substitution rates than mainland lineages, as predicted by the nearly neutral theory, although overall substitution rates do not differ significantly.

The approximately 30-million-year-old ERVPb1 envelope gene is evolutionarily conserved among hominoids and Old World monkeys

Most human endogenous retroviruses (HERVs) are ancient and their genes are rendered nonfunctional by debilitating mutations. One exception is a recently discovered envelope gene located on chromosome 14. This envelope protein was also recently shown to be expressed in various human tissues and to mediate cell-cell fusion ex vivo. In this study, we demonstrate that this locus (designated ERVPb1) is preserved in Old World monkeys and that the reading frame is maintained. This is congruent with the entry of the HERV-P(b) group between 27 and 36 million years ago as suggested by long terminal repeat divergence. Although the coding capacity is generally lost in the HERV-IP supergroup, the analysis of nucleotide substitutions, lack of stop codons, and single-nucleotide polymorephisms strongly indicates a selective advantage of the ERVPb1 envelope genes during primate evolution. The purifying selection and tissue-specific expression of the human ERVPb1 envelope gene provide strong evidence of a beneficial role for the host.

Potential photosynthesis gene recombination between Prochlorococcus and Synechococcus via viral intermediates

Genes (psbA and psbD) encoding for photosynthetically important proteins were recently found in a number of cultured cyanophage genomes. This phenomenon may be a beneficial trait to the viruses or their photosynthetic cyanobacterial hosts, or may represent an untapped pool of genes involved in the formation of the photosynthetic apparatus that are prone to lateral gene transfer. Here we show analyses of psbA genes from uncultured environmental viruses and prophage populations. We observe a statistically significant separation between viral genes and their potential Synechococcus hosts' genes, and statistical analyses under models of codon evolution indicate that the psbA genes of viruses are evolving under levels of purifying selection that are virtually indistinguishable from their hosts. Furthermore, our data also indicate the possible exchange and reshuffling of psbA genes between Synechococcus and Prochlorococcus via phage intermediates. Overall, these observations raise the possibility that marine viruses serve as a potential genetic pool in shaping the evolution of cyanobacterial photosynthesis.

Relaxed molecular clocks for dating historical plant dispersal events

Age estimation from molecular sequences has emerged as a powerful tool for inferring when a plant lineage arrived in a particular area. Knowing the tenure of lineages within a region is key to understanding the evolution of traits, the evolution of biotic interactions, and the evolution of floras. New analytical methods model change in substitution rates along individual branches of a phylogenetic tree by combining molecular data with time constraints, usually from fossils. These "relaxed clock" approaches can be applied to several gene regions that need not all have the same substitution rates, and they can also incorporate multiple simultaneous fossil calibrations. Since 1995, at least 100 plant biogeographic studies have used molecular-clock dating, and about a fifth has used relaxed clocks. Many of these report evidence of long-distance dispersal. Meta-analyses of studies from the same geographic region can identify directional biases because of prevailing wind or water currents and the relative position and size of landmasses.

Independent intrachromosomal recombination events underlie the pericentric inversions of chimpanzee and gorilla chromosomes homologous to human chromosome 16

Analyses of chromosomal rearrangements that have occurred during the evolution of the hominoids can reveal much about the mutational mechanisms underlying primate chromosome evolution. We characterized the breakpoints of the pericentric inversion of chimpanzee chromosome 18 (PTR XVI), which is homologous to human chromosome 16 (HSA 16). A conserved 23-kb inverted repeat composed of satellites, LINE and Alu elements was identified near the breakpoints and could have mediated the inversion by bringing the chromosomal arms into close proximity with each other, thereby facilitating intrachromosomal recombination. The exact positions of the breakpoints may then have been determined by local DNA sequence homologies between the inversion breakpoints, including a 22-base pair direct repeat. The similarly located pericentric inversion of gorilla (GGO) chromosome XVI, was studied by FISH and PCR analysis. The p- and q-arm breakpoints of the inversions in PTR XVI and GGO XVI were found to occur at slightly different locations, consistent with their independent origin. Further, FISH studies of the homologous chromosomal regions in macaque and orangutan revealed that the region represented by HSA BAC RP11-696P19, which spans the inversion breakpoint on HSA 16q11-12, was derived from the ancestral primate chromosome homologous to HSA 1. After the divergence of orangutan from the other great apes approximately 12 million years ago (Mya), a duplication of the corresponding region occurred followed by its interchromosomal transposition to the ancestral chromosome 16q. Thus, the most parsimonious interpretation is that the gorilla and chimpanzee homologs exhibit similar but nonidentical derived pericentric inversions, whereas HSA 16 represents the ancestral form among hominoids.

Fitting background-selection predictions to levels of nucleotide variation and divergence along the human autosomes

The roles of positive directional selection (selective sweeps) and negative selection (background selection) in shaping the genome-wide distribution of genetic variation in humans remain largely unknown. Here, we optimize the parameter values of a model of the removal of deleterious mutations (background selection) to observed levels of human polymorphism, controlling for mutation rate heterogeneity by using interspecific divergence. A point of "best fit" was found between background-selection predictions and estimates of human effective population sizes, with reasonable parameter estimates whose uncertainty was assessed by bootstrapping. The results suggest that the purging of deleterious alleles has had some influence on shaping levels of human variation, although the effects may be subtle over the majority of the human genome. A significant relationship was found between background-selection predictions and measures of skew in the allele frequency distribution. The genome-wide action of selection (positive and/or negative) is required to explain this observation.

Phylogeny and age of diversification of the planitibia species group of the Hawaiian Drosophila

The Hawaiian Drosophila offer a unique opportunity to examine evolutionary questions because of the known ages of the Hawaiian Islands and the large number of species endemic to this archipelago. One of the more well studied groups of Hawaiian Drosophila is the planitibia species group, a long-standing population genetic model system. Here we present a molecular phylogenetic hypothesis of all 17 taxa in the planitibia group based on nucleotide sequences from two mitochondrial (16S and COII) and four nuclear (Adh, Gpdh, Yp1, and Yp2) loci, accounting for over 4kb of sequence per taxon. We use these data to estimate major divergence times within this group. Our results suggest that the basal diversification within this group, calculated at 6.1 +/- 0.47 MY, predates the oldest high island of Kauai. The older diversifications in this group took place on Kauai, with subsequent colonization and speciation events occurring as new islands became available to Drosophila. Understanding of the phylogenetic relationships of this important group will place the existing population genetic work in a macroevolutionary context and stimulate additional work, particularly on those taxa endemic to the Maui Nui complex of islands.

Bayesian Estimation of Species Divergence Times Under a Molecular Clock Using Multiple Fossil Calibrations with Soft Bounds

We implement a Bayesian Markov chain Monte Carlo algorithm for estimating species divergence times that uses heterogeneous data from multiple gene loci and accommodates multiple fossil calibration nodes. A birth-death process with species sampling is used to specify a prior for divergence times, which allows easy assessment of the effects of that prior on posterior time estimates. We propose a new approach for specifying calibration points on the phylogeny, which allows the use of arbitrary and flexible statistical distributions to describe uncertainties in fossil dates. In particular, we use soft bounds, so that the probability that the true divergence time is outside the bounds is small but nonzero. A strict molecular clock is assumed in the current implementation, although this assumption may be relaxed. We apply our new algorithm to two data sets concerning divergences of several primate species, to examine the effects of the substitution model and of the prior for divergence times on Bayesian time estimation. We also conduct computer simulation to examine the differences between soft and hard bounds. We demonstrate that divergence time estimation is intrinsically hampered by uncertainties in fossil calibrations, and the error in Bayesian time estimates will not go to zero with increased amounts of sequence data. Our analyses of both real and simulated data demonstrate potentially large differences between divergence time estimates obtained using soft versus hard bounds and a general superiority of soft bounds. Our main findings are as follows. (1) When the fossils are consistent with each other and with the molecular data, and the posterior time estimates are well within the prior bounds, soft and hard bounds produce similar results. (2) When the fossils are in conflict with each other or with the molecules, soft and hard bounds behave very differently; soft bounds allow sequence data to correct poor calibrations, while poor hard bounds are impossible to overcome by any amount of data. (3) Soft bounds eliminate the need for "safe" but unrealistically high upper bounds, which may bias posterior time estimates. (4) Soft bounds allow more reliable assessment of estimation errors, while hard bounds generate misleadingly high precisions when fossils and molecules are in conflict.

Taxon sampling effects in molecular clock dating: An example from the African Restionaceae

Three commonly used molecular dating methods for correction of variable rates (non-parametric rate smoothing, penalized likelihood, and Bayesian rate correction) as well as the assumption of a global molecular clock were tested for sensitivity to taxon sampling. The test dataset of 6854 basepairs for 300 terminals includes a nearly complete sample of the Restio-clade of the African Restionaceae (272 of the 288 species), as well as 26 outgroup species. Of this, nested subsets of 35, 51, 80, 120, 150, and the full 300 species were used. Molecular dating experiments with these datasets showed that all methods are sensitive to undersampling, but that this effect is more severe in analyses that use more extreme rate smoothing. Additionally, the undersampling effect is positively related to distance from the calibration node. The combined effect of undersampling and distance from the calibration node resulted in up to threefold differences in the age estimation of nodes from the same dataset with the same calibration point. We suggest that the most suitable methods are penalized likelihood and Bayesian when a global clock assumption has been rejected, as these methods are more successful at finding optimal levels of smoothing to correct for rate heterogeneity, and are less sensitive to undersampling.

Fungi evolution revisited: application of the penalized likelihood method to a Bayesian fungal phylogeny provides a new perspective on phylogenetic relationships and divergence dates of Ascomycota groups

The depiction of evolutionary relationships within phylum Ascomycota is still controversial because of unresolved branching orders in the radiation of major taxa. Here we generated a dataset of 166 small subunit (18S) rDNA sequences, representative of all groups of Fungi and used as input in a Bayesian phylogenetic analysis. This phylogeny suggests that Discomycetes are a basal group of filamentous Ascomycetes and probably maintain ancestor characters since their representatives are intermingled among other filamentous fungi. Also, we show that the evolutionary rate heterogeneity within Ascomycota precludes the assumption of a global molecular clock. Accordingly, we used the penalized likelihood method, and for calibration we included a 400 million-year-old Pyrenomycete fossil considering two distinct scenarios found in the literature, one with an estimated date of 1576 Myr for the plant-animal-fungus split and the other with an estimated date of 965 Myr for the animal-fungus split. Our data show that the current classification of the fossil as a Pyrenomycete is not compatible with the second scenario. Estimates under the first scenario are older than dates proposed in previous studies based on small subunit rDNA sequences but support estimates based on multiprotein analysis, suggesting that the radiation of the major Ascomycota groups occurred into the Proterozoic era.

Molecular characterization of the pericentric inversion of chimpanzee chromosome 11 homologous to human chromosome 9

In addition to the fusion of human chromosome 2, nine pericentric inversions are the most conspicuous karyotype differences between humans and chimpanzees. In this study we identified the breakpoint regions of the pericentric inversion of chimpanzee chromosome 11 (PTR 11) homologous to human chromosome 9 (HSA 9). The break in homology between PTR 11p and HSA 9p12 maps to pericentromeric segmental duplications, whereas the breakpoint region orthologous to 9q21.33 is located in intergenic single-copy sequences. Close to the inversion breakpoint in PTR 11q, large blocks of alpha satellites are located, which indicate the presence of the centromere. Since G-banding analysis and the comparative BAC analyses performed in this study imply that the inversion breaks occurred in the region homologous to HSA 9q21.33 and 9p12, but not within the centromere, the structure of PTR 11 cannot be explained by a single pericentric inversion. In addition to this pericentric inversion of PTR 11, further events like centromere repositioning or a second smaller inversion must be assumed to explain the structure of PTR 11 compared with HSA 9.

Evolution of the Class 2 cytokines and receptors, and discovery of new friends and relatives

The sequencing of a wide variety of genomes and their transcripts has allowed researchers to determine how proteins or protein families evolved and how strongly during evolution a protein has been conserved. In this report, we analyze the evolution of the Class 2 ligands and their cognate receptors by analyzing Class 2 ligand and receptor chain gene sequences from a variety of DNA sequence databases. Both the Class 2 cytokines and receptor chains appear to have developed during the evolution of the chordate phyla: distant homologues of type I interferon (IFN) receptors are the only Class 2 cytokine receptors identified in the Ciona genomes, while a wide variety of Class 2 ligands and receptor chains are encoded in the currently available genomes of bony vertebrates (teleost fish, amphibians, reptiles, birds, mammals). Phylogenetic trees of ligands and ligand-binding receptor chains demonstrate that proteins involved in conferring antiviral activity diverged before those involved in adaptive immunity. Genes encoding IFNs and IFN receptors duplicated multiple times during chordate evolution, suggesting that duplication of genes encoding IFN activity conveyed an evolutionary advantage. Altogether, these data support a model whereby the original Class 2 cytokines and receptors evolved and duplicated during the evolution of the chordate innate immune response system; new receptor and ligand duplications evolved into signaling molecules to fulfill communication requirements of a highly specialized and differentiated vertebrate immune system. In addition, the genomic analysis led to the discovery of some new members of this family.

Primate phylogeny, evolutionary rate variations, and divergence times: a contribution from the nuclear gene IRBP

The first third (ca. 1200 bp) of exon 1 of the nuclear gene encoding the interstitial retinoid-binding protein (IRBP) has been sequenced for 12 representative primates belonging to Lemuriformes, Lorisiformes, Tarsiiformes, Platyrrhini, and Catarrhini, and combined with available data (13 other primates, 11 nonprimate placentals, and 2 marsupials). Phylogenetic analyses using maximum likelihood on nucleotides and amino acids robustly support the monophyly of primates, Strepsirrhini, Lemuriformes, Lorisiformes, Anthropoidea, Catarrhini, and Platyrrhini. It is interesting to note that 1) Tarsiidae grouped with Anthropoidea, and the support for this node depends on the molecular characters considered; 2) Cheirogaleidae grouped within Lemuriformes; and 3) Daubentonia was the sister group of all other Lemuriformes. Study of the IRBP evolutionary rate shows a high heterogeneity within placentals and also within primates. Maximum likelihood local molecular clocks were assigned to three clades displaying significantly contrasted evolutionary rates. Paenungulata were shown to evolve 2.5-3 times faster than Perissodactyla and Lemuriformes. Six independent calibration points were used to estimate splitting ages of the main primate clades, and their compatibility was evaluated. Divergence ages were obtained for the following crown groups: 13.8-14.2 MY for Lorisiformes, 26.5-27.2 MY for Lemuroidea, 39.6-40.7 MY for Lemuriformes, 45.4-46.7 MY for Strepsirrhini, and 56.7-58.4 MY for Haplorrhini. The incompatibility between some paleontological and molecular estimates may reflect the incompleteness of the placental fossil record, and/or indicate that the variable IRBP evolutionary rates are not fully accommodated by local molecular clocks.

Incorporating gene-specific variation when inferring and evaluating optimal evolutionary tree topologies from multilocus sequence data

Because of the increase of genomic data, multiple genes are often available for the inference of phylogenetic relationships. The simple approach for combining multiple genes from the same taxon is to concatenate the sequences and then ignore the fact that different positions in the concatenated sequence came from different genes. Here, we discuss two criteria for inferring the optimal tree topology from data sets with multiple genes. These criteria are designed for multigene data sets where gene-specific evolutionary features are too important to ignore. One criterion is conventional and is obtained by taking the sum of log-likelihoods over all genes. The other criterion is obtained by dividing the log-likelihood for a gene by its sequence length and then taking the arithmetic mean over genes of these ratios. A similar strategy could be adopted with parsimony scores. The optimal tree is then declared to be the one for which the sum or the arithmetic mean is maximized. These criteria are justified within a two-stage hierarchical framework. The first level of the hierarchy represents gene-specific evolutionary features, and the second represents site-specific features for given genes. For testing significance of the optimal topology, we suggest a two-stage bootstrap procedure that involves resampling genes and then resampling alignment columns within resampled genes. An advantage of this procedure over concatenation is that it can effectively account for gene-specific evolutionary features. We discuss the applicability of the two-stage bootstrap idea to the Kishino-Hasegawa test and the Shimodaira-Hasegawa test.

Molecular evolution of recombination hotspots and highly recombining pseudoautosomal regions in hominoids

We examined the effects of recombination on the molecular evolution of noncoding regions in pseudoautosomal regions (PARs) and recombination hotspots in hominoids. The PAR-linked regions analyzed had on average longer branch lengths than those of the recombination hotspots. Moreover, contrary to previous observations, we found no correlation between recombination rate and silent site divergence in our data set and little change in the GC content during recent hominoid evolution. This suggests that the current rate of recombination is not a good indicator of the past rates of recombination for these highly recombining regions. Furthermore, human recombination hotspots show increased AT to GC substitutions in the human lineage, while no such pattern is detected for PAR-linked regions. Together, these observations suggest that recombination hotspots in hominoids are transient in the evolutionary time-scale. Interestingly, the 16p13.3 recombination hotspot locus violates a local molecular clock, though the locus appears to be noncoding and should evolve neutrally. We hypothesize that sudden changes in recombination rate have caused the changes in substitution rate at this locus.

The Effect of Paralogous Lineages on the Application of Reconciliation Analysis by Cophylogeny Mapping

Paralogy defines similarity caused by duplication rather than common descent and is well known in the case of paralogous gene copies within a single genome. The term is here extended to paralogous lineages of associates within a single host. The phylogenies of four genera within the Herpesviridae were reconciled with host phylogenies using cophylogenetic mapping. The observed correspondence for each pair of phylogenies was evaluated through randomization of the viral phylogeny and demonstrated to be greater than expected by chance. A simulation study was then carried out to assess the influence of paralogous lineages on the efficacy of reconciliation analysis. Combining viral taxa from different genera that infected common hosts introduced incongruence into the cophylogenies and reduced both the minimum and maximum observed number of codivergence events relative to the initial analysis of orthologous clades. However, at an average sample size this did not alter the fundamental significance of observed correspondence. With smaller sample sizes, the number of orthologous taxa selected at random from the pool of taxa was reduced. False-negative results then increased in proportion from 0.02 to 0.33. These results demonstrated that reconciliation analysis is robust under conditions of paralogy at "normal" sample sizes but is adversely affected by a combination of paralogy and low sample size. Consideration of phylogenies for Papillomavirus, Atadenovirus, and Mastadenovirus suggest that paralogous lineages may be a widespread phenomenon among DNA viruses and that duplication irrespective of host speciation is an important cause of viral diversification.