Study protocol: Identification and validation of integrative biomarkers of physical activity level and health in children and adolescents (INTEGRActiv)
Catalina Picó1,2,3 
Empar Lurbe3,4 
Jaap Keijer5 
Jan Kopecky6 
Jean-François Landrier7 
Julio Álvarez-Pitti3,4 
Jean Charles Martin7 
Paula Oliver1,2,3* 
Andreu Palou1,2,3 
Mariona Palou1,2,3 
Petr Zouhar6 
Joan Ribot1,2,3 
Ana M. Rodríguez1,2,3 
Juana Sánchez1,2,3 Francisca Serra1,2,3 
M. Luisa Bonet1,2,3
- 1Laboratory of Molecular Biology, Nutrition and Biotechnology (Group of Nutrigenomics, Biomarkers and Risk Evaluation), University of the Balearic Islands (UIB), Palma, Spain
- 2Health Research Institute of the Balearic Islands (IdISBa), Palma, Spain
- 3CIBER de Fisiopatología de la Obesidad y Nutrición (CIBEROBN), Madrid, Spain
- 4Department of Pediatrics (Innovation in Paediatrics and Technologies-iPEDITEC- research group), Fundación de Investigación, Consorcio Hospital General, University of Valencia, Valencia, Spain
- 5Human and Animal Physiology, Wageningen University, Wageningen, Netherlands
- 6Laboratory of Adipose Tissue Biology, Institute of Physiology of the Czech Academy of Sciences, Prague, Czech Republic
- 7C2VN, INSERM, INRAE, Aix-Marseille University, Marseille, France
Background: Physical activity (PA) provides health benefits across the lifespan and improves many established cardiovascular risk factors that have a significant impact on overall mortality. However, discrepancies between self-reported and device-based measures of PA make it difficult to obtain consistent results regarding PA and its health effects. Moreover, PA may produce different health effects depending on the type, intensity, duration, and frequency of activities and individual factors such as age, sex, body weight, early life conditions/exposures, etc. Appropriate biomarkers relating the degree of PA level with its effects on health, especially in children and adolescents, are required and missing. The main objective of the INTEGRActiv study is to identify novel useful integrative biomarkers of PA and its effects on the body health in children and adolescents, who represent an important target population to address personalized interventions to improve future metabolic health.
Methods/design: The study is structured in two phases. First, biomarkers of PA and health will be identified at baseline in a core cohort of 180 volunteers, distributed into two age groups: prepubertal (n = 90), and postpubertal adolescents (n = 90). Each group will include three subgroups (n = 30) with subjects of normal weight, overweight, and obesity, respectively. Identification of new biomarkers will be achieved by combining physical measures (PA and cardiorespiratory and muscular fitness, anthropometry) and molecular measures (cardiovascular risk factors, endocrine markers, cytokines and circulating miRNA in plasma, gene expression profile in blood cells, and metabolomics profiling in plasma). In the second phase, an educational intervention and its follow-up will be carried out in a subgroup of these subjects (60 volunteers), as a first validation step of the identified biomarkers.
Discussion: The INTEGRActiv study is expected to provide the definition of PA and health-related biomarkers (PA-health biomarkers) in childhood and adolescence. It will allow us to relate biomarkers to factors such as age, sex, body weight, sleep behavior, dietary factors, and pubertal status and to identify how these factors quantitatively affect the biomarkers’ responses. Taken together, the INTEGRActiv study approach is expected to help monitor the efficacy of interventions aimed to improve the quality of life of children/adolescents through physical activity.
Clinical Trial Registration: ClinicalTrials.gov, Identifier NCT05907785.
1. Introduction
Physical activity (PA) provides significant health benefits across the lifespan and improves many of the established cardiovascular risk factors, such as hypertension, obesity, and diabetes, which have a significant negative impact on overall mortality (1, 2). Although there is strong evidence of the positive relationship between PA and health, it is not very clear what kind and how much PA is necessary, nor can precise individual PA recommendations for the maintenance of health be made. There has been a high heterogeneity both in the definition of PA and its measurement (3). PA level is generally determined from self-report questionnaires, which are subject to potential bias, or through more objective measures, such as pedometers and accelerometers, but there are differences in the way intensity, duration, and frequency of activities are reported that make it difficult to obtain clear, consistent results (2, 3). There is an important need for markers to measure PA more precisely and objectively and, more specifically, markers allowing to relate the degree of PA with effects on the body and health status. PA may produce different health effects depending on type, intensity, duration, and frequency of activities, and also depending on individual factors such as age, sex, body weight, pre-existing health conditions, lifestyle, early life conditions/exposures, etc (4, 5). Not only may such missing biomarkers be a useful tool to measure specific health effects of PA, but they may also be essential for a deeper insight into the complex biological relationships that exist between physical fitness and health determinants.
Most studies dealing with the identification of biomarkers of PA and its relation to health have been performed in adults and athletes [e.g., studies (6, 7)]. Although diseases such as coronary artery disease, hypertension, and type 2 diabetes are generally manifested in adult life, the processes underlying the development of these diseases generally originate in early ages, and childhood and adolescence are critical periods (8). The general WHO recommendation of an average of 60 min/day of moderate-to-vigorous PA in children and adolescents (9) is inversely associated with biomarkers of cardiometabolic risk (10), including lower obesity rates (11). However, globally, it is recognized that children spend an insufficient amount of time in physical activities and too much time in sedentary activities (11). The availability of biomarkers indicative of PA and its relationship with health status (PA-health biomarkers) in these sensitive age groups, will be very useful to make precise, more personalized recommendations to improve health and prevent disease development. The development of PA-health biomarkers will be addressed by the INTEGRActiv study.
Because of the multifaceted nature of homeostasis, omics technologies, such as transcriptomics and metabolomics that broadly analyse functional genomic responses may be particularly valuable for the identification and characterization of biomarkers of PA on health. Peripheral blood cells (PBC) or the subpopulation of mononuclear cells within them (PBMC, including lymphocytes and monocytes) are an attractive source of biomarkers, particularly transcriptomic-based biomarkers because they can be easily obtained by minimally invasive techniques (12, 13). These cells travel through the blood, respond to external and internal factors, and their gene expression profile can partly reflect the expression profile of other tissues and be indicative of metabolic, physiological, and pathological states, as suggested by studies in animal models (14–17) and in humans (18–20). Moreover, PBMC were recently shown to be able to functionally reflect differences in fitness in adults (21). Hence, changes in gene expression in blood cells may inform on physiopathological states and may also have a predictive component (18, 22, 23). The use of PBC as a source of transcriptome markers concerning PA in children/adolescents is almost unexplored so far. In a recent study in children, we have described changes in PBC gene expression associated with low PA, which may be related to the cardiometabolic health effects of PA (24). The INTEGRActiv study will help to validate levels of our previously identified transcripts as indicators of low PA and to identify other potential biomarkers by using targeted and untargeted analyses. Functional links of the encoded transcripts/proteins with cardiovascular and endocrine markers may give information on the effects of PA on body health and will be further explored through literature and data mining/integration.
Among the transcript-based biomarkers, microRNAs (miRNAs) represent specific species of non-coding RNAs that post-transcriptionally regulate gene expression. They are found inside blood cells, but also in the circulation, where they remain stable and can be easily measured (25). PA has been shown to affect levels of specific circulating miRNAs in adult subjects (26), but further research is needed to define them as useful biomarkers. Therefore, circulating miRNAs constitute a new potential regulatory component that may play a role in exercise-induced adaptations and be of interest as PA biomarkers (27, 28). Interestingly, circulating miRNAs also have a great potential as cardiac biomarkers. Therefore, changes in circulating miRNAs could provide further insight into the molecular physiological response to PA, as well as valuable information on the link between PA and cardiovascular health (26, 29, 30) and may thus serve as PA-health biomarkers.
Metabolomics is also a promising approach to exploring the relationship between PA and health, a relationship that can differ between groups and individuals. Extensive plasma metabolite profiling can provide an overview of the metabolism with a level of description that goes beyond genetic information and more closely reflects the phenotype, thus helping to connect genotype to phenotype at the molecular level (31). There is no question that PA affects the metabolome [e.g., (7)], but it is necessary to better understand when these changes occur and for how long they persist. It would also be desirable to identify those biomarkers of PA that indicate and predict metabolic health, including cardiovascular health and propensity to obesity (32). To identify such biomarkers, we will conduct comprehensive untargeted metabolomics and lipidomics analysis on a wide range of low-molecular-weight analytes. In our previous studies in mice within the European project BIOCLAIMS, we learned that plasma acylcarnitines (markers of lipid catabolism) could also serve as a very early and sex-specific biomarker of obesity propensity (33). Here, we will characterize these and other analytes, such as the novel biomarker of PA, N-lactoyl-phenylalanine (34).
Secreted regulatory molecules are also of interest as potential biomarkers and mediators of PA responses. Among them, myokines are signaling molecules released by muscles and involved in the crosstalk of muscle with other organs/tissues (35). Circulating myokines have been associated with muscle contraction and with multiple metabolic effects, such as cardiovascular and anti-inflammatory effects, and the promotion of energy metabolism and “browning” in white adipose tissue, most of them related to the beneficial effects of PA (36). Interestingly, adipose tissue is targeted by myokines and can respond to myokines by changing the production of specific adipose-born cytokines, the so-called adipokines. This myokine/adipokine crosstalk is very relevant in modulating molecular events to ensure whole-body metabolic homeostasis (37). Hence, the study of selected myokines and adipokines is planned as potential PA-health biomarkers.
Thus, the INTREGRActiv study focuses on children and adolescents and will assess PA and relate PA to cardiorespiratory and muscular fitness, measures of body weight and body fat, cardiovascular risk factors and endocrine markers, cytokines, circulating miRNA, gene expression profile in blood cells, and metabolomics profile in plasma to derive new, useful PA-health biomarkers (Integrative biomarkers of PA and health) in children and adolescents. Data regarding perinatal factors, food consumption frequency, sleep behavior, and other lifestyle and social determinants may also have a great impact on health and PA responses and will also be considered.
2. Methods and analysis
2.1. Aim and setting of the study
The main objective of the INTEGRActiv study is to identify novel useful integrative biomarkers of PA level and its effects on body health (PA-health biomarkers) in children and adolescents and to carry out a first validation stage. The identification of novel candidate biomarkers will be achieved by combining measures of PA and cardiorespiratory and muscular fitness assessment with anthropometric measures, measures of cardiovascular risk factors, endocrine markers, cytokines, and circulating miRNA in plasma, gene expression profile in blood cells, and metabolomics profile in plasma in prepubertal and pubertal children/adolescents in different categories of BMI z-score. Besides the influence of BMI, age, pubertal state, and gender, possible confounding factors such as diet and lifestyle-related factors including sleep behaviour will be considered. It is expected that newly identified biomarkers reflecting PA level and its relation with health may guide clinical and nutritional/lifestyle advice and public policies related to endorsing personalized PA, with better adherence and response, to promote health and prevent disease risk factors since early stages of life.
2.2. Design of the study
The INTEGRActiv study is structured in two phases. The first baseline phase (phase 1) is to identify biomarkers of PA and health in a new core cohort of 180 volunteers to be established. In the second phase (phase 2), an educational intervention and its follow-up will be carried out in a subgroup of these subjects (60 volunteers), to allow a first validation step based on the evolution of the identified biomarkers. The general flowchart of the study is illustrated in Figure 1.
Figure 1. The general flowchart of the INTEGRActiv study.
2.3. Establishment of the core cohort, data collection, and procedures
2.3.1. Participants
A total of 180 volunteers (core cohort) will be recruited and distributed into two age groups: prepubertal individuals (n = 90), including children, from 8 to 10 years, and postpubertal adolescents (n = 90) from 13 to 15 years. Each group will include three subgroups (n = 30) of subjects with normal weight (NW), overweight (OW), and obesity (OB) as defined by WHO criteria (38). Volunteers will be equally distributed between both sexes in each group and subgroup. The sample size of the core cohort is based on our previous expertise in different projects aimed at identifying biomarkers (mainly gene expression biomarkers analysed in blood cells by omics techniques), in which we used even a lower number of volunteers [e.g., studies (24, 39)]. Moreover, the combination/integration of different types of biomarkers that will be performed (based not only on gene expression, but also considering lipidomic, metabolic and miRNA profiles) makes its more likely to be able to detect significant metabolic differences and more robust biomarkers. Additionally, the number of volunteers of the core cohort is also expected to be sufficient to achieve the required sample size for the intervention study, which was calculated as explained in Section 2.4.
Recruitment of children/adolescents will be performed through public advertisement, among the patients of the Cardiovascular Risk Unit of the Pediatric Department of the General University Hospital of Valencia (Spain), among patients of primary care pediatricians, as well as in recruitment campaigns in sports academies and summer camps. Children/adolescents will be informed of the study, and those who express interest and meet the specific requirements of the different groups in the core cohort will be selected for recruitment. Priority will be given to individuals expressing interest in the subsequent intervention to guarantee the minimum number of volunteers required. Pubertal assessment will be performed by experts in pediatric endocrinology. Children/adolescents with regular medication, reported pathologies, and/or the presence of eating disorders will be excluded to avoid potential confounding factors. Individuals with clinical recommendations against the performance of PA will be also excluded. Screening of eating disorders will be done using the SCOFF and DEBQ questionnaires (40, 41). During recruitment, all volunteers will be informed that they will have the possibility to participate in a 6-month duration educational intervention (see below).
2.3.2. Data collection and procedures
Each volunteer will be given three initial appointments to obtain relevant information and blood samples, and to perform different procedures (detailed information of the procedures developed in each appointment is provided in Table 1).
Table 1. Information of the procedures developed in each appointment. .
2.3.2.1. Anthropometric parameters
Weight, height, and waist circumference will be measured by trained nurses to derive BMI, the corresponding standard deviation (SD) and BMI Z score, and waist-to-height values as useful indicators of abdominal obesity and cardiometabolic risk in children. Subjects with a BMI Z score below 1 will be classified as presenting normal-weight; those ranging from 1 to 2 will be defined as subjects with overweight, while obesity is defined by a BMI + 2 SD value (38). Waist circumference will be measured at the midpoint between the iliac crest and the costal margin at the midaxillary line in the standing position at the end of a gentle expiration, and waist percentiles will be calculated according to growth curves (42).
2.3.2.2. Sociodemographic and clinical variables/information
Data regarding perinatal factors (e.g., birth weight and birth complications), breastfeeding duration, weight gain during the first 5 years of life, eating habits during childhood and feeding difficulties, assessment of pubertal stage, and date of menarche and menstrual cycles (in women) will be collected through a clinical interview. Personal and family health history, leisure time activities, family structure and peers, school setting, family affluence, social inequality, and parental weight/height and age will also be recorded.
2.3.2.3. Questionnaires on dietary intake, PA, sedentary and sleep habits
Dietary intake will be assessed by an expert nutritionist following the Guidance on the EU Menu methodology (43); adherence to the Mediterranean diet will be evaluated using the KIDMED questionnaire (44, 45). PA level will be evaluated using the PA Questionnaire for Older Children and Adolescents (PAQ-C, PAQ-A), which has been previously used to classify them into different activity levels (46, 47) and to investigate the PA and health outcomes relationship (46, 48). Sedentary and screen time (49) and sleep habits (50) will be assessed using validated questionnaires. All questionnaires are in Spanish. For each participant, questionnaire responses will be completed online (with the support of parents and investigators) and automatically recorded in the software database (Lyme Survey).
2.3.2.4. Blood sample collection
PBC will be obtained via antecubital fossa venipuncture in the early morning, under fasting conditions. Blood samples will be collected in PAX tubes, to allow optimal conservation and extraction of RNA from PBC for gene expression analyses, and in tubes containing or not anticoagulants (heparin) to obtain plasma and serum, respectively, for analyses of biochemical parameters, miRNA, cytokines, and metabolomics. Circulating biochemical parameters analyzed will include those classically related to metabolic syndrome; abnormal values will be defined from normative data (51).
2.3.2.5. Blood pressure (BP) measurements
BP will be measured three consecutive times by nurses using a validated monitor (Omron-M705) oscillometric recorder, according to the European Society of Hypertension (ESH) Guidelines (52). Volunteers will be classified as normotensive, high-normal, stage 1, or stage 2 hypertensive according to the criteria of the ESH Guidelines, in which the fourth report normative values were used for interpretation in adolescents <16 years (52).
2.3.2.6. Cardiorespiratory fitness (CRF) assessment
CRF refers to the capacity of the circulatory and respiratory systems to supply oxygen to skeletal muscle mitochondria for energy production needed during PA. Low or unhealthy CRF is a strong, independent predictor of cardiovascular disease (CVD) and all-cause mortality in adults (53). In youth, CRF is a predictor of several health indicators, including cardiometabolic health (54, 55), premature CVD (56), academic achievement (57), and mental health (58). CRF will be measured as maximal oxygen consumption (V˙O2max) from cardiopulmonary exercise testing, using a treadmill test and the Balke protocol (59) modified to enhance a secure environment, i.e., avoiding running in steep slopes and taking subjects to exhaustion, especially for the youngest population and those volunteers with moderate or severe obesity (60).
2.3.2.7. Muscular fitness assessment
Muscular fitness will be assessed following the recommendations of the ALPHA Health-Related Fitness Test Battery for Children and Adolescents. The handgrip strength and the standing broad jump tests will be used to assess musculoskeletal fitness (61).
2.3.2.8. Objective assessment of gait and PA
Triaxial accelerometers (Actigraph model wGT3X, Manufacturing Technology Inc, Fort Walton Beach, USA) will be used for objective measurement of PA and to detect gait disturbances. The accelerometer will be fixed at the waist above the right iliac crest. (1) In the laboratory, the subject will be asked to walk on a treadmill for 10 min adjusting the speed to walking pace, without interfering with his/her usual walking pattern. (2) PA will be monitored for 7 consecutive days. The device will only be removed in the case of swimming (this will be recorded through questionnaires). During this period, adolescents who own a mobile phone will be advised to use the Google Fit App to validate their step count accuracy. The software (ActiLife Software) will be used to analyze the signals. Among other variables, the total PA time will be selected, segmenting it into different degrees of intensity (light, moderate and vigorous).
2.4. Educational intervention
To validate the biomarkers identified in the first stage of the project, an intervention study will be carried out in a subset of the core cohort. All cohort volunteers and their parents will be invited at enrolment to attend a 1 h educational session conducted by two pediatricians at the hospital. Topics covered will include the importance of weight loss in overweight/obesity and its maintenance, a therapeutic nutritional approach to childhood obesity, and the role of PA in cardiovascular fitness. A subgroup of the core cohort encompassing 25% of the participants with normal weight (NW), 25% with overweight (OW), and 50% with obesity (OB), i.e., 60 volunteers (30 for prepubertal and 30 for postpubertal groups) will be formally included in the intervention. Each volunteer participating in the intervention will be given two additional appointments to perform the different proofs/questionnaires and obtain the required biological samples (Table 1). Each individual will be his/her control, so no additional control group is required, and the results will be analyzed individually by performing a pre-post intervention statistical analysis. The number of participants to be included in the educational intervention study has been calculated to achieve a significant 0.2-point BMI z-score reduction between the estimated mean and the sampling mean, with a statistical power of 80% and an alpha risk of 0.05 in a one-sided test. This reduction in BMI z-score has been associated with a clinically significant improvement in children and adolescents with overweight (62). Potential abandonments will be compensated by incorporating more subjects from the core cohort, which is feasible since educational intervention is offered to all subjects from the core cohort (n = 180). A similar sample size has been used with success in previous intervention studies (63, 64).
The educational intervention proposed is based on the PAIDO Programme (www.programapaido.es), an outpatient family-based multidisciplinary programme that combines PA, education on nutrition, and behavioral therapy (63). It is a 6-month intervention aimed to improve PA and dietary habits and thus an improvement of metabolic health and body weight loss is expected, particularly in the OW and OB groups. This intervention has demonstrated efficacy in improving body composition and cardiovascular fitness (63, 64). The dietary intervention will be focused on the promotion of the Mediterranean diet as a healthy nutritional recommendation, with follow-up and advice from a dietitian-nutritionist. Participants will be encouraged to reduce sedentary behavior (watching television, playing computer games, playing board games). The performance of aerobic and strength physical exercises will also be scheduled, progressively increasing the objectives. During the intervention, children, adolescents, and their parents will participate in 6 (online and on-site) sessions (1 per month). At 3 months after the initial evaluation, volunteers of the intervention study will be re-evaluated (body composition) and at 6 months they will be re-sampled (PBC, plasma/serum) and re-evaluated for general parameters (body composition, PA by accelerometry and questionnaires, and KIDMED questionnaire), as above described.
2.5. Laboratory procedures—omic analyses
General metabolic assessment and myokines/adipokines, transcriptome, miRNAs, and metabolome analyses will be performed in plasma and PBC samples from subjects of the core cohort during the biomarker generation phase (phase 1) to provide data for subsequent integrative analysis. In particular:
- Glucose, insulin, lipid profile (total, LDL, and HDL-cholesterol; triglycerides), uric acid and creatinine levels will be evaluated in plasma samples of all subjects of the core cohort using standard automated laboratory enzymatic methods. Insulin resistance will be assessed by the homeostatic model assessment (HOMA) index.
- Plasma levels of circulating myokines and adipokines of special relevance, such as contraction-induced myokines, myokines involved in muscle-adipose tissue crosstalk, and main regulatory adipokines (Irisin, Interleukin-6 and 15, myostatin, FGF21, leptin, adiponectin, etc.) will be measured using commercial assay kits based on the Luminex® xMAP® technology, and ELISA kits in all subjects of the core cohort.
- Whole-genome transcriptomic analysis of PBC samples will be performed by next-generation RNA Sequencing (RNA-Seq) in a representative subset of the core cohort (90 subjects). RT-qPCR analysis of the expression of the most remarkable genes (i.e., those most physiologically plausible and affected) identified in the unbiased transcriptomic study will be used to validate the results in all subjects of the core cohort.
- Combined untargeted and targeted metabolomics and lipidomics analysis of plasma will be performed in all subjects of the core cohort by a LIMeX-4D workflow with 4 different liquid chromatography–mass spectrometry (LC-MS) based methods to annotate a wide range of complex lipids, polar metabolites, and exposome compounds (mainly food components and drugs) (65, 66). Promising biomarkers originating from untargeted metabolomics/lipidomics will be quantified using the reference standardization approach (67).
- For circulating miRNA analysis, RNA will be extracted from plasma samples using the RNAeasy serum/plasma Kit (Qiagen). A miRNome analysis (around 2000 miRNAs) will be assessed in plasma samples of a representative subset of the core cohort (15 of each of the 6 groups), using miRNA sequencing. The most remarkable miRNAs identified (up to 20) will be analysed in the whole cohort by RT-qPCR.
Candidate biomarkers emerging from phase 1 will be measured in plasma and PBC of the subjects in the educational intervention (phase 2) study, at the end of it, using ELISA, RT-qPCR, or biochemical assays, depending on the marker, as above described.
2.6. Data analysis and integration
Data will be expressed as mean ± standard deviation (SD). Data will be checked for normality using Shapiro–Wilks normality test and log-transformation will be applied when required. For multi-group comparisons depending on PA levels, body composition, or other health-related parameters available, Levene's test will be performed to assess whether the variance is equal between groups; if the variance is heterogeneous, data will be log-transformed before analysis. Regression models will be applied to analyse the complex data. Comparisons between two groups will be assessed by the appropriate test (T-test or the non-parametric Mann—Whitney U-test). When appropriate, two-way ANOVA will be used to further determine the effects of different factors, in the case of independent categorical variables. ANCOVA analyses including variables such as BMI, age and gender as covariates will be also used, when appropriate. The threshold of significance will be defined at P < 0.05. Correlation analyses are also envisaged to study the possible relationships between the candidate biomarkers and PA levels and other health-related parameters. SPSS and R software will be used for general statistical analyses.
Omics data from transcriptomics and miRNomics (subset of the core cohort) will be subjected to an exploratory analysis of data through principal component analysis (PCA) and partial least squares discriminant analysis (PLS-DA) (when appropriate). Analysis of differential gene expression will be performed through T-tests, such as Student's or limma, ANOVA, and false discovery rate (FDR) correction. Gene set enrichment analysis and functional and pathway analyses will be performed. Metabolomics datasets will also be processed using PCA and PLS-DA, and metabolite network mapping using KEGG reactant pairs along with Tanimoto similarity.
2.6.1. Integrated analysis
Data from the PA endpoints (PA, muscular fitness, and CRF describing variables) collected for each individual will be mapped in a principal component analysis plot. Hierarchical cluster analysis will be applied to the PCA scores in order to extract 3 clusters corresponding to low, medium, and high PA scores. The significant differences between the group of individuals will be challenged in a supervised partial-least square discriminant analysis (cross-validation ANOVA and permutation test). All the collected biological data, whether clinical or omics, will be examined in each of the 3 groups of individuals defined therefrom (low, medium, and high PA outcomes). Variables selection to fit the class assignment will be performed independently for each kind of variable with a sight of sparsity. PLS-discriminant analysis (orthogonalized or not), sparse-PLS, and machine-learning methods such as random forest will be applied and resilient variables, selected with all the methods, retained. Statistical models will be built on 4/5 of the population, and their performance assessed by predicting the 1/5 unknowns. This will be iterated 10 times. Complementary AUROC analyses will confirm biomarkers selection robustness. The most resilient variables (expected ∼20 variables from any origin) will be saved, and computed in a final PLS-discriminant analysis. From there, an equation calculated from the beta score of the PLS algorithm and the variable values will be established (using the NIPALS algorithm) to give each individual a score value (with 95% confidence intervals) that fits its PA class (e.g., low, medium, high). This equation will be challenged in the subgroup of the core cohort after the educational intervention. Also, the predictive robustness of the individual score value will be challenged to identified confounding factors (BMI, age, pubertal state, and gender, diet and lifestyle-related factors including sleep behaviour) in a logistic regression model.
After biomarkers discovery, the same methodology as above will be used to identify the biological variables that are related to each PA group, but with much less sparsity. This will allow to examine the molecular regulations that occur at the pathway level. Multi-block PLS method will determine the most influencing kind of variables (miR, transcriptomics, lipidomics, metabolomics, clinics) on group assignment. Pathway enrichment analyses will be also used to provide a biological explanation output. Some focus on specific pathways using in silico reconstituted pathway gathering several omics levels could be also performed (ingenuity pathway analysis, MetExplore, and/or Metscape plugin of cytoscape).
3. Discussion
The INTEGRActiv study will address the scientific and societal challenge of identifying integrated markers reflecting both PA level and health (PA-health biomarkers) in children and adolescents, who represent an important target population to address personalized interventions to improve future metabolic health at a societal level. Implementation of the use of the biomarkers relating PA level and health to be identified in the INTEGRActiv study will help the objective assessment of the actual benefits of public health education interventions and PA programs addressed to children and adolescents that are often developed through schools and the primary health care systems. It can also help identify non-responders requiring special attention so that PA programs more adapted to their characteristics can be recommended. Such public programs are urgently required because unhealthy PA habits among children and adolescents are very prevalent in many countries, going hand-in-hand with a high prevalence of overweight and obesity, and there is the awareness that good physical (and nutritional) habits need to be promoted especially in children and youths for their incorporation into their lifestyle, hence ensuring better adherence throughout life. Furthermore, appropriate levels of PA in these critical stages of development might also be required to ensure a proper biological metabolic programming of long-term health.
All in all, the INTEGRActiv study is expected to encompass relevant progress beyond the state-of-the-art for (a) the definition of potential biomarkers for PA level in childhood and adolescence; (b) the provision of relevant mechanistic information for the link between PA and metabolic health in young subjects; (c) the identification of factors such as age, gender, body weight, sleep behaviour and pubertal status that quantitatively affect biomarker responses; (d) the use and development of new tools in biomarker research, including integrative analysis; (e) further assessment and first-step validation of promising candidate biomarkers in an intervention study. It is expected that newly identified robust biomarkers reflecting PA level and its relation with health (PA-health biomarkers) will guide nutritional/lifestyle and clinical advice and public policies related to endorsing biomarker-based personalised PA, with better adherence and response, to promote health and prevent disease risk factors, starting in the early stages of life.
Ethics statement
The studies involving humans were approved by Ethics Committees of the Consorcio Hospital General Universitario de Valencia (registry number 2/2022). The studies will be conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation in this study will be provided by the participants’ legal guardians/next of kin.
Author contributions
CP, EL, JKe, JKo, and J-FL are the principal investigators of the INTEGRActiv project. All authors contributed to the study design. CP, PO, and MB wrote the manuscript. All authors revised and approved the final manuscript. All authors contributed to the article and approved the submitted version.
Funding
The study is supported by the project AC21_2/00033, funded by Instituto de Salud Carlos III (ISCIII) and co-funded by Next Generation EU, Mecanismo para la Recuperación y la Resiliencia (MRR); the Ministry of Education, Youth and Sports of the Czech Republic; and the French National Research Agency in the joint funding activity HDHL-INTIMIC: Standardised measurement, monitoring and/or biomarkers to study food intake, physical activity and health (STAMIFY 2021).
Acknowledgments
We are grateful for the support of the UIB and IdISBa (to CP, PO, AP, MP, JR, AR, JS, FS, and MB), the CIBEROBN (to CP, EL, JA-P, PO, AP, MP, JR, AR, JS, FS, and MB), the Fundación de Investigación, and Hospital General Universitario de Valencia (to EL and JA-P), the Wageningen University (to JKe), the Institute of Physiology of the Czech Academy of Sciences (to JKo and PZ), and to INRAE, INSERM and AMU (to J-FL and JM).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1. Fogelholm M. Physical activity, fitness and fatness: relations to mortality, morbidity and disease risk factors. A systematic review. Obes Rev. (2010) 11(3):202–21. doi: 10.1111/j.1467-789X.2009.00653.x
2. Reiner M, Niermann C, Jekauc D, Woll A. Long-term health benefits of physical activity–a systematic review of longitudinal studies. BMC Public Health. (2013) 13:813. doi: 10.1186/1471-2458-13-813
3. Daskalopoulou C, Stubbs B, Kralj C, Koukounari A, Prince M, Prina AM. Physical activity and healthy ageing: a systematic review and meta-analysis of longitudinal cohort studies. Ageing Res Rev. (2017) 38:6–17. doi: 10.1016/j.arr.2017.06.003
4. Zakari M, Alsahly M, Koch LG, Britton SL, Katwa LC, Lust RM. Are there limitations to exercise benefits in peripheral arterial disease? Front Cardiovasc Med. (2018) 5:173. doi: 10.3389/fcvm.2018.00173
5. Galmes-Panades AM, Bennasar-Veny M, Oliver P, Garcia-Coll N, Chaplin A, Fresneda S, et al. Efficacy of different modalities and frequencies of physical exercise on glucose control in people with prediabetes (GLYCEX randomised trial). Metabolites. (2022) 12(12):1286. doi: 10.3390/metabo12121286
6. Zieker D, Fehrenbach E, Dietzsch J, Fliegner J, Waidmann M, Nieselt K, et al. cDNA microarray analysis reveals novel candidate genes expressed in human peripheral blood following exhaustive exercise. Physiol Genomics. (2005) 23(3):287–94. doi: 10.1152/physiolgenomics.00096.2005
7. Janssen JJE, Lagerwaard B, Nieuwenhuizen AG, Escote X, Canela N, Del Bas JM, et al. Single and joined behaviour of circulating biomarkers and metabolic parameters in high-fit and low-fit healthy females. Int J Mol Sci. (2023) 24(4):4202. doi: 10.3390/ijms24044202
8. Magnussen CG, Smith KJ, Juonala M. When to prevent cardiovascular disease? As early as possible: lessons from prospective cohorts beginning in childhood. Curr Opin Cardiol. (2013) 28(5):561–8. doi: 10.1097/HCO.0b013e32836428f4
9. WHO Guidelines on physical activity and sedentary behaviour. Geneva: World Health Organization 2020. Licence: CC BY-NC-SA 3.0 IGO.
10. Gutin B, Owens S. The influence of physical activity on cardiometabolic biomarkers in youths: a review. Pediatr Exerc Sci. (2011) 23(2):169–85. doi: 10.1123/pes.23.2.169
11. Katzmarzyk PT, Barreira TV, Broyles ST, Champagne CM, Chaput JP, Fogelholm M, et al. Physical activity, sedentary time, and obesity in an international sample of children. Med Sci Sports Exerc. (2015) 47(10):2062–9. doi: 10.1249/MSS.0000000000000649
12. de Mello VD, Kolehmanien M, Schwab U, Pulkkinen L, Uusitupa M. Gene expression of peripheral blood mononuclear cells as a tool in dietary intervention studies: what do we know so far? Mol Nutr Food Res. (2012) 56(7):1160–72. doi: 10.1002/mnfr.201100685
13. Reynes B, Priego T, Cifre M, Oliver P, Palou A. Peripheral blood cells, a transcriptomic tool in nutrigenomic and obesity studies: current state of the art. Compr Rev Food Sci Food Saf. (2018) 17(4):1006–20. doi: 10.1111/1541-4337.12363
14. Caimari A, Oliver P, Rodenburg W, Keijer J, Palou A. Slc27a2 expression in peripheral blood mononuclear cells as a molecular marker for overweight development. Int J Obes (Lond). (2010) 34(5):831–9. doi: 10.1038/ijo.2010.17
15. Oliver P, Reynes B, Caimari A, Palou A. Peripheral blood mononuclear cells: a potential source of homeostatic imbalance markers associated with obesity development. Pflugers Arch. (2013) 465(4):459–68. doi: 10.1007/s00424-013-1246-8
16. Konieczna J, Sanchez J, van Schothorst EM, Torrens JM, Bunschoten A, Palou M, et al. Identification of early transcriptome-based biomarkers related to lipid metabolism in peripheral blood mononuclear cells of rats nutritionally programmed for improved metabolic health. Genes Nutr. (2014) 9(1):366. doi: 10.1007/s12263-013-0366-2
17. Diaz-Rua R, Keijer J, Caimari A, van Schothorst EM, Palou A, Oliver P. Peripheral blood mononuclear cells as a source to detect markers of homeostatic alterations caused by the intake of diets with an unbalanced macronutrient composition. J Nutr Biochem. (2015) 26(4):398–407. doi: 10.1016/j.jnutbio.2014.11.013
18. Sanchez J, Priego T, Pico C, Ahrens W, De Henauw S, Fraterman A, et al. Blood cells as a source of transcriptional biomarkers of childhood obesity and its related metabolic alterations: results of the IDEFICS study. J Clin Endocrinol Metab. (2012) 97(4):E648–52. doi: 10.1210/jc.2011-2209
19. Jang K, Tong T, Lee J, Park T, Lee H. Altered gene expression profiles in peripheral blood mononuclear cells in obese subjects. Obes Facts. (2020) 13(3):375–85. doi: 10.1159/000507817
20. Costa A, Reynes B, Konieczna J, Martin M, Fiol M, Palou A, et al. Use of human PBMC to analyse the impact of obesity on lipid metabolism and metabolic status: a proof-of-concept pilot study. Sci Rep. (2021) 11(1):18329. doi: 10.1038/s41598-021-96981-6
21. Janssen JJE, Lagerwaard B, Porbahaie M, Nieuwenhuizen AG, Savelkoul HFJ, van Neerven RJJ, et al. Extracellular flux analyses reveal differences in mitochondrial PBMC metabolism between high-fit and low-fit females. Am J Physiol Endocrinol Metab. (2022) 322(2):E141–E53. doi: 10.1152/ajpendo.00365.2021
22. Sanchez J, Pico C, Ahrens W, Foraita R, Fraterman A, Moreno LA, et al. Transcriptome analysis in blood cells from children reveals potential early biomarkers of metabolic alterations. Int J Obes (Lond). (2017) 41(10):1481–8. doi: 10.1038/ijo.2017.132
23. Szostaczuk N, van Schothorst EM, Sanchez J, Priego T, Palou M, Bekkenkamp-Grovenstein M, et al. Identification of blood cell transcriptome-based biomarkers in adulthood predictive of increased risk to develop metabolic disorders using early life intervention rat models. FASEB J. (2020) 34(7):9003–17. doi: 10.1096/fj.202000071RR
24. Galmes S, Ruperez AI, Sanchez J, Moreno LA, Foraita R, Hebestreit A, et al. KLB and NOX4 expression levels as potential blood-based transcriptional biomarkers of physical activity in children. Sci Rep. (2023) 13(1):5563. doi: 10.1038/s41598-023-31537-4
25. Kosaka N, Iguchi H, Yoshioka Y, Takeshita F, Matsuki Y, Ochiya T. Secretory mechanisms and intercellular transfer of microRNAs in living cells. J Biol Chem. (2010) 285(23):17442–52. doi: 10.1074/jbc.M110.107821
26. Zhou Q, Shi C, Lv Y, Zhao C, Jiao Z, Wang T. Circulating microRNAs in response to exercise training in healthy adults. Front Genet. (2020) 11:256. doi: 10.3389/fgene.2020.00256
27. Silva GJJ, Bye A, El Azzouzi H, Wisloff U. MicroRNAs as important regulators of exercise adaptation. Prog Cardiovasc Dis. (2017) 60(1):130–51. doi: 10.1016/j.pcad.2017.06.003
28. Wang L, Lv Y, Li G, Xiao J. MicroRNAs in heart and circulation during physical exercise. J Sport Health Sci. (2018) 7(4):433–41. doi: 10.1016/j.jshs.2018.09.008
29. Landrier JF, Derghal A, Mounien L. MicroRNAs in obesity and related metabolic disorders. Cells. (2019) 8(8):859. doi: 10.3390/cells8080859
30. Dimassi S, Karkeni E, Laurant P, Tabka Z, Landrier JF, Riva C. Microparticle miRNAs as biomarkers of vascular function and inflammation response to aerobic exercise in obesity? Obesity (Silver Spring). (2018) 26(10):1584–93. doi: 10.1002/oby.22298
31. Pico C, Serra F, Rodriguez AM, Keijer J, Palou A. Biomarkers of nutrition and health: new tools for new approaches. Nutrients. (2019) 11(5):1092. doi: 10.3390/nu11051092
32. Kelly RS, Kelly MP, Kelly P. Metabolomics, physical activity, exercise and health: a review of the current evidence. Biochim Biophys Acta Mol Basis Dis. (2020) 1866(12):165936. doi: 10.1016/j.bbadis.2020.165936
33. Horakova O, Hansikova J, Bardova K, Gardlo A, Rombaldova M, Kuda O, et al. Plasma acylcarnitines and amino acid levels as an early complex biomarker of propensity to high-fat diet-induced obesity in mice. PLoS One. (2016) 11(5):e0155776. doi: 10.1371/journal.pone.0155776
34. Li VL, He Y, Contrepois K, Liu H, Kim JT, Wiggenhorn AL, et al. An exercise-inducible metabolite that suppresses feeding and obesity. Nature. (2022) 606(7915):785–90. doi: 10.1038/s41586-022-04828-5
35. Karstoft K, Pedersen BK. Skeletal muscle as a gene regulatory endocrine organ. Curr Opin Clin Nutr Metab Care. (2016) 19(4):270–5. doi: 10.1097/MCO.0000000000000283
36. Leal LG, Lopes MA, Batista ML Jr. Physical exercise-induced myokines and muscle-adipose tissue crosstalk: a review of current knowledge and the implications for health and metabolic diseases. Front Physiol. (2018) 9:1307. doi: 10.3389/fphys.2018.01307
37. Li F, Li Y, Duan Y, Hu CA, Tang Y, Yin Y. Myokines and adipokines: involvement in the crosstalk between skeletal muscle and adipose tissue. Cytokine Growth Factor Rev. (2017) 33:73–82. doi: 10.1016/j.cytogfr.2016.10.003
38. de Onis M, Onyango AW, Borghi E, Siyam A, Nishida C, Siekmann J. Development of a WHO growth reference for school-aged children and adolescents. Bull World Health Organ. (2007) 85(9):660–7. doi: 10.2471/BLT.07.043497
39. Costa A, van der Stelt I, Reynes B, Konieczna J, Fiol M, Keijer J, et al. Whole-genome transcriptomics of PBMC to identify biomarkers of early metabolic risk in apparently healthy people with overweight-obesity and in normal-weight subjects. Mol Nutr Food Res. (2023) 67(4):e2200503. doi: 10.1002/mnfr.202200503
40. Morgan JF, Reid F, Lacey JH. The SCOFF questionnaire: a new screening tool for eating disorders. West J Med. (2000) 172(3):164–5. doi: 10.1136/ewjm.172.3.164
41. Banos RM, Cebolla A, Etchemendy E, Felipe S, Rasal P, Botella C. Validation of the Dutch eating behavior questionnaire for children (DEBQ-C) for use with spanish children. Nutr Hosp. (2011) 26(4):890–8. doi: 10.1590/S0212-16112011000400032
42. Cook S, Auinger P, Huang TT. Growth curves for cardio-metabolic risk factors in children and adolescents. J Pediatr. (2009) 155(3):S6.e15–26. doi: 10.1016/j.jpeds.2009.04.051
43. European Food Safety Authority, 2014. Guidance on the EU menu methodology. EFSA Journal. (2014) 12(12):3944. 80 pp.
44. Serra-Majem L, Ribas L, Ngo J, Ortega RM, Garcia A, Perez-Rodrigo C, et al. Food, youth and the Mediterranean diet in Spain. Development of KIDMED, Mediterranean diet quality index in children and adolescents. Public Health Nutr. (2004) 7(7):931–5. doi: 10.1079/PHN2004556
45. Altavilla C, Caballero-Perez P. An update of the KIDMED questionnaire, a Mediterranean diet quality Index in children and adolescents. Public Health Nutr. (2019) 22(14):2543–7. doi: 10.1017/S1368980019001058
46. Mackelvie KJ, McKay HA, Khan KM, Crocker PR. A school-based exercise intervention augments bone mineral accrual in early pubertal girls. J Pediatr. (2001) 139(4):501–8. doi: 10.1067/mpd.2001.118190
47. Kowalski CK, Crocker PR, Kowalski NP. Convergent validity of the physical activity questionnaire for adolescents. Pediatr Exerc Sci. (1997) 9:342–52. doi: 10.1123/pes.9.4.342
48. Bailey DA, McKay HA, Mirwald RL, Crocker PR, Faulkner RA. A six-year longitudinal study of the relationship of physical activity to bone mineral accrual in growing children: the university of Saskatchewan bone mineral accrual study. J Bone Miner Res. (1999) 14(10):1672–9. doi: 10.1359/jbmr.1999.14.10.1672
49. Rey-Lopez JP, Ruiz JR, Ortega FB, Verloigne M, Vicente-Rodriguez G, Gracia-Marco L, et al. Reliability and validity of a screen time-based sedentary behaviour questionnaire for adolescents: the HELENA study. Eur J Public Health. (2012) 22(3):373–7. doi: 10.1093/eurpub/ckr040
50. Chervin RD, Hedger K, Dillon JE, Pituch KJ. Pediatric sleep questionnaire (PSQ): validity and reliability of scales for sleep-disordered breathing, snoring, sleepiness, and behavioral problems. Sleep Med. (2000) 1(1):21–32. doi: 10.1016/S1389-9457(99)00009-X
51. Expert Panel on Integrated Guidelines for Cardiovascular Health and Risk Reduction in Children and Adolescent, National Heart, Lung and Blood Institute. Expert panel on integrated guidelines for cardiovascular health and risk reduction in children and adolescents: summary report. Pediatrics. (2011) 128 (Suppl 5):S213–56. doi: 10.1542/peds.2009-2107C
52. Lurbe E, Agabiti-Rosei E, Cruickshank JK, Dominiczak A, Erdine S, Hirth A, et al. 2016 European society of hypertension guidelines for the management of high blood pressure in children and adolescents. J Hypertens. (2016) 34(10):1887–920. doi: 10.1097/HJH.0000000000001039
53. Ross R, Blair SN, Arena R, Church TS, Despres JP, Franklin BA, et al. Importance of assessing cardiorespiratory fitness in clinical practice: a case for fitness as a clinical vital sign: a scientific statement from the American heart association. Circulation. (2016) 134(24):e653–99. doi: 10.1161/CIR.0000000000000461
54. Lang JJ, Belanger K, Poitras V, Janssen I, Tomkinson GR, Tremblay MS. Systematic review of the relationship between 20 m shuttle run performance and health indicators among children and youth. J Sci Med Sport. (2018) 21(4):383–97. doi: 10.1016/j.jsams.2017.08.002
55. Ortega FB, Ruiz JR, Castillo MJ, Sjostrom M. Physical fitness in childhood and adolescence: a powerful marker of health. Int J Obes (Lond). (2008) 32(1):1–11. doi: 10.1038/sj.ijo.0803774
56. Hogstrom G, Nordstrom A, Nordstrom P. High aerobic fitness in late adolescence is associated with a reduced risk of myocardial infarction later in life: a nationwide cohort study in men. Eur Heart J. (2014) 35(44):3133–40. doi: 10.1093/eurheartj/eht527
57. Santana CCA, Azevedo LB, Cattuzzo MT, Hill JO, Andrade LP, Prado WL. Physical fitness and academic performance in youth: a systematic review. Scand J Med Sci Sports. (2017) 27(6):579–603. doi: 10.1111/sms.12773
58. Lubans D, Richards J, Hillman C, Faulkner G, Beauchamp M, Nilsson M, et al. Physical activity for cognitive and mental health in youth: a systematic review of mechanisms. Pediatrics. (2016) 138(3):e20161642. doi: 10.1542/peds.2016-1642
59. Balke B, Ware RW. An experimental study of physical fitness of air force personnel. U S Armed Forces Med J. (1959) 10(6):675–88.13659732
60. Guixeres J, Redon P, Saiz J, Alvarez J, Torro MI, Cantero L, et al. Cardiovascular fitness in youth: association with obesity and metabolic abnormalities. Nutr Hosp. (2014) 29(6):1290–7. doi: 10.3305/nh.2014.29.6.7383
61. Ruiz JR, Castro-Pinero J, Espana-Romero V, Artero EG, Ortega FB, Cuenca MM, et al. Field-based fitness assessment in young people: the ALPHA health-related fitness test battery for children and adolescents. Br J Sports Med. (2011) 45(6):518–24. doi: 10.1136/bjsm.2010.075341
62. Markert J, Herget S, Petroff D, Gausche R, Grimm A, Kiess W, et al. Telephone-based adiposity prevention for families with overweight children (T.A.F.F.-study): one year outcome of a randomized, controlled trial. Int J Environ Res Public Health. (2014) 11(10):10327–44. doi: 10.3390/ijerph111010327
63. Lison JF, Real-Montes JM, Torro I, Arguisuelas MD, Alvarez-Pitti J, Martinez-Gramage J, et al. Exercise intervention in childhood obesity: a randomized controlled trial comparing hospital-versus home-based groups. Acad Pediatr. (2012) 12(4):319–25. doi: 10.1016/j.acap.2012.03.003
64. Bruño A, Escobar P, Cebolla A, Alvarez-Pitti J, Guixeres J, Lurbe E, et al. Home-exercise childhood obesity intervention: a randomized clinical trial comparing print versus web-based (move it) platforms. J Pediatr Nurs. (2018) 42:e79–84. doi: 10.1016/j.pedn.2018.04.008
65. Paluchova V, Oseeva M, Brezinova M, Cajka T, Bardova K, Adamcova K, et al. Lipokine 5-PAHSA is regulated by adipose triglyceride lipase and primes adipocytes for de novo lipogenesis in mice. Diabetes. (2020) 69(3):300–12. doi: 10.2337/db19-0494
66. Sistilli G, Kalendova V, Cajka T, Irodenko I, Bardova K, Oseeva M, et al. Krill oil supplementation reduces exacerbated hepatic steatosis induced by thermoneutral housing in mice with diet-induced obesity. Nutrients. (2021) 13(2):437. doi: 10.3390/nu13020437
Keywords: integration analysis, cardiorespiratory fitness, metabolomics, transcriptomics, miRNome, adipokines, cytokines
Citation: Picó C, Lurbe E, Keijer J, Kopecky J, Landrier J-F, Álvarez-Pitti J, Martin JC, Oliver P, Palou A, Palou M, Zouhar P, Ribot J, Rodríguez AM, Sánchez J, Serra F and Bonet ML (2023) Study protocol: Identification and validation of integrative biomarkers of physical activity level and health in children and adolescents (INTEGRActiv). Front. Pediatr. 11:1250731. doi: 10.3389/fped.2023.1250731
Received: 30 June 2023; Accepted: 30 August 2023;
Published: 12 September 2023.
Edited by:
Alessandro Leone, University of Milan, ItalyReviewed by:
Tiago Jeronimo Dos Santos, University of Almeria, SpainIrene Rutigliano, IRCCS Casa Sollievo della Sofferenza Hospital, Italy
© 2023 Picó, Lurbe, Keijer, Kopecky, Landrier, Álvarez-Pitti, Martin, Oliver, Palou, Palou, Zouhar, Ribot, Rodríguez, Sánchez, Serra and Bonet. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Paula Oliver paula.oliver@uib.es