Skip Navigation

Systematic Biology 2005 54(4):575-594; doi:10.1080/10635150590947131
This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (68)
Right arrowRequest Permissions
Google Scholar
Right arrow Articles by Lavin, M.
Right arrow Articles by Wojciechowski, M. F.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Lavin, M.
Right arrow Articles by Wojciechowski, M. F.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

© 2005 Society of Systematic Biologists

Evolutionary Rates Analysis of Leguminosae Implicates a Rapid Diversification of Lineages during the Tertiary

Edited by Peter Linder: Associate Editor

Matt Lavin1, Patrick S. Herendeen2 and Martin F. Wojciechowski3

1 Department of Plant Sciences, Montana State University Bozeman, Montana, 59717, USA E-mail: mlavin{at}montana.edu
2 Department of Biological Sciences, The George Washington University 2023 G Street NW, Washington, DC, 20052, USA
3 School of Life Sciences, Arizona State University Tempe, Arizona, 85287–4501, USA


    Abstract
 Top
 Abstract
 Materials and Methods
 Results
 Discussion
 Appendix 1. The 13...
 References
 
Tertiary macrofossils of the flowering plant family Leguminosae (legumes) were used as time constraints to estimate ages of the earliest branching clades identified in separate plastid matK and rbcL gene phylogenies. Penalized likelihood rate smoothing was performed on sets of Bayesian likelihood trees generated with the AIC-selected GTR+ {Gamma} +I substitution model. Unequivocal legume fossils dating from the Recent continuously back to about 56 million years ago were used to fix the family stem clade at 60 million years (Ma), and at 1-Ma intervals back to 70 Ma. Specific fossils that showed distinctive combinations of apomorphic traits were used to constrain the minimum age of 12 specific internal nodes. These constraints were placed on stem rather than respective crown clades in order to bias for younger age estimates. Regardless, the mean age of the legume crown clade differs by only 1.0 to 2.5 Ma from the fixed age of the legume stem clade. Additionally, the oldest caesalpinioid, mimosoid, and papilionoid crown clades show approximately the same age range of 39 to 59 Ma. These findings all point to a rapid family-wide diversification, and predict few if any legume fossils prior to the Cenozoic. The range of the matK substitution rate, 2.1–24.6 x 10–10 substitutions per site per year, is higher than that of rbcL, 1.6–8.6 x 10–10, and is accompanied by more uniform rate variation among codon positions. The matK and rbcL substitution rates are highly correlated across the legume family. For example, both loci have the slowest substitution rates among the mimosoids and the fastest rates among the millettioid legumes. This explains why groups such as the millettioids are amenable to species-level phylogenetic analysis with these loci, whereas other legume groups are not.

Keywords: Age estimation; Bayesian phylogenetics; Fabaceae; Leguminosae; matK; penalized likelihood rate smoothing; rbcL; substitution rates

Received June 3, 2004; Revised August 31, 2004; Accepted November 18, 2004


The flowering plant family Leguminosae (Fabaceae) contains over 18,000 species distributed throughout the world in many ecological settings, from deserts of high latitudes to seasonally dry or wet tropical forests of equatorial regions (Lewis et al., 2005). Legumes diversified during the Early Tertiary (Herendeen et al., 1992) to become a ubiquitous feature of modern terrestrial biotas, similar to the timing of diversification of other prevailing terrestrial groups, such as the dominant modern families of angiosperms, polypod ferns, mammals, teleost fishes, birds, and insects (e.g., Wilf and Labandeira, 1999; Engel, 2001; Manos and Standford, 2001; Schneider et al., 2004). Second only to the grass family (Poaceae) in agricultural and economic importance, the legumes includes many species harvested as crops, or used for oils, fiber, fuel, timber, medicinals, chemicals, and horticultural varieties. Legumes play an important role in the terrestrial nitrogen cycle regardless of whether they form root nodules with symbiotic rhizobia (Sprent, 2001).

The first definitive legumes appear during the Late Paleocene, about 56 million years ago (Mya) (Herendeen, 2001; Herendeen and Wing, 2001; Wing et al., 2004). All three traditionally recognized subfamilies of legumes, the caesalpinioids, mimosoids, and papilionoids (Polhill et al., 1981), as well as other taxonomically large clades within these subfamilies (e.g., genistoids), are recorded from the fossil record soon afterward, beginning around 50 to 55 Mya (e.g., Herendeen et al., 1992). A prediction derived from the legume fossil record is that there should be little difference between the estimated age of the origin of legumes and their subsequent diversification. Now that legumes and close relatives are well sampled for molecular data (e.g., Kajita et al., 2001; Herendeen et al., 2003a; Wojciechowski et al., 2004), the goal of this article is to sample diverse legumes and immediate outgroups such that the age of the legume stem (origin) and crown (diversification) clades can be confidently estimated.

A hypothesis that posits a rapid diversification of legumes also predicts that the oldest caesalpinioid, mimosoid, and papilionoid crown clades should be equivalent in age and nearly as old as the legume crown. This finding would counter conventional wisdom that implicitly suggests that caesalpinioids harbor many ancestral traits and comprise the oldest lineages (e.g., Polhill et al., 1981). The "basal" characterization of caesalpinioids (e.g., Bruneau et al., 2001) is possibly one reason why Tucker and Douglas (1994) choose outgroups that shared radial floral symmetry with caesalpinioids (e.g., Connaraceae, Cunoniaceae, and Sapindaceae). The term "basal," however, is vague at best (e.g., Krell and Cranston, 2004). Although caesalpinioids form a paraphyletic grade from which stem the mimosoids and papilionoids (Kajita et al., 2001; Wojciechowski et al., 2004), attributing "basal" to basally branching clades would mean nothing regarding ages of extant diversifications.

With the exception of the few cursory reports in Wojciechowski (2003), no molecular age estimates have yet been detailed for crown groups within the legume family. This rate and age analysis is not only comprehensive, but incorporates multiple fossil constraints and ages estimated from multiple trees, approaches that increase the reliability of rate and age estimates (e.g., Yang and Yoder, 2003; Yoder and Yang, 2004). The rate and age analysis detailed in this study therefore can reveal not only the tempo of legume diversification, but also provide credible age estimates for crown clades with no fossil record.


    Materials and Methods
 Top
 Abstract
 Materials and Methods
 Results
 Discussion
 Appendix 1. The 13...
 References
 
Taxon Sampling
Only two studies have extensively sampled molecular sequence data from across legumes and close relatives (Kajita et al., 2001; Wojciechowski et al., 2004). The sequences sampled from these two studies for this analysis represent 38 of the 40 legume tribes recognized by Polhill (1994), all of the well-supported legume crown clades recently delineated by molecular phylogenetic studies (reviewed in Wojciechowski, 2003), and the designated outgroups representing the two other families of Fabales (sensu Angiosperm Phylogeny Group, 2003), Polygalaceae and Surianaceae, as well as the rosaceous genus Quillaja.

DNA Sequence Data
Sequence data from Wojciechowski et al. (2004) and Kajita et al. (2001) represent two plastid loci, matK and rbcL, commonly used in phylogenetic studies of plants. The matK data set contained 335 sequences with an aligned length of 1674 sites, 1430 of which were used in the rates analysis after gapped regions were omitted (the flanking non-coding portion of the trnK intron was not included). The rbcL data set included 241 sequences with an aligned length of 1399 sites, all of which were used in the rates analysis. The chloroplast trnL-F region (e.g., Pennington et al., 2001; Bruneau et al., 2001; Hu et al., 2002) was too invariable and estimates at one node broadly overlapped with those of adjacent nodes (Lavin, unpublished data). The nuclear ribosomal ITS/5.8S region (Baldwin et al., 1995) is unalignable for family-wide comparisons. DNA isolations, primer specifications, PCR amplification conditions, DNA sequencing, methods of data management, and voucher specimen information follow Wojciechowski et al. (2004) and Kajita et al. (2001). Data matrices are deposited with TreeBase (http://www.treebase.org/) study accession number S1968 for matK (Wojciechowski et al., 2004) and S578 for rbcL (Kajita et al., 2001).

Evolutionary Rates Analysis
Branch lengths were estimated using a Bayesian approach (Huelsenbeck and Ronquist, 2001; Huelsenbeck et al., 2001). The GTR+SS, GTR+ {Gamma} +I, and other nucleotide substitution models for both the matK and rbcL data sets were evaluated with Akaike information criterion (AIC) and Schwarz criterion (SC, or Bayesian information criterion), as implemented manually (e.g., Burnham and Anderson, 2002; Johnson and Omland, 2004) and in ModelTest (Posada and Crandall, 1998). Selected models were subjected to multiple runs of Metropolis-coupled Markov Chain Monte Carlo chains, each initiated with random starting trees and default substitution parameters.

Each Bayesian run comprised four chains of default temperatures and 5–20 x 106 generations. Tree parameters were sampled every 5 x 104 generations after likelihood stationarity was attained. One hundred Bayesian trees at stationarity were then systematically sampled (e.g., once every 50,000 trees) for the rates analysis in order to estimate the mean, variance, minimum, and maximum rate and age for each of 82 nodes representing the oldest crown clades within the legume family. Relative substitution rates and ages were converted to absolute rates and ages by enforcing age constraints derived from the fossil record.

A likelihood ratio test (Felsenstein, 1988; Huelsenbeck and Crandall, 1997; Posada and Crandall, 1998; Johnson and Omland, 2004) using either the GTR+ {Gamma} +I or GTR+SS nucleotide substitution models rejected a molecular clock for both the matK (LR = 2372, df = 333, P = 0.00000) and rbcL (LR = 1411, df = 239, P = 0.00000) data sets. The penalized likelihood method (Sanderson, 2002) in the program r8s (Sanderson, 2003) was thus used to estimate nucleotide substitution rates and ages of selected stem and crown clades. Via a cross-validation procedure, the penalized likelihood method finds an optimum rate-smoothing parameter for the transition of substitution rate between ancestor and descendants. The optimum level of smoothing should generally lie between the single-parameter (e.g., molecular clock; Langley and Fitch, 1974) and the parameter-rich model (e.g., typically nonparametric rate-smoothing; Sanderson, 1997).

Fossil Constraints
Age constraints derived from the legume fossil record were imposed on each of the matK and rbcL Bayesian consensus trees. The age of the root node was fixed and minimum ages were assigned to 12 internal nodes. Many fossils are potentially available, but our selection was limited to those with a unique set of apomorphic characters that could be assigned to a specific stem clade of legumes. Although either stem or crown clades may be constrained (e.g., Magallón and Sanderson, 2001), assigning fossil constraints to stem clades allowed more branch length (nucleotide substitutions) to be attributed to the minimum time frame defined by the fossil. This resulted in a faster rate of substitution and a bias toward younger age estimates. Because the objective of our study is to test the hypothesis that the oldest crown clades within legumes are indeed old relative to the legume stem clade, placing constraints on stem clades will not bias the results in favor of the hypothesis being tested. Of the 13 age constraints (A to M in Appendix 1) imposed on the matK phylogeny, only 9 (A, B, D, F, G, H, I, J, and K; see Appendix 1) could be placed on the rbcL phylogeny because of less exhaustive sampling. Absolute ages of relevant geological periods follow Berggren et al. (1996).


    Results
 Top
 Abstract
 Materials and Methods
 Results
 Discussion
 Appendix 1. The 13...
 References
 
Rate-Smoothing Parameters
The optimal rate-smoothing parameter selected during cross-validation of the penalized likelihood analysis was 101.5 for matK and 102 for rbcL. In both instances, these smoothing parameters resulted in a better goodness of fit (estimated terminal branch lengths versus the actual terminal branch lengths) than that determined for the globally rate constant (Langley-Fitch) analysis. Nonparametric rate-smoothing consistently yielded the worst fit.

Nucleotide Substitution Estimates for the matK Locus
Both AIC and SC revealed the best fit of the GTR+ {Gamma} +I model compared to all other nested and non-nested models (e.g., AIC of GTR+SS = 105608 and of GTR+ {Gamma} +I = 103094 for a {Delta} AIC of 2514; or SC = 105677 for GTR+SS and 103158 for GTR+ {Gamma} +I for a {Delta} SC of 2519). Parameters estimated during the Bayesian analysis using the GTR+ {Gamma} +I and GTR+SS models (Table 1) reveal, among other things, a fairly uniform among-site rate variation. For example, the substitution rate at the 3rd codon position is only about twice that of the 2nd position.


View this table:
[in this window]
[in a new window]

 
Table 1 Nucleotide substitution parameters for the GTR+ {Gamma} +I and GTR+SS models for each of the matK and rbcL loci. Parameters were estimated by Bayesian inference and sampling at likelihood stationarity (i.e., 100 trees). These models involve six substitution parameters (r), four base frequency parameters (pi), the alpha shape parameter, the proportion of invariant sites parameter [p(I)], and three relative substitution rates among codon positions, 1st (m{1}), 2nd (m{2}), and 3rd (m{3}). The GTR+ {Gamma} +I model for both loci had the best fit to the data, as revealed by either AIC or SC.

 
The mean penalized likelihood (PL) substitution rate across 100 Bayesian trees for the 82 nodes identified in the matK phylogeny (Figs. 1, 2, 3, Table 2) ranges from 2.1 to 24.6 x 10–10 substitutions per site per year. When the 12 minimum age constraints (nodes B to M, Appendix 1) were not imposed during the PL rates analysis, a modest trend toward a faster substitution rates was observed (Fig. 4a). Three points lie distinctly above the trend line (Fig. 4a) and include two constrained nodes, K (MRCA of Machaerium falciforme and Aeschynomene purpusii; Table 2, Fig. 2) and M (MRCA of Robinia pseudoacacia and Coursetia axillaris; Table 2, Fig. 3), and one adjacent to K (node 45; Table 2, Fig. 2).


Figure 1
View larger version (79K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 1 Bayesian consensus matK phylogeny (i.e., with averaged branch lengths) showing outgroups, caesalpinioid, mimosoid, and some of the papilionoid crown clades (e.g., the Swartzia node and the Cladrastis crown). Estimated substitution parameters for this Bayesian analysis are given in Table 1. Detailed age and rate estimates for all nodes labeled with letters or numbers are presented in Table 2. Age estimates are reported for the older crown clades. Ma = million years. The nodes labeled with diamonds are those with fixed (node A) or minimum age constraints (nodes B to M).

 


Figure 2
View larger version (72K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 2 Continuation of the Bayesian consensus matK phylogeny from Figure 1 showing the genistoid, dalbergioid, and mirbelioid crown clades. Estimated substitution parameters for this Bayesian analysis are given in Table 1. Detailed age and rate estimates for all nodes labeled with letters or numbers are presented in Table 2. Age estimates are reported for the older crown clades. Ma = million years. The nodes labeled with diamonds are those with fixed (node A) or minimum fossil constraints (nodes B to M).

 


Figure 3
View larger version (72K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 3 Continuation of the Bayesian consensus matK phylogeny from Figure 2 showing the indigoferoid, millettioid, and Hologalegina crown clades. Estimated substitution parameters for this Bayesian analysis are given in Table 1. Detailed age and rate estimates for all nodes labeled with letters or numbers are presented in Table 2. Age estimates are reported for the older crown clades. Ma = million years. The nodes labeled with diamonds are those with fixed (node A) or minimum fossil constraints (nodes B to M).

 


Figure 4
View larger version (39K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 4 Descriptive statistical comparisons of mean rate estimates (see Tables 2 and 3) derived from the matK and rbcL data sets. (a) Rate estimates derived from matK using penalized likelihood (PL) for the 94 nodes (B to M, 1 to 82) listed in Table 1 compared to PL estimates for the same nodes derived from not invoking the 12 fossil-calibrated minimum age constraints; subs/site/Ma = substitution per site per million years. (b) Rate estimates derived from matK using PL with the 12 minimum age constraints imposed and compared to the rate constant Langley-Fitch (LF) and rate variable nonparametric rate smoothing (NPRS) methods. The optimum substitution rate estimated with LF using minimal age constraints is indicated by the line with an intercept of 0.00081 substitutions/site/Ma (diamonds). Without minimum age constraints, the rate constant estimate is 0.00144 substitutions/site/Ma. (c) Comparison of rate estimates derived from each of the matK and rbcL Bayesian consensus phylogenies for the 28 comparable nodes listed in Table 3 (not including root node A) using PL rate smoothing (see Figs. 1 and 3 for the identification of the labeled nodes).

 


View this table:
[in this window]
[in a new window]

 
Table 2 Penalized likelihood estimates of ages for fossil-constrained and unconstrained nodes across 100 matK Bayesian trees sampled at likelihood stationarity (for node assignment, see the Bayesian consensus, Figs. 1, 2, 3). The first node (A) is assigned a fixed age, whereas the next 12 nodes (B to M) are assigned minimum age constraints (indicated in parentheses) as derived from the fossil record. The letter and numeric codes for each node correspond to those shown in Figures 1, 2, 3. Ma = million years; SD = standard deviation; s/s = substitutions per site.

 
An average of 8.1 x 10–10 substitutions per site per year was estimated with the rate constant Langley-Fitch (LF) model (Fig. 4b), which is well within the range of estimates derived from the penalized likelihood method (Fig. 4a). This rate was increased to 14.4 x 10–10 substitutions per site per year when the 12 minimum fossil constraints were removed (data not shown). The nonparametric rate-smoothing (NPRS) approach resulted in highly variable rate estimates that showed no correlation with rates estimated using the penalized likelihood approach (Fig. 4b).

Nucleotide Substitution Estimates for the rbcL Locus
Both AIC and SC revealed the best fit of the GTR+ {Gamma} +I model compared to all other nested and non-nested models (e.g., AIC of GTR+SS = 49008 and of GTR+ {Gamma} +I = 44614 for a {Delta} AIC of 4395; or SC = 49077 for GTR+SS and 44677 for GTR+ {Gamma} +I for a {Delta} SC of 4400). Parameters estimated during the Bayesian analysis using the GTR+ {Gamma} +I and GTR+SS models (Table 1) reveal that substitution variation at the 3rd codon position was almost 10-fold that of the 2nd position.

The mean PL substitution rate across 100 Bayesian trees for the 29 nodes identified in the rbcL phylogeny ranges from 1.6 to 8.6 x 10–10 substitutions per site per year (Table 3). If minimum fossil constraints are not imposed, a modest trend toward faster substitution rates is detected. Only two nodes show a significant deviation from the general trend (graph not shown) and include node K (MRCA of Machaerium lunatum and Dalbergia hupeana; Table 3, Fig. 2), and node G (MRCA of Acacia farnesiana and Mimosa speggazzinii; Table 3, Fig. 1). The rate constant LF model showed the same trend in estimating faster substitution rates when the minimum fossil constraints were not imposed (means of 4.8 versus 7.0 x 10–10 substitutions per site per year).


View this table:
[in this window]
[in a new window]

 
Table 3 Penalized likelihood estimates of ages for fossil-constrained and unconstrained nodes across 100 rbcL Bayesian trees sampled at stationarity (not shown but see parsimony consensus in Kajita et al., 2001). The legume stem clade (A) is fixed at 60 Ma, and the next eight nodes listed (B, D, F, G, H, J, K, and M) are assigned minimum age constraints (indicated in parentheses) as derived from the fossil record. These lettered and numbered nodes in the rbcL phylogeny correspond to a node in the matK phylogeny and are labeled to identify this association (see Figs. 1, 2, 3). Ma = million years; SD = standard deviation; s/s = substitutions per site.

 
Comparisons of Substitution Rates Estimated for the matK and rbcL Loci
Although penalized likelihood estimates derived from the matK locus have a substitution rate that is on average about twice that of rbcL, some lineages show a similar rate between these two loci (Fig. 4c). The mimosoid stem clade (node F; Fig. 1) and the Cladrastis crown clade (node 17; Fig. 1) distinctly show the slowest mean rates for both loci (Fig. 4c). For example, the Cladrastis crown is estimated at 2.6 x 10–10 substitution per site per year for matK and 1.6 x 10–10 for rbcL (node 17, Figs. 1, 4b; Tables 2, 3). The estimated substitution rate for the inverted-repeat-loss clade (IRLC; node 75, Figs. 3, 4c) is 7.8 x 10–10 substitutions per site per year for the matK locus and 7.2 x 10–10 for rbcL. In general, the matK locus has a much faster rate of evolution, with the most extreme rates observed in lineages traditionally classified in the papilionoid tribe Phaseoleae (nodes 63 to 68 in Table 2; Fig. 3) with mean rates of about 20 x 10–10 substitutions per site per year or greater. The substitution rate estimated for node 65 in particular (the MRCA of Pueraria montana and Phaseolus vulgaris in Fig. 2; Tables 2, 3) is distinctly the highest for both matK and rbcL (Fig. 4c).

Estimated Ages of Old Legume Crown Clades
The ages estimated for the 28 nodes that are comparable between the matK and rbcL phylogenies show a positive linear relationship (Tables 2, 3; Fig. 5a). In both phylogenies, the estimated age of the legume crown is 59 million years (Ma) (node 1 in Tables 2, 3). The estimated ages derived from the penalized likelihood rates smoothing for older legume crown clades are mostly less than 30 Ma (Figs. 1, 2, 3, 6, 7). These age estimates show a strong positive relationship with the rate constant LF and the rate variable NPRS estimates, although NPRS tends to yield the oldest ages of the three methods (Fig. 5b). Imposing the 12 minimum age constraints on nodes B to M (Figs. 1, 2, 3) generally resulted in age estimates that are older than estimates without the fossil constraints (Fig. 5c).


Figure 5
View larger version (41K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 5 Descriptive statistical comparisons of mean age estimates derived from the matK and rbcL data sets (see Tables 2 and 3). (a) Comparison of age estimates derived from the matK and rbcL Bayesian phylogenies for the 28 comparable nodes listed in Table 3 (not including root node A) using the penalized likelihood (PL) rate smoothing approach. Ma = million years. (b) Age estimates derived from matK imposing minimal age constraints. PL estimates for the 94 nodes listed in Table 1 (not including root node A) are compared with those derived from rate constant Langley-Fitch (LF; diamonds) and rate variable nonparametric rate smoothing (NPRS; squares) methods. (c) Age estimates derived from matK with and without minimal age constraints. PL estimates using minimum age constraints for the 94 nodes listed in Table 1 (not including root node A) are compared with those derived from PL (diamonds) and LF (squares) methods without minimum age constraints.

 


Figure 6
View larger version (75K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 6 Phylogenies derived from the matK data set. (a) Bayesian consensus phylogram. Scale bar equals 0.005 substitutions per site. (b) Penalized likelihood rate-smoothed chronogram with the 12 minimum age constraints imposed. Scale bar equals 10 Ma. The vertical lines (b) divide the 60-Ma duration of legume evolution into 15-Ma segments. The oldest and taxonomically large crown clades are alternately shaded black and gray. The shapes of both phylogenies depict generally short internal branch lengths indicative of a rapid family-wide diversification.

 


Figure 7
View larger version (67K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 7 Phylogenies derived from the rbcL data set. (a) Bayesian consensus phylogram. Scale bar equals 0.005 substitutions per site. (b) Penalized likelihood rate-smoothed chronogram with the 12 minimum time constraints enforced. Scale bar equals 10 Ma. The vertical lines (b) divide the 60-Ma duration of legume evolution into 15-Ma segments. The shapes of both phylogenies depict generally short internal branch lengths indicative of a rapid family-wide diversification. The oldest and taxonomically large crown clades are alternately shaded black and gray, with the top gray shaded lines representing the outgroups and the following paraphyletic grade of black lines representing the caesalpinioid legumes. Otherwise, the groups are labeled. The arrows indicate the legume stem clade.

 
The ages of the oldest caesalpinioid crown clades range from approximately 54 to 30 Ma (nodes B to E, Table 2; and nodes 2 to 9, Fig. 1). The oldest of these, nodes 5 and D (the Umtiza crown node), have estimated mean ages of about 54 Ma (Fig. 1). In contrast, the extant mimosoids trace back to a most recent common ancestor at just over 40 Ma (node 10, Fig. 1). The older papilionoid crown clades show ages comparable to the oldest extant caesalpinioid diversifications. For example, the age of the genistoid diversification is estimated at about 56 Ma (node 21, Fig. 2), that of the dalbergioids at about 55 Ma (node 34, Fig. 2), and Hologalegina at about 50 Ma (node 70, Fig. 3). Indeed, nearly all of the remaining oldest papilionoid crown clades (e.g., nodes 15, 17 in Fig. 1; node 55 in Fig. 2; node 58 in Fig. 3) are around 45 Ma or older. The ages of extant papilionoid diversifications relative to those of caesalpinioids do not change when minimum age constraints are not imposed. In this case, all of the estimated ages are younger by about 1 to 5 Ma (data not shown), the difference being greater with caesalpinioid clades because this is where most of the minimum fossil constraints are imposed (see positions of constraints B to F in Fig. 1).

Effects of Varying the Age of the Root Node
The analyses presented above are derived from fixing the age of the legume stem clade at 60 Ma. Moving this age at 1-Ma intervals back to 70 Ma has little qualitative effect on the results. At the 60-Ma fixed age, the estimated ages of the legume and papilionoid crown clades lay distinctly above the trend line to the upper right (Fig. 8a and b). These two age estimates, when plotted at 1-Ma intervals (Fig. 8c; the legume crown is illustrated for both matK and rbcL), are exceptional in showing successively older ages as the fixed age of the root node is moved back to 70 Ma. All other nodes for which ages are estimated have a flat trend line over the 60- to 70-Ma interval, as exemplified by those of the caesalpinioid and mimosoid crowns (Fig. 8). Regardless, the difference in mean age estimates between the legume stem and crown clades vary between 1.0 and 2.5 Ma within this 60- to 70-Ma interval.


Figure 8
View larger version (42K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 8 Descriptive statistical comparisons of mean age estimates when the fixed age of the legume stem clade is scaled between 60 and 70 Ma. (a) Age estimates derived from the matK data set for the 94 nodes identified in Table 2 (B to M, 1 to 82; not including the root node). (b) Age estimates derived from the rbcL data set for the 28 nodes identified in Table 3 that are shared with the matK phylogeny (not including root node). (c) Ages estimates derived from the matK data set or rbcL data set (where noted) for selected old crown clades.

 

    Discussion
 Top
 Abstract
 Materials and Methods
 Results
 Discussion
 Appendix 1. The 13...
 References
 
The matK tree is more resolved and better supported than the rbcL tree, a finding consistent with the parsimony analyses of essentially these same data sets (cf. Wojciechowski et al., 2004; Kajita et al., 2001). The present study suggests this may be the result of the rate of substitution at the matK locus, which is more evenly distributed among the three codon positions and generally about twice as fast as that of rbcL (matK range of 2.1 to 24.6 x 10–10 versus rbcL range of 1.6 to 8.6 x 10–10). These estimated rates of substitution in matK and rbcL, the rate differences between the two loci, and the among-site rate variation, are very similar to those reported for other plant groups (e.g., Albert et al., 1994; Steele and Vilgalys, 1994; Wang et al., 1999; Osaloo and Kawano, 1999; Young and dePamphilis, 2000; Xiang et al., 2000; Zhang and Renner, 2003; Neel and Cummings, 2004).

The rapid diversification of the legume family soon after its origin is revealed by both the matK and rbcL phylogenies (Figs. 6, 7). An estimated 1.0 to 2.5 Ma distinguishes the mean ages of the legume stem and crown clade. This may explain why clade support values for the monophyly of Leguminosae are sometimes lower compared to support levels for the entire Fabales crown (i.e., the MRCA of Quillaja, Surianaceae, Polygalaceae, and Leguminosae; sensu Angiosperm Phylogeny Group, 2003) or many of the older crown clades within legumes (Bruneau et al., 2001; Kajita et al., 2001; Wojciechowski, 2003; Wojciechowski et al., 2004).

The total duration of legume evolution was fixed at 60 Ma in our analyses, just older than the oldest definitive legume fossil at 56 Ma (Herendeen, 2001). Several of the earliest minimal time constraints are in the 55- to 56-Ma range (nodes F, H, and J, Table 2; Figs. 1, 2). Regardless, the legume crown could have been estimated closer to a 56-Ma age instead of about 59 Ma (node 1, Fig. 1). If all minimum fossil constraints are removed, the age of the legume crown is estimated at 55.9 ± 1.0 Ma, also a short time span relative to the total duration of legume evolution. Notably, increasing the fixed age of the legume stem to 70 Ma still resolves a difference on average of up to 2.5 Ma between the estimated age of the legume stem and crown (Fig. 8). These findings support the overall conclusion that an Early Tertiary legume diversification immediately followed the origin of the family.

The fossil record of legumes predicts this genetic finding of a rapid diversification of extant lineages. Many diverse legumes, especially representatives from each of the three traditionally recognized subfamilies, have continuous fossil records from the Recent back to the Late Paleocene (e.g., Herendeen and Dilcher, 1992). Moreover, this record comes from North and South America, Europe, Africa, and Asia (e.g., Axelrod, 1992; Herendeen et al., 1992). The caesalpinioids, with an estimated 161 genera and about 3000 species (Lewis et al., 2005), are most diverse in tropical regions throughout the New World, Africa, and southeast Asia. Although caesalpinioids form a basal grade within which are nested the other two subfamilies (e.g., Polhill, 1981, 1994; Herendeen et al., 2003a, 2003b), the estimated ages of the older caesalpinioid crown clades differ little from those of the mimosoids, and especially the papilionoids (compare marked crown clades in Figs. 1, 2, 3).

The papilionoid crown clade, and not a caesalpinioid one, is most sensitive to the effects of moving the fixed root age from 60 back to 70 Ma (Fig. 8). This implies that the oldest legume fossils are as likely to represent papilionoids as caesalpinioids. This is actually predicted by the fossil record. The oldest known legume fossils, with an age of 56 Ma (Herendeen, 2001), are fruits similar to Diplotropis and Bowdichia, representatives of the genistoid clade of papilionoids (Crisp et al., 2000). Fossil flowers assigned to Barnebyanthus buchananensis are close in age to these oldest legume fossils and in overall morphology to extant genistoid genera such as Bowringia and Clathrotropis (Crepet and Herendeen, 1992). Fossil leaves and fruits with affinities to the genistoid Ormosia come from strata at least 56 Ma (Wing et al., 2004, and unpublished).

The dalbergioids (Lavin et al., 2001), Hologalegina (Wojciechowski et al., 2000), and mirbelioids (Crisp et al., 2000), in contrast to the genistoids, lack an Early Eocene fossil record. Although these three large crown clades include papilionoids traditionally considered derived (Polhill et al., 1981; Polhill, 1994), they all have age estimates in the 50-Ma time frame or older (Figs. 2, 3). Indeed, Hologalegina contains many of the well-known temperate and herbaceous species of legumes, such as alfalfa, clovers, peas, and vetches. Remarkably, the mean age of this crown diversification is estimated at 50.6 Ma (Fig. 3, Table 2).

Subfamily Mimosoideae, with an estimated 80 genera and some 3000 species (Lewis et al., 2005), are mostly tropical to subtropical in distribution, and are abundant in arid and semiarid regions throughout the world. Of all the old legume crown clades, the mimosoids have a paucity of nucleotide substitution variation for both matK and rbcL (nodes F and 10 in Tables 2, 3; Fig. 4c). The apparent slow down in the rate of nucleotide substitution in this group may explain the findings in recent molecular systematic studies (Luckow et al., 2003; Miller et al., 2003) of many short and poorly supported branch lengths within the mimosoid crown. However, the mimosoids have a relatively long stem branch leading to the crown and accordingly a significant difference between the age of the origin and diversification of extant members the subfamily (about 15 Ma; compare nodes F and 10 in Tables 2, 3). These estimated ages derived independently from matK and rbcL are highly congruent and suggest that the lack of nucleotide variation in mimosoids is a function of both a rate slow-down and a relatively recent extant diversification, beginning about 40 Ma (see also Figs. 7b, 8b).

Numerous mimosoid leaf fossils come from the Early Tertiary, such as the 46 million year old Acacia mahenge fossil (Appendix 1). However, none can be assigned unequivocally to a modern subgroup. Given that no fossil constraints older than that of a 15-Ma-old fossil Acacia flower can be confidently assigned to an internal node within the mimosoid crown (node G, Tables 2, 3; Fig. 1), the combined molecular and fossil evidence suggests not only a relatively recent diversification of mimosoids, but also that many Early Tertiary mimosoid fossils belong to extinct lineages.

In contrast to the mimosoids, the tribe Phaseoleae (stemming from node 63 of the millettioid clade, Fig. 3) shows a very high rate of nucleotide substitution for both the matK and rbcL loci (Figs. 3, 4c, 5b; nodes 63 to 68 in Tables 2, 3). These two plastid loci contribute high levels of phylogenetic resolution and clade support in studies that sample at the genus and even species level within these clades (e.g., Kajita et al., 2001; Delgado-Salinas et al., 2003; Riley-Hulting et al., 2004; Thulin et al., 2004). That we detect much greater rates of nucleotide substitution in the Phaseoleae compared to the mimosoid crown is not an artifact of denser sampling in the former, as might be inferred from the findings of Webster et al. (2003, 2004). Both clades are approximately the same size in terms of actual numbers (i.e., approximately 3270 species of mimosoids and 2060 species of Phaseoleae; Lewis et al., 2005) and sampled numbers of species (i.e., 35 mimosoids and 25 Phaseoleae; Figs. 1, 3).

This study provides estimated ages that can serve as calibration points for rates analyses of legume clades lacking a fossil record (e.g., Schrire et al., 2003; Thulin et al., 2004; Lavin et al., 2004; Percy et al., 2004). Also, we have examined rate variation at the matK and rbcL loci across all of legumes, and have thus revealed that some groups such as the millettioids have a sufficiently fast rate of substitution so that these two plastid loci can be used in phylogenetic analyses to successfully resolve relationships even at the species level (e.g., Delgado-Salinas et al., 2003; Riley-Hulting et al., 2004; Thulin et al., 2004). In addition, the profusion of basally branching and morphologically diverse crown clades within legumes reveals that the paraphyletic grade of caesalpinioid lineages harbor neither the oldest diversifications nor some other quality of legume antiquity. Finally, the integration of fossil constraints with molecular phylogenetic data allows the identification of the rapid diversification early in the history of legume evolution. This new perspective would not have been possible with fossil data alone because as it turns out many of the oldest crown clades in legumes are unknown from the fossil record.


    Appendix 1. The 13 Fossil Constraints Imposed During This Study
 Top
 Abstract
 Materials and Methods
 Results
 Discussion
 Appendix 1. The 13...
 References
 
The legume family represents an ideal group for conducting an analysis of evolutionary rates and ages because the family has a high detectability in the fossil record. A rapid Paleocene appearance means that the duration of legume evolution probably includes just the last 60 to 65 Ma. The high detectability of legumes in the fossil record is postulated by the following evidence: (1) the characteristic nitrogen-rich metabolism predisposes legumes to high biomass production (McKey, 1994); (2) legumes are colonizers of disturbed sites close to depositional environments; (3) woody legume plants often produce abundant deciduous leaves, leaflets, flowers, or fruits; (4) these deciduous parts are often coriaceous or woody; and (5) legumes are taxonomically easily recognized in the fossil record (e.g., the cross-hatchings evident on the pulvinules of leaflets or pulvinus of leaves, compound leaves, leaflets often with a smooth margin and brochidodromous venation, and pod with single placental suture bearing a characteristic alternating branching pattern of the funiculi). The relationship of this high detectability to the continuous Cenozoic fossil record that abruptly begins by the Middle to Late Paleocene is being evaluated in a separate study with mark-recapture models (Rotella and Lavin, in preparation).

  1. The Leguminosae stem clade was fixed initially at 60 Ma, but also at every one million year interval back to 70 Ma. Only this age was fixed because a fixed time span facilitates rates and age estimation in all relaxed clock methods (Sanderson, personal communication). Furthermore, the temporally continuous and spatially extensive distribution of diverse fossil legumes from the Recent back to the Late Paleocene supports a hypothesized Early to Middle Paleocene origin of the family. The legume stem clade is defined as the most recent common ancestor (MRCA) of Polygala californica and Cercis occidentalis in the matK phylogeny (Fig. 1) and Polygala cruciata and Glycine tabacina in the rbcL phylogeny (not shown, but see Kajita et al., 2001). The most important pattern derived from the study of fossil legumes during the past several decades is that all three subfamilies have an abundant fossil record in North America, Europe, and Asia, and to some degree in South America and Africa (due to relatively poor sampling). This abundant record is represented continuously by at least fossil fruits and leaves from the Recent back to the Late Paleocene (e.g., Brown, 1962; Herendeen and Dilcher, 1992; Herendeen, 2001; Jacobs, 2003; and Wing et al., 2004). Earlier reports of legume fossils from as early as the Maastrichtian are tenuous and include only sporadic pollen and wood specimens that lack definitive legume apomorphies. Such legume apomorphies would include a prolate pollen grain with narrow colpi and large operculi, or wood with vestured pits (summarized in Herendeen et al., 1992).
    Wikström et al. (2001) posited a 74 to 79-Ma age (Late Cretaceous) for the legumes, and generally old estimates for other angiosperm families. These authors used only nonparametric rate smoothing, and constrained just a single node in the angiosperm phylogeny, the split between the Cucurbitales and the Fagales, which are two of three potential sister groups to Fabales (sensu Angiosperm Phylogeny Group, 2003). Implementation of a single age constraint combined with nonparametric rate-smoothing is bound to yield estimates highly biased toward older ages (Sanderson, 2002). Magallón and Sanderson (2001) used 59.9 Ma for the Fabales crown clade, which is equivalent to the legume stem. To accommodate uncertainty, we fixed the age of the legume stem at 60 Ma, but then compared results using this fixed age with those derived from fixed root ages at every one million year interval between 61 and 70 Ma.
  2. The Cercis stem clade was set to a minimum of 34 Ma (Late Eocene). This node is defined as the MRCA of Cercis occidentalis and Bauhinia tomentosa in the matK phylogeny (Fig. 1) and Cercis canadensis and Bauhinia purpurea in the rbcL phylogeny (Kajita et al., 2001). Fossils showing apomorphic traits of Cercis come from the Late Eocene of Sheep Rock Creek, Oregon (Herendeen and S. Manchester, unpublished data; site information provided in Wolfe and Tanai, 1987; McClain and Manchester, 2001). Fossils from this locality reveal a combination of apomorphies allowing assignment to Cercis, including unifoliolate leaves (i.e., with a lower pulvinus and upper pulvini), a thickish texture to the orbiculate-acuminate lamina, and a palmate primary venation combined with a pinnate secondary venation. Also found are thin-walled pods bearing a narrow non-vascularized placental wing suggestive of the fruits of Cercis. Additional leaf and fruit fossils also assigned to Cercis come from the Late Eocene Florissant Beds, Colorado (MacGinite, 1953; Manchester, 2001). Several types of legume fruits at Florissant have been identified as Cercis, however; and further study is required before these fossils can be definitively utilized here.
  3. The Hymenaea stem clade was set to a minimum of 34 Ma. This node is defined as the MRCA of Hymenaea courbaril and Tessmannia lescrauwaetii in the matK phylogeny (Fig. 1) and is undefined in the rbcL phylogeny. Poinar and Brown (2002) described Hymenaea mexicana from flowers preserved in amber with an estimated age of 22.5 to 26.0 Ma from Chiapas, Mexico. Hueber and Langenheim (1986) and Poinar (1991) place the fossil Hymenaea flower from La Toca Mines amber (Dominican Republic), formally named Hymenaea protera, as sister to the African Hymenaea verrucosa. Graham (1992) validated H. protera as bona fide Hymenaea with a minimum age of Late Eocene or 34 Ma. Poinar et al. (1993) further validated the generic assignment of H. protera with "fossil" rbcL sequence data. Only the Hymenaea crown clade, however, is well supported in their rbcL study. As it turns out, the Hymenaea rbcL sequences are small partial fragments, and contain certain sequence anomalies that need to be verified. In addition, proper outgroups were not sampled in Poinar et al. (1993), so the distinction of the Hymenaea stem and crown cannot be made. The above assignment, therefore, is the most precise possible given inadequate molecular sampling of Hymenaea and closely related genera.
  4. The Arcoa stem (or Umtiza crown; Herendeen et al., 2003a, 2003b) clade was set to a minimum of 34 Ma. This node is defined as the MRCA of Arcoa gonavensis and Gleditsia triacanthos in the matK phylogeny (Fig. 1) and Acrocarpus sp. (Genbank AF308699 [GenBank] ) and Gleditsia triacanthos in the rbcL phylogeny (Kajita et al., 2001). Fossil leaves from the Florissant Beds, Colorado (MacGinitie, 1953) referred to as "Prosopis linearifolia" reveal a mix of pinnate leaves and a bipinnate leaf where the pinnate leaves are larger than the bipinnate one. The leaflets are distinctly linear and asymmetric. Furthermore, in the single bipinnate leaf, the terminal group of three pinnae results from a sessile terminal pinna (e.g., specimen USNM40563). These leaf characteristic are very similar to those of Arcoa (now endemic to southern Hispaniola). The bipinnate leaf bearing a distinctively sessile terminal pinna is also found in Acrocarpus and Tetrapterocarpon, although leaflet morphology is considerably different in the latter two genera. These three genera plus Ceratonia, Gleditsia, Gymnocladus, and Umtiza constitute the Umtiza clade (Herendeen et al., 2003b). Notably, dehiscent pods winged along both sutures, very similar to those of Acrocarpus, are found in the Middle Eocene Bovay Clay pit, Tennessee (Herendeen, 1992). The combination of two wings and pod dehiscence, however, is not definitive for constraining the Umtiza crown to the earlier age of 40 Ma.
  5. The Mezoneuron stem clade was set to a minimum of 45 Ma. This node is defined as the MRCA of Caesalpinia (subgenus Mezoneuron) andamanica and Haematoxylum brasiletto in the matK phylogeny (Fig. 1) and is undefined in the rbcL phylogeny (Kajita et al., 2001). Several species of Mezoneuron (or Caesalpinia subg. Mezoneuron) are reported from North America and England as Middle Eocene to Miocene in age, even though Mezoneuron is presently confined to the Paleotropics (Herendeen and Zarucchi, 1990; Herendeen and Dilcher, 1991). The apomorphic traits that allow assignment of these fossils to the Mezoneuron stem include a membranous indehiscent pod with multiple ovules, and a broad placental wing with looping to longitudinal venation.
  6. The mimosoid stem clade was set to a minimum of 55 Ma. This node is defined as the MRCA of Inga punctata and Cercidium floridum in the matK phylogeny (Fig. 1) and Erythrophleum ivorense and Albizia saman in the rbcL phylogeny (Kajita et al., 2001). The fossil flower Protomimosoidea buchananensis (Crepet and Taylor, 1985, 1986) comes from the Wilcox Formation of western Tennessee, which has an age estimated at the Paleocene-Eocene boundary (55 Ma). The apomorphic characters that allow placement of this fossil to the mimosoid stem include a spicate inflorescence, actinomorphic bisexual flowers, valvate petals, tubular stigmas, and 10 exserted free stamens. Crepet and Taylor (1986) suggested an affinity to the tribe Mimoseae, but this was explicitly because "the fossil flowers lack the distinctive taxonomic characters that define the other tribes of mimosoid legumes...." In essence, the apomorphically diagnostic traits of these fossil flowers represent generalized mimosoid attributes and thus Protomimosoidea buchananensis constrains the mimosoid stem. Regardless, there are numerous other fossil leaves and fruits that collectively represent generalized mimosoid legumes from the earliest Eocene (e.g., Gros, 1992; Herendeen and Dilcher, 1990a; Herendeen and Jacobs, 2000). Notably, leaves, inflorescences, and fruits of a fossil species very similar to Dinizia are widespread from the southeastern USA during the Eocene (Herendeen and Dilcher, 1990a). Although Dinizia is no longer considered a mimosoid legume (Luckow et al., 2003), the Eocene fossils potentially referred to this genus could impose the same time constraints as the Protomimosoideae fossils. This is because the Dinizia stem clade is coeval with the mimosoid stem clade in the matK phylogeny (see Fig. 1).
  7. The Acacia stem clade was set to a minimum of 15 Ma (middle Miocene). This node is defined as the MRCA of Acacia hindsii and Inga punctata in the matK phylogeny (Fig. 1) and Acacia hindsii and Mimosa speggazzini in the rbcL phylogeny (Kajita et al., 2001). Dilcher et al. (1992) detailed the floral and leaflet morphology of specimens preserved in Dominican Republic amber from the Palo Alto mine. These mimosoid flowers come from sediments dated at the Oligocene-Miocene boundary. Dilcher et al. (1992) argue that the amber has been weathered out of older sediments and redeposited in this Oligocene-Miocene sediment. In contrast, Iturralde-Vinent and MacPhee (1996) presented evidence for the deposition of Dominican Amber during an interval 15 to 20 million years ago. Thus, the minimum age of the Acacia stem clade is 15 Ma. The assignment of the amber fossils to the Acacia stem (i.e., the first branching Acacia-containing lineage within the mimosoid crown) can be made because these fossils clearly show diagnostic apomorphic characters, such as numerous stamens with free filaments each with an anther bearing a stalked gland. Associated leaf pinnae bearing numerous small leaflets by themselves are not diagnostic, but they bear hairs identical to those on the flowers. The amber Acacia flowers are the only mimosoid fossils that unequivocally constrain a node nested within the mimosoid crown.
    Herendeen and Jacobs (2000) described Acacia mahengense from Middle Eocene (46 Ma) fossil leaves in Tanzania. This taxonomic assignment was derived from morphological similarity rather than definitive apomorphies. This is problematic from the perspective of this study given that Acacia, as traditionally circumscribed, contains at least several unrelated mimosoid lineages (Luckow et al., 2003; Miller et al., 2003).
  8. The Papilionoideae stem clade was set to a minimum of 55 Ma (Late Paleocene). This node is defined as the MRCA of Amburana cearensis and Cercidium floridum in the matK phylogeny (Fig. 1) and Peltogyne confertiflora and Glycine tabacina in the rbcL phylogeny (Kajita et al., 2001). The flowers of Barnebyanthus buchananensis are the earliest generalized papilionoid fossils from the Late Paleocene to Early Eocene sediments of the Buchanan Clay Pit, western Tennessee (Crepet and Herendeen, 1992). In addition to bilateral symmetry, these fossil flowers clearly show a connate calyx with five small calyx lobes, an adaxial median sepal, a differentiated standard petal that is positioned outside the wing petals, and distinct wing and keel petals, the latter of which are not fused. Crepet and Herendeen included Barnebyanthus in a cladistic analysis of traits scored also for selected genera of the tribe Sophoreae and concluded that Barnebyanthus was sister to the genistoid genera Bowringia and Clathrotropis. Although these fossil flowers show synapomorphies associated with the basic papilionoid zygomorphic flower, they actually lack other floral synapomorphies (e.g., fused keel petals, connate filaments) that could serve to unequivocally constrain their position within the papilionoid crown clade. Thus, the fossils are consistent with several lineages within "Sophoreae" and other papilionoids with free petals and filaments. This analysis of Crepet and Herendeen illustrates that legume fossils, because they include mainly disarticulated leaves, leaflets, flowers, and fruits, generally comprise few apomorphic characters.
    The minimum age of 55 Ma for the papilionoid stem clade is substantiated by a set of fossils that are suggestive of specific papilionoid groups. For example, fossils representing the lineage containing Diplotropis and Bowdichia (see node J, below) have an estimated age at around 56 Ma. Fossil leaves and fruits possibly representing Ormosia, are found in Colombia with an estimated age somewhere between 55 and 60 Ma (Wing et al., 2004).
  9. The Styphnolobium stem clade was set to a minimum of 40 Ma (Middle Eocene). This node is defined as the MRCA of Styphnolobium japonicum and Pickeringia montana in the matK phylogeny (Fig. 1) and is undefined in the rbcL phylogeny. Because both Styphnolobium and Cladrastis subgenus Cladrastis (here represented by Cladrastis delavayi, C. lutea, and C. sinensis) have a Middle Eocene fossil record from Tennessee (Herendeen, 1992), the stem leading to these two clades is constrained. Apomorphic traits that strongly suggest an affinity of these fossils with subgenus Cladrastis include a stipitate indehiscent fruit with thin pod walls and an attenuated nonplacental margin, a narrow wing along the placental suture, and seeds that are oriented longitudinally. Leaves and leaflets highly similar to those of subgenus Cladrastis co-occur with these same legume fruits in the Puryear Clay pit. Apomorphic traits that show an affinity of a single pod from the Miller Clay pit with Styphnolobium include indehiscent moniliform valves in which the constrictions between seeds narrow to less than half the diameter of the widest part of the pod. Most importantly, the seed-bearing portion is rounded in outline and longitudinally wrinkled, which suggests a diagnostic fleshy pericarp. The fossil fruit of Cladrastis subgenus Platycarpus from the Oligocene Bridge Creek flora (Brown, 1937) is superfluous to the Cladrastis subgenus Cladrastis and Styphnolobium fossils.
  10. The Diplotropis stem clade (essentially the Genistoid crown clade) was set to a minimum of 56 Ma, which is the oldest record for a definitive legume fossil. This node is defined as the MRCA of Diplotropis brasiliensis and Spartium junceum in the matK phylogeny (Fig. 2) and Acosmium dasycarpum and Spartium junceum in the rbcL phylogeny (Kajita et al., 2001). Fossil leaves and pods very similar to Bowdichia and Diplotropis are known from the Late Paleocene of western Wyoming (Herendeen and Wing, 2001), as well as the Middle Eocene of southeastern USA (Herendeen and Dilcher, 1990b). The apomorphic traits that suggest an affinity of these fossils to Bowdichia and Diplotropis include membranous pod valves, numerous seeds that are transversally oriented, a narrow wing on the placental suture, strong lateral nerves that demarcate the body (seed-containing region) from the placental and lower margins, and leaflets variably alternate to opposite (Calpurnia and Maackia of this same clade also have similar fruits). In the matK phylogeny (Fig. 2), the difference between the Bowdichia-Diplotropis stem and the Genistoid crown clade is trivial such that assignment of these fossils to one or the other node makes little difference with respect to rate and age estimation.
  11. The Machaerium stem clade was set to a minimum of 40 Ma. This node is defined as the MRCA of Aeschynomene purpusii and Machaerium falciforme in the matK phylogeny (Fig. 2) and Dalbergia hupeana and Machaerium lunatum in the rbcL phylogeny (Kajita et al., 2001). A series of fossil leaflets from the Puryear, Bolden, and Bovay Clay Pits (Herendeen, 1992) have numerous closely spaced craspedodromous secondary veins, a strong marginal vein, and poorly organized higher order venation that are diagnostic of a majority of species in the genus Machaerium. Epidermal cell structure fits the range of variation also observed in Machaerium.
  12. The Tipuana stem clade was set to a minimum of 10 Ma. This node is defined as the MRCA of Tipuana tipu and Maraniona lavinii in the matK phylogeny (Fig. 2) and is undefined in the rbcL phylogeny. Numerous winged fruits come from Loja and Cuenca in southern Ecuador, including Tipuana ecuatoriana (Burnham, 1995). This fossil fruit has a diagnostic long stipe, a unilateral wing emanating from mostly the style, parallel and distally spreading venation on wing, weak venation on stylar region of the wing, and seeds that are round in cross section (Burnham, 1995). Hughes et al. (2004) point out that the sessile fruits of Maraniona, which lack a pronounced wing venation, are highly similar to those of Tipuana, but distinctly lack the apomorphies shared between the fossil and extant species of Tipuana.
  13. The Robinia stem clade was set to a minimum of 34 Ma. This node is defined as the MRCA of Robinia pseudoacacia and Coursetia axillaris in the matK phylogeny (Fig. 3) and is undefined in the rbcL phylogeny. The fossil woods of Robinia zirkelii and the assignment of these to the Robinia stem clade are discussed in Lavin et al. (2003). This fossil wood species was widespread over North America and Europe and possesses such apomorphic traits as vestured intervessel pits, storied axial parenchyma and vessel elements, numerous thin-walled tyloses, ring-porous wood, homocellular rays, and spiral sculpturing in narrow vessels, which make assignment to the Robinia stem clade unequivocal.

Fossil Record of Fabales Other Than Leguminosae
With respect to the Fabales other than legumes (i.e., Polygalaceae, Surianaceae, and the genus Quillaja; sensu Angiosperm Phylogeny Group, 2003), fossil representation comes primarily from Polygalaceae, including several pollen records, fossil leaves and fruits currently under study (Kathleen Pigg, Arizona State University, personal communication). The fruits ("non-schizocarpic samaras") from the Late Paleocene Almont and Beicegel Creek (North Dakota) floras show similarities to Securidaca (Crane et al., 1990; Pigg et al., 2004). Polygalaceae has distinctive pollen diagnosed by more than four colporate apertures in which endoapertures may fuse. The exine stratification is obscure and surface is almost always smooth (Muller, 1981; Kathleen Pigg, Arizona State University, personal communications). Fossil pollen similar to the extant Monnina (Heusser, 1971) is described as Psilastephanocolporites fissilis from the Paleocene of Chile (Doubinger and Chotin, 1975; Muller, 1981) and from the Eocene of Guyana (Leidelmeyer, 1966). Polygalacidites calrus pollen is reported from the lower and upper Eocene of Assam, India (Sah and Dutta, 1968). Securidaca bombacopsis (Leopold and MacGinitie, 1972) is known from the Middle Eocene of North America, and pollen of Securidaca and a taxon similar to Bredemeyera are known from the upper Miocene of Mexico (Graham, 1976). Additional reports of polygalaceous pollen come from the Eocene of southeastern North America and California (Frederiksen, 1980, 1983; Taylor, 1990; Pigg et al., 2004).

These fossils are assigned to their taxonomic categories by using overall morphological similarity. Additionally, we did not sample Polygalaceae with respect to these fossils, so adding minimum constraints to nodes in the outgroups was not possible. Regardless, a minimum Late Paleocene age (Almont and Beicegel Creek; Crane et al., 1990; Manchester et al., in press; Pigg et al., 2004) assigned to the Polygalaceae stem clade (i.e., 55 Ma) does not affect the results presented because this clade is already constrained by older fixed ages (60 to 70 Ma) that are assigned to the coeval legume stem clade.

Surianaceae is now restricted to areas with subtropical to tropical coastal climates, but fossil wood with affinities to Suriana has been reported from the Eocene of Wyoming (Kruse, 1954). The tenuous taxonomic assignment of this fossil renders any age constraint dubious.

Potential fossil constraints may eventually come from groups related to Fabales. Wikstrom et al. (2001) assigned a 84-Ma fossil constraint to the split between Fagales and Cucurbitales (i.e., the Fagales stem clade). The three gene (Soltis et al., 2000) and matK (Hilu et al., 2003) phylogenies of angiosperms provide little resolution of the clades involving Fabales, Fagales, Cucurbitales, and Rosales. Fagales is sister to Cucurbitales in Soltis et al. (2000), but this larger clade is unresolved with respect to Rosales and Fabales. The Cucurbitales is sister to the other three orders in the matK strict consensus but with no clade support (Hilu et al., 2003). The uncertainty of the interrelationships among Fabales and related eurosid clades render assignments of the fossil constraint by Wikstrom et al. (2001) uncertain and in need of a separate study detailing the morphologies of extinct and extant eurosid taxa.


    Acknowledgements
 
Anne Bruneau, Alfonso Delgado-Salinas, Colin Hughes, R. Toby Pennington, and Mats Thulin provided unpublished matK sequences for selected taxa. Brian Schrire and Gwilym Lewis assisted with taxon sampling, organized input files involving tree descriptions and node identifications, and provided taxonomic details taken from their unpublished Legumes of the World. Scott Wing willingly showed us legume fossils from Cerrajón, Colombia. Kathleen Pigg kindly provided information on the fossil record of Fabales other than Leguminosae. Mike Crisp, Colin Hughes, Peter Linder, Rod Page, Kathleen Pigg, and Mike Sanderson provided critical comments that improved the manuscript. Jay Rotella guided us in dealing with the intricacies of model selection. This study was supported in part by grants from the National Science Foundation (DEB-9407824 to MFW and M. J. Sanderson, and DEB-0075202 to ML, and DEB-0316375 to PH) and funds provided by Arizona State University (to MFW). Useful discussion of this study was provided at the "Deep Time" sponsored symposium on dating methods (Botany 2003, Mobile, Alabama; NSF RCN-0090283).


    References
 Top
 Abstract
 Materials and Methods
 Results
 Discussion
 Appendix 1. The 13...
 References
 

    Albert V. A., Backlund A., Bremer K., Chase M. W., Manhart J., Mishler B. D., Nixon K. C. Functional constraints and rbcL evidence for land plant phylogeny. Ann. Mo. Bot. Garden (1994) 81:534–567.[CrossRef]

    Angiosperm Phylogeny Group, An update of the Angiosperm Phylogeny Group classification for the orders and families of flowering plants: APG II. Bot. J. Linn. Soc. (2003) 141:399–436. Angiosperm Phylogeny Group.[CrossRef][Web of Science]

    Axelrod D. I. Climatic pulses, a major factor in legume evolution. In: Advances in legume systematics, part 4, the fossil record—Herendeen P. S., Dilcher D. L., eds. (1992) Kew, UK: Royal Botanic Gardens. 259–279.

    Baldwin B. G., Sanderson M. J., Porter J. M., Wojciechowski M. F., Campbell C. S., Donoghue M. J. The ITS region of nuclear ribosomal DNA: A valuable source of evidence on angiosperm phylogeny. Ann. Mo. Bot. Garden (1995) 82:247–277.[CrossRef]

    Berggren W. A., Kent D. V., Swisher C. C. III, Aubry M.-P. A revised Cenozoic geochronology and chronostratigraphy. In: Geochronology, time scales and global stratigraphic correlation—Berggren W. A., Kent D. V., Aubry M.-P., Hardenbol J., eds. (1996) 129–212. SEPM Special Publication No 54.

    Brown R. W. Fossil legumes from Bridge Creek, Oregon. J. Wa. Acad. Sci. (1937) 27:414–418.

    Brown R. W. Paleocene flora of the Rocky Mountains and Great Plains (1962) Washington, D.C.: U. S. Government Printing Office. Geological Survey Professional Paper 375.

    Bruneau A., Forest F., Herendeen P. S., Klitgaard B. B., Lewis G. P. Phylogenetic relationships in the Caesalpinioideae (Leguminosae) as inferred from chloroplast trn L intron sequences. Syst. Bot. (2001) 26:487–514.

    Burnham K. P., Anderson D. R. Model selection and multimodel inference: A practical information-theoretic approach (2002) 2nd. New York: Springer-Verlag.

    Burnham R. A new species of winged fruit from the Miocene of Ecuador: Tipuana ecuatoriana (Leguminosae). Am. J. Bot. (1995) 82:1599–1607.[CrossRef][Web of Science]

    Crane P. R., Manchester S. R., Dilcher D. L. A preliminary survey of fossil leaves and well-preserved reproductive structures from the Sentinal Butte Formation (Paleocene) near Almont, North Dakota. Fieldiana Geol. New Ser. (1990) 20:1–63.

    Crepet W. L., Herendeen P. S. Papilionoid flowers from the early Eocene of southeastern North America. In: Advances in legume systematics, part 4, the fossil record—Herendeen P. S., Dilcher D. L., eds. (1992) Kew, UK: Royal Botanic Gardens. 43–55.

    Crepet W. L., Taylor D. W. The diversification of the Leguminosae: First fossil evidence of the Mimosoideae and Papilionoideae. Science (1985) 288:1087–1089.

    Crepet W. L., Taylor D. W. Primitive mimosoid flowers from the Paleocene-Eocene and their systematic and evolutionary implications. Am. J. Bot. (1986) 73:548–563.[CrossRef][Web of Science]

    Crisp M. D., Gilmore S., van Wyk B.-E. Molecular phylogeny of the genistoid tribes of papilionoid legumes. In: Advances in legume systematics, part 9—Herendeen P. S., Bruneau A., eds. (2000) Kew, UK: Royal Botanic Garden. 249–276.

    Delgado-Salinas A., Bibler R., Lavin M. Molecular phylogeny of Phaseolus. Abstract. (2003) The Seventeenth Biennial Meeting of The Bean Improvement Cooperative: SacramentoCalifornia.

    Dilcher D. L., Herendeen P. S., Hueber F. Fossil Acacia flowers with attached anther glands from Dominican Republic amber. In: Advances in legume systematics, part 4, the fossil record—Herendeen P. S., Dilcher D. L., eds. (1992) Kew, UK: Royal Botanic Gardens. 33–42.

    Doubinger J., Chotin P. Etude palynologique de lignites Tertiaires du basin d'Arauco-Concepción (Chile). Rev. Esp. Micropaleont. (1975) 7:549–565.

    Engel M. S. Monophyly and extensive extinction of advanced eusocial bees: Insights from an unexpected Eocene diversity. Proc. Nat. Acad. Sci. USA (2001) 98:1661–1664.[Abstract/Free Full Text]

    Felsenstein J. Phylogenies from molecular sequences: Inference and reliability. Ann. Rev. Genet. (1988) 22:521–565.[CrossRef][Web of Science][Medline]

    Frederiksen N. O. Sporomorphs from the Jackson Group (Upper Eocene) and adjacent strata of Mississippi and western Alabama. U.S. Geological Survey Professional Paper (1980) 1084:1–75.

    Frederiksen N. O. Middle Eocene palynomorphs from San Diego, California. Part II. Angiosperm pollen and miscellanea. American Association Stratigraphic Palynological Contributions Series (1983) 12:32–155.

    Graham A. Studies in Neotropical paleobotany. II. The Miocene communities of Veracruz, Mexico. Ann. Mo. Botanical Garden (1976) 63:787–842.

    Graham A. The current status of the legume fossil record in the Caribbean region. In: Advances in legume systematics, part 4, the fossil record—Herendeen P. S., Dilcher D. L., eds. (1992) Kew, UK: Royal Botanic Gardens. 161–167.

    Gros J. P. A synopsis of the fossil record of mimosoid legume wood. In: Advances in legume systematics, part 4, the fossil record—Herendeen P. S., Dilcher D. L., eds. (1992) Kew, UK: Royal Botanic Gardens. 69–83.

    Herendeen P. S. The fossil record of the Leguminosae: Recent advances. (2001) Legumes down under: The Fourth International Legume Conference, Abstracts, 34–35. Canberra, Australia: Australian National University.

    Herendeen P. S., Bruneau A., Lewis G. P. Phylogenetic relationships in caesalpinioid legumes: A preliminary analysis based on morphological and molecular data. In: Advances in legume systematics, part 10, higher level systematics—Klitgaard B. B., Bruneau A., eds. (2003a) Kew, UK: Royal Botanic Gardens. 37–62.

    Herendeen P. S., Crepet W. L., Dilcher D. L. The fossil history of the Leguminosae: Phylogenetic and biogeographic implications. In: Advances in legume systematics, part 4, the fossil record—Herendeen P. S., Dilcher D. L., eds. (1992) Kew, UK: Royal Botanic Gardens. 303–316.

    Herendeen P. S., Dilcher D. L. Fossil mimosoid legumes from the Eocene and Oligocene of southeastern North America. Rev. Palaeobot. Palynol. (1990a) 62:339–361.[CrossRef][Web of Science]

    Herendeen P. S., Dilcher D. L. Diplotropis (Leguminosae, Papilionoideae) from the Middle Eocene of southeastern North America. Syst. Bot. (1990b) 15:526–533.[CrossRef]

    Herendeen P. S., Dilcher D. L. Caesalpinia subgenus Mezoneuron (Leguminosae, Caesalpinioideae) from the Tertiary of North America. Am. J. Bot. (1991) 78:1–12.[CrossRef][Web of Science]

    Herendeen P. S., Dilcher D. L. Advances in legume systematics, part 4. The fossil record (1992) Kew, UK: Royal Botanic Gardens.

    Herendeen P. S., Jacobs B. F. Fossil legumes from the Middle Eocene (460 Ma) Mahenge Flora of Singida, Tanzania. Am. J. Bot. (2000) 87:1358–1366.[Abstract/Free Full Text]

    Herendeen P. S., Lewis G. P., Bruneau A. Floral morphology in caesalpinioid legumes: Testing the monophyly of the "Umtiza clade". Int. J. Plant Sci. (2003b) 164:S393–S407.[CrossRef][Web of Science]

    Herendeen P. S., Wing S. Papilionoid legume fruits and leaves from the Paleocene of northwestern Wyoming (2001) Botanical Society of America. http://www.botany2001.org/Botany 2001 Abstracts, published by.

    Herendeen P. S., Zarucchi J. L. Validation of Caesalpinia subgenus Mezoneuron (Desf.) Vidal and new combinations in Caesalpinia for two species of Mezoneuron from Africa. Ann. Mo. Bot. Garden (1990) 77:854–855.[CrossRef]

    Heusser C. J. Pollen and spores of Chile (1971) Tucson: University of Arizona Press.

    Hilu K. W., Borsch T., Müller K., Soltis D. E., Soltis P. S., Savolainen V., Chase M. W., Powell M. P., Alice L. A., Evans R., Sauquet H., Neinhuis C., Slotta T. A. B., Rohwer J. G., Campbell C. S., Chatrou L. W. Angiosperm phylogeny based on matK sequence information. Am. J. Bot. (2003) 90:1758–1776.[Abstract/Free Full Text]

    Hu J.-M., Lavin M., Wojciechowski M. F., Sanderson M. J. Phylogenetic analysis of nuclear ribosomal ITS/5.8 S sequences in the tribe Millettieae (Fabaceae): Poecilanthe - Cyclolobium, the core Millettieae, and the Callerya group. Sys. Bot. (2002) 27:722–733.

    Hueber F. M., Langenheim J. Dominican amber tree had African ancestors. Geotimes (1986) 31:8–10.

    Huelsenbeck J. P., Crandall K. A. Phylogeny estimation and hypothesis testing using maximum likelihood. Ann. Rev. Ecol. Sys. (1997) 28:437–466.[CrossRef]

    Huelsenbeck J. P., Ronquist F. R. MrBayes: Bayesian inference of phylogeny. Bioinformatics (2001) 17:754.[Abstract/Free Full Text]

    Huelsenbeck J. P., Ronquist R., Nielson R., Bollback J. P. Bayesian inference of phylogeny and its impact on evolutionary biology. Science (2001) 294:2310–2314.[Abstract/Free Full Text]

    Hughes C. E., Lewis G. P., Daza Yamona A., Reynel C. Maraniona, a new dalbergioid legume genus (Leguminosae, Papilionoideae) from Peru. Sys. Bot. (2004) 29:366–374.[CrossRef]

    Iturralde-Vinent M. A., MacPhee R. D. E. Age and paleogeographical origin of Dominican amber. Science (1996) 273:1850–1852.[Abstract/Free Full Text]

    Jacobs B. F. The plant fossil record and implications for phytogeography in tropical Africa. XVIIth AETFAT Congress Abstracts: 47 (2003) Addis Ababa: Addis Ababa University.

    Johnson J. B., Omland K. S. Model selection in ecology and evolution. Trends Ecol. Evol. (2004) 19:101–108.[CrossRef][Medline]

    Kajita T., Ohashi H., Tateishi Y., Bailey C. D., Doyle J. J. rbcL and legume phylogeny, with particular reference to Phaseoleae, Millettieae, and Allies. Sys. Bot. (2001) 26:515–536.

    Krell F.-T., Cranston P. S. Which side of the tree is more basal? Sys. Entomol. (2004) 29:279–281.[CrossRef]

    Kruse H. O. Some Eocene dicotyledonous woods from the Eden Valley, Wyoming. Ohio J. Sci. (1954) 54:243–268.

    Langley C. H., Fitch W. An estimation of the constancy of the rate of molecular evolution. J. Mol. Evol. (1974) 3:161–177.[CrossRef][Web of Science][Medline]

    Lavin M., Pennington R. T., Klitgaard B. B., Sprent J. I., De Lima H. C., Gasson P. E. The dalbergioid legumes (Fabaceae): Delimitation of a monophyletic pantropical clade. Am. J. Bot. (2001) 88:503–533.[Abstract/Free Full Text]

    Lavin M., Schrire B. D., Lewis G. P., Pennington R. T., Delgado-Salinas A., Thulin M., Hughes C. E., Beyra Matos A., Wojciechowski M. F. Metacommunity process rather than continental tectonic history better explains geographically structured phylogenies in legumes. Phil. Trans. R. Soc. Biol. Sci. (2004) 359(1450):1509–1522.[CrossRef]

    Lavin M., Wojciechowski M. F., Gasson P., Hughes C. E., Wheeler E. Phylogeny of robinioid legumes (Fabaceae) revisited: Coursetia and Gliricidia recircumscribed, and a biogeographical appraisal of the Caribbean endemics. Sys. Bot. (2003) 28:387–409.

    Leidelmeyer P. The Paleocene and Lower Eocene pollen flora of Guyana. Leidse Geologische Medelingen (1966) 36:49–70.

    Leopold E. B., MacGinitie H. D. Development and affinities of Tertiary floras in the Rocky Mountains. In: Floristics and paleofloristics of Asia and eastern North America—Graham A., ed. (1972) Amsterdam: Elsevier Publishing Company. 147–200.

    Lewis G., Schrire B., Mackinder B., Lock M., eds. Legumes of the world (2005) Kew: Royal Botanic Gardens. 592.

    Luckow M., Miller J. T., Murphy D. J., Livshultz T. A phylogenetic analysis of the Mimosoideae (Leguminosae) based on chloroplast DNA sequence data. In: Advances in legume systematics, part 10, higher level systematics—Klitgaard B. B., Bruneau A., eds. (2003) Kew, UK: Royal Botanic Gardens. 197–220.

    MacGinitie H. D. Fossil plants of the Florissant Beds (1953) 599 i–iii. Colorado: Carnegie Institution of Washington Publications. 1–198.

    Magallón S., Sanderson M. J. Absolute diversification rates in angiosperms clades. Evolution (2001) 55:1762–1780.[CrossRef][Web of Science][Medline]

    Manchester S. R. Update on the megafossil flora of Florissant, Colorado. Denver Museum of Nature & Science, Series 4, (2001) 1:137–161.

    Manchester S. R., Pigg K. B., Crane P. R. Palaeocarpinus dakotaensis sp. nov. (Betulaceae: Coryloideae) and associated staminate catkins, pollen and leaves from the Paleocene of North Dakota. Int. J. Plant Sci. (2004) 165:1135–1148.[CrossRef][Web of Science]

    Manos P. S., Stanford A. M. The biogeography of Fagaceae: Tracking the Tertiary history of temperate and subtropical forests of the Northern Hemisphere. Int. J. Plant Sci. (2001) 162:S77–S93.[CrossRef][Web of Science]

    McClain A. M., Manchester S. R. Dipteronia (Sapindaceae) from the Tertiary of North America and implications for the phytogeographic history of the Aceroideae. Am. J. Bot. (2001) 88:1316–1325.[Abstract/Free Full Text]

    McKey D. Legumes and nitrogen: The evolutionary ecology of a nitrogen-demanding lifestyle. In: Advances in legume systematics, part 5, the nitrogen factor—Sprent J. I., McKey D., eds. (1994) Kew, UK: Royal Botanic Gardens. 211–228.

    Miller J. T., Grimes J. W., Murphy D. J., Bayer R. J., Ladiges P. Y. A phylogenetic analysis of the Acacieae and Ingeae (Mimosoideae: Fabaceae) based on trnK, matK, psbA, trnH, and trnL trnF sequence data. Syst. Bot. (2003) 28:558–566.

    Muller J. Fossil pollen records of extant angiosperms. Bot. Rev. (1981) 47:1–142.

    Neel M. C., Cummings M. P. Section-level relationships of North American Agalinis (Orobanchaceae) based on DNA sequence analysis of three chloroplast gene regions. BMC Evol. Biol. (2004) 4:15.[CrossRef][Medline]

    Osaloo S. K., Kawano S. Molecular systematics of Trilliaceae II. Phylogenetic analyses of Trillium and its allies using sequences of rbcL and matK genes of cpDNA and internal transcribed spacers of 18S–26S nrDNA. Plant Species Bio. (1999) 14:75–94.

    Pennington R. T., Lavin M., Ireland H., Klitgaard B., Preston J., Hu J.-M. Phylogenetic relationships of basal papilionoid legumes based upon sequence of the chloroplast trn L intron. Syst. Bot. (2001) 26:537–556.

    Percy D. M., Page R. D. M., Cronk Q. C. B. Plant-insect interactions: Double-dating associated insect and plant lineage reveals asynchrous radiations. Syst. Biol. (2004) 53:120–127.[Abstract/Free Full Text]

    Pigg K. B., Wojciechowski M. F., DeVore M. L. Samaras from the Late Paleocene Almont and Beicegel Creek floras of North Dakota, U.S.A. with potential affinities to Securidaca (Polygalaceae). Botany 2004 (2004) abstract.

    Poinar G. O. Jr. Hymenaea protera sp.n. (Leguminosae, Caesalpinioideae) from Dominican amber has African affinities. Experientia (1991) 47:1075–1082.[CrossRef][Web of Science]

    Poinar G. O. Jr., Brown A. E. Hymenaea mexicana sp. nov. (Leguminosae: Caesalpinioideae) from Mexican amber indicates Old World connections. Bot. J. Linn. Soc. (2002) 139:125–132.[CrossRef][Web of Science]

    Poinar H. N., Cano R. J., Poinar G. O. Jr. DNA from an extinct plant. Nature (1993) 363:677.[CrossRef]

    Polhill R. M. Classification of the Leguminosae. In: Phytochemical dictionary of the Leguminosae—Bisby F. A., Buckingham J., Harborne J. B., eds. (1994) New York: Chapman and Hall. xxxv–xlviii.

    Polhill R. M., Raven P. H., Stirton C. H. Evolution and systematics of the Leguminosae. In: Advances in legume systematics, part 1—Polhill R. M., Raven P. H., eds. (1981) Kew, UK: Royal Botanic Gardens. 1–26.

    Posada D., Crandall K. A. ModelTest: Testing the model of DNA substitution. Bioinformatics (1998) 14:817–818.[Abstract/Free Full Text]

    Riley-Hulting E., Delgado-Salinas A., Lavin M. Phylogenetic systematics of Strophostyles (Fabaceae): A North American temperate genus within a Neotropical diversification. Syst. Bot. (2004) 29:627–653.[CrossRef]

    Sah S. C. D., Dutta K. Palynostratigraphy of the Tertiary sedimentary formations of Assam. 2. Stratigraphic significance of spores and pollen in the Tertiary succession of Assam. Palaeobotanist (1968) 16:177–195.

    Sanderson M. A nonparametric approach to estimating divergence times in the absence of rate constancy. Mol. Biol. Evol. (1997) 14:1218–1231.[Web of Science]

    Sanderson M. Estimating absolute rates of molecular evolution and divergence times: A penalized likelihood approach. Mol. Biol. Evol. (2002) 19:101–109.[Abstract/Free Full Text]

    Sanderson M. (2003) Davis: University of California. r8s, version 1.6, User's Manual (April 2003). Distributed by the author (http://ginger.ucdavis.edu/r8s/).

    Schneider H., Schuettpelz E., Pryer K. M., Cranfill R., Magallón S., Lupia R. Ferns diversified in the shadow of angiosperms. Nature (2004) 428:553–557.[CrossRef][Medline]

    Schrire B. D., Lavin M., Barker N. P., Cortes-Burns H., von Senger I., Kim J.-H. Towards a phylogeny of Indigofera (Leguminosae-Papilionoideae): Identification of major clades and relative ages. In: Advances in legume systematics, part 10, higher level systematics—Klitgaard B. B., Bruneau A., eds. (2003) Kew, UK: Royal Botanic Gardens. 269–302.

    Soltis D. E., Soltis P. S., Chase M. W., Mort M. E., Albach D. C., et al. Angiosperm phylogeny inferred from 18S rDNA, rbcL, and atpB sequences. Bot. J. Linn. Soc. (2000) 133:381–461.[CrossRef][Web of Science]

    Sprent J. I. Nodulation in legumes (2001) Kew, UK: Royal Botanic Gardens.

    Steele K. P., Vilgalys R. Phylogenetic analyses of Polemoniaceae using nucleotide sequences of the plastid gene matK. Syst. Bot. (1994) 19:126–142.[CrossRef]

    Taylor D. W. Paleobiogeographic relationships from the Cretaceous and Early Tertiary of the North American area. Bot. Rev. (1990) 56:279–417.[CrossRef]

    Thulin M., Lavin M., Pasquet R., Delgado-Salinas A. Phylogeny and biogeography of Wajira (Leguminosae): A monophyletic segregate of Vigna centered in the Horn of Africa region. Syst. Bot. (2004) 29:903–920.[CrossRef]

    Tucker S. C., Douglas A. W. Ontogenetic evidence and phylogenetic relationships among basal taxa of legumes. In: Advances in legume systematics, part 6, structural botany—Ferguson I. K., Tucker S. C., eds. (1994) Kew, UK: Royal Botanic Gardens. 11–32.

    Wang X.-R., Tsumura Y., Yoshimaru H., Nagasaka K., Szmidt A. E. Phylogenetic relationships of Eurasian Pines ( Pinus, Pinaceae) based on chloroplast rbcL, matK, rpl20-rps18 spacer, and trnV intron sequences. Am. J. Bot. (1999) 86:1742–1753.[Abstract/Free Full Text]

    Webster A. J., Payne R. J. H., Pagel M. Molecular phylogenies link rates of evolution and speciation. Science (2003) 301:478.[Free Full Text]

    Webster A. J., Payne R. J. H., Pagel M. Response to comments on "Molecular phylogenies link rates of evolution and speciation". Science (2004) 303:173d.[Free Full Text]

    Wikström N., Savolainen V., Chase M. W. Evolution of angiosperms: Calibrating the family tree. Proc. R. Soc. B (2001) 268:2211–2220.[Medline]

    Wilf P., Labandeira C. C. Response of plant-insect associations to Paleocene-Eocene warming. Science (1999) 284:2153–2156.[Abstract/Free Full Text]

    Wing S. L., Herrera F., Jaramillo C. A Paleocene flora from the Cerrajón Formation, Guajíra Peninsula, northeastern Colombia. (2004) VII International Organization of Paleobotany Conference Abstracts. Trelew, Argentina: Museo Egidio Feruglio. 146–147.

    Wojciechowski M. F. Reconstructing the phylogeny of legumes (Leguminosae): An early 21st century perspective. In: Advances in legume systematics, part 10, Higher level systematics—Klitgaard B. B., Bruneau A., eds. (2003) Kew, UK: Royal Botanic Gardens. 5–35.

    Wojciechowski M. F., Lavin M., Sanderson M. J. A phylogeny of legumes (Leguminosae) based on analysis of the plastid matK gene resolves many well-supported subclades within the family. Am. J. Bot. (2004) 91:1846–1862.[Abstract/Free Full Text]

    Wojciechowski M. F., Sanderson M. J., Steele K. P., Liston A. Molecular phylogeny of the "temperate herbaceous tribes" of papilionoid legumes: A supertree approach. In: Advances in legume systematics, part 9—Herendeen P. S., Bruneau A., eds. (2000) Kew, UK: Royal Botanic Gardens. 277–298.

    Wolfe J. A., Tanai T. Systematics, phylogeny, and distribution of Acer (maples) in the Cenozoic of western North America. Journal of the Faculty of Science of Hokkaido University Series IV (1987) 22:1–10.

    Xiang Q.-Y., Soltis D. E., Soltis P. S., Manchester S. R., Crawford D. J. Timing the eastern Asian–eastern North American floristic disjunction: Molecular clock corroborates paleontological estimates. Mol. Phylogenet. Evol. (2000) 15:462–472.[CrossRef][Web of Science][Medline]

    Yang Z., Yoder A. D. Comparison of likelihood and Bayesian methods for estimating divergence times using multiple gene loci and calibration points, with application to a radiation of cute-looking mouse legume species. Syst. Biol. (2003) 15:705–716.

    Yoder A. D., Yang Z. Divergence dates for Malagasy lemurs estimated from multiple gene loci: Geological and evolutionary context. Mol. Ecol. (2004) 13:757–773.[CrossRef][Medline]

    Young N. D., dePamphilis C. W. Purifying selection detected in the plastid gene matK and flanking ribozyme regions within a group II intron of nonphotosynthetic plants. Mol. Biol. Evol. (2000) 17:1933–1941.[Abstract/Free Full Text]

    Zhang L.-B., Renner S. The deepest splits in Chloranthaceae as resolved by chloroplast sequences. Int. J. Plant Sci. (2003) 164(5 Suppl.):S383–S392.[CrossRef][Web of Science]


Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us    What's this?


This article has been cited by other articles:


Home page
Proc. Natl. Acad. Sci. USAHome page
S. L. Wing, F. Herrera, C. A. Jaramillo, C. Gomez-Navarro, P. Wilf, and C. C. Labandeira
Late Paleocene fossils from the Cerrejon Formation, Colombia, are the earliest record of Neotropical rainforest
PNAS, November 3, 2009; 106(44): 18627 - 18632.
[Abstract] [Full Text] [PDF]


Home page
Plant Physiol.Home page
S. B. Cannon, G. D. May, and S. A. Jackson
Three Sequenced Legume Genomes and Many Crop Species: Rich Opportunities for Translational Genomics
Plant Physiology, November 1, 2009; 151(3): 970 - 977.
[Full Text] [PDF]


Home page
Plant Physiol.Home page
N. Gill, S. Findley, J. G. Walling, C. Hans, J. Ma, J. Doyle, G. Stacey, and S. A. Jackson
Molecular and Chromosomal Evidence for Allopolyploidy in Soybean
Plant Physiology, November 1, 2009; 151(3): 1167 - 1174.
[Abstract] [Full Text] [PDF]


Home page
Plant Physiol.Home page
P. David, N. W.G. Chen, A. Pedrosa-Harand, V. Thareau, M. Sevignac, S. B. Cannon, D. Debouck, T. Langin, and V. Geffroy
A Nomadic Subtelomeric Disease Resistance Gene Cluster in Common Bean
Plant Physiology, November 1, 2009; 151(3): 1048 - 1065.
[Abstract] [Full Text] [PDF]


Home page
Proc. Natl. Acad. Sci. USAHome page
J. Peccoud, J.-C. Simon, H. J. McLaughlin, and N. A. Moran
Post-Pleistocene radiation of the pea aphid complex revealed by rapidly evolving endosymbionts
PNAS, September 22, 2009; 106(38): 16315 - 16320.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Bot.Home page
B. D. Schrire, M. Lavin, N. P. Barker, and F. Forest
Phylogeny of the tribe Indigofereae (Leguminosae-Papilionoideae): Geographically structured more in succulent-rich and temperate settings than in grass-rich environments
Am. J. Botany, April 1, 2009; 96(4): 816 - 852.
[Abstract] [Full Text] [PDF]


Home page
Proc. Natl. Acad. Sci. USAHome page
T. Kawashima, X. Wang, K. F. Henry, Y. Bi, K. Weterings, and R. B. Goldberg
Identification of cis-regulatory sequences that activate transcription in the suspensor of plant embryos
PNAS, March 3, 2009; 106(9): 3627 - 3632.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Bot.Home page
D. E. Soltis, V. A. Albert, J. Leebens-Mack, C. D. Bell, A. H. Paterson, C. Zheng, D. Sankoff, C. W. dePamphilis, P. K. Wall, and P. S. Soltis
Polyploidy and angiosperm diversification
Am. J. Botany, January 1, 2009; 96(1): 336 - 348.
[Abstract] [Full Text] [PDF]


Home page
Plant Physiol.Home page
R. W. Innes, C. Ameline-Torregrosa, T. Ashfield, E. Cannon, S. B. Cannon, B. Chacko, N. W.G. Chen, A. Couloux, A. Dalwani, R. Denny, et al.
Differential Accumulation of Retroelements and Diversification of NB-LRR Disease Resistance Genes in Duplicated Regions following Polyploidy in the Ancestor of Soybean
Plant Physiology, December 1, 2008; 148(4): 1740 - 1759.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Bot.Home page
G. Prenner and B. B. Klitgaard
Towards unlocking the deep nodes of Leguminosae: Floral development and morphology of the enigmatic Duparquetia orchidacea (Leguminosae, Caesalpinioideae)
Am. J. Botany, November 1, 2008; 95(11): 1349 - 1365.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Bot.Home page
C. Saslis-Lagoudakis, M. W. Chase, D. N. Robinson, S. J. Russell, and B. B. Klitgaard
Phylogenetics of neotropical Platymiscium (Leguminosae: Dalbergieae): systematics, divergence times, and biogeography inferred from nuclear ribosomal and plastid DNA sequence data
Am. J. Botany, October 1, 2008; 95(10): 1270 - 1286.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Bot.Home page
R. A. Scherson, R. Vidal, and M. J. Sanderson
Phylogeny, biogeography, and rates of diversification of New World Astragalus (Leguminosae) with an emphasis on South American radiations
Am. J. Botany, August 1, 2008; 95(8): 1030 - 1039.
[Abstract] [Full Text] [PDF]


Home page
Syst BiolHome page
I. J. Maureira-Butler, B. E. Pfeil, A. Muangprom, T. C. Osborn, and J. J. Doyle
The Reticulate History of Medicago (Fabaceae)
Syst Biol, June 1, 2008; 57(3): 466 - 482.
[Abstract] [Full Text] [PDF]


Home page
DNA ResHome page
K. Van, D. H. Kim, C. M. Cai, M. Y. Kim, J. H. Shin, M. A. Graham, R. C. Shoemaker, B.-S. Choi, T.-J. Yang, and S.-H. Lee
Sequence Level Analysis of Recently Duplicated Regions in Soybean [Glycine max (L.) Merr.] Genome
DNA Res, April 1, 2008; 15(2): 93 - 102.
[Abstract] [Full Text] [PDF]


Home page
J HeredHome page
J. A. Schlueter, B. E. Scheffler, S. Jackson, and R. C. Shoemaker
Fractionation of Synteny in a Genomic Region Containing Tandemly Duplicated Genes across Glycine max, Medicago truncatula, and Arabidopsis thaliana
J. Hered., March 2, 2008; (2008) esn010v1.
[Abstract] [Full Text] [PDF]


Home page
DNA ResHome page
H. Hisano, S. Sato, S. Isobe, S. Sasamoto, T. Wada, A. Matsuno, T. Fujishiro, M. Yamada, S. Nakayama, Y. Nakamura, et al.
Characterization of the Soybean Genome Using EST-derived Microsatellite Markers
DNA Res, January 11, 2008; (2008) dsm025v1.
[Abstract] [Full Text] [PDF]


Home page
DNA ResHome page
Z. Xia, Y. Tsubokura, M. Hoshi, M. Hanawa, C. Yano, K. Okamura, T. A. Ahmed, T. Anai, S. Watanabe, M. Hayashi, et al.
An Integrated High-density Linkage Map of Soybean with RFLP, SSR, STS, and AFLP Markers Using A Single F2 Population
DNA Res, January 11, 2008; (2008) dsm027v1.
[Abstract] [Full Text] [PDF]


Home page
Plant Physiol.Home page
C. Ameline-Torregrosa, B.-B. Wang, M. S. O'Bleness, S. Deshpande, H. Zhu, B. Roe, N. D. Young, and S. B. Cannon
Identification and Characterization of Nucleotide-Binding Site-Leucine-Rich Repeat Genes in the Model Plant Medicago truncatula
Plant Physiology, January 1, 2008; 146(1): 5 - 21.
[Abstract] [Full Text] [PDF]


Home page
Plant CellHome page
C. E.M. Champagne, T. E. Goliber, M. F. Wojciechowski, R. W. Mei, B. T. Townsley, K. Wang, M. M. Paz, R. Geeta, and N. R. Sinha
Compound Leaf Development and Evolution in the Legumes
PLANT CELL, November 1, 2007; 19(11): 3369 - 3378.
[Abstract] [Full Text] [PDF]


Home page
Syst BiolHome page
P. S. Manos, P. S. Soltis, D. E. Soltis, S. R. Manchester, S.-H. Oh, C. D. Bell, D. L. Dilcher, and D. E. Stone
Phylogeny of Extant and Fossil Juglandaceae Inferred from the Integration of Molecular and Morphological Data Sets
Syst Biol, June 1, 2007; 56(3): 412 - 430.
[Abstract] [Full Text] [PDF]


Home page
Plant Physiol.Home page
B. H. Le, J. A. Wagmaister, T. Kawashima, A. Q. Bui, J. J. Harada, and R. B. Goldberg
Using Genomics to Study Legume Seed Development
Plant Physiology, June 1, 2007; 144(2): 562 - 574.
[Full Text] [PDF]


Home page
Crop Sci.Home page
S. A. Jackson, D. Rokhsar, G. Stacey, R. C. Shoemaker, J. Schmutz, and J. Grimwood
Toward a Reference Sequence of the Soybean Genome: A Multiagency Effort
Crop Sci., November 1, 2006; 46(Supplement_1): S-55 - S-61.
[Abstract] [Full Text] [PDF]


Home page
Proc. Natl. Acad. Sci. USAHome page
S. B. Cannon, L. Sterck, S. Rombauts, S. Sato, F. Cheung, J. Gouzy, X. Wang, J. Mudge, J. Vasdewani, T. Schiex, et al.
Legume genome evolution viewed through the Medicago truncatula and Lotus japonicus genomes
PNAS, October 3, 2006; 103(40): 14959 - 14964.
[Abstract] [Full Text] [PDF]


Home page
Syst BiolHome page
C. Yesson and A. Culham
Phyloclimatic Modeling: Combining Phylogenetics and Bioclimatic Modeling
Syst Biol, October 1, 2006; 55(5): 785 - 802.
[Abstract] [Full Text] [PDF]


Home page
Proc. Natl. Acad. Sci. USAHome page
C. Hughes and R. Eastwood
From the Cover: Island radiation on a continental scale: Exceptional rates of plant diversification after uplift of the Andes
PNAS, July 5, 2006; 103(27): 10334 - 10339.
[Abstract] [Full Text] [PDF]


Home page
Phil Trans R Soc BHome page
S. Conway Morris
Darwin's dilemma: the realities of the Cambrian 'explosion'
Phil Trans R Soc B, June 29, 2006; 361(1470): 1069 - 1083.
[Abstract] [Full Text] [PDF]


Home page
Proc. Natl. Acad. Sci. USAHome page
Q. C. B. Cronk
Legume flowers bear fruit.
PNAS, March 28, 2006; 103(13): 4801 - 4802.
[Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (68)
Right arrowRequest Permissions
Google Scholar
Right arrow Articles by Lavin, M.
Right arrow Articles by Wojciechowski, M. F.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Lavin, M.
Right arrow Articles by Wojciechowski, M. F.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?