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

Spatial Distribution of Genetic, Ploidy, and Morphological Variation of the Edaphic Steno-Endemic Alyssum moellendorfianum (Brassicaceae) from the Western Balkans

Polyploidy is a powerful mechanism driving genetic, physiological, and phenotypic changes among cytotypes of the same species across both large and small geographic scales. These changes can significantly shape population structure and increase the evolutionary and adaptation potential of cytotypes. Alyssum moellendorfianum, an edaphic steno-endemic species with a narrow distribution in the Balkan Peninsula, serves as an intriguing case study. We conducted a comprehensive analysis of genetic diversity and population structure across the species' range, employing an array of genetic techniques (nuclear microsatellites, amplified fragment length polymorphisms, and plastid DNA sequences), flow cytometry (FCM), morphometry, and pollen analysis. The study reveals two genetic lineages: spatially distributed diploid and tetraploid cytotypes. Clear divergence between diploids and tetraploids was shown by AFLP, while plastid DNA sequences confirmed private haplotypes in each of the studied populations. Higher genetic diversity and allelic richness following the north-south pattern were documented in tetraploids compared to diploids, as indicated by nuclear microsatellites. Morphometric analysis via principal component analysis (PCA) and canonical discriminant analysis (CDA) did not reveal any divergence between diploid and tetraploid cytotypes. Nonetheless, a distinction in pollen size was clearly observed. The results suggest an autopolyploid origin of tetraploids from diploid ancestors. Despite the population fragmentation in a very small geographic range, these populations harbour high genetic diversity, which would allow them to remain stable if natural processes remain undisturbed.

Patterns of Genetic and Morphological Variability of Teucrium montanum sensu lato (Lamiaceae) on the Balkan Peninsula

The Balkan Peninsula represents an important center of plant diversity, exhibiting remarkable ecological heterogeneity that renders it an optimal region for studying the diversification patterns of complex taxa such as Teucrium montanum. In the Balkan Peninsula, T. montanum is a highly plastic and morphologically variable species with unresolved taxonomic status. To ascertain the patterns of genetic and morphological diversification, a comparative genetic and morphological analysis was conducted. In total, 57 populations were subjected to analysis using AFLP and a multivariate morphometric approach. A Bayesian analysis of population structure distinguished two main genetic clusters, labelled A and B. Cluster B was found to be geographically restricted to the northwestern Dinarides, while cluster A occurred in the rest of the Balkans. Genetic cluster A was further subdivided into four subclusters that were spatially separated from each other. The contact populations between the subclusters exhibited a mixed genetic structure. There was a partial correlation between genetic and morphological diversification. The peripheral populations of the genetic clusters displayed morphological differences, while both genetic and morphological differences decreased in the contact zones. The observed genetic structure can be attributed to the reproductive biology of this species and the complex geological history of the Balkan Peninsula.

Patterns of Genetic Variation of Nepeta nuda L. from the Central Balkans: Understanding Drivers of Chemical Diversity

Nepeta nuda L., a notable medicinal species in the tradition of the Balkan region, is a rich source of bioactive iridoids and phenolics previously described as high-resolution taxonomical classifiers for the genus Nepeta. However, their potential in investigating intra-species differentiation is here described for the first time. The aim was to recognize the sources of natural chemical diversity and their association with the genetic variability both within and among N. nuda populations in the Central Balkans. Chemical diversity was assessed from methanol extracts and essential oils through untargeted and targeted metabolomics using state-of-the-art analytical tools, covering a broad spectrum of compounds that represent the N. nuda metabolome. We found that chemodiversity primarily resides within populations of N. nuda, and similar results were obtained at the DNA level using microsatellite markers. The low genetic and chemical differentiation of the studied N. nuda populations implies that their metabolomic profiles may be less influenced by geographic distance and variable environmental conditions within the Central Balkans, as they are under the pivotal control of their genetic backgrounds. Screening the distribution of the major bioactive compounds belonging to phenolics (phenolic acids and flavonoids) and iridoids (both aglycones and glycosylated forms), within and among N. nuda populations, is able to guarantee mass spectrometry-based tools for the selection of elite representative genotypes with practical importance. The knowledge acquired will allow us to delve deeper into the molecular background of N. nuda chemical diversity, which is the course of our further work.

Phylogeographic structure of common sage (Salvia officinalis L.) reveals microrefugia throughout the Balkans and colonizations of the Apennines

Studying the population-genetic and phylogeographic structures of a representative species of a particular geographical region can not only provide us with information regarding its evolutionary history, but also improve our understanding of the evolutionary processes underlying the patterns of species diversity in that area. By analysing eight highly polymorphic microsatellite loci and two chloroplast DNA regions, we have investigated the influence of Pleistocene climate fluctuations on the evolutionary history of Salvia officinalis L. (common sage). The populations with the highest genetic diversity were located in the central parts of the Balkan distribution range. A large group of closely related haplotypes was distributed throughout the Balkans and the central Apennines, while the private lineage occupied the southern Apennines. In addition, two highly differentiated lineages were scattered only over the Balkans. The results suggest that a single refugium of the studied species from the last glacial period was located in the central part of the range in the Balkans. Numerous microrefugia, probably spanning several glaciation cycles, were scattered across the Balkans, while colonisation of the Apennines from the Balkans occurred at least on two occasions.

The Renaissance of Wild Food Plants: Insights from Tuscany (Italy)

This paper provides an overview of wild food plants traditionally used in the gastronomy of Tuscany, an Italian region with high biological diversity and whose cultural heritage is well known. Forty-nine bibliographic sources, including five unpublished studies, were reviewed. A list of species with ecological characteristics, plant parts used, use category (food, liquor, or seasoning), methods of preparation (raw or cooked), and recipes is presented. The use of 357 taxa (3711 use reports, URs), was recorded, belonging to 215 genera and 72 botanical families. Over the total taxa, 12 are new for Tuscany, 52 seem not to be present in other Italian regions, and 54 were not detected in the consulted European ethnobotanical literature. Of these taxa, 324 (3117 URs) were used as food, while 49 (178 URs) and 81 (416 URs) were used for liquor and seasoning, respectively. Of the 17 different food recipes, cooked vegetables constituted the largest group, followed by salads, omelets, snacks, and fillings. The chemical composition of the recorded food plants and the possible safety risks associated to their consumption, as well as their traditional medicinal use, are also shown. This review highlights the richness of ethnobotanical knowledge in Tuscany. Such biocultural heritage can be a “source of inspiration” for agriculture. As a reservoir of potential new crops, wild edible flora may contribute to the development of emerging horticultural sectors such as vertical farming and microgreens production. Moreover, the nutrient content and healthy properties of many wild food plants reported in this study has the ability to meet consumer demand for functional foods.

Population structure and adaptive variation of Helichrysum italicum (Roth) G. Don along eastern Adriatic temperature and precipitation gradient

Immortelle (Helichrysum italicum (Roth) G. Don; Asteraceae) is a perennial plant species native to the Mediterranean region, known for many properties with wide application mainly in perfume and cosmetic industry. A total of 18 wild H. italicum populations systematically sampled along the eastern Adriatic environmental gradient were studied using AFLP markers to determine genetic diversity and structure and to identify loci potentially responsible for adaptive divergence. Results showed higher levels of intrapopulation diversity than interpopulation diversity. Genetic differentiation among populations was significant but low, indicating extensive gene flow between populations. Bayesian analysis of population structure revealed the existence of two genetic clusters. Combining the results of F_ST - outlier analysis (Mcheza and BayeScan) and genome-environment association analysis (Samβada, LFMM) four AFLP loci strongly associated with the bioclimatic variables Bio03 Isothermality, Bio08 Mean temperature of the wettest quarter, Bio15 Precipitation seasonality, and Bio17 Precipitation of driest quarter were found to be the main variables driving potential adaptive genetic variation in H. italicum along the eastern Adriatic environmental gradient. Redundancy analysis revealed that the partitioning of genetic variation was mainly associated with the adaptation to temperature oscillations. The results of the research may contribute to a clearer understanding of the importance of local adaptations for the genetic differentiation of Mediterranean plants and allow the planning of appropriate conservation strategies. However, considering that the identified outlier loci may be linked to genes under selection rather than being the target of natural selection, future studies must aim at their additional analysis.

Exploration of genetic, morphological and essential oil variation reveals tools for the authentication and breeding of Salvia pomifera subsp. calycina (Sm.) Hayek

Salvia pomifera subsp. calycina (Sm.) Hayek (Lamiaceae), is an Eastern Mediterranean element, which is used in traditional medicine and cuisine in the same manner as S. fruticosa Mill. and S. officinalis L.. The essential oil (EO) and the extracts of S. pomifera possess bioactive compounds with anti-proliferative, anticholinesterase, antioxidant, antiviral and antifungal properties. In this study, the chemical (EO), genetic (DNA microsatellites, SSRs) and morphological diversity of forty-nine individuals of Salvia pomifera subsp. calycina, originating from five natural populations of the Peloponnese (Greece) were determined, in order to explore the potential for successful breeding and to reveal tools and biomarkers for identification and authentication. Chemical and genetic analyses revealed high levels of variation both within and among populations, while morphological analysis mainly within populations. Essential oil yield ranged from 1.79 to 5.79 ml 100 g^-1 dry wt, among individuals while β-thujone ranged from 6.04 to 64.75%. Consistency was found in the EO yield and composition of specific individuals, when sampled during the same period, for three consecutive years, while the analysis during spring and summer months showed differentiation that still retained individual's discrimination. Genetic analysis using SSRs showed that the observed population heterozygosity (H_o) ranged from 0.48 to 0.67, while high number of private alleles were revealed in all populations. Considerable genetic differentiation was observed among the three Salvia taxa (S. pomifera subsp. calycina, S. fruticosa, S. officinalis) (Fst values ranged from 0.27 to 0.48) and lower among S. pomifera subsp. calycina populations (Fst values ranged from 0.06 to 0.13). The great variation that was revealed in all measured traits, in combination with the demonstrated, genetically based, consistency of their EO yield and composition, advocates to a successful breeding, whereas SSR genotyping presents a strong identification and authentication tool.

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.

Morphological, genetic and epigenetic aspects of homoploid hybridization between Salvia officinalis L. and Salvia fruticosa Mill

The inheritance of phenotypic, genetic and epigenetic traits in hybridization events is difficult to predict, as numerous evolutionary, ecological, and genetic factors can play a crucial role in the process of hybridization. In the middle Adriatic island of Vis, we investigated hybridization between Salvia officinalis and S. fruticosa at morphological, genetic and epigenetic levels. SSR results revealed that hybrid individuals were characterized by diploid set of chromosomes suggesting homoploid hybridization. A well-defined group that mostly comprised of F1 generation individuals was detected. For the majority of analysed morphological characteristics, hybrids were placed in-between parental taxa, while at the same time, values of different genetic parameters were mostly higher in hybrids than in parental species. The results revealed a high contrast in the levels of phenotypic variability and epigenetic excitation between parental taxa. Environmental niche modelling confirmed that in the studied location S. officinalis experiences optimal climatological conditions, while S. fruticosa struggles with unsuitable conditions. Very low levels of gene flow between the parental species were detected. In addition, contrasting levels of epigenetic excitation in the studied groups clearly demonstrated the importance of an epigenetic response to an altered environment and confirmed the trans-generational nature of the epigenetic changes.

Understanding biogeographical patterns in the western Balkan Peninsula using environmental niche modelling and geostatistics in polymorphic Edraianthus tenuifolius

The Balkan Peninsula represents one of the three southern European glacial refugia where biodiversity persisted throughout the climatically unstable Quaternary. This study considered the 'refugia within refugia' model, which assumes the environmental heterogeneity over time and space in larger refugia. To better understand patterns shaped during the Quaternary climatic oscillations, the present and last glacial maximum (LGM) environmental conditions and current morphological variability of Edraianthus tenuifolius, an endemic plant of the western Balkans with a well-known therphical structure, were considered. Potential present and LGM distributions were studied through environmental niche modelling using 161 data of occurrences and six bioclimatic variables, hindcasting the model to LGM conditions using three different global circulation models. To test the geographical variability of the reproductive region, 41 characters of 667 inflorescences from 35 populations within the distribution range were measured. Geographical patterns, using geostatistics together with univariate and multivariate statistical approaches, were analysed. The environmental niche model suggested the current potential distribution in correspondence to its known occurrences. The hindcast to LGM conditions suggested two separate areas of environmental suitability, one in the present-day northern Adriatic coast of Croatia (Istrian Peninsula, Kvarner) and another on the present-day south-eastern Adriatic coast (southern Dalmatia, Montenegro and northern Albania). Morphological variability showed a similar pattern, where southern populations separated from northern populations according to a major split in the central part of its distribution range (central Dalmatia). On other levels, stronger barriers were predicted to separate northern populations from the eastern Istrian Peninsula and the Kvarner area. The results suggested congruent biogeographical patterns to the already known phylogeographical structure. Both environmental niche modelling and the geographical variability of morphological characters suggested spatial partitioning, indicating the potential presence of two separate refugia during the LGM.

Genetic diversity and metabolic profile of Salvia officinalis populations: implications for advanced breeding strategies

As a result of this work, we were able to characterize seven indigenous to Greece Salvia officinalis populations using genetic and metabolomic tools. These tools can be used to select the most promising genotypes, capable to design future breeding programs for high valuable varieties. An initial investigation was carried out to compare the genetic and metabolic diversity in S. officinalis grown in Greece and to discern the relationship between the two sets of data. Analysis of inter-simple sequence repeats (ISSR) revealed significant genetic differences among seven sage populations, which were grouped into three main clusters according to an UPGMA ISSR data-based dendrogram and Principle Coordinate Analysis. 80 loci were scored of which up to 90% were polymorphic at species level. According to the composition of their essential oil, the populations were classified into two chemotypes: 1.8 cineole/α-thujone and α-thujone/1.8 cineole. Additionally, a targeted ultra performance liquid chromatography (UPLC-MS/MS) method was used to qualify and quantify phenolic compounds in methanolic extracts of the seven sage genotypes according to which they were districted in six clusters among the sage populations. The main compounds characterizing the seven genotypes were rosmarinic acid and carnosol, followed by apigenin-7-O-glucoside (Ap7glc), and luteolin-7-O-glucoside (Lu7glc). The correlation between matrices obtained from ISSR data and metabolic profiles was non-significant. However, based on the differences in metabolic fingerprint, we aimed to define populations using as main selection criteria the high polyphenol content and desired essential oil composition, using state to the art analytical tools for the identification of parent lines for breeding programs.