- Original investigation
- Open Access
Longitudinal association of C-reactive protein and Haemoglobin A1c over 13 years: the European Prospective Investigation into Cancer - Norfolk study
Cardiovascular Diabetologyvolume 14, Article number: 61 (2015)
Type-2 diabetes is associated with systemic inflammation and higher C-reactive protein (CRP) levels. However, the longitudinal association of CRP and haemoglobin-A1c (HbA1c) has not been described in large prospective studies. Understanding such associations may shed light on the role of inflammation in development of type-2 diabetes and its complications such as cardiovascular diseases.
EPIC-Norfolk is a cohort study of men and women aged 40–79 years at time of recruitment (1993–1997). Serum CRP (mg/l) was measured using a high-sensitivity assay at baseline and 13-years follow-up. HbA1c (%) was measured at baseline, 4, and 13 years. Participants were excluded if they were diagnosed with diabetes or were taking diabetes medication. Data on at least one measurement of CRP and HbA1c was available for 14228 participants (55 % of the cohort).
In the cross-sectional analysis of baseline data, a 1-SD higher loge-CRP (about three-fold higher CRP) was associated with 0.06 (95 % CI 0.04, 0.08) higher HbA1c (%) adjusted for potential confounders. In longitudinal analysis using multivariable linear mixed models, change in CRP over 13 years was to a similar extent positively associated with increase in HbA1c, such that 1-SD higher longitudinal change in loge-CRP was associated with 0.04 (95 % CI 0.02, 0.05) increase in HbA1c.
In this study we found longitudinal observational evidence suggesting that increase in systemic inflammation is associated with an increase in HbA1c and thus systemic inflammation may have a role in development of type-2 diabetes and its complications.
Serum C-reactive protein (CRP) is a marker of systemic inflammation and has been shown to be associated with incident type 2 diabetes [1,2]. Patients with type 2 diabetes are also at risk of developing micro-vascular and macro-vascular complications such as diabetic retinopathy, nephropathy, neuropathy, atherosclerosis, and cardiovascular diseases. The potential etiologic role of chronic systemic inflammation in the pathogenesis of type 2 diabetes and in the link between diabetes and micro-and macro-vascular complications is not well known.
Haemoglobin A1c (HbA1c) concentrations in patients with diabetes are associated with risk of developing micro- and macro-vascular complications, and is thus used as the “gold standard” for monitoring glucose control in patients with type 2 diabetes . HbA1c has also been advocated to be used as a marker for diagnosis of type 2 diabetes . Therefore, understanding the association between markers of systemic inflammation (such as CRP) and haemoglobin A1c may help better understand the potential role of inflammation in the development of type 2 diabetes and in the mechanisms involved in development of complications in patients with diabetes. A positive cross-sectional association between CRP and HbA1c has been reported previously [4-7]. However, the temporal association of CRP and HbA1c over time has not been described. Exploring the longitudinal association between changes in CRP and changes in HbA1c over time in large-scale observational studies should be more informative of the nature of any temporal associations than in cross-sectional studies. In this study we aimed to examine the longitudinal association of CRP and HbA1c in a large prospective cohort study with repeated measurements of CRP, HbA1c, and other confounders, over 13 years of follow up.
Study participants were recruited in the European Prospective Investigation into Cancer in Norfolk (EPIC-Norfolk) study. Details of this study have been previously described . Briefly, 25,639 men and women, aged 40–79, resident in Norfolk, England were recruited to the study using general practice age and sex registers. The age-sex structure of the participants recruited in the study was approximately similar to the UK population age-sex structure according to registers in the United Kingdom. Study participants attended the baseline health examination between 1993 and 1997, a second health examination at around 4 years, and a third health examination around 13 years after recruitment. The EPIC-Norfolk study was approved by the Norfolk Local Research Ethics Committee, and all volunteers gave written informed consent.
At the clinic visits, extensive demographic, medical, lifestyle, family history and dietary data were collected by asking participants to complete a health and lifestyle questionnaire. Smoking status was categorized into current-, former-, and never-smokers, and alcohol consumption was computed as units per week. Medication use (including use of corticosteroids and postmenopausal hormone replacement therapy (HRT)) was assessed by self-report on the health and lifestyle questionnaire and by examining medications brought by the participants to the clinic visit. Habitual physical activity was assessed using the EPIC-validated short physical activity questionnaire by combining levels of occupational and leisure-time physical activity into 4 categories .
Participants were excluded from the present analysis if they were diagnosed with diabetes by a physician at baseline, or were taking diabetes medication. Prevalent and incident cases of type-2 diabetes were either ascertained by self-report of diabetes diagnosed by a physician, as well as by self-report of new administration of diabetes medication at follow up provided that it was verified by medical records, or by record linkage with external sources independent of individuals’ participation in follow up health examinations or postal questionnaires. External sources used to ascertain diabetes cases include listing with general practice diabetes registers, local hospital outpatient diabetes registers, hospital admissions information that screened for diabetes related conditions, and Office for National Statistics mortality data with coding for diabetes.
At the clinic visits, trained nurses took anthropometric measurements on individuals in light clothing without shoes. Height was measured to the nearest millimetre using a free standing stadiometer and weight was measured to the nearest 100 g using digital scales. Body mass index (BMI) was calculated as weight in kilograms divided by the square of height in metres. Waist circumference was measured to the nearest millimetre using a D-loop tape measure.
HbA1c (%) was measured in fresh EDTA blood samples by using high-performance liquid chromatography at baseline and second health examination using the BioRad Diamat Automated Glycosylated Haemoglobin Analyser (Hemel Hempstead, United Kingdom), and at the third health examination using the Tosoh Automated Glycohemoglobin Analyzer.
Blood samples used for measurement of CRP at baseline were centrifuged at 2,100 g for 15 min at 4 °C and then kept frozen in −80 °C freezers until being thawed in 2008 for assaying CRP. Serum high-sensitivity CRP (mg/l) was measured using the Olympus AU640 chemistry analyzer (Olympus Diagnostics, United Kingdom). CRP was not measured at the second health examination, but was measured at the third health examination in fresh blood samples using the Siemens Dimension clinical chemistry analyzer (Siemens Dimension clinical chemistry analyzer, Newark, Delaware, US).
Serum CRP and HbA1c levels were analysed as continuous and categorical variables in separate analyses. Serum CRP had a right skewed distribution and for continuous analysis was log-transformed to obtain a normal distribution. For categorical analyses CRP was divided into 4 clinically relevant categories using cut-points 1, 3, and 10 mg/l.
In cross-sectional analysis using baseline data, multiple linear regression analysis was used to examine mean difference in HbA1c across CRP categories. Analyses were adjusted for age, sex, BMI, waist circumference, smoking, physical activity, alcohol intake, corticosteroid medication, and for women only, menopausal status and HRT use. Missing values for categorical variables were coded as such and thus individuals with missing data were not excluded from the analyses.
In longitudinal analyses, linear mixed models were used to examine the rate of change of HbA1c across baseline CRP categories. The rate of change of HbA1c across categories of CRP was modelled as a fixed-effect interaction between follow up time (years) and CRP categories, allowing for subject-specific random effects for the intercept and time coefficients. Statistical significance of differences in longitudinal rates of change in HbA1c levels by baseline CRP levels was assessed based on P-value for continuous interaction between baseline loge-CRP and time, with subject-specific random effects for the intercept and time. The magnitude of longitudinal change in HbA1c associated with a standard deviation (SD) longitudinal increase in loge-CRP (corresponding to a three-fold increase in original CRP (mg/l) values) over time was calculated by fitting a linear mixed model allowing for random subject specific slopes for loge-CRP and random subject specific intercepts. All covariates were treated as time-varying in the longitudinal analysis to account for changes in life-style and anthropometric factors over time.
Linear mixed models treat missing data as “missing at random”. However, loss to follow-up in our study was likely to be due to higher morbidity and mortality (N = 1,074 (14 %) of participants included in the analysis at baseline died during follow-up), both of which are associated with CRP and HbA1c measures. Therefore, to minimize selection bias, we performed a complete case analysis in which we restricted the longitudinal analyses to participants that had adequate data for analysis at baseline and at least one follow-up health examination. The analyses were also repeated stratified by sex and BMI. Since the associations were similar in sex stratified analyses, all subsequent results are presented as sex-combined. All analyses were performed using Stata 12.0 (StataCorp. College Station, TX).
A total of 14,228 participants who were free of diabetes (as ascertained from self-report or medical records), had data on CRP at baseline, and at least one measurement of HbA1c, were included in the present analysis. The mean age of participants at baseline was 58 years (SD = 9) and 44 % were male (Table 1). The mean (SD) duration of follow up was 3.7 (0.7) and 13.1 (1.8) years at the second and third examinations respectively. Among participants that attended the baseline and follow-up health examinations, median (inter-quartile range) CRP (mg/l) at baseline and 13 years was 1.4 (0.7-3.0) and 2.0 (1.3-3.4) respectively. Mean (standard deviation (SD)) HbA1c (%) at baseline, 4, and 13 years were 5.3 (0.7), 5.5 (0.5), and 5.8 (0.6) respectively.
At baseline assessment, participants with higher serum CRP levels tended to be older, had higher HbA1c levels, waist circumference, and BMI, were less likely to be physically active, more likely to smoke, have a pre-existing history of disease, and more likely to take corticosteroid medication or HRT (Table 1).
In cross-sectional multiple linear regression analyses (Table 2), a 1-SD higher baseline loge-CRP (about three-fold higher CRP on the original scale (mg/L)) was associated with 0.09 (95 % CI 0.07, 0.10) higher HbA1c (%) adjusted for age and sex. The corresponding value was 0.06 (95 % CI 0.04, 0.08) after adjustment for other covariates. Multivariable adjusted mean HbA1c levels were higher in participants with higher baseline CRP levels (Fig. 1). No significant interaction by sex or BMI was observed and the positive associations persisted in stratified analyses (Fig. 1).
HbA1c increased by about 0.05 (SD = 0.05) annually. In longitudinal analyses, the annualised rates of change of HbA1c was significantly associated with baseline CRP levels adjusted for age, sex and follow-up time (P = 0.02 for continuous baseline loge CRP*time interaction, Table 2). Annualized rates of change of HbA1c were no longer significantly associated with baseline CRP levels in multivariable analysis. The results were not changed when analysis was restricted to individuals who had data for CRP, HbA1c, and other covariates at baseline and third health examination (data not shown).
In addition to baseline levels of CRP, change in CRP over time was also significantly and positively associated with longitudinal change in HbA1c (Fig. 2). A 1-SD increment in loge CRP (corresponding to 3-fold increase in CRP on the original mg/l scale) over 13 years of follow-up was associated with an increase of 0.06 (95 % CI 0.04, 0.07) in HbA1c adjusted for age, sex, and follow up time. After adjustment for changes in other covariates in multivariable models, the associations attenuated but remained statistically significant, such that a three-fold increase in CRP over 13 years was associated with 0.04 (95 % CI 0.02, 0.05) increase in HbA1c over the same period.
In this study we found a cross-sectional positive association between serum CRP and HbA1c. A 1-SD higher loge-CRP was cross-sectionally associated with approximately 0.06 higher HbA1c (%). Changes in levels of serum CRP over time were also positively associated with change in HbA1c levels, but to a lesser extent to that found in cross-sectional analysis; such that 1-SD increase in loge-CRP levels over time was associated with ~ 0.04 increase in HbA1c (%) over the same period. These associations were largely independent of age, sex, BMI, waist circumference, smoking status, physical activity, alcohol intake, medical history of diseases, corticosteroid medication, and HRT use.
The cross-sectional positive association between CRP and HbA1c levels has been previously shown in studies involving participants with diabetes  and without diabetes [4, 6, 7] among different populations and age groups. The multivariable adjusted cross-sectional associations of HbA1c and log-CRP observed in the present analysis were similar to those in previous studies [4, 6, 7]. CRP has also been shown to be positively associated with incident type 2 diabetes in previous studies [1, 2]. However, to the best of our knowledge, the longitudinal association between CRP and HbA1c has not been previously reported.
Strengths and limitations
An important consideration in investigating the role of inflammation in development of diabetes and abnormal glucose metabolism is confounding by obesity. Obesity and high body fat are associated with diabetes and abnormal glucose metabolism, as well as with elevated circulating levels of CRP and other pro-inflammatory cytokines [10-12]. A strength of the present study is the large sample size and the availability of repeated measurements of all variables including BMI and waist circumference. Availability of repeated measurements reduces the possibility of measurement error and also enables us to better account for concurrent changes in confounding factors (such as BMI and waist circumference) over time and thus enabling better understanding of the association between inflammation and glucose metabolism independent of obesity.
We excluded participants who were diagnosed with diabetes by a physician or were taking diabetes medication at baseline because treatment for diabetes might affect the association between CRP and HbA1c. We also restricted the analyses to individuals who attended both the baseline and follow up health examinations, and had adequate data required for analysis, to avoid potential bias due to lower survival in those with higher CRP and HbA1c levels.
A limitation of the present study was that CRP at baseline was measured in serum samples that were kept frozen for about 15 years . However, previous studies have shown that CRP is relatively stable in frozen plasma . CRP at the third health examination was measured in fresh serum samples using a different chemistry analyzer. Difference in the methods and assays used for measurement of CRP and HbA1c across different time points may have increased the measurement error. However, the intra-class correlation coefficient between measurements of HbA1c (0.62) was considerably higher than that of systolic blood pressure (0.39) and LDL-cholesterol (0.28), and the intra-class correlation coefficient between the two measurements of log CRP (0.39) was similar to that of systolic blood pressure, and higher than that of LDL cholesterol measured at the same time points in the EPIC-Norfolk study (unpublished data). Similar results were reported in other studies [10, 15].
Although we observed a temporal positive association between change in CRP and change in HbA1c, it is not possible from an observational study to infer causality. Knowledge of biological mechanisms justifies the plausibility of association. Tumour necrosis factor alpha (TNF-α), and interleukin-6 (IL-6) are pro-inflammatory cytokines over expressed in adipose tissue. IL-6 and TNF-α in turn stimulate hepatic production of CRP . The link between obesity inflammation, and abnormal glucose metabolism may in part be related to increased production of pro-inflammatory cytokines such as IL-6 and TNF-α in adipose tissue. TNF-α is suggested to induce insulin resistance . Interleukin-6 and CRP are also shown in animal studies to impair insulin signalling and insulin sensitivity mainly through insulin receptor substrate phosphorylation [18-22]. Further investigation of the association between CRP and markers of insulin resistance such as HOMA-IR is required to underpin this hypothesis.
Higher levels of CRP and leptin were shown to be associated with increased metabolic syndrome components . On the other hand, plasma leptin levels are associated with vascular endothelial function in overweight patients with type 2 diabetes  and leptin was found to be directly related to CRP independently of BMI and other confounding factors especially in men . These findings are consistent with and support a role of CRP and inflammation in abnormal glucose metabolism.
CRP is a marker of an age-related pro-inflammatory process described as “inflamm-aging” [26-28]. According to this theory, long-term exposure to a variety of antigens (infection, food, etc.) results in progressive filling of the immune system by activated immune cells and pro-inflammatory cytokines. This process may lead to a reduced capacity to respond to infectious and stress factors later in life. It may additionally lead to a pro-inflammatory state resulting from increased activated immune cells and pro-inflammatory cytokines. The pro-inflammatory state is suggested to cause long-term tissue damage and an increased susceptibility to non-infectious diseases, such as diabetes and impairments in glucose metabolism, where immunity and inflammation play a role [26-29]. The studies exploring the association between CRP related genes and HbA1c levels or incident diabetes do not support a causal role for CRP . Nevertheless, CRP genotype polymorphisms are associated with a small difference (about 60 %) in serum CRP levels, and of insufficient magnitude for an association with incident diabetes to be detected. Moreover, the negative findings of the studies regarding CRP related genes and diabetes risk do not preclude a causal role for the inflammatory process (of which CRP is an indicator) in development of diabetes.
It remains unclear whether systemic inflammation has a causal role in development of diabetes and related conditions, or whether it is a case of reverse causation, i.e., disturbance in glucose metabolism being accountable for causing systemic inflammation. Previous studies have shown that hyperglycaemia induces secretion of pro-inflammatory cytokines and acute phase reactants by adipocytes and monocytes [31,32]. Moreover, the possibility of the temporal association between CRP and HbA1c being due to residual confounding by unmeasured genetic, lifestyle, or environmental factors cannot be ruled out.
The associations observed in the present study provide observational evidence suggesting that systemic inflammation may potentially have an etiologic role in development of diabetes and may be involved in the link between elevated HbA1c levels and development of diabetic complications such as cardiovascular diseases. Further biological and interventional studies are required to elucidate these hypotheses.
Lee CC, Adler AI, Sandhu MS, Sharp SJ, Forouhi NG, Erqou S, et al. Association of C-reactive protein with type 2 diabetes: prospective analysis and meta-analysis. Diabetologia. 2009;52:1040–7.
Wang X, Bao W, Liu J, Ouyang YY, Wang D, Rong S, et al. Inflammatory markers and risk of type 2 diabetes: a systematic review and meta-analysis. Diabetes Care. 2013;36:166–75.
Lenters-Westra E, Schindhelm RK, Bilo HJ, Slingerland RJ. Haemoglobin A1c: Historical overview and current concepts. Diabetes Res Clin Pract. 2013;99:75–84.
Rogowski O, Shapira I, Peretz H, Berliner S. Glycohaemoglobin as a determinant of increased fibrinogen concentrations and low-grade inflammation in apparently healthy nondiabetic individuals. Clin Endocrinol (Oxf). 2008;68:182–9.
Fukuhara M, Matsumura K, Wakisaka M, Takata Y, Sonoki K, Fujisawa K, et al. Hyperglycemia promotes microinflammation as evaluated by C-reactive protein in the very elderly. Intern Med. 2007;46:207–12.
Temelkova-Kurktschiev T, Siegert G, Bergmann S, Henkel E, Koehler C, Jaross W, et al. Subclinical inflammation is strongly related to insulin resistance but not to impaired insulin secretion in a high risk population for diabetes. Metabolism. 2002;51:743–9.
Wu T, Dorn JP, Donahue RP, Sempos CT, Trevisan M. Associations of serum C-reactive protein with fasting insulin, glucose, and glycosylated hemoglobin: the Third National Health and Nutrition Examination Survey, 1988–1994. Am J Epidemiol. 2002;155:65–71.
Day N, Oakes S, Luben R, Khaw KT, Bingham S, Welch A, et al. EPIC-Norfolk: study design and characteristics of the cohort. European Prospective Investigation of Cancer. Br J Cancer. 1999;80 Suppl 1:95–103. 95–103.
Wareham NJ, Jakes RW, Rennie KL, Schuit J, Mitchell J, Hennings S, et al. Validity and repeatability of a simple index derived from the short physical activity questionnaire used in the European Prospective Investigation into Cancer and Nutrition (EPIC) study. Public Health Nutr. 2003;6:407–13.
Ahmadi-Abhari S, Luben RN, Wareham NJ, Khaw KT. Distribution and determinants of C-reactive protein in the older adult population: European Prospective Investigation into Cancer-Norfolk study. Eur J Clin Invest. 2013;43:899–911.
Arcari A, Zito F, Di CA, De CA, Dirckx C, Arnout J, et al. C reactive protein and its determinants in healthy men and women from European regions at different risk of coronary disease: the IMMIDIET Project. J Thromb Haemost. 2008;6:436–43.
Thorand B, Baumert J, Doring A, Herder C, Kolb H, Rathmann W, et al. Sex differences in the relation of body composition to markers of inflammation. Atherosclerosis. 2006;184:216–24.
Ahmadi-Abhari S, Kaptoge S, Luben RN, Wareham NJ, Khaw KT: Longitudinal Association of C-Reactive Protein and Lung Function Over 13 Years: The EPIC-Norfolk Study. Am J Epidemiol 2013.
Nilsson TK, Boman K, Jansson JH, Thogersen AM, Berggren M, Broberg A, et al. Comparison of soluble thrombomodulin, von Willebrand factor, tPA/PAI-1 complex, and high-sensitivity CRP concentrations in serum, EDTA plasma, citrated plasma, and acidified citrated plasma (Stabilyte) stored at −70 degrees C for 8–11 years. Thromb Res. 2005;116:249–54.
Kaptoge S, Di AE, Lowe G, Pepys MB, Thompson SG, Collins R, et al. C-reactive protein concentration and risk of coronary heart disease, stroke, and mortality: an individual participant meta-analysis. Lancet. 2010;375:132–40.
Marnell L, Mold C, Du Clos TW. C-reactive protein: ligands, receptors and role in inflammation. Clin Immunol. 2005;117:104–11.
Hotamisligil GS. Inflammation and metabolic disorders. Nature. 2006;444:860–7.
Xu JW, Morita I, Ikeda K, Miki T, Yamori Y. C-reactive protein suppresses insulin signaling in endothelial cells: role of spleen tyrosine kinase. Mol Endocrinol. 2007;21:564–73.
Xi L, Xiao C, Bandsma RH, Naples M, Adeli K, Lewis GF. C-reactive protein impairs hepatic insulin sensitivity and insulin signaling in rats: role of mitogen-activated protein kinases. Hepatology. 2011;53:127–35.
D'Alessandris C, Lauro R, Presta I, Sesti G. C-reactive protein induces phosphorylation of insulin receptor substrate-1 on Ser307 and Ser 612 in L6 myocytes, thereby impairing the insulin signalling pathway that promotes glucose transport. Diabetologia. 2007;50:840–9.
Weigert C, Hennige AM, Lehmann R, Brodbeck K, Baumgartner F, Schauble M, et al. Direct cross-talk of interleukin-6 and insulin signal transduction via insulin receptor substrate-1 in skeletal muscle cells. J Biol Chem. 2006;281:7060–7.
Tzeng TF, Liu IM, Cheng JT. Activation of opioid mu-receptors by loperamide to improve interleukin-6-induced inhibition of insulin signals in myoblast C2C12 cells. Diabetologia. 2005;48:1386–92.
Abu-Farha M, Behbehani K, Elkum N: Comprehensive analysis of circulating adipokines and hsCRP association with cardiovascular disease risk factors and metabolic syndrome in Arabs. Cardiovasc Diabetol 2014, 13:76. doi:10.1186/1475-2840-13-76.: 76–13.
Morioka T, Emoto M, Yamazaki Y, Kawano N, Imamura S, Numaguchi R et al.: Leptin is associated with vascular endothelial function in overweight patients with type 2 diabetes. Cardiovasc Diabetol 2014, 13:10. doi:10.1186/1475-2840-13-10.: 10–13.
Nakamura Y, Ueshima H, Okuda N, Miura K, Kita Y, Okamura T et al.: Relation of Serum Leptin and Adiponectin Level to Serum C-Reactive Protein: The INTERLIPID Study. Int J Vasc Med 2013, 2013:601364. doi:10.1155/2013/601364. Epub;%2013 Nov 24.: 601364.
De MM, Franceschi C, Monti D, Ginaldi L. Inflammation markers predicting frailty and mortality in the elderly. Exp Mol Pathol. 2006;80:219–27.
De MM, Franceschi C, Monti D, Ginaldi L. Inflamm-ageing and lifelong antigenic load as major determinants of ageing rate and longevity. FEBS Lett. 2005;579:2035–9.
Goto M. Inflammaging (inflammation + aging): A driving force for human aging based on an evolutionarily antagonistic pleiotropy theory? Biosci Trends. 2008;2:218–30.
Franceschi C, Bonafe M, Valensin S. Human immunosenescence: the prevailing of innate immunity, the failing of clonotypic immunity, and the filling of immunological space. Vaccine. 2000;18:1717–20.
Brunner EJ, Kivimaki M, Witte DR, Lawlor DA, Davey SG, Cooper JA, et al. Inflammation, insulin resistance, and diabetes–Mendelian randomization using CRP haplotypes points upstream. PLoS Med. 2008;5, e155.
Morohoshi M, Fujisawa K, Uchimura I, Numano F. Glucose-dependent interleukin 6 and tumor necrosis factor production by human peripheral blood monocytes in vitro. Diabetes. 1996;45:954–9.
Lin Y, Rajala MW, Berger JP, Moller DE, Barzilai N, Scherer PE. Hyperglycemia-induced production of acute phase reactants in adipose tissue. J Biol Chem. 2001;276:42077–83.
We thank all study participants, general practitioners and the EPIC-Norfolk study team for their contribution. The EPIC-Norfolk study is supported by funding from the Medical Research Council and Cancer Research UK. SA is supported by the Gates Cambridge scholarship. Funding sources did not have a role in the design and conduct of the study, the collection, management, analysis, and interpretation of the data or the preparation, review, approval, or decision to submit the manuscript.
The authors declare that they have no competing interests.
SA analyzed the data and wrote the manuscript with co-authors. SK advised on the statistical analysis. RNL performed all data management and record linkage. KTK and NJW are principal investigators in the EPIC-Norfolk population study. All authors provided detailed comments on the manuscript, reviewed and edited the manuscript and contributed to the discussion. All authors read and approved the final manuscript.