Available online at www.sciencedirect.com

South African Journal of Botany 75 (2009) 560–572 www.elsevier.com/locate/sajb

Nuclear and chloroplast DNA-based phylogenies of Tourn. ex Medik. (; ) reveal extensive incongruence and generic paraphyly, but support the recognition of infraspecific taxa in C. monilifera ⁎ N.P. Barker a, , S. Howis a, B. Nordenstam b, M. Källersjö b, P. Eldenäs b, C. Griffioen c, H.P. Linder c, 1

a Molecular Ecology and Systematics Group, Department of Botany, Rhodes University, Grahamstown 6140, South Africa b Swedish Museum of Natural History, Box 50007, S-104 05 Stockholm, Sweden c Botany Department, University of Cape Town, Private Bag, Rondebosch, 7700, South Africa Received 3 December 2008; received in revised form 28 May 2009; accepted 29 May 2009

Abstract

The small Chrysanthemoides comprises two species within which a number of infraspecific taxa have been recognized, some of which are invasive aliens in Australia and New Zealand. Here we investigate the relationships of the species and infraspecific taxa using both chloroplast and nuclear non-coding DNA sequence data. Results of the analyses of the plastid and nuclear data sets are incongruent, and neither Chry- santhemoides nor is resolved as monophyletic, although there is some support for the recognition of infraspecific taxa. Analyses of the separate and combined data sets resolve two clades within Chrysanthemoides (which include some species of Osteospermum), and these appear to have a geographic basis, one being restricted to the mainly winter rainfall region, the other the eastern bi-seasonal rainfall area. Our results suggest that there is evidence of past or ongoing hybridization within and possibly between these two lineages. © 2009 SAAB. Published by Elsevier B.V. All rights reserved.

Keywords: Biocontrol; Chrysanthemoides; Genetic diversity; Hybridization; Incongruence; Intraspecific variation; Osteospermum; Phylogeny

1. Introduction included in Osteospermum by Norlindh (1943, 1977). Norden- stam (1996) also questioned the monophyly of the reduced Os- The tribe Calenduleae comprises some 120 species and 12 teospermum, noting that despite these taxonomic changes, it was genera (Nordenstam, 2006, 2007; Nordenstam and Källersjö, still heterogeneous. Further splits from Osteospermum were 2009). Earlier assessments of the tribe recognized seven described as the new genera Norlindhia B. Nord. and (Norlindh, 1977) or ten genera (Nordenstam, 1994). Since 1994 B. Nord. (Nordenstam, 2006), and a species of Cass. there have been several taxonomic re-arrangements. Nordenstam was recognized as a distinct new genus Nephrotheca B. Nord. & (1994, 1996) considered the previously recognized genus Cas- Källersjö (Nordenstam et al., 2006). This move left Gibbaria as a talis Cass. and sect. Blaxium (Cass.) T. Norl. of Osteospermum monotypic genus until a second species was recently added by a L. to belong to Vaill. and revived the genera transfer from Osteospermum (Nordenstam and Källersjö, 2009). Less. and Less., both of which had been Norlindh (1943) established the genus Chrysanthemoides Tourn. ex Medik. for two species of Osteospermum that had fleshy drupe-like fruits. However, fleshy or semi-drupaceous fruits have later also been reported in some species of Osteospermum such as ⁎ Corresponding author. E-mail address: [email protected] (N.P. Barker). O. junceum Berg., O. asperulum (DC.) T. Norl., O. corymbosum 1 Current address: Institute for Systematic Botany, Zollikerstrasse 107, CH L. and O. triquetrum L. f. (Wood and Nordenstam, 2003). The 8008 Zürich, Switzerland. evolution of drupes or fruits with a fleshy exocarp is extremely rare

0254-6299/$ - see front matter © 2009 SAAB. Published by Elsevier B.V. All rights reserved. doi:10.1016/j.sajb.2009.05.006 N.P. Barker et al. / South African Journal of Botany 75 (2009) 560–572 561

Table 1 Comparison of the species and infraspecific taxon names recognized by Norlindh (1943) and the unpublished entities recognized by Griffioen (1995), with key characteristics of each entity provided. Norlindh (1943) Griffioen (1995) Morphological characteristics Habit Involucral scales Drupes Ecology C. monilifera “C. monilifera Shrub, 1.5 m in Elliptic or slightly Inner involucral Fleshy, orange-red Foundindisturbedsites subsp. subsp. height×2.0 m in obovate, 0.4–0.6 mm scales lanceolate, when mature, on margins of monilifera monilifera” diameter, stems thick, 7 –30 mm× outer involucral globose to climax vegetation in a tanniferous, 20–50 mm, glabrous, scales linear, sub-globose, range of soils pubescence absent. tannins absent, margins slightly pubescent, 5–8 mm in diameter, including TMS, toothed, petioles 4–7 mm long. length: breadth ratio limestone and 6–2-mm long. 1.1–1.2, Achene loams in the SW Cape, not ridged. on mountain slopes from Piketberg in the north to Worcester and Hermanus in the east. “C. monilifera Small to large erect Obovate (6–40 mm× Inner and outer Fleshy, purple-black Coastal sand dunes, subsp. floribunda bush, 1–3.5 m in 10–50 mm), margins involucral scales when mature, limestones and along form 1” height×1–6m spinescent or scalloped, narrowly ovoid, 4–6mm roadsides from diameter, stems lamina leathery or fleshy. lanceolate to ovate, long, achene ridged. Langebaan to Knysna. tanniferous, Young leaves and capitula glabrous of pubescence absent clustered at branch pubescent, from young tissues. terminals, covered with 4–7mmlong. loose pubescence. “C. monilifera Small to large erect Narrowly to broadly Broadly ovate, Ovoid, 4–6 mm long, Costal dunes from subsp. floribunda bush, 1–3.5 m in elliptic, size variable covered with a achene not ridged. Knysna to Hiumansdorp, form 2” height×1–6m (15–60 mm×40– loose pubescence. spreading into forest diameter. Stems 100 mm), margins edges in Knysna and tanniferous, scalloped, lamina thin. George area. pubescence absent from young tissues. C. monilifera “C. monilifera subsp. Large bushes or 15–30 mm×25–45 mm, Lanceolate, Ovate, 4–6mm Grahamstown eastwards subsp. pisifera var. pisifera small trees, 2 m in elliptic, glabrous, margins 3–7 mm long, long. into Transkei region, in pisifera form 1” diameter, 2.5 m 3–5 dentate. sparsely fynbos or false fynbos, high. pubescent. shallow but fertile soils at altitudes between 1000 and 1500 m a.s.l. “C. monilifera subsp. Small shrub to Variable, 5–30 mm× Lanceolate, Ovoid, white when Bredasdorp to Uniondale pisifera var. pisifera 1.5 m high. 8–45 mm, elliptic, 4–10 mm, immature, turning and Humansdorp, in form 2” glabrous, leathery, little or no purple-black when disturbed ground, in deeply dentate, pubescence. ripe, 4–6 mm long. sands, clays and rich 2–4 teeth soils of river gullies at per margin. altitudes between 20 and 900 m a.s.l. size increases westwards. “C. monilifera subsp. Stems without Somewhat leathery, Lanceolate, Ovate. Kamiesberg and pisifera var. borealis” tannins, narrowly elliptic, 4–9 mm long, surrounds in pubescence absent. 10–20 mm×25–45 mm, glabrous, Namakwaland. glabrous, margins spinescent. “C. monilifera subsp. Erect bush, Narrowly elliptic, Ovate, partly Ovoid, Montagu, Swellendam, pisifera var. 1–2 m tall, 10–25 mm×35–60 mm, pubescent. purple-black Bredasdorp region in angustifolia” 1 –2m glabrous, margins when ripe. Resnosterbosveld diameter. scalloped. and false fynbos, restricted to moist Bokkeveld shales and sandy acid flats. C. monilifera “C. monilifera Stems pubescent, Elliptic to broadly Inner scales Obovoid, purple- Associated with Protea subsp. subsp. canescens” tanins absent. elliptic, 0.3–0.5 mm slightly ovate, black when mature, and grassland vegetation canescens thick, 10–35 mm× outer scales 2–6 mm long, in the Drakensberg and 25–55 mm, pubescent, lanceolate, achene not ridged. Suikerbosrand. tannins absent, margins pubescent, toothed, petiole 5–7 mm long. 5–20 mm long. (continued on next page) 562 N.P. Barker et al. / South African Journal of Botany 75 (2009) 560–572

Table 1 (continued) Norlindh (1943) Griffioen (1995) Morphological characteristics Habit Leaves Involucral scales Drupes Ecology C. monilifera “C. monilifera subsp. Erect shrub or Elliptic to narrowly Inner scales ovate, Obovoid, purple- Tropical montane subsp. septentrionalis” bush, stems elliptic, 7–25 mm× outer scales black when mature, habitats between 1500 septentrionalis without tannins, 10–40 mm, pubescent or broadly lanceolate, 3–6 mm long, and 2400 m a.s.l. pubescent. glabrous, tannins absent, slightly pubescent, achene ridged. in Zimbabwe, margins entire or 3–7 mm long. Tanzania and Kenya, spinescent, petiole montane grasslands and 4–12 mm long. woodlands, become pubescent above 2000 m a.s.l. prefer dolerite derived soils. C. monilifera “C. monilifera subsp. Scrambling shrub Obovate to broadly so, Inner scales Fleshy, obovoid, Coastal sands from Port subsp. rotundata” or small tree, 15–60 mm×30–70 mm, lanceolate or purple-black Elizabeth to Inharrime rotundata stems tanniferous, glabrous when mature, narrowly ovate, when mature, (Mozambique). glabrous. tanniferous, margins outer scales 3–6 mm long, entire or scalloped, lanceolate, slightly achene ridged. petiole 6–20 mm long. pubescent or glabrous, 2.5–5 mm long. C. monilifera “C. incana subsp. Bushes up to 2.5 m Narrowly elliptic, Inner scales ovate, Ovoid, 4–6mm Scattered, infrequent. subsp. subcanescens” high, spreading, 3 –10 mm×10–32 mm, outer scales long, achene Robertson, Oudtshoorn, subcanescens young tissue margins minutely lanceolate, not ridged. Middelberg, pubescent, stems spinescent or entire, 2.5–7 mm long. Mossel Bay and up to 2 m long, petiole 2–10 mm long. Graaff Reinet region. glabrous. C. incana “C. incana subsp. Low, spreading Narrowly elliptic, Pubescent. Ovoid, 3.5–7mm Clanwilliam, Calvinia incana var. gracilis” bush, stems 3–12 mm×10–30 mm, long. and Hondeklip Bay slender, pubescent. pubescent, somewhat leathery. “C. incana subsp. Prostrate, Broadly obovate, Pubescent, Found only at Cape incana var. hirsuta” spinescent bush, 7–25 mm×10–30 mm, 3–5 mm long. Agulhas on coastal 0.5 m high, leathery, pubescent, sand dunes. stems pubescent. margins entire or dentate. “C. incana subsp. Prostrate Alternate, broadly obovate, Ovate, 2–5mm Ovoid, 5–m7m Coastal sands from incana var. incana” spinescent bushes or 10–35 mm×20–40 mm, long, pubescent. long, achenes Gouritz River mouth small shrubs up pubescence peeling not ridged. to Saldanha and to 1 m tall, 1 –4m or glabrous, margins northwards to in diamater, stems toothed. Vredendal. with peeling pubescence. “C. incana subsp. Prostrate, Alternate, obovate, Ovate, 2–5mm Ovoid, Mainly localised to incana var. spinescent bushes 5–15 mm×10–20 mm, long, pubescent. achenes not Saldanha area, in microphylla” or small shrubs, pubescent, margin ridged. Strandveld, but also 04.–2 m tall, toothed. known from 1 –3 m diamater. Bredasdorp and Stems pubescent. Laingsberg. “C. incana Prostrate shrub 0.5 m Obovate, 4–5mm× Pubescent. Ovoid, Coastal, from Kleinsee subsp. incana high, stems 5–25 mm, leathery, 4.5–7mm to Spencer Bay var. rangei” pubescent, bearing pubescent, margins long. (Namibia). leaves at tips. toothed; clustered at brannh tips. in the Asteraceae, and as far as is known, no other Old World analysis of morphological characters indicated that C. incana members of this have drupes (Norlindh, 1977). could be distinguished from C. monilifera on the basis of Norlindh (1943) recognized two species in Chrysanthe- spinescence, prostrate growth form, and the distribution of moides,viz.C. incana (Burm. f.) T. Norl. and C. monilifera pubescence on the stems, leaves and receptacles in the former (L.) T. Norl., and further divided the latter into five subspecies, species (Griffioen, 1995). These differences are accompanied by some of which were noted to be highly variable. In an various geographic and ecological attributes. Furthermore, unpublished thesis, Griffioen (1995) used morphological data isozyme electrophoresis indicated that the two species are supplemented with isozyme and ecological data to re-assess the distinct, and do not hybridize when co-occurring. Within infraspecific of both species in the genus. A phenetic C. incana, Griffioen (1995) recognized six infraspecific taxa, N.P. Barker et al. / South African Journal of Botany 75 (2009) 560–572 563 and within C. monilifera 10 taxa were recognized (at subspecies, 2000; Syring et al., 2007; Howis et al., in press; Ramdhani et al., variety and form rank; Table 1). in press), and when noted requires careful analysis and Norlindh (1943) noted that the drupes of Chrysanthemoides explanation. Monophyly at infraspecific ranks is likely to be are edible and most likely bird dispersed, and he ascribed the compromised by hybridization and incomplete lineage sorting. presence of the species on St. Helena since before 1839 to bird However, as noted by Holder et al. (2001), distinguishing dispersal. Chrysanthemoides fruits are eaten by a variety of between these two processes is difficult, but a phylogeographic birds in South Africa (Rowan, 1967; Keath et al., 1992; Joffe, approach may enable us to identify monophyletic infraspecific 2001). In St. Helena dispersal by the introduced Indian Myna taxa. If so, then we suggest that this be viewed as evidence has been reported (Ashmole and Ashmole, 2000). In Australia, favouring their recognition as valid infraspecific taxa, and where Chrysanthemoides is an invasive weed, various mam- certainly as “Evolutionary Significant Units” (ESUs sensu mals have been recorded as dispersers as well as emus and Ryder, 1986; cf. Fraser and Bernatchez, 2001). flying frugivorous birds (Weiss, 1986; Meek, 1998). Man has Two widely utilised chloroplast spacer regions were selected thus also acted as a dispersal agent, and both C. monilifera for this investigation: the psbA-trnH spacer, and the trnL intron subsp. rotundata (DC.) T. Norl. and C. monilifera subsp. mo- in conjunction with the associated trnL-trnF spacer (hereafter nilifera are legislated as “weeds of national significance” in termed the trnL-F region). The psbA-trnH region is being Australia (http://www.weeds.org.au/WoNS/bitoubush/), and increasingly used in phylogenetic studies at the intrageneric there has been over two decades of biocontrol research in level (Gielly et al., 1996; Sang et al., 1997; Kim et al., 1999; Australia on C. monilifera (Downie et al., 2007). This species is Chandler et al., 2001; Pelser et al., 2003; McKenzie et al., 2006; also a problem in New Zealand (Roy et al., 2004). The study of McKenzie and Barker, 2008) as well as the intraspecific level these taxa in their native region is thus of vital importance if the (Štorchová and Olson, 2004; Yamashiro et al., 2004). The trnL- spread of these species as weeds is to be controlled (Scott, trnF intergenic spacer region is one of the most commonly used 1996). However, the fact that these taxa are able to establish non-coding regions of cpDNA in phylogenetic studies at the easily has resulted in their use in rehabilitation efforts following intrageneric and species level, and has occasionally been found mining activities (Hälbich, 2003) and to stabilize coastal dunes to be sufficiently variable for use below the species level in urban areas of South Africa (Nichols, 1996). (Barker et al., 2005). Chrysanthemoides is distributed across a number of biomes While the use of the Internal Transcribed Spacer (ITS) region and vegetation types in southern Africa, ranging from the for phylogenetic purposes is considered controversial by some Fynbos of the South Western Cape, to the montane Grassland of (see for example Alvarez and Wendel, 2003; Bailey et al., 2003; the Drakensberg, Chimanimani and mountains of eastern Africa Small et al., 2004 for critiques), it is still the most tractable (Griffioen, 1995). Of the two recognized species, C. incana is nuclear region for molecular systematics at the species and mostly restricted to the South Western Cape but extends mainly genus level (e.g. Feliner and Rosselló, 2007; Mort et al., 2007), along the coast northwards to Namaqualand and southern and has been widely used in phylogenetics of many groups of Namibia (Angra Pequena) and along the southern coast to about Asteraceae. Baldwin (1993) was probably the first to notice Cape Agulhas. C. monilifera is more widespread, with intraspecific variability in ITS sequences in the Asteraceae, and C. monilifera subsp. canescens (DC.) T. Norl. extending from this region has been used in phylogeographic studies in the the Eastern Cape mountains and the Drakensberg into northern Asteraceae (for example Comes and Abbott 2001; Simurda Transvaal, and C. monilifera subsp. septentrionalis T. Norl. et al., 2005; Zachariades et al., 2004; Pelser et al., 2007). distributed in the montane regions of Tropical East Africa from However, at least some of these studies use ITS data in Zimbabwe north to Tanzania (Norlindh, 1943). conjunction with additional data such as AFLP's or plastid Thus, while Chrysanthemoides in current taxonomy comprises DNA data so as to address the potentially problematic issues or only two species, there is ample evidence (published and paralogy and lineage sorting. unpublished) of considerable variation within each species. Here we use both nuclear and chloroplast DNA sequence data to 2. Methods undertake a phylogenetic analysis of the genus to test not only the monophyly of the genus, but also to assess if there is any genetic 2.1. Sampling, DNA extraction, amplification and sequencing evidence to support the recognition of the infraspecific taxa as recognized by either or both Norlindh (1943) and Griffioen (1995), Multiple samples from each of the species of Chrysanthemoides especially within C. monilifera. We emphasise that while we are were obtained, and in particular we focused on ensuring coverage using the taxa as recognized by Griffioen, they are not validly of the various subspecific entities within C. monilifera as published, and hence the names in Table 1 and the text that follows recognized by Norlindh (1943) and Griffioen (1995).Where appear in quotation marks to indicate this status. possible, we collected material at sites that had been sampled for In order to test taxon monophyly, we adopt a multiple Griffioen's phenetic analysis (Griffioen, 1995). All samples were exemplar approach, including two or more specimens repre- identified to subspecies, variety, and in some cases form, using sentative of each taxonomic entity. This approach is important, Griffioen's (1995) key. A total of 35 Chrysanthemoides samples as some molecular studies have indicated that species non- were used (Table 2). Additional sequence data for several species of monophyly (as determined by molecular data) can be quite Osteospermum (sections Homocarpa T. Norl. and Coriacea common (e.g. Crisp and Chandler, 1996; Ohsako and Ohnishi, T. Norl.) were selected in order to test generic monophyly. Only 564 N.P. Barker et al. / South African Journal of Botany 75 (2009) 560–572

Table 2 Taxon, voucher, locality and GenBank details of samples used in this study. Taxon (sensu Griffioen 1995) Voucher Locality GenBank numbers trnL-trnF psbA-trnH ITS 1 “C. incana subsp. incana” AW 20 WC: Darling FJ861600 FJ861556 FJ861513 2 “C. incana subsp. incana” SH 75 Ex. Cult. Kirstenbosch Botanic Gardens FJ861601 FJ861557 FJ861512 (origin unknown) 3 “C. incana subsp. incana” SH 77 WC: North of Bloubergstrand FJ861602 FJ861558 FJ861511 4 “C. monilifera subsp. canescens” NPB 1818 KZN: Bushmans Nek FJ861591 FJ861547 FJ861499 5 “C. monilifera subsp. canescens” NPB 1820 KZN: Sani pass FJ861570 FJ861526 FJ861498 6 “C. monilifera subsp. floribunda form 1” SH 58 WC: 54 km west of Albertinia FJ861577 FJ861533 FJ861502 7 “C. monilifera subsp. floribunda form 1” SH 68 WC: Struisbaai FJ861596 FJ861552 FJ861510 8 “C. monilifera subsp. floribunda form 1” SH 70 WC: Baardskeerdersbos junction FJ861587 FJ861543 FJ861497 to Gansbaai 9 “C. monilifera subsp. floribunda form 1” SH 71 WC:13 km east of Gansbaai FJ861586 FJ861542 FJ861496 10 “C. monilifera subsp. floribunda form 1” SH 94 WC: Swartberg pass FJ861583 FJ861539 FJ861503 11 “C. monilifera subsp. floribunda form 1” SH 105 WC: De Hoop FJ861598 FJ861554 FJ861509 12 “ C. monilifera subsp. floribunda form 1” NPB 1882 EC: Grahamstown FJ861590 FJ861546 FJ861495 13 “C. monilifera subsp. floribunda form 1” SH 57 WC: 19.5 km west of Albertinia FJ861576 FJ861532 FJ861508 14 “C. monilifera subsp. floribunda form 2”⁎ SH 50 WC: Knysna FJ861574 FJ861530 FJ861501 15 “C. monilifera subsp. floribunda form 2”⁎ SH 52 WC: Outeniqua Pass FJ861575 FJ861531 FJ861507 16 “C. monilifera subsp. monilifera” SH 88 WC: Gydo pass FJ861579 FJ861535 FJ861506 17 “C. monilifera subsp. monilifera” SH 62 WC: Between Swellendam FJ861588 FJ861544 FJ861494 and Riviersonderend 18 “C. monilifera subsp. monilifera” SH 86 WC: Piketberg Pass FJ861585 FJ861541 FJ861505 19 “C. monilifera subsp. pisifera var. angustifolia” SH 72 WC: Hermanus FJ861578 FJ861534 FJ861500 20 “C. monilifera subsp. pisifera var. angustifolia” SH 90 WC: Barrydale, Tradouw's pass FJ861584 FJ861540 FJ861504 21 “C. monilifera subsp. pisifera var. pisifera form 1” SH 111 WC: Joubertina FJ861572 FJ861528 FJ861489 22 “C. monilifera subsp. pisifera var. pisifera form 1” NPB 1831 EC: Van Stadens wild flower reserve FJ861571 FJ861527 FJ861488 23 “C. monilifera subsp. pisifera var. pisifera form 1” SH 43 EC: Kareedouw FJ861597 FJ861553 FJ861486 24 “C. monilifera subsp. pisifera var. pisifera form 1” SR 173 EC: Jeffrey's Bay FJ861580 FJ861536 FJ861492 25 “C. monilifera subsp. pisifera var. pisifera form 1” SH 40 EC: Langkloof road (N2/R62 junction) FJ861573 FJ861529 FJ861491 26 “C. monilifera subsp. pisifera var. pisifera form 2” SH 60 WC: Tradouw Pass FJ861589 FJ861545 FJ861514 27 “C. monilifera subsp. pisifera var. pisifera form 2”⁎ SH 100 WC: Nature's Valley FJ861594 FJ861550 FJ861484 28 “C. monilifera subsp. pisifera var. pisifera form 2”⁎ SR 189 EC: Kareedouw FJ861582 FJ861538 FJ861481 29 “C. monilifera subsp. pisifera var. pisifera form 2”⁎ SR 178 WC: Storm's river FJ861581 FJ861537 FJ861490 30 “C. monilifera subsp. rotundata” SH “Pe2” EC: Port Elizabeth FJ861593 FJ861549 FJ861482 31 “C. monilifera subsp. rotundata” SH “RKOS” EC: Kenton-on-Sea FJ861595 FJ861551 FJ861485 32 “C. monilifera subsp. rotundata” CP 480 KZN: Kozi mouth FJ861592 FJ861548 FJ861483 33 “C. monilifera subsp. rotundata” CP 491 KZN: Umkomaas FJ861568 FJ861524 FJ861480 34 “C. monilifera subsp. rotundata” NPB 1832 EC: Sardinia Bay FJ861599 FJ861555 FJ861487 “C. monilifera subsp. septentrionalis” M&P 452 Tanzania FJ861569 FJ861525 FJ861493 Osteospermum aciphyllum AW 409 (S) WC: Swartberg FJ861566 FJ861522 FJ861478 Osteospermum asperulum AW 435 (S) WC: Swartberg Pass FJ861567 FJ861523 FJ861477 Osteospermum ciliatum AW 428 (S) WC: Greyton FJ861564 FJ861520 FJ861476 Osteospermum corymbosum AW 448 (S) WC: Grootvadersbosch FJ861563 FJ861519 FJ861475 Osteospermum junceum AW 411 (S) WC: Jonkershoek FJ861560 FJ861516 FJ861479 Osteospermum pyrifolium AW 447 (S) WC: Grootvadersbosch FJ861562 FJ861518 FJ861474 Osteospermum subulatum AW 458 (S) WC: Agulhas FJ861561 FJ861517 FJ861473 Osteospermum triquetrum AW 454 (S) WC: Langeberg FJ861565 FJ861521 FJ861472 Tripteris microcarpa A. and B. Strid NC: Garies. FJ861559 FJ861515 FJ861471 37658 (S) Unless otherwise indicated, all vouchers are housed in the Selmar Schonland herbarium (GRA). An ⁎ indicates those samples which are incongruent between the ITS and cpDNA phyologenies (i.e. they swap positions between Clades 1 and 2). Numbers in the left column correspond to localities numbered on the map provided in Fig. 2. Key to collectors: SH = Seranne Howis, NPB = Nigel Barker, AW = Alan Wood, SR = Syd Ramdhani, CP = Craig Peter, M&P = Massawe and Phillipson. Key to province codes in locality column: EC = Eastern Cape; KZN = Kwa Zulu Natal; NC = Northern Cape; WC = Western Cape. species from these sections were considered, as other studies have The more distantly related Tripteris microcarpa Harv. was used as indicated a close relationship between Chrysanthemoides and these outgroup. sections of Osteospermum (Wood and Nordenstam, 2003; Leaf samples were collected into silica gel according to the Nordenstam and Källersjö, 2009). It should be noted that method of Chase and Hills (1991), and the DNA was extracted O. subulatum DC. (of sect. Trialata in Norlindh, 1943)and using a modified CTAB DNA extraction protocol (Doyle and O. triquetrum (unassigned to section in Norlindh, 1943)were Doyle, 1987). The ITS region was amplified using “ITS-5” transferredtosect.Homocarpa by Wood and Nordenstam (2003). (White et al., 1990) and a modified “ITS-4” primer (“ITS- N.P. Barker et al. / South African Journal of Botany 75 (2009) 560–572 565

Chrys-4”;5′-TCCTCCGCTTATGGATATGC-3′). The psbA- the TBR branch swapping algorithm. A strict consensus tree trnH spacer was amplified using the primers “psbA” and “trnH” was then calculated from the set of equally parsimonious trees. (Sang et al., 1997), and the trnL-F region amplified using the A FULL HEURISTIC Bootstrap analysis was conducted on primers “c” and “f” (Taberlet et al., 1991). each data set using 1000 replicates, with MAXTREES set to PCR amplifications were carried out using either a Thermo- 2000. Hybaid PCRSprint Temperature Cycling System or a Corbett Research PC-960G Microplate Gradient Thermal Cycler, with 3. Results 35-40 cycles of amplification. Successful PCR amplification was confirmed by electrophoresing the PCR products on a 1% MrModelTest identified the best DNA substitution model for agarose gel. PCR products were cleaned using either the the ITS data as the HKY + G (Hasegawa et al., 1985) model QIAGEN QIAquick PCR purification kit or the PROMEGA with variable sites assumed to follow a discrete gamma Wizard SV Gel and PCR purification kit and resuspended in distribution. For the cpDNA data, the best DNA substitution 30 µl of dH2O before being sequenced using the ABI prism model was identified as GTR+G. The Bayesian Inference BigDye Terminator v3.0 or v3.1 Ready Reaction Cycle consensus topologies for the ITS and cpDNA data sets are sequencing kit (Applied Biosystems) according to the manufac- presented in Fig. 1, with the parsimony bootstrap values also turer's instructions. indicated. The parsimony analyses produce poorly resolved The ITS PCR product was sequenced in both directions with the consensus trees, and the nodes that are retained in common with primers “ITS-1” (White et al., 1990), “ITS-Chrys-4”, “Chromo- the BI trees are indicated on Fig. 1. Table 3 lists the details of the 5.8-R” (5′-GATTCTGCAATTCACACC-3′; Barker et al., 2005), parsimony analysis of each data set. Support for the various and the novel internal primer “Chrysanth-5.8-F” (5′-GACTCTCG- nodes in both trees is generally lacking, if one is to accept that GCAACGGATATC-3′). The psbA-trnH spacer was sequenced any Bayesian posterior probability (PP) valueb0.95 and usingthesameprimersasusedinthePCRprocess,andthetrnL-F bootstrap percentagesb70% is weak (Alfaro and Holder, 2006). region was sequenced using the primers “c”, “d”, “e” and “f” The ITS region is twice as variable and contains approxi- (Taberlet et al., 1991). mately three times as many parsimony-informative sites as the Sequence data was checked and edited using SEQUENCHER cpDNA regions (Table 3). Despite this lower variability, the (Version 3.1.1; Gene Code Corporation). Assembled sequences psbA-trnH region included three synapomorphic insertion– were imported into MACLADE (Version 4.06; Maddison and deletion events (indels) that, when mapped on the plastid tree, Maddison, 2000) and aligned manually. provide additional support for three of the main nodes retrieved (Fig. 1; note the exception of a subsequent loss of one of these 2.2. Phylogenetic analyses insertions in specimen SH100 – “C. monilifera subsp. pisifera var. pisifera form 2”–within Clade 1). However, indel data As the psbA-trnH and trnL-F are both found on the chloroplast were not coded and included in phylogenetic analyses. genome which is inherited as a single unit, these data sets were It must be noted that all three samples of C. incana possessed combined to form what we term the cpDNA data set. Prior to multiple ITS paralogues of different lengths, which meant that not analysis, the incomplete 5′ and 3′ ends of the psbA-trnH and all the data from this region could be used, as the trace files became trnL-F regions were excluded. The ITS nrDNA data set was unreadable at the point where the paralogues diverged. This analysed separately. Data sets were analysed by means of problem was particularly severe in the Wood 20 specimen, where Bayesian inference (BI) and maximum parsimony (MP). As 361 aligned sites were affected, and coded as unknown. In the Bayesian analysis is based on explicit models of DNA evolution, remaining two C. incana samples, this problem affected less than the program MrModelTest (Nylander, 2004) was used to select 10 sites. Reasons for the presence of multiple paralogues in these the model of DNA substitution that best fit the data. The Bayesian samples include hybridization followed by incomplete lineage analysis was run using MrBayes v3.1.2 (Huelsenbeck and sorting, or the presence of pseudogenes. However, as the 5.8S Ronquist, 2001) as follows: four Markov chains, three heated regions of these sequences are highly conserved, it seems unlikely and one cold, were run simultaneously for 2,000,000 generations that the sequences used here are pseudogenes (Razafimandimbison and trees were saved every 100 generations. The starting tree was et al., 2004). The obvious solution to this problem would be to clone random, the branch lengths were saved and the first 4000 trees the PCR product and sequence the different paralogues. While this were discarded as burn-in following a visual inspection of a plot has been done for a larger study based on ITS data only (Howis of the likelihood values to ensure stationarity had been reached et al., in prep.), we used the data obtained from directly sequenced well within this limit. The remaining trees were combined and PCR product in this study. This finding is nonetheless interesting, used to generate a 50% majority rule consensus tree and to and this species needs detailed molecular and morphological determine the PP for each node. investigation. Parsimony analyses were conducted on the two data sets as follows: Using PAUP version 4.0 b10, one hundred random 3.1. Topological comparisons and incongruence input analyses were conducted using the parsimony-informative characters, keeping one tree (TKEEP=1) from each of the It is immediately apparent that the ITS and cpDNA topologies analyses. A heuristic search was then conducted on all the are not congruent, and in some instances this incongruence is well shortest trees in memory (with MAXTREES set 20,000) using supported. The genus Chrysanthemoides is indicated by both 566 N.P. Barker et al. / South African Journal of Botany 75 (2009) 560–572 N.P. Barker et al. / South African Journal of Botany 75 (2009) 560–572 567

Table 3 phyletic. Other taxa such as “C. monilifera subsp.pisiferavar. Details of the data sets and results from the parsimony analyses of each data set. pisifera form 2” are paraphyletic. cpDNA ITS Combined These results indicate some (but not complete) agreement in Alignment length 1136 746 1882 terms of membership of each of the two clades. This incomplete No. variable characters 107 123 230 agreement is further exacerbated by five samples of C. monilifera (9.4%) (16.5%) (12.2%) that swap clades between the two phylogenies. Three samples in No. potentially 20 (1.76%) 47 (6.3%) 67 (3.6%) Clade 1 of the cpDNA phylogeny are placed in Clade 2 of the ITS parsimony-informative characters Number of most parsimonious trees 336 20,000 20,000 phylogeny (indicated by solid dots in Fig. 1), and two samples Tree length 28 75 129 placed in Clade 2 of the cpDNA phylogeny are in Clade 1 of the Consistency index 0.786 0.733 0.597 ITS phylogeny (indicated by squares in Fig. 1). Retention index 0.963 0.885 0.846 Osteospermum ciliatum Berg. and O. aciphyllum are placed near the base of the tree in both phylogenies, while O. subulatum is in an unsupported clade of four Osteospermum taxa collectively datasets to be non-monophyletic, as in each tree there are various sister to Clade 1 in the cpDNA phylogeny, but placed within Clade species of Osteospermum embedded within the Chrysanthemoides 1 of the ITS topology, with moderate to good support (Fig. 1). clade. Furthermore, the two additional lineages, representing the O. junceum and O. asperulum are sister taxa in the ITS topology, Drakensberg and East African taxa (“C. monilifera subsp. canes- and are in turn sister to the “C. monilifera subsp. canescens” clade, cens” and “C. monilifera subsp. septentrionalis”),areplacedina a relationship that lacks support. However, in the cpDNA topology, more basal position, among species of Osteospermum in the ITS O. junceum is sister to part of Clade 1 that comprises samples of analysis. Both data sets resolve “C. incana subsp. incana” as a C. monilifera (with moderate) and O. asperulum is sister to the rest monophyletic species which is embedded within C. monilifera. of Clade 2 (with moderate support at best). The two main clades of samples on the plastid phylogeny O. pyrifolium T. Norl. and O. triquetrum are members of Clade receive at best moderate support (pp=0.93 and BS=84% for 2 of the cpDNA phylogeny, and placed in a well supported clade Clade 2). The ITS analysis shows higher levels of support (pp=0.97) with “C. monilifera subsp. canescens” samples. (pp=0.99, BS=84% for Clade 1 and pp=0.92 for Clade 2, which However, in the ITS phylogeny they placed as members of has no BS support). Furthermore, each of these clades corresponds separate clades (O. pyrifolium is part of Clade 1 and O. triquetrum reasonably well with the Griffioen's taxonomy. Within Clade 1, the is sister to Clade 2). samples of “C. incana subsp. incana” are sister to the clade containing the remaining samples of “C. monilifera subsp. pisifera 3.2. Combined data set var. pisifera form 1” and “C. monilifera subsp. rotundata”.Clade2 comprises samples of “C. monilifera subsp. floribunda form 1”, Wiens (1998) argues that conflicting data sets can be combined, “C. monilifera subsp. monilifera”, “C. monilifera subsp. pisifera and that the accuracy of the recovered topology as being a reflection var. pisifera form 2” and “C. monilifera subsp. pisifera var. an- of the true phylogeny is enhanced. We thus combined the data in an gustifolia”. In addition, the cpDNA places the two montane attempt to enhance the phylogenetic signal in the data. Parsimony subspecies (“C. monilifera subsp. canescens”d an “C. monilifera analysis of the combined data resultedinapoorlyresolvedtree subsp. septentrionalis”) in this clade as well. However, the ITS data (nodes retained in the strict consensus tree are indicated in Fig. 3), a places these taxa in a more basal position, among species result typical of instances where hybridization has occurred and of Osteospermum, but with no support. reduced resolution in parsimony analyses (McDade, 1992). How- In the cpDNA phylogeny, none of the infraspecific taxa apart ever, the results of the BI analysis are encouraging, in that most from “C. monilifera subsp. canescens” are resolved as mono- samples of C. monilifera formed monophyletic lineages corre- phyletic, the two-sample clade of which receives pp=1.0 and sponding to their taxonomic identity (Fig. 3), although BI posterior BS=88. In contrast, the ITS data resolves samples of four of the probability support for most nodes was still lacking. taxa as recognized by Griffioen (1995) as monophyletic, but (with one exception) with insignificant support. These are “C. monilifera 4. Discussion subsp.pisiferavar.pisiferaform 1” (Clade 1), “C. monilifera subsp. floribunda form 2” (Clade 1), “C. monilifera subsp. monilifera” 4.1. Support for infraspecific taxa within C. monilifera (Clade 2, with good PP and BS support) and “C. monilifera subsp. pisifera var. angustifolia” (Clade 2). With one exception, samples While our sample size here is limited, the ITS sequence data of “C. monilifera subsp. rotundata” (Clade 1) are also mono- does provide some evidence to support the recognition of the

Fig. 1. Bayesian consensus trees from analysis of cpDNA (right) and nrDNA (left) data sets of 31 Chrysanthemoides samples and eight Osteospermum species. Numbers above the branches indicate posterior probabilities, and numbers below the branches are parsimony bootstrap values. Those nodes that are retained in the parsimony consensus trees are indicated as thick branches. The vertical bars indicate the two main clades of Chrysanthemoides samples discussed in the text. Osteospermum species are in bold, and the samples outlined in solid or dashed lines are the Afromontane subspecies of C. monilifera. The solid circles and squares indicate samples that are incongruent within the two main clades of Chrysanthemoides. The grey vertical bar at the base of Clade 1 in the cpDNA tree represents a 10 base pair insertion in the psbA-trnH data (with a subsequent loss in one specimen, indicated by the grey circle), the open bar at the base of Clade two in the cpDNA tree represent a 4-base pair deletion in the psbA-trnH data, and the thin black bar at the node above this represents a 6-base pair deletion in the psbA-trnH data. Key to taxon abbreviations: C. m = C. monilifera;C.m.p.=C. monilifera subsp. pisifera;C.i.=C. incana. 568 N.P. Barker et al. / South African Journal of Botany 75 (2009) 560–572 multitude of infraspecific taxa recognized by Griffioen (1995), are placed within a clade that comprises most samples from Clade although few nodes receive good support. It is unfortunate that 1 with good support (pp=0.95), reflecting the results of cpDNA the plastid data is too conservative to provide independent analysis. The other two samples are placed basal to a clade that support for these lineages. It should also be borne in mind that comprises [Clade 1; O. subulatum], a position that approximates gene trees do not equate to species trees (Doyle, 1992), and that the ITS result. The geographic patterns based on the analysis of the infraspecific relationships may be reticulate rather than combined data are thus even more striking, with the enlarged Clade hierarchical in nature. If so, then our results represent efforts 1 (which includes O. pyrifolium and O. subulatum and having a to place square pegs (reticulating taxonomic entities) into round pp=0.98) now assuming a distinct Eastern distribution, while the holes (a hierarchical representation of relationships). remaining samples form a Western lineage (Figs. 2 and 3). Our results (especially the combined analysis; Fig. 3)do however indicate some genetic support for the infraspecific entities 4.3. Relationship of Chrysanthemoides to Osteospermum recognized by Griffioen (1995) on the basis of morphology (which is a phenotypic representation of the nuclear genotype). This lends On the basis of both ITS and cpDNA gene trees as well as the weight to their validity, and we thus feel that it is important to combined analysis, both Osteospermum and Chrysanthemoides encourage their use by the end-users of taxonomies. We thus are paraphyletic, bearing in mind that many placements are not present the key features of these entities in Table 1.Thisresult well supported. If these gene trees are taken to be straightforward highlights the value of careful phenetic and ecological studies in reflections of evolutionary history, the simplest nomenclatural and variable species (Griffioen, 1995), which should accompany taxonomic solution to this problem is to subsume all Chry- molecular phylogeographic studies, as all sources of data need to santhemoides taxa as well as Osteospermum sect. Coriacea be considered in assessing species limits and the recognition of (O. junceum)intoOsteospermum sect. Homocarpa, which would ESUs within species. The fact that other of Griffioen's (1995) then be monophyletic. This would mean that the defining taxonomic entities are not monophyletic should not be viewed morphological characteristic of Chrysanthemoides (drupaceous negatively — it is entirely possible that paraphyletic taxa (such as fruit) would then have to be interpreted as having originated “C. monilifera subsp. rotundata”, “C. monilifera subsp. pisifera several times, followed by losses/reversals in those species of form 2” and “C. monilifera subsp. floribunda form 1”)represent Osteospermum shown to be embedded within the lineages of surviving ancestral lineages, out of which some of the mono- Chrysanthemoides. The number of gains and losses depends on phyletic lineages have recently evolved, and/or lineage sorting has the phylogeny used: cpDNA or nrDNA. However, it would be not yet reached completion in these taxa. Another possibility is that most profitable to re-examine the fruit morphology and its this is a side effect of relatively rare hybridization. This should, in ontogeny for all species in Osteospermum section Homocarpa,as particular, affect the plastid phylogenies, as they tend to give more there is recent evidence (Wood and Nordenstam, 2003) that fruits categorical results, and so can be expected to be categorically with a fleshy exocarp are not restricted to Chrysanthemoides as wrong. claimed by Norlindh (1943) and that exocarp features (including colour) may have evolved independently and repeatedly as 4.2. Correlation to distribution adaptations to different dispersal strategies (bird vs. ant dispersal; cf. Wood and Nordenstam, 2003). When the geographic distribution of the Chrysanthemoides samples in each clade identified in Fig. 1 is mapped, it becomes 4.4. The role of hybridization in Osteospermum apparent that Clades 1 and 2 are correlated to geographic and Chrysanthemoides distribution: Clade 1 comprises samples from the eastern portion of the distribution range, and Clade 2 the western Irrespective of proposed taxonomic and nomenclatural (Fig. 2). The western clade (Clade 2) comprises specimens from changes required to preserve generic monophyly, the issue of the predominantly winter rainfall region, whereas Clade 1 conflicting species-level relationships within this assemblage covers the bi-seasonal (but predominantly summer) rainfall remains. This could be caused by several processes. Generally, region. Interestingly, the five samples that exchange clade such conflict is attributed to hybridization and/or incomplete membership between the ITS and cpDNA phylogenies are all lineage sorting as well as the questionable utility of ITS as a from the geographic region intermediate between the main suitable nuclear marker (see for example Vriesendorp and Bakker, distribution areas of the two main clades; the Tsitsikamma – 2005; Mort et al., 2008). Given the numerous instances of Nature's Valley – Oudtshoorn area, suggesting that hybridiza- incongruence in this group, we would have to invoke multiple tion between the two main clades in this region cannot be ruled instances of hybridization (and possibly subsequent introgres- out. It is thus possible that gene flow between these lineages has sion), such that the general “Chrysanthemoides” form or mor- taken place, and that the different positions of these five samples phology is retained throughout in disparate lineages, such as in the ITS phylogeny reflects the fixation of one particular set of “C. monilifera subsp. canescens” and “C. monilifera subsp. ro- parental paralogues via concerted evolution — most likely tundata”. Certainly the intermediate geographic distribution of those from the paternal source which are incongruent with the samples that swap clade membership may be considered as maternal cpDNA topology. evidence for past or ongoing hybridization. The report of a natural In the analysis of the combined data, three of the five samples hybrid between Osteospermum potbergense A. R. Wood & that swap clade membership (those indicated in Fig. 1 by circles) B. Nord. and Chrysanthemoides monilifera is noteworthy in this ..Bre ta./SuhArcnJunlo oay7 20)560 (2009) 75 Botany of Journal African South / al. et Barker N.P. – 572

Fig. 2. Map showing sample sites for South African specimens of Chrysanthemoides. Key to shapes: squares indicate localities of C. monilifera samples in Clade 1; circles indicate localities of C. monilifera samples in Clade 2; stars indicate localities of C. monilifera samples which swap clade membership between cpDNA and ITS data sets; triangles indicate localities of C. monilifera subsp. canescens; diamonds indicate localities of samples of C. incana. The grey outline indicates samples that fall into the Eastern clade in the analysis of the combined data set. Numbers within symbols indicate samples as listed in Table 1. 569 570 N.P. Barker et al. / South African Journal of Botany 75 (2009) 560–572

Fig. 3. Bayesian consensus tree from analysis of the combined (cpDNA and nrDNA) data set of 31 Chrysanthemoides samples and eight Osteospermum species. Numbers above the branches indicate posterior probabilities, and numbers below the branches are parsimony bootstrap values. Those nodes that are retained in the parsimony consensus tree are indicated as thick branches. The vertical bars indicate the two main clades of Chrysanthemoides samples discussed in the text. Osteospermum species are in bold, and the samples outlined in solid or dashed lines are the Afromontane subspecies of C. monilifera. The solid circles and squares indicate samples that are incongruent between the cpDNA and ITS data sets. Key to taxon abbreviations: C. m = C. monilifera;C.m.p.=C. monilifera subsp. pisifera;C.i.=C. incana. context (Wood and Nordenstam, 2003). Alternatively, the entire 5. Conclusion Osteospermum section Homocarpa is of very recent origin, and the results here represent incomplete lineage sorting, a scenario Our results indicate the existence of substantial intraspecific that cannot be examined fully until a more complete phylogeny is variation within C. monilifera, corresponding in part to the obtained and subjected to some form of dating analysis. various morphotypes recognized by Griffioen (1995) at different N.P. Barker et al. / South African Journal of Botany 75 (2009) 560–572 571 taxonomic ranks. Unfortunately, the taxonomic recognition and Baldwin, B.G., 1993. Molecular phylogenetics of Calycadenia (Compositae) disentangling of these entities, along with the relationships of based on ITS sequences of nuclear ribosomal DNA. American Journal of Botany 80, 222–238. these to each other and other Osteospermum species is bedeviled Barker, N.P., Von Senger, I., Howis, S., Zachariades, C., Ripley, B.S., 2005. by incongruence between nuclear and plastid markers, suggesting phylogeography based on rDNA ITS sequence data: two examples a history of hybridization events and reticulation or incomplete from the Asteraceae. In: Bakker, F.T., Chatrou, L.W., Gravendeel, B., Pelser, lineage sorting further confused by the distinct possibility of P.B. (Eds.), Plant Species-level Systematics: New Perspectives on Patterns – multiple paralogues of the ITS region. Despite the slight im- and Process. Gantner Verlag, Ruggell, pp. 217 244. Chandler, G.T., Bayer, R.J., Crisp, M.D., 2001. A molecular phylogeny of the provement in the phylogeny obtained when the data are endemic Australian genus Gastrolobium (Fabaceae: Mirbelieae) and allied combined, we refrain from nomenclatural validation of Griffoen's genera using chloroplast and nuclear markers. American Journal of Botany 88, unpublished infraspecific taxa. 1675–1687. Chase, M.W., Hills, H.H., 1991. Silica gel: an ideal material for field preservation of leaf samples for DNA studies. Taxon 40, 215–220. 5.1. Significance to users of taxonomies Comes, H.P., Abbott, R.J., 2001. Molecular phylogeography, reticulation, and lineage sorting in Mediterranean Senecio sect. Senecio (Asteraceae). Evolution Our results are of considerable significance for biocontrol 55, 1943–1962. – scientists working to control Chrysanthemoides in regions Crisp, M.D., Chandler, G.T., 1996. Paraphyletic species. Telopea 6, 813 844. Downie, P.O., Holtkamp, R.H., Ireson, J.E., Kwong, R.M., Swirepek, A.E., where it is invasive. As the current taxonomy does not 2007. A review of the Chrysanthemoides monilifera biological control adequately address the genetic diversity within species of program in Australia: 1987–2005. Plant Protection Quarterly 22, 24–32. Chrysanthemoides, it is essential for biocontrol scientists to Doyle, J.J., 1992. Gene trees and species trees: molecular systematics as one- ascertain which genetic lineage and geographic area the character taxonomy. Systematic Botany 17, 144–163. invasive plants are from e.g. Zachariades et al. (2004) and Doyle, J.J., Doyle, J.L., 1987. A rapid DNA isolation procedure for small quantities of fresh leaf material. Phytochemical Bulletin 19, 11–15. Paterson et al. (2009). Once this is known, it may be possible to Feliner, G.N., Rosselló, J.A., 2007. Better the devil you know? Guidelines for obtain more effective biocontrol organisms from natural insightful utilization of nrDNA ITS in species-level evolutionary studies in populations of that specific genotype, as for example shown plants. Molecular Phylogenetics and Evolution 44, 911–919. for the control of the invasive fern Lygodium microphyllum by Fraser, D.J., Bernatchez, L., 2001. Adaptive evolutionary conservation: towards phytophagous mites (Goolsby et al., 2006). a unified concept for defining conservation units. Molecular Ecology 10, 2741–2752. Our findings are also relevant to the horticultural and Gielly, L., Yuan, Y.-M., Kupper, P., Taberlet, P., 1996. Phylogenetic use of landscaping industry, as C. monilifera is used for a range of noncoding regions in the genus Gentiana L.: chloroplast trnL(UAA) intron horticultural practices (including landscaping and restoration, versus nuclear ribosomal internal transcribed spacer sequences. Molecular noted above). As we show that genetic diversity is correlated Phylogenetics and Evolution 5, 460–466. with geographic origin, it is clearly important to ensure that only Goolsby, J.A., De Barro, P.J., Makinson, J.R., Pemberton, R.W., Hartley, D.M., Frohlich, D.R., 2006. Matching the origin of an invasive weed for selection locally adapted lineages are used for these activities, otherwise of a herbivore haplotype for a biological control programme. Molecular it is possible that foreign genotypes will be imported and may Ecology 15, 287–297. interbreed with local genotypes. This can have disastrous Griffioen, R.C., 1995. A taxonomic study of Chrysanthemoides Tourn. ex Medik. genetic consequences such as outbreeding depression and Compositae). MSc Thesis, University of Cape Town. genetic swamping (Hufford and Mazer, 2003). Hälbich, T.J.R., 2003. Mine rehabilitation in the arid Succulent Karoo vegetation zone of the South African west coast, Namakwa Sands — case study. Heavy Minerals 2003, Johannesburg, South African Institute of Mining and Metalurgy. Acknowledgements Hasegawa, M., Kishino, H., Yano, T., 1985. Dating the human–ape split by a molecular clock of mitochondrial DNA. Journal of Molecular Evolution 22, 160–174. We thank those collectors who contributed to the samples listed Holder, M.T., Anderson, J.A., Holloway, A.K., 2001. Difficulties in detecting in Table 2. S. Howis acknowledges financial support from the hybridization. Systematic Biology 50, 978–982. National Research Foundation of South Africa (NRF), who also Howis, S., Barker, N.P., Mucina, L., (in press). Globally grown, but poorly known: provided support through Grant Unique Numbers (GUN) species delimitation, relationships and biogeography in Gazania Gaertn. 2042600, 2046932 and 2053645. Two anonymous reviewers are (Asteraceae). Taxon.(http://caliban.ingentaselect.co.uk/fstemp/ a459e2fb6892d0d6a7b7f00922c3f593.pdf). thanked for their insights and suggestions which improved this Huelsenbeck, J.P., Ronquist, F., 2001. MrBayes: Bayesian inference of paper. phylogenetic trees. Bioinformatics 17, 754–755. Hufford, K.M., Mazer, S.J., 2003. Plant ecotypes: genetic differentiation in the age of ecological restoration. Trends in Ecology and Evolution 18, 147–155. References Joffe, P., 2001. Creative Gardening with Indigenous Plants: A South African Guide. Briza, Pretoria. Alfaro, M.E., Holder, M.T., 2006. The posterior and the prior in Bayesian Keath, S., Urban, E.K., Fry, C.H., 1992. The Birds of Africa, vol. 4. Academic phylogenetics. Annual Review of Ecology and Systematics 37, 19–42. Press, London. Alvarez, I., Wendel, J.F., 2003. Ribosomal ITS sequences and plant phylogenetic Kim, S.-C., Crawford, D.J., Jansen, R.K., Arnoldo, S.-G., 1999. The use of a non- inference. Molecular Phylogenetics and Evolution 29, 417–434. coding region of chloroplast DNA in phylogenetic studies of the subtribe Ashmole, P., Ashmole, M., 2000. St. Helena and Ascension Island: A Natural Sonchinae (Asteraceae: Lactuceae). Plant Systematics and Evolution 215, 85–99. History. Anthony Nelson, Oswestry. Maddison, W.P., Maddison, D.R., 2000. MacClade 4: Analysis of Phylogeny Bailey, C.D., Carr, T.G., Harris, S.A., Hughes, C.E., 2003. Characterization of and Character Evolution. Sinauer Associates, Sunderland. angiosperm nrDNA polymorphism, paralogy and pseudogenes. Molecular McDade, L.A., 1992. Hybrids and phylogenetic systematics II. The impact of Phylogenetics and Evolution 29, 435–455. hybrids on cladistic analysis. Evolution 46, 1329–1346. 572 N.P. Barker et al. / South African Journal of Botany 75 (2009) 560–572

McKenzie, R.J., Barker, N.P., 2008. Radiation of southern African daisies: Razafimandimbison, S.G., Kellogg, E.A., Bremer, B., 2004. Recent origin and biogeographic inferences for subtribe Arctotidinae (Asteraceae, Arctotideae). phylogenetic utility of divergent ITS putative pseudogenes: a case study Molecular Phylogenetics and Evolution 49, 1–16. from Naucleeae (Rubiaceae). Systematic Biology 53, 177–192. McKenzie, R.J., Barker, N.P., Muller, E.M., Karis, P.O., 2006. Phylogenetic Rowan, M.K., 1967. A study of the colies of southern Africa. The Ostrich 38, relationships and generic delimitation in subtribe Arctotidinae (Asteraceae: 63–115. Arctotideae) inferred by DNA sequence data from ITS and five chloroplast Roy, B., Popay, I., Champion, P., James, T., Rahman, A., 2004. An illustrated regions. American Journal of Botany 93, 1222–1235. guide to the common weeds of New Zealand. New Zealand Plant Protection Meek, P.D., 1998. ‘Weed seeds and whoopsie daisies’: viability of bitou bush Society. Manaaki-Whenua Press. Chrysanthemoides monilifera seeds in fox (Vulpes vulpes) scats. Plant Protection Ryder, O.A., 1986. Species conservation and systematics: the dilemma of Quarterly 13, 21–26. subspecies. Trends in Ecology and Evolution 1, 9–10. Mort, M.E., Archibald, J.K., Randle, C.P., Levsen, N.D., O'Leary, T.R., Sang, T., Crawford, D.J., Stuessy, T.F., 1997. Chloroplast DNA phylogeny, Topalov, K., Wiegand, C.M., Crawford, D.J., 2007. Inferring phylogeny at reticulate evolution and biogeography of Paeonia (Paeoniaceae). American low taxonomic levels: utility of rapidly evolving cpDNA and nuclear ITS Journal of Botany 84, 1120–1136. loci. American Journal of Botany 94, 173–183. Scott, J.K., 1996. Population ecology of Chrysanthemoides monilifera in South Mort, M.E., Randle, C.P., Kimball, R.T., Mesfin Tadesse, Crawford, D.J., 2008. Africa: implications for its control in Australia. Journal of Applied Ecology Phylogeny of Coreopsideae (Asteraceae) inferred from nuclear and plastid 33, 1496–1508. DNA sequences. Taxon 57, 109–120. Simurda, M.C., Marshall, D.C., Knox, J.S., 2005. Phylogeography of the narrow Nichols, G.R., 1996. Preliminary observations on managing and reclaiming frontal endemic, virginicum (Asteraceae), based upon ITS sequence dunes within the Durban municipal area. Landscape and Urban Planning 34, comparisons. Systematic Botany 30, 887–898. 383–388. Small, R.L., Cronn, R.C., Wendel, J.F., 2004. Use of nuclear genes for phylogeny Nordenstam, B., 1994. Tribe Calenduleae. In: Bremer, K. (Ed.), Asteraceae: Cladistics reconstruction in plants. Australian Systematic Botany 17, 145–170. and Classification. Timber Press, Portland, Oregon, USA, pp. 365–376. Štorchová, H., Olson, M.S., 2004. Comparison between mitochondrial and Nordenstam, B., 1996. Recent revision of and Calenduleae chloroplast DNA variation in the native range of Silene vulgaris. Molecular systematics. In: Hind, D.J.N., Beentje, H.J. (Eds.), Compositae: Systematics. Ecology 13, 2909–2919. Proceedings of the International Compositae Conference, Kew. Royal Syring, J., Farrell, K., Businský, R., Cronn, R., Liston, A., 2007. Widespread Botanic Gardens, Kew, UK, pp. 591–596. genealogical nonmonophyly in species of Pinus subgenus Strobus.Systematic Nordenstam, B., 2006. Generic revisions in the tribe Calenduleae (Compositae). Biology 56, 163–181. Compositae Newsletter 44, 38–49. Taberlet, P., Gielly, L., Pautou, G., Bouvet, J., 1991. Universal primers for Nordenstam, B., 2007. Tribe Calenduleae Cass. In: Kadereit, J.W., Jeffrey, C. amplification of three non-coding regions of chloroplast DNA. Plant (Eds.), . Vol. 8 of Kubitzki, K. (Ed.), The Families and Genera of Molecular Biology 17, 1105–1109. Vascular Plants. Springer, Berlin, pp. 241–245. Vriesendorp, B., Bakker, F.T., 2005. Reconstructing patterns of reticulate Nordenstam, B., Källersjö, M., 2009. Calenduleae. In: Funk, V.A., Susanna, A., evolution in angiosperms: what can we do? Taxon 54, 593–604. Stuessy, T., Bayer, R.B. (Eds.), Systematics, Evolution and Biogeography of Weiss, P.W., 1986. The biology of Australian weeds 14. Chrysanthemoides the Compositae. IAPT, Vienna, pp. 527–538. monilifera (L.) T. Norl. The Journal of the Australian Institute of Agricultural Nordenstam, B., Källersjö, M., Eldenäs, P., 2006. Nephrotheca, a new Science 52, 127–134. monotypic genus of the Compositae-Calenduleae from the southwestern White, T.J., Bruns, T., Lee, S., Taylor, J., 1990. Amplification and direct Cape Province. Compositae Newsletter 44, 32–37. sequencing of fungal ribosomal RNA genes for phylogenetics. In: Innis, M.A., Norlindh, T., 1943. Studies in the Calenduleae I. Monograph of the genera Di- Gelfand, D.H., Sninsky, J.J., White, T.J. (Eds.), PCR Protocols: A Guide to morphotheca, Castalis, Osteospermum, Gibbaria and Chrysanthemoides. Methods and Applications. Academic Press, San Diego, California, USA, Lund. Carl Blom Boktryckeri, CWK Gleerup, p. 432. pp. 315–324. Norlindh, T., 1977. Calenduleae — systematic review. In: Heywood, V.H., Wiens, J.J., 1998. Combining data sets with different phylogenetic histories. Harborne, J.B., Turner, B.L. (Eds.), The Biology and Chemistry of the Systematic Biology 47, 568–581. Compositae, vol. 2. Academic Press, London, pp. 961–988. Wood, A.R., Nordenstam, B., 2003. An interesting new species of Osteosper- Nylander, J.A.A., 2004. MrModeltest v2. Evolutionary Biology Centre, Uppsala mum (Asteraceae-Calenduleae) from the Western Cape Province, South University. Program distributed by the author. Africa, providing a link to the genus Chrysanthemoides. South African Ohsako, T., Ohnishi, O., 2000. Intra- and interspecific phylogeny of wild Fago- Journal of Botany 69, 572–578. pyrum (Polygonaceae) species based on nucleotide sequences of noncoding Yamashiro, T., Fukuda, T., Yokoyama, J., Maki, M., 2004. Molecular phylogeny regions in chloroplast DNA. American Journal of Botany 87, 573–582. of Vincetoxicum (Apocynaceae-Asclepiadoideae) based on the nucleotide Paterson, I.D., Downie, D.A., Hill, M.P., 2009. Using molecular methods to sequences of cpDNA and nrDNA. Molecular Phylogenetics and Evolution determine the origin of weed populations of Pereskia aculeata in South 31, 689–700. Africa and its relevance to biological control. Biological Control 48, 84–91. Zachariades, C., Von Senger, I., Barker, N.P., 2004. Evidence for a northern Pelser, P.B., Gravendeel, B., Van der Meijden, R., 2003. Phylogeny reconstruction Caribbean origin for the southern African biotype of Chromolaena odorata. in the gap between too little and too much divergence: the closest relatives of In: Day, M.D., McFadyen, R.E. (Eds.), Chromolaena in the Asia-Pacific Senecio jacobaea (Asteraceae) according to DNA sequences and AFLPs. region. Proceedings of the 6th International Workshop on Biological Control Molecular Phylogenetics and Evolution 29, 613–628. and Management of Chromolaena, held in Cairns, Australia, May 6–9, Pelser, P.B., Nordenstam, B., Kadereit, J.W., Watson, L.E., 2007. An ITS 2003, pp. 25–27. phylogeny of tribe Senecioneae (Asteraceae) and a new delimitation of Senecio L. Taxon 56, 1077–1104. Ramdhani, S., Barker, N.P., Baijnath, H., (in press). Molecular data indicates rampant non-monophyly of species in Kniphofia Moench (Asphodelaceae): a consequence of a recent Afromontane radiation? Taxon.

Edited by JC Manning