Time flies, a new molecular time-scale for brachyceran fly evolution without a clock

The insect order Diptera, the true flies, contains one of the four largest Mesozoic insect radiations within its suborder Brachycera. Estimates of phylogenetic relationships and divergence dates among the major brachyceran lineages have been problematic or vague because of a lack of consistent evidence and the rarity of well-preserved fossils. Here, we combine new evidence from nucleotide sequence data, morphological reinterpretations, and fossils to improve estimates of brachyceran evolutionary relationships and ages. The 28S ribosomal DNA (rDNA) gene was sequenced for a broad diversity of taxa, and the data were combined with recently published morphological scorings for a parsimony-based phylogenetic analysis. The phylogenetic topology inferred from the combined 28S rDNA and morphology data set supports brachyceran monophyly and the monophyly of the four major brachyceran infraorders and suggests relationships largely consistent with previous classifications. Weak support was found for a basal brachyceran clade comprising the infraorders Stratiomyomorpha (soldier flies and relatives), Xylophagomorpha (xylophagid flies), and Tabanomorpha (horse flies, snipe flies, and relatives). This topology and similar alternative arrangements were used to obtain Bayesian estimates of divergence times, both with and without the assumption of a constant evolutionary rate. The estimated times were relatively robust to the choice of prior distributions. Divergence times based on the 28S rDNA and several fossil constraints indicate that the Brachycera originated in the late Triassic or earliest Mesozoic and that all major lower brachyceran fly lineages had near contemporaneous origins in the mid-Jurassic prior to the origin of flowering plants (angiosperms). This study provides increased resolution of brachyceran phylogeny, and our revised estimates of fly ages should improve the temporal context of evolutionary inferences and genomic comparisons between fly model organisms.     

Unraveling the evolutionary radiation of the thoracican barnacles using molecular and morphological evidence: a comparison of several divergence time estimation approaches

The Thoracica includes the ordinary barnacles found along the sea shore and is the most diverse and well-studied superorder of Cirripedia. However, although the literature abounds with scenarios explaining the evolution of these barnacles, very few studies have attempted to test these hypotheses in a phylogenetic context. The few attempts at phylogenetic analyses have suffered from a lack of phylogenetic signal and small numbers of taxa. We collected DNA sequences from the nuclear 18S, 28S, and histone H3 genes and the mitochondrial 12S and 16S genes (4,871 bp total) and data for 37 adult and 53 larval morphological characters from 43 taxa representing all the extant thoracican suborders (except the monospecific Brachylepadomorpha). Four Rhizocephala (highly modified parasitic barnacles) taxa and a Rhizocephala + Acrothoracica (burrowing barnacles) hypothetical ancestor were used as the outgroup for the molecular and morphological analyses, respectively. We analyzed these data separately and combined using maximum likelihood (ML) under "hill-climbing" and genetic algorithm heuristic searches, maximum parsimony procedures, and Bayesian inference coupled with Markov chain Monte Carlo techniques under mixed and homogeneous models of nucleotide substitution. The resulting phylogenetic trees answered key questions in barnacle evolution. The four-plated Iblomorpha were shown as the most primitive thoracican, and the plateless Heteralepadomorpha were placed as the sister group of the Lepadomorpha. These relationships suggest for the first time in an invertebrate that exoskeleton biomineralization may have evolved from phosphatic to calcitic. Sessilia (nonpedunculate) barnacles were depicted as monophyletic and appear to have evolved from a stalked (pedunculate) multiplated (5+) scalpelloidlike ancestor rather than a five-plated lepadomorphan ancestor. The Balanomorpha (symmetric sessile barnacles) appear to have the following relationship: (Chthamaloidea(Coronuloidea(Tetraclitoidea, Balanoidea))). Thoracican divergence times were estimated under ML-based local clock, Bayesian, and penalized likelihood approaches using an 18S data set and three calibration points: Heteralepadomorpha = 530 million years ago (MYA), Scalpellomorpha = 340 MYA, and Verrucomorpha = 120 MYA. Estimated dates varied considerably within and between approaches depending on the calibration point. Highly parameterized local clock models that assume independent rates (r > or = 15) for confamilial or congeneric species generated the most congruent estimates among calibrations and agreed more closely with the barnacle fossil record. Reasonable estimates were also obtained under the Bayesian procedure of Kishino et al. (2001, Mol. Biol. Evol. 18:352-361) but using multiple calibrations. Most of the dates estimated under the Bayesian procedure of Aris-Brosou and Yang (2002, Syst. Biol. 51:703-714) and the penalized likelihood method using single and/or multiple calibrations were inconsistent among calibrations and did not fit the fossil record.     

Phylogenetic analysis of the Brachyura (Crustacea, Decapoda) based on characters of the foregut with establishment of a new taxon

The Brachyura, within the decapod crustaceans, is one of the most species-rich taxa with up to 10 000 species. However, its phylogenetic history, evolution and fossil record remain subjects of controversy. In our study, we examined the phylogenetic relationships of the Brachyura based on morphological characters of the foregut. The cladistic analysis supports a monophyletic Brachyura including the Dromiidae and Raninidae. A clade comprising Dromiidae and Dynomenidae forms the most basal assemblage within the Brachyura, followed by the Homolidae and Latreilliidae. As a result, neither Podotremata nor Archaeobrachyura form a clade. In contrast, foregut data suggest that the classical taxon Oxystomata, comprising Calappidae, Parthenopidae, Dorippidae, Leucosiidae, Cymonomidae and Raninidae, is monophyletic. This makes the Heterotremata paraphyletic or polyphyletic. A newly established taxon, Neobrachyura, embraces some representatives of the Heterotremata and the monophyletic Thoracotremata.     

Phylogeny of Decapoda using two nuclear protein-coding genes: origin and evolution of the Reptantia

The phylogeny of Decapoda is contentious and many hypotheses have been proposed based on morphological cladistic analyses. Recent molecular studies, however, yielded contrasting results despite their use of similar data (nuclear and mitochondrial rDNA). Here we present the first application of two nuclear protein-coding genes, phosphoenolpyruvate carboxykinase and sodium-potassium ATPase alpha-subunit, to reconstruct the phylogeny of major infraorders within Decapoda. A total of 64 species representing all infraorders of Pleocyemata were analyzed with five species from Dendrobranchiata as outgroups. Maximum likelihood and Bayesian inference reveal that the Reptantia and all but one infraorder are monophyletic. Thalassinidea, however, is polyphyletic. The nodal support for most of the infraordinal and inter-familial relationships is high. Stenopodidea and Caridea form a clade sister to Reptantia, which comprises two major clades. The first clade, consisting of Astacidea, Achelata, Polychelida and three thalassinidean families (Axiidae, Calocarididae and Eiconaxiidae), corresponds essentially to the old taxon suborder Macrura Reptantia. Polychelida nests within Macrura Reptantia instead of being the most basal reptant as suggested in previous studies. The high level of morphological and genetic divergence of Polychelida from Achelata and Astacidea justifies its infraorder status. The second major reptant clade consists of Anomura, Brachyura and two thalassindean families (Thalassinidae and Upogebiidae). Anomura and Brachyura form Meiura, with moderate support. Notably thalassinidean families are sister to both major reptant clades, suggesting that the stem lineage reptants were thalassinidean-like. Moreover, some families (e.g. Nephropidae, Diogenidae, Paguridae) are paraphyletic, warranting further studies to evaluate their status. The present study ably demonstrates the utility of nuclear protein-coding genes in phylogenetic inference in decapods. The topologies obtained are robust and the two molecular markers are informative across a wide range of taxonomic levels. We propose that nuclear protein-coding genes should constitute core markers for future phylogenetic studies of decapods, especially for higher systematics.     

Molecular clocks do not support the Cambrian explosion

The fossil record has long supported the view that most animal phyla originated during a brief period approximately 520 MYA known as the Cambrian explosion. However, molecular data analyses over the past 3 decades have found deeper divergences among animals (approximately 800 to 1,200 MYA), with and without the assumption of a global molecular clock. Recently, two studies have instead reported time estimates apparently consistent with the fossil record. Here, we demonstrate that methodological problems in these studies cast doubt on the accuracy and interpretations of the results obtained. In the study by Peterson et al., young time estimates were obtained because fossil calibrations were used as maximum limits rather than as minimum limits, and not because invertebrate calibrations were used. In the study by Aris-Brosou and Yang, young time estimates were obtained because of problems with rate models and other methods specific to the study, and not because Bayesian methods were used. This also led to many anomalous findings in their study, including a primate-rodent divergence at 320 MYA. With these results aside, molecular clocks continue to support a long period of animal evolution before the Cambrian explosion of fossils.     

Are pollen fossils useful for calibrating relaxed molecular clock dating of phylogenies? A comparative study using Myrtaceae

The identification and application of reliable fossil calibrations represents a key component of many molecular studies of evolutionary timescales. In studies of plants, most paleontological calibrations are associated with macrofossils. However, the pollen record can also inform age calibrations if fossils matching extant pollen groups are found. Recent work has shown that pollen of the myrtle family, Myrtaceae, can be classified into a number of morphological groups that are synapomorphic with molecular groups. By assembling a data matrix of pollen morphological characters from extant and fossil Myrtaceae, we were able to measure the fit of 26 pollen fossils to a molecular phylogenetic tree using parsimony optimisation of characters. We identified eight Myrtaceidites fossils as appropriate for calibration based on the most parsimonious placements of these fossils on the tree. These fossils were used to inform age constraints in a Bayesian phylogenetic analysis of a sequence alignment comprising two sequences from the chloroplast genome (matK and ndhF) and one nuclear locus (ITS), sampled from 106 taxa representing 80 genera. Three additional analyses were calibrated by placing pollen fossils using geographic and morphological information (eight calibrations), macrofossils (five calibrations), and macrofossils and pollen fossils in combination (12 calibrations). The addition of new fossil pollen calibrations led to older crown ages than have previously been found for tribes such as Eucalypteae and Myrteae. Estimates of rate variation among lineages were affected by the choice of calibrations, suggesting that the use of multiple calibrations can improve estimates of rate heterogeneity among lineages. This study illustrates the potential of including pollen-based calibrations in molecular studies of divergence times.     

Higher level phylogeny and the first divergence time estimation of Heteroptera (Insecta: Hemiptera) based on multiple genes

Heteroptera, or true bugs, are the largest, morphologically diverse and economically important group of insects with incomplete metamorphosis. However, the phylogenetic relationships within Heteroptera are still in dispute and most of the previous studies were based on morphological characters or with single gene (partial or whole 18S rDNA). Besides, so far, divergence time estimates for Heteroptera totally rely on the fossil record, while no studies have been performed on molecular divergence rates. Here, for the first time, we used maximum parsimony (MP), maximum likelihood (ML) and Bayesian inference (BI) with multiple genes (18S rDNA, 28S rDNA, 16S rDNA and COI) to estimate phylogenetic relationships among the infraorders, and meanwhile, the Penalized Likelihood (r8s) and Bayesian (BEAST) molecular dating methods were employed to estimate divergence time of higher taxa of this suborder. Major results of the present study included: Nepomorpha was placed as the most basal clade in all six trees (MP trees, ML trees and Bayesian trees of nuclear gene data and four-gene combined data, respectively) with full support values. The sister-group relationship of Cimicomorpha and Pentatomomorpha was also strongly supported. Nepomorpha originated in early Triassic and the other six infraorders originated in a very short period of time in middle Triassic. Cimicomorpha and Pentatomomorpha underwent a radiation at family level in Cretaceous, paralleling the proliferation of the flowering plants. Our results indicated that the higher-group radiations within hemimetabolous Heteroptera were simultaneously with those of holometabolous Coleoptera and Diptera which took place in the Triassic. While the aquatic habitat was colonized by Nepomorpha already in the Triassic, the Gerromorpha independently adapted to the semi-aquatic habitat in the Early Jurassic.     

Phylogenetic systematics of the reptantian Decapoda (Crustacea, Malacostraca)

Although the biology of the reptantian Decapoda has been much studied, the last comprehensive review of reptantian systematics was published more than 80 years ago. We have used cladistic methods to reconstruct the phylogenetic system of the reptantian Decapoda. We can show that the Reptantia represent a monophyletic taxon. The classical groups, the 'Palinura', 'Astacura' and 'Anomura' are paraphyletic assemblages. The Polychelida is the sister-group of all other reptantians. The Astacida is not closely related to the Homarida, but is part of a large monophyletic taxon which also includes the Thalassinida, Anomala and Brachyura. The Anomala and Brachyura are sister-groups and the Thalassinida is the sister-group of both of them. Based on our reconstruction of the sister-group relationships within the Reptantia, we discuss alternative hypotheses of reptantian interrelationships, the systematic position of the Reptantia within the decapods, and draw some conclusions concerning the habits and appearance of the reptantian stem species.     

Phylogeny of the parasitic microgastroid subfamilies (Hymenoptera: Braconidae) based on sequence data from seven genes, with an improved time estimate of the origin of the lineage

The microgastroid complex of braconid wasps is a widely recognized and biologically coherent lineage of endoparasitoids of lepidopteran larvae (caterpillars). The complex has received significant phylogenetic attention in recent years due in part to the taxons' association with mutualistic polydnaviruses, with which they compromise host immune systems. A number of previous attempts using a variety of morphological and molecular approaches have not unequivocally resolved relationships amongst the main subfamilies. This work represents a more extensive attempt to resolve the microgastroid relationships, using seven genes (16S rRNA, cytochrome oxidase I (CO1), 28S rRNA, arginine kinase (ArgK), long wavelength rhodopsin (Ops), elongation factor 1 alpha (EF1a) and wingless (Wg)) and a greater taxonomic representation. Bayesian, likelihood and parsimony phylogenetic reconstructions of this improved data set has determined that the chelonines diverged first from the remainder of the microgastroids, however the relationships amongst the other subfamilies are still unclear, suggesting a greater nucleotide sample is required to resolve them. Examination of the contribution of individual gene trees to the phylogeny demonstrates why the relationships between subfamilies are still unclear, with not all groups monophyletic for all trees. Filtered supernetworks demonstrate that monophyly of all subfamilies is only recovered when splits found in only one or two genes are excluded, but this also results in little remaining structure left in the deep nodes to resolve inter-subfamily relationships. By increasing the breadth of the study we were also able to re-evaluate previous attempts at dating the lineage and, therefore the origin of the polydnavirus association. Previous attempts used a much reduced data set and fewer fossil calibrations. Thorough literature searches have revealed a substantial increase in the fossil calibrations and these, combined with more sophisticated molecular dating analysis, have substantially increased the age of the microgastroid lineage from previous estimates of approximately 73MYA to approximately 100MYA. Examination of the resultant linearized clock tree also allows an insight into the evolution of the more species rich subfamilies. The chelonines appear to have had a steady rate of evolution, whilst the microgastrines and cardiochilines appear to have undergone a more significant "burst" of evolution. It is hypothesized that the different parasitism strategies of subfamilies (Chelonines are egg parasitoids and the remainder are larval parasitioids) may have influenced the evolutionary rates of the groups.     

Mitogenomic analysis of decapod crustacean phylogeny corroborates traditional views on their relationships

Phylogenetic relationships within decapod crustaceans are highly controversial. Even recent analyses based on molecular datasets have shown largely contradictory results. Previous studies using mitochondrial genomes are promising but suffer from a poor and unbalanced taxon sampling. To fill these gaps we sequenced the (nearly) complete mitochondrial genomes of 13 decapod species: Stenopus hispidus, Polycheles typhlops, Panulirus versicolor, Scyllarides latus, Enoplometopus occidentalis, Homarus gammarus, Procambarus fallax f. virginalis, Upogebia major, Neaxius acanthus, Calocaris macandreae, Corallianassa coutierei, Cryptolithodes sitchensis, Neopetrolisthes maculatus, and add that of Dromia personata. Our new data allow for comprehensive analyses of decapod phylogeny using the mitochondrial genomes of 50 species covering all major taxa of the Decapoda. Five species of Stomatopoda and one species of Euphausiacea serve as outgroups. Most of our analyses using Maximum Likelihood (ML) and Bayesian inference (BI) of nucleotide and amino acid datasets revealed congruent topologies for higher level decapod relationships: (((((((Anomala, Brachyura), Thalassinida: Gebiidea), Thalassinida: Axiidea), (Astacidea, Polychelida), Achelata), Stenopodidea), Caridea), Dendrobranchiata). This result corroborates several traditional morphological views and adds new perspectives. In particular, the position of Polychelida is surprising. Nevertheless, some problems can be identified. In a minority of analyses the basal branching of Reptantia is not fully resolved, Thalassinida are monophyletic; Polychelida are the sister group to Achelata, and Stenopodidea are resolved as sister group to Caridea. Despite this and although some nodal supports are low in our phylogenetic trees, we think that the largely stable topology of the trees regardless of different types of analyses suggests that mitochondrial genomes show good potential to resolve the relationship within Decapoda.     

Molecular phylogenetics and biogeography of Lepus in Eastern Asia based on mitochondrial DNA sequences

In spite of several classification attempts among taxa of the genus Lepus, phylogenetic relationships still remain poorly understood. Here, we present molecular genetic evidence that may resolve some of the current incongruities in the phylogeny of the leporids. The complete mitochondrial cytb, 12S genes, and parts of ND4 and control region fragments were sequenced to examine phylogenetic relationships among Chinese hare taxa and other leporids throughout the World using maximum parsimony, maximum likelihood, and Bayesian phylogenetic reconstruction approaches. Using reconstructed phylogenies, we observed that the Chinese hare is not a single monophyletic group as originally thought. Instead, the data infers that the genus Lepus is monophyletic with three unique species groups: North American, Eurasian, and African. Ancestral area analysis indicated that ancestral Lepus arose in North America and then dispersed into Eurasia via the Bering Land Bridge eventually extending to Africa. Brooks Parsimony analysis showed that dispersal events followed by subsequent speciation have occurred in other geographic areas as well and resulted in the rapid radiation and speciation of Lepus. A Bayesian relaxed molecular clock approach based on the continuous autocorrelation of evolutionary rates along branches estimated the divergence time between the three major groups within Lepus. The genus appears to have arisen approximately 10.76 MYA (+/-0.86 MYA), with most speciation events occurring during the Pliocene epoch (5.65+/-1.15 MYA approximately 1.12 +/- 0.47 MYA).     

Evolutionary relationships among scyphozoan jellyfish families based on complete taxon sampling and phylogenetic analyses of 18S and 28S ribosomal DNA

A stable phylogenetic hypothesis for families within jellyfish class Scyphozoa has been elusive. Reasons for the lack of resolution of scyphozoan familial relationships include a dearth of morphological characters that reliably distinguish taxa and incomplete taxonomic sampling in molecular studies. Here, we address the latter issue by using maximum likelihood and Bayesian methods to reconstruct the phylogenetic relationships among all 19 currently valid scyphozoan families, using sequence data from two nuclear genes: 18S and 28S rDNA. Consistent with prior morphological hypotheses, we find strong evidence for monophyly of subclass Discomedusae, order Coronatae, rhizostome suborder Kolpophorae and superfamilies Actinomyariae, Kampylomyariae, Krikomyariae, and Scapulatae. Eleven of the 19 currently recognized scyphozoan families are robustly monophyletic, and we suggest recognition of two new families pending further analyses. In contrast to long-standing morphological hypotheses, the phylogeny shows coronate family Nausithoidae, semaeostome family Cyaneidae, and rhizostome suborder Daktyliophorae to be nonmonophyletic. Our analyses neither strongly support nor strongly refute monophyly of order Rhizostomeae, superfamily Inscapulatae, and families Ulmaridae, Catostylidae, Lychnorhizidae, and Rhizostomatidae. These taxa, as well as familial relationships within Coronatae and within rhizostome superfamily Inscapulatae, remain unclear and may be resolved by additional genomic and taxonomic sampling. In addition to clarifying some historically difficult taxonomic questions and highlighting nodes in particular need of further attention, the molecular phylogeny presented here will facilitate more robust study of phenotypic evolution in the Scyphozoa, including the evolution characters associated with mass occurrences of jellyfish.     

Fossil calibration of molecular divergence infers a moderate mutation rate and recent radiations for pinus

Silent mutation rate estimates for Pinus vary 50-fold, ranging from angiosperm-like to among the slowest reported for plants. These differences either reflect extraordinary genomic processes or inconsistent fossil calibration, and they have important consequences for population and biogeographical inferences. Here we estimate mutation rates from 4 Pinus species that represent the major lineages using 11 nuclear and 4 chloroplast loci. Calibration was tested at the divergence of Pinus subgenera with the oldest leaf fossil from subg. Strobus (Eocene; 45 MYA) or a recently published subg. Strobus wood fossil (Cretaceous; 85 MYA). These calibrations place the origin of Pinus 190-102 MYA and give absolute silent rate estimates of 0.70-1.31x10(-9) and 0.22-0.42x10(-9).site-1.year-1 for the nuclear and chloroplast genomes, respectively. These rates are approximately 4- to 20-fold slower than angiosperms, but unlike many previous estimates, they are more consistent with the high per-generation deleterious mutation rates observed in pines. Chronograms from nuclear and chloroplast genomes show that the divergence of subgenera accounts for about half of the time since Pinus diverged from Picea, with subsequent radiations occurring more recently. By extending the sampling to encompass the phylogenetic diversity of Pinus, we predict that most extant subsections diverged during the Miocene. Moreover, subsect. Australes, Ponderosae, and Contortae, containing over 50 extant species, radiated within a 5 Myr time span starting as recently as 18 MYA. An Eocene divergence of pine subgenera (using leaf fossils) does not conflict with fossil-based estimates of the Pinus-Picea split, but a Cretaceous divergence using wood fossils accommodates Oligocene fossils that may represent modern subsections. Because homoplasy and polarity of character states have not been tested for fossil pine assignments, the choice of fossil and calibration node represents a significant source of uncertainty. Based on several lines of evidence (including agreement with ages inferred using calibrations outside of Pinus), we conclude that the 85 MYA calibration at the divergence of pine subgenera provides a reasonable lower bound and that further refinements in age and mutation rate estimates will require a synthetic examination of pine fossil history.     

Phylogenetic analysis of Myobia musculi (Schranck, 1781) by using the 18S small ribosomal subunit sequence

We used high-fidelity PCR to amplify 2 overlapping regions of the ribosomal gene complex from the rodent fur mite Myobia musculi. The amplicons encompassed a large portion of the mite's ribosomal gene complex spanning 3128 nucleotides containing the entire 18S rRNA, internal transcribed spacer (ITS) 1,5.8S rRNA, ITS2, and a portion of the 5'-end of the 28S rRNA. M. musculi's 179-nucleotide 5.8S rRNA nucleotide sequence was not conserved, so this region was identified by conservation of rRNA secondary structure. Maximum likelihood and Bayesian inference phylogenetic analyses were performed by using multiple sequence alignment consisting of 1524 nucleotides of M. musculi 18S rRNA and homologous sequences from 42 prostigmatid mites and the tick Dermacentor andersoni. The phylograms produced by both methods were in agreement regarding terminal, secondary, and some tertiary phylogenetic relationships among mites. Bayesian inference discriminated most infraordinal relationships between Eleutherengona and Parasitengona mites in the suborder Anystina. Basal relationships between suborders Anystina and Eupodina historically determined by comparing differences in anatomic characteristics were less well-supported by our molecular analysis. Our results recapitulated similar 18S rRNA sequence analyses recently reported. Our study supports M. musculi as belonging to the suborder Anystina, infraorder Eleutherenona, and superfamily Cheyletoidea.     

Calibration choice, rate smoothing, and the pattern of tetrapod diversification according to the long nuclear gene RAG-1

A phylogeny of tetrapods is inferred from nearly complete sequences of the nuclear RAG-1 gene sampled across 88 taxa encompassing all major clades, analyzed via parsimony and Bayesian methods. The phylogeny provides support for Lissamphibia, Theria, Lepidosauria, a turtle-archosaur clade, as well as most traditionally accepted groupings. This tree allows simultaneous molecular clock dating for all tetrapod groups using a set of well-corroborated calibrations. Relaxed clock (PLRS) methods, using the amniote = 315 Mya (million years ago) calibration or a set of consistent calibrations, recovers reasonable divergence dates for most groups. However, the analysis systematically underestimates divergence dates within archosaurs. The bird-crocodile split, robustly documented in the fossil record as being around approximately 245 Mya, is estimated at only approximately 190 Mya, and dates for other divergences within archosaurs are similarly underestimated. Archosaurs, and particulary turtles have slow apparent rates possibly confounding rate modeling, and inclusion of calibrations within archosaurs (despite their high deviances) not only improves divergence estimates within archosaurs, but also across other groups. Notably, the monotreme-therian split ( approximately 210 Mya) matches the fossil record; the squamate radiation ( approximately 190 Mya) is younger than suggested by some recent molecular studies and inconsistent with identification of approximately 220 and approximately 165 Myo (million-year-old) fossils as acrodont iguanians and approximately 95 Myo fossils colubroid snakes; the bird-lizard (reptile) split is considerably older than fossil estimates (< or = 285 Mya); and Sphenodon is a remarkable phylogenetic relic, being the sole survivor of a lineage more than a quarter of a billion years old. Comparison with other molecular clock studies of tetrapod divergences suggests that the common practice of enforcing most calibrations as minima, with a single liberal maximal constraint, will systematically overestimate divergence dates. Similarly, saturation of mitochondrial DNA sequences, and the resultant greater compression of basal branches means that using only external deep calibrations will also lead to inflated age estimates within the focal ingroup.     

Phylogenetic divergences of the true bugs (Insecta: Hemiptera: Heteroptera), with emphasis on the aquatic lineages: the last piece of the aquatic insect jigsaw originated in the Late Permian/Early Triassic

Heteroptera are among the most diverse hemimetabolous insects. Seven infraorders have been recognized within this suborder of Hemiptera. Apart from the well‐established sister‐group relationship between Cimicomorpha and Pentatomomorpha (= Terheteroptera), the two terminal lineages, the relationships among the other five infraorders are still controversial, of which three (Gerromorpha, Nepomorpha and Leptopodomorpha) are intimately connected to aquatic environments. However, the various and often conflicting available phylogeny hypotheses do not offer a clear background for a connection between diversification and palaeoenvironments. In this study, a molecular data set representing 79 taxa and 10 149 homologous sites is used to infer the phylogenetic relationships within Heteroptera. Bayesian inference, maximum‐likelihood and maximum parsimony analyses were employed. The results of phylogenetic inferences largely confirm the widely accepted phylogenetic context. Estimation of the divergence time based on the phylogenetic results revealed that Gerromorpha, Nepomorpha and Leptopodomorpha originated successively during the period from the Late Permian to Early Triassic (269–246 Ma). This timescale is consistent with the origin and radiation time of various aquatic holometabolans. Our results indicate that the aquatic and semi‐aquatic true bugs evolved under environmental conditions of high air temperature and humidity in an evolutionary scenario similar to that of the aquatic holometabolans.     

Phylogenetics and temporal diversification of the earliest true flies (Insecta: Diptera) based on multiple nuclear genes

Abstract Relationships among families of the lower Diptera (formerly suborder ‘Nematocera’) have been exceptionally difficult to resolve. Multiple hypotheses based on morphology have been proposed to identify the earliest lineages of flies and place the phylogenetic origin of the higher flies (Brachycera), but convincing support is limited. Here we resolve relationships among the major groups of lower Diptera using sequence data from four nuclear markers, including both ribosomal (28S rDNA) and protein‐coding (CAD, TPI and PGD) genes. Our results support both novel and traditional arrangements. Most unexpectedly, the small, highly‐specialized family Deuterophlebiidae appears to be sister to all remaining Diptera. Other results include the resolution of the traditional infra‐orders Culicomorpha (including a novel superfamily Simulioidea = Thaumaleidae + Simuliidae), Tipulomorpha (Tipulidae sensu lato + Trichoceridae) and Bibionomorpha sensu lato. We find support for a limited Psychodomorpha (Blephariceridae, Tanyderidae and Psychodidae) and Ptychopteromorpha (Ptychopteridae), whereas the placement of several enigmatic families (Nymphomyiidae, Axymyiidae and Perissommatidae) remains ambiguous. According to genetic data, the infra‐order Bibionomorpha is sister to the Brachycera. Much of the phylogenetic signal for major lineages was found in the 28S rDNA gene, whereas protein‐coding genes performed variably at different levels. In addition to elucidating relationships, we also estimate the age of major lower dipteran clades, based on molecular divergence time estimates using relaxed‐clock Bayesian methods and fossil calibration points.     

The phylogenetic relationships among infraorders and superfamilies of Diptera based on morphological evidence

Members of the megadiverse insect order Diptera (flies) have successfully colonized all continents and nearly all habitats. There are more than 154 000 described fly species, representing 10–12% of animal species. Elucidating the phylogenetic relationships of such a large component of global biodiversity is challenging, but significant advances have been made in the last few decades. Since Hennig first discussed the monophyly of major groupings, Diptera has attracted much study, but most researchers have used non‐numerical qualitative methods to assess morphological data. More recently, quantitative phylogenetic methods have been used on both morphological and molecular data. All previous quantitative morphological studies addressed narrower phylogenetic problems, often below the suborder or infraorder level. Here we present the first numerical analysis of phylogenetic relationships of the entire order using a comprehensive morphological character matrix. We scored 371 external and internal morphological characters from larvae, pupae and adults for 42 species, representing all infraorders selected from 42 families. Almost all characters were obtained from previous studies but required revision for this ordinal‐level study, with homology assessed beyond their original formulation and across all infraorders. We found significant support for many major clades (including the Diptera, Culicomorpha, Bibionomorpha, Brachycera, Eremoneura, Cyclorrhapha, Schizophora, Calyptratae and Oestroidea) and we summarize the character evidence for these groups. We found low levels of support for relationships between the infraorders of lower Diptera, lower Brachycera and major lineages of lower Cyclorrhapha, and this is consistent with findings from molecular studies. These poorly supported areas of the tree may be due to periods of rapid radiation that left few synapomorphies in surviving lineages.     

Gondwanan radiation of the Southern Hemisphere crayfishes (Decapoda: Parastacidae): evidence from fossils and molecules

Aim The sequential break‐up of Gondwana is thought to be a dominant process in the establishment of shared biota across landmasses of the Southern Hemisphere. Yet similar distributions are shared by taxa whose radiations clearly post‐date the Gondwanan break‐up. Thus, determining the contribution of vicariance versus dispersal to seemingly Gondwanan biota is complex. The southern freshwater crayfishes (family Parastacidae) are distributed on Australia and New Guinea, South America, Madagascar and New Zealand and are unlikely to have dispersed via oceans, owing to strict freshwater limitations. We test the hypotheses that the break‐up of Gondwana has led to (1) a predominately east–west (((Australia, New Zealand: 80 Ma) Madagascar: 160–121 Ma) South America: 165–140 Ma), or (2) a southern (((Australia, South America: 52–35 Ma) New Zealand: 80 Ma) Madagascar: 160–121 Ma) pattern for parastacid crayfish. Further, we examine the evidence for a complete drowning of New Zealand and subsequent colonization by freshwater crayfish. Location Southern Hemisphere. Methods The evolutionary relationships among the 15 genera of Parastacidae were reconstructed using mitochondrial [16S, cytochrome c oxidase subunit I (COI)] and nuclear (18S, 28S) sequence data and maximum likelihood and Bayesian methods of phylogenetic reconstruction. A Bayesian (multidivtime ) molecular dating method using six fossil calibrations and phylogenetic inference was used to estimate divergence time among crayfish clades on Gondwanan landmasses. Results The South American crayfish are monophyletic and a sister group to all other southern crayfish. Australian crayfish are not monophyletic, with two Tasmanian genera, Spinastacoides and Ombrastacoides, forming a clade with New Zealand and Malagasy crayfish (both monophyletic). Divergence of crayfish among southern landmasses is estimated to have occurred around the Late Jurassic to Early Cretaceous (109–178 Ma). Main conclusions The estimated phylogenetic relationships and time of divergence among the Southern Hemisphere crayfishes were consistent with an east–west pattern of Gondwanan divergence. The divergence between Australia and New Zealand (109–160 Ma) pre‐dated the rifting at around 80 Ma, suggesting that these lineages were established prior to the break‐up. Owing to the age of the New Zealand crayfish, we reject the hypothesis that there was a complete drowning of New Zealand crayfish habitat.     

The tempo and mode of barnacle evolution

Previous phylogenetic attempts at resolving barnacle evolutionary relationships are few and have relied on limited taxon sampling. Here we combine DNA sequences from three nuclear genes (18S, 28S and H3) and 44 morphological characters collected from 76 thoracican (ingroup) and 15 rhizocephalan (outgroup) species representing almost all the Thoracica families to assess the tempo and mode of barnacle evolution. Using phylogenetic methods of maximum parsimony, maximum likelihood, and Bayesian inference and 14 fossil calibrations, we found that: (1) Iblomorpha form a monophyletic group; (2) pedunculated barnacles without shell plates (Heteralepadomorpha) are not ancestral, but have evolved, at least twice, from plated forms; (3) the ontogenetic pattern with 5-->6-->8-->12+ plates does not reflect Thoracica shell evolution; (4) the traditional asymmetric barnacles (Verrucidae) and the Balanomorpha are each monophyletic and together they form a monophyletic group; (5) asymmetry and loss of a peduncle have evolved twice in the Thoracica, resulting in neither the Verrucomorpha nor the Sessilia forming monophyletic groups in their present definitions; (6) the Scalpellomorpha are not monophyletic; (7) the Thoracica suborders evolved since the Early Carboniferous (340mya) with the final radiation of the Sessilia in the Upper Jurassic (147mya). These results, therefore, reject many of the underlying hypotheses about character evolution in the Cirripedia Thoracica, stimulate a variety of new thoughts on thoracican radiation, and suggest the need for a major rearrangement in thoracican classification based on estimated phylogenetic relationships.