Noradrenaline release from the locus coeruleus shapes stress-induced hippocampal gene expression

Exposure to an acute stressor triggers a complex cascade of neurochemical events in the brain. However, deciphering their individual impact on stress-induced molecular changes remains a major challenge. Here, we combine RNA sequencing with selective pharmacological, chemogenetic, and optogenetic manipulations to isolate the contribution of the locus coeruleus-noradrenaline (LC-NA) system to the acute stress response in mice. We reveal that NA release during stress exposure regulates a large and reproducible set of genes in the dorsal and ventral hippocampus via β-adrenergic receptors. For a smaller subset of these genes, we show that NA release triggered by LC stimulation is sufficient to mimic the stress-induced transcriptional response. We observe these effects in both sexes, and independent of the pattern and frequency of LC activation. Using a retrograde optogenetic approach, we demonstrate that hippocampus-projecting LC neurons directly regulate hippocampal gene expression. Overall, a highly selective set of astrocyte-enriched genes emerges as key targets of LC-NA activation, most prominently several subunits of protein phosphatase 1 (Ppp1r3c, Ppp1r3d, Ppp1r3g) and type II iodothyronine deiodinase (Dio2). These results highlight the importance of astrocytic energy metabolism and thyroid hormone signaling in LC-mediated hippocampal function and offer new molecular targets for understanding how NA impacts brain function in health and disease.

On the identification of differentially-active transcription factors from ATAC-seq data

ATAC-seq has emerged as a rich epigenome profiling technique, and is commonly used to identify Transcription Factors (TFs) underlying given phenomena. A number of methods can be used to identify differentially-active TFs through the accessibility of their DNA-binding motif, however little is known on the best approaches for doing so. Here we benchmark several such methods using a combination of curated datasets with various forms of short-term perturbations on known TFs, as well as semi-simulations. We include both methods specifically designed for this type of data as well as some that can be repurposed for it. We also investigate variations to these methods, and identify three particularly promising approaches (chromVAR-limma with critical adjustments, monaLisa and a combination of GC smooth quantile normalization and multivariate modeling). We further investigate the specific use of nucleosome-free fragments, the combination of top methods, and the impact of technical variation. Finally, we illustrate the use of the top methods on a novel dataset to characterize the impact on DNA accessibility of TRAnscription Factor TArgeting Chimeras (TRAFTAC), which can deplete TFs - in our case NFkB - at the protein level.

Harnessing PROTAC technology to combat stress hormone receptor activation

Counteracting the overactivation of glucocorticoid receptors (GR) is an important therapeutic goal in stress-related psychiatry and beyond. The only clinically approved GR antagonist lacks selectivity and induces unwanted side effects. To complement existing tools of small-molecule-based inhibitors, we present a highly potent, catalytically-driven GR degrader, KH-103, based on proteolysis-targeting chimera technology. This selective degrader enables immediate and reversible GR depletion that is independent of genetic manipulation and circumvents transcriptional adaptations to inhibition. KH-103 achieves passive inhibition, preventing agonistic induction of gene expression, and significantly averts the GR's genomic effects compared to two currently available inhibitors. Application in primary-neuron cultures revealed the dependency of a glucocorticoid-induced increase in spontaneous calcium activity on GR. Finally, we present a proof of concept for application in vivo. KH-103 opens opportunities for a more lucid interpretation of GR functions with translational potential.

Sperm chromatin accessibility's involvement in the intergenerational effects of stress hormone receptor activation

Dexamethasone is a stress hormone receptor agonist used widely in clinics. We and others previously showed that paternal administration of dexamethasone in mice affects the phenotype of their offspring. The substrate of intergenerational transmission of environmentally induced effects often involves changes in sperm RNA, yet other epigenetic modifications in the germline can be affected and are also plausible candidates. First, we tested the involvement of altered sperm RNAs in the transmission of dexamethasone induced phenotypes across generations. We did this by injecting sperm RNA into naïve fertilized oocytes, before performing metabolic and behavioral phenotyping of the offspring. We observed phenotypic changes in discordance with those found in offspring generated by in vitro fertilization using sperm from dexamethasone exposed males. Second, we investigated the effect of dexamethasone on chromatin accessibility using ATAC sequencing and found significant changes at specific genomic features and gene regulatory loci. Employing q-RT-PCR, we show altered expression of a gene in the tissue of offspring affected by accessibility changes in sperm. Third, we establish a correlation between specific DNA modifications and stress hormone receptor activity as a likely contributing factor influencing sperm accessibility. Finally, we independently investigated this dependency by genetically reducing thymine-DNA glycosylase levels and observing concomitant changes at the level of chromatin accessibility and stress hormone receptor activity.

Reuniting philosophy and science to advance cancer research

Cancers rely on multiple, heterogeneous processes at different scales, pertaining to many biomedical fields. Therefore, understanding cancer is necessarily an interdisciplinary task that requires placing specialised experimental and clinical research into a broader conceptual, theoretical, and methodological framework. Without such a framework, oncology will collect piecemeal results, with scant dialogue between the different scientific communities studying cancer. We argue that one important way forward in service of a more successful dialogue is through greater integration of applied sciences (experimental and clinical) with conceptual and theoretical approaches, informed by philosophical methods. By way of illustration, we explore six central themes: (i) the role of mutations in cancer; (ii) the clonal evolution of cancer cells; (iii) the relationship between cancer and multicellularity; (iv) the tumour microenvironment; (v) the immune system; and (vi) stem cells. In each case, we examine open questions in the scientific literature through a philosophical methodology and show the benefit of such a synergy for the scientific and medical understanding of cancer.

The vascular gene Apold1 is dispensable for normal development but controls angiogenesis under pathological conditions

The molecular mechanisms of angiogenesis have been intensely studied, but many genes that control endothelial behavior and fate still need to be described. Here, we characterize the role of Apold1 (Apolipoprotein L domain containing 1) in angiogenesis in vivo and in vitro. Single-cell analyses reveal that - across tissues - the expression of Apold1 is restricted to the vasculature and that Apold1 expression in endothelial cells (ECs) is highly sensitive to environmental factors. Using Apold1-/- mice, we find that Apold1 is dispensable for development and does not affect postnatal retinal angiogenesis nor alters the vascular network in adult brain and muscle. However, when exposed to ischemic conditions following photothrombotic stroke as well as femoral artery ligation, Apold1-/- mice display dramatic impairments in recovery and revascularization. We also find that human tumor endothelial cells express strikingly higher levels of Apold1 and that Apold1 deletion in mice stunts the growth of subcutaneous B16 melanoma tumors, which have smaller and poorly perfused vessels. Mechanistically, Apold1 is activated in ECs upon growth factor stimulation as well as in hypoxia, and Apold1 intrinsically controls EC proliferation but not migration. Our data demonstrate that Apold1 is a key regulator of angiogenesis in pathological settings, whereas it does not affect developmental angiogenesis, thus making it a promising candidate for clinical investigation.

Meta-analysis of (single-cell method) benchmarks reveals the need for extensibility and interoperability

Computational methods represent the lifeblood of modern molecular biology. Benchmarking is important for all methods, but with a focus here on computational methods, benchmarking is critical to dissect important steps of analysis pipelines, formally assess performance across common situations as well as edge cases, and ultimately guide users on what tools to use. Benchmarking can also be important for community building and advancing methods in a principled way. We conducted a meta-analysis of recent single-cell benchmarks to summarize the scope, extensibility, and neutrality, as well as technical features and whether best practices in open data and reproducible research were followed. The results highlight that while benchmarks often make code available and are in principle reproducible, they remain difficult to extend, for example, as new methods and new ways to assess methods emerge. In addition, embracing containerization and workflow systems would enhance reusability of intermediate benchmarking results, thus also driving wider adoption.

enrichMiR predicts functionally relevant microRNAs based on target collections

MicroRNAs (miRNAs) are small non-coding RNAs that are among the main post-transcriptional regulators of gene expression. A number of data collections and prediction tools have gathered putative or confirmed targets of these regulators. It is often useful, for discovery and validation, to harness such collections to perform target enrichment analysis in given transcriptional signatures or gene-sets in order to predict involved miRNAs. While several methods have been proposed to this end, a flexible and user-friendly interface for such analyses using various approaches and collections is lacking. enrichMiR (https://ethz-ins.org/enrichMiR/) addresses this gap by enabling users to perform a series of enrichment tests, based on several target collections, to rank miRNAs according to their likely involvement in the control of a given transcriptional signature or gene-set. enrichMiR results can furthermore be visualised through interactive and publication-ready plots. To guide the choice of the appropriate analysis method, we benchmarked various tests across a panel of experiments involving the perturbation of known miRNAs. Finally, we showcase enrichMiR functionalities in a pair of use cases.

Doublet identification in single-cell sequencing data using scDblFinder

Doublets are prevalent in single-cell sequencing data and can lead to artifactual findings. A number of strategies have therefore been proposed to detect them. Building on the strengths of existing

approaches, we developed scDblFinder, a fast, flexible and accurate Bioconductor-based doublet detection method. Here we present the method, justify its design choices, demonstrate its performance on both single-cell RNA and accessibility (ATAC) sequencing data, and provide some observations on doublet formation, detection, and enrichment analysis. Even in complex datasets, scDblFinder can accurately identify most heterotypic doublets, and was already found by an independent benchmark to outcompete alternatives.

Multiomic profiling of the acute stress response in the mouse hippocampus

The acute stress response mobilizes energy to meet situational demands and re-establish homeostasis. However, the underlying molecular cascades are unclear. Here, we use a brief swim exposure to trigger an acute stress response in mice, which transiently increases anxiety, without leading to lasting maladaptive changes. Using multiomic profiling, such as proteomics, phospho-proteomics, bulk mRNA-, single-nuclei mRNA-, small RNA-, and TRAP-sequencing, we characterize the acute stress-induced molecular events in the mouse hippocampus over time. Our results show the complexity and specificity of the response to acute stress, highlighting both the widespread changes in protein phosphorylation and gene transcription, and tightly regulated protein translation. The observed molecular events resolve efficiently within four hours after initiation of stress. We include an interactive app to explore the data, providing a molecular resource that can help us understand how acute stress impacts brain function in response to stress.

Acute stress can help individuals to respond to challenging events, although chronic stress leads to maladaptive changes. Here, the authors present a multi omic analysis profiling acute stress-induced changes in the mouse hippocampus, providing a resource for the scientific community.

microRNA-138 controls hippocampal interneuron function and short-term memory in mice

The proper development and function of neuronal circuits rely on a tightly regulated balance between excitatory and inhibitory (E/I) synaptic transmission, and disrupting this balance can cause neurodevelopmental disorders, for example, schizophrenia. MicroRNA-dependent gene regulation in pyramidal neurons is important for excitatory synaptic function and cognition, but its role in inhibitory interneurons is poorly understood. Here, we identify miR138-5p as a regulator of short-term memory and inhibitory synaptic transmission in the mouse hippocampus. Sponge-mediated miR138-5p inactivation specifically in mouse parvalbumin (PV)-expressing interneurons impairs spatial recognition memory and enhances GABAergic synaptic input onto pyramidal neurons. Cellular and behavioral phenotypes associated with miR138-5p inactivation are paralleled by an upregulation of the schizophrenia (SCZ)-associated Erbb4, which we validated as a direct miR138-5p target gene. Our findings suggest that miR138-5p is a critical regulator of PV interneuron function in mice, with implications for cognition and SCZ. More generally, they provide evidence that microRNAs orchestrate neural circuit development by fine-tuning both excitatory and inhibitory synaptic transmission.

scanMiR: a biochemically-based toolkit for versatile and efficient microRNA target prediction

MOTIVATION

microRNAs are important post-transcriptional regulators of gene expression, but the identification of functionally relevant targets is still challenging. Recent research has shown improved prediction of microRNA-mediated repression using a biochemical model combined with empirically-derived k-mer affinity predictions, however these findings are not easily applicable.

RESULTS

We translate this approach into a flexible and user-friendly bioconductor package, scanMiR, also available through a web interface. Using lightweight linear models, scanMiR efficiently scans for binding sites, estimates their affinity, and predicts aggregated transcript repression. Moreover, flexible 3'-supplementary alignment enables the prediction of unconventional interactions, such as bindings potentially leading to target-directed microRNA degradation or slicing. We showcase scanMiR through a systematic scan for such unconventional sites on neuronal transcripts, including lncRNAs and circRNAs. Finally, in addition to the main bioconductor package implementing these functions, we provide a user-friendly web application enabling the scanning of sequences, the visualization of predicted bindings, and the browsing of predicted target repression.

AVAILABILITY

scanMiR and companion packages are implemented in R, released under the GPL-3 and accessible on Bioconductor (https://bioconductor.org/packages/release/bioc/html/scanMiR.html) as well as through a shiny web server (https://ethz-ins.org/scanMiR/).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

From cohorts to molecules: Adverse impacts of endocrine disrupting mixtures

Convergent evidence associates exposure to endocrine disrupting chemicals (EDCs) with major human diseases, even at regulation-compliant concentrations. This might be because humans are exposed to EDC mixtures, whereas chemical regulation is based on a risk assessment of individual compounds. Here, we developed a mixture-centered risk assessment strategy that integrates epidemiological and experimental evidence. We identified that exposure to an EDC mixture in early pregnancy is associated with language delay in offspring. At human-relevant concentrations, this mixture disrupted hormone-regulated and disease-relevant regulatory networks in human brain organoids and in the model organisms Xenopus leavis and Danio rerio, as well as behavioral responses. Reinterrogating epidemiological data, we found that up to 54% of the children had prenatal exposures above experimentally derived levels of concern, reaching, for the upper decile compared with the lowest decile of exposure, a 3.3 times higher risk of language delay.

Doublet identification in single-cell sequencing data using scDblFinder

Doublets are prevalent in single-cell sequencing data and can lead to artifactual findings. A number of strategies have therefore been proposed to detect them. Building on the strengths of existing approaches, we developed scDblFinder, a fast, flexible and accurate Bioconductor-based doublet detection method. Here we present the method, justify its design choices, demonstrate its performance on both single-cell RNA and accessibility (ATAC) sequencing data, and provide some observations on doublet formation, detection, and enrichment analysis. Even in complex datasets, scDblFinder can accurately identify most heterotypic doublets, and was already found by an independent benchmark to outcompete alternatives.

Pervasive compartment-specific regulation of gene expression during homeostatic synaptic scaling

Synaptic scaling is a form of homeostatic plasticity which allows neurons to adjust their action potential firing rate in response to chronic alterations in neural activity. Synaptic scaling requires profound changes in gene expression, but the relative contribution of local and cell‐wide mechanisms is controversial. Here we perform a comprehensive multi‐omics characterization of the somatic and process compartments of primary rat hippocampal neurons during synaptic scaling. We uncover both highly compartment‐specific and correlating changes in the neuronal transcriptome and proteome. Whereas downregulation of crucial regulators of neuronal excitability occurs primarily in the somatic compartment, structural components of excitatory postsynapses are mostly downregulated in processes. Local inhibition of protein synthesis in processes during scaling is confirmed for candidate synaptic proteins. Motif analysis further suggests an important role for trans‐acting post‐transcriptional regulators, including RNA‐binding proteins and microRNAs, in the local regulation of the corresponding mRNAs. Altogether, our study indicates that, during synaptic scaling, compartmentalized gene expression changes might co‐exist with neuron‐wide mechanisms to allow synaptic computation and homeostasis.

The Acute Stress Response in the Multiomic Era

Studying the stress response is a major pillar of neuroscience research not only because stress is a daily reality but also because the exquisitely fine-tuned bodily changes triggered by stress are a neuroendocrinological marvel. While the genome-wide changes induced by chronic stress have been extensively studied, we know surprisingly little about the complex molecular cascades triggered by acute stressors, the building blocks of chronic stress. The acute stress (or fight-or-flight) response mobilizes organismal energy resources to meet situational demands. However, successful stress coping also requires the efficient termination of the stress response. Maladaptive coping-particularly in response to severe or repeated stressors-can lead to allostatic (over)load, causing wear and tear on tissues, exhaustion, and disease. We propose that deep molecular profiling of the changes triggered by acute stressors could provide molecular correlates for allostatic load and predict healthy or maladaptive stress responses. We present a theoretical framework to interpret multiomic data in light of energy homeostasis and activity-dependent gene regulation, and we review the signaling cascades and molecular changes rapidly induced by acute stress in different cell types in the brain. In addition, we review and reanalyze recent data from multiomic screens conducted mainly in the rodent hippocampus and amygdala after acute psychophysical stressors. We identify challenges surrounding experimental design and data analysis, and we highlight promising new research directions to better understand the stress response on a multiomic level.

treeclimbR pinpoints the data-dependent resolution of hierarchical hypotheses

treeclimbR is for analyzing hierarchical trees of entities, such as phylogenies or cell types, at different resolutions. It proposes multiple candidates that capture the latent signal and pinpoints branches or leaves that contain features of interest, in a data-driven way. It outperforms currently available methods on synthetic data, and we highlight the approach on various applications, including microbiome and microRNA surveys as well as single-cell cytometry and RNA-seq datasets. With the emergence of various multi-resolution genomic datasets, treeclimbR provides a thorough inspection on entities across resolutions and gives additional flexibility to uncover biological associations.

Streamlining differential exon and 3' UTR usage with diffUTR

Background

Despite the importance of alternative poly-adenylation and 3′ UTR length for a variety of biological phenomena, there are limited means of detecting UTR changes from standard transcriptomic data.

Results

We present the diffUTR Bioconductor package which streamlines and improves upon differential exon usage (DEU) analyses, and leverages existing DEU tools and alternative poly-adenylation site databases to enable differential 3′ UTR usage analysis. We demonstrate the diffUTR features and show that it is more flexible and more accurate than state-of-the-art alternatives, both in simulations and in real data.

Conclusions

diffUTR enables differential 3′ UTR analysis and more generally facilitates DEU and the exploration of their results.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04114-7.

Isoform-resolved correlation analysis between mRNA abundance regulation and protein level degradation

Profiling of biological relationships between different molecular layers dissects regulatory mechanisms that ultimately determine cellular function. To thoroughly assess the role of protein post‐translational turnover, we devised a strategy combining pulse stable isotope‐labeled amino acids in cells (pSILAC), data‐independent acquisition mass spectrometry (DIA‐MS), and a novel data analysis framework that resolves protein degradation rate on the level of mRNA alternative splicing isoforms and isoform groups. We demonstrated our approach by the genome‐wide correlation analysis between mRNA amounts and protein degradation across different strains of HeLa cells that harbor a high grade of gene dosage variation. The dataset revealed that specific biological processes, cellular organelles, spatial compartments of organelles, and individual protein isoforms of the same genes could have distinctive degradation rate. The protein degradation diversity thus dissects the corresponding buffering or concerting protein turnover control across cancer cell lines. The data further indicate that specific mRNA splicing events such as intron retention significantly impact the protein abundance levels. Our findings support the tight association between transcriptome variability and proteostasis and provide a methodological foundation for studying functional protein degradation.

muscat detects subpopulation-specific state transitions from multi-sample multi-condition single-cell transcriptomics data

Single-cell RNA sequencing (scRNA-seq) has become an empowering technology to profile the transcriptomes of individual cells on a large scale. Early analyses of differential expression have aimed at identifying differences between subpopulations to identify subpopulation markers. More generally, such methods compare expression levels across sets of cells, thus leading to cross-condition analyses. Given the emergence of replicated multi-condition scRNA-seq datasets, an area of increasing focus is making sample-level inferences, termed here as differential state analysis; however, it is not clear which statistical framework best handles this situation. Here, we surveyed methods to perform cross-condition differential state analyses, including cell-level mixed models and methods based on aggregated pseudobulk data. To evaluate method performance, we developed a flexible simulation that mimics multi-sample scRNA-seq data. We analyzed scRNA-seq data from mouse cortex cells to uncover subpopulation-specific responses to lipopolysaccharide treatment, and provide robust tools for multi-condition analysis within the muscat R package.

Alterations in sperm long RNA contribute to the epigenetic inheritance of the effects of postnatal trauma

Psychiatric diseases have a strong heritable component known to not be restricted to DNA sequence-based genetic inheritance alone but to also involve epigenetic factors in germ cells. Initial evidence suggested that sperm RNA is causally linked to the transmission of symptoms induced by traumatic experiences. Here, we show that alterations in long RNA in sperm contribute to the inheritance of specific trauma symptoms. Injection of long RNA fraction from sperm of males exposed to postnatal trauma recapitulates the effects on food intake, glucose response to insulin and risk-taking in adulthood whereas the small RNA fraction alters body weight and behavioural despair. Alterations in long RNA are maintained after fertilization, suggesting a direct link between sperm and embryo RNA.

pipeComp, a general framework for the evaluation of computational pipelines, reveals performant single cell RNA-seq preprocessing tools

We present pipeComp ( https://github.com/plger/pipeComp ), a flexible R framework for pipeline comparison handling interactions between analysis steps and relying on multi-level evaluation metrics. We apply it to the benchmark of single-cell RNA-sequencing analysis pipelines using simulated and real datasets with known cell identities, covering common methods of filtering, doublet detection, normalization, feature selection, denoising, dimensionality reduction, and clustering. pipeComp can easily integrate any other step, tool, or evaluation metric, allowing extensible benchmarks and easy applications to other fields, as we demonstrate through a study of the impact of removal of unwanted variation on differential expression analysis.

A complete pupillometry toolbox for real-time monitoring of locus coeruleus activity in rodents

The locus coeruleus (LC) is a region in the brainstem that produces noradrenaline and is involved in both normal and pathological brain function. Pupillometry, the measurement of pupil diameter, provides a powerful readout of LC activity in rodents, primates and humans. The protocol detailed here describes a miniaturized setup that can screen LC activity in rodents in real-time and can be established within 1-2 d. Using low-cost Raspberry Pi computers and cameras, the complete custom-built system costs only ~300 euros, is compatible with stereotaxic surgery frames and seamlessly integrates into complex experimental setups. Tools for pupil tracking and a user-friendly Pupillometry App allow quantification, analysis and visualization of pupil size. Pupillometry can discriminate between different, physiologically relevant firing patterns of the LC and can accurately report LC activation as measured by noradrenaline turnover. Pupillometry provides a rapid, non-invasive readout that can be used to verify accurate placement of electrodes/fibers in vivo, thus allowing decisions about the inclusion/exclusion of individual animals before experiments begin.

Dosage analysis of the 7q11.23 Williams region identifies BAZ1B as a major human gene patterning the modern human face and underlying self-domestication

We undertook a functional dissection of chromatin remodeler BAZ1B in neural crest (NC) stem cells (NCSCs) from a uniquely informative cohort of typical and atypical patients harboring 7q11.23 copy number variants. Our results reveal a key contribution of BAZ1B to NCSC in vitro induction and migration, coupled with a crucial involvement in NC-specific transcriptional circuits and distal regulation. By intersecting our experimental data with new paleogenetic analyses comparing modern and archaic humans, we found a modern-specific enrichment for regulatory changes both in BAZ1B and its experimentally defined downstream targets, thereby providing the first empirical validation of the human self-domestication hypothesis and positioning BAZ1B as a master regulator of the modern human face. In so doing, we provide experimental evidence that the craniofacial and cognitive/behavioral phenotypes caused by alterations of the Williams-Beuren syndrome critical region can serve as a powerful entry point into the evolution of the modern human face and prosociality.

KMT2B Is Selectively Required for Neuronal Transdifferentiation, and Its Loss Exposes Dystonia Candidate Genes

Transdifferentiation of fibroblasts into induced neuronal cells (iNs) by the neuron-specific transcription factors Brn2, Myt1l, and Ascl1 is a paradigmatic example of inter-lineage conversion across epigenetically distant cells. Despite tremendous progress regarding the transcriptional hierarchy underlying transdifferentiation, the enablers of the concomitant epigenome resetting remain to be elucidated. Here, we investigated the role of KMT2A and KMT2B, two histone H3 lysine 4 methylases with cardinal roles in development, through individual and combined inactivation. We found that Kmt2b, whose human homolog's mutations cause dystonia, is selectively required for iN conversion through suppression of the alternative myocyte program and induction of neuronal maturation genes. The identification of KMT2B-vulnerable targets allowed us, in turn, to expose, in a cohort of 225 patients, 45 unique variants in 39 KMT2B targets, which represent promising candidates to dissect the molecular bases of dystonia.

Multi-omic measurements of heterogeneity in HeLa cells across laboratories

Reproducibility in research can be compromised by both biological and technical variation, but most of the focus is on removing the latter. Here we investigate the effects of biological variation in HeLa cell lines using a systems-wide approach. We determine the degree of molecular and phenotypic variability across 14 stock HeLa samples from 13 international laboratories. We cultured cells in uniform conditions and profiled genome-wide copy numbers, mRNAs, proteins and protein turnover rates in each cell line. We discovered substantial heterogeneity between HeLa variants, especially between lines of the CCL2 and Kyoto varieties, and observed progressive divergence within a specific cell line over 50 successive passages. Genomic variability has a complex, nonlinear effect on transcriptome, proteome and protein turnover profiles, and proteotype patterns explain the varying phenotypic response of different cell lines to Salmonella infection. These findings have implications for the interpretation and reproducibility of research results obtained from human cultured cells.

Exploratory rearing: a context- and stress-sensitive behavior recorded in the open-field test

Stressful experiences are linked to anxiety disorders in humans. Similar effects are observed in rodent models, where anxiety is often measured in classic conflict tests such as the open-field test. Spontaneous rearing behavior, in which rodents stand on their hind legs to explore, can also be observed in this test yet is often ignored. We define two forms of rearing, supported rearing (in which the animal rears against the walls of the arena) and unsupported rearing (in which the animal rears without contacting the walls of the arena). Using an automated open-field test, we show that both rearing behaviors appear to be strongly context dependent and show clear sex differences, with females rearing less than males. We show that unsupported rearing is sensitive to acute stress, and is reduced under more averse testing conditions. Repeated testing and handling procedures lead to changes in several parameters over varying test sessions, yet unsupported rearing appears to be rather stable within a given animal. Rearing behaviors could therefore provide an additional measure of anxiety in rodents relevant for behavioral studies, as they appear to be highly sensitive to context and may be used in repeated testing designs.

The European politics of animal experimentation: From Victorian Britain to 'Stop Vivisection'

This paper identifies a common political struggle behind debates on the validity and permissibility of animal experimentation, through an analysis of two recent European case studies: the Italian implementation of the European Directive 2010/63/EC regulating the use of animals in science, and the recent European Citizens' Initiative (ECI) 'Stop Vivisection'. Drawing from a historical parallel with Victorian antivivisectionism, we highlight important threads in our case studies that mark the often neglected specificities of debates on animal experimentation. From the representation of the sadistic scientist in the XIX century, to his/her claimed capture by vested interests and evasion of public scrutiny in the contemporary cases, we show that animals are not simply the focus of the debate, but also a privileged locus at which much broader issues are being raised about science, its authority, accountability and potential misalignment with public interest. By highlighting this common socio-political conflict underlying public controversies around animal experimentation, our work prompts the exploration of modes of authority and argumentation that, in establishing the usefulness of animals in science, avoid reenacting the traditional divide between epistemic and political fora.

YY1 Haploinsufficiency Causes an Intellectual Disability Syndrome Featuring Transcriptional and Chromatin Dysfunction

Yin and yang 1 (YY1) is a well-known zinc-finger transcription factor with crucial roles in normal development and malignancy. YY1 acts both as a repressor and as an activator of gene expression. We have identified 23 individuals with de novo mutations or deletions of YY1 and phenotypic features that define a syndrome of cognitive impairment, behavioral alterations, intrauterine growth restriction, feeding problems, and various congenital malformations. Our combined clinical and molecular data define "YY1 syndrome" as a haploinsufficiency syndrome. Through immunoprecipitation of YY1-bound chromatin from affected individuals' cells with antibodies recognizing both ends of the protein, we show that YY1 deletions and missense mutations lead to a global loss of YY1 binding with a preferential retention at high-occupancy sites. Finally, we uncover a widespread loss of H3K27 acetylation in particular on the YY1-bound enhancers, underscoring a crucial role for YY1 in enhancer regulation. Collectively, these results define a clinical syndrome caused by haploinsufficiency of YY1 through dysregulation of key transcriptional regulators.

RNAontheBENCH: computational and empirical resources for benchmarking RNAseq quantification and differential expression methods

RNA sequencing (RNAseq) has become the method of choice for transcriptome analysis, yet no consensus exists as to the most appropriate pipeline for its analysis, with current benchmarks suffering important limitations. Here, we address these challenges through a rich benchmarking resource harnessing (i) two RNAseq datasets including ERCC ExFold spike-ins; (ii) Nanostring measurements of a panel of 150 genes on the same samples; (iii) a set of internal, genetically-determined controls; (iv) a reanalysis of the SEQC dataset; and (v) a focus on relative quantification (i.e. across-samples). We use this resource to compare different approaches to each step of RNAseq analysis, from alignment to differential expression testing. We show that methods providing the best absolute quantification do not necessarily provide good relative quantification across samples, that count-based methods are superior for gene-level relative quantification, and that the new generation of pseudo-alignment-based software performs as well as established methods, at a fraction of the computing time. We also assess the impact of library type and size on quantification and differential expression analysis. Finally, we have created a R package and a web platform to enable the simple and streamlined application of this resource to the benchmarking of future methods.

From replica to instruments: animal models in biomedical research

The ways in which other animal species can be informative about human biology are not exhausted by the traditional picture of the animal model. In this paper, I propose to distinguish two roles which laboratory organisms can have in biomedical research. In the more traditional case, organisms act as surrogates for human beings, and as such are expected to be more manageable replicas of humans. However, animal models can inform us about human biology in a much less straightforward way, by being used as measuring devices-what I call their instrumental role. I first characterize this role and provide criteria for it, before illustrating it with some examples from biomedical research, especially cancer research. In such an instrumental role, phenotypes are not expected to phenocopy human phenomena, but instead have the purely instrumental value of detecting or measuring differences. I argue that the instrumental role is more prevalent than might first be suspected, and that some characteristics of contemporary biomedical research are increasingly shifting the use of laboratory organisms to the instrumental role. Finally, in light of the distinction proposed, I discuss the meaning of the expression "animal model".

Germinal center dysregulation by histone methyltransferase EZH2 promotes lymphomagenesis

Protection against deadly pathogens requires the production of high-affinity antibodies by B cells, which are generated in germinal centers (GCs). Alteration of the GC developmental program is common in many B cell malignancies. Identification of regulators of the GC response is crucial to develop targeted therapies for GC B cell dysfunctions, including lymphomas. The histone H3 lysine 27 methyltransferase enhancer of zeste homolog 2 (EZH2) is highly expressed in GC B cells and is often constitutively activated in GC-derived non-Hodgkin lymphomas (NHLs). The function of EZH2 in GC B cells remains largely unknown. Herein, we show that Ezh2 inactivation in mouse GC B cells caused profound impairment of GC responses, memory B cell formation, and humoral immunity. EZH2 protected GC B cells against activation-induced cytidine deaminase (AID) mutagenesis, facilitated cell cycle progression, and silenced plasma cell determinant and tumor suppressor B-lymphocyte-induced maturation protein 1 (BLIMP1). EZH2 inhibition in NHL cells induced BLIMP1, which impaired tumor growth. In conclusion, EZH2 sustains AID function and prevents terminal differentiation of GC B cells, which allows antibody diversification and affinity maturation. Dysregulation of the GC reaction by constitutively active EZH2 facilitates lymphomagenesis and identifies EZH2 as a possible therapeutic target in NHL and other GC-derived B cell diseases.

What mechanisms can't do: explanatory frameworks and the function of the p53 gene in molecular oncology

What has been called the new mechanistic philosophy conceives of mechanisms as the main providers of biological explanation. We draw on the characterization of the p53 gene in molecular oncology, to show that explaining a biological phenomenon (cancer, in our case) implies instead a dynamic interaction between the mechanistic level-rendered at the appropriate degree of ontological resolution-and far more general explanatory tools that perform a fundamental epistemic role in the provision of biological explanations. We call such tools "explanatory frameworks". They are called frameworks to stress their higher level of generality with respect to bare mechanisms; on the other hand, they are called explanatory because, as we show in this paper, their importance in explaining biological phenomena is not secondary with respect to mechanisms. We illustrate how explanatory frameworks establish selective and local criteria of causal relevance that drive the search for, characterisation and usage of biological mechanisms. Furthermore, we show that explanatory frameworks allow for changes of scientific perspective on the causal relevance of mechanisms going beyond the account provided by the new mechanistic philosophy.

Cell reprogramming requires silencing of a core subset of polycomb targets

Transcription factor (TF)–induced reprogramming of somatic cells into induced pluripotent stem cells (iPSC) is associated with genome-wide changes in chromatin modifications. Polycomb-mediated histone H3 lysine-27 trimethylation (H3K27me3) has been proposed as a defining mark that distinguishes the somatic from the iPSC epigenome. Here, we dissected the functional role of H3K27me3 in TF–induced reprogramming through the inactivation of the H3K27 methylase EZH2 at the onset of reprogramming. Our results demonstrate that surprisingly the establishment of functional iPSC proceeds despite global loss of H3K27me3. iPSC lacking EZH2 efficiently silenced the somatic transcriptome and differentiated into tissues derived from the three germ layers. Remarkably, the genome-wide analysis of H3K27me3 in Ezh2 mutant iPSC cells revealed the retention of this mark on a highly selected group of Polycomb targets enriched for developmental regulators controlling the expression of lineage specific genes. Erasure of H3K27me3 from these targets led to a striking impairment in TF–induced reprogramming. These results indicate that PRC2-mediated H3K27 trimethylation is required on a highly selective core of Polycomb targets whose repression enables TF–dependent cell reprogramming.

Multicellular organisms are composed of a variety of cell types. Over the last years we have learned that cell differentiation is fully reversible and that it takes few specific transcription factors (proteins that bind to DNA and regulate gene expression) to convert one cell type into another. The most dramatic example is the reprogramming of somatic cells into stem cells that reacquire the potential to give rise to all cell types of the body. This process entails the resetting of the gene expression program of the somatic cells necessary to acquire a pluripotent state but remains poorly understood. Here, we defined the role in cell reprogramming of the Polycomb axis, one of the key effectors of gene silencing that operates through a specific chemical modification of histones (the trimethylation of lysine-27 on histone H3) that is stably inherited from one cell generation to the next. Surprisingly, we uncovered a selective requirement for Polycomb silencing during somatic cell reprogramming. Upon inactivation of the essential Polycomb protein EZH2, histones trimethylated on lysine-27 were retained on a selected subset of physiological Polycomb target genes whose products control normal development, and erasure of the histone mark from these critical genes abolished reprogramming.