Serum anti-AP3D1 antibodies are risk factors for acute ischemic stroke related with atherosclerosis

Atherosclerosis has been considered as the main cause of morbidity, mortality, and disability worldwide. The first screening for antigen markers was conducted using the serological identification of antigens by recombinant cDNA expression cloning, which has identified adaptor-related protein complex 3 subunit delta 1 (AP3D1) as an antigen recognized by serum IgG antibodies of patients with atherosclerosis. Serum antibody levels were examined using the amplified luminescent proximity homogeneous assay-linked immunosorbent assay (AlphaLISA) using a recombinant protein as an antigen. It was determined that the serum antibody levels against AP3D1 were higher in patients with acute ischemic stroke (AIS), transient ischemic attack, diabetes mellitus (DM), cardiovascular disease, chronic kidney disease (CKD), esophageal squamous cell carcinoma (ESCC), and colorectal carcinoma than those in the healthy donors. The area under the curve values of DM, nephrosclerosis type of CKD, and ESCC calculated using receiver operating characteristic curve analysis were higher than those of other diseases. Correlation analysis showed that the anti-AP3D1 antibody levels were highly associated with maximum intima-media thickness, which indicates that this marker reflected the development of atherosclerosis. The results of the Japan Public Health Center-based Prospective Study indicated that this antibody marker is deemed useful as risk factors for AIS.


Methods
Patients and HDs' sera. This study was approved by the Local Ethical Review Board of Chiba University, Graduate School of Medicine (Chiba, Japan) and the review boards of the participating hospitals. All experimental procedures were performed in accordance with the Declaration of Helsinki.
Blood samples were collected from patients who had provided their informed consent. Each serum sample was centrifuged at 2000×g for 10 min at 4°C, and the supernatant was stored at −80°C until use. Repeated thawing and the freezing of samples were avoided.
Serum samples of patients with AIS, transient ischemic attack (TIA), deep and subcortical white matter hyperintensity (DSWMH), asymptomatic cerebral infarction (asymptCI), and chronic-phase cerebral infarction (cCI) were provided by the Chiba Prefectural Sawara Hospital, Chiba Rosai Hospital, and Chiba Aoba Municipal Hospital, and sera of patients with DM and CVD were obtained from Chiba University Hospital. The stroke subtype of each patient was determined according to the criteria of the Trial of Org 10,172 in Acute Stroke Treatment classification system 24 . In this study, large-artery atherosclerosis or small-artery occlusion (lacune) were included as AIS or cerebral infarction. Serum samples associated with AIS, TIA, and CVD were obtained within 2 weeks after disease onset. Sera of CKD patients were obtained from the Kumamoto cohort 25,26 , and Chiba University Hospital provided the serum samples of patients with ESCC and colorectal carcinoma (CRC). Sera of HDs were obtained from Chiba University Hospital, Port Square Kashiwado Clinic, and Chiba Prefectural Sawara Hospital. Sample of HDs from Port Square Kashiwado Clinic and Chiba Prefectural Sawara Hospital were selected from individuals who exhibited no abnormalities in cranial magnetic resonance imaging.
Sera of patients with atherosclerosis were used to search for antigens that could be recognized by the serum IgG antibody. The library used was a Uni-ZAP XR cDNA phage library containing a human microvascular endothelial cell cDNA library (Stratagene; Agilent Technologies, Inc., La Jolla, CA), which was infected into Escherichia coli (E. coli) XL1-Blue MRF′. Proteins were then transferred onto nitrocellulose membranes [NitroBind, Osmonics Inc., Minnetonka, MN)], which were pretreated with 10 mM isopropyl-β-D-thiogalactoside (IPTG) (Wako Pure Chemicals, Osaka, Japan) for 30 min. The membranes were blocked for 1 h with 1% protease-free bovine serum albumin (Wako Pure Chemicals), 20 mM Tris-HCl (pH 7.5), 0.15 M NaCl, and 0.05% Tween-20, treated overnight to 1:2000 diluted sera from the patients, and then incubated for 1 h with 1:5000 diluted alkaline phosphatase-conjugated goat anti-human IgG (Jackson ImmunoResearch Laboratories, Inc., West Grove, PA). Finally, using a color development solution [0.3 mg/ml nitroblue tetrazolium (Wako Pure Chemicals), 0.15 mg/ ml 5-bromo-4-chloro-3-indolyl-phosphate (Wako Pure Chemicals), 100 mM Tris-HCl (pH 9.5), 100 mM NaCl, and 5 mM MgCl 2 ], the positive reactions could be identified. The resulting positive antibody were then re-cloned twice in order to obtain monoclonality 13-17,27-29 . Sequence analysis of identified clones. Monoclonal phage cDNA clones were converted to pBluescript phagemids by in vivo excision using the ExAssist helper phage (Stratagene; Agilent Technologies, La Jolla, CA). Plasmid DNA was obtained from the E. coli SOLR strains transformed by the phagemids. Homology search of the inserted and sequenced cDNAs using a public database provided by the the National Center for Biotechnology Information (https:// blast. ncbi. nlm. nih. gov/ Blast. cgi) identified the genes. www.nature.com/scientificreports/ AP3D1 protein expression, extraction, and purification. The region of 2490-4347 of the AP3D1 gene was isolated and was recombined into the EcoRI/NotI site of pGEX-4T-1 (GE Healthcare Life Sciences, Pittsburgh, PA), followed by confirmation by DNA sequencing. The E. coli BL-21 transfected with pGEX-4T-1-AP3D1 was then treated with 0.1 mM IPTG at 37°C for 3 h to induce the expression of cDNA products. Then, the cells were lysed in BugBuster Master Mix (Merck KGaA, Darmstadt, Germany). With the Glutathione-Sepharose (GE Healthcare Life Sciences, Pittsburgh, PA) column chromatography according to the manufacturer's instructions, glutathione S-transferase (GST)-fused-AP3D1 protein was purified, as has been described previously 17,19,20 . AlphaLISA (amplified luminescence proximity homogeneous assay-linked immunosorbent assay). To evaluate the serum antibody levels, AlphaLISA was used. Initially, AlphaLISA was performed in 384-well microtiter plates (white opaque OptiPlate, PerkinElmer, Waltham, MA, USA) containing either 2.5 µl of 1:100 diluted serum or 2.5 µl of 10 µg/ml of GST and GST-AP3D1 protein (25 mM HEPES, pH 7.4, 0.1% casein, 0.5% Triton X-100, 1 mg/ml dextran-500 and 0.05% Proclin-300). The reaction mixture was then incubated at room temperature for 6-8 h. Secondly, after adding anti-human IgG-conjugated acceptor beads (2.5 µl at 40 µg/ml) and glutathione-conjugated donor beads (2.5 µl at 40 µg/ml), the mixture was further incubated at room temperature in the dark for 7-14 days. Chemical emissions were read on an EnSpire Alpha microplate reader (PerkinElmer) as previously described [18][19][20][21][22][23]30 . Specific reactions were then calculated by subtracting the alpha photon counts of the GST and buffer control from the counts of the GST-AP3D1 protein.
JPHC cohort analysis. The longitudinal association between plasma AP3D1 levels (using the above AlphaLISA detection antibody levels) and incident AIS were examined in the Japan Public Health Center-based prospective Study (JPHC). The study nested within JPHC cohort 31,32 , involving approximately 30,000 Japanese individuals aged 40-69 years at the baseline period of 1990-1994 whose plasma were stored. The antibody levels of AP3D1 protein were measured in 202 cases of AIS in the cohort developed between the baseline and 2008, and in 202 controls whose sex, age (within 2 years), date of blood sampling (within 3 months), time since last meal (within 4 h) and study location (Public Health Center area) were matched with the cases. We used a conditional logistic regression model to estimate the odds ratios (ORs) and 95% confidence intervals (CIs) for AIS with respect to the antibody levels of AP3D1 protein.
Statistical analysis. All statistical analyses were conducted using GraphPad Prism 5 (GraphPad Software, La Jolla, CA) and EZR software 33 . The Kruskal-Wallis test (Mann-Whitney U test with Bonferroni's correction applied) was used to evaluate differences among > 3 groups, and the Mann-Whitney U test was employed to determine significant differences between the two groups. Correlations were calculated using Spearman's correlation analysis and logistic regression analysis. The predictive values of the putative disease markers were assessed using a receiver operating characteristic (ROC) curve analysis, and the cutoff values were set to maximize the sums of sensitivity and specificity. All tests were two-tailed, and P values lower than 0.05 were considered to be statistically significant. The power calculation was performed using G-Power 3.1 software (Heinrich-Heine-Universität Düsseldorf).
Ethics declarations. The present study was approved by the Local Ethical Review Board of Chiba University, Graduate School of Medicine (Chiba, Japan) as well as the review boards of co-operating hospitals or institutes. Serum or plasma was collected from participants who had provided informed consent by following the protocols approved by their institutional ethical committees.

Results
Initial screening of AP3D1 antigens using SEREX. As an initial SEREX screening, sera of patients with atherosclerosis were used to search for antigens that could be recognized by serum IgG antibodies, one of which was AP3D1 (accession no. NM_003938.8). The region of 2490-4347 of the AP3D1 gene was then isolated and recombined into the EcoRI/NotI site of pGEX-4T-1, followed by confirmation by DNA sequencing. The cDNA was then expressed in E. coli, purified by affinity chromatography, and employed as an antigen in order to examine the serum antibody levels.
The levels of anti-AP3D1 antibodies were elevated in patients with AIS and TIA. . The levels of s-AP3D1-Abs were determined to be significantly higher in patients with AIS and TIA than those in HDs (Fig. 1a). When the ages of the subjects were matched to 65 years, the s-AP3D1-Ab levels were still significantly higher in patients with AIS than those in HDs ( Supplementary Fig. S1a). At a cutoff value equivalent to the average plus two SDs of the HD values, the s-AP3D1-Ab-positive rates in HDs and patients with AIS and TIA were 2.4%, 10.1%, and 10.4%, respectively (Table 1 lower  Elevation of s-AP3D1-Abs levels in patients with DM. The levels of s-AP3D1-Abs were also examined for DM. Sera of HD were obtained from Chiba University, whereas the sera of patients with DM were provided by the Chiba University Hospital. The average ages (± SDs) of the HDs and patients with DM were 45.20 ± 10.95 and 63.12 ± 12.04 years, respectively. The AlphaLISA results revealed that s-AP3D1-Ab levels were significantly higher in patients with DM than in the HDs (Fig. 2a). When the levels of s-AP3D1-Abs were compared between age-matched (60 years) HDs and patients with DM, the levels were also significantly higher in patients with DM than in HDs ( Supplementary Fig. S1b). When the positive samples for which the AlphaLISA counts exceeded the cutoff value were scored, the positive rates of s-AP3D1-Abs in the HDs and the patients with DM were 3.7% and 41.8%, respectively ( Table 2). The AUC value of s-AP3D1-Abs vs. DM was as high as 0.791 (Fig. 2b). Therefore, it can be concluded that the s-AP3D1-Ab levels were closely associated with DM.
The s-AP3D1-Abs levels were associated with CVD. For the next step, the antibody levels in samples from CVD patients were examined. The samples of CVD patients were obtained from Chiba University Hospital, and those in HDs were from Chiba University, Port Square Kashiwado Clinic, and Chiba Prefectural Sawara Hospital. The average ages (± SDs) of the HDs and CVD patients were 45.27 ± 11.20 and 66.07 ± 11.32 years, respectively. Compared with HDs, s-AP3D1-Abs levels were significantly higher in patients with CVD (Fig. 3a), and the s-AP3D1-Ab positivity rates in HDs and patients with CVD were 5.1% and 24.0%, respectively (Table 3). ROC curve analysis revealed that AUC of s-AP3D1-Abs for CVD was 0.758 (Fig. 3b).
The s-AP3D1-Ab levels were closely related to CKD. The antibody levels in the sera of CKD patients were examined, which were assumed to be closely related to atherosclerosis. The sera of the CKD patients were obtained from the Kumamoto cohort 22,23 , including 145 from patients with diabetic kidney disease (type 1 CKD), 32 from patients with nephrosclerosis (type 2 CKD), and 123 from patients with glomerulonephritis (type 3 CKD). The sera of HDs (82 specimens) were obtained from Chiba University, Chiba Prefectural Sawara www.nature.com/scientificreports/ Hospital, and the National Hospital Organization of Shimoshizu Hospital. Patients from all three groups of CKD were found to have significantly higher levels of s-AP3D1-Abs compared to that in HDs (Fig. 4a). The s-AP3D1-Ab-positive rates in HDs and patients with type 1, type 2, and type 3 CKD were 4.9%, 27.6%, 37.5%, and 22.8%, respectively (Table 4). ROC curve analysis revealed AUC of s-AP3D1-Abs of type 1, type 2, and type 3 CKD to be 0.791, 0.874, and 0.735, respectively ( Fig. 4b-d). Type 2 CKD showed the highest AUC value among all diseases examined.
Association of s-AP3D1-Ab levels with ESCC or CRC . The s-AP3D1-Abs levels were also measured in serum samples from the HDs and patients with ESCC or CRC. AlphaLISA results revealed that s-AP3D1Ab levels were significantly higher in patients with ESCC and CRC than in HDs (Fig. 5a). The positivity rates of s-AP3D1-Abs in HDs and patients with ESCC and CRC were 3.1%, 42.2%, and 15.6%, respectively ( Table 5). The AUC values were 0.872 and 0.743 for ESCC and CRC, respectively (Fig. 5b,c).

Correlation analysis.
Correlation analysis of s-AP3D1-Ab levels and subject data was performed using 633 specimens from Chiba Prefectural Sawara Hospital, including 139 samples from HDs, 121 from patients with DSWMH, 17 from patients with asymptCI, 43 from patients with TIA, 226 from patients with AIS, 57 from patients with cCI, and 30 from other diseases. Baseline characteristics of the study subjects in Sawara Hospital cohort are summarized in Supplementary Table S1. Using the Mann-Whitney U test, the s-AP3D1-Ab levels were compared between participants with body mass index (BMI) < 25 and BMI ≥ 25; participants with or without diseases of DM, hypertension (HT), CVD, and dyslipidemia; and between those patients who were smokers or nonsmokers and those who consumed alcohol or not. The analysis showed that the s-AP3D1-Ab levels were significantly higher in the subjects with HT than in those without HT and those with DM than without DM (Table 6). Conversely, no significant differences in s-AP3D1-Ab levels were observed in the other categories. Performing a logistic regression analysis of the predictors for AIS using the results of the Sawara Hospital was also considered, which included 139 samples from HDs and 226 from patients with AIS. An elevated s-AP3D1-Ab level was associated with an increased risk of AIS as shown by the univariate logistic regression analysis (P < 0.0001). A multivariate logistic regression analysis has identified age, HT, and DM, but not s-AP3D1-Ab, as independent predictors of AIS (Supplementary Table S2).  www.nature.com/scientificreports/  Table 2. ROC curves to assess the ability of s-AP3D1-Abs to predict DM are shown (b). The numbers in the figures are the same as those shown in Fig. 1. The results of the power calculation were shown in Supplementary Table S3. Table 2. Comparing the anti-AP3D1 antibody levels between HDs and patients with diabetes mellitus (DM). The upper panel indicates the number of total samples, samples from male and female participants, and ages (average ± SD). The lower panel summarizes the s-AP3D1-Ab levels examined using AlphaLISA. Numbers are as shown in Table 1; P-values of < 0.05 and positive rates of > 10% are marked in bold font. The plots for these data are shown in Fig. 2b. www.nature.com/scientificreports/ Next, correlation analysis was performed using a Spearman application in order to determine the correlation between s-AP3D1-Ab levels and subject parameters, including general information such as age, body height, weight, BMI, and the degree of artery stenosis [the maximum intima-media thickness (max IMT)]. The serum s-AP3D1-Ab levels were determined to be closely correlated with age (P < 0.001), max IMT (P < 0.001), blood pressure (BP) (P < 0.001), and smoking period (P < 0.001) ( Table 7). Conversely, inverse correlation was observed between s-AP3D1-Ab levels and height, weight, Ca, and low-density lipoprotein cholesterol. Blood glucose and glycated hemoglobin, which are identified as DM markers, were not significantly correlated with  Table 3. ROC curves to assess the ability of s-AP3D1-Abs to predict CVD are shown (b). The numbers in the figures are the same as those shown in Fig. 1. The results of the power calculation were shown in Supplementary Table S3. Table 3. Comparing the anti-AP3D1 antibody levels between HDs and cardiovascular disease (CVD) patients. The upper panel indicates the number of total samples, samples from male and female participants, and ages (average ± SD). The lower panel summarizes the s-AP3D1-Ab levels examined using AlphaLISA. Numbers are as shown in Table 1; P-values of < 0.05 and positive rates of > 10% are marked in bold font. The plots for these data are shown in Fig. 2. www.nature.com/scientificreports/ the s-AP3D1-Abs levels. These results suggest that s-AP3D1-Ab reflected atherosclerosis and its causal HT and smoking.

JPHC cohort analysis.
To determine whether s-AP3D1-Ab marker can be applied to predict the onset of AIS, JPHC cohort samples were examined. The antibody level of AP3D1 protein was positively and strongly associated with the risk of AIS: the ORs (95% CI) were 1.40 (0.75-2.63), 1.97 (1.07-3.65), and 2.28 (1.26-4.13) for the samples with the second, third, and highest quartiles of antibody level, respectively, compared with the lowest quartile (Table 8). These results indicate that the antibody markers against AP3D1 are useful in predicting the onset of AIS.

Discussion
Through the initial SEREX screening, AP3D1 was identified as an antigen recognized by serum IgG in patients with atherosclerosis. The s-AP3D1-Ab levels were higher in patients with AIS, TIA, DM, CVD, CKD, ESCC, and CRC than in the HDs (Figs. 1-5 and Tables 1-5). Meanwhile, to reduce the effect of age, we compared the serum antibody levels of HDs and age-matched patients. The results showed that the s-AP3D1-Ab levels in patients with AIS and DM were significantly higher than those in HDs (Supplementary Fig. S1). Among these diseases, the highest positive rates were observed for ESCC, DM, and type 2 CKD (Tables 1-5). The AUC values  Table 4. The total average values, SDs, cutoff values, positive numbers, positive rates (%), and P-values are summarized in Table 4. The ability of s-AP3D1-Abs to (b) predict type 1, (c) type 2, and (d) type 3 CKD was also evaluated using the ROC curve analysis. The results of the power calculation were shown in Supplementary www.nature.com/scientificreports/ for nephrosclerosis type 2 CKD and ESCC, diabetic type 1 CKD, and DM were 0.874 and 0.872, 0.791, and 0.791, respectively, which were higher than other diseases. The comparison using the Mann-Whitney U test revealed that the s-AP3D1-Ab levels were significantly higher in the subjects with DM than in those without DM (Table 6). In contrast, no significant correlation was found between the s-AP3D1-Ab levels and DM markers, including blood glucose and glycated hemoglobin (Table 7). Consequently, the s-AP3D1-Ab levels do not directly reflect DM, but are associated with DM-induced atherosclerotic disorders, which are also related to CKD and cancer. Consistently, Spearman correlation analysis revealed a significant association between s-AP3D1-Ab levels and max IMT (P < 0.001), which reflects arterial stenosis, namely, atherosclerosis ( Table 7). The antibody levels significantly correlated with HT (P < 0.0001) ( Table 6), which are well-known risk factors for atherosclerosis 34 . A univariate logistic regression analysis revealed that an elevated AP3D1-Ab level was associated with an increased risk of AIS (P < 0.0001). A multivariate logistic regression analysis has also identified age, HT, and DM, but not AP3D1-Ab as independent predictors of AIS (Supplementary Table S2). Therefore, s-AP3D1-Ab marker may discriminate a certain type, if not all, of atherosclerosis caused by HT or DM, leading to the development of AIS and CVD. AP3D1 is a subunit of the AP3 adaptor-like complex 35 , which is expressed in the ubiquitous AP3 complex and also in the neuronal form 36 . AP3 is a heterotetrameric adaptor protein involved in the biogenesis of lysosomerelated organelles, such as platelet-dense bodies. Mouse mutants of the null AP3D1 allele were reported to show abnormal bleeding due to the absence of a storage pool of dense platelet particles, raising the possibility that AP3D1 functions as a thrombogenic regulator through platelet function 37 . Platelets have been identified to play a significant role in hemostatic and thrombotic processes, where abnormal platelet adhesion/activation can lead to the formation of clots (thrombosis) 38 . Thrombosis is known to be closely associated with atherosclerosis 39 . Hirokawa et al. 40 reported that AP3D1-DOT1L-SF3A2 was identified as a new susceptibility locus for myocardial infarction (MI) by European genome-wide association studies (GWAS), which is consistent with our observation that s-AP3D1-Abs levels were significantly higher in patients with CVD including MI (Fig. 3a, Table 3). On the other hand, Xiao et al. reported that angiotensin II facilitated the binding of AP3D1 with beta-arrestin, which resulted in the activation of AP3D1 as a scaffold protein 41 . Angiotensin II plays a key role in the pathogenesis of HT 42-45 , leading to endothelial dysfunction and atherosclerosis 46,47 . AP3D1 can possibly mediate angiotensin Table 4. Comparing the s-AP3D1-Ab levels between HDs and patients with chronic kidney disease (CKD). The numbers shown are as described in Table 1. CKD was divided into three groups as follows: type 1, diabetic kidney disease; type 2, nephrosclerosis; and type 3, glomerulonephritis. P values of < 0.05 and positive rates of > 10% are marked in bold font. The plots for these data are shown in Fig. 4a.
Numbers for the total samples, samples from male and female participants and ages (average ± SD) www.nature.com/scientificreports/  Table 5. ROC curves to assess the ability of s-AP3D1-Abs to predict (b) ESCC and (c) CRC are shown. The numbers in the figures are the same as those shown in Fig. 1. The results of the power calculation were shown in Supplementary Table S3. Table 5. Comparing the serum anti-AP3D1 antibody levels of HDs versus those of patients with esophageal squamous cell carcinoma (ESCC) and colorectal carcinoma (CRC). The s-AP3D1-Ab levels examined using AlphaLISA in HDs and patients with ESCC and CRC are shown. Purified AP3D1-GST proteins were used as antigens. The numbers shown are as described in Table 1. P-values of < 0.05 and positive rates of > 10% are marked in bold font. The plots for these data are shown in Fig. 5a. www.nature.com/scientificreports/ www.nature.com/scientificreports/ www.nature.com/scientificreports/ II-induced HT and atherosclerosis. This was compatible with the results that the s-AP3D1-Ab levels were significantly associated with HT (P < 0.0001) ( Table 6) and BP (P < 0.001) ( Table 7). HT is also known as a risk factor for various types of cancers [48][49][50][51] . Previous reports have proved that angiotensin II is associated with cancer development 52,53 . Angiotensin II is a major upstream regulator of cancer cachexia 54 and can further stimulate angiogenesis and tumor growth of breast cancer 55,56 . Thus, AP3D1 could also mediate angiotensin II-induced carcinogenesis. Based on the results, s-AP3D1-Abs levels were elevated in cancer patients (Fig. 5).

Patient group
The progression of atherosclerosis and cancer often takes several years or more and, in the early stages, is sometimes accompanied by low levels of tissue destruction, which can lead to leakage of proteins from the cells. During this repeated leakage of proteins, even low levels of antigens can induce amplified expression of the antibodies 57 . Thus, antibody markers are deemed more sensitive than antigen markers and may be useful for the early diagnosis of solid cancers including ESCC and CRC. s-AP3D1-Ab marker was closely associated with max IMT, an index of atherosclerosis, which then leads to the onset of AIS and AMI. Thus, predicting these onsets using s-AP3D1-Ab marker is possible. In fact, the results of JPHC cohort analysis indicated that the antibody marker against AP3D1 is useful in predicting the onset of AIS (Table 8). Atherosclerotic AIS and AMI and cancer have been identified as the leading causes of death; thus, the s-AP3D1-Ab marker would be highly useful to reduce its mortality.
As antihypertensive agents, statins, and antiplatelet agents are generally known to prevent the pathogenesis of atherosclerosis 58-60 , the potential modulatory effects of these drugs on s-AP3D1-Ab levels must be considered. Second, as the controls were healthy volunteer donors, potential confounding factors between patients with atherosclerosis, cancers, and controls (age, BMI, HT, DM, and hyperlipidemia) were not adjusted in the analysis of this study. Lastly, physiological testing, such as baPWV, or coronary artery calcification, was not performed to evaluate atherosclerosis in subjects subjected to the s-AP3D1-Ab analysis. Nevertheless, these tests might be expected to confirm the results of this study. The study population included only Japanese patients; thus, further studies are required in patients who are not taking drugs that can affect atherosclerosis and in other ethnic groups. Development of more biomarkers for the early diagnosis of atherosclerotic diseases and early detection of the development of tumors may improve the quality of life.

Data availability
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request. www.nature.com/scientificreports/