Molecular evidence for the age, origin, and evolutionary history of the American desert plant genus Tiquilia (Boraginaceae)

Although the deserts of North America are of very recent origin, their characteristic arid-adapted endemic plant lineages have been suggested to be much older. Earlier researchers have hypothesized that the ancestors of many of these modern desert lineages first adapted to aridity in highly localized arid or semi-arid sites as early as the late Cretaceous or early Tertiary, and that these lineages subsequently spread and diversified as global climate became increasingly arid during the Cenozoic. No study has explicitly examined these hypotheses for any North American arid-adapted plant group. The current paper tests these hypotheses using the genus Tiquilia (Boraginaceae), a diverse North American desert plant group. A strongly supported phylogeny of the genus is estimated using combined sequence data from three chloroplast markers (matK, ndhF, and rps16) and two nuclear markers (ITS and waxy). Ages of divergence events within the genus are estimated using penalized likelihood and a molecular clock approach on the ndhF tree for Tiquilia and representative outgroups, including most of the major lineages of Boraginales. The dating analysis suggests that the stem lineage of Tiquilia split from its nearest extant relative in the Paleocene or Eocene ( approximately 59-48 Ma). This was followed by a relatively long period before the first divergence in the crown group near the Eocene/Oligocene boundary ( approximately 33-29 Ma), shortly after the greatest Cenozoic episode of rapid aridification. Divergence of seven major lineages of Tiquilia is dated to the early-to-mid Miocene ( approximately 23-13 Ma). Several major lineages show a marked increase in diversification concomitant with the onset of more widespread semi-arid and then arid conditions beginning in the late Miocene ( approximately 7 Ma). This sequence of divergence events in Tiquilia agrees well with earlier researchers' ideas concerning North American desert flora assembly.     

Monophyly of Kelloggia Torrey ex Benth. (Rubiaceae) and evolution of its intercontinental disjunction between western North America and eastern Asia

Kelloggia Torrey ex Bentham (Rubiaceae) consists of two species disjunctly distributed in western North America (K. galioides Torrey) and the western part of eastern Asia (K. chinensis Franch.). The two species exhibit a high level of morphological divergence. To test its monophyly and to infer its biogeographic history, we estimated the phylogeny of Kelloggia and its relatives from sequences of three chloroplast DNA regions (rbcL gene, atpB-rbcL spacer, and rps16 intron). The monophyly of Kelloggia was strongly supported, and it forms a sister relationship with the tribe Rubieae. The divergence time between the two disjunct species of Kelloggia was estimated to be 5.42 ± 2.32 million years ago (mya) using the penalized likelihood method based on rbcL sequence data with fossil calibration. Our result does not support the Madrean-Tethyan hypothesis, which assumes an earlier divergence time of 20-25 mya. Ancestral area analysis, as well as dispersal-vicariance (DIVA) analysis, suggests the Asian origin of Kelloggia and the importance of Eurasia in the diversification of its close relatives in the Rubieae-Theligoneae-Paederieae group. The intercontinental disjunction in Kelloggia is suggested to have evolved via long-distance dispersal from Asia into western North America.     

Phylogeny and Phytogeography of Eupatorium (Eupatorieae, Asteraceae): Insights from Sequence Data of the nrDNA ITS Regions and cpDNA RFLP

Eupatorium were examined by sequencing the internal transcribed spacers (ITS) of nuclear ribosomal DNA and restriction site analysis of chloroplast DNA. Molecular data provided strong evidence that (1) this genus originated in North America, (2) the genus diverged into three morphological species groups, Eutrochium, Traganthes and Uncasia in North America, and (3) one of the North American Uncasia lineages migrated into temperate Europe and eastern Asia over the Bering land bridge. The estimated divergence times support a late Miocene to early Pliocene migration from North America to Eurasia via the Bering land bridge. A European species was sister to all of the eastern Asian species examined. The disjunct distribution pattern of the genus Eupatorium is incongruent with the classical Arcto-Tertiary geoflora concept. Received 13 September 1999/ Accepted in revised form 4 January 2000     

Phylogeography of the Northern Alligator Lizard (Squamata,Anguidae): Hidden diversity in a western endemic

Western North America includes the California Floristic Province and the Pacific Northwest, biologically diverse regions highlighted by a complex topography, geology, climate and history. A number of animals span these regions and show distinctive patterns of dispersal, vicariance and lineage diversification. Examining phylogeographic patterns in the fauna of this area aids in our understanding of the forces that have contributed to the generation and maintenance of regional biodiversity. Here, we investigate the biogeography and population structure of the Northern Alligator Lizard (Elgaria coerulea), a wide‐ranging anguid endemic to western North America. We sequenced two mtDNA fragments (ND2 and ND4) for 181 individuals across the range of the species and analysed these data with phylogenetic approaches to infer population and biogeographic history, and date major divergences within the taxon. We further used Bayesian clustering methods to assess major patterns of population structure and performed ecological niche modelling (ENM) to aid in our interpretation of geographic structure and diversification of E. coerulea lineages. Our phylogeographic examination of E. coerulea uncovered surprising diversity and structure, recovering 10 major lineages, each with substantial geographic substructure. While some divergences within the species are relatively old (Pliocene, 5.3–2.6 mya), most intraspecific variation appears to be of more recent origin (Pleistocene, 2.6 mya‐11,700 ya). Current diversity appears to have arisen in the Sierra Nevada Mountains and spread west and north since the Pliocene. Finally, our ENMs suggest that much of the Coast Ranges in California provided ideal habitat during the Last Glacial Maxima (LGM) that has since contracted dramatically and shifted northwards, whereas significant portions of the Sierra Nevada were unsuitable during the LGM and have since become more suitable. Interestingly, E. coerulea shares a number of genetic boundaries with other sympatric taxa, suggesting common historical events and geomorphological features have shaped the biota of this region.     

Phylogeny and biogeography of Alyssum (Brassicaceae) based on nuclear ribosomal ITS DNA sequences

The genus Alyssum consists of about 195 species native to Europe, Asia and northern Africa. All species were assigned to six sections. Previous molecular phylogeny studies indicate that Alyssum is polyphyletic. However, the divergence time and dispersal of the genus are not well studied. In this study, the phylogenetic relationships within the genus Alyssum were studied with nrDNA ITS sequences obtained from five sections. The divergence time was estimated by fossil calibration and the biogeography was examined by spread analysis. The phylogeny indicated two main lineages: lineage 1 includes the section of Alyssum, Gamosepalum and Psilonema; lineage 2 includes the section of Odontarrhena, Meniocus and Clypeola. The phylogenetic relationship was not congruent with the previous sectional classifications. The age of Alyssum was dated to the upper Miocene. Molecular data suggested the diversification of Alyssum in Mediterranean areas and wide-ranging distribution such as North Africa, eastward into Central Asia and immigration into North America. Climatic aridification and arid/semiarid areas established in the Pliocene/Pleistocene could have provided favourable conditions for the migration and diversification of Alyssum.     

Phylogenetic assessment and biogeographic analyses of tribe <Emphasis Type="Italic">Peracarpeae</Emphasis> (Campanulaceae)

Peracarpeae is a small tribe consisting of three genera: Homocodon, Heterocodon and Peracarpa, with a disjunct distribution between eastern Asia and western North America. Homocodon is endemic to southwestern China and was previously placed in the western North American genus Heterocodon. Our phylogenetic analysis using four plastid markers (matK, atpB, rbcL and trnL-F) suggests the polyphyly of Peracarpeae. Homocodon is sister to a clade consisting of the eastern Asian Adenophora, Hanabusaya and species of Campanula from the Mediterranean region and North America, rather than forming a clade with Heterocodon. Homocodon and its Eurasia relatives are estimated to have diverged in the early Miocene (16.84 mya, 95% HPD 13.35–21.45 mya). The eastern Asian Peracarpa constitutes a clade with the North American Heterocodon, Githopsis and three species of Campanula, supporting a disjunction between eastern Asia and North America in Campanulaceae. The Asian-North American disjunct lineages diverged in the early Miocene (16.17 mya, 95% HPD 13.12–20.9 mya). The biogeographic analyses suggest that Homocodon may be a relict of an early radiation in eastern Asia, and that Peracarpa and its closest North American relatives most likely originated from a Eurasian ancestor.     

Temporal patterns of diversification and microendemism in Eastern Highland endemic barcheek darters (Percidae: Etheostomatinae)

Eastern North America is the location of the world's most species-rich temperate freshwater fish fauna. Hypotheses regarding the geographic and temporal scale of teleost diversification in this region have not been broadly investigated using absolute divergence time estimates among the constituent lineages. This study used time-calibrated molecular phylogenies estimated from mitochondrial and nuclear genes to investigate the temporal and geographic signatures of diversification within barcheek darters, a clade of allopatrically distributed species endemic to the Eastern Highlands. Results from divergence time estimates using an uncorrelated lognormal model suggest that the barcheek darters are an ancient group with a crown node estimate of 16.3 mya, 95% highest posterior density (HPD): [12.4, 20.5], and the clade is characterized by substantial intraspecific divergence times within several species. In particular, the Caney Fork endemic Etheostoma basilare comprises five strongly supported and deeply divergent clades with a most recent common ancestor estimated at 8.0 mya, 95% HPD: [5.6, 10.7]. These results are concordant with the hypothesis that geologically stable areas of eastern North America have facilitated both the generation and preservation of lineages across a substantial breadth of evolutionary time, and that allopatric speciation in darters has occurred at much smaller spatial scales than previously realized.     

Phylogeny and biogeography of Rhododendron subsection Pontica, a group with a tertiary relict distribution

Rhododendron subgenus Hymenanthes subsection Pontica is exceptional among Tertiary relict groups in having a high proportion of species (4 of 11) native to SW Eurasia. A phylogeny based on cpDNA matK and trnL-F indicated that multiple Pontica lineages colonised each of SW Eurasia, SE North America, and NE Asia, with little or no speciation within regions thereafter. Therefore, multiple (3-4) Pontica lineages survived the Quaternary in SW Eurasia, in contrast to other Tertiary relict genera. Pontica comprises two major clades, one of which is wholly Eurasian, and paraphyletic with respect to at least some of the remaining 200 species of subgenus Hymenanthes, which are all distributed in SE Asia. The other clade has species from W and SE North America, SW Eurasia, and NE Asia. According to synonymous matK substitution data, the two clades diverged 9-6 million years ago (mya), whereas most divergence within them happened 5-3 mya. Although the phylogeny indicates probable trans-Atlantic migration for one of two America-Eurasia disjunctions in Pontica, the timing supports migration via Beringia for both.     

A biogeographical analysis of Muhlenbergia (Poaceae: Chloridoideae: Cynodonteae: Muhlenbergiinae)

A phylogeny based on the analysis of six DNA sequence markers (ITS, ndhA intron, rpl32-trnL, rps3, rps16 intron, and rps16-trnK) is used to infer ancestral areas and divergence times, and reconstruct the biogeographical history and evolution of 150 of the 183 (82%) species of Muhlenbergia. Our results suggest that the genus originated 9.3 mya in the Sierra Madre (Occidental and Oriental) in Mexico, splitting into six lineages: M. ramulosa diverging 8.2 mya, M. subg. Muhlenbergia at 5.9 mya, M. subg. Pseudosporobolus at 5.9 mya, M. subg. Clomena at 5.4 mya, M. subg. Bealia at 4.3 mya, and M. subg. Trichochloa at 1 mya, each of these with a high probability of Sierra Madrean origin. Our results further suggest that founder-event speciation from Sierra Madre to South America occurred independently multiple times in all five subgenera during the Pleistocene and late Pliocene. One long-distance dispersal event most likely originating from Central or Eastern North America to East and Central Asia occurred 1.6–1 mya in M. subg. Muhlenbergia. In our cladogram, members of M. subg. Trichochloa show little genetic resolution, suggesting very low levels of divergence among the species, and this may be a consequence of rapid radiation.     

The origins and evolution of the genus Myosotis L. (Boraginaceae)

Although morphologically well defined, the phylogeny and taxonomy of Myosotis has been uncertain. In particular it has been unclear whether the genus had a Northern Hemisphere or Australasian origin. However, separate analyses of the ITS and the 3' region of matK, as well as a combined analysis of ITS, 3' matK, the psbA-trnA spacer, and 3' ndhF regions indicate that several distinct lineages exist within Myosotis and strongly support a Northern Hemisphere origin for the genus. Further, the observed transoceanic distributions and levels of genetic divergence between lineages indicate that long distance dispersal has been important for establishing the current geographic range expansion of Myosotis. Our molecular data also suggest that the diversification of Australasian Myosotis has occurred since the late Tertiary and is largely due to radiation within and from New Zealand. This inference is consistent with the findings of recent phylogenetic studies on other New Zealand alpine genera. Our results highlight the important role played by late Tertiary and Quaternary climate change in explaining current floristic diversity. The genetic relationships reported here also suggest that the current infrageneric taxonomy of Myosotis does not fully reflect the evolution of the genus.     

Diversification in North American arid lands: Niche conservatism, divergence and expansion of habitat explain speciation in the genus Ephedra

A lineage of 12 arid land shrubby species in the gymnosperm genus Ephedra (Gnetales) from North America is used to evaluate the influence of climate on speciation. With a long evolutionary history, and a well documented fossil record this lineage is an ideal model for understanding the process of speciation under a niche conservatism scenario. Using seven DNA molecular markers, Bayesian inference is carried out to uncover sister species and to estimate time of divergence of the lineages. Ecological niche models are generated for four parapatric and sympatric sister species and two analyses of niche evolution are performed, one based on ecological niche models and another using raw data and multivariate analysis. As previous analyses suggest, the diversification of North America Ephedra species may be the result of a recent secondary radiation. Both parapatric and sympatric species diverged mostly in a scenario of climatic niche conservatism. However, we also found strong evidence for niche divergence for one of the sister species pairs (E. californica-E. trifurca). Moreover, the multivariate analysis found environmental differences for some variables between sister species. The estimated divergence time of three pairs of sister species distributed in southwestern North America (E. cutleri-E. aspera, E. californica-E. trifurca and E. torreyana-E. viridis) is inferred to have occurred in the Late Miocene to Pliocene and for the sister species pair E. antisyphilitica-E. coryi distributed in the southern United States and northeastern Mexico, it was inferred from the Pliocene to Pleistocene. The orogenetic and climatic changes documented for these regions related to expansion of arid lands, may have contributed to the diversification in North American Ephedra, rather than adaptations to new climatic conditions.     

Phylogenetic relationships and biogeography of Fuchsia (Onagraceae) based on noncoding nuclear and chloroplast DNA data

To examine relationships and test previous sectional delimitations within Fuchsia, this study used parsimony and maximum likelihood analyses with nuclear ITS and chloroplast trnL-F and rpl16 sequence data for 37 taxa representing all sections of Fuchsia and four outgroup taxa. Results support previous sectional delimitations, except for F. verrucosa, which is related to a Central American clade rather than to section Fuchsia and is described here as a new section Verrucosa. The basal relationships within Fuchsia are poorly resolved, suggesting an initial rapid diversification of the genus. Among the species sampled, there is strong support for a single South Pacific lineage, a southern South American/southern Brazilian lineage, a tropical Andean lineage, and one or two Central American and Mexican lineages. There is no clear support for an austral origin of the genus, as previously proposed, which is more consistent with Fuchsia's sister group relationship with the boreal Circaea. An ultrametric molecular clock analysis (all minimal dates) places the split between Fuchsia and Circaea at 41 million years ago (mya), with the diversification of the modern-day lineages of Fuchsia beginning at 31 mya. The South Pacific Fuchsia lineage branches off around 30 mya, consistent with fossil records from Australia and New Zealand. The large Andean section Fuchsia began to diversify around 22 mya, preceded by the divergence of the Caribbean F. triphylla at 25 mya. The Brazilian members of section Quelusia separated from the southern Andean F. magellanica around 13 mya, and the ancestor of the Tahitian F. cyrtandroides split off from the New Zealand species of section Skinnera approximately 8 mya.     

Phylogenetic relationships within the New World endemic Zeltnera (Gentianaceae-Chironiinae) inferred from molecular and karyological data

The New World endemic genus Zeltnera consists of 25 species mainly distributed in the western part of the United States and Mexico. Chromosome counts performed on 113 populations (24 species) reveal extensive congruence between chromosomal groups and the assemblages obtained from analyses of nuclear ribosomal DNA (ITS) and chloroplast DNA (trnL intron and trnL-F intergenic spacer) sequences. Karyological and molecular data sets support three main biogeographic groups for Zeltnera. A first and mainly unresolved cluster (n = 17 and n = 20) occurs in California, whereas two other clades are centered in the Texas region (n = 20 and n = 21) and in Mexico (n = 21 and n = 22). Under the assumption of a molecular clock, and using both dispersal and vicariance explanations for the current distribution of the respective species, the genus is thought to have a North American origin with considerable diversification in the early Pliocene (ca. 5 million years ago). Geological events, such as desert formation and mountain orogenies, have created insuperable barriers that today separate the three major and likely vicariant groups.     

Adaptive radiation of shrubby tarweeds (Deinandra) in the California Islands parallels diversification of the Hawaiian silversword alliance (Compositae-Madiinae)

Phylogenetic analyses of nuclear rDNA transcribed spacers and cytogenetic studies of interspecific hybrids reported here uphold Carlquist's hypothesis (1965, Island Biology) that shrubby tarweeds (Deinandra) of Guadalupe Island, Mexico, are products of in situ radiation in the California Islands, where evidence of plant diversification has been equivocal. Based on the rDNA findings, the Guadalupe Island endemics (D. frutescens, D. greeneana subsp. greeneana, and D. palmeri) constitute a clade that arose since the late Pliocene, well after the origin of Guadalupe Island and diversification of annual, mainland Californian lineages of Deinandra. High interfertility and normal meiosis in F(1) hybrids between the three endemics contrast with reduced interfertility (to complete intersterility) and meiotic irregularities in F(1) hybrids between other, mostly mainland species of Deinandra. Cloned rDNA sequences provided no convincing evidence of introgression among the Guadalupe Island deinandras; morphological, phenological, and/or habitat differences among those taxa indicate ecological barriers to gene flow and a probable role for ecological divergence in diversification. Biosystematic and molecular phylogenetic data for shrubby tarweeds of Guadalupe Island and another secondarily woody, oceanic-island tarweed lineage, the Hawaiian silversword alliance, reveal strikingly similar evolutionary histories. Both groups violate Baker's Rule by stemming from self-incompatible ancestors in western North America, and each has undergone within-island diversification without evolution of strong sterility barriers among lineages. Evolutionary parallels between these Hawaiian and California Island lineages of Madiinae, first suggested by Carlquist, may reflect characteristics of tarweeds that facilitate insular colonization and adaptive radiation.     

Phylogeographic diversification of antelope squirrels (Ammospermophilus) across North American deserts

We investigated the biogeographic history of antelope squirrels, genus Ammospermophilus, which are widely distributed across the deserts and other arid lands of western North America. We combined range‐wide sampling of all currently recognized species of Ammospermophilus with a multilocus data set to infer phylogenetic relationships. We then estimated divergence times within identified clades of Ammospermophilus using fossil‐calibrated and rate‐calibrated molecular clocks. Lastly, we explored generalized distributional changes of Ammospermophilus since the last glacial maximum using species distribution models, and assessed responses to Quaternary climate change by generating demographic parameter estimates for the three wide‐ranging clades of A. leucurus. From our phylogenetic estimates we inferred strong phylogeographic structure within Ammospermophilus and the presence of three well‐supported major clades. Initial patterns of historical divergence were coincident with dynamic alterations in the landscape of western North America, and the formation of regional deserts during the Late Miocene and Pliocene. Species distribution models and demographic parameter estimates revealed patterns of recent population expansion in response to glacial retreat. When combined with evidence from co‐distributed taxa, the historical biogeography of Ammospermophilus provides additional insight into the mechanisms that impacted diversification of arid‐adapted taxa across the arid lands of western North America. We propose species recognition of populations of the southern Baja California peninsula to best represent our current understanding of evolutionary relationships among genetic units of Ammospermophilus. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2013, 109 , 949–967.     

DNA barcoding of European Herbertus (Marchantiopsida, Herbertaceae) and the discovery and description of a new species

DNA barcoding of a group of European liverwort species from the genus Herbertus was undertaken using three plastid (matK, rbcL and trnH-psbA) and one nuclear (ITS) marker. The DNA barcode data were effective in discriminating among the sampled species of Herbertus and contributed towards the detection of a previously overlooked European Herbertus species, described here as H. norenus sp. nov. This species shows clear-cut differences in DNA sequence for multiple barcode regions and is also morphologically distinct. The DNA barcode data were also useful in clarifying taxonomic relationships of the European species with some species from Asia and North America. In terms of the discriminatory power of the different barcode markers, ITS was the most informative region, followed closely by matK. All species were distinguishable by ITS alone, rbcL + matK and various other multimarker combinations.     

Evolution of biogeographic disjunction between eastern Asia and eastern North America in Phryma (Phrymaceae)

This study examines molecular and morphological differentiation in Phryma L., which has only one species with a well-known classic intercontinental disjunct distribution between eastern Asia (EA) and eastern North America (ENA). Phylogenetic analysis of nuclear ribosomal ITS and chloroplast rps16 and trnL-F sequences revealed two highly distinct clades corresponding to EA and ENA. The divergence time between the intercontinental populations was estimated to be 3.68 ± 2.25 to 5.23 ± 1.37 million years ago (mya) based on combined chloroplast data using Bayesian and penalized likelihood methods. Phylogeographic and dispersal-vicariance (DIVA) analysis suggest a North American origin of Phryma and its migration into EA via the Bering land bridge. Multivariate analysis based on 23 quantitative morphological characters detected no geographic groups at the intercontinental level. The intercontinental populations of Phryma thus show distinct molecular divergence with little morphological differentiation. The discordance of the molecular and morphological patterns may be explained by morphological stasis due to ecological similarity in both continents. The divergence of Phryma from its close relatives in the Phrymaceae was estimated to be at least 32.32 ± 4.46 to 49.35 ± 3.18 mya.     

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.     

Rapid speciation and the evolution of hummingbird pollination in neotropical Costus subgenus Costus (Costaceae): evidence from nrDNA ITS and ETS sequences

We estimate phylogenetic relationships and the biogeographic and pollination history of Costus subgenus Costus (Costaceae) using sequence data from the internal and external transcribed spacer (ITS and ETS) regions of 18S-26S nuclear ribosomal DNA. The African members of the subgenus form a series of lineages basal to a monophyletic neotropical species radiation. The neotropical species have large, showy flowers visited by either euglossine bees or hummingbirds. The hummingbird pollination syndrome is supported as a derived character state from the bee pollination syndrome, and we estimate that it has evolved independently seven or more times in the neotropics. A molecular clock approach suggests that diversification of the neotropical clade has been recent and rapid and that it coincides with dramatic climatic and geologic changes, Andean orogeny, and the closing of the Panama isthmus that occurred in the Pliocene and Pleistocene epochs. We propose a scenario for the diversification of Costus, in which rapid floral adaptation in geographic isolation and range shifts in response to environmental changes contribute to reproductive isolation among close relatives. We suggest that these processes may be common in other recently diversified plant lineages centered in Central America or the Northern Andean phytogeographic region.     

Molecular relationships in Encephalartos (Zamiaceae, Cycadales) based on nucleotide sequences of nuclear ITS 1&2, rbcL, and genomic ISSR fingerprinting

The cycad genus Encephalartos is restricted to Africa and is threatened with extinction in most of its range. Total DNA was extracted from 51, i.e., 78 %, of the described species of Encephalartos. The accessions were sampled from the furthest western occurrence of the genus in Nigeria, via Sudan and Uganda, to southern South Africa. The sequences of nuclear ribosomal internal transcribed spacer regions 1 and 2 (ITS 1&2), the chloroplast encoded rbcL gene, and ISSR genomic fingerprinting were employed to resolve the molecular history and the relationships within the genus. Sequence alignment, as well as ISSR fingerprinting, data show low genetic variation among all analysed accessions, indicating diversification within the Pliocene/Pleistocene. ITS 1&2 data agree well with morphological and geographical characters and resolved three major genetic clusters with overlapping distribution ranges in eastern South Africa. This area, that contains the largest diversity of genotypes of Encephalartos, may have served as a Pliocene/Pleistocene refugium.