The application of multiple reaction monitoring and multi-analyte profiling to HDL proteins
Lipids in Health and Diseasevolume 13, Article number: 8 (2014)
HDL carries a rich protein cargo and examining HDL protein composition promises to improve our understanding of its functions. Conventional mass spectrometry methods can be lengthy and difficult to extend to large populations. In addition, without prior enrichment of the sample, the ability of these methods to detect low abundance proteins is limited. Our objective was to develop a high-throughput approach to examine HDL protein composition applicable to diabetes and cardiovascular disease (CVD).
We optimized two multiplexed assays to examine HDL proteins using a quantitative immunoassay (Multi-Analyte Profiling- MAP) and mass spectrometric-based quantitative proteomics (Multiple Reaction Monitoring-MRM). We screened HDL proteins using human xMAP (90 protein panel) and MRM (56 protein panel). We extended the application of these two methods to HDL isolated from a group of participants with diabetes and prior cardiovascular events and a group of non-diabetic controls.
We were able to quantitate 69 HDL proteins using MAP and 32 proteins using MRM. For several common proteins, the use of MRM and MAP was highly correlated (p < 0.01). Using MAP, several low abundance proteins implicated in atherosclerosis and inflammation were found on HDL. On the other hand, MRM allowed the examination of several HDL proteins not available by MAP.
MAP and MRM offer a sensitive and high-throughput approach to examine changes in HDL proteins in diabetes and CVD. This approach can be used to measure the presented HDL proteins in large clinical studies.
Recent findings suggest that HDL carries over 80 proteins involved in lipid metabolism, complement regulation, acute phase response and proteinase inhibition . This protein cargo appears to be remodeled in metabolic syndrome , cardiovascular disease [1, 3] and after cholesterol therapies . However, measurement of proteins associated with lipids and that span several orders of magnitude in abundance can be challenging. Conventional LC-MS/MS assays and label-free quantitation methods using spectral counting or extracted ion chromatograms (XIC) used in the previous studies [1–4] can be lengthy and are limited to small sample sizes. Multiple Reaction Monitoring (MRM), on the other hand, is a tandem MS (MS/MS) scan mode unique to triple quadrupole MS instrumentation that is capable of rapid, sensitive, and specific quantitation of peptides in highly complex sample matrices, such as plasma [5, 6]. MRM is a targeted approach that requires knowledge of the molecular weight the peptide of interest and its fragmentation pattern, leading to the generation of target “transitions” for monitoring protein levels. When coupled with stable isotope peptide standards (SIS peptides), quantitation using MRM can be highly reproducible . MRM quantitation has been successfully applied to plasma proteins [5, 6] and, more recently, to a limited subset of HDL proteins . Without pre analytical depletion or enrichment, MRM sensitivity can be compromised by the existence of high abundance proteins. Although the sensitivity of MRM allows analyses at levels down to one ng/mL, the accuracy at the lower end of the abundance range can be problematic. In this case, the use of methods that deplete the higher abundance proteins is often needed. These methods, however, can lead to the unintentional removal of proteins that are attached to the depleted proteins. An alternative approach is the use of immunoassays that offer much improved sensitivity. One example of these assays is multiplexed immunoassay panels (MAP). MAP compliments MRM at this low end and helps ensure that proteins are fully and correctly quantitated at the low range.
There has been considerable interest in understanding HDL functions in light of the strong inverse association of HDL cholesterol and cardiovascular disease (CVD) and the lack of improved CVD outcomes after raising HDL’s cholesterol content in three recent trials [9–11]. Knowledge of HDL protein composition promises to improve our understanding of HDL’s multiple functions in CVD. In this study, we optimized two approaches to measure HDL proteins using multiplexed immunoassays (Multi-Analyte Profiling-MAP) and quantitative proteomics (MRM) platforms. We extended the application of these two techniques to HDL isolated from participants with and without diabetes and CVD. Our findings suggest that MAP can be used to monitor low abundance proteins on HDL, whereas MRM allows examining higher abundance HDL proteins without the need for pre-existing antibodies.
Materials and methods
The study was approved by the University of Arizona Institutional Review Board, and all patients provided written informed consent prior to testing. The samples were collected at University of Arizona Medicine Clinics. Participants reported to the Center for Clinical and Translational Sciences (CaTS) after an overnight fast. Samples were collected in EDTA containing tubes. The presence of cardiovascular disease was defined by a prior history of coronary artery bypass surgery (CABG), percutaneous transluminal angioplasty (PTCA), prior MI, or thrombotic stroke. Patients were older than 18 years of age. New diagnosis of diabetes was based on 2 hour oral glucose tolerance test (OGTT) > 200 mg/dl or glycated hemoglobin (HbA1c > 6.5%). Established diabetes was defined by clinical history. The study included disease free participants and participants with advanced kidney disease, diabetes and CVD. The study excluded subjects if they met any of the following criteria: had type had 1 diabetes, were on an active weight loss program, history of cancer, HIV, or steroid use.
HDL Isolation technique
HDL isolation by centrifugation was based on a modification of a previously published protocol . In brief, potassium bromide (KBr ~55 mg) was added to 310 μl of plasma samples to create a density of 1.21 g/mL. The sample was overlaid with 200 μL of 1.21 g/mL density solution for a total volume was 500 μL. Samples were then spun at 120,000 rpm, at 16°C for 2 hours (Beckman TLX ultracentrifuge with a type 120.1 fixed angle rotor using thick-walled 500 μL Polycarbonate tubes, item 343776). The upper 125 μL solution that had a density of less than 1.21 g/mL was removed and 150 μL of NaCl/EDTA solution (0.9% (w/v) NaCl, 0.1% (w/v) EDTA, pH 7.4) was added to each tube for a final density of 1.063 g/mL. Subsequently, 225 μL of 1.06 KBr solutions in NaCl/EDTA was underlaid with a final volume of 500 μL for a second 2 hour spin at the same parameters listed. The bottom 125 μL (HDL fraction) of solution was removed for further analysis. Four HDL samples were sent to Myriad RBM to externally validate our measurements using an immunoassay in a CLIA certified laboratory. To confirm depletion of albumin and apoB proteins from HDL samples, HDL from a CVD pool was isolated with a second approach that involved long centrifugations at two sequential spins each of 10 hours duration using the above technique. Albumin levels were then measured using a commercial ELISA (Assaypro).
Measurements of HDL proteins
We screened HDL using the human MAP panel (90 proteins) and MRM panel (56 proteins). The proteins that were detected on HDL are summarized in Table 1. Four HDL samples were used to compare protein measures by MAP and MRM. Subsequently, two pools of HDL samples (each pool 500 μL combined from 10 HDL isolates) - one pooled from ten non-diabetic subjects defined as the control group and a second pooled from ten subjects with both diabetes and CVD–defined as the disease group were run on the MAP platform using the HumanMAP panel and by MRM.
The samples were analyzed at Myriad RBM that uses Luminex xMAP. Luminex xMAP is a well-established particle array system that is based on beads with unique fluorescent signatures with proprietary multi-analyte panel targets assessed for cross-reactivity. This technology has been used for the detection of numerous targets, such as cytokines , cancer markers , and indicators for various disease states . The assay details of this platform are well described in Myriad RBM website (http://www.myriadrbm.com). Here we applied MAP to HDL samples obtained by ultracentrifugation in both 2X2 and 10X10 isolation techniques, calculating the limit of detection and the recovery analysis after 12 dilutions. The data is presented in Additional file 1: Table S1. MAP was then applied to the control and disease HDL sample pool.
The samples were analyzed at the University of Victoria - Genome BC Proteomics Centre with a dedicated core service for MRMs with the capacity of high throughput stable isotope peptide production for absolute quantification. We screened HDL protein using existing published transitions that were previously applied to plasma proteins  or to newly developed transitions from proteins that are associated with HDL (LCAT, CETP, PLTP, PON1, Apolipoprotein D, Apolipoprotein M, Apolipoprotein L1, Apolipoprotein CII) or of interest to CVD (Macrophage migration inhibitory factor, Matrix Gla protein). The transition list used is provided in Additional file 1: Table S2. For the generation of CV data, the samples were injected four times per method, with two methods in total. The total number of transitions per sample was 142 (method 1-low abundant specific targets) + 88 (method 2-general high abundance plasma protein targets). There were between 1–5 transitions monitored for each peptide. Proteins with new transitions, or those that were low in abundance were monitored using either multiple peptides (as CETP, MIF, PON1) or up to 5 transitions per peptide. The methods were not scheduled as the retention times shifted between HDL isolation procedures that can result in non-quantifiable data. Based on these replicate runs, we selected one transition for quantitation analysis based on the lowest coefficient of variation by the 4 replicate “technical” runs. These transitions are summarized in Table 2. The selection process of these transitions (to exclude interferences or amino acid modifications such as phosphorylation and glycosylation) was previously described [6, 15] in detail and included in the (Additional file 1: MRM methods). Representative chromatograms of the healthy pool HDL transitions and disease pool HDL transitions are also included in the figure Additional file 1: Figure S1.
We used the statistical program R2.1 (R core development team). The characteristics of the two pools were compared using an independent sample test or a proportion test. The correlation coefficient and the p value generated between plasma proteins and HDL proteins was obtained using spearman correlation test.
Results and discussion
Description of study participants
The goal of this project was to use sensitive and high throughout approach to analyze HDL proteins in health and vascular disease such as diabetes and CVD. The samples were selected to detect differences in HDL proteins between healthy controls and patients with chronic vascular disease. The study demographics and biochemical measurements are listed in Table 3. The two groups were age and sex matched. All diseased individuals had a history of CVD events prior to participation. As shown in Table 3, diseased subjects were obese, with lower HDL cholesterol, elevated triglyceride levels, uncontrolled diabetes, elevated inflammation (CRP), and evidence of chronic kidney disease as revealed by the elevated plasma creatinine.
HDL and Plasma proteins
One challenge in HDL proteomics is differentiating whether proteins are actually present within HDL or are predominately carried within plasma and have become loosely associated with HDL before or during the isolation procedure. We initially isolated HDL using a two sequential spins (2 hrs × 2) and we are able to detect several low abundance proteins commonly associated with atherogenesis and inflammation using MAP (Table 1) of participants with both diabetes and CVD. To confirm whether these proteins are part of the HDL fraction or loosely attached plasma proteins, a longer centrifugation process was employed to isolate HDL (10 hrs × 2). Plasma albumin levels were determined to be 45 mg/mL. After the (2 hrs × 2) centrifugation, albumin levels within the HDL fraction were 0.1 mg/mL. After prolonged sequential centrifugation step (10 hrs × 2), albumin levels were below the detection limits of this assay. Proteins detectable on HDL (by both centrifugation techniques) assessed by MAP, along with plasma protein concentrations are summarized in Table 4. We then correlated the concentrations of these proteins between the different fractions. We found no correlation between plasma and HDL proteins identified in either the (2 hrs ×2) or (10 hrs ×2) HDL isolation preparations. For example, Apo A-I concentrations of 0.36, 0.27, 0.18 mg/mL in plasma HDL (2 hrs × 2) and HDL (10 hrs × 2), respectively shows clear retention on HDL compared to fibrinogen with concentrations of 4.7, 0.0021, 0.0012 mg/mL. These findings confirm that non-HDL proteins were efficiently depleted with the longer centrifugation. In contrast, the HDL proteins isolated with the (2 hrsx2) and (10 hrs × 10) methods were highly correlated (Figure 1, r = 0.95, p < 0.001). As expected the concentrations of HDL proteins were greater in the (2 hrs × 2) than the (10 hrs × 2) technique (Table 3). The concentration of Apo A-I in plasma was low in this sample pool from diseased participants compared to healthy controls. Published Apo A -I levels are in the 0.9-1.6 mg/mL range. This may suggest significant HDL remodeling in advanced vascular disease.
Correlation of MAP and MRM
MRM’s performance for low abundance proteins is limited without prior enrichment. Thus, MAP allows the measurements of low abundance proteins on HDL. There were 10 proteins common to both MAP and MRM, however, one protein (Apo E) was not detected in our samples in 3 out of 4 samples by analyzed by MRM. We correlated 9 proteins measured by both approaches in 4 HDL samples. The results are summarized in Table 5. We did not observe a good correlation for lipoprotein (a), Apo D and Apo CIII between the two methods. Most of Lipoprotein (a) is associated with LDL and present in low abundance on HDL  reflecting a challenge for measuring this protein on HDL by MRM. Apo D on the other hand, has a high degree of homology to retinol binding protein or other members of the alpha 2 microglobulin protein superfamily . This might potentially pose a challenge for the MAP technique. Apo CIII was another protein that was measured by MAP and MRM with a weak correlation between the two measurement platforms. It is not clear why these assays were discordant. It is possible that the peptide used is not a good surrogate of Apo CIII levels. Given that Apo CIII exists in multiple isoforms , both assays could be measuring different forms.
HDL Proteome in diabetes and CVD
There is strong inverse association of HDL cholesterol and cardiovascular disease (CVD). However, recent studies suggested that niacin or CETP inhibition designed to raise HDL-C content did not lead to improved CVD outcomes [9–11]. This discrepancy highlights the need to revise our view of HDL and its functions. One previous study demonstrated that HDL acquires an inflammatory phenotype in acute coronary syndrome . Our results suggest that examining HDL by MAP and MRM can reveal important insights into its biology. One significant finding in our study is the ability to detect for the first time important proposed mediators of atherosclerosis (ENRAGE, MPO, and PAI-1) on HDL using MAP. We also optimized an MRM transition library for the high throughput monitoring of HDL peptides. To investigate the HDL proteome in CVD, HDL-containing samples from 10 non-diabetic participants and HDL from subjects with both diabetes and prior CVD events were pooled into a control and a DM/CVD sample pool. 69 out of the 90 humanMAP proteins showed acceptable recovery in response to multiple dilutions (Additional file 1: Table S1) and were above the lower limit of detection in both sample pools. Figure 2 represents the ratio of the 69 proteins on HDL between disease and control participants. Our findings suggested that several atherosclerosis-associated and inflammatory proteins were increased on HDL of diseased individuals. In addition, MPO, PAI-1, IL1beta, and ENRAGE were only detected on HDL of diseased individuals. We acknowledge that by pooling the samples, the clinical utility of this approach is limited. A larger study of subjects with and without disease is needed to fully characterize the distribution of these low abundance proteins on HDL. Unlike MAP, the sensitivity of MRM is limited without prior sample enrichment. Hence, the accuracy of the assay is compromised when sample concentrations fall below the linear range of the assay as evidenced by the increased coefficient of variation (CV) in several MRM proteins. By MRM, we targeted transitions of 56 proteins. We were able to detect 32 proteins of which 11 proteins had CVs less than 20% on the replicate runs (Tables 1 and 2). Clusterin is a high abundant protein that is associated with HDL. We were able to quantitate clusterin accurately in plasma (CV < 5% ); however we were not able to reliably measure it in HDL as was previously demonstrated  reflecting the challenges of using MRMs for protein quantification when levels of a given target falls out of the dynamic range in diseased states. Comparing HDL proteins in the CVD versus control pool by MRM revealed the depletion of several lipid metabolism proteins such as Apo CI, CII and PON-1 in CVD. In contrast, the concentrations of several acute phase proteins such as clusterin, Complement C9, alpha-1-antitrypsin were increased in CVD (Figure 3). These findings would facilitate conducting larger studies using the presented approach to examine the effect of inflammation and atherosclerosis associated proteins on HDL composition and function in individuals with diabetes and CVD.
Challenges of measuring lipid associated proteins
Measurement of lipid associated proteins has been successful using immunoassays based on nephelometry and mass spectrometry [1, 4, 8]. Nephelometry is not widely available and not amenable to multiplexing. Conventional mass spectrometry techniques using spectral counting or extracted ion chromatograms can be lengthy and challenging in large sample sizes. Here, we present two alternatives (MAP and MRM) that are amenable to multiplexing, and are high throughput. MAP is limited by the availability and quality of pre-exiting antibody panels. MRM, on the other hand, can lose sensitivity at the lower end of abundance. Although the performance of MRM in lower abundance HDL proteins was inferior to immune based assays, the performance of MRM is likely to improve with the development of more sensitive mass spectrometry analyzers and better sample fractionation methods.
There are several limitations of this study. First, the present study was done with two pools of clinical samples. As such, the present work describes an alternative analytical tool for the measurement of HDL proteins. Application of these assays to larger clinical data sets is feasible and will allow determination of their clinical utility. Second, the digestion efficiency of each protein monitored by MRM was not assessed. However, we have previously shown  that if the digestion procedure is consistent and reproducible, then the ability to compare quantitative values between samples is maintained. Thus, the results presented are better described as “relative accurate abundances”. Despite this limitation, several peptides quantitated by MRM were correlated with measurement using MAP (Table 5) affirming our claim that the MRM assay can provide absolute quantitation. A third limitation of this study is that we did not use an independent method of HDL isolation (such as reciprocal co-immunoprecipitation of a given protein with Apo A-I, or gel filtration as suggested by Davidson et al. ) to validate that presence of these low abundance proteins on HDL. However, HDL isolated by a longer centrifugation step had successfully depleted albumin with several of the low abundance proteins still detectable on HDL. In addition, the proteins on HDL and the corresponding plasma concentrations did not correlate, whereas the two HDL fractions were highly correlated. This finding suggests that these low abundance proteins were not contaminant plasma proteins. Our findings however, need to be replicated in a larger study group.
This study suggests the feasibility of measuring HDL proteins using MRM and MAP. The application of MAP and MRM to the HDL proteome offers the potential to improve our understanding of HDL functions and help direct interventions aimed at remodeling the HDL phenotype in diabetes and CVD.
High density lipoprotein
- Apo A-I:
Multiple reaction monitoring
Vaisar T, Pennathur S, Green PS, Gharib SA, Hoofnagle AN, Cheung MC, Byun J, Vuletic S, Kassim S, Singh P: Shotgun proteomics implicates protease inhibition and complement activation in the antiinflammatory properties of HDL. J Clin Invest. 2007, 117: 746-756. 10.1172/JCI26206
Hoofnagle AN, Wu M, Gosmanova AK, Becker JO, Wijsman EM, Brunzell JD, Kahn SE, Knopp RH, Lyons TJ, Heinecke JW: Low clusterin levels in high-density lipoprotein associate with insulin resistance, obesity, and dyslipoproteinemia. Arterioscler Thromb Vasc Biol. 2010, 30: 2528-2534. 10.1161/ATVBAHA.110.212894
Alwaili K, Bailey D, Awan Z, Bailey SD, Ruel I, Hafiane A, Krimbou L, Laboissiere S, Genest J: The HDL proteome in acute coronary syndromes shifts to an inflammatory profile. Biochim Biophys Acta. 1821, 2012: 405-415.
Green PS, Vaisar T, Pennathur S, Kulstad JJ, Moore AB, Marcovina S, Brunzell J, Knopp RH, Zhao XQ, Heinecke JW: Combined statin and niacin therapy remodels the high-density lipoprotein proteome. Circulation. 2008, 118: 1259-1267. 10.1161/CIRCULATIONAHA.108.770669
Anderson L, Hunter CL: Quantitative mass spectrometric multiple reaction monitoring assays for major plasma proteins. Mol Cell Proteomics. 2006, 5: 573-588.
Kuzyk MA, Smith D, Yang J, Cross TJ, Jackson AM, Hardie DB, Anderson NL, Borchers CH: Multiple reaction monitoring-based, multiplexed, absolute quantitation of 45 proteins in human plasma. Mol Cell Proteomics. 2009, 8: 1860-1877. 10.1074/mcp.M800540-MCP200
Addona TA, Abbatiello SE, Schilling B, Skates SJ, Mani DR, Bunk DM, Spiegelman CH, Zimmerman LJ, Ham AJ, Keshishian H: Multi-site assessment of the precision and reproducibility of multiple reaction monitoring-based measurements of proteins in plasma. Nat Biotechnol. 2009, 27: 633-641. 10.1038/nbt.1546
Hoofnagle AN, Becker JO, Oda MN, Cavigiolio G, Mayer P, Vaisar T: Multiple-reaction monitoring-mass spectrometric assays can accurately measure the relative protein abundance in complex mixtures. Clin Chem. 2012, 58: 777-781. 10.1373/clinchem.2011.173856
Boden WE, Probstfield JL, Anderson T, Chaitman BR, Desvignes-Nickens P, Koprowicz K, McBride R, Teo K, Weintraub W: Niacin in patients with low HDL cholesterol levels receiving intensive statin therapy. N Engl J Med. 2011, 365: 2255-2267.
Nissen SE, Tardif JC, Nicholls SJ, Revkin JH, Shear CL, Duggan WT, Ruzyllo W, Bachinsky WB, Lasala GP, Tuzcu EM: Effect of torcetrapib on the progression of coronary atherosclerosis. N Engl J Med. 2007, 356: 1304-1316. 10.1056/NEJMoa070635
Schwartz GG, Olsson AG, Abt M, Ballantyne CM, Barter PJ, Brumm J, Chaitman BR, Holme IM, Kallend D, Leiter LA: Effects of dalcetrapib in patients with a recent acute coronary syndrome. N Engl J Med. 2012, 367: 2089-2099. 10.1056/NEJMoa1206797
Carson RT, Vignali DA: Simultaneous quantitation of 15 cytokines using a multiplexed flow cytometric assay. J Immunol Methods. 1999, 227: 41-52. 10.1016/S0022-1759(99)00069-1
Arellano‒Garcia M, Hu S, Wang J, Henson B, Zhou H, Chia D, Wong D: Multiplexed immunobead‒based assay for detection of oral cancer protein biomarkers in saliva. Oral diseases. 2008, 14: 705-712. 10.1111/j.1601-0825.2008.01488.x
Morgan E, Varro R, Sepulveda H, Ember JA, Apgar J, Wilson J, Lowe L, Chen R, Shivraj L, Agadir A: Cytometric bead array: a multiplexed assay platform with applications in various areas of biology. Clin Immunol. 2004, 110: 252-266. 10.1016/j.clim.2003.11.017
Yassine H, Borges CR, Schaab MR, Billheimer D, Stump C, Reaven P, Lau SS, Nelson R: Mass spectrometric immunoassay and MRM as targeted MS-based quantitative approaches in biomarker development: Potential applications to cardiovascular disease and diabetes. Proteomics Clin Appl. 2013, 7: 528-540. 10.1002/prca.201200028
Scanu AM, Fless GM: Lipoprotein (a). Heterogeneity and biological relevance. J Clin Investig. 1990, 85: 1709- 10.1172/JCI114625
DRAYNA DT, McLEAN JW, WION KL, TRENT JM, DRABKIN HA, LAWN RM: Human Apolipoprotein D Gene: Gene Sequence, Chromosome Localization, and Homology to the α2u-Globulin Superfamily. Dna. 1987, 6: 199-204. 10.1089/dna.1987.6.199
Niederkofler EE, Tubbs KA, Kiernan UA, Nedelkov D, Nelson RW: Novel mass spectrometric immunoassays for the rapid structural characterization of plasma apolipoproteins. J Lipid Res. 2003, 44: 630-639. 10.1194/jlr.D200034-JLR200
Proc JL, Kuzyk MA, Hardie DB, Yang J, Smith DS, Jackson AM, Parker CE, Borchers CH: A quantitative study of the effects of chaotropic agents, surfactants, and solvents on the digestion efficiency of human plasma proteins by trypsin. J Proteome Res. 2010, 9: 5422-5437. 10.1021/pr100656u
Davidson WS, Silva RA, Chantepie S, Lagor WR, Chapman MJ, Kontush A: Proteomic analysis of defined HDL subpopulations reveals particle-specific protein clusters: relevance to antioxidative function. Arterioscler Thromb Vasc Biol. 2009, 29: 870-876. 10.1161/ATVBAHA.109.186031
Sources of funding
Dr. Yassine was supported by K23HL107389, AHA12CRP11750017 and USC CTSI pilot UL1TR000130. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health. MRM mass spectrometric data were acquired by the Arizona Proteomics Consortium supported by NIEHS grant P30ES06694 to the Southwest Environmental Health Sciences Center (SWEHSC to Dr. Lau), NIH/NCI grant P30CA023074 to the Arizona Cancer Center (AZCC), and by the BIO5 Institute of the University of Arizona.
We would like to thank Genome Canada and Genome BC for Science and Technology Innovation Centre funding support. We thank Carol Parker for her assistance with writing the manuscript. We would also like to recognize Tyra J. Cross and Suping Zhang for the synthesis of all SIS peptides and Juncong Yang for exemplary technical support. We also thank Dr. George Tsaprailis with his assistance in running MRMs at the Arizona Proteomics Consortium. Finally, we would like to thank Jeffrey Freiser for his assistance with Myriad RBM HDL runs.
The author’s declare that they have no competing interests.
Participated in research design: HY, SL, DS, CB, DB; conducted experiments: HY, HK, AJ; performed data analysis and interpretation: HY, DB; contributed to the writing of the manuscript: PR, CB; critically revised the manuscript: PR; all authors read and approved the final manuscript.