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.

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.

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.

Quantitatively Unraveling Hierarchy of Factors Impacting Virgin Olive Oil Phenolic Profile and Oxidative Stability.

A single phenolic group and even a compound play different roles in the sensory properties and stability of virgin olive oil (VOO), which in turn are strongly influenced by several factors. Understanding the causes of differences in phenolic compound composition and oxidative stability (OS) in VOOs is essential for targeted and timely harvest and processing while maintaining desired oil quality. The phenolic profile and OS of two monocultivar VOOs (Oblica and Leccino) grown in two geographical sites of different altitudes (coastal plain and hilly hinterland) were analyzed throughout the ripening period over two years. Concentration of secoiridoids was 30% higher in the Oblica than in the Leccino VOOs, which in turn had significantly higher values of OS. Both cultivars had more than twice as high concentrations of the two most abundant phenolic compounds, the dialdehyde form of decarboxymethyl oleuropein aglycone and the dialdehyde form of decarboxymethyl ligstroside aglycone, and OS values in a colder growing site of higher altitude. Among the studied monocultivar VOOs, the secoiridoid group did not behave equally during ripening. The hierarchy of different influencing factors was investigated using multivariate statistics and revealed: cultivar > geographical site > harvest period > growing season. In addition, the possibility of traceability of VOO using molecular markers was investigated by establishing SSR profiles of oils of the studied cultivars and comparing them with SSR profiles of leaves.

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.

Morphological and Chemical Diversity and Antioxidant Capacity of the Service Tree (Sorbus domestica L.) Fruits from Two Eco-Geographical Regions.

Service tree, Sorbus domestica L., is a rare and neglected wild fruit tree species of southern and central Europe. Being distributed in different eco-geographical regions, with fragmented and low-density populations, S. domestica represents an interesting model case for investigating patterns of within- and between-population diversity at geographical and environmental scales. This study aimed to analyze the proximate composition, antioxidant activity, and morphometric fruit characteristics. We examined the diversity and population divergences of 49 S. domestica individuals originating from seven populations across continental and Mediterranean eco-geographical regions. In addition, tests of isolation by distance and environment were performed to detect the magnitude of divergence explained by geographic and environmental variables. Significant differences between the studied populations were found in almost all of the studied morphometric and chemical fruit characteristics. The studied service tree populations were characterized by high phenotypic variation despite the low number of trees per population. Model-based population structure analysis using morphometric and chemical fruit characteristics revealed three groups of service tree populations. We concluded that non-effective pollen and seed dispersal along with genetic drift and specific environmental factors resulted in a distinct phenotype with a specific chemical composition in the isolated island population. In addition, a pattern of isolation by the environment was revealed. We infer that morphological and chemical differences between the studied populations in the true service tree from different eco-geographical regions were mediated by adaptation to the specific environmental conditions.

Towards the Well-Tempered Chloroplast DNA Sequences.

With the development of next-generation sequencing technology and bioinformatics tools, the process of assembling DNA sequences has become cheaper and easier, especially in the case of much shorter organelle genomes. The number of available DNA sequences of complete chloroplast genomes in public genetic databases is constantly increasing and the data are widely used in plant phylogenetic and biotechnological research. In this work, we investigated possible inconsistencies in the stored form of publicly available chloroplast genome sequence data. The impact of these inconsistencies on the results of the phylogenetic analysis was investigated and the bioinformatic solution to identify and correct inconsistencies was implemented. The whole procedure was demonstrated using five plant families (Apiaceae, Asteraceae, Campanulaceae, Lamiaceae and Rosaceae) as examples.

Assessment of the Origin and Diversity of Croatian Common Bean Germplasm Using Phaseolin Type, SSR and SNP Markers and Morphological Traits.

Landraces represent valuable genetic resources for breeding programmes to produce high-yielding varieties adapted to stressful environmental conditions. Although the common bean (Phaseolus vulgaris L.) is an economically important food legume for direct human consumption worldwide, common bean production in Croatia is based almost exclusively on landraces and there is no common bean breeding program. Information on phaseolin type and results of population structure and genetic diversity obtained by analysis of SSR and SNP markers, in combination with the morphological characterization of 174 accessions of 10 common bean landraces (morphotypes), enabled thorough classification of accessions. The accessions were classified into phaseolin type H1 ("S") of Mesoamerican origin and phaseolin types H2 ("H" or "C") and H3 ("T") of Andean origin. By applying distance- and model-based clustering methods to SSR markers, the accessions were classified into two clusters at K = 2 separating the accessions according to the centres of origin, while at K = 3, the accessions of Andean origin were further classified into two clusters of accessions that differed in phaseolin type (H2 and H3). Using SNP markers, model-based analysis of population structure was performed, the results of which were consistent with those of SSR markers. In addition, 122 accessions were assigned to 14 newly formed true-type morphogenetic groups derived from three different domestication events: (1) Mesoamerican (H1A) ("Biser", "Kukuruzar", "Tetovac", "Tresnjevac"), (2) Andean-indeterminate type (H2B1) ("Dan noc", "Sivi", "Puter", "Sivi prosarani", "Tresnjevac") and (3) Andean-determinate type (H3B2) ("Bijeli", "Dan noc", "Puter", "Tresnjevac", "Zelencec"). The rest of the accessions could represent putative hybrids between morphogenetic groups. The differences between the true-type groups of accessions were further analysed based on nine quantitative traits, and the subsets of traits that best distinguish among centres of origin (A: Mesoamerican, B: Andean) and genetic groups (H1A, H2B1, H3B2) were proposed.

Genome-Wide Association Studies of Mineral Content in Common Bean.

Micronutrient malnutrition is one of the main public health problems in many parts of the world. This problem raises the attention of all valuable sources of micronutrients for the human diet, such as common bean (Phaseolus vulgaris L.). In this research, a panel of 174 accessions representing Croatian common bean landraces was phenotyped for seed content of eight nutrients (N, P, K, Ca, Mg, Fe, Zn, and Mn), and genotyped using 6,311 high-quality DArTseq-derived SNP markers. A genome-wide association study (GWAS) was then performed to identify new genetic sources for improving seed mineral content. Twenty-two quantitative trait nucleotides (QTN) associated with seed nitrogen content were discovered on chromosomes Pv01, Pv02, Pv03, Pv05, Pv07, Pv08, and Pv10. Five QTNs were associated with seed phosphorus content, four on chromosome Pv07, and one on Pv08. A single significant QTN was found for seed calcium content on chromosome Pv09 and for seed magnesium content on Pv08. Finally, two QTNs associated with seed zinc content were identified on Pv06 while no QTNs were found to be associated with seed potassium, iron, or manganese content. Our results demonstrate the utility of GWAS for understanding the genetic architecture of seed nutritional traits in common bean and have utility for future enrichment of seed with macro- and micronutrients through genomics-assisted breeding.

The Influence of a Seedling Recruitment Strategy and a Clonal Architecture on a Spatial Genetic Structure of a Salvia brachyodon (Lamiaceae) Population.

By performing a high-resolution spatial-genetic analysis of a partially clonal Salvia brachyodon population, we elucidated its clonal architecture and seedling recruitment strategy. The sampling of the entire population was based on a 1 × 1 m grid and each sampled individual was genotyped. Population-genetic statistics were combined with geospatial analyses. On the population level, the presence of both sexual and clonal reproduction and repeated seedling recruitment as the prevailing strategy of new genets establishment were confirmed. On the patch level, a phalanx clonal architecture was detected. A significant negative correlation between patches' sizes and genotypic richness was observed as young plants were not identified within existing patches of large genets but almost exclusively in surrounding areas. The erosion of the genetic variability of older patches is likely caused by the inter-genet competition and resulting selection or by a random die-off of individual genets accompanied by the absence of new seedlings establishment. This study contributes to our understanding of how clonal architecture and seedling recruitment strategies can shape the spatial-genetic structure of a partially clonal population and lays the foundation for the future research of the influence of the population's clonal organization on its sexual reproduction.

Divergent selection and genetic structure of Sideritis scardica populations from southern Balkan Peninsula as revealed by AFLP fingerprinting.

Sideritis scardica Giseb. is a subalpine/alpine plant species endemic to the central part of the Balkan Peninsula. In this study, we combined Amplified Fragment Length Polymorphism (AFLP) and environmental data to examine the adaptive genetic variations in S. scardica natural populations sampled in contrasting environments. A total of 226 AFLP loci were genotyped in 166 individuals from nine populations. The results demonstrated low gene diversity, ranging from 0.095 to 0.133 and significant genetic differentiation ranging from 0.115 to 0.408. Seven genetic clusters were revealed by Bayesian clustering methods as well as by Discriminant Analysis of Principal Components and each population formed its respective cluster. The exception were populations P02 Mt. Shara and P07 Mt. Vermio, that were admixed between two clusters. Both landscape genetic methods Mcheza and BayeScan identified a total of seven (3.10%) markers exhibiting higher levels of genetic differentiation among populations. The spatial analysis method Samßada detected 50 individual markers (22.12%) associated with bioclimatic variables, among them seven were identified by both Mcheza and BayeScan as being under directional selection. Four bioclimatic variables associated with five out of seven outliers were related to precipitation, suggesting that this variable is the key factor affecting the adaptive variation of S. scardica.