BACKGROUND
divergence times are of major interest for studies of evolutionary biology and historical biogeography, but also to researchers who focus on understanding various types of trait evolution, such as the development of chemical and genetic pathways, climatic niche and geographic range sizes, and morphological, ecological and behavioural characters.
with the recent publication of fully annotated genomes in the solanaceae  <cit> , genomic tools now exist for unravelling genetic mechanisms that control such traits and their development.
what is lacking, however, is a robust phylogenetic framework that encompasses species and generic diversity across the family in order to maximise the potential of these new data sources in a wider evolutionary context.
although several studies have focused on understanding evolution of particular characteristics in solanaceae in a phylogenetic context, including analyses of genome and chromosome evolution  <cit> , life history and polyploidy  <cit> , floral and fruit morphology  <cit> , gene family evolution and sub-functionalization  <cit> , and broad-scale biogeographic patterns  <cit> , only a single study has examined character evolution through time  <cit> .
a central problem has been the lack of a robust, densely sampled, dated molecular phylogeny for the entire family.
solanaceae are a particularly interesting angiosperm family, not only because they include many major crop species  and ornamentals , but also a number of taxa used as biological model systems   <cit> .
a taxonomic and phylogenetic framework is now available for the family at the generic level   <cit> .
separate studies have explored relationships at tribal  <cit>  and generic levels  <cit> , and many phylogenetic studies have focused on solanum, a genus that comprises nearly half of the species within the family  <cit> ; see  <cit>  for references prior to 2006].
some of these studies have resulted in taxonomic changes, such as re-circumscription of the previously distinct genera lycopersicon, cyphomandra, normania, and triguera as parts of a monophyletic solanum <cit> .
such changes, although sometimes disruptive in the short term, have helped to both stabilise names and provide a better evolutionary context for future studies.
although a relatively robust understanding of the major clades within solanaceae exists, a densely sampled species-level phylogeny is still lacking.
the most recent molecular systematic study focused on establishing major relationships within the family, but lacked depth in terms of sampling as it only included  <dig>  of a total c.  <dig>  solanaceae species  <cit> .
a larger phylogenetic analysis with  <dig> species by goldberg and colleagues  <cit>  focused on the evolution of breeding systems, and did not discuss details of topology or implications for family-wide systematics.
for solanum itself, the most recent phylogeny included only  <dig>  of the total c.  <dig>  species in the genus  <cit> .
because several phylogenetic studies at various taxonomic levels across solanaceae have since been published , there is now a large quantity of new sequence data available for a wider family-level analysis.
molecular divergence time analyses do not only depend on the availability of a robust, well-sampled phylogeny, but also require robust fossil calibration points  <cit> .
the solanaceae fossil record has never been fully reviewed, and only a few fossils have been used in molecular studies  <cit> .
these studies used fossils as calibration points without a careful comparison of fossil morphology in relation to extant diversity.
a recent survey of the earliest fossil record of the asterid clade, including solanaceae, highlighted the need to re-assess the earliest putative solanaceae fossils that could provide robust calibration points for the crown or stem node of the family  <cit> .
this study is part of a collaborative approach to studying the taxonomy and phylogeny of the solanaceae.
here we present a densely sampled phylogenetic study of the family coupled with a molecular dating analysis with fossil calibrations.
we review all known seed fossils in the family, and assess them for identity, age, and phylogenetic position.
we then use all available sequences for seven dna loci found in genbank with nearly all genera and  <dig>  species represented.
a dating analysis is run using an uncorrelated relaxed molecular clock model within a bayesian framework with direct fossil calibrations.
the resulting time-calibrated phylogeny offers important insights into the evolution of the family at different taxonomic levels, and a robust platform for future evolutionary studies.
RESULTS
fossil review
a total of  <dig> fossil records previously assigned to solanaceae were found in the literature .
these included  <dig> seed fossils, one leaf fossil, five flowers, two wood and three pollen fossils .
none of the leaf or flower fossils showed any distinct morphological characters that allowed us to definitely assign them to the family.
of the two wood fossils, solanumxylum paranensis can be clearly assigned to solanaceae based on a large number of anatomical characters such as para- and apotracheal axial parenchyma that is diffuse in aggregates, simple perforation plates, bordered and alternate intervessel pits, homocellular rays, fibres that are polygonal and quandrangular in section, and the presence of septate fibers   <cit> .
the other wood fossil shows no specific characters of solanaceae except those common to solanaceae and asteraceae and lacks axial parenchyma; we do not consider this a member of solanaceae   <cit> .
of the two pollen records, the classification of datura cf.
discolor awaits further examination, since no description or illustration of the fossil was provided in the original publication .
a pollen fossil-taxon from california, based on two poorly preserved specimens of 3-colporate, 5-colpate, prolate shaped grains with striate ornamentation, resembles pollen grains of lycium, nolana, and hyoscyamus <cit> .
similar characters appear in the pollen of the unrelated genera brucea  and skimmia   <cit> , and hence we have not assigned this pollen fossil-taxon to solanaceae for our analysis .
list of all known fossil records of solanaceae.
fossils highlighted in bold were scanned using high-resolution x-ray computed tomography .
for full details of the fossil and references to primary sources, see additional file  <dig>
<dig> no description nor illustration or proper reference to specimen.
<dig> reference not located.
the putative solanaceae seed fossils were analysed using a combination of characteristics known from clades within the family  <cit> :  seeds flattened,  circular to reniform in shape,  hilum sub-laterally or laterally positioned, and  testa cells sinuate-margined.
we assigned seeds with all four of these characters to the subfamily solanoideae , while those with some but not all of these were assigned to the family as a whole    <cit> .
the currently recognised earliest fossils assignable to solanaceae include two seed fossils from eocene europe: solanispermum reniforme recorded from various beds from southern england  <cit> , and solanum arnense, a fossil-taxon described based on a few specimens found from the lower bagshot   <cit> .
neither of these shows the combination of flat seeds with sinuate-margined testa cells, a unique combination that could tie them to the tribe solanoideae.
the flattened seeds of solanispermum reniforme lack sinuate margined testa cells, and solanum arnense seeds show the characteristic testa cells but seeds are round rather than flattened.
hence, we consider these fossils as earliest evidence of solanaceae and the presence of the family in eocene europe, but do not assign them to any particular clade within the family.
seeds of the fossil-taxon cantisolanum daturoides <cit>  have previously been cited as the oldest known solanaceae fossil by some authors  <cit>  but doubtfully a member of the family by others  <cit> .
results from a ct-scanning study have shown that this cantisolanum seed is anatropous and does not belong to solanaceae, but has likely affinities to the monocot family philydraceae .
solanaceae phylogeny
our final supermatrix had a taxon coverage density of  <dig> , and included  <dig>  species of solanaceae, representing all but two genera  and 40% of total species within the family, including 34% sampling of species within the large genus solanum.
two plastid regions, ndhf and trnl-f, were available for all genera except darcyanthus and capsicophysalis and the plastid and nuclear regions its, waxy and trnl-f were the most densely sampled regions at the species level .
the matrix included a total of  <dig>  variable characters, with an aligned length of  <dig>  bp .
a total of  <dig>  bp were excluded from analyses due to ambiguous alignment  resulting in a matrix of  <dig>  bp .
proportionately, waxy  and ndhf  contributed most pi  characters .
the relatively little-used plastid region trns-g showed a surprising number of pi characters , considering it had relatively poor taxon coverage density , compared to trnl-f which had a coverage of  <dig>  but only  <dig> % of total pi characters .
the final matrix included  <dig> % missing data .
at the species level, there was an average of  <dig> % missing data, as measured by number of base pairs, but only  <dig> % when measured in terms of pi characters expected from the missing regions.
list of the seven regions used for building the solanaceae mega-phylogeny, showing level of sampling at specific and generic levels, aligned length, and proportions of parsimony informative  and variable sites per region.
values have been calculated without outgroups.
1ingroup only, excluded regions omitted.
the resolved maximum likelihood topology shows strong support for all previously identified major clades within solanaceae  <cit> , and increased node support is observed particularly within solanum .
only major clades and their relationships are discussed here due to the fact that our analyses only accounted for incongruence issues amongst data sets between major clades rather than at shallow taxonomic levels.
we encourage readers to refer back to available clade-specific studies for detailed species-level phylogenies ; these studies have incorporated larger sets of markers than used here, incorporate methods that test/account for gene tree – species tree incongruence, and discuss issues that could have led to any detected incongruences between gene trees such as polyploidy and/or hybridisation, and incomplete lineage sorting.
the branching order at the base of solanaceae is not well defined, similar to the findings of olmstead et al.
<cit> , and four groups are identified as the first branching taxa: schizanthus, duckeodendron, the previously unplaced reyesia, and the tribe goetzeoideae .
reyesia has been previously associated with salpiglossis <cit> , but is here placed with goetzeoideae and duckeodendron .
the previously unsampled genera heteranthia, trianaea and schraderanthus are placed within schwenckieae, juanulloeae, and physalinae, respectively .
the informally named x =  <dig> clade is here recovered with strong support and nicotianoideae is resolved as sister to the rest of the clade .
within the physalinae, work is clearly needed to delimit monophyletic genera .
two closely related genera, larnax and deprea, are resolved as sister to withaninae, in agreement with morphology .
these genera have been linked with iochrominae in some molecular analyses  <cit> , but considered distant outgroups of iochrominae by others  <cit> .
the molecular data support the treatment of schraderanthus as distinct from leucophysalis <cit> , and schraderanthus is here found as sister to brachistus + witheringia .
within solanum, all  <dig> major clades identified by weese & bohs  <cit>  are recovered, with nearly fully resolved relationships among them .
the thelopodium clade is resolved as the first branching group, and the remaining solanum species are divided into two strongly supported clades.
clade i comprises all non-spiny, often herbaceous  species without stellate hairs, but also includes woody climbers  and some shrubby species .
clade ii comprises species that are often shrubs or small trees , many with prickles and/or stellate hairs .
within clade i, which includes a total of c.  <dig> known species, two clear clades are resolved:  the potato clade, with regmandra clade as the first branching group, and  clade m, including morelloid, dulcamaroid, archaesolanum, normania, and the african non-spiny clades .
relationships within clade m are well resolved and highly supported, revealing the position of the african non-spiny clade as distinct from and not closely related to the dulcamaroid clade, despite their morphological similarities such as a twining habit and twisting petioles  <cit> .
within the potato clade, relationships are equally well resolved: section petota is resolved as sister to a group comprising the tomato clade plus a set of smaller early-branching clades .
the regmandra clade, a group of  <dig> species whose centre of diversity is the hyper-arid atacama desert, is here resolved as part of the potato clade for the first time , a result supported by morphology  <cit> .
relationships within clade ii are less well-resolved.
the clade consists of c.  <dig> mostly woody species, and includes the large leptostemonum clade known as "spiny solanums".
there is moderate support for s. clandestinum + s. mapiriense as sister to the rest of clade ii .
relationships within the large leptostemonum clade remain relatively unresolved, but all  <dig> major clades found in previous analyses  <cit>  are supported.
a set of previously unplaced species, s. crotonoides, s. hayesii, and s. multispinum, are resolved sequentially as sister to the torva clade , although on morphological grounds s. hayesii would be a member of the torva clade.
molecular dating
the general topology of the bayesian maximum clade credibility tree matched that of the best scoring maximum likelihood tree with similar levels of support for major clades .
the only topological difference, although not a hard incongruence, was observed at the base of solanaceae: bayesian analyses resolved schwenckieae as the first branching group within the family, while the base of the tree remained largely unresolved in the maximum likelihood topology.
results from pathd <dig> gave generally similar ages as those from the beast analysis .
a notable trend is that beast ages were consistently younger especially towards the early-branching nodes .
the younger ages obtained from the beast analysis reflect that diversification rates across solanaceae have been non-linear especially towards the base of the tree, and/or that extinction and speciation rates have varied across the tree.
we will focus our discussion on the beast results, which we consider to be more robust due to the more realistic model assumptions used, including the relaxed molecular clock model that accounts for rate variation across lineages, as well as birth-death tree model accounting for extinction  <cit> .
minimum age estimates produced based on both pathd <dig> and beast analyses for solanaceae and its major clades.
clades refer to crown nodes unless otherwise indicated.
ages in brackets are 95% confidence intervals for beast.
a densely sampled, time calibrated solanaceae mega-phylogeny with associated confidence intervals can be found in additional file  <dig>
the beast results place the stem age of solanaceae at c.  <dig> million years ago : 46–54 ma), and the crown node at c. 30 ma  .
the crown node of the x =  <dig> clade, which is the split between nicotiana and solanum, was estimated to be c.  <dig> my old  .
the solanoideae began diversifying c. 21 ma  .
solanum, a genus which includes nearly half of the total species diversity in the family, split from jaltomata c. 17 ma  and started diversifying c. 16 ma  .
the solanum – capsicum split, corresponding to the most common ancestor of solanum & physalis, occurred c. 19 ma  .
within solanum, major splits include tomato – potato c. 8 ma  and the eggplant –tomato/potato lineages corresponding to the clade i – clade ii split c. 14 ma  .
crown node age estimates show that section petota, which includes all cultivated potatoes, started diversifying c. 7 ma , section lycopersicon, which includes the cultivated tomato, c. 2 ma , the group containing all species of cultivated eggplants  c. 3 ma , and the group  including all cultivated pepper species c. 3 ma  .
discussion
phylogenetic relationships within solanaceae
although individual studies have contributed significantly to a better understanding of the systematics and evolution of the family at generic and tribal levels, our results bring together data from a large number of studies into a single analysis, and present a coherent view on the current systematic knowledge of this diverse family and its major clades.
our analyses support all of the major clades previously identified within solanaceae  <cit> , solanum <cit>  and the leptostemonum clade of solanum <cit> .
all of these major clades within the family are now strongly supported, and furthermore, our results reveal strongly-supported relationships between the major clades of the mega-diverse genus solanum, strengthening the backbone.
the increased resolution in the current phylogeny can be attributed to both the increased sampling of markers as well as species.
in the quest for better resolved phylogenies, studies often seek large amounts of sequence data, but it is now well established that increased species sampling can have an equally positive effect on phylogenetic resolution and accuracy  <cit> .
the addition of more species to a data set has the effect of splitting long branches and detecting multiple substitutions, as well as resolving phylogenetic conflict, improving parameter estimation, and making inferences less dependent on particular evolutionary models  <cit> .
in our approach we chose to maximise species sampling, while minimising missing data by choosing only the most densely sampled markers available.
this approach generally boosted resolution without introducing any of the significant negative effects that large amounts of missing data can have on phylogeny estimation.
our study presented here is a significant step forward in working towards a fully sampled species-level phylogeny for solanaceae.
a previous study by goldberg et al.
<cit>  included  <dig> species but did not present a fully annotated molecular phylogeny that would allow an analysis of systematic relationships within the family.
with >  <dig>  species now covered, the current phylogeny includes 40% of known species and all genera of solanaceae, except the monospecific and recently segregated darcyanthus and capsicophysalis.
this is a substantial improvement on previous studies, and our current phylogeny is one of the best sampled family-level studies in angiosperms e.g.,  <cit> .
the sampling is now adequate to test for generic monophyly in previously poorly sampled groups.
although the number of genera is becoming stable with  <dig> currently recognised genera in solanaceae , our analyses support previous results in identifying a set of groups where generic re-evaluation will be necessary, including lycianthes/capsicum, the genera in the physalineae   <cit> , deprea/larnax, the iochrominae  <cit> , and the australian endemics in the anthocercideae ; many of these clusters of generic problems have been identified by previous authors.
broader level relationships within solanaceae and solanum, as well as generic delimitations and problems identified in previous studies are supported by our species-rich dataset.
relationships between some of the major clades remain unresolved, however, most notably those at the base of the family and within the solanoideae, and the leptostemonum clade of solanum.
resolving these nodes will be a priority in order to better understand evolution of some particularly complex traits, such as chromosome evolution.
for example, resolving the sister group to the x =  <dig> clade, as well as the first branching taxa within solanaceae, would allow us to determine the ancestral base chromosome number in the family and to fully understand directionality of chromosome evolution.
despite the increased resolution introduced by the use of more sequence data and higher species-level sampling, our results do not show any improvement in the resolution in these critical nodes.
more genes will be needed to resolve these relationships, but the question remains which genes should be used.
highly variable nuclear loci, such as cosii markers already used in solanaceae  <cit> , and the ppr genes used in families within the related asterid order lamiales  <cit> , present the most promising candidates.
the widely sequenced regions its, waxy, ndhf, and trnsg are the most variable across the solanaceae and species-level sampling using these regions should be increased.
the traditionally used plastid marker trnt-f, which is relatively slowly evolving within solanaceae, is known to include pseudogenes in solanum <cit>  and care should be taken when using this region in phylogenetic studies.
solanaceae fossil record
a few fossils have been used in previous molecular dating studies of solanaceae, but without re-evaluation of fossil morphology and hence their placement within the phylogeny  <cit> .
as revealed by our literature review, a relatively large record for the family exists.
the most usable evidence comes from fossil solanaceae seeds, the oldest of which are from eocene europe , with a sharp increase in the number of seed morphotypes observed towards the pleistocene.
the fossil seeds can be divided into two sets:  seeds showing four morphological characters present in the extant members of the solanoideae, and  seeds that bear resemblance to the family in general but cannot be assigned to more specific clades within it because they lack the unique combination of seed flattening and presence of sinuate-margined testa cells.
although some of these fossils have been described with names associated with extant species and/or genera, our morphological review shows that none of them show unique morphological characters that can be used to place them to any extant genera.
we consider the placement of these fossils on terminal nodes as has been done by previous authors  <cit>  unjustified.
all of the fossils we were able to unambiguously identify as solanaceae are from eocene europe, where none of the first branching lineages of the family occur.
south america is the centre of diversity of extant solanaceae, and all of the early diverging lineages are exclusively found in the new world.
this suggests that the fossil record of the family is still far from complete, and that further studies on south american fossils might reveal crucial evidence with respect to the timing of diversification in solanaceae.
a promising avenue for future fossil studies would be to carefully evaluate wood fossil records, especially cretaceous-eocene material from the area in which the early-branching lineages all now occur  <cit> .
dates for solanaceae
our study is the largest bayesian molecular dating analysis executed to date in terms of taxon sampling.
most previous studies have used bayesian dating methods after pruning their original, large phylogenetic datasets largely due to an a priori assumption that bayesian methods cannot cope with datasets with > <dig> terminals e.g.,  <cit> .
our study with  <dig>  species and > <dig>  bp of sequence data demonstrates that large matrices with > <dig> terminals can be analysed using bayesian dating methods.
further studies are needed, however, to fully explore best methods for analysing large datasets with the currently available dating methods that implement relaxed molecular clock models required for analyses of diverse clades where rates are expected to vary  <cit> .
such studies should focus on exploring trade-offs between number of taxa, complexity of models and partitions used in order to fully understand limitations and potential error sources in large scale analyses.
in our dating analysis, we followed the recent recommended best practice guidelines for fossil calibration  <cit>  and placed fossil calibrations at stem nodes of the most inclusive extant groups using apomorphy-based morphological assignment.
morphological evidence from the seed fossils only allowed assignment to the broad groups solanoideae or solanaceae as a whole.
fossils provide only minimum age estimates for the nodes they are assigned, and hence results from our dating analysis where fossil calibrations were used should be considered as minimum age estimates.
we further biased our results towards younger ages by assigning the oldest known fossils of solanaceae to the stem node of the family rather than to more specific nodes within solanaceae due to lack of morphological and anatomical characters that could be used to assign them to more specific nodes.
there is always a possibility, however, that these seeds represent more specific clades within solanaceae, which would push back age estimates for the family.
currently, the earliest fossil evidence for the family comes from eocene europe, but based on biogeographic analyses, the crown group of solanaceae is thought to have originated and first diversified in south america  <cit> .
total evidence analysis, where fossils are placed as terminal taxa in the dating analysis using both molecular and morphological data matrix, could help in exploring the robustness of fossil placement  <cit> , but as pointed above, the lack of characters in the solanaceae seed fossils does not currently permit such analyses.
the most promising avenue in strengthening the dating analysis would be in finding further fossil records .
this would increase the number of fossil calibration points and allow the use of cross-validation methods  <cit> .
the rate of molecular evolution in plants has been found to correlate with life history traits, whereby longer living species show consistently lower substitution rates compared to shorter living species  <cit> .
molecular clock models should incorporate such rate variation, especially in groups such as solanaceae which include a range of growth and life forms.
our dating analyses did not incorporate such models, although the model used in our bayesian analysis allows rates to vary between lineages independently.
the lack of such models in our analyses implies that the age of herbaceous, shorter lived plants  will be systematically overestimated, while ages in dominantly woody clades  will be consistently underestimated.
future studies should explore how molecular clock models that account for rate variation due to life history traits could be implemented.
previous studies have produced a wide range of estimates for the stem age of the family, ranging from 34–85 ma  <cit> , but none of these studies included dense sampling within the family nor used robust solanaceae-specific fossil calibrations.
paape et al.
<cit>  analysed divergence times within solanaceae but with a small dataset consisting of  <dig> species only.
this study was based on three fossil calibration points without re-assessment or morphological study of the original fossils, and estimated solanaceae stem age to have diverged 62 ma   <cit> .
the oldest estimates for the family stem node age come from earlier molecular studies which used calibration points with more simplistic dating methods   <cit> , while the most recent molecular dating study of angiosperms by bell et al.
<cit>  who used  <dig> fossil calibrations across the tree and a relaxed molecular clock model, estimated the solanaceae stem node to have diverged c. 59 ma .
our results, which we consider as minimum ages, are broadly consistent with bell et al.
<cit>  in estimating the stem node of solanaceae to date back to c. 49 ma .
the age of the major splits within the family has been of interest to various fields, including studies on chromosomal  <cit>  and genome evolution  <cit> .
our minimum age estimates for the major splits between tomato – potato , eggplant – tomato/potato , solanum – capsicum , and solanum – nicotiana  are consistent with the age estimates produced in previous studies without fossil calibrations using much sparser sampling and more simplistic molecular clock models  <cit> .
our results for the nicotiana – symonanthus split  corroborate results obtained using island age   <cit>  and those calculated using paralogy-free subtree analysis   <cit> .
our results presented here suggest that the rate of chromosomal and genome evolution within solanaceae has been marginally slower at least within particular lineages than previously thought.
with the densely sampled chronogram presented in this study, a more detailed analysis of chromosomal evolution at the species level could now be performed in the solanaceae to study rate differences and drivers of chromosomal changes such as environmental or life history factors.
similarly, morphological characters such as fruit type  <cit>  could be analysed in relation to diversification rates to identify whether particular morphological traits are associated with speciation rate shifts in solanaceae.
CONCLUSIONS
despite much focus on character and trait evolution within solanaceae, little has been known about the origin of traits in the family in terms of time.
we present here minimum age estimates and associated confidence intervals for the entire solanaceae using a species-rich dataset comprising almost half of the species diversity within the family.
this densely sampled chronogram will provide the basis for unravelling the tempo and mode of evolution of many of the much-studied and complex traits in this diverse and economically important family such as self-incompatibility, fruit type, cold and salt tolerance, disease resistance, chromosomal re-arrangements, genome size, and gene sub-functionalization.

 divergence times are of major interest for studies of evolutionary biology and historical biogeography, but also to researchers who focus on understanding various types of trait evolution, such as the development of chemical and genetic pathways, climatic niche and geographic range sizes, and morphological, ecological and behavioural characters. with the recent publication of fully annotated genomes in the solanaceae, genomic tools now exist for unravelling genetic mechanisms that control such traits and their development. what is lacking, however, is a robust phylogenetic framework that encompasses species and generic diversity across the family in order to maximise the potential of these new data sources in a wider evolutionary context. although several studies have focused on understanding evolution of particular characteristics in solanaceae in a phylogenetic context, including analyses of genome and chromosome evolution, life history and polyploidy, floral and fruit morphology, gene family evolution and sub-functionalization, and broad-scale biogeographic patterns, only a single study has examined character evolution through time. a central problem has been the lack of a robust, densely sampled, dated molecular phylogeny for the entire family. solanaceae are a particularly interesting angiosperm family, not only because they include many major crop species and ornamentals, but also a number of taxa used as biological model systems. a taxonomic and phylogenetic framework is now available for the family at the generic level. separate studies have explored relationships at tribal and generic levels, and many phylogenetic studies have focused on solanum, a genus that comprises nearly half of the species within the family; see for references prior to 2006]. some of these studies have resulted in taxonomic changes, such as re-circumscription of the previously distinct genera lycopersicon, cyphomandra, normania, and triguera as parts of a monophyletic solanum. such changes, although sometimes disruptive in the short term, have helped to both stabilise names and provide a better evolutionary context for future studies. although a relatively robust understanding of the major clades within solanaceae exists, a densely sampled species-level phylogeny is still lacking. the most recent molecular systematic study focused on establishing major relationships within the family, but lacked depth in terms of sampling as it only included of a total c. solanaceae species. a larger phylogenetic analysis with species by goldberg and colleagues focused on the evolution of breeding systems, and did not discuss details of topology or implications for family-wide systematics. for solanum itself, the most recent phylogeny included only of the total c. species in the genus. because several phylogenetic studies at various taxonomic levels across solanaceae have since been published, there is now a large quantity of new sequence data available for a wider family-level analysis. molecular divergence time analyses do not only depend on the availability of a robust, well-sampled phylogeny, but also require robust fossil calibration points. the solanaceae fossil record has never been fully reviewed, and only a few fossils have been used in molecular studies. these studies used fossils as calibration points without a careful comparison of fossil morphology in relation to extant diversity. a recent survey of the earliest fossil record of the asterid clade, including solanaceae, highlighted the need to re-assess the earliest putative solanaceae fossils that could provide robust calibration points for the crown or stem node of the family. this study is part of a collaborative approach to studying the taxonomy and phylogeny of the solanaceae. here we present a densely sampled phylogenetic study of the family coupled with a molecular dating analysis with fossil calibrations. we review all known seed fossils in the family, and assess them for identity, age, and phylogenetic position. we then use all available sequences for seven dna loci found in genbank with nearly all genera and species represented. a dating analysis is run using an uncorrelated relaxed molecular clock model within a bayesian framework with direct fossil calibrations. the resulting time-calibrated phylogeny offers important insights into the evolution of the family at different taxonomic levels, and a robust platform for future evolutionary studies.  fossil review a total of fossil records previously assigned to solanaceae were found in the literature. these included seed fossils, one leaf fossil, five flowers, two wood and three pollen fossils. none of the leaf or flower fossils showed any distinct morphological characters that allowed us to definitely assign them to the family. of the two wood fossils, solanumxylum paranensis can be clearly assigned to solanaceae based on a large number of anatomical characters such as para- and apotracheal axial parenchyma that is diffuse in aggregates, simple perforation plates, bordered and alternate intervessel pits, homocellular rays, fibres that are polygonal and quandrangular in section, and the presence of septate fibers. the other wood fossil shows no specific characters of solanaceae except those common to solanaceae and asteraceae and lacks axial parenchyma; we do not consider this a member of solanaceae. of the two pollen records, the classification of datura cf. discolor awaits further examination, since no description or illustration of the fossil was provided in the original publication. a pollen fossil-taxon from california, based on two poorly preserved specimens of 3-colporate, 5-colpate, prolate shaped grains with striate ornamentation, resembles pollen grains of lycium, nolana, and hyoscyamus. similar characters appear in the pollen of the unrelated genera brucea and skimmia, and hence we have not assigned this pollen fossil-taxon to solanaceae for our analysis. list of all known fossil records of solanaceae. fossils highlighted in bold were scanned using high-resolution x-ray computed tomography. for full details of the fossil and references to primary sources, see additional file no description nor illustration or proper reference to specimen. reference not located. the putative solanaceae seed fossils were analysed using a combination of characteristics known from clades within the family: seeds flattened, circular to reniform in shape, hilum sub-laterally or laterally positioned, and testa cells sinuate-margined. we assigned seeds with all four of these characters to the subfamily solanoideae, while those with some but not all of these were assigned to the family as a whole. the currently recognised earliest fossils assignable to solanaceae include two seed fossils from eocene europe: solanispermum reniforme recorded from various beds from southern england, and solanum arnense, a fossil-taxon described based on a few specimens found from the lower bagshot. neither of these shows the combination of flat seeds with sinuate-margined testa cells, a unique combination that could tie them to the tribe solanoideae. the flattened seeds of solanispermum reniforme lack sinuate margined testa cells, and solanum arnense seeds show the characteristic testa cells but seeds are round rather than flattened. hence, we consider these fossils as earliest evidence of solanaceae and the presence of the family in eocene europe, but do not assign them to any particular clade within the family. seeds of the fossil-taxon cantisolanum daturoides have previously been cited as the oldest known solanaceae fossil by some authors but doubtfully a member of the family by others. results from a ct-scanning study have shown that this cantisolanum seed is anatropous and does not belong to solanaceae, but has likely affinities to the monocot family philydraceae. solanaceae phylogeny our final supermatrix had a taxon coverage density of, and included species of solanaceae, representing all but two genera and 40% of total species within the family, including 34% sampling of species within the large genus solanum. two plastid regions, ndhf and trnl-f, were available for all genera except darcyanthus and capsicophysalis and the plastid and nuclear regions its, waxy and trnl-f were the most densely sampled regions at the species level. the matrix included a total of variable characters, with an aligned length of  bp. a total of  bp were excluded from analyses due to ambiguous alignment resulting in a matrix of  bp. proportionately, waxy and ndhf contributed most pi characters. the relatively little-used plastid region trns-g showed a surprising number of pi characters, considering it had relatively poor taxon coverage density, compared to trnl-f which had a coverage of but only % of total pi characters. the final matrix included % missing data. at the species level, there was an average of % missing data, as measured by number of base pairs, but only % when measured in terms of pi characters expected from the missing regions. list of the seven regions used for building the solanaceae mega-phylogeny, showing level of sampling at specific and generic levels, aligned length, and proportions of parsimony informative and variable sites per region. values have been calculated without outgroups. 1ingroup only, excluded regions omitted. the resolved maximum likelihood topology shows strong support for all previously identified major clades within solanaceae, and increased node support is observed particularly within solanum. only major clades and their relationships are discussed here due to the fact that our analyses only accounted for incongruence issues amongst data sets between major clades rather than at shallow taxonomic levels. we encourage readers to refer back to available clade-specific studies for detailed species-level phylogenies; these studies have incorporated larger sets of markers than used here, incorporate methods that test/account for gene tree – species tree incongruence, and discuss issues that could have led to any detected incongruences between gene trees such as polyploidy and/or hybridisation, and incomplete lineage sorting. the branching order at the base of solanaceae is not well defined, similar to the findings of olmstead et al., and four groups are identified as the first branching taxa: schizanthus, duckeodendron, the previously unplaced reyesia, and the tribe goetzeoideae. reyesia has been previously associated with salpiglossis, but is here placed with goetzeoideae and duckeodendron. the previously unsampled genera heteranthia, trianaea and schraderanthus are placed within schwenckieae, juanulloeae, and physalinae, respectively. the informally named x =  clade is here recovered with strong support and nicotianoideae is resolved as sister to the rest of the clade. within the physalinae, work is clearly needed to delimit monophyletic genera. two closely related genera, larnax and deprea, are resolved as sister to withaninae, in agreement with morphology. these genera have been linked with iochrominae in some molecular analyses, but considered distant outgroups of iochrominae by others. the molecular data support the treatment of schraderanthus as distinct from leucophysalis, and schraderanthus is here found as sister to brachistus + witheringia. within solanum, all major clades identified by weese & bohs are recovered, with nearly fully resolved relationships among them. the thelopodium clade is resolved as the first branching group, and the remaining solanum species are divided into two strongly supported clades. clade i comprises all non-spiny, often herbaceous species without stellate hairs, but also includes woody climbers and some shrubby species. clade ii comprises species that are often shrubs or small trees, many with prickles and/or stellate hairs. within clade i, which includes a total of c. known species, two clear clades are resolved: the potato clade, with regmandra clade as the first branching group, and clade m, including morelloid, dulcamaroid, archaesolanum, normania, and the african non-spiny clades. relationships within clade m are well resolved and highly supported, revealing the position of the african non-spiny clade as distinct from and not closely related to the dulcamaroid clade, despite their morphological similarities such as a twining habit and twisting petioles. within the potato clade, relationships are equally well resolved: section petota is resolved as sister to a group comprising the tomato clade plus a set of smaller early-branching clades. the regmandra clade, a group of species whose centre of diversity is the hyper-arid atacama desert, is here resolved as part of the potato clade for the first time, a result supported by morphology. relationships within clade ii are less well-resolved. the clade consists of c. mostly woody species, and includes the large leptostemonum clade known as "spiny solanums". there is moderate support for s. clandestinum + s. mapiriense as sister to the rest of clade ii. relationships within the large leptostemonum clade remain relatively unresolved, but all major clades found in previous analyses are supported. a set of previously unplaced species, s. crotonoides, s. hayesii, and s. multispinum, are resolved sequentially as sister to the torva clade, although on morphological grounds s. hayesii would be a member of the torva clade. molecular dating the general topology of the bayesian maximum clade credibility tree matched that of the best scoring maximum likelihood tree with similar levels of support for major clades. the only topological difference, although not a hard incongruence, was observed at the base of solanaceae: bayesian analyses resolved schwenckieae as the first branching group within the family, while the base of the tree remained largely unresolved in the maximum likelihood topology. results from pathd gave generally similar ages as those from the beast analysis. a notable trend is that beast ages were consistently younger especially towards the early-branching nodes. the younger ages obtained from the beast analysis reflect that diversification rates across solanaceae have been non-linear especially towards the base of the tree, and/or that extinction and speciation rates have varied across the tree. we will focus our discussion on the beast results, which we consider to be more robust due to the more realistic model assumptions used, including the relaxed molecular clock model that accounts for rate variation across lineages, as well as birth-death tree model accounting for extinction. minimum age estimates produced based on both pathd and beast analyses for solanaceae and its major clades. clades refer to crown nodes unless otherwise indicated. ages in brackets are 95% confidence intervals for beast. a densely sampled, time calibrated solanaceae mega-phylogeny with associated confidence intervals can be found in additional file the beast results place the stem age of solanaceae at c. million years ago: 46–54 ma), and the crown node at c. 30 ma. the crown node of the x =  clade, which is the split between nicotiana and solanum, was estimated to be c. my old. the solanoideae began diversifying c. 21 ma. solanum, a genus which includes nearly half of the total species diversity in the family, split from jaltomata c. 17 ma and started diversifying c. 16 ma. the solanum – capsicum split, corresponding to the most common ancestor of solanum & physalis, occurred c. 19 ma. within solanum, major splits include tomato – potato c. 8 ma and the eggplant –tomato/potato lineages corresponding to the clade i – clade ii split c. 14 ma. crown node age estimates show that section petota, which includes all cultivated potatoes, started diversifying c. 7 ma, section lycopersicon, which includes the cultivated tomato, c. 2 ma, the group containing all species of cultivated eggplants c. 3 ma, and the group including all cultivated pepper species c. 3 ma. discussion phylogenetic relationships within solanaceae although individual studies have contributed significantly to a better understanding of the systematics and evolution of the family at generic and tribal levels, our results bring together data from a large number of studies into a single analysis, and present a coherent view on the current systematic knowledge of this diverse family and its major clades. our analyses support all of the major clades previously identified within solanaceae, solanum and the leptostemonum clade of solanum. all of these major clades within the family are now strongly supported, and furthermore, our results reveal strongly-supported relationships between the major clades of the mega-diverse genus solanum, strengthening the backbone. the increased resolution in the current phylogeny can be attributed to both the increased sampling of markers as well as species. in the quest for better resolved phylogenies, studies often seek large amounts of sequence data, but it is now well established that increased species sampling can have an equally positive effect on phylogenetic resolution and accuracy. the addition of more species to a data set has the effect of splitting long branches and detecting multiple substitutions, as well as resolving phylogenetic conflict, improving parameter estimation, and making inferences less dependent on particular evolutionary models. in our approach we chose to maximise species sampling, while minimising missing data by choosing only the most densely sampled markers available. this approach generally boosted resolution without introducing any of the significant negative effects that large amounts of missing data can have on phylogeny estimation. our study presented here is a significant step forward in working towards a fully sampled species-level phylogeny for solanaceae. a previous study by goldberg et al. included species but did not present a fully annotated molecular phylogeny that would allow an analysis of systematic relationships within the family. with >  species now covered, the current phylogeny includes 40% of known species and all genera of solanaceae, except the monospecific and recently segregated darcyanthus and capsicophysalis. this is a substantial improvement on previous studies, and our current phylogeny is one of the best sampled family-level studies in angiosperms e.g.,. the sampling is now adequate to test for generic monophyly in previously poorly sampled groups. although the number of genera is becoming stable with currently recognised genera in solanaceae, our analyses support previous results in identifying a set of groups where generic re-evaluation will be necessary, including lycianthes/capsicum, the genera in the physalineae, deprea/larnax, the iochrominae, and the australian endemics in the anthocercideae; many of these clusters of generic problems have been identified by previous authors. broader level relationships within solanaceae and solanum, as well as generic delimitations and problems identified in previous studies are supported by our species-rich dataset. relationships between some of the major clades remain unresolved, however, most notably those at the base of the family and within the solanoideae, and the leptostemonum clade of solanum. resolving these nodes will be a priority in order to better understand evolution of some particularly complex traits, such as chromosome evolution. for example, resolving the sister group to the x =  clade, as well as the first branching taxa within solanaceae, would allow us to determine the ancestral base chromosome number in the family and to fully understand directionality of chromosome evolution. despite the increased resolution introduced by the use of more sequence data and higher species-level sampling, our results do not show any improvement in the resolution in these critical nodes. more genes will be needed to resolve these relationships, but the question remains which genes should be used. highly variable nuclear loci, such as cosii markers already used in solanaceae, and the ppr genes used in families within the related asterid order lamiales, present the most promising candidates. the widely sequenced regions its, waxy, ndhf, and trnsg are the most variable across the solanaceae and species-level sampling using these regions should be increased. the traditionally used plastid marker trnt-f, which is relatively slowly evolving within solanaceae, is known to include pseudogenes in solanum and care should be taken when using this region in phylogenetic studies. solanaceae fossil record a few fossils have been used in previous molecular dating studies of solanaceae, but without re-evaluation of fossil morphology and hence their placement within the phylogeny. as revealed by our literature review, a relatively large record for the family exists. the most usable evidence comes from fossil solanaceae seeds, the oldest of which are from eocene europe, with a sharp increase in the number of seed morphotypes observed towards the pleistocene. the fossil seeds can be divided into two sets: seeds showing four morphological characters present in the extant members of the solanoideae, and seeds that bear resemblance to the family in general but cannot be assigned to more specific clades within it because they lack the unique combination of seed flattening and presence of sinuate-margined testa cells. although some of these fossils have been described with names associated with extant species and/or genera, our morphological review shows that none of them show unique morphological characters that can be used to place them to any extant genera. we consider the placement of these fossils on terminal nodes as has been done by previous authors unjustified. all of the fossils we were able to unambiguously identify as solanaceae are from eocene europe, where none of the first branching lineages of the family occur. south america is the centre of diversity of extant solanaceae, and all of the early diverging lineages are exclusively found in the new world. this suggests that the fossil record of the family is still far from complete, and that further studies on south american fossils might reveal crucial evidence with respect to the timing of diversification in solanaceae. a promising avenue for future fossil studies would be to carefully evaluate wood fossil records, especially cretaceous-eocene material from the area in which the early-branching lineages all now occur. dates for solanaceae our study is the largest bayesian molecular dating analysis executed to date in terms of taxon sampling. most previous studies have used bayesian dating methods after pruning their original, large phylogenetic datasets largely due to an a priori assumption that bayesian methods cannot cope with datasets with > terminals e.g.,. our study with species and >  bp of sequence data demonstrates that large matrices with > terminals can be analysed using bayesian dating methods. further studies are needed, however, to fully explore best methods for analysing large datasets with the currently available dating methods that implement relaxed molecular clock models required for analyses of diverse clades where rates are expected to vary. such studies should focus on exploring trade-offs between number of taxa, complexity of models and partitions used in order to fully understand limitations and potential error sources in large scale analyses. in our dating analysis, we followed the recent recommended best practice guidelines for fossil calibration and placed fossil calibrations at stem nodes of the most inclusive extant groups using apomorphy-based morphological assignment. morphological evidence from the seed fossils only allowed assignment to the broad groups solanoideae or solanaceae as a whole. fossils provide only minimum age estimates for the nodes they are assigned, and hence results from our dating analysis where fossil calibrations were used should be considered as minimum age estimates. we further biased our results towards younger ages by assigning the oldest known fossils of solanaceae to the stem node of the family rather than to more specific nodes within solanaceae due to lack of morphological and anatomical characters that could be used to assign them to more specific nodes. there is always a possibility, however, that these seeds represent more specific clades within solanaceae, which would push back age estimates for the family. currently, the earliest fossil evidence for the family comes from eocene europe, but based on biogeographic analyses, the crown group of solanaceae is thought to have originated and first diversified in south america. total evidence analysis, where fossils are placed as terminal taxa in the dating analysis using both molecular and morphological data matrix, could help in exploring the robustness of fossil placement, but as pointed above, the lack of characters in the solanaceae seed fossils does not currently permit such analyses. the most promising avenue in strengthening the dating analysis would be in finding further fossil records. this would increase the number of fossil calibration points and allow the use of cross-validation methods. the rate of molecular evolution in plants has been found to correlate with life history traits, whereby longer living species show consistently lower substitution rates compared to shorter living species. molecular clock models should incorporate such rate variation, especially in groups such as solanaceae which include a range of growth and life forms. our dating analyses did not incorporate such models, although the model used in our bayesian analysis allows rates to vary between lineages independently. the lack of such models in our analyses implies that the age of herbaceous, shorter lived plants will be systematically overestimated, while ages in dominantly woody clades will be consistently underestimated. future studies should explore how molecular clock models that account for rate variation due to life history traits could be implemented. previous studies have produced a wide range of estimates for the stem age of the family, ranging from 34–85 ma, but none of these studies included dense sampling within the family nor used robust solanaceae-specific fossil calibrations. paape et al. analysed divergence times within solanaceae but with a small dataset consisting of species only. this study was based on three fossil calibration points without re-assessment or morphological study of the original fossils, and estimated solanaceae stem age to have diverged 62 ma. the oldest estimates for the family stem node age come from earlier molecular studies which used calibration points with more simplistic dating methods, while the most recent molecular dating study of angiosperms by bell et al. who used fossil calibrations across the tree and a relaxed molecular clock model, estimated the solanaceae stem node to have diverged c. 59 ma. our results, which we consider as minimum ages, are broadly consistent with bell et al. in estimating the stem node of solanaceae to date back to c. 49 ma. the age of the major splits within the family has been of interest to various fields, including studies on chromosomal and genome evolution. our minimum age estimates for the major splits between tomato – potato, eggplant – tomato/potato, solanum – capsicum, and solanum – nicotiana are consistent with the age estimates produced in previous studies without fossil calibrations using much sparser sampling and more simplistic molecular clock models. our results for the nicotiana – symonanthus split corroborate results obtained using island age and those calculated using paralogy-free subtree analysis. our results presented here suggest that the rate of chromosomal and genome evolution within solanaceae has been marginally slower at least within particular lineages than previously thought. with the densely sampled chronogram presented in this study, a more detailed analysis of chromosomal evolution at the species level could now be performed in the solanaceae to study rate differences and drivers of chromosomal changes such as environmental or life history factors. similarly, morphological characters such as fruit type could be analysed in relation to diversification rates to identify whether particular morphological traits are associated with speciation rate shifts in solanaceae.  despite much focus on character and trait evolution within solanaceae, little has been known about the origin of traits in the family in terms of time. we present here minimum age estimates and associated confidence intervals for the entire solanaceae using a species-rich dataset comprising almost half of the species diversity within the family. this densely sampled chronogram will provide the basis for unravelling the tempo and mode of evolution of many of the much-studied and complex traits in this diverse and economically important family such as self-incompatibility, fruit type, cold and salt tolerance, disease resistance, chromosomal re-arrangements, genome size, and gene sub-functionalization. 
the minimum age estimates for major splits within the family provided here correspond well with results from previous studies, indicating splits between tomato & potato around million years ago with a 95% highest posterior density 7–10 ma, solanum & capsicum c. 19 ma, and solanum & nicotiana c. 24 ma. we also review the fossil record of the family to derive robust calibration points, and estimate a chronogram using an uncorrelated relaxed molecular clock. our densely-sampled phylogeny shows strong support for all previously identified clades of solanaceae and strongly supported relationships between the major clades, particularly within solanum. the tomato clade is shown to be sister to section petota, and the regmandra clade is the first branching member of the potato clade. despite a wealth of phylogenetic work on individual clades and a deep knowledge of particular cultivated species such as tomato and potato, a robust evolutionary framework with a dated molecular phylogeny for the family is still lacking. our large time-calibrated phylogeny provides a significant step towards completing a fully sampled species-level phylogeny for solanaceae, and provides age estimates for the whole family. here we investigate molecular divergence times for solanaceae using a densely-sampled species-level phylogeny.
 the solanaceae is a plant family of great economic importance. despite a wealth of phylogenetic work on individual clades and a deep knowledge of particular cultivated species such as tomato and potato, a robust evolutionary framework with a dated molecular phylogeny for the family is still lacking. here we investigate molecular divergence times for solanaceae using a densely-sampled species-level phylogeny. we also review the fossil record of the family to derive robust calibration points, and estimate a chronogram using an uncorrelated relaxed molecular clock.  our densely-sampled phylogeny shows strong support for all previously identified clades of solanaceae and strongly supported relationships between the major clades, particularly within solanum. the tomato clade is shown to be sister to section petota, and the regmandra clade is the first branching member of the potato clade. the minimum age estimates for major splits within the family provided here correspond well with results from previous studies, indicating splits between tomato & potato around million years ago with a 95% highest posterior density 7–10 ma, solanum & capsicum c. 19 ma, and solanum & nicotiana c. 24 ma.  our large time-calibrated phylogeny provides a significant step towards completing a fully sampled species-level phylogeny for solanaceae, and provides age estimates for the whole family. the chronogram now includes 40% of known species and all but two monotypic genera, and is one of the best sampled angiosperm family phylogenies both in terms of taxon sampling and resolution published thus far. the increased resolution in the chronogram combined with the large increase in species sampling will provide much needed data for the examination of many biological questions using solanaceae as a model system. 