Peripheral blood circular RNA circ-0008102 may serve as a novel clinical biomarker in beta-thalassemia patients

Circular RNA circ-0008102 has previously been found dysregulated in β-thalassemia (β-thal) in circRNAs microarray (GSE196682 and GSE241141). Our study is aimed at identifying whether circ-0008102 could be a novel biomarker in β-thal. The peripheral blood of pediatric β-thal patients with (n = 39) or without (n = 20) blood transfusion and healthy controls (n = 30) was selected. qRT-PCR, ROC curve analysis, Spearman correlation analysis, and FISH were used to analyze clinical value of circ-0008102. qRT-PCR confirmed that circ-0008102 expression in pediatric β-thal patients without blood transfusion was significantly higher. ROC curves analysis showed that the AUC of circ-0008102 for differentiating patients without blood transfusion from patients with blood transfusion and healthy controls with an AUC of 0.733 and 0.711. Furthermore, circ-0008102 expression was positively correlated with the levels of RBC, HbF, β-globin, and γ-globin mRNA, but was negatively corrected with the levels of HbA and Cr. circ-0008102 was mainly located in the cytoplasm. circ-0008102 could induce the activation of γ-globin and negatively regulate the expression of the five highest-ranking candidate miRNAs (miR-372-3p, miR-329-5p, miR-198, miR-152-5p, and miR-627-3p) in K562 cells. Conclusion: We demonstrate that peripheral blood upregulated circ-0008102 may serve as a novel clinical biomarker for pediatric β-thal without blood transfusion. What is known: • CircRNAs are known to be involved in various human diseases, and several circRNAs are regarded as a class of promising blood-based biomarkers for detection of β-thal. • CircRNAs exert biological functions by epigenetic modification and gene expression regulation, and dysregulated circRNAs in β-thal might be involved in the induction of HbF in β-thal. What is new: • Peripheral blood circ-0008102 maybe serve as a novel clinical biomarker for detection of pediatric β-thal without blood transfusion. • Circ-0008102 participates in the pathogenesis of β-thal through regulating γ-globin expression, and negatively regulates the expression of miR-372-3p, miR-329-5p, miR-198, miR-152-5p and miR-627-3p. Supplementary Information The online version contains supplementary material available at 10.1007/s00431-023-05398-y.


Background
β-Thalassemia (β-thal) is an autosomal recessive hereditary anemia characterized by reduced or absent β-globin chain synthesis, which is most widespread among Mediterranean countries, Central Asia, Middle East, India, South America, countries alongside the North coast of Africa, and Southern China.It has been assessed that approximately 68,000 children are born with various thalassemia syndromes each year worldwide [1].According to clinical manifestations, there are three main forms of β-thal, including β-thal major (TM), β-thal intermedia (TI), and β-thal carrier (TC, also called thalassemia minor).Moreover, β-thal can be clinically classified as transfusion-dependent thalassemia (TDT) and non-transfusion-dependent thalassemia (NTDT) according to the severity of the phenotype, caused by a wide spectrum of mutations in a homozygous or compound heterozygous state.The β-globin gene, located on chromosome 11p15.5, is responsible for the coding of β-globin chain.Patients with TM or TI are chiefly homozygotes or compound heterozygotes for β 0 or β + mutation, and patients with TC are heterozygotes for β 0 or β + mutation.Although no treatment is needed for individuals who are diagnosed as TC, patients with TM are mainly treated by long-term blood transfusion and its high cost has brought a heavy economic burden to the family and society [2].Currently, bone marrow or stem cell transplantation is a curable treatment for patients with TM, but it is difficult to popularize clinically due to the limitation of matching and marrow source [3].Thus, an increased knowledge about the pathophysiology of β-thal may contribute to identify novel treatment methods for this disease.
As research progresses, more and more biomarkers associated with β-thal have been identified, including traditional proteins (alpha-hemoglobin pool, transferrin receptor 1, and heat shock protein 70) and non-coding RNAs (ncRNAs) [4,5].At present, studies on these biomarkers are still at the stage of expression analysis and the relevant functional mechanisms have not yet been elucidated.As a specific type of ncRNAs, circular RNAs (circRNAs) are commonly regarded to be derived from precursor mRNAs, with a covalently closed loop structure and non-poly-adenylated tail.Because of this unique structure, circRNAs are stable in the various body fluids and are not easily degraded by Ribonuclease R (RNase R).In recent years, with the development of highthroughput sequencing analysis and bioinformatic methods, thousands of circRNAs have been identified [6].circRNAs are evolutionarily conserved, exhibit cell and tissue-type specific expression across species, and have shown great potential in diagnosis, therapy, and prognosis for a variety of human diseases, including hematological diseases [7,8].For example, circ-MYBL2, a circRNA from MYBL2 gene, is reported to facilitate the progression of acute myeloid leukemia by regulating fms-related receptor tyrosine kinase 3 (FLT3) translation through recruiting polypyrimidine tract binding protein 1 (PTBP1) [9].circ-BUB1B_544aa, generated from BUB1 mitotic checkpoint serine/threonine kinase B (BUB1B) gene, has the potential to aggravate multiple myeloma malignancy through evoking chromosomal instability [10].Besides, in myelodysplastic syndrome, peripheral blood circ-100352, circ-104056, and circ-102817 are valuable for diagnosis and prognosis and are linked with some important cancer-related functions and pathways [11].Nowadays, increasing study has revealed that large amounts of circRNAs are abnormally expressed in β-thal and play important roles in epigenetic modification and gene expression regulation [12].However, the roles and clinical value of many circRNAs in β-thal still remain to be explored.
Here, circ-0008102, also known as circ-100654, is derived from ligand-dependent nuclear receptor corepressor (LCOR) gene and is largely functionally unknown.A previous expression profile of circRNAs by microarray from the GSE196682 dataset has revealed that circ-0008102 expression is upregulated in the peripheral blood from β-thal carriers with high HbF.However, the circRNA microarray data (GSE241141) from our lab have shown that circ-0008102 level was downregulated in β-thal patients.These results suggested that circ-0008102 is dysregulated in β-thal.As a new circRNA, the role and mechanism of circ-0008102 in β-thal remain unknown.Our pre-experiments revealed that circ-0008102 expression in peripheral blood of pediatric β-thal patients was associated with HbF level.Therefore, whether increased circ-0008102 expression is involved in the development of β-thal is worth investigating.The present study is aimed at confirming the expression levels of circ-0008102 in large cohort of pediatric β-thal patients and healthy controls and analyzing its relationship with clinical characteristics.The receiver operating characteristic (ROC) curve analysis was performed to evaluate the diagnostic value of circ-0008102 for pediatric β-thal.Moreover, based on bioinformatics analysis, the circ-0008102/miRNA/ mRNA interaction network was established to help us better understand its biological roles and underlying mechanisms in pediatric β-thal.

Study participants
Fifty-nine pediatric β-thal patients (age 4-13 years, average age 8.12 ± 2.12 years) (20 cases without any blood transfusion for at least the past month prior to sample collection and 39 cases with blood transfusion history before the past month prior to sample collection) were prospective collected in Fujian Province and Guangxi Province, China, between January 2019 and December 2020.Thirty age-and gender-matched healthy controls (age 4-11 years, average age 7.27 ± 2.08 years) were obtained from Fujian Maternity and Child Health Hospital College of Clinical Medicine for Obstetrics and Gynecology and Pediatrics, Fujian Medical University.Including criteria: (1) β-thal patients have anemia symptoms and diagnosed with β-thal by genetic diagnosis; (2) healthy controls have no symptoms of anemia and no pathogenic mutation in the β-globin gene by genetic diagnosis.Exclusion criteria: (1) patients with other chronic and congenital diseases; (2) patients with acute and chronic lung infection; (3) patients with abnormal blood coagulation; (4) patients were under treatment with hydroxyurea; (5) patients with the three common deletional mutations (-SEA, -α3.7, and -α4.2) and three non-deletional mutations (αQSα, αCSα, and αWSα).All subjects had no genetic relationship.This study was approved by the Ethics Review Committee of Fujian Maternity and Child Health Hospital College of Clinical Medicine for Obstetrics and Gynecology and Pediatrics, Fujian Medical University (no: 073,2019) and was conducted in conformity to the Declaration of Helsinki.Written informed consent was obtained from all participants or their parents following a detailed description of the purpose of the study.

Peripheral blood samples
The peripheral blood (5 ml each) was collected into PAXgene Blood RNA Tubes (Qiagen, Hilden, Germany) and was stored immediately at − 80 °C for further RNA extraction.After mixing by inversion of the tubes, approximately 2 ml peripheral blood samples anticoagulated with EDTA-K2 were collected for further hematological analysis.Hematological parameters, including red blood cells (RBC), hemoglobin (Hb), mean corpuscular volume (MCV), and mean corpuscular hemoglobin (MCH), were measured using a Sysmex XN-3000 automatic hematology analyzer (Sysmex Corporation; Shanghai, China).The hemoglobin components and levels, including hemoglobin A (HbA), hemoglobin A 2 (HbA 2 ), and fetal hemoglobin (HbF), were analyzed using an automated capillary S2 electrophoresis system version 6.2 (Sebia, Paris, France).

RNA extraction and reverse transcription
Total RNA from the peripheral blood samples was isolated using a PAXgene Blood RNA Kit (Qiagen) according to the manufacturer's protocol.The purity of the RNA samples was assessed using a Bio Photometer MULTISKAN GO, and the OD260/280 ratio at 1.9 to 2.0 indicated good quality.The RNA integrity was determined by 1% formaldehyde-denatured gel electrophoresis.After calculating the concentration, RNA was frozen in − 80 °C refrigerator.RNA was reverse transcribed into cDNA using a SuperScript™ III Reverse Transcriptase (Invitrogen, Carlsbad, CA, USA) with random primers, according to the manufacturer's protocol.The cDNA was preserved in − 20 °C for further use.

Quantitative real-time PCR (qRT-PCR)
The cDNA of miRNAs was synthesized with a Mir-XTM First Strand Synthesis kit (Takara Bio, Inc.,), according to the manufacturer's instructions.The StepOnePlus™ Real-Time PCR System (Applied Biosystems, Foster City, CA, USA) was adopted to perform qRT-PCR with a PrimeScript™ RT-PCR kit and Mir-X miRNA qRT-PCR TB Green® Kit (Takara Bio, Inc.,), according to the manufacturer's instruction.The primers of circ-0008102, LCOR, β-globin, γ-globin, miR-372-3p, miR-329-5p, miR-198, miR-152-5p, and miR-627-3p were shown in Supplementary Table S1.The β-actin and U6 were used as the internal control for quantitative analysis.qRT-PCR was performed in triplicate for both target and internal control genes, and the negative control (non-cDNA) was run with every experimental plate to assess specificity and to rule out contamination.The relative fold-change of each target gene with respect to β-actin and U6 was calculated by the 2 −ΔΔCt method.

RNAse R treatment
The RNAse R reagent (Epicentre Biotechnologies, Madison, WI, USA) was used to digest linear RNAs.Owing to circular structure, circRNAs could avoid RNase degradation.The total RNA was divided into 2 parts (2 μg each): one was for the digestion of 3 U/μg RNase R at 37 °C for 20 min; another was treated with equal amount of RNAse-free water under the same conditions.Subsequently, the RNA was reversetranscribed into cDNA and qRT-PCR helped to analyze the detected LCOR mRNA and circ-0008102 following the above experimental methods.

Sanger sequencing
Sanger sequencing was used to validate the back-splice junction sequences of circ-0008102.Briefly, the partial sequence fragments of circ-0008102 that contained the back-splice junction were amplified from cDNA using a PrimeSTAR® GXL Premix (Takara Bio, Inc.), according to the manufacturer's recommendation.The amplification primers were as follows: sense, 5′-CCT GCC GAA AGC ATC TCC AGT-3′ and antisense, 5′-TCG CTG CAT GAT CTC ACG G-3′.The amplification primers were also used as primers for sequencing.Sanger sequencing was performed on an ABI 3730XL DNA analyzer (Applied Biosystems) and was analyzed by using a PC ChromasPro software version 2.1.5(Applied Biosystems).

Cell culture and cell transfection
Human K562 cells were purchased from the Cell Center of Shanghai Institutes for Biological Sciences (Shanghai, China).K562 cells were cultured in RPMI-1640 medium (Gibco, Grand Island, NY, USA), supplemented with 10% FBS (Gibco) and 1% penicillin-streptomycin (Sigma-Aldrich, Louis, MO, USA).All cells were maintained at 37 °C in a humidified atmosphere containing 5% CO 2 .

Fluorescence in situ hybridization (FISH)
The FISH kits of circ-0008102 and negative control were designed and purchased from GenePharma (Shanghai, China).Approximately 6 × 10 4 K562 cells were seeded on cell climbing in 24-well plates overnight.Then, the cells were washed with PBS and fixed in 4% paraformaldehyde for 15 min at 37 °C and followed by permeabilized them for 15 min in 0.1% buffer A at 37 °C.After being washed with PBS, 2 × buffer C was added and maintained for 30 min at 37 °C, and then buffer E was added and rehydrated for 30 min at 73 °C.Subsequently, 1 μL of biotin-labeled circ-0008102-probe (5′-AGT GTC CCC GTT CCC TTG AGT ACT G-3′) or negative control-probe (5′-TGC TTT GCA CGG TAA CGC CTG TTT T-3′), 1 μL of SA-Cy3, and 8 μL of PBS were added into the cells and incubated at 37 °C overnight by avoiding light.The next day, the cells were sequentially washed with 0.1% buffer F for 10 min at 37 °C, 2 × buffer C for 10 min (repeat 3 times) at 60 °C, and 2 × buffer C for 10 min (repeat 3 times) at 37 °C.Finally, the nucleus of cells was dyed in 4′,6′-diamidino-2-phenylindole (DAPI), and the subcellular localization of circ-0008102 was detected using a Leica TCS SP8 CARS Confocal Microscope (Leica, Wetzlar, Germany).

Statistical analysis
GraphPad Prism software (version 9.0; GraphPad Software Inc., San Diego, CA, USA) and SPSS software (version 25.0; SPSS, Chicago, IL, USA) were used for statistical analyses.Measurement data were expressed as the mean ± standard deviation (SD), and count data was expressed as case (%).The normality of the sample distributions was calculated by the Kolmogorov-Smirnoff test.The Mann-Whitney U test or independent t-test was used to assess the differences between the pediatric β-thal patients and healthy controls and the differences between pediatric β-thal patients with blood transfusion and pediatric β-thal patients without blood transfusion.The correlations between the expression level of circ-0008102 and hematological parameters, biochemical indicators, and β-globin and γ-globin mRNA expression were evaluated by Spearman correlation analysis.The ROC curve analysis and the area under curve (AUC) were carried on to detect circ-0008102 as biomarker for differentiating pediatric β-thal with or without blood transfusion.P values less than 0.05 were considered statistically significant.

Comparison of baseline characteristics in pediatric β-thal patients and healthy controls
The genotypes of pediatric β-thal patients were shown in Supplementary Table S2, and a total of 19 genotypes were found in the pediatric β-thal patients (n = 59), and the pathogenic mutations of β-globin gene were mainly concentrated in the βCD41-42(-TCTT) (HBB:c.There was no statistical difference of age and sex distribution between pediatric β-thal patients and healthy controls.As expected, compared with healthy controls, pediatric β-thal patients displayed significantly decreased  1).But the levels of HbA 2 (P = 0.008) and HbF (P < 0.001) in pediatric β-thal patients were significantly higher than that in healthy controls.In addition, pediatric β-thal patients showed abnormal biochemical indicators, including TBIL, DBIL, ALT, AST, and SF, and the levels of these parameters were markedly enhanced compared with that in healthy controls (P < 0.001); however, the levels of Cr and TP in pediatric β-thal patients were decreased compared with healthy controls (P < 0.01).There was no statistical difference in MCV, BUN, UA, ALB, and GGT between pediatric β-thal patients and healthy controls.These data confirmed that pediatric β-thal patients had phenotypes of anemia, abnormal liver function, and iron metabolism.Moreover, based on blood transfusion or not, the pediatric β-thal patients were further categorized into two groups: without blood transfusion group (n = 20) and blood transfusion group (n = 39).The hematological parameters and biochemical indicators of the two groups were compared and shown in Supplementary Table S3.Compared with that in patients without blood transfusion, patients with blood transfusion showed higher levels of Hb, MCV, MCH, and HbA and lower level of HbF (P < 0.05).However, no significant difference was observed in age, sex distribution, RBC, HbA 2 , BUN, Cr, UA, TP, ALB, TBIL, DBIL, ALT, AST, GGT, and SF levels in the two groups (P > 0.05).These data demonstrated that blood transfusion could improve anemia indicators and was an effective treatment option for β-thal.

Identification and characterization of circ-0008102
Based on the annotation of circBase, circ-0008102 with a length of 515 nt was produced by back-splicing from exons 5, 6, and 7 of the LCOR gene locus (NM_032440), which was located at human chromosome 10: 98,703,869-98711953 position (strand +, GRCh37/hg19) (Fig. 1A).The junction sequence of circ-0008102 was validated by Sanger sequencing (Fig. 1B).Subsequently, to verify that circ-0008102 was generated from the back-splicing instead of trans-splicing or genomic rearrangements, the divergent and convergent primers were designed to amplify circ-0008102 and LCOR mRNA, respectively.RT-PCR revealed that both cDNA and gDNA could amplify LCOR mRNA using the convergent primers, whereas only cDNA could amplify circ-0008102 using the divergent primers (Fig. 1C).Additionally, qRT-PCR showed that circ-0008102 was more stable than LCOR mRNA in K562 cells under treatment with RNAse R (Fig. 1D, P < 0.01).To summarize, circ-0008102 was confirmed to be a circular RNA.

Validation of upregulated circ-0008102 expression in pediatric β-thal patients without blood transfusion
We explored the expression levels of circ-0008102 in peripheral blood of large cohort of subjects, and results showed that there was no statistical difference of circ-0008102 expression between pediatric β-thal patients and healthy controls (Fig. 2A, P > 0.05).Surprisingly, the expression levels of circ-0008102 in peripheral blood of pediatric β-thal patients without blood transfusion were significantly higher than that in pediatric β-thal patients with blood transfusion and healthy controls (Fig. 2B, P < 0.05).These data demonstrated that circ-0008102 expression was upregulated in pediatric β-thal patients without blood transfusion.

Evaluation of circ-0008102 in peripheral blood as a novel diagnostic biomarker for pediatric β-thal patients without blood transfusion
We then evaluated the diagnostic ability of circ-0008102 to distinguish pediatric β-thal patients from healthy controls through ROC curves analysis.Figure 3A showed that the AUC value was 0.552, and 95% confidence interval (CI) was 0.424-0.679,and the cut-off value was 1.11, with sensitivity was 0.491 and specificity was 0.690 (Fig. 3A, P = 0.441), hinting that circ-0008102 had low clinical diagnostic value for the diagnosis of pediatric β-thal.As shown in Fig. 3B, peripheral blood circ-0008102 expression discriminated pediatric β-thal patients without blood transfusion from pediatric β-thal patients with blood transfusion with an AUC of 0.733 (95% CI 0.590-0.877,sensitivity 0.722, and specificity 0.686) (P = 0.006).Similarly, the AUC of circ-0008102 for differentiating pediatric β-thal patients without blood transfusion from healthy controls was 0.711 (95% CI 0.554-0.868),with 0.722 sensitivity and 0.690 specificity (Fig. 3C, P = 0.016).These data demonstrated peripheral blood circ-0008102 might be a novel clinical biomarker in pediatric β-thal without blood transfusion.

Analysis of correlation between circ-0008102 expression and clinical characteristics in pediatric β-thal patients without blood transfusion
As shown in Supplementary Table S4, the association between circ-0008102 expression and hematological parameters in pediatric β-thal patients without blood transfusion showed that the expression level of circ-0008102 was negatively correlated with the level of HbA (r = − 0.494, P = 0.037), but positively correlated with the levels of RBC (r = 0.521, P = 0.026) and HbF level (r = 0.480, P = 0.044).However, no association was found between circ-0008102 expression and Hb, MCV, MCH, and HbA 2 (all P > 0.05) .Similarly, Spearman correlation analysis revealed that the expression level of circ-0008102 was negatively correlated with the level of Cr (r = − 0.726, P = 0.018) in pediatric β-thal patients without blood transfusion.However, no association was found in BUN, UA, TP, ALB, TBIL, DBIL, ALT, AST, GGT, and SF (all P > 0.05) .Together, these data suggested that upregulated circ-0008102 expression was Fig. 2 The upregulated expression of circ-0008102 in pediatric β-thal patients without blood transfusion.A qRT-PCR analysis of the expression levels of circ-0008102 (normalized to β-actin) in peripheral blood of pediatric healthy controls and pediatric β-thal patients.
B The circ-0008102 expression was upregulated in pediatric β-thal patients without blood transfusion relative to pediatric β-thal patients with blood transfusion and healthy controls.Data were presented as the means ± SD. *P < 0.05, ns P > 0.05 associated with the pathogenesis of pediatric β-thal without blood transfusion.

Analysis of β-globin and γ-globin expression in pediatric β-thal patients without blood transfusion and its association with circ-0008102 expression
Compared with healthy controls, the mRNA expression levels of β-globin in the peripheral blood from pediatric β-thal patients without blood transfusion were significantly downregulated with 0.411 fold changes (Fig. 4A, P < 0.05); however, the mRNA expression levels of γ-globin were significantly upregulated with 1087.09fold changes (Fig. 4B, P < 0.05).A positive correlation between the mRNA expression levels of β-globin and γ-globin was shown in Fig. 4C (r = 0.804, P < 0.001).Furthermore, in pediatric β-thal patients without blood transfusion, circ-0008102 expression was positively correlated with the mRNA expression levels of β-globin (r = 0.534) and γ-globin (r = 0.781) (Fig. 4D, E, P < 0.05).Moreover, the mRNA expression levels of γ-globin were regulated when circ-0008102 overexpression or knockdown (Supplementary Fig. S1, P < 0.05).All data stated that circ-0008102 might be involved in the regulation of γ-globin expression in pediatric β-thal without blood transfusion.

GO and KEGG enrichment analyses of circ-0008102/ miRNA/mRNA network
The top 10 significantly enriched GO terms in the biological process (BP), cellular components (CC), and molecular function (MF) were shown in Supplementary Fig. S5A-C.In terms of BP, the top three enriched items were "regulation of DNA-templated transcription," "regulation of nucleic acidtemplated transcription," and "regulation of RNA biosynthetic process."In terms of CC, the top three enriched items were "transcription regulator complex," "chromation," and "chromosome."In terms of MF, "sequence-specific DNA binding," "transcriptional cis-regulatory region binding," and "transcription regulatory region nucleic acid binding" were enriched for the first three places.Besides, the top 10 significantly enriched pathways of KEGG enrichment analysis were shown in Supplementary Fig. S5D.Among these pathways, "TGF-beta signaling pathway," "stem cells pluripotency," "Th17 cell differentiation," and "cellular senescence" were closely correlated with the pathogenesis of pediatric β-thal without blood transfusion.

Discussion
With the deepening of the research on the underlying molecular basis of β-thal, its clinical treatment methods have been updated and gradually improved, and more and more biomarkers related to this disease have been reported, including erythrocyte transferrin receptor 1, alpha-hemoglobin stabilizer protein, heat shock protein 70, and CD35 and CD55 [16][17][18].In recent years, ncRNAs, including miRNAs and long non-coding RNAs (lncRNAs), are regarded as a class of promising blood-based biomarkers for detection of human diseases including β-thal [19].For example, Leecharoenkiat et al. reported that plasma upregulated miR-451 expression was a novel hemolytic marker in β 0 -thal/HbE patients and was associated with the degree of hemolysis and promoted erythropoiesis [20].Kuno et al. demonstrated that elevated miR-125b expression from activated phagocytic monocytes was a genetic modifier related to anemia severity in β-thal patients [21].El-Khazragy et al. suggested that plasma-dysregulated miR-let-7d, miR-122, and miR-200 served as novel and non-invasive predictor biomarkers for cellular damage under condition of tissue iron excess in transfusion-dependent β-thal patients [22].Recently, a cross-sectional study concluded that circulating MALAT1 (metastasis-associated lung adenocarcinoma transcript 1) and MIAT (myocardial infarction-associated transcript) were promising predictors for blood transfusion status in β-thal patients [23].Nevertheless, the potential biomarkers of β-thal remain largely unknown.
Accumulated evidence has suggested that the expression of circRNAs is tissue-and cell type-specific and specifically regulates disease progression [7].Compared with proteins, miRNAs, lncRNAs, and circRNAs are ideal candidate biomarkers in human diseases, as they are abundant and stable in body fluids, including peripheral blood, serum, cerebrospinal fluid, urine, gastric fluid, and joint fluid [8].Therefore, it is of great significance to identify circRNA-related biomarkers in β-thal.In this study, we found upregulated circ-0008102 expression was associated with the pathogenesis of pediatric β-thal and firstly demonstrated that circ-0008102 might be an effective biomarker for detection of pediatric β-thal without blood transfusion.
LCOR, as a master transcription factor, is involved in converting hemogenic endothelium into hematopoietic stem and progenitor cells [24].It is also reported that LCOR negatively regulates early adipogenesis by inhibiting C/EBPβ transcriptional activity [25].circ-0008102 is derived from LCOR gene locus, which is an exonic circRNA with 515 nt in length.Recently, dysregulated circRNAs have received more attention as diagnostic or predictive biomarkers for malignant tumors, nervous system diseases, cardiovascular disease, and endocrine diseases [26,27].However, the role of circRNAs in β-thal is poorly known.Herein, 59 pediatric β-thal patients were collected and their hematological parameters and biochemical indicators were consistent with the phenotypes of TM or TI.Encouragingly, we found that the expression levels of circ-0008102 in peripheral blood from pediatric β-thal patients without blood transfusion were significantly higher than that in pediatric β-thal patients with blood transfusion and healthy controls.The Fig. 5 The subcellular localization of circ-0008102.Fluorescence in situ hybridization (FISH) found that circ-0008102 was mainly located in the cytoplasm of K562 cells, scale bar = 20 μm data were consistent with upregulated circ-0008102 expression in a circRNA microarray profile of GSE196682 dataset (adult β-thal patients vs. healthy controls), indicating high credibility of our data.Furthermore, we found circ-0008102 expression was positively correlated with the levels of RBC and HbF, but negatively corrected with the levels of HbA and Cr.Blood transfusion affected the expression of circ-0008102 in pediatric β-thal patients probably for the following reasons: (1) due to the low expression of circ-0008102 in healthy controls, blood transfusion led to a decrease in the expression of circ-0008102 in the peripheral blood of pediatric β-thal patients; (2) due to the positive correlation between circ-0008102 and HbF, HbF decreased after blood transfusion, which caused a decrease in circ-0008102 expression as a result of feedback regulation; (3) due to the negative correlation between circ-0008102 and HbA, an increase in HbA after blood transfusion led to a decrease in circ-0008102.All these strong unraveled upregulated circ-0008102 expression was associated with the pathogenesis of pediatric β-thal without blood transfusion.Therefore, the detailed biological role of circ-0008102 in pediatric β-thal deserved further study.
circRNAs display good diagnostic performance, not only because of their high specificity and sensitivity but also because of their great stability in peripheral blood [28].Currently, the predictive role of circRNAs in β-thal remains to be blank.To our knowledge, it is the first time to identify the diagnostic power of circ-0008102 in pediatric β-thal patients in our study.The peripheral blood circ-0008102 revealed a significant diagnostic value on pediatric β-thal patients without blood transfusion when discriminated with pediatric β-thal patients with blood transfusion and healthy controls (AUC 0.733 and 0.711, respectively).These results suggested that circ-0008102 had potential as a novel biomarker in the detection of pediatric β-thal without blood transfusion.Without doubt, future multi-center, large-scale study is needed to verify the diagnostic value of circ-0008102 in pediatric β-thal.
Increasing evidence showed that exonic circRNAs served as sponges of miRNAs to regulate the expression of target mRNAs [38].We also found circ-0008102 was mainly located in cytoplasm of K562 cells by FISH assay, confirming circ-0008102 functioned as miRNA sponges.To better understand the regulatory mechanism of circ-0008102 in this study, we constructed the circ-0008102/miRNA/ mRNA network based on bioinformatics analysis and found circ-0008102 might act as ceRNAs to capture miR-372-3p, miR-329-5p, miR-198, miR-152-5p, and miR-627-3p and subsequently regulated the 651 mRNA expression.Previous studies had demonstrated that the five miRNAs could be detected in peripheral blood [39][40][41][42][43], suggesting that they might be aberrantly expressed in pediatric β-thal.The negative regulatory effect of circ-0008102 on the five miRNAs was verified in K562 cells by overexpression and knockdown of circ-0008102, confirming circ-0008102/miRNA/ mRNA network.The GO and pathway enrichment analyses indicated that these mRNAs were involved in DNA binding and transcription regulatory region binding and were associated with Th17 cell differentiation and stem cell pluripotency signaling pathways.All of these data indicated that circ-0008102 might be an effective biomarker for pediatric β-thal without blood transfusion.
(5) In addition to regulating the activation of γ-globin, the role of circ-0008102 in the progression of pediatric β-thal remains largely unknown.Future research will use in vitro and in vivo experiments to address these issues.
In conclusion, our data provide preliminary evidence that the expression level of circ-0008102 is upregulated in peripheral blood of pediatric β-thal without blood transfusion, and circ-0008102 may serve as a novel clinical biomarker for detection of pediatric β-thal without blood transfusion.circ-0008102 participates in the pathogenesis of β-thal through regulating γ-globin expression, which needs to be investigated further.

Fig. 1
Fig. 1 Identification of the circular structure of circ-0008102.A The genomic locus and backspliced junction of circ-0008102 were shown in schematic diagram.B The junction sequences of circ-0008102 were validated by Sanger sequencing.C RT-PCR analysis of circ-0008102 and LCOR mRNA expression in cDNA and gDNA from pediatric β-thal using divergent and convergent primers, respectively.D qRT-PCR analysis of circ-0008102 and LCOR mRNA levels in K562 cells under treatment with RNAse R. RT-PCR, reverse transcription PCR; LCOR, ligand dependent nuclear receptor corepressor; gDNA, genomic DNA; qRT-PCR, quantitative real-time PCR; RNase R, Ribonuclease R. Error bars represented the means ± SD. **P < 0.01, ns P > 0.05