Spanish Guidelines for the use of targeted deep sequencing in myelodysplastic syndromes and chronic myelomonocytic leukaemia.

The landscape of medical sequencing has rapidly changed with the evolution of next generation sequencing (NGS). These technologies have contributed to the molecular characterization of the myelodysplastic syndromes (MDS) and chronic myelomonocytic leukaemia (CMML), through the identification of recurrent gene mutations, which are present in >80% of patients. These mutations contribute to a better classification and risk stratification of the patients. Currently, clinical laboratories include NGS genomic analyses in their routine clinical practice, in an effort to personalize the diagnosis, prognosis and treatment of MDS and CMML. NGS technologies have reduced the cost of large-scale sequencing, but there are additional challenges involving the clinical validation of these technologies, as continuous advances are constantly being made. In this context, it is of major importance to standardize the generation, analysis, clinical interpretation and reporting of NGS data. To that end, the Spanish MDS Group (GESMD) has expanded the present set of guidelines, aiming to establish common quality standards for the adequate implementation of NGS and clinical interpretation of the results, hoping that this effort will ultimately contribute to the benefit of patients with myeloid malignancies.

Identification of Recurrent Mutations in the microRNA-Binding Sites of B-Cell Lymphoma-Associated Genes in Follicular Lymphoma.

Follicular lymphoma (FL) is a common indolent B-cell lymphoma that can transform into the more aggressive transformed FL (tFL). However, the molecular process driving this transformation is uncertain. In this work, we aimed to identify microRNA (miRNA)-binding sites recurrently mutated in follicular lymphoma patients, as well as in transformed FL patients. Using whole-genome sequencing data from FL tumors, we discovered 544 mutations located in bioinformatically predicted microRNA-binding sites. We then studied these specific regions using targeted sequencing in a cohort of 55 FL patients, found 16 recurrent mutations, and identified a further 69 variants. After filtering for QC, we identified 21 genes with mutated miRNA-binding sites that were also enriched for B-cell-associated genes by Gene Ontology. Over 40% of mutations identified in these genes were present exclusively in tFL patients. We validated the predicted miRNA-binding sites of five of the genes by luciferase assay and demonstrated that the identified mutations in BCL2 and EZH2 genes impaired the binding efficiency of miR-5008 and miR-144 and regulated the endogenous levels of messenger RNA (mRNA).

The First Dose of Fingolimod Affects Circulating Extracellular Vesicles in Multiple Sclerosis Patients.

Extracellular vesicles (EVs) are membrane-bound particles involved in intercellular communication. They carry proteins, lipids, and nucleotides such as microRNAs (miRNAs) from the secreting cell that can modulate target cells. We and others have previously described the presence of EVs in peripheral blood of multiple sclerosis (MS) patients and postulated them as novel biomarkers. However, their immune function in MS pathogenesis and the effect during the onset of new immunomodulatory therapies on EVs remain elusive. Here, we isolated plasma EVs from fingolimod-treated MS patients in order to assess whether EVs are affected by the first dose of the treatment. We quantified EVs, analyzed their miRNA cargo, and checked their immune regulatory function. Results showed an elevated EV concentration with a dramatic change in their miRNA cargo 5 h after the first dose of fingolimod. Besides, EVs obtained prior to fingolimod treatment showed an increased immune regulatory activity compared to EVs obtained 5 h post-treatment. This work suggests that EVs are implicated in the mechanism of action of immunomodulatory treatments from the initial hours and opens a new avenue to explore a potential use of EVs for early treatment monitoring.

New Concepts in Cancer Biomarkers: Circulating miRNAs in Liquid Biopsies.

The effective and efficient management of cancer patients relies upon early diagnosis and/or the monitoring of treatment, something that is often difficult to achieve using standard tissue biopsy techniques. Biological fluids such as blood hold great possibilities as a source of non-invasive cancer biomarkers that can act as surrogate markers to biopsy-based sampling. The non-invasive nature of these "liquid biopsies" ultimately means that cancer detection may be earlier and that the ability to monitor disease progression and/or treatment response represents a paradigm shift in the treatment of cancer patients. Below, we review one of the most promising classes of circulating cancer biomarkers: microRNAs (miRNAs). In particular, we will consider their history, the controversy surrounding their origin and biology, and, most importantly, the hurdles that remain to be overcome if they are really to become part of future clinical practice.

The circulating transcriptome as a source of non-invasive cancer biomarkers: concepts and controversies of non-coding and coding RNA in body fluids.

The gold standard for cancer diagnosis remains the histological examination of affected tissue, obtained either by surgical excision, or radiologically guided biopsy. Such procedures however are expensive, not without risk to the patient, and require consistent evaluation by expert pathologists. Consequently, the search for non-invasive tools for the diagnosis and management of cancer has led to great interest in the field of circulating nucleic acids in plasma and serum. An additional benefit of blood-based testing is the ability to carry out screening and repeat sampling on patients undergoing therapy, or monitoring disease progression allowing for the development of a personalized approach to cancer patient management. Despite having been discovered over 60 years ago, the clear clinical potential of circulating nucleic acids, with the notable exception of prenatal diagnostic testing, has yet to translate into the clinic. The recent discovery of non-coding (nc) RNA (in particular micro(mi)RNAs) in the blood has provided fresh impetuous for the field. In this review, we discuss the potential of the circulating transcriptome (coding and ncRNA), as novel cancer biomarkers, the controversy surrounding their origin and biology, and most importantly the hurdles that remain to be overcome if they are really to become part of future clinical practice.

Targeted next-generation sequencing and non-coding RNA expression analysis of clear cell papillary renal cell carcinoma suggests distinct pathological mechanisms from other renal tumour subtypes.

Clear cell tubulopapillary renal cell carcinoma (CCPRCC) is a recently described rare renal malignancy that displays characteristic gross, microscopic and immunohistochemical differences from other renal tumour types. However, CCPRCC remains a very poorly understood entity. We therefore sought to elucidate some of the molecular mechanisms involved in this neoplasm by carrying out targeted next-generation sequencing (NGS) to identify associated mutations, and in addition examined the expression of non-coding (nc) RNAs. We identified multiple somatic mutations in CCPRCC cases, including a recurrent [3/14 cases (21%)] non-synonymous T992I mutation in the MET proto-oncogene, a gene associated with epithelial-to-mesenchymal transition (EMT). Using a microarray approach, we found that the expression of mature (n = 1105) and pre-miRNAs (n = 1105), as well as snoRNA and scaRNAs (n = 2214), in CCPRCC cases differed from that of clear cell renal cell carcinoma (CCRCC) or papillary renal cell carcinoma (PRCC) tumours. Surprisingly, and unlike other renal tumour subtypes, we found that all five members of the miR-200 family were over-expressed in CCPRCC cases. As these miRNAs are intimately involved with EMT, we stained CCPRCC cases for E-cadherin, vimentin and ß-catenin and found that the tumour cells of all cases were positive for all three markers, a combination rarely reported in other renal tumours that could have diagnostic implications. Taken together with the mutational analysis, these data suggest that EMT in CCPRCC tumour cells is incomplete or blocked, consistent with the indolent clinical course typical of this malignancy. In summary, as well as describing a novel pathological mechanism in renal carcinomas, this study adds to the mounting evidence that CCPRCC should be formally considered a distinct entity. Microarray data have been deposited in the GEO database [GEO accession number (GSE51554)].

MicroRNAs in Lymphoma: Regulatory Role and Biomarker Potential.

Although it is now evident that microRNAs (miRNAs) play a critical regulatory role in many, if not all, pathological and physiological processes, remarkably they have only formally been recognized for less than fifteen years. These endogenously produced short non-coding RNAs have created a new paradigm of gene control and have utility as both novel biomarkers of cancer and as potential therapeutics. In this review we consider the role of miRNAs in lymphoid biology both under physiological (i.e. lymphopoiesis) and malignant (i.e. lymphomagenesis) conditions. In addition to the functional significance of aberrant miRNA expression in lymphomas we discuss their use as novel biomarkers, both as a in situ tumour biomarker and as a non-invasive surrogate for the tumour by testing miRNAs in the blood of patients. Finally we consider the use of these molecules as potential therapeutic agents for lymphoma (and other cancer) patients and discuss some of the hurdles yet to be overcome in order to translate this potential into clinical practice.

Silencing of ASXL1 impairs the granulomonocytic lineage potential of human CD34⁺ progenitor cells.

The ASXL1 gene encodes a chromatin-binding protein involved in epigenetic regulation in haematopoietic cells. Loss-of-function ASXL1 mutations occur in patients with a range of myeloid malignancies and are associated with adverse outcome. We have used lentiviral-based shRNA technology to investigate the effects of ASXL1 silencing on cell proliferation, apoptosis, myeloid differentiation and global gene expression in human CD34(+) cells differentiated along the myeloid lineage in vitro. ASXL1-deficient cells showed a significant decrease in the generation of CD11b(+) and CD15(+) cells, implicating impaired granulomonocytic differentiation. Furthermore, colony-forming assays showed a significant increase in the number of multipotent mixed lineage colony-forming unit (CFU-GEMM) colonies and a significant decrease in the numbers of granulocyte-macrophage CFU (CFU-GM) and granulocyte CFU (CFU-G) colonies in ASXL1-deficient cells. Our data suggests that ASXL1 knockdown perturbs human granulomonocytic differentiation. Gene expression profiling identified many deregulated genes in the ASXL1-deficient cells differentiated along the granulomonocytic lineage, and pathway analysis showed that the most significantly deregulated pathway was the LXR/RXR activation pathway. ASXL1 may play a key role in recruiting the polycomb repressor complex 2 (PRC2) to specific loci, and we found over-representation of PRC2 targets among the deregulated genes in ASXL1-deficient cells. These findings shed light on the functional role of ASXL1 in human myeloid differentiation.

Mutations in SETBP1 are recurrent in myelodysplastic syndromes and often coexist with cytogenetic markers associated with disease progression.

Whole exome sequencing was performed in a patient with myelodysplastic syndrome before and after progression to acute myeloid leukaemia. Mutations in several genes, including SETBP1, were identified following leukaemic transformation. Screening of 328 patients with myeloid disorders revealed SETBP1 mutations in 14 patients (4·3%), 7 of whom had -7/del(7q) and 3 had i(17)(q10), cytogenetic markers associated with shortened overall survival and increased risk of leukaemic evolution. SETBP1 mutations were frequently acquired at the time of leukaemic evolution, coinciding with increase of leukaemic blasts. These data suggest that SETBP1 mutations may play a role in MDS and chronic myelomonocytic leukaemia disease progression.

Targeted re-sequencing analysis of 25 genes commonly mutated in myeloid disorders in del(5q) myelodysplastic syndromes.

Interstitial deletion of chromosome 5q is the most common chromosomal abnormality in myelodysplastic syndromes. The catalogue of genes involved in the molecular pathogenesis of myelodysplastic syndromes is rapidly expanding and next-generation sequencing technology allows detection of these mutations at great depth. Here we describe the design, validation and application of a targeted next-generation sequencing approach to simultaneously screen 25 genes mutated in myeloid malignancies. We used this method alongside single nucleotide polymorphism-array technology to characterize the mutational and cytogenetic profile of 43 cases of early or advanced del(5q) myelodysplastic syndromes. A total of 29 mutations were detected in our cohort. Overall, 45% of early and 66.7% of advanced cases had at least one mutation. Genes with the highest mutation frequency among advanced cases were TP53 and ASXL1 (25% of patients each). These showed a lower mutation frequency in cases of 5q- syndrome (4.5% and 13.6%, respectively), suggesting a role in disease progression in del(5q) myelodysplastic syndromes. Fifty-two percent of mutations identified were in genes involved in epigenetic regulation (ASXL1, TET2, DNMT3A and JAK2). Six mutations had allele frequencies <20%, likely below the detection limit of traditional sequencing methods. Genomic array data showed that cases of advanced del(5q) myelodysplastic syndrome had a complex background of cytogenetic aberrations, often encompassing genes involved in myeloid disorders. Our study is the first to investigate the molecular pathogenesis of early and advanced del(5q) myelodysplastic syndromes using next-generation sequencing technology on a large panel of genes frequently mutated in myeloid malignancies, further illuminating the molecular landscape of del(5q) myelodysplastic syndromes.

Assessment of Minimal Residual Disease by Next Generation Sequencing in Peripheral Blood as a Complementary Tool for Personalized Transplant Monitoring in Myeloid Neoplasms.

Patients with myeloid neoplasms who relapsed after allogenic hematopoietic stem cell transplant (HSCT) have poor prognosis. Monitoring of chimerism and specific molecular markers as a surrogate measure of relapse is not always helpful; therefore, improved systems to detect early relapse are needed. We hypothesized that the use of next generation sequencing (NGS) could be a suitable approach for personalized follow-up post-HSCT. To validate our hypothesis, we analyzed by NGS, a retrospective set of peripheral blood (PB) DNA samples previously evaluated by high-sensitive quantitative PCR analysis using insertion/deletion polymorphisms (indel-qPCR) chimerism engraftment. Post-HCST allelic burdens assessed by NGS and chimerism status showed a similar time-course pattern. At time of clinical relapse in 8/12 patients, we detected positive NGS-based minimal residual disease (NGS-MRD). Importantly, in 6/8 patients, we were able to detect NGS-MRD at time points collected prior to clinical relapse. We also confirmed the disappearance of post-HCST allelic burden in non-relapsed patients, indicating true clinical specificity. This study highlights the clinical utility of NGS-based post-HCST monitoring in myeloid neoplasia as a complementary specific analysis to high-sensitive engraftment testing. Overall, NGS-MRD testing in PB is widely applicable for the evaluation of patients following HSCT and highly valuable to personalized early treatment intervention when mixed chimerism is detected.

Strategy for identification of a potential inherited leukemia predisposition in a 299 patient's cohort with tumor-only sequencing data.

Myeloid neoplasms (MN) are usually sporadic late-onset cancers; nevertheless, growing evidence suggests that ∼5% of the cases could emerge as a consequence of inherited predisposition. Distinguishing somatic from germline variants is of vital importance, in order to establish an appropriate individualized management and counsel the patients and their relatives. Since many of the genes associated with myeloid neoplasm germline predisposition (MNGP) are also affected in sporadic MN, we intended to design a strategy to identify potentially inherited variants in a tumor only NGS panel in a cohort of 299 patients with a variety of MN. We considered as indicative of potential inherited origin, variants detected in BM sample at a ∼50% VAF classified as pathogenic, likely pathogenic or of unknown significance detected in MNGP-related genes. A total of 104 suspicious variants from 90 patients were filtered-in in tumor samples. Mutational patterns, follow-up data, and sequencing of a range of non-myeloid tissues were used for narrowing down the list of suspicious variants, and ultimately discriminate their nature. Our data supports the importance of considering variants found upon tumor-only sequencing as potentially of germline origin, and we offer a pipeline to define the nature of the variants.

Assessment of the clinical utility of four NGS panels in myeloid malignancies. Suggestions for NGS panel choice or design.

The diagnosis of myeloid neoplasms (MN) has significantly evolved through the last few decades. Next Generation Sequencing (NGS) is gradually becoming an essential tool to help clinicians with disease management. To this end, most specialized genetic laboratories have implemented NGS panels targeting a number of different genes relevant to MN. The aim of the present study is to evaluate the performance of four different targeted NGS gene panels based on their technical features and clinical utility. A total of 32 patient bone marrow samples were accrued and sequenced with 3 commercially available panels and 1 custom panel. Variants were classified by two geneticists based on their clinical relevance in MN. There was a difference in panel's depth of coverage. We found 11 discordant clinically relevant variants between panels, with a trend to miss long insertions. Our data show that there is a high risk of finding different mutations depending on the panel of choice, due both to the panel design and the data analysis method. Of note, CEBPA, CALR and FLT3 genes, remains challenging the use of NGS for diagnosis of MN in compliance with current guidelines. Therefore, conventional molecular testing might need to be kept in place for the correct diagnosis of MN for now.

The Circulating Transcriptome as a Source of Biomarkers for Melanoma.

The circulating transcriptome is a valuable source of cancer biomarkers, which, with the exception of microRNAs (miRNAs), remains relatively unexplored. To elucidate which RNAs are present in plasma from melanoma patients and which could be used to distinguish cancer patients from healthy individuals, we used next generation sequencing (NGS), and validation was carried out by qPCR and/or ddPCR. We identified 442 different microRNAs in samples, eleven of which were differentially expressed (p < 0.05). Levels of miR-134-5p and miR-320a-3p were significantly down-regulated (p < 0.001) in melanoma samples (n = 96) compared to healthy controls (n = 28). Differentially expressed protein-encoding mRNA 5'-fragments were enriched for the angiopoietin, p21-activated kinase (PAK), and EIF2 pathways. Levels of ATM1, AMFR, SOS1, and CD109 gene fragments were up-regulated (p < 0.001) in melanoma samples (n = 144) compared to healthy controls (n = 41) (AUC = 0.825). Over 40% of mapped reads were YRNAs, a class of non-coding RNAs that to date has been little explored. Expression levels of RNY3P1, RNY4P1, and RNY4P25 were significantly higher in patients with stage 0 disease than either healthy controls or more advanced stage disease (p < 0.001). In conclusion, we have identified a number of novel RNA biomarkers, which, most importantly, we validated in multi-center retrospective and prospective cohorts, suggesting potential diagnostic use of these RNA species.

Noncoding RNA Expression and Targeted Next-Generation Sequencing Distinguish Tubulocystic Renal Cell Carcinoma (TC-RCC) from Other Renal Neoplasms.

Tubulocystic renal cell carcinoma (TC-RCC) is a rare recently described renal neoplasm characterized by gross, microscopic, and immunohistochemical differences from other renal tumor types and was recently classified as a distinct entity. However, this distinction remains controversial particularly because some genetic studies suggest a close relationship with papillary RCC (PRCC). The molecular basis of this disease remains largely unexplored. We therefore performed noncoding (nc) RNA/miRNA expression analysis and targeted next-generation sequencing mutational profiling on 13 TC-RCC cases (11 pure, two mixed TC-RCC/PRCC) and compared with other renal neoplasms. The expression profile of miRNAs and other ncRNAs in TC-RCC was distinct and validated 10 differentially expressed miRNAs by quantitative RT-PCR, including miR-155 and miR-34a, that were significantly down-regulated compared with PRCC cases (n = 22). With the use of targeted next-generation sequencing we identified mutations in 14 different genes, most frequently (>60% of TC-RCC cases) in ABL1 and PDFGRA genes. These mutations were present in <5% of clear cell RCC, PRCC, or chromophobe RCC cases (n > 600) of The Cancer Genome Atlas database. In summary, this study is by far the largest molecular study of TC-RCC cases and the first to investigate either ncRNA expression or their genomic profile. These results add molecular evidence that TC-RCC is indeed a distinct entity from PRCC and other renal neoplasms.

Mutational profiling can identify laryngeal dysplasia at risk of progression to invasive carcinoma.

Early diagnosis of laryngeal squamous cell carcinoma (LSCC) at the stage of dysplasia could greatly improve the outcome of affected patients. For the first time we compared the mutational landscape of non-progressing dysplasia (NPD; n = 42) with progressing dysplasia (PD; n = 24), along with patient-matched LSCC biopsies; a total of 90 samples. Using targeted next-generation sequencing identified non-synonymous mutations in six genes (PIK3CA, FGFR3, TP53, JAK3, MET, FBXW7), and mutations were validated by Sanger sequencing and/or qPCR. Analysis was extended in silico to 530 head and neck (HNSCC) cases using TCGA data. Mutations in PIK3CA and FGFR3 were detected in PD and LSCC cases, as well as other HNSCC cases, but absent in NPD cases. In contrast, mutations in JAK3, MET and FBXW7 were found in NPD cases but not PD, LSCC or other HNSCC cases. TP53 was the most frequently mutated gene in both PD and NPD cases. With the exception of R248W, mutations were mutually exclusive. Moreover, five of seven PD mutations were located in motif H2 of p53, whereas none of the NPD mutations were. In summary, we propose that the mutational profile of laryngeal dysplasia has utility for the early detection of patients at risk of progression.

Spatial intratumoural heterogeneity in the expression of GIT1 is associated with poor prognostic outcome in oestrogen receptor positive breast cancer patients with synchronous lymph node metastases.

Background: The outcome for oestrogen receptor positive (ER+) breast cancer patients has improved greatly in recent years largely due to targeted therapy. However, the presence of involved multiple synchronous lymph nodes remains associated with a poor outcome. Consequently, these patients would benefit from the identification of new prognostic biomarkers and therapeutic targets. The expression of G-protein-coupled receptor kinase-interacting protein 1 (GIT1) has recently been shown to be an indicator of advanced stage breast cancer. Therefore, we investigated its expression and prognostic value of GIT1 in a cohort of 140 ER+ breast cancer with synchronous lymph node involvement. Methods: Immunohistochemistry was employed to assess GIT1 expression in a tissue microarray (TMA) containing duplicate non-adjacent cores with matched primary tumour and lymph node tissue (n=140). GIT1 expression in tumour cells was scored and statistical correlation analyses were carried out. Results: The results revealed a sub-group of patients that displayed discordant expression of GIT1 between the primary tumour and the lymph nodes (i.e. spatial intratumoural heterogeneity). We observed that loss of GIT1 expression in the metastasis was associated with a shorter time to recurrence, poorer overall survival, and a shorter median survival time. Moreover, multivariate analysis demonstrated that GIT1 expression was an independent prognostic indicator. Conclusions: GIT1 expression enabled the identification of a sub-class of ER+ patients with lymph node metastasis that have a particularly poor prognostic outcome. We propose that this biomarker could be used to further stratify ER+ breast cancer patients with synchronous lymph node involvement and therefore facilitate adjuvant therapy decision making.

Why some tumours trigger neovascularisation and others don't: the story thus far.

BACKGROUND: Angiogenesis is not essential for tumours to develop and expand, as cancer can also grow in a non-angiogenic fashion, but why this type of growth occurs is unknown. Surprisingly, our data from mRNA transcription profiling did not show any differences in the classical angiogenic pathways, but differences were observed in mitochondrial metabolic pathways, suggesting a key role for metabolic reprogramming. We then validated these results with mRNA profiling by investigating differential protein expression via immunohistochemistry in angiogenic and non-angiogenic non-small cell lung cancers (NSCLCs). METHODS: Immunohistochemical staining for 35 angiogenesis- and hypoxia-related biomarkers were performed on a collection of 194 angiogenic and 73 non-angiogenic NSCLCs arranged on tissue microarrays. Sequencing of P53 was performed with frozen tissue samples of NSCLC. RESULTS: The non-angiogenic tumours were distinguished from the angiogenic ones by having higher levels of proteins associated with ephrin pathways, mitochondria, cell biogenesis, and hypoxia-inducible factor 1 (HIF1) regulation by oxygen and transcription of HIF-controlled genes but lower levels of proteins involved in the stroma, cell-cell signaling and adhesion, integrins, and Delta-Notch and epidermal growth factor (EGF)-related signaling. However, proteins classically associated with angiogenesis were present in both types of tumours at very comparable levels. Cytoplasmic expression of P53 was strongly associated with non-angiogenic tumours. A pilot investigation showed that P53 mutations were observed in 32.0% of angiogenic cases but in 71.4% of non-angiogenic tumours. CONCLUSIONS: Our observations thus far indicate that both angiogenic and non-angiogenic tumours experience hypoxia/HIF and vascular endothelial growth factor (VEGF) pathway protein expression in a comparable fashion. However, angiogenesis does not ensue in the non-angiogenic tumours. Surprisingly, metabolic reprogramming seems to distinguish these two types of neoplastic growth. On the basis of these results, we raise the hypothesis that in some, but not in all cases, initial tissue remodeling and/or inflammation could be one of the secondary steps necessary to trigger angiogenesis. In the non-angiogenic tumours, in which neovascularisation fails to occur, HIF pathway activation could be the driving force toward metabolic reprogramming.

The transporter ABCB7 is a mediator of the phenotype of acquired refractory anemia with ring sideroblasts.

Refractory anemia with ring sideroblasts (RARS) is characterized by mitochondrial ferritin (FTMT) accumulation and markedly suppressed expression of the iron transporter ABCB7. To test the hypothesis that ABCB7 is a key mediator of ineffective erythropoiesis of RARS, we modulated its expression in hematopoietic cells. ABCB7 up and downregulation did not influence growth and survival of K562 cells. In normal bone marrow, ABCB7 downregulation reduced erythroid differentiation, growth and colony formation, and resulted in a gene expression pattern similar to that observed in intermediate RARS erythroblasts, and in the accumulation of FTMT. Importantly, forced ABCB7 expression restored erythroid colony growth and decreased FTMT expression level in RARS CD34+ marrow cells. Mutations in the SF3B1 gene, a core component of the RNA splicing machinery, were recently identified in a high proportion of patients with RARS and 11 of the 13 RARS patients in this study carried this mutation. Interestingly, ABCB7 exon usage differed between normal bone marrow and RARS, as well as within the RARS cohort. In addition, SF3B1 silencing resulted in downregulation of ABCB7 in K562 cells undergoing erythroid differentiation. Our findings support that ABCB7 is implicated in the phenotype of acquired RARS and suggest a relation between SF3B1 mutations and ABCB7 downregulation.