Cryptic repeated genomic recombination during speciation in Gossypium gossypioides

The Mexican cotton Gossypium gossypioides is a perplexing entity, with conflicting morphological, cytogenetic, and molecular evidence of its phylogenetic affinity to other American cottons. We reevaluated the evolutionary history of this enigmatic species using 16.4 kb of DNA sequence. Phylogenetic analyses show that chloroplast DNA (7.3 kb), nuclear ribosomal internal transcribed spacers (ITS; 0.69 kb), and unique nuclear genes (8.4 kb) yield conflicting resolutions for G. gossypioides. Eight low-copy nuclear genes provide a nearly unanimous resolution of G. gossypioides as the basalmost American diploid cotton, whereas cpDNA sequences resolve G. gossypioides deeply nested within the American diploid clade sister to Peruvian G. raimondii, and ITS places G. gossypioides in an African (rather than an American) clade. These data, in conjunction with previous evidence from the repetitive fraction of the genome, implicate a complex history for G. gossypioides possibly involving temporally separated introgression events from genetically divergent cottons that are presently restricted to different hemispheres. Based on repetitive nuclear DNA, it appears that G. gossypioides experienced nuclear introgression from an African species shortly after divergence from the remainder of the American assemblage. More recently, hybridization with a Mexican species may have resulted in cpDNA introgression, and possibly a second round of cryptic nuclear introgression. Gossypium gossypioides provides a striking example of the previously unsuspected chimeric nature of some plant genomes and the resulting phylogenetic complexity produced by multiple historical reticulation events.     

Analysis of SSR dynamics in chloroplast genomes of Brassicaceae family

Simple sequence repeats (SSRs) are present abundantly in most eukaryotic genomes. They affect several cellular processes like chromatin organization, regulation of gene activity, DNA repair, DNA recombination, etc. Though considerable data exists on using nuclear SSRs to infer phylogenetic relationships, the potential of chloroplast microsatellites (cpSSR), in this regard, remains largely unexplored. In the present study we probe various nucleotide repeat motifs (NRMs) / types of SSRs present in chloroplast genomes (cpDNA) of 12 species belonging to Brassicaceae family. NRMs show a non-random distribution in coding and non-coding compartments of cpDNA. As expected, trinucleotide repeats are more common in coding regions while other repeat motifs are prominent in non-coding DNA. Total numbers of SSRs in coding region show little variation between species while considerable variation is exhibited by SSRs in non-coding regions. Finally, we have designed universal primers that yield polymorphic amplicons from all 12 species. Our analysis also suggests that amplicon length polymorphism shows no significant relationship with sequence based phylogeny of SSRs in cpDNA of Brassicaceae family.     

Reticulate evolution in Thuja inferred from multiple gene sequences: implications for the study of biogeographical disjunction between eastern Asia and North America

The eastern Asia-North America disjunction is one of the most interesting biogeographical patterns, but its formation is still in much debate. Here nucleotide sequences of five cpDNA regions, nrDNA ITS and two low-copy nuclear genes (LEAFY, 4CL) were employed to reconstruct the phylogeny and to explore the historical biogeography of Thuja, a typical eastern Asia-North America disjunct genus. High topological discordance was observed between chloroplast and nuclear gene trees, even between different nuclear gene trees, suggesting that Thuja could have a reticulate evolutionary history due to multiple interspecific hybridization events. The eastern Asian species Thuja koraiensis might have obtained its chloroplast genome from the eastern North American species T. occidentalis by chloroplast capture, while the western North American species T. plicata is very likely to have inherited a recombinant cpDNA. Based on the phylogenetic analysis of multiple genes, DIVA-reconstruction of the distribution history, molecular clock estimation and fossil data, we inferred that Thuja could have originated from the high-latitude areas of North America in the Paleocene or earlier with subsequent expansion into eastern Asia through the Bering Land Bridge. The two eastern Asia species T. standishii and T. sutchuenensis have a sister relationship, and their split could have occurred in the Oligocene or early Miocene. In the present study, the selection of molecular markers in biogeographic studies was also discussed. Since most previous studies on the eastern Asia and North America disjunction are based on uniparentally inherited cpDNA and (or) directly sequenced nrDNA ITS data, the historical reticulate evolution in the studied groups might have been underestimated. Therefore, we suggest that multiple genes from different genomes, especially low-copy nuclear genes, be used in this research area in the future.     

Origin and evolution of North American polyploid Silene (Caryophyllaceae)

Nuclear DNA sequences from introns of the low-copy nuclear gene family encoding the second largest subunit of RNA polymerases and the ribosomal internal transcribed spacer (ITS) regions, combined with the psbE-petL spacer and the rps16 intron from the chloroplast genome were used to infer origins and phylogenetic relationships of North American polyploid Silene species and their closest relatives. Although the vast majority of North American Silene species are polyploid, which contrasts to the diploid condition dominating in other parts of the world, the phylogenetic analyses rejected a single origin of the North American polyploids. One lineage consists of tetraploid Silene menziesii and its diploid allies. A second lineage, Physolychnis s.l., consists of Arctic, European, Asian, and South American taxa in addition to the majority of the North American polyploids. The hexaploid S. hookeri is derived from an allopolyploidization between these two lineages. The tetraploid S. nivea does not belong to any of these lineages, but is closely related to the European diploid S. baccifera. The poor resolution within Physolychnis s.l. may be attributed to rapid radiation, recombination among homoeologues, homoplasy, or any combination of these factors. No extant diploid donors could be identified in Physolychnis s.l.     

Phylogenetic reticulation in subtribe Helianthinae

Incongruence between phylogenetic estimates based on nuclear and chloroplast DNA (cpDNA) markers was used to infer that there have been at least two instances of chloroplast transfer, presumably through wide hybridization, in subtribe Helianthinae. One instance involved Simsia dombeyana, which exhibited a cpDNA restriction site phenotype that was markedly divergent from all of the other species of the genus that were surveyed but that matched the restriction site pattern previously reported for South American species of Viguiera. In contrast, analysis of sequence data from the nuclear ribosomal DNA internal transcribed spacer (ITS) region showed Simsia to be entirely monophyletic and placed samples of S. dombeyana as the sister group to the relatively derived S. foetida, a result concordant with morphological information. A sample of a South American species of Viguiera was placed by ITS sequence data as the sister group to a member of V. subg. Amphilepis, which was consistent with cpDNA restriction site data. Samples of Tithonia formed a single monophyletic clade based on ITS sequence data, whereas they were split between two divergent clades based on cpDNA restriction site analysis. The results suggested that cpDNA transfer has occurred between taxa diverged to the level of morphologically distinct genera, and highlight the need for careful and complete assessment of molecular data as a source of phylogenetic information.     

Discordant mtDNA and cpDNA phylogenies indicate geographic speciation and reticulation as driving factors for the diversification of the genus Picea

The phylogeny of the genus Picea was investigated by sequencing three loci from the paternally inherited chloroplast genome (trnK, rbcL and trnTLF) and the intron 2 of the maternally transmitted mitochondrial gene nad1 for 35 species. Significant topological differences were found between the trnK tree and the rbcL and trnTLF phylogenetic trees, and between cpDNA and mtDNA phylogenies. None of the phylogenies matched morphological classifications. The mtDNA phylogeny was geographically more structured than cpDNA phylogenies, reflecting the different inheritance of the two cytoplasmic genomes in the Pinaceae and their differential dispersion by seed only and seed and pollen, respectively. Most North American taxa formed a monophyletic group on the mtDNA tree, with topological patterns suggesting geographic speciation by range fragmentation or by dispersal and isolation. Similar patterns were also found among Asian taxa. Such a trend towards geographic speciation is anticipated in other Pinaceae genera with similar life history, autecology and reproductive system. Incongruences between organelle phylogenies suggested the occurrence of mtDNA capture by invading cpDNA. Incongruences between cpDNA partitions further suggested heterologous recombination presumably also linked to ancient reticulate evolution. Whilst cpDNA appears potentially valuable for molecular taxonomy and systematics purposes, these results emphasize the reduced value of cpDNA to infer vertical descent and the speciation history for plants with paternal transmission and high dispersal of their chloroplast genome.     

Evidence of chloroplast capture in South American Nothofagus (subgenus Nothofagus,Nothofagaceae)

Subgenus Nothofagus, although geographically restricted at present to temperate areas of South America, has captured much attention in discussions of plant biogeography due to its widespread distribution through Gondwanan continents during the Tertiary. However, phylogenetic relationships within the subgenus Nothofagus have not yet been resolved. We examined geographic patterns of intraspecific and interspecific genetic variation to detect whether incongruences in nuclear or plastid DNA phylogenies occur, in order to better understand the evolutionary history of the subgenus Nothofagus. We conducted spatially-explicit sampling at 10 distinct locations throughout the range of austral South American forests and sampled all present Nothofagus species. We used ITS and chloroplast DNA sequences to estimate phylogenetic relationships. A phylogeny constructed from nuclear genes resolved the subgenus Nothofagus as monophyletic. We found that N. antarctica was a sister to a clade of evergreen species (N. betuloides, N. dombeyi, and N. nitida), while N. pumilio likely diverged earlier. Nine cpDNA haplotypes were distinguished in the subgenus Nothofagus which were associated to geographic locations rather than to taxonomic relationships. This species-independent cpDNA phylogeographic structures within the subgenus Nothofagus may be related to repeated chloroplast capture events over geological time in Patagonia.     

Phylogenetic relationships in Gleditsia (Leguminosae) based on ITS sequences

We used nucleotide sequences from the internal transcribed spacers and 5.8S gene of nuclear ribosomal DNA to test competing phylogenetic and biogeographic hypotheses in Gleditsia. Eleven of 13 Gleditsia species were sampled, along with two species of its sister genus, Gymnocladus. Analyses of ITS data and of a combined data set that included sequences of ITS and two chloroplast genes supported several conclusions that were interpreted in light of fossil data and current legume phylogeny. Gleditsia and Gymnocladus appear to have originated in eastern Asia during the Eocene. Eastern North American species of both genera most likely evolved from ancestors that migrated across the Bering land bridge, but the eastern Asian/eastern North American disjunction appears to be much older in Gymnocladus than in Gleditsia. Gleditsia amorphoides, from temperate South America, is sister to the rest of the genus, suggesting early long-distance dispersal from Asia. The remainder of Gleditsia is divided into three unresolved clades, possibly indicating a split early in the evolution of the genus. Two of those clades contain only Asian species, and one contains Asian and North American species. The North American species, Gleditsia triacanthos and Gleditsia aquatica, are polymorphic and paraphyletic with respect to their ITS and cpDNA sequences, which suggests recent diversification.     

Are nuclear loci ideal for barcoding plants? A case study of genetic delimitation of two sister species using multiple loci and multiple intraspecific individuals

Considerable debate remains as to which DNA region should be used to barcode plants. Several different chloroplast (cp) DNA regions (rbcL, matK, and trnH-psbA) and nuclear ribosomal internal transcribed regions (ITS) have been suggested as suitable barcodes in plants. Recently, low-copy nuclear loci were also suggested to be potentially ideal barcode regions. The aim of the present study was to test the effectiveness of these proposed DNA fragments and five additional low-copy loci (CHS, DETl, COPl, PGICl, and RPS2; comprising both coding and non-coding regions) in barcoding closely related species. We examined the divergences within and between two species of Pugioniun (Brassicaceae). We failed to find any interspecific variation from three cpDNA fragments with which to discriminate the two species. However, a single base mutation in the internal transcribed spacer (ITS) could discriminate between the two species consistently. We found more variations among all individuals of the two species using each of the other five low-copy nuclear loci. However, only alleles from one locus (DET1) of the five low-copy loci related to flowering regulations was able to distinguish the sampled individuals into two species. We failed to amplify the corresponding fragments out of Brassicaceae using the designed DETl primers. We further discussed the discrimination power of different loci due to incomplete lineage sorting, gene flow, and species-specific evolution. Our results highlight the possibility of using the nuclear ITS as a core or complementary fragment to barcode recent diverged species.     

Molecular phylogenetics of the clover genus (Trifolium--Leguminosae)

Trifolium, the clover genus, is one of the largest genera of the legume family. We conducted parsimony and Bayesian phylogenetic analyses based on nuclear ribosomal DNA internal transcribed spacer and chloroplast trnL intron sequences obtained from 218 of the ca. 255 species of Trifolium, representatives from 11 genera of the vicioid clade, and an outgroup Lotus. We confirm the monophyly of Trifolium, and propose a new infrageneric classification of the genus based on the phylogenetic results. Incongruence between the nrDNA and cpDNA results suggests five to six cases of apparent hybrid speciation, and identifies the putative progenitors of the allopolyploids T. dubium, a widespread weed, and T. repens, the most commonly cultivated clover species. Character state reconstructions confirm 2n=16 as the ancestral chromosome number in Trifolium, and infer a minimum of 19 instances of aneuploidy and 22 of polyploidy in the genus. The ancestral life history is hypothesized to be annual in subgenus Chronosemium and equivocal in subgenus Trifolium. Transitions between the annual and perennial habit are common. Our results are consistent with a Mediterranean origin of the genus, probably in the Early Miocene. A single origin of all North and South American species is hypothesized, while the species of sub-Saharan Africa may originate from three separate dispersal events.     

Topological incongruence between nuclear and chloroplast DNA trees suggesting hybridization in the urophyllum group of the genus Fagopyrum (Polygonaceae)

We performed phylogenetic analyses of Fagopyrum species in the urophyllum group based on nucleotide sequences of two nuclear genes, FLORICAULA/LEAFY (FLO/LFY) and AGAMOUS (AG), and three segments of chloroplast DNA (cpDNA), rbcL-accD, trnK intron, and trnC-rpoB spacer. The FLO/LFY and AG sequences turned out to be phylogenetically more informative at the intrageneric level than the cpDNA sequences. Congruence among these gene trees, inferred by a maximum-likelihood (ML) method, demonstrated that topologies were partially incongruent between the nuclear and chloroplast DNA phylogenies. The nuclear DNA sequence data supported a monophyletic relation of F. statice, F. gilesii, and F. jinshaense, whereas the former two species formed another monophyletic relation with the F. capillatum-F. gracilipes-F. gracilipedoides-F. rubifolium clade excluding F. jinshaense in the synthetic cpDNA phylogeny. In addition, two divergent sequences of FLO/LFY were found in F. rubifolium (tetraploid). One of these was sister to F. gracilipedoides and another was sister to F. statice, and a monophyletic relation of these two genes was rejected by a bootstrap analysis. These results suggest that hybridization may have occurred during diversification of Fagopyrum species in the urophyllum group, and that F. rubifolium is possibly allotetraploid species.     

Are nuclear loci ideal for barcoding plants? A case study of genetic delimitation of two sister species using multiple loci and multiple intraspecific individuals

Considerable debate remains as to which DNA region should be used to barcode plants. Several different chloroplast (cp) DNA regions (rbcL, matK, and trnH–psbA) and nuclear ribosomal internal transcribed regions (ITS) have been suggested as suitable barcodes in plants. Recently, low-copy nuclear loci were also suggested to be potentially ideal barcode regions. The aim of the present study was to test the effectiveness of these proposed DNA fragments and five additional low-copy loci (CHS, DET1, COP1, PGIC1, and RPS2; comprising both coding and non-coding regions) in barcoding closely related species. We examined the divergences within and between two species of Pugionium (Brassicaceae). We failed to find any interspecific variation from three cpDNA fragments with which to discriminate the two species. However, a single base mutation in the internal transcribed spacer (ITS) could discriminate between the two species consistently. We found more variations among all individuals of the two species using each of the other five low-copy nuclear loci. However, only alleles from one locus (DET1) of the five low-copy loci related to flowering regulations was able to distinguish the sampled individuals into two species. We failed to amplify the corresponding fragments out of Brassicaceae using the designed DET1 primers. We further discussed the discrimination power of different loci due to incomplete lineage sorting, gene flow, and species-specific evolution. Our results highlight the possibility of using the nuclear ITS as a core or complementary fragment to barcode recent diverged species.     

Rapid diversification of the cotton genus (Gossypium: Malvaceae) revealed by analysis of sixteen nuclear and chloroplast genes

Previous molecular phylogenetic studies have failed to resolve the branching order among the major cotton (Gossypium) lineages, and it has been unclear whether this reflects actual history (rapid radiation) or sampling properties of the genes evaluated. In this paper, we reconsider the phylogenetic relationships of diploid cotton genome groups using DNA sequences from 11 single-copy nuclear loci (10?293 base pairs [bp]), nuclear ribosomal DNA (695 bp), and four chloroplast loci (7370 bp). Results from individual loci and combined nuclear and chloroplast DNA partitions reveal that the cotton genome groups radiated in rapid succession following the formation of the genus. Maximum likelihood analysis of nuclear synonymous sites shows that this radiation occurred within a time span equivalent to 17% of the time since the separation of Gossypium from its nearest extant relatives in the genera Kokia and Gossypioides. Chloroplast and nuclear phylogenies differ significantly with respect to resolution of the basal divergence in the genus and to interrelationships among African cottons. This incongruence is due to limited character evolution in cpDNA and either previously unsuspected hybridization or unreliable phylogenetic performance of the cpDNA characters. This study highlights the necessity of using multiple, independent data sets for resolving phylogenetic relationships of rapidly diverged lineages.     

Phylogeny, biogeography, and rates of diversification of New World Astragalus (Leguminosae) with an emphasis on South American radiations

This study uses phylogenetic relationships of New World representatives of the species-rich genus Astragalus (Leguminosae; Papilionoideae) to follow up on recent evidence pointing to rapid and recent plant diversification patterns in the Andes. Bayesian and maximum likelihood phylogenetic analyses were done using nuclear rDNA ITS and chloroplast spacers trnD-trnT and trnfM-trnS1, either separately or in combination. The effect of using partitioned vs. nonpartitioned analyses in a Bayesian approach was evaluated. Highest resolution was obtained when the data were combined in partitioned or nonpartitioned Bayesian analyses. All phylogenies support two clades of South American species nested within the North American species, implying two separate invasions from North to South America. These two clades correspond to the original morphological classification of Johnston (1947 Journal of the Arnold Arboretum 28: 336-409). The mean ages of the South American clades were very recent but still significantly different (1.89 and 0.98 Ma). Upper and lower bounds on rates of diversification varied between 2.01 and 0.65 species/Ma for the older clade and 2.06 and 1.24 species/Ma for the younger clade. Even the lower bounds are still very high, reasserting Neo-Astragalus in the growing list of recent rapid radiations of plants, especially in areas with a high physiographic diversity, such as the Andes.     

Molecular phylogenetics of Houstonia (Rubiaceae): descending aneuploidy and breeding system evolution in the radiation of the lineage across North America

Nuclear and chloroplast DNA sequence variation was used to infer evolutionary relationships within and among members of Houstonia (Rubiaceae) and other closely related genera in North America. Sequences from the internal transcribed spacer (ITS) of the nrDNA and a cpDNA intron in the trnL gene were used to reconstruct phylogenetic relationships of 30 species of Houstonia and closely related genera. The data suggest that the North American species of Houstonia are not monophyletic, but belong to the same lineage as Stenaria. The radiation of this lineage has been accompanied by changes in the basic chromosome number of the major clades through descending aneuploidy. This loss of chromosomes was also associated with northward colonization of North America. However, other characters, such as an annual versus perennial habit, heterostylous and homostylous breeding systems, and the evolution of self fertilization, seem to be labile throughout the lineage, originating multiple times throughout the evolutionary history of the lineage.     

Chloroplast DNA phylogeny, reticulate evolution, and biogeography of Paeonia (Paeoniaceae)

The coding region of the mat K gene and two intergenic spacers, psb A-trn H and trn L(UAA)-trn F(GAA), of cpDNA were sequenced to study phylogenetic relationships of 32 Paeonia species. In the psb A-trn H intergenic spacer, short sequences bordered by long inverted repeats have undergone inversions that are often homoplasious mutations. Insertions/deletions found in the two intergenic spacers, mostly resulting from slipped-strand mispairing, provided relatively reliable phylogenetic information. The mat K coding region, evolving more rapidly than the trnL-trn F spacer and more slowly than the psb A-trn H spacer, produced the best resolved phylogenetic tree. The mat K phylogeny was compared with the phylogeny obtained from sequences of internal transcribed spacers (ITS) of nuclear ribosomal DNA. A refined hypothesis of species phylogeny of section Paeonia was proposed by considering the discordance between the nuclear and cpDNA phylogenies to be results of hybrid speciation followed by inheritance of cpDNA of one parent and fixation of ITS sequences of another parent. The Eurasian and western North American disjunct distribution of the genus may have resulted from interrruption of the continuous distribution of ancestral populations of extant peony species across the Bering land bridge during the Miocene. Pleistocene glaciation may have played an important role in triggering extensive reticulate evolution within section Paeonia and shifting distributional ranges of both parental and hybrid species.     

Evolution and phylogenetic utility of the PHOT gene duplicates in the Verbena complex (Verbenaceae): dramatic intron size variation and footprint of ancestral recombination

A well-resolved species level phylogeny is critically important in studying organismal evolution (e.g., hybridization, polyploidization, adaptive speciation). Lack of appropriate molecular markers that give sufficient resolution to gene trees is one of the major impediments to inferring species level phylogenies. In addition, sampling multiple independent loci is essential to overcome the lineage sorting problem. The availability of nuclear loci has often been a limiting factor in plant species-level phylogenetic studies. Here the two PHOT loci were developed as new sources of nuclear gene trees. The PHOT1 and PHOT2 gene trees of the Verbena complex (Verbenaceae) are well resolved and have good clade support. These gene trees are consistent with each other and previously generated chloroplast and nuclear waxy gene trees in most of the phylogenetic backbone as well as some terminal relationships, but are incongruent in some other relationships. Locus-specific primers were optimized for amplifying and sequencing these two loci in all Lamiales. Comparing intron size in the context of the gene trees shows dramatic variation within the Verbena complex, particularly at the PHOT1 locus. These variations are largely caused by invasions of short transposable elements and frequent long deletions and insertions of unknown causes. In addition, inspection of DNA sequences and phylogenetic analyses unmask a clear footprint of ancestral recombination in one species.     

Mitochondrial DNA variation and reticulate evolution of the genus Abies

The sequences of three regions of mitochondrial DNA (mtDNA) of a total length of 5226 bp were used to study the phylogeography of the genus Abies. The mtDNA haplotype network, comprising 36 studied Abies taxa, consisted of two branches; the first represented all American species plus two Asian, and the second included the remaining Eurasian species. Within these clusters, the haplotypes formed nine major groups, generally corresponding to the clades of the previously obtained phylogeny based on chloroplast DNA (cpDNA), but the relationships of these groups were significantly different; species assignment to the particular mtDNA haplotype group was more in line with its geographical distribution. In addition, the mtDNA haplotype network contains cycles indicating the recombination. It is assumed that the incongruence of cpDNA and mtDNA phylogenies is caused by the introgression capture of alien mtDNA during species hybridization and thus contains information about past migrations. The cases of incongruence of mitochondrial and chloroplast DNA suggesting a migration of Abies between Asia and North America are discussed.     

Utility of Low-Copy Nuclear Gene Sequences in Plant Phylogenetics

Low-copy nuclear genes in plants are a rich source of phylogenetic information. They hold a great potential to improve the robustness of phylogenetic reconstruction at all taxonomic levels, especially where universal markers such as cpDNA and nrDNA are unable to generate strong phylogenetic hypotheses. Low-copy nuclear genes, however, remain underused in plant phylogenetic studies due to practical and theoretical complications in unraveling the evolutionary dynamics of nuclear gene families. The lack of the universal markers or universal PCR primers of low-copy nuclear genes has also hampered their phylogenetic utility. It has recently become clear that low-copy nuclear genes are particularly helpful in resolving close interspecific relationships and in reconstructing allopolyploidization in plants. Gene markers that are widely, if not universally, useful have begun to emerge. Although utilizing low-copy nuclear genes usually requires extra lab work such as designing PCR primers, PCR-cloning, and/or Southern blotting, rapid accumulation of gene sequences in the databases and advances in cloning techniques have continued to make such studies more feasible. With the growing number of theoretical studies devoted to the gene tree and species tree problem, a solid foundation for reconstructing complex plant phylogenies based on multiple gene trees began to build. It is also realized increasingly that fast evolving introns of the low-copy nuclear genes will provide much needed phylogenetic information around the species boundary and allow us to address fundamental questions concerning processes of plant speciation. Phylogenetic and molecular evolutionary analyses of developmentally important genes will add a new dimension to systematic and evolutionary studies of plant diversity.     

Frequent chloroplast capture among Isodon (Lamiaceae) species in Japan revealed by phylogenies based on variation in chloroplast and nuclear DNA

A recent phylogenetic study based only on chloroplast DNA (cpDNA) variation revealed that populations of an Isodon species are frequently embedded paraphyletically among other Isodon species. This phylogenetic discrepancy between species taxonomy and molecular phylogeny was considered to have resulted from chloroplast DNA captures and/or incomplete lineage sorting. To elucidate which of these factors was mainly responsible for the observed phylogenetic pattern, we performed phylogenetic analyses of multiple populations of Isodon species in Japan using cpDNA variation, three single-copy nuclear genes, and double-digest restriction-site-associated DNA sequencing (ddRAD-seq). Although a species often shared chlorotypes with other species, our phylogenetical analyses based on variation in the three single-copy nuclear genes and the ddRAD-seq data showed that most populations belonging to the same species were monophyletic at the species level, suggesting that chloroplast capture may have frequently occurred between Isodon species. Some populations of an intraspecific taxon were embedded paraphyletically within the species, regardless of the large amount of phylogenetic information in nuclear DNA; this incongruity may have resulted from incomplete lineage sorting.