Establishment of a zebrafish inbred strain, M-AB, capable of regular breeding and genetic manipulation

Inbred strains of organisms are genetically highly uniform and thus useful for life science research. We have previously reported the ongoing generation of the zebrafish IM strain from the India (IND) strain through full sib-pair mating for 16 generations. However, the IM fish laid a small number of offspring and had a short lifespan, implying the need for discreet care in breeding. Here, we report the subsequent establishment of IM strain as well as the generation of a new inbred zebrafish strain, Mishima-AB (M-AB). M-AB was derived from the *AB strain by full sib-pair mating for over 20 generations, which fulfills the general criterion for the establishment of an inbred strain. In contrast to the IM case, maintenance of the M-AB strain by sib-pair mating required almost no special handling. Genome sequencing of IM individuals from the 47th generation and M-AB individuals from the 27th generation revealed that SNP-based genomic heterogeneity across whole-genome nucleotides was 0.008% and 0.011%, respectively. These percentages were much lower than those of the parental IND (0.197%) and *AB (0.086%) strains. These results indicate that the genomes of these inbred strains were highly homogenous. We also demonstrated the successful microinjection of antisense morpholinos, CRISPR/Cas9, and foreign genes into M-AB embryos at the 1-cell stage. Overall, we report the establishment of a zebrafish inbred strain, M-AB, which is capable of regular breeding and genetic manipulation. This strain will be useful for the analysis of genetically susceptible phenotypes such as behaviors, microbiome features and drug susceptibility.

Congenital anaemia associated with loss-of-function variants in DNA polymerase epsilon 1

DNA polymerase epsilon (Pol ε), a component of the core replisome, is involved in DNA replication. Although genetic defects of Pol ε have been reported to cause immunodeficiency syndromes, its role in haematopoiesis remains unknown. Here, we identified compound heterozygous variants (p.[Asp1131fs];[Thr1891del]) in POLE, encoding Pol ε catalytic subunit A (POLE1), in siblings with a syndromic form of severe congenital transfusion-dependent anaemia. In contrast to Diamond-Blackfan anaemia, marked reticulocytopenia or marked erythroid hypoplasia was not found. Their bone marrow aspirates during infancy revealed erythroid dysplasia with strongly positive TP53 in immunostaining. Repetitive examinations demonstrated trilineage myelodysplasia within 2 years from birth. They had short stature and facial dysmorphism. HEK293 cell-based expression experiments and analyses of patient-derived induced pluripotent stem cells (iPSCs) disclosed a reduced mRNA level of Asp1131fs-POLE1 and defective nuclear translocation of Thr1891del-POLE1. Analysis of iPSCs showed compensatory mRNA upregulation of the other replisome components and increase of the TP53 protein, both suggesting dysfunction of the replisome. We created Pole-knockout medaka fish and found that heterozygous fishes were viable, but with decreased RBCs. Our observations expand the phenotypic spectrum of the Pol ε defect in humans, additionally providing unique evidence linking Pol ε to haematopoiesis.

A resource of single-cell gene expression profiles in a planarian Dugesia japonica

The freshwater planarian Dugesia japonica maintains an abundant heterogeneous cell population called neoblasts, which include adult pluripotent stem cells. Thus, it is an excellent model organism for stem cell and regeneration research. Recently, many single-cell RNA sequencing (scRNA-seq) databases of several model organisms, including other planarian species, have become publicly available; these are powerful and useful resources to search for gene expression in various tissues and cells. However, the only scRNA-seq dataset for D. japonica has been limited by the number of genes detected. Herein, we collected D. japonica cells, and conducted an scRNA-seq analysis. A novel, automatic, iterative cell clustering strategy produced a dataset of 3,404 cells, which could be classified into 63 cell types based on gene expression profiles. We introduced two examples for utilizing the scRNA-seq dataset in this study using D. japonica. First, the dataset provided results consistent with previous studies as well as novel functionally relevant insights, that is, the expression of DjMTA and DjP2X-A genes in neoblasts that give rise to differentiated cells. Second, we conducted an integrative analysis of the scRNA-seq dataset and time-course bulk RNA-seq of irradiated animals, demonstrating that the dataset can help interpret differentially expressed genes captured via bulk RNA-seq. Using the R package "Seurat" and GSE223927, researchers can easily access and utilize this dataset.

Compensatory Enlargement of the Liver after Proton Beam Therapy for Hepatocellular Carcinoma

PURPOSE/OBJECTIVE(S)

Charged particle therapy (CPT) has been applied as a safe and effective treatment option for hepatocellular carcinoma (HCC). Although most HCC patients have cirrhosis, favorable treatment outcome has been achieved with CPT preserving liver function. After proton beam therapy (PBT) for patients with HCC, the liver volume in the non-irradiated area is often enlarged. Here, we evaluated whether enlargement of the non-irradiated liver affects preserving hepatic function and prognosis in HCC patients treated with PBT.

MATERIALS/METHODS

Among consecutive patients with HCC treated with PBT between April 2011 and July 2017, we retrospectively identified patients who fulfilled the following criteria: (i) receiving PBT to the right hepatic lobe, (ii) the left lateral segment was not irradiated, (iii) no local treatment was performed for liver within 12 months after PBT, and (iv) the albumin-bilirubin (ALBI) score was evaluable at 12 months after PBT. The left lateral segment was defined as the non-irradiated region and measured by contrast-enhanced CT just before and 3 months after PBT. ALBI scores just before and 12 months after PBT were compared to evaluate changes of hepatic function. Overall survival rate was estimated using the Kaplan-Meier method, and differences in survival between subgroups were examined using the log-rank test. The ALBI scores were compared using the Wilcoxon signed-rank test.

RESULTS

We identified 40 patients (male/female = 32/8). The median age at the start date of PBT was 72 (range, 54-87) years. The prescribed dose was 66.0-76.0 Gy (relative biological effectiveness) delivered in 10-38 fractions. The median follow-up was 61 (range, 12-126) months. The 5-year overall survival rates were 79.0% (95% CI: 60.4-100.0%) in the larger enlargement group (n = 16, enlarged volume of non-irradiated region 3 months after PBT ≥75 cm3) and 53.7% (95% CI: 36.0-79.9%) in the smaller enlargement group (n = 24, as above, <75 cm3), respectively (p = 0.21). The median ALBI scores just before and 12 months after PBT were -3.14 (95% CI: -3.22- (-2.53)) and -2.74 (95% CI: -3.07- (-2.42)) in the larger enlargement group (p = 0.09), and -2.91 (95% CI: -3.15- (-2.40)) and -2.59 (95% CI: -2.74- (-2.12)) in the smaller enlargement group (p = 0.006), respectively.

CONCLUSION

Our study suggests that larger enlargement of the non-irradiated liver after PBT is related to well-preserved liver function at 1 year and modestly associated with a favorable prognosis.

Search for a τ^{+}τ^{-} Resonance in e^{+}e^{-}→μ^{+}μ^{-}τ^{+}τ^{-} Events with the Belle II Experiment

We report the first search for a nonstandard-model resonance decaying into τ pairs in e^{+}e^{-}→μ^{+}μ^{-}τ^{+}τ^{-} events in the 3.6-10  GeV/c^{2} mass range. We use a 62.8  fb^{-1} sample of e^{+}e^{-} collisions collected at a center-of-mass energy of 10.58 GeV by the Belle II experiment at the SuperKEKB collider. The analysis probes three different models predicting a spin-1 particle coupling only to the heavier lepton families, a Higgs-like spin-0 particle that couples preferentially to charged leptons (leptophilic scalar), and an axionlike particle, respectively. We observe no evidence for a signal and set exclusion limits at 90% confidence level on the product of cross section and branching fraction into τ pairs, ranging from 0.7 to 24 fb, and on the couplings of these processes. We obtain world-leading constraints on the couplings for the leptophilic scalar model for masses above 6.5  GeV/c^{2} and for the axionlike particle model over the entire mass range.

Activation of XBP1 but not ATF6α rescues heart failure induced by persistent ER stress in medaka fish

The unfolded protein response is triggered in vertebrates by ubiquitously expressed IRE1α/β (although IRE1β is gut-specific in mice), PERK, and ATF6α/β, transmembrane-type sensor proteins in the ER, to cope with ER stress, the accumulation of unfolded and misfolded proteins in the ER. Here, we burdened medaka fish, a vertebrate model organism, with ER stress persistently from fertilization by knocking out the AXER gene encoding an ATP/ADP exchanger in the ER membrane, leading to decreased ATP concentration-mediated impairment of the activity of Hsp70- and Hsp90-type molecular chaperones in the ER lumen. ER stress and apoptosis were evoked from 4 and 6 dpf, respectively, leading to the death of all AXER-KO medaka by 12 dpf because of heart failure (medaka hatch at 7 dpf). Importantly, constitutive activation of IRE1α signaling-but not ATF6α signaling-rescued this heart failure and allowed AXER-KO medaka to survive 3 d longer, likely because of XBP1-mediated transcriptional induction of ER-associated degradation components. Thus, activation of a specific pathway of the unfolded protein response can cure defects in a particular organ.

Observation of e^{+}e^{-}→ωχ_{bJ}(1P) and Search for X_{b}→ωϒ(1S) at sqrt[s] near 10.75 GeV

We study the processes e^{+}e^{-}→ωχ_{bJ}(1P) (J=0, 1, or 2) using samples at center-of-mass energies sqrt[s]=10.701, 10.745, and 10.805 GeV, corresponding to 1.6, 9.8, and 4.7  fb^{-1} of integrated luminosity, respectively. These data were collected with the Belle II detector during special operations of the SuperKEKB collider above the ϒ(4S) resonance. We report the first observation of ωχ_{bJ}(1P) signals at sqrt[s]=10.745  GeV. By combining Belle II data with Belle results at sqrt[s]=10.867  GeV, we find energy dependencies of the Born cross sections for e^{+}e^{-}→ωχ_{b1,b2}(1P) to be consistent with the shape of the ϒ(10753) state. These data indicate that the internal structures of the ϒ(10753) and ϒ(10860) states may differ. Including data at sqrt[s]=10.653  GeV, we also search for the bottomonium equivalent of the X(3872) state decaying into ωϒ(1S). No significant signal is observed for masses between 10.45 and 10.65  GeV/c^{2}.

MED27, SLC6A7, and MPPE1 Variants in a Complex Neurodevelopmental Disorder with Severe Dystonia

BACKGROUND

Despite advances in next generation sequencing technologies, the identification of variants of uncertain significance (VUS) can often hinder definitive diagnosis in patients with complex neurodevelopmental disorders.

OBJECTIVE

The objective of this study was to identify and characterize the underlying cause of disease in a family with two children with severe developmental delay associated with generalized dystonia and episodic status dystonicus, chorea, epilepsy, and cataracts.

METHODS

Candidate genes identified by autozygosity mapping and whole-exome sequencing were characterized using cellular and vertebrate model systems.

RESULTS

Homozygous variants were found in three candidate genes: MED27, SLC6A7, and MPPE1. Although the patients had features of MED27-related disorder, the SLC6A7 and MPPE1 variants were functionally investigated. SLC6A7 variant in vitro overexpression caused decreased proline transport as a result of reduced cell-surface expression, and zebrafish knockdown of slc6a7 exhibited developmental delay and fragile motor neuron morphology that could not be rescued by L-proline transporter-G396S RNA. Lastly, patient fibroblasts displayed reduced cell-surface expression of glycophosphatidylinositol-anchored proteins linked to MPPE1 dysfunction.

CONCLUSIONS

We report a family harboring a homozygous MED27 variant with additional loss-of-function SLC6A7 and MPPE1 gene variants, which potentially contribute to a blended phenotype caused by multilocus pathogenic variants. © 2022 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

DeLTa-Seq: direct-lysate targeted RNA-Seq from crude tissue lysate

BACKGROUND

Quantification of gene expression such as RNA-Seq is a popular approach to study various biological phenomena. Despite the development of RNA-Seq library preparation methods and sequencing platforms in the last decade, RNA extraction remains the most laborious and costly step in RNA-Seq of tissue samples of various organisms. Thus, it is still difficult to examine gene expression in thousands of samples.

RESULTS

Here, we developed Direct-RT buffer in which homogenization of tissue samples and direct-lysate reverse transcription can be conducted without RNA purification. The DTT concentration in Direct-RT buffer prevented RNA degradation but not RT in the lysates of several plant tissues, yeast, and zebrafish larvae. Direct reverse transcription on these lysates in Direct-RT buffer produced comparable amounts of cDNA to those synthesized from purified RNA. To maximize the advantage of the Direct-RT buffer, we integrated Direct-RT and targeted RNA-Seq to develop a cost-effective, high-throughput quantification method for the expressions of hundreds of genes: DeLTa-Seq (Direct-Lysate reverse transcription and Targeted RNA-Seq). The DeLTa-Seq method could drastically improve the efficiency and accuracy of gene expression analysis. DeLTa-Seq analysis of 1056 samples revealed the temperature-dependent effects of jasmonic acid and salicylic acid in Arabidopsis thaliana.

CONCLUSIONS

The DeLTa-Seq method can realize large-scale studies using thousands of animal, plant, and microorganism samples, such as chemical screening, field experiments, and studies focusing on individual variability. In addition, Direct-RT is also beneficial for gene expression analysis in small tissues from which it is difficult to purify enough RNA for the experiments.

The p.Thr395Met missense variant of NFIA found in a patient with intellectual disability is a defective variant

Nuclear factor one A (NFIA) is a transcription factor that regulates the development of the central nervous system. Haploinsufficiency of the NFIA gene causes NFIA-related disorder, which includes brain abnormalities and intellectual disability, with or without urinary tract defects. Intragenic deletions, nonsense variants, frameshift variants, and missense variants in one allele of the NFIA gene have been reported to cause various neurological and urogenital symptoms. Here we report a 10-year-old male patient with developmental delay, coarctation of the aorta, and distinctive facial features. Exome analysis identified a rare de novo heterozygous missense variant p.Thr395Met in NFIA. We employed zebrafish as a model organism in our NFIA analysis and found that nfia^-/- zebrafish initially showed a loss of commissural axons in the brain, and eventually underwent growth retardation resulting in premature death. Impairment of the commissural neurons in nfia^-/- zebrafish embryos could be restored by the expression of wild-type human NFIA protein, but not of mutant human protein harboring the p.Thr395Met substitution, indicating that this variant affects the function of NFIA protein. Taken together, we suggest that the p.Thr395Met allele in the NFIA gene is relevant to the pathogenesis of NFIA-related disorder.

Precise Measurement of the D^{0} and D^{+} Lifetimes at Belle II

We report a measurement of the D^{0} and D^{+} lifetimes using D^{0}→K^{-}π^{+} and D^{+}→K^{-}π^{+}π^{+} decays reconstructed in e^{+}e^{-}→cc[over ¯] data recorded by the Belle II experiment at the SuperKEKB asymmetric-energy e^{+}e^{-} collider. The data, collected at center-of-mass energies at or near the ϒ(4S) resonance, correspond to an integrated luminosity of 72  fb^{-1}. The results, τ(D^{0})=410.5±1.1(stat)±0.8(syst)  fs and τ(D^{+})=1030.4±4.7(stat)±3.1(syst)  fs, are the most precise to date and are consistent with previous determinations.

Intracellular and Intercellular Gene Regulatory Network Inference From Time-Course Individual RNA-Seq

Gene regulatory network (GRN) inference is an effective approach to understand the molecular mechanisms underlying biological events. Generally, GRN inference mainly targets intracellular regulatory relationships such as transcription factors and their associated targets. In multicellular organisms, there are both intracellular and intercellular regulatory mechanisms. Thus, we hypothesize that GRNs inferred from time-course individual (whole embryo) RNA-Seq during development can reveal intercellular regulatory relationships (signaling pathways) underlying the development. Here, we conducted time-course bulk RNA-Seq of individual mouse embryos during early development, followed by pseudo-time analysis and GRN inference. The results demonstrated that GRN inference from RNA-Seq with pseudo-time can be applied for individual bulk RNA-Seq similar to scRNA-Seq. Validation using an experimental-source-based database showed that our approach could significantly infer GRN for all transcription factors in the database. Furthermore, the inferred ligand-related and receptor-related downstream genes were significantly overlapped. Thus, the inferred GRN based on whole organism could include intercellular regulatory relationships, which cannot be inferred from scRNA-Seq based only on gene expression data. Overall, inferring GRN from time-course bulk RNA-Seq is an effective approach to understand the regulatory relationships underlying biological events in multicellular organisms.

Search for B^{+}→K^{+}νν[over ¯] Decays Using an Inclusive Tagging Method at Belle II

A search for the flavor-changing neutral-current decay B^{+}→K^{+}νν[over ¯] is performed at the Belle II experiment at the SuperKEKB asymmetric energy electron-positron collider. The data sample corresponds to an integrated luminosity of 63  fb^{-1} collected at the ϒ(4S) resonance and a sample of 9  fb^{-1} collected at an energy 60 MeV below the resonance. Because the measurable decay signature involves only a single charged kaon, a novel measurement approach is used that exploits not only the properties of the B^{+}→K^{+}νν[over ¯] decay, but also the inclusive properties of the other B meson in the ϒ(4S)→BB[over ¯] event, to suppress the background from other B meson decays and light-quark pair production. This inclusive tagging approach offers a higher signal efficiency compared to previous searches. No significant signal is observed. An upper limit on the branching fraction of B^{+}→K^{+}νν[over ¯] of 4.1×10^{-5} is set at the 90% confidence level.

Direct-TRI: High-throughput RNA-extracting Method for All Stages of Zebrafish Development

Recent popularization of next-generation sequencing enables conducting easy transcriptome analysis. Nevertheless, substantial RNA isolation work prior to RNA sequencing, as well as the high cost involved, still makes the routine use of large-scale transcriptome analysis difficult. For example, conventional phenol-chloroform RNA extraction cannot be easily applied to hundreds of samples. Therefore, we developed Direct-TRI, a new cost-effective and high throughput RNA-extraction method that uses a commercial guanidine-phenol-based RNA extraction reagent and a 96-well silica column plate. We applied Direct-TRI to zebrafish whole larvae and juvenile samples and obtained comparable RNA qualities by several different homogenization methods such as vortexing, manual homogenizing, and freezing/crushing. Direct-TRI enabled the extraction of 192 RNA samples in an hour with a cost of less than a dollar per sample. Direct-TRI is useful for large-scale transcriptome studies, manipulating hundreds of zebrafish individuals, and may be used with other animal samples.

Loss of αklotho causes reduced motor ability and short lifespan in zebrafish

The klotho gene encodes a transmembrane protein αKlotho that interacts with a fibroblast growth factor (FGF) receptor in renal tubular epithelial cells and functions as a co-receptor for FGF23, which is an osteocytes-derived hormone. This bone-to-kidney signal promotes urinary phosphate excretion. Interestingly, αKlotho knockout mice show an accelerated aging and a shortened life span. Similarly, C. elegans lacking the αklotho homologue showed a short life span. However, the physiological basis of aging-related function of αklotho remain unclear. The αklotho-deficient vertebrate animals other than mice have been awaited as an alternative model of premature aging. We here employed zebrafish in our study and revealed that αklotho mutant zebrafish appeared to be normal at 3 months postfertilization (mpf) but eventually underwent premature death by 9 mpf, while normal zebrafish is known to survive for 42 months. We also assessed the motor ability of zebrafish in a forced swimming assay and found that αklotho mutant zebrafish displayed reduced swimming performance before their survival declined. A recent study also reported a similar finding that αklotho-deficient zebrafish exhibited a short life span and reduced spontaneous movements. Taken together, these results suggest that αKlotho mutant zebrafish show premature aging and are useful to investigate aging in vertebrates.

Recurrent NFIA K125E substitution represents a loss-of-function allele: Sensitive in vitro and in vivo assays for nontruncating alleles

Nuclear factor I A (NFIA) is a transcription factor that belongs to the NFI family. Truncating variants or intragenic deletion of the NFIA gene are known to cause the human neurodevelopmental disorder known as NFIA‐related disorder, but no patient heterozygous for a missense mutation has been reported. Here, we document two unrelated patients with typical phenotypic features of the NFIA‐related disorder who shared a missense variant p.Lys125Glu (K125E) in the NFIA gene. Patient 1 was a 6‐year‐old female with global developmental delay, corpus callosum anomaly, macrocephaly, and dysmorphic facial features. Patient 2 was a 14‐month‐old male with corpus callosum anomaly and macrocephaly. By using Drosophila and zebrafish models, we functionally evaluated the effect of the K125E substitution. Ectopic expression of wild‐type human NFIA in Drosophila caused developmental defects such as eye malformation and premature death, while that of human NFIA K125E variant allele did not. nfia‐deficient zebrafish embryos showed defects of midline‐crossing axons in the midbrain/hindbrain boundary. This impairment of commissural neurons was rescued by expression of wild‐type human NFIA, but not by that of mutant variant harboring K125E substitution. In accordance with these in vivo functional analyses, we showed that the K125E mutation impaired the transcriptional regulation of HES1 promoter in cultured cells. Taken together, we concluded that the K125E variant in the NFIA gene is a loss‐of‐function mutation.

Characterization of zebrafish GABAA receptor subunits

γ-Aminobutyric acid (GABA), the major inhibitory neurotransmitter in the central nervous system, exerts its effect through the activation of GABA receptors. GABA_A receptors are ligand-gated chloride channels composed of five subunit proteins. Mammals have 19 different GABA_A receptor subunits (α1–6, β1–3, γ1–3, δ, ε, π, θ, and ρ1–3), the physiological properties of which have been assayed by electrophysiology. However, the evolutionary conservation of the physiological characteristics of diverged GABA_A receptor subunits remains unclear. Zebrafish have 23 subunits (α1, α2a, α2b, α3–5, α6a, α6b, β1–4, γ1–3, δ, π, ζ, ρ1, ρ2a, ρ2b, ρ3a, and ρ3b), but the electrophysiological properties of these subunits have not been explored. In this study, we cloned the coding sequences for zebrafish GABA_A receptor subunits and investigated their expression patterns in larval zebrafish by whole-mount in situ hybridization. We also performed electrophysiological recordings of GABA-evoked currents from Xenopus oocytes injected with one or multiple zebrafish GABA_A receptor subunit cRNAs and calculated the half-maximal effective concentrations (EC50s) for each. Our results revealed the spatial expressions and electrophysiological GABA sensitivities of zebrafish GABA_A receptors, suggesting that the properties of GABA_A receptor subunits are conserved among vertebrates.

Search for Axionlike Particles Produced in e^{+}e^{-} Collisions at Belle II

We present a search for the direct production of a light pseudoscalar a decaying into two photons with the Belle II detector at the SuperKEKB collider. We search for the process e^{+}e^{-}→γa, a→γγ in the mass range 0.2<m_{a}<9.7  GeV/c^{2} using data corresponding to an integrated luminosity of (445±3)  pb^{-1}. Light pseudoscalars interacting predominantly with standard model gauge bosons (so-called axionlike particles or ALPs) are frequently postulated in extensions of the standard model. We find no evidence for ALPs and set 95% confidence level upper limits on the coupling strength g_{aγγ} of ALPs to photons at the level of 10^{-3}  GeV^{-1}. The limits are the most restrictive to date for 0.2<m_{a}<1  GeV/c^{2}.

A Reproducible Protocol to Measure the Critical Swimming Speed of Adult Zebrafish

The quantitative measurement of water flow-induced swimming of fish species using a swimmill is a powerful method to evaluate motor ability of individual fish. Zebrafish is a commonly used vertebrate that enables the study of morphological, physiological and behavioral characteristics associated with genes. We here established a reproducible method that allows to measure the body length and the critical swimming speed of adult zebrafish using a swimmill.

InMeRF: prediction of pathogenicity of missense variants by individual modeling for each amino acid substitution

In predicting the pathogenicity of a nonsynonymous single-nucleotide variant (nsSNV), a radical change in amino acid properties is prone to be classified as being pathogenic. However, not all such nsSNVs are associated with human diseases. We generated random forest (RF) models individually for each amino acid substitution to differentiate pathogenic nsSNVs in the Human Gene Mutation Database and common nsSNVs in dbSNP. We named a set of our models ‘Individual Meta RF’ (InMeRF). Ten-fold cross-validation of InMeRF showed that the areas under the curves (AUCs) of receiver operating characteristic (ROC) and precision–recall curves were on average 0.941 and 0.957, respectively. To compare InMeRF with seven other tools, the eight tools were generated using the same training dataset, and were compared using the same three testing datasets. ROC-AUCs of InMeRF were ranked first in the eight tools. We applied InMeRF to 155 pathogenic and 125 common nsSNVs in seven major genes causing congenital myasthenic syndromes, as well as in VANGL1 causing spina bifida, and found that the sensitivity and specificity of InMeRF were 0.942 and 0.848, respectively. We made the InMeRF web service, and also made genome-wide InMeRF scores available online (https://www.med.nagoya-u.ac.jp/neurogenetics/InMeRF/).

Swimming capability of zebrafish is governed by water temperature, caudal fin length and genetic background

Several zebrafish strains such as AB, Tübingen (TU), Wild India Kolkata (WIK) and Tupfel long fin (TL) have been established for genetic study. Each strain has its morphological and behavioral traits. Motor traits, however, have not been explored in zebrafish strains. We here applied a treadmill for fish (swimmill) and measured swimming capability of adult zebrafish by critical swimming speed, which is the maximum water velocity in which fish can keep swimming. First, we confirmed that swimming capability does not vary between female and male. Second, we found that the appropriate water temperature for swimming was between 16 and 30 °C. Third, our fin clip experiments using long-finned zebrafish revealed that they can exhibit high swimming capability when the caudal fin length was set between 3 and 10 mm, implying that long-finned zebrafish are unfavorable for fast swimming. Finally, we compared swimming capability of several zebrafish strains and demonstrated that WIK fish was significantly less capable of swimming despite that they have short caudal fin (~9 mm). The offspring of WIK fish were less capable of swimming, while hybrids of WIK and TU showed high swimming performance comparable to TU. Thus, lower swimming capability of WIK strain is inheritable as a motor trait.

Mechanical regulation of bone homeostasis through p130Cas-mediated alleviation of NF-κB activity

Mechanical loading plays an important role in bone homeostasis. However, molecular mechanisms behind the mechanical regulation of bone homeostasis are poorly understood. We previously reported p130Cas (Cas) as a key molecule in cellular mechanosensing at focal adhesions. Here, we demonstrate that Cas is distributed in the nucleus and supports mechanical loading-mediated bone homeostasis by alleviating NF-κB activity, which would otherwise prompt inflammatory processes. Mechanical unloading modulates Cas distribution and NF-κB activity in osteocytes, the mechanosensory cells in bones. Cas deficiency in osteocytes increases osteoclastic bone resorption associated with NF-κB-mediated RANKL expression, leading to osteopenia. Upon shear stress application on cultured osteocytes, Cas translocates into the nucleus and down-regulates NF-κB activity. Collectively, fluid shear stress-dependent Cas-mediated alleviation of NF-κB activity supports bone homeostasis. Given the ubiquitous expression of Cas and NF-κB together with systemic distribution of interstitial fluid, the Cas-NF-κB interplay may also underpin regulatory mechanisms in other tissues and organs.

Deleterious de novo variants of X-linked ZC4H2 in females cause a variable phenotype with neurogenic arthrogryposis multiplex congenita

Pathogenic variants in the X-linked gene ZC4H2, which encodes a zinc-finger protein, cause an infrequently described syndromic form of arthrogryposis multiplex congenita (AMC) with central and peripheral nervous system involvement. We present genetic and detailed phenotypic information on 23 newly identified families and simplex cases that include 19 affected females from 18 families and 14 affected males from nine families. Of note, the 15 females with deleterious de novo ZC4H2 variants presented with phenotypes ranging from mild to severe, and their clinical features overlapped with those seen in affected males. By contrast, of the nine carrier females with inherited ZC4H2 missense variants that were deleterious in affected male relatives, four were symptomatic. We also compared clinical phenotypes with previously published cases of both sexes and provide an overview on 48 males and 57 females from 42 families. The spectrum of ZC4H2 defects comprises novel and recurrent mostly inherited missense variants in affected males, and de novo splicing, frameshift, nonsense, and partial ZC4H2 deletions in affected females. Pathogenicity of two newly identified missense variants was further supported by studies in zebrafish. We propose ZC4H2 as a good candidate for early genetic testing of males and females with a clinical suspicion of fetal hypo-/akinesia and/or (neurogenic) AMC.

The AAA+ ATPase/ubiquitin ligase mysterin stabilizes cytoplasmic lipid droplets

Mutations in mysterin cause cerebrovascular moyamoya, but the mechanism of pathogenesis is unknown. Sugihara et al. report that mysterin stabilizes cytoplasmic lipid droplets through the activity of ATPase and ubiquitin ligase. Disease-associated mutations in mysterin impair this process, suggesting a potential link between moyamoya disease and metabolism.

Mysterin, also known as RNF213, is an intracellular protein that forms large toroidal oligomers. Mysterin was originally identified in genetic studies of moyamoya disease (MMD), a rare cerebrovascular disorder of unknown etiology. While mysterin is known to exert ubiquitin ligase and putative mechanical ATPase activities with a RING finger domain and two adjacent AAA+ modules, its biological role is poorly understood. Here, we report that mysterin is targeted to lipid droplets (LDs), ubiquitous organelles specialized for neutral lipid storage, and markedly increases their abundance in cells. This effect was exerted primarily through specific elimination of adipose triglyceride lipase (ATGL) from LDs. The ubiquitin ligase and ATPase activities of mysterin were both important for its proper LD targeting. Notably, MMD-related mutations in the ubiquitin ligase domain of mysterin significantly impaired its fat-stabilizing activity. Our findings identify a unique new regulator of cytoplasmic LDs and suggest a potential link between the pathogenesis of MMD and fat metabolism.

Integrin β1 activity is required for cardiovascular formation in zebrafish

Integrins are transmembrane molecules that facilitate cell-to-cell and cell-to-extracellular matrix (ECM) interactions. Integrin molecules are heterodimers that consist of α- and β-subunits. The integrin β1 gene is widely expressed in vivo and is the major β molecule in many tissues; however, tissue-specific roles of integrin β1 are still elusive. In this study, we investigated integrin β1 function in endothelial cells of zebrafish. An integrin β1b mutant zebrafish exhibited morphological abnormalities in blood vessel formation, cephalic hemorrhage and a decreased responsiveness to tactile stimulation during development. To determine the role of integrin β1b in vascular formation, we developed a Gal4/UAS-mediated conditional inactivation of integrin β1 by expressing the cytoplasmic region of integrin β1 that acts as a dominant-negative (DN) isoform. Expression of integrin β1 DN in endothelial cells induced blood vessel abnormalities as in integrin β1b mutants. These results show that endothelial cells require integrin activity for the formation and/or maintenance of blood vessels in zebrafish. Furthermore, our time-lapse recording visualized the breakpoint of cephalic vessels and the hemorrhage onset. Taken together, our tissue-specific inactivation of integrin β1 in zebrafish is powerful tools for functional analysis of integrin β1 in developing tissues.

Characterization of the Zebrafish Glycine Receptor Family Reveals Insights Into Glycine Receptor Structure Function and Stoichiometry

To study characterization of zebrafish glycine receptors (zGlyRs), we assessed expression and function of five α- and two ß-subunit encoding GlyR in zebrafish. Our qPCR analysis revealed variable expression during development, while in situ hybridizations uncovered expression in the hindbrain and spinal cord; a finding consistent with the reported expression of GlyR subunits in these tissues from other organisms. Electrophysiological recordings using Xenopus oocytes revealed that all five α subunits form homomeric receptors activated by glycine, and inhibited by strychnine and picrotoxin. In contrast, ß subunits only formed functional heteromeric receptors when co-expressed with α subunits. Curiously, the second transmembranes of both ß subunits were found to lack a phenylalanine at the sixth position that is commonly associated with conferring picrotoxin resistance to heteromeric receptors. Consistent with the absence of phenylalanines at the sixth position, heteromeric zGlyRs often lacked significant picrotoxin resistance. Subsequent efforts revealed that resistance to picrotoxin in both zebrafish and human heteromeric GlyRs involves known residues within transmembrane 2, as well as previously unknown residues within transmembrane 3. We also found that a dominant mutation in human GlyRα1 that gives rise to hyperekplexia, and recessive mutations in zebrafish GlyRßb that underlie the bandoneon family of motor mutants, result in reduced receptor function. Lastly, through the use of a concatenated construct we demonstrate that zebrafish heteromeric receptors assemble with a stoichiometry of 3α:2ß. Collectively, our findings have furthered our knowledge regarding the assembly of heteromeric receptors, and the molecular basis of ß subunit-conferred picrotoxin resistance. These results should aid in future investigations of glycinergic signaling in zebrafish and mammals.

Functional Consequences of the Postnatal Switch From Neonatal to Mutant Adult Glycine Receptor α1 Subunits in the Shaky Mouse Model of Startle Disease

Mutations in GlyR α1 or β subunit genes in humans and rodents lead to severe startle disease characterized by rigidity, massive stiffness and excessive startle responses upon unexpected tactile or acoustic stimuli. The recently characterized startle disease mouse mutant shaky carries a missense mutation (Q177K) in the β8-β9 loop within the large extracellular N-terminal domain of the GlyR α1 subunit. This results in a disrupted hydrogen bond network around K177 and faster GlyR decay times. Symptoms in mice start at postnatal day 14 and increase until premature death of homozygous shaky mice around 4–6 weeks after birth. Here we investigate the in vivo functional effects of the Q177K mutation using behavioral analysis coupled to protein biochemistry and functional assays. Western blot analysis revealed GlyR α1 subunit expression in wild-type and shaky animals around postnatal day 7, a week before symptoms in mutant mice become obvious. Before 2 weeks of age, homozygous shaky mice appeared healthy and showed no changes in body weight. However, analysis of gait and hind-limb clasping revealed that motor coordination was already impaired. Motor coordination and the activity pattern at P28 improved significantly upon diazepam treatment, a pharmacotherapy used in human startle disease. To investigate whether functional deficits in glycinergic neurotransmission are present prior to phenotypic onset, we performed whole-cell recordings from hypoglossal motoneurons (HMs) in brain stem slices from wild-type and shaky mice at different postnatal stages. Shaky homozygotes showed a decline in mIPSC amplitude and frequency at P9-P13, progressing to significant reductions in mIPSC amplitude and decay time at P18-24 compared to wild-type littermates. Extrasynaptic GlyRs recorded by bath-application of glycine also revealed reduced current amplitudes in shaky mice compared to wild-type neurons, suggesting that presynaptic GlyR function is also impaired. Thus, a distinct, but behaviorally ineffective impairment of glycinergic synapses precedes the symptoms onset in shaky mice. These findings extend our current knowledge on startle disease in the shaky mouse model in that they demonstrate how the progression of GlyR dysfunction causes, with a delay of about 1 week, the appearance of disease symptoms.

Thumb boutonnière deformity without rheumatoid arthritis or trauma

Boutonnière deformity of the thumb without rheumatoid arthritis or trauma is not widely recognised. This study aimed to investigate its prevalence, relation to sex and age, and identifying factors associated with the extensor mechanism using ultrasonography. We examined 248 thumbs from 127 participants who were asymptomatic for hand disorders and had no history of hand injury. Boutonnière deformity was identified in 20 thumbs of 17 participants and was significantly more prevalent among elderly participants than among young participants. No significant differences in sex or hand dominance were noted. The deformity had a significant effect on range of motion and grip and pinch strengths. The extensor pollicis brevis (EPB) tendon was significantly narrower in affected thumbs than in non-affected thumbs. The prevalence of boutonnière deformity without rheumatoid arthritis or trauma was 13%, and the deformity was associated with advanced age and a narrow EPB tendon.

LEVEL OF EVIDENCE

III.

Relationships between pathologic subjective halitosis, olfactory reference syndrome, and social anxiety in young Japanese women

BACKGROUND

Pathologic subjective halitosis is known as a halitosis complaint without objective confirmation of halitosis by others or by halitometer measurements; it has been reported to be associated with social anxiety disorder. Olfactory reference syndrome is a preoccupation with the false belief that one emits a foul and offensive body odor. Generally, patients with olfactory reference syndrome are concerned with multiple body parts. However, the mouth is known to be the most common source of body odor for those with olfactory reference syndrome, which could imply that the two conditions share similar features. Therefore, we investigated potential causal relationships among pathologic subjective halitosis, olfactory reference syndrome, social anxiety, and preoccupations with body part odors.

METHODS

A total of 1360 female students (mean age 19.6 ± 1.1 years) answered a self-administered questionnaire regarding pathologic subjective halitosis, olfactory reference syndrome, social anxiety, and preoccupation with odors of body parts such as mouth, body, armpits, and feet. The scale for pathologic subjective halitosis followed that developed by Tsunoda et al.; participants were divided into three groups based on their scores (i.e., levels of pathologic subjective halitosis). A Bayesian network was used to analyze causal relationships between pathologic subjective halitosis, olfactory reference syndrome, social anxiety, and preoccupations with body part odors.

RESULTS

We found statistically significant differences in the results for olfactory reference syndrome and social anxiety among the various levels of pathologic subjective halitosis (P < 0.001). Residual analyses indicated that students with severe levels of pathologic subjective halitosis showed greater preoccupations with mouth and body odors (P < 0.05). Bayesian network analysis showed that social anxiety directly influenced pathologic subjective halitosis and olfactory reference syndrome. Preoccupations with mouth and body odors also influenced pathologic subjective halitosis.

CONCLUSIONS

Social anxiety may be a causal factor of pathologic subjective halitosis and olfactory reference syndrome.

Activation-Dependent Rapid Postsynaptic Clustering of Glycine Receptors in Mature Spinal Cord Neurons

Inhibitory synapses are established during development but continue to be generated and modulated in strength in the mature nervous system. In the spinal cord and brainstem, presynaptically released inhibitory neurotransmitter dominantly switches from GABA to glycine during normal development in vivo. While presynaptic mechanisms of the shift of inhibitory neurotransmission are well investigated, the contribution of postsynaptic neurotransmitter receptors to this shift is not fully elucidated. Synaptic clustering of glycine receptors (GlyRs) is regulated by activation-dependent depolarization in early development. However, GlyR activation induces hyperpolarization after the first postnatal week, and little is known whether and how presynaptically released glycine regulates postsynaptic receptors in a depolarization-independent manner in mature developmental stage. Here we developed spinal cord neuronal culture of rodents using chronic strychnine application to investigate whether initial activation of GlyRs in mature stage could change postsynaptic localization of GlyRs. Immunocytochemical analyses demonstrate that chronic blockade of GlyR activation until mature developmental stage resulted in smaller clusters of postsynaptic GlyRs that could be enlarged upon receptor activation for 1 h in the mature stage. Furthermore, live cell-imaging techniques show that GlyR activation decreases its lateral diffusion at synapses, and this phenomenon is dependent on PKC, but neither Ca2+ nor CaMKII activity. These results suggest that the GlyR activation can regulate receptor diffusion and cluster size at inhibitory synapses in mature stage, providing not only new insights into the postsynaptic mechanism of shifting inhibitory neurotransmission but also the inhibitory synaptic plasticity in mature nervous system.

Efficacy and safety of collagenase Clostridium histolyticum injection for Dupuytren's contracture in non-Caucasian Japanese patients (CORD-J Study): the first clinical trial in a non-Caucasian population

To assess the efficacy, safety and pharmacokinetics of 0.58 mg collagenase Clostridium histolyticum injections for the treatment of Dupuytren's contracture in Japanese patients, we conducted a phase III, multicentre, uncontrolled, open-label clinical study in patients with Dupuytren's contracture. Of the 77 patients, 66 achieved clinical success in the primary treated joint (86%; 95% confidence interval: 76% to 93%), confirming the efficacy of collagenase Clostridium histolyticum injections. More improvement was seen in the metacarpophalangeal joints than in the proximal interphalangeal joints (94% versus 73%). The main adverse reaction was a local reaction in the injected hand. No tendon rupture or anaphylactic reactions were seen. The concentrations of collagenase Clostridium histolyticum were below the lower limit of quantification in plasma samples at all time points. As seen in global studies in Caucasian patients, a corrective effect on Dupuytren's contracture and good tolerance were observed in most non-Caucasian (Asian) Japanese patients.

LEVEL OF EVIDENCE

Level 3.

Loss of ift122, a Retrograde Intraflagellar Transport (IFT) Complex Component, Leads to Slow, Progressive Photoreceptor Degeneration Due to Inefficient Opsin Transport

In the retina, aberrant opsin transport from cell bodies to outer segments leads to retinal degenerative diseases such as retinitis pigmentosa. Opsin transport is facilitated by the intraflagellar transport (IFT) system that mediates the bidirectional movement of proteins within cilia. In contrast to functions of the anterograde transport executed by IFT complex B (IFT-B), the precise functions of the retrograde transport mediated by IFT complex A (IFT-A) have not been well studied in photoreceptor cilia. Here, we analyzed developing zebrafish larvae carrying a null mutation in ift122 encoding a component of IFT-A. ift122 mutant larvae show unexpectedly mild phenotypes, compared with those of mutants defective in IFT-B. ift122 mutants exhibit a slow onset of progressive photoreceptor degeneration mainly after 7 days post-fertilization. ift122 mutant larvae also develop cystic kidney but not curly body, both of which are typically observed in various ciliary mutants. ift122 mutants display a loss of cilia in the inner ear hair cells and nasal pit epithelia. Loss of ift122 causes disorganization of outer segment discs. Ectopic accumulation of an IFT-B component, ift88, is observed in the ift122 mutant photoreceptor cilia. In addition, pulse-chase experiments using GFP-opsin fusion proteins revealed that ift122 is required for the efficient transport of opsin and the distal elongation of outer segments. These results show that IFT-A is essential for the efficient transport of outer segment proteins, including opsin, and for the survival of retinal photoreceptor cells, rendering the ift122 mutant a unique model for human retinal degenerative diseases.