Impairments of cerebellar structure and function in a zebrafish KO of neuropsychiatric risk gene znf536.

Genetic variants in ZNF536 contribute to the risk for neuropsychiatric disorders such as schizophrenia, autism, and others. The role of this putative transcriptional repressor in brain development and function is, however, largely unknown. We generated znf536 knockout (KO) zebrafish and studied their behavior, brain anatomy, and brain function. Larval KO zebrafish showed a reduced ability to compete for food, resulting in decreased total body length and size. This phenotype can be rescued by segregating the homozygous KO larvae from their wild-type and heterozygous siblings, enabling studies of adult homozygous KO animals. In adult KO zebrafish, we observed significant reductions in anxiety-like behavior and social interaction. These znf536 KO zebrafish have decreased cerebellar volume, corresponding to decreased populations of specific neuronal cells, especially in the valvular cerebelli (Va). Finally, using a Tg[mbp:mgfp] line, we identified a previously undetected myelin structure located bilaterally within the Va, which also displayed a reduction in volume and disorganization in KO zebrafish. These findings indicate an important role for ZNF536 in brain development and implicate the cerebellum in the pathophysiology of neuropsychiatric disorders.

Mitochondrial proteins encoded by the 22q11.2 neurodevelopmental locus regulate neural stem and progenitor cell proliferation.

Microdeletion of a 3Mb region encompassing 45 protein-coding genes at chromosome 22q11.2 (22q11.2DS) predisposes individuals to multiple neurodevelopmental disorders and is one of the greatest genetic risk factors for schizophrenia. Defective mitochondrial function has been hypothesized to contribute to 22q11.2DS pathogenesis; however, which of the six mitochondrial genes contribute to neurodevelopmental phenotypes and their underlying mechanisms remain unresolved. To systematically test 22q11.2DS genes for functional roles in neurodevelopment and behavior, we generated genetic mutants for each of the 37 conserved zebrafish orthologs and performed high throughput behavioral phenotyping using seven behavioral assays. Through this unbiased approach, we identified five single-gene mutants with partially overlapping behavioral phenotypes. Two of these genes, mrpl40 and prodha, encode for mitochondrial proteins and, similar to what we observed in mrpl40 and prodha mutants, pharmacologic inhibition of mitochondrial function during development results in microcephaly. Single mutant analysis shows that both mrpl40 and prodha mutants display aberrant neural stem and progenitor cell proliferation, with each gene regulating distinct cell populations. Finally, double mutants for both mrpl40 and prodha display aggravated behavioral phenotypes and neural stem and progenitor cell analysis reveals a previously unrecognized partially redundant role for mrpl40 and prodha in regulating radial glia-like cell proliferation. Combined, our results demonstrate a critical role for mitochondrial function in neural stem and progenitor cell populations in the developing vertebrate brain and provide compelling evidence that mitochondrial dysfunction during neurodevelopment is linked to brain volume and behavioral phenotypes observed in models of 22q11.2DS.

Behavioral analysis through the lifespan of disc1 mutant zebrafish identifies defects in sensorimotor transformation.

DISC1 is a genetic risk factor for multiple psychiatric disorders. Compared to the dozens of murine Disc1 models, there is a paucity of zebrafish disc1 models-an organism amenable to high-throughput experimentation. We conducted the longitudinal neurobehavioral analysis of disc1 mutant zebrafish across key stages of life. During early developmental stages, disc1 mutants exhibited abrogated behavioral responses to sensory stimuli across multiple testing platforms. Moreover, during exposure to an acoustic sensory stimulus, loss of disc1 resulted in the abnormal activation of neurons in the pallium, cerebellum, and tectum-anatomical sites involved in the integration of sensory perception and motor control. In adulthood, disc1 mutants exhibited sexually dimorphic reduction in anxiogenic behavior in novel paradigms. Together, these findings implicate disc1 in sensorimotor processes and the genesis of anxiogenic behaviors, which could be exploited for the development of novel treatments in addition to investigating the biology of sensorimotor transformation in the context of disc1 deletion.

Multiple adaptations underly co-option of a CRISPR surveillance complex for RNA-guided DNA transposition.

CRISPR-associated transposons (CASTs) are natural RNA-directed transposition systems. We demonstrate that transposon protein TniQ plays a central role in promoting R-loop formation by RNA-guided DNA-targeting modules. TniQ residues, proximal to CRISPR RNA (crRNA), are required for recognizing different crRNA categories, revealing an unappreciated role of TniQ to direct transposition into different classes of crRNA targets. To investigate adaptations allowing CAST elements to utilize attachment sites inaccessible to CRISPR-Cas surveillance complexes, we compared and contrasted PAM sequence requirements in both I-F3b CAST and I-F1 CRISPR-Cas systems. We identify specific amino acids that enable a wider range of PAM sequences to be accommodated in I-F3b CAST elements compared with I-F1 CRISPR-Cas, enabling CAST elements to access attachment sites as sequences drift and evade host surveillance. Together, this evidence points to the central role of TniQ in facilitating the acquisition of CRISPR effector complexes for RNA-guided DNA transposition.

A single-cell transcriptome atlas of the maturing zebrafish telencephalon.

The zebrafish telencephalon is composed of highly specialized subregions that regulate complex behaviors such as learning, memory, and social interactions. The transcriptional signatures of the neuronal cell types in the telencephalon and the timeline of their emergence from larva to adult remain largely undescribed. Using an integrated analysis of single-cell transcriptomes of approximately 64,000 cells obtained from 6-day-postfertilization (dpf), 15-dpf, and adult telencephalon, we delineated nine main neuronal cell types in the pallium and eight in the subpallium and nominated novel marker genes. Comparing zebrafish and mouse neuronal cell types revealed both conserved and absent types and marker genes. Mapping of cell types onto a spatial larval reference atlas created a resource for anatomical and functional studies. Using this multiage approach, we discovered that although most neuronal subtypes are established early in the 6-dpf fish, some emerge or expand in number later in development. Analyzing the samples from each age separately revealed further complexity in the data, including several cell types that expand substantially in the adult forebrain and do not form clusters at the larval stages. Together, our work provides a comprehensive transcriptional analysis of the cell types in the zebrafish telencephalon and a resource for dissecting its development and function.

Mitochondrial genes in the 22q11.2 deleted region regulate neural stem and progenitor cell proliferation.

Microdeletion of a 3Mbp region encompassing 45 protein-coding genes at chromosome 22q11.2 (22q11.2DS) predisposes to multiple neurodevelopmental disorders and is one of the greatest genetic risk factors for schizophrenia. Defective mitochondrial function has been hypothesized to contribute to 22q11.2DS pathogenesis; however, which of the six mitochondrial genes contribute to neurodevelopmental phenotypes and their underlying mechanisms remain unresolved. To systematically test 22q11.2DS genes for functional roles in neurodevelopment and behavior, we generated genetic mutants for each of the 37 conserved zebrafish orthologs and performed high throughput behavioral phenotyping using seven behavioral assays. Through this unbiased approach, we identified five single-gene mutants with partially overlapping behavioral phenotypes. Two of these genes, mrpl40 and prodha , encode for mitochondrial proteins and, similar to what we observed in mrpl40 and prodha mutants, pharmacologic inhibition of mitochondrial function during development results in microcephaly. Finally, we show that both mrpl40 and prodha mutants display neural stem and progenitor cell phenotypes, with each gene regulating different neural stem cell populations. Combined, our results demonstrate a critical role for mitochondrial function in neural stem and progenitor cell populations in the developing vertebrate brain and provide compelling evidence that mitochondrial dysfunction during neurodevelopment is linked to brain volume and behavioral phenotypes observed in models of 22q11.2DS.

Mosaic and non-mosaic protocadherin 19 mutation leads to neuronal hyperexcitability in zebrafish.

Epilepsy is one of the most common neurological disorders. The X-linked gene PCDH19 is associated with sporadic and familial epilepsy in humans, typically with early-onset clustering seizures and intellectual disability in females but not in so-called 'carrier' males, suggesting that mosaic PCDH19 expression is required to produce epilepsy. To characterize the role of loss of PCDH19 function in epilepsy, we generated zebrafish with truncating pcdh19 variants. Evaluating zebrafish larvae for electrophysiological abnormalities, we observed hyperexcitability phenotypes in both mosaic and non-mosaic pcdh19+/- and pcdh19-/- mutant larvae. Thus, we demonstrate that the key feature of epilepsy-network hyperexcitability-can be modeled effectively in zebrafish, even though overt spontaneous seizure-like swim patterns were not observed. Further, zebrafish with non-mosaic pcdh19 mutation displayed reduced numbers of inhibitory interneurons suggesting a potential cellular basis for the observed hyperexcitability. Our findings in both mosaic and non-mosaic pcdh19 mutant zebrafish challenge the prevailing theory that mosaicism governs all PCDH19-related phenotypes and point to interneuron-mediated mechanisms underlying these phenotypes.

ZebraShare: a new venue for rapid dissemination of zebrafish mutant data.

BACKGROUND: In the past decade, the zebrafish community has widely embraced targeted mutagenesis technologies, resulting in an abundance of mutant lines. While many lines have proven to be useful for investigating gene function, many have also shown no apparent phenotype, or phenotypes not of interest to the originating lab. In order for labs to document and share information about these lines, we have created ZebraShare as a new resource offered within ZFIN. METHODS: ZebraShare involves a form-based submission process generated by ZFIN. The ZebraShare interface (https://zfin.org/action/zebrashare) can be accessed on ZFIN under "Submit Data". Users download the Submission Workbook and complete the required fields, then submit the completed workbook with associated images and captions, generating a new ZFIN publication record. ZFIN curators add the submitted phenotype and mutant information to the ZFIN database, provide mapping information about mutations, and cross reference this information across the appropriate ZFIN databases. We present here examples of ZebraShare submissions, including phf21aa, kdm1a, ctnnd1, snu13a, and snu13b mutant lines. RESULTS: Users can find ZebraShare submissions by searching ZFIN for specific alleles or line designations, just as for alleles submitted through the normal process. We present several potential examples of submission types to ZebraShare including a phenotypic mutants, mildly phenotypic, and early lethal mutants. Mutants for kdm1a show no apparent skeletal phenotype, and phf21aa mutants show only a mild skeletal phenotype, yet these genes have specific human disease relevance and therefore may be useful for further studies. The p120-catenin encoding gene, ctnnd1, was knocked out to investigate a potential role in brain development or function. The homozygous ctnnd1 mutant disintegrates during early somitogenesis and the heterozygote has localized defects, revealing vital roles in early development. Two snu13 genes were knocked out to investigate a role in muscle formation. The snu13a;snu13b double mutant has an early embryonic lethal phenotype, potentially related to a proposed role in the core splicing complex. In each example, the mutants submitted to ZebraShare display phenotypes that are not ideally suited to their originating lab's project directions but may be of great relevance to other researchers. CONCLUSION: ZebraShare provides an opportunity for researchers to directly share information about mutant lines within ZFIN, which is widely used by the community as a central database of information about zebrafish lines. Submissions of alleles with a phenotypic or unexpected phenotypes is encouraged to promote collaborations, disseminate lines, reduce redundancy of effort and to promote efficient use of time and resources. We anticipate that as submissions to ZebraShare increase, they will help build an ultimately more complete picture of zebrafish genetics and development.

A Customizable Low-Cost System for Massively Parallel Zebrafish Behavioral Phenotyping.

High-throughput behavioral phenotyping is critical to genetic or chemical screening approaches. Zebrafish larvae are amenable to high-throughput behavioral screening because of their rapid development, small size, and conserved vertebrate brain architecture. Existing commercial behavioral phenotyping systems are expensive and not easily modified for new assays. Here, we describe a modular, highly adaptable, and low-cost system. Along with detailed assembly and operation instructions, we provide data acquisition software and a robust, parallel analysis pipeline. We validate our approach by analyzing stimulus response profiles in larval zebrafish, confirming prepulse inhibition phenotypes of two previously isolated mutants, and highlighting best practices for growing larvae prior to behavioral testing. Our new design thus allows rapid construction and streamlined operation of many large-scale behavioral setups with minimal resources and fabrication expertise, with broad applications to other aquatic organisms.

Phenotypic Landscape of Schizophrenia-Associated Genes Defines Candidates and Their Shared Functions.

Genomic studies have identified hundreds of candidate genes near loci associated with risk for schizophrenia. To define candidates and their functions, we mutated zebrafish orthologs of 132 human schizophrenia-associated genes. We created a phenotype atlas consisting of whole-brain activity maps, brain structural differences, and profiles of behavioral abnormalities. Phenotypes were diverse but specific, including altered forebrain development and decreased prepulse inhibition. Exploration of these datasets identified promising candidates in more than 10 gene-rich regions, including the magnesium transporter cnnm2 and the translational repressor gigyf2, and revealed shared anatomical sites of activity differences, including the pallium, hypothalamus, and tectum. Single-cell RNA sequencing uncovered an essential role for the understudied transcription factor znf536 in the development of forebrain neurons implicated in social behavior and stress. This phenotypic landscape of schizophrenia-associated genes prioritizes more than 30 candidates for further study and provides hypotheses to bridge the divide between genetic association and biological mechanism.

Kctd13 deletion reduces synaptic transmission via increased RhoA.

Copy-number variants of chromosome 16 region 16p11.2 are linked to neuropsychiatric disorders and are among the most prevalent in autism spectrum disorders. Of many 16p11.2 genes, Kctd13 has been implicated as a major driver of neurodevelopmental phenotypes. The function of KCTD13 in the mammalian brain, however, remains unknown. Here we delete the Kctd13 gene in mice and demonstrate reduced synaptic transmission. Reduced synaptic transmission correlates with increased levels of Ras homolog gene family, member A (RhoA), a KCTD13/CUL3 ubiquitin ligase substrate, and is reversed by RhoA inhibition, suggesting increased RhoA as an important mechanism. In contrast to a previous knockdown study, deletion of Kctd13 or kctd13 does not increase brain size or neurogenesis in mice or zebrafish, respectively. These findings implicate Kctd13 in the regulation of neuronal function relevant to neuropsychiatric disorders and clarify the role of Kctd13 in neurogenesis and brain size. Our data also reveal a potential role for RhoA as a therapeutic target in disorders associated with KCTD13 deletion.

Internal guide RNA interactions interfere with Cas9-mediated cleavage.

The CRISPR/Cas system uses guide RNAs (gRNAs) to direct sequence-specific DNA cleavage. Not every gRNA elicits cleavage and the mechanisms that govern gRNA activity have not been resolved. Low activity could result from either failure to form a functional Cas9-gRNA complex or inability to recognize targets in vivo. Here we show that both phenomena influence Cas9 activity by comparing mutagenesis rates in zebrafish embryos with in vitro cleavage assays. In vivo, our results suggest that genomic factors such as CTCF inhibit mutagenesis. Comparing near-identical gRNA sequences with different in vitro activities reveals that internal gRNA interactions reduce cleavage. Even though gRNAs containing these structures do not yield cleavage-competent complexes, they can compete with active gRNAs for binding to Cas9. These results reveal that both genomic context and internal gRNA interactions can interfere with Cas9-mediated cleavage and illuminate previously uncharacterized features of Cas9-gRNA complex formation.

CHOPCHOP v2: a web tool for the next generation of CRISPR genome engineering.

In just 3 years CRISPR genome editing has transformed biology, and its popularity and potency continue to grow. New CRISPR effectors and rules for locating optimum targets continue to be reported, highlighting the need for computational CRISPR targeting tools to compile these rules and facilitate target selection and design. CHOPCHOP is one of the most widely used web tools for CRISPR- and TALEN-based genome editing. Its overarching principle is to provide an intuitive and powerful tool that can serve both novice and experienced users. In this major update we introduce tools for the next generation of CRISPR advances, including Cpf1 and Cas9 nickases. We support a number of new features that improve the targeting power, usability and efficiency of CHOPCHOP. To increase targeting range and specificity we provide support for custom length sgRNAs, and we evaluate the sequence composition of the whole sgRNA and its surrounding region using models compiled from multiple large-scale studies. These and other new features, coupled with an updated interface for increased usability and support for a continually growing list of organisms, maintain CHOPCHOP as one of the leading tools for CRISPR genome editing. CHOPCHOP v2 can be found at http://chopchop.cbu.uib.no.

Polq-Mediated End Joining Is Essential for Surviving DNA Double-Strand Breaks during Early Zebrafish Development.

Error-prone repair of DNA double-strand breaks (DSBs) has been postulated to occur through classical non-homologous end joining (NHEJ) in systems ranging from nematode somatic tissues to zebrafish embryos. Contrary to this model, we show that zebrafish embryos mutant for DNA polymerase theta (Polq), a critical component of alternative end joining (alt-EJ), cannot repair DSBs induced by CRISPR/Cas9 or ionizing radiation. In the absence of DSBs, polq mutants are phenotypically normal, but they do not survive mutagenesis and display dramatic differences in the mutation profiles compared with the wild-type. These results show that alt-EJ repair is essential and dominant during the early development of a vertebrate.

Computational Design of DNA-Binding Proteins.

Predicting the outcome of engineered and naturally occurring sequence perturbations to protein-DNA interfaces requires accurate computational modeling technologies. It has been well established that computational design to accommodate small numbers of DNA target site substitutions is possible. This chapter details the basic method of design used in the Rosetta macromolecular modeling program that has been successfully used to modulate the specificity of DNA-binding proteins. More recently, combining computational design and directed evolution has become a common approach for increasing the success rate of protein engineering projects. The power of such high-throughput screening depends on computational methods producing multiple potential solutions. Therefore, this chapter describes several protocols for increasing the diversity of designed output. Lastly, we describe an approach for building comparative models of protein-DNA complexes in order to utilize information from homologous sequences. These models can be used to explore how nature modulates specificity of protein-DNA interfaces and potentially can even be used as starting templates for further engineering.

Massively parallel determination and modeling of endonuclease substrate specificity.

We describe the identification and characterization of novel homing endonucleases using genome database mining to identify putative target sites, followed by high throughput activity screening in a bacterial selection system. We characterized the substrate specificity and kinetics of these endonucleases by monitoring DNA cleavage events with deep sequencing. The endonuclease specificities revealed by these experiments can be partially recapitulated using 3D structure-based computational models. Analysis of these models together with genome sequence data provide insights into how alternative endonuclease specificities were generated during natural evolution.

Efficient mutagenesis by Cas9 protein-mediated oligonucleotide insertion and large-scale assessment of single-guide RNAs.

The CRISPR/Cas9 system has been implemented in a variety of model organisms to mediate site-directed mutagenesis. A wide range of mutation rates has been reported, but at a limited number of genomic target sites. To uncover the rules that govern effective Cas9-mediated mutagenesis in zebrafish, we targeted over a hundred genomic loci for mutagenesis using a streamlined and cloning-free method. We generated mutations in 85% of target genes with mutation rates varying across several orders of magnitude, and identified sequence composition rules that influence mutagenesis. We increased rates of mutagenesis by implementing several novel approaches. The activities of poor or unsuccessful single-guide RNAs (sgRNAs) initiating with a 5' adenine were improved by rescuing 5' end homogeneity of the sgRNA. In some cases, direct injection of Cas9 protein/sgRNA complex further increased mutagenic activity. We also observed that low diversity of mutant alleles led to repeated failure to obtain frame-shift mutations. This limitation was overcome by knock-in of a stop codon cassette that ensured coding frame truncation. Our improved methods and detailed protocols make Cas9-mediated mutagenesis an attractive approach for labs of all sizes.

Progressive engineering of a homing endonuclease genome editing reagent for the murine X-linked immunodeficiency locus.

LAGLIDADG homing endonucleases (LHEs) are compact endonucleases with 20-22 bp recognition sites, and thus are ideal scaffolds for engineering site-specific DNA cleavage enzymes for genome editing applications. Here, we describe a general approach to LHE engineering that combines rational design with directed evolution, using a yeast surface display high-throughput cleavage selection. This approach was employed to alter the binding and cleavage specificity of the I-Anil LHE to recognize a mutation in the mouse Bruton tyrosine kinase (Btk) gene causative for mouse X-linked immunodeficiency (XID)-a model of human X-linked agammaglobulinemia (XLA). The required re-targeting of I-AniI involved progressive resculpting of the DNA contact interface to accommodate nine base differences from the native cleavage sequence. The enzyme emerging from the progressive engineering process was specific for the XID mutant allele versus the wild-type (WT) allele, and exhibited activity equivalent to WT I-AniI in vitro and in cellulo reporter assays. Fusion of the enzyme to a site-specific DNA binding domain of transcription activator-like effector (TALE) resulted in a further enhancement of gene editing efficiency. These results illustrate the potential of LHE enzymes as specific and efficient tools for therapeutic genome engineering.