The insect molecular clock

Genomic analysis reveals strong population structure in the Giant Sydney Crayfish (Euastacus spinifer (Heller, 1865))

Australia is home to over 140 species of freshwater crayfish (Decapoda: Parastacidae), representing a centre of diversity for this group in the Southern Hemisphere. Species delimitation in freshwater crayfish is difficult because many species show significant variation in colouration and morphology. This is particularly evident in the genus Euastacus, which exhibits large variations in colour and spination throughout its putative range. To understand this variation, we investigated the genetic diversity, population structure, phylogeny, and evolutionary timescale of the Giant Sydney Crayfish (Euastacus spinifer (Heller, 1865)). Our data set is sampled from over 70 individuals from across the ∼600 km range of the species, and includes a combination of two mitochondrial markers and more than 7000 single-nucleotide polymorphisms (SNPs) from the nuclear genome. Data were also obtained for representatives of the close relative, Euastacus vesper McCormack and Ahyong, 2017. Genomic SNP analyses revealed strong population structure, with multiple distinct populations showing little evidence of gene flow or migration. Phylogenetic analyses of mitochondrial data revealed similar structure between populations. Taken together, our analyses suggest that E. spinifer, as currently understood, represents a species complex, of which E. vesper is a member. Molecular clock estimates place the divergences within this group during the Pleistocene. The isolated and highly fragmented populations identified in our analyses probably represent relict populations of a previously widespread ancestral species. Periodic flooding events during the Pleistocene are likely to have facilitated the movement of these otherwise restricted freshwater crayfish within and between drainage basins, including the Murray-Darling and South East Coast Drainages. We present evidence supporting the recognition of populations in the southern parts of the range of E. spinifer as one or two separate species, which would raise the number of species within the E. spinifer complex to at least three. Our results add to the growing body of evidence that many freshwater crayfish exhibit highly fragmented, range-restricted distributions. In combination with the life-history traits of these species, the restricted distributions exacerbate the threats already placed on freshwater crayfish, which are among the five most endangered animal groups globally.

The Impact of Fast Radiation on the Phylogeny of Bactrocera Fruit Flies as Revealed by Multiple Evolutionary Models and Mutation Rate-Calibrated Clock

Several true fruit flies (Tephritidae) cause major damage to agriculture worldwide. Among them, species of the genus Bactrocera are extensively studied to understand the traits associated with their invasiveness and ecology. Comparative approaches based on a reliable phylogenetic framework are particularly effective, but several nodes of the Bactrocera phylogeny are still controversial, especially concerning the reciprocal affinities of the two major pests B. dorsalis and B. tryoni. Here, we analyzed a newly assembled genomic-scaled dataset using different models of evolution to infer a phylogenomic backbone of ten representative Bactrocera species and two outgroups. We further provide the first genome-scaled inference of their divergence by calibrating the clock using fossil records and the spontaneous mutation rate. The results reveal a closer relationship of B. dorsalis with B. latifrons than to B. tryoni, contrary to what was previously supported by mitochondrial-based phylogenies. By employing coalescent-aware and heterogeneous evolutionary models, we show that this incongruence likely derives from a hitherto undetected systematic error, exacerbated by incomplete lineage sorting and possibly hybridization. This agrees with our clock analysis, which supports a rapid and recent radiation of the clade to which B. dorsalis, B. latifrons and B. tryoni belong. These results provide a new picture of Bactrocera phylogeny that can serve as the basis for future comparative analyses.

Complex co-evolutionary relationships between cicadas and their symbionts

Previous evidence suggests that cicadas lacking Hodgkinia may harbour the yeast-like fungal symbionts (YLS). Here, we reinforce an earlier conclusion that the pathogenic ancestor of YLS independently infected different cicada lineages instead of the common ancestor of Cicadidae. Five independent replacement events in the loss of Hodgkinia/acquisition of YLS and seven other replacement events of YLS (from an Ophiocordyceps fungus to another Ophiocordyceps fungus) are hypothesised to have occurred within the sampled cicada taxa. The divergence time of YLS lineages was later than that of corresponding cicada lineages. The rapid shift of diversification rates of YLS and related cicada-parasitizing Ophiocordyceps began at approximately 32.94 Ma, and the diversification rate reached the highest value at approximately 24.82 Ma, which corresponds to the cooling climate changes at the Eocene-Oligocene boundary and the Oligocene-Miocene transition respectively. Combined with related acquisition/replacement events of YLS occurred during the cooling-climate periods, we hypothesise that the cooling-climate changes impacted the interactions between cicadas and related Ophiocordyceps, which coupled with the unusual life cycle and the differentiation of cicadas may finally led to the diversification of YLS in Cicadidae. Our results contribute to a better understanding of the evolutionary transition of YLS from entomopathogenic fungi in insects.

Evolutionary Rates are Correlated Between Buchnera Endosymbionts and the Mitochondrial Genomes of Their Aphid Hosts

The evolution of bacterial endosymbiont genomes is strongly influenced by host-driven selection. Factors affecting host genome evolution will potentially affect endosymbiont genomes in similar ways. One potential outcome is correlations in molecular rates between the genomes of the symbiotic partners. Recently, we presented the first evidence of such correlations between the mitochondrial genomes of cockroaches and the genomes of their endosymbiont (Blattabacterium cuenoti). Here we investigate whether similar patterns are found in additional host-symbiont partners. We use partial genome data from multiple strains of the bacterial endosymbionts Buchnera aphidicola and Sulcia muelleri, and the mitochondrial genomes of their sap-feeding insect hosts. Both endosymbionts show phylogenetic congruence with the mitochondria of their hosts, a result that is expected due to their identical mode of inheritance. We compared root-to-tip distances and branch lengths of phylogenetically independent species pairs. Both analyses showed a highly significant correlation of molecular rates between the genomes of Buchnera and the mitochondrial genomes of their hosts. A similar correlation was detected between Sulcia and their hosts, but was not statistically significant. Our results indicate that evolutionary rate correlations between hosts and long-term symbionts may be a widespread phenomenon.

Chronological Incongruences between Mitochondrial and Nuclear Phylogenies of Aedes Mosquitoes

One-third of all mosquitoes belong to the Aedini, a tribe comprising common vectors of viral zoonoses such as Aedes aegypti and Aedes albopictus. To improve our understanding of their evolution, we present an updated multigene estimate of Aedini phylogeny and divergence, focusing on the disentanglement between nuclear and mitochondrial phylogenetic signals. We first show that there are some phylogenetic discrepancies between nuclear and mitochondrial markers which may be caused by wrong taxa assignment in samples collections or by some stochastic effect due to small gene samples. We indeed show that the concatenated dataset is model and framework dependent, indicating a general paucity of signal. Our Bayesian calibrated divergence estimates point toward a mosquito radiation in the mid-Jurassic and an Aedes radiation from the mid-Cretaceous on. We observe, however a strong chronological incongruence between mitochondrial and nuclear data, the latter providing divergence times within the Aedini significantly younger than the former. We show that this incongruence is consistent over different datasets and taxon sampling and that may be explained by either peculiar evolutionary event such as different levels of saturation in certain lineages or a past history of hybridization throughout the genus. Overall, our updated picture of Aedini phylogeny, reveal a strong nuclear-mitochondrial incongruence which may be of help in setting the research agenda for future phylogenomic studies of Aedini mosquitoes.

Molecular phylogeny and trait evolution in an ancient terrestrial arthropod lineage: Systematic revision and implications for ecological divergence (Collembola, Tomocerinae)

Phylogenetic assessments of functional traits are important for mechanistically understanding the interactions between organisms and environments, but such practices are strongly limited by the availability of phylogenetic frameworks. The tomocerin springtails are an ancient, widespread and ecologically important group of terrestrial arthropods, whereas their phylogeny and trait evolution remained unaddressed. In the present study, we conducted the first comprehensive phylogenetic reconstruction of Tomocerinae, based on a multi-loci molecular dataset covering all major lineages within the subfamily, using Bayesian inference (BI), maximum-likelihood (ML) and maximum-parsimony (MP) approaches. Divergence time was estimated and ancestral character state reconstruction (ACSR) was performed to trace the evolutionary history of five ecomorphological traits correlated with sensory and locomotory functions. Our results support the monophyly of Tomocerinae, and indicate that current classification of Tomocerinae only partially reflects evolutionary relationships, notably the commonest and speciose genus Tomocerus is polyphyletic. The subfamily probably originated in Early Cretaceous and diversified in two Cretaceous and one Eocene radiation events. As indicated by the evolutionary patterns of functional traits, multiple ecological divergences took place during the diversification of Tomocerinae. The study suggests a potential underestimation of ecological divergence and functional diversity in terrestrial arthropods, calls for an update of present trait databases, and demonstrates the value of macroevolutionary knowledge for improving the trait-based ecology. In addition, Tomocerus, Tomocerina and Tritomurus are redefined, a new genus Yoshiicerusgen. n. and new subgenera Coloratomurussubgen. n., Ciliatomurussubgen. n., Striatomurussubgen. n. and Ocreatomurussubgen. n. are described in the appendix.

Extreme genetic diversity among springtails (Collembola) in subterranean calcretes of arid Australia

The subterranean islands hypothesis for calcretes of the Yilgarn region in Western Australia applies to many stygobitic (subterranean-aquatic) species that are "trapped" evolutionarily within isolated aquifers due to their aquatic lifestyles. In contrast, little is known about the distribution of terrestrial-subterranean invertebrates associated with the calcretes. We used subterranean Collembola from the Yilgarn calcretes to test the hypothesis that troglobitic species, those inhabiting the subterranean unsaturated (non-aquatic) zone of calcretes, are also restricted in their distribution and represent reciprocally monophyletic and endemic lineages. We used the barcoding fragment of the mtDNA cytochrome c oxidase subunit 1 (COI) gene from 183 individuals to reconstruct the phylogenetic history of the genus Pseudosinella Schäffer (Collembola, Lepidocyrtidae) from 10 calcretes in the Yilgarn. These calcretes represent less than 5% of the total possible calcretes in this region, yet we show that their diversity for subterranean Collembola comprises a minimum of 25 new species. Regionally, multiple levels of diversity exist in Pseudosinella, indicative of a complex evolutionary history for this genus in the Yilgarn. These species have probably been impacted by climatic oscillations, facilitating their dispersal across the landscape. The results represent a small proportion of the undiscovered diversity in Australia's arid zone.

Bayesian Inference of Evolutionary Histories under Time-Dependent Substitution Rates

Many factors complicate the estimation of time scales for phylogenetic histories, requiring increasingly complex evolutionary models and inference procedures. The widespread application of molecular clock dating has led to the insight that evolutionary rate estimates may vary with the time frame of measurement. This is particularly well established for rapidly evolving viruses that can accumulate sequence divergence over years or even months. However, this rapid evolution stands at odds with a relatively high degree of conservation of viruses or endogenous virus elements over much longer time scales. Building on recent insights into time-dependent evolutionary rates, we develop a formal and flexible Bayesian statistical inference approach that accommodates rate variation through time. We evaluate the novel molecular clock model on a foamy virus cospeciation history and a lentivirus evolutionary history and compare the performance to other molecular clock models. For both virus examples, we estimate a similarly strong time-dependent effect that implies rates varying over four orders of magnitude. The application of an analogous codon substitution model does not implicate long-term purifying selection as the cause of this effect. However, selection does appear to affect divergence time estimates for the less deep evolutionary history of the Ebolavirus genus. Finally, we explore the application of our approach on woolly mammoth ancient DNA data, which shows a much weaker, but still important, time-dependent rate effect that has a noticeable impact on node age estimates. Future developments aimed at incorporating more complex evolutionary processes will further add to the broad applicability of our approach.

Evolutionary rates are correlated between cockroach symbionts and mitochondrial genomes

Bacterial endosymbionts evolve under strong host-driven selection. Factors influencing host evolution might affect symbionts in similar ways, potentially leading to correlations between the molecular evolutionary rates of hosts and symbionts. Although there is evidence of rate correlations between mitochondrial and nuclear genes, similar investigations of hosts and symbionts are lacking. Here, we demonstrate a correlation in molecular rates between the genomes of an endosymbiont (Blattabacterium cuenoti) and the mitochondrial genomes of their hosts (cockroaches). We used partial genome data for multiple strains of B. cuenoti to compare phylogenetic relationships and evolutionary rates for 55 cockroach/symbiont pairs. The phylogenies inferred for B. cuenoti and the mitochondrial genomes of their hosts were largely congruent, as expected from their identical maternal and cytoplasmic mode of inheritance. We found a correlation between evolutionary rates of the two genomes, based on comparisons of root-to-tip distances and on comparisons of the branch lengths of phylogenetically independent species pairs. Our results underscore the profound effects that long-term symbiosis can have on the biology of each symbiotic partner.

Host and geography together drive early adaptive radiation of Hawaiian planthoppers

The interactions between insects and their plant host have been implicated in driving diversification of both players. Early arguments highlighted the role of ecological opportunity, with the idea that insects "escape and radiate" on new hosts, with subsequent hypotheses focusing on the interplay between host shifting and host tracking, coupled with isolation and fusion, in generating diversity. Because it is rarely possible to capture the initial stages of diversification, it is particularly difficult to ascertain the relative roles of geographic isolation versus host shifts in initiating the process. The current study examines genetic diversity between populations and hosts within a single species of endemic Hawaiian planthopper, Nesosydne umbratica (Hemiptera, Delphacidae). Given that the species was known as a host generalist occupying unrelated hosts, Clermontia (Campanulaceae) and Pipturus (Urticaceae), we set out to determine the relative importance of geography and host in structuring populations in the early stages of differentiation on the youngest islands of the Hawaiian chain. Results from extensive exon capture data showed that N. umbratica is highly structured, both by geography, with discrete populations on each volcano, and by host plant, with parallel radiations on Clermontia and Pipturus leading to extensive co-occurrence. The marked genetic structure suggests that populations can readily become established on novel hosts provided opportunity; subsequent adaptation allows monopolization of the new host. The results support the role of geographic isolation in structuring populations and with host shifts occurring as discrete events that facilitate subsequent parallel geographic range expansion.

Variation across Species and Levels: Implications for Model Species Research

The selection of model species tends to involve two typically unstated assumptions, namely: (1) that the similarity between species decreases steadily with phylogenetic distance, and (2) that similarities are greater at lower levels of biological organization. The first assumption holds on average, but species similarities tend to decrease with the square root of divergence time, rather than linearly, and lineages with short generation times (which includes most model species) tend to diverge faster than average, making the decrease in similarity non-monotonic. The second assumption is more difficult to test. Comparative molecular research has traditionally emphasized species similarities over differences, whereas comparative research at higher levels of organization frequently highlights the species differences. However, advances in comparative genomics have brought to light a great variety of species differences, not just in gene regulation but also in protein coding genes. Particularly relevant are cases in which homologous high-level characters are based on non-homologous genes. This phenomenon of non-orthologous gene displacement, or "deep non-homology," indicates that species differences at the molecular level can be surprisingly large. Given these observations, it is not surprising that some findings obtained in model species do not generalize across species as well as researchers had hoped, even if the research is molecular.

Origin and macroevolution of micro-moths on sunken Hawaiian Islands

The origins and evolution of Hawaiian biodiversity are a matter of controversy, and the mechanisms of lineage diversification for many organisms on this remote archipelago remain unclear. Here we focus on the poorly known endemic leaf-mining moth genus Philodoria (Lepidoptera, Gracillariidae), whose species feed on a diversity of Hawaiian plant lineages, many of which are critically endangered. We use anchored hybrid enrichment to assemble the first phylogenomic dataset (507 loci) for any Hawaiian animal taxon. To uncover the timing and pattern of diversification of these moths, we apply two frequently used dating calibration strategies, biogeographic calibrations and secondary calibrations. Island calibrations on their own resulted in much younger and unrealistic dates compared to strategies that relied on secondary calibrations. Philodoria probably originated on the now partially sunken islands of Laysan or Lisianski, approximately 21 Ma, and were associated with host plants in the families Ebenaceae, Malvaceae or Primulaceae. Major feeding groups associated with specific host-plant families originated soon after the plants colonized the islands. Allopatric isolation and host shifts, in concert and independently, probably play major roles in the diversification of Philodoria Our dating results indicate that Philodoria is among the oldest known Hawaiian arthropod lineages, and that island calibrations alone can lead to unrealistically young dates.

Reconstructing the phylogeny of Blattodea: robust support for interfamilial relationships and major clades

Cockroaches are among the most recognizable of all insects. In addition to their role as pests, they play a key ecological role as decomposers. Despite numerous studies of cockroach phylogeny in recent decades, relationships among most major lineages are yet to be resolved. Here we examine phylogenetic relationships among cockroaches based on five genes (mitochondrial 12S rRNA, 16S rRNA, COII; nuclear 28S rRNA and histone H3), and infer divergence times on the basis of 8 fossils. We included in our analyses sequences from 52 new species collected in China, representing 7 families. These were combined with data from a recent study that examined these same genes from 49 species, resulting in a significant increase in taxa analysed. Three major lineages, Corydioidea, Blaberoidea, and Blattoidea were recovered, the latter comprising Blattidae, Tryonicidae, Lamproblattidae, Anaplectidae, Cryptocercidae and Isoptera. The estimated age of the split between Mantodea and Blattodea ranged from 204.3 Ma to 289.1 Ma. Corydioidea was estimated to have diverged 209.7 Ma (180.5-244.3 Ma 95% confidence interval [CI]) from the remaining Blattodea. The clade Blattoidea diverged from their sister group, Blaberoidea, around 198.3 Ma (173.1-229.1 Ma). The addition of the extra taxa in this study has resulted in significantly higher levels of support for a number of previously recognized groupings.

Time-dependent estimates of molecular evolutionary rates: evidence and causes

We are writing in response to a recent critique by Emerson & Hickerson (2015), who challenge the evidence of a time-dependent bias in molecular rate estimates. This bias takes the form of a negative relationship between inferred evolutionary rates and the ages of the calibrations on which these estimates are based. Here, we present a summary of the evidence obtained from a broad range of taxa that supports a time-dependent bias in rate estimates, with a consideration of the potential causes of these observed trends. We also describe recent progress in improving the reliability of evolutionary rate estimation and respond to the concerns raised by Emerson & Hickerson (2015) about the validity of rates estimated from time-structured sequence data. In doing so, we hope to dispel some misconceptions and to highlight several research directions that will improve our understanding of time-dependent biases in rate estimates.

Inflation of Molecular Clock Rates and Dates: Molecular Phylogenetics, Biogeography, and Diversification of a Global Cicada Radiation from Australasia (Hemiptera: Cicadidae: Cicadettini)

Dated phylogenetic trees are important for studying mechanisms of diversification, and molecular clocks are important tools for studies of organisms lacking good fossil records. However, studies have begun to identify problems in molecular clock dates caused by uncertainty of the modeled molecular substitution process. Here we explore Bayesian relaxed-clock molecular dating while studying the biogeography of ca. 200 species from the global cicada tribe Cicadettini. Because the available fossils are few and uninformative, we calibrate our trees in part with a cytochrome oxidase I (COI) clock prior encompassing a range of literature estimates for arthropods. We show that tribe-level analyses calibrated solely with the COI clock recover extremely old dates that conflict with published estimates for two well-studied New Zealand subclades within Cicadettini. Additional subclade analyses suggest that COI relaxed-clock rates and maximum-likelihood branch lengths become inflated relative to EF-1[Formula: see text] intron and exon rates and branch lengths as clade age increases. We present corrected estimates derived from: (i) an extrapolated EF-1[Formula: see text] exon clock derived from COI-calibrated analysis within the largest New Zealand subclade; (ii) post hoc scaling of the tribe-level chronogram using results from subclade analyses; and (iii) exploitation of a geological calibration point associated with New Caledonia. We caution that considerable uncertainty is generated due to dependence of substitution estimates on both the taxon sample and the choice of model, including gamma category number and the choice of empirical versus estimated base frequencies. Our results suggest that diversification of the tribe Cicadettini commenced in the early- to mid-Cenozoic and continued with the development of open, arid habitats in Australia and worldwide. We find that Cicadettini is a rare example of a global terrestrial animal group with an Australasian origin, with all non-Australasian genera belonging to two distal clades. Within Australia, we show that Cicadettini is more widely distributed than any other cicada tribe, diverse in temperate, arid and monsoonal habitats, and nearly absent from rainforests. We comment on the taxonomic implications of our findings for thirteen cicada genera.

Molecular Phylogeny of Grassland Caterpillars (Lepidoptera: Lymantriinae: Gynaephora) Endemic to the Qinghai-Tibetan Plateau

Gynaephora (Lepidoptera Erebidae: Lymantriinae) is a small genus, consisting of 15 nominated species, of which eight species are endemic to the Qinghai-Tibetan Plateau (QTP). In this study, we employed both mitochondrial and nuclear loci to infer a molecular phylogeny for the eight QTP Gynaephora spp. We used the phylogeny to estimate divergence dates in a molecular dating analysis and to delimit species. This information allowed us to investigate associations between the diversification history of the eight QTP species and geological and climatic events. Phylogenetic analyses indicated that the eight QTP species formed a monophyletic group with strong supports in both Bayesian and maximum likelihood analyses. The low K2P genetic distances between the eight QTP species suggested that diversification occurred relatively quickly and recently. Out of the eight species, five species were highly supported as monophyletic, which were also recovered by species delimitation analyses. Samples of the remaining three species (G. aureata, G. rouergensis, and G. minora) mixed together, suggesting that further studies using extensive population sampling and comprehensive morphological approaches are necessary to clarify their species status. Divergence time estimation results demonstrated that the diversification and speciation of Gynaephora on the QTP began during the late Miocene/early Pliocene and was potentially affected by the QTP uplift and associated climate changes during this time.

Prolonged decay of molecular rate estimates for metazoan mitochondrial DNA

Evolutionary timescales can be estimated from genetic data using the molecular clock, often calibrated by fossil or geological evidence. However, estimates of molecular rates in mitochondrial DNA appear to scale negatively with the age of the clock calibration. Although such a pattern has been observed in a limited range of data sets, it has not been studied on a large scale in metazoans. In addition, there is uncertainty over the temporal extent of the time-dependent pattern in rate estimates. Here we present a meta-analysis of 239 rate estimates from metazoans, representing a range of timescales and taxonomic groups. We found evidence of time-dependent rates in both coding and non-coding mitochondrial markers, in every group of animals that we studied. The negative relationship between the estimated rate and time persisted across a much wider range of calibration times than previously suggested. This indicates that, over long time frames, purifying selection gives way to mutational saturation as the main driver of time-dependent biases in rate estimates. The results of our study stress the importance of accounting for time-dependent biases in estimating mitochondrial rates regardless of the timescale over which they are inferred.

Genome evolution in the primary endosymbiont of whiteflies sheds light on their divergence

Whiteflies are important agricultural insect pests, whose evolutionary success is related to a long-term association with a bacterial endosymbiont, Candidatus Portiera aleyrodidarum. To completely characterize this endosymbiont clade, we sequenced the genomes of three new Portiera strains covering the two extant whitefly subfamilies. Using endosymbiont and mitochondrial sequences we estimated the divergence dates in the clade and used these values to understand the molecular evolution of the endosymbiont coding sequences. Portiera genomes were maintained almost completely stable in gene order and gene content during more than 125 Myr of evolution, except in the Bemisia tabaci lineage. The ancestor had already lost the genetic information transfer autonomy but was able to participate in the synthesis of all essential amino acids and carotenoids. The time of divergence of the B. tabaci complex was much more recent than previous estimations. The recent divergence of biotypes B (MEAM1 species) and Q (MED species) suggests that they still could be considered strains of the same species. We have estimated the rates of evolution of Portiera genes, synonymous and nonsynonymous, and have detected significant differences among-lineages, with most Portiera lineages evolving very slowly. Although the nonsynonymous rates were much smaller than the synonymous, the genomic dN/dS ratios were similar, discarding selection as the driver of among-lineage variation. We suggest variation in mutation rate and generation time as the responsible factors. In conclusion, the slow evolutionary rates of Portiera may have contributed to its long-term association with whiteflies, avoiding its replacement by a novel and more efficient endosymbiont.

Beyond fossil calibrations: realities of molecular clock practices in evolutionary biology

Molecular-based divergence dating methods, or molecular clocks, are the primary neontological tool for estimating the temporal origins of clades. While the appropriate use of vertebrate fossils as external clock calibrations has stimulated heated discussions in the paleontological community, less attention has been given to the quality and implementation of other calibration types. In lieu of appropriate fossils, many studies rely on alternative sources of age constraints based on geological events, substitution rates and heterochronous sampling, as well as dates secondarily derived from previous analyses. To illustrate the breadth and frequency of calibration types currently employed, we conducted a literature survey of over 600 articles published from 2007 to 2013. Over half of all analyses implemented one or more fossil dates as constraints, followed by geological events and secondary calibrations (15% each). Vertebrate taxa were subjects in nearly half of all studies, while invertebrates and plants together accounted for 43%, followed by viruses, protists and fungi (3% each). Current patterns in calibration practices were disproportionate to the number of discussions on their proper use, particularly regarding plants and secondarily derived dates, which are both relatively neglected in methodological evaluations. Based on our survey, we provide a comprehensive overview of the latest approaches in clock calibration, and outline strengths and weaknesses associated with each. This critique should serve as a call to action for researchers across multiple communities, particularly those working on clades for which fossil records are poor, to develop their own guidelines regarding selection and implementation of alternative calibration types. This issue is particularly relevant now, as time-calibrated phylogenies are used for more than dating evolutionary origins, but often serve as the backbone of investigations into biogeography, diversity dynamics and rates of phenotypic evolution.

Is agriculture driving the diversification of the Bemisia tabaci species complex (Hemiptera: Sternorrhyncha: Aleyrodidae)?: Dating, diversification and biogeographic evidence revealed

BACKGROUND

Humans and insect herbivores are competing for the same food crops and have been for thousands of years. Despite considerable advances in crop pest management, losses due to insects remain considerable. The global homogenisation of agriculture has supported the range expansion of numerous insect pests and has been driven in part by human-assisted dispersal supported through rapid global trade and low-cost air passenger transport. One of these pests, is the whitefly, Bemisia tabaci, a cryptic species complex that contains some of the world's most damaging pests of agriculture. The complex shows considerable genetic diversity and strong phylogeographic relationships. One consequence of the considerable impact that members of the B. tabaci complex have on agriculture, is the view that human activity, particularly in relation to agricultural practices, such as use of insecticides, has driven the diversification found within the species complex. This has been particularly so in the case of two members of the complex, Middle East-Asia Minor 1 (MEAM1) and Mediterranean (MED), which have become globally distributed invasive species. An alternative hypothesis is that diversification is due to paleogeographic and paleoclimatological changes.

RESULTS

The idea that human activity is driving speciation within the B. tabaci complex has never been tested, but the increased interest in fossil whiteflies and the growth in molecular data have enabled us to apply a relaxed molecular clock and so estimate divergence dates for the major lineages within the B. tabaci species complex. The divergence estimates do not support the view that human activity has been a major driver of diversification.

CONCLUSIONS

Our analysis suggests that the major lineages within the complex arose approximately 60-30 mya and the highly invasive MED and MEAM1 split from the rest of the species complex around 12 mya well before the evolution of Homo sapiens and agriculture. Furthermore, the divergence dates coincide with a period of global diversification that occurred broadly across the plant and animal kingdoms and was most likely associated with major climatic and tectonic events.