Long noncoding RNA PCED1B-AS1 promotes erythroid differentiation coordinating with GATA1 and chromatin remodeling

Junwi Zhu,Yunxio Rn,c,Yunyun Hn,Tingting Jin,Ynming Li,Xiuyn Run,Hongzhu Qu,Shngwn Hung,Zhojun Zhng,c,*,Xingdong Fng,c,d,*

aCAS Key Laboratory of Genome Science & Information,Beijing Institute of Genomics,Chinese Academy of Sciences,Beijing,China; bInstitute for Stem Cell and Regeneration,Chinese Academy of Sciences,Beijing,China; cSino-Danish College,University of Chinese Academy of Sciences,Beijing,China; dCollege of Life Sciences,University of Chinese Academy of Sciences,Beijing,China; eGuizhou University,Medical College,Guiyang,Guizhou,China

Abstract

Keywords:Chromatin accessibility,Erythroid differentiation,Long non-coding RNA, PCED1B-AS1

1.INTRODUCTION

Hematopoietic stem cells(HSCs)initiate erythroid differentiation into common myeloid progenitor cells (CMPs),further differentiation into erythroblasts (EBs),and eventually maturation into red blood cells—a complex and dynamic process.1,2HSCs have self-renewal and pluripotency differentiation abilities,and can differentiate into CMPs and lymphoid progenitor cells(CLPs).3CMPs can be further differentiated into EBs,which approach the characteristics of mature red blood cells with morphological and structural changes,including a decrease in cell size and nuclear condensation.4After denucleation,the cells form reticulocytes,which are released into the blood and eventually mature into red blood cells,playing the role of transporting oxygen.3,5

Erythroid differentiation is regulated by a variety of factors,including cytokines,signaling pathways,and transcription factors.3,6In particular,a variety of transcription factors playing a significant role in erythroid differentiation have been discovered,such as GATA1—one of the important transcription factors during erythroid differentiation.7,8The effects ofalncRNA-EC3,alncRNA-EC7,andlincRNA-EPSon erythroid differentiation have been verified by many years of lncRNA research.9–11A number of studies have shown that lncRNA has tissue specificity and plays an important role in stem cell differentiation and other biological process.12–14Chromatin conformation also showed dynamic changes in erythroid differentiation,15,16during the late stages of which the nucleus was concentrated and chromatin might change dramatically.17However,the association between lncRNA in erythroid differentiation regulation and chromatin status remains to be studied.

In order to explore the association of lncRNA and chromatin status in erythroid differentiation,we first screened the lncRNAs 2kb upstream and downstream with GATA1 binding sites based on ChIP-Seq in the ENCODE database,and further identified lncRNA candidates in the Blueprint database containing transcriptome and epigenome data relating to the hemotopoietic system.We found that lncRNAPCED1B-AS1may have an effect on erythroid differentiation.Interestingly,we found that erythroid differentiation regulated byPCED1B-AS1could be associated with GATA1 and chromatin remodeling in erythroblast cells.

2.METHODS

2.1.Cell culture

We cultured TF-1 cells in RPMI 1640 medium containing 10%fetal bovine serum,1%penicillin and streptomycin,and 2μg/ml–8μg/ml granulocyte-macrophage colony stimulating factor.We cultured K562 cells in RPMI 1640 medium containing 10%fetal bovine serum,1% penicillin,and streptomycin,and induced k562 erythroid differentiation with 50μM hemin.We cultured 293T cells in DMEM containing 10% fetal calf serum,1%penicillin and streptomycin.We isolated CD34+ cells from human cord blood using the EasySep Hu CB CD34 Pos Sel Kit II kit and tested its purity by flow cytometry.CD34+cells were first cultured at 105cells/ml for 6 days in SFEM II supplemented with 10%FBS,10ng/ml SCF,1ng/ml IL-3,and 3IU/ml EPO at 37°C in 5% CO2,and then cultured for 8 more days in the above complete medium with the presence of 30%FBS and the absence of IL-3 and SCF.We cultured all of the above cells a constant temperature incubator at 37°C with 5% CO2.The subject has signed an informed consent form,and the sample collection was approved by the ethics committee of Beijing Institute of Genomics(BIG),CAS.

2.2.GATA1 knock down in TF-1 cells

We purchased three GATA1 siRNA sequences from Gemma gene company and diluted them to 20μM with ddH2O(GATA1-siRNA-1:CAUUGCUCAACUGUAUGGATT,GATA1-siRNA-2: CCAAGCUUCGUGGAACUCUTT,GATA1-siRNA-3:GGUACUC AGUGCACCAACUTT).We added 0.8×106TF-1 cells to each well of a six-well plates.We transferred cells of each well with 7.5μl of 20μM of siRNA via lipofectamine 2000 reagent,and cultured the cells in incubator for 48h.We collected 1×106cells and extracted total RNA using TRIzol.We detected knockdown efficiency by RT-PCR after reverse transcription.

2.3.Overexpression of PCED1B-AS1 in K562 cells

We cloned the sequence ofPCED1B-AS1into lentiviral plasmid pHAGE-fEF-1a-IRES-ZsGreen-2.We transfected the overexpressing recombinant construct into 293T cells,together with the packaging plasmids psPAX2 and pMD2.G with Lipofectamine 2000 regents,and collected the cell supernatant medium for 48h and 72h,respectively,as a virus stock solution.After filtration through a 0.45μm filter,we concentrated the virus stock to 50× virus using 5× PEG8000.We added the concentrated virus solution and polybrene at a final concentration of 6μg/ml to a medium in which we cultured K562 cells for infection,and isolated GFP-positive cells by flow cytometry and cultured them to obtain cells stably expressingPCED1B-AS1.

2.4.Quantitative real-time PCR

We extracted total RNA from the cells using the TRIzol Reagent kit,and reverse transcribed total RNA into cDNA using the PrimeScript RT reagent Kit with gDNA Eraser kit.We performed RT-PCR in a CFX96 Real-Time PCR detection system using Maxima SYBR Green/ROX qPCR Master Mix.

2.5.Flow cytometry

We collected 1×106cells,washed them twice with 1ml of 1×PBS,and then resuspended them in 100μl of 1×PBS.We added 0.1μl of phycoerythrin-conjugated anti-CD235a and 0.5μl of allophycocyanin-conjugated anti-CD71 in a dark area,and incubated them in a refrigerator at 4°C for 10min in the dark.We washed the cells twice with 1ml of 1× PBS,and discarded the supernatant to obtain the cell pellet,which we resuspended in 500μl of 1×PBS and detected on a FACSCalibur flow cytometer(BD Biosciences).

2.6.Benzidine staining

We dissolved 10mg benzidine hydrochloride with 1 ml of 0.5 M acetic acid,and added 1μl 30% hydrogen peroxide to 50μl benzidine solution immediately before use.We collected 0.3 to 0.5M cells in 100μl 1×PBS for 2 washes,then resuspended the cells with 10μl of 1×PBS.We added 1μl benzidine solution to the cell suspension and incubated it for 5min at room temperature.We spread the cells on a glass slide,and observed the staining results under a microscope after 30min,counting the proportion of positive cells in three different fields.

2.7.ChIP-seq analysis of K562 cells in ENCODE

We downloaded the GATA1 ChIP-seq data in K562 erythroid cells from ENCODE.18We screened for the GATA1 binding sites within 2kb in the presence of the lncRNAs.The lncRNA is regarded to be potentially regulated by GATA1 in erythroid cells.

2.8.Chromatin immunoprecipitation

Appropriate 50μl magnetic beads (Dynabeads M-280 Sheep anti-Rabbit IgG) were respectively mixed with 2.5μg GATA1 antibody and IgG antibody,and rotated overnight at 4°C.A total of 107K562 cells were cross-linked by 9% formaldehyde solution,and the chromatin was sheared to a fragment of 300 to 500bp in size.The sheared chromatin was divided into two equal parts,in addition to 20μl input and 20μl to identify the ultrasound efficiency,and was added to the beads-GATA1-antibody and beads-IgG-antibody for immunoprecipitation overnight at 4°C.The DNA obtained by immunoprecipitation was purified by Agencourt AMPure XP Kit and detected by qPCR.Primers used are as follows:

Rank1-F: GCCATCTAGTCCCCAGACAA,Rank1-R:CCTCCAAGCTGTCCTTGTTC; Rank2-F: CAGTCCCACTGAGGAAAGGA,Rank2-R:AAGGAAGGGAGGTCTGTGGT;Rank3-F:ACCTGCTGGTTTCAAGTGGA,Rank3-R:CCACTGACGTCTTGATGGAA; Rank4-F: TCCGCTTCCTAAAGCAAGAA,Rank4-R: AGCTTGGGCTATAGGGGAAA.

2.9.Chromatin accessibility analysis

We downloaded the ATAC-seq raw data of HSC,CMP,and EB from the GSE74912 dataset.15We based the data processing on a previous study.15In brief,we trimmed reads by TrimGalore and aligned them to the human genome (GRCh37/hg19) using Bowtie 2.19We identified peaks using MACS2 software.20

2.10.Motif analysis

We selected the peaks near thePCED1B-AS1transcription start site of HSC,CMP,and EB,and performed motif analysis using Homer software with the findmotifsgenome.pl command module.

2.11.The potential regulation factor of PCED1B-AS1

We first identified the regulation factors ofPCED1B-AS1at 100kb using the Cistrom Data Browser Tookit(http://dbtoolkit.cistrome.org/)21,22and then retained the factors present in EB cells with a regulatory potential(RP)score greater than 0.5.We set regulation potential(RP)scores of ＞0.8,0.6–0.8,and ＜0.6 as high,medium,and low,respectively.

2.12.Target gene prediction of PCED1B-AS1 and correlation analysis

We selected GATA1 as a potential regulation factor ofPCED1B-AS1based on the previous step analysis,which we also performed using the Cistrom Data Browser Tookit.We then obtained the top 10 putative targets ofPCED1B-AS1at the bottom of the results page according to public data.We also predicted target genes ofPCED1B-AS1using Starbase 3.0.23We performed correlation analysis using the rcorr function in the Hmisc package and visualized it using the corrplot R package.

2.13.Statistical analysis

We performed each experiment three times or more,and tested the statistical difference between samples for significance using the Kruskal-Wallis test in R.We consideredP≤.05 to be statistically significant.

3.RESULTS

3.1.LncRNA PCED1B-AS1 was identified as a potential regulator in promoting erythroid differentiation

In erythroid differentiation,GATA1 is an essential transcription factor involved in a considerable number of target genes that regulate the proliferation,differentiation,and survival of hematopoietic progenitors.24Therefore,lncRNA regulated by GATA1 is more likely to play a role in erythropoiesis.To screen for potential GATA1-regulated lncRNAs,we first analyzed the ChIP-seq data from K562 erythroid cells in the ENCODE database and identified 107 lncRNAs with GATA1 binding sites within 2kb of their genes.40We extracted the expression data from these 107 lncRNAs in the following three stages—HSCs,CMPs,and EBs,considered to be representative stages in the process of erythropoiesis,and analyzed them for differences in expression.

Significant changes in the expression of lncRNAs might indicate that they have the potential to regulate erythroid differentiation.We utilized DESeq2 to analyze lncRNA expression in three stages of HSCs,CMPs,and EBs.With padj ＜.05 and|log2FoldChange|＞1 as thresholds,we finally identified one significantly differentially expressed lncRNA—PCED1B-AS1—expressed in EBs at a far higher level than those of the other two stages,HSCs and CMPs(Fig.1A).The variant ofPCED1B-AS1is ENST00000552269 that is the only transcript ofPCED1BAS1screened in this study.

Next,to further confirm the potential ofPCED1B-AS1in erythroid differentiation,we observed the endogenous expression pattern ofPCED1B-AS1during erythroid differentiation of cord blood-derived hematopoietic stem and progenitor cells(HSPCs),and from day 7 to day 15,the result showed a continuous upward expression trend ofPCED1B-AS1(Fig.1B),suggesting its potential role in promoting erythroid differentiation.

3.2. PCED1B-AS1 promotes erythroid differentiation in K562 cells

To verify whetherPCED1B-AS1promotes erythroid differentiation,we conducted phenotypic analysis of K562 cells withPCED1B-AS1overexpression (Fig.1C) and control cells that were induced toward erythroid differentiation for 3 days by hemin.We first examined the expression changes in globin genes,finding a slight increase in the expression of ε- and β-globin before and after induction in the cells overexpressingPCED1BAS1compared with the control cells.However,the expression of γ-globin significantly increases after induction compared with the control cells(Fig.1D).We analyzed the expression of cell surface markers CD235a and CD71 by flow cytometry.We found that the proportion of CD235a and CD71 double positive cells in the overexpression group was higher than that in the control cells after induction (Fig.1F),but there was no significant difference.Benzidine staining results also showed the proportion of stained positive cells in total cells to be far higher than that of the control cells(Fig.1G,H).These results demonstrated the role ofPCED1BAS1in promoting erythroid differentiation of K562 cells.

3.3.GATA1-regulated PCED1B-AS1 coordinating with accessible chromatin participates in erythroid differentiation

To further explore the mechanism behindPCED1B-AS1participation in erythroid differentiation,we analyzed the dynamic changes in chromatin accessibility aroundPCED1BAS1loci during erythroid differentiation.Interestingly,we found that the chromatin was tightly condensed at thePCED1B-AS1promoter region in the EB compared with those at the HSCs and CMPs,while the dense structure of the nucleosome was newly opened at the distal region of thePCED1B-AS1in EB,and the chromatin accessibility increased compared to that at other periods(Fig.2A).Moreover,from the public ChIP-seq data and histone modification data of ENCODE,we observed strong GATA1 binding signals and H3K27ac signals at the distal region ofPCED1B-AS1(Fig.2A).We verified the physical binding of GATA1 in that region by ChIP-qPCR(Fig.2B).H3K27ac histone modification at the distal region ofPCED1B-AS1indicated the presence of an enhancer signal.

We then performed motif analysis of chromatin accessibility peaks aroundPCED1B-AS1loci,including the promoter region and the distal region.We found that the transcription factor enriched in HSCs was mainly from the RUNX family (Fig.2C),which play a significant role in hematopoiesis regulation—particularly in the important function of Runx1 in the generation of HSCs.25–27The transcription factor CLOCK,a well known factor regulating circadian rhythms was enriched in CMPs(Fig.2C).28,29Transcription factors enriched at the EB stage belonged mainly to the GATA family,such as GATA1(Fig.2C),suggesting that accessibility of the chromatin was complete in preparation for the entrance of the GATA family TFs to regulate erythroid differentiation in EBs byPCED1B-AS1.

To further clarify which transcription factor may regulatePCED1B-AS1at the distal region,we first used the published data of ChIP-seq(the protein factor and histone marker)and the chromatin accessibility(DNase-Seq and ATAC-seq)of EB cells to predict the potential regulatory factors within 100kb upstream and downstream ofPCED1B-AS1.The regulatory potential score is reflected by the standardized values of 0–1.We observed that GATA1 has the highest regulatory potential forPCED1BAS1(Fig.2D).Moreover,we found that the expression ofPCED1B-AS1also decreased with GATA1 knockdown in the TF-1 cell line (Fig.2E).These results suggested that—when chromatin was accessible—GATA1 regulatedPCED1B-AS1at the distal region ofPCED1B-AS1.

In order to further explore the target genes regulated byPCED1B-AS1,we predicted the potential target genes ofPCED1B-AS1through Cistrom Data Browser and starBase 3.0 database and calculated the correlation analysis based on their expression.We found thatPCED1B-AS1has the strongest correlation withPPT2(Fig.2F).The palmitoyl protein thioesterase-2 (PPT2) gene encodes a lysosomal thioesterase homologous to PPT1.30The correlation betweenPCED1B-AS1and bothHAUS1andACVRL1is also high(Fig.2F).HAUS1is involved mainly in the assembly of mitotic spindles,which facilitates the maintenance of centrosome integrity and completion of cytokinesis.ACVRL1is a type I receptor of TGF-beta family ligands BMP9/GDF2 and BMP10,and is an important regulator of normal vascular development.

FIGURE1.Identification and functional verification of lncRNA PCED1B-AS1.(A) Dynamic expression of lncRNA PCED1B-AS1 during erythroid differentiation in Blueprint database.(B)Relative expression of PCED1B-AS1 during late stage of erythroid differentiation of CD34+cells.The horizontal axis represents the days after EPO induction.(C)The confirmation of PCED1B-AS1 overexpression in K562 cells by RT-PCR.D3 represents the cells induced with hemin for 3 days.(D)The analysis for globin gene expression changes before and after induction of PCED1B-AS1-overexpressed K562 cells compared with that in control cells.HBE,HBG and HBB are β-like genes that respectively represent ε(embryonic)-,γ(fetal)-and β(adult)-globin and are specifically expressed at different developmental stages during erythroid differentiation.(E)The expression changes of CD235a and CD71 was detected by flow cytometry before and after induction.The percentage in the upper right of each figure represents the proportion of CD235a+-CD71+double positive cells.(F)Percentage of CD235a+-CD71+double positive cells among the K562 cells before and after induction.Data are the means±SD of three experimental replicates.N.S.=no significance between samples.(G,H)Benzidine-stained positive rate changes of K562 cells before and after induction with hemin.All relative mRNA expression was normalized to GAPDH.Statistical results were analyzed by student t-test and Kruskal-Wallis test,*P＜.05,**P＜.01,***P＜.001,nsP＞.05,ns:no significance.D0:Day0,D3:Day3.Ctr:Control,OE:overexpression of PCED1B-AS1.

FIGURE2.The association between lncRNA PCED1B-AS1 and chromatin accessibility during erythroid differentiation.(A) Dynamic profile of chromatin accessibility and histone modification of H3K27ac around PCED1B-AS1 loci during erythroid differentiation.(B)Confirmation of physical binding of GATA1 around PCED1B-AS1 loci by ChIP-qPCR in K562 cells.Rank1 and Rank2 are from the red box on the left side in(A)picture,and Rank3 and Rank4 are from the red box on the right side in(A).The bar graphs represent the average of percentage input(mean±SD)from three independent ChIP experiments.(C)Known motif identification of peaks within 20kb of PCED1B-AS1.(D)Regulation factors of PCED1B-AS1 within its 100kb genomic regions in EB.(E)Relative expression of PCED1B-AS1 in TF-1 cells with GATA1 knockdown by siRNA.(F) Correlation between PCED1B-AS1 with its potential target genes.(G) Model of PCED1B-AS1 participates in erythroid differentiation via dynamic chromatin remodeling involving GATA1.

In summary,we observed that enhanced chromatin accessibility in EBs allowed the transcription factor GATA1 and other GATA family members to bind to the distal region,thereby regulating the transcription ofPCED1B-AS1,which in turn affected erythroid differentiation (Fig.2G).

4.DISCUSSION

LncRNAs are abundant,large in number,and participate in many biological processes,including the process of erythropoiesis.Although several studies have identified multiple lncRNAs that function during erythroid differentiation,the function mechanism of lncRNA remains unclear.In our study,we screened antisense lncRNAPCED1B-AS1from public ENCODE databases and Blueprint,and verified its role in promoting erythroid differentiation in the K562 cell model widely used for phenotypic analysis of the role played by unknown regulators in erythroid differentiation.

lncRNA is characterized by specificity at the cellular stage,a lower expression level than that of the coding gene,14and a low conservation between species.31Different lncRNAs play roles at specific cell stages during erythroid differentiation.For example,H19,highly expressed in long-term HSCs (LT-HSCs),regulates the maintenance and self-renewal of LT-HSCs via the Igf2 and Igf1r pathways,32while lncRNA-EC6,highly expressed in red blood cells in the late stage of erythropoiesis,regulates the denucleation of mouse red blood cells through the RAC1/PIP5K signaling pathway.33The expression level ofPCED1B-AS1gradually increases in the late stage of erythroid differentiation,indicating its important role in the late stage of red blood cell differentiation.The recent study ofPCED1B-AS1found that its expression was decreased in patients with active tuberculosis compared with that in healthy people,and that the apoptosis of monocytes was attenuated and autophagy was enhanced.34Because,in the late stage of erythropoiesis,red blood cell islands are formed by macrophages and cells at different stages of terminal differentiation,and the role of macrophages is to help red blood cells denucleate by phagocytosis.35,36Due to the rapidly increased expression ofPCED1B-AS1at day 14 of the culture,we suspect thatPCED1B-AS1may play a role in this process of denucleation.

The chromatin of eukaryotes is a high-level complex structure formed by dense compaction of nucleosomes.When the dense nucleosome structure is destroyed,the chromatin DNA can be approached by cis-regulatory elements such as a promoter,an enhancer,an insulator,a silencer,and a trans-acting factor.In the past,the study of lncRNA associated with erythroid differentiation was concerned mainly with the regulation of lncRNA by direct cis-regulation or trans-regulation with other factors,37and few explore the changes of chromatin accessibility in this process.By inhibiting Fas-mediated cell death,lncNRA Fas-antisense 1(Fas-AS1 or Saf) participates in the process of erythrocyte maturation,and binding sites of GATA1,KLF1,and NF-κB transcription factors were found in its promoter region.Moreover,in the early stage of EB expansion,NF-kB signaling inhibits the expression of Saf,followed by increased expression of Saf in the late stage of RBC maturation,and the concurrent elevation of GATA1 and KLF1 expression,indicating that Saf may be regulated by these transcription factors.38This suggests that sequential actions of different transcription factors at different times on the same loci lead to changes in gene expression,which in turn regulate the differentiation process.In this study,we found that,at different stages of erythroid differentiation,chromatin accessibility changed in the vicinity ofPCED1B-AS1and was accompanied by successive binding of different transcription factors (RUNX/PU.1/GATA1).After knocking down GATA1 in TF-1 erythroid cells,we detected a decrease in the expression ofPCED1B-AS1,indicating that GATA1 may be the factor that affectsPCED1B-AS1,directly or indirectly,and further regulates erythroid differentiation.

We found there is higher correlation betweenPCED1B-AS1with its potential target genes.For example,PPT2,one of candidate of α-Globin poly(C)binding proteins(αCPs)complex involved in posttranscriptional regulation,in addition to maintaining the stability of RNA,it also participates in translation control during erythroid differentiation.39

The link between lncRNA and chromatin has been reported in other research areas.LncRNA controls the chromatin structure and accessibility of genetic information through interaction with chromatin-modifying enzymes and nucleosome-remodeling factors.40,41For instance,lncRNA TARID acts by coupling with GADD45A to direct DNA demethylation to specific loci in cancer cells to regulate gene expression.41,42The relationship between lncRNA and chromatin in the process of erythroid differentiation remains to be revealed,and in our study we have carried out preliminary explorations.

ACKNOWLEDGMENTS

This research was supported by the Strategic Priority Research Program of the Chinese Academy of Sciences (XDA16010602)and the National Natural Science Foundation of China(81670109,81700097,81870097,81700116).