Impact of environmental and individual factors on COVID-19 mortality in children and adolescents in Mexico: An observational study.

BACKGROUND: During the Covid-19 pandemic, children and adolescents faced poverty, potentially dying from preventable causes, or missing out essential vaccines. The aim of this study was to assess potential environmental and individual factors associated to COVID-19 mortality in children and adolescents in Mexico.
METHODS: This cross-sectional study analysed the official data of 131,001 children under 10 years of age and adolescents between 10 and 19 years of age with COVID-19 disease, in Mexico. Participants were diagnosticated between March 2020 and June 13, 2021. The environmental variables such as malnutrition, vaccination coverage and social lag index were evaluated at the state level. Multilevel models were fitted to evaluate the association between environmental and individual factors and COVID-19 mortality.
FINDINGS: A total of 773 (0.6%) children and adolescents died due to COVID-19. Younger age (OR = 0.878, 95%CI: 0.869-0.888), diabetes (OR = 3.898, 95%CI: 2.596-5.851), immunosuppression (OR = 5.410, 95%CI: 4.088-7.158), obesity (OR = 1.876, 95%CI: 1.397-2.521), hypertension (OR = 1.906, 95%CI: 1.239-2.932), cardiovascular disease (OR = 2.288, 95%CI: 1.482-3.531), and chronic kidney disease (OR = 13.250, 95%CI: 9.066-19.350) were associated with mortality. COVID-19 mortality was directly associated with social lag index and malnutrition (ORvery high = 2.939, 95%CI: 1.111-7.775, and OR = 1.390, 95%CI: 1.073-1.802, respectively), and inversely associated with population density (OR = 0.374, 95%CI: 0.204-0.688). Finally, children and adolescents living in areas with a higher percentage of people with incomplete education (OR = 1.045, 95%CI: 1.011-1.081), of children of school age of 6-14 years who do not attend school (OR = 1.266, 95%CI: 1.032-1.554), and of illiterate population aged 15 and over (OR = 1.086, 95%CI: 0.999-1.179) were associated with a higher risk of COVID-19 mortality.
INTERPRETATION: Malnutrition, social lag index and population density are key factors to understand COVID-19 mortality in children and adolescents. Also, age and pre-existing comorbidities were also associated with worse COVID-19 prognosis. FUNDING: No funding was secured for this study.

Evidence before this study

Preliminary evidence suggests the relevance of health inequities in COVID-19 outcomes. In Mexico, COVID-19 mortality in children and adolescents was higher compared to high-income countries. Malnutrition and social lag index have been associated with worse health outcomes in other diseases. Previous studies evaluated the impact of environmental factors on COVID-19 mortality in older adults or the general population. Also, it is known that children younger than four years of age with cardiovascular risk or immunosuppressed had an increased risk of mortality from COVID-19. However, no data are reported on the association between environmental and individual factors and COVID-19 mortality in children and adolescents.

Added value of this study

We analyzed COVID-19 mortality and risk factors in 131,001 children and adolescents with COVID-19 in Mexico. To our best knowledge, this is the largest study of COVID-19 in children and adolescents in Mexico aimed to evaluate environmental and individual factors associated with a worse COVID-19 prognosis. This study confirms the impact of factors such as malnutrition, social lag index or density population on COVID-19 mortality.

Implications of all the available evidence

Our results seem to show that COVID-19 is a mild disease in most of the paediatric cases. However, it is recommendable to consider the impact of malnutrition and social lag on COVID-19 outcomes, specially, in a low-middle-income country. Our findings point to the need to consider environmental factors to face future challenges regarding the pandemic. Likewise, this evidence on the impact of environmental characteristics on mortality in children and adolescents should contribute to rethinking public policies that improve the health of children and adolescents, and even help prepare for new pandemics.

Alt-text: Unlabelled box

Introduction

Less than 10% of COVID-19 cases involve children up to the age of 16.1, 2, 3 Although the clinical course of the symptomatic SARS-CoV-2 infection (COVID-19) is usually mild among young patients, children affected by the disease can suffer a moderate-serious infection. , The main COVID-19 complication in children is the multisystem inflammatory syndrome (MIS). COVID-19 mortality in young people is low, varing between 0 and 13%, however, recent information has revealed that the impact of this disease in pediatric populations is likely to vary between and within countries.

The World Health Organization (WHO) highlighted the association between a high risk of COVID-19 infection and poverty. Factors such as reduced access to healthcare, malnutrition and environmental and living conditions can increase the risk of a worse outcome. Differences in health outcomes according to income increased during the pandemic. An estimated 47 million children younger than 5 years old were suffering malnutrition previous to the COVID-19 pandemic. This condition was associated with prolonged hospitalization due to COVID-19 infection. Moreover, patients with COVID‐19 disease are prone to develop significant weight loss, malnutrition, and cachexia and, therefore, malnutrition could also be a prognostic factor for COVID-19–related increased mortality.

Comorbidities are associated with worst outcomes in pediatric COVID-19 patients. Obesity is associated with a worsened prognosis of the infection in both adults and children. Although diabetes, hypertension, malignancies and chronic respiratory disease were previously shown to be associated with COVID-19 prognosis in adults, evidence on the role of these comorbidities in children is still scarce. ,

This analysis evaluated the impact of environmental and individual factors on COVID-19 mortality in Mexican children and adolescents during the first two years of the pandemic. We hypothesized that children and adolescents living in areas with less advantageous conditions would experience worse COVID-19 outcomes.

Methods

This is a retrospective analysis of the open COVID-19 registry dataset collected by the Mexican government in its official coronavirus web page and updated daily (available at https://www.gob.mx/salud/documentos/datos-abiertos-152127). The study sample included children (0–9 years old) and adolescents (10–19 years old) with COVID-19 diagnosis from February 27, 2020 to June 13, 2021 (Figure 1 ). The WHO defines adolescence as the phase of life between childhood and adulthood, from ages 10–19; early adolescents are those aged 10–14 years, and late adolescents are those aged 15–19 years. COVID-19 diagnosis was made according to the Diagnosis and Treatment Guideline for COVID-19 published by the Mexican Secretary of Health (available at https://coronavirus.gob.mx/). A positive COVID-19 case was considered if the child or adolescent presented a positive RT-PCR test, a positive rapid Ag-T test, or a diagnosis defined by patients with a clinical-epidemiological relation compatible with SARS-CoV-2 infection, as well as deceased persons that were assigned the diagnosis by a specialized committee. The average time between onset of symptoms and death from COVID-19 was 13 days.

Figure 1

Flow diagram for study participants.

We used a harmonized dataset of individual and state-level data for 32 states (Mexico is a federal republic composed of 32 states). Based on the patient's place of residence, we linked individual level with state-level data using unique state codes.

Outcomes

The dependent variable was death in patients with COVID-19, which was assessed using date that the person died (alive=0; dead=1).

Clinical covariates

For this study, socio-demographic information on the patient was used, such as age, sex, state of residence, and ethnicity (indigenous/not indigenous). Clinical information included the presence of pneumonia, type of patient care (outpatient vs. inpatient), admission at intensive care unit (ICU), and treatment with intubation. Pre-existing comorbidities were determined based on the patient's self-report on admission, and included obesity, diabetes, hypertension, chronic obstructive pulmonary disease (COPD), asthma, immunosuppression, cardiovascular disease, chronic kidney disease (CKD), and other comorbidities.

Exposure variables of the environment

The influence of environmental factors on mortality in children and adolescents with COVID-19 was assessed according to different data sources. Exposure variables were all measured at the entity level.

Social lag index

The 2020 social lag index is an indicator that measures social development in Mexico, which is calculated by the National Council for the Evaluation of Social Development Policy (CONEVAL) for each federal entity. The index contains 11 indicators:  the percentage of the population (1) aged 15 and over that is illiterate, (2) with children of school age of 6–14 years who do not attend school, (3) aged 15 and over with incomplete basic education, and (4) not affiliated to any health service, and (5–11) percentage of inhabited private housing units without the following basic infrastructure: drainage (5), floor covering (6), toilets or sanitary services (7), piped water from the public network (8), electricity services (9), washing machine (10), or refrigerator (11). The social lag index was categorized as very low, low, medium, high and very high.

Malnutrition

The number of people with malnutrition by area in 2020 was obtained from the Mexican Secretary of Health. These persons were identified using the International Classification of Diseases, Tenth Revision (ICD-10) diagnosis codes E40-E46. ,

Vaccination

Information on vaccination coverage was captured from the National Survey on Health and Nutrition (ENSANUT) 2018. Two categories of vaccination schedule were analyzed as “complete schedule” and “schedule with 4 vaccines” for children of one year of age and up to two years old. A vaccination scheme was considered complete if there was a record that the child received a dose of BCG, three doses of hepatitis B, three doses of pentavalent, two doses of pneumococcal vaccine, two doses of anti-rotavirus, and one dose of MMR (triple viral). The 4-vaccine schedule completely excludes rotavirus and pneumococcal vaccines.

Population density

This variable was extracted from the National Institute of Statistics and Geography (INEGI, for its Spanish acronym) (inhabitants per km2).

Statistical analysis

Individual and environmental characteristics according to whether the patient was alive or dead were described, and means for continuous variables and percentages for categorical variables were used. Bivariate analyses were conducted to determine the association between individual and environmental variables and COVID-19 mortality. Finally, multilevel multivariable logistic regression models were performed with individuals nested within states, to test the association between each variable and COVID-19 mortality.

Models were adjusted for sex, age, diabetes, COPD, asthma, immunosuppression, hypertension, cardiovascular disease and obesity. Furthermore, the model with the malnutrition variable was also adjusted for population density, since this variable is in absolute numbers. Therefore, a set of adjusted models were fitted separately for each environmental variable, adjusting for individual covariates. Also, a model was fitted for each indicator related to social lag index. Models where the environmental variables were associated significantly with mortality at the 0.05 level are shown.

The overall model is showed bellow

∼ Binomial(1, )

Where πij is the probability of dying for individual i in state j, Xi is the set of explanatory variables at the individual level (sex, age, diabetes, obesity, immunosuppression, hypertension, cardiovascular disease and CKD), Zj is the set of explanatory variables defined for the states (social lag index or its indicators, vaccination coverage in children, people with malnutrition or population density), and uj are the residuals of level 2, for which it is assumed that they are independent and follow a normal distribution with mean 0 and variance .

To quantify the random effect, the Median Odds Ratio (MOR) was used, which can be interpreted as the increase in the median risk of worse control of a patient if he/she changes from one state to another with a higher risk. All statistical analyses were carried out using STATA version 16.0 (StataCorp, Stata Statistical Software, 2019). P-value < 0.05 was considered as statistically significant.

Ethics approval and consent to participate

The conduct of this study did not require the approval of an institutional ethics committee, because it is a secondary analysis of public data available on the platform of the Mexican Secretary of Health under the Mexico's open government data (OGD) policy, in order to facilitate the access, use, reuse and redistribution of information.

Role of the funding source

This research received no external funding. All authors had full accesses to all data in the study and had final responsibility for the decision to submit for publication.

Results

Table 1 summarizes characteristics of 131,001 children and adolescents with COVID-19 by outcome. A total of 773 (0.60%) patients died due to COVID-19: 421 were male (54.46%) and 352 female (45.54%), 419 (54.20%) patients were less than 10 years old. The mean age at the moment of COVID-19 diagnosis was 13.5 years (13.6 years for alive cases versus 8.2 years for deceased cases). Among indigenous children with COVID-19, observed mortality was 2.25% (17/757).] Also, 905 (0.69%) patients required ICU admission and 210 (27.17%) of those patients died. Comorbidities were more frequent in patients who died: 10.31% immunosuppression, 7.76% hypertension, 7.28% CKD, 6.64% diabetes, 5.07% cardiovascular disease, 3.99% obesity, 1.56% asthma, and 0.52% COPD. A total of 127 (16.43%) deceased cases presented more than one comorbidity, and the percentage of deceased patients with pneumonia was 66.36%.

Table 1

Characteristics of children and adolescents according to COVID-19 outcome.

		Alive (n = 130,228)	Dead (n = 773)	Total (n = 131,001)	p-value
Sex, n (%)	Female	65551 (50.34%)	352 (45.54%)	65903 (50.31%)	0.007
	Male	64677 (49.66%)	421 (54.46%)	65098 (49.69%)
Edad, mean (SD)	13.56 (5.27)	8.23 (7.46)	13.53 (5.3)	<0.001
Age group (years)	< 10	27203 (20.89%)	419 (54.20%)	27622 (21.09%)	<0.001
	10 - 14	31143 (23.91%)	104 (13.45%)	31247 (23.85%)
	≥ 15	71882 (55.20%)	250 (32.34%)	72132 (55.06%)
Intubated, n (%)	No	129813 (99.71%)	465 (60.63%)	130278 (99.48%)	<0.001
	Yes	378 (0.29%)	302 (39.37%)	680 (0.52%)
Pneumonia, n (%)	No	127137 (97.63%)	260 (33.64%)	127397 (97.25%)	<0.001
	Yes	3091 (2.37%)	513 (66.36%)	3604 (2.75%)
Indigenous, n (%)	No	121488 (99.32%)	740 (97.75%)	122228 (99.31%)	<0.001
	Yes	835 (0.68%)	17 (2.25%)	852 (0.69%)
Diabetes, n (%)	No	129294 (99.40%)	717 (93.36%)	130011 (99.36%)	<0.001
	Yes	780 (0.60%)	51 (6.64%)	831 (0.64%)
Obesity, n (%)	No	124933 (96.04%)	698 (90.53%)	125631 (96.01%)	<0001
	Yes	5154 (3.96%)	73 (9.47%)	5227 (3.99%)
COPD, n (%)	No	129970 (99.91%)	765 (99.48%)	130735 (99.90%)	<0.001
	Yes	121 (0.09%)	4 (0.52%)	125 (0.10%)
Asthma, n (%)	No	126275 (97.07%)	757 (98.44%)	127032 (97.08%)	0.025
	Yes	3806 (2.93%)	12 (1.56%)	3818 (2.92%)
Immunosuppression	No	129258 (99.37%)	687 (89.69%)	129945 (99.31%)	<0.001
	Yes	820 (0.63%)	79 (10.31%)	899 (0.69%)
Hypertension, n (%)	No	129409 (99.48%)	717 (93.24%)	130126 (99.45%)	<0.001
	Yes	674 (0.52%)	52 (6.76%)	726 (0.55%)
Cardiovascular disease, n (%)	No	129471 (99.53%)	730 (94.93%)	130201 (99.50%)	<0.001
	Yes	615 (0.47%)	39 (5.07%)	654 (0.50%)
CKD, n (%)	No	129753 (99.75%)	713 (92.72%)	130466 (99.70%)	<0.001
	Yes	331 (0.25%)	56 (7.28%)	387 (0.30%)
Patient type, n (%)	Outpatients	124094 (95.29%)	50 (6.47%)	124144 (94.77%)	<0.001
	Inpatients	5439 (4.18%)	513 (66.36%)	5952 (4.54%)
	ICU	695 (0.53%)	210 (27.17%)	905 (0.69%)
Comorbidities, n (%)	0	117708 (90.39%)	400 (51.75%)	118108 (90.16%)	<0.001
	1	10860 (8.34%)	246 (31.82%)	11106 (8.48%)
	Mayor de 1	1660 (1.27%)	127 (16.43%)	1787 (1.36%)

COPD, chronic obstructive pulmonary disease; CKD, chronic kidney disease; ICU, intensive care unit.

Table 2 presents the environmental characteristics where children and adolescents resided at the time of diagnosis of SARS-CoV-2 infection. A total of 198 (25.61%) deaths occurred in children and adolescents living in areas with high and very high social lag index. COVID-19 mortality was significantly different among the indicators related to social lag index, except for population not affiliated to any health service. Deceased children and adolescents lived in areas with a greater number of people with malnutrition compared to children and adolescents who survived (an average of 2022.68 vs 1872.03 inhabitants). Finally, children and adolescents living in areas with higher rates of complete vaccination coverage who were one year old and up to two years old showed lower mortality percentage than those with “schedule with 4 vaccines”.

Table 2

Neighborhood characteristics for the entities where children and adolescents resided according to COVID-19 outcome.

		Alive (n = 130,228)	Dead (n = 773)	Total (n = 131,001)	p-value*
Social lag index, n (%)	Very low	65564 (50.35%)	148 (19.15%)	65712 (50.16%)	<0.001
	Low	39436 (30.28%)	375 (48.51%)	39811 (30.39%)
	Medium	10842 (8.33%)	52 (6.73%)	10894 (8.32%)
	High	8943 (6.87%)	95 (12.29%)	9038 (6.90%)
	Very high	5443 (4.18%)	103 (13.32%)	5546 (4.23%)
Indicators related to social lag, mean (SD)	Percentage of population aged 15 and over that is illiterate	3.05 (2.41)	4.42 (3.21)	3.05 (2.42)	<0.001
	Percentage of population with children of school age of 6–14 years who do not attend school	5.48 (0.77)	5.84 (0.97)	5.48 (0.78)	<0.001
	Percentage of population with incomplete basic education	23.76 (7.34)	28.74 (7.78)	23.79 (7.35)	<0.001
	Percentage of population not affiliated to any health service	25.78 (4.88)	25.69 (5.88)	25.78 (4.89)	0.626
	Percentage of inhabited private housing units without floor covering	2.04 (2.47)	3.4 (3.56)	2.05 (2.48)	<0.001
	Percentage of inhabited private housing units without toilets or sanitary services	1.29 (1.45)	1.96 (1.88)	1.29 (1.46)	<0.001
	Percentage of inhabited private housing units without piped water from the public network	2.31 (2.22)	3.29 (2.94)	2.32 (2.23)	<0.001
	Percentage of inhabited private housing units without drainage	2.35 (3.31)	4.18 (4.26)	2.36 (3.32)	<0.001
	Percentage of inhabited private housing units without electricity services	0.44 (0.49)	0.72 (0.56)	0.44 (0.49)	<0.001
	Percentage of inhabited private housing units without washing machine	21.83 (8.47)	26.56 (12.16)	21.86 (8.5)	<0.001
	Percentage of inhabited private housing units without refrigerator	9.06 (5.51)	12.12 (7.74)	9.08 (5.53)	<0.001
Indicators related to malnutrition, mean (SD)	Number of people with malnutrition	1872.03 (1105.02)	2022.68 (1431.78)	1872.92 (1107.29)	<0.001
Vaccination coverage for children of one year of age, mean (SD)	Complete schedule	21.22 (7.04)	17.69 (8.46)	21.2 (7.06)	<0.001
	Schedule with 4 vaccines	22.01 (6.93)	18.9 (8.65)	21.99 (6.94)	<0.001
Vaccination coverage for children up to age two, mean (SD)	Complete schedule	34.75 (6.46)	32.84 (9.23)	34.74 (6.48)	<0.001
	Schedule with 4 vaccines	37.61 (6.51)	35.08 (9.44)	37.6 (6.54)	<0.001
Population density, mean (SD)	2782.28 (2960.73)	674.12 (1585.66)	2769.84 (2958.9)	<0.001

Significant value (p < 0.05).

Table 3 shows results of the multilevel model to estimate association between COVID-19 mortality and individual characteristics, adjusted for population density. From this model, it was estimated that the median odds of mortality due to COVID-19 in an area with a high mortality would be approximately twice that of one with a low mortality (MOR = 1.797, 95% CI: 1.482–2.113). This analysis demonstrated a negative association between COVID-19 mortality and age (OR = 0.878, 95% CI: 0.869–0.888). Children and adolescents with diabetes (OR = 3.898, 95% CI: 2.596–5.851), immunosuppression (OR = 5.410, 95% CI: 4.088–7.158), obesity (OR = 1.876, 95% CI: 1.397–2.521), hypertension (OR = 1.906, 95% CI: 1.239–2.932), cardiovascular disease (OR = 2.288, 95% CI: 1.482–3.531), and chronic kidney disease (OR = 13.250, 95% CI: 9.066–19.350) showed a higher risk of mortality.

Table 3

Factors associated with COVID-19 mortality in children and adolescents in Mexico. Results from the multilevel logistic regression models adjusted for sex, age and comorbidities.

	OR	IC 95%	P-value*
Sex (Ref. Female)	1.115	0.962 - 1.292	0.147
Age	0.878	0.869 - 0.888	0.000
Diabetes	3.898	2.596 - 5.851	0.000
Obesity	1.876	1.397 - 2.521	0.000
Immunosuppression	5.410	4.088 - 7.158	0.000
Hypertension	1.906	1.239 - 2.932	0.003
Cardiovascular disease	2.288	1.482 - 3.531	0.000
CKD	13.250	9.066 - 19.350	0.000
Population density	0.374	0.204 - 0.688	0.002
MOR	1.797	1.482 - 2.113	0.000

Significant value (p <  0.05).

CKD, chronic kidney disease.

Regarding environmental factors, social lag index and malnutrition were positively associated with COVID-19 mortality (ORvery high = 2.939, 95%CI: 1.111-7.775, and OR = 1.390, 95%CI: 1.073-1.802, respectively). Population density presented a negative association with COVID-19 mortality (OR = 0.374, 95%CI: 0.204–0.688). The models were adjusted by age, sex, diabetes, obesity, immunosuppression, hypertension, and cardiovascular disease. In addition, the malnutrition model was also adjusted for population density (Figure 2 ).

Figure 2

Adjusted associations between environmental factors and COVID-19 mortality.

Odds Ratio (OR); 95% Confidence interval (CI). Values were obtained using multivariable multilevel logistic regression models adjusted for sex, age, diabetes, obesity, immunosuppression, hypertension, cardiovascular disease and CKD. The models with the malnutrition variables were also adjusted for population density.

Figure 3 presents the odds of the indicators that make up the social lag index and COVID-19 mortality adjusted by state and individual characteristics. A higher percentage of population aged 15 and over with incomplete education was associated with a higher risk of COVID-19 mortality (OR = 1.045, 95%CI: 1.011–1.081), as well as a higher percentage of population with children of school age of 6-14 years who do not attend school (OR = 1.266, 95%CI 1.032–1.554), and a higher percentage of population aged 15 and over that are illiterate (OR = 1.086, 95%CI 0.999-1.179).

Figure 3

Adjusted associations between indicators relates to social lag index and COVID-19 mortality.

Odds Ratio (OR); 95% Confidence interval (CI). Values were obtained using multivariable multilevel logistic regression models adjusted for sex, age, diabetes, obesity, immunosuppression, hypertension, cardiovascular disease and CKD.

Discussion

In this study, we assessed the influence of environmental and individual factors on COVID-19 mortality in children and adolescents in Mexico during the first two years of the pandemic. The most relevant findings of this work can be summarized as follows: (i) Children and adolescents living in areas with higher prevalence of malnutrition and very high social lag indexes (especially its indicators on access to education) showed greater odds of dying from COVID-19; and (ii) Regarding individual characteristics, age and pre-existing comorbidities such as diabetes, obesity, immunosuppression, hypertension, cardiovascular disease, and CKD were associated with an increased risk of mortality.

Child malnutrition is a serious public health issue in Mexico. , Malnutrition was associated with socially disadvantaged populations and the pandemic could have worsened this problem. Malnutrition is actually a risk factor associated with mortality in some viral infection such as influenza A (H1N1). The impact of undernutrition in COVID-19 outcomes became a subject of interest. Our findings highlight the relevance of malnutrition in COVID-19 mortality among children and adolescents in Mexico. Mertens and Peñalvo described a potential relationship between COVID-19 mortality and areas with an elevated burden of undernutrition. Although this analysis used country-level data, its conclusions were robust enough to consider the importance of malnutrition in the prognosis of the COVID-19 disease. A study conducted in India concluded that malnutrition among under-five children was a risk factor of COVID-19. Alec Kurtz et al. showed that malnutrition was correlated with COVID-19 severity in an age-dependent way, and found the worst odds for severe COVID-19 in children between ages 6 and 17 with a history of malnutrition.

Social inequalities in COVID-19 mortality were previously studied. In the present cohort, a very high category in social lag index was associated with a worse COVID-19 disease prognosis in children. Previous research found similar results in other age groups. Other studies from Brazil and Mexico showed that the poorest population groups have higher mortality from COVID-19. , Lower education level showed a strong association with COVID-19 mortality in recent studies in the United States , Ribeiro et al. showed similar findings in Brazil: Lower education level was the strongest association with COVID-19 fatalities.

Additionally, our findings showed that living in an area with a high population density is associated with lower odds of dying from COVID-19. Some studies conducted by researchers at the Johns Hopkins Bloomberg School of Public Health found that the death rate is inversely related to population density, since dense counties had a significantly lower rate of COVID-19 mortality. , This is in contrast to others assumptions, but density could lead to better access to health care facilities and infrastuctures, as well as achieving a greater impact with other measures like physical distancing and lockdown.

Our findings also pointed to a connection between obesity and a worse COVID-19 prognosis. Several studies found a strong association between this variable and short-term mortality caused by COVID-19 in the general population.42, 43, 44 Tripathi et al. focused on the impact of previous conditions on COVID-19 outcomes in children. This work, based on a large cohort in the United States, found that obesity was a key factor in understanding the evolution of COVID-19 in this group.

Other comorbidities were associated with COVID-19 mortality in children and adolescents in our analysis. According to previous studies, other underlying comorbidities may lead to a worse outcome in children and adolescents. , , A recent meta-analysis showed that excessive adipose tissue, deficit in lean mass, insulin resistance, dyslipidemia, hypertension, high levels of proinflammatory cytokines, and low intake of essential nutrients are associated with damage to immune, cardiovascular, respiratory, and urinary systems, as well as with modification of the intestinal microbiota (dysbiosis). Moreira et al. found that children with pre-existing conditions, such as asthma, autoimmune disease, cardiovascular disease, chronic lung disease, gastrointestinal/liver disease or hypertension, had an increased risk of hospitalization and death. Differences between baseline characteristics in these cohorts and access to healthcare systems should be considered when interpreting these findings.

Regarding age, our analysis showed that younger patients showed a higher risk of COVID-19 death. This finding is controversial considering previous studies. A large cohort in the United States found that younger age increases the risk of hospitalization, but not death. However, the review by Kitano et al. found that over 90% of pediatric COVID-19 deaths were from Low-Middle-Income Countries and that mortality risk was highest in children under 1-year old. Potential explanations to the higher risk of COVID-19 mortality Mexican younger patients could involve individual factors as immune function differences and lack of prior exposure to other coronaviruses. Also, the important role of social determinants of health and the quality of health care system could be related to the negative results.

This study has several strengths. Data from a relevant number of Mexican patients were analyzed. To our knowledge, this is the largest pediatric evaluation assessing factors associated with COVID-19 mortality in children and adolescents in Mexico. Nonetheless, some limitations must be considered. First, the presence of comorbidities was self-reported (or reported by a caregiver), which could lead to misclassification bias. In addition, there is a small group of patients who reported other comorbidities that are nonspecific and unrelated to an increased risk of severe infection and poor prognosis. Second, we used cross-sectional information, limiting our ability to establish a clear causal link. Third, the use of the sentinel surveillance system to detect and report COVID-19 cases in Mexico likely skews detection towards more severe cases. Also, since we rely on data from the population who had symptoms and sought health care, we do not have any information on patients who were asymptomatic or mildly symptomatic and did not seek health care. Another limitation of our study is that we cannot exclude the possibility that the number of deaths in patients who had COVID-19 could be underreported. Finally, the environmental factors were measured at the state level, since our ability to geocode the patients was restricted to that level. States are fairly large and heterogeneous, composed by a variable mix of rich and poor neighborhoods. Nevertheless, despite its limitations, we used a large national registry of COVID-19 cases manteined by the Mexican government and many of the cases were tested using both Rapid Ag-T and RT-PCR.

Our findings show the impact of malnutrition, poverty and pre-existing comorbidities in the prognosis of COVID-19 pandemic. The relevance of malnutrition seems to be associated with population density. This finding could be related to the scarce access to health care services in rural areas. Our study can have relevant implications for reaching a more comprehensive understanding of the factors associated with COVID-19 mortality in children and adolescents.

Contributors

CSP, and FJPG conceptualised the analysis. CSP, and FJPG accessed and verified the data, completed the formal analysis, and drafted the original manuscript. All authors had access to all data and contributed to study design, data collection, and manuscript editing. All authors had full access to all the data in the study and had final responsibility for the decision to submit for publication.

Data sharing

The datasets generated and/or analyzed during the current study are available in: https://datos.gob.mx/busca/dataset/informacion-referente-a-casos-covid-19-en-mexico/resource/f9418a4c-10e6-4099-8036-f51234b0e81b

Declaration of interests

The authors declare no competing interests.