Search

IntroductionInherited retinal diseases (IRDs) affect ∼4.5 million people worldwide. Elusive pathogenic variants in over 280 genes are associated with one or more clinical forms of IRDs. It is necessary to understand the complex interaction among retinal cell types and pathogenic genes by constructing a regulatory network. In this study, we attempt to establish a panoramic expression view of the cooperative work in retinal cells to understand the clinical manifestations and pathogenic bases underlying IRDs.MethodsSingle-cell RNA sequencing (scRNA-seq) data on the retinas from 35 retina samples of 3 species (human, mouse, and zebrafish) including 259,087 cells were adopted to perform a comparative analysis across species. Bioinformatic tools were used to conduct weighted gene co-expression network analysis (WGCNA), single-cell regulatory network analysis, cell–cell communication analysis, and trajectory inference analysis.ResultsThe cross-species comparison revealed shared or species-specific gene expression patterns at single-cell resolution, such as the stathmin family genes, which were highly expressed specifically in zebrafish Müller glias (MGs). Thirteen gene modules were identified, of which nine were associated with retinal cell types, and Gene Ontology (GO) enrichment of module genes was consistent with cell-specific highly expressed genes. Many IRD genes were identified as hub genes and cell-specific regulons. Most IRDs, especially the retinitis pigmentosa (RP) genes, were enriched in rod-specific regulons. Integrated expression and transcription regulatory network genes, such as congenital stationary night blindness (CSNB) genes GRK1, PDE6B, and TRPM1, showed cell-specific expression and transcription characteristics in either rods or bipolar cells (BCs). IRD genes showed evolutionary conservation (GNAT2, PDE6G, and SAG) and divergence (GNAT2, MT-ND4, and PDE6A) along the trajectory of photoreceptors (PRs) among species. In particular, the Leber congenital amaurosis (LCA) gene OTX2 showed high expression at the beginning of the trajectory of both PRs and BCs.ConclusionWe identified molecular pathways and cell types closely connected with IRDs, bridging the gap between gene expression, genetics, and pathogenesis. The IRD genes enriched in cell-specific modules and regulons suggest that these diseases share common etiological bases. Overall, mining of interspecies transcriptome data reveals conserved transcriptomic features of retinas across species and promising applications in both normal retina anatomy and retina pathology.

Summary: DEAH-Box Helicase 38 (DHX38) is a pre-mRNA splicing factor and also a disease-causing gene of autosomal recessive retinitis pigmentosa (arRP). The role of DHX38 in the development and maintenance of the retina remains largely unknown. In this study, by using the dhx38 knockout zebrafish model, we demonstrated that Dhx38 deficiency causes severe differentiation defects and apoptosis of retinal progenitor cells (RPCs) through disrupted mitosis and increased DNA damage. Furthermore, we found a significant accumulation of R-loops in the dhx38-deficient RPCs and human cell lines. Finally, we found that DNA replication stress is the prerequisite for R-loop-induced DNA damage in the DHX38 knockdown cells. Taken together, our study demonstrates a necessary role of DHX38 in the development of retina and reveals a DHX38/R-loop/replication stress/DNA damage regulatory axis that is relatively independent of the known functions of DHX38 in mitosis control.

Abstract Mutations in TUB-like protein 1 (TULP1) are associated with severe early-onset retinal degeneration in humans. However, the pathogenesis remains largely unknown. There are two homologous genes of TULP1 in zebrafish, namely tulp1a and tulp1b. Here, we generated the single knockout (tulp1a −/− and tulp1b −/− ) and double knockout (tulp1-dKO) models in zebrafish. Knockout of tulp1a resulted in the mislocalization of UV cone opsins and the degeneration of UV cones specifically, while knockout of tulp1b resulted in mislocalization of rod opsins and rod-cone degeneration. In the tulp1-dKO zebrafish, mislocalization of opsins was present in all types of photoreceptors, and severe degeneration was observed at a very early age, mimicking the clinical manifestations of TULP1 patients. Photoreceptor cilium length was significantly reduced in the tulp1-dKO retinas. RNA-seq analysis showed that the expression of tektin2 (tekt2), a ciliary and flagellar microtubule structural component, was downregulated in the tulp1-dKO zebrafish. Dual-luciferase reporter assay suggested that Tulp1a and Tulp1b transcriptionally activate the promoter of tekt2. In addition, ferroptosis might be activated in the tulp1-dKO zebrafish, as suggested by the up-regulation of genes related to the ferroptosis pathway, the shrinkage of mitochondria, reduction or disappearance of mitochondria cristae, and the iron and lipid droplet deposition in the retina of tulp1-dKO zebrafish. In conclusion, our study establishes an appropriate zebrafish model for TULP1-associated retinal degeneration and proposes that loss of TULP1 causes defects in cilia structure and opsin trafficking through the downregulation of tekt2, which further increases the death of photoreceptors via ferroptosis. These findings offer insight into the pathogenesis and clinical treatment of early-onset retinal degeneration.

Introduction: As a congenital and genetically related disease, many single nucleotide polymorphisms (SNPs) have been reported to be associated with the risk of HSCR. Our previous research showed that SNP rs2439302 (NRG1) interacted with rs2435357 (RET) to increase the risk of HSCR development. However, the underlying molecular mechanism is still not well understood.Methods: SNP rs2439302 (NRG1) and rs2435357 (RET) were genotyped in 470 HSCR cases. The expression of NRG1 and RET was investigated in the colon of HSCR patients. Knockdown of the NRG1 and RET homologs was performed in zebrafish to investigate their synergistic effect on ENS development. The effect of SNP rs2439302 and rs2435357 polymorphism on neuron proliferation, migration, and differentiation were investigated in SHSY-5Y cells and IPSCs.Results: Significant downregulation of NRG1 and RET expression was noticed in the aganglionic segment of HSCR patients and SHSY-5Y cells with rs2439302 GG/rs2435357 TT genotype. NRG1 and RET double mutants caused the most severe reduction in enteric neuron numbers than NRG1 single mutant or RET single mutant in the hindgut of zebrafish. SHSY-5Y cells and IPSCs with rs2439302 GG/rs2435357 TT genotype exhibited a decreased proliferative, migration, and differentiative capacity. CTCF showed a considerably higher binding ability to SNP rs2439302 CC than GG. NRG1 reduction caused a further decrease in SOX10 expression via the PI3K/Akt pathway, which regulates RET expression by directly binding to rs2435357.Discussion: SNP rs2439302 (NRG1) GG increases the risk of developing HSCR by affecting the binding of transcription factor CTCF and interacting with rs2435357 (RET) to regulate RET expression via the PI3K/Akt/SOX10 pathway.

Neural retina leucine zipper (NRL) is an essential gene for the fate determination and differentiation of the precursor cells into rod photoreceptors in mammals. Mutations in NRL are associated with the autosomal recessive enhanced S-cone syndrome and autosomal dominant retinitis pigmentosa. However, the exact role of Nrl in regulating the development and maintenance of photoreceptors in the zebrafish (Danio rerio), a popular animal model used for retinal degeneration and regeneration studies, has not been fully determined. In this study, we generated an nrl knockout zebrafish model via the CRISPR-Cas9 technology and observed a surprising phenotype characterized by a reduced number, but not the total loss, of rods and over-growth of green cones. We discovered two waves of rod genesis, nrl-dependent and -independent at the embryonic and post-embryonic stages, respectively, in zebrafish by monitoring the rod development. Through bulk and single-cell RNA sequencing, we characterized the gene expression profiles of the whole retina and each retinal cell type from the wild type and nrl knockout zebrafish. The over-growth of green cones and mis-expression of green-cone-specific genes in rods in nrl mutants suggested that there are rod/green-cone bipotent precursors, whose fate choice between rod versus green-cone is controlled by nrl. Besides, we identified the mafba gene as a novel regulator of the nrl-independent rod development, based on the cell-type-specific expression patterns and the retinal phenotype of nrl/mafba double-knockout zebrafish. Gene collinearity analysis revealed the evolutionary origin of mafba and suggested that the function of mafba in rod development is specific to modern fishes. Furthermore, the altered photoreceptor composition and abnormal gene expression in nrl mutants caused progressive retinal degeneration and subsequent regeneration. Accordingly, this study revealed a novel function of the mafba gene in rod development and established a working model for the developmental and regulatory mechanisms regarding the rod and green-cone photoreceptors in zebrafish.

Zebrafish is an excellent model for exploring the development of the inner ear. Its inner ear has similar functions to that of humans, specifically in the maintenance of hearing and balance. Mafba is a component of the Maf transcription factor family. It participates in multiple biological processes, but its role in inner-ear development remains poorly understood. In this study, we constructed a mafba knockout (mafba−/−) zebrafish model using CRISPR/Cas9 technology. The mafba−/− mutant inner ear displayed severe impairments, such as enlarged otocysts, smaller or absent otoliths, and insensitivity to sound stimulation. The proliferation of p63+ epidermal stem cells and dlc+ ionocyte progenitors was inhibited in mafba−/− mutants. Moreover, the results showed that mafba deletion induces the apoptosis of differentiated K+-ATPase-rich (NR) cells and H+-ATPase-rich (HR) cells. The activation of p53 apoptosis and G0/G1 cell cycle arrest resulted from DNA damage in the inner-ear region, providing a mechanism to account for the inner ear deficiencies. The loss of homeostasis resulting from disorders of ionocyte progenitors resulted in structural defects in the inner ear and, consequently, loss of hearing. In conclusion, the present study elucidated the function of ionic channel homeostasis and inner-ear development using a zebrafish Mafba model and clarified the possible physiological roles.

Pre-mRNA splicing is a precise regulated process and is crucial for system development and homeostasis maintenance. Mutations in spliceosomal components have been found in various hematopoietic malignancies (HMs) and have been considered as oncogenic derivers of HMs. However, the role of spliceosomal components in normal and malignant hematopoiesis remains largely unknown. Pre-mRNA processing factor 31 (PRPF31) is a constitutive spliceosomal component, which mutations are associated with autosomal dominant retinitis pigmentosa. PRPF31 was found to be mutated in several HMs, but the function of PRPF31 in normal hematopoiesis has not been explored. In our previous study, we generated a prpf31 knockout (KO) zebrafish line and reported that Prpf31 regulates the survival and differentiation of retinal progenitor cells by modulating the alternative splicing of genes involved in mitosis and DNA repair. In this study, by using the prpf31 KO zebrafish line, we discovered that prpf31 KO zebrafish exhibited severe defects in hematopoietic stem and progenitor cell (HSPC) expansion and its sequentially differentiated lineages. Immunofluorescence results showed that Prpf31-deficient HSPCs underwent malformed mitosis and M phase arrest during HSPC expansion. Transcriptome analysis and experimental validations revealed that Prpf31 deficiency extensively perturbed the alternative splicing of mitosis-related genes. Collectively, our findings elucidate a previously undescribed role for Prpf31 in HSPC expansion, through regulating the alternative splicing of mitosis-related genes. [ABSTRACT FROM AUTHOR]

Background There are many subtypes of major depressive disorder (MDD) . Evidence about the differentiation of them is mostly based on symptomatological characteristics, but rarely by objective biological indicators. Objective To explore the differences in event-related potential (ERP) N400 between patients with different subtypes of depression, providing an objective electrophysiological basis for subtyping depression. Methods Two hundred and nine outpatients and inpatients who met the diagnostic criteria of MDD in Diagnostic and Statistical Manual of Mental Disorders, 5th Edition were selected from First Hospital of Shanxi Medical University from February 2019 to December 2021, including 78 with anxious depression (ASD+MD-) , 46 with melancholic depression (ASD-MD+) , 61 with anxious and melancholic depression (ASD+MD+) , and 24 with non-anxious and melancholic depression (ASD-MD-) subtyped using the 30-item Inventory of Depressive Symptomatology (IDS-30) and 17-item Hamilton Depression Rating Scale (HAMD-17) . Thirty-five community-living health volunteers were recruited as controls at the same time.The factor and total scores of the HAMD-17 and the Repeatable Battery for the Assessment of Neuropsychological Status (RBANS) of all subjects were collected. And the latency and amplitude of N400 component of all subjects were collected by inducing N400 components through understanding ambiguous words at the end of sentences. Results (1) The total score and factor scores of HAMD-17 differed significantly across MDD subgroups and control group (P< 0.05) . The total score of HAMD-17 in each of the four MDD subgroups was much higher than that in control group (P0.05) .There were no significant differences in latency and amplitude of N400 component across different electrode sites (P>0.05) , and there were no interaction between no-depression or each of the four subtypes of MDD and electrode sites (P>0.05) . (4) The N400 amplitude value at Fz electrode was negatively correlated with the cognitive impairment factor score (r=-0.170, P=0.016) .The N400 amplitude value at Cz electrode was positively correlated with the block factor score (r=0.151, P=0.033) .The N400 amplitude valueat Pz electrode was positively correlated with the block factor score (r=0.174, P=0.014) .The N400 amplitude value of the Fz electrode was positively correlated with the immediate memory score (r=0.138, P=0.050) .The N400 latency value at Cz electrode was negatively correlated with the delayed memory score (r=-0.155, P=0.028) . Conclusion Extensive cognitive impairment was found in MDD patients, including impaired immediate memory, speech function, attention and delayed memory. ASD+MD+ and ASD+ MD-patients had more severe symptoms, and ASD+MD+ patients had the most severe cognitive impairment. N400 amplitude value was negatively correlated with cognitive impairment and positively correlated with block factor in MDD patients. But these patients demonstrated no impaired function of speech integration induced by ambiguous words at the end of sentences.