Expression Atlas update: insights from sequencing data at both bulk and single cell level

Expression Atlas (www.ebi.ac.uk/gxa) and its newest counterpart the Single Cell Expression Atlas (www.ebi.ac.uk/gxa/sc) are EMBL-EBI's knowledgebases for gene and protein expression and localisation in bulk and at single cell level. These resources aim to allow users to investigate their expression in normal tissue (baseline) or in response to perturbations such as disease or changes to genotype (differential) across multiple species. Users are invited to search for genes or metadata terms across species or biological conditions in a standardised consistent interface. Alongside these data, new features in Single Cell Expression Atlas allow users to query metadata through our new cell type wheel search. At the experiment level data can be explored through two types of dimensionality reduction plots, t-distributed Stochastic Neighbor Embedding (tSNE) and Uniform Manifold Approximation and Projection (UMAP), overlaid with either clustering or metadata information to assist users' understanding. Data are also visualised as marker gene heatmaps identifying genes that help confer cluster identity. For some data, additional visualisations are available as interactive cell level anatomograms and cell type gene expression heatmaps.

Expression Atlas update: gene and protein expression in multiple species

The EMBL-EBI Expression Atlas is an added value knowledge base that enables researchers to answer the question of where (tissue, organism part, developmental stage, cell type) and under which conditions (disease, treatment, gender, etc) a gene or protein of interest is expressed. Expression Atlas brings together data from >4500 expression studies from >65 different species, across different conditions and tissues. It makes these data freely available in an easy to visualise form, after expert curation to accurately represent the intended experimental design, re-analysed via standardised pipelines that rely on open-source community developed tools. Each study's metadata are annotated using ontologies. The data are re-analyzed with the aim of reproducing the original conclusions of the underlying experiments. Expression Atlas is currently divided into Bulk Expression Atlas and Single Cell Expression Atlas. Expression Atlas contains data from differential studies (microarray and bulk RNA-Seq) and baseline studies (bulk RNA-Seq and proteomics), whereas Single Cell Expression Atlas is currently dedicated to Single Cell RNA-Sequencing (scRNA-Seq) studies. The resource has been in continuous development since 2009 and it is available at https://www.ebi.ac.uk/gxa.

From ArrayExpress to BioStudies

ArrayExpress (https://www.ebi.ac.uk/arrayexpress) is an archive of functional genomics data at EMBL-EBI, established in 2002, initially as an archive for publication-related microarray data and was later extended to accept sequencing-based data. Over the last decade an increasing share of biological experiments involve multiple technologies assaying different biological modalities, such as epigenetics, and RNA and protein expression, and thus the BioStudies database (https://www.ebi.ac.uk/biostudies) was established to deal with such multimodal data. Its central concept is a study, which typically is associated with a publication. BioStudies stores metadata describing the study, provides links to the relevant databases, such as European Nucleotide Archive (ENA), as well as hosts the types of data for which specialized databases do not exist. With BioStudies now fully functional, we are able to further harmonize the archival data infrastructure at EMBL-EBI, and ArrayExpress is being migrated to BioStudies. In future, all functional genomics data will be archived at BioStudies. The process will be seamless for the users, who will continue to submit data using the online tool Annotare and will be able to query and download data largely in the same manner as before. Nevertheless, some technical aspects, particularly programmatic access, will change. This update guides the users through these changes.

Ensembl Genomes 2020-enabling non-vertebrate genomic research

Ensembl Genomes (http://www.ensemblgenomes.org) is an integrating resource for genome-scale data from non-vertebrate species, complementing the resources for vertebrate genomics developed in the context of the Ensembl project (http://www.ensembl.org). Together, the two resources provide a consistent set of interfaces to genomic data across the tree of life, including reference genome sequence, gene models, transcriptional data, genetic variation and comparative analysis. Data may be accessed via our website, online tools platform and programmatic interfaces, with updates made four times per year (in synchrony with Ensembl). Here, we provide an overview of Ensembl Genomes, with a focus on recent developments. These include the continued growth, more robust and reproducible sets of orthologues and paralogues, and enriched views of gene expression and gene function in plants. Finally, we report on our continued deeper integration with the Ensembl project, which forms a key part of our future strategy for dealing with the increasing quantity of available genome-scale data across the tree of life.

Expression Atlas update: from tissues to single cells

Expression Atlas is EMBL-EBI's resource for gene and protein expression. It sources and compiles data on the abundance and localisation of RNA and proteins in various biological systems and contexts and provides open access to this data for the research community. With the increased availability of single cell RNA-Seq datasets in the public archives, we have now extended Expression Atlas with a new added-value service to display gene expression in single cells. Single Cell Expression Atlas was launched in 2018 and currently includes 123 single cell RNA-Seq studies from 12 species. The website can be searched by genes within or across species to reveal experiments, tissues and cell types where this gene is expressed or under which conditions it is a marker gene. Within each study, cells can be visualized using a pre-calculated t-SNE plot and can be coloured by different features or by cell clusters based on gene expression. Within each experiment, there are links to downloadable files, such as RNA quantification matrices, clustering results, reports on protocols and associated metadata, such as assigned cell types.

New Drosophila Circadian Clock Mutants Affecting Temperature Compensation Induced by Targeted Mutagenesis of Timeless

Drosophila melanogaster has served as an excellent genetic model to decipher the molecular basis of the circadian clock. Two key proteins, PERIOD (PER) and TIMELESS (TIM), are particularly well explored and a number of various arrhythmic, slow, and fast clock mutants have been identified in classical genetic screens. Interestingly, the free running period (tau, τ) is influenced by temperature in some of these mutants, whereas τ is temperature-independent in other mutant lines as in wild-type flies. This, so-called "temperature compensation" ability is compromised in the mutant timeless allele "ritsu" (tim rit ), and, as we show here, also in the tim blind allele, mapping to the same region of TIM. To test if this region of TIM is indeed important for temperature compensation, we generated a collection of new mutants and mapped functional protein domains involved in the regulation of τ and in general clock function. We developed a protocol for targeted mutagenesis of specific gene regions utilizing the CRISPR/Cas9 technology, followed by behavioral screening. In this pilot study, we identified 20 new timeless mutant alleles with various impairments of temperature compensation. Molecular characterization revealed that the mutations included short in-frame insertions, deletions, or substitutions of a few amino acids resulting from the non-homologous end joining repair process. Our protocol is a fast and cost-efficient systematic approach for functional analysis of protein-coding genes and promoter analysis in vivo. Interestingly, several mutations with a strong temperature compensation defect map to one specific region of TIM. Although the exact mechanism of how these mutations affect TIM function is as yet unknown, our in silico analysis suggests they affect a putative nuclear export signal (NES) and phosphorylation sites of TIM. Immunostaining for PER was performed on two TIM mutants that display longer τ at 25°C and complete arrhythmicity at 28°C. Consistently with the behavioral phenotype, PER immunoreactivity was reduced in circadian clock neurons of flies exposed to elevated temperatures.

Role of co-repressor genomic landscapes in shaping the Notch response

Repressors are frequently deployed to limit the transcriptional response to signalling pathways. For example, several co-repressors interact directly with the DNA-binding protein CSL and are proposed to keep target genes silenced in the absence of Notch activity. However, the scope of their contributions remains unclear. To investigate co-repressor activity in the context of this well defined signalling pathway, we have analysed the genome-wide binding profile of the best-characterized CSL co-repressor in Drosophila, Hairless, and of a second CSL interacting repressor, SMRTER. As predicted there was significant overlap between Hairless and its CSL DNA-binding partner, both in Kc cells and in wing discs, where they were predominantly found in chromatin with active enhancer marks. However, while the Hairless complex was widely present at some Notch regulated enhancers in the wing disc, no binding was detected at others, indicating that it is not essential for silencing per se. Further analysis of target enhancers confirmed differential requirements for Hairless. SMRTER binding significantly overlapped with Hairless, rather than complementing it, and many enhancers were apparently co-bound by both factors. Our analysis indicates that the actions of Hairless and SMRTER gate enhancers to Notch activity and to Ecdysone signalling respectively, to ensure that the appropriate levels and timing of target gene expression are achieved.

The communication between cells that occurs during development, as well as in disease contexts, involves a small number of signalling pathways of which the Notch pathway is one. One outstanding question is how these pathways can bring about different gene responses in different contexts. As gene expression is co-ordinated by a mixture of activators and repressors, we set out to investigate whether the distribution of repressors across the genome is important in shaping whether genes are able to respond to Notch activity. Our results from analyzing the binding profile of two repressors, Hairless and SMRTER, show that, in many cases, they are not essential for preventing a gene from responding. Instead they are deployed at a limited number of genetic loci where they gate the response, helping to set a threshold for gene activation. Perturbations to their function lead to enhanced gene expression in limited territories rather than to new programmes of gene expression. Their main role therefore is to restrict the time or levels of signal that a gene needs to receive before it will respond.

Juvenile hormone signaling during reproduction and development of the linden bug, Pyrrhocoris apterus

Juvenile hormone (JH), a sesquiterpenoid produced by the insect corpus allatum gland (CA), prevents metamorphosis in larvae and stimulates vitellogenesis in adult females. Whether the same JH signaling pathway regulates both processes is presently unknown. Here, we employ the robust JH response during reproduction and development of the linden bug, Pyrrhocoris apterus, to compare the function of key JH-signaling genes encoding the JH receptor, Methoprene-tolerant (Met), its binding partner Taiman (Tai), and a JH-inducible protein, Krüppel-homolog 1 (Kr-h1). RNA interference (RNAi) with Met or Tai, but not Kr-h1, blocked ovarian development and suppressed vitellogenin gene expression in the fat body of females raised under reproduction-inducing conditions. Loss of Met and Tai matched the effects of CA ablation or the natural absence of JH during reproductive diapause. Stimulation of vitellogenesis by treatment of diapausing females with a JH mimic methoprene also required both Met and Tai in the fat body, whereas Kr-h1 RNAi had no effect. Therefore, the Met-Tai complex likely functions as a JH receptor during vitellogenesis. In contrast to Met and Kr-h1 that are both required for JH to prevent precocious metamorphosis in P. apterus larvae, removal of Tai disrupted larval ecdysis without causing premature adult development. Our results show that while Met operates during metamorphosis in larvae and reproduction in adult females, its partner Tai is only required for the latter. The diverse functions of JH thus likely rely on a common receptor whose actions are modulated by distinct components.