|
|
||||||||
|
What's this? |
Systemics and Phytogeography |
2Missouri Botanical Garden, Graduate Studies Program, P.O. Box 299, St. Louis, Missouri 63166 USA; 3Department of Biology, Washington University, One Brookings Drive, Campus Box 1229, St. Louis, Missouri 63130-4899 USA; 4Research Center for Biodiversity, Academia Sinica, Nankang, Taipei 115, Taiwan; 5Department of Plant Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801 USA; 6Department of Plant Systematics and Geography, Warsaw University, Aleje Ujazdowskie 4, PL-00-478 Warszawa, Poland
Received for publication March 18, 2005. Accepted for publication August 22, 2005.
ABSTRACT
The alpine ecosystem is the only terrestrial biogeographic unit that is distributed globally. Studying phylogenetics of the plant species in this widespread ecosystem can provide insights into the historical biogeographic processes that have shaped the global biodiversity. The trans-Pacific disjunct alpine genus Oreomyrrhis (Apiaceae) was investigated using nrDNA ITS sequences to test the taxonomic and biogeographic hypotheses. Phylogenetic analyses using maximum parsimony, maximum likelihood, and Bayesian inference revealed that species of Oreomyrrhis form a weakly supported monophyletic clade that is nested within Chaerophyllum sect. Chaerophyllum (subtribe Scandicinae, tribe Scandiceae). The optimal solutions of dispersal-vicariance analysis indicate that the ancestor of Chaerophyllum sect. Chaerophyllum (including Oreomyrrhis) was distributed in Eurasia and subsequently dispersed to North America and southern Pacific Rim. Based on dating using ITS sequence variation, these dispersal events were most likely recent, probably during late Tertiary to Quaternary. The structure of the ITS haplotype network suggests that a rapid range expansion via long-distance dispersal had been crucial in generating the trans-Pacific disjunction of Oreomyrrhis. Furthermore, evolution toward smaller mericarp size and a transition from outcrossing to selfing during Oreomyrrhis's evolution might have increased the chances for long-distance dispersal, facilitating its range expansion and occupation on alpine environments.
Key Words: alpine plants Apiaceae Chaerophyllum internal transcribed spacer (ITS) long-distance dispersal Oreomyrrhis Pleistocene glaciations South Pacific disjunction
The alpine ecosystem is the only terrestrial biogeographic unit that is distributed globally (Körner, 2003
). Studying phylogenetics of the plant species in this widespread ecosystem presents a unique opportunity to understand the historical biogeographic processes that have shaped global biodiversity. Recent studies in Europe suggest that dispersal dynamics, Pleistocene glaciations, and landscape heterogeneity have been the most influential factors shaping the geographic distribution and genetic variation in alpine plants (Tribsch and Stuessy, 2003
). However, few studies have attempted to understand how these historical, geological, and biological factors have affected the alpine biodiversity outside of Europe, on a global scale. The trans-Pacific disjunct genus, Oreomyrrhis, found in isolated alpine life zones across different southern latitudes around the southern Pacific Basin, presents an ideal system to test these biogeographic hypotheses.
Oreomyrrhis Endl. (Apiaceae subfamily Apioideae), as circumscribed in the seminal work of Mathias and Constance (1955)
and with novelties described subsequently by Mathias and Constance (1977)
and Chen and Wang (2001)
, consists of ca. 25 species distributed in alpine and subalpine areas and sub-Antarctic islands around the southern Pacific Rim (Mexico, Central American highlands, northern South American Andes, Tierra del Fuego, Falkland Islands, New Zealand, SE Australia, Tasmania, New Guinea, Borneo, and Taiwan; Fig. 1). It is one of the few genera in Apioideae with a mainly Southern Hemisphere distribution (Mathias, 1965
). Known as an "exemplary demonstration of the sub-Antarctic-South Pacific pattern of distribution" (van Steenis, 1963
), its unique geographic pattern has intrigued generations of phytogeographers (Hayata, 1911
; Merrill, 1918
; Mathias and Constance, 1955
; van Steenis, 1962
, 1963
, 1964
; Raven, 1973
; Melville, 1981
).
|
The taxonomic history of Oreomyrrhis began with the description of Myrrhis andicola by Kunth in 1821, based on a Humboldt and Bonpland collection from Ecuador (Humboldt et al., 1821
). Later that year, a new genus/species Caldasia chaerophylloides Lag. was described from Peru (La Gasca, 1821
). De Candolle (1829
, 1830
) accepted this new genus and transferred M. andicola to Caldasia. However, because Caldasia Lag. was a later homonym, Endlicher (1839)
proposed to replace it with the new generic name Oreomyrrhis, which has been accepted since that time (Mathias and Constance, 1955
).
Morphologically, Oreomyrrhis is characterized by an inflorescence consisting of a simple umbel borne terminally on peduncles that arise from the bases of sheathing and often rosetted leaves (Mathias and Constance, 1955
). Such a simple umbel is extremely rare in Apioideae, appearing consistently only in Oreomyrrhis, Lilaeopsis Greene, and Neogoezia Hemsl. (Constance, 1987
). The first tribal classification that included Oreomyrrhis placed the genus (as Caldasia) in tribe Scandiceae Spreng. (de Candolle, 1830
), which was challenged by Hooker (1856)
but followed by Bentham (1867)
and Koso-Poljansky (1916)
. Emphasizing the presence or absence of vittae (secretory canals) and secondary ridges on the mericarps, Bentham (1866)
suggested that, among the genera of Scandiceae, Oreomyrrhis is most similar to Chaerophyllum L. of the subtribe Scandicinae Tausch.
In one of the most comprehensive and influential taxonomic works on Apiaceae, Drude (1898)
, whose taxonomy emphasized microscopic structure of the mericarps, removed Oreomyrrhis from Scandiceae on the basis of mericarp shape and the absence of crystal druses around the carpophore (Mathias and Constance, 1955
). Drude placed Oreomyrrhis in tribe Smyrnieae Spreng., that included other simple or irregularly compound-umbelled genera, such as Neogoezia, Apiastrum Nutt. ex Torr. et A. Gray, Erigenia Nutt., and Orogenia S. Watson (Mathias and Constance, 1955
). Drude's treatment has been followed by many subsequent authors (e.g., Shan and Sheh, 1979
) and was adopted with minor changes in the most recent synopsis of Apiaceae (Pimenov and Leonov, 1993
). However, the phylogenetic position of Oreomyrrhis was left unanswered by Mathias and Constance (1955
, p. 349) as they concluded "... We are nevertheless inclined to regard Oreomyrrhis as having a vague relationship to Chaerophyllum and probably no close affinity with other genera bearing simple umbels. Although Oreomyrrhis is clearly referable to the subfamily Apioideae, we have been unable to recognize any close relative among the living Umbelliferae of either the northern or southern hemisphere."
All 13 species of Oreomyrrhis investigated to date have a base chromosome number of n = 6 (Mathias and Constance, 1955
; Moore, 1971
; Chen and Wang, 2001
; Pimenov et al., 2003
). This number is rare in Apioideae, representing only 4% of the 1461 species with published chromosome counts in the subfamily (Moore, 1971
; Pimenov et al., 2003
). The rarity of its base number has stimulated much speculation regarding Oreomyrrhis's phylogenetic position (Moore, 1971
; Raven, 1975
; Pimenov et al., 2003
). Recently, serological (Shneyer et al., 1992
) and molecular phylogenetic studies (Downie and Katz-Downie, 1996
; Downie et al., 2001
) have shown that Drude's tribe Smyrnieae is a heterogeneous and polyphyletic assemblage. Unfortunately, Oreomyrrhis has not been included in any modern systematic study and its phylogenetic affinity remains puzzling.
In addition to its dubious tribal position, the infrageneric taxonomy of Oreomyrrhis has been equally contentious (Mathias and Constance, 1955
; Stevens, 1990
). The majority of these controversies have centered on the delimitation of O. andicola (Kunth) Endl. ex Hook. f. (Mathias and Constance, 1955
). First described from the Ecuadorian Andes (as Myrrhis andicola), subsequent collections extending from the northern Andes north to the Mexican highlands indicate that O. andicola is an extremely polymorphic complex. These populations display continuous morphological variation in habit, size, leaf shape, color, and pubescence that nearly spans the range of morphological diversity of the entire genus (Mathias and Constance, 1955
; Melville, 1981
). While there have been attempts to divide the Andean Oreomyrrhis into multiple species (e.g., Urban, 18801882
; Johnston, 1938
), some authors have extended the species boundary of O. andicola beyond South America to include populations in New Zealand (Kirk, 1899
), Australia (Bentham, 1866
), and New Guinea (Buwalda, 1951
).
By restricting the name Oreomyrrhis andicola to encompass only the Andean populations and by accepting most regional endemics, Mathias and Constance (1955)
recognized 23 species of Oreomyrrhis (Fig. 1). Similar species recognition criteria were later adopted in Mathias and Constance (1977)
and Chen and Wang (2001)
to describe new species from New Guinea and Taiwan, respectively. However, the species delimitation by Mathias and Constance (1955)
was later criticized by Kern (1958)
as being too narrow. Kern (1958)
surmised that Oreomyrrhis comprises only a few species, with O. andicola being a species with a circum-Pacific distribution from Mexico to Taiwan. Despite Kern's criticism, Mathias and Constance's delimitations have been widely adopted in most local floras (e.g., Allan, 1961
; van Royen, 1983
; Walsh and Entwisle, 1999
). Based on the work of Mathias and Constance (1955)
, Melville (1981)
reconstructed the infrageneric relationships of Oreomyrrhis and, based on his own hypothesis of Oreomyrrhis's phylogeny, proposed a vicariance scenario that is linked with ancient geological events. Melville's phylogenetic hypothesis was later adopted by Crisci et al. (1991)
for a historical biogeographic analysis of southern South America.
In recent years, DNA sequence data have proven to be a promising and powerful tool for elucidating relationships within Apiaceae that have been difficult to resolve with morphological data alone (Downie et al., 2001
). In particular, the large quantity of data in GenBank from the nuclear ribosomal DNA (nrDNA) internal transcribed spacer (ITS) region for numerous taxa of Apiaceae are a valuable asset for testing taxonomic hypotheses surrounding Oreomyrrhis. Moreover, ITS sequence data are useful for reconstructing infrageneric relationships (e.g., Downie et al., 2004
; Neves and Watson, 2004
) and for inferring the biogeographic histories of Apiaceae genera (e.g., Wen et al., 2002
).
In an ongoing project of global alpine biogeography, Oreomyrrhis was investigated using a combination of phylogenetic, phylogeographic, and population genetic approaches. In this study, ITS sequence data were analyzed by various phylogenetic methods to address the following questions: (1) Is Oreomyrrhis monophyletic? (2) What are its closest relatives? (3) Does vicariance or long-distance dispersal better explain the present distribution of Oreomyrrhis?
MATERIALS AND METHODS
Sampling strategy for Oreomyrrhis
Species of Oreomyrrhis were sampled to optimize coverage of the geographic distribution of the genus (Fig. 1) and intraspecific morphological variation. Leaves were collected and preserved on silica gel by the first author during fieldwork in Ecuador (O. andicola), Taiwan (O. involucrata var. pubescens, O. nanhuensis, and O. taiwaniana), SE Australia (O. brevipes, O. ciliata, O. eriopoda, and O. pulvinifica), Tasmania (O. argentea, O. ciliata, O. eriopoda, and O. sessiliflora), New Zealand (O. colensoi, O. ramosa, and O. rigida), and Guatemala (O. daucifolia) and through the courtesy of local researchers in Taiwan (O. involucrata and O. taiwaniana), Tierra del Fuego (O. hookeri), the Falkland Islands (O. hookeri), and Mexico (O. orizabae). DNA was also extracted from herbarium specimens for some taxa (O. azorellacea, O. borneensis, O. gunnii, O. involucrata, O. linearis, O. papuana, O. pumila, and O. sessiliflora), including two of the most unique New Guinea species, O. azorellacea (moss-like cushion plant) and O. linearis (characterized by grass-like linear leaves). The species not sampled in this study include the rare New Guinean species O. buwaldiana and O. plicata, that are known only by type collections, and the Mexican high alpine species O. tolucana for which no ITS sequence was amplified after several trials using herbarium material.
To evaluate the monophyly of species with broad distributions and substantial morphological variation, multiple individuals (number in parentheses) from different localities were sampled for O. andicola (5), O. argentea (2), O. ciliata (5), O. colensoi (4), O. eriopoda (6), O. hookeri (2), O. involucrata (4), O. linearis (2), O. ramosa (4), O. rigida (4), O. sessiliflora (2), and O. taiwaniana (4). Oreomyrrhis involucrata var. pubescens, a taxon not recognized by Mathias and Constance (1955)
and Chen and Wang (2001)
but readily recognizable in the field (K. Chung, personal observation), was also surveyed. One individual from Awahokomo karstland, Otago, New Zealand, previously identified as O. rigida (Molloy et al., 1999
), but likely a new species (P. Heenan, Landcare Research, personal communication), was also investigated. In total, 57 accessions of Oreomyrrhis, representing 22 of the 25 species recognized and/or later described by Mathias and Constance (1955
, 1977
) and Chen and Wang (2001)
, were included. The taxa sampled, taxon authorities, voucher information, and Genbank accession numbers are listed in the Appendix. Throughout this study, Mathias and Constance's (1955)
species circumscription was closely followed for plant identification, with further verification based on examination of type specimens and/or specimens annotated by Mathias and Constance (1955)
.
DNA extraction, PCR amplification, and sequencing
Total genomic DNA was isolated from silica-dried leaf material using the CTAB procedure of Doyle and Doyle (1987)
. For herbarium specimens, DNA was isolated using a VIOGENE Plant Genomic DNA extraction miniprep system (Viogen U.S.A., Sunnyvale, California) following the manufacturer's protocol. The entire ITS region (ITS-1, 5.8S, and ITS-2) was amplified using primers ITS1 and ITS4 described in White et al. (1990)
. For older herbarium material, the internal primers ITS2 and ITS3 (Downie and Katz-Downie, 1996
) or primers ITS-C and ITS-D (Blattner, 1999
), in combination with the external primers ITS1 and ITS4, were used to amplify the entire ITS region. All amplifications were performed in a 50-µL reaction. Polymerase chain reaction (PCR) conditions were set for an initial 4 min denaturation at 94°C, followed by 35 cycles of 94°C for 30 s, 50°C for 30 s, and 72°C for 1 min, with a final extension at 72°C for 10 min and a final incubation at 30°C for 10 min. PCR products (templates) were visualized via agarose gel electrophoresis and were further purified using the QIAquick Gel Extraction Kit (Qiagen, Valencia, California). Sequence reactions were carried out in both directions for each purified PCR template using the ABI Prism dGTP BigDye Terminator Ready Reaction kit (Applied Biosystems, Foster City, California, USA) following the manufacturer's protocol. Excessive Bigdye in the sequence reactions was removed using CENTRI-SEP columns (Princeton Separations, Adelphia, New Jersey, USA). Sequence reaction products were visualized using either an Applied Biosystem 373 automated DNA sequencer or a BaseStation DNA Fragment analyzer (MJ Research, South San Francisco, California, USA). Contigs were assembled using SEQMAN (Swindell and Plasterer, 1997
). All newly acquired sequences have been archived in Genbank (Appendix).
Molecular cloning
To examine the extent of sequence homogenization among reiterated ITS copies, molecular cloning was conducted on 10 accessions of nine Oreomyrrhis species (O. andicola [Chung 1493], O. daucifolia [Chung 1650], O. eriopoda [Chung 1554 and Chung 1612], O. hookeri [D. Broughton, no voucher], O. involucrata [Juan 53], O. nanhuensis [Chung 1510], O. orizabae [Avendaño R. 5349], O. ramosa [Chung 1630], and O. rigida [Chung 1618]). The purified ITS templates were ligated to the pGEM-T Easy vector (Promega, Madison, Wisconsin, USA) following the manufacturer's protocol and subsequently transformed into competent cells. After an overnight culture at 37°C on the LB ampicillin/IPTG/X-gal selective plate, colonies carrying the ITS insert were identified by color (white) and further verified by PCR using the T7 and SP6 promoter primer pairs (Promega, Madison, Wisconsin, USA). For each accession, up to 20 white colonies were checked for inserts. The amplified PCR products were purified and subsequently sequenced using the ITS1 and ITS4 primers as described in the previous section.
Tests of sequence evolution
The rate constancy of the ITS sequence evolution of Oreomyrrhis and its closest relatives Chaerophyllum procumbens and C. tainturieri (see Figs. 3 and 4) were inferred from relative rate tests (Wu and Li, 1985
) based on Kimura two-parameter distance (K2P; Kimura, 1980
), implemented by the program K2WuLi (Jermiin, 1996
). Chaerophyllum temulum was used as a reference taxon. Given a constant ITS evolution rate, approximate times of diversification with 95% confidence intervals between Oreomyrrhis and its outgroup and among Oreomyrrhis species were estimated as ± 2 SE of the K2P distances, calibrated with an ITS sequence divergence rate of 0.79% per million years (Sang et al., 1994
; Kropf et al., 2002
). This more conservative ITS divergence rate was chosen because we were mainly interested in testing competing biogeographic scenarios related to recent (late Tertiary-Quaternary) vs. old (Jurassic or early Tertiary) diversification of Oreomyrrhis, as predicted by long-distance dispersal (e.g., Dawson, 1971
; Raven, 1973
) vs. vicariance hypotheses (Melville, 1981
; van Steenis, 1964
), respectively. Pairwise K2P distances and their SEs were calculated using MEGA2 (Kumar et al., 2001
).
|
|
Analysis 1
To elucidate Oreomyrrhis's phylogenetic position within Apioideae, the data matrix of Downie and Katz-Downie (1996)
was adopted with some modifications. This matrix was composed of ITS-1 and ITS-2 sequences of representatives of tribe Smyrnieae, subtribes Scandicinae, Daucinae, and Torilidinae of tribe Scandiceae, and several major clades of Apioideae identified by recent molecular systematic studies (Downie et al., 2001
). Pleurospermum foetens Franch. (tribe Symrnieae, sensu Drude [1898]
) was chosen as outgroup (cf. Downie et al., 2001
). ITS sequences of eight Oreomyrrhis species from different geographic regions [O. orizabae (Mexico), O. andicola (South American Andes), O. hookeri (Tierra del Fuego), O. colensoi (New Zealand), O. eriopoda (Australia), O. papuana (New Guinea), O. borneensis (Borneo), and O. taiwaniana (Taiwan)] were added to the matrix. Chaerophyllum temulum L., the type species of the largest genus of Scandiceae (Spalik and Downie, 2001
), was added to test Bentham's (1866)
proposition that Oreomyrrhis is most similar to Chaerophyllum. Other additions included selected taxa of Drude's (1898)
tribe Smyrnieae that possess irregularly compound umbels (Erigenia, Neogoezia, and Orogenia), and Lilaeopsis, another predominantly Southern Hemispheric distributed genus also possessing a simple umbel (Constance, 1987
; Petersen et al., 2001
). Sequences were aligned using DAMBE (Xia and Xie, 2001
) and were manually adjusted in the alignment editor Se-Al (Rambaut, 1996
). Only those positions that were aligned unambiguously were used in the phylogenetic analyses; the aligned data matrix is available in TreeBASE (study accession number S1350, matrix accession number M2384).
Phylogenetic analyses were conducted using maximum parsimony (MP) algorithm implemented in PAUP* version 4.0b10 (Swofford, 2003
). All characters were treated as unordered, and all character transformations were weighted equally. A heuristic search with 1000 random addition replicates was performed with the following options in effect: MULTREES, TBR branch swapping, and gaps treated as missing data. Clade support was assessed by bootstrap (Felsenstein, 1985
) and decay (Bremer, 1988
) analyses. One hundred bootstrap replicates were analyzed by PAUP* using the heuristic search option with 10 random sequence additions. Decay indices were calculated using the program TreeRot.v2 (Sorenson, 1999
). To evaluate the phylogenetic content of the data, the g1 statistic of tree-length distribution skewness was calculated from 10 000 random MP trees generated by PAUP*, and was subsequently compared to the critical values provided in Hillis and Huelsenbeck (1992)
.
Analysis 2
The results of analysis 1 unambiguously placed Oreomyrrhis in subtribe Scandicinae of tribe Scandiceae, with Chaerophyllum being its closest relative (Fig. 2). To further elucidate their relationships within Scandicinae, ITS sequences of 57 Oreomyrrhis accessions (22 species) were combined with all available ITS sequences (ITS-1 and ITS-2) of the genus Chaerophyllum (22 species) and representatives of other major clades in the subtribe (Anthriscus, Kozlovia, Myrrhis, Osmorhiza, Scandix, Sphallerocarpus; see Spalik and Downie, 2001
). Athamanta cretensis and Conopodium ramosum were selected as outgroups based on recent phylogenetic analysis (Downie et al., 2000
). Voucher information of sampled taxa is shown in the Appendix.
|
The gap-excluding data matrix was then analyzed using maximum likelihood (ML). Modeltest 3.5 (Posada and Crandall, 1998
) was used to identify the appropriate evolutionary model of nucleotide substitution. Based on the model selected by Modeltest 3.5, 10 heuristic searches were performed in PAUP* using random addition sequence and TBR branch swapping. To assess clade support, a bootstrap analysis of 1000 replicates was conducted using the neighbor-joining method with ML distance (e.g., Downie et al., 2004
; Sun et al., 2004
).
A Bayesian phylogenetic analysis was also performed using the program MrBayes 3 (Ronquist and Huelsenbeck, 2003
) based on the same nucleotide substitution model used in ML analysis. In MrBayes, parameters were set as follows: random starting tree, 2 000 000 generations runs with sampling occurring every 100 generations, and a 4000-tree "burn-in" period. PAUP* was used to generate a 50% majority consensus tree of the output tree file from MrBayes.
Analysis 3
Results of analysis 2 revealed that Oreomyrrhis accessions formed a poorly resolved polytomous clade probably resulting from low sequence divergences. Such levels of sequence variation suggest ongoing speciation processes in Oreomyrrhis. To better elucidate relationships among Oreomyrrhis accessions, a haplotype network was constructed. Theoretically, haplotype networks can straddle populationspecies interfaces where speciation occurs, providing a powerful tool to investigate evolutionary processes (Templeton et al., 2000
; Schaal and Leverich, 2001
). We compiled an ITS sequence (ITS-1, 5.8S, and ITS-2) matrix of 57 Oreomyrrhis accessions. Based on the results of analysis 2 (see Figs. 3 and 4), Chaerophyllum procumbens and C. tainturieri were selected as outgroups. The haplotype network was calculated using statistical parsimony (Templeton et al., 1992
) in the program TCS (Clement et al., 2000
), with alignment gaps treated as a fifth state.
Biogeographic analysis
A dispersal-vicariance analysis (DIVA; Ronquist, 1997
) was performed to infer biogeographic history of Oreomyrrhis and subtribe Scandicinae. Given a phylogeny, DIVA reconstructs distribution of ancestral areas that minimize dispersal and extinction events under a parsimony criterion (Ronquist, 1997
), presenting an ideal tool to identify sources and directions of dispersal for Oreomyrrhis and its close relatives. For implementing DIVA, a fully resolved phylogeny is demanded. Because parts of the Scandicinae phylogeny are unresolved, clades composed of taxa distributed in the same general unit area were compressed to generate a simplified bifurcating phylogeny that is required by DIVA. Information of taxon distribution of subtribe Scandicinae was gathered from Spalik and Downie (2001
; Fig. 3). Athamanta and Conopodium were chosen as outgroups (Spalik et al., 2001a
). Since we were mainly interested in identifying the ancestral source of dispersal for Oreomyrrhis, not the exact locality of initial diversification, the unresolved clade of Oreomyrrhis (Fig. 4A) was also compressed to a single terminal and designated with a Pacific distribution (D). Three additional general unit areas were defined to cover the distributional range of Scandicinae: Eurasia (including North Africa; A), North America (B), and South America (C).
RESULTS
Characteristics and evolution of Oreomyrrhis ITS sequences
From 57 accessions of 22 Oreomyrrhis species, 26 distinct ITS sequences (haplotypes) were obtained. The alignment of the polymorphic sites, haplotype designations, and their frequencies are shown in Table 1. The length of the entire ITS region for the 26 Oreomyrrhis haplotypes ranged from 602 to 607 bp (ITS-1: 217218 bp; 5.8S: 163164 bp; ITS-2: 221226 bp). Alignment of the 26 haplotypes resulted in a matrix of 610 positions (ITS-1: 219 bp; 5.8S: 164 bp; ITS-2: 227 bp), including two 1-bp gaps in ITS-1, a single 1-bp gap in 5.8S, and three 1-bp gaps and a 4-bp gap in ITS-2.
|
Within the 378 relative rate tests among C. procumbens, C. tainturieri, and the 26 Oreomyrrhis ITS haplotypes, there were only two tests where rate constancy of ITS sequence evolution was rejected significantly (P < 0.05; z score > |1.96|). This figure is much lower than the expected number by chance alone (one of every 20 comparisons will be rejected significantly at the 5% level by chance; Sytsma and Schaal, 1985
; Abbott and Comes, 2004
). Given the apparent rate constancy of ITS evolution, times of speciation and diversification were approximated based on pairwise sequence divergences, calibrated with an ITS sequence divergence rate of 0.79% per million years. The average pairwise sequence divergence between C. temulum and Oreomyrrhis and its closest relatives (C. procumbens and C. tainturieri) was 6.45 ± 1.18%, corresponding to a separation time of 8.16 (95% confidence intervals = 5.1811.15) million years ago (mya). Between Oreomyrrhis and its closest relatives (C. procumbens and C. tainturieri), the average pairwise sequence divergence is 1.76 ± 0.49%, indicating that Oreomyrrhis diverged from its most recent common ancestor at 2.23 (95% confidence intervals = 0.993.47) mya. Within Oreomyrrhis, the average pairwise sequence divergence is 0.80 ± 0.18%, corresponding to a diversification time of 1.01 (95% confidence intervals = 0.56 1.47) mya.
Analysis 1
Alignment of the ITS-1 and ITS-2 sequences of the seven Oreomyrrhis ITS sequences (O. hookeri and O. orizabae are identical) of eight Oreomyrrhis species and 44 representatives of Apioideae resulted in a data matrix of 475 positions. Due to ambiguous alignment, 42 characters (7074, 102104, 249258, 268285, 388389, and 409412) were excluded. The final 433-bp matrix contained 264 parsimony-informative characters. Parsimony analysis resulted in 144 equally most parsimonious trees (MPTs) of 1218 steps. The uncorrected consistency index (CI) was 0.48 and the corrected CI was 0.44, the latter of which was higher than the expected CI (0.31) for 52 taxa (Sanderson and Donoghue, 1989
). The retention index (RI) was 0.72. The strict consensus tree of the 144 MPTs with bootstrap support values (BV) and decay indices (DI) for each clade is presented in Fig. 2. The g1 statistic for 10 000 random trees was 0.39, which is significantly more skewed than random data (g1 = 0.09 for 250 variable position and 25 or more taxa; P < 0.01; Hillis and Huelsenbeck, 1992
).
Compared to the results of Downie and Katz-Downie (1996)
, whose sampling strategy was basically adopted here, relationships among taxa that were both sampled by previous and current studies are highly concordant (Fig. 2). One of the major differences was the placement of Smyrnium olusatrum, which was placed sister to tribe Scandiceae with strong support in Downie and Katz-Downie (1996)
but was placed next to tribe Oenantheae with weak support in the current study (Fig. 2).
The eight Oreomyrrhis species formed a well-supported monophyletic group (BV = 93%; DI = 3) with low internal resolution (Fig. 2). This Oreomyrrhis clade was sister to Chaerophyllum temulum, and together they formed a strongly supported clade (BV = 99%; DI = 8) within subtribe Scandicinae. None of the genera that were assigned to tribe Smyrnieae by Drude (1898)
and other genera with simple and/or irregular compound umbels (Fig. 2) are closely related to Oreomyrrhis.
Analysis 2
Alignment of 26 Oreomyrrhis haplotypes, 22 Chaerophyllum species (sequences of Chaerophyllum elegans/ C. hirsutum and C. hakkiaricum/C. macrospermum are identical, respectively), and 11 outgroups resulted in a matrix of 57 terminals and 452 characters. Twenty-nine gaps of various sizes (nineteen 1-bp, four 2-bp, one 3-bp, two 4-bp, one 8-bp, one 11-bp, and one 16-bp) were introduced to facilitate alignment. Ten of these gaps were parsimony informative (gap A: position 68; B: 72; C: 101; D: 237, E: 255; F: 375377; G: 411; H: 247262; I: 266269; J: 344345; Fig. 3).
Treating gaps as missing data, MP analysis of the ITS data matrix resulted in 2680 MPTs of 463 steps (CI = 0.70, corrected CI = 0.62; RI = 0.84). The g1 statistic for 10 000 random trees was 0.46, which is significantly more skewed than that of random data (Hillis and Huelsenbeck, 1992
). To test the phylogenetic utility of the alignment gaps, parsimony-informative gaps were scored as binary characters, resulting in a second matrix of 462 aligned positions. Maximum parsimony analysis revealed 1580 MPTs of 490 steps with slightly lower consistency and retention indices (CI = 0.68, corrected CI = 0.60; RI = 0.83), suggesting an increase in homoplasy when gaps were taken into account. The strict consensus trees derived from the two matrices (with and without gap characters), however, were almost identical, except for two additional clades (denoted by asterisks in Fig. 3) uncovered in the gap-including matrix. The strict consensus tree derived from the gap-including matrix is presented in Fig. 3. Mapping the 10 parsimony informative gaps onto the consensus tree indicated that only three gaps were free of homoplasy (Fig. 3).
Modeltest 3.5 selected the GTR + G model (base frequencies: 0.2403, A; 0.2182, C; 0.2480, G; 0.2935, T; estimates of substitution rates A
C: 1.3354; A
G: 2.3532; A
T: 1.4904; C
G: 0.4394; C
T: 3.1987; G
T: 1; proportion of invariable sites = 0; gamma distribution shape parameter = 0.9001). Based on these parameters, a single maximum likelihood (ML) tree was uncovered (Ln likelihood score = 3153.10211; Fig. 4A). The same model and parameters were used to perform a Bayesian analysis. The 50% majority rule consensus tree with Bayesian posterior probabilities (PP) is shown in Fig. 4B.
Topologies uncovered from the three phylogenetic inference methods were very similar, with only minor differences in resolution between some terminal taxa (Figs. 3 and 4). In all analyses, monophyly of the Chaerophyllum-Oreomyrrhis clade was well to strongly supported (BV = 85% in MP and 91% in ML; DI = 3; PP = 99%). Within this clade, four strongly supported clades were uncovered. Three of these corresponded perfectly to Chaerophyllum sect. Chrysocarpum Spalik & S. R. Downie, sect. Dasypetalon Neilr., and sect. Physocaulis DC. [a monotypic section composed of only C. nodosum (L.) Crantz], respectively, that were identified by recent molecular systematic studies (Spalik and Downie, 2001
; Spalik et al., 2001a
). Section Chaerophyllum, however, was rendered paraphyletic by Oreomyrrhis.
In all analyses, Oreomyrrhis was placed within Chaerophyllum sect. Chaerophyllum with moderate to strong support (BV = 76% in MP and 88% in ML; DI = 2; PP = 99%). Within this clade, all 26 Oreomyrrhis accessions and the two North American Chaerophyllum species (C. procumbens and C. tainturieri) formed a strongly supported subclade (BV = 88% in MP and 93% in ML; DI = 2; PP = 90%) that was sister to C. temulum. Relationships within the subclade, however, were polytomous, with four moderate to strongly supported clades nested within it. These four clades are the North American Chaerophyllum clade (BV = 99% in MP and 98% in ML; DI = 5; PP = 100%), the New Zealand Oreomyrrhis clade (BV = 80% in MP and 88% in ML; DI = 2; PP = 99%), the East Asian Oreomyrrhis clade (BV = 74% in MP and 76% in ML; DI = 1; PP = 96%), and the Australian endemic clade consisting of several Australian Oreomyrrhis taxa (BV = 60% in MP and 59% in ML; DI = 1; PP = 89%).
In the ML analysis (Fig. 4A), all Oreomyrrhis accessions were grouped together in a weakly supported monophyletic group (BV = 57%) that was sister to the North American Chaerophyllum clade. However, supports for the monophyletic Oreomyrrhis were weak in both MP and Bayesian analysis (BV = 36% in MP; PP = 38%). Whether or not the Oreomyrrhis species that were represented by multiple accessions were monophyletic was inconclusive.
Analysis 3
Alignment of the 57 Oreomyrrhis accessions and the two North American Chaerophyllum species resulted in a matrix of 611 positions (Table 1). Alignment position 568 was excluded due to its potential homoplasy (Table 1; see Discussion). With alignment gaps treated as a fifth state, a single, fully resolved network was reconstructed (Fig. 5). The ITS network has a starlike topology, within which the most frequent haplotype (common Oreomyrrhis haplotype) is linked to the outgroup taxa (Chaerophyllum clade) and connected to 12 relatively low frequency haplotypes and clades, including the New Zealand clade, the East Asian clade, and the Australian clade that were also recovered in analysis 2 (Figs. 3 and 4). For species represented by multiple haplotypes (O. argentea, O. ciliata, O. colensoi, O. eriopoda, and O. ramosa), these constituent haplotypes were not grouped into separate haplotype lineages (Fig. 5).
|
|
Utilization of nrDNA ITS region for phylogenetic inference
The properties of ITS loci that have attracted and facilitated their application in plant phylogenetic studies (Baldwin et al., 1995
) were recently reexamined by Álvarez and Wendel (2003)
. Although ITS loci are generally assumed to be homogenized by concerted evolution, Álvarez and Wendel (2003)
cautioned that, under circumstances such as hybridization, polyploidization, or pseudogene formation, separate ITS sequences might be present and persistent in a genome. This could interfere with sequencing, distort assessment of sequence homology, and subsequently mislead phylogenetic inference.
The ITS region was chosen primarily because it is the only molecular marker for which a comprehensive sampling of Apiaceae (Downie et al., 2000
, 2001
) is available for testing taxonomic hypotheses regarding the tribal placement of Oreomyrrhis. The ITS region is also the fastest evolving sequence known in Apiaceae (Downie et al., 2001
) and therefore is more likely to resolve infrageneric relationships. In Oreomyrrhis, a low base chromosome number precludes its origin through polyploidization. The autogamous nature (discussed later) and currently isolated distribution of most Oreomyrrhis species assure that interspecific hybridization and introgression are probably rare. Plant species characterized by these biological attributes are frequently assumed to have had their ITS paralogues homogenized by concerted evolution (Denduangboripant and Cronk, 2000
).
However, molecular cloning of all 10 selected Oreomyrrhis accessions revealed intra-individual ITS polymorphisms that were imperceptible from the electropherogram of direct sequencing. The persistence of intra-individual ITS polymorphism in Oreomyrrhis suggests that concerted evolution might not have proceeded fast enough to homogenize different ITS sequence types (Denduangboripant and Cronk, 2000
; Hughes et al., 2002
). One possible scenario for the failure of concerted evolution in this non-hybrid diploid genus is the existence of multiple nrDNA arrays in the genome (O'Kane et al., 1996
), which can be detected by fluorescent in situ hybridization (FISH) and genomic in situ hybridization (GISH) (Álvarez and Wendel, 2003
). Alternatively, the existence of different ITS paralogues within each accession of Oreomyrrhis may suggest that intra-individual ITS polymorphism might be more widespread than previously assumed (Hughes et al., 2002
). Nevertheless, since the revealed intra-individual ITS polymorphisms present shallow paralogy that does not impede phylogenetic inference, ITS sequences remain a credible molecular marker for phylogenetic reconstruction in Oreomyrrhis.
Phylogenetic affinities of Oreomyrrhis
The ITS data presented here unambiguously placed Oreomyrrhis in the subtribe Scandicinae of tribe Scandiceae (Fig. 2), supporting the taxonomic hypotheses of de Candolle (1830)
, Bentham (1867)
, and Koso-Poljansky (1916)
, and rejecting Drude's (1898)
placement of Oreomyrrhis in tribe Smyrnieae. Drude's tribe Smyrnieae (encompassing Arracacia, Conium, Erigenia, Neogoezia, Orogenia, Oreomyrrhis, Pleurospermum, and Smyrnium) is polyphyletic in the present analysis (Fig. 2), confirming relationships suggested by previous studies using serological (Shneyer et al., 1992
) and molecular data (Downie and Katz-Downie, 1996
; Katz-Downie et al., 1999
). Our ITS phylogenies also indicate that taxa characterized by simple or irregular compound umbels (e.g., Oreomyrrhis, Orogenia, Neogoezia, Lilaeopsis, and Erigenia) do not form a monophyletic group (Fig. 2), indicating that this trait has had multiple independent origins in Apioideae and that its phylogenetic utility has been overemphasized.
The close relationship between Oreomyrrhis and Chaerophyllum proposed by Bentham (1866)
was verified by our ITS data (Figs. 2 4). Chaerophyllum is the largest (ca. 30 species) and most diverse genus in subtribe Scandicinae. It is distributed mainly in Eurasia (including North Africa) and is most diversified in the Mediterranean and the Caucasus regions, with two species native to North America (Spalik and Downie, 2001
). Recent molecular studies support the monophyly of Chaerophyllum and reveal four well-supported clades within it (Spalik and Downie, 2001
). Lacking morphological synapomorphies, these four clades were each treated as section by Spalik and Downie (2001)
.
Based on the ITS data, Oreomyrrhis is nested within the sect. Chaerophyllum (Figs. 3 and 4) that comprises one species native to Europe (C. temulum) and two to North America (C. procumbens and C. tainturieri). Within this section, Oreomyrrhis is most closely related to the two North American species. The phylogenetic relationships revealed by the current data indicate the paraphyletic nature of Chaerophyllum if Oreomyrrhis is treated as a separate genus. Because C. temulum is the type species of the genus, preservation of the generic status of Oreomyrrhis would simultaneously mandate the erection of four small genera (the North American Chaerophyllum clade, sect. Chrysocarpum, sect. Dasypetalon, and the monotypic sect. Physocaulis), leaving Chaerophyllum a monotypic genus composed of only C. temulum. Before proposing formal nomenclatural changes, further study by employing additional molecular markers of chloroplast genome will be conducted to test the result reached by ITS data.
Monophyly of Oreomyrrhis
The monophyly of Oreomyrrhis was only weakly supported in the ML analysis (BV = 57%; Fig. 4A). The bootstrap support in MP analysis and Bayesian posterior probability for this clade are both less than 50%. Examination of the aligned ITS matrix revealed a likely homoplasious character at position 409 (position 568 in Table 1 and analysis 3). To evaluate how this character had affected phylogeny reconstruction, we excluded this character and performed a MP analysis. Using the same parameters for the gap-excluding matrix in analysis 2, a heuristic search resulted in 536 equal length MPTs of 458 steps, two steps shorter than the original trees (460 steps). We used the "fast stepwise addition" option of PAUP* to perform a bootstrap analysis with 100 000 replicates (e.g., Sun et al., 2004
). The strict consensus of the 536 MPTs was almost identical to the ML tree (Fig. 4A), and the support for the monophyly of Oreomyrrhis accessions increased to 66%, compared to 36% when character 568 was included.
Although the ITS data provide only weak supports, the monophyly of Oreomyrrhis is unambiguously defined by its morphological synapomorphy (simple umbel) and apparent cytological and ecological distinctiveness. The monophyly indicates a single origin of Oreomyrrhis, suggesting that the trans-Pacific disjunction of the genus reflects a genuine biogeographic process, not a taxonomic artifact.
Biogeographic implications of the ITS phylogenies
Because a great proportion of the distribution range of Oreomyrrhis parallels the Gondwanan disjunction in the South Pacific (e.g., Sanmartín and Ronquist, 2004
), this genus had often been assumed to have a Southern Hemispheric and/or ancient origin (e.g., van Steenis, 1963
; Dawson, 1963
; Raven, 1973
; Melville, 1981
; Smith, 1986
; Wu, 1998
). However, our ITS data strongly support that Oreomyrrhis is closely related to Chaerophyllum of subtribe Scandicinae, that is mainly distributed in the Northern Hemisphere. Within Chaerophyllum, Oreomyrrhis is nested within sect. Chaerophyllum and, together with the two North American species of the section (C. procumbens and C. tainturieri), forms a clade that is sister to the European species C. temulum (Figs. 3 and 4). Spalik and Downie (2001)
hypothesized that Chaerophyllum arrived in North America by the incidental dispersion of seeds from Europe, probably by vagrant birds. Based on the optimization of DIVA, the ancestral area for sect. Chaerophyllum (including Oreomyrrhis) was suggested to be Eurasia, and the disjunction in North American and Pacific involved two dispersal events (Fig. 6). Because the distribution range of Oreomyrrhis and Chaerophyllum do not overlap, DIVA optimization indicates that the current distribution could have resulted from dispersal from Europe to North America followed by a second dispersal to Pacific, to Pacific first and second dispersal to North America, or two independent dispersals to Europe and Pacific (Fig. 6).
Based on our time estimates, these dispersal events probably took place around late Miocene to early Pliocene ([5.18] 8.16 [11.15] mya), when dispersals from Europe to North America through island hopping via the North Atlantic Land Bridge (Tiffney and Manchester, 2001
) were still possible. However, since the basal lineage of Oreomyrrhis cannot be identified confidently by our data, the initial dispersal event that gave rise to subsequent Oreomyrrhis diversification remains unclear.
The estimated times for the separation of the North American Chaerophyllum and Oreomyrrhis (2.23 mya [95% confidence intervals = 0.993.47 mya]) and diversification within Oreomyrrhis (1.06 mya [95% confidence intervals = 0.56 1.47 mya]) are relatively recent, corresponding to a late Tertiary to Quaternary evolution of the group. Note that the 0.79% sequence divergence rate is at the lower end of the range reported for flowering plants (Zhang et al., 2001
), indicating that the dates inferred here might have been overestimated and the actual times for Oreomyrrhis diversification might be more recent. Even under this more conservative estimate, our results strongly support a late Tertiary to Quaternary origin and diversification of Oreomyrrhis, and most likely long-distance dispersal as the major biogeographic mechanism underlying the trans-Pacific disjunction of Oreomyrrhis (Dawson, 1971
; Thorne, 1972
; Raven, 1973
). The vicariant scenarios that invoke early Tertiary land-bridge (van Steenis, 1962
, 1964
) or plate tectonic movements dating back to the Jurassic (Melville, 1981
) can be rejected confidently because their predicted temporal contexts significantly predate the origin of Oreomyrrhis inferred from the ITS data.
Recent phylogenetic studies suggest that long-distance dispersal underlies the distributions of many trans-Pacific disjunct plant groups (Winkworth et al., 2002
; Sanmartín and Ronquist, 2004
). Comparisons of these Pacific taxa have revealed concordant phylogenetic relationships, suggesting that recent long-distance dispersal events might not occur randomly but have followed certain dispersal routes, probably associated with ancient vegetation types that may have served as migratory corridors (Winkworth et al., 2002
; Sanmartín and Ronquist, 2004
). Unfortunately, the ITS data contain only limited information regarding the pattern of dispersal that have shaped the trans-Pacific disjunction of Oreomyrrhis.
Biogeographic implications of the ITS haplotype network
To better understand genealogical processes that underlying the Pacific disjunction, an ITS haplotype network was constructed (Fig. 5). The predominant feature of the network is a single widespread and high frequency haplotype (common Oreomyrrhis haplotype) linked to numerous low-frequency haplotypes by a few mutational steps, a pattern that is commonly seen in temperate taxa in Europe and North America that have experienced drastic distributional changes during the Pleistocene (Hewitt, 1996
; Schaal et al., 1998
; Milá et al., 2000
). This type of network structure could have resulted from (1) a postglacial range expansion from refugium populations that had experienced genetic bottleneck, (2) gene flow, or (3) a selective sweep (Milá et al., 2000
). In Oreomyrrhis, trans-Pacific gene flow seems an unlikely explanation for the wide distribution of the common Oreomyrrhis haplotype (Abbott and Comes, 2004
). Because the ITS region is known to subject little functional constraint (Baldwin et al., 1995
), unless this region is closely linked to a gene favored by natural selection, a selective sweep is a less likely explanation than range expansion (Milá et al., 2000
).
However, unlike low-elevation temperate taxa whose distribution ranges expanded after the retreat of the ice sheet (Hewitt, 1996
), alpine plants such as Oreomyrrhis would have actually experienced postglacial range contraction, as rising global temperatures forced alpine plants to migrate up toward mountain tops (Raven, 1973
; Larena et al., 2002
). During the glacial maxima, alpine plant populations would have actually expanded as the alpine zone moved downward. This reverse scenario of glacial population expansion and interglacial population contraction is known as the "displacement refugia model" (Kropf et al., 2003
). Under this model, the common Oreomyrrhis haplotype likely achieved its broad distribution by "glacial" range expansion as more suitable alpine habitats became available during glacial maxima. Since the common Oreomyrrhis haplotype is the only haplotype distributed in more than one geographic region, its glacial range expansion probably had occurred before other haplotypes arose.
Species circumscription and speciation in Oreomyrrhis
Although many of the species circumscribed by Mathias and Constance (1955)
and Chen and Wang (2001)
are characterized by distinct ITS haplotypes (e.g., O. azorellacea, O. borneensis, O. brevipes, O. gunnii, O. linearis, O. nanhuensis, and O. papuana), suggesting that they have evolved as separate evolutionary lineages (Templeton et al., 2000
), such genetic distinctness is absent in the majority of the species investigated (Fig. 5). For instance, the common Oreomyrrhis haplotype is carried by nine species that are distributed throughout most of the range of Oreomyrrhis (Fig. 5). The sharing of haplotypes by recently differentiated species could result from interspecific hybridization or incomplete lineage sorting of ancestral haplotypes that had been present prior to the divergence of the species (Olsen and Schaal, 1999
). In the current study, the ancestral position of the common Oreomyrrhis haplotype suggested by its high mutational connectedness and high frequency (Castelloe and Templeton, 1994
) strongly favors lineage sorting, for interspecific hybridizations between disjunct taxa across Pacific seem unlikely (Abbott and Comes, 2004
). On the other hand, gene flows among species in the same geographic area (e.g., O. argentea A, O. eriopoda B, and O. pulvinifica in SE Australia; O. involucrata and O. taiwaniana in Taiwan; O. colensoi A, O. ramosa A, and O. rigida in New Zealand) are possible for haplotype sharing, especially during glacial maxima when alpine plant populations underwent range expansion into lower elevations (Larena et al., 2002
). Furthermore, the recency of species divergences within Oreomyrrhis suggests that there probably has not been enough time for achieving reciprocal monophyly of the newly evolved lineages (Avise, 2000
). The preliminary results from ITS data provide an exciting first step toward an in-depth understanding of the speciation in Oreomyrrhis. In the ongoing phylogeographic and population genetic studies, we are employing chloroplast sequence data and AFLP markers to explore the intriguing evolutionary history and biogeography of Oreomyrrhis.
Cytological and morphological considerations
The base chromosome number for the genus Chaerophyllum is n = 11, with deviant and intraspecific variable counts only being reported in sect. Chaerophyllum (C. temulum [2n = 14 and 22] and C. procumbens [2n = 12 and 22]) (Pimenov et al., 2003
). Mapping base chromosome numbers onto the ITS phylogeny (Fig. 3) indicates that the numbers n = 6 and n = 7 occur only in the sect. Chaerophyllum clade (including Oreomyrrhis).
Although mericarps of most species of Oreomyrrhis are not equipped with structures such as hooks or nutrient tissues (Fig. 7) that could facilitate dispersal (Mathias and Constance, 1955
; van Steenis, 1963
), Raven (1973)
suggested that, during glacial maxima when high latitude alpine plants were forced to migrate to lower elevations closer to sea levels, chances for seed dispersals by transoceanic wanderers (e.g., albatross) would greatly increase. Compared to other species in Chaerophyllum, members of sect. Chaerophyllum (including Oreomyrrhis) possess relatively smaller mericarps (Spalik et al., 2001b
), with a further size reduction in Oreomyrrhis (Fig. 7). The evolution toward smaller mericarps in Oreomyrrhis might have been advantageous in enhancing their dispersibility (Howe and Smallwood, 1982
), especially when carried in mud adhering to the legs of wanderers (Salisbury, 1970
; Raven, 1973
).
|
Despite inbreeding depression, the genetic benefits (e.g., reproductive assurance and purge of genetic load) associated with the shift from outcrossing to selfing promote successful long-distance colonization in plants, an observation known as Baker's Rule (Baker, 1955
). Recent studies suggest that a rapid transition from outcrossing to self-pollination in the early evolution of Arabidopsis thaliana could have contributed to its rapid postglacial range expansion around 17 000 years ago (Shimizu et al., 2004
). In conjunction with these lines of evidence, our results suggest that the evolution in chromosome number and transition in reproductive syndrome, along with reduction in mericarp size, probably played pivotal roles during the early evolution of sect. Chaerophyllum, followed by further diversification and occupation of cold and alpine environments by Oreomyrrhis.
Downie, S. R., and D. S. Katz-Downie. 1996. A molecular phylogeny of Apiaceae subfamily Apioideae: evidence from nuclear ribosomal DNA internal transcribed spacer sequences. American Journal of Botany 83: 243 251.
Downie, S. R., D. S. Katz-Downie, and K. Spalik. 2000. A phylogeny of Apiaceae tribe Scandiceae: evidence from nuclear ribosomal DNA internal transcribed spacer sequences. American Journal of Botany 87: 7695.
Downie, S. R., R. L. Hartman, F.-J. Sun, and D. S. Katz-Downie. 2002. Polyphyly of the spring-parsleys (Cymopterus): molecular and morphological evidence suggests complex relationships among the perennial endemic genera of western North American Apiaceae. Canadian Journal of Botany 80: 12951324.
Downie, S. R., S. Ramanath, D. S. Katz-Downie, and E. Llanas. 1998. Molecular systematics of Apiaceae subfamily Apioideae: phylogenetic analyses of nuclear ribosomal DNA internal transcribed spacer and plastid rpoC1 intron sequences. American Journal of Botany 85: 563591.
Hardway, T. M., K. Spalik, M. F. Watson, D. S. Katz-Downie, and S. R. Downie. 2004. Circumscription of Apiaceae tribe Oenantheae. South African Journal of Botany 70: 393406.
Katz-Downie, D. S., C. M. Valiejo-Roman, E. I. Teretieva, A. V. Troitsky, M. G. Pimenov, B. Lee, and S. R. Downie. 1999. Towards a molecular phylogeny of Apiaceae subfamily Apioideae: additional information from nuclear ribosomal DNA ITS sequences. Plant Systematics and Evolution 216: 167195.
Sun, F.-J., S. R. Downie, and R. L. Hartman. 2004. An ITS-based phylogenetic analysis of the perennial endemic Apiaceae subfamily Apioideae of western North America. Systematic Botany 29: 419431.
FOOTNOTES
The authors thank P. Raven, M. Richardson, and A. Larson for guidance; J. Beck, B. Torke, P. Stevens, and members of the Schaal lab at Washington University for valuable discussions and comments on the manuscript; C.-H. Chen for sharing unpublished ITS data; S. Avendaño R., N. Goodall, D. Broughton, H.-M. Chang, L.-J. Juan, A. Mark, and G. Plunkett for plant materials; A, B, CANB, F, G, GH, HO, K, L, LA, LAE, NSW, NY, P, UC, US, VCU, W, and WELT for loans; A, CANB, HAST, HO, MO, and UC for permission to sample herbarium material; H.-F. Kuo for technical assistance in cloning; and M. Veith for instruction on SEM. KFC is deeply grateful for the hospitalities, assistances, and logistic supports from friends and institutions worldwide: Vargas family, and QCNE staff in Ecuador; C.-I Huang, Y.-Y. Huang, W.-C. Leong, C.-H. Liu, and HAST staff in Taiwan; K. Wilson (NSW), J. Bruhl (UNE), B. Lepschi (CANB), M. Fagg (CANB), L. Halasz (CANB), N. Walsh (MEL), A. Buchanan (HO), and M. Song in Australia; P. Heenan (CHR), CHR staff, S. Halloy (Otago), and A. Anderberg in New Zealand; and A. MacVean (UVAL) in Guatemala. Permission for fieldwork from the following authorities is sincerely appreciated: Ministry of Environment of Ecuador; Yushan, Taroko, and Shei-pa National Parks of Taiwan; Environment ACT, NSW National Parks and Wildlife Service, Victoria Department of Sustainability and Environment; Tasmania Department of Environment, Forestry Tasmania, Hydro Tasmania; New Zealand Department of Conservation; and Consejo Nacional de Areas Protegidas (Guatemala). Field and laboratory work was supported by the West Pacific Phylogeography Project of Academia Sinica, Taiwan (CIP and KFC), a National Science Foundation Doctoral Dissertation Improvement Grant (DEB-0408105), BSA Karling Award, ASPT Graduate Student Research Grant, a grant from the Mellon Foundation to the Missouri Botanical Garden, and the DBBS at Washington University. ![]()
7 Author for correspondence (kchung{at}artsci.wustl.edu
) ![]()
LITERATURE CITED
Abbott R. J. H. P. Comes 2004 Evolution in the Arctic: a phylogeographic analysis of the circumarctic plant, Saxifraga oppositifolia (purple saxifrage). New Phytologist 161: 211-224[CrossRef][Web of Science]
Allan H. H. 1961 Flora of New Zealand, vol. 1. R. E. Owen Government Printer, Wellington, New Zealand
Álvarez I. J. F. Wendel 2003 Ribosomal ITS sequences and plant phylogenetic inference. Molecular Phylogenetics and Evolution 29: 417-434[CrossRef][Web of Science][Medline]
Avise C. J. 2000 Phylogeography: the history and formation of species. Harvard University Press, Cambridge, Massachusetts, USA
Bailey C. D. T. G. Carr S. A. Harris C. E. Hughes 2003 Characterization of angiosperm nrDNA polymorphism, paralogy, and pseudogenes. Molecular Phylogenetics and Evolution 29: 435-455[CrossRef][Web of Science][Medline]
Baldwin B. G. M. J. Sanderson J. M. Porter M. F. Wojciechowski C. S. Campbell M. J. Donoghue 1995 The ITS region of nuclear ribosomal DNA: a valuable source of evidence on angiosperm phylogeny. Annals of the Missouri Botanical Garden 82: 247-277[CrossRef][Web of Science]
Baker H. G. 1955 Self-compatibility and establishment after "long-distance" dispersal. Evolution 9: 347-349[CrossRef][Web of Science]
Baskin C. C. T. S. Hawkins J. M. Baskin 2004 Ecological life cycle of Chaerophyllum procumbens variety shortii (Apiaceae), a winter annual of the North American eastern deciduous forest. Bulletin of the Torrey Botanical Club 131: 126-139
Baskin J. M. C. C. Baskin 1990 Germination ecophysiology of seeds of the winter annual Chaerophyllum tainturieri: a new type of morphophysiological dormancy. Journal of Ecology 78: 993-1004[CrossRef]
Bentham G. 1866 Flora australiensis: a description of the plants of the Australian territory, vol. 3. Lovell Reeve, London, UK
Bentham G. 1867 Umbelliferae. In G. Bentham and J. D. Hooker [eds.], Genera plantarum, vol. 1, 859931. Lovell Reeve, London, UK
Blattner F. R. 1999 Direct amplification of the entire ITS region from poorly preserved plant material using recombinant PCR. BioTechniques 29: 1180-1186
Bremer K. 1988 The limits of amino acid sequence data in angiosperm phylogenetic reconstruction. Evolution 42: 795-803[CrossRef][Web of Science]
Buwalda P. 1951 Notes on New Guinea Umbelliferae. Journal of Arnold Arboretum 32: 59-66
Castelloe J. A. R. Templeton 1994 Root probabilities for intraspecific gene trees under neutral coalescent theory. Molecular Phylogenetics and Evolution 3: 102-113[CrossRef][Medline]
Chen C.-H. J.-C. Wang 2001 Revision of the genus Oreomyrrhis Endl. (Apiaceae) in Taiwan. Botanical Bulletin of Academia Sinica 42: 303-312[Web of Science]
Clement M. D. Posada K. A. Crandall 2000 TCS: a computer program to estimate gene genealogies. Molecular Ecology 9: 1657-1659[CrossRef][Medline]
Constance L. 1987 Neogoezia (Apiaceae), a very distinct and elegant genus of Mexican Umbelliferae. Opera Botanica 92: 59-71
Crisci J. V. M. M. Cigliano J. J. Morrone S. Roig-Junent 1991 Historical biogeography of southern South America. Systematic Zoology 40: 152-171[CrossRef]
Dawson J. W. 1963 Origins of New Zealand alpine flora. New Zealand Journal of Ecology 10: 12-15
Dawson J. W. 1971 Relationships of the New Zealand Umbelliferae. In V. H. Heywood [ed.], The biology and chemistry of the Umbelliferae, 43 61. Academic Press, New York, New York, USA
de Candolle A. P. 1829 Mémoire sur la famille des Ombellifères. Treuttel et Würtz, Paris, France
de Candolle A. P. 1830 Umbelliferae. In A. P. de Candolle [ed.], Prodromus systematis naturalis regni vegetabilis, vol. 4, 55250. Treuttel et Würtz, Paris, France
Denduangboripant J. Q. C. B. Cronk 2000 High intraindividual variation in internal transcribed spacer sequences in Aeschynanthus (Gesneriaceae): implication for phylogenetics. Proceedings of the Royal Society of London, series B 267: 1407-1415[Medline]
Downie S. R. D. S. Katz-Downie 1996 A molecular phylogeny of Apiaceae subfamily Apioideae: evidence from nuclear ribosomal DNA internal transcribed spacer sequences. American Journal of Botany 83: 243-251
Downie S. R. D. S. Katz-Downie K. Spalik 2000 A phylogeny of Apiaceae tribe Scandiceae: evidence from nuclear ribosomal DNA internal transcribed spacer sequences. American Journal of Botany 87: 76-95
Downie S. R. G. M. Plunkett M. F. Watson K. Spalik D. S. Katz-Downie C. M. Valiejo-Roman E. I. Terentieva A. V. Troitsky B.-Y. Lee J. Lahham A. El-Oqlah 2001 Tribes and clades within Apiaceae subfamily Apioideae: the contribution of molecular data. Edinburgh Journal of Botany 58: 301-330[CrossRef]
Downie S. R. F.-J. Sun D. S. Katz-Downie G. J. Colletti 2004 A phylogenetic study of Perideridia (Apiaceae) based on nuclear ribosomal DNA ITS sequences. Systematic Botany 29: 737-751[CrossRef][Web of Science]
Doyle J. J. J. L. Doyle 1987 A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochemistry Bulletin 19: 11-15
Drude O. 1898 Umbelliferae. In A. Engler and K. Prantl [eds.], Die natürlichen Pflanzenfamilien, vol. 3, 63250. Wilhelm Engelmann, Leipzig, Germany
Endlicher S. 1839 Genera plantarum secundum ordines naturals disposita. Vindobonae, Apud Fr. Beck Universitatis Bibliopolam, Wien, Austria
Felsenstein J. 1985 Confidence limits on phylogenies: an approach using the bootstrap. Evolution 39: 783-791[CrossRef][Web of Science]
Hayata B. 1911 Materials for a flora of Formosa. Journal of the College of Science, Imperial University of Tokyo 30: 1-471
Hershkovitz M. A. E. A. Zimmer W. J. Hahn 1999 Ribosomal DNA sequences and angiosperm systematics. In P. M. Hollingsworth, R. M. Bateman, and R. J. Gornall [eds.], Molecular systematics and plant evolution, 268326. Taylor & Francis, London, UK
Hewitt G. M. 1996 Some genetic consequences of ice ages and their role in divergence and speciation. Biological Journal of the Linnean Society 58: 247-276[CrossRef]
Hillis D. M. J. P. Huelsenbeck 1992 Signal, noise, and reliability in molecular phylogenetic analyses. Journal of Heredity 83: 189-195
Holmgren P. K. N. H. Holmgren L. C. Barnett 1990 Index herbariorum. New York Botanical Garden, New York, New York, USA
Hooker J. D. 1856 The botany of the Antarctic voyage of H. M. Discovery Ships Erebus and Terror in the year 18391843. II. Flora Novae-Zelandiae. Lovell-Reeve, London, UK
Howe H. J. J. Smallwood 1982 Ecology of seed dispersal. Annual Review of Ecology and Evolution 13: 201-228
Hughes C. E. C. D. Bailey S. A. Harris 2002 Divergent and reticulate species relationships in Leucaena (Fabaceae) inferred from multiple data sources: insights into polyploidy origins and nrDNA polymorphism. American Journal of Botany 87: 1057-1073
Humboldt A. A. Bonpland C. S. Kunth 1821 Nova genera et species plantarum, vol. 5. Maze, Paris, France
Jermiin L. S. 1996 K2WuLi, version 1.0. Australian National University, Canberra, Australia. Available at website, http://jcsmr.anu.edu.au/dicb/humgen/lars/k2wulisub.htm/
Johnston I. M. 1938 Some undescribed species from Mexico and Guatemala. Journal of the Arnold Arboretum 19: 117-128
Katz-Downie D. S. C. M. Valiejo-Roman E. I. Teretieva A. V. Troitsky M. G. Pimenov B. Lee S. R. Downie 1999 Towards a molecular phylogeny of Apiaceae subfamily Apioideae: additional information from nuclear ribosomal DNA ITS sequences. Plant Systematics and Evolution 216: 167-195
Kern J. H. 1958 Umbelliferae. In C. G. G. J. van Steenis [ed.], Flora Malesiana, series I, vol. 5, part 4, 555557. Noordhoff-Kolff N. V., Djakarta, Indonesia
Kimura M. 1980 A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. Journal of Molecular Evolution 16: 111-120[CrossRef][Web of Science][Medline]
Kirk T. 1899 The student flora of New Zealand. J. Mackay, Government Printer, Wellington, New Zealand
Körner C. 2003 Alpine plant life: functional plant ecology of high mountain ecosystems, 2nd ed. Springer-Verlag, Berlin, Germany
Koso-Poljansky B. M. 1916 Sciadophytorum systematis lineamenta. Bulletin de la Société Impériale des Naturalistes (Moscou) 29: 93-222
Kropf M. J. W. Kadereit H. P. Comes 2002 Late Quaternary distributional stasis in the submediterranean mountain plant Anthyllis montana L. (Fabaceae) inferred from ITS sequences and amplified fragment length polymorphism markers. Molecular Ecology 11: 447-463[CrossRef][Medline]
Kropf M. J. W. Kadereit H. P. Comes 2003 Differential cycles of range contraction and expansion in European high mountain plants during the late Quaternary: insights from Pritzelago alpina (L.) O. Kuntze (Brassicaceae). Molecular Ecology 12: 931-949[CrossRef][Medline]
Kumar S. K. Tamura I. B. Jakobsen M. Nei 2001 MEGA2: molecular evolutionary genetics analysis software. Bioinformatics 17: 1244-1245
La Gasca M. 1821 Amenidades naturals de las Españas, vol. I, part 2. Ibarra, Madrid, Spain
Larena B. G. J. F. Aguilar G. N. Feliner 2002 Glacial-induced altitudinal migration in Armeria (Plumbaginaceae) inferred from patterns of chloroplast DNA haplotype sharing. Molecular Ecology 11: 1965-1974[CrossRef][Medline]
Mathias M. E. 1965 Distribution patterns of certain Umbelliferae. Annals of the Missouri Botanical Garden 52: 387-398[CrossRef]
Mathias M. E. L. Constance 1955 The genus Oreomyrrhis (Umbelliferae), a problem in South Pacific distribution. University of California Publications in Botany 27: 347-416
Mathias M. E. L. Constance 1977 A new species of Oreomyrrhis (Umbelliferae, Apiaceae) from New Guinea. Journal of Arnold Arboretum 58: 190-192
Melville R. 1981 Vicarious plant distributions and paleogeography of the Pacific region. In G. Nelson and D. E. Rosen [eds.], Vicariance biogeography: a critique, 238274. Columbia University Press, New York, New York, USA
Merrill E. D. 1918 Oreomyrrhis borneensis Merr. sp. nov., an interesting addition to our knowledge of the Malayan flora. American Journal of Botany 5: 514-515[CrossRef][Web of Science]
Milá B. D. J. Girman M. Kimura T. B. Smith 2000 Genetic evidence for the effect of postglacial population range expansion on the phylogeography of a North American songbird. Proceedings of the Royal Society of London, series B 267: 1033-1040[Medline]
Molloy B. P. E. Edgar P. B. Heenan P. J. de Lange 1999 New species of Poa (Gramineae) and Ischnocarpus (Brassicaceae) from limestone, North Otago, South Island, New Zealand. New Zealand Journal of Botany 37: 41-50[Web of Science]
Moore D. M. 1968 The vascular flora of the Falkland Islands. British Antarctic Survey Scientific Reports, no. 60. British Antarctic Survey, Natural Environment Research Council, London, UK
Moore D. M. 1971 Chromosome studies in Umbelliferae. In V. H. Heywood [ed.], The biology and chemistry of the Umbelliferae, 233255. Academic Press, New York, New York, USA
Neves S. S. M. F. Watson 2004 Phylogenetic relationships in Bupleurum (Apiaceae) based on nuclear ribosomal DNA ITS sequence data. Annals of Botany 93: 379-398
Olsen K. O. B. A. Schaal 1999 Evidence on the origin of cassava: phylogeography of Manihot esculenta. Proceedings of the National Academy of Sciences, USA 96: 5586-5591
O'Kane S. L. B. A. Schaal I. Al-Shehbaz 1996 The origins of Arabidopsis suecica (Brassicaceae) as indicated by rDNA sequences. Systematic Botany 21: 559-566
Petersen G. O. Seberg S. Larsen 2001 The phylogenetic and taxonomic position of Lilaeopsis (Apiaceae), with notes on the applicability of ITS sequence data for phylogenetic reconstruction. Australian Systematic Botany 15: 181-191[CrossRef][Web of Science]
Pimenov M. G. M. V. Leonov 1993 The genera of the Umbelliferae. Royal Botanic Gardens, Kew, UK
Pimenov M. M. G. Vasil'eva M. V. Leonov J. V. Daushkevich 2003 Karyotaxonomical analysis in the Umbelliferae. Science Publishers, Enfield, New Hampshire, USA
Posada D. K. A. Crandall 1998 Modeltest: testing the model of DNA substitution. Bioinformatics 14: 817-818
Rambaut A. 1996 Se-Al: sequence alignment editor, version 1.dl. Department of Zoology, University of Oxford, Oxford, UK. Available at website, http://evolve.zoo.ox.ac.uk/
Raven P. H. 1973 Evolution of subalpine and alpine plant groups in New Zealand. New Zealand Journal of Botany 11: 177-200
Raven P. H. 1975 The bases of angiosperm phylogeny: cytology. Annals of the Missouri Botanical Garden 62: 724-764[CrossRef][Web of Science]
Ronquist F. 1997 Dispersal-vicariance analysis: a new approach to the quantification of historical biogeography. Systematic Biology 46: 195-203[CrossRef][Web of Science]
Ronquist F. J. P. Huelsenbeck 2003 MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19: 1572-1574
Salisbury E. 1970 The pioneer vegetation of exposed muds and its biological features. Philosophical Transactions of the Royal Society of London, B, Biological Sciences 259: 207-255[CrossRef]
Sanderson M. J. M. J. Donoghue 1989 Patterns of variation in levels of homoplasy. Evolution 43: 1781-1795[CrossRef][Web of Science]
Sang T. D. J. Crawford S.-C. Kim T. F. Stuessy 1994 Radiation of the endemic genus Dendroseris (Asteraceae) on the Juan Fernandez Islands: evidence from sequences of the ITS regions of nuclear ribosomal DNA. American Journal of Botany 81: 1494-1501[CrossRef][Web of Science]
Sanmartín I. F. Ronquist 2004 Southern hemisphere biogeography inferred by event-based models: plant versus animal patterns. Systematic Biology 53: 216-243[CrossRef][Web of Science][Medline]
Schaal B. A. W. J. Leverich 2001 Plant population biology and systematics. Taxon 50: 679-695[CrossRef][Web of Science]
Schaal B. A. D. A. Hayworth K. M. Olsen J. T. Rauscher W. A. Smith 1998 Phylogeographic studies in plants: problems and prospects. Molecular Ecology: 7: 465-474[CrossRef]
Shan R. H. M. L. Sheh [eds.] 1979 Umbelliferae. Flora Republicae Popularis Sinicae, vol. 55, part 1. Science Press, Beijing, China (in Chinese)
Shimizu K. K. J. M. Cork A. L. Caicedo C. A. Mays R. C. Moore K. M. Olsen S. Ruzsa G. Coop C. D. Bustamante P. Awadalla M. D. Purugganan 2004 Darwinian selection on a selfing locus. Science 306: 2081-2084
Shneyer V. S. G. P. Borschtschenko M. G. Pimenov 1992 The tribe Smyrnieae (Umbelliferae) in the light of serotaxonomical analysis. Plant Systematics and Evolution 182: 135-148[CrossRef][Web of Science]
Simmons M. P. H. Ochoterena 2000 Gaps as characters in sequence-based phylogenetic analyses. Systematic Biology 49: 369-381[CrossRef][Web of Science][Medline]
Smith J. M. B. 1986 Origins of Australasian tropicalpine and alpine floras. In B. Barlow [ed.], Flora and fauna of alpine Australasia: ages and origins, 109128. CSIRO, Melbourne, Australia
Sorenson M. D. 1999 TreeRot, version 2. Boston University, Boston, Massachusetts, USA. Available at website, http://people.bu.edu/msoren/TreeRot.html/
Spalik K. S. R. Downie 2001 The utility of morphological characters for inferring phylogeny in Scandiceae subtribe Scandicinae (Apiaceae). Annals of the Missouri Botanical Garden 88: 270-301[CrossRef][Web of Science]
Spalik K. A. Wojewódzka S. R. Downie 2001a Delimitation of genera in Apiaceae with examples from Scandiceae subtribe Scandicinae. Edinburgh Journal of Botany 58: 331-346[CrossRef]
Spalik K. A. Wojewódzka S. R. Downie 2001b The evolution of fruit in Scandiceae subtribe Scandicinae (Apiaceae). Canadian Journal of Botany 79: 1358-1374[CrossRef][Web of Science]
Spalik K. J.-P. Reduron S. R. Downie 2004 The phylogenetic position of Peucedanum sensu lato and allied genera and their placement in tribe Selineae (Apiaceae, subfamily Apioideae). Plant Systematics and Evolution 243: 189-210[CrossRef][Web of Science]
Stevens P. F. 1990 Nomenclatural stability, taxonomic instinct, and flora writinga recipe for disaster?. In P. Baas, K. Kalkman, and R. Geesink [eds.], The plant diversity of Malesia, 387410. Kluwer, Dordrecht, Netherlands
Sun F.-J. S. R. Downie R. L. Hartman 2004 An ITS-based phylogenetic analysis of the perennial endemic Apiaceae subfamily Apioideae of western North America. Systematic Botany 29: 419-431[CrossRef][Web of Science]
Swindell S. R. T. N. Plasterer 1997 SEQMAN. Contig assembly. In S. R. Swindell [ed.] Sequence data analysis guidebook, 7589. Human Press, Totowa, New Jersey, USA
Swofford D. L. 2003 PAUP*: phylogenetic analysis using parsimony (*and other methods). Sinauer, Sunderland, Massachusetts, USA
Sytsma K. J. B. A. Schaal 1985 Phylogenetics of the Lisianthius skinneri (Gentianaceae) species complex in Panama utilizing DNA restriction fragment analysis. Evolution 39: 594-608[CrossRef][Web of Science]
Templeton A. R. K. A. Crandall C. F. Sing 1992 A cladistic analysis of phenotypic associations with haplotypes inferring from restriction endonuclease mapping and DNA sequence data. III. Cladogram estimation. Genetics 132: 619-633[Abstract]
Templeton A. R. S. D. Maskas M. B. Cruzan 2000 Gene trees: a powerful tool for exploring the evolutionary biology of species and speciation. Plant Species Biology 15: 211-222
Thorne R. F. 1972 Major disjunctions in the geographic ranges of seed plants. Quarterly Review of Biology 47: 365-411[CrossRef]
Tiffney B. H. S. R. Manchester 2001 The use of geological and paleontological evidence in evaluating plant phylogeographic hypotheses in the Northern Hemisphere Tertiary. International Journal of Plant Sciences 162: (Supplement) S3-S17[CrossRef]
Tribsch A. T. F. Stuessy 2003 Evolution and phylogeography of arctic and alpine plants in Europe: introduction. Taxon 52: 415-416[CrossRef][Web of Science]
Urban I. 18801882 Zur flora Südamerikas, besonders Brasiliens. Linnaea 43: 859-931
van Royen P. 1983 The alpine flora of New Guinea, vol. 4. J. Cramer, Vaduz, Liechtenstein
van Steenis C. G. G. J. 1962 The land-bridge theory in botany with particular reference to tropical plants. Blumea 11: 235-372
van Steenis C. G. G. J. 1963 Pacific plant areas, vol. 1. Bureau of Printing, Manila, Philippines
van Steenis C. G. G. J. 1964 Plant geography of the mountain flora of Mt. Kinabalu. Proceedings of the Royal Society of London, B, Biological Sciences 161: 7-38
Walsh N. G. T. J. Entwisle 1999 Flora of Victoria, vol. 4. Inkata Press, Melbourne, Australia
Wen J. P. P. Lowry II J. L. Walck K.-O. Yoo 2002 Phylogenetic and biogeographic diversification in Osmorhiza (Apiaceae). Annals of the Missouri Botanical Garden 89: 414-428[CrossRef][Web of Science]
White T. J. T. Bruns S. Lee J. Taylor 1990 Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. In M. A. Innis, D. H. Gelfand, J. J. Sninsky, and T. J. White [eds.], PCR protocols: a guide to methods and applications, 315322. Academic Press, San Diego, California, USA
Winkworth R. C. S. J. Wagstaff D. Glenny P. J. Lockhart 2002 Plant dispersal N.E.W.S from New Zealand. Trends in Ecology and Evolution 17: 514-520[CrossRef]
Wu C.-I. W.-H. Li 1985 Evidence for higher rates of nucleotide substitution in rodents than in man. Proceedings of the National Academy of Sciences, USA 82: 1741-1745
Wu C.-Y. 1998 Delineation and unique features of the Sino-Japanese floristic region. In D. E. Boufford and H. Ohba [eds.], Sino-Japanese flora: its characteristics and diversification. University Museum, University of Tokyo, Bulletin Number 37, 135. Tokyo University Press, Tokyo, Japan
Xia X. Z. Xie 2001 DAMBE: data analysis in molecular biology and evolution. Journal of Heredity 92: 371-373
Zhang L.-B. H. P. Comes J. W. Kadereit 2001 Phylogeny and Quaternary history of the European montane/alpine endemic Soldanella (Primulaceae) based on ITS and AFLP variation. American Journal of Botany 88: 2331-2345
![]()
CiteULike
Complore
Connotea
Del.icio.us
Digg
Facebook
Reddit
Technorati
Twitter What's this?
This article has been cited by other articles:
![]() |
F. Vandelook, N. Bolle, and J. A. Van Assche Seed Dormancy and Germination of the European Chaerophyllum temulum (Apiaceae), a Member of a Trans-Atlantic Genus Ann. Bot., August 1, 2007; 100(2): 233 - 239. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. I. Calvino, P. M. Tilney, B.-E. v. Wyk, and S. R. Downie A molecular phylogenetic study of southern African Apiaceae Am. J. Botany, December 1, 2006; 93(12): 1828 - 1847. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |