© 2004 Society of Systematic Biologists
Southern Hemisphere Biogeography Inferred by Event-Based Models: Plant versus Animal Patterns
Edited by Cliff Cunningham: Associate Editor
Department of Systematic Zoology, Evolutionary Biology Centre Norbyvägen 18D, SE-752 36 Uppsala, Sweden; E-mail: isabel.sanmartin{at}ebc.uu.se (I.S.)
| Abstract |
|---|
|
|
|---|
The Southern Hemisphere has traditionally been considered as having a fundamentally vicariant history. The common trans-Pacific disjunctions are usually explained by the sequential breakup of the supercontinent Gondwana during the last 165 million years, causing successive division of an ancestral biota. However, recent biogeographic studies, based on molecular estimates and more accurate paleogeographic reconstructions, indicate that dispersal may have been more important than traditionally assumed. We examined the relative roles played by vicariance and dispersal in shaping Southern Hemisphere biotas by analyzing a large data set of 54 animal and 19 plant phylogenies, including marsupials, ratites, and southern beeches (1,393 terminals). Parsimony-based tree fitting in conjunction with permutation tests was used to examine to what extent Southern Hemisphere biogeographic patterns fit the breakup sequence of Gondwana and to identify concordant dispersal patterns. Consistent with other studies, the animal data are congruent with the geological sequence of Gondwana breakup: (Africa(New Zealand(southern South America, Australia))). Trans-Antarctic dispersal (Australia
southern South America) is also significantly more frequent than any other dispersal event in animals, which may be explained by the long period of geological contact between Australia and South America via Antarctica. In contrast, the dominant pattern in plants, (southern South America(Australia, New Zealand)), is better explained by dispersal, particularly the prevalence of trans-Tasman dispersal between New Zealand and Australia. Our results also confirm the hybrid origin of the South American biota: there has been surprisingly little biotic exchange between the northern tropical and the southern temperate regions of South America, especially for animals.
Keywords: Concordant dispersal; Gondwana; historical biogeography; parsimony-based tree fitting; Southern Hemisphere; vicariant history; West Wind Drift
Received December 18, 2002; Revised May 14, 2003; Accepted November 15, 2003
Dispersal and vicariance are often considered competing hypotheses in historical biogeography. Disjunct distribution patterns can be explained either by fragmentation of widespread ancestors by vicariant (isolating) events or by dispersal across a preexisting barrier. Current methods of biogeographic analysis are based on the vicariance model because nearly any distribution pattern can be explained by dispersal, making dispersal hypotheses resilient to falsification (Morrone and Crisci, 1995). Vicariance hypotheses, however, can be tested by concordance between the phylogenetic and distribution patterns of different organisms. Taxa exhibiting the same phylogenetic and distribution pattern are assumed to share a common history of geological isolation events. The resulting pattern can be represented in the form of a hierarchical area cladogram, presumably reflecting the biotic relationships among the areas analyzed (Nelson and Platnick, 1981).
The biogeographic history of the Southern Hemisphere is considered a prime example of the vicariance scenario: disjunct trans-Pacific distributions resulted from the sequential breakup of the southern supercontinent Gondwana during the last 165 million years (MY) (Fig. 1), causing vicariant division of its ancestral biota (Brundin, 1966; Raven and Axelrod, 1972). Long-distance dispersal is usually discarded as the primary causal explanation of Southern Hemisphere patterns. Most Gondwanan groups (e.g., marsupials, ratites) are presumed to be poor dispersers, unable to move across oceanic barriers. Moreover, dispersal is thought to be a random process, and hence it could not have given rise to the type of congruent or concordant patterns found in so many different groups (Croizat et al., 1974; Craw, 1982).
|
Paleogeographic reconstructions, however, indicate that the biogeographic history of the Southern Hemisphere cannot be entirely reduced to a simple sequence of vicariance events. For instance, some of the southern landmasses (e.g., New Guinea) were formed by the accretion of multiple terranes (Pigram and Davies, 1987), whereas others (India) became isolated only to join new landmasses later on. Although recent molecular estimates indicate that the Cretaceous breakup of Gondwana was an important mechanism in the ordinal diversification of birds and mammals (Cracraft, 2001; Eizirik et al., 2001), molecular studies of other groups increasingly indicate that dispersal has played a more important role in the development of the Southern Hemisphere biota than traditionally assumed. Molecular clock estimates of divergence times show that many classic Gondwanan groups combine an ancient vicariance pattern with relatively recent dispersal events (Waters et al., 2000b; Cooper et al., 2001), whereas other groups originated after continental breakup, and their distribution can be explained only by long-distance dispersal across oceanic barriers (Baum et al., 1998; Waters et al., 2000a; Buckley et al., 2002). The notion that Gondwanan groups are poor dispersers has also been questioned by their occurrence on young volcanic islands (Nothofagus; Pole, 1994) or their presence in Gondwana landmasses that were isolated at the time of cladogenesis (e.g., ratite kiwis in New Zealand; Cooper et al., 2001). The flora of New Zealand and New Caledonia has traditionally been considered a relict of the late Cretaceous Gondwanan flora, isolated since the time of vicariance. However, according to Pole (1994), paleobotanical records indicate that most of New Zealand's present flora came by long-distance dispersal from Australia during the late Tertiary after land connections with the rest of Gondwana were broken. Trans-Tasman dispersal via the West Wind Drift is often argued to explain the frequency of relationships ((Australia, New Zealand)South America) in austral plants, a pattern in conflict with the accepted sequence of Gondwana breakup (Linder and Crisp, 1995).
Despite numerous studies of single groups, a comprehensive analysis of Southern Hemisphere biogeographic patterns has not been attempted yet. The shortcomings of studies based on only a few groups and the necessity of large data sets in separating general biogeographic patterns from distributional noise are widely recognized (Crisci et al., 1991; Seberg, 1991; Linder and Crisp, 1995; Sanmartín et al., 2001). At the same time, recent methodological developments have made it possible to study complex biogeographic patterns in much more detail than previously. Thus, the time is ripe for a new look at Southern Hemisphere biogeography.
Traditional biogeographic methods are characterized by a relatively lose relation between pattern and process. This causes a number of problems, including difficulties in inferring and testing biogeographic hypotheses that include both dispersal and vicariance patterns. In contrast, recent event-based methods (Page, 1995; Ronquist, 1997, 1998) are derived from explicit (although not strict probabilistic) process models and reconstruct both ancestral distributions and biogeographic events. Each event (e.g., vicariance, duplication, dispersal, and extinction) is associated with a cost that should be inversely related to its likelihood: the more likely the event, the less costly it should be (Ronquist and Nylin, 1990; Ronquist, 1998). The most-parsimonious solution is the one minimizing the total cost of the implied events (Ronquist, 1995, 1998, 2002a); it can be searched for using modifications of standard heuristic and exhaustive parsimony search algorithms (Ronquist, 1995, 1998, 2002a).
Dispersal–vicariance analysis (DIVA; Ronquist, 1997) is an event-based biogeographic inference method that does not force vicariance events to conform to a hierarchical pattern. DIVA is particularly appropriate when inferring biogeographic patterns in complex geological settings where areas have split and fused again (Sanmartín et al., 2001), but its generality causes unnecessary loss of statistical power when there is a strong underlying hierarchical pattern of vicariance events, as one might expect in the Southern Hemisphere. In this situation, it is more appropriate to use event-based tree fitting methods, which only allow vicariance events predicted by a hierarchical area cladogram (Page, 1995; Ronquist, 1998, 2002a). Tree-fitting methods have been used extensively for studying host–parasite systems but are still relatively unexplored in the biogeographic context. These methods are good for finding significantly supported area cladograms, and if these cladograms exist tree fitting will help to reduce the noise level in the data when looking at other patterns. For instance, tree fitting allows one to study dispersal patterns while accounting for the vicariance events predicted by a geologically based area cladogram. In the absence of geologically predicted vicariance events, parsimony-based tree fitting is essentially the same as DIVA with ancestral distributions restricted to two unit areas.
We analyzed Southern Hemisphere biogeography based on the published phylogenies of 19 plant and 54 nonmarine animal taxa. In total, the phylogenies included 1,393 terminals, which makes this the largest biogeographic analysis of the Southern Hemisphere attempted so far. We use parsimony-based tree fitting in conjunction with permutation tests to examine the relative role played by dispersal and vicariance in shaping the Southern Hemisphere biota. In particular, we wanted to test to what extent southern biogeographic patterns fit the breakup sequence of Gondwana and to search for the existence of concordant dispersal patterns in the southern flora and fauna.
Geological Scenario
There are several conflicting hypotheses on earth history that explain the evolution of continents by continental drift in the Mesozoic and Cenozoic. The sequential breakup of the supercontinent Pangaea (Scotese et al., 1988) is the most accepted and the one followed here, but the expanding earth theory (Shields, 1979) and the idea of the lost supercontinent of Pacifica (Kamp, 1980) have also been considered by biogeographers in explaining biogeographic patterns incongruent with the Pangaea model (Crisci et al., 1991; Seberg, 1991).
The following account is a synthesis of ideas presented by Scotese et al. (1988), Veevers (1991), Veevers et al. (1991), Lawver et al. (1992), and McLoughlin (2001) (see http://www.ebc.uu/systzoo/staff/sanmartin.html for a more complete account). They are summarized in terms of a geological area cladogram in Figure 1.
During the Triassic, Gondwana formed the southern portion of the Pangaea supercontinent. The climate of Gondwana was not uniform, and paleobotanists and biogeographers usually recognize two climatically different biotic provinces: a Northern Tropical Gondwana province, including northern South America, Africa, Madagascar, India, northern Australia, and New Guinea, with a tropical climate; and a Southern Temperate Gondwana province, including southern South America, Australia, Antarctica, New Zealand, and New Caledonia (Craw, 1982; Karol et al., 2000) and southern temperate Africa (Amorim and Tozoni, 1994).
Gondwana started to break up in the Jurassic (165–150 MY ago [MYA]), when rifting began between India and Australia–east Antarctica about 165 MYA. Shortly after, the continental block of Madagascar and India, which was then adjacent to Somalia, broke away from Africa and began moving southeast, attaining its present position in front of Mozambique in the early Cretaceous, 121 MYA (Rabinowitz et al., 1983). India separated from Madagascar in the late Cretaceous (88–84 MYA), with the opening of the Mascarene Basin (Storey et al., 1995), and began drifting northward, eventually to collide with Asia at about 50 MYA. South America began to separate from Africa in the Early Cretaceous (135 MYA) with the opening of the South Atlantic Ocean at the latitude of Argentina and Chile. Northern South America and Africa remained connected until the mid-late Cretaceous (110–95 MYA), when a transform fault opened between Brazil and Guinea. As a result, Africa started drifting northeast and collided with Eurasia in the Paleocene (60 MYA), whereas southern South America drifted southwest into contact with Antarctica.
New Zealand, Australia, South America, and Antarctica remained connected until the late Cretaceous: east Antarctica was adjacent to southern Australia, and New Zealand and southern South America were in contact with west Antarctica (Fig. 2a). Antarctica had at that time a warm-temperate climate. New Zealand was the first to break away from Antarctica in the late Cretaceous. About 80 MYA the continental block of Tasmantis, including New Zealand and New Caledonia, broke away from west Antarctica and moved northwest, opening the Tasman Sea (Fig. 2b). During most of the Paleogene, New Zealand and New Caledonia were progressively submerged under a marine transgression, with probably > 80% of New Zealand's present area beneath water by the mid-Oligocene (Brothers and Lillie, 1988; Cooper and Milliner, 1993). McLoughlin (2001), however, suggested that New Zealand was not completely submerged but was divided into smaller islands intermittently separated by shallow seas.
|
New Zealand and New Caledonia were finally separated in the mid-Tertiary (40–30 MYA), when the Norfolk Ridge foundered, opening the New Caledonian Basin (Walley and Ross, 1991). Australia and South America remained in contact across Antarctica until the Eocene. Australia began to separate from Antarctica in the late Cretaceous (90 MYA), but both continents remained in contact along Tasmania, and complete separation did not occur until the late Eocene (35 MYA) with the opening of the South Tasman Sea. Some authors, however, have proposed an earlier biotic separation (50 MYA), when a shallow marine seaway formed between Australia and Antarctica (Woodburne and Case, 1996). Southern South America and Antarctica remained in contact through the Antarctic Peninsula until the Oligocene (30–28 MYA), when the Drake Passage opened between these continents, allowing the establishment of the Antarctic Circumpolar Current and the onset of the first Antarctic glaciation. Following its separation from Antarctica, Australia began to drift rapidly toward Asia. New Guinea was then joined to the northern margin of the Australian Plate, although only the southern margin of New Guinea was emergent at that time. The collision of the Australian and Pacific plates in the Oligocene (30 MYA) initiated the tectonic uplift of New Guinea, but by the early Miocene, much of southern New Guinea was again submerged. Subsequent episodes of uplift began in the Miocene after the collision of the Australian and Asian plates and led to the accretion of numerous terranes to the northern margin of New Guinea (Pigram and Davies, 1987). A second phase of collision between Australia and the Asian Plate in the late Miocene–Pliocene (10 MYA) resulted in the rise of mountain chains in Southeast Asia, northern Australia, New Guinea, and New Zealand (e.g., the Southern Alps).
Biogeographic Patterns
As a background for the discussion of distribution patterns, we identified the most commonly discussed area relationships in the Southern Hemisphere biogeographic literature (Fig. 3, Table 1); several of these patterns previously had not been named or studied in a comparative framework. The most famous is the southern Gondwana pattern (SGP; Fig. 3a), which refers to Brundin's (1966) classic pattern based on his studies of chironomid midges (Diptera). This pattern, found in many other animal groups, shows African taxa diverging basally, followed by a New Zealand clade and then a southern South America–Australian clade. The relative position of these areas in the cladogram (Fig. 3a) is congruent with the sequential breakup of these continents from the Gondwanan landmass (Fig. 1): ((Africa...)(New Zealand(southern South America, Australia))). In some groups such as marsupials (Krajewski et al., 2000), New Guinea appears as sister area to Australia, a placement supported by geological evidence (Fig. 1).
|
|
The SGP is usually explained as the result of the sequential breakup of southern temperate Gondwana during the Mesozoic–Tertiary coupled with extinction or primitive absence in some areas (Brundin, 1966; Edmunds, 1981; Cranston and Edward, 1992). The tropical Gondwana pattern (TGP; Fig. 3e) is also explained by vicariance: it follows the breakup sequence of the western part of the northern Gondwana province (Amorim and Tozoni, 1994). In contrast, the plant southern pattern (PSP; Fig. 3b) and the inverted southern pattern (ISP; Fig. 3c) are incongruent with the geological scenario (Fig. 1) and are typically explained by dispersal. The northern gondwana pattern (NGP; Fig. 3d) has been explained by either vicariance or dispersal. It shows a track between the landmasses that once formed northeast tropical Gondwana and is therefore complementary to the TGP. Some groups show a combination of the northern and tropical Gondwana patterns (e.g., Aplocheiloidei) or an NGP in which northern South America branches close to Africa (e.g., Eucnemidae). In most cases, however, groups support either one pattern or the other. The transamerican pattern (AMP; Fig. 3f) follows earlier biogeographic studies (Crisci et al., 1991; Lopretto and Morrone, 1998), which showed South America divided into two biotic provinces with different biogeographic affinities.
| Material and Methods |
|---|
|
|
|---|
Selection of Groups for Analysis
We searched the literature, mainly based on tips from colleagues, for data on groups meeting the following criteria: (1) monophyletic, distributed almost exclusively in the Southern Hemisphere, and claimed to be of Gondwanan origin; (2) represented in at least two (preferably three or more) of the areas studied; (3) including four or more species; and (4) represented by a well-resolved phylogeny with no more than two trichotomies or one tetratomy (unless the extra polytomy occurred within a clade with the same distribution, in which case the clade was summarized in one taxon or one of the possible topologies was arbitrarily chosen for the analysis). In several cases, phylogenies were not resolved to the species level (genera or species groups as terminals) or did not include all the species in the group (i.e., taxon sampling). The Biological Abstracts and Zoological Record databases were carefully searched to ensure that the included phylogenies were the most recent ones for the selected groups. Groups were separated into animal (54 phylogenies) and plant (19 phylogenies) data sets. Insects dominate the animal data set (Fig. 4a) and were also analyzed separately from the other animals. The most well-represented group in the animal data set is the beetle family Eucnemidae (Muona, 1991), with 22 of 54 phylogenies but only 17% of the total number of animal species. The large number of eucnemid groups potentially biased some analyses, which were therefore run both with and without eucnemids included. Vertebrates were represented by iguanid lizards and classic Gondwanan groups such as marsupials, ratites, and cichlid and galaxiid fishes. The plant data set (Fig. 4b) included all main groups in the angiosperm classification: monocots, eudicots, and basal angiosperms (APG, 2003). Following Linder and Crisp (1995), we also included Cyttaria, a fungus parasitic on Nothofagus, in the plant data set because despite being phylogenetically closer to animals, fungi and plants probably have similar means of dispersal (spores/seeds). Table 2 lists the selected groups and references and the oldest appearance of the group in the fossil record. The complete data sets (Appendix A) and supporting references for fossil occurrences and reference taxa (Appendix B) are available at the Systematic Biology website (http://systematicbiology.org/). Fossil evidence indicates that many of the selected groups are or could be old enough to display Gondwanan patterns (Table 2; Appendix B).
|
|
It has been customary in biogeographic studies to exclude groups occurring in less than three areas because they are not considered to provide evidence about area relationships. However, in event-based tree fitting, groups present in two areas do provide potential evidence concerning both dispersal patterns and vicariance events. Therefore, we included in the animal data set three groups that are present in only two areas: Arachnochampinae (NZ-AUS), Botryocladius (SSA-AUS), and characiform fishes (NSA-AFR) (Table 2). These groups have in common that they present a repeated pattern of sister-group relationships between the two areas they occupy. For these particular groups, a vicariance explanation is also supported in the original studies because either the fossil record or the sequence divergences between the taxa indicate that the group originated before the geological separation of the two landmasses (Ortí and Meyer, 1997; Cranston and Edward, 1999).
Not all potential groups are included in our sample. It was brought to our attention after the analyses were completed that the austral biota was undersampled and that there are data for several more Australian–New Zealand plant groups that we were not aware of. Had these groups been included, however, we predict that they would have strengthened the plant dispersal patterns found in the present study. An even more important sampling bias, which is very difficult to correct for, is that certain regions and certain organism groups have been extensively studied than others and therefore are significantly overrepresented in our data sets.
Delimitation of Areas
The areas used in the analysis correspond to historically persistent Gondwanan landmasses according to paleogeographic reconstructions. Composite areas, such as New Zealand (McLoughlin, 2001), were treated as single landmasses because, following accretion, they broke away from Gondwana as one unit. Small landmasses or oceanic archipelagos were included within a larger area when they shared a similar biota or recent land connections during the Pleistocene. Nine areas were finally identified.
- AFR: Africa, excluding the region north of the Saharan belt. We consider subsaharan Africa a single unit because the division between tropical and temperate regions was often not clear from the distribution of the terminal taxa in many of our study groups. In other groups (Karol et al., 2000), the fossil record indicates that the African taxa were part of the tropical African biota in the early Cretaceous but became restricted to South Africa in more recent times. We would have needed more African taxa to be able to differentiate between the histories of tropical Central Africa and temperate South Africa, if they are different.
- MAD: Madagascar, including the Seychelles and the near Indian Ocean islands (Reunion, Mauritius, etc) with faunistically similar biotas (Johanson, 1998).
- IND: India, including Nepal, Tibet, and Sri Lanka. Sri Lanka was connected to India during the Pleistocene glaciations (McLoughlin, 2001).
- AUS: Australia and Tasmania. Tasmania was connected to Australia until 9,000–6,500 years ago (McLoughlin, 2001).
- NZ: New Zealand, including the Kermadec and the subantarctic islands, which lie on the same continental block as New Zealand and probably have been colonized by New Zealand's biota via long-distance dispersal (Pole, 1994).
- NC: New Caledonia.
- NG: New Guinea, including the Solomon and New Hebrides islands, which show close biotic affinities (Raven and Axelrod, 1972; Muona, 1991).
Following previous biogeographic studies (Crisci et al., 1991), we divided South America into two separate regions.
- SSA: southern South America, defined as the southern temperate region (Argentina, Chile, Paraguay, Bolivia, and the Falkland and Antarctic islands), and north-central South America west of the Andes.
- NSA: northern South America, defined as north-central South America east of the Andes (Brazil, Colombia, and Venezuela).
Three other areas were recognized but considered as outside occurrences (not part of Gondwana).
- SEA: Southeast Asia, including the Malaysian Peninsula, Philippines, Sumatra, Borneo, and the Inner Banda Arc. Southeast Asia was formed by the accretion of numerous small terranes that rifted away from Gondwana and accreted to Asia at different times during the Mesozoic–Cenozoic (Metcalfe, 1996). The region thus has a composite biogeographic history (Burret et al., 1991).
- SWP: Southwest Pacific. In our analyses, this area includes the Melanesian archipelagos formed after the Oligocene collision of the Pacific and Australian plates (Fiji, Samoa, Tonga, and Vanuatu), and Hawaii but not New Zealand and New Caledonia, which have always been part of the Gondwanan continent.
- HOL: The Holarctic region comprises the Palearctic (Europe, Asia, and North Africa) and Nearctic (North America and Mexico) infraregions. Central American distributions usually have been considered Nearctic except when the group is also present in northern South America, in which case they are coded as NSA.
Parsimony-Based Tree Fitting
For all analyses, we used parsimony-based tree fitting under the four-event model (Page, 1995) as implemented in the program TreeFitter 1.2 (Ronquist, 2002b). The general methodology, also known as tree reconciliation and implemented in TreeMap (Page, 1995; Page and Charleston, 1998), has been used widely in coevolutionary studies (e.g., Page, 2002, and references therein) but is relatively new to biogeographic analyses (but see Sanmartín et al., 2001). We used several important but not widely employed improvements of the original method (Ronquist, 1998; Sanmartín et al., 2001; Ronquist, 2002a; Sanmartín and Ronquist, 2002), described here (see Ronquist, 2002a, for more detailed coverage).
In the biogeographic context, the four events in the model are (1) vicariance: allopatric speciation caused by the origin of a dispersal barrier affecting many organisms simultaneously; (2) duplication (speciation within an area): usually allopatric and associated with a local or temporary dispersal barrier within an area; (3) dispersal: occurring between isolated areas and associated with speciation; and (4) extinction: leading to the disappearance of a lineage from an area where it is predicted to occur. The intricacies of this model and its relation to alternative models are explained in detail elsewhere (Ronquist, 1998, 2002a).
Once each event in the model is associated with a cost, it is possible to fit an organism phylogeny and the associated geographic distributions of the terminals to any given area cladogram (Fig. 5a), producing a reconstruction of the biogeographic history of the organisms according to the area cladogram. Each reconstruction specifies a set of ancestral distributions and biogeographic events that could have produced the observed distribution patterns on the given area cladogram. The reconstruction can most easily be illustrated using a trackogram displaying the organism phylogeny on top of the area cladogram, with symbols denoting the four kinds of events (Fig. 5b). The cost of the reconstruction is simply the sum of the costs of the implied events. Of particular interest in event-based tree fitting is the reconstruction with the minimum cost, the most-parsimonious reconstruction (MPR).
|
Searching for the optimal area cladogram
The cost of the MPR fitting one (or more) group(s) of organisms to a given area cladogram represents a natural measure of the fit of that area cladogram to the organism data. The area cladogram with the lowest cost, i.e., the most-parsimonious area cladogram, best explains the taxon distributions in the phylogeny. This optimal area cladogram can be found, for instance, by explicit enumeration of all possible area cladograms or by heuristic search for the best area cladogram using algorithms similar to those used in ordinary parsimony analysis of phylogeny (e.g., stepwise addition followed by branch swapping). In our searches for optimal area cladograms, we used different hold settings in TreeFitter (determining how many trees are kept in each step of the stepwise addition algorithm) and different swapping neighborhoods (reconnection limits), determining how exhaustive the branch swapping is (Ronquist, 2002b).
Randomization tests
Any aspect of the MPR calculated from the observed data can be tested by comparing it to MPRs derived from random data sets. Depending on the property of interest, the random data sets can have the terminals in the organism phylogeny, the terminals in the area cladogram, or both randomly permuted (Siddall, 1996; Ronquist, 2002a). Alternatively, one can randomize the topology of the area cladogram or the organism phylogeny instead of permuting the terminals (Page, 1995). However, this approach can be slightly biased if the random trees are not generated in the appropriate way (Ronquist, 2002a). Therefore, we consistently used permutation of terminals.
To examine whether organism distributions are phylogenetically constrained in general according to a given area cladogram, we permuted both areas and organism terminals. The null hypothesis of no phylogenetic constraint was rejected when few of the permutations fell below the observed MPR cost. For instance, the result was considered significant at the 0.05 level when < 5% of the permutations had a cost less than or equal to that for the observed data. Often, one is interested in a particular event or process, in which case other types of permutations are more appropriate. For instance, to examine whether an area cladogram is successful in predicting vicariance events, touching the organism phylogenies will only cause unnecessary noise. To obtain maximum sensitivity in these tests, we used permutations of the terminals in the area cladogram when focusing on vicariance and extinction events and permutation of the organism terminals when focusing on dispersal and duplication events.
When the area cladogram is inferred from the organism phylogenies, it is inappropriate to randomize the area cladogram because the error rate will be inflated (Siddall, 1996). Instead, the organism phylogenies must be randomized and the best area cladogram inferred for each random data set. This approach was taken here for the area cladograms based on organism data and not on geological evidence.
Event costs
Setting the event costs is the central problem in parsimony-based tree fitting. Using the same cost for all events is not a viable approach because then we can guarantee that no events are fitted to internal nodes and internal branches in the area cladogram (Ronquist and Nylin, 1990). As a possible solution, Page (1995) suggested maximizing codivergence events (MC), which is equivalent to maximizing vicariance events in biogeographic problems. Although MC has been popular in coevolutionary studies, it has serious drawbacks. In particular, MC is sensitive to dispersal and duplication events (Ronquist, 2002a). Even a moderate number of such events can make it impossible to identify a strong underlying vicariance pattern. Therefore, it is not possible to use MC in inferring predominantly dispersal patterns or mixed vicariance–dispersal patterns (Ronquist, 2002a).
Assigning a positive cost to all events is more in line with the parsimony principle in general and produces a more powerful analytical tool. However, it does leave the difficulty of setting the exact costs to be used. Ronquist (2002a) proposed a procedure for selecting the optimal event costs based on maximizing the chances of finding significance in permutation tests focusing on total cost. The individual event costs are simply varied until a cost combination is found that minimizes the probability of permuted data sets having a total cost lower than that of the observed data (Fig. 5c). Using this approach in examining a variety of hypothetical data sets, Ronquist (2002a) found a set of event costs that performed well under a wide range of conditions. This solution has an arbitrarily low cost for vicariance and duplication events and a higher cost for dispersal than for extinction events. It performs well because both vicariance and duplication generate phylogenetically constrained distribution patterns, whereas dispersal and extinction do not. We used the default cost assignments of TreeFitter 1.2 (0.01 for vicariance and duplication, 1.0 for extinction, 2.0 for dispersal), which conform to this scheme. The primary results given here (fitting the entire animal or plant data set to the geological area cladogram) indicate that these event costs are in the optimal range for our data, given 1,000 permutations of organism area cladogram terminals.
An alternative to setting specific event costs is to enumerate all possible solutions under a range of event costs (the jungle approach; Charleston, 1998; see also Ronquist, 2002a). This approach is still too computationally complex to be employed in a large study such as the one attempted here.
Algorithms
When the event costs are fixed, it is possible to use fast polynomial algorithms for finding the upper and lower bounds on the total cost of the true MPR (Ronquist, 1995, 1998, 2002a). In practice, lower-bound solutions are almost always identical to the true MPR (differences may arise because of incompatible dispersals; see Ronquist, 1995). Here, we consistently used the lower-bound algorithm implemented in TreeFitter 1.2.
Widespread terminals
Widespread taxa (terminal taxa distributed in more than one area) pose a problem in biogeographic reconstruction because they introduce ambiguity in the data set (Morrone and Crisci, 1995). This problem has traditionally been dealt with using Assumptions 0, 1, and 2 of Nelson and Platnick (1981), but these assumptions are inapplicable to event-based methods for several reasons. Each of them specifies a set of possible solutions but not the relative cost of these solutions. This approach may result in nonoverlapping solution sets (Van Veller et al., 1999) and other problems, but more critically it is incompatible with the fact that all area cladograms must have a relative cost in an event-based method (Ronquist, 2002a). Therefore, Sanmartín and Ronquist (2002; see also Sanmartín et al., 2001; Ronquist, 2002a) introduced three novel event-based options to treat widespread taxa. Under the Recent option, the widespread distribution is assumed to be the result of recent dispersal. Thus, only one of the areas in the distribution is considered ancestral. The Ancient option assumes that the widespread distribution is the result of a failure to vicariate and explains any mismatch between the distribution and the area cladogram by extinction. Under this option, all of the areas in the distribution are considered ancestral. The Free option treats the widespread terminal as an unresolved higher taxon consisting of one lineage occurring in each area and permits any combination of events and any resolution of the terminal polytomy in explaining the widespread distribution. Each of these options can be easily implemented in algorithmic form (Sanmartín and Ronquist, 2002) and are available in TreeFitter (Ronquist, 2002b).
Comparing the event-based options with the traditional assumptions is like comparing apples with oranges, but there are some similarities. The Ancient option is, in some sense, a combination of Assumption 0 and Assumption 1, and the Free option is reminiscent of Assumption 2. The Recent option, however, is a different concept altogether; it is similar to the way polymorphic taxa are usually treated in ordinary parsimony analysis of phylogeny. Using empirical data, Sanmartín et al. (2001) showed that the Recent option, unlike the Free and Ancient options, forces vicariance events onto the ancestral nodes in the phylogeny, and this option is therefore more powerful in separating phylogenetically constrained distribution patterns from random data in permutation tests. For this reason, we consistently used the Recent option in treating widespread terminals. When examining dispersal patterns, we excluded dispersals within widespread terminals and only looked at inferred ancestral dispersals, unless explicitly stated otherwise. By focusing exclusively on ancestral events, we reduced the inherent problem of comparing phylogenies that have widespread higher taxa as terminals with phylogenies that have endemic species as terminals.
Analyses
Best area cladogram
We first searched for the best area cladogram fitting the animal and the plant phylogenies, respectively. We used the heuristic tree search of TreeFitter because the exhaustive search was too time consuming. An initial tree was constructed by stepwise addition, holding 100 trees in each step, followed by tree bisection-reconnection (TBR) branch swapping with a reconnection limit (neighborhood) of 20 nodes (essentially identical to a full standard TBR search). Using the same approach, we searched for the best area cladogram for each of three different subgroups of the animal data set (animals excluding insects, insects alone, and insects excluding eucnemids).
Randomization tests implemented in TreeFitter 1.2 were used to assess the statistical significance of the fit between the best area cladograms and the organism phylogenies (Ronquist, 2002a). We generated 1,000 random data sets for which the terminal distributions in the organism phylogenies had been permuted and inferred the best area cladogram for each (hold = 5, neighborhood = 20). We then compared the total cost and the number of vicariance events obtained when fitting each of the random data sets to its best area cladogram, with the corresponding values obtained when fitting the observed data to its best area cladogram obtained under the same search conditions. The significance value (P) was calculated as the percentage of random data sets with a lower cost or a higher number of vicariance events than the original data set.
Biogeographic patterns
We also searched for the optimal area cladogram for each individual group in the data sets (Table 2) and then compared that area cladogram with each of the biogeographic patterns (Fig. 3). To calculate the degree of congruence (similarity) between the area cladogram and the biogeographic patterns, we used a modification of the method of Robinson and Foulds (1981) for calculating the metric distance between two phylogenetic trees (see Fig. 6).
|
The biogeographic pattern with the highest positive score for a particular organism group was assigned a support value of 1.0 for that group, except when the score was lower than 0.1, in which case the group was considered as not supporting any biogeographic pattern. When a group supported two compatible biogeographic patterns such as TGP and NGP (e.g., Aplocheiloidei, Table 2), the pattern with the highest score received a support value of 1.0 and the pattern with the lowest score received a value of 0.5. The total support for each biogeographic pattern was computed by summing across all individual groups in the animal and plant data sets separately.
Fit to the geological scenario
To address whether the biogeographic history of the Southern Hemisphere reflects the geological breakup of Gondwana, we fitted the geological area cladogram (Fig. 1) to each of the organism phylogenies and used randomization tests in TreeFitter to assess the significance of the fit. When the group fitted the geological scenario significantly (P < 0.05), we calculated the frequency of each of the four different biogeographic events required to explain the fit (vicariance, duplication, dispersal, and extinction) and tested whether the observed frequency of the events departed from what would be expected by chance, as indicated by 1,000 randomly permuted data sets.
Dispersal analysis
We also examined the role played by dispersal in shaping Southern Hemisphere biogeographic patterns by fitting the organism phylogenies to the geological area cladogram, i.e., we looked at the dispersals remaining after geologically predicted events, such as vicariance, had been removed.
First, we compared the overall frequency of dispersal events in animals and plants. For this analysis, we divided dispersal events into terminal dispersals, occurring within widespread terminal taxa, and ancestral dispersals, associated with interior nodes in the organism phylogenies. Ancestral and terminal dispersals were considered both separately and combined. To test whether the frequency of these dispersals differed significantly between the animal and plant data sets, we used a permutation test. The distribution of expected values was obtained by generating 100 data sets in which each group was categorized as animal or plant at random while maintaining the original number of groups in each class (54 and 19), as in a partition homogeneity test. The null hypothesis, i.e., no association between the frequency of dispersals and the type of organism involved (plant or animal) would be rejected (P < 0.05) if the difference in dispersal frequency between the original animal and plant data sets were larger than the difference in > 95 % of the random data sets.
Second, we calculated the frequency of dispersal events between the major southern landmasses and tested the significance of the dispersal patterns using the randomization tests in TreeFitter. In these tests, the observed frequency of the predicted dispersal patterns was compared with that expected by chance by using a modified chi-square test in which the expected frequency and the reference distribution were both calculated from 100 randomly permuted data sets. Permuting the terminals of the area cladogram is not possible in these tests because such permutations would change the set of biogeographic events. Therefore, only the terminals in the organism phylogenies were permuted.
Before determining the direction of dispersal asymmetries, observed values were normalized against the number of expected events from the randomizations. This step is important because not all the southern landmasses are equally represented in our data sets (Figs. 4c, 4d), and this inequality would heavily bias the number of dispersal events (Wares et al., 2001). The corrected frequencies compared here represent the number of dispersal events in each direction that we would have seen if all areas had been equally well represented in the data set.
| Results |
|---|
|
|
|---|
Area Cladograms
Figure 7 shows the optimal area cladograms for the animal and plant data sets found using the heuristic search in TreeFitter. All cladograms were significant in terms of overall cost (P < 0.01 in all cases) and had a higher number of vicariance events than expected by chance (P = 0.01 for the plant data set; P < 0.01 for all others). The animal cladogram (Fig. 7a) has a pattern in agreement with the SGP (Fig. 3a), except for the position of Africa, which appears closely related to other areas: (Africa, ...)((New Zealand(S. South America, Australia))). The same area cladogram was obtained when insects were analyzed separately, but when the Eucnemidae were excluded from the insect data set the SGP became complete in that it now also included Africa (Fig. 7b). Separate analysis of the animal groups excluding insects (Fig. 7c) also reflected the SGP, in this case including New Guinea as sister group to Australia: (New Zealand (S. South America (Australia, New Guinea))). In contrast to all animal area cladograms, the plant area cladogram (Fig. 7d) supports the PSP (Fig. 3b): (S. South America(New Zealand, Australia)).
|
Both the animal and plant area cladograms support a hybrid origin for the South American biota (Crisci et al., 1991): northern South America is consistently grouped with the Holarctic region (mainly North American distributions), whereas southern South America appears more closely related to Australia or New Zealand (Fig. 7). None of the animal or plant area cladograms show the sister-group relationship between New Zealand and New Caledonia predicted by the geological scenario (Fig. 1); instead, biotic relationships separate these areas in the area cladogram.
Biogeographic Patterns
Figure 8 shows a comparison between the animal and plant data sets in terms of the number of groups supporting each of the six major biogeographic patterns (cf. Fig. 3). The animal data set (Fig. 8a) is dominated by the SGP and by the NGP. Support for the NGP, however, is almost exclusively due to the eucnemid beetles. They constitute a small percentage of the animal data set in terms of number of species (17%), but they represent a significant proportion of the groups (22 of 54; 41%). When the Eucnemidae are removed from the analysis, the SGP becomes the overwhelmingly dominant pattern in animals and in insects (Fig. 8b). Similarly, most of the noninsect animal groups support the SGP or the TGP (Fig. 8c). In contrast, the majority of groups in the plant data set support either the PSP or the ISP, with almost no support for the SGP (Fig. 8d). Several groups did not support any clear pattern (multiple area cladograms) or supported a pattern different from those represented in Figure 3 (cf. Table 2).
|
Fit to the Geological Scenario
Both the animal and plant data sets display highly significant, phylogenetically conserved distribution patterns (P < 0.001) when fit to the geological scenario (Fig. 1). However, the processes responsible for the distribution patterns differ between plants and animals (Table 3). The animal data set as a whole shows a higher frequency of vicariance, extinction, and duplication events than expected by chance (P < 0.05), whereas dispersals are significantly rare. The plant data set also exhibits a higher frequency of duplication and a lower frequency of dispersal than expected, but the frequencies of vicariance and extinction events do not depart significantly from expected values. This difference can also be observed when individual groups are fitted to the geological area cladogram: with the exception of Strelitziaceae, none of the plant groups shows a significant frequency of vicariance or extinction, whereas many of the animal groups do (Table 3).
|
Dispersal Analysis
Comparison of the frequency of dispersal events between the animal and plant data sets shows that the frequency of terminal dispersals is significantly higher in plants than in animals (P = 0.02), but the frequencies of ancestral and total dispersals, although higher in plants, do not differ from expected values (P = 0.42 and 0.96, respectively). Thus, there is more dispersal in the plant data set, but it appears to be solely due to dispersal within widespread terminals.
Table 4 shows comparisons of the observed and expected frequencies of dispersal between the southern landmasses in the animal, insect, and plant data sets (the number of events in the noninsect animal data set was too small for statistical comparisons). Figures 9 and 10 represent the same comparisons for the animal and plant data sets, but the dispersal frequencies are reported as a ratio against the number of expected events from the randomizations to correct for the unequal representation of the areas in the data set. We found several significant concordant dispersal patterns in the data sets, but these patterns were different in animals and plants (the insect data set showed patterns identical to those observed for animals; see Table 4).
|
|
|
In animals, trans-Antarctic dispersal between Australia and southern South America (AUS
SSA) is significantly more frequent than any of the other dispersal events involving the austral landmasses (P = 0.01, Table 4). In contrast, trans-Tasman dispersal (NZ
AUS) is significantly more frequent than trans-Antarctic (AUS
SSA) and trans-Pacific (NZ
SSA) dispersals in plants (P = 0.01). New Zealand is more frequently the source area of trans-Tasman dispersal than is Australia in animals (Fig. 9), whereas the opposite pattern is found in plants (Fig. 10), but none of these differences are significant (P = 0.27, Table 4). There is, however, highly significant directional asymmetry in trans-Pacific biotic exchange (NZ
SSA), but it follows opposite directions in animals and plants: New Zealand is the most frequent source area of dispersal in the animal data set (P = 0.01, Fig. 9), whereas southern South America is almost the only source of dispersal in the plant data set (P < 0.01, Fig. 10).
Trans-American exchange (SSA
NSA) in animals is significantly lower than the biotic exchange between northern South America and North America (NSA
HOL) (P < 0.01), and southern South America seems to be the main source area of these dispersals (P = 0.04). This comparison was not significant in plants.
Not surprisingly, dispersal between Madagascar and Africa (AFR
MAD) is significantly more frequent than trans-Atlantic dispersal between northern South America and Africa (AFR
NSA) in animals (P = 0.01) and plants (P < 0.01). Northern South America is more frequently the source area of trans-Atlantic dispersal than is Africa (AFR
NSA) in plants (P < 0.01), but the values are very small.
| Discussion |
|---|
|
|
|---|
Hierarchical Vicariance Patterns
Sanmartín et al. (2001) found that a maximum-vicariance analysis of 57 Holarctic animal groups failed to support any hierarchical area cladogram, and they concluded that this occurs because the northern landmasses were joined to each other in various combinations over time, creating a reticulate biogeographic pattern. In contrast, if current distribution patterns in the Southern Hemisphere were shaped to a large extent by the sequential breakup of Gondwana and the associated vicariance events, we should be able to retrieve a strongly supported hierarchical pattern reflecting the geological area cladogram (Fig. 1).
Our results indicate that the breakup of Gondwana has indeed played an important role in molding the biogeography of our animal groups (presumably Mesozoic). The best area cladograms for the animal data set as a whole and for different animal subgroups (Figs. 7a–c) all contain important components of the geological area cladogram (Fig. 1). Geologically predicted vicariance events are significantly more common than expected by chance in animals (Table 3), and many of the species-rich animal groups support the SGP (Fig. 8b), which agrees well with the hypothesized geological breakup of southern temperate Gondwana (Fig. 1).
The congruence between zoogeographic patterns and southern temperate Gondwana breakup implies that many of the animal groups studied here could be old. Although the precise timing of the loss of land connections among southern Gondwana landmasses is unclear, most authors agree on a date of 80 MYA for complete continental separation between New Zealand and Australia, and a date of 35 MYA is generally agreed for the final breakup of terrestrial connections between South America and Australia (Veevers et al., 1991), although some authors maintain that biotic exchange was terminated already in the early Eocene, 50 MYA (Woodburne and Case, 1996).
Somewhat surprisingly, our results suggest that biogeographic patterns in our plants (of presumably Mesozoic or early Cenozoic age) have not been significantly influenced by the Gondwanan breakup. The best plant area cladogram supports the PSP (S. South America(New Zealand, Australia)), which has often been found in plant biogeographic studies (Seberg, 1991; Linder and Crisp, 1995) but is in conflict with continental fragmentation patterns (Fig. 1). The dominance among individual plant groups of the PSP and the ISP (Fig. 8d), two patterns typically explained by dispersal (Table 1), further corroborates the departure of plant phylogenies from geological vicariance. When the plant data set was fitted to the geological area cladogram, none of the plant groups (except Strelitziaceae) showed a higher frequency of vicariance and extinction events than expected by chance. The same result was obtained for the entire plant data set. If plants had responded to the breakup of Gondwana with speciation (vicariance) followed by local extinction in some areas, both vicariance and extinction events should be more common than expected by chance.
Why do plant distributions show no traces of Gondwanan breakup? There seems to be two possibilities: either the studied plant groups are too young to have been affected or they once fragmented in response to Gondwanan splits but the original patterns were subsequently lost because of dispersal and extinction events. Several lines of evidence support the latter hypothesis for many of the plant groups studied here. First, fossil evidence (Table 2, Appendix B) indicates that many plant groups are old enough to have been influenced by at least the more recent Gondwanan geological events (Fig. 1). Second, for at least one group, the southern beeches (Nothofagus), adding fossil distributions makes the biogeographic pattern congruent with the breakup sequence of Gondwana (Manos, 1997; Swenson et al., 2001). Thus, the distributions of southern beeches were initially shaped by geologically induced vicariance events, but subsequent extinctions have resulted in the extant species showing a pattern in conflict with geology, the PSP (Fig. 3b). Third, colonization of new areas might be expected to be successful more often in plants than in animals, given that plants have a resistant dispersal stage (seeds/spores), that should facilitate long-distance dispersal. Not surprisingly, there is also more specific evidence for some plant groups that suggests that recent dispersal has indeed been superimposed onto an ancient vicariance pattern, for example in Atherospermataceae (Renner et al., 2000). We found significantly more dispersal events among plant terminals than among animal terminals, consistent with the hypothesis that vicariance patterns in plants have been obscured by later dispersal. However, the number of terminal dispersals is also influenced by the proportion of terminals that are higher taxa, which tend to be more widespread and therefore imply more dispersal events. In our data set, there were more higher-taxon terminals in plants than in animals (12% vs. 5%), so it is possible that the difference in terminal dispersals was caused entirely due to this effect.
If plant biogeography has not been shaped by Gondwanan breakup, why do we find strongly supported hierarchical patterns, such as the PSP, among plants? It is difficult to exclude the possibility of other sets of geological events underlying these alternative hierarchical patterns, but a more likely explanation is probably that they are the result of recurring combinations of biotic events other than geologically induced vicariance. Consider, for instance, the close relationship between New Zealand and Australia in the PSP (Fig. 3b). This type of pattern could easily arise in old groups that were originally present on New Zealand, then went extinct there, and finally recolonized New Zealand from Australia. Unless the colonizing lineage was deeply nested within an Australian clade, it would be difficult or impossible to distinguish such a pattern from a true hierarchical vicariance pattern. The same pattern could have arisen in groups that were too young to be in New Zealand at the time of rifting (80 MYA) and instead dispersed there from Australia later. This appears to be the case in, e.g., ratite birds (Van Tuinen et al., 1998; Cooper et al., 2001). Thus, hierarchical biogeographic patterns could be created both by a classical vicariance scenario and by recurring combinations of dispersal and extinction events that generate hierarchical patterns.
Concordant Dispersal Patterns
To be able to infer with any confidence dispersal events from phylogenetic evidence, these events must be rare in comparison with the number of lineage-splitting events. That is, dispersal events must have occurred against a background of phylogenetically constrained events, such as vicariance or local speciation (duplication). This was indeed the case in both animals and plants: dispersals were significantly less common and vicariance or duplication events were significantly more common than expected by chance (Table 3). Once it has been determined that dispersals are phylogenetically constrained, so that the individual events can be inferred with reasonable confidence, it becomes possible to search for general dispersal patterns.
In biogeography, dispersals are often considered to be random, but our analyses reveal strong concordant dispersal patterns in the Southern Hemisphere, both in animals and in plants. In other words, dispersals tend to occur between certain areas and in certain directions.
Before discussing the individual patterns, it is useful to consider the way in which the four kinds of biogeographic events relate to each other when biogeographic patterns are looked upon from different perspectives. For example, consider local speciation (duplication) events. Even in groups with Gondwanan vicariance patterns, our results typically indicate that local speciation was significantly more common than expected (Table 3). However, if the same patterns had been analyzed in more geographic detail, with more groups represented, we might have discovered that some of the local speciation events correspond to vicariance or dispersal events within the broadly defined regions used in the present study. Thus, the high frequency of local speciation seen here undoubtedly reflects to some extent the lack of geographic detail.
A more subtle example of the effect of different perspectives concerns dispersal and vicariance events. A dispersal event is a range expansion across an existing dispersal barrier followed by allopatric speciation along the barrier. Thus, a dispersal event is just a vicariance event that is not predicted by the area cladogram. Under certain conditions, concordant dispersal events, possibly in combination with extinctions, may occur in a sequence that generates a hierarchical pattern, which is then indistinguishable from a sequential vicariance pattern. Fitting organism phylogenies to such a dispersal-generated area cladogram will result in the concordant dispersal events being identified as predicted vicariance events. It is only when the organism phylogenies are fitted to a geologically based area cladogram that these dispersal events can be identified as true dispersals. The situation is analogous to the mapping of ordinary characters onto phylogenetic trees. A feature that is indicated to be a homoplasy (dispersal) on one tree may appear as a synapomorphy (vicariance) on another tree. The best way to accurately separate homoplasies from homologies is to map characters onto the true phylogeny. However, most concordant dispersal patterns do not generate hierarchical distribution patterns (just like sets of convergences rarely conform to a hierarchical pattern), and the identification of them is hence more robust to the area cladogram used. The concordant dispersal patterns discussed here (Figs. 9, 10; Table 4) are all identified against the background of the geological area cladogram (Fig. 1), and this approach should be borne in mind when interpreting them.
Trans-Antarctic exchange
In animals, trans-Antarctic dispersal between Australia and southern South America was significantly more frequent than any of the dispersal events involving New Zealand and the other austral landmasses (Fig. 9; Table 4). This pattern is probably an effect of time and past continental connections. Given that the isolation of New Zealand (80 MYA) occurred before the vicariance between southern South America and Australia (35 MYA), we would expect extinction to have erased the traces of the old dispersal events involving New Zealand, making them rarer than Australian–South American dispersals among the events documented by extant organisms and their phylogenies. The dominance of trans-Antarctic dispersal also emphasizes the importance of the Antarctic route in southern faunal exchange. This route was established during the late Cretaceous–early Tertiary, when the Southern Andes became connected through the Scotia Arc to the Antarctic Peninsula, which in turn was connected to Australia through east Antarctica, and it was not broken until the late Eocene (35 MYA), when the South Tasman Sea opened between Australia and east Antarctica (Lawver et al., 1992). Throughout most of this time period, the connected landmasses had a temperate climate, and Antarctica was covered by an angiosperm-rich flora dominated by Nothofagus. The route thus allowed a long period of biotic exchange between Australia and southern South America. Paleogene fossils of marsupials (Woodburne and Case, 1996) and ratite birds (Van Tuinen et al., 1998) from Antarctica further document the importance of the trans-Antarctic dispersal route.
Trans-Tasman dispersal
The dominance of trans-Tasman dispersal in plants (Fig. 10) is associated with the high frequency of sister-group relationships, (New Zealand, Australia)S. South America), in austral plant groups (Linder and Crisp, 1995). Although the PSP has been sometimes explained by an alternative sequence of breakup events involving the lost Pacifica continent (Crisci et al., 1991; Seberg, 1991), a more common hypothesis is that this pattern is the result of trans-Tasman dispersal between Australia and New Zealand after terrestrial connections were broken 80 MYA (Pole, 1994, 2001). The opening of the Drake Passage between South America and Antarctica in the Oligocene (30 MYA) led to the establishment in the region of the easterly directed West Wind Drift (Veevers, 1991). Historical records show that birds and insects are regularly blown across the Tasman Sea from Australia to New Zealand (Michaux, 1991). Marine dispersal in currents associated with the West Wind Drift is also invoked to explain sister-group relationships between New Zealand–Australian marine organisms (Fell, 1962; Waters et al., 2000a). Pole (1994, 2001) went even further in arguing that most of New Zealand's present flora came by long-distance dispersal from Australia during the late Tertiary because the old Gondwanan lineages, dating back to New Zealand's isolation in the late Cretaceous, went extinct during New Zealand's submergence in the Paleogene (but see McLoughlin, 2001). The direction of the West Wind Drift, however, predicts that dispersal should be significantly more frequent from Australia to New Zealand than in the other direction. Although Australia was the main source of trans-Tasman dispersals in the plant data set (Fig. 10), this difference is not significant in permutation tests. In the animal data set, trans-Tasman dispersal seems to have been against the direction of the prevailing winds (Fig. 9), but the difference is again not significant. More data are needed to document whether the expected directional asymmetry due to the West Wind Drift exists. However, other studies on plants (Linder, 1999; Winkworth et al., 2002) have also documented trans-Tasman dispersal in both directions, contradicting the West Wind Drift hypothesis.
An alternative hypothesis is that much of the trans-Tasman dispersal has been over stepping-stone islands, making the predominant wind direction less important. This dispersal would have involved interconnecting islands and shallow water bridges across the Lord Howe Rise and Norfolk Ridge, which according to some reconstructions remained close to Australia until the Middle Tertiary (Walley and Ross, 1991). A new uplift in the late Tertiary would have permitted biotic dispersal as late as the Miocene (Herzer et al., 1997). However, the existence of these biotic routes is not clear because alternative paleogeographic reconstructions show the Lord Howe Rise and Norfolk Ridge covered by shallow seas since the end of the Cretaceous, 70 MYA (McLoughlin, 2001).
Trans-Pacific dispersal
A sister-group relationship, (southern South America, New Zealand)Australia), in conflict with the geologically predicted vicariance pattern, is often found in the Southern Hemisphere flora and fauna (Swenson and Bremer, 1997; Renner et al., 2000; Winkworth et al., 2002). Although this disjunction has been explained sometimes by a different sequence of vicariance events (Glasby and Alvarez, 1999), a more convincing explanation is that it is the result of long-distance dispersal between New Zealand and South America after New Zealand drifted away from west Antarctica, 80 MYA. Wind dispersal in the South Pacific is presumed to be predominantly eastward, from New Zealand to southern South America, following the direction of the West Wind Drift (Winkworth et al., 2002). The dominance of eastward dispersals in the animal data set, with New Zealand acting as the main source area of trans-Pacific dispersal (Fig. 9), seems to provide evidence of the importance of the West Wind Drift in forming Southern Hemisphere biogeographic patterns. However, the highly directional dispersal asymmetry documented here in plants (Fig. 10) suggests exclusively westward long-distance dispersal (S. South America
New Zealand), against the prevailing direction of winds and ocean currents (see also Winkworth et al., 2002). A possible explanation is that dispersal from southern South America to New Zealand did not involve direct jumps but rather stepping-stone dispersal along the coastline of Antarctica (Swenson and Bremer, 1997; Renner et al., 2000). This hypothesis is supported by the fact that temperate forests existed in the Antarctic coastal areas until at least the Pliocene (Swenson and Bremer, 1997). This dispersal could have been mediated by the west-flowing East Wind Drift, which runs close to the Antarctic coast, or could have followed the West Wind Drift around Antarctica, involving dispersal first to the subantarctic islands (and/or Australia) and from there to New Zealand (Swenson and Bremer, 1997).
South American Biotic Divide
North America and South America became connected during the middle Cretaceous (100 MYA) across the proto-Caribbean archipelago, located further to the west than it is today. This connection was broken in the early Eocene, about 49 MYA (Hay et al., 1999). Significant faunal exchange presumably took place across the Caribbean connection and involved a first invasion of placental and marsupial mammals (Springer et al., 1998). North and South America became connected again during the late Tertiary, first via the Panama Island Arc (15 MYA) and later across the Panama Isthmus in the late Pliocene (3.5 MYA). This connection established a second period of biotic exchange, the Great American Exchange, with the replacement of the southern marsupial fauna by a second invasion of placental mammals. The importance of this faunal exchange is supported by the sister-group relationship between northern South America and the Holarctic (North America) in the animal and insect area cladograms (Figs. 7a, 7b) and the close relationship in the noninsect cladogram (Fig. 7c) and the plant cladogram (Fig. 7d). In contrast, the biotas of northern South America and southern South America do not appear to be closely related. Southern South America appears more closely related to the austral landmasses (Australia, New Zealand) than to northern South America in the plant and all animal area cladograms. The frequency of trans-American dispersal between northern and southern South America is significantly lower than that between northern South America and the Holarctic in the animal and insect data sets (Table 4). This result confirms previous biogeographic results suggesting that the South American biota is formed by a northern tropical component and a southern temperate component, each with different biogeographic affinities (Crisci et al., 1991; Amorim and Tozoni, 1994; Lopretto and Morrone, 1998). It is difficult to explain this biotic divide on geological grounds. Northern and southern South America have always been connected, except during a brief period during the Cretaceous (100–80 MYA) when an epicontinental sea separated both halves. The most likely explanation is ecological: the tropical forests of the Amazonian Basin have probably constituted a long-lasting, powerful barrier to dispersal for temperate organisms. Dispersal of South American animals was significantly more frequent from the southern to the northern region than in the other direction, in contrast to the direction of the Great American Exchange, which is thought to have been predominantly southward. However, the presumed direction of the trans-American exchange is based on fossil evidence from mammals (Woodburne and Case, 1996), whereas most of the (few) South American dispersals are in the insect phylogenies (Table 4), the main component of our animal data set.
Northern Gondwana Pattern
Biogeographic relationships among the landmasses that form the southeastern tropics (Africa, Madagascar, India, Southeast Asia/Southwest Pacific, northern Australia, and New Guinea) are typically explained as the result of recent dispersal along the coasts of the Indian Ocean. These areas share a tropical biota, and exchange between the eastern landmasses has been frequent since the collision of the Australian and Asian plates in the Miocene (Raven and Axelrod, 1972). Long-distance dispersal over water barriers is probably the likeliest explanation for angiosperm groups such as the baobabs (Baum et al., 1998; see Table 2), which radiated long after the separation of Africa and Australia in the early Cretaceous, 135 MYA.
However, for other groups this pattern could be much older. The Kerguelen Plateau, a continental block situated off northeastern India that foundered at the end of the Cretaceous (McLoughlin, 2001), probably acted as a land bridge between Australia–Antarctica and India–Madagascar–Africa during the late Cretaceous, after the separation of Africa and Australia. The existence of this land link is supported by the pan-Gondwanan distribution of many late Cretaceous terrestrial floras (Karol et al., 2000; McLoughlin, 2001) and faunas (Cranston and Hare, 1995; Cracraft, 2001). Macey et al. (2000) propose an even older chain of dispersal events to explain the NGP observed in acrodont lizards (Table 2). Small continental fragments could have drifted off from Gondwana during the late Paleozoic–Mesozoic and accreted to the Eurasian Plate during the Mesozoic–Cenozoic, e.g., India or the small terranes forming Southeast Asia (Metcalfe, 1996). These fragments would have acted as vessels transporting taxa across the Tethys Sea long before the mid-Miocene collision between the Australian and Asian plates.
Tropical Gondwana Pattern
Madagascar diverged from Africa in the early Cretaceous (121 MYA) with the opening of the Somali Basin, whereas Africa and northern South America remained connected until the mid-Cretaceous, about 100 MYA (Fig. 1). However, Africa and Madagascar have always remained relatively close, whereas Africa and tropical South America quickly drifted away from each other once they started to separate. The narrow Mozambique Channel is probably less of a barrier to dispersal for plants and animals than the South Atlantic Ocean. This difference is obvious when comparing vicariance and dispersal patterns. Some of the plant and animal groups (e.g., Strelitziaceae, Aplocheiloid fishes) support in their area cladogram the old vicariance pattern shown in Figure 3e: (Madagascar(Africa, N. South America)). However, the animal area cladogram (Fig. 7a) and one of the two plant area cladograms (Fig. 7d) show a sister-group relationship between Madagascar and Africa, with northern South America closer to the Holarctic. More significantly, the dispersal analysis shows that dispersal between Africa and Madagascar is significantly more frequent than that between Africa and northern South America (Table 4). It is therefore likely that the more recent dispersal pattern between Madagascar and Africa has obscured the basal vicariance pattern. Also, some of the studied groups are probably not old enough to have been present before the vicariance of Madagascar in the early Cretaceous (e.g., cichlid fishes, Table 2). When the direction of dispersals is compared, northern South America appears to have acted as the main source of dispersal in the trans-Atlantic exchange, but the difference is not significant in animals, although it is significant or close to significant in plants and insects, respectively (Table 4). If the asymmetry exists, it could be explained by differential extinction during the post-Miocene aridification, which was more severe in Africa than in northern South America (Farias et al., 1999).
New Zealand and New Caledonia
New Zealand and New Caledonia were linked by continental crust along the Norfolk Ridge from the Permian (McLoughlin, 2001) until the foundering of the Norfolk Ridge and the opening of the New Caledonia Basin in the early Oligocene (30 MYA). This long-lasting connection is presumably reflected in the presence of many common relict groups, which date back to the isolation of New Zealand and New Caledonia from Gondwana in the late Cretaceous (Raven and Axelrod, 1972). One would therefore expect to see the sister-group relationship in the geological area cladogram (Fig. 1) to be reflected in the biogeographic patterns of the southern flora and fauna. However, none of the animal or plant area cladograms in Figure 7 show New Zealand and New Caledonia as sister areas. On the contrary, New Zealand is grouped with other austral areas (Australia and South America), whereas New Caledonia usually appears more closely related to New Guinea or the Southwest Pacific. Our results, therefore, do not support the notion of a common relict late Cretaceous Gondwanan biota in the two landmasses (Raven and Axelrod, 1972). Pole (1994) argued that the old Gondwanan lineages went extinct during the Paleogene submergence of New Zealand and New Caledonia and that the biotic similarity between these areas is explained by post-Eocene long-distance dispersal. Another possibility is that this dispersal took place via stepping stones across the Norfolk Ridge (Cracraft, 2001). According to Herzer et al. (1997), a new uplift of the Norfolk Ridge and the Reinga Ridge during the Miocene could have permitted at least intermittent land connections between New Zealand and New Caledonia until the late Tertiary. Molecular studies seem to support this explanation for New Zealand cicadas (Arensburger et al., 2004).
| Conclusions |
|---|
|
|
|---|
Unlike in the Northern Hemisphere, biogeography in the Southern Hemisphere is characterized by significant hierarchical distribution patterns. However, these patterns are different in the animals and plants we studied. In the animal groups, the distribution patterns largely conform to the geological breakup of Gondwana, particularly the SGP, and appear to have been generated by geologically induced vicariance. In the plants we studied, however, the hierarchical patterns, e.g., the PSP and the ISP, are incongruent with the commonly accepted sequence of geological events and have apparently been shaped by concordant dispersal and extinction events. In many of the plant groups, it appears that old, geologically induced vicariance patterns have been obscured by more recent dispersal. The resistant dispersal stage of plants (the seed) may facilitate colonization events, making plants less constrained by geological history than are many groups of animals.
When organism distributions are fitted onto a geologically based area cladogram, several types of concordant dispersal patterns emerge. The animal data set is dominated by trans-Antarctic dispersals, i.e., faunal links between South America and Australia that cannot be explained by vicariance. This pattern is presumably due to the long period of geological contact between Australia and South America via Antarctica. In plants, however, there has been a lot of exchange between landmasses that were not connected at the time, for instance, trans-Tasman dispersal between Australia and New Zealand. Both animal and plant data sets point to the hybrid nature of the South American biota, with surprisingly few dispersals between tropical and temperate regions of the continent. New Caledonia and New Zealand appear to have retained few old Gondwanan lineages, particularly those of plants.
Although our results clarify some points concerning Southern Hemisphere biogeography, many questions remain to be answered. Further studies should provide additional resolution by further development of the analytical tools, by expanding the data sets to include more groups and more completely resolved phylogenies, by incorporating the fossil distributions into the phylogenies (Manos, 1997; Swenson et al., 2001), and by dating the pylogenetic events in order to discriminate between different biogeographic scenarios (Sanmartín et al., 2001).
|
| Acknowledgment |
|---|
|
|
|---|
We thank Henrik Enghoff, Juan. J. Morrone, Ulf Swenson, and Lars Vilhelmsem for comments on the first drafts of this manuscript. Special thanks are extended to Henrik Enghoff for providing us with numerous animal cladograms used in this study. We are also grateful to Kåre Bremer, James Carpenter, Kjell-Arne Johanson, Peter Linder, Ewa Krzeminska, and Jyrki Muona for their help with the data sets and to Cliff Cuningham, Matt Lavin, Peter Linder, and Chris Simon for their constructive criticisms of the manuscript. We are especially indebted to Associate Editor Cliff Cuningham, who greatly helped to improve this manuscript with many useful suggestions. This research was supported through a European Community Marie Curie Fellowship (to I.S.) under the Improving Human Potential Program (Project MCFI-2000-00794) and by the Swedish Science Research Council (grants B10712-308 and 621-2001-2963 to F.R.).
| References |
|---|
|
|
|---|
-
Amorim D. S., Tozoni S. H. S. Phylogenetic and biogeographic analysis of the Anisopodoidea (Diptera: Bibionomorpha), with an area cladogram for intercontinental relationships. Rev. Bras. Entomol. (1994) 38:517–543.
Andersen N. M. Phylogeny and taxonomy of water striders, genus Aquarius Schellenberg (Insecta, Hemiptera, Gerridae), with a new species from Australia. Steenstrupia (1990) 16:37–81.
An update of the Angiosperm Phylogeny Group classification for the orders and families of flowering plants: APG II. Bot. J. Linn. Soc.—Angiosperm Phylogeny Group (APG), ed. (2003) 141:399–436.[CrossRef][Web of Science]
Arensburger P., Buckley T. R., Simon C., Moulds M., Holsinger K. E. Biogeography and phylogeny of the New Zealand cicada genera (Hemiptera: Cicadidae) based on nuclear and mitochondrial DNA data. J. Biogeogr. (2004) 31:1–13.
Barron E. J., Harrison C. G. A., Sloan J. L., Hay W. W. Paleogeography, 180 million years ago to the present. Eclogae Geol. Helv. (1981) 74:443–470.
Baum D. A., Small R. L., Wendel J. F. Biogeography and floral evolution of baobabs (Adansonia, Bombaceae) as inferred from multiple data sets. Syst. Biol. (1998) 47:181–207.
Bradford J. C. Molecular phylogenetics and morphological evolution in Cunonieae (Cunoniaceae). Ann. Mo. Bot. Gard. (2002) 89:491–503.[CrossRef]
Bradford J. C., Barnes R. W. Phylogenetics and classification of Cunoniaceae (Oxalidales) using chloroplast DNA sequences and morphology. Syst. Bot. (2001) 26:354–385.
Breure A. S. H. Systematics, phylogeny and zoogeography of Bulimulinae (Mollusca). Zool. Verh. Leiden (1979) 168:1–215.
Brothers R. N., Lillie A. R. Regional geology of New Caledonia. The ocean basin and margins. In: The Pacific Ocean, Volume 7B—Nairn A. E. M., Stehli F. G., Uyeda S., eds. (1988) New York: Plenum Press. 325–374. Pages.
Brundin L. Transantarctic relationships and their significance, as evidenced by chironomid midges with a monograph of the subfamilies Podonominae and Aphroteniinae and the austral Heptagynae. K. Sven. Vetenskapsakad. Handl. (1966) 11:1–472.
Burret C., Duhig N., Varne R. Asian and South-western Pacific continental terranes derived from Gondwana, and their biogeographic significance. Aust. Syst. Bot. (1991) 4:13–24.[CrossRef]
Buckley T. R., Arensburger P., Simon C., Chambers G. K. Combined data, Bayesian phylogenetics, and the origin of the New Zealand cicada genera. Syst. Biol. (2002) 51:4–18.
Carpenter J. Biogeographic patterns in the Vespidae (Hymenoptera): Two views of Africa and South America. In: Biological relationships between Africa and South America—Goldblatt P., ed. (1993) New Haven, Connecticut: Yale Univ. Press. 139–155. Pages.
Charleston M. A. Jungles: A new solution to the host/parasite phylogeny reconciliation problem. Math. Biosci. (1998) 149:191–223.[CrossRef][Web of Science][Medline]
Coode M. J. E. Aristotelia and Vallea, closely related in Elaeocarpaceae. Kew Bull. (1985) 40:479–507.[CrossRef]
Cooper A., Lalueza-Fox C., Anderson S., Rambaut A., Austin J., Ward R. Complete mitochondrial genome sequences of two extinct moas clarify ratite evolution. Nature (2001) 409:704–707.[CrossRef]
Cooper R. A., Milliner P. R. The New Zealand biota: Historical background and new research. Trends Ecol. Evol. (1993) 8:429–433.[CrossRef]
Cracraft J. Avian evolution, Gondwana biogeography and the Cretaceous–Tertiary mass extinction event. Proc. R. Soc. Lond. B (2001) 268:459–469.
Cranston P. S., Edward D. H. D. A systematic reappraisal of the Australian Aphroteniinae (Diptera: Chironomidae) with dating from vicariance biogeography. Syst. Entomol. (1992) 17:41–54.[CrossRef]
Cranston P. S., Edward D. H. D. Afrochlus Freeman: An African Gondwanan midge and the phylogeny of the Podonominae (Diptera: Chironomidae). Syst. Entomol. (1998) 23:77–90.[CrossRef]
Cranston P. S., Edward D. H. D. Botryocladius gen. n.: A new transantarctic genus of orthocladiine midge (Diptera: Chironomidae). Syst. Entomol. (1999) 24:305–333.[CrossRef]
Cranston P. S., Edward D. H. D., Cook L. G. New status, species, distribution records and phylogeny for Australian mandibulate Chironomidae (Diptera). Aust. J. Entomol. (2002) 41:357–366.[CrossRef]
Cranston P. S., Hare L. Conochironomus Freeman: An Afro-Australian Chironomini genus revised (Diptera: Chironomidae). Syst. Entomol. (1995) 20:247–264.[CrossRef]
Craw R. Phylogenetics, areas, geology and the biogeography of Croizat: A radical view. Syst. Zool. (1982) 31:304–316.
Crisci J. V., Cigliano M. M., Morrone J. J., Roig-Juñent S. Historical biogeography of southern South America. Syst. Zool. (1991) 40:152–171.
Crisci J. V., Gamundi I. J., Cabello M. N. A cladistic analysis of the genus Cyttaria (Fungi—Ascomycotina). Cladistics (1988) 4:279–290.[CrossRef][Web of Science]
Croizat L. Panbiogeography, Volumes 1, 2a, 2b. Privately published by the author. (1958) Caracas, Venezuela.
Croizat L., Nelson L. G., Rosen D. E. Centres of origin and related concepts. Syst. Zool. (1974) 23:265–287.
Damgaard J., Zettel H. Genetic diversity, species phylogeny and historical biogeography of the Aquarius paludum group (Heteroptera: Gerridae). Insect Syst. Evol. (2003) 34:313–328.
Edmunds G. F. Discussion. In: Vicariance biogeography—Nelson G., Rosen D. E., eds. (1981) New York: Columbia Univ. Press. 287–297. Pages.
Eizirik E., Murphy W. J., O'Brien S. J. Molecular dating and biogeography of the early placental mammal radiation. J. Hered. (2001) 92:212–219.
Farias I. P., Ortí G., Sampaio I., Schneider H., Meyer A. Mitochondrial DNA phylogeny of the family Cichlidae: Monophyly and fast molecular evolution of the Neotropical assemblage. J. Mol. Evol. (1999) 48:703–711.[CrossRef][Web of Science][Medline]
Fell H. B. West-Wind-Drift dispersal of echinoderms in the Southern Hemisphere. Nature (1962) 24:759–761.
Gauld I. D., Wahl D. B. The Labeninae (Hymenoptera: Ichneumonidae): A study in phylogenetic reconstruction and evolutionary biology. Zool. J. Linn. Soc. (2000) 129:271–347.[CrossRef]
Glasby C. J., Alvarez B. Distribution patterns and biogeographic analysis of austral Polychaeta (Annelida). J. Biogeogr. (1999) 26:507–533.[CrossRef]
Goldblatt P. Biological relationships between Africa and South America: An overview. In: Biological relationships between Africa and South America—Goldblatt P., ed. (1993) New Haven, Connecticut: Yale Univ. Press. 3–14. Pages.
Haddrath O., Baker A. J. Complete mitochondrial DNA genome sequences of extinct birds: Ratite phylogenetics and the vicariance biogeography hypothesis. Proc. R. Soc. Lond. B (2001) 268:939–945.
Hay W. W., DeConto R. M., Wold C. N., Wilson K. M., Voigt S., Schulz M., Wold-Rossby A., Dullo W. C., Ronov A. B., Balukhovsky A. N., Soeding E. Alternative global cretaceous paleogeography. In: Evolution of Cretaceous ocean-climate systems—Barrera E., Johnson C., eds. (1999) Boulder, Colorado: Geological Society of America. 1–47. Special Paper 332.
Henderson I. M. Biogeography without area? Aust. Syst. Bot. (1991) 4:59–71.[CrossRef]
Herzer R. H., Chaproniere G. C. H., Edwards A. R., Hollis C. J., Pelletier B., Raine J. I., Scott G. H., Stagpoole V., Strong C. P., Symonds P., Wilson G. J., Zhu H. Seismic stratigraphy and structural history of the Reinga Basin and its margins, southern Norfolk Ridge system. N.Z.J. Geol. Geophys. (1997) 40:425–451.
Hopper S. D., Fay M. F., Rosetto M., Chase M. W. A molecular phylogenetic analysis of the bloodroot and kangaroo paw family, Haemodoraceae: Taxonomic, biogeographic and conservation implications. Bot. J. Linn. Soc. (1999) 131:285–299.[CrossRef][Web of Science]
Hunt G. S., Cokendolpher J. C. Ballarinae, a new family of harvestmen from the Southern Hemisphere (Arachnida, Opiliones, Neopilinidae). Rec. Austr. Mus. (1991) 43:131–169.[CrossRef]
Johanson K. A. Phylogenetic and biogeographic analysis of the family Helicopsychidae (Insecta: Trichoptera). Entomol. Scan. Suppl. (1998) 53:1–172.
Kamp P. J. J. Pacifica and New Zealand, proposed eastern elements in Gondwanaland's history. Nature (1980) 288:659–664.[CrossRef]
Kårehed J., Lundberg J., Bremer B., Bremer K. Evolution of the Australasian families Alseuosmiaceae, Argophyllaceae, and Phellinaceae. Syst. Bot. (1999) 24:660–682.[CrossRef]
Karol K. G., Suh Y., Schaltz G. E., Zimmer E. A. Molecular evidence for the phylogenetic position of Takhtajania in the Winteraceae: Inference from nuclear ribosomal and chloroplast gene spacer sequences. Ann. M. Bot. Gard. (2000) 87:414–432.[CrossRef]
Krajewski C., Wroe S., Westerman M. Molecular evidence for the pattern and timing of cladogenesis in dasyurid marsupials. Zool. J. Linn. Soc. (2000) 130:375–404.[CrossRef]
Kress W. J., Schatz G. E., Andrianifahanana M., Morland H. S. Pollination of Ravenala madagascariensis (Strelitziaceae) by lemurs in Madagascar: Evidence for an archaic coevolution system? Am. J. Bot. (1994) 81:542–551.[CrossRef][Web of Science]
Lawver L. A., Gahagan L. M., Coffin M. F. The development of paleoseaways around Antarctica. The Antarctic paleoenvironment: A perspective on global change. Am. Geophys. Union Antarctic Res. Ser. (1992) 56:7–30.
Linder H. P. Rytidosperma vickeryae—A new danthonioid grass from Kosciuszko (New South Wales, Australia): Morphology, phylogeny and biogeography. Aust. Syst. Bot. (1999) 12:743–755.[CrossRef]
Linder H. P. Vicariance, climate change, anatomy and phylogeny of Restionaceae. Bot. J. Linn. Soc. (2000) 134:159–177.[Web of Science]
Linder H. P., Barker N. P. Biogeography of Danthonieae. In: Grasses: Systematics and evolution—Jacobs S. W. L., Everett J., eds. (2000) Melbourne: CSIRO. 231–238. Pages.
Linder H. P., Crisp M. D. Nothofagus and Pacific biogeography. Cladistics (1995) 11:5–32.[CrossRef][Web of Science]
Lopretto E. C., Morrone J. J. Anaspidacea, Bathynellacea (Crustacea, Syncarida), generalised tracks, and the biogeographic relationships of South America. Zool. Scr. (1998) 27:311–318.[CrossRef][Web of Science]
Macey J. R., Schulte J. A. II, Ananjeva N. B., Larson A., Wang Y., Pethiyagoda R., Rastegar-Pouyani N., Papenfuss T. J. Evaluating the trans-Tethys migration: An example using acrodont lizard phylogenetics. Syst. Biol. (2000) 49:233–256.
Manos P. S. Systematics of Nothofagus (Nothofagaceae) based on rDNA spacer sequences (ITS): Taxonomic congruence with morphology and plastid sequences. Am. J. Bot. (1997) 84:1137–1155.[Abstract]
Mathews E. G. Classification, phylogeny and biogeography of the genera of Adeliini (Coleoptera: Tenebrionidae). Invertebr. Taxon. (1998) 12:685–824.[CrossRef]
Matile L. Recherches sur la systemátique et l'évolution des Keroplatidae (Diptera, Mycetophiloidea). Mém. Mus. Natl. Hist. Nat. Paris Sér. A Zool. (1990) 148:1–682.
McLoughlin S. The breakup history of Gondwana and its impact on pre-Cenozoic floristic provincialism. Aust. J. Bot. (2001) 49:271–300.[CrossRef][Web of Science]
Metcalfe I. Pre-Cretaceous evolution of SE Asian terranes. In: Tectonic evolution of South-east Asia—Hall R., Blundell D. J., eds. (1996) Boulder, Colorado: Geological Society of America. 97–122. Special Publication 106.
Michaux B. Distributional patterns and tectonic development in Indonesia: Wallace reinterpreted. Aust. Syst. Bot. (1991) 4:25–36.[CrossRef]
Morrone J. J., Crisci J. V. Historical biogeography: Introduction to methods. Annu. Rev. Ecol. Syst. (1995) 26:373–401.[Web of Science]
Munroe D. D. The systematics, phylogeny, and zoogeography of Symmerus and Australosymmerus (Diptera: Micetophilidae, Ditomyiinae). Mém. Entomol. Soc. Can. (1974) 92:1–183.
Muona J. The Eucnemidae of South-east Asia and the Western-Pacific: A biogeographic study. Aust. J. Bot. (1991) 4:165–182.[CrossRef]
Murphy W. J., Collier G. E. A molecular phylogeny for apocheiloid fishes (Atherinomorpha, Cyprinodontiformes): The role of vicariance and the origins of annualism. Mol. Biol. Evol. (1997) 14:790–799.[Abstract]
Nelson G., Platnick N. I. Systematics and biogeography: Cladistics and vicariance (1981) New York: Columbia Univ. Press.
Ortí G., Meyer A. The radiation of characiform fishes and the limits of resolution of mitochondrial ribosomal DNA sequences. Syst. Biol. (1997) 46:75–100.
Page R. D. M. Parallel phylogenies: Reconstructing the history of host–parasite assemblages. Cladistics (1995) 10:155–173.[CrossRef][Web of Science]
Page R. D. M. Introduction. In: Cospeciation—Page R. D. M., ed. (2002) Chicago: Chicago Univ. Press. 1–22. Pages.
Page R. D. M., Charleston M. A. Trees within trees: Phylogeny and historical associations. Trends Ecol. Evol. (1998) 13:356–359.[CrossRef]
Palma R. E., Spotorno A. E. Molecular systematics of marsupials based on the rRNA 12S mitochondrial gene: The phylogeny of Dideldimorphia and of the living fossil microbiotheriid Dromiciops gliroides Thomas. Mol. Phylogenet. Evol. (1999) 13:525–535.[CrossRef][Web of Science][Medline]
Palmer K. R., Geissman J. Geologic time scale. (1999) Geological Society of America. Available at http://www.geosociety.org/science/timescale/timescl.pdf.
Patterson C. Methods of paleobiogeography. In: Vicariance biogeography: A critique—Nelson G., Rosen D. E., eds. (1981) New York: Columbia Univ. Press. 490–497. Pages.
Pigram C. J., Davies H. L. Terranes and the accretion history of the New Guinean orogen. Bur. Miner. Resour. J. Aust. Geol. Geophys. (1987) 10:193–211.
Pole M. The New Zealand flora–Entirely long-distance dispersal? J. Biogeogr. (1994) 21:625–635.[CrossRef]
Pole M. Can long-distance dispersal be inferred from the New Zealand plant fossil record? Aust. J. Bot. (2001) 49:357–366.
Pollock D. A. Classification, reconstructed phylogeny and geographical history of genera of Pilipalpinae (Coleoptera: Tenebrionoidea: Pyrochroidae). Invertebr. Taxon. (1995) 9:563–708.[CrossRef]
Qin T. K., Gullan P. J. A cladistic analysis of wax scales (Hemiptera: Coccoidea: Coccidae: Ceroplastinae). Syst. Entomol. (1995) 20:289–308.[CrossRef]
Qin T. K., Gullan P. J., Beattie A. C. Biogeography of the wax scales (Insecta: Hemiptera: Coccidae: Ceroplastinae). J. Biogeogr. (1998) 25:37–45.[CrossRef]
Rabinowitz P. D., Coffin M. F., Falvey D. The separation of Madagascar and Africa. Science (1983) 220:67–69.
Raven P. H., Axelrod D. I. Plate tectonics and Australasian paleobiogeography. Science (1972) 176:1379–1386.
Reeves G., Chase M. W., Goldblatt P., Rudall P., Fay M. F., Cox A. V., Lejeune B., Souza–Chies T. Molecular systematics of Iridaceae: Evidence from four plastid DNA regions. Am. J. Bot. (2001) 88:2074–2087.
Reid A. L. Review of the Peripaptosidae (Onicophora) in Australia, with comments on peripaptosid relationships. Invertebr. Taxon. (1996) 10:663–936.[CrossRef]
Reid A. L., Tait N. N., Briscoe D. A., Rowell D. M. Morphological, cytogenetic and allozymic variation within Cephalofovea (Onycophora, Peripaptosidae) with descriptions of three new species. Zool. J. Linn. Soc. (1995) 114:115–138.[CrossRef]
Renner S. S., Foreman D. B., Murray D. Timing transantarctic disjunctions in the Atherospermataceae (Laurales): Evidence from coding and noncoding chloroplast sequences. Syst. Biol. (2000) 49:579–591.
Robinson D. F., Foulds L. R. Comparison of phylogenetic trees. Math. Biosci. (1981) 53:131–147.[CrossRef][Web of Science]
Ronquist F. Reconstructing the history of host–parasite associations using generalized parsimony. Cladistics (1995) 11:73–89.[CrossRef][Web of Science]
Ronquist F. Dispersal–vicariance analysis: A new biogeographic approach to the quantification of historical biogeography. Syst. Biol. (1997) 46:195–203.
Ronquist F. Phylogenetic approaches in coevolution and biogeography. Zool. Scr. (1998) 26:313–322.[CrossRef][Web of Science]
Ronquist F. Parsimony analysis of coevolving species associations. In: Cospeciation—Page R. D. M., ed. (2002a) Chicago: Chicago Univ. Press. 22–64. Pages.
Ronquist F. TreeFitter, version 1.2. (2002b) Software available from http://morphbank.ebc.uu.se/TreeFitter.
Ronquist F., Nylin S. Process and pattern in the evolution of species associations. Syst. Zool. (1990) 39:323–344.
Sanmartín I., Enghoff H., Ronquist F. Patterns of animal dispersal, vicariance and diversification in the Holarctic. Biol. J. Linn. Soc. (2001) 73:345–390.[CrossRef][Web of Science]
Sanmartín I., Ronquist F. New solutions to old problems: Widespread taxa, redundant distributions and missing areas in event-based biogeography. Anim. Biodivers. Conserv. (2002) 25:75–93.
Schmidt F. Les Hydrobiosides (Trichoptera, Annulipalpia). Bull. Inst. R. Sci. Nat. Belg. Suppl. (1989) 59:1–154.
Scotese C. R., Gahagan L. M., Larson R. L. Plate tectonic reconstructions of the Cretaceous and Cenozoic ocean basins. Tectonophysics (1988) 155:27–48.[CrossRef][Web of Science]
Seberg O. Taxonomy, phylogeny and biogeography of the genus Oreobolus R. Br. (Cyperaceae), with comments on the biogeography of the South Pacific continents. Bot. J. Linn. Soc. (1988) 96:119–195.
Seberg O. Biogeographic congruence in the South Pacific. Aust. Syst. Bot. (1991) 4:127–136.[CrossRef]
Shaw S. R. Phylogeny and biogeography of the parasitoid wasp family Megalyridae (Hymenoptera). J. Biogeogr. (1990) 17:569–581.[CrossRef]
Shields O. Evidence for initial opening of the Pacific Ocean in the Jurassic. Palaeogeogr. Palaeoclimatol. Paleoecol. (1979) 26:181–220.[CrossRef]
Siddall M. E. Phylogenetic covariance probability: Confidence and historical associations. Syst. Biol. (1996) 45:48–66.
Springer M. S., Westerman M., Kavanagh J. R., Burk A., Woodburne M. O., Kao D. J., Krajewski C. The origin of the Australasian marsupial fauna and the phylogenetic affinities of the enigmatic monito del monte and marsupial mole. Proc. R. Soc. Lond. B. (1998) 265:2381–2386.
Storey M., Mahoney J. J., Saunders A. D., Duncan R. A., Kelley S. P., Coffin M. F. Timing of hot spot-related volcanism and the breakup of Madagascar and India. Science (1995) 267:852–855.
Swenson U., Bremer K. Pacific biogeography of the Asteraceae genus Abrotanella (Senecioneae, Blennospermatinae). Syst. Bot. (1997) 22:493–508.[CrossRef]
Swenson U., Hill R. S., McLoughlin S. Biogeography of Nothofagus supports the sequence of Gondwana break-up. Taxon (2001) 50:1–17.
Thayer M. K. Revision, phylogeny and biogeography of the austral genus Metacorneolabium Steel (Coleoptera: Staphylinidae: Omaliinae). In: Taxonomy, phylogeny and zoogeography of beetles and ants.—Ball G. E., ed. (1985) Dordrecht: Dr. W. Junk Publishers. 113–236. Pages.
Van Tuinen M., Sibley C. G., Hedges S. B. Phylogeny and biogeography of ratite birds inferred from DNA sequences of the mithocondrial ribosomal genes. Mol. Biol. Evol. (1998) 15:370–376.[Abstract]
Van Veller M. G. P., Kornet D. J., Zandee M. Two requirements for obtaining valid common patterns under assumptions 0, 1 and 2 in vicariance biogeography. Cladistics (1999) 15:393–406.[Web of Science]
Veevers J. J. Phanerozoic Australia in the changing configuration of Proto-Pangea through Gondwanaland and Pangea to the present dispersed continents. Aust. Syst. Bot. (1991) 4:1–11.
Veevers J. J., Powell C. M., Roots S. R. Review of sea floor spreading around Australia. I. Synthesis of the patterns of spreading. Aust. J. Earth Sci. (1991) 38:373–389.[CrossRef]
Walley A. M., Ross M. I. Preliminary reconstructions for the Cretaceous to Cainozoic of the New Zealand–New Caledonian Region (1991) Canberra: Bureau of Mineral Resources.
Wares J. P., Gaines S. D., Cunningham C. W. A comparative study of asymmetric migration events across a marine biogeographic boundary. Evolution (2001) 55:295–306.[Web of Science][Medline]
Waters J. M., Dijkstra L. H., Wallis G. P. Biogeography of a southern hemisphere freshwater fish: How important is marine dispersal? Mol. Ecol. (2000a) 9:1815–1821.[Medline]
Waters J. M., López J. A., Wallis G. P. Molecular phylogenetics and biogeography of galaxiid fishes (Osteichthyes: Galaxiidae): Dispersal, vicariance, and the position of Lepidogalaxias salamandroides. Syst. Biol. (2000b) 49:777–795.
Weston P. H., Crisp M. D. Cladistic biogeography of waratahs (Proteaceae: Embothriedae) and their allies across the Pacific. Aust. Syst. Bot. (1994) 7:225–249.[CrossRef]
Winkworth R. C., Wagstaff S. J., Glenny D., Lockhart P. Plant dispersal N.E.W.S. from New Zealand. Trends Ecol. Evol. (2002) 17:514–520.[CrossRef]
Woodburne M. O., Case J. A. Dispersal, vicariance, and the late Cretaceous to early Tertiary land mammal biogeography from South America to Australia. J. Mammal. Evol. (1996) 3:121–161.[CrossRef]
Yeates D. K., Irwin M. E. Apioceridae (Insecta: Diptera): Cladistic reappraisal and biogeography. Zool. J. Linn. Soc. (1996) 116:247–301.[CrossRef]
This article has been cited by other articles:
![]() |
M. J. Phillips, G. C. Gibb, E. A. Crimp, and D. Penny Tinamous and Moa Flock Together: Mitochondrial Genome Sequence Analysis Reveals Independent Losses of Flight among Ratites Syst Biol, November 13, 2009; (2009) syp079v1. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Balke, I. Ribera, L. Hendrich, M. A. Miller, K. Sagata, A. Posman, A. P. Vogler, and R. Meier New Guinea highland origin of a widespread arthropod supertramp Proc R Soc B, July 7, 2009; 276(1666): 2359 - 2367. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. E.H Jones, A. J.D Tennyson, J. P Worthy, S. E Evans, and T. H Worthy A sphenodontine (Rhynchocephalia) from the Miocene of New Zealand and palaeobiogeography of the tuatara (Sphenodon) Proc R Soc B, April 7, 2009; 276(1660): 1385 - 1390. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Goldberg, S. A Trewick, and A. M Paterson Evolution of New Zealand's terrestrial fauna: a review of molecular evidence Phil Trans R Soc B, October 27, 2008; 363(1508): 3319 - 3334. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. D Smith, P. J Makovicky, F. L Agnolin, M. D Ezcurra, D. F Pais, and S. W Salisbury A Megaraptor-like theropod (Dinosauria: Tetanurae) in Australia: support for faunal exchange across eastern and western Gondwana in the Mid-Cretaceous Proc R Soc B, September 22, 2008; 275(1647): 2085 - 2093. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. R. Mast, C. L. Willis, E. H. Jones, K. M. Downs, and P. H. Weston A smaller Macadamia from a more vagile tribe: inference of phylogenetic relationships, divergence times, and diaspore evolution in Macadamia and relatives (tribe Macadamieae; Proteaceae) Am. J. Botany, July 1, 2008; 95(7): 843 - 870. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Huttunen, L. Hedenas, M. S. Ignatov, N. Devos, and A. Vanderpoorten Origin and evolution of the northern hemisphere disjunction in the moss genus Homalothecium (Brachytheciaceae) Am. J. Botany, June 1, 2008; 95(6): 720 - 730. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. A. A. Nylander, U. Olsson, P. Alstrom, and I. Sanmartin Accounting for Phylogenetic Uncertainty in Biogeography: A Bayesian Approach to Dispersal-Vicariance Analysis of the Thrushes (Aves: Turdus) Syst Biol, April 1, 2008; 57(2): 257 - 268. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. D. Crisp Ghosts of Gondwana: The History of Life in New Zealand.--George Gibbs. 2006, reprinted 2007. Craig Potton Publishing, Nelson, New Zealand. 232 pp. ISBN 978-1-877333-48-4 (ISBN-10 1-877333-48-4). NZ$49.99 (hardcover) Syst Biol, April 1, 2008; 57(2): 329 - 332. [Full Text] [PDF] |
||||
![]() |
A. D. Rogers Evolution and biodiversity of Antarctic organisms: a molecular perspective Phil Trans R Soc B, December 29, 2007; 362(1488): 2191 - 2214. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. M. Waters Malte C. Ebach and Raymond S. Tangney (editors), Biogeography in a Changing World Syst Biol, October 1, 2007; 56(5): 871 - 873. [Full Text] [PDF] |
||||
![]() |
J. M. Waters, D. L. Rowe, S. Apte, T. M. King, G. P. Wallis, L. Anderson, R. J. Norris, D. Craw, and C. P. Burridge Geological Dates and Molecular Rates: Rapid Divergence of Rivers and Their Biotas Syst Biol, April 1, 2007; 56(2): 271 - 282. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. Wahlberg That Awkward Age for Butterflies: Insights from the Age of the Butterfly Subfamily Nymphalinae (Lepidoptera: Nymphalidae) Syst Biol, October 1, 2006; 55(5): 703 - 714. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. Bossuyt, R. M. Brown, D. M. Hillis, D. C. Cannatella, and M. C. Milinkovitch Phylogeny and Biogeography of a Cosmopolitan Frog Radiation: Late Cretaceous Diversification Resulted in Continent-Scale Endemism in the Family Ranidae Syst Biol, August 1, 2006; 55(4): 579 - 594. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Poux, P. Chevret, D. Huchon, W. W. de Jong, and E. J. P. Douzery Arrival and Diversification of Caviomorph Rodents and Platyrrhine Primates in South America Syst Biol, April 1, 2006; 55(2): 228 - 244. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. M. Waters and D. Craw Goodbye Gondwana? New Zealand Biogeography, Geology, and the Problem of Circularity Syst Biol, April 1, 2006; 55(2): 351 - 356. [Full Text] [PDF] |
||||
![]() |
M. P. Schwarz, S. Fuller, S. M. Tierney, and S. J. B. Cooper Molecular Phylogenetics of the Exoneurine Allodapine Bees Reveal an Ancient and Puzzling Dispersal from Africa to Australia Syst Biol, February 1, 2006; 55(1): 31 - 45. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. G Cook and M. D Crisp Not so ancient: the extant crown group of Nothofagus represents a post-Gondwanan radiation Proc R Soc B, December 7, 2005; 272(1580): 2535 - 2544. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. T Monaghan, J.-L. Gattolliat, M. Sartori, J.-M. Elouard, H. James, P. Derleth, O. Glaizot, F. de Moor, and A. P Vogler Trans-oceanic and endemic origins of the small minnow mayflies (Ephemeroptera, Baetidae) of Madagascar Proc R Soc B, September 7, 2005; 272(1574): 1829 - 1836. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. M. Waters Historical Bioegeography: An Introduction.--J. V. Crisci, L. Katinas, and P. Posadas. 2003. Harvard University Press, Cambridge MA. 264 pp. ISBN 0-674-01059-0. $45.00 hardcover. Syst Biol, April 1, 2005; 54(2): 338 - 340. [Full Text] [PDF] |
||||
![]() |
Y.-M. Yuan, S. Wohlhauser, M. Moller, J. Klackenberg, M. W. Callmander, and P. Kupfer Phylogeny and Biogeography of Exacum (Gentianaceae): A Disjunctive Distribution in the Indian Ocean Basin Resulted from Long Distance Dispersal and Extensive Radiation Syst Biol, February 1, 2005; 54(1): 21 - 34. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||














