Data source and study population
A cross-sectional study of adult residents was conducted in Catalonia, a Mediterranean region of southern Europe with 7,434,632 inhabitants (16 % of the Spanish population, 2010 census). In Catalonia, 358 primary health care teams (PHCTs) comprised of general practitioners (GPs), nurses, social workers and support personnel are assigned by geographical area and responsible for the health care of the population in their areas. The Catalan Health Institute manages 274 PHCTs (76.5 %), serving a population of 4,859,725 adults; the remaining PHCTs are managed by other providers. Primary care professionals systematically use electronic health records (EHR) to record diagnoses, prescriptions and other clinical information, patient management, and administrative activities. The Catalan Health Institute Information System for the Development of Research in Primary Care (SIDIAP) [11] compiles coded clinical information from the EHR system based on data from its 274 PHCTs. A subset of SIDIAP records meeting the highest quality criteria for clinical data (SIDIAP-Q) includes 40 % of the SIDIAP population (1,936,443 patients), attended by the 1,319 GPs assigned to 251 PHCTs whose data recording scored highest in a validated comparison process [12]. SIDIAP has been shown to be highly representative of the general Catalan population in terms of geography, age and sex distribution, according to the official 2010 census [12].
The study sample, was selected from the SIDIAP-Q database, included 1,749,710 patients aged 19 years or older, assigned to 251 PHCTs during the period of study (1 January- 31 December 2010); 186,733 individuals were excluded because they were younger than 19 years.
Coding of diseases
International Classification of Diseases (ICD-10) codes were mapped to the International Classification of Primary Care (ICPC-2e-v.4.2, available at: http://www.kith.no/templates/kith_WebPage____1111.aspx). R codes (symptoms, signs and abnormal clinical and laboratory findings, not elsewhere classified) and Z codes (factors influencing health status and contact with health services) were excluded, resulting in 686 included codes. Each diagnosis was then classified using O’Halloran criteria for chronic disease [13]: (i) have a duration that has lasted, or is expected to last, at least 6 months; (ii) have a pattern of recurrence or deterioration; (iii) have a poor prognosis and (iv) produce consequences, or sequel, that have an impact on the individual’s quality of life [13, 14]. All results were described by ICPC-2 codes and 146 chronic diseases were included in the analysis.
Cardiovascular morbidity was studied in 24 selected cardiovascular chronic diseases from chapter K of ICPC- 2, based on active diagnoses recorded in EHR as of December 31, 2010: Rheumatic fever/heart disease (K71); Neoplasm, cardiovascular (K72); Congenital anomaly, cardiovascular (K73); Ischaemic heart disease with angina (K74); Acute myocardial infarction (K75); Ischaemic heart disease without angina (K76); Heart failure (K77); Atrial fibrillation/flutter (K78); Paroxysmal tachycardia (K79); Cardiac arrhythmia NOS (K80); Heart/arterial murmur NOS (K81); Pulmonary heart disease (K82); Heart valve disease NOS (K83); Heart disease, other (K84); Hypertension, uncomplicated (K86); Hypertension, complicated (K87); Postural hypotension (K88); Transient cerebral ischaemia (K89); Stroke/cerebrovascular accident (K90); Cerebrovascular disease (K91); Atherosclerosis/peripheral vascular disease (K92); Pulmonary embolism (K93); Phlebitis/thrombophlebitis (K94) and Varicose veins of leg (K95).
Outcomes and variables
The main outcome was cardiovascular morbidity burden defined as the coexistence of one or more chronic cardiovascular diseases in patients with MM. MM was defined as the coexistence of 2 or more chronic diseases.
Two study groups were defined, with and without MM. The MM group was further divided into subgroups, constituting a cardiovascular morbidity group (MM-CMG) –i.e., MM patients with one or more chronic cardiovascular disease– and non-cardiovascular morbidity group with other chronic diseases (MM-NCMG).
Secondary outcomes included CVR and CVRF profile. CVR was assessed by the REGICOR (Registre Gironí del Cor) score, which evaluates the 10-year risk of a coronary event (angina, myocardial infarct with/without symptoms, fatal or non-fatal), with four categories of severity: low, <5 %; moderate, 5–9.9 %; high, 10–14.9 %; and very high, ≥ 15 % [15]. This score is only applicable to individuals aged 35 to 74 years.
Modifiable major CVRFs registered in the EHR and the sum of these major factors were analysed: smoking, hypertension, hypercholesterolaemia and diabetes. Other CVRFs evaluated were: hypertriglyceridaemia, obesity, and alcoholism (in average units of weekly consumption, classified as: low risk consumption [17–28 units in men and 11–17 units in women]; risky consumption [>28 and > 17, respectively]) [16, 17]. Five additional variables were considered in the analysis: sex (female/male), age (years), age groups (19–24, 25–44, 45–64, 65–79, and 80+ years), number of chronic diseases and setting (urban/rural). Physical examinations yielded values for body mass index (BMI, kg/cm2) and blood pressure (mm Hg) and included laboratory tests: glycated haemoglobin (%), creatinine (mg/dl), uric acid (mg/dl), total cholesterol (mg/dl) and triglycerides (mg/dl), along with glomerular filtration rate < 60 ml/min/1.73 m2 to determine a decrease in renal function [18].
Statistical analysis
Descriptive statistics were used to summarize overall information. Categorical variables were expressed as frequencies (percentage) and continuous as mean (Standard deviation, SD) or median (interquartile range, IQR). The MM crude prevalence and 95 % confidence intervals (CI) were calculated. The differences between MM and non-MM groups were tested using Student t, Mann–Whitney or Chi-square for unadjusted comparison, as appropriate. The crude prevalence (95 % CI) of MM-CMG was calculated. Prevalence estimates of each CV chronic condition and 95 % CI were obtained. The 95 % CI for the prevalences was calculated using the continuity-corrected Wilson score interval.
For comparison of MM-CMG, MM-NCMG and non-MM groups, ANOVA, Kruskal-Wallis or Chi square tests were used as appropriate. To determine the most prevalent MM patterns in MM-CMG, all possible combinations of each CV disease, one and two chronic conditions and their frequencies were calculated for each sex and age group. CVR score (REGICOR) distribution was compared within the three groups studied..We assumed that missing data were missing completely at random (MCAR), and so we performed a complete case analysis to handle missing data. We had sufficient power for our analysis, even though we lost part of our data set.
Spearman correlation between number of chronic diseases and chronic cardiovascular diseases and age was assessed.
All statistical tests were two-sided at the 5 % significance level. The analyses were performed using SPSS for Windows, version 18 and R version 3.2.3.