Extracellular protein interactions are crucial to the development of multicellular organisms because they initiate signaling pathways and enable cellular recognition cues. Despite their importance, extracellular protein interactions are often under-represented in large scale protein interaction data sets because most high throughput assays are not designed to detect low affinity extracellular interactions. Due to the lack of a comprehensive data set, the evolution of extracellular signaling pathways has remained largely a mystery. We investigated this question using a combined data set of physical pairwise interactions between zebrafish extracellular proteins, mainly from the immunoglobulin superfamily and leucine-rich repeat families, and their spatiotemporal expression profiles. We took advantage of known homology between proteins to estimate the relative rates of changes of four parameters after gene duplication, namely extracellular protein interaction, expression pattern, and the divergence of extracellular and intracellular protein sequences. We showed that change in expression profile is a major contributor to the evolution of signaling pathways followed by divergence in intracellular protein sequence, whereas extracellular sequence and interaction profiles were relatively more conserved. Rapidly evolving expression profiles will eventually drive other parameters to diverge more quickly because differentially expressed proteins get exposed to different environments and potential binding partners. This allows homologous extracellular receptors to attain specialized functions and become specific to tissues and/or developmental stages.

The dorsal habenular nuclei of the zebrafish epithalamus have become a valuable model for studying the development of left-right (L-R) asymmetry and its function in the vertebrate brain. The bilaterally paired dorsal habenulae exhibit striking differences in size, neuroanatomical organization, and molecular properties. They also display differences in their efferent connections with the interpeduncular nucleus (IPN) and in their afferent input, with a subset of mitral cells distributed on both sides of the olfactory bulb innervating only the right habenula. Previous studies have implicated the dorsal habenulae in modulating fear/anxiety responses in juvenile and adult zebrafish. It has been suggested that the asymmetric olfactory-habenula pathway (OB-Ha), revealed by selective labeling from an lhx2a:YFP transgene, mediates fear behaviors elicited by alarm pheromone. Here we show that expression of the fam84b gene demarcates a unique region of the right habenula that is the site of innervation by lhx2a:YFP-labeled olfactory axons. Upon ablation of the parapineal, which normally promotes left habenular identity; the fam84b domain is present in both dorsal habenulae and lhx2a:YFP-labeled olfactory bulb neurons form synapses on the left and the right side. To explore the relevance of the asymmetric olfactory projection and how it might influence habenular function, we tested activation of this pathway using odorants known to evoke behaviors. We find that alarm substance or other aversive odors, and attractive cues, activate fos expression in subsets of cells in the olfactory bulb but not in the lhx2a:YFP expressing population. Moreover, neither alarm pheromone nor chondroitin sulfate elicited fos activation in the dorsal habenulae. The results indicate that L-R asymmetry of the epithalamus sets the directionality of olfactory innervation, however, the lhx2a:YFP OB-Ha pathway does not appear to mediate fear responses to aversive odorants.

The cellular retinoic acid-binding protein 2 (CRABP2) is believed to be involved in regulating access of retinoic acid to nuclear retinoic acid receptors. We have determined the cDNA sequence and the genomic organization of the duplicated crabp2 gene (crabp2b) in zebrafish. The crabp2b cDNA was 522bp in length and encodes a polypeptide consisting of 146 amino acids. Radiation hybrid mapping assigned the crabp2b gene to zebrafish linkage group 19. The comparison of the mapped human CRABP2 gene, zebrafish crabp2a and zebrafish crabp2b genes revealed that human chromosome 1 has a syntenic relationship to zebrafish linkage groups 16 and 19. Reverse transcription-polymerase chain reaction (RT-PCR) detected crabp2b mRNA in total RNA extracted from whole adult zebrafish, but not in any of the adult zebrafish tissues examined. The crabp2a mRNA was detected in total RNA extracted from whole adult zebrafish, adult zebrafish muscle, testes, and skin and to a lesser extent in heart, ovary and brain. No crabp2a mRNA-specific product was detected in kidney, liver or intestine of the adult zebrafish. Whole mount in situ hybridization detected crabp2b and crabp2a mRNA in a number of structures known to require retinoic acid signaling during embryonic development. The crabp2b mRNA was detected in the central nervous system, branchial arches, pectoral fins, retina (dorsal to the lens), epidermis and otic vesicle of the developing zebrafish. The crabp2a transcripts were detected by whole mount in situ hybridization in the central nervous system, epidermis, proliferative zone of the retina, intestinal bulb, oesophagus, pectoral fins and branchial arches during zebrafish embryonic development.

Mutations in the SEPN1 gene encoding the selenoprotein N (SelN) have been described in different congenital myopathies. Here, we report the first mutation in the selenocysteine insertion sequence (SECIS) of SelN messenger RNA, a hairpin structure located in the 3' untranslated region, in a patient presenting a classical although mild form of rigid spine muscular dystrophy. We detected a significant reduction in both mRNA and protein levels in the patient's skin fibroblasts. The SECIS element is crucial for the insertion of selenocysteine at the reprogrammed UGA codon by recruiting the SECIS-binding protein 2 (SBP2), and we demonstrated that this mutation abolishes SBP2 binding to SECIS in vitro, thereby preventing co-translational incorporation of selenocysteine and SelN synthesis. The identification of this mutation affecting a conserved base in the SECIS functional motif thereby reveals the structural basis for a novel pathological mechanism leading to SEPN1-related myopathy.

In many vertebrates, establishment of Left-Right (LR) asymmetry results from the activity of a ciliated organ functioning as the LR Organizer (LRO). While regulation of the formation of this structure by major signaling pathways has been described, the transcriptional control of LRO formation is poorly understood. Using the zebrafish model, we show that the transcription factors and cofactors mediating or regulating the transcriptional outcome of the Hippo signaling pathway play a pivotal role in controlling the expression of genes essential to the formation of the LRO including ligands and receptors of signaling pathways involved in this process and most genes required for motile ciliogenesis. Moreover, the transcription cofactor, Vgll4l regulates epigenetic programming in LRO progenitors by controlling the expression of writers and readers of DNA methylation marks. Altogether, our study uncovers a novel and essential role for the transcriptional effectors and regulators of the Hippo pathway in establishing LR asymmetry.

Selenium is important for embryogenesis in vertebrates but little is known about the expression patterns and biological functions of most selenoprotein genes. Taking advantage of the zebrafish model, systematic analysis of selenoprotein gene expression was performed by in situ hybridization on whole-mount embryos at different developmental stages. Twenty-one selenoprotein mRNAs were analyzed and all of them exhibited expression patterns restricted to specific tissues. Moreover, we demonstrated that highly similar selenoprotein paralogs were expressed within distinct territories. Therefore, tissue- and development-specific expression patterns provided new information for selenoproteins of unknown function.

To gain an understanding of molecular events that underlie pattern formation in the retina, we evaluated the expression profiles of over 8000 transcripts randomly selected from an embryonic zebrafish library. Detailed analysis of cDNAs that display restricted expression patterns revealed factors that are specifically expressed in single cell classes and are potential regulators of neurogenesis. These cDNAs belong to numerous molecular categories and include cell surface receptors, cytoplasmic enzymes, and transcription factors. To test whether expression patterns that we have uncovered using this approach are indicative of function in neurogenesis, we used morpholino-mediated knockdown approach. The knockdown of soxp, a transcript expressed in the vicinity of the inner plexiform layer, revealed its role in cell type composition of amacrine and ganglion cell layers. Blocking the function of cxcr4b, a chemokine receptor specifically expressed in ganglion cells, suggests a role in ganglion cell survival. These experiments demonstrate that in situ hybridization-based reverse genetic screens can be applied to isolate genetic regulators of neurogenesis. This approach very well complements forward genetic mutagenesis studies previously used to study retinal neurogenesis in zebrafish.

Nuclear receptors (NRs) are transcription factors that are implicated in several biological processes such as embryonic development, homeostasis, and metabolic diseases. To study the role of NRs in development, it is critically important to know when and where individual genes are expressed. Although systematic expression studies using reverse transcriptase PCR and/or DNA microarrays have been performed in classical model systems such as Drosophila and mouse, no systematic atlas describing NR involvement during embryonic development on a global scale has been assembled. Adopting a systems biology approach, we conducted a systematic analysis of the dynamic spatiotemporal expression of all NR genes as well as their main transcriptional coregulators during zebrafish development (101 genes) using whole-mount in situ hybridization. This extensive dataset establishes overlapping expression patterns among NRs and coregulators, indicating hierarchical transcriptional networks. This complete developmental profiling provides an unprecedented examination of expression of NRs during embryogenesis, uncovering their potential function during central nervous system and retina formation. Moreover, our study reveals that tissue specificity of hormone action is conferred more by the receptors than by their coregulators. Finally, further evolutionary analyses of this global resource led us to propose that neofunctionalization of duplicated genes occurs at the levels of both protein sequence and RNA expression patterns. Altogether, this expression database of NRs provides novel routes for leading investigation into the biological function of each individual NR as well as for the study of their combinatorial regulatory circuitry within the superfamily.

Mutations in the gene coding for selenoprotein N (SelN), a selenium containing protein of unknown function, cause different forms of congenital muscular dystrophy in humans. These muscular diseases are characterized by early onset of hypotonia which predominantly affect in axial muscles. We used zebrafish as a model system to understand the function of SelN in muscle formation during embryogenesis. Zebrafish SelN is highly homologous to its human counterpart and amino acids corresponding to the mutated positions in human muscle diseases are conserved in the zebrafish protein. The sepn1 gene is highly expressed in the somites and notochord during early development. Inhibition of the sepn1 gene by injection of antisense morpholinos does not alter the fate of the muscular tissue, but causes muscle architecture disorganization and greatly reduced motility. Ultrastructural analysis of the myotomes reveals defects in muscle sarcomeric organization and in myofibers attachment, as well as altered myoseptum integrity. These studies demonstrate the important role of SelN for muscle organization during early development. Moreover, alteration of myofibrils architecture and tendon-like structure in embryo deficient for SelN function provide new insights into the pathological mechanism of SelN-related myopathy.

Sialyltransferases are responsible for the synthesis of a diverse range of sialoglycoconjugates predicted to be pivotal to deuterostomes' evolution. In this work, we reconstructed the evolutionary history of the metazoan α2,3-sialyltransferases family (ST3Gal), a subset of sialyltransferases encompassing six subfamilies (ST3Gal I-ST3Gal VI) functionally characterized in mammals. Exploration of genomic and expressed sequence tag databases and search of conserved sialylmotifs led to the identification of a large data set of st3gal-related gene sequences. Molecular phylogeny and large scale sequence similarity network analysis identified four new vertebrate subfamilies called ST3Gal III-r, ST3Gal VII, ST3Gal VIII, and ST3Gal IX. To address the issue of the origin and evolutionary relationships of the st3gal-related genes, we performed comparative syntenic mapping of st3gal gene loci combined to ancestral genome reconstruction. The ten vertebrate ST3Gal subfamilies originated from genome duplication events at the base of vertebrates and are organized in three distinct and ancient groups of genes predating the early deuterostomes. Inferring st3gal gene family history identified also several lineage-specific gene losses, the significance of which was explored in a functional context. Toward this aim, spatiotemporal distribution of st3gal genes was analyzed in zebrafish and bovine tissues. In addition, molecular evolutionary analyses using specificity determining position and coevolved amino acid predictions led to the identification of amino acid residues with potential implication in functional divergence of vertebrate ST3Gal. We propose a detailed scenario of the evolutionary relationships of st3gal genes coupled to a conceptual framework of the evolution of ST3Gal functions.

A gene for the zebrafish brain-type fatty acid-binding protein (fabp7b) was identified and its structure defined. The zebrafish fabp7b gene spans 1479 bp and consists of four exons encoding 24, 58, 34 and 16 amino acids, respectively, which is identical to the structure of the fabp7a gene previously described. The complete fabp7b cDNA was isolated by 5' and 3' RACE and its nucleotide sequence determined. The deduced amino acid sequence of FABP7B encoded by the zebrafish fabp7b gene shares 82% identity with that of FABP7A encoded by the zebrafish fabp7a gene. A single transcription start site for the fabp7b gene was mapped by 5' RNA ligase-mediated RACE. Phylogenetic analysis indicated that the duplication of the fabp7 genes occurred in the fish lineage after their divergence from mammals. The zebrafish fabp7b gene was assigned to linkage group 20 by radiation hybrid mapping. Reverse transcription-polymerase chain reaction detected fabp7b transcripts in the same adult tissues as fabp7a transcripts. In the brain, levels of fabp7b transcripts were lower than fabp7a transcripts. Whole-mount in situ hybridization showed that the zebrafish fabp7a transcripts were distributed in the early developing central nervous system. In addition to being expressed in the developing brain and retina, zebrafish fabp7b mRNA was also detected in the swim bladder and pharynx during the embryonic to larval transitory phase.

Selenium is an essential trace element and selenoprotein N (SelN) was the first selenium-containing protein shown to be directly involved in human inherited diseases. Mutations in the SEPN1 gene, encoding SelN, cause a group of muscular disorders characterized by predominant affection of axial muscles. SelN has been shown to participate in calcium and redox homeostasis, but its pathophysiological role in skeletal muscle remains largely unknown. To address SelN function in vivo, we generated a Sepn1-null mouse model by gene targeting. The Sepn1(-/-) mice had normal growth and lifespan, and were macroscopically indistinguishable from wild-type littermates. Only minor defects were observed in muscle morphology and contractile properties in SelN-deficient mice in basal conditions. However, when subjected to challenging physical exercise and stress conditions (forced swimming test), Sepn1(-/-) mice developed an obvious phenotype, characterized by limited motility and body rigidity during the swimming session, as well as a progressive curvature of the spine and predominant alteration of paravertebral muscles. This induced phenotype recapitulates the distribution of muscle involvement in patients with SEPN1-Related Myopathy, hence positioning this new animal model as a valuable tool to dissect the role of SelN in muscle function and to characterize the pathophysiological process.

Zebrafish (Danio rerio) can serve as a model system to study heritable skin diseases. The skin is rapidly developed during the first 5-6 days of embryonic growth, accompanied by expression of skin-specific genes. Transmission electron microscopy (TEM) of wild-type zebrafish at day 5 reveals a two-cell-layer epidermis separated from the underlying collagenous stroma by a basement membrane with fully developed hemidesmosomes. Scanning electron microscopy (SEM) reveals an ordered surface contour of keratinocytes with discrete microridges. To gain insight into epidermal morphogenesis, we have employed morpholino-mediated knockdown of the abca12 and snap29 genes, which are crucial for secretion of lipids and intracellular trafficking of lamellar granules, respectively. Morpholinos, when placed on exon-intron junctions, were >90% effective in preventing the corresponding gene expression when injected into one- to four-cell-stage embryos. By day 3, TEM of abca12 morphants showed accumulation of lipid-containing electron-dense lamellar granules, whereas snap29 morphants showed the presence of apparently empty vesicles in the epidermis. Evaluation of epidermal morphogenesis by SEM revealed similar perturbations in both cases in the microridge architecture and the development of spicule-like protrusions on the surface of keratinocytes. These morphological findings are akin to epidermal changes in harlequin ichthyosis and CEDNIK syndrome, autosomal recessive keratinization disorders due to mutations in the ABCA12 and SNAP29 genes, respectively. The results indicate that interference of independent pathways involving lipid transport in the epidermis can result in phenotypically similar perturbations in epidermal morphogenesis, and that these fish mutants can serve as a model to study the pathomechanisms of these keratinization disorders.

Phosphopantothenate is a precursor to synthesis of coenzyme A, a molecule essential to many metabolic pathways. Organisms of the archaeal phyla were shown to utilize a different phosphopantothenate biosynthetic pathway from the eukaryotic and bacterial one. In this study, we report that symbiotic bacteria from the group Candidatus poribacteria present enzymes of the archaeal pathway, namely pantoate kinase and phosphopantothenate synthetase, mirroring what was demonstrated for Picrophilus torridus, an archaea partially utilizing the bacterial pathway. Our results not only support the ancient origin of the coenzyme A pathway in the three domains of life but also highlight its complex and dynamic evolution. Importantly, this study helps to improve protein annotation for this pathway in the C. poribacteria group and other related organisms.

Despite the fundamental importance of patterning along the dorsal-ventral (DV) and anterior-posterior (AP) axes during embryogenesis, uncertainty exists in the orientation of these axes for the mesoderm. Here we examine the origin and formation of the zebrafish kidney, a ventrolateral mesoderm derivative, and show that AP patterning of the non-axial mesoderm occurs across the classic gastrula stage DV axis while DV patterning aligns along the animal-vegetal pole. We find that BMP signalling acts early to establish broad anterior and posterior territories in the non-axial mesoderm while retinoic acid (RA) functions later, but also across the classic DV axis. Our data support a model in which RA on the dorsal side of the embryo induces anterior kidney fates while posterior kidney progenitors are protected ventrally by the RA-catabolizing enzyme Cyp26a1. This work clarifies our understanding of vertebrate axis orientation and establishes a new paradigm for how the kidney and other mesodermal derivatives arise during embryogenesis.

Genomic and cDNA sequences coding for two cellular retinol-binding proteins (rbp) in zebrafish were retrieved from DNA sequence databases. Phylogenetic analysis revealed that these proteins were most similar to mammalian RBP7/Rbp7 proteins. Hence, the genes coding for these proteins were named rbp7a and rbp7b. Using a radiation hybrid panel, rbp7a and rbp7b were mapped to the zebrafish chromosomes 23 and 6, respectively. Conserved gene synteny indicated that these genes most likely arose as a result of a fish-specific whole-genome duplication event that had occurred 230-400 million years ago. Whole-mount in situ hybridization to zebrafish embryos detected rbp7a transcripts from the sphere stage (4h post-fertilization (hpf)) in the forerunner cells and the yolk syncytial layer, as well as in Kuppfer's vesicle and the periderm at 12 hpf. The transcripts of rbp7b were seen primarily in the somite stages (10-24 hpf) of zebrafish embryos, but also in the floor plate and hypochord, and did not overlap with the distribution of rbp7a transcripts in embryos. The hybridization signal for rbp7a and rbp7b transcripts was not detected in embryos after 12 hpf and 24 hpf, respectively. While transcripts for both rbp7a and rbp7b were found in all adult tissues assayed by RT-qPCR, the steady-state level of rbp7a transcripts were significantly higher than that of rbp7b transcripts in gill and ovary, whereas rbp7b transcripts were significantly higher than rbp7a transcripts in muscle and brain. The distribution of rbp7a and rbp7b transcripts in embryos and adult zebrafish indicate that the cis-elements that control the transcriptional regulation of the rbp7a and rbp7b genes have diverged considerably since their duplication.

Extracellular interactions involving both secreted and membrane-tethered receptor proteins are essential to initiate signaling pathways that orchestrate cellular behaviors within biological systems. Because of the biochemical properties of these proteins and their interactions, identifying novel extracellular interactions remains experimentally challenging. To address this, we have recently developed an assay, AVEXIS (avidity-based extracellular interaction screen) to detect low affinity extracellular interactions on a large scale and have begun to construct interaction networks between zebrafish receptors belonging to the immunoglobulin and leucine-rich repeat protein families to identify novel signaling pathways important for early development. Here, we expanded our zebrafish protein library to include other domain families and many more secreted proteins and performed our largest screen to date totaling 16,544 potential unique interactions. We report 111 interactions of which 96 are novel and include the first documented extracellular ligands for 15 proteins. By including 77 interactions from previous screens, we assembled an expanded network of 188 extracellular interactions between 92 proteins and used it to show that secreted proteins have twice as many interaction partners as membrane-tethered receptors and that the connectivity of the extracellular network behaves as a power law. To try to understand the functional role of these interactions, we determined new expression patterns for 164 genes within our clone library by using whole embryo in situ hybridization at five key stages of zebrafish embryonic development. These expression data were integrated with the binding network to reveal where each interaction was likely to function within the embryo and were used to resolve the static interaction network into dynamic tissue- and stage-specific subnetworks within the developing zebrafish embryo. All these data were organized into a freely accessible on-line database called ARNIE (AVEXIS Receptor Network with Integrated Expression; www.sanger.ac.uk/arnie) and provide a valuable resource of new extracellular signaling interactions for developmental biology.

The endoplasmic reticulum (ER) is the reservoir for calcium in cells. Luminal calcium levels are determined by calcium-sensing proteins that trigger calcium dynamics in response to calcium fluctuations. Here we report that Selenoprotein N (SEPN1) is a type II transmembrane protein that senses ER calcium fluctuations by binding this ion through a luminal EF-hand domain. In vitro and in vivo experiments show that via this domain, SEPN1 responds to diminished luminal calcium levels, dynamically changing its oligomeric state and enhancing its redox-dependent interaction with cellular partners, including the ER calcium pump sarcoplasmic/endoplasmic reticulum calcium ATPase (SERCA). Importantly, single amino acid substitutions in the EF-hand domain of SEPN1 identified as clinical variations are shown to impair its calcium-binding and calcium-dependent structural changes, suggesting a key role of the EF-hand domain in SEPN1 function. In conclusion, SEPN1 is a ER calcium sensor that responds to luminal calcium depletion, changing its oligomeric state and acting as a reductase to refill ER calcium stores.

Generating properly differentiated embryonic structures in vitro from pluripotent stem cells remains a challenge. Here we show that instruction of aggregates of mouse embryonic stem cells with an experimentally engineered morphogen signalling centre, that functions as an organizer, results in the development of embryo-like entities (embryoids). In situ hybridization, immunolabelling, cell tracking and transcriptomic analyses show that these embryoids form the three germ layers through a gastrulation process and that they exhibit a wide range of developmental structures, highly similar to neurula-stage mouse embryos. Embryoids are organized around an axial chordamesoderm, with a dorsal neural plate that displays histological properties similar to the murine embryo neuroepithelium and that folds into a neural tube patterned antero-posteriorly from the posterior midbrain to the tip of the tail. Lateral to the chordamesoderm, embryoids display somitic and intermediate mesoderm, with beating cardiac tissue anteriorly and formation of a vasculature network. Ventrally, embryoids differentiate a primitive gut tube, which is patterned both antero-posteriorly and dorso-ventrally. Altogether, embryoids provide an in vitro model of mammalian embryo that displays extensive development of germ layer derivatives and that promises to be a powerful tool for in vitro studies and disease modelling.

PTEN is a tumor suppressor gene associated with multiple tumor types. PTEN function is essential for early embryonic development and is involved in the regulation of cell size, number, and survival. By dephosphorylating PIP(3), PTEN normally acts to inhibit the PI3-Kinase/AKT pathway. Here we have identified two zebrafish orthologs, ptena and ptenb, of the single mammalian PTEN gene and analyzed the role of these genes in zebrafish development. Ptena transcripts were expressed throughout the embryo at early somitogenesis. By 24 hpf, expression was predominant in the central nervous system, axial vasculature, retina, branchial arches, ear, lateral line primordium, and pectoral fin bud. Ptenb was also ubiquitously expressed early in somitogenesis, but transcripts became more restricted to the somites and central nervous system as development progressed. By 48 hpf, ptena and ptenb were expressed predominantly in the central nervous system, branchial arches, pectoral fins, and eye. Antisense morpholinos were used to knock down translation of ptena and ptenb mRNA in zebrafish embryos. Knockdown of either pten gene caused increased levels of phosphorylated Akt in morphant embryos, indicating that Ptena and Ptenb each possess PIP(3) lipid phosphatase activity. Ptena morphants had irregularities in notochord shape (73%), vasculogenesis (83%), head shape (72%), and inner ear development (59%). The most noticeable defects in ptenb morphants were upward hooked tails (73%), domed heads (83%), and reduced yolk extensions (90%). These results indicate that ptena and ptenb encode functional enzymes and that each pten gene plays a distinct role during zebrafish embryogenesis.