Abstract
Introduction: Unlike homozygous hemoglobin SS (HbSS) disease, stroke is a rare complication in hemoglobin SC (HbSC) disease. However, recent studies have demonstrated a high prevalence of silent stroke in HbSC disease. The factors associated with stroke and cerebral vasculopathy in the HbSC population are unknown. Methods: We conducted a retrospective study of all patients with sickle cell disease treated at the University of Missouri, Columbia, over an 18-year period (2000–2018). The goal of the study was to characterize the silent, overt stroke, and cerebral vasculopathy in HbSC patients and compare them to patients with HbSS and HbS/β thalassemia1 (thal) in this cohort. We also analyzed the laboratory and clinical factors associated with stroke and cerebral vasculopathy in the HbSC population. Results: Of the 34 HbSC individuals, we found that the overall prevalence of stroke and cerebral vasculopathy was 17.7%. Only females had evidence of stroke or cerebral vasculopathy in our HbSC cohort (33.3%, p = 0.019). Time-averaged means of maximum velocities were lower in the HbSC group than the HbSS group and did not correlate with stroke outcome. Among HbSC individuals, those with stroke and cerebral vasculopathy had a marginally higher serum creatinine than those without these complications (0.77 mg/dL vs. 0.88 mg/dL, p = 0.08). Stroke outcome was associated with recurrent vaso-occlusive pain crises (Rec VOCs) (75 vs. 25%, p = 0.003) in HbSC patients. The predominant cerebrovascular lesions in HbSC included microhemorrhages and leukoencephalopathy. Conclusion: There is a distinct subset of individuals with HbSC who developed overt, silent stroke, and cerebral vasculopathy. A female predominance and association with Rec VOCs were identified in our cohort; however, larger clinical trials are needed to identify and confirm specific clinical and laboratory markers associated with stroke and vasculopathy in HbSC disease.
Introduction
Children and adults with sickle cell disease (SCD) suffer from overt stroke, the prevalence of which has been reported to be 11% in individuals <20 years of age [1]. More commonly, they suffer from silent cerebral infarction (SCI) that has a reported prevalence of 37.1% in children <4 years of age [2]. SCI is associated with neurocognitive deficits in children and is associated with subsequent risk of overt stroke [3-5].
Unlike homozygous sickle cell disease (HbSS), the reported prevalence of overt and silent stroke in hemoglobin SC (HbSC) disease is low, at 0 and 5.8%, respectively, in Cooperative Study of Sickle Cell Disease (CSSD) [6]. However, in a recent study, the prevalence of SCI in HbSC was found to be 13.5% [7].
Given the association of SCI with both subsequent stroke [3] and overall decline in neurocognitive function [5], the detection and institution of comprehensive interventions, including judicious use of disease-modifying agents, are crucial [4]. Despite recent advances, both causal and progression factors associated with silent stroke and cerebral vasculopathy in HbSS disease are poorly understood. Comparatively lesser number of studies has been conducted looking into the prevalence, laboratory, and clinical factors associated with overt, silent stroke, and cerebral vasculopathy in HbSC disease, the second most common form of SCD. Currently, most of our understanding of HbSC is from extrapolation of data from HbSS disease [8].
The aim of our study was to characterize the silent and overt stroke, cerebral vasculopathy, and neurovascular complications in children and adults with HbSC and compare them to HbSS and HbS/β thalassemia (thal) in the Central Missouri Cohort (University of Missouri [MU]-SCD). Additionally, we examined the SCD-related complications and laboratory factors associated with cerebrovascular complications in this population.
Methods
In this retrospective cross-sectional study, pediatric and adult SCD patients treated at MU between 2000 (first availability of the electronic medical record) and 2018 were included. Retrospective chart review was conducted on all HbSC, HbSS, and HbS/β thal patients treated at MU medical center. The Institutional Review Board for MU, Columbia, MO, approved this study. Hemoglobin electrophoresis and, in some cases, β-globin gene sequencing was used to identify patients with HbSS, HbSC, and HbS/β thal. The classification of silent or overt stroke was done based on reported neurological events on chart review. Other SCD-related complications, such as pain crisis, acute chest syndrome, cutaneous ulceration, gallbladder disease, pulmonary hypertension, renal complications, and priapism, were reviewed throughout the course of the patient’s medical history. SCD-related vaso-occlusive pain crises (VOC) were defined as an acute episode of pain that resulted in a visit to our medical facility and treatment with oral or parenteral narcotic medication or with a parenteral nonsteroidal anti-inflammatory drug.
Recurrent VOC (Rec VOC) was defined as 2 or more episodes of VOCs warranting admission to the hospital for control of pain in 12 months preceding the identification of stroke or cerebral vasculopathy. Elevation in serum unconjugated (indirect) fraction of the bilirubin was defined as unconjugated hyperbilirubinemia. Indirect hyperbilirubinemia was defined as unconjugated bilirubin >1.0 mg/dL at our institution.
The use of disease-modifying treatments with hydroxyurea (HU) and red blood cell exchange transfusion (RBCx) was assessed in all SCD individuals in the cohort. Children with SCD between ages 10–20 years with history of surveillance magnetic resonance imaging (MRI) scans were assessed for indications, including headache, seizures, and neurological deficits. All patients between the ages of 2 years–16 years also received annual transcranial doppler ultrasound (TCD) per institutional guidelines. TCD velocities were analyzed and verified by 2 independent radiologists on separate occasions. Three consecutive time-averaged means of maximum velocities (TAMMV) were obtained bilaterally at 3 locations in the middle cerebral artery (MCA) and 1 location in the anterior cerebral artery (ACA), as described previously [22]. Per institutional protocol, MRI was performed in SCD individuals if they had a history of headaches, neurological deficits, or an abnormal TCD. Clinical, laboratory, TCD, and MRI/ magnetic resonance angiography data for clinically documented overt stroke, SCI, and vasculopathy were analyzed.
Indication for MRI scanning for adults included transient or focal neurological events, severe headaches, or seizures. Routine surveillance MRI was not performed in adults at our institution. All MRIs were reviewed by a board-certified radiologist and re-reviewed in a blinded manner by a neuroradiologist. The study neuroradiologist assessment was compared to the original report for agreement.
Imaging Acquisition and Analysis
Patients underwent MRI of the brain on Siemens Aera 1.5T or Skyra 3T MR scanners (Siemens Healthineers, Malvern, PA) with a standard 12-channel head coil. Sagittal three-dimensional magnetization prepared rapid acquisition gradient echo (MP-RAGE) T1-weighted, coronal T2-weighted, and axial fluid-attenuated inversion recovery (FLAIR), susceptibility-weighted, and diffusion-weighted images were acquired as per the clinical protocols at the time of the examination. The neuroradiologist reviewed the MR images to classify the participants as having deep lacunar infarcts, encephalomalacia, brain atrophy, leukoencephalopathy/microvascular changes, parenchymal intracranial hemorrhage, chronic microhemorrhages, and superficial hemosiderosis.
We considered lacunes as small focal <15 mm hyperintense/CSF signal lesions in the basal ganglia or deep white matter of the brain. Such lacunes and encephalomalacia that were not associated with documented neurological deficits were termed as silent cerebral infarction [9]. Microhemorrhages were defined as punctate foci of hypo-intensity on susceptibility-weighted image and distinct from vascular flow voids and leptomeningeal hemosiderosis. The latter was considered when linear/curvilinear magnetic susceptibility within the subarachnoid space.
Statistical Analysis
Descriptive analysis was used and expressed as means +/− standard deviation (SD), medians (range), and numbers (%). Fisher’s exact tests for categorical data and Mann-Whitney U test for continuous data analysis were used. A p value <0.05 was considered to be statistically significant. The analysis was conducted using PRISM and SAS software (9.3).
Results
Demographic Characteristics of the MU-SCD Cohort
We analyzed charts of 131 children and adults with SCD, with an age range of 1–65 years. On average, SCD patients received follow-up at MU for 1.7 decades of time at the medical center. Of the 131 individuals with SCD who sought care at MU, 28 had incomplete electronic medical data for assessment, and therefore, 103 patients were included in the statistical analysis.
We first determined the cerebrovascular changes in the HbSC group and compared them to the HbSS and HbS/β thal group. Of the 38 patients with HbSC, 34 patients had evaluable clinical data and 9 patients had prior brain MRI. In patients with HbSS disease, 60 patients had evaluable data and 30 patients had prior MRI. Of the 9 patients with HbS/β subtype, 5 patients had prior MRI of the brain.
No overt strokes were observed in children <10 years of age in the MU-SCD cohort (shown in Table 1). It was observed that 32% of children between 10 and 20 years had a cerebrovascular lesion, and the prevalence was highest compared to other age-groups. Among individuals with SCD surviving into the sixth decade of life, 62.5% had evidence of MRI lesions in the brain consistent with stroke. A significantly higher number of patients in the stroke category were on HU (76.9 vs. 23.1%, p = 0.001) and RBCx therapy (61.5 vs. 38.5%, p = 0.02).
Stroke and Vasculopathy in HbSC Disease
In the MU-SCD cohort, stroke and vasculopathy were observed in 17.7% (n = 6) in the HbSC group. In comparison, stroke and vasculopathy were observed in 30% (n = 18) of the HbSS and 22.2% (n = 2) of the HbS/β thal patients (shown in Table 2).
Silent stroke and cerebral vasculopathy were more frequently observed than overt stroke in the HbSC group (shown in Table 2). As expected, an overall higher prevalence of both silent and overt stroke was observed in the HbSS group than patients with HbSC and HbS/β thal. The percentage of cerebral vasculopathy was observed to be higher in HbSC than HbSS and HbS/β thal group in our cohort; however, these results were not statistically significant (shown in Table 2). Stroke and cerebral vasculopathy were identified only in women with HbSC in the MU-SCD cohort (33 vs. 0%, p = 0.019) (shown in Tables 3, 5). All the patients with stroke or vasculopathy were of African-American ethnicity (shown in Table 3).
TCD Velocities and Stroke in MU-SCD
As a part of MU institutional policy, TCD was routinely performed as part of annual surveillance in all children with SCD (HbSS, HbSC, and HbS/β thal) <16 years of age (shown in online supplementary. Table 1; see www.karger.com/doi/10.1159/000519360 for all online suppl. material).
TCD was performed in 25 (65.7%) patients with HbSC as a part of annual surveillance program in children. Three consecutive results were used to calculate the mean velocities as shown (shown in Fig. 1). The TAMMV for ACA, MCA, and posterior cerebral artery is shown (shown in Fig. 1). The TAMMV was 80. 16 cm/s [+/− 28.95 cm/s, SD and 4.46 cm/s, standard error of mean (SEM)] in the right MCA; 47.34 cm/s (+/− 20.61 cm/s, 3.072 cm/s; SD and SEM, respectively), in the right ACA and 43.86 cm/s (+/− 16.71 cm/s and 2.74 cm/s; SD and SEM, respectively), in the right posterior cerebral artery. On the left side, the mean velocities were 79.12 cm/s (+/− 24.17 cm/s, 3.56 cm/s; SD and SEM, respectively) in MCA, 50.13 cm/s (+/− 16.4 cm/s and 2.63 cm/s; SD and SEM, respectively) in the ACA and 39.53 cm/s (+/− 14.65 cm/s, 2.35 cm/s; SD and SEM, respectively).
TCD velocities in HbSC disease. TAMMV, timed averaged mean of maximum velocity; mRMCA, mid-right middle cerebral artery; RACA, right anterior cerebral artery, rPCA, right posterior cerebral artery; mLMCA, mid-left middle cerebral artery; LACA, left anterior cerebral artery, LPCA, left posterior cerebral artery; HbSC, hemoglobin SC; TCD, transcranial Doppler ultrasound.
TCD velocities in HbSC disease. TAMMV, timed averaged mean of maximum velocity; mRMCA, mid-right middle cerebral artery; RACA, right anterior cerebral artery, rPCA, right posterior cerebral artery; mLMCA, mid-left middle cerebral artery; LACA, left anterior cerebral artery, LPCA, left posterior cerebral artery; HbSC, hemoglobin SC; TCD, transcranial Doppler ultrasound.
No patient with HbSC had TAMMV above 170 cm/s in this cohort (shown in Fig. 1). The left ACA was not visualized in 1 patient, but TAMMV was <150 cm/s in all other intracranial vessels. No elevation in TAMMV was observed in the preceding 12 months in 1 child with HbSC who developed overt stroke.
Laboratory Parameters Associated with Stroke and Cerebral Vasculopathy in HbSC Disease
The mean hemoglobin at steady state was above 10 g/dL in both groups of patients with and without stroke in HbSC disease (shown in Table 4). Although not statistically significant, HbF was found to be lower in the group with stroke than those without a history of stroke. Higher serum creatinine was observed in individuals with stroke and vasculopathy than those who did not have stroke (0.77 mg/dL vs. 0.88 mg/dL, p = 0.08), but these findings were not statistically significant. In accordance with this finding, there was evidence of lower estimated glomerular filtration rate in the stroke and cerebral vasculopathy group than patients without stroke and cerebral vasculopathy in HbSC disease. Other laboratory parameters, including platelet count, white blood cell count, and aspartate aminotransaminases, were not statistically different in both groups.
Sickle Cell-Related Complications Associated with Stroke in HbSC Disease
We then proceeded to analyze SCD-related complications associated with stroke outcome in the HbSC individuals (shown in Fig. 2). We observed that some laboratory markers of hemolysis (anemia, reticulocytosis, unconjugated hyperbilirubinemia, and increased urobilinogen) were not significantly different in the overt and silent stroke patients with HbSC compared to those without stroke (shown in Table 4). Complications associated with the hemolytic phenotype, such as pulmonary hypertension, priapism, cutaneous ulcer, and gallbladder stones were not statistically different in both groups. Osteonecrosis and proliferative retinopathy (62.5 vs. 37.5%, p = 0.9) of SCD were not associated with stroke outcome in HbSC disease (shown in Fig. 2d). Rec VOC in HbSC individuals was associated with stroke outcome compared to those without Rec VOC (75 vs. 25%, p = 0.003) (shown in Fig. 2e).
Sickle cell-related complications associated with stroke. HbSC, hemoglobin SC; VOC, vaso-occlusive crisis; ACS, acute chest syndrome; UCB, unconjugated hyperbilirubinemia; RP, retinitis proliferans; Rec VOEs, recurrent vaso-occlusive episodes.
Sickle cell-related complications associated with stroke. HbSC, hemoglobin SC; VOC, vaso-occlusive crisis; ACS, acute chest syndrome; UCB, unconjugated hyperbilirubinemia; RP, retinitis proliferans; Rec VOEs, recurrent vaso-occlusive episodes.
MRI in HbSC Disease
We analyzed the cerebrovascular lesions in patients in accordance with the criteria stated in the methods section. As brain MRI is not routinely performed in patients with HbSC, clinical documentation was reviewed for determination of need for MRI. Of the 9 patients who obtained MRI, 8 had a history of severe, recurrent headaches, and 1 patient had a focal neurological deficit.
Among the HbSC patients with MRI, 1 (11.1%) had lacunar infarcts (shown in Fig. 3a), 1 (11.1%) showed focal encephalomalacia (shown in Fig. 3b), and 1 (11.1%) cerebral atrophy (shown in Fig. 3c). The brain MRI revealed leukoencephalopathy in 3 (33.3%) (shown in Fig. 3d) and microhemorrhages in 3 (33.3%) (shown in Fig. 3e) subjects. No hemosiderosis was observed in HbSC patients in the cohort (online suppl. Table 2).
Spectrum of MRI findings in patients with HbSC disease. MRI, magnetic resonance imaging; HbSC, hemoglobin SC; FLAIR, fluid-attenuated inversion recovery. a Axial T2-weighted image shows punctate hyperintense foci within the left caudate head and right medial thalamus, compatible with chronic lacunar infarcts. b Axial T2-weighted shows encephalomalacia and hemosiderin staining within the left periventricular white matter. c Axial T2-weighted image shows global cerebral volume loss with mild ventriculomegaly. d Axial FLAIR image shows scattered punctate hyperintense foci within the white matter of the bilateral cerebral hemispheres, most compatible with chronic microvascular changes. e Axial susceptibility-weighted image shows multiple scattered foci of magnetic susceptibility within the subcortical white matter of the bilateral cerebral hemispheres consistent with chronic microhemorrhages.
Spectrum of MRI findings in patients with HbSC disease. MRI, magnetic resonance imaging; HbSC, hemoglobin SC; FLAIR, fluid-attenuated inversion recovery. a Axial T2-weighted image shows punctate hyperintense foci within the left caudate head and right medial thalamus, compatible with chronic lacunar infarcts. b Axial T2-weighted shows encephalomalacia and hemosiderin staining within the left periventricular white matter. c Axial T2-weighted image shows global cerebral volume loss with mild ventriculomegaly. d Axial FLAIR image shows scattered punctate hyperintense foci within the white matter of the bilateral cerebral hemispheres, most compatible with chronic microvascular changes. e Axial susceptibility-weighted image shows multiple scattered foci of magnetic susceptibility within the subcortical white matter of the bilateral cerebral hemispheres consistent with chronic microhemorrhages.
Discussion
The prevalence of SCI in HbSC has been reported to be around 5.8% compared to 21.8% in HbSS children between 6 and 18 years of age in the CSSD cohort [1, 6]. More recent studies have shown that the prevalence of SCI in HbSC might be higher than previously reported [7]. As the number of individuals with HbSC obtaining routine brain MRI is low, it has been difficult to quantify the overall prevalence of silent stroke and cerebral vasculopathy in this population. Additionally, the clinical and laboratory variables associated with stroke in HbSC have not been previously studied. To address these gaps in our understanding, we retrospectively analyzed the clinical characteristics and laboratory and clinical parameters associated with both overt and silent stroke and vasculopathy in HbSC individuals.
Our study identified a higher than previously reported prevalence of SCI and cerebral vasculopathy in the HbSC population compared to previous CSSD reports [6]. SCI was more prevalent than overt stroke in the HbSC population, and interestingly similar white matter lesions have been reported in 1.3% of healthy children as incidental findings [10]. In this analysis, the overall prevalence of SCI and cerebral vasculopathy was around 14.7% whereas the prevalence of similar lesions was estimated to be 13.3% in the HbSS group which is in accordance with the prevalence that has been reported in other prospective studies [6]. The reported prevalence of cerebrovascular accidents by 20, 30, and 45 years of age was estimated to be 11, 15, and 24%, respectively, in HbSS individuals but was significantly lower in the HbSC individuals at 2, 4, and 10%, respectively [11]. Similar to our findings, some recent studies have also shown a higher incidence of SCI in HbSC population [7, 12, 13]. It is unclear if the wider use of high-resolution MRI is the cause of higher detection of these SCIs.
There were no dominant CNS lesions in HbSC disease compared to HbSS disease in our cohort. There was evidence of a slightly higher prevalence of encephalomalacia, atrophy, and microhemorrhage in HbSC than HbSS individuals (online suppl. Table 2). The most common MRI lesions reported in the brain of adults with HbSS disease has been leukoencephalopathy, followed by atrophy and encephalomalacia [14].
Our preliminary data indicate that the pathogenesis of stroke in patients with HbSC may be different from that seen in HbSS individuals. Patients with stroke or cerebral vasculopathy had mean hemoglobin above 10 g/dL, indicating that anemia is less likely a causative factor in these patients. Based on our observations, it is likely that hyperviscosity may play a role in the pathogenesis of both silent and overt strokes in HbSC individuals [15-18].
Although not statistically significant, higher concentration of HbF and lower creatinine were detected in the HbSC individuals without stroke than patients with stroke. Other indicators of hemolysis, including reticulocyte count, aspartate aminotransferase, and total bilirubin, were not significantly different between the 2 groups. We did not routinely perform LDH testing on our patients with SCD, and therefore, these data were not available to us. Due to small size of study population, it is difficult to draw a definitive conclusion from these observations and more studies are needed.
In HbSS patients, both ACS and pain crises were previously observed to be negatively associated with silent stroke [19]. Anemia and hemolysis were found to be associated with silent stroke in HbSS, in previous studies [20]. In another study in HbSS individuals, it was found that SCI was associated with lower Hb and abnormal magnetic resonance angiography whereas VOC or acute chest syndrome did not positively correlate with stroke [20].
In our HbSC cohort, we found that patients with a previous history of Rec VOCs had a positive correlation with stroke outcome. However, the other sickle cell-related complications, like retinopathy, priapism, pulmonary hypertension, and cutaneous ulcer did not correlate with stroke outcome. Even though a single documented vaso-occlusive episode did not correlate with stroke outcome, recurrent vaso-occlusive episodes were significantly associated with stroke and cerebral vasculopathy.
Factors known to play a role in HbSC disease pathology could in part play a role in HbSC stroke pathology. These include abnormal micro-and macrovascular function [21] and unique hemorheological properties of HbSC red blood cells (RBC) [22]. Impaired RBC deformability and increased RBC aggregation have been distinctly observed in HbSC disease [22].
TCD results did not correlate with stroke in HbSC, and the individual who developed overt stroke had normal TCD in the preceding year. This lack of TCD correlation was a notable difference observed in HbSC patients compared to HbSS and HbS/β thal patients from historical trials [23, 24]. Similar lack of correlation between TCD and SCI has been reported in other SCD cohorts [20, 25]. On a broader perspective, TCDs have not been standardized in HbSC disease. TAMMV has been found to be significantly lower in HbSC (<130 cm/s) than HbSS individuals from the Brazilian cohort [26, 27]. Further studies are warranted to assess if higher TAMMV is associated with overt or silent stroke in HbSC disease.
Stroke and vasculopathy were seen only in women in our HbSC cohort (shown in Table 5). This female predominance has not been noted in our HbSS cohort or in the historical HbSS controls [6]. Due to small numbers, it is difficult to ascertain at this time if female gender in itself contributes to the risk of stroke in HbSC disease.
Treatment considerations for overt and silent stroke in HbSC have mostly followed HbSS disease guidelines, even though the etiopathology of the disease might be different. Transfusion, HU, and phlebotomy have been undertaken in HbSC disease in small clinical trials [28-31]. Currently, there is an ongoing debate on the optimal management of SCI in SCD [4, 32]. Furthermore, there are no available guidelines regarding the optimal management of SCI and overt stroke in HbSC individuals, and a multidisciplinary approach is warranted looking into all options, including stem cell transplantation, transfusions, phlebotomy, and HU [33].
Limitations of our study include the retrospective nature and the small number of patients who underwent MRI scans. There was also a selection bias as symptomatic patients were more likely to obtain brain MRI. Significant causal inferences would not be possible with small numbers and need to be tested in a larger prospective trial. Additionally, accurate prevalence cannot be ascertained by retrospective study and can be determined only by prospective screening. The incidence of stroke in our HbSC cohort was higher than what has been previously reported, but it is difficult to compare the incidence rate to historical controls by statistical analyses. There is also a possibility of ascertainment bias, with our radiologist using more advanced MRI techniques than historical cohorts, leading to higher detection rate of stroke and vasculopathy. Higher HU and RBCx use in the stroke category is likely to be confounded as an indication rather than cause-effect. More longitudinal studies are required to understand the onset and evolution of these cerebrovascular lesions in HbSC individuals.
In conclusion, there is a recognizable subset of patients with HbSC who suffered from both silent and overt strokes. Stroke in HbSC has distinct differences compared to stroke in HbSS individuals. More studies are needed to understand the pathophysiology and risk factors associated with stroke in HbSC individuals. Earlier identification of at-risk individuals through clinical and laboratory biomarkers and institution of preventive and treatment strategies is warranted.
Acknowledgments
We would like to thank Ward Hagar, MD, for his editorial comments and Shalini Sathi, for editing the manuscript.
Statement of Ethics
This study was approved by the Institutional Review Board, University of Missouri- Columbia, MO, USA (IRB#2004777 HS). Health Insurance Portability and Accountability Act (HIPAA) waiver of consent and assent was obtained for this study (MU IRB contact: cassanij@health.missouri.edu).
Conflict of Interest Statement
Speaker’s Bureau (V.B.) for Novartis. Authors (B.K.S., Y.Y., M.R.W., L.N., B.G., T.A., C.L.S.) declare no competing conflict of interest.
Funding Sources
This project was funded by Missouri State Grant for Hemoglobinopathy (AOC17380120-B.K.S.) and Leda Sears Trust Fund, University of Missouri (B.K.S.). Y.Y. was partly supported by NIH grant (K12HDO43451) from Eunice Kennedy Shriver National Institute of Child and Human Development of Building Interdisciplinary Research Careers in Women’s Health (BIRCWH) scholar. Some of the findings of this study were presented in abstract form at Annual American Society of Hematology Meeting, 2019.
Author Contributions
All authors contributed significantly to the article. B.K.S. and C.L.S. conceptualized the study, collected the data, and wrote the manuscript; B.G., V.B., and T.A. wrote and edited the manuscript; M.R.W. and L.N. collected the data and edited the manuscript; Y.Y. analyzed the data and edited the manuscript.