Phylogenetic analysis of the internal transcribed spacer region of Japanese Lilium species

 The DNA from 16 Lilium species and one variety endemic to or naturalized in Japan were obtained and their internal transcribed spacer regions of nuclear ribosomal DNA (nrDNA) were amplified by PCR and sequenced by cycle sequencing. Phylogenetic analysis of the ITS sequences supported the validity of Comber’s classification system. It has also provided molecular evidence for the transfer of Lilium dauricum to sect. Sinomartagon. The phylogenetic relationships revealed by ITS DNA analysis were supported by previously published crossability data. The molecular phylogeny of Japanese Lilium species was discussed with reference to the putative migration routes of these species. Received: 30 June 1998 / Accepted: 19 October 1998     

Phylogeny of Allium L. subg. Melanocrommyum (Webb et Berth.) Rouy based on DNA sequence analysis of the internal transcribed spacer region of nrDNA

 Sequence analysis of the ITS region of nuclear ribosomal DNA from subgeneric representatives of Allium L. produced phylogenetic trees which concurred with previous conclusions based on classical taxonomy. Phylogenetic analysis revealed a closer relationship between Nectaroscordum siculum and Allium cernuum (representing Amerallium) than between A. cernuum and the rest of the Allium species employed in this study. The phylogeny of subg. Melanocrommyum based on ITS sequences largely agreed with inferences made by previous researchers based on morphology or a restriction analysis of chloroplast DNA. However, the phylogenetic positions of Allium protensum and Allium macleanii based on ITS sequences did not correspond to their morphological similarity with Allium schubertii and Allium giganteum, respectively. Received: 15 February 1998 / Accepted: 12 March 1998     

The taxonomic relationships within the genus Excoecaria L. (Euphorbiaceae) based on leaf morphology and rDNA sequence data

The tropical Indo-Pacific genus Excoecaria L. (Euphorbiaceae) has several closely related species in Australia whose taxonomic relationships are unclear. The most widely reported species in Australia is the mangrove species Excoecaria agallocha L. (type species), whose taxonomic and geographic limits are difficult to define from its closely related species or sub-species. Two additional taxa have also been described but not clearly differentiated from the type species: Excoecaria dallachyana Baillon and Excoecaria ovalis Endl. This project aimed to determine the taxonomic relationships of the Australian Excoecaria species using both leaf morphological data and DNA sequence data from the internal transcribed spacer (ITS) region of ribosomal genes. The nucleotide differences in the examined ITS1 region show that E. agallocha from eastern Australia and E. ovalis from Western Australia respectively, are genetically uniform within species but differ from each other consistently, thus supporting species status. The leaf morphological data also support this view: single factor analysis of variance consistently separated E. ovalis from E. agallocha on the basis of leaf width, leaf length and length of petiole. In contrast, E. ovalis from the Gulf of Carpentaria differs only slightly from E. ovalis in Western Australia, but no evidence was found to suggest any leaf morphological differentiation within this species. The analysis also suggests that E. dallachyana is not closely related to either mangrove species E. agallocha or E. ovalis, despite superficial morphological similarities.     

A PCR based SNPs marker for specific characterization of English walnut (<Emphasis Type="Italic">Juglans regia</Emphasis> L.) cultivars

English walnut (Juglans regia L.) is the most economically important species from all the 21 species belonging to the genus Juglans and is an important and healthy food as well as base material for timber industry. The aim of this study was to develop a simple technique for specific characterization of English walnut using DNA method. The first and second internal transcribed spacers (ITS1 and ITS2) as well as the intervening 5.8S coding region of the rRNA gene for 18 cultivars of J. regia L. isolated from different geographic origins were characterized. The size of the spacers sequences ranged from 257 to 263 bases for ITS1 and from 217 to 219 bases for ITS2. Variation of GC contents has also been observed and scored as 55–56.7 and 57.1–58.9% for ITS1 and ITS2, respectively. This data exhibited the presence of polymorphism among cultivars. Alignment of the ITS1-5.8S-ITS2 sequences from 18 walnut cultivars showed that there were 244 single nucleotide polymorphisms (SNPs) and 1 short insertion–deletion (indel) at 5′ end ITS1. Amplification refractory mutation system strategy was successfully applied to the SNP markers of the ITS1 and ITS2 sequences for the fingerprinting analysis of 17 on 18 walnut cultivars. The prediction of ITS1 and ITS2 RNA secondary structure from each cultivar was improved by detecting key functional elements shared by all sequences in the alignments. Phylogenetic analysis of the ITS1-5.8S-ITS2 region clearly separated the isolated sequences into two clusters. The results showed that ITS1 and ITS2 region could be used to discriminate these walnut cultivars.     

A Molecular Phylogeny of Lilium in the Internal Transcribed Spacer Region of Nuclear Ribosomal DNA

Phylogenetic relationships among 55 species of Lilium, Cardiocrinum giganteum, and Nomocharis saluenensis were inferred from nucleotide sequence variations in the internal transcribed spacer (ITS) regions of 18S–25S nuclear ribosomal DNA. The phylogeny derived from ITS sequences estimated using maximum-likelihood methods indicated that (1) most of the species construct their own clade according to the classification based on morphological features at the section level; (2) section Daurolirion is not independent of Sinomartagon, and it is appropriate to integrate two sections as Sinomartagon; (3) it is appropriate that L. henryi and L. bulbiferum are classified into subsection 6a and Sinomartagon–Daurolirion, respectively; (4) subsection 6b is much closer to Sinomartagon than subsection 6a and Archelirion, and it arose directly from Sinomartagon; and (5) Lilium is much closer to Nomocharis than Cardiocrinum. Phylogenetic estimation using sequences of the ITS region is suitable at the levels of genus, section, and most of subsection. Received: 18 December 1998 / Accepted: 14 March 1999     

DNA barcoding of Pedicularis L. (Orobanchaceae): Evaluating four universal barcode loci in a large and hemiparasitic genus

Abstract One application of DNA barcoding is species identification based on sequences of a short and standardized DNA region. In plants, various DNA regions, alone or in combination, have been proposed and investigated, but consensus on a universal plant barcode remains elusive. In this study, we tested the utility of four candidate barcoding regions (rbcL, matK, trnH‐psbA, and internal transcribed spacer (ITS)) as DNA barcodes for discriminating species in a large and hemiparasitic genus Pedicularis (Orobanchaceae). Amplification and sequencing was successful using single primer pairs for rbcL, trnH‐psbA, and ITS, whereas two primer pairs were required for matK. Patterns of sequence divergence commonly showed a “barcoding gap”, that is, a bimodal frequency distribution of pairwise distances representing genetic diversity within and between species, respectively. Considering primer universality, ease of amplification and sequencing, and performance in discriminating species, we found the most effective single‐region barcode for Pedicularis to be ITS, and the most effective two‐region barcode to be rbcL + ITS. Both discriminated at least 78% of the 88 species and correctly identified at least 89% of the sequences in our sample, and were effective in placing unidentified samples in known species groups. Our results suggest that DNA barcoding has the potential to aid taxonomic research in Pedicularis, a species‐rich cosmopolitan clade much in need of revision, as well as ecological studies in its center of diversity, the Hengduan Mountains region of China.     

Genetic differentiation of Colletotrichum gloeosporioides and C. truncatum associated with Anthracnose disease of papaya (Carica papaya L.) and bell pepper (Capsium annuum L.) based on ITS PCR-RFLP fingerprinting

Members of the genus Colletotrichum include some of the most economically important fungal pathogens in the world. Accurate diagnosis is critical to devising disease management strategies. Two species, Colletotrichum gloeosporioides and C. truncatum, are responsible for anthracnose disease in papaya (Carica papaya L.) and bell pepper (Capsicum annuum L.) in Trinidad. The ITS1–5.8S–ITS2 region of 48 Colletotrichum isolates was sequenced, and the ITS PCR products were analyzed by PCR-RFLP analysis. Restriction site polymorphisms generated from 11 restriction enzymes enabled the identification of specific enzymes that were successful in distinguishing between C. gloeosporioides and C. truncatum isolates. Species-specific restriction fragment length polymorphisms generated by the enzymes AluI, HaeIII, PvuII, RsaI, and Sau3A were used to consistently resolve C. gloeosporioides and C. truncatum isolates from papaya. AluI, ApaI, PvuII, RsaI, and SmaI reliably separated isolates of C. gloeosporioides and C. truncatum from bell pepper. PvuII, RsaI, and Sau3A were also capable of distinguishing among the C. gloeosporioides isolates from papaya based on the different restriction patterns that were obtained as a result of intra-specific variation in restriction enzyme recognition sites in the ITS1–5.8S–ITS2 rDNA region. Of all the isolates tested, C. gloeosporioides from papaya also had the highest number of PCR-RFLP haplotypes. Cluster analysis of sequence and PCR-RFLP data demonstrated that all C. gloeosporioides and C. truncatum isolates clustered separately into species-specific clades regardless of host species. Phylograms also revealed consistent topologies which suggested that the genetic distances for PCR-RFLP-generated data were comparable to that of ITS sequence data. ITS PCR-RFLP fingerprinting is a rapid and reliable method to identify and differentiate between Colletotrichum species.     

Phylogenetic relationships among species of <Emphasis Type="Italic">Leotia</Emphasis> (Leotiales) based on ITS and RPB2 sequences

Thirty-three collections of Leotia were used to investigate inter-and infra-specific relationships in the genus. Collections were obtained from various parts of the world and represent the ascomatal color forms typical in species of the genus. The ITS rDNA and a variable region of the RNA polymerase II (RPB2) gene were sequenced and analyzed using parsimony and maximum likelihood methods. Although ITS and RPB2 tree topologies differed in regard to the position of two clades of L. lubrica and L. atrovirens, no significant conflict between ITS and RPB2 data or trees was found as determined by the partition homogeneity test. RPB2 sequences in general gave results comparable to ITS; the RPB2 sequences were more easily aligned. Phylogenetic analysis of the sequence data indicates that L. viscosa, L. lubrica and L. atrovirens are polyphyletic species. This suggests that ascomatal color in fresh specimens is not a reliable character alone for determining species in this group. Four major well-supported groups were found; these do not fully correspond to the commonly recognized species. Stipe color, in both fresh and dry condition, seems to correlate with the major recognized groups but features of the ascospores, asci and paraphyses prove too variable to be informative. The most basal group of Leotia species, identified as L. atrovirens, differ from all others by having stipes without gel tissue in their outer layers.     

RFLP marker analysis supports tetrasonic inheritance in Lotus corniculatus L.

Lotus corniculatus is a tetraploid (2n=4x=24) perennial forage legume and has been reported to have tetrasomic inheritance for several traits, although it has also been reported to show disomic inheritance. Molecular markers were used to clarify whether tetrasomic inheritance, disomic inheritance, or a combination of both, was found within an F₂ population arising from a cross between two diverse L. corniculatus accessions. The inheritance of ”tetra-allelic” RFLP markers (markers with four segregating bands) indicated that disomic inheritance could not account for the phenotypic F₂ classes observed, and that only tetrasomic inheritance would explain the observed results. Goodness of fit tests for ”tetra-allelic” and ”tri-allelic” (three segregating bands) RFLP marker data suggested support for chromosomal-type tetrasomic inheritance. RFLP genotypes interpreted from autoradiographic signal intensity provided additional support for tetrasomic inheritance and the occurrence of preferential pairing between parental chromosomes. Bivalent pairing was predominant in the two parental lines and their F₁ hybrid in cytological analyses. L. corniculatus has been classified as both an autotetraploid and an allotetraploid species. RFLP evidence of tetrasomic inheritance gives support for L. corniculatus being classified as an autotetraploid species. Even though bivalent pairing occurs, as seen in other autotetraploid species, pairing between any of the four homologous chromosomes is possible. Preferential pairing in the F₁ hybrid suggests that genome differentiation appears to be minimal between homologs within an accession, while genome differentiation is greater between homologs from different accessions of this genetically diverse species. Received: 16 November 1999 / Accepted: 14 July 2000     

Phylogenetic relationships of coffee-tree species (Coffea L.) as inferred from ITS sequences of nuclear ribosomal DNA

 Phylogenetic relationships of Coffea species were estimated from the sequences of the internal transcribed spacer (ITS 2) region of nuclear ribosomal DNA. The ITS 2 region of 37 accessions belonging to 26 Coffea taxa and to three Psilanthus species was directly sequenced from polymerase chain reaction (PCR)-amplified DNA fragments. The level of variation was high enough to make the ITS 2 a useful tool for phylogenetic reconstruction. However, an unusual level of intraspecific variation was observed leading to some difficulty in interpreting rDNA sequence divergences. Sequences were analysed using Wagner parsimony as well as the neighbour-joining distance method. Coffea taxa were divided into several major groups which present a strong geographical correspondence (i.e. Madagascar, East Africa, Central Africa and West Africa). This organisation is well supported by cytogenetic evidence. On the other hand, the results were in contradiction with the present classification of coffee-tree taxa into two genera, namely Coffea and Psilanthus. Furthermore, additivity of parental rDNA types was not observed in the allotetraploid species C. arabica. Received: 25 July 1996 / Accepted: 18 October 1996     

Internal transcribed spacer region evolution in Larix and Pseudotsuga (Pinaceae)

The nuclear ribosomal DNA (nrDNA) internal transcribed spacer (ITS) region has been characterized in the sister genera Larix and Pseudotsuga (Pinaceae). Complete sequences were obtained for seven species of Larix from North America and Eurasia and five species of Pseudotsuga from western North America and eastern Asia. ITS region lengths ranged from 1759 to 1770 bp in Larix and from 1564 to 1571 bp in Pseudotsuga. In both genera, ITS1 is three times as long as the 5.8S plus ITS2 and contains subrepeats as observed in other genera of Pinaceae. Secondary structure models predicted that the subrepeats fold into terminal stem and loop domains. ITS polymorphism detected within individuals of Larix and Pseudotsuga suggests a slow rate of concerted evolution among nrDNA loci. Except for the placement of L. sibirica, phylogenetic analyses of the ITS region agreed with previously reported restriction site analyses of Larix and Pseudotsuga. The data were not consistent with phylogenetic hypotheses for Larix based primarily upon ovulate cone characters, failing to support a derivation of the North American L. laricina from a short-bracted Eurasian lineage. The phylogenetic hypothesis did not conflict with a stepping stone model of evolution for Pseudotsuga, but a basal lineage could not be inferred for either genus.     

Two Erysiphe species associated with recent outbreak of soybean powdery mildew: results of molecular phylogenetic analysis based on nuclear rDNA sequences

 Serious outbreaks of powdery mildew by a fungus belonging to the mitosporic genus Oidium subgenus Pseudoidium have been reported on soybean (Glycine max) in a wide area of eastern Asia since 1998. The taxonomic and phylogenetic placement of the causal fungus has not yet been determined because of lack of the perfect stage. We found ascomata having mycelioid appendages on a single leaf of soybean infested by powdery mildew. Molecular phylogenetic analysis was conducted based on a total of 14 sequences of the rDNA internal transcribed spacer (ITS) region from 13 soybean and wild soybean (Glycine soja) materials collected in Japan, Korea, Vietnam, and the United States, combined with 47 sequence data obtained from the DNA databases. It was revealed that two Erysiphe species were associated with the outbreak of soybean powdery mildew. There was 16% difference between the two species in genetic divergence of the ITS sequence. One species with perfect stage has an ITS sequence identical to that of Erysiphe glycines on Amphicarpaea and is identified as Erysiphe glycines based on the ITS sequence and morphology of ascomata. The second species, without the perfect stage, is likely to be Erysiphe diffusa (= Microsphaera diffusa), known as the fungus causing soybean powdery mildew in the United States, because the ITS sequences are identical to those from materials collected in the United States. However, we need materials having ascomata of E. diffusa to confirm the species name. Received: March 15, 2002 / Accepted: May 22, 2002     

A technique for identifying the roots of different species in mixed samples using nuclear ribosomal DNA

Abstract. Question: Is it possible to determine the species composition of root samples containing multiple species, without first disentangling individual roots? Methods: The internal transcribed spacer (ITS) region of nuclear ribosomal DNA was amplified and sequenced from four California annual grassland species (two Poaceae and two Asteraceae). Restriction enzymes that cut the ITS region of each species into uniquely sized fragments were identified based on DNA sequence variation of the ITS regions. Mixed root samples were analysed to test the ability of the method to identify the presence or absence of each species in multi‐species samples. Results: The technique successfully identified species present in multi‐species samples. ITS regions were shorter in Poaceae than in Asteraceae, so size differences alone were sufficient to distinguish these taxonomic groups. At the species level, digestion of ITS regions with the appropriate restriction enzymes yielded at least one uniquely sized fragment for each species. Conclusions: This method is the first to identify the species composition of mixed root samples. It should be applicable to most plant species because the ITS region is flanked by universal primers and most species have unique ITS sequences. The ability to determine species‐specific rooting distributions has broad applications in vegetation science.     

Taxonomy and phylogeny of Leptopilina species (Hymenoptera: Cynipoidea: Figitidae) attacking frugivorous drosophilid flies in Japan,with description of three new species

Despite the intensive use of the Leptopilina genus and its drosophilid hosts as model systems in the study of host–parasitoid interactions, the diversity and distribution of the species occurring in the Asian region remain elusive. Here we report the phylogeny of Japanese Leptopilina species attacking frugivorous drosophilid flies, based on COI, ITS1 and ITS2 sequences. Consistent with molecular data, hybridization experiments and morphological examination, five species were recorded in Japan: Leptopilina heterotoma, L. victoriae and three new species, two occurring in the Ryukyu archipelago, L. ryukyuensis and L. pacifica, and another species, L. japonica, distributed in Honshu and Hokkaido. Leptopilina japonica is further divided into two subspecies, L. j. japonica occurring in Japan, and L. j. formosana occurring in Taiwan. According to these results, we discuss the evolution, speciation and colonization history of Japanese Leptopilina species.     

Do floral odor profiles geographically vary with the degree of specificity for pollinators? Investigation in two sapromyophilous Arum species (Araceae)

We compared floral odour profiles among populations of two Arum species which show different degrees of specificity for their fly pollinators. Insects were collected from inflorescences in four populations of Arum italicum and two populations of Arum maculatum. In six Arum populations, we compared inflorescences odour profiles collected by Solid Phase Micro Extraction (SPME) and analysed by gas chromatography. We confirmed that from a pollination point of view, A. italicum is an opportunist species, as it is mainly pollinated by insects of the families Psychodidae, Chironomidae and Sciaridae, whereas A. maculatum is a specialist species, as it is 90% pollinated by Psychodidae. In all populations, Arum italicum was less attractive to pollinators than Arum maculatum. Floral odour profiles of A. italicum were not geographically structured among populations, suggesting a high gene flow or adaptation to a fluctuant guild of pollinators. On the contrary, odour profiles of A. maculatum varied between the two populations studied suggesting a lower gene flow or adaptation to different local pollinator preferences     

Molecular characterisation and origin of the Coffea arabica L. genome

Restriction fragment length polymorphism (RFLP) markers were used in combination with genomic in situ hybridisation (GISH) to investigate the origin of the allotetraploid species Coffea arabica (2n = 44). By comparing the RFLP patterns of potential diploid progenitor species with those of C. arabica, the sources of the two sets of chromosomes, or genomes, combined in C. arabica were identified. The genome organisation of C. arabica was confirmed by GISH using simultaneously labelled total genomic DNA from the two putative genome donor species as probes. These results clearly suggest that C. arabica is an amphidiploid formed by hybridisation between C. eugenioides and C. canephora, or ecotypes related to these diploid species. Our results also indicate low divergence between the two constituent genomes of C. arabica and those of its progenitor species, suggesting that the speciation of C. arabica took place relatively recently. Precise localisation in Central Africa of the site of the speciation of C. arabica, based on the present distribution of the coffee species, appears difficult, since the constitution and extent of tropical forest has varied considerably during the late Quaternary period. Received: 6 June 1998 / Accepted: 10 November 1998     

Examining morphological and molecular diagnostic character states of Cichorium intybus L. (Asteraceae) and C. spinosum L.

The inter- and intraspecific variability of Cichorium intybus L. was examined to evaluate potential morphological and molecular diagnostic character states. Two diagnostic and one overlapping morphological character clearly delimit the two species C. intybus and C. spinosum. All applied molecular methods (ITS, AFLP, Microsatellites) failed to significantly discriminate between these accepted species. As the morphological traits are fixed and heritable, criteria for species delimitation are fulfilled. These traits, however, are apparently due to mutations of a few crucial loci affecting the morphological diagnostic character states. Intraspecific variability within C. intybus revealed to be highly influenced by plastic response to local environmental factors and subspecific delimitation cannot be supported.     

Hybridization between the threatened plant, <Emphasis Type="Italic">Lespedeza leptostachya</Emphasis> Englem. and its co-occurring congener <Emphasis Type="Italic">Lespedeza capitata</Emphasis> Michx.: morphological and molecular evidence

Natural hybridization is common in the genus Lespedeza. No hybrids between Lespedeza leptostachya Englem. and Lespedeza capitata Michx. are formally recognized in any of the current floras, however observations in the field suggest that hybridization might occur in many of their shared habitats. Putative hybrids were compared to L. leptostachya and L. capitata using morphological measurements and screened for the presence of species-specific trnL-F gene region (cpDNA) and the ITS gene region (nrDNA). A discriminate analysis of 10 morphological measurements identified the hybrids as intermediate to both parents with two PCA axes explaining 99% of the variation between taxa. The presence of hybrids was confirmed by genetic markers with individuals morphologically identified as hybrids having cpDNA trnL-F genotypes identical to L. leptostachya and the ITS (nrDNA) phenotypes in most cases contain the ITS genotype of both parents, however, some putative hybrid individuals contained the ITS genotype of only one parents. Those individuals with L. leptostachya ITS and trnL-F could be a case of misclassification, but the presence of both L. capitata ITS genotypes and L. leptostachya trnL-F genotypes suggest segregation has occurred, which may result from either selfing or backcrossing.     

Molecular phylogenetics of Hippophae L. (Elaeagnaceae) based on the internal transcribed spacer (ITS) sequences of nrDNA

 The genus Hippophae comprises 7 species and 8 subspecies according to the latest classification, and has shown enormous ecological, nutrient and medicinal values. Here we analyzed the phylogenetic relationships among 15 taxa of the genus by comparing sequences of the internal transcribed spacer (ITS) region of nuclear ribosomal DNA (nrDNA). ITS sequences in Hippophae varied in length from 651 bp to 666 bp. The aligned sequences were 690 bp in length and 269 (39.0%) were variable sites with 150 being parsimony-informative. The amount of polymorphism observed within a taxon was extremely low in most taxa except for two putative hybrid species. The aligned sequences were analyzed by maximum parsimony (MP) and neighbor-joining (NJ) methods. In the strict consensus trees of parsimony analysis, the monophyly of Hippophae was supported by 100% bootstrap value. H. tibetana was at the basal position of the genus, and the remaining taxa formed two clades with high bootstrap support. The first clade included subspecies of H.␣rhamnoides and the other one consisted of remaining species. Parsimony analysis also suggested that the species H. tibetana, H. neurocarpa and H.␣salicifolia were all distinct. Although the sequence divergence among subspecies of H. rhamnoides was also remarkably high, the molecular data supported the monophyly of H. rhamnoides when H. rhamnoides subsp. gyantsensis Rousi was excxluded. The NJ trees showed essentially the same topology. The taxonomical arrangement that divided the genus into two sections was not supported based on the ITS sequences. However, the hybrid origin of H. goniocarpa and H. litangensis proposed previously was supported by the present ITS data. Received January 7, 2002; accepted May 10, 2002 Published online: November 22, 2002 Addresses of the authors: Kun Sun, Xuelin Chen, Ruijun Ma, Qin Wang, Institute of Botany, Northwest Normal University, Lanzhou 730070, China. Changbao Li, Song Ge (e-mail: gesong@ns.ibcas.ac.cn or song_ge@hotmail.com), Laboratory of Systematic and Evolutionary Botany, Institute of Botany, Chinese Academy of Sciences, Beijing 100093, China.     

Phylogeny of korean rubus (rosaceae) based on its (nrDNA) and trnL/F intergenic region (cpdna)

We examined the phylogeny of the genusRubus in Korea using an internal transcribed spacer (ITS) of the nuclear ribosomal DNA and a trnL-trnF (trnL/F) intergenic region of the chloroplast DNA. In all, 21 ingroup species (1.2 kb for each species) were analyzed using parsimony, resulting in 672 aligned sequences from ITS, and 502 bases with trnL/F. Individual and combined analysis of ITS and trnL/F data proved that the genusRubus is a monophyletic group. This phylogeny also substantiated a previous sectional classification scheme rather than a subgerius classification scheme. However, our results did not support the earlier sectional classification by Focke (sect.Corchorifolii), but did support the sectional classification of Nakai: sect.Crataegifolii (R. crataegifolius, R. takesimensis andR. trifidus), and sect.Villosii (R. corchorifolius). Most of these species, which are found in Korea and belong to subg.Idaeobatus, appeared in two different groups in all data sets. This suggests that this subgenus is a polyphyletic group that has gone through at least two independent evolutionary processes. The taxa, when mapped onto the combined tree, showed that the occurrence of their morphological characters of simple and compound leaves was concurrent in KoreanRubus. ITS sequence data were consistent overall with the geographical distribution of each species. Furthermore, the trnL/F sequence data provided phylogenetic information within closely related species.