Search

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.

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.

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.

Severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) first emerged in December 2019 and spread quickly causing the coronavirus disease 2019 (COVID‐19) pandemic. Recent single cell RNA‐Seq analyses have shown the presence of SARS‐CoV‐2 entry factors in the human corneal, limbal, and conjunctival superficial epithelium, leading to suggestions that the human ocular surface may serve as an additional entry gateway and infection hub for SARS‐CoV‐2. In this article, we review the ocular clinical presentations of COVID‐19 and the features of the ocular surface that may underline the overall low ocular SARS‐CoV‐2 infection. We critically evaluate the studies performed in nonhuman primates, ex vivo organ culture ocular models, stem cell derived eye organoids and the differences in infection efficiency observed in different parts of human ocular surface epithelium. Finally, we highlight the additional work that needs to be carried out to understand the immune response of the ocular surface to SARS‐CoV‐2 infection, which can be translated into prophylactic treatments that may be applied to other organ systems., Schematic showing potential SARS‐CoV‐2 entry into the eye via infection of ocular surface epithelium (left hand side panel) or the systemic routes (right hand side panel).

Background & Aims Serum microRNA (miRNA) levels are known to change in non-alcoholic fatty liver disease (NAFLD) and may serve as useful biomarkers. This study aimed to profile miRNAs comprehensively at all NAFLD stages. Methods We profiled 2,083 serum miRNAs in a discovery cohort (183 cases with NAFLD representing the complete NAFLD spectrum and 10 population controls). miRNA libraries generated by HTG EdgeSeq were sequenced by Illumina NextSeq. Selected serum miRNAs were profiled in 372 additional cases with NAFLD and 15 population controls by quantitative reverse transcriptase PCR. Results Levels of 275 miRNAs differed between cases and population controls. Fewer differences were seen within individual NAFLD stages, but miR-193a-5p consistently showed increased levels in all comparisons. Relative to NAFL/non-alcoholic steatohepatitis (NASH) with mild fibrosis (stage 0/1), 3 miRNAs (miR-193a-5p, miR-378d, and miR378d) were increased in cases with NASH and clinically significant fibrosis (stages 2–4), 7 (miR193a-5p, miR-378d, miR-378e, miR-320b, miR-320c, miR-320d, and miR-320e) increased in cases with NAFLD activity score (NAS) 5–8 compared with lower NAS, and 3 (miR-193a-5p, miR-378d, and miR-378e) increased but 1 (miR-19b-3p) decreased in steatosis, activity, and fibrosis (SAF) activity score 2–4 compared with lower SAF activity. The significant findings for miR-193a-5p were replicated in the additional cohort with NAFLD. Studies in Hep G2 cells showed that following palmitic acid treatment, miR-193a-5p expression decreased significantly. Gene targets for miR-193a-5p were investigated in liver RNAseq data for a case subgroup (n = 80); liver GPX8 levels correlated positively with serum miR-193a-5p. Conclusions Serum miR-193a-5p levels correlate strongly with NAFLD activity grade and fibrosis stage. MiR-193a-5p may have a role in the hepatic response to oxidative stress and is a potential clinically tractable circulating biomarker for progressive NAFLD. Lay summary MicroRNAs (miRNAs) are small pieces of nucleic acid that may turn expression of genes on or off. These molecules can be detected in the blood circulation, and their levels in blood may change in liver disease including non-alcoholic fatty liver disease (NAFLD). To see if we could detect specific miRNA associated with advanced stages of NAFLD, we carried out miRNA sequencing in a group of 183 patients with NAFLD of varying severity together with 10 population controls. We found that a number of miRNAs showed changes, mainly increases, in serum levels but that 1 particular miRNA miR-193a-5p consistently increased. We confirmed this increase in a second group of cases with NAFLD. Measuring this miRNA in a blood sample may be a useful way to determine whether a patient has advanced NAFLD without an invasive liver biopsy., Graphical abstract, Highlights • Serum miRNA was sequenced in 183 NAFLD cases of varying severity and 10 population controls. • Plasma levels of miR-193a-5p were significantly increased in patients with advanced fibrosis, high NAS scores, or high SAF scores. • Other miRNAs including miR378d and miR378e were also significantly increased in certain comparisons. • The findings for miR-193a-5p were replicated in a cohort of 372 additional NAFLD cases.

Single-cell sequencing technologies have emerged as a revolutionary tool with transformative new methods to profile genetic, epigenetic, spatial, and lineage information in individual cells. Single-cell RNA sequencing (scRNA-Seq) allows researchers to collect large datasets detailing the transcriptomes of individual cells in space and time and is increasingly being applied to reveal cellular heterogeneity in retinal development, normal physiology, and disease, and provide new insights into cell-type specific markers and signaling pathways. In recent years, scRNA-Seq datasets have been generated from retinal tissue and pluripotent stem cell-derived retinal organoids. Their cross-comparison enables staging of retinal organoids, identification of specific cells in developing and adult human neural retina and provides deeper insights into cell-type sub-specification and geographical differences. In this article, we review the recent rapid progress in scRNA-Seq analyses of retina and retinal organoids, the questions that remain unanswered and the technical challenges that need to be overcome to achieve consistent results that reflect the complexity, functionality, and interactions of all retinal cell types.

Although tropism of SARS-CoV-2 for respiratory tract epithelial cells is well established, an open question is whether the conjunctival epithelium is also a target for SARS-CoV-2. Conjunctival epithelial cells, which express viral entry receptors ACE2 and 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 progenitor, basal and superficial epithelial cells and fibroblasts, which could be maintained successfully up to day 75 of differentiation. Using single-cell RNA Seq, with complementary imaging and virological assays, we observed that while all conjunctival cell types were permissive to 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-Kβ 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 transplants.

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). A single-cell atlas of human fetal bone marrow in healthy fetuses and fetuses with Down syndrome provides insight into developmental haematopoiesis in humans and the transcription and functional differences that occur in Down syndrome.

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.

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.

Throughout postnatal life, haematopoiesis in the bone marrow (BM) maintains blood and immune cell production. Haematopoiesis first emerges in human BM at 12 post conception weeks while fetal liver (FL) haematopoiesis is still expanding. Yet, almost nothing is known about how fetal BM evolves to meet the highly specialised needs of the fetus and newborn infant. Here, we detail the development of fetal BM including stroma using single cell RNA-sequencing. We find that the full blood and immune cell repertoire is established in fetal BM in a short time window of 6-7 weeks early in the second trimester. Fetal BM promotes rapid and extensive diversification of myeloid cells, with granulocytes, eosinophils and dendritic cell (DC) subsets emerging for the first time. B-lymphocyte expansion occurs, in contrast with erythroid predominance in FL at the same gestational age. We identify transcriptional and functional differences that underlie tissue-specific identity and cellular diversification in fetal BM and FL. Finally, we reveal selective disruption of B-lymphocyte, erythroid and myeloid development due to cell intrinsic differentiation bias as well as extrinsic regulation through an altered microenvironment in the fetal BM from constitutional chromosome anomaly Down syndrome during this crucial developmental time window.

Retinoblastoma (Rb) is a childhood cancer of the developing retina in response to bi-allelic inactivation of RB1 or MYCN amplification, accounting for up to 17% of all tumours in infancy. To gain insights into the transcriptional events of cell state transitions during Rb development, we developed two models via retinal organoid differentiation of a RB1 depleted human embryonic stem cell line (RB1-null hESCs) and a RB1 patient-specific induced pluripotent (iPSC) line harbouring a RB1 biallelic mutation (c.2082delC). Both models were characterised by RB1 depletion and a significant increase in the fraction of proliferating cone precursors (RXRγ+Ki67+), which were defined as the cell of origin for Rb by single cell RNA-Seq analysis. The RB1 depleted retinal organoids displayed similar features to Rb tumours, 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. The patient-specific iPSC model displayed an enhanced reduction of amacrine, horizontal and retinal ganglion cells and an accelerated loss of cone photoreceptor markers during transition towards Rb, compared to the RB1-null hESC model. In both models, the Rb cells of origin, intermediate retinoma and/or Rb cells expressed retinal ganglion and horizontal cell markers in addition to cone markers, a novel finding, which could help to better characterise these tumours with possible therapeutic implications. Application of Melphalan, Topotecan and TW-37 led to a significant reduction in the fraction of proliferating cone precursors, validating the suitability of these in vitro models for testing novel therapeutics for Rb.

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., Graphical abstract Schematic presentation of main techniques and findings presented in this manuscript.Image 1

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.

International audience; The mechanisms that drive nonalcoholic fatty liver disease (NAFLD) remain incompletely understood. This large multicenter study characterized the transcriptional changes that occur in liver tissue across the NAFLD spectrum as disease progresses to cirrhosis to identify potential circulating markers. We performed high-throughput RNA sequencing on a discovery cohort comprising histologically characterized NAFLD samples from 206 patients. Unsupervised clustering stratified NAFLD on the basis of disease activity and fibrosis stage with differences in age, aspartate aminotransferase (AST), type 2 diabetes mellitus, and carriage of PNPLA3 rs738409 , a genetic variant associated with NAFLD. Relative to early disease, we consistently identified 25 differentially expressed genes as fibrosing steatohepatitis progressed through stages F2 to F4. This 25-gene signature was independently validated by logistic modeling in a separate replication cohort ( n = 175), and an integrative analysis with publicly available single-cell RNA sequencing data elucidated the likely relative contribution of specific intrahepatic cell populations. Translating these findings to the protein level, SomaScan analysis in more than 300 NAFLD serum samples confirmed that circulating concentrations of proteins AKR1B10 and GDF15 were strongly associated with disease activity and fibrosis stage. Supporting the biological plausibility of these data, in vitro functional studies determined that endoplasmic reticulum stress up-regulated expression of AKR1B10 , GDF15 , and PDGFA , whereas GDF15 supplementation tempered the inflammatory response in macrophages upon lipid loading and lipopolysaccharide stimulation. This study provides insights into the pathophysiology of progressive fibrosing steatohepatitis, and proof of principle that transcriptomic changes represent potentially tractable and clinically relevant markers of disease progression.

Summary 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., Graphical Abstract, Highlights • Distinct development trajectories of DC2 and DC3 underpin human cDC2 heterogeneity • pDC, cDC1, and DC2 (classical DCs) develop from LMPPs along a CD123+ IRF8high pathway • DC3 and monocytes develop from CD33+ GMPs along an IRF8low SIRPA+ pathway • IRF8 deficiency causes gene dose-dependent loss of IRF8high then IRF8low pathway DCs, Heterogeneity of human CD1c+ dendritic cells (cDC2s) is described, but how this arises is unknown. Cytlak and colleagues demonstrate that the cDC2 subsets, DC2 and DC3, develop along distinct hematopoietic trajectories, defined by differential IRF8 expression. DC2s develop from LMPPs along an IRF8hi pathway, while DC3 differentiation follows an IRF8low trajectory.

SummaryTo study the development and composition of human ocular surface, we performed single cell (sc) RNA-Seq at key embryonic, fetal and adult stages and generated the first atlas of the corneal cell types from development to adulthood. Our data indicate 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 proliferative epithelial progenitors, which predate the formation of limbal niche by a few weeks. Bioinformatic comparison of adult cell clusters identified GPHA2, a novel cell-surface marker for quiescent limbal stem cells (qLSCs), whose function is to maintain qLSCs self-renewal. Combining scRNA- and ATAC-Seq analysis, we identified multiple upstream regulators for qLSCs and transit amplifying (TA) cells and demonstrated a close interaction between the immune cells and epithelial stem and progenitor cells in the cornea. RNA-Seq analysis indicated loss of qLSCs and acquisition of proliferative limbal basal epithelial progenitor markers duringex vivolimbal epithelial cell 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 TA cells in the limbal epithelium as two key changes underlying the disease phenotype. Our scRNA- and ATAC-Seq data of developing and adult cornea in steady state and disease conditions provide a unique resource for defining pathways/genes that can lead to improvement inex vivoexpansion and differentiation methods for cell based replacement therapies and better understanding and treatment of ocular surface disorders.Key findingsscRNA-Seq of adult human cornea and conjunctiva reveals the signature of various ocular surface cell populationsscRNA-Seq of human developing cornea identifies stage-specific definitions of corneal epithelial, stromal and endothelial layersscRNA-Seq analysis results in identification of novel markers for qLSCs and TA cellsCombined scRNA- and ATAC-Seq analysis reveals key transcriptional networks in qLSCs and TA cells and close interactions with immune cellsExpansion of limbal epithelium results in downregulation of qLSCs and acquisition of proliferative limbal epithelial progenitor markersscRNA-Seq of keratoconus corneas reveals activation of collagenase in the corneal stroma and a reduced pool of TA cells in the limbal epitheliumGraphical abstractSchematic presentation of main techniques and findings presented in this manuscript.

Purpose The high infection rate of SARS-CoV-2 necessitates the need for multiple studies identifying the molecular mechanisms that facilitate the viral entry and propagation. Currently the potential extra-respiratory transmission routes of SARS-CoV-2 remain unclear. Methods Using single-cell RNA Seq and ATAC-Seq datasets and immunohistochemical analysis, we investigated SARS-CoV-2 tropism in the embryonic, fetal and adult human ocular surface. Results The co-expression of ACE2 receptor and entry protease TMPRSS2 was detected in the human adult conjunctival, limbal and corneal epithelium, but not in the embryonic and fetal ocular surface up to 21 post conception weeks. These expression patterns were corroborated by the single cell ATAC-Seq data, which revealed a permissive chromatin in ACE2 and TMPRSS2 loci in the adult conjunctival, limbal and corneal epithelium. Co-expression of ACE2 and TMPRSS2 was strongly detected in the superficial limbal, corneal and conjunctival epithelium, implicating these as target entry cells for SARS-CoV-2 in the ocular surface. Strikingly, we also identified the key pro-inflammatory signals TNF, NFKβ and IFNG as upstream regulators of the transcriptional profile of ACE2+TMPRSS2+ cells in the superficial conjunctival epithelium, suggesting that SARS-CoV-2 may utilise inflammatory driven upregulation of ACE2 and TMPRSS2 expression to enhance infection in ocular surface. Conclusions Together our data indicate that the human ocular surface epithelium provides an additional entry portal for SARS-CoV-2, which may exploit inflammatory driven upregulation of ACE2 and TMPRSS2 entry factors to enhance infection.

The COVID-19 pandemic, caused by the novel coronavirus SARS-CoV-2, creates an urgent need for identifying molecular mechanisms that mediate viral entry, propagation, and tissue pathology. Cell membrane bound angiotensin-converting enzyme 2 (ACE2) and associated proteases, transmembrane protease serine 2 (TMPRSS2) and Cathepsin L (CTSL), were previously identified as mediators of SARS-CoV2 cellular entry. Here, we assess the cell type-specific RNA expression of ACE2, TMPRSS2, and CTSL through an integrated analysis of 107 single-cell and single-nucleus RNA-Seq studies, including 22 lung and airways datasets (16 unpublished), and 85 datasets from other diverse organs. Joint expression of ACE2 and the accessory proteases identifies specific subsets of respiratory epithelial cells as putative targets of viral infection in the nasal passages, airways, and alveoli. Cells that co-express ACE2 and proteases are also identified in cells from other organs, some of which have been associated with COVID-19 transmission or pathology, including gut enterocytes, corneal epithelial cells, cardiomyocytes, heart pericytes, olfactory sustentacular cells, and renal epithelial cells. Performing the first meta-analyses of scRNA-seq studies, we analyzed 1,176,683 cells from 282 nasal, airway, and lung parenchyma samples from 164 donors spanning fetal, childhood, adult, and elderly age groups, associate increased levels of ACE2, TMPRSS2, and CTSL in specific cell types with increasing age, male gender, and smoking, all of which are epidemiologically linked to COVID-19 susceptibility and outcomes. Notably, there was a particularly low expression of ACE2 in the few young pediatric samples in the analysis. Further analysis reveals a gene expression program shared by ACE2+TMPRSS2+ cells in nasal, lung and gut tissues, including genes that may mediate viral entry, subtend key immune functions, and mediate epithelial-macrophage cross-talk. Amongst these are IL6, its receptor and co-receptor, IL1R, TNF response pathways, and complement genes. Cell type specificity in the lung and airways and smoking effects were conserved in mice. Our analyses suggest that differences in the cell type-specific expression of mediators of SARS-CoV-2 viral entry may be responsible for aspects of COVID-19 epidemiology and clinical course, and point to putative molecular pathways involved in disease susceptibility and pathogenesis.

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. Electronic supplementary material The online version of this article (10.1007/s00239-018-9855-7) contains supplementary material, which is available to authorized users.

Spatial Transcriptomics allows the sequencing of the complete transcriptomes from barcoded regions of intact tissue. The technology has the potential to answer a wide range of biological questions concerning cellular function, but analysis of the data presents a number of challenges which are not met by existing analysis tools. Here we present Spaniel, an R package providing a framework for analysing and sharing Spatial Transcriptomics data.

Mitochondrial disorders are heterogeneous, showing variable presentation and penetrance. Over the last three decades, our ability to recognize mitochondrial patients and diagnose these mutations, linking genotype to phenotype, has greatly improved. However, it has become increasingly clear that these strides in diagnostics have not benefited all population groups. Recent studies have demonstrated that patients from genetically understudied populations, in particular those of black African heritage, are less likely to receive a diagnosis of mtDNA disease. It has been suggested that haplogroup context might influence the presentation and penetrance of mtDNA disease; thus, the spectrum of mutations that are associated with disease in different populations. However, to date there is only one well‐established example of such an effect: the increased penetrance of two Leber's hereditary optic neuropathy mutations on a haplogroup J background. This paper conducted the most extensive investigation to date into the importance of haplogroup context on the pathogenicity of mtDNA mutations. We searched for proven human point mutations across 726 multiple sequence alignments derived from 33 non‐human species absent of disease. A total of 58 pathogenic point mutations arise in the sequences of these species. We assessed the sequence context and found evidence of population variants that could modulate the phenotypic expression of these point mutations masking the pathogenic effects seen in humans. This supports the theory that sequence context is influential in the presentation of mtDNA disease and has implications for diagnostic practices. We have shown that our current understanding of the pathogenicity of mtDNA point mutations, primarily built on studies of individuals with haplogroups HVUKTJ, will not present a complete picture. This will have the effect of creating a diagnostic inequality, whereby individuals who do not belong to these lineages are less likely to receive a genetic diagnosis.

The formation of human dendritic cells (DC) is controlled by multiple hematopoietic transcription factors, including IRF8. Loss of IRF8 differentially affects plasmacytoid DC (pDC) and the classical DC lineages, cDC1 and cDC2. Additionally, cDC2 are heterogeneous, comprising DC-like DC2 and monocyte-like DC3. Using high dimensional analysis, in vitro differentiation, and an allelic series of human IRF8 deficiency, we show that CD1c+DC heterogeneity originates from two distinct pathways of DC development. The lymphoid-primed IRF8high pathway carries pDC, cDC1, and DC2 trajectories, while the common myeloid IRF8low pathway forms DC3 and monocytes. We trace distinct trajectories of cDC1, DC2 and DC3 through the GMP showing that AXL+SIGLEC6+ pre-DC map exclusively to the DC2 pathway. The preservation of human CD1c+DC in partial IRF8 deficiency is explained by the replacement of DC2 by an expanded DC3 population. These observations are congruent with lineage-primed models of hematopoiesis, illustrating new links between immune cell function and ontogeny.

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

Mitochondrial DNA (mtDNA) mutations are well recognized as an important cause of inherited disease. Diseases caused by mtDNA mutations exhibit a high degree of clinical heterogeneity with a complex genotype-phenotype relationship, with many such mutations exhibiting incomplete penetrance. There is evidence that the spectrum of mutations causing mitochondrial disease might differ between different mitochondrial lineages (haplogroups) seen in different global populations. This would point to the importance of sequence context in the expression of mutations. To explore this possibility, we looked for mutations which are known to cause disease in humans, in animals of other species unaffected by mtDNA disease. The mt-tRNA genes are the location of many pathogenic mutations, with the m.3243A>G mutation on the mt-tRNA-Leu(UUR) being the most frequently seen mutation in humans. This study looked for the presence of m.3243A>G in 2784 sequences from 33 species, as well as any of the other mutations reported in association with disease located on mt-tRNA-Leu(UUR). We report a number of disease associated variations found on mt-tRNA-Leu(UUR) in other chordates, as the major population variant, with m.3243A>G being seen in 6 species. In these, we also found a number of mutations which appear compensatory and which could prevent the pathogenicity associated with this change in humans. This work has important implications for the discovery and diagnosis of mtDNA mutations in non-European populations. In addition, it might provide a partial explanation for the conflicting results in the literature that examines the role of mtDNA variants in complex traits.

Angiotensin-converting enzyme 2 (ACE2) and accessory proteases (TMPRSS2 and CTSL) are needed for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) cellular entry, and their expression may shed light on viral tropism and impact across the body. We assessed the cell-type-specific expression of ACE2, TMPRSS2 and CTSL across 107 single-cell RNA-sequencing studies from different tissues. ACE2, TMPRSS2 and CTSL are coexpressed in specific subsets of respiratory epithelial cells in the nasal passages, airways and alveoli, and in cells from other organs associated with coronavirus disease 2019 (COVID-19) transmission or pathology. We performed a meta-analysis of 31 lung single-cell RNA-sequencing studies with 1,320,896 cells from 377 nasal, airway and lung parenchyma samples from 228 individuals. This revealed cell-type-specific associations of age, sex and smoking with expression levels of ACE2, TMPRSS2 and CTSL. Expression of entry factors increased with age and in males, including in airway secretory cells and alveolar type 2 cells. Expression programs shared by ACE2+TMPRSS2+ cells in nasal, lung and gut tissues included genes that may mediate viral entry, key immune functions and epithelial-macrophage cross-talk, such as genes involved in the interleukin-6, interleukin-1, tumor necrosis factor and complement pathways. Cell-type-specific expression patterns may contribute to the pathogenesis of COVID-19, and our work highlights putative molecular pathways for therapeutic intervention.