Barbara R. Holland

Outgroup Misplacement and Phylogenetic Inaccuracy Under a Molecular Clock—A Simulation Study

We conducted a simulation study of the phylogenetic methods UPGMA, neighbor joining, maximum parsimony, and maximum likelihood for a five-taxon tree under a molecular clock. The parameter space included a small region where maximum parsimony is inconsistent, so we tested inconsistency correction for parsimony and distance correction for neighbor joining. As expected, corrected parsimony was consistent. For these data, maximum likelihood with the clock assumption outperformed each of the other methods tested. The distance-based methods performed marginally better than did maximum parsimony and maximum likelihood without the clock assumption. Data correction was generally detrimental to accuracy, especially for short sequence lengths. We identified another region of the parameter space where, although consistent for a given method, some incorrect trees were each selected with up to twice the frequency of the correct (generating) tree for sequences of bounded length. These incorrect trees are those where the outgroup has been incorrectly placed. In addition to this problem, the placement of the outgroup sequence can have a confounding effect on the ingroup tree, whereby the ingroup is correct when using the ingroup sequences alone, but with the inclusion of the outgroup the ingroup tree becomes incorrect.

Optimizing Automated AFLP Scoring Parameters to Improve Phylogenetic Resolution

The amplified fragment length polymorphism (AFLP) technique is an increasingly popular component of the phylogenetic toolbox, particularly for plant species. Technological advances in capillary electrophoresis now allow very precise estimates of DNA fragment mobility and amplitude, and current AFLP software allows greater control of data scoring and the production of the binary character matrix. However, for AFLP to become a useful modern tool for large data sets, improvements to automated scoring are required. We design a procedure that can be used to optimize AFLP scoring parameters to improve phylogenetic resolution and demonstrate it for two AFLP scoring programs (GeneMapper and GeneMarker). In general, we found that there was a trade-off between getting more characters of lower quality and fewer characters of high quality. Conservative settings that gave the least error did not give the best phylogenetic resolution, as too many useful characters were discarded. For example, in GeneMapper, we found that bin width was a crucial parameter, and that although reducing bin width from 1.0 to 0.5 base pairs increased the error rate, it nevertheless improved resolution due to the increased number of informative characters. For our 30-taxon data sets, moving from default to optimized parameter settings gave between 3 and 11 extra internal edges with >50% bootstrap support, in the best case increasing the number of resolved edges from 14 to 25 out of a possible 27. Nevertheless, improvements to current AFLP software packages are needed to (1) make use of replicate profiles to calibrate the data and perform error calculations and (2) perform tests to optimize scoring parameters in a rigorous and automated way. This is true not only when AFLP data are used for phylogenetics, but also for other applications, including linkage mapping and population genetics.

Using supernetworks to distinguish hybridization from lineage-sorting

BackgroundA simple and widely used approach for detecting hybridization in phylogenies is to reconstruct gene trees from independent gene loci, and to look for gene tree incongruence. However, this approach may be confounded by factors such as poor taxon-sampling and/or incomplete lineage-sorting.ResultsUsing coalescent simulations, we investigated the potential of supernetwork methods to differentiate between gene tree incongruence arising from taxon sampling and incomplete lineage-sorting as opposed to hybridization. For few hybridization events, a large number of independent loci, and well-sampled taxa across these loci, we found that it was possible to distinguish incomplete lineage-sorting from hybridization using the filtered Z-closure and Q-imputation supernetwork methods. Moreover, we found that the choice of supernetwork method was less important than the choice of filtering, and that count-based filtering was the most effective filtering technique.ConclusionFiltered supernetworks provide a tool for detecting and identifying hybridization events in phylogenies, a tool that should become increasingly useful in light of current genome sequencing initiatives and the ease with which large numbers of independent gene loci can be determined using new generation sequencing technologies.

Molecular Epidemiology of Campylobacter jejuni Isolates from Wild-Bird Fecal Material in Children's Playgrounds

ABSTRACT In many countries relatively high notification rates of campylobacteriosis are observed in children under 5 years of age. Few studies have considered the role that environmental exposure plays in the epidemiology of these cases. Wild birds inhabit parks and playgrounds and are recognized carriers of Campylobacter, and young children are at greater risk of ingesting infective material due to their frequent hand-mouth contact. We investigated wild-bird fecal contamination in playgrounds in parks in a New Zealand city. A total of 192 samples of fresh and dried fecal material were cultured to determine the presence of Campylobacter spp. Campylobacter jejuni isolates were also characterized by pulsed-field gel electrophoresis (PFGE) and multilocus sequence typing (MLST), and the profiles obtained were compared with those of human isolates. C. jejuni was isolated from 12.5% of the samples. MLST identified members of clonal complexes ST-45, ST-682, and ST-177; all of these complexes have been recovered from wild birds in Europe. PFGE of ST-45 isolates resulted in profiles indistinguishable from those of isolated obtained from human cases in New Zealand. Members of the ST-177 and ST-682 complexes have been found in starlings (Sturnus vulgaris) in the United Kingdom, and these birds were common in playgrounds investigated in New Zealand in this study. We suggest that feces from wild birds in playgrounds could contribute to the occurrence of campylobacteriosis in preschool children. Further, the C. jejuni isolates obtained in this study belonged to clonal complexes associated with wild-bird populations in the northern hemisphere and could have been introduced into New Zealand in imported wild garden birds in the 19th century.

Consensus Networks: A Method for Visualising Incompatibilities in Collections of Trees

We present a method for summarising collections of phylogenetic trees that extends the notion of consensus trees. Each branch in a phylogenetic tree corresponds to a bipartition or split of the set of taxa labelling its leaves. Given a collection of phylogenetic trees, each labelled by the same set of taxa, all those splits that appear in more than a predefined threshold proportion of the trees are displayed using a median network. The complexity of this network is bounded as a function of the threshold proportion. We demonstrate the method for a collection of 5000 trees resulting from a Monte Carlo Markov Chain analysis of 37 mammal mitochondrial genomes, and also for a collection of 80 equally parsimonious trees resulting from a heuristic search on 53 human mitochondrial sequences.

Heterozygosity and functional allelic variation in the Candida albicans efflux pump genes CDR1 and CDR2

Elevated expression of the plasma membrane drug efflux pump proteins Cdr1p and Cdr2p was shown to accompany decreased azole susceptibility in Candida albicans clinical isolates. DNA sequence analysis revealed extensive allelic heterozygosity, particularly of CDR2. Cdr2p alleles showed different abilities to transport azoles when individually expressed in Saccharomyces cerevisiae. Loss of heterozygosity, however, did not accompany decreased azole sensitivity in isogenic clinical isolates. Two adjacent non‐synonymous single nucleotide polymorphisms (NS‐SNPs), G1473A and I1474V in the putative transmembrane (TM) helix 12 of CDR2, were found to be present in six strains including two isogenic pairs. Site‐directed mutagenesis showed that the TM‐12 NS‐SNPs, and principally the G1473A NS‐SNP, contributed to functional differences between the proteins encoded by the two Cdr2p alleles in a single strain. Allele‐specific PCR revealed that both alleles were equally frequent among 69 clinical isolates and that the majority of isolates (81%) were heterozygous at the G1473A/I1474V locus, a significant (P < 0.001) deviation from the Hardy–Weinberg equilibrium. Phylogenetic analysis by maximum likelihood (Paml) identified 33 codons in CDR2 in which amino acid allelic changes showed a high probability of being selectively advantageous. In contrast, all codons in CDR1 were under purifying selection. Collectively, these results indicate that possession of two functionally different CDR2 alleles in individual strains may confer a selective advantage, but that this is not necessarily due to azole resistance.

Visualizing Conflicting Evolutionary Hypotheses in Large Collections of Trees: Using Consensus Networks to Study the Origins of Placentals and Hexapods

Many phylogenetic methods produce large collections of trees as opposed to a single tree, which allows the exploration of support for various evolutionary hypotheses. However, to be useful, the information contained in large collections of trees should be summarized; frequently this is achieved by constructing a consensus tree. Consensus trees display only those signals that are present in a large proportion of the trees. However, by their very nature consensus trees require that any conflicts between the trees are necessarily disregarded. We present a method that extends the notion of consensus trees to allow the visualization of conflicting hypotheses in a consensus network. We demonstrate the utility of this method in highlighting differences amongst maximum likelihood bootstrap values and Bayesian posterior probabilities in the placental mammal phylogeny, and also in comparing the phylogenetic signal contained in amino acid versus nucleotide characters for hexapod monophyly.

Untangling Long Branches: Identifying Conflicting Phylogenetic Signals Using Spectral Analysis, Neighbor-Net, and Consensus Networks

Long-branch attraction is a well-known source of systematic error that can mislead phylogenetic methods; it is frequently invoked post hoc, upon recovering a different tree from the one expected based on prior evidence. We demonstrate that methods that do not force the data onto a single tree, such as spectral analysis, Neighbor-Net, and consensus networks, can be used to detect conflicting signals within the data, including those caused by long-branch attraction. We illustrate this approach using a set of taxa from three unambiguously monophyletic families within the Pelecaniformes: the darters, the cormorants and shags, and the gannets and boobies. These three families are universally acknowledged as forming a monophyletic group, but the relationship between the families remains contentious. Using sequence data from three mitochondrial genes (12S, ATPase 6, and ATPase 8) we demonstrate that the relationship between these three families is difficult to resolve because they are separated by a short internal branch and there are conflicting signals due to long-branch attraction, which are confounded with nonhomogeneous sequence evolution across the different genes. Spectral analysis, Neighbor-Net, and consensus networks reveal conflicting signals regarding the placement of one of the darters, with support found for darter monophyly, but also support for a conflicting grouping with the outgroup, pelicans. Furthermore, parsimony and maximum-likelihood analyses produced different trees, with one of the two most parsimonious trees not supporting the monophyly of the darters. Monte Carlo simulations, however, were not sensitive enough to reveal long-branch attraction unless the branches are longer than those actually observed. These results indicate that spectral analysis, Neighbor-Net, and consensus networks offer a powerful approach to detecting and understanding the source of conflicting signals within phylogenetic data.

Imputing Supertrees and Supernetworks from Quartets

Inferring species phylogenies is an important part of understanding molecular evolution. Even so, it is well known that an accurate phylogenetic tree reconstruction for a single gene does not always necessarily correspond to the species phylogeny. One commonly accepted strategy to cope with this problem is to sequence many genes; the way in which to analyze the resulting collection of genes is somewhat more contentious. Supermatrix and supertree methods can be used, although these can suppress conflicts arising from true differences in the gene trees caused by processes such as lineage sorting, horizontal gene transfer, or gene duplication and loss. In 2004, Huson et al. (IEEE/ACM Trans. Comput. Biol. Bioinformatics 1:151-158) presented the Z-closure method that can circumvent this problem by generating a supernetwork as opposed to a supertree. Here we present an alternative way for generating supernetworks called Q-imputation. In particular, we describe a method that uses quartet information to add missing taxa into gene trees. The resulting trees are subsequently used to generate consensus networks, networks that generalize strict and majority-rule consensus trees. Through simulations and application to real data sets, we compare Q-imputation to the matrix representation with parsimony (MRP) supertree method and Z-closure, and demonstrate that it provides a useful complementary tool.

Mosaic Genomes of the Six Major Primate Lentivirus Lineages Revealed by Phylogenetic Analyses

ABSTRACT To clarify the origin and evolution of the primate lentiviruses (PLVs), which include human immunodeficiency virus types 1 and 2 as well as their simian relatives, simian immunodeficiency viruses (SIVs), isolated from several host species, we investigated the phylogenetic relationships among the six supposedly nonrecombinant PLV lineages for which the full genome sequences are available. Employing bootscanning as an exploratory tool, we located several regions in the PLV genome that seem to have uncertain or conflicting phylogenetic histories. Phylogeny reconstruction based on distance and maximum-likelihood algorithms followed by a number of statistical tests confirms the existence of at least five putative recombinant fragments in the PLV genome with different clustering patterns. Split decomposition analysis also shows that phylogenetic relationships among PLVs may be better represented by network-based graphs, such as the ones produced by SplitsTree. Our findings not only imply that the six so-called pure PLV lineages have in fact mosaic genomes but also make more unlikely the hypothesis of cospeciation of SIVs and their simian hosts.