Prevalence of analgesic use in patients with chronic pain referred to a multidisciplinary pain centre and its correlation with patient-reported outcome measures: A cross-sectional, registry-based study.

BACKGROUND: Chronic pain is prevalent in Sweden, nearing 20% in the adult population. Treatment often requires a multimodal approach, with medication, physical therapy and psychological interventions. However, the frequency of medication in patients with chronic pain in Sweden, and its correlation with patient-reported outcome measures (PROMs), are currently unknown. OBJECTIVES: To investigate the frequency of use of analgesics and other medication in patients with chro-nic pain referred to a multidisciplinary pain centre, and how opioid treatment relates to PROMs. DESIGN: Cross-sectional, registry-based study. PATIENTS: New referral visits (n = 1,275) to the Pain and Rehabilitation Center in Linköping, Sweden in 2015. 441 patients had complete medication and PROM data. METHODS: Patient-reported analgesic and other medications were matched with patient PROM data from the Swedish Quality Registry for Pain Rehabilitation. Univariate analysis was conducted with IBM Statistical Package for the Social Sciences (SPSS; IBM Corporation, Somers, NY, USA) version 24.0, and multivariate analysis with SIMCA-P+ (version 13, Umetrics AB, Umeå, Sweden), with a special emphasis on opioids. RESULTS: n = 132 (30%) patients used opioids daily, and this group differed from other patients on many PROMs, with medium effect sizes for pain severity, interference, health-related quality of life, activity engagement, and satisfaction with social life. Multivariate analysis identified four groups and showed that daily use of opioids was significantly correlated with high pain intensity and low physical functioning. CONCLUSION: Prevalence of daily opioid use was 30% and daily opioid use did not correlate with better outcome of PROMs. Longitudinal studies are warranted (e.g. on the clinical effect of tapering), as are studies that can better explain the medication variability in patients with complex chronic pain.

Salivary beta-endorphin and substance P are not biomarkers of neuropathic chronic pain propensity.

OBJECTIVE: The pathophysiology of chronic pain is complex, with most of our knowledge being derived from preclinical studies. The search for biomarkers mirroring the pathophysiology of chronic pain is ongoing, and there is an increasing interest in saliva as a diagnostic tool. Given what is known about salivary substance P and salivary gland innervation, we hypothesized that salivary substance P and/or beta-endorphin might reflect the basal activity of these neuropeptides in the central nervous system, thereby perhaps mirroring a general propensity to chronic pain. Based on this overall hypothesis, our aim was to compare salivary levels of these neuropeptides in chronic neuropathic pain patients with healthy controls. An additional aim was to relate salivary levels to plasma levels. MATERIALS AND METHODS: We compared salivary concentrations of beta-endorphin and substance P in 14 chronic neuropathic pain patients with concentrations in 18 healthy controls using a Luminex technology kit. Salivary-to-plasma quotients were also calculated. RESULTS: We found no significant difference between the groups' salivary concentrations of substance P and beta-endorphin. No correlation was found between salivary and plasma concentrations of each neuropeptide, which we hypothesize might point to local production of beta-endorphin and/or substance P in the salivary glands. Given high substance P salivary-to-plasma quotients, such a local production seems more likely for substance P than for beta-endorphin. CONCLUSIONS: Propensity to neuropathic chronic pain was not substantiated by our analysis of salivary levels of substance P and/or beta-endorphin. However, we report salivary-to-plasma quotients that give potentially important physiological insight about these neuropeptides.

Transcriptional and Epigenetic Changes Influencing Skeletal Muscle Metabolism in Women With Polycystic Ovary Syndrome.

Context: Insulin resistance in skeletal muscle is a major risk factor for the development of type 2 diabetes in women with polycystic ovary syndrome (PCOS). Despite this, the mechanisms underlying insulin resistance in PCOS are largely unknown. Objective: To investigate the genome-wide DNA methylation and gene expression patterns in skeletal muscle from women with PCOS and controls and relate them to phenotypic variations. Design/Participants: In a case-control study, skeletal muscle biopsies from women with PCOS (n = 17) and age-, weight-, and body mass index‒matched controls (n = 14) were analyzed by array-based DNA methylation and mRNA expression profiling. Results: Eighty-five unique transcripts were differentially expressed in muscle from women with PCOS vs controls, including DYRK1A, SYNPO2, SCP2, and NAMPT. Furthermore, women with PCOS had reduced expression of genes involved in immune system pathways. Two CpG sites showed differential DNA methylation after correction for multiple testing. However, an mRNA expression of ∼30% of the differentially expressed genes correlated with DNA methylation levels of CpG sites in or near the gene. Functional follow-up studies demonstrated that KLF10 is under transcriptional control of insulin, where insulin promotes glycogen accumulation in myotubes of human muscle cells. Testosterone downregulates the expression levels of COL1A1 and MAP2K6. Conclusion: PCOS is associated with aberrant skeletal muscle gene expression with dysregulated pathways. Furthermore, we identified specific changes in muscle DNA methylation that may affect gene expression. This study showed that women with PCOS have epigenetic and transcriptional changes in skeletal muscle that, in part, can explain the metabolic abnormalities seen in these women.

SWI/SNF regulates half of its targets without the need of ATP-driven nucleosome remodeling by Brahma.

BACKGROUND: Brahma (BRM) is the only catalytic subunit of the SWI/SNF chromatin-remodeling complex of Drosophila melanogaster. The function of SWI/SNF in transcription has long been attributed to its ability to remodel nucleosomes, which requires the ATPase activity of BRM. However, recent studies have provided evidence for a non-catalytic function of BRM in the transcriptional regulation of a few specific genes. RESULTS: Here we have used RNA-seq and ChIP-seq to identify the BRM target genes in S2 cells, and we have used a catalytically inactive BRM mutant (K804R) that is unable to hydrolyze ATP to investigate the magnitude of the non-catalytic function of BRM in transcription regulation. We show that 49% of the BRM target genes in S2 cells are regulated through mechanisms that do not require BRM to have an ATPase activity. We also show that the catalytic and non-catalytic mechanisms of SWI/SNF regulation operate on two subsets of genes that differ in promoter architecture and are linked to different biological processes. CONCLUSIONS: This study shows that the non-catalytic role of SWI/SNF in transcription regulation is far more prevalent than previously anticipated and that the genes that are regulated by SWI/SNF through ATPase-dependent and ATPase-independent mechanisms have specialized roles in different cellular and developmental processes.

Evaluation of the ISL1 gene in the pathogenesis of bladder exstrophy in a Swedish cohort.

Bladder exstrophy is a congenital closure defect of the urinary bladder with a profound effect on morbidity. Although the malformation is usually sporadic, a genetic background is supported by an increased recurrence risk in relatives, higher concordance rates in monozygotic twins and several associated chromosomal aberrations. Recently, the ISL1 gene was presented as a candidate gene for bladder exstrophy and epispadias complex (BEEC) development in two different studies. In our study, we screened for genetic variants in the ISL1 gene in DNA from 125 Swedish patients using Sanger sequencing and array-CGH analysis. In addition, we evaluated ISL1 expression in RNA of human bladder during embryonic and fetal weeks 5-10 relative to that in lung tissue (week 9). In total, 21 single-nucleotide variants were identified, including a potentially novel missense variant, c.137C>G p.(Ala46Gly), substituting a conserved amino acid. This variant was inherited from an unaffected mother. No structural variants were identified. RNA sequencing revealed ISL1 mRNA expression during the critical time frame of human bladder development. In conclusion, we did not detect any known or likely pathogenic variants in the ISL1 gene in 125 Swedish BEEC patients, indicating that variation in the ISL1 gene is not a common genetic mechanism of BEEC development in the Swedish population.

Gene expression profiling across ontogenetic stages in the wood white (Leptidea sinapis) reveals pathways linked to butterfly diapause regulation.

In temperate latitudes, many insects enter diapause (dormancy) during the cold season, a period during which developmental processes come to a standstill. The wood white (Leptidea sinapis) is a butterfly species distributed across western Eurasia that shows photoperiod-induced diapause with variation in critical day-length across populations at different latitudes. We assembled transcriptomes and estimated gene expression levels at different developmental stages in experimentally induced directly developing and diapausing cohorts of a single Swedish population of L. sinapis to investigate the regulatory mechanisms underpinning diapause initiation. Different day lengths resulted in expression changes of developmental genes and affected the rate of accumulation of signal molecules, suggesting that diapause induction might be controlled by increased activity of monoamine neurotransmitters in larvae reared under short-day light conditions. Expression differences between light treatment groups of two monoamine regulator genes (DDC and ST) were observed already in instar III larvae. Once developmental pathways were irreversibly set at instar V, a handful of genes related to dopamine production were differentially expressed leading to a significant decrease in expression of global metabolic genes and increase in expression of genes related to fatty acid synthesis and sequestration. This is in line with a time-dependent (hour-glass) model of diapause regulation where a gradual shift in the concentration of monoamine neurotransmitters and their metabolites during development of larvae under short-day conditions leads to increased storage of fat, decreased energy expenditures, and ultimately developmental stasis at the pupal stage.

A Single Bout of Electroacupuncture Remodels Epigenetic and Transcriptional Changes in Adipose Tissue in Polycystic Ovary Syndrome.

A single bout of electroacupuncture results in muscle contractions and increased whole body glucose uptake in women with polycystic ovary syndrome (PCOS). Women with PCOS have transcriptional and epigenetic alterations in the adipose tissue and we hypothesized that electroacupuncture induces epigenetic and transcriptional changes to restore metabolic alterations. Twenty-one women with PCOS received a single bout of electroacupuncture, which increased the whole body glucose uptake. In subcutaneous adipose tissue biopsies, we identified treatment-induced expression changes of 2369 genes (Q < 0.05) and DNA methylation changes of 7055 individual genes (Q = 0.11). The largest increase in expression was observed for FOSB (2405%), and the largest decrease for LOC100128899 (54%). The most enriched pathways included Acute phase response signaling and LXR/RXR activation. The DNA methylation changes ranged from 1-16%, and 407 methylation sites correlated with gene expression. Among genes known to be differentially expressed in PCOS, electroacupuncture reversed the expression of 80 genes, including PPARγ and ADIPOR2. Changes in the expression of Nr4a2 and Junb are reversed by adrenergic blockers in rats demonstrating that changes in gene expression, in part, is due to activation of the sympathetic nervous system. In conclusion, low-frequency electroacupuncture with muscle contractions remodels epigenetic and transcriptional changes that elicit metabolic improvement.

RNA-sequence data normalization through in silico prediction of reference genes: the bacterial response to DNA damage as case study.

BACKGROUND: Measuring how gene expression changes in the course of an experiment assesses how an organism responds on a molecular level. Sequencing of RNA molecules, and their subsequent quantification, aims to assess global gene expression changes on the RNA level (transcriptome). While advances in high-throughput RNA-sequencing (RNA-seq) technologies allow for inexpensive data generation, accurate post-processing and normalization across samples is required to eliminate any systematic noise introduced by the biochemical and/or technical processes. Existing methods thus either normalize on selected known reference genes that are invariant in expression across the experiment, assume that the majority of genes are invariant, or that the effects of up- and down-regulated genes cancel each other out during the normalization. RESULTS: Here, we present a novel method, moose2 , which predicts invariant genes in silico through a dynamic programming (DP) scheme and applies a quadratic normalization based on this subset. The method allows for specifying a set of known or experimentally validated invariant genes, which guides the DP. We experimentally verified the predictions of this method in the bacterium Escherichia coli, and show how moose2 is able to (i) estimate the expression value distances between RNA-seq samples, (ii) reduce the variation of expression values across all samples, and (iii) to subsequently reveal new functional groups of genes during the late stages of DNA damage. We further applied the method to three eukaryotic data sets, on which its performance compares favourably to other methods. The software is implemented in C++ and is publicly available from http://grabherr.github.io/moose2/. CONCLUSIONS: The proposed RNA-seq normalization method, moose2 , is a valuable alternative to existing methods, with two major advantages: (i) in silico prediction of invariant genes provides a list of potential reference genes for downstream analyses, and (ii) non-linear artefacts in RNA-seq data are handled adequately to minimize variations between replicates.

Cerium oxide nanoparticles inhibit differentiation of neural stem cells.

Cerium oxide nanoparticles (nanoceria) display antioxidant properties and have shown cytoprotective effects both in vitro and in vivo. Here, we explored the effects of nanoceria on neural progenitor cells using the C17.2 murine cell line as a model. First, we assessed the effects of nanoceria versus samarium (Sm) doped nanoceria on cell viability in the presence of the prooxidant, DMNQ. Both particles were taken up by cells and nanoceria, but not Sm-doped nanoceria, elicited a temporary cytoprotective effect upon exposure to DMNQ. Next, we employed RNA sequencing to explore the transcriptional responses induced by nanoceria or Sm-doped nanoceria during neuronal differentiation. Detailed computational analyses showed that nanoceria altered pathways and networks relevant for neuronal development, leading us to hypothesize that nanoceria inhibits neuronal differentiation, and that nanoceria and Sm-doped nanoceria both interfere with cytoskeletal organization. We confirmed that nanoceria reduced neuron specific ß3-tubulin expression, a marker of neuronal differentiation, and GFAP, a neuroglial marker. Furthermore, using super-resolution microscopy approaches, we could show that both particles interfered with cytoskeletal organization and altered the structure of neural growth cones. Taken together, these results reveal that nanoceria may impact on neuronal differentiation, suggesting that nanoceria could pose a developmental neurotoxicity hazard.

Persistent Effects of Developmental Exposure to 17α-Ethinylestradiol on the Zebrafish (Danio rerio) Brain Transcriptome and Behavior.

The synthetic estrogen 17α-ethinylestradiol (EE2) is an endocrine disrupting compound of concern due to its persistence and widespread presence in the aquatic environment. Effects of developmental exposure to low concentrations of EE2 in fish on reproduction and behavior not only persisted to adulthood, but have also been observed to be transmitted to several generations of unexposed progeny. To investigate the possible biological mechanisms of the persistent anxiogenic phenotype, we exposed zebrafish embryos for 80 days post fertilization to 0, 3, and 10 ng/L EE2 (measured concentrations 2.14 and 7.34 ng/L). After discontinued exposure, the animals were allowed to recover for 120 days in clean water. Adult males and females were later tested for changes in stress response and shoal cohesion, and whole-brain gene expression was analyzed with RNA sequencing. The results show increased anxiety in the novel tank and scototaxis tests, and increased shoal cohesion in fish exposed during development to EE2. RNA sequencing revealed 34 coding genes differentially expressed in male brains and 62 in female brains as a result of EE2 exposure. Several differences were observed between males and females in differential gene expression, with only one gene, sv2b, coding for a synaptic vesicle protein, that was affected by EE2 in both sexes. Functional analyses showed that in female brains, EE2 had significant effects on pathways connected to the circadian rhythm, cytoskeleton and motor proteins and synaptic proteins. A large number of non-coding sequences including 19 novel miRNAs were also differentially expressed in the female brain. The largest treatment effect in male brains was observed in pathways related to cholesterol biosynthesis and synaptic proteins. Circadian rhythm and cholesterol biosynthesis, previously implicated in anxiety behavior, might represent possible candidate pathways connecting the transcriptome changes to the alterations to behavior. Further the observed alteration in expression of genes involved in synaptogenesis and synaptic function may be important for the developmental modulations resulting in an anxiety phenotype. This study represents an initial survey of the fish brain transcriptome by RNA sequencing after long-term recovery from developmental exposure to an estrogenic compound.

Local Adaptation in European Firs Assessed through Extensive Sampling across Altitudinal Gradients in Southern Europe.

BACKGROUND: Local adaptation is a key driver of phenotypic and genetic divergence at loci responsible for adaptive traits variations in forest tree populations. Its experimental assessment requires rigorous sampling strategies such as those involving population pairs replicated across broad spatial scales. METHODS: A hierarchical Bayesian model of selection (HBM) that explicitly considers both the replication of the environmental contrast and the hierarchical genetic structure among replicated study sites is introduced. Its power was assessed through simulations and compared to classical 'within-site' approaches (FDIST, BAYESCAN) and a simplified, within-site, version of the model introduced here (SBM). RESULTS: HBM demonstrates that hierarchical approaches are very powerful to detect replicated patterns of adaptive divergence with low false-discovery (FDR) and false-non-discovery (FNR) rates compared to the analysis of different sites separately through within-site approaches. The hypothesis of local adaptation to altitude was further addressed by analyzing replicated Abies alba population pairs (low and high elevations) across the species' southern distribution range, where the effects of climatic selection are expected to be the strongest. For comparison, a single population pair from the closely related species A. cephalonica was also analyzed. The hierarchical model did not detect any pattern of adaptive divergence to altitude replicated in the different study sites. Instead, idiosyncratic patterns of local adaptation among sites were detected by within-site approaches. CONCLUSION: Hierarchical approaches may miss idiosyncratic patterns of adaptation among sites, and we strongly recommend the use of both hierarchical (multi-site) and classical (within-site) approaches when addressing the question of adaptation across broad spatial scales.

Different Alleles of a Gene Encoding Leucoanthocyanidin Reductase (PaLAR3) Influence Resistance against the Fungus Heterobasidion parviporum in Picea abies.

Despite the fact that fungal diseases are a growing menace for conifers in modern silviculture, only a very limited number of molecular markers for pathogen resistance have been validated in conifer species. A previous genetic study indicated that the resistance of Norway spruce (Picea abies) to Heterobasidion annosum s.l., a pathogenic basidiomycete species complex, is linked to a quantitative trait loci that associates with differences in fungal growth in sapwood (FGS) that includes a gene, PaLAR3, which encodes a leucoanthocyanidin reductase. In this study, gene sequences showed the presence of two PaLAR3 allelic lineages in P. abies. Higher resistance was associated with the novel allele, which was found in low frequency in the four P. abies populations that we studied. Norway spruce plants carrying at least one copy of the novel allele showed a significant reduction in FGS after inoculation with Heterobasidion parviporum compared to their half-siblings carrying no copies, indicating dominance of this allele. The amount of (+) catechin, the enzymatic product of PaLAR3, was significantly higher in bark of trees homozygous for the novel allele. Although we observed that the in vitro activities of the enzymes encoded by the two alleles were similar, we could show that allele-specific transcript levels were significantly higher for the novel allele, indicating that regulation of gene expression is responsible for the observed effects in resistance, possibly caused by differences in cis-acting elements that we observe in the promoter region of the two alleles.

Identifying Genetic Signatures of Natural Selection Using Pooled Population Sequencing in Picea abies.

The joint inference of selection and past demography remain a costly and demanding task. We used next generation sequencing of two pools of 48 Norway spruce mother trees, one corresponding to the Fennoscandian domain, and the other to the Alpine domain, to assess nucleotide polymorphism at 88 nuclear genes. These genes are candidate genes for phenological traits, and most belong to the photoperiod pathway. Estimates of population genetic summary statistics from the pooled data are similar to previous estimates, suggesting that pooled sequencing is reliable. The nonsynonymous SNPs tended to have both lower frequency differences and lower FST values between the two domains than silent ones. These results suggest the presence of purifying selection. The divergence between the two domains based on synonymous changes was around 5 million yr, a time similar to a recent phylogenetic estimate of 6 million yr, but much larger than earlier estimates based on isozymes. Two approaches, one of them novel and that considers both FST and difference in allele frequencies between the two domains, were used to identify SNPs potentially under diversifying selection. SNPs from around 20 genes were detected, including genes previously identified as main target for selection, such as PaPRR3 and PaGI.

The extent and meaning of hybridization and introgression between Siberian spruce (Picea obovata) and Norway spruce (Picea abies): cryptic refugia as stepping stones to the west?

Boreal species were repeatedly exposed to ice ages and went through cycles of contraction and expansion while sister species alternated periods of contact and isolation. The resulting genetic structure is consequently complex, and demographic inferences are intrinsically challenging. The range of Norway spruce (Picea abies) and Siberian spruce (Picea obovata) covers most of northern Eurasia; yet their geographical limits and histories remain poorly understood. To delineate the hybrid zone between the two species and reconstruct their joint demographic history, we analysed variation at nuclear SSR and mitochondrial DNA in 102 and 88 populations, respectively. The dynamics of the hybrid zone was analysed with approximate Bayesian computation (ABC) followed by posterior predictive structure plot reconstruction and the presence of barriers across the range tested with estimated effective migration surfaces. To estimate the divergence time between the two species, nuclear sequences from two well-separated populations of each species were analysed with ABC. Two main barriers divide the range of the two species: one corresponds to the hybrid zone between them, and the other separates the southern and northern domains of Norway spruce. The hybrid zone is centred on the Urals, but the genetic impact of Siberian spruce extends further west. The joint distribution of mitochondrial and nuclear variation indicates an introgression of mitochondrial DNA from Norway spruce into Siberian spruce. Overall, our data reveal a demographic history where the two species interacted frequently and where migrants originating from the Urals and the West Siberian Plain recolonized northern Russia and Scandinavia using scattered refugial populations of Norway spruce as stepping stones towards the west.

A Picea abies linkage map based on SNP markers identifies QTLs for four aspects of resistance to Heterobasidion parviporum infection.

A consensus linkage map of Picea abies, an economically important conifer, was constructed based on the segregation of 686 SNP markers in a F1 progeny population consisting of 247 individuals. The total length of 1889.2 cM covered 96.5% of the estimated genome length and comprised 12 large linkage groups, corresponding to the number of haploid P. abies chromosomes. The sizes of the groups (from 5.9 to 9.9% of the total map length) correlated well with previous estimates of chromosome sizes (from 5.8 to 10.8% of total genome size). Any locus in the genome has a 97% probability to be within 10 cM from a mapped marker, which makes the map suited for QTL mapping. Infecting the progeny trees with the root rot pathogen Heterobasidion parviporum allowed for mapping of four different resistance traits: lesion length at the inoculation site, fungal spread within the sapwood, exclusion of the pathogen from the host after initial infection, and ability to prevent the infection from establishing at all. These four traits were associated with two, four, four and three QTL regions respectively of which none overlapped between the traits. Each QTL explained between 4.6 and 10.1% of the respective traits phenotypic variation. Although the QTL regions contain many more genes than the ones represented by the SNP markers, at least four markers within the confidence intervals originated from genes with known function in conifer defence; a leucoanthocyanidine reductase, which has previously been shown to upregulate during H. parviporum infection, and three intermediates of the lignification process; a hydroxycinnamoyl CoA shikimate/quinate hydroxycinnamoyltransferase, a 4-coumarate CoA ligase, and a R2R3-MYB transcription factor.

Patterns of nucleotide diversity at photoperiod related genes in Norway spruce [Picea abies (L.) Karst].

The ability of plants to track seasonal changes is largely dependent on genes assigned to the photoperiod pathway, and variation in those genes is thereby important for adaptation to local day length conditions. Extensive physiological data in several temperate conifer species suggest that populations are adapted to local light conditions, but data on the genes underlying this adaptation are more limited. Here we present nucleotide diversity data from 19 genes putatively involved in photoperiodic response in Norway spruce (Picea abies). Based on similarity to model plants the genes were grouped into three categories according to their presumed position in the photoperiod pathway: photoreceptors, circadian clock genes, and downstream targets. An HKA (Hudson, Kreitman and Aquade) test showed a significant excess of diversity at photoreceptor genes, but no departure from neutrality at circadian genes and downstream targets. Departures from neutrality were also tested with Tajima's D and Fay and Wu's H statistics under three demographic scenarios: the standard neutral model, a population expansion model, and a more complex population split model. Only one gene, the circadian clock gene PaPRR3 with a highly positive Tajima's D value, deviates significantly from all tested demographic scenarios. As the PaPRR3 gene harbours multiple non-synonymous variants it appears as an excellent candidate gene for control of photoperiod response in Norway spruce.

Clinal variation at phenology-related genes in spruce: parallel evolution in FTL2 and Gigantea?

Parallel clines in different species, or in different geographical regions of the same species, are an important source of information on the genetic basis of local adaptation. We recently detected latitudinal clines in SNPs frequencies and gene expression of candidate genes for growth cessation in Scandinavian populations of Norway spruce (Picea abies). Here we test whether the same clines are also present in Siberian spruce (P. obovata), a close relative of Norway spruce with a different Quaternary history. We sequenced nine candidate genes and 27 control loci and genotyped 14 SSR loci in six populations of P. obovata located along the Yenisei river from latitude 56°N to latitude 67°N. In contrast to Scandinavian Norway spruce that both departs from the standard neutral model (SNM) and shows a clear population structure, Siberian spruce populations along the Yenisei do not depart from the SNM and are genetically unstructured. Nonetheless, as in Norway spruce, growth cessation is significantly clinal. Polymorphisms in photoperiodic (FTL2) and circadian clock (Gigantea, GI, PRR3) genes also show significant clinal variation and/or evidence of local selection. In GI, one of the variants is the same as in Norway spruce. Finally, a strong cline in gene expression is observed for FTL2, but not for GI. These results, together with recent physiological studies, confirm the key role played by FTL2 and circadian clock genes in the control of growth cessation in spruce species and suggest the presence of parallel adaptation in these two species.

No time for spruce: rapid dampening of circadian rhythms in Picea abies (L. Karst).

The identification and cloning of full-length homologs of circadian clock genes from Picea abies represent a first step to study the function and evolution of the circadian clock in gymnosperms. Phylogenetic analyses suggest that the sequences of key circadian clock genes are conserved between angiosperms and gymnosperms. though fewer homologous copies were found for most gene families in P. abies. We detected diurnal cycling of circadian clock genes in P. abies using quantitative real-time PCR; however, cycling appeared to be rapidly dampened under free-running conditions. Given the unexpected absence of transcriptional cycling during constant conditions, we employed a complementary method to assay circadian rhythmic outputs and measured delayed fluorescence in seedlings of Norway spruce. Neither of the two approaches to study circadian rhythms in Norway spruce could detect robust ∼24 h cycling behavior under constant conditions. These data suggest gene conservation but fundamental differences in clock function between gymnosperms and other plant taxa.

Conserved function of core clock proteins in the gymnosperm Norway spruce (Picea abies L. Karst).

From studies of the circadian clock in the plant model species Arabidopsis (Arabidopsis thaliana), a number of important properties and components have emerged. These include the genes CIRCADIAN CLOCK ASSOCIATED 1 (CCA1), GIGANTEA (GI), ZEITLUPE (ZTL) and TIMING OF CAB EXPRESSION 1 (TOC1 also known as PSEUDO-RESPONSE REGULATOR 1 (PRR1)) that via gene expression feedback loops participate in the circadian clock. Here, we present results from ectopic expression of four Norway spruce (Picea abies) putative homologs (PaCCA1, PaGI, PaZTL and PaPRR1) in Arabidopsis, their flowering time, circadian period length, red light response phenotypes and their effect on endogenous clock genes were assessed. For PaCCA1-ox and PaZTL-ox the results were consistent with Arabidopsis lines overexpressing the corresponding Arabidopsis genes. For PaGI consistent results were obtained when expressed in the gi2 mutant, while PaGI and PaPRR1 expressed in wild type did not display the expected phenotypes. These results suggest that protein function of PaCCA1, PaGI and PaZTL are at least partly conserved compared to Arabidopsis homologs, however further studies are needed to reveal the protein function of PaPRR1. Our data suggest that components of the three-loop network typical of the circadian clock in angiosperms were present before the split of gymnosperms and angiosperms.

A significant fraction of 21-nucleotide small RNA originates from phased degradation of resistance genes in several perennial species.

Small RNAs (sRNAs), including microRNA (miRNA) and short-interfering RNA (siRNA), are important in the regulation of diverse biological processes. Comparative studies of sRNAs from plants have mainly focused on miRNA, even though they constitute a mere fraction of the total sRNA diversity. In this study, we report results from an in-depth analysis of the sRNA population from the conifer spruce (Picea abies) and compared the results with those of a range of plant species. The vast majority of sRNA sequences in spruce can be assigned to 21-nucleotide-long siRNA sequences, of which a large fraction originate from the degradation of transcribed sequences related to nucleotide-binding site-leucine-rich repeat-type resistance genes. Over 90% of all genes predicted to contain either a Toll/interleukin-1 receptor or nucleotide-binding site domain showed evidence of siRNA degradation. The data further suggest that this phased degradation of resistance-related genes is initiated from miRNA-guided cleavage, often by an abundant 22-nucleotide miRNA. Comparative analysis over a range of plant species revealed a huge variation in the abundance of this phenomenon. The process seemed to be virtually absent in several species, including Arabidopsis (Arabidopsis thaliana), rice (Oryza sativa), and nonvascular plants, while particularly high frequencies were observed in spruce, grape (Vitis vinifera), and poplar (Populus trichocarpa). This divergent pattern might reflect a mechanism to limit runaway transcription of these genes in species with rapidly expanding nucleotide-binding site-leucine-rich repeat gene families. Alternatively, it might reflect variation in a counter-counter defense mechanism between plant species.