Phylogeny of Henicopidae (Chilopoda: Lithobiomorpha): a combined analysis of morphology and five molecular loci

Abstract. Relationships in Henicopidae, the dominant southern temperate clade of Lithobiomorpha, are appraised based on parsimony analysis of forty-nine morphological characters and sequence data from five loci (nuclear ribosomal RNAs 18S and 28S, mitochondrial ribosomal RNAs 12S and 16S, protein-coding mitochondrial cytochrome c oxidase I). A combined analysis of these data used direct character optimization, and tested stability of hypotheses through parameter-sensitivity analysis. The morphology dataset highlighted the mandibles as a source of new characters. Morphology, as well as the most congruent parameters for the sequence data and combined analysis, resolved Zygethobiini within Henicopini. Groups retrieved by combined analysis of the sequences and combination with morphology for all parameters include Anopsobiinae/Henicopinae, Lamyctes + Lamyctinus + Henicops , Paralamyctes ( Paralamyctes ), and a clade that groups the southeastern Australian/New Zealand Paralamyctes ( Haasiella ) and P . ( Thingathinga ). Paralamyctes (including Haasiella ) is a Gondwanan clade in the most congruent cladograms based on all molecular data and combination with morphological data. Biogeographic analysis of subtrees for Paralamyctes resolved the interrelationships of Gondwana as (Patagonia ((New South Wales + southeastern Queensland) ((Tasmania) (southern Africa + India) (New Zealand + north Queensland)))).     

Cladistic analysis of self‐grooming indicates a single origin of eusociality in corbiculate bees (Hymenoptera: Apidae)

Behavioural traits have been used extensively in recent years as an important character source for making phylogenetic inferences. The phylogenetic positions of the members of the Apini subtribe are increasingly being debated, and new characters must be examined. We analysed the presence and absence of certain behavioural patterns, as well as the sequences of some of these patterns, to generate 79 characters. Eleven species comprised the ingroup, and Xylocopini comprised the outgroup. Parsimony analysis showed that the most parsimonious tree was (Euglossina(Bombina(Apina+Meliponina))). This topology is consistent with most studies that use morphological data and the few that use behavioural data, which suggests that advanced eusociality arose only once in a common ancestor of the clade Apina plus Meliponina; however, this hypothesis is inconsistent with our molecular data. Thus we considered behavioural, molecular, and morphological data and recovered the same topology, in which eusociality has a single origin in corbiculate bees.     

Phylogeny of sipunculan worms: A combined analysis of four gene regions and morphology

The intra-phyletic relationships of sipunculan worms were analyzed based on DNA sequence data from four gene regions and 58 morphological characters. Initially we analyzed the data under direct optimization using parsimony as optimality criterion. An implied alignment resulting from the direct optimization analysis was subsequently utilized to perform a Bayesian analysis with mixed models for the different data partitions. For this we applied a doublet model for the stem regions of the 18S rRNA. Both analyses support monophyly of Sipuncula and most of the same clades within the phylum. The analyses differ with respect to the relationships among the major groups but whereas the deep nodes in the direct optimization analysis generally show low jackknife support, they are supported by 100% posterior probability in the Bayesian analysis. Direct optimization has been useful for handling sequences of unequal length and generating conservative phylogenetic hypotheses whereas the Bayesian analysis under mixed models provided high resolution in the basal nodes of the tree.     

A phylogenetic analysis of Coelopidae (Diptera) based on morphological and DNA sequence data

The phylogenetic relationships of 22 species of Coelopidae are reconstructed based on a data matrix consisting of morphological and DNA sequence characters (16S rDNA, EF-1alpha). Optimal gap and transversion costs are determined via a sensitivity analysis and both equal weighting and a transversion cost of 2 are found to perform best based on taxonomic congruence, character incongruence, and tree support. The preferred phylogenetic hypothesis is fully resolved and well-supported by jackknife, bootstrap, and Bremer support values, but it is in conflict with the cladogram based on morphological characters alone. Most notably, the Coelopidae and the genus Coelopa are not monophyletic. However, partitioned Bremer Support and an analysis of node stability under different gap and transversion costs reveal that the critical clades rendering these taxa non-monophyletic are poorly supported. Furthermore, the monophyly of Coelopidae and Coelopa is not rejected in analyses using 16S rDNA that was manually aligned. The resolution of the tree based on this reduced data sets is, however, lower than for the tree based on the full data sets. Partitioned Bremer support values reveal that 16S rDNA characters provide the largest amount of tree support, but the support values are heavily dependent on analysis conditions. Problems with direct comparison of branch support values for trees derived using fixed alignments with those obtained under optimization alignment are discussed. Biogeographic history and available behavioral and genetic data are also discussed in light of this first cladogram for Coelopidae based on a quantitative phylogenetic analysis.     

Phylogenetic relationships of Steinernema Travassos, 1927 (Nematoda: Cephalobina: Steinernematidae) based on nuclear, mitochondrial and morphological data

Entomopathogenic nematodes of the genus Steinernema are lethal parasites of insects that are used as biological control agents of several lepidopteran, dipteran and coleopteran pests. Phylogenetic relationships among 25 Steinernema species were estimated using nucleotide sequences from three genes and 22 morphological characters. Parsimony analysis of 28S (LSU) sequences yielded a well-resolved phylogenetic hypothesis with reliable bootstrap support for 13 clades. Parsimony analysis of mitochondrial DNA sequences (12S rDNA and cox 1 genes) yielded phylogenetic trees with a lower consistency index than for LSU sequences, and with fewer reliably supported clades. Combined phylogenetic analysis of the 3-gene dataset by parsimony and Bayesian methods yielded well-resolved and highly similar trees. Bayesian posterior probabilities were high for most clades; bootstrap (parsimony) support was reliable for approximately half of the internal nodes. Parsimony analysis of the morphological dataset yielded a poorly resolved tree, whereas total evidence analysis (molecular plus morphological data) yielded a phylogenetic hypothesis consistent with, but less resolved than trees inferred from combined molecular data. Parsimony mapping of morphological characters on the 3-gene trees showed that most structural features of steinernematids are highly homoplastic. The distribution of nematode foraging strategies on these trees predicts that S. hermaphroditum, S. diaprepesi and S. longicaudum (US isolate) have cruise forager behaviours.     

Phylogeny,biogeography and the stepwise evolutionary colonization of intertidal habitat in the Liparocephalini based on morphological and molecular characters (Coleoptera: Staphylinidae: Aleocharinae)

A phylogenetic analysis of the tribe Liparocephalini Fenyes is presented based on morphological and molecular characters. The data set comprised 50 adult morphological characters, partial COI (907 bp), COII (366 bp) and 12S rDNA (325–355 bp), and nearly complete sequences of 18S rDNA (1768–1902 bp) for 21 species. Eighteen species of liparocephaline beetles from all eight genera and three outgroups, are included. The sequences were analysed separately and simultaneously with morphological characters by direct optimization in the program POY4 and by partitioned Bayesian analysis for the combined data. The direct optimization (DO) tree for the combined data under equal weighting, which also shows a minimum incongruence length difference value, resulted in a monophyletic Liparocephalini with the following patterns of phylogenetic relationships (outgroup ((Baeostethus, Ianmoorea) (Paramblopusa ((Amblopusa, Halorhadinus) (Liparocephalus, Diaulota))))). A sensitivity analysis using 16 different parameter sets for the combined data shows the monophyly of the liparocephalines and all its genera under all parameter sets. Bayesian analysis resulted in topological differences in comparison with the DO tree under equal weighting only in the position of the genus Paramblopusa and clade (Amblopusa + Halorhadinus), which were reversed. Historical biogeography and the stepwise evolutionary colonization of intertidal habitat in the Liparocephalini are discussed. Based on the biogeographical analyses, we hypothesize that the ancestor of the Liparocephalini occurred along the Panthallassan Ocean, the direct antecedent of the Pacific Ocean, followed by repeated dispersals to the Nearctic from the Palearctic. We also hypothesize that ancestors of the Liparocephalini appear to have arisen in the littoral zone of beaches and then colonized rocky reef areas in the low tidal zone later through high‐ to mid‐tide zones. © The Willi Hennig Society 2009.     

Adding mitochondrial sequence data (16S rRNA and cytochrome c oxidase subunit I) to the phylogeny of centipedes (Myriapoda: Chilopoda): an analysis of morphology and four molecular loci

Relationships within Chilopoda (centipedes) are assessed based on 222 morphological characters, complete 18S rRNA sequences for 70 chilopod terminals, the D3 region of 28S rRNA for 65 terminals, 16S rRNA sequences for 54 terminals and cytochrome c oxidase subunit I sequences for 45 terminals. Morphological and molecular data for seven orders of Diplopoda are used to root cladograms for Chilopoda. Analyses use direct character optimization for 15 gap and substitution models. The Pleurostigmophora and Epimorpha s.l. hypotheses are largely stable to parameter variation for the combined data; the latter clade is formalized as the new taxon Phylactometria. The combined data include parameter sets that support either the monophyly of Epimorpha s.str. (=Scolopendromorpha + Geophilomorpha) or Craterostigmus + Geophilomorpha; the former derives its support from morphology and the nuclear ribosomal genes. Monophyly of Lithobiomorpha and the sister group relationship between Lithobiidae and Henicopidae are stable for morphological and combined data, and are also resolved for the molecular data for 14 of 15 parameter sets. The fundamental split in Scolopendromorpha is between Cryptopidae and Scolopendridae sensu Attems. Blind scolopendromorphs unite as a clade in most molecular and combined analyses, including those that minimize incongruence between data partitions. Geophilomorpha divides into Placodesmata and Adesmata under nine of 15 explored parameter sets.     

Phylogenetic relationships of Hamadryas (Nymphalidae: Biblidinae) based on the combined analysis of morphological and molecular data

A new phylogenetic hypothesis for the Neotropical butterfly genus Hamadryas, based on the combination of a morphological matrix, one mitochondrial (COI) and four nuclear markers (CAD, RpS5, EF1a, and Wingless), is presented. Results from analyses of the molecular evidence are compared with a previously published morphological phylogeny. Molecular data and the analysis of the complete dataset support the monophyly of Hamadryas and most sister groups suggested by morphological data alone. The addition of DNA sequences to the morphological matrix helped define species groups for which no morphological synapomorphies were found. Partitioned Bremer support indicates that COI, CAD, and morphology were consistently in agreement with the combined evidence tree. In contrast, signal from the nuclear markers Rps5, EF1a, and Wingless showed indifference at most levels of the tree, and minor conflict at nodes solving the relationships between species groups. Though resolved, the combined evidence tree shows low resample values, particularly among species groups whose relationships were characterized by short internodes. A reassessment about the pattern of character change for sound production is presented and discussed.     

Phylogeny of Glaucosomatidae inferred from molecular evidence

Most species of glaucosomatids (Teleostei: Glaucosomatidae) are endemic to Australia, except Glaucosoma buergeri that is widely distributed from Australia to Japan. This study elucidated phylogenetic relationships among glaucosomatids based on the morphological characters of the saccular‐otolith sagitta, in addition to molecular evidence of mitochondrial 16S rDNA, cytochrome oxidase I (COI) and cytochrome b (cyt b) sequences, and nuclear rhodopsin sequences. The topologies of individuals' phylogenetic trees, based on 16S rDNA, COI and cyt b sequences, were statistically indistinguishable from one another, and were only slightly different from a tree based on rhodopsin sequences. These molecular tree topologies, however, differed from species relationships in morphology‐based phylogenetic hypothesis proposed in previous studies. Specimens of G. buergeri from Australia and Taiwan showed differences in the sagitta and molecular differentiation at the four genes, suggesting a possible speciation event. Both molecular and morphological evidences indicate that Glaucosoma magnificum is the plesiomorphic sister species of other glaucosomatid species. Glaucosoma hebraicum is the sister species of a clade composed of G. buergeri and Glaucosoma scapulare. Molecular and morphological evidences also support the species status of G. hebraicum.     

Phylogenetic relationships of Steinernema Travassos, 1927 (Nematoda: Cephalobina: Steinernematidae) based on nuclear, mitochondrial and morphological data

Entomopathogenic nematodes of the genus Steinernema are lethal parasites of insects that are used as biological control agents of several lepidopteran, dipteran and coleopteran pests. Phylogenetic relationships among 25 Steinernema species were estimated using nucleotide sequences from three genes and 22 morphological characters. Parsimony analysis of 28S (LSU) sequences yielded a well-resolved phylogenetic hypothesis with reliable bootstrap support for 13 clades. Parsimony analysis of mitochondrial DNA sequences (12S rDNA and cox 1 genes) yielded phylogenetic trees with a lower consistency index than for LSU sequences, and with fewer reliably supported clades. Combined phylogenetic analysis of the 3-gene dataset by parsimony and Bayesian methods yielded well-resolved and highly similar trees. Bayesian posterior probabilities were high for most clades; bootstrap (parsimony) support was reliable for approximately half of the internal nodes. Parsimony analysis of the morphological dataset yielded a poorly resolved tree, whereas total evidence analysis (molecular plus morphological data) yielded a phylogenetic hypothesis consistent with, but less resolved than trees inferred from combined molecular data. Parsimony mapping of morphological characters on the 3-gene trees showed that most structural features of steinernematids are highly homoplastic. The distribution of nematode foraging strategies on these trees predicts that S. hermaphroditum, S. diaprepesi and S. longicaudum (US isolate) have cruise forager behaviours.     

Coding behavioural data for cladistic analysis: using dynamic homology without parsimony

Many of the controversies around the concept of homology rest on the subjectivity inherent to primary homology propositions. Dynamic homology partially solves this problem, but there has been up to now scant application of it outside of the molecular domain. This is probably because morphological and behavioural characters are rich in properties, connections and qualities, so that there is less space for conflicting character delimitations. Here we present a new method for the direct optimization of behavioural data, a method that relies on the richness of this database to delimit the characters, and on dynamic procedures to establish character state identity. We use between‐species congruence in the data matrix and topological stability to choose the best cladogram. We test the methodology using sequences of predatory behaviour in a group of spiders that evolved the highly modified predatory technique of spitting glue onto prey. The cladogram recovered is fully compatible with previous analyses in the literature, and thus the method seems consistent. Besides the advantage of enhanced objectivity in character proposition, the new procedure allows the use of complex, context‐dependent behavioural characters in an evolutionary framework, an important step towards the practical integration of the evolutionary and ecological perspectives on diversity. © The Willi Hennig Society 2010.     

Morphology,molecules, and the phylogenetics of cetaceans

Recent phylogenetic analyses of cetacean relationships based on DNA sequence data have challenged the traditional view that baleen whales (Mysticeti) and toothed whales (Odontoceti) are each monophyletic, arguing instead that baleen whales are the sister group of the odontocete family Physeteridae (sperm whales). We reexamined this issue in light of a morphological data set composed of 207 characters and molecular data sets of published 12S, 16S, and cytochrome b mitochondrial DNA sequences. We reach four primary conclusions: (1) Our morphological data set strongly supports the traditional view of odontocete monophyly; (2) the unrooted molecular and morphological trees are very similar, and most of the conflict results from alternative rooting positions; (3) the rooting position of the molecular tree is sensitive to choice of artiodactyls outgroup taxa and the treatment of two small but ambiguously aligned regions of the 12S and 16S sequences, whereas the morphological root is strongly supported; and (4) combined analyses of the morphological and molecular data provide a well-supported phylogenetic estimate consistent with that based on the morphological data alone (and the traditional view of toothed-whale monophyly) but with increased bootstrap support at nearly every node of the tree.     

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

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

Towards an 18S phylogeny of hexapods: accounting for group-specific character covariance in optimized mixed nucleotide/doublet models

The phylogenetic diversification of Hexapoda is still not fully understood. Morphological and molecular analyses have resulted in partly contradicting hypotheses. In molecular analyses, 18S sequences are the most frequently employed, but it appears that 18S sequences do not contain enough phylogenetic signals to resolve basal relationships of hexapod lineages. Until recently, character interdependence in these data has never been treated seriously, though possibly accounting for the occurrence of biased results. However, software packages are readily available which can incorporate information on character interdependence within a Bayesian approach. Accounting for character covariation derived from a hexapod consensus secondary structure model and applying mixed DNA/RNA substitution models, our Bayesian analysis of 321 hexapod sequences yielded a partly robust tree that depicts many hexapod relationships congruent with morphological considerations. It appears that the application of mixed DNA/RNA models removes many of the anomalies seen in previous studies. We focus on basal hexapod relationships for which unambiguous results are missing. In particular, the strong support for a “Chiastomyaria” clade (Ephemeroptera+Neoptera) obtained in Kjer's [2004. Aligned 18S and insect phylogeny. Syst. Biol. 53, 1–9] study of 18S sequences could not be confirmed by our analysis. The hexapod tree can be rooted with monophyletic Entognatha but not with a clade Ellipura (Collembola+Protura). Compared to previously published contributions, accounting for character interdependence in analyses of rRNA data presents an improvement of phylogenetic resolution. We suggest that an integration of explicit clade-specific rRNA structural refinements is not only possible but an important step in the optimization of substitution models dealing with rRNA data.     

Phylogeny of Steinernema travassos, 1927 (Cephalobina: Steinernematidae) inferred from ribosomal DNA sequences and morphological characters

Entomopathogenic nematodes in Steinernema, together with their symbiont bacteria Xenorhabdus, are obligate and lethal parasites of insects that can provide effective biological control of some important lepidopteran, dipteran, and coleopteran pests of commercial crops. Phylogenetic relationships among 21 Steinernema species were estimated using 28S ribosomal DNA (rDNA) sequences and morphological characters. Sequences of the rDNA internal transcribed spacers were obtained to provide additional molecular characters to resolve relationships among Steinernema carpocapsae, Steinernema scapterisci, Steinernema siamkavai, and Steinernema monticolum. Four equally parsimonious trees resulted from combined analysis of 28S sequences and 22 morphological characters. Clades inferred from analyses of molecular sequences and combined datasets were primarily reliably supported as assessed by bootstrap resampling, whereas those inferred from morphological data alone were not. Although partially consistent with some traditional expectations and previous phylogenetic studies, the hypotheses inferred from molecular evidence, and those from combined analysis of morphological and molecular data, provide a new and comprehensive framework for evaluating character evolution of steinernematids. Interpretation of morphological character evolution on 6 trees inferred from sequence data and combined evidence suggests that many structural features of these nematodes are highly homoplastic, and that some structures previously used to hypothesize relationships represent ancestral character states.     

Phylogeny of Tunicata inferred from molecular and morphological characters

The phylogeny of the Tunicata was reconstructed using molecular and morphological characters. Mitochondrial cytochrome oxidase I (cox1) and 18S rDNA sequences were obtained for 14 and 4 tunicate species, respectively. 18S rDNA sequences were aligned with gene sequences obtained from GenBank of 57 tunicates, a cephalochordate, and a selachian craniate. Cox1 sequences were aligned with the sequence of two ascidians and a cephalochordate obtained from GenBank. Traditional, morphological, life history, and biochemical characters of larval and adult stages were compiled from the literature and analyzed cladistically. Separate and simultaneous parsimony analyses of molecular and morphological data were carried out. Aplousobranch ascidians from three different families were included in a molecular phylogenetic analysis for the first time. Analysis of the morphological, life history, and biochemical characters results in a highly unresolved tree. Aplousobranchiata form a strongly supported monophylum in the analysis of the 18S rDNA data, the morphological data, and the combined data set. Cionidae is not included in the Aplousobranchiata but nests within the Phlebobranchiata. Appendicularia (=Larvacea) nest within the 'Ascidiacea' as the sister taxon of Aplousobranchiata in the parsimony analysis of the 18S rDNA data and the combined analysis. A potential morphological synapomorphy of Aplousobranchiata plus Appendicularia is the horizontal orientation of the larval tail. In the 18S rDNA and the combined analysis, Thaliacea is included in the 'Ascidiacea' as the sister group to Phlebobranchiata. Pyrosomatida is found to be the sister taxon to the Salpidae in analyses of 18S rDNA and combined data, whereas the analysis of the morphological data recovers a sister group relationship between Doliolidae and Salpidae. Results of cox1 analyses are incongruent with both the 18S rDNA and the morphological phylogenies. Cox1 sequences may evolve too rapidly to resolve relationships of higher tunicate taxa. However, the cox1 data may be useful at lower taxonomic levels.     

From kissing to belly stridulation: comparative analysis reveals surprising diversity,rapid evolution,and much homoplasy in the mating behaviour of 27 species of sepsid flies (Diptera: Sepsidae)

Our understanding of how fast mating behaviour evolves in insects is rather poor due to a lack of comparative studies among insect groups for which phylogenetic relationships are known. Here, we present a detailed study of the mating behaviour of 27 species of Sepsidae (Diptera) for which a well‐resolved and supported phylogeny is available. We demonstrate that mating behaviour is extremely diverse in sepsids with each species having its own mating profile. We define 32 behavioural characters and document them with video clips. Based on sister species comparisons, we provide several examples where mating behaviour evolves faster than all sexually dimorphic morphological traits. Mapping the behaviours onto the molecular tree reveals much homoplasy, comparable to that observed for third positions of mitochondrial protein‐encoding genes. A partitioned Bremer support (PBS) analysis reveals conflict between the molecular and behavioural data, but behavioural characters have higher PBS values per parsimony‐informative character than DNA sequence characters.     

Evolution and phylogenetic utility of alignment gaps within intron sequences of three nuclear genes in bumble bees (Bombus)

To test whether gaps resulting from sequence alignment contain phylogenetic signal concordant with those of base substitutions, we analyzed the occurrence of indel mutations upon a well-resolved, substitution-based tree for three nuclear genes in bumble bees (Bombus, Apidae: Bombini). The regions analyzed were exon and intron sequences of long-wavelength rhodopsin (LW Rh), arginine kinase (ArgK), and elongation factor-1alpha (EF-1alpha) F2 copy genes. LW Rh intron had only a few uninformative gaps, ArgK intron had relatively long gaps that were easily aligned, and EF-1alpha intron had many short gaps, resulting in multiple optimal alignments. The unambiguously aligned gaps within ArgK intron sequences showed no homoplasy upon the substitution-based tree, and phylogenetic signals within ambiguously aligned regions of EF-1alpha intron were highly congruent with those of base substitutions. We further analyzed the contribution of gap characters to phylogenetic reconstruction by incorporating them in parsimony analysis. Inclusion of gap characters consistently improved support for nodes recovered by substitutions, and inclusion of ambiguously aligned regions of EF-1alpha intron resolved several additional nodes, most of which were apical on the phylogeny. We conclude that gaps are an exceptionally reliable source of phylogenetic information that can be used to corroborate and refine phylogenies hypothesized by base substitutions, at least at lower taxonomic levels. At present, full use of gaps in phylogenetic reconstruction is best achieved in parsimony analysis, pending development of well-justified and generally applicable methods for incorporating indels in explicitly model-based methods.