Umeå University's logo

Long non-coding RNAs (lncRNAs) have emerged as important regulators of many bio- logical processes, although their regulatory roles remain poorly characterized in woody plants, especially in gymnosperms. A major challenge of working with lncRNAs is to assign functional annotations, since they have a low coding potential and low cross-species conservation.

We utilised an existing RNA-Sequencing resource and performed short RNA sequencing of somatic embryogenesis developmental stages in Norway spruce (Picea abies L. Karst). We implemented a pipeline to identify lncRNAs located within the intergenic space (lincRNAs) and generated a co-expression network including protein coding, lincRNA and miRNA genes.

To assign putative functional annotation, we employed a guilt-by-association approach using the co-expression network and integrated these results with annota- tion assigned using semantic similarity and co-expression. Moreover, we evaluated the relationship between lincRNAs and miRNAs, and identified which lincRNAs are conserved in other species. We identified lincRNAs with clear evidence of differential expression during somatic embryogenesis and used network connectivity to identify those with the greatest regulatory potential.

This work provides the most comprehensive view of lincRNAs in Norway spruce and is the first study to perform global identification of lincRNAs during somatic embryogen- esis in conifers. The data have been integrated into the expression visualisation tools at the PlantGenIE.org web resource to enable easy access to the community. This will facilitate the use of the data to address novel questions about the role of lincRNAs in the regulation of embryogenesis and facilitate future comparative genomics studies.

Aspen (Populus tremula L.) is a keystone species and a model system for forest tree genomics. We present an updated resource comprising a chromosome-scale assem- bly, population genetics and genomics data. Using the resource, we explore the genetic basis of natural variation in leaf size and shape, traits with complex genetic architecture.

We generated the genome assembly using long-read sequencing, optical and high-density genetic maps. We conducted whole-genome resequencing of the Umeå Aspen (UmAsp) collection. Using the assembly and re-sequencing data from the UmAsp, Swedish Aspen (SwAsp) and Scottish Aspen (ScotAsp) collections we performed genome-wide association analyses (GWAS) using Single Nucleotide Polymorphisms (SNPs) for 26 leaf physiognomy phenotypes. We conducted Assay of Transposase Accessible Chromatin sequencing (ATAC-Seq), identified genomic regions of accessible chromatin, and subset SNPs to these regions, improving the GWAS detection rate. We identified candidate long non-coding RNAs in leaf samples, quantified their expression in an updated co-expression network, and used this to explore the functions of candidate genes identified from the GWAS.

A GWAS found SNP associations for seven traits. The associated SNPs were in or near genes annotated with developmental functions, which represent candidates for further study. Of particular interest was a !177-kbp region harbouring associations with several leaf phenotypes in ScotAsp.

We have incorporated the assembly, population genetics, genomics, and GWAS data into the PlantGenIE.org web resource, including updating existing genomics data to the new genome version, to enable easy exploration and visualisation. We provide all raw and processed data to facilitate reuse in future studies.

Leaf senescence can be induced by stress or aging, sometimes in a synergistic manner. It is generally acknowledged that the ability to withstand senescence-inducing conditions can provide plants with stress resilience. Although the signaling and transcriptional networks responsible for a delayed senescence phenotype, often referred to as a functional stay-green trait, have been actively investigated, very little is known about the subsequent metabolic adjustments conferring this aptitude to survival. First, using the individually darkened leaf (IDL) experimental setup, we compared IDLs of wild-type (WT) Arabidopsis (Arabidopsis thaliana) to several stay-green contexts, that is IDLs of two functional stay-green mutant lines, oresara1-2 (ore1-2) and an allele of phytochrome-interacting factor 5 (pif5), as well as to leaves from a WT plant entirely darkened (DP). We provide compelling evidence that arginine and ornithine, which accumulate in all stay-green contexts—likely due to the lack of induction of amino acids (AAs) transport—can delay the progression of senescence by fueling the Krebs cycle or the production of polyamines (PAs). Secondly, we show that the conversion of putrescine to spermidine (SPD) is controlled in an age-dependent manner. Thirdly, we demonstrate that SPD represses senescence via interference with ethylene signaling by stabilizing the ETHYLENE BINDING FACTOR1 and 2 (EBF1/2) complex. Taken together, our results identify arginine and ornithine as central metabolites influencing the stress- and age-dependent progression of leaf senescence. We propose that the regulatory loop between the pace of the AA export and the progression of leaf senescence provides the plant with a mechanism to fine-tune the induction of cell death in leaves, which, if triggered unnecessarily, can impede nutrient remobilization and thus plant growth and survival.

Boreal conifers possess a tremendous ability to survive and remain evergreen during harsh winter conditions and resume growth during summer. This is enabled by coordinated regulation of major cellular functions at the level of gene expression, metabolism, and physiology. Here we present a comprehensive characterization of the annual changes in the global transcriptome of Norway spruce (Picea abies) needles as a resource to understand needle development and acclimation processes throughout the year. In young, growing needles (May 15 until June 30), cell walls, organelles, etc., were formed, and this developmental program heavily influenced the transcriptome, explained by over-represented Gene Ontology (GO) categories. Later changes in gene expression were smaller but four phases were recognized: summer (July–August), autumn (September–October), winter (November–February), and spring (March–April), where over-represented GO categories demonstrated how the needles acclimated to the various seasons. Changes in the seasonal global transcriptome profile were accompanied by differential expression of members of the major transcription factor families. We present a tentative model of how cellular activities are regulated over the year in needles of Norway spruce, which demonstrates the value of mining this dataset, accessible in ConGenIE together with advanced visualization tools.

Summary: Since its introduction, RNA-Seq technology has been used extensively in studies of pathogenic bacteria to identify and quantify differences in gene expression across multiple samples from bacteria exposed to different conditions. With some exceptions, tools for studying gene expression, determination of differential gene expression, downstream pathway analysis and normalization of data collected in extreme biological conditions is still lacking. Here, we describe ProkSeq, a user-friendly, fully automated RNA-Seq data analysis pipeline designed for prokaryotes. ProkSeq provides a wide variety of options for analysing differential expression, normalizing expression data and visualizing data and results.

Availability and implementation: ProkSeq is implemented in Python and is published under the MIT source license. The pipeline is available as a Docker container https://hub.docker.com/repository/docker/snandids/prokseq-v2.0, or can be used through Anaconda: https://anaconda.org/snandiDS/prokseq. The code is available on Github: https://github.com/snandiDS/prokseq and a detailed user documentation, including a manual and tutorial can be found at https://prokseqV20.readthedocs.io.

Background: Translation is a tightly regulated process, controlling the rate of protein synthesis in cells. Ribosome sequencing (Ribo-Seq) is a recently developed tool for studying actively translated mRNA and can thus directly address translational regulation. Ribo-Seq libraries need to be sequenced to a great depth due to high contamination by rRNA and other contaminating nucleic acid fragments. Deep sequencing is expensive, and it generates large volumes of data, making data analysis complicated and time consuming.

Methods and results: Here we developed a platform for Ribo-Seq library construction and data analysis to enable rapid quality assessment of Ribo-Seq libraries with the help of a small-scale sequencer. Our data show that several qualitative features of a Ribo-Seq library, such as read length distribution, P-site distribution, reading frame and triplet periodicity, can be effectively evaluated using only the data generated by a benchtop sequencer with a very limited number of reads.

Conclusion: Our pipeline enables rapid evaluation of Ribo-Seq libraries, opening up possibilities for optimization of Ribo-Seq library construction from difficult samples, and leading to better decision making prior to more costly deep sequencing.

PIRIN2 (PRN2) was earlier reported to suppress syringyl (S)-type lignin accumulation of xylem vessels of Arabidopsis thaliana. In the present study, we report yeast two-hybrid results supporting the interaction of PRN2 with HISTONE MONOUBIQUITINATION2 (HUB2) in Arabidopsis. HUB2 has been previously implicated in several plant developmental processes, but not in lignification. Interaction between PRN2 and HUB2 was verified by β-galactosidase enzymatic and co-immunoprecipitation assays. HUB2 promoted the deposition of S-type lignin in the secondary cell walls of both stem and hypocotyl tissues, as analysed by pyrolysis-GC/MS. Chemical fingerprinting of individual xylem vessel cell walls by Raman and Fourier transform infrared microspectroscopy supported the function of HUB2 in lignin deposition. These results, together with a genetic analysis of the hub2 prn2 double mutant, support the antagonistic function of PRN2 and HUB2 in deposition of S-type lignin. Transcriptome analyses indicated the opposite regulation of the S-type lignin biosynthetic gene FERULATE-5-HYDROXYLASE1 by PRN2 and HUB2 as the underlying mechanism. PRN2 and HUB2 promoter activities co-localized in cells neighbouring the xylem vessel elements, suggesting that the S-type lignin-promoting function of HUB2 is antagonized by PRN2 for the benefit of the guaiacyl (G)-type lignin enrichment of the neighbouring xylem vessel elements.