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.     

Phylogenetic and biogeographic diversification of Rhus (Anacardiaceae) in the Northern Hemisphere

Sequences of internal transcribed spacers (ITS) of nuclear ribosomal DNA, the chloroplast ndhF gene, and chloroplast trnL-F regions (trnL intron, and trnL [UAA] 3' exon-trnF [GAA] intergenic spacer) were used for phylogenetic analyses of Rhus, a genus disjunctly distributed in Asia, Europe, Hawaii, North America, and Northern Central America. Both ITS and cpDNA data sets support the monophyly of Rhus. The monophyly of subgenus Rhus was suggested by the combined cpDNA and ITS data, and largely supported in the cpDNA data except that Rhus microphylla of subgenus Lobadium was nested within it. The monophyly of subgenus Lobadium was strongly supported in the ITS data, whereas the cpDNA data revealed two main clades within the subgenus, which formed a trichotomy with the clade of subgenus Rhus plus R. microphylla. The ITS and cpDNA trees differ in the positions of Rhus michauxii, R. microphylla, and Rhus rubifolia, and hybridization may have caused this discordance. Fossil evidence indicates that Rhus dates back to the early Eocene. The penalized likelihood method was used to estimate divergence times, with fossils of Rhus subgenus Lobadium, Pistacia and Toxicodendron used for age constraints. Rhus diverged from its closest relative at 49.1+/-2.1 million years ago (Ma), the split of subgenus Lobadium and subgenus Rhus was at 38.1+/-3.0 Ma. Rhus most likely migrated from North America into Asia via the Bering Land Bridge during the Late Eocene (33.8+/-3.1 Ma). Rhus coriaria from southern Europe and western Asia diverged from its relatives in eastern Asia at 24.4+/-3.2 Ma. The Hawaiian Rhus sandwicensis diverged from the Asian Rhus chinensis at 13.5+/-3.0 Ma. Subgenus Lobadium was inferred to be of North American origin. Taxa of subgenus Lobadium then migrated southward to Central America. Furthermore, we herein make the following three nomenclatural combinations: (1) Searsia leptodictya (Diels) T. S. Yi, A. J. Miller and J. Wen, comb. nov., (2) Searsia pyroides (A. Rich.) T. S. Yi, A. J. Miller and J. Wen, comb. nov., and (3) Searsia undulata (Jacq.) T. S. Yi, A. J. Miller and J. Wen, because our analyses support the segregation of Searsia from Rhus.     

A phylogenetic analysis ofPanax (Araliaceae): Integrating cpDNA restriction site and nuclear rDNA ITS sequence data

A phylogenetic analysis ofPanax was performed using restriction site variations of eight PCR-amplified chloroplast regions. Twenty populations were examined, representing 13 of the 14 species ofPanax. Aralia cordata was used as the outgroup. The 11 restriction endonucleases produced a total of 105 restriction sites and length variations from the large single-copy region of cpDNA. Forty restriction variations are polymorphic. The cpDNA tree is largely congruent with the nuclear ribosomal ITS phylogeny. Similar to the ITS tree, the cpDNA dataset suggests the following relationships: (1)P. trifolius from eastern North America is sister to the clade consisting of all otherPanax species; (2)P. ginseng andP. japonicus from eastern Asia form a clade withP. quinquefolius from eastern North America; (3) the HimalayanP. pseudoginseng is most closely related toP. stipuleanatus of southwestern China; and (4) the medicinally importantP. notoginseng forms a clade with the closely relatedP. bipinnatifidus, P. ginseng, P. japonicus, P. major, P. quinquefolius, P. sinensis, P. wangianus, andP. zingiberensis. Two biogeographic disjunctions are detectable withinPanax. One is the connection of the eastern North AmericanP. trifolius with the rest ofPanax species. The other is the more recent disjunction between the North AmericanP. quinquefolius and the eastern AsianP. ginseng andP. japonicus. The active orogenies caused by the collision of the Indian Plate with Asia may have facilitated the diversification ofPanax taxa in Asia in the late Tertiary.     

Phylogeny and biogeography of the flowering plant genus Styrax (Styracaceae) based on chloroplast DNA restriction sites and DNA sequences of the internal transcribed spacer region

Phylogenetic relationships within the flowering plant genus Styrax were investigated with DNA sequence data from the internal transcribed spacer (ITS) region of nuclear ribosomal DNA (nrDNA) and with chloroplast DNA restriction site data from the genes trnK, rpoC1, and rpoC2. The data sets from each genome were analyzed separately and in combination with parsimony methods. The results strongly support the monophyly of each of the four series of the genus but provide little phylogenetic resolution among them. Reticulate evolution may at least partly explain discordance between the molecular phylogenetic estimates and a prior morphological estimate within series Cyrta. The historical biogeography of the genus was inferred with unweighted parsimony character optimization of trees recovered from a combined ITS and morphological data set, after a series of combinability tests for data set congruence was conducted. The results are consistent with the fossil record in supporting a Eurasian origin of Styrax. The nested phylogenetic position of the South American members of the genus within those from southern North America and Eurasia suggests that the boreotropics hypothesis best explains the amphi-Pacific tropical disjunct distribution occurring within section Valvatae. The pattern of relationship recovered among the species of section Styrax ((western North America + western Eurasia) (eastern North America + eastern Eurasia)) is rare among north-temperate Tertiary forest relicts. The monophyly of the group of species from western North America and western Eurasia provides qualified support for the Madrean-Tethyan hypothesis, which posits a Tertiary floristic connection among the semiarid regions in which these taxa occur. A single vicariance event between eastern Asia and eastern North America accounts for the pattern of relationship among intercontinental disjuncts in series Cyrta.     

Phylogeny and Biogeography of Thuja L. (Cupressaceae), an Eastern Asian and North American Disjunct Genus

In order to develop better insights into biogeographic patterns of eastern Asian and North American disjunct plant genera, sequences of nuclear ribosomal DNA internal transcribed spacer (nr DNA ITS) region were used to estimate interspecific relationships of Thuja L. (Cupressaceae) and infer its biogeography based on the phylogeny. According to the phylogenetic analysis, two clades were recognized. The first clade included Thuja plicata D. Don (western North America) and T. koraiensis Nakai (northeastern Asia), and the second one contained T. occidentalis (Gord.) Carr. (Japan). The ancestral area of Thuja was inferred to be eastern Asia, and two dispersal events were responsible for the modern distribution of Thuja in North America. Both the North Atlantic land bridge and Bering land bridge were possible routes for the migration of ancestral populations to North America.     

Phylogenetic utility of nuclear nitrate reductase: A multi-locus comparison of nuclear and chloroplast sequence data for inference of relationships among American Lycieae (Solanaceae)

There has been considerable interest and research into the benefits of multiple low-copy nuclear regions for phylogenetic studies at low taxonomic levels. In this study, the phylogenetic utility of DNA sequence data from two low-copy nuclear genes, nitrate reductase (NIA) and granule-bound starch synthase I (GBSSI), was compared with data from nuclear ribosomal ITS and the cpDNA spacers trnT–trnF and trnD–trnT across 33 closely related taxa in tribe Lycieae (Solanaceae). The NIA data (introns 1–2) had the most parsimony-informative characters, with over twice the number provided by GBSSI, but NIA also had greater homoplasy. Although gene trees were generally concordant across the four regions, there were some notable areas of incongruence, suggesting both incomplete lineage sorting as well as possible reticulate origins.     

Phylogenetics and biogeography of eastern Asian-North American disjunct genus Pachysandra （Buxaceae） inferred from nucleotide sequences

Pachysandra is an eastern Asian-North American disjtunct genus with three species, two in eastern Asia （Pachysandra axillaris and Pachysandra terminalis） and one in eastern North America （Pachysandra procurnbens）. Although morphological and cytological studies suggest a close affinity of Pprocumbens with P axillaris, molecular data from nuclear and chloroplast DNA regions have provided conflicting signals. In this study, we tested previous phylogenetic hypotheses using sequences of nuclear ribosomal DNA internal transcribed spacers and chloroplast ndhF gene from multiple individuals of each of the three species. We also estimated the time of divergence between eastem Asia and eastern North America. Our results support the morphological and cytological conclusion that P procumbens is more closely related to P axillaris than to P terminalis. The estimated time of divergence of P axillaris and P procumbens was 14.6±5.5 mya, consistent with estimates from many other eastern Asian-North American disjunct genera. The migration of Pachysandra populations from eastern Asia to North America might have occurred by way of the North Atlantic land bridge.     

Polyphyly of the Padus group of Prunus (Rosaceae) and the evolution of biogeographic disjunctions between eastern Asia and eastern North America

Prunus subgenus Padus is a group with a wide distribution in temperate eastern Asia and eastern North America with one species extending to Europe and one to Central America. Phylogenetic relationships of subgenus Padus were reconstructed using sequences of nuclear ribosomal ITS, and plastid ndhF gene, and rps16 intron and rpl16 intron. Prunus subgenus Padus is shown to be polyphyletic. Taxa of subgenus Padus and subgenus Laurocerasus are highly intermixed in both the ITS and the plastid trees. The results support two disjunctions between eastern North America and Eurasia within the Padus group. One disjunction is between Prunus virginiana of eastern North America and P. padus of Eurasia, estimated to have diverged at 2.99 (95 % HPD 0.59–6.15)–4.1 (95 % HPD 0.63–8.59) mya. The other disjunction is between P. serotina and its Asian relatives. The second disjunction may have occurred earlier than the former one, but the age estimate is difficult due to the unresolved phylogenetic position of the P. serotina complex.     

Phylogenetic Use of Noncoding Regions in the GenusGentianaL.: ChloroplasttrnL (UAA) Intron versus Nuclear Ribosomal Internal Transcribed Spacer Sequences

Sequence divergence was estimated within noncoding sequences of both chloroplast DNA (cpDNA)trnL (UAA) intron and nuclear ribosomal DNA (nrDNA) internal transcribed spacer sequences (ITS1 and ITS2) for 10 species of the genusGentianaL. (Gentianaceae). Comparisons of evolutionary rates among these sequences (cpDNA versus nrDNA, ITS1 versus ITS2) were performed. It appears that sequence divergence is on average two to three times higher in ITSs than in thetrnL intron sequences and higher in ITS1 than in ITS2. Both the cpDNA intron and ITSs of nrDNA give concordant phylogenetic trees. However, the ITS-based phylogeny displays higher bootstrap values. At the intrageneric level, at least inGentiana,ITSs (especially ITS2) sequences seem to be more appropriate in the assessment of plant phylogenies. Nevertheless, the cpDNAtrnL intron seems to be preferable at the intergeneric level.     

Divergent and reticulate species relationships in Leucaena (Fabaceae) inferred from multiple data sources: insights into polyploid origins and nrDNA polymorphism

Previous analyses of species relationships and polyploid origins in the mimosoid legume genus Leucaena have used chloroplast DNA (cpDNA) restriction site data and morphology. Here we present an analysis of a new DNA sequence data set for the nuclear ribosomal DNA (nrDNA) 5.8S subunit and flanking ITS 1 and ITS 2 spacers, a simultaneous analysis of the morphology, ITS and cpDNA data sets for the diploid species, and a detailed comparison of the cpDNA and ITS gene trees, which include multiple accessions of all five tetraploid species. Significant new insights into species relationships and polyploid origins, including that of the economically important tropical forage tree L. leucocephala, are discussed. Heterogeneous ITS copy types, including 26 putative pseudogene sequences, were found within individuals of four of the five tetraploid and one diploid species. Potential pseudogenes were identified using two pairwise comparison approaches as well as a tree-based method that compares observed and expected proportions of total ITS variation contributed by the 5.8S subunit optimized onto branches of one of the ITS gene trees. Inclusion of putative pseudogene sequences in the analysis provided evidence that some pseudogenes in allopolyploid L. leucocephala are not the result of post-allopolyploidization gene silencing, but were inherited from its putative diploid maternal progenitor L. pulverulenta.     

Phylogenetic relationships and biogeographic patterns in North American members of Carex section Acrocystis (Cyperaceae) using nrDNA ITS and ETS sequence data

Carex section Acrocystis currently includes 27 taxa in North America. Recent phylogenetic studies have suggested that the North American and some but not all of the Eurasian species form a clade. Relationships and biogeographic patterns among species in this core-Acrocystis group are explored here using nuclear ribosomal (nrDNA) internal transcribed spacer region (ITS) and nrDNA external transcribed spacer region (ETS) sequence data. While maximum parsimony analysis of the ITS and ETS data provides only a moderately resolved branching structure for species relationships within the core-Acrocystis clade, maximum likelihood analysis provides a more resolved hypothesis of relationships in the section. The core-Acrocystis clade consists of a grade of Eurasian and primarily western North American species, with a well-supported clade of only eastern North American species nested within this grade. ITS and ETS types do not coalesce within many species or species complexes. Possible explanations for the non-coalescent nature of ITS and ETS copies in Acrocystis are explored, including lineage sorting, hybridization, and cryptic species.     

Incongruence among mitochondrial,chloroplast and nuclear gene trees in <Emphasis Type="Italic">Pinus</Emphasis> subgenus <Emphasis Type="Italic">Strobus</Emphasis> (Pinaceae)

Introgression has been considered to be one of main factors leading to phylogenetic incongruence among different datasets at lower taxonomic levels. In the plants of Pinaceae, the mtDNA, cpDNA, and nuclear DNA (nrDNA) may have different evolutionary histories through introgression because they are inherited maternally, paternally and biparentally, respectively. We compared mtDNA, cpDNA, and two low-copy nrDNA phylogenetic trees in the genus Pinus subgenus Strobus, in order to detect unknown past introgression events in this group. nrDNA trees were mostly congruent with the cpDNA tree, and supported the recent sectional and subsectional classification system. In contrast, mtDNA trees split the members of sect. Quinquefoliae into two groups that were not observed in the other gene trees. The factors constituting incongruence may be divided into the following two categories: the different splits within subsect. Strobus, and the non-monophyly of subsect. Gerardianae. The former was hypothesized to have been caused by the past introgression of cpDNA, mtDNA or both between Eurasian and North American species through Beringia. The latter was likely caused by the chimeric structure of the mtDNA sequence of P. bungeana, which might have originated through past hybridization, or through a horizontal transfer event and subsequent recombination. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

A species-level phylogenetic study of the Verbena complex (Verbenaceae) indicates two independent intergeneric chloroplast transfers

Two major impediments to infer plant phylogenies at inter- or intra- species level include the lack of appropriate molecular markers and the gene tree/species tree discordance. Both of these problems require more extensive investigations. One of the foci of this study is examining the phylogenetic utility of a combined chloroplast DNA dataset (>5.0kb) of seven non-coding regions, in comparison with that of a large fragment (ca. 3.0kb) of a low-copy nuclear gene (waxy), in a recent, rapidly diversifying group, the Verbena complex. The complex includes three very closely related genera, Verbena (base chromosome number x=7), Glandularia (x=5), and Junellia (x=10), comprising some 150 species distributed predominantly in South and North America. Our results confirm the inadequacy of non-coding cpDNA in resolving relationships among closely related species due to lack of variation, and the great potential of low-copy nuclear gene as source of variation. However, this study suggests that when both cpDNA and nuclear DNA are employed in low-level phylogenetic studies, cpDNA might be very useful to infer organelle evolutionary history (e.g., chloroplast transfer) and more comprehensively understand the evolutionary history of organisms. The phylogenetic framework of the Verbena complex resulted from this study suggests that Junellia is paraphyletic and most ancestral among the three genera; both Glandularia and Verbena are monophyletic and have been derived from within Junellia. Implications of this phylogenetic framework to understand chromosome number evolution and biogeography are discussed. Most interestingly, the comparison of the cpDNA and nuclear DNA phylogenies indicates two independent intergeneric chloroplast transfers, both from Verbena to Glandularia. One is from a diploid North American Verbena species to a polyploid North American Glandularia species. The other is more ancient, from the South American Verbena group to the common ancestor of a major Glandularia lineage, which has radiated subsequently in both South and North America. The commonly assumed introgressive hybridization may not explain the chloroplast transfers reported here. The underlying mechanism remains uncertain.     

The role of disjunction and postglacial population expansion on phylogeographical history and genetic diversity of the circumboreal plant Chamaedaphne calyculata

In the last decade a number of studies has illustrated quite different phylogeographical patterns amongst plants with a northern present‐day geographical distribution, spanning the entire circumboreal region and/or circumarctic region and southern mountains. These works, employing several marker systems, have brought to light the complex evolutionary histories of this group. Here I focus on one circumboreal plant species, Chamaedaphne calyculata (leatherleaf), to unravel its phylogeographical history and patterns of genetic diversity across its geographical range. A survey of 29 populations with combined analyses of chloroplast DNA (cpDNA), internal transcribed spacer (ITS) and AFLP markers revealed structuring into two groups: Eurasian/north‐western North American, and north‐eastern North American. The present geographical distribution of C. calyculata has resulted from colonization from two putative refugial areas: east Beringia and south‐eastern North America. The variation of chloroplast DNA (cpDNA) and ITS sequences strongly indicated that the evolutionary histories of the Eurasian/north‐western North American and the north‐eastern North American populations were independent of each other because of a geographical disjunction in the distribution area and ice‐sheet history between north‐eastern and north‐western North America. Mismatch analysis using ITS confirmed that the present‐day population structure is the result of rapid expansion, probably since the last glacial maximum. The AFLP data revealed low genetic diversity of C. calyculata (P = 19.5%, H = 0.085) over the whole geographical range, and there was no evidence of loss of genetic diversity within populations in the continuous range, either at the margins or in formerly glaciated and nonglaciated regions. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 105 , 761–775.     

Molecular phylogeny and phylogeography ofLupinus (Leguminosae) inferred from nucleotide sequences of therbcL gene and ITS 1 + 2 regions of rDNA

Total DNA was extracted from 55 species of theLeguminosae (including 29 species ofLupinus). The chloroplast generbcL and the ITS 1 + 2 regions of nuclear RNA genes were amplified by polymerase chain reaction (PCR) and sequenced directly. The sequences obtained were evaluated with character state (Maximum Parsimony) and distance methods (Neighbour Joining). Phylogenetic trees obtained with both data sets and methods are mostly congruent.Genisteae andCrotalarieae are sister groups and share ancestry with theThermopsideae/Podalyrieae. The genusLupinus, which forms a monophyletic clade within theGenisteae, shows a distinct Old-New World disjunction and appears to be divided into several more or less distinct groups: (1) The species from the eastern part of South America. (2) The homogeneous rough-seeded group (Scabrispermae) of the Old World species which is well distinguished from the smooth-seeded group (Malacospermae). (3) Within the rather heterogeneous smooth-seeded lupins a smaller subgroup withL. angustifolius, L. hispanicus andL. luteus is recognized. (4) Also separated are North American lupins and South American species with a western distribution. Genetic distances imply that the genusLupinus evolved during the last 12–14 million years, ruling out the hypothesis that the present Old-New World disjunction can be interpreted as a result of the continental drift. The genetic data suggest an origin in the Old World and an independant colonisation of the Eastern parts of South America as opposed to North America and the Western parts of South America.     

Out of the Qinghai-Tibet Plateau: evidence for the origin and dispersal of Eurasian temperate plants from a phylogeographic study of Hippophaë rhamnoides (Elaeagnaceae)

Numerous temperate plants now distributed across Eurasia are hypothesized to have originated and migrated from the Qinghai-Tibet Plateau (QTP) and adjacent regions. However, this hypothesis has never been tested through a phylogeographic analysis of a widely distributed species. Here, we use Hippopha? rhamnoides as a model to test this hypothesis. We collected 635 individuals from 63 populations of the nine subspecies of H. rhamnoides. We sequenced two maternally inherited chloroplast (cp) DNA fragments and also the bi-paternally inherited nuclear ribosomal ITS. We recovered five major clades in phylogenetic trees constructed from cpDNA and internal transcribed spacer (ITS) sequence variation. Most sampled individuals of six subspecies that are distributed in northern China, central Asia and Asia Minor/Europe, respectively, comprised monophyletic clades (or subclades) nested within those found in the QTP. Two subspecies in the QTP were paraphyletic, while the placement of another subspecies from the Mongolian Plateau differed between the ITS and cpDNA phylogenetic trees. Our phylogeographic analyses supported an 'out-of-QTP' hypothesis for H. rhamnoides followed by allopatric divergence, hybridization and introgression. These findings highlight the complexity of intraspecific evolutions and the importance of the QTP as a center of origin for many temperate plants.     

Genetic evidence for host specificity in the hemi-parasitic Viscum album L. (Viscaceae)

Nuclear ribosomal DNA (nrDNA) ITS sequences and partial sequences of three non-coding chloroplast DNA (cpDNA) introns and spacers were used to assess genetic variation within and among three presumed host races of the hemi-parasite Viscum album L. Currently, identification of host races occurs via the host trees, and morphological differences are minute at best. cpDNA and nrDNA ITS sequences revealed little sequence variation, but the variation found consistently supported the distinction of three host races. cpDNA and ITS sequences were not incongruent, as assessed by the incongruence length difference test. A combined analysis supported the sister group relationship between mistletoes from deciduous trees and fir.     

Molecular phylogenetic systematics and biogeography of tribe Neillieae (Rosaceae) using DNA sequences of cpDNA, rDNA, and LEAFY

A phylogeny of the tribe Neillieae (Rosaceae), which comprises Neillia, Stephanandra, and Physocarpus, was reconstructed based on nucleotide sequences of several regions of cpDNA, the ITS and ETS regions of rDNA, and the second intron of LEAFY, to elucidate relationships among genera and species in Neillieae and to assess the historical biogeography of the tribe. Phylogenetic analyses indicated that Physocarpus and Neillia-Stephanandra were strongly supported as monophyletic and suggested that Stephanandra may have originated by hybridization between two lineages of Neillia. Dispersal-vicariance analyses suggested that the most recent common ancestor of Neillieae may have occupied eastern Asia and western North America and that Physocarpus and Neillia-Stephanandra may have been split by an intercontinental vicariance event in the early Miocene. The biogeographic analyses also suggested that species of Neillia and Stephanandra diversified in eastern Asia, whereas in Physocarpus one dispersal event from western North America to eastern Asia occurred. Two divergent types of LEAFY sequences were found in the eastern North American species, P. opulifolius, but only one type was present in each plant. The two types of sequences may represent homeologous genes that originated by hybridization between P. capitatus and P. monogynus, both western North American species.     

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.     

Recent vicariance and the origin of the rare, edaphically specialized Sandhills lily, Lilium pyrophilum (Liliaceae): evidence from phylogenetic and coalescent analyses

Establishing the phylogenetic and demographic history of rare plants improves our understanding of mechanisms that have led to their origin and can lead to valuable insights that inform conservation decisions. The Atlantic coastal plain of eastern North America harbours many rare and endemic species, yet their evolution is poorly understood. We investigate the rare Sandhills lily (Lilium pyrophilum), which is endemic to seepage slopes in a restricted area of the Atlantic coastal plain of eastern North America. Using phylogenetic evidence from chloroplast, nuclear internal transcribed spacer and two low-copy nuclear genes, we establish a close relationship between L. pyrophilum and the widespread Turk's cap lily, L. superbum. Isolation-with-migration and coalescent simulation analyses suggest that (i) the divergence between these two species falls in the late Pleistocene or Holocene and almost certainly post-dates the establishment of the edaphic conditions to which L. pyrophilum is presently restricted, (ii) vicariance is responsible for the present range disjunction between the two species, and that subsequent gene flow has been asymmetrical and (iii) L. pyrophilum harbours substantial genetic diversity in spite of its present rarity. This system provides an example of the role of edaphic specialization and climate change in promoting diversification in the Atlantic coastal plain.