Urotensin II-related peptides, Urp1 and Urp2, control zebrafish spine morphology

  1. Elizabeth A Bearce
  2. Zoe H Irons
  3. Johnathan R O'Hara-Smith
  4. Colin J Kuhns
  5. Sophie I Fisher
  6. William E Crow
  7. Daniel T Grimes  Is a corresponding author
  1. Institute of Molecular Biology, Department of Biology, University of Oregon, United States

Abstract

The spine provides structure and support to the body, yet how it develops its characteristic morphology as the organism grows is little understood. This is underscored by the commonality of conditions in which the spine curves abnormally such as scoliosis, kyphosis, and lordosis. Understanding the origin of these spinal curves has been challenging in part due to the lack of appropriate animal models. Recently, zebrafish have emerged as promising tools with which to understand the origin of spinal curves. Using zebrafish, we demonstrate that the urotensin II-related peptides (URPs), Urp1 and Urp2, are essential for maintaining spine morphology. Urp1 and Urp2 are 10-amino acid cyclic peptides expressed by neurons lining the central canal of the spinal cord. Upon combined genetic loss of Urp1 and Urp2, adolescent-onset planar curves manifested in the caudal region of the spine. Highly similar curves were caused by mutation of Uts2r3, an URP receptor. Quantitative comparisons revealed that urotensin-associated curves were distinct from other zebrafish spinal curve mutants in curve position and direction. Last, we found that the Reissner fiber, a proteinaceous thread that sits in the central canal and has been implicated in the control of spine morphology, breaks down prior to curve formation in mutants with perturbed cilia motility but was unaffected by loss of Uts2r3. This suggests a Reissner fiber-independent mechanism of curvature in urotensin-deficient mutants. Overall, our results show that Urp1 and Urp2 control zebrafish spine morphology and establish new animal models of spine deformity.

Editor's evaluation

This is a beautifully executed study on the role of Urp signaling in spine morphogenesis in zebrafish. This work challenges the model that Urp1/ 2 controls the extension and straightening of the body axis of the zebrafish embryos, and thus makes a significant contribution to the literature.

https://doi.org/10.7554/eLife.83883.sa0

eLife digest

The backbone, or spine, is an integral part of the human body, providing support to our torsos so that we can sit, stand, bend and twist. If this structure does not form correctly, it can lead to pain, neurologic problems, and mobility issues. The spine normally has curves, but these can become deformed for many reasons, including genetic and muscular factors. There are also cases in which the cause of a spine distortion is unknown, such as in scoliosis (where the spine twists to the side), lordosis (where the lower part of the spine curves excessively), and kyphosis (where the upper part of the spine shows extreme curvature).

The structure of the spine is laid out during embryonic development and maintained throughout life. Experiments in zebrafish have shown that a crucial element in preserving the shape of the spine is the flow of cerebrospinal fluid or CSF. Propelled by the movement of little ‘hairs’ at the surface of specialized cells, this liquid runs through our central nervous system along a cavity lined with neurons. These nerve cells produce Urp1 and Urp2, two short molecules (or peptides) built from the same components as proteins. In zebrafish embryos, lowering the levels of these peptides had previously been shown to cause early body deformities. But what role, if any, do Urp1 and Urp2 play in maintaining the shape of the spine in adult zebrafish?

Bearce et al. set out to answer this question. First, they generated mutant zebrafish which did not carry either Urp1, Urp2 or both peptides. Contrary to previous findings, all three of these mutants developed normally as embryos. Once they were adults, zebrafish lacking Urp1 exhibited normal spines, while those lacking Urp2 had slightly deformed curves. However, zebrafish lacking both peptides had prominent curves in the tail-region of their spines, somewhat akin to lordosis in humans. This indicates that both peptides are necessary for adult spine structure, but work in a semi-redundant manner. Interestingly, the defects observed first appeared in adolescent fish and gradually worsened as they grew; many forms of human spinal abnormalities follow a similar trajectory.

Bearce et al. also tested the role of the protein Uts2r3, a receptor for peptides which belong to the urotensin family (such as Urp1 and Urp2). Fish lacking this protein showed normal spine structure as embryos, but distorted spinal curves as adults, suggesting that Urp1 and Urp2 might control spine morphology by signaling via the Uts2r3 receptor.

Together, Bearce et al.’s observations show that disturbing urotensin signaling leads to a lordosis-like condition in adult zebrafish, with evident deformities in the tail-region of the spine. Considering the broad similarities in structures between the zebrafish and the human spine, these results point to a possible involvement of urotensin signaling in spine distortion in humans. More studies using zebrafish will likely provide further insights into the principles that control the shape of the spine and what goes wrong when it breaks down.

Introduction

Understanding how the shape of organisms is acquired is a central goal of developmental biology. The chordate body axis forms during embryonic development, when it is based around the rod-like notochord (Stemple, 2005). Later, the vertebrate axis comprises a column of repeating vertebrae which grows during juvenile and adolescent phases and is then maintained during adult life for up to several decades in some species (Bagnat and Gray, 2020). While a great deal has been learned about how the body axis emerges during embryogenesis, less is known about how spine morphology is maintained during growth and adulthood.

A breakdown of spine morphology occurs in scoliosis, lordosis, and kyphosis. Scoliosis is medically defined as lateral curvatures of the spine greater than 10° (Cheng et al., 2015; Mesiti, 2021; Wise et al., 2008) and can be caused by congenital defects of vertebral patterning or as a secondary consequence of neuromuscular disease (Pourquié, 2011; Wishart and Kivlehan, 2021). However, most cases of scoliosis are idiopathic in nature, with no known etiology: approximately 3% of children are afflicted by idiopathic scoliosis, which most often onsets during adolescence (Cheng et al., 2015; Labrom et al., 2021). By contrast, kyphosis and lordosis occur when there is excessive curvature of the thoracic and lumbar regions of the vertebral column, respectively, resulting in a hunched upper back (kyphosis) or a concave lower back (lordosis) (Ogura et al., 2021) without vertebral structural defects. Since these categories of curves can co-occur, there are likely to be overlapping as well as distinct causes.

A challenge to understanding the origin of spinal curvature has been the dearth of suitable animal models recapitulating disease states. Recently, teleost fishes, in particular zebrafish (Danio rerio), have emerged as prominent animal models of spinal deformity (Bagnat and Gray, 2020; Bearce and Grimes, 2021; Boswell and Ciruna, 2017; Gorman and Breden, 2009; Roy, 2021). Using zebrafish, it was found that motile cilia-generated cerebrospinal fluid (CSF) flow is essential for maintaining body and spine morphology (Grimes et al., 2016). Mutants with defective motile cilia failed to undergo axial straightening during embryogenesis and so developed a misshapen early embryonic body axis called ‘curly tail down’ (CTD; Brand et al., 1996). If rescued during this early stage, mutants went on to develop three-dimensional spinal curves that recapitulated some features of idiopathic scoliosis, including adolescent-stage onset in the absence of vertebral patterning defects (Grimes et al., 2016; Marie-Hardy et al., 2021; Wang et al., 2022). Precisely how motile cilia and CSF flow maintain spine morphology during growth is not understood, but it is known that during early larval stages cilia motility is essential for the assembly of the Reissner fiber (RF), an extracellular thread-like structure composed predominantly of the large glycoprotein SCOspondin (encoded by sspo) which sits in the CSF in brain ventricles and the central canal (Cantaut-Belarif et al., 2018; Rodríguez et al., 1998). Zebrafish sspo mutants exhibited CTD as embryos while hypomorphic mutants which can survive beyond embryonic stages also manifested spinal curves (Cantaut-Belarif et al., 2018; Lu et al., 2020; Rose et al., 2020; Troutwine et al., 2020).

The URPs, Urp1 and Urp2, may also function downstream of motile cilia in the central canal. Urp1 and Urp2 are 10-amino acid cyclic peptides previously linked to heart disease and mental illness (Sugo et al., 2003; Konno et al., 2013; Nobata et al., 2011; Parmentier et al., 2011; Quan et al., 2021; Tostivint et al., 2006; Vaudry et al., 2010). In zebrafish, Urp1 and Urp2 are expressed in CSF-contacting neurons (CSF-cNs), flow sensory neurons in the central canal, and their expression is increased by motile cilia function and the RF (Cantaut-Belarif et al., 2020; Lu et al., 2020; Quan et al., 2015; Zhang et al., 2018). Morpholino knockdown of Urp1/Urp2 results in embryonic CTD phenotypes while addition of Urp1/Urp2 peptides can rescue the CTD of cilia motility- and RF-deficient mutants (Lu et al., 2020; Zhang et al., 2018). This suggested that Urp1 and Urp2 act downstream of cilia motility to promote early axial straightening (Grimes, 2019; Lu et al., 2020; Zhang et al., 2018).

Here, we set out to address whether Urp1 and Urp2 function beyond embryogenesis in maintaining body and spine morphology during growth and adulthood. By generating zebrafish mutants lacking Urp1 and Urp2 peptides, we found that they are essential, in a semi-redundant fashion, for adult spine morphology. Loss of Urp1 and Urp2 together led to the onset of spinal curves during adolescent stages and, by adulthood, resulted in planar curves in the caudal region of the spine that occurred without vertebral patterning defects or significant structural malformations. A similar phenotype was present upon mutation of the urotensin receptor (UT) gene, uts2r3, suggesting that Urp1 and Urp2 signal via Uts2r3 to maintain spine morphology. Urotensin-associated curves were quantitatively distinct from the curves displayed by cfap298 mutants, which lack cilia motility, and pkd2l1 mutants in which a CSF-cN-localized ion channel is mutated, suggestive of overlapping but distinct roles of these components. Moreover, RF breakdown preceded curve formation in cfap298 mutants while RF structure was maintained before and after curves appeared in uts2r3 mutants. Overall, this demonstrates that Urp1 and Urp2 peptides control the morphology of the zebrafish spine. We suggest that urotensin-deficient zebrafish model human spinal deformities and will be important tools for deciphering how the spine is maintained and how this process goes wrong in disease.

Results

Urp1 and Urp2 peptides are dispensable for embryonic axial straightening

To determine whether Urp1 and Urp2 are required for spine morphology, we used CRISPR/Cas9 to generate zebrafish mutant lines. Urp1 and Urp2 are encoded by 5-exon genes with the final exon coding for the 10-amino acid peptides that are released by cleavage from the pro-domain (Figure 1A–B, Figure 1—figure supplement 1A–B). We used pairs of guide RNAs to induce deletions across the genetic region coding for the peptides (Figure 1A, Figure 1—figure supplement 2A–B). We refer to the resulting mutant lines as urp1∆P and urp2∆P because they lack the peptide coding sequence. In addition, mRNA quantitation revealed downregulation of urp1 and urp2 in their respective mutant backgrounds, indicating transcript decay (Figure 1—figure supplement 2E).

Figure 1 with 3 supplements see all
Urp1 and Urp2 are dispensable for axial straightening.

(A) urp1 and urp2 are 5-exon genes (gray boxes). The final exon codes for the 10-amino acid peptides produced after cleavage from the prodomain at a dibasic site (KR). Pairs of gRNAs were used to induce deletions of Urp1 and Urp2 peptide coding sequences, resulting in urp1∆P and urp2∆P mutants, respectively. SP – signal peptide. (B) Urp1 and Urp2 peptide sequences with identical hexacyclic regions. (C) Zebrafish posterior axial straightening, the morphogenetic process which straightens the embryonic body. (D) Fluorescence in situ hybridization based on hybridization chain reaction analysis of pkd2l1, urp1, and urp2 expression in the central canal at 28 hpf. pkd2l1 expression marks CSF-cNs. urp1 expression is restricted to ventral CSF-cNs while urp2 is expressed in all CSF-cNs. Both urp1 and urp2 are expressed in cfap298tm304 and sspob1446 mutants, though comparison of expression between samples was non-quantitative. (i) Shows the zebrafish trunk with the yolk stalk labeled (*). (ii) Shows zoomed regions taken at the rostro-caudal level at the end of the yolk stalk. Scale bars: 150 µm (i), 10 µm (ii). (E) Lateral views of 28–30 hpf germline mutants (i) and crispants (ii). The urp1∆P;urp2∆P double mutants are maternal zygotic (MZ) mutants. Sibling (sib) and Cas9-only injected embryos served as controls. All embryos were incubated at 28°C, which is a restrictive temperature for cfap298tm304. (F) Quantitative reverse transcriptase PCR (qRT-PCR) analysis of urp1 and urp2 mRNA expression levels in cfap298tm304 and sspob1446 mutants at 28 hpf. n>3 biologically independent samples. Bars represent mean ± s.e.m. Two-tailed student’s t test used to calculate p-values. (G) Schematic of crispant generation and body curve analysis. (H) Quantitation of crispant body curves where bars represent mean ± s.d. for at least three independent clutches and injection mixes. The total number of embryos analyzed is given. *p<0.0001, student’s t test applied. UI – uninjected.

Figure 1—source data 1

Raw data for qRT-PCR and crispant body angle measurements.

https://cdn.elifesciences.org/articles/83883/elife-83883-fig1-data1-v2.xlsx

We first assessed urp1∆P and urp2∆P mutants for embryonic phenotypes. A previous morpholino-based knockdown study concluded that Urp1 and Urp2 are required for axial straightening, the process by which the ventrally curved zebrafish embryo straightens as the trunk elongates and detaches from the yolk (Figure 1C). Urp1/Urp2 morphants failed to undergo straightening and therefore displayed CTD (Zhang et al., 2018). Surprisingly, both urp1∆P and urp2∆P mutants underwent normal axial straightening and did not exhibit CTD (Figure 1Ei). By contrast, we observed CTD in both cfap298tm304 mutants that lack cilia motility in the central canal (Bearce et al., 2022) and sspob1446 mutants in which the RF constituent SCOspondin is mutated, as expected (Figure 1Ei, Figure 1—figure supplement 3A–B). Notably, cfap298tm304 and sspob1446 mutants maintained urp1 and urp2 expression in CSF-cNs, central canal neurons marked by pkd2l1 expression (Figure 1D). However, urp1 and urp2 transcripts were quantitatively reduced in cfap298tm304 and sspob1446 mutants (Figure 1F). We reasoned that the absence of CTD in urp1∆P and urp2∆P mutants might reflect redundancy, since Urp1 and Urp2 peptides are highly similar, with identical hexacyclic regions (Figure 1B, Figure 1—figure supplement 1A–B). Alternatively, maternally derived urp1 and/or urp2 transcripts may function to prevent phenotypes from manifesting. However, maternal zygotic urp1∆P;urp2∆P double mutants also exhibited linear body axes (Figure 1Ei), ruling out redundant or maternal gene product function. This demonstrates that Urp1 and Urp2 peptide-null mutants undergo axial straightening.

To confirm this finding, we performed additional Urp1 and Urp2 loss-of-function experiments. By injecting four guide RNAs (gRNAs) along with Cas9 into wild-type embryos at the one-cell stage, we generated mosaic mutants, called crispants, that were then assessed for body shape phenotypes (Figure 1G). In positive control experiments, cfap298 and sspo crispants exhibited robust CTD, phenocopying germline cfap298tm304 and sspob1446 mutants (Figure 1Eii). Quantitation of body curvature revealed that crispant generation was highly efficient, with CTD penetrance being close to 100% (Figure 1H). By contrast, urp1 and urp2 single and double crispants exhibited straight body axes that were not different to uninjected embryos or embryos injected with Cas9 only (Figure 1Eii and H). Using T7 endonuclease assays, we confirmed that high levels of insertion-deletion mutations were generated at gRNA sites in crispants (Figure 1—figure supplement 2C–D). We used the AB genetic background for the majority of our work, but we also generated and phenotyped urp1;urp2 double crispants on WIK and TU backgrounds to test for potential background effects. Normal axial straightening upon mutation of urp1 and urp2 was also observed on these backgrounds (Figure 1—figure supplement 2F). Overall, crispant results confirmed germline mutant findings. We conclude that Urp1 and Urp2 peptides are dispensable for axial straightening in embryonic zebrafish.

Urp1 and Urp2 function semi-redundantly to maintain spine morphology

Next, we determined the impact of Urp1 and Urp2 loss on adult spine morphology. Outwardly, urp1∆P mutant adults at 3 months post fertilization (mpf) appeared normal whereas urp2∆P mutants exhibited minor body dysmorphologies and kinked tails (n=72 for urp1∆P mutants and n=92 for urp1∆P mutants, Figure 2—figure supplement 1A). To assess spine morphology directly, we imaged bone by X-ray microcomputed tomography (µCT). Three-dimensional reconstitutions of µCT data from 3 mpf fish showed that urp1∆P mutants indeed exhibited overtly normal skeletal morphology (n=7) while urp2∆P mutants showed slight sagittal curves (n=4; Figure 2—figure supplement 2A, Figure 2—videos 1–3). By contrast to these absent or mild deformities in single mutants, urp1∆P;urp2∆P double mutants exhibited prominent curves, with significant dorsal-ventral Cobb angles, a measure of deviation from straightness, especially in the caudal region of the spine (Figure 2B and D–F, Figure 2—figure supplements 1A and 2A, Figure 2—videos 1–4). These data indicate that Urp1 and Urp2 are essential for adult spine morphology, and that they function in a semi-redundant fashion in this context.

Figure 2 with 11 supplements see all
Urp1 and Urp2 are required for proper adult spine morphology.

(A–C) Lateral views of microcomputed tomography reconstitutions of wild-type (A), urp1∆P;urp2∆P (B) and uts2r3b1436 (C) mutants at 3 mpf. (D) Cobb angle measurements for individual fish in the sagittal plane for urp1∆P;urp2∆P and uts2r3b1436 mutants. Circles represent angles for individual curves. (E-E’) Total Cobb angles with each circle representing an individual fish. The mean ± s.d. is shown. (G’) is the data from G parsed for sex. p-Values are given from two-tailed unpaired student’s t tests. (F) The position of curve apex is plotted and shows that most curves are in caudal vertebrae. n=9 and 8 for urp1∆P;urp2∆P and uts2r3b1436 mutants, respectively.

Figure 2—source data 1

Raw data from spinal curve phenotypic measurements.

https://cdn.elifesciences.org/articles/83883/elife-83883-fig2-data1-v2.xlsx

To assess the long-term maintenance of spine morphology in Urp1- and Urp2-deficient conditions, we aged urp1∆P and urp2∆P single mutants to 12 mpf then performed µCT. At this later time point, urp1∆P and urp2∆P mutants exhibited mild kyphosis-like curves though urp2∆P mutants were more severe (Figure 2—figure supplement 2B, Figure 2—videos 5–7). These degenerative phenotypes demonstrate that Urp1 and Urp2 are essential for maintenance of spine morphology throughout adulthood and aging, and suggest that Urp2 plays a larger role than Urp1.

Urp1 and Urp2 signal through the Uts2r3 receptor to control spine morphology

Urp1 and Urp2 peptides engage G-protein-coupled receptors (Ames et al., 1999; Chatenet et al., 2004; Elshourbagy et al., 2002; Labarrère et al., 2003; Liu et al., 1999; Nothacker et al., 1999). While a single urotensin II receptor (UT) gene is found in humans and has recently been linked to abnormal spinal curvature (Dai et al., 2021), the zebrafish genome encodes five such receptors. One of those, Uts2r3, was previously implicated in spine morphology (Zhang et al., 2018). To systematically compare the effects of Uts2r3 receptor mutation with loss of Urp1 and Urp2, we generated a uts2r3 mutant line harboring a 178-amino acid deletion after the third amino acid, significantly disrupting the protein (Figure 2—figure supplement 3). Like urp1∆P;urp2∆P double mutants, these uts2r3b1436 mutants underwent normal axial straightening as embryos (Figure 1Ei) and went on to exhibit spinal curves as adults (Figure 2C, Figure 2—figure supplement 1A, Figure 2—video 8). Cobb angle measurements showed that uts2r3b1436 mutants and urp1∆P;urp2∆P mutants were similar, though curves in urp1∆P;urp2∆P mutants were slightly more severe (Figure 2D–E). Like urp1∆P;urp2∆P mutants, uts2r3b1436 mutants showed mostly caudally located curves, especially in the most rostral of the caudal vertebrae (Figure 2F). Thus, although we cannot rule out minor roles for other UT receptors, these data suggest that Urp1 and Urp2 control spine morphology largely by signaling through Uts2r3.

Urotensin pathway mutants display adolescent-onset spinal curves in the absence of structural vertebral defects

Next, we determined whether urotensin pathway mutants recapitulated any signs of disease present in patients. Several types of human spinal curves onset during adolescent growth (Cheng et al., 2015). To discern the stage of onset of curves in urotensin pathway mutants, we monitored urp1∆P;urp2∆P double mutant cohorts as they grew. Subtle curves first became apparent between 9 and 11 days pf (dpf), corresponding to a standard length between 3.9±0.7 mm and 5.9±0.4 mm (mean ± s.d.; Figure 3A–B, Figure 3—video 1), a stage when adolescents were rapidly growing (Figure 3C). By 13 dpf (standard length 6.2±0.3 mm), curves were evident in all urp1∆P;urp2∆P mutants and progressively worsened up to 17 dpf (8.3±0.4 mm) when we ended this analysis (Figure 3A–B). At 1 mpf, we assessed urp1∆P;urp2∆P mutants by µCT and found variability in curve position and amplitude (Figure 3D). Notably, at this stage, five of seven mutants exhibited a significant curve in the pre-caudal vertebrae, in addition to a caudal curve (Figure 3D, Figure 3—figure supplement 1Ai-Bii). Since pre-caudal curves were rare in mutants at 3 mpf (Figure 2B and F), this suggested that curve location is dynamic and that pre-caudal curves form then resolve or shift in some mutants as they grow to adulthood.

Figure 3 with 3 supplements see all
urp1∆P;urp2∆P mutants exhibit adolescent-onset spinal curves without significant structural vertebral defects.

(A) Lateral views of control fish and age-matched urp1∆P;urp2∆P mutants. Arrows point to forming body curves. (B) Traces of body shape every 2 days for 5 fish per time point from 3 to 17 dpf. (C) Growth curves for control and urp1∆P;urp2∆P mutants were indistinguishable. Arrow shows time of curve onset. Mean ± s.d. is plotted. (D) Microcomputed tomography (µCT) reconstitutions of spines at 1 mpf with heads, fins, and ribs digitally removed. Scale bar: 1 mm. (E) µCT reconstitutions of three pre-caudal and two caudal vertebrae including frontal and lateral views of the highlighted vertebra. No major structural defects such as fusions were observed in urp1∆P;urp2∆P mutants. (F) Vertebral body rostral-caudal length (L1/L2) and dorsal-ventral height (H1/H2) aspect ratios for six vertebrae of wild type (n=4) and urp1∆P;urp2∆P mutants (n=4). Length aspect ratios were significantly more variable in mutants, but height aspect ratios were unchanged (p=0.022 and 0.745, respectively, Bartlett’s test for equal variances). (G) Calcein staining revealed well-structured vertebrae forming in control (standard length 5.7 mm) and urp1∆P;urp2∆P mutant (standard length 5.7 mm) fish. n>30 fish per condition.

Figure 3—source data 1

Raw data from larval growth measurements and vertebral quantitation.

https://cdn.elifesciences.org/articles/83883/elife-83883-fig3-data1-v2.xlsx

Next, we assessed whether spinal curves in urotensin pathway mutants were caused by congenital defects of vertebral patterning or structure. Staining of juveniles with the vital dye calcein revealed no defects in vertebral patterning or spacing in urp1∆P;urp2∆P mutants at 10 dpf (4.5–6 mm standard length [Figure 3C]) (Figure 3G). We then assessed µCT data from 3 mpf fish and quantified vertebral body shape for vertebrae 11–16, where curves occurred in urp1∆P;urp2∆P mutants (Figure 2F). Vertebral body length and height aspect ratios were 1.00±0.007 (mean ± s.d.) and 0.99±0.007, respectively, for wild type (n=4; Figure 3F). By contrast, urp1∆P;urp2∆P mutants exhibited more variable vertebral length aspect ratios (0.97±0.03, p=0.022; Bartlett’s test for equal variances, Figure 3F). Vertebral height aspect ratios in mutants (1.00±0.006) were not significantly different to controls (p=0.745, Figure 3F). These data are consistent with subtle vertebral shape defects at the points of curvature in urp1∆Purp2∆P mutants. However, we do not observe vertebral fusions or missing or transformed appendages, suggesting that vertebral defects do not underlie spinal curves; instead, small changes in vertebral shape are likely due to the presence of curves themselves rather than causative of curves in the first place. This is similar to what occurs in non-congenital forms of human scoliosis (Cheng et al., 2015).

Additionally, we parsed our phenotypic data for sex since spinal curves often show sex bias in severity in humans (Cheng et al., 2015), something which has been recapitulated in some zebrafish spinal curve mutants (Marie-Hardy et al., 2021). However, in both urp1∆P;urp2∆P and uts2r3b1436 mutants, we found no significant differences in curve penetrance or severity between males and females (Figure 2E’).

Urotensin-deficient mutants are phenotypically distinct from cfap298tm304 and pkd2l1icm02 mutants

We next compared the phenotypes of urotensin pathway mutants to other mutant lines that exhibit spinal curves. The cfap298tm304 line harbors a temperature-sensitive mutation in cfap298, a gene required for cilia motility in several organisms including humans (Austin-Tse et al., 2013; Bearce et al., 2022; Jaffe et al., 2016). cfap298tm304 mutants exhibit reduced cilia motility in the central canal and, if the resulting CTD is embryonically rescued by temperature shifts, develop adolescent-onset spinal curves (Grimes et al., 2016). These curves were argued to model an adolescent idiopathic scoliosis (AIS)-like condition (Figure 4A–B, Figure 4—video 1; Grimes et al., 2016; Marie-Hardy et al., 2021). Importantly, both urp1 and urp2 transcripts were significantly downregulated in cfap298tm304 mutants (Figure 1F) suggesting that spinal curves in cfap298tm304 might be the result of reduced Urp1/Urp2 expression. To systematically compare cfap298tm304 mutants with urotensin-deficient mutants, we raised cohorts of urp1∆P;urp2∆P mutants, uts2r3b1436 mutants, and temperature-shift-rescued cfap298tm304 mutants alongside one another in the same aquatics facility after backcrossing all lines to the AB strain for multiple generations. At 3 mpf, we performed µCT scanning and three-dimensional reconstitutions. First, we calculated dorso-ventral Cobb angles, which revealed that cfap298tm304 mutants were more severely curved (average total Cobb angle: 260.4±32.7°) than either urp1∆P;urp2∆P mutants (197.5±38.9°) or uts2r3b1436 mutants (146.1±47.4°) (Figure 4E, Figure 2E, Figure 2—figure supplement 1D–E). Second, cfap298tm304 mutants showed prominent dorsal-ventral curves in the pre-caudal as well as caudal vertebrae, a distinct pattern compared with the predominantly caudal curves in urp1∆P;urp2∆P and uts2r3b1436 mutants (Figure 4A–C, Figure 2—figure supplement 1C). Third, cfap298tm304 mutants exhibited significant lateral curvature of the spine, often with spinal twisting, a hallmark of AIS-like curves (Figure 4Fi–Fii, Figure 4—figure supplement 1). By contrast, urp1∆P;urp2∆P and uts2r3b1436 mutants showed planar curves, with very minor or no lateral deviations (Figure 4Fi–Fii, Figure 4—figure supplement 1). These results demonstrated that cilia motility mutants and urotensin-deficient mutants exhibit distinct spinal curve phenotypes. As such, the causes of spinal curves in cfap298tm304 mutants can only be partially explained by reduced Urp1/Urp2 expression.

Figure 4 with 5 supplements see all
urp1∆P;urp2∆P mutants and cfap298tm304 mutants are phenotypically distinct.

(A–D) Lateral views of microcomputed tomography (µCT) reconstitutions of wild-type (A), cfap298tm304 mutants (B), urp1∆P;urp2∆P double mutants (C), and cfap298tm304;urp1∆P;urp2∆P triple mutants (D). All fish shown are female. Scale bar: 10 mm. (E) Total Cobb angles with each circle representing an individual fish. The mean ± s.d. is shown. p-Values are given from two-tailed unpaired student’s t tests. U — urp1∆P;urp2∆P double mutants; C — cfap298tm304 mutants; UC — cfap298tm304;urp1∆P;urp2∆P triple mutants. (Fi) Dorsal views of µCT reconstitutions with ribs and fins removed. Scale bar: 5 mm. (Fii) Quantitation of degree of lateral curvature for wild type (n=5) and urp1∆P;urp2∆P (n=8), uts2r3b1436 (n=3), cfap298tm304 (n=3), and cfap298tm304;urp1∆Purp2∆P (n=6) mutants. y-axis is the arbitrary units.

To further explore the relationship between cilia motility and urotensin peptides, we generated adult cfap298tm304;urp1∆P;urp2∆P triple mutants that were embryonically rescued by temperature shifts. Triple mutants exhibited significant curves, similar to cfap298tm304 single mutants (Figure 4D–E, Figure 2—figure supplement 1C–E, Figure 4—video 2; average total Cobb angle: 284.4±33.5°). The curves of triple mutants were three-dimensional in nature, with both dorso-ventral and medio-lateral deviations, as well as incidences of spinal torsion (Figure 4Fi–Fii, Figure 4—figure supplement 1); curves were present in both pre-caudal and caudal vertebrae, as opposed to being more restricted to caudal vertebrae in urp1∆P;urp2∆P double mutants (Figure 2—figure supplement 1C). Overall, this suggests that motile cilia contribute to urotensin-dependent and urotensin-independent pathways controlling spine morphology.

Pkd2l1 is a polycystin family ion channel expressed in CSF-cNs, the same cell type which expresses Urp1 and Urp2 (Figure 1D; Quan et al., 2015). Pkd2l1 is responsible for flow-induced Ca2+ signaling in CSF-cNs (Böhm et al., 2016; Sternberg et al., 2018). While pkd2l1icm02 mutants exhibited normal early axis development, they went on to develop mild kyphosis-like curves upon aging (Sternberg et al., 2018), a result we recapitulated after raising pkd2l1icm02 mutants on the same genetic background and under the same conditions as urotensin-deficient mutants for direct comparison (Figure 2—figure supplement 2A–B). Notably, pkd2l1icm02 mutants showed very subtle curves in the pre-caudal vertebrae and absence of lateral deviation at 12 mpf but no obvious curves of any kind at 3 mpf (Figure 2—figure supplement 2A-B, Figure 4—videos 3–4). This mild kyphosis-like phenotype was therefore also highly distinct from urp1∆P;urp2∆P and uts2r3b1436 mutants.

Overall, our phenotypic data show that urp1∆P;urp2∆P mutants develop adolescent-onset curves without significant vertebral structural defects or sex bias. Coupled to the consistently caudal location of curves at 3 mpf as well as the lack of lateral deviation, we suggest that urotensin pathway mutants most closely reflect a lordosis-like condition. In agreement, urotensin mutants were phenotypically distinct from established AIS-like (cfap298tm304) and kyphosis-like (pkd2l1icm02) models.

RF breakdown precedes AIS-like curves in cfap298tm304 mutants

Given the links between motile cilia, the RF and urp1 and urp2 expression (Figure 1D and F; Cantaut-Belarif et al., 2020; Lu et al., 2020; Zhang et al., 2018) as well as the requirement for proper RF assembly to prevent spinal curves (Rose et al., 2020; Troutwine et al., 2020), we assessed Sspo, the major component of the RF, in spinal curve mutants. To visualize Sspo localization, we used the sspo-GFPut24 line in which GFP coding sequence is fused to the endogenous sspo locus, producing Sspo-GFP protein (Troutwine et al., 2020). First, we assessed Sspo localization in the central canal of cfap298tm304 mutants at 28 hr pf (hpf). In sibling controls, Sspo localized into an RF throughout the central canal (Figure 5A). By contrast, cfap298tm304 mutants raised at restrictive temperatures, which exhibited reduced central canal cilia motility and CTD (Figure 1E; Bearce et al., 2022), lacked RF. Instead, Sspo was diffusely localized in the central canal in these mutants (Figure 5B), in agreement with previous work showing that cilia motility is required for RF assembly in embryos (Cantaut-Belarif et al., 2018).

Reissner fiber (RF) breakdown in cfap298tm304 mutants but not urotensin-deficient mutants.

(A–F) Grayscale maximal intensity projection of Sspo-GFP localization in the central canal in 28 hpf embryos (A–C) and 12 dpf adolescents (D–F). RF is denoted by arrow heads in D and F. Arrows point to structures along the central canal that become GFP-positive in cfap298tm304 and uts2r3b1436 mutants. Scale bar: 10 µm. (G) Schematic of temperature shift experiment in which cfap298tm304 mutants are initially raised at permissive temperatures before being shifted to restrictive temperatures at 6 dpf, then imaged at 12 dpf. (H–I) Lateral views of cfap298tm304 (H) and uts2r3b1436 (I) mutants at 12 dpf when Sspo-GFP imaging took place. The white box in H shows the location imaged in D–F.

Next, we took advantage of the temperature-sensitive nature of the cfap298tm304 mutation to determine whether RF is also disrupted at later timepoints, during adolescent stages when true spinal curves begin to develop. To do so, we initially raised cfap298tm304 mutants at permissive temperatures, allowing RF to correctly form and embryos to fully straighten. At 6 dpf, we transitioned larvae to restrictive temperatures (Figure 5G), which led to the appearance of curves by 11–14 dpf (standard length 6.5±0.2 mm), then assessed Sspo localization at 12 dpf. As in embryos, Sspo localized into a defined RF in the central canal in sibling controls (arrow head in Figure 5D) but was diffuse in temperature upshifted cfap298tm304 mutants, both in mutants that had yet to develop curves and those exhibiting subtle curves (Figure 5E, Figure 5H). This demonstrates (1) that cilia motility is not only required for the initial formation of the RF but also for its maintenance; and (2) breakdown of the RF precedes curve onset in cilia motility-deficient mutants. This supports a model in which continued cilia motility maintains the RF structure, and loss of the RF causes AIS-like curves to develop in cilia motility mutants.

The RF remains intact in Uts2r3-deficient mutants both before and after curve formation

Next, we imaged Sspo-GFP localization in the central canal of uts2r3b1436 mutants to determine whether RF breakdown could be causative in urotensin-associated spinal curves. In agreement with that lack of axial straightening phenotype in uts2r3b1436 mutant embryos (Figure 1Ei), we found that Sspo localized normally into the RF in the central canal at 28 hpf in uts2r3b1436 mutants (Figure 5C). At 12 dpf, as spinal curves were beginning to manifest (Figure 5I), Sspo still formed an intact RF in uts2r3b1436 mutants (arrow head in Figure 5F). The fact that RF is present in uts2r3b1436 mutants as curves form suggests that urotensin-associated curves are not caused by defective RF formation. This result coheres with the distinct spinal phenotypes exhibited by cfap298tm304 mutants and uts2r3b1436 mutants.

While imaging Sspo-GFP at 12 dpf, we noted that in both cfap298tm304 mutants and uts2r3b1436 mutants, large central canal cells with the appearance of CSF-cNs became GFP-positive (arrows in Figure 5E–F), something we rarely observed in control fish (Figure 5D). Indeed, the GFP-positive central canal cells in uts2r3b1436 give the effect of making the RF appear comparatively smaller/dimmer (compare RF in Figure 5F and Figure 5D). We suggest that CSF-cNs may endocytose Sspo-GFP monomers in cfap298tm304 mutants where the RF has broken down and in uts2r3b1436 mutants where RF is likely to be making increasing numbers of contacts with CSF-cNs (Orts-Del’Immagine et al., 2020) owing to the onset of spinal curves.

Discussion

Urotensin II (UII) is a cyclic peptide that was first identified from the teleost urophysis (Pearson et al., 1980) and subsequently found to exist in amphibians (Conlon et al., 1992) and mammals (Coulouarn et al., 1998). A highly similar peptide, called URP was then isolated from the brains of rodents (Sugo et al., 2003). UII and URP both signal via the UT, a G-protein-coupled receptor (Ames et al., 1999; Liu et al., 1999; Mori et al., 1999; Nothacker et al., 1999). UIIs, URPs, and UTs have been linked to cardiovascular function and inflammation, but their roles in the development of morphology are little understood.

In this study, we discovered a role for two of the URPs, Urp1 and Urp2, in zebrafish spine morphology. To do so, we generated urp1∆P and urp2∆P mutants that lacked the genetic region coding for the Urp1 and Urp2 peptides, respectively, then phenotyped skeletal morphology by µCT. This revealed that Urp1 and Urp2 function semi-redundantly to control spine morphology, with double mutants, as well as Uts2r3 receptor mutants, developing spinal curves in the caudal region during adolescent growth. By contrast, single urp1∆P and urp2∆P mutants developed more subtle curves that nevertheless worsened with age. The lack of major vertebral defects, the location and direction of the curves coupled with phenotypic differences compared with mutants that model AIS and kyphosis, suggested that urotensin-deficient mutants model a lordosis-like condition.

The RF, a long proteinaceous thread-like structure which sits in the central canal and is mostly made from SCOspondin (encoded by sspo), has been implicated in controlling body axis and spine morphology (Cantaut-Belarif et al., 2018; Lu et al., 2020; Rose et al., 2020; Troutwine et al., 2020). We find that RF breaks down prior to curve formation in the cilia motility cfap298tm304 mutant. Given other studies linking presence of the RF with a linear body axis, this strongly suggests that RF breakdown is a major factor driving spinal curves in cfap298tm304 mutants. By contrast, Uts2r3-deficient mutants exhibited an intact RF, demonstrating that curves are not formed by RF breakdown in urotensin pathway mutants and nor do the presence of curves significantly disrupt RF structure. This coheres with a model in which urotensin signals act downstream of RF function in controlling spine morphology. Similarly, urp1 and urp2 expression are known to be controlled by RF function during embryonic phases (Figure 1F; Cantaut-Belarif et al., 2020; Lu et al., 2020; Rose et al., 2020; Zhang et al., 2018).

Intriguingly, motile cilia mutants and sspo mutants exhibit three-dimensional spinal curves, with dorso-ventral and medio-lateral curvature (Grimes et al., 2016; Lu et al., 2020; Rose et al., 2020; Troutwine et al., 2020). By contrast, we find that urotensin-deficient mutants exhibit largely planar curves, only in the dorso-ventral direction, thereby potentially uncoupling two systems controlling posture (Picton et al., 2021). Urp1 and Urp2 are expressed in CSF-cNs (Figure 1D; Quan et al., 2015), a cell type which consists of both dorsal and ventral subpopulations. While Urp1 and Urp2 are co-expressed in ventral CSF-cNs, only Urp2 is expressed in dorsal CSF-cNs (Figure 1D; Quan et al., 2015). It is therefore tempting to speculate that dorso-ventral spine shape is mediated specifically by ventral CSF-cNs which express higher amounts of Urp1/Urp2 peptides, resulting in dorso-ventral curves in urotensin-deficient mutants. This also suggests, as above, that while Urp1/Urp2 expression is in part controlled by upstream cilia motility and RF function, decreased urotensin signaling cannot account fully for the spinal curve phenotypes that occur upon loss of cilia motility or the RF. Indeed, the RF is required for both dorsal and ventral CSF-cN function (Orts-Del’Immagine et al., 2020), which may explain why dorso-ventral and medio-lateral curves occur when RF is disrupted either by cilia motility mutations or mutations to SCOspondin. This scenario is further complicated by the finding that increased urp1 and urp2 expression occurs upon mutation of rpgrip1l, a gene encoding a component of the ciliary transition zone (Vesque et al., 2019 — preprint). As such, various cilia-dependent signals likely control the precise levels of Urp1/Urp2 peptides, and controlling those levels appears critical for maintaining the shape of the spine.

In addition to exhibiting phenotypic differences in terms of curve direction compared with cilia motility and RF mutants, urotensin-deficient mutants also showed gradual worsening of spinal phenotypes upon aging. urp1∆P mutants exhibited no obvious phenotypes at 3 mpf but by 12 mpf showed abnormal curves, while urp2∆P showed mild curves at 3 mpf and more severe deformity at 12 mpf. This implies that Urp1/Urp2 function throughout adulthood and aging to maintain spine morphology. Since double mutants were more severe at 3 mpf than either single mutant, partial redundancy between Urp1 and Urp2 appears to occur while overall dose likely sets how early phenotypes manifest. In contrast to this long-term role for urotensin peptides, temperature-shift experiments in which motile cilia were inactivated after 34 dpf showed no role for motile cilia beyond this stage in maintaining the spine (Grimes et al., 2016). It is worth noting, however, that this result does not preclude a role for CSF flow or the RF during long-term spine morphostasis because it has not been determined if motile cilia are required for CSF flow or RF formation in adult fish.

While the study of zebrafish spinal curve mutants holds great promise for understanding the basic science of spine morphology, it will also be important for the field to grapple with the question of how closely zebrafish spinal deformity mutants truly recapitulate human spinal curve diseases. The spines of humans and zebrafish are broadly similar, and it has been suggested that spinal loads are comparable (Gorman and Breden, 2009). Moreover, zebrafish spines, like human spines, seem predisposed to curvature, with high levels of scoliosis-like curves naturally developing with age (Bearce and Grimes, 2021; Gorman and Breden, 2009). The overall shape of the zebrafish spine is also similar to human, with a natural kyphotic curve in the pre-caudal (rib-bearing) vertebrae and a compensatory, albeit very minor, lordotic curve in the most anterior caudal vertebrae. However, these curves are not as pronounced as in humans. Moreover, zebrafish also exhibit some fish-specific structures such as the Weberian apparatus (Dietrich et al., 2021). Nevertheless, zebrafish cilia motility mutants appear to model several features of AIS including the three-dimensional nature of curves, lack of vertebral patterning defects or significant vertebral structural malformations, adolescent-onset and, in some cases, sex bias (Grimes et al., 2016; Marie-Hardy et al., 2021; Wang et al., 2022). The curves of urp1∆P;urp2∆P mutants display some of these features as well but, importantly, are not three-dimensional. Instead, urotensin pathway-deficient mutants display primarily planar curves, with little or no lateral deviation. This is more similar to what occurs in hyper-kyphosis and hyper-lordosis in humans, when natural curves are accentuated. Given this planarity, and since curves are mostly present in the caudal vertebrae, we suggest that urp1∆P;urp2∆P mutants model some aspects of lordosis and so refer to this phenotype as lordosis-like. However, we note that human lumbar vertebrae and zebrafish caudal vertebrae are structurally distinct (Boswell and Ciruna, 2017), and humans have a significant natural lordotic curve that allows for an efficient upright walking gait, whereas zebrafish do not. Thus, urotensin-deficient mutants recapitulate some aspects of lordosis but clearly cannot mimic human-specific aspects of hyper-lordotic curves.

A surprising finding from our work was that Urp1 and Urp2 peptides are genetically dispensable for embryonic axial straightening. This interpretation is challenged by morpholino knockdown of Urp1/Urp2, which does result in failure of axial straightening in some individuals, resulting in a CTD phenotype (Zhang et al., 2018). One possibility is that the CTD of morphants results from morpholino off-target effects. However, this seems unlikely for three reasons: (1) adding exogenous Urp1/Urp2 peptides to the central canal can rescue the CTD phenotype of cilia motility- and RF-deficient mutants (Lu et al., 2020; Zhang et al., 2018), suggesting the involvement of Urp1/Urp2 in axial straightening, at least in gain-of-function experiments; (2) since motile cilia and the RF are required for both axial straightening during embryogenesis and for the maintenance of spine morphology during adolescence, it seems parsimonious that Urp1/Urp2 peptides would also function across these two life stages; and (3) urp1 and urp2 transcript levels are reduced in motile cilia and RF mutants (Figure 1F; Cantaut-Belarif et al., 2020; Lu et al., 2020; Zhang et al., 2018), suggestive of a link between Urp1 and Urp2 upregulation and axial straightening.

Nevertheless, the lack of CTD in our mutants, in which the Urp1 and Urp2 peptide coding sequences were entirely removed, is clear: single, double, and maternal-zygotic urp1∆P and urp2∆P mutants all underwent normal axial straightening. This strongly argues that Urp1 and Urp2 are dispensable for straightening. Since the deletions were induced toward the end of the protein, it does leave open the possibility that the pro-domain sequences are required for straightening. However, this seems unlikely for three reasons: (1) urp1∆P and urp2∆P mutants showed urp1 and urp2 transcript downregulation, respectively, in addition to deletion of the peptide coding regions; (2) urp1 and urp2 single and double crispants, in which gRNAs targeted several regions of the gene, also showed normal straightening; and (3) exogenous addition of Urp1 and Urp2 peptides, without pro-domains, rescued CTD of a cilia motility and RF mutant (Lu et al., 2020; Zhang et al., 2018), suggesting that it is the peptide itself and not some other region which is functional.

One possibility is that genetic compensation explains the mutant/morphant phenotypic differences (Rossi et al., 2015; Sztal and Stainier, 2020). In this putative scenario, a feedback response in the cell buffers otherwise harmful mutations, preventing their effects from manifesting phenotypically. A recently discovered compensation mechanism is transcriptional adaptation in which mutant mRNA is decayed and the products of that decay are, via a sequence-dependent mechanism, recruited to genes with similar sequences where they promote transcriptional upregulation (El-Brolosy et al., 2019; Ma et al., 2019). The upregulation of adapting genes then masks phenotypes in mutants but not morphants. We did find slight upregulation of urp2 in urp1∆P mutants (Figure 1—figure supplement 2E), which may indicate transcriptional adaptation, although this alone cannot explain the lack of CTD phenotypes in urp1∆P mutants because urp1∆P;urp2∆P double mutants also lacked CTD. It will be informative in the future to determine whether other urotensin II-encoding peptides (Figure 1—figure supplement 1A–B) are able to compensate, during embryonic phases, for loss of urp1 and urp2 or if other factors explain the mutant/morphant discrepancies.

While we were in the final stages of preparing this manuscript, a study was released which made several complementary findings, also concluding that Urp1 and Urp2 function redundantly to maintain spine shape (Gaillard et al., 2022 — preprint). Importantly, and in agreement with our work, Gaillard and colleagues observed no embryonic axial defects upon genetic loss of Urp1 and Urp2. Moreover, they found that other urotensin II-encoding genes were not upregulated in urp1 and urp2 mutants, arguing against phenotypic masking by genetic compensation. Our findings and those of Gaillard and colleagues together therefore strongly argue that Urp1 and Urp2 are not essential for axial straightening during embryogenesis but are instead required for the maintenance of the body axis during growth and adulthood. As such, other mechanisms, currently unknown, likely operate downstream of cilia motility and RF function to mediate embryonic axial straightening.

Future efforts will be required to discern which tissues respond to Urp1/Urp2 signals during the control of spine morphology. During embryonic stages, Uts2r3 is expressed in dorsal muscle (Zhang et al., 2018), but it remains unclear how Urp1/Urp2 peptides released by CSF-cNs could signal to effect muscle during adolescent stages. Moreover, based on single-cell RNA-sequencing gene expression atlases (Farnsworth et al., 2020), uts2r3 is expressed in several other cell types in addition to muscle. Tissue-specific ablations and rescue experiments should be used to untangle precisely where Uts2r3-dependent Urp1/Urp2 signaling occurs to control spine morphology. It will also be critical to determine the timing of action of urotensin signaling. While we observe phenotypes first appearing between 9 and 11 dpf in urp1∆P;urp2∆P mutants, it is possible that the underlying defect in mutants is caused by an earlier event which only phenotypically manifests later. One candidate is subtle disruptions to the notochord, which may later result in spinal curves (Bagwell et al., 2020). On the other hand, our experiments in which we aged urp1∆P and urp2∆P single mutants argue that Urp1 and Urp2 peptides function throughout life, and not only during early growth, to maintain spine morphology.

Another major question is whether the function of urotensin signaling in spine morphology is conserved in other species, and whether our findings are directly relevant to humans. These questions will require future work to answer, but two recent studies shed some light on these matters. First, in the frog Xenopus laevis, disruption of Utr4, a counterpart of Uts2r3, causes abnormal curvature of the body axis (Alejevski et al., 2021). Second, a human genetics study reported that rare mutations in UTS2R are significantly associated with spinal curvature, being discovered within AIS patient cohorts (Dai et al., 2021). Thus, a deeper understanding of the role of urotensin signaling in maintaining spine shape will not only provide insight into principles of morphogenesis but potentially also human disease.

Materials and methods

Zebrafish

AB, TU, and WIK strains of D. rerio were used. Zebrafish lines generated were sspob1446, urp1b1420 (called urp1∆P), urp2b1421 (called urp2∆P), uts2r3b1436 as well as previously published lines: cfap298tm304 (Jaffe et al., 2016), pkd2l1icm02 (Sternberg et al., 2018), and sspo-gfput24 (Troutwine et al., 2020). Experiments were undertaken in accordance with research guidelines of the International Association for Assessment and Accreditation of Laboratory Animal Care and approved by the University of Oregon Institutional Animal Care and Use Committee (Protocol number 21–45). All zebrafish strains and other materials are available upon request.

Generation and genotyping of mutant lines sspob1446, urp1b1420, urp2b1421, and uts2r3b1436 were generated using CRISPR/Cas9. gRNA oligos were designed using CRISPRscan (Moreno-Mateos et al., 2015). gRNA templates (IDT) were assembled by annealing and extension with Bottom strand ultramer_1 (Key resource table) using Taq Polymerase (NEB, M0273) with cycling parameters of 95°C (3 min), 95°C (30 s), 45°C (30 s), 72°C (30 s), and 72oC (10 min) with 30 cycles of the middle three steps. PCR product was purified (Zymo DNA Clean and Concentrator Kit, D4013) and then used for in vitro RNA synthesis using a MEGAshortscript T7 Transcription Kit (ThermoFisher, AM1354). Synthesized gRNAs were purified (Zymo RNA Clean and Concentrator Kit, R1013), then 150 pg along with 320 pg/nl Cas9 (IDT, 1081058) were injected into one-cell stage fertilized eggs. The mosaic mutant fish resulting from these injections (F0 fish) were raised and outcrossed to AB wild-types, and DNA was extracted from the resulting F1 embryos. Mutant alleles were screened by PCR coupled with restriction enzyme digestion and/or Sanger sequencing (GeneWiz). Embryos from F1 clutches harboring mutations were raised to adulthood and outcrossed to AB wild-types to generate F2 families which were screened for mutations and raised. The nature of mutations was identified by sequencing DNA extracted from adult fin clips of F2 heterozygous fish using CRISP-ID to deconvolute (Dehairs et al., 2016) and confirmed by sequencing DNA of F3 homozygous embryos.

urp1b1420 mutants contain a 279 bp deletion and 1 bp insertion that were genotyped by PCR amplification with urp1_geno_1 and urp1_geno_2 primers which generate a 460 bp band from wild-type DNA and a 184 bp band from mutant DNA. urp2b1421 mutants harbor a 61 bp deletion and were genotyped by PCR amplification with urp2_geno_1 and urp2_geno_2 primers followed by gel electrophoresis to distinguish the 283 bp wild-type band and the 226 bp mutant band. uts2r3b1436 mutants contain a 534 bp deletion and were also genotyped by PCR, using uts2r3_geno_1 and uts2r3_geno_2, in which wild-type sequence led to an 832 bp band and mutant sequence a 298 bp band.

The nature of the sspob1446 mutation was determined by whole genome sequencing. DNA was extracted from mutant embryos using a phenol/chloroform procedure. Libraries were prepared using the FS DNA Library Prep Kit for Illumina sequencing (NEB, E7805). DNA was digested into 150 bp fragments, and paired-end sequencing was performed using a NovaSeq 6000 Sequencing System. Trimmomatic (version 0.36; ILLUMNIACLIP: TruSeq3-PE-2.fa:2:30:10:1:true LEADING:3 TRAILING:3 SLIDINGWINDOW:5:20 MINLEN:42 AVGQUAL:30) was used to remove Illumina adaptor sequences from paired-end reads. Illumina short-read sequences were then aligned to the GRCz11 reference sequence of chromosome 24 using BWA-MEM (version 0.7.01). SAMtools (version 1.8) was used to sort and index reads. Aligned reads in BAM format were analyzed in IGV (version 2.13.1). Mutants were routinely genotyped by PCR amplification with oligos sspo_geno_1 and sspo_geno_2, followed by BsaI-HFv2 (NEB, R2733) restriction digestion to produce 300 bp and 99 bp bands from wild-type DNA and a protected 399 bp band from mutant DNA.

Generation of somatic mosaic F0 mutants (crispants)

Request a detailed protocol

Four gRNA oligos per gene for cfap298, sspo, urp1, and urp2 were chosen from a look-up table (Key resource appendix; Wu et al., 2018). gRNAs were synthesized from oligos in multiplex. After being pooled at 10 µM, oligos were annealed and extended with Bottom strand ultramer_2 using Phusion High-Fidelity PCR Mastermix (NEB, M0531) with Phusion High-Fidelity DNA Polymerase (NEB, M05030) using incubations: 98oC (2 min), 50°C (10 min), and 72oC (10 min). Assembled oligos were purified and used as templates for in vitro RNA synthesis, as described in ‘Generation and genotyping of mutant lines’ section. For mutagenesis, 1000 pg of gRNAs along with 1600 pg/nl Cas9 (IDT, 1081058) were injected into one-cell-stage embryos. To assess rates of mutagenesis, DNA was extracted from 1 dpf crispants and subjected to T7 endonuclease I assays (NEB, E3321).

Quantitation of body curvature at 1–2 dpf

Request a detailed protocol

Zebrafish larvae at 28–30 hpf were imaged using a Leica S9i stereomicroscope with integrated 10-megapixel camera. Body angles were calculated using ImageJ (Schindelin et al., 2012) as described in Bearce et al., 2022.

Quantitative reverse transcriptase PCR (qRT-PCR)

Request a detailed protocol

RNA was extracted using a Zymo Direct-Zol RNA Miniprep kit (Zymo Research, R2051). cDNA was prepared from 25 ng of RNA using oligoDT primers in a 20 µl reaction using a High Capacity cDNA Reverse Transcription Kit (ThermoFisher, 4368814). qRT-PCR reactions were performed in real time using 5 µl PowerUp SYBR Green Master Mix (ThermoFisher, A25741), 0.8 µl of 10 µM forward and reverse primers, 1.4 µl of nuclease-free water, and 2 µl of diluted cDNA. PCR was performed using a QuantStudio Real Time PCR System (Applied Biosystems) with cycling parameters: 50°C (2 min), 95°C (10 min) then 40 cycles of 95°C (15 s), and 60°C (1 min). Each reaction was performed in quadruplicate. Quantitation was relative to rpl13 and used the ∆∆CT relative quantitation method in which fold changes are calculated as 2−∆∆CT. The efficiency of amplification was verified to be close to 100% with a standard curve of RNA dilutions.

Calcein staining

Request a detailed protocol

Larvae were incubated in water containing 0.2% calcein (Sigma-Aldrich, C0875) for 10 min then rinsed two to three times in water (5 min per rinse). Larvae were immobilized with 0.005% tricaine, mounted in 0.8% low melt agarose, and imaged with a Leica THUNDER stereoscope.

Multiplex fluorescent in situ hybridization chain reaction (in situ HCR)

Request a detailed protocol

Embryos were fixed in 4% paraformaldehyde at 4°C overnight, washed with phosphate buffered saline (PBS) then serially dehydrated to 100% methanol, and stored at –20°C. Embryos were rehydrated, washed with PBS containing 0.1% Tween-20, incubated in hybridization buffer (Molecular Instruments), then incubated in 2 pmol of probes at 37°C overnight in a total volume of 500 µl of hybridization buffer. Embryos were washed in wash buffer (Molecular Instruments), washed twice in 5× SSCT (sodium chloride sodium citrate with 0.1% tween-20), and then incubated in amplification buffer (Molecular Instruments) for 1 hr. RNA hairpins designed to bind either pkd2l1, urp1, or urp2 were prepared by heating 10 pmol of each to 95°C for 90 s then snap-cooled in the dark for 30 min. Embryos were then incubated overnight in 500 µl of amplification buffer containing 30 pmol hairpins at room temperature in the dark. Embryos were washed five times in 5× SSCT, stored at 4oC, and then mounted for confocal microscopy. Images were acquired using a Zeiss LSM880 using either a ×20 air or ×40 water objective. Acquisition settings were derived using wild type embryos and then applied to all embryos. Images were exported to IMARIS 9.5.0 (Oxford Instruments). A Gaussian filter of width 0.42 µm (×20) or 0.21 µm (×40), and a rolling ball background subtraction of 10 µm was applied.

X-ray microcomputed tomography

Request a detailed protocol

Scans were performed using a vivaCT 80 (Scanco Medical) at 18.5-µm voxel resolution (for 3 mpf and 12 mpf fish) or 10-µm voxel resolution (1 mpf fish) as previously described (Bearce et al., 2022). Digital dissections of the spinal column were performed in 3D Slicer (Kikinis et al., 2013) using the Segmentation Editor. A threshold of 3200 was used to mask 1 mm tube (Draw Tube function) around the spine in the axial slice view, beginning between the otic vesicles rostral to the first vertebrae and ending at the split of the tail. The center of the tube was set at the narrowest ‘hollow”’ within the lumen of centra.

Quantitation of lateral spine curvature

Request a detailed protocol

Quantitation of lateral spine curvature was performed in ImageJ by orienting isolated spine images in a dorsal view with heads to the left and with the otic vesicles and first Weberian vertebrae parallel to the x-axis. Landmarks were assigned to the narrowest point of each centrum rostral to caudal; where maximum projection resulted in hidden or overlapping vertebrae, the appropriate number of points was added in closest approximation. The y-value from each landmark was subtracted from the point rostral to it, resulting in a map of local deflections where positive values indicate rightward displacement, and negative values indicate leftward displacement.

Live imaging of SCOspondin-GFP

Request a detailed protocol

Embryos (28 hpf) and larvae (12 dpf) were anesthetized in tricaine until touch response was abolished and then embedded in 0.8% low-melt agarose laced with tricaine in inverted imaging chambers (14 mm #1.5 coverslips, VWR cat no. 10810–054). In larvae, care was taken to align the posterior body close to the coverslip to the RF within the working distance of the objective. A Nikon Ti2 inverted microscope equipped with Plan Apo ×40 and ×60 WI DIC (1.2 NA) objectives, a Yokogawa Spinning Disk and pco.edge sCMOS camera were used to capture 512×256 images in time series. Exposure time varied with age (100 ms–300 ms) as Sspo-GFP brightened in intensity over time; exposure, camera settings, and laser power were kept constant between age-matched individuals. Images were cropped, rotated, and intensity-adjusted in ImageJ (Schindelin et al., 2012).

Appendix 1

Appendix 1—key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
Commercial assay, kitDNA Clean and Concentrator KitZymo ResearchCat no: D4013
Commercial assay, kitRNA Clean and Concentrator KitZymo ResearchCat no: R1016
Commercial assay, kitDirect-zol RNA MiniPrep KitZymo ResearchCat no: R2050
Commercial assay, kitGeneJET Gel Extraction KitThermo Fisher ScientificCat no: K0691
Commercial assay, kitHigh Capacity RNA-to-cDNA KitThermo Fisher ScientificCat no: 4387406
Commercial assay, kitMEGAshortscript T7 Transcription KitThermo Fisher ScientificCat no: AM1354
Commercial assay, kitHiScribe T7 High Yield
RNA Synthesis Kit
New England BiolabsCat no: E2040
Commercial assay, kitFS DNA Library Prep
Kit for Illumina
New England BiolabsCat no: E7805
Commercial assay, kitHCR-RNA FISH
Hybridization Buffer
Molecular Instruments
Commercial assay, kitT7 Endonuclease INew England BiolabsCat no: E3321
Commercial assay, kitHCR-RNA FISH
Amplification Buffer
Molecular Instruments
Commercial assay, kitAlexaFluor-647 HairpinsMolecular Instruments
Commercial assay, kitAlexaFluor-546 HairpinsMolecular Instruments
Commercial assay, kitAlexaFluor-488 HairpinsMolecular Instruments
Commercial assay, kiturp1 HCR-RNA FISH probeMolecular InstrumentsProject-specific design
Commercial assay, kiturp2 HCR-RNA FISH probeMolecular InstrumentsProject-specific design
Commercial assay, kitpkd2l1 HCR-RNA FISH probeMolecular InstrumentsProject-specific design
Commercial assay, kitTaq PolymeraseNew England BiolabsCat no: M0273
Commercial assay, kitTURBO DNaseThermo Fisher ScientificCat no: AM2238
Commercial assay, kitPhusion High-Fidelity
DNA Polymerase
New England BiolabsCat no: M0530
Commercial assay, kitPhusion High-Fidelity
PCR Master Mix
New England BiolabsCat no: M0531
Commercial assay, kitSYBR Green PCR Master MixThermo Fisher ScientificCat no: 4309155
Strain, strain background (Danio rerio)AB strainUniversity of Oregon
Strain, strain background (D. rerio)WIK strainUniversity of Oregon
Strain, strain background (D. rerio)TU strainUniversity of Oregon
Strain, strain background (D. rerio)cfap298tm304 lineJaffe et al., 2016ZDB-FISH-150901–23024
Strain, strain background (D. rerio)pkd2l1icm02 lineSternberg et al., 2018ZDB-FISH-160811–9
Strain, strain background (D. rerio)sspob1446 lineThis studyFigure 1—figure supplement 3
Strain, strain background (D. rerio)sspo-GFPut24 lineTroutwine et al., 2020ZDB-FISH-190313–20
Strain, strain background (D. rerio)urp1b1420 lineThis studyFigure 1—figure supplement 2
Strain, strain background (D. rerio)urp2b1421 lineThis studyFigure 1—figure supplement 2
Strain, strain background (D. rerio)uts2r3b1436 lineThis studyFigure 2—figure supplement 3
Software, algorithm3D SlicerFedorov et al., 2012
Software, algorithmIMARIS 9.9Oxford Instruments
Software, algorithmNIS-ElementsNikon Instruments Inc
Software, algorithmZEN SoftwareCarl Zeiss AG
Software, algorithmQuantStudio Design
and Analysis Software
Applied Biosystems
Software, algorithmIntegrated Genome
Viewer (version 2.13.1)
Robinson et al., 2011
Sequence-based reagentBottom strand ultramer_1This studyTail ultramer for generating singe gRNA oligosAAAAGCACCGACTCG
GTGCCACTTTTTC
AAGTTGATAACGGACT
AGCCTTATTTT
AACTTGCTAT
Sequence-based reagenturp1_gRNA_1This studygRNA_1 oligo for generating urp1b1420 linetaatacgactcactataGGCGT
TGGTCAGCCTGACAT
gttttagagctagaa
Sequence-based reagenturp1_gRNA_2This studygRNA_2 oligo for generating urp1b1420 linetaatacgactcactataGGG
TCCTCTGTCCATCTCCG
gttttagagctagaa
Sequence-based reagenturp2_gRNA_1This studygRNA_1 oligo for generating urp2b1421 linetaatacgactcactataGGCA
GATGGAGAAAGATTGA
gttttagagctagaa
Sequence-based reagenturp2_gRNA_2This studygRNA_2 oligo for generating urp2b1421 linetaatacgactcactataGGC
GTTTGCAGAAATCAGCG
gttttagagctagaa
Sequence-based reagentuts2r3_gRNAThis studygRNA oligo for generating uts2r3b1436 linetaatacgactcactataGGG
TGAAGGGGAAGAGAAGA
gttttagagctagaa
Sequence-based reagentsspo_gRNA_1This studygRNA_1 oligo for generating sspob1446 linetaatacgactcactataGGT
CCCCAGTGGTCCGCG
GTgttttagagctagaa
Sequence-based reagentsspo_gRNA_2This studygRNA_2 oligo for generating sspob1446 linetaatacgactcactataGGCAC
AGTGTGTGAGACCAG
gttttagagctagaa
Sequence-based reagentBottom strand ultramer_2This studyTail ultramer for generating multiplexed gRNA oligosAAAAGCACCGACTCGG
TGCCACTTTTTC
AAGTTGATAA
CGGACTAGCCTTATT
TTAACTTGCTATTTC
TAGCTCTAAAAC
Sequence-based reagenturp1_F0_gRNA_1Wu et al., 2018gRNA_1 oligo for generating urp1 F0 embryosTAATACGACTCACTAT
AGGAAAGTGAAGAT
CGCGGCCGTTTTAG
AGCTAGAAATAGC
Sequence-based reagenturp1_F0_gRNA_2Wu et al., 2018gRNA_2 oligo for generating urp1 F0 embryosTAATACGACTCACT
ATAGGACACGG
CTCTGCC
ACAACGTTTTAGA
GCTAGAAATAGC
Sequence-based reagenturp1_F0_gRNA_3Wu et al., 2018gRNA_3 oligo for generating urp1 F0 embryosTAATACGACTCACTA
TAGGTTCAGAAGC
TGGTAGCAGGTTTTA
GAGCTAGAAATAGC
Sequence-based reagenturp1_F0_gRNA_4Wu et al., 2018gRNA_4 oligo for generating urp1 F0 embryosTAATACGACTCACT
ATAGGGAAAATAAAT
AACATGGTGTTTTA
GAGCTAGAAATAGC
Sequence-based reagenturp2_F0_gRNA_1Wu et al., 2018gRNA_1 oligo for generating urp2 F0 embryosTAATACGACTCACTA
TAGGTGACTGTCGC
TTCAATCGGTTTTAG
AGCTAGAAATAGC
Sequence-based reagenturp2_F0_gRNA_2Wu et al., 2018gRNA_2 oligo for generating urp2 F0 embryosTAATACGACTCACT
ATAGGGACATTTCCT
GACGGAGAGTTTTA
GAGCTAGAAATAGC
Sequence-based reagenturp2_F0_gRNA_3Wu et al., 2018gRNA_3 oligo for generating urp2 F0 embryosTAATACGACTCACTA
TAGGTGGACACGA
GGAGACCGAGTTTT
AGAGCTAGAAATAGC
Sequence-based reagenturp2_F0_gRNA_4Wu et al., 2018gRNA_4 oligo for generating urp2 F0 embryosTAATACGACTCACT
ATAGGTCACCAGGTAG
TGACGGAGTTTTAG
AGCTAGAAATAGC
Sequence-based reagentsspo_F0_gRNA_1Wu et al., 2018gRNA_1 oligo for generating sspo F0 embryosTAATACGACTCACT
ATAGGTTCGTCCCC
AGTGGTCCGGTTTT
AGAGCTAGAAATAGC
Sequence-based reagentsspo_F0_gRNA_2Wu et al., 2018gRNA_2 oligo for generating sspo F0 embryosTAATACGACTCACT
ATAGGAAACGG
CCGTCAGTGTCGGT
TTTAGAGCTAGAAATAGC
Sequence-based reagentsspo_F0_gRNA_3Wu et al., 2018gRNA_3 oligo for generating sspo F0 embryosTAATACGACTCAC
TATAGGTGTTGC
AACACCAACCGGGT
TTTAGAGCTAGAAATAGC
Sequence-based reagentsspo_F0_gRNA_4Wu et al., 2018gRNA_4 oligo for generating sspo F0 embryosTAATACGACTCACT
ATAGGAGCCTAGACC
TGCTCACGGTTTTA
GAGCTAGAAATAGC
Sequence-based reagentcfap298_F0_gRNA_1Wu et al., 2018gRNA_1 oligo for generating cfap298 F0 embryosTAATACGACTCAC
TATAGGTTCTCTT
CAACACTACGGGT
TTTAGAGCTAGAAATAGC
Sequence-based reagentcfap298_F0_gRNA_2Wu et al., 2018gRNA_2 oligo for generating cfap298 F0 embryosTAATACGACTCAC
TATAGGGCTCC
ACAATCTGATCATG
TTTTAGAGCTA
GAAATAGC
Sequence-based reagentcfap298_F0_gRNA_3Wu et al., 2018gRNA_3 oligo for generating cfap298 F0 embryosTAATACGACTCAC
TATAGGCATTC
TTATTGGATCATGG
TTTTAGAGCTA
GAAATAGC
Sequence-based reagentcfap298_F0_gRNA_4Wu et al., 2018gRNA_4 oligo for generating cfap298 F0 embryosTAATACGACTCACT
ATAGGTCTCTGG
CAGGTGCGCCCGTT
TTAGAGCTAGAAATAGC
Sequence-based reagentpkd2l1_geno_1Sternberg et al., 2018pkd2l1icm02 genotyping oligo 1TGTGTGCTAGG
ACTGTGGGG
Sequence-based reagentpkd2l1_geno_2Sternberg et al., 2018pkd2l1icm02 genotyping oligo 2AGGGCAAGAGAA
TGGCAAGACG
Sequence-based reagenturp1_geno_1This Studyurp1b1420 genotyping oligo 1GCACCCAAAAT
CCAACGACT
Sequence-based reagenturp1_geno_2This Studyurp1b1420 genotyping oligo 2TGTATGGGGAA
AACAAAGGCA
Sequence-based reagenturp2_geno_1This Studyurp2b1421 genotyping oligo 1TTGGGGTTGT
AACAGGTAGTG
Sequence-based reagenturp2_geno_2This Studyurp2b1421 genotyping oligo 2AACAAGGAAGA
CGCTGCAAG
Sequence-based reagentuts2r3_geno_1This Studyuts2r3b1436 genotyping oligo 1ATGGATCCCC
TGATGTCCTG
Sequence-based reagentuts2r3_geno_2This Studyuts2r3b1436 genotyping oligo 2TCGAACTCTGC
TCATCCCAG
Sequence-based reagentsspo_geno_1This Studysspob1446 genotyping oligo 1CGCAAACACTT
CCACTTCCA
Sequence-based reagentsspo_geno_2This Studysspob1446 genotyping oligo 2TTGAAGCCAGATGT
AAAGGATGAGTGT
Sequence-based reagenturp1_gRNA1+2_T7E1_FThis StudyForward primer to amplify genomic DNA for T7E1 assay in crispantsGACAGCGCAC
CCTTAATTGT
Sequence-based reagenturp1_gRNA1+2_T7E1_RThis StudyReverse primer to amplify genomic DNA for T7E1 assay in crispantsACATTTAGCCTT
AACAAGCACAA
Sequence-based reagenturp1_gRNA3+4_T7E1_FThis StudyForward primer to amplify genomic DNA for T7E1 assay in crispantsCAGACAAGGG
AACAGAGAGGA
Sequence-based reagenturp1_gRNA3+4_T7E1_RThis StudyReverse primer to amplify genomic DNA for T7E1 assay in crispantsCCACTGCTTTTA
AATCATCCACC
Sequence-based reagenturp2_gRNA1_T7E1_FThis StudyForward primer to amplify genomic DNA for T7E1 assay in crispantsATCTTAGAGG
CGCATTGGTG
Sequence-based reagenturp2_gRNA1_T7E1_RThis StudyReverse primer to amplify genomic DNA for T7E1 assay in crispantsGCATGAGGCG
GTTTGTTTTG
Sequence-based reagenturp2_gRNA2+3_T7E1_FThis StudyForward primer to amplify genomic DNA for T7E1 assay in crispantsTGAAGCAACT
GAGGAGCAAA
Sequence-based reagenturp2_gRNA2+3_T7E1_RThis StudyReverse primer to amplify genomic DNA for T7E1 assay in crispantsACAGTACAGTT
CAGCACACCT
Sequence-based reagenturp2_gRNA4_T7E1_FThis StudyForward primer to amplify genomic DNA for T7E1 assay in crispantsTGACCTATACATC
AAAGCCAAGG
Sequence-based reagenturp2_gRNA4_T7E1_RThis StudyReverse primer to amplify genomic DNA for T7E1 assay in crispantsCCTGGGCTGA
TCATACCTCT
Sequence-based reagentrpl13_qPCR_FThis StudyForward primer for quantitative RT-PCRTAAGGACGGAG
TGAACAACCA
Sequence-based reagentrpl13_qPCR_RThis StudyReverse primer for quantitative RT-PCRCTTACGTCTGC
GGATCTTTCTG
Sequence-based reagenturp1_qPCR_FThis StudyForward primer for quantitative RT-PCRACATTCTGG
CTGTGGTTTG
Sequence-based reagenturp1_qPCR_RThis StudyReverse primer for quantitative RT-PCRGTCCGTCTTCA
ACCTCTGCTAC
Sequence-based reagenturp2_qPCR_FThis StudyForward primer for quantitative RT-PCRAGAGGAAACA
GCAATGGACG
Sequence-based reagenturp2_qPCR_RThis StudyReverse primer for quantitative RT-PCRTGTTGGTTTTG
GTTGACG

Data availability

All data generated or analysed during this study are included in the manuscript and supporting file.

References

    1. Mesiti BL
    (2021)
    Scoliosis: an overview
    Radiologic Technology 93:55–72.
    1. Wishart BD
    2. Kivlehan E
    (2021) Neuromuscular scoliosis
    Physical Medicine and Rehabilitation Clinics of North America 32:547–556.
    https://doi.org/10.1016/j.pmr.2021.02.007

Decision letter

  1. Marianne E Bronner
    Senior and Reviewing Editor; California Institute of Technology, United States

In the interests of transparency, eLife publishes the most substantive revision requests and the accompanying author responses.

[Editors' note: this paper was reviewed by Review Commons.]

https://doi.org/10.7554/eLife.83883.sa1

Author response

General Statements

We were immensely pleased that the reviewers considered our conclusions “well supported” and our study “beautifully executed”. Reviewers also recognized the significance of our work. Reviewer 1 stated that “building a model that describes one of these pathways will allow us to begin to test therapies to treat or prevent scoliosis” then noted that we “help to build a larger model of normal spine morphogenesis” and that this is “important”. Reviewer 2 called our work an “exciting advance in our understanding of one of the essential signaling pathways that help regulate body axis straightening and spine morphogenesis in zebrafish” and mentioned that our work “may also help to further our understanding of the etiology and pathophysiology of multiple forms of neuromuscular scoliosis in humans”. Reviewer 3 agreed, stating that our work “adds important information on the role of urotensin signaling in spine formation” and noted that it is timely: “findings are of special significance in the light of recent reports that mutations in UTS2R3 show association with spinal curvature in patients with adolescent idiopathic scoliosis”.

We thank the three reviewers for reading our research and providing feedback. In all cases, we have incorporated their suggestions, and we believe this has made our manuscript much stronger. Indeed, reviewers had only a small number of “major points”, and all have now been addressed. We have also addressed all “minor points” raised by reviewers.

Major points

Reviewer 1

The key conclusions are well supported, see below for my two major issues.

1. Please don't call this lordosis. Lordosis or hyperlordosis effects lumbar vertebra. The curve in the lumbar region shifts body weight so that human gait is more efficient that in the great apes, or so the story goes. Zebrafish do not have lumbar vertebra equivalents or a natural curve in the caudal region. Similarly, fish do not have the equivalent vertebra to generate kyphosis, which is again a hyper flexion of a normal human spinal curve. Instead zebrafish have Weberian, precaudal and caudal vertebra. It would be so much more useful for the field if the authors used these terms and specified ranges, i.e. numbered vertebrae, that are effected so we can directly and accurately compare regions of defects between zebrafish mutants. It would help to make the point that the uts2r3 mutant has more caudally located curves than urp1/2 double mutants.

We agree with the reviewer and have made changes in line with their suggestions. For example, at the end of the abstract, we have changed “new animal models of lordosis-like curves” to “new animal models of spine deformity”; at the end of the Introduction, we say “mutants model human spinal deformity” rather than “mutants model a lordosis-like condition”.

We have also added a new paragraph to the Discussion about whether zebrafish can truly model scoliosis, lordosis and kyphosis which, among other points, mentions some of what the reviewer brought up about lumbar vs. caudal vertebral differences.

When we do mention lordosis (less often than previously) in the text, we are careful to say “lordosis-like” and we fully explain what we mean by this, including the aspects of the zebrafish phenotypes that are not lordosis-like. Last, when describing our phenotypes, we now use the appropriate names for zebrafish vertebrae, i.e., caudal rather than lumbar.

2. The observation that urp1/2 double mutants have curves only in the D/V plane and almost completely lack side-to-side curves is noteworthy. Does the urp1-/-urp2-/- mutant uncouple two systems for posture? If this separate a DV from side-to-side postural control system, that would be very interesting. It is particularly important to describe how penetrant the phenotype is and how many times it was observed. See 9 minor comments. It would help the reader if the authors explicitly described the features that they see in the cfap298 mutant that constitute lateral curves and that are lacking in urp1/2 (e.g. in figure 4E).

To address this comment, we have supplied more data on the extent of lateral curvature in cfap298tm304 and urp1∆P;urp2∆P mutants and, in the new version, also uts2r3b1436 and cfap298tm304;urp1∆P;urp2∆P triple mutants. These new data are:

New data showing the extent of lateral curvature in various mutant conditions (manuscript Figure 4Fi-Fii).

Additionally, in the interests of making all of our data available, we provide a supplementary figure showing all scans performed that were used to analyze lateral curvature where we also note the sex of the fish. See new Figure 4—figure supplement 1.

In agreement with our previous observations, cfap298tm304 mutants indeed exhibit significantly more lateral curvature than urp1∆P;urp2∆P or uts2r3b1436 mutants. Thus, these mutants may uncouple two systems for posture, as the reviewer suggests. We now note this possibility in the Discussion.

Reviewer 2

1. Need to show that the CRISPANT targeting was effective for mutagenesis at each loci screened in the work presented in Figure 1E.

We have now performed T7 endonuclease assays (as also suggested by Reviewer 3) to assess mutagenesis in crispants. These assays revealed significant insertion-deletion mutations at all gRNA targeting sites as well as cases of large deletions between sites.

New data demonstrating that urp1 and urp2 targeting was effective for mutagenesis (see also Figure 1—figure supplement 2C-D):

These data are referenced in the main text in the following way: “Using T7 endonuclease assays, we confirmed that high levels of insertion-deletion mutations were generated at gRNA sites in crispants (Figure 1—figure supplement 2C-D).”

Overall, the evidence we have that Urp1 and Urp2 peptides do not contribute to axial straightening is:

  • urp1∆P and urp2∆P single and double mutants (and maternal-zygotic mutants) show normal axial straightening. These mutants completely lack the coding region for the Urp1 and Urp2 peptides.

  • urp1 and urp2 crispants with significant indels (as assessed by T7E1 assays) also showed normal straightening. By contrast, all other positive control crispants showed expected axial phenotypes.

  • A preprint produced independently of our work also generated urp1 and urp2 mutants and also reported lack of axial defects (Gaillard et al., 2022 – BioRxiv).

Thus, we conclude that Urp1 and Urp2 are dispensable for axial straightening in zebrafish. We devote several paragraphs in our Discussion to this somewhat unexpected finding.

Reviewer 3

1. The addition of the F0 crispant experiment to show that the pro-peptide of urp1/2 does not have a function and is responsible for the difference between the observed morpholino and the crispr phenotype was important. However, since no phenotype was observed in crispants it is important to add evidence of induced cuts for all guide RNAs used in the crispant experiment. These control experiments might have been done already. If not, they can easily be done in a short period of time by performance of T7 assays on injected fish and would not require additional reagents.

Please see response to Reviewer 2 above including Figure1-figure supplement 2. In summary, we have now performed the suggested T7 assays and they indicated that indels were indeed generated in urp1 and urp2 crispants (in the manuscript, now described in Figure 1—figure supplement 2C-D).

Overall, the evidence we have that Urp1 and Urp2 peptides do not contribute to axial straightening is:

  • urp1∆P and urp2∆P single and double mutants (and maternal-zygotic mutants) show normal axial straightening. These mutants completely lack the coding region for the Urp1 and Urp2 peptides.

  • urp1 and urp2 crispants with significant indels (as assessed by T7E1 assays) also showed normal straightening. By contrast, all other positive control crispants showed expected axial phenotypes.

  • A preprint produced independently of our work also generated urp1 and urp2 mutants and also reported lack of axial defects (Gaillard et al., 2022 – BioRxiv).

2. The authors claim that there were no structural defects observed in urp1/2 double mutants. However, the hemal arch in figure 3 E seems to be deformed. This could be normal variance or a phenotype. This can be addressed by simple reinspection of the scans.

Reinspection of the scans shows no significant hemal arch defects in urp1∆P;urp2∆P mutants. We now explicitly state this in the text. We also note that we do not see significant vertebral defects such as fusions or hemivertebrae and nor do we see any patterning defects at 10 dpf based on calcein staining. Additionally, we have added vertebral body shape quantitation (length and height aspect ratios):

New data showing vertebral body shape quantitation (see also Figure 3F)

Since these structural changes are subtle, and we do not observe vertebral patterning changes based on calcein staining at earlier time points, we conclude that the vertebral structural changes are not causative in spinal curvature but rather are a consequence of the presence of curves.

Minor points

Reviewer 1

3. Supplementary FigS3B How to measure the Cobb Angle is unclear. Why is the first curve not counted? I count 3 curves. First a ventral displacement, then a dorsal to ventral return, then a sharp flex before the tail. How to measure Cobb angle might be easier to explain if the figure is expanded into steps. Identify the apical vertebra, then showing how the lines are drawn parallel to those vertebrae, then where the measured angle forms between the lines perpendicular to the drawn parallel lines.

We have now redrawn our schematic figure more thoroughly explained how Cobb angle is measured with the following figure and figure legend.

5a. I think we (zebrafish biologists) need be explicit about what we mean with "without vertebral defects." What do we count as defects? Vertebrae can be fused, bent, shortened or the growing edges can be slanted. In Figure 3E, and movie7, it is clear that the highlighted mutant vertebrae are shorter than WT. The growing ends of normal vertebra are perpendicular to the long axis of the vertebra. In the mutants the ends are slanted. Please define in the text what you consider a relevant vertebral defect, because these vertebrae have defects. Or are you only considering the calcein stained centra at 10dpf?

We have made textual edits to clearly state that we do not see vertebral fusions or missing appendages. We have also added vertebral body length and height aspect ratios to quantify vertebrae shape and find subtle changes in the variance of vertebral body length aspect ratios in urp1∆P;urp2∆P mutants (see Figure 3F in the new manuscript). Since these structural changes are subtle, and we do not observe vertebral patterning changes based on calcein staining at earlier time points, we conclude that the vertebral structural changes are not causative in spinal curvature but rather are a consequence of the presence of curves.

5b. Do you want to base your patterning conclusion on primarily the calcein data as these are closer to the notochord patterning time window. Please anchor this conclusion to a specific time or standard length e.g. 10dpf/5.6mm.

We now note the time window when calcein was performed and provide the standard lengths of fish analyzed. As such, our conclusion about vertebral patterning is now anchored to the standard length, as the reviewer suggested:

“Staining of juveniles with the vital dye calcein revealed no defects in vertebral patterning or spacing in urp1∆P;urp2∆P mutants at 10 dpf (4.5-6 mm standard length [Figure 3C]) (Figure 3F).”

6. "At 30 dpf… several mutants exhibited a significant curve in the pre-caudal vertebrae, in addition to a caudal curve (Figure 3D and S3C). Since pre-caudal curves were rare in mutants at 3-months, this suggested that curve location is dynamic". The frequency of this observation is important. Does it effect all or a fraction of mutants? Can you provide some numbers to anchor these observations? Maybe fractions e.g.. 3 of 4 fish had precaudal curves at 30pdf, and 0 of 10 fish had precaudal curves by 3 mpf?

We have now analyzed more 1 mpf fish by µCT and we observe that 5 of 7 exhibit pre-caudal curves:

New data showing spinal curves are variable in 1 mpf urp1∆P;urp2∆P mutants (see also Figure 3—figure supplement 1).

These data are also mentioned in our revised text. By contrast to the significant pre-caudal curves at 1 mpf, 3 mpf fish showed predominantly caudal curves (Figure 2F). Only 2 of 9 fish analyzed by µCT at 3 mpf were scored as exhibiting an “apex of curve” vertebrae in the precaudal vertebrae, and in both of those cases the curve was at the very end of the precaudals. This reinforces our conclusion that curve position is dynamic through growth.

7. The description of the pkd2l1 mutant, instead of terming it kyphosis can you tell the reader the vertebra number at the peak of the curve. The authors say the pkd2l1 mutant is highly distinct from urp1/urp2-/-, but the reader needs to hear exactly what is distinct. For example, does this mutant have both lateral and D/V curves?

We have now scanned several pkd2l1icm02 mutant fish include images of pkd2l1icm02 mutants at two different timepoints (3 mpf and 12 mpf; see new Figure 2—figure supplement 2). Our results agreed with those previously published for the pkd2l1icm02 mutant line (Sternberg et al., 2018) but we believe it is important to include our data so readers can see side-by-side images of the various mutant conditions using the same skeleton visualization technique.

At 3 mpf, pkd2l1icm02 mutants essentially appeared wild-type but by 12 mpf they had developed very subtle D/V curves in the pre-caudal vertebrae. They do not exhibit M/L curves at either stage.

We called the phenotype displayed by pkd2l1icm02 mutants “kyphosis” to be in line with a previous publication describing these mutants (Sternberg et al., 2018). We have now added new wording in the Discussion about whether or not zebrafish can truly model kyphosis and lordosis, and we now make clear in our Results that the “mutants went on to develop very subtle kyphosis-like curves” rather than “is kyphosis”.

It is intriguing that pkd2l1icm02 mutants do not exhibit any curves until much later in life than urp1∆P;urp2∆P and uts2r3b1436 mutants. Inspired by this finding, we aged urp1∆P and urp2∆P single mutants and found that they go on to develop D/V curves by 12 mpf, with urp1∆P mutants being less severe than urp2∆P mutants (Figure 2—figure supplement 2). To summarize all of our findings:

Author response table 1
3-months12-monthsPosition of curve
urp1∆Pno curvesmild D/V curvesMostly caudal
urp2∆Pmild D/V curvesintermediate D/V curvesMostly caudal
urp1∆P;urp2∆Psevere D/V curvessevere D/V curvesMostly caudal
uts2r3b1436severe D/V curvessevere D/V curvesMostly caudal
cfap298tm304severe 3D curvessevere 3D curvesCaudal and pre-caudal
pkd2l1icm02no curvesvery mild D/V curvesMostly pre-caudal

Phenotypes in urp1∆P and urp2∆P single mutants upon aging shows: (1) Urp1 and Urp2 are not entirely redundant in long-term spine maintenance and (2) proper Urp1/Urp2 dose is essential. We have now included these new insights.

This is an interesting question. To address it, we imaged more uts2r3b1436 mutant spines and reconstructed views from dorsal aspect, including these in new Figure 4Fi and Figure 4—figure supplement 1. We also quantified the degree of lateral curvature (Figure 4Fii). The reviewer’s suggestion is correct – there are minimal side-to-side curves in uts2r3b1436 mutants, highly similar to what we found for urp1∆P;urp2∆P double mutants.

9. One finding that deserves more discussion is the observation that urp1/urp2 double mutants have almost no side-to-side defects and all the obvious bends are in the D/V plane. Does this uncouple two systems for posture? Please consider the following paper. It shows a proprioception system that maintains normal side-to-side posture. A spinal organ of proprioception for integrated motor action feedback.

Picton LD, Bertuzzi M, Pallucchi I, Fontanel P, Dahlberg E, Björnfors ER, Iacoviello F, Shearing PR, El Manira A. Neuron. 2021 Apr 7;109(7):1188-1201.e7. doi: 10.1016/j.neuron.2021.01.018. Epub 2021 Feb 11. PMID: 33577748

Thank you for pointing out this manuscript. We have now mentioned this potential “uncoupling” in our Discussion and included this citation.

Reviewer 2

1. Figure 3F: might be improved by making the images black and white and possibly inverted. It is not easy to clearly see the vertebrae as is.

Thanks for the suggestion, we made this change.

https://doi.org/10.7554/eLife.83883.sa2

Article and author information

Author details

  1. Elizabeth A Bearce

    Institute of Molecular Biology, Department of Biology, University of Oregon, Eugene, United States
    Contribution
    Formal analysis, Funding acquisition, Investigation, Visualization, Methodology, Writing - review and editing
    Competing interests
    No competing interests declared
  2. Zoe H Irons

    Institute of Molecular Biology, Department of Biology, University of Oregon, Eugene, United States
    Contribution
    Funding acquisition, Validation, Methodology
    Competing interests
    No competing interests declared
  3. Johnathan R O'Hara-Smith

    Institute of Molecular Biology, Department of Biology, University of Oregon, Eugene, United States
    Contribution
    Investigation, Visualization, Methodology, Writing - review and editing
    Competing interests
    No competing interests declared
  4. Colin J Kuhns

    Institute of Molecular Biology, Department of Biology, University of Oregon, Eugene, United States
    Contribution
    Validation, Investigation
    Competing interests
    No competing interests declared
  5. Sophie I Fisher

    Institute of Molecular Biology, Department of Biology, University of Oregon, Eugene, United States
    Contribution
    Investigation
    Competing interests
    No competing interests declared
  6. William E Crow

    Institute of Molecular Biology, Department of Biology, University of Oregon, Eugene, United States
    Contribution
    Investigation
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-2991-3076
  7. Daniel T Grimes

    Institute of Molecular Biology, Department of Biology, University of Oregon, Eugene, United States
    Contribution
    Conceptualization, Supervision, Funding acquisition, Investigation, Writing - original draft, Project administration, Writing - review and editing
    For correspondence
    dtgrimes@uoregon.edu
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-0173-1887

Funding

National Institutes of Health (R00AR70905)

  • Daniel T Grimes

National Institutes of Health (F32AR078002)

  • Elizabeth A Bearce

National Institutes of Health (F31HD105435)

  • Zoe H Irons

National Institutes of Health (R35GM142949)

  • Daniel T Grimes

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We thank Judy Peirce, Tim Mason, and the Aquatics Facility for zebrafish husbandry, the GC3F Biological Imaging Facility, and the X-Ray Imaging Core, all at the University of Oregon. We thank Claire Wyart, John Postlethwait, and Ryan Gray for sharing zebrafish lines, Zac Bush for help with genome sequencing, Mike Harms and Ron Kwon for discussions, and Katie Fisher for proof reading. Funding: This study was supported by a National Institutes of Health grants R00AR70905 and R35GM142949 (to DTG), F32AR078002 (to EAB), F31HD105435 and T32HD007348 (to ZHI).

Ethics

Experiments were undertaken in accordance with research guidelines of the International Association for Assessment and Accreditation of Laboratory Animal Care and approved by the University of Oregon Institutional Animal Care and Use Committee (# AUP-21-45).

Senior and Reviewing Editor

  1. Marianne E Bronner, California Institute of Technology, United States

Publication history

  1. Preprint posted: August 15, 2022 (view preprint)
  2. Received: October 2, 2022
  3. Accepted: November 24, 2022
  4. Accepted Manuscript published: December 1, 2022 (version 1)
  5. Version of Record published: January 12, 2023 (version 2)

Copyright

© 2022, Bearce et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 462
    Page views
  • 107
    Downloads
  • 1
    Citations

Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Elizabeth A Bearce
  2. Zoe H Irons
  3. Johnathan R O'Hara-Smith
  4. Colin J Kuhns
  5. Sophie I Fisher
  6. William E Crow
  7. Daniel T Grimes
(2022)
Urotensin II-related peptides, Urp1 and Urp2, control zebrafish spine morphology
eLife 11:e83883.
https://doi.org/10.7554/eLife.83883

Further reading

    1. Developmental Biology
    Marianne E Emmert, Parul Aggarwal ... Roger Cornwall
    Research Article Updated

    Neonatal brachial plexus injury (NBPI) causes disabling and incurable muscle contractures that result from impaired longitudinal growth of denervated muscles. This deficit in muscle growth is driven by increased proteasome-mediated protein degradation, suggesting a dysregulation of muscle proteostasis. The myostatin (MSTN) pathway, a prominent muscle-specific regulator of proteostasis, is a putative signaling mechanism by which neonatal denervation could impair longitudinal muscle growth, and thus a potential target to prevent NBPI-induced contractures. Through a mouse model of NBPI, our present study revealed that pharmacologic inhibition of MSTN signaling induces hypertrophy, restores longitudinal growth, and prevents contractures in denervated muscles of female but not male mice, despite inducing hypertrophy of normally innervated muscles in both sexes. Additionally, the MSTN-dependent impairment of longitudinal muscle growth after NBPI in female mice is associated with perturbation of 20S proteasome activity, but not through alterations in canonical MSTN signaling pathways. These findings reveal a sex dimorphism in the regulation of neonatal longitudinal muscle growth and contractures, thereby providing insights into contracture pathophysiology, identifying a potential muscle-specific therapeutic target for contracture prevention, and underscoring the importance of sex as a biological variable in the pathophysiology of neuromuscular disorders.

    1. Developmental Biology
    2. Genetics and Genomics
    Ankit Sabharwal, Mark D Wishman ... Stephen C Ekker
    Research Advance Updated

    The clinical and largely unpredictable heterogeneity of phenotypes in patients with mitochondrial disorders demonstrates the ongoing challenges in the understanding of this semi-autonomous organelle in biology and disease. Previously, we used the gene-breaking transposon to create 1200 transgenic zebrafish strains tagging protein-coding genes (Ichino et al., 2020), including the lrpprc locus. Here, we present and characterize a new genetic revertible animal model that recapitulates components of Leigh Syndrome French Canadian Type (LSFC), a mitochondrial disorder that includes diagnostic liver dysfunction. LSFC is caused by allelic variations in the LRPPRC gene, involved in mitochondrial mRNA polyadenylation and translation. lrpprc zebrafish homozygous mutants displayed biochemical and mitochondrial phenotypes similar to clinical manifestations observed in patients, including dysfunction in lipid homeostasis. We were able to rescue these phenotypes in the disease model using a liver-specific genetic model therapy, functionally demonstrating a previously under-recognized critical role for the liver in the pathophysiology of this disease.