CytoTree: an R/Bioconductor package for analysis and visualization of flow and mass cytometry data.

BACKGROUND: The rapidly increasing dimensionality and throughput of flow and mass cytometry data necessitate new bioinformatics tools for analysis and interpretation, and the recently emerging single-cell-based algorithms provide a powerful strategy to meet this challenge. RESULTS: Here, we present CytoTree, an R/Bioconductor package designed to analyze and interpret multidimensional flow and mass cytometry data. CytoTree provides multiple computational functionalities that integrate most of the commonly used techniques in unsupervised clustering and dimensionality reduction and, more importantly, support the construction of a tree-shaped trajectory based on the minimum spanning tree algorithm. A graph-based algorithm is also implemented to estimate the pseudotime and infer intermediate-state cells. We apply CytoTree to several examples of mass cytometry and time-course flow cytometry data on heterogeneity-based cytology and differentiation/reprogramming experiments to illustrate the practical utility achieved in a fast and convenient manner. CONCLUSIONS: CytoTree represents a versatile tool for analyzing multidimensional flow and mass cytometry data and to producing heuristic results for trajectory construction and pseudotime estimation in an integrated workflow.

The short isoform of the long-type PML-RARA fusion gene in acute promyelocytic leukaemia lacks sensitivity to all-trans-retinoic acid.

Alternative splicing is associated with human disease. In acute promyelocytic leukaemia (APL) patients with the long (L)-type promyelocytic leukaemia-retinoic acid receptor α fusion gene (PML-RARA), three alternative splicing isoforms can be detected: E5(+)E6(+), E5(-)E6(+), and E5(-)E6(-). This study is the first to demonstrate that alternative splicing of L-type PML-RARA is associated with time to achieve complete remission (CR) in APL. Higher expression of the E5(-)E6(-) isoform, the short isoform, was related to longer time to achieve CR. Each isoform was constructed into recombinant lentiviral vector and transfected into U937 cells. Compared with the E5(-)E6(+) and E5(+)E6(+) groups, the U937 cells with E5(-)E6(-) showed lower sensitivity to all-trans-retinoic acid treatment.

A systematic survey of LU domain-containing proteins reveals a novel human gene, LY6A, which encodes the candidate ortholog of mouse Ly-6A/Sca-1 and is aberrantly expressed in pituitary tumors.

The Ly-6 and uPAR (LU) domain-containing proteins represent a large family of cell-surface markers. In particular, mouse Ly-6A/Sca-1 is a widely used marker for various stem cells; however, its human ortholog is missing. In this study, based on a systematic survey and comparative genomic study of mouse and human LU domain-containing proteins, we identified a previously unannotated human gene encoding the candidate ortholog of mouse Ly-6A/Sca-1. This gene, hereby named LY6A, reversely overlaps with a lncRNA gene in the majority of exonic sequences. We found that LY6A is aberrantly expressed in pituitary tumors, but not in normal pituitary tissues, and may contribute to tumorigenesis. Similar to mouse Ly-6A/Sca-1, human LY6A is also upregulated by interferon, suggesting a conserved transcriptional regulatory mechanism between humans and mice. We cloned the full-length LY6A cDNA, whose encoded protein sequence, domain architecture, and exon-intron structures are all well conserved with mouse Ly-6A/Sca-1. Ectopic expression of the LY6A protein in cells demonstrates that it acts the same as mouse Ly-6A/Sca-1 in their processing and glycosylphosphatidylinositol anchoring to the cell membrane. Collectively, these studies unveil a novel human gene encoding a candidate biomarker and provide an interesting model gene for studying gene regulatory and evolutionary mechanisms.

Detection of MPL exon10 mutations in 103 Chinese patients with JAK2V617F-negative myeloproliferative neoplasms.

JAK2V617F mutation has been reported in 90% of patients with polycythemia vera (PV) and about 50% of patients with essential thromobocythemia (ET) and primary myelofibrosis (PMF). Recently, acquired mutations in the transmembrane-juxtamembrane region of MPL (MPLW515 mutations) have been reported in approximately 5% of JAK2V617F-negative PMF and about 1% of all cases of ET. MPL is the receptor for thrombopoietin that regulates the production of platelets by bone marrow. It is likely that some mutations more closely related to ET in MPL exon10 may have been missed by current assays. We inferred that there might be other mutations in MPL exon10 for MPN patients in addition to MPLW515 mutations. To investigate its mutation types and prevalence in Chinese patients with myeloproliferative neoplasms (MPN), we performed mutation detection on MPL exon10 in 103 JAK2V617F-negative MPN patients by single strand conformation polymorphism (SSCP) and allele-specific PCR (AS-PCR) combined with sequencing. As a result, one previously unrecognized MPL mutation (12-bp in-frame insertion) was identified in one patient with ET in addition to an MPLW515K mutation identified in one PMF patient. This confirms our hypothesis that BCR/ABL negative and JAK2V617F-negative MPN patients have other mutations besides W515 mutation in MPL exon10 and mutations other than single nucleotide exchange also exist. In addition, MPL mutation was associated with Chinese MPN patients.

Selective and competitive functions of the AAR and UPR pathways in stress-induced angiogenesis.

The amino acid response (AAR) and unfolded protein response (UPR) pathways converge on eIF2α phosphorylation, which is catalyzed by Gcn2 and Perk, respectively, under different stresses. This close interconnection makes it difficult to specify different functions of AAR and UPR. Here, we generated a zebrafish model in which loss of threonyl-tRNA synthetase (Tars) induces angiogenesis dependent on Tars aminoacylation activity. Comparative transcriptome analysis of the tars-mutant and wild-type embryos with/without Gcn2- or Perk-inhibition reveals that only Gcn2-mediated AAR is activated in the tars-mutants, whereas Perk functions predominantly in normal development. Mechanistic analysis shows that, while a considerable amount of eIF2α is normally phosphorylated by Perk, the loss of Tars causes an accumulation of uncharged tRNAThr, which in turn activates Gcn2, leading to phosphorylation of an extra amount of eIF2α. The partial switchover of kinases for eIF2α largely overwhelms the functions of Perk in normal development. Interestingly, although inhibition of Gcn2 and Perk in this stress condition both can reduce the eIF2α phosphorylation levels, their functional consequences in the regulation of target genes and in the rescue of the angiogenic phenotypes are dramatically different. Indeed, genetic and pharmacological manipulations of these pathways validate that the Gcn2-mediated AAR, but not the Perk-mediated UPR, is required for tars-deficiency induced angiogenesis. Thus, the interconnected AAR and UPR pathways differentially regulate angiogenesis through selective functions and mutual competitions, reflecting the specificity and efficiency of multiple stress response pathways that evolve integrally to enable an organism to sense/respond precisely to various types of stresses.

Anti-Tyro3 IgG Associates with Disease Activity and Reduces Efferocytosis of Macrophages in New-Onset Systemic Lupus Erythematosus.

BACKGROUND: Systemic lupus erythematosus (SLE) is a disease characterized by the production of a large number of autoantibodies. Defected phagocytosis of macrophage plays an important role in innate immunity in the pathogenesis of SLE. Tyro3 is a receptor responsible for the recognition of apoptotic cells during efferocytosis by macrophages. To investigate the role of Tyro3 receptor in macrophages' efferocytosis of apoptotic cells in SLE, we aimed to reveal the clinical relevance and impact of Tyro3 autoantibody on SLE. METHODS: The serum levels of IgG-type autoantibody against Tyro3 receptor were detected in new-onset, treatment-naïve SLE patients (n = 70), rheumatoid arthritis (RA) (n = 24), primary Sjögren's Syndrome (pSS) (n = 21), and healthy controls (HCs) (n = 70) using enzyme-linked immunosorbent assay (ELISA). The effects of purified Tyro3 autoantibody from SLE patients on the efferocytosis of human monocyte-derived macrophages were measured by flow cytometry and immunofluorescence. RESULTS: The serum levels of IgG-type autoantibody against Tyro3 receptor were significantly elevated in patients with SLE compared to RA, pSS, and HCs (all p < 0.0001). The levels of anti-Tyro3 IgG were positively associated with the SLE disease activity index (SLEDAI) score (r = 0.254, p = 0.034), erythrocyte sedimentation rate (ESR) (r = 0.430, p < 0.001), C-reactive protein (CRP) (r = 0.246, p = 0.049), and immunoglobulin G (IgG) (r = 0.408, p = 0.001) and negatively associated with haemoglobin (Hb) (r = -0.294, p = 0.014). ROC curves illustrated that the anti-Tyro3 antibody could differentiate patients with SLE from HCs. Furthermore, flow cytometry and immunofluorescence demonstrated that purified anti-Tyro3 IgG inhibited the efferocytosis of macrophages (p = 0.004 and 0.044, respectively) compared with unconjugated human IgG. CONCLUSIONS: These observations indicated that autoantibody against Tyro3 was associated with disease activity and could impair efferocytosis of macrophages. It might be a potential novel disease biomarker and might be involved in the pathogenesis of SLE.

Setd2 deficiency impairs hematopoietic stem cell self-renewal and causes malignant transformation.

The histone H3 lysine 36 methyltransferase SETD2 is frequently mutated in various cancers, including leukemia. However, there has not been any functional model to show the contribution of SETD2 in hematopoiesis or the causal role of SETD2 mutation in tumorigenesis. In this study, using a conditional Setd2 knockout mouse model, we show that Setd2 deficiency skews hematopoietic differentiation and reduces the number of multipotent progenitors; although the number of phenotypic hematopoietic stem cells (HSCs) in Setd2-deleted mice is unchanged, functional assays, including serial BM transplantation, reveal that the self-renewal and competitiveness of HSCs are impaired. Intriguingly, Setd2-deleted HSCs, through a latency period, can acquire abilities to overcome the growth disadvantage and eventually give rise to hematopoietic malignancy characteristic of myelodysplastic syndrome. Gene expression profile of Setd2-deleted hematopoietic stem/progenitor cells (HSPCs) partially resembles that of Dnmt3a/Tet2 double knockout HSPCs, showing activation of the erythroid transcription factor Klf1-related pathway, which plays an important role in hematopoietic malignant transformation. Setd2 deficiency also induces DNA replication stress in HSCs, as reflected by an activated E2F gene regulatory network and repressed expression of the ribonucleotide reductase subunit Rrm2b, which results in proliferation and cell cycle abnormalities and genomic instability, allowing accumulation of secondary mutation(s) that synergistically contributes to tumorigenesis. Thus, our results demonstrate that Setd2 is required for HSC self-renewal, and provide evidence supporting the causal role of Setd2 deficiency in tumorigenesis. The underlying mechanism shall advance our understanding of epigenetic regulation of cancer and provide potential new therapeutic targets.

Chemical transdifferentiation: closer to regenerative medicine.

Cell transdifferentiation, which directly switches one type of differentiated cells into another cell type, is more advantageous than cell reprogramming to generate pluripotent cells and differentiate them into functional cells. This process is crucial in regenerative medicine. However, the cell-converting strategies, which mainly depend on the virus-mediated expression of exogenous genes, have clinical safety concerns. Small molecules with compelling advantages are a potential alternative in manipulating cell fate conversion. In this review, we briefly retrospect the nature of cell transdifferentiation and summarize the current developments in the research of small molecules in promoting cell conversion. Particularly, we focus on the complete chemical compound-induced cell transdifferentiation, which is closer to the clinical translation in cell therapy. Despite these achievements, the mechanisms underpinning chemical transdifferentiation remain largely unknown. More importantly, identifying drugs that induce resident cell conversion in vivo to repair damaged tissue remains to be the end-goal in current regenerative medicine.

[Establishment of diagnostic model to monitor minimal residual disease of acute promyelocytic leukemia by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry].

OBJECTIVE: To screen the potential protein biomarkers in minimal residual disease (MRD) of the acute promyelocytic leukemia (APL) by comparison of differentially expressed serum protein between APL patients at diagnosis and after complete remission (CR) and healthy controls, and to establish and verify a diagnostic model. METHODS: Serum proteins from 36 cases of primary APL, 29 cases of APL during complete remission and 32 healthy controls were purified by magnetic beads and then analyzed by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS). The spectra were analyzed statistically using FlexAnalysis(TM) and ClinProt(TM) software. RESULTS: Two prediction model of primary APL/healthy control, primary APL/APL CR were developed. Thirty four statistically significant peptide peaks were obtained with the m/z value ranging from 1000 to 10 000 (P < 0.001) in primary APL/healthy control model. Seven statistically significant peptide peaks were obtained in primary APL/APL CR model (P < 0.001). Comparison of the protein profiles between the two models, three peptides with m/z 4642, 7764 and 9289 were considered as the protein biomarker of APL MRD. A diagnostic pattern for APL CR using m/z 4642 and 9289 was established. Blind validation yielded correct classification of 6 out of 8 cases. CONCLUSIONS: The MALDI-TOF MS analysis of APL patients serum protein can be used as a promising dynamic method for MRD detection and the two peptides with m/z 4642 and 9289 may be better biomarkers.

Identification of a novel circularized transcript of the AML1 gene.

The AML1 gene is an essential transcription factor regulating the differentiation of hematopoietic stem cells into mature blood cells. Though at least 12 different alternatively spliced AML1 mRNAs are generated, three splice variants (AML1a, AML1b and AML1c) have been characterized. Here, using the reverse transcription-polymerase chain reaction with outward-facing primers, we identified a novel non-polyadenylated transcript from the AML1 gene, with exons 5 and 6 scrambled. The novel transcript resisted RNase R digestion, indicating it is a circular RNA structure that may originate from products of mRNA alternative splicing. The expression of the novel transcript in different cells or cell lines of human and a number of other species matched those of the canonical transcripts. The discovery provides additional evidence that circular RNA could stably exist in vivo in human, and may also help to understand the mechanism of the regulation of the AML1 gene transcription.