Global gene expression reveals a set of new genes involved in the modification of cells during erythroid differentiation

Differentiation is accompanied by a progressive loss in transcriptional memory

BACKGROUND

Cell differentiation requires the integration of two opposite processes, a stabilizing cellular memory, especially at the transcriptional scale, and a burst of gene expression variability which follows the differentiation induction. Therefore, the actual capacity of a cell to undergo phenotypic change during a differentiation process relies upon a modification in this balance which favors change-inducing gene expression variability. However, there are no experimental data providing insight on how fast the transcriptomes of identical cells would diverge on the scale of the very first two cell divisions during the differentiation process.

RESULTS

In order to quantitatively address this question, we developed different experimental methods to recover the transcriptomes of related cells, after one and two divisions, while preserving the information about their lineage at the scale of a single cell division. We analyzed the transcriptomes of related cells from two differentiation biological systems (human CD34+ cells and T2EC chicken primary erythrocytic progenitors) using two different single-cell transcriptomics technologies (scRT-qPCR and scRNA-seq).

CONCLUSIONS

We identified that the gene transcription profiles of differentiating sister cells are more similar to each other than to those of non-related cells of the same type, sharing the same environment and undergoing similar biological processes. More importantly, we observed greater discrepancies between differentiating sister cells than between self-renewing sister cells. Furthermore, a progressive increase in this divergence from first generation to second generation was observed when comparing differentiating cousin cells to self renewing cousin cells. Our results are in favor of a gradual erasure of transcriptional memory during the differentiation process.

Peroxiredoxins in erythrocytes: far beyond the antioxidant role

The red blood cells (RBCs) are essential to transport oxygen (O2) and nutrients throughout the human body. Changes in the structure or functioning of the erythrocytes can lead to several deficiencies, such as hemolytic anemias, in which an increase in reactive oxidative species generation is involved in the pathophysiological process, playing a significant role in the severity of several clinical manifestations. There are important lines of defense against the damage caused by oxidizing molecules. Among the antioxidant molecules, the enzyme peroxiredoxin (Prx) has the higher decomposition power of hydrogen peroxide, especially in RBCs, standing out because of its abundance. This review aimed to present the recent findings that broke some paradigms regarding the three isoforms of Prxs found in RBC (Prx1, Prx2, and Prx6), showing that in addition to their antioxidant activity, these enzymes may have supplementary roles in transducing peroxide signals, as molecular chaperones, protecting from membrane damage, and maintenance of iron homeostasis, thus contributing to the overall survival of human RBCs, roles that seen to be disrupted in hemolytic anemia conditions.

Thirty Years with ERH: An mRNA Splicing and Mitosis Factor Only or Rather a Novel Genome Integrity Protector?

ERH is a 100 to about 110 aa nuclear protein with unique primary and three-dimensional structures that are very conserved from simple eukaryotes to humans, albeit some species have lost its gene, with most higher fungi being a noteworthy example. Initially, studies on Drosophila melanogaster implied its function in pyrimidine metabolism. Subsequently, research on Xenopus laevis suggested that it acts as a transcriptional repressor. Finally, studies in humans pointed to a role in pre-mRNA splicing and in mitosis but further research, also in Caenorhabditis elegans and Schizosaccharomyces pombe, demonstrated its much broader activity, namely involvement in the biogenesis of mRNA, and miRNA, piRNA and some other ncRNAs, and in repressive heterochromatin formation. ERH interacts with numerous, mostly taxon-specific proteins, like Mmi1 and Mei2 in S. pombe, PID-3/PICS-1, TOST-1 and PID-1 in C. elegans, and DGCR8, CIZ1, PDIP46/SKAR and SAFB1/2 in humans. There are, however, some common themes in this wide range of processes and partners, such as: (a) ERH homodimerizes to form a scaffold for several complexes involved in the metabolism of nucleic acids, (b) all these RNAs are RNA polymerase II transcripts, (c) pre-mRNAs, whose splicing depends on ERH, are enriched in transcripts of DNA damage response and DNA metabolism genes, and (d) heterochromatin is formed to silence unwanted transcription, e.g., from repetitive elements. Thus, it seems that ERH has been adopted for various pathways that serve to maintain genome integrity.

Role of Cytokine-Inducible SH2 Domain-Containing (CISH) Protein in the Regulation of Erythropoiesis

The cytokine-inducible SH2 domain-containing (CISH) protein was the first member of the suppressor of cytokine signaling (SOCS) family of negative feedback regulators discovered, being identified in vitro as an inducible inhibitor of erythropoietin (EPO) signaling. However, understanding of the physiological role played by CISH in erythropoiesis has remained limited. To directly assess the function of CISH in this context, mice deficient in CISH were characterized with respect to developmental, steady-state, and EPO-induced erythropoiesis. CISH was strongly expressed in the fetal liver, but CISH knockout (KO) mice showed only minor disruption of primitive erythropoiesis. However, adults exhibited mild macrocytic anemia coincident with subtle perturbation particularly of bone marrow erythropoiesis, with EPO-induced erythropoiesis blunted in the bone marrow of KO mice but enhanced in the spleen. Cish was expressed basally in the bone marrow with induction following EPO stimulation in bone marrow and spleen. Overall, this study indicates that CISH participates in the control of both basal and EPO-induced erythropoiesis in vivo.

Identification of key pathways, genes and immune cell infiltration in hypoxia of high-altitude acclimatization via meta-analysis and integrated bioinformatics analysis

Background: For individuals acutely exposed to high-altitude regions, environmental hypobaric hypoxia induces several physiological or pathological responses, especially immune dysfunction. Therefore, hypoxia is a potentially life-threatening factor, which has closely related to high-altitude acclimatization. However, its specific molecular mechanism is still unclear. Methods: The four expression profiles about hypoxia and high altitude were downloaded from the Gene Expression Omnibus database in this study. Meta-analysis of GEO datasets was performed by NetworkAnalyst online tool. Kyoto Encyclopedia of Genes and Genomes (KEGG), Gene ontology (GO) enrichment analysis, and visualization were performed using R (version 4.1.3) software, respectively. The CIBERSORT analysis was conducted on GSE46480 to examine immune cell infiltration. In addition, we experimentally verified the bioinformatics analysis with qRT-PCR. Results: The meta-analysis identified 358 differentially expressed genes (DEGs), with 209 upregulated and 149 downregulated. DEGs were mostly enriched in biological processes and pathways associated with hypoxia acclimatization at high altitudes, according to both GO and KEGG enrichment analyses. ERH, VBP1, BINP3L, TOMM5, PSMA4, and POLR2K were identified by taking intersections of the DEGs between meta-analysis and GSE46480 and verified by qRT-PCR experiments, which were inextricably linked to hypoxia. Immune infiltration analysis showed significant differences in immune cells between samples at sea level and high altitudes. Conclusion: Identifying the DEGs and pathways will improve our understanding of immune function during high-altitude hypoxia at a molecular level. Targeting hypoxia-sensitive pathways in immune cells is interesting in treating high-altitude sickness. This study provides support for further research on high-altitude acclimatization.

ERH Gene and Its Role in Cancer Cells

Cancer is a major public health problem worldwide. Studies on oncogenes and tumor-targeted therapies have become an important part of cancer treatment development. In this review, we summarize and systematically introduce the gene enhancer of rudimentary homolog (ERH), which encodes a highly conserved small molecule protein. ERH mainly exists as a protein partner in human cells. It is involved in pyrimidine metabolism and protein complexes, acts as a transcriptional repressor, and participates in cell cycle regulation. Moreover, it is involved in DNA damage repair, mRNA splicing, the process of microRNA hairpins as well as erythroid differentiation. There are many related studies on the role of ERH in cancer cells; however, there are none on tumor-targeted therapeutic drugs or related therapies based on the expression of ERH. This study will provide possible directions for oncologists to further their research studies in this field.

LSD1-mediated repression of GFI1 super-enhancer plays an essential role in erythroleukemia

Super-enhancers (SEs) consist of enhancer clusters with abundant binding of transcription factors (TFs) and cofactors. LSD1 is a histone modifier that eliminates SE activity. However, whether SE suppression by LSD1 is associated with leukemogenesis remains unknown. In erythro-megakaryocyte lineage leukemia cells, activation of the SE of GFI1 (GFI1-SE) is related to induction of myeloid differentiation by LSD1 inhibitors NCD38 and NCD25 and to their antileukemia effect. Although functional TF-motifs were concentrated in an evolutionally conserved area, NCD38 barely induced additional TF recruitment. Instead, the transcription cofactors including LSD1, CoREST, HDAC1, and HDAC2 were evicted from GFI1-SE. Deletion of GFI1-SE impaired induction of myeloid differentiation by NCD38 and NCD25 in erythroleukemia cells. Gene set enrichment analysis revealed that the GFI1-SE deletion impaired NCD38-induced programs related to granulocyte differentiation and the CEBPA network, but restored NCD38-suppressed programs related to erythroid development, GATA1 targets, and acute myeloid leukemia (AML) clusters including FAB subtype M6 and AML with myelodysplastic syndrome-related chromosomal abnormalities. Ontologies of genes whose expression changes by NCD38 were canceled due to the GFI1-SE deletion showed enrichment in AML and neutropenia signatures. Collectively, our data suggest that sustainable repression of GFI1-SE by LSD1 is essential for sustenance of erythroleukemia cells.

A systems biology pipeline identifies regulatory networks for stem cell engineering

A major challenge for stem cell engineering is achieving a holistic understanding of the molecular networks and biological processes governing cell differentiation. To address this challenge, we describe a computational approach that combines gene expression analysis, previous knowledge from proteomic pathway informatics and cell signaling models to delineate key transitional states of differentiating cells at high resolution. Our network models connect sparse gene signatures with corresponding, yet disparate, biological processes to uncover molecular mechanisms governing cell fate transitions. This approach builds on our earlier CellNet and recent trajectory-defining algorithms, as illustrated by our analysis of hematopoietic specification along the erythroid lineage, which reveals a role for the EGF receptor family member, ErbB4, as an important mediator of blood development. We experimentally validate this prediction and perturb the pathway to improve erythroid maturation from human pluripotent stem cells. These results exploit an integrative systems perspective to identify new regulatory processes and nodes useful in cell engineering.

Global analysis of erythroid cells redox status reveals the involvement of Prdx1 and Prdx2 in the severity of beta thalassemia

β-thalassemia is a worldwide distributed monogenic red cell disorder, characterized by an absent or reduced beta globin chain synthesis. The unbalance of alpha-gamma chain and the presence of pathological free iron promote severe oxidative damage, playing crucial a role in erythrocyte hemolysis, exacerbating ineffective erythropoiesis and decreasing the lifespan of red blood cells (RBC). Catalase, glutathione peroxidase and peroxiredoxins act together to protect RBCs from hydrogen peroxide insult. Among them, peroxiredoxins stand out for their overall abundance and reactivity. In RBCs, Prdx2 is the third most abundant protein, although Prdxs 1 and 6 isoforms are also found in lower amounts. Despite the importance of these enzymes, Prdx1 and Prdx2 may have their peroxidase activity inactivated by hyperoxidation at high hydroperoxide concentrations, which also promotes the molecular chaperone activity of these proteins. Some studies have demonstrated the importance of Prdx1 and Prdx2 for the development and maintenance of erythrocytes in hemolytic anemia. Now, we performed a global analysis comparatively evaluating the expression profile of several antioxidant enzymes and their physiological reducing agents in patients with beta thalassemia intermedia (BTI) and healthy individuals. Furthermore, increased levels of ROS were observed not only in RBC, but also in neutrophils and mononuclear cells of BTI patients. The level of transcripts and the protein content of Prx1 were increased in reticulocyte and RBCs of BTI patients and the protein content was also found to be higher when compared to beta thalassemia major (BTM), suggesting that this peroxidase could cooperate with Prx2 in the removal of H₂O₂. Furthermore, Prdx2 production is highly increased in RBCs of BTM patients that present high amounts of hyperoxidized species. A significant increase in the content of Trx1, Srx1 and Sod1 in RBCs of BTI patients suggested protective roles for these enzymes in BTI patients. Finally, the upregulation of Nrf2 and Keap1 transcription factors found in BTI patients may be involved in the regulation of the antioxidant enzymes analyzed in this work.

RPL27A is a target of miR-595 and may contribute to the myelodysplastic phenotype through ribosomal dysgenesis

We investigated the functional consequences following deletion of a microRNA (miR) termed miR-595 which resides on chromosome 7q and is localised within one of the commonly deleted regions identified for Myelodysplasia (MDS) with monosomy 7 (−7)/isolated loss of 7q (7q-). We identified several targets for miR-595, including a large ribosomal subunit protein RPL27A. RPL27A downregulation induced p53 activation, apoptosis and inhibited proliferation. Moreover, p53-independent effects were additionally identified secondary to a reduction in the ribosome subunit 60s. We confirmed that RPL27A plays a pivotal role in the maintenance of nucleolar integrity and ribosomal synthesis/maturation. Of note, RPL27A overexpression, despite showing no significant effects on p53 mRNA levels, did in fact enhance cellular proliferation. In normal CD34+ cells, RPL27A knockdown preferentially blocked erythroid proliferation and differentiation. Lastly, we show that miR-595 expression appears significantly downregulated in the majority of primary samples derived from MDS patients with (−7)/(7q-), in association with RPL27A upregulation. This significant downregulation of miR-595 is also apparent when higher risk MDS cases are compared to lower risk cases. The potential clinical importance of these findings requires further validation.

The gene expression signature of anagrelide provides an insight into its mechanism of action and uncovers new regulators of megakaryopoiesis

BACKGROUND

Anagrelide is a cytoreductive agent used to lower platelet counts in essential thrombocythemia. Although the drug has been known to selectively inhibit megakaryopoiesis for many years, the molecular mechanism accounting for this activity is still unclear.

OBJECTIVES AND METHODS

To address this issue we have compared the global gene expression profiles of human hematopoietic cells treated ex-vivo with and without anagrelide while growing under megakaryocyte differentiation conditions, using high-density oligonucleotide microarrays. Gene expression data were validated by the quantitative polymerase chain reaction and mined to identify functional subsets and regulatory pathways.

RESULTS

We identified 328 annotated genes differentially regulated by anagrelide, including many genes associated with platelet functions and with the control of gene transcription. Prominent among the latter was TRIB3, whose expression increased in the presence of anagrelide. Pathway analysis revealed that anagrelide up-regulated genes that are under the control of the transcription factor ATF4, a known TRIB3 inducer. Notably, immunoblot analysis demonstrated that anagrelide induced the phosphorylation of eIF2α, which is an upstream regulator of ATF4, and increased ATF4 protein levels. Furthermore, salubrinal, an inhibitor of eIF2α dephosphorylation, increased the expression of ATF4-regulated genes and blocked megakaryocyte growth.

CONCLUSIONS

These findings link signaling through eIF2α/ATF4 to the anti-megakaryopoietic activity of anagrelide and identify new potential modulators of megakaryopoiesis.

Molecular pathways of early CD105-positive erythroid cells as compared with CD34-positive common precursor cells by flow cytometric cell-sorting and gene expression profiling

Special attention has recently been drawn to the molecular network of different genes that are responsible for the development of erythroid cells. The aim of the present study was to establish in detail the immunophenotype of early erythroid cells and to compare the gene expression profile of freshly isolated early erythroid precursors with that of the CD34-positive (CD34(+)) compartment. Multiparameter flow cytometric analyses of human bone marrow mononuclear cell fractions (n=20) defined three distinct early erythroid stages. The gene expression profile of sorted early erythroid cells was analyzed by Affymetrix array technology. For 4524 genes, a differential regulation was found in CD105-positive erythroid cells as compared with the CD34(+) progenitor compartment (2362 upregulated genes). A highly significant difference was observed in the expression level of genes involved in transcription, heme synthesis, iron and mitochondrial metabolism and transforming growth factor-β signaling. A comparison with recently published data showed over 1000 genes that as yet have not been reported to be upregulated in the early erythroid lineage. The gene expression level within distinct pathways could be illustrated directly by applying the Ingenuity software program. The results of gene expression analyses can be seen at the Gene Expression Omnibus repository.