Ross Kettleborough

A systematic genome-wide analysis of zebrafish protein-coding gene function

Since the publication of the human reference genome, the identities of specific genes associated with human diseases are being discovered at a rapid rate. A central problem is that the biological activity of these genes is often unclear. Detailed investigations in model vertebrate organisms, typically mice, have been essential for understanding the activities of many orthologues of these disease-associated genes. Although gene-targeting approaches and phenotype analysis have led to a detailed understanding of nearly 6,000 protein-coding genes, this number falls considerably short of the more than 22,000 mouse protein-coding genes. Similarly, in zebrafish genetics, one-by-one gene studies using positional cloning, insertional mutagenesis, antisense morpholino oligonucleotides, targeted re-sequencing, and zinc finger and TAL endonucleases have made substantial contributions to our understanding of the biological activity of vertebrate genes, but again the number of genes studied falls well short of the more than 26,000 zebrafish protein-coding genes. Importantly, for both mice and zebrafish, none of these strategies are particularly suited to the rapid generation of knockouts in thousands of genes and the assessment of their biological activity. Here we describe an active project that aims to identify and phenotype the disruptive mutations in every zebrafish protein-coding gene, using a well-annotated zebrafish reference genome sequence, high-throughput sequencing and efficient chemical mutagenesis. So far we have identified potentially disruptive mutations in more than 38% of all known zebrafish 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 1,000 alleles. All mutant alleles and data are available to the community and our phenotyping scheme is adaptable to phenotypic analysis beyond embryogenesis.

Exome sequencing identifies NBEAL2 as the causative gene for gray platelet syndrome

Gray platelet syndrome (GPS) is a predominantly recessive platelet disorder that is characterized by 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 individuals and identified NBEAL2 as the causative gene; it has no previously known function but is a member of a gene family that is involved in granule development. Silencing of nbeal2 in zebrafish abrogated thrombocyte formation.

Identification of the zebrafish maternal and paternal transcriptomes

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.

High-throughput target-selected gene inactivation in zebrafish

There is an increasing requirement for efficient reverse genetics in the zebrafish, Here we describe a method that takes advantage of conventional mutagenized libraries (identical to ones used in forward screens) and re-sequencing to identify ENU-induced mutations in genes of interest. The efficiency of TILLING (Targeting Induced Local Legions IN Genomes) depends on the rate of mutagenesis in the library being screened, the amount of base pairs screened, and the ability to effectively identify and retrieve mutations on interest. Here we show that by improving the mutagenesis protocol, using in silico methods to predict codon changes for target selection, efficient PCR and re-sequencing, and accurate mutation detection we can vastly improve current TILLING protocols. Importantly it is also possible to use this method for screening for splice and mis-sense mutations, and with even a relatively small library, there is a high chance of identifying mutations across any given gene.

Identification of a plant isoflavonoid that causes biliary atresia

A plant toxin with a previously uncharacterized chemical structure causes biliary atresia in zebrafish and mammals. Fishing for a toxin Biliary atresia is a childhood disease of unknown etiology, which causes rapid destruction of the biliary tree. Lorent et al. used a naturally occurring outbreak of biliary atresia in Australian lambs to discover that one cause of this disease may be a plant toxin. During a period of drought, the pregnant ewes’ diet included unusual species of plants growing on land that is normally covered by water. The researchers tested extracts from these plants in a zebrafish model and identified a compound that caused biliary atresia–like changes in the fish as well as in mammalian cells. Although these plants are not eaten by humans, the findings suggest that plant toxins may contribute to the development of biliary atresia in people as well. 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.

A radiation hybrid transcript map of the mouse genome.

Expressed-sequence tag (EST) maps are an adjunct to sequence-based analytical methods of gene detection and localization for those species for which such data are available, and provide anchors for high-density homology and orthology mapping in species for which large-scale sequencing has yet to be done. Species for which radiation hybrid–based transcript maps have been established include human, rat, mouse, dog, cat and zebrafish. We have established a comprehensive first-generation–placement radiation hybrid map of the mouse consisting of 5,904 mapped markers (3,993 ESTs and 1,911 sequence-tagged sites (STSs)). The mapped ESTs, which often originate from small-EST clusters, are enriched for genes expressed during early mouse embryogenesis and are probably different from those localized in humans. We have confirmed by in situ hybridization that even singleton ESTs, which are usually not retained for mapping studies, may represent bona fide transcribed sequences. Our studies on mouse chromosomes 12 and 14 orthologous to human chromosome 14 show the power of our radiation hybrid map as a predictive tool for orthology mapping in humans.

Multi-allelic phenotyping--a systematic approach for the simultaneous analysis of multiple induced mutations.

The zebrafish mutation project (ZMP) aims to generate a loss of function allele for every protein-coding gene, but importantly to also characterise the phenotypes of these alleles during the first five days of development. Such a large-scale screen requires a systematic approach both to identifying phenotypes, and also to linking those phenotypes to specific mutations. This phenotyping pipeline simultaneously assesses the consequences of multiple alleles in a two-step process. First, mutations that do not produce a visible phenotype during the first five days of development are identified, while a second round of phenotyping focuses on detailed analysis of those alleles that are suspected to cause a phenotype. Allele-specific PCR single nucleotide polymorphism (SNP) assays are used to genotype F2 parents and individual F3 fry for mutations known to be present in the F1 founder. With this method specific phenotypes can be linked to induced mutations. In addition a method is described for cryopreserving sperm samples of mutagenised males and their subsequent use for in vitro fertilisation to generate F2 families for phenotyping. Ultimately this approach will lead to the functional annotation of the zebrafish genome, which will deepen our understanding of gene function in development and disease.

Alternative haplotypes of antigen processing genes in zebrafish diverged early in vertebrate evolution

Significance Antigen presentation genes are exceptionally polymorphic, enhancing immune defense. Polymorphism within additional components of the MHC pathway, particularly the antigen processing genes, may also shape immune responses. Using transcriptome, exome, and whole-genome sequencing to examine immune gene variation in zebrafish, we uncovered several antigen processing genes not found in the reference genome clustered within a deeply divergent haplotype of the core MHC locus. Our data provide evidence that these previously undescribed antigen processing genes retain ancient alternative sequence lineages, likely derived during the formation of the adaptive immune system, and represent the most divergent collection of antigen processing and presentation genes yet identified. These findings offer insights into the evolution of vertebrate adaptive immunity. 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.

Erratum: The deubiquitinase USP9X suppresses pancreatic ductal adenocarcinoma (Nature (2012) 486 (266-270) DOI: 10.1038/nature11114)

This corrects the article DOI: 10.1038/nature11114

Erratum: Corrigendum: The deubiquitinase USP9X suppresses pancreatic ductal adenocarcinoma

This corrects the article DOI: 10.1038/nature11114