Cell cycle gene regulation dynamics revealed by RNA velocity and deep-learning.

Despite the fact that the cell cycle is a fundamental process of life, a detailed quantitative understanding of gene regulation dynamics throughout the cell cycle is far from complete. Single-cell RNA-sequencing (scRNA-seq) technology gives access to these dynamics without externally perturbing the cell. Here, by generating scRNA-seq libraries in different cell systems, we observe cycling patterns in the unspliced-spliced RNA space of cell cycle-related genes. Since existing methods to analyze scRNA-seq are not efficient to measure cycling gene dynamics, we propose a deep learning approach (DeepCycle) to fit these patterns and build a high-resolution map of the entire cell cycle transcriptome. Characterizing the cell cycle in embryonic and somatic cells, we identify major waves of transcription during the G1 phase and systematically study the stages of the cell cycle. Our work will facilitate the study of the cell cycle in multiple cellular models and different biological contexts.

Regulation of transcription reactivation dynamics exiting mitosis.

Proliferating cells experience a global reduction of transcription during mitosis, yet their cell identity is maintained and regulatory information is propagated from mother to daughter cells. Mitotic bookmarking by transcription factors has been proposed as a potential mechanism to ensure the reactivation of transcription at the proper set of genes exiting mitosis. Recently, mitotic transcription and waves of transcription reactivation have been observed in synchronized populations of human hepatoma cells. However, the study did not consider that mitotic-arrested cell populations progressively desynchronize leading to measurements of gene expression on a mixture of cells at different internal cell-cycle times. Moreover, it is not well understood yet what is the precise role of mitotic bookmarking on mitotic transcription as well as on the transcription reactivation waves. Ultimately, the core gene regulatory network driving the precise transcription reactivation dynamics remains to be identified. To address these questions, we developed a mathematical model to correct for the progressive desynchronization of cells and estimate gene expression dynamics with respect to a cell-cycle pseudotime. Furthermore, we used a multiple linear regression model to infer transcription factor activity dynamics. Our analysis allows us to characterize waves of transcription factor activities exiting mitosis and predict a core gene regulatory network responsible of the transcription reactivation dynamics. Moreover, we identified more than 60 transcription factors that are highly active during mitosis and represent new candidates of mitotic bookmarking factors which could be relevant therapeutic targets to control cell proliferation.

A Model-Driven Quantitative Analysis of Retrotransposon Distributions in the Human Genome.

Retrotransposons, DNA sequences capable of creating copies of themselves, compose about half of the human genome and played a central role in the evolution of mammals. Their current position in the host genome is the result of the retrotranscription process and of the following host genome evolution. We apply a model from statistical physics to show that the genomic distribution of the two most populated classes of retrotransposons in human deviates from random placement, and that this deviation increases with time. The time dependence suggests a major role of the host genome dynamics in shaping the current retrotransposon distributions. Focusing on a neutral scenario, we show that a simple model based on random placement followed by genome expansion and sequence duplications can reproduce the empirical retrotransposon distributions, even though more complex and possibly selective mechanisms can have contributed. Besides the inherent interest in understanding the origin of current retrotransposon distributions, this work sets a general analytical framework to analyze quantitatively the effects of genome evolutionary dynamics on the distribution of genomic elements.

Temporal specificity and heterogeneity of Drosophila immune cells.

Immune cells provide defense against non-self and have recently been shown to also play key roles in diverse processes such as development, metabolism, and tumor progression. The heterogeneity of Drosophila immune cells (hemocytes) remains an open question. Using bulk RNA sequencing, we find that the hemocytes display distinct features in the embryo, a closed and rapidly developing system, compared to the larva, which is exposed to environmental and metabolic challenges. Through single-cell RNA sequencing, we identify fourteen hemocyte clusters present in unchallenged larvae and associated with distinct processes, e.g., proliferation, phagocytosis, metabolic homeostasis, and humoral response. Finally, we characterize the changes occurring in the hemocyte clusters upon wasp infestation, which triggers the differentiation of a novel hemocyte type, the lamellocyte. This first molecular atlas of hemocytes provides insights and paves the way to study the biology of the Drosophila immune cells in physiological and pathological conditions.

A rare codon-based translational program of cell proliferation.

BACKGROUND: The speed of translation elongation is primarily determined by the abundance of tRNAs. Thus, the codon usage influences the rate with which individual mRNAs are translated. As the nature of tRNA pools and modifications can vary across biological conditions, codon elongation rates may also vary, leading to fluctuations in the protein production from individual mRNAs. Although it has been observed that functionally related mRNAs exhibit similar codon usage, presumably to provide an effective way to coordinate expression of multiple proteins, experimental evidence for codon-mediated translation efficiency modulation of functionally related mRNAs in specific conditions is scarce and the associated mechanisms are still debated. RESULTS: Here, we reveal that mRNAs whose expression increases during cell proliferation are enriched in rare codons, poorly adapted to tRNA pools. Ribosome occupancy profiling and proteomics measurements show that upon increased cell proliferation, transcripts enriched in rare codons undergo a higher translation boost than transcripts with common codons. Re-coding of a fluorescent reporter with rare codons increased protein output by ~ 30% relative to a reporter re-coded with common codons. Although the translation capacity of proliferating cells was higher compared to resting cells, we did not find evidence for the regulation of individual tRNAs. Among the models that were proposed so far to account for codon-mediated translational regulation upon changing conditions, the one that seems most consistent with our data involves a global upregulation of ready-to-translate tRNAs, which we show can lead to a higher increase in the elongation velocity at rare codons compared to common codons. CONCLUSIONS: We propose that the alleviation of translation bottlenecks in rapidly dividing cells enables preferential upregulation of pro-proliferation proteins, encoded by mRNAs that are enriched in rare codons.

Protein synthesis rates and ribosome occupancies reveal determinants of translation elongation rates.

Although protein synthesis dynamics has been studied both with theoretical models and by profiling ribosome footprints, the determinants of ribosome flux along open reading frames (ORFs) are not fully understood. Combining measurements of protein synthesis rate with ribosome footprinting data, we here inferred translation initiation and elongation rates for over a 1,000 ORFs in exponentially growing wild-type yeast cells. We found that the amino acid composition of synthesized proteins is as important a determinant of translation elongation rate as parameters related to codon and transfer RNA (tRNA) adaptation. We did not find evidence of ribosome collisions curbing the protein output of yeast transcripts, either in high translation conditions associated with exponential growth, or in strains in which deletion of individual ribosomal protein (RP) genes leads to globally increased or decreased translation. Slow translation elongation is characteristic of RP-encoding transcripts, which have markedly lower protein output compared with other transcripts with equally high ribosome densities.

Terminal exon characterization with TECtool reveals an abundance of cell-specific isoforms.

Sequencing of RNA 3' ends has uncovered numerous sites that do not correspond to the termination sites of known transcripts. Through their 3' untranslated regions, protein-coding RNAs interact with RNA-binding proteins and microRNAs, which regulate many properties, including RNA stability and subcellular localization. We developed the terminal exon characterization (TEC) tool ( http://tectool.unibas.ch ), which can be used with RNA-sequencing data from any species for which a genome annotation that includes sites of RNA cleavage and polyadenylation is available. We discovered hundreds of previously unknown isoforms and cell-type-specific terminal exons in human cells. Ribosome profiling data revealed that many of these isoforms were translated. By applying TECtool to single-cell sequencing data, we found that the newly identified isoforms were expressed in subpopulations of cells. Thus, TECtool enables the identification of previously unknown isoforms in well-studied cell systems and in rare cell types.

Explicit Modeling of siRNA-Dependent On- and Off-Target Repression Improves the Interpretation of Screening Results.

RNAi is broadly used to map gene regulatory networks, but the identification of genes that are responsible for the observed phenotypes is challenging, as small interfering RNAs (siRNAs) simultaneously downregulate the intended on targets and many partially complementary off targets. Additionally, the scarcity of publicly available control datasets hinders the development and comparative evaluation of computational methods for analyzing the data. Here, we introduce PheLiM (https://github.com/andreariba/PheLiM), a method that uses predictions of siRNA on- and off-target downregulation to infer gene-specific contributions to phenotypes. To assess the performance of PheLiM, we carried out siRNA- and CRISPR/Cas9-based genome-wide screening of two well-characterized pathways, bone morphogenetic protein (BMP) and nuclear factor κB (NF-κB), and we reanalyzed publicly available siRNA screens. We demonstrate that PheLiM has the overall highest accuracy and most reproducible results compared to other available methods. PheLiM can accommodate various methods for predicting siRNA off targets and is broadly applicable to the identification of genes underlying complex phenotypes.

Hybridization of powdery mildew strains gives rise to pathogens on novel agricultural crop species.

Throughout the history of agriculture, many new crop species (polyploids or artificial hybrids) have been introduced to diversify products or to increase yield. However, little is known about how these new crops influence the evolution of new pathogens and diseases. Triticale is an artificial hybrid of wheat and rye, and it was resistant to the fungal pathogen powdery mildew (Blumeria graminis) until 2001 (refs. 1,2,3). We sequenced and compared the genomes of 46 powdery mildew isolates covering several formae speciales. We found that B. graminis f. sp. triticale, which grows on triticale and wheat, is a hybrid between wheat powdery mildew (B. graminis f. sp. tritici) and mildew specialized on rye (B. graminis f. sp. secalis). Our data show that the hybrid of the two mildews specialized on two different hosts can infect the hybrid plant species originating from those two hosts. We conclude that hybridization between mildews specialized on different species is a mechanism of adaptation to new crops introduced by agriculture.

Interplay of microRNA and epigenetic regulation in the human regulatory network.

The expression of protein-coding genes is controlled by a complex network of regulatory interactions. It is becoming increasingly appreciated that post-transcriptional repression by microRNAs, a class of small non-coding RNAs, is a key layer of regulation in several biological processes. In this contribution, we discuss the interplay between microRNAs and epigenetic regulators. Among the mixed genetic circuits composed by these two different kinds of regulation, it seems that a central role is played by double-negative feedback loops in which a microRNA inhibits an epigenetic regulator and in turn is controlled at the epigenetic level by the same regulator. We discuss a few relevant properties of this class of network motifs and their potential role in cell differentiation. In particular, using mathematical modeling we show how this particular circuit can exhibit a switch-like behavior between two alternative steady states, while being robust to stochastic transitions between these two states, a feature presumably required for circuits involved in cell fate decision. Finally, we present a list of putative double-negative feedback loops from a literature survey combined with bioinformatic analysis, and discuss in detail a few examples.

A combination of transcriptional and microRNA regulation improves the stability of the relative concentrations of target genes.

It is well known that, under suitable conditions, microRNAs are able to fine tune the relative concentration of their targets to any desired value. We show that this function is particularly effective when one of the targets is a Transcription Factor (TF) which regulates the other targets. This combination defines a new class of feed-forward loops (FFLs) in which the microRNA plays the role of master regulator. Using both deterministic and stochastic equations, we show that these FFLs are indeed able not only to fine-tune the TF/target ratio to any desired value as a function of the miRNA concentration but also, thanks to the peculiar topology of the circuit, to ensure the stability of this ratio against stochastic fluctuations. These two effects are due to the interplay between the direct transcriptional regulation and the indirect TF/Target interaction due to competition of TF and target for miRNA binding (the so called "sponge effect"). We then perform a genome wide search of these FFLs in the human regulatory network and show that they are characterized by a very peculiar enrichment pattern. In particular, they are strongly enriched in all the situations in which the TF and its target have to be precisely kept at the same concentration notwithstanding the environmental noise. As an example we discuss the FFL involving E2F1 as Transcription Factor, RB1 as target and miR-17 family as master regulator. These FFLs ensure a tight control of the E2F/RB ratio which in turns ensures the stability of the transition from the G0/G1 to the S phase in quiescent cells.