Reconstructing plumage evolution in orioles (Icterus): repeated convergence and reversal in patterns

Several empirical studies suggest that sexually selected characters, including bird plumage, may evolve rapidly and show high levels of convergence and other forms of homoplasy. However, the processes that might generate such convergence have not been explored theoretically. Furthermore, no studies have rigorously addressed this issue using a robust phylogeny and a large number of signal characters. We scored the appearance of 44 adult male plumage characters that varied across New World orioles (Icterus). We mapped the plumage characters onto a molecular phylogeny based on two mitochondrial genes. Reconstructing the evolution of these characters revealed evidence of convergence or reversal in 42 of the 44 plumage characters. No plumage character states are restricted to any groups of species higher than superspecies in the oriole phylogeny. The high frequency of convergence and reversal is reflected in the low overall retention index (RI = 0.66) and the low overall consistency index (CI = 0.28). We found similar results when we mapped plumage changes onto a total evidence tree. Our findings reveal that plumage patterns and colors are highly labile between species of orioles, but highly conserved within the oriole genus. Furthermore, there are at least two overall plumage types that have convergently evolved repeatedly in the three oriole clades. This overall convergence leads to significant conflict between the molecular and plumage data. It is not clear what evolutionary processes lead to this homoplasy in individual characters or convergence in overall pattern. However, evolutionary constraints such as developmental limitations and genetic correlations between characters are likely to play a role. Our results are consistent with the belief that avian plumage and other sexually selected characters may evolve rapidly and may exhibit high homoplasy. The overall convergence in oriole plumage patterns is an interesting evolutionary phenomenon, but it cautions against heavy reliance on plumage characters for constructing phylogenies.     

Recombination confounds interpretations of Wolbachia evolution.

Wolbachia are vertically transmitted bacteria known from arthropods and nematode worms, which are maintained in host populations because they either physiologically benefit infected individuals or parasitically manipulate their reproduction. The different manipulation phenotypes are scattered across the Wolbachia phylogeny, suggesting that there have been multiple evolutions of similar phenotypes. This conclusion relies on the assumption of an absence of recombination between bacterial strains, so that the gene used to reconstruct the phylogeny reflects the evolutionary history of the genes involved in the trait. We tested for recombination by reconstructing the phylogeny of two Wolbachia genes from seven B-subdivision strains. The two genes produced mutually incompatible topologies, indicating that these lineages are subject to genetic recombination. This means that many evolutionary patterns inferred from Wolbachia phylogenies must be re-evaluated. Furthermore, recombination may be an important feature both in the evolution of the manipulation phenotypes and avoidance of Müller's ratchet. Finally, we discuss the implications of recombination for attempts to genetically engineer Wolbachia for use in the control of crop pests and human pathogens.     

A phylogenomic analysis of Escherichia coli / Shigella group: implications of genomic features associated with pathogenicity and ecological adaptation

ABSTRACT: BACKGROUND: The Escherichia coli species contains a variety of commensal and pathogenic strains, and its intraspecific diversity is extraordinarily high. With the availability of an increasing number of E. coli strain genomes, a more comprehensive concept of their evolutionary history and ecological adaptation can be developed using phylogenomic analyses. In this study, we constructed two types of whole-genome phylogenies based on 34 E. coli strains using collinear genomic segments. The first phylogeny was based on the concatenated collinear regions shared by all of the studied genomes, and the second phylogeny was based on the variable collinear regions that are absent from at least one genome. Intuitively, the first phylogeny is likely to reveal the lineal evolutionary history among these strains (i.e., an evolutionary phylogeny), whereas the latter phylogeny is likely to reflect the whole-genome similarities of extant strains (i.e., a similarity phylogeny). RESULTS: Within the evolutionary phylogeny, the strains were clustered in accordance with known phylogenetic groups and phenotypes. When comparing evolutionary and similarity phylogenies, a concept emerges that Shigella may have originated from at least three distinct ancestors and evolved into a single clade. By scrutinizing the properties that are shared amongst Shigella strains but missing in other E. coli genomes, we found that the common regions of the Shigella genomes were mainly influenced by mobile genetic elements, implying that they may have experienced convergent evolution via horizontal gene transfer. Based on an inspection of certain key branches of interest, we identified several collinear regions that may be associated with the pathogenicity of specific strains. Moreover, by examining the annotated genes within these regions, further detailed evidence associated with pathogenicity was revealed. CONCLUSIONS: Collinear regions are reliable genomic features used for phylogenomic analysis among closely related genomes while linking the genomic diversity with phenotypic differences in a meaningful way. The pathogenicity of a strain may be associated with both the arrival of virulence factors and the modification of genomes via mutations. Such phylogenomic studies that compare collinear regions of whole genomes will help to better understand the evolution and adaptation of closely related microbes and E. coli in particular.     

猫科动物的分类——分子系统学框架和化石证据

猫科动物的演化关系传统上一直是富有争议的,这很大程度上造成了猫科在属一级的分类众说纷纭。新兴的分子生物学技术为我们理解猫科动物的系统发育提供了新的视角,本文回顾了近些年来的猫科分子系统学研究,以Johnson等（2006）和李钢等（2016）的研究结果为核心,同时结合确切的化石记录进行佐证,推断了猫科的演化历史。我们发现,分子系统学提出的现代猫科动物于晚中新世开始适应辐射、先后分化出了8个支系等观点,与化石证据基本吻合;但是分子系统学对部分支系起源的推断与化石证据存在差异。结合化石记录的分析,我们认为现代猫科8个支系中除狞猫支系和虎猫支系以外,其余均最有可能起源于亚洲,现代猫科动物在演化过程中至少经历了30次洲际迁徙。结合演化历史和形态学研究,我们遵循邱占祥等（2004）的观点,将现存的猫科动物全部归于猫亚科（Felinae）,并将现代猫亚科动物划分为15个属40个物种。     

Correlated evolution in fig pollination

This is the first comparative study of correlated evolution between figs (Ficus species, Moraceae) and their pollinators (Hymenoptera: Agaoninae) based on molecular phylogenies of both lineages. Fig relationships based on the internal transcribed spacer region (ITS) of nuclear ribosomal DNA and pollinator relationships inferred from mitochondrial cytochrome oxidase I (COI) sequences enabled the study of correlated evolution based on molecular phylogenies for the largest set of interacting species ever compared. Comparative methods have been applied to tests of adaptation, but the application of these methods in tests of coadaptation, defined as reciprocal evolutionary change in interacting lineages, has received less attention. I have extended tests of correlated evolution between two traits along a phylogeny to the case of interacting lineages, where two traits may or may not share a common phylogenetic history. Independent contrasts and phylogenetic autocorrelation rejected the null hypothesis that trait correlations within lineages are stronger than trait correlations between interacting lineages. Fig style lengths and pollinator ovipositor lengths, for example, were more highly correlated than were pollinator body size and ovipositor length. Mutualistic interactions between figs and their pollinators illustrate the novel ways in which phylogenies and comparative methods can detect patterns of correlated evolution. The most outstanding evidence of correlated evolution between these obligate mutualists is that interacting trait correlations are stronger than within-lineage allometric relationships.     

Ecotypes of the Mycobacterium tuberculosis complex

A phylogeny of the Mycobacterium tuberculosis complex has recently shown that the animal-adapted strains are found in a single lineage marked by the deletion of chromosomal region 9 (RD9) [Brosch et al., 2002. A new evolutionary scenario for the Mycobacterium tuberculosis complex. Proc. Natl Acad. Sci. USA 99 (6), 3684-3689]. We have obtained the spoligotype patterns of the RD9 deleted strains used to generate this new evolutionary scenario and we show that the presence of spoligotype spacers 3, 9, 16, 39, and 40-43 is phylogenetically informative in this lineage. We have used the phylogenetically informative spoligotype spacers to screen a database of spoligotype patterns and have identified further members of a group of strains apparently host-adapted to antelopes. The presence of the spoligotype spacers is congruent with the phylogeny generated by chromosomal deletions, suggesting that recombination is rare or absent between strains of this lineage. The phylogenetically informative spacers, in concert with the previously identified single nucleotide mutations and chromosomal deletions, can be used to identify a series of clades in the RD9 deleted lineage each with a separate host preference. Finally, we discuss the application of the ecotype concept to this series of clades and suggest that the M. tuberculosis complex may best be described as a series of host-adapted ecotypes.     

Genomic Diversity and Evolution of Mycobacterium ulcerans Revealed by Next-Generation Sequencing

Mycobacterium ulcerans is the causative agent of Buruli ulcer, the third most common mycobacterial disease after tuberculosis and leprosy. It is an emerging infectious disease that afflicts mainly children and youths in West Africa. Little is known about the evolution and transmission mode of M. ulcerans, partially due to the lack of known genetic polymorphisms among isolates, limiting the application of genetic epidemiology. To systematically profile single nucleotide polymorphisms (SNPs), we sequenced the genomes of three M. ulcerans strains using 454 and Solexa technologies. Comparison with the reference genome of the Ghanaian classical lineage isolate Agy99 revealed 26,564 SNPs in a Japanese strain representing the ancestral lineage. Only 173 SNPs were found when comparing Agy99 with two other Ghanaian isolates, which belong to the two other types previously distinguished in Ghana by variable number tandem repeat typing. We further analyzed a collection of Ghanaian strains using the SNPs discovered. With 68 SNP loci, we were able to differentiate 54 strains into 13 distinct SNP haplotypes. The average SNP nucleotide diversity was low (average 0.06–0.09 across 68 SNP loci), and 96% of the SNP locus pairs were in complete linkage disequilibrium. We estimated that the divergence of the M. ulcerans Ghanaian clade from the Japanese strain occurred 394 to 529 thousand years ago. The Ghanaian subtypes diverged about 1000 to 3000 years ago, or even much more recently, because we found evidence that they evolved significantly faster than average. Our results offer significant insight into the evolution of M. ulcerans and provide a comprehensive report on genetic diversity within a highly clonal M. ulcerans population from a Buruli ulcer endemic region, which can facilitate further epidemiological studies of this pathogen through the development of high-resolution tools.     

Advances in the use of DNA barcodes to build a community phylogeny for tropical trees in a Puerto Rican forest dynamics plot

Background

Species number, functional traits, and phylogenetic history all contribute to characterizing the biological diversity in plant communities. The phylogenetic component of diversity has been particularly difficult to quantify in species-rich tropical tree assemblages. The compilation of previously published (and often incomplete) data on evolutionary relationships of species into a composite phylogeny of the taxa in a forest, through such programs as Phylomatic, has proven useful in building community phylogenies although often of limited resolution. Recently, DNA barcodes have been used to construct a robust community phylogeny for nearly 300 tree species in a forest dynamics plot in Panama using a supermatrix method. In that study sequence data from three barcode loci were used to generate a well-resolved species-level phylogeny.

Methodology/Principal Findings

Here we expand upon this earlier investigation and present results on the use of a phylogenetic constraint tree to generate a community phylogeny for a diverse, tropical forest dynamics plot in Puerto Rico. This enhanced method of phylogenetic reconstruction insures the congruence of the barcode phylogeny with broadly accepted hypotheses on the phylogeny of flowering plants (i.e., APG III) regardless of the number and taxonomic breadth of the taxa sampled. We also compare maximum parsimony versus maximum likelihood estimates of community phylogenetic relationships as well as evaluate the effectiveness of one- versus two- versus three-gene barcodes in resolving community evolutionary history.

Conclusions/Significance

As first demonstrated in the Panamanian forest dynamics plot, the results for the Puerto Rican plot illustrate that highly resolved phylogenies derived from DNA barcode sequence data combined with a constraint tree based on APG III are particularly useful in comparative analysis of phylogenetic diversity and will enhance research on the interface between community ecology and evolution.     

Modeling body size evolution in Felidae under alternative phylogenetic hypotheses

The use of phylogenetic comparative methods in ecological research has advanced during the last twenty years, mainly due to accurate phylogenetic reconstructions based on molecular data and computational and statistical advances. We used phylogenetic correlograms and phylogenetic eigenvector regression (PVR) to model body size evolution in 35 worldwide Felidae (Mammalia, Carnivora) species using two alternative phylogenies and published body size data. The purpose was not to contrast the phylogenetic hypotheses but to evaluate how analyses of body size evolution patterns can be affected by the phylogeny used for comparative analyses (CA). Both phylogenies produced a strong phylogenetic pattern, with closely related species having similar body sizes and the similarity decreasing with increasing distances in time. The PVR explained 65% to 67% of body size variation and all Moran's I values for the PVR residuals were non-significant, indicating that both these models explained phylogenetic structures in trait variation. Even though our results did not suggest that any phylogeny can be used for CA with the same power, or that "good" phylogenies are unnecessary for the correct interpretation of the evolutionary dynamics of ecological, biogeographical, physiological or behavioral patterns, it does suggest that developments in CA can, and indeed should, proceed without waiting for perfect and fully resolved phylogenies.     

Global phylogeny of Mycobacterium tuberculosis based on single nucleotide polymorphism (SNP) analysis: insights into tuberculosis evolution, phylogenetic accuracy of other DNA fingerprinting systems, and recommendations for a minimal standard SNP set

We analyzed a global collection of Mycobacterium tuberculosis strains using 212 single nucleotide polymorphism (SNP) markers. SNP nucleotide diversity was high (average across all SNPs, 0.19), and 96% of the SNP locus pairs were in complete linkage disequilibrium. Cluster analyses identified six deeply branching, phylogenetically distinct SNP cluster groups (SCGs) and five subgroups. The SCGs were strongly associated with the geographical origin of the M. tuberculosis samples and the birthplace of the human hosts. The most ancestral cluster (SCG-1) predominated in patients from the Indian subcontinent, while SCG-1 and another ancestral cluster (SCG-2) predominated in patients from East Asia, suggesting that M. tuberculosis first arose in the Indian subcontinent and spread worldwide through East Asia. Restricted SCG diversity and the prevalence of less ancestral SCGs in indigenous populations in Uganda and Mexico suggested a more recent introduction of M. tuberculosis into these regions. The East African Indian and Beijing spoligotypes were concordant with SCG-1 and SCG-2, respectively; X and Central Asian spoligotypes were also associated with one SCG or subgroup combination. Other clades had less consistent associations with SCGs. Mycobacterial interspersed repetitive unit (MIRU) analysis provided less robust phylogenetic information, and only 6 of the 12 MIRU microsatellite loci were highly differentiated between SCGs as measured by GST. Finally, an algorithm was devised to identify two minimal sets of either 45 or 6 SNPs that could be used in future investigations to enable global collaborations for studies on evolution, strain differentiation, and biological differences of M. tuberculosis.     

Molecular phylogeny and character evolution in terete-stemmed Andean opuntias (Cactaceae-Opuntioideae)

The cacti of tribe Tephrocacteae (Cactaceae-Opuntioideae) are adapted to diverse climatic conditions over a wide area of the southern Andes and adjacent lowlands. They exhibit a range of life forms from geophytes and cushion-plants to dwarf shrubs, shrubs or small trees. To confirm or challenge previous morphology-based classifications and molecular phylogenies, we sampled DNA sequences from the chloroplast trnK/matK region and the nuclear low copy gene phyC and compared the resulting phylogenies with previous data gathered from nuclear ribosomal DNA sequences. The here presented chloroplast and nuclear low copy gene phylogenies were mutually congruent and broadly coincident with the classification based on gross morphology and seed micro-morphology and anatomy. Reconstruction of hypothetical ancestral character states suggested that geophytes and cushion-forming species probably evolved several times from dwarf shrubby precursors. We also traced an increase of embryo size at the expense of the nucellus-derived storage tissue during the evolution of the Tephrocacteae, which is thought to be an evolutionary advantage because nutrients are then more rapidly accessible for the germinating embryo. In contrast to these highly concordant phylogenies, nuclear ribosomal DNA data sampled by a previous study yielded conflicting phylogenetic signals. Secondary structure predictions of ribosomal transcribed spacers suggested that this phylogeny is strongly influenced by the inclusion of paralogous sequence probably arisen by genome duplication during the evolution of this plant group.     

Phylogeny of Mycobacterium tuberculosis Beijing strains constructed from polymorphisms in genes involved in DNA replication, recombination and repair

Background

The Beijing family is a successful group of M. tuberculosis strains, often associated with drug resistance and widely distributed throughout the world. Polymorphic genetic markers have been used to type particular M. tuberculosis strains. We recently identified a group of polymorphic DNA repair replication and recombination (3R) genes. It was shown that evolution of M. tuberculosis complex strains can be studied using 3R SNPs and a high-resolution tool for strain discrimination was developed. Here we investigated the genetic diversity and propose a phylogeny for Beijing strains by analyzing polymorphisms in 3R genes.

Methodology/Principal Findings

A group of 3R genes was sequenced in a collection of Beijing strains from different geographic origins. Sequence analysis and comparison with the ones of non-Beijing strains identified several SNPs. These SNPs were used to type a larger collection of Beijing strains and allowed identification of 26 different sequence types for which a phylogeny was constructed. Phylogenetic relationships established by sequence types were in agreement with evolutionary pathways suggested by other genetic markers, such as Large Sequence Polymorphisms (LSPs). A recent Beijing genotype (Bmyc10), which included 60% of strains from distinct parts of the world, appeared to be predominant.

Conclusions/Significance

We found SNPs in 3R genes associated with the Beijing family, which enabled discrimination of different groups and the proposal of a phylogeny. The Beijing family can be divided into different groups characterized by particular genetic polymorphisms that may reflect pathogenic features. These SNPs are new, potential genetic markers that may contribute to better understand the success of the Beijing family.     

Exon-intron structure and evolution of the Lipocalin gene family

The Lipocalins are an ancient protein family whose expression is currently confirmed in bacteria, protoctists, plants, arthropods, and chordates. The evolution of this protein family has been assessed previously using amino acid sequence phylogenies. In this report we use an independent set of characters derived from the gene structure (exon-intron arrangement) to infer a new lipocalin phylogeny. We also present the novel gene structure of three insect lipocalins. The position and phase of introns are well preserved among lipocalin clades when mapped onto a protein sequence alignment, suggesting the homologous nature of these introns. Because of this homology, we use the intron position and phase of 23 lipocalin genes to reconstruct a phylogeny by maximum parsimony and distance methods. These phylogenies are very similar to the phylogenies derived from protein sequence. This result is confirmed by congruence analysis, and a consensus tree shows the commonalities between the two source trees. Interestingly, the intron arrangement phylogeny shows that metazoan lipocalins have more introns than other eukaryotic lipocalins, and that intron gains have occurred in the C-termini of chordate lipocalins. We also analyze the relationship of intron arrangement and protein tertiary structure, as well as the relationship of lipocalins with members of the proposed structural superfamily of calycins. Our congruence analysis validates the gene structure data as a source of phylogenetic information and helps to further refine our hypothesis on the evolutionary history of lipocalins.     

分子进化研究中系统发生树的重建

在现代分子进化研究中,根据现有生物基因或物种多样性来重建生物的进化史是一个非常重要的问题。一个可靠的系统发生的推断,将揭示出有关生物进化过程的顺序,有助于我们了解生物进化的历史和进化机制。本文简要介绍了系统发生树推断的几个重要问题：建树方法、数据转换、树的可靠性及目前使用较多的几种分析软件。     

Reconstructing labroid evolution with single-copy nuclear DNA.

Fifteen per cent of all living fishes are united in a single suborder (Labroidei) and display a dazzling array of behavioural and ecological traits. The labroids are considered monophyletic and members share a pharyngeal jaw apparatus (PJA) modified for crushing and processing prey. Outside of the explicitly functional PJA, there is no corroborative evidence for a monophyletic Labroidei. Here, we report the first molecular phylogenetic analysis of the suborder. Contrary to morphology-based phylogenies, our single-copy nuclear DNA data do not support labroid families as a natural group. Our data indicate that pharyngognathy has evolved independently among labroid families and that characters of the PJA are not reliable markers of perciform evolution. This work ''crushes'' conventional views of fish phylogeny and should engender novel concepts of piscine life history evolution.     

耐药结核分枝杆菌的适应性代价与补偿性进化

结核病耐药率的攀升是目前全球结核病防控面临的重大挑战。结核分枝杆菌主要通过其基因组中耐药相关基因发生点突变而获得耐药性。由于耐药相关基因通常具有重要的生理功能,其突变往往会导致结核分枝杆菌自身适应性下降,即产生“适应性代价”。然而,耐药结核分枝杆菌可通过进一步积累其他特定突变来回复其适应性,这种能使其适应性上升的突变称为“补偿性突变”。耐药结核分枝杆菌的补偿性进化被认为是耐药结核病广泛传播与流行的生物学基础。近年来,在结核病分子流行病学和基础研究领域,针对耐药结核分枝杆菌的补偿性进化开展了大量研究。本文从结核分枝杆菌的耐药分子机制、耐药突变的适应性代价与补偿性进化,以及补偿性进化如何影响耐药结核病传播等方面,综述耐药结核分枝杆菌补偿性进化的研究进展。     

FACTORS DETERMINING THE ACCURACY OF CLADOGRAM ESTIMATION: EVALUATION USING COMPUTER SIMULATION

We developed a simulation model of phylogenesis with which we generated a large number of phylogenies and associated data matrices. We examined the characteristics of these and evaluated the success of three taxonomic methods (Wagner parsimony, character compatibility, and UPGMA clustering) as estimators of phylogeny, paying particular attention to the consequences of changes in certain evolutionary assumptions: relative rate of evolution in three different evolutionary contexts (phyletic, parent lineage, and daughter lineage); relative rate of evolution in different directions (novel forward, convergent forward, or reverse); variation of evolutionary rates; and topology of the phylogenetic tree. Except for variation of evolutionary rates, all the evolutionary parameters that we controlled had significant effects on accuracy of phylogenetic reconstructions. Unexpectedly, the topology of the phylogeny was the most important single factor affecting accuracy; some phylogenies are more readily estimated than others for simply historical reasons. We conclude that none of the three estimation methods is very accurate, that the differences in accuracy among them are rather small, and that historical effects (the branching pattern of a phylogeny) may outweigh biological effects in determining the accuracy with which a phylogeny can be reconstructed.     

Hominoid cytogenetics and evolution

Much of the literature on the chromosomes of the Hominoidea exists in virtual isolation from both evolutionary theory and physical anthropology. Several unjustified speculations about hominoid affinities in the literature of cytogenetics may be attributed to the effects of this isolation. In this paper, the literature of comparative hominoid cytogenetics is reviewed, and that on chromosomal band patterns and repetitive DNA distributions relative to current evolutionary theory is discussed. These data are critically analyzed and shown to be more consistent with an orthodox hominoid phylogeny than with heterodox phylogenies. Rates and modes of karyotypic evolution are also discussed in an attempt to begin to assimilate the study of hominoid chromosomes within the framework of physical anthropology.     

Assessing phylogenetic dependence of morphological traits using co-inertia prior to investigate character evolution in Loricariinae catfishes

With the increase of laboratory facilities, molecular phylogenies are playing a predominant role in evolutionary analyses. However, understanding the evolution of morphological traits remains essential for a comprehensive view of the evolution of a group. Here we present a new approach based on co-inertia analysis for identifying characters which variations are dependent to the phylogeny, a prerequisite for analyzing the evolution of characters. Our approach has the advantage of treating the full data set at once, including qualitative and quantitative variables. It provides a graphical output giving the contribution of each variable to the co-structure, allowing a direct discrimination among phylogenetically dependent and independent variables. We have implemented this approach in deciphering the evolution of morphological traits in a highly specialized group of Neotropical catfishes: the Loricariinae. We have first inferred a molecular phylogeny of this group based on the 12S and 16S mitochondrial genes. The resulting phylogeny indicated that the subtribe Harttiini was restricted to the single genus Harttia, and within the subtribe Loricariini, two sister subtribes were distinguished, Sturisomina (new subtribe), and Loricariina. Among Loricariina, the morphological groups Loricariichthys and Loricaria+Pseudohemiodon were confirmed. The co-inertia analysis highlighted a strong relationship between the morphological and the genetic data sets, and identified three quantitative and eight qualitative variables linked to the phylogeny. The evolution of quantitative variables was assessed using the orthogram method and showed a major punctual event in the evolution of the number of caudal-fin rays, and a more gradual pattern of evolution of the number of teeth along the phylogeny. The evolution of qualitative variables was inferred using ancestral states reconstructions and highlighted parallel patterns of evolution in characters linked to the mouth, suggesting co-evolution of the traits for adapting to divergent substrates.     

A cophylogenetic perspective of RNA-virus evolution

The extent to which viruses and their hosts codiverge remains an open question, given that numerous cases of both "cospeciation" and horizontal switching have recently been documented. DNA viruses that form persistent infections are thought to be the most likely candidates for phylogenetic congruence. Phylogenetic reconciliation analysis was used to compare established phylogenies for four RNA viruses and their hosts. The analysis employs a cophylogeny mapping technique, implemented in TreeMap v2.0, to find the most parsimonious combinations of evolutionary events able to reconcile any incongruence. This technique is guaranteed to recover all potentially optimal solutions to the reconciled tree and specifically tests the null hypothesis that an associate phylogeny is no more congruent with a host phylogeny than would be a random tree with the same taxon set. Phylogenies for Hantavirus, Spumavirus, and avian sarcoma leukosis virus were found to be significantly similar to their host trees, whereas Lyssavirus and Arenavirus displayed no significant congruence. These results demonstrate that RNA viruses are able to form stable associations with their hosts over evolutionary time scales and that the details of such associations are consistent with persistent infection being a necessary but not sufficient precondition.