Metabolomics is an emerging ‘omics’ science involving the characterization of metabolites and metabolism in biological systems. Few bioinformatic tools have been developed for the visualization, exploration and analysis of metabolomic data within the context of metabolic pathways: some of them became rapidly obsolete and are no longer supported, others are based on a single database. A systematic collection of existing annotations has the potential of considerably boosting the investigation and contextualization of metabolomic measurements. We have released a major update of our Bioconductor package graphite which explicitly tracks small molecules within pathway topologies and their interactions with proteins. The package gathers the information stored in eight major databases, oriented both at genes and at metabolites, across 14 different species. Depending on user preferences, all pathways can be retrieved as gene-only, genemetabolite or metabolite-only networks. The new graphite version (1.24) is available on Bioconductor. Supplementary data are available at Bioinformatics online.

It is still continuously debated whether the low-dose/dose-rate (LDR) of ionizing radiation represents a hazard for humans. Model organisms, such as fruit flies, are considered valuable systems to reveal insights into this issue. We found that, in wild-type Drosophila melanogaster larval neuroblasts, the frequency of Chromosome Breaks (CBs), induced by acute γ-irradiation, is considerably reduced when flies are previously exposed to a protracted dose of 0.4 Gy delivered at a dose rate of 2.5 mGy/h. This indicates that this exposure, which is associated with an increased expression of DNA damage response proteins, induces a radioadaptive response (RAR) that protects Drosophila from extensive DNA damage. Interestingly, the same exposure reduces the frequency of telomere fusions (TFs) from Drosophila telomere capping mutants suggesting that the LDR can generally promote a protective response on chromatin sites that are recognized as DNA breaks. Deep RNA sequencing revealed that RAR is associated with a reduced expression of Loquacious D (Loqs-RD) gene that encodes a well-conserved dsRNA binding protein required for esiRNAs biogenesis. Remarkably, loss of Loqs mimics the LDR-mediated chromosome protection as it decreases the IR-induced CBs and TFs frequency. Thus, our molecular characterization of RAR identifies Loqs as a key factor in the cellular response to LDR and in the epigenetic routes involved in radioresistance.

Time-course gene expression experiments are useful tools for exploring biological processes. In this type of experiments, gene expression changes are monitored along time. Unfortunately, replication of time series is still costly and usually long time course do not have replicates. Many approaches have been proposed to deal with this data structure, but none of them in the field of pathway analysis. Pathway analyses have acquired great relevance for helping the interpretation of gene expression data. Several methods have been proposed to this aim: from the classical enrichment to the more complex topological analysis that gains power from the topology of the pathway. None of them were devised to identify temporal variations in time course data.

The krill species Euphausia superba plays a critical role in the food chain of the Antarctic ecosystem. Significant changes in climate conditions observed in the Antarctic Peninsula region in the last decades have already altered the distribution of krill and its reproductive dynamics. A deeper understanding of the adaptation capabilities of this species is urgently needed. The availability of a large body of RNA-seq assays allowed us to extend the current knowledge of the krill transcriptome. Our study covered the entire developmental process providing information of central relevance for ecological studies. Here we identified a series of genes involved in different steps of the krill moulting cycle, in the reproductive process and in sexual maturation in accordance with what was already described in previous works. Furthermore, the new transcriptome highlighted the presence of differentially expressed genes previously unknown, playing important roles in cuticle development as well as in energy storage during the krill life cycle. The discovery of new opsin sequences, specifically rhabdomeric opsins, one onychopsin, and one non-visual arthropsin, expands our knowledge of the krill opsin repertoire. We have collected all these results into the KrillDB² database, a resource combining the latest annotation of the krill transcriptome with a series of analyses targeting genes relevant to krill physiology. KrillDB² provides in a single resource a comprehensive catalog of krill genes; an atlas of their expression profiles over all RNA-seq datasets publicly available; a study of differential expression across multiple conditions. Finally, it provides initial indications about the expression of microRNA precursors, whose contribution to krill physiology has never been reported before.

Genomic imprinting and X chromosome inactivation (XCI) are two prototypical epigenetic mechanisms whereby a set of genes is expressed mono-allelically in order to fine-tune their expression levels. Defects in genomic imprinting have been observed in several neurodevelopmental disorders, in a wide range of tumours and in induced pluripotent stem cells (iPSCs). Single Nucleotide Variants (SNVs) are readily detectable by RNA-sequencing allowing the determination of whether imprinted or X-linked genes are aberrantly expressed from both alleles, although standardised analysis methods are still missing. We have developed a tool, named BrewerIX, that provides comprehensive information about the allelic expression of a large, manually-curated set of imprinted and X-linked genes. BrewerIX does not require programming skills, runs on a standard personal computer, and can analyze both bulk and single-cell transcriptomes of human and mouse cells directly from raw sequencing data. BrewerIX confirmed previous observations regarding the bi-allelic expression of some imprinted genes in naive pluripotent cells and extended them to preimplantation embryos. BrewerIX also identified misregulated imprinted genes in breast cancer cells and in human organoids and identified genes escaping XCI in human somatic cells. We believe BrewerIX will be useful for the study of genomic imprinting and XCI during development and reprogramming, and for detecting aberrations in cancer, iPSCs and organoids. Due to its ease of use to non-computational biologists, its implementation could become standard practice during sample assessment, thus raising the robustness and reproducibility of future studies.

• Regulation of gene expression through non-coding RNAs at single myofiber and nucleus resolution. • Reinterpretation of KEGG pathways with microRNA and long non-coding RNA activities. • miR-149, -214, and let-7e alter mitochondrial shape. • The long non-coding RNA Pvt1 is a sponge for miR-27a. • miR-208b regulates Sox6; miR-214 regulates both Sox6 and Slc16a3.