Epigenetic stress memory in gymnosperms.

Gymnosperms are long-lived, cone-bearing seed plants that include some of the most ancient extant plant species. These relict land plants have evolved to survive in habitats marked by chronic or episodic stress. Their ability to thrive in these environments is partly due to their phenotypic flexibility, and epigenetic regulation likely plays a crucial part in this plasticity. We review the current knowledge on abiotic and biotic stress memory in gymnosperms and the possible epigenetic mechanisms underlying long-term phenotypic adaptations. We also discuss recent technological improvements and new experimental possibilities that likely will advance our understanding of epigenetic regulation in these ancient and hard-to-study plants.

1,2,3,4-dithiadiazole derivatives as a novel class of calcium signaling modulators.

Aberrant calcium signaling is associated with a diverse range of pathologies, including cardiovascular and neurodegenerative diseases, diabetes, cancer, etc… So, therapeutic strategies based on the correction of pathological calcium signaling are becoming extremely in demand. Thus, the development of novel calcium signaling modulators remains highly actual. Previously we found that 1,2,3,4-dithiadiazole derivative 3-(4-nitrophenyl)-5-phenyl-3H-1,2,3,4-dithiadiazole-2-oxide can strongly reduce calcium uptake through store-operated calcium (SOC) channels. Here we tested several structurally related compounds and found that most of them can effectively affect SOC channels and attenuate calcium content in the endoplasmic reticulum, thus, establishing 1,2,3,4-dithiadiazoles as a novel class of SOC channel inhibitors. Comparing different 1,2,3,4-dithiadiazole derivatives we showed that previously published 3-(4-nitrophenyl)-5-phenyl-3H-1,2,3,4-dithiadiazole-2-oxide and newly tested 3-(3,5-difluorophenyl)-5-phenyl-3H-1,2,3,4-dithiadiazole 2-oxide demonstrated the highest efficacy of SOC entry reduction, supposing the important role of electron-withdrawing substituents to realize the inhibitory activity of 1,2,3,4-dithiadiazoles.

Parametric modeling and analysis of intake phases for side-ported Wankel rotary engines.

This research presents a novel port parametric modeling technique using three-dimensional computational fluid dynamics for the design and optimization of intake and exhaust phases in side-ported Wankel rotary engines (WREs). Definitions for the port phases encompass parameters such as port start opening, port full opening, port start closing, and port full closing timings. The four port phase control arcs are obtained by translating and rotating the rotor flank to satisfy the high control accuracy. Further, the shape of the port is further smoothed and varied by four auxiliary circular arcs. Moreover, the influence of port full closing timing and the size of auxiliary circular arcs (R1, and R3) on the intake characteristics is studied. The results show that the novel method can flexibly and effectively control the phases and shapes. The early port full closing timing reduces fluid backflow and improves volumetric efficiency (VE) but increases intake loss (IL). The small size of R1 facilitates to increase the VE and reduce IL. A larger or smaller size of R3 is not conducive to reducing IL, and the smaller size of R3 improves the VE. The novel generation method proposed in this paper provides a theoretical basis to optimize the design of various sizes of side-ported WREs and guidance for practical manufacturing.

Warmer temperature during asexual reproduction induce methylome, transcriptomic, and lasting phenotypic changes in Fragaria vesca ecotypes.

Plants must adapt with increasing speed to global warming to maintain their fitness. One rapid adaptation mechanism is epigenetic memory, which may provide organisms sufficient time to adapt to climate change. We studied how the perennial Fragaria vesca adapted to warmer temperatures (28°C vs. 18°C) over three asexual generations. Differences in flowering time, stolon number, and petiole length were induced by warmer temperature in one or more ecotypes after three asexual generations and persisted in a common garden environment. Induced methylome changes differed between the four ecotypes from Norway, Iceland, Italy, and Spain, but shared methylome responses were also identified. Most differentially methylated regions (DMRs) occurred in the CHG context, and most CHG and CHH DMRs were hypermethylated at the warmer temperature. In eight CHG DMR peaks, a highly similar methylation pattern could be observed between ecotypes. On average, 13% of the differentially methylated genes between ecotypes also showed a temperature-induced change in gene expression. We observed ecotype-specific methylation and expression patterns for genes related to gibberellin metabolism, flowering time, and epigenetic mechanisms. Furthermore, we observed a negative correlation with gene expression when repetitive elements were found near (±2 kb) or inside genes. In conclusion, lasting phenotypic changes indicative of an epigenetic memory were induced by warmer temperature and were accompanied by changes in DNA methylation patterns. Both shared methylation patterns and transcriptome differences between F. vesca accessions were observed, indicating that DNA methylation may be involved in both general and ecotype-specific phenotypic variation.

Epitype-inducing temperatures drive DNA methylation changes during somatic embryogenesis in the long-lived gymnosperm Norway spruce.

An epigenetic memory of the temperature sum experienced during embryogenesis is part of the climatic adaptation strategy of the long-lived gymnosperm Norway spruce. This memory has a lasting effect on the timing of bud phenology and frost tolerance in the resulting epitype trees. The epigenetic memory is well characterized phenotypically and at the transcriptome level, but to what extent DNA methylation changes are involved have not previously been determined. To address this, we analyzed somatic epitype embryos of Norway spruce clones produced at contrasting epitype-inducing conditions (18 and 28°C). We screened for differential DNA methylation in 2744 genes related mainly to the epigenetic machinery, circadian clock, and phenology. Of these genes, 68% displayed differential DNA methylation patterns between contrasting epitype embryos in at least one methylation context (CpG, CHG, CHH). Several genes related to the epigenetic machinery (e.g., DNA methyltransferases, ARGONAUTE) and the control of bud phenology (FTL genes) were differentially methylated. This indicates that the epitype-inducing temperature conditions induce an epigenetic memory involving specific DNA methylation changes in Norway spruce.

Major transcriptomic differences are induced by warmer temperature conditions experienced during asexual and sexual reproduction in Fragaria vesca ecotypes.

A major challenge for plants in a rapidly changing climate is to adapt to rising temperatures. Some plants adapt to temperature conditions by generating an epigenetic memory that can be transmitted both meiotically and mitotically. Such epigenetic memories may increase phenotypic variation to global warming and provide time for adaptation to occur through classical genetic selection. The goal of this study was to understand how warmer temperature conditions experienced during sexual and asexual reproduction affect the transcriptomes of different strawberry (Fragaria vesca) ecotypes. We let four European F. vesca ecotypes reproduce at two contrasting temperatures (18 and 28°C), either asexually through stolon formation for several generations, or sexually by seeds (achenes). We then analyzed the transcriptome of unfolding leaves, with emphasis on differential expression of genes belonging to the epigenetic machinery. For asexually reproduced plants we found a general transcriptomic response to temperature conditions but for sexually reproduced plants we found less significant responses. We predicted several splicing isoforms for important genes (e.g. a SOC1, LHY, and SVP homolog), and found significantly more differentially presented splicing event variants following asexual vs. sexual reproduction. This difference could be due to the stochastic character of recombination during meiosis or to differential creation or erasure of epigenetic marks during embryogenesis and seed development. Strikingly, very few differentially expressed genes were shared between ecotypes, perhaps because ecotypes differ greatly both genetically and epigenetically. Genes related to the epigenetic machinery were predominantly upregulated at 28°C during asexual reproduction but downregulated after sexual reproduction, indicating that temperature-induced change affects the epigenetic machinery differently during the two types of reproduction.

Methylome, transcriptome, and phenotype changes induced by temperature conditions experienced during sexual reproduction in Fragaria vesca.

Temperature conditions experienced during embryogenesis and seed development may induce epigenetic changes that increase phenotypic variation in plants. Here we investigate if embryogenesis and seed development at two different temperatures (28 vs. 18°C) result in lasting phenotypic effects and DNA methylation changes in woodland strawberry (Fragaria vesca). Using five European ecotypes from Spain (ES12), Iceland (ICE2), Italy (IT4), and Norway (NOR2 and NOR29), we found statistically significant differences between plants from seeds produced at 18 or 28°C in three of four phenotypic features investigated under common garden conditions. This indicates the establishment of a temperature-induced epigenetic memory-like response during embryogenesis and seed development. The memory effect was significant in two ecotypes: in NOR2 flowering time, number of growth points and petiole length were affected, and in ES12 number of growth points was affected. This indicates that genetic differences between ecotypes in their epigenetic machinery, or other allelic differences, impact this type of plasticity. We observed statistically significant differences between ecotypes in DNA methylation marks in repetitive elements, pseudogenes, and genic elements. Leaf transcriptomes were also affected by embryonic temperature in an ecotype-specific manner. Although we observed significant and lasting phenotypic change in at least some ecotypes, there was considerable variation in DNA methylation between individual plants within each temperature treatment. This within-treatment variability in DNA methylation marks in F. vesca progeny may partly be a result of allelic redistribution from recombination during meiosis and subsequent epigenetic reprogramming during embryogenesis.

Plant-made vaccines against viral diseases in humans and farm animals.

Plants provide not only food and feed, but also herbal medicines and various raw materials for industry. Moreover, plants can be green factories producing high value bioproducts such as biopharmaceuticals and vaccines. Advantages of plant-based production platforms include easy scale-up, cost effectiveness, and high safety as plants are not hosts for human and animal pathogens. Plant cells perform many post-translational modifications that are present in humans and animals and can be essential for biological activity of produced recombinant proteins. Stimulated by progress in plant transformation technologies, substantial efforts have been made in both the public and the private sectors to develop plant-based vaccine production platforms. Recent promising examples include plant-made vaccines against COVID-19 and Ebola. The COVIFENZ® COVID-19 vaccine produced in Nicotiana benthamiana has been approved in Canada, and several plant-made influenza vaccines have undergone clinical trials. In this review, we discuss the status of vaccine production in plants and the state of the art in downstream processing according to good manufacturing practice (GMP). We discuss different production approaches, including stable transgenic plants and transient expression technologies, and review selected applications in the area of human and veterinary vaccines. We also highlight specific challenges associated with viral vaccine production for different target organisms, including lower vertebrates (e.g., farmed fish), and discuss future perspectives for the field.

Chromone-Containing Allylmorpholines Influence Ion Channels in Lipid Membranes via Dipole Potential and Packing Stress.

Herein, we report that chromone-containing allylmorpholines can affect ion channels formed by pore-forming antibiotics in model lipid membranes, which correlates with their ability to influence membrane boundary potential and lipid-packing stress. At 100 µg/mL, allylmorpholines 1, 6, 7, and 8 decrease the boundary potential of the bilayers composed of palmitoyloleoylphosphocholine (POPC) by about 100 mV. At the same time, the compounds do not affect the zeta-potential of POPC liposomes, but reduce the membrane dipole potential by 80-120 mV. The allylmorpholine-induced drop in the dipole potential produce 10-30% enhancement in the conductance of gramicidin A channels. Chromone-containing allylmorpholines also affect the thermotropic behavior of dipalmytoylphosphocholine (DPPC), abolishing the pretransition, lowering melting cooperativity, and turning the main phase transition peak into a multicomponent profile. Compounds 4, 6, 7, and 8 are able to decrease DPPC's melting temperature by about 0.5-1.9 °C. Moreover, derivative 7 is shown to increase the temperature of transition of palmitoyloleoylphosphoethanolamine from lamellar to inverted hexagonal phase. The effects on lipid-phase transitions are attributed to the changes in the spontaneous curvature stress. Alterations in lipid packing induced by allylmorpholines are believed to potentiate the pore-forming ability of amphotericin B and gramicidin A by several times.

Dithiadiazole derivative 3-(4-nitrophenyl)-5-phenyl-3H-1,2,3,4-dithiadiazole-2-oxide - Novel modulator of store-operated calcium entry.

Pathological calcium homeostasis accompanies the development of a large number of different diseases, therefore, the search for new modulators of calcium signaling remains highly actual. Last decades store-operated calcium channels have been repeatedly postulated as a therapeutic target, so the compounds acting on them can be considered promising drug prototypes. Here, we tested several derivatives of 1,2,3,4-dithiadiazole, 1,3-thiazine, pyrazolopyrimidine and thiohydrazides for the ability to affect the thapsigargin-induced calcium response. Using calcium imaging and the patch-clamp technique we found that dithiadiazole derivative3-(4-nitrophenyl)-5-phenyl-3H-1,2,3,4-dithiadiazole-2-oxidehad a strong inhibitory effect on store-operated calcium entry at the micromolar concentration in HEK293 cells. Moreover, incubation of the cells with this compound also resulted in the decrease of ER calcium content. Thus, we have postulated 3-(4-nitrophenyl)-5-phenyl-3H-1,2,3,4-dithiadiazole-2-oxide as a novel inhibitor of store-operated calcium entry and suggested the derivatives of 1,2,3,4-dithiadiazole as a prospective class of compounds for searching new calcium modulators.

Quantification of a New Antifungal Drug Based on 1,3,4-Thiadiazole by HPLC-ESI/MS: Method Development, Validation and Pharmacokinetic Application.

The high sensitive HPLC-ESI/MS method for quantitative determination of a new antifungal drug - 2-[(1Z)-1-(3,5-diphenyl-1,3,4-thiadiazol-2(3Н)-ylidene)methyl]-3,5-diphenyl-1,3,4-thiadiazol-3-ium chloride (TDZ) - was developed and fully validated. TDZ was separated from plasma and urine samples by acetonitrile deproteinization and extraction without time-consuming sample preparation. The reversed-phase high-performance liquid chromatography on Kromasil 100-3.5 C8 column of TDZ in isocratic elution mode using 0.03% trifluoroacetic acid : acetonitrile (65:35, v/v) at a flow rate of 0.2 mL min-1 was performed. Determination of TDZ was carried out by a positive electrospray ionization in a selected ion monitoring mode for [M+]=489 m/z. The method of absolute calibration was used for quantification of TDZ in two concentrations ranges: 100-2500 pg mL-1 and 2500-30 000 pg mL-1. The established method showed a good linearity (R=0.999 for both ranges), the limits of determination and quantification were 50 and 100 pg mL-1, respectively. The Intra- and Inter-day precision values were measured by t-Distribution and Fisher's Exact Test and were in accordance with the regulatory guidance. Low matrix effects and good recovery were found for TDZ. The present method was successfully applied to determine the pharmacokinetic parameters of TDZ by means of intravenous and oral administrations to rats at 5.0 mg kg-1 and 10.0 mg kg-1, respectively.

A Flexible Synthetic Approach to Fluorescent Chromeno[4,3-b]pyridines and Pyrano[3,2-c]chromenes from Electron-Deficient 3-Vinylchromones.

We report a flexible approach to the synthesis of phenanthrene-like heterocycles through organocatalytic ANRORC (Addition of the Nucleophile, Ring Opening, and Ring Closure) reaction of electron-deficient 3-vinylchromones with cyanoacetamide. Addition of highly basic DBU (1,8-diazabicyclo[5.4.0]undec-7-ene) or tetramethylguanidine (TMG) at 80 °C leads to chromeno[4,3-b]pyridines in good yields, whereas Et3 N at 20 °C made it possible to obtain the less accessible pyrano[3,2-c]chromenes and their 2-imines. The synthesis proceeds in mild conditions (EtOH, 20-80 °C), is versatile and applicable for a wide scope of reactants. The obtained compounds show bright fluorescence in the range 460-595 nm with high quantum yields (up to 0.84) in various solvents (MeCN, DMSO, EtOH, H2 O).

Plant-Produced Vaccines: Future Applications in Aquaculture.

Aquaculture has undergone rapid development in the past decades. It provides a large part of high-quality protein food for humans, and thus, a sustainable aquaculture industry is of great importance for the worldwide food supply and economy. Along with the quick expansion of aquaculture, the high fish densities employed in fish farming increase the risks of outbreaks of a variety of aquatic diseases. Such diseases not only cause huge economic losses, but also lead to ecological hazards in terms of pathogen spread to marine ecosystems causing infection of wild fish and polluting the environment. Thus, fish health is essential for the aquaculture industry to be environmentally sustainable and a prerequisite for intensive aquaculture production globally. The wide use of antibiotics and drug residues has caused intensive pollution along with risks for food safety and increasing antimicrobial resistance. Vaccination is the most effective and environmentally friendly approach to battle infectious diseases in aquaculture with minimal ecological impact and is applicable to most species of farmed fish. However, there are only 34 fish vaccines commercially available globally to date, showing the urgent need for further development of fish vaccines to manage fish health and ensure food safety. Plant genetic engineering has been utilized to produce genetically modified crops with desirable characteristics and has also been used for vaccine production, with several advantages including cost-effectiveness, safety when compared with live virus vaccines, and plants being capable of carrying out posttranslational modifications that are similar to naturally occurring systems. So far, plant-derived vaccines, antibodies, and therapeutic proteins have been produced for human and animal health. However, the development of plant-made vaccines for animals, especially fish, is still lagging behind the development of human vaccines. The present review summarizes the development of fish vaccines currently utilized and the suitability of the plant-production platform for fish vaccine and then addresses considerations regarding fish vaccine production in plants. Developing fish vaccines by way of plant biotechnology are significant for the aquaculture industry, fish health management, food safety, and human health.

Immunochemical Detection of Modified Species of Cytosine in Plant Tissues.

Methylated cytosine (5-methylcytosine) is the most studied epigenetic mark involved in the regulation of gene expression. Although it displays highly variable dynamics during plant ontogenesis, it is possible to gain a fine spatial perspective with immunohistochemistry techniques that use specific antibodies and fluorochromes. Besides, there are other cytosine modifications described in plants, although their biological significance is still unknown (i.e., 5-hydroxymethylcytosine, 5-formylcytosine and 5-carboxylcytosine). Here we present a standardized protocol to detect cytosine modifications in plant tissues.

1,3-Thiazine, 1,2,3,4-Dithiadiazole, and Thiohydrazide Derivatives Affect Lipid Bilayer Properties and Ion-Permeable Pores Induced by Antifungals.

Over the past decade, thiazines, thiadiazoles, and thiohydrazides have attracted increasing attention due to their sedative, antimicrobial, antiviral, antifungal, and antitumor activities. The clinical efficacy of such drugs, as well as the possibility of developing resistance to antimicrobials, will depend on addressing a number of fundamental problems, including the role of membrane lipids during their interaction with plasma membranes. The effects of the eight 1,3- thiazine-, 1,2,3,4- dithiadiazole-, and thiohydrazide-related compounds on the physical properties of model lipid membranes and the effects on reconstituted ion channels induced by the polyene macrolide antimycotic nystatin and antifungal cyclic lipopeptides syringomycin E and fengycin were observed. We found that among the tested agents, the fluorine-containing compound N'-(3,5-difluorophenyl)-benzenecarbothiohydrazide (C6) was the most effective at increasing the electric barrier for anion permeation into the hydrophobic region of the membrane and reducing the conductance of anion-permeable syringomycin pores. A decrease in the membrane boundary potential with C6 adsorption also facilitated the immersion of positively charged syringomycin molecules into the lipid bilayer and increases the pore-forming ability of the lipopeptide. Using differential scanning microcalorimetry, we showed that C6 led to disordering of membrane lipids, possibly by potentiating positive curvature stress. Therefore, we used C6 as an agonist of antifungals forming the pores that are sensitive to membrane curvature stress and lipid packing, i.e., nystatin and fengycin. The dramatic increase in transmembrane current induced by syringomycin E, nystatin, and fengycin upon C6 treatment suggests its potential in combination therapy for treating invasive fungal infections.

Transcriptome profiling and in silico detection of the antimicrobial peptides of red king crab Paralithodes camtschaticus.

Endogenous antimicrobial peptides (AMPs) are evolutionarily ancient factors of innate immunity, which are produced by all multicellular organisms and play a key role in their protection against infection. Red king crab (Paralithodes camtschaticus), also called Kamchatka crab, is widely distributed and the best known species of all king crabs belonging to the family Lithodidae. Despite their economic importance, the genetic resources of king crabs are scarcely known and no full-genome sequences are available to date. Therefore, analysis of the red king crab transcriptome and identification and characterization of its AMPs could potentially contribute to the development of novel antimicrobial drug candidates when antibiotic resistance has become a global health threat. In this study, we sequenced the P. camtschaticus transcriptomes from carapace, tail flap and leg tissues using an Illumina NGS platform. Libraries were systematically analyzed for gene expression profiles along with AMP prediction. By an in silico approach using public databases we defined 49 cDNAs encoding for AMP candidates belonging to diverse families and functional classes, including buforins, crustins, paralithocins, and ALFs (anti-lipopolysaccharide factors). We analyzed expression patterns of 27 AMP genes. The highest expression was found for Paralithocin 1 and Crustin 3, with more than 8,000 reads. Other paralithocins, ALFs, crustins and ubiquicidins were among medium expressed genes. This transcriptome data set and AMPs provide a solid baseline for further functional analysis in P. camtschaticus. Results from the current study contribute also to the future application of red king crab as a bio-resource in addition to its being a known seafood delicacy.

Mass spectrometry reveals the presence of specific set of epigenetic DNA modifications in the Norway spruce genome.

5-Methylcytosine (5mC) is an epigenetic modification involved in regulation of gene expression in metazoans and plants. Iron-(II)/α-ketoglutarate-dependent dioxygenases can oxidize 5mC to 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC) and 5-carboxylcytosine (5caC). Although these oxidized forms of 5mC may serve as demethylation intermediates or contribute to transcriptional regulation in animals and fungi, experimental evidence for their presence in plant genomes is ambiguous. Here, employing reversed-phase HPLC coupled with sensitive mass spectrometry, we demonstrated that, unlike 5caC, both 5hmC and 5fC are detectable in non-negligible quantities in the DNA of a conifer, Norway spruce. Remarkably, whereas 5hmC content of spruce DNA is approximately 100-fold lower relative to human colorectal carcinoma cells, the levels of both - 5fC and a thymine base modification, 5-hydroxymethyluracil, are comparable in these systems. We confirmed the presence of modified DNA bases by immunohistochemistry in Norway spruce buds based on peroxidase-conjugated antibodies and tyramide signal amplification. Our results reveal the presence of specific range of noncanonical DNA bases in conifer genomes implying potential roles for these modifications in plant development and homeostasis.

Comparative sensitivity to gamma radiation at the organismal, cell and DNA level in young plants of Norway spruce, Scots pine and Arabidopsis thaliana.

MAIN CONCLUSION: Persistent DNA damage in gamma-exposed Norway spruce, Scots pine and Arabidopsis thaliana, but persistent adverse effects at the organismal and cellular level in the conifers only. Gamma radiation emitted from natural and anthropogenic sources may have strong negative impact on plants, especially at high dose rates. Although previous studies implied different sensitivity among species, information from comparative studies under standardized conditions is scarce. In this study, sensitivity to gamma radiation was compared in young seedlings of the conifers Scots pine and Norway spruce and the herbaceous Arabidopsis thaliana by exposure to 60Co gamma dose rates of 1-540 mGy h-1 for 144 h, as well as 360 h for A. thaliana. Consistent with slightly less prominent shoot apical meristem, in the conifers growth was significantly inhibited with increasing dose rate ≥ 40 mGy h-1. Post-irradiation, the conifers showed dose-rate-dependent inhibition of needle and root development consistent with increasingly disorganized apical meristems with increasing dose rate, visible damage and mortality after exposure to ≥ 40 mGy h-1. Regardless of gamma duration, A. thaliana showed no visible or histological damage or mortality, only delayed lateral root development after ≥ 100 mGy h-1 and slightly, but transiently delayed post-irradiation reproductive development after ≥ 400 mGy h-1. In all species dose-rate-dependent DNA damage occurred following ≥ 1-10 mGy h-1 and was still at a similar level at day 44 post-irradiation. In conclusion, the persistent DNA damage (possible genomic instability) following gamma exposure in all species may suggest that DNA repair is not necessarily mobilized more extensively in A. thaliana than in Norway spruce and Scots pine, and the far higher sensitivity at the organismal and cellular level in the conifers indicates lower tolerance to DNA damage than in A. thaliana.

In Silico Analysis of Small RNAs Suggest Roles for Novel and Conserved miRNAs in the Formation of Epigenetic Memory in Somatic Embryos of Norway Spruce.

Epigenetic memory in Norway spruce affects the timing of bud burst and bud set, vitally important adaptive traits for this long-lived forest species. Epigenetic memory is established in response to the temperature conditions during embryogenesis. Somatic embryogenesis at different epitype inducing (EpI) temperatures closely mimics the natural processes of epigenetic memory formation in seeds, giving rise to epigenetically different clonal plants in a reproducible and predictable manner, with respect to altered bud phenology. MicroRNAs (miRNAs) and other small non-coding RNAs (sRNAs) play an essential role in the regulation of plant gene expression and may affect this epigenetic mechanism. We used NGS sequencing and computational in silico methods to identify and profile conserved and novel miRNAs among small RNAs in embryogenic tissues of Norway spruce at three EpI temperatures (18, 23 and 28°C). We detected three predominant classes of sRNAs related to a length of 24 nt, followed by a 21-22 nt class and a third 31 nt class of sRNAs. More than 2100 different miRNAs within the prevailing length 21-22 nt were identified. Profiling these putative miRNAs allowed identification of 1053 highly expressed miRNAs, including 523 conserved and 530 novels. 654 of these miRNAs were found to be differentially expressed (DEM) depending on EpI temperature. For most DEMs, we defined their putative mRNA targets. The targets represented mostly by transcripts of multiple-repeats proteins, like TIR, NBS-LRR, PPR and TPR repeat, Clathrin/VPS proteins, Myb-like, AP2, etc. Notably, 124 DE miRNAs targeted 203 differentially expressed epigenetic regulators. Developing Norway spruce embryos possess a more complex sRNA structure than that reported for somatic tissues. A variety of the predicted miRNAs showed distinct EpI temperature dependent expression patterns. These putative EpI miRNAs target spruce genes with a wide range of functions, including genes known to be involved in epigenetic regulation, which in turn could provide a feedback process leading to the formation of epigenetic marks. We suggest that TIR, NBS and LRR domain containing proteins could fulfill more general functions for signal transduction from external environmental stimuli and conversion them into molecular response. Fine-tuning of the miRNA production likely participates in both developmental regulation and epigenetic memory formation in Norway spruce.

The epigenetic memory of temperature during embryogenesis modifies the expression of bud burst-related genes in Norway spruce epitypes.

MAIN CONCLUSION: Epigenetic memory affects the timing of bud burst phenology and the expression of bud burst-related genes in genetically identical Norway spruce epitypes in a manner usually associated with ecotypes. In Norway spruce, a temperature-dependent epigenetic memory established during embryogenesis affects the timing of bud burst and bud set in a reproducible and predictable manner. We hypothesize that the clinal variation in these phenological traits, which is associated with adaptation to growth under frost-free conditions, has an epigenetic component. In Norway spruce, dehydrins (DHNs) have been associated with extreme frost tolerance. DHN transcript levels decrease gradually prior to flushing, a time when trees are highly sensitive to frost. Furthermore, EARLY BUD BREAK 1 genes (EBB1) and the FT-TFL1-LIKE 2-gene (PaFTL2) were previously suggested to be implied in control of bud phenology. Here we report an analysis of transcript levels of 12 DHNs, 3 EBB1 genes and FTL2 in epitypes of the same genotype generated at different epitype-inducing temperatures, before and during spring bud burst. Earlier flushing of epitypes originating from embryos developed at 18 °C as compared to 28 °C, was associated with differential expression of these genes between epitypes and between buds and last year's needles. The majority of these genes showed significantly different expressions between epitypes in at least one time point. The general trend in DHN expression pattern in buds showed the expected reduction in transcript levels when approaching flushing, whereas, surprisingly, transcript levels peaked later in needles, mainly at the moment of bud burst. Collectively, our results demonstrate that the epigenetic memory of temperature during embryogenesis affects bud burst phenology and expression of the bud burst-related DHN, EBB1 and FTL2 genes in genetically identical Norway spruce epitypes.