Search

BACKGROUND: Although there is clear evidence of genetic contributions to susceptibility to nicotine addiction, it has proved difficult to identify causal alleles and pathways from studies in humans. Mutagenesis in model species generates strong phenotypes not present in wildtype populations and can be used to identify biological mechanisms underlying quantifiable behaviours. We tested the hypothesis that a forward genetic screen for nicotine preference in zebrafish can predict loci and biological mechanisms influencing human smoking behaviour. METHODS: A population-based forward genetic screen of ethylnitrosurea- mutagenized zebrafish was used to identify lines of fish showing altered nicotine preference. Immunohistochemical, behavioral and quantitative PCR analyses were used to characterize mutant larvae. Focussed SNP analysis of the homologous human locus in cohorts from the UK and a Finnish Twin study assessed the predictive validity of the zebrafish data for human smoking behavior. RESULTS: We show nicotine preference is heritable in fish as in humans and identify loss-of-function mutations in the zebrafish Slit3 gene as associated with increased nicotine preference. Slit3 mutant larval fish showed altered sensitivity of habituation to acoustic startle to dopaminergic antagonists and increased Drd2 and Drd3 mRNA expression. Dopaminergic neuronal pathfinding was unaffected. Analysis of the SLIT3 locus in two independent human cohorts identified 2 genetic markers that predict level of cigarette consumption and likelihood of cessation. CONCLUSION: These findings suggest a role for SLIT3 signaling in development of dopaminergic pathways affecting behaviours associated with nicotine dependence and confirm the translational relevance of the zebrafish model in exploring complex human behaviors.

Dystroglycan, an extracellular matrix receptor, has essential functions in various tissues. Loss of α‐dystroglycan‐laminin interaction due to defective glycosylation of α‐dystroglycan underlies a group of congenital muscular dystrophies often associated with brain malformations, referred to as dystroglycanopathies. The lack of isogenic human dystroglycanopathy cell models has limited our ability to test potential drugs in a human‐ and neural‐specific context. Here, we generated induced pluripotent stem cells (iPSCs) from a severe dystroglycanopathy patient with homozygous FKRP (fukutin‐related protein gene) mutation. We showed that CRISPR/Cas9‐mediated gene correction of FKRP restored glycosylation of α‐dystroglycan in iPSC‐derived cortical neurons, whereas targeted gene mutation of FKRP in wild‐type cells disrupted this glycosylation. In parallel, we screened 31,954 small molecule compounds using a mouse myoblast line for increased glycosylation of α‐dystroglycan. Using human FKRP‐iPSC‐derived neural cells for hit validation, we demonstrated that compound 4‐(4‐bromophenyl)‐6‐ethylsulfanyl‐2‐oxo‐3,4‐dihydro‐1H‐pyridine‐5‐carbonitrile (4BPPNit) significantly augmented glycosylation of α‐dystroglycan, in part through upregulation of LARGE1 glycosyltransferase gene expression. Together, isogenic human iPSC‐derived cells represent a valuable platform for facilitating dystroglycanopathy drug discovery and therapeutic development.

The Deciphering the Mechanisms of Developmental Disorders programme has analysed the morphological and molecular phenotypes of embryonic and perinatal lethal mouse mutant lines in order to investigate the causes of embryonic lethality. Here we show that individual whole-embryo RNA-seq of 73 mouse mutant lines (>1000 transcriptomes) identifies transcriptional events underlying embryonic lethality and associates previously uncharacterised genes with specific pathways and tissues. For example, our data suggest that Hmgxb3 is involved in DNA-damage repair and cell-cycle regulation. Further, we separate embryonic delay signatures from mutant line-specific transcriptional changes by developing a baseline mRNA expression catalogue of wild-type mice during early embryogenesis (4–36 somites). Analysis of transcription outside coding sequence identifies deregulation of repetitive elements in Morc2a mutants and a gene involved in gene-specific splicing. Collectively, this work provides a large scale resource to further our understanding of early embryonic developmental disorders., The transcriptional signature of embryonic lethality has not been defined. Here, the authors, as part of the Deciphering the Mechanisms of Developmental Disorders programme, define genes causing murine embryonic lethality around E9.5 and identify developmental delay transcriptional signatures.

Skin pigment patterns are important, being under strong selection for multiple roles including camouflage and UV protection. Pigment cells underlying these patterns form from adult pigment stem cells (APSCs). In zebrafish, APSCs derive from embryonic neural crest cells, but sit dormant until activated to produce pigment cells during metamorphosis. The APSCs are set-aside in an ErbB signaling dependent manner, but the mechanism maintaining quiescence until metamorphosis remains unknown. Mutants for a pigment pattern gene, parade, exhibit ectopic pigment cells localised to the ventral trunk, but also supernumerary cells restricted to the Ventral Stripe. Contrary to expectations, these melanocytes and iridophores are discrete cells, but closely apposed. We show that parade encodes Endothelin receptor Aa, expressed in the blood vessels, most prominently in the medial blood vessels, consistent with the ventral trunk phenotype. We provide evidence that neuronal fates are not affected in parade mutants, arguing against transdifferentiation of sympathetic neurons to pigment cells. We show that inhibition of BMP signaling prevents specification of sympathetic neurons, indicating conservation of this molecular mechanism with chick and mouse. However, inhibition of sympathetic neuron differentiation does not enhance the parade phenotype. Instead, we pinpoint ventral trunk-restricted proliferation of neural crest cells as an early feature of the parade phenotype. Importantly, using a chemical genetic screen for rescue of the ectopic pigment cell phenotype of parade mutants (whilst leaving the embryonic pattern untouched), we identify ErbB inhibitors as a key hit. The time-window of sensitivity to these inhibitors mirrors precisely the window defined previously as crucial for the setting aside of APSCs in the embryo, strongly implicating adult pigment stem cells as the source of the ectopic pigment cells. We propose that a novel population of APSCs exists in association with medial blood vessels, and that their quiescence is dependent upon Endothelin-dependent factors expressed by the blood vessels., Author summary Pigment patterns are crucial for the many aspects of animal biology, for example, providing camouflage, enabling mate selection and protecting against UV irradiation. These patterns are generated by one or more pigment cell-types, localised in the skin, but derived from specialised stem cells (adult pigment stem cells, APSCs). In mammals, such as humans, but also in birds and fish, these APSCs derive from a transient population of multipotent progenitor cells, the neural crest. Formation of the adult pigment pattern is perhaps best studied in the zebrafish, where the adult pigment pattern is formed during a metamorphosis beginning around 21 days of development. The APSCs are set-aside in the embryo around 1 day of development, but then remain inactive until that metamorphosis, when they become activated to produce the adult pigment cells. We know something of how the cells are set-aside, but what signals maintain them in an inactive state is a mystery. Here we study a zebrafish mutant, called parade, which shows ectopic pigment cells in the embryo. We clone the parade gene, identifying it as ednraa encoding a component of a cell-cell communication process, which is expressed in blood vessels. By characterising the changes in the neural crest and in the pigment cells formed, and by combining this with an innovative assay identifying drugs that prevent the ectopic cells from forming, we deduce that the ectopic cells in the larva derive from precocious activation of APSCs to form pigment cells. We propose that a novel population of APSCs are associated with the blood vessels, that these are held in a quiescent state by signals coming from these vessels, and that these signals depend upon ednraa. Together this opens up an exciting opportunity to identify the signals maintaining APSC quiescence in zebrafish.

The neural crest (NC) is a vertebrate-specific cell type that contributes to a wide range of different tissues across all three germ layers. The gene regulatory network (GRN) responsible for the formation of neural crest is conserved across vertebrates. Central to the induction of the NC GRN are AP-2 and SoxE transcription factors. NC induction robustness is ensured through the ability of some of these transcription factors to compensate loss of function of gene family members. However the gene regulatory events underlying compensation are poorly understood. We have used gene knockout and RNA sequencing strategies to dissect NC induction and compensation in zebrafish. We genetically ablate the NC using double mutants of tfap2a;tfap2c or remove specific subsets of the NC with sox10 and mitfa knockouts and characterise genome-wide gene expression levels across multiple time points. We find that compensation through a single wild-type allele of tfap2c is capable of maintaining early NC induction and differentiation in the absence of tfap2a function, but many target genes have abnormal expression levels and therefore show sensitivity to the reduced tfap2 dosage. This separation of morphological and molecular phenotypes identifies a core set of genes required for early NC development. We also identify the 15 somites stage as the peak of the molecular phenotype which strongly diminishes at 24 hpf even as the morphological phenotype becomes more apparent. Using gene knockouts, we associate previously uncharacterised genes with pigment cell development and establish a role for maternal Hippo signalling in melanocyte differentiation. This work extends and refines the NC GRN while also uncovering the transcriptional basis of genetic compensation via paralogues., Author summary Embryonic development is an intricate process that requires genes to be active at the right time and place. Organisms have evolved mechanisms that ensure faithful execution of developmental programmes even if genes fail to function. For example, in a process called genetic compensation, one or more genes become activated in response to loss of function of another. In this work we use the zebrafish model to investigate how two related genes, tfap2a and tfap2c, interact to ensure establishment of the neural crest, a vertebrate-specific cell type that contributes to many different tissues. Losing tfap2a activity causes mild morphological defects and losing tfap2c has no visible effect. Yet when both are inactive, embryos are severely abnormal due to lack of neural crest-derived tissues. Here we show that loss of tfap2a triggers upregulation of tfap2c which prevents the loss of neural crest tissue. However, the genes normally regulated by tfap2a respond differently to tfap2c allowing us to identify the first tier of the Ap2 network and new players in neural crest biology. Our work demonstrates that the expression signature of partial, but morphologically sufficient, genetic compensation provides an opportunity to dissect gene regulatory networks.

We have produced an mRNA expression time course of zebrafish development across 18 time points from 1 cell to 5 days post-fertilisation sampling individual and pools of embryos. Using poly(A) pulldown stranded RNA-seq and a 3′ end transcript counting method we characterise temporal expression profiles of 23,642 genes. We identify temporal and functional transcript co-variance that associates 5024 unnamed genes with distinct developmental time points. Specifically, a class of over 100 previously uncharacterised zinc finger domain containing genes, located on the long arm of chromosome 4, is expressed in a sharp peak during zygotic genome activation. In addition, the data reveal new genes and transcripts, differential use of exons and previously unidentified 3′ ends across development, new primary microRNAs and temporal divergence of gene paralogues generated in the teleost genome duplication. To make this dataset a useful baseline reference, the data can be browsed and downloaded at Expression Atlas and Ensembl.

Background: Mutations in proteins involved in telomere maintenance lead to a range of human diseases, including dyskeratosis congenita, idiopathic pulmonary fibrosis and cancer. Telomerase functions to add telomeric repeats back onto the ends of chromosomes, however non-canonical roles of components of telomerase have recently been suggested. Methods: Here we use a zebrafish telomerase mutant which harbours a nonsense mutation in tert to investigate the adult phenotypes of fish derived from heterozygous parents of different ages. Furthermore we use whole genome sequencing data to estimate average telomere lengths. Results: We show that homozygous offspring from older heterozygotes exhibit signs of body wasting at a younger age than those of younger parents, and that offspring of older heterozygous parents weigh less irrespective of genotype. We also demonstrate that tert homozygous mutant fish have a male sex bias, and that clutches from older parents also have a male sex bias in the heterozygous and wild-type populations. Telomere length analysis reveals that the telomeres of younger heterozygous parents are shorter than those of older heterozygous parents. Conclusions: These data indicate that the phenotypes observed in offspring from older parents cannot be explained by telomere length. Instead we propose that Tert functions outside of telomere length maintenance in an age-dependent manner to influence the adult phenotypes of the next generation.

Large-scale forward genetic screens have been instrumental for identifying genes that regulate development, homeostasis, and regeneration, as well as the mechanisms of disease. The zebrafish, Danio rerio, is an established genetic and developmental model used in genetic screens to uncover genes necessary for early development. However, the regulation of postembryonic development has received less attention as these screens are more labor intensive and require extensive resources. The lack of systematic interrogation of late development leaves large aspects of the genetic regulation of adult form and physiology unresolved. To understand the genetic control of postembryonic development, we performed a dominant screen for phenotypes affecting the adult zebrafish. In our screen, we identified 72 adult viable mutants showing changes in the shape of the skeleton as well as defects in pigmentation. For efficient mapping of these mutants and mutation identification, we devised a new mapping strategy based on identification of mutant-specific haplotypes. Using this method in combination with a candidate gene approach, we were able to identify linked mutations for 22 out of 25 mutants analyzed. Broadly, our mutational analysis suggests that there are key genes and pathways associated with late development. Many of these pathways are shared with humans and are affected in various disease conditions, suggesting constraint in the genetic pathways that can lead to change in adult form. Taken together, these results show that dominant screens are a feasible and productive means to identify mutations that can further our understanding of gene function during postembryonic development and in disease.

We have produced an mRNA expression time course of zebrafish development across 18 time points from 1-cell to 5 days post-fertilisation sampling individual and pools of embryos. Using poly(A) pulldown stranded RNA-seq and a 3′ end transcript counting method we characterise the temporal expression profiles of 23,642 genes. We identify temporal and functional transcript co-variance that associates 5,024 unnamed genes with distinct developmental time points. Specifically, a class of over 100 previously uncharacterised zinc finger domain containing genes, located on the long arm of chromosome 4, is expressed in a sharp peak during zygotic genome activation. The data reveal complex and widespread differential use of exons and previously unidentified 3′ ends across development, new primary microRNA transcripts and temporal divergence of gene paralogues generated in the teleost genome duplication. To make this dataset a useful baseline reference, the data are accessible to browse and download at Expression Atlas and Ensembl.

SummaryMutations in genes required for functional glycosylation of α-dystroglycan cause a group of congenital muscular dystrophies associated with brain malformations, referred to as dystroglycanopathies. The lack of isogenic, physiology-relevant human cellular models has limited our understanding of the cortical abnormalities in dystroglycanopathies. Here we generate induced pluripotent stem cells (iPSCs) from a severe dystroglycanopathy patient with homozygous mutations in the ribitol-5-phosphate transferase gene, FKRP. We carry out targeted gene correction in FKRP-iPSCs using CRISPR/Cas9-mediated genome editing. We characterise the directed differentiation of FKRP- and corrected-iPSCs to neural stem cells, cortical progenitors and cortical neurons. Importantly, we show that targeted gene correction of FKRP restores functional glycosylation of α-dystroglycan in iPSC-derived cortical neurons. We independently validate this result by showing targeted gene mutation of FKRP disrupts functional glycosylation of α-dystroglycan. This work demonstrates the feasibility of using CRISPR/Cas9-engineered human iPSCs for modelling dystroglycanopathies and provides a foundation for therapeutic development.HighlightsGeneration of FKRP-iPSCs for modelling cortical abnormalities in dystroglycanopathiesPrecise gene correction by CRISPR/Cas9-mediated genome editingDirected differentiation of isogenic control and FKRP-iPSC to cortical neuronsFunctional glycosylation of α-dystroglycan is restored in cortical neurons derived from CRISPR/Cas9-corrected iPSCsTargeted gene mutation of FKRP disrupts functional glycosylation of α-dystroglycan in cortical neurons

KDM2A is a histone demethylase associated with transcriptional silencing, however very little is known about its in vivo role in development and disease. Here we demonstrate that loss of the orthologue kdm2aa in zebrafish causes widespread transcriptional disruption and leads to spontaneous melanomas at a high frequency. Fish homozygous for two independent premature stop codon alleles show reduced growth and survival, a strong male sex bias, and homozygous females exhibit a progressive oogenesis defect. kdm2aa mutant fish also develop melanomas from early adulthood onwards which are independent from mutations in braf and other common oncogenes and tumour suppressors as revealed by deep whole exome sequencing. In addition to effects on translation and DNA replication gene expression, high-replicate RNA-seq in morphologically normal individuals demonstrates a stable regulatory response of epigenetic modifiers and the specific de-repression of a group of zinc finger genes residing in constitutive heterochromatin. Together our data reveal a complex role for Kdm2aa in regulating normal mRNA levels and carcinogenesis. These findings establish kdm2aa mutants as the first single gene knockout model of melanoma biology., Author summary Epigenetic modifications of DNA and histones, the major components of chromatin, play a central role in transcriptional regulation. KDM2A is a histone demethylase that integrates DNA and histone modification signatures and is involved in transcriptional silencing through heterochromatin maintenance. Here we show that adult zebrafish homozygous for the orthologue kdm2aa develop melanomas, a malignant form of skin cancer, independently from oncogenes known to drive melanoma formation. We observe that transcript abundance is widely affected in kdm2aa mutants and find that gene expression of several DNA- and histone- modifying enzymes is stably altered. We furthermore demonstrate a specific de-repression of a group of genes encoding zinc finger-containing proteins that has the potential to be involved in transcriptional regulation. We suggest that these molecular disruptions underlie the melanoma formation, as well as the other observed phenotypes such as reduced growth and survival, a male sex bias and an oogenesis defect. This work demonstrates in vivo a role for Kdm2aa as a tumour suppressor and establishes, to our knowledge, kdm2aa-deficient fish as the first single gene knockout vertebrate model of melanoma.

SummaryThe retinal anterior homeobox (rax) gene encodes a transcription factor necessary for vertebrate eye development. rax transcription is initiated at the end of gastrulation in Xenopus, and is a key part of the regulatory network specifying anterior neural plate and retina. We describe here a Xenopus tropicalis rax mutant, the first mutant analyzed in detail from a reverse genetic screen. As in other vertebrates, this nonsense mutation results in eyeless animals, and is lethal peri-metamorphosis. Tissue normally fated to form retina in these mutants instead forms tissue with characteristics of diencephalon and telencephalon. This implies that a key role of rax, in addition to defining the eye field, is in preventing alternative forebrain identities. Our data highlight that brain and retina regions are not determined by the mid-gastrula stage but are by the neural plate stage. An RNA-Seq analysis and in situ hybridization assays for early gene expression in the mutant revealed that several key eye field transcription factors (e.g. pax6, lhx2 and six6) are not dependent on rax activity through neurulation. However, these analyses identified other genes either up- or down-regulated in mutant presumptive retinal tissue. Two neural patterning genes of particular interest that appear up-regulated in the rax mutant RNA-seq analysis are hesx1 and fezf2. These genes were not previously known to be regulated by rax. The normal function of rax is to partially repress their expression by an indirect mechanism in the presumptive retina region in wildtype embryos, thus accounting for the apparent up-regulation in the rax mutant. Knock-down experiments using antisense morpholino oligonucleotides directed against hesx1 and fezf2 show that failure to repress these two genes contributes to transformation of presumptive retinal tissue into non-retinal forebrain identities in the rax mutant.

The initial phases of embryonic development occur in the absence of de novo transcription and are instead controlled by maternally inherited mRNAs and proteins. During this initial period, cell cycles are synchronous and lack gap phases. Following this period of transcriptional silence, zygotic transcription begins, the maternal influence on development starts to decrease, and dramatic changes to the cell cycle take place. Here, we discuss recent work that is shedding light on the maternal to zygotic transition and the interrelated but distinct mechanisms regulating the onset of zygotic transcription and changes to the cell cycle during early embryonic development.

The need for scalable strategies to probe the biological consequences of candidate cancer genes has never been more pressing. The zebrafish, with its capacity for high-throughput transgenesis, in vivo imaging and chemical/genetic screening, has ideal features for undertaking this task. Unique biological insights from zebrafish have already led to the identification of novel oncogenic drivers and small molecules being used to treat the human cancer. This review summarizes the recent main findings and describes pertinent areas where the zebrafish can greatly contribute to our understanding of cancer biology and treatment.

The zebrafish research community is celebrating! The zebrafish genome has recently been sequenced, the Zebrafish Mutation Project (launched by the Wellcome Trust Sanger Institute) has published the results of its first large-scale ethylnitrosourea (ENU) mutagenesis screen, and a host of new techniques, such as the genome editing technologies TALEN and CRISPR-Cas, are enabling specific mutations to be created in model organisms and investigated in vivo. The zebrafish truly seems to be coming of age. These powerful resources invoke the question of whether zebrafish can be increasingly used to model human disease, particularly common, chronic diseases of metabolism such as obesity and type 2 diabetes. In recent years, there has been considerable success, mainly from genomic approaches, in identifying genetic variants that are associated with these conditions in humans; however, mechanistic insights into the role of implicated disease loci are lacking. In this Review, we highlight some of the advantages and disadvantages of zebrafish to address the organism’s utility as a model system for human metabolic diseases.

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.

The blood group Vel was discovered 60 years ago(1), but the underlying gene is unknown. Individuals negative for the Vel antigen are rare and are required for the safe transfusion of patients with antibodies to Vel(2). To identify the responsible gene, we sequenced the exomes of five individuals negative for the Vel antigen and found that four were homozygous and one was heterozygous for a low-frequency 17-nucleotide frameshift deletion in the gene encoding the 78-amino-acid transmembrane protein SMIM1. A follow-up study showing that 59 of 64 Vel-negative individuals were homozygous for the same deletion and expression of the Vel antigen on SMIM1-transfected cells confirm SMIM1 as the gene underlying the Vel blood group. An expression quantitative trait locus (eQTL), the common SNP rs1175550 contributes to variable expression of the Vel antigen (P = 0.003) and influences the mean hemoglobin concentration of red blood cells (RBCs; P = 8.6 x 10(-15))(3). In vivo, zebrafish with smim1 knockdown showed a mild reduction in the number of RBCs, identifying SMIM1 as a new regulator of RBC formation. Our findings are of immediate relevance, as the homozygous presence of the deletion allows the unequivocal identification of Vel-negative blood donors.

Autosomal-recessive albinism is a hypopigmentation disorder with a broad phenotypic range. A substantial fraction of individuals with albinism remain genetically unresolved, and it has been hypothesized that more genes are to be identified. By using homozygosity mapping of an inbred Faroese family, we identified a 3.5 Mb homozygous region (10q22.2–q22.3) on chromosome 10. The region contains five protein-coding genes, and sequencing of one of these, C10orf11, revealed a nonsense mutation that segregated with the disease and showed a recessive inheritance pattern. Investigation of additional albinism-affected individuals from the Faroe Islands revealed that five out of eight unrelated affected persons had the nonsense mutation in C10orf11. Screening of a cohort of autosomal-recessive-albinism-affected individuals residing in Denmark showed a homozygous 1 bp duplication in C10orf11 in an individual originating from Lithuania. Immunohistochemistry showed localization of C10orf11 in melanoblasts and melanocytes in human fetal tissue, but no localization was seen in retinal pigment epithelial cells. Knockdown of the zebrafish (Danio rerio) homolog with the use of morpholinos resulted in substantially decreased pigmentation and a reduction of the apparent number of pigmented melanocytes. The morphant phenotype was rescued by wild-type C10orf11, but not by mutant C10orf11. In conclusion, we have identified a melanocyte-differentiation gene, C10orf11, which when mutated causes autosomal-recessive albinism in humans.

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.

Multiple glycosyltransferases are essential for the proper modification of alpha-dystroglycan, as mutations in the encoding genes cause congenital/limb-girdle muscular dystrophies. Here we elucidate further the structure of an O-mannose-initiated glycan on alpha-dystroglycan that is required to generate its extracellular matrix-binding polysaccharide. This functional glycan contains a novel ribitol structure that links a phosphotrisaccharide to xylose. ISPD is a CDP-ribitol (ribose) pyrophosphorylase that generates the reduced sugar nucleotide for the insertion of ribitol in a phosphodiester linkage to the glycoprotein. TMEM5 is a UDP-xylosyl transferase that elaborates the structure. We demonstrate in a zebrafish model as well as in a human patient that defects in TMEM5 result in muscular dystrophy in combination with abnormal brain development. Thus, we propose a novel structure—a ribitol in a phosphodiester linkage—for the moiety on which TMEM5, B4GAT1, and LARGE act to generate the functional receptor for ECM proteins having LG domains. DOI: http://dx.doi.org/10.7554/eLife.14473.001

Background The CRISPR/Cas9 system is a prokaryotic immune system that infers resistance to foreign genetic material and is a sort of 'adaptive immunity'. It has been adapted to enable high throughput genome editing and has revolutionised the generation of targeted mutations. Results We have developed a scalable analysis pipeline to identify CRISPR/Cas9 induced mutations in hundreds of samples using next generation sequencing (NGS) of amplicons. We have used this system to investigate the best way to screen mosaic Zebrafish founder individuals for germline transmission of induced mutations. Screening sperm samples from potential founders provides much better information on germline transmission rates and crucially the sequence of the particular insertions/deletions (indels) that will be transmitted. This enables us to combine screening with archiving to create a library of cryopreserved samples carrying known mutations. It also allows us to design efficient genotyping assays, making identifying F1 carriers straightforward. Conclusions The methods described will streamline the production of large numbers of knockout alleles in selected genes for phenotypic analysis, complementing existing efforts using random mutagenesis. Electronic supplementary material The online version of this article (doi:10.1186/s12864-016-2563-z) contains supplementary material, which is available to authorized users.

Chordoma is a rare malignant tumour of bone, the molecular marker of which is the expression of the transcription factor, brachyury. Having recently demonstrated that silencing brachyury induces growth arrest in a chordoma cell line, we now seek to identify its downstream target genes. Here we use an integrated functional genomics approach involving shRNA-mediated brachyury knockdown, gene expression microarray, ChIP-seq experiments, and bioinformatics analysis to achieve this goal. We confirm that the T-box binding motif of human brachyury is identical to that found in mouse, Xenopus, and zebrafish development, and that brachyury acts primarily as an activator of transcription. Using human chordoma samples for validation purposes, we show that brachyury binds 99 direct targets and indirectly influences the expression of 64 other genes, thereby acting as a master regulator of an elaborate oncogenic transcriptional network encompassing diverse signalling pathways including components of the cell cycle, and extracellular matrix components. Given the wide repertoire of its active binding and the relative specific localization of brachyury to the tumour cells, we propose that an RNA interference-based gene therapy approach is a plausible therapeutic avenue worthy of investigation.

As vertebrate transcriptomes continue to be scrutinized they reveal ever-increasing levels of complexity. Deciphering the transcribed regions from the nontranscribed regions and presenting a comprehensive, yet artifact free, gene set is a significant undertaking. Annotating the gene content of a genomic reference sequence is fundamental to understanding the biological processes of the organism. Moreover, for a model organism like zebrafish, being able to link a human gene of medical interest to the zebrafish ortholog is key to elucidating human gene function. Gene annotation methods rely on cDNA sequencing or ab initio prediction, or a combination of both. Vertebrate methods tend to rely on sequence information. The Vertebrate Genome Annotation Database (VEGA) (Wilming et al. 2008) provides a manually curated set of gene models derived from transcript sequence data, while the RefSeq Database (Pruitt et al. 2009) provides a predicted and manually curated gene set based on genomic and transcript sequence data. In contrast, Ensembl provides an automated annotation system for vertebrate and other eukaryote species (Flicek et al. 2011). An Ensembl genebuild predicts gene models using a publicly available species-specific and orthologous sequence aligned to the reference genome. Some organisms have a lot of species-specific cDNA sequence data, while others have little. As more information enters the public databases, particularly species-specific data, the Ensembl annotation improves in both quantity and quality. The advent of high-throughput transcriptome sequencing, RNA-seq (Wang et al. 2009), marks a major advance in species-specific experimentally derived sequence data and has the potential to make a significant impact on Ensembl genebuilds. In particular, RNA-seq can add untranslated regions of predicted orthologous transcript and provide proof of transcription. Recently, new sequencing technologies, such as the Illumina Genome Analyzer (Bentley et al. 2008) used in this study, have increased sequencing capacity and enabled alternative strategies for many sequence-based investigations. RNA-seq, for example, has allowed extremely deep sequencing of complementary DNA to an extent that was not possible using cDNA libraries and capillary sequencing. Unlike traditional directed cDNA-sequencing strategies (Temple et al. 2009), multitranscript sampling and the depth of the sequence in RNA-seq reduce noise caused by occasional misspliced mRNA. RNA-seq also assays the frequency of alternative splice forms and their spatial and temporal expression patterns, giving a comprehensive snapshot of the transcription of the whole sample. Studies have used RNA-seq to look at the transcriptome of yeast (Nagalakshmi et al. 2008), fission yeast (Wilhelm et al. 2008), mouse (Mortazavi et al. 2008), and human (Cloonan et al. 2008; Sultan et al. 2008; Wang et al. 2008) and found new levels of complexity. Methods have been published that use RNA-seq sequence reads alone to build gene models in a reference genome guided manner (Denoeud et al. 2008; Yassour et al. 2009; Guttman et al. 2010; Trapnell et al. 2010) or assemble the transcriptome independently of a reference sequence (Zerbino and Birney 2008; Robertson et al. 2010), in one case using a related proteome to aid assembly (Surget-Groba and Montoya-Burgos 2010). RNA-seq data is just one piece of information that can be used to predict a gene, but it is a cost-effective way to increase the amount of species-specific cDNA data for the less well-studied genomes. Recently, zebrafish gene collections derived using high-throughput sequencing technologies have been published, which look at the transcriptome of early development and lincRNA (Ulitsky et al. 2011; Pauli et al. 2012). Our aim was to distill a large quantity of zebrafish RNA-seq short reads into a useful gene annotation and incorporate this into the Ensembl genebuild pipeline. RNA-seq assays all polyA transcriptional activity in the sample, and we present our initial endeavor at a concise interpretation. In an attempt to minimize the loss of biologically significant data while excluding artifacts, we have incorporated numerous filtering steps throughout the genebuilding process. Our RNA-seq-only genebuild and Ensembl/VEGA genebuild bring complementary elements to the final gene set. Incorporating the RNA-seq data into the Ensembl/VEGA genebuilds adds a large amount of new species-specific cDNA sequence to current publicly available gene data providing a comprehensive gene annotation. In this analysis we have concentrated on protein-coding genes, although there is scope for further modifications to annotate all of the transcribed sequence. We present an advance in the annotation of the zebrafish genome and a template for improving Ensembl annotation in other species.

Graphical abstract Highlights ► Automatized characterization of thrombus formation in time-lapse microscopy. ► Novel energy model for segmentation of multiple dynamic regions. ► Novel algorithm for the joint segmentation of thrombus and aortic regions. ► Exhaustive validation on synthetic and real microscopic data., The characterization of thrombus formation in time-lapse DIC microscopy is of increased interest for identifying genes which account for atherothrombosis and coronary artery diseases (CADs). In particular, we are interested in large-scale studies on zebrafish, which result in large amount of data, and require automatic processing. In this work, we present an image-based solution for the automatized extraction of parameters quantifying the temporal development of thrombotic plugs. Our system is based on the joint segmentation of thrombotic and aortic regions over time. This task is made difficult by the low contrast and the high dynamic conditions observed in vivo DIC microscopic scenes. Our key idea is to perform this segmentation by distinguishing the different motion patterns in image time series rather than by solving standard image segmentation tasks in each image frame. Thus, we are able to compensate for the poor imaging conditions. We model motion patterns by energies based on the idea of dynamic textures, and regularize the model by two prior energies on the shape of the aortic region and on the topological relationship between the thrombus and the aorta. We demonstrate the performance of our segmentation algorithm by qualitative and quantitative experiments on synthetic examples as well as on real in vivo microscopic sequences.

Contains fulltext : 108772.pdf (Publisher’s version ) (Open Access) Walker-Warburg syndrome (WWS) is an autosomal recessive multisystem disorder characterized by complex eye and brain abnormalities with congenital muscular dystrophy (CMD) and aberrant a-dystroglycan glycosylation. Here we report mutations in the ISPD gene (encoding isoprenoid synthase domain containing) as the second most common cause of WWS. Bacterial IspD is a nucleotidyl transferase belonging to a large glycosyltransferase family, but the role of the orthologous protein in chordates is obscure to date, as this phylum does not have the corresponding non-mevalonate isoprenoid biosynthesis pathway. Knockdown of ispd in zebrafish recapitulates the human WWS phenotype with hydrocephalus, reduced eye size, muscle degeneration and hypoglycosylated a-dystroglycan. These results implicate ISPD in a-dystroglycan glycosylation in maintaining sarcolemma integrity in vertebrates.

Model organisms have played a huge part in the history of studies of human genetic disease, both in identifying disease genes and characterizing their normal and abnormal functions. But is the importance of model organisms diminishing? The direct discovery of disease genes and variants in humans has been revolutionized, first by genome-wide association studies and now by whole-genome sequencing. Not only is it now much easier to directly identify potential disease genes in humans, but the genetic architecture that is being revealed in many cases is hard to replicate in model organisms. Furthermore, disease modelling can be done with increasing effectiveness using human cells. Where does this leave non-human models of disease?

Allelic mutations in putative glycosyltransferase genes, fukutin and fukutin-related protein (fkrp), lead to a wide range of muscular dystrophies associated with hypoglycosylation of a-dystroglycan, commonly referred to as dystroglycanopathies. Defective glycosylation affecting dystroglycan‐ligand interactions is considered to underlie the disease pathogenesis. We have modelled dystroglycanopathies in zebrafish using a novel loss-of-function dystroglycan allele and by inhibition of Fukutin family protein activities. We show that muscle pathology in embryos lacking Fukutin or FKRP is different from loss of dystroglycan. In addition to hypoglycosylated a-dystroglycan, knockdown of Fukutin or FKRP leads to a notochord defect and a perturbation of laminin expression before muscle degeneration. These are a consequence of endoplasmic reticulum stress and activation of the unfolded protein response (UPR), preceding loss of dystroglycan‐ ligand interactions. Together, our results suggest that Fukutin family proteins may play important roles in protein secretion and that the UPR may contribute to the phenotypic spectrum of some dystroglycanopathies in humans.

Background: Mutations in proteins involved in telomere maintenance lead to a range of human diseases, including dyskeratosis congenita, idiopathic pulmonary fibrosis and cancer. Telomerase functions to add telomeric repeats back onto the ends of chromosomes, however non-canonical roles of components of telomerase have recently been suggested.Methods: Here we use a zebrafish telomerase mutant which harbours a nonsense mutation intertto investigate the adult phenotypes of fish derived from heterozygous parents of different ages. Furthermore we use whole genome sequencing data to estimate average telomere lengths.Results: We show that homozygous offspring from older heterozygotes exhibit signs of body wasting at a younger age than those of younger parents, and that offspring of older heterozygous parents weigh less irrespective of genotype. We also demonstrate thatterthomozygous mutant fish have a male sex bias, and that clutches from older parents also have a male sex bias in the heterozygous and wild-type populations. Telomere length analysis reveals that the telomeres of younger heterozygous parents are shorter than those of older heterozygous parents.Conclusions: These data indicate that the phenotypes observed in offspring from older parents cannot be explained by telomere length. Instead we propose that Tert functions outside of telomere length maintenance in an age-dependent manner to influence the adult phenotypes of the next generation.

In humans, OFD1 is mutated in oral-facial-digital type I syndrome leading to prenatal death in hemizygous males and dysmorphic faces and brain malformations, with polycystic kidneys presenting later in life in heterozygous females. To elucidate the function of Ofd1, we have studied its function during zebrafish embryonic development. In wild-type embryos, ofd1 mRNA is widely expressed and Ofd1-green fluorescent protein (GFP) fusion localizes to the centrosome/basal body. Disrupting Ofd1 using antisense morpholinos (MOs) led to bent body axes, hydrocephalus and oedema. Laterality was randomized in the brain, heart and viscera, likely a consequence of shorter cilia with disrupted axonemes and perturbed intravesicular fluid flow in Kupffer's vesicle. Embryos injected with ofd1 MOs also displayed convergent extension (CE) defects, which were enhanced by loss of Slb/Wnt11 or Tri/Vangl2, two proteins functioning in a non-canonical Wnt/Planar Cell Polarity (PCP) pathway. Pronephric glomerular midline fusion was compromised in vangl2 and ofd1 loss of function embryos and we suggest this anomaly may be a novel CE defect. Thus, Ofd1 is required for ciliary motility and function in zebrafish, supporting data showing that Ofd1 is essential for primary cilia function in mice. In addition, our data show that Ofd1 is important for CE during gastrulation, consistent with data linking primary cilia and non-canonical Wnt/PCP signalling. © 2008 The Author(s).

Ribosome biogenesis is a ubiquitous and essential process in cells. Defects in ribosome biogenesis and function result in a group of human disorders, collectively known as ribosomopathies. In this study, we describe a zebrafish mutant with a loss-of-function mutation in nol9, a gene that encodes a non-ribosomal protein involved in rRNA processing. nol9 sa1022/sa1022 mutants have a defect in 28S rRNA processing. The nol9 sa1022/sa1022 larvae display hypoplastic pancreas, liver and intestine and have decreased numbers of hematopoietic stem and progenitor cells (HSPCs), as well as definitive erythrocytes and lymphocytes. In addition, ultrastructural analysis revealed signs of pathological processes occurring in endothelial cells of the caudal vein, emphasizing the complexity of the phenotype observed in nol9 sa1022/sa1022 larvae. We further show that both the pancreatic and hematopoietic deficiencies in nol9 sa1022/sa1022 embryos were due to impaired cell proliferation of respective progenitor cells. Interestingly, genetic loss of Tp53 rescued the HSPCs but not the pancreatic defects. In contrast, activation of mRNA translation via the mTOR pathway by L-Leucine treatment did not revert the erythroid or pancreatic defects. Together, we present the nol9 sa1022/sa1022 mutant, a novel zebrafish ribosomopathy model, which recapitulates key human disease characteristics. The use of this genetically tractable model will enhance our understanding of the tissue-specific mechanisms following impaired ribosome biogenesis in the context of an intact vertebrate., Author Summary The production of ribosomes, the protein-synthesizing machines, is fundamental in all cells. It is a very complex process that requires the coordinated actions of ribosomal and non-ribosomal proteins. Impairment of ribosome formation and function leads to a class of disorders known as “ribosomopathies”. Here, we describe the identification and characterization of a zebrafish mutant in nol9, a gene encoding a non-ribosomal protein involved in ribosome biogenesis. The nol9 sa1022/sa1022 mutants show defects in the exocrine pancreas and erythrocytes due to impaired cell proliferation. These phenotypic features of nol9 sa1022/sa1022 mutants are reminiscent of the clinical symptoms of Shwachman-Diamond syndrome, a ribosomopathy characterized by exocrine pancreatic insufficiency and hematopoietic defects. Interestingly, we found that hematopoiesis but not pancreas morphogenesis in nol9 sa1022/sa1022 larvae is Tp53-dependent, highlighting that the consequences of impaired ribosome biogenesis differ between tissues within the same organism. This study provides novel insight into the function of the ribosome biogenesis protein Nol9 in zebrafish development and presents a novel model that will help to decipher the tissue-specific mechanisms of ribosomopathies.

Background We present a genome-wide messenger RNA (mRNA) sequencing technique that converts small amounts of RNA from many samples into molecular phenotypes. It encompasses all steps from sample preparation to sequence analysis and is applicable to baseline profiling or perturbation measurements. Results Multiplex sequencing of transcript 3′ ends identifies differential transcript abundance independent of gene annotation. We show that increasing biological replicate number while maintaining the total amount of sequencing identifies more differentially abundant transcripts. Conclusions This method can be implemented on polyadenylated RNA from any organism with an annotated reference genome and in any laboratory with access to Illumina sequencing. Electronic supplementary material The online version of this article (doi:10.1186/s12864-015-1788-6) contains supplementary material, which is available to authorized users.

The RAB5 gene family is the best characterised of all human RAB families and is essential for in vitro homotypic fusion of early endosomes. In recent years, the disruption or activation of Rab5 family proteins has been used as a tool to understand growth factor signal transduction in whole animal systems such as Drosophila melanogaster and zebrafish. In this study we have examined the functions for four rab5 genes in zebrafish. Disruption of rab5ab expression by antisense morpholino oligonucleotide (MO) knockdown abolishes nodal signalling in early zebrafish embryos, whereas overexpression of rab5ab mRNA leads to ectopic expression of markers that are normally downstream of nodal signalling. By contrast MO disruption of other zebrafish rab5 genes shows little or no effect on expression of markers of dorsal organiser development. We conclude that rab5ab is essential for nodal signalling and organizer specification in the developing zebrafish embryo., Highlights • We have examined the activities of each of the zebrafish Rab5 genes using morpholino knockdowns. • Loss of one Rab5 isoform, Rab5ab, affects formation of the dorsal organizer. • Rab5ab overexpression leads to ectopic expression of dorsal markers.

Research using Xenopus laevis has made enormous contributions to our understanding of vertebrate development, control of the eukaryotic cell cycle and the cytoskeleton. One limitation, however, has been the lack of systematic genetic studies in Xenopus to complement molecular and cell biological investigations. Work with the closely related diploid frog Xenopus tropicalis is beginning to address this limitation. Here, we review the resources that will make genetic studies using X. tropicalis a reality.

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.

The coatomer vesicular coat complex is essential for normal Golgi and secretory activities in eukaryotic cells. Through positional cloning of genes controlling zebrafish notochord development, we found that the sneezy, happy, and dopey loci encode the alpha, beta, and beta' subunits of the coatomer complex. Export from mutant endoplasmic reticulum is blocked, Golgi structure is disrupted, and mutant embryos eventually degenerate due to widespread apoptosis. The early embryonic phenotype, however, demonstrates that despite its "housekeeping" functions, coatomer activity is specifically and cell autonomously required for normal chordamesoderm differentiation, perinotochordal basement membrane formation, and melanophore pigmentation. Hence, differential requirements for coatomer activity among embryonic tissues lead to tissue-specific developmental defects. Moreover, we note that the mRNA encoding alpha coatomer is strikingly upregulated in notochord progenitors, and we present data suggesting that alpha coatomer transcription is tuned to activity- and cell type-specific secretory loads.

The availability of the whole genome sequence of many model organisms, combined with well-established chemical mutagenesis methods and cost-effective high-throughput DNA genotyping, allows mutations to be identified for virtually any gene. Recently dubbed TILLING (for targeting induced local lesions in genomes), this general method is gaining popularity. In this article, I discuss some of the TILLING methods that are available, the successes that have been reported for several organisms and the future outlook for such methods.

Activities of a variety of signaling proteins that regulate embryogenesis are limited by endogenous antagonists. The zebrafish Nodal-related ligands, Squint and Cyclops, and their antagonists, Lefty1 and Lefty2, belong to the TGFβ-related protein superfamily, whose members have widespread biological activities [1]. Among other activities, Nodals direct the formation of most mesendoderm [2]. By inducing their own transcription and that of the Lefties, Nodal signals establish positive and negative autoregulatory loops [3, 4]. To investigate how these autoregulatory pathways regulate development, we depleted zebrafish embryos of Lefty1 and/or Lefty2 by using antisense morpholino oligonucleotides. Loss of Lefty1 causes aberrations during somitogenesis stages, including left-right patterning defects, whereas Lefty2 depletion has no obvious consequences. Depletion of both Lefty1 and Lefty2, by contrast, causes unchecked Nodal signaling, expansion of mesendoderm, and loss of ectoderm. The expansion of mesendoderm correlates with an extended period of rapid cellular internalization and a failure of deep-cell epiboly. The gastrulation defects of embryos depleted of Lefty1 and Lefty2 result from the deregulation of Squint signaling. In contrast, deregulation of Cyclops does not affect morphology or the transcription of Nodal target genes during gastrulation. Furthermore, we find that Cyclops is specifically required for the maintenance of lefty1 and lefty2 transcription.

Muscular dystrophy is frequently caused by disruption of the dystrophin-glycoprotein complex (DGC), which links muscle cells to the extracellular matrix. Dystroglycan, a central component of the DGC, serves as a laminin receptor via its extracellular α subunit, and interacts with dystrophin (and thus the actin cytoskeleton) through its integral membrane β subunit. We have removed the function of dystroglycan in zebrafish embryos. In contrast to mouse, where dystroglycan mutations lead to peri-implantation lethality, dystroglycan is dispensable for basement membrane formation during early zebrafish development. At later stages, however, loss of dystroglycan leads to a disruption of the DGC, concurrent with loss of muscle integrity and necrosis. In addition, we find that loss of the DGC leads to loss of sarcomere and sarcoplasmic reticulum organisation. The DGC is required for long-term survival of muscle cells in zebrafish, but is dispensable for muscle formation. Dystroglycan or the DGC is also required for normal sarcomere and sarcoplasmic reticulum organisation. Because zebrafish embryos lacking dystroglycan share several characteristics with human muscular dystrophy, they should serve as a useful model for the disease. In addition, knowing the dystroglycan null phenotype in zebrafish will facilitate the isolation of other molecules involved in muscular dystrophy pathogenesis.

Zebrafish embryos homozygous for the masterblind(mbl) mutation exhibit a striking phenotype in which the eyes and telencephalon are reduced or absent and diencephalic fates expand to the front of the brain. Here we show that mbl−/−embryos carry an amino-acid change at a conserved site in the Wnt pathway scaffolding protein, Axin1. The amino-acid substitution present in the mbl allele abolishes the binding of Axin to Gsk3 and affects Tcf-dependent transcription. Therefore, Gsk3 activity may be decreased in mbl−/− embryos and in support of this possibility, overexpression of either wild-type Axin1 or Gsk3β can restore eye and telencephalic fates to mbl−/−embryos. Our data reveal a crucial role for Axin1-dependent inhibition of the Wnt pathway in the early regional subdivision of the anterior neural plate into telencephalic, diencephalic, and eye-forming territories.