Factor analysis for the clustering of cardiometabolic risk factors and sedentary behavior, a cross-sectional study

Abstract

Background

Few studies have reported on the clustering pattern of CVD risk factors, including sedentary behavior, systemic inflammation, and cadiometabolic components in the general population.

Objective

We aimed to explore the clustering pattern of CVD risk factors using exploratory factor analysis to investigate the underlying relationships between various CVD risk factors.

Methods

A total of 5606 subjects (3157 male, 51.5±11.7 y/o) were enrolled, and 14 cardiovascular risk factors were analyzed in an exploratory group (n = 3926) and a validation group (n = 1676), including sedentary behaviors.

Results

Five factor clusters were identified to explain 69.4% of the total variance, including adiposity (BMI, TG, HDL, UA, and HsCRP; 21.3%), lipids (total cholesterol and LDL-cholesterol; 14.0%), blood pressure (SBP and DBP; 13.3%), glucose (HbA1C, fasting glucose; 12.9%), and sedentary behavior (MET and sitting time; 8.0%). The inflammation biomarker HsCRP was clustered with only adiposity factors and not with other cardiometabolic risk factors, and the clustering pattern was verified in the validation group.

Conclusion

This study confirmed the clustering structure of cardiometabolic risk factors in the general population, including sedentary behavior. HsCRP was clustered with adiposity factors, while physical inactivity and sedentary behavior were clustered with each other.

About the Authors:

Tsung-Ying Tsai

Roles Conceptualization, Formal analysis, Investigation, Methodology, Writing – original draft

Affiliation: Division of Cardiology, Department of Medicine, Taipei Veterans General Hospital, Taipei, Taiwan

Pai-Feng Hsu

Roles Conceptualization, Data curation, Formal analysis, Methodology, Writing – review & editing

Affiliations Division of Cardiology, Department of Medicine, Taipei Veterans General Hospital, Taipei, Taiwan, Healthcare and Management Center, Taipei Veterans General Hospital, Taipei, Taiwan, Department of Public Health, National Yang-Ming University, Taipei, Taiwan, School of Medicine, National Yang-Ming University, Taipei, Taiwan

Chung-Chi Lin

Roles Data curation, Formal analysis, Methodology, Supervision

Affiliations Division of Cardiology, Department of Medicine, Taipei Veterans General Hospital, Taipei, Taiwan, School of Medicine, National Yang-Ming University, Taipei, Taiwan

Yuan-Jen Wang

Roles Methodology, Writing – review & editing

Affiliations Division of Cardiology, Department of Medicine, Taipei Veterans General Hospital, Taipei, Taiwan, School of Medicine, National Yang-Ming University, Taipei, Taiwan

Yaw-Zon Ding

Roles Supervision, Writing – review & editing

Affiliations Division of Cardiology, Department of Medicine, Taipei Veterans General Hospital, Taipei, Taiwan, School of Medicine, National Yang-Ming University, Taipei, Taiwan

Teh-Ling Liou

Roles Supervision, Writing – review & editing

Affiliations Division of Cardiology, Department of Medicine, Taipei Veterans General Hospital, Taipei, Taiwan, School of Medicine, National Yang-Ming University, Taipei, Taiwan

Ying-Wen Wang

Roles Supervision, Writing – review & editing

Affiliations Division of Cardiology, Department of Medicine, Taipei Veterans General Hospital, Taipei, Taiwan, School of Medicine, National Yang-Ming University, Taipei, Taiwan

Shao-Sung Huang

Roles Methodology, Supervision, Writing – review & editing

Affiliations Division of Cardiology, Department of Medicine, Taipei Veterans General Hospital, Taipei, Taiwan, Healthcare and Management Center, Taipei Veterans General Hospital, Taipei, Taiwan, School of Medicine, National Yang-Ming University, Taipei, Taiwan

Wan-Leong Chan

Roles Supervision, Writing – review & editing

Affiliations Division of Cardiology, Department of Medicine, Taipei Veterans General Hospital, Taipei, Taiwan, Healthcare and Management Center, Taipei Veterans General Hospital, Taipei, Taiwan, School of Medicine, National Yang-Ming University, Taipei, Taiwan

Shing-Jong Lin

Roles Supervision, Writing – review & editing

Affiliations Division of Cardiology, Department of Medicine, Taipei Veterans General Hospital, Taipei, Taiwan, Healthcare and Management Center, Taipei Veterans General Hospital, Taipei, Taiwan, School of Medicine, National Yang-Ming University, Taipei, Taiwan, Cardiovascular Research Center, National Yang-Ming University, Taipei, Taiwan, Institute of Clinical Medicine, National Yang-Ming University, Taipei, Taiwan

Jaw-Wen Chen

Roles Supervision, Writing – review & editing

Affiliations Division of Cardiology, Department of Medicine, Taipei Veterans General Hospital, Taipei, Taiwan, Healthcare and Management Center, Taipei Veterans General Hospital, Taipei, Taiwan, School of Medicine, National Yang-Ming University, Taipei, Taiwan

Hsin-Bang Leu

Roles Conceptualization, Investigation, Methodology, Supervision, Writing – review & editing

* E-mail: [email protected]

Affiliations Division of Cardiology, Department of Medicine, Taipei Veterans General Hospital, Taipei, Taiwan, Healthcare and Management Center, Taipei Veterans General Hospital, Taipei, Taiwan, School of Medicine, National Yang-Ming University, Taipei, Taiwan, Cardiovascular Research Center, National Yang-Ming University, Taipei, Taiwan, Institute of Clinical Medicine, National Yang-Ming University, Taipei, Taiwan

ORCID logo https://orcid.org/0000-0001-6832-3234

Abstract

Background

Few studies have reported on the clustering pattern of CVD risk factors, including sedentary behavior, systemic inflammation, and cadiometabolic components in the general population.

Objective

We aimed to explore the clustering pattern of CVD risk factors using exploratory factor analysis to investigate the underlying relationships between various CVD risk factors.

Methods

A total of 5606 subjects (3157 male, 51.5±11.7 y/o) were enrolled, and 14 cardiovascular risk factors were analyzed in an exploratory group (n = 3926) and a validation group (n = 1676), including sedentary behaviors.

Results

Five factor clusters were identified to explain 69.4% of the total variance, including adiposity (BMI, TG, HDL, UA, and HsCRP; 21.3%), lipids (total cholesterol and LDL-cholesterol; 14.0%), blood pressure (SBP and DBP; 13.3%), glucose (HbA1C, fasting glucose; 12.9%), and sedentary behavior (MET and sitting time; 8.0%). The inflammation biomarker HsCRP was clustered with only adiposity factors and not with other cardiometabolic risk factors, and the clustering pattern was verified in the validation group.

Conclusion

This study confirmed the clustering structure of cardiometabolic risk factors in the general population, including sedentary behavior. HsCRP was clustered with adiposity factors, while physical inactivity and sedentary behavior were clustered with each other.

Figures

Table 5

Table 1

Table 2

Table 3

Table 4

Fig 1

Fig 2

Table 5

Table 1

Table 2

Table 3

Citation: Tsai T-Y, Hsu P-F, Lin C-C, Wang Y-J, Ding Y-Z, Liou T-L, et al. (2020) Factor analysis for the clustering of cardiometabolic risk factors and sedentary behavior, a cross-sectional study. PLoS ONE 15(11): e0242365. https://doi.org/10.1371/journal.pone.0242365

Editor: Mauro Lombardo, San Raffaele Roma Open University, ITALY

Received: May 16, 2020; Accepted: October 30, 2020; Published: November 16, 2020

Copyright: © 2020 Tsai et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: Data cannot be shared publicly due to sensitive patient information. Data are available from the Taipei Veterans General Hospital Institutional Data Access / Ethics Committee ([email protected]) or via the corresponding author ([email protected]) for researchers who meet the criteria for access to confidential data.

Funding: This study is funded by the Research foundation of cardiovascular medicine. The primary recipient of research fund is Dr. Hsin-Bang Leu. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Background

Despite significant advances over the past decades, atherosclerotic cardiovascular diseases (CVD) remain the most significant cause of mortality worldwide [1]. CVD is a complex disease, and decades of research has recognized that CVD risk factors are clustered in specific patterns imply common underlying disease processes [2]. Metabolic abnormalities such as central obesity, insulin resistance, dyslipidemia, and high blood pressure are highly involved in the pathogenesis of CVD and are recognized as metabolic syndrome. Metabolic syndrome has variable presentation and is considered as a high-risk factor for CVD [3]. However, with the evolving understanding of the pathophysiology of atherosclerosis, risk factors are being added to the ever-expanding battery of CV risk factors.

Sedentary behavior (SB) and physical activity (PA) has been established as an independent risk factor for CVD on its own [4–7]. The current professional guidelines recommend avoiding SB and maintaining adequate PA, which are considered major goals in public health policy [8]. SB has also been linked to reduced triglyceride metabolism, insufficient antioxidant production, and glucose intolerance in several animal studies [9]. The underlying pathophysiologic link between SB and other CVD risk factors are deeply intertwined. There has been a great number of studies demonstrating the association between SB and tradition CVD risk factors. For example, in the landmark NHANES 2003–2006 study. total sedentary time was detrimentally associated with several biomarkers including waist circumference, HDL-cholesterol, C-reactive protein, triglycerides, insulin, and insulin resistance. However, while this study demonstrated the close relationship between SB and inflammation, the information of other important risk factors including total cholesterol, UA, HbA1C, were missing [10]. In a later meta-analysis, Edwardson et al. demonstrated that patients with longer sedentary time have greater odds of having metabolic syndrome [11]. However, it is impossible to delineate the association between SB and a particular component of metabolic syndrome. Hence, although the importance of SB is undeniable, there have only been a few studies have investigated the clustering relationship between SB and other risk factors, and there is insufficient information to explain the variances of metabolic abnormalities observed [10–14]. Furthermore, the definition of SB has not been well established in previous studies [15].

Exploratory factor analysis is a statistical method of data reduction that allows investigators to overcome the analytical difficulty posed by the vast number of risk variables by demonstrating the underlying relationships between different risk factors [16]. Factor analysis has been performed on many different populations but PA and SB have not been included in most major studies [17–19]. In addition, CVD has long been considered a systemic inflammatory disease. Inflammatory biomarkers such as high-sensitive C-reactive protein (HsCRP) are used to reflect disease severity and guide treatment strategies [20, 21]. Therefore, the current study investigates the clustering relationship of SB, cardiometabolic components, and inflammatory biomarkers in the general population of Taiwan.

Methods

Study population

We derived our data from the VGH-HEALTHCARE study, which is a prospective cohort study to evaluate the impact of PA and SB on long-term outcomes [22]. Adult subjects who received a comprehensive health examination at the Healthcare Center of Taipei Veterans General Hospital from February 2015 to July 2019 were invited to join this study. In brief, the healthcare center of Taipei Veterans general hospital provides elective, self-paid health examination services to all individuals who wish to receive health examination to identify undetected conditions for primary prevention purposes. We included patients without significant symptoms or illness and excluded those who refused to participate, whose exam revealed an acute illness, or had a chronic condition that require regular follow up such as active cancer, heart failure, coronary artery disease or stroke. We believe that data from our study population could provide vital information on primary prevention from the general population perspective. Demographic data, biochemical blood tests, and information on PA and SB were collected.

VGH-HEALTHCARE is an ongoing prospective study, so information on long-term outcomes is currently not available. The present study serves as a cross-section analysis to investigate the relationship between baseline cardiometabolic factors, SB, and inflammation biomarkers. The enrolled subjects were divided into two groups with similar baseline characteristics, with 70% in the exploratory group and 30% in the validation group. All participants provided a written consent. This study was conducted in concordance with the Declaration of Helsinki and was approved by the Internal Research Board of Taipei Veterans General Hospital. All information was obtained after receiving informed consent from the study participants.

Clinical assessments and biochemical parameters

The baseline information included age, sex, body height, body weight, body mass index (BMI), waist circumference, monthly income, education level, alcohol drinking behavior, and smoking status. The collected medical history included hypertension, type 2 diabetes, and hyperlipidemia. After an overnight fast, a TBA-c16000 automatic analyzer (Toshiba Medical Systems, Tochigi, Japan) was used to measure biochemical parameters, including fasting glucose, hemoglobin A1C (HbA1C), total cholesterol (TC), triglycerides (TG), low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol (HDL-C), aspartate aminotransferase (AST), alanine aminotransferase (ALT), total bilirubin, uric acid (UA), and high-sensitivity C-reactive protein (hsCRP) [23, 24].

Physical activity and sedentary behavior

PA levels and sedentary status were assessed using the validated Chinese version of the International Physical Activity Questionnaire–Short Form (IPAQ-SF) [25, 26]. All of the included patients filled out the questionnaires within one hour at noon with the assistance of specially trained nurses. The IPAQ-SF includes the number of days and the duration of vigorous, moderate, and walking activities during the previous week [27, 28]. The IPAQ-SF enables the calculation of metabolic equivalents (MET minutes per week), which are derived by assigning standardized MET values of 3.3, 4, and 8 for walking, moderate-intensity activity, and vigorous-intensity activity, respectively. These data were quantified, and an estimated metabolic equivalent of a task for each individual was classified as high, moderate, or low PA according to the IPAQ-SF score. The total daily sitting time was also collected for all participants.

Statistical analysis

The data are expressed as the mean ± standard deviation for normally distributed continuous variables and as number(percentage) for categorical variables. Demographic characteristics and biochemical variables were compared using student’s t-tests and the Mann–Whitney U test for the comparison of continuous variables, while chi-squared tests were used for categorical variables. Statistical significance was considered as P < 0.05. All statistical analyses were carried out using SPSS 20.0 software (IBM, Inc. Chicago, USA).

We performed an exploratory factor analysis to determine the clustering of cardiovascular risk factors, PA, and SB. The detailed statistical method has been described previously [29]. In brief, exploratory factor analysis is a statistical method designed to reveal the inter-correlations between the analyzed variables by reducing the collected variables into summary factors while retaining as much of the variance in the original variables as possible. There are three main steps of factor analysis: (1) extraction of factors (principal component analysis (PCA)); (2) rotation of factors to obtain a simple structure that can be easily interpreted; and (3) naming and interpreting each factor based on estimated values for factor loadings.

We used PCA to identify the principal components that reflect a group of variables that act together on a common pathophysiological process. We used an eigenvalue >1 as the extraction threshold, which was calculated as the sum of the squared factor loadings, which is a measurement of the amount of variation in the total sample accounted for by each factor. We then used orthogonal rotation (varimax rotation) to obtain factor loadings. We used an absolute loading value of >0.4 to interpret the factor pattern, which has been used by previous major factor analysis studies [30, 31].

In the primary analysis, factors were derived from 14 potential CVD risk factors: BMI, waist circumference, systolic blood pressure (SBP), diastolic blood pressure (DBP), fasting blood glucose, HbA1C, TC, HDL-C, LDL-C, TG, UA, HsCRP, MET, and daily sitting time. We first tested pair-wise correlations. We defined a KMO value greater than 0.6 and a significant Bartlett’s test of sphericity (P<0.001) as an indication for sampling adequacy and a lack of an identity matrix. Because it is difficult to find another population with the same low risk who had also completed the questionnaires as well as a biochemical study including CRP, we consider that it was reasonable to divide the whole population into two groups, one for training and one for validation to avoid overfitting the developed model. We randomly assigned 70% of the total population to the training group while the remaining 30% were assigned to the validation group. In statistics for data mining, it is a common method that training data is often a subset of the total data set, and the test set is a subset of the test-trained model. Similar analyses have been reported previously. For example, in Hravnak’s machine learning algorithms data were divided into Block 1 for the ML training/cross-validation set and Block 2 for the test set [32]. Furthermore, Goodman et al. used factor analysis for the cardiovascular clustering risk, and they selected 20% of the cases as an exploratory sample and the remaining 80% cases were used as a validation sample [30]. Because our current study aimed to investigate the cluster risk including inflammation makers and SB information, it is reasonable to test and validate sample using the same number of variables.

Results

A total of 5606 subjects (3157 male, age 51.5±11.7) were enrolled in this study. The enrollees were divided into a training group (N = 3926) and a validation group (N = 1676). The participants’ baseline characteristics, biochemical data, PA, and SB are shown in Table 1. The enrolled population consists of healthy middle age Taiwanese people. There was no significant difference between the demographic and biochemical data of the two groups.

[Figure omitted. See PDF.]

Table 1. Baseline characteristics, comorbidities, and physical activity of the two groups.

https://doi.org/10.1371/journal.pone.0242365.t001

The correlation among the 14 variables among all subjects is shown in Table 2. Sitting time demonstrated a negative correlation with MET (Pearson correlation coefficient (r) = -0.165, P < 0.001). MET value was significantly negatively correlated with waist circumference, blood pressure, cholesterol, triglyceride, LDL and sitting time, but positively correlate with HDL. We have further categories subjects into low, moderate, and high PA group and we found that subjects with low activity have unfavorable lipid profiles and higher baseline inflammatory makers, supporting the connection between higher cardiovascular risk to low PA (Table 3). However, although sitting time is negatively correlated with MET values, the correlation between lipid profiles only existed in HDL. It is not surprising because sitting time is only one factor among the definition of SB and total MET activity estimation was generated from more activity’s information. To explore the effect of gender and age of factor clustering pattern, we re-analyzed the correlation, taking into consideration gender and age which showed similar factor clustering in all subgroups (S1–S4 Tables).

[Figure omitted. See PDF.]

Table 2. Correlation of cardiovascular risk factors.

https://doi.org/10.1371/journal.pone.0242365.t002

[Figure omitted. See PDF.]

Table 3. Baseline characteristics, comorbidities, and physical activity of different physical activity groups.

https://doi.org/10.1371/journal.pone.0242365.t003

Tables 3 and 4 show the factor analysis results of the training group and the validation group. In the training group (Table 3), PCA identified five factors with an eigenvalue >1. The combined factors explained 69.43% of the variance among the original 14 factors. BMI, WC, HsCRP, TG, UA, and HDL-C were grouped together in the first common factor, the adiposity factor, which is similar to the metabolic syndrome criteria put forth by the NCEP/ATP III guidelines [33]. The inflammation biomarker hsCRP, which represents the underlying inflammation status of a subject, was shown within the adiposity group. This factor explained 21.30% of the total variance. The second common factor was the lipid factor, which contained LDL-C and total cholesterol and explained approximately 13.97% of the total variance.

[Figure omitted. See PDF.]

Table 4. Factor analysis of the training group.

https://doi.org/10.1371/journal.pone.0242365.t004

The third common factor, the blood pressure factor, consisted of SBP and DBP and accounted for 13.03% of the variance. The fourth common factor, the glucose factor, included fasting blood glucose and HbA1C and accounted for 12.91% of the variance. The final common factor, the SB factor, consisted of both the MET and the daily sitting time, which explained 7.96% of the variance.

For the results of the validation samples (Table 5), the PCA also identified five factors. The combined factors explained 69.48% of the variance in the original 14 factors. The adiposity factor, lipid factor, blood pressure factor, glucose factor, and activity factor explained 20.26%, 13.96%, 13.81%, 13.03, and 8.43% of the variance, respectively. Figs 1 and 2 show the component plots with factor diagrams from the PCA with varimax rotation. Both groups demonstrated a consistent clustering of risk factors.

[Figure omitted. See PDF.]

Fig 1. Component plots of the training group with factor diagrams from principle component analysis with varimax rotation.

BMI = body mass index; GLU = serum glucose; HDL = high density lipoprotein; LDL = low density lipoprotein; SBP = systolic blood pressure; DBP = diastolic blood pressure; MET = metabolic equivalent; HbA1C = hemoglobin A1C.

https://doi.org/10.1371/journal.pone.0242365.g001

[Figure omitted. See PDF.]

Fig 2. Component plots of the validation group with factor diagrams from principle component analysis with varimax rotation.

BMI = body mass index; GLU = serum glucose; HDL = high density lipoprotein; LDL = low density lipoprotein; SBP = systolic blood pressure; DBP = diastolic blood pressure; MET = metabolic equivalent; HbA1C = hemoglobin A1C.

https://doi.org/10.1371/journal.pone.0242365.g002

[Figure omitted. See PDF.]

Table 5. Factor analysis of the validation group.

https://doi.org/10.1371/journal.pone.0242365.t005

To explore the composition of factor clustering in different gender and age, we performed age and sex-stratified factor analysis. The results showed a similar clustering pattern of risk factors in both gender and age groups (S1–S4 Tables).

Discussion

This single-center cross-sectional analysis examined 5606 healthy Asian adults and demonstrated that complex clustering cardiometabolic factors can be divided into five factor clusters: the adiposity factor (waist circumference, BMI, TG, HDL, and UA), the blood pressure factor (SBP and DBP), the lipid factor (TC and LDL), the glucose factor (fasting glucose and HbA1C), and the PA factor. These factors explained 21.97%, 13.97%, 13.30%, 12.91%, and 7.96% of the total variance, respectively. Systemic inflammation was linked to the adiposity factor, while SB and PA were clustered together and formed an independent CVD risk factor. Multiple cardiometabolic factors were involved in the development and progression of atherosclerotic cardiovascular disease.

This study is the first to investigate the relationship of cardiometabolic factors, systemic inflammation, and sedentary information simultaneously in the general population of Taiwan. We also demonstrated that the PA factors do not cluster with traditional CV risk factors. Our results suggest that physical inactivity may exert its effect on cardiovascular disease in an independent and unique way. This result may prompt future researchers to explore the possible pathophysiologic mechanism behind the independent effect of the PA level. Our findings could also provide important evidence that adiposity is linked to baseline inflammation, which explained nearly 20% of the variance in subjects without CVD.

Systemic inflammation and cardiovascular risk factors

Increased baseline inflammation is believed to play a crucial role in all stages of the artherothrombotic disease process, and treatment strategies to reduce inflammation have been demonstrated to lower the residual risk of CVD in addition to lipid-lowering medications [20, 21, 34]. We demonstrated that HsCRP was a part of the adiposity factor but was independent of the lipid factor (TC and LDL). The correlation between baseline CRP and LDL has been a controversial topic in previous observations [35]. The JUPITER trial demonstrated that risk reduction with statin therapy is related to the level of CRP, but no such relationship was observed for LDL-C. This suggests that statins reduced hsCRP independently of LDL and that LDL was not related to hsCRP [36, 37].

The current study extended the understanding that inflammation is linked to the adiposity factor and contributes to cardiometabolic clustering variance in the general population of Taiwan. Nicklas et al. showed that the change in hsCRP after weight loss is significantly correlated with changes in total body fat, abdominal adiposity, visceral adiposity, and lipid profiles, especially TG and HDL, but not LDL. This supports the close relationship between baseline inflammation and adiposity [38]. It is also worth mentioning that uric acid was also linked with adiposity in our study. This result is consistent with a previous factor analysis study of 2945 adults from the FIBER study [39]. Uric acid has been recognized as a maker for inflammation and oxidative stress, which underlie the disease processes of gout and CVD [40]. The inclusion of uric acid further supports our finding that insulin resistance and inflammation are features of adiposity.

Clustering of sedentary behavior and physical inactivity

SB and PA can explain some of the variance in our study, but the proportion of explained variance was smaller than that of other factors. We also demonstrated that daily sitting time has a negative correlation with weekly PA, suggesting a close relationship between them. The correlation between sitting time and PA has not been consistent in previous studies, with some demonstrating that PA and SB are interdependent, while others show an independent relationship [7, 41–43]. This discrepancy can be explained by the fact that a person can be both sedentary and physically active (the Active Couch Potato phenomenon, describes someone who meets the recommendations for physical activity but still sits around for long periods of the day) [44]. Another explanation is the inconsistent definition of SB across previous studies. Many studies used the lowest MET as a definition of SB, which does not require the subject to be actually sitting [15]. In this study, we used the IPAQ questionnaire, which is a widely validated questionnaire that is recommended by the WHO, to evaluate the association between sedentary behavior and weekly physical activity [45, 46]. Our result is compatible with a recent large study with robust methodology by Stamatakis et al., who showed that the daily sitting time and PA of 149,077 middle age adults were well correlated [42].

However, MET and daily sitting time were not significantly correlated with other cardiometabolic risk factors in our study. To evaluate whether PA and SB are truly independent of other risk factors, we divided our subjects into low, mediate, and high-intensity PA groups according to the definition of the IPAQ questionnaire. Table 3 shows the baseline characteristics of all subjects according to PA. Subjects with low PA have longer sitting times, larger waist circumferences, unfavorable lipid profiles (higher TG, higher LDL-C, and lower HDL-C) and higher baseline CRP. This indicates that low PA was still associated with risk factors of cardiovascular disease. Indeed, as shown in Table 3, patients with low PA have multiple factors at the same time, making it difficult to evaluate the importance and correlation of individual risk factors. Factor analysis of this complex data has provided more information to dissect the complex clustering of risk factors.

Table 3 have also showed that nearly 43% of our participants belonged to the low PA categories. Although there might be a selection bias that our population may have higher socioeconomic status, but lower daily activity, the high prevalence of low PA is not unique to our population. Both global and Taiwanese survey have showed equally alarming rate of insufficiently physically activity [47, 48]. Moreover, the WHO report have shown that over the past 15 years, the levels of insufficient activity did not improve [49].

Study limitations

There were some limitations to the present study. First, the information on sitting time was self-reported and is thus susceptible to reporting error. However, IPAQ was developed to measure health-related PA in the general population, has been tested extensively, and is now used in many international studies [45]. The sitting time was defined as the time subjects spend sitting while at work, at home, and during leisure time. This may include time spent sitting at a desk, socializing with friends, reading, sitting, or lying down watching television. The screen time spent on individual activities such as TV watching or computer/mobile usage was not recorded individually. Second, all of the participants were recruited from those who received a comprehensive annual or biennial examination at a healthcare center, and all subjects were asymptomatic and had few risk factors. A detailed history of medication was not available, and the possibility of selection bias cannot be excluded. Larger studies are also still needed to investigate whether our findings could be applied to other populations. Third, the study was cross-sectional in nature and provided no outcome data. Because the VGH-HEALTHCARE study is still ongoing, the long-term outcome data will be available when the VGH- HEALTHCARE study completes. Fourth, our study did not include emerging risk factors for CVD such as hematological factors and liver function. However Although these factors are considered to be new cardiovascular markers in some studies, there is little evidence supporting their role in the general population compared with the factors analyzed in the current study. Fifth, we did not apply machine learning method in the analysis of our data. We plan to use machine learning method for the analysis of outcome data after the VGH-HEALTHCARE study is completed. However, we did not consider using machine learning in our current study, for two reasons. First, using machine learning for factor analysis has not been widely accepted for lack of long-term outcome information. Second. To make our results more convincing to the general audience, we chose to investigate the percentage of variance of clustering cardiometabolic risk factors and the association between sedentary behavior and other associated risk factors with a well validated statistical method.

Conclusion

The complex relationship of cardiometabolic factors, inflammation, and sedentary information among the general population of Taiwan can be divided into five factor clusters: the adiposity factor (waist circumference, BMI, TG, HDL, and UA), the blood pressure factor (SBP and DBP), the lipid factor (TC and LDL), the glucose factor (fasting glucose and HbA1C), and the PA factor, which explained 21.97%, 13.97%, 13.30%, 12.91%, and 7.96% of the total variance, respectively. Our results suggest that systemic inflammation shares the same underlying disease process with metabolic syndrome, while the independent role of PA warrants exploration with future studies.

Supporting information

[Figure omitted. See PDF.]

S1 Table. Factor analysis in patients <65 years of age.

https://doi.org/10.1371/journal.pone.0242365.s001

(DOCX)

S2 Table. Factor analysis in patients >65 years of age.

https://doi.org/10.1371/journal.pone.0242365.s002

(DOCX)

S3 Table. Factor analysis in male patients.

https://doi.org/10.1371/journal.pone.0242365.s003

(DOCX)

S4 Table. Factor analysis in female patients.

https://doi.org/10.1371/journal.pone.0242365.s004

(DOCX)

Citation: Tsai T-Y, Hsu P-F, Lin C-C, Wang Y-J, Ding Y-Z, Liou T-L, et al. (2020) Factor analysis for the clustering of cardiometabolic risk factors and sedentary behavior, a cross-sectional study. PLoS ONE 15(11): e0242365. https://doi.org/10.1371/journal.pone.0242365

References

1. Roth GA, Johnson C, Abajobir A, Abd-Allah F, Abera SF, Abyu G, et al. Global, Regional, and National Burden of Cardiovascular Diseases for 10 Causes, 1990 to 2015. J Am Coll Cardiol. 2017;70: 1–25. pmid:28527533

2. Reaven GM. Banting lecture 1988. Role of insulin resistance in human disease. Diabetes. 1988;37: 1595–607. pmid:3056758

3. Kassi E, Pervanidou P, Kaltsas G, Chrousos G. Metabolic syndrome: definitions and controversies. BMC Med. 2011;9: 48. pmid:21542944

4. Sattelmair J, Pertman J, Ding EL, Kohl HW, Haskell W, Lee I-M. Dose response between physical activity and risk of coronary heart disease: a meta-analysis. Circulation. 2011;124: 789–95. pmid:21810663

5. Young DR, Hivert M-F, Alhassan S, Camhi SM, Ferguson JF, Katzmarzyk PT, et al. Sedentary Behavior and Cardiovascular Morbidity and Mortality: A Science Advisory From the American Heart Association. Circulation. 2016;134: e262–79. pmid:27528691

6. Pandey A, Salahuddin U, Garg S, Ayers C, Kulinski J, Anand V, et al. Continuous dose-response association between sedentary time and risk for cardiovascular disease a meta-analysis. JAMA Cardiol. 2016;1: 575–583. pmid:27434872

7. Ekelund U, Bauman A, Lee IM. Effects of early physical exercise on later health–Authors’ reply. The Lancet. Lancet Publishing Group; 2017. p. 801. pmid:28248173

8. Piepoli MF, Hoes AW, Agewall S, Albus C, Brotons C, Catapano AL, et al. 2016 European Guidelines on cardiovascular disease prevention in clinical practice: The Sixth Joint Task Force of the European Society of Cardiology and Other Societies on Cardiovascular Disease Prevention in Clinical Practice (constituted by representatives of 10 societies and by invited experts)Developed with the special contribution of the European Association for Cardiovascular Prevention & Rehabilitation (EACPR). Eur Heart J. 2016;37: 2315–2381. pmid:27222591

9. Matthews CE. Minimizing Risk Associated With Sedentary Behavior: Should We Focus on Physical Activity, Sitting, or Both? Journal of the American College of Cardiology. Elsevier USA; 2019. pp. 2073–2075. pmid:31023431

10. Healy GN, Matthews CE, Dunstan DW, Winkler EAH, Owen N. Sedentary time and cardio-metabolic biomarkers in US adults: NHANES 200306. Eur Heart J. 2011;32: 590–597. pmid:21224291

11. Edwardson CL, Gorely T, Davies MJ, Gray LJ, Khunti K, Wilmot EG, et al. Association of Sedentary Behaviour with Metabolic Syndrome: A Meta-Analysis. O’Connor KA, editor. PLoS One. 2012;7: e34916. pmid:22514690

12. Sungwacha JN, Tyler J, Longo-Mbenza B, Lasi On’Kin JBK, Gombet T, Erasmus RT. Assessing clustering of metabolic syndrome components available at primary care for Bantu Africans using factor analysis in the general population. BMC Res Notes. 2013;6: 228. pmid:23758878

13. Esteghamati A, Zandieh A, Khalilzadeh O, Morteza A, Meysamie A, Nakhjavani M, et al. Clustering of leptin and physical activity with components of metabolic syndrome in Iranian population: An exploratory factor analysis. Endocrine. 2010;38: 206–213. pmid:21046482

14. Kelishadi R, Ardalan G, Adeli K, Motaghian M, Majdzadeh R, Mahmood-Arabi MS, et al. Factor Analysis of Cardiovascular Risk Clustering in Pediatric Metabolic Syndrome: CASPIAN Study. Ann Nutr Metab. 2007;51: 208–215. pmid:17587791

15. Tremblay MS, Aubert S, Barnes JD, Saunders TJ, Carson V, Latimer-Cheung AE, et al. Sedentary Behavior Research Network (SBRN)–Terminology Consensus Project process and outcome. Int J Behav Nutr Phys Act. 2017;14: 75. pmid:28599680

16. Meigs JB. Invited commentary: insulin resistance syndrome? Syndrome X? Multiple metabolic syndrome? A syndrome at all? Factor analysis reveals patterns in the fabric of correlated metabolic risk factors. Am J Epidemiol. 2000;152: 908–11; discussion 912. pmid:11092432

17. Choi KM, Lee J, Kim KB, Kim DR, Kim SK, Shin DH, et al. Factor analysis of the metabolic syndrome among elderly Koreans—the South-west Seoul Study. Diabet Med. 2003;20: 99–104. pmid:12581260

18. H A J, K A J, F A, D R, W LE, S P, et al. Factor Analysis of Metabolic Syndrome Using Directly Measured Insulin Sensitivity: The Insulin Resistance Atherosclerosis Study. Diabetes. 2002;51. pmid:12145182

19. Tsai C-H, Li T-C, Lin C-C, Tsay H-S. Factor analysis of modifiable cardiovascular risk factors and prevalence of metabolic syndrome in adult Taiwanese. Endocrine. 2011;40: 256–64. pmid:21499818

20. Emerging Risk Factors Collaboration K, Kaptoge S, DiAngelantonio E, Lowe G, Pepys MB, Thompson SG, et al. C-reactive protein concentration and risk of coronary heart disease, stroke, and mortality: an individual participant meta-analysis. Lancet (London, England). 2010;375: 132–40. pmid:20031199

21. Ridker PM, Everett BM, Thuren T, MacFadyen JG, Chang WH, Ballantyne C, et al. Antiinflammatory Therapy with Canakinumab for Atherosclerotic Disease. N Engl J Med. 2017;377: 1119–1131. pmid:28845751

22. Yang H-C, Liang Y, Hsu H-C, Shu J-H, Chou R-H, Hsu P-F, et al. InVestiGation of the Association of Physical Activity and Sedentary Behavior with tHe Occurrence of Future Cardiovascular Disease and Long Term Outcome in General Population Using the HEALTHCARE Database (VGH-HEALTHCARE). Acta Cardiol Sin. 2019;35: 534–541. pmid:31571803

23. Chen S-C, Lin C-P, Hsu H-C, Shu J-H, Liang Y, Hsu P-F, et al. Serum bilirubin improves the risk predictions of cardiovascular and total death in diabetic patients. Clin Chim Acta. 2019;488: 1–6. pmid:30359587

24. Chiang C-H, Huang C-C, Chan W-L, Chen J-W, Leu H-B. The severity of non-alcoholic fatty liver disease correlates with high sensitivity C-reactive protein value and is independently associated with increased cardiovascular risk in healthy population. Clin Biochem. 2010;43: 1399–404. pmid:20846522

25. Liou YM, Jwo CJC, Yao KG, Chiang L-C, Huang L-H. Selection of appropriate Chinese terms to represent intensity and types of physical activity terms for use in the Taiwan version of IPAQ. J Nurs Res. 2008;16: 252–63. pmid:19061172

26. Hwang A-C, Zhan Y-R, Lee W-J, Peng L-N, Chen L-Y, Lin M-H, et al. Higher Daily Physical Activities Continue to Preserve Muscle Strength After Mid-Life, But Not Muscle Mass After Age of 75. Medicine (Baltimore). 2016;95: e3809. pmid:27258519

27. Guthold R, Ono T, Strong KL, Chatterji S, Morabia A. Worldwide variability in physical inactivity a 51-country survey. Am J Prev Med. 2008;34: 486–94. pmid:18471584

28. Brugnara L, Murillo S, Novials A, Rojo-Martínez G, Soriguer F, Goday A, et al. Low Physical Activity and Its Association with Diabetes and Other Cardiovascular Risk Factors: A Nationwide, Population-Based Study. Vina J, editor. PLoS One. 2016;11: e0160959. pmid:27532610

29. Williams B, Onsman A, Brown T, Brown T. Exploratory factor analysis: A five-step guide for novices. Australas J Paramed. 2010;8.

30. Goodman E, Dolan LM, Morrison JA, Daniels SR. Factor analysis of clustered cardiovascular risks in adolescence: obesity is the predominant correlate of risk among youth. Circulation. 2005;111: 1970–7. pmid:15837951

31. Stoner L, Weatherall M, Skidmore P, Castro N, Lark S, Faulkner J, et al. Cardiometabolic Risk Variables in Preadolescent Children: A Factor Analysis. J Am Heart Assoc. 2017;6. pmid:29021277

32. Hravnak M, Chen L, Dubrawski A, Bose E, Clermont G, Pinsky MR. Real alerts and artifact classification in archived multi-signal vital sign monitoring data: implications for mining big data. J Clin Monit Comput. 2016;30: 875–888. pmid:26438655

33. Expert Panel on Detection, Evaluation, and Treatment of High Blood Cholesterol in Adults. Executive Summary of The Third Report of The National Cholesterol Education Program (NCEP) Expert Panel on Detection, Evaluation, And Treatment of High Blood Cholesterol In Adults (Adult Treatment Panel III). JAMA. 2001;285: 2486–97. pmid:11368702

34. Libby P, Ridker PM, Maseri A. Inflammation and Atherosclerosis. Circulation. 2002;105: 1135–1143. pmid:11877368

35. Christian AH, Mochari H, Mosca LJ. Waist Circumference, Body Mass Index, and Their Association With Cardiometabolic and Global Risk. J Cardiometab Syndr. 2009;4: 12–19. pmid:19245511

36. Ridker PM. Moving Beyond JUPITER: Will Inhibiting Inflammation Reduce Vascular Event Rates? Curr Atheroscler Rep. 2013;15: 295. pmid:23225175

37. Ridker PM, Danielson E, Fonseca FAH, Genest J, Gotto AM, Kastelein JJP, et al. Rosuvastatin to Prevent Vascular Events in Men and Women with Elevated C-Reactive Protein. N Engl J Med. 2008;359: 2195–2207. pmid:18997196

38. Nicklas JM, Sacks FM, Smith SR, LeBoff MS, Rood JC, Bray GA, et al. Effect of dietary composition of weight loss diets on high-sensitivity c-reactive protein: The Randomized POUNDS LOST trial. Obesity. 2013;21: 681–689. pmid:23712970

39. Mannucci E, Monami M, Rotella CM. How many components for the metabolic syndrome? Results of exploratory factor analysis in the FIBAR study. Nutr Metab Cardiovasc Dis. 2007;17: 719–726. pmid:17387006

40. Maruhashi T, Hisatome I, Kihara Y, Higashi Y. Hyperuricemia and endothelial function: From molecular background to clinical perspectives. Atherosclerosis. 2018;278: 226–231. pmid:30326405

41. Schuna JM, Johnson WD, Tudor-Locke C. Adult self-reported and objectively monitored physical activity and sedentary behavior: NHANES 2005–2006. Int J Behav Nutr Phys Act. 2013;10: 126. pmid:24215625

42. Stamatakis E, Gale J, Bauman A, Ekelund U, Hamer M, Ding D. Sitting Time, Physical Activity, and Risk of Mortality in Adults. J Am Coll Cardiol. 2019;73: 2062–2072. pmid:31023430

43. Silfee V, Lemon S, Lora V, Rosal M. Sedentary behavior and cardiovascular disease risk factors among latino adults. J Health Care Poor Underserved. 2017;28: 798–811. pmid:28529225

44. Owen N, Healy GN, Matthews CE, Dunstan DW. Too much sitting: The population health science of sedentary behavior. Exerc Sport Sci Rev. 2010;38: 105–113. pmid:20577058

45. Craig CL, Marshall AL, Sjöström M, Bauman AE, Booth ML, Ainsworth BE, et al. International physical activity questionnaire: 12-Country reliability and validity. Med Sci Sports Exerc. 2003;35: 1381–1395. pmid:12900694

46. Boon RM, Hamlin MJ, Steel GD, Ross JJ. Validation of the New Zealand physical activity questionnaire (NZPAQ-LF) and the international physical activity questionnaire (IPAQ-LF) with accelerometry. Br J Sports Med. 2010;44: 741–746. pmid:18981036

47. CDC maps show variance in adult physical activity levels by state, territory | AHA News. [cited 11 Aug 2020]. https://www.healthycommunities.org/news/headline/2020-01-22-cdc-maps-show-variance-adult-physical-activity-levels-state-territory

48. Wu X, Tsai SP, Tsao CK, Chiu ML, Tsai MK, Lu PJ, et al. Cohort Profile: The Taiwan MJ Cohort: Half a million Chinese with repeated health surveillance data. Int J Epidemiol. 2017;46: 1744–1744g. pmid:28204597

49. WHO | Prevalence of insufficient physical activity. WHO. 2018.