Novel cardiovascular gene functions revealed via systematic phenotype prediction in zebrafish

Distinct roles of core autophagy-related genes in zebrafish definitive hematopoiesis

Despite the well-described discrepancy between ATG (macroautophagy/autophagy-related) genes in the regulation of hematopoiesis, varying essentiality of core ATG proteins in vertebrate definitive hematopoiesis remains largely unclear. Here, we employed zebrafish (Danio rerio) to compare the functions of six core atg genes, including atg13, becn1 (beclin1), atg9a, atg2a, atg5, and atg3, in vertebrate definitive hematopoiesis via clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9 ribonucleoprotein and morpholino targeting. Zebrafish with various atg mutations showed autophagic deficiency and presented partially consistent hematopoietic abnormalities during early development. All six atg mutations led to a declined number of spi1b+ (Spi-1 proto-oncogene b) myeloid progenitor cells. However, only becn1 mutation resulted in the expansion of myb+ (v-myb avian myeloblastosis viral oncogene homolog) hematopoietic stem and progenitor cells (HSPCs) and transiently increased coro1a+ (coronin, actin binding protein, 1A) leukocytes, whereas atg3 mutation decreased the number of HSPCs and leukocytes. Proteomic analysis of caudal hematopoietic tissue identified sin3aa (SIN3 transcription regulator family member Aa) as a potential modulator of atg13- and becn1-regulated definitive hematopoiesis. Disruption of sin3aa rescued the expansion of HSPCs and leukocytes in becn1 mutants and exacerbated the decrease of HSPCs in atg13 mutants. Double mutations were also performed to examine alternative functions of various atg genes in definitive hematopoiesis. Notably, becn1 mutation failed to induce HSPCs expansion with one of the other five atg mutations. These findings demonstrated the distinct roles of atg genes and their interplays in zebrafish definitive hematopoiesis, thereby suggesting that the vertebrate definitive hematopoiesis is regulated in an atg gene-dependent manner.Abbreviations: AGM: aorta-gonad-mesonephros; AO: acridine orange; atg: autophagy related; becn1: beclin 1, autophagy related; CHT: caudal hematopoietic tissue; CKO: conditional knockout; coro1a: coronin, actin binding protein, 1A; CQ: chloroquine; CRISPR: clustered regularly interspaced short palindromic repeats; dpf: days post fertilization; FACS: fluorescence-activated cell sorting; hbae1.1: hemoglobin, alpha embryonic 1.1; HSCs: hematopoietic stem cells; HSPCs: hematopoietic stem and progenitor cells; KD: knockdown; KO: knockout; map1lc3/lc3: microtubule-associated protein 1 light chain 3; MO: morpholino; mpeg1.1: macrophage expressed 1, tandem duplicate 1; mpx: myeloid-specific peroxidase; myb: v-myb avian myeloblastosis viral oncogene homolog; PE: phosphatidylethanolamine; p-H3: phospho-H3 histone; PtdIns3K: class 3 phosphatidylinositol 3-kinase; rag1: recombination activating 1; rb1cc1/fip200: RB1-inducible coiled-coil 1; RFLP: restriction fragment length polymorphism; RNP: ribonucleoprotein; sin3aa: SIN3 transcription regulator family member Aa; spi1b: Spi-1 proto-oncogene b; ulk: unc-51 like autophagy activating kinase; vtg1: vitellogenin 1; WISH: whole-mount in situ hybridization.

Zebrafish as a model organism - can a fish mimic human?

From pre-historic era, all scientific discoveries have evolved around a concept - THINK BIG but for a change zebrafish as a model organism in research had managed to halt the entire medical community and made us realize that it's time to think small. From a barely imagined being in research few years ago to around 4,000 publications in just last year, zebrafish has definitely come a long way. Through these tiny fish, scientists have managed to find genes that caused human diseases and have also developed various specific models to know more about the pathology behind such diseases. This review will focus on zebrafish as a model organism from the time it was introduced to the most novel targets with particular emphasis on central nervous system (CNS) as it is rapidly evolving branch in zebrafish research these days. This review will try to shed light on the early stages of zebrafish as a model organism and will try to cover the journey of it developing as a successful model organism to map many diseases like diabetes, Alzheimer's and autism describing the rationale for using this specific model and briefly the techniques under each category and finally will summarize the pros and cons of the model with its expected future directions.

Zebrafish as a Mainstream Model for In Vivo Systems Pharmacology and Toxicology

Pharmacology and toxicology are part of a much broader effort to understand the relationship between chemistry and biology. While biomedicine has necessarily focused on specific cases, typically of direct human relevance, there are real advantages in pursuing more systematic approaches to characterizing how health and disease are influenced by small molecules and other interventions. In this context, the zebrafish is now established as the representative screenable vertebrate and, through ongoing advances in the available scale of genome editing and automated phenotyping, is beginning to address systems-level solutions to some biomedical problems. The addition of broader efforts to integrate information content across preclinical model organisms and the incorporation of rigorous analytics, including closed-loop deep learning, will facilitate efforts to create systems pharmacology and toxicology with the ability to continuously optimize chemical biological interactions around societal needs. In this review, we outline progress toward this goal.

The structural basis of histone modifying enzyme specificity and promiscuity: Implications for metabolic regulation and drug design

Histone modifying enzymes regulate chromatin architecture through covalent modifications and ultimately control multiple aspects of cellular function. Disruption of histone modification leads to changes in gene expression profiles and may lead to disease. Both small molecule inhibitors and intermediary metabolites have been shown to modulate histone modifying enzyme activity although our ability to identify successful drug candidates or novel metabolic regulators of these enzymes has been limited. Using a combination of large scale in silico screens and in vivo phenotypic analysis, we identified several small molecules and intermediary metabolites with distinctive HME activity. Our approach using unsupervised learning identifies the chemical fingerprints of both small molecules and metabolites that facilitate recognition by the enzymes active sites which can be used as a blueprint to design novel inhibitors. Furthermore, this work supports the idea that histone modifying enzymes sense intermediary metabolites integrating genes, environment and cellular physiology.

In vivo identification and validation of novel potential predictors for human cardiovascular diseases

Genetics crucially contributes to cardiovascular diseases (CVDs), the global leading cause of death. Since the majority of CVDs can be prevented by early intervention there is a high demand for the identification of predictive causative genes. While genome wide association studies (GWAS) correlate genes and CVDs after diagnosis and provide a valuable resource for such causative candidate genes, often preferentially those with previously known or suspected function are addressed further. To tackle the unaddressed blind spot of understudied genes, we particularly focused on the validation of human heart phenotype-associated GWAS candidates with little or no apparent connection to cardiac function. Building on the conservation of basic heart function and underlying genetics from fish to human we combined CRISPR/Cas9 genome editing of the orthologs of human GWAS candidates in isogenic medaka with automated high-throughput heart rate analysis. Our functional analyses of understudied human candidates uncovered a prominent fraction of heart rate associated genes from adult human patients impacting on the heart rate in embryonic medaka already in the injected generation. Following this pipeline, we identified 16 GWAS candidates with potential diagnostic and predictive power for human CVDs.

GeneMANIA update 2018

GeneMANIA (http://genemania.org) is a flexible user-friendly web site for generating hypotheses about gene function, analyzing gene lists and prioritizing genes for functional assays. Given a query gene list, GeneMANIA finds functionally similar genes using a wealth of genomics and proteomics data. In this mode, it weights each functional genomic dataset according to its predictive value for the query. Another use of GeneMANIA is gene function prediction. Given a single query gene, GeneMANIA finds genes likely to share function with it based on their interactions with it. Enriched Gene Ontology categories among this set can point to the function of the gene. Nine organisms are currently supported (Arabidopsis thaliana, Caenorhabditis elegans, Danio rerio, Drosophila melanogaster, Escherichia coli, Homo sapiens, Mus musculus, Rattus norvegicus and Saccharomyces cerevisiae). Hundreds of data sets and hundreds of millions of interactions have been collected from GEO, BioGRID, IRefIndex and I2D, as well as organism-specific functional genomics data sets. Users can customize their search by selecting specific data sets to query and by uploading their own data sets to analyze. We have recently updated the user interface to GeneMANIA to make it more intuitive and make more efficient use of visual space. GeneMANIA can now be used effectively on a variety of devices.

Macrophages target Salmonella by Lc3-associated phagocytosis in a systemic infection model

Innate immune defense against intracellular pathogens, like Salmonella, relies heavily on the autophagy machinery of the host. This response is studied intensively in epithelial cells, the target of Salmonella during gastrointestinal infections. However, little is known of the role that autophagy plays in macrophages, the predominant carriers of this pathogen during systemic disease. Here we utilize a zebrafish embryo model to study the interaction of S. enterica serovar Typhimurium with the macroautophagy/autophagy machinery of macrophages in vivo. We show that phagocytosis of live but not heat-killed Salmonella triggers recruitment of the autophagy marker GFP-Lc3 in a variety of patterns labeling tight or spacious bacteria-containing compartments, also revealed by electron microscopy. Neutrophils display similar GFP-Lc3 associations, but genetic modulation of the neutrophil/macrophage balance and ablation experiments show that macrophages are critical for the defense response. Deficiency of atg5 reduces GFP-Lc3 recruitment and impairs host resistance, in contrast to atg13 deficiency, indicating that Lc3-Salmonella association at this stage is independent of the autophagy preinitiation complex and that macrophages target Salmonella by Lc3-associated phagocytosis (LAP). In agreement, GFP-Lc3 recruitment and host resistance are impaired by deficiency of Rubcn/Rubicon, known as a negative regulator of canonical autophagy and an inducer of LAP. We also found strict dependency on NADPH oxidase, another essential factor for LAP. Both Rubcn and NADPH oxidase are required to activate a Salmonella biosensor for reactive oxygen species inside infected macrophages. These results identify LAP as the major host protective autophagy-related pathway responsible for macrophage defense against Salmonella during systemic infection. Abbreviations: ATG: autophagy related gene; BECN1: Beclin 1; CFU: colony forming units; CYBA/P22PHOX: cytochrome b-245, alpha chain; CYBB/NOX2: cytochrome b-245 beta chain; dpf: days post fertilization; EGFP: enhanced green fluorescent protein; GFP: green fluorescent protein; hfp: hours post fertilization; hpi: hours post infection; IRF8: interferon regulatory factor 8; Lcp1/L-plastin: lymphocyte cytosolic protein 1; LAP: LC3-associated phagocytosis; MAP1LC3/LC3: microtubule-associated protein 1A/1B-light chain 3; mCherry: red fluorescent protein; mpeg1: macrophage expressed gene 1; mpx: myeloid specific peroxidase; NADPH oxidase: nicotinamide adenine dinucleotide phosphate oxidase; NCF4/P40PHOX: neutrophil cytosolic factor 4; NTR-mCherry: nitroreductase-mCherry fusion; PTU: phenylthiourea; PtdIns3K: class III phosphatidylinositol 3-kinase; PtdIns3P: phosphatidylinositol 3-phosphate; RB1CC1/FIP200: RB-1 inducible coiled coin 1; ROS: reactive oxygen species; RT-PCR: reverse transcriptase polymerase chain reaction; RUBCN/RUBICON: RUN and cysteine rich domain containing BECN1-interacting protein; SCV: Salmonella-containing vacuole; S. Typhimurium/S.T: Salmonella enterica serovar Typhimurium; TEM: transmission electron microscopy; Tg: transgenic; TSA: tyramide signal amplification; ULK1/2: unc-51-like autophagy activating kinase 1/2; UVRAG: UVRAG: UV radiation resistance associated; wt: wild type.

Cardiac injections of AntagomiRs as a novel tool for knockdown of miRNAs during heart development

BACKGROUND

Studying microRNA networks during heart development is essential to obtain a better understanding of developmental defects and diseases associated with the heart and to identify novel opportunities for therapeutics. Here we highlight the advantages of chicken embryos as a vertebrate model for studying intermediate processes of heart development. Avians develop a four-chambered heart closely resembling human anatomy and they develop ex utero, which makes them easily accessible. Furthermore, embryos are available all year with a steady supply.

RESULTS

In this report we established a novel method for the knockdown of microRNA function by microinjecting AntagomiRs into the chicken heart in ovo. Our approach enables the targeted delivery of antagomirs into a locally restricted area and is not impacted by circulation. After further embryo development the successful miRNA knockdown was confirmed. Loss of function phenotypes can be evaluated rapidly, compared to more time-consuming genetic ablation experiments. The local application avoids potential systemic effects of microRNA knockdown, therefore allowing the assessment of impacts on heart development only. The method can be adjusted for different stages of chicken embryos (HH13-HH18) as well as for knockdown or targeted overexpression of coding genes.

CONCLUSION

In conclusion our method allows targeted and locally restricted delivery of Antagomirs to the heart leading to successful knockdown of microRNA function. This method enables rapid phenotypic assessment, for example by gene expression analysis of multiple cardiac genes.

Prediction of Human Phenotype Ontology terms by means of hierarchical ensemble methods

BACKGROUND

The prediction of human gene-abnormal phenotype associations is a fundamental step toward the discovery of novel genes associated with human disorders, especially when no genes are known to be associated with a specific disease. In this context the Human Phenotype Ontology (HPO) provides a standard categorization of the abnormalities associated with human diseases. While the problem of the prediction of gene-disease associations has been widely investigated, the related problem of gene-phenotypic feature (i.e., HPO term) associations has been largely overlooked, even if for most human genes no HPO term associations are known and despite the increasing application of the HPO to relevant medical problems. Moreover most of the methods proposed in literature are not able to capture the hierarchical relationships between HPO terms, thus resulting in inconsistent and relatively inaccurate predictions.

RESULTS

We present two hierarchical ensemble methods that we formally prove to provide biologically consistent predictions according to the hierarchical structure of the HPO. The modular structure of the proposed methods, that consists in a "flat" learning first step and a hierarchical combination of the predictions in the second step, allows the predictions of virtually any flat learning method to be enhanced. The experimental results show that hierarchical ensemble methods are able to predict novel associations between genes and abnormal phenotypes with results that are competitive with state-of-the-art algorithms and with a significant reduction of the computational complexity.

CONCLUSIONS

Hierarchical ensembles are efficient computational methods that guarantee biologically meaningful predictions that obey the true path rule, and can be used as a tool to improve and make consistent the HPO terms predictions starting from virtually any flat learning method. The implementation of the proposed methods is available as an R package from the CRAN repository.

The Cyprinodon variegatus genome reveals gene expression changes underlying differences in skull morphology among closely related species

BACKGROUND

Understanding the genetic and developmental origins of phenotypic novelty is central to the study of biological diversity. In this study we identify modifications to the expression of genes at four developmental stages that may underlie jaw morphological differences among three closely related species of pupfish (genus Cyprinodon) from San Salvador Island, Bahamas. Pupfishes on San Salvador Island are trophically differentiated and include two endemic species that have evolved jaw morphologies unlike that of any other species in the genus Cyprinodon.

RESULTS

We find that gene expression differs significantly across recently diverged species of pupfish. Genes such as Bmp4 and calmodulin, previously implicated in jaw diversification in African cichlid fishes and Galapagos finches, were not found to be differentially expressed among species of pupfish. Instead we find multiple growth factors and cytokine/chemokine genes to be differentially expressed among these pupfish taxa. These include both genes and pathways known to affect craniofacial development, such as Wnt signaling, as well as novel genes and pathways not previously implicated in craniofacial development. These data highlight both shared and potentially unique sources of jaw diversity in pupfish and those identified in other evolutionary model systems such as Galapagos finches and African cichlids.

CONCLUSIONS

We identify modifications to the expression of genes involved in Wnt signaling, Igf signaling, and the inflammation response as promising avenues for future research. Our project provides insight into the magnitude of gene expression changes contributing to the evolution of morphological novelties, such as jaw structure, in recently diverged pupfish species.

Comparison of Zebrafish tmem88a mutant and morpholino knockdown phenotypes

Tmem88a is a transmembrane protein that is thought to be a negative regulator of the Wnt signalling pathway. Several groups have used antisense morpholino oligonucleotides in an effort to characterise the role of tmem88a in zebrafish cardiovascular development, but they have not obtained consistent results. Here, we generate an 8 bp deletion in the coding region of the tmem88a locus using TALENs, and we have gone on to establish a viable homozygous tmem88aΔ8 mutant line. Although tmem88aΔ8 mutants have reduced expression of some key haematopoietic genes, differentiation of erythrocytes and neutrophils is unaffected, contradicting our previous study using antisense morpholino oligonucleotides. We find that expression of the tmem88a paralogue tmem88b is not significantly changed in tmem88aΔ8 mutants and injection of the tmem88a splice-blocking morpholino oligonucleotide into tmem88aΔ8 mutants recapitulates the reduction of erythrocytes observed in morphants using o-Dianisidine. This suggests that there is a partial, but inessential, requirement for tmem88a during haematopoiesis and that morpholino injection exacerbates this phenotype in tmem88a morpholino knockdown embryos.

Keratin mediated attachment of stem cells to augment cardiomyogenic lineage commitment

The objective of this work was to develop a simple surface modification technique using keratin derived from human hair for efficient cardiomyogenic lineage commitment of human mesenchymal stem cells (hMSCs). Keratin was extracted from discarded human hair containing both the acidic and basic components along with the heterodimers. The extracted keratin was adsorbed to conventional tissue culture polystyrene surfaces at different concentration. Keratin solution of 500μg/ml yielded a well coated layer of 12±1nm thickness with minimal agglomeration. The keratin coated surfaces promoted cell attachment and proliferation. Large increases in the mRNA expression of known cardiomyocyte genes such as cardiac actinin, cardiac troponin and β-myosin heavy chain were observed. Immunostaining revealed increased expression of sarcomeric α-actinin and tropomyosin whereas Western blots confirmed higher expression of tropomyosin and myocyte enhancer factor 2C in cells on the keratin coated surface than on the non-coated surface. Keratin promoted DNA demethylation of the Atp2a2 and Nkx2.5 genes thereby elucidating the importance of epigenetic changes as a possible molecular mechanism underlying the increased differentiation. A global gene expression analysis revealed a significant alteration in the expression of genes involved in pathways associated in cardiomyogenic commitment including cytokine and chemokine signaling, cell-cell and cell-matrix interactions, Wnt signaling, MAPK signaling, TGF-β signaling and FGF signaling pathways among others. Thus, adsorption of keratin offers a facile and affordable yet potent route for inducing cardiomyogenic lineage commitment of stem cells with important implications in developing xeno-free strategies in cardiovascular regenerative medicine.

Function-driven discovery of disease genes in zebrafish using an integrated genomics big data resource

Whole exome sequencing (WES) accelerates disease gene discovery using rare genetic variants, but further statistical and functional evidence is required to avoid false-discovery. To complement variant-driven disease gene discovery, here we present function-driven disease gene discovery in zebrafish (Danio rerio), a promising human disease model owing to its high anatomical and genomic similarity to humans. To facilitate zebrafish-based function-driven disease gene discovery, we developed a genome-scale co-functional network of zebrafish genes, DanioNet (www.inetbio.org/danionet), which was constructed by Bayesian integration of genomics big data. Rigorous statistical assessment confirmed the high prediction capacity of DanioNet for a wide variety of human diseases. To demonstrate the feasibility of the function-driven disease gene discovery using DanioNet, we predicted genes for ciliopathies and performed experimental validation for eight candidate genes. We also validated the existence of heterozygous rare variants in the candidate genes of individuals with ciliopathies yet not in controls derived from the UK10K consortium, suggesting that these variants are potentially involved in enhancing the risk of ciliopathies. These results showed that an integrated genomics big data for a model animal of diseases can expand our opportunity for harnessing WES data in disease gene discovery.

An extended set of yeast-based functional assays accurately identifies human disease mutations

We can now routinely identify coding variants within individual human genomes. A pressing challenge is to determine which variants disrupt the function of disease-associated genes. Both experimental and computational methods exist to predict pathogenicity of human genetic variation. However, a systematic performance comparison between them has been lacking. Therefore, we developed and exploited a panel of 26 yeast-based functional complementation assays to measure the impact of 179 variants (101 disease- and 78 non-disease-associated variants) from 22 human disease genes. Using the resulting reference standard, we show that experimental functional assays in a 1-billion-year diverged model organism can identify pathogenic alleles with significantly higher precision and specificity than current computational methods.

Metallothionein 2 regulates endothelial cell migration through transcriptional regulation of vegfc expression

Analysis of developmental angiogenesis can help to identify regulatory networks, which also contribute to disease-related vascular growth. Vascular endothelial growth factors (Vegf) drive angiogenic processes such as sprouting, endothelial cell (EC) migration and proliferation. However, how Vegf expression is regulated during development is not well understood. By analyzing developmental zebrafish angiogenesis, we have identified Metallothionein 2 (Mt2) as a novel regulator of vegfc expression. While Metallothioneins (Mts) have been extensively analyzed for their capability of regulating homeostasis and metal detoxification, we demonstrate that Mt2 is required for EC migration, proliferation and angiogenic sprouting upstream of vegfc expression. We further demonstrate that another Mt family member cannot compensate Mt2 deficiency and therefore postulate that Mt2 regulates angiogenesis independent of its canonical Mt function. Our data not only reveal a non-canonical function of Mt2 in angiogenesis, but also propose Mt2 as a novel regulator of vegfc expression.

Cardiac transcriptome and dilated cardiomyopathy genes in zebrafish

BACKGROUND

Genetic studies of cardiomyopathy and heart failure have limited throughput in mammalian models. Adult zebrafish have been recently pursued as a vertebrate model with higher throughput, but genetic conservation must be tested.

METHODS AND RESULTS

We conducted transcriptome analysis of zebrafish heart and searched for fish homologues of 51 known human dilated cardiomyopathy-associated genes. We also identified genes with high cardiac expression and genes with differential expression between embryonic and adult stages. Among tested genes, 30 had a single zebrafish orthologue, 14 had 2 homologues, and 5 had ≥3 homologues. By analyzing the expression data on the basis of cardiac abundance and enrichment hypotheses, we identified a single zebrafish gene for 14 of 19 multiple-homologue genes and 2 zebrafish homologues of high priority for ACTC1. Of note, our data suggested vmhc and vmhcl as functional zebrafish orthologues for human genes MYH6 and MYH7, respectively, which are established molecular markers for cardiac remodeling.

CONCLUSIONS

Most known genes for human dilated cardiomyopathy have a corresponding zebrafish orthologue, which supports the use of zebrafish as a conserved vertebrate model. Definition of the cardiac transcriptome and fetal gene program will facilitate systems biology studies of dilated cardiomyopathy in zebrafish.

Selecting causal genes from genome-wide association studies via functionally coherent subnetworks

Genome-wide association (GWA) studies have linked thousands of loci to human diseases, but the causal genes and variants at these loci generally remain unknown. Although investigators typically focus on genes closest to the associated polymorphisms, the causal gene is often more distal. Reliance on published work to prioritize candidates is biased toward well-characterized genes. We describe a 'prix fixe' strategy and software that uses genome-scale shared-function networks to identify sets of mutually functionally related genes spanning multiple GWA loci. Using associations from ∼100 GWA studies covering ten cancer types, our approach outperformed the common alternative strategy in ranking known cancer genes. As more GWA loci are discovered, the strategy will have increased power to elucidate the causes of human disease.

Prioritizing causal disease genes using unbiased genomic features

BACKGROUND

Cardiovascular disease (CVD) is the leading cause of death in the developed world. Human genetic studies, including genome-wide sequencing and SNP-array approaches, promise to reveal disease genes and mechanisms representing new therapeutic targets. In practice, however, identification of the actual genes contributing to disease pathogenesis has lagged behind identification of associated loci, thus limiting the clinical benefits.

RESULTS

To aid in localizing causal genes, we develop a machine learning approach, Objective Prioritization for Enhanced Novelty (OPEN), which quantitatively prioritizes gene-disease associations based on a diverse group of genomic features. This approach uses only unbiased predictive features and thus is not hampered by a preference towards previously well-characterized genes. We demonstrate success in identifying genetic determinants for CVD-related traits, including cholesterol levels, blood pressure, and conduction system and cardiomyopathy phenotypes. Using OPEN, we prioritize genes, including FLNC, for association with increased left ventricular diameter, which is a defining feature of a prevalent cardiovascular disorder, dilated cardiomyopathy or DCM. Using a zebrafish model, we experimentally validate FLNC and identify a novel FLNC splice-site mutation in a patient with severe DCM.

CONCLUSION

Our approach stands to assist interpretation of large-scale genetic studies without compromising their fundamentally unbiased nature.

Small heat shock proteins are necessary for heart migration and laterality determination in zebrafish

Small heat shock proteins (sHsps) regulate cellular functions not only under stress, but also during normal development, when they are expressed in organ-specific patterns. Here we demonstrate that two small heat shock proteins expressed in embryonic zebrafish heart, hspb7 and hspb12, have roles in the development of left-right asymmetry. In zebrafish, laterality is determined by the motility of cilia in Kupffer's vesicle (KV), where hspb7 is expressed; knockdown of hspb7 causes laterality defects by disrupting the motility of these cilia. In embryos with reduced hspb7, the axonemes of KV cilia have a 9+0 structure, while control embyros have a predominately 9+2 structure. Reduction of either hspb7 or hspb12 alters the expression pattern of genes that propagate the signals that establish left-right asymmetry: the nodal-related gene southpaw (spaw) in the lateral plate mesoderm, and its downstream targets pitx2, lefty1 and lefty2. Partial depletion of hspb7 causes concordant heart, brain and visceral laterality defects, indicating that loss of KV cilia motility leads to coordinated but randomized laterality. Reducing hspb12 leads to similar alterations in the expression of downstream laterality genes, but at a lower penetrance. Simultaneous reduction of hspb7 and hspb12 randomizes heart, brain and visceral laterality, suggesting that these two genes have partially redundant functions in the establishment of left-right asymmetry. In addition, both hspb7 and hspb12 are expressed in the precardiac mesoderm and in the yolk syncytial layer, which supports the migration and fusion of mesodermal cardiac precursors. In embryos in which the reduction of hspb7 or hspb12 was limited to the yolk, migration defects predominated, suggesting that the yolk expression of these genes rather than heart expression is responsible for the migration defects.