Opinions on multiple sequence alignment, and an empirical comparison of repeatability and accuracy between POY and structural alignment

Phenotypic characters of static homology increase phylogenetic stability under direct optimization of otherwise dynamic homology characters

Direct optimization of unaligned sequence characters provides a natural framework to explore the sensitivity of phylogenetic hypotheses to variation in analytical parameters. Phenotypic data, when combined into such analyses, are typically analyzed with static homology correspondences unlike the dynamic homology sequence data. Static homology characters may be expected to constrain the direct optimization and thus, potentially increase the similarity of phylogenetic hypotheses under different cost sets. However, whether a total-evidence approach increases the phylogenetic stability or not remains empirically largely unexplored. Here, I studied the impact of static homology data on sensitivity using six empirical data sets composed of several molecular markers and phenotypic data. The inclusion of static homology phenotypic data increased the average stability of phylogenetic hypothesis in five out of the six data sets. To investigate if any static homology characters would have similar effect, the analyses were repeated with randomized phenotypic data, and with one of the molecular markers fixed as static homology characters. These analyses had, on average, almost no effect on the phylogenetic stability, although the randomized phenotypic data sometimes resulted in even higher stability than empirical phenotypic data. The impact was related to the strength of the phylogenetic signal in the phenotypic data: higher average jackknife support of the phenotypic tree correlated with stronger stabilizing effect in the total-evidence analysis. Phenotypic data with a strong signal made the total-evidence trees topologically more similar to the phenotypic trees, thus, they constrained the dynamic homology correspondences of the sequence data. Characters that increase phylogenetic stability are particularly valuable for phylogenetic inference. These results indicate an important role and additive value of phenotypic data in increasing the stability of phylogenetic hypotheses in total-evidence analyses.

LMAP_S: Lightweight Multigene Alignment and Phylogeny eStimation

Background

Recent advances in genome sequencing technologies and the cost drop in high-throughput sequencing continue to give rise to a deluge of data available for downstream analyses. Among others, evolutionary biologists often make use of genomic data to uncover phenotypic diversity and adaptive evolution in protein-coding genes. Therefore, multiple sequence alignments (MSA) and phylogenetic trees (PT) need to be estimated with optimal results. However, the preparation of an initial dataset of multiple sequence file(s) (MSF) and the steps involved can be challenging when considering extensive amount of data. Thus, it becomes necessary the development of a tool that removes the potential source of error and automates the time-consuming steps of a typical workflow with high-throughput and optimal MSA and PT estimations.

Results

We introduce LMAP_S (Lightweight Multigene Alignment and Phylogeny eStimation), a user-friendly command-line and interactive package, designed to handle an improved alignment and phylogeny estimation workflow: MSF preparation, MSA estimation, outlier detection, refinement, consensus, phylogeny estimation, comparison and editing, among which file and directory organization, execution, manipulation of information are automated, with minimal manual user intervention. LMAP_S was developed for the workstation multi-core environment and provides a unique advantage for processing multiple datasets. Our software, proved to be efficient throughout the workflow, including, the (unlimited) handling of more than 20 datasets.

Conclusions

We have developed a simple and versatile LMAP_S package enabling researchers to effectively estimate multiple datasets MSAs and PTs in a high-throughput fashion. LMAP_S integrates more than 25 software providing overall more than 65 algorithm choices distributed in five stages. At minimum, one FASTA file is required within a single input directory. To our knowledge, no other software combines MSA and phylogeny estimation with as many alternatives and provides means to find optimal MSAs and phylogenies. Moreover, we used a case study comparing methodologies that highlighted the usefulness of our software. LMAP_S has been developed as an open-source package, allowing its integration into more complex open-source bioinformatics pipelines. LMAP_S package is released under GPLv3 license and is freely available at https://lmap-s.sourceforge.io/.

Fully automated sequence alignment methods are comparable to, and much faster than, traditional methods in large data sets: an example with hepatitis B virus

Aligning sequences for phylogenetic analysis (multiple sequence alignment; MSA) is an important, but increasingly computationally expensive step with the recent surge in DNA sequence data. Much of this sequence data is publicly available, but can be extremely fragmentary (i.e., a combination of full genomes and genomic fragments), which can compound the computational issues related to MSA. Traditionally, alignments are produced with automated algorithms and then checked and/or corrected "by eye" prior to phylogenetic inference. However, this manual curation is inefficient at the data scales required of modern phylogenetics and results in alignments that are not reproducible. Recently, methods have been developed for fully automating alignments of large data sets, but it is unclear if these methods produce alignments that result in compatible phylogenies when compared to more traditional alignment approaches that combined automated and manual methods. Here we use approximately 33,000 publicly available sequences from the hepatitis B virus (HBV), a globally distributed and rapidly evolving virus, to compare different alignment approaches. Using one data set comprised exclusively of whole genomes and a second that also included sequence fragments, we compared three MSA methods: (1) a purely automated approach using traditional software, (2) an automated approach including by eye manual editing, and (3) more recent fully automated approaches. To understand how these methods affect phylogenetic results, we compared resulting tree topologies based on these different alignment methods using multiple metrics. We further determined if the monophyly of existing HBV genotypes was supported in phylogenies estimated from each alignment type and under different statistical support thresholds. Traditional and fully automated alignments produced similar HBV phylogenies. Although there was variability between branch support thresholds, allowing lower support thresholds tended to result in more differences among trees. Therefore, differences between the trees could be best explained by phylogenetic uncertainty unrelated to the MSA method used. Nevertheless, automated alignment approaches did not require human intervention and were therefore considerably less time-intensive than traditional approaches. Because of this, we conclude that fully automated algorithms for MSA are fully compatible with older methods even in extremely difficult to align data sets. Additionally, we found that most HBV diagnostic genotypes did not correspond to evolutionarily-sound groups, regardless of alignment type and support threshold. This suggests there may be errors in genotype classification in the database or that HBV genotypes may need a revision.

Enabling precision medicine via standard communication of HTS provenance, analysis, and results

A personalized approach based on a patient's or pathogen’s unique genomic sequence is the foundation of precision medicine. Genomic findings must be robust and reproducible, and experimental data capture should adhere to findable, accessible, interoperable, and reusable (FAIR) guiding principles. Moreover, effective precision medicine requires standardized reporting that extends beyond wet-lab procedures to computational methods. The BioCompute framework (https://w3id.org/biocompute/1.3.0) enables standardized reporting of genomic sequence data provenance, including provenance domain, usability domain, execution domain, verification kit, and error domain. This framework facilitates communication and promotes interoperability. Bioinformatics computation instances that employ the BioCompute framework are easily relayed, repeated if needed, and compared by scientists, regulators, test developers, and clinicians. Easing the burden of performing the aforementioned tasks greatly extends the range of practical application. Large clinical trials, precision medicine, and regulatory submissions require a set of agreed upon standards that ensures efficient communication and documentation of genomic analyses. The BioCompute paradigm and the resulting BioCompute Objects (BCOs) offer that standard and are freely accessible as a GitHub organization (https://github.com/biocompute-objects) following the “Open-Stand.org principles for collaborative open standards development.” With high-throughput sequencing (HTS) studies communicated using a BCO, regulatory agencies (e.g., Food and Drug Administration [FDA]), diagnostic test developers, researchers, and clinicians can expand collaboration to drive innovation in precision medicine, potentially decreasing the time and cost associated with next-generation sequencing workflow exchange, reporting, and regulatory reviews.

This Community Page article presents a communication standard for the provenance of high-throughput sequencing data; a BioCompute Object (BCO) can serve as a history of what was computed, be used as part of a validation process, or provide clarity and transparency of an experimental process to collaborators.

Topological support and data quality can only be assessed through multiple tests in reviewing Blattodea phylogeny

Assessing support for molecular phylogenies is difficult because the data is heterogeneous in quality and overwhelming in quantity. Traditionally, node support values (bootstrap frequency, Bayesian posterior probability) are used to assess confidence in tree topologies. Other analyses to assess the quality of phylogenetic data (e.g. Lento plots, saturation plots, trait consistency) and the resulting phylogenetic trees (e.g. internode certainty, parameter permutation tests, topological tests) exist but are rarely applied. Here we argue that a single qualitative analysis is insufficient to assess support of a phylogenetic hypothesis and relate data quality to tree quality. We use six molecular markers to infer the phylogeny of Blattodea and apply various tests to assess relationship support, locus quality, and the relationship between the two. We use internode-certainty calculations in conjunction with bootstrap scores, alignment permutations, and an approximately unbiased (AU) test to assess if the molecular data unambiguously support the phylogenetic relationships found. Our results show higher support for the position of Lamproblattidae, high support for the termite phylogeny, and low support for the position of Anaplectidae, Corydioidea and phylogeny of Blaberoidea. We use Lento plots in conjunction with mutation-saturation plots, calculations of locus homoplasy to assess locus quality, identify long branch attraction, and decide if the tree's relationships are the result of data biases. We conclude that multiple tests and metrics need to be taken into account to assess tree support and data robustness.

Aligning evidence: concerns regarding multiple sequence alignments in estimating the phylogeny of the Nudibranchia suborder Doridina

Molecular estimates of phylogenetic relationships rely heavily on multiple sequence alignment construction. There has been little consensus, however, on how to properly address issues pertaining to the alignment of variable regions. Here, we construct alignments from four commonly sequenced molecular markers (16S, 18S, 28S and cytochrome c oxidase subunit I) for the Nudibranchia using three different methodologies: (i) strict mathematical algorithm; (ii) exclusion of variable or divergent regions and (iii) manually curated, and examine how different alignment construction methods can affect phylogenetic signal and phylogenetic estimates for the suborder Doridina. Phylogenetic informativeness (PI) profiles suggest that the molecular markers tested lack the power to resolve relationships at the base of the Doridina, while being more robust at family-level classifications. This supports the lack of consistent resolution between the 19 families within the Doridina across all three alignments. Most of the 19 families were recovered as monophyletic, and instances of non-monophyletic families were consistently recovered between analyses. We conclude that the alignment of variable regions has some effect on phylogenetic estimates of the Doridina, but these effects can vary depending on the size and scope of the phylogenetic query and PI of molecular markers.

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.

Phylogenetic study of Class Armophorea (Alveolata, Ciliophora) based on 18S-rDNA data

The 18S rDNA phylogeny of Class Armophorea, a group of anaerobic ciliates, is proposed based on an analysis of 44 sequences (out of 195) retrieved from the NCBI/GenBank database. Emphasis was placed on the use of two nucleotide alignment criteria that involved variation in the gap-opening and gap-extension parameters and the use of rRNA secondary structure to orientate multiple-alignment. A sensitivity analysis of 76 data sets was run to assess the effect of variations in indel parameters on tree topologies. Bayesian inference, maximum likelihood and maximum parsimony phylogenetic analyses were used to explore how different analytic frameworks influenced the resulting hypotheses. A sensitivity analysis revealed that the relationships among higher taxa of the Intramacronucleata were dependent upon how indels were determined during multiple-alignment of nucleotides. The phylogenetic analyses rejected the monophyly of the Armophorea most of the time and consistently indicated that the Metopidae and Nyctotheridae were related to the Litostomatea. There was no consensus on the placement of the Caenomorphidae, which could be a sister group of the Metopidae + Nyctorheridae, or could have diverged at the base of the Spirotrichea branch or the Intramacronucleata tree.

Parsimony and model-based analyses of indels in avian nuclear genes reveal congruent and incongruent phylogenetic signals

Insertion/deletion (indel) mutations, which are represented by gaps in multiple sequence alignments, have been used to examine phylogenetic hypotheses for some time. However, most analyses combine gap data with the nucleotide sequences in which they are embedded, probably because most phylogenetic datasets include few gap characters. Here, we report analyses of 12,030 gap characters from an alignment of avian nuclear genes using maximum parsimony (MP) and a simple maximum likelihood (ML) framework. Both trees were similar, and they exhibited almost all of the strongly supported relationships in the nucleotide tree, although neither gap tree supported many relationships that have proven difficult to recover in previous studies. Moreover, independent lines of evidence typically corroborated the nucleotide topology instead of the gap topology when they disagreed, although the number of conflicting nodes with high bootstrap support was limited. Filtering to remove short indels did not substantially reduce homoplasy or reduce conflict. Combined analyses of nucleotides and gaps resulted in the nucleotide topology, but with increased support, suggesting that gap data may prove most useful when analyzed in combination with nucleotide substitutions.

Phylogeny as a proxy for ecology in seagrass amphipods: which traits are most conserved?

Increasingly, studies of community assembly and ecosystem function combine trait data and phylogenetic relationships to gain novel insight into the ecological and evolutionary constraints on community dynamics. However, the key to interpreting these two types of information is an understanding of the extent to which traits are phylogenetically conserved. In this study, we develop the necessary framework for community phylogenetics approaches in a system of marine crustacean herbivores that play an important role in the ecosystem functioning of seagrass systems worldwide. For 16 species of amphipods and isopods, we (1) reconstructed phylogenetic relationships using COI, 16S, and 18S sequences and Bayesian analyses, (2) measured traits that are potentially important for assembling species between and within habitats, and (3) compared the degree to which each of these traits are evolutionarily conserved. Despite poor phylogenetic resolution for the order Amphipoda as a whole, we resolved almost all of the topology for the species in our system, and used a sampling of ultrametric trees from the posterior distribution to account for remaining uncertainty in topology and branch lengths. We found that traits varied widely in their degree of phylogenetic signal. Body mass, fecundity, and tube building showed very strong phylogenetic signal, and temperature tolerance and feeding traits showed much less. As such, the degree of signal was not predictable based on whether the trait is related to environmental filtering or to resource partitioning. Further, we found that even with strong phylogenetic signal in body size, (which may have large impacts on ecosystem function), the predictive relationship between phylogenetic diversity and ecosystem function is not straightforward. We show that patterns of phylogenetic diversity in communities of seagrass mesograzers could lead to a variety of interpretations and predictions, and that detailed study of trait similarities and differences will be necessary to interpret these patterns.

Stable phylogenetic patterns in scutigeromorph centipedes (Myriapoda : Chilopoda : Scutigeromorpha): dating the diversification of an ancient lineage of terrestrial arthropods

Although stable and well-supported relationships are in place for the three main clades (families) of Scutigeromorpha, the interrelationships of particular taxa within the most diverse family, Scutigeridae, are less clearly resolved. Novel molecular data for taxa from Mesoamerica, the Caribbean, southern Africa, New Guinea and previously unsampled parts of the Pacific are incorporated into phylogenetic analyses. Relationships across the tree are stable under variable analytical conditions, whether these are homology-based (multiple sequence alignment versus implied alignment; untrimmed versus trimmed datasets) or method-based (parsimony versus maximum likelihood). Hypervariable regions, contrary to common belief, add phylogenetic structure to the data, as measured by the increased support for many nodes when compared with the same alignments trimmed with Gblocks. Our analyses show that a Yule-3-rate model best explained the diversification of Scutigeromorpha during their 400 million years of history. More complete molecular data for the New Guinea genus Ballonema stabilise its position as sister group to Thereuoneminae. To reconcile scutigeromorph systematics with the phylogeny, the monotypic genus Madagassophora Verhoeff, 1936, is placed in synonymy with Scutigerina Silvestri, 1901 (n. syn.), its type species M. hova becoming Scutigerina hova (de Saussure & Zehntner, 1902) new comb. (from Scutigera), and Lassophora Verhoeff, 1905, is re-established for an Afro-Malagasy clade containing Lassophora nossibei (de Saussure & Zehntner, 1902) new comb. (from Scutigera) and a newly sequenced species from Mozambique that diverged at the base of the lineage to Thereuoneminae. The dated phylogeny of Scutigeromorpha is more consistent with ancient vicariant splits between Madagascar–southern Africa and Australia–New Caledonia than with younger dispersal scenarios, though some geologically young Pacific islands that harbour lineages dating to the Cretaceous demonstrate the potential for trans-oceanic dispersal.

A method of alignment masking for refining the phylogenetic signal of multiple sequence alignments

Inaccurate inference of positional homologies in multiple sequence alignments and systematic errors introduced by alignment heuristics obfuscate phylogenetic inference. Alignment masking, the elimination of phylogenetically uninformative or misleading sites from an alignment before phylogenetic analysis, is a common practice in phylogenetic analysis. Although masking is often done manually, automated methods are necessary to handle the much larger data sets being prepared today. In this study, we introduce the concept of subsplits and demonstrate their use in extracting phylogenetic signal from alignments. We design a clustering approach for alignment masking where each cluster contains similar columns-similarity being defined on the basis of compatible subsplits; our approach then identifies noisy clusters and eliminates them. Trees inferred from the columns in the retained clusters are found to be topologically closer to the reference trees. We test our method on numerous standard benchmarks (both synthetic and biological data sets) and compare its performance with other methods of alignment masking. We find that our method can eliminate sites more accurately than other methods, particularly on divergent data, and can improve the topologies of the inferred trees in likelihood-based analyses. Software available upon request from the author.

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.

Phylogeny of the Asparagales based on three plastid and two mitochondrial genes

PREMISE OF THE STUDY

The Asparagales, with ca. 40% of all monocotyledons, include a host of commercially important ornamentals in families such as Orchidaceae, Alliaceae, and Iridaceae, and several important crop species in genera such as Allium, Aloe, Asparagus, Crocus, and Vanilla. Though the order is well defined, the number of recognized families, their circumscription, and relationships are somewhat controversial.

METHODS

Phylogenetic analyses of Asparagales were based on parsimony and maximum likelihood using nucleotide sequence variation in three plastid genes (matK, ndhF, and rbcL) and two mitochondrial genes (atp1 and cob). Branch support was assessed using both jackknife analysis implementing strict-consensus (SC) and bootstrap analysis implementing frequency-within-replicates (FWR). The contribution of edited sites in the mitochondrial genes to topology and branch support was investigated.

KEY RESULTS

The topologies recovered largely agree with previous results, though some clades remain poorly resolved (e.g., Ruscaceae). When the edited sites were included in the analysis, the plastid and mitochondrial genes were highly incongruent. However, when the edited sites were removed, the two partitions became congruent.

CONCLUSIONS

Some deeper nodes in the Asparagales tree remain poorly resolved or unresolved as do the relationships of certain monogeneric families (e.g., Aphyllanthaceae, Ixioliriaceae, Doryanthaceae), whereas support for many families increases. However, the increased support is dominated by plastid data, and the potential influence of mitochondrial and biparentially inherited single or low-copy nuclear genes should be investigated.

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.

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.

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.

AlexSys: a knowledge-based expert system for multiple sequence alignment construction and analysis

Multiple sequence alignment (MSA) is a cornerstone of modern molecular biology and represents a unique means of investigating the patterns of conservation and diversity in complex biological systems. Many different algorithms have been developed to construct MSAs, but previous studies have shown that no single aligner consistently outperforms the rest. This has led to the development of a number of ‘meta-methods’ that systematically run several aligners and merge the output into one single solution. Although these methods generally produce more accurate alignments, they are inefficient because all the aligners need to be run first and the choice of the best solution is made a posteriori. Here, we describe the development of a new expert system, AlexSys, for the multiple alignment of protein sequences. AlexSys incorporates an intelligent inference engine to automatically select an appropriate aligner a priori, depending only on the nature of the input sequences. The inference engine was trained on a large set of reference multiple alignments, using a novel machine learning approach. Applying AlexSys to a test set of 178 alignments, we show that the expert system represents a good compromise between alignment quality and running time, making it suitable for high throughput projects. AlexSys is freely available from http://alnitak.u-strasbg.fr/∼aniba/alexsys.

The effects of fossil placement and calibration on divergence times and rates: an example from the termites (Insecta: Isoptera)

Among insects, eusocial behavior occurs in termites, ants, some bees and wasps. Isoptera and Hymenoptera convergently share social behavior, and for both taxa its evolution remains poorly understood. While dating analyses provide researchers with the opportunity to date the origin of eusociality, fossil calibration methodology may mislead subsequent ecological interpretations. Using a comprehensive termite dataset, we explored the effect of fossil placement and calibration methodology. A combined molecular and morphological dataset for 42 extant termite lineages was used, and a second dataset including these 42 taxa, plus an additional 39 fossil lineages for which we had only morphological data. MrBayes doublet-model analyses recovered similar topologies, with one minor exception (Stolotermitidae is sister to the Hodotermitidae, s.s., in the 42-taxon analysis but is in a polytomy with Hodotermitidae and (Kalotermitidae + Neoisoptera) in the 81-taxon analysis). Analyses using the r8s program on these topologies were run with either minimum/maximum constraints (analysis a = 42-taxon and analysis c = 81-taxon analyses) or with the fossil taxon ages fixed (ages fixed to be the geological age of the deposit from which they came, analysis b = 81-taxon analysis). Confidence intervals were determined for the resulting ultrametric trees, and for most major clades there was significant overlap between dates recovered for analyses A and C (with exceptions, such as the nodes Neoisoptera, and Euisoptera). With the exception of isopteran and eusiopteran node ages, however, none of the major clade ages overlapped when analysis B is compared with either analysis A or C. Future studies on Dictyoptera should note that the age of Kalotermitidae was underestimated in absence of kalotermitid fossils with fixed ages.

Arthropod phylogeny revisited, with a focus on crustacean relationships

Higher-level arthropod phylogenetics is an intensely active field of research, not least as a result of the hegemony of molecular data. However, not all areas of arthropod phylogenetics have so far received equal attention. The application of molecular data to infer a comprehensive phylogeny of Crustacea is still in its infancy, and several emerging results are conspicuously at odds with morphology-based studies. In this study, we present a series of molecular phylogenetic analyses of 88 arthropods, including 57 crustaceans, representing all the major lineages, with Onychophora and Tardigrada as outgroups. Our analyses are based on published and new sequences for two mitochondrial markers, 16S rDNA and cytochrome c oxidase subunit I (COI), and the nuclear ribosomal gene 18S rDNA. We designed our phylogenetic analyses to assess the effects of different strategies of sequence alignment, alignment masking, nucleotide coding, and model settings. Our comparisons show that alignment optimization of ribosomal markers based on secondary structure information can have a radical impact on phylogenetic reconstruction. Trees based on optimized alignments recover monophyletic Arthropoda (excluding Onychophora), Pancrustacea, Malacostraca, Insecta, Myriapoda and Chelicerata, while Maxillopoda and Hexapoda emerge as paraphyletic groups. Our results are unable to resolve the highest-level relationships within Arthropoda, and none of our trees supports the monophyly of Myriochelata or Mandibulata. We discuss our results in the context of both the methodological variations between different analyses, and of recently proposed phylogenetic hypotheses. This article offers a preliminary attempt to incorporate the large diversity of crustaceans into a single molecular phylogenetic analysis, assessing the robustness of phylogenetic relationships under varying analysis parameters. It throws into sharp relief the relative strengths and shortcomings of the combined molecular data for assessing this challenging phylogenetic problem, and thereby provides useful pointers for future studies.

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.

Molecular dating and biogeography of fig-pollinating wasps

Figs and fig-pollinating wasps are obligate mutualists that have coevolved for over 60 million years. But when and where did pollinating fig wasps (Agaonidae) originate? Some studies suggest that agaonids arose in the Late Cretaceous and the current distribution of fig-wasp faunas can be explained by the break-up of the Gondwanan landmass. However, recent molecular-dating studies suggest divergence time estimates that are inconsistent with the Gondwanan vicariance hypothesis and imply that long distance oceanic dispersal could have been an important process for explaining the current distribution of both figs and fig wasps. Here, we use a combination of phylogenetic and biogeographical data to infer the age, the major period of diversification, and the geographic origin of pollinating fig wasps. Age estimates ranged widely depending on the molecular-dating method used and even when using the same method but with slightly different constraints, making it difficult to assess with certainty a Gondwanan origin of agaonids. The reconstruction of ancestral areas suggests that the most recent common ancestor of all extant fig-pollinating wasps was most likely Asian, although a southern Gondwana origin cannot be rejected. Our analysis also suggests that dispersal has played a more important role in the development of the fig-wasp biota than previously assumed.