BACKGROUND
ever since darwin, biologists have been interested in using phylogenetic trees to investigate the pattern and tempo of lineage diversification.
we see that the branching patterns in phylogenies contain information about the processes of speciation and extinction.
molecular phylogenies of species rich groups supply hypotheses of evolutionary relationships among higher taxa, and also allow us to test hypotheses about the tempo and pattern of diversification  <cit> , revealing both the topology of ancestor-descendant relationships and the tempo of descent among members of a clade.
it is clear that organismal diversification rates  have been variable across lineages and through time, with some lineages showing slowing rates of lineage accumulation  <cit>  and others showing rapid diversification  <cit>  or radiations.
radiations are generally defined in a way that includes rapid cladogenesis from a common ancestor  <cit> , yielding taxon rich clades.
the tempo and pattern of diversification of many major lineages of organisms, including insects, remain controversial  <cit>  - does macroevolution proceed at a relatively constant rate, accumulating lineages exponentially, or does it proceed through bursts of speciation triggered by processes such as adaptive radiation but otherwise remaining relatively constant?
rapid or explosive diversifications are characterized by lineages that have diverged in rapid succession within a relatively short time span.
molecular phylogenies can display a pattern that characterizes ancient rapid radiations as its signature  <cit> .
the signature can be described as close temporal spacing of a number cladogenetic or lineage-splitting events in a phylogeny, such that the internal branches  that link taxa together are very short.
although such a phylogenetic pattern is to be expected from an ancient rapid radiation, the pattern can also be observed in the case of other factors involved in phylogenetic reconstruction such as inadequate data, conflict within or among data sets, loss of phylogenetic signal over time, or data with inappropriate evolutionary rates  <cit> .
phylogenies of major lineages of insects based on morphological and/or molecular data have sometimes been contentious, often lacking the data to distinguish between alternative evolutionary relationships  <cit> .
the difficulty in confidently resolving basal splits in these phylogenies may be a result of insufficient time during cladogenesis to accumulate strong phylogenetic signal in the data because of rapid lineage splitting.
the ground beetle  subfamily harpalinae may display such a signature of rapid diversification.
harpalinae is the largest subfamily of carabid beetles and includes about  <dig>  species  <cit> , the bulk of the family's species-level diversity.
some harpalines exhibit a number of unusual morphological forms such as an elongated body in the genus agra, an ant-like form in calybe, snail-shell cracking mandibles in licinus  <cit> , and an extremely dorso-ventrally flattened body as in mormolyce.
not only are members of harpalinae diverse in morphological form, but also they display a variety of unusual lifestyles including granivory  <cit> , ovoviviparity  <cit> , symbiosis with ants and termites  <cit> , ectoparasitism of other insects  <cit> , specialized host mimicry by ectoparasites  <cit> , and arboreality  <cit> .
resolving the tribal relationships within the harpalines has been difficult, in part, because branch lengths at the base of the harpaline clade are very short in molecular based trees  and morphological phylogenetic analyses have not been able to resolve basal relationships within harpalinae  <cit> .
evidence from the fossil record that the subfamily harpalinae underwent an explosive radiation in the cretaceous period  <cit> , and the very short branches at the base of the harpaline clade in the molecular phylogeny of three nuclear genes point to the possibility that this speciose group represents a rapid radiation within carabids, and the tribes may have originated almost simultaneously.
harpalines contain more than  <dig> times the number of species as brachininae , their sister group and almost  <dig> times more than the next closest clade, the austral psydrines   <cit> .
rapid, ancient radiations are a difficult challenge to phylogenetic inference  <cit>  because short interior branches limit the historical record of early diversification.
in harpalinae, there are also long terminal branches  <cit>  that are prone to the analytical artifact of long-branch attraction  <cit> .
using the preferred phylogeny of harpalines and close relatives  <cit> , molecular sequence data  <cit> , and fossil data, we investigated the timing and tempo of the diversification of harpalines.
the aim of this work was to test the hypothesis that harpalines evolved their incredible diversity through an explosive radiation shortly after their evolutionary origin.
our goals were to determine  the age of divergence of the major lineages  of harpalines,  whether the evolutionary history of diversification of the subfamily harpalinae and its brachinine sister group depart significantly from a constant rate model, and  whether harpalines and/or brachinines have experienced changes in diversification rates during their evolutionary history indicative of an explosive radiation of lineages.
in this study, we report on the timing and patterns of evolution in the largest subfamily of ground beetles and its sister group, the brachinine bombardier beetles.
methods
phylogenetic tree and molecular sequence data
using the 28sbe+wg maximum likelihood tree topology from ober and maddison  <cit>  and dna sequence data sets for 18s from ober  <cit>  and 28s and wingless from ober and maddison  <cit> , we estimated branch lengths under the maximum likelihood gtr+Γ+i model in paup*  <cit> .
branch lengths for the all data combined tree were estimated from the 28s data due to the large amount of missing data for many taxa.
the 28sbe+wg likelihood tree  <cit>  represents the best phylogenetic hypothesis available for harpalinae.
it has the most extensive taxon sampling for the subfamily and uses a model of molecular evolution to estimate the phylogeny from two nuclear genes.
monophyly of harpalinae is based on molecular data  <cit> , morphological data  <cit> , defensive chemical data  <cit> , and cytological data  <cit> .
most tribes within the subfamily are monophyletic, but the relationships among tribes are not strongly supported  <cit> .
the effects of phylogenetic uncertainty on divergence dates were not explored in depth in this study.
one beast analysis was done allowing tree topology to be estimated along with other parameter values.
the topology was largely congruent with the 28sbe+wg maximum likelihood tree topology from ober and maddison  <cit>  with a few exceptions noted in figure  <dig>  a complete list of taxa and dna sequence data used in this study from ober  <cit>  and ober and maddison  <cit>  is in additional file  <dig>  table s <dig>  we estimated divergence dates from combined 18s, 28s and wg data and also from each gene separately.
for separate gene analyses, taxa for which no molecular data was available were pruned from the tree yielding a 18s tree with  <dig> harpaline taxa and  <dig> outgroups, a wg tree with  <dig> harpalines and  <dig> outgroups, and a 28s tree with  <dig> harpalines and  <dig> outgroups.
we used a likelihood ratio test to look for evidence of a molecular clock in the molecular phylogenies.
divergence time estimation
divergence times were estimated using the penalized likelihood approach provided by the program r8s v. <dig>   <cit>  and a bayesian relaxed clock method in beast v. <dig> .
<dig>  <cit>  for all data combined and separate gene datasets.
for r8s and beast analyses, we calibrated several internal nodes  based on the fossil record .
greek symbols correspond to fossil-based minimum age constraints placed on the phylogeny in figure  <dig>
the r8s penalized likelihood method with the tn algorithm uses semiparametric rate-smoothing to optimize rate changes among branches.
smoothing parameters were estimated by using the fossil-based model cross-validation feature in r8s  <cit> .
all  <dig> fossils passed the cross-validation test and were included in the analysis.
the outermost outgroup was pruned from the tree.
two r8s analyses were done for each data set, one with the younger fossil dates for each calibration point and one with the older dates.
all nodes ages were free to be estimated.
all nodes for which there was fossil information were constrained to be no older than a maximum age of  <dig>  or  <dig>  my ago, the fossil age of carabidae  <cit> , except the root node of the tree which was constrained with a maximum age of  <dig> my ago, the age of adephaga  <cit> , with no minimum age constraint.
minimum age constraints were applied to nodes according to fossil information .
we also inferred divergence dates using a bayesian relaxed clock uncorrelated lognormal method in beast for all data combined and the three genes separately.
for the all data combined analysis, we concatenated the three data sets for all taxa, but created three data partitions for 18s, 28s, and wg.
we chose a separate gtr+i+Γ model with four gamma categories, based on mrmodeltest  <cit> ; a yule process of speciation, and an uncorrelated lognormal relaxed clock model of rate variation for each partition.
model parameters were unlinked across partitions.
for the separate gene analyses, we applied a gtr+i+Γ model with four gamma categories, a yule process of speciation, and an uncorrelated lognormal relaxed clock model of rate.
the tree topology prior was fixed for all beast analyses except one all-data-combined analysis where the tree topology was estimated along with other parameter values.
the same calibration points as in the r8s analyses were used .
node constraints were assigned a normal prior distribution with means equal to the midpoint of the fossil date and the standard deviations encompassing the youngest and oldest age of each fossil.
a normal distribution was chosen because it allows uncertainty in the calibration estimates  <cit> .
after an initial period of fine-tuning the operators, two separate mcmc analyses were run for  <dig> million generations for each data set with parameters sampled every  <dig> generations.
independent runs were combined using logcombiner <dig> .
<dig>  <cit> , and the first 20% of the generations from each run was discarded as burnin.
convergence of the chains was checked using tracer  <dig>   <cit> .
the searches achieved adequate mixing as assessed by the high effective sampling size  values for all parameters of  <dig> or greater.
node ages and upper and lower bounds of the 95% highest posterior density interval for divergence times was calculated using treeannotator  <dig> .
<dig> and visualized using figtree  <dig> .
<dig>  <cit> .
biodiversity estimate
to test if the harpalinae clade  is more diverse than its sister clade, the brachininae , we determined if the two groups had significantly different number of extant species using a conservative test of slowinsky and guyer  <cit> .
diversification rates
several diversification statistics were run on the chronograms for each clade  and for the overall tree.
all diversification tests were performed in r using ape  <cit> , geiger  <cit> , and laser  <cit> .
the overall rates of net diversification for the harpalinae clade and the brachininae clade were estimated from the pruned all-data-combined beast and r8s chronograms using geiger based on a pure-birth process of diversification and also with a speciation:extinction rate of 2: <dig> using the method of magallón and sanderson  <cit> .
the rate estimate took into account the extant taxa missing from each clade not sampled in the chronograms.
the relative cladogenesis  statistic, which uses a broken-stick distribution to identify branches in the tree ancestral to a greater proportion of extant descendant lineages than expected by chance  <cit> , was used to calculate the probability that a particular lineage at time t will have k descendants given the total number of descendants at time  <dig> .
this test detects and locates unusually rapid shifts in diversification rates by looking for unusually lineage-rich clades in the cohort of clades all originating at one slice of time in the tree.
we estimated the rc statistic with bonferroni correction as implemented in geiger for the all-data-combined r8s and beast chronograms.
rates of cladogenesis through time were investigated in the harpalinae and brachininae all data combined beast chronograms using constant rate  test of pybus and harvey  <cit>  that estimates the gamma statistic  of a given chronogram.
under a yule  process, γ values of completely sampled phylogenies have been shown to fit a standard normal distribution with mean =  <dig>  <cit> .
significantly negative values of γ  are indicative of decreasing rates of cladogenesis through time .
however, incomplete taxon sampling has been shown to inflate the type i error of the cr test  <cit> .
to correct for the under sampling in our analysis , we used monte carlo constant rate  test  <cit>  in which full topologies for each clade were simulated under the yule process and then randomly sub sampled in mesquite  <cit>  to generate a corrected null distribution.
we compared our observed γ for brachinines and harpalines to the appropriate null distribution based on  <dig> simulated trees.
we calculated γ using the laser package in r for each of these simulated trees and then used the resulting set of values as our null distribution under incomplete taxon sampling  <cit> .
for both the mccr tests, we performed a one-tailed test with a critical value of p =  <dig> .
rates of diversification were evaluated using lineage through time  plots, which illustrate the accumulation of lineages over time for the all data combined beast chronograms for the ingroup harpalinae and its sister group brachininae using ape.
to evaluate the effects of incomplete taxon sampling on the slope of the ltt plot we used  <dig> subsampled simulated yule process brachinine and harpaline trees to construct mean ltt curves for comparison with the empirical ltt curves.
we evaluate the fit of six models of diversification  <cit>  in laser using akaike information criteria  and/or hierarchical likelihood ratio tests  for the full all-data-combined beast chronogram, the harpaline clade and the brachinine clade.
model  <dig>  assumes a pure birth or yule model of diversification, model  <dig>  assumes a constant birth death model of diversification, model  <dig>  specifies a variable speciation rate and constant extinction rate, model  <dig>  specifies a variable extinction rate and constant speciation rate, model  <dig>  assumes an abrupt change in diversification rate at some breakpoint in the past, model  <dig>  assumes three different diversification rates with two breakpoints in the past.
for models  <dig> and  <dig> we evaluated the hypothesized shifts in diversification rate every one million years.
RESULTS
divergence time estimations of harpalinae
all three molecular datasets rejected a strict molecular clock for the tree .
phylogenetic relationships of tribes within harpalines estimated by beast from all data combined were largely congruent with the 28sbe+wg maximum likelihood tree topology from ober and maddison  <cit> .
for the most part, tribes and supertribes were monophyletic except for perigonini.
however, some of the sistergroup relationships of tribes differed between the ober and maddison  <cit>  and beast tree .
the difference in tree topology did not affect the estimate of the age of the origin of harpalinae , but the origin of most of the tribes was estimated to be slightly older than in the ober and maddison  <cit>  tree.
nonetheless, by the eocene,  <dig> my after the origin of the subfamily, all tribes are present, similar to the timing of diversification estimated from the ober and maddison  <cit>  tree.
our estimate for the origin of harpalinae ranged from about  <dig> my ago to about  <dig> my ago  depending on the type of analysis and the gene used to estimate the divergence times.
r8s estimates of divergence times were generally older than the beast estimates, outside the 95% confidence intervals for the beast estimates.
r8s divergence date estimates for the origin of harpalines from all dna data sets were much older  than to the fossil record .
results of divergence date estimates from individual genes showed that the oldest dates for all nodes were inferred from the wg gene and the youngest dates for all nodes were inferred from the 18s gene with r8s .
the results of the beast analyses did not show a clear bias of particular genes for older or younger date estimates .
estimated age of stem lineage origins within harpalinae from all-data-combined and18s rdna, wingless, and 28s rdna data sets from r8s relaxed clock analyses using the youngest hypothesized age of the fossil calibration  and the oldest hypothesized age of the fossil calibration  from table  <dig>  young and old ages resulted in essentially identical estimates for the wg gene.
not all nodes were present in the 18s data set due to limited taxon sampling  <cit> .
the capital letters correspond to clades in figure  <dig>
estimated age of stem lineage origins within harpalinae from the all-data-combined and separate 18s rdna, wingless, and 28s rdna data sets from beast analyses.
not all nodes were present in the 18s data set due to limited taxon sampling  <cit> .
the capital letters correspond to clades in figure  <dig>
the combined molecular data suggest harpalines originated approximately  <dig> my ago , but may not have diversified much until ~ <dig> my after the origin of the subfamily.
our analyses suggest that many harpaline tribes probably originated ~ <dig> my ago with most of the extant genera evolving  <dig> -  <dig> my ago.
the average age estimates of most of the tribes imply that they originated within a  <dig> my time window  about  <dig> my after the origin of harpalinae .
approximately  <dig> my elapsed from the origin of harpalines to the time when all major lineages and tribes are present.
tempo of diversification
although the sister groups are the same age and descended from a common speciation event, there is a significant difference  in species number between harpalinae and brachininae, indicating an increase in net diversification  within the harpaline clade.
the clade that contains harpalines is also significantly more diverse  than the sister group to brachinines + harpalines, the austral psydrines, which contains only  <dig> extant species .
a conservative slowinsk and guyer test, to determine if clades differ in number of species, shows harpalinae with significantly more species than brachininae .
the estimated overall net diversification rate with a pure birth model of diversification in brachinines was  <dig>  per my  with the beast all-data-combined chronogram to  <dig>  per my with the r8s all data combined chronogram, and in harpalines the estimated overall net diversification rate with a pure birth model was  <dig>  per my  with the beast chronogram to  <dig>  per my with the r8s chronogram.
the rc test indicated a significant difference in cladogenesis rate between harpalines and brachinines.
we located unusually rapid shifts in diversification in the node subtending harpalinae and along at least four basal branches of harpalines .
there were also significant increases in diversification within and the tribe lebiini, especially the subtribe calleidina .
results of the cr and mccr tests showed the observed value of γ in the harpalinae clade was negative , which rejects the hypothesis that rates of lineage accumulation remained constant over time, however, when compared to the distribution of γ statistics of the simulations, the harpaline γ statistic was not statistically significant .
the brachinine γ statistic was not significantly different than zero, and the cr and mccr tests do not reject a constant rate of diversification .
the semilogarithmic ltt plots for the harpaline and brachinine clades are shown in figure  <dig>  the accumulation of harpaline lineages did not differ significantly from a constant rate of diversification early in harpaline evolution when compared to  <dig> simulated yule process trees .
however, there was a steep increase in the number of harpaline lineages early after the origin of the subfamily, approximately at  <dig> my ago, and a shift in rate at about  <dig> my ago .
this was the time period when most harpaline tribes were diversifying according to the beast divergence time estimates.
brachinines showed a fairly slow but constant rate of lineage accumulation , within or slightly below the range of simulated trees.
for the harpaline clade, a model of two abrupt rate shifts, model 6: yule3rate , was chosen over other constant and variable rate models .
a diversification rate decrease  was estimated to have happened first at  <dig>  my ago and then again at  <dig>  my ago  .
statistical analyses of the ltt plot for brachinines indicated that a pure birth model of constant rate of diversification  could not be rejected .
a model of two abrupt rate shifts, model 6: yule3rate  best fit the overall tree with decreases in diversification rates occurring at  <dig>  my ago and  <dig>  my ago .
pure birth assumes a pure birth or yule model of diversification, birth death assumes a constant birth death model of diversification, spvar specifies a variable speciation rate and constant extinction rate, exvar specifies a variable extinction rate and constant speciation rate, yule2rate assumes an abrupt change in diversification rate at some breakpoint in the past, yule3rate assumes three different diversification rates with two breakpoints in the past.
a breakpoint at  <dig>  my ago
b breakpoints at  <dig>  and  <dig>  my ago
c breakpoint at  <dig>  my ago
d breakpoints at  <dig>  and  <dig>  my ago
e breakpoint at  <dig>  my ago
f breakpoints at  <dig>  and  <dig>  my ago
discussion
divergence time estimations of harpalinae
we found a mid-cretaceous origin of harpalinae with most tribes arising by the paleocene.
crown diversification of tribes appears to have commenced by the late-middle eocene, and nearly all tribes appeared by the late eocene.
in a relatively short amount of time , harpalines diversified from just a few lineages into almost all the major lineages extant today.
from the divergence date estimates, it appears this rapid diversification was not immediately after the origin of harpalines, but about  <dig> my later.
beast estimates of divergence dates from all molecular data combined were slightly older than the fossil record.
r8s estimated divergence dates much older that the fossil record  and the beast estimates .
differences between the r8s and beast estimates may due to the relaxed clock method used to infer dates and the prior distributions for the divergence time of calibration points in beast versus the minimum and maximum age constraints in r8s.
we prefer the dates from the all data combined beast analyses as a more accurate estimate of harpaline evolution based on the partitioned model of molecular evolution and bayesian framework.
estimates of divergence dates from molecular data are subject to several sources of error.
estimates may be subject to errors in fossil calibration points or topology, inappropriate priors, or invalid clock assumptions, including much more rate heterogeneity than can be accommodated by the programs we used.
the error introduced by finite sequence data or model misspecification, such as neglect of multiple hits and across site rate variation  <cit>  can cause spurious results in divergence time estimates.
sources of error have been reviewed elsewhere .
we did not explicitly explore the effect of differences in tree topology on divergence time estimates.
we assumed the ober and maddison  <cit>  28s+wg maximum likelihood tree was the best phylogenetic hypothesis for harpalines available, but results from the beast analysis of all data combined where topology was also estimated, did not change our conclusions about the timing of harpaline origin or the rapidness of diversification of tribes.
our estimates of harpaline divergence dates are a working hypothesis subject to the addition of more taxa, molecular data, improved phylogenetic hypotheses, and fossils yet to be discovered.
tempo of diversification
harpalines have had an increase in net diversification compared to their sister group, the brachinines.
it is difficult to tease apart whether this is due to increased speciation or reduced extinction, because the net diversification rate is simply the difference between speciation and extinction rates.
we tested for a change in extinction rate through time in harpalines and brachinines, and a change in extinction rate does not explain the difference in species number between the sister groups because he estimated extinction rate for both is the same .
the net diversification rate for harpalines  is higher than the rate estimated for explosive diversification of angiosperms   <cit>  and on par with the rapid mammalian radiation   <cit> .
an increase in diversification is evident in the rc test where several harpaline clades have experienced higher than usual cladogenesis.
the time of increased cladogenesis occurs approximately  <dig> -  <dig> my ago, during a time of incredible faunal and floral turnover in many plant and animal groups.
in contrast, there were not shifts in diversification in brachinines or austral psydrines.
the rc test is sensitive to temporal depth, phylogenetic scope, and nonindependence of diversification rate shifts  <cit> .
therefore rate shifts identified by rc should be interpreted with caution, especially nested shifts.
if a clade is especially diverse, then one or more of its parent clades will also show significant diversification through a "trickle-down" effect  <cit> .
this may be what is being observed along the basal braches of harpalinae, and a conservative interpretation of the rc results would be to attribute the increase in diversification rate to the shallowest significant branch .
in this case, there are significant increases in diversification at the base of the harpaline subfamily and within the tribe lebiini, the largest tribe of harpalines with over  <dig>  species.
alternatively, a more complex series of radiations may have taken place sequentially.
an increase in diversification could have taken place over a longer period of time and shifts in diversification could have occurred over several nodes or branches.
evidence of shifts in diversification rates in harpalines is harder to interpret from the cr and mccr tests and the ltt plots.
older lineages have a higher risk of being extinct, and thus not sampled, than younger ones .
incomplete taxon sampling can give the appearance of nodes near the root of the tree giving rise to more extant descendants than nodes near the tips and are therefore more likely to be included in a small random sample   <cit> .
the mccr test for harpalines did not indicate that nodes were closer to the root than expected under a constant rates model after correcting for incomplete taxon sampling.
the taxa included in our phylogenetic tree did not represent a random subsample of the total harpaline species, as assumed by the mccr test with yule process simulations.
instead, they most likely represent an oversampling of deeper lineages and an undersampling of closely related lineages, which may have affected the conclusions drawn from this test.
furthermore the mccr test suffers from low statistical power  <cit>  and has been shown to have decreased sensitivity in detecting an initial high rate of diversification followed by a shift to a lower rate  <cit> .
the harpaline ltt plot with simulated constant rate trees agreed with the results of the mccr test in that it did not show significantly different rate shifts from trees simulated under a yule process, but the harpaline ltt plot showed a decrease in diversification rate after about  <dig> my ago, the period after the origin of tribes.
although the simulated trees diversified at a constant rate , randomly removing taxa generated an apparent deceleration of lineage accumulation .
as with the mccr test, our tree is not a random subsample of harpaline taxa.
this bias may affect how the ltt plot is interpreted.
the ltt plot and associated confidence interval generated by simulations is just one method by which timing and tempo of diversification can be visualized and evaluated, other lines of evidence suggest a fairly rapid accumulation of harpaline lineages in the paleocene and eocene.
the shift in diversification rate from a high rate initially to lower rates later in harpaline evolutionary history are evident in the results of the tests of models of diversification .
an increase in extinction rate later in harpaline evolutionary history does not explain the pattern of high diversification early and declining diversification through time.
the pattern of lineage accumulation through time rejected a model of changing extinction rate.
the best fitting model for harpalines shows a decrease in rate at  <dig> my ago and again at about  <dig> my ago.
the first rate decrease corresponds well to the end of tribe diversification estimated in beast, and perhaps the second shift, which is much more subtle signals the end of a burst of generic diversification, but more taxa will need to be included to test this hypothesis.
estimated dates of divergence, the rc test, evidence from the ltt plot, and models of diversification tests, seem to suggest a late cretaceous through middle eocene period of rapid diversification in harpalines.
the paleogene was a time of profound reorganization of biota, perhaps in part due to a period of warm "greenhouse" climate during the paleocene/eocene thermal maximum or continental fission during the late cretaceous.
the late cretaceous to middle eocene is the same time period when many other plant and animal taxa are experiencing tremendous changes in evolution from mass extinctions at the k-t boundary to explosive radiations in birds  <cit> , mammals  <cit>  and some plants  <cit>  and insects  <cit> .
harpaline diversification occurred just after the angiosperm cretaceous radiation  <dig> -  <dig> my ago, when flowering plants achieved widespread floristic dominance for the first time  <cit> .
it does not appear that harpalines necessarily co-evolved with angiosperms, but the angiosperm radiation, along with other biotic and abiotic factors, during this period provided harpalines with the possibility of novel ecological interactions.
although the vast majority of harpalines are predatory, the evolving complex angiosperm habitats may have afforded harpalines with new niches in which to diversify, especially because some major groups of harpalines, like the species rich tribe lebiini, evolved arboreal lifestyles instead of ground dwelling  <cit> .
recent studies have proposed that another species rich predatory insect group, ants, may have also radiated primarily in the late cretaceous  <cit>  or paleogene  <cit> .
CONCLUSIONS
harpaline species diversity is remarkable compared to its relatives and other carabid clades.
their high rates of net diversification have made them an evolutionary success story and an important part of many modern ecosystems.
from the dates of divergence estimates, there was a large amount of rapid speciation in a short amount of time in the late cretaceous and early paleogene, but it does not appear that it was a particularly explosive radiation immediately after their evolutionary origin.
instead, harpalines underwent a somewhat elevated bursts of diversification about  <dig> my after their origin.
this suggests that harpalines may have had a long "phylogenetic fuse"  <cit>  before tribes and major lineages began to diversify as seen in birds  <cit>  and mammals  <cit> .
relationships of tribes and within the carabid subfamily harpalinae are not strongly supported  <cit> .
the lack of support from the data for deep branches in the harpaline molecular phylogeny could be explained by a rapid radiation or cladogenesis during the evolutionary history of this group.
short internal branches created by rapid cladogenesis leave a meager record of diversification that is further potentially obscured by long terminal branches.
this explanation was supported in a study attempting to resolve early metazoan evolution  <cit>  and supported by fiala and sokal's  <cit>  simulations.
they suggested that short internal branches and long terminal branches reduces the accuracy with which a phylogeny can be estimated, hence the resolution that can be obtained.
the pattern seen in the molecular phylogenetic trees of short internal nodes at the base of harpalines with longer terminal branches  <cit>  can be explained in several ways.
it may be the case that harpalines did undergo a burst of diversification during the late cretaceous and early paleogene, but the taxon sampling in the molecular phylogeny was not complete enough to distinguish the pattern clearly.
incomplete taxon sampling in such a large group makes it difficult to interpret the results of the constant rate tests and ltt plots as a significant burst of diversification in harpalinae evolutionary history.
it may be the case that harpalines did diversify rapidly, but it was not as explosive as carabidologists have assumed.
on the other hand, the tempo of harpaline evolution may have been fairly constant though time with a high , but steady rate of diversification and lineage accumulation, leading to the large and speciose clade we see today.
alternatively, the molecular phylogenies inferred from the three nuclear genes with short basal branches may be due to the information content and phylogenetic utility of the genes or systematic errors in selecting models of evolution for phylogenetic inference.
in any case, the inclusion of additional harpaline and brachinine taxa and new molecular data would help shed light on the evolutionary history of the harpalinae and the patterns and tempo of diversification of this group of beetles.
authors' contributions
kao conceived of the study, carried out some of the analyses, and drafted the manuscript.
tnh participated in the design of the study and performed many of the analyses.
both authors read and approved the final manuscript.
supplementary material
additional file 1
table s <dig>  taxa and genbank numbers used in this study.
click here for file
acknowledgements
duane mckenna, jessica thomas, and three anonymous reviewers generously provided comments that greatly improved the manuscript.
the authors thank the college of the holy cross and richard b. fisher summer research fellowship for funding for this project.

 ever since darwin, biologists have been interested in using phylogenetic trees to investigate the pattern and tempo of lineage diversification. we see that the branching patterns in phylogenies contain information about the processes of speciation and extinction. molecular phylogenies of species rich groups supply hypotheses of evolutionary relationships among higher taxa, and also allow us to test hypotheses about the tempo and pattern of diversification, revealing both the topology of ancestor-descendant relationships and the tempo of descent among members of a clade. it is clear that organismal diversification rates have been variable across lineages and through time, with some lineages showing slowing rates of lineage accumulation and others showing rapid diversification or radiations. radiations are generally defined in a way that includes rapid cladogenesis from a common ancestor, yielding taxon rich clades. the tempo and pattern of diversification of many major lineages of organisms, including insects, remain controversial - does macroevolution proceed at a relatively constant rate, accumulating lineages exponentially, or does it proceed through bursts of speciation triggered by processes such as adaptive radiation but otherwise remaining relatively constant? rapid or explosive diversifications are characterized by lineages that have diverged in rapid succession within a relatively short time span. molecular phylogenies can display a pattern that characterizes ancient rapid radiations as its signature. the signature can be described as close temporal spacing of a number cladogenetic or lineage-splitting events in a phylogeny, such that the internal branches that link taxa together are very short. although such a phylogenetic pattern is to be expected from an ancient rapid radiation, the pattern can also be observed in the case of other factors involved in phylogenetic reconstruction such as inadequate data, conflict within or among data sets, loss of phylogenetic signal over time, or data with inappropriate evolutionary rates. phylogenies of major lineages of insects based on morphological and/or molecular data have sometimes been contentious, often lacking the data to distinguish between alternative evolutionary relationships. the difficulty in confidently resolving basal splits in these phylogenies may be a result of insufficient time during cladogenesis to accumulate strong phylogenetic signal in the data because of rapid lineage splitting. the ground beetle subfamily harpalinae may display such a signature of rapid diversification. harpalinae is the largest subfamily of carabid beetles and includes about species, the bulk of the family's species-level diversity. some harpalines exhibit a number of unusual morphological forms such as an elongated body in the genus agra, an ant-like form in calybe, snail-shell cracking mandibles in licinus, and an extremely dorso-ventrally flattened body as in mormolyce. not only are members of harpalinae diverse in morphological form, but also they display a variety of unusual lifestyles including granivory, ovoviviparity, symbiosis with ants and termites, ectoparasitism of other insects, specialized host mimicry by ectoparasites, and arboreality. resolving the tribal relationships within the harpalines has been difficult, in part, because branch lengths at the base of the harpaline clade are very short in molecular based trees and morphological phylogenetic analyses have not been able to resolve basal relationships within harpalinae. evidence from the fossil record that the subfamily harpalinae underwent an explosive radiation in the cretaceous period, and the very short branches at the base of the harpaline clade in the molecular phylogeny of three nuclear genes point to the possibility that this speciose group represents a rapid radiation within carabids, and the tribes may have originated almost simultaneously. harpalines contain more than times the number of species as brachininae, their sister group and almost times more than the next closest clade, the austral psydrines. rapid, ancient radiations are a difficult challenge to phylogenetic inference because short interior branches limit the historical record of early diversification. in harpalinae, there are also long terminal branches that are prone to the analytical artifact of long-branch attraction. using the preferred phylogeny of harpalines and close relatives, molecular sequence data, and fossil data, we investigated the timing and tempo of the diversification of harpalines. the aim of this work was to test the hypothesis that harpalines evolved their incredible diversity through an explosive radiation shortly after their evolutionary origin. our goals were to determine the age of divergence of the major lineages of harpalines, whether the evolutionary history of diversification of the subfamily harpalinae and its brachinine sister group depart significantly from a constant rate model, and whether harpalines and/or brachinines have experienced changes in diversification rates during their evolutionary history indicative of an explosive radiation of lineages. in this study, we report on the timing and patterns of evolution in the largest subfamily of ground beetles and its sister group, the brachinine bombardier beetles. methods phylogenetic tree and molecular sequence data using the 28sbe+wg maximum likelihood tree topology from ober and maddison and dna sequence data sets for 18s from ober and 28s and wingless from ober and maddison, we estimated branch lengths under the maximum likelihood gtr+Γ+i model in paup*. branch lengths for the all data combined tree were estimated from the 28s data due to the large amount of missing data for many taxa. the 28sbe+wg likelihood tree represents the best phylogenetic hypothesis available for harpalinae. it has the most extensive taxon sampling for the subfamily and uses a model of molecular evolution to estimate the phylogeny from two nuclear genes. monophyly of harpalinae is based on molecular data, morphological data, defensive chemical data, and cytological data. most tribes within the subfamily are monophyletic, but the relationships among tribes are not strongly supported. the effects of phylogenetic uncertainty on divergence dates were not explored in depth in this study. one beast analysis was done allowing tree topology to be estimated along with other parameter values. the topology was largely congruent with the 28sbe+wg maximum likelihood tree topology from ober and maddison with a few exceptions noted in figure a complete list of taxa and dna sequence data used in this study from ober and ober and maddison is in additional file table s we estimated divergence dates from combined 18s, 28s and wg data and also from each gene separately. for separate gene analyses, taxa for which no molecular data was available were pruned from the tree yielding a 18s tree with harpaline taxa and outgroups, a wg tree with harpalines and outgroups, and a 28s tree with harpalines and outgroups. we used a likelihood ratio test to look for evidence of a molecular clock in the molecular phylogenies. divergence time estimation divergence times were estimated using the penalized likelihood approach provided by the program r8s v. and a bayesian relaxed clock method in beast v.. for all data combined and separate gene datasets. for r8s and beast analyses, we calibrated several internal nodes based on the fossil record. greek symbols correspond to fossil-based minimum age constraints placed on the phylogeny in figure the r8s penalized likelihood method with the tn algorithm uses semiparametric rate-smoothing to optimize rate changes among branches. smoothing parameters were estimated by using the fossil-based model cross-validation feature in r8s. all fossils passed the cross-validation test and were included in the analysis. the outermost outgroup was pruned from the tree. two r8s analyses were done for each data set, one with the younger fossil dates for each calibration point and one with the older dates. all nodes ages were free to be estimated. all nodes for which there was fossil information were constrained to be no older than a maximum age of or my ago, the fossil age of carabidae, except the root node of the tree which was constrained with a maximum age of my ago, the age of adephaga, with no minimum age constraint. minimum age constraints were applied to nodes according to fossil information. we also inferred divergence dates using a bayesian relaxed clock uncorrelated lognormal method in beast for all data combined and the three genes separately. for the all data combined analysis, we concatenated the three data sets for all taxa, but created three data partitions for 18s, 28s, and wg. we chose a separate gtr+i+Γ model with four gamma categories, based on mrmodeltest; a yule process of speciation, and an uncorrelated lognormal relaxed clock model of rate variation for each partition. model parameters were unlinked across partitions. for the separate gene analyses, we applied a gtr+i+Γ model with four gamma categories, a yule process of speciation, and an uncorrelated lognormal relaxed clock model of rate. the tree topology prior was fixed for all beast analyses except one all-data-combined analysis where the tree topology was estimated along with other parameter values. the same calibration points as in the r8s analyses were used. node constraints were assigned a normal prior distribution with means equal to the midpoint of the fossil date and the standard deviations encompassing the youngest and oldest age of each fossil. a normal distribution was chosen because it allows uncertainty in the calibration estimates. after an initial period of fine-tuning the operators, two separate mcmc analyses were run for million generations for each data set with parameters sampled every generations. independent runs were combined using logcombiner., and the first 20% of the generations from each run was discarded as burnin. convergence of the chains was checked using tracer. the searches achieved adequate mixing as assessed by the high effective sampling size values for all parameters of or greater. node ages and upper and lower bounds of the 95% highest posterior density interval for divergence times was calculated using treeannotator. and visualized using figtree.. biodiversity estimate to test if the harpalinae clade is more diverse than its sister clade, the brachininae, we determined if the two groups had significantly different number of extant species using a conservative test of slowinsky and guyer. diversification rates several diversification statistics were run on the chronograms for each clade and for the overall tree. all diversification tests were performed in r using ape, geiger, and laser. the overall rates of net diversification for the harpalinae clade and the brachininae clade were estimated from the pruned all-data-combined beast and r8s chronograms using geiger based on a pure-birth process of diversification and also with a speciation:extinction rate of 2: using the method of magallón and sanderson. the rate estimate took into account the extant taxa missing from each clade not sampled in the chronograms. the relative cladogenesis statistic, which uses a broken-stick distribution to identify branches in the tree ancestral to a greater proportion of extant descendant lineages than expected by chance, was used to calculate the probability that a particular lineage at time t will have k descendants given the total number of descendants at time. this test detects and locates unusually rapid shifts in diversification rates by looking for unusually lineage-rich clades in the cohort of clades all originating at one slice of time in the tree. we estimated the rc statistic with bonferroni correction as implemented in geiger for the all-data-combined r8s and beast chronograms. rates of cladogenesis through time were investigated in the harpalinae and brachininae all data combined beast chronograms using constant rate test of pybus and harvey that estimates the gamma statistic of a given chronogram. under a yule process, γ values of completely sampled phylogenies have been shown to fit a standard normal distribution with mean =. significantly negative values of γ are indicative of decreasing rates of cladogenesis through time. however, incomplete taxon sampling has been shown to inflate the type i error of the cr test. to correct for the under sampling in our analysis, we used monte carlo constant rate test in which full topologies for each clade were simulated under the yule process and then randomly sub sampled in mesquite to generate a corrected null distribution. we compared our observed γ for brachinines and harpalines to the appropriate null distribution based on simulated trees. we calculated γ using the laser package in r for each of these simulated trees and then used the resulting set of values as our null distribution under incomplete taxon sampling. for both the mccr tests, we performed a one-tailed test with a critical value of p =. rates of diversification were evaluated using lineage through time plots, which illustrate the accumulation of lineages over time for the all data combined beast chronograms for the ingroup harpalinae and its sister group brachininae using ape. to evaluate the effects of incomplete taxon sampling on the slope of the ltt plot we used subsampled simulated yule process brachinine and harpaline trees to construct mean ltt curves for comparison with the empirical ltt curves. we evaluate the fit of six models of diversification in laser using akaike information criteria and/or hierarchical likelihood ratio tests for the full all-data-combined beast chronogram, the harpaline clade and the brachinine clade. model assumes a pure birth or yule model of diversification, model assumes a constant birth death model of diversification, model specifies a variable speciation rate and constant extinction rate, model specifies a variable extinction rate and constant speciation rate, model assumes an abrupt change in diversification rate at some breakpoint in the past, model assumes three different diversification rates with two breakpoints in the past. for models and we evaluated the hypothesized shifts in diversification rate every one million years.  divergence time estimations of harpalinae all three molecular datasets rejected a strict molecular clock for the tree. phylogenetic relationships of tribes within harpalines estimated by beast from all data combined were largely congruent with the 28sbe+wg maximum likelihood tree topology from ober and maddison. for the most part, tribes and supertribes were monophyletic except for perigonini. however, some of the sistergroup relationships of tribes differed between the ober and maddison and beast tree. the difference in tree topology did not affect the estimate of the age of the origin of harpalinae, but the origin of most of the tribes was estimated to be slightly older than in the ober and maddison tree. nonetheless, by the eocene, my after the origin of the subfamily, all tribes are present, similar to the timing of diversification estimated from the ober and maddison tree. our estimate for the origin of harpalinae ranged from about my ago to about my ago depending on the type of analysis and the gene used to estimate the divergence times. r8s estimates of divergence times were generally older than the beast estimates, outside the 95% confidence intervals for the beast estimates. r8s divergence date estimates for the origin of harpalines from all dna data sets were much older than to the fossil record. results of divergence date estimates from individual genes showed that the oldest dates for all nodes were inferred from the wg gene and the youngest dates for all nodes were inferred from the 18s gene with r8s. the results of the beast analyses did not show a clear bias of particular genes for older or younger date estimates. estimated age of stem lineage origins within harpalinae from all-data-combined and18s rdna, wingless, and 28s rdna data sets from r8s relaxed clock analyses using the youngest hypothesized age of the fossil calibration and the oldest hypothesized age of the fossil calibration from table young and old ages resulted in essentially identical estimates for the wg gene. not all nodes were present in the 18s data set due to limited taxon sampling. the capital letters correspond to clades in figure estimated age of stem lineage origins within harpalinae from the all-data-combined and separate 18s rdna, wingless, and 28s rdna data sets from beast analyses. not all nodes were present in the 18s data set due to limited taxon sampling. the capital letters correspond to clades in figure the combined molecular data suggest harpalines originated approximately my ago, but may not have diversified much until ~ my after the origin of the subfamily. our analyses suggest that many harpaline tribes probably originated ~ my ago with most of the extant genera evolving - my ago. the average age estimates of most of the tribes imply that they originated within a my time window about my after the origin of harpalinae. approximately my elapsed from the origin of harpalines to the time when all major lineages and tribes are present. tempo of diversification although the sister groups are the same age and descended from a common speciation event, there is a significant difference in species number between harpalinae and brachininae, indicating an increase in net diversification within the harpaline clade. the clade that contains harpalines is also significantly more diverse than the sister group to brachinines + harpalines, the austral psydrines, which contains only extant species. a conservative slowinsk and guyer test, to determine if clades differ in number of species, shows harpalinae with significantly more species than brachininae. the estimated overall net diversification rate with a pure birth model of diversification in brachinines was per my with the beast all-data-combined chronogram to per my with the r8s all data combined chronogram, and in harpalines the estimated overall net diversification rate with a pure birth model was per my with the beast chronogram to per my with the r8s chronogram. the rc test indicated a significant difference in cladogenesis rate between harpalines and brachinines. we located unusually rapid shifts in diversification in the node subtending harpalinae and along at least four basal branches of harpalines. there were also significant increases in diversification within and the tribe lebiini, especially the subtribe calleidina. results of the cr and mccr tests showed the observed value of γ in the harpalinae clade was negative, which rejects the hypothesis that rates of lineage accumulation remained constant over time, however, when compared to the distribution of γ statistics of the simulations, the harpaline γ statistic was not statistically significant. the brachinine γ statistic was not significantly different than zero, and the cr and mccr tests do not reject a constant rate of diversification. the semilogarithmic ltt plots for the harpaline and brachinine clades are shown in figure the accumulation of harpaline lineages did not differ significantly from a constant rate of diversification early in harpaline evolution when compared to simulated yule process trees. however, there was a steep increase in the number of harpaline lineages early after the origin of the subfamily, approximately at my ago, and a shift in rate at about my ago. this was the time period when most harpaline tribes were diversifying according to the beast divergence time estimates. brachinines showed a fairly slow but constant rate of lineage accumulation, within or slightly below the range of simulated trees. for the harpaline clade, a model of two abrupt rate shifts, model 6: yule3rate, was chosen over other constant and variable rate models. a diversification rate decrease was estimated to have happened first at my ago and then again at my ago. statistical analyses of the ltt plot for brachinines indicated that a pure birth model of constant rate of diversification could not be rejected. a model of two abrupt rate shifts, model 6: yule3rate best fit the overall tree with decreases in diversification rates occurring at my ago and my ago. pure birth assumes a pure birth or yule model of diversification, birth death assumes a constant birth death model of diversification, spvar specifies a variable speciation rate and constant extinction rate, exvar specifies a variable extinction rate and constant speciation rate, yule2rate assumes an abrupt change in diversification rate at some breakpoint in the past, yule3rate assumes three different diversification rates with two breakpoints in the past. a breakpoint at my ago b breakpoints at and my ago c breakpoint at my ago d breakpoints at and my ago e breakpoint at my ago f breakpoints at and my ago discussion divergence time estimations of harpalinae we found a mid-cretaceous origin of harpalinae with most tribes arising by the paleocene. crown diversification of tribes appears to have commenced by the late-middle eocene, and nearly all tribes appeared by the late eocene. in a relatively short amount of time, harpalines diversified from just a few lineages into almost all the major lineages extant today. from the divergence date estimates, it appears this rapid diversification was not immediately after the origin of harpalines, but about my later. beast estimates of divergence dates from all molecular data combined were slightly older than the fossil record. r8s estimated divergence dates much older that the fossil record and the beast estimates. differences between the r8s and beast estimates may due to the relaxed clock method used to infer dates and the prior distributions for the divergence time of calibration points in beast versus the minimum and maximum age constraints in r8s. we prefer the dates from the all data combined beast analyses as a more accurate estimate of harpaline evolution based on the partitioned model of molecular evolution and bayesian framework. estimates of divergence dates from molecular data are subject to several sources of error. estimates may be subject to errors in fossil calibration points or topology, inappropriate priors, or invalid clock assumptions, including much more rate heterogeneity than can be accommodated by the programs we used. the error introduced by finite sequence data or model misspecification, such as neglect of multiple hits and across site rate variation can cause spurious results in divergence time estimates. sources of error have been reviewed elsewhere. we did not explicitly explore the effect of differences in tree topology on divergence time estimates. we assumed the ober and maddison 28s+wg maximum likelihood tree was the best phylogenetic hypothesis for harpalines available, but results from the beast analysis of all data combined where topology was also estimated, did not change our conclusions about the timing of harpaline origin or the rapidness of diversification of tribes. our estimates of harpaline divergence dates are a working hypothesis subject to the addition of more taxa, molecular data, improved phylogenetic hypotheses, and fossils yet to be discovered. tempo of diversification harpalines have had an increase in net diversification compared to their sister group, the brachinines. it is difficult to tease apart whether this is due to increased speciation or reduced extinction, because the net diversification rate is simply the difference between speciation and extinction rates. we tested for a change in extinction rate through time in harpalines and brachinines, and a change in extinction rate does not explain the difference in species number between the sister groups because he estimated extinction rate for both is the same. the net diversification rate for harpalines is higher than the rate estimated for explosive diversification of angiosperms and on par with the rapid mammalian radiation. an increase in diversification is evident in the rc test where several harpaline clades have experienced higher than usual cladogenesis. the time of increased cladogenesis occurs approximately - my ago, during a time of incredible faunal and floral turnover in many plant and animal groups. in contrast, there were not shifts in diversification in brachinines or austral psydrines. the rc test is sensitive to temporal depth, phylogenetic scope, and nonindependence of diversification rate shifts. therefore rate shifts identified by rc should be interpreted with caution, especially nested shifts. if a clade is especially diverse, then one or more of its parent clades will also show significant diversification through a "trickle-down" effect. this may be what is being observed along the basal braches of harpalinae, and a conservative interpretation of the rc results would be to attribute the increase in diversification rate to the shallowest significant branch. in this case, there are significant increases in diversification at the base of the harpaline subfamily and within the tribe lebiini, the largest tribe of harpalines with over species. alternatively, a more complex series of radiations may have taken place sequentially. an increase in diversification could have taken place over a longer period of time and shifts in diversification could have occurred over several nodes or branches. evidence of shifts in diversification rates in harpalines is harder to interpret from the cr and mccr tests and the ltt plots. older lineages have a higher risk of being extinct, and thus not sampled, than younger ones. incomplete taxon sampling can give the appearance of nodes near the root of the tree giving rise to more extant descendants than nodes near the tips and are therefore more likely to be included in a small random sample. the mccr test for harpalines did not indicate that nodes were closer to the root than expected under a constant rates model after correcting for incomplete taxon sampling. the taxa included in our phylogenetic tree did not represent a random subsample of the total harpaline species, as assumed by the mccr test with yule process simulations. instead, they most likely represent an oversampling of deeper lineages and an undersampling of closely related lineages, which may have affected the conclusions drawn from this test. furthermore the mccr test suffers from low statistical power and has been shown to have decreased sensitivity in detecting an initial high rate of diversification followed by a shift to a lower rate. the harpaline ltt plot with simulated constant rate trees agreed with the results of the mccr test in that it did not show significantly different rate shifts from trees simulated under a yule process, but the harpaline ltt plot showed a decrease in diversification rate after about my ago, the period after the origin of tribes. although the simulated trees diversified at a constant rate, randomly removing taxa generated an apparent deceleration of lineage accumulation. as with the mccr test, our tree is not a random subsample of harpaline taxa. this bias may affect how the ltt plot is interpreted. the ltt plot and associated confidence interval generated by simulations is just one method by which timing and tempo of diversification can be visualized and evaluated, other lines of evidence suggest a fairly rapid accumulation of harpaline lineages in the paleocene and eocene. the shift in diversification rate from a high rate initially to lower rates later in harpaline evolutionary history are evident in the results of the tests of models of diversification. an increase in extinction rate later in harpaline evolutionary history does not explain the pattern of high diversification early and declining diversification through time. the pattern of lineage accumulation through time rejected a model of changing extinction rate. the best fitting model for harpalines shows a decrease in rate at my ago and again at about my ago. the first rate decrease corresponds well to the end of tribe diversification estimated in beast, and perhaps the second shift, which is much more subtle signals the end of a burst of generic diversification, but more taxa will need to be included to test this hypothesis. estimated dates of divergence, the rc test, evidence from the ltt plot, and models of diversification tests, seem to suggest a late cretaceous through middle eocene period of rapid diversification in harpalines. the paleogene was a time of profound reorganization of biota, perhaps in part due to a period of warm "greenhouse" climate during the paleocene/eocene thermal maximum or continental fission during the late cretaceous. the late cretaceous to middle eocene is the same time period when many other plant and animal taxa are experiencing tremendous changes in evolution from mass extinctions at the k-t boundary to explosive radiations in birds, mammals and some plants and insects. harpaline diversification occurred just after the angiosperm cretaceous radiation - my ago, when flowering plants achieved widespread floristic dominance for the first time. it does not appear that harpalines necessarily co-evolved with angiosperms, but the angiosperm radiation, along with other biotic and abiotic factors, during this period provided harpalines with the possibility of novel ecological interactions. although the vast majority of harpalines are predatory, the evolving complex angiosperm habitats may have afforded harpalines with new niches in which to diversify, especially because some major groups of harpalines, like the species rich tribe lebiini, evolved arboreal lifestyles instead of ground dwelling. recent studies have proposed that another species rich predatory insect group, ants, may have also radiated primarily in the late cretaceous or paleogene.  harpaline species diversity is remarkable compared to its relatives and other carabid clades. their high rates of net diversification have made them an evolutionary success story and an important part of many modern ecosystems. from the dates of divergence estimates, there was a large amount of rapid speciation in a short amount of time in the late cretaceous and early paleogene, but it does not appear that it was a particularly explosive radiation immediately after their evolutionary origin. instead, harpalines underwent a somewhat elevated bursts of diversification about my after their origin. this suggests that harpalines may have had a long "phylogenetic fuse" before tribes and major lineages began to diversify as seen in birds and mammals. relationships of tribes and within the carabid subfamily harpalinae are not strongly supported. the lack of support from the data for deep branches in the harpaline molecular phylogeny could be explained by a rapid radiation or cladogenesis during the evolutionary history of this group. short internal branches created by rapid cladogenesis leave a meager record of diversification that is further potentially obscured by long terminal branches. this explanation was supported in a study attempting to resolve early metazoan evolution and supported by fiala and sokal's simulations. they suggested that short internal branches and long terminal branches reduces the accuracy with which a phylogeny can be estimated, hence the resolution that can be obtained. the pattern seen in the molecular phylogenetic trees of short internal nodes at the base of harpalines with longer terminal branches can be explained in several ways. it may be the case that harpalines did undergo a burst of diversification during the late cretaceous and early paleogene, but the taxon sampling in the molecular phylogeny was not complete enough to distinguish the pattern clearly. incomplete taxon sampling in such a large group makes it difficult to interpret the results of the constant rate tests and ltt plots as a significant burst of diversification in harpalinae evolutionary history. it may be the case that harpalines did diversify rapidly, but it was not as explosive as carabidologists have assumed. on the other hand, the tempo of harpaline evolution may have been fairly constant though time with a high, but steady rate of diversification and lineage accumulation, leading to the large and speciose clade we see today. alternatively, the molecular phylogenies inferred from the three nuclear genes with short basal branches may be due to the information content and phylogenetic utility of the genes or systematic errors in selecting models of evolution for phylogenetic inference. in any case, the inclusion of additional harpaline and brachinine taxa and new molecular data would help shed light on the evolutionary history of the harpalinae and the patterns and tempo of diversification of this group of beetles. authors' contributions kao conceived of the study, carried out some of the analyses, and drafted the manuscript. tnh participated in the design of the study and performed many of the analyses. both authors read and approved the final manuscript. supplementary material additional file 1 table s taxa and genbank numbers used in this study. click here for file acknowledgements duane mckenna, jessica thomas, and three anonymous reviewers generously provided comments that greatly improved the manuscript. the authors thank the college of the holy cross and richard b. fisher summer research fellowship for funding for this project.
here we investigate the tempo of evolution in harpalines using molecular divergence dating techniques and explore the rates of lineage accumulation in harpalines and their sister group. we did not see a significant decrease in diversification rate through time in the mccr test, but a model of diversification with two shift points to lower diversification rates fit the harpaline lineage accumulation through time the best. how this diversity evolved has been difficult to reconstruct, perhaps because harpalines underwent a rapid burst of diversification early in their evolutionary history. in a relatively small window of time, harpalines underwent rapid speciation. harpalinae is a species rich clade of carabid beetles with many unusual morphological forms and ecological interactions. instead of an immediate burst of explosive diversification, harpalines may have had a long "fuse" before major lineages diversified during the early paleogene when other taxa such as mammals, birds, and some flowering plants were also rapidly diversifying. harpalines have a relative high net diversification rate and increased cladogenesis in some regions of the clade. our results indicate harpalines are significantly more diverse and have higher diversification than their sistergroup. according to molecular divergence date estimates, harpalines originated in the mid cretaceous but did not diversify extensively until the late cretaceous or early paleogene about million years after their origin.
 harpalinae is a species rich clade of carabid beetles with many unusual morphological forms and ecological interactions. how this diversity evolved has been difficult to reconstruct, perhaps because harpalines underwent a rapid burst of diversification early in their evolutionary history. here we investigate the tempo of evolution in harpalines using molecular divergence dating techniques and explore the rates of lineage accumulation in harpalines and their sister group.  according to molecular divergence date estimates, harpalines originated in the mid cretaceous but did not diversify extensively until the late cretaceous or early paleogene about million years after their origin. in a relatively small window of time, harpalines underwent rapid speciation. harpalines have a relative high net diversification rate and increased cladogenesis in some regions of the clade. we did not see a significant decrease in diversification rate through time in the mccr test, but a model of diversification with two shift points to lower diversification rates fit the harpaline lineage accumulation through time the best.  our results indicate harpalines are significantly more diverse and have higher diversification than their sistergroup. instead of an immediate burst of explosive diversification, harpalines may have had a long "fuse" before major lineages diversified during the early paleogene when other taxa such as mammals, birds, and some flowering plants were also rapidly diversifying. 