Corneal nonmyelinating Schwann cells illuminated by single-cell transcriptomics and visualized by protein biomarkers

Single cell analysis of short-term dry eye induced changes in cornea immune cell populations

Background

Dry eye causes corneal inflammation, epitheliopathy and sensorineural changes. This study evaluates the hypothesis that dry eye alters the percentages and transcriptional profiles of immune cell populations in the cornea.

Methods

Desiccating stress (DS) induced dry eye was created by pharmacologic suppression of tear secretion and exposure to drafty low humidity environment. Expression profiling of corneal immune cells was performed by single-cell RNA sequencing (scRNA-seq). Cell differentiation trajectories and cell fate were modeled through RNA velocity analysis. Confocal microscopy was used to immunodetect corneal immune cells. Irritation response to topical neurostimulants was assessed.

Results

Twelve corneal immune cell populations based on their transcriptional profiles were identified at baseline and consist of monocytes, resident (rMP) and MMP12/13 high macrophages, dendritic cells (cDC2), neutrophils, mast cells, pre T/B cells, and innate (γDT, ILC2, NK) and conventional T and B lymphocytes. T cells and resident macrophages (rMP) were the largest populations in the normal cornea comprising 18.6 and 18.2 percent, respectively. rMP increased to 55.2% of cells after 5 days of DS. Significant changes in expression of 1,365 genes (adj p < 0.0001) were noted in rMP with increases in cytokines and chemokines (Tnf, Cxcl1, Ccl12, Il1rn), inflammatory markers (Vcam, Adam17, Junb), the TAM receptor (Mertk), and decreases in complement and MHCII genes. A differentiation trajectory from monocytes to terminal state rMP was found. Phagocytosis, C-type lectin receptor signaling, NF-kappa B signaling and Toll-like receptor signaling were among the pathways with enhanced activity in these cells. The percentage of MRC1+ rMPs increased in the cornea and they were observed in the basal epithelium adjacent to epithelial nerve plexus. Concentration of the chemokine CXCL1 increased in the cornea and it heightened irritation/pain responses to topically applied hypertonic saline.

Conclusion

These findings indicate that DS recruits monocytes that differentiate to macrophages with increased expression of inflammation associated genes. The proximity of these macrophages to cornea nerves and their expression of neurosensitizers suggests they contribute to the corneal sensorineural changes in dry eye.

Modulating Growth Factor Receptor Signaling to Promote Corneal Epithelial Homeostasis

The corneal epithelium is the first anatomical barrier between the environment and the cornea; it is critical for proper light refraction onto the retina and prevents pathogens (e.g., bacteria, viruses) from entering the immune-privileged eye. Trauma to the highly innervated corneal epithelium is extremely painful and if not resolved quickly or properly, can lead to infection and ultimately blindness. The healthy eye produces its own growth factors and is continuously bathed in tear fluid that contains these proteins and other nutrients to maintain the rapid turnover and homeostasis of the ocular surface. In this article, we review the roles of growth factors in corneal epithelial homeostasis and regeneration and some of the limitations to their use therapeutically.

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.

Leprosy-specific subsets of macrophages and Schwann cells identified by single-cell RNA-sequencing

In Mycobacterium leprae (M. leprae)-infection, inflammatory cells' subsets and dynamics as well as the interactions with Schwann cells have remained elusive. We investigated individual cells in M. leprae-inoculated nude mice by single-cell RNA-sequencing (scRNA-seq). For macrophages, we dissected two M1-like subsets and five M2-like subsets, where lipid-associated signatures were pervasive in both M1-like and M2-like subsets. There were four macrophage trajectories showing: (i) pro-inflammatory (M1), (ii) lipid metabolism-related (M2), (iii) anti-inflammatory (M2), and (iv) interferon-stimulated gene-related (M2) fates. They displayed early divergence without ever rejoining along the paths, suggesting simultaneous or continuous stimuli for macrophage activation in leprosy. The scRNA-seq predicted Schwann cell-macrophage interactions (Notch1-Jag1, Plxnb1-Sema4d interactions). An immature Schwann cell subset showing Tfap2a expression was identified, indicating Schwann cell dedifferentiation in leprosy tissues. Expressions of Notch1, Jag1, Plxnb1, Sema4d, and Tfap2a were validated in mouse or human leprosy tissues by immunohistochemistry. We identified both pro-inflammatory and inflammation-resolution signatures, where lipid-associated signatures were pervasive to the macrophages, representing leprosy-specific macrophage states for prolonged and repeated episodes of inflammation and resolution. Our study identified refined molecular states and interactions of macrophages and Schwann cells, suggesting novel insights into the pathogenesis of unhealed inflammation with neuropathy and potential therapeutic targets for leprosy.

The Role of Sensory Innervation in Homeostatic and Injury-Induced Corneal Epithelial Renewal

The cornea is the window through which we see the world. Corneal clarity is required for vision, and blindness occurs when the cornea becomes opaque. The cornea is covered by unique transparent epithelial cells that serve as an outermost cellular barrier bordering between the cornea and the external environment. Corneal sensory nerves protect the cornea from injury by triggering tearing and blink reflexes, and are also thought to regulate corneal epithelial renewal via unknown mechanism(s). When protective corneal sensory innervation is absent due to infection, trauma, intracranial tumors, surgery, or congenital causes, permanent blindness results from repetitive epithelial microtraumas and failure to heal. The condition is termed neurotrophic keratopathy (NK), with an incidence of 5:10,000 people worldwide. In this report, we review the currently available therapeutic solutions for NK and discuss the progress in our understanding of how the sensory nerves induce corneal epithelial renewal.

Dry eye disease in mice activates adaptive corneal epithelial regeneration distinct from constitutive renewal in homeostasis

Many epithelial compartments undergo constitutive renewal in homeostasis but activate unique regenerative responses following injury. The clear corneal epithelium is crucial for vision and is renewed from limbal stem cells (LSCs). Using single-cell RNA sequencing, we profiled the mouse corneal epithelium in homeostasis, aging, diabetes, and dry eye disease (DED), where tear deficiency predisposes the cornea to recurrent injury. In homeostasis, we capture the transcriptional states that accomplish continuous tissue turnover. We leverage our dataset to identify candidate genes and gene networks that characterize key stages across homeostatic renewal, including markers for LSCs. In aging and diabetes, there were only mild changes with <15 dysregulated genes. The constitutive cell types that accomplish homeostatic renewal were conserved in DED but were associated with activation of cell states that comprise "adaptive regeneration." We provide global markers that distinguish cell types in homeostatic renewal vs. adaptive regeneration and markers that specifically define DED-elicited proliferating and differentiating cell types. We validate that expression of SPARC, a marker of adaptive regeneration, is also induced in corneal epithelial wound healing and accelerates wound closure in a corneal epithelial cell scratch assay. Finally, we propose a classification system for LSC markers based on their expression fidelity in homeostasis and disease. This transcriptional dissection uncovers the dramatically altered transcriptional landscape of the corneal epithelium in DED, providing a framework and atlas for future study of these ocular surface stem cells in health and disease.

The transcriptional profile of keloidal Schwann cells

Recently, a specific Schwann cell type with profibrotic and tissue regenerative properties that contributes to keloid formation has been identified. In the present study, we reanalyzed published single-cell RNA sequencing (scRNA-seq) studies of keloids, healthy skin, and normal scars to reliably determine the specific gene expression profile of keloid-specific Schwann cell types in more detail. We were able to confirm the presence of the repair-like, profibrotic Schwann cell type in the datasets of all three studies and identified a specific gene-set for these Schwann cells. In contrast to keloids, in normal scars, the number of Schwann cells was not increased, nor was their gene expression profile distinctly different from that of Schwann cells of normal skin. In addition, our bioinformatics analysis provided evidence for a role of transcription factors of the AP1, STAT, and KLF families, and members of the IER genes in the dedifferentiation process of keloidal Schwann cells. Together, our analysis strengthens the role of the profibrotic Schwann cell type in the formation of keloids. Knowledge of the exact gene expression profile of these Schwann cells will facilitate their identification in other organs and diseases.

Single-Cell RNA-Seq Reveals the Heterogeneity of Cell Communications between Schwann Cells and Fibroblasts within the Microenvironment in Vestibular Schwannoma

Vestibular schwannomas (VSs), which develop from Schwann cells (SCs) of the vestibular nerve, are the most prevalent benign tumors of the cerebellopontine angle and internal auditory canal. Despite advances in treatment, the cellular components and mechanisms of VS tumor progression remain unclear. Here, we performed single-cell RNA-sequencing on clinically surgically isolated VS samples and clarified their cellular composition in a sophisticated manner, including the heterogeneous SC subtypes. Advanced bioinformatics analysis also revealed the associated biological functions, pseudotime trajectory, and transcriptional network of the SC subgroups. We also found that there was tight intercellular communication between SCs and tumor-associated fibroblasts via integrin and growth factor signaling and that the gene expression differences in SCs and fibroblasts determined the heterogeneity of cellular communication in different individuals. Our findings suggest a microenvironmental mechanism underlying the development of VS.

Nimodipine Exerts Beneficial Effects on the Rat Oligodendrocyte Cell Line OLN-93

Multiple sclerosis (MS) is a chronic autoimmune disease of the central nervous system (CNS). Therapy is currently limited to drugs that interfere with the immune system; treatment options that primarily mediate neuroprotection and prevent neurodegeneration are not available. Here, we studied the effects of nimodipine on the rat cell line OLN-93, which resembles young mature oligodendrocytes. Nimodipine is a dihydropyridine that blocks the voltage-gated L-type calcium channel family members Cav1.2 and Cav1.3. Our data show that the treatment of OLN-93 cells with nimodipine induced the upregulation of myelin genes, in particular of proteolipid protein 1 (Plp1), which was confirmed by a significantly greater expression of PLP1 in immunofluorescence analysis and the presence of myelin structures in the cytoplasm at the ultrastructural level. Whole-genome RNA sequencing additionally revealed the upregulation of genes that are involved in neuroprotection, remyelination, and antioxidation pathways. Interestingly, the observed effects were independent of Cav1.2 and Cav1.3 because OLN-93 cells do not express these channels, and there was no measurable response pattern in patch-clamp analysis. Taking into consideration previous studies that demonstrated a beneficial effect of nimodipine on microglia, our data support the notion that nimodipine is an interesting drug candidate for the treatment of MS and other demyelinating diseases.

A single cell atlas of human cornea that defines its development, limbal progenitor cells and their interactions with the immune cells

Purpose

Single cell (sc) analyses of key embryonic, fetal and adult stages were performed to generate a comprehensive single cell atlas of all the corneal and adjacent conjunctival cell types from development to adulthood.

Methods

Four human adult and seventeen embryonic and fetal corneas from 10 to 21 post conception week (PCW) specimens were dissociated to single cells and subjected to scRNA- and/or ATAC-Seq using the 10x Genomics platform. These were embedded using Uniform Manifold Approximation and Projection (UMAP) and clustered using Seurat graph-based clustering. Cluster identification was performed based on marker gene expression, bioinformatic data mining and immunofluorescence (IF) analysis. RNA interference, IF, colony forming efficiency and clonal assays were performed on cultured limbal epithelial cells (LECs).

Results

scRNA-Seq analysis of 21,343 cells from four adult human corneas and adjacent conjunctivas revealed the presence of 21 cell clusters, representing the progenitor and differentiated cells in all layers of cornea and conjunctiva as well as immune cells, melanocytes, fibroblasts, and blood/lymphatic vessels. A small cell cluster with high expression of limbal progenitor cell (LPC) markers was identified and shown via pseudotime analysis to give rise to five other cell types representing all the subtypes of differentiated limbal and corneal epithelial cells. A novel putative LPCs surface marker, GPHA2, expressed on the surface of 0.41% ± 0.21 of the cultured LECs, was identified, based on predominant expression in the limbal crypts of adult and developing cornea and RNAi validation in cultured LECs. Combining scRNA- and ATAC-Seq analyses, we identified multiple upstream regulators for LPCs and demonstrated a close interaction between the immune cells and limbal progenitor cells. RNA-Seq analysis indicated the loss of GPHA2 expression and acquisition of proliferative limbal basal epithelial cell markers during ex vivo LEC expansion, independently of the culture method used. Extending the single cell analyses to keratoconus, we were able to reveal activation of collagenase in the corneal stroma and a reduced pool of limbal suprabasal cells as two key changes underlying the disease phenotype. Single cell RNA-Seq of 89,897 cells obtained from embryonic and fetal cornea indicated that during development, the conjunctival epithelium is the first to be specified from the ocular surface epithelium, followed by the corneal epithelium and the establishment of LPCs, which predate the formation of limbal niche by a few weeks.

Conclusions

Our scRNA-and ATAC-Seq data of developing and adult cornea in steady state and disease conditions provide a unique resource for defining genes/pathways that can lead to improvement in ex vivo LPCs expansion, stem cell differentiation methods and better understanding and treatment of ocular surface disorders.

Corneal nonmyelinating Schwann cells illuminated by single-cell transcriptomics and visualized by protein biomarkers

The cornea is the most innervated tissue in the human body. Myelinated axons upon inserting into the peripheral corneal stroma lose their myelin sheaths and continue into the central cornea wrapped by only nonmyelinating corneal Schwann cells (nm-cSCs). This anatomical organization is believed to be important for central vision. Here we employed single-cell RNA sequencing (scRNA-seq), microscopy, and transgenics to characterize these nm-cSCs of the central cornea. Using principal component analysis, uniform manifold approximation and projection, and unsupervised hierarchal cell clustering of scRNA-seq data derived from central corneal cells of male rabbits, we successfully identified several clusters representing different corneal cell types, including a unique cell cluster representing nm-cSCs. To confirm protein expression of cSC genes, we performed cross-species validation, employing corneal whole-mount immunostaining with confocal microscopy in mouse corneas. The expression of several representative proteins of nm-cSCs were validated. As the proteolipid protein 1 (PLP1) gene was also expressed in nm-cSCs, we explored the Plp1-eGFP transgenic reporter mouse line to visualize cSCs. Specific and efficient eGFP expression was observed in cSCs in adult mice of different ages. Of several putative cornea-specific SC genes identified, Dickkopf-related protein 1 was shown to be present in nm-cSCs. Taken together, our findings, for the first time, identify important insights and tools toward the study nm-cSCs in isolated tissue and adult animals. We expect that our results will advance the future study of nm-cSCs in applications of nerve repair, and provide a resource for the study of corneal sensory function.