The hagfish genome and the evolution of vertebrates.

As the only surviving lineages of jawless fishes, hagfishes and lampreys provide a crucial window into early vertebrate evolution1-3. Here we investigate the complex history, timing and functional role of genome-wide duplications4-7 and programmed DNA elimination8,9 in vertebrates in the light of a chromosome-scale genome sequence for the brown hagfish Eptatretus atami. Combining evidence from syntenic and phylogenetic analyses, we establish a comprehensive picture of vertebrate genome evolution, including an auto-tetraploidization (1RV) that predates the early Cambrian cyclostome-gnathostome split, followed by a mid-late Cambrian allo-tetraploidization (2RJV) in gnathostomes and a prolonged Cambrian-Ordovician hexaploidization (2RCY) in cyclostomes. Subsequently, hagfishes underwent extensive genomic changes, with chromosomal fusions accompanied by the loss of genes that are essential for organ systems (for example, genes involved in the development of eyes and in the proliferation of osteoclasts); these changes account, in part, for the simplification of the hagfish body plan1,2. Finally, we characterize programmed DNA elimination in hagfish, identifying protein-coding genes and repetitive elements that are deleted from somatic cell lineages during early development. The elimination of these germline-specific genes provides a mechanism for resolving genetic conflict between soma and germline by repressing germline and pluripotency functions, paralleling findings in lampreys10,11. Reconstruction of the early genomic history of vertebrates provides a framework for further investigations of the evolution of cyclostomes and jawed vertebrates.

De novo genome assembly and in natura epigenomics reveal salinity-induced DNA methylation in the mangrove tree Bruguiera gymnorhiza.

Mangroves are adapted to harsh environments, such as high ultraviolet (UV) light, low nutrition, and fluctuating salinity in coastal zones. However, little is known about the transcriptomic and epigenomic basis of the resilience of mangroves due to limited available genome resources. We performed a de novo genome assembly and in natura epigenome analyses of the mangrove Bruguiera gymnorhiza, one of the dominant mangrove species. We also performed the first genome-guided transcriptome assembly for mangrove species. The 309 Mb of the genome is predicted to encode 34 403 genes and has a repeat content of 48%. Depending on its growing environment, the natural B. gymnorhiza population showed drastic morphological changes associated with expression changes in thousands of genes. Moreover, high-salinity environments induced genome-wide DNA hypermethylation of transposable elements (TEs) in the B. gymnorhiza. DNA hypermethylation was concurrent with the transcriptional regulation of chromatin modifier genes, suggesting robust epigenome regulation of TEs in the B. gymnorhiza genome under high-salinity environments. The genome and epigenome data in this study provide novel insights into the epigenome regulation of mangroves and a better understanding of the adaptation of plants to fluctuating, harsh natural environments.

Transcriptional regulation analysis reveals the complexity of metamorphosis in the Pacific oyster (Crassostrea gigas).

Many marine invertebrate phyla are characterized by indirect development. These animals transit from planktonic larvae to benthic spats via settlement and metamorphosis, which contributes to their adaption to the marine environment. Studying the biological process of metamorphosis is, thus, key to understanding the origin and evolution of indirect development. Although numerous studies have been conducted on the relationship between metamorphosis and the marine environment, microorganisms, and neurohormones, little is known about gene regulation network (GRN) dynamics during metamorphosis. Metamorphosis-competent pediveligers of the Pacific oyster Crassostrea gigas were assayed in this study. By assaying gene expression patterns and open chromatin region changes of different samples of larvae and spats, the dynamics of molecular regulation during metamorphosis were examined. The results indicated significantly different gene regulation networks before, during and post-metamorphosis. Genes encoding membrane-integrated receptors and those related to the remodeling of the nervous system were upregulated before the initiation of metamorphosis. Massive biogenesis, e.g., of various enzymes and structural proteins, occurred during metamorphosis as inferred from the comprehensive upregulation of the protein synthesis system post epinephrine stimulation. Hierarchical downstream gene networks were then stimulated. Some transcription factors, including homeobox, basic helix-loop-helix and nuclear receptors, showed different temporal response patterns, suggesting a complex GRN during the transition stage. Nuclear receptors, as well as their retinoid X receptor partner, may participate in the GRN controlling oyster metamorphosis, indicating an ancient role of the nuclear receptor regulation system in animal metamorphosis. Supplementary Information: The online version contains supplementary material available at 10.1007/s42995-023-00204-y.

The hagfish genome and the evolution of vertebrates.

As the only surviving lineages of jawless fishes, hagfishes and lampreys provide a critical window into early vertebrate evolution. Here, we investigate the complex history, timing, and functional role of genome-wide duplications in vertebrates in the light of a chromosome-scale genome of the brown hagfish Eptatretus atami. Using robust chromosome-scale (paralogon-based) phylogenetic methods, we confirm the monophyly of cyclostomes, document an auto-tetraploidization (1RV) that predated the origin of crown group vertebrates ~517 Mya, and establish the timing of subsequent independent duplications in the gnathostome and cyclostome lineages. Some 1RV gene duplications can be linked to key vertebrate innovations, suggesting that this early genomewide event contributed to the emergence of pan-vertebrate features such as neural crest. The hagfish karyotype is derived by numerous fusions relative to the ancestral cyclostome arrangement preserved by lampreys. These genomic changes were accompanied by the loss of genes essential for organ systems (eyes, osteoclast) that are absent in hagfish, accounting in part for the simplification of the hagfish body plan; other gene family expansions account for hagfishes' capacity to produce slime. Finally, we characterise programmed DNA elimination in somatic cells of hagfish, identifying protein-coding and repetitive elements that are deleted during development. As in lampreys, the elimination of these genes provides a mechanism for resolving genetic conflict between soma and germline by repressing germline/pluripotency functions. Reconstruction of the early genomic history of vertebrates provides a framework for further exploration of vertebrate novelties.

Creation of an albino squid line by CRISPR-Cas9 and its application for in vivo functional imaging of neural activity.

Cephalopods are remarkable among invertebrates for their cognitive abilities, adaptive camouflage, novel structures, and propensity for recoding proteins through RNA editing. Due to the lack of genetically tractable cephalopod models, however, the mechanisms underlying these innovations are poorly understood. Genome editing tools such as CRISPR-Cas9 allow targeted mutations in diverse species to better link genes and function. One emerging cephalopod model, Euprymna berryi, produces large numbers of embryos that can be easily cultured throughout their life cycle and has a sequenced genome. As proof of principle, we used CRISPR-Cas9 in E. berryi to target the gene for tryptophan 2,3 dioxygenase (TDO), an enzyme required for the formation of ommochromes, the pigments present in the eyes and chromatophores of cephalopods. CRISPR-Cas9 ribonucleoproteins targeting tdo were injected into early embryos and then cultured to adulthood. Unexpectedly, the injected specimens were pigmented, despite verification of indels at the targeted sites by sequencing in injected animals (G0s). A homozygote knockout line for TDO, bred through multiple generations, was also pigmented. Surprisingly, a gene encoding indoleamine 2,3, dioxygenase (IDO), an enzyme that catalyzes the same reaction as TDO in vertebrates, was also present in E. berryi. Double knockouts of both tdo and ido with CRISPR-Cas9 produced an albino phenotype. We demonstrate the utility of these albinos for in vivo imaging of Ca2+ signaling in the brain using two-photon microscopy. These data show the feasibility of making gene knockout cephalopod lines that can be used for live imaging of neural activity in these behaviorally sophisticated organisms.

Analysis of the P. lividus sea urchin genome highlights contrasting trends of genomic and regulatory evolution in deuterostomes.

Sea urchins are emblematic models in developmental biology and display several characteristics that set them apart from other deuterostomes. To uncover the genomic cues that may underlie these specificities, we generated a chromosome-scale genome assembly for the sea urchin Paracentrotus lividus and an extensive gene expression and epigenetic profiles of its embryonic development. We found that, unlike vertebrates, sea urchins retained ancestral chromosomal linkages but underwent very fast intrachromosomal gene order mixing. We identified a burst of gene duplication in the echinoid lineage and showed that some of these expanded genes have been recruited in novel structures (water vascular system, Aristotle's lantern, and skeletogenic micromere lineage). Finally, we identified gene-regulatory modules conserved between sea urchins and chordates. Our results suggest that gene-regulatory networks controlling development can be conserved despite extensive gene order rearrangement.

Evidence from oyster suggests an ancient role for Pdx in regulating insulin gene expression in animals.

Hox and ParaHox genes encode transcription factors with similar expression patterns in divergent animals. The Pdx (Xlox) homeobox gene, for example, is expressed in a sharp spatial domain in the endodermal cell layer of the gut in chordates, echinoderms, annelids and molluscs. The significance of comparable gene expression patterns is unclear because it is not known if downstream transcriptional targets are also conserved. Here, we report evidence indicating that a classic transcriptional target of Pdx1 in vertebrates, the insulin gene, is a likely direct target of Pdx in Pacific oyster adults. We show that one insulin-related gene, cgILP, is co-expressed with cgPdx in oyster digestive tissue. Transcriptomic comparison suggests that this tissue plays a similar role to the vertebrate pancreas. Using ATAC-seq and ChIP, we identify an upstream regulatory element of the cgILP gene which shows binding interaction with cgPdx protein in oyster hepatopancreas and demonstrate, using a cell culture assay, that the oyster Pdx can act as a transcriptional activator through this site, possibly in synergy with NeuroD. These data argue that a classic homeodomain-target gene interaction dates back to the origin of Bilateria.

Macrophages directly contribute collagen to scar formation during zebrafish heart regeneration and mouse heart repair.

Canonical roles for macrophages in mediating the fibrotic response after a heart attack include extracellular matrix turnover and activation of cardiac fibroblasts to initiate collagen deposition. Here we reveal that macrophages directly contribute collagen to the forming post-injury scar. Unbiased transcriptomics shows an upregulation of collagens in both zebrafish and mouse macrophages following heart injury. Adoptive transfer of macrophages, from either collagen-tagged zebrafish or adult mouse GFPtpz-collagen donors, enhances scar formation via cell autonomous production of collagen. In zebrafish, the majority of tagged collagen localises proximal to the injury, within the overlying epicardial region, suggesting a possible distinction between macrophage-deposited collagen and that predominantly laid-down by myofibroblasts. Macrophage-specific targeting of col4a3bpa and cognate col4a1 in zebrafish significantly reduces scarring in cryoinjured hosts. Our findings contrast with the current model of scarring, whereby collagen deposition is exclusively attributed to myofibroblasts, and implicate macrophages as direct contributors to fibrosis during heart repair.

Reconstruction of the Global Neural Crest Gene Regulatory Network In Vivo.

Precise control of developmental processes is encoded in the genome in the form of gene regulatory networks (GRNs). Such multi-factorial systems are difficult to decode in vertebrates owing to their complex gene hierarchies and dynamic molecular interactions. Here we present a genome-wide in vivo reconstruction of the GRN underlying development of the multipotent neural crest (NC) embryonic cell population. By coupling NC-specific epigenomic and transcriptional profiling at population and single-cell levels with genome/epigenome engineering in vivo, we identify multiple regulatory layers governing NC ontogeny, including NC-specific enhancers and super-enhancers, novel trans-factors, and cis-signatures allowing reverse engineering of the NC-GRN at unprecedented resolution. Furthermore, identification and dissection of divergent upstream combinatorial regulatory codes has afforded new insights into opposing gene circuits that define canonical and neural NC fates early during NC ontogeny. Our integrated approach, allowing dissection of cell-type-specific regulatory circuits in vivo, has broad implications for GRN discovery and investigation.

A genome-wide assessment of the ancestral neural crest gene regulatory network.

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.

Generation of a double binary transgenic zebrafish model to study myeloid gene regulation in response to oncogene activation in melanocytes.

A complex network of inflammatory genes is closely linked to 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 regulatory programming of the myeloid cells as they respond to oncogene-activated melanocytes to be explored, focussing on the initial phase when cells become the precursors of cancer. A hormone-inducible binary system allows for temporal control of expression of different Ras oncogenes (NRasQ61K, HRasG12V and KRasG12V) 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.

From Pioneer to Repressor: Bimodal foxd3 Activity Dynamically Remodels Neural Crest Regulatory Landscape In Vivo.

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.

Active nuclear transcriptome analysis reveals inflammasome-dependent mechanism for early neutrophil response to Mycobacterium marinum.

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.

Biotagging of Specific Cell Populations in Zebrafish Reveals Gene Regulatory Logic Encoded in the Nuclear Transcriptome.

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.

BRG1-SWI/SNF-dependent regulation of the Wt1 transcriptional landscape mediates epicardial activity during heart development and disease.

Epicardium-derived cells (EPDCs) contribute cardiovascular cell types during development and in adulthood respond to Thymosin ß4 (Tß4) and myocardial infarction (MI) by reactivating a fetal gene programme to promote neovascularization and cardiomyogenesis. The mechanism for epicardial gene (re-)activation remains elusive. Here we reveal that BRG1, the essential ATPase subunit of the SWI/SNF chromatin-remodelling complex, is required for expression of Wilms' tumour 1 (Wt1), fetal EPDC activation and subsequent differentiation into coronary smooth muscle, and restores Wt1 activity upon MI. BRG1 physically interacts with Tß4 and is recruited by CCAAT/enhancer-binding protein ß (C/EBPß) to discrete regulatory elements in the Wt1 locus. BRG1-Tß4 co-operative binding promotes optimal transcription of Wt1 as the master regulator of embryonic EPDCs. Moreover, chromatin immunoprecipitation-sequencing reveals BRG1 binding at further key loci suggesting SWI/SNF activity across the fetal epicardial gene programme. These findings reveal essential functions for chromatin-remodelling in the activation of EPDCs during cardiovascular development and repair.

High-throughput spatial mapping of single-cell RNA-seq data to tissue of origin.

Understanding cell type identity in a multicellular organism requires the integration of gene expression profiles from individual cells with their spatial location in a particular tissue. Current technologies allow whole-transcriptome sequencing of spatially identified cells but lack the throughput needed to characterize complex tissues. Here we present a high-throughput method to identify the spatial origin of cells assayed by single-cell RNA-sequencing within a tissue of interest. Our approach is based on comparing complete, specificity-weighted mRNA profiles of a cell with positional gene expression profiles derived from a gene expression atlas. We show that this method allocates cells to precise locations in the brain of the marine annelid Platynereis dumerilii with a success rate of 81%. Our method is applicable to any system that has a reference gene expression database of sufficiently high resolution.

Thyrotroph embryonic factor regulates light-induced transcription of repair genes in zebrafish embryonic cells.

Numerous responses are triggered by light in the cell. How the light signal is detected and transduced into a cellular response is still an enigma. Each zebrafish cell has the capacity to directly detect light, making this organism particularly suitable for the study of light dependent transcription. To gain insight into the light signalling mechanism we identified genes that are activated by light exposure at an early embryonic stage, when specialised light sensing organs have not yet formed. We screened over 14,900 genes using micro-array GeneChips, and identified 19 light-induced genes that function primarily in light signalling, stress response, and DNA repair. Here we reveal that PAR Response Elements are present in all promoters of the light-induced genes, and demonstrate a pivotal role for the PAR bZip transcription factor Thyrotroph embryonic factor (Tef) in regulating the majority of light-induced genes. We show that tefbeta transcription is directly regulated by light while transcription of tefalpha is under circadian clock control at later stages of development. These data leads us to propose their involvement in light-induced UV tolerance in the zebrafish embryo.