Abstract

Background. Information on meal patterns (type, number, relative contribution to energy/day (%), time, and location of meals) is limited or nonexistent. Design. Cross-sectional, nationally representative surveys. Setting. Colombia. Participants. n = 26,115 from 3 to 64 years old. The sample analyzed included 3,127 children between 3 and 4 years old, 13,384 children between 5 and 17 years old, and 9,604 adults between 18 and 64 years old. Data Analysis. Meal patterns were described by age group. Through multiple linear regression, crude and adjusted differences in the categories of the covariates studied were estimated. The number/day of meals was the dependent variable. Results. The number of meals/day (mean ± SD) was 4.4 ± 0.0, without differences by sex , current weight , or wealth index , but there were differences in the level of education of the head of the household and the level of food security of the household . A total of 96.8% of the population eats 3 or more meals/day (95% CI: 96.2, 97.2). The consumption frequency (mean ± SD) of the three main meals was 0.95 ± 0.0 “times/day,” 1.0 ± 0.0 and 0.95 ± 0.0, for breakfast, lunch, and dinner, respectively. Lunch is the meal that makes the greatest relative contribution to the total energy consumed (energy/day), 33.9% (95% CI: 32.7, 35.1). Breakfast is eaten outside the home by 13.0% of the subjects, lunch by 26.0%, and dinner by 3.8%. The minimum fasting interval is 9 hours and the maximum is 10 hours and 30 minutes. The meal pattern is equivalent to type “A,” with three main meals and two or three intermediate meals (midmorning and midafternoon) taken during 15 hours of the day. Conclusions. All age groups had more than four meals/day. The number is directly related to socioeconomic level. Lunch is the main meal.

1. Introduction

The dietary transition as a complementary concept to the nutritional transition is based on the establishment of dietary patterns and on the study of the dietary changes that occur in such patterns [14]. It is important to identify dietary consumption patterns (Traditional, Western, Prudent, etc.), as well as the type and number of the meals that compose them (meal patterns or eating patterns) [57]. It is possible to establish a pattern of meals and associate it with cultural, political, social, economic, and biological changes [5, 811]. Meal patterns are also associated with the time in which they are made (meal timing) and with the quality of the diet [5, 1214]. A key aspect is to establish the definition of what the term meal means. There are many definitions, but four are the most relevant: the one based on the time of the day when it is performed, the one based on the quantity and quality of the foods that compose it, and the one, that is, self-reported by the subjects and/or based on a pre elaborated list and another one called neutral, which presupposes a post classification by convenience by the researcher in order to establish comparisons [5, 15]. The used definition influences the results by relating the pattern of meals with the consumption of nutrients, the quality of the diet, the frequency of consumption, or even health-disease events [5]. The main sources of data to establish a meal pattern are the reminders of the dietary consumption of the last 24 hours (24-hour dietary recall, (24 HR)) and the food frequency questionnaires (FFQ) [5, 16].

In Colombia, based on the two most recent National Surveys of the Nutritional Situation (ENSIN), which were conducted in 2010 and 2015, three dietary patterns were identified and the stages of alimentary and nutritional transitions were established [2, 17]. In addition, adherence to dietary consumption patterns by geographic area [4], consumption of macro and micronutrients, their dietary sources, and the quality of the diet throughout the life cycle was characterized, all associated with biological and socioeconomic characteristics [18, 19]. In Colombia, as in many countries, the formal study of meal patterns is limited. The study of meal patterns is a developing field, and the type of food, its frequency of consumption [20], the number of meals/day, and the time in which they are made have been found to be associated with circadian rhythm [12, 13], metabolic changes [21], body weight and its control [22, 23], mediating events of chronic disease such as high levels of cholesterol and hypertension, and the development of diabetes and cardiovascular disease [1234]. Meal patterns are also associated with glycemic alterations in pregnant women and preterm births [35].

Meal patterns can function as a potential cause and consequence. As a cause, they are linked to biological and health-disease events. Complementarily, they are the consequence of complex interactions between culture, economic and political crisis, climate change, food marketing, employment, income, migrations, educational level, and food insecurity; additionally, in Colombia, they are the consequence of internal conflict, violence and forced internal displacement of populations that are associated with territorial control by violent groups to develop illegal economies [3638].

Based on information collected in 2015 in Colombia by the National Nutritional Situation Survey (ENSIN-2015) [39, 40], the objectives of this study were: (a) to estimate the number of meals/day associated with some biological and socioeconomic variables of interest, (b) to estimate crude and adjusted differences in the number of meals in the categories of these variables, (c) to estimate the absolute kilocalories and relative (%) contributions of each type of meal to the total energy intake/day (kilocalories), (d) to establish the prevalence of meal intake (%), the frequency of consumption (times/day) and the place where the different types of meal are made, and (e) to estimate and describe the time at which they are taken.

2. Methods

Colombia is a developing country with medium income and widespread social inequalities, that is, located in the northeastern corner of South America [41, 42]. In Colombia during 2015, under the leadership of the Colombian Institute of Family Welfare (ICBF) and the Ministry of Health (MinSalud), ENSIN-2015 was conducted. In ENSIN-2015, 44,202 households were surveyed, representing 4,739 groups of 295 strata representing 99% of the country’s population. The methods, the populations studied, and the scope and limitations of ENSIN-2015 have already been published [39]. Minor children between 3 and 17 years old and adults between 18 and 64 years old, excluding pregnant women and girls, were the target population of this analysis. The study population answered, among others, a FFQ, a 24 HR, and a sociodemographic survey. The ENSIN-2015 included 151,343 people, of which 28,902 answered questions about the type of food consumed (yes/no), the frequency, and the place, where they were made. After excluding pregnant women (n = 1,939) and subjects outside the age range or with incomplete information (n = 848), the final sample comprised 26,115 subjects including 3,127 children between 3 and 4 years of age, 13,384 children between 5 and 17 years of age, and 9,604 adults between 18 and 64 years of age.

2.1. Data Sources

Based on the FFQ, we established the number of meals/day, the prevalence (%), the frequency of consumption (times/day) of each meal, and the place, where it was consumed (not where they were prepared). Based on the 24 HR, we established the time of consumption of each type of meal (Hour:Minutes), the usual intake/day (kilocalories), and the relative contribution (%) of each meal to the total intake.

The main outcome variables were: (a) the number of meals per day (meals/day), (b) the prevalence (%) of each type of meal, (c) the frequency of each meal (times/day), (d) the predominant place (%) where they are taken, and (e) the time of day at which they are taken (Median, P25, and P75). Based on the FFQ, the number of meals was estimated by investigating whether eight specific meals from a preestablished list were usually made (before breakfast, breakfast, mid-morning, lunch, mid-afternoon, dinner, after dinner, and another) and the prevalence of realization (%) through a dichotomous response (yes/no) on the usual intake in the last week. The frequency of each type of meal was established based on five categories referring to the past week (every day, between five and six days, between three and four days, between two and three days, and between one and two days). After translating the previous response options to a continuous time variable, the frequency of each type of meal was expressed as times per day (times/day). In addition, the FFQ inquired about the predominant location of each type of meal, also referring to the last week. For this variable, there were ten response options (home, school/university, family or friend’s house, restaurant, cafeteria, street-mobile stand, institutional program, work, transportation, and others). Based on 24 HR, the average energy (kilocalories) provided by each type of meal and the relative contribution (%) that each of them makes to the total energy consumed per day (kilocalories/day) were estimated. ENSIN-2015 applied a repeated 24 HR in subsample [39], following the methodology developed in 1999 by the United States Department of Agriculture (USDA) [43, 44]. The 24 HR had a response rate of 92%. The distribution of the intake of kilocalories/day was normalized and corrected by incorporating the intraindividual variability with the methods proposed by the University of Iowa and the PC-Side software, v1.0 [45]; therefore, the reported intake (kilocalories) and the relative contribution (%) correspond to the usual intake of the subjects. In summary, the usual intake is a value derived from the two 24 HR (but not as a simple average), to represent the long-term intake, rather than the intake of the last 24 hours [45].

In addition, based on 24 HR, the time at which each type of meal was eaten was estimated (Median, P25, and P75).

Eight biological and sociodemographic variables were the covariates of interest; these covariates, along with the FFQ and the 24-hour-recall (24 HR), were assessed by direct interviews by dietitians, nutritionists, and trained personnel. Household food insecurity was established with the Latin American and Caribbean Scale (ELCSA) adopted by the Food and Agriculture Organization of the United Nations (FAO) [46]. The ELCSA classifies households as safe and unsafe and the level of food insecurity as mild, moderate, or severe [46]. The wealth index is a continuous indicator of the socioeconomic level of the household, and it was categorized into quartiles. It was determined in the total sample of the survey (ENSIN-2015). It was established based on the analysis of the main components of a set of physical characteristics of the household, the goods in the household, and the availability of services, all of which were previously established in international surveys of demography and health [39, 47]. The level of education of the head of the household was established based on the approved years of study. The geographical region is a variable that represents the territory and the conditions of structural, economic, and cultural development of the subjects. Colombia has five geographical regions [39]. Bogotá is the capital of the country and, together with the central region, has the largest human development index, while the Pacific region along with the Amazonia-Orinoquia has the worst human development index [42]. The most prevalent population in the Pacific region is Afro-descendant individuals, followed by indigenous and mestizo individuals. The level of urbanism was established based on the concentration of the population in three categories: urban centers have more than one million inhabitants and large cities; small towns have between one hundred thousand and one million inhabitants; dispersed populations have fewer than one hundred thousand inhabitants [39].

Anthropometric measurements were obtained from all household members with the use of calibrated standards, techniques, and instruments. In ENSIN-2015, height was measured with the use of stadiometers (ShorrBoard), with a sensitivity of 1 mm. Weight was measured with SECA 874 scales, with a precision of 100 g [39]. In minors, Z-scores were established for the height/age and Body Mass Index (BMI) indicators according to the growth standards of the World Health Organization (WHO) [48].

2.2. Data Analysis

All analyses were conducted under the consideration of the complex sampling design and using Stata v14.1 [49] to, (a) estimate and describe the average number of meals/day by the categories of the covariates of interest, (b) estimate the crude and adjusted differences in the number of meals/day between categories of the covariates studied and by age group, (c) describe the prevalence (%) of the type of meal made by age, the frequency/day (times/day), the usual energy provided by each meal (kilocalories) and the relative contribution of each meal (%) to the total energy consumed in the day (kilocalories/day), (d) establish the prevalence (%) of the place where each type of meal is made by age group, and (e) establish the time (Median, P25, P50) at which each meal was taken. The variables were represented with means or proportions ± standard error (SE) and 95% confidence intervals (95% CI). The adjusted differences were obtained through multiple linear regression models with number of meals as a continuous result and biological and sociodemographic variables as predictors except the height-for-age or height- and BMI-for-age or BMI. The estimates for education come from a model that excludes the wealth index and food security, which could be on the causal path. The wealth index estimates excluded food security.

The study was conducted according to the guidelines laid down in the Declaration of Helsinki [50]. Consent for participation in the survey was obtained by the Colombian Institute of Family Welfare prior to enrollment [39]. The ethics committee in health research of the Universidad Industrial de Santander determined that analyses of these anonymized data were exempt from review.

3. Results

A total of 49.4% of the participants were men. A total of 8.2% of the children presented stunting. A total of 13.6% (95% CI: 11.0, 16.8) of children between 3 and 4 years of age, 15% (95% CI: 13.2, 16.9) of children between 5 and 17 years of age, and 14.5% (95% CI: 12.8, 16.4) of adults were obese (Z-score >2 or BMI ≥30). A total of 24.7% (95% CI: 22.4, 27.2) of adults had primary education or less, and 9% (95% CI: 6.9, 11.6) had higher education or more. 38.0% (95% CI: 29.5, 49.3) of the participants were in Q1 of the wealth index, and 13.7% (95% CI: 10.8, 17.1) were in Q4. 36% (95% CI: 26.6, 46.5) live in the four main cities of the country. The number of meals/day (mean ± SD) in the population aged 3 to 64 years was 4.4 ± 0.0; and there were no differences based on sex, value = 0.068, current weight, value = 0.336 or wealth index value = 0.480, but there were differences based on the level of education of the head of the household, value < 0.0001 and the household food security, value < 0.0001. The number of meals/day in households with heads with less than primary school was 4.3 ± 0.0, 0.4 meals/day less than in households, where the head has a complete university education, 4.7 ± 0.1, for trend -value<0.0001.

In Tables 14 show the mean number of meals/day and the crude and adjusted differences between the categories of the socioeconomic variables studied. There are significant differences in favor of men, and the relationship between the number of meals/day with the level of education, height-for-age (Z-score), and household food security is direct. In addition to the above, an inverse relationship between the number of meals/day and body mass index is observed in adults. Figure 9(s) shows the mean number of meals/day by age.

The point prevalence of having one or fewer meals/day was 0.00% (95% CI: 0.00, 0.00), that of having two meals was 0.03% (95% CI: 0.03, 0.04), that of having three meals/day was 0.23% (95% CI: 0.22, 0.25), that of having four meals/day was 0.27% (95% CI: 0.25, 0.28), that of having five meals/day was 0.28% (95% CI: 0.27, 0.29), that of having six meals/day was 0.14% (95% CI: 0.13, 0.15), and that of having seven meals was 0.04% (95% CI: 0.04, 0.05). Finally, the prevalence of eating eight meals/day was 0.00% (95% CI: 0.00, 0.00). A total of 96.8% of the population had 3 or more meals/day (95% CI: 96.2, 97.2). Figure 2(s) shows the probability of having a specific type of meal/day according to the level of food security of the household.

The prevalence of before breakfast was 25.0% (95% CI: 22.8, 27.3), that of breakfast was 95.6% (95% CI: 94.9, 96.2), that of the midmorning was 46.8% (95% CI: 44.7, 48.9), the lunch was 98.7% (95% CI: 98.1, 99.0), the midafternoon was 51.9% (95% CI: 50.0, 53.7), the dinner was 94.8% (95% CI: 93.7, 95.6), and finally, the prevalence of after dinner was 25.3% (95% CI: 23.2, 27.5). The frequency of consumption, times/day (mean ± SD) of the three main meals was 0.95 ± 0.0, 1.0 ± 0.0, and 0.95 ± 0.0 for breakfast, lunch, and dinner, respectively. Lunch was the meal that made the greatest relative contribution to the total energy consumed (energy/day), 33.9% (95% CI: 32.7, 35.1), followed by dinner, with a relative contribution of 24.7% (95% CI: 23.8, 25.6), and breakfast, with a relative contribution of 22.7% (95% CI: 21.4, 23.9). Table 5 presents the prevalence and frequency (times/day) by type of meal, in addition to the absolute kilocalories and relative (%) contribution to the total energy/day by the age groups studied. The relative contribution (%) that lunch makes to the total energy/day is directly related to the level of education of the head of the household, value = 0.012. The relationship between the relative contribution (%) and the wealth index is direct with the midmorning, value = 0.001, lunch, , the midafternoon, and with after-dinner value 0.020. Figure 3(s) shows the relative contribution (%) of each type of food to the total energy/day according to age.

All meals were mostly eaten at home; breakfast was done outside the home by 13.0% of the subjects, lunch by 26.0%, and dinner by 3.8%. Table 6 presents in detail the place where each type of meal is made by the age groups studied. Finally, the time (Hour:Minutes) at which each meal is done (Median, P25, P75) is 6:30 (6:00, 7:20) before breakfast, 8:00 (7:00, 9:00) for breakfast, 10:00 (9:30, 11:00) for midmorning, 12:30 (12:00, 13:00) for lunch, 16:00 (15:00, 17:00) for midafternoon, 19:00 (18:00, 19:30) for dinner and 20:20 (19:30, 21:00) for after dinner. Figure 4(s) shows, in addition to the time of completion, the relative density by the number of food records at each meal time. Lunch is the most concentrated meal regarding the time in which it is eaten, and breakfast and midafternoon are the most dispersed at the time of eating.

4. Discussion

The findings of this study established the pattern of meals in the Colombian population between 3 and 64 years of age. In 2015, all age groups had more than four meals/day. The number is directly related to variables that represent socioeconomic level (wealth index) or proximal to it (household food security and level of education). The home continues to be the main place, where all meals are made. Lunch and midafternoon are the two meals that are most often eaten outside the home. The lower prevalence of midmorning and midafternoon is what differentiates insecure households from those with food security. Contrary to popular belief in this population, breakfast is not the main meal in terms of absolute kilocalories or relative (%) contribution of energy/day. Lunch is the main meal, with an advantage of nine percentage points in the relative contribution over dinner and breakfast. According to the time in which the meals are made, lunch is the most concentrated, and breakfast and midafternoon are the most dispersed during the day. The meal pattern reported in this study is classified as type A [13]. Since there is no reason to suppose that this pattern changed in the last decade in Colombia, it is reasonable to state that it is a potential risk factor for the development of chronic diseases. The “A” meal pattern is a common meal pattern in which the obesity epidemic developed and in which diabetes and associated chronic diseases have emerged [13].

As already noted, the definition and classification of the term “meal” influence subsequent results [5], which makes it difficult to compare our findings. A study with noninstitutionalized adults in the USA (>45 y), where the three main meals were investigated and the rest were called “snacks” showed that at least 67% ate three or more meals/day (herein, the prevalence was 96.8%), and dinner had the highest prevalence (36%) [6]. A study conducted among U.S Adventists (≥30 y) that measured meals with six categories found that they eat 4 ± 1 meals/day (mean ± SD) and that lunch is the meal with the least variability in the time of day it is eaten [30]. A study of Puerto Rican children aged 10–17 years, in which eight types of meals were established, showed that they consumed between 4.3 ± 1.4 and 5.0 ± 1.1 meals/day (mean ± SD). Younger children had more meals/day and those with inadequate weight less meals/day [16]. The report of fewer meals with a higher body mass index (BMI) in Puerto Rican children, as also reported in this study, is likely to be a reporting bias, which is well known when studying dietary intake [51]. A study in Swedish women (>45 y), where the number of meals was established based on the time of day at which they were made, showed that they consumed between 6.1 vs. 5.2 meals/day, with an advantage for obese women compared to those who had adequate weight. In addition, it was established that obese women ate more meals/day in the evening (after 16:00 hours) than women with adequate weight [24]. In adults residing in the USA (20–70 y), based on eight preestablished meal times, it was estimated that they ate an average of 3.9 ± 0.8 meals/day. In addition, the percentage of meals outside the home was 29.7%; 18.9% for breakfast, 53.5% for lunch, and 19.6% for dinner [51]. A study conducted with an adult population in Europe (45–65 y), which preestablished 11 meal times, found that the number of meals/day was (mean ± SD) 4.3 ± 0.8 and 4.5 ± 1.0 for men and women in France, respectively, 5.2 ± 1.1 and 5.0 ± 0.9 for men and women in Norway, 5.7 ± 1.3 and 5.4 ± 1.0 for men and women in Belgium, 5.5 ± 31.3 and 5.7 ± 1.0 for men and women in the Czech Republic and 7.0 ± 1.6 and 7.1 ± 1.3 for men and women in the Netherlands. In Belgium and the Netherlands, subjects with a lower level of education were found to eat more meals/day than those with a higher educational level, and as reported here, the prevalence of the three main meals ranges from 79.8% (dinner in Norway) to 98.1% (breakfast in the Czech Republic) [20].

Comparisons with previous studies are difficult not only due to the different definitions and classifications that are made of the number of meals but also due to the methods used in the estimation of the prevalence, the frequency/day, and in the analytical approach, where the biological approaches are more common than socioeconomic approaches. In Colombia, the inability of households to physically access food measured through the ELCSA [37, 52] is one of the main determinants of alimentary and nutritional transitions [1720, 2529]. However, fewer meals/day, especially midmorning and midafternoon meals, do not necessarily translate into a lower quality diet if, as has been evidenced in the Colombian case, in population terms, there are no deficiencies in the consumption of energy and proteins [19, 28]. There are often such deficiencies among vulnerable groups (indigenous, displaced population, extremely poor, migrants, etc.), but this issue is beyond the scope of the current study.

Two or three meals/day, in particular the main ones and especially breakfast and in addition, fasting between 12 and 16 hours, are an ideal reported in some studies to maintain body weight, increase insulin sensitivity, decrease serum cholesterol and the sensation of hunger, improve circadian rhythm, decrease inflammation, increase resistance to stress, and improve intestinal microbiota and other risk markers for the development of cardiovascular and chronic disease in general [12]. However, the above statement is only a hypothesis that deserves further study and should be carefully observed. Based on the hours in which each of the meals that make up the pattern of meals in the Colombian population are made, the minimum fasting interval is 9 hours, and the maximum is 10 hours and 30 minutes. We established that the pattern of meals is equivalent to the pattern of type “A,” with three main meals and two or three intermediate snacks, all taken at an interval of 15 hours of the day [13]. In addition, two of them–dinner and a smaller snack–are made at night (Figure 4(s)). In Colombia, ENSIN data show that the average value of the BMI continues to increase. This shift to the right of the weight distribution has occurred under the meal pattern reported here (Type A) [13, 17].

Based on the prevalence and the place, where meals are made, it can be hypothesized that Colombia is far from being an industrialized society, where basic preindustrial and cultural relationships with food prevail over economic ones. Making the main meals mostly at home shows both a high level of economic dependence as well as limited participation in the formal production apparatus and the industry. The unemployment rate in 2015 was 8.9%, and the informal employment rate was 48.5%, thus supporting this hypothesis [53]. While in France, 53.5% of the subjects eat lunch outside the home [51], only 26.0% do so in Colombia.

4.1. Scope and Limitations of the Study

The main strength of this study is that the data come from a national survey, with high-quality data obtained through the 24 HR and FFQ methods. However, given the cross-sectional nature of the data, it is not possible to establish causal relationships. Another limitation is that, given the intentionality of the analysis, the number of meals/day, their prevalence or their frequency/day, or even the usual absolute intake in each type of meal (kilocalories) or their relative contribution (%) to the total energy intake (kilocalories/day) cannot be equated with the quality of the diet. This study is another piece of a complex puzzle to try to characterize different aspects of what has been conceived as the alimentary transition in Colombia. The methods used, such as the questionnaires applied in ENSIN-2015, can be used to search for causal relationships with biological variables and, in addition, to establish baselines in comparisons with future ENSINs and with studies that aim to establish the impact of food insecurity, wealth index, poverty, employment rate, income, economic or political crises, etc., and other socioeconomic variables on the number of meals/day. It is well known that the coronavirus SARS-CoV-2 pandemic and its social, economic, and political consequences have negatively impacted, first acutely and now chronically, all the indicators that potentially modify meal patterns [54, 55].

Finally, the results can be incorporated into public policy to overcome the vision focused on nutrients or foods and complement it with eating patterns and meals, closer to the daily life of the subjects, and that allow improvements in consumption habits and related practices.

Data Availability

To access the ENSIN 2015 public database, you must register in the repository of the Ministry of Public Health: [email protected] and make the request through the format available at: https://www.minsalud.gov.co/sites/rid/paginas/freesearchresultsf.aspx?k=Base%20de%20datos%20Encuesta%20Nacional%20de%20la%20Situaci%C3%B3n%20Nutricional%20ENSIN%202015.

Additional Points

The lead author affirms that this manuscript is an honest, accurate, and transparent account of the study being reported. The reporting of this work is compliant with STROBE guidelines. The lead author affirms that no important aspects of the study have been omitted and that any discrepancies from the study as planned have been explained.

Disclosure

Funding sources did not play any role in the design, conduct, or interpretation of the study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

The authors’ contributions were as follows: M. P. Z. and O. F. H. designed the research. O. F. H led the statistical analysis. M. P. Z. and O. F. H. wrote the paper. O. F. H. has primary responsibility for the final content. All authors have read and approved the final version of the manuscript.

Supplementary Materials

Figure 1(s): Number of meals per day (mean and 95% CI) according to age in the Colombian population (3–64 y, non-pregnant women). National Survey of Nutritional Situation in Colombia (ENSIN-2015). Figure 2(s): Probability of meal type per day according to food security at the home in the Colombian population (3–64 y, non-pregnant women). National Survey of Nutritional Situation in Colombia (ENSIN-2015). Figure 3(s): Relative contribution (%) to total kilocalories/day according to the type of meal in the Colombian population (3–64 y, non-pregnant women). National Survey of Nutritional Situation in Colombia (ENSIN-2015). Figure 4(s): Time of day in which meal times are made (Median, P25, P75) and density by the number of foods in them in the Colombian population (3–64 y, non-pregnant women). National Survey of Nutritional Situation in Colombia (ENSIN-2015). Appendix. Appendix 1-(a) to Appendix 1-(g) present the contributing food categories for each type of meal (proportion, %). Due to the volume of presenting this information in detail by subcategories that are part of each category or foods that make up each subcategory, this can be requested from the authors with the intention of collaborative work. There are four supplementary figures and one appendix. (Supplementary Materials)