A "Total Evidence" Analysis of the Phylogenetic Relationships among the Photosynthetic Stramenopiles

Phylogenetic relationships among the nine major autotrophic stramenopile taxa were inferred in a combined analysis of the rbcL, SSU rDNA, partial LSU rRNA, carotenoid, and ultrastructural data sets. The structure of the shortest combined tree is: (Outgroup, ((((Bacillariophyceae, (Pelagophyceae, Dictyochophyceae)),((Phaeophyceae, Xanthophyceae), Raphidophyceae)), Eustigmatophyceae),(Chrysophyceae, Synurophyceae))). The Synurophyceae/Chrysophyceae is the best supported group followed by the Phaeophyceae/Xanthophyceae and the Pelagophyceae/Dictyochophyceae clades. The monophyletic groups composed of Bacillariophyceae/Pelagophyceae/Dictyochophyceae and Phaeophyceae/Xanthophyceae/Raphidophyceae received the lowest Bremer support values. The optimal combined tree suggests that the diatom frustule is derived from the siliceous "skeleton" in Dictyochophyceae, that the reduced flagellar apparatus arose once in the Bacillariophyceae/Dictyochophyceae/Pelagophyceae clade, and that the specific photoreceptor-eyespot apparatus in Chrysophyceae and the Phaeophyceae/Xantophyceae clade originated independently within the autotrophic stramenopiles. Despite conflicts in tree structure between the most-parsimonious combined phylogeny and the optimal tree(s) of each data partition, it cannot be concluded that extensive incongruence exists between the data sets.     

Sensitivity of phylogeny estimation to taxonomic sampling

Recent studies have shown that addition or deletion of taxa from a data matrix can change the estimate of phylogeny. I used 29 data sets from the literature to examine the effect of taxon sampling on phylogeny estimation within data sets. I then used multiple regression to assess the effect of number of taxa, number of characters, homoplasy, strength of support, and tree symmetry on the sensitivity of data sets to taxonomic sampling. Sensitivity to sampling was measured by mapping characters from a matrix of culled taxa onto optimal trees for that reduced matrix and onto the pruned optimal tree for the entire matrix, then comparing the length of the reduced tree to the length of the pruned complete tree. Within-data-set patterns can be described by a second-order equation relating fraction of taxa sampled to sensitivity to sampling. Multiple regression analyses found number of taxa to be a significant predictor of sensitivity to sampling; retention index, number of informative characters, total support index, and tree symmetry were nonsignificant predictors. I derived a predictive regression equation relating fraction of taxa sampled and number of taxa potentially sampled to sensitivity to taxonomic sampling and calculated values for this equation within the bounds of the variables examined. The length difference between the complete tree and a subsampled tree was generally small (average difference of 0-2.9 steps), indicating that subsampling taxa is probably not an important problem for most phylogenetic analyses using up to 20 taxa.     

A Phylogeny of the Gulls (Aves: Larinae) Inferred from Osteological and Integumentary Characters

Abstract — Gulls (Aves: Larinae) are among the best-studied of birds, yet prior attempts to reconstruct gull relationships have met with little success. In the present study I use 117 characters from the skeleton and 64 from the integument to test gull monophyly and estimate gull phylogeny. One shortest tree, requiring 9747 unweighted changes and having a CI of 0.267, wasLarusis polyphyletic. Although the tree is fully resolved, support for many of the inferred clades is poor. In a comparison of osteological and integumentary evidence, I found that incongruence between the osteological and integumentary character sets accounts for only a minority of the total incongruence observed, and suggest that low between-set incongruence may be a consequence of the low signal-to-noise ratio in each set of characters. I also found that osteological evidence is particularly important for determining higher-level structure, whereas integumentary evidence is important for resolving lower-level relationships within the gull group. Finally, I found that integumentary characters are not dramatically more homoplasious than osteological characters, and argue that casual dismissal of integumentary characters as “too labile” is unwarranted.     

Data congruence, paedomorphosis and salamanders

Background

The retention of ancestral juvenile characters by adult stages of descendants is called paedomorphosis. However, this process can mislead phylogenetic analyses based on morphological data, even in combination with molecular data, because the assessment if a character is primary absent or secondary lost is difficult. Thus, the detection of incongruence between morphological and molecular data is necessary to investigate the reliability of simultaneous analyses. Different methods have been proposed to detect data congruence or incongruence. Five of them (PABA, PBS, NDI, LILD, DRI) are used herein to assess incongruence between morphological and molecular data in a case study addressing salamander phylogeny, which comprises several supposedly paedomorphic taxa. Therefore, previously published data sets were compiled herein. Furthermore, two strategies ameliorating effects of paedomorphosis on phylogenetic studies were tested herein using a statistical rigor. Additionally, efficiency of the different methods to assess incongruence was analyzed using this empirical data set. Finally, a test statistic is presented for all these methods except DRI.

Results

The addition of morphological data to molecular data results in both different positions of three of the four paedomorphic taxa and strong incongruence, but treating the morphological data using different strategies ameliorating the negative impact of paedomorphosis revokes these changes and minimizes the conflict. Of these strategies the strategy to just exclude paedomorphic character traits seem to be most beneficial. Of the three molecular partitions analyzed herein the RAG1 partition seems to be the most suitable to resolve deep salamander phylogeny. The rRNA and mtDNA partition are either too conserved or too variable, respectively. Of the different methods to detect incongruence, the NDI and PABA approaches are more conservative in the indication of incongruence than LILD and PBS.

Conclusion

Paedomorphosis induces strong conflicts and can mislead the phylogenetic analyses even in combined analyses. However, different strategies are efficiently minimizing these problems. Though the exploration of different methods to detect incongruence is preferable NDI and PABA are more conservative than the others and NDI is computational less extensive than PABA.     

Sources of character conflict in a clade of water striders (Heteroptera: Gerridae)

Incongruence among trees reconstructed with different data may stem from historical (gene tree‐species tree conflict) or process (character change biases) phenomena. Regardless of the source, incongruent data, as determined with “global” measures of homoplasy, have often been excluded from parsimony analysis of the combined data. Recent studies suggest that these homoplasy measures do not predict the contribution of each character to overall tree structure. Branch support measures identify, on a character to node basis, sources of support and conflict resulting from a simultaneous analysis of the data. We implement these branch support measures to identify sources of character conflict in a clade of water striders consisting of Gerris Fabricius, Aquarius Schellenberg, and Limnoporus Stål species. Separate analyses of morphology, mitochondrial cytochrome oxidase I (COI), large mitochondrial ribosomal subunit (16SrRNA), and elongation factor‐1α (EF‐1α) data resulted in cladograms that varied in resolution and topological concordance. Simultaneous analysis of the data resulted in two trees that were unresolved for one node in a strict consensus. The topology agreed with current classification except for the placements of Aquarius chilensis and the Aquarius remigis species group closer to Gerris than to congeneric species. Branch support measures indicated that support derived from each data set varied among nodes, but COI had an overall negative effect on branch support. However, Spearman rank correlation of partitioned branch support values indicated no negative associations of branch support between any data sets and a positive association between EF‐1α and 16SrRNA. Thus incongruence among data sets was not drastic and the gene‐tree versus species tree phenomenon was not implicated. Biases in character change were a more likely reason for incongruence, although saturation curves and incongruence length difference for COI indicated little potential for homoplasy. However, a posteriori inspection of COI nucleotide change with reference to the simultaneous analysis tree revealed AT and codon biases. These biases were not associated with branch support measures. Therefore, it is difficult to predict incongruence or identify its cause. Exclusion of data is ill advised because every character is potentially parsimony informative.     

Exploring data interaction and nucleotide alignment in a multiple gene analysis of Ips (Coleoptera: Scolytinae)

The possibility of gene tree incongruence in a species-level phylogenetic analysis of the genus Ips (Coleoptera: Scolytidae) was investigated based on mitochondrial 16S rRNA (16S) and nuclear elongation factor-1 alpha (EF-1 alpha) sequences, and existing cytochrome oxidase I (COI) and nonmolecular data sets. Separate cladistic analyses of the data partitions resulted in partially discordant most-parsimonious trees but revealed only low conflict of the phylogenetic signal. Interactions among data partitions, which differed in the extent of sequence divergence (COI > 16S > EF-1 alpha), base composition, and homoplasy, revealed that much of the branch support emerges only in the simultaneous analysis, particularly for deeper nodes in the tree, which are almost entirely supported through "hidden support" (sensu Gatesy et al., Cladistics 15:271-313, 1999). Apparent incongruence between data partitions is in part due to suboptimal alignments and bias of character transformations, but little evidence supports invoking incongruent phylogenetic histories of genetic loci. There is also no justification for eliminating or downweighting gene partitions on the basis of their apparent homoplasy or incongruence with other partitions, because the signal emerges only in the interaction of all data. In comparison with traditional taxonomy, the pini, plastographus, and perturbatus groups are polyphyletic, whereas the grandicollis group is monophyletic except for inclusion of the (monophyletic) calligraphus group. The latidens group and some European species are distantly related and closer to other genera within Ipini. Our robust cladogram was used to revise the classification of Ips. We provide new diagnoses for Ips and four subgeneric taxa.     

A priori assumptions about characters as a cause of incongruence between molecular and morphological hypotheses of primate interrelationships

When molecules and morphology produce incongruent hypotheses of primate interrelationships, the data are typically viewed as incompatible, and molecular hypotheses are often considered to be better indicators of phylogenetic history. However, it has been demonstrated that the choice of which taxa to include in cladistic analysis as well as assumptions about character weighting, character state transformation order, and outgroup choice all influence hypotheses of relationships and may positively influence tree topology, so that relationships between extant taxa are consistent with those found using molecular data. Thus, the source of incongruence between morphological and molecular trees may lie not in the morphological data themselves but in assumptions surrounding the ways characters evolve and their impact on cladistic analysis. In this study, we investigate the role that assumptions about character polarity and transformation order play in creating incongruence between primate phylogenies based on morphological data and those supported by multiple lines of molecular data. By releasing constraints imposed on published morphological analyses of primates from disparate clades and subjecting those data to parsimony analysis, we test the hypothesis that incongruence between morphology and molecules results from inherent flaws in morphological data. To quantify the difference between incongruent trees, we introduce a new method called branch slide distance (BSD). BSD mitigates many of the limitations attributed to other tree comparison methods, thus allowing for a more accurate measure of topological similarity. We find that releasing a priori constraints on character behavior often produces trees that are consistent with molecular trees. Case studies are presented that illustrate how congruence between molecules and unconstrained morphological data may provide insight into issues of polarity, transformation order, homology, and homoplasy.     

Incongruence of mitochondrial and nuclear gene trees in the Carabid beetles Ohomopterus

We studied the molecular phylogeny of the carabid subgenus Ohomopterus (genus Carabus), using two mitochondrial (mt) DNA regions (16SrRNA and NADH dehydrogenase subunit 5) and three nuclear DNA regions (wingless, phosphoenolpyruvate carboxykinase, and an anonymous locus). We revisited the previously reported incongruence between the distribution of mtDNA markers and morphologically defined species (Su et al., 1996; J. Mol. Evol. 43:662-671), which those authors attributed to "type switching", a concerted change in many morphological characters that results in the repeated evolution of a particular morphological type. Our mtDNA gene tree obtained from 44 individuals representing all 15 currently recognized species of Ohomopterus revealed that haplotypes isolated from individuals of a single "species" were frequently separated into distant clades, confirming the previous report. The three nuclear markers generally conformed better-with the morphologically defined species than did the mitochondrial markers. The phylogenetic signal in mtDNA and nuclear DNA data differed strongly, and these two partitions were significantly incongruent with each other according to the incongruence length difference test of Farris et al. (1994; Cladistics 10:315-320), although the three nuclear partitions were not homogeneous either. Our results did not support the type-switching hypothesis that had been proposed to fit the morphological data to the mitochondrial gene tree: The incongruence of the mtDNA tree with other nuclear markers indicates that the mtDNA-based tree does not reflect species history any better than the morphological data do. Incongruence of gene trees in Ohomopterus may have been promoted by the complex processes of geographic isolation and hybridization in the Japanese Archipelago that have led to occasional gene flow and recombination between separated entities. The occurrence of reticulate patterns in this group is intriguing, because species of Ohomopterus exhibit extremely divergent genitalic structures that represent a highly efficient reproductive isolation mechanism.     

Molecular and morphological phylogenetic analysis of an insular radiation in Pacific black flies (Simulium)

Ecological adaptation within islands may have figured prominently in the insular radiation of black flies (subgenus Inseliellum) in the Society Islands, French Polynesia. To aid in understanding the sequence of ecological shifts in this group, we have constructed a phylogeny by using morphology, the cytochrome oxidase I (COI) gene, and the small ribosomal subunit (12S) gene. The strong influence of COI on the combined analysis tree was evident from its contribution to the partitioned Bremer support (62%). The net effect of including 12S was to reduce overall tree support. Different character sets resolved different portions of the combined analysis tree, with COI resolving recent lineages, 12S resolving basal relationships, and morphology supporting the monophyly of taxa having smaller larval feeding fans (oviceps group). The Partition Homogeneity and Kashino-Hasegawa tests indicated significant incongruence between morphological and mitochondrial data. The Templeton test revealed that morphology and the combined (COI + 12S) mitochondrial data were incongruent. This conflict stems primarily from disagreement over the monophyly of taxa having much smaller larval feeding fans. Either convergence in a subset of morphological characters, low phylogenetic signal among mitochondrial sequences, or lineage-sorting causing the mitochondrial data to track an incorrect evolutionary history may be responsible for these results.     

Resolving the relationships of apid bees (Hymenoptera: Apidae) through a direct optimization sensitivity analysis of molecular,morphological, and behavioural characters

Phylogenetic analyses that incorporate the most character information also provide the most explanatory power. Here I demonstrate the value of such an approach through a direct optimization sensitivity analysis of apid bee phylogeny. Whereas prior studies have relied solely on one class of data or the other, this analysis combines previously published molecular, morphological, and behavioural characters into a single supermatrix. The final dataset includes 191 ingroup and 30 outgroup taxa, and includes data from seven unaligned gene sequences (18S, 28S, wingless, EF1‐α, polII, Nak, LW rhodopsin), 209 adult and larval morphological characters, and two behavioural characters. Nine different sets of transformation cost parameters are evaluated, along with their relative degrees of character incongruence. The preferred parameter set returns a strict consensus tree somewhat similar to, but more resolved than, a previous parsimony tree based on molecules alone. I also describe the effects of including EF1‐α and LW rhodopsin intron sequences on the outcome of the direct optimization analysis. By accounting for more evidence, this study provides the most comprehensive treatment yet of apid phylogenetic relationships.     

Mitochondrial DNA,morphology, and the phylogenetic relationships of Antarctic icefishes (Notothenioidei: Channichthyidae)

The Channichthyidae is a lineage of 16 species in the Notothenioidei, a clade of fishes that dominate Antarctic near-shore marine ecosystems with respect to both diversity and biomass. Among four published studies investigating channichthyid phylogeny, no two have produced the same tree topology, and no published study has investigated the degree of phylogenetic incongruence between existing molecular and morphological datasets. In this investigation we present an analysis of channichthyid phylogeny using complete gene sequences from two mitochondrial genes (ND2 and 16S) sampled from all recognized species in the clade. In addition, we have scored all 58 unique morphological characters used in three previous analyses of channichthyid phylogenetic relationships. Data partitions were analyzed separately to assess the amount of phylogenetic resolution provided by each dataset, and phylogenetic incongruence among data partitions was investigated using incongruence length difference (ILD) tests. We utilized a parsimony-based version of the Shimodaira-Hasegawa test to determine if alternative tree topologies are significantly different from trees resulting from maximum parsimony analysis of the combined partition dataset. Our results demonstrate that the greatest phylogenetic resolution is achieved when all molecular and morphological data partitions are combined into a single maximum parsimony analysis. Also, marginal to insignificant incongruence was detected among data partitions using the ILD. Maximum parsimony analysis of all data partitions combined results in a single tree, and is a unique hypothesis of phylogenetic relationships in the Channichthyidae. In particular, this hypothesis resolves the phylogenetic relationships of at least two species (Channichthys rhinoceratus and Chaenocephalus aceratus), for which there was no consensus among the previous phylogenetic hypotheses. The combined data partition dataset provides substantial statistical power to discriminate among alternative hypotheses of channichthyid relationships. These findings suggest the optimal strategy for investigating the phylogenetic relationships of channichthyids is one that uses all available phylogenetic data in analyses of combined data partitions.     

A data based parsimony method of cophylogenetic analysis

Phylogenies of closely interacting groups, such as hosts and parasites, are seldom completely congruent. Incongruence can arise from biologically meaningful differences in the histories of the two groups, or can be generated by artifactual differences that are merely the result of incorrect phylogenies with weakly supported nodes. We present a method that distinguishes between these sources of incongruence and identifies lineages that are responsible for significant differences between phylogenies. We use the logic of conditional combination in that we first test for statistically significant incongruence using the partition homogeneity test. Then we remove all possible combinations of taxa until a non-significant result of this test is achieved. Finally, we construct a 'combined evidence' phylogeny and then reposition the incongruent taxa. This method produces trees for final comparison using reconciliation methods, but it includes only as many incongruence events as can be statistically justified from the data sets. We apply this method to a host–parasite (gopher–louse) data set and identify many fewer incongruence events than do topology based analyses alone. Our method is broadly applicable to comparisons of phylogenies of interacting taxa, such as hosts and parasites, or mutualists. The method should also be useful for other problems involving comparisons of phylogenies, such as multiple gene trees or cladistic biogeography.     

Resolution of the Hylobates phylogeny: congruence of mitochondrial D-loop sequences with molecular, behavioral, and morphological data sets

Gibbons of the genus Hylobates likely speciated very rapidly following isolation by rising sea levels during the Pleistocene. We sequenced the hypervariable region I (HV-I) of the mitochondrial D-loop to reconstruct the phylogeny of this group. Although the results clearly supported monophyly of each of the six species, the relationships among them were not clearly resolved by these data alone. A homogeneity test against published data sets of a coding mitochondrial locus (ND3-ND4 region), behavioral characters (vocalizations), and morphological traits (including skeletal and soft tissue anatomy) revealed no significant incongruence, and combining them resulted in a phylogenetic tree with much stronger support. The Kloss's gibbon (H. klossii), long considered a primitive taxon based on morphology, shares many molecular and vocal characteristics with the Javan gibbon (H. moloch), and appear as the most recently derived species. The northernmost species (H. lar and H. pileatus) are the most basal taxa. These data suggest that ancestral gibbons radiated from north to south. Unlike other markers, the HV-I region can accurately identify members of different gibbon species much like a DNA barcode, with potential applications to conservation.     

Genome-scale phylogenetics: inferring the plant tree of life from 18,896 gene trees

Phylogenetic analyses using genome-scale data sets must confront incongruence among gene trees, which in plants is exacerbated by frequent gene duplications and losses. Gene tree parsimony (GTP) is a phylogenetic optimization criterion in which a species tree that minimizes the number of gene duplications induced among a set of gene trees is selected. The run time performance of previous implementations has limited its use on large-scale data sets. We used new software that incorporates recent algorithmic advances to examine the performance of GTP on a plant data set consisting of 18,896 gene trees containing 510,922 protein sequences from 136 plant taxa (giving a combined alignment length of >2.9 million characters). The relationships inferred from the GTP analysis were largely consistent with previous large-scale studies of backbone plant phylogeny and resolved some controversial nodes. The placement of taxa that were present in few gene trees generally varied the most among GTP bootstrap replicates. Excluding these taxa either before or after the GTP analysis revealed high levels of phylogenetic support across plants. The analyses supported magnoliids sister to a eudicot + monocot clade and did not support the eurosid I and II clades. This study presents a nuclear genomic perspective on the broad-scale phylogenic relationships among plants, and it demonstrates that nuclear genes with a history of duplication and loss can be phylogenetically informative for resolving the plant tree of life.     

Comparable disparity in the appendicular skeleton across the fish–tetrapod transition,and the morphological gap between fish and tetrapod postcrania

Appendicular skeletal traits are used to quantify changes in morphological disparity and morphospace occupation across the fish–tetrapod transition and to explore the informativeness of different data partitions in phylogeny reconstruction. Anterior appendicular data yield trees that differ little from those built from the full character set, whilst posterior appendicular data result in considerable loss of phylogenetic resolution and tree branch rearrangements. Overall, there is a significant incongruence in the signals associated with pectoral and pelvic data. The appendicular skeletons of fish and tetrapods attain similar levels of morphological disparity (at least when data are rarefied at the maximum sample size for fish in our study) and occupy similarly sized regions of morphospace. However, fish appear more dispersed in morphospace than tetrapods do. All taxa show a heterogeneous distribution in morphospace, and there is a clear separation between fish and tetrapods despite the presence of several evolutionarily intermediate taxa.     

The relative lability of floral vs. non-floral characters and a morphological phylogenetic analysis of Cobaea (Polemoniaceae)

The relative levels of lability in floral vs. vegetative characters have been suggested to give insight into the mechanisms of adaptation and speciation. Cobaea (Polemoniaceae) exhibits a remarkable diversity in floral form. A morphology-based phylogeny was constructed and is congruent in most regards with a previously published molecular phylogeny, but significant incongruence was found in the placement of two taxa. A combined analysis was performed, excluding C. aquatoriensis and C. lutea, and was nearly identical to the molecular analysis. The phylogenies are compatible with a recently published classification. Incongruencies between the phylogenies within section Rosenbergia may have profound implications for floral evolution. A simple method for testing the levels of homoplasy, as an indicator of lability between two classes of data, is developed and used to test for differences between floral and non-floral classes of characters. The method is based on repeated randomization of the data into partitions and recalculation of a test statistic defined as the difference between consistency index (CI) values of the data in each partition over a given tree. The null hypothesis that the floral characters are no more homoplasious than the non-floral characters cannot be rejected (P=0.09).     

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.     

Phylogenetic incongruence among oncogenic genital alpha human papillomaviruses

The human papillomaviruses (HPVs) have long been thought to follow a monophyletic pattern of evolution with little if any evidence for recombination between genomes. On the basis of this model, both oncogenicity and tissue tropism appear to have evolved once. Still, no systematic statistical analyses have shown whether monophyly is the rule across all HPV open reading frames (ORFs). We conducted a taxonomic analysis of 59 mucosal/genital HPVs using whole-genome and sliding-window similarity measures; maximum-parsimony, neighbor-joining, and Bayesian phylogenetic analyses; and localized incongruence length difference (LILD) analyses. The algorithm for the LILD analyses localized incongruence by calculating the tree length differences between constrained and unconstrained nodes in a total-evidence tree across all HPV ORFs. The process allows statistical evaluation of every ORF/node pair in the total-evidence tree. The most significant incongruence was observed at the putative high-risk (i.e., cancer-associated) node, the common oncogenic ancestor for alpha HPV species 9 (e.g., HPV type 16 [HPV16]), 11, 7 (e.g., HPV18), 5, and 6. Although these groups share early-gene homology, including high degrees of similarity among E6 and E7, groups 9 and 11 diverge from groups 7, 5, and 6 with respect to L2 and L1. The HPV species groups primarily associated with cervical and anogenital cancers appear to follow two distinct evolutionary paths, one conferred by the early genes and another by the late genes. The incongruence in the genital HPV phylogeny could have occurred from an early recombination event, an ecological niche change, and/or asymmetric genome convergence driven by intense selection. These data indicate that the phylogeny of the oncogenic HPVs is complex and that their evolution may not be monophyletic across all genes.     

Canonical phylogenetic ordination

A phylogenetic comparative method is proposed for estimating historical effects on comparative data using the partitions that compose a cladogram, i.e., its monophyletic groups. Two basic matrices, Y and X, are defined in the context of an ordinary linear model. Y contains the comparative data measured over t taxa. X consists of an initial tree matrix that contains all the xj monophyletic groups (each coded separately as a binary indicator variable) of the phylogenetic tree available for those taxa. The method seeks to define the subset of groups, i.e., a reduced tree matrix, that best explains the patterns in Y. This definition is accomplished via regression or canonical ordination (depending on the dimensionality of Y) coupled with Monte Carlo permutations. It is argued here that unrestricted permutations (i.e., under an equiprobable model) are valid for testing this specific kind of groupwise hypothesis. Phylogeny is either partialled out or, more properly, incorporated into the analysis in the form of component variation. Direct extensions allow for testing ecomorphological data controlled by phylogeny in a variation partitioning approach. Currently available statistical techniques make this method applicable under most univariate/multivariate models and metrics; two-way phylogenetic effects can be estimated as well. The simplest case (univariate Y), tested with simulations, yielded acceptable type I error rates. Applications presented include examples from evolutionary ethology, ecology, and ecomorphology. Results showed that the new technique detected previously overlooked variation clearly associated with phylogeny and that many phylogenetic effects on comparative data may occur at particular groups rather than across the entire tree.