How many taxa? Spatiotemporal evolution and taxonomy of Amphoricarpos (Asteraceae, Carduoideae) on the Balkan Peninsula

Phytochemicals and Their Correlation with Molecular Data in Micromeria and Clinopodium (Lamiaceae) Taxa

A study of the phytochemical and molecular characteristics of ten Micromeria and six Clinopodium taxa (family Lamiaceae) distributed in the Balkan Peninsula was carried out. The phytochemicals detected in essential oils by gas chromatography, mass spectrometry, and molecular data amplified fragment length polymorphism were used to study the taxonomic relationships among the taxa and the correlations between phytochemical and molecular data. STRUCTURE analysis revealed three genetic groups, while Bayesian Analysis of Population Structure grouped the studied taxa into 11 clusters nested in the groups obtained by STRUCTURE. Principal components analysis performed with the 21 most represented compounds in the essential oils yielded results that were partly consistent with those obtained by STRUCTURE and neighbour-joining analyses. However, their geographic distributions did not support the genetic grouping of the studied taxa and populations. The Mantel test showed a significant correlation between the phytochemical and genetic data (r = 0.421, p < 0.001). Genetic distance explained 17.8% of the phytochemical distance between populations. The current taxonomic position of several of the studied taxa is yet to be satisfactorily resolved, and further studies are needed. Such future research should include nuclear and plastid DNA sequences from a larger sample of populations and individuals.

Genome Size Variation in Dianthus sylvestris Wulfen sensu lato (Caryophyllaceae)

Genome size (GS) is an important characteristic that may be helpful in delimitation of taxa, and multiple studies have shown correlations between intraspecific GS variation and morphological or environmental factors, as well as its geographical segregation. We estimated a relative GS (RGS) of 707 individuals from 162 populations of Dianthus sylvestris with a geographic focus on the Balkan Peninsula, but also including several populations from the European Alps. Dianthus sylvestris is morphologically variable species thriving in various habitats and six subspecies have been recognized from the Balkan Peninsula. Our RGS data backed-up with chromosome counts revealed that the majority of populations were diploid (2n = 30), but ten tetraploid populations have been recorded in D. sylvestris subsp. sylvestris from Istria (Croatia, Italy). Their monoploid RGS is significantly lower than that of the diploids, indicating genome downsizing. In addition, the tetraploids significantly differ from their diploid counterparts in an array of morphological and environmental characteristics. Within the diploid populations, the RGS is geographically and only partly taxonomically correlated, with the highest RGS inferred in the southern Balkan Peninsula and the Alps. We demonstrate greater RGS variation among the Balkan populations compared to the Alps, which is likely a result of more pronounced evolutionary differentiation within the Balkan Peninsula. In addition, a deep RGS divergence within the Alps likely points to persistence of the alpine populations in different Pleistocene refugia.

Allele Sorting as a Novel Approach to Resolving the Origin of Allotetraploids Using Hyb-Seq Data: A Case Study of the Balkan Mountain Endemic Cardamine barbaraeoides

Mountains of the Balkan Peninsula are significant biodiversity hotspots with great species richness and a large proportion of narrow endemics. Processes that have driven the evolution of the rich Balkan mountain flora, however, are still insufficiently explored and understood. Here we focus on a group of Cardamine (Brassicaceae) perennials growing in wet, mainly mountainous habitats. It comprises several Mediterranean endemics, including those restricted to the Balkan Peninsula. We used target enrichment with genome skimming (Hyb-Seq) to infer their phylogenetic relationships, and, along with genomic in situ hybridization (GISH), to resolve the origin of tetraploid Cardamine barbaraeoides endemic to the Southern Pindos Mts. (Greece). We also explored the challenges of phylogenomic analyses of polyploid species and developed a new approach of allele sorting into homeologs that allows identifying subgenomes inherited from different progenitors. We obtained a robust phylogenetic reconstruction for diploids based on 1,168 low-copy nuclear genes, which suggested both allopatric and ecological speciation events. In addition, cases of plastid-nuclear discordance, in agreement with divergent nuclear ribosomal DNA (nrDNA) copy variants in some species, indicated traces of interspecific gene flow. Our results also support biogeographic links between the Balkan and Anatolian-Caucasus regions and illustrate the contribution of the latter region to high Balkan biodiversity. An allopolyploid origin was inferred for C. barbaraeoides, which highlights the role of mountains in the Balkan Peninsula both as refugia and melting pots favoring species contacts and polyploid evolution in response to Pleistocene climate-induced range dynamics. Overall, our study demonstrates the importance of a thorough phylogenomic approach when studying the evolution of recently diverged species complexes affected by reticulation events at both diploid and polyploid levels. We emphasize the significance of retrieving allelic and homeologous variation from nuclear genes, as well as multiple nrDNA copy variants from genome skim data.

Genetic Diversity and Demographic History of Wild and Cultivated/Naturalised Plant Populations: Evidence from Dalmatian Sage (Salvia officinalis L., Lamiaceae)

Dalmatian sage (Salvia officinalis L., Lamiaceae) is a well-known aromatic and medicinal Mediterranean plant that is native in coastal regions of the western Balkan and southern Apennine Peninsulas and is commonly cultivated worldwide. It is widely used in the food, pharmaceutical and cosmetic industries. Knowledge of its genetic diversity and spatiotemporal patterns is important for plant breeding programmes and conservation. We used eight microsatellite markers to investigate evolutionary history of indigenous populations as well as genetic diversity and structure within and among indigenous and cultivated/naturalised populations distributed across the Balkan Peninsula. The results showed a clear separation between the indigenous and cultivated/naturalised groups, with the cultivated material originating from one restricted geographical area. Most of the genetic diversity in both groups was attributable to differences among individuals within populations, although spatial genetic analysis of indigenous populations indicated the existence of isolation by distance. Geographical structuring of indigenous populations was found using clustering analysis, with three sub-clusters of indigenous populations. The highest level of gene diversity and the greatest number of private alleles were found in the central part of the eastern Adriatic coast, while decreases in gene diversity and number of private alleles were evident towards the northwestern Adriatic coast and southern and eastern regions of the Balkan Peninsula. The results of Ecological Niche Modelling during Last Glacial Maximum and Approximate Bayesian Computation suggested two plausible evolutionary trajectories: 1) the species survived in the glacial refugium in southern Adriatic coastal region with subsequent colonization events towards northern, eastern and southern Balkan Peninsula; 2) species survived in several refugia exhibiting concurrent divergence into three genetic groups. The insight into genetic diversity and structure also provide the baseline data for conservation of S. officinalis genetic resources valuable for future breeding programmes.

Pleistocene phylogeography and cryptic diversity of a tiger beetle, Calomera littoralis, in North-Eastern Mediterranean and Pontic regions inferred from mitochondrial COI gene sequences

Background. Calomera littoralis is a Palearctic species, widely distributed in Europe; inhabiting predominantly its Atlantic, Mediterranean and Black Sea coastlines. Methods. Its phylogeography on the Balkan Peninsula and on the north-western Black Sea coast was inferred using a 697 bp long portion of the mitochondrial COI gene, amplified from 169 individuals collected on 43 localities. Results. The results revealed two genetically divergent groups/lineages, the southern one inhabiting both the Balkan Peninsula and the Pontic Region and the northern one found exclusively in the Pontic Region. Species delimitation based on DNA barcoding gap suggested an interspecific level of divergence between these groups. Multivariate analysis of eight male and female morphometric traits detected no difference between the groups, implying they may represent cryptic species. The Bayesian time-calibrated reconstruction of phylogeny suggested that the lineages diverged ca. 2.3 Ma, in early Pleistocene. Discussion. The presence of the two genetically divergent groups results most likely from contemporary isolation of the Pontic basin from the Mediterranean that broke the continuous strip of coastal habitats inhabited by C. littoralis. Demographic analyses indicated that both lineages have been in demographic and spatial expansion since ca. 0.15 Ma. It coincides with the terminal stage of MIS-6, i.e., Wartanian/Saalian glaciation, and beginning of MIS-5e, i.e., Eemian interglacial, during which, due to eustatic sea level rise, a wide connection between Mediterranean and the Pontic basin was re-established. This, along with re-appearance of coastal habitats could initiate north-east expansion of the southern lineage and its secondary contact with the northern one. The isolation of the Pontic basin from the Mediterranean during the Weichselian glaciation most likely did not have any effect on their phylogeography.

Heteroploid Knautia drymeia includes K. gussonei and cannot be separated into diagnosable subspecies

PREMISE OF THE STUDY

Knautia drymeia is a morphologically variable, diploid and tetraploid temperate forest understory species distributed in southeastern Europe and adjacent areas. The species is an excellent system to explore the influence of polypoidy on taxonomic delineations, the role of hybridization among genetically distant populations in polyploid evolution, and the impact of glacial refugia on the evolution of polyploids.

METHODS

Amplified fragment length polymorphism fingerprinting and multivariate analyses of morphological characters were performed on 57 populations spanning the distribution area of K. drymeia. K-means clustering, comparison of in-silico tetraploids and observed tetraploids, and a phylogeographic analysis using relaxed random walks were used to explore the genetic structure within the diploids, to infer the origin of the tetraploids and to reconstruct range expansion through time. Further, we contrasted the morphology and genetic groups with current taxonomy and evaluated the status of the tetraploid Apennine endemic K. gussonei and the intraspecific taxa of K. drymeia.

KEY RESULTS

The genetic structure was strongly geographically correlated and yielded four genetic groups; K. gussonei was inseparable from K. drymeia. Distributions of diploid lineages are suggestive of glacial refugia in the northwesternmost and southeastern Balkan Peninsula. Polyploids originated at least two times, as autopolyploids and probably additionally also as allopolyploids. Morphological divergence corresponded with neither genetic groups nor current taxonomy.

CONCLUSIONS

Genetic and morphometric data confirmed neither divergence of K. gussonei nor recognition of subspecies within K. drymeia. We therefore propose treating K. drymeia as a morphologically and genetically variable species without infraspecific taxa.