Multi-species single-cell transcriptomic analysis of ocular compartment regulons

scQA: A dual-perspective cell type identification model for single cell transcriptome data

Single-cell RNA sequencing technologies have been pivotal in advancing the development of algorithms for clustering heterogeneous cell populations. Existing methods for utilizing scRNA-seq data to identify cell types tend to neglect the beneficial impact of dropout events and perform clustering focusing solely on quantitative perspective. Here, we introduce a novel method named scQA, notable for its ability to concurrently identify cell types and cell type-specific key genes from both qualitative and quantitative perspectives. In contrast to other methods, scQA not only identifies cell types but also extracts key genes associated with these cell types, enabling bidirectional clustering for scRNA-seq data. Through an iterative process, our approach aims to minimize the number of landmarks to approximately a dozen while maximizing the inclusion of quasi-trend-preserved genes with dropouts both qualitatively and quantitatively. It then clusters cells by employing an ingenious label propagation strategy, obviating the requirement for a predetermined number of cell types. Validated on 20 publicly available scRNA-seq datasets, scQA consistently outperforms other salient tools. Furthermore, we confirm the effectiveness and potential biological significance of the identified key genes through both external and internal validation. In conclusion, scQA emerges as a valuable tool for investigating cell heterogeneity due to its distinctive fusion of qualitative and quantitative facets, along with bidirectional clustering capabilities. Furthermore, it can be seamlessly integrated into border scRNA-seq analyses. The source codes are publicly available at https://github.com/LD-Lyndee/scQA.

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.

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.

A Non-canonical Excitatory PV RGC-PV SC Visual Pathway for Mediating the Looming-evoked Innate Defensive Response

Parvalbumin-positive retinal ganglion cells (PV+ RGCs) are an essential subset of RGCs found in various species. However, their role in transmitting visual information remains unclear. Here, we characterized PV+ RGCs in the retina and explored the functions of the PV+ RGC-mediated visual pathway. By applying multiple viral tracing strategies, we investigated the downstream of PV+ RGCs across the whole brain. Interestingly, we found that the PV+ RGCs provided direct monosynaptic input to PV+ excitatory neurons in the superficial layers of the superior colliculus (SC). Ablation or suppression of SC-projecting PV+ RGCs abolished or severely impaired the flight response to looming visual stimuli in mice without affecting visual acuity. Furthermore, using transcriptome expression profiling of individual cells and immunofluorescence colocalization for RGCs, we found that PV+ RGCs are predominant glutamatergic neurons. Thus, our findings indicate the critical role of PV+ RGCs in an innate defensive response and suggest a non-canonical subcortical visual pathway from excitatory PV+ RGCs to PV+ SC neurons that regulates looming visual stimuli. These results provide a potential target for intervening and treating diseases related to this circuit, such as schizophrenia and autism.

iPSC-derived cells stimulate ABCG2+ /NES+ endogenous trabecular meshwork cell proliferation and tissue regeneration

A major risk factor for glaucoma, the first leading cause of irreversible blindness worldwide, is the decellularisation of the trabecular meshwork (TM) in the conventional outflow pathway. Stem cell-based therapy, particularly the utilisation of induced pluripotent stem cells (iPSCs), presents an enticing potential for tissue regeneration and intraocular pressure (IOP) maintenance in glaucoma. We have previously observed that differentiated iPSCs can stimulate endogenous cell proliferation in the TM, a pivotal factor in TM regeneration and aqueous humour outflow restoration. In this study, we investigated the response of TM cells in vivo after interacting with iPSC-derived cells and identified two subpopulations responsible for this relatively long-term tissue regeneration: ATP Binding Cassette Subfamily G Member 2 (ABCG2)-positive cells and Nestin (NES)-positive cells. We further uncovered that alterations of these responsive cells are linked to ageing and different glaucoma etiologies, suggesting that ABCG2+ subpopulation decellularization could serve as a potential risk factor for TM decellularization in glaucoma. Taken together, our findings illustrated the proliferative subpopulations in the conventional outflow pathway when stimulated with iPSC-derived cells and defined them as TM precursors, which may be applied to develop novel therapeutic approaches for glaucoma.

The single-cell transcriptomic atlas and RORA-mediated 3D epigenomic remodeling in driving corneal epithelial differentiation

Proper differentiation of corneal epithelial cells (CECs) from limbal stem/progenitor cells (LSCs) is required for maintenance of ocular homeostasis and clear vision. Here, using a single-cell transcriptomic atlas, we delineate the comprehensive and refined molecular regulatory dynamics during human CEC development and differentiation. We find that RORA is a CEC-specific molecular switch that initiates and drives LSCs to differentiate into mature CECs by activating PITX1. RORA dictates CEC differentiation by establishing CEC-specific enhancers and chromatin interactions between CEC gene promoters and distal regulatory elements. Conversely, RORA silences LSC-specific promoters and disrupts promoter-anchored chromatin loops to turn off LSC genes. Collectively, our work provides detailed and comprehensive insights into the transcriptional dynamics and RORA-mediated epigenetic remodeling underlying human corneal epithelial differentiation.

NRG1 Regulates Proliferation, Migration and Differentiation of Human Limbal Epithelial Stem Cells

Limbal epithelial stem/progenitor cells (LESCs) proliferate, migrate and differentiate into mature corneal epithelium cells (CECs) that cover the ocular surface. LESCs play a crucial role in the maintenance and regeneration of the corneal epithelium, and their dysfunction can lead to various corneal diseases. Neuregulin 1 (NRG1) is a member of the epidermal growth factor family that regulates the growth and differentiation of epithelial tissues. Here, we depicted the dynamic transcriptomic profiles during human CEC differentiation, identifying six gene co-expression modules that were specific to different differentiation stages. We found that the expression of NRG1 was high in human LESCs and decreased dramatically upon differentiation. Knockdown of NRG1 significantly inhibited LESC proliferation and upregulated the expression of the terminal differentiation marker genes KRT3, KRT12 and CLU. In addition, the scratch wound closure assay showed that knockdown of NRG1 attenuated wound closure of LESCs over 24 h. Together, we dissected the transcriptional regulatory dynamics during CEC differentiation and identified NRG1 as a key regulator that promoted LESC proliferation and migration and maintained the undifferentiated state.

Single-cell transcriptomics of the ocular anterior segment: a comprehensive review

Elucidating the cellular and genetic composition of ocular tissues is essential for uncovering the pathophysiology of ocular diseases. Since the introduction of single-cell RNA sequencing (scRNA-seq) in 2009, vision researchers have performed extensive single-cell analyses to better understand transcriptome complexity and heterogeneity of ocular structures. This technology has revolutionized our ability to identify rare cell populations and to make cross-species comparisons of gene expression in both steady state and disease conditions. Importantly, single-cell transcriptomic analyses have enabled the identification of cell-type specific gene markers and signalling pathways between ocular cell populations. While most scRNA-seq studies have been conducted on retinal tissues, large-scale transcriptomic atlases pertaining to the ocular anterior segment have also been constructed in the past three years. This timely review provides vision researchers with an overview of scRNA-seq experimental design, technical limitations, and clinical applications in a variety of anterior segment-related ocular pathologies. We review open-access anterior segment-related scRNA-seq datasets and illustrate how scRNA-seq can be an indispensable tool for the development of targeted therapeutics.

Liquid-biopsy proteomics combined with AI identifies cellular drivers of eye aging and disease in vivo

Single-cell analysis in living humans is essential for understanding disease mechanisms, but it is impractical in non-regenerative organs, such as the eye and brain, because tissue biopsies would cause serious damage. We resolve this problem by integrating proteomics of liquid biopsies with single-cell transcriptomics from all known ocular cell types to trace the cellular origin of 5,953 proteins detected in the aqueous humor. We identified hundreds of cell-specific protein markers, including for individual retinal cell types. Surprisingly, our results reveal that retinal degeneration occurs in Parkinson's disease, and the cells driving diabetic retinopathy switch with disease stage. Finally, we developed artificial intelligence (AI) models to assess individual cellular aging and found that many eye diseases not associated with chronological age undergo accelerated molecular aging of disease-specific cell types. Our approach, which can be applied to other organ systems, has the potential to transform molecular diagnostics and prognostics while uncovering new cellular disease and aging mechanisms.

CellAnn: a comprehensive, super-fast, and user-friendly single-cell annotation web server

MOTIVATION

Single-cell sequencing technology has become a routine in studying many biological problems. A core step of analyzing single-cell data is the assignment of cell clusters to specific cell types. Reference-based methods are proposed for predicting cell types for single-cell clusters. However, the scalability and lack of preprocessed reference datasets prevent them from being practical and easy to use.

RESULTS

Here, we introduce a reference-based cell annotation web server, CellAnn, which is super-fast and easy to use. CellAnn contains a comprehensive reference database with 204 human and 191 mouse single-cell datasets. These reference datasets cover 32 organs. Furthermore, we developed a cluster-to-cluster alignment method to transfer cell labels from the reference to the query datasets, which is superior to the existing methods with higher accuracy and higher scalability. Finally, CellAnn is an online tool that integrates all the procedures in cell annotation, including reference searching, transferring cell labels, visualizing results, and harmonizing cell annotation labels. Through the user-friendly interface, users can identify the best annotation by cross-validating with multiple reference datasets. We believe that CellAnn can greatly facilitate single-cell sequencing data analysis.

AVAILABILITY AND IMPLEMENTATION

The web server is available at www.cellann.io, and the source code is available at https://github.com/Pinlyu3/CellAnn_shinyapp.

PLAE Web App Enables Powerful Searching and Multiple Visualizations Across One Million Unified Single-Cell Ocular Transcriptomes

Purpose

To create a high-performance reactive web application to query single-cell gene expression data across cell type, species, study, and other factors.

Methods

We updated the content and structure of the underlying data (single cell Eye in a Disk [scEiaD]) and wrote the web application PLAE (https://plae.nei.nih.gov) to visualize and explore the data.

Results

The new portal provides quick visualization of over a million individual cells from vertebrate eye and body transcriptomes encompassing four species, 60 cell types, six ocular tissues, and 23 body tissues across 35 publications. To demonstrate the value of this unified pan-eye dataset, we replicated known neurogenic and cone macula markers in addition to proposing six new cone human region markers.

Conclusions

The PLAE web application offers the eye community a powerful and quick means to test hypotheses related to gene expression across a highly diverse, community-derived database.

Translational Relevance

The PLAE resource enables any researcher or clinician to study and research gene expression patterning across a wide variety of curated ocular cell types with a responsive web app.

Single-nucleus multiomic mapping of m6A methylomes and transcriptomes in native populations of cells with sn-m6A-CT

N6-methyladenosine (m6A) RNA modification plays important roles in the governance of gene expression and is temporally regulated in different cell states. In contrast to global m6A profiling in bulk sequencing, single-cell technologies for analyzing m6A heterogeneity are not extensively established. Here, we developed single-nucleus m6A-CUT&Tag (sn-m6A-CT) for simultaneous profiling of m6A methylomes and transcriptomes within a single nucleus using mouse embryonic stem cells (mESCs). m6A-CT is capable of enriching m6A-marked RNA molecules in situ, without isolating RNAs from cells. We adapted m6A-CT to the droplet-based single-cell omics platform and demonstrated high-throughput performance in analyzing nuclei isolated from thousands of cells from various cell types. We show that sn-m6A-CT profiling is sufficient to determine cell identity and allows the generation of cell-type-specific m6A methylome landscapes from heterogeneous populations. These indicate that sn-m6A-CT provides additional dimensions to multimodal datasets and insights into epitranscriptomic landscape in defining cell fate identity and states.

Transcriptomic analysis of the ocular posterior segment completes a cell atlas of the human eye

Although the visual system extends through the brain, most vision loss originates from defects in the eye. Its central element is the neural retina, which senses light, processes visual signals, and transmits them to the rest of the brain through the optic nerve (ON). Surrounding the retina are numerous other structures, conventionally divided into anterior and posterior segments. Here, we used high-throughput single-nucleus RNA sequencing (snRNA-seq) to classify and characterize cells in six extraretinal components of the posterior segment: ON, optic nerve head (ONH), peripheral sclera, peripapillary sclera (PPS), choroid, and retinal pigment epithelium (RPE). Defects in each of these tissues are associated with blinding diseases-for example, glaucoma (ONH and PPS), optic neuritis (ON), retinitis pigmentosa (RPE), and age-related macular degeneration (RPE and choroid). From ~151,000 single nuclei, we identified 37 transcriptomically distinct cell types, including multiple types of astrocytes, oligodendrocytes, fibroblasts, and vascular endothelial cells. Our analyses revealed a differential distribution of many cell types among distinct structures. Together with our previous analyses of the anterior segment and retina, the data presented here complete a "Version 1" cell atlas of the human eye. We used this atlas to map the expression of >180 genes associated with the risk of developing glaucoma, which is known to involve ocular tissues in both anterior and posterior segments as well as the neural retina. Similar methods can be used to investigate numerous additional ocular diseases, many of which are currently untreatable.

Low-dose radiation induces unstable gene expression in developing human iPSC-derived retinal ganglion organoids

The effects of low-dose radiation on undifferentiated cells carry important implications. However, the effects on developing retinal cells remain unclear. Here, we analyzed the gene expression characteristics of neuronal organoids containing immature human retinal cells under low-dose radiation and predicted their changes. Developing retinal cells generated from human induced pluripotent stem cells (iPSCs) were irradiated with either 30 or 180 mGy on days 4-5 of development for 24 h. Genome-wide gene expression was observed until day 35. A knowledge-based pathway analysis algorithm revealed fluctuations in Rho signaling and many other pathways. After a month, the levels of an essential transcription factor of eye development, the proportion of paired box 6 (PAX6)-positive cells, and the proportion of retinal ganglion cell (RGC)-specific transcription factor POU class 4 homeobox 2 (POU4F2)-positive cells increased with 30 mGy of irradiation. In contrast, they decreased after 180 mGy of irradiation. Activation of the "development of neurons" pathway after 180 mGy indicated the dedifferentiation and development of other neural cells. Fluctuating effects after low-dose radiation exposure suggest that developing retinal cells employ hormesis and dedifferentiation mechanisms in response to stress.

Postmortem Ultrastructural Analysis of the Retina from COVID-19 Deceased Patients

PURPOSE

COVID-19 (coronavirus disease 2019) is an infectious disease caused by SARS-CoV-2, first reported in 2019 in Wuhan, China. Among the common complications is a pro-inflammatory and hypercoagulative response that compromises the vasculature among various organs.

METHODS

In this report, we present the postmortem retinal findings of five patients observed by means of optical microscopy and transmission and scanning electron microscopy techniques.

RESULTS

Clinical manifestations such as retinal hemorrhages and exacerbated inflammatory infiltrate, altered ultra structure with swollen mitochondria and pyknotic cells in both layers of the retina were observed in all analyzed eyes.

CONCLUSION

Our data point to the fragility of this tissue in cases of severe COVID-19.

Single-Cell RNA Sequencing: Opportunities and Challenges for Studies on Corneal Biology in Health and Disease

The structure and major cell types of the multi-layer human cornea have been extensively studied. However, various cell states in specific cell types and key genes that define the cell states are not fully understood, hindering our comprehension of corneal homeostasis, related diseases, and therapeutic discovery. Single-cell RNA sequencing is a revolutionary and powerful tool for identifying cell states within tissues such as the cornea. This review provides an overview of current single-cell RNA sequencing studies on the human cornea, highlighting similarities and differences between them, and summarizing the key genes that define corneal cell states reported in these studies. In addition, this review discusses the opportunities and challenges of using single-cell RNA sequencing to study corneal biology in health and disease.

Single cell transcriptome analyses of the developing zebrafish eye- perspectives and applications

Within a relatively short period of time, single cell transcriptome analyses (SCT) have become increasingly ubiquitous with transcriptomic research, uncovering plentiful details that boost our molecular understanding of various biological processes. Stemming from SCT analyses, the ever-growing number of newly assigned genetic markers increases our understanding of general function and development, while providing opportunities for identifying genes associated with disease. SCT analyses have been carried out using tissue from numerous organisms. However, despite the great potential of zebrafish as a model organism, other models are still preferably used. In this mini review, we focus on eye research as an example of the advantages in using zebrafish, particularly its usefulness for single cell transcriptome analyses of developmental processes. As studies have already shown, the unique opportunities offered by zebrafish, including similarities to the human eye, in combination with the possibility to analyze and extract specific cells at distinct developmental time points makes the model a uniquely powerful one. Particularly the practicality of collecting large numbers of embryos and therefore isolation of sufficient numbers of developing cells is a distinct advantage compared to other model organisms. Lastly, the advent of highly efficient genetic knockouts methods offers opportunities to characterize target gene function in a more cost-efficient way. In conclusion, we argue that the use of zebrafish for SCT approaches has great potential to further deepen our molecular understanding of not only eye development, but also many other organ systems.

A multi-omics atlas of the human retina at single-cell resolution

Cell classes in the human retina are highly heterogeneous with their abundance varying by several orders of magnitude. Here, we generated and integrated a multi-omics single-cell atlas of the adult human retina, including more than 250,000 nuclei for single-nuclei RNA-seq and 137,000 nuclei for single-nuclei ATAC-seq. Cross-species comparison of the retina atlas among human, monkey, mice, and chicken revealed relatively conserved and non-conserved types. Interestingly, the overall cell heterogeneity in primate retina decreases compared with that of rodent and chicken retina. Through integrative analysis, we identified 35,000 distal cis-element-gene pairs, constructed transcription factor (TF)-target regulons for more than 200 TFs, and partitioned the TFs into distinct co-active modules. We also revealed the heterogeneity of the cis-element-gene relationships in different cell types, even from the same class. Taken together, we present a comprehensive single-cell multi-omics atlas of the human retina as a resource that enables systematic molecular characterization at individual cell-type resolution.

Single-cell RNA sequencing in cornea research: Insights into limbal stem cells and their niche regulation

The corneal epithelium is composed of stratified squamous epithelial cells on the outer surface of the eye, which acts as a protective barrier and is critical for clear and stable vision. Its continuous renewal or wound healing depends on the proliferation and differentiation of limbal stem cells (LSCs), a cell population that resides at the limbus in a highly regulated niche. Dysfunction of LSCs or their niche can cause limbal stem cell deficiency, a disease that is manifested by failed epithelial wound healing or even blindness. Nevertheless, compared to stem cells in other tissues, little is known about the LSCs and their niche. With the advent of single-cell RNA sequencing, our understanding of LSC characteristics and their microenvironment has grown considerably. In this review, we summarized the current findings from single-cell studies in the field of cornea research and focused on important advancements driven by this technology, including the heterogeneity of the LSC population, novel LSC markers and regulation of the LSC niche, which will provide a reference for clinical issues such as corneal epithelial wound healing, ocular surface reconstruction and interventions for related diseases.

Transcriptomic Analysis of the Ocular Posterior Segment Completes a Cell Atlas of the Human Eye

Although the visual system extends through the brain, most vision loss originates from defects in the eye. Its central element is the neural retina, which senses light, processes visual signals, and transmits them to the rest of the brain through the optic nerve (ON). Surrounding the retina are numerous other structures, conventionally divided into anterior and posterior segments. Here we used high-throughput single nucleus RNA sequencing (snRNA-seq) to classify and characterize cells in the extraretinal components of the posterior segment: ON, optic nerve head (ONH), peripheral sclera, peripapillary sclera (PPS), choroid, and retinal pigment epithelium (RPE). Defects in each of these tissues are associated with blinding diseases - for example, glaucoma (ONH and PPS), optic neuritis (ON), retinitis pigmentosa (RPE), and age-related macular degeneration (RPE and choroid). From ∼151,000 single nuclei, we identified 37 transcriptomically distinct cell types, including multiple types of astrocytes, oligodendrocytes, fibroblasts, and vascular endothelial cells. Our analyses revealed a differential distribution of many cell types among distinct structures. Together with our previous analyses of the anterior segment and retina, the new data complete a "Version 1" cell atlas of the human eye. We used this atlas to map the expression of >180 genes associated with the risk of developing glaucoma, which is known to involve ocular tissues in both anterior and posterior segments as well as neural retina. Similar methods can be used to investigate numerous additional ocular diseases, many of which are currently untreatable.

A temporal single cell transcriptome atlas of zebrafish anterior segment development

Anterior segment dysgenesis (ASD), resulting in vision impairment, stems from maldevelopment of anterior segment (AS) tissues. Incidence of ASD has been linked to malfunction of periocular mesenchyme cells (POM). POM cells specify into anterior segment mesenchyme (ASM) cells which colonize and produce AS tissues. In this study we uncover ASM developmental trajectories associated with formation of the AS. Using a transgenic line of zebrafish that fluorescently labels the ASM throughout development, Tg[foxc1b:GFP], we isolated GFP+ ASM cells at several developmental timepoints (48-144 hpf) and performed single cell RNA sequencing. Clustering analysis indicates subdifferentiation of ASM as early as 48 hpf and subsequent diversification into corneal epithelium/endothelium/stroma, or annular ligament (AL) lineages. Tracking individual clusters reveals common developmental pathways, up to 72 hpf, for the AL and corneal endothelium/stroma and distinct pathways for corneal epithelium starting at 48 hpf. Spatiotemporal validation of over 80 genes found associated with AS development demonstrates a high degree of conservation with mammalian trabecular meshwork and corneal tissues. In addition, we characterize thirteen novel genes associated with annular ligament and seven with corneal development. Overall, the data provide a molecular verification of the long-standing hypothesis that POM derived ASM give rise to AS tissues and highlight the high degree of conservation between zebrafish and mammals.

Quantitative systems pharmacology of the eye: Tools and data for ocular QSP

Good eyesight belongs to the most-valued attributes of health, and diseases of the eye are a significant healthcare burden. Case numbers are expected to further increase in the next decades due to an aging society. The development of drugs in ophthalmology, however, is difficult due to limited accessibility of the eye, in terms of drug administration and in terms of sampling of tissues for drug pharmacokinetics (PKs) and pharmacodynamics (PDs). Ocular quantitative systems pharmacology models provide the opportunity to describe the distribution of drugs in the eye as well as the resulting drug-response in specific segments of the eye. In particular, ocular physiologically-based PK (PBPK) models are necessary to describe drug concentration levels in different regions of the eye. Further, ocular effect models using molecular data from specific cellular systems are needed to develop dose-response correlations. We here describe the current status of PK/PBPK as well as PD models for the eyes and discuss cellular systems, data repositories, as well as animal models in ophthalmology. The application of the various concepts is highlighted for the development of new treatments for postoperative fibrosis after glaucoma surgery.

Whole Exome Sequencing Reveals Novel Candidate Genes in Familial Forms of Glaucomatous Neurodegeneration

Glaucoma is the largest cause of irreversible blindness with a multifactorial genetic etiology. This study explores novel genes and gene networks in familial forms of primary open angle glaucoma (POAG) and primary angle closure glaucoma (PACG) to identify rare mutations with high penetrance. Thirty-one samples from nine MYOC-negative families (five POAG and four PACG) underwent whole-exome sequencing and analysis. A set of prioritized genes and variations were screened in an independent validation cohort of 1536 samples and the whole-exome data from 20 sporadic patients. The expression profiles of the candidate genes were analyzed in 17 publicly available expression datasets from ocular tissues and single cells. Rare, deleterious SNVs in AQP5, SRFBP1, CDH6 and FOXM1 from POAG families and in ACACB, RGL3 and LAMA2 from PACG families were found exclusively in glaucoma cases. AQP5, SRFBP1 and CDH6 also revealed significant altered expression in glaucoma in expression datasets. Single-cell expression analysis revealed enrichment of identified candidate genes in retinal ganglion cells and corneal epithelial cells in POAG; whereas for PACG families, retinal ganglion cells and Schwalbe's Line showed enriched expression. Through an unbiased exome-wide search followed by validation, we identified novel candidate genes for familial cases of POAG and PACG. The SRFBP1 gene found in a POAG family is located within the GLC1M locus on Chr5q. Pathway analysis of candidate genes revealed enrichment of extracellular matrix organization in both POAG and PACG.

Towards a clinically-based common coordinate framework for the human gut cell atlas: the gut models

BACKGROUND

The Human Cell Atlas resource will deliver single cell transcriptome data spatially organised in terms of gross anatomy, tissue location and with images of cellular histology. This will enable the application of bioinformatics analysis, machine learning and data mining revealing an atlas of cell types, sub-types, varying states and ultimately cellular changes related to disease conditions. To further develop the understanding of specific pathological and histopathological phenotypes with their spatial relationships and dependencies, a more sophisticated spatial descriptive framework is required to enable integration and analysis in spatial terms.

METHODS

We describe a conceptual coordinate model for the Gut Cell Atlas (small and large intestines). Here, we focus on a Gut Linear Model (1-dimensional representation based on the centreline of the gut) that represents the location semantics as typically used by clinicians and pathologists when describing location in the gut. This knowledge representation is based on a set of standardised gut anatomy ontology terms describing regions in situ, such as ileum or transverse colon, and landmarks, such as ileo-caecal valve or hepatic flexure, together with relative or absolute distance measures. We show how locations in the 1D model can be mapped to and from points and regions in both a 2D model and 3D models, such as a patient's CT scan where the gut has been segmented.

RESULTS

The outputs of this work include 1D, 2D and 3D models of the human gut, delivered through publicly accessible Json and image files. We also illustrate the mappings between models using a demonstrator tool that allows the user to explore the anatomical space of the gut. All data and software is fully open-source and available online.

CONCLUSIONS

Small and large intestines have a natural "gut coordinate" system best represented as a 1D centreline through the gut tube, reflecting functional differences. Such a 1D centreline model with landmarks, visualised using viewer software allows interoperable translation to both a 2D anatomogram model and multiple 3D models of the intestines. This permits users to accurately locate samples for data comparison.

Ribosomal proteins regulate 2-cell-stage transcriptome in mouse embryonic stem cells

A rare sub-population of mouse embryonic stem cells (mESCs), the 2-cell-like cell, is defined by the expression of MERVL and 2-cell-stage-specific transcript (2C transcript). Here, we report that the ribosomal proteins (RPs) RPL14, RPL18, and RPL23 maintain the identity of mESCs and regulate the expression of 2C transcripts. Disregulation of the RPs induces DUX-dependent expression of 2C transcripts and alters the chromatin landscape. Mechanically, knockdown (KD) of RPs triggers the binding of RPL11 to MDM2, an interaction known to prevent P53 protein degradation. Increased P53 protein upon RP KD further activates its downstream pathways, including DUX. Our study delineates the critical roles of RPs in 2C transcript activation, ascribing a novel function to these essential proteins.

NAD salvage pathway machinery expression in normal and glaucomatous retina and optic nerve

Glaucoma is the leading cause of irreversible blindness and is a major health and economic burden. Current treatments do not address the neurodegenerative component of glaucoma. In animal models of glaucoma, the capacity to maintain retinal nicotinamide adenine dinucleotide (NAD) pools declines early during disease pathogenesis. Treatment with nicotinamide, an NAD precursor through the NAD salvage pathway, robustly protects against neurodegeneration in a number of glaucoma models and improves vision in existing glaucoma patients. However, it remains unknown in humans what retinal cell types are able to process nicotinamide to NAD and how these are affected in glaucoma. To address this, we utilized publicly available RNA-sequencing data (bulk, single cell, and single nucleus) and antibody labelling in highly preserved enucleated human eyes to identify expression of NAD synthesizing enzyme machinery. This identifies that the neural retina favors expression of the NAD salvage pathway, and that retinal ganglion cells are particularly enriched for these enzymes. NMNAT2, a key terminal enzyme in the salvage pathway, is predominantly expressed in retinal ganglion cell relevant layers of the retina and declines in glaucoma. These findings suggest that human retinal ganglion cells can directly utilize nicotinamide and could maintain a capacity to do so in glaucoma, showing promise for ongoing clinical trials.

ETS1-HMGA2 Axis Promotes Human Limbal Epithelial Stem Cell Proliferation

Purpose

This study aimed to investigate the role and molecular mechanism of ETS1 in the proliferation and differentiation of human limbal epithelial stem cells (LESCs).

Methods

RNA-seq and quantitative real-time PCR were used to determine gene expression changes when ETS1 and HMGA2 was knocked down using short-hairpin RNAs or overexpressed by lentivirus. Immunofluorescence and flow cytometry experiments were performed to assess the roles of ETS1 and HMGA2 in LESC proliferation. ETS1-bound cis-regulatory elements and target genes in LESCs were identified using chromatin immunoprecipitation sequencing. The epigenetic features of ETS1-binding sites were assessed by the published histone modification and chromatin accessibility profiles.

Results

ETS1 was robustly expressed in LESCs but dramatically reduced on differentiation into corneal epithelial cells (CECs). ETS1 knockdown in LESCs inhibited cellular proliferation and activated CEC markers (KRT3, KRT12, CLU, and ALDH3A1). When ETS1 was overexpressed during CEC differentiation, LESC-associated genes were upregulated while CEC-associated genes were downregulated. The genome-wide binding profile of ETS1 was identified in LESCs. ETS1 occupied H3K4me3-marked promoters and H3K27ac/H3K4me1-marked enhancers. ETS1-binding sites were also enriched for chromatin accessibility signal. HMGA2 showed a consistent expression pattern with ETS1. ETS1 activates HMAG2 by binding to its promoter. Knockdown and overexpression experiments suggested that HMGA2 can promote LESC proliferation and inhibits its differentiation.

Conclusions

ETS1 promotes LESC proliferation and inhibits its differentiation via activating HMGA2.

Single-cell omics: A new direction for functional genetic research in human diseases and animal models

Over the past decade, with the development of high-throughput single-cell sequencing technology, single-cell omics has been emerged as a powerful tool to understand the molecular basis of cellular mechanisms and refine our knowledge of diverse cell states. They can reveal the heterogeneity at different genetic layers and elucidate their associations by multiple omics analysis, providing a more comprehensive genetic map of biological regulatory networks. In the post-GWAS era, the molecular biological mechanisms influencing human diseases will be further elucidated by single-cell omics. This review mainly summarizes the development and trend of single-cell omics. This involves single-cell omics technologies, single-cell multi-omics technologies, multiple omics data integration methods, applications in various human organs and diseases, classic laboratory cell lines, and animal disease models. The review will reveal some perspectives for elucidating human diseases and constructing animal models.

COVID-19 signalome: Pathways for SARS-CoV-2 infection and impact on COVID-19 associated comorbidity

The COVID-19 pandemic has been the focus of research the past two years. The major breakthrough was made by discovering pathways related to SARS-CoV-2 infection through cellular interaction by angiotensin-converting enzyme (ACE2) and cytokine storm. The presence of ACE2 in lungs, intestines, cardiovascular tissues, brain, kidneys, liver, and eyes shows that SARS-CoV-2 may have targeted these organs to further activate intracellular signalling pathways that lead to cytokine release syndrome. It has also been reported that SARS-CoV-2 can hijack coatomer protein-I (COPI) for S protein retrograde trafficking to the endoplasmic reticulum-Golgi intermediate compartment (ERGIC), which, in turn, acts as the assembly site for viral progeny. In infected cells, the newly synthesized S protein in endoplasmic reticulum (ER) is transported first to the Golgi body, and then from the Golgi body to the ERGIC compartment resulting in the formation of specific a motif at the C-terminal end. This review summarizes major events of SARS-CoV-2 infection route, immune response following host-cell infection as an important factor for disease outcome, as well as comorbidity issues of various tissues and organs arising due to COVID-19. Investigations on alterations of host-cell machinery and viral interactions with multiple intracellular signaling pathways could represent a major factor in more effective disease management.

Transcriptomic profiling of iris tissue highlights LCK signaling and T cell-mediated immunity in Behcet's uveitis

Uveitis is the most common form of ocular lesions in Behcet's disease, severely affecting visual function. Molecular pathological changes of ocular lesions in patients with Behcet's uveitis (BU) are largely unknown. In this study, we performed the first comprehensive transcriptomic profiling of iris specimens from BU patients and healthy donors to provide an insight into intraocular immunopathogenesis. The mRNA sequencing identified 1633 differentially expressed genes (DEGs) between the BU group and healthy controls. GO functional enrichment analysis on DEGs showed that T cell activation was the most significantly enriched biological process. KEGG analysis of DEGs also revealed several prominently enriched T cell-related pathways, including the T cell receptor signaling pathway, Th17 cell differentiation, and Th1 and Th2 cell differentiation. The lymphocyte-specific protein tyrosine kinase (LCK) was identified as the key hub gene in the protein interaction network of DEGs. Western blot analysis further showed increased expression of active LCK in the BU group, suggesting activation of LCK signaling. Using publicly accessible single-cell RNA-sequencing data of the healthy iris, LCK was found to be expressed in clusters of activated T cells but not in other iris cell clusters, suggesting an overt association between LCK upregulation and T cell-mediated immune dysregulation. Additionally, 16 drugs were predicted to be potential inhibitors of LCK. Overall, these findings not only highlighted the central role of T cell-mediated immunity and previously unreported LCK signaling in intraocular immunopathogenesis but also revealed the potential value of LCK as a new therapeutic target for BU patients.

Cellular and molecular profiles of larval and adult Xenopus corneal epithelia resolved at the single-cell level

Corneal Epithelial Stem Cells (CESCs) and their proliferative progeny, the Transit Amplifying Cells (TACs), are responsible for homeostasis and maintaining corneal transparency. Owing to our limited knowledge of cell fates and gene activity within the cornea, the search for unique markers to identify and isolate these cells remains crucial for ocular surface reconstruction. We performed single-cell RNA sequencing of corneal cells from larval and adult stages of Xenopus. Our results indicate that as the cornea develops and matures, there is an increase in cellular diversity, which is accompanied by a substantial shift in transcriptional profile, gene regulatory network and cell-cell communication dynamics. Our data also reveals several novel genes expressed in corneal cells and changes in gene expression during corneal differentiation at both developmental time-points. Importantly, we identify specific basal cell clusters in both the larval and adult cornea that comprise a relatively undifferentiated cell type and express distinct stem cell markers, which we propose are the putative larval and adult CESCs, respectively. This study offers a detailed atlas of single-cell transcriptomes in the frog cornea. In the future, this work will be useful to elucidate the function of novel genes in corneal epithelial homeostasis, wound healing and regeneration.

Evaluation of choroidal thickness with OCT in COVID-19 patients with high D-dimer levels

To evaluate retinal and choroidal thickness with optical coherence tomography (OCT) to detect retinal and choroidal pathologies in coronavirus disease 2019 (COVID-19) patients with high D-dimer levels. Thirty patients who were hospitalized in the intensive care unit due to severe acute respiratory syndrome coronavirus 2 and whose D-dimer levels were high during this period, who applied to the internal medicine outpatient clinic between 15 and 30 days after discharge, and 30 healthy volunteers with similar age and gender as the control group was included in the study. After full ophthalmological examination, central foveal and choroidal thicknesses were evaluated using optical coherence tomography. Statistical analysis of the study data demonstrated that there was no significant difference between the groups in terms of age or gender (p > 0.05). There was also no statistically significant difference between the groups in terms of central foveal thickness, central choroidal thickness, or nasal 500, nasal 1500, temporal 500, or temporal 500 micron distances (p > 0.05 for all parameters). Choroidal and retinal vascular thicknesses were not affected in the short term in COVID-19 patients hospitalized in the intensive care unit.

The Essential Role of the Choriocapillaris in Vision: Novel Insights from Imaging and Molecular Biology

The choriocapillaris, a dense capillary network located at the posterior pole of the eye, is essential for supporting normal vision, supplying nutrients, and removing waste products from photoreceptor cells and the retinal pigment epithelium. The anatomical location, heterogeneity, and homeostatic interactions with surrounding cell types make the choroid complex to study both in vivo and in vitro. Recent advances in single-cell RNA sequencing, in vivo imaging, and in vitro cell modeling are vastly improving our knowledge of the choroid and its role in normal health and in age-related macular degeneration (AMD). Histologically, loss of endothelial cells (ECs) of the choriocapillaris occurs early in AMD concomitant with elevated formation of the membrane attack complex of complement. Advanced imaging has allowed us to visualize early choroidal blood flow changes in AMD in living patients, supporting histological findings of loss of choroidal ECs. Single-cell RNA sequencing is being used to characterize choroidal cell types transcriptionally and discover their altered patterns of gene expression in aging and disease. Advances in induced pluripotent stem cell protocols and 3D cultures will allow us to closely mimic the in vivo microenvironment of the choroid in vitro to better understand the mechanism leading to choriocapillaris loss in AMD.

Screening of primary open-angle glaucoma diagnostic markers based on immune-related genes and immune infiltration

Purpose

Primary open-angle glaucoma (POAG) continues to be a poorly understood disease. Although there were multiple researches on the identification of POAG biomarkers, few studies systematically revealed the immune-related cells and immune infiltration of POAG. Bioinformatics analyses of optic nerve (ON) and trabecular meshwork (TM) gene expression data were performed to further elucidate the immune-related genes of POAG and identify candidate target genes for treatment.

Methods

We performed a gene analysis of publicly available microarray data, namely, the GSE27276-GPL2507, GSE2378-GPL8300, GSE9944-GPL8300, and GSE9944-GPL571 datasets from the Gene Expression Omnibus database. The obtained datasets were used as input for parallel pathway analyses. Based on random forest and support vector machine (SVM) analysis to screen the key genes, significantly changed pathways were clustered into functional categories, and the results were further investigated. CIBERSORT was used to evaluate the infiltration of immune cells in POAG tissues. A network visualizing the differences between the data in the POAG and normal groups was created. GO and KEGG enrichment analyses were performed using the Metascape database. We divided the differentially expressed mRNAs into upregulated and downregulated groups and predicted the drug targets of the differentially expressed genes through the Connectivity Map (CMap) database.

Results

A total of 49 differentially expressed genes, including 19 downregulated genes and 30 upregulated genes, were detected. Five genes ((Keratin 14) KRT14, (Hemoglobin subunit beta) HBB, (Acyl-CoA Oxidase 2) ACOX2, (Hephaestin) HEPH and Keratin 13 (KRT13)) were significantly changed. The results showed that the expression profiles of drug disturbances, including those for avrainvillamide-analysis-3, cytochalasin-D, NPI-2358, oxymethylone and vinorelbine, were negatively correlated with the expression profiles of disease disturbances. This finding indicated that these drugs may reduce or even reverse the POAG disease state.

Conclusion

This study provides an overview of the processes involved in the molecular pathogenesis of POAG in the ON and TM. The findings provide a new understanding of the molecular mechanism of POAG from the perspective of immunology.

FOXD1 Is a Transcription Factor Important for Uveal Melanocyte Development and Associated with High-Risk Uveal Melanoma

Simple Summary

Despite successful treatment of primary uveal melanoma (UM), metastases still occur in approximately 50% of the patients. Unfortunately, little is known about the mechanism behind metastasized UM. By reanalyzing publicly available single-cell RNA sequencing data of embryonic zebrafish larvae and validating the results with UM data, we have identified five transcription regulators of interest: ELL2, KDM5B, REXO4, RBFOX2 and FOXD1. The most significant finding is FOXD1, which is nearly exclusively expressed in high-risk UM and is associated with poor survival. FOXD1 is a novel gene which could be involved in the metastatic capability of UM. Elucidating its function and role in metastatic UM could help to understand and develop treatment for UM.

Abstract

Uveal melanoma (UM) is a deadly ocular malignancy, originating from uveal melanocytes. Although much is known regarding prognostication in UM, the exact mechanism of metastasis is mostly unknown. Metastatic tumor cells are known to express a more stem-like RNA profile which is seen often in cell-specific embryonic development to induce tumor progression. Here, we identified novel transcription regulators by reanalyzing publicly available single cell RNA sequencing experiments. We identified five transcription regulators of interest: ELL2, KDM5B, REXO4, RBFOX2 and FOXD1. Our most significant finding is FOXD1, as this gene is nearly exclusively expressed in high-risk UM and its expression is associated with a poor prognosis. Even within the BAP1-mutated UM, the expression of FOXD1 is correlated with poor survival. FOXD1 is a novel factor which could potentially be involved in the metastatic capacity of high-risk UM. Elucidating the function of FOXD1 in UM could provide insight into the malignant transformation of uveal melanocytes, especially in high-risk UM.

Cell atlas of the human ocular anterior segment: Tissue-specific and shared cell types

The anterior segment of the eye consists of the cornea, iris, ciliary body, crystalline lens, and aqueous humor outflow pathways. Together, these tissues are essential for the proper functioning of the eye. Disorders of vision have been ascribed to defects in all of them; some disorders, including glaucoma and cataract, are among the most prevalent causes of blindness in the world. To characterize the cell types that compose these tissues, we generated an anterior segment cell atlas of the human eye using high-throughput single-nucleus RNA sequencing (snRNAseq). We profiled 195,248 nuclei from nondiseased anterior segment tissues of six human donors, identifying >60 cell types. Many of these cell types were discrete, whereas others, especially in the lens and cornea, formed continua corresponding to known developmental transitions that persist in adulthood. Having profiled each tissue separately, we performed an integrated analysis of the entire anterior segment, revealing that some cell types are unique to a single structure, whereas others are shared across tissues. The integrated cell atlas was then used to investigate cell type-specific expression patterns of more than 900 human ocular disease genes identified through either Mendelian inheritance patterns or genome-wide association studies.

Uncovering Pharmacological Opportunities for Cancer Stem Cells-A Systems Biology View

Cancer stem cells (CSCs) represent a small fraction of the total cancer cell population, yet they are thought to drive disease propagation, therapy resistance and relapse. Like healthy stem cells, CSCs possess the ability to self-renew and differentiate. These stemness phenotypes of CSCs rely on multiple molecular cues, including signaling pathways (for example, WNT, Notch and Hedgehog), cell surface molecules that interact with cellular niche components, and microenvironmental interactions with immune cells. Despite the importance of understanding CSC biology, our knowledge of how neighboring immune and tumor cell populations collectively shape CSC stemness is incomplete. Here, we provide a systems biology perspective on the crucial roles of cellular population identification and dissection of cell regulatory states. By reviewing state-of-the-art single-cell technologies, we show how innovative systems-based analysis enables a deeper understanding of the stemness of the tumor niche and the influence of intratumoral cancer cell and immune cell compositions. We also summarize strategies for refining CSC systems biology, and the potential role of this approach in the development of improved anticancer treatments. Because CSCs are amenable to cellular transitions, we envision how systems pharmacology can become a major engine for discovery of novel targets and drug candidates that can modulate state transitions for tumor cell reprogramming. Our aim is to provide deeper insights into cancer stemness from a systems perspective. We believe this approach has great potential to guide the development of more effective personalized cancer therapies that can prevent CSC-mediated relapse.