Introduction: Myelodysplastic syndromes (MDS) with SF3B1 mutation represent a subtype of hematopoietic stem cell disorders with bone marrow (BM) erythroid dysplasia and ring sideroblasts. This mutation in a splicing factor gene causes multiple alterations of mRNA maturation. Other mutated splicing factors SRSF2 or U2AF1, produce splicing alterations and formation of unscheduled R-loops, which trigger DNA damage response and genetic instability. Mechanistically, mutant SRSF2 impairs the transcription pause release downstream of TSS allowing invasion of relaxed DNA double helix by nascent RNA. Although previous studies suggested that R-loops inappropriately form in SF3B1 mutant cells, investigations are required to elucidate why MDS with SF3B1 mutation are less prone to genetic instability. Methods: A total of 96 patients including 53 MDS with SF3B1 mutation and 43 other low risk MDS was enrolled in this study. RNA-seq of BM mononuclear cells (MNC) was performed in 21 MDS with and 6 MDS w/o SF3B1 mutation. BM CD34+-derived erythroblasts were expanded from 69 MDS and 22 healthy controls to analyze transcriptome, proteome, R-loop landscape, and DNA replication stress. For global proteome analysis, after protein digestion, peptides were fractionated on SCX column and analyzed by mass spectrometry. R-loops were detected by DNA-RNA immunoprecipitation (DRIP-seq) using S9.6 antibody. Immunofluorescence labeling of pRPA32, 53BP1 and γH2AX, DNA combing, BrdU assay and metabolomics were used to investigate DNA replication stress and DNA biosynthesis in human primary erythroblasts or in murine proerythroblastic (proE) G1E-ER4 Crispr-cas9 Sf3b1K700E/+ cell line. expressing GATA1 under the control of estrogen receptor. Results: In SF3B1-mutated BM MNC samples, DeSeq2 analysis of RNA-seq data revealed 1,058 upregulated genes (log2(FC)>0.6, BH-adj P-value<0.05) involved in DNA replication, cell cycle, nucleotide biosynthesis, DNA repair, and erythroid differentiation. Analysis of splicing identified 3,937 differential splicing events (ΔPSI>0.10, BH-adj P-value<0.05) including 1,256 events of intron retention (IR). Among IR events, 384 loss of IR events correlated with altered expression of 293 genes. In SF3B1-mutated erythroblasts, loss of IR events were more frequent in polychromatophilic erythroblasts (polyE: 383) than in proE/early basophilic erythroblasts (eBasoE: 171) suggesting that gene extinction mediated by IR during normal differentiation was altered. Quantitative proteomics revealed 924 significantly deregulated proteins with little overlap with IR transcripts, consistently involved in DNA damage, DNA repair, and nucleotide biosynthesis pathways. DRIP-seq was performed in proE/eBasoE of 8 MDS w/o or with SF3B1 mutation to investigate R-loop formation. Stringent peak calling identified 2,208 shared peaks in wild-type samples compared to 159 shared peaks in SF3B1-mutated samples, suggesting a loss of R-loops. Integrative analysis and DRIP-qPCR demonstrated that R-loop losses overlapped IR loss events, consistent with the notion that intron excision prevents R-loop formation. Compared to SRSF2 or U2AF1-mutated erythroblasts, SF3B1-mutated cells exhibited a significant increase of DNA synthesis by BrdU labeling and of replication fork speed by DNA combing. Immunofluorescence analysis detected pRPA32 foci marking single strand DNA but no evidence of double strand breaks in contrast to SRSF2 mutant cells in which 53BP1 or γH2AX foci were detected. Compared to isogenic cell lines, Sf3b1K700E/+ proE G1E-ER4 cells demonstrated higher proliferation rate and fork speed in proliferation and in estradiol-induced differentiation conditions, and significantly lower amounts of deoxynucleotides (dATP, dCTP, TTP) suggesting persistence of DNA synthesis during differentiation. Targeting of DNA polymerase by aphidicolin or ribonucleotide reductase by low doses of hydroxyurea partially rescued erythroid differentiation. Conclusion: Our data demonstrated a highly specific signature of DNA replication stress at transcriptomic and proteomic level in SF3B1-mutant cells. High replication fork speed reflects the loss of R-loops in these cells, which may reduce the risk of genetic instability at the expense of normal differentiation. Pharmacological targeting of the DNA replication stress could improve erythropoiesis.