Highly Efficient CRISPR-Cas9-Based Methods for Generating Deletion Mutations and F0 Embryos that Lack Gene Function in Zebrafish

Endothelial Cell Flow-Mediated Quiescence Is Temporally Regulated and Utilizes the Cell Cycle Inhibitor p27

BACKGROUND

Endothelial cells regulate their cell cycle as blood vessels remodel and transition to quiescence downstream of blood flow-induced mechanotransduction. Laminar blood flow leads to quiescence, but how flow-mediated quiescence is established and maintained is poorly understood.

METHODS

Primary human endothelial cells were exposed to laminar flow regimens and gene expression manipulations, and quiescence depth was analyzed via time-to-cell cycle reentry after flow cessation. Mouse and zebrafish endothelial expression patterns were examined via scRNA-seq (single-cell RNA sequencing) analysis, and mutant or morphant fish lacking p27 were analyzed for endothelial cell cycle regulation and in vivo cellular behaviors.

RESULTS

Arterial flow-exposed endothelial cells had a distinct transcriptome, and they first entered a deep quiescence, then transitioned to shallow quiescence under homeostatic maintenance conditions. In contrast, venous flow-exposed endothelial cells entered deep quiescence early that did not change with homeostasis. The cell cycle inhibitor p27 (CDKN1B) was required to establish endothelial flow-mediated quiescence, and expression levels positively correlated with quiescence depth. p27 loss in vivo led to endothelial cell cycle upregulation and ectopic sprouting, consistent with loss of quiescence. HES1 and ID3, transcriptional repressors of p27 upregulated by arterial flow, were required for quiescence depth changes and the reduced p27 levels associated with shallow quiescence.

CONCLUSIONS

Endothelial cell flow-mediated quiescence has unique properties and temporal regulation of quiescence depth that depends on the flow stimulus. These findings are consistent with a model whereby flow-mediated endothelial cell quiescence depth is temporally regulated downstream of p27 transcriptional regulation by HES1 and ID3. The findings are important in understanding endothelial cell quiescence misregulation that leads to vascular dysfunction and disease.

CRISPR-RfxCas13d screening uncovers Bckdk as a post-translational regulator of the maternal-to-zygotic transition in teleosts

The Maternal-to-Zygotic transition (MZT) is a reprograming process encompassing zygotic genome activation (ZGA) and the clearance of maternally-provided mRNAs. While some factors regulating MZT have been identified, there are thousands of maternal RNAs whose function has not been ascribed yet. Here, we have performed a proof-of-principle CRISPR-RfxCas13d maternal screening targeting mRNAs encoding protein kinases and phosphatases in zebrafish and identified Bckdk as a novel post-translational regulator of MZT. Bckdk mRNA knockdown caused epiboly defects, ZGA deregulation, H3K27ac reduction and a partial impairment of miR-430 processing. Phospho-proteomic analysis revealed that Phf10/Baf45a, a chromatin remodeling factor, is less phosphorylated upon Bckdk depletion. Further, phf10 mRNA knockdown also altered ZGA and Phf10 constitutively phosphorylated rescued the developmental defects observed after bckdk mRNA depletion. Altogether, our results demonstrate the competence of CRISPR-RfxCas13d screenings to uncover new regulators of early vertebrate development and shed light on the post-translational control of MZT mediated by protein phosphorylation.

Functional role of myosin-binding protein H in thick filaments of developing vertebrate fast-twitch skeletal muscle

Myosin-binding protein H (MyBP-H) is a component of the vertebrate skeletal muscle sarcomere with sequence and domain homology to myosin-binding protein C (MyBP-C). Whereas skeletal muscle isoforms of MyBP-C (fMyBP-C, sMyBP-C) modulate muscle contractility via interactions with actin thin filaments and myosin motors within the muscle sarcomere "C-zone," MyBP-H has no known function. This is in part due to MyBP-H having limited expression in adult fast-twitch muscle and no known involvement in muscle disease. Quantitative proteomics reported here reveal MyBP-H is highly expressed in prenatal rat fast-twitch muscles and larval zebrafish, suggesting a conserved role in muscle development, and promoting studies to define its function. We take advantage of the genetic control of the zebrafish model and a combination of structural, functional, and biophysical techniques to interrogate the role of MyBP-H. Transgenic, FLAG-tagged MyBP-H or fMyBP-C both localize to the C-zones in larval myofibers, whereas genetic depletion of endogenous MyBP-H or fMyBP-C leads to increased accumulation of the other, suggesting competition for C-zone binding sites. Does MyBP-H modulate contractility from the C-zone? Globular domains critical to MyBP-C's modulatory functions are absent from MyBP-H, suggesting MyBP-H may be functionally silent. However, our results suggest an active role. Small angle x-ray diffraction of intact larval tails revealed MyBP-H contributes to the compression of the myofilament lattice accompanying stretch or contraction, while in vitro motility experiments indicate MyBP-H shares MyBP-C's capacity as a molecular "brake". These results provide new insights and raise questions about the role of the C-zone during muscle development.

Genetic tools for the study of the mangrove killifish, Kryptolebias marmoratus, an emerging vertebrate model for phenotypic plasticity

Kryptolebias marmoratus (Kmar), a teleost fish of the order Cyprinodontiformes, has a suite of unique phenotypes and behaviors not observed in other fishes. Many of these phenotypes are discrete and highly plastic-varying over time within an individual, and in some cases reversible. Kmar and its interfertile sister species, K. hermaphroditus, are the only known self-fertile vertebrates. This unusual sexual mode has the potential to provide unique insights into the regulation of vertebrate sexual development, and also lends itself to genetics. Kmar is easily adapted to the lab and requires little maintenance. However, its internal fertilization and small clutch size limits its experimental use. To support Kmar as a genetic model, we compared alternative husbandry techniques to maximize recovery of early cleavage-stage embryos. We find that frequent egg collection enhances yield, and that protease treatment promotes the greatest hatching success. We completed a forward mutagenesis screen and recovered several mutant lines that serve as important tools for genetics in this model. Several will serve as useful viable recessive markers for marking crosses. Importantly, the mutant kissylips lays embryos at twice the rate of wild-type. Combining frequent egg collection with the kissylips mutant background allows for a substantial enhancement of early embryo yield. These improvements were sufficient to allow experimental analysis of early development and the successful mono- and bi-allelic targeted knockout of an endogenous tyrosinase gene with CRISPR/Cas9 nucleases. Collectively, these tools will facilitate modern developmental genetics in this fascinating fish, leading to future insights into the regulation of plasticity.

Differentiation and functioning of the lateral line organ in zebrafish require Smpx activity

The small muscle protein, X-linked (SMPX) gene encodes a cytoskeleton-associated protein, highly expressed in the inner ear hair cells (HCs), possibly regulating auditory function. In the last decade, several mutations in SMPX have been associated with X-chromosomal progressive non syndromic hearing loss in humans and, in line with this, Smpx-deficient animal models, namely zebrafish and mouse, showed significant impairment of inner ear HCs development, maintenance, and functioning. In this work, we uncovered smpx expression in the neuromast mechanosensory HCs of both Anterior and Posterior Lateral Line (ALL and PLL, respectively) of zebrafish larvae and focused our attention on the PLL. Smpx was subcellularly localized throughout the cytoplasm of the HCs, as well as in their primary cilium. Loss-of-function experiments, via both morpholino-mediated gene knockdown and CRISPR/Cas9 F0 gene knockout, revealed that the lack of Smpx led to fewer properly differentiated and functional neuromasts, as well as to a smaller PLL primordium (PLLp), the latter also Smpx-positive. In addition, the kinocilia of Smpx-deficient neuromast HCs appeared structurally and numerically altered. Such phenotypes were associated with a significant reduction in the mechanotransduction activity of the neuromast HCs, in line with their positivity for Smpx. In summary, this work highlights the importance of Smpx in lateral line development and, specifically, in proper HCs differentiation and/or maintenance, and in the mechanotransduction process carried out by the neuromast HCs. Because lateral line HCs are both functionally and structurally analogous to the cochlear HCs, the neuromasts might represent an invaluable-and easily accessible-tool to dissect the role of Smpx in HCs development/functioning and shed light on the underlying mechanisms involved in hearing loss.

Foxp and Skor family proteins control differentiation of Purkinje cells from Ptf1a- and Neurog1-expressing progenitors in zebrafish

Cerebellar neurons, such as GABAergic Purkinje cells (PCs), interneurons (INs) and glutamatergic granule cells (GCs) are differentiated from neural progenitors expressing proneural genes, including ptf1a, neurog1 and atoh1a/b/c. Studies in mammals previously suggested that these genes determine cerebellar neuron cell fate. However, our studies on ptf1a;neurog1 zebrafish mutants and lineage tracing of ptf1a-expressing progenitors have revealed that the ptf1a/neurog1-expressing progenitors can generate diverse cerebellar neurons, including PCs, INs and a subset of GCs in zebrafish. The precise mechanisms of how each cerebellar neuron type is specified remains elusive. We found that genes encoding the transcriptional regulators Foxp1b, Foxp4, Skor1b and Skor2, which are reportedly expressed in PCs, were absent in ptf1a;neurog1 mutants. foxp1b;foxp4 mutants showed a strong reduction in PCs, whereas skor1b;skor2 mutants completely lacked PCs, and displayed an increase in immature GCs. Misexpression of skor2 in GC progenitors expressing atoh1c suppressed GC fate. These data indicate that Foxp1b/4 and Skor1b/2 function as key transcriptional regulators in the initial step of PC differentiation from ptf1a/neurog1-expressing neural progenitors, and that Skor1b and Skor2 control PC differentiation by suppressing their differentiation into GCs.

Distinct roles of core autophagy-related genes in zebrafish definitive hematopoiesis

Despite the well-described discrepancy between ATG (macroautophagy/autophagy-related) genes in the regulation of hematopoiesis, varying essentiality of core ATG proteins in vertebrate definitive hematopoiesis remains largely unclear. Here, we employed zebrafish (Danio rerio) to compare the functions of six core atg genes, including atg13, becn1 (beclin1), atg9a, atg2a, atg5, and atg3, in vertebrate definitive hematopoiesis via clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9 ribonucleoprotein and morpholino targeting. Zebrafish with various atg mutations showed autophagic deficiency and presented partially consistent hematopoietic abnormalities during early development. All six atg mutations led to a declined number of spi1b+ (Spi-1 proto-oncogene b) myeloid progenitor cells. However, only becn1 mutation resulted in the expansion of myb+ (v-myb avian myeloblastosis viral oncogene homolog) hematopoietic stem and progenitor cells (HSPCs) and transiently increased coro1a+ (coronin, actin binding protein, 1A) leukocytes, whereas atg3 mutation decreased the number of HSPCs and leukocytes. Proteomic analysis of caudal hematopoietic tissue identified sin3aa (SIN3 transcription regulator family member Aa) as a potential modulator of atg13- and becn1-regulated definitive hematopoiesis. Disruption of sin3aa rescued the expansion of HSPCs and leukocytes in becn1 mutants and exacerbated the decrease of HSPCs in atg13 mutants. Double mutations were also performed to examine alternative functions of various atg genes in definitive hematopoiesis. Notably, becn1 mutation failed to induce HSPCs expansion with one of the other five atg mutations. These findings demonstrated the distinct roles of atg genes and their interplays in zebrafish definitive hematopoiesis, thereby suggesting that the vertebrate definitive hematopoiesis is regulated in an atg gene-dependent manner.Abbreviations: AGM: aorta-gonad-mesonephros; AO: acridine orange; atg: autophagy related; becn1: beclin 1, autophagy related; CHT: caudal hematopoietic tissue; CKO: conditional knockout; coro1a: coronin, actin binding protein, 1A; CQ: chloroquine; CRISPR: clustered regularly interspaced short palindromic repeats; dpf: days post fertilization; FACS: fluorescence-activated cell sorting; hbae1.1: hemoglobin, alpha embryonic 1.1; HSCs: hematopoietic stem cells; HSPCs: hematopoietic stem and progenitor cells; KD: knockdown; KO: knockout; map1lc3/lc3: microtubule-associated protein 1 light chain 3; MO: morpholino; mpeg1.1: macrophage expressed 1, tandem duplicate 1; mpx: myeloid-specific peroxidase; myb: v-myb avian myeloblastosis viral oncogene homolog; PE: phosphatidylethanolamine; p-H3: phospho-H3 histone; PtdIns3K: class 3 phosphatidylinositol 3-kinase; rag1: recombination activating 1; rb1cc1/fip200: RB1-inducible coiled-coil 1; RFLP: restriction fragment length polymorphism; RNP: ribonucleoprotein; sin3aa: SIN3 transcription regulator family member Aa; spi1b: Spi-1 proto-oncogene b; ulk: unc-51 like autophagy activating kinase; vtg1: vitellogenin 1; WISH: whole-mount in situ hybridization.

EVOLUTIONARY CO-OPTION OF AN ANCESTRAL CLOACAL REGULATORY LANDSCAPE DURING THE EMERGENCE OF DIGITS AND GENITALS

The transition from fins to limbs has been a rich source of discussion for more than a century. One open and important issue is understanding how the mechanisms that pattern digits arose during vertebrate evolution. In this context, the analysis of Hox gene expression and functions to infer evolutionary scenarios has been a productive approach to explain the changes in organ formation, particularly in limbs. In tetrapods, the transcription of Hoxd genes in developing digits depends on a well-characterized set of enhancers forming a large regulatory landscape1,2. This control system has a syntenic counterpart in zebrafish, even though they lack bona fide digits, suggestive of deep homology3 between distal fin and limb developmental mechanisms. We tested the global function of this landscape to assess ancestry and source of limb and fin variation. In contrast to results in mice, we show here that the deletion of the homologous control region in zebrafish has a limited effect on the transcription of hoxd genes during fin development. However, it fully abrogates hoxd expression within the developing cloaca, an ancestral structure related to the mammalian urogenital sinus. We show that similar to the limb, Hoxd gene function in the urogenital sinus of the mouse also depends on enhancers located in this same genomic domain. Thus, we conclude that the current regulation underlying Hoxd gene expression in distal limbs was co-opted in tetrapods from a preexisting cloacal program. The orthologous chromatin domain in fishes may illustrate a rudimentary or partial step in this evolutionary co-option.

A Simple and Scalable Zebrafish Model of Sonic Hedgehog Medulloblastoma

Medulloblastoma (MB) is the most common malignant brain tumor in children and is stratified into three major subgroups. The Sonic hedgehog (SHH) subgroup represents ~30% of all MB cases and has significant survival disparity depending upon TP53 status. Here, we describe the first zebrafish model of SHH MB using CRISPR to mutate ptch1, the primary genetic driver in human SHH MB. These tumors rapidly arise adjacent to the valvula cerebelli and resemble human SHH MB by histology and comparative genomics. In addition, ptch1-deficient MB tumors with loss of tp53 have aggressive tumor histology and significantly worse survival outcomes, comparable to human patients. The simplicity and scalability of the ptch1 MB model makes it highly amenable to CRISPR-based genome editing screens to identify genes required for SHH MB tumor formation in vivo, and here we identify the grk3 kinase as one such target.

Scalable CRISPR Screens in Zebrafish Using MIC-Drop

CRISPR-Cas9 is a powerful tool to interrogate gene function in a targeted and systematic manner. Although the technology has been scaled up for large-scale genetic screens in cell culture, similar scale screens in vivo have been extremely challenging due to the cost, labor, and time required to generate and keep track of thousands of mutant animals. We reported the development of Multiplexed Intermixed CRISPR Droplets (MIC-Drop), a platform that makes large-scale reverse genetic screens possible in zebrafish. In this chapter, we provide a detailed protocol to conduct large-scale genetic screens using this novel platform.

Integrating non-mammalian model organisms in the diagnosis of rare genetic diseases in humans

Next-generation sequencing technology has rapidly accelerated the discovery of genetic variants of interest in individuals with rare diseases. However, showing that these variants are causative of the disease in question is complex and may require functional studies. Use of non-mammalian model organisms - mainly fruitflies (Drosophila melanogaster), nematode worms (Caenorhabditis elegans) and zebrafish (Danio rerio) - enables the rapid and cost-effective assessment of the effects of gene variants, which can then be validated in mammalian model organisms such as mice and in human cells. By probing mechanisms of gene action and identifying interacting genes and proteins in vivo, recent studies in these non-mammalian model organisms have facilitated the diagnosis of numerous genetic diseases and have enabled the screening and identification of therapeutic options for patients. Studies in non-mammalian model organisms have also shown that the biological processes underlying rare diseases can provide insight into more common mechanisms of disease and the biological functions of genes. Here, we discuss the opportunities afforded by non-mammalian model organisms, focusing on flies, worms and fish, and provide examples of their use in the diagnosis of rare genetic diseases.

Genetic mechanisms of multiciliated cell development: from fate choice to differentiation in zebrafish and other models

Multiciliated cells (MCCS) form bundles of cilia and their activities are essential for the proper development and physiology of many organ systems. Not surprisingly, defects in MCCs have profound consequences and are associated with numerous disease states. Here, we discuss the current understanding of MCC formation, with a special focus on the genetic and molecular mechanisms of MCC fate choice and differentiation. Furthermore, we cast a spotlight on the use of zebrafish to study MCC ontogeny and several recent advances made in understanding MCCs using this vertebrate model to delineate mechanisms of MCC emergence in the developing kidney.

Proteostasis governs differential temperature sensitivity across embryonic cell types

Embryonic development is remarkably robust, but temperature stress can degrade its ability to generate animals with invariant anatomy. Phenotypes associated with environmental stress suggest that some cell types are more sensitive to stress than others, but the basis of this sensitivity is unknown. Here, we characterize hundreds of individual zebrafish embryos under temperature stress using whole-animal single-cell RNA sequencing (RNA-seq) to identify cell types and molecular programs driving phenotypic variability. We find that temperature perturbs the normal proportions and gene expression programs of numerous cell types and also introduces asynchrony in developmental timing. The notochord is particularly sensitive to temperature, which we map to a specialized cell type: sheath cells. These cells accumulate misfolded protein at elevated temperature, leading to a cascading structural failure of the notochord and anatomic defects. Our study demonstrates that whole-animal single-cell RNA-seq can identify mechanisms for developmental robustness and pinpoint cell types that constitute key failure points.

Embryo-scale reverse genetics at single-cell resolution

The maturation of single-cell transcriptomic technologies has facilitated the generation of comprehensive cellular atlases from whole embryos1-4. A majority of these data, however, has been collected from wild-type embryos without an appreciation for the latent variation that is present in development. Here we present the 'zebrafish single-cell atlas of perturbed embryos': single-cell transcriptomic data from 1,812 individually resolved developing zebrafish embryos, encompassing 19 timepoints, 23 genetic perturbations and a total of 3.2 million cells. The high degree of replication in our study (eight or more embryos per condition) enables us to estimate the variance in cell type abundance organism-wide and to detect perturbation-dependent deviance in cell type composition relative to wild-type embryos. Our approach is sensitive to rare cell types, resolving developmental trajectories and genetic dependencies in the cranial ganglia neurons, a cell population that comprises less than 1% of the embryo. Additionally, time-series profiling of individual mutants identified a group of brachyury-independent cells with strikingly similar transcriptomes to notochord sheath cells, leading to new hypotheses about early origins of the skull. We anticipate that standardized collection of high-resolution, organism-scale single-cell data from large numbers of individual embryos will enable mapping of the genetic dependencies of zebrafish cell types, while also addressing longstanding challenges in developmental genetics, including the cellular and transcriptional plasticity underlying phenotypic diversity across individuals.

Investigating the role of lipid genes in liver disease using fatty liver models of alcohol and high fat in zebrafish (Danio rerio)

BACKGROUND

Accumulation of lipid in the liver is the first hallmark of both alcohol-related liver disease (ALD) and non-alcohol-related fatty liver disease (NAFLD). Recent studies indicate that specific mutations in lipid genes confer risk and might influence disease progression to irreversible liver cirrhosis. This study aimed to understand the function/s of lipid risk genes driving disease development in zebrafish genetic models of alcohol-related and non-alcohol-related fatty liver.

METHODS

We used zebrafish larvae to investigate the effect of alcohol and high fat to model fatty liver and tested the utility of this model to study lipid risk gene functions. CRISPR/Cas9 gene editing was used to create knockdowns in 5 days post-fertilisation zebrafish larvae for the available orthologs of human cirrhosis risk genes (pnpla3, faf2, tm6sf2). To establish fatty liver models, larvae were exposed to ethanol and a high-fat diet (HFD) consisting of chicken egg yolk. Changes in morphology (imaging), survival, liver injury (biochemical tests, histopathology), gene expression (qPCR) and lipid accumulation (dye-specific live imaging) were analysed across treatment groups to test the functions of these genes.

RESULTS

Exposure of 5-day post-fertilisation (dpf) WT larvae to 2% ethanol or HFD for 48 h developed measurable hepatic steatosis. CRISPR-Cas9 genome editing depleted pnpla3, faf2 and tm6sf2 gene expression in these CRISPR knockdown larvae (crispants). Depletion significantly increased the effects of ethanol and HFD toxicity by increasing hepatic steatosis and hepatic neutrophil recruitment ≥2-fold in all three crispants. Furthermore, ethanol or HFD exposure significantly altered the expression of genes associated with ethanol metabolism (cyp2y3) and lipid metabolism-related gene expression, including atgl (triglyceride hydrolysis), axox1, echs1 (fatty acid β-oxidation), fabp10a (transport), hmgcra (metabolism), notch1 (signalling) and srebp1 (lipid synthesis), in all three pnpla3, faf2 and tm6sf2 crispants. Nile Red staining in all three crispants revealed significantly increased lipid droplet size and triglyceride accumulation in the livers following exposure to ethanol or HFD.

CONCLUSIONS

We identified roles for pnpla3, faf2 and tm6sf2 genes in triglyceride accumulation and fatty acid oxidation pathways in a zebrafish larvae model of fatty liver.

Pioneer statoacoustic neurons guide neuroblast behaviour during otic ganglion assembly

Cranial ganglia are aggregates of sensory neurons that mediate distinct types of sensation. The statoacoustic ganglion (SAG) develops into several lobes that are spatially arranged to connect appropriately with hair cells of the inner ear. To investigate the cellular behaviours involved in the 3D organization of the SAG, we use high-resolution confocal imaging of single-cell, labelled zebrafish neuroblasts (NBs), photoconversion, photoablation, and genetic perturbations. We show that otic NBs delaminate out of the otic epithelium in an epithelial-mesenchymal transition-like manner, rearranging apical polarity and primary cilia proteins. We also show that, once delaminated, NBs require RhoGTPases in order to perform active migration. Furthermore, tracking of recently delaminated NBs revealed their directed migration and coalescence around a small population of pioneer SAG neurons. These pioneer SAG neurons, not from otic placode origin, populate the coalescence region before otic neurogenesis begins and their ablation disrupts delaminated NB migratory pathways, consequentially affecting SAG shape. Altogether, this work shows for the first time the role of pioneer SAG neurons in orchestrating SAG development.

Hand2 represses non-cardiac cell fates through chromatin remodeling at cis- regulatory elements

Developmental studies have revealed the importance of the transcription factor Hand2 in cardiac development. Hand2 promotes cardiac progenitor differentiation and epithelial maturation, while repressing other tissue types. The mechanisms underlying the promotion of cardiac fates are far better understood than those underlying the repression of alternative fates. Here, we assess Hand2-dependent changes in gene expression and chromatin remodeling in cardiac progenitors of zebrafish embryos. Cell-type specific transcriptome analysis shows a dual function for Hand2 in activation of cardiac differentiation genes and repression of pronephric pathways. We identify functional cis- regulatory elements whose chromatin accessibility are increased in hand2 mutant cells. These regulatory elements associate with non-cardiac gene expression, and drive reporter gene expression in tissues associated with Hand2-repressed genes. We find that functional Hand2 is sufficient to reduce non-cardiac reporter expression in cardiac lineages. Taken together, our data support a model of Hand2-dependent coordination of transcriptional programs, not only through transcriptional activation of cardiac and epithelial maturation genes, but also through repressive chromatin remodeling at the DNA regulatory elements of non-cardiac genes.

CRISPR/Cas9-induced asap1a and asap1b co-knockout mutant zebrafish displayed abnormal embryonic development and impaired neutrophil migration

ASAP1 (Arf-GAP with SH3 domain, the ankyrin repeat and the PH domain) is the GTPase activating protein of the small G protein Arf. To understand more about the physiological functions of ASAP1 in vivo, we chose to use the zebrafish as an animal model, and analyzed the characterization of asap1 using loss-of-function studies. Here, two isoforms in zebrafish, asap1a and asap1b, were found to be homologous to human ASAP1, and the gene knockout zebrafish lines for asap1a and asap1b were established using the CRISPR/Cas9 technique with different insertions and deletions of bases. Zebrafish with asap1a and asap1b co-knockout showed a significant reduction in survival and hatching rates, as well as an increase in malformation rates during the early stages of development, while the asap1a or asap1b single knockout mutants did not affect the growth and development of individual zebrafish. Exploring the gene expression compensation between asap1a and asap1b using qRT-PCR, we found that asap1b had increased expression when asap1a was knocked out, showing a clear compensatory effect against asap1a knockout; In turn, asap1a did not have detectable compensating expression after asap1b knockout. Furthermore, the co-knockout homozygous mutants displayed impaired neutrophil migration to Mycobacterium marinum infection, and showed an increased bacterial load. Together, these are the first inherited asap1a and/or asap1b mutant zebrafish lines by the CRISPR/Cas9 gene editing approach, and by serving as useful models, they can significantly contribute to better annotation and follow-up physiological studies of human ASAP1.

Fluorescent PCR-based Screening Methods for Precise Knock-in of Small DNA Fragments and Point Mutations in Zebrafish

Generation of zebrafish (Danio rerio) models with targeted insertion of epitope tags and point mutations is highly desirable for functional genomics and disease modeling studies. Currently, CRISPR/Cas9-mediated knock-in is the method of choice for insertion of exogeneous sequences by providing a repair template for homology-directed repair (HDR). A major hurdle in generating knock-in models is the labor and cost involved in screening of injected fish to identify the precise knock-in events due to low efficiency of the HDR pathway in zebrafish. Thus, we developed fluorescent PCR-based high-throughput screening methods for precise knock-in of epitope tags and point mutations in zebrafish. Here, we provide a step-by-step guide that describes selection of an active sgRNA near the intended knock-in site, design of single-stranded oligonucleotide (ssODN) templates for HDR, quick validation of somatic knock-in using injected embryos, and screening for germline transmission of precise knock-in events to establish stable lines. Our screening method relies on the size-based separation of all fragments in an amplicon by fluorescent PCR and capillary electrophoresis, thus providing a robust and cost-effective strategy. Although we present the use of this protocol for insertion of epitope tags and point mutations, it can be used for insertion of any small DNA fragments (e.g., LoxP sites, in-frame codons). Furthermore, the screening strategy described here can be used to screen for precise knock-in of small DNA sequences in any model system, as PCR amplification of the target region is its only requirement. Key features This protocol expands the use of fluorescent PCR and CRISPR-STAT for screening of precise knock-in of small insertions and point mutations in zebrafish. Allows validation of selected sgRNA and HDR template within two weeks by somatic knock-in screening. Allows robust screening of point mutations by combining restriction digest with CRISPR-STAT. Graphical overview Overview of the three-phase knock-in pipeline in zebrafish (created with BioRender.com).

CRISPR-based knockout and base editing confirm the role of MYRF in heart development and congenital heart disease

High-throughput DNA sequencing studies increasingly associate DNA variants with congenital heart disease (CHD). However, functional modeling is a crucial prerequisite for translating genomic data into clinical care. We used CRISPR-Cas9-mediated targeting of 12 candidate genes in the vertebrate model medaka (Oryzias latipes), five of which displayed a novel cardiovascular phenotype spectrum in F0 (crispants): mapre2, smg7, cdc42bpab, ankrd11 and myrf, encoding a transcription factor recently linked to cardiac-urogenital syndrome. Our myrf mutant line showed particularly prominent embryonic cardiac defects recapitulating phenotypes of pediatric patients, including hypoplastic ventricle. Mimicking human mutations, we edited three sites to generate specific myrf single-nucleotide variants via cytosine and adenine base editors. The Glu749Lys missense mutation in the conserved intramolecular chaperon autocleavage domain fully recapitulated the characteristic myrf mutant phenotype with high penetrance, underlining the crucial function of this protein domain. The efficiency and scalability of base editing to model specific point mutations accelerate gene validation studies and the generation of human-relevant disease models.

Single cell transcriptome analyses of the developing zebrafish eye- perspectives and applications

Within a relatively short period of time, single cell transcriptome analyses (SCT) have become increasingly ubiquitous with transcriptomic research, uncovering plentiful details that boost our molecular understanding of various biological processes. Stemming from SCT analyses, the ever-growing number of newly assigned genetic markers increases our understanding of general function and development, while providing opportunities for identifying genes associated with disease. SCT analyses have been carried out using tissue from numerous organisms. However, despite the great potential of zebrafish as a model organism, other models are still preferably used. In this mini review, we focus on eye research as an example of the advantages in using zebrafish, particularly its usefulness for single cell transcriptome analyses of developmental processes. As studies have already shown, the unique opportunities offered by zebrafish, including similarities to the human eye, in combination with the possibility to analyze and extract specific cells at distinct developmental time points makes the model a uniquely powerful one. Particularly the practicality of collecting large numbers of embryos and therefore isolation of sufficient numbers of developing cells is a distinct advantage compared to other model organisms. Lastly, the advent of highly efficient genetic knockouts methods offers opportunities to characterize target gene function in a more cost-efficient way. In conclusion, we argue that the use of zebrafish for SCT approaches has great potential to further deepen our molecular understanding of not only eye development, but also many other organ systems.

Large-scale F0 CRISPR screens in vivo using MIC-Drop

The zebrafish is a powerful model system for studying animal development, for modeling genetic diseases, and for large-scale in vivo functional genetics. Because of its ease of use and its high efficiency in targeted gene perturbation, CRISPR-Cas9 has recently gained prominence as the tool of choice for genetic manipulation in zebrafish. However, scaling up the technique for high-throughput in vivo functional genetics has been a challenge. We recently developed a method, Multiplexed Intermixed CRISPR Droplets (MIC-Drop), that makes large-scale CRISPR screening in zebrafish possible. Here, we outline the step-by-step protocol for performing functional genetic screens in zebrafish by using MIC-Drop. MIC-Drop uses multiplexed single-guide RNAs to generate biallelic mutations in injected zebrafish embryos, allowing genetic screens to be performed in F0 animals. Combining microfluidics and DNA barcoding enables simultaneous targeting of tens to hundreds of genes from a single injection needle, while also enabling retrospective and rapid identification of the genotype responsible for an observed phenotype. The primary target audiences for MIC-Drop are developmental biologists, zebrafish geneticists, and researchers interested in performing in vivo functional genetic screens in a vertebrate model system. MIC-Drop will also prove useful in the hands of chemical biologists seeking to identify targets of small molecules that cause phenotypic changes in zebrafish. By using MIC-Drop, a typical screen of 100 genes can be conducted within 2-3 weeks by a single user.

Nr2f1a maintains atrial nkx2.5 expression to repress pacemaker identity within venous atrial cardiomyocytes of zebrafish

Maintenance of cardiomyocyte identity is vital for normal heart development and function. However, our understanding of cardiomyocyte plasticity remains incomplete. Here, we show that sustained expression of the zebrafish transcription factor Nr2f1a prevents the progressive acquisition of ventricular cardiomyocyte (VC) and pacemaker cardiomyocyte (PC) identities within distinct regions of the atrium. Transcriptomic analysis of flow-sorted atrial cardiomyocytes (ACs) from nr2f1a mutant zebrafish embryos showed increased VC marker gene expression and altered expression of core PC regulatory genes, including decreased expression of nkx2.5, a critical repressor of PC differentiation. At the arterial (outflow) pole of the atrium in nr2f1a mutants, cardiomyocytes resolve to VC identity within the expanded atrioventricular canal. However, at the venous (inflow) pole of the atrium, there is a progressive wave of AC transdifferentiation into PCs across the atrium toward the arterial pole. Restoring Nkx2.5 is sufficient to repress PC marker identity in nr2f1a mutant atria and analysis of chromatin accessibility identified an Nr2f1a-dependent nkx2.5 enhancer expressed in the atrial myocardium directly adjacent to PCs. CRISPR/Cas9-mediated deletion of the putative nkx2.5 enhancer leads to a loss of Nkx2.5-expressing ACs and expansion of a PC reporter, supporting that Nr2f1a limits PC differentiation within venous ACs via maintaining nkx2.5 expression. The Nr2f-dependent maintenance of AC identity within discrete atrial compartments may provide insights into the molecular etiology of concurrent structural congenital heart defects and associated arrhythmias.

STAT5b is a key effector of NRG-1/ERBB4-mediated myocardial growth

The growth factor Neuregulin-1 (NRG-1) regulates myocardial growth and is currently under clinical investigation as a treatment for heart failure. Here, we demonstrate in several in vitro and in vivo models that STAT5b mediates NRG-1/EBBB4-stimulated cardiomyocyte growth. Genetic and chemical disruption of the NRG-1/ERBB4 pathway reduces STAT5b activation and transcription of STAT5b target genes Igf1, Myc, and Cdkn1a in murine cardiomyocytes. Loss of Stat5b also ablates NRG-1-induced cardiomyocyte hypertrophy. Dynamin-2 is shown to control the cell surface localization of ERBB4 and chemical inhibition of Dynamin-2 downregulates STAT5b activation and cardiomyocyte hypertrophy. In zebrafish embryos, Stat5 is activated during NRG-1-induced hyperplastic myocardial growth, and chemical inhibition of the Nrg-1/Erbb4 pathway or Dynamin-2 leads to loss of myocardial growth and Stat5 activation. Moreover, CRISPR/Cas9-mediated knockdown of stat5b results in reduced myocardial growth and cardiac function. Finally, the NRG-1/ERBB4/STAT5b signaling pathway is differentially regulated at mRNA and protein levels in the myocardium of patients with pathological cardiac hypertrophy as compared to control human subjects, consistent with a role of the NRG-1/ERBB4/STAT5b pathway in myocardial growth.

Application of an F0-based genetic assay in adult zebrafish to identify modifier genes of an inherited cardiomyopathy

Modifier genes contribute significantly to our understanding of pathophysiology in human diseases; however, effective approaches to identify modifier genes are still lacking. Here, we aim to develop a rapid F0-based genetic assay in adult zebrafish using the bag3 gene knockout (bag3e2/e2) cardiomyopathy model as a paradigm. First, by utilizing a classic genetic breeding approach, we identified dnajb6b as a deleterious modifier gene for bag3 cardiomyopathy. Next, we established an F0-based genetic assay in adult zebrafish through injection of predicted microhomology-mediated end joining (MMEJ)-inducing single guide RNA/Cas9 protein complex. We showed that effective gene knockdown is maintained in F0 adult fish, enabling recapitulation of both salutary modifying effects of the mtor haploinsufficiency and deleterious modifying effects of the dnajb6b gene on bag3 cardiomyopathy. We finally deployed the F0-based genetic assay to screen differentially expressed genes in the bag3 cardiomyopathy model. As a result, myh9b was identified as a novel modifier gene for bag3 cardiomyopathy. Together, these data prove the feasibility of an F0 adult zebrafish-based genetic assay that can be effectively used to discover modifier genes for inherited cardiomyopathy.

Loss of RREB1 in pancreatic beta cells reduces cellular insulin content and affects endocrine cell gene expression

AIMS/HYPOTHESIS

Genome-wide studies have uncovered multiple independent signals at the RREB1 locus associated with altered type 2 diabetes risk and related glycaemic traits. However, little is known about the function of the zinc finger transcription factor Ras-responsive element binding protein 1 (RREB1) in glucose homeostasis or how changes in its expression and/or function influence diabetes risk.

METHODS

A zebrafish model lacking rreb1a and rreb1b was used to study the effect of RREB1 loss in vivo. Using transcriptomic and cellular phenotyping of a human beta cell model (EndoC-βH1) and human induced pluripotent stem cell (hiPSC)-derived beta-like cells, we investigated how loss of RREB1 expression and activity affects pancreatic endocrine cell development and function. Ex vivo measurements of human islet function were performed in donor islets from carriers of RREB1 type 2 diabetes risk alleles.

RESULTS

CRISPR/Cas9-mediated loss of rreb1a and rreb1b function in zebrafish supports an in vivo role for the transcription factor in beta cell mass, beta cell insulin expression and glucose levels. Loss of RREB1 also reduced insulin gene expression and cellular insulin content in EndoC-βH1 cells and impaired insulin secretion under prolonged stimulation. Transcriptomic analysis of RREB1 knockdown and knockout EndoC-βH1 cells supports RREB1 as a novel regulator of genes involved in insulin secretion. In vitro differentiation of RREB1KO/KO hiPSCs revealed dysregulation of pro-endocrine cell genes, including RFX family members, suggesting that RREB1 also regulates genes involved in endocrine cell development. Human donor islets from carriers of type 2 diabetes risk alleles in RREB1 have altered glucose-stimulated insulin secretion ex vivo, consistent with a role for RREB1 in regulating islet cell function.

CONCLUSIONS/INTERPRETATION

Together, our results indicate that RREB1 regulates beta cell function by transcriptionally regulating the expression of genes involved in beta cell development and function.

New CRISPR Technology for Creating Cell Models of Lipoprotein Assembly and Secretion

PURPOSE OF REVIEW

This review is aimed at providing an overview of new developments in gene editing technology, including examples of how this technology has been used to develop cell models for studying the effects of gene ablation or missense mutations on lipoprotein assembly and secretion.

RECENT FINDINGS

CRISPR/Cas9-mediated gene editing is superior to other technologies because of its ease, sensitivity, and low off-target effects. This technology has been used to study the importance of microsomal triglyceride transfer protein in the assembly and secretion of apolipoprotein B-containing lipoproteins, as well as to establish causal effects of APOB gene missense mutations on lipoprotein assembly and secretion. CRISPR/Cas9 technology is anticipated to provide unprecedented flexibility in studying protein structure and function in cells and animals and to yield mechanistic insights into variants in the human genome.

A CRISPR-Cas9-mediated F0 screen to identify pro-regenerative genes in the zebrafish retinal pigment epithelium

Ocular diseases resulting in death of the retinal pigment epithelium (RPE) lead to vision loss and blindness. There are currently no FDA-approved strategies to restore damaged RPE cells. Stimulating intrinsic regenerative responses within damaged tissues has gained traction as a possible mechanism for tissue repair. Zebrafish possess remarkable regenerative abilities, including within the RPE; however, our understanding of the underlying mechanisms remains limited. Here, we conducted an F0 in vivo CRISPR-Cas9-mediated screen of 27 candidate RPE regeneration genes. The screen involved injection of a ribonucleoprotein complex containing three highly mutagenic guide RNAs per target gene followed by PCR-based genotyping to identify large intragenic deletions and MATLAB-based automated quantification of RPE regeneration. Through this F0 screening pipeline, eight positive and seven negative regulators of RPE regeneration were identified. Further characterization of one candidate, cldn7b, revealed novel roles in regulating macrophage/microglia infiltration after RPE injury and in clearing RPE/pigment debris during late-phase RPE regeneration. Taken together, these data support the utility of targeted F0 screens for validating pro-regenerative factors and reveal novel factors that could regulate regenerative responses within the zebrafish RPE.

Establishment of an Integrated CRISPR/Cas9 Plasmid System for Simple and Efficient Genome Editing in Medaka In Vitro and In Vivo

Although CRISPR/Cas9 has been used in gene manipulation of several fish species in vivo, its application in fish cultured cells is still challenged and limited. In this study, we established an integrated CRISPR/Cas9 plasmid system and evaluated its efficiency of gene knock-out or knock-in at a specific site in medaka (Oryzias latipes) in vitro and in vivo. By using the enhanced green fluorescent protein reporter plasmid pGNtsf1, we demonstrate that pCas9-U6sgRNA driven by endogenous U6 promoter (pCas9-mU6sgRNA) mediated very high gene editing efficiency in medaka cultured cells, but not by exogenous U6 promoters. After optimizing the conditions, the gene editing efficiencies of eight sites targeting for four endogenous genes were calculated, and the highest was up to 94% with no detectable off-target. By one-cell embryo microinjection, pCas9-mU6sgRNA also mediated efficient gene knock-out in vivo. Furthermore, pCas9-mU6sgRNA efficiently mediated gene knock-in at a specific site in medaka cultured cells as well as embryos. Collectively, our study demonstrates that the genetic relationship of U6 promoter is critical to gene editing efficiency in medaka cultured cells, and a simple and efficient system for medaka genome editing in vitro and in vivo has been established. This study provides an insight into other fish genome editing and promotes gene functional analysis.

Improved Genome Editing in the Ascidian Ciona with CRISPR/Cas9 and TALEN

The ascidian Ciona intestinalis type A (or Ciona robusta) is an important organism for elucidating the mechanisms that make the chordate body plan. CRISPR/Cas9 and TAL effector nuclease (TALEN) are widely used to quickly address genetic functions in Ciona. Our previously reported method of CRISPR/Cas9-mediated mutagenesis in this animal has inferior mutation rates compared to those of TALENs. We here describe an updated way to effectively mutate genes with CRISPR/Cas9 in Ciona. Although the construction of TALENs is much more laborious than that of CRISPR/Cas9, this technique is useful for tissue-specific knockouts that are not easy even by the optimized CRISPR/Cas9 method.

Transcription factors underlying photoreceptor diversity

During development, retinal progenitors navigate a complex landscape of fate decisions to generate the major cell classes necessary for proper vision. Transcriptional regulation is critical to generate diversity within these major cell classes. Here, we aim to provide the resources and techniques required to identify transcription factors necessary to generate and maintain diversity in photoreceptor subtypes, which are critical for vision. First, we generate a key resource: a high-quality and deep transcriptomic profile of each photoreceptor subtype in adult zebrafish. We make this resource openly accessible, easy to explore, and have integrated it with other currently available photoreceptor transcriptomic datasets. Second, using our transcriptomic profiles, we derive an in-depth map of expression of transcription factors in photoreceptors. Third, we use efficient CRISPR-Cas9 based mutagenesis to screen for null phenotypes in F0 larvae (F0 screening) as a fast, efficient, and versatile technique to assess the involvement of candidate transcription factors in the generation of photoreceptor subtypes. We first show that known phenotypes can be easily replicated using this method: loss of S cones in foxq2 mutants and loss of rods in nr2e3 mutants. We then identify novel functions for the transcription factor Tbx2, demonstrating that it plays distinct roles in controlling the generation of all photoreceptor subtypes within the retina. Our study provides a roadmap to discover additional factors involved in this process. Additionally, we explore four transcription factors of unknown function (Skor1a, Sall1a, Lrrfip1a, and Xbp1), and find no evidence for their involvement in the generation of photoreceptor subtypes. This dataset and screening method will be a valuable way to explore the genes involved in many other essential aspects of photoreceptor biology.

Oxoeicosanoid signaling mediates early antimicrobial defense in zebrafish

5-oxoETE is a bioactive lipid derived from arachidonic acid generated when phospholipase A₂ activation coincides with oxidative stress. Through its G protein-coupled receptor OXER1, pure 5-oxoETE is a potent leukocyte chemoattractant. Yet, its physiological function has remained elusive owing to the unusual OXER1 conservation pattern. OXER1 is conserved from fish to primates but not in rodents, precluding genetic loss-of-function studies in mouse. To determine its physiological role, we combine transcriptomic, lipidomic, and intravital imaging assays with genetic perturbations of the OXER1 ortholog hcar1-4 in zebrafish. Pseudomonas aeruginosa infection induces the synthesis of 5-oxoETE and its receptor, along with other inflammatory pathways. Hcar1-4 deletion attenuates neutrophil recruitment and decreases post-infection survival, which could be rescued by ectopic expression of hcar1-4 or human OXER1. By revealing 5-oxoETE as dominant lipid regulator of the early antimicrobial response in a non-rodent vertebrate, our work expands the current, rodent-centric view of early inflammation.

Morphological and sensorimotor phenotypes in a zebrafish CHARGE syndrome model are domain-dependent

CHARGE syndrome is a heterogeneous disorder characterized by a spectrum of defects affecting multiple tissues and behavioral difficulties such as autism, attention-deficit/hyperactivity disorder, obsessive-compulsive disorder, anxiety, and sensory deficits. Most CHARGE cases arise from de novo, loss-of-function mutations in chromodomain-helicase-DNA-binding-protein-7 (CHD7). CHD7 is required for processes such as neuronal differentiation and neural crest cell migration, but how CHD7 affects neural circuit function to regulate behavior is unclear. To investigate the pathophysiology of behavioral symptoms in CHARGE, we established a mutant chd7 zebrafish line that recapitulates multiple CHARGE phenotypes including ear, cardiac, and craniofacial defects. Using a panel of behavioral assays, we found that chd7 mutants have specific auditory and visual behavior deficits that are independent of defects in sensory structures. Mauthner cell-dependent short-latency acoustic startle responses are normal in chd7 mutants, while Mauthner-independent long-latency responses are reduced. Responses to sudden decreases in light are also reduced in mutants, while responses to sudden increases in light are normal, suggesting that the retinal OFF pathway may be affected. Furthermore, by analyzing multiple chd7 alleles we observed that the penetrance of morphological and behavioral phenotypes is influenced by genetic background but that it also depends on the mutation location, with a chromodomain mutation causing the highest penetrance. This pattern is consistent with analysis of a CHARGE patient dataset in which symptom penetrance was highest in subjects with mutations in the CHD7 chromodomains. These results provide new insight into the heterogeneity of CHARGE and will inform future work to define CHD7-dependent neurobehavioral mechanisms.

De novo PAM generation to reach initially inaccessible target sites for base editing

Base editing by CRISPR crucially depends on the presence of a protospacer adjacent motif (PAM) at the correct distance from the editing site. Here, we present and validate an efficient one-shot approach termed 'inception' that expands the editing range. This is achieved by sequential, combinatorial base editing: de novo generated synonymous, non-synonymous or intronic PAM sites facilitate subsequent base editing at nucleotide positions that were initially inaccessible, further opening the targeting range of highly precise editing approaches. We demonstrate the applicability of the inception concept in medaka (Oryzias latipes) in three settings: loss of function, by introducing a pre-termination STOP codon in the open reading frame of oca2; locally confined multi-codon changes to generate allelic variants with different phenotypic severity in kcnh6a; and the removal of a splice acceptor site by targeting intronic sequences of rx3. Using sequentially acting base editors in the described combinatorial approach expands the number of accessible target sites by 65% on average. This allows the use of well-established tools with NGG PAM recognition for the establishment of thus far unreachable disease models, for hypomorphic allele studies and for efficient targeted mechanistic investigations in a precise and predictable manner.

Myelination generates aberrant ultrastructure that is resolved by microglia

To enable rapid propagation of action potentials, axons are ensheathed by myelin, a multilayered insulating membrane formed by oligodendrocytes. Most of the myelin is generated early in development, resulting in the generation of long-lasting stable membrane structures. Here, we explored structural and dynamic changes in central nervous system myelin during development. To achieve this, we performed an ultrastructural analysis of mouse optic nerves by serial block face scanning electron microscopy (SBF-SEM) and confocal time-lapse imaging in the zebrafish spinal cord. We found that myelin undergoes extensive ultrastructural changes during early postnatal development. Myelin degeneration profiles were engulfed and phagocytosed by microglia using exposed phosphatidylserine as one "eat me" signal. In contrast, retractions of entire myelin sheaths occurred independently of microglia and involved uptake of myelin by the oligodendrocyte itself. Our findings show that the generation of myelin early in development is an inaccurate process associated with aberrant ultrastructural features that require substantial refinement.

A Critical Review of Zebrafish Neurological Disease Models-1. The Premise: Neuroanatomical, Cellular and Genetic Homology and Experimental Tractability

The last decade has seen a dramatic rise in the number of genes linked to neurological disorders, necessitating new models to explore underlying mechanisms and to test potential therapies. Over a similar period, many laboratories adopted zebrafish as a tractable model for studying brain development, defining neural circuits and performing chemical screens. Here we discuss strengths and limitations of using the zebrafish system to model neurological disorders. The underlying premise for many disease models is the high degree of homology between human and zebrafish genes, coupled with the conserved vertebrate Bauplan and repertoire of neurochemical signaling molecules. Yet, we caution that important evolutionary divergences often limit the extent to which human symptoms can be modeled meaningfully in zebrafish. We outline advances in genetic technologies that allow human mutations to be reproduced faithfully in zebrafish. Together with methods that visualize the development and function of neuronal pathways at the single cell level, there is now an unprecedented opportunity to understand how disease-associated genetic changes disrupt neural circuits, a level of analysis that is ideally suited to uncovering pathogenic changes in human brain disorders.

CRISPR-Cas9-Mediated Genome Modifications in Zebrafish

CRISPR-Cas9 genome editing technology has been successfully applied to generate various genetic modifications in zebrafish. The CRISPR-Cas9 system, which originally consisted of three components, CRISPR RNA (crRNA), trans-activating crRNA (tracrRNA), and Cas9, efficiently induces DNA double-strand breaks (DSBs) at targeted genomic loci, often resulting in frameshift-mediated target gene disruption (knockout). However, it remains difficult to perform the targeted integration of exogenous DNA fragments (knock-in) with CRISPR-Cas9. DSBs can be restored through DNA repair mechanisms, such as nonhomologous end joining (NHEJ), microhomology-mediated end joining (MMEJ), and homology-directed repair (HDR). One of our two research groups established a method for the precise MMEJ-mediated targeted integrations of exogenous genes containing homologous microhomology sequences flanking a targeted genomic locus in zebrafish. The other group recently developed a method for knocking in ~200 nt sequences encoding composite tags using long single-stranded DNA (ssDNA) donors. This chapter summarizes the CRISPR-Cas9-mediated genome modification strategy in zebrafish.

Targeted Manipulation of Histone Modification in Medaka Embryos

Recent development of targeted manipulation of histone modification enables us to experimentally and directly test the functional relevance of histone modifications accumulated at specific genomic regions. In particular, dCas9 epigenome editing has been widely used for site-specific manipulation of epigenetic modification. Here, we describe how to apply dCas9 epigenome editing in fish (medaka, Oryzias latipes) embryos and how to analyze induced changes in histone modification.

Rapid Testing of Gene Function in Axonal Regeneration After Spinal Cord Injury Using Larval Zebrafish

Larval zebrafish show axonal regrowth over a complex spinal injury site and recovery of function within days after injury. Here we describe a simple protocol to disrupt gene function in this model using acute injections of highly active synthetic gRNAs to rapidly detect loss-of-function phenotypes without the need for breeding.

Dissecting cell identity via network inference and in silico gene perturbation

Cell identity is governed by the complex regulation of gene expression, represented as gene-regulatory networks¹. Here we use gene-regulatory networks inferred from single-cell multi-omics data to perform in silico transcription factor perturbations, simulating the consequent changes in cell identity using only unperturbed wild-type data. We apply this machine-learning-based approach, CellOracle, to well-established paradigms-mouse and human haematopoiesis, and zebrafish embryogenesis-and we correctly model reported changes in phenotype that occur as a result of transcription factor perturbation. Through systematic in silico transcription factor perturbation in the developing zebrafish, we simulate and experimentally validate a previously unreported phenotype that results from the loss of noto, an established notochord regulator. Furthermore, we identify an axial mesoderm regulator, lhx1a. Together, these results show that CellOracle can be used to analyse the regulation of cell identity by transcription factors, and can provide mechanistic insights into development and differentiation.

A minimally invasive fin scratching protocol for fast genotyping and early selection of zebrafish embryos

Current genetic modification and phenotyping methods in teleost fish allow detailed investigation of vertebrate mechanisms of development, modeling of specific aspects of human diseases and efficient testing of drugs at an organ/organismal level in an unparalleled fast and large-scale mode. Fish-based experimental approaches have boosted the in vivo verification and implementation of scientific advances, offering the quality guaranteed by animal models that ultimately benefit human health, and are not yet fully replaceable by even the most sophisticated in vitro alternatives. Thanks to highly efficient and constantly advancing genetic engineering as well as non-invasive phenotyping methods, the small zebrafish is quickly becoming a popular alternative to large animals' experimentation. This approach is commonly associated to invasive procedures and increased burden. Here, we present a rapid and minimally invasive method to obtain sufficient genomic material from single zebrafish embryos by simple and precise tail fin scratching that can be robustly used for at least two rounds of genotyping already from embryos within 48 h of development. The described protocol betters currently available methods (such as fin clipping), by minimizing the relative animal distress associated with biopsy at later or adult stages. It allows early selection of embryos with desired genotypes for strategizing culturing or genotype-phenotype correlation experiments, resulting in a net reduction of "surplus" animals used for mutant line generation.

Deficiency of Adipose Triglyceride Lipase Induces Metabolic Syndrome and Cardiomyopathy in Zebrafish

Lipid metabolism dysfunction is related to clinical disorders including obesity, cancer, liver steatosis, and cardiomyopathy. Impaired lipolytic enzymes result in altered release of free fatty acids. The dramatic change in dyslipidemia is important in lipotoxic cardiomyopathy. Adipose triglyceride lipase (ATGL) catalyzes the lipolysis of triacylglycerol to reduce intramyocardial triglyceride levels in the heart and improve myocardial function. We examined the role of ATGL in metabolic cardiomyopathy by developing an Atgl knockout (ALKO) zebrafish model of metabolic cardiomyopathy disease by continuously expressing CRISPR/Cas9 protein and atgl gene guide RNAs (gRNAs). The expressed Cas9 protein bound to four gRNAs targeting the atgl gene locus, facilitating systemic gene KO. Ablation of Atgl interfered with lipid metabolism, which induced hyperlipidemia and hyperglycemia. ALKO adults and embryos displayed hypertrophic hearts. ALKO presented a typical dilated cardiomyopathy profile with a remarkable reduction in four sarcomere genes (myosin heavy chain 7-like, actin alpha cardiac muscle 1b, myosin binding protein C3, and troponin T type 2a) and two Ca²⁺ handling regulator genes (tropomyosin 4b and ATPase sarcoplasmic/endoplasmic reticulum Ca²⁺ transporting 2b). Immune cell infiltration in cardiac tissue of ALKO provided direct evidence of advanced metabolic cardiomyopathy. The presently described model could become a powerful tool to clarify the underlying mechanism between metabolic disorders and cardiomyopathies.

Modeling leukemia with zebrafish (Danio rerio): Towards precision medicine

Leukemia is a type of blood cancer characterized by high genetic heterogeneity and fatality. While chemotherapy remains the primary form of treatment for leukemia, its effectiveness was profoundly diminished by the genetic heterogeneity and cytogenetic abnormalities of leukemic cells. Therefore, there is an unmet need to develop precision medicine for leukemia with distinct genetic backgrounds. Zebrafish (Danio rerio), a freshwater fish with exceptional feasibility in genome editing, is a powerful tool for rapid human cancer modeling. In the past decades, zebrafish have been adopted in modeling human leukemia, exploring the molecular mechanisms of underlying genetic abnormalities, and discovering novel therapeutic agents. Although many recurrent mutations of leukemia have been modeled in zebrafish for pathological study and drug discovery, its great potential in leukemia modeling was not yet fully exploited, particularly in precision medicine. In this review, we evaluated the current zebrafish models of leukemia/pre-leukemia and genetic techniques and discussed the potential of zebrafish models with novel techniques, which may contribute to the development of zebrafish as a disease model for precision medicine in treating leukemia.

Macrophage NFATC2 mediates angiogenic signaling during mycobacterial infection

During mycobacterial infections, pathogenic mycobacteria manipulate both host immune and stromal cells to establish and maintain a productive infection. In humans, non-human primates, and zebrafish models of infection, pathogenic mycobacteria produce and modify the specialized lipid trehalose 6,6'-dimycolate (TDM) in the bacterial cell envelope to drive host angiogenesis toward the site of forming granulomas, leading to enhanced bacterial growth. Here, we use the zebrafish-Mycobacterium marinum infection model to define the signaling basis of the host angiogenic response. Through intravital imaging and cell-restricted peptide-mediated inhibition, we identify macrophage-specific activation of NFAT signaling as essential to TDM-mediated angiogenesis in vivo. Exposure of cultured human cells to Mycobacterium tuberculosis results in robust induction of VEGFA, which is dependent on a signaling pathway downstream of host TDM detection and culminates in NFATC2 activation. As granuloma-associated angiogenesis is known to serve bacterial-beneficial roles, these findings identify potential host targets to improve tuberculosis disease outcomes.

A robust pipeline for efficient knock-in of point mutations and epitope tags in zebrafish using fluorescent PCR based screening

BACKGROUND

Genome editing using CRISPR/Cas9 has become a powerful tool in zebrafish to generate targeted gene knockouts models. However, its use for targeted knock-in remains challenging due to inefficient homology directed repair (HDR) pathway in zebrafish, highlighting the need for efficient and cost-effective screening methods.  RESULTS: Here, we present our fluorescent PCR and capillary electrophoresis based screening approach for knock-in using a single-stranded oligodeoxynucleotide donor (ssODN) as a repair template for the targeted insertion of epitope tags, or single nucleotide changes to recapitulate pathogenic human alleles. For the insertion of epitope tags, we took advantage of the expected change in size of the PCR product. For point mutations, we combined fluorescent PCR with restriction fragment length polymorphism (RFLP) analysis to distinguish the fish with the knock-in allele. As a proof-of-principle, we present our data on the generation of fish lines with insertion of a FLAG tag at the tcnba locus, an HA tag at the gata2b locus, and a point mutation observed in Gaucher disease patients in the gba gene. Despite the low number of germline transmitting founders (1-5%), combining our screening methods with prioritization of founder fish by fin biopsies allowed us to establish stable knock-in lines by screening 12 or less fish per gene.

CONCLUSIONS

We have established a robust pipeline for the generation of zebrafish models with precise integration of small DNA sequences and point mutations at the desired sites in the genome. Our screening method is very efficient and easy to implement as it is PCR-based and only requires access to a capillary sequencer.

Rab35 governs apicobasal polarity through regulation of actin dynamics during sprouting angiogenesis

In early blood vessel development, trafficking programs, such as those using Rab GTPases, are tasked with delivering vesicular cargo with high spatiotemporal accuracy. However, the function of many Rab trafficking proteins remain ill-defined in endothelial tissue; therefore, their relevance to blood vessel development is unknown. Rab35 has been shown to play an enigmatic role in cellular behaviors which differs greatly between tissue-type and organism. Importantly, Rab35 has never been characterized for its potential contribution in sprouting angiogenesis; thus, our goal was to map Rab35’s primary function in angiogenesis. Our results demonstrate that Rab35 is critical for sprout formation; in its absence, apicobasal polarity is entirely lost in vitro and in vivo. To determine mechanism, we systematically explored established Rab35 effectors and show that none are operative in endothelial cells. However, we find that Rab35 partners with DENNd1c, an evolutionarily divergent guanine exchange factor, to localize to actin. Here, Rab35 regulates actin polymerization through limiting Rac1 and RhoA activity, which is required to set up proper apicobasal polarity during sprout formation. Our findings establish that Rab35 is a potent brake of actin remodeling during blood vessel development.

The promiscuous GTPase Rab35 has been shown to be involved in many important cellular functions. In this article, Francis et al. illustrate how Rab35 acts as a critical brake to actin remodeling during sprouting angiogenesis and how it is necessary for proper blood vessel development.

Teleost Fish and Organoids: Alternative Windows Into the Development of Healthy and Diseased Brains

The variety in the display of animals’ cognition, emotions, and behaviors, typical of humans, has its roots within the anterior-most part of the brain: the forebrain, giving rise to the neocortex in mammals. Our understanding of cellular and molecular events instructing the development of this domain and its multiple adaptations within the vertebrate lineage has progressed in the last decade. Expanding and detailing the available knowledge on regionalization, progenitors’ behavior and functional sophistication of the forebrain derivatives is also key to generating informative models to improve our characterization of heterogeneous and mechanistically unexplored cortical malformations. Classical and emerging mammalian models are irreplaceable to accurately elucidate mechanisms of stem cells expansion and impairments of cortex development. Nevertheless, alternative systems, allowing a considerable reduction of the burden associated with animal experimentation, are gaining popularity to dissect basic strategies of neural stem cells biology and morphogenesis in health and disease and to speed up preclinical drug testing. Teleost vertebrates such as zebrafish, showing conserved core programs of forebrain development, together with patients-derived in vitro 2D and 3D models, recapitulating more accurately human neurogenesis, are now accepted within translational workflows spanning from genetic analysis to functional investigation. Here, we review the current knowledge of common and divergent mechanisms shaping the forebrain in vertebrates, and causing cortical malformations in humans. We next address the utility, benefits and limitations of whole-brain/organism-based fish models or neuronal ensembles in vitro for translational research to unravel key genes and pathological mechanisms involved in neurodevelopmental diseases.

In vivo dissection of Rhoa function in vascular development using zebrafish

The small monomeric GTPase RHOA acts as a master regulator of signal transduction cascades by activating effectors of cellular signaling, including the Rho-associated protein kinases ROCK1/2. Previous in vitro cell culture studies suggest that RHOA can regulate many critical aspects of vascular endothelial cell (EC) biology, including focal adhesion, stress fiber formation, and angiogenesis. However, the specific in vivo roles of RHOA during vascular development and homeostasis are still not well understood. In this study, we examine the in vivo functions of RHOA in regulating vascular development and integrity in zebrafish. We use zebrafish RHOA-ortholog (rhoaa) mutants, transgenic embryos expressing wild type, dominant negative, or constitutively active forms of rhoaa in ECs, pharmacological inhibitors of RHOA and ROCK1/2, and Rock1 and Rock2a/b dgRNP-injected zebrafish embryos to study the in vivo consequences of RHOA gain- and loss-of-function in the vascular endothelium. Our findings document roles for RHOA in vascular integrity, developmental angiogenesis, and vascular morphogenesis in vivo, showing that either too much or too little RHOA activity leads to vascular dysfunction.

Massively targeted evaluation of therapeutic CRISPR off-targets in cells

Methods for sensitive and high-throughput evaluation of CRISPR RNA-guided nucleases (RGNs) off-targets (OTs) are essential for advancing RGN-based gene therapies. Here we report SURRO-seq for simultaneously evaluating thousands of therapeutic RGN OTs in cells. SURRO-seq captures RGN-induced indels in cells by pooled lentiviral OTs libraries and deep sequencing, an approach comparable and complementary to OTs detection by T7 endonuclease 1, GUIDE-seq, and CIRCLE-seq. Application of SURRO-seq to 8150 OTs from 110 therapeutic RGNs identifies significantly detectable indels in 783 OTs, of which 37 OTs are found in cancer genes and 23 OTs are further validated in five human cell lines by targeted amplicon sequencing. Finally, SURRO-seq reveals that thermodynamically stable wobble base pair (rG•dT) and free binding energy strongly affect RGN specificity. Our study emphasizes the necessity of thoroughly evaluating therapeutic RGN OTs to minimize inevitable off-target effects.

The embryonic zebrafish brain is seeded by a lymphatic-dependent population of mrc1+ microglia precursors

Microglia are the resident macrophages of the CNS that serve critical roles in brain construction. Although human brains contain microglia by 4 weeks gestation, an understanding of the earliest microglia that seed the brain during its development remains unresolved. Using time-lapse imaging in zebrafish, we discovered a mrc1a+ microglia precursor population that seeds the brain before traditionally described microglia. These early microglia precursors are dependent on lymphatic vasculature that surrounds the brain and are independent of pu1+ yolk sac-derived microglia. Single-cell RNA-sequencing datasets reveal Mrc1+ microglia in the embryonic brains of mice and humans. We then show in zebrafish that these early mrc1a+ microglia precursors preferentially expand during pathophysiological states in development. Taken together, our results identify a critical role of lymphatics in the microglia precursors that seed the early embryonic brain.