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.

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.

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.

Tracing the oomycete pathogen Saprolegnia parasitica in aquaculture and the environment.

Saprolegnia parasitica causes saprolegniosis, a disease responsible for significant economic losses in aquaculture and declines of fish populations in the wild, but the knowledge of its distribution and prevalence in the environment is limited. We developed a fast, sensitive and specific S. parasitica droplet digital PCR (ddPCR) assay and demonstrated its applicability for the detection and quantification of the pathogen in environmental samples: swab DNA collected from the host (trout skin, surface of eggs) and environmental DNA extracted from water. The developed assay was used to assess how abiotic (i.e. physico-chemical parameters of the water) and biotic (health status of the host) factors influence the S. parasitica load in the environment. The pathogen load in water samples was positively correlated with some site-specific abiotic parameters such as electrical conductivity (EC) and calcium, while fluorides were negatively correlated, suggesting that physico-chemical parameters are important for determining S. parasitica load in natural waters. Furthermore, skin swabs of injured trout had significantly higher pathogen load than swabs collected from healthy fish, confirming that S. parasitica is a widespread opportunistic pathogen. Our results provide new insights into various environmental factors that influence the distribution and abundance of S. parasitica.

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.

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.

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.

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.

Octahedral-Tetrahedral Systems [Co( dppm O, O)3]2+[CoX4]2- Showing Slow Magnetic Relaxation with Two Relaxation Modes.

Three compounds with octahedral-tetrahedral Co(II) moieties of [Co( dppm O, O)3][CoX4] type, where X = SCN (1), Cl (2), or I (4) have been synthesized and characterized by the X-ray structure analysis (1 and 4), and spectroscopic methods. The dc magnetic measurements show high magnetic anisotropy for octahedral centers whereas tetrahedral sites possess moderate D values. These results are confirmed by the ab initio calculations. The ac susceptibility data reveals a slow magnetic relaxation for 2 and 4, similar to that of the X = Br analogue (3), whereas 1 displays no ac-absorption signal. There are two relaxation channels; the slower for 2 (4) possesses a relaxation time as long as τLF= 178 (588) ms at T = 1.9 K and Bdc = 0.7 T. Also, the half-Zn analogue, [Co( dppm O, O)3][ZnI4], shows slow magnetic relaxation with two relaxation channels conditioned by the cationic unit [Co( dppm O, O)3]2+.