Mark Gerstein

Numerous statistical methods have emerged for inferring DNA motifs for transcription factors (TFs) from genomic regions. However, the process of selecting informative regions for motif inference remains understudied. Current approaches select regions with strong ChIP-seq signal for a given TF, assuming that such strong signal primarily results from specific interactions between the TF and its motif. Additionally, these selection approaches do not account for non-target motifs, i.e. motifs of other TFs; they presume the occurrence of these non-target motifs infrequent compared to that of the target motif, and thus assume these have minimal interference with the identification of the target. Leveraging extensive ChIP-seq datasets, we introduced the concept of TF signal ‘crowdedness’, referred to as C-score, for each genomic region. The C-score helps in highlighting TF signals arising from non-specific interactions. Moreover, by considering the C-score (and adjusting for the length of genomic regions), we can effectively mitigate interference of non-target motifs. Using these tools, we find that in many instances, strong ChIP-seq signal stems mainly from non-specific interactions, and the occurrence of non-target motifs significantly impacts the accurate inference of the target motif. Prioritizing genomic regions with reduced crowdedness and short length markedly improves motif inference. This ‘less-is-more’ effect suggests that ChIP-seq region selection warrants more attention.

Essential functions of mitogen-activated protein kinases (MAPKs) depend on their capacity to selectively phosphorylate a limited repertoire of substrates. MAPKs harbor a conserved groove located outside of the catalytic cleft that binds to short linear sequence motifs found in substrates and regulators. However, the weak and transient nature of these “docking” interactions poses a challenge to defining MAPK interactomes and associated sequence motifs. Here, we describe a yeast-based genetic screening pipeline to evaluate large collections of MAPK docking sequences in parallel. Using this platform, we analyzed a combinatorial library based on the docking sequences from the MAPK kinases MKK6 and MKK7, defining features critical for binding to the stress-activated MAPKs JNK1 and p38α. Our screen of a library consisting of ~12,000 sequences from the human proteome revealed multiple MAPK-selective interactors, including many that did not conform to previously defined docking motifs. Analysis of p38α/JNK1 exchange mutants identified specific docking groove residues that mediate selective binding. Last, we verified that docking sequences identified in the screen functioned in substrate recruitment in vitro and in cultured cells. Together, these studies establish an approach to characterize MAPK docking sequences and provide a resource for future investigation of signaling downstream of p38 and JNK.

Recent advances in single-cell sequencing technologies have provided unprecedented opportunities to measure the gene expression profile and RNA velocity of individual cells. However, modeling transcriptional dynamics is computationally challenging because of the high-dimensional, sparse nature of the single-cell gene expression measurements and the nonlinear regulatory relationships. Here, we present DeepVelo , a neural network–based ordinary differential equation that can model complex transcriptome dynamics by describing continuous-time gene expression changes within individual cells. We apply DeepVelo to public datasets from different sequencing platforms to (i) formulate transcriptome dynamics on different time scales, (ii) measure the instability of cell states, and (iii) identify developmental driver genes via perturbation analysis. Benchmarking against the state-of-the-art methods shows that DeepVelo can learn a more accurate representation of the velocity field. Furthermore, our perturbation studies reveal that single-cell dynamical systems could exhibit chaotic properties. In summary, DeepVelo allows data-driven discoveries of differential equations that delineate single-cell transcriptome dynamics.