Evaluating nuclear protein-coding genes for phylogenetic utility in beetles

Contrasting patterns of phylogeographic relationships in sympatric sister species of ironclad beetles (Zopheridae: Phloeodes spp.) in California's Transverse Ranges

BACKGROUND

Comparative phylogeography of sympatric sibling species provides an opportunity to isolate the effects of geography and demographics on the evolutionary history of two lineages over the same, known time scale. In the current study, we investigated the phylogeographic structure of two zopherid beetle species, Phloeodes diabolicus and P. plicatus, where their ranges overlap in California's Transverse Ranges.

RESULTS

Although P. diabolicus and P. plicatus share similar habitats with largely overlapping distributions, the results of this study revealed different evolutionary histories for each species since divergence from their most recent common ancestor. In general, P. plicatus had higher genetic diversity, and more among population isolation than P. diabolicus. The mismatch distributions indicated that one major difference between the two species was the timing of population expansion. This result was consistent with genetic patterns revealed by the Phist values and genetic diversity. Lastly, there were no parallel genetic breaks at similar geographic barriers between the species.

CONCLUSIONS

Our data revealed that differential demographics rather than geography were responsible for the genetic patterns of the two species.

Phylogeny and biogeography of dolichoderine ants: effects of data partitioning and relict taxa on historical inference

Ants (Hymenoptera: Formicidae) are conspicuous organisms in most terrestrial ecosystems, often attaining high levels of abundance and diversity. In this study, we investigate the evolutionary history of a major clade of ants, the subfamily Dolichoderinae, whose species frequently achieve ecological dominance in ant communities. This group has also produced some of the world's most successful invasive ants. We use an extensive molecular data set ( approximately 9 kb of sequence data from 10 nuclear genes, covering 48 dolichoderine species and 6 outgroup taxa) to infer the phylogenetic relationships, divergence dates, and biogeographic history of these ants. We evaluate the effects of data partitioning and outgroup composition on phylogenetic inference by estimating relationships under a series of increasingly partitioned data sets and by running analyses both with and without Aneuretus simoni, a rare and localized species that is the nearest living relative of Dolichoderinae. We also examine the effects of excluding 2 data partitions with significant base composition heterogeneity. Our results reveal 4 well-supported and mutually exclusive clades of dolichoderines, corresponding to 4 newly defined tribes: Bothriomyrmecini (B), Dolichoderini (D), Leptomyrmecini (L), and Tapinomini (T). All Bayesian and likelihood analyses yield the same unrooted (ingroup-only) topology, ((D,L),(B,T)), with the outgroups attaching either on the Dolichoderini branch or on the Tapinomini branch. Placement of the root is highly sensitive to choice of model partition and to inclusion/exclusion of Aneuretus. Bayes' factors strongly favor the more partitioned models, and in these Tapinomini is recovered as sister to the remaining dolichoderines, but only if Aneuretus is included. Exclusion of Aneuretus precludes recovery of this topology in all but the most highly partitioned Bayesian analyses and then only with nonsignificant support, underscoring the importance of relict, taxonomically isolated taxa for phylogenetic inference. Removal of 2 partitions with heterogeneous base composition also markedly increases support for placement of the root on the Tapinomini branch. Our divergence date estimates and biogeographic analyses indicate that crown-group dolichoderines arose about 65 million years ago (Ma), although this was preceded by a substantial period (30 million years) of stem group evolution. The 4 extant tribes are estimated to have crown-group origins in the late Paleocene or Eocene (40-60 Ma). Tapinomini and Bothriomyrmecini originated in the Paleotropics and subsequently dispersed to other biogeographic regions. Crown-group Leptomyrmecini arose and diversified in the Neotropics, but they also gave rise to one clade that colonized Australia about 30 Ma and subsequently experienced a massive radiation on that continent. This event occurred later than the diversification of dolichoderines in the northern hemisphere, so that by the time dolichoderines came to dominate the Australian fauna they had already declined in abundance in the Holarctic region.

Placement of Conopidae (Diptera) within Schizophora based on mtDNA and nrDNA gene regions

The first attempt to phylogenetically place Conopidae using molecular characters, as well as the largest molecular analysis of relationships within Schizophora (Diptera) to date, is presented. Twenty-eight taxa from 11 acalyptrate families and seven acalyptrate superfamilies are represented. Nearly 12,800 bp of sequence data from 10 genes representing both mitochondrial (cytochrome oxidase I (COI), cytochrome b (cytB), and 12S) and nuclear genes (28S, the carbamoyl phosphate synthetase region of CAD (CAD), elongation factor-1alpha (EF-1alpha), white, alanyl-tRNA synthetase (AATS), triose phosphate isomerase (TPI), and phosphogluconate dehydrogenase (PGD)) are analysed. Parsimony and Bayesian analyses strongly support the monophyly of both Conopidae and Schizophora. While in the parsimony analysis, Conopidae are placed as sister to the remaining Schizophora, the Bayesian analysis recovers a Conopidae+Lauxaniidae clade. The value of nuclear, mitochondrial, ribosomal, and protein-coding gene sequence data for answering phylogenetic questions at different levels of divergence is evaluated.

Phylogenetic analysis and temporal diversification of mosquitoes (Diptera: Culicidae) based on nuclear genes and morphology

Background

Phylogenetic analyses provide a framework for examining the evolution of morphological and molecular diversity, interpreting patterns in biogeography, and achieving a stable classification. The generic and suprageneric relationships within mosquitoes (Diptera: Culicidae) are poorly resolved, making these subjects difficult to address.

Results

We carried out maximum parsimony and maximum likelihood, including Bayesian, analyses on a data set consisting of six nuclear genes and 80 morphological characters to assess their ability to resolve relationships among 25 genera. We also estimated divergence times based on sequence data and fossil calibration points, using Bayesian relaxed clock methods. Strong support was recovered for the basal position and monophyly of the subfamily Anophelinae and the tribes Aedini and Sabethini of subfamily Culicinae. Divergence times for major culicid lineages date to the early Cretaceous.

Conclusions

Deeper relationships within the family remain poorly resolved, suggesting the need for additional taxonomic sampling. Our results support the notion of rapid radiations early in the diversification of mosquitoes.

Polyphyly of Xylosandrus Reitter inferred from nuclear and mitochondrial genes (Coleoptera: Curculionidae: Scolytinae)

The Xyleborina ambrosia beetle genus Xylosandrus contains 54 species, several of which are of economic importance. The monophyly of the genus was tested using a data set comprised of multiple gene loci: 28S rDNA; the mitochondrial gene cytochrome oxidase I (COI); and the nuclear genes arginine kinase (ArgK), rudimentary (CAD), and Elongation Factor 1alpha (EF-1alpha). The nuclear protein-coding genes CAD and ArgK were used for the first time in phylogenetics of Scolytinae. Analyses were performed using Parsimony and Bayesian optimality criteria. Our analyses included 43 specimens representing 15 Xylosandrus species and 20 species from Amasa, Anisandrus, Cnestus, Euwallacea and Xyleborus, and two species from the outgroup genus Coccotrypes. All analyses recovered a polyphyletic Xylosandrus. Several species of Xylosandrus were consistently placed in clades with the genera Anisandrus and Cnestus with high support values (100% bootstrap support). Among these, was the economically important invasive species X. mutilatus, which was consistently recovered as part of the "Cnestus" clade. In our analyses, both CAD and ArgK demonstrated phylogenetic utility across varying nodal depths. Despite the selection of genes with signals at complementary phylogenetic depths, the data set used herein did not resolve the phylogeny of Xylosandrus and related genera. Since the taxon sample available for molecular work represents only a fraction of Xylosandrus species, a complete revision that combines molecular and morphological data in a total evidence approach is recommended for the genus.

Testing the monophyly of Calocidae (Insecta: Trichoptera) based on multiple molecular data

Calocidae constitute a hypothesised monophyletic group of caddisflies (Trichoptera) being geographically restricted to New Zealand (one genus) and Australia (five genera). This analysis tests the monophyly of the family based on sequences from five different molecular genes. The complete data set includes 29 species and covers a complete genus representation of the Calocidae as well as representatives of other families in which one or more calocid genera have been classified. Sequences from two mitochondrial (cytochrome oxidase I and 16S) and three nuclear (elongation factor 1-alpha, RNA polymerase-II, and Cadherin) genes were used, resulting in a 3958bp data set and 37.1% parsimony informative characters. The Cadherin (CAD) and RNA polymerase-II (POL-II) genes are used for the first time for revealing Trichoptera phylogenies. The character matrix was analyzed by using maximum parsimony (MP) and Bayesian criteria, the latter by applying three different partition strategies for comparison. Two most parsimonious trees were found, differing in the position of one clade within the sister-group to a monophyletic Calocidae. The Bayesian tree based on the maximum number of partitions differs from trees based on a reduced partition analysis with respect to taxa outside the current circumscription of Calocidae. Both the MP and Bayesian analyses left Calocidae monophyletic, with a monophyletic clade of all Australian genera being sister-group to the New Zealand genus. The results from the agreement subtree analysis demonstrates that CAD performs well both separately and in combination with other genes and adds substantial resolution to the calocid phylogeny in a combined MP analysis.

Monophyly of terrestrial adephagan beetles as indicated by three nuclear genes (Coleoptera: Carabidae and Trachypachidae)

The beetle suborder Adephaga is traditionally divided into two sections on the basis of habitat, terrestrial Geadephaga and aquatic Hydradephaga. Monophyly of both groups is uncertain, and the relationship of the two groups has implications for inferring habitat transitions within Adephaga. Here we examine phylogenetic relationships of these groups using evidence provided by DNA sequences from all four suborders of beetles, including 60 species of Adephaga, four Archostemata, three Myxophaga, and ten Polyphaga. We studied 18S ribosomal DNA and 28S ribosomal DNA, aligned with consideration of secondary structure, as well as the nuclear protein-coding gene wingless. Independent and combined Bayesian, likelihood, and parsimony analyses of all three genes supported placement of Trachypachidae in a monophyletic Geadephaga, although for analyses of 28S rDNA and some parsimony analyses only if Coleoptera is constrained to be monophyletic. Most analyses showed limited support for the monophyly of Hydradephaga. Outside of Adephaga, there is support from the ribosomal genes for a sister group relationship between Adephaga and Polyphaga. Within the small number of sampled Polyphaga, analyses of 18S rDNA, wingless, and the combined matrix supports monophyly of Polyphaga exclusive of Scirtoidea. Unconstrained analyses of the evolution of habitat suggest that Adephaga was ancestrally aquatic with one transition to terrestrial. However, in analyses constrained to disallow changes from aquatic to terrestrial habitat, the phylogenies imply two origins of aquatic habit within Adephaga.

Multigene analysis of phylogenetic relationships and divergence times of primate sucking lice (Phthiraptera: Anoplura)

Cospeciation between hosts and parasites offers a unique opportunity to use information from parasites to infer events in host evolutionary history. Although lice (Insecta: Phthiraptera) are known to cospeciate with their hosts and have frequently served as important markers to infer host evolutionary history, most molecular studies are based on only one or two markers. Resulting phylogenies may, therefore, represent gene histories (rather than species histories), and analyses of multiple molecular markers are needed to increase confidence in the results of phylogenetic analyses. Herein, we phylogenetically examine nine molecular markers in primate sucking lice (Phthiraptera: Anoplura) and we use these markers to estimate divergence times among louse lineages. Individual and combined analyses of these nine markers are, for the most part, congruent, supporting relationships hypothesized in previous studies. Only one marker, the nuclear protein-coding gene Histone 3, has a significantly different tree topology compared to the other markers. The disparate evolutionary history of this marker, however, has no significant effect on topology or nodal support in the combined phylogenetic analyses. Therefore, phylogenetic results from the combined data set likely represent a solid hypothesis of species relationships. Additionally, we find that simultaneous use of multiple markers and calibration points provides the most reliable estimates of louse divergence times, in agreement with previous studies estimating divergences among species. Estimates of phylogenies and divergence times also allow us to verify the results of [Reed, D.L., Light, J.E., Allen, J.M., Kirchman, J.J., 2007. Pair of lice lost or parasites regained: the evolutionary history of anthropoid primate lice. BMC Biol. 5, 7.]; there was probable contact between gorilla and archaic hominids roughly 3 Ma resulting in a host switch of Pthirus lice from gorillas to archaic hominids. Thus, these results provide further evidence that data from cospeciating organisms can yield important information about the evolutionary history of their hosts.

Phylogeny, scaling, and the generation of extreme forces in trap-jaw ants

Trap-jaw ants of the genus Odontomachus produce remarkably fast predatory strikes. The closing mandibles of Odontomachus bauri, for example, can reach speeds of over 60 m s–1. They use these jaw strikes for both prey capture and locomotion – by striking hard surfaces, they can launch themselves into the air. We tested the hypothesis that morphological variation across the genus is correlated with differences in jaw speeds and accelerations. We video-recorded jaw-strikes at 70 000–100 000 frames s–1 to measure these parameters and to model force production. Differences in mean speeds ranged from 35.9±7.7 m s–1 for O. chelifer, to 48.8±8.9 m s–1 for O. clarus desertorum. Differences in species' accelerations and jaw sizes resulted in maximum strike forces in the largest ants (O. chelifer) that were four times those generated by the smallest ants (O. ruginodis). To evaluate phylogenetic effects and make statistically valid comparisons, we developed a phylogeny of all sampled Odontomachus species and seven outgroup species (19 species total) using four genetic loci. Jaw acceleration and jaw-scaling factors showed significant phylogenetic non-independence, whereas jaw speed and force did not. Independent contrast (IC) values were used to calculate scaling relationships for jaw length, jaw mass and body mass, which did not deviate significantly from isometry. IC regression of angular acceleration and body size show an inverse relationship, but combined with the isometric increase in jaw length and mass results in greater maximum strike forces for the largest Odontomachus species. Relatively small differences (3%) between IC and species-mean based models suggest that any deviation from isometry in species' force production may be the result of recent selective evolution, rather than deep phylogenetic signal.