Aligned 18S and insect phylogeny

The nuclear small subunit rRNA (18S) has played a dominant role in the estimation of relationships among insect orders from molecular data. In previous studies, 18S sequences have been aligned by unadjusted automated approaches (computer alignments that are not manually readjusted), most recently with direct optimization (simultaneous alignment and tree building using a program called "POY"). Parsimony has been the principal optimality criterion. Given the problems associated with the alignment of rRNA, and the recent availability of the doublet model for the analysis of covarying sites using Bayesian MCMC analysis, a different approach is called for in the analysis of these data. In this paper, nucleotide sequence data from the 18S small subunit rRNA gene of insects are aligned manually with reference to secondary structure, and analyzed under Bayesian phylogenetic methods with both GTR+I+G and doublet models in MrBayes. A credible phylogeny of Insecta is recovered that is independent of the morphological data and (unlike many other analyses of 18S in insects) not contradictory to traditional ideas of insect ordinal relationships based on morphology. Hexapoda, including Collembola, are monophyletic. Paraneoptera are the sister taxon to a monophyletic Holometabola but weakly supported. Ephemeroptera are supported as the sister taxon of Neoptera, and this result is interpreted with respect to the evolution of direct sperm transfer and the evolution of flight. Many other relationships are well-supported but several taxa remain problematic, e.g., there is virtually no support for relationships among orthopteroid orders. A website is made available that provides aligned 18S data in formats that include structural symbols and Nexus formats.     

Multilocus ribosomal RNA phylogeny of the leaf beetles (Chrysomelidae)

Basal relationships in the Chrysomelidae (leaf beetles) were investigated using two nuclear (small and partial large subunits) and mitochondrial (partial large subunit) rRNA (≈ 3000 bp total) for 167 taxa covering most major lineages and relevant outgroups. Separate and combined data analyses were performed under parsimony and model‐based tree building algorithms from dynamic (direct optimization) and static (Clustal and BLAST) sequence alignments. The performance of methods differed widely and recovery of well established nodes was erratic, in particular when using single gene partitions, but showed a slight advantage for Bayesian inferences and one of the fast likelihood algorithms (PHYML) over others. Direct optimization greatly gained from simultaneous analysis and provided a valuable hypothesis of chrysomelid relationships. The BLAST‐based alignment, which removes poorly aligned sequence segments, in combination with likelihood and Bayesian analyses, resulted in highly defensible trees obtained in much shorter time than direct optimization, and hence is a viable alternative when data sets grow. The main taxonomic findings include the recognition of three major lineages of Chrysomelidae, including a basal “sagrine” clade (Criocerinae, Donaciinae, Bruchinae), which was sister to the “eumolpine” (Spilopyrinae, Eumolpinae, Cryptocephalinae, Cassidinae) plus “chrysomeline” (Chrysomelinae, Galerucinae) clades. The analyses support a broad definition of subfamilies (i.e., merging previously separated subfamilies) in the case of Cassidinae (cassidines + hispines) and Cryptocephalinae (chlamisines + cryptocephalines + clytrines), whereas two subfamilies, Chrysomelinae and Eumolpinae, were paraphyletic. The surprising separation of monocot feeding Cassidinae (associated with the eumolpine clade) from the other major monocot feeding groups in the sagrine clade was well supported. The study highlights the need for thorough taxon sampling, and reveals that morphological data affected by convergence had a great impact when combined with molecular data in previous phylogenetic analyses of Chrysomelidae. © The Willi Hennig Society 2007.     

Phylogenetics of barbets (Aves: Piciformes) based on nuclear and mitochondrial DNA sequence data

A combination of nuclear (beta-fibrinogen intron 7; 938 bp) and mitochondrial (cytochrome b; 1045 bp) DNA sequence data analyzed with model-based phylogenetic methods yields a hypothesis of barbet and toucan relationships supported by high Bayesian posterior probabilities and several synapomorphic indels in the nuclear intron data. The basal topology differs from previous morphology and mitochondrial DNA sequence based hypotheses, but is consistent with DNA-DNA hybridization results. The mitochondrial DNA sequence data provide resolution at the tips of the tree, but when analyzed alone, yield a different basal topology than the combined data. However, the basal nodes in the mitochondrial-based tree have little statistical support. Combined data analysis produced strong statistical support at basal nodes and a relatively simple geographic structure. Barbets from each of the three tropical regions are monophyletic, but the Old World barbets are paraphyletic. The African and New World clades are sister taxa, and the Asian clade is basal. This geographic structure indicates that similarities in plumage, voice, and behavior between Gymnobucco and Calorhamphus are convergent. The data are not conclusive, but suggest paraphyly of the New World barbets. Only 4% of the Bayesian posterior distribution unites Semnornis with the other New World barbets.     

A phylogeny of anisopterous dragonflies (Insecta,Odonata) using mtRNA genes and mixed nucleotide/doublet models

The application of mixed nucleotide/doublet substitution models has recently received attention in RNA‐based phylogenetics. Within a Bayesian approach, it was shown that mixed models outperformed analyses relying on simple nucleotide models. We analysed an mt RNA data set of dragonflies representing all major lineages of Anisoptera plus outgroups, using a mixed model in a Bayesian and parsimony (MP) approach. We used a published mt 16S rRNA secondary consensus structure model and inferred consensus models for the mt 12S rRNA and tRNA valine. Secondary structure information was used to set data partitions for paired and unpaired sites on which doublet or nucleotide models were applied, respectively. Several different doublet models are currently available of which we chose the most appropriate one by a Bayes factor test. The MP reconstructions relied on recoded data for paired sites in order to account for character covariance and an application of the ratchet strategy to find most parsimonious trees. Bayesian and parsimony reconstructions are partly differently resolved, indicating sensitivity of the reconstructions to model specification. Our analyses depict a tree in which the damselfly family Lestidae is sister group to a monophyletic clade Epiophlebia + Anisoptera, contradicting recent morphological and molecular work. In Bayesian analyses, we found a deep split between Libelluloidea and a clade ‘Aeshnoidea’ within Anisoptera largely congruent with Tillyard’s early ideas of anisopteran evolution, which had been based on evidently plesiomorphic character states. However, parsimony analysis did not support a clade ‘Aeshnoidea’, but instead, placed Gomphidae as sister taxon to Libelluloidea. Monophyly of Libelluloidea is only modestly supported, and many inter‐family relationships within Libelluloidea do not receive substantial support in Bayesian and parsimony analyses. We checked whether high Bayesian node support was inflated owing to either: (i) wrong secondary consensus structures; (ii) under‐sampling of the MCMC process, thereby missing other local maxima; or (iii) unrealistic prior assumptions on topologies or branch lengths. We found that different consensus structure models exert strong influence on the reconstruction, which demonstrates the importance of taxon‐specific realistic secondary structure models in RNA phylogenetics.     

Modeling nucleotide evolution at the mesoscale: the phylogeny of the neotropical pitvipers of the Porthidium group (viperidae: crotalinae)

We analyzed the phylogeny of the Neotropical pitvipers within the Porthidium group (including intra-specific through inter-generic relationships) using 1.4 kb of DNA sequences from two mitochondrial protein-coding genes (ND4 and cyt-b). We investigated how Bayesian Markov chain Monte-Carlo (MCMC) phylogenetic hypotheses based on this 'mesoscale' dataset were affected by analysis under various complex models of nucleotide evolution that partition models across the dataset. We develop an approach, employing three statistics (Akaike weights, Bayes factors, and relative Bayes factors), for examining the performance of complex models in order to identify the best-fit model for data analysis. Our results suggest that: (1) model choice may have important practical effects on phylogenetic conclusions even for mesoscale datasets, (2) the use of a complex partitioned model did not produce widespread increases or decreases in nodal posterior probability support, and (3) most differences in resolution resulting from model choice were concentrated at deeper nodes. Our phylogenetic estimates of relationships among members of the Porthidium group (genera: Atropoides, Cerrophidion, and Porthidium) resolve the monophyly of the three genera. Bayesian MCMC results suggest that Cerrophidion and Porthidium form a clade that is the sister taxon to Atropoides. In addition to resolving the intra-specific relationships among a majority of Porthidium group taxa, our results highlight phylogeographic patterns across Middle and South America and suggest that each of the three genera may harbor undescribed species diversity.     

基于稳定氧同位素确定植物水分来源不同方法的比较

利用稳定同位素技术确定植物水分来源,对提高生态水文过程的认识和对干旱半干旱区的生态管理至关重要。目前基于稳定同位素技术确定植物水分来源的方法众多,但不同方法之间对比的研究较少。本研究基于原位样品采集,室内实验测试,利用直接对比法、多元线性混合模型(IsoSource)、贝叶斯混合模型(MixSIR、MixSIAR)和吸水深度模型分析植物水分来源,并对比各方法的优缺点。结果表明:相对于多元线性混合模型(IsoSource)而言,贝叶斯混合模型(MixSIR、MixSIAR)具有更好的水源区分性能,但对数据要求较高,且植物木质部水和潜在水源同位素组成的标准差越小,模型运行结果的可信度更高。本研究中贝叶斯混合模型(MixSIR)为最优解。在利用稳定氢氧同位素技术确定植物水分来源时,可先通过直接对比法定性判断植物可能利用的潜在水源,然后再用多元线性混合模型(IsoSource)、贝叶斯混合模型(MixSIR、MixSIAR)计算出各潜在水源对植物的贡献率和贡献范围,必要时可评估模型性能,选择出最优模型,定量分析植物的水分来源。若植物主要吸收利用不同土层深度的土壤水,可结合吸水深度模型计算出植物吸收土壤水的平均深度。本研究为干旱半干旱地区利用同位素技术确定植物水分来源方法的选择提供了理论依据。     

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.     

The devil in the details: interactions between the branch-length prior and likelihood model affect node support and branch lengths in the phylogeny of the Psoraceae

In popular use of Bayesian phylogenetics, a default branch-length prior is almost universally applied without knowing how a different prior would have affected the outcome. We performed Bayesian and maximum likelihood (ML) inference of phylogeny based on empirical nucleotide sequence data from a family of lichenized ascomycetes, the Psoraceae, the morphological delimitation of which has been controversial. We specifically assessed the influence of the combination of Bayesian branch-length prior and likelihood model on the properties of the Markov chain Monte Carlo tree sample, including node support, branch lengths, and taxon stability. Data included two regions of the mitochondrial ribosomal RNA gene, the internal transcribed spacer region of the nuclear ribosomal RNA gene, and the protein-coding largest subunit of RNA polymerase II. Data partitioning was performed using Bayes' factors, whereas the best-fitting model of each partition was selected using the Bayesian information criterion (BIC). Given the data and model, short Bayesian branch-length priors generate higher numbers of strongly supported nodes as well as short and topologically similar trees sampled from parts of tree space that are largely unexplored by the ML bootstrap. Long branch-length priors generate fewer strongly supported nodes and longer and more dissimilar trees that are sampled mostly from inside the range of tree space sampled by the ML bootstrap. Priors near the ML distribution of branch lengths generate the best marginal likelihood and the highest frequency of "rogue" (unstable) taxa. The branch-length prior was shown to interact with the likelihood model. Trees inferred under complex partitioned models are more affected by the stretching effect of the branch-length prior. Fewer nodes are strongly supported under a complex model given the same branch-length prior. Irrespective of model, internal branches make up a larger proportion of total tree length under the shortest branch-length priors compared with longer priors. Relative effects on branch lengths caused by the branch-length prior can be problematic to downstream phylogenetic comparative methods making use of the branch lengths. Furthermore, given the same branch-length prior, trees are on average more dissimilar under a simple unpartitioned model compared with a more complex partitioned models. The distribution of ML branch lengths was shown to better fit a gamma or Pareto distribution than an exponential one. Model adequacy tests indicate that the best-fitting model selected by the BIC is insufficient for describing data patterns in 5 of 8 partitions. More general substitution models are required to explain the data in three of these partitions, one of which also requires nonstationarity. The two mitochondrial ribosomal RNA gene partitions need heterotachous models. We found no significant correlations between, on the one hand, the amount of ambiguous data or the smallest branch-length distance to another taxon and, on the other hand, the topological stability of individual taxa. Integrating over several exponentially distributed means under the best-fitting model, node support for the family Psoraceae, including Psora, Protoblastenia, and the Micarea sylvicola group, is approximately 0.96. Support for the genus Psora is distinctly lower, but we found no evidence to contradict the current classification.     

Phylogenetic relationships of the Dactyloa clade of Anolis lizards based on nuclear and mitochondrial DNA sequence data

The Dactyloa clade, one of two major subgroups of mainland Anolis lizards, is distributed from Costa Rica to Peru, including the Amazon region and the southern Lesser Antilles. We estimated the phylogenetic relationships within Dactyloa based on mitochondrial (ND2, five transfer-RNAs, COI) and nuclear (RAG1) gene regions using likelihood and Bayesian methods under different partition strategies. In addition, we tested the monophyly of five previously recognized groups within Dactyloa. The data strongly support the monophyly of Dactyloa and five major clades: eastern, latifrons, Phenacosaurus, roquet and western, each of which exhibits a coherent geographic range. Relationships among the five major clades are less clear: support for basal nodes within Dactyloa is weak and some contradictory relationships are supported by different datasets and/or phylogenetic methods. Of the previously recognized subgroups within Dactyloa, only the roquet series consistently passed the topology tests applied. The monophyly of the aequatorialis, latifrons (as traditionally circumscribed) and punctatus series was strongly rejected, and the monophyly of Phenacosaurus (as traditionally circumscribed) yielded mixed results. The results of the phylogenetic analyses suggest the need for a revised taxonomy and have implications for the biogeography and tempo of the Dactyloa radiation.     

First molecular data on the phylum Loricifera: an investigation into the phylogeny of ecdysozoa with emphasis on the positions of Loricifera and Priapulida

Recent progress in molecular techniques has generated a wealth of information for phylogenetic analysis. Among metazoans all but a single phylum have been incorporated into some sort of molecular analysis. However, the minute and rare species of the phylum Loricifera have remained elusive to molecular systematists. Here we report the first molecular sequence data (nearly complete 18S rRNA) for a member of the phylum Loricifera, Pliciloricus sp. from Korea. The new sequence data were analyzed together with 52 other ecdysozoan sequences, with all other phyla represented by three or more sequences. The data set was analyzed using parsimony as an optimality criterion under direct optimization as well as using a Bayesian approach. The parsimony analysis was also accompanied by a sensitivity analysis. The results of both analyses are largely congruent, finding monophyly of each ecdysozoan phylum, except for Priapulida, in which the coelomate Meiopriapulus is separate from a clade of pseudocoelomate priapulids. The data also suggest a relationship of the pseudocoelomate priapulids to kinorhynchs, and a relationship of nematodes to tardigrades. The Bayesian analysis placed the arthropods as the sister group to a clade that includes tardigrades and nematodes. However, these results were shown to be parameter dependent in the sensitivity analysis. The position of Loricifera was extremely unstable to parameter variation, and support for a relationship of loriciferans to any particular ecdysozoan phylum was not found in the data.     

Bayesian phylogenetic analysis of combined data

The recent development of Bayesian phylogenetic inference using Markov chain Monte Carlo (MCMC) techniques has facilitated the exploration of parameter-rich evolutionary models. At the same time, stochastic models have become more realistic (and complex) and have been extended to new types of data, such as morphology. Based on this foundation, we developed a Bayesian MCMC approach to the analysis of combined data sets and explored its utility in inferring relationships among gall wasps based on data from morphology and four genes (nuclear and mitochondrial, ribosomal and protein coding). Examined models range in complexity from those recognizing only a morphological and a molecular partition to those having complex substitution models with independent parameters for each gene. Bayesian MCMC analysis deals efficiently with complex models: convergence occurs faster and more predictably for complex models, mixing is adequate for all parameters even under very complex models, and the parameter update cycle is virtually unaffected by model partitioning across sites. Morphology contributed only 5% of the characters in the data set but nevertheless influenced the combined-data tree, supporting the utility of morphological data in multigene analyses. We used Bayesian criteria (Bayes factors) to show that process heterogeneity across data partitions is a significant model component, although not as important as among-site rate variation. More complex evolutionary models are associated with more topological uncertainty and less conflict between morphology and molecules. Bayes factors sometimes favor simpler models over considerably more parameter-rich models, but the best model overall is also the most complex and Bayes factors do not support exclusion of apparently weak parameters from this model. Thus, Bayes factors appear to be useful for selecting among complex models, but it is still unclear whether their use strikes a reasonable balance between model complexity and error in parameter estimates.     

First combined cladistic analysis of marsupial mammal interrelationships

We combine osteological, dental, and soft tissue data with sequences from three nuclear and five mitochondrial genes, sampling all major living clades of marsupials plus several extinct taxa, to perform a simultaneous analysis of marsupial interrelationships. These data were analyzed using direct optimization and sensitivity analysis on a parallel supercomputing cluster, and compared with trees produced with conventional parsimony and likelihood algorithms using a static alignment. A major issue in marsupial phylogeny is the relationships among australidelphians. Optimal analyses using direct optimization and those based on the static alignment support the basal positions of peramelians (bandicoots) and Dromiciops ('monito del monte') within Australidelphia, and in all but one case these analyses support a monophyletic Eometatheria, a group consisting of all australidelphians excluding peramelians. Dromiciops is basal within Eometatheria in analyses that maximize congruence across partitions, including the equally weighted parameter set. The topologies resulting from direct optimization under all parameter sets show some differences, but all show a high degree of resolution. Direct optimization supports high-level clades supported by analyses of partitioned molecular (e.g., Notoryctes as sister group of Dasyuromorphia) and morphological (e.g., Diprotodontia) data.     

Phylogenetic relationships and heterogeneous evolutionary processes among phrynosomatine sand lizards (Squamata, Iguanidae) revisited

Phylogenetic analyses of DNA sequences were conducted to evaluate four alternative hypotheses of phrynosomatine sand lizard relationships. Sequences comprising 2871 aligned base pair positions representing the regions spanning ND1-COI and cyt b-tRNA(Thr) of the mitochondrial genome from all recognized sand lizard species were analyzed using unpartitioned parsimony and likelihood methods, likelihood methods with assumed partitions, Bayesian methods with assumed partitions, and Bayesian mixture models. The topology (Uma, (Callisaurus, (Cophosaurus, Holbrookia))) and thus monophyly of the "earless" taxa, Cophosaurus and Holbrookia, is supported by all analyses. Previously proposed topologies in which Uma and Callisaurus are sister taxa and those in which Holbrookia is the sister group to all other sand lizard taxa are rejected using both parsimony and likelihood-based significance tests with the combined, unparitioned data set. Bayesian hypothesis tests also reject those topologies using six assumed partitioning strategies, and the two partitioning strategies presumably associated with the most powerful tests also reject a third previously proposed topology, in which Callisaurus and Cophosaurus are sister taxa. For both maximum likelihood and Bayesian methods with assumed partitions, those partitions defined by codon position and tRNA stem and nonstems explained the data better than other strategies examined. Bayes factor estimates comparing results of assumed partitions versus mixture models suggest that mixture models perform better than assumed partitions when the latter were not based on functional characteristics of the data, such as codon position and tRNA stem and nonstems. However, assumed partitions performed better than mixture models when functional differences were incorporated. We reiterate the importance of accounting for heterogeneous evolutionary processes in the analysis of complex data sets and emphasize the importance of implementing mixed model likelihood methods.     

The effects of nucleotide substitution model assumptions on estimates of nonparametric bootstrap support

The use of parameter-rich substitution models in molecular phylogenetics has been criticized on the basis that these models can cause a reduction both in accuracy and in the ability to discriminate among competing topologies. We have explored the relationship between nucleotide substitution model complexity and nonparametric bootstrap support under maximum likelihood (ML) for six data sets for which the true relationships are known with a high degree of certainty. We also performed equally weighted maximum parsimony analyses in order to assess the effects of ignoring branch length information during tree selection. We observed that maximum parsimony gave the lowest mean estimate of bootstrap support for the correct set of nodes relative to the ML models for every data set except one. For several data sets, we established that the exact distribution used to model among-site rate variation was critical for a successful phylogenetic analysis. Site-specific rate models were shown to perform very poorly relative to gamma and invariable sites models for several of the data sets most likely because of the gross underestimation of branch lengths. The invariable sites model also performed poorly for several data sets where this model had a poor fit to the data, suggesting that addition of the gamma distribution can be critical. Estimates of bootstrap support for the correct nodes often increased under gamma and invariable sites models relative to equal rates models. Our observations are contrary to the prediction that such models cause reduced confidence in phylogenetic hypotheses. Our results raise several issues regarding the process of model selection, and we briefly discuss model selection uncertainty and the role of sensitivity analyses in molecular phylogenetics.     

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.     

Molecular phylogenetic analysis of nycteribiid and streblid bat flies (Diptera: Brachycera, Calyptratae): implications for host associations and phylogeographic origins

Bat flies are a small but diverse group of highly specialized ectoparasitic, obligatory bloodsucking Diptera. For the first time, the phylogenetic relationships of 26 species and five subfamilies were investigated using four genes (18S rDNA, 16S rDNA, COII, and cytB) under three optimality criteria (maximum parsimony (MP), maximum likelihood (ML), and Bayesian inference). Tree topology tests of previous hypotheses were conducted under likelihood (Shimodaira-Hasegawa test). Major findings include the non-monophyly of the Streblidae and the recovery of an Old World- and a New World-Clade of bat flies. These data ambiguously resolve basal relationships between Hippoboscidae, Glossinidae, and bat flies. Recovered phylogenies resulted in either monophyly (Bayesian approach) or paraphyly (MP/ML topologies) of the bat flies, thus obscuring the potential number of possible associations with bats throughout the history of this group. Dispersal-vicariance analysis suggested the Neotropical region as the possible ancestral distribution area of the New World Streblidae and the Oriental region for the Old World bat flies. The genes examined show conflicting support across the nodes of the tree, particularly in the basal positions. Additionally, there is poor character support among all genes for the nodes associated with early hippoboscoid diversification. This results in extremely short basal branches, adding support to the idea of a rapid radiation among the four major groups of Hippoboscoidea.     

Secondary structure alignment and direct optimization of 28S rDNA sequences provide limited phylogenetic resolution in bark and ambrosia beetles (Curculionidae: Scolytinae)

Phylogenetic relationships in Scolytinae were reconstructed from 107 DNA sequences that spanned the D2 and D3 expansion segments, and related core regions of the nuclear large ribosomal subunit (28S). Sequences were analysed by parsimony and Bayesian analyses of aligned sequences aided by a new secondary structure model for the D2–D3 domains. Direct optimization was performed on ambiguous alignment regions in combination with fixed states optimization of unambiguous regions, but performed poorly compared to the Bayesian and parsimony analyses. Generally, the phylogenetic signal mainly resolved relationships within tribes, while deeper divergences were either not resolved or received marginal support. In addition to confirming several previously established clades, we found that Micracini formed the sister group to Cactopinus , a group of mainly cactus feeding scolytine beetles. Furthermore, Ipini was monophyletic with Pseudips and Acanthotomicus subtending to the most basal node of that clade. The monophyly of Corthylini, which consists of the bark and cone feeding Pityophtorina and the ambrosia fungus-feeding Corthylina, was supported in some of the analyses. A close relationship was found between Phloeotribus and the two Phloeosinini genera Chramesus and Pseudochramesus , suggesting an evolutionary trajectory for the origin of a lamellate antennal club in Phloeotribus .     

Reliability of Bayesian posterior probabilities and bootstrap frequencies in phylogenetics

Many empirical studies have revealed considerable differences between nonparametric bootstrapping and Bayesian posterior probabilities in terms of the support values for branches, despite claimed predictions about their approximate equivalence. We investigated this problem by simulating data, which were then analyzed by maximum likelihood bootstrapping and Bayesian phylogenetic analysis using identical models and reoptimization of parameter values. We show that Bayesian posterior probabilities are significantly higher than corresponding nonparametric bootstrap frequencies for true clades, but also that erroneous conclusions will be made more often. These errors are strongly accentuated when the models used for analyses are underparameterized. When data are analyzed under the correct model, nonparametric bootstrapping is conservative. Bayesian posterior probabilities are also conservative in this respect, but less so.     

Phylogenetic relationships of the dwarf boas and a comparison of Bayesian and bootstrap measures of phylogenetic support

Four New World genera of dwarf boas (Exiliboa, Trachyboa, Tropidophis, and Ungaliophis) have been placed by many systematists in a single group (traditionally called Tropidophiidae). However, the monophyly of this group has been questioned in several studies. Moreover, the overall relationships among basal snake lineages, including the placement of the dwarf boas, are poorly understood. We obtained mtDNA sequence data for 12S, 16S, and intervening tRNA-val genes from 23 species of snakes representing most major snake lineages, including all four genera of New World dwarf boas. We then examined the phylogenetic position of these species by estimating the phylogeny of the basal snakes. Our phylogenetic analysis suggests that New World dwarf boas are not monophyletic. Instead, we find Exiliboa and Ungaliophis to be most closely related to sand boas (Erycinae), boas (Boinae), and advanced snakes (Caenophidea), whereas Tropidophis and Trachyboa form an independent clade that separated relatively early in snake radiation. Our estimate of snake phylogeny differs significantly in other ways from some previous estimates of snake phylogeny. For instance, pythons do not cluster with boas and sand boas, but instead show a strong relationship with Loxocemus and Xenopeltis. Additionally, uropeltids cluster strongly with Cylindrophis, and together are embedded in what has previously been considered the macrostomatan radiation. These relationships are supported by both bootstrapping (parametric and nonparametric approaches) and Bayesian analysis, although Bayesian support values are consistently higher than those obtained from nonparametric bootstrapping. Simulations show that Bayesian support values represent much better estimates of phylogenetic accuracy than do nonparametric bootstrap support values, at least under the conditions of our study.     

Data partitions and complex models in Bayesian analysis: the phylogeny of Gymnophthalmid lizards

Phylogenetic studies incorporating multiple loci, and multiple genomes, are becoming increasingly common. Coincident with this trend in genetic sampling, model-based likelihood techniques including Bayesian phylogenetic methods continue to gain popularity. Few studies, however, have examined model fit and sensitivity to such potentially heterogeneous data partitions within combined data analyses using empirical data. Here we investigate the relative model fit and sensitivity of Bayesian phylogenetic methods when alternative site-specific partitions of among-site rate variation (with and without autocorrelated rates) are considered. Our primary goal in choosing a best-fit model was to employ the simplest model that was a good fit to the data while optimizing topology and/or Bayesian posterior probabilities. Thus, we were not interested in complex models that did not practically affect our interpretation of the topology under study. We applied these alternative models to a four-gene data set including one protein-coding nuclear gene (c-mos), one protein-coding mitochondrial gene (ND4), and two mitochondrial rRNA genes (12S and 16S) for the diverse yet poorly known lizard family Gymnophthalmidae. Our results suggest that the best-fit model partitioned among-site rate variation separately among the c-mos, ND4, and 12S + 16S gene regions. We found this model yielded identical topologies to those from analyses based on the GTR+I+G model, but significantly changed posterior probability estimates of clade support. This partitioned model also produced more precise (less variable) estimates of posterior probabilities across generations of long Bayesian runs, compared to runs employing a GTR+I+G model estimated for the combined data. We use this three-way gamma partitioning in Bayesian analyses to reconstruct a robust phylogenetic hypothesis for the relationships of genera within the lizard family Gymnophthalmidae. We then reevaluate the higher-level taxonomic arrangement of the Gymnophthalmidae. Based on our findings, we discuss the utility of nontraditional parameters for modeling among-site rate variation and the implications and future directions for complex model building and testing.