Towards inferring causal gene regulatory networks from single cell expression Measurements, bioRxiv, 2018-09-27

AbstractSingle-cell transcriptome sequencing now routinely samples thousands of cells, potentially providing enough data to reconstruct causal gene regulatory networks from observational data. Here, we present Scribe, a toolkit for detecting and visualizing causal regulatory interactions between genes and explore the potential for single-cell experiments to power network reconstruction. Scribe employs Restricted Directed Information to determine causality by estimating the strength of information transferred from a potential regulator to its downstream target. We apply Scribe and other leading approaches for causal network reconstruction to several types of single-cell measurements and show that there is a dramatic drop in performance for pseudotime” ordered single-cell data compared to true time series data. We demonstrate that performing causal inference requires temporal coupling between measurements. We show that methods such as “RNA velocity” restore some degree of coupling through an analysis of chromaffin cell fate commitment. These analyses therefore highlight an important shortcoming in experimental and computational methods for analyzing gene regulation at single-cell resolution and point the way towards overcoming it.

biorxiv genomics 100-200-users 2018

beditor A computational workflow for designing libraries of guide RNAs for CRISPR-mediated base editing, bioRxiv, 2018-09-26

ABSTRACTCRISPR-mediated base editors have opened unique avenues for scar-free genome-wide mutagenesis. Here, we describe a comprehensive computational workflow called beditor that can be broadly adapted for designing guide RNA libraries with a range of CRISPR-mediated base editors, PAM recognition sequences and genomes of many species. Additionally, in order to assist users in selecting the best sets of guide RNAs for their experiments, a priori estimates, called beditor scores are calculated. These beditor scores are intended to select guide RNAs that conform to requirements for optimal base editing the editable base falls within maximum activity window of the CRISPR-mediated base editor and produces non-confounding mutational effects with minimal predicted off-target effects. We demonstrate the utility of the software by designing guide RNAs for base-editing to create or remove thousands of clinically important human disease mutations.

biorxiv synthetic-biology 0-100-users 2018

Efficient generation of endogenous fluorescent reporters by Nested CRISPR in Caenorhabditis elegans, bioRxiv, 2018-09-26

AbstractCRISPR-based genome editing methods in model organisms are evolving at an extraordinary speed. Whereas the generation of deletion or missense mutants is quite straightforward, the production of endogenous fluorescent reporters is still inefficient. The use of plasmids with selection markers is an effective methodology, but often requires laborious and complicated cloning steps. We have established a cloning-free ribonucleoprotein-driven Nested CRISPR method that robustly produces endogenous fluorescent reporters. This methodology is based on the division of the GFP and mCherry sequences in three fragments. In the first step we use ssDNA donors (≤200 bp) to insert 5’ and 3’ fragments in the place of interest. In the second step, we use these sequences as homology regions for Homology Directed Repair (HDR) with a dsDNA donor (PCR product, ≈700 bp) including the middle fragment, thus completing the fluorescent protein sequence. This method is advantageous because the first step with ssDNA donors is known to be very efficient, and the second step, uses universal reagents, including validated PCR products and crRNAs, to create fluorescent reporters reaching reliable editing efficiencies as high as 40%. We have also used Nested CRISPR in a non-essential gene to produce a deletion mutant in the first step and a transcriptional reporter in the second step.In the search of modifications to optimize the method, we tested synthetic sgRNAs, but we did not observe a significant increase in the efficacy compared to independently adding tracrRNA and crRNA to the injection mix. Conveniently, we also found that both steps of Nested CRISPR could be performed in a single injection. Finally, we discuss the utility of Nested CRISPR for targeted insertion of long DNA fragments in other systems and prospects of this method in the future.

biorxiv genetics 0-100-users 2018

Exploiting selection at linked sites to infer the rate and strength of adaptation, bioRxiv, 2018-09-26

Genomic data encodes past evolutionary events and has the potential to reveal the strength, rate, and biological drivers of adaptation. However, robust estimation of adaptation rate (alpha) and adaptation strength remains a challenging problem because evolutionary processes such as demography, linkage, and non-neutral polymorphism can confound inference. Here, we exploit the influence of background selection to reduce the fixation rate of weakly-beneficial alleles to jointly infer the strength and rate of adaptation. We develop a novel MK-based method (ABC-MK) to infer adaptation rate and strength, and estimate alpha=0.135 in human protein-coding sequences, 72% of which is contributed by weakly adaptive variants. We show that in this adaptation regime alpha is reduced approx. 25% by linkage genome-wide. Moreover, we show that virus-interacting proteins (VIPs) undergo adaptation that is both stronger and nearly twice as frequent as the genome average (alpha=0.224, 56% due to strongly-beneficial alleles). Our results suggest that while most adaptation in human proteins is weakly-beneficial, adaptation to viruses is often strongly-beneficial. Our method provides a robust framework for estimating adaptation rate and strength across species.

biorxiv evolutionary-biology 0-100-users 2018

Slow transcriptional elongation causes embryonic lethality and perturbs kinetic coupling of long neural genes, bioRxiv, 2018-09-26

The rate of RNA Polymerase II (RNAPII) elongation has an important role in the control of Alternative splicing (AS); however, the in vivo consequences of an altered elongation rate are unknown. Here, we generated mouse embryonic stem cells (ESCs) knocked-in for a slow elongating form of RNAPII. We show that a reduced transcriptional elongation rate results in early embryonic lethality in mice and impairs the differentiation of ESCs into the neural lineage. This is accompanied by changes in splicing and in gene expression in ESCs and along the pathway of neuronal differentiation. In particular, we found a crucial role for RNAPII elongation rate in transcription and splicing of long neuronal genes involved in synapse signaling. The impact of the kinetic coupling of RNAPII elongation rate with AS is more predominant in ESC-differentiated neurons than in pluripotent cells. Our results demonstrate the requirement for an appropriate transcriptional elongation rate to ensure proper gene expression and to regulate AS during development.

biorxiv molecular-biology 0-100-users 2018

Stronger and higher proportion of beneficial amino acid changing mutations in humans compared to mice and flies, bioRxiv, 2018-09-26

ABSTRACTQuantifying and comparing the amount of adaptive evolution among different species is key to understanding evolutionary processes. Previous studies have shown differences in adaptive evolution across species, however their specific causes remain elusive. Here, we use improved modeling of weakly deleterious mutations and the demographic history of the outgroup species and estimate that 30–34% of nonsynonymous substitutions between humans and outgroup species have been fixed by positive selection. This estimate is much higher than previous estimates, which did not account for the population size of the outgroup species. Next, we directly estimate the proportion and selection coefficients of newly arising strongly beneficial nonsynonymous mutations in humans, mice, and D. melanogaster by examining patterns of polymorphism and divergence. We develop a novel composite likelihood framework to test whether these parameters differ across species. Overall, we reject a model with the same proportion and the same selection coefficients of beneficial mutations across species, and estimate that humans have a higher proportion of beneficial mutations compared to Drosophila and mice. We demonstrate that this result cannot be attributed to biased gene conversion. In summary, we find the proportion of beneficial mutations is higher in humans than in D. melanogaster or mice, suggesting that organismal complexity, which increases the number of steps required in adaptive walks, may be a key predictor of the amount of adaptive evolution within a species.

biorxiv evolutionary-biology 0-100-users 2018