A point mutation decouples the lipid transfer activities of microsomal triglyceride transfer protein.

Apolipoprotein B-containing lipoproteins (B-lps) are essential for the transport of hydrophobic dietary and endogenous lipids through the circulation in vertebrates. Zebrafish embryos produce large numbers of B-lps in the yolk syncytial layer (YSL) to move lipids from yolk to growing tissues. Disruptions in B-lp production perturb yolk morphology, readily allowing for visual identification of mutants with altered B-lp metabolism. Here we report the discovery of a missense mutation in microsomal triglyceride transfer protein (Mtp), a protein that is essential for B-lp production. This mutation of a conserved glycine residue to valine (zebrafish G863V, human G865V) reduces B-lp production and results in yolk opacity due to aberrant accumulation of cytoplasmic lipid droplets in the YSL. However, this phenotype is milder than that of the previously reported L475P stalactite (stl) mutation. MTP transfers lipids, including triglycerides and phospholipids, to apolipoprotein B in the ER for B-lp assembly. In vitro lipid transfer assays reveal that while both MTP mutations eliminate triglyceride transfer activity, the G863V mutant protein unexpectedly retains ~80% of phospholipid transfer activity. This residual phospholipid transfer activity of the G863V mttp mutant protein is sufficient to support the secretion of small B-lps, which prevents intestinal fat malabsorption and growth defects observed in the mttpstl/stl mutant zebrafish. Modeling based on the recent crystal structure of the heterodimeric human MTP complex suggests the G865V mutation may block triglyceride entry into the lipid-binding cavity. Together, these data argue that selective inhibition of MTP triglyceride transfer activity may be a feasible therapeutic approach to treat dyslipidemia and provide structural insight for drug design. These data also highlight the power of yolk transport studies to identify proteins critical for B-lp biology.

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.

Melanoma brain metastases: is it time to eliminate radiotherapy?

PURPOSE: Immunotherapy has demonstrated efficacy in treatment of intracranial metastasis from melanoma, calling into question the role of intracranial radiotherapy (RT). Herein, we assessed the utilization patterns of intracranial RT in patients with melanoma brain metastasis and compared outcomes in patients treated with immunotherapy alone versus immunotherapy in addition to intracranial RT. METHODS: We queried the National Cancer Database (NCDB) for patients with melanoma brain metastases treated with immunotherapy and intracranial RT or immunotherapy alone. Multivariable logistic regression identified variables associated with increased likelihood of receiving immunotherapy alone. Multivariable Cox regression was used to identify co-variates predictive of overall survival (OS). Propensity matching was used to account for indication bias. RESULTS: We identified 528 and 142 patients that were treated with combination therapy and immunotherapy alone, respectively. Patients with lower comorbidity score were more likely to receive immunotherapy alone. Extracranial disease and treatment at a non-academic center were associated with worse OS. Median OS for all patients was 11.0 months. Treatment with stereotactic radiosurgery (SRS) in addition to immunotherapy was superior to immunotherapy alone, median OS of 19.0 versus 11.5 months (p = 0.006). Whole brain radiation therapy (WBRT) in combination with immunotherapy performed worse than immunotherapy alone, median OS of 7.7 versus 11.5 months (p = 0.0255). CONCLUSIONS: For melanoma patients requiring WBRT, immunotherapy alone may be reasonable in asymptomatic patients. For those eligible for SRS, combination therapy may provide better outcomes. Results of ongoing prospective studies will help provide guidance regarding the use of radioimmunotherapy in this population.

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.

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.

KDM2A integrates DNA and histone modification signals through a CXXC/PHD module and direct interaction with HP1.

Functional genomic elements are marked by characteristic DNA and histone modification signatures. How combinatorial chromatin modification states are recognized by epigenetic reader proteins and how this is linked to their biological function is largely unknown. Here we provide a detailed molecular analysis of chromatin recognition by the lysine demethylase KDM2A. Using biochemical approaches we identify a nucleosome interaction module within KDM2A consisting of a CXXC type zinc finger, a PHD domain and a newly identified Heterochromatin Protein 1 (HP1) interaction motif that mediates direct binding between KDM2A and HP1. This nucleosome interaction module enables KDM2A to decode nucleosomal H3K9me3 modification in addition to CpG methylation signals. The multivalent engagement with DNA and HP1 results in a nucleosome binding circuit in which KDM2A can be recruited to H3K9me3-modified chromatin through HP1, and HP1 can be recruited to unmodified chromatin by KDM2A. A KDM2A mutant deficient in HP1-binding is inactive in an in vivo overexpression assay in zebrafish embryos demonstrating that the HP1 interaction is essential for KDM2A function. Our results reveal a complex regulation of chromatin binding for both KDM2A and HP1 that is modulated by DNA- and H3K9-methylation, and suggest a direct role for KDM2A in chromatin silencing.

Efficient identification of CRISPR/Cas9-induced insertions/deletions by direct germline screening in zebrafish.

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.

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.

Flaminal enzyme alginogel: a novel approach to the control of wound exudate, bioburden and debridement.

A multidisciplinary panel of Woundcare experts of international repute was assembled to review the clinical evidence and advise on the classification of the Flaminal products. This is based on their exact role in wound management and is to be defined on the basis of clinical efficacy, evidence and utility. Experts of international repute from Australia, Belgium, Czech Republic, France, Germany, Italy, Netherlands, and UK participated in this exercise.

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.

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.

Zebrafish Fukutin family proteins link the unfolded protein response with dystroglycanopathies.

Allelic mutations in putative glycosyltransferase genes, fukutin and fukutin-related protein (fkrp), lead to a wide range of muscular dystrophies associated with hypoglycosylation of α-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 α-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.

Bilateral adrenal haemorrhage in antiphospholipid syndrome and a short review of the literature.

Antiphospholipid syndrome (APS) is an uncommon autoantibody-mediated condition characterised by acquired thrombophilia resulting in recurrent arterial and venous thrombosis. An inciting factor allows for the exposure of endothelial phospholipids, causing antigen formation and subsequent creation of antibodies. A woman in her 70s presented after vehicular trauma, suffering broken ribs, pneumothorax and incidentally discovered left adrenal haemorrhage. Two weeks later she presented with acute-onset abdominal pain and was found to have a right adrenal gland haemorrhage on CT imaging without interval trauma occurring. The patient had antiphospholipid antibody laboratory studies drawn and was given intravenous heparin with a bridge to warfarin at discharge. Laboratory results returned positive for lupus anticoagulant, beta-2 glycoprotein and anticardiolipin antibodies indicating triple positivity, with repeated laboratory tests positive in 12 weeks' time, confirming the diagnosis. Bilateral adrenal haemorrhage, rather than traditional venous thromboembolism, was the presenting pathology in this patient's diagnosis of APS.

Total Nucleic Acid Extraction from Single Zebrafish Embryos for Genotyping and RNA-seq.

RNA sequencing allows for the quantification of the transcriptome of embryos to investigate transcriptional responses to various perturbations (e.g., mutations, infections, drug treatments). Previous protocols either lack the option to genotype individual samples, or are laborious and therefore difficult to do at a large scale. We have developed a protocol to extract total nucleic acid from individual zebrafish embryos. Individual embryos are lysed in 96-well plates and nucleic acid is extracted using SPRI beads. The total nucleic acid can be genotyped and then DNase I treated to produce RNA samples for sequencing. This protocol allows for processing large numbers of individual samples, with the ability to genotype each sample, which makes it possible to undertake transcriptomic analysis on mutants at timepoints before the phenotype is visible. Graphic abstract: Extraction of total nucleic acid from individual zebrafish embryos for genotyping and RNA-seq.

Allele-specific gene expression can underlie altered transcript abundance in zebrafish mutants.

In model organisms, RNA-sequencing (RNA-seq) is frequently used to assess the effect of genetic mutations on cellular and developmental processes. Typically, animals heterozygous for a mutation are crossed to produce offspring with different genotypes. Resultant embryos are grouped by genotype to compare homozygous mutant embryos to heterozygous and wild-type siblings. Genes that are differentially expressed between the groups are assumed to reveal insights into the pathways affected by the mutation. Here we show that in zebrafish, differentially expressed genes are often over-represented on the same chromosome as the mutation due to different levels of expression of alleles from different genetic backgrounds. Using an incross of haplotype-resolved wild-type fish, we found evidence of widespread allele-specific expression, which appears as differential expression when comparing embryos homozygous for a region of the genome to their siblings. When analysing mutant transcriptomes, this means that the differential expression of genes on the same chromosome as a mutation of interest may not be caused by that mutation. Typically, the genomic location of a differentially expressed gene is not considered when interpreting its importance with respect to the phenotype. This could lead to pathways being erroneously implicated or overlooked due to the noise of spurious differentially expressed genes on the same chromosome as the mutation. These observations have implications for the interpretation of RNA-seq experiments involving outbred animals and non-inbred model organisms.

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.

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.

Complex regulation of cyp26a1 creates a robust retinoic acid gradient in the zebrafish embryo.

Positional identities along the anterior-posterior axis of the vertebrate nervous system are assigned during gastrulation by multiple posteriorizing signals, including retinoic acid (RA), fibroblast growth factors (Fgfs), and Wnts. Experimental evidence has suggested that RA, which is produced in paraxial mesoderm posterior to the hindbrain by aldehyde dehydrogenase 1a2 (aldh1a2/raldh2), forms a posterior-to-anterior gradient across the hindbrain field, and provides the positional information that specifies the locations and fates of rhombomeres. Recently, alternative models have been proposed in which RA plays only a permissive role, signaling wherever it is not degraded. Here we use a combination of experimental and modeling tools to address the role of RA in providing long-range positional cues in the zebrafish hindbrain. Using cell transplantation and implantation of RA-coated beads into RA-deficient zebrafish embryos, we demonstrate that RA can directly convey graded positional information over long distances. We also show that expression of Cyp26a1, the major RA-degrading enzyme during gastrulation, is under complex feedback and feedforward control by RA and Fgf signaling. The predicted consequence of such control is that RA gradients will be both robust to fluctuations in RA synthesis and adaptive to changes in embryo length during gastrulation. Such control also provides an explanation for the fact that loss of an endogenous RA gradient can be compensated for by RA that is provided in a spatially uniform manner.