- Research article
- Open Access
Have giant lobelias evolved several times independently? Life form shifts and historical biogeography of the cosmopolitan and highly diverse subfamily Lobelioideae (Campanulaceae)
BMC Biology volume 7, Article number: 82 (2009)
The tendency of animals and plants to independently develop similar features under similar evolutionary pressures - convergence - is a widespread phenomenon in nature. In plants, convergence has been suggested to explain the striking similarity in life form between the giant lobelioids (Campanulaceae, the bellflower family) of Africa and the Hawaiian Islands. Under this assumption these plants would have developed the giant habit from herbaceous ancestors independently, in much the same way as has been suggested for the giant senecios of Africa and the silversword alliance of Hawaii.
Phylogenetic analyses based on plastid (rbcL, trnL-F) and nuclear (internal transcribed spacer [ITS]) DNA sequences for 101 species in subfamily Lobelioideae demonstrate that the large lobelioids from eastern Africa the Hawaiian Islands, and also South America, French Polynesia and southeast Asia, form a strongly supported monophyletic group. Ancestral state reconstructions of life form and distribution, taking into account phylogenetic uncertainty, indicate their descent from a woody ancestor that was probably confined to Africa. Molecular dating analyses using Penalized Likelihood and Bayesian relaxed clock approaches, and combining multiple calibration points, estimate their first diversification at ~25-33 million years ago (Ma), shortly followed by several long-distance dispersal events that resulted in the current pantropical distribution.
These results confidently show that lobelioid species, commonly called 'giant', are very closely related and have not developed their giant form from herbaceous ancestors independently. This study, which includes the hitherto largest taxon sampling for subfamily Lobelioideae, highlights the need for a broad phylogenetic framework for testing assumptions about morphological development in general, and convergent evolution in particular.
When the great 18th-century naturalist Carl Linnaeus took a closer look at a the whale during the preparation of his natural classification system, he realized for the first time that it was not a fish but a mammal . Its form and fins were certainly like the fish, but it also had mammary glands and lungs and should, thus, be classified as a mammal. Ever since Linnaeus, what we today refer to as convergent evolution (see, for example, [2, 3]) has been advocated as a common and widespread phenomenon in nature. From single features (such as the wings of birds and insects or the cones of cycads and conifers) to whole organisms (shrimps and krill, cacti and euphorbias), superficial similarity has been repeatedly demonstrated to evolve independently in distantly related evolutionary lineages.
In plants, the giant lobelias from the Hawaiian Islands and tropical Africa have been cited as remarkable examples of morphological convergence in the family Campanulaceae (the bellflower family). According to some earlier authors [4–6], convergence from herbaceous plants into tall treelets would have occurred independently in different mountain systems in response to similar tropical alpine climates consisting of nightly frosts and rapid temperature fluctuations. The idea that many morphological features in giant lobelias represent adaptations to proximal environmental factors  is based on the observation that the leaf rosettes of many giant species provide insulation for the central axis, protecting the shoot apex and the young leaves and inflorescences from extreme temperatures. According to Hedberg , a similar evolutionary pattern can be observed between the giant senecios of Africa and the silversword alliance of the Hawaiian Islands (both in Asteraceae), which Hedberg cited as an example of convergent evolution alongside the giant lobelias (Figure 1).
Since these pioneering works, many studies have addressed specific aspects of giant lobelioids, not only providing novel insights into the conspicuous habit of the giant forms  but also their detailed morphology [8–10], chromosome numbers [11, 12], conservation status  and physiology [14–17]. Based on a phylogenetic analysis of 17 species in subfamily Lobelioideae (including 13 giant forms and four herbaceous taxa), Knox et al.  proposed that the giant lobelioids consisted of a Chilean hexaploid group and a pantropical tetraploid group, all derived from herbaceous ancestors. Later, using the Chilean giant species as outgroup and analysing the relationships and biogeographic history among 21 African and one Brazilian giant lobelia, Knox and Palmer  suggested that the eastern African giant lobelias arrived in Africa from the Asian/Pacific region. More recently, Givnish et al.  used a larger data set (including 38 lobelioid species, of which two are herbaceous) to demonstrate that all Hawaiian lobelioids constitute a monophyletic clade, thus corroborating previous results by Givnish and co-workers [8, 19]
Despite these major advances, we still need to verify whether the biogeographic and phylogenetic patterns obtained so far for Lobelioideae will stand the inclusion of increased taxon sampling. In order to confidently test whether the giant forms worldwide evolved independently from herbaceous ancestors, it seems vital to include a considerably larger representation of herbaceous taxa and to root the Lobelioideae tree with representative taxa from outside the subfamily.
Here I build on previous studies to reconstruct the life form and biogeographic evolution of Lobelioideae. I focus on the following questions: (i) have giant lobelioids evolved from herbaceous ancestors several times independently?; (ii) when and where did the giant habit evolve?; and (iii) is the life form and geographic distribution phylogenetically conservative? To address these questions, I generated fast-evolving nuclear sequence data for several giant species and, in combination with all available DNA sequences for selected markers, I performed phylogenetic, biogeographic, character evolution and molecular dating analyses. The results reported here demonstrate that the giant lobelioids have a very different evolutionary history than the giant Asteraceae.
Results and discussion
The maximum likelihood trees based on all Campanulaceae sequences available at GenBank and complemented by novel sequences generated in this study for trnL-F, ITS and rbcL are shown in Figure 2a, c. Lobelioideae is retrieved as monophyletic in all three analyses of each marker separately, providing an indication of the choice of taxa for the subsequent analyses. The results from the Bayesian and bootstrap analyses on the Lobelioideae data set are shown together in Figure 3. Generally, strongly supported relationships were in agreement with, or at least did not contradict, previous phylogenetic analyses [7, 8, 10, 18–21]. The few differences obtained here (especially in clade N4, Figure 3) are discussed below.
The results from the molecular dating analyses are summarized in Table 1. Generally, node age estimations exhibited relatively large confidence intervals, both under the Penalized Likelihood (Figure 4) and the Bayesian relaxed clock (Figure 5) methods. The tree with highest sum of clade credibilities generated under the Bayesian evolutionary analysis by sampling trees (BEAST) analysis (Figure 5; effective sample size = 322) did not exhibit any strongly supported phylogenetic conflict with the consensus cladogram of the MrBayes analysis (Figure 3). Although there was great variation between the age estimates obtained in Penalized Likelihood as compared to BEAST, in all key nodes listed in Table 1 there was considerable overlap of the 95% confidence intervals and highest probability densities. The mean covariance between parental and child lineages in the BEAST analysis was 0.157 (Figure 6a). Since a covariance close to zero indicates no significant autocorrelation of rates (the fundamental assumption made by Penalized Likelihood), this indicates that the BEAST results reported here are probably more realistic than the ages obtained in Penalized Likelihood . For comparison, both results are reported throughout the paper.
Life form shifts
The character state reconstruction of life form is shown in Figure 7 (see Table 2 for the definitions). The reconstruction of life form on the Bayesian cladogram required 68 steps and the character is phylogenetically conservative (P < 0.001; Figure 8a, b). The ancestral life form of the most recent common ancestor (MRCA) of Lobelioideae (N1 in Figure 3) is ambiguous, due to an early split between a chamaephyte lineage (leading to the clade Lobelia physaloides - L. tomentosa) and a hemicryptophyte lineage (the ancestor to all the others).
There is great life form variation in the Lobelioideae at the intergeneric, intrageneric and intraspecific levels (Figure 7, Table 3). Only in a few cases is there a homogeneous life form correlating with a certain clade, as indicated by the low Consistency Index of life form overall (0.27 on the tree shown in Figure 7). One example of local high consistency is in the clade Lobelia oligophylla - L. arenaria (N11 in Figure 3), in which all exhibit the hemicryptophyte habit, at least facultatively. Although the overall low consistency is partly biased by an uneven taxon sampling (for example, there are only two species of Lysipomia sampled but all 30 species are strictly chamaephytes), it does reflect a general pattern of difficulties in identifying natural taxa in Lobelioideae based on morphology (see, for example, [21, 23]).
Evolution of the giant lobelioids
The nanophanerophyte habit (Table 2) developed early in the history of the Lobelioideae, most probably several times (at nodes N4, N10 and N12 in Figure 3; with 81.6%, 69.1% and 97.2% of all reconstructions with uniquely best states, respectively), or, less likely, just once (N3; 45.7% of all reconstructions; Figure 7, Table 3). Present day lineages inhabiting the Hawaiian Islands, French Polynesia, southeast Asia and eastern Brazil (N4 in Figure 3) all derive from a single ancestor (Bayesian posterior probability, Bpp = 1.00, Bootstrap support, Bs = 98). Ancestral state reconstructions indicate that this ancestor was most likely confined to Africa and that it was a nanophanerophyte (in 950 and 816 of 1000 reconstructions, respectively; the other reconstructions being ambiguous rather than a different state; see Figures 7 and 9). Indeed, it is in this clade that the truly giant habit occurs, as exhibited by L. gloria-montis and L. rhynchopetalum (Figure 1). This corroborates previous results, based on fewer species, that giant lobelias are closely related and are ultimately derived from herbaceous ancestors [7, 18, 21] but contradicts some earlier suggestions the giant habit would be plesiomorphic in the family .
Nanophanerophyte ancestors (with buds below 3 m high; see Table 2) have given rise to phanerophytes (buds above 3 m) several times independently. Although several species are strictly phanerophytic (for example, L. giberroa, L. petiolata, L. stricklandiae, Delissea undulata and Cyanea leptostegia), others often vary in height (for example, Sclerotheca jayorum, Brighamia insignis, Trematolobelia macrostachys, Cyanea angustifolia and Siphocampylus giganteus). These results indicate that, in the Lobelioideae, height seems to be a more labile state than woodiness.
Geographic range development of the Lobelioideae is depicted in Figure 9. As with life form, distribution is phylogenetically conservative (P < 0.001; Figure 8c-d), but even more so given its higher consistency index (0.40) and fewer parsimony steps required for reconstructing the Bayesian cladogram (25).
The MRCA of the Lobelioideae, as well as those of most of the early diverging splits (N1 - N3 in Figure 3), are all reconstructed to Africa with confidence. Range shifts are unambiguously inferred from: Africa to Oceania (Lobelia physaloides); Africa to the Neotropics (Lobelia aquatica); the Neotropics to Oceania (Lobelia macrodon - L. arenaria); and Oceania to southeast Asia (Lobelia chinensis). That the Brazilian giant lobelias are together sister to an African species is strongly supported, as previously suggested, by Knox et al. .
There is no support here for the hypothesis that the giant lobelias arrived in eastern Africa from Asia or the Pacific region . However, the resolution at the base of clade N4 is very poor and the relationships among the strongly supported clades within it - and possibly their biogeographic reconstructions - may substantially change with the addition of more species and sequence data. Indeed, in the analysis by Givnish et al. , L. nicotianifolia (a southeast Asian species) was inferred to be sister to a clade comprising French Polynesian, African and Hawaiian subclades. In the present analysis, that 'basal' position is occupied by African species (L. columnaris, L. gibberoa and L. gregoriana), although the support for this placement is weak (Bpp = 0.81, Bs < 50). As taxon sampling, sequence regions and phylogenetic methods have all varied considerably between different studies, further investigation is clearly needed in order to reconstruct a solid biogeographic scenario for the large lobelioids.
Dispersal versus vicariance
Although the exact timing for transcontinental range shifts inferred here are prone to large error margins associated with topological and branch length uncertainties (see the credibility bars at nodes in Figures 4 and 5), all range shifts are estimated to have occurred in the last 50 Ma and, in several cases, presumably much more recently. Given that the final break-up of Gondwanaland took place around ~100 Ma , these results provide evidence that long-distance dispersals must have played a major role in shaping the present-day distribution of taxa, a scenario also corroborated by the molecular dating analysis of Givnish et al. . This contrasts to early suggestions  that distribution patterns may be the result of vicariance caused by continental drift.
Long-distance dispersal seems plausible for lobelioids, since most species produce large amounts of minuscule seeds. From herbarium specimens, I calculated that 1 g of seed of the Hawaiian endemic Clermontia kakeana contains about 36,000 seeds (Figure 10). Although seeds of Clermontia and some other Hawaiian genera (Cyanea and Delissea) are today contained in bird-dispersed berries, other Hawaiian taxa (Lobelia, Trematolobelia and Brighamia) still possess wind-dispersed capsules, a condition also inferred for the MRCA of the Hawaiian clade by Givnish et al. . It is therefore reasonable to conclude that wind dispersal of tiny seeds was the dispersal mode for the ancestor of the Hawaiian lobelioids, as suggested by Givnish et al.  and that, presumably, the same applies for other lobelioid dispersal events. Millions of tons of dust are transported annually from the Sahara into South America , so it is conceivable that, given enough time, lobelioid seeds could also be carried over the 15,000 km separating the Hawaiian Islands and Africa and, subsequently, succeed in germinating and establishing a new population.
Colonization of the Hawaiian Islands
The earliest diversification of lineages endemic to the Hawaiian Islands corresponds to the crown age of node N6 retrieved in the MrBayes analysis (Figure 3). In the PL analysis, that clade had a mean of 22.4 Ma (16.9 - 29.0 Ma, 95% confidence intervals, CI), whereas in the BEAST analysis the mean was 12.2 Ma (5.74 - 16.2 Ma, 95% Highest Posterior Density, HPD; Table 1). However, in the BEAST analysis, where the Hawaiian taxa had to be constrained as monophyletic (see Methods), the crown group of the whole Hawaiian clade had a mean of 16.4 Ma (9.43 - 23.0 Ma, 95% HPD; see Figure 6b). These ages are only slightly older than the estimates by Givnish et al. , which varied between 13.6 ± 3.11 Ma and 13.0 ± 1.00 depending on calibration methodology.
The relatively old node ages estimated here may seem surprising, given that the oldest of the modern Hawaiian Islands with well-developed vegetation only dates back to 5.1 Ma . However, they are fully explainable when considering that these islands are part of a much older archipelago formed by the movement of the Pacific plate over a fixed hot spot. A continuous chain of islands have been elevated since 29.8 Ma  and lobelioid lineages could have continuously colonized the rising islands and ceased to exist in the subsiding ones. Island hopping by a plant adapted to long-distance dispersal is biologically feasible and has been suggested for the Hawaiian genus Hillebrandia in the plant family Begoniaceae  and the fern genus Diellia .
Correlates of lobelioid diversification
In order to identify possible drivers of biotic diversification, it is essential to infer when and where lineages diversified. Ages and ancestral areas reconstructed here may, therefore, provide a suitable starting point for further examination of particular clades in the Lobelioideae.
One appealing example of how these reconstructions can be useful concerns the SCBL clade (Siphocampylus, Centropogon, Burmeistera and Lysipomia - clade N12 in Figure 3). This group of mainly nanophanerophytes and chamaephytes comprises over 580 species, which is about half of all Lobelioideae species [21, 23, 29–31]. The SCBL clade is entirely confined to the Neotropical region, and is particularly rich in species in the Northern Andes. The uplift of the Northern Andes began at ~31 Ma and intensified in the last ~20 Ma, with several discrete phases of uplift which especially affected the Eastern Cordilleras . The temporal match between the Andean uplift and the ages inferred here is striking: onset of diversification was inferred at ~27 Ma (PL) or ~19 Ma (BEAST) (Table 1) and the MRCA of this clade was unambiguously inferred as the Neotropics (Figure 9). Although a correlation in time and space does not necessarily imply causation, there is an emerging consensus that the Andean uplift has played a major role in Neotropical diversification [32–35]. Further studies are clearly needed in order to disentangle the relative role of competing hypotheses of drivers of diversification in this clade, be they mainly abiotic (for example, geotectonic events, climatic fluctuations) or biotic (for example, pollinator interactions, such as those demonstrated for certain species of Centropogon and Burmeistera [36–40]).
The temporal framework for the Lobelioideae provided here may also contribute to the discussion on the diversification of particular Lobelioideae groups, for which a fine-scale biogeographic analysis has been performed but for which no absolute divergence times have been estimated. One of these is the giant lobelias of eastern Africa, which have been thoroughly studied by Knox and co-workers [7, 13, 18, 41]. Crown ages for the eastern African clade, comprising Lobelia columnaris, L. gibberoa, and L. gregoriana, estimated here (in PL: mean 22.3 Ma; 13.7 - 30.4 Ma, 95% CI; in BEAST: median 11.2 Ma; 3.09 - 20.8 Ma, 95% HPD) clearly predate the age of most of the mountains in the Eastern and Western African Rifts . These findings influence the interpretation of the tempo and possible processes underlying diversification in the region (see  for a detailed account on competing hypotheses and their predictions).
Let us recapitulate the thought-provoking hypothesis that the giant lobelias of Africa and the Hawaiian Islands converged into the giant life form from herbaceous ancestors in the same way as the giant senecios and the silversword alliance have (pages 84-85 in ; Figure 1). The distant relationship between Dendrosenecio and the silversword alliance, and their separate diversification, has been generally accepted for more than a century. This is reflected in their taxonomic placement in different tribes based on conspicuous floral dissimilarities and was recently corroborated by molecular phylogenetic analyses [42–45]. In contrast, the results presented here show a very different evolutionary history for the giant form in the subfamily Lobelioideae. The giant lobelioids of Africa and the Hawaiian Islands, together with similarly large species from eastern Brazil, French Polynesia and southeast Asia, are all derived from a single ancestor that was woody and probably African.
Detailed comparative studies are needed in order to evaluate the morphological similarities and differences in habits among and within clades N4, N10 and N12 (Figure 3). A different way of coding for the morphology and habits of species in those clades might affect the reconstruction of the ancestral states for the Lobelioideae. Eventually, more phylogenetically informative sequence data and denser taxon sampling may help elucidate the evolution of this exceedingly diverse plant group.
Choice of taxa and molecular markers
So far, some 2600 sequences classified as belonging to family Campanulaceae have been deposited in GenBank (as of August 2009). In order to assess which of these were correctly placed in subfamily Lobelioideae, I first downloaded all Campanulaceae accessions for the most widely used sequence regions: rbcL, trnL-F and ITS. Two criteria were then used for choosing the ingroup taxa in order to infer relationships within subfamily Lobelioideae: (i) that the species could be confidently assigned to the Lobelioideae based on the results from the large-scale analyses for the Campanulaceae; and (ii) that the trnL-F region was sequenced for that species, since this region has been shown to contain considerably more phylogenetically informative characters than rbcL (41.7% as compared to 25.9%, respectively ). Based on these criteria, a total of 108 species (including seven outgroup species) were then selected for inferring the phylogeny of the Lobelioideae. This represents an increase of 33 species compared to the high-level phylogeny of the Lobelioideae recently presented by Antonelli . In that analysis, sequences for several taxa exhibiting a large/giant habit had not yet been made available on GenBank. These include representatives from the genera Apetahia and Sclerotheca (from French Polynesia), Lobelia boninensis (from the Bonin Islands), Lobelia sect. Galeatella (from the Hawaiian Islands) and Lobelia nicotianifolia (from southeast Asia), as well as many other herbaceous species.
In an attempt to increase the phylogenetic resolution among representatives of the giant lobelioids, sequences of the ITS (ITS 1 - 5.8S - ITS 2) of nuclear ribosomal DNA were generated de novo for 10 giant species plus Lobelia dortmanna. The region was amplified and sequenced with the primers ITS1, ITS2, ITS3, ITS4 , and ITS10 , following the amplification and sequencing techniques described by Antonelli . Since ITS is one of the fastest evolving molecular markers available today [47–49], its use was expected to add resolution for identifying the relationships between the main groups of giant lobelioids. Table 3 lists all species used in the analyses of subfamily Lobelioideae and indicates the species sequenced for this study.
Alignment and phylogenetic estimation
Sequences were aligned using the L-INS-I algorithm implemented in the software MAFFT v. 6 . Aligned matrices were inspected manually and all unreliable sequences excluded, iteratively, until all sequences were deemed homologous. The final aligned matrices for Campanulaceae comprised: for trnL-F, 452 sequences and 1541 characters; for ITS, 445 sequences and 1197 characters; and for rbcL, 438 sequences and 1401 characters.
Maximum Likelihood trees were inferred in the software Garli 0.960  by performing 10 independent runs with the default settings. The analyses were performed in the CIPRES cluster at the San Diego Supercomputer Center http://www.phylo.org/portal2 and the results visualized using the Interactive Tree of Life facility http://itol.embl.de. MrModeltest [52, 53] was used to select the best-fitting evolutionary model for each sequence region. The Partition Homogeneity Test  was applied to test for conflicting phylogenetic signal among these regions by performing a heuristic search with 5000 replicates, 100 random addition sequences, TBR branch swapping and saving up to 50 trees per replicate in PAUP* version 4.0b10 . Since the test did not approach significance (P = 0.962), the three markers were combined for the subsequent analyses. A bootstrap analysis was run in PAUP under the maximum parsimony criterion, by sampling 10000 replicates, with 100 random addition sequences and saving one tree per replicate. A Bayesian phylogenetic analysis was then conducted with the software MrBayes v. 3.1 , performing two parallel runs of 20 million generations each, using four chains, sampling every 1000 generations and saving branch lengths. The performance of the analyses was evaluated using the software Tracer v.1.4.1  and node frequency statistics calculated on 20000 post burn-in trees.
Several methods are currently available for estimating divergence times in a phylogeny [58, 59]. In order to compare the results from two widely used methods, divergence times were estimated using the Bayesian relaxed clock and PL approaches [60, 61]. Both methods have the advantage of enabling direct calibration on one or more nodes of a phylogeny. A major difference between them, however, is that PL assumes that rates are auto-correlated (inheritable), whereas in Bayesian relaxed clock dating each branch is allowed to evolve at its own rate.
PL estimations were done using the software r8s [61, 62]. An automated cross-validation algorithm was run to identify the optimal smoothing value for the final analysis, with log10 increments of 0.1 and using the Truncated Newton method implemented in r8s. For estimating the effects of phylogenetic uncertainty on node age estimations, 1000 trees from the stationary sample of Bayesian trees were independently dated and their statistics (mean and 95% confidence interval values) computed for each node of the Bayesian consensus tree using the software TreeAnnotator .
For the Bayesian dating analysis, five runs of 10 million generations each were performed in the software BEAST v.1.4.8  at the Computational Biology Service Unit hosted by Cornell University, USA http://cbsuapps.tc.cornell.edu. The analysis assumed a pure birth (Yule) process, since this tree prior is most suitable for inferring relationships between individuals from different species . The performance of the analysis (convergence of the independent runs, effective sample sizes) was evaluated using Tracer v.1.4.1 . Mean and 95% HPD intervals of ages were then calculated from 40,000 post burn-in trees using the software TreeAnnotator v.1.4.8 , and visualized using FigTree v1.2.2 .
Three simultaneous calibration points were applied in both analyses. (i) Since the fossil record of Campanulaceae is exceptionally scarce , the root of the Asterales tree (see ) was calibrated at 93 Ma, as estimated under a major study including 83 asterid families and based on six plastid markers and six fossil calibrations . Reliability of this age was recently corroborated by a study that employed the same set of fossils, but a different set of sequence regions and taxa (mean 94 Ma; ). In the PL analysis, 93 Ma was set as a fixed calibration, whereas in the BEAST analyses it constituted a normally distributed prior with mean 93 Ma and standard deviation of 1.0. (ii) An unchallenged fossil Campanula from the Miocene of Poland  was used to impose a minimal age constraint of 5.33 Ma for subfamily Campanuloideae. (iii) A 29.8 Ma maximum age was imposed for the radiation of the endemic Hawaiian taxa. This corresponds to the age of the oldest island of the Hawaiian Ridge (Kure), after which a continuous chain of islands has been available for propagules of the highland Hawaiian biota. In the PL analysis, this provided minimal ages for the five strongly supported (Bpp = 1.00) crown groups of Hawaiian species identified in the MrBayes analysis (see Figure 4, Results). In order to be able to use the latter as a calibration point in the BEAST analyses, it first was necessary to constrain all Hawaiian species to be monophyletic (failure to do so caused the program to crash, apparently due to a failure in generating starting trees). Since the monophyly of the Hawaiian lobelioids was demonstrably strongly supported [Bpp = 1.00, 8, 10], this a priori constraint was not expected to significantly influence the results.
Biogeographic and character evolution analyses
One problematic aspect in referring to 'giant lobelioids' is how we should characterize these plants morphologically. Giant lobelias have been described as:
'pachycaul plants; often long-lived and pliestesial; of large stature and bulk; with dense apical rosettes of typically sessile leaves, which may close every night and are retained even after they die for insulation and fire protection; with nightly closing, large racemose inflorescences; with a large apical bud of developing leaves to protect young tissue; which secrete ice-nucleating polysaccharide fluids that prevent frost damage; and a hollow pith for internal water storage that narrows basipetally'. (see  and references therein [23, 69, 70]).
These features correctly describe several of the eastern African species but they do not apply to all of them, let alone to other large lobelioids elsewhere (for example, the Andean species Siphocampylus giganteus, Figure 11). To overcome this problem, life form was coded here according to the recent checklist of the Campanulaceae by Lammers  (Table 2). Besides being a comprehensive and up-to-date account on life form and distribution for all species, this work provides a consistent framework for evaluating evolutionary traits in the family. The current distribution of all species was compiled from the same source , but simplified into eight larger operational units: Africa, French Polynesia, the Hawaiian Islands, temperate North America, the Neotropical region including temperate South America, Oceania, southeast Asia and temperate Eurasia. For both characters, the outgroup taxon (Helianthus) was coded as having a different state to any of the ingroup taxa, to prevent biasing the reconstruction of basal nodes on the tree (since Helianthus is used here to root the Campanulaceae tree but it is certainly not the sister species to the Campanulaceae). Life form and distribution for all species analysed are listed in Table 3.
The most widely used methods of ancestral area reconstruction that are able to take into account phylogenetic uncertainty and multistate characters are Fitch Parsimony (FP), implemented in the software Mesquite v. 2.7  and Bayes-DIVA [72, 73]. These two methods are based on very different biogeographic evolutionary models. FP constrains ancestors to be monomorphic (that is, restricted to single areas) and models changes in distribution from ancestor to descendant as a change in character state, equivalent to dispersal between single areas. It thus implements a dispersalist explanation. In contrast, DIVA allows widespread distributions at ancestral nodes. Although the maximum number of areas can be constrained in DIVA, single-area ancestors are not allowed and widespread distributions are always divided at speciation events by vicariance. Thus DIVA favours a more vicariant explanation . Based on the molecular dating results obtained (see Results and Discussion), and taking into consideration that the Hawaiian archipelago has never been connected to any land mass , it is clear that long-distance dispersals have played a crucial role in shaping the distribution of the Lobelioideae. This indicated that Fitch optimization was a more suitable method for inferring ancestral ranges of lobelioid nodes than DIVA. Mesquite v. 2.7  was therefore used for reconstructing both ancestral areas and ancestral life forms. To take into account phylogenetic uncertainty, reconstructions were performed on 1000 trees from the stationary Bayesian tree sample, using the Maximum Parsimony criterion and counting all trees with uniquely best states. The results were summarized by computing the relative frequencies of ancestral area reconstructions for each node of the Bayesian consensus cladogram of the Lobelioideae.
Tests of phylogenetic conservatism
In order to test whether life form and geographic distribution were phylogenetically conservative, I first calculated the number of parsimony steps necessary to explain the occurrence of each character on a particular tree, repeating this for 1000 trees randomly chosen from the post burn-in Bayesian sample (the same sample used for the dating and character state reconstruction analyses). I then generated 1000 trees with the same topology as the Bayesian 50% majority-rule consensus but with character states randomly shuffled among the tips of each tree, while keeping their relative frequency constant. The values obtained after reconstructing the characters on the simulated data set were used to compute 99% credibility intervals, to which the observed numbers of parsimony steps for each respective character could be compared. The analyses were performed in Mesquite v. 2.7 .
Bayesian evolutionary analysis by sampling trees
Bayesian posterior probability
Highest Posterior Density
internal transcribed spacer
Mega-annum (million years)
most recent common ancestor
Siphocampylus, Centropogon, Burmeistera, Lysipomia.
Linnaeus C: Systema Naturae. 1758, Laurentii Salvii, Holmiae, 1: 824-10
Givnish TJ, Sytsma KJ: Molecular Evolution and Adaptive Radiation. 1997, Cambridge: Cambridge University Press
Futuyma DJ, Edwards SV, John R: Evolution. 2005, Massachusetts: Sinauer Associates
Hedberg O: Features of afroalpine plant ecology. 1964, Stockholm: Almqvist & Wiksell
Fries RE, Fries TCE: Die riesen-lobelien Africas. Svensk Botanisk Tidskrift. 1922, 16: 383-416.
Mabberley DJ: The giant lobelias: pachycauly, biogeography, ornithophily and continental drift. New Phytologist. 1975, 74: 365-374. 10.1111/j.1469-8137.1975.tb02623.x.
Knox EB, Palmer JD: Chloroplast DNA evidence on the origin and radiation of the giant lobelias in eastern Africa. Systematic Botany. 1998, 23: 109-149. 10.2307/2419583.
Givnish TJ: Adaptive radiation, dispersal, and diversification of the Hawaiian lobeliads. The biology of biodiversity. 1999, Tokyo: Springer, 67-90.
Givnish TJ, Sytsma KJ, Smith JF, Hahn WJ: Thorn-like prickles and heterophylly in Cyanea: Adaptations to extinct avian browsers on Hawaii?. Proceedings of the National Academy of Sciences of the United States of America. 1994, 91: 2810-2814. 10.1073/pnas.91.7.2810.
Givnish TJ, Millam KC, Mast AR, Paterson TB, Theim TJ, Hipp AL, Henss JM, Smith JF, Wood KR, Sytsma KJ: Origin, adaptive radiation and diversification of the Hawaiian lobeliads (Asterales: Campanulaceae). Proceedings of the Royal Society B: Biological Sciences. 2009, 276: 407-416. 10.1098/rspb.2008.1204.
Lammers TG: Chromosome numbers of Campanulaceae. III. Review and integration of data for subfamily Lobelioideae. Amer J Bot. 1993, 80: 660-675. 10.2307/2445437.
Lammers TG, Hensold N: Chromosome numbers of Campanulaceae. II. The Lobelia tupa complex of Chile. Amer J Bot. 1992, 79: 585-588. 10.2307/2444871.
Knox EB: The conservation status of the giant senecios and giant lobelias in eastern Africa. Opera Botanica. 1993, 121: 195-216.
Embuscado ME, BeMiller JN, Knox EB: A survey and partial characterization of ice-nucleating fluids secreted by giant-rosette (Lobelia and Dendrosenecio) plants of the mountains of eastern Africa. Carbohydrate Polymers. 1996, 31: 1-9. 10.1016/S0144-8617(96)00120-8.
Givnish TJ, Montgomery RA, Goldstein G: Adaptive radiation of photosynthetic physiology in the Hawaiian lobeliads: Light regimes, static light responses, and whole-plant compensation points. American Journal of Botany. 2004, 91: 228-246. 10.3732/ajb.91.2.228.
Montgomery RA, Givnish TJ: Adaptive radiation of photosynthetic physiology in the Hawaiian lobeliads: dynamic photosynthetic responses. Oecologia. 2008, 155: 455-467. 10.1007/s00442-007-0936-3.
Montgomery RA, Goldstein G, Givnish TJ: Photoprotection of PSII in Hawaiian lobeliads from diverse light environments. Functional Plant Biology. 2008, 35: 595-605. 10.1071/FP08031.
Knox EB, Downie SR, Palmer JD: Chloroplast genome rearrangements and the evolution of giant lobelias from herbaceous ancestors. Molecular Biology and Evolution. 1993, 10: 414-430.
Givnish TJ, Knox E, Patterson TB, Hapeman JR, Palmer JD, Sytsma KJ: The Hawaiian lobelioids are monophyletic and underwent a rapid initial radiation roughly 15 million years ago. American Journal of Botany. 1996, 83: 159-
Givnish TJ, Sytsma KJ, Smith JF, Hahn WJ: Molecular evolution, adaptive radiation, and geographic speciation in Cyanea (Campanulaceae, Lobelioideae). Hawaiian Biogeography: Evolution on a Hot Spot Archipelago. 1995, 288-337.
Antonelli A: Higher level phylogeny and evolutionary trends in Campanulaceae subfam. Lobelioideae: Molecular signal overshadows morphology. Molecular Phylogenetics and Evolution. 2008, 46: 1-18. 10.1016/j.ympev.2007.06.015.
Drummond AJ, Rambaut A: BEAST: Bayesian evolutionary analysis by sampling trees. Bmc Evolutionary Biology. 2007, 7: 8-10.1186/1471-2148-7-214.
Lammers TG: World Checklist and Bibliography of Campanulaceae. 2007, Kew: Royal Botanical Gardens
Scotese CR: Atlas of Earth History. Paleogeography. 2001, Arlington:PALEOMAP Project, 1:
Koren I, Kaufman YJ, Washington R, Todd MC, Rudich Y, Martins JV, Rosenfeld D: The Bodélé depression: A single spot in the Sahara that provides most of the mineral dust to the Amazon forest. Environmental Research Letters. 2006, 1: 10.1088/1748-9326/1/1/014005.
Clague DA: The growth and subsidence of the Hawaiian-Emperor volcanic chain. The origin and evolution of Pacific Island Biotas, New Guinea to Eastern Polynesia: pattern and processes. 1996
Clement WL, Tebbitt MC, Forrest LL, Blair JE, Brouillet L, Eriksson T, Swensen SM: Phylogenetic position and biogeography of Hillebrandia sandwicensis (Begoniaceae): a rare Hawaiian relict. American Journal of Botany. 2004, 91: 905-917. 10.3732/ajb.91.6.905.
Schneider H, Ranker TA, Russell SJ, Cranfill R, Geiger JMO, Aguraiuja R, Wood KR, Grundmann M, Kloberdanz K, Vogel JC: Origin of the endemic fern genus Diellia coincides with the renewal of Hawaiian terrestrial life in the Miocene. Proc R Soc B-Biol Sci. 2005, 272: 455-460. 10.1098/rspb.2004.2965.
Lammers TG: Review of the neotropical endemics Burmeistera, Centropogon, and Siphocampylus (Campanulaceae: Lobelioideae), with description of 18 new species and a new section. Brittonia. 1998, 50: 233-262. 10.2307/2807855.
Lammers TG: Seventeen new species of Lobelioideae (Campanulaceae) from South America. Novon. 2002, 12: 206-233. 10.2307/3392958.
Knox EB, Muasya AM, Muchhala N: The predominantly South American clade of Lobeliaceae. Systematic Botany. 2008, 33: 462-468. 10.1600/036364408784571590.
Wesselingh FP, Hoorn C, Kroonenberg SB, Antonelli A, Lundberg JG, Vonhof HB, Hooghiemstra H: On the origin of Amazonian landscapes and biodiversity: a synthesys. Amazonia, Landscape and Species Evolution. Edited by: Hoorn, C, Wesselingh, FP. 2010, Oxford: Blackwell, 1
Antonelli A: Spatiotemporal Evolution of Neotropical Organisms: New Insights into an Old Riddle. Doctoral Thesis. 2008, University of Gothenburg, Department of Plant and Environmental Sciences, Göteborg, Sweden
Antonelli A, Nylander JAA, Persson C, Sanmartín I: Tracing the impact of the Andean uplift on Neotropical plant evolution. Proceedings of the National Academy of Sciences. 2009, 106: 9749-9754. 10.1073/pnas.0811421106.
Antonelli A, Quijada-Mascareñas A, Crawford AJ, Bates JM, Velazco PM, Wüster W: Molecular studies and phylogeography of Amazonian tetrapods and their relation to geological and climatic models. Amazonia, Landscape and Species Evolution. Edited by: Hoorn C, Wesselingh F. 2010, Oxford: Blackwell, 386-404. 1
Muchhala N: Exploring the boundary between pollination syndromes: Bats and hummingbirds as pollinators of Burmeistera cyclostigmata and B. tenuiflora (Campanulaceae). Oecologia. 2003, 134: 373-380.
Muchhala N: The pollination biology of Burmeistera (Campanulaceae): Specialization and syndromes. American Journal of Botany. 2006, 93: 1081-1089. 10.3732/ajb.93.8.1081.
Muchhala N: Adaptive trade-off in floral morphology mediates specialization for flowers pollinated by bats and hummingbirds. American Naturalist. 2007, 169: 494-504. 10.1086/512047.
Muchhala N: Nectar bat stows huge tongue in its rib cage. Nature. 2006, 444: 701-702. 10.1038/444701a.
Muchhala N, Potts MD: Character displacement among bat-pollinated flowers of the genus Burmeistera: analysis of mechanism, process and pattern. Proceedings of the Royal Society B: Biological Sciences. 2007, 274: 2731-2737. 10.1098/rspb.2007.0670.
Knox EB, Luke Q, Thulin M: A new giant Lobelia from the Eastern Arc Mts, Tanzania. Kew Bulletin. 2004, 59: 189-194. 10.2307/4115849.
Knox EB, Palmer JD: Chloroplast DNA variation and the recent radiation of the giant senecios (Asteraceae) on the tall mountains of eastern Africa. Proceedings of the National Academy of Sciences of the United States of America. 1995, 92: 10349-10353. 10.1073/pnas.92.22.10349.
Knox EB, Palmer JD: The origin of Dendrosenecio within the Senecioneae (Asteraceae) based on chloroplast DNA evidence. American Journal of Botany. 1995, 82: 1567-1573. 10.2307/2446185.
Panero JL, Funk VA: The value of sampling anomalous taxa in phylogenetic studies: major clades of the Asteraceae revealed. Molecular Phylogenetics and Evolution. 2008, 47: 757-782. 10.1016/j.ympev.2008.02.011.
Goertzen LR, Cannone JJ, Gutell RR, Jansen RK: ITS secondary structure derived from comparative analysis: implications for sequence alignment and phylogeny of the Asteraceae. Molecular Phylogenetics and Evolution. 2003, 29: 216-234. 10.1016/S1055-7903(03)00094-0.
White TJ, Bruns T, Lee S, Taylor J: Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. PCR Protocols: a Guide to Methods and Applications. 1990, London and New York: Academic Press, 315-322.
Hughes CE, Eastwood RJ, Bailey CD: From famine to feast? Selecting nuclear DNA sequence loci for plant species-level phylogeny reconstruction. Philosophical Transactions of the Royal Society B: Biological Sciences. 2006, 361: 211-225. 10.1098/rstb.2005.1735.
Shaw J, Lickey EB, Beck JT, Farmer SB, Liu W, Miller J, Siripun KC, Winder CT, Schilling EE, Small RL: The tortoise and the hare II: Relative utility of 21 noncoding chloroplast DNA sequences for phylogenetic analysis. American Journal of Botany. 2005, 92: 142-166. 10.3732/ajb.92.1.142.
Shaw J, Lickey EB, Schilling EE, Small RL: Comparison of whole chloroplast genome sequences to choose noncoding regions for phylogenetic studies in angiosperms: the tortoise and the hare III. American Journal of Botany. 2007, 94: 275-10.3732/ajb.94.3.275.
Katoh K, Kuma K-i, Toh H, Miyata T: MAFFT version 5: improvement in accuracy of multiple sequence alignment. Nucleic Acids Research. 2005, 33: 511-518. 10.1093/nar/gki198.
Zwickl DJ: Genetic algorithm approaches for the phylogenetic analysis of large biological sequence datasets under the maximum likelihood criterion. 2006, PhD dissertation, The University of Texas at Austin, [http://www.bio.utexas.edu/faculty/antisense/garli/Garli.html]
Nylander JAA: MrModeltest v2. Program distributed by the author. 2004, Evolutionary Biology Centre, Uppsala University
Posada D, Crandall KA: MODELTEST: Testing the model of DNA substitution. Bioinformatics. 1998, 14: 817-818. 10.1093/bioinformatics/14.9.817.
Farris JS, Kallersjo M, Kluge AG, Bult C: Testing significance of incongruence. Cladistics. 1994, 10: 315-319. 10.1111/j.1096-0031.1994.tb00181.x.
Swofford DL: PAUP*: phylogenetic analysis using parsimony (* and other methods). Version. 2002
Huelsenbeck JP, Ronquist F: MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001, 17: 754-755. 10.1093/bioinformatics/17.8.754.
Rambaut A, Drummond AJ: Tracer v14. 2007
Renner SS: Relaxed molecular clocks for dating historical plant dispersal events. Trends in Plant Science. 2005, 10: 550-558. 10.1016/j.tplants.2005.09.010.
Ricklefs RE: Estimating diversification rates from phylogenetic information. Trends in Ecology and Evolution. 2007, 22: 601-610. 10.1016/j.tree.2007.06.013.
Drummond AJ, Ho SYW, Phillips MJ, Rambaut A: Relaxed phylogenetics and dating with confidence. PLoS Biol. 2006, 4: e88-10.1371/journal.pbio.0040088.
Sanderson MJ: Estimating Absolute Rates of Molecular Evolution and Divergence Times: A Penalized Likelihood Approach. Molecular Biology and Evolution. 2002, 19: 101-109.
Sanderson MJ: r8s: Inferring absolute rates of molecular evolution and divergence times in the absence of a molecular clock. Bioinformatics. 2003, 19: 301-302. 10.1093/bioinformatics/19.2.301.
Rambaut A: FigTree v1.1. 2008, [http://tree.bio.ed.ac.uk/software/figtree]
Lammers TG: Campanulaceae. The families and genera of vascular plants. Edited by: Kubitzki K. 2007, Berlin: Springer-Verlag, 8: 26-56. full_text.
Winkworth RC, Lundberg J, Donoghue MJ: Toward a resolution of Campanulid phylogeny, with special reference to the placement of Dipsacales. Taxon. 2008, 57: 53-65.
Bremer K, Friis EM, Bremer B: Molecular phylogenetic dating of asterid flowering plants shows Early Cretaceous diversification. Systematic Biology. 2004, 53: 496-505. 10.1080/10635150490445913.
Janssens SB, Knox EB, Huysmans S, Smets EF, Merckx V: Rapid radiation of Impatiens (Balsaminaceae) during Pliocene and Pleistocene: Result of a global climate change. Molecular Phylogenetics and Evolution. 2009, 52: 806-824. 10.1016/j.ympev.2009.04.013.
Lancucka-Srodoniowa M: Macroscopic plant remains from the freswater miocene of the Nowy Sacz Basin (West Carpathians, Poland). Acta Palaeobotanica. 1979, XX: 3-117.
Lammers TG: Revision of Lobelia sect. Galeatella (Campanulaceae: Lobelioideae). Journal of the Botanical Research Institute of Texas. 2007, 1: 789-810.
Lammers TG: Revision of the endemic Hawaiian genus Trematolobelia (Campanulaceae: Lobelioideae). Brittonia. 2009, 61: 126-143.
Maddison WP, Maddison DR: Mesquite: a modular system for evolutionary analysis. Version 2.01. 2007, [http://mesquiteproject.org]
Ronquist F: Dispersal-vicariance analysis: a new approach to the quantification of historical biogeography. Systematic Biology. 1997, 46: 195-203.
Nylander JAA, Olsson U, Alström P, Sanmartín I: Accounting for Phylogenetic Uncertainty in Biogeography: A Bayesian Approach to Dispersal-Vicariance Analysis of the Thrushes (Aves: Turdus). Systematic Biology. 2008, 57: 257-268. 10.1080/10635150802044003.
Sanmartin I: Event-based biogeography: Integrating patterns, processes, and time. Biogeography in a Changing World. Edited by: Ebach MC, Tangney R. 2007, Boca Raton: Crc Press-Taylor & Francis, 135-159.
I would like to thank Roger Eriksson, Claes Persson, Aelys Humphreys, Bengt Oxelman, Thomas Lammers, Eric Knox, Uno Eliasson, Derrick Zwickl and Toby Pennington for fruitful discussions and advice, Vivian Aldén for laboratory assistance, and Gerald Carr, Frederick Warshauer and Christian Puff for kindly providing the photographs in Figure 1. Thomas Givnish and three anonymous reviewers are kindly acknowledged for reviewing and providing constructive criticism on this manuscript. This research has been financially supported by grants from the Royal Swedish Academy of Sciences, the Swedish Research Council to Claes Persson, the Royal Society of Arts and Sciences in Göteborg, Kungliga och Hvitfeldtska Stiftelsen, Helge Ax:son Johnsons Stiftelse, Carl Tryggers Stiftelse and University of Gothenburg.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.