Genetic diversity and structure of Dalmatian pyrethrum (Tanacetum cinerariifolium Trevir./Sch./Bip., Asteraceae) within the Balkan refugium

Local adaptation to the native environment affects pyrethrin variability in Dalmatian pyrethrum populations

The insecticidal compound pyrethrin is synthesized in Dalmatian pyrethrum (Tanacetum cinerariifolium (Trevis.) Sch.Bip.; Asteraceae), a plant species endemic to the eastern Mediterranean. Pyrethrin is a mixture of six compounds, pyrethrin I and II, cinerin I and II, and jasmolin I and II. For this study we sampled 15 natural Dalmatian pyrethrum populations covering the entire natural distribution range of the species; Croatian coastal regions and the islands, inland Bosnia and Herzegovina and Montenegro. The plants were grown in a field experiment under uniform growing conditions to exclude a short-term response to environmental factors and instead observe variation in pyrethrin content and composition among and within populations due to genetic adaptation to the native environment. The drivers of local adaptation were explored by examining the role of bioclimatic factors as a cause of population differentiation. Pyrethrins were extracted by ultrasound-assisted extraction, and the extracts were analyzed by HPLC-UV-DAD. The populations differed significantly in the content and composition of pyrethrins. The highest levels of total pyrethrins (1.27% flower DW), were found in population P14 Budva and the significantly highest levels of pyrethrin I in population P14 Vranjske Njive, Podgorica (66.47% of total pyrethrin). Based on bioclimatic conditions of the sampling sites, populations were grouped into five bioclimatic groups (A, B, C, D, and E), which showed qualitative and quantitative variability in pyrethrin content. The most abundant bioclimatic group was bioclimatic group E, which was characterized by the highest average values for pyrethrin I (53.87% of total pyrethrin), total pyrethrin content (1.06% flower DW) and the ratio of pyrethrin I and II (1.85). The correlation analysis between the pyrethrin compounds and some of the bioclimatic variables (e. g., BIO03 Isothermality and BIO04 Temperature seasonality) showed their significant contribution in explaining the variation of pyrethrins in T. cinerariifolium. The differences in pyrethrin content and composition may be partly due to genetic adaptation to the ecological conditions of the native environment. The obtained data would enable the selection of source populations for breeding programs aimed at producing cultivars with desirable biochemical properties and adaptation to different bioclimatic conditions.

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.

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.

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.

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 Pelješac (0.22%) and the highest in the island populations of Mali Lošinj and Zlarin (0.87%). The population of Mali Lošinj 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 Pelješac 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 Pelješac 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.

Xylem cavitation isolates leaky flowers during water stress in pyrethrum

Flowers underpin plant evolution, genetic legacy and global food supply. They are exposed to similar evaporative conditions as leaves, yet floral physiology is a product of different selective forces. We used Tanacetum cinerariifolium, a perennial daisy, to examine the response of flowers to whole-plant water stress, determining if flowers constitute a liability during drought, and how this species has adapted to minimize risk associated with reproduction. We determined the relative transpiration cost of flowers and leaves and confirmed that flowers in this species are xylem-hydrated. The relative water stress tolerance of leaves and flowers then was compared using xylem vulnerability measurements linked with observed tissue damage during an acute drought treatment. Flowers were a major source of water loss during drought but the xylem supplying them was much more vulnerable to cavitation than leaves. This xylem vulnerability segmentation was confirmed by observations that most flowers died whereas leaves were minimally affected during drought. Early cavitation and hydraulic isolation of flowers during drought benefits the plant by slowing the dehydration of perennial vegetative organs and delaying systemic xylem damage. Our results highlight the need to understand flower xylem vulnerability as a means of predicting plant reproductive failure under future drought.

Genetic diversity and population structure of the sweet leaf herb, Stevia rebaudiana B., cultivated and landraces germplasm assessed by EST-SSRs genotyping and steviol glycosides phenotyping

BACKGROUND

Stevia rebaudiana (Asteraceae), native from Paraguay, accumulates steviol glycosides (SGs) into its leaves. These compounds exhibit acaloric intense sweet taste which answers to consumer demands for reducing daily sugar intake. Despite the developpement of S. rebaudiana cultivation all over the world, the development of new cultivars is very recent, in particular due to a colossal lack of (1) germplasm collection and breeding, (2) studies on genetic diversity and its structuring, (3) genomic tools.

RESULTS

In this study, we developped 18 EST-SSR from 150,258 EST from The Compositae Genome Project of UC Davis ( http://compgenomics.ucdavis.edu/data/ ). We genotyped 145 S. rebaudiana individuals, issued from thirty-one cultivars and thirty-one landraces of various origins worldwide. Markers polymorphic information content (PIC) ranged between 0.60 and 0.84. An average of 12 alleles per locus and a high observed heterozygoty of 0.69 could be observed. The landraces revealed twice as many private alleles as cultivars. The genotypes could be clustered into 3 genetic populations. The landraces were grouped in the same cluster in which the oldest cultivars "Eirete" and "MoritaIII" type are also found. The other two clusters only include cultivated genotypes. One of them revealed an original genetic variability. SG phenotypes could not discriminate the three genetic clusters but phenotyping showed a wide range of composition in terms of bitter to sweet SGs.

CONCLUSION

This is the first study of genetic diversity in Stevia rebaudiana involving 145 genotypes, including known cultivars as well as landrace populations of different origin. This study pointed out the structuration of S. rebaudiana germplasm and the resource of the landrace populations for genetic improvement, even on the trait of SG's composition.

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