Aberrant gene expression yet undiminished retinal ganglion cell genesis in iPSC-derived models of optic nerve hypoplasia.

BACKGROUND: Optic nerve hypoplasia (ONH), the leading congenital cause of permanent blindness, is characterized by a retinal ganglion cell (RGC) deficit at birth. Multifactorial developmental events are hypothesized to underlie ONH and its frequently associated neurologic and endocrine abnormalities; however, environmental influences are unclear and genetic underpinnings are unexplored. This work investigates the genetic contribution to ONH RGC production and gene expression using patient induced pluripotent stem cell (iPSC)-derived retinal organoids (ROs). MATERIALS AND METHODS: iPSCs produced from ONH patients and controls were differentiated to ROs. RGC genesis was assessed using immunofluorescence and flow cytometry. Flow-sorted BRN3+ cells were collected for RNA extraction for RNA-Sequencing. Differential gene expression was assessed using DESeq2 and edgeR. PANTHER was employed to identify statistically over-represented ontologies among the differentially expressed genes (DEGs). DEGs of high interest to ONH were distinguished by assessing function, mutational constraint, and prior identification in ONH, autism and neurodevelopmental disorder (NDD) studies. RESULTS: RGC genesis and survival were similar in ONH and control ROs. Differential expression of 70 genes was identified in both DESeq2 and edgeR analyses, representing a ~ 4-fold higher percentage of DEGs than in randomized study participants. DEGs showed trends towards over-representation of validated NDD genes and ONH exome variant genes. Among the DEGs, RAPGEF4 and DMD had the greatest number of disease-relevant features. CONCLUSIONS: ONH genetic background was not associated with impaired RGC genesis but was associated with DEGs exhibiting disease contribution potential. This constitutes some of the first evidence of a genetic contribution to ONH.

Episodic live imaging of cone photoreceptor maturation in GNAT2-EGFP retinal organoids.

Fluorescent reporter pluripotent stem cell-derived retinal organoids are powerful tools to investigate cell type-specific development and disease phenotypes. When combined with live imaging, they enable direct and repeated observation of cell behaviors within a developing retinal tissue. Here, we generated a human cone photoreceptor reporter line by CRISPR/Cas9 genome editing of WTC11-mTagRFPT-LMNB1 human induced pluripotent stem cells (iPSCs) by inserting enhanced green fluorescent protein (EGFP) coding sequences and a 2A self-cleaving peptide at the N-terminus of guanine nucleotide-binding protein subunit alpha transducin 2 (GNAT2). In retinal organoids generated from these iPSCs, the GNAT2-EGFP alleles robustly and exclusively labeled immature and mature cones. Episodic confocal live imaging of hydrogel immobilized retinal organoids allowed tracking of the morphological maturation of individual cones for >18 weeks and revealed inner segment accumulation of mitochondria and growth at 12.2 µm3 per day from day 126 to day 153. Immobilized GNAT2-EGFP cone reporter organoids provide a valuable tool for investigating human cone development and disease.

Episodic live imaging of cone photoreceptor maturation in GNAT2-EGFP retinal organoids.

Fluorescent reporter pluripotent stem cell (PSC) derived retinal organoids are powerful tools to investigate cell type-specific development and disease phenotypes. When combined with live imaging, they enable direct and repeated observation of cell behaviors within a developing retinal tissue. Here, we generated a human cone photoreceptor reporter line by CRISPR/Cas9 genome editing of WTC11-mTagRFPT-LMNB1 human induced pluripotent stem cells (iPSCs) by inserting enhanced green fluorescent protein (EGFP) coding sequences and a 2A self-cleaving peptide at the N-terminus of Guanine Nucleotide-Binding Protein Subunit Alpha Transducin 2 (GNAT2). In retinal organoids generated from these iPSCs, the GNAT2-EGFP allele robustly and exclusively labeled both immature and mature cones starting at culture day 34. Episodic confocal live imaging of hydrogel immobilized retinal organoids allowed tracking of morphological maturation of individual cones for >18 weeks and revealed inner segment accumulation of mitochondria and growth at 12.2 cubic microns per day from day 126 to day 153. Immobilized GNAT2-EGFP cone reporter organoids provide a valuable tool for investigating human cone development and disease.

Non-synonymous, synonymous, and non-coding nucleotide variants contribute to recurrently altered biological processes during retinoblastoma progression.

Retinoblastomas form in response to biallelic RB1 mutations or MYCN amplification and progress to more aggressive and therapy-resistant phenotypes through accumulation of secondary genomic changes. Progression-related changes include recurrent somatic copy number alterations and typically non-recurrent nucleotide variants, including synonymous and non-coding variants, whose significance has been unclear. To determine if nucleotide variants recurrently affect specific biological processes, we identified altered genes and over-represented variant gene ontologies in 168 exome or whole-genome-sequenced retinoblastomas and 12 tumor-matched cell lines. In addition to RB1 mutations, MYCN amplification, and established retinoblastoma somatic copy number alterations, the analyses revealed enrichment of variant genes related to diverse biological processes including histone monoubiquitination, mRNA processing (P) body assembly, and mitotic sister chromatid segregation and cytokinesis. Importantly, non-coding and synonymous variants increased the enrichment significance of each over-represented biological process term. To assess the effects of such mutations, we examined the consequences of a 3' UTR variant of PCGF3 (a BCOR-binding component of Polycomb repressive complex I), dual 3' UTR variants of CDC14B (a regulator of sister chromatid segregation), and a synonymous variant of DYNC1H1 (a regulator of P-body assembly). One PCGF3 and one of two CDC14B 3' UTR variants impaired gene expression whereas a base-edited DYNC1H1 synonymous variant altered protease sensitivity and stability. Retinoblastoma cell lines retained only ~50% of variants detected in tumors and enriched for new variants affecting p53 signaling. These findings reveal potentially important differences in retinoblastoma cell lines and tumors and implicate synonymous and non-coding variants, along with non-synonymous variants, in retinoblastoma oncogenesis.

Characterization and staging of outer plexiform layer development in human retina and retinal organoids.

The development of the first synapse of the visual system between photoreceptors and bipolar cells in the outer plexiform layer (OPL) of the human retina is crucial for visual processing but poorly understood. By studying the maturation state and spatial organization of photoreceptors, depolarizing bipolar cells and horizontal cells in the human fetal retina, we establish a pseudo-temporal staging system for OPL development that we term OPL-Stages 0 to 4. This was validated through quantification of increasingly precise subcellular localization of bassoon to the OPL with each stage (P<0.0001). By applying these OPL staging criteria to human retinal organoids (HROs) derived from human embryonic and induced pluripotent stem cells, we observed comparable maturation from OPL-Stage 0 at day 100 in culture up to OPL-Stage 3 by day 160. Quantification of presynaptic protein localization confirmed progression from OPL-Stage 0 to 3 (P<0.0001). Overall, this study defines stages of human OPL development through mid-gestation and establishes HROs as a model system that recapitulates key aspects of human photoreceptor-bipolar cell synaptogenesis in vitro.

Structural and Functional Characterization of Human Stem-Cell-Derived Retinal Organoids by Live Imaging.

Purpose: Human pluripotent stem cell (hPSC)-derived retinal organoids are a platform for investigating retinal development, pathophysiology, and cellular therapies. In contrast to histologic analysis in which multiple specimens fixed at different times are used to reconstruct developmental processes, repeated analysis of the same living organoids provides a more direct means to characterize changes. New live imaging modalities can provide insights into retinal organoid structure and metabolic function during in vitro growth. This study employed live tissue imaging to characterize retinal organoid development, including metabolic changes accompanying photoreceptor differentiation. Methods: Live hPSC-derived retinal organoids at different developmental stages were examined for microanatomic organization and metabolic function by phase contrast microscopy, optical coherence tomography (OCT), fluorescence lifetime imaging microscopy (FLIM), and hyperspectral imaging (HSpec). Features were compared to those revealed by histologic staining, immunostaining, and microcomputed tomography (micro-CT) of fixed organoid tissue. Results: We used FLIM and HSpec to detect changes in metabolic activity as organoids differentiated into organized lamellae. FLIM detected increased glycolytic activity and HSpec detected retinol and retinoic acid accumulation in the organoid outer layer, coinciding with photoreceptor genesis. OCT enabled imaging of lamellae formed during organoid maturation. Micro-CT revealed three-dimensional structure, but failed to detect lamellae. Conclusions: Live imaging modalities facilitate real-time and nondestructive imaging of retinal organoids as they organize into lamellar structures. FLIM and HSpec enable rapid detection of lamellar structure and photoreceptor metabolism. Live imaging techniques may aid in the continuous evaluation of retinal organoid development in diverse experimental and cell therapy settings.

A Rapid and Sensitive Next-Generation Sequencing Method to Detect RB1 Mutations Improves Care for Retinoblastoma Patients and Their Families.

Retinoblastoma is a childhood eye malignancy that can lead to the loss of vision, eye(s), and sometimes life. The tumors are initiated by inactivating mutations in both alleles of the tumor-suppressor gene, RB1, or, rarely, by MYCN amplification. Timely identification of a germline RB1 mutation in blood samples or either somatic RB1 mutation or MYCN amplification in tumors is important for effective care and management of retinoblastoma patients and their families. However, current procedures to thoroughly test RB1 mutations are complicated and lengthy. Herein, we report a next-generation sequencing-based method capable of detecting point mutations, small indels, and large deletions or duplications across the entire RB1 gene and amplification of MYCN gene on a single platform. From DNA extraction to clinical interpretation requires only 3 days, enabling early molecular diagnosis of retinoblastoma and optimal treatment outcomes. This method can also detect low-level mosaic mutations in blood samples that can be missed by routine Sanger sequencing. In addition, it can differentiate between RB1 mutation- and MYCN amplification-driven retinoblastomas. This rapid, comprehensive, and sensitive method for detecting RB1 mutations and MYCN amplification can readily identify RB1 mutation carriers and thus improve the management and genetic counseling for retinoblastoma patients and their families.

ChIP-Seq to Analyze the Binding of Replication Proteins to Chromatin.

Chromatin immunoprecipitation (ChIP) is a widely used method to study interactions between proteins and discrete chromosomal loci in vivo. ChIP was originally developed for in vivo analysis of protein associations with candidate DNA sequences known or suspected to bind the protein of interest. The advent of DNA microarrays enabled the unbiased, genome-scale identification of all DNA sequences enriched by ChIP, providing a genomic map of a protein's chromatin binding. This method, termed ChIP-chip, is broadly applicable and has been particularly valuable in DNA replication studies to map potential replication origins in Saccharomyces cerevisiae and other organisms based on the specific association of certain replication proteins with these chromosomal elements, which are distributed throughout the genome. More recently, high-throughput sequencing (HTS) technologies have replaced microarrays as the preferred method for genomic analysis of ChIP experiments, and this combination is termed ChIP-Seq. We present a detailed ChIP-Seq protocol for S. cerevisiae that can be adapted for different HTS platforms and for different organisms. We also outline general schemes for data analysis; however, HTS data analyses usually must be tailored specifically for individual studies, depending on the experimental design, data characteristics, and the genome being analyzed.

ChIP-chip to analyze the binding of replication proteins to chromatin using oligonucleotide DNA microarrays.

Chromatin immunoprecipitation (ChIP) is a widely used method to study the interactions between proteins and discrete chromosomal loci in vivo. Originally, ChIP was developed for analysis of protein associations with DNA sequences known or suspected to bind the protein of interest. The advent of DNA microarrays has enabled the identification of all DNA sequences enriched by ChIP, providing a genomic view of protein binding. This powerful approach, termed ChIP-chip, is broadly applicable and has been particularly valuable in DNA replication studies to map replication origins in Saccharomyces cerevisiae based on the association of replication proteins with these chromosomal elements. We present a detailed ChIP-chip protocol for S. cerevisiae that uses oligonucleotide DNA microarrays printed on polylysine-coated glass slides and can also be easily adapted for commercially available high-density tiling microarrays from NimbleGen. We also outline general protocols for data analysis; however, microarray data analyses usually must be tailored specifically for individual studies, depending on experimental design, microarray format, and data quality.

Rad53 regulates replication fork restart after DNA damage in Saccharomyces cerevisiae.

Replication fork stalling at a DNA lesion generates a damage signal that activates the Rad53 kinase, which plays a vital role in survival by stabilizing stalled replication forks. However, evidence that Rad53 directly modulates the activity of replication forks has been lacking, and the nature of fork stabilization has remained unclear. Recently, cells lacking the Psy2-Pph3 phosphatase were shown to be defective in dephosphorylation of Rad53 as well as replication fork restart after DNA damage, suggesting a mechanistic link between Rad53 deactivation and fork restart. To test this possibility we examined the progression of replication forks in methyl-methanesulfonate (MMS)-damaged cells, under different conditions of Rad53 activity. Hyperactivity of Rad53 in pph3Delta cells slows fork progression in MMS, whereas deactivation of Rad53, through expression of dominant-negative Rad53-KD, is sufficient to allow fork restart during recovery. Furthermore, combined deletion of PPH3 and PTC2, a second, unrelated Rad53 phosphatase, results in complete replication fork arrest and lethality in MMS, demonstrating that Rad53 deactivation is a key mechanism controlling fork restart. We propose a model for regulation of replication fork progression through damaged DNA involving a cycle of Rad53 activation and deactivation that coordinates replication restart with DNA repair.

Genome-wide mapping of ORC and Mcm2p binding sites on tiling arrays and identification of essential ARS consensus sequences in S. cerevisiae.

BACKGROUND: Eukaryotic replication origins exhibit different initiation efficiencies and activation times within S-phase. Although local chromatin structure and function influences origin activity, the exact mechanisms remain poorly understood. A key to understanding the exact features of chromatin that impinge on replication origin function is to define the precise locations of the DNA sequences that control origin function. In S. cerevisiae, Autonomously Replicating Sequences (ARSs) contain a consensus sequence (ACS) that binds the Origin Recognition Complex (ORC) and is essential for origin function. However, an ACS is not sufficient for origin function and the majority of ACS matches do not function as ORC binding sites, complicating the specific identification of these sites. RESULTS: To identify essential origin sequences genome-wide, we utilized a tiled oligonucleotide array (NimbleGen) to map the ORC and Mcm2p binding sites at high resolution. These binding sites define a set of potential Autonomously Replicating Sequences (ARSs), which we term nimARSs. The nimARS set comprises 529 ORC and/or Mcm2p binding sites, which includes 95% of known ARSs, and experimental verification demonstrates that 94% are functional. The resolution of the analysis facilitated identification of potential ACSs (nimACSs) within 370 nimARSs. Cross-validation shows that the nimACS predictions include 58% of known ACSs, and experimental verification indicates that 82% are essential for ARS activity. CONCLUSION: These findings provide the most comprehensive, accurate, and detailed mapping of ORC binding sites to date, adding to the emerging picture of the chromatin organization of the budding yeast genome.

Diminished S-phase cyclin-dependent kinase function elicits vital Rad53-dependent checkpoint responses in Saccharomyces cerevisiae.

Cyclin-dependent kinase (CDK) is required for the initiation of chromosomal DNA replication in eukaryotes. In Saccharomyces cerevisiae, the Clb5 and Clb6 cyclins activate Cdk1 and drive replication origin firing. Deletion of CLB5 reduces initiation of DNA synthesis from late-firing origins. We have examined whether checkpoints are activated by loss of Clb5 function and whether checkpoints are responsible for the DNA replication defects associated with loss of Clb5 function. We present evidence for activation of Rad53 and Ddc2 functions with characteristics suggesting the presence of DNA damage. Deficient late origin firing in clb5Delta cells is not due to checkpoint regulation, but instead, directly reflects the decreased abundance of S-phase CDK, as Clb6 activates late origins when its dosage is increased. Moreover, the viability of clb5Delta cells depends on Rad53. Activation of Rad53 by either Mrc1 or Rad9 contributes to the survival of clb5Delta cells, suggesting that both DNA replication and damage pathways are responsive to the decreased origin usage. These results suggest that reduced origin usage leads to stress or DNA damage at replication forks, necessitating the function of Rad53 in fork stabilization. Consistent with the notion that decreased S-CDK function creates stress at replication forks, deletion of RRM3 helicase, which facilitates replisome progression, greatly diminished the growth of clb5Delta cells. Together, our findings indicate that deregulation of S-CDK function has the potential to exacerbate genomic instability by reducing replication origin usage.

The Rpd3-Sin3 histone deacetylase regulates replication timing and enables intra-S origin control in Saccharomyces cerevisiae.

The replication of eukaryotic genomes follows a temporally staged program, in which late origin firing often occurs within domains of altered chromatin structure(s) and silenced genes. Histone deacetylation functions in gene silencing in some late-replicating regions, prompting an investigation of the role of histone deacetylation in replication timing control in Saccharomyces cerevisiae. Deletion of the histone deacetylase Rpd3 or its interacting partner Sin3 caused early activation of late origins at internal chromosomal loci but did not alter the initiation timing of early origins or a late-firing, telomere-proximal origin. By delaying initiation relative to the earliest origins, Rpd3 enables regulation of late origins by the intra-S replication checkpoint. RPD3 deletion suppresses the slow S phase of clb5Delta cells by enabling late origins to fire earlier, suggesting that Rpd3 modulates the initiation timing of many origins throughout the genome. Examination of factors such as Ume6 that function together with Rpd3 in transcriptional repression indicates that Rpd3 regulates origin initiation timing independently of its role in transcriptional repression. This supports growing evidence that for much of the S. cerevisiae genome transcription and replication timing are not linked.