Digenic inheritance involving a muscle-specific protein kinase and the giant titin protein causes a skeletal muscle myopathy.

In digenic inheritance, pathogenic variants in two genes must be inherited together to cause disease. Only very few examples of digenic inheritance have been described in the neuromuscular disease field. Here we show that predicted deleterious variants in SRPK3, encoding the X-linked serine/argenine protein kinase 3, lead to a progressive early onset skeletal muscle myopathy only when in combination with heterozygous variants in the TTN gene. The co-occurrence of predicted deleterious SRPK3/TTN variants was not seen among 76,702 healthy male individuals, and statistical modeling strongly supported digenic inheritance as the best-fitting model. Furthermore, double-mutant zebrafish (srpk3-/-; ttn.1+/-) replicated the myopathic phenotype and showed myofibrillar disorganization. Transcriptome data suggest that the interaction of srpk3 and ttn.1 in zebrafish occurs at a post-transcriptional level. We propose that digenic inheritance of deleterious changes impacting both the protein kinase SRPK3 and the giant muscle protein titin causes a skeletal myopathy and might serve as a model for other genetic diseases.

Transcriptional profiling of zebrafish identifies host factors controlling susceptibility to Shigella flexneri.

Shigella flexneri is a human-adapted pathovar of Escherichia coli that can invade the intestinal epithelium, causing inflammation and bacillary dysentery. Although an important human pathogen, the host response to S. flexneri has not been fully described. Zebrafish larvae represent a valuable model for studying human infections in vivo. Here, we use a Shigella-zebrafish infection model to generate mRNA expression profiles of host response to Shigella infection at the whole-animal level. Immune response-related processes dominate the signature of early Shigella infection (6 h post-infection). Consistent with its clearance from the host, the signature of late Shigella infection (24 h post-infection) is significantly changed, and only a small set of immune-related genes remain differentially expressed, including acod1 and gpr84. Using mutant lines generated by ENU, CRISPR mutagenesis and F0 crispants, we show that acod1- and gpr84-deficient larvae are more susceptible to Shigella infection. Together, these results highlight the power of zebrafish to model infection by bacterial pathogens and reveal the mRNA expression of the early (acutely infected) and late (clearing) host response to Shigella infection.

Loss of slc39a14 causes simultaneous manganese hypersensitivity and deficiency in zebrafish.

Manganese neurotoxicity is a hallmark of hypermanganesemia with dystonia 2, an inherited manganese transporter defect caused by mutations in SLC39A14. To identify novel potential targets of manganese neurotoxicity, we performed transcriptome analysis of slc39a14-/- mutant zebrafish that were exposed to MnCl2. Differentially expressed genes mapped to the central nervous system and eye, and pathway analysis suggested that Ca2+ dyshomeostasis and activation of the unfolded protein response are key features of manganese neurotoxicity. Consistent with this interpretation, MnCl2 exposure led to decreased whole-animal Ca2+ levels, locomotor defects and changes in neuronal activity within the telencephalon and optic tectum. In accordance with reduced tectal activity, slc39a14-/- zebrafish showed changes in visual phototransduction gene expression, absence of visual background adaptation and a diminished optokinetic reflex. Finally, numerous differentially expressed genes in mutant larvae normalised upon MnCl2 treatment indicating that, in addition to neurotoxicity, manganese deficiency is present either subcellularly or in specific cells or tissues. Overall, we assembled a comprehensive set of genes that mediate manganese-systemic responses and found a highly correlated and modulated network associated with Ca2+ dyshomeostasis and cellular stress. This article has an associated First Person interview with the first author of the paper.

Translational relevance of forward genetic screens in animal models for the study of psychiatric disease.

Psychiatric disorders represent a significant burden in our societies. Despite the convincing evidence pointing at gene and gene-environment interaction contributions, the role of genetics in the etiology of psychiatric disease is still poorly understood. Forward genetic screens in animal models have helped elucidate causal links. Here we discuss the application of mutagenesis-based forward genetic approaches in common animal model species: two invertebrates, nematodes (Caenorhabditis elegans) and fruit flies (Drosophila sp.); and two vertebrates, zebrafish (Danio rerio) and mice (Mus musculus), in relation to psychiatric disease. We also discuss the use of large scale genomic studies in human populations. Despite the advances using data from human populations, animal models coupled with next-generation sequencing strategies are still needed. Although with its own limitations, zebrafish possess characteristics that make them especially well-suited to forward genetic studies exploring the etiology of psychiatric disorders.

Behavioral and Gene Regulatory Responses to Developmental Drug Exposures in Zebrafish.

Developmental consequences of prenatal drug exposure have been reported in many human cohorts and animal studies. The long-lasting impact on the offspring-including motor and cognitive impairments, cranial and cardiac anomalies and increased prevalence of ADHD-is a socioeconomic burden worldwide. Identifying the molecular changes leading to developmental consequences could help ameliorate the deficits and limit the impact. In this study, we have used zebrafish, a well-established behavioral and genetic model with conserved drug response and reward pathways, to identify changes in behavior and cellular pathways in response to developmental exposure to amphetamine, nicotine or oxycodone. In the presence of the drug, exposed animals showed altered behavior, consistent with effects seen in mammalian systems, including impaired locomotion and altered habituation to acoustic startle. Differences in responses seen following acute and chronic exposure suggest adaptation to the presence of the drug. Transcriptomic analysis of exposed larvae revealed differential expression of numerous genes and alterations in many pathways, including those related to cell death, immunity and circadian rhythm regulation. Differential expression of circadian rhythm genes did not correlate with behavioral changes in the larvae, however, two of the circadian genes, arntl2 and per2, were also differentially expressed at later stages of development, suggesting a long-lasting impact of developmental exposures on circadian gene expression. The immediate-early genes, egr1, egr4, fosab, and junbb, which are associated with synaptic plasticity, were downregulated by all three drugs and in situ hybridization showed that the expression for all four genes was reduced across all neuroanatomical regions, including brain regions implicated in reward processing, addiction and other psychiatric conditions. We anticipate that these early changes in gene expression in response to drug exposure are likely to contribute to the consequences of prenatal exposure and their discovery might pave the way to therapeutic intervention to ameliorate the long-lasting deficits.

The gene regulatory basis of genetic compensation during neural crest induction.

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.

Common and distinct transcriptional signatures of mammalian embryonic lethality.

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.

Compensatory growth renders Tcf7l1a dispensable for eye formation despite its requirement in eye field specification.

The vertebrate eye originates from the eye field, a domain of cells specified by a small number of transcription factors. In this study, we show that Tcf7l1a is one such transcription factor that acts cell-autonomously to specify the eye field in zebrafish. Despite the much-reduced eye field in tcf7l1a mutants, these fish develop normal eyes revealing a striking ability of the eye to recover from a severe early phenotype. This robustness is not mediated through genetic compensation at neural plate stage; instead, the smaller optic vesicle of tcf7l1a mutants shows delayed neurogenesis and continues to grow until it achieves approximately normal size. Although the developing eye is robust to the lack of Tcf7l1a function, it is sensitised to the effects of additional mutations. In support of this, a forward genetic screen identified mutations in hesx1, cct5 and gdf6a, which give synthetically enhanced eye specification or growth phenotypes when in combination with the tcf7l1a mutation.

A high-resolution mRNA expression time course of embryonic development in zebrafish.

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.

Loss of the chromatin modifier Kdm2aa causes BrafV600E-independent spontaneous melanoma in zebrafish.

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.

The age of heterozygous telomerase mutant parents influences the adult phenotype of their offspring irrespective of genotype in zebrafish.

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.

gEVAL - a web-based browser for evaluating genome assemblies.

MOTIVATION: For most research approaches, genome analyses are dependent on the existence of a high quality genome reference assembly. However, the local accuracy of an assembly remains difficult to assess and improve. The gEVAL browser allows the user to interrogate an assembly in any region of the genome by comparing it to different datasets and evaluating the concordance. These analyses include: a wide variety of sequence alignments, comparative analyses of multiple genome assemblies, and consistency with optical and other physical maps. gEVAL highlights allelic variations, regions of low complexity, abnormal coverage, and potential sequence and assembly errors, and offers strategies for improvement. Although gEVAL focuses primarily on sequence integrity, it can also display arbitrary annotation including from Ensembl or TrackHub sources. We provide gEVAL web sites for many human, mouse, zebrafish and chicken assemblies to support the Genome Reference Consortium, and gEVAL is also downloadable to enable its use for any organism and assembly. AVAILABILITY AND IMPLEMENTATION: Web Browser: http://geval.sanger.ac.uk, Plugin: http://wchow.github.io/wtsi-geval-plugin CONTACT: kj2@sanger.ac.uk SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Genome-wide mapping of Hif-1α binding sites in zebrafish.

BACKGROUND: Hypoxia Inducible Factor (HIF) regulates a cascade of transcriptional events in response to decreased oxygenation, acting from the cellular to the physiological level. This response is evolutionarily conserved, allowing the use of zebrafish (Danio rerio) as a model for studying the hypoxic response. Activation of the hypoxic response can be achieved in zebrafish by homozygous null mutation of the von Hippel-Lindau (vhl) tumour suppressor gene. Previous work from our lab has focused on the phenotypic characterisation of this mutant, establishing the links between vhl mutation, the hypoxic response and cancer. To further develop fish as a model for studying hypoxic signalling, we examine the transcriptional profile of the vhl mutant with respect to Hif-1α. As our approach uses embryos consisting of many cell types, it has the potential to uncover additional HIF regulated genes that have escaped detection in analogous mammalian cell culture studies. RESULTS: We performed high-density oligonucleotide microarray analysis of the gene expression changes in von Hippel-Lindau mutant zebrafish, which identified up-regulation of well-known hypoxia response genes and down-regulation of genes primarily involved in lipid processing. To identify the dependency of these transcriptional changes on HIF, we undertook Chromatin Immunoprecipitation linked next generation sequencing (ChIP-seq) for the transcription factor Hypoxia Inducible Factor 1α (HIF-1α). We identified HIF-1α binding sites across the genome, with binding sites showing enrichment for an RCGTG motif, showing conservation with the mammalian hypoxia response element. CONCLUSIONS: Transcriptome analysis of vhl mutant embryos detected activation of key hypoxia response genes seen in human cell models of hypoxia, but also suppression of many genes primarily involved in lipid processing. ChIP-seq analysis of Hif-1α binding sites unveiled an unprecedented number of loci, with a high proportion containing a canonical hypoxia response element. Whether these sites are functional remains unknown, nevertheless their frequent location near transcriptional start sites suggests functionality, and will allow for investigation into the potential hypoxic regulation of genes in their vicinity. We expect that our data will be an excellent starting point for analysis of both fish and mammalian gene regulation by HIF.

High-throughput and quantitative genome-wide messenger RNA sequencing for molecular phenotyping.

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.

Identification of a plant isoflavonoid that causes biliary atresia.

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.

Clinical findings associated with homozygous sickle cell disease in the Barbadian population--do we need a national SCD registry?

BACKGROUND: Comprehensive care in homozygous sickle cell disease (HbSS) entails universal neonatal screening and subsequent monitoring of identified patients, a process which has been streamlined in the neighbouring island of Jamaica. In preparation for a similar undertaking in Barbados, we have developed a database of persons with known HbSS, and have piloted processes for documenting clinical manifestations. We now present a brief clinical profile of these findings with comparisons to the Jamaican cohort. METHODS: HbSS participants were recruited from clinics and support groups. A history of select clinical symptoms was taken and blood and urine samples and echocardiograms were analysed. A re-analysis of data from a previous birth cohort was completed. RESULTS: Forty-eight persons participated (32 F/16 M); age range 10-62 yrs. 94% had a history of ever having a painful crisis. In the past year, 44% of participants had at least one crisis. There were >69 crises in 21 individuals; 61% were self-managed at home and the majority of the others were treated and discharged from hospital; few were admitted. The prevalence of chronic leg ulceration was 27%. Forty-two persons had urinalysis, 44% were diagnosed with albuminuria (urinary protein/creatinine ratio ≥30 mg/g). Thirty-two participants had echocardiography, 28% had a TRJV ≥ 2.5 m/s. Re-analysis of the incidence study revealed a sickle gene frequency (95% CI) of 2.01% (0.24 to 7.21). CONCLUSION: Although we share a common ancestry, it is thought that HbSS is less common and less severe in Barbados compared to Jamaica. The Jamaican studies reported a sickle gene frequency of 3.15 (2.81 to 3.52); the prevalence of chronic leg ulcers and albuminuria was 29.5% and 42.5% respectively. These comparisons suggest that our initial thoughts may be speculative and that HbSS may be an underestimated clinical problem in Barbados. A prospective neonatal screening programme combined with centralized, routine monitoring of HbSS morbidity and outcomes will definitively answer this question and will improve the evidence-based care and management of HbSS in Barbados.

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.

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.

The zebrafish reference genome sequence and its relationship to the human genome.

Zebrafish have become a popular organism for the study of vertebrate gene function. The virtually transparent embryos of this species, and the ability to accelerate genetic studies by gene knockdown or overexpression, have led to the widespread use of zebrafish in the detailed investigation of vertebrate gene function and increasingly, the study of human genetic disease. However, for effective modelling of human genetic disease it is important to understand the extent to which zebrafish genes and gene structures are related to orthologous human genes. To examine this, we generated a high-quality sequence assembly of the zebrafish genome, made up of an overlapping set of completely sequenced large-insert clones that were ordered and oriented using a high-resolution high-density meiotic map. Detailed automatic and manual annotation provides evidence of more than 26,000 protein-coding genes, the largest gene set of any vertebrate so far sequenced. Comparison to the human reference genome shows that approximately 70% of human genes have at least one obvious zebrafish orthologue. In addition, the high quality of this genome assembly provides a clearer understanding of key genomic features such as a unique repeat content, a scarcity of pseudogenes, an enrichment of zebrafish-specific genes on chromosome 4 and chromosomal regions that influence sex determination.

Genomic regulatory blocks encompass multiple neighboring genes and maintain conserved synteny in vertebrates.

We report evidence for a mechanism for the maintenance of long-range conserved synteny across vertebrate genomes. We found the largest mammal-teleost conserved chromosomal segments to be spanned by highly conserved noncoding elements (HCNEs), their developmental regulatory target genes, and phylogenetically and functionally unrelated "bystander" genes. Bystander genes are not specifically under the control of the regulatory elements that drive the target genes and are expressed in patterns that are different from those of the target genes. Reporter insertions distal to zebrafish developmental regulatory genes pax6.1/2, rx3, id1, and fgf8 and miRNA genes mirn9-1 and mirn9-5 recapitulate the expression patterns of these genes even if located inside or beyond bystander genes, suggesting that the regulatory domain of a developmental regulatory gene can extend into and beyond adjacent transcriptional units. We termed these chromosomal segments genomic regulatory blocks (GRBs). After whole genome duplication in teleosts, GRBs, including HCNEs and target genes, were often maintained in both copies, while bystander genes were typically lost from one GRB, strongly suggesting that evolutionary pressure acts to keep the single-copy GRBs of higher vertebrates intact. We show that loss of bystander genes and other mutational events suffered by duplicated GRBs in teleost genomes permits target gene identification and HCNE/target gene assignment. These findings explain the absence of evolutionary breakpoints from large vertebrate chromosomal segments and will aid in the recognition of position effect mutations within human GRBs.