Abstract
Background Circulating human red blood cells (RBCs) consist of mature erythrocytes and immature reticulocytes. Being anucleated, RBCs lack typical transcriptomes, but are known to contain small amounts of diverse long transcripts and microRNAs. However, the exact role and importance of these RNAs is lacking. Shedding of extracellular vesicles (EVs) from the plasma membrane constitutes an integral mechanism of RBC homeostasis, by which RBCs remove unnecessary cytoplasmic content and cell membrane.Results To study this further, we explored the transcriptomes of RBCs and extracellular vesicles (EVs) of RBCs using next-generation sequencing. Furthermore, to understand the dynamics of the RBC transcriptome, we performed single-cell RNA sequencing on RBCs. An analysis of the single-cell transcriptomes revealed that approximately 10% of the cells contained detectable levels of mRNA and fell into three subpopulations based on their transcriptomes. Decrease in the mRNA quantity was observed across the populations. Qualitative changes included the differences in the globin transcripts and changes in the expression of ribosomal genes. A specific short splice form of a long non-coding RNA, Metastasis Associated Lung Adenocarcinoma Transcript 1 (MALAT1), was the most enriched marker in one subpopulation of RBCs, co-expressing with ribosomal structural transcripts. MALAT1 expression was confirmed by qPCR in CD71-enriched reticulocytes, which were also characterized with imaging flow cytometry at single cell level.Conclusions Analysis of the RBC transcriptome shows enrichment of pathways and functional categories required for the maturation of reticulocytes and erythrocyte functions. The RBC transcriptome was detected in their EVs, making these transcripts available for intercellular communication in blood. Our experiments on single cell level revealed that lncRNA MALAT1 is the marker for one of the three RBC populations co-expressing with a group of ribosomal protein transcripts.