- Research article
- Open Access
- Open Peer Review
Impact of multimorbidity: acute morbidity, area of residency and use of health services across the life span in a region of south Europe
BMC Family Practice volume 15, Article number: 55 (2014)
Concurrent diseases, multiple pathologies and multimorbidity patterns are topics of increased interest as the world’s population ages. To explore the impact of multimorbidity on affected patients and the consequences for health services, we designed a study to describe multimorbidity by sex and life-stage in a large population sample and to assess the association with acute morbidity, area of residency and use of health services.
A cross-sectional study was conducted in Catalonia (Spain). Participants were 1,749,710 patients aged 19+ years (251 primary care teams). Primary outcome: Multimorbidity (≥2 chronic diseases). Secondary outcome: Number of new events of each acute disease. Other variables: number of acute diseases per patient, sex, age group (19–24, 25–44, 45–64, 65–79, and 80+ years), urban/rural residence, and number of visits during 2010.
Multimorbidity was present in 46.8% (95% CI, 46.7%-46.8%) of the sample, and increased as age increased, being higher in women and in rural areas. The most prevalent pair of chronic diseases was hypertension and lipid disorders in patients older than 45 years. Infections (mainly upper respiratory infection) were the most common acute diagnoses. In women, the highest significant RR of multimorbidity vs. non-multimorbidity was found for teeth/gum disease (aged 19–24) and acute upper respiratory infection. In men, this RR was only positive and significant for teeth/gum disease (aged 65–79). The adjusted analysis showed a strongly positive association with multimorbidity for the oldest women (80+ years) with acute diseases and women aged 65–79 with 3 or more acute diseases, compared to patients with no acute diseases (OR ranged from 1.16 to 1.99, p < 0.001). Living in a rural area was significantly associated with lower probability of having multimorbidity. The odds of multimorbidity increased sharply as the number of visits increased, reaching the highest probability in those aged 65–79 years.
Multimorbidity is related to greater use of health care services and higher incidence of acute diseases, increasing the burden on primary care services. The differences related to sex and life-stage observed for multimorbidity and acute diseases suggest that further research on multimorbidity should be stratified according to these factors.
Concurrent diseases, multiple pathologies and multimorbidity patterns are topics of increased interest as the world’s population ages . Multimorbidity is the coexistence of two or more chronic health problems in the same person at one point in time , and multimorbidity patterns are any combination of chronic diseases . Both considerations have important consequences for the individual and for health services . Multimorbidity is a challenge for industrialized countries and can jeopardize the viability of national health systems.
Traditionally, the construct of multimorbidity has been inherently associated with persistent or chronic disease. Methods to measure multimorbidity include disease scores, case-mix systems, indexes and disease counts, the latter being the common method . Far too little attention has been paid to the use of health services and the role of urban or rural residency in patients with multimorbidity [6, 7]. Furthermore, the literature lacks comparisons by sex of acute morbidity in patients with multimorbidity in a large population sample [8, 9].
The classification of acute and chronic disease remains controversial. Acute disease is characterized by a single or repeated episode of relatively rapid onset and short duration with a recovery to previous stage of activity . Nevertheless, some diseases fall into a grey area. Knowledge of specific acute diseases that may occur more frequently than expected and of the underlying vulnerabilities  could help to focus attention on the patients with multimorbidity rather than emphasizing the diseases.
To explore the impact of multimorbidity on affected patients and the consequences for health services, we designed a study to describe multimorbidity by sex and life-stage in a large population sample and to assess the association with acute morbidity, area of residency and use of health services.
Data source and study population
Cross-sectional study of adults resident in Catalonia, a Mediterranean region of southern European with 7,434,632 inhabitants (2010 census), 16% of the population of Spain. In Catalonia, 358 primary health care teams (PHCT) comprised of doctors, nurses, social workers and support staff are assigned by geographical area and responsible for the health care of the population in their areas. The Catalan Health Institute (CHI) manages 274 PHCT (76.5%), serving a population of 4,859,725 adults; the remaining PHCT are managed by other providers. Doctors and nurses systematically use electronic health records (EHR) to record diagnoses, prescriptions and other clinical, patient management and administrative activities. The CHI Information System for the Development of Research in Primary Care (SIDIAP)  compiles coded clinical information from the EHR system. A subset of SIDIAP records meeting the highest quality criteria for clinical data (SIDIAP-Q) includes 40% of the SIDIAP population (1,936,443 patients), attended by the 1,319 general practitioners (GP) whose data recording scored highest in a validated comparison process. The sample is representative of the general Catalan population in terms of geography, age and sex distributions, according to the official 2010 census .
A sample of 1,749,710 patients aged 19 years or older, assigned to 251 PHCT during the period of study (1 January- 31 December 2010), was selected from the SIDIAP-Q database.
Coding of diseases
International Classification of Diseases (ICD-10) codes were mapped to the International Classification of Primary Care (ICPC-2e-v.4.2, available at: http://www.kith.no/templates/kith_WebPage____1111.aspx). R codes (symptoms, signs and abnormal clinical and laboratory findings, not elsewhere classified) and Z codes (factors influencing health status and contact with health services) were excluded, resulting in 686 included codes. Each diagnosis was then classified using O’Halloran criteria for chronic disease . We included all 146 diagnoses considered as chronic diseases by O’Halloran criteria: (i) have a duration that has lasted, or is expected to last, at least 6 months; (ii) have a pattern of recurrence or deterioration; (iii) have a poor prognosis and (iv) produce consequences, or sequelae, that have an impact on the individual’s quality of life [14, 15]. Any disease not meeting the O’Halloran criteria was considered an acute disease.
All results were described with ICPC-2 codes. Diseases were classified as acute if diagnosed during the study period and chronic if recorded as such in EHR as of 31 December 2010.
Outcomes and variables
The main outcome was multimorbidity, defined as the coexistence of 2 or more chronic diseases. Secondary outcome was the number of new events of each acute disease. Other variables recorded for each patient were the following: number of all acute diseases (0, 1, 2, > = 3), sex (male, female), age (young adult, 19 to 24; adult, 25–44; older adult, 45–64; elderly, 65–79; and oldest adult, 80+), number of visits during the study period (0, 1–2, 3–5, 6–10, ≥11), and area of residence (rural if <10,000 inhabitants and/or population density <150 people/km2, otherwise urban) . Number of all acute diseases (or 0 diseases) and visits (or 0 visits) were categorized as quartiles of the study population. Number of visits was used as a proxy of use of health services and included visits recorded in EHR by GP, nurses or social workers, either at the primary care centre or as home health care.
Analysis was stratified by sex and age group. Descriptive statistics were used to summarize overall information. Categorical variables were expressed as frequencies (percentage) and continuous as mean (Standard deviation, SD) or median (interquartile range, IQR).
Cumulative incidence of acute morbidity events was calculated as the number of new acute events during the study period divided by the at-risk population in the sample (e.g., if a patient had bronchitis twice in the one-year study period, the total number of events accounted for was 2). We took into account the five acute diseases with the highest cumulative incidence within each stratum. Risk ratios (RR) of multimorbidity vs. non-multimorbidity were calculated for the number of events for each acute disease, using Poisson, negative binomial (if overdispersion was present) or zero inflated (when data had an excess of zero counts) equations, as appropriate. All models were adjusted for number of visits and area of residency.
To determine the most prevalent multimorbidity patterns, all possible combinations of any two chronic diseases and their frequencies were calculated. Observed (O) and expected (E) prevalence of those two chronic diseases with each acute disease was then computed. Expected co-occurrence of diseases was obtained as the product of these prevalences, assuming statistical independence of the diseases concerned. The overlapping of those combinations that presented the highest O/E ratio was reported.
Logistic regression was used to assess the association between multimorbidity and the variables listed above.
All statistical tests were two-sided at the 5% significance level. The analyses were performed using SPSS for Windows, version 18 (SPSS Inc., Chicago, IL, USA), Stata/SE version 11 for Windows (Stata Corp. LP, College Station, TX, USA) and R version 2.15.2 (R Foundation for Statistical Computing, Vienna, Austria).
The study protocol was approved by the Committee on the Ethics of Clinical Research, Institut Universitari d’Investigació en Atenció Primària (IDIAP) Jordi Gol (Protocol No: P12/28). All data were anonymized and the confidentiality of EHR was respected at all times in accordance with international law.
We included 1,749,710 patients; mean age was 47.4 years (SD: 17.8), 50.7% were female and 16% lived in rural areas. Multimorbidity (≥ 2 diseases) was present in 46.8% (95%CI, 46.7%-46.8%) of the sample, being higher in female (52.3%) than in male (41.1%) and in rural areas (47.6%) than in urban areas (46.6%).
The prevalence of the most common chronic diseases differed by sex below 45: anxiety disorder/anxiety (women aged 19–44); acne (men aged 19–24) and lipid disorder (men aged 25–44). After 45 both sex groups first chronic disease was lipid disorder in 45–64 and uncomplicated hypertension in 65-80+ (Table 1). Upper respiratory infection acute is the most incident acute disease in all age groups (except in 80+).
Multimorbidity prevalence increased as age increased, being higher in female (ranged from 19.0% to 92.1%) than male (ranged from 12.9% to 92.0%). In patients with multimorbidity, the number of acute diseases was higher in female than male and decreased as age increased, except in male older than 65. In addition, the number of visits increased as age increased, and was higher for female than male in all age groups except 80+ (Table 2).
Patients with multimorbidity had a higher incidence of acute diseases and number of visits in all age strata than non-multimorbidity patients; in both cases, the incidence was higher for female than male (Table 2). Overall, the median (IQR) of number of visits was 8(4–14) in patients with multimorbidity vs. 1(0–4) in the non-multimorbidity group.
The two most prevalent combinations of two chronic diseases were hypertension and lipid disorders in patients older than 45 years. The only acute disease that appeared in both sexes was “bursitis/tendinitis/synovitis NOS” in the oldest age group (80+). In the other age groups, the acute health disease varied by sex (Figure 1).
The five acute diseases with the highest cumulative incidence were similar by sex in any age group. Infections were the most common diagnosis. Cystitis/urinary infection was present among the five most prevalent acute conditions only in women and in the oldest men. In women, the highest significant RR of multimorbidity vs. non-multimorbidity was found for teeth/gum disease (aged 19–24) and upper respiratory infection, acute (80+). In men, this RR was only positive and significant for teeth/gum disease (aged 65–79) (Table 3).
The adjusted analysis of factors associated with multimorbidity showed that the oldest patients (80+ years) with acute diseases and women aged 65–79 with 3 or more acute diseases were more likely to have multimorbidity than patients with no acute diseases. This positive association was only significant in women. Living in a rural area was significantly associated with lower probability of having multimorbidity. Patients who visited a GP more often were more likely than those without visits to have multimorbidity, reaching the highest probability in those aged 65–79 years (Table 4).
Statement of principal findings
Almost half of the study population had multimorbidity, with infections (mainly acute upper respiratory infection) the most common acute disease in both sexes and all age groups. The most frequent multimorbidity pattern of chronic diseases was the combination of hypertension and dyslipidemia in adults over 45 years of age.
We observed a decrease in the number of acute diseases recorded as age increased. Nonetheless, in adjusted models female older than 65 who had acute diseases were more likely to have multimorbidity.
Finally, the use of health services was positively associated with a diagnosis of multimorbidity. Living in a rural area decreased the probability of multimorbidity.
Strengths and weaknesses of the study
A major strength of this study is the analysis of a large, high-quality database of primary-care records, representative of a large population. In the context of a national health system with universal coverage, EHR data have been shown to yield more reliable and representative conclusions than those derived from survey-based studies . Another important strength was the inclusion of all chronic and acute diagnoses registered in EHR, which contributed to a more accurate analysis of the association between acute and chronic diseases and of the disease combinations present in multimorbidity in this population. To synthesize the results, we present here only the most frequent combinations. Finally, few studies have incorporated acute diseases in the study of multimorbidity patterns  and none analyzed the relationship between multimorbidity and acute morbidity.
Some possible biases could have influenced our results. First, diseases could be underreported, especially for male of normal workforce age who tend to see their doctors less often than other strata of patients. This effect would diminish in the two oldest age groups because of the retirement age (65 years) in Spain. In patients with multimorbidity, the true incidence of acute disease could be underreported because the GP would place higher priority on the chronic problems in these patients. On the other hand, there could be an over-representation of chronic diagnoses (e.g., hypertension, diabetes, hyperlipidemia, etc.) that are included in the goals/incentives contracts of Catalan PHCT.
The diseases that form part of the CatSalut treatment objectives may be more carefully recorded than other conditions. However, these same diseases are the most prevalent (high blood pressure, diabetes, hypercholesterolemia, smoking, dyslipidemia, ischemic heart disease, atrial fibrillation) and therefore have the greatest impact on population health. The quality-recorders database (SIDIAP-Q) used in this study minimizes the under-reporting of diseases not included in the CatSalut objectives.
Furthermore, the stratified analysis allows more accurate estimation within each age-sex strata and universal access to free health care and medications makes it more likely that patients seeking care will acquire a diagnosis, either acute or chronic [17, 18]. Second, there is no universally accepted criterion for consensus classification of acute and chronic disease. This lack of accurate case definitions impedes the establishment of the true incidence/prevalence of a disease . Finally, a residual confounding cannot be completely excluded, and could occur because of epidemiological factors not considered in this study, such as patients’ socioeconomic status .
Strengths and weaknesses in relation to other studies
The estimated multimorbidity prevalence in our sample is higher than in other European studies [21–24], perhaps because of the analysis of a greater number of diseases in our study than in most other published studies . Nonetheless, the patterns of multimorbidity observed were similar to those observed in other studies .
As in other studies, multimorbidity was more prevalent among female [21, 22, 24]. This could be due to the longer female life expectancy and worse health status, compared to male, differences that are due to both biological and social factors . In addition, sex is a social determinant that influences health status, health behaviours and the use of health services [27–30]. Recent studies, however, suggest a dismantling of this paradigm based on sex-stratified analysis of consultations for common symptoms .
Acute problems are time-consuming for health professional , and therefore should be considered part of the primary care workload. Although current health policy, health care services, and research are all heavily focused on chronic diseases, we must not forget that 41% of primary care visits are motivated by an acute disease . The incidence of acute diseases observed in our study concurs with other reports of acute upper respiratory infection and other health problems of infectious aetiology (acute tonsillitis, cystitis) as the primary reasons for seeking primary care, along with non-infectious diseases such as dorsalgia [33, 34].
Our study observed a higher prevalence of multimorbidity in rural settings. Other studies conducted in rural areas have reported only a greater prevalence of multimorbidity in elderly people [6, 7]. Nonetheless, our adjusted analysis showed that living in a rural area is negatively associated with multimorbidity. This phenomenon could be due to the environmental and sociocultural context and access to both public and private services .
Implications for clinicians and policymakers
Our study considered multimorbidity in patients who received primary medical care, considering all visits and diagnoses (acute and chronic diseases). This approach allowed the identification of vulnerable subgroups in our population. A major advantage of our methodology is the use of data obtained directly from standard clinical practice. Knowing the distribution of acute and chronic diseases by life-stage and sex will help the clinician faced with a particular patient to anticipate disease patterns based on the patient’s sex and stage of life, recognizing that these vary with age. This will encourage the implementation of personalized disease prevention and health promotion activities.
At the level of health policy and health care administration, the organization of services should be reviewed to ensure that continuity and coordination of patient care are guaranteed; current evidence suggests the potential for improvement in this regard .
Unanswered questions and future research
The classification of chronic and acute disease remains unresolved, and there is no consensus on the type and number of chronic diseases that define multimorbidity. A personalized measure to determine the severity of diagnosed multimorbidity is also needed.
If longitudinal studies confirm a higher incidence of morbidity in patients with multimorbidity, evidence-based interventions will be needed to prevent the onset of acute disease. Further studies are needed to study possible genetic and pathophysiological explanations that corroborate the observed multimorbidity patterns.
In-depth analysis of other contextual factors related to multimorbidity is also required, along with studies of the relationship between area of residence and multimorbidity and of the differences in health status that may exist between different territories. Finally, there is a need for the implementation and evaluation of health literacy and self-management interventions to improve patient competence in resolving routine acute diseases, which in turn will decrease the care burden in primary care systems.
Multimorbidity is related to more use of health services and higher incidence of acute diseases, which increases the burden on primary care services. Residence in urban vs rural settings is a factor for future in-depth study.
The association of acute morbidity, area of residency or use of health services with multimorbidity differs according to life-stage and sex. Therefore, the study of multimorbidity should be stratified by life-stage and sex. Understanding these trends across life-stages will allow health systems to adjust their clinical and management models to adapt and prioritize interventions.
Primary health care teams
Catalan Health Institute
Electronic health records
Information System for the Development of Research in Primary Care
International Classification of Diseases
International Classification of Primary Care
Institut Universitari d’Investigació en Atenció Primària
Bodenheimer T, Wagner EH, Grumbach K: Improving primary care for patients with chronic illness. JAMA. 2002, 288: 1775-1779.
Van den Akker M, Buntinx F, Knottnerus JA: Comorbidity or multimorbidity: what’s in a name? A review of literature. Eur J Gen Pract. 1996, 2: 65-70.
Freund T, Kunz CU, Ose D, Szecsenyi J, Peters-Klimm F: Patterns of multimorbidity in primary care patients at high risk of future hospitalization. Popul Health Manag. 2012, 15: 119-124.
Gijsen R, Hoeymans N, Schellevis FG, Ruwaard D, Satariano WA, van den Bos GA: Causes and consequences of comorbidity: a review. J Clin Epidemiol. 2001, 54: 661-674.
van den Akker M, Buntinx F, Roos S, Knottnerus JA: Problems in determining occurrence rates of multimorbidity. J Clin Epidemiol. 2001, 54: 675-679.
Khanam MA, Streatfield PK, Kabir ZN, Qiu C, Cornelius C, Wahlin A: Prevalence and patterns of multimorbidity among elderly people in rural Bangladesh: a cross-sectional study. J Health Popul Nutr. 2011, 29: 406-414.
John R, Kerby DS, Hennessy CH: Patterns and impact of comorbidity and multimorbidity among community-resident American Indian elders. Gerontologist. 2003, 43: 649-660.
Saltman DC, Sayer GP, Whicker SD: Co-morbidity in general practice. Postgrad Med J. 2005, 81: 474-480.
Broemeling A, Watson D, Black C: Chronic conditions. Co-morbidity among residents of British Columbia. 2005, Vancouver, BC: University of British Columbia
World Health Organization: A Glossary of Terms for Community Health Care and Services for Older Persons. 2004, Kobe: WHO
Wolff JL, Starfield B, Anderson G: Prevalence, expenditures, and complications of multiple chronic conditions in the elderly. Arch Intern Med. 2002, 162: 2269-2276.
Information system for the development of research in primary care (SIDIAP data base). [http://www.sidiap.org/]
García-Gil Mdel M, Hermosilla E, Prieto-Alhambra D, Fina F, Rosell M, Ramos R, Rodriguez J, Williams T, Van Staa T, Bolibar B: Construction and validation of a scoring system for the selection of high-quality data in a Spanish population primary care database (SIDIAP). Inform Prim Care. 2011, 19 (3): 135-145.
O’Halloran J, Miller GC, Britt H: Defining chronic conditions for primary care with ICPC-2. Fam Pract. 2004, 21: 381-386.
Defining Chronic Conditions for Primary Care Using ICPC-2. Available in: http://www.fmrc.org.au/Download/DefiningChronicConditions.pdf
Violán C, Foguet-Boreu Q, Hermosilla-Pérez E, Valderas JM, Bolíbar B, Fàbregas-Escurriola M, Brugulat-Guiteras P, Muñoz-Pérez MÁ: Comparison of the information provided by electronic health records data and a population health survey to estimate prevalence of selected health conditions and multimorbidity. BMC Public Health. 2013, 13: 251-
Jordan K, Porcheret M, Croft P: Quality of morbidity coding in general practice computerized medical records: a systematic review. Fam Pract. 2004, 21: 396-412.
Valderas JM, Starfield B, Sibbald B, Salisbury C, Roland M: Defining comorbidity: implications for understanding health and health services. Ann Fam Med. 2009, 7: 357-363.
de Lusignan S, Tomson C, Harris K, van Vlymen J, Gallagher H: Creatinine fluctuation has a greater effect than the formula to estimate glomerular filtration rate on the prevalence of chronic kidney disease. Nephron Clin Pract. 2011, 117: c213-c224.
Barnett K, Mercer SW, Norbury M, Watt G, Wyke S, Guthrie B: Epidemiology of multimorbidity and implications for health care, research, and medical education: a cros-sectional study. Lancet. 2012, 380: 37-43.
Prados-Torres A, Poblador-Plou B, Calderón-Larrañaga A, Gimeno-Feliu LA, González-Rubio F, Poncel-Falcó A, Sicras-Mainar A, Alcalá-Nalvaiz JT: Multimorbidity patterns in primary care: interactions among chronic diseases using factor analysis. PLoS One. 2012, 7: e32190-
Glynn LG, Valderas JM, Healy P, Burke E, Newell J, Gillespie P, Murphy AW: The prevalence of multimorbidity in primary care and its effect on health care utilization and cost. Fam Pract. 2011, 28: 516-523.
Salisbury C, Johnson L, Purdy S, Valderas JM, Montgomery AA: Epidemiology and impact of multimorbidity in primary care: a retrospective cohort study. Br J Gen Pract. 2011, 61: e12-e21.
García-Olmos L, Salvador CH, Alberquilla Á, Lora D, Carmona M, García-Sagredo P, Pascual M, Muñoz A, Monteagudo JL, García-López F: Comorbidity patterns in patients with chronic diseases in general practice. PLoS One. 2012, 7: e32141-
Marengoni A, Angleman S, Melis R, Mangialasche F, Karp A, Garmen A, Meinow B, Fratiglioni L: Aging with multimorbidity: a systematic review of the literature. Ageing Res Rev. 2011, 10: 430-439.
Laux G, Kuehlein T, Rosemann T, Szecsenyi J: Co- and multimorbidity patterns in primary care based on episodes of care: results from the German CONTENT project. BMC Health Serv Res. 2008, 8: 14-
Lee JH, Sadana R: Improving Equity in Health by Addressing Social Determinants. The Commission on Social Determinants of Health Knowledge Networks. 2011, World Health Organization, [http://whqlibdoc.who.int/publications/2011/9789241503037_eng.pdf]
Marmot M, Friel S, Bell R, Houweling TA, Taylor S: Commission on social determinants of health: closing the gap in a generation: health equity through action on the social determinants of health. Lancet. 2008, 372: 1661-1669.
Bertakis KD, Azari R, Helms LJ, Callahan EJ, Robbins JA: Gender differences in the utilization of health care services. J Fam Pract. 2000, 49: 147-152.
Redondo-Sendino A, Guallar-Castillón P, Banegas JR, Rodríguez-Artalejo F: Gender differences in the utilization of health-care services among the older adult population of Spain. BMC Public Health. 2006, 6: 155-
Hunt K, Adamson J, Hewitt C, Nazareth J: Do women consult more than men? A review of gender and consultation for back pain and headache. J Health Serv Res Policy. 2011, 16: 108-117.
Stange KC, Zyzanski SJ, Jaén CR, Callahan EJ, Kelly RB, Gillanders WR, Shank JC, Chao J, Medalie JH, Miller WL, Crabtree BF, Flocke SA, Gilchrist VJ, Langa DM, Goodwin MA: Illuminating the ‘black box’. A description of 4454 patient visits to 138 family physicians. J Fam Pract. 1998, 46: 377-389.
Woodwell DA, Cherry DK: National ambulatory medical care survey: 2002 summary. Adv Data. 2004, 346: 1-44.
Wändell P, Carlsson AC, Wettermark B, Lord G, Cars T, Ljunggren G: Most common diseases diagnosed in primary care in Stockholm, Sweden, in 2011. Fam Pract. 2013, 506: 13-
Smith SM, Soubhi H, Fortin M, Hudon C, O’Dowd T: Managing patients with multimorbidity: systematic review of interventions in primary care and community settings. BMJ. 2012, 345: e5205-
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2296/15/55/prepub
We thank the Catalan Health Institute and especially the SIDIAP Unit, which provided the database for the study. The authors also appreciate the English language review by Elaine Lilly, PhD.
This study was funded by the Ministry of Science and Innovation through the Instituto Carlos III (ISCiii) as part of the Primary Care Prevention and Health Promotion Research Network (rediAPP), by ISCiii-RETICS (RD12/0005), by a grant for research projects ISCiii (PI12/00427), and by a 2011–2013 scholarship that aims to promote research in Primary Health Care by health professionals who have completed their specialty training, awarded by Institut Universitari d’Investigació en Atenció Primària Jordi Gol (IDIAP Jordi Gol). The funders had no role in the study design, collection, analysis and interpretation of data, writing of the manuscript and decision to submit for publication.
The authors declare that they have no competing interests.
All authors contributed to the design of the study, revised the article, and approved the final version. CV, QFB, JMV, MMP, ARL drafted the study protocol and obtained the funding. TRB, CV, QFB, JMV, ARL contributed to the analysis and interpretation of data. CV, QFB, JMV, TRB, ARL, MPV, MMP, and EPR wrote the first draft, and all authors contributed ideas, interpreted the findings and reviewed rough drafts of the manuscript. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.