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.

Single-cell insights into immune dysregulation in rheumatoid arthritis flare versus drug-free remission.

Immune-mediated inflammatory diseases (IMIDs) are typically characterised by relapsing and remitting flares of inflammation. However, the unpredictability of disease flares impedes their study. Addressing this critical knowledge gap, we use the experimental medicine approach of immunomodulatory drug withdrawal in rheumatoid arthritis (RA) remission to synchronise flare processes allowing detailed characterisation. Exploratory mass cytometry analyses reveal three circulating cellular subsets heralding the onset of arthritis flare - CD45RO+PD1hi CD4+ and CD8+ T cells, and CD27+CD86+CD21- B cells - further characterised by single-cell sequencing. Distinct lymphocyte subsets including cytotoxic and exhausted CD4+ memory T cells, memory CD8+CXCR5+ T cells, and IGHA1+ plasma cells are primed for activation in flare patients. Regulatory memory CD4+ T cells (Treg cells) increase at flare onset, but with dysfunctional regulatory marker expression compared to drug-free remission. Significant clonal expansion is observed in T cells, but not B cells, after drug cessation; this is widespread throughout memory CD8+ T cell subsets but limited to the granzyme-expressing cytotoxic subset within CD4+ memory T cells. Based on our observations, we suggest a model of immune dysregulation for understanding RA flare, with potential for further translational research towards novel avenues for its treatment and prevention.

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

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

Release of Histone H3K4-reading transcription factors from chromosomes in mitosis is independent of adjacent H3 phosphorylation.

Histone modifications influence the recruitment of reader proteins to chromosomes to regulate events including transcription and cell division. The idea of a histone code, where combinations of modifications specify unique downstream functions, is widely accepted and can be demonstrated in vitro. For example, on synthetic peptides, phosphorylation of Histone H3 at threonine-3 (H3T3ph) prevents the binding of reader proteins that recognize trimethylation of the adjacent lysine-4 (H3K4me3), including the TAF3 component of TFIID. To study these combinatorial effects in cells, we analyzed the genome-wide distribution of H3T3ph and H3K4me2/3 during mitosis. We find that H3T3ph anti-correlates with adjacent H3K4me2/3 in cells, and that the PHD domain of TAF3 can bind H3K4me2/3 in isolated mitotic chromatin despite the presence of H3T3ph. Unlike in vitro, H3K4 readers are still displaced from chromosomes in mitosis in Haspin-depleted cells lacking H3T3ph. H3T3ph is therefore unlikely to be responsible for transcriptional downregulation during cell division.

Programmed cell death-1 receptor-mediated regulation of Tbet+NK1.1- innate lymphoid cells within the tumor microenvironment.

Innate lymphoid cells (ILCs) play a key role in tissue-mediated immunity and can be controlled by coreceptor signaling. Here, we define a subset of ILCs that are Tbet+NK1.1- and are present within the tumor microenvironment (TME). We show programmed death-1 receptor (PD-1) expression on ILCs within TME is found in Tbet+NK1.1- ILCs. PD-1 significantly controlled the proliferation and function of Tbet+NK1.1- ILCs in multiple murine and human tumors. We found tumor-derived lactate enhanced PD-1 expression on Tbet+NK1.1- ILCs within the TME, which resulted in dampened the mammalian target of rapamycin (mTOR) signaling along with increased fatty acid uptake. In line with these metabolic changes, PD-1-deficient Tbet+NK1.1- ILCs expressed significantly increased IFNγ and granzyme B and K. Furthermore, PD-1-deficient Tbet+NK1.1- ILCs contributed toward diminished tumor growth in an experimental murine model of melanoma. These data demonstrate that PD-1 can regulate antitumor responses of Tbet+NK1.1- ILCs within the TME.

Tracheostomy in children is associated with neutrophilic airway inflammation.

BACKGROUND: Tracheostomies in children are associated with significant morbidity, poor quality of life, excess healthcare costs and excess mortality. The underlying mechanisms facilitating adverse respiratory outcomes in tracheostomised children are poorly understood. We aimed to characterise airway host defence in tracheostomised children using serial molecular analyses. METHODS: Tracheal aspirates, tracheal cytology brushings and nasal swabs were prospectively collected from children with a tracheostomy and controls. Transcriptomic, proteomic and metabolomic methods were applied to characterise the impact of tracheostomy on host immune response and the airway microbiome. RESULTS: Children followed up serially from the time of tracheostomy up to 3 months postprocedure (n=9) were studied. A further cohort of children with a long-term tracheostomy were also enrolled (n=24). Controls (n=13) comprised children without a tracheostomy undergoing bronchoscopy. Long-term tracheostomy was associated with airway neutrophilic inflammation, superoxide production and evidence of proteolysis when compared with controls. Reduced airway microbial diversity was established pre-tracheostomy and sustained thereafter. CONCLUSIONS: Long-term childhood tracheostomy is associated with a inflammatory tracheal phenotype characterised by neutrophilic inflammation and the ongoing presence of potential respiratory pathogens. These findings suggest neutrophil recruitment and activation as potential exploratory targets in seeking to prevent recurrent airway complications in this vulnerable group of patients.

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

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

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

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

Regulation of CHK1 inhibitor resistance by a c-Rel and USP1 dependent pathway.

Previously, we discovered that deletion of c-Rel in the Eµ-Myc mouse model of lymphoma results in earlier onset of disease, a finding that contrasted with the expected function of this NF-κB subunit in B-cell malignancies. Here we report that Eµ-Myc/cRel-/- cells have an unexpected and major defect in the CHK1 pathway. Total and phospho proteomic analysis revealed that Eµ-Myc/cRel-/- lymphomas highly resemble wild-type (WT) Eµ-Myc lymphomas treated with an acute dose of the CHK1 inhibitor (CHK1i) CCT244747. Further analysis demonstrated that this is a consequence of Eµ-Myc/cRel-/- lymphomas having lost expression of CHK1 protein itself, an effect that also results in resistance to CCT244747 treatment in vivo. Similar down-regulation of CHK1 protein levels was also seen in CHK1i resistant U2OS osteosarcoma and Huh7 hepatocellular carcinoma cells. Further investigation revealed that the deubiquitinase USP1 regulates CHK1 proteolytic degradation and that its down-regulation in our model systems is responsible, at least in part, for these effects. We demonstrate that treating WT Eµ-Myc lymphoma cells with the USP1 inhibitor ML323 was highly effective at reducing tumour burden in vivo. Targeting USP1 activity may thus be an alternative therapeutic strategy in MYC-driven tumours.

Up-regulation of the PI3K/AKT and RHO/RAC/PAK signalling pathways in CHK1 inhibitor resistant Eµ-Myc lymphoma cells.

The development of resistance and the activation of bypass pathway signalling represents a major problem for the clinical application of protein kinase inhibitors. While investigating the effect of either a c-Rel deletion or RelAT505A phosphosite knockin on the Eµ-Myc mouse model of B-cell lymphoma, we discovered that both NF-κB subunit mutations resulted in CHK1 inhibitor resistance, arising from either loss or alteration of CHK1 activity, respectively. However, since Eµ-Myc lymphomas depend on CHK1 activity to cope with high levels of DNA replication stress and consequent genomic instability, it was not clear how these mutant NF-κB subunit lymphomas were able to survive. To understand these survival mechanisms and to identify potential compensatory bypass signalling pathways in these lymphomas, we applied a multi-omics strategy. With c-Rel-/- Eµ-Myc lymphomas we observed high levels of Phosphatidyl-inositol 3-kinase (PI3K) and AKT pathway activation. Moreover, treatment with the PI3K inhibitor Pictilisib (GDC-0941) selectively inhibited the growth of reimplanted c-Rel-/- and RelAT505A, but not wild type (WT) Eµ-Myc lymphomas. We also observed up-regulation of a RHO/RAC pathway gene expression signature in both Eµ-Myc NF-κB subunit mutation models. Further investigation demonstrated activation of the RHO/RAC effector p21-activated kinase (PAK) 2. Here, the PAK inhibitor, PF-3758309 successfully overcame resistance of RelAT505A but not WT lymphomas. These findings demonstrate that up-regulation of multiple bypass pathways occurs in CHK1 inhibitor resistant Eµ-Myc lymphomas. Consequently, drugs targeting these pathways could potentially be used as either second line or combinatorial therapies to aid the successful clinical application of CHK1 inhibitors.

Strain-specific impacts of probiotics are a significant driver of gut microbiome development in very preterm infants.

The development of the gut microbiome from birth plays important roles in short- and long-term health, but factors influencing preterm gut microbiome development are poorly understood. In the present study, we use metagenomic sequencing to analyse 1,431 longitudinal stool samples from 123 very preterm infants (<32 weeks' gestation) who did not develop intestinal disease or sepsis over a study period of 10 years. During the study period, one cohort had no probiotic exposure whereas two cohorts were given different probiotic products: Infloran (Bifidobacterium bifidum and Lactobacillus acidophilus) or Labinic (B. bifidum, B. longum subsp. infantis and L. acidophilus). Mothers' own milk, breast milk fortifier, antibiotics and probiotics were significantly associated with the gut microbiome, with probiotics being the most significant factor. Probiotics drove microbiome transition into different preterm gut community types (PGCTs), each enriched in a different Bifidobacterium sp. and significantly associated with increased postnatal age. Functional analyses identified stool metabolites associated with PGCTs and, in preterm-derived organoids, sterile faecal supernatants impacted intestinal, organoid monolayer, gene expression in a PGCT-specific manner. The present study identifies specific influencers of gut microbiome development in very preterm infants, some of which overlap with those impacting term infants. The results highlight the importance of strain-specific differences in probiotic products and their impact on host interactions in the preterm gut.

hiPSC-derived bone marrow milieu identifies a clinically actionable driver of niche-mediated treatment resistance in leukemia.

Leukemia cells re-program their microenvironment to augment blast proliferation and enhance treatment resistance. Means of clinically targeting such niche-driven treatment resistance remain ambiguous. We develop human induced pluripotent stem cell (hiPSC)-engineered niches to reveal druggable cancer-niche dependencies. We reveal that mesenchymal (iMSC) and vascular niche-like (iANG) hiPSC-derived cells support ex vivo proliferation of patient-derived leukemia cells, affect dormancy, and mediate treatment resistance. iMSCs protect dormant and cycling blasts against dexamethasone, while iANGs protect only dormant blasts. Leukemia proliferation and protection from dexamethasone-induced apoptosis is dependent on cancer-niche interactions mediated by CDH2. Consequently, we test CDH2 antagonist ADH-1 (previously in Phase I/II trials for solid tumors) in a very aggressive patient-derived xenograft leukemia mouse model. ADH-1 shows high in vivo efficacy; ADH-1/dexamethasone combination is superior to dexamethasone alone, with no ADH-1-conferred additional toxicity. These findings provide a proof-of-concept starting point to develop improved, potentially safer therapeutics targeting niche-mediated cancer dependencies in blood cancers.

Single-cell DNA sequencing identifies risk-associated clonal complexity and evolutionary trajectories in childhood medulloblastoma development.

We reconstructed the natural history and temporal evolution of the most common childhood brain malignancy, medulloblastoma, by single-cell whole-genome sequencing (sc-WGS) of tumours representing its major molecular sub-classes and clinical risk groups. Favourable-risk disease sub-types assessed (MBWNT and infant desmoplastic/nodular MBSHH) typically comprised a single clone with no evidence of further evolution. In contrast, highest risk sub-classes (MYC-amplified MBGroup3 and TP53-mutated MBSHH) were most clonally diverse and displayed gradual evolutionary trajectories. Clinically adopted biomarkers (e.g. chromosome 6/17 aberrations; CTNNB1/TP53 mutations) were typically early-clonal/initiating events, exploitable as targets for early-disease detection; in analyses of spatially distinct tumour regions, a single biopsy was sufficient to assess their status. Importantly, sc-WGS revealed novel events which arise later and/or sub-clonally and more commonly display spatial diversity; their clinical significance and role in disease evolution post-diagnosis now require establishment. These findings reveal diverse modes of tumour initiation and evolution in the major medulloblastoma sub-classes, with pathogenic relevance and clinical potential.

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

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

Heteroplasmic mitochondrial DNA variants in cardiovascular diseases.

Mitochondria are implicated in the pathogenesis of cardiovascular diseases (CVDs) but the reasons for this are not well understood. Maternally-inherited population variants of mitochondrial DNA (mtDNA) which affect all mtDNA molecules (homoplasmic) are associated with cardiometabolic traits and the risk of developing cardiovascular disease. However, it is not known whether mtDNA mutations only affecting a proportion of mtDNA molecules (heteroplasmic) also play a role. To address this question, we performed a high-depth (~1000-fold) mtDNA sequencing of blood DNA in 1,399 individuals with hypertension (HTN), 1,946 with ischemic heart disease (IHD), 2,146 with ischemic stroke (IS), and 723 healthy controls. We show that the per individual burden of heteroplasmic single nucleotide variants (mtSNVs) increases with age. The age-effect was stronger for low-level heteroplasmies (heteroplasmic fraction, HF, 5-10%), likely reflecting acquired somatic events based on trinucleotide mutational signatures. After correcting for age and other confounders, intermediate heteroplasmies (HF 10-95%) were more common in hypertension, particularly involving non-synonymous variants altering the amino acid sequence of essential respiratory chain proteins. These findings raise the possibility that heteroplasmic mtSNVs play a role in the pathophysiology of hypertension.

Blood mRNA Expression in Alzheimer's Disease and Dementia With Lewy Bodies.

OBJECTIVES: The objective of this study was to investigate the expression of genes in Alzheimer's disease (AD) and dementia with Lewy bodies (DLB), both at the mild cognitive impairment (MCI) and dementia stages, to improve our understanding of disease pathophysiology and investigate the potential for diagnostic and prognostic biomarkers based on mRNA expression. DESIGN: Cross-sectional observational study. SETTING: University research center. PARTICIPANTS: People with MCI with Lewy bodies (MCI-LB, n=55), MCI-AD (n=19), DLB (n=38), AD (n=24) and a cognitively unimpaired comparison group (n=28). MEASUREMENTS: Ribonucleic acid sequencing of whole blood. Differentially expressed genes (DEGs) were identified and gene set enrichment analysis was carried out. RESULTS: Compared with the cognitively unimpaired group, there were 22 DEGs in MCI-LB/DLB and 61 DEGs in MCI-AD/AD. DEGS were also identified when comparing the two disease groups. Expression of ANP32A was associated with more rapid cognitive decline in MCI-AD/AD. Gene set enrichment analysis identified downregulation in gene sets including MYC targets and oxidative phosphorylation in MCI-LB/DLB; upregulation of immune and inflammatory responses in MCI-AD/AD; and upregulation of interferon-α and -γ responses in MCI-AD/AD compared with MCI-LB/DLB. CONCLUSION: This study identified multiple DEGs in MCI-LB/DLB and MCI-AD/AD. One of these DEGs, ANP32A, may be a prognostic marker in AD. Genes related to mitochondrial function were downregulated in MCI-LB/DLB. Previously reported upregulation of genes associated with inflammation and immune responses in MCI-AD/AD was confirmed in this cohort. Differences in interferon responses between MCI-AD/AD and MCI-LB/DLB suggest that there are key differences in peripheral immune responses between these diseases.

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

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