Transfer learning for versatile and training free high content screening analyses.

High content screening (HCS) is a technology that automates cell biology experiments at large scale. A High Content Screen produces a high amount of microscopy images of cells under many conditions and requires that a dedicated image and data analysis workflow be designed for each assay to select hits. This heavy data analytic step remains challenging and has been recognized as one of the burdens hindering the adoption of HCS. In this work we propose a solution to hit selection by using transfer learning without additional training. A pretrained residual network is employed to encode each image of a screen into a discriminant representation. The deep features obtained are then corrected to account for well plate bias and misalignment. We then propose two training-free pipelines dedicated to the two main categories of HCS for compound selection: with or without positive control. When a positive control is available, it is used alongside the negative control to compute a linear discriminant axis, thus building a classifier without training. Once all samples are projected onto this axis, the conditions that best reproduce the positive control can be selected. When no positive control is available, the Mahalanobis distance is computed from each sample to the negative control distribution. The latter provides a metric to identify the conditions that alter the negative control's cell phenotype. This metric is subsequently used to categorize hits through a clustering step. Given the lack of available ground truth in HCS, we provide a qualitative comparison of the results obtained using this approach with results obtained with handcrafted image analysis features for compounds and siRNA screens with or without control. Our results suggests that the fully automated and generic pipeline we propose offers a good alternative to handcrafted dedicated image analysis approaches. Furthermore, we demonstrate that this solution select conditions of interest that had not been identified using the primary dedicated analysis. Altogether, this approach provides a fully automated, reproducible, versatile and comprehensive alternative analysis solution for HCS encompassing compound-based or downregulation screens, with or without positive controls, without the need for training or cell detection, or the development of a dedicated image analysis workflow.

Comprehensive mapping of exon junction complex binding sites reveals universal EJC deposition in Drosophila.

BACKGROUND: The exon junction complex (EJC) is involved in most steps of the mRNA life cycle, ranging from splicing to nonsense-mediated mRNA decay (NMD). It is assembled by the splicing machinery onto mRNA in a sequence-independent manner. A fundamental open question is whether the EJC is deposited onto all exon‒exon junctions or only on a subset of them. Several previous studies have made observations supportive of the latter, yet these have been limited by methodological constraints. RESULTS: In this study, we sought to overcome these limitations via the integration of two different approaches for transcriptome-wide mapping of EJCs. Our results revealed that nearly all, if not all, internal exons consistently harbor an EJC in Drosophila, demonstrating that EJC presence is an inherent consequence of the splicing reaction. Furthermore, our study underscores the limitations of eCLIP methods in fully elucidating the landscape of RBP binding sites. Our findings highlight how highly specific (low false positive) methodologies can lead to erroneous interpretations due to partial sensitivity (high false negatives). CONCLUSIONS: This study contributes to our understanding of EJC deposition and its association with pre-mRNA splicing. The universal presence of EJC on internal exons underscores its significance in ensuring proper mRNA processing. Additionally, our observations highlight the need to consider both specificity and sensitivity in RBP mapping methodologies.

Choroid plexuses carry nodal-like cilia that undergo axoneme regression from early adult stage.

Choroid plexuses (ChPs) produce cerebrospinal fluid and sense non-cell-autonomous stimuli to control the homeostasis of the central nervous system. They are mainly composed of epithelial multiciliated cells, whose development and function are still controversial. We have thus characterized the stepwise order of mammalian ChP epithelia cilia formation using a combination of super-resolution-microscopy approaches and mouse genetics. We show that ChP ciliated cells are built embryonically on a treadmill of spatiotemporally regulated events, starting with atypical centriole amplification and ending with the construction of nodal-like 9+0 cilia, characterized by both primary and motile features. ChP cilia undergo axoneme resorption at early postnatal stages through a microtubule destabilization process controlled by the microtubule-severing enzyme spastin and mitigated by polyglutamylation levels. Notably, this phenotype is preserved in humans, suggesting a conserved ciliary resorption mechanism in mammals.

Revealing invisible cell phenotypes with conditional generative modeling.

Biological sciences, drug discovery and medicine rely heavily on cell phenotype perturbation and microscope observation. However, most cellular phenotypic changes are subtle and thus hidden from us by natural cell variability: two cells in the same condition already look different. In this study, we show that conditional generative models can be used to transform an image of cells from any one condition to another, thus canceling cell variability. We visually and quantitatively validate that the principle of synthetic cell perturbation works on discernible cases. We then illustrate its effectiveness in displaying otherwise invisible cell phenotypes triggered by blood cells under parasite infection, or by the presence of a disease-causing pathological mutation in differentiated neurons derived from iPSCs, or by low concentration drug treatments. The proposed approach, easy to use and robust, opens the door to more accessible discovery of biological and disease biomarkers.

The model diatom Phaeodactylum tricornutum provides insights into the diversity and function of microeukaryotic DNA methyltransferases.

Cytosine methylation is an important epigenetic mark involved in the transcriptional control of transposable elements in mammals, plants and fungi. The Stramenopiles-Alveolate-Rhizaria (SAR) lineages are a major group of ecologically important marine microeukaryotes, including the phytoplankton groups diatoms and dinoflagellates. However, little is known about their DNA methyltransferase diversity. Here, we performed an in-silico analysis of DNA methyltransferases found in marine microeukaryotes and showed that they encode divergent DNMT3, DNMT4, DNMT5 and DNMT6 enzymes. Furthermore, we found three classes of enzymes within the DNMT5 family. Using a CRISPR/Cas9 strategy we demonstrated that the loss of the DNMT5a gene correlates with a global depletion of DNA methylation and overexpression of young transposable elements in the model diatom Phaeodactylum tricornutum. The study provides a view of the structure and function of a DNMT family in the SAR supergroup using an attractive model species.

Evolution is not Uniform Along Coding Sequences.

Amino acids evolve at different speeds within protein sequences, because their functional and structural roles are different. Notably, amino acids located at the surface of proteins are known to evolve more rapidly than those in the core. In particular, amino acids at the N- and C-termini of protein sequences are likely to be more exposed than those at the core of the folded protein due to their location in the peptidic chain, and they are known to be less structured. Because of these reasons, we would expect that amino acids located at protein termini would evolve faster than residues located inside the chain. Here we test this hypothesis and found that amino acids evolve almost twice as fast at protein termini compared with those in the center, hinting at a strong topological bias along the sequence length. We further show that the distribution of solvent-accessible residues and functional domains in proteins readily explain how structural and functional constraints are weaker at their termini, leading to the observed excess of amino acid substitutions. Finally, we show that the specific evolutionary rates at protein termini may have direct consequences, notably misleading in silico methods used to infer sites under positive selection within genes. These results suggest that accounting for positional information should improve evolutionary models.

Cell painting transfer increases screening hit rate.

Drug discovery uses high throughput screening to identify compounds that interact with a molecular target or that alter a phenotype favorably. The cautious selection of molecules used for such a screening is instrumental and is tightly related to the hit rate. In this work, we wondered if cell painting, a general-purpose image-based assay, could be used as an efficient proxy for compound selection, thus increasing the success rate of a specific assay. To this end, we considered cell painting images with 30,000 molecules treatments, and selected compounds that produced a visual effect close to the positive control of an assay, by using the Frechet Inception Distance. We then compared the hit rates of such a preselection with what was actually obtained in real screening campaigns. As a result, cell painting would have permitted a significant increase in the success rate and, even for one of the assays, would have allowed to reach 80% of the hits with 10 times fewer compounds to test. We conclude that images of a cell painting assay can be directly used for compound selection prior to screening, and we provide a simple quantitative approach in order to do so.

A medullary centre for lapping in mice.

It has long been known that orofacial movements for feeding can be triggered, coordinated, and often rhythmically organized at the level of the brainstem, without input from higher centers. We uncover two nuclei that can organize the movements for ingesting fluids in mice. These neuronal groups, IRtPhox2b and Peri5Atoh1, are marked by expression of the pan-autonomic homeobox gene Phox2b and are located, respectively, in the intermediate reticular formation of the medulla and around the motor nucleus of the trigeminal nerve. They are premotor to all jaw-opening and tongue muscles. Stimulation of either, in awake animals, opens the jaw, while IRtPhox2b alone also protracts the tongue. Moreover, stationary stimulation of IRtPhox2b entrains a rhythmic alternation of tongue protraction and retraction, synchronized with jaw opening and closing, that mimics lapping. Finally, fiber photometric recordings show that IRtPhox2b is active during volitional lapping. Our study identifies one of the subcortical nuclei underpinning a stereotyped feeding behavior.

Genome wide natural variation of H3K27me3 selectively marks genes predicted to be important for cell differentiation in Phaeodactylum tricornutum.

In multicellular organisms, Polycomb Repressive Complex2 (PRC2) is known to deposit tri-methylation of lysine 27 of histone H3 (H3K27me3) to establish and maintain gene silencing, critical for developmentally regulated processes. The PRC2 complex is absent in both widely studied model yeasts, which initially suggested that PRC2 arose with the emergence of multicellularity. However, its discovery in several unicellular species including microalgae questions its role in unicellular eukaryotes. Here, we use Phaeodactylum tricornutum enhancer of zeste E(z) knockouts and show that P. tricornutum E(z) is responsible for di- and tri-methylation of lysine 27 of histone H3. H3K27me3 depletion abolishes cell morphology in P. tricornutum providing evidence for its role in cell differentiation. Genome-wide profiling of H3K27me3 in fusiform and triradiate cells further revealed genes that may specify cell identity. These results suggest a role for PRC2 and its associated mark in cell differentiation in unicellular species, and highlight their ancestral function in a broader evolutionary context than currently is appreciated.

PhaeoNet: A Holistic RNAseq-Based Portrait of Transcriptional Coordination in the Model Diatom Phaeodactylum tricornutum.

Transcriptional coordination is a fundamental component of prokaryotic and eukaryotic cell biology, underpinning the cell cycle, physiological transitions, and facilitating holistic responses to environmental stress, but its overall dynamics in eukaryotic algae remain poorly understood. Better understanding of transcriptional partitioning may provide key insights into the primary metabolism pathways of eukaryotic algae, which frequently depend on intricate metabolic associations between the chloroplasts and mitochondria that are not found in plants. Here, we exploit 187 publically available RNAseq datasets generated under varying nitrogen, iron and phosphate growth conditions to understand the co-regulatory principles underpinning transcription in the model diatom Phaeodactylum tricornutum. Using WGCNA (Weighted Gene Correlation Network Analysis), we identify 28 merged modules of co-expressed genes in the P. tricornutum genome, which show high connectivity and correlate well with previous microarray-based surveys of gene co-regulation in this species. We use combined functional, subcellular localization and evolutionary annotations to reveal the fundamental principles underpinning the transcriptional co-regulation of genes implicated in P. tricornutum chloroplast and mitochondrial metabolism, as well as the functions of diverse transcription factors underpinning this co-regulation. The resource is publically available as PhaeoNet, an advanced tool to understand diatom gene co-regulation.

Unraveling spatial cellular pattern by computational tissue shuffling.

Cell biology relies largely on reproducible visual observations. Unlike cell culture, tissues are heterogeneous, making difficult the collection of biological replicates that would spotlight a precise location. In consequence, there is no standard approach for estimating the statistical significance of an observed pattern in a tissue sample. Here, we introduce SET (for Synthesis of Epithelial Tissue), a method that can accurately reconstruct the cell tessellation formed by an epithelium in a microscopy image as well as thousands of alternative synthetic tessellations made of the exact same cells. SET can build an accurate null distribution to statistically test if any local pattern is necessarily the result of a process, or if it could be explained by chance in the given context. We provide examples in various tissues where visible, and invisible, cell and subcellular patterns are unraveled in a statistically significant manner using a single image and without any parameter settings.

Coordination of transcriptional and translational regulations in human epithelial cells infected by Listeria monocytogenes.

The invasion of mammalian cells by intracellular bacterial pathogens reshuffles their gene expression and functions; however, we lack dynamic insight into the distinct control levels that shape the host response. Here, we have addressed the respective contribution of transcriptional and translational regulations during a time-course of infection of human intestinal epithelial cells by an epidemic strain of Listeria monocytogenes, using transcriptome analysis paralleled with ribosome profiling. Upregulations were dominated by early transcriptional activation of pro-inflammatory genes, whereas translation inhibition appeared as the major driver of downregulations. Instead of a widespread but transient shutoff, translation inhibition affected specifically and durably transcripts encoding components of the translation machinery harbouring a 5'-terminal oligopyrimidine motif. Pre-silencing the most repressed target gene (PABPC1) slowed down the intracellular multiplication of Listeria monocytogenes, suggesting that the infected host cell can benefit from the repression of genes involved in protein synthesis and thereby better control infection.

In vivo large-scale analysis of Drosophila neuronal calcium traces by automated tracking of single somata.

How does the concerted activity of neuronal populations shape behavior? Impediments to address this question are primarily due to critical experimental barriers. An integrated perspective on large scale neural information processing requires an in vivo approach that can combine the advantages of exhaustively observing all neurons dedicated to a given type of stimulus, and simultaneously achieve a resolution that is precise enough to capture individual neuron activity. Current experimental data from in vivo observations are either restricted to a small fraction of the total number of neurons, or are based on larger brain volumes but at a low spatial and temporal resolution. Consequently, fundamental questions as to how sensory information is represented on a population scale remain unanswered. In Drosophila melanogaster, the mushroom body (MB) represents an excellent model to analyze sensory coding and memory plasticity. In this work, we present an experimental setup coupled with a dedicated computational method that provides in vivo measurements of the activity of hundreds of densely packed somata uniformly spread in the MB. We exploit spinning-disk confocal 3D imaging over time of the whole MB cell body layer in vivo while it is exposed to olfactory stimulation. Importantly, to derive individual signal from densely packed somata, we have developed a fully automated image analysis procedure that takes advantage of the specificities of our data. After anisotropy correction, our approach operates a dedicated spot detection and registration over the entire time sequence to transform trajectories to identifiable clusters. This enabled us to discard spurious detections and reconstruct missing ones in a robust way. We demonstrate that this approach outperformed existing methods in this specific context and made possible high-throughput analysis of approximately 500 single somata uniformly spread over the MB in various conditions. Applying this approach, we find that learned experiences change the population code of odor representations in the MB. After long-term memory (LTM) formation, we quantified an increase in responsive somata count and a stable single neuron signal. We predict that this method, which should further enable studying the population pattern of neuronal activity, has the potential to uncover fine details of sensory processing and memory plasticity.

Artificially decreasing cortical tension generates aneuploidy in mouse oocytes.

Human and mouse oocytes' developmental potential can be predicted by their mechanical properties. Their development into blastocysts requires a specific stiffness window. In this study, we combine live-cell and computational imaging, laser ablation, and biophysical measurements to investigate how deregulation of cortex tension in the oocyte contributes to early developmental failure. We focus on extra-soft cells, the most common defect in a natural population. Using two independent tools to artificially decrease cortical tension, we show that chromosome alignment is impaired in extra-soft mouse oocytes, despite normal spindle morphogenesis and dynamics, inducing aneuploidy. The main cause is a cytoplasmic increase in myosin-II activity that could sterically hinder chromosome capture. We describe here an original mode of generation of aneuploidies that could be very common in oocytes and could contribute to the high aneuploidy rate observed during female meiosis, a leading cause of infertility and congenital disorders.

PySpacell: A Python Package for Spatial Analysis of Cell Images.

Technologies such as microscopy, sequential hybridization, and mass spectrometry enable quantitative single-cell phenotypic and molecular measurements in situ. Deciphering spatial phenotypic and molecular effects on the single-cell level is one of the grand challenges and a key to understanding the effects of cell-cell interactions and microenvironment. However, spatial information is usually overlooked by downstream data analyses, which usually consider single-cell read-out values as independent measurements for further averaging or clustering, thus disregarding spatial locations. With this work, we attempt to fill this gap. We developed a toolbox that allows one to test for the presence of a spatial effect in microscopy images of adherent cells and estimate the spatial scale of this effect. The proposed Python module can be used for any light microscopy images of cells as well as other types of single-cell data such as in situ transcriptomics or metabolomics. The input format of our package matches standard output formats from image analysis tools such as CellProfiler, Fiji, or Icy and thus makes our toolbox easy and straightforward to use, yet offering a powerful statistical approach for a wide range of applications. © 2019 International Society for Advancement of Cytometry.

Active Fluctuations of the Nuclear Envelope Shape the Transcriptional Dynamics in Oocytes.

Nucleus position in cells can act as a developmental cue. Mammalian oocytes position their nucleus centrally using an F-actin-mediated pressure gradient. The biological significance of nucleus centering in mammalian oocytes being unknown, we sought to assess the F-actin pressure gradient effect on the nucleus. We addressed this using a dedicated computational 3D imaging approach, biophysical analyses, and a nucleus repositioning assay in mouse oocytes mutant for cytoplasmic F-actin. We found that the cytoplasmic activity, in charge of nucleus centering, shaped the nucleus while promoting nuclear envelope fluctuations and chromatin motion. Off-centered nuclei in F-actin mutant oocytes were misshaped with immobile chromatin and modulated gene expression. Restoration of F-actin in mutant oocytes rescued nucleus architecture fully and gene expression partially. Thus, the F-actin-mediated pressure gradient also modulates nucleus dynamics in oocytes. Moreover, this study supports a mechano-transduction model whereby cytoplasmic microfilaments could modulate oocyte transcriptome, essential for subsequent embryo development.

m5C Methylation Guides Systemic Transport of Messenger RNA over Graft Junctions in Plants.

In plants, transcripts move to distant body parts to potentially act as systemic signals regulating development and growth. Thousands of messenger RNAs (mRNAs) are transported across graft junctions via the phloem to distinct plant parts. Little is known regarding features, structural motifs, and potential base modifications of transported transcripts and how these may affect their mobility. We identified Arabidopsis thaliana mRNAs harboring the modified base 5-methylcytosine (m5C) and found that these are significantly enriched in mRNAs previously described as mobile, moving over graft junctions to distinct plant parts. We confirm this finding with graft-mobile methylated mRNAs TRANSLATIONALLY CONTROLLED TUMOR PROTEIN 1 (TCTP1) and HEAT SHOCK COGNATE PROTEIN 70.1 (HSC70.1), whose mRNA transport is diminished in mutants deficient in m5C mRNA methylation. Together, our results point toward an essential role of cytosine methylation in systemic mRNA mobility in plants and that TCTP1 mRNA mobility is required for its signaling function.

ALFA: annotation landscape for aligned reads.

BACKGROUND: The last 10 years have seen the rise of countless functional genomics studies based on Next-Generation Sequencing (NGS). In the vast majority of cases, whatever the species, whatever the experiment, the two first steps of data analysis consist of a quality control of the raw reads followed by a mapping of those reads to a reference genome/transcriptome. Subsequent steps then depend on the type of study that is being made. While some tools have been proposed for investigating data quality after the mapping step, there is no commonly adopted framework that would be easy to use and broadly applicable to any NGS data type. RESULTS: We present ALFA, a simple but universal tool that can be used after the mapping step on any kind of NGS experiment data for any organism with available genomic annotations. In a single command line, ALFA can compute and display distribution of reads by categories (exon, intron, UTR, etc.) and biotypes (protein coding, miRNA, etc.) for a given aligned dataset with nucleotide precision. We present applications of ALFA to Ribo-Seq and RNA-Seq on Homo sapiens, CLIP-Seq on Mus musculus, RNA-Seq on Saccharomyces cerevisiae, Bisulfite sequencing on Arabidopsis thaliana and ChIP-Seq on Caenorhabditis elegans. CONCLUSIONS: We show that ALFA provides a powerful and broadly applicable approach for post mapping quality control and to produce a global overview using common or dedicated annotations. It is made available to the community as an easy to install command line tool and from the Galaxy Tool Shed.

Adult Neural Stem Cells and Multiciliated Ependymal Cells Share a Common Lineage Regulated by the Geminin Family Members.

Adult neural stem cells and multiciliated ependymal cells are glial cells essential for neurological functions. Together, they make up the adult neurogenic niche. Using both high-throughput clonal analysis and single-cell resolution of progenitor division patterns and fate, we show that these two components of the neurogenic niche are lineally related: adult neural stem cells are sister cells to ependymal cells, whereas most ependymal cells arise from the terminal symmetric divisions of the lineage. Unexpectedly, we found that the antagonist regulators of DNA replication, GemC1 and Geminin, can tune the proportion of neural stem cells and ependymal cells. Our findings reveal the controlled dynamic of the neurogenic niche ontogeny and identify the Geminin family members as key regulators of the initial pool of adult neural stem cells.

Monitored eCLIP: high accuracy mapping of RNA-protein interactions.

CLIP-seq methods provide transcriptome-wide snapshots of RNA-protein interactions in live cells. Reverse transcriptases stopping at cross-linked nucleotides sign for RNA-protein binding sites. Reading through cross-linked positions results in false binding site assignments. In the 'monitored enhanced CLIP' (meCLIP) method, a barcoded biotinylated linker is ligated at the 5' end of cross-linked RNA fragments to purify RNA prior to the reverse transcription. cDNAs keeping the barcode sequence correspond to reverse transcription read-throughs. Read through occurs in unpredictable proportions, representing up to one fourth of total reads. Filtering out those reads strongly improves reliability and precision in protein binding site assignment.