Single-cell analysis reveals the continuum of human lympho-myeloid progenitor cells.

The hierarchy of human hemopoietic progenitor cells that produce lymphoid and granulocytic-monocytic (myeloid) lineages is unclear. Multiple progenitor populations produce lymphoid and myeloid cells, but they remain incompletely characterized. Here we demonstrated that lympho-myeloid progenitor populations in cord blood - lymphoid-primed multi-potential progenitors (LMPPs), granulocyte-macrophage progenitors (GMPs) and multi-lymphoid progenitors (MLPs) - were functionally and transcriptionally distinct and heterogeneous at the clonal level, with progenitors of many different functional potentials present. Although most progenitors had the potential to develop into only one mature cell type ('uni-lineage potential'), bi- and rarer multi-lineage progenitors were present among LMPPs, GMPs and MLPs. Those findings, coupled with single-cell expression analyses, suggest that a continuum of progenitors execute lymphoid and myeloid differentiation, rather than only uni-lineage progenitors' being present downstream of stem cells.

A polymorphic p53 response element in KIT ligand influences cancer risk and has undergone natural selection.

The ability of p53 to regulate transcription is crucial for tumor suppression and implies that inherited polymorphisms in functional p53-binding sites could influence cancer. Here, we identify a polymorphic p53 responsive element and demonstrate its influence on cancer risk using genome-wide data sets of cancer susceptibility loci, genetic variation, p53 occupancy, and p53-binding sites. We uncover a single-nucleotide polymorphism (SNP) in a functional p53-binding site and establish its influence on the ability of p53 to bind to and regulate transcription of the KITLG gene. The SNP resides in KITLG and associates with one of the largest risks identified among cancer genome-wide association studies. We establish that the SNP has undergone positive selection throughout evolution, signifying a selective benefit, but go on to show that similar SNPs are rare in the genome due to negative selection, indicating that polymorphisms in p53-binding sites are primarily detrimental to humans.

Canonical Notch signaling is dispensable for adult steady-state and stress myelo-erythropoiesis.

Although an essential role for canonical Notch signaling in generation of hematopoietic stem cells in the embryo and in thymic T-cell development is well established, its role in adult bone marrow (BM) myelopoiesis remains unclear. Some studies, analyzing myeloid progenitors in adult mice with inhibited Notch signaling, implicated distinct roles of canonical Notch signaling in regulation of progenitors for the megakaryocyte, erythroid, and granulocyte-macrophage cell lineages. However, these studies might also have targeted other pathways. Therefore, we specifically deleted, in adult BM, the transcription factor recombination signal-binding protein J κ (Rbpj), through which canonical signaling from all Notch receptors converges. Notably, detailed progenitor staging established that canonical Notch signaling is fully dispensable for all investigated stages of megakaryocyte, erythroid, and myeloid progenitors in steady state unperturbed hematopoiesis, after competitive BM transplantation, and in stress-induced erythropoiesis. Moreover, expression of key regulators of these hematopoietic lineages and Notch target genes were unaffected by Rbpj deficiency in BM progenitor cells.

Impact of spliceosome mutations on RNA splicing in myelodysplasia: dysregulated genes/pathways and clinical associations.

SF3B1, SRSF2, and U2AF1 are the most frequently mutated splicing factor genes in the myelodysplastic syndromes (MDS). We have performed a comprehensive and systematic analysis to determine the effect of these commonly mutated splicing factors on pre-mRNA splicing in the bone marrow stem/progenitor cells and in the erythroid and myeloid precursors in splicing factor mutant MDS. Using RNA-seq, we determined the aberrantly spliced genes and dysregulated pathways in CD34+ cells of 84 patients with MDS. Splicing factor mutations result in different alterations in splicing and largely affect different genes, but these converge in common dysregulated pathways and cellular processes, focused on RNA splicing, protein synthesis, and mitochondrial dysfunction, suggesting common mechanisms of action in MDS. Many of these dysregulated pathways and cellular processes can be linked to the known disease pathophysiology associated with splicing factor mutations in MDS, whereas several others have not been previously associated with MDS, such as sirtuin signaling. We identified aberrantly spliced events associated with clinical variables, and isoforms that independently predict survival in MDS and implicate dysregulation of focal adhesion and extracellular exosomes as drivers of poor survival. Aberrantly spliced genes and dysregulated pathways were identified in the MDS-affected lineages in splicing factor mutant MDS. Functional studies demonstrated that knockdown of the mitosis regulators SEPT2 and AKAP8, aberrantly spliced target genes of SF3B1 and SRSF2 mutations, respectively, led to impaired erythroid cell growth and differentiation. This study illuminates the effect of the common spliceosome mutations on the MDS phenotype and provides novel insights into disease pathophysiology.

Cell-intrinsic depletion of Aml1-ETO-expressing pre-leukemic hematopoietic stem cells by K-Ras activating mutation.

Somatic mutations in acute myeloid leukemia are acquired sequentially and hierarchically. First, pre-leukemic mutations, such as t(8;21) that encodes AML1-ETO, are acquired within the hematopoietic stem cell (HSC) compartment, while signaling pathway mutations, including KRAS activating mutations, are late events acquired during transformation of leukemic progenitor cells and are rarely detectable in HSC. This raises the possibility that signaling pathway mutations are detrimental to clonal expansion of pre-leukemic HSC. To address this hypothesis, we used conditional genetics to introduce Aml1-ETO and K-RasG12D into murine HSC, either individually or in combination. In the absence of activated Ras, Aml1-ETO-expressing HSC conferred a competitive advantage. However, activated K-Ras had a marked detrimental effect on Aml1-ETO-expressing HSC, leading to loss of both phenotypic and functional HSC. Cell cycle analysis revealed a loss of quiescence in HSC co-expressing Aml1-ETO and K-RasG12D, accompanied by an enrichment in E2F and Myc target gene expression and depletion of HSC self-renewal-associated gene expression. These findings provide a mechanistic basis for the observed absence of KRAS signaling mutations in the pre-malignant HSC compartment.

CTAS: a CT score to quantify disease activity in pulmonary sarcoidosis.

BACKGROUND: A major gap in the management of sarcoidosis is the lack of accessible and objective methods to measure disease activity. Since 90% of patients have pulmonary involvement, we explored if a disease activity score based on thoracic CT scans could address this clinical issue. METHODS: High-resolution CT scans from 100 consecutive patients with sarcoidosis at a regional sarcoidosis service were scored for extent of CT abnormalities known to relate to granuloma or lymphocytic infiltration from published CT-pathological studies. These individual abnormality scores were then correlated against serum ACE, sIL-2R and change in FVC to identify CT abnormalities that reflect contemporaneous disease activity. The sum of these scores, or CT Activity Score (CTAS), was then validated against FVC response to treatment. FINDINGS: CT extent scores for nodularity, ground-glass opacification, interlobular septal thickening and consolidation correlated significantly with at least one of the disease activity parameters and were used to form CTAS. CTAS was found to predict FVC response to treatment at 1 year and was highly reproducible between radiologists. An abbreviated CTAS (aCTAS), constructed from presence or absence of the four CT abnormalities, also showed significant correlation with FVC response to treatment. CTAS and aCTAS also correlated with response to treatment in the fibrotic subgroup. INTERPRETATION: CTAS provides a concept for an objective and reproducible CT scoring method to quantify disease activity in sarcoidosis. The score can potentially be used to stratify patients according to disease activity, determine response to treatment and establish if fibrotic sarcoidosis is active.

Single cell spatial analysis reveals inflammatory foci of immature neutrophil and CD8 T cells in COVID-19 lungs.

Single cell spatial interrogation of the immune-structural interactions in COVID -19 lungs is challenging, mainly because of the marked cellular infiltrate and architecturally distorted microstructure. To address this, we develop a suite of mathematical tools to search for statistically significant co-locations amongst immune and structural cells identified using 37-plex imaging mass cytometry. This unbiased method reveals a cellular map interleaved with an inflammatory network of immature neutrophils, cytotoxic CD8 T cells, megakaryocytes and monocytes co-located with regenerating alveolar progenitors and endothelium. Of note, a highly active cluster of immature neutrophils and CD8 T cells, is found spatially linked with alveolar progenitor cells, and temporally with the diffuse alveolar damage stage. These findings offer further insights into how immune cells interact in the lungs of severe COVID-19 disease. We provide our pipeline [Spatial Omics Oxford Pipeline (SpOOx)] and visual-analytical tool, Multi-Dimensional Viewer (MDV) software, as a resource for spatial analysis.

Effect of five genetic variants associated with lung function on the risk of chronic obstructive lung disease, and their joint effects on lung function.

RATIONALE: Genomic loci are associated with FEV1 or the ratio of FEV1 to FVC in population samples, but their association with chronic obstructive pulmonary disease (COPD) has not yet been proven, nor have their combined effects on lung function and COPD been studied. OBJECTIVES: To test association with COPD of variants at five loci (TNS1, GSTCD, HTR4, AGER, and THSD4) and to evaluate joint effects on lung function and COPD of these single-nucleotide polymorphisms (SNPs), and variants at the previously reported locus near HHIP. METHODS: By sampling from 12 population-based studies (n = 31,422), we obtained genotype data on 3,284 COPD case subjects and 17,538 control subjects for sentinel SNPs in TNS1, GSTCD, HTR4, AGER, and THSD4. In 24,648 individuals (including 2,890 COPD case subjects and 13,862 control subjects), we additionally obtained genotypes for rs12504628 near HHIP. Each allele associated with lung function decline at these six SNPs contributed to a risk score. We studied the association of the risk score to lung function and COPD. MEASUREMENTS AND MAIN RESULTS: Association with COPD was significant for three loci (TNS1, GSTCD, and HTR4) and the previously reported HHIP locus, and suggestive and directionally consistent for AGER and TSHD4. Compared with the baseline group (7 risk alleles), carrying 10-12 risk alleles was associated with a reduction in FEV1 (ß = -72.21 ml, P = 3.90 × 10(-4)) and FEV1/FVC (ß = -1.53%, P = 6.35 × 10(-6)), and with COPD (odds ratio = 1.63, P = 1.46 × 10(-5)). CONCLUSIONS: Variants in TNS1, GSTCD, and HTR4 are associated with COPD. Our highest risk score category was associated with a 1.6-fold higher COPD risk than the population average score.

No evidence for age-related differences in mitochondrial RNA quality in the female germline.

Abstract: Mitochondrial quality is implicated as a contributor to declining fertility with aging. We investigated mitochondrial transcripts in oocytes and their associated cumulus cells from mice of different ages using RNA-seq. Mice aged 3 weeks, 9 weeks, and 1 year were superovulated, and 48 h later, oocyte cumulus complexes were collected by follicle puncture. We did not detect any major differences that could be attributed to aging. However, mitochondrial RNA transcripts which deviated from the consensus sequence were found at a higher frequency in cumulus cells than in their corresponding oocyte. Previous investigations have shown that variation in the sequence of mtRNA transcripts is substantial, and at least some of this can be accounted for by post-transcriptional modifications which impact base calling during sequencing. Our data would be consistent with either less post-transcriptional modification in mitochondrial RNA from oocytes than cumulus cells or with lower mtDNA mutational load. Lay summary: Women become less fertile as they age. Shortage of energy contributes to this, caused by a decline in the quality of mitochondria (the powerhouses of the cell) in the egg. Genes are the blueprint for the cell. They are made of DNA which is copied into an RNA message, or instructions, for making proteins. We counted differences in the RNA message of developing eggs and the cells that support them during development (cumulus cells). We compared the number of these differences in mice of different ages. These age groups represent mice had not reached puberty, those of prime reproductive age, and old mothers. We did not find any differences linked to the age of the mice. However, we did find differences between the egg and the cumulus cells. In most cases, there were lower levels of mutations in eggs than there were in cumulus cells.

The chromatin remodeller ATRX facilitates diverse nuclear processes, in a stochastic manner, in both heterochromatin and euchromatin.

The chromatin remodeller ATRX interacts with the histone chaperone DAXX to deposit the histone variant H3.3 at sites of nucleosome turnover. ATRX is known to bind repetitive, heterochromatic regions of the genome including telomeres, ribosomal DNA and pericentric repeats, many of which are putative G-quadruplex forming sequences (PQS). At these sites ATRX plays an ancillary role in a wide range of nuclear processes facilitating replication, chromatin modification and transcription. Here, using an improved protocol for chromatin immunoprecipitation, we show that ATRX also binds active regulatory elements in euchromatin. Mutations in ATRX lead to perturbation of gene expression associated with a reduction in chromatin accessibility, histone modification, transcription factor binding and deposition of H3.3 at the sequences to which it normally binds. In erythroid cells where downregulation of α-globin expression is a hallmark of ATR-X syndrome, perturbation of chromatin accessibility and gene expression occurs in only a subset of cells. The stochastic nature of this process suggests that ATRX acts as a general facilitator of cell specific transcriptional and epigenetic programmes, both in heterochromatin and euchromatin.

BET inhibition disrupts transcription but retains enhancer-promoter contact.

Enhancers are DNA sequences that enable complex temporal and tissue-specific regulation of genes in higher eukaryotes. Although it is not entirely clear how enhancer-promoter interactions can increase gene expression, this proximity has been observed in multiple systems at multiple loci and is thought to be essential for the maintenance of gene expression. Bromodomain and Extra-Terminal domain (BET) and Mediator proteins have been shown capable of forming phase condensates and are thought to be essential for super-enhancer function. Here, we show that targeting of cells with inhibitors of BET proteins or pharmacological degradation of BET protein Bromodomain-containing protein 4 (BRD4) has a strong impact on transcription but very little impact on enhancer-promoter interactions. Dissolving phase condensates reduces BRD4 and Mediator binding at enhancers and can also strongly affect gene transcription, without disrupting enhancer-promoter interactions. These results suggest that activation of transcription and maintenance of enhancer-promoter interactions are separable events. Our findings further indicate that enhancer-promoter interactions are not dependent on high levels of BRD4 and Mediator, and are likely maintained by a complex set of factors including additional activator complexes and, at some sites, CTCF and cohesin.

IL-6 effector function of group 2 innate lymphoid cells (ILC2) is NOD2 dependent.

Cutaneous group 2 innate lymphoid cells (ILC2) are spatially and epigenetically poised to respond to barrier compromise and associated immunological threats. ILC2, lacking rearranged antigen-specific receptors, are primarily activated by damage-associated cytokines and respond with type 2 cytokine production. To investigate ILC2 potential for direct sensing of skin pathogens and allergens, we performed RNA sequencing of ILC2 derived from in vivo challenged human skin or blood. We detected expression of NOD2 and TLR2 by skin and blood ILC2. Stimulation of ILC2 with TLR2 agonist alone not only induced interleukin-5 (IL-5) and IL-13 expression but also elicited IL-6 expression in combination with Staphylococcus aureus muramyl dipeptide (MDP). Heat-killed skin-resident bacteria provoked an IL-6 profile in ILC2 in vitro that was notably impaired in ILC2 derived from patients with nucleotide-binding oligomerization domain-containing protein 2 (NOD2) mutations. In addition, we show that NOD2 signaling can stimulate autophagy in ILC2, which was also impaired in patients with NOD2 mutations. Here, we have identified a role for ILC2 NOD2 signaling in the differential regulation of ILC2-derived IL-6 and have reported a previously unrecognized pathway of direct ILC2 bacterial sensing.

Multi-Modal Characterization of Monocytes in Idiopathic Pulmonary Fibrosis Reveals a Primed Type I Interferon Immune Phenotype.

Idiopathic pulmonary fibrosis (IPF) is the most severe form of chronic lung fibrosis. Circulating monocytes have been implicated in immune pathology in IPF but their phenotype is unknown. In this work, we determined the immune phenotype of monocytes in IPF using multi-colour flow cytometry, RNA sequencing and corresponding serum factors, and mapped the main findings to amount of lung fibrosis and single cell transcriptomic landscape of myeloid cells in IPF lungs. We show that monocytes from IPF patients displayed increased expression of CD64 (FcγR1) which correlated with amount of lung fibrosis, and an amplified type I IFN response ex vivo. These were accompanied by markedly raised CSF-1 levels, IL-6, and CCL-2 in serum of IPF patients. Interrogation of single cell transcriptomic data from human IPF lungs revealed increased proportion of CD64hi monocytes and "transitional macrophages" with higher expression of CCL-2 and type I IFN genes. Our study shows that monocytes in IPF patients are phenotypically distinct from age-matched controls, with a primed type I IFN pathway that may contribute to driving chronic inflammation and fibrosis. These findings strengthen the potential role of monocytes in the pathogenesis of IPF.

Single-cell atlas of colonic CD8+ T cells in ulcerative colitis.

Colonic antigen-experienced lymphocytes such as tissue-resident memory CD8+ T cells can respond rapidly to repeated antigen exposure. However, their cellular phenotypes and the mechanisms by which they drive immune regulation and inflammation remain unclear. Here we compiled an unbiased atlas of human colonic CD8+ T cells in health and ulcerative colitis (UC) using single-cell transcriptomics with T-cell receptor repertoire analysis and mass cytometry. We reveal extensive heterogeneity in CD8+ T-cell composition, including expanded effector and post-effector terminally differentiated CD8+ T cells. While UC-associated CD8+ effector T cells can trigger tissue destruction and produce tumor necrosis factor (TNF)-α, post-effector cells acquire innate signatures to adopt regulatory functions that may mitigate excessive inflammation. Thus, we identify colonic CD8+ T-cell phenotypes in health and UC, define their clonal relationships and characterize terminally differentiated dysfunctional UC CD8+ T cells expressing IL-26, which attenuate acute colitis in a humanized IL-26 transgenic mouse model.

Reconstruction of the Global Neural Crest Gene Regulatory Network In Vivo.

Precise control of developmental processes is encoded in the genome in the form of gene regulatory networks (GRNs). Such multi-factorial systems are difficult to decode in vertebrates owing to their complex gene hierarchies and dynamic molecular interactions. Here we present a genome-wide in vivo reconstruction of the GRN underlying development of the multipotent neural crest (NC) embryonic cell population. By coupling NC-specific epigenomic and transcriptional profiling at population and single-cell levels with genome/epigenome engineering in vivo, we identify multiple regulatory layers governing NC ontogeny, including NC-specific enhancers and super-enhancers, novel trans-factors, and cis-signatures allowing reverse engineering of the NC-GRN at unprecedented resolution. Furthermore, identification and dissection of divergent upstream combinatorial regulatory codes has afforded new insights into opposing gene circuits that define canonical and neural NC fates early during NC ontogeny. Our integrated approach, allowing dissection of cell-type-specific regulatory circuits in vivo, has broad implications for GRN discovery and investigation.

DOT1L inhibition reveals a distinct subset of enhancers dependent on H3K79 methylation.

Enhancer elements are a key regulatory feature of many important genes. Several general features including the presence of specific histone modifications are used to demarcate potentially active enhancers. Here we reveal that putative enhancers marked with H3 lysine 79 (H3K79) di or trimethylation (me2/3) (which we name H3K79me2/3 enhancer elements or KEEs) can be found in multiple cell types. Mixed lineage leukemia gene (MLL) rearrangements (MLL-r) such as MLL-AF4 are a major cause of incurable acute lymphoblastic leukemias (ALL). Using the DOT1L inhibitor EPZ-5676 in MLL-AF4 leukemia cells, we show that H3K79me2/3 is required for maintaining chromatin accessibility, histone acetylation and transcription factor binding specifically at KEEs but not non-KEE enhancers. We go on to show that H3K79me2/3 is essential for maintaining enhancer-promoter interactions at a subset of KEEs. Together, these data implicate H3K79me2/3 as having a functional role at a subset of active enhancers in MLL-AF4 leukemia cells.

Transcriptomic profiling of the myeloma bone-lining niche reveals BMP signalling inhibition to improve bone disease.

Multiple myeloma is an incurable, bone marrow-dwelling malignancy that disrupts bone homeostasis causing skeletal damage and pain. Mechanisms underlying myeloma-induced bone destruction are poorly understood and current therapies do not restore lost bone mass. Using transcriptomic profiling of isolated bone lining cell subtypes from a murine myeloma model, we find that bone morphogenetic protein (BMP) signalling is upregulated in stromal progenitor cells. BMP signalling has not previously been reported to be dysregulated in myeloma bone disease. Inhibition of BMP signalling in vivo using either a small molecule BMP receptor antagonist or a solubilized BMPR1a-FC receptor ligand trap prevents trabecular and cortical bone volume loss caused by myeloma, without increasing tumour burden. BMP inhibition directly reduces osteoclastogenesis, increases osteoblasts and bone formation, and suppresses bone marrow sclerostin levels. In summary we describe a novel role for the BMP pathway in myeloma-induced bone disease that can be therapeutically targeted.

Does osteogenic potential of clonal human bone marrow mesenchymal stem/stromal cells correlate with their vascular supportive ability?

BACKGROUND: Human bone marrow-derived mesenchymal stem/stromal cells (hBM MSCs) have multiple functions, critical for skeletal formation and function. Their functional heterogeneity, however, represents a major challenge for their isolation and in developing potency and release assays to predict their functionality prior to transplantation. Additionally, potency, biomarker profiles and defining mechanisms of action in a particular clinical setting are increasing requirements of Regulatory Agencies for release of hBM MSCs as Advanced Therapy Medicinal Products for cellular therapies. Since the healing of bone fractures depends on the coupling of new blood vessel formation with osteogenesis, we hypothesised that a correlation between the osteogenic and vascular supportive potential of individual hBM MSC-derived CFU-F (colony forming unit-fibroblastoid) clones might exist. METHODS: We tested this by assessing the lineage (i.e. adipogenic (A), osteogenic (O) and/or chondrogenic (C)) potential of individual hBM MSC-derived CFU-F clones and determining if their osteogenic (O) potential correlated with their vascular supportive profile in vitro using lineage differentiation assays, endothelial-hBM MSC vascular co-culture assays and transcriptomic (RNAseq) analyses. RESULTS: Our results demonstrate that the majority of CFU-F (95%) possessed tri-lineage, bi-lineage or uni-lineage osteogenic capacity, with 64% of the CFU-F exhibiting tri-lineage AOC potential. We found a correlation between the osteogenic and vascular tubule supportive activity of CFU-F clones, with the strength of this association being donor dependent. RNAseq of individual clones defined gene fingerprints relevant to this correlation. CONCLUSIONS: This study identified a donor-dependent correlation between osteogenic and vascular supportive potential of hBM MSCs and important gene signatures that support these functions that are relevant to their bone regenerative properties.

From Pioneer to Repressor: Bimodal foxd3 Activity Dynamically Remodels Neural Crest Regulatory Landscape In Vivo.

The neural crest (NC) is a transient embryonic stem cell-like population characterized by its multipotency and broad developmental potential. Here, we perform NC-specific transcriptional and epigenomic profiling of foxd3-mutant cells in vivo to define the gene regulatory circuits controlling NC specification. Together with global binding analysis obtained by foxd3 biotin-ChIP and single cell profiles of foxd3-expressing premigratory NC, our analysis shows that, during early steps of NC formation, foxd3 acts globally as a pioneer factor to prime the onset of genes regulating NC specification and migration by re-arranging the chromatin landscape, opening cis-regulatory elements and reshuffling nucleosomes. Strikingly, foxd3 then gradually switches from an activator to its well-described role as a transcriptional repressor and potentially uses differential partners for each role. Taken together, these results demonstrate that foxd3 acts bimodally in the neural crest as a switch from "permissive" to "repressive" nucleosome and chromatin organization to maintain multipotency and define cell fates.

SCL/TAL1 cooperates with Polycomb RYBP-PRC1 to suppress alternative lineages in blood-fated cells.

During development, it is unclear if lineage-fated cells derive from multilineage-primed progenitors and whether active mechanisms operate to restrict cell fate. Here we investigate how mesoderm specifies into blood-fated cells. We document temporally restricted co-expression of blood (Scl/Tal1), cardiac (Mesp1) and paraxial (Tbx6) lineage-affiliated transcription factors in single cells, at the onset of blood specification, supporting the existence of common progenitors. At the same time-restricted stage, absence of SCL results in expansion of cardiac/paraxial cell populations and increased cardiac/paraxial gene expression, suggesting active suppression of alternative fates. Indeed, SCL normally activates expression of co-repressor ETO2 and Polycomb-PRC1 subunits (RYBP, PCGF5) and maintains levels of Polycomb-associated histone marks (H2AK119ub/H3K27me3). Genome-wide analyses reveal ETO2 and RYBP co-occupy most SCL target genes, including cardiac/paraxial loci. Reduction of Eto2 or Rybp expression mimics Scl-null cardiac phenotype. Therefore, SCL-mediated transcriptional repression prevents mis-specification of blood-fated cells, establishing active repression as central to fate determination processes.