BACKGROUND
the legume family  is the third largest plant family on earth with an estimated  <dig>  species  <cit> .
a study of diversification across angiosperms found fabales  as one of  <dig> angiosperm orders to exhibit higher than expected diversification rates  <cit> .
within fabales, fabaceae exhibits greater diversification than the other relatively species-poor families.
furthermore, within fabaceae, the subfamily papilionoideae contains the bulk of species  compared with mimosoideae and caesalpinioideae.
this uneven distribution of species diversity continues down to lower taxonomic levels in papilionoideae with astragalus containing upwards of  <dig> species  <cit> .
although astragalus does not exhibit higher diversification rates relative to other related genera, the astragalean clade, comprising astragalus and a subset of the galegeae tribe, does  <cit> .
other studies of fabaceae have examined diversification rates and timing of divergence, focusing on generic or higher levels in the family  <cit> .
the collection of research discussed above provides precedent for uneven diversification and shifts in diversification rates within fabaceae, phenomena relatively unexplored at the species level .
studies based on chloroplast coding genes matk and rbcl found elevated nucleotide substitution rates for the tribe psoraleeae relative to other legume clades, suggestive of rapid evolution or diversification  <cit> .
in addition, research on transcontinental clades estimated the age of the divergence between two psoraleeae genera, cullen of australia, and rupertia, endemic to the western united states, at approximately  <dig>  million years ago   <cit> .
this age estimate illustrates the recent diversification of this group, especially as the north american  clade of psoraleeae diversified after this transcontinental split.
evidence suggests that psoraleeae represents a recent, rapid radiation, making this group well suited to a species-level diversification study within fabaceae.
current circumscription of the psoraleeae tribe includes a monophyletic new world group of six genera: hoita , orbexilum , pediomelum , psoralidium , and rupertia  from north america and otholobium  in south america  <cit> .
four old world genera also exist: cullen, bituminaria, psoralea, and otholobium, the latter of which spans the old and new worlds  <cit> .
when considering the nam clade, the skewed distribution of species richness towards pediomelum suggests a shift in diversification rate within nam psoraleeae.
diversification is the net speciation and extinction of a group.
when coupled with evolutionary and ecological processes, diversification is central to the creation of species richness.
differential species richness across taxonomic and geographic levels has been widely studied, with many hypotheses proposed to explain this variation .
most hypotheses concerning species richness focus on the impact of climatic factors on the diversification of species across geographic regions, endeavoring to find and explain global patterns of biodiversity.
hypotheses concerning the cause of diversification rate shifts, or net changes in speciation or extinction rate that bear on relative species diversity, have also been suggested including geographic range fragmentation, ecology, competition, key innovations, and climate change  <cit> .
climatic variables such as water and temperature account for much of the variation in species richness seen globally  <cit> .
indeed, climatic variables have been implicated in studies of species richness and diversification of a wide range of organisms, including angiosperms  <cit> , amphibians  <cit> , reptiles  <cit> , vertebrates  <cit> , insects  <cit> , and birds  <cit> , among others  <cit> .
hypotheses of climate change are particularly relevant to nam psoraleeae as climatic variables may affect habitat and range variation past and present.
the five nam genera span a variety of ecological habitats: hoita is riparian, living strictly in moist soils along streams and springs in california; rupertia resides on moderately moist slopes in oak-pine communities of california; orbexilum resides in the semi-moist soils of forest communities of the southeastern united states; psoralidium lives in the semi-dry soils of midwestern and intermountain grasslands and deserts; and pediomelum has adapted to well-drained soils in dry, arid ecosystems of the deserts of the intermountain west and the grasslands and rock outcrops of the midwestern united states  <cit> .
many habitats of endemic pediomelum species are characterized by arid, hot environments such as the cedar glades of tennessee, the sandy deserts of the southwestern united states, or the limestone outcrops of the edwards plateau in texas.
observation of nam psoraleeae in light of evolutionary relationships suggests a shift in climate regimes with pediomelum, the largest genus, and psoralidium shifting to more xeric habitats, perhaps leading to increased diversification in these genera relative to other nam psoraleeae.
this leads to the postulate that adaptation to xeric habitats may have acted as a key innovation, a trait whose origin has spurred diversification within a group  <cit> .
the young age estimate of the tribe suggests that nam psoraleeae may have diversified in part in the quaternary era, a period replete with climate oscillations brought on by glacial cycles.
these shifts greatly influenced the range and distribution of species across north america, as ranges were contracted and expanded repeatedly.
the onset of climate change required plant species to move, adapt, or go extinct.
upon glacial recession, new environments were open to colonization or expansion, often requiring adaptation of existing flora  <cit> .
these cycles could influence speciation either through genetic differentiation spurred by range fragmentation into glacial refugia or through admixture or introgression as previously isolated entities mixed and created new and differing gene pools which contributed to speciation  <cit> .
these events, coupled with selection pressures caused by new environments or competitors, may have spurred genetic and/or morphological change, perhaps influencing the tempo of diversification.
the idea of quaternary climate shifts impacting genetic divergence and diversification is not new  <cit> .
while most studies involving the influence of this time period investigate intraspecific and population-level genetic divergence , studies concerning the impact of this era on speciation also exist but are more rare .
if nam psoraleeae diversified during this time, the tempo of diversification of this group may have been influenced by these climate shifts.
in this study, we present nam psoraleeae as a model study system for investigating the impact of climate change on diversification within a recent, rapid radiation.
we test the hypothesis of a diversification rate shift within nam psoraleeae against the null hypothesis of no shift in diversification rate within the group.
we then explore diversification in light of various climate-related hypotheses.
in particular, we test the hypothesis that a shift to xeric habitats acted as a key innovation, thus spurring a diversification rate increase in pediomelum and psoralidium as opposed to the rest of nam psoraleeae against the null hypothesis of no net increase in diversification rate of these genera.
we also investigate the impact of quaternary climate oscillations on diversification of nam psoraleeae.
to do so, we estimate divergence times of nam psoraleeae to determine if the group was diversifying during this period.
we then correlate diversification rate shifts to geologic time.
the young age of this group provides a glimpse into the emerging evolution of nam psoraleeae and the impact of changing environments and climes on its radiation.
in addition, the knowledge gained herein may provide insight into other species radiations that have endured similar conditions, and serves as a model for the study of diversification rates and ecological fluctuations in the robust phylogenetic context of a recent, rapid radiation.
RESULTS
divergence dates
a molecular clock was rejected for each data set , justifying the use of a relaxed molecular clock for divergence dating.
the mixed model combined analysis from beast yielded a maximum clade credibility tree very similar to previously published phylogenies  <cit> , suggesting a monophyletic nam psoraleeae .
individual gene analyses differed somewhat in their divergence date estimates of key nodes across psoraleeae, illustrating the variation in rates of molecular evolution and phylogenetic information content inherent in each gene region .
for the most part, date estimates determined using nuclear markers are similar to or younger than dates based on chloroplast markers.
two exceptions to this pattern exist: the node representing the most recent common ancestor  of rupertia has the minimum date  and the maximum date  estimated by nuclear markers, while the dates estimated for the mrca of glycine and psoraleeae by nuclear genes exceed those estimated by chloroplast markers.
mrca is most recent common ancestor.
sd is standard deviation.
95% upper and lower hpd is high posterior density credibility interval.
ave is the average of the dates across individual gene estimates with accompanying standard deviation below.
total is the date based on the combined dataset using mixed models.
considering the gene regions as a whole, the combined mixed model analysis estimated the divergence date of calibrated nodes closer to the date expected than an analysis that averaged date estimates across all gene regions.
for example, the calibrated node of the mrca of glycine and phaseolus should have had a mean of  <dig>  mya; the combined analysis estimated this divergence date at  <dig>  mya while the average across the individual genes for this node was  <dig>  mya .
with the exception of glycine and the mrca of glycine and phaseolus, the mean divergence date summarized across individual genes is older than the mean divergence date estimated under the combined mixed model analysis.
the minimum and maximum divergence date estimates across the individual genes could be taken as a range of uncertainty surrounding the divergence date of respective nodes, much the way we utilize the credibility intervals from bayesian analyses.
the ranges suggested across individual gene regions differ somewhat from the 95% high posterior density  credibility interval determined in the combined mixed model analysis.
with the exception of the mrcas of glycine/psoraleeae, orbexilum, cullen, and glycine, the range obtained from the minimum and maximum estimates across individual gene regions is more conservative  than the bayesian credibility interval from the combined analysis .
while divergence estimates from the combined mixed model analysis differed from the average across individual gene estimates, the difference is slight.
furthermore, the divergence ranges suggested by utilizing the minimum and maximum estimates across individual genes or the range suggested by the 95% confidence interval  inferred using the mean across individual genes both exhibit considerable overlap with the range suggested by the 95% hpd credibility interval.
we therefore confine our discussion to the dates and ranges observed under the combined mixed model framework.
as estimated in the combined analysis, psoraleeae began diversifying in north america around  <dig>  mya with the split of the genus hoita from the remaining taxa .
orbexilum split from the ancestor of pediomelum, psoralidium, and rupertia shortly thereafter, with the genus beginning to diversify around  <dig>  mya .
around  <dig>  mya, the pediomelum lineage split from the lineage leading to psoralidium and rupertia, two genera which would later split around  <dig>  mya.
rupertia began diversifying around  <dig>  mya while psoralidium speciated around  <dig>  mya.
pediomelum, the largest genus, began diversifying around  <dig>  mya.
differential diversification rates
the results of the seven whole-tree statistics computed in symmetree using the  <dig> maximum parsimony tree topologies suggest a significant shift in diversification rates within nam psoraleeae with strong support .
the delta statistics, those that help to localize the diversification rate shift to a certain area of the phylogeny, did not find statistical support for a single shift point of diversification rate.
however, the test does support that a shift took place, either gradually over the whole phylogeny or in multiple places over time.
the birth-death likelihood  analysis based on the maximum clade credibility tree from beast determined the maximum likelihood estimate of the speciation rate under the purebirth model to be  <dig>  speciation events per million years for the combined data.
bdl chose purebirth as the best rate-constant model .
it was therefore used as the null to calculate Δaicrc for the rate variable models.
the critical value for Δaicrc, as compared with the best rate-constant model, to maintain a significance level of α =  <dig>  is  <dig> .
bdl found yule3rate as the best rate-variable model .
we therefore reject the null hypothesis of rate-constancy because the observed Δaicrc >  <dig> .
according to the scenario suggested by the yule3rate model, nam psoraleeae began diversifying with a net diversification rate of  <dig>  speciation events per million years.
a shift in net diversification rate took place  <dig>  mya in which the rate shifts dramatically to  <dig>  speciation events per million years.
the net diversification rate shifted again around  <dig>  mya, decreasing to  <dig>  speciation events per million years .
a lineage through time plot illustrates the definitive 3-rate model as outlined in the bdl analysis .
combined evidence from the bdl and divergence date analyses places the diversification rate shift at the origin of the main pediomelum radiation.
the maximum clade credibility chronogram was calibrated such that the divergence of hoita and pediomelum was scaled to  <dig>  mya.
r = net diversification rate ; a = extinction fraction; st = time of rate shift ; k = carrying capacity parameter; x = rate change parameter; ln = log-likelihood; aic = akaike information criterion; Δaic = change in aic relative to purebirth.
testing of macroevolutionary hypotheses
the hypothesis of xeric habitats acting as a key innovation in the diversification of nam psoraleeae based on the  <dig> posterior distribution trees from the beast analysis was rejected  using ree's key innovation test.
the hypothesis that quaternary climate shifts influenced diversification of nam psoraleeae is supported by the acceptance of the yule3rate model having a diversification rate shift, st, between  <dig>  and  <dig>  mya, the time span equating to glaciation-induced climate oscillations during the quaternary  <cit> .
this translates to a  <dig> -fold increase in net diversification rate around  <dig> mya .
discussion
recent, rapid radiations
the legume family is somewhat remarkable in that age estimates of this incredibly diverse group are relatively young, with transcontinental crown clade ages ranging from  <dig> to  <dig> mya and intracontinental crown clade ages ranging from  <dig> to  <dig> mya, providing a mostly quaternary to neogene age distribution  <cit> .
while astragalus is the legume poster child for rapid radiations, other recent, rapid radiations exist within fabaceae  <cit> .
such radiations are often characterized by short branch lengths, little sequence variation, and diverse morphology and ecological habitats  <cit> .
nam psoraleeae, especially pediomelum, exhibits all of these characteristics  <cit> .
this, coupled with evidence of increased rates of molecular evolution and recent age estimates of transcontinental clades, suggests that tribe psoraleeae is a recent, rapid radiation  <cit> .
divergence dates
our divergence date estimates place the origin of nam psoraleeae between  <dig>  and  <dig>  mya with a mean estimate of  <dig>  mya, concomitant with many other legume genera  <cit> .
of particular interest in terms of divergence dates is the largest genus, pediomelum.
grimes' circumscription segregates pediomelum into three subgenera based on branching pattern and dehiscence of inflorescences: subgenus leucocraspedon includes p. rhombifolium and p. palmeri with remaining taxa split between subgenera disarticulatum and pediomelum  <cit> .
the diversification of this genus began around  <dig> mya  with the split of subgenus leucocraspedon from the remainder of pediomelum.
the two species in subgenus leucocraspedon differ greatly from the rest of the genus both morphologically and genetically  <cit> , this evidenced in the divergence of the two main pediomelum lineages.
the remainder of pediomelum began diversifying around  <dig> mya, constituting the main body of species richness of nam psoraleeae.
to place this in perspective with other rapid radiations, the hawaiian silverswords, the model for adaptive radiations, diverged around  <dig> mya  <cit> .
within fabaceae, astragalus was estimated to have diverged from oxytropis around  <dig> to  <dig> mya with neo-astragalus radiating approximately  <dig>  mya  <cit> , while the andean radiation of lupinus is dated around  <dig>  mya  <cit> .
these dates place the diversification of pediomelum on level ground with other well-known examples of recent radiations.
diversification rates
a rapid, family-wide diversification is strongly supported in the legume family  <cit> .
psoraleeae is no exception.
estimates of diversification rates vary depending on the model of diversification used in bdl on the maximum clade credibility tree from beast .
however, all our estimates of nam psoraleeae diversification are as fast or faster than rates determined for other notable rapid radiations: under the purebirth model, a diversification rate of  <dig>  for nam psoraleeae is similar to other rapid radiations including the hawaiian silverswords  <cit> , mosses  <cit> , squamate reptiles  <cit> , and salamanders  <cit> , all of which range between  <dig>  and  <dig>  speciation events per million years.
skewed species richness to pediomelum suggests a shift in net diversification rates within nam psoraleeae.
our analyses, including both topological and temporal methods, confirmed this hypothesis.
whole-tree statistics based on topology suggested a shift in diversification rate mid-depth in the phylogeny of nam psoraleeae.
while a shift was strongly supported, the topological method was unable to decipher between the presence of multiple shifts in diversification rates along the phylogeny or a gradual change in diversification rate throughout the phylogeny as a whole.
in addition, the topological method could not determine whether the shift or shifts were net increases or decreases in diversification rate.
this lack of information is due to limitations of the topological method stemming from the fact that it relies solely on the order of branching events as opposed to the timing of branching events.
the greater power of bdl, a temporal method, found two shifts in net diversification rate along the phylogeny.
according to diversification rates suggested by the yule3rate model and bdl, nam psoraleeae diversification rates shifted dramatically with a net increase from  <dig>  to  <dig>  speciation events per million years approximately  <dig> mya.
the increased rate is one of the fastest diversification rates published to date, rivaling those of ice plants in south africa  <cit>  and the andean radiation of lupinus  <cit> .
the increase in net diversification rate did indeed occur at mid-depth in the phylogeny and represents a large increase in the diversification of the tribe within north america.
the concordance of dates between the divergence date and bdl analyses confirms that the diversification rate shift took place as the main radiation of pediomelum , began diversifying around  <dig> mya .
in addition to a net increase, analyses also provide evidence for a net decrease in diversification rate from  <dig>  to  <dig>  speciation events per million years around  <dig>  years ago.
this rate is far lower than the previous rates, suggesting a decrease in the amount of speciation or an increase in the amount of extinction or both.
the low speciation rate seen at this recent time point could represent the speciation/population boundary.
diversification from this point forward to the present could be seen as existing on a population-level scale, as genetic differentiation between populations not sufficiently diverged to be distinguished as species.
furthermore, this decrease in diversification rate may be indirectly influenced by recent events, including those of anthropogenic origin.
indeed, we know for certain that human influence has been the cause of extinction within nam psoraleeae: orbexilum stipulatum, once endemic to rock island, kentucky, at falls of the ohio river, is confirmed extinct after the island was destroyed for construction of a dam; orbexilum macrophyllum was collected from a single site in north carolina and is presumed extinct  <cit> .
still other species are severely threatened, such as pediomelum pentaphyllum and p. humile, due to the elimination of extant populations through urban development, leaving only one or two known populations alive.
estimates of more recent diversification shifts are especially dependent on the number of extant species.
if the number of species is depleted through recent events such as human influence or through the lack of species recognition, it will more likely be reflected in diversification rate shifts taking place at more recent times, this due to the fact that we are estimating the diversification rate based on the number of species, instead of the number of clades.
diversification shifts – causes?
uneven diversification is evident in both taxonomic and geographical realms, exemplified in the ideas of species richness and the latitudinal diversity gradient.
a great number of hypotheses have been proposed as explanations for uneven diversification with many centering on key innovations or climate.
while the pattern of the latitudinal diversity gradient, meaning the existence of greater species numbers centered around lower latitudes, is well documented, the mechanism behind the phenomenon is unclear, leading to a plethora of hypotheses  <cit> .
it is possible and even probable that no single hypothesis can explain such patterns  <cit> .
furthermore, these and other patterns could be the results of purely stochastic processes.
however, evidence has shown that climatic factors are especially influential on the evolutionary history of many organismal groups .
as nam psoraleeae is within temporal and geographical range of quaternary climate oscillations and exhibits a wide variety of climate-based habitats, we chose to address hypotheses involving climate change as possible effectors of diversification rate shifts.
habitat shifts
past climate change has impacted current climate regimes, creating a spectrum of climate-based habitats across the earth.
movement of ancestral species from one habitat to another may also have spurred diversification, relegating climate-based habitat as a key innovation.
climate regime shifts have been found to have impacted diversification across various organismal groups, including plants  <cit>  and lizards  <cit> .
hoita, orbexilum, and rupertia all reside in more mesic environments characterized by warm, moist, humid climates whereas psoralidium and pediomelum reside in more xeric habitats characterized by hot, dry environments such as the deserts of the intermountain and southwest united states .
in testing habitat shifts from mesic to xeric habitats, our analyses suggest that past climate change affected diversification rates in nam psoraleeae more than shifts in present climate regimes as defined by current habitats.
the shift to xeric habitats was not considered a key innovation in terms of increasing diversification rates within the group.
instead, current habitats could be seen as relicts from the effects of past climate change on its influence on the diversification of north american psoraleeae.
quaternary climate oscillations
the recent origin of the main radiation of pediomelum around  <dig> mya places its diversification in the middle of the pleistocene epoch of the quaternary era, a time period replete with dramatic glaciation-induced climate oscillations.
these glacial fluxes include milankovitch cycles which consist of repeated formation and recession of glacial ice sheets having periodicities of  <dig>  to  <dig>  years, resulting in wide-ranging climate fluxes  <cit> .
it has been well documented that much of the earth's biodiversity and many species' evolutionary trajectories were drastically influenced by the climate fluxes brought on by quaternary glaciations .
in north america, the sheer presence of ice sheets shifted climates and strongly influenced the range and distribution of species  <cit> .
ice formation may have influenced speciation via genetic differentiation through range fragmentations such as glacial refugia.
upon glacial recession, new environments were formed and subsequently colonized through species expansion and adaptation  <cit> .
given the geographical and temporal distribution of nam psoraleeae, these cycles may have influenced its evolutionary history.
indeed, the shift point for the increase in diversification along the phylogeny at  <dig> mya and the lineage through time plot from our analyses support the hypothesis that quaternary climate fluxes, occurring between  <dig>  and  <dig>  mya, influenced the diversification of nam psoraleeae.
geography
the effects of these oscillations are not homogenous, however.
the quaternary ice ages and climate oscillations affected geographic regions across north america differently.
the southeast provided several refugia in which species could survive  <cit>  while the montane regions of the west were characterized by more varied environments such as ice, tundra, pluvial lakes, and deserts  <cit> .
in addition, elevation played an important role in species distributions as elevational ranges of plants were  <dig> to 1200+ meters below modern ranges, an effect felt more readily in the topographically diverse regions of the west  <cit> .
upon glacial recession, the intermountain west offered a greater variety of habitats for colonization and adaptation as a result of topography.
this is in line with cracraft's lithospheric complexity hypothesis, stating that geographic separation spurred by topographical changes promotes speciation  <cit> .
in addition to orogeny, glacial cycles contributed to topographical change, and thus are linked with the lithospheric hypothesis.
geographical trends are seen in the diversity of psoraleeae.
nam psoraleeae ranges throughout the united states: orbexilum is mostly confined to the southeast united states, with one species reaching into mexico, pediomelum is mostly found in central and intermountain states, hoita is endemic to california, rupertia ranges along the west coast, and psoralidium is spread over the great plains with one species endemic to utah.
of interest is the seeming correlation of lineage age with southerly distribution.
orbexilum is one of the oldest nam psoraleeae genera, diversifying earlier than most other lineages.
this could be demonstrative of long-established ranges of orbexilum through pleistocene refugia in the southeast united states.
of import also is the age of pediomelum as a whole, beginning its diversification around  <dig> mya with the break off of subgenus leucocraspedon.
the two leucocraspedon species exhibit the most southern distribution of all pediomelum species with p. rhombifolium being distributed throughout texas and into mexico while p. palmeri is almost exclusively found in mesoamerica.
this southern distribution perhaps allowed this lineage to persist through the quaternary climate oscillations as the effects of glacial climate change were moderated toward the equator.
the main radiation of pediomelum exhibits the greatest number of species with the lowest intraspecific variation, findings in accord with the lithospheric complexity hypothesis, geographic separation, ecological diversification due to habitat diversity, and range fragmentation brought on by quaternary climate shifts.
the pattern of less diversity in the southeast  and greater diversity in the southwest  is also found in cicindela beetles  <cit>  and agelenopsis spiders  <cit> , among others.
CONCLUSIONS
nam psoraleeae has proved useful for studying evolutionary divergence and the impact of climate change on diversification within a recent, rapid radiation.
this work effectively demonstrates differential rates of diversification at the species and generic levels, trends manifested across taxonomic levels within fabaceae.
an estimation of divergence dates within psoraleeae found that the main radiation of pediomelum, comprising the bulk of species richness in nam psoraleeae, originated  <dig> to  <dig> mya during the quaternary era, concurrent with the milankovitch climate oscillations brought on by glacial cycles.
the hypothesis of a diversification rate shift in nam psoraleeae was confirmed by both topological and temporal methods.
we demonstrated a large increase in diversification rate near the beginning of the pleistocene epoch.
our findings suggest that the diversification rate increase in nam psoraleeae is not due to current climate regimes as represented by habitat, but rather imply that quaternary climate oscillations impacted the diversification of the group.
we show nonrandom spatial diversification as well, with higher species diversity in the american southwest.
nam psoraleeae diversification is a good example of how earthly dynamics including global climate change and topography work together to shape biodiversity.
in a time when the issue of global climate change pervades society and science, this study aptly illustrates the control that global climate change can have over the evolutionary trajectories and fates of lineages and species.

 the legume family is the third largest plant family on earth with an estimated species. a study of diversification across angiosperms found fabales as one of angiosperm orders to exhibit higher than expected diversification rates. within fabales, fabaceae exhibits greater diversification than the other relatively species-poor families. furthermore, within fabaceae, the subfamily papilionoideae contains the bulk of species compared with mimosoideae and caesalpinioideae. this uneven distribution of species diversity continues down to lower taxonomic levels in papilionoideae with astragalus containing upwards of species. although astragalus does not exhibit higher diversification rates relative to other related genera, the astragalean clade, comprising astragalus and a subset of the galegeae tribe, does. other studies of fabaceae have examined diversification rates and timing of divergence, focusing on generic or higher levels in the family. the collection of research discussed above provides precedent for uneven diversification and shifts in diversification rates within fabaceae, phenomena relatively unexplored at the species level. studies based on chloroplast coding genes matk and rbcl found elevated nucleotide substitution rates for the tribe psoraleeae relative to other legume clades, suggestive of rapid evolution or diversification. in addition, research on transcontinental clades estimated the age of the divergence between two psoraleeae genera, cullen of australia, and rupertia, endemic to the western united states, at approximately million years ago. this age estimate illustrates the recent diversification of this group, especially as the north american clade of psoraleeae diversified after this transcontinental split. evidence suggests that psoraleeae represents a recent, rapid radiation, making this group well suited to a species-level diversification study within fabaceae. current circumscription of the psoraleeae tribe includes a monophyletic new world group of six genera: hoita, orbexilum, pediomelum, psoralidium, and rupertia from north america and otholobium in south america. four old world genera also exist: cullen, bituminaria, psoralea, and otholobium, the latter of which spans the old and new worlds. when considering the nam clade, the skewed distribution of species richness towards pediomelum suggests a shift in diversification rate within nam psoraleeae. diversification is the net speciation and extinction of a group. when coupled with evolutionary and ecological processes, diversification is central to the creation of species richness. differential species richness across taxonomic and geographic levels has been widely studied, with many hypotheses proposed to explain this variation. most hypotheses concerning species richness focus on the impact of climatic factors on the diversification of species across geographic regions, endeavoring to find and explain global patterns of biodiversity. hypotheses concerning the cause of diversification rate shifts, or net changes in speciation or extinction rate that bear on relative species diversity, have also been suggested including geographic range fragmentation, ecology, competition, key innovations, and climate change. climatic variables such as water and temperature account for much of the variation in species richness seen globally. indeed, climatic variables have been implicated in studies of species richness and diversification of a wide range of organisms, including angiosperms, amphibians, reptiles, vertebrates, insects, and birds, among others. hypotheses of climate change are particularly relevant to nam psoraleeae as climatic variables may affect habitat and range variation past and present. the five nam genera span a variety of ecological habitats: hoita is riparian, living strictly in moist soils along streams and springs in california; rupertia resides on moderately moist slopes in oak-pine communities of california; orbexilum resides in the semi-moist soils of forest communities of the southeastern united states; psoralidium lives in the semi-dry soils of midwestern and intermountain grasslands and deserts; and pediomelum has adapted to well-drained soils in dry, arid ecosystems of the deserts of the intermountain west and the grasslands and rock outcrops of the midwestern united states. many habitats of endemic pediomelum species are characterized by arid, hot environments such as the cedar glades of tennessee, the sandy deserts of the southwestern united states, or the limestone outcrops of the edwards plateau in texas. observation of nam psoraleeae in light of evolutionary relationships suggests a shift in climate regimes with pediomelum, the largest genus, and psoralidium shifting to more xeric habitats, perhaps leading to increased diversification in these genera relative to other nam psoraleeae. this leads to the postulate that adaptation to xeric habitats may have acted as a key innovation, a trait whose origin has spurred diversification within a group. the young age estimate of the tribe suggests that nam psoraleeae may have diversified in part in the quaternary era, a period replete with climate oscillations brought on by glacial cycles. these shifts greatly influenced the range and distribution of species across north america, as ranges were contracted and expanded repeatedly. the onset of climate change required plant species to move, adapt, or go extinct. upon glacial recession, new environments were open to colonization or expansion, often requiring adaptation of existing flora. these cycles could influence speciation either through genetic differentiation spurred by range fragmentation into glacial refugia or through admixture or introgression as previously isolated entities mixed and created new and differing gene pools which contributed to speciation. these events, coupled with selection pressures caused by new environments or competitors, may have spurred genetic and/or morphological change, perhaps influencing the tempo of diversification. the idea of quaternary climate shifts impacting genetic divergence and diversification is not new. while most studies involving the influence of this time period investigate intraspecific and population-level genetic divergence, studies concerning the impact of this era on speciation also exist but are more rare. if nam psoraleeae diversified during this time, the tempo of diversification of this group may have been influenced by these climate shifts. in this study, we present nam psoraleeae as a model study system for investigating the impact of climate change on diversification within a recent, rapid radiation. we test the hypothesis of a diversification rate shift within nam psoraleeae against the null hypothesis of no shift in diversification rate within the group. we then explore diversification in light of various climate-related hypotheses. in particular, we test the hypothesis that a shift to xeric habitats acted as a key innovation, thus spurring a diversification rate increase in pediomelum and psoralidium as opposed to the rest of nam psoraleeae against the null hypothesis of no net increase in diversification rate of these genera. we also investigate the impact of quaternary climate oscillations on diversification of nam psoraleeae. to do so, we estimate divergence times of nam psoraleeae to determine if the group was diversifying during this period. we then correlate diversification rate shifts to geologic time. the young age of this group provides a glimpse into the emerging evolution of nam psoraleeae and the impact of changing environments and climes on its radiation. in addition, the knowledge gained herein may provide insight into other species radiations that have endured similar conditions, and serves as a model for the study of diversification rates and ecological fluctuations in the robust phylogenetic context of a recent, rapid radiation.  divergence dates a molecular clock was rejected for each data set, justifying the use of a relaxed molecular clock for divergence dating. the mixed model combined analysis from beast yielded a maximum clade credibility tree very similar to previously published phylogenies, suggesting a monophyletic nam psoraleeae. individual gene analyses differed somewhat in their divergence date estimates of key nodes across psoraleeae, illustrating the variation in rates of molecular evolution and phylogenetic information content inherent in each gene region. for the most part, date estimates determined using nuclear markers are similar to or younger than dates based on chloroplast markers. two exceptions to this pattern exist: the node representing the most recent common ancestor of rupertia has the minimum date and the maximum date estimated by nuclear markers, while the dates estimated for the mrca of glycine and psoraleeae by nuclear genes exceed those estimated by chloroplast markers. mrca is most recent common ancestor. sd is standard deviation. 95% upper and lower hpd is high posterior density credibility interval. ave is the average of the dates across individual gene estimates with accompanying standard deviation below. total is the date based on the combined dataset using mixed models. considering the gene regions as a whole, the combined mixed model analysis estimated the divergence date of calibrated nodes closer to the date expected than an analysis that averaged date estimates across all gene regions. for example, the calibrated node of the mrca of glycine and phaseolus should have had a mean of mya; the combined analysis estimated this divergence date at mya while the average across the individual genes for this node was mya. with the exception of glycine and the mrca of glycine and phaseolus, the mean divergence date summarized across individual genes is older than the mean divergence date estimated under the combined mixed model analysis. the minimum and maximum divergence date estimates across the individual genes could be taken as a range of uncertainty surrounding the divergence date of respective nodes, much the way we utilize the credibility intervals from bayesian analyses. the ranges suggested across individual gene regions differ somewhat from the 95% high posterior density credibility interval determined in the combined mixed model analysis. with the exception of the mrcas of glycine/psoraleeae, orbexilum, cullen, and glycine, the range obtained from the minimum and maximum estimates across individual gene regions is more conservative than the bayesian credibility interval from the combined analysis. while divergence estimates from the combined mixed model analysis differed from the average across individual gene estimates, the difference is slight. furthermore, the divergence ranges suggested by utilizing the minimum and maximum estimates across individual genes or the range suggested by the 95% confidence interval inferred using the mean across individual genes both exhibit considerable overlap with the range suggested by the 95% hpd credibility interval. we therefore confine our discussion to the dates and ranges observed under the combined mixed model framework. as estimated in the combined analysis, psoraleeae began diversifying in north america around mya with the split of the genus hoita from the remaining taxa. orbexilum split from the ancestor of pediomelum, psoralidium, and rupertia shortly thereafter, with the genus beginning to diversify around mya. around mya, the pediomelum lineage split from the lineage leading to psoralidium and rupertia, two genera which would later split around mya. rupertia began diversifying around mya while psoralidium speciated around mya. pediomelum, the largest genus, began diversifying around mya. differential diversification rates the results of the seven whole-tree statistics computed in symmetree using the maximum parsimony tree topologies suggest a significant shift in diversification rates within nam psoraleeae with strong support. the delta statistics, those that help to localize the diversification rate shift to a certain area of the phylogeny, did not find statistical support for a single shift point of diversification rate. however, the test does support that a shift took place, either gradually over the whole phylogeny or in multiple places over time. the birth-death likelihood analysis based on the maximum clade credibility tree from beast determined the maximum likelihood estimate of the speciation rate under the purebirth model to be speciation events per million years for the combined data. bdl chose purebirth as the best rate-constant model. it was therefore used as the null to calculate Δaicrc for the rate variable models. the critical value for Δaicrc, as compared with the best rate-constant model, to maintain a significance level of α = is. bdl found yule3rate as the best rate-variable model. we therefore reject the null hypothesis of rate-constancy because the observed Δaicrc >. according to the scenario suggested by the yule3rate model, nam psoraleeae began diversifying with a net diversification rate of speciation events per million years. a shift in net diversification rate took place mya in which the rate shifts dramatically to speciation events per million years. the net diversification rate shifted again around mya, decreasing to speciation events per million years. a lineage through time plot illustrates the definitive 3-rate model as outlined in the bdl analysis. combined evidence from the bdl and divergence date analyses places the diversification rate shift at the origin of the main pediomelum radiation. the maximum clade credibility chronogram was calibrated such that the divergence of hoita and pediomelum was scaled to mya. r = net diversification rate; a = extinction fraction; st = time of rate shift; k = carrying capacity parameter; x = rate change parameter; ln = log-likelihood; aic = akaike information criterion; Δaic = change in aic relative to purebirth. testing of macroevolutionary hypotheses the hypothesis of xeric habitats acting as a key innovation in the diversification of nam psoraleeae based on the posterior distribution trees from the beast analysis was rejected using ree's key innovation test. the hypothesis that quaternary climate shifts influenced diversification of nam psoraleeae is supported by the acceptance of the yule3rate model having a diversification rate shift, st, between and mya, the time span equating to glaciation-induced climate oscillations during the quaternary. this translates to a -fold increase in net diversification rate around mya. discussion recent, rapid radiations the legume family is somewhat remarkable in that age estimates of this incredibly diverse group are relatively young, with transcontinental crown clade ages ranging from to mya and intracontinental crown clade ages ranging from to mya, providing a mostly quaternary to neogene age distribution. while astragalus is the legume poster child for rapid radiations, other recent, rapid radiations exist within fabaceae. such radiations are often characterized by short branch lengths, little sequence variation, and diverse morphology and ecological habitats. nam psoraleeae, especially pediomelum, exhibits all of these characteristics. this, coupled with evidence of increased rates of molecular evolution and recent age estimates of transcontinental clades, suggests that tribe psoraleeae is a recent, rapid radiation. divergence dates our divergence date estimates place the origin of nam psoraleeae between and mya with a mean estimate of mya, concomitant with many other legume genera. of particular interest in terms of divergence dates is the largest genus, pediomelum. grimes' circumscription segregates pediomelum into three subgenera based on branching pattern and dehiscence of inflorescences: subgenus leucocraspedon includes p. rhombifolium and p. palmeri with remaining taxa split between subgenera disarticulatum and pediomelum. the diversification of this genus began around mya with the split of subgenus leucocraspedon from the remainder of pediomelum. the two species in subgenus leucocraspedon differ greatly from the rest of the genus both morphologically and genetically, this evidenced in the divergence of the two main pediomelum lineages. the remainder of pediomelum began diversifying around mya, constituting the main body of species richness of nam psoraleeae. to place this in perspective with other rapid radiations, the hawaiian silverswords, the model for adaptive radiations, diverged around mya. within fabaceae, astragalus was estimated to have diverged from oxytropis around to mya with neo-astragalus radiating approximately mya, while the andean radiation of lupinus is dated around mya. these dates place the diversification of pediomelum on level ground with other well-known examples of recent radiations. diversification rates a rapid, family-wide diversification is strongly supported in the legume family. psoraleeae is no exception. estimates of diversification rates vary depending on the model of diversification used in bdl on the maximum clade credibility tree from beast. however, all our estimates of nam psoraleeae diversification are as fast or faster than rates determined for other notable rapid radiations: under the purebirth model, a diversification rate of for nam psoraleeae is similar to other rapid radiations including the hawaiian silverswords, mosses, squamate reptiles, and salamanders, all of which range between and speciation events per million years. skewed species richness to pediomelum suggests a shift in net diversification rates within nam psoraleeae. our analyses, including both topological and temporal methods, confirmed this hypothesis. whole-tree statistics based on topology suggested a shift in diversification rate mid-depth in the phylogeny of nam psoraleeae. while a shift was strongly supported, the topological method was unable to decipher between the presence of multiple shifts in diversification rates along the phylogeny or a gradual change in diversification rate throughout the phylogeny as a whole. in addition, the topological method could not determine whether the shift or shifts were net increases or decreases in diversification rate. this lack of information is due to limitations of the topological method stemming from the fact that it relies solely on the order of branching events as opposed to the timing of branching events. the greater power of bdl, a temporal method, found two shifts in net diversification rate along the phylogeny. according to diversification rates suggested by the yule3rate model and bdl, nam psoraleeae diversification rates shifted dramatically with a net increase from to speciation events per million years approximately mya. the increased rate is one of the fastest diversification rates published to date, rivaling those of ice plants in south africa and the andean radiation of lupinus. the increase in net diversification rate did indeed occur at mid-depth in the phylogeny and represents a large increase in the diversification of the tribe within north america. the concordance of dates between the divergence date and bdl analyses confirms that the diversification rate shift took place as the main radiation of pediomelum, began diversifying around mya. in addition to a net increase, analyses also provide evidence for a net decrease in diversification rate from to speciation events per million years around years ago. this rate is far lower than the previous rates, suggesting a decrease in the amount of speciation or an increase in the amount of extinction or both. the low speciation rate seen at this recent time point could represent the speciation/population boundary. diversification from this point forward to the present could be seen as existing on a population-level scale, as genetic differentiation between populations not sufficiently diverged to be distinguished as species. furthermore, this decrease in diversification rate may be indirectly influenced by recent events, including those of anthropogenic origin. indeed, we know for certain that human influence has been the cause of extinction within nam psoraleeae: orbexilum stipulatum, once endemic to rock island, kentucky, at falls of the ohio river, is confirmed extinct after the island was destroyed for construction of a dam; orbexilum macrophyllum was collected from a single site in north carolina and is presumed extinct. still other species are severely threatened, such as pediomelum pentaphyllum and p. humile, due to the elimination of extant populations through urban development, leaving only one or two known populations alive. estimates of more recent diversification shifts are especially dependent on the number of extant species. if the number of species is depleted through recent events such as human influence or through the lack of species recognition, it will more likely be reflected in diversification rate shifts taking place at more recent times, this due to the fact that we are estimating the diversification rate based on the number of species, instead of the number of clades. diversification shifts – causes? uneven diversification is evident in both taxonomic and geographical realms, exemplified in the ideas of species richness and the latitudinal diversity gradient. a great number of hypotheses have been proposed as explanations for uneven diversification with many centering on key innovations or climate. while the pattern of the latitudinal diversity gradient, meaning the existence of greater species numbers centered around lower latitudes, is well documented, the mechanism behind the phenomenon is unclear, leading to a plethora of hypotheses. it is possible and even probable that no single hypothesis can explain such patterns. furthermore, these and other patterns could be the results of purely stochastic processes. however, evidence has shown that climatic factors are especially influential on the evolutionary history of many organismal groups. as nam psoraleeae is within temporal and geographical range of quaternary climate oscillations and exhibits a wide variety of climate-based habitats, we chose to address hypotheses involving climate change as possible effectors of diversification rate shifts. habitat shifts past climate change has impacted current climate regimes, creating a spectrum of climate-based habitats across the earth. movement of ancestral species from one habitat to another may also have spurred diversification, relegating climate-based habitat as a key innovation. climate regime shifts have been found to have impacted diversification across various organismal groups, including plants and lizards. hoita, orbexilum, and rupertia all reside in more mesic environments characterized by warm, moist, humid climates whereas psoralidium and pediomelum reside in more xeric habitats characterized by hot, dry environments such as the deserts of the intermountain and southwest united states. in testing habitat shifts from mesic to xeric habitats, our analyses suggest that past climate change affected diversification rates in nam psoraleeae more than shifts in present climate regimes as defined by current habitats. the shift to xeric habitats was not considered a key innovation in terms of increasing diversification rates within the group. instead, current habitats could be seen as relicts from the effects of past climate change on its influence on the diversification of north american psoraleeae. quaternary climate oscillations the recent origin of the main radiation of pediomelum around mya places its diversification in the middle of the pleistocene epoch of the quaternary era, a time period replete with dramatic glaciation-induced climate oscillations. these glacial fluxes include milankovitch cycles which consist of repeated formation and recession of glacial ice sheets having periodicities of to years, resulting in wide-ranging climate fluxes. it has been well documented that much of the earth's biodiversity and many species' evolutionary trajectories were drastically influenced by the climate fluxes brought on by quaternary glaciations. in north america, the sheer presence of ice sheets shifted climates and strongly influenced the range and distribution of species. ice formation may have influenced speciation via genetic differentiation through range fragmentations such as glacial refugia. upon glacial recession, new environments were formed and subsequently colonized through species expansion and adaptation. given the geographical and temporal distribution of nam psoraleeae, these cycles may have influenced its evolutionary history. indeed, the shift point for the increase in diversification along the phylogeny at mya and the lineage through time plot from our analyses support the hypothesis that quaternary climate fluxes, occurring between and mya, influenced the diversification of nam psoraleeae. geography the effects of these oscillations are not homogenous, however. the quaternary ice ages and climate oscillations affected geographic regions across north america differently. the southeast provided several refugia in which species could survive while the montane regions of the west were characterized by more varied environments such as ice, tundra, pluvial lakes, and deserts. in addition, elevation played an important role in species distributions as elevational ranges of plants were to 1200+ meters below modern ranges, an effect felt more readily in the topographically diverse regions of the west. upon glacial recession, the intermountain west offered a greater variety of habitats for colonization and adaptation as a result of topography. this is in line with cracraft's lithospheric complexity hypothesis, stating that geographic separation spurred by topographical changes promotes speciation. in addition to orogeny, glacial cycles contributed to topographical change, and thus are linked with the lithospheric hypothesis. geographical trends are seen in the diversity of psoraleeae. nam psoraleeae ranges throughout the united states: orbexilum is mostly confined to the southeast united states, with one species reaching into mexico, pediomelum is mostly found in central and intermountain states, hoita is endemic to california, rupertia ranges along the west coast, and psoralidium is spread over the great plains with one species endemic to utah. of interest is the seeming correlation of lineage age with southerly distribution. orbexilum is one of the oldest nam psoraleeae genera, diversifying earlier than most other lineages. this could be demonstrative of long-established ranges of orbexilum through pleistocene refugia in the southeast united states. of import also is the age of pediomelum as a whole, beginning its diversification around mya with the break off of subgenus leucocraspedon. the two leucocraspedon species exhibit the most southern distribution of all pediomelum species with p. rhombifolium being distributed throughout texas and into mexico while p. palmeri is almost exclusively found in mesoamerica. this southern distribution perhaps allowed this lineage to persist through the quaternary climate oscillations as the effects of glacial climate change were moderated toward the equator. the main radiation of pediomelum exhibits the greatest number of species with the lowest intraspecific variation, findings in accord with the lithospheric complexity hypothesis, geographic separation, ecological diversification due to habitat diversity, and range fragmentation brought on by quaternary climate shifts. the pattern of less diversity in the southeast and greater diversity in the southwest is also found in cicindela beetles and agelenopsis spiders, among others.  nam psoraleeae has proved useful for studying evolutionary divergence and the impact of climate change on diversification within a recent, rapid radiation. this work effectively demonstrates differential rates of diversification at the species and generic levels, trends manifested across taxonomic levels within fabaceae. an estimation of divergence dates within psoraleeae found that the main radiation of pediomelum, comprising the bulk of species richness in nam psoraleeae, originated to mya during the quaternary era, concurrent with the milankovitch climate oscillations brought on by glacial cycles. the hypothesis of a diversification rate shift in nam psoraleeae was confirmed by both topological and temporal methods. we demonstrated a large increase in diversification rate near the beginning of the pleistocene epoch. our findings suggest that the diversification rate increase in nam psoraleeae is not due to current climate regimes as represented by habitat, but rather imply that quaternary climate oscillations impacted the diversification of the group. we show nonrandom spatial diversification as well, with higher species diversity in the american southwest. nam psoraleeae diversification is a good example of how earthly dynamics including global climate change and topography work together to shape biodiversity. in a time when the issue of global climate change pervades society and science, this study aptly illustrates the control that global climate change can have over the evolutionary trajectories and fates of lineages and species. 
diversification rate shifts in nam psoraleeae are not due to current climate regimes as represented by habitat, but instead to past global climate change resulting from quaternary glaciations. nam psoraleeae represents a recent, rapid radiation with several genera originating during the pleistocene, to million years ago. this study outlines nam psoraleeae as a model example of a recent, rapid radiation. nam psoraleeae diversification is a good example of how earthly dynamics including global climate change and topography work together to shape biodiversity. we investigate the impact of climate change on diversification in this group by testing the hypothesis that a shift from mesic to xeric habitats acted as a key innovation and investigating diversification rate shifts along geologic time, discussing the impact of quaternary climate oscillations on diversification. past studies in the legume family have uncovered several evolutionary trends including differential mutation and diversification rates across varying taxonomic levels. the hypothesis that a climate regime shift from mesic to xeric habitats drove increased diversification in affected taxa was not supported. nonrandom spatial diversification also exists, with greater species richness in the american southwest. this group includes a sizable diversification within north america dated at approximately million years ago with skewed species distribution to the most recently derived genus, pediomelum, suggesting a diversification rate shift. we also test the hypothesis of a diversification rate shift within nam psoraleeae using topological and temporal methods. a shift in diversification rate is supported by both methods with a -fold increase suggested around million years ago followed by a -fold decrease years ago. timing of the diversification rate increase supports the hypothesis that glaciation-induced climate changes during the quaternary influenced diversification of the group.
 past studies in the legume family have uncovered several evolutionary trends including differential mutation and diversification rates across varying taxonomic levels. the legume tribe psoraleeae is shown herein to exemplify these trends at the generic and species levels. this group includes a sizable diversification within north america dated at approximately million years ago with skewed species distribution to the most recently derived genus, pediomelum, suggesting a diversification rate shift. we estimate divergence dates of north american psoraleeae using bayesian mcmc sampling in beast based on eight dna regions. we also test the hypothesis of a diversification rate shift within nam psoraleeae using topological and temporal methods. we investigate the impact of climate change on diversification in this group by testing the hypothesis that a shift from mesic to xeric habitats acted as a key innovation and investigating diversification rate shifts along geologic time, discussing the impact of quaternary climate oscillations on diversification.  nam psoraleeae represents a recent, rapid radiation with several genera originating during the pleistocene, to million years ago. a shift in diversification rate is supported by both methods with a -fold increase suggested around million years ago followed by a -fold decrease years ago. the hypothesis that a climate regime shift from mesic to xeric habitats drove increased diversification in affected taxa was not supported. timing of the diversification rate increase supports the hypothesis that glaciation-induced climate changes during the quaternary influenced diversification of the group. nonrandom spatial diversification also exists, with greater species richness in the american southwest.  this study outlines nam psoraleeae as a model example of a recent, rapid radiation. diversification rate shifts in nam psoraleeae are not due to current climate regimes as represented by habitat, but instead to past global climate change resulting from quaternary glaciations. nam psoraleeae diversification is a good example of how earthly dynamics including global climate change and topography work together to shape biodiversity. 