Smart-RRBS for single-cell methylome and transcriptome analysis

The integration of DNA methylation and transcriptional state within single cells is of broad interest. Several single-cell dual- and multi-omics approaches have been reported that enable further investigation into cellular heterogeneity, including the discovery and in-depth study of rare cell populations. Such analyses will continue to provide important mechanistic insights into the regulatory consequences of epigenetic modifications. We recently reported a new method for profiling the DNA methylome and transcriptome from the same single cells in a cancer research study. Here, we present details of the protocol and provide guidance on its utility. Our Smart-RRBS (reduced representation bisulfite sequencing) protocol combines Smart-seq2 and RRBS and entails physically separating mRNA from the genomic DNA. It generates paired epigenetic promoter and RNA-expression measurements for ~24% of protein-coding genes in a typical single cell. It also works for micro-dissected tissue samples comprising hundreds of cells. The protocol, excluding flow sorting of cells and sequencing, takes ~3 d to process up to 192 samples manually. It requires basic molecular biology expertise and laboratory equipment, including a PCR workstation with UV sterilization, a DNA fluorometer and a microfluidic electrophoresis system

Epigenetic evolution and lineage histories of chronic lymphocytic leukaemia

Genetic and epigenetic intra-tumoral heterogeneity cooperate to shape the evolutionary course of cancer1. Chronic lymphocytic leukaemia (CLL) is a highly informative model for cancer evolution as it undergoes substantial genetic diversification and evolution after therapy2,3. The CLL epigenome is also an important disease-defining feature4,5, and growing populations of cells in CLL diversify by stochastic changes in DNA methylation known as epimutations6. However, previous studies using bulk sequencing methods to analyse the patterns of DNA methylation were unable to determine whether epimutations affect CLL populations homogeneously. Here, to measure the epimutation rate at single-cell resolution, we applied multiplexed single-cell reduced-representation bisulfite sequencing to B cells from healthy donors and patients with CLL. We observed that the common clonal origin of CLL results in a consistently increased epimutation rate, with low variability in the cell-to-cell epimutation rate. By contrast, variable epimutation rates across healthy B cells reflect diverse evolutionary ages across the trajectory of B cell differentiation, consistent with epimutations serving as a molecular clock. Heritable epimutation information allowed us to reconstruct lineages at high-resolution with single-cell data, and to apply this directly to patient samples. The CLL lineage tree shape revealed earlier branching and longer branch lengths than in normal B cells, reflecting rapid drift after the initial malignant transformation and a greater proliferative history. Integration of single-cell bisulfite sequencing analysis with single-cell transcriptomes and genotyping confirmed that genetic subclones mapped to distinct clades, as inferred solely on the basis of epimutation information. Finally, to examine potential lineage biases during therapy, we profiled serial samples during ibrutinib-associated lymphocytosis, and identified clades of cells that were preferentially expelled from the lymph node after treatment, marked by distinct transcriptional profiles. The single-cell integration of genetic, epigenetic and transcriptional information thus charts the lineage history of CLL and its evolution with therapy

In Vitro Megakaryocyte Differentiation and Proplatelet Formation in Ph-Negative Classical Myeloproliferative Neoplasms: Distinct Patterns in the Different Clinical Phenotypes

Background: Ph-negative myeloproliferative neoplasms (MPNs) are clonal disorders that include primary myelofibrosis (PMF), polycythemia vera (PV) and essential thrombocythemia (ET). Although the pathogenesis of MPNs is still incompletely understood, an involvement of the megakaryocyte lineage is a distinctive feature. Methodology/Principal Findings: We analyzed the in vitro megakaryocyte differentiation and proplatelet formation in 30 PMF, 8 ET, 8 PV patients, and 17 healthy controls (CTRL). Megakaryocytes were differentiated from peripheral blood CD34+ or CD45+ cells in the presence of thrombopoietin. Megakaryocyte output was higher in MPN patients than in CTRL with no correlation with the JAK2 V617F mutation. PMF-derived megakaryocytes displayed nuclei with a bulbous appearance, were smaller than ET- or PV-derived megakaryocytes and formed proplatelets that presented several structural alterations. In contrast, ET- and PV-derived megakaryocytes produced more proplatelets with a striking increase in bifurcations and tips compared to both control and PMF. Proplatelets formation was correlated with platelet counts in patient peripheral blood. Patients with pre-fibrotic PMF had a pattern of megakaryocyte proliferation and proplatelet formation that was similar to that of fibrotic PMF and different from that of ET. Conclusions/Significance: In conclusion, MPNs are associated with high megakaryocyte proliferative potential. Profound differences in megakaryocyte morphology and proplatelet formation distinguish PMF, both fibrotic and prefibrotic, from ET and PV

Functional Studies on the IBD Susceptibility Gene IL23R Implicate Reduced Receptor Function in the Protective Genetic Variant R381Q

Genome-wide association studies (GWAS) in several populations have demonstrated significant association of the IL23R gene with IBD (Crohn's disease (CD) and ulcerative colitis (UC)) and psoriasis, suggesting that perturbation of the IL-23 signaling pathway is relevant to the pathophysiology of these diseases. One particular variant, R381Q (rs11209026), confers strong protection against development of CD. We investigated the effects of this variant in primary T cells from healthy donors carrying IL23RR381 and IL23RQ381 haplotypes. Using a proprietary anti-IL23R antibody, ELISA, flow cytometry, phosphoflow and real-time RT-PCR methods, we examined IL23R expression and STAT3 phosphorylation and activation in response to IL-23. IL23RQ381 was associated with reduced STAT3 phosphorylation upon stimulation with IL-23 and decreased number of IL-23 responsive T-cells. We also observed slightly reduced levels of proinflammatory cytokine secretion in IL23RQ381 positive donors. Our study shows conclusively that IL23RQ381 is a loss-of-function allele, further strengthening the implication from GWAS results that the IL-23 pathway is pathogenic in human disease. This data provides an explanation for the protective role of R381Q in CD and may lead to the development of improved therapeutics for autoimmune disorders like CD

Novel mutations and their functional and clinical relevance in myeloproliferative neoplasms: JAK2, MPL, TET2, ASXL1, CBL, IDH and IKZF1

Myeloproliferative neoplasms (MPNs) originate from genetically transformed hematopoietic stem cells that retain the capacity for multilineage differentiation and effective myelopoiesis. Beginning in early 2005, a number of novel mutations involving Janus kinase 2 (JAK2), Myeloproliferative Leukemia Virus (MPL), TET oncogene family member 2 (TET2), Additional Sex Combs-Like 1 (ASXL1), Casitas B-lineage lymphoma proto-oncogene (CBL), Isocitrate dehydrogenase (IDH) and IKAROS family zinc finger 1 (IKZF1) have been described in BCR-ABL1-negative MPNs. However, none of these mutations were MPN specific, displayed mutual exclusivity or could be traced back to a common ancestral clone. JAK2 and MPL mutations appear to exert a phenotype-modifying effect and are distinctly associated with polycythemia vera, essential thrombocythemia and primary myelofibrosis; the corresponding mutational frequencies are ∼99, 55 and 65% for JAK2 and 0, 3 and 10% for MPL mutations. The incidence of TET2, ASXL1, CBL, IDH or IKZF1 mutations in these disorders ranges from 0 to 17% these latter mutations are more common in chronic (TET2, ASXL1, CBL) or juvenile (CBL) myelomonocytic leukemias, mastocytosis (TET2), myelodysplastic syndromes (TET2, ASXL1) and secondary acute myeloid leukemia, including blast-phase MPN (IDH, ASXL1, IKZF1). The functional consequences of MPN-associated mutations include unregulated JAK-STAT (Janus kinase/signal transducer and activator of transcription) signaling, epigenetic modulation of transcription and abnormal accumulation of oncoproteins. However, it is not clear as to whether and how these abnormalities contribute to disease initiation, clonal evolution or blastic transformation

The S505A thrombopoietin receptor mutation in childhood hereditary thrombocytosis and essential thrombocythemia is S505N: single letter amino acid code matters.

Dear Editor, Myeloproliferative neoplasms (MPNs) are rare diseases in children [1] when compared to adults [2]. Targeted sequencing revealed differences in the mutational landscape of pediatric and adult MPNs [1]. Children exhibit a lower frequency of mutations in the MPN driver genes JAK2, MPL, and CALR, and a much higher proportion of children exhibit no mutation among the 104 classically involved genes; these children exhibiting a trend towards essential thrombocythemia (ET) [1]. We were surprised to read about an activating thrombopoietin receptor (TpoR/Mpl) S505A transmembrane domain (TM) mutation in childhood hereditary thrombocytosis (HT) and ET [3–5]. [...

Ordered chromatin changes and human X chromosome reactivation by cell fusion-mediated pluripotent reprogramming.

Erasure of epigenetic memory is required to convert somatic cells towards pluripotency. Reactivation of the inactive X chromosome (Xi) has been used to model epigenetic reprogramming in mouse, but human studies are hampered by Xi epigenetic instability and difficulties in tracking partially reprogrammed iPSCs. Here we use cell fusion to examine the earliest events in the reprogramming-induced Xi reactivation of human female fibroblasts. We show that a rapid and widespread loss of Xi-associated H3K27me3 and XIST occurs in fused cells and precedes the bi-allelic expression of selected Xi-genes by many heterokaryons (30-50%). After cell division, RNA-FISH and RNA-seq analyses confirm that Xi reactivation remains partial and that induction of human pluripotency-specific XACT transcripts is rare (1%). These data effectively separate pre- and post-mitotic events in reprogramming-induced Xi reactivation and reveal a complex hierarchy of epigenetic changes that are required to reactivate the genes on the human Xi chromosome

Ordered chromatin changes and human X chromosome reactivation by cell fusion-mediated pluripotent reprogramming

Erasure of epigenetic memory is required to convert somatic cells towards pluripotency. Reactivation of the inactive X chromosome (Xi) has been used to model epigenetic reprogramming in mouse, but human studies are hampered by Xi epigenetic instability and difficulties in tracking partially reprogrammed iPSCs. Here we use cell fusion to examine the earliest events in the reprogramming-induced Xi reactivation of human female fibroblasts. We show that a rapid and widespread loss of Xi-associated H3K27me3 and XIST occurs in fused cells and precedes the bi-allelic expression of selected Xi-genes by many heterokaryons (30-50%). After cell division, RNA-FISH and RNA-seq analyses confirm that Xi reactivation remains partial and that induction of human pluripotency-specific XACT transcripts is rare (1%). These data effectively separate pre- and post-mitotic events in reprogramming-induced Xi reactivation and reveal a complex hierarchy of epigenetic changes that are required to reactivate the genes on the human Xi chromosome