Asymmetric Introgression and Cryptic Natural Hybridization between Two Species of Teucrium Section Polium (Lamiaceae) on the Balkan Peninsula.

In this work, we analyzed the morphology and genetic structure of Teucrium montanum, T. capitatum and their hybrid T. × rohlenae from three syntopic populations. A morphometric study showed that the parents and their hybrids exhibited continuous morphological variation, with the hybrid positioned exactly between the parents. Genetic analysis revealed that plants morphologically identified as T. × rohlenae are fertile hybrids that produce hybrid swarms dominated by later-generation hybrids. This suggests that introgression, rather than speciation, is the more likely outcome of hybridization between these plant species. The extent and direction of gene flow between the two species differed markedly between the three syntopic localities. At the Trilj locality, it was clearly unidirectional, with T. capitatum playing the dominant role. At the Sicevo locality, gene flow was slightly asymmetric, favoring the genetic background of T. capitatum, while at the Sliven site, it was completely asymmetric in the opposite direction. The extreme case of unidirectional gene flow was observed at the Trilj locality where plants morphologically identified as T. montanum could not be genetically distinguished from T. capitatum. This suggests that interspecific hybridization occurred long ago, leading to introgression and cryptic hybrids, blurring of species boundaries and generating evolutionary noise.

Dig-up Primers: A Pipeline for Identification of Polymorphic Microsatellites Loci within Assemblies of Related Species.

Simple sequence repeats (SSRs) have become one of the most popular molecular markers and are used in numerous fields, including conservation genetics, population genetic studies, and genetic mapping. Advances in next-generation sequencing technology and the growing amount of genomic data are driving the development of bioinformatics tools for SSR marker design. These tools work with different combinations of input data, which can be raw reads or assemblies, and with one or more input datasets. We present here a new strategy and implementation of a simple standalone pipeline that utilizes more than one assembly for the in silico design of PCR primers for microsatellite loci in more than one species. Primers are tested in silico to determine if they are polymorphic, eliminating the need to test time-consuming cross-species amplification in the laboratory. The end result is a set of markers that are in silico polymorphic in all analyzed species and have great potential for the identification of interspecies hybrids. The efficiency of the tool is demonstrated using two examples at different taxonomic levels and with different numbers of input assemblies to generate promising, high-quality SSR markers.

Revealing the Diversity and Complex Relationships of Croatian Olive Germplasm.

Olive growing in Croatia has a long tradition and is of great economic and social impact. The present study includes a set of 108 tree samples (88 samples corresponding to 60 presumed cultivars and 20 trees of unnamed ones) collected from 27 groves in the entire olive growing area, and is the most comprehensive survey to be conducted in Croatia. The genetic diversity, relationships, and structures of olive plants were studied using eight microsatellite loci. All loci were polymorphic and revealed a total of 90 alleles. A total of 74 different genotypes were detected that were subjected to further diversity and genetic relationship studies. The Fitch-Margoliash tree and Bayesian analysis of population structure revealed a complex relationship between the identified olive genotypes, which were clustered into three gene pools, indicating different origins of Croatian olive germplasms. Excluding the redundant germplasms, 44 different genotypes among the sampled trees of well-known cultivars and 16 new local germplasms were identified. In addition, we provide the etymology of 46 vernacular names, which confirms that the vast majority of traditional Croatian cultivars have common and widespread names. The results presented herein underline the importance of safeguarding local cultivars and conducting continuous surveys.

Lack of pollinators selects for increased selfing, restricted gene flow and resource allocation in the rare Mediterranean sage Salvia brachyodon.

Range contraction and habitat fragmentation can cause biodiversity loss by creating conditions that directly or indirectly affect the survival of plant populations. Fragmented habitats can alter pollinator guilds and impact their behavior, which may result in pollen/pollinator limitation and selection for increased selfing as a mechanism for reproductive assurance. We used Salvia brachyodon, a narrowly distributed and endangered sage from eastern Adriatic, to test the consequences of range contraction and habitat fragmentation. Molecular data indicate a severe and relatively recent species range reduction. While one population is reproductively almost completely isolated, moderate gene flow has been detected between the remaining two populations. The high pollen-to-ovule ratio and the results of controlled hand pollination indicate that S. brachyodon has a mixed mating system. Quantitative and qualitative differences in the community and behaviour of flower visitors resulted in limited pollination services in one population where no effective pollinator other than pollen and nectar robbers were observed. In this population, self-pollination predominated over cross-pollination. Various environmental factors, in which plant-pollinator interactions play a pivotal role, have likely created selection pressures that have led to genetic and phenotypic differentiation and different resource allocation strategies among populations.

Patterns of Leaf and Fruit Morphological Variation in Marginal Populations of Acer tataricum L. subsp. tataricum.

Marginal populations are usually smaller and more isolated and grow in less favourable conditions than those at the distribution centre. The variability of these populations is of high importance, as it can support the adaptations needed for the conditions that they grow in. In this research, the morphological variability of eight Tatar maple (Acer tataricum L. subsp. tataricum) populations was analysed. Tatar maple is an insect-pollinated and wind-dispersed shrub/tree, whose northwestern distribution edge is in southeastern Europe. Morphometric methods were used to analyse the variability of the populations using leaf and fruit morphology. The research revealed significant differences between and within populations. Furthermore, differences in the distribution of the total variability were noted, which suggest that different evolutionarily factors affect different plant traits. Correlation analysis confirmed a weak dependency between the vegetative and generative traits. In addition, no evidence was found for the presence of isolation by environment (IBE). However, the Mantel test for isolation by distance (IBD) was significant for the leaf morphometric traits and non-significant for the fruit morphometric traits. Being the marginal leading-edge populations, they are younger and were less likely to have had time for adaptation to local environments, which would have resulted in the development of IBE. Overall, edge populations of Tatar maple were characterised by great morphological variability, which helps these populations in their response to the intensive selective pressures they face in their environment.

Genetic diversity and selection signatures in a gene bank panel of maize inbred lines from Southeast Europe compared with two West European panels.

Southeast Europe (SEE) is a very important maize-growing region, comparable to the Corn belt region of the United States, with similar dent germplasm (dent by dent hybrids). Historically, this region has undergone several genetic material swaps, following the trends in the US, with one of the most significant swaps related to US aid programs after WWII. The imported accessions used to make double-cross hybrids were also mixed with previously adapted germplasm originating from several more distant OPVs, supporting the transition to single cross-breeding. Many of these materials were deposited at the Maize Gene Bank of the Maize Research Institute Zemun Polje (MRIZP) between the 1960s and 1980s. A part of this Gene Bank (572 inbreds) was genotyped with Affymetrix Axiom Maize Genotyping Array with 616,201 polymorphic variants. Data were merged with two other genotyping datasets with mostly European flint (TUM dataset) and dent (DROPS dataset) germplasm. The final pan-European dataset consisted of 974 inbreds and 460,243 markers. Admixture analysis showed seven ancestral populations representing European flint, B73/B14, Lancaster, B37, Wf9/Oh07, A374, and Iodent pools. Subpanel of inbreds with SEE origin showed a lack of Iodent germplasm, marking its historical context. Several signatures of selection were identified at chromosomes 1, 3, 6, 7, 8, 9, and 10. The regions under selection were mined for protein-coding genes and were used for gene ontology (GO) analysis, showing a highly significant overrepresentation of genes involved in response to stress. Our results suggest the accumulation of favorable allelic diversity, especially in the context of changing climate in the genetic resources of SEE.

Variation in Chloroplast Genome Size: Biological Phenomena and Technological Artifacts.

The development of bioinformatic solutions is guided by biological knowledge of the subject. In some cases, we use unambiguous biological models, while in others we rely on assumptions. A commonly used assumption for genomes is that related species have similar genome sequences. This is even more obvious in the case of chloroplast genomes due to their slow evolution. We investigated whether the lengths of complete chloroplast sequences are closely related to the taxonomic proximity of the species. The study was performed using all available RefSeq sequences from the asterid and rosid clades. In general, chloroplast length distributions are narrow at both the family and genus levels. In addition, clear biological explanations have already been reported for families and genera that exhibit particularly wide distributions. The main factors responsible for the length variations are parasitic life forms, IR loss, IR expansions and contractions, and polyphyly. However, the presence of outliers in the distribution at the genus level is a strong indication of possible inaccuracies in sequence assembly.

New insights in the Spanish gene pool of olive (Olea europaea L.) preserved ex situ and in situ based on high-throughput molecular markers.

In Spain, several local studies have highlighted the likely presence of unknown olive cultivars distinct from the approximately 260 ones previously described in the literature. Furthermore, recent advancements in identification techniques have significantly enhanced in terms of efficacy and precision. This scenario motivated a new nationwide prospecting effort aimed at recovering and characterizing new cultivated germplasm using high-throughput molecular markers. In the present study, the use of 96 EST-SNP markers allowed the identification of a considerable amount of new material (173 new genotypes) coming from areas with low intensification of production in different regions of Spain. As a result, the number of distinct national genotypes documented in the World Olive Germplasm Bank of IFAPA, Córdoba (WOGBC-ESP046) increased to 427. Likewise, 65 and 24 new synonymy and homonymy cases were identified, respectively. This rise in the number of different national cultivars allowed to deepen the knowledge about the underlying genetic structure. The great genetic variability of Spanish germplasm was confirmed, and a new hot spot of diversity was identified in the northern regions of La Rioja and Aragon. Analysis of the genetic structure showed a clear separation between the germplasm of southern and northern-northeastern Spain and indicated a significantly higher level of admixture in the latter. Given the expansion of modern olive cultivation with only a few cultivars, this cryptic germplasm is in great danger of disappearing. This underlines the fact that maintaining as many cultivars as possible will increase the genetic variability of the olive gene pool to meet the future challenges of olive cultivation.

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.

Gene flow between wild trees and cultivated varieties shapes the genetic structure of sweet chestnut (Castanea sativa Mill.) populations.

Gene flow between cultivated and wild gene pools is common in the contact zone between agricultural lands and natural habitats and can be used to study the development of adaptations and selection of novel varieties. This is likely the case in the northern Adriatic region, where centuries-old cultivated orchards of sweet chestnut (Castanea sativa Mill.) are planted within the natural distribution area of the species. Thus, we investigated the population structure of several orchards of sweet chestnuts. Furthermore, the genetic background of three toponymous clonal varieties was explored. Six genomic simple sequence repeat (gSSR) and nine EST-derived SSR (EST-SSR) loci were utilized in this research, and both grafted and non-grafted individuals were included in this study. Five closely related clones were identified, which represent a singular, polyclonal marron variety, found in all three cultivation areas. Furthermore, many hybrids, a result of breeding between cultivated and wild chestnuts, have been found. Analyzed semi-wild orchards defined by a diverse genetic structure, represent a hotspot for further selection and could result in creation of locally adapted, high-yielding varieties.

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.

Chloroplast Genome Annotation Tools: Prolegomena to the Identification of Inverted Repeats.

The development of next-generation sequencing technology and the increasing amount of sequencing data have brought the bioinformatic tools used in genome assembly into focus. The final step of the process is genome annotation, which works on assembled genome sequences to identify the location of genome features. In the case of organelle genomes, specialized annotation tools are used to identify organelle genes and structural features. Numerous annotation tools target chloroplast sequences. Most chloroplast DNA genomes have a quadripartite structure caused by two copies of a large inverted repeat. We investigated the strategies of six annotation tools (Chloë, Chloroplot, GeSeq, ORG.Annotate, PGA, Plann) for identifying inverted repeats and analyzed their success using publicly available complete chloroplast sequences of taxa belonging to the asterid and rosid clades. The annotation tools use two different approaches to identify inverted repeats, using existing general search tools or implementing stand-alone solutions. The chloroplast sequences studied show that there are different types of imperfections in the assembled data and that each tool performs better on some sequences than the others.

Development of Microsatellite Markers for Tanacetum cinerariifolium (Trevis.) Sch. Bip., a Plant with a Large and Highly Repetitive Genome.

Dalmatian pyrethrum (Tanacetum cinerariifolium (Trevis.) Sch. Bip.) is an outcrossing plant species (2n = 18) endemic to the eastern Adriatic coast and source of the natural insecticide pyrethrin. Due to the high repeatability and large genome (1C-value = 9.58 pg) our previous attempts to develop microsatellite markers using the traditional method were unsuccessful. Now we have used Illumina paired-end whole genome sequencing and developed a specific procedure to obtain useful microsatellite markers. A total of 796,130,142 high-quality reads (approx. 12.5× coverage) were assembled into 6,909,675 contigs using two approaches (de novo assembly and joining of overlapped pair-end reads). A total of 31,380 contigs contained one or more microsatellite sequences, of which di-(59.7%) and trinucleotide (25.9%) repeats were the most abundant. Contigs containing microsatellites were filtered according to various criteria to achieve better yield of functional markers. After two rounds of testing, 17 microsatellite markers were developed and characterized in one natural population. Twelve loci were selected for preliminary genetic diversity analysis of three natural populations. Neighbor-joining tree, based on the proportion of shared alleles distances, grouped individuals into clusters according to population affiliation. The availability of codominant SSR markers will allow analysis of genetic diversity and structure of natural Dalmatian pyrethrum populations as well as identification of breeding lines and cultivars.

The complete chloroplast genome of dalmatian pyrethrum (Tanacetum cinerariifolium (Trevir.) Sch. Bip. (Asteraceae)), the source of the natural insecticide pyrethrin.

Tanacetum cinerariifolium is an endemic species of the eastern Adriatic coast that synthesizes the natural insecticide pyrethrin. We have characterized the complete chloroplast genome of the species and analyzed its phylogeny within the Asteraceae family. The complete chloroplast genome of T. cinerariifolium has a size of 150,136 bp, including a large single-copy (LSC) region of 82,717 bp, a small single-copy (SSC) region of 18,411bp, and a pair of inverted repeats (IRs) of 24,504 bp. The chloroplast genome of T. cinerariifolium encodes 108 genes, including 77 protein-coding genes (PCGs), 27 tRNA genes, and 4 rRNA genes. Phylogenetic analyses based on the complete chloroplast genomes placed T. cinerariifolium in a sister position to species of the genera Artemisia and Chrysanthemum.

Crop breeding for a changing climate in the Pannonian region: towards integration of modern phenotyping tools.

The Pannonian Plain, as the most productive region of Southeast Europe, has a long tradition of agronomic production as well as agronomic research and plant breeding. Many research institutions from the agri-food sector of this region have a significant impact on agriculture. Their well-developed and fruitful breeding programmes resulted in productive crop varieties highly adapted to the specific regional environmental conditions. Rapid climatic changes that occurred during the last decades led to even more investigations of complex interactions between plants and their environments and the creation of climate-smart and resilient crops. Plant phenotyping is an essential part of botanical, biological, agronomic, physiological, biochemical, genetic, and other omics approaches. Phenotyping tools and applied methods differ among these disciplines, but all of them are used to evaluate and measure complex traits related to growth, yield, quality, and adaptation to different environmental stresses (biotic and abiotic). During almost a century-long period of plant breeding in the Pannonian region, plant phenotyping methods have changed, from simple measurements in the field to modern plant phenotyping and high-throughput non-invasive and digital technologies. In this review, we present a short historical background and the most recent developments in the field of plant phenotyping, as well as the results accomplished so far in Croatia, Hungary, and Serbia. Current status and perspectives for further simultaneous regional development and modernization of plant phenotyping are also discussed.

Utility of EST-SNP Markers for Improving Management and Use of Olive Genetic Resources: A Case Study at the Worldwide Olive Germplasm Bank of Córdoba.

Olive, the emblematic Mediterranean fruit crop, owns a great varietal diversity, which is maintained in ex situ field collections, such as the World Olive Germplasm Bank of Córdoba (WOGBC), Spain. Accurate identification of WOGBC, one of the world's largest collections, is essential for efficient management and use of olive germplasm. The present study is the first report of the use of a core set of 96 EST-SNP markers for the fingerprinting of 1273 accessions from 29 countries, including both field and new acquired accessions. The EST-SNP fingerprinting made possible the accurate identification of 668 different genotypes, including 148 detected among the new acquired accessions. Despite the overall high genetic diversity found at WOGBC, the EST-SNPs also revealed the presence of remarkable redundant germplasm mostly represented by synonymy cases within and between countries. This finding, together with the presence of homonymy cases, may reflect a continuous interchange of olive cultivars, as well as a common and general approach for their naming. The structure analysis revealed a certain geographic clustering of the analysed germplasm. The EST-SNP panel under study provides a powerful and accurate genotyping tool, allowing for the foundation of a common strategy for efficient safeguarding and management of olive genetic resources.

Phenotypic Variation in European Wild Pear (Pyrus pyraster (L.) Burgsd.) Populations in the North-Western Part of the Balkan Peninsula.

Leaves play a central role in plant fitness, allowing efficient light capture, gas exchange and thermoregulation, ensuring optimal growing conditions for the plant. Phenotypic variability in leaf shape and size has been linked to environmental heterogeneity and habitat characteristics. Therefore, the study of foliar morphology in plant populations can help us to identify the environmental factors that may have influenced the process of species diversification. In this study, we used European wild pear (Pyrus pyraster (L.) Burgsd., Rosaceae) as a model species to investigate the phenotypic variability of leaves under different environmental conditions. Using leaf morphometric data from 19 natural populations from the north-western part of the Balkan Peninsula, a high level of variability among and within populations were found. Leaf traits related to leaf size were more variable compared to leaf shape traits, with both influenced by geographic and environmental factors. Consequently, patterns of isolation by environment (IBE) and distance (IBD) were identified, with IBE showing a stronger influence on leaf variability. Multivariate statistical analysis revealed that European wild pear populations from the north-western part of the Balkan Peninsula can be divided into two morphological clusters, consistent with their geographical distance and environmental conditions. Our results confirm a high level of phenotypic variability in European wild pear populations, providing additional data on this poorly studied species, emphasizing phenotypic plasticity as a major driver in the adaptation of this noble hardwood species to rapid climate change.

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 FST - 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.

Grass pea natural variation reveals oligogenic resistance to Fusarium oxysporum f. sp. pisi.

Grass pea (Lathyrus sativus L.) is an annual legume species, phylogenetically close to pea (Pisum sativum L.), that may be infected by Fusarium oxysporum f. sp. pisi (Fop), the causal agent of fusarium wilt in peas with vast worldwide yield losses. A range of responses varying from high resistance to susceptibility to this pathogen has been reported in grass pea germplasm. Nevertheless, the genetic basis of that diversity of responses is still unknown, hampering its breeding exploitation. To identify genomic regions controlling grass pea resistance to fusarium wilt, a genome-wide association study approach was applied on a grass pea worldwide collection of accessions inoculated with Fop race 2. Disease responses were scored in this collection that was also subjected to high-throughput based single nucleotide polymorphisms (SNP) screening through genotyping-by-sequencing. A total of 5,651 high-quality SNPs were considered for association mapping analysis, performed using mixed linear models accounting for population structure. Because of the absence of a fully assembled grass pea reference genome, SNP markers' genomic positions were retrieved from the pea's reference genome v1a. In total, 17 genomic regions were associated with three fusarium wilt response traits in grass pea, anticipating an oligogenic control. Seven of these regions were located on pea chromosomes 1, 6, and 7. The candidate genes underlying these regions were putatively involved in secondary and amino acid metabolism, RNA (regulation of transcription), transport, and development. This study revealed important fusarium wilt resistance favorable grass pea SNP alleles, allowing the development of molecular tools for precision disease resistance breeding.

High diversity of natural Dalmatian pyrethrum based on pyrethrin composition at intra- and interpopulation level.

Seeds collected from 10 Dalmatian pyrethrum (Tanacetum cinerariifolium /Trevir./ Sch. Bip. Asteraceae) populations naturally occurring along the Adriatic coast and islands, were used to establish a common garden with the aim of estimating the variation of pyrethrin compounds between and within populations. A recently optimized matrix solid phase dispersion (MSPD) method was used for the extraction of six pyrethrin compounds from dry flower heads of 200 Dalmatian pyrethrum individuals. Separation and quantification of pyrethrin compounds were performed by high performance liquid chromatography with diode array detector (HPLC-DAD). High variability in pyrethrin content and composition was observed within and between populations. Total pyrethrin content in individual samples varied from 0.10% to 1.35% of flower dry weight. On average, the lowest total pyrethrin content was observed in the population of Peljesac (0.22%) and the highest in the island populations of Mali Losinj and Zlarin (0.87%). The population of Mali Losinj had the most favourable pyrethrin profile, including the lowest variability in total pyrethrin content, the highest pyrethrin I (PI) content (up to 60.47%) and PI/PII ratio (up to 5.88), as well as the lowest pyrethrin II (PII) content. In contrast, the populations from Biokovo and Peljesac represented pyrethrin profiles with the lowest content of total pyrethrin, PI and PI/PII ratio. The Biokovo population also had the highest PII content (43.18% on average), while the Peljesac population had the highest jasmolin I (JI) and jasmolin II (JII) content. Four different chemotypes were determined by cluster analysis. Chemotype P1 is the most promising for future breeding programs and is characterized by the highest total pyrethrin content, PI content and PI/PII ratio, as important measures of insecticidal activity. Correlation analysis of each pyrethrin compound with 25 environmental variables and Principal Component Analysis (PCA) were carried out and the importance of optimum temperature and precipitation on pyrethrin quantity and quality was suggested. Total pyrethrin content was significantly correlated with six temperature and five precipitation parameters. Pyrethrin I, PII and cinerin II (CII) were significantly correlated mainly with temperature parameters; JII and cinerin I (CI) with precipitation parameters, while JI was correlated with temperature range and precipitation parameters. The mountainous Biokovo population was an outlier in the PCA analysis, mainly due to the extremely low values of temperature parameters at a higher altitude. Precipitation variables separated the central Adriatic islands inhabiting drier habitats from populations inhabiting habitats with greater precipitation. The high-resolution data from this study provide the opportunity for individual-based selection for breeding plants with the best pyrethrin profiles, as well as for increased adaptation to broader climatic conditions.