Genetic similarities versus morphological resemblance: Unraveling a polyploid complex in a Mediterranean biodiversity hotspot.

The Balkan Peninsula is recognized as one of the hotspots of biodiversity in Europe. This area has shown since the Last Glacial Maximum appropriate conditions for species diversification and hybridization, which has led to the existence of numerous taxonomically unresolved entities. Here, we focus on the Western Balkans and explore the genetic structure and relationships among species belonging to the V. austriaca - V. orbiculata diploid-polyploid complex, including populations showing intermediate morphologies. A combination of nuclear markers (microsatellites), plastid DNA regions (trnH-psbA, ycf6-psbM) and ploidy level estimations using flow cytometry are employed to assess the genetic structure and evolutionary dynamics of this polyploid complex. To reconstruct the evolutionary history, an approximate Bayesian computation approach is combined with projections of the species distribution models onto the climatic scenarios of the Mid-Holocene (6 ka BP) and Last Glacial Maximum (22 ka BP). Four main groups were found: one well-established entity within the diploid level, V. dalmatica, a second diploid-tetraploid group which corresponds to V. orbiculata, a hexaploid cluster harboring V. austriaca subsp. jacquinii individuals, and an enigmatic tetraploid group. According to the molecular data obtained, this latter cluster represents an allopolyploid cryptic lineage −with V. orbiculata and V. dalmatica as putative parents− morphologically similar to V. orbiculata, but genetically more related to V. austriaca subsp. jacquinii. Veronica dalmatica and this "uncertain tetraploid" group are involved in the formation of the hexaploid taxon V. austriaca subsp. jacquinii, with the possibility of recent gene flow among different cytotypes. The present study supports a scenario of diversification from a diploid common ancestor leading to two different but interrelated lineages. The first one would correspond with the diploid V. orbiculata plus tetraploid individuals of this species arising through allo- and autopolyploidization, and the second one would involve all ploidy levels with allopolyploidization being prevalent.

Divide and conquer! Data-mining tools and sequential multivariate analysis to search for diagnostic morphological characters within a plant polyploid complex (Veronica subsect. Pentasepalae, Plantaginaceae).

This study exhaustively explores leaf features seeking diagnostic characters to aid the classification (assigning cases to groups, i.e. populations to taxa) in a polyploid plant-species complex. A challenging case study was selected: Veronica subsection Pentasepalae, a taxonomically intricate group. The "divide and conquer" approach was implemented-that is, a difficult primary dataset was split into more manageable subsets. Three techniques were explored: two data-mining tools (artificial neural networks and decision trees) and one unsupervised discriminant analysis. However, only the decision trees and discriminant analysis were finally used to select diagnostic traits. A previously established classification hypothesis based on other data sources was used as a starting point. A guided discriminant analysis (i.e. involving manual character selection) was used to produce a grouping scheme fitting this hypothesis so that it could be taken as a reference. Sequential unsupervised multivariate analysis enabled the recognition of all species and infraspecific taxa; however, a suboptimal classification rate was achieved. Decision trees resulted in better classification rates than unsupervised multivariate analysis, but three complete taxa were misidentified (not present in terminal nodes). The variable selection led to a different grouping scheme in the case of decision trees. The resulting groups displayed low misclassification rates when analyzed using artificial neural networks. The decision trees as well as the discriminant analysis are recommended in the search of diagnostic characters. Due to the high sensitivity that artificial neural networks have to the combination of input/output layers, they are proposed as evaluation tools for morphometric studies. The "divide and conquer" principle is a promising strategy, providing success in the present case study.

The challenge of species delimitation in the diploid-polyploid complex Veronica subsection Pentasepalae.

A reliable taxonomic framework and the identification of evolutionary lineages are essential for effective decisions in conservation biodiversity programs. However, phylogenetic reconstruction becomes extremely difficult when polyploidy and hybridization are involved. Veronica subsection Pentasepalae is a diploid-polyploid complex of ca. 20 species with ploidy levels ranging from 2x to 10x. Here, DNA-ploidy level estimations and AFLP fingerprinting were used to determine the evolutionary history, and species boundaries were reviewed in an integrated approach including also previous data (mainly morphology and sequence-based phylogenetic reconstructions). Molecular analyses were performed for 243 individuals from 95 populations, including for the first time all taxa currently recognized within the subsection. Phylogenetic reconstruction identified four main groups corresponding almost completely to the four clusters identified by genetic structure analyses. Multiple autopolyploidization events have occurred in the tetraploid V. satureiifolia giving rise to octoploid entities in central Europe and north of Spain, whereas hybridization is demonstrated to have occurred in several populations from the Balkan Peninsula. Furthermore, our study has established the taxonomic status of taxa, for the most part recovered as monophyletic. Cryptic taxa within the group have been identified, and a new species, Veronica dalmatica, is fully described. This study highlights the implications of polyploidy in species delimitation, and illustrates the importance to conserve polyploid populations as potential sources of diversification due to evolutionary significance of genome duplications in plant evolution.

Degradation of sexual reproduction in Veronica filiformis after introduction to Europe.

BACKGROUND: Baker's law predicts that self-incompatible plant species are generally poor colonizers because their mating system requires a high diversity of genetically differentiated individuals and thus self-compatibility should develop after long-distance dispersal. However, cases like the introduction of the self-incompatible Veronica filiformis (Plantaginaceae) to Europe constitute an often overlooked alternative to this rule. This species was introduced from subalpine areas of the Pontic-Caucasian Mountains and colonized many parts of Central and Western Europe in the last century, apparently without producing seeds. To investigate the consequences of the absence of sexual reproduction in this obligate outcrosser since its introduction, AFLP fingerprints, flower morphology, pollen and ovule production and seed vitality were studied in introduced and native populations. RESULTS: Interpopulation crossings of 19 introduced German populations performed in the greenhouse demonstrated that introduced populations are often unable to reproduce sexually. These results were similar to intrapopulation crossings, but this depended on the populations used for crossings. Results from AFLP fingerprinting confirmed a lack of genetic diversity in the area of introduction, which is best explained by the dispersal of clones. Flower morphology revealed the frequent presence of mutations affecting the androecium of the flower and decreasing pollen production in introduced populations. The seeds produced in our experiments were smaller, had a lower germination rate and had lower viability than seeds from the native area. CONCLUSIONS: Taken together, our results demonstrate that V. filiformis was able to spread by vegetative means in the absence of sexual reproduction. This came at the cost of an accumulation of phenotypically observable mutations in reproductive characters, i.e. Muller's ratchet.

Safer DNA extraction from plant tissues using sucrose buffer and glass fiber filter.

For some plant species, DNA extraction and downstream experiments are inhibited by various chemicals such as polysaccharides and polyphenols. This short communication proposed an organic-solvent free (except for ethanol) extraction method. This method consists of an initial washing step with STE buffer (0.25 M sucrose, 0.03 M Tris, 0.05 M EDTA), followed by DNA extraction using a piece of glass fiber filter. The advantages of this method are its safety and low cost. The purity of the DNA solution obtained using this method is not necessarily as high as that obtained using the STE/CTAB method, but it is sufficient for PCR experiments. These points were demonstrated empirically with two species, Japanese speedwell and common dandelion, for which DNA has proven difficult to amplify via PCR in past studies.

Species delimitation using dominant and codominant multilocus markers.

We propose a method for delimiting species based on dominant or codominant multilocus data using Gaussian clustering with a noise component for outliers. Case studies show that provisional species delimited using Gaussian clustering based on dominant multilocus data correspond well with provisional species delimited based on other data. However, the performance of Gaussian clustering in delimiting species based on few codominant markers was only moderate. Species represented by few individuals are usually included in the noise component because clusters are difficult to recognize with limited data. As alternative methods, we evaluated two model-based clustering methods originally proposed to infer population structure and assign individuals to populations based on the assumption of Hardy-Weinberg equilibrium within populations, namely STRUCTURE and STRUCTURAMA, as well as the "fields for recombination" approach. The latter resulted in lumping all individuals of each data set with codominant markers together, and whereas STRUCTURE often provides no decision about the number of clusters, STRUCTURAMA usually yields correct or almost correct numbers of clusters. The classification success of STRUCTURAMA analyses based on codominant markers was very good, but its performance with dominant markers was less consistent. Based on the classification success of the different methods for delimiting species with dominant and codominant multilocus markers in the case studies, we recommend using Gaussian clustering for data sets with dominant markers and STRUCTURAMA for data sets with codominant markers.

Disentangling phylogeography, polyploid evolution and taxonomy of a woodland herb (Veronica chamaedrys group, Plantaginaceae s.l.) in southeastern Europe.

Southeastern Europe is a centre of European biodiversity, but very little is known about factors causing the observed richness. Here, we contribute to fill this gap by reconstructing the spatio-temporal diversification of the cytologically variable and taxonomically intricate complex of Veronica chamaedrys (Plantaginaceae s.l.), growing in open forests, forest edges and grasslands, with flow cytometry, molecular markers (AFLPs, plastid DNA sequences) and morphometry. Our results show that both diploid and tetraploid cytotypes are widespread, but diploids predominate on the southern Balkan Peninsula. Plastid sequences suggest a first split into three main lineages in the mid-Pleistocene and a continuous diversification during the last 0.4 my. Two of the identified plastid lineages coincide with geographically distinct AFLP clusters. Altogether, the genetic data suggest forest refugia on the southern-most Balkan Peninsula (Greece), in Bulgaria, Istria (Croatia and Slovenia) and maybe the southeastern Carpathians (Romania). Morphometric and genetic data show little congruence with current taxonomy.

Phylogeny of Veronica in the Southern and Northern Hemispheres based on plastid, nuclear ribosomal and nuclear low-copy DNA.

The cosmopolitan and ecologically diverse genus Veronica with approximately 450 species is the largest genus of the newly circumscribed Plantaginaceae. Previous analyses of Veronica DNA sequences were in stark contrast to traditional systematics. However, analyses did not allow many inferences regarding the relationship between major groups identified, hindering further analysis of diversification and evolutionary trends in the genus. To resolve the backbone relationships of Veronica, we added sequences from additional plastid DNA regions to existing data and analyzed matching data sets for 78 taxa and more than 5000 aligned characters from nuclear ribosomal DNA and plastid DNA regions. The results provide the best resolved and supported estimate of relationships among major groups in the Northern (Veronica s. str.) and Southern Hemisphere (hebes). We present new informal names for the five main species groups within the Southern Hemisphere sect. Hebe. Furthermore, in two instances we provide morphological and karyological characters supporting these relationships. Finally, we present the first evidence from nuclear low-copy CYCLOIDEA2-region to compare results from the plastid genome with the nuclear genome.

Evolutionary rates in Veronica L. (Plantaginaceae): disentangling the influence of life history and breeding system.

The evolutionary rate at which DNA sequences evolve is known to differ between different groups of organisms. However, the reasons for these different rates are seldom known. Among plants, the generation-time hypothesis, which states that organisms that reproduce faster also have more DNA substitutions per time, has gained most popularity. We evaluate the generation-time hypothesis using 131 DNA sequences from the plastid trnLF region and the nuclear ribosomal ITS region of the genus Veronica (Plantaginaceae). We also examine the alternative hypothesis that a higher substitution rate is correlated with selfing breeding system. Selfing is associated with annual life history in many organisms and may thus often be the underlying reason for observed correlations of annual life history with other characters. We provide evidence that annual life history is more likely to be the responsible factor for higher substitution rates in Veronica than a selfing breeding system. Nevertheless, the way in which annual life history may influence substitution rate in detail remains unknown, and some possibilities are discussed.

Phylogeographic patterns in the Australasian genus Chionohebe (Veronica s.l., Plantaginaceae) based on AFLP and chloroplast DNA sequences.

The alpine genus Chionohebe is one of seven genera in the southern hemisphere Hebe complex. The main aims of this study were to infer the evolutionary relationships and assess phylogeographic patterns among the six species of Chionohebe, determine the origin of the two species with trans-Tasman distributions, and test species delimitations and specimen identifications based on morphology. Analyses of AFLP data recovered five major lineages within Chionohebe, some of which corresponded to species and varieties as currently circumscribed. Although the cushion chionohebes were strongly supported as monophyletic, the sole non-cushion species, C. densifolia, was sister to Parahebe trifida, and thus the AFLP data do not support a monophyletic Chionohebe as usually circumscribed. Strong north/south and west/east phylogeographic patterns were found among and within the main AFLP lineages in New Zealand. Analyses of chloroplast DNA (cpDNA) revealed eight haplotypes in Chionohebe, but these did not correspond to current taxonomy or geography due to widespread interspecific haplotype sharing. Based on both AFLP and cpDNA results, the two trans-Tasman species are shown to have originated in New Zealand and dispersed to Australia independently.

Comparative phylogeography of the Veronica alpina complex in Europe and North America.

The Veronica alpina complex comprises eight species of alpine habitats over a wide range of mountain systems in the Northern Hemisphere. The occurrence of sympatric species in the European and North American mountain systems allowed us not only to investigate the effect of the ice ages on intraspecific phylogeographical patterns and genetic diversity in different continents of the Northern Hemisphere, but also to compare these patterns in closely related species. Plastid DNA trnL-F sequences and AFLP (amplified fragment length polymorphism) fingerprints were used to infer the phylogenetic history of the group and phylogeographical patterns within species. Hybrid origin of tetraploid eastern North American V. wormskjoldii from western North American V. nutans (= V. wormskjoldii s.l.) and Eurasian V. alpina is suggested. A number of phylogeographical groups have been found both in V. alpina from Europe and in V. nutans from western North America. Phylogeographical substructuring in the Alps is inferred for V. alpina but not for V. bellidioides, which is moreover characterized by an overall very low genetic diversity. Western North American V. cusickii is much more genetically diverse than its sympatric relative, V. nutans, an effect that is likely due to differences in the breeding system. Populations of V. nutans are differentiated into three groups, those from the Cascades and from the southern and the northern Rocky Mountains. Genetic diversity seems to be higher in the North American V. nutans than in the morphologically and ecologically similar European V. alpina. A possible scenario to explain this pattern is suggested.

Phylogeny of Veronica--a combination of molecular and chemical evidence.

A molecular phylogenetic tree of the genus Veronica, based on sequences of the ITS region of nuclear ribosomal DNA, provided a framework to test the use of iridoid chemical data in the classification of the genus. The distribution of ten iridoid glucosides, determined in 27 species of Veronica, supported the molecular results. The nine clades of the molecular tree appeared to correlate with a specific pattern of catalpol esters and chromosome base numbers. Some relationships within Veronica, previously suggested on the basis of morphological and karyological characters, were supported by our data. Neither the molecular, nor the chemical features were in favour of the classical circumscription and classification of the genus.

Genome size variation and evolution in Veronica.

BACKGROUND AND AIMS: The amount of DNA per chromosome set is known to be a fairly constant characteristic of a species. Its interspecific variation is enormous, but the biological significance of this variation is little understood. Some of the characters believed to be correlated with DNA amount are alpine habitat, life history and breeding system. In the present study, the aim is to distinguish between direct causal connections and chance correlation of the amount of DNA in the genus Veronica. METHODS: Estimates of DNA amount were analysed for 42 members of Veroniceae in connection with results from a phylogenetic analysis of plastid trnL-F DNA sequences and tested correlations using standard statistical tests, phylogenetically independent contrasts and a model-based generalized least squares method to distinguish the phylogenetic effect on the results. KEY RESULTS: There appears to be a lower upper limit for DNA amount in annuals than in perennials. Most DNAC-values in Veroniceae are below the mean DNA C-value for annuals in angiosperms as a whole. However, the long-debated correlation of low genome size with annual life history is not significant (P = 0.12) using either standard statistical tests or independent contrasts, but it is significant with the generalized least squares method (P < 0.01). CONCLUSIONS: The correlation of annual life history and low genome size found in earlier studies could be due to the association of annual life history and selfing, which is significantly correlated with low genome size using any of the three tests applied. This correlation can be explained by models showing a reduction in transposable elements in selfers. A significant correlation of higher genome sizes with alpine habitats was also detected.

Iridoid glucosides of Paederota lutea and the relationships between Paederota and Veronica.

In a chemical investigation of the water soluble compounds in Paederota lutea eight known iridoids were isolated together with a new one with a 8,9-double bond, namely paederotoside (10-O-benzoyl-6'-O-alpha-arabino(1-->6)-beta-glucopyranosyl arborescosidic acid) and the 6-hydroxy-flavone glucoside 4'-O-methylscutellarein 7-O-beta-glucopyranoside. The known iridoid glucosides were 8-epiloganic acid, gardoside, aucubin, catalpol and the 6-O-esters of catalpol: veronicoside, catalposide, amphicoside and verproside. The compounds isolated show that Paederota has a glycoside content almost identical to that of a general Veronica species, and this is in good agreement with the results from recent investigations of the phylogeny of Veronica and its closest relatives, where Paederota is placed as a sister-group next to Veronica. In an analysis of the iridoid glucosides from some of these relatives, it is shown that Veronica, Paederota, Picrorhiza and Veronicastrum are all characterized by containing 6-O-esters of catalpol. Some less closely related taxa namely: Lagotis, Wulfenia, Plantago, Aragoa and Globularia instead contain 10-O-esters of catalpol or aucubin.