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

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

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.

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

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

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

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

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.

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.

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.

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.

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.

Msx1 haploinsufficiency modifies the Pax9-deficient cardiovascular phenotype.

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

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

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

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

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

Extreme-Depth Re-sequencing of Mitochondrial DNA Finds No Evidence of Paternal Transmission in Humans.

Recent reports have questioned the accepted dogma that mammalian mitochondrial DNA (mtDNA) is strictly maternally inherited. In humans, the argument hinges on detecting a signature of inter-molecular recombination in mtDNA sequences sampled at the population level, inferring a paternal source for the mixed haplotypes. However, interpreting these data is fraught with difficulty, and direct experimental evidence is lacking. Using extreme-high depth mtDNA re-sequencing up to ~1.2 million-fold coverage, we find no evidence that paternal mtDNA haplotypes are transmitted to offspring in humans, thus excluding a simple dilution mechanism for uniparental transmission of mtDNA present in all healthy individuals. Our findings indicate that an active mechanism eliminates paternal mtDNA which likely acts at the molecular level.

Universal heteroplasmy of human mitochondrial DNA.

Mammalian cells contain thousands of copies of mitochondrial DNA (mtDNA). At birth, these are thought to be identical in most humans. Here, we use long read length ultra-deep resequencing-by-synthesis to interrogate regions of the mtDNA genome from related and unrelated individuals at unprecedented resolution. We show that very low-level heteroplasmic variance is present in all tested healthy individuals, and is likely to be due to both inherited and somatic single base substitutions. Using this approach, we demonstrate an increase in mtDNA mutations in the skeletal muscle of patients with a proofreading-deficient mtDNA polymerase γ due to POLG mutations. In contrast, we show that OPA1 mutations, which indirectly affect mtDNA maintenance, do not increase point mutation load. The demonstration of universal mtDNA heteroplasmy has fundamental implications for our understanding of mtDNA inheritance and evolution. Ostensibly de novo somatic mtDNA mutations, seen in mtDNA maintenance disorders and neurodegenerative disease and aging, will partly be due to the clonal expansion of low-level inherited variants.

Molecular pathogenesis of polymerase γ-related neurodegeneration.

OBJECTIVE: Polymerase gamma (POLG) mutations are a common cause of mitochondrial disease and have also been linked to neurodegeneration and aging. We studied the molecular mechanisms underlying POLG-related neurodegeneration using postmortem tissue from a large number of patients. METHODS: Clinical information was available from all subjects. Formalin-fixed and frozen brain tissue from 15 patients and 23 controls was studied employing a combination of histopathology, immunohistochemistry, and molecular studies of microdissected neurons. RESULTS: The primary consequence of POLG mutation in neurons is mitochondrial DNA depletion. This was already present in infants with little evidence of neuronal loss or mitochondrial dysfunction. With longer disease duration, we found an additional, progressive accumulation of mitochondrial DNA deletions and point mutations accompanied by increasing numbers of complex I-deficient neurons. Progressive neurodegeneration primarily affected the cerebellar systems and dopaminergic cells of the substantia nigra. Superimposed on this chronic process were acute, focal cortical lesions that correlated with epileptogenic foci and that showed massive neuronal loss. INTERPRETATION: POLG mutations appear to compromise neuronal respiration via a combination of early and stable depletion and a progressive somatic mutagenesis of the mitochondrial genome. This leads to 2 distinct but overlapping biological processes: a chronic neurodegeneration reflected clinically by progressive ataxia and cognitive impairment, and an acute focal neuronal necrosis that appears to be related to the presence of epileptic seizures. Our findings offer an explanation of the acute-on-chronic clinical course of this common mitochondrial encephalopathy.

Post-transcriptional exon shuffling events in humans can be evolutionarily conserved and abundant.

In silico analyses have established that transcripts from some genes can be processed into RNAs with rearranged exon order relative to genomic structure (post-transcriptional exon shuffling, or PTES). Although known to contribute to transcriptome diversity in some species, to date the structure, distribution, abundance, and functional significance of human PTES transcripts remains largely unknown. Here, using high-throughput transcriptome sequencing, we identify 205 putative human PTES products from 176 genes. We validate 72 out of 112 products analyzed using RT-PCR, and identify additional PTES products structurally related to 61% of validated targets. Sequencing of these additional products reveals GT-AG dinucleotides at >95% of the splice junctions, confirming that they are processed by the spliceosome. We show that most PTES transcripts are expressed in a wide variety of human tissues, that they can be polyadenylated, and that some are conserved in mouse. We also show that they can extend into 5' and 3' UTRs, consistent with formation via trans-splicing of independent pre-mRNA molecules. Finally, we use real-time PCR to compare the abundance of PTES exon junctions relative to canonical exon junctions within the transcripts from seven genes. PTES exon junctions are present at <0.01% to >90% of the levels of canonical junctions, with transcripts from MAN1A2, PHC3, TLE4, and CDK13 exhibiting the highest levels. This is the first systematic experimental analysis of PTES in human, and it suggests both that the phenomenon is much more widespread than previously thought and that some PTES transcripts could be functional.