Prdm13 forms a feedback loop with Ptf1a and is required for glycinergic amacrine cell genesis in the Xenopus Retina

Temporal Transcriptomic Profiling of the Developing Xenopus laevis Eye

Retinal progenitor cells (RPCs) are a multipotent and highly proliferative population that give rise to all retinal cell types during organogenesis. Defining their molecular signature is a key step towards identifying suitable approaches to treat visual impairments. Here, we performed RNA-sequencing of whole eyes from Xenopus at three embryonic stages and used differential expression analysis to define the transcriptomic profiles of optic tissues containing proliferating and differentiating RPCs during retinogenesis. Gene Ontology and KEGG pathway analyses showed that genes associated with developmental pathways (including Wnt and Hedgehog signaling) were upregulated during the period of active RPC proliferation in early retinal development (Nieuwkoop Faber st. 24 and 27). Developing eyes had dynamic expression profiles and shifted to enrichment for metabolic processes and phototransduction during RPC progeny specification and differentiation (st. 35). Furthermore, conserved adult eye regeneration genes were also expressed during early retinal development including sox2, pax6, nrl, and Notch signaling components. The eye transcriptomic profiles presented here span RPC proliferation to retinogenesis and included regrowth-competent stages. Thus, our dataset provides a rich resource to uncover molecular regulators of RPC activity and will allow future studies to address regulators of RPC proliferation during eye repair and regrowth.

Revisiting the development of cerebellar inhibitory interneurons in the light of single-cell genetic analyses

The present review aims to provide a short update of our understanding of the inhibitory interneurons of the cerebellum. While these cells constitute but a minority of all cerebellar neurons, their functional significance is increasingly being recognized. For one, inhibitory interneurons of the cerebellar cortex are now known to constitute a clearly more diverse group than their traditional grouping as stellate, basket, and Golgi cells suggests, and this diversity is now substantiated by single-cell genetic data. The past decade or so has also provided important information about interneurons in cerebellar nuclei. Significantly, developmental studies have revealed that the specification and formation of cerebellar inhibitory interneurons fundamentally differ from, say, the cortical interneurons, and define a mode of diversification critically dependent on spatiotemporally patterned external signals. Last, but not least, in the past years, dysfunction of cerebellar inhibitory interneurons could also be linked with clinically defined deficits. I hope that this review, however fragmentary, may stimulate interest and help focus research towards understanding the cerebellum.

Gene targeting in amyotrophic lateral sclerosis using causality-based feature selection and machine learning

BACKGROUND

Amyotrophic lateral sclerosis (ALS) is a rare progressive neurodegenerative disease that affects upper and lower motor neurons. As the molecular basis of the disease is still elusive, the development of high-throughput sequencing technologies, combined with data mining techniques and machine learning methods, could provide remarkable results in identifying pathogenetic mechanisms. High dimensionality is a major problem when applying machine learning techniques in biomedical data analysis, since a huge number of features is available for a limited number of samples. The aim of this study was to develop a methodology for training interpretable machine learning models in the classification of ALS and ALS-subtypes samples, using gene expression datasets.

METHODS

We performed dimensionality reduction in gene expression data using a semi-automated preprocessing systematic gene selection procedure using Statistically Equivalent Signature (SES), a causality-based feature selection algorithm, followed by Boosted Regression Trees (XGBoost) and Random Forest to train the machine learning classifiers. The SHapley Additive exPlanations (SHAP values) were used for interpretation of the machine learning classifiers. The methodology was developed and tested using two distinct publicly available ALS RNA-seq datasets. We evaluated the performance of SES as a dimensionality reduction method against: (a) Least Absolute Shrinkage and Selection Operator (LASSO), and (b) Local Outlier Factor (LOF).

RESULTS

The proposed methodology achieved 85.18% accuracy for the classification of cerebellum or frontal cortex samples as C9orf72-related familial ALS, sporadic ALS or healthy samples. Importantly, the genes identified as the most determinative have also been reported as disease-associated in ALS literature. When tested in the evaluation dataset, the methodology achieved 88.89% accuracy for the classification of sporadic ALS motor neuron samples. When LASSO was used as feature selection method instead of SES, the accuracy of the machine learning classifiers ranged from 74.07 to 96.30%, depending on tissue assessed, while LOF underperformed significantly (77.78% accuracy for the classification of pooled cerebellum and frontal cortex samples).

CONCLUSIONS

Using SES, we addressed the challenge of high dimensionality in gene expression data analysis, and we trained accurate machine learning ALS classifiers, specific for the gene expression patterns of different disease subtypes and tissue samples, while identifying disease-associated genes.

Disease-causing mutations in genes encoding transcription factors critical for photoreceptor development

Photoreceptor development of the vertebrate visual system is controlled by a complex transcription regulatory network. OTX2 is expressed in the mitotic retinal progenitor cells (RPCs) and controls photoreceptor genesis. CRX that is activated by OTX2 is expressed in photoreceptor precursors after cell cycle exit. NEUROD1 is also present in photoreceptor precursors that are ready to specify into rod and cone photoreceptor subtypes. NRL is required for the rod fate and regulates downstream rod-specific genes including the orphan nuclear receptor NR2E3 which further activates rod-specific genes and simultaneously represses cone-specific genes. Cone subtype specification is also regulated by the interplay of several transcription factors such as THRB and RXRG. Mutations in these key transcription factors are responsible for ocular defects at birth such as microphthalmia and inherited photoreceptor diseases such as Leber congenital amaurosis (LCA), retinitis pigmentosa (RP) and allied dystrophies. In particular, many mutations are inherited in an autosomal dominant fashion, including the majority of missense mutations in CRX and NRL. In this review, we describe the spectrum of photoreceptor defects that are associated with mutations in the above-mentioned transcription factors, and summarize the current knowledge of molecular mechanisms underlying the pathogenic mutations. At last, we deliberate the outstanding gaps in our understanding of the genotype-phenotype correlations and outline avenues for future research of the treatment strategies.

Multi-omics approach dissects cis-regulatory mechanisms underlying North Carolina macular dystrophy, a retinal enhanceropathy

North Carolina macular dystrophy (NCMD) is a rare autosomal-dominant disease affecting macular development. The disease is caused by non-coding single-nucleotide variants (SNVs) in two hotspot regions near PRDM13 and by duplications in two distinct chromosomal loci, overlapping DNase I hypersensitive sites near either PRDM13 or IRX1. To unravel the mechanisms by which these variants cause disease, we first established a genome-wide multi-omics retinal database, RegRet. Integration of UMI-4C profiles we generated on adult human retina then allowed fine-mapping of the interactions of the PRDM13 and IRX1 promoters and the identification of eighteen candidate cis-regulatory elements (cCREs), the activity of which was investigated by luciferase and Xenopus enhancer assays. Next, luciferase assays showed that the non-coding SNVs located in the two hotspot regions of PRDM13 affect cCRE activity, including two NCMD-associated non-coding SNVs that we identified herein. Interestingly, the cCRE containing one of these SNVs was shown to interact with the PRDM13 promoter, demonstrated in vivo activity in Xenopus, and is active at the developmental stage when progenitor cells of the central retina exit mitosis, suggesting that this region is a PRDM13 enhancer. Finally, mining of single-cell transcriptional data of embryonic and adult retina revealed the highest expression of PRDM13 and IRX1 when amacrine cells start to synapse with retinal ganglion cells, supporting the hypothesis that altered PRDM13 or IRX1 expression impairs interactions between these cells during retinogenesis. Overall, this study provides insight into the cis-regulatory mechanisms of NCMD and supports that this condition is a retinal enhanceropathy.

Ptf1a function and transcriptional cis-regulation, a cornerstone in vertebrate pancreas development

Vertebrate pancreas organogenesis is a stepwise process regulated by a complex network of signaling and transcriptional events, progressively steering the early endoderm toward pancreatic fate. Many crucial players of this process have been identified, including signaling pathways, cis‐regulatory elements, and transcription factors (TFs). Pancreas‐associated transcription factor 1a (PTF1A) is one such TF, crucial for pancreas development. PTF1A mutations result in dramatic pancreatic phenotypes associated with severe complications, such as neonatal diabetes and impaired food digestion due to exocrine pancreatic insufficiency. Here, we present a brief overview of vertebrate pancreas development, centered on Ptf1a function and transcriptional regulation, covering similarities and divergences in three broadly studied organisms: human, mouse and zebrafish.

North Carolina Macular Dystrophy: Phenotypic Variability and Computational Analysis of Disease-Associated Noncoding Variants

Purpose

North Carolina macular dystrophy (NCMD) is an autosomal dominant, congenital disorder affecting the central retina. Here, we report clinical and genetic findings in three families segregating NCMD and use epigenomic datasets from human tissues to gain insights into the effect of NCMD-implicated variants.

Methods

Clinical assessment and genetic testing were performed. Publicly available transcriptomic and epigenomic datasets were analyzed and the activity-by-contact method for scoring enhancer elements and linking them to target genes was used.

Results

A previously described, heterozygous, noncoding variant upstream of the PRDM13 gene was detected in all six affected study participants (chr6:100,040,987G>C [GRCh37/hg19]). Interfamilial and intrafamilial variability were observed; the visual acuity ranged from 0.0 to 1.6 LogMAR and fundoscopic findings ranged from visually insignificant, confluent, drusen-like macular deposits to coloboma-like macular lesions. Variable degrees of peripheral retinal spots (which were easily detected on widefield retinal imaging) were observed in all study subjects. Notably, a 6-year-old patient developed choroidal neovascularization and required treatment with intravitreal bevacizumab injections. Computational analysis of the five single nucleotide variants that have been implicated in NCMD revealed that these noncoding changes lie within two putative enhancer elements; these elements are predicted to interact with PRDM13 in the developing human retina. PRDM13 was found to be expressed in the fetal retina, with greatest expression in the amacrine precursor cell population.

Conclusions

We provide further evidence supporting the role of PRDM13 dysregulation in the pathogenesis of NCMD and highlight the usefulness of widefield retinal imaging in individuals suspected to have this condition.

Single-Cell Transcriptomic Comparison of Human Fetal Retina, hPSC-Derived Retinal Organoids, and Long-Term Retinal Cultures

To study the development of the human retina, we use single-cell RNA sequencing (RNA-seq) at key fetal stages and follow the development of the major cell types as well as populations of transitional cells. We also analyze stem cell (hPSC)-derived retinal organoids; although organoids have a very similar cellular composition at equivalent ages as the fetal retina, there are some differences in gene expression of particular cell types. Moreover, the inner retinal lamination is disrupted at more advanced stages of organoids compared with fetal retina. To determine whether the disorganization in the inner retina is due to the culture conditions, we analyze retinal development in fetal retina maintained under similar conditions. These retinospheres develop for at least 6 months, displaying better inner retinal lamination than retinal organoids. Our single-cell RNA sequencing (scRNA-seq) comparisons of fetal retina, retinal organoids, and retinospheres provide a resource for developing better in vitro models for retinal disease.

Investigation of PRDM10 and PRDM13 Expression in Developing Mouse Embryos by an Optimized PACT-Based Embryo Clearing Method

Recent developments in tissue clearing methods have significantly advanced the three-dimensional analysis of biological structures in whole, intact tissue, providing a greater understanding of spatial relationships and biological circuits. Nonetheless, studies have reported issues with maintaining structural integrity and preventing tissue disintegration, limiting the wide application of these techniques to fragile tissues such as developing embryos. Here, we present an optimized passive tissue clearing technique (PACT)-based embryo clearing method, initial embedding PACT (IMPACT)-Basic, that improves tissue rigidity without compromising optical transparency. We also present IMPACT-Advance, which is specifically optimized for thin slices of mouse embryos past E13.5. We demonstrate proof-of-concept by investigating the expression of two relatively understudied PR domain (PRDM) proteins, PRDM10 and PRDM13, in intact cleared mouse embryos at various stages of development. We observed strong PRDM10 and PRDM13 expression in the developing nervous system and skeletal cartilage, suggesting a functional role for these proteins in these tissues throughout embryogenesis.

Sox2 knockdown in the neonatal retina causes cell fate to switch from amacrine to bipolar

Transcription factor Sox2 is widely recognized for its critical roles in the nervous system, including the neural retina. Here, we aimed to reveal the function of Sox2 in the process of mouse postnatal development. After the suppression of Sox2 at P0, there was an increase number in bipolar cells but a decrease in amacrine cells. Inhibited Sox2 expression also led to decreased visual function. Furthermore, we found a distinctive type of retinal cells expressing the characteristic proteins of both bipolar cells and amacrine cells at P6, which may be an intermediate state in which amacrine cells were transforming into bipolar cells. Transcription factors associated with the development of bipolar cells and amacrine cells also support those changes. Our work indicated that inhibition of Sox2 could change cell fate by affecting transcription factors in the development of bipolar cells and amacrine cells, may provide new directions for the study and treatment of retinal genetic diseases and retinal dysplasia.

Multimodal Imaging and Functional Testing in a North Carolina Macular Disease Family: Toxoplasmosis, Fovea Plana, and Torpedo Maculopathy Are Phenocopies

PURPOSE

To describe multimodal imaging and corresponding functional studies in a newly found family with North Carolina macular dystrophy (NCMD). To our knowledge, this is an original report using OCT angiography to evaluate persons with NCMD.

DESIGN

A descriptive, retrospective study of a family with NCMD.

PARTICIPANTS

A total of 3 participants, representing 3 generations of a single family, each demonstrating a different grade of NCMD, underwent clinical and genetic testing.

METHODS

Diagnostic multimodal imaging and functional testing of the retina included color fundus photography, fundus autofluorescence, intravenous fluorescein angiography, spectral-domain OCT and OCT angiography, multifocal electroretinography, full-field electroretinography, and microperimetry. DNA sequencing was performed using Sanger sequencing techniques.

MAIN OUTCOME MEASURES

Spectral-domain OCT images, fundus photographs, fundus autofluorescence images, fluorescein angiograms, OCT angiograms, multifocal electroretinography images, full-field electroretinography images, and microperimetry maps. Sanger sequencing was performed for molecular diagnosis.

RESULTS

Multimodal imaging helped to demonstrate the nature of the retinal and choroidal lesions in each participant and the extent of visual function. Genetic testing demonstrated the variant 2 point mutation (chromosome 6: 99593111) in the deoxyribonuclease 1 hypersensitivity binding site on chromosome 6q16 causing overexpression of the retinal transcription factor PRDM13.

CONCLUSIONS

NCMD has great phenotypic variability, which can be appreciated only by examining multiple family members. To our knowledge, this is an original report that shows a correlation of functional studies with multimodal imaging. These findings are consistent with NCMD being a developmental abnormality of the macula. All layers of the retina and choroid demonstrate maldevelopment and varying degrees of malfunction. Although PRDM13 is expressed in the amacrine cells, we have yet to identify an abnormality specific to this cellular layer. The retinal vasculature appears to be surprisingly well preserved or intact by OCT angiogram compared with that shown in intravenous fluorescein angiograms. OCT angiograms suggest that foveal hypoplasia is a phenocopy of grade 1 NCMD, torpedo maculopathy is a phenocopy of grade 2 NCMD, and in this single family, congenital toxoplasmosis is a phenocopy of grade 3 NCMD.

Unique noncoding variants upstream of PRDM13 are associated with a spectrum of developmental retinal dystrophies including progressive bifocal chorioretinal atrophy

The autosomal dominant progressive bifocal chorioretinal atrophy (PBCRA) disease locus has been mapped to chromosome 6q14-16.2 that overlaps the North Carolina macular dystrophy (NCMD) locus MCDR1. NCMD is a nonprogressive developmental macular dystrophy, in which variants upstream of PRDM13 have been implicated. Whole genome sequencing was performed to interrogate structural variants (SVs) and single nucleotide variants (SNVs) in eight individuals, six affected individuals from two families with PBCRA, and two individuals from an additional family with a related developmental macular dystrophy. A SNV (chr6:100,046,804T>C), located 7.8 kb upstream of the PRDM13 gene, was shared by all PBCRA-affected individuals in the disease locus. Haplotype analysis suggested that the variant arose independently in the two families. The two affected individuals from Family 3 were screened for rare variants in the PBCRA and NCMD loci. This revealed a de novo variant in the proband, 21 bp from the first SNV (chr6:100,046,783A>C). This study expands the noncoding variant spectrum upstream of PRDM13 and suggests altered spatio-temporal expression of PRDM13 as a candidate disease mechanism in the phenotypically distinct but related conditions, NCMD and PBCRA.

Expression of the inactivating deiodinase, Deiodinase 3, in the pre-metamorphic tadpole retina

Thyroid hormone (TH) orchestrates amphibian metamorphosis. Thus, this developmental phase is often used to study TH-dependent responses in specific tissues. However, TH signaling appears early in development raising the question of the control of TH availability in specific cell types prior to metamorphosis. TH availability is under strict temporal and tissue-specific control by deiodinases. We examined the expression of the TH-inactivating enzyme, deiodinase type 3 (D3), during early retinal development. To this end we created a Xenopus laevis transgenic line expressing GFP from the Xenopus dio3 promoter region (pdio3) and followed pdio3-GFP expression in pre-metamorphic tadpoles. To validate retinal GFP expression in the transgenic line as a function of dio3 promoter activity, we used in situ hybridization to compare endogenous dio3 expression to reporter-driven GFP activity. Retinal expression of dio3 increased during pre-metamorphosis through stages NF41, 45 and 48. Both sets of results show dio3 to have cell-specific, dynamic expression in the pre-metamorphic retina. At stage NF48, dio3 expression co-localised with markers for photoreceptors, rods, Opsin-S cones and bipolar neurons. In contrast, in post-metamorphic juveniles dio3 expression was reduced and spatially confined to certain photoreceptors and amacrine cells. We compared dio3 expression at stages NF41 and NF48 with TH-dependent transcriptional responses using another transgenic reporter line: THbZIP-GFP and by analyzing the expression of T3-regulated genes in distinct TH availability contexts. At stage NF48, the majority of retinal cells expressing dio3 were negative for T3 signaling. Notably, most ganglion cells were virtually both dio3-free and T3-responsive. The results show that dio3 can reduce TH availability at the cellular scale. Further, a reduction in dio3 expression can trigger fine-tuned T3 action in cell-type specific maturation at the right time, as exemplified here in photoreceptor survival in the pre-metamorphic retina.

Prdm13 is required for Ebf3+ amacrine cell formation in the retina

Amacrine interneurons play a critical role in the processing of visual signals within the retina. They are highly diverse, representing 30 or more distinct subtypes. Little is known about how amacrine subtypes acquire their unique gene expression and morphological features. We characterized the gene expression pattern of the zinc-finger transcription factor Prdm13 in the mouse. Consistent with a developmental role, Prdm13 was expressed by Ptf1a+ amacrine and horizontal precursors. Over time, Prdm13 expression diverged from the transiently expressed Ptf1a and marked just a subset of amacrine cells in the adult retina. While heterogeneous, we show that most of these Prdm13+ amacrine cells express the transcription factor Ebf3 and the calcium binding protein calretinin. Loss of Prdm13 did not affect the number of amacrine cells formed during development. However, we observed a modest loss of amacrine cells and increased apoptosis that correlated with the onset timing of Ebf3 expression. Adult Prdm13 loss-of-function mice had 25% fewer amacrine cells, altered calretinin expression, and a lack of Ebf3+ amacrines. Forcing Prdm13 expression in retinal progenitor cells did not significantly increase amacrine cell formation, Ebf3 or calretinin expression, and appeared detrimental to the survival of photoreceptors. Our data show that Prdm13 is not required for amacrine fate as a class, but is essential for the formation of Ebf3+ amacrine cell subtypes. Rather than driving subtype identity, Prdm13 may act by restricting competing fate programs to maintain identity and survival.