COPB2 loss of function causes a coatopathy with osteoporosis and developmental delay.

Coatomer complexes function in the sorting and trafficking of proteins between subcellular organelles. Pathogenic variants in coatomer subunits or associated factors have been reported in multi-systemic disorders, i.e., coatopathies, that can affect the skeletal and central nervous systems. We have identified loss-of-function variants in COPB2, a component of the coatomer complex I (COPI), in individuals presenting with osteoporosis, fractures, and developmental delay of variable severity. Electron microscopy of COPB2-deficient subjects' fibroblasts showed dilated endoplasmic reticulum (ER) with granular material, prominent rough ER, and vacuoles, consistent with an intracellular trafficking defect. We studied the effect of COPB2 deficiency on collagen trafficking because of the critical role of collagen secretion in bone biology. COPB2 siRNA-treated fibroblasts showed delayed collagen secretion with retention of type I collagen in the ER and Golgi and altered distribution of Golgi markers. copb2-null zebrafish embryos showed retention of type II collagen, disorganization of the ER and Golgi, and early larval lethality. Copb2+/- mice exhibited low bone mass, and consistent with the findings in human cells and zebrafish, studies in Copb2+/- mouse fibroblasts suggest ER stress and a Golgi defect. Interestingly, ascorbic acid treatment partially rescued the zebrafish developmental phenotype and the cellular phenotype in Copb2+/- mouse fibroblasts. This work identifies a form of coatopathy due to COPB2 haploinsufficiency, explores a potential therapeutic approach for this disorder, and highlights the role of the COPI complex as a regulator of skeletal homeostasis.

Heterozygous loss-of-function variants significantly expand the phenotypes associated with loss of GDF11.

PURPOSE: Growth differentiation factor 11 (GDF11) is a key signaling protein required for proper development of many organ systems. Only one prior study has associated an inherited GDF11 variant with a dominant human disease in a family with variable craniofacial and vertebral abnormalities. Here, we expand the phenotypic spectrum associated with GDF11 variants and document the nature of the variants. METHODS: We present a cohort of six probands with de novo and inherited nonsense/frameshift (4/6 patients) and missense (2/6) variants in GDF11. We generated gdf11 mutant zebrafish to model loss of gdf11 phenotypes and used an overexpression screen in Drosophila to test variant functionality. RESULTS: Patients with variants in GDF11 presented with craniofacial (5/6), vertebral (5/6), neurological (6/6), visual (4/6), cardiac (3/6), auditory (3/6), and connective tissue abnormalities (3/6). gdf11 mutant zebrafish show craniofacial abnormalities and body segmentation defects that match some patient phenotypes. Expression of the patients' variants in the fly showed that one nonsense variant in GDF11 is a severe loss-of-function (LOF) allele whereas the missense variants in our cohort are partial LOF variants. CONCLUSION: GDF11 is needed for human development, particularly neuronal development, and LOF GDF11 alleles can affect the development of numerous organs and tissues.

A fish with no sex: gonadal and adrenal functions partition between zebrafish NR5A1 co-orthologs.

People with NR5A1 mutations experience testicular dysgenesis, ovotestes, or adrenal insufficiency, but we do not completely understand the origin of this phenotypic diversity. NR5A1 is expressed in gonadal soma precursor cells before expression of the sex-determining gene SRY. Many fish have two co-orthologs of NR5A1 that likely partitioned ancestral gene subfunctions between them. To explore ancestral roles of NR5A1, we knocked out nr5a1a and nr5a1b in zebrafish. Single-cell RNA-seq identified nr5a1a-expressing cells that co-expressed genes for steroid biosynthesis and the chemokine receptor Cxcl12a in 1-day postfertilization (dpf) embryos, as does the mammalian adrenal-gonadal (interrenal-gonadal) primordium. In 2dpf embryos, nr5a1a was expressed stronger in the interrenal-gonadal primordium than in the early hypothalamus but nr5a1b showed the reverse. Adult Leydig cells expressed both ohnologs and granulosa cells expressed nr5a1a stronger than nr5a1b. Mutants for nr5a1a lacked the interrenal, formed incompletely differentiated testes, had no Leydig cells, and grew far larger than normal fish. Mutants for nr5a1b formed a disorganized interrenal and their gonads completely disappeared. All homozygous mutant genotypes lacked secondary sex characteristics, including male breeding tubercles and female sex papillae, and had exceedingly low levels of estradiol, 11-ketotestosterone, and cortisol. RNA-seq showed that at 21dpf, some animals were developing as females and others were not, independent of nr5a1 genotype. By 35dpf, all mutant genotypes greatly under-expressed ovary-biased genes. Because adult nr5a1a mutants form gonads but lack an interrenal and conversely, adult nr5a1b mutants lack a gonad but have an interrenal, the adrenal, and gonadal functions of the ancestral nr5a1 gene partitioned between ohnologs after the teleost genome duplication, likely owing to reciprocal loss of ancestral tissue-specific regulatory elements. Identifying such elements could provide hints to otherwise unexplained cases of Differences in Sex Development.

BICRA, a SWI/SNF Complex Member, Is Associated with BAF-Disorder Related Phenotypes in Humans and Model Organisms.

SWI/SNF-related intellectual disability disorders (SSRIDDs) are rare neurodevelopmental disorders characterized by developmental disability, coarse facial features, and fifth digit/nail hypoplasia that are caused by pathogenic variants in genes that encode for members of the SWI/SNF (or BAF) family of chromatin remodeling complexes. We have identified 12 individuals with rare variants (10 loss-of-function, 2 missense) in the BICRA (BRD4 interacting chromatin remodeling complex-associated protein) gene, also known as GLTSCR1, which encodes a subunit of the non-canonical BAF (ncBAF) complex. These individuals exhibited neurodevelopmental phenotypes that include developmental delay, intellectual disability, autism spectrum disorder, and behavioral abnormalities as well as dysmorphic features. Notably, the majority of individuals lack the fifth digit/nail hypoplasia phenotype, a hallmark of most SSRIDDs. To confirm the role of BICRA in the development of these phenotypes, we performed functional characterization of the zebrafish and Drosophila orthologs of BICRA. In zebrafish, a mutation of bicra that mimics one of the loss-of-function variants leads to craniofacial defects possibly akin to the dysmorphic facial features seen in individuals harboring putatively pathogenic BICRA variants. We further show that Bicra physically binds to other non-canonical ncBAF complex members, including the BRD9/7 ortholog, CG7154, and is the defining member of the ncBAF complex in flies. Like other SWI/SNF complex members, loss of Bicra function in flies acts as a dominant enhancer of position effect variegation but in a more context-specific manner. We conclude that haploinsufficiency of BICRA leads to a unique SSRIDD in humans whose phenotypes overlap with those previously reported.

Correction: yippee like 3 (ypel3) is a novel gene required for myelinating and perineurial glia development.

[This corrects the article DOI: 10.1371/journal.pgen.1008841.].

yippee like 3 (ypel3) is a novel gene required for myelinating and perineurial glia development.

Hypomyelination, a neurological condition characterized by decreased production of myelin sheets by glial cells, often has no known etiology. Elucidating the genetic causes of hypomyelination provides a better understanding of myelination, as well as means to diagnose, council, and treat patients. Here, we present evidence that YIPPEE LIKE 3 (YPEL3), a gene whose developmental role was previously unknown, is required for central and peripheral glial cell development. We identified a child with a constellation of clinical features including cerebral hypomyelination, abnormal peripheral nerve conduction, hypotonia, areflexia, and hypertrophic peripheral nerves. Exome and genome sequencing revealed a de novo mutation that creates a frameshift in the open reading frame of YPEL3, leading to an early stop codon. We used zebrafish as a model system to validate that YPEL3 mutations are causative of neuropathy. We found that ypel3 is expressed in the zebrafish central and peripheral nervous system. Using CRISPR/Cas9 technology, we created zebrafish mutants carrying a genomic lesion similar to that of the patient. Our analysis revealed that Ypel3 is required for development of oligodendrocyte precursor cells, timely exit of the perineurial glial precursors from the central nervous system (CNS), formation of the perineurium, and Schwann cell maturation. Consistent with these observations, zebrafish ypel3 mutants have metabolomic signatures characteristic of oligodendrocyte and Schwann cell differentiation defects, show decreased levels of Myelin basic protein in the central and peripheral nervous system, and develop defasciculated peripheral nerves. Locomotion defects were observed in adult zebrafish ypel3 mutants. These studies demonstrate that Ypel3 is a novel gene required for perineurial cell development and glial myelination.

Bi-allelic Variants in TONSL Cause SPONASTRIME Dysplasia and a Spectrum of Skeletal Dysplasia Phenotypes.

SPONASTRIME dysplasia is an autosomal-recessive spondyloepimetaphyseal dysplasia characterized by spine (spondylar) abnormalities, midface hypoplasia with a depressed nasal bridge, metaphyseal striations, and disproportionate short stature. Scoliosis, coxa vara, childhood cataracts, short dental roots, and hypogammaglobulinemia have also been reported in this disorder. Although an autosomal-recessive inheritance pattern has been hypothesized, pathogenic variants in a specific gene have not been discovered in individuals with SPONASTRIME dysplasia. Here, we identified bi-allelic variants in TONSL, which encodes the Tonsoku-like DNA repair protein, in nine subjects (from eight families) with SPONASTRIME dysplasia, and four subjects (from three families) with short stature of varied severity and spondylometaphyseal dysplasia with or without immunologic and hematologic abnormalities, but no definitive metaphyseal striations at diagnosis. The finding of early embryonic lethality in a Tonsl-/- murine model and the discovery of reduced length, spinal abnormalities, reduced numbers of neutrophils, and early lethality in a tonsl-/- zebrafish model both support the hypomorphic nature of the identified TONSL variants. Moreover, functional studies revealed increased amounts of spontaneous replication fork stalling and chromosomal aberrations, as well as fewer camptothecin (CPT)-induced RAD51 foci in subject-derived cell lines. Importantly, these cellular defects were rescued upon re-expression of wild-type (WT) TONSL; this rescue is consistent with the hypothesis that hypomorphic TONSL variants are pathogenic. Overall, our studies in humans, mice, zebrafish, and subject-derived cell lines confirm that pathogenic variants in TONSL impair DNA replication and homologous recombination-dependent repair processes, and they lead to a spectrum of skeletal dysplasia phenotypes with numerous extra-skeletal manifestations.

Grxcr1 Promotes Hair Bundle Development by Destabilizing the Physical Interaction between Harmonin and Sans Usher Syndrome Proteins.

Morphogenesis and mechanoelectrical transduction of the hair cell mechanoreceptor depend on the correct assembly of Usher syndrome (USH) proteins into highly organized macromolecular complexes. Defects in these proteins lead to deafness and vestibular areflexia in USH patients. Mutations in a non-USH protein, glutaredoxin domain-containing cysteine-rich 1 (GRXCR1), cause non-syndromic sensorineural deafness. To understand the deglutathionylating enzyme function of GRXCR1 in deafness, we generated two grxcr1 zebrafish mutant alleles. We found that hair bundles are thinner in homozygous grxcr1 mutants, similar to the USH1 mutants ush1c (Harmonin) and ush1ga (Sans). In vitro assays showed that glutathionylation promotes the interaction between Ush1c and Ush1ga and that Grxcr1 regulates mechanoreceptor development by preventing physical interaction between these proteins without affecting the assembly of another USH1 protein complex, the Ush1c-Cadherin23-Myosin7aa tripartite complex. By elucidating the molecular mechanism through which Grxcr1 functions, we also identify a mechanism that dynamically regulates the formation of Usher protein complexes.

A Recurrent De Novo Heterozygous COG4 Substitution Leads to Saul-Wilson Syndrome, Disrupted Vesicular Trafficking, and Altered Proteoglycan Glycosylation.

The conserved oligomeric Golgi (COG) complex is involved in intracellular vesicular transport, and is composed of eight subunits distributed in two lobes, lobe A (COG1-4) and lobe B (COG5-8). We describe fourteen individuals with Saul-Wilson syndrome, a rare form of primordial dwarfism with characteristic facial and radiographic features. All affected subjects harbored heterozygous de novo variants in COG4, giving rise to the same recurrent amino acid substitution (p.Gly516Arg). Affected individuals' fibroblasts, whose COG4 mRNA and protein were not decreased, exhibited delayed anterograde vesicular trafficking from the ER to the Golgi and accelerated retrograde vesicular recycling from the Golgi to the ER. This altered steady-state equilibrium led to a decrease in Golgi volume, as well as morphologic abnormalities with collapse of the Golgi stacks. Despite these abnormalities of the Golgi apparatus, protein glycosylation in sera and fibroblasts from affected subjects was not notably altered, but decorin, a proteoglycan secreted into the extracellular matrix, showed altered Golgi-dependent glycosylation. In summary, we define a specific heterozygous COG4 substitution as the molecular basis of Saul-Wilson syndrome, a rare skeletal dysplasia distinct from biallelic COG4-CDG.

Usherin defects lead to early-onset retinal dysfunction in zebrafish.

Mutations in USH2A are the most frequent cause of Usher syndrome and autosomal recessive nonsyndromic retinitis pigmentosa. To unravel the pathogenic mechanisms underlying USH2A-associated retinal degeneration and to evaluate future therapeutic strategies that could potentially halt the progression of this devastating disorder, an animal model is needed. The available Ush2a knock-out mouse model does not mimic the human phenotype, because it presents with only a mild and late-onset retinal degeneration. Using CRISPR/Cas9-technology, we introduced protein-truncating germline lesions into the zebrafish ush2a gene (ush2armc1: c.2337_2342delinsAC; p.Cys780GlnfsTer32 and ush2ab1245: c.15520_15523delinsTG; p.Ala5174fsTer). Homozygous mutants were viable and displayed no obvious morphological or developmental defects. Immunohistochemical analyses with antibodies recognizing the N- or C-terminal region of the ush2a-encoded protein, usherin, demonstrated complete absence of usherin in photoreceptors of ush2armc1, but presence of the ectodomain of usherin at the periciliary membrane of ush2ab1245-derived photoreceptors. Furthermore, defects of usherin led to a reduction in localization of USH2 complex members, whirlin and Adgrv1, at the photoreceptor periciliary membrane of both mutants. Significantly elevated levels of apoptotic photoreceptors could be observed in both mutants when kept under constant bright illumination for three days. Electroretinogram (ERG) recordings revealed a significant and similar decrease in both a- and b-wave amplitudes in ush2armc1 as well as ush2ab1245 larvae as compared to strain- and age-matched wild-type larvae. In conclusion, this study shows that mutant ush2a zebrafish models present with early-onset retinal dysfunction that is exacerbated by light exposure. These models provide a better understanding of the pathophysiology underlying USH2A-associated RP and a unique opportunity to evaluate future therapeutic strategies.

Mutation of POC1B in a severe syndromic retinal ciliopathy.

We describe a consanguineous Iraqi family with Leber congenital amaurosis (LCA), Joubert syndrome (JBTS), and polycystic kidney disease (PKD). Targeted next-generation sequencing for excluding mutations in known LCA and JBTS genes, homozygosity mapping, and whole-exome sequencing identified a homozygous missense variant, c.317G>C (p.Arg106Pro), in POC1B, a gene essential for ciliogenesis, basal body, and centrosome integrity. In silico modeling suggested a requirement of p.Arg106 for the formation of the third WD40 repeat and a protein interaction interface. In human and mouse retina, POC1B localized to the basal body and centriole adjacent to the connecting cilium of photoreceptors and in synapses of the outer plexiform layer. Knockdown of Poc1b in zebrafish caused cystic kidneys and retinal degeneration with shortened and reduced photoreceptor connecting cilia, compatible with the human syndromic ciliopathy. A recent study describes homozygosity for p.Arg106ProPOC1B in a family with nonsyndromic cone-rod dystrophy. The phenotype associated with homozygous p.Arg106ProPOC1B may thus be highly variable, analogous to homozygous p.Leu710Ser in WDR19 causing either isolated retinitis pigmentosa or Jeune syndrome. Our study indicates that POC1B is required for retinal integrity, and we propose POC1B mutations as a probable cause for JBTS with severe PKD.

The cone-dominant retina and the inner ear of zebrafish express the ortholog of CLRN1, the causative gene of human Usher syndrome type 3A.

Clarin-1 (CLRN1) is the causative gene in Usher syndrome type 3A, an autosomal recessive disorder characterized by progressive vision and hearing loss. CLRN1 encodes Clarin-1, a glycoprotein with homology to the tetraspanin family of proteins. Previous cell culture studies suggest that Clarin-1 localizes to the plasma membrane and interacts with the cytoskeleton. Mouse models demonstrate a role for the protein in mechanosensory hair bundle integrity, but the function of Clarin-1 in hearing remains unclear. Even less is known of its role in vision, because the Clrn1 knockout mouse does not exhibit a retinal phenotype and expression studies in murine retinas have provided conflicting results. Here, we describe cloning and expression analysis of the zebrafish clrn1 gene, and report protein localization of Clarin-1 in auditory and visual cells from embryonic through adult stages. We detect clrn1 transcripts as early as 24h post-fertilization, and expression is maintained through adulthood. In situ hybridization experiments show clrn1 transcripts enriched in mechanosensory hair cells and supporting cells of the inner ear and lateral line organ, photoreceptors, and cells of the inner retina. In mechanosensory hair cells, Clarin-1 is polarized to the apical cell body and the synapses. In the retina, Clarin-1 localizes to lateral cell contacts between photoreceptors and is associated with the outer limiting membrane and subapical processes emanating from Müller glial cells. We also find Clarin-1 protein in the outer plexiform, inner nuclear and ganglion cell layers of the retina. Given the importance of Clarin-1 function in the human retina, it is imperative to find an animal model with a comparable requirement. Our data provide a foundation for exploring the role of Clarin-1 in retinal cell function and survival in a diurnal, cone-dominant species.

Identification of differentially expressed genes in the zebrafish hypothalamic-pituitary axis.

The vertebrate hypothalamic-pituitary axis (HP) is the main link between the central nervous system and endocrine system. Although several signal pathways and regulatory genes have been implicated in adenohypophysis ontogenesis, little is known about hypothalamic-neurohypophysial development or when the HP matures and becomes functional. To identify markers of the HP, we constructed subtractive cDNA libraries between adult zebrafish hypothalamus and pituitary. We identified previously published genes, ESTs and novel zebrafish genes, some of which were predicted by genomic database analysis. We also analyzed expression patterns of these genes and found that several are expressed in the embryonic and larval hypothalamus, neurohypophysis, and/or adenohypophysis. Expression at these stages makes these genes useful markers to study HP maturation and function.

Caveolin-1 is required for lateral line neuromast and notochord development.

Caveolae have been linked to diverse cellular functions and to many disease states. In this study we have used zebrafish to examine the role of caveolin-1 and caveolae during early embryonic development. During development, expression is apparent in a number of tissues including Kupffer's vesicle, tailbud, intersomite boundaries, heart, branchial arches, pronephric ducts and periderm. Particularly strong expression is observed in the sensory organs of the lateral line, the neuromasts and in the notochord where it overlaps with expression of caveolin-3. Morpholino-mediated downregulation of Cav1alpha caused a dramatic inhibition of neuromast formation. Detailed ultrastructural analysis, including electron tomography of the notochord, revealed that the central regions of the notochord has the highest density of caveolae of any embryonic tissue comparable to the highest density observed in any vertebrate tissue. In addition, Cav1alpha downregulation caused disruption of the notochord, an effect that was enhanced further by Cav3 knockdown. These results indicate an essential role for caveolin and caveolae in this vital structural and signalling component of the embryo.

Zebrafish as a model for caveolin-associated muscle disease; caveolin-3 is required for myofibril organization and muscle cell patterning.

Caveolae are an abundant feature of many animal cells. However, the exact function of caveolae remains unclear. We have used the zebrafish, Danio rerio, as a system to understand caveolae function focusing on the muscle-specific caveolar protein, caveolin-3 (Cav3). We have identified caveolin-1 (alpha and beta), caveolin-2 and Cav3 in the zebrafish. Zebrafish Cav3 has 72% identity to human CAV3, and the amino acids altered in human muscle diseases are conserved in the zebrafish protein. During embryonic development, cav3 expression is apparent by early segmentation stages in the first differentiating muscle precursors, the adaxial cells and slightly later in the notochord. cav3 expression appears in the somites during mid-segmentation stages and then later in the pectoral fins and facial muscles. Cav3 and caveolae are located along the entire sarcolemma of late stage embryonic muscle fibers, whereas beta-dystroglycan is restricted to the muscle fiber ends. Down-regulation of Cav3 expression causes gross muscle abnormalities and uncoordinated movement. Ultrastructural analysis of isolated muscle fibers reveals defects in myoblast fusion and disorganized myofibril and membrane systems. Expression of the zebrafish equivalent to a human muscular dystrophy mutant, CAV3P104L, causes severe disruption of muscle differentiation. In addition, knockdown of Cav3 resulted in a dramatic up-regulation of eng1a expression resulting in an increase in the number of muscle pioneer-like cells adjacent to the notochord. These studies provide new insights into the role of Cav3 in muscle development and demonstrate its requirement for correct intracellular organization and myoblast fusion.