SiliFish: A desktop application to model swimming behavior in developing zebrafish (Danio rerio).

SiliFish is an open-source desktop application to model and study zebrafish swimming. Here, we explain how to define the general parameters of the model, define cell populations, place them within the spinal cord, and define their projections. We explain how to run a simulation and how to visualize the network output and single-cell activity. The choice of C# as the programming language allows higher speed performance, simulating models with larger spinal circuits in less time. For complete details on the use and execution of this protocol, please refer to Roussel et al. (2021).1.

Mapping of morpho-electric features to molecular identity of cortical inhibitory neurons.

Knowledge of the cell-type-specific composition of the brain is useful in order to understand the role of each cell type as part of the network. Here, we estimated the composition of the whole cortex in terms of well characterized morphological and electrophysiological inhibitory neuron types (me-types). We derived probabilistic me-type densities from an existing atlas of molecularly defined cell-type densities in the mouse cortex. We used a well-established me-type classification from rat somatosensory cortex to populate the cortex. These me-types were well characterized morphologically and electrophysiologically but they lacked molecular marker identity labels. To extrapolate this missing information, we employed an additional dataset from the Allen Institute for Brain Science containing molecular identity as well as morphological and electrophysiological data for mouse cortical neurons. We first built a latent space based on a number of comparable morphological and electrical features common to both data sources. We then identified 19 morpho-electrical clusters that merged neurons from both datasets while being molecularly homogeneous. The resulting clusters best mirror the molecular identity classification solely using available morpho-electrical features. Finally, we stochastically assigned a molecular identity to a me-type neuron based on the latent space cluster it was assigned to. The resulting mapping was used to derive inhibitory me-types densities in the cortex.

A method to estimate the cellular composition of the mouse brain from heterogeneous datasets.

The mouse brain contains a rich diversity of inhibitory neuron types that have been characterized by their patterns of gene expression. However, it is still unclear how these cell types are distributed across the mouse brain. We developed a computational method to estimate the densities of different inhibitory neuron types across the mouse brain. Our method allows the unbiased integration of diverse and disparate datasets into one framework to predict inhibitory neuron densities for uncharted brain regions. We constrained our estimates based on previously computed brain-wide neuron densities, gene expression data from in situ hybridization image stacks together with a wide range of values reported in the literature. Using constrained optimization, we derived coherent estimates of cell densities for the different inhibitory neuron types. We estimate that 20.3% of all neurons in the mouse brain are inhibitory. Among all inhibitory neurons, 18% predominantly express parvalbumin (PV), 16% express somatostatin (SST), 3% express vasoactive intestinal peptide (VIP), and the remainder 63% belong to the residual GABAergic population. We find that our density estimations improve as more literature values are integrated. Our pipeline is extensible, allowing new cell types or data to be integrated as they become available. The data, algorithms, software, and results of our pipeline are publicly available and update the Blue Brain Cell Atlas. This work therefore leverages the research community to collectively converge on the numbers of each cell type in each brain region.

Modeling spinal locomotor circuits for movements in developing zebrafish.

Many spinal circuits dedicated to locomotor control have been identified in the developing zebrafish. How these circuits operate together to generate the various swimming movements during development remains to be clarified. In this study, we iteratively built models of developing zebrafish spinal circuits coupled to simplified musculoskeletal models that reproduce coiling and swimming movements. The neurons of the models were based upon morphologically or genetically identified populations in the developing zebrafish spinal cord. We simulated intact spinal circuits as well as circuits with silenced neurons or altered synaptic transmission to better understand the role of specific spinal neurons. Analysis of firing patterns and phase relationships helped to identify possible mechanisms underlying the locomotor movements of developing zebrafish. Notably, our simulations demonstrated how the site and the operation of rhythm generation could transition between coiling and swimming. The simulations also underlined the importance of contralateral excitation to multiple tail beats. They allowed us to estimate the sensitivity of spinal locomotor networks to motor command amplitude, synaptic weights, length of ascending and descending axons, and firing behavior. These models will serve as valuable tools to test and further understand the operation of spinal circuits for locomotion.

The spinal cord is a column of nerve tissue that connects the brain to the rest of the body in vertebrate animals. Nerve cells in the spinal cord, called neurons, help to control and coordinate the body's movements. As the spinal cord develops, new neurons are born and new connections are made between neurons and muscles, resulting in more coordinated and skillful movements as time goes on. Zebrafish, for example, display body-bending maneuvers called coils within 24 hours of the egg being fertilized. Next, bursts of swimming movements emerge, which are driven by sporadic tail beats. These tail maneuvers become more consistent as the fish develops, and eventually result in smooth movements called beat-and-glide swimming. The groups of spinal cord neurons that appear at each stage of zebrafish development have been characterized, but it remains unclear how newly formed circuits (groups of neurons recently connected to each other) work together to produce swimming maneuvers. To answer this question, Roussel et al. simulated changes in the spinal cord that help zebrafish acquire new swimming movements as they grow. The computer models encoded neural circuits based on cell populations identified in experimental studies, and replicated swimming behaviors that emerge during the first few days of zebrafish development. Simulations tested how specific neural circuits generate the characteristic swimming movements that represent key developmental milestones in zebrafish. The results showed that adding new neurons and more cell-to-cell connections led to increasingly sophisticated swimming maneuvers. As the zebrafish spinal cord matured, the fish were better able to control the pace and duration of their swimming movements. Roussel et al. also identified specific patterns of neural activity linked to particular maneuvers. For example, tail beats switch direction when neurons on one side of the spinal cord excite neurons on the opposite side. This activity, which becomes more rhythmic, also needs to be exquisitely timed to produce and coordinate the right motion. Roussel et al.'s modelling of developmental milestones in growing zebrafish provides insights into how neural networks control movement. The computer models are among the first to accurately reproduce swimming behaviors in developing zebrafish. More experimental data could be added to the models to capture the full range of early zebrafish movements, and to further investigate how maturing spinal cord circuits control swimming. Since zebrafish and mammals have many spinal neurons in common, further research may aid our understanding of movement disorders in humans.

Spinal V1 neurons inhibit motor targets locally and sensory targets distally.

Rostro-caudal coordination of spinal motor output is essential for locomotion. Most spinal interneurons project axons longitudinally to govern locomotor output, yet their connectivity along this axis remains unclear. In this study, we use larval zebrafish to map synaptic outputs of a major inhibitory population, V1 (Eng1+) neurons, which are implicated in dual sensory and motor functions. We find that V1 neurons exhibit long axons extending rostrally and exclusively ipsilaterally for an average of 6 spinal segments; however, they do not connect uniformly with their post-synaptic targets along the entire length of their axon. Locally, V1 neurons inhibit motor neurons (both fast and slow) and other premotor targets, including V2a, V2b, and commissural premotor neurons. In contrast, V1 neurons make robust long-range inhibitory contacts onto a dorsal horn sensory population, the commissural primary ascending neurons (CoPAs). In a computational model of the ipsilateral spinal network, we show that this pattern of short-range V1 inhibition to motor and premotor neurons underlies burst termination, which is critical for coordinated rostro-caudal propagation of the locomotor wave. We conclude that spinal network architecture in the longitudinal axis can vary dramatically, with differentially targeted local and distal connections, yielding important consequences for function.

Spatiotemporal Transition in the Role of Synaptic Inhibition to the Tail Beat Rhythm of Developing Larval Zebrafish.

Significant maturation of swimming in zebrafish (Danio rerio) occurs within the first few days of life when fish transition from coiling movements to burst swimming and then to beat-and-glide swimming. This maturation occurs against a backdrop of numerous developmental changes - neurogenesis, a transition from predominantly electrical to chemical-based neurotransmission, and refinement of intrinsic properties. There is evidence that spinal locomotor circuits undergo fundamental changes as the zebrafish transitions from burst to beat-and-glide swimming. Our electrophysiological recordings confirm that the operation of spinal locomotor circuits becomes increasingly reliant on glycinergic neurotransmission for rhythmogenesis governing the rhythm of tail beats. This transition occurred at the same time that we observed a change in rhythmicity of synaptic inhibition to spinal motoneurons (MNs). When we examined whether the transition from weakly to strongly glycinergic dependent rhythmogenesis occurred at a uniform pace across the length of the spinal cord, we found that this transition occurred earlier at caudal segments than at rostral segments of the spinal cord. Furthermore, while this rhythmogenic transition occurred when fish transition from burst swimming to beat-and-glide swimming, these two transitions were not interdependent. These results suggest that there is a developmental transition in the operation of spinal locomotor circuits that is gradually set in place in the spinal cord in a caudo-rostral temporal sequence.

PLPHP deficiency: clinical, genetic, biochemical, and mechanistic insights.

Biallelic pathogenic variants in PLPBP (formerly called PROSC) have recently been shown to cause a novel form of vitamin B6-dependent epilepsy, the pathophysiological basis of which is poorly understood. When left untreated, the disease can progress to status epilepticus and death in infancy. Here we present 12 previously undescribed patients and six novel pathogenic variants in PLPBP. Suspected clinical diagnoses prior to identification of PLPBP variants included mitochondrial encephalopathy (two patients), folinic acid-responsive epilepsy (one patient) and a movement disorder compatible with AADC deficiency (one patient). The encoded protein, PLPHP is believed to be crucial for B6 homeostasis. We modelled the pathogenicity of the variants and developed a clinical severity scoring system. The most severe phenotypes were associated with variants leading to loss of function of PLPBP or significantly affecting protein stability/PLP-binding. To explore the pathophysiology of this disease further, we developed the first zebrafish model of PLPHP deficiency using CRISPR/Cas9. Our model recapitulates the disease, with plpbp-/- larvae showing behavioural, biochemical, and electrophysiological signs of seizure activity by 10 days post-fertilization and early death by 16 days post-fertilization. Treatment with pyridoxine significantly improved the epileptic phenotype and extended lifespan in plpbp-/- animals. Larvae had disruptions in amino acid metabolism as well as GABA and catecholamine biosynthesis, indicating impairment of PLP-dependent enzymatic activities. Using mass spectrometry, we observed significant B6 vitamer level changes in plpbp-/- zebrafish, patient fibroblasts and PLPHP-deficient HEK293 cells. Additional studies in human cells and yeast provide the first empirical evidence that PLPHP is localized in mitochondria and may play a role in mitochondrial metabolism. These models provide new insights into disease mechanisms and can serve as a platform for drug discovery.

Pyridoxine-Dependent Epilepsy in Zebrafish Caused by Aldh7a1 Deficiency.

Pyridoxine-dependent epilepsy (PDE) is a rare disease characterized by mutations in the lysine degradation gene ALDH7A1 leading to recurrent neonatal seizures, which are uniquely alleviated by high doses of pyridoxine or pyridoxal 5'-phosphate (vitamin B6 vitamers). Despite treatment, neurodevelopmental disabilities are still observed in most PDE patients underlining the need for adjunct therapies. Over 60 years after the initial description of PDE, we report the first animal model for this disease: an aldh7a1-null zebrafish (Danio rerio) displaying deficient lysine metabolism and spontaneous and recurrent seizures in the larval stage (10 days postfertilization). Epileptiform electrographic activity was observed uniquely in mutants as a series of population bursts in tectal recordings. Remarkably, as is the case in human PDE, the seizures show an almost immediate sensitivity to pyridoxine and pyridoxal 5'-phosphate, with a resulting extension of the life span. Lysine supplementation aggravates the phenotype, inducing earlier seizure onset and death. By using mass spectrometry techniques, we further explored the metabolic effect of aldh7a1 knockout. Impaired lysine degradation with accumulation of PDE biomarkers, B6 deficiency, and low γ-aminobutyric acid levels were observed in the aldh7a1-/- larvae, which may play a significant role in the seizure phenotype and PDE pathogenesis. This novel model provides valuable insights into PDE pathophysiology; further research may offer new opportunities for drug discovery to control seizure activity and improve neurodevelopmental outcomes for PDE.