Robustness of single-cell RNA-seq for identifying differentially expressed genes.

BACKGROUND: A common feature of single-cell RNA-seq (scRNA-seq) data is that the number of cells in a cell cluster may vary widely, ranging from a few dozen to several thousand. It is not clear whether scRNA-seq data from a small number of cells allow robust identification of differentially expressed genes (DEGs) with various characteristics. RESULTS: We addressed this question by performing scRNA-seq and poly(A)-dependent bulk RNA-seq in comparable aliquots of human induced pluripotent stem cells-derived, purified vascular endothelial and smooth muscle cells. We found that scRNA-seq data needed to have 2,000 or more cells in a cluster to identify the majority of DEGs that would show modest differences in a bulk RNA-seq analysis. On the other hand, clusters with as few as 50-100 cells may be sufficient for identifying the majority of DEGs that would have extremely small p values or transcript abundance greater than a few hundred transcripts per million in a bulk RNA-seq analysis. CONCLUSION: Findings of the current study provide a quantitative reference for designing studies that aim for identifying DEGs for specific cell clusters using scRNA-seq data and for interpreting results of such studies.

The miR-23a∼27a∼24-2 microRNA Cluster Promotes Inflammatory Polarization of Macrophages.

Macrophages are critical for regulating inflammatory responses. Environmental signals polarize macrophages to either a proinflammatory (M1) state or an anti-inflammatory (M2) state. We observed that the microRNA (miRNA) cluster mirn23a, coding for miRs-23a, -27a, and -24-2, regulates mouse macrophage polarization. Gene expression analysis of mirn23a-deficient myeloid progenitors revealed a decrease in TLR and IFN signaling. Mirn23a -/- bone marrow-derived macrophages (BMDMs) have an attenuated response to LPS, demonstrating an anti-inflammatory phenotype in mature cells. In vitro, mirn23a-/- BMDMs have decreased M1 responses and an enhanced M2 responses. Overexpression of mirn23a has the opposite effect, enhancing M1 and inhibiting M2 gene expression. Interestingly, expression of mirn23a miRNAs goes down with inflammatory stimulation and up with anti-inflammatory stimulation, suggesting that its regulation prevents locking macrophages into polarized states. M2 polarization of tumor-associated macrophages (TAMs) correlates with poor outcome for many tumors, so to determine if there was a functional consequence of mirn23a loss modulating immune cell polarization, we assayed syngeneic tumor growth in wild-type and mirn23a -/- mice. Consistent with the increased anti-inflammatory/immunosuppressive phenotype in vitro, mirn23a -/- mice inoculated with syngeneic tumor cells had worse outcomes compared with wild-type mice. Coinjecting tumor cells with mirn23a -/- BMDMs into wild-type mice phenocopied tumor growth in mirn23a -/- mice, supporting a critical role for mirn23a miRNAs in macrophage-mediated tumor immunity. Our data demonstrate that mirn23a regulates M1/M2 polarization and suggests that manipulation of mirn23a miRNA can be used to direct macrophage polarization to drive a desired immune response.

Genome editing demonstrates that the -5 kb Nanog enhancer regulates Nanog expression by modulating RNAPII initiation and/or recruitment.

Transcriptional enhancers have been defined by their ability to operate independent of distance and orientation in plasmid-based reporter assays of gene expression. At present, histone marks are used to identify and define enhancers but do not consider the endogenous role of an enhancer in the context of native chromatin. We employed a combination of genomic editing, single cell analyses, and sequencing approaches to investigate a Nanog-associated cis-regulatory element, which has been reported by others to be either an alternative promoter or a super-enhancer. We first demonstrate both distance and orientation independence in native chromatin, eliminating the issues raised with plasmid-based approaches. We next demonstrate that the dominant super-enhancer modulates Nanog globally and operates by recruiting and/or initiating RNA Polymerase II. Our studies have important implications to how transcriptional enhancers are defined and how they regulate gene expression.

Decitabine and vorinostat with FLAG chemotherapy in pediatric relapsed/refractory AML: Report from the therapeutic advances in childhood leukemia and lymphoma (TACL) consortium.

Survival outcomes for relapsed/refractory pediatric acute myeloid leukemia (R/R AML) remain dismal. Epigenetic changes can result in gene expression alterations which are thought to contribute to both leukemogenesis and chemotherapy resistance. We report results from a phase I trial with a dose expansion cohort investigating decitabine and vorinostat in combination with fludarabine, cytarabine, and G-CSF (FLAG) in pediatric patients with R/R AML [NCT02412475]. Thirty-seven patients enrolled with a median age at enrollment of 8.4 (range, 1-20) years. There were no dose limiting toxicities among the enrolled patients, including two patients with Down syndrome. The recommended phase 2 dose of decitabine in combination with vorinostat and FLAG was 10 mg/m2 . The expanded cohort design allowed for an efficacy evaluation and the overall response rate among 35 evaluable patients was 54% (16 complete response (CR) and 3 complete response with incomplete hematologic recovery (CRi)). Ninety percent of responders achieved minimal residual disease (MRD) negativity (<0.1%) by centralized flow cytometry and 84% (n = 16) successfully proceeded to hematopoietic stem cell transplant. Two-year overall survival was 75.6% [95%CI: 47.3%, 90.1%] for MRD-negative patients vs. 17.9% [95%CI: 4.4%, 38.8%] for those with residual disease (p < .001). Twelve subjects (34%) had known epigenetic alterations with 8 (67%) achieving a CR, 7 (88%) of whom were MRD negative. Correlative pharmacodynamics demonstrated the biologic activity of decitabine and vorinostat and identified specific gene enrichment signatures in nonresponding patients. Overall, this therapy was well-tolerated, biologically active, and effective in pediatric patients with R/R AML, particularly those with epigenetic alterations.

Rap1A, Rap1B, and ß-Adrenergic Signaling in Autologous HCT: A Randomized Controlled Trial of Propranolol.

Successful hematopoietic cell transplantation (HCT) depends on rapid engraftment of the progenitor and stem cells that will reestablish hematopoiesis. Rap1A and Rap1B are two closely related small GTPases that may affect platelet and neutrophil engraftment during HCT through their roles in cell adhesion and migration. ß-adrenergic signaling may regulate the participation of Rap1A and Rap1B in engraftment through their inhibition or activation. We conducted a correlative study of a randomized controlled trial evaluating the effects of the nonselective ß-antagonist propranolol on expression and prenylation of Rap1A and Rap1B during neutrophil and platelet engraftment in 25 individuals receiving an autologous HCT for multiple myeloma. Propranolol was administered for 1 week prior to and 4 weeks following HCT. Blood was collected 7 days (baseline) and 2 days (Day -2) before HCT, and 28 days after HCT (Day +28). Circulating polymorphonuclear cells (PMNC) were isolated and analyzed via immunoblotting to determine levels of prenylated and total Rap1A versus Rap1B. Twelve participants were randomized to the intervention and 13 to the control. Rap1A expression significantly correlated with Rap1B expression. Rap1B expression significantly correlated with slower platelet engraftment; however, this association was not observed in the propranolol-treated group. There were no significant associations between neutrophil engraftment and Rap1A or Rap1B expression. Post hoc exploratory analyses did not reveal an association between social health variables and Rap1A or Rap1B expression. This study identifies a greater regulatory role for Rap1B than Rap1A in platelet engraftment and suggests a possible role for ß-adrenergic signaling in modulating Rap1B function during HCT.

Degenerate minigene library analysis enables identification of altered branch point utilization by mutant splicing factor 3B1 (SF3B1).

Cancer-associated mutations of the core splicing factor 3 B1 (SF3B1) result in selection of novel 3' splice sites (3'SS), but precise molecular mechanisms of oncogenesis remain unclear. SF3B1 stabilizes the interaction between U2 snRNP and branch point (BP) on the pre-mRNA. It has hence been speculated that a change in BP selection is the basis for novel 3'SS selection. Direct quantitative determination of BP utilization is however technically challenging. To define BP utilization by SF3B1-mutant spliceosomes, we used an overexpression approach in human cells as well as a complementary strategy using isogenic murine embryonic stem cells with monoallelic K700E mutations constructed via CRISPR/Cas9-based genome editing and a dual vector homology-directed repair methodology. A synthetic minigene library with degenerate regions in 3' intronic regions (3.4 million individual minigenes) was used to compare BP usage of SF3B1K700E and SF3B1WT. Using this model, we show that SF3B1K700E spliceosomes utilize non-canonical sequence variants (at position -1 relative to BP adenosine) more frequently than wild-type spliceosomes. These predictions were confirmed using minigene splicing assays. Our results suggest a model of BP utilization by mutant SF3B1 wherein it is able to utilize non-consensus alternative BP sequences by stabilizing weaker U2-BP interactions.

Cyclin-dependent kinase 1 (CDK1) and CDK2 have opposing roles in regulating interactions of splicing factor 3B1 with chromatin.

Splicing factor 3B1 (SF3B1) is a core splicing protein that stabilizes the interaction between the U2 snRNA and the branch point in the mRNA target during splicing. SF3B1 is heavily phosphorylated at its N terminus and a substrate of cyclin-dependent kinases (CDKs). Although SF3B1 phosphorylation coincides with splicing catalysis, the functional significance of SF3B1 phosphorylation is largely undefined. Here, we show that SF3B1 phosphorylation follows a dynamic pattern during cell cycle progression that depends on CDK activity. SF3B1 is known to interact with chromatin, and we found that SF3B1 maximally interacts with nucleosomes during G1/S and that this interaction requires CDK2 activity. In contrast, SF3B1 disassociates from nucleosomes at G2/M, coinciding with a peak in CDK1-mediated SF3B1 phosphorylation. Thus, CDK1 and CDK2 appear to have opposing roles in regulating SF3B1-nucleosome interactions. Importantly, these interactions were modified by the presence and phosphorylation status of linker histone H1, particularly the H1.4 isoform. Performing genome-wide analysis of SF3B1-chromatin binding in synchronized cells, we observed that SF3B1 preferentially bound exons. Differences in SF3B1 chromatin binding to specific sites, however, did not correlate with changes in RNA splicing, suggesting that the SF3B1-nucleosome interaction does not determine cell cycle-dependent changes to mRNA splicing. Our results define a cell cycle stage-specific interaction between SF3B1 and nucleosomes that is mediated by histone H1 and depends on SF3B1 phosphorylation. Importantly, this interaction does not seem to be related to SF3B1's splicing function and, rather, points toward its potential role as a chromatin modifier.

Deletion of Tet proteins results in quantitative disparities during ESC differentiation partially attributable to alterations in gene expression.

BACKGROUND: The Tet protein family (Tet1, Tet2, and Tet3) regulate DNA methylation through conversion of 5-methylcytosine to 5-hydroxymethylcytosine which can ultimately result in DNA demethylation and play a critical role during early mammalian development and pluripotency. While multiple groups have generated knockouts combining loss of different Tet proteins in murine embryonic stem cells (ESCs), differences in genetic background and approaches has made it difficult to directly compare results and discern the direct mechanism by which Tet proteins regulate the transcriptome. To address this concern, we utilized genomic editing in an isogenic pluripotent background which permitted a quantitative, flow-cytometry based measurement of pluripotency in combination with genome-wide assessment of gene expression and DNA methylation changes. Our ultimate goal was to generate a resource of large-scale datasets to permit hypothesis-generating experiments. RESULTS: We demonstrate a quantitative disparity in the differentiation ability among Tet protein deletions, with Tet2 single knockout exhibiting the most severe defect, while loss of Tet1 alone or combinations of Tet genes showed a quantitatively intermediate phenotype. Using a combination of transcriptomic and epigenomic approaches we demonstrate an increase in DNA hypermethylation and a divergence of transcriptional profiles in pluripotency among Tet deletions, with loss of Tet2 having the most profound effect in undifferentiated ESCs. CONCLUSIONS: We conclude that loss of Tet2 has the most dramatic effect both on the phenotype of ESCs and the transcriptome compared to other genotypes. While loss of Tet proteins increased DNA hypermethylation, especially in gene promoters, these changes in DNA methylation did not correlate with gene expression changes. Thus, while loss of different Tet proteins alters DNA methylation, this change does not appear to be directly responsible for transcriptome changes. Thus, loss of Tet proteins likely regulates the transcriptome epigenetically both through altering 5mC but also through additional mechanisms. Nonetheless, the transcriptome changes in pluripotent Tet2-/- ESCs compared to wild-type implies that the disparities in differentiation can be partially attributed to baseline alterations in gene expression.

Nanog Expression in Embryonic Stem Cells - An Ideal Model System to Dissect Enhancer Function.

Embryonic stem cells (ESCs) are derived from the preimplantation embryo and can differentiate into virtually any other cell type (termed pluripotency), which is governed by lineage specific transcriptions factors (TFs) binding to cis regulatory elements (CREs) to mediate changes in gene expression. The reliance on transcriptional regulation to maintain pluripotency makes ESCs a valuable model to study the role of distal CREs such as enhancers in modulating gene expression to affect cell fate decisions. This review will highlight recent advance on transcriptional enhancers, focusing on studies performed in ESCs. In addition, we argue that the Nanog locus, which encodes for an ESC-critical TF, is particularly informative because it contains multiple co-regulated genes and enhancers in close proximity to one another. The unique landscape at Nanog permits the study of ongoing questions including whether multiple enhancers function additively versus synergistically, determinants of gene specificity, and cell-to-cell variability in gene expression.

CRISPR/Cas9-mediated conversion of human platelet alloantigen allotypes.

Human platelet alloantigens (HPAs) reside on functionally important platelet membrane glycoproteins and are caused by single nucleotide polymorphisms in the genes that encode them. Antibodies that form against HPAs are responsible for several clinically important alloimmune bleeding disorders, including fetal and neonatal alloimmune thrombocytopenia and posttransfusion purpura. The HPA-1a/HPA-1b alloantigen system, also known as the Pl(A1)/Pl(A2) polymorphism, is the most frequently implicated HPA among whites, and a single Leu33Pro amino acid polymorphism within the integrin ß3 subunit is responsible for generating the HPA-1a/HPA-1b alloantigenic epitopes. HPA-1b/b platelets, like those bearing other low-frequency platelet-specific alloantigens, are relatively rare in the population and difficult to obtain for purposes of transfusion therapy and diagnostic testing. We used CRISPR/Cas9 (clustered regularly interspaced short palindromic repeats/CRISPR associated protein 9) gene-editing technology to transform Leu33 (+) megakaryocytelike DAMI cells and induced pluripotent stem cells (iPSCs) to the Pro33 allotype. CD41(+) megakaryocyte progenitors derived from these cells expressed the HPA-1b (Pl(A2)) alloantigenic epitope, as reported by diagnostic NciI restriction enzyme digestion, DNA sequencing, and western blot analysis using HPA-1b-specific human maternal alloantisera. Application of CRISPR/Cas9 technology to genetically edit this and other clinically-important HPAs holds great potential for production of designer platelets for diagnostic, investigative, and, ultimately, therapeutic use.

A cell surfaceome map for immunophenotyping and sorting pluripotent stem cells.

Induction of a pluripotent state in somatic cells through nuclear reprogramming has ushered in a new era of regenerative medicine. Heterogeneity and varied differentiation potentials among induced pluripotent stem cell (iPSC) lines are, however, complicating factors that limit their usefulness for disease modeling, drug discovery, and patient therapies. Thus, there is an urgent need to develop nonmutagenic rapid throughput methods capable of distinguishing among putative iPSC lines of variable quality. To address this issue, we have applied a highly specific chemoproteomic targeting strategy for de novo discovery of cell surface N-glycoproteins to increase the knowledge-base of surface exposed proteins and accessible epitopes of pluripotent stem cells. We report the identification of 500 cell surface proteins on four embryonic stem cell and iPSCs lines and demonstrate the biological significance of this resource on mouse fibroblasts containing an oct4-GFP expression cassette that is active in reprogrammed cells. These results together with immunophenotyping, cell sorting, and functional analyses demonstrate that these newly identified surface marker panels are useful for isolating iPSCs from heterogeneous reprogrammed cultures and for isolating functionally distinct stem cell subpopulations.

Cohesin loss and MLL-AF9 are not synthetic lethal in murine hematopoietic stem and progenitor cells.

Objective As cohesin mutations are rarely found in MLL-rearranged acute myeloid leukemias, we investigated the potential synthetic lethality between cohesin mutations and MLL-AF9 using murine hematopoietic stem and progenitor cells. Results Contrary to our hypothesis, a complete loss of Stag2 or haploinsufficiency of Smc3 were well tolerated in MLL-AF9-expressing hematopoietic stem and progenitor cells. Minimal effect of cohesin subunit loss on the in vitro self-renewal of MLL-AF9-expressing cells was observed. Despite the differing mutational landscapes of cohesin-mutated and MLL fusion AMLs, previous studies showed that cohesin and MLL fusion mutations similarly drive abnormal self-renewal through HOXA gene upregulation. The utilization of a similar mechanism suggests that little selective pressure exists for the acquisition of cohesin mutations in AMLs expressing MLL fusions, explaining their lack of co-occurrence. Our results emphasize the importance of using genetic models to test suspected synthetic lethality and suggest that a lack of co-occurrence may instead point to a common mechanism of action between two mutations.

Transcriptomic diversity of innate lymphoid cells in human lymph nodes compared to BM and spleen.

Innate lymphoid cells (ILCs) are largely tissue-resident, mostly described within the mucosal tissues. However, their presence and functions in the human draining lymph nodes (LNs) are unknown. Our study unravels the tissue-specific transcriptional profiles of 47,287 CD127+ ILCs within the human abdominal and thoracic LNs. LNs contain a higher frequency of CD127+ ILCs than in BM or spleen. We define independent stages of ILC development, including EILP and pILC in the BM. These progenitors exist in LNs in addition to naïve ILCs (nILCs) that can differentiate into mature ILCs. We define three ILC1 and four ILC3 sub-clusters in the LNs. ILC1 and ILC3 subsets have clusters with high heat shock protein-encoding genes. We identify previously unrecognized regulons, including the BACH2 family for ILC1 and the ATF family for ILC3. Our study is the comprehensive characterization of ILCs in LNs, providing an in-depth understanding of ILC-mediated immunity in humans.

Histone Modifications and Their Contributions to Hypertension.

Essential hypertension, a multifaceted disorder, is a worldwide health problem. A complex network of genetic, epigenetic, physiological, and environmental components regulates blood pressure (BP), and any dysregulation of this network may result in hypertension. Growing evidence suggests a role for epigenetic factors in BP regulation. Any alterations in the expression or functions of these epigenetic regulators may dysregulate various determinants of BP, thereby promoting the development of hypertension. Histone posttranslational modifications are critical epigenetic regulators that have been implicated in hypertension. Several studies have demonstrated a clear association between the increased expression of some histone-modifying enzymes, especially HDACs (histone deacetylases), and hypertension. In addition, treatment with HDAC inhibitors lowers BP in hypertensive animal models, providing an excellent opportunity to design new drugs to treat hypertension. In this review, we discuss the potential contribution of different histone modifications to the regulation of BP.

N-MYC regulates cell survival via eIF4G1 in inv(16) acute myeloid leukemia.

N-MYC (encoded by MYCN) is a critical regulator of hematopoietic stem cell function. While the role of N-MYC deregulation is well established in neuroblastoma, the importance of N-MYC deregulation in leukemogenesis remains elusive. Here, we demonstrate that N-MYC is overexpressed in acute myeloid leukemia (AML) cells with chromosome inversion inv(16) and contributes to the survival and maintenance of inv(16) leukemia. We identified a previously unknown MYCN enhancer, active in multiple AML subtypes, essential for MYCN mRNA levels and survival in inv(16) AML cells. We also identified eukaryotic translation initiation factor 4 gamma 1 (eIF4G1) as a key N-MYC target that sustains leukemic survival in inv(16) AML cells. The oncogenic role of eIF4G1 in AML has not been reported before. Our results reveal a mechanism whereby N-MYC drives a leukemic transcriptional program and provides a rationale for the therapeutic targeting of the N-MYC/eIF4G1 axis in myeloid leukemia.