The architecture of the chloroplast psbA-trnH non-coding region in angiosperms

The psbA-trnH intergenic region is among the most variable regions in the angiosperm chloroplast genome. It is a popular tool for plant population genetics and species level phylogenetics and has been proposed as suitable for DNA barcoding studies. This region contains two parts differing in their evolutionary conservation: 1) the psbA 3′UTR (untranslated region) and 2) the psbA-trnH intergenic non-transcribed spacer. We compared the sequence and RNA secondary structure of the psbA 3′ UTR across angiosperms and found consensus motifs corresponding to the stem portions of the RNA stem-loop structures and a consensus TTAGTGTATA box. The psbA-trnH spacer exhibited patterns that can be explained by the independent evolution of large inversions in the psbA 3′UTR and mutational hot spots in the remaining portion of the psbA-trnH spacer. We conclude that a comparison of chloroplast UTRs across angiosperms offer clues to the identity of putative regulatory elements and information about selective constraints imposed on the chloroplast non-coding regions.     

Phylogeny of Carpinus and subfamily Coryloideae (Betulaceae) based on chloroplast and nuclear ribosomal sequence data

Phylogenetic studies were conducted for Carpinus and the subfamily Coryloideae (Betulaceae) using sequences of the chloroplast matK gene, the trnL-trnF region (trnL intron, and trnL [UAA] 3' exon-trnF [GAA] intergenic spacer) and the psbA-trnH intergenic spacer, and the nuclear ribosomal ITS regions. The combined analyses of the three chloroplast regions suggest that Coryloideae is monophyletic; Ostryopsis is sister to the Carpinus - Ostrya clade; Corylus is monophyletic and sister to the Ostrya - Carpinus - Ostryopsis clade; Ostrya is paraphyletic; and within Carpinus, species of sect. Carpinus from eastern Asia form a monophyletic group, whereas the positions of C. betulus from Europe and C. caroliniana from eastern North America are unresolved within the Carpinus clade. The cpDNA tree generated in this study is largely congruent with the previously published ITS results, but the ITS tree places Carpinus sect. Distegocarpus as sister to the Ostrya - Carpinus sect. Carpinus clade. Future work is needed to examine the relationships within the Ostrya - Carpinus clade, evaluate the generic status of Ostrya, and test the phylogenetic position of Ostryopsis.     

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

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

A molecular phylogeny of the genus Gagea (Liliaceae) in Germany inferred from non-coding chloroplast and nuclear DNA sequences

The systematics of the mainly yellow flowered Gagea species complex (Liliaceae) has long been considered difficult because only a few phenotypic features within this genus and as a result of hypothesized interspecific hybridisation. A molecular phylogenetic study of seven Gagea species (G. bohemica, G. lutea, G. minima, G. pomeranica, G. pratenis, G. spathacea and G. villosa) from Germany has been undertaken, based on plastid DNA sequences (trnL(UAA)-trnF(GAA), psbA-trnH) and on the nuclear ribosomal internal transcribed spacer (ITS). Sequence divergence within the Gagea species ranges up to 15.5% for psbA-trnH, 22.0% for trnL-trnF and 23.7% for ITS (ITS1 + 5.8S rRNA + ITS2). Two subspecies of Gagea bohemica: G. bohemica subsp. saxatilis and G. bohemica subsp. bohemica are characterized by trnL-trnF data and morphological features. Analysis of the ITS region shows that G. pomeranica represents a hybrid of G. pratensis and G. lutea. Lloydia serotina was initially used as an outgroup species, but it was placed within the investigated Gagea species in the psbA-trnH and the trnL-trnF phylogenetic tree.     

Phylogeny and biogeography of Sassafras (Lauraceae) disjunct between eastern Asia and eastern North America

Sassafras (Lauraceae) consists of three species disjunct between eastern Asia (S. tzumu and S. randaiense) and eastern North America (S. albidum). Phylogenetic analysis based on sequences of nuclear ribosomal ITS and three chloroplast non-coding regions (rpl16, trnL-F, and psbA-trnH) showed that Sassafras is monophyletic and that the eastern North American S. albidum is sister to the clade of its two eastern Asian counterparts. Their intercontinental divergence was estimated to be 13.80 ± 2.29−16.69 ± 2.52 million years ago (mya) using the penalized likelihood method with the ITS and three chloroplast markers. Biogeographic analyses combined with fossil evidence suggest that Sassafras has a relict distribution in the Northern Hemisphere without a Gondwanan link. The divergence time of the two eastern Asian species (the continental Chinese Sassafras tzumu and S. randaiense endemic to Taiwan) is estimated to be 0.61 ± 0.75−2.23 ± 0.76 mya. Sassafras randaiense from Taiwan was most likely derived from an ancestor from continental China.     

山麦冬及其近缘种DNA条形码研究

通过分析山麦冬及其近缘种cpDNA trnL-F、psbA-trnH间隔区序列特点,探讨trnL-F、psbA-trnH序列作为山麦冬及其近缘种DNA条形码的潜力。分别对两者进行PCR扩增后,并进行纯化测序。其中山麦冬及其近缘种9个物种18个样品的trnL-F序列长度为355~356 bp,在山麦冬属内序列完全一致,但在沿阶草属中存在特异变异位点,可对麦冬与沿阶草两物种进行鉴别。而psbA-trnH序列长度为543~544 bp,仅在麦冬中存在特异鉴别位点。结果表明,trnL-F、psbA-trnH序列由于进化速率低、保守性强,仅适合麦冬类植物属间鉴别的DNA条形码,而在属下水平的应用有一定的局限。     

History and evolution of alpine plants endemic to the Qinghai-Tibetan Plateau: Aconitum gymnandrum (Ranunculaceae)

How Quaternary climatic oscillations affected range distributions and intraspecific divergence of alpine plants on the Qinghai‐Tibetan Plateau (QTP) remains largely unknown. Here, we report a survey of chloroplast DNA (cpDNA) and nuclear ribosomal internal transcribed spacer (ITS) DNA variation aimed at exploring the phylogeographical history of the QTP alpine endemic Aconitum gymnandrum. We sequenced three cpDNA fragments (rpl20–rps12 intergenic spacer, the trnV intron and psbA‐trnH spacer) and also the nuclear (ITS) region in 245 individuals from 23 populations sampled throughout the species’ range. Two distinct lineages, with eastern and western geographical distributions respectively, were identified from a phylogenetic analysis of ITS sequence variation. Based on a fast substitution rate, these were estimated to have diverged from each other in the early Pleistocene approximately 1.45 Ma. The analysis of cpDNA variation identified nine chlorotypes that clustered into two major clades that were broadly congruent in geographical distribution with the two ITS lineages. The east–west split of cpDNA divergence was supported by an amova which partitioned approximately half of the total variance between these two groups of populations. Analysis of the spatial distribution of chlorotypes showed that each clade was subdivided into two groups of populations such that a total of four population groups existed in the species. It is suggested that these different groups derive from four independent glacial refugia that existed during the Last Glacial Maximum (LGM), and that three of these refugia were located at high altitude on the QTP platform itself at that time. Coalescent simulation of chlorotype genealogies supported both an early Pleistocene origin of the two main cpDNA clades and also the ‘four‐refugia’ hypothesis during the LGM. Two previous phylogeographical studies of QTP alpine plants indicated that such plants retreated to refugia at the eastern/south‐eastern plateau edge during the LGM and/or previous glacial maxima. However, the results for A. gymnandrum suggest that at least some of these cold‐tolerant species may have also survived centrally on the QTP platform throughout the Quaternary.     

Multilocus phylogenetic reconstruction informing polyploid relationships of <Emphasis Type="Italic">Aconitum</Emphasis> subgenus <Emphasis Type="Italic">Lycoctonum</Emphasis> (Ranunculaceae) in China

Polyploidization has long been recognized as one of the most important driving forces of plant evolution. Aconitum subgenus Lycoctonum (Ranunculaceae) has a wide distribution range and well-known background of polyploidy, thereby providing a potentially valuable model to explore polyploid origin and evolutionary history. However, the phylogeny of subg. Lycoctonum has not yet been completely resolved. In the current study, 29 species including diploid, tetraploid and hexaploid species were sampled in subg. Lycoctonum. Using four cpDNA regions (ndhF-trnL, psbA-trnH, psbD-trnT and trnT-L) and two nrDNA regions (internal transcribed spacer, ITS, and external transcribed spacer, ETS), phylogenetic relationship was first reconstructed for the polyploid species within subg. Lycoctonum. In combination with nuclear diversification rate estimation, cpDNA haplotype network, ancestral area reconstruction as well as morphological and karyotypic evidence, potential origin and divergence time were further assessed among the polyploid species. Hybridization was inferred for A. angustius and A. finetianum was suggested to be the potential maternal progenitor, due to their close phylogenetic relationship, highly similar morphologies and overlapping distribution range. Local origin was inferred for tetraploids in the Hengduan Mountains (HDM) with eight groups of chromosomes of four homeologous, which diverged approximately 3.00 Ma in the same period of the orogeny of the HDM. The hexaploid A. apetalum was inferred to suffer from geographical isolation due to the formation of the Qinghai–Tibetan Plateau (QTP) and the HDM. Hybridization and heterogeneous habitats in the HDM were suggested to play an important role in the polyploidization in subg. Lycoctonum.     

Chloroplast DNA-relationship in palaeoaustralLopidium concinnum (Hypopterygiaceae, Musci). An example of stenoevolution in mosses Studies in austral temperate rain forest bryophytes 2

Evolutionary relationship between disjunct populations of the palaeoaustral moss taxonLopidium concinnum (Hypopterygiaceae) from New Zealand and southern South America were studied using non-coding chloroplast DNA sequences. No or only slight changes could be observed within the sequences oftrnT_UGU —trnL_UAA 5exon intergenic spacer,trnL_UAA intron andtrnL_UAA 3exon —trnF_GAA intergenic spacer. This indicates nearly no genetic divergence between extant New Zealand and Chilean populations, i.e. no significant differing pathways of evolution within the 80–60 million years of disrupted areas with interrupted gene flow. Molecular data support the idea of an old Gondwanan relict species of stenoevolutionary character. Ecological data on short-range dispersal strengthen this assessment.     

Comparisons of phylogenetic hypotheses among different data sets in dwarf dandelions (Krigia,Asteraceae): Additional information from internal transcribed spacer sequences of nuclear ribosomal DNA

The internal transcribed spacer (ITS) region of the 18 S–25 S nuclear ribosomal DNA repeat was sequenced from 19 populations of the tribeLactuceae, including all species of dwarf dandelion (Krigia) and five outgroup genera. The incidence of length changes and base substitutions was at least two times higher for ITS 1 than ITS 2. Interspecific sequence divergence withinKrigia averaged 9.62% (1.61%–15.19%) and 4.26% (0%–6.64%) in ITS 1 and ITS 2, respectively. Intergeneric sequence divergence ranged from 15.6% to 44.5% in ITS 1 and from 8.0% to 28.6% in ITS 2. High sequence divergence and homoplasy among genera of tribeLactuceae suggest that the phylogenetic utility of ITS sequence data is limited to interspecific studies or comparisons among closely related genera. Trees generated from ITS sequences are essentially identical to those from restriction site comparisons of the entire nuclear ribosomal (nr) DNA region. The degree of tree resolution differed depending on how gaps were treated in phylogenetic analyses. The ITS trees were congruent with the chloroplast DNA and morphological phylogenies in three major ways: 1) the sister group relationship betweenKrigia andPyrrhopappus; 2) the recognition of two monophyletic sections,Krigia andCymbia, in genusKrigia; and 3) the monophyly of theK. occidentalis-K. cespitosa clade in sect.Cymbia. However, the two nrDNA-based trees are not congruent with morphology/chloroplast DNA-based trees for the interspecific relationships in sect.Krigia. An average of 22.5% incongruence was observed among fourKrigia data sets. The relatively high degree of incongruence among data sets is due primarily to conflict between trees based on nrDNA and morphological/cpDNA data. The incongruence is probably due to the concerted evolution of nrDNA repeating units. The results fromKrigia and theLactuceae suggest that nrDNA data may have limited utility in phylogenetic studies of plants, especially in groups which exhibit high levels of sequence divergence. Our combined phylogenetic analysis as a total evidence shows the least conflict to each of the individual data sets.     

Molecular and morphological analyses of EuropeanAconitum species (Ranunculaceae)

The systematics of the yellow-floweredAconitum lycoctonum species complex (Ranunculaceae) has long been considered difficult because of high morphological variability and hypothesized hybridization. To elucidate the systematic value of these taxa, and to study their phylogeny in a broader context, we studied all species of the subgenusLycoctonum from Europe and the Caucasus Mountains and the European species of the subgen.Aconitum. We sequenced the ITS region and the chloroplast intergenic spacerpsbA-trnH and conducted parsimony analyses. We also measured morphological characters from 19A. lycoctonum populations. The two subgenera in Europe are clearly separated but hybridization between subgenera might have taken place in the past. Sequence variation withinA. lycoctonum and among all taxa of subgenusLycoctonum was very low and probably the result of recent speciation, leading to one morphologically highly variable species in which color morphs could be seen as subspecies. There was high morphological variability within and among populations and the morphological characters have no value as systematic characters.     

Adaptive radiation and genetic differentiation in the woodySonchus alliance (Asteraceae:Sonchinae) in the Canary Islands

The woodySonchus alliance consists of 19 species ofSonchus subg.Dendrosonchus, one species ofSonchus subg.Sonchus and species of five genera (i.e.Babcockia, Sventenia, Taeckholmia, Lactucosonchus, Prenanthes), and is restricted primarily to the archipelago of the Canaries in the Macaronesian phytogeographical region. An enzyme electrophoretic study, including 13 loci, was conducted to assess genetic diversity within and divergence among species of the alliance. Nei's genetic identities (distances) between genera and/or subgenera range from 0.490 (0.714) to 0.980 (0.013), and pairwise comparisons of all populations show relatively high genetic identities, with a mean of 0.804. The high identities further support the genetic cohesiveness of the alliance and its single origin on the Macaronesian islands. Species in the alliance also show about 50%; higher total genetic diversity (H_T) than the mean for other oceanic endemics. There is greater divergence between endemics or species on older islands compared to those on younger islands, which suggests that time is a factor for divergence at allozyme loci. Furthermore, populations on older islands have higher total genetic diversities and lower identities than conspecific populations on younger islands. These results imply early colonization, radiation, and divergence of the woodySonchus alliance on older islands followed by subsequent colonization to younger islands. The taxonomic distribution of alleles in the alliance indicates lineage sorting also played a role in divergence among species. Lineage sorting may also produce nonconcordance with either taxonomic designation or the pattern of variation obtained from other molecular markers such as ITS sequences of nrDNA. Timing for the origin and radiation of the alliance agrees with the estimate based on ITS sequences, and suggests that the early divergence and rapid radiation took place during the Late Tertiary on either Gran Canaria or Tenerife.     

ITS phylogeny of West Asian Heracleum species and related taxa of Umbelliferae–Tordylieae W.D.J.Koch, with notes on evolution of their psbA-trnH sequences

Summary Heracleum is a large and taxonomically complex genus of the Umbelliferae–Tordylieae. The phylogenetic relationships of West Asian Heracleum species and related taxa were explored using data from sequences of the internal transcribed spacer (ITS1 and ITS2) region of nuclear ribosomal DNA. The data set consists of 56 species, of which 47 were analyzed for the first time; it represents all subdivisions of the genus Heracleum, as well as some representatives of Pastinaca complex. Heracleum was shown to be a polyphyletic genus, as its species fall into two different clades, one of which comprises also Symphyoloma and Mandenovia. Section Pubescentia was confirmed, in contrast to the sections Villosa and Heracleum being polyphyletic. A separate position of the section Wendia was supported. H. marashicum was shown to be a member of a clade comprising Pastinaca and related genera. The sequences of chloroplast psbA-trnH intergenic spacer, the region recently proposed for DNA barcoding in plants, were also analyzed for 33 species, representing all principal clades within Heracleum and its relatives. They have been proven to be very similar and not suitable for DNA barcoding in this group. However, some sequence variation was revealed. This variation could be explained by the combination of such evolutionary events as inversion and duplication. It was shown that these events are rather common in Tordylieae and can occur independently in different lineages. The evolution patterns of psbA-trnH spacer are hypothesized.     

Genetic divergence within the monotypic tree genus <Emphasis Type="Italic">Platycarya</Emphasis> (Juglandaceae) and its implications for species’ past dynamics in subtropical China

Subtropical East Asia harbours a large plant diversity that is often attributed to allopatric speciation in this topographically complex region characterized by a relative climate stability. Here, we use observations of Platycarya, a widespread subtropical Asian tree genus, to explore the consequences of past climate stability on species’ evolutionary history in subtropical China. This genus has a controversial taxonomy: while it is now prevailingly treated as monotypic, two species have been originally described, Platycarya strobilacea and P. longipes. Previous information from species distribution models, fossil pollen data and genetic data based on chloroplast DNA (cpDNA) were integrated with newly obtained genetic data from the two putative species. We used both cpDNA (psbA-trnH and trnL-F intergenic spacers, including a partial trnL gene sequence) and nuclear markers. The latter included sequences of the internal transcribed spacer region (ITS1–5.8S–ITS2) of the nuclear ribosomal DNA and random genomic single nucleotide polymorphisms. Using these nuclear genetic markers, we found interspecific genetic divergence fitting with the ‘two species’ scenario and geographically structured intraspecific variation. Using cpDNA markers, we also found geographically structured intraspecific variation. Despite deep inter- and intraspecific genetic divergence, we detected genetic admixture in southwest China. Overall, our findings of genetic divergence within Platycarya support the hypothesis of allopatric speciation. However, episodes of population interconnection were identified, at least in southwest China, suggesting that the genus has had a dynamic population history.     

The phylogenetic relationships of Cynoglottis (Boraginaceae- Boragineae) inferred from ITS, 5.8S and trnL sequences

A molecular phylogenetic analysis of Cynoglottis was performed to evaluate previous hypotheses based on non-molecular evidence concerning the position of this genus within Boraginaceae tribe Boragineae. ITS-5.8S and trnL_UAA sequences from the nuclear and chloroplast non-coding genomes were obtained for four Cynoglottis taxa and selected members of the related genera Anchusa, Anchusella, Gastrocotyle, Brunnera and Pentaglottis. Cynoglottis is monophyletic, but neither trnL nor ITS support a close relationship with Brunnera, unlike previously supposed on morphological grounds. Brunnera is, instead, related to the southwestern European monotypic genus Pentaglottis, with which it forms a basal clade. ITS-5.8S sequences show that Anchusa thessala, a southeastern European annual species of Anchusa subg. Buglossellum, is sister to Cynoglottis and that the two taxa form a clade which also includes the Balkan endemic Gastrocotyle macedonica. Species of Anchusa subg. Anchusa form a separate lineage with high bootstrap support, suggesting that this heterogeneous genus is paraphyletic with respect to Cynoglottis. ITS sequences also discriminate between the Balkan-Apenninic diploid C. barrelieri and the Anatolian tetraploid C. chetikiana, albeit with low support. The molecular results are discussed in the light of karyological, morphological and chorological aspects.This work has been supported by M.I.U.R. 40% 2003 and the University of Firenze.     

Historical biogeography and phylogenetic relationships of the genus Chamaecyparis (Cupressaceae) inferred from chloroplast DNA polymorphism

The genus Chamaecyparis comprises five species and one variety native to Taiwan, Japan, Canada, and USA, which demonstrates a classical eastern Asian, western North American, and eastern North American disjunct distributional pattern. The phylogenetic relationships of the species of Chamaecyparis were inferred by comparing 1130 bp of the combined data set of chloroplast trnV intron and petG-trnP intergenic spacer. The phylogenetic tree shows that Chamaecyparis nootkatensis (Cupressus nootkatensis or Xanthocyparis nootkatensis) is clearly diverged from other Chamaecyparis species. For Chamaecyparis species, C. thyoides is sister to C. formosensis and C. pisifera and these together form a monophyletic group. C. lawsoniana is sister to C. obtusa and C. taiwanensis; and these form another monophyletic group. Homogeneity in evolutionary rates was found among species in these two monophyletic groups. Results indicate the divergent evolution of C. taiwanensis and C. formosensis and molecular evidence in this investigation supports C. taiwanensis as a variety of C. obtusa. Utility of cpDNA intergenic spacer petG-trnP in Chamaecyparis is also discussed. Several biogeographical implications were inferred: (1) at least two divergence events have produced the eastern Asian, and both western and eastern North American disjunct distribution in Chamaecyparis; (2) intercontinental sister species pairs are found in Chamaecyparis; (3) cpDNA divergence between two intercontinental sister pairs of C. thyoides and C. pisifera, and C. lawsoniana and C. obtusa is 2.8% and 1.1%, which suggest an estimated divergence time of 14 and 5.5 million years ago during middle and late Miocene, respectively; (4) cpDNA divergence of two Asian Chamaecyparis groups between C. obtusa and C. taiwanensis, and between C. pisifera and C. formosensis is 0.25% and 0.57%, which suggest an estimated divergence time of 1.3 and 2.9 million years ago during Pleistocene and late Pliocene, respectively; these estimated divergence times suggest a relatively recent migration of Chamaecyparis to Taiwan from the Japanese Archipelago; (5) that climatic deterioration caused the disappearance of Chamaecyparis in continental Asia is probable.     

Phylogeny of Rhododendron section Vireya (Ericaceae) based on two non-coding regions of cpDNA

Section Vireya (Blume) Copel.f. is one of the most morphologically diverse groups of the genus Rhododendron. Vireyas have a unique distribution for the genus, being predominantly found throughout the Malesian Archipelago. The alpha taxonomy of section Vireya is relatively well understood and taxa are easily distinguished from other rhododendrons by their general appearance. Defining characteristics of the section the possession of seeds with tailed appendages at both ends, the twisting of capsule valves after opening, placentas that separate as thread-like structures from the central axis as the capsule opens are all subject to exceptions. Phylogenetic analyses of two cpDNA regions, psbA-trnH and trnT-trnL intergenic spacers are reported. The results of each analysis were generally congruent, with clades relating strongly to geographic areas. Section Vireya is monophyletic with the inclusion of at least one species of section Rhododendron. Only two of the seven subsections currently circumscribed are monophyletic: Malayovireya Sleumer and Siphonovireya Sleumer. ``Euvireya', all subsections excluding the paraphyletic subsection Pseudovireya, is monophyletic and includes two major clades, one restricted to eastern Malesia and the other to the western and middle Malesia.     

The position ofCastrilanthemum Vogt & Oberprieler and the phylogeny of Mediterranean Anthemideae (Compositae) as inferred from nrDNA ITS and cpDNAtrnL/trnF IGS sequence variation

Sequences of the nrDNA internal transcribed spacer (ITS) region and the cpDNAtrnL/trnF intergenic spacer (IGS) region were analysed for 41 Mediterranean and Eurasian representatives of the Anthemideae (Compositae) to ascertain the systematic position of the unispecific genusCastrilanthemum Vogt & Oberprieler and to study the phylogeny of the Anthemideae in the Mediterranean region. Both markers consistently point to the close relationship ofCastrilanthemum with the W Mediterranean generaLeucanthemopsis (Giroux) Heywood andProlongoa Boiss., forming a strongly supported monophyletic group (theLeucanthemopsis-group which also comprises the unispecific genusHymenostemma (Kunze) Willk.). Results also demonstrate that subtribes Achilleinae, Leucantheminae, and Matricariinae sensu Bremer & Humphries are non-monophyletic. Besides results from sequence variation, a marked 5bp-deletion intrnL/trnF IGS divides all these subtribes into more basal subgroups which are related to each other and to the large Eurasian generaTanacetum L. andAnthemis L., and a monophyletic group of closely related and more advanced subgroups which also contain the monophyletic Chrysantheminae sensu Bremer & Humphries. For this second group a W Mediterranean centre of diversification is suggested, however, its sister-group relationships within the basal grade of generic groups remain unclear.