Clark, Karl J.

The zebrafish has emerged as an excellent model for cardiovascular research thanks to its ex-utero and rapid embryonic development, its embryonic transparent nature, and its capacity to survive in the absence of a functional cardiovascular system during the first week of development, which enables functional characterization of mutations that would otherwise induce lethality in traditional murine models. The aim of this project is to identify and characterize novel genes involved in zebrafish cardiovascular development using mutagenic gene trapping, a technique that generates random insertional mutations across the genome. We use the well-defined RP2 gene-breaking transposon system, which not only mutates, but also fluorescently tags the trapped gene product(s) (Clark, Nat Methods, 2011). RP2 also introduces loxP sites into the mutant locus, which can be used for Cre-mediated phenotype rescue by microinjection of Cre recombinase, or by crossing to tissue-specific Cre lines. From over 3000 RP2-injected embryos, 141 fish showed germline transmission. Among them, 51 expressed strong fluorescence in different tissues, including the heart, vessels, notochord, central nervous system (CNS), and eyes. Three cardiovascular lines were selected for phenotypic characterization and functional studies. RP2#C2 strain expressed fluorescence in the heart valves, CNS, eyes and pectoral fin buds. Inverse PCR in RP2#C2 demonstrated a trapped gene at meis4.1a, which encodes a homeobox transcription factor that has not been previously studied in zebrafish. RP2#121 strain expressed strong fluorescence in cardiac and skeletal muscles. RP2#91 strain showed expression in the vasculature and demonstrated a trapped gene at pdgfra. Homozygous RP2#91 mutants showed severe defects in the heart, blood flow, and other body parts including the head and musculature. We are currently creating a panel of tissue-specific Cre lines targeting tissues such as cardiomyocytes, endothelial cells and smooth muscle cells for spatiotemporal rescuing of the mutant phenotype. The generated zebrafish protein-trap lines are invaluable tools to annotate gene function, dissect the molecular mechanisms of cardiovascular development, and potentially serve as disease models.

Abstract Background GNB1 encodes a subunit of a heterotrimeric G‐protein complex that transduces intracellular signaling cascades. Disruptions to the gene have previously been shown to be embryonic lethal in knockout mice and to cause complex neurodevelopmental disorders in humans. To date, the majority of variants associated with disease in humans have been missense variants in exons 5‐7. Methods Genetic sequencing was performed on two patients presenting with complex neurological phenotypes including intellectual disability, hypotonia, and in one patient seizures. Reported variants were assessed using RNA sequencing and functional BRET/BiFC assays. Results A splice variant reported in patient 1 was confirmed to cause usage of a cryptic splice site leading to a truncated protein product. Patient 2 was reported to have a truncating variant. BRET and BiFC assays of both patient variants confirmed both were deficient in inducing GPCR‐induced G protein activation due to lack of dimer formation with the Gγ subunit. Conclusion Here, we report two patients with functionally confirmed loss of function variants in GNB1 and neurodevelopmental phenotypes including intellectual disability, hypotonia, and seizures in one patient. These results suggest haploinsufficiency of GNB1 is a mechanism for neurodevelopmental disorders in humans.

One key bottleneck in understanding the human genome is the relative under-characterization of 90% of protein coding regions. We report a collection of 1200 transgenic zebrafish strains made with the gene-break transposon (GBT) protein trap to simultaneously report and reversibly knockdown the tagged genes. Protein trap-associated mRFP expression shows previously undocumented expression of 35% and 90% of cloned genes at 2 and 4 days post-fertilization, respectively. Further, investigated alleles regularly show 99% gene-specific mRNA knockdown. Homozygous GBT animals in ryr1b , fras1 , tnnt2a , edar and hmcn1 phenocopied established mutants. 204 cloned lines trapped diverse proteins, including 64 orthologs of human disease-associated genes with 40 as potential new disease models. Severely reduced skeletal muscle Ca 2+ transients in GBT ryr1b homozygous animals validated the ability to explore molecular mechanisms of genetic diseases. This GBT system facilitates novel functional genome annotation towards understanding cellular and molecular underpinnings of vertebrate biology and human disease.

Efficient precision genome engineering requires high frequency and specificity of integration at the genomic target site. Here, we describe a set of resources to streamline reporter gene knock-ins in zebrafish and demonstrate the broader utility of the method in mammalian cells. Our approach uses short homology of 24–48 bp to drive targeted integration of DNA reporter cassettes by homology-mediated end joining (HMEJ) at high frequency at a double strand break in the targeted gene. Our vector series, pGTag (plasmids for Gene Tagging), contains reporters flanked by a universal CRISPR sgRNA sequence which enables in vivo liberation of the homology arms. We observed high rates of germline transmission (22–100%) for targeted knock-ins at eight zebrafish loci and efficient integration at safe harbor loci in porcine and human cells. Our system provides a straightforward and cost-effective approach for high efficiency gene targeting applications in CRISPR and TALEN compatible systems.

CRISPR and CRISPR-Cas effector proteins enable the targeting of DNA double-strand breaks to defined loci based on a variable length RNA guide specific to each effector. The guide RNAs are generally similar in size and form, consisting of a ∼20 nucleotide sequence complementary to the DNA target and an RNA secondary structure recognized by the effector. However, the effector proteins vary in protospacer adjacent motif requirements, nuclease activities, and DNA binding kinetics. Recently, ErCas12a, a new member of the Cas12a family, was identified in Eubacterium rectale . Here, we report the first characterization of ErCas12a activity in zebrafish and expand on previously reported activity in human cells. Using a fluorescent reporter system, we show that CRISPR-ErCas12a elicits strand annealing mediated DNA repair more efficiently than CRISPR-Cas9. Further, using our previously reported gene targeting method that utilizes short homology, GeneWeld, we demonstrate the use of CRISPR-ErCas12a to integrate reporter alleles into the genomes of both zebrafish and human cells. Together, this work provides methods for deploying an additional CRISPR-Cas system, thus increasing the flexibility researchers have in applying genome engineering technologies.

AbstractDiffuse idiopathic skeletal hyperostosis (DISH) is a disorder principally characterized by calcification and ossification of spinal ligaments and entheses. Fibrodysplasia ossificans progressiva (FOP) is a rare autosomal dominant disabling disorder characterized by progressive ossification of skeletal muscle, fascia, tendons, and ligaments. These conditions manifest phenotypic overlap in the ossification of tendons and ligaments. We describe herein a patient with DISH, exhibiting heterotopic ossification of the posterior longitudinal ligament where clinical whole exome sequencing identified a variant within ACVR1, a gene implicated in FOP. This variant, p.K400E, is a novel variant, not identified previously, and occurs in a highly conserved region across orthologs. We used sequence‐based predicative algorithms, molecular modeling, and molecular dynamics simulations, to test the potential for p.K400E to alter the structure and dynamics of ACVR1. We applied the same modeling and simulation methods to established FOP variants, to identify the detailed effects that they have on the ACVR1 protein, as well as to act as positive controls against which the effects of p.K400E could be evaluated. Our in silico molecular analyses support p.K400E as altering the behavior of ACVR1. In addition, functional testing to measure the effect of this variant on BMP‐pSMAD 1/5/8 target genes was carried out which revealed this variant to cause increased ID1 and Msx2 expression compared with the wild‐type receptor. This analysis supports the potential for the variant of uncertain significance to contribute to the patient's phenotype.

Abstract The discovery and development of DNA-editing nucleases (Zinc Finger Nucleases, TALENs, CRISPR/Cas systems) has given scientists the ability to precisely engineer or edit genomes as never before. Several different platforms, protocols and vectors for precision genome editing are now available, leading to the development of supporting web-based software. Here we present the Gene Sculpt Suite (GSS), which comprises three tools: (i) GTagHD, which automatically designs and generates oligonucleotides for use with the GeneWeld knock-in protocol; (ii) MEDJED, a machine learning method, which predicts the extent to which a double-stranded DNA break site will utilize the microhomology-mediated repair pathway; and (iii) MENTHU, a tool for identifying genomic locations likely to give rise to a single predominant microhomology-mediated end joining allele (PreMA) repair outcome. All tools in the GSS are freely available for download under the GPL v3.0 license and can be run locally on Windows, Mac and Linux systems capable of running R and/or Docker. The GSS is also freely available online at www.genesculpt.org.

When vertebrates face acute stressors, their bodies rapidly undergo a repertoire of physiological and behavioral adaptations, which is termed the stress response. Rapid changes in heart rate and blood glucose levels occur via the interaction of glucocorticoids and their cognate receptors following hypothalamic‐pituitary‐adrenal axis activation. These physiological changes are observed within minutes of encountering a stressor and the rapid time domain rules out genomic responses that require gene expression changes. Although behavioral changes corresponding to physiological changes are commonly observed, it is not clearly understood to what extent hypothalamic‐pituitary‐adrenal axis activation dictates adaptive behavior. We hypothesized that rapid locomotor response to acute stressors in zebrafish requires hypothalamic‐pituitary‐interrenal (HPI) axis activation. In teleost fish, interrenal cells are functionally homologous to the adrenocortical layer. We derived eight frameshift mutants in genes involved in HPI axis function: two mutants in exon 2 ofmc2r(adrenocorticotropic hormone receptor), five in exon 2 or 5 ofnr3c1(glucocorticoid receptor [GR]) and two in exon 2 ofnr3c2(mineralocorticoid receptor [MR]). Exposing larval zebrafish to mild environmental stressors, acute changes in salinity or light illumination, results in a rapid locomotor response. We show that this locomotor response requires a functioning HPI axis via the action ofmc2rand the canonical GR encoded bynr3c1gene, but not MR (nr3c2). Our rapid behavioral assay paradigm based on HPI axis biology can be used to screen for genetic and environmental modifiers of the hypothalamic‐pituitary‐adrenal axis and to investigate the effects of corticosteroids and their cognate receptor interactions on behavior.