The neural crest (NC) is a transient embryonic stem cell-like population characterized by its multipotency and broad developmental potential. Here, we perform NC-specific transcriptional and epigenomic profiling of foxd3-mutant cells in vivo to define the gene regulatory circuits controlling NC specification. Together with global binding analysis obtained by foxd3 biotin-ChIP and single cell profiles of foxd3-expressing premigratory NC, our analysis shows that, during early steps of NC formation, foxd3 acts globally as a pioneer factor to prime the onset of genes regulating NC specification and migration by re-arranging the chromatin landscape, opening cis-regulatory elements and reshuffling nucleosomes. Strikingly, foxd3 then gradually switches from an activator to its well-described role as a transcriptional repressor and potentially uses differential partners for each role. Taken together, these results demonstrate that foxd3 acts bimodally in the neural crest as a switch from ‘‘permissive’’ to ‘‘repressive’’ nucleosome and chromatin organization to maintain multipotency and define cell fates.

The neural crest (NC) is an embryonic cell population that contributes to key vertebrate-specific features including the craniofacial skeleton and peripheral nervous system. Here we examine the transcriptional and epigenomic profiles of NC cells in the sea lamprey, in order to gain insight into the ancestral state of the NC gene regulatory network (GRN). Transcriptome analyses identify clusters of co-regulated genes during NC specification and migration that show high conservation across vertebrates but also identify transcription factors (TFs) and cell-adhesion molecules not previously implicated in NC migration. ATAC-seq analysis uncovers an ensemble of cis-regulatory elements, including enhancers of Tfap2B, SoxE1 and Hox-α2 validated in the embryo. Cross-species deployment of lamprey elements identifies the deep conservation of lamprey SoxE1 enhancer activity, mediating homologous expression in jawed vertebrates. Our data provide insight into the core GRN elements conserved to the base of the vertebrates and expose others that are unique to lampreys.

Interrogation of gene regulatory circuits in complex organisms requires precise and robust methods to label cell-types for profiling of target proteins in a tissue-specific fashion as well as data analysis to understand interconnections within the circuits. There are several strategies for obtaining cell-type and subcellular specific genome-wide data. We have developed a methodology, termed “biotagging” that uses tissue-specific, genetically encoded components to biotinylate target proteins, enabling in depth genome-wide profiling in zebrafish. We have refined protocols to use the biotagging approach that led to enhanced isolation of coding and non-coding RNAs from ribosomes and nuclei of genetically defined cell-types. The ability to study both the actively translated and transcribed transcriptome in the same cell population, coupled to genomic accessibility assays has enabled the study of cell-type specific gene regulatory circuits in zebrafish due to the high signal-to-noise achieved via its stringent purification protocol. Here, we provide detailed methods to isolate, profile and analyze cell-type specific polyribosome and nuclear transcriptome in zebrafish.

Due to its genetic amenability coupled with recent advances in genome editing, the zebrafish serves as an excellent model to examine the function of both coding and non-coding elements. Recently, the non-coding genome has gained prominence due to its critical role in development and disease. Here, we have re-purposed the Ac/Ds maize transposition system to reliably screen and efficiently characterise zebrafish enhancers, with or without germline propagation. We further utilised the system to stably express guide RNAs in microinjected embryos enabling tissue-specific CRISPR/dCas9-interference (CRISPRi) knockdown of lncRNA and enhancer activity without disrupting the underlying genetic sequence. Our study highlights the utility of Ac/Ds transposition for transient epigenome modulation of non-coding elements in zebrafish.

A complex network of inflammation succeeds somatic cell transformation and malignant disease. Immune cells and their associated molecules are responsible for detecting and eliminating cancer cells as they establish themselves as the precursors of a tumour. By the time a patient has a detectable solid tumour, cancer cells have escaped the initial immune response mechanisms. Here, we describe the development of a double binary zebrafish model that enables exploring regulatory programming of the myeloid cells as they respond to oncogene-activated melanocytes, focussing on the initial phase when cells become the precursors of cancer. A hormone-inducible binary system allows for temporal control of different Ras-oncogenes (NRasQ61K, HRasG12V, KRasG12V) expression in melanocytes leading to proliferation and changes in morphology of the melanocytes. This model was coupled to binary cell-specific biotagging models allowing in vivo biotinylation and subsequent isolation of macrophage or neutrophil nuclei for regulatory profiling of their active transcriptomes. Nuclear transcriptional profiling of neutrophils, performed as they respond to the earliest precursors of melanoma in vivo, revealed an intricate landscape of regulatory factors that may promote progression to melanoma including serpinb1l4, fgf1, fgf6, cathepsin H, galectin 1 and galectin 3. The model presented here provides a powerful platform to study the myeloid response to the earliest precursors of melanoma.

CRISPR-Cas9 genome engineering has revolutionised all aspects of biological research, with epigenome engineering transforming gene regulation studies. Here, we present an optimised, adaptable toolkit enabling genome and epigenome engineering in the chicken embryo, and demonstrate its utility by probing gene regulatory interactions mediated by neural crest enhancers. First, we optimise novel efficient guide-RNA mini expression vectors utilising chick U6 promoters, provide a strategy for rapid somatic gene knockout and establish a protocol for evaluation of mutational penetrance by targeted next generation sequencing. We show that CRISPR/Cas9-mediated disruption of transcription factors causes a reduction in their cognate enhancer-driven reporter activity. Next, we assess endogenous enhancer function using both enhancer deletion and nuclease-deficient Cas9 (dCas9) effector fusions to modulate enhancer chromatin landscape, thus providing the first report of epigenome engineering in a developing embryo. Finally, we use the synergistic activation mediator (SAM) system to activate an endogenous target promoter. The novel genome and epigenome engineering toolkit developed here enables manipulation of endogenous gene expression and enhancer activity in chicken embryos, facilitating high-resolution analysis of gene regulatory interactions in vivo.

Circular RNAs are a puzzling class of RNAs with covalently closed loop structure, lacking 5′ and 3′ ends. Present in all cells, no unifying theme has emerged regarding their function. Here, we use transcriptional data obtained by biotagging in zebrafish to uncover a high-resolution cohort of embryonic circRNAs expressed in nuclear and polysomal subcellular compartments in three developmental cell types. The ample presence of embryonic circRNAs on polysomes contradicts previous reports suggesting predominant nuclear localisation. We uncover a novel class of circRNAs, significantly enriched at tandem duplicated genes. Using newly-developed NGS-based approach, we simultaneously resolve the full sequence of putative "tandem-circRNAs" and their underlying mRNAs. These form by long-range cis-splicing events often between distant tandem duplicated genes, resulting in chimaeric mRNA transcripts and circRNAs containing their supernumerary excised exons, integrated from multiple tandem loci. Taken together, our results suggest that circularisation events rather than circRNAs themselves are functionally important.

Interrogation of gene regulatory circuits in complex organisms requires precise tools for the selection of individual cell types and robust methods for biochemical profiling of target proteins. We have developed a versatile, tissue-specific binary in vivo biotinylation system in zebrafish termed biotagging that uses genetically encoded components to biotinylate target proteins, enabling in-depth genome- wide analyses of their molecular interactions. Using tissue-specific drivers and cell compartment-specific effector lines, we demonstrate the specificity of the biotagging toolkit at the biochemical, cellular, and transcriptional levels. We use biotagging to characterize the in vivo transcriptional landscape of migratory neural crest and myocardial cells in different cellular compartments (ribosomes and nucleus). These analyses reveal a comprehensive network of coding and non-coding RNAs and cis-regulatory modules, demonstrating that tissue-specific identity is embedded in the nuclear transcriptomes. By eliminating background inherent to complex embryonic environments, biotagging allows analyses of molecular interactions at high resolution.

The mechanisms governing neutrophil response to Mycobacterium tuberculosis remain poorly understood. In this study we utilise biotagging, a novel genome-wide profiling approach based on cell type-specific in vivo biotinylation in zebrafish to analyse the initial response of neutrophils to Mycobacterium marinum, a close genetic relative of M. tuberculosis used to model tuberculosis. Differential expression analysis following nuclear RNA-seq of neutrophil active transcriptomes reveals a significant upregulation in both damage-sensing and effector components of the inflammasome, including caspase b, NLRC3 ortholog (wu: fb15h11) and il1β. Crispr/Cas9-mediated knockout of caspase b, which acts by proteolytic processing of il1β, results in increased bacterial burden and less infiltration of macrophages to sites of mycobacterial infection, thus impairing granuloma development. We also show that a number of immediate early response genes (IEGs) are responsible for orchestrating the initial neutrophil response to mycobacterial infection. Further perturbation of the IEGs exposes egr3 as a key transcriptional regulator controlling il1β transcription.

Current screening methods for ovarian cancer can only detect advanced disease. Earlier detection has proved difficult because the molecular precursors involved in the natural history of the disease are unknown. To identify early driver mutations in ovarian cancer cells, we used dense whole genome sequencing of micrometastases and microscopic residual disease collected at three time points over three years from a single patient during treatment for high-grade serous ovarian cancer (HGSOC). The functional and clinical significance of the identified mutations was examined using a combination of population-based whole genome sequencing, targeted deep sequencing, multi-center analysis of protein expression, loss of function experiments in an in-vivo reporter assay and mammalian models, and gain of function experiments in primary cultured fallopian tube epithelial (FTE) cells. We identified frequent mutations involving a 40 kb distal repressor region for the key stem cell differentiation gene SOX2. In the apparently normal FTE, the region was also mutated. This was associated with a profound increase in SOX2 expression (p < 2−16), which was not found in patients without cancer (n = 108). Importantly, we show that SOX2 overexpression in FTE is nearly ubiquitous in patients with HGSOCs (n = 100), and common in BRCA1-BRCA2 mutation carriers (n = 71) who underwent prophylactic salpingo-oophorectomy. We propose that the finding of SOX2 overexpression in FTE could be exploited to develop biomarkers for detecting disease at a premalignant stage, which would reduce mortality from this devastating disease.