Comparison of alternative approaches for analysing multi-level RNA-seq data

RNA sequencing (RNA-seq) is widely used for RNA quantification in the environmental, biological and medical sciences. It enables the description of genome-wide patterns of expression and the identification of regulatory interactions and networks. The aim of RNA-seq data analyses is to achieve rigorous quantification of genes/transcripts to allow a reliable prediction of differential expression (DE), despite variation in levels of noise and inherent biases in sequencing data. This can be especially challenging for datasets in which gene expression differences are subtle, as in the behavioural transcriptomics test dataset from D. melanogaster that we used here. We investigated the power of existing approaches for quality checking mRNA-seq data and explored additional, quantitative quality checks. To accommodate nested, multi-level experimental designs, we incorporated sample layout into our analyses. We employed a subsampling without replacement-based normalization and an identification of DE that accounted for the hierarchy and amplitude of effect sizes within samples, then evaluated the resulting differential expression call in comparison to existing approaches. In a final step to test for broader applicability, we applied our approaches to a published set of H. sapiens mRNA-seq samples, The dataset-tailored methods improved sample comparability and delivered a robust prediction of subtle gene expression changes. The proposed approaches have the potential to improve key steps in the analysis of RNA-seq data by incorporating the structure and characteristics of biological experiments

Inhibitory top-down projections from zona incerta mediate neocortical memory

Top-down projections convey a family of signals encoding previous experiences and current aims to the sensory neocortex, where they converge with external bottom-up information to enable perception and memory. Whereas top-down control has been attributed to excitatory pathways, the existence, connectivity, and information content of inhibitory top-down projections remain elusive. Here, we combine synaptic two-photon calcium imaging, circuit mapping, cortex-dependent learning, and chemogenetics in mice to identify GABAergic afferents from the subthalamic zona incerta as a major source of top-down input to the neocortex. Incertocortical transmission undergoes robust plasticity during learning that improves information transfer and mediates behavioral memory. Unlike excitatory pathways, incertocortical afferents form a disinhibitory circuit that encodes learned top-down relevance in a bidirectional manner where the rapid appearance of negative responses serves as the main driver of changes in stimulus representation. Our results therefore reveal the distinctive contribution of long-range (dis)inhibitory afferents to the computational flexibility of neocortical circuits

Microguards and micromessengers of the genome

The regulation of gene expression is of fundamental importance to maintain organismal function and integrity and requires a multifaceted and highly ordered sequence of events. The cyclic nature of gene expression is known as ‘transcription dynamics’. Disruption or perturbation of these dynamics can result in significant fitness costs arising from genome instability, accelerated ageing and disease. We review recent research that supports the idea that an important new role for small RNAs, particularly microRNAs (miRNAs), is in protecting the genome against short-term transcriptional fluctuations, in a process we term ‘microguarding’. An additional emerging role for miRNAs is as ‘micromessengers’—through alteration of gene expression in target cells to which they are trafficked within microvesicles. We describe the scant but emerging evidence that miRNAs can be moved between different cells, individuals and even species, to exert biologically significant responses. With these two new roles, miRNAs have the potential to protect against deleterious gene expression variation from perturbation and to themselves perturb the expression of genes in target cells. These interactions between cells will frequently be subject to conflicts of interest when they occur between unrelated cells that lack a coincidence of fitness interests. Hence, there is the potential for miRNAs to represent both a means to resolve conflicts of interest, as well as instigate them. We conclude by exploring this conflict hypothesis, by describing some of the initial evidence consistent with it and proposing new ideas for future research into this exciting topic

MicroRNA regulation of the paired-box transcription factor Pax3 confers robustness to developmental timing of myogenesis

Commitment of progenitors in the dermomyotome to myoblast fate is the first step in establishing the body musculature. Pax3 is a crucial transcription factor, important for skeletal muscle development and expressed in myogenic progenitors in the dermomyotome of developing somites and in migratory muscle progenitors that populate the limb buds. Down-regulation of Pax3 is essential to ignite the myogenic program, including up-regulation of myogenic regulators, Myf-5 and MyoD. MicroRNAs (miRNAs) confer robustness to developmental timing by posttranscriptional repression of genetic programs that are related to previous developmental stages or to alternative cell fates. Here we demonstrate that the muscle-specific miRNAs miR-1 and miR-206 directly target Pax3. Antagomir-mediated inhibition of miR-1/miR-206 led to delayed myogenic differentiation in developing somites, as shown by transient loss of myogenin expression. This correlated with increased Pax3 and was phenocopied using Pax3-specific target protectors. Loss of myogenin after antagomir injection was rescued by Pax3 knockdown using a splice morpholino, suggesting that miR-1/miR-206 control somite myogenesis primarily through interactions with Pax3. Our studies reveal an important role for miR-1/miR-206 in providing precision to the timing of somite myogenesis. We propose that posttranscriptional control of Pax3 downstream of miR-1/miR-206 is required to stabilize myoblast commitment and subsequent differentiation. Given that mutually exclusive expression of miRNAs and their targets is a prevailing theme in development, our findings suggest that miRNA may provide a general mechanism for the unequivocal commitment underlying stem cell differentiation

ClimBar : An Integrated Approach to Evaluate and Utilize Genetic Diversity

European agriculture anticipates an unprecedented combination of stress factors, production threats and quality needs due to climate change. Various regions of Europe will be affected differently. Barley & wheat domestication, and landrace formation in Europe, were under very different climates than those emerging now. Alleles needed for sustainable, resilient, quality yields in a changed climate are likely not combined in current haplotypes of elite barley cultivars. These alleles are likely found in diverse landraces and wild relatives in the Mediterranean basin and Fertile Crescent -- areas that prefigure expected climate change. New precision, high-throughput phenotyping tools are essential to find trait-allele associations needed for future-climate breeding. Combining genetics, genomics, modelling, molecular biology, morphology, and physiology, ClimBar takes an interdisciplinary approach to develop a strategy for breeding an increased resilience to climate change in barley. ClimBar, a new project under the framework of FACCE ERA-NET Plus Joint Programming Initiative on Climate Smart Agriculture, will identify genome regions, genes, and alleles conferring the traits needed to breed resilient barley varieties adapted to the climatic conditions predicted for 2070 in different European environments. Adapted, resilient germplasm created using ClimBar data, tools and models will provide food-chain security, economic stability and environmental sustainability. Website: http://plen.ku.dk/english/research/plant_soil/breeding/quality/climbar

Detecting new microRNAs in human osteoarthritic chondrocytes identifies miR-3085 as a human, chondrocyte-selective, microRNA

Objective: To use deep sequencing to identify novel microRNAs in human osteoarthritic cartilage which have a functional role in chondrocyte phenotype or function. Design: A small RNA library was prepared from human osteoarthritic primary chondrocytes using in-house adaptors and analysed by Illumina sequencing. Novel candidate microRNAs were validated by northern blot and qRT-PCR. Expression was measured in cartilage models. Targets of novel candidates were identified by microarray and computational analysis, validated using 3’-UTR-luciferase reporter plasmids. Protein levels were assessed by western blot and functional analysis by cell adhesion. Results: We identified 990 known microRNAs and 1621 potential novel microRNAs in human osteoarthritic chondrocytes, 60 of the latter were expressed in all samples assayed. MicroRNA-140-3p was the most highly expressed microRNA in osteoarthritic cartilage. Sixteen novel candidate microRNAs were analysed further, of which 6 remained after northern blot analysis. Three novel microRNAs were regulated across models of chondrogenesis, chondrocyte differentiation or cartilage injury. One sequence (novel #11), annotated in rodents as microRNA-3085-3p, was preferentially expressed in cartilage, dependent on chondrocyte differentiation and, in man, is located in an intron of the cartilage-expressed gene CRTAC-1. This microRNA was shown to target the ITGA5 gene directly (which encodes integrin alpha5) and inhibited adhesion to fibronectin (dependent on alpha5beta1 integrin). Conclusion: Deep sequencing has uncovered many potential microRNA candidates expressed in human cartilage. At least three of these show potential functional interest in cartilage homeostasis and osteoarthritis. Particularly, novel #11 (microRNA-3085-3p) which has been identified for the first time in man

Targeting the MAPK7/MMP9 axis for metastasis in primary bone cancer

Metastasis is the leading cause of cancer-related death. This multistage process involves contribution from both tumour cells and the tumour stroma to release metastatic cells into the circulation. Circulating tumour cells (CTCs) survive circulatory cytotoxicity, extravasate and colonise secondary sites effecting metastatic outcome. Reprogramming the transcriptomic landscape is a metastatic hallmark, but detecting underlying master regulators that drive pathological gene expression is a key challenge, especially in childhood cancer. Here we used whole tumour plus single-cell RNA-sequencing in primary bone cancer and CTCs to perform weighted gene co-expression network analysis to systematically detect coordinated changes in metastatic transcript expression. This approach with comparisons applied to data collected from cell line models, clinical samples and xenograft mouse models revealed mitogen-activated protein kinase 7/matrix metallopeptidase 9 (MAPK7/MMP9) signalling as a driver for primary bone cancer metastasis. RNA interference knockdown of MAPK7 reduces proliferation, colony formation, migration, tumour growth, macrophage residency/polarisation and lung metastasis. Parallel to these observations were reduction of activated interleukins IL1B, IL6, IL8 plus mesenchymal markers VIM and VEGF in response to MAPK7 loss. Our results implicate a newly discovered, multidimensional MAPK7/MMP9 signalling hub in primary bone cancer metastasis that is clinically actionable

Transfer RNA-derived small RNAs in the cancer transcriptome

The cellular lifetime includes stages such as differentiation, proliferation, division, senescence and apoptosis.These stages are driven by a strictly ordered process of transcription dynamics. Molecular disruption to RNA polymerase assembly, chromatin remodelling and transcription factor binding through to RNA editing, splicing, post-transcriptional regulation and ribosome scanning can result in significant costs arising from genome instability. Cancer development is one example of when such disruption takes place. RNA silencing is a term used to describe the effects of post-transcriptional gene silencing mediated by a diverse set of small RNA molecules. Small RNAs are crucial for regulating gene expression and microguarding genome integrity.RNA silencing studies predominantly focus on small RNAs such as microRNAs, short-interfering RNAs and piwi-interacting RNAs. We describe an emerging renewal of inter-est in a‘larger’small RNA, the transfer RNA (tRNA).Precisely generated tRNA-derived small RNAs, named tRNA halves (tiRNAs) and tRNA fragments (tRFs), have been reported to be abundant with dysregulation associated with cancer. Transfection of tiRNAs inhibits protein translation by displacing eukaryotic initiation factors from messenger RNA (mRNA) and inaugurating stress granule formation.Knockdown of an overexpressed tRF inhibits cancer cell proliferation. Recovery of lacking tRFs prevents cancer metastasis. The dual oncogenic and tumour-suppressive role is typical of functional small RNAs. We review recent reports on tiRNA and tRF discovery and biogenesis, identification and analysis from next-generation sequencing data and a mechanistic animal study to demonstrate their physiological role in cancer biology. We propose tRNA-derived small RNA-mediated RNA silencing is an innate defence mechanism to prevent oncogenic translation. We expect that cancer cells are percipient to their ablated control of transcription and attempt to prevent loss of genome control through RNA silencing

Assessment of genetically modified maize MIR604 for renewal authorisation under Regulation (EC) No 1829/2003 (application EFSA‐GMO‐RX‐013)

Following the submission of application EFSA‐GMO‐RX‐013 under Regulation (EC) No 1829/2003 from Syngenta Crop Protection NV/SA, the EFSA Panel on Genetically Modified Organisms (GMO) was asked to deliver a scientific risk assessment on the data submitted in the context of the renewal of authorisation application for the insect‐resistant genetically modified maize MIR604, for food and feed uses, excluding cultivation within the EU. The data received in the context of this renewal application contained post‐market environmental monitoring reports, a systematic search and evaluation of literature, updated bioinformatic analyses, and additional documents or studies performed by or on behalf of the applicant. The GMO Panel assessed these data for possible new hazards, modified exposure or new scientific uncertainties identified during the authorisation period and not previously assessed in the context of the original application. Under the assumption that the DNA sequence of the event in maize MIR604 considered for renewal is identical to the corrected sequence of the originally assessed event, the GMO Panel concludes that there is no evidence in renewal application EFSA‐GMO‐RX‐013 for new hazards, modified exposure or scientific uncertainties that would change the conclusions of the original risk assessment on maize MIR604