Introduction

Coronavirus disease 2019 (COVID-19), a respiratory disease caused by the severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2), was rapidly expanded worldwide within a short period. The new SARS-CoV-2 coronavirus infection is clinically characterized by mild upper respiratory tract infection [1]; however, some cases may progress to aggressive viral pneumonia, acute respiratory failure and death [2].

SARS-CoV-2 coronavirus infection has become challenging for public health systems because of the high infection rates. It is crucial to understand the epidemiology and pathobiology of COVID-19 to reduce its impact. The rapid evolving pandemic requires sharing all information and instant experiences among healthcare workers in order to assess prognosis, risk stratification, allocation of healthcare resources and designing of clinical trials. Control strategies and isolation are necessary for prevention while innovative and costly therapies, including hospital beds and resources, are required for management. Implementing these measures requires considerable efforts by healthcare systems; thus, quantifying the influence of risk factors on clinical outcomes at the population level is critical. As for the literature revised, to date few analyses of this sort have been performed.

Since COVID-19 is a new infectious disease, its relationship with chronic comorbidities is still unclear. According to four systematic reviews and meta-analyses, chronic underlying diseases are highly prevalent among COVID-19 patients and are associated with severe cases [3,4,5,6]. Elderly people and those with chronic comorbidities such as respiratory, cardiovascular, metabolic, renal, and cancer have been associated with adverse clinical outcomes [6,7,8,9,10,11]. This study was aimed to determine the effects of chronic comorbidities on clinical outcomes of patients with and without COVID-19.

Methods

Design, setting and patients

The study was conducted in accordance with The REporting of studies Conducted using Observational Routinely-collected health Data (RECORD) Statement. This is a retrospective analysis of the national administrative healthcare open data from the General Directorate of Epidemiology of Mexico [12]. National administrative healthcare open data routinely collects information regarding origin, demography, epidemiology, chronic comorbidities and outcomes of patients with suspicion of viral respiratory disease, defined as cough, fever or headache last week more at least one of the following: dyspnoea, arthralgia, myalgia, sore throat, rhinorrhoea, conjunctivitis, or chest pain [13]. A complete list of codes used to classify chronic comorbidities, outcomes, and confounders is provided in Online Resource 1. The sentinel healthcare units of the System of Epidemiological Surveillance of Viral Respiratory Diseases, located at the first, second or third level of healthcare, transmitted the data to the institutional/ regional health authorities and then recorded in the SISVER platform once patients were registered in the medical unit. Quality of data was verified at each step. Patients with suspicion of viral respiratory disease were isolated and combined pharyngeal/ nasopharyngeal (for all patients), bronchoalveolar lavage (for intubated patients only), or lung biopsy (for deceased patients only) samples were collected for reverse-transcriptase–polymerase-chain-reaction (RT-PCR) testing. RT-PCR testing was carried out by the National Network of Epidemiological Surveillance Laboratories, namely the Institute of Epidemiological Diagnosis and Reference “Dr Manuel Martínez Báez”, State Laboratory of Public Health, and Support Laboratories for Epidemiological Surveillance; results were recorded in the SISVER platform [13]. Finally, data were summarized and published by the General Directorate of Epidemiology.

A total of 65,535 patients with suspicion of viral respiratory disease from 01 January to 13 May 2020 were analysed; of them, 38,324 were SARS-CoV-2 positive and 27 211 were SARS-CoV-2 negative. Flowchart of participants is depicted in Online Resource 2. Distribution of patients according to the sectors of the National Health System and regions of Mexico involved in the study is illustrated in Online Resource 3.

Reverse-transcriptase–polymerase-chain-reaction testing

RT-PCR testing for detecting SARS-CoV-2 virus was carried out as stated by the guidelines of the Institute of Epidemiological Diagnosis and Reference “Dr Manuel Martínez Báez” and endorsed by the General Directorate of Epidemiology of Mexico [20]. All types of RT-PCR tests had an analytical sensitivity of ≤ 250 copies per reaction and analytical specificity of 100%. As recommended by the World Health Organization, RT-PCR assays targeting the RNA-dependent RNA polymerase, envelope, and nucleocapsid genes of SARS-CoV-2 virus were used for 93.59% of patients (analytical sensitivity: 10 copies per reaction; analytical specificity: 100%) [14]; other RT-PCR tests used by private laboratories (and listed in Online Resource 4) accounted for 5.45%; unspecified type of RT-PCR tests was 0.96%.

General variables, chronic comorbidities and outcomes

The following general data were examined: age, child/adult, biological sex, smoking habit, pregnancy, time from symptoms onset to medical contact, time from medical contact to death, and time from symptoms onset to death. Medical contact was defined as the date of the patient´s registration into the medical record. The analysed chronic comorbidities were the following: chronic obstructive pulmonary disease (COPD), bronchial asthma, hypertension, cardiovascular disease, diabetes mellitus, obesity, chronic kidney disease, and others. The adverse clinical outcomes of interest were hospital admission, pneumonia, intensive care unit (ICU) admission, endotracheal intubation, and mortality.

Statistical analysis

Bias in detecting SARS-CoV-2 virus infection was unlikely since RT-PCR was achieved for all patients and was carried out by validated methods [13, 14]. Study objective was unknown for healthcare physicians and patients, which allowed to minimize the following sources of biases: (1) selective report of chronic comorbidities and adverse clinical outcomes; and (2) selective mistake or unconscious omission in the report of adverse clinical outcomes. In this regard, the likelihood of mistake/ omission was the same for patients with and without chronic comorbidities. Finally, as national administrative healthcare open data might be another source of selection biases, we could expect that patients with SARS-CoV-2 infection were more likely to be reported; however, since patients with and without SARS-CoV-2 infection were separately analysed, this sort of bias was absolutely controlled.

The total of 38,324 SARS-CoV-2 positive and 27,211 SARS-CoV-2 negative patients reported in the national administrative healthcare open data were used for analysis. Variables registered in the dataset as “ignored” or “unspecified” were considered as missing data. As a small number of absent data was found (Online Resource 5), no treatment for missing data was carried out. Description of missing data was reported in the result section and not included in statistical analyses.

Categorical variables are shown as absolute numbers with percentage, whereas numerical variables are represented as means with standard deviation (SD) or medians with 25th – 75th intercuartil range (IQR) depending on whether the population was normally distributed (assessed by Kolmogorov–Smirnov test and Q-Q plot). Differences between groups were assessed using the chi-square test (χ2) for categorical variables; t-test or Mann–Whitney U test were used for continuous variables depending on whether the population was normally distributed or not, respectively.

Simple and multivariate logistic regression analyses were used to identify the factors associated with adverse clinical outcome. For each multivariate model, age, sex, smoking habit, time from symptoms onset to medical contact, and all the comorbidities were introduced. Classification tree with exhaustive CHAID method was used for categorizing numerical variables. Results were summarized as odds ratio (OR), 99% confidence interval (CI), and p value.

Two analyses were carried out. In the first, the total cohort of patients was examined. In the second, a subgroup analysis (for SARS-CoV-2 positive and for SARS-CoV-2 negative patients) was made.

Statistical tests with a two-tailed p value < 0.01 were considered as significant. Data were analysed using IBM®SPSS® Statistics 23.0 (IBM, Chicago, IL, USA). Microsoft Excel software was also used for graphs.

Results

General characteristics of patients

The mean age of the population studied was 44.0 years (SD 16.9 years), and 10,872 (16.6%) patients had > 60 years old. Male and female patients accounted for 53.2% and 46.8%, respectively. Smoking habit was identified in 6,029 (9.2%) patients. Pregnant women accounted for 1.0% of the total number of cases (2.2% of the female patients). Children comprised 4.4% of patients. The median time from symptoms onset to medical contact was 3.0 days (IQR 1.0–6.0 days). Hospital and ICU admission was required in 32.3% and 2.9% of patients, respectively. Pneumonia was diagnosed in 22.7%, endotracheal intubation was used in 2.6%, and overall mortality rate was 6.7%. The median time from symptoms onset to death and time from medical contact to death was 9.0 days (IQR 6.0–13.0 days) and 5.0 days (IQR 2.0–9.0 days), respectively. Table 1 illustrates the general characteristics of patients.

Table 1 General characteristics of participants

Chronic comorbidities

Chronic comorbidities were found in 28,668 patients (43.7%); of them, 7,780 (27.1%) had two comorbidities and 4,122 (14.4%) had three or more comorbidities. The prevalence of chronic comorbidities in both SARS-CoV-2 positive and negative patients was 46.9% and 39.4%, respectively. The most common diseases were hypertension, obesity and diabetes mellitus (Fig. 1a and b).

Fig. 1
figure 1

Distribution of chronic comorbidities in SARS-CoV-2 positive (a) and negative (b) patients. COPD chronic obstructive pulmonary disease

Among SARS-CoV-2 positive patients, age in patients with chronic comorbidities was higher than that in patients without comorbidities (52.2 vs. 42.0 years; p < 0.0001); additionally, differences in distribution of time from symptoms onset to medical contact (p = 0.003), time from medical contact to death (p < 0.0001) and time from symptoms onset to death (p < 0.0001) were observed between patients with and without chronic comorbidities. In SARS-CoV-2 negative group, patients with chronic comorbidities were older than those without chronic comorbidities (mean 47.6 vs. 34.8 years; p < 0.0001); however, time from symptoms onset to medical contact (p = 0.149), time from symptoms onset to death (p = 0.021) and time from medical contact to death (p = 0.184) in patients with chronic comorbidities were similar to those in patients without chronic comorbidities.

Factors associated with adverse clinical outcomes in the total cohort

Univariate analysis of factors associated with adverse clinical outcomes in the total cohort is depicted in Online Resource 6. Multivariate analyses demonstrated that age, sex male, time from symptoms onset to medical contact, numerous chronic comorbidities, and SARS-CoV-2 positive were independently associated with several adverse clinical outcomes (Table 2). Noteworthy, SARS-CoV-2 positive and age ≥ 80 years had the greatest size of effect.

Table 2 Multivariate logistic regression analyses of factors associated with adverse clinical outcomes in the total cohort

Factors associated with adverse clinical outcomes by subgroup analysis

The rates of adverse clinical outcomes were significantly higher in patients with chronic comorbidities than those in patients without chronic comorbidities for both SARS-CoV-2 positive and negative subgroups, as shown by univariate analysis (Online Resource 7, 8, 9, 10, 11). Multivariate logistic regression analysis demonstrated that age > 40 years, sex male, time from symptoms onset to medical contact > 3 days, and chronic comorbidities such as COPD, hypertension, diabetes mellitus, obesity and chronic kidney disease were independent risk factors for adverse clinical outcomes in SARS-CoV-2 positive patients (Table 3). For SARS-CoV-2 negative patients, age > 60 years, sex male and chronic comorbidities such as cardiovascular disease, diabetes mellitus and chronic kidney disease were the risk factors associated with an increased risk of adverse clinical outcomes (Table 4).

Table 3 Multivariate logistic regression analyses of factors associated with adverse clinical outcomes in the SARS-CoV-2 positive patients
Table 4 Multivariate logistic regression analyses of factors associated with adverse clinical outcomes in the SARS-CoV-2 negative patients

Discussion

Since December 2019 until present, more than 1 million people have died globally as a result of COVID-19 [15]; furthermore, COVID-19 pandemic has strained the medical systems worldwide and new outbreaks have been reported in several regions including Europe [15]. We found that SARS-CoV-2 infection and age ≥ 80 years had the greatest size of effect, demonstrating that COVID-19 and oldest patients had the most powerful association with worse clinical outcomes at the population level. These results support that the persistent cries of alarm is imperative because the world remains in danger.

Age was related to adverse clinical outcomes in both SARS-CoV-2 positive and negative patients, but the age at risk was lower for COVID-19 patients (40 vs. 60 years). Although age has been repeatedly related to adverse clinical outcomes in previous studies of COVID-19 [8, 16,17,18,19], the present results show that the strength of association could be greater than believed. Immune response to viral infections is reduced over time because of a limited cross-reactivity between anti-seasonal coronavirus and anti-SARS antibodies. In previously exposed individuals to seasonal coronaviruses but waning antibody titres (e.g. the elderly), the recall antibody production after SARS-CoV-2 infection can results in a pro-inflammatory state [20].

Among patients with SARS-CoV-2 infection, the time from symptoms onset to medical contact was found as a factor related to adverse clinical outcomes, highlighting the importance of an early medical seeking. General availability of health coverage and population-educating programs should be implemented in order to avoid delayed medical care. We found that sex male increased the risk of adverse clinical outcomes. Other studies reported similar results [18, 19, 21]. Immunity system is weaker in male than in female sex, therefore men are at higher risk of viral infection than women [22]. Because the low prevalence of smoking in patients hospitalized with COVID-19 [23, 24], the protective effect of tobacco is largely discussed [25]; however, several meta-analyses of observational studies have confirmed that smoking habit is related to severe COVID-19 and adverse outcomes [19, 21, 26,27,28]. Tobacco smoking is associated with upregulation of the angiotensin-converting-enzyme-2 (ACE2), immune system disturbance, endothelial injury and hypercoagulable state [28, 29]; these factors and pre-existing comorbidities could explain adverse outcomes observed in smokers with COVID-19. We found no effect of smoking habit on clinical outcomes, which could be partially due to the aforementioned low prevalence of smoking in COVID-19 patients, misclassification of smokers as nonsmokers, and the composition of the study participants: most previous studies have been focused on hospitalized patients while we carried out a population-based study including out-hospital patients. Our results are placed in the middle region of the risk spectrum: no protective or harmful effect, which increases the controversy on this topic. According to the currently contradictory knowledge, further large prospective studies are required to clarify the relationship between smoking habit and COVID-19.

We found the prevalence of respiratory diseases, hypertension, cardiovascular diseases, obesity and chronic kidney disease in COVID-19 patients was similar to those described in previous studies [3, 6, 8, 10]; nonetheless, the prevalence of diabetes mellitus was slightly higher than those reported by other authors [6, 9], which may be explained by the variation in the prevalence of chronic comorbidities among different countries. As previous studies in COVID-19 patients [4, 6, 8,9,10, 16,17,18], we observed a relationship between chronic comorbidities (e.g. COPD, hypertension, diabetes mellitus, obesity and chronic kidney disease) and adverse clinical outcomes. This association could be expected because patients with chronic diseases are usually aged, immunosenescent, or have impaired respiratory tract function. Particularly, COVID-19 pandemic have been challenged in countries with a high prevalence of obesity [30]; morbidity and mortality rates are higher in individuals with obesity and COVID-19 [31, 32]. Explanations for this association include the chronic pro-inflammatory state, the excessive oxidative stress response, and the impaired immunity; cytokines, mammalian target of rapamycin, and altered natural killer cell also play a central role [33]. Preventive weight-loss measures should be stimulated in obese patients.

Bronchial asthma was not associated with clinical outcomes in our analysis; other studies have shown similar results [34,35,36]. Type 2 immune response and its related cytokines such as IL-4 and IL-13, accumulation of eosinophils, and conventional therapeutics (e.g. inhaled corticosteroids, allergen immunotherapy, and anti-IgE monoclonal antibody) might reduce the risks of asthmatics suffering viral infection through alleviating inflammation or enhancing antiviral defence, which provides potential protective effects against COVID-19 [37]. Because the effect of bronchial asthma, diabetes mellitus and chronic kidney disease on clinical outcomes was similar for both SARS-CoV-2 positive and negative patients, it may represent the impact of these chronic comorbidities on the health status of general population; indeed, some meta-analyses have failed to demonstrate a pooled effect of chronic kidney disease [4] and diabetes mellitus on severity of COVID-19 [6]. Cardiovascular disease has been linked to poor clinical outcomes in COVID-19 patients [19, 21], but such association was not found in the present study. Potential reasons for this unexpected result are the study population (other researches have been centred in hospitalized or critically ill patients while we carried out a population-based study including out-hospital patients), the number and type of covariates included in multivariate analysis (we analysed the covariates with all chronic comorbidities, therefore all patients´ information was explored and taken as a whole, depicting the real clinical context), the definition of cardiovascular disease used in specific studies, and the random characteristics of the study population. Conversely, COVID-19 patients with COPD, hypertension, or obesity seem to have a real increased risk for adverse clinical outcomes. These results suggest a different impact of the specific type of chronic comorbidity on clinical outcomes between patients with and without SARS-CoV-2 infection, which requires further researches.

According to the current knowledge, the most important pathophysiological factors in COVID-19 seem to be the viral load and the patient´s immune response as described before [38, 39]; then chronic comorbidities might modify the magnitude of the response to viral aggression by amplifying the inflammatory pathways [40]. The loss of ACE2-mediated protection against the pro-inflammatory, vasoconstrictor and pro-fibrotic effects of angiotensin II [41]; endothelial dysfunction and pro-thrombotic state [42]; and single nucleotide polymorphisms in the host DNA might contribute to worse clinical outcomes among COVID-19 patients with chronic comorbidities [43].

COVID-19 patients with chronic underlying diseases are high-risk patients requiring strict adherence to preventive and treatment protocols. This approach supports those interventions aimed to reduce the viral load (e.g. preventive strategies including vaccination program, antiviral therapies, plasma from convalescent patients) and immunological reinforcement (e.g. type I interferon, innate immunity enhancement approaches). Since the COVID-19 pandemic remains as a real world threat [15], ICU practitioners, hospital administrators, governments, and policy makers must prepare for an increase in hospital and ICU beds requirements; quantifying the level of risk related to chronic underlying diseases is vital for these purposes. This information may help the healthcare systems to guide vulnerable populations, assess the risk of deterioration, organize the human and material health resources, and optimize the health cost.

There are several strengths in this study. First, this is an adequately powered study which includes a large number of patients. Second, this is a large population-based study assessing the influence of a variety of specific chronic underlying diseases on several clinical outcomes in patients with COVID-19, which depicts the real clinical context. Most previous studies used small sample sizes, single centre observations, heterogeneous data, and usually focused on just one clinical outcome [6,7,8,9,10,11]; these limitations may lead to biases and uncertain statistical estimations, which were overcome through this large population-based study. Third, we also studied the clinical impact of chronic comorbidities among patients without COVID-19, which has not been performed before. These data account for a holistic analysis about the real impact of chronic underlying diseases for healthcare systems in the current setting of respiratory diseases and COVID-19 pandemic.

Limitations of this study should be addressed. First, this is a retrospective analysis where chronic underlying diseases were self-reported. Underreport of chronic comorbidities could modify the estimation of the strength of association with clinical outcomes; however, the prevalence of chronic comorbidities was similar to those informed in other cohorts of COVID-19 patients [3, 6, 9, 10]. Second, some patients with pneumonia in the SARS-CoV-2 negative group might actually have a falsely negative RT-PCR in the upper respiratory tract sample. Bronchoalveolar lavage by bronchoscopy, a method with higher diagnostic yield [44], is an aerosol-generating procedure used in selected patients (e.g. intubated patients) in order to minimise exposure of healthcare workers to SARS-CoV-2 [45]. Nevertheless, considering the false negative rate in COVID-19 pneumonia reported in previous studies [46, 47] and the number of SARS-CoV-2 negative patients with pneumonia in our study, a little proportion of patients were misclassified (≈ 3.1%). According to the high sample size of the study, misclassification had slight effect on statistical estimations. Third, in nonCOVID-19 patients, the aetiology of the acute respiratory disease was not investigated. Patients with chronic comorbidities infected with nonSARS-CoV-2 respiratory viruses have a poorer prognosis [48, 49]. Fourth, since we used an open-data dataset, adverse clinical outcomes could develop in patients after updating the web page; however, patients with and without comorbidities had the same chance to develop outcomes, with little or no impact on estimations. Fifth, we only used the chronic comorbidities recorded in the administrative healthcare open data; therefore other chronic underlying diseases likely associated with adverse clinical outcomes were not explored; for instances, malignancies, cerebrovascular diseases and liver diseases have also linked to poorer outcome in COVID-19 patients [3,4,5, 8]. Finally, the study was carried out in Mexican population, which could be a limitation for generalizing results. Because some limitations could affect statistical estimations, it should be considered before extrapolating results to a specific population. Overcoming these limitations is required for designing future researches.

Conclusions

Chronic underlying diseases are commonly found in patients with and without COVID-19, which are significant risk factors for adverse clinical outcomes. This knowledge can better define those COVID-19 patients at higher risk, who require more personalized and specific approaches to prevent complications. The different impact of the specific type of chronic comorbidity on clinical outcomes in patients with and without SARS-CoV-2 infection requires further researches. Although these findings need confirmation using other data sources, it is a useful contribution in our understanding of COVID-19.