Non-equilibrium dynamics and floral trait interactions shape extant angiosperm diversity.

Why are some traits and trait combinations exceptionally common across the tree of life, whereas others are vanishingly rare? The distribution of trait diversity across a clade at any time depends on the ancestral state of the clade, the rate at which new phenotypes evolve, the differences in speciation and extinction rates across lineages, and whether an equilibrium has been reached. Here we examine the role of transition rates, differential diversification (speciation minus extinction) and non-equilibrium dynamics on the evolutionary history of angiosperms, a clade well known for the abundance of some trait combinations and the rarity of others. Our analysis reveals that three character states (corolla present, bilateral symmetry, reduced stamen number) act synergistically as a key innovation, doubling diversification rates for lineages in which this combination occurs. However, this combination is currently less common than predicted at equilibrium because the individual characters evolve infrequently. Simulations suggest that angiosperms will remain far from the equilibrium frequencies of character states well into the future. Such non-equilibrium dynamics may be common when major innovations evolve rarely, allowing lineages with ancestral forms to persist, and even outnumber those with diversification-enhancing states, for tens of millions of years.

Consistent phenological shifts in the making of a biodiversity hotspot: the Cape flora.

BACKGROUND: The best documented survival responses of organisms to past climate change on short (glacial-interglacial) timescales are distributional shifts. Despite ample evidence on such timescales for local adaptations of populations at specific sites, the long-term impacts of such changes on evolutionary significant units in response to past climatic change have been little documented. Here we use phylogenies to reconstruct changes in distribution and flowering ecology of the Cape flora--South Africa's biodiversity hotspot--through a period of past (Neogene and Quaternary) changes in the seasonality of rainfall over a timescale of several million years. RESULTS: Forty-three distributional and phenological shifts consistent with past climatic change occur across the flora, and a comparable number of clades underwent adaptive changes in their flowering phenology (9 clades; half of the clades investigated) as underwent distributional shifts (12 clades; two thirds of the clades investigated). Of extant Cape angiosperm species, 14-41% have been contributed by lineages that show distributional shifts consistent with past climate change, yet a similar proportion (14-55%) arose from lineages that shifted flowering phenology. CONCLUSIONS: Adaptive changes in ecology at the scale we uncover in the Cape and consistent with past climatic change have not been documented for other floras. Shifts in climate tolerance appear to have been more important in this flora than is currently appreciated, and lineages that underwent such shifts went on to contribute a high proportion of the flora's extant species diversity. That shifts in phenology, on an evolutionary timescale and on such a scale, have not yet been detected for other floras is likely a result of the method used; shifts in flowering phenology cannot be detected in the fossil record.

Floral organogenesis and the developmental basis for pollinator deception in the asiatic dayflower, Commelina communis (Commelinaceae).

The upper half of flowers in Commelina communis deceptively lures potential pollinators with its showy petals and staminodes on the false promise of abundant pollen. This paper presents evidence that staminodization in the upper half is associated with a severe retardation of the entire upper floral hemisphere early in development. Possible consequences of this developmental retardation are seen also in the gynoecium, where the upper carpel of the three-carpellate ovary is underdeveloped and sterile at maturity. Only late in development do the upper petals and staminodes expand and acquire pigments necessary for their attractive function. We surmise that retardations of this severity are unlikely to be found for functionally fertile organs such as stamens and ovule-producing carpels, because key preparatory events preceding sporogenesis might otherwise be disrupted. Such differential growth about the floral apex resembles that known in some eudicots to be regulated by the TCP gene family; thus, future molecular developmental studies in Commelina may help to extend our understanding of the evolutionary genetics of floral monosymmetry to monocots.

Intraspecific variability and timing in ancestral ecology reconstruction: a test case from the cape flora.

Thamnochortus (ca. 32 species) is an ecologically diverse genus of Restionaceae. Restionaceae comprise a major component of the southern African Cape flora, wherein eco-diversification might have been important in the generation of high levels of species richness. In an attempt to reconstruct the macroecological history of Thamnochortus, it was found that standard procedures for character state optimization make two inappropriate assumptions. The first is that ancestors are monomorphic (i.e., ecologically uniform) and the second is that eco-diversification follows, or is slower than, lineage diversification. We demonstrate a variety of coding schemes with which the assumption of monomorphy can be avoided. For unordered discrete ecological characters, presence coding and generalized frequency coding (GFC) are suboptimal because they occasionally yield illogical assignments of no state to ancestors. Polymorphism coding or use of the program DIVA are preferable in this respect but are applicable only with parsimony. For continuous eco-characters (e.g., a rainfall gradient, where individual species occur in ranges), GFC and MaxMin coding provide equally valid solutions to optimizing ranges with parsimony. However, MaxMin can be extended to likelihood approaches and is therefore preferable. With respect to rates and timing, all algorithms currently employed for ancestral ecology reconstruction bias toward slow rates of eco-diversification relative to lineage diversification. An alternative to this bias is provided by DIVA, which biases toward accelerated rates of eco-diversification and thus inferences of ecology-driven speciation. We see no way of choosing between these biases; however, phylogeneticists should be aware of them. Applying these methods to Thamnochortus, we find there to be important differences in details, yet general congruence, regarding the historical ecology of this clade. We infer the most recent common ancestor of Thamnochortus to have been a post-fire resprouting species distributed on rocky, well-drained, sandstone-derived soils at lower-middle elevations, in regions of moderate levels of yearly (primarily winter) rainfall. This species would have been distributed in habitats much like those of the southwestern Cape mountains today. Major ecological trends include shifts to lower rainfall regimes and shifts from sandstone to limestone-derived alkaline soils at lower altitudes.

Taxon sampling effects in molecular clock dating: an example from the African Restionaceae.

Three commonly used molecular dating methods for correction of variable rates (non-parametric rate smoothing, penalized likelihood, and Bayesian rate correction) as well as the assumption of a global molecular clock were tested for sensitivity to taxon sampling. The test dataset of 6854 basepairs for 300 terminals includes a nearly complete sample of the Restio-clade of the African Restionaceae (272 of the 288 species), as well as 26 outgroup species. Of this, nested subsets of 35, 51, 80, 120, 150, and the full 300 species were used. Molecular dating experiments with these datasets showed that all methods are sensitive to undersampling, but that this effect is more severe in analyses that use more extreme rate smoothing. Additionally, the undersampling effect is positively related to distance from the calibration node. The combined effect of undersampling and distance from the calibration node resulted in up to threefold differences in the age estimation of nodes from the same dataset with the same calibration point. We suggest that the most suitable methods are penalized likelihood and Bayesian when a global clock assumption has been rejected, as these methods are more successful at finding optimal levels of smoothing to correct for rate heterogeneity, and are less sensitive to undersampling.