Am. J. Bot. Download to a Citation Manager
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


  Add to CiteULike   Add to Complore   Add to Connotea   Add to Del.icio.us   Add to Digg   Add to Facebook   Add to Reddit   Add to Technorati   Add to Twitter
What's this?
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplemental Data
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Van den heede, C. J.
Right arrow Articles by Chase, M. W.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Van den heede, C. J.
Right arrow Articles by Chase, M. W.
Agricola
Right arrow Articles by Van den heede, C. J.
Right arrow Articles by Chase, M. W.
Social Bookmarking
 Add to CiteULike   Add to Complore   Add to Connotea   Add to Del.icio.us   Add to Digg   Add to Facebook   Add to Reddit   Add to Technorati   Add to Twitter  
What's this?
(American Journal of Botany. 2003;90:481-495.)
© 2003 Botanical Society of America, Inc.


Systematics

Phylogenetic analysis of Asplenium subgenus Ceterach (Pteridophyta: Aspleniaceae) based on plastid and nuclear ribosomal ITS DNA sequences1

Caroline J. Van den heede2, Ronald L. L. Viane2,4 and Mark W. Chase3

2Pteridological Section, Department of Biology, Ghent University, K. L. Ledeganckstraat 35, B-9000 Ghent, Belgium; 3Jodrell Laboratory, Royal Botanic Gardens, Richmond, Surrey, TW9 3DS, UK

Received for publication May 23, 2002. Accepted for publication October 31, 2002.


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Phylogenetic relationships among 20 taxa of the fern genus Asplenium subgenus Ceterach (Filicopsida, represented by 73 accessions) were investigated using DNA sequence data from the nuclear ribosomal internal transcribed spacers (ITS nDNA) and plastid trnL-F intergenic spacer. In addition, a single sample per taxon was used in an analysis of the plastid rbcL gene. Chromosome counts were determined for all the samples, and these demonstrated a range from diploid to octoploid. Analyses of the DNA sequence data indicated that Asplenium subgenus Ceterach is polyphyletic, implicating homoplasy in the characters previously used to circumscribe this taxon. Plastid trnL-F and rbcL analyses resulted in identical tree topologies. The trees produced from the separate plastid and nuclear matrices agree in (1) the recognition of identical groups of accessions corresponding to A. dalhousiae, A. ceterach, A. aureum, A. cordatum, A. phillipsianum, and A. haughtonii; (2) the division of A. subg. Ceterach into two subclades, a Eurasian-Macaronesian and a strictly African alliance; (3) the position of A. dalhousiae as a member of the former subclade; (4) the lack of genetic variation in A. cordatum despite its morphological variability; and (5) the clustering of each autopolyploid with their diploid ancestor. However, the plastid and nuclear trees differ in their placement of A. haughtonii and A. dalhousiae, which might be due to different evolutionary histories of nuclear and plastid genomes, and is possibly an indication of ancient hybridization. The analyses confirm the existence of several strictly African taxa. Asplenium phillipsianum and A. cordatum each form species complexes of diploid and autopolyploid taxa, from which a third, morphologically intermediate, allotetraploid species has originated. Asplenium haughtonii is a distinct endemic species from Saint Helena. The maternally inherited plastid sequences support the hypothesis that A. aureum is an ancestor of A. lolegnamense and of A. octoploideum. Because gene conversion did not eliminate divergent ITS alleles in the allopolyploids, their reticulate ancestry could be demonstrated. Biparentally inherited nrITS sequences support the allopolyploid status of A. aureum, A. lolegnamense, and A. punjabense, indicating they share the ancestral A. javorkeanum genome.

Key Words: Africa • Aspleniaceae • Asplenium subgenus Ceterach • gene conversion • ITS • molecular phylogenetics • Pteridophyta • rbcLtrnL-trnF intergenic spacer


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Asplenium subgenus Ceterach (Willd.) Bir et al. is a small group of about 10 fern taxa within the large, subcosmopolitan genus Asplenium L. (720 species; Kramer and Viane, 1990 ). Although characterized by a variable gross morphology, Aspleniaceae are anatomically and cytologically homogeneous. The basic chromosome number (x) is 36. Both auto- and allopolyploidy are common driving forces of evolution in Aspleniaceae (Wagner, 1954 ; Lovis, 1977 ; Reichstein, 1981 ). Thus, in the well-known Holarctic fern flora, Asplenium consists of about 53% (ancestral) diploids, 24% autopolyploids, 20% allopolyploids, and 3% apogamous taxa. Polyploidy is at least equally common in the tropical taxa, but their status and ancestry remain largely unexplored.

There is no modern monograph of Aspleniaceae including a phylogenetic assessment, and the monophyly of A. subgenus Ceterach was never investigated. The only worldwide overview of Asplenium is that of Mettenius (1859) . Cytological studies have led to the present well-known phylogeny of northern temperate taxa (Viane et al., 1993 ; Wagner et al., 1993 ), but relationships among tropical and southern temperate taxa remain largely unknown. Chromosome numbers were compiled by Löve et al. (1977) and Reichstein (1981) .

Asplenium subgenus Ceterach contains xerophytic rock ferns with the dorsal side of the lamina densely covered with scales (Fig. 1). The indusium is reduced or absent, and the veins are anastomosing (Eurasian group) or free (strictly African group). Ploidy ranges from diploid to octoploid. The number of species is still unclear and varies from five or seven to 14, depending on the weight given to certain morphological characters and whether different cytotypes are considered separate species (Van den heede et al., 2002 ). Most taxa are morphologically similar and difficult to distinguish in the field. Representatives of Ceterach occur in Europe, Macaronesia, throughout southwestern Asia to the Himalayas, and in the drier parts of North, East, and southern Africa. Islands like Madeira and Saint Helena have endemic species.



View larger version (51K):
[in this window]
[in a new window]
 
Fig. 1. Pinnae with venation pattern and sori (A, C), and dorsal laminal scales (B, D). Asplenium ceterach (A, B): left pinna with scales removed to show veins with marginal anastomoses and almost naked sori; scales with penta- or hexagonal cells. Asplenium dalhousiae (C, D): pinna with free veins and sori with distinct dorsal indusium; scales with elongated cells. Scale bar (A, C) = 1 cm, (B, D) = 0.5 cm. (A) from Luerssen, 1889

 
Rasbach et al. (1987) have shown that A. ceterach L. is an autotetraploid, which must have originated via chromosome doubling in the morphologically identical diploid, A. javorkeanum Vida (synonym: A. ceterach subsp. bivalens (D.E.Mey.) Greuter and Burdet). Asplenium xmantoniae Váróczy and Vida is the triploid hybrid between A. javorkeanum and A. ceterach. Within the same complex, the autohexaploid A. cyprium Viane and Van den heede described from Cyprus is also known from Greece and Sicily (Viane et al., 1996 ; Van den heede et al., 2002 ).

To date, three cytologically different endemic species are known from Macaronesia. Asplenium aureum Cav. was found to be tetraploid by Manton (1950) and suggested to be allotetraploid by Vida and Reichstein (Viane and Reichstein, 1992 ). Asplenium octoploideum Viane and Van den heede (synonym: A. parvifolium sensu Vida and Reichstein) is the allo-octoploid species that probably formed by chromosome doubling in the tetraploid hybrid between A. aureum and A. ceterach (Viane and Reichstein, 1992 ; Van den heede and Viane, 2002 ). Manton et al. (1986) and Gibby and Lovis (1989) suggested that the Madeiran allohexaploid A. lolegnamense (Gibby and Lovis) Viane originated by hybridization between A. aureum and A. octoploideum. Though at present their ancestors only occur together in this archipelago, there seems to be no conclusive molecular evidence that A. octoploideum and A. lolegnamense were formed on the Canary Islands, thus they may show a relict distribution pattern. Asplenium xchasmophilum Van den heede and Viane is the natural, sterile hexaploid hybrid between A. ceterach and A. octoploideum, at present only known from its type locality on Tenerife (Van den heede and Viane, 2002 ). Chromosome pairing behavior in this hybrid confirms the allotetraploid status of A. aureum (C. J. Van den heede and R. L. L. Viane, unpublished data).

For the African A. capense-cordatum-phillipsianum complex, some chromosome counts were published by Bir et al. (1985) , who reported one diploid and three tetraploid specimens without mentioning auto- or allopolyploidy. We were able to confirm (Bir et al., 1985 ) that Asplenium haughtonii (Hook.) Bir et al. from Saint Helena is tetraploid.

Within Aspleniaceae, A. dalhousiae Hook., sometimes put into subg. Ceterachopsis (J.Sm. ex Ching) Bir on the sole basis of its frond morphology, represents the putatively closest relative of A. subg. Ceterach. Copeland (1947) already pointed out that it is impossible to define Asplenium so as to exclude the Ceterachopsis group. According to Pichi Sermolli (1979 , p. 191), "There are good grounds for thinking that Ceterach derived from Asplenium through the intermediate genus Ceterachopsis." Four characters distinguish A. dalhousiae from the European Ceterach group: absence of a dense layer of dorsal scales, structure of the laminal paleae, open venation pattern, and presence of a true, well-developed indusium (Fig. 1). In fact, these are all putative plesiomorphies within Asplenium, from which Ceterachopsis only differs by the single apomorphy: a pinnatipartite to pinnatisect (De Candolle, 1813 ; Pichi Sermolli, 1972 ) lamina. In most enzyme systems studied, A. dalhousiae shows unique zymograms different from those of the members of A. subg. Ceterach (C. J. Van den heede, unpublished data). Asplenium dalhousiae has a disjunct range with a wide distribution in southwestern Asia and smaller ones in Ethiopia and in southwestern North America. We investigated one specimen from each of these areas of distribution.

We also included Asplenium punjabense Bir et al., endemic to the western Himalayas, in the analysis. This fertile allohexaploid has been previously documented to have originated by chromosome doubling of a triploid hybrid between diploid A. dalhousiae and tetraploid A. ceterach (Bir et al., 1985 ) and can be used to evaluate which molecular method traces reticulate evolution best.

While preparing a biosystematic revision of A. subg. Ceterach (C. J. Van den heede, unpublished data), we generated hypotheses about possible relationships between members (and hybrids) of this complex based on (micro)morphological and cytological studies. We felt the need to obtain independent estimates for our phylogenetic ideas with molecular data, such as isozymes and DNA. During our fern studies for the Flora of Tropical East Africa, it also became evident that further research was necessary to find out how many Ceterach species exist in Africa. Schelpe (1970) and Schelpe and Anthony (1986) recognized only a single, variable species, Ceterach cordatum Sw., because they could not find characters to distinguish additional taxa.

The first DNA sequencing of pteridophytes was done by Hasebe et al. (1993) , who analyzed nucleotide variation in the plastid rbcL gene from several leptosporangiate fern families. At present, most molecular analyses of Pteridophyta are based on rbcL sequences (e.g., Hasebe et al., 1994 , 1995 ; Wolf et al., 1994 , 1999 ; Murakami, 1995 ; Pryer et al., 1995 ; Murakami et al., 1998 , 1999a , b ; Yatabe et al., 1999 ; Schulze et al., 2001 ). Besides rbcL, some fern studies used sequences of other regions of plastid DNA, such as 16S rDNA (Manhart, 1995 ), internal transcribed spacers (plastid ITS; Pahnke et al., 1996 ) and trnL-F (Hauk et al., 1996 ; Haufler et al., 2000 ). The first nuclear gene used for phylogenetic research in pteridophytes was 18S rRNA (Raubeson and Stein, 1995 ; Wolf, 1995 , 1996 ; Kranz and Huss, 1996 ), which is useful for examining relationships among the major groups of Pteridophyta, although the resulting trees are similar but less resolved than those for rbcL. Wolf (1996) also tried to sequence two other nuclear regions of fern genomes, 26S and the internal transcribed spacers (ITS), but concluded that it was difficult or impossible to get "clean" sequences and that fine tuning of techniques and cloning might be necessary. Additionally, the few sequences he obtained were difficult to align because of the apparent rapid accumulation of insertions and deletions. Wolf (1996) suggested that nucleotide variation in nrITS could be taxonomically useful within fern genera, species, and populations. The first analysis of fern ITS regions was on cheilanthoid ferns by Gastony and Rollo (1998) , who did not report any difficulties similar to those encountered by Wolf (1996) .

In this study, three DNA sequence regions (plastid trnL-trnF intergenic spacer, the rbcL gene, and nrITS) are used to examine interspecific relationships within A. subgenera Ceterach and Ceterachopsis. We also compare sequences of specimens from different populations of the same species for trnL-F and ITS. Intraspecific variation in sequences of plastid DNA (reviewed by Harris and Ingram, 1991 ; Soltis et al., 1992 ; Fujii et al., 1997 ) and ITS of nrDNA have been reported (Baldwin et al., 1995 ; Fujii et al., 2001 ). According to Harris and Ingram (1991) , intra-individual sampling should be considered for studies at the species level. Plastid DNA polymorphism is quite common in perennial taxa, and even in the absence of morphological variation they may show variation due to substitutions and insertions/deletions (indels). In wild taxa, plastid DNA variation may become established by hybridization and "introgression" and thus increase genetic diversity (Harris and Ingram, 1991 ).

The trnL-F region is the non-coding, intergenic spacer between the trnL (UAA)3' exon and trnF (GAA) of plastid DNA, for which Taberlet et al. (1991) designed "universal" primers. Taberlet et al. (1991) suggested that sequences of intergenic spacers of plastid DNA might be useful for evolutionary studies of closely related species and populations of the same species. In Iridaceae, the trnL-F region evolves three times faster than rbcL (Soltis and Soltis, 1998 ). However, in other groups of plants, the trnL-F intergenic spacer is as, or more, conserved than rbcL (Asmussen and Chase, 2001 ). In Ophioglossaceae, the number of variable positions in trnL-F is 3–5 times higher than in rbcL despite the fact that the trnL-F spacer is about one-third the length of rbcL; Hauk et al. (1996) found that rbcL and trnL-F analyses produced nearly identical tree topologies.

Sequencing of the ITS region of 18S–26S nuclear ribosomal DNA has been useful for phylogenetic reconstruction at specific and generic levels in many angiosperms (reviewed in Baldwin et al., 1995 ). This region consists of two spacers, ITS1 and ITS2, separated by the conserved 5.8S gene. Plant ribosomal DNA consists of thousands of repeats and often paralogues (Soltis and Soltis, 1998 ). Due to concerted evolution, the sequence similarity between different, functional ribosomal copies within an organism is extremely high in most plant species (Dover, 1982 ; Arnheim, 1983 ; King, 1993 ; Baldwin et al., 1995 ). Concerted evolution is not well understood, but includes several mechanisms such as unequal crossing over, gene conversion, and transposition (Dover, 1982 ; Soltis and Soltis, 1998 ; Graur and Li, 2000 ). However, concerted evolution has not homogenized rDNA repeats in all taxa. Some plant genomes exhibit ribosomal DNA diversity including multiple functional ITS copies (divergent paralogues), nonfunctional ITS copies ("pseudogenes"), and recombinants (Buckler et al., 1997 ; Soltis and Soltis, 1998 ; Hershkovitz et al., 1999 ). According to Hershkovitz et al. (1999) , the 5.8S gene provides an indicator of the functionality of cloned ITS copies. An aligned 5.8S sequence showing transversions and/or indels that are absent in the other sampled taxa might indicate that the particular rDNA copy is not expressed. Within-individual ITS polymorphism may result from mutation or hybridization between genomes with different, functional ITS sequences. Sang et al. (1995) were able to document reticulate evolution in Paeonia and detect hybridization using nrITS. Both parental ITS sequences have been maintained in allopolyploid species of Krigia and Tragopogon (Asteraceae), as well as in Paeonia species of about one million years of age (Kim and Jansen, 1994 ; Sang et al., 1995 ; Soltis and Soltis, 1998 ). Similarly, Campbell et al. (1997) uncovered extensive ITS sequence polymorphism within Amelanchier (Rosaceae) taxa, which are suspected to be relatively old. They postulated that this polymorphism is mainly created by hybridization and that polyploidy and/or agamospermy could retard concerted evolution.

Comparisons of plastid and nuclear DNA analyses can significantly improve understanding of the origin of polyploid species (Soltis and Soltis, 1991 ; Soltis et al., 1992 ; Wendel and Doyle, 1998 ). Because several of the ferns in this study were thought to be autopolyploid, it was necessary to know the exact ploidy of each individual to interpret the molecular results. Therefore we determined the cytological status of each accession used in this study. The main objectives of this study were to investigate the interspecific relationships within A. subg. Ceterach, evaluate earlier hypotheses about the phylogeny in this group, look for evidence of reticulate evolution, and determine how many Ceterach species occur in Africa.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Taxon sampling
Apart from MO04926887, collected by A. Reina et al. (MO), vouchers are deposited at GENT (taxa, localities, vouchers, ploidies, genome formulae, and GenBank accession numbers have been archived at the Botanical Society of America website [http://ajbsupp.botany.org/v90/]). Localities of DNA vouchers are shown on the map of Fig. 2. Material for this study has been cultivated in Ghent Botanical Garden (Belgium) from spores collected by C. J. Van den heede and R. L. L. Viane in the wild or from herbarium specimens. From the living Ceterach(opsis) collection of about 1100 specimens, 73 accessions previously examined in chromosome and isozyme studies were selected for our DNA work (http://ajbsupp.botany.org/v90/). Because plant health influences isozyme results (Stebbins, 1989 ; Acquaah, 1992 ; Van den heede et al., 2002 ), we included only vigorously growing plants in this study.



View larger version (28K):
[in this window]
[in a new window]
 
Fig. 2. Localities of vouchers used for DNA sequencing (see http://ajbsupp.botany.org/v90/).

 
Cytological preparation
For chromosome counts, immature spore mother cells were fixed in freshly prepared 3 : 1 absolute ethanol : glacial acetic acid and stored frozen until required. Acetocarmine squash preparations were manually made as described by Heitz (1925 , 1950 ) and Manton (1950) . Photographs were taken with an Olympus BH2 microscope. Preparations were made permanent by dehydrating the cover slip and slide in graded mixtures of acetic acid and absolute ethanol, followed by mounting in Euparal (Manton, 1950 ). All permanent preparations are retained in the Pteridological Section of the Department of Biology at Ghent University, Ghent, East Flanders, Belgium. Counted plants were used as standards to compare the nuclear DNA content of additional specimens by a flow cytometer (Partec PA-1), using the manufacturer's protocol (Partec, Münster, Nordrein-Westfalen, Germany).

DNA extraction
Total DNA was extracted from fresh leaf material (0.5–1.0 g) using the modified 2 x cetyltrimethyl ammonium bromide (CTAB) protocol of Doyle and Doyle (1987) . The DNA was precipitated in cold (4°C) absolute ethanol for at least 5 d, pelleted, and purified via cesium chloride/ethidium bromide gradient (1.55 g/mL). Subsequently, CsCl and other residual salts were removed by dialysis. Ethidium bromide was removed with butanol.

Amplification and cloning
The trnL-trnF intergenic spacer was amplified using the universal primers e and f of Taberlet et al. (1991) . Polymerase chain reaction (PCR) amplification of rbcL fragments followed Hasebe et al. (1994) using the same forward and reverse primers. The thermal cycling protocol for plastid DNA consisted of 28–30 cycles, each comprising 1 min denaturation at 94°C, 30 s annealing at 48°C, and 1 min chain extension at 72°C. The final extension was at 72°C for 7 min.

The ITS region (ITS1, the 5.8 gene, and ITS2) was amplified with a primer designed to exclude a frequent bryophyte contaminant (epiphyte) (5'-CCTGCGGAAGGATACTTGTCG-3'). We used the reverse primer, ITS4, described by White et al. (1990) . Thermal cycling of ITS involved an initial premelt at 94°C for 2 min followed by 30 cycles, each consisting of 1 min denaturation at 94°C, 1 min annealing at 50°C, and 2.5 min chain extension at 72°C. The final stage was a 7-min extension at 72°C. Addition of 2% dimethylsulfoxide (DMSO) greatly improved both PCR and cycle sequencing reactions of the guanine and cylosine (GC)-rich ITS region (Varadaraj and Skinner, 1994 ); in PCR its addition helped to guarantee that the amplified product contained mostly functional copies of ITS, whereas in cycle sequencing DMSO helped to overcome premature terminations caused by GC-rich regions.

The ITS PCR products of samples with multiple sequence signals were cloned using Promega's pGEM-T Easy Vector System (Crawley, West Sussex, UK). Ligations and transformations were prepared according to the technical manual provided by the manufacturer. Transformed E. coli cells were spread onto agar plates and incubated at 37°C for 17 h. Between 9 and 12 white colonies per plate (selected at random) were touched with a sterile 1 µL tip and subsequently used as template for reamplification with the same primers. Double-stranded PCR products were purified with Concert Rapid PCR Purification System using the manufacturer's protocols (Invitrogen Life Technologies, Paisley, Scotland, UK).

Sequencing and alignment
Cycle sequencing reactions of 10 µL were carried out using the ABI Taq Dye-Deoxy Terminator Cycle Sequencing Kit (Applied Biosystems, Warrington, Cheshire, UK) with 1 µL of the dye mixture. The PCR amplification primers were used also as sequencing primers. For trnL-F/rbcL, both strands were sequenced on an ABI 377 automated sequencer (according to the manufacturer's protocols). Sequencing was usually straightforward with the exception of A. phillipsianum (Kümmerle) Bir et al. specimens due to a GC-rich area near the middle of the trnL-F intergenic spacer. All ITS clones from PCR products with heterogeneity were first sequenced with the forward primer. From identical sequences, we chose one clone for sequencing the reverse strand. Sequences were edited and assembled using Sequencher software (Gene Codes Corporation, Ann Arbor, Michigan, USA). All ITS sequence versions were submitted to BLAST search in GenBank to exclude contaminants (e.g., fungi; Soltis and Soltis, 1998 ), which resulted in complete nrITS sequences for 42 taxa. Putative allopolyploids showed two sequence types. Assembled sequences were manually aligned prior to analysis. No gaps were needed to align any of the rbcL and 5.8S sequences. Because ITS sequences for other ferns were not available in the GenBank, boundaries of ITS1 and ITS2 were determined by comparison with sequences of Asteraceae (Baldwin, 1992 ) and of Nyssa ogeche (Nyssaceae; Fan and Xiang, 2001 ). Complete sequences for the trnL-F intergenic spacer were obtained for 72 samples. For rbcL, we randomly selected a single accession from each species and produced 11 new rbcL sequences (http://ajbsupp.botany.org/v90/).

Data analysis
Based on rbcL evidence (Hasebe et al., 1995 ; Pryer et al., 1995 ), trnL-F sequences of the following leptosporangiate ferns were used as outgroups: a putatively more primitive species of Dennstaedtia (R. Cranfill, University of California, Berkeley, unpublished data); two Polypodium species (AF159183–159184) published by Haufler et al. (2000) ; and A. nidus L. (AF425118), A. scolopendrium L., and A. unilaterale Lam. sequences provided by R. Cranfill (unpublished data). For rbcL, sequences of 27 Asplenium species taken from GenBank (http://ajbsupp.botany.org/v90/) were included. Dennstaedtia samoensis (Brack.) T. Moore (U18637) and Microlepia strigosa (Thunb.) C. Presl (U05931) from the more distantly related Dennstaedtiaceae were designated as the outgroup for the rbcL analysis.

Because A. dalhousiae came out as a member of the ingroup in our trnL-F/rbcL trees and due to the fact that no ITS sequences of other ferns were available in GenBank, we used midpoint rooting for the ITS tree. This procedure gave a similar pattern of relationships as the plastid trees, for which we could use outgroups. The matrices were analyzed with PAUP 4.0b10 (Swofford, 2002 ), using a heuristic search and maximum parsimony methods under the Fitch criterion (unordered and equally weighted characters; Fitch, 1971 ). Shortest trees were obtained with 100 replicates of random taxon addition and tree bisection-reconnection (TBR) branch swapping. MulTrees was on (saving multiple, equally parsimonious trees), and no tree limit was used. All other settings were the standard defaults. Gaps (hypothesized insertions/deletions [indels]) were treated as missing values and not separately scored; thus gap information was not included in the analyses (bootstrap percentages were already high without including them). Internal support was assessed using 1000 bootstrap replicates, saving five trees per replicate to reduce time spent in swapping large numbers of trees, and TBR swapping. Percentages of less than 50% are not reported here; 85–100% is considered as strong support. The trees of Figs. 3 and 5 are shown with branch lengths proportional to the number of estimated substitutions (ACCTRAN optimization). The sequences are deposited in GenBank (http://ajbsupp.botany.org/v90/). Alternative topologies among the equally most parsimonious trees identified affect only groupings of accessions within species, and such patterns are not issues in this study.



View larger version (54K):
[in this window]
[in a new window]
 
Fig. 3. One randomly selected tree of the 100 shortest trees of Asplenium subg. Ceterach and A. dalhousiae resulting from parsimony analysis of sequences of the trnL-F intergenic spacer; length = 446 steps, CI = 0.88, and RI = 0.97. Numbers above the branches are lengths (ACCTRAN optimization); those below the branches indicate bootstrap percentages greater than 50% (1000 replicates). Dennstaedtia was specified as outgroup. The arrow indicates the single node that collapses in the strict consensus of all shortest trees. Next to ploidy (2x, 3x, etc.), putative ancestral (anc), autoploid (auto), or alloploid (allo) status is indicated. The aureum, ceterach, and African clades are indicated by a vertical bar

 


View larger version (39K):
[in this window]
[in a new window]
 
Fig. 5. One randomly selected of the 96 most parsimonious, midpoint rooted trees of Asplenium subg. Ceterach and A. dalhousiae based on analysis of nrITS sequence data; length = 335 steps, CI = 0.85, and RI = 0.98. Fitch branch lengths (ACCTRAN optimization) are shown above and bootstrap percentages greater than 50% (1000 replicates) below the branches. The tree shown is identical to the strict consensus tree. Next to ploidy (2x, 3x, etc.), putative ancestral (anc), autoploid (auto), or alloploid (allo) status is indicated. The aureum, ceterach, and African clades are indicated by a vertical bar

 

    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Cytology
Chromosomes were counted for 25 of the 73 accessions used in this study. These samples were subsequently used as standards to determine by flow cytometry the ploidy of the remaining 48 specimens (http://ajbsupp.botany.org/v90/).

Analysis of trnL-trnF intergenic spacer data
The length of the trnL-trnF intergenic spacer varied from 338 base pairs (bp) in A. cordatum (CV814a) to 404 bp in A. dalhousiae (TR7634). The aligned matrix was 525 bp long with the outgroups and 450 bp without. The matrix contained 238 variable sites of which 148 were potentially parsimony informative, whereas the Ceterach matrix (without the outgroups) included only 76 potentially parsimony informative characters (80 variable positions in all). Thirteen gaps varying from 1 to 20 bp in length were introduced to align the Eurasian-Macaronesian Ceterach group with the African ceterach group. All accessions of Asplenium dalhousiae differ from those of the other taxa in the possession of a unique duplication of 19 bp. Furthermore, A. punjabense (TR7675) and A. ceterach from Italy (CV494) share another unique duplication of 17 bp. Sequence alignment of the complete matrix including the outgroups required 10 additional gaps of 1–20 bp.

Analysis resulted in 100 most parsimonious trees of 327 steps with a consistency index (CI) = 0.88 and a retention index (RI) = 0.97. One of the shortest trees (selected at random) is presented in Fig. 3 with the branch lengths (i.e., the numbers of substitutions; ACCTRAN optimization) above and the bootstrap percentages below each branch (groups that received less than 50% are not indicated). There is little variation among the trees due to the fact that PAUP calculates the number of trees based on arbitrary dichotomies of accessions that have no differences in their sequences. A strict consensus tree of all equally most parsimonious trees is consistent with the interspecific relationships of the ingroup species shown in this single tree. In all shortest trees, A. subg. Ceterach is polyphyletic. Within Aspleniaceae, A. unilaterale is sister to all other Asplenium species (bootstrap percentage [BP] 100), which in turn are composed of two well-supported subgroups. The first group (BP 83) corresponds to A. nidus, A. scolopendrium, and the Eurasian Ceterach group including A. dalhousiae; the second group (BP 100) consists of the strictly African ceterach species as currently recognized by C. J. Van den heede and R. L. L. Viane. Within the first subclade, the relationships of A. nidus, A. scolopendrium, and the Eurasian Ceterach group are unresolved in the strict consensus tree (see arrow; Fig. 3).

The African ceterach clade is divided into an A. cordatum group (BP 90) and an A. phillipsianum-haughtonii group (BP 93). The first group (BP 90) clusters all A. cordatum samples ranging from diploid (2x) to hexaploid (6x). The second group (BP 87) bifurcates into the A. phillipsianum (2x–6x; BP 94) and A. haughtonii (4x–6x; BP 99) clades. Intermediate forms (labeled as A. sp.) between A. cordatum and A. phillipsianum are found in either the A. cordatum or A. phillipsianum clade.

The well-supported (BP 97) non-African Ceterach clade includes two strongly supported subclades. The first (BP 100) consists of A. punjabense and all specimens belonging to the A. ceterach complex from the Mediterranean (Cyprus, Italy, Slovenia) and Tenerife. All species (2x, 4x, 6x) and hybrids (3x and 6x) have identical trnL-F sequences except for the morphologically identical tetraploid CV494 (Italy), which differs by three substitutions. In the second subclade (BP 100), A. dalhousiae is sister to the A. aureum clade. The A. dalhousiae clade forms a well-supported group (BP 100) in which the Ethiopian and Pakistani specimens differ from the Mexican specimen by a single substitution. The aureum clade is poorly supported (BP 64) and groups three Macaronesian species, A. aureum, A. lolegnamense, and A. octoploideum.

Analysis of rbcL data
To evaluate the putative polyphyly of A. subg. Ceterach shown by the trnL-F results, we assembled an rbcL matrix (sequences of 11 "ceterach" taxa and 27 other Asplenium species) using Dennstaedtia and Microlepia as outgroups. This rbcL matrix had 1301 characters (no indels present), of which 391 were variable and 260 potentially parsimony informative. A heuristic search performed under the same conditions as for the trnL-F analysis resulted in four most parsimonious trees of 788 steps with CI = 0.60 and RI = 0.75. One of the shortest trees (selected at random), Fig. 4, is identical to that produced with just the trnL-F data.



View larger version (25K):
[in this window]
[in a new window]
 
Fig. 4. One of four most parsimonious trees for Aspleniaceae based on our 11 rbcL sequences of taxa of Asplenium subg. Ceterach (bold) and 27 rbcL sequences of other Asplenium species presently available in GenBank (http://ajbsupp.botany.org/v90/); length = 788 steps, CI = 0.60, and RI = 0.75. Dennstaedtia and Microlepia were specified as outgroups. Numbers are bootstrap percentages greater than 50% (1000 replicates). Arrows indicate nodes that collapse in the strict consensus of all shortest trees. The polyphyly of A. subg. Ceterach (branches with thick lines) is confirmed by rbcL data

 
Within Aspleniaceae, the Hymenasplenium group (Mitui et al., 1989 ; Kramer and Viane, 1990 ; Murakami and Moran, 1993 ; Murakami, 1995 ; Cheng and Murakami, 1998 ) is sister to the other Asplenium species (BP 99), which agrees with the results of Murakami (1995) . Due to limited taxon sampling and low sequence divergence, several BP are low, and relationships of A. marinum, A. hemionitis, and the strictly African ceterach group are unresolved in the strict consensus tree. However, the polyphyly of A. subg. Ceterach is clearly supported by the rbcL data. Increased taxon sampling is needed to identify the sister group of the well-supported clade (BP 100) of strictly African ceterach taxa. Sister to the Eurasian Ceterach group (BP 99) is the strongly supported clade (BP 100) of Mediterranean A. sagittatum (DC.) Bange and Holarctic A. scolopendrium. Again, the non-African Ceterach group includes two fully supported subclades (BP 100). In the first, A. dalhousiae is sister to the strongly supported (BP 95) A. aureum clade clustering the three Macaronesian species. The second subclade groups all taxa of the A. ceterach complex, A. hybridum (Milde) Bange, and A. punjabense. Asplenium hybridum, endemic to the Dalmatian Islands (Croatia), is a fertile allotetraploid that originated by chromosome doubling in a diploid hybrid between diploid A. javorkeanum and diploid A. sagittatum (Vida, 1963 ; Reichstein, 1981 ). The position of A. hybridum indicates that A. javorkeanum served as the maternal parent in this event.

Analysis of ITS data
Internal transcribed spacer 1 ranged in length from 237 to 279 bp; ITS2 varied from 191 to 225 bp. The shortest ITS1 and ITS2 regions were found in A. dalhousiae. The GC content in ITS1 ranged from 47.7% in A. haughtonii (CV915) to 61.0% in A. aureum (Cav.) Buch (CV712) and for ITS2 from 50.5% in A. haughtonii (CV924) to 68.1% in A. javorkeanum (CV14). With its 164 bp, the 5.8S subunit was uniform in size although slightly variable in GC content (53.0–57.3%). Internal transcribed spacer 2 includes more (127) potentially parsimony informative characters than ITS1 (110). The ten variable positions of the 5.8S rDNA gene are all potentially parsimony-informative. In addition, the flanking portion of the 26S rDNA gene (48 bp) contains six informative sites, whereas the flanking part (only 2 bp) of the 18S gene is not informative. Sequence alignment required 23 and 15 gaps (varying from 1 to 40 bp in length) in ITS1 and ITS2, respectively. Indels were absent in the 5.8S rRNA gene.

The ITS1 + 5.8 S + ITS2 region had an aligned length of 878 bp. However, the total length of the region analyzed, including the short regions of 18S and 26S, consisted of 927 sites. Analysis resulted in 96 most parsimonious trees of 335 steps with CI = 0.85 and RI = 0.98. One of the shortest, midpoint-rooted trees (selected at random) is shown in Fig. 5 with branch lengths (ACCTRAN optimization) indicated above the branches and bootstrap percentages below. If ITS clones from a single plant were positioned in two different clades, we mark them as (copy type) I and II. Although PAUP calculated that there were 96 equally parsimonious shortest trees, the strict consensus tree is identical to the single tree illustrated.

Midpoint rooting of the ITS tree again divides subgenus Ceterach in two main groups: the Eurasian-Macaronesian species plus A. dalhousiae and the strictly African taxa. The former, maximally supported (BP 100) clade consists of the Eurasian and Macaronesian Ceterach species and bifurcates into two strongly supported subclades: the aureum-genome clade (BP 100) and the ceterach-genome clade (BP 90). The aureum clade consists of the ITS copy type I of both A. aureum (CV712) and (CV174). The ceterach clade splits into two highly supported groups: one (BP 89) with A. xmantoniae and two diploid A. javorkeanum specimens from Italy and Slovenia and the other (BP 100) with tetraploid A. ceterach from Tenerife (CV187), copy type II of A. punjabense (TR7634), A. aureum (CV712), and A. lolegnamense (CV174).

The A. dalhousiae clade (BP 100) is most closely related to the Eurasian-Macaronesian clade (this is of course dependent on the topology created by midpoint rooting) and includes sequences from three taxa. Copy type I of the Pakistani allohexaploid A. punjabense (TR7675) has an identical sequence to that of the diploid A. dalhousiae (TR7634) from Pakistan, whereas the third sequence of the Mexican plant differs by two substitutions. Copy type II of the allohexaploid A. punjabense (TR7675) clusters with the ceterach clade and differs by only three substitutions from A. ceterach (autotetraploid) from Tenerife.

The African group is divided in two subgroups (each with BP 100). One subgroup corresponds to A. haughtonii in which two clusters are present, characterized by variation of up to 24 substitutions. The second subgroup bifurcates into two taxa: A. cordatum (BP 100) and A. phillipsianum (BP 81). The intermediate forms (A. sp.) turn up either in one of the sister taxa (A. cordatum or A. phillipsianum) or in both species (CV398b). The cordatum clade is homogenous; there are no cytological or geographical differences between the 15 specimens from different regions in eastern and southern Africa. The phillipsianum clade shows more variation, with a subcluster (BP 100) grouping the specimens from Ethiopia and Socotra.


    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Cytological results
Although new ploidies within several taxa are described elsewhere (e.g., Van den heede and Viane, 2002 ; Van den heede et al., 2002 ), we want to stress our cytological proof of the occurrence of tetraploid A. ceterach on Gran Canaria and Tenerife (C. J. Van den heede et al., unpublished data). Our study of micromorphology, cytology, and isozymes has shown that besides A. aureum and A. octoploideum, A. lolegnamense and A. ceterach occur in the Canarian Archipelago (Van den heede and Viane, 2002 ; C. J. Van den heede and R. L. L. Viane, unpublished data), and this has important implications for the assessment of the ancestry of tetraploid A. aureum.

Patterns for trnL-F/rbcL in Asplenium subgenus Ceterach
Intraspecific variation in trnL-F DNA was low and limited to single substitutions within some accessions in most species. Because double signals were never found in our trnL-F intergenic spacer and rbcL sequences, plastid inheritance must be at least uniparental (Soltis et al., 1992 ; Sewell et al., 1993 ). Maternal plastid inheritance has generally been recorded in ferns (Stein and Barrington, 1990 ; Gastony and Yatskievych, 1992 ; Vogol et al., 1998 ; Van den heede et al., 2002 ). Consequently, reticulate patterns of evolution cannot be detected. Reticulate evolution in Asplenium was documented in detail by e.g., Wagner (1954) , Lovis (1977) , and Reichstein (1981 , 1984 ).

Internal transcribed spacer in Asplenium subgenus Ceterach
The ITS characteristics of Asplenium fall within the range of variation noted for other land plants, but the length of each of the spacers in these taxa was shorter than those reported by Gastony and Rollo (1998) . The total length of the ITS regions (including the 5.8S gene) varies from 600 to 700 bp, which agrees with the range found across the angiosperms (Soltis and Soltis, 1998 ). Similarly, the GC content does not exceed the variation (45–75%) found within the angiosperms (Hershkovitz et al., 1999 ).

Although generally rare, intraspecific nrITS variability was reported by Sytsma and Schaal (1990) , who detected differences of up to 100 bp among specimens of populations of Lisianthus skinneri (Gentianaceae). According to Smith (1976) large indels may be due to unequal crossing over (believed to be one of the driving forces of concerted evolution) of rDNA repeat units. Baldwin (1993) found intraspecific ITS sequence divergence in populations of Calycadenia (Asteraceae). Similarly, we found intraspecific variation (indels) in the two nuclear ribosomal spacers, ITS1 and ITS2, but not in the conserved 5.8S gene. Limited interspecific, substitutional variation in the 5.8S gene was observed among the three African species but not among the Eurasian-Macaronesian species.

Because gene conversion has not eliminated one allele of the biparentally inherited nuclear ribosomal ITS regions of the allopolyploids, their reticulate origins can be demonstrated. Autopolyploids show a single copy type, whereas allopolyploids are characterized by two, and according to A. Coleman (Brown University, personal communication) these are all functional ITS copies, which is also indicated by the uniform length and low numbers of substitutions observed in the 5.8S gene.

Gastony and Rollo (1998) did not report heterogeneity among ITS copies produced from a single accession. They do not give cytological data but most taxa used in their ITS analysis were diploids or autopolyploids, and gene conversion may have eliminated all but a single ITS allele in any allopolyploids. However, our results are concordant with those of Hoot and Taylor (2001) , showing different copy types in allopolyploids and interspecific hybrids in Isoetes.

Polyphyly of Asplenium subgenus Ceterach
Only Moore (1857) excluded the strictly African species from A. subgenus Ceterach because of their free venation pattern, which contrasts with the anastomosing veins in the Eurasian-Macaronesian complex. Among "modern" pteridologists (e.g., Copeland, 1947 ; Pichi Sermolli, 1979 ), it is generally accepted that A. subg. Ceterach is natural (monophyletic), although this was never evaluated in a phylogenetic context. However, in our plastid analyses, all trees using Dennstaedtia as outgroup show that A. subg. Ceterach is polyphyletic. Moreover, in both our trnL-F and ITS matrices, species within each main group were easily alignable (i.e., required few gaps), but the Eurasian-Macaronesian group was more time-consuming to align with the African group and necessitated many more gaps, some of which were ambiguous in position.

Most comparative studies of rbcL and trnL-F data have produced highly similar but not always identical trees (e.g., Chase et al., 2000 ; Lledo et al., 2000 ; Richardson et al., 2000 ). If we accept the rbcL and trnL-F trees as an accurate reflection of phylogenetic patterns, then the following characters, present in most (99%) Asplenium species, are likely to be plesiomorphic: a pinnate lamina, subglabrous fronds, free venation, and an indusium. Because the two major groups of Ceterach taxa are not each other's closest relatives, two synapomorphies must have originated in parallel: (1) a pinnatisect lamina and (2) a dense scale cover. Adaptation to xeric habitats most likely resulted in the multiple origin of a dense scale cover. The presence of these extra scales made a protective indusium redundant; this consequently became reduced (in the European-Macaronesian species) or totally disappeared (in the African taxa). There is, however, no explanation why the typical pinnatisect lamina (Fig. 6) should have been evolved twice within the family.



View larger version (23K):
[in this window]
[in a new window]
 
Fig. 6. Concept of relationships based on nuclear ITS data of the Asplenium taxa studied. Putative ancestral diploids, "A. semi-aureum" and "A. proto-haughtonii," may be extinct (they are at least undiscovered) and gave rise to the extant taxa via polyploidization and/or hybridization followed by chromosome doubling (sterile hybrids not shown). Thick lines and arrows indicate chromosome doubling

 
The topology of the ITS tree is more in agreement with previous morphological interpretations than the plastid tree, and placement of A. dalhousiae as sole sister species to A. aureum in the plastid trees seems unreasonable to us on morphological grounds. It would imply the secondary loss in A. dalhousiae of four synapomorphies shared by A. ceterach and A. aureum: anastomosing veins, a reduced indusium, dense scale cover, and identical ceterachoid scales. It is more plausible that there was an ancient transfer of plastid DNA from an unknown diploid, implicated in the origin of A. aureum, to A. dalhousiae. This unknown progenitor of A. aureum is hereafter provisionally called "A. semi-aureum" (Fig. 6).

All molecular data collected to date reveal the polyphyly of Asplenium subgenus Ceterach, implicating homoplasy in the lamina shape and the dense scale cover, characters previously used to circumscribe this "taxon." The taxonomic consequences are that either A. subg. Ceterach should be limited to the Eurasian and Macaronesian species (excluding A. dalhousiae) and a new subgenus be created to accommodate the African taxa or that distinction at subgeneric level is abandoned altogether.

Comparison of trnL-F/rbcL and ITS phylogenetic trees
Because the ITS trees could not be rooted with an outgroup, they are not directly comparable to the rooted plastid trees. However, the plastid trees are congruent with the ITS tree in five fundamental ways: (1) recognition of six identical groups, i.e., A. dalhousiae, A. ceterach, "A. semi-aureum," A. cordatum, A. phillipsianum, and A. haughtonii; (2) division of the former subgenus Ceterach into two distinct groups, a Eurasian-Macaronesian and a strictly African alliance; (3) the position of A. dalhousiae in the ingroup; (4) a lack of molecular divergence in A. cordatum despite its morphological variability; and (5) clustering of autopolyploids with their diploid progenitor in the A. ceterach complex, in A. cordatum, and in A. phillipsianum.

The plastid and ITS trees differ in their placements of A. haughtonii and A. dalhousiae. These are "hard incongruences" (Seelanan et al., 1997 , p. 286), and presumably such incongruence is due to different evolutionary histories of these regions, which evolved at their own rates (Page and Holmes, 1998 ; Wendel and Doyle, 1998 ). Various phenomena may lead to discordant patterns and include processes that operate at the organismal level as well as at the molecular level (reviewed in Wendel and Doyle, 1998 ).

The cause of the incongruence between taxa previously included in A. subgenus Ceterach remains unclear. However, it is possible that ancient hybridization including plastid DNA capture unaccompanied by nrDNA introgression might have produced this conflict (Rieseberg and Brunsfeld, 1992 ; Wendel and Doyle, 1998 ). This cytoplasmic "introgression" involves ancient hybridization among lineages that presently have separate ranges and/or are even incompatible. Evidence for the existence of this phenomenon is seen in Gossypium gossypioides (Wendel et al., 1995 ) and Allium section Cepa (van Raamsdonk et al., 1997 ). The occurrence of plastid capture is generally seen as a disadvantage for the use of plastid DNA to study phylogenetic relationships at lower taxonomic levels (Soltis and Soltis, 1998 ), but it also provides additional valuable evidence if hybridization is known to have occurred on the basis of other information.

Further results are discussed by clade or taxon.

Asplenium aureum clade
The aureum clade contains three Macaronesian species, A. aureum, A. lolegnamense, and A. octoploideum, sharing identical nuclear ("semi-aureum") and plastid sequences. Tetraploid A. aureum, endemic to the Canary Islands, was suggested previously (not fully confirmed) to be allotetraploid (Viane and Reichstein, 1992 ). From the ITS pattern, we know that one of its diploid progenitors was A. javorkeanum: one of the ITS copies of A. aureum retrieved by cloning is nearly identical to the ITS copy of the autotetraploid A. ceterach from Tenerife. However, we still lack a diploid that carries the other ITS allele found in A. aureum. For trnL-F/rbcL, A. aureum has a sequence type that is similar to that of A. dalhousiae but not identical. We thus hypothesize the existence of a now extinct or at least unknown diploid ("A. semi-aureum"; Fig. 6) that donated its ITS allele and plastid genome to A. aureum. Our results from plastid trnL-F/rbcL also support the hypothesis that A. aureum is an ancestor of both allohexaploid A. lolegnamense (Gibby and Lovis, 1989 ) and allo-octoploid A. octoploideum (Manton et al., 1986 ). Asplenium lolegnamense produced two ITS alleles, one from "A. semi-aureum" and another nearly identical to A. ceterach from Tenerife.

These results not only document reticulate evolution within subgenus Ceterach sensu stricto, but also support the allopolyploid status of A. aureum and A. lolegnamense (Fig. 5). This confirms our interpretation of morphological, cytological, and isozyme data after the discovery of true A. ceterach on the Canary Islands (Van den heede and Viane, 2002 ).

Asplenium dalhousiae clade
In addition to the large population in southwestern Asia, diploid A. dalhousiae is also known from a few localities in Ethiopia, northwestern Mexico, and southwestern USA (Arizona). This disjunct distribution, which is similar to that of Asplenium exiguum Bedd. and once led Clute (1911) to hypothesize a multiple origin for this species, is now best explained by long-distance dispersal (Moran and Smith, 2001 ; Parris, 2001 ). According to our analyses of the trnL-F/rbcL data, the A. dalhousiae clade is sister to the A. aureum clade. This topology is unexpected because it was generally accepted among pteridologists (e.g., Copeland, 1947 ; Pichi Sermolli, 1979 ) that the members of A. subg. Ceterach originated from an ancestor close to A. dalhousiae. The trnL-F/rbcL results show that A. ceterach, A. dalhousiae, and A. aureum form a clade with A. dalhousiae sister to "A. semi-aureum," the unknown diploid ancestor of A. aureum, A. lolegnamense, and A. octoploideum. Such a phylogeny would imply the secondary loss (reversal) of four apomorphic morphological ceterachoid characters in A. dalhousiae: (1) anastomosing venation, (2) indusium, (3) the dense scale cover, and (4) scale type. These synapomorphies are shared by A. ceterach and A. aureum. Given that the typical ceterach condition prevails in the allotetraploid A. aureum, we assume that the diploid "A. semi-aureum" that gave rise to it had the typical ceterachoid characters, but it most likely had much wider fronds than in either A. ceterach or A. dalhousiae (Fig. 6). In contrast, according to the ITS results, "A. semi-aureum" is sister to A. ceterach, which has at least two divergent alleles, and this group is then a distant sister to A. dalhousiae. This relationship makes more sense morphologically because it does not involve the four reversals mentioned above. The topological and divergence inconsistencies in the trnL-F/rbcL and ITS trees indicate that the evolutionary histories of the nuclear and plastid genomes in this clade are different: putative ancient hybridization and plastid capture, in which the unknown diploid progenitor of A. aureum hybridized with A. dalhousiae, may have left the latter with plastid sequences similar to those of "A. semi-aureum." Additional nuclear markers should be evaluated to determine if this incongruence relates to more of the nuclear genome than ITS; it could also be that the ITS copy obtained from A. dalhousiae is a paralogue to that retrieved from the A. ceterach group.

Asplenium ceterach clade
Asplenium ceterach is an autotetraploid that originated via chromosome doubling in the morphologically identical diploid A. javorkeanum (see introduction). Asplenium ceterach occurs mainly in Europe, southwestern Asia to the Himalayas, and more rarely in northern Africa. The allohexaploid A. punjabense originated by chromosome doubling of a triploid hybrid between diploid A. dalhousiae and tetraploid A. ceterach and is endemic to the Punjab area in India and Pakistan (Bir et al., 1985 ). Such an origin is congruent with the two ITS alleles obtained from cloning (Fig. 5). The trnL-F sequence of A. punjabense (TR7675) is similar to that of a tetraploid A. ceterach specimen (CV494) from Italy. We were unable to investigate any A. ceterach from the Punjab where A. punjabense is endemic, so the match might be better with a specimen from there. For trnL-F, the "A. ceterach clade" includes all investigated taxa of the A. ceterach complex (2x, 3x, 4x, 6x), A. xchasmophilum, and A. punjabense. For both the hexaploid hybrids (CV682a and CV 690a) between tetraploid A. ceterach and octoploid A. octoploideum from Tenerife (Van den heede and Viane, 2002 ) and A. punjabense, autotetraploid A. ceterach is probably the maternal parent. We found almost no intraspecific trnL-F variation in A. ceterach, whereas ITS variation was much higher.

Asplenium cordatum
Lovis (in Bir et al., 1985 ) found A. cordatum to be tetraploid but did not mention either auto- or allopolyploidy. The species is relatively common in rocky habitats in the drier parts of eastern and southern Africa. Both trnL-F and ITS data cluster all A. cordatum samples from all 18 localities in Africa and one intermediate plant (CV961) in which A. cordatum putatively served as maternal parent. Although this species shows considerable (micro)morphological (scale structure and leaf silhouette) variability, intraspecific trnL-F intergenic spacer variation is limited to 1–3 substitutions, whereas no intraspecific ITS divergence was detected, and neither cytological, nor geographical subdivisions are present. Most trnL-F variation (three substitutions; Fig. 3) is found among specimens from Ethiopia. The absence of length and sequence variation in the ITS region of the tetraploid specimens compared to that of their diploid progenitors supports the hypothesis of an autopolyploid origin.

Asplenium phillipsianum
Bir et al. (1985) reported diploid and tetraploid specimens within this taxon, although its exact polyploid status was not investigated. In East and South Africa, A. phillipsianum is less common than A. cordatum. However, it apparently is the only species occurring on the Indian Ocean islands of Socotra, Madagascar, and Réunion. The trnL-F clade includes 13 A. phillipsianum specimens and three plants intermediate in morphology between A. cordatum and A. phillipsianum (CV357a, CV398b, CV605), in which it is likely that A. phillipsianum was the female parent. In the ITS tree, A. phillipsianum shows some geographic variation with a subcluster grouping specimens from Ethiopia and Socotra. The Ethiopian specimens are the most divergent for both trnL-F and ITS data.

African intermediates
In the ITS tree, several morphologically intermediate specimens have ITS copies that turn up in A. cordatum (CV961), A. phillipsianum (CV357a, CV605), or in both species clusters (CV398b). Based on the ITS and chromosome results, the Kenyan tetraploid specimen (CV398b) must be allotetraploid and probably originated via hybridization between A. phillipsianum and A. cordatum. It thus can be considered to represent a new taxon on the African continent. This allotetraploid may have originated either by chromosome doubling in the sterile diploid hybrid between diploids of A. phillipsianum and A. cordatum (conventional alloploidy) or by direct hybridization between the autotetraploids (delayed alloploidy or allohomoploidy). Lovis (1977) introduced the term "delayed allotetraploidy" to denote the phenomenon whereby a fertile allotetraploid hybrid originates from a cross between two different autotetraploid species. At least three European taxa have originated this way (Reichstein, 1984 ). More recently, Haufler (1997) used the term "allohomoploidy" to define hybrid speciation that does not involve chromosome doubling.

Within the African A. cordatum/phillipsianum complex both conventional and delayed allotetraploidy could have taken place because the diploid as well as tetraploid parents are still extant. However, because diploid A. phillipsianum is rare in (South) Africa whereas the tetraploid cytotypes of both A. cordatum and A. phillipsianum are common in Africa, we consider hybridization by delayed allotetraploidy the most probable explanation (Fig. 6). The subtle gross morphological differences between A. cordatum, A. phillipsianum, and intermediate forms (Fig. 6) are supported by the molecular results. Additional micromorphological studies to distinguish these taxa are in progress.

Asplenium haughtonii
The Saint Helenian endemic, A. haughtonii, was reported to be tetraploid by Lovis in Bir et al. (1985) . Several authors (Hooker, 1860 , 1861 ; Sim, 1915 ; Tardieu-Blot, 1958 ; Burrows, 1990 ) considered this plant conspecific with A. cordatum, although both species are easy to distinguish morphologically (Fig. 6). The 12 specimens studied have identical trnL-F sequences. According to the trnL-F and rbcL results, plastid DNA of A. haughtonii is more similar to that of A. phillipsianum than the latter is to that of A. cordatum. The ITS results show a tree with A. haughtonii sister to the pair A. phillipsianum/A. cordatum. This pattern is better supported by morphological and isozyme data (C. J. Van den heede, unpublished data), in which A. haughtonii comes out as a species distinct from its continental African relatives. In any event, the three African entities form distinct clusters, which should be recognized as distinct species.

General conclusions
Using DNA sequences, we were able to clarify interspecific relationships within A. subg. Ceterach (Fig. 6), evaluate earlier hypotheses about the phylogenetic patterns in this group, and document reticulate evolution in the production of some taxa. Morphological, cytological, and isozyme studies (C. J. Van den heede and R. L. L. Viane, unpublished data) provided hypotheses about the number of taxa and the ancestry of species within A. "subg. Ceterach." DNA sequencing has corroborated these hypotheses independently, although we emphasize that cytological knowledge of individual plants was needed to accurately interpret the molecular results.

Our DNA results support the existence of several strictly African taxa. Asplenium cordatum and A. phillipsianum are distinct species or species complexes consisting of diploid and autopolyploid taxa. A third undescribed allotetraploid species has originated from A. cordatum and A. phillipsianum. These three taxa can also be recognized on morphological grounds. Furthermore, our results show that A. haughtonii from Saint Helena is a distinct endemic species apart from its African relatives, A. phillipsianum and A. cordatum.

Reticulate evolution could be demonstrated because gene conversion apparently did not eliminate the divergent nrITS alleles in the allopolyploids. Uniparentally inherited plastid DNA only discloses the maternal lineage, whereas in this case biparentally inherited nuclear ribosomal ITS sequences detected the allopolyploid status of A. aureum, A. lolegnamense, and A. punjabense and indicated that A. ceterach was one of their progenitors. The trnL-L and rbcL results support the hypotheses that A. aureum is an ancestor of both A. lolegnamense and A. octoploideum. Independently of our other studies, these DNA analyses demonstrate that A. ceterach occurs in the Canary Islands.

This study explored the usefulness of DNA sequence data for resolving interspecific relationships within Aspleniaceae. Because A. subg. Ceterach in its current circumscription is polyphyletic, increased taxon sampling is needed to reveal the position within the family of these two clades and identify the sister group of the strictly African taxa. Future work will include additional molecular data, both plastid and nuclear, of many more species, in combination with study of cytology, isozymes, and (micro)morphological characters. More effort needs to be focused on the apparent disagreement of the plastid and ITS patterns; in particular additional evidence from the nuclear genome should be sought. Nomenclatural changes indicated by our results will be published elsewhere.


    FOOTNOTES
 
1 The authors thank R. Cranfill for providing trnL-F sequences of some outgroups; A. Coleman for checking the secondary structure stability of the different ITS copy types; N. Murakami for giving information about rbcL sequencing; N. Villarroel for his help at Ghent University; R. Johns and the "Flora of Tropical East Africa" committee for hosting the visit of the first author at the Jodrell Laboratory; the Royal Botanic Gardens Kew (UK) for funding the lab work; Ghent University for financial support; and the management of Ghent Botanical Garden (Belgium) for providing greenhouse facilities. We appreciate the constructive comments on an earlier draft of this manuscript, given by the reviewers. A portion of the fieldwork in Africa was sponsored by the Belgian FSR "Indian Ocean project" G.0024.96 N. Back

4 Author for reprint requests (ronnie.viane{at}rug.ac.be ) Back


    LITERATURE CITED
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Acquaah G. 1992 Practical protein electrophoresis for genetic research. Dioscorides Press, Portland, Oregon, USA

Arnheim N. 1983 Concerted evolution in multigene families. In M. Nei and R. Koehn [eds.], Evolution of genes and proteins, 38–61. Sinauer, Sunderland, Massachusetts, USA

Asmussen C. B. M. W. Chase 2001 Coding and noncoding plastid DNA in palm systematics. American Journal of Botany 88: 1103-1117[Abstract/Free Full Text]

Baldwin B. G. 1992 Phylogenetic utility of the internal transcribed spacers of nuclear ribosomal DNA in plants: an example from the Compositae. Molecular Phylogenetics and Evolution 1: 3-16[CrossRef][Medline]

Baldwin B. G. 1993 Molecular phylogenetics of Calycadenia (Compositae) based on ITS sequences of nuclear ribosomal DNA: chromosomal and morphological evolution reexamined. American Journal of Botany 80: 222-238[CrossRef][Web of Science]

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]

Bir S. S. C. R. Fraser-Jenkins J. D. Lovis 1985 Asplenium punjabense sp. nov. and its significance for the status of Ceterach and Ceterachopsis. Fern Gazette 13: 54-63

Buckler E. S. A. Ippolito T. P. Holtsford 1997 The evolution of ribosomal DNA: divergent paralogues and phylogenetic implications. Genetics 145: 821-832[Abstract]

Burrows J. E. 1990 Southern African ferns and fern allies. Frandsen, Sandton, South Africa

Campbell C. S. M. F. Wojciechowski B. G. Baldwin L. A. Alice M. J. Donoghue 1997 Persistent nuclear ribosomal DNA sequence polymorphism in the Amelanchier agamic complex (Rosaceae). Molecular Biology and Evolution 14: 81-90[Abstract]

Chase M. W. A. De Bruijn A. V. Cox G. Reeves P. J. Rudall M. A. T. Johnson L. E. Eguiarte 2000 Phylogenetics of Asphodelaceae (Aparagales): an analysis of plastid rbcL and trnL-F DNA sequences. Annals of Botany 86: 935-951[Abstract/Free Full Text]

Cheng X. N. Murakami 1998 Cytotaxonomic study of genus Hymenasplenium (Aspleniaceae) in Xishuangbanna, southwestern China. Journal of Plant Research 111: 495-500[CrossRef][Web of Science]

Clute W. N. 1911 The distribution of Asplenium alternans. Fern Bulletin 19: 38-42

Copeland E. B. 1947 Genera Filicum. Chronica Botanica, Waltham, Massachusetts, USA

De Candolle A. P. 1813 Théorie élémentaire de la botanique, ou exposition des principes de la classification naturelle et de l'art de décrire et d'étudier les végétaux. Déterville, Paris, France

Dover G. 1982 Molecular drive: a cohesive mode of species evolution. Nature 299: 111-117[CrossRef][Medline]

Doyle J. J. J. L. Doyle 1987 A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochemical Bulletin 19: 11-15

Fan C. Q.-Y. Xiang 2001 Phylogenetic relationships within Cornus (Cornaceae) based on rDNA sequences. American Journal of Botany 88: 1131-1138[Abstract/Free Full Text]

Fitch W. M. 1971 Toward defining the course of evolution: minimum change for a specific tree topology. Systematic Zoology 20: 406-416[CrossRef][Web of Science]

Fujii N. K. Ueda Y. Watano T. Shimizu 1997 Intraspecific sequence variation of chloroplast DNA in Pedicularis chamissonis Steven (Scrophulariaceae) and geographic structuring of the Japanese "Alpine" plants. Journal of Plant Research 110: 195-207[CrossRef][Web of Science]

Fujii N. K. Ueda Y. Watano T. Shimizu 2001 Two genotypes of Pedicularis chamissonis (Scrophulariaceae) distributed at Mt. Gassan, Japan: additional genetic and morphological studies. Journal of Plant Research 114: 133-140[CrossRef][Web of Science]

Gastony G. J. D. R. Rollo 1998 Cheilanthoid ferns (Pteridaceae: Cheilanthoideae) in the southwestern United States and adjacent Mexico—a molecular phylogenetic reassessment of generic lines. Aliso 17: 131-144

Gastony G. J. G. Yatskievych 1992 Maternal inheritance of the chloroplast and mitochondrial genomes in Cheilanthoid ferns. American Journal of Botany 79: 716-722[CrossRef][Web of Science]

Gibby M. J. D. Lovis 1989 New ferns of Madeira. Fern Gazette 13: 285-290

Graur D. W.-H. Li 2000 Fundamentals of molecular evolution, 2nd ed. Sinauer, Sunderland, Massachusetts, USA

Harris S. A. R. Ingram 1991 Chloroplast DNA and biosystematics: the effects of intraspecific diversity and plastid transmission. Taxon 40: 393-412

Hasebe M. M. Ito R. Kofuji K. Ueda K. Iwatsuki 1993 Phylogenetic relationships of ferns deduced from rbcL gene sequence. Journal of Molecular Evolution 37: 476-482[CrossRef][Web of Science][Medline]

Hasebe M. T. Omori M. Nakazawa T. Sano M. Kato K. Iwatsuki 1994 rbcL gene sequences provide evidence for the evolutionary lineages of leptosporangiate ferns. Proceedings of the National Academy of Sciences, USA 91: 5730-5734[Abstract/Free Full Text]

Hasebe M. et al., 1995 Fern phylogeny based on rbcL nucleotide sequences. American Fern Journal 85: 134-181[CrossRef]

Haufler C. H. 1997 Modes and mechanisms of speciation in pteridophytes. In K. Iwatsuki and P. Raven [eds.], Evolution and diversification of land plants, 291–307. Springer-Verlag, Tokyo, Japan

Haufler C. H. E. A. Hooper J. P. Therrien 2000 Modes and mechanisms of speciation in pteridophytes: implications of contrasting patterns in ferns representing temperate and tropical habitats. Plant Species Biology 15: 223-236[CrossRef]

Hauk W. D. C. R. Parks M. W. Chase 1996 A comparison between trnL-F intergenic spacer and rbcL DNA sequence data: an example from Ophioglossaceae. American Journal of Botany 83: 126 (Abstract)

Heitz E. 1925 Der Nachweis der Chromosomen. Vergleichende Studien über ihre Zahl, Größe und Form im Pflanzenreich I. Zeitschrift für Botanik 18: 625-681

Heitz E. 1950 Elemente der Botanik. Eine Anleitung zum Studium der Pflanze durch Beobachtungen und Versuche an Crepis capillaris (L.) Wallr. Springer-Verlag, Wien, Austria

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, 268–326. Taylor and Francis, London, UK

Hooker W. J. 1860 Species Filicum, vol. III, Lomaria-Actiniopteris. W. Pamplin, London, UK

Hooker W. J. 1861 A second century of ferns. W. Pamplin, London, UK

Hoot S. B. W. C. Taylor 2001 The utility of nuclear ITS, a LEAFY homolog intron, and chloroplast atpB-rbcL spacer region data in phylogenetic analyses and species delimitation in Isoetes. American Fern Journal 91: 166-177[CrossRef]

Kim K.-J. R. K. Jansen 1994 Comparisons of phylogenetic hypotheses among different data sets in dwarf dandelions (Krigia, Asteraceae): additional information from internal transcribed spacer sequences of nuclear ribosomal DNA. Plant Systematics and Evolution 190: 157-185[CrossRef][Web of Science]

King M. 1993 Species evolution, the role of chromosome change. Cambridge University Press, Cambridge, UK

Kramer K. U. R. L. L. Viane 1990 Aspleniaceae. In K. U. Kramer and P. S. Green [eds.], Pteridophytes and gymnosperms, 52–56. Springer-Verlag, Berlin, Germany

Kranz H. D. V. A. R. Huss 1996 Molecular evolution of pteridophytes and their relationship to seed plants: evidence from complete 18S rRNA gene sequences. Plant Systematics and Evolution 202: 1-11[CrossRef][Web of Science]

Lledo M. D. M. B. Crespo A. V. Cox M. F. Fay M. W. Chase 2000 Polyphyly of Limoniastrum (Plumbaginaceae): evidence from DNA sequences of plastid rbcL, trnL intron and trnL-F intergene spacer. Botanical Journal of the Linnean Society 132: 175-191[CrossRef]

Löve A. D. Löve R. E. G. Pichi Sermolli 1977 Cytotaxonomical atlas of the pteridophytes. J. Cramer, Vaduz, Liechtenstein

Lovis J. D. 1977 Evolutionary patterns and processes in ferns. In R. D. Preston and H. W. Woolhouse [eds.], Advances in botanical research, vol. 4, 229–415. Academic Press, London, UK

Luerssen C. 1889 Die Farnpflanzen oder Gefässbündelkryptogamen (Pteridopyta). In A. Grunow, A. Fischer, F. Hauck, G. Limpricht, C. Luerssen, W. Migula, H. Rehm, P. Richter, and G. Winter [eds.], Dr. L. Rabenhorst's Kryptogamen-Flora von Deutschland, Oesterreich und der Schweiz, 2nd ed. vol. 3, 1–906. Eduard Kummer, Leipzig, Germany

Manhart J. R. 1995 Chloroplast 16S rDNA sequences and phylogenetic relationships of fern allies and ferns. American Fern Journal 85: 182-192[CrossRef][Web of Science]

Manton I. 1950 Problems of cytology and evolution in the Pteridophyta. Cambridge University Press, Cambridge, UK

Manton I. J. D. Lovis G. Vida M. Gibby 1986 Cytology of the fern flora of Madeira. Bulletin of the British Museum (Natural History) Botany series 15: 123-161

Mettenius G. 1859 Über einige Farngattungen. VI. Asplenium. H.L. Brönner, Frankfurt am Main, Germany

Mitui K. N. Murakami K. Iwatsuki 1989 Chromosomes and systematics of Asplenium sect. Hymenasplenium (Aspleniaceae). American Journal of Botany 76: 1689-1697[CrossRef][Web of Science]

Moore T. 1857 Index Filicum: a synopsis, with characters, of the genera, and an enumeration of the species of ferns. W. Pamplin, London, UK

Moran R. C. A. R. Smith 2001 Phytogeographic relationships between neotropical and African-Madagascan pteridophytes. Brittonia 53: 304-351[Web of Science]

Murakami N. 1995 Systematics and evolutionary biology of the fern genus Hymenasplenium (Aspleniaceae). Journal of Plant Research 108: 257-268[CrossRef][Web of Science]

Murakami N. R. C. Moran 1993 Monograph of the neotropical species of Asplenium sect. Hymenasplenium (Aspleniaceae). Annals of the Missouri Botanical Garden 80: 1-38

Murakami N. S. Nogami M. Watanabe K. Iwatsuki 1999a Phylogeny of Aspleniaceae inferred from rbcL nucleotide sequences. American Fern Journal 89: 232-243[CrossRef]

Murakami N. M. Watanabe J. Yokoyama Y. Yatabe H. Iwasaki S. Serizawa 1999b Molecular taxonomic study and revision of the three Japanese species of Asplenium sect. Thamnopteris. Journal of Plant Research 112: 15-25[CrossRef][Web of Science]

Murakami N. J. Yokoyama X. Cheng H. Iwasaki R. Imaichi K. Iwatsuki 1998 Molecular {alpha}-taxonomy of Hymenasplenium obliquissimum complex (Aspleniaceae) based on rbcL sequence comparisons. Plant Species Biology 13: 51-56

Page R. D. M. E. C. Holmes 1998 Molecular evolution. A phylogenetic approach. Blackwell Science, Oxford, UK

Pahnke J. V. Goremykin V. Bobrova A. Troitsky A. Antonov W. Martin 1996 Utility of rDNA internal transcribed spacer sequences from the inverted repeat of chloroplast DNA in pteridophyte molecular phylogenetics. In J. M. Camus, M. Gibby, and R. J. Johns [eds.], Pteridology in perspective, 217–230. Royal Botanic Gardens, Kew, UK

Parris B. S. 2001 Circum-Antarctic continental distribution patterns in pteridophyte species. Brittonia 53: 270-283[Web of Science]

Pichi Sermolli R. E. G. 1972 Fragmenta Pteridologiae III. Webbia 27: 389-460

Pichi Sermolli R. E. G. 1979 A survey of the pteridological flora of the Mediterranean Region. Webbia 34: 175-242

Pryer K. M. A. R. Smith J. E. Skog 1995 Phylogenetic relationships of extant ferns based on evidence from morphology and rbcL sequences. American Fern Journal 85: 205-282[CrossRef][Web of Science]

Rasbach H. K. Rasbach H. W. Bennert 1987 xAsplenoceterach barrancense Bennert et Meyer (Aspleniaceae, Pteridophyta)—Neufunde und cytologische Untersuchungen. Farnblätter 17: 3-16

Raubeson L. A. D. B. Stein 1995 Insights into fern evolution from mapping chloroplast genomes. American Fern Journal 85: 193-204

Reichstein T. 1981 Hybrids in European Aspleniaceae (Pteridophyta). Botanica Helvetica 91: 89-139

Reichstein T. 1984 Aspleniaceae. In H. J. Conert, U. Hamann, W. Schultze-Motel, and G. Wagenitz [eds.], Gustav Hegi Illustrierte Flora von Mitteleuropa. Band 1, Pteridophyta. Teil 1, 211–275. P. Parey, Berlin, Hamburg, Germany

Richardson J. E. M. F. Fay Q. C. B. Cronk D. Bowman M. W. Chase 2000 A phylogenetic analysis of Rhamnaceae using rbcL and trnL-F plastid DNA sequences. American Journal of Botany 87: 1309-1324[Abstract/Free Full Text]

Rieseberg L. H. S. J. Brunsfeld 1992 Molecular evidence and plant introgression. In D. E. Soltis, P. S. Soltis, and J. J. Doyle [eds.], Molecular systematics of plants, 151–176. Chapman and Hall, New York, New York, USA

Sang T. D. J. Crawford T. F. Stuessy 1995 Documentation of reticulate evolution in peonies (Paeonia) using internal transcribed spacer sequences of nuclear ribosomal DNA: implications for biogeograpgy and concerted evolution. Proceedings of the National Academy of Sciences, USA 92: 6813-6817[Abstract/Free Full Text]

Schelpe E. A. C. L. E. 1970 Pteridophyta. In A. W. Exell and E. Launert [eds.], Flora Zambesiaca. Crown Agents, London, UK

Schelpe E. A. C. L. E. N. C. Anthony 1986 Pteridophyta. In O. A. Leistner [ed.], Flora of Southern Africa. Government Printer, Pretoria, South Africa

Schulze G. J. Treutlein M. Wink 2001 Phylogenetic relationships between Asplenium bourgaei (Boiss.) Milde and A. jahandiezii (Litard.) Rouy inferred from morphological characters and rbcL sequences. Plant Biology 3: 364-371[CrossRef][Web of Science]

Seelanan T. H. Schnabel J. F. Wendel 1997 Congruence and consensus in the cotton tribe. Systematic Botany 22: 259-290[CrossRef][Web of Science]

Sewell M. M. Y.-L. Qiu C. R. Parks M. W. Chase 1993 Genetic evidence for trace paternal transmission of plastids in Liriodendron and Magnolia (Magnoliaceae). American Journal of Botany 80: 854-858[CrossRef][Web of Science]

Sim T. R. 1915 The ferns of South Africa containing descriptions and figures of the ferns and fern allies of South Africa. University Press, Cambridge, UK

Smith G. P. 1976 Evolution of repeated DNA sequences by unequal crossover. Science 191: 528-535[Abstract/Free Full Text]

Soltis D. E. P. S. Soltis 1998 Choosing an approach and an appropriate gene for phylogenetic analysis. In D. E. Soltis, P. S. Soltis, and J. J. Doyle [eds.], Molecular systematics of plants II. DNA sequencing, 1–42. Kluwer, Boston, Massachusetts, USA

Soltis D. E. P. S. Soltis B. G. Milligan 1992 Intraspecific chloroplast DNA variation: systematic and phylogenetic implications. In D. E. Soltis, P. S. Soltis, and J. J. Doyle [eds.], Molecular systematics of plants, 117–150. Chapman and Hall, New York, New York, USA

Soltis P. S. J. J. Doyle D. E. Soltis 1992 Molecular data and polyploid evolution in plants. In D. E. Soltis, P. S. Soltis, and J. J. Doyle [eds.], Molecular systematics of plants, 177–201. Chapman and Hall, New York, New York, USA

Soltis P. S. D. E. Soltis 1991 Multiple origins of the allotetraploid Tragopogon mirus (Compositae): rDNA evidence. Systematic Botany 16: 407-413[CrossRef][Web of Science]

Stebbins G. L. 1989 Introduction. In D. E. Soltis and P. S. Soltis [eds.], Isozymes in plant biology, 1–4. Dioscorides Press, Portland, Oregon, UK

Stein D. B. D. S. Barrington 1990 Recurring hybrid formation in a population of Polystichum x potteri: evidence from chloroplast DNA comparisons. Annals of the Missouri Botanical Garden 77: 334-339[CrossRef][Web of Science]

Swofford D. L. 2002 PAUP*: phylogenetic analysis using parsimony (*and other methods), version 4.0b10. Sinauer, Sunderland, Massachusetts, USA

Sytsma K. J. B. A. Schaal 1990 Ribosomal DNA variation within and among individuals of Lisianthus (Gentianaceae) populations. Plant Systematics and Evolution 170: 97-106[CrossRef][Web of Science]

Taberlet P. L. Gielly G. Pautou J. Bouvet 1991 Universal primers for amplification of three non-coding regions of chloroplast DNA. Plant Molecular Biology 17: 1105-1109[CrossRef][Web of Science][Medline]

Tardieu-Blot M.-L. 1958 Polypodiacées (sensu lato). In H. Humbert [ed.], Flore de Madagascar et des Comores (plantes vasculaires). Tome I. Gouvernement général de Madagascar, Muséum National d'Histoire Naturelle, Paris, France

Van den heede C. J. R. L. L. Viane 2002 New species and new hybrids in Asplenium subgenus Ceterach (Aspleniaceae). GEP News 9: 1-4

Van den heede C. J. S. Pajarón E. Pangua R. L. L. Viane 2002 A new species and a new hybrid of Asplenium (Aspleniaceae) from Cyprus and evidence of their origin. Belgian Journal of Botany 135: 92-116[Web of Science]

van Raamsdonk L. W. D. M. P. Smiech J. M. Sandbrink 1997 Introgression explains incongruence between nuclear and chloroplast DNA-based phylogenies in Allium section Cepa. Botanical Journal of the Linnean Society 123: 91-108[CrossRef]

Varadaraj K. D. M. Skinner 1994 Denaturants or cosolvents improve the specificity of PCR amplification of a GC-rich DNA using genetically engineered DNA polymerases. Gene 140: 1-5[CrossRef][Web of Science][Medline]

Viane R. L. L. A. C. Jermy J. D. Lovis 1993 Asplenium. In T. G. Tutin, N. A. Burges, A. O. Chater, J. R. Edmondson, V. H. Heywood, D. M. Moore, D. H. Valentine, S. M. Walters, and D. A. Webb [eds.], Flora Europaea, 2nd ed., vol. 1, 18–23. Cambridge University Press, Cambridge, UK

Viane R. L. L. H. Rasbach K. Rasbach T. Reichstein 1996 Observations on some ferns of Poros and the adjacent parts of the Peloponnesus (Greece). Bocconea 5: 279-300

Viane R. L. L. T. Reichstein 1992 Notes about Asplenium II: some new names and combinations in Asplenium L. (Aspleniaceae, Pteridophyta). Biologisch Jaarboek Dodonaea 59: 157-165

Vida G. 1963 A new Asplenium (Sectio Ceterach) species and the problem of the origin of Phyllitis hybrida (Milde) C. Christ. Acta Botanica Academiae Scientiarum Hungaricae 9: 197-215

Vogol J. C. S. J. Russel F. J. Rumsey J. A. Barrett M. Gibby 1998 Evidence for maternal transmission of chloroplast DNA in the genus Asplenium (Aspleniaceae, Pteridophyta). Botanica Acta 111: 247-249[Web of Science]

Wagner W. H., Jr. 1954 Reticulate evolution in the Appalachian Aspleniums. Evolution 7: 103-118

Wagner W. H., Jr. R. C. Moran C. R. Werth 1993 Aspleniaceae. In N. R. Morin et al. [eds.], Flora of North America north of Mexico, vol. 2, 228–245. Oxford University Press, New York, New York, USA

Wendel J. F. J. J. Doyle 1998 Phylogenetic incongruence: window into genome history and molecular evolution. In D. E. Soltis, P. S. Soltis, and J. J. Doyle [eds.], Molecular systematics of plants. II. DNA sequencing, 265–296. Kluwer, Boston, Massachusetts, USA

Wendel J. F. A. Schnabel T. Seelanan 1995 An unusual ribosomal DNA sequence from Gossypium gossypioides reveals ancient, cryptic, intergenomic introgression. Molecular Phylogenetics and Evolution 4: 298-313[CrossRef][Web of Science][Medline]

White T. J. T. Bruns S. Lee J. Taylor 1990 Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. In M. Innis, D. Gelfand, J. Sninsky, and T. White [eds.], PCR protocols: a guide to methods and applications, 315–322. Academic Press, San Diego, California, USA

Wolf P. G. 1995 Phylogenetic analyses of rbcL and nuclear ribosomal RNA gene sequences in Dennstaedtiaceae. American Fern Journal 85: 306-327[CrossRef]

Wolf P. G. 1996 Pteridophyte phylogenies based on analysis of DNA sequences: a multiple gene approach. In J. M. Camus, M. Gibby, and R. J. Johns [eds.], Pteridology in perspective, 203–215. Royal Botanic Gardens, Kew, UK

Wolf P. G. S. D. Sipes M. R. White M. L. Martines K. M. Pryer A. R. Smith K. Ueda 1999 Phylogenetic relationships of the enigmatic fern families Hymenophyllopsidaceae and Lophosoraceae: evidence from rbcL nucleotide sequences. Plant Systematics and Evolution 219: 263-270[CrossRef][Web of Science]

Wolf P. G. P. S. Soltis D. S. Soltis 1994 Phylogenetic relationships of Dennstaedtioid ferns: evidence from rbcL sequences. Molecular Phylogenetics and Evolution 3: 383-392[CrossRef][Medline]

Yatabe Y. H. Nishida N. Murakami 1999 Phylogeny of Osmundaceae inferred from rbcL nucleotide sequences and comparison of the fossil evidences. Journal of Plant Research 112: 397-404[CrossRef][Web of Science]


Add to CiteULike CiteULike   Add to Complore Complore   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us   Add to Digg Digg   Add to Facebook Facebook   Add to Reddit Reddit   Add to Technorati Technorati   Add to Twitter Twitter    What's this?


This article has been cited by other articles:


Home page
Am. J. Bot.Home page
A. Ebihara, H. Ishikawa, S. Matsumoto, S.-J. Lin, K. Iwatsuki, M. Takamiya, Y. Watano, and M. Ito
Nuclear DNA, chloroplast DNA, and ploidy analysis clarified biological complexity of the Vandenboschia radicans complex (Hymenophyllaceae) in Japan and adjacent areas
Am. J. Botany, September 1, 2005; 92(9): 1535 - 1547.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplemental Data
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Van den heede, C. J.
Right arrow Articles by Chase, M. W.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Van den heede, C. J.
Right arrow Articles by Chase, M. W.
Agricola
Right arrow Articles by Van den heede, C. J.
Right arrow Articles by Chase, M. W.
Social Bookmarking
 Add to CiteULike   Add to Complore   Add to Connotea   Add to Del.icio.us   Add to Digg   Add to Facebook   Add to Reddit   Add to Technorati   Add to Twitter  
What's this?


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS