Correlates of hyperdiversity in southern African ice plants (Aizoaceae).

The exceptionally high plant diversity of the Greater Cape Floristic Region (GCFR) comprises a combination of ancient lineages and young radiations. A previous phylogenetic study of Aizoaceae subfamily Ruschioideae dated the radiation of this clade of > 1500 species in the GCFR to 3.8-8.7 Mya, establishing it as a flagship example of a diversification event triggered by the onset of a summer-arid climate in the region. However, a more recent analysis found an older age for the Ruschioideae lineage (17 Mya), suggesting that the group may in fact have originated much before the aridification of the region 10-15 Mya. Here, we reassess the tempo of radiation of ice plants by using the most complete generic-level phylogenetic tree for Aizoaceae to date, a revised calibration age and a new dating method. Our estimates of the age of the clade are even younger than initially thought (stem age 1.13-6.49 Mya), supporting the hypothesis that the radiation post-dates the establishment of an arid environment in the GCFR and firmly placing the radiation among the fastest in angiosperms (diversification rate of 4.4 species per million years). We also statistically examine environmental and morphological correlates of richness in ice plants and find that diversity is strongly linked with precipitation, temperature, topographic complexity and the evolution of highly succulent leaves and wide-band tracheids.

Introduction

The Greater Cape Floristic Region (GCFR) of southern Africa is the most biodiverse temperate region of the globe in terms of plant species (Born, Linder & Desmet, 2007; Kreft & Jetz, 2007). Considerable research has been conducted with the aim of disentangling the major forces that have driven diversification in the region (Linder, 2003; Verboom ; Schnitzler ). The current consensus is that high diversity in the GCFR is due to a combination of the gradual accumulation of species from old plant lineages (Linder, 2008; Verboom ; Valente , 2011), and of recent and rapid radiations that are thought to have been triggered by the establishment of a summer-dry climate in the south-western tip of southern Africa in the Miocene (Richardson ; Linder, 2003; Verboom ; Klak, Reeves & Hedderson, 2004).

The most charismatic and perhaps most spectacular of the recent ‘explosive’ radiations in the GCFR is that of subfamily Ruschioideae of Aizoaceae (ice plant family). Represented by 1585 species in 112 genera, Ruschioideae are one of the most species-rich and diverse clades of angiosperms in southern Africa (Smith, 1998; Goldblatt & Manning, 2002). With subfamily Mesembryanthemoideae, species of Ruschioideae are commonly referred to as the ‘mesembs’. Using phylogenetic methods, Klak dated the radiation of a subclade of the ruschioids (core Ruschioideae) at 3.8–8.7 Mya, with a diversification rate of 0.77–1.75 species Myr−1. This recent and fast radiation is considered one of the most rapid recorded in angiosperms (Valente, Savolainen & Vargas, 2010b) and rivals some of the fastest radiations in the world, such as that of cichlid fish (Verheyen ).

Nearly a decade since the publication of the ice plant dating analysis of Klak , little progress has been made with regard to understanding the precise tempo and causes of the radiation of core Ruschioideae. A key question remains whether the radiation pre-dates or post-dates the establishment of a summer-arid climate in the south-western tip of southern Africa in the mid-Miocene (10–15 Mya, Cowling, Procheş & Partridge, 2009). Aridification of the region was associated with the establishment of the Benguela upwelling system and is thought to have led to the extinction of moist-adapted lineages and the opening of new niches (Zachos ; Dupont ), which may have triggered the radiation of mesembs. However, a study by Arakaki suggested that the radiation of Ruschioideae may in fact be much older than previously thought (approximately 17 Mya), challenging the view of Klak that the radiation post-dates the establishment of the semi-arid regime in the GCFR. In their dating analysis, Klak used a limited sampling of the clade (< 50% of the genera), derived their calibration point from an angiosperm-wide chronogram (Wikström, Savolainen & Chase, 2001) that has been often criticized owing to its estimates often being much older than the fossil record (Anderson, Bremer & Friis, 2005) and did not have access to the recent advances in relaxed molecular-clock Bayesian dating methods (Drummond ). Therefore, the question remains as to whether the hypothesis of Klak et al. that the ice plant radiation was triggered by climate change in the GCFR is valid, or whether a more ancient origin, as proposed by Arakaki , is preferred. An older age could also imply that rates of diversification of mesembs are not as spectacular as previously thought (Valente ).

Another pending question is whether morphological innovations may have aided the radiation of the core ruschioids. Klak proposed that two key characters (leaf shape and tracheid cell type) have been linked with evolutionary success of mesembs. First, core ruschioids show great leaf shape diversity, often with succulent trigonous shapes or stones that allow increased storage of water and prevent excess water loss, in contrast to earlier diverging species of Aizoaceae, which tend to have flat, non-succulent leaves (Smith, 1998). Secondly, several ice plant species possess wide-band tracheids (WBTs) that have larger secondary cell walls than normal tracheids and which are thought to be an adaptation to aridity by better withstanding water stress (Mauseth ; Landrum, 2001, 2006, 2008). However, to date, the evolution of these traits has never been examined using phylogenetic ancestral state reconstruction methods and therefore little is known regarding their evolutionary history and whether they may have indeed been linked to diversification, as proposed by Klak .

A final question that remains unresolved is what were the drivers of diversification of ruschioids in southern Africa. Studies in southern African Aizoaceae have suggested a positive link between mesemb diversity and environmental factors associated with aridity (Ihlenfeldt, 1994; Ellis, Weis & Brandon, 2006). Mesembs form a major component of the most arid region of the GCFR (Fig. 1), the Succulent Karoo (SK) eco-region, comprising nearly 20% of its species (Goldblatt & Manning, 2002; Born ). This suggests that there may be a link between the arid climatic conditions of the SK and mesemb species richness, but this has never formally been tested. We also know little about the relationship between mesemb evolutionary success and other key environmental factors that have commonly been thought to have triggered diversification of plant clades in the GCFR, namely topography and soil type (Linder, 2003; Schnitzler ). Both topographic complexity and edaphic diversity have been hypothesized to have played a role in the radiation of mesembs, through the generation of increased opportunities for, respectively, allopatry and ecological divergence (Ellis ).

Figure 1

Map of southern African eco-regions and mesemb genus richness per quarter degree square (QDS).

Here, we reconstruct the largest generic-level phylogenetic tree for Aizoaceae in order to re-examine the evolutionary history of mesembs using the most representative and complete sampling of the clade to date (nearly 90% of the genera), in combination with a revised calibration age and the latest Bayesian relaxed-clock dating methods. Our main aim is to test whether the radiation of core Ruschioideae does indeed post-date the establishment of the summer-arid climate in southern Africa 10–15 Mya as widely accepted and proposed by Klak , or whether a more ancient origin as proposed by Arakaki is favoured. We also examine whether diversification rates in the clade, as derived from our new dating analysis, still rate among the highest recorded in angiosperms. Secondly, we take advantage of our new phylogenetic tree for Aizoaceae to reconstruct the evolution of leaf shape and tracheid type, two morphological traits proposed by Klak to have been linked with the evolutionary success of the core ruschioids. Finally, we use detailed information on geographical patterns of mesemb genus richness to test formally the importance of key environmental factors such as rainfall and temperature in determining ice plant diversity in the SK and the GCFR as a whole.

Material and Methods

Taxon sampling

Aizoaceae are composed of four subfamilies with 135 genera and 1830 species (Smith, 1998; Klak ; Klak, Bruyns & Hedderson, 2007). Two of these subfamilies, Aizooideae and Sesuvioideae, are slightly succulent shrubs with a worldwide distribution of 145 species in 12 genera (Smith, 1998; Klak ). The remaining 1685 species and 123 genera are in two subfamilies, Mesembryanthemoideae and Ruschioideae, and are commonly referred to as the ‘mesembs’ (the succulent members of Aizoaceae; Smith, 1998). Mesembryanthemoideae, the weedy mesembs, are a relatively small group with around 100 species in 11 genera. In contrast, Ruschioideae contains about 1585 species in 112 genera, with the early diverging group consisting of ten genera and 22 species, and the species-rich core group consisting of 1563 species in 102 genera (Klak ; Table 1).

Table 1

Genera and species numbers for the three ‘mesemb’ groups (data modified from Klak , 2007)

Group	No. of genera	No. of species
Mesembryanthemoids	11	100
Early-diverging ruschioids	10	22
Core ruschioids	102	1563
Total	123	1685

Representatives of 106 out of the 123 currently recognized mesemb genera were sampled in this study; the majority were collected in South Africa during 2008 and 2009, and stored in silica gel. To provide as complete a sampling as possible, accessions were taken from GenBank for genera for which we did not have silica-dried material, in addition to outgroup species (Appendix 12009).

DNA sequencing

Total cellular DNA isolation was performed using a modified CTAB procedure (Doyle & Doyle, 1987; Savolainen ). 2× CTAB lysis buffer (500 μL) was added to 50–100 mg of ground plant material. An equal volume of an active protein denaturant (SEVAG; 24 parts chloroform to one part isoamyl alcohol) was added to allow the DNA phase to be isolated. The pellet was precipitated and washed in ethanol, then resuspended in 50 μL TE buffer for freezer storage. Two non-coding plastid gene regions, trnL-F and psbA-trnH, were amplified to expand upon the previous phylogenetic tree produced by Klak , 2004). The trnL intron and trnL-F intergenic spacer were amplified, using primers c, d, e and f (Taberlet, Gielly & Bouvet, 1991), in two separate reactions. The psbA–trnH intergenic spacer, which has been recommended as a putative DNA barcode due to it being one of the most variable non-coding regions of the plastid genome (Kress ; Shaw ), was amplified using primers psbAF and trnHR (Sang, Crawford & Stuessy, 1997).

For each 1 μL of DNA template, 22.5 μL of ReddyMix master mix (ABgene), 4 μL MgCl2 (2.5 mm), 0.5 μL bvine serum albumin (0.4%), and 0.5 μL forward and 0.5 μL reverse primer (10 μm) were added, to make a 29-μL reaction. For all three regions the following PCR conditions were used: an initial denaturation at 94 °C for 3 min to melt the double strands of DNA, followed by 28 cycles of 1 min denaturation at 94 °C, 1 min annealing at 48 °C and an extension of 72 °C for 1 min, with a final additional extension at 72 °C for 10 min. Success of PCR was verified by 1% agarose gel electrophoresis and successful reactants were purified either with the QIAquick PCR Purification Kit (Qiagen), the products being eluted in EB elution buffer (Qiagen), or with ExoSAP-IT [Exonuclease I and Shrimp Alkaline Phosphatase Recombinant (rSAP); USB]. Cycle sequencing reactions were carried out in 10-μL reactions consisting of: 40 ng cleaned amplification product, 0.5 μL BigDye Terminator Cycle Sequencing Ready Reaction kit v3.1 (Applied Biosystems), 0.75 μL primer (0.1 ng μL−1; PCR primers used as sequencing primers), 3.0 μL sequencing buffer prepared according to the manufacturer's instructions, and sterile distilled water to make up a final volume of 10 μL. The cycle sequencing thermal profile consisted of 26 cycles of 10 s denaturation at 96 °C, 5 s annealing at 50 °C and 4 min at 60 °C. Complementary strands were sequenced on an ABI 3130xl automated DNA sequencer (Applied Biosystems), following the manufacturer's protocols.

Phylogenetic inference, divergence time estimation and diversification rates

Complementary strands were edited and assembled in Sequencher v4.5 (Genes Codes Corp.). Sequences were aligned by eye in Geneious Pro v5.6.3, with missing genera obtained from GenBank added to the matrix (Appendix 12009). Both regions were plastid and non-coding and were expected to produce congruent results, and were therefore combined for further analyses. Bayesian analysis was performed in BEAST v1.7.4 (Drummond & Rambaut, 2007). ModelTest (Posada & Crandall, 2010) was used to select the most appropriate model of sequence evolution, based on the lowest Akaike information criterion (AIC) score. Two independent BEAST analyses were run for 20 million generations, sampled every 2000. Tracer v1.5 was used to check the progress of the Bayesian analysis, and TreeAnnotator v1.7.4 obtained the consensus tree excluding the first three million generations (‘burn in’ phase). Posterior probabilities (PPs) were assigned in FigTree v1.3.1.

Due to the absence of fossil data, molecular dating was conducted using the previously estimated divergence time between Aizoaceae/Phytolaccaceae and Nyctaginaceae of 21 Mya (Wikström ; Forest & Chase, 2009). This age was estimated from an angiosperm-wide tree using non-parametric rate smoothing (NPRS) applied to branch lengths obtained using accelerated transformation optimization in maximum parsimony (ACCTRAN). Klak used a calibration age for the same node of 26 Mya, opting for the age obtained based on maximum-likelihood branch lengths from the Wikström tree. However, the analysis of Wikström has often been criticized for disagreeing with the fossil record and overestimating node ages (Anderson ; Forest & Chase, 2009), and we therefore opted for the youngest estimate from their study, which was that based on the ACCTRAN optimization. The prior for the stem age of Aizoaceae/Phytolaccaceae was set using a normal distribution (mean = 21, SD = 1). We repeated the analysis using the same calibration age as Klak , who used an age of 26 Mya for this split (normal distributed prior, mean = 26, SD = 1).

We estimated net rates of diversification for Aizoaceae and the core Ruschioideae clade using the whole-clade estimator of Magallón & Sanderson (2001). Rates were calculated assuming no extinction (E = 0) or a high rate of extinction relative to speciation (E = 0.9) for both crown and stem groups using the R package Geiger (Harmon ). We repeated analyses for the ages obtained using the two alternative calibration ages.

Ancestral trait reconstruction of morphological traits

We conducted character optimization analyses to reconstruct the evolution of two traits that have been hypothesized to have played a role in diversification of mesembs: leaf shape and tracheid type (Klak ). Leaf shape was scored as: cylindrical, flat, trigonal and stone. WBTs were scored as present or absent. Data on leaf shape and tracheid type were extracted from Smith (1998) and Landrum (2001, 2008). In addition, leaf shape data were supplemented with personal communications from P. Burgoyne, based on field observations. Character states are mostly conserved within genera (Smith, 1998), and when that is not the case the characters were scored as polymorphic. The trait data are given in Appendix 22005.

The Bayesian maximum clade credibility tree was used in Mesquite v2.7.2 (Maddison & Maddison, 2009) for unordered parsimony ancestral state reconstruction. To account for uncertainty in tree topology and branch lengths, character optimizations were repeated for each of 1000 trees from the BEAST output.

Environmental analysis

We tested whether four environmental factors (precipitation, temperature, topographic complexity and soil type diversity) are associated with mesemb genus richness. Genus diversity may be decoupled from species diversity if species richness is unevenly distributed among genera. However, genus richness has been shown to provide an excellent proxy for species richness in hyperdiverse plant clades and biodiversity hotspots that are not yet amenable to species-level analyses due to their exceptionally high numbers of species (Villaseñor ; Mazaris ), as is the case of mesembs (> 1500 species).

South African quarter degree squared (QDS) generic distribution data were obtained from the PRECIS database at SANBI Pretoria. The mesemb-specific data were extracted in R v2.12.1 (http://www.r-project.org), and ArcGIS v9.2 (Esri) was used to map genus richness per QDS. ArcGIS was also used to map ecological data using the statistics function to quantify values per QDS. Data on the mean annual precipitation (mm, 2.5-min resolution), mean temperature of the driest quarter (°C × 10, 2.5-min resolution), and topographical complexity (measured as altitudinal standard deviation; m, 2.5-min resolution) were gathered from Bioclim (http://www.worldclim.org/bioclim). Soil type diversity (number of soil types per QDS) was obtained from Schnitzler and based on the SOTER-based soil parameter estimates for southern Africa (version 1.a; Batjes, 2004). To test for environmental predictors of generic richness (of mesembs overall and of the main clades separately: Mesembryanthemoideae, the early diverging ruschioids and the core ruschioids), linear regressions of generic richness per QDS and environmental variables were performed in R.

Results

The total number of species included in the combined dataset was 143, representing 108 genera of Aizoaceae, plus two outgroup taxa; this included 80 species generated for trnL-F and 100 species for psbA–trnH. Sequences for the remaining species were taken from GenBank, resulting in a total of 106 mesemb genera, 89 of which were from the core ruschioid group (Appendix 12009). Topology of the strict consensus tree was inferred with Bayesian methods, using a general time-reversible model with gamma distributed rate variation, and topological results were congruent with previous studies (Klak ; Fig. 2). The analysis shows strong support at subfamily level for Aizooideae, Mesembryanthemoideae and Ruschioideae (PPs = 1); resolution is also strongly supported between the early diverging and core groups of Ruschioideae (PP = 1). Core Ruschioideae are recovered as a monophyletic group (PP = 1), although the lack of resolution within the clade is apparent, and hence the relationships among those genera remain largely unresolved. However, some subclades of core Ruschioideae are strongly supported: for example, Hallianthus H.E.K.Hartmann and Leipoldtia L.Bolus (PP = 0.97), and Bijlia N.E.Br. and Argyroderma N.E.Br. (PP = 0.88). In some cases for which multiple accessions were analysed per genus, the genera were recovered as monophyletic, such as Malephora N.E.Br. (PP = 0.97), Orthopterum L.Bolus (PP = 1), Lithops N.E.Br. (PP = 0.88) and Drosanthemum Schwantes (PP = 0.62). However, in the majority of cases where multiple accessions were sequenced, accessions Acrodon N.E.Br., Delosperma N.E.Br., Faucaria Schwantes, Lampranthus N.E.Br., Oscularia Schwantes and Ruschia Schwantes were distributed across core Ruschioideae, highlighting the lack of resolution or, alternatively, providing evidence for non-monophyly of these genera.

Fig 2

Maximum-clade credibility tree of all 143 taxa based on analysis of the trnL-F and psbA–trnH matrix. Dark branches indicate nodes with support greater than PP = 0.8. Bars to right shows species of: green, outgroup; grey, Aizooideae; light blue, Mesembryanthemoids; dark blue, early-diverging Ruschioideae; red, core Ruschioideae. Mean ages (Mya) obtained for key nodes in the phylogenetic tree are indicated with arrows. The 95% highest posterior density intervals for the node ages are shown in blue horizontal bars at each node. The stem node of the core ruschioid radiation is indicated with a red circle and the crown node with an orange circle. Time scale bar shown at the bottom (Mya).

The dating analysis (Fig. 2, Table 2) using the new calibration age estimated the split between Aizooideae and the mesembs to have occurred 7.88 Mya (3.01–14.93; 95% highest posterior density interval), with Mesembryanthemoideae splitting from Ruschiodeae 6.02 Mya (2.18–11.82). The radiation of the core ruschioids has a stem age of 3.31 Mya (1.13–6.49; split with early diverging Dorotheantheae clade) and a crown age of 1.50 Mya (0.35–3.14; Fig. 2). The analysis using the calibration age of Klak produced slightly older estimates (Table 2). However, both analyses in this study produced younger ages for all nodes when compared with the study of Klak , which dated the age of core ruschioids at 3.8–8.7 Mya.

Table 2

Clade ages (Mya) and net diversification rates (r; species Myr–1) for selected clades

		Aizoaceae	Core ruschioids
	No. of species	1860	1563
Calibration age of 21 Mya	Crown age	7.88 (3.01–14.93)	1.50 (0.35–3.14)
		Stem age	13.75 (5.43–20.56)	3.31 (1.13–6.49)
E = 0	r (crown)	0.86 (0.45–2.27)	4.44 (2.12–19.03)
		r (stem)	0.55 (0.37–1.39)	2.22 (1.13–6.51)
E = 0.9	r (crown)	0.66 (0.35–1.72)	3.34 (1.59–14.30)
		r (stem)	0.38 (0.25–0.96)	1.53 (0.78–4.48)
Calibration age of 26 Mya	Crown age	9.48 (3.31–15.80)	2.02 (0.56–3.76)
		Stem age	15.45 (7.09–24.92)	3.94 (1.15–6.56)
E = 0	r (crown)	0.72 (0.43–2.07)	3.30 (1.77–11.90)
		r (stem)	0.49 (0.30–1.06)	1.87 (1.16–6.40)
E = 0.9	r (crown)	0.55 (0.33–1.56)	2.48 (1.81–8.94)
		r (stem)	0.34 (0.21–0.74)	1.28 (0.80–4.40)

Ages were obtained from the maximum-clade credibility trees of the Bayesian divergence dating analyses using the calibration age of 21 Mya (Wikström ; Forest & Chase, 2009) or, alternatively, 26 Mya, the same calibration age used by Klak . Diversification rates were estimated using the whole-clade method of Magallón & Sanderson (2001), assuming no extinction (E = 0) or high rate of extinction relative to speciation (E = 0.9).

We obtained a net diversification rate for the crown group of core ruschioids of 4.44 (2.12–19.03) sp Myr−1 using the younger calibration age, or 3.30 (1.77–11.90) sp Myr−1 using the older calibration age; the net rate for the core ruschioids has been at least four times faster than the background rate for the family (Table 2).

Ancestral state reconstruction: role of key innovations in ruschioid diversification

Morphological ancestral state reconstruction reveals that flat leaves are the most likely ancestral leaf state (Fig. 3A) for Aizoaceae (100% of the 1000 trees analysed had flat leaf as the state optimized to the root of Aizoaceae). The two species of Aizooideae sampled, Aizoon canariense L. and an unidentified species of Galenia L., have flat leaves, as do all the early-diverging ruschioids, with the exception of Conicosia N.E.Br., which has cylindrical leaves. Just under half of the genera of Mesembryanthemoideae sampled have flat leaves, with the remainder having cylindrical or a combination of flat and cylindrical leaves within a genus. The character optimization analyses support a major shift in succulence of leaf type at the origin of the core ruschioids, with 82.1% of the trees sampled presenting a shift from flat leaves to highly succulent leaves at the most recent common ancestor of core Ruschioideae, whereas 17.8% of the trees presented an equivocal reconstruction of leaf shape at that node, meaning that character state could not be confidently identified in those trees. The majority of ruschioid genera sampled have trigonous leaves. Across the 1000 posterior trees, cylindrical leaves have evolved an average of ten times (minimum eight, maximum 17) in core Ruschioideae, whereas stone leaves have evolved an average of 13 times (minimum nine, maximum 17) in core Ruschioideae. Reversal to non-succulence has never occurred in the radiation.

Fig 3

Character state optimization of (A) leaf shape and (B) tracheid type (WBT, wide-band tracheid). Characters mapped onto the maximum clade credibility tree from the Bayesian analysis.

The ancestral tracheid cell type was recovered as unmodified in 100% of the trees, with no members of Aizooideae, Mesembryanthemoideae or early-diverging Ruschioideae having WBTs, whereas the vast majority (61 out of 68 sampled) of core ruschioids have them (Fig. 3B). The most recent common ancestor of the core ruschioids was optimized to have WBTs in 100% of the trees. According to the posterior distribution of trees, WBTs have been lost an average of eight times in the core ruschioids and are absent in Carpobrotus N.E.Br., Chasmatophyllum Dinter & Schwantes, Conicosia, Conophytum N.E.Br., Erepsia N.E.Br., Gibbaeum Haw., Glottiphyllum Haw. and Rabiea N.E.Br.

The analysis of environmental and genus richness data shows that core ruschioid genus diversity is higher in areas of low precipitation, such as the Little Karoo and Richtersveld Mountains (Fig. 4A, Table 3). The opposite is true of genera of early-diverging ruschioids, which are present in highest numbers in areas of higher precipitation, such as coastal fynbos habitats in the South Western Cape (Fig. 4B). The linear models (Table 3) demonstrate that precipitation, temperature and topographical complexity all have a significant effect (P < 0.0001) on overall mesemb genus richness, whereas soil type does not. Increasing precipitation has a negative effect on genus richness, whereas topographic heterogeneity and increasing temperature both have positive effects. These explanatory variables have a similar effect on core ruschioid genus richness (Table 3), but there are no significant correlates with generic richness of the species-poor early-diverging ruschioids. Precipitation and temperature have significant effects on mesembryanthemoid genus richness but topographical complexity does not.

Fig 4

Map and correlation graphs of mean annual precipitation (mm) against core (A) and early-diverging (B) ruschioid genus richness (per QDS), showing a map of genus richness (size of dots indicates genus richness) against a background of mean annual precipitation.

Table 3

Results of linear models with generic richness as a response variable and environmental parameters as predictor variables for the different groups

	Effect size	Standard error	t	P value
All mesembs
Topographic complexity	0.009	0.001	6.364	< 0.0001
Precipitation	−0.005	0.0006	−8.120	< 0.0001
Dry-season temperature	0.0370	0.003	14.413	< 0.0001
Soil type diversity	−0.060	0.101	−0.594	0.5530
Mesembryanthemoids
Topographic complexity	0.001	0.0009	1.478	0.1401
Precipitation	−0.002	0.0003	−6.204	< 0.0001
Dry-season temperature	0.004	0.001	2.738	0.0065
Soil type diversity	0.052	0.057	0.945	0.3453
Early diverging ruschioids
Topographic complexity	0.001	0.0009	1.393	0.1670
Precipitation	0.0007	0.0004	1.543	0.1260
Dry-season temperature	0.007	0.004	1.593	0.1140
Soil type diversity	0.019	0.0669	0.284	0.7770
Core ruschioids
Topographic complexity	0.006	0.001	6.093	< 0.0001
Precipitation	−0.003	0.0004	−7.718	< 0.0001
Dry-season temperature	0.023	0.0019	12.545	< 0.0001
Soil type diversity	−0.054	0.0735	−0.746	0.4561

Significant P values are given in bold.

Discussion

Evolutionary relationships

This is the largest generic-level phylogenetic tree of Aizoaceae reconstructed to date, including 87% of mesemb genera. It upholds previous molecular studies with regard to good subfamilial-level support (Klak , 2004); however, within the core ruschioids, there is a lack of sequence divergence that is a common occurrence when attempting to reconstruct phylogenies of young radiations (Valente ). Young, closely related taxa often exhibit little variation in their sequences, and this has been found to be the case in many phylogenetic studies of southern African plant groups (e.g. Richardson ; Schnitzler ; Valente ). In Aizoaceae, previous work has shown that the trnL-F and psbA–trnH regions have proved efficient at resolving subfamily-level relationships, but the difficulty of resolving closely related sister taxa remains (Klak ). Recent studies have used amplified fragment length polymorphism (AFLP) loci and microsatellite markers to explain the evolutionary relationships between and within rapidly radiated GCFR clades (e.g. Buys ; Prunier & Holsinger, 2010; Rymer ; Valente ). In mesembs, an AFLP study of Argyroderma by Ellis found low support for the monophyly of individuals of the same species in many cases, again suggesting that the young age of the radiation hampers the ability to identify evolutionary units in core Ruschioideae, even when more sensitive markers are used.

Timing of the radiation

Our divergence time estimates for the core ruschioids (Fig. 2, mean crown age 1.50 Mya; mean stem age 3.31 Mya) are younger than those estimated by Klak et al. (3.8–8.7 Mya; Klak ) and considerably younger than those obtained by Arakaki et al. (17 Mya; Arakaki ). The new dates are clearly consistent with the hypothesis that the radiation of the core ruschioids post-dates the onset of a summer-arid climate in south-western southern Africa 10–15 Mya. Our new dating analysis for Aizoaceae therefore rejects the alternative hypothesis that the radiation is ancient as recently proposed by Arakaki . A young age for ice plants in the GCFR, as favoured by our dating analysis, is also consistent with a recent palynological study that found that Aizoaceae pollen records were absent before 8 Mya (Dupont ).

Coupled with the younger divergence time estimated for the core ruschioids, we find that the diversification rates we estimate for this group (Table 2, younger calibration – 4.44 sp Myr−1; older calibration 3.30 sp Myr−1) are even higher than formerly postulated (0.77–1.75 sp Myr−1; Klak ). These findings uphold the view that this represents one of the most recent and rapid radiations known in angiosperms (Valente ). In fact, the confidence intervals of our estimates for diversification rate of core Ruschioideae surpass those for any other plant group, including Dianthus L. (Caryophyllaceae), the most rapid plant radiation document to date (Valente ), although we acknowledge that limited sequence divergence may bias the dating calculations to a degree.

There are several reasons why our dating estimates are younger than those of Klak . First, we chose to use a younger calibration age for the same node that was used by Klak because the study from which their age was derived (Wikström ) has been criticized for producing ages that disagree with the fossil record by overestimating the ages of nodes (Anderson ). By choosing the minimum estimate from that study (21 versus 26 Mya), as proposed by Forest & Chase (2009), we aimed to reduce the reported bias towards older estimates. However, even when the same calibration age of 26 Mya was used, our estimate for the age of the core ruschioids was also younger than that obtained by Klak ; mean crown age 2.02 Mya; mean stem age 3.94 Mya). In this case, the younger age obtained is therefore potentially due the use of Bayesian divergence dating, which is known to produce younger ages than those estimated by NPRS (Linder, Hardy & Rutschmann, 2005), the dating method used by Klak and colleagues. The Bayesian relaxed-clock divergence dating method implemented in BEAST (Drummond & Rambaut, 2007) is generally thought to produce more accurate ages than NPRS (Sanderson, 1997), because it models molecular rate among lineages as varying in an autocorrelated manner and incorporates phylogenetic uncertainty into the dating process. It therefore does not require as much prior information about rate variation within different clades in the tree (Drummond ). The use of much denser genus sampling, as is the case in our phylogenetic tree, can also affect dating estimates, although it usually does so by increasing rather than decreasing node ages (Linder ).

Our dating analysis lacks a direct fossil calibration, and we did not include large error estimates around the calibration age (other than the standard deviation of the node age normal prior). Therefore, our results must be viewed with caution. The study that found an ancient age of Ruschioideae (Arakaki ) did use fossils to calibrate a wider tree of worldwide succulent clades that included core Ruschioideae, but the fossils were distantly related to Aizoaceae, and the authors used a much sparser sampling of mesembs. Therefore, it is difficult at this stage to assess whether our young estimates are more accurate, as both approaches present drawbacks.

Morphology

The majority of posterior trees used in our analyses of character optimization presented a shift from non-WBT to WBT and from flat leaves to highly succulent leaves at the stem node of core Ruschioideae, suggesting that WBTs and high leaf succulence evolved around the same time in the most recent common ancestor of the radiation (Fig. 3). Ancestral trait reconstruction of leaf shape revealed that all core ruschioids sampled in this study evolved from the ancestral state of slightly succulent flat leaves to highly succulent leaves (Fig. 3A). These highly succulent leaves vary from cylindrical (e.g. Trichodiadema Schwantes) and trigonal (e.g. Erepsia and Faucaria) to compact stones (e.g. Lithops, Conophytum and Argyroderma; in ten out of 11 species; Ellis ). All have a decreased surface area that can account for much lower water loss than flat-leaved ice plants (van Jaarsveld, 1987) and they also act as water storage organs, allowing the core ruschioids to survive in areas of much lower rainfall (Smith, 1998). The most extreme level of succulence seen, namely the miniature succulents with reduced stone-shaped leaves, appears to have independently evolved multiple times and is evident in 15 genera, suggesting convergent evolution of a character state that may be beneficial in semi-arid conditions. Most miniature succulent genera are either monospecific or contain only a few species, with the exceptions of Conophytum (88 species), Lithops (37), Gibbaeum (16) and Argyroderma (11), indicating that this trait may not have been a universal driver of diversification in the core ruschioids (to test this hypothesis an analysis of the effect of the trait on speciation/extinction rates would be required; see below). Species in the early-diverging ruschioid group all have the ancestral state of flat leaves, with the exception of Conocosia, which has cylindrical leaves. Approximately half of the genera in Mesembryanthemoideae have flattened mesomorphic leaves, with the other half having cylindrical leaves (Klak ).

The presence of WBTs with wider secondary walls that prevent collapse under water stress development is purely a core ruschioid adaptation, and it evolved in mesembs at the most recent common ancestor of the core ruschioids (Fig. B). It is likely that WBTs could have evolved in tandem with leaf succulence to improve water storage ability. A few core ruschioid genera appear to have lost WBTs but, as suggested by Landrum (2001), this may be because these structures were not yet developed in the examined specimens. WBTs require high light levels and low water availability to be initiated, and therefore the possibility exists that they are present in these genera, but were not detected. If the evolution of WBTs is tied to diversification, then it could be expected that those genera lacking them (Carpobrotus, Chasmatophyllum, Conicosia, Conophytum, Erepsia, Gibbaeum, Glottiphyllum and Rabiea) would have, on average, fewer species than those with WBTs. Species delimitation is problematic in the core ruschioids due to their recent radiation (Klak, Hedderson & Linder, 2003a; Ellis ), but the current taxonomic status of these seven genera indicates that Conophytum is one of the most diverse genera in Aizoaceae, with 88 species, whereas five of the remaining six genera have < 20 species (Erepsia has 27 species) and thus are not considered species-rich (Smith, 1998). The remainder of the genera in the core ruschioids have WBTs, including some of the most species-rich genera [e.g. Delosperma (c. 163 species), Drosanthemum (120), Lampranthus (> 220) and Ruschia (220); Smith, 1998; Klak ]. Importantly, and contrary to what would be expected, the loss of WBTs is not associated with emigration out of the arid zone, as several genera of core Ruschioideae that lack them are present in the SK.

We were not able to conduct formal tests of a link between traits and diversification rates in mesembs for two main reasons. First, the lack of resolution in our genus-level tree, with low PPs for most nodes within the radiation, leads to many of the genera not being retrieved as monophyletic, thus preventing the use of diversification methods that allow the assignment of unsampled species richness to well-defined tree terminals (e.g. Medusa; Alfaro ). Second, although we sampled the majority of mesemb genera, our species-level sampling was still low (< 10%), and we were therefore not able to use diversification analyses that allow the identification of key innovations by estimating the effect of a character on speciation and extinction rates, such as in the Bisse framework (Maddison, Midford & Otto, 2007). Such methods require high levels of sampling (> 80% of species), the use of characters for which no state occurs in < 10% of species (which does not apply to any of the characters in our study, all of which are rare, e.g. absence of WBTs and flat/cylindrical leaves) and trees with > 300 terminals (Davis, Midford & Maddison, 1998). Nevertheless, such diversification analyses would be the ideal new direction to explore in further studies that aim to detect a link between high rates of diversification and morphological adaptations to an arid environment in ice plants.

Environment

Lack of rainfall is the main defining character of an arid environment (Hopkins & Jones, 1983), but high temperatures are also important, notably in the GCFR. Both these features are significant predictors of both mesembryanthemoid and core ruschioid richness (Table 3), highlighting the link between mesemb success and aridity. The core ruschioids differ from the early-diverging ruschioids in terms of areas of occurrence, with the core group typically occurring in drier regions of southern Africa, particularly in the SK. The significant negative relationship (P < 0.0001) between core ruschioid genus richness and precipitation suggests their morphological adaptations to arid environments (e.g. WBTs and highly succulent leaves) have allowed them to take advantage of niches other plants cannot occupy. Higher genus richness is evident in the arid SK, particularly in the dry Little Karoo and Richtersveld Mountain areas where they have radiated in high numbers (Fig. 4A). On the other hand, we found no relationship between genus richness of the early-diverging ruschioid group and precipitation (Table 3). Early-diverging ruschioids possess flat leaves and unmodified tracheid cells and display habitat preference for the fynbos biome of the south-western Cape (Fig. 4B). They are therefore less well adapted to arid conditions, and may show higher rates of extinction in the more arid SK, which may explain the lack of a relationship with precipitation. The distribution of Mesembryanthemoideae has more in common with that of the core ruschioids, with respect to higher genus richness in areas of lower precipitation (Table 3). This could be a result of the development of more succulent, cylindrical leaves, the conspicuous bladder cells on their leaves or their weedy generalist habit (Smith, 1998; Klak ).

We found a significant relationship between core ruschioid genus richness and increasing topographical complexity but no relationship with edaphic diversity. The high topographical complexity of areas where core ruschioid richness is high, such as the Richtersveld Mountains, could have provided ecological opportunities that are known to drive diversification in the absence of novel traits (Hughes & Eastwood, 2006). The rugged mountainous quartz fields found in the SK harbour many endemic plant species, with fine-scale discrimination of species between patches (Ellis & Weis, 2006; Ellis ). The topographical complexity of the quartz habitat seems to have selected for habitat-specific, short-lived drought-resistant flowering stones (Cowling ; Schmiedel & Jürgens, 1999; Ellis ), and hence core ruschioid species occur there in high abundance. However, topography alone cannot explain why core ruschioid richness is higher, given that the south-western Cape, where the early-diverging ruschioids are most diverse, is also topographically complex. It appears that topography may only be positively associated with genus-richness when interacting with other factors, as was previously found in other groups [e.g. Lupinus L. (Fabaceae), Hughes & Eastwood, 2006]. We hypothesize that the reflective ability of quartz in the mountainous quartz fields of the SK could increase levels of UV radiation, which has been known to increase mutation rates. The combination of high opportunities for allopatry associated with topographical complexity and the potentially higher rates of mutation in the quartz habitat could have driven high speciation rates (Rozema ; Rothschild, 1999), an area which warrants further study.

Conclusions

Our new dating analysis of Aizoaceae using denser sampling, a revised calibration age and more powerful dating methods corroborates the long-held hypothesis that the radiation post-dates the establishment of the contemporary summer-arid climate in the GCFR, and therefore rejects the more recent hypothesis that core Ruschioideae are an ancient clade that was already present before the significant aridification process of the GCFR. The new dating analysis also clearly places the radiation of mesembs among the most rapid angiosperm diversification events documented to date. In addition, this study provides strong new evidence suggesting a link between mesemb evolutionary success and arid conditions. We showed that two morphological adaptations to aridity evolved at the origin of the core ruschioid clade, and our results corroborate the hypothesis that high ice plant diversity is associated with arid conditions such as low precipitation and high temperatures. In addition, we also found a significant link between topographical complexity and mesemb diversity, suggesting a role of allopatry in promoting reproductive isolation in this rapidly radiating clade. The morphological innovations that have evolved in members of core Ruschioideae were likely to be crucial for their survival in an arid environment while other lineages became extinct. Whether these adaptations drove speciation as key innovations, at the same time resulting in low extinction rates, cannot be confirmed with the current data, but they are certainly linked to extant patterns of generic richness.

Species	Voucher	Source	GenBank accession numbers
trnL-F	psbA–trnH
Aloinopsis spathulata	PMB10422	This study	KC834485	KC834404
Acrodon bellidiflorus	–	This study	–	AM230592.1
Acrodon purpurea	PMB9850	This study	–	KC834403
Aizoon canariense	Goldblatt & Manning 11 708 (MO)	GenBank	AJ558042.1	–
Amphibolia laevis	PMB10389	This study	KC834486	KC834405
Antegibbaeum fissoides	PMB10721	This study	KC834487	KC834406
Antimima crassifolia	PMB10361b	This study	KC834488	KC834407
Antimima ventricosa	–	GenBank	AJ439015.1	AJ532896.1
Apatesia helianthoides	Klak 800 (BOL)	GenBank	AJ558064.1	–
Apatesia sabulosa	PMB12063	This study	KC834489	KC834408
Aptenia geniculiflora	PMB8859	This study	KC834490	KC834409
Argyroderma pearsonii	PMB10387	This study	KC834491	KC834410
Aridaria brevicarpa	Bruyns 9469 (BOL)	GenBank	AM161375.1	–
Aridaria vespertina	PMB10359	This study	KC834492	KC834411
Aspazoma amplectens	PMB10338	This study	KC834493	KC834412
Astridia longifolia	PMB33	This study	KC834494	KC834413
Bergeranthus scapiger	PMB9123	This study	KC834495	KC834414
Bijlia tugwelliae	Bruyns 2762 (BOL)	GenBank	AJ558093.1	AJ532874.1
Braunsia geminata	–	GenBank	AJ439018.1	AJ532884.1
Braunsia stayneri	PMB10439	This study	–	KC834415
Brianhuntleya intrusa	PMB1333	This study	KC834496	KC834416
Brownanthus neglectus	–	GenBank	AY993973.1	AY996734.1
Brownanthus pubescens	–	GenBank	AJ438998.1	AY996737.1
Carpanthea pomeridiana	Klak 801 (BOL)	GenBank	AJ558065.1	–
Carpobrotus deliciosus	PMB9139	This study	KC834497	KC834417
Carpobrotus muirii	–	GenBank	AJ439021.1	AJ532885.1
Carruanthus ringens	Bruyns 8173a (BOL)	GenBank	AJ558094.1	AJ532872.1
Caryotophora skiatophytoides	Klak 805 (BOL)	GenBank	AJ558066.1	–
Caryotophora sp.	none	This study	KC834498	–
Caulipsolon rapaceum	Klak 750 (BOL)	GenBank	AJ558053.1	–
Cephalophyllum pillansii	Klak 785 (BOL)	GenBank	AJ558100.1	AJ532895.1
Cephalophyllum sp.	PMB11815	This study	KC834499	–
Cerochlamys pachyphylla	PMB19455	This study	KC834500	KC834418
Chasmatophyllum musculinum	PMB11398	This study	KC834501	KC834419
Cheiridopsis turbita	PMB10325	This study	KC834502	KC834420
Cleretum papulosum	Bruyns 8825a (BOL)	GenBank	AJ558070.1	–
Cleretum sp.	PMB11728	This study	KC834503	KC834421
Conicosia elongata	PMB11723	This study	KC834504	KC834422
Conophytum bilobum	PMB10357	This study	KC834505	KC834423
Conophytum bruynsii	Bruyns 6784 (BOL)	GenBank	AJ558090.1	AJ532869.1
Cylindrophyllum comptonii	PMB10429	This study	–	KC834424
Dactylopsis digitata	–	This study	–	AY996740.1
Deilanthe peersii	PMB9234	This study	KC834506	KC834425
Delosperma echinatum	–	GenBank	AJ439001.1	AJ532848.1
Delosperma esterhuyensiae	–	GenBank	AJ439002.1	AJ532849.1
Dicrocaulon sp.	PMB11983	This study	KC834507	KC834426
Didymaotus lapidiformis	PMB11507	This study	KC834508	KC834427
Dinteranthus puberulus	PMB7227	This study	KC834509	KC834428
Diplosoma retroversum	Klak 835 (BOL)	GenBank	AJ558071.1	AJ532845.1
Disphyma dunsdonii	Klak 808 (BOL)	GenBank	AJ558072.1	AJ532846.1
Dorotheanthus bellidiformis	–	GenBank	AJ439000.1	AJ532843.1
Dracophilus dealbatus	PMB20	This study	–	KC834429
Drosanthemum schoenlandianum	–	GenBank	AJ439003.1	AJ532852.1
Drosanthemum sp.	PMB11917	This study	KC834510	KC834430
Eberlanzia dichotoma	–	GenBank	AJ439014.1	AJ532889.1
Eberlanzia gravida	PMB10325	This study	KC834511	KC834431
Ebracteola fulleri	PMB11572	This study	KC834512	KC834432
Enarganthe octonaria	PMB10358	This study	KC834513	KC834433
Erepsia inclaudens	PMB10399b	This study	KC834514	KC834434
Erepsia pillansii	–	GenBank	AJ439027.1	–
Esterhuysenia drepanophylla	–	GenBank	AJ439028.1	–
Faucaria britteniae	PMB8935	This study	KC834515	KC834435
Faucaria felina	Klak 338 (BOL)	GenBank	AJ558085.1	AJ532864.1
Fenestraria rhopalophylla	PMB7371	This study	KC834516	KC834436
Frithia pulchra	PMB1	This study	KC834517	–
Galenia sp.	PMB8898	This study	–	KC834437
Gibbaeum pachypodium	Klak 380 (BOL)	GenBank	AJ558082.1	AJ532859.1
Gibbaeum pubescens	PMB10452	This study	KC834518	KC834438
Glottiphyllum depressum	PMB10435	This study	KC834519	KC834439
Hallianthus planus	PMB7375	This study	KC834520	KC834440
Hammeria sp.	PMB11527	This study	KC834521	KC834441
Hartmanthus sp.	PMB16	This study	KC834522	KC834442
Hereroa crassa	PMB10400	This study	–	KC834443
Hymenogene conica	Klak 802 (BOL)	GenBank	AJ558068.1	–
Jacobsenia sp.	PMB11963	This study	KC834523	–
Jensenobotrya lossowiana	PMB215	This study	KC834524	KC834444
Juttadinteria simpsonii	–	GenBank	AJ439009.1	KC834445
Khadia carolinensis	PMB4542	This study	KC834525	KC834446
Lampranthus amphibolius	–	GenBank	AJ439045.1	AJ532878.1
Lampranthus bicolor	–	GenBank	AJ439042.1	AJ532876.1
Leipoldtia weigangiana	PMB10350	This study	KC834526	KC834447
Lithops julii	–	GenBank	AJ439007.1	AJ532866.1
Lithops sp.	PMB11542	This study	KC834527	KC834448
Machairophyllum albidum	Klak 182 (BOL)	GenBank	AJ558096.1	AJ532875.1
Machairophyllum sp.	PMB8485	This study	KC834528	KC834449
Malephora crassa	PMB11525	This study	–	KC834450
Malephora lutea	Klak 664 (BOL)	GenBank	AJ558083.1	AJ532860.1
Mesembryanthemum crystallinum	PMB11964	This study	KC834529	KC834451
Mestoklema sp.	PMB11863	This study	KC834530	KC834452
Meyerophytum meyeri	PMB10332	This study	KC834531	KC834453
Mitrophyllum grande	PMB10344	This study	–	KC834454
Monilaria moniliformis	Klak787 (BOL)	GenBank	AJ558074.1	AJ532844.1
Monilaria pissiformis	PMB10386	This study	KC834532	KC834455
Mossia intervallaris	PMB8890	This study	–	KC834456
Namaquanthus vanheerdei		GenBank	AJ439049.1	AJ532879.1
Namibia sp.	PMB8480b	This study	KC834533	KC834457
Nananthus aloides	PMB10494	This study	KC834534	KC834458
Neea floribunda		GenBank	FJ039169.1	FJ039025.2
Nelia pillansii	Klak777(BOL)	GenBank	AJ558092.1	AJ532871.1
Nelia schlechteri	PMB10340	This study	KC834535	KC834459
Neohenricia sibbettii	PMB11358	This study	KC834536	KC834460
Neohenricia spiculata	Bruyns 7289 (BOL)	GenBank	AJ558087.1	AJ532862.1
Octopoma sp.	PMB11795	This study	KC834537	KC834461
Odontophorus marlothii	Klak862 (BOL)	GenBank	AJ558101.1	AJ532898.1
Odontophorus nanus	PMB10315b	This study	KC834538	KC834462
Oophytum nanum	PMB10387b	This study	–	KC834463
Orthopterum coeganum	Klak350 (BOL)	GenBank	AJ558088.1	AJ532865.1
Orthopterum waltoniae	PMB8936	This study	KC834539	KC834464
Oscularia comptonii	PMB1139	This study	KC834540	KC834465
Oscularia deltoides	–	GenBank	AJ439004.1	AJ532861.1
Phiambolia unca	PMB7850	This study	KC834541	KC834466
Phyllobolus deciduus	PMB11804	This study	KC834542	KC834467
Phytolacca sp.		GenBank	AJ558037.1	DQ006209.1
Pleiospilos simulans	Klak4988 (BOL)	GenBank	AJ558102.1	AJ532897.1
Pleiospilos sp.	PMB3736	This study	KC834543	KC834468
Polymita steenbokensis	Bruyns 8267 (BOL)	GenBank	AJ558097.1	AJ532893.1
Prepodesma orpenii	PMB10254	This study	KC834544	KC834469
Psammophora modesta	PMB8244	This study	KC834545	KC834470
Psilocaulon dinteri	Bruyns 9511 (BOL)	GenBank	AM161435.1	–
Psilocaulon parviflorum	Klak699 (BOL)	GenBank	AJ558062.1	AYg996741.1
Rabiea albinota	PMB8553	This study	KC834546	KC834471
Rhinephyllum sp.	PMB11485	This study	KC834547	KC834472
Ruschia crassa	PMB8208	This study	KC834548	KC834473
Ruschia strubeniae	Klak318 (BOL)	GenBank	AJ558099.1	AJ532892.1
Saphesia flaccida	Klak799 (BOL)	GenBank	AJ558069.1	–
Schlechteranthus hallii	–	This study	–	AM230589.1
Schlechteranthus maximilliani	PMB10364	This study	KC834549	KC834474
Schwantesia herrei	PMB10297	This study	KC834550	KC834475
Scopelogena bruynseii	–	GenBank	AJ439050.1	AJ532882.1
Sesuvium sesuvioides	Bruyns8876(BOL)	GenBank	AJ558038.1	–
Skiatophytum tripolium	Klak 1030 (BOL)	GenBank	AM161451.1	–
Smicrostigma viride	–	GenBank	AJ439051.1	AJ532881.1
Stoeberia beetzii	PMB11906	This study	KC834551	KC834476
Stomatium alboroseum	PMB10409	This study	KC834552	KC834477
Synaptophyllum juttae	PMB8481	This study	KC834553	KC834478
Tanquana prismatica	PMB10401	This study	–	KC834479
Titanopsis calcarea	none	This study	KC834554	KC834480
Titanopsis hugo schlechteri	–	GenBank	AJ439008.1	AJ532867.1
Trichodiadema emarginatum	Klak817 (BOL)	GenBank	AJ558084.1	AJ532851.1
Trichodiadema rogerii	PMB10434	This study	KC834554	KC834481
Vanzijlia annulata	PMB1390	This study	KC834556	KC834482
Vlokia ater	–	GenBank	AJ439052.1	–
Wooleya farinosa	PMB11919	This study	KC834557	KC834483
Zeuktophyllum sp.	PMB6151	This study	KC834558	KC834484
Zeuktophyllum suppositum	–	GenBank	AJ439054.1	–

Species	Wide-band tracheids: present (0); absent (1); or unknown (?)	Leaf shape: 0 = flat; 1 = cylindrical; 2 = trigonous; 3 = stone; 4 = flat/cylindrical; 5 = trigonous/cylindrical
Acrodon bellidiflorus	0	2
Acrodon purpurea	0	2
Aizoon canariense	1	0
Aloinopsis spathulata	0	2
Amphibolia leavis	?	5
Antegibbaeum fissoides	0	5
Antimima crassa	0	2
Antimima ventricosa	0	2
Apatesia helianthoides	?	0
Apatesia sabulosa	?	0
Aptenia genicula	1	0
Argyroderma pearsonii	0	3
Aridaria brevicarpa	1	1
Aridaria vespertina	1	1
Aspazoma amplectens	?	1
Astridia longifolia	0	5
Bergeranthus scapiger	0	2
Bijlia tugwelliae	0	2
Braunsia geminata	0	2
Braunsia stayneri	0	2
Brianhuntleya intrusa	?	?
Brownanthus neglectus	?	1
Brownanthus pubescens	?	1
Carpanthea pomeridiana	1	0
Carpobrotus nuirii	1	2
Carpobrotus deliciosus	?	?
Carruanthus ringens	0	2
Caryotophora skiatophytoides	?	0
Caryotophora sp.	?	0
Caulipsolon rapaceum	?	1
Cephalophyllum pillansii	?	?
Cephalophyllum sp.	?	2
Cerochlamys pachyphylla	?	2
Chasmatophyllum musculinum	1	2
Cheiridopsis turbita	0	2
Cleretum papulosum	1	0
Cleretum sp.	1	0
Conicosia elongata	1	5
Conophytum bilobum	1	3
Conophytum bruynsii	1	3
Cylindrophyllum comptonii	0	1
Dactylopsis digitata	?	4
Deilanthe peersii	?	2
Delosperma echinatum	0	5
Delosperma esterhuyensiae	0	5
Dicraulon sp.	0	2
Didymaotus lapidiformis	0	3
Dinteranthus puberulus	0	3
Diplosoma retroversum	?	3
Disphyma dunsdonii	0	5
Dorotheanthus bellidiformis	1	0
Dracophilus dealbatus	0	2
Drosanthemum schoenlandianum	0	5
Drosanthemum sp.	0	5
Eberlanzia dichotoma	0	1
Eberlanzia gravida	0	1
Ebracteola fulleri	0	2
Enarganthe octonaria	?	2
Erepsia inclaudens	1	2
Erepsia pillansii	1	2
Esterhuysenia drepanophylla	?	2
Faucaria britteniae	0	2
Faucaria felina	0	2
Fenestraria rhopalophylla	0	3
Frithia pulchra	0	3
Galenia sp.	1	0
Gibbaeum pachypodium	1	3
Gibbaeum pubescens	1	3
Glottiphyllum depressum	1	2
Hallianthus planus	0	2
Hammeria sp.	?	1
Hartmanthus sp.	?	2
Hereroa crassa	0	1
Hymenogene conica	?	0
Jensenobotrya lossowiana	0	3
Juttadinteria simpsonii	0	5
Khadia carolinensis	?	5
Lampranthus amphibolius	0	5
Lampranthus bicolor	0	5
Leipoldtia weigangiana	0	2
Lithops julii	0	3
Lithops sp.	0	3
Machairophyllum albidum	0	2
Machairophyllum sp.	0	2
Malephora crassa	0	5
Malephora lutea	0	5
Mesembryanthemum crystallinum	1	4
Mestoklema sp.	0	5
Meyerophytum meyeri	0	1
Mitrophyllum grande	0	1
Monilaria moniliformis	0	1
Monilaria pissiformis	0	1
Mossia intervallaris	0	2
Namaquanthus vanheerdei	0	1
Namibia sp.	0	3
Nananthus aloides	0	2
Neea floribunda	?	?
Nelia pillansii	?	2
Nelia schlechteri	?	2
Neohenricia sibbettii	0	2
Neohenricia spiculata	0	2
Octopoma sp.	0	2
Odontophorus marlothii	0	2
Odontophorus nanus	0	2
Oophytum nanum	?	3
Orthopterum coeganum	0	2
Orthopterum waltoniae	0	2
Oscularia comptonii	0	2
Oscularia deltoides	0	2
Phiambolia unca	?	?
Phyllobolus deciduus	?	4
Phytolacca sp.	?	?
Pleiospilos simulans	0	3
Pleiospilos sp.	0	3
Polymita steenbokensis	?	2
Prepodesma orpenii	0	2
Psammophora modesta	0	5
Psilocaulon dinteri	?	1
Psilocaulon parviflorum	?	1
Rabiea albinota	1	2
Rhinephyllum sp.	0	5
Ruschia crassa	0	2
Ruschia strubeniae	0	2
Saphesia flaccida	?	0
Schlechteranthus hallii	?	2
Schlechteranthus maximilliani	?	2
Schwantesia herrei	0	2
Scopelogena bruynseii	?	5
Skiatophytum tripolium	?	0
Smicrostigma viride	0	2
Stoberia beetzi	0	2
Stomatium alboroseum	0	2
Synaptophyllum juttae	?	0
Tanquana prismatica	0	3
Titanopsis calcarea	0	2
Titanopsis hugo-schlechteri	0	2
Trichodiadema emarginatum	0	1
Trichodiadema rogerii	0	1
Vanzijlia annulata	0	5
Vlokia ater	?	3
Wooleya farinosa	?	5
Zeuktophyllum sp.	?	2
Zeuktophyllum suppositum	?	2