Molecular phylogenetics reveal multiple tertiary vicariance origins of the African rain forest trees
BMC Biology volume 6, Article number: 54 (2008)
Tropical rain forests are the most diverse terrestrial ecosystems on the planet. How this diversity evolved remains largely unexplained. In Africa, rain forests are situated in two geographically isolated regions: the West-Central Guineo-Congolian region and the coastal and montane regions of East Africa. These regions have strong floristic affinities with each other, suggesting a former connection via an Eocene pan-African rain forest. High levels of endemism observed in both regions have been hypothesized to be the result of either 1) a single break-up followed by a long isolation or 2) multiple fragmentation and reconnection since the Oligocene. To test these hypotheses the evolutionary history of endemic taxa within a rain forest restricted African lineage of the plant family Annonaceae was studied. Molecular phylogenies and divergence dates were estimated using a Bayesian relaxed uncorrelated molecular clock assumption accounting for both calibration and phylogenetic uncertainties.
Our results provide strong evidence that East African endemic lineages of Annonaceae have multiple origins dated to significantly different times spanning the Oligocene and Miocene epochs. Moreover, these successive origins (c. 33, 16 and 8 million years – Myr) coincide with known periods of aridification and geological activity in Africa that would have recurrently isolated the Guineo-Congolian rain forest from the East African one. All East African taxa were found to have diversified prior to Pleistocene times.
Molecular phylogenetic dating analyses of this large pan-African clade of Annonaceae unravels an interesting pattern of diversification for rain forest restricted trees co-occurring in West/Central and East African rain forests. Our results suggest that repeated reconnections between the West/Central and East African rain forest blocks allowed for biotic exchange while the break-ups induced speciation via vicariance, enhancing the levels of endemicity. These results provide an explanation for present day distribution patterns and origins of endemicity for African rain forest trees. Moreover, given the pre-Pleistocene origins of all the studied endemic East African genera and species, these results also offer important insights for setting conservation priorities in these highly diversified but threatened ecosystems.
Tropical rain forests harbor outstanding levels of biodiversity. Numerous studies have focused on how such high levels of species richness might be maintained through time [1–3]. Much less is known about the evolutionary origins that have led to the generation of these hyper-rich ecosystems [4, 5]. Understanding the evolutionary history and biogeographic processes behind the origin of such areas plays a central role in setting conservation priorities [6, 7]. In Africa, rain forests of the Guineo-Congolian floristic region  occur in a coastal band in West Africa and occupy vast areas in Central Africa (Figure 1A). In East Africa, rain forests occur in small, highly fragmented patches near the coast and in the adjacent Eastern Arc Mountains (Figure 1A), covering a small surface area (c. 7000 km2). Despite the relatively small size of the area, the East African rain forests harbor an exceptional density of endemic plant species [9, 10], which is one of the highest on the planet , and have been identified as one of 25 global biodiversity hotspots . The Guineo-Congolian and East African regions are geographically isolated by a c. 1000 km-wide North-South arid corridor creating an effective barrier to dispersal for rain forest-restricted taxa. Nevertheless, floristic comparisons indicate strong affinities between the two rain forest regions within numerous distantly related plant families [8, 9, 11, 12]. Indeed, many endemic species are representatives of widespread genera co-occurring in the Guineo-Congolian and East African rain forest regions. In addition, the endemic East African genera are thought to be closely related to Guineo-Congolian genera. These strong affinities are generally explained by the presence of a continuous rain forest across tropical Africa during the Eocene [9–16].
One hypothesis for explaining the high levels of endemicity suggests that this pan-African rain forest broke up only once at the onset of East African aridification during the Oligocene-Early Miocene (c. 33-20 Myr) [11–13, 15, 17] and that the c. 30 Myr of isolation of the East African rain forests coupled with a highly fragmented habitat led to the high levels of endemicity observed in East Africa. In contrast, an alternative hypothesis suggests that from the mid-Tertiary onwards (c. 33-2 Myr) African rain forests expanded and contracted on multiple occasions after the initial break-up [9, 12, 14, 16]. Expansion and reconnection would have allowed new lineages to spread from West-Central Africa to East Africa or vice versa. Subsequent contraction and isolation would have resulted in diversification by vicariance or allopatric speciation, thus enhancing the levels of biodiversity and endemicity.
In the context of studying the origins of East African dry-adapted endemic species in the genus Acridocarpus (Malpighiaceae), dated molecular phylogenies have proven to be a powerful tool . Here we use a similar approach on a rain forest-restricted clade with disjunct West-Central and East African distributions, which provides a way to test the above hypotheses. If a single break-up of the pan-African rain forest is responsible for the observed distribution, then all the splits between Guineo-Congolian and East African lineages should be dated to the same period (Figure 1B). In contrast, in an oscillating forest scenario we expect to find a chronological sequence of multiple vicariance events, with splits dated at significantly different time periods (Figure 1C).
We studied the evolutionary history of the largest monophyletic African clade of rain forest trees  within the diverse pan-tropical family Annonaceae. Annonaceae is one of the best examples of a tropical plant family for which strong positive correlations were found between abundance and species richness, and temperature and rainfall, respectively [see ]). This implies that Annonaceae are mainly restricted to tropical rain forests and thus provide an ideal model family to test hypotheses on the evolution of rain forests. The studied African clade is composed of three genera endemic to the East African rain forests, two endemic to the Guineo-Congolian region, and six disjunct genera with numerous species endemic to either East or West-Central African rain forests . These genera are not found in arid parts of Africa and are restricted to forests with high rainfall regimes and relatively short dry seasons. Moreover, the comparatively large fruits found within this group are not adapted to long-distance dispersal. Within this clade, monocarps are sessile in genera with apocarpous fruits  and differ from the stalked monocarps that are more common in Annonaceae. So, even though individual fruitlets are free, in the context of dispersal these genera have fruits that functionally are equivalent to syncarpous fruits, far too large for effective long-distance dispersers such as birds. Rather, they are dispersed by mammals that are restricted to rain forests, such as gorillas , and their crossing wide arid corridors is a highly improbable event. This low dispersal capacity as well as the restriction and abundance of this clade in both the Guineo-Congolian and East African rain forests [11, 20, 22] renders this clade suitable for testing the above hypotheses.
Two dated, well-resolved molecular phylogenies were generated using plastid DNA sequences, focusing specifically on two different taxonomic levels within the African clade. The first phylogeny contained representatives of all eleven African genera and was used to date the deeper splits within this clade. Second, a species-level phylogeny was estimated for two sister genera, Isolona and Monodora, each of which has endemics in both Guineo-Congolian and East African rain forests [20, 23]. Both trees were dated under a Bayesian relaxed uncorrelated molecular clock assumption.
Taxon sampling and sequence data
For this study two different chloroplast DNA sequence datasets were used. First, in order to date the origin of the East African endemic genera, a previously published DNA sequence data matrix (dataset A in Additional file 1) of the African clade was used (Couvreur et al., 2008, treeBASE SN3554). This dataset, totaling 64 taxa, included all the genera of the African clade (11 genera) and presented a thorough sampling of the genera within the long branch clade . The family Eupomatiaceae was used as the outgroup [Eupomatiaceae has been recovered as sister to Annonaceae, [25, 26]]). The dataset was composed of six plastid markers, three non-coding (trnL-trnF, trnSG, and psbA-trnH) and three coding (ndhF, rbcL, and partial matK), totaling 7945 characters.
Second, in order to date the splits within genera, a species-level phylogeny of the two sister genera Isolona (15 out of 20 species) and Monodora (13 out of 14 species) was generated (dataset B, see Additional file 1, treeBASE number SN3633). For both genera all the known East African species were included.
A modified cetyl trimethyl ammonium bromide (CTAB) protocol of Doyle and Doyle  following Bakker et al.  was used for DNA extraction. The universal primers C/D and E/F  were used to amplify and sequence the trnL intron and trnL-trnF spacer. The psbA-trnH intergenic spacer was amplified and sequenced using primers psbA and trnH (GUG) . The trnS-trnG intergenic spacer was amplified and sequenced using primers trnS (GCU) and trnG (UCC)  and the trnD-trnT marker was also amplified and sequenced using primers trnTGGU and trnDGUC . In addition, the second part of the ndhF gene was sequenced as it has been shown to be more variable than the more conserved 5' region and is thus more appropriate for species-level analyses . The Annonaceae specific primer LBC-intF  was used in combination with the usual 2110R primer  amplifying a region of c. 620 bp. PCR reactions were performed with 30–50 ng of genomic DNA, 0.4% of BSA, 0.2 μM of each primer, 0.2 mM dNTP PCR mix (Promega, Madison, WI), 3 μM MgCl2, 1× PCR buffer (Promega, Madison, WI), and 0.5 U of Taq DNApolymerase (Promega, Madison, WI) in a total volume of 50 μl. The PCR program was as follows: 35 thermal cycles at 94°C for 1 minute, 50–55°C for 50 seconds, 72°C for 50 seconds and a final extension at 72°C for 3 minutes.
For both datasets, sequences were edited using the program Staden  and aligned manually using PAUP* [version 4.10b; ]. Gaps were coded following the simple coding model of Simmons and Ochoterena . Microsatellites were excluded from the analysis, as these structures probably originate through slipped-strand mispairing  and are likely to be highly homoplasious.
Phylogenetic analysis and divergence date estimates
In order to test whether the estimated dates of the origin of the East African endemic lineages are significantly different the molecular trees were dated under a Bayesian framework using the software BEAST v. 1.4.7 [37, 38]. The method implemented in BEAST simultaneously estimates divergence times, tree topology, and rates, thereby providing a clear advantage over previous relaxed clock methods  that estimate tree topology and divergence dates separately [e.g. [40–42]]). Both datasets were partitioned into the number of markers used by directly editing of the XLM file and following Couvreur . The best performing evolutionary model for each marker was identified under two different model selection criteria, the hierarchical likelihood ratio test (hLRT) and the Akaike information criterion [AIC; ]) as implemented in MrModelTest . For both datasets two independent analyses were undertaken to check for convergence of the MCMC chains. Analyses were undertaken by sampling every 1000th generation, and were considered complete once the effective sampling size (ESS) of each parameter was above 200, as suggested in the BEAST manual. In order to assess that the MCMC chain reached stationarity we examined the lnL plots using Tracer v. 1.3 . In particular, we searched for evidence that the model likelihood and parameter estimates had reached stationarity after a burn-in period.
The testing of the two alternative hypotheses presented in Figure 1 can be viewed under a relative time scale , i.e. with no reference to any absolute timing. However, fossil information is available and was used in order to calibrate trees in an absolute time frame. The fossil taxon Archaeanthus [98 Myr, ]) is characterized by distinctive derived stipules, an elongated receptacle, and fruits with numerous well spaced follicles, all of which are interpreted as synapomorphies shared with the family Magnoliaceae [47–49]. Thus, this fossil provides a minimum age for the stem of Magnoliaceae of 98 Myr. Even though the correct placement of this fossil within the Magnoliales has been questioned in the past , recent studies [24, 48, 51] have largely relied on Archaeanthus to date this family. Fossil information within Annonaceae is limited ([e.g. [52, 53]]) and is thought [24, 48] to be unreliable with respect to their exact placement and was thus not used in this study. One other relevant fossil is the Early Cretaceous Endressinia taxon (112 Myr ). This fossil is characterized by broad staminodes with lateral glands, which implies a possible sister relationship to Eupomatiaceae and would thus provide an older minimum age of 112 Myr for the split between Eupomatiaceae and Annonaceae than with Archaeanthus. However, the method used to assign Endressinia to Eupomatiaceae (see ) has been criticized and is thought to be unreliable . Thus, we deemed Archaeanthus the most reliable fossil calibration point and this was used to assign a minimum age to the Magnoliaceae stem node.
Using the age of Archaeanthus as well as a wider taxon sampling within Annonaceae (80 genera out of c. 110, 205 OTUs) and rbcL and trnLF sequences, Pirie  provided an age estimate of 91 Myr (± 1.5 Myr) for the stem of Annonaceae (split between Eupomatiaceae and Annonaceae) when using the Penalized Likelihood method of Sanderson . This estimated date is in accordance with previously published date estimates for the origin of Annonaceae [24, 60]. Thus, 91 Myr was used as a calibration point for dataset A. Secondary calibration (calibrating a node with a date provided by a previous analysis) is a commonly used alternative given the absence of a method for direct calibration (e.g. fossil or geological [61–63]). However, it has been shown that, unless particular care is taken, secondary calibration can generate internal inconsistencies leading to unreliable dates . BEAST accommodates for calibration uncertainty by applying a prior probability distribution on the age, e.g. a prior distribution defined in terms of its mean and standard deviation . A wide variety of prior probability distributions are available . For this study a normal probability distribution was used as it is thought to better reflect uncertainty related to secondary calibration points . For dataset A, a mean of 91 and a standard deviation of 1.5 were specified. This effectively encloses a range of possible ages from c. 89 to 93 Myr. In dataset B, a normal distribution was also used for the age of the stem node of Isolona and Monodora, with the mean taken from the analysis of dataset A. The standard deviation was set to contain the lower and higher boundaries of the 95% highest posterior density values effectively accommodating for age uncertainty.
Both sequence datasets (A and B) deviated from a strict molecular clock and rates between adjacent branches were uncorrelated as shown by the values of the parameters 'coefficient of variation' and 'covariance', respectively  (Table 1). Thus divergence times were estimated under a lognormal non-correlated relaxed clock method and using the Yule model of speciation as implemented in BEAST 1.4.7. In all cases user-specified chronogram trees were used as starting trees and were obtained using the r8s program . Finally, taxon subsets were specified for each clade of interest allowing recording of the mean time of the most recent common ancestor (t MRCA), the 95% highest posterior density intervals (HPD), and the effective sampling size (ESS). Dates were considered significantly different only when the 95% HPDs were not overlapping.
Our results provide strong support for the multiple connection-reconnection hypothesis with relative ages of splits occurring at three significantly different periods in time (Figure 2). The estimated absolute ages for all of the major clades recognized within Annonaceae are in perfect agreement with recent studies [24, 50]. Within the African clade five different origins of endemic East African lineages were identified (Figure 2A and 2B, nodes a-e). Node 'a' represents the origin of the Sanrafaelia/Ophrypetalum clade with a mean age of 32.9 Myr (95% HPD: 42.9–23.6 Myr) and node 'b' that of Asteranthe with a mean age of 16.8 (HPD: 23.4–10.5). Node 'c' indicates the origin of the East African endemic species Uvariodendron kirkii, with a mean of 8.4 Myr (HPD: 13.2–3.7). Two further nodes are indicated in Figure 2B, with the first one corresponding to the origin of the East African species of Monodora (node 'd'), with a mean of 8.4 Myr (HPD: 12.2–4.7). The lack of phylogenetic resolution makes it impossible to distinguish the origin of the East African species of Isolona (Figure 2B). Therefore, only a minimum age for this potential split can be provided corresponding to node 'e' with a mean of 5.4 Myr (HPD: 8.4–2.5). Each age estimate (Figure 2C) for these nodes is dated to three significantly different periods of time as the 95% HPDs do not overlap. For node 'c', 'd', and 'e' the 95% HPDs do largely overlap, indicating no significant differences in the age estimates.
The chronograms and phylogenetic relationships resulting from both analyses are presented in Figure 2A and 2B. At the generic level (dataset A) the resulting phylogenetic relationships as well as corresponding branch support were identical to those found in Couvreur et al. . The species-level chronogram (dataset B) of the two sister genera Isolona and Monodora is represented in Figure 2B. Monodora has two well-supported sister clades (Bayesian posterior probability – PP = 1.00): a West African clade, and an East African clade. Relationships within the West African clade are weakly supported except for the sister position of Monodora angolensis (PP = 1.00) to the remaining West African species. In the East African clade, Monodora grandidieri is strongly supported as sister to the rest of the taxa. Moreover, Monodora globiflora with Monodora carolinae and Monodora stenopetala, as well as Monodora hastipetala with Monodora junodii form well-supported subclades (PP = 1.00). The Malagasy species of Isolona form a highly supported, nested clade. The East African taxa are also strongly supported as nested within Isolona; however, their relationship with the other clades remains unresolved.
Our results demonstrate that the splits between the Guineo-Congolian and East African Annonaceae lineages have multiple origins at significantly different times (Figure 2), which is congruent with the hypothesis presented in Figure 1C. Given that the majority of species found within this clade are restricted to rain forests and present no long-range dispersal potential, we interpret the observed pattern as a likely consequence of a series of connection-isolation events between the East African and Guineo-Congolian forests. Our results were obtained with just one plant group, but the biogeographic patterns observed in Annonaceae are common amongst many other tree families in Africa in terms of ecology (taxa restricted to rain forests) and distribution (disjunct between West/Central and East Africa) [8, 11, 12]. Therefore, we predict that the evolutionary process of divergence by vicariance on multiple occasions could be a common pattern for numerous other plant taxa in Africa, thus explaining much of the plant distribution and endemicity within African rain forest tree families. Interestingly, multiple splits between East and West-Central rain forest dwelling lineages at different times were also found to have occurred in African caecilian amphibians . Thus, it would appear that this process could also be the main factor behind present day distribution patterns for a wide range of African rain forest restricted organisms, including animals.
Our results also provide important insights into the absolute time scale of these events, especially for the origin of the East African endemic lineages. The dates inferred from our study were obtained based on time calibration using the fossil taxon Archaeanthus (see Methods) as it is generally accepted to represent a minimum age constraint for the stem of Magnoliaceae [24, 49, 50, 53, 58, 59]. Even acknowledging the complex history of vegetation change in East Africa  and the uncertainty in age estimates, the estimated origins of these splits using the Archaeanthus fossil calibration coincide well with known periods of renewed aridity and/or continental uplift in East Africa. These phenomena are generally thought to have played a role in the break-up of the pan-African rain forest, thus providing an explanation for the observed splits, which can be looked upon as further support of their correctness. The first split is dated to the Oligocene (mean: 32.9 Myr; Figure 2A – a), which corresponds to a period of drastic global cooling resulting, in part, from the development of permanent continental ice-sheets in Antarctica . This cooling induced a drier climate at equatorial levels in Africa, fragmenting the pan-African Eocene rain forest . The second split, dated to the Late-Early Miocene (mean: 16.8, Figure 2A – b), post-dates an Early Miocene warmer and moister climate during which rain forest is thought to have extended again from coast to coast [13, 15], and coincides with the Miocene climatic optimum (17-15 Myr) characterized by high global temperatures . This period also witnessed the closure of the Tethys sea when the African plate collided with the Eurasian one (c. 18-17 Myr), bringing an end to the moist influence of the latitudinal oceanic circulation system [15, 69]. These drier conditions, as well as higher global temperatures, induced a new period of marked aridity, allowing the spread of savannas at the expense of rain forests [15, 69]. This is visible at the paleobotanical level by an increase in grass abundance across Africa by 16 Myr . The third significant split (mean: 8.4–5.4 Myr) occurred in three lineages, Uvariodendron (Figure 2A – c), Monodora (Figure 2B – d), and Isolona (Figure 2B – e), and took place after the initiation of geological activity in the western East African Rift System that uplifted the central Tanganyikan plateau (c. 10 Myr) [12, 70]. Such uplifting has been shown to have played a significant role in the aridification of eastern Africa . Moreover, this estimated date also coincides with a period of significant extension of savannas in East Africa [16, 69], as judged by an increase in biomass of plants using C4 photosynthetic pathways (8-6 Myr) .
Furthermore, even with our fossil-based minimum age estimates most speciation events in East African Isolona and Monodora (Figure 2B) have a pre-Pleistocene origin indicating a lack of diversification in the East African rain forests during the climatic fluctuations of the Pleistocene (e.g. ). For East Africa, this pattern was also documented in clawed frogs where all extant East African taxa were found to have originated before the Pleistocene . Paleoecological records from the Eastern Arcs in Tanzania spanning the last 38,000 years provided evidence of ecosystem stability during the last glacial maximum period . This lack of diversification during the Pleistocene is generally thought to be due to the stabilizing moist influence of the Indian Ocean on tropical rain forests in parts of coastal and montane East Africa, allowing them to persist throughout glacial periods without any important fluctuations in size [12, 75]. Therefore, our results provide strong evidence that a larger Oligocene-Miocene continental-scale fragmentation between, rather than within, East and West-Central rain forest regions has played a major role in generating the patterns in plant distribution and endemicity that are currently observed across Africa. However, it may be that other, possibly herbaceous, groups diversified in East Africa as a result of Pleistocene climatic change but this is certainly not the case for Annonaceae and is likely also not true for other woody representatives of the East African rain forest flora.
Finally, our results indicate that the highly fragmented and threatened rain forests of Eastern Africa contain large amounts of pre-Pleistocene derived endemic lineages having evolutionary histories that date back to the Early Oligocene. A large majority of these lineages are species-poor , making them vulnerable to extinction. Additionally, most of the endemic species or genera are highly vulnerable to extinction because of small distributional ranges . Subsequently, they have been red-listed with some level of threat to their survival (an estimated 22–25% of the total number of endemic species [Roy Gereau, personal communication]). Thus, because of their ancient origins and species-poor composition, the extinction of just a handful of species could result in the loss of significant phylogenetic diversity . This information should be carefully and urgently taken into consideration in future conservation planning within the East African region as knowledge of the presence of such ancient lineages allows us to prioritize areas with high phylogenetic diversity [6, 76, 77].
Our molecular phylogenetic study on the evolutionary history of this pan-African Annonaceae rain forest clade provided numerous important insights into the evolutionary origins of African rain forest trees. We show that the East African lineages within this large clade have significantly different temporal origins. This would suggest that present day distribution and endemicity in African rain forest is a likely consequence of a series of connection-isolation events between the East African and Guineo-Congolian forests. Moreover, all the East African taxa appear to have diversified before the Pleistocene climatic oscillations. This provides evidence for the increasingly probable hypothesis that East African rain forests have been ecologically stable during those times. Finally, the old pre-Pleistocene origins of the endemic East African taxa suggest that these biodiversity hotspot regions contain high levels of phylogenetic diversity, an important factor to take into account for conservation priorities.
- PP :
Bayesian posterior probability
highest posterior distributions.
Hubbell SP: A unified theory of biogeography and relative species abundance and its application to tropical rain forests and coral reefs. Coral Reefs. 1997, 16: S9-S21. 10.1007/s003380050237.
Molino JF, Sabatier D: Tree diversity in tropical rain forests: A validation of the intermediate disturbance hypothesis. Science. 2001, 294: 1702-1704. 10.1126/science.1060284.
Wills C, Harms KE, Condit R, King D, Thompson J, He F, Muller-Landau HC, Ashton P, Losos E, Comita L, Hubbell S, LaFrankie J, Bunyavejchewin S, Dattaraja HS, Davies S, Esufali S, Foster R, Gunatilleke N, Gunatilleke S, Hall P, Itoh A, John R, Kiratiprayoon S, Loo de Lao S, Massa M, Nath C, Supardi Noor MN, Rahman Kassim A, Sukumar R, Suresh HS, Sun I, Tan S, Yamakura T, Zimmerman J: Nonrandom processes maintain diversity in tropical forests. Science. 2006, 311: 527-531. 10.1126/science.1117715.
Moritz C, Patton JL, Schneider CJ, Smith TB: Diversification of rainforest faunas: An integrated molecular approach. Annu Rev Ecol Syst. 2000, 31: 533-563. 10.1146/annurev.ecolsys.31.1.533.
Pennington RT, Richardson JE, Lavin M: Insights into the historical construction of species-rich biomes from dated plant phylogenies, neutral ecological theory and phylogenetic community structure. New Phytol. 2006, 172: 605-616. 10.1111/j.1469-8137.2006.01902.x.
Mace GM, Gittleman JL, Purvis A: Preserving the Tree of Life. Science. 2003, 300: 1707-1709. 10.1126/science.1085510.
Myers N, Mittermeier RA, Mittermeier CG, da Fonseca GAB, Kent J: Biodiversity hotspots for conservation priorities. Nature. 2000, 403: 853-858. 10.1038/35002501.
White F: The Guineo-Congolian region and its relationships to other phytochoria. Bull Jard Bot Nat Belg. 1979, 49: 11-55. 10.2307/3667815.
Burgess ND, Clarke GP, Rodgers WA: Coastal forests of eastern Africa: status, endemism patterns and their potential causes. Biol J Linn Soc Lond. 1998, 64: 337-10.1111/j.1095-8312.1998.tb00337.x.
Burgess ND, Butynski TM, Cordeiro NJ, Doggart NH, Fjeldså J, Howell KM, Kilahama FB, Loader SP, Lovett JC, Mbilinyi B, Menegon M, Moyer DC, Nashanda E, Perkin A, Rovero F, Stanley WT, Stuart SN: The biological importance of the Eastern Arc Mountains of Tanzania and Kenya. Biol Conserv. 2007, 134: 209-231. 10.1016/j.biocon.2006.08.015.
Burgess ND, Clarke GP: Coastal Forests of Eastern Africa. 2000, Gland: IUCN
Wasser SK, Lovett JC: Biogeography and Ecology of the Rainforests of Eastern Africa. 1993, Cambridge: Cambridge University Press
Morley RJ: Origin and Evolution of Tropical Rain Forests. 2000, New York: John Wiley & Sons
Coetzee JA: African flora since the terminal Jurassic. Biological relationships between Africa and South America. Edited by: Goldblatt P. 1993, New Haven: Yale University Press, 37-61.
Axelrod DI, Raven PH: Late Cretaceous and Tertiary vegetation history of Africa. Biogeography and Ecology of Southern Africa. Edited by: Werger MJA. 1978, The Hague: W Junk bv Publishers, 77-130.
Jacobs BF, Kingston JD, Jacobs LL: The origin of grass-dominated ecosystems. Ann Missouri Bot Gard. 1999, 86: 590-643. 10.2307/2666186.
Davis CC, Bell CD, Fritsch PW, Mathews S: Phylogeny of Acridocarpus-Brachylophon (Malpighiaceae): implications for tertiary tropical floras and Afroasian biogeography. Evolution. 2002, 56: 2395-2405.
Couvreur TLP, Richardson JE, Sosef MSM, Erkens RHJ, Chatrou LW: Evolution of syncarpy and other morphological characters in African Annonaceae: a posterior mapping approach. Mol Phylogenet Evol. 2008, 47: 302-318. 10.1016/j.ympev.2008.01.018.
Punyasena SW, Eshel G, McElwain JC: The influence of climate on the spatial patterning of Neotropical plant families. J Biogeogr. 2008, 35: 117-130.
Couvreur TLP, Gereau RE, Wieringa JJ, Richardson JE: Description of four new species of Monodora and Isolona (Annonaceae) from Tanzania and an overview of Tanzanian Annonaceae diversity. Adansonia (Paris). 2006, 28: 243-266.
Rogers ME, Abernethy K, Bermejo M, Cipolletta C, Doran D, McFarland K, Nishihara T, Remis M, Tutin CEG: Western gorilla diet: A synthesis from six sites. Am J Primatol. 2004, 64: 173-192. 10.1002/ajp.20071.
Gentry A: Diversity and floristic composition of lowland tropical forest in Africa and South America. Biological Relationships between Africa and South America. Edited by: Goldblatt P. 1993, New Haven: Yale University Press, 500-547.
Couvreur TLP: Monograph of the syncarpous African genera Isolona and Monodora(Annonaceae). Syst Bot Monogr.
Richardson JE, Chatrou LW, Mols JB, Erkens RHJ, Pirie MD: Historical biogeography of two cosmopolitan families of flowering plants: Annonaceae and Rhamnaceae. Philos Trans R Soc Lond B Biol Sci. 2004, 359: 1495-1508. 10.1098/rstb.2004.1537.
Qiu YL, Lee J, Bernasconi-Quadroni F, Soltis DE, Soltis PS, Zanis M, Zimmer EA, Chen Z, Savolainen V, Chase MW: Phylogeny of basal angiosperms: Analyses of five genes from three genomes. Int J Plant Sci. 2000, 161: S3-S27. 10.1086/317584.
Sauquet H, Doyle JA, Scharaschkin T, Borsch T, Hilu KW, Chatrou LW, Le Thomas A: Phylogenetic analysis of Magnoliales and Myristicaceae based on multiple data sets: implications for character evolution. Bot J Linn Soc. 2003, 142: 125-186. 10.1046/j.1095-8339.2003.00171.x.
Doyle JJ, Doyle JL: A rapid DNA isolation procedure from small quantities of fresh leaf tissue. Phyt Bull. 1987, 19: 11-15.
Bakker FT, Hellbriigge D, Culham A, Gibby M: Phylogenetic relationships within Pelargonium sect. Peristera (Geraniaceae) inferred from nrDNA and cpDNA sequence comparisons. Plant Syst Evol. 1998, 211: 273-287. 10.1007/BF00985364.
Taberlet P, Gielly L, Pautou G, Bouvet J: Universal primers for amplification of 3 noncoding regions of chloroplast DNA. Plant Mol Biol. 1991, 17: 1105-1109. 10.1007/BF00037152.
Hamilton MB: Four primer pairs for the amplification of chloroplast intergenic regions with intraspecific variation. Mol Ecol. 1999, 8: 521-523.
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. Am J Bot. 2005, 92: 142-166. 10.3732/ajb.92.1.142.
Erkens RHJ, Chatrou LW, Maas JW, Pirie MD: Phylogenetic relationships, saturation and marker-use in the Long Branch Clade of Annonaceae. PhD thesis, From morphological nightmare to molecular conundrum. Phylogenetic, evolutionary and taxonomic studies on Guatteria (Annonaceae). Edited by: Erkens RHJ. 2007, Utrecht University, Utrecht, 25-41.
Olmstead RG, Sweere JA: Combining data in phylogenetic systematics: an empirical approach using three molecular data sets in the Solanaceae. Syst Biol. 1994, 43: 467-481. 10.2307/2413546.
Staden Package WWW site. [http://staden.sourceforge.net]
Swofford DL: PAUP* Phylogenetic Analysis Using Parsimony (* and other methods), v. 4.0 beta 10. 2002, Sunderland: Sinauer Associates
Simmons MP, Ochoterena H: Gaps as characters in sequence-based phylogenetic analyses. Syst Biol. 2000, 49: 369-381. 10.1080/10635159950173889.
Levinson G, Gutman G: Slipped-strand mispairing: a major mechanism for DNA sequence evolution. Mol Biol Evol. 1987, 4: 203-221.
Drummond AJ, Ho SYW, Phillips MJ, Rambaut A: Relaxed phylogenetics and dating with confidence. PLoS Biol. 2006, 4: 699-710. 10.1371/journal.pbio.0040088.
Drummond AJ, Rambaut A: BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007, 7: 214-10.1186/1471-2148-7-214.
Ho SYW, Phillips MJ, Drummond AJ, Cooper A: Accuracy of rate estimation using relaxed-clock models with a critical focus on the early metazoan radiation. Mol Biol Evol. 2005, 22: 1355-1363. 10.1093/molbev/msi125.
Sanderson M: A nonparametric approach to estimating divergence times in the absence of rate constancy. Mol Biol Evol. 1997, 14: 1218-1231.
Thorne JL, Kishino H, Painter IS: Estimating the rate of evolution of the rate of molecular evolution. Mol Biol Evol. 1998, 15: 1647-1657.
Sanderson MJ: Estimating absolute rates of molecular evolution and divergence times: A penalized likelihood approach. Mol Biol Evol. 2002, 19: 101-109.
BEAST partitioning. [http://tlpcouvreur.googlepages.com/beastpartitioning]
Akaike H: Information theory as an extension of the maximum likelihood principle. Second International Symposium on Information Theory. Edited by: Petrov BN, Csaki F. 1973, Budapest: Akademiai Kiado, 267-281.
Nylander JAA: MrModeltest v2. Program distributed by the author. 2004, Uppsala: Evolutionary Biology Centre, Uppsala University
Rambaut A, Drummond AJ: Tracer. Version 1.4. 2003, [http://evolve.zoo.ox.ac.uk/]
Loader SP, Pisani D, Cotton JA, Gower DJ, Day JJ, Wilkinson M: Relative time scales reveal multiple origins of parallel disjunct distributions of African caecilian amphibians. Biol Lett. 2007, 3: 505-508. 10.1098/rsbl.2007.0266.
Dilcher DL, Crane PR: Archaeanthus – an early angiosperm from the Cenomanian of the western interior of North-America. Ann Missouri Bot Gard. 1984, 71: 351-383. 10.2307/2399030.
Doyle JA, Sauquet H, Scharaschkin T, Le Thomas A: Phylogeny, molecular and fossil dating, and biogeographic history of Annonaceae and Myristicaceae (Magnoliales). Int J Plant Sci. 2004, 165: S55-S67. 10.1086/421068.
Friis EM, Pedersen KR, Crane PR: Cretaceous anglosperin flowers: Innovation and evolution in plant reproduction. Palaeogeogr Palaeoclimatol Palaeoecol. 2006, 232: 251-293. 10.1016/j.palaeo.2005.07.006.
Crepet WL, Nixon KC, Gandolfo MA: Fossil evidence and phylogeny: The age of major angiosperm clades based on mesofossil and macrofossil evidence from cretaceous deposits. Am J Bot. 2004, 91: 1666-1682. 10.3732/ajb.91.10.1666.
Pirie MD, Chatrou LW, Mols JB, Erkens RHJ, Oosterhof J: 'Andean-centred' genera in the short-branch clade of Annonaceae: testing biogeographical hypotheses using phylogeny reconstruction and molecular dating. J Biogeogr. 2006, 33: 31-46. 10.1111/j.1365-2699.2005.01388.x.
Chesters KIM: Some plant remains from the Upper Cretaceous and Tertiary of West Africa: Maastrichtian seeds of Annonaceae. Ann Mag Nat Hist Ser. 1955, 12: 498-504.
Kvacek J, Eklund H: A report on newly recovered reproductive structures from the Cenomanian of Bohemia (Central Europe). Int J Plant Sci. 2003, 164: 1021-1039. 10.1086/378824.
Mohr BAR, Bernardes-de-Oliveira MEC: Endressinia brasiliana, a magnolialean angiosperm from the Lower Cretaceous Crato Formation (Brazil). Int J Plant Sci. 2004, 165: 1121-1133. 10.1086/423879.
Mohr BAR, Eklund H: Araripia florifera, a magnoliid angiosperm from the Lower Cretaceous Crato Formation (Brazil). Rev Palaeobot Palynol. 2003, 126: 279-292. 10.1016/S0034-6667(03)00092-7.
Crepet WL: The fossil record of angiosperms: Requiem or renaissance?. Ann Missouri Bot Gard. 2008, 95: 3-33. 10.3417/2007065.
Pirie MD: The age of clades in Annonaceae: progress in molecular dating. The Natural History of the Annonaceae. Edited by: Chatrou LW, Erkens RHJ, Richardson JE. Regnum Vegetabile
Wikström N, Savolainen V, Chase MW: Evolution of the angiosperms: calibrating the family tree. Proc Biol Sci. 2001, 268 (1482): 2211-2220. 10.1098/rspb.2001.1782.
Zhou SL, Renner SS, Wen J: Molecular phylogeny and intra- and intercontinental biogeography of Calycanthaceae. Mol Phylogenet Evol. 2006, 39: 1-15. 10.1016/j.ympev.2006.01.015.
Berry PE, Hahn WJ, Sytsma KJ, Hall JC, Mast A: Phylogenetic relationships and biogeography of Fuchsia (Onagraceae) based on noncoding nuclear and chloroplast DNA data. Am J Bot. 2004, 91: 601-614. 10.3732/ajb.91.4.601.
Renner SS: Relaxed molecular clocks for dating historical plant dispersal events. Trends Plant Sci. 2005, 10: 550-558. 10.1016/j.tplants.2005.09.010.
Shaul S, Graur D: Playing chicken (Gallus gallus): methodological inconsistencies of molecular divergence date estimates due to secondary calibration points. Gene. 2002, 300: 59-61. 10.1016/S0378-1119(02)00851-X.
Ho SYW: Calibrating molecular estimates of substitution rates and divergence times in birds. J Avian Biol. 2007, 38: 409-414.
Peng ZG, Ho SYW, Zhang YG, He SP: Uplift of the Tibetan plateau: Evidence from divergence times of glyptosternoid catfishes. Mol Phylogenet Evol. 2006, 39: 568-572. 10.1016/j.ympev.2005.10.016.
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.
Zachos J, Pagani M, Sloan L, Thomas E, Billups K: Trends, rhythms, and aberrations in global climate 65 Ma to present. Science. 2001, 292: 686-693. 10.1126/science.1059412.
Jacobs BF: Palaeobotanical studies from tropical Africa: relevance to the evolution of forest, woodland and savannah biomes. Philos Trans R Soc Lond B Biol Sci. 2004, 359: 1573-1583. 10.1098/rstb.2004.1533.
Sepulchre P, Ramstein G, Fluteau F, Schuster M, Tiercelin JJ, Brunet M: Tectonic uplift and Eastern Africa aridification. Science. 2006, 313: 1419-1423. 10.1126/science.1129158.
Cerling TE, Harris JM, MacFadden BJ, Leakey MG, Quade J, Eisenmann V, Ehleringer JR: Global vegetation change through the Miocene/Pliocene boundary. Nature. 1997, 389: 153-158. 10.1038/38229.
Haffer J: Speciation in Amazonian forest birds. Science. 1969, 165: 131-137. 10.1126/science.165.3889.131.
Evans BJ, Kelley DB, Tinsley RC, Melnick DJ, Cannatella DC: A mitochondrial DNA phylogeny of African clawed frogs: phylogeography and implications for polyploid evolution. Mol Phylogenet Evol. 2004, 33: 197-213. 10.1016/j.ympev.2004.04.018.
Mumbi CT, Marchant R, Hooghiemstra H, Wooller MJ: Late Quaternary vegetation reconstruction from the Eastern Arc Mountains, Tanzania. Quaternary Research. 2008, 69: 326-341. 10.1016/j.yqres.2007.10.012.
Marchant R, Mumbi C, Behera S, Yamagata T: The Indian Ocean dipole – the unsung driver of climatic variability in East Africa. Afr J Ecol. 2007, 45: 4-16. 10.1111/j.1365-2028.2006.00707.x.
Yamagata T, Balmford A, Manning JC, Proches S, Bank van der M, Reeves G, Hedderson TAJ, Savolainen V: Preserving the evolutionary potential of floras in biodiversity hotspots. Nature. 2007, 445: 757-760. 10.1038/nature05587.
Sechrest W, Brooks TM, da Fonseca GAB, Konstant WR, Mittermeier RA, Purvis A, Rylands AB, Gittleman JL: Hotspots and the conservation of evolutionary history. Proc Natl Acad Sci USA. 2002, 99: 2067-2071. 10.1073/pnas.251680798.
We thank JJ Wieringa, BF Jacobs, JC Lovett, Q Luke, DJ Harris, E Haston, RT Pennington, and T van der Niet for useful comments and discussions. We also acknowledge A Robertson, D Johnson, and the Missouri Botanical Garden for providing plant DNA material. We are grateful to FM Mbago, RE Gereau, L Ngok Balak, Y Issembe, R Niangadouma, and M Botermans for their excellent help and assistance in the field. The governmental authorities of Gabon and Tanzania (COSTECH) as well as national park directors are thanked for granting collection permits. Funding for fieldwork came from the Netherlands Organization for Scientific Research (NWO), National Geographic Society, Alberta Mennega Stichting, Hugo de Vries Fonds, and Air France-KLM.
TLPC carried out the molecular genetic studies, the sequence alignment, phylogenetic and dating analysis, and drafted the manuscript. LWC participated in the sequence alignment, supervised the study and helped to draft the manuscript. MSM participated in the design of the study. JER conceived the study, and participated in its design and coordination, and helped to draft the manuscript. All authors critically read and improved earlier drafts and approved the final manuscript.
Electronic supplementary material
Additional file 1: Vouchers and genebank info. Taxon sampling, voucher information and GenBank accession numbers for each of the six (dataset A) and five (dataset B) chloroplast markers used in this study. UUCB: University Utrecht Botanical Garden; DRC: Democratic Republic of Congo. (XLS 45 KB)
About this article
Cite this article
Couvreur, T.L., Chatrou, L.W., Sosef, M.S. et al. Molecular phylogenetics reveal multiple tertiary vicariance origins of the African rain forest trees . BMC Biol 6, 54 (2008). https://doi.org/10.1186/1741-7007-6-54