Alignment and Topological Accuracy of the Direct Optimization approach via POY and Traditional Phylogenetics via ClustalW + PAUP*

Evidence of absence treated as absence of evidence: The effects of variation in the number and distribution of gaps treated as missing data on the results of standard maximum likelihood analysis

Although numerous studies have demonstrated the theoretical and empirical importance of treating gaps as insertion/deletion (indel) events in phylogenetic analyses, the standard approach to maximum likelihood (ML) analysis employed in the vast majority of empirical studies codes gaps as nucleotides of unknown identity ("missing data"). Therefore, it is imperative to understand the empirical consequences of different numbers and distributions of gaps treated as missing data. We evaluated the effects of variation in the number and distribution of gaps (i.e., no base, coded as IUPAC "." or "-") treated as missing data (i.e., any base, coded as "?" or IUPAC "N") in standard ML analysis. We obtained alignments with variable numbers and arrangements of gaps by aligning seven diverse empirical datasets under different gap opening costs using MAFFT. We selected the optimal substitution model for each alignment using the corrected Akaike Information Criterion in jModelTest2 and searched for optimal trees using GARLI. We also employed a Monte Carlo approach to randomly replace nucleotides with gaps (treated as missing data) in an empirical dataset to understand more precisely the effects of varying their number and distribution. To compare alignments, we developed four new indices and used several existing measures to quantify the number and distribution of gaps in all alignments. Our most important finding is that ML scores correlate negatively with gap opening costs and the amount of missing data. However, this negative relationship is not due to the increase in missing data per se-which increases ML scores-but instead to the effect of gaps on nucleotide homology. These variables also cause significant but largely unpredictable effects on tree topology.

Differences between hard and soft phylogenetic data

When building the tree of life, variability of phylogenetic signal is often accounted for by partitioning gene sequences and testing for differences. The same considerations, however, are rarely applied to morphological data, potentially undermining its use in evolutionary contexts. Here, we apply partition heterogeneity tests to 59 animal datasets to demonstrate that significant differences exist between the phylogenetic signal conveyed by 'hard' and 'soft' characters (bones, teeth and shells versus myology, integument etc). Furthermore, the morphological partitions differ significantly in their consistency relative to independent molecular trees. The observed morphological differences correspond with missing data biases, and as such their existence presents a problem not only for phylogeny reconstruction, but also for interpretations of fossil data. Evolutionary inferences drawn from clades in which hard, readily fossilizable characters are relatively less consistent and different from other morphology (mammals, bivalves) may be less secure. More secure inferences might be drawn from the fossil record of clades that exhibit fewer differences, or exhibit more consistent hard characters (fishes, birds). In all cases, it will be necessary to consider the impact of missing data on empirical data, and the differences that exist between morphological modules.

Unique Strain of Rickettsia parkeri Associated with the Hard Tick Dermacentor parumapertus Neumann in the Western United States

In 1953, investigators at the Rocky Mountain Laboratories in Hamilton, MT, described the isolation of a spotted fever group Rickettsia (SFGR) species from Dermacentor parumapertus ticks collected from black-tailed jackrabbits (Lepus californicus) in northern Nevada. Several decades later, investigators characterized this SFGR (designated the parumapertus agent) by using mouse serotyping methods and determined that it represented a distinct rickettsial serotype closely related to Rickettsia parkeri; nonetheless, the parumapertus agent was not further characterized or studied. To our knowledge, no isolates of the parumapertus agent remain in any rickettsial culture collection, which precludes contemporary phylogenetic placement of this enigmatic SFGR. To rediscover the parumapertus agent, adult-stage D. parumapertus ticks were collected from black-tailed jackrabbits shot or encountered as roadkills in Arizona, Utah, or Texas from 2011 to 2016. A total of 339 ticks were collected and evaluated for infection with Rickettsia species. Of 112 D. parumapertus ticks collected in south Texas, 16 (14.3%) contained partial ompA sequences with the closest identity (99.6%) to Rickettsia sp. strain Atlantic rainforest Aa46, an SFGR that is closely related or identical to an SFGR species that causes a mild rickettsiosis in several states of Brazil. A pure isolate, designated strain Black Gap, was cultivated in Vero E6 cells, and sequence analysis of the rrs, gltA, sca0, sca5, and sca4 genes also revealed the closest genetic identity to Rickettsia sp. Atlantic rainforest Aa46. Phylogenetic analysis of the five concatenated rickettsial genes place Rickettsia sp. strain Black Gap and Rickettsia sp. Atlantic rainforest Aa46 with R. parkeri in a distinct and well-supported clade.IMPORTANCE We suggest that Rickettsia sp. Black Gap and Rickettsia sp. Atlantic rainforest Aa46 represent nearly identical strains of R. parkeri and that Rickettsia sp. Black Gap or a very similar strain of R. parkeri represents the parumapertus agent. The close genetic relatedness among these taxa, as well as the response of guinea pigs infected with the Black Gap strain, suggests that R. parkeri Black Gap could cause disease in humans. The identification of this organism could also account, at least in part, for the remarkable differences in severity ascribed to Rocky Mountain spotted fever (RMSF) among various regions of the American West during the early 20th century. We suggest that the wide variation in case fatality rates attributed to RMSF could have occurred by the inadvertent inclusion of cases of milder disease caused by R. parkeri Black Gap.

Progress, pitfalls and parallel universes: a history of insect phylogenetics

The phylogeny of insects has been both extensively studied and vigorously debated for over a century. A relatively accurate deep phylogeny had been produced by 1904. It was not substantially improved in topology until recently when phylogenomics settled many long-standing controversies. Intervening advances came instead through methodological improvement. Early molecular phylogenetic studies (1985-2005), dominated by a few genes, provided datasets that were too small to resolve controversial phylogenetic problems. Adding to the lack of consensus, this period was characterized by a polarization of philosophies, with individuals belonging to either parsimony or maximum-likelihood camps; each largely ignoring the insights of the other. The result was an unfortunate detour in which the few perceived phylogenetic revolutions published by both sides of the philosophical divide were probably erroneous. The size of datasets has been growing exponentially since the mid-1980s accompanied by a wave of confidence that all relationships will soon be known. However, large datasets create new challenges, and a large number of genes does not guarantee reliable results. If history is a guide, then the quality of conclusions will be determined by an improved understanding of both molecular and morphological evolution, and not simply the number of genes analysed.

The fine structure of the rectal pads of Zorotypus caudelli Karny (Zoraptera, Insecta)

The rectal pads of a species of the controversial polyneopteran order Zoraptera were examined using histological sections and TEM micrographs. Six pads are present along the thin rectal epithelium. Each pad consists of a few large principal cells surrounded by flattened junctional cells, which extend also beneath the principal cells. The cells are lined by a thin apical cuticle. No basal cells and no cavity have been observed beneath the pad. Principal cells have a regular layer of apical microvilli and are joined by intercellular septate junctions, which are interrupted by short dilatations of the intercellular space. At these levels the two adjacent plasma membranes are joined by short zonulae adhaerentes. In the cytoplasm, a rich system of strict associations between lateral plasma membranes and mitochondria forms scalariform junctions. Rectal pads share ultrastructural features with similar excretory organs of several neopteran groups, in particular with Blattodea (roaches and termites) and Thysanoptera, and are involved in fluid reabsorption and ion regulation.

Comparative morphology of spermatozoa and reproductive systems of zorapteran species from different world regions (Insecta, Zoraptera)

The male and female reproductive apparatus of Zorotypus magnicaudelli (Malaysia), Zorotypus huxleyi (Ecuador) and Zorotypus weidneri (Brazil) were examined and documented in detail. The genital apparatus and sperm of the three species show only minor differences. The testes are larger in Z. magnicaudelli. Z. huxleyi lacks the helical appendage in the accessory glands. A long cuticular flagellum is present in Z. magnicaudelli and in the previously studied Zorotypus caudelli like in several other species, whereas it is absent in Z. weidneri, Z. huxleyi, Zorotypus hubbardi, Zorotypus impolitus and Zorotypus guineensis. Characteristic features of the very similar sperm are the presence of: a) two dense arches above the axoneme; b) a 9 + 9+2 axoneme with detached subtubules A and B of doublets 1 and 6; c) the axonemal end degenerating with enlarging accessory tubules; d) accessory tubules with 17 protofilaments; e) three accessory bodies beneath the axoneme; and f) two mitochondrial derivatives of equal shape. The first characteristic (a) is unknown outside of Zoraptera and possibly autapomorphic. The sperm structure differs distinctly in Z. impolitus and Z. hubbardi, which produce giant sperm and possess a huge spermatheca. The presence of the same sperm type in species either provided with a sclerotized coiled flagellum in males or lacking this structure indicates that a different organization of the genital apparatus does not necessarily affect the sperm structure. The flagellum and its pouch has probably evolved within Zoraptera, but it cannot be excluded that it is a groundplan feature and was reduced several times. The fossil evidence and our findings suggest that distinct modifications in the genital apparatus occurred before the fragmentation of the Gondwanan landmass in the middle Cretaceous.

PCR-based molecular characterization of Toxocara spp. using feces of stray cats: a study from Southwest Iran

Feces of stray cat are potential sources of gastrointestinal parasites and play a crucial role in spreading and transmitting parasite eggs, larvae, and oocysts through contamination of soil, food, or water. In this study, we investigated the prevalence of Toxocara spp. infection in stray cats in Ahvaz city, southwest Iran. Eggs of Toxocara spp. in feces of stray cats were detected by the sucrose flotation method, and identification was conducted by polymerase chain reaction (PCR) and DNA sequencing. Of the 140 fecal samples that were randomly collected from public environments during the months of January to May 2012, 45% were found to harbour Toxocara spp. eggs. The highest prevalence of Toxocara spp. eggs was found in the central area of Ahvaz city (28.6%). T. canis eggs were found in 4 (6.34%) of the 63 positive samples. Stray cats are found in parks, playgrounds, and other public places and may be a potential contamination risk. Identification of Toxocara spp. using molecular methods is sufficiently sensitive to detect low levels of parasites and identify the different Toxocara spp. in feces. The relatively high prevalence of Toxocara spp. infection may continue to increase due to lack of effective environmental hygiene control in Iran. Consequently, there is a need to plan adequate programs to detect, identify, and control this infection as well as stray cats in the region.

Local search for the generalized tree alignment problem

BACKGROUND

A phylogeny postulates shared ancestry relationships among organisms in the form of a binary tree. Phylogenies attempt to answer an important question posed in biology: what are the ancestor-descendent relationships between organisms? At the core of every biological problem lies a phylogenetic component. The patterns that can be observed in nature are the product of complex interactions, constrained by the template that our ancestors provide. The problem of simultaneous tree and alignment estimation under Maximum Parsimony is known in combinatorial optimization as the Generalized Tree Alignment Problem (GTAP). The GTAP is the Steiner Tree Problem for the sequence edit distance. Like many biologically interesting problems, the GTAP is NP-Hard. Typically the Steiner Tree is presented under the Manhattan or the Hamming distances.

RESULTS

Experimentally, the accuracy of the GTAP has been subjected to evaluation. Results show that phylogenies selected using the GTAP from unaligned sequences are competitive with the best methods and algorithms available. Here, we implement and explore experimentally existing and new local search heuristics for the GTAP using simulated and real data.

CONCLUSIONS

The methods presented here improve by more than three orders of magnitude in execution time the best local search heuristics existing to date when applied to real data.

Re-mind the gap! Insertion - deletion data reveal neglected phylogenetic potential of the nuclear ribosomal internal transcribed spacer (ITS) of fungi

Rapidly evolving, indel-rich phylogenetic markers play a pivotal role in our understanding of the relationships at multiple levels of the tree of life. There is extensive evidence that indels provide conserved phylogenetic signal, however, the range of phylogenetic depths for which gaps retain tree signal has not been investigated in detail. Here we address this question using the fungal internal transcribed spacer (ITS), which is central in many phylogenetic studies, molecular ecology, detection and identification of pathogenic and non-pathogenic species. ITS is repeatedly criticized for indel-induced alignment problems and the lack of phylogenetic resolution above species level, although these have not been critically investigated. In this study, we examined whether the inclusion of gap characters in the analyses shifts the phylogenetic utility of ITS alignments towards earlier divergences. By re-analyzing 115 published fungal ITS alignments, we found that indels are slightly more conserved than nucleotide substitutions, and when included in phylogenetic analyses, improved the resolution and branch support of phylogenies across an array of taxonomic ranges and extended the resolving power of ITS towards earlier nodes of phylogenetic trees. Our results reconcile previous contradicting evidence for the effects of data exclusion: in the case of more sophisticated indel placement, the exclusion of indel-rich regions from the analyses results in a loss of tree resolution, whereas in the case of simpler alignment methods, the exclusion of gapped sites improves it. Although the empirical datasets do not provide to measure alignment accuracy objectively, our results for the ITS region are consistent with previous simulations studies alignment algorithms. We suggest that sophisticated alignment algorithms and the inclusion of indels make the ITS region and potentially other rapidly evolving indel-rich loci valuable sources of phylogenetic information, which can be exploited at multiple taxonomic levels.

The tree alignment problem

BACKGROUND

The inference of homologies among DNA sequences, that is, positions in multiple genomes that share a common evolutionary origin, is a crucial, yet difficult task facing biologists. Its computational counterpart is known as the multiple sequence alignment problem. There are various criteria and methods available to perform multiple sequence alignments, and among these, the minimization of the overall cost of the alignment on a phylogenetic tree is known in combinatorial optimization as the Tree Alignment Problem. This problem typically occurs as a subproblem of the Generalized Tree Alignment Problem, which looks for the tree with the lowest alignment cost among all possible trees. This is equivalent to the Maximum Parsimony problem when the input sequences are not aligned, that is, when phylogeny and alignments are simultaneously inferred.

RESULTS

For large data sets, a popular heuristic is Direct Optimization (DO). DO provides a good tradeoff between speed, scalability, and competitive scores, and is implemented in the computer program POY. All other (competitive) algorithms have greater time complexities compared to DO. Here, we introduce and present experiments a new algorithm Affine-DO to accommodate the indel (alignment gap) models commonly used in phylogenetic analysis of molecular sequence data. Affine-DO has the same time complexity as DO, but is correctly suited for the affine gap edit distance. We demonstrate its performance with more than 330,000 experimental tests. These experiments show that the solutions of Affine-DO are close to the lower bound inferred from a linear programming solution. Moreover, iterating over a solution produced using Affine-DO shows little improvement.

CONCLUSIONS

Our results show that Affine-DO is likely producing near-optimal solutions, with approximations within 10% for sequences with small divergence, and within 30% for random sequences, for which Affine-DO produced the worst solutions. The Affine-DO algorithm has the necessary scalability and optimality to be a significant improvement in the real-world phylogenetic analysis of sequence data.

The morphology and evolution of the female postabdomen of Holometabola (Insecta)

In the present article homology issues, character evolution and phylogenetic implications related to the female postabdomen of the holometabolan insects are discussed, based on an earlier analysis of a comprehensive morphological data set. Hymenoptera, the sistergroup of the remaining Holometabola, are the only group where the females have retained a fully developed primary ovipositor of the lepismatid type. There are no characters of the female abdomen supporting a clade Coleopterida + Neuropterida. The invagination of the terminal segments is an autapomorphy of Coleoptera. The ovipositor is substantially modified in Raphidioptera and distinctly reduced in Megaloptera and Neuroptera. The entire female abdomen is extremely simplified in Strepsiptera. The postabdomen is tapering posteriorly in Mecopterida and retractile in a telescopic manner (oviscapt). The paired ventral sclerites of segments VIII and IX are preserved, but valvifers and valvulae are not distinguishable. In Amphiesmenoptera sclerotizations derived from the ventral appendages VIII are fused ventromedially, forming a solid plate, and the appendages IX are reduced. The terminal segments are fused and form a terminal unit which bears the genital opening subapically. The presence of two pairs of apophyses and the related protraction of the terminal unit by muscle force are additional autapomorphies, as is the fusion of the rectum with the posterior part of the genital chamber (cloaca). Antliophora are supported by the presence of a transverse muscle between the ventral sclerites of segment VIII. Secondary egg laying tubes have evolved independently within Boreidae (absent in Caurinus) and in Tipulomorpha. The loss of two muscle associated with the genital chamber are likely autapomorphies of Diptera. The secondary loss of the telescopic retractability of the postabdomen is one of many autapomorphies of Siphonaptera.

Treelength optimization for phylogeny estimation

The standard approach to phylogeny estimation uses two phases, in which the first phase produces an alignment on a set of homologous sequences, and the second phase estimates a tree on the multiple sequence alignment. POY, a method which seeks a tree/alignment pair minimizing the total treelength, is the most widely used alternative to this two-phase approach. The topological accuracy of trees computed under treelength optimization is, however, controversial. In particular, one study showed that treelength optimization using simple gap penalties produced poor trees and alignments, and suggested the possibility that if POY were used with an affine gap penalty, it might be able to be competitive with the best two-phase methods. In this paper we report on a study addressing this possibility. We present a new heuristic for treelength, called BeeTLe (Better Treelength), that is guaranteed to produce trees at least as short as POY. We then use this heuristic to analyze a large number of simulated and biological datasets, and compare the resultant trees and alignments to those produced using POY and also maximum likelihood (ML) and maximum parsimony (MP) trees computed on a number of alignments. In general, we find that trees produced by BeeTLe are shorter and more topologically accurate than POY trees, but that neither POY nor BeeTLe produces trees as topologically accurate as ML trees produced on standard alignments. These findings, taken as a whole, suggest that treelength optimization is not as good an approach to phylogenetic tree estimation as maximum likelihood based upon good alignment methods.

Standard maximum likelihood analyses of alignments with gaps can be statistically inconsistent

BackgroundMost statistical methods for phylogenetic estimation in use today treat a gap (generally representing an insertion or deletion, i.e., indel) within the input sequence alignment as missing data. However, the statistical properties of this treatment of indels have not been fully investigated.ResultsWe prove that maximum likelihood phylogeny estimation, treating indels as missing data, can be statistically inconsistent for a general (and rather simple) model of sequence evolution, even when given the true alignment. Therefore, accurate phylogeny estimation cannot be guaranteed for maximum likelihood analyses, even given arbitrarily long sequences, when indels are present and treated as missing data.ConclusionsOur result shows that the standard statistical techniques used to estimate phylogenies from sequence alignments may have unfavorable statistical properties, even when the sequence alignment is accurate and the assumed substitution model matches the generation model. This suggests that the recent research focus on developing statistical methods that treat indel events properly is an important direction for phylogeny estimation.

The fine structure of the female reproductive system of Zorotypus caudelli Karny (Zoraptera)

The general structure of the female genital system of Zorotypus caudelli is described. The ovarioles are of the panoistic type. Due to the reduction of the envelope (tunica externa) the ovarioles are in direct contact with the hemolymph like in some other insect groups, Plecoptera included. The calices are much larger in Z. caudelli then in Zorotypus hubbardi and their epithelial cells produce large amounts of secretions, probably protecting the surface of the eggs deposited on the substrate. Eggs taken from the calyx bear a series of long fringes, which are missing in the eggs found in the ovariole, and in other zorapteran species. The long sperm of Z. caudelli and the long spermathecal duct are likely related to a sexual isolating mechanism (cryptic female choice), impeding female re-mating. The apical receptacle and the spermathecal duct - both of ectodermal origin - consist of three cell types. In addition to the cells beneath the cuticle lining the lumen, two other cell types are visible: secretory and canal cells. The cytoplasm of the former is rich in rough endoplasmic reticulum cisterns and Golgi complexes, which produce numerous discrete dense secretory bodies. These products are released into the receiving canal crossing the extracellular cavity of secretory cells, extending over a series of long microvilli. The secretion is transported towards the lumen of the apical receptacle of the spermatheca or to that of the spermathecal duct by a connecting canal formed by the canal cells. It is enriched by material produced by the slender canal cells. Before mating, the sperm cells are enveloped by a thick glycocalyx produced at the level of the male accessory glands, but it is absent when they have reached the apical receptacle, and also in the spermathecal duct lumen. It is likely removed by secretions of the spermatheca. The eggs are fertilized at the level of the common oviduct where the spermathecal duct opens. Two micropyles at the dorsal side of the equator level possibly facilitate fertilization. The presence of these two micropyles is a presumably derived feature shared with Phasmatodea. The fine structure of the female reproductive system of Z. caudelli does not allow to assess the phylogenetic position at the present stage of knowledge. The enlarged calyx and the temporary presence of long fringes on the eggs are potential autapomorphies of Z. caudelli or may indicate relationships with other Zorotypus species.

BitterDB: a database of bitter compounds

Basic taste qualities like sour, salty, sweet, bitter and umami serve specific functions in identifying food components found in the diet of humans and animals, and are recognized by proteins in the oral cavity. Recognition of bitter taste and aversion to it are thought to protect the organism against the ingestion of poisonous food compounds, which are often bitter. Interestingly, bitter taste receptors are expressed not only in the mouth but also in extraoral tissues, such as the gastrointestinal tract, indicating that they may play a role in digestive and metabolic processes. BitterDB database, available at http://bitterdb.agri.huji.ac.il/bitterdb/, includes over 550 compounds that were reported to taste bitter to humans. The compounds can be searched by name, chemical structure, similarity to other bitter compounds, association with a particular human bitter taste receptor, and so on. The database also contains information on mutations in bitter taste receptors that were shown to influence receptor activation by bitter compounds. The aim of BitterDB is to facilitate studying the chemical features associated with bitterness. These studies may contribute to predicting bitterness of unknown compounds, predicting ligands for bitter receptors from different species and rational design of bitterness modulators.

The impact of multiple protein sequence alignment on phylogenetic estimation

Multiple sequence alignment is typically the first step in estimating phylogenetic trees, with the assumption being that as alignments improve, so will phylogenetic reconstructions. Over the last decade or so, new multiple sequence alignment methods have been developed to improve comparative analyses of protein structure, but these new methods have not been typically used in phylogenetic analyses. In this paper, we report on a simulation study that we performed to evaluate the consequences of using these new multiple sequence alignment methods in terms of the resultant phylogenetic reconstruction. We find that while alignment accuracy is positively correlated with phylogenetic accuracy, the amount of improvement in phylogenetic estimation that results from an improved alignment can range from quite small to substantial. We observe that phylogenetic accuracy is most highly correlated with alignment accuracy when sequences are most difficult to align, and that variation in alignment accuracy can have little impact on phylogenetic accuracy when alignment error rates are generally low. We discuss these observations and implications for future work.

PICS-Ord: unlimited coding of ambiguous regions by pairwise identity and cost scores ordination

BACKGROUND

We present a novel method to encode ambiguously aligned regions in fixed multiple sequence alignments by 'Pairwise Identity and Cost Scores Ordination' (PICS-Ord). The method works via ordination of sequence identity or cost scores matrices by means of Principal Coordinates Analysis (PCoA). After identification of ambiguous regions, the method computes pairwise distances as sequence identities or cost scores, ordinates the resulting distance matrix by means of PCoA, and encodes the principal coordinates as ordered integers. Three biological and 100 simulated datasets were used to assess the performance of the new method.

RESULTS

Including ambiguous regions coded by means of PICS-Ord increased topological accuracy, resolution, and bootstrap support in real biological and simulated datasets compared to the alternative of excluding such regions from the analysis a priori. In terms of accuracy, PICS-Ord performs equal to or better than previously available methods of ambiguous region coding (e.g., INAASE), with the advantage of a practically unlimited alignment size and increased analytical speed and the possibility of PICS-Ord scores to be analyzed together with DNA data in a partitioned maximum likelihood model.

CONCLUSIONS

Advantages of PICS-Ord over step matrix-based ambiguous region coding with INAASE include a practically unlimited number of OTUs and seamless integration of PICS-Ord codes into phylogenetic datasets, as well as the increased speed of phylogenetic analysis. Contrary to word- and frequency-based methods, PICS-Ord maintains the advantage of pairwise sequence alignment to derive distances, and the method is flexible with respect to the calculation of distance scores. In addition to distance and maximum parsimony, PICS-Ord codes can be analyzed in a Bayesian or maximum likelihood framework. RAxML (version 7.2.6 or higher that was developed for this study) allows up to 32-state ordered or unordered characters. A GTR, MK, or ORDERED model can be applied to analyse the PICS-Ord codes partition, with GTR performing slightly better than MK and ORDERED.

AVAILABILITY

An implementation of the PICS-Ord algorithm is available from http://scit.us/projects/ngila/wiki/PICS-Ord. It requires both the statistical software, R http://www.r-project.org and the alignment software Ngila http://scit.us/projects/ngila.

Recent trends in molecular phylogenetic analysis: where to next?

The acquisition of large multilocus sequence data is providing researchers with an unprecedented amount of information to resolve difficult phylogenetic problems. With these large quantities of data comes the increasing challenge regarding the best methods of analysis. We review the current trends in molecular phylogenetic analysis, focusing specifically on the topics of multiple sequence alignment and methods of tree reconstruction. We suggest that traditional methods are inadequate for these highly heterogeneous data sets and that researchers employ newer more sophisticated search algorithms in their analyses. If we are to best extract the information present in these data sets, a sound understanding of basic phylogenetic principles combined with modern methodological techniques are necessary.

Phylogenetic position of the acariform mites: sensitivity to homology assessment under total evidence

BACKGROUND

Mites (Acari) have traditionally been treated as monophyletic, albeit composed of two major lineages: Acariformes and Parasitiformes. Yet recent studies based on morphology, molecular data, or combinations thereof, have increasingly drawn their monophyly into question. Furthermore, the usually basal (molecular) position of one or both mite lineages among the chelicerates is in conflict to their morphology, and to the widely accepted view that mites are close relatives of Ricinulei.

RESULTS

The phylogenetic position of the acariform mites is examined through employing SSU, partial LSU sequences, and morphology from 91 chelicerate extant terminals (forty Acariformes). In a static homology framework, molecular sequences were aligned using their secondary structure as guide, whereby regions of ambiguous alignment were discarded, and pre-aligned sequences analyzed under parsimony and different mixed models in a Bayesian inference. Parsimony and Bayesian analyses led to trees largely congruent concerning infra-ordinal, well-supported branches, but with low support for inter-ordinal relationships. An exception is Solifugae + Acariformes (P. P = 100%, J. = 0.91). In a dynamic homology framework, two analyses were run: a standard POY analysis and an analysis constrained by secondary structure. Both analyses led to largely congruent trees; supporting a (Palpigradi (Solifugae Acariformes)) clade and Ricinulei as sister group of Tetrapulmonata with the topology (Ricinulei (Amblypygi (Uropygi Araneae))). Combined analysis with two different morphological data matrices were run in order to evaluate the impact of constraining the analysis on the recovered topology when employing secondary structure as a guide for homology establishment. The constrained combined analysis yielded two topologies similar to the exclusively molecular analysis for both morphological matrices, except for the recovery of Pedipalpi instead of the (Uropygi Araneae) clade. The standard (direct optimization) POY analysis, however, led to the recovery of trees differing in the absence of the otherwise well-supported group Solifugae + Acariformes.

CONCLUSIONS

Previous studies combining ribosomal sequences and morphology often recovered topologies similar to purely morphological analyses of Chelicerata. The apparent stability of certain clades not recovered here, like Haplocnemata and Acari, is regarded as a byproduct of the way the molecular homology was previously established using the instrumentalist approach implemented in POY. Constraining the analysis by a priori homology assessment is defended here as a way of maintaining the severity of the test when adding new data to the analysis. Although the strength of the method advocated here is keeping phylogenetic information from regions usually discarded in an exclusively static homology framework; it still has the inconvenience of being uninformative on the effect of alignment ambiguity on resampling methods of clade support estimation. Finally, putative morphological apomorphies of Solifugae + Acariformes are the reduction of the proximal cheliceral podomere, medial abutting of the leg coxae, loss of sperm nuclear membrane, and presence of differentiated germinative and secretory regions in the testis delivering their products into a common lumen.

Using the Multiple Analysis Approach to Reconstruct Phylogenetic Relationships among Planktonic Foraminifera from Highly Divergent and Length-polymorphic SSU rDNA Sequences

The high sequence divergence within the small subunit ribosomal RNA gene (SSU rDNA) of foraminifera makes it difficult to establish the homology of individual nucleotides across taxa. Alignment-based approaches so far relied on time-consuming manual alignments and discarded up to 50% of the sequenced nucleotides prior to phylogenetic inference. Here, we investigate the potential of the multiple analysis approach to infer a molecular phylogeny of all modern planktonic foraminiferal taxa by using a matrix of 146 new and 153 previously published SSU rDNA sequences. Our multiple analysis approach is based on eleven different automated alignments, analysed separately under the maximum likelihood criterion. The high degree of congruence between the phylogenies derived from our novel approach, traditional manually homologized culled alignments and the fossil record indicates that poorly resolved nucleotide homology does not represent the most significant obstacle when exploring the phylogenetic structure of the SSU rDNA in planktonic foraminifera. We show that approaches designed to extract phylogenetically valuable signals from complete sequences show more promise to resolve the backbone of the planktonic foraminifer tree than attempts to establish strictly homologous base calls in a manual alignment.

Phylogenetic support values are not necessarily informative: the case of the Serialia hypothesis (a mollusk phylogeny)

BACKGROUND

Molecular phylogenies are being published increasingly and many biologists rely on the most recent topologies. However, different phylogenetic trees often contain conflicting results and contradict significant background data. Not knowing how reliable traditional knowledge is, a crucial question concerns the quality of newly produced molecular data. The information content of DNA alignments is rarely discussed, as quality statements are mostly restricted to the statistical support of clades. Here we present a case study of a recently published mollusk phylogeny that contains surprising groupings, based on five genes and 108 species, and we apply new or rarely used tools for the analysis of the information content of alignments and for the filtering of noise (masking of random-like alignment regions, split decomposition, phylogenetic networks, quartet mapping).

RESULTS

The data are very fragmentary and contain contaminations. We show that that signal-like patterns in the data set are conflicting and partly not distinct and that the reported strong support for a "rather surprising result" (monoplacophorans and chitons form a monophylum Serialia) does not exist at the level of primary homologies. Split-decomposition, quartet mapping and neighbornet analyses reveal conflicting nucleotide patterns and lack of distinct phylogenetic signal for the deeper phylogeny of mollusks.

CONCLUSION

Even though currently a majority of molecular phylogenies are being justified with reference to the 'statistical' support of clades in tree topologies, this confidence seems to be unfounded. Contradictions between phylogenies based on different analyses are already a strong indication of unnoticed pitfalls. The use of tree-independent tools for exploratory analyses of data quality is highly recommended. Concerning the new mollusk phylogeny more convincing evidence is needed.

Simultaneous phylogeny reconstruction and multiple sequence alignment

Background

A phylogeny is the evolutionary history of a group of organisms. To date, sequence data is still the most used data type for phylogenetic reconstruction. Before any sequences can be used for phylogeny reconstruction, they must be aligned, and the quality of the multiple sequence alignment has been shown to affect the quality of the inferred phylogeny. At the same time, all the current multiple sequence alignment programs use a guide tree to produce the alignment and experiments showed that good guide trees can significantly improve the multiple alignment quality.

Results

We devise a new algorithm to simultaneously align multiple sequences and search for the phylogenetic tree that leads to the best alignment. We also implemented the algorithm as a C program package, which can handle both DNA and protein data and can take simple cost model as well as complex substitution matrices, such as PAM250 or BLOSUM62. The performance of the new method are compared with those from other popular multiple sequence alignment tools, including the widely used programs such as ClustalW and T-Coffee. Experimental results suggest that this method has good performance in terms of both phylogeny accuracy and alignment quality.

Conclusion

We present an algorithm to align multiple sequences and reconstruct the phylogenies that minimize the alignment score, which is based on an efficient algorithm to solve the median problems for three sequences. Our extensive experiments suggest that this method is very promising and can produce high quality phylogenies and alignments.

Barking up the wrong treelength: the impact of gap penalty on alignment and tree accuracy

Several methods have been developed for simultaneous estimation of alignment and tree, of which POY is the most popular. In a 2007 paper published in Systematic Biology, Ogden and Rosenberg reported on a simulation study in which they compared POY to estimating the alignment using ClustalW and then analyzing the resultant alignment using maximum parsimony. They found that ClustalW+MP outperformed POY with respect to alignment and phylogenetic tree accuracy, and they concluded that simultaneous estimation techniques are not competitive with two-phase techniques. Our paper presents a simulation study in which we focus on the NP-hard optimization problem that POY addresses: minimizing treelength. Our study considers the impact of the gap penalty and suggests that the poor performance observed for POY by Ogden and Rosenberg is due to the simple gap penalties they used to score alignment/tree pairs. Our study suggests that optimizing under an affine gap penalty might produce alignments that are better than ClustalW alignments, and competitive with those produced by the best current alignment methods. We also show that optimizing under this affine gap penalty produces trees whose topological accuracy is better than ClustalW+MP, and competitive with the current best two-phase methods.

Secondary structure of ITS2 rRNA provides taxonomic characters for systematic studies--a case in Lycoperdaceae (Basidiomycota)

The secondary structure of the ITS2 rDNA transcript (pre-rRNA) could provide information for identifying homologous nucleotide characters useful for cladistic inference of relationships. Such structure data could become taxonomic characters. This work compares the effect of several modern nucleotide alignment strategies, including those making use of structure data, on phylogenetic inference. From both the phylogenetic analyses and comparative secondary structure, implications for taxonomy and evolution of puffball fungi are discussed. Lycoperdaceae remain insufficiently resolved with present taxon and data sampling. Neither alignment allows statistically robust phylogenetic hypotheses under any current optimality criterion. The secondary structure data at this time are best used as accessory taxonomic characters as their phylogenetic resolving power and confidence in validity is limited compared with underlying nucleotide characters. We introduce a preliminary nomenclature convention to describe secondary structure for defining consensus features. These consensus structures are illustrated for the clades /Calvatia, /Handkea-Echinatum, /Vascellum, /Morganella, and /Plumbea-Paludosa (Bovista).