A pre‐Miocene Irano‐Turanian cradle: Origin and diversification of the species‐rich monocot genus Gagea (Liliaceae)

The Irano‐Turanian (IT) floristic region is considered an important center of origin for many taxa. However, there is a lack of studies dealing with typical IT genera that also occur in neighboring areas. The species‐rich monocot genus Gagea Salisb. shows a center of diversity in IT region and a distribution in adjacent regions, therefore representing a good study object to investigate spatial and temporal relationships among IT region and its neighboring areas (East Asia, Euro‐Siberia, Himalaya, and Mediterranean). We aimed at (a) testing the origin of the genus and of its major lineages in the IT region, (b) reconstructing divergence times, and (c) reconstructing colonization events. To address these problems, sequences of the ribosomal DNA internal transcribed spacer (ITS) region of 418 individuals and chloroplast intergenic spacers sequences (psbA‐trnH, trnL‐trnF) of 497 individuals, representing 116 species from all sections of the genus and nearly its entire distribution area were analyzed. Divergence times were estimated under a random molecular clock based on nrITS phylogeny, which was the most complete data set regarding the representation of species and distribution areas. Ancestral distribution ranges were estimated for the nrITS data set as well as for a combined data set, revealing that Gagea most likely originated in southwestern Asia. This genus first diversified there starting in the Early Miocene. In the Middle Miocene, Gagea migrated to the Mediterranean and to East Asia, while migration into Euro‐Siberia took place in the Late Miocene. During the Pleistocene, the Arctic was colonized and Gagea serotina, the most widespread species, reached North America. The Mediterranean basin was colonized multiple times from southwestern Asia or Euro‐Siberia. Most of the currently existing species originated during the last 3 Ma.     

Evaluation of six candidate DNA barcoding loci in Ficus (Moraceae) of China

Ficus, with about 755 species, diverse habits and complicated co‐evolutionary history with fig wasps, is a notoriously difficult group in taxonomy. DNA barcoding is expected to bring light to the identification of Ficus but needs evaluation of candidate loci. Based on five plastid loci (rbcL, matK, trnH‐psbA, psbK‐psbI, atpF‐atpH) and a nuclear locus [internal transcribed spacer (ITS)], we calculated genetic distances and DNA barcoding gaps individually and in combination and constructed phylogenetic trees to test their ability to distinguish the species of the genus. A total of 228 samples representing 63 putative species in Ficus (Moraceae) of China were included in this study. The results demonstrated that ITS has the most variable sites, greater intra‐ and inter‐specific divergences, the highest species discrimination rate (72%) and higher primer universality among the single loci. It is followed by psbK‐psbI and trnH‐psbA with moderate variation and considerably lower species discrimination rates (about 19%), whereas matK, rbcL and atpF‐atpH could not effectively separate the species. Among the possible combinations of loci, ITS + trnH‐psbA performed best but only marginally improved species resolution over ITS alone (75% vs. 72%). Therefore, we recommend using ITS as a single DNA barcoding locus in Ficus.     

Molecular phylogeny of Aerides (Orchidaceae) based on one nuclear and two plastid markers: a step forward in understanding the evolution of the Aeridinae

Phylogenetic relationships of the orchid genus Aerides (Epidendroideae, Vandeae, Aeridinae) from Southeast Asia were inferred from DNA sequences of one nuclear (nrITS) and two plastid (matK, trnL–trnL–F) regions of 48 taxa (21 Aerides, 25 other Aeridinae, 2 outgroup). Analyses of the combined datasets with parsimony, maximum likelihood and Bayesian methods revealed that Aerides is monophyletic and consists of three well-supported subclades which are only partly in accordance with previous sectional delimitations based on floral characters. The two different flower types in Aerides (hidden versus open spur entrance) seem to have evolved at least twice in geographically distinct areas. The phylogeny presented here is yet another example in Orchidaceae where floral morphology cannot be relied on to reconstruct phylogenetic history but rather is the result of pollinator-driven selection. The Aerides subclades are characterized by three different length classes of the mutation-rich P8 region in the trnL intron. To our knowledge, this is the first time that the P8 region was studied in orchids. The matK gene has been assumed to be a pseudogene in orchids due to occasional occurrence of frameshift indels, low transition/transversion (ts:tv) ratios and low substitution rates at the 3rd codon position. However, matK does not appear to be a pseudogene in Aerides and a comparison with data from other angiosperms suggests that ts:tv ratios and low substitution rates have been overestimated as arguments for a pseudogene status of matK in orchids.     

Molecular identification of species in Prunus sect. Persica (Rosaceae), with emphasis on evaluation of candidate barcodes for plants

Abstract Species of Prunus L. sect. Persica are not only important fruit trees, but also popular ornamental and medicinal plants. Correct identification of seedlings, barks, or fruit kernels is sometimes required, but no reliable morphological characters are available. Nowadays, the technique of DNA barcoding has the potential to meet such requirements. In this study, we evaluated the suitability of 11 DNA loci (atpB‐rbcL, trnH‐psbA, trnL‐F, trnS‐G, atpF‐H, rbcL, matK, rpoB, rpoC1, nad1, and internal transcribed spacer [ITS]) as candidate DNA barcodes for peaches, using samples from 38 populations, covering all the species in sect. Persica. On the whole, the primers worked well in this group and sequencing difficulties were met only in the case of ITS locus. Five loci (rbcL, matK, rpoB, rpoC, and nad1) have very low variation rates, whereas atpB‐rbcL, atpF‐H, trnH‐psbA, trnL‐F and trnS‐G show more variability. The most variable loci, atpB‐rbcL and trnH‐psbA, can distinguish three of the five species. Two two‐locus combinations, atpB‐rbcL+trnL‐F and atpB‐rbcL+atpF‐H, can resolve all five species. We also find that identification powers of the loci are method‐dependent. The NeighborNet method shows higher species identification power than maximum parsimony, neighbor joining, and unweighted pair group method with arithmetic mean methods.     

Phylogenetic relationships between disjunctly occurring groups of Tristicha trifaria (Podostemaceae)

Aim To reveal the phylogeographic relationship of disjunct specimens of Tristicha trifaria (Bory ex Willd.) Spreng., a member of the Podostemaceae river‐weed family, which is distributed exceptionally widely, but disjunctly, in Africa and the Americas. Location Brazil, Mexico, Ghana, Tanzania and Madagascar. Methods The chloroplast matK and rbcL genes, a trnK intron, the trnS‐trnG intergenic spacer (IGS), the two IGSs of trnT‐trnL‐trnF, a trnL intron, and nuclear ribosomal ITS regions were sequenced and analysed. Phylogenetic analyses were conducted using maximum likelihood and maximum parsimony methods. Results The T. trifaria samples analysed were separated into two groups in a rooted tree based on a combined matK/rbcL/ITS dataset; one contained the West African and all of the American samples, and the other contained the East African and Madagascan samples. An unrooted tree obtained from a combined analysis of all the chloroplast DNA and nuclear ITS data showed that a sample from West Africa was sister to an American T. trifaria group. Main conclusions The American and West African T. trifaria are closely related, despite the great distance between their locations. This observation, along with a tree of the whole Tristichoideae subfamily and estimated divergence times, suggests that an ancestor of T. trifaria migrated from Asia to Africa during the early Tertiary, and that this was followed by further westward migration to the Americas at the end of the Miocene or in the early Pliocene.     

Potential use of DNA barcoding for the identification of Salvia based on cpDNA and nrDNA sequences

An effective DNA marker for authenticating the genus Salvia was screened using seven DNA regions (rbcL, matK, trnL–F, and psbA–trnH from the chloroplast genome, and ITS, ITS1, and ITS2 from the nuclear genome) and three combinations (rbcL + matK, psbA–trnH + ITS1, and trnL–F + ITS1). The present study collected 232 sequences from 27 Salvia species through DNA sequencing and 77 sequences within the same taxa from the GenBank. The discriminatory capabilities of these regions were evaluated in terms of PCR amplification success, intraspecific and interspecific divergence, DNA barcoding gaps, and identification efficiency via a tree-based method. ITS1 was superior to the other marker for discriminating between species, with an accuracy of 81.48%. The three combinations did not increase species discrimination. Finally, we found that ITS1 is a powerful barcode for identifying Salvia species, especially Salvia miltiorrhiza.     

Phylogenetic relationships of Iranian Allium sect. Allium (Amaryllidaceae,Allioideae) as inferred from nrDNA ITS,cpDNA rps16 and trnL–F sequences

Allium is a particularly species rich (more than 800 species) and economically important genus, with numerous taxonomic problems at all levels of classification. In this study, we try to uncover the phylogenetic relationships of the common leek Allium ampeloprasum based on selected samples of this species and its putative relatives in A. sect. Allium from Iran. The silica‐dried leaf samples of 56 accessions representing 23 species of Allium were sequenced; 53 sequences of nrDNA ITS, 35 sequences of plastid rps16 and 52 sequences of trnL–F were generated and additional accessions were extracted from GenBank in order to cover all recognized main lineages in the genus. Maximum Parsimony and Bayesian Inference generated similar trees, but the placement of A. ampeloprasum and its relatives differed slightly between the nuclear and the plastid phylogenies. In the nrITS tree, A. ampeloprasum is retrieved as a highly supported clade with A. iranicum, while in the combined plastid tree A. ampeloprasum formed a highly supported clade with A. vineale. This supports the hypothesis of a possible hybrid origin of A. ampeloprasum. Allium iranicum formed a clade in the plastid tree, but was resolved as paraphyletic in the nrITS tree, probably due to presence of multiple non‐concerted copies of nrITS. Close relationships are suggested between the following species: A. aznavense and A. wendelboi with A. talyschense, A. erubescens and A. rotundum with A. scorodoprasum and A. abbasii with A. phanerantherum.     

Induction of multiple shoots in a monopodial orchid hybrid (Aerides vandarum Reichb.f × Vanda stangeana Reichb.f) using thidiazuron and analysis of their genetic stability

The effect of thidiazuron (TDZ) on direct multiple shoot induction in axenic seedlings of a monopodial orchid hybrid Aerides vandarum × Vanda stangeana, using a dual phase culture medium was studied. The culture system consisted of a basal agar solidified half-strength Murashige and Skoog medium overlaid by a liquid fraction of the same composition. Highest regeneration of multiple shoots (15.8 shoots per seedling) was obtained in the medium containing 2% sucrose (w/v) supplemented with 2 mgl⁻¹ TDZ. Genetic stability of the regenerated shoots was assessed using random amplified polymorphic DNA (RAPD), inter simple sequence repeat (ISSR), and restriction fragment length polymorphism of the PCR amplified (PCR-RFLP) nrITS region, as well as those of the coding (matK) and non-coding (trnL-F) regions of the cpDNA. Across the randomly selected mother plant and nine of its regenerated shoots, 2,680 bands were generated by 19 RAPD and 12 ISSR primers, exhibiting monomorphic banding profiles. Homogenous PCR-RFLP profiles were generated for nrITS using four restriction enzymes (REs), matK using five REs, and trnL-F using six REs. These molecular analyses showed no genomic alterations in all regenerated shoots obtained on medium containing 2 mgl⁻¹ TDZ.     

Development of high‐resolution DNA barcodes for Dioscorea species discrimination and phylogenetic analysis

The genus Dioscorea is widely distributed in tropical and subtropical regions, and is economically important in terms of food supply and pharmaceutical applications. However, DNA barcodes are relatively unsuccessful in discriminating between Dioscorea species, with the highest discrimination rate (23.26%) derived from matK sequences. In this study, we compared genic and intergenic regions of three Dioscorea chloroplast genomes and found that the density of SNPs and indels in intergenic sites was about twice and seven times higher than that of SNPs and indels in the genic regions, respectively. A total of 52 primer pairs covering highly variable regions were designed and seven pairs of primers had 80%–100% PCR success rate. PCR amplicons of 73 Dioscorea individuals and assembled sequences of 47 Dioscorea SRAs were used for estimating intraspecific and interspecific divergence for the seven loci: The rpoB‐trnC locus had the highest interspecific divergence. Automatic barcoding gap discovery (ABGD), Poisson tree processes (PTP), and generalized mixed Yule coalescence (GMYC) analysis were applied for species delimitation based on the seven loci and successfully identified the majority of species, except for species in the Enantiophyllum section. Phylogenetic analysis of 51 Dioscorea individuals (28 species) showed that most individuals belonging to the same species tended to cluster in the same group. Our results suggest that the variable loci derived from comparative analysis of plastid genome sequences could be good DNA barcode candidates for taxonomic analysis and species delimitation.     

中国葫芦科1新记录种——厚叶棒锤瓜

报道了葫芦科厚叶棒锤瓜[Neoalsomitra sarcophylla(Wall.)Hutch.]在中国的分布新记录。该物种形态上与N.balansae(Gagnep.)Hutch.近似,但蒴果长3~4cm;种子长6~7mm,宽3~4mm,两端呈短角状。本研究利用DNA条形码技术对该物种进行测序,获得matK、rbcL、psbA-trnH 3个基因序列。应用Blast法为该物种的分类处理提供佐证。凭证标本存放于广西药用植物园标本馆(GXMG)。     

Phylogeny of Carpinus and subfamily Coryloideae (Betulaceae) based on chloroplast and nuclear ribosomal sequence data

Phylogenetic studies were conducted for Carpinus and the subfamily Coryloideae (Betulaceae) using sequences of the chloroplast matK gene, the trnL-trnF region (trnL intron, and trnL [UAA] 3' exon-trnF [GAA] intergenic spacer) and the psbA-trnH intergenic spacer, and the nuclear ribosomal ITS regions. The combined analyses of the three chloroplast regions suggest that Coryloideae is monophyletic; Ostryopsis is sister to the Carpinus - Ostrya clade; Corylus is monophyletic and sister to the Ostrya - Carpinus - Ostryopsis clade; Ostrya is paraphyletic; and within Carpinus, species of sect. Carpinus from eastern Asia form a monophyletic group, whereas the positions of C. betulus from Europe and C. caroliniana from eastern North America are unresolved within the Carpinus clade. The cpDNA tree generated in this study is largely congruent with the previously published ITS results, but the ITS tree places Carpinus sect. Distegocarpus as sister to the Ostrya - Carpinus sect. Carpinus clade. Future work is needed to examine the relationships within the Ostrya - Carpinus clade, evaluate the generic status of Ostrya, and test the phylogenetic position of Ostryopsis.     

New insights into DNA barcoding of seagrasses

Taxonomists find some plant genera challenging because of the few morphological differences or unclear characters among closely related species, which leads to the misidentification of taxa. DNA barcoding is an approach to identify species by using short orthologous DNA sequences, known as ‘DNA barcodes’. Concatenated rbcL and matK sequences are considered DNA barcodes for seagrasses. However, these markers are not applicable to all members of seagrasses at the species level, especially within the genus Halophila. Our previous studies indicated that the internal transcribed spacer (ITS) showed higher species resolution than the concatenated rbcL and matK sequences in the case of Halophila ovalis and closely related species. In this study, 26 ITS, two rbcL and two matK consensus sequences from 18 seagrass taxa belonging to four families collected in India, Vietnam, Germany, Croatia and Egypt were processed. Molecular ITS analysis resolved five clades. The results also indicate that the Cymodoceaceae family might be a non-monophyletic group. In conclusion, ITS could be applied as a DNA barcode for seagrasses instead of the rbcL/matK system previously proposed.     

Molecular phylogeny of Juglans (Juglandaceae): a biogeographic perspective

The eastern Asian and eastern North American disjunction in Juglans offers an opportunity to estimate the time since divergence of the Eurasian and American lineages and to compare it with paleobotanical evidence. Five chloroplast DNA noncoding spacer (NCS) sequences: trnT−trnF, psbA−trnH, atpB−rbcL, trnV-16S rRNA, and trnS-trnfM and data from earlier studies (matK, ITS, and nuclear RFLP) were used to reconstruct phylogeny and to estimate the divergence time of major lineages. Seventeen taxa from four sections of Juglans and two outgroup taxa, Pterocarya stenoptera and Carya illinoiensis were included. NCS data was congruent only with matK data. Both maximum parsimony (MP) and maximum likelihood (ML) cladograms were concordant at the sectional level and revealed three well-supported monophyletic clades corresponding to sections Juglans, Cardiocaryon, and Rhysocaryon in both NCS and combined analyses. The single extant American butternut, Juglans cinerea was placed within the poorly resolved, but well-supported Rhysocaryon. Placement of taxa within Rhysocaryon and Cardiocaryon were inconsistent between NCS and combined analyses. Overall, the results suggest that: (1) the NCS sequence divergence observed within and between sections of Juglans is low and the addition of matK data only marginally improved resolution within Rhysocaryon; (2) the early divergence of section Juglans in both MP and ML analyses of NCS and combined data implies its ancient origin in contrast to fossil evidence, which suggests the earliest divergence of sections Rhysocaryon and Cardiocaryon; and (3) the extant taxa may not hold the footprints to unravel the evolutionary history of the genus.     

Application of DNA Barcodes in Asian Tropical Trees – A Case Study from Xishuangbanna Nature Reserve,Southwest China

Background

Within a regional floristic context, DNA barcoding is more useful to manage plant diversity inventories on a large scale and develop valuable conservation strategies. However, there are no DNA barcode studies from tropical areas of China, which represents one of the biodiversity hotspots around the world.

Methodology and Principal Findings

A DNA barcoding database of an Asian tropical trees with high diversity was established at Xishuangbanna Nature Reserve, Yunnan, southwest China using rbcL and matK as standard barcodes, as well as trnH–psbA and ITS as supplementary barcodes. The performance of tree species identification success was assessed using 2,052 accessions from four plots belonging to two vegetation types in the region by three methods: Neighbor-Joining, Maximum-Likelihood and BLAST. We corrected morphological field identification errors (9.6%) for the three plots using rbcL and matK based on Neighbor-Joining tree. The best barcode region for PCR and sequencing was rbcL (97.6%, 90.8%), followed by trnH–psbA (93.6%, 85.6%), while matK and ITS obtained relative low PCR and sequencing success rates. However, ITS performed best for both species (44.6–58.1%) and genus (72.8–76.2%) identification. With trnH–psbA slightly less effective for species identification. The two standard barcode rbcL and matK gave poor results for species identification (24.7–28.5% and 31.6–35.3%). Compared with other studies from comparable tropical forests (e.g. Cameroon, the Amazon and India), the overall performance of the four barcodes for species identification was lower for the Xishuangbanna Nature Reserve, possibly because of species/genus ratios and species composition between these tropical areas.

Conclusions/Significance

Although the core barcodes rbcL and matK were not suitable for species identification of tropical trees from Xishuangbanna Nature Reserve, they could still help with identification at the family and genus level. Considering the relative sequence recovery and the species identification performance, we recommend the use of trnH–psbA and ITS in combination as the preferred barcodes for tropical tree species identification in China.     

Chloroplast DNA phylogeography of Fagus crenata (Fagaceae) in Japan

 CpDNA variation in Japanese beech, Fagus crenata (Fagaceae), was studied in 45 populations distributed throughout the species' range. Two cpDNA regions were sequenced: the non-coding region between the trnL (UAA) 5′exon and trnF (GAA), and the trnK region (including matK). Thirteen distinct cpDNA haplotypes were recognized and each haplotype was found to be geographically structured. Two major clades (I and II+III) were revealed in phylogenetic analyses among the haplotypes using F. sylvatica as an outgroup. The haplotypes of Clade I were distributed mainly along the Japan Sea side of the Japanese Archipelago, while those of Clade II+III occurred chiefly along the Pacific Ocean side. Consequently, the distribution of the two major cpDNA clades suggests that there were two migration routes in the history of F. crenata; one along the Japan Sea and the other along the Pacific Ocean side of the Japanese Islands. Received March 19, 2001 Accepted November 22, 2001     

Phylogeography and population genetics of Acanthophyllum squarrosum complex (Caryophyllaceae) in the Irano-Turanian region

Acanthophyllum squarrosum and two closely related species, A. heratense and A. laxiusculum (Caryophyllaceae), form a complex that covers parts of subalpine steppes of the Irano-Turanian (IT) region. In this study, we explored the genetic structure and phylogeography of this complex based on partial sequences of two chloroplasts (psbA–trnH and rpl32–trnL (UAG)) and two nuclear (EST24 and nrITS) DNA regions. We analysed 80 individuals from eight populations and detected 12 chloroplast haplotypes, 16 and eight nuclear alleles in EST24 and nrITS sequences, respectively. Phylogenetic trees and haplotype networks did not show distinct genetic groups in the complex and this could be explained by incomplete lineage sorting or introgression between species. Divergence time analysis revealed a Quaternary origin for A. squarrosum complex at approximately 1.8 million years ago (Mya) and the neutrality test results indicated that this complex experienced a recent population expansion. AMOVA analysis of the chloroplast regions showed a significant genetic differentiation among populations and low genetic differentiation within populations, but opposite results were found with nuclear markers, implying introgression between A. squarrosum complex populations.     

Promise and Challenge of DNA Barcoding in Venus Slipper (Paphiopedilum)

Orchidaceae are one of the largest families of flowering plants, with over 27,000 species described and all orchids are listed in CITES. Moreover, the seedlings of orchid species from the same genus are similar. The objective of DNA barcoding is rapid, accurate, and automated species identification, which may be used to identify illegally traded endangered species from vegetative specimens of Paphiopedilum (Venus slipper), a flagship group for plant conservation with high ornamental and commercial values. Here, we selected eight chloroplast barcodes and nrITS to evaluate their suitability in Venus slippers. The results indicate that all tested barcodes had no barcoding gap and the core plant barcodes showed low resolution for the identification of Venus slippers (18.86%). Of the single-locus barcodes, nrITS is the most efficient for the species identification of the genus (52.27%), whereas matK + atpF-atpH is the most efficient multi-locus combination (28.97%). Therefore, we recommend the combination of matK + atpF-atpH + ITS as a barcode for Venus slippers. Furthermore, there is an upper limit of resolution of the candidate barcodes, and only half of the taxa with multiple samples were identified successfully. The low efficiency of these candidate barcodes in Venus slippers may be caused by relatively recent speciation, the upper limit of the barcodes, and/or the sampling density. Although the discriminatory power is relatively low, DNA barcoding may be a promising tool to identify species involved in illegal trade, which has broad applications and is valuable for orchid conservation.     

Phylogenetics and taxonomic delimitation of the genus Guizotia (Asteraceae) based on sequences derived from various chloroplast DNA regions

Parsimony-based phylogenetic analyses of the genus Guizotia were undertaken based on DNA sequence data from the following chloroplast DNA (cpDNA) regions: trnT-trnL, trnL-trnF, trnY-rpoB, trnC-petN, psbM-trnD and rps16-trnQ intergenic spacers, trnL, rps16 and matK-5′trnK introns and matK gene. Out of the 26 primers used in this study, 14 were newly designed. The study was conducted to determine (1) the closest relative of Guizotia abyssinica, (2) the taxonomic status of some Guizotia taxa and (3) the subtribal placement of Guizotia in the tribe Heliantheae. The analyses of the sequence data showed that G. abyssinica, G. scabra ssp. scabra, G. scabra ssp. schimperi and G. villosa are phylogenetically closely related. However, G. scabra ssp. schimperi appeared as the most closely related taxon to G. abyssinica. Based on this phylogenetic analysis, we suggest that the two subspecies of G. scabra are better treated as separate species. The analysis also clearly demonstrated that “Chelelu” and “Ketcha” are distinct Guizotia species. The trnT-trnL and trnL-trnF intergenic spacer-based phylogenetic analysis of various subtribes of the tribe Heliantheae strongly supports the placement of the genus Guizotia within the subtribe Milleriinae.     

A reassessment of tribal affinities of Cratystylis and Haegiela (Asteraceae) based on three chloroplast DNA sequences

 The tribal affinities of Cratystylis and Haegiela were assessed using three chloroplast DNA sequences, the trnL/F spacer, the trnL intron and the matK coding region. The outgroup was represented by two species of the subfamily Barnadesioideae, whereas one to seven genera (45 species including Cratystylis and Haegiela) of the tribes of the Asteroideae [Anthemideae (6 genera), Astereae (7), Calenduleae (2), Gnaphalieae (7), Heliantheae s.l. (5), Inuleae s.str. (3), Plucheeae (3), Senecioneae (4)] and Cichorioideae, [Arctotideae (1), Cardueae (2), Lactuceae (2), Liabeae (1), Mutisieae (1) and Vernonieae (1)] comprise the ingroup. Phylogenetic analysis indicates that Cratystylis has strong support as a member of the tribe Plucheeae, whereas Haegiela is a member of Gnaphalieae. At some point in their taxonomic history, both genera have been placed in these tribes and there are good morphological and chemical characters that justify these placements. The monotypic Haegiela was once included in Epaltes (Plucheeae) and this study has confirmed the need for the separation of the two genera. The genus Cratystylis appears to be monophyletic. Received August 26, 2002; accepted September 19, 2002 Published online: February 7, 2003