Skip to main content

Comparative epidemiology of five waves of COVID-19 in Mexico, March 2020–August 2022

Abstract

Background

The Mexican Institute of Social Security (IMSS) is the largest health care provider in Mexico, covering about 48% of the Mexican population. In this report, we describe the epidemiological patterns related to confirmed cases, hospitalizations, intubations, and in-hospital mortality due to COVID-19 and associated factors, during five epidemic waves recorded in the IMSS surveillance system.

Methods

We analyzed COVID-19 laboratory-confirmed cases from the Online Epidemiological Surveillance System (SINOLAVE) from March 29th, 2020, to August 27th, 2022. We constructed weekly epidemic curves describing temporal patterns of confirmed cases and hospitalizations by age, gender, and wave. We also estimated hospitalization, intubation, and hospital case fatality rates. The mean days of in-hospital stay and hospital admission delay were calculated across five pandemic waves. Logistic regression models were employed to assess the association between demographic factors, comorbidities, wave, and vaccination and the risk of severe disease and in-hospital death.

Results

A total of 3,396,375 laboratory-confirmed COVID-19 cases were recorded across the five waves. The introduction of rapid antigen testing at the end of 2020 increased detection and modified epidemiological estimates. Overall, 11% (95% CI 10.9, 11.1) of confirmed cases were hospitalized, 20.6% (95% CI 20.5, 20.7) of the hospitalized cases were intubated, and the hospital case fatality rate was 45.1% (95% CI 44.9, 45.3). The mean in-hospital stay was 9.11 days, and patients were admitted on average 5.07 days after symptoms onset. The most recent waves dominated by the Omicron variant had the highest incidence. Hospitalization, intubation, and mean hospitalization days decreased during subsequent waves. The in-hospital case fatality rate fluctuated across waves, reaching its highest value during the second wave in winter 2020. A notable decrease in hospitalization was observed primarily among individuals ≥ 60 years. The risk of severe disease and death was positively associated with comorbidities, age, and male gender; and declined with later waves and vaccination status.

Conclusion

During the five pandemic waves, we observed an increase in the number of cases and a reduction in severity metrics. During the first three waves, the high in-hospital fatality rate was associated with hospitalization practices for critical patients with comorbidities.

Peer Review reports

Background

SARS-CoV-2 has spread worldwide in the form of a series of waves of infection with varying levels of transmissibility and severity [1, 2] and modulated by multiple factors, including the intrinsic transmission properties of evolving variants, social distancing measures (e.g., face mask wearing), access to medical care and vaccination campaigns. Investigating the epidemiological factors that have shaped the multi-wave COVID-19 patterns of infection, hospitalization, and death at different spatial scales could improve public health policy to save lives and minimize economic disruptions for future pandemics.

The COVID-19 pandemic has disproportionately affected many Latin American countries, with 176,521,913 cases, including 2,822,535 deaths reported as of September 8th, 2022 [3]. Mexico, one of North America’s most highly populated countries with high poverty levels [4] and low COVID-19 testing rates, documented the first imported COVID-19 case on February 27th, 2020. Mexico currently ranks 7th in the world regarding COVID-19 reported deaths, with a total of 329,652 recorded deaths and 7,046,220 reported cases (4.68% of confirmed case fatality rate) as of September 8th, 2022 [5].

In this study, we investigated the epidemiological patterns of COVID-19 infections, hospitalizations, deaths, and factors associated with severe disease outcomes during the five epidemic waves in Mexico. To that end, we used data reported by the Mexican Institute of Social Security (IMSS) surveillance system, the most extensive Latin-American social security system, and Mexico’s leading health institution [6].

Methods

Study design

To understand the COVID-19 epidemic in Mexico, we analyzed epidemiological surveillance data from IMSS, a tripartite Mexican health system covering active and retired workers and their families from the private sector. This group comprises roughly 60 million people, 48% of the Mexican population (126 million individuals), with a network of 1530 primary health care units and 251 2nd care level hospitals, and 25 high specialization hospitals nationwide. Overall, the IMSS population is representative of the urban Mexican population in the low-medium deciles of income. However, it is not representative of rural areas and those populations in the lowest deciles of income [6].

Epidemiological data from the Online Epidemiological Surveillance System (SINOLAVE) was used to describe trends. SINOLAVE is an epidemiological surveillance platform implemented by the IMSS during the 2009 influenza pandemic in México to collect information about acute respiratory infection cases seeking care at IMSS health care facilities [7, 8]. For the COVID-19 pandemic, SINOLAVE was adapted for COVID-19 suspected patients that fulfill the clinical signs and symptoms of the national operational case definition (during a timeframe of 10 days, presence of at least one of the following significant symptoms: fever, cough, headache, or dyspnea; accompanied by at least one of the following minor symptoms: arthralgias, myalgias odynophagia, runny nose, conjunctivitis, chest pain, anosmia, dysgeusia, and chill).

For each case, SINOLAVE collects data relating to demographic characteristics (age, gender), underlying conditions (self-reported chronic diseases such as diabetes, hypertension, COPD, cardiovascular disease, among others), symptoms at the time of testing, COVID-19 laboratory results, COVID-19 vaccination, as well as dates of symptoms onset (self-reported). Regarding vaccination status, patients were classified into four groups: (1) unvaccinated, (2) with incomplete schema (with a single dose of vaccine AstraZeneca, Sputnik, Moderna, Novavax, Pfizer BioNTech, Sinopharma, or Sinovac), (3) with complete schema (with two doses of the previous vaccines or one dose of CanSino or Janssen vaccines) and (4) with booster schema if the subject has a complete schema and an additional dose of any COVID-19 vaccine.

From the beginning and throughout the COVID-19 epidemic, the Mexican surveillance system created by the Ministry of Health established mandatory testing for 10% of symptomatic outpatient cases and 100% of hospitalized patients and deaths. Additionally, IMSS tested symptomatic active workers to issue temporal total disability certificates. However, testing coverage varied as the epidemic evolved. During the first and second waves, PCR testing was prioritized for people at high risk of severe disease (defined by age and comorbidities). Testing was available for all symptomatic cases seeking care after the introduction of rapid tests at the end of December 2020.

SINOLAVE information was supplemented with hospitalization data from the Virtual Center in Emergencies and Disasters (CVOED) platform, using an IMSS unique identification number. CVOED was used to complement case records with hospitalization and intubation status information. The latter system collects information from health care settings equipped to care for COVID-19 hospitalizations, including isolation areas and mechanical ventilation. A COVID-19 hospitalization was registered as such on SINOLAVE or CVOED. Criteria for hospitalization were based on medical judgment. Hospitalized cases frequently present with comorbidities, pneumonia, low oxygenation saturation, or respiratory failure. Additionally, a score on the scale CURB-65 (Confusion, Urea level, Respiratory rate, Blood pressure, and age ≥ 65), Pneumonia Severity Index (PSI), National Early Warning Score (NEWS) or qSOFA helped guide hospitalization decisions [9]. Since space in intensive care units was limited, hospitals were adapted to care for critical patients in general hospitalization wards in designated COVID-19 hospitals. These areas were equipped with oxygen supply, mechanical ventilation, and a health care team for every 24 patients. The team was led by a medical specialist in internal medicine, pneumology, infectiology, intensive medicine, or anesthesiology with help from other specialists or general practitioners’ doctors and by a group of nurses. These designated hospitals grew their capacity as the demand increased and have maintained a varying number of beds over time.

We included laboratory-confirmed COVID-19 cases by RT-PCR or Rapid Antigen tests (RAT) with the onset of symptoms from March 29th, 2020, to August 27th, 2022. Testing coverage varied as the epidemic evolved. During the first and second waves, PCR testing was prioritized for people at high risk of severe disease (defined by age and comorbidities). Testing was available for all suspected cases after the introduction of rapid tests at the end of 2020.

Days of admission delay were defined as the number of days elapsed from symptoms onset to hospitalization. In-hospital days were defined as the number of days from hospital admission to discharge or death.

Additionally, to describe the temporal distribution of SARS-CoV-2 variants, we analyzed a subset of 18,839 genetically sequenced samples. This subset with variant information included cases with onset of symptoms from April 4th, 2021, to August 20, 2022.

Statistical analysis

Epidemic curves construction

Epidemic curves of the weekly rates of total confirmed and hospitalized cases were constructed from March 29th, 2020 (week 2020-14) to August 27th, 2022 (week 2022-34). We also displayed the weekly hospitalization percentage among confirmed cases, mechanical ventilation (intubated) cases among those hospitalized, and the in-hospital case fatality rate. Weekly periods were defined based on the dates of onset of symptoms. Epidemic waves correspond to the following onset of symptoms periods: the first wave from week 2020-14 to week 2020-40 (from March 29th, 2020 to October 3rd, 2020); the second wave from week 2020-41 until week 2021-21 (from October 4th, 2020 to May 29th, 2021); the third wave from week 2021-22 to week 2021-50 (from May 30th, 2021 until December 18th, 2021); the fourth wave from week 2021-51 to 2022-17 (from December 19th, 2021 to April 30th, 2022); and fifth wave from week 2022-18 to week 2022-34 (from May 1st, 2022 until August 27th, 2022).

To estimate confirmed incidence rates, we used the IMSS-affiliated population in June 2022 (n = 60,293,938) as the denominator (Additional file 1: Table S1).

Epidemic wave summary statistics

We analyzed the epidemiological characteristics across five epidemic waves occurring during summer and winter 2020, summer and winter 2021, and summer 2022, respectively (Additional file 1: Table S2). Hospitalization percentage among confirmed cases, intubation percentages among hospitalized, hospital case fatality rate, mean days of hospital admission delay, and mean days of in-hospital stay were calculated in each epidemic wave with a 95% confidence interval (95% CI). Stratified analyses by gender and age group (< 20, 20 to 39, 40 to 50, and 60 or more years old) were also carried out.

Factors associated with COVID-19 severe outcomes

To investigate the influence of hospitalization practices, healthcare resources-related factors, the vulnerability of the patients, and vaccination coverage on the in-hospital case fatality rate observed during the pandemic, we evaluated temporal patterns of PCR and RAT tests, the proportion of hospitalizations with comorbidities, and vaccination coverage. Multiple logistic regression models were used to assess the association of age, sex, comorbidities, vaccination status, mechanical ventilation, ICU admission, and indicator variables for the number of waves with the risk of hospitalization among all positive cases and the risk of death among hospitalized patients.

Statistical analyses were carried out using Stata version 14. Confidence intervals were calculated for exact proportions and means, assuming a normal distribution. To compare the distribution of variables among comparison groups, we used a t-test, one-way ANOVA, and Chi-square trend test.

Results

IMSS confirmed the first case at the end of March 2020. The epidemiological system reported symptomatic and laboratory-confirmed cases. The number of cases and estimates of hospitalization and in-hospital case fatality rates are presented in Table 1.

Table 1 Summary of COVID-19 cases and estimates by epidemic wave

The five waves comprised 8,783,565 suspected cases, of which 7,502,195 were tested (85.4%) and 3,396,375 laboratory-confirmed COVID-19 cases with onset of symptoms from March 29th, 2020, to August 27th, 2022 were analyzed (positivity rate of 45.3% among tested cases). Overall, 11% cases were hospitalized (373,646 of 3,396,375 confirmed cases), 20.6% were intubated (76,929 of 373,646 hospitalized cases), and 45.1% was the overall hospital case fatality rate (168,499 of 373,646 hospitalized cases). Below we describe the epidemiological patterns related to case incidence rates of confirmed cases, hospitalizations, in-hospital mortality, hospital admission delays, and in-hospital stays across waves, age, and gender.

Lab-confirmed COVID-19 cases

The number of weekly cases showed five peaks during the study period. The peak that occurred during January 2022 had the highest incidence, followed by the fifth peak in summer 2022. The first peak in 2020 had the lowest incidence (Figs. 1A and 2A).

Fig. 1
figure 1

COVID-19 epidemic curve and estimates. A Shows the number of COVID-19 confirmed cases (orange bars for total cases and dark green bars for hospitalized cases). The middle panel shows the number of hospitalization cases in dark green bars. B Contains the weekly trend of the following estimations: hospitalization percentage among confirmed cases (green line), intubation percentage among hospitalized (yellow line) and hospital case fatality rate (CFR) (red line). Figures include the onset of symptoms week period from 2020-14 (April 1st, 2020) to 2022-34 (August 27th, 2022). Dotted vertical lines represent the separation of the five epidemic waves. Epidemic waves correspond to the following onset of symptoms periods: the first wave from week 2020-14 to week 2020-40 (from March 29th, 2020 to October 3rd, 2020); the second wave from week 2020-41 until week 2021-21 (from October 4th, 2020 to May 29th, 2021); the third wave from week 2021-22 to week 2021-50 (from May 30th, 2021 until December 18th, 2021); the fourth wave from week 2021-51 to 2022-17 (from December 19th, 2021 to April 30th, 2022); and fifth wave from week 2022-18 to week 2022-34 (from May 1st, 2022 until August 27th, 2022)

Fig. 2
figure 2

COVID-19 estimates according to the five epidemic waves. The figure shows the following estimations with a 95% confidence interval: incidence rate among the population (orange bars), hospitalization rate among the population (dark green bars), hospitalization percentage among confirmed cases (green bars), hospital case fatality rate (red bars), mean hospital admission delay (purple bars) intubation percentage (yellow bars), and mean hospitalization days (blue bars). Epidemic waves correspond to the following onset of symptoms periods: the first wave from week 2020-14 to week 2020-40 (from March 29th, 2020 to October 3rd, 2020); the second wave from week 2020-41 until week 2021-21 (from October 4th, 2020 to May 29th, 2021); the third wave from week 2021-22 to week 2021-50 (from May 30th, 2021 until December 18th, 2021); the fourth wave from week 2021-51 to 2022-17 (from December 19th, 2021 to April 30th, 2022); and fifth wave from week 2022-18 to week 2022-34 (from May 1st, 2022 until August 27th, 2022)

Based on a subset of cases with variant information available since spring 2021, we found that the second, third, fourth, and fifth waves were dominated by the Alpha-Gamma, Delta, and Omicron variants, respectively (Additional file 1: Figure S1).

The multi-wave pattern was also evident when the epidemic curves were stratified by gender, and age groups, especially the curves for those younger than 60 years. Also, the incidence rate among males was higher during the first three waves than among females. In contrast, during the fourth and fifth waves, the incidence among females was higher than that for males. Across age groups, we found that individuals 20–39 and 40–59 years old showed the highest incidence rates during the last two waves and the first two epidemic waves, respectively. Individuals ≥ 60 years exhibited their highest incidence during the second wave in the winter of 2020 (Additional file 1: Tables S3 and S4, and Figures S2A, S3A, S4A, and S5A).

Hospitalization rates among the population

Regarding hospitalization, the second wave showed the highest rate, followed by the first wave. Gradually, the hospitalization rates declined from the third to fifth waves (Figs. 1 and 2). Compared to females, males showed higher hospitalization rates than females. Individuals ≥ 60 years exhibited higher hospitalization rates than younger groups having their highest value during the second wave in the winter of 2020 (Additional file 1: Figures S2A, S3B, S4A, and S5B).

Hospitalization among confirmed cases

Overall, 11% (95% CI 10.9%, 11.1%) of the cases were hospitalized. The weekly hospitalization percentage declined across waves from 39.9% (95% CI 39.7%, 40.1%) in the first wave to 1.8% (95% CI 1.7%, 1.8%) in the fifth wave (Figs. 1B and 2C). Males and older individuals (≥ 60 years) were more likely to require hospitalization than females and younger groups. The hospitalization percentage by gender and age groups decreased across subsequent waves. A substantial decrease in hospitalization percentage was observed among older individuals (≥ 60 years.) during the third and fourth waves (Fig. 3, Additional file 1: Tables S3 and S4, and Figures S2B, S3C, S4B, and S5C).

Fig. 3
figure 3

Severe COVID-19 outcomes according to age group. The figure shows the weekly trend of the following severe COVID-19 outcomes by age group: hospitalization percentage among confirmed cases (left panel), intubation percentage among hospitalized (middle panel) and hospital case fatality rate (right panel). Age groups are represented as follows: age below 20 years old in the green line, 20 to 39 years old in the blue line, 40 to 59 years in the orange line, and red line for the group aged 60 years old and over. Figures include the onset of symptoms from week 2020-14 to week 2022-34 (from April 1st, 2020, to August 27th, 2022). Dotted vertical lines represent the separation of the five epidemic waves. Epidemic waves correspond to the following onset of symptoms periods: the first wave from week 2020-14 to week 2020-40 (from March 29th, 2020 to October 3rd, 2020); the second wave from week 2020-41 until week 2021-21 (from October 4th, 2020 to May 29th, 2021); the third wave from week 2021-22 to week 2021-50 (from May 30th, 2021 until December 18th, 2021); the fourth wave from week 2021-51 to 2022-17 (from December 19th, 2021 to April 30th, 2022); and fifth wave from week 2022-18 to week 2022-34 (from May 1st, 2022 until August 27th, 2022)

Intubation among hospitalized cases

The overall Intubation percentage was 20.6% (95% CI 20.5, 20.7) during the study period. It decreased from 26.0% (95% CI 25.7%, 26.2%) during the first wave to 5.3% (95% CI 4.9%, 5.7%) in the fifth wave (Figs. 1B and 2D). Intubation among hospitalized cases remained stable at ~ 20% during the second and third waves. Similarly, those ≥ 60 years and males were more likely to undergo intubation. The intubation percentage declined over subsequent waves across age groups and for males and females (Fig. 3, Additional file 1: Tables S3 and S4, and Figures S2B, S3D, S4B, and S5D).

Hospital case fatality rates

The overall hospital case fatality rate was 45.1% (95% CI 44.9, 45.3) and fluctuated over time, having its highest value during the second winter 2020 epidemic wave (50.8%, 95% CI 50.6%, 51.1%) followed by a significant decrease during the third (43.6%, 95% CI 43.3%, 44.0%), fourth (34.7%, 95% CI 34.2%, 25.2%), and fifth waves (17.7%, 95% CI 17.1%, 18.4%). It correlated with confirmed cases' epidemic peaks (Figs. 1B and 2E).

The hospital case fatality rate was higher in males than females and increased with age. Individuals ≥ 60 years exhibited the highest case fatality rate (59.3%), and those aged < 20 years showed the lowest hospital case fatality rate (Fig. 3, Additional file 1: Tables S3 and S4, and Figures S2B, S3E, S4B, and S5E).

Hospital admission delay and length of stay

The overall hospital admission delay was 5.07 (95% CI 5.06, 5.09) days, and the mean in-hospital stay was 9.11 (95% CI 9.09, 9.13) days. While the mean in-hospital stay showed a sustained decrease over time, the hospital admission delay increased from the first wave to the third wave and decreased in the last wave (Fig. 2F and G). Males had a slightly longer hospital admission delay than females (5.22 vs. 4.87 days, t-test − 24.7, P < 0.001). Similarly, the in-hospital stay tended to be longer among males than females (9.26 vs. 8.91 days, t-test − 14.8 P < 0.001). The age group 40 to 59 years old showed the highest mean of these metrics (Additional file 1: Tables S3 and S4, and Figures S3F, G, S5F, and G). Compared with survival patients, subjects that died had lower in-hospital stay (8.1 vs. 9.9 days; t-test = 75.1, P < 0.001) and slightly higher hospital admission delay (5.2 vs. 4.9 days; t-test = 27.3, P < 0.001).

Factors associated with COVID-19 severe outcomes

During the COVID-19 epidemic, we observed a decrease in the frequency of previous medical conditions among confirmed cases. We also observed an increase in RAT during the second wave and in COVID-19 vaccination from the third to the fifth wave (Additional file 1: Figure S6).

The number of persons who reported vaccination before the onset of illness is shown in Table 2. Ambulatory confirmed cases reported previous complete vaccination in 20.9%, while 5.4% of the hospitalized cases reported being fully vaccinated. Among those hospitalized, survivors reported 5.8% being fully vaccinated, while 4.9% of non-survivors reported being fully vaccinated.

Table 2 COVID-19 vaccination according to the hospitalization and in-hospital death status

In a multivariate logistic regression model, the risk of hospitalization among confirmed cases was positively associated with comorbidities, age, and male gender; and inversely associated with wave progression and complete vaccination before the onset of symptoms. The risk of death during hospitalization was reduced for those patients previously vaccinated, the fourth and fifth waves, and those admitted to an intensive care unit. As expected, this risk was increased for men and patients with comorbidities and was positively associated with age and intubation (Table 3).

Table 3 Multivariate logistic models for severe COVID-19 outcomes

Discussion

The COVID-19 pandemic has been characterized by a multi-wave pattern driven by multiple factors, including social control measurements, SARS-CoV-2 variants, testing practices, vaccination, and health care infrastructure. In this report, we described the epidemiological patterns of COVID-19 during five epidemic waves from surveillance data reported by IMSS, the most extensive Mexican and Latin-American social security system. Roughly, we observed an increase in reported incidence and a decrease in the severity of COVID-19 confirmed cases over time. Additionally, based on a subset of samples, we found that the second, third, fourth, and fifth waves were dominated by the Alpha-Gamma, Delta, and Omicron SARS-CoV-2 variants, respectively. During the first three waves, Mexico experienced a high in-hospital mortality rate associated with hospitalization criteria for critical patients with comorbidities. Further research is needed to improve case management during future pandemics.

Although surveillance systems do not represent all disease cases, they provide time patterns and valuable information for policy decision-making. Confirmed cases represent only a sample of all cases because not all ill persons seek medical care and have a specimen collected [10]. Additionally, the test’s sensitivity and specificity depend on the collection timing and the specimen’s quality. This under-ascertainment of COVID-19 ill persons has been evaluated with seroprevalence and excess mortality studies [11,12,13]. A seroprevalence survey just before vaccination began in December 2020 estimated the extent of SARS-CoV-2 infections in Mexico, regardless of disease severity and test availability, at as much as 25%, from a national representative sample [11]. The excess mortality percentage in Mexico has been estimated at 32% across the 2 pandemic years, and the number of death certificates associated with COVID-19 has been ~ 1.5 times the number of laboratory-confirmed deaths reported in the epidemiologic surveillance system [12].

In this case, trends of positive cases are representative of testing recommendations, adequately represent critically ill patients, and underrepresent ambulatory symptomatic cases. Although the recommendations for testing were similar during the five waves, the introduction of rapid testing presented a differential practice after the second wave (Additional file 1: Figure S7). Therefore, the magnitude of the first and second waves could be underestimated, and the comparison with subsequent waves may be affected by under ascertainment bias.

During the first wave in Mexico, the hospitalization percentage among 264,792 confirmed cases was 39.9% and declined to 1.8% in the fifth wave. The first figure was similar to the hospitalization percentage of 31.9% reported in March–August 2020 in the United States [14]. We also observed a higher hospitalization percentage in males and older people. In this respect, the Palmer et al. study reported that the exponential increase in COVID-19 hospitalization rates with age could be associated inversely with thymic T-cell production, which declines in males and older age [15]. In line with our results, a Denmark study observed a reduction in hospitalization rates for Omicron cases compared with Delta variant cases [16]. Our higher hospitalization rates compared with this previous study could also be explained by the Mexican testing strategy, which has been focused on severe COVID-19 cases.

Our in-hospital case fatality rate estimates were higher than other reports worldwide [17] but similar to those reported for critically ill hospitalized patients [18, 19]. A combination of factors could explain these elevated in-hospital case fatality rate estimates. These include the hospitalization practices characterized by strict admission criteria, admitting mostly critically ill cases, healthcare resources-related factors, the high vulnerability of the sick population, and the low COVID-19 vaccination coverage during the first waves (Additional file 1: Figures S6, S8, and S9). The contribution of organizational factors, the low proportion of intubated patients during hospitalization, and the reduced number of patients admitted to ICU deserve further research.

The multivariate logistic model results confirm the independent contribution of gender, age, previous medical conditions, vaccination status, and the reduction of disease severity over time. Similar to results from Milan, Italy, reported by Giacomelli et al., the hospital case fatality rate had its highest value during the second wave [20]. De la Peña reported that the case fatality rate increased with diabetes, hypertension, obesity, immunosuppression, and end-stage kidney disease. The population attributable fraction due to obesity in outpatients was 16.8% [21]. In another study, Hong et al. [22] reported that Mexico was among the countries with the highest CFRs. They also mentioned that medical resources might be important in preventing COVID-19-related deaths; however, considering the slight variation in fatality among the elderly, preventive measures such as vaccination are more critical, especially for the elderly population, to minimize the mortality rates. Olivas-Martinez et al. [23] found that the mortality rate over time was related to the availability of ICU beds, indirectly suggesting that overcrowding contributed to hospital mortality. In contrast with other studies that report an increase in COVID-19 severity with delta strain [24, 25], our study found that the third wave, mainly driven by this variant, had a lower in-hospital case fatality rate than the first two waves. These differences could be explained by the abovementioned hospitalization criteria and testing practices during the first two waves and the COVID-19 vaccination strategy that initially targeted the most vulnerable groups, leading to a decline in hospitalization rates for high-risk groups.

As expected, we found an important contribution of vaccination in preventing severe disease and reducing the risk of hospitalization. The OR for hospitalization among confirmed cases was 0.41 (95% CI 0.40, 0.41) compared to unvaccinated or incomplete vaccination. The impact of vaccination on the population level requires additional information for the general population. This information was not available in this study, restricted to laboratory-confirmed cases. However, the impact of vaccination on the COVID-19 pandemic has recently been reported with promising results [26,27,28].

A strength of our study is the large sample size, including records for nearly 3.4 million laboratory-confirmed COVID-19 cases during a 2-year pandemic period stemming from a population of 60 million individuals, representing the Mexican Population. Nevertheless, our findings need to be contextualized, considering some limitations. The analyzed data derives from passive database systems, which collect information from patients that seek and receive medical care at IMSS health care facilities. In this respect, milder COVID-19 cases would not be fully represented. Additionally, specific information about administering specific therapeutics before or during hospitalization was unavailable. With our data, we could not distinguish between hospitalization due to COVID-19 and hospitalizations for reasons other than COVID-19. Moreover, only vaccination coverage at the Mexican population level was available, with no information on specific coverage for IMSS affiliates.

Conclusions

In summary, our findings indicate that reported incidence tended to increase during the five COVID-19 epidemic waves experienced by the Mexican population. Still, severity tended to decline over subsequent waves, with an important effect of vaccination. The high in-hospital mortality rate was associated with hospitalization criteria for critical patients with comorbidities during the first three waves.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author upon reasonable request.

Abbreviations

IMSS:

Mexican Institute of Social Security

IRR:

Incidence rate ratio

COPD:

Chronic obstructive pulmonary disease

COVID-19:

Coronavirus disease 2019

CVOED:

Virtual Center in Emergencies and Disasters

RR:

Relative risk

RT-PCR:

Real-Time Reverse transcription polymerase chain reaction

SARS-CoV-2:

Severe acute respiratory syndrome coronavirus 2

SINOLAVE:

Online Epidemiological Surveillance System

References

  1. World Health Organization. Coronavirus disease (COVID-19) situation reports. https://www.who.int/emergencies/diseases/novel-coronavirus-2019/situation-reports. Accessed 8th Sept 2022.

  2. Dong E, Du H, Gardner L. An interactive web-based dashboard to track COVID-19 in real-time. Lancet Infect Dis. 2020;20(5):533–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Pan American Health Organization. Geo-Hub COVID-19-information system for the region of the Americas. https://paho-covid19-response-who.hub.arcgis.com/. Accessed 8th Sept 2022.

  4. World Bank. Atlas of sustainable development goals 2018: from world development indicators. https://openknowledge.worldbank.org/handle/10986/29788. Accessed 8th Sept 2022.

  5. Our World in Data. Total confirmed deaths due to COVID-19. https://ourworldindata.org/grapher/cumulative-deaths-and-cases-covid-19. Accessed 8th Sept 2022.

  6. Instituto Mexicano del Seguro Social. Informe al Ejecutivo Federal y al Congreso de la Unión sobre la Situación Financiera y los Riesgos del Instituto Mexicano del Seguro Social 2020–2021. https://www.imss.gob.mx/conoce-al-imss/informe-2020-2021. Accessed 8th Sept 2022.

  7. Echevarria-Zuno S, Mejia-Arangure JM, Mar-Obeso AJ, Grajales-Muniz C, Robles-Perez E, Gonzalez-Leon M, Ortega-Alvarez MC, Gonzalez-Bonilla C, Rascon-Pacheco RA, Borja-Aburto VH. Infection and death from influenza A H1N1 virus in Mexico: a retrospective analysis. Lancet. 2009;374(9707):2072–9.

    Article  PubMed  Google Scholar 

  8. Chowell G, Echevarria-Zuno S, Viboud C, Simonsen L, Tamerius J, Miller MA, Borja-Aburto VH. Characterizing the epidemiology of the 2009 influenza A/H1N1 pandemic in Mexico. PLoS Med. 2011;8(5):e1000436.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Secretaria de Salud. Lineamiento para la atención de pacientes por COVID-19. http://cvoed.imss.gob.mx/lineamiento-para-la-atencion-de-pacientes-por-covid-19/. Accessed 8th Sept 2022.

  10. Reed C, Angulo FJ, Swerdlow DL, Lipsitch M, Meltzer MI, Jernigan D, Finelli L. Estimates of the prevalence of pandemic (H1N1) 2009, United States, April-July 2009. Emerg Infect Dis. 2009;15(12):2004–7.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Basto-Abreu A, Carnalla M, Torres-Ibarra L, Romero-Martinez M, Martinez-Barnetche J, Lopez-Martinez I, Aparicio-Antonio R, Shamah-Levy T, Alpuche-Aranda C, Rivera JA, et al. Nationally representative SARS-CoV-2 antibody prevalence estimates after the first epidemic wave in Mexico. Nat Commun. 2022;13(1):589.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Gobierno de México. Exceso de Mortalidad en México. https://coronavirus.gob.mx/exceso-de-mortalidad-en-mexico/. Accessed 8th Sept 2022.

  13. Dahal S, Luo R, Swahn MH, Chowell G. Geospatial variability in excess death rates during the COVID-19 pandemic in Mexico: examining socio demographic, climate and population health characteristics. Int J Infect Dis. 2021;113:347–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Vahidy FS, Pan AP, Ahnstedt H, Munshi Y, Choi HA, Tiruneh Y, Nasir K, Kash BA, Andrieni JD, McCullough LD. Sex differences in susceptibility, severity, and outcomes of coronavirus disease 2019: cross-sectional analysis from a diverse US metropolitan area. PLoS ONE. 2021;16(1):e0245556.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Palmer S, Cunniffe N, Donnelly R. COVID-19 hospitalization rates rise exponentially with age, inversely proportional to thymic T-cell production. J R Soc Interface. 2021;18(176):20200982.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Bager P, Wohlfahrt J, Bhatt S, Stegger M, Legarth R, Moller CH, Skov RL, Valentiner-Branth P, Voldstedlund M, Fischer TK, et al. Risk of hospitalisation associated with infection with SARS-CoV-2 omicron variant versus delta variant in Denmark: an observational cohort study. Lancet Infect Dis. 2022;22(7):967–976.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Zimmermann IR, Sanchez MN, Frio GS, Alves LC, Pereira CCA, Lima RTS, Machado C, Santos LMP, Silva END. Trends in COVID-19 case-fatality rates in Brazilian public hospitals: a longitudinal cohort of 398,063 hospital admissions from 1st March to 3rd October 2020. PLoS ONE. 2021;16(7):e0254633.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Lim ZJ, Subramaniam A, Ponnapa Reddy M, Blecher G, Kadam U, Afroz A, Billah B, Ashwin S, Kubicki M, Bilotta F, et al. Case fatality rates for patients with COVID-19 requiring invasive mechanical ventilation. A meta-analysis. Am J Respir Crit Care Med. 2021;203(1):54–66.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Macedo A, Goncalves N, Febra C. COVID-19 fatality rates in hospitalized patients: systematic review and meta-analysis. Ann Epidemiol. 2021;57:14–21.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Giacomelli A, Ridolfo AL, Pezzati L, Oreni L, Carrozzo G, Beltrami M, Poloni A, Caloni B, Lazzarin S, Colombo M, et al. Mortality rates among COVID-19 patients hospitalised during the first three waves of the epidemic in Milan, Italy: a prospective observational study. PLoS ONE. 2022;17(4):e0263548.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Pena JE, Rascon-Pacheco RA, Ascencio-Montiel IJ, Gonzalez-Figueroa E, Fernandez-Garate JE, Medina-Gomez OS, Borja-Bustamante P, Santillan-Oropeza JA, Borja-Aburto VH. Hypertension, diabetes and obesity, major risk factors for death in patients with COVID-19 in Mexico. Arch Med Res. 2021;52(4):443–9.

    Article  CAS  PubMed  Google Scholar 

  22. Hong D, Lee S, Choi YJ, Moon S, Jang Y, Cho YM, Lee H, Min S, Park H, Hahn S, et al. The age-standardized incidence, mortality, and case fatality rates of COVID-19 in 79 countries: a cross-sectional comparison and their correlations with associated factors. Epidemiol Health. 2021;43:e2021061.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Olivas-Martinez A, Cardenas-Fragoso JL, Jimenez JV, Lozano-Cruz OA, Ortiz-Brizuela E, Tovar-Mendez VH, Medrano-Borromeo C, Martinez-Valenzuela A, Roman-Montes CM, Martinez-Guerra B, et al. In-hospital mortality from severe COVID-19 in a tertiary care center in Mexico City; causes of death, risk factors and the impact of hospital saturation. PLoS ONE. 2021;16(2):e0245772.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Taylor CA, Patel K, Pham H, Whitaker M, Anglin O, Kambhampati AK, Milucky J, Chai SJ, Kirley PD, Alden NB, et al. Severity of disease among adults hospitalized with laboratory-confirmed COVID-19 before and during the period of SARS-CoV-2 B.1.617.2 (Delta) Predominance—COVID-NET, 14 States, January–August 2021. MMWR Morb Mortal Wkly Rep. 2021;70(43):1513–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Twohig KA, Nyberg T, Zaidi A, Thelwall S, Sinnathamby MA, Aliabadi S, Seaman SR, Harris RJ, Hope R, Lopez-Bernal J, et al. Hospital admission and emergency care attendance risk for SARS-CoV-2 delta (B.1.617.2) compared with alpha (B.1.1.7) variants of concern: a cohort study. Lancet Infect Dis. 2022;22(1):35–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. McNamara LA, Wiegand RE, Burke RM, Sharma AJ, Sheppard M, Adjemian J, Ahmad FB, Anderson RN, Barbour KE, Binder AM, et al. Estimating the early impact of the US COVID-19 vaccination programme on COVID-19 cases, emergency department visits, hospital admissions, and deaths among adults aged 65 years and older: an ecological analysis of national surveillance data. Lancet. 2022;399(10320):152–60.

    Article  CAS  PubMed  Google Scholar 

  27. Vilches TN, Moghadas SM, Sah P, Fitzpatrick MC, Shoukat A, Pandey A, Galvani AP. Estimating COVID-19 infections, hospitalizations, and deaths following the US vaccination campaigns during the pandemic. JAMA Netw Open. 2022;5(1):e2142725.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Suthar AB, Wang J, Seffren V, Wiegand RE, Griffing S, Zell E. Public health impact of covid-19 vaccines in the US: observational study. BMJ. 2022;377:e069317.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

IMSS-affiliated authors were supported by intramural funding. GC is partially supported by NSF Grants 1610429 and 1633381 and R01 GM 130900.

Author information

Authors and Affiliations

Authors

Contributions

GC conceived the study question. IJAM carried out statistical analyses and interpreted the data. RARP, ODOL contributed to the analysis plan. VHBA, IJAM and GC drafted and reviewed the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Ramón Alberto Rascón-Pacheco.

Ethics declarations

Ethics approval and consent to participate

Our study followed ethical standards, institutional regulations, national laws, and the Helsinki declaration. This study was based on available databases with no personal identifiers. The Scientific Research National Committee of Mexican Institute of Social Security (IMSS) approved the research with the register number R-2020-785-165. As the nature of the study is retrospective, the Scientific Research National Committee of Mexican Institute of Social Security (IMSS) has waived informed consent for the study participants.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Additional Tables S1–S4 and Figures S1–S9.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ascencio-Montiel, I.J., Ovalle-Luna, O.D., Rascón-Pacheco, R.A. et al. Comparative epidemiology of five waves of COVID-19 in Mexico, March 2020–August 2022. BMC Infect Dis 22, 813 (2022). https://doi.org/10.1186/s12879-022-07800-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12879-022-07800-w

Keywords