Cis-regulatory code for determining the action of Foxd as both an activator and a repressor in ascidian embryos

Dynamics of Chromatin Opening across Larval Development in the Urochordate Ascidian Ciona savignyi

Ascidian larvae undergo tail elongation and notochord lumenogenesis, making them an ideal model for investigating tissue morphogenesis in embryogenesis. The cellular and mechanical mechanisms of these processes have been studied; however, the underlying molecular regulatory mechanism remains to be elucidated. In this study, assays for transposase-accessible chromatin using sequencing (ATAC-seq) and RNA sequencing (RNA-seq) were applied to investigate potential regulators of the development of ascidian Ciona savignyi larvae. Our results revealed 351 and 138 differentially accessible region genes through comparisons of ATAC-seq data between stages 21 and 24 and between stages 24 and 25, respectively. A joint analysis of RNA-seq and ATAC-seq data revealed a correlation between chromatin accessibility and gene transcription. We further verified the tissue expression patterns of 12 different genes. Among them, Cs-matrix metalloproteinase 24 (MMP24) and Cs-krüppel-like factor 5 (KLF5) were highly expressed in notochord cells. Functional assay results demonstrated that both genes are necessary for notochord lumen formation and expansion. Finally, we performed motif enrichment analysis of the differentially accessible regions in different tailbud stages and summarized the potential roles of these motif-bearing transcription factors in larval development. Overall, our study found a correlation between gene expression and chromatin accessibility and provided a vital resource for understanding the mechanisms of the development of ascidian embryos.

Transcription factor interactions explain the context-dependent activity of CRX binding sites

The effects of transcription factor binding sites (TFBSs) on the activity of a cis-regulatory element (CRE) depend on the local sequence context. In rod photoreceptors, binding sites for the transcription factor (TF) Cone-rod homeobox (CRX) occur in both enhancers and silencers, but the sequence context that determines whether CRX binding sites contribute to activation or repression of transcription is not understood. To investigate the context-dependent activity of CRX sites, we fit neural network-based models to the activities of synthetic CREs composed of photoreceptor TFBSs. The models revealed that CRX binding sites consistently make positive, independent contributions to CRE activity, while negative homotypic interactions between sites cause CREs composed of multiple CRX sites to function as silencers. The effects of negative homotypic interactions can be overcome by the presence of other TFBSs that either interact cooperatively with CRX sites or make independent positive contributions to activity. The context-dependent activity of CRX sites is thus determined by the balance between positive heterotypic interactions, independent contributions of TFBSs, and negative homotypic interactions. Our findings explain observed patterns of activity among genomic CRX-bound enhancers and silencers, and suggest that enhancers may require diverse TFBSs to overcome negative homotypic interactions between TFBSs.

A digital twin reproducing gene regulatory network dynamics of early Ciona embryos indicates robust buffers in the network

How gene regulatory networks (GRNs) encode gene expression dynamics and how GRNs evolve are not well understood, although these problems have been studied extensively. We created a digital twin that accurately reproduces expression dynamics of 13 genes that initiate expression in 32-cell ascidian embryos. We first showed that gene expression patterns can be manipulated according to predictions by this digital model. Next, to simulate GRN rewiring, we changed regulatory functions that represented their regulatory mechanisms in the digital twin, and found that in 55 of 100 cases, removal of a single regulator from a conjunctive clause of Boolean functions did not theoretically alter qualitative expression patterns of these genes. In other words, we found that more than half the regulators gave theoretically redundant temporal or spatial information to target genes. We experimentally substantiated that the expression pattern of Nodal was maintained without one of these factors, Zfpm, by changing the upstream regulatory sequence of Nodal. Such robust buffers of regulatory mechanisms may provide a basis of enabling developmental system drift, or rewiring of GRNs without changing expression patterns of downstream genes, during evolution.

Active learning of enhancer and silencer regulatory grammar in photoreceptors

Cis-regulatory elements (CREs) direct gene expression in health and disease, and models that can accurately predict their activities from DNA sequences are crucial for biomedicine. Deep learning represents one emerging strategy to model the regulatory grammar that relates CRE sequence to function. However, these models require training data on a scale that exceeds the number of CREs in the genome. We address this problem using active machine learning to iteratively train models on multiple rounds of synthetic DNA sequences assayed in live mammalian retinas. During each round of training the model actively selects sequence perturbations to assay, thereby efficiently generating informative training data. We iteratively trained a model that predicts the activities of sequences containing binding motifs for the photoreceptor transcription factor Cone-rod homeobox (CRX) using an order of magnitude less training data than current approaches. The model's internal confidence estimates of its predictions are reliable guides for designing sequences with high activity. The model correctly identified critical sequence differences between active and inactive sequences with nearly identical transcription factor binding sites, and revealed order and spacing preferences for combinations of motifs. Our results establish active learning as an effective method to train accurate deep learning models of cis-regulatory function after exhausting naturally occurring training examples in the genome.