Background Many serum biochemical indices have been found to be associated with coronary heart disease (CHD); however, few studies have evaluated the value on screening CHD of the integrated serum biochemical indices.
Methods In this study, 627 healthy controls and 1049 patients with CHD were recruited to develop CHD screening models for males and females using unconditional logistic regression. The performance of the screening models was evaluated by areas under the receiver operating characteristic (ROC) curves (AUCs), and externally validated in another population comprised of 190 healthy controls and 246 patients with CHD.
Results Backward stepwise variable selection showed that increasing age, total cholesterol (TC), logarithm-transformed homocysteine (lnHCY), logarithm-transformed γ-glutamyl transpeptidase (lnGGT), and decreasing uric acid, logarithm-transformed triglyceride, apolipoprotein A (apoA) and apolipoprotein B (apoB), increased the detection of CHD in males. In comparison, increasing age, TC, lnHCY, lnGGT and high-density lipoprotein cholesterol versus low-density lipoprotein cholesterol, and decreasing apoA, apoB, logarithm transformed lipoprotein (A) and logarithm transformed total bilirubin, increased the detection of CHD in females. The AUCs for the screening models for males and females were 0.958 (95% CI 0.946 to 0.969) and 0.986 (95% CI 0.977 to 0.994), respectively. The performance of the screening models was further evaluated in external validation samples, the AUCs for males and females were 0.907 and 0.992, respectively.
Conclusions Our study suggests that integrated serum biochemical indices may be used to screen for suspected CHD in participants.
- Cardiovascular disease
- CHD/CORONORARY HEART
- Clinical epidemiology
This is an Open Access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/
Statistics from Altmetric.com
Coronary heart disease (CHD) remains a persistent public health burden worldwide and is the second leading cause of cardiovascular death in the Chinese population.1 ,2 In CHD, the supply of blood, oxygen and nutrients to the heart muscle is reduced because of the build-up of a waxy substance, known as a plaque, (ie, atherosclerosis), in the coronary arteries.3 CHD develops when the coronary arteries become damaged or diseased; over time, plaques can harden or rupture, narrowing the coronary arteries due to the formation of blood clots, and reducing or blocking the flow of blood to the heart.4 ,5 The decrease in blood flow may cause angina, shortness of breath, or other signs and symptoms of CHD.6 Complete blockage can cause a myocardial infarction or even sudden cardiac death.7
Atherosclerosis is now recognised as a complex inflammatory disorder of the coronary arteries.8 When the arterial endothelium encounters certain endovascular substances such as proinflammatory cytokines, the products of glycoxidation, dyslipidaemia or vasoconstrictor hormones, endothelial cells augment the expression of adhesion molecules that promote the adherence of blood leucocytes to the inner surface of the arterial wall.9 ,10 Adherent mononuclear phagocytes and T lymphocytes communicate with the cells of the arterial wall endothelium and smooth muscle cells (SMCs), and trigger an exchange of messages among the various cell types, leading to inflammation in the early atheroma.6 ,11 Lipid mediators such as prostanoids, leucotrienes, histamine, cytokines and complement components participate in this pathophysiological process.6 Subsequently, SMCs migrate from the tunica media to the intima.12 These cells proliferate and elaborate a rich and complex extracellular matrix in response to various oxidative, haemodynamic, inflammatory and autoimmune signals.6 Certain constituents of the extracellular matrix (eg, proteoglycans) bind lipoproteins, prolonging their residence in the intima, and rendering them more susceptible to oxidative modification and glycation.13 These products of lipoprotein modification, including oxidised phospholipids and advanced glycation end products, sustain and propagate the inflammatory response.14 As the lesion progresses, calcification may occur.15 In addition to proliferation, apoptotic cell death occurs in the established atherosclerotic lesion.16 The death of lipid–laden macrophages leads to the extracellular deposition of tissue factor, some in particulate form.6 ,16 The accumulating extracellular lipids in the intima can coalesce, forming the core of the atherosclerotic plaque.6
Epidemiological studies have identified a series of risk factors associated with CHD.17 Age, sex, family history and ethnicity are classic unmodifiable risk factors;18 of these, age is the strongest predictor of adverse cardiovascular outcomes.19 Although cardiovascular disease is the biggest cause of mortality in women, men have about twice the total morbidity and mortality, compared to women.20 Family history is a significant contributing factor to CHD.21 Ethnicity is a well-recognised risk factor, with a higher prevalence of cardiovascular disease in South Asian and black populations.22 Hypertension, dyslipidaemia and diabetes are classic modifiable risk factors.17 Systolic and diastolic blood pressures are strongly positively associated with cardiovascular disease; patients with CHD can benefit from antihypertensive therapy.23 ,24 Cholesterol and triglyceride abnormalities constitute almost half the population-attributable risk.25 Patients with diabetes mellitus are considered to be at high risk of CHD regardless of other factors;26 increasing glycaemic levels are positively correlated with vascular complications.27 Smoking, alcohol abuse, obesity and high dietary intake of high-fat, high-carbohydrate and high-calorie diets are behavioural risk factors associated with CHD.17 Recently, a number of new biomarkers involved in the pathophysiology of CHD have been described.28 Genome-wide association studies have also identified a number of genetic cardiovascular risk biomarkers.17
Awareness of the importance of blood biochemical markers is growing.29 CHD may change the estimated glomerular filtration rate derived from serum creatinine measurement.30 Other blood biochemical markers, including homocysteine (HCY),31 uric acid (UA),32 microalbuminuria and cystatin C,33 ,34 are considered to change in patients with CHD as well; however, they are confounded by the relationship with other risk factors such as smoking.17 As research on CHD has progressed, scientists have started to combine the identified markers in order to produce an index with high screening value.17 Specialist societies recommend initiating preventive treatment of cardiovascular disease on the basis of a person's 10-year risk for CHD. Well-known prognostic models to estimate this risk originate from the Framingham Heart Study,35 the Women's Health Study,36 the PROCAM (Prospective Cardiovascular Münster) study (37) and the SCORE (Systematic Coronary Risk Evaluation) project (38). These risk-estimation systems were derived from large-scale prospective cohorts and focused on the primary prevention of CHD worldwide.35–38 In this study, we aimed to develop a model utilising biochemical markers to screen for suspected CHD in patients with symptoms.
This study was conducted according to the principles of the World Medical Association Declaration of Helsinki. The study was approved by the Review Board of the People's Hospital of Yuxi City (Yuxi, China) (approval number: YNYXH2010-0012). All participants provided written informed consent prior to participating in this study. The participants were simultaneously informed of their right to repeal consent by themselves or their kin, caretakers or guardians.
Model developing population
Patients with CHD were consecutively recruited from Department of Cardiology, People's Hospital of Yuxi City, between September 2010 and July 2013. Patients suspected of having CHD were initially screened based on their symptoms, medical and family histories, risk factors for CHD and a physical examination. Patients’ symptoms included (1) heaviness, tightness, pressure and/or pain in the chest behind the breastbone; (2) pain radiating in the arms, shoulders, jaw, neck and/or back; (3) shortness of breath; (4) weakness and fatigue. Patients suspected of having CHD were further diagnosed according to American Heart Association guidelines. All patients with CHD were confirmed to have arteriostenosis (>50%) of at least one large epicardial coronary artery, caused by an atheromatous plaque, using coronary angiography.
Unrelated consecutive inpatients, in the department of cardiology for other reasons, were recruited as controls during the same time period. These individuals had no history of CHD and were further evaluated by two or more cardiological tests to exclude latent CHD.
Participants who met the following exclusion criteria were excluded from the study: alcohol abuse; diabetes; history of smoking; evidence of non-coronary atherosclerotic disease; a chronic infectious disease, chronic or acute liver disease, chronic lung disease, chronic kidney disease, chronic wasting disease, malignant tumour, autoimmune disease or xanthelasma.
External validation population
To test the value of the screening strategy derived from the above modelling population, another group of unrelated consecutive inpatients were recruited from the same department of the hospital in the following year. Inclusion and exclusion criteria for participants with CHD and controls were the same as those of the model deviation population.
Laboratory-based biochemical testing and questionnaire
A fasting blood sample was collected from each participant via the antecubital vein, in the morning. Serum lipid profiles, including total cholesterol (TC), triglycerides (TG), high-density lipoprotein-C (HDL-C), low-density lipoprotein-C (LDL-C), apolipoprotein (apo)A, apoB, apoE and lipoprotein (a) Lp(a); indicators of liver function, including γ-glutamyl transpeptidase (GGT), total bilirubin (TBil), indirect bilirubin (IBil) and direct bilirubin (DBil); indicators of kidney function, including UA; and HCY levels, were measured in the Hospital's laboratory according to routine procedures.
A questionnaire was self-administered by all participants; the items included age, sex, ethnicity, medical history, family history and manifestations of CHD, alongside behavioural risk factors such as smoking, alcohol abuse, obesity and high dietary intake of high-fat, high-carbohydrate and high-calorie diets.
All analyses were performed using SPSS for Windows (V.18). Continuous variables that were normally or approximately normally distributed are presented as the mean±SD, while those with a skewed distribution are presented as the median (1st and 3rd quantiles). Univariable analyses were performed using independent t tests or Wilcoxon rank tests, as appropriate. To eliminate multicollinearity between or among variables, collinearity diagnostics were performed before further statistical analysis. Continuous variables that were not normally distributed were natural logarithm-transformed prior to the analysis.
Multiple unconditional logistic regression models were developed for risk estimation. Potentially relevant clinical predictors according to the univariable analysis (p<0.10) were introduced into the starting model and then eliminated manually using the backward step-by-step approach, depending on the largest p value. All items showing statistical significance at p<0.05 were retained in the final screening model. Internal cross-validation was used to minimise overfitting. The performance of discrimination was evaluated by an area under receiver operating characteristic (ROC) curve (AUC). The screening model was externally validated in the validation population, with AUC calculated. In view of the discrepancies between males and females in previous risk score models, we developed and validated our models by sex.
Participants included in the model development population
Data collected from the participants between 1 October 2010, and 31 July 2013, were used to develop the screening model. The characteristics of the participants used for model development are listed in table 1.
The model development population included 1049 patients with CHD, with 664 males and 385 females, aged 63.3±11.2 and 62.5±9.0 years, respectively (table 1). The model development population also included 627 controls, with 400 males and 227 females, who had no history of CHD and were further evaluated by two or more cardiological tests to exclude latent CHD (table 1). The mean age (in years) of the control participants was 45.5±12.1 for males and 43.3±10.4 for females (table 1). The laboratory test variables were not always normally distributed (eg, TG, Lpa, HCY, TBil, DBil, IBil and GGT); these are presented as medians (1st and 3rd quantiles) rather than means±SD (table 1).
Deprivation of the screening models
Sex-specific univariate analyses were performed on all available data, to primarily identify variables associated with an increased detection of CHD. These variables included age, UA, TC, logarithm-transformed TG (lnTG), HDL-C, LDL-C, apoA, apoB, logarithm-transformed Lp(a) (lnLPa), logarithm-transformed HCY (lnHCY), logarithm-transformed TBil (lnTBil), logarithm-transformed DBIL (lnDBil), logarithm-transformed IBil (lnIBil) and logarithm-transformed GGT (lnGGT). Our univariate analysis results are given in table 2. Of the male population, age, UA, TC, lnTG, HDL-C, LDL-C, apoA, apoB, lnLp, lnHCY, lnTBil, lnDBil and lnIBil, are the variables associated with an increased diagnosis of CHD (table 2). Of the female population, age, UA, TC, lnTG, HDL-C, LDL-C, apoA, apoB, lnLp, lnHCY, lnTBil, lnIBil and lnGGT, are the variables associated with an increased diagnosis of CHD (table 2). At this stage, the main differences between males and females are variables lnGGT and lnDBil, lnGGT is significantly associated with CHD in females but not males and, conversely, lnDBil is significantly associated with CHD in males but not females (table 2). Variables with p<0.10 were included in a subsequent multiple logistic regression.
In view of the strong multicollinearity among three serum bilirubin indices (lnTBil, lnDBil and lnIBil), we used lnTBil exclusively to represent the serum bilirubin level and introduced the value into the multiple logistic regression model. There was also strong multicollinearity among three serum cholesterol indices (TC, HDL-C and LDL-C); thus, we computed the HDL-C and LDL-C ratio (HDL-C vs LDL-C) as a composite indicator for the multiple regression analysis.
Age, UA, TC, lnTG and HDL-C versus LDL-C, apoA, apoB, lnLp(a), lnHCY and lnTBil, were used to build the male model, while age, UA, TC, lnTG and HDL-C versus LDL-C, apoA, apoB, lnLp(a), lnHCY, lnTBil and lnGGT, were used to build the female model. The results of our backward stepwise selection process are listed in table 3; we found that increasing age, TC, lnHCY and lnGGT, and decreasing UA, lnTG, apoA and apoB, increased the diagnosis of CHD in males. In comparison, increasing age, TC, lnHCY, lnGGT and HDL-C versus LDL-C and decreasing apoA, apoB, lnLp(a) and lnTBil, increased the diagnosis of CHD in females.
Probabilities of diagnosis of CHD for males and females were as following:and
Discrimination of the screening models
We evaluated the performance of the prediction models by drawing ROC curves of the CHD probabilities predicted by the above models and calculating the AUC. The AUC was 0.958 (95% CI 0.946 to 0.969) for males and 0.986 (95% CI 0.977 to 0.994) for females (figure 1, upper). The ROC curves showed that the sensitivity and specificity of the models were near optimal for males and females.
Subjects included in the external validation population
The data used for model validation were collected the following year; namely, 1 August 2013, to 31 July 2014. The model validation population included 246 patients with CHD, with 176 males and 70 females, aged 61.8±11.1 and 64.5±9.3 years, respectively (table 4). The mean age (in years) of the control participants was 52.5±9.2 for males and 51.0±7.4 for females.
External validation of the screening models
To validate the above screening models, a separate population, identified the following year, was used. As described above, the model validation participants included 176 males and 70 females with CHD confirmed by coronary angiography (table 4). We obtained each participant's probability of diagnosing with CHD using the above screening models. Next, ROC curves were plotted for both males and females (figure 1, Lower). The AUC for males (0.907, 95% CI 0.870 to 0.943) declined very slightly, while that for females (0.992, 95% CI 0.984 to 1.000) increased slightly, with no statistical significance.
To explore whether routine blood biochemical indices, together with age and sex, might be used to screen for suspected CHD in participants, we developed and externally validated screening models for males and females. Statistical analyses revealed that increasing age, TC, lnHCY and lnGGT, and decreasing UA, lnTG, apoA and apoB, increased the detection probability of CHD in males. In comparison, increasing age, TC, lnHCY, lnGGT and HDL-C versus LDL-C, and decreasing apoA, apoB, lnLp(a) and lnTBil, increased the diagnosis of CHD in females. ROC curves showed that the sensitivity and specificity of the models approached the optimum level in males and females. Both models were further externally validated in a separate population, with high performance. Several well-known risk factors for CHD, including age, TC, TG, apoA, apoB, HDL-C versus LDL-C, and HCY, were established as useful predictors of CHD in our screening model, while UA and TBil were also found to be predictors of CHD. In terms of UA, some studies have identified it as a risk factor for CHD,32 while other studies have shown that it is not a risk factor for CHD.39 These conflicting conclusions might be the result of a small sample size, unreasonable stratification, or inadequate statistical analysis. Although GGT, in clinical practice, is considered to be a biomarker for liver function, recent research has indicated that GGT is correlated with cardiovascular disease.40 Considering this fact, together with our data, the role of GGT in cardiovascular disease should be considered. Age, TC, HCY and GGT increased the risk of CHD, while apoA and apoB decreased the risk of CHD, although these factors had different coefficients among males and females in our screening models. UA and TG were found to be male-specific factors that decreased the risk of CHD; TBil and Lp(a) were found to be female-specific factors that decreased the risk of CHD, and HDL-C versus LDL-C was found to be a dominant factor increasing the risk of CHD. Further study is needed to determine why some risk factors are common to both sexes and why others are sex-specific.
Age, sex, family history and ethnicity are classic un-modifiable risk factors; hypertension, dyslipidaemia and diabetes, are classic modifiable risk factors; and smoking, alcohol abuse, obesity and high dietary intake of high-fat, high-carbohydrate and high-calorie diets, are behavioural risk factors.17 Considering the many biomarkers identified thus far,17 the risk factors for CHD form an intricate network. To date, it has been impossible to integrate all of the complicated relationships existing among these risk factors to produce a single screening strategy with high performance in clinical practice. This difficulty is attributed to the present technological threshold of computational biology, the intrinsic complexity of the pathophysiological process, and interactions between host and risk factors. Also, our understanding of the physiological processes underlying the metabolic balance in humans remains limited. Although UA, TC, TG, HDL-C, LDL-C, apoA, apoB, Lp(a), HCY, TBil, DBil, IBil and GGT were found to be associated with CHD,17–40 most studies have focused on blood fat and HCY; few studies have tried to integrate these biomarkers to make a comprehensive correlation between them and CHD. In our statistical analysis, we did not include other known risk factors for CHD, including family history, hypertension, diabetes, smoking, alcohol abuse, obesity and high dietary intake of high-fat, high-carbohydrate and high-calorie diets.17 The reasons for this omission are as follows: the use of a questionnaire to obtain certain information (eg, family history and high dietary intake of high-fat, high-carbohydrate and high-calorie diets) was not considered to be reliable enough; smoking and alcohol abuse are difficult to quantify; the causal relationship between some risk factors and CHD remains undefined (eg, whether obesity is an independent risk factor for CHD, a confounding factor, or a compensatory factor)17; and current two-dimensional and three-dimensional analysis strategies cannot integrate all of the information effectively (ie, having more factors is not necessarily better).
What is already known on this subject
Current cardiovascular risk-estimation systems include gender, age, smoking, cholesterol, systolic blood pressure etc, but neglect serum biochemical indices. Many serum biochemical indices, in addition to cholesterol, have been found to have changed with occurrence and developing of coronary heart disease; however, few studies have evaluated value of integrated biochemical indices on screening of coronary heart disease in suspected participants.
What this study adds
Statistical analyses revealed that increasing age, total cholesterol, homocysteine and γ-glutamyl transpeptidase, and decreasing uric acid, triglyceride, apoA and apoB, increased the diagnosis probability of coronary heart disease in males. Increasing age, TC total cholesterol, homocysteine, γ-glutamyl transpeptidase and high-density lipoprotein cholesterol versus low-density lipoprotein cholesterol, and decreasing apoA, apoB, lipoprotein (a) and total bilirubin, increased the diagnosis of coronary heart disease in females. These models were further validated in an external validation population, with high performance.
LF, CZ and SN contributed equally.
Contributors LF, JZ and LL designed the research; CZ and SN performed the statistical analysis; LF and CZ wrote the paper; YH, WX and XZ carried out the research.
Funding Yunnan Provincial Science and Technology Department of Basic Research on the Application of Self-Financing Projects 2013FZ257; the Regional Fund Project of the National Natural Science Foundation of China 81460326; the Scientific Research Fund of Yunnan Provincial Education Department 2011C083; Joint Special Funds from the Yunnan Province Science and Technology Department and Department of Applied Basic Research of Kunming Medical University 2013FZ283.
Competing interests None declared.
Patient consent Obtained.
Ethics approval The Review Board of the People's Hospital of Yuxi City.
Provenance and peer review Not commissioned; externally peer reviewed.
Data sharing statement All data are included in the paper.