Since the publication of the human reference genome, the identities of specific genes associated with human diseases are being discovered at an enormous rate. A central problem is that the biological activity of these genes is often unclear. Detailed investigations in vertebrate model organisms, typically mice, have been essential for understanding the activities of many orthologues of these disease-associated genes. Although gene-targeting approaches1-3 and phenotype analysis have led to a detailed understanding of nearly 6,000 protein-coding genes3,4, this number falls significantly short of all >22,000 mouse protein-coding genes5. Similarly, in zebrafish genetics, one-by-one gene studies using positional cloning6, insertional mutagenesis7-9, antisense morpholino oligonucleotides10, targeted re-sequencing11-13 and zinc finger and TAL endonucleases14-17 have made significant contributions to our understanding of the biological activity of vertebrate genes, but the number of genes studied again falls well short of the >26,000 zebrafish protein-coding genes18. Importantly, for both mice and zebrafish, none of these strategies is particularly suited to the rapid generation of knockouts in thousands of genes and the assessment of their biological activity. Enabled by a well-annotated zebrafish reference genome sequence18,19, high-throughput sequencing and efficient chemical mutagenesis, we describe an active project that aims to identify and phenotype disruptive mutations in every zebrafish protein-coding gene. Thus far we have identified potentially disruptive mutations in more than 38% of all known protein coding genes. We have developed a multi-allelic phenotyping scheme to efficiently assess the effects of each allele during embryogenesis and have analysed the phenotypic consequences of over 1000 alleles. All mutant alleles and data are available to the community and our phenotyping scheme is adaptable to phenotypic analysis beyond embryogenesis.

Pancreatic ductal adenocarcinoma (PDA) remains a lethal malignancy despite tremendous progress in its molecular characterization. Indeed, PDA tumors harbor four signature somatic mutations1–4, and a plethora of lower frequency genetic events of uncertain significance5. Here, we used Sleeping Beauty (SB) transposon-mediated insertional mutagenesis6,7 in a mouse model of pancreatic ductal preneoplasia8 to identify genes that cooperate with oncogenic KrasG12D to accelerate tumorigenesis and promote progression. Our screen revealed new candidates and confirmed the importance of many genes and pathways previously implicated in human PDA. Interestingly, the most commonly mutated gene was the X-linked deubiquitinase Usp9x, which was inactivated in over 50% of the tumors. Although prior work had attributed a pro-survival role to USP9X in human neoplasia9, we found instead that loss of Usp9x enhances transformation and protects pancreatic cancer cells from anoikis. Clinically, low USP9X protein and mRNA expression in PDA correlates with poor survival following surgery, and USP9X levels are inversely associated with metastatic burden in advanced disease. Furthermore, chromatin modulation with trichostatin A or 5-aza-2′-deoxycytidine elevates USP9X expression in human PDA cell lines to suggest a clinical approach for certain patients. The conditional deletion of Usp9x cooperated with KrasG12D to rapidly accelerate pancreatic tumorigenesis in mice, validating their genetic interaction. Therefore, we propose USP9X as a major new tumor suppressor gene with prognostic and therapeutic relevance in PDA.

Mouse ENU mutagenesis programmes have yielded a series of independent mutations on proximal chromosome 4 leading to dominant head-bobbing and circling behaviour due to truncations of the lateral semicircular canal of the inner ear. Here, we report the identification of mutations in the Chd7 gene in nine of these mutant alleles including six nonsense and three splice site mutations. The human CHD7 gene is known to be involved in CHARGE syndrome, which also shows inner ear malformations and a variety of other features with varying penetrance and appears to be due to frequent de novo mutation. We found widespread expression of Chd7 in early development of the mouse in organs affected in CHARGE syndrome including eye, olfactory epithelium, inner ear and vascular system. Closer inspection of heterozygous mutant mice revealed a range of defects with reduced penetrance, such as cleft palate, choanal atresia, septal defects of the heart, haemorrhages, prenatal death, vulva and clitoral defects and keratoconjunctivitis sicca. Many of these defects mimic the features of CHARGE syndrome. There were no obvious features of the gene that might make it more mutable than other genes. We conclude that the large number of mouse mutants and human de novo mutations may be due to the combination of the Chd7 gene being a large target and the fact that many heterozygous carriers of the mutations are viable individuals with a readily detectable phenotype.

Transcription is an essential component of basic cellular and developmental processes. However, early embryonic development occurs in the absence of transcription and instead relies upon maternal mRNAs and proteins deposited in the egg during oocyte maturation. Although the early zebrafish embryo is competent to transcribe exogenous DNA, factors present in the embryo maintain genomic DNA in a state that is incompatible with transcription. The cell cycles of the early embryo titrate out these factors, leading to zygotic transcription initiation, presumably in response to a change in genomic DNA chromatin structure to a state that supports transcription. To understand the molecular mechanisms controlling this maternal to zygotic transition, it is important to distinguish between the maternal and zygotic transcriptomes during this period. Here we use exome sequencing and RNA-seq to achieve such discrimination and in doing so have identified the first zygotic genes to be expressed in the embryo. Our work revealed different profiles of maternal mRNA post-transcriptional regulation prior to zygotic transcription initiation. Finally, we demonstrate that maternal mRNAs are required for different modes of zygotic transcription initiation, which is not simply dependent on the titration of factors that maintain genomic DNA in a transcriptionally incompetent state.

Gray platelet syndrome (GPS) is a predominantly recessive platelet disorder characterized by a mild thrombocytopenia with large platelets and a paucity of α-granules; these abnormalities cause mostly moderate but in rare cases severe bleeding. We sequenced the exomes of four unrelated cases and identified as the causative gene NBEAL2, a gene with previously unknown function but a member of a gene family involved in granule development. Silencing of nbeal2 in zebrafish abrogated thrombocyte formation.

Biliary atresia (BA) is a rapidly progressive and destructive fibrotic disorder of unknown etiology affecting the extrahepatic biliary tree of neonates. Epidemiological studies suggest that an environmental factor, such as a virus or toxin, is the cause of the disease, although none have been definitively established. Several naturally occurring outbreaks of BA in Australian livestock have been associated with the ingestion of unusual plants by pregnant animals during drought conditions. We used a biliary secretion assay in zebrafish to isolate a previously undescribed isoflavonoid, biliatresone, from Dysphania species implicated in a recent BA outbreak. This compound caused selective destruction of the extrahepatic, but not intrahepatic, biliary system of larval zebrafish. A mutation that enhanced biliatresone toxicity mapped to a region of the zebrafish genome that has conserved synteny with an established human BA susceptibility locus. The toxin also caused loss of cilia in neonatal mouse extrahepatic cholangiocytes in culture and disrupted cell polarity and monolayer integrity in cholangiocyte spheroids. Together, these findings provide direct evidence that BA could be initiated by perinatal exposure to an environmental toxin.

Laminins are major constituents of basement membranes and have wide ranging functions during development and in the adult. They are a family of heterotrimeric molecules created through association of an alpha, beta and gamma chain. We previously reported that two zebrafish loci, grumpy (gup) and sleepy (sly), encode laminin beta1 and gamma1, which are important both for notochord differentiation and for proper intersegmental blood vessel (ISV) formation. In this study we show that bashful (bal) encodes laminin alpha1 (lama1). Although the strongest allele, bal(m190), is fully penetrant, when compared to gup or sly mutant embryos, bal mutants are not as severely affected, as only anterior notochord fails to differentiate and ISVs are unaffected. This suggests that other alpha chains, and hence other isoforms, act redundantly to laminin 1 in posterior notochord and ISV development. We identified cDNA sequences for lama2, lama4 and lama5 and disrupted the expression of each alone or in mutant embryos also lacking laminin alpha1. When expression of laminin alpha4 and laminin alpha1 are simultaneously disrupted, notochord differentiation and ISVs are as severely affected as sly or gup mutants. Moreover, live imaging of transgenic embryos expressing enhanced green fluorescent protein in forming ISVs reveals that the vascular defects in these embryos are due to an inability of ISV sprouts to migrate correctly along the intersegmental, normally laminin-rich regions.

Antigen processing and presentation genes found within the MHC are among the most highly polymorphic genes of vertebrate genomes, providing populations with diverse immune responses to a wide array of pathogens. Here, we describe transcriptome, exome, and whole-genome sequencing of clonal zebrafish, uncovering the most extensive diversity within the antigen processing and presentation genes of any species yet examined. Our CG2 clonal zebrafish assembly provides genomic context within a remarkably divergent haplotype of the core MHC region on chromosome 19 for six expressed genes not found in the zebrafish reference genome: mhc1uga, proteasome-β 9b (psmb9b), psmb8f, and previously unknown genes psmb13b, tap2d, and tap2e. We identify ancient lineages for Psmb13 within a proteasome branch previously thought to be monomorphic and provide evidence of substantial lineage diversity within each of three major trifurcations of catalytic-type proteasome subunits in vertebrates: Psmb5/Psmb8/Psmb11, Psmb6/ Psmb9/Psmb12, and Psmb7/Psmb10/Psmb13. Strikingly, nearby tap2 and MHC class I genes also retain ancient sequence lineages, indicating that alternative lineages may have been preserved throughout the entire MHC pathway since early diversification of the adaptive immune system ∼500 Mya. Furthermore, polymorphisms within the three MHC pathway steps (antigen cleavage, transport, and presentation) are each predicted to alter peptide specificity. Lastly, comparative analysis shows that antigen processing gene diversity is far more extensive than previously realized (with ancient coelacanth psmb8 lineages, shark psmb13, and tap2t and psmb10 outside the teleost MHC), implying distinct immune functions and conserved roles in shaping MHC pathway evolution throughout vertebrates.comparative genomics | proteasome and TAP evolution | major histocompatibility | MHC class I pathway | CG2 clonal zebrafish