Article Text

other Versions

Original research
Seroprevalence of the SARS-CoV-2 antibody in healthcare workers: a multicentre cross-sectional study in 10 Colombian cities
  1. Jeadran Nevardo Malagón-Rojas1,2,
  2. Marcela Mercado-Reyes3,
  3. Yezith G Toloza-Pérez2,
  4. Eliana L Parra Barrera2,
  5. Marien Palma2,
  6. Esperanza Muñoz4,
  7. Ronald López2,
  8. Julia Almentero2,
  9. Vivian V Rubio3,
  10. Edgar Ibáñez5,
  11. Eliana Téllez2,6,
  12. Lucy G Delgado-Murcia7,
  13. Claudia P Jimenez8,9,
  14. Diego Viasus-Pérez10,
  15. Marisol Galindo3,
  16. Luisa Lagos2
  17. Interinstitutional Support Group
    1. 1Doctorado en Salud Pública, El Bosque University Faculty of Medicine, Bogota, Colombia
    2. 2Grupo de investigación en Salud Ambiental y Laboral, Instituto Nacional de Salud, Bogotá, Colombia
    3. 3Research Directorate, National Institute of Health, Bogota, Colombia
    4. 4Faculty of Nursing, National University of Colombia, Bogota, Colombia
    5. 5El Bosque University Faculty of Medicine, Bogota, Colombia
    6. 6Red de Investigación, Innovación y Desarrollo en Seguridad y Salud en el Trabajo, Bogota, Colombia
    7. 7Immunotoxicology Research Group, National University of Colombia—Bogota Headquarters, Bogota, Colombia
    8. 8National Open and Distance University, Bucaramanga, Colombia
    9. 9Red de Investigación, Innovación y Desarrollo en Seguridad y Salud en el Trabajo, Bogotá, Colombia
    10. 10School of Medicine, Universidad del Norte, Barranquilla, Colombia
    1. Correspondence to Dr Jeadran Nevardo Malagón-Rojas, Grupo de investigación en Salud Ambiental y Laboral, Instituto Nacional de Salud, Bogota, Colombia; jnmalagon{at}unbosque.edu.co

    Abstract

    Background Healthcare workers are at increased risk of infection due to occupational exposure to SARS-CoV-2-infected patients. The objective of this study was to determine the seroprevalence of SARS-CoV-2 in healthcare workers in Colombia.

    Methods This study is a cross-sectional study focused on estimating the seroprevalence of SARS-CoV-2 antibodies in healthcare workers from 65 hospitals in 10 cities in Colombia during the second semester of 2020. The seroprevalence was determined using an automated immunoassay (Abbott SARS-CoV-2 CLIA IgG). The study included a survey to establish the sociodemographic variables and the risk of infection. A multivariate model was used to evaluate the association between the results of seroprevalence and risk factors.

    Results The global seroprevalence of antibodies against SARS-CoV-2 was 35% (95% Bayesian CI 33% to 37%). All the personnel reported the use of protective equipment. General services personnel and nurses presented the highest ratios of seroprevalence among the healthcare workers. Low socioeconomic strata have shown a strong association with seropositivity.

    Conclusion This study estimates the prevalence of SARS-CoV-2 infection among healthcare workers. Even though all the personnel reported the use of protective equipment, the seroprevalence in the general services personnel and nurses was high. Also, a significant difference by cities was observed.

    • COVID-19
    • public health
    • occupational health
    • health personnel
    • cross-sectional studies

    Data availability statement

    Data are available upon reasonable request. Data may be obtained from a third party and are not publicly available.

    http://creativecommons.org/licenses/by-nc/4.0/

    This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.

    Statistics from Altmetric.com

    Request Permissions

    If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

    Key messages

    What is already known about this subject?

    • Worldwide, healthcare workers (HCWs) have been mainly affected due to the direct contact with patients with COVID-19.

    • Serological studies are essential to understand the infection in healthcare and to establish control strategies.

    • Studies have estimated that the seroprevalence of antibodies against SARS-COV-2 in health workers is highly variable and is correlated with the proportion of infections at the local level.

    • Latin America has one of the highest rates of COVID-19 infection among HCWs; except for some reports from specific hospitals, there is a lack of studies of seroprevalence on this group of workers in the region.

    What are the new findings?

    • The present study was carried out during the second semester of 2020.

    • We reported a seroprevalence rate of SARS-CoV-2 of 35% in HCWs, one of the higher ratios of HCWs reported worldwide.

    • The SARS-CoV-2 seroprevalence in HCWs varied in the Colombian cities from 21% to 71%.

    • Workers from general services and nurses showed greater seroprevalence ratios.

    How might this impact on policy or clinical practice in the foreseeable future?

    • Guidance for prevention of the transmission of SARS-CoV-2 may include extra labour factors, such as recommendations for the protection of close relatives and the promotion of protective behaviours out of working hours.

    Introduction

    SARS-CoV-2 has generated multiple and diverse challenges worldwide in all areas of work. One of the work environments that has attracted the most attention is the provision of health services, given the relevance of healthcare workers (HCWs) and their role during the pandemic. HCWs possibly have the highest exposure and risk of infection because they are in direct contact with infected patients. The WHO estimates that around 14% of reported COVID-19 cases correspond to HCWs even reaching 35% in some countries.1 By September 2020, in the Americas region, almost 570 000 HCWs were reported with COVID-19, in addition to 2500 deaths from SARS-CoV-2.2 Around the world, it is estimated that the COVID-19 infections among this specific group range between 1% and 45.3%, being higher in countries from the northern hemisphere.3–10 It was indicated that HCWs who are male and from ethnic minorities resulted in higher seroprevalence level.11

    By 1 July 2021, nearly 4.4 million cases and more than 110 000 COVID-19-related deaths have been reported in Colombia.12 The country has experienced three peaks: the first between July and August 2020, the second between December 2020 and January 2021, and the last one between April and June 2021. The average effective reproductive number (R0) for the country during 2020 was 1.09, while the estimated average R0 between January and July 2021 was 0.92.12 During the lockdown period, the Colombian health system increased its capability of intensive care unit (ICU) beds number and boosted HCWs’ response, as well as the number of labs able to carry out RT-PCR in the country.

    Since the pandemic declared by the WHO began in March 2020,13 various active surveillance strategies have been implemented, such as self-reporting of symptoms through mobile applications, the usage of RT-PCR tests regularly and performing serological tests to the identification of antibodies against SARS-CoV-2.14 Due to the costs of surveillance strategies through the identification of nucleic acids, this type of intervention has been poorly adopted among health institutions. Given that up to 50% of SARS-CoV-2 cases correspond to asymptomatic infections,15 the incidence of SARS-CoV-2 cases based on the notification of confirmed SARS-CoV-2 cases implies a considerable underestimation of the incidence of this virus infection.16 In this sense, the carrying out of studies allowed us to estimate the proportion of HCWs who have antibodies against SARS-CoV-2. Besides, this type of analysis also provides information about the immune response to the virus, natural susceptibility, as well as useful information when prioritising the application of the vaccine.

    HCWs are the most exposed to the risk of infection with the new SARS-CoV-2. This study aimed to determine the seroprevalence of antibodies against SARS-CoV-2 among HCWs in Colombia, along with describing the associations between seroprevalence and occupational exposure to SARS-CoV-2 in 10 Colombian cities.

    Materials and methods

    This was a cross-sectional study of HCWs across medical services in 10 cities of Colombia from September to November 2020. The study was designed following the recommendations from the Strengthening the Reporting of Observational Studies in Epidemiology statement for observational studies.17

    Sample

    A cross-sectional study with non-probability sampling was designed for health workers from public and private hospitals in 10 cities in Colombia. To calculate the sample size, an expected seroprevalence of 30% (p=0.30), (q=1–0.30=0.70) was taken with a marginal sampling error of ±2% (δ=0.02) and a confidence level=95% (α=0.05, Zα=1.96). A 10% loss percentage was established to calculate a minimum sample size of 2241 participants. We used the definition of HCW designated by the WHO providers of healthcare attention.18 The public and private hospitals (IPS) were chosen using the municipal records, choosing the clinics and hospitals that concentrated 80% of the attention of COVID-19 cases in the municipality. Health workers were invited to participate through the personnel office. Participants in the study were selected from a list of volunteers in each IPS. The selection was performed using a random number generator in Excel. We included personnel either directly or indirectly involved in the healthcare attention: doctors, nurses, pharmacists, physiotherapists, respiratory therapists, bacteriologists, healthcare technicians, admission assistants, general services (catering and cleaning staff), and security personnel.18

    Serum samples and serological test

    Serum samples obtained from 6 to 7 mL of venous blood were collected. Samples were refrigerated and transported to a local laboratory in each city. Later, samples were centrifuged to separate the serum and were stored at −30°C to −80°C until processing at the reference laboratory at the Instituto Nacional de Salud. The detection of total antibodies was made by the chemiluminescence technique (CLIA) ‘SARS-CoV-2 Total Advia Centaur-Siemens’. The Advia Centaur-Siemen’s test detects serum total antibodies against SARS-CoV-2. According to the manufacturer, the range of index values oscillates between 0.05 and 10 (cut-off point of reactive=>10).

    The CLIA test was selected after performing a secondary validation with samples from the Colombian population. The sensitivity and specificity of the test were 86% (95% CI 79% to 91%) and 99% (95% CI 96% to 100%), respectively.19

    Collection of samples and field surveys

    An electronic questionnaire was applied online using Google Forms. The questions were based on the guidelines from the WHO.18 The questionnaire included questions for sociodemographic characterisation, usage of personal protective equipment, characterisation of work conditions and dwelling, and previous exposure to COVID-19. Also, participants who declared having had COVID-19 were asked whether their insurance recognised their COVID-19 episode as being work related.

    Five experts were asked to perform a virtual judgement to validate the content of the instrument.20 The criteria included clarity, coherence, relevance and sufficiency. The criteria were evaluated on a scale of 1–5 by each expert. Finally, the questionnaire was validated by 300 HCWs from Bogota. A Spanish version of the questionnaire is available (online supplemental material 1).

    Supplemental material

    Statistical analysis

    Sociodemographic characteristics of HCWs were described for each city. For the quantitative variables, means and SD were estimated. Subsequently, a bivariate analysis was performed comparing the nominal or ordinal variables regarding the presence or absence of antibodies against SARS-CoV-2, analysed using Pearson’s Χ2 test with Yates correction. In the case of quantitative variables, the Spearman correlation was used. The level of statistical significance established was p<0.05.

    A Poisson regression model was applied to test the relationship between the results of the CLIA tests (dependent variable) and the independent variables. Independent variables were selected from the bivariate results. Variables with a p<0.2 were included in the multivariate model. Associations were presented in prevalence ratios with 95% CIs. The Akaike information criterion was used as an estimator of the quality of statistical models.21

    The statistical analysis was conducted using R (V.4.0.3). The overall crude frequencies of seropositivity tests were estimated. Later, the crude seroprevalence was stratified by age, sex, ethnicity and role at the IPS. For both cases, crude seroprevalence was adjusted using the Bayesian method in R V.2.21.2 (pack RStan)22 using the data of sensitivity and specificity reported in previous studies carried out using CLIA in Colombian populations.19 Also, 95% Bayesian credibility intervals (BCIs) were obtained.

    The model was run according to the following equation:

    Real prevalence=adjusted prevalence×sensitivity+(1−adjusted prevalence)×(1−specificity) (see online supplemental file 2).

    Supplemental material

    Results

    Sociodemographic description

    A total of 4042 HCWs were included in the study in 65 health centres from 10 cities. The questionnaire was completed by 82% (n=3294) of the participants. The mean age was 36.45±10.5 years old. Most of the participants belonged to socioeconomic strata 2 and 3 (63.4%). Most of the participants were nurses (34.9%) and general practitioners (10.1%). The distribution of participants by cities were Bogotá (n=677; 16.7%), Bucaramanga (n=508; 12.6%), Cali (n=500; 12.4%), Medellín (n=470; 11.6%), Barranquilla (n=434; 10.7%), Cucuta (n=423; 10.5%), Villavicencio (n=395; 9.8%), Ipiales (n=388; 9.6%), Leticia (n=176; 4.4%) and Guapi (n=71; 1.8%) (figure 1).

    Figure 1

    Distribution of seroprevalence by city. The figure presents in greyscale the prevalence of antibodies against SARS-CoV-2 in each city included in the study. Cities are marked with a white dot. The greater the intensity of the grey, the greater the seroprevalence was found in the city.

    SARS-CoV-2 prevalence

    The prevalence of self-declared SARS-CoV-2 infection was 30.06% (95% CI 29.01% to 32.22%; n=995). The percentage of workers who declared having been diagnosed with the PCR test was 89.54% (95% CI 87.49% to 91.3%; n=891). The proportion of seropositivity among the workers who declared the COVID-19 infection was 81% (95% CI 77.50% to 84.0%). The proportion of workers who received legal recognition of COVID-19 infection as a work-related disease was 40.60% (95% CI 37.50% to 43.70%; n=404).

    Protective equipment

    The adherence to using protective equipment such as face masks (disposable surgical and N95 masks) was 100% among HCWs. It was observed that male participants tend to perform a lower number of protective elements (Χ2=44.69, p=0.00001). Participants who declared having been vaccinated against influenza during the last year have a lower proportion of seroconversion than those who did not receive the vaccine (Χ2=9.7425, p=0.0001).

    Seroprevalence of SARS-CoV-2 in HCWs

    The global seroprevalence was 35% (95% BCI 33.0% to 37.0%). The highest seroprevalence by cities was found in Guapi (71%), Villavicencio (54%) and Barranquilla (44%) (table 1). No significant difference was found between male and female HCWs concerning seropositivity. In the bivariate analysis, there was an association between seroprevalence and occupation, age, socioeconomic strata and educational level (p<0.05) (table 1). There was an association between seropositivity and families with two or more members (Χ2=7.74; p=0.005).

    Table 1

    Demographics, symptomatology, occupational exposure and seroprevalence in healthcare workers in Colombia

    The seropositivity was higher among personnel from general services (48%; 95% CI 37% to 59%) and nurses (46%; 95% CI 42% to 49%) (figure 2). The occupation with the lowest seropositivity rate was physiotherapist (7%; 95% CI 0% to 18%). A reversed social gradient was found out between the presence of antibodies against SARS-CoV-2 and socioeconomic level (Χ2=100.87; p trend=0.0000001) (figure 3).

    Figure 2

    Seroprevalence among healthcare workers in Colombia by occupation. *Other: nutritionists, Rx technician, physiotherapist and clinical laboratory technician. The figure shows the crude and adjusted seroprevalence levels among different healthcare occupations. Seroprevalence levels were higher among nurses and general services workers.

    Figure 3

    Seroprevalence and socioeconomic stratum in Colombia. The figure presents the differences in the seroprevalence of antibodies against SARS-CoV-2 among healthcare workers from different socioeconomic levels. There is a linear trend showing higher seroprevalence levels in lower socioeconomic strata (1–3) and lower seroprevalence in higher socioeconomic strata (4–6).

    The multivariate model showed that participants from lower socioeconomic strata have more chance of having a reactive CLIA test (table 2). Besides, HCWs with blood type AB compared with type O were 68% more likely to have a reactive CLIA test. People who work in the emergency room and hospitalisation were more likely to have a reactive CLIA test (by 57% and 37%, respectively). Participants who worked in the ICU and COVID-19 services had not increased their risk of having a reactive test and were not significant in the multivariate model (table 2).

    Table 2

    Multivariate analysis of seroprevalence by demographic characteristics

    Discussion

    HCWs are a population with a high-risk of acquiring SARS-CoV-2 infection due to direct contact with patients.7 23 24 We conducted a study to assess the seroprevalence of SARS-CoV-2 infection associated with characteristic demographics and the occupation of HCWs from 65 hospitals and medical centres in 10 Colombian cities. We observed an overall seroprevalence of 32%. Except for Bucaramanga (26%), the seroprevalence was higher among workers in cities with less than 1.5 million inhabitants: Guapi (71%), Villavicencio (54%), Leticia (43%) and Ipiales (37%). The seroprevalence in cities with larger populations was lower: Bogotá (34%), Cali (35%), Cúcuta (27%) and Medellín (22%).

    Comparing our findings with the reported seroprevalence of SARS-CoV-2 antibodies in the general population in Colombia (September–December 2020),25 except for Bucaramanga (32%), the seroprevalence in the HCWs tended to be lower in cities located in the North region (Barranquilla 55%, Cúcuta 40% and Medellín 27%). Nevertheless, the seroprevalence of SARS-CoV-2 antibodies in HCWs in cities from the Central, Southern and Western region of the country was higher than that reported in the general population (Bogota 30%, Leticia 59%, Villavicencio 34%, Guapi 68% and Ipiales 35%).25

    The seroprevalence in HCWs from Bogota was higher than in a previous study in one hospital in the city carried out in August 2020 (8.26%)26 but similar to the seroprevalence reported in a cohort of airport workers in Bogotá (September 2020).27

    The reported seroprevalence of antibodies against SARS-CoV-2 in HCWs was greater than that reported in the studies conducted during the second semester of 2020 in North America (12.7%), Africa (8.2%) and Asia (4%).28 Besides, reports from studies carried out in European countries are lower than our results: Denmark (4.04%), England (24.4%), Germany (4.36%),29 Greece (1.26%),30 Italy (14.4%) and Switzerland (1%).31

    We have found no differences in the distribution of seroprevalence between men and women. This issue has been approached in several seroprevalence studies. A recent meta-analysis has reported that seroprevalence levels were higher among male HCWs.11 Another systematic review observed higher seroprevalence ratios among men.28 This association may be correlated with the fact that men tend to show less adherence to protective protocols compared with women.32 In our study, we observed that men seem to be less willing to use all personal protection items compared with women.

    Also, we observed that occupations that are performed mostly by women were associated with a higher risk of infection. It has been stated that gender is a social determinant of health, linked to the health disparities among the COVID-19 pandemic.33 Also, it was highlighted that personal protective equipment does not protect female HCWs as well as their male colleagues. It has been pointed out, for example, that the glasses do not fit their faces, the gloves are too long and the face shields collide with the chest, making it uncomfortable to perform procedures.34 35 These conditions constitute a relevant concern considering that the COVID-19 pandemic has highlighted the extent to which society depends on women, both in the first line of response in the health sector, as well as at homes. Women constitute the majority of the workforce in the health sector36 and in Colombia, more than 70% of HCWs are women.37 Nevertheless, these statistics do not include personnel involved in activities of cleaning and catering. Women have an increased risk of contracting SARS-CoV-2 given the close interaction with patients and visitors amid shortages of personal protective equipment.33 Also, women are concentrated in roles requiring the closest, prolonged contact with patients.38

    We did not observe differences in seroprevalence levels regarding ethnicity. Nevertheless, several authors have reported that exposition levels are higher among Afro-Americans and Hispanics.11 28 39

    Also, an association between the presence of antibodies in serum and the social stratum in which the worker resides was found. It was observed that there is a reverse gradient in the seroprevalence proportion; as far as the stratum increases, the proportion of people with a reactive test for antibodies decreases. These findings have been previously described at the community level, showing the association between socioeconomic aspects and COVID-19 transmission,40 COVID-19 severity41 and antibody presence.42 Nevertheless, to our knowledge, this is the first report showing the links between seroprevalence and socioeconomic strata among the HCWs. A study carried out among a group of workers with high risk of SARS-CoV-2 transmission has stated that prevention programmes should include extra labour factors such as including recommendations for biological protection at home, supermarket and other places.43 In this sense, various authors have indicated that most of the COVID-19 cases took place in places such as homes.44 45 Considering that only 40% of the self-declared COVID-19 cases were recognised as linked to occupational activities, extra occupational risk factors such as positive close contact at home, family size and house conditions should be studied in detail to understand the SARS-CoV-2 transmission in HCWs.46

    Additionally, the regression model found out that participants who work in the emergency room and hospitalisation were more likely to have antibodies against SARS-CoV-2. Conversely, workers from the ICU and COVID-19 services had not increased their risk of having a reactive test. These findings differ from the literature reports which have stated that people working on COVID-19 units have an increased risk of having a positive SARS-CoV-2 antibody test.8 47 48 Nevertheless, our results have shown that participants who work in emergency services showed a significant increase in the risk of having an antibody reactive test. A&E characteristics vary depending on the location, clinical specialties and availability of technology. Most of the emergency services in Colombia have not divided the attention of respiratory cases from other emergencies. Also, patients could stay for longer periods waiting for diagnosis, treatment or transference. These work conditions may increase the exposure risk of HCWs in emergency areas.

    Studies have shown a slightly increased infection among non-O types. Also, the risk of intubation decreased in type A.49 Here, we observed an increased risk in the AB blood types. However, this result may contribute to the knowledge of blood type and the relationship with the role of the infection with COVID-19.

    The present study has limitations. First, the aspects related to the design of the research. A cross-sectional study was formulated and carried out between September and December 2020. Also, the characteristics of the sampling may introduce a selection bias among the workers who have been previously infected. Second, the HCWs answered the survey in different moments of the pandemic, even several months after having the COVID-19 episode, which may introduce a recall bias. Likewise, we do not control or implement so methodology to determine if the exposure was occupational or if the disease was caused outside the medical centres. In addition, the test used to identify antibodies against SARS-COV-2 did not allow quantifying the differences in antibody titres among the reactive participants. In this sense, we could not perform an analysis including the role that the job position played in the generation or not of antibodies against SARS-CoV-2. Since the study used the declared information related to COVID-19 antecedent, it was not possible to determine how many of those who tested positive at CLIA had also a previous molecular test positive for SARS-CoV-2 infection. Finally, the study did not evaluate the inadequate use of personal protective equipment that has been associated with an increased risk of COVID-19.50–52 Moreover, the source of contagion was not determined in the study, though it could be a main key issue to the protection and to ensure treatment and recovery of the health workers.

    In conclusion, to our knowledge, this is the first national study to quantify the level of seropositivity to SARS-CoV-2 in HCWs in the Andean region. The impact of the transmission on HCWs varies significantly from one city to another. Our findings have important implications for understanding the spread of SARS-CoV-2 and for planning control programmes in this population, as it could be the information of seroprevalence before the introduction of the SARS-CoV-2 vaccine.

    Abstract translation

    This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.

    Data availability statement

    Data are available upon reasonable request. Data may be obtained from a third party and are not publicly available.

    Ethics statements

    Patient consent for publication

    Ethics approval

    The study proposal and protocol were approved by the ethics committee of Instituto Nacional de Salud (CEMIN 012/2020) and the Research Ethics Committee from Universidad del Norte (223/10/2020). We have obtained written informed consent from each participant in the study.

    Acknowledgments

    The authors thank the hundreds of volunteers, community leaders, health institutions and the thousands of participants of the study. Special thanks to Secretarias de Salud Distritales and Secretarias de Salud Departamentales from the cities involved in the study. Special thanks to the public institutions, Armada Nacional de Colombia, Policía Nacional de Colombia, Fontur, SATENA, and civil organisations that facilitated the transfers of personnel in the field and the logistics of the study in each city. Special recognition for the support of the field epidemiologists’ team of the Field-Epidemiology-Training-Program (FETP) from Colombia. Also, a special thanks to healthcare workers of all medical centres in the country who have treated patients with COVID-19 since the beginning of the pandemic. Recognition and honours to the health personnel who lost their lives in this noble fight and to those who continue in this tough challenge. Authors want to thank Owen Harrison for his valuable comments in the final version of the manuscript. Finally, JNM-R would like to express his sincere gratitude to Carolina Garrote, Manuel Cruz and Victoria Malagón.

    References

    Supplementary materials

    Footnotes

    • Twitter @jeadranmalagon

    • Collaborators Interinstitutional Support Group: Magdalena Weisner, Gloria Puerto, Lyda Muñoz-Galindo, María Teresa Herrera, Jhonantan Reales, Edwin Cárdenas-Villamil, Jessica Ortíz, Ligia Ovideo and Liliana Serrano, Stephany Botero, Helena María Rodríguez-Perea, Norma, Celis-Cruz and Juan Felipe Bedoya.

    • Funding acquisition

    • Resources

    • Writing – review and editing

    • Contributors JNM-R, MM-R, ELPB, ET, CPJ, LGD-M, DV-P, LL and EM designed the study. JNM-R, MM-R, JA, VVR, MG, EI, YGT-P, RL, MP and ELPB conducted the analysis and developed the first draft of the manuscript, and edited the paper according to coauthors’ suggestions. DV-P, EM and LGD-M provided their technical input to improve the final manuscript. All authors contributed to the final draft of the manuscript. JNM-R and MM-R accepts full responsibility for the finished article, had access to any data, and controlled the decision to publish.

    • Funding The present research was financed by the Instituto Nacional de Salud.

    • Map disclaimer The depiction of boundaries on the map(s) in this article does not imply the expression of any opinion whatsoever on the part of BMJ (or any member of its group) concerning the legal status of any country, territory, jurisdiction or area or of its authorities. The map(s) are provided without any warranty of any kind, either express or implied.

    • Competing interests None declared.

    • Provenance and peer review Not commissioned; externally peer reviewed.

    • Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.