Mapping the transcriptome: Realizing the full potential of spatial data analysis.

RNA sequencing in situ allows for whole-transcriptome characterization at high resolution, while retaining spatial information. These data present an analytical challenge for bioinformatics-how to leverage spatial information effectively? Properties of data with a spatial dimension require special handling, which necessitate a different set of statistical and inferential considerations when compared to non-spatial data. The geographical sciences primarily use spatial data and have developed methods to analye them. Here we discuss the challenges associated with spatial analysis and examine how we can take advantage of practice from the geographical sciences to realize the full potential of spatial information in transcriptomic datasets.

Differential IRF8 Transcription Factor Requirement Defines Two Pathways of Dendritic Cell Development in Humans.

The formation of mammalian dendritic cells (DCs) is controlled by multiple hematopoietic transcription factors, including IRF8. Loss of IRF8 exerts a differential effect on DC subsets, including plasmacytoid DCs (pDCs) and the classical DC lineages cDC1 and cDC2. In humans, cDC2-related subsets have been described including AXL+SIGLEC6+ pre-DC, DC2 and DC3. The origin of this heterogeneity is unknown. Using high-dimensional analysis, in vitro differentiation, and an allelic series of human IRF8 deficiency, we demonstrated that cDC2 (CD1c+DC) heterogeneity originates from two distinct pathways of development. The lymphoid-primed IRF8hi pathway, marked by CD123 and BTLA, carried pDC, cDC1, and DC2 trajectories, while the common myeloid IRF8lo pathway, expressing SIRPA, formed DC3s and monocytes. We traced distinct trajectories through the granulocyte-macrophage progenitor (GMP) compartment showing that AXL+SIGLEC6+ pre-DCs mapped exclusively to the DC2 pathway. In keeping with their lower requirement for IRF8, DC3s expand to replace DC2s in human partial IRF8 deficiency.

Blood and immune development in human fetal bone marrow and Down syndrome.

Haematopoiesis in the bone marrow (BM) maintains blood and immune cell production throughout postnatal life. Haematopoiesis first emerges in human BM at 11-12 weeks after conception1,2, yet almost nothing is known about how fetal BM (FBM) evolves to meet the highly specialized needs of the fetus and newborn. Here we detail the development of FBM, including stroma, using multi-omic assessment of mRNA and multiplexed protein epitope expression. We find that the full blood and immune cell repertoire is established in FBM in a short time window of 6-7 weeks early in the second trimester. FBM promotes rapid and extensive diversification of myeloid cells, with granulocytes, eosinophils and dendritic cell subsets emerging for the first time. The substantial expansion of B lymphocytes in FBM contrasts with fetal liver at the same gestational age. Haematopoietic progenitors from fetal liver, FBM and cord blood exhibit transcriptional and functional differences that contribute to tissue-specific identity and cellular diversification. Endothelial cell types form distinct vascular structures that we show are regionally compartmentalized within FBM. Finally, we reveal selective disruption of B lymphocyte, erythroid and myeloid development owing to a cell-intrinsic differentiation bias as well as extrinsic regulation through an altered microenvironment in Down syndrome (trisomy 21).

Spatial transcriptomics reveals novel genes during the remodelling of the embryonic human arterial valves.

Abnormalities of the arterial valves, including bicuspid aortic valve (BAV) are amongst the most common congenital defects and are a significant cause of morbidity as well as predisposition to disease in later life. Despite this, and compounded by their small size and relative inaccessibility, there is still much to understand about how the arterial valves form and remodel during embryogenesis, both at the morphological and genetic level. Here we set out to address this in human embryos, using Spatial Transcriptomics (ST). We show that ST can be used to investigate the transcriptome of the developing arterial valves, circumventing the problems of accurately dissecting out these tiny structures from the developing embryo. We show that the transcriptome of CS16 and CS19 arterial valves overlap considerably, despite being several days apart in terms of human gestation, and that expression data confirm that the great majority of the most differentially expressed genes are valve-specific. Moreover, we show that the transcriptome of the human arterial valves overlaps with that of mouse atrioventricular valves from a range of gestations, validating our dataset but also highlighting novel genes, including four that are not found in the mouse genome and have not previously been linked to valve development. Importantly, our data suggests that valve transcriptomes are under-represented when using commonly used databases to filter for genes important in cardiac development; this means that causative variants in valve-related genes may be excluded during filtering for genomic data analyses for, for example, BAV. Finally, we highlight "novel" pathways that likely play important roles in arterial valve development, showing that mouse knockouts of RBP1 have arterial valve defects. Thus, this study has confirmed the utility of ST for studies of the developing heart valves and broadens our knowledge of the genes and signalling pathways important in human valve development.

Single-cell RNA sequencing reveals transcriptional changes of human choroidal and retinal pigment epithelium cells during fetal development, in healthy adult and intermediate age-related macular degeneration.

Age-related macular degeneration (AMD) is the most prevalent cause of blindness in the developed world. Vision loss in the advanced stages of the disease is caused by atrophy of retinal photoreceptors, overlying retinal pigment epithelium (RPE) and choroidal endothelial cells. The molecular events that underline the development of these cell types from in utero to adult as well as the progression to intermediate and advanced stages AMD are not yet fully understood. We performed single-cell RNA-sequencing (RNA-Seq) of human fetal and adult RPE-choroidal tissues, profiling in detail all the cell types and elucidating cell type-specific proliferation, differentiation and immunomodulation events that occur up to midgestation. Our data demonstrate that progression from the fetal to adult state is characterized by an increase in expression of genes involved in the oxidative stress response and detoxification from heavy metals, suggesting a better defence against oxidative stress in the adult RPE-choroid tissue. Single-cell comparative transcriptional analysis between a patient with intermediate AMD and an unaffected subject revealed a reduction in the number of RPE cells and melanocytes in the macular region of the AMD patient. Together these findings may suggest a macular loss of RPE cells and melanocytes in the AMD patients, but given the complex processing of tissues required for single-cell RNA-Seq that is prone to technical artefacts, these findings need to be validated by additional techniques in a larger number of AMD patients and controls.

Incorporating microglia-like cells in human induced pluripotent stem cell-derived retinal organoids.

Microglia are the primary resident immune cells in the retina. They regulate neuronal survival and synaptic pruning making them essential for normal development. Following injury, they mediate adaptive responses and under pathological conditions they can trigger neurodegeneration exacerbating the effect of a disease. Retinal organoids derived from human induced pluripotent stem cells (hiPSCs) are increasingly being used for a range of applications, including disease modelling, development of new therapies and in the study of retinogenesis. Despite many similarities to the retinas developed in vivo, they lack some key physiological features, including immune cells. We engineered an hiPSC co-culture system containing retinal organoids and microglia-like (iMG) cells and tested their retinal invasion capacity and function. We incorporated iMG into retinal organoids at 13 weeks and tested their effect on function and development at 15 and 22 weeks of differentiation. Our key findings showed that iMG cells were able to respond to endotoxin challenge in monocultures and when co-cultured with the organoids. We show that retinal organoids developed normally and retained their ability to generate spiking activity in response to light. Thus, this new co-culture immunocompetent in vitro retinal model provides a platform with greater relevance to the in vivo human retina.

Msx1 haploinsufficiency modifies the Pax9-deficient cardiovascular phenotype.

BACKGROUND: Successful embryogenesis relies on the coordinated interaction between genes and tissues. The transcription factors Pax9 and Msx1 genetically interact during mouse craniofacial morphogenesis, and mice deficient for either gene display abnormal tooth and palate development. Pax9 is expressed specifically in the pharyngeal endoderm at mid-embryogenesis, and mice deficient for Pax9 on a C57Bl/6 genetic background also have cardiovascular defects affecting the outflow tract and aortic arch arteries giving double-outlet right ventricle, absent common carotid arteries and interruption of the aortic arch. RESULTS: In this study we have investigated both the effect of a different genetic background and Msx1 haploinsufficiency on the presentation of the Pax9-deficient cardiovascular phenotype. Compared to mice on a C57Bl/6 background, congenic CD1-Pax9-/- mice displayed a significantly reduced incidence of outflow tract defects but aortic arch defects were unchanged. Pax9-/- mice with Msx1 haploinsufficiency, however, have a reduced incidence of interrupted aortic arch, but more cases with cervical origins of the right subclavian artery and aortic arch, than seen in Pax9-/- mice. This alteration in arch artery defects was accompanied by a rescue in third pharyngeal arch neural crest cell migration and smooth muscle cell coverage of the third pharyngeal arch arteries. Although this change in phenotype could theoretically be compatible with post-natal survival, using tissue-specific inactivation of Pax9 to maintain correct palate development whilst inducing the cardiovascular defects was unable to prevent postnatal death in the mutant mice. Hyoid bone and thyroid cartilage formation were abnormal in Pax9-/- mice. CONCLUSIONS: Msx1 haploinsufficiency mitigates the arch artery defects in Pax9-/- mice, potentially by maintaining the survival of the 3rd arch artery through unimpaired migration of neural crest cells to the third pharyngeal arches. With the neural crest cell derived hyoid bone and thyroid cartilage also being defective in Pax9-/- mice, we speculate that the pharyngeal endoderm is a key signalling centre that impacts on neural crest cell behaviour highlighting the ability of cells in different tissues to act synergistically or antagonistically during embryo development.

Deconstructing Retinal Organoids: Single Cell RNA-Seq Reveals the Cellular Components of Human Pluripotent Stem Cell-Derived Retina.

The rapid improvements in single cell sequencing technologies and analyses afford greater scope for dissecting organoid cultures composed of multiple cell types and create an opportunity to interrogate these models to understand tissue biology, cellular behavior and interactions. To this end, retinal organoids generated from human embryonic stem cells (hESCs) were analyzed by single cell RNA-sequencing (scRNA-Seq) at three time points of differentiation. Combinatorial data from all time points revealed the presence of nine clusters, five of which corresponded to key retinal cell types: retinal pigment epithelium (RPE), retinal ganglion cells (RGCs), cone and rod photoreceptors, and Müller glia. The remaining four clusters expressed genes typical of mitotic cells, extracellular matrix components and those involved in homeostasis. The cell clustering analysis revealed the decreasing presence of mitotic cells and RGCs, formation of a distinct RPE cluster, the emergence of cone and rod photoreceptors from photoreceptor precursors, and an increasing number of Müller glia cells over time. Pseudo-time analysis resembled the order of cell birth during retinal development, with the mitotic cluster commencing the trajectory and the large majority of Müller glia completing the time line. Together, these data demonstrate the feasibility and potential of scRNA-Seq to dissect the inherent complexity of retinal organoids and the orderly birth of key retinal cell types. Stem Cells 2019;37:593-598.

CRX Expression in Pluripotent Stem Cell-Derived Photoreceptors Marks a Transplantable Subpopulation of Early Cones.

Death of photoreceptors is a common cause of age-related and inherited retinal dystrophies, and thus their replenishment from renewable stem cell sources is a highly desirable therapeutic goal. Human pluripotent stem cells provide a useful cell source in view of their limitless self-renewal capacity and potential to not only differentiate into cells of the retina but also self-organize into tissue with structure akin to the human retina as part of three-dimensional retinal organoids. Photoreceptor precursors have been isolated from differentiating human pluripotent stem cells through application of cell surface markers or fluorescent reporter approaches and shown to have a similar transcriptome to fetal photoreceptors. In this study, we investigated the transcriptional profile of CRX-expressing photoreceptor precursors derived from human pluripotent stem cells and their engraftment capacity in an animal model of retinitis pigmentosa (Pde6brd1), which is characterized by rapid photoreceptor degeneration. Single cell RNA-Seq analysis revealed the presence of a dominant cell cluster comprising 72% of the cells, which displayed the hallmarks of early cone photoreceptor expression. When transplanted subretinally into the Pde6brd1 mice, the CRX+ cells settled next to the inner nuclear layer and made connections with the inner neurons of the host retina, and approximately one-third of them expressed the pan cone marker, Arrestin 3, indicating further maturation upon integration into the host retina. Together, our data provide valuable molecular insights into the transcriptional profile of human pluripotent stem cells-derived CRX+ photoreceptor precursors and indicate their usefulness as a source of transplantable cone photoreceptors. Stem Cells 2019;37:609-622.

CD200 Expression Marks a Population of Quiescent Limbal Epithelial Stem Cells with Holoclone Forming Ability.

One of the main challenges in limbal stem cell (LSC) biology and transplantation is the lack of definitive cell surface markers which can be used to identify and enrich viable LSCs. In this study, expression of 361 cell surface proteins was assessed in ex vivo expanded limbal epithelial cells. One marker, CD200 was selected for further characterization based on expression in a small subset of limbal epithelial cells (2.25% ± 0.69%) and reduced expression through consecutive passaging and calcium induced differentiation. CD200 was localized to a small population of cells at the basal layer of the human and mouse limbal epithelium. CD200+ cells were slow cycling and contained the majority of side population (SP) and all the holoclone forming progenitors. CD200+ cells displayed higher expression of LSCs markers including PAX6, WNT7A, CDH3, CK14, CK15, and ABCB5 and lower expression of Ki67 when compared to CD200- . Downregulation of CD200 abrogated the ability of limbal epithelial cells to form holoclones, suggesting an important function for CD200 in the maintenance and/or self-renewal of LSCs. A second marker, CD109, which was expressed in 56.29% ± 13.96% of limbal epithelial cells, was also found to co-localize with ΔNp63 in both human and mouse cornea, albeit more abundantly than CD200. CD109 expression decreased slowly through calcium induced cell differentiation and CD109+ cells were characterized by higher expression of Ki67, when compared to CD109- subpopulation. Together our data suggest that CD200 expression marks a quiescent population of LSCs with holoclone forming potential, while CD109 expression is associated with a proliferative progenitor phenotype. Stem Cells 2018;36:1723-1735.

Haplogroup Context is Less Important in the Penetrance of Mitochondrial DNA Complex I Mutations Compared to mt-tRNA Mutations.

Mitochondrial diseases are a highly complex, heterogeneous group of disorders. Mitochondrial DNA variants that are linked to disease can exhibit variable expression and penetrance. This has an implication for mitochondrial diagnostics as variants that cause disease in one individual may not in another. It has been suggested that the sequence context in which a variant arises could influence the genotype-phenotype relationship. However, the consequence of sequence variation between different haplogroups on the expression of disease is not well understood. European haplogroups are the most widely studied. To ensure accurate diagnostics for patients globally, we first need to understand how, if at all, the sequence context in which a variant arises contributes to the manifestion of disease. To help us understand this, we used 2752 sequences from 33 non-human species that do not have disease. We searched for variants in the seven complex I genes that are associated with disease in humans. Our findings indicate that only three reported pathogenic complex I variants have arisen in these species. More importantly, only one of these, m.3308T>C, has arisen with its associated amino acid change in the studied non-human species. With the status of m.3308T>C as a disease causing variant being a matter of debate. This is a stark contrast to previous findings in the mitochondrial tRNA genes and suggests that sequence context may be less important in the complex I genes. This information will help us improve the identification and diagnosis of mitochondrial DNA variants in non-European populations.

Deciphering the spatiotemporal transcriptional and chromatin accessibility of human retinal organoid development at the single-cell level.

Molecular information on the early stages of human retinal development remains scarce due to limitations in obtaining early human eye samples. Pluripotent stem cell-derived retinal organoids (ROs) provide an unprecedented opportunity for studying early retinogenesis. Using a combination of single cell RNA-seq and spatial transcriptomics we present for the first-time a single cell spatiotemporal transcriptome of RO development. Our data demonstrate that ROs recapitulate key events of retinogenesis including optic vesicle/cup formation, presence of a putative ciliary margin zone, emergence of retinal progenitor cells and their orderly differentiation to retinal neurons. Combining the scRNA- with scATAC-seq data, we were able to reveal cell-type specific transcription factor binding motifs on accessible chromatin at each stage of organoid development, and to show that chromatin accessibility is highly correlated to the developing human retina, but with some differences in the temporal emergence and abundance of some of the retinal neurons.

Single-cell analyses reveal transient retinal progenitor cells in the ciliary margin of developing human retina.

The emergence of retinal progenitor cells and differentiation to various retinal cell types represent fundamental processes during retinal development. Herein, we provide a comprehensive single cell characterisation of transcriptional and chromatin accessibility changes that underline retinal progenitor cell specification and differentiation over the course of human retinal development up to midgestation. Our lineage trajectory data demonstrate the presence of early retinal progenitors, which transit to late, and further to transient neurogenic progenitors, that give rise to all the retinal neurons. Combining single cell RNA-Seq with spatial transcriptomics of early eye samples, we demonstrate the transient presence of early retinal progenitors in the ciliary margin zone with decreasing occurrence from 8 post-conception week of human development. In retinal progenitor cells, we identified a significant enrichment for transcriptional enhanced associate domain transcription factor binding motifs, which when inhibited led to loss of cycling progenitors and retinal identity in pluripotent stem cell derived organoids.

PRPF8-mediated dysregulation of hBrr2 helicase disrupts human spliceosome kinetics and 5´-splice-site selection causing tissue-specific defects.

The carboxy-terminus of the spliceosomal protein PRPF8, which regulates the RNA helicase Brr2, is a hotspot for mutations causing retinitis pigmentosa-type 13, with unclear role in human splicing and tissue-specificity mechanism. We used patient induced pluripotent stem cells-derived cells, carrying the heterozygous PRPF8 c.6926 A > C (p.H2309P) mutation to demonstrate retinal-specific endophenotypes comprising photoreceptor loss, apical-basal polarity and ciliary defects. Comprehensive molecular, transcriptomic, and proteomic analyses revealed a role of the PRPF8/Brr2 regulation in 5'-splice site (5'SS) selection by spliceosomes, for which disruption impaired alternative splicing and weak/suboptimal 5'SS selection, and enhanced cryptic splicing, predominantly in ciliary and retinal-specific transcripts. Altered splicing efficiency, nuclear speckles organisation, and PRPF8 interaction with U6 snRNA, caused accumulation of active spliceosomes and poly(A)+ mRNAs in unique splicing clusters located at the nuclear periphery of photoreceptors. Collectively these elucidate the role of PRPF8/Brr2 regulatory mechanisms in splicing and the molecular basis of retinal disease, informing therapeutic approaches.

Spatial Gene Expression Changes in the Mouse Heart After Base-Targeted Irradiation.

PURPOSE: Radiation cardiotoxicity (RC) is a clinically significant adverse effect of treatment for patients with thoracic malignancies. Clinical studies in lung cancer have indicated that heart substructures are not uniformly radiosensitive, and that dose to the heart base drives RC. In this study, we aimed to characterize late changes in gene expression using spatial transcriptomics in a mouse model of base regional radiosensitivity. METHODS AND MATERIALS: An aged female C57BL/6 mouse was irradiated with 16 Gy delivered to the cranial third of the heart using a 6 × 9 mm parallel opposed beam geometry on a small animal radiation research platform, and a second mouse was sham-irradiated. After echocardiography, whole hearts were collected at 30 weeks for spatial transcriptomic analysis to map gene expression changes occurring in different regions of the partially irradiated heart. Cardiac regions were manually annotated on the capture slides and the gene expression profiles compared across different regions. RESULTS: Ejection fraction was reduced at 30 weeks after a 16 Gy irradiation to the heart base, compared with the sham-irradiated controls. There were markedly more significant gene expression changes within the irradiated regions compared with nonirradiated regions. Variation was observed in the transcriptomic effects of radiation on different cardiac base structures (eg, between the right atrium [n = 86 dysregulated genes], left atrium [n = 96 dysregulated genes], and the vasculature [n = 129 dysregulated genes]). Disrupted biological processes spanned extracellular matrix as well as circulatory, neuronal, and contractility activities. CONCLUSIONS: This is the first study to report spatially resolved gene expression changes in irradiated tissues. Examination of the regional radiation response in the heart can help to further our understanding of the cardiac base's radiosensitivity and support the development of actionable targets for pharmacologic intervention and biologically relevant dose constraints.

pRB-Depleted Pluripotent Stem Cell Retinal Organoids Recapitulate Cell State Transitions of Retinoblastoma Development and Suggest an Important Role for pRB in Retinal Cell Differentiation.

Retinoblastoma (Rb) is a childhood cancer of the developing retina, accounting for up to 17% of all tumors in infancy. To gain insights into the transcriptional events of cell state transitions during Rb development, we established 2 disease models via retinal organoid differentiation of a pRB (retinoblastoma protein)-depleted human embryonic stem cell line (RB1-null hESCs) and a pRB patient-specific induced pluripotent (iPSC) line harboring a RB1 biallelic mutation (c.2082delC). Both models were characterized by pRB depletion and accumulation of retinal progenitor cells at the expense of amacrine, horizontal and retinal ganglion cells, which suggests an important role for pRB in differentiation of these cell lineages. Importantly, a significant increase in the fraction of proliferating cone precursors (RXRγ+Ki67+) was observed in both pRB-depleted organoid models, which were defined as Rb-like clusters by single-cell RNA-Seq analysis. The pRB-depleted retinal organoids displayed similar features to Rb tumors, including mitochondrial cristae aberrations and rosette-like structures, and were able to undergo cell growth in an anchorage-independent manner, indicative of cell transformation in vitro. In both models, the Rb cones expressed retinal ganglion and horizontal cell markers, a novel finding, which could help to better characterize these tumors with possible therapeutic implications. Application of Melphalan, Topotecan, and TW-37 led to a significant reduction in the fraction of Rb proliferating cone precursors, validating the suitability of these in vitro models for testing novel therapeutics for Rb.

Human Retinal Organoids Provide a Suitable Tool for Toxicological Investigations: A Comprehensive Validation Using Drugs and Compounds Affecting the Retina.

Retinal drug toxicity screening is essential for the development of safe treatment strategies for a large number of diseases. To this end, retinal organoids derived from human pluripotent stem cells (hPSCs) provide a suitable screening platform due to their similarity to the human retina and the ease of generation in large-scale formats. In this study, two hPSC cell lines were differentiated to retinal organoids, which comprised all key retinal cell types in multiple nuclear and synaptic layers. Single-cell RNA-Seq of retinal organoids indicated the maintenance of retinal ganglion cells and development of bipolar cells: both cell types segregated into several subtypes. Ketorolac, digoxin, thioridazine, sildenafil, ethanol, and methanol were selected as key compounds to screen on retinal organoids because of their well-known retinal toxicity profile described in the literature. Exposure of the hPSC-derived retinal organoids to digoxin, thioridazine, and sildenafil resulted in photoreceptor cell death, while digoxin and thioridazine additionally affected all other cell types, including Müller glia cells. All drug treatments caused activation of astrocytes, indicated by dendrites sprouting into neuroepithelium. The ability to respond to light was preserved in organoids although the number of responsive retinal ganglion cells decreased after drug exposure. These data indicate similar drug effects in organoids to those reported in in vivo models and/or in humans, thus providing the first robust experimental evidence of their suitability for toxicological studies.

Conjunctival epithelial cells resist productive SARS-CoV-2 infection.

Conjunctival epithelial cells, which express viral-entry receptors angiotensin-converting enzyme 2 (ACE2) and transmembrane protease serine type 2 (TMPRSS2), constitute the largest exposed epithelium of the ocular surface tissue and may represent a relevant viral-entry route. To address this question, we generated an organotypic air-liquid-interface model of conjunctival epithelium, composed of basal, suprabasal, and superficial epithelial cells, and fibroblasts, which could be maintained successfully up to day 75 of differentiation. Using single-cell RNA sequencing (RNA-seq), with complementary imaging and virological assays, we observed that while all conjunctival cell types were permissive to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genome expression, a productive infection did not ensue. The early innate immune response to SARS-CoV-2 infection in conjunctival cells was characterised by a robust autocrine and paracrine NF-κB activity, without activation of antiviral interferon signalling. Collectively, these data enrich our understanding of SARS-CoV-2 infection at the human ocular surface, with potential implications for the design of preventive strategies and conjunctival transplantation.

Retinal pigment epithelium extracellular vesicles are potent inducers of age-related macular degeneration disease phenotype in the outer retina.

Age-related macular degeneration (AMD) is a leading cause of blindness. Vision loss is caused by the retinal pigment epithelium (RPE) and photoreceptors atrophy and/or retinal and choroidal angiogenesis. Here we use AMD patient-specific RPE cells with the Complement Factor H Y402H high-risk polymorphism to perform a comprehensive analysis of extracellular vesicles (EVs), their cargo and role in disease pathology. We show that AMD RPE is characterised by enhanced polarised EV secretion. Multi-omics analyses demonstrate that AMD RPE EVs carry RNA, proteins and lipids, which mediate key AMD features including oxidative stress, cytoskeletal dysfunction, angiogenesis and drusen accumulation. Moreover, AMD RPE EVs induce amyloid fibril formation, revealing their role in drusen formation. We demonstrate that exposure of control RPE to AMD RPE apical EVs leads to the acquisition of AMD features such as stress vacuoles, cytoskeletal destabilization and abnormalities in the morphology of the nucleus. Retinal organoid treatment with apical AMD RPE EVs leads to disrupted neuroepithelium and the appearance of cytoprotective alpha B crystallin immunopositive cells, with some co-expressing retinal progenitor cell markers Pax6/Vsx2, suggesting injury-induced regenerative pathways activation. These findings indicate that AMD RPE EVs are potent inducers of AMD phenotype in the neighbouring RPE and retinal cells.

Dissecting the Transcriptional and Chromatin Accessibility Heterogeneity of Proliferating Cone Precursors in Human Retinoblastoma Tumors by Single Cell Sequencing-Opening Pathways to New Therapeutic Strategies?

Purpose: Retinoblastoma (Rb) is a malignant neoplasm arising during retinal development from mutations in the RB1 gene. Loss or inactivation of both copies of RB1 results in initiation of retinoblastoma tumors; however, additional genetic changes are needed for the continued growth and spread of the tumor. Ex vivo research has shown that in humans, retinoblastoma may initiate from RB1-depleted cone precursors. Notwithstanding, it has not been possible to assess the full spectrum of clonal types within the tumor itself in vivo and the molecular changes occurring at the cells of origin, enabling their malignant conversion. To overcome these challenges, we have performed the first single cell (sc) RNA- and ATAC-Seq analyses of primary tumor tissues, enabling us to dissect the transcriptional and chromatin accessibility heterogeneity of proliferating cone precursors in human Rb tumors. Methods: Two Rb tumors each characterized by two pathogenic RB1 mutations were dissociated to single cells and subjected to scRNA-Seq and scATAC-Seq using the 10× Genomics platform. In addition, nine human embryonic and fetal retina samples were dissociated to single cells and subjected to scRNA- and ATAC-Seq analyses. The scRNA- and ATAC-Seq data were embedded using Uniform Manifold Approximation and Projection and clustered with Seurat graph-based clustering. Integrated scATAC-Seq analysis of Rb tumors and human embryonic/fetal retina samples was performed to identify Rb cone enriched subclusters. Pseudo time analysis of proliferating cones in the Rb samples was performed with Monocle. Ingenuity Pathway Analysis was used to identify the signaling pathway and upstream regulators in the Rb cone-enriched subclusters. Results: Our single cell analyses revealed the predominant presence of cone precursors at different stages of the cell cycle in the Rb tumors and among those identified the G2/M subset as the cell type of origin. scATAC-Seq analysis identified two Rb enriched cone subclusters, each characterized by activation of different upstream regulators and signaling pathways, enabling proliferating cone precursors to escape cell cycle arrest and/or apoptosis. Conclusions: Our study provides evidence of Rb tumor heterogeneity and defines molecular pathways that can be targeted to define new treatment strategies.