Malnutrition is the top risk factor for death and disability in Ethiopia, responsible for 24 % of all deaths(1). The country faces a high double burden of malnutrition, with 22 % of women and 33 % of men underweight, and 8 % of women and 3 % of men overweight or obese. In Ethiopia’s capital city, Addis Ababa, the prevalence of overweight and obesity (29 % of women and 20 % of men) has overtaken that of underweight (13 % of women and 18 % of men)(2). Overweight and obesity, along with unhealthy lifestyles, are important risk factors for non-communicable diseases such as CVD, type 2 diabetes and some cancers(3). Obesity is particularly harmful in low-income settings with high rates of childhood malnutrition, as exposure to undernutrition in early life exacerbates the relationship between obesity and risk of non-communicable diseases(Reference Wells, Sawaya and Wibaek4). As in other urban areas in sub-Saharan Africa (SSA), non-communicable diseases such as diabetes and CVD have become major causes of morbidity and mortality in Addis Ababa, responsible for 31 % of deaths reported by hospitals(Reference Tebekaw, Teller and Colón-Ramos5).
Innovative strategies are needed to address the double burden of malnutrition, especially in low-income, urban settings. Increasing the proportion of energy intake from protein could be a potential intervention to simultaneously prevent or treat underweight and overweight. Protein deficiency, contributing to protein-energy malnutrition, is a major cause of disability in Ethiopia(Reference Misganaw, Melaku and Tessema6). In 2013, an estimated 49 % of women and 69 % of men in Ethiopia were below the estimated average requirement for protein intake(7). Protein intake in low-income countries such as Ethiopia tends to be limited, as the majority of energy intake comes from cereal-based staple foods, with on average only 3 % of energy from meat, 11 % from roots and tubers and 6 % from pulses, nuts and oilseeds(Reference Schönfeldt and Gibson Hall8). Low protein intake can lead to stunting and wasting, especially in childhood(Reference Schönfeldt and Gibson Hall8). An analysis of 180 countries found that higher national-level estimates of total protein intake were associated with lower prevalence of child stunting(Reference Ghosh, Suri and Uauy9). Higher utilisable protein estimates (which take into account the digestibility of the protein sources and the essential amino acid composition) were associated with lower prevalence of child stunting, independent of energy intake(Reference Ghosh, Suri and Uauy9). One study estimated a potential 1·19-year increase in life expectancy at birth if the adult population of Addis Ababa ate enough protein to meet the required daily amounts and a 0·42–2·0 percentage point decrease in child stunting, with animal-source foods more efficacious than plant-source foods in reducing stunting(Reference Blakstad, Danaei and Tadesse10).
While evidence from low-income settings is limited, studies in Western populations suggest that increased protein intake may have favourable effects on body weight and cardiometabolic health. In the USA and Canada, protein intake above the RDA has been associated with lower body weight, more favourable body composition and lower waist circumference(Reference Green, Shea and Vasdev11,Reference Pasiakos, Lieberman and Fulgoni12) . A meta-analysis of short-term randomised trials suggests that high protein diets (20–35 % of energy from protein) may lower cardiometabolic risk through changes in body composition and/or weight(Reference Wycherley, Moran and Clifton13). Higher-protein diets often have a relatively low energy density, which can aid appetite suppression, reduce overall food intake, preserve lean body mass and promote weight management(Reference Liu, Silvestre and Poppitt14).
While evidence from short-term studies indicates that overall protein intake may be beneficial for improving cardiometabolic health, it is important to consider potential divergent health effects of distinct protein sources. Traditionally, animal-source protein has been thought to be beneficial for the reduction of undernutrition(Reference Willett, Rockström and Loken15,Reference Krasevec, An and Kumapley16) , but the relationship between animal protein, plant protein and cardiometabolic health is unclear. In the USA, diets higher in both plant and animal protein, independent of other dietary factors, have been associated with cardiometabolic benefits including decreased BMI and waist circumference(Reference Berryman, Agarwal and Lieberman17). In Korea, increased amounts of both plant and animal protein have been found to decrease BMI and waist circumference(Reference Park, Park and Kang18). However, in Belgium, plant protein intake was associated with lower BMI and waist circumference in males and females, while animal protein was associated with higher BMI and higher waist circumference in males(Reference Lin, Bolca and Vandevijvere19). Plant protein intake has also been associated with favourable changes in waist circumference in the USA(Reference Hruby and Jacques20) and lower risk of metabolic syndrome in Australia(Reference Shang, Scott and Hodge21). In two large US cohorts, total protein intake was not associated with total mortality, but animal protein was associated with higher mortality and plant protein was associated with lower mortality. However, these findings may reflect in part the types of fat in animal and plant foods and other unhealthy components of animal-source foods commonly consumed in Western settings (such as sodium, nitrates and nitrites found in processed red meat) rather than just the type of protein(Reference Song, Fung and Hu22).
Nearly all evidence on protein intake and cardiometabolic health is from high-income countries, which may not be generalisable to low-income settings where animal protein consumption is much lower and food sources of protein differ. While shifts from animal to plant sources of protein are encouraged in Western populations for health and environmental benefits, promotion of animal-source foods in place of carbohydrate in low-resource settings in SSA, where the majority of energy is from starchy carbohydrates, may improve dietary quality, micronutrient intake and have metabolic benefits through reducing glycaemic load(Reference Willett, Rockström and Loken15,23) . Ethiopia is currently creating dietary guidelines to address the double burden of malnutrition(Reference Bekele, de Vries and Trijsburg24); however, there is a lack of robust dietary data on protein consumption in Ethiopia and a need to clarify how distinct protein sources are related to underweight and overweight in order to form evidence-based guidelines. Therefore, we aimed to examine the associations between proportion of energy from total, animal and plant protein, as well as food sources of protein, with underweight and overweight/obesity in Ethiopian adults.
Methods
Household survey
Between January and March 2018, we conducted a cross-sectional study among 1050 urban households in Addis Ababa including 1050 adult women and 635 adult men 18–49 years of age to characterise patterns in protein-source food production, access and consumption. We used multi-stage random sampling to select households from five sub-cities of Addis Ababa, Ethiopia. Eligible households included a woman of reproductive age (18–49) and at least one child between 6 and 59 months of age. If there was an adult male age 18–49 years available at the household, they were also included in the sample. If more than one adult male or adult female was living in the household and present at the time of the interview, one was randomly selected. If the woman of reproductive age was pregnant, she was included in the survey; however, pregnant women were excluded from this analysis (n 56). We invited 1083 eligible households to participate and 1050 provided informed consent (response rate of 97 %). The woman of reproductive age completed a household survey including questions on socio-economic status and demographic characteristics. Both the adult man and woman completed a seventy-three-item non-quantitative FFQ assessing the number of days each food was consumed out of the past 7 d, which was administered at the participant’s household by a trained interviewer. Portion sizes and frequencies of intake were not collected. This FFQ was adapted from a semi-quantitative FFQ previously validated for use against two 24 h diet recalls among urban adults in Dar es Salaam, Tanzania(Reference Zack, Irema and Kazonda25) and was modified to reflect foods commonly consumed in Ethiopia. Height was measured for all men and women to the nearest 0·1 cm with the subject barefoot using a stadiometer. Weight was measured to the nearest 0·1 kg with the subject barefoot and in light clothes using a SECA digital scale. BMI was calculated as weight in kilograms divided by height in metres squared. Based on BMI, individuals were classified using standard cut-offs as underweight (< 18·5 kg/m2), healthy weight (18·5–24·9 kg/m2), overweight (25–29·9 kg/m2) or obese (≥ 30 kg/m2). All participants had complete height, weight and dietary data.
Assessment of dietary intake
We calculated energy, protein, carbohydrate and fat intakes using (1) FFQ dietary intake data, (2) nutrient contents from the Ethiopian Food Composition Table(26) and (3) portion sizes for men and women obtained from previous nutrition surveys in Addis Ababa(Reference Tadesse, Hemler and Andersen27). We used the Tanzanian Food Composition Table(Reference Lukmanji, Hertzmark and Mlingi28) where Ethiopian estimates were missing. No participants had implausibly high energy intakes (greater than 4500 kcal/18 828 kJ per day), but five participants with energy intakes less than 500 kcal/2092 kJ per day were excluded from the analysis. We estimated the prevalence of inadequate protein intake by applying the estimated average requirements for protein intake for each individual (0·66 g/kg/d for adult men and non-lactating women and 1·05 g/kg/d for lactating women)(29).
Covariate Measurements
Covariates were collected in the socio-demographic household survey administered to women and included age, highest level of school completed, religion, marital status, employment, household size, indicators of household living standards and household food insecurity. The women’s covariates were used for the men living in their households. All households had complete covariate data. To measure household wealth, participants were categorised into quartiles based on a wealth index that was constructed using principal components analysis and a set of twelve indicators describing household living standards and possessions owned(Reference Vyas and Kumaranayake30). Household food insecurity was assessed using the Household Food Insecurity Access Scale, which is a standardised scale ranging from 0 to 27 based on a nine-item questionnaire including the following food security domains: uncertainty about the household food supply, insufficient quality of food and insufficient food intake. This scale has been used across cultural contexts and categorises households and populations as food secure, mildly food insecure, moderately food insecure or severely food insecure(Reference Coates, Swindale and Bilinsky31).
Statistical analysis
Descriptive statistics were presented as means and standard deviations; medians and 25th/75th percentiles or proportions. Relationships between socio-demographic characteristics and protein intake were assessed by testing for linear trends in socio-demographic characteristics across quartiles of percentage energy from total protein intake using linear regression (for continuous variables), logistic regression (for binary variables) and ordered multinomial logistic regression (for categorical variables with more than two levels). Linear models adjusting for socio-demographic variables were used to examine associations between percentage energy from total protein, plant protein and animal protein and continuous BMI. Percentage energy from total protein, plant protein and animal protein were examined in linear models as continuous exposures and also in quartiles, with the first quartile representing the lowest percentage energy intake from total, plant or animal protein. We assessed potential nonlinearity in the associations between BMI and percentage energy from total, plant and animal protein by modelling each exposure in multivariate-adjusted models as cubic splines with 2 df.
Logistic regression models were used to examine the associations of percentage energy from total, plant and animal protein (both in quartiles and as continuous exposures) with a binary variable for overweight/obese (defined as BMI ≥ 25 kg/m2) and a binary variable for underweight (BMI < 18·5 kg/m2) to examine effects of protein sources on underweight and overweight/obesity independently. Overweight/obese and underweight models were run separately excluding underweight participants from models with an outcome of overweight/obesity and vice versa to maintain a comparison group of those with a healthy weight. Analyses were also performed separately by sex because women and men may have different risks of underweight and overweight(Reference Hoque, Hasan and Rahman32,Reference Finucane, Stevens and Cowan33) . In addition, the top ten food sources of protein intake in the sample were identified and logistic regression models were used to investigate the associations between percentage energy from each protein food with overweight and underweight.
All dietary factors were expressed as nutrient densities (percentage of total energy intake) to examine diet composition without constraining total energy intake, as this may be the primary mediator of how dietary exposures affect outcomes of body weight(Reference Willett34,Reference Mozaffarian, Hao and Rimm35) . We used these nutrient density models to estimate the effect of substituting percentage energy from total, animal and plant protein for an equal percentage energy from carbohydrate as well as substituting animal for plant protein. In all analyses, multivariable-adjusted models controlled for female’s age, sex, female’s educational attainment (never attended school, completed primary or less, completed higher than primary), if female is married or living with partner (yes, no), wealth index (poorest, poor, wealthy and wealthiest), female’s religion (orthodox Christian or not orthodox Christian), household size (continuous) and female employed (yes, no). Carbohydrate substitution models controlled for the same covariates as the multivariable-adjusted models plus percentage of calories from total fat (continuous). Carbohydrate substitution models with plant protein as the exposure additionally adjusted for percentage energy from animal protein and models with animal protein as the exposure additionally adjusted for percentage energy from plant protein. Animal and plant protein substitution models were also included to estimate the effect of substituting animal for plant protein (or vice versa) by adjusting for all socio-demographic covariates in the multivariable-adjusted models, plus percentage of calories from carbohydrate (continuous) and percentage of calories from fat (continuous).
Secondary analyses examined the associations between percentage energy from total protein, plant protein and animal protein and continuous BMI separately excluding underweight participants and excluding overweight participants as the relationship between protein intake and BMI may differ in the context of overfeeding and underfeeding(Reference Westerterp-Plantenga, Lemmens and Westerterp36). Secondary analyses were also performed adjusting for total energy intake as a continuous variable in the linear models examining percentage energy from total protein, plant protein and animal protein with continuous BMI; as well as the logistic models with binary outcomes of overweight and underweight. In addition, because 29 % of participants reported currently lactating, secondary analyses were conducted adjusting for lactation status to rule out potential confounding due to postpartum weight retention. To confirm our estimation of substitution effects using the nutrient density models, we also estimated substitutions using the energy partition model(Reference Willett34) by including percentage of energy from fat, carbohydrate and protein as continuous variables in the same multivariable model. For this secondary analysis, the difference in the coefficients for each macronutrient plus their covariance was used to estimate substitution effects(Reference Willett34).
A two-sided probability value < 0·05 was considered to indicate a statistically significant difference. Statistical analyses were conducted using Stata 16·1 (StataCorp LP).
Results
Characteristics of the survey population are presented in Table 1. After excluding participants with implausible energy intakes (n 5) and pregnant women (n 56), there were 632 men and 992 women included in the final analysis. Most participants were Orthodox Christian, had completed primary school or higher and were married. According to the Household Food Insecurity Access Scale, 38 % of households were food insecure(Reference Coates, Swindale and Bilinsky31). Female’s age, male sex, Orthodox Christian (v. not), married (v. not), household size, household wealth quartile, household food security and BMI were all positively associated with percentage energy intake from total protein in bivariate analyses (P < 0·05; Table 1). Median energy intake (P25, P75) was 1695 (1353, 1986) kcal/d for men and 1556 (1223, 1890) kcal/d for women. Mean protein intake was 14 % of total energy for men and 14 % of total energy for women. Mean percentage energy from carbohydrates and fat was 69 % and 17 %, respectively. The proportion of participants who did not meet the estimated average requirement for protein intake was 21 % of men, 23 % of non-lactating women and 62 % of lactating women. Plant protein made up 83 % of total protein intake in men and 84 % in women. Major sources of dietary protein for both men and women included teff (36 % of protein), beef meat (12 %), peas (10 %), lentils (9 %) and wheat (5 %). In this sample, 10 % of women and 8 % of men were underweight, 24 % of women and 23 % of men were overweight and 10 % of women and 3 % of men were obese.
* P < 0·05.
** P < 0·001.
† Demographic characteristics are presented by quartiles of percent energy from protein intake, with the first quartile representing the lowest percent energy intake from protein.
‡ Data are presented as mean ± sd or median (P25, P75) for continuous measures and % (n) for categorical measures.
§ EAR for protein intake are 0·66 g/kg/d for adult men and women aged 19–70 years.
|| EAR for protein intake are 1.05 g/kg/d for lactating women.
In linear models adjusted for socio-demographic covariates, percentage energy from total protein was not associated with continuous BMI (Table 2). A greater percentage energy from plant protein was associated with lower BMI in linear models adjusted for age and sex, but this was NS after controlling for covariates. In age- and sex-adjusted linear models, a greater percentage energy from animal protein was positively and significantly associated with increased BMI (β for 1 % energy increment 0·15; 95 % CI: 0·07, 0·23; P < 0·001). However, this was NS after adjusting for socio-demographic characteristics. Of all covariates, wealth and education were strong confounders of this relationship and inclusion of either in the model nullified the positive association between animal protein and BMI. Additionally adjusting for percentage energy from fat (which changes the model into a substitution of animal protein for carbohydrate) flipped the association to inverse, although this was NS (β for 1 % energy increment –0·09; 95 % CI: –0·22, 0·05; P = 0·20) (Table 2). There was no evidence of nonlinearity for the associations between BMI and percentage energy from total protein, plant protein and animal protein in multivariate-adjusted models (see online supplementary material, Supplemental Fig. 1) or multivariate-adjusted carbohydrate substitution models (see online supplementary material, Supplemental Fig. 2).
* P < 0·05.
** P < 0·001.
† Percentage energy from total protein, plant protein and animal protein were examined in linear models as continuous exposures and also in quartiles; Q1 represents the lowest percent energy intake from total, plant or animal protein.
‡ Beta coefficients are in kg/m2 per 1 % difference in percent energy
§ The age- and sex-adjusted model is adjusted for female’s age (years, continuous) and sex (female, male)
|| The multivariable-adjusted model is adjusted for female’s age, sex, female’s educational attainment (never attended school, completed primary or less, completed higher than primary), if female is married or living with partner (yes, no), wealth index (poorest, poor, wealthy, wealthiest), female’s religion (orthodox Christian or not orthodox Christian), household size (continuous) and female employed (yes, no).
¶ The multivariable-adjusted carbohydrate substitution model is the multivariable-adjusted model additionally adjusted for percentage of calories from total fat (continuous). When plant protein is the exposure, the model is additionally adjusted for percentage of calories from animal protein. When animal protein is the exposure, the model is additionally adjusted for percentage of calories from plant protein.
†† The multivariable-adjusted animal protein substitution model is adjusted for female’s age, sex, percentage of calories from fat, percentage of calories from carbohydrate (continuous), female’s educational attainment, if female is married or living with partner, wealth index, female’s religion, household size and female employed.
‡‡ The multivariable-adjusted plant protein substitution model is adjusted for the same covariates as the multivariable-adjusted animal protein substitution model.
Total protein intake was associated with lower odds of overweight/obesity (OR per 1 % energy increment 0·92; 95 % CI: 0·86, 0·99; P = 0·02) compared with healthy weight, after adjusting for age, sex, percentage energy from total fat, educational attainment, marital status, wealth index, religion, household size and employment status (Table 3). Participants in the third quartile of percentage energy from total protein had a lower odds of overweight/obesity compared with those in the lowest quartile (OR 0·69; 95 % CI: 0·49, 0·97; P = 0·04) when substituting protein for carbohydrate, but no significant associations were observed when comparing the highest and lowest quartiles (OR 0·89; 95 % CI: 0·62, 1·28; P = 0·54). These results were consistent after adjustment for lactation status. In fully adjusted models, total protein was associated with a lower odds of overweight/obesity when substituting for carbohydrate (OR per 1 % energy increment 0·88; 95 % CI: 0·81, 0·97; P = 0·009) in women, but no significant associations were observed in men (OR per 1 % energy increment 0·99; 95 % CI: 0·88, 1·12; P = 0·91) (see online supplementary material, Supplemental Table 1). In multivariable-adjusted carbohydrate substitution models, total protein intake was not significantly associated with odds of underweight (compared with healthy weight) in the overall sample (OR per 1 % energy increment 0·95; 95 % CI: 0·85, 1·07; P = 0·40), in women (OR per 1 % energy increment 0·92; 95 % CI: 0·80, 1·06; P = 0·24) or in men (OR per 1 % energy increment 1·09; 95 % CI: 0·89, 1·32; P = 0·39) (Table 3; see online supplementary material, Supplemental Table 1).
* P < 0·05
† Percentage energy intake from total protein, plant protein and animal protein were examined in logistic regression models as continuous exposures and also in quartiles; Q1 represents the lowest percentage energy intake from total, plant or animal protein.
‡ Overweight/obese is defined as BMI ≥ 25 and is compared with a reference category of non-overweight/obese (BMI < 25). All underweight participants (n 148) were excluded from the overweight analysis. Underweight is defined as BMI < 18·5 and is compared with a reference category of non-underweight (BMI ≥ 18·5). All overweight participants (n 501) were excluded from the underweight analysis.
§ The age- and sex-adjusted model is adjusted for female’s age (years, continuous) and sex (female, male)
|| The multivariable-adjusted model is adjusted for female’s age, sex, female’s educational attainment (never attended school, completed primary or less, completed higher than primary), if female is married or living with partner (yes, no), wealth index (poorest, poor, wealthy, wealthiest), female’s religion (orthodox Christian or not orthodox Christian), household size (continuous) and female employed (yes, no).
¶ The multivariable-adjusted carbohydrate substitution model is the multivariable-adjusted model additionally adjusted for percentage of calories from total fat (continuous). When plant protein is the exposure, the model is additionally adjusted for percentage of calories from animal protein. When animal protein is the exposure, the model is additionally adjusted for percentage of calories from plant protein.
** The multivariable-adjusted animal protein substitution model is adjusted for female’s age, sex, percentage of calories from fat, percentage of calories from carbohydrate (continuous), female’s educational attainment, if female is married or living with partner, wealth index, female’s religion, household size and female employed.
†† The multivariable-adjusted plant protein substitution model is adjusted for the same covariates as the multivariable-adjusted animal protein substitution model.
Among men and women combined, plant protein intake was not significantly associated with odds of underweight or overweight/obesity compared with healthy weight (Table 3). However, plant protein intake was associated with increased odds of overweight/obesity in women when substituted for animal protein (OR per 1 % energy increment 1·11; 95 % CI: 1·01, 1·22; P = 0·03), but not when substituted for carbohydrate (OR per 1 % energy increment 0·96; 95 % CI: 0·86, 1·07; P = 0·50) (see online supplementary material, Supplemental Table 1). Greater animal protein intake was significantly associated with lower odds of overweight in fully adjusted models when substituted for total carbohydrate (OR per 1 % energy increment 0·89; 95 % CI: 0·82, 0·96; P = 0·004) and for plant protein (OR per 1 % energy increment 0·91; 95 % CI: 0·85, 0·98; P = 0·008) (Table 3). Participants in the highest quartile of animal protein intake had a 50 % (95 % CI: 19, 69; P = 0·005) lower odds of overweight/obesity compared with those in the lowest quartile, substituting animal protein for carbohydrate. These results were consistent after adjustment for lactation status. Women had 68 % (95 % CI: 41, 83; P < 0·001) lower odds of overweight/obesity comparing the highest to lowest quartile of animal protein intake, substituted for carbohydrate, but this association was not seen in men. No associations between animal protein intake and underweight were observed in the fully adjusted models overall, in women, or in men (Table 3; see online supplementary material, Supplemental Table 1).
In food-based analyses, percentage energy from most commonly consumed protein sources was not related to underweight or overweight/obesity (Table 4). In fully adjusted models, a greater percentage energy from milk was associated with lower odds of underweight but the CI was wide (OR per 5 % energy increment 0·54; 95 % CI: 0·30, 1·00; P = 0·05) and this finding was NS after adjusting for lactation status (OR per 5 % energy increment 0·55; 95 % CI: 0·30, 1·01; P = 0·06). Although not statistically significant, a greater percentage energy intake from barley was related to lower odds of overweight/obesity (OR per 5 % energy increment 0·85; 95 % CI: 0·72, 1·01; P = 0·06).
* P < 0·05.
† Protein source foods were included in models as percentage of total energy intake from each food. Models are examining a 5 % increase in percentage of total energy intake of each food.
‡ Overweight/obese is defined as BMI ≥ 25. The overweight/obese analysis was conducted among 1476 healthy weight, overweight and obese participants. 148 underweight participants were excluded from this analysis.
§ Underweight is defined as BMI < 18·5. The underweight analysis was conducted among 1123 healthy weight and underweight participants. 501 overweight and obese participants were excluded from this analysis.
|| Age- and sex-adjusted model is adjusted for female’s age (years, continuous) and sex (female, male)
¶ Multivariable-adjusted model is adjusted for female’s age, sex, female’s educational attainment (never attended school, completed primary or less, completed higher than primary), if female is married or living with partner (yes, no), wealth index (poorest, poor, wealthy, wealthiest), female’s religion (orthodox Christian or not orthodox Christian), household size (continuous) and female employed (yes, no).
For linear models with BMI as an outcome, results from secondary analyses adjusting for total energy intake (see online supplementary material, Supplemental Table 2) and excluding overweight participants (n 501) (see online supplementary material, Supplemental Table 3) were largely consistent with the primary analyses. When excluding all underweight participants (n 148), total protein was weakly associated with lower BMI when substituted for carbohydrate comparing the third and lowest quartiles of total protein intake (β –0·55; 95 % CI: –1·08, –0·01; P = 0·05), but not when comparing the highest and lowest quartiles (β –0·07; 95 % CI: –0·19, 0·04; P = 0·22) (see online supplementary material, Supplemental Table 4). Animal protein intake was related to lower BMI when substituted for plant protein in analyses excluding underweight participants, but this association was marginally significant (β per 1 % energy increment –0·11; 95 % CI: –0·23, 0·00; P = 0·05; P trend across Q1–Q4 = 0·04).
In secondary analyses adjusting for total energy intake with overweight/obesity as an outcome, the association between total protein and overweight/obesity was NS in the fully adjusted carbohydrate substitution model (OR per 1 % energy increment: 0·93; 95 % CI: 0·87, 1·00; P = 0·06) (see online supplementary material, Supplemental Table 5). However, when adjusting for total energy intake, greater percentage energy from animal protein was still associated with lower odds of overweight/obesity in fully adjusted models substituting animal protein for carbohydrate (OR per 1 % energy increment: 0·90; 95 % CI: 0·83, 0·98; P = 0·01) and for plant protein (OR per 1 % energy increment: 0·91; 95 % CI: 0·85, 0·98; P = 0·01).
Results from secondary analyses using multivariable-adjusted energy partition models to examine substitution effects, which did not control for total energy intake, were consistent with the primary analyses. Among all participants, substituting fat for carbohydrate was associated with higher BMI, but substituting protein for carbohydrate was not. Among overweight/obese and healthy weight participants, substituting percentage energy from fat for carbohydrate was associated with higher odds of overweight/obesity, while substitution of protein for carbohydrate and protein for fat was associated with lower odds of obesity (see online supplementary material, Supplemental Table 6).
Discussion
In this cross-sectional analysis among 1624 adults in Addis Ababa, Ethiopia, we found that percentage energy intake from total protein and animal protein were both significantly associated with lower odds of overweight/obesity when substituted for carbohydrate, especially among women. This study population had a very high intake of carbohydrates and very low intake of animal protein, which is consistent with other studies in low-income countries(Reference Schönfeldt and Gibson Hall8). In addition, one-third of the study population was overweight or obese (34 % of women and 25 % of men), which is similar to previous estimates in Addis Ababa(2,7) .
To our knowledge, this is the first study in SSA to examine the associations between total, animal and plant protein intake in relation to underweight and overweight/obesity. Our findings are consistent with some previous studies in Western and Asian populations that have found that increased total protein as well as animal protein were associated with improved waist circumference(Reference Park, Park and Kang18) and body composition(Reference Berryman, Agarwal and Lieberman17). Some randomised trials comparing higher protein diets to lower protein diets have reported decreases in body weight(Reference Green, Shea and Vasdev11,Reference Foster, Wyatt and Hill37–Reference Dong, Zhang and Wang40) and greater fat loss with higher protein diets(Reference Green, Shea and Vasdev11,Reference Samaha, Iqbal and Seshadri38,Reference Halton and Hu41) . However, other trials have found that higher protein diets did not result in increased weight loss compared with other diets, especially over the long term(Reference Halton and Hu41–Reference Noakes, Keogh and Foster43).A meta-analysis of thirty-two trials with greater than 12 months of follow up found that although higher protein diets showed benefits for weight loss in the short term, these benefits persisted to a smaller degree in the long term. However, greater benefits were observed long term in those with better compliance to higher-protein diets(Reference Clifton, Condo and Keogh44).
Although additional evidence is needed from long-term trials, existing evidence from shorter-term trials suggests that higher-protein diets may have beneficial metabolic effects when protein is used to partially replace carbohydrates (especially carbohydrates from refined sources)(Reference Halton and Hu41). Replacement of carbohydrate with protein may promote diet-induced thermogenesis, increase satiety which can lead to reduced subsequent energy intake and lower glycaemic load(Reference Halton and Hu41,Reference Weigle, Breen and Matthys45) . Protein intake may also protect against loss of lean body mass(Reference Liu, Silvestre and Poppitt14). In the Ethiopian context, increasing protein intake in place of carbohydrate may be beneficial because carbohydrate intake is very high and food sources of carbohydrates are often refined(7). In our sample, 69 % of calories were derived from carbohydrates and the major sources were teff (a cereal common in Ethiopia), refined wheat, maize and pasta.
Interestingly, we also found percentage energy from animal protein was associated with lower odds of overweight/obesity when substituted for plant protein, which conflicts with some previous studies in Western populations that have found beneficial effects of plant protein and harmful effects of animal protein on waist circumference, BMI and metabolic syndrome(Reference Lin, Bolca and Vandevijvere19–Reference Shang, Scott and Hodge21). Food sources of animal protein tend to be higher in cholesterol, energy and saturated fats while food sources of plant protein could help to control body weight and improve body composition because of their associations with lower intakes of energy, total fat, cholesterol and saturated fat(Reference Lin, Bolca and Vandevijvere19). However, these findings and mechanisms in Western populations are likely not generalisable to low-income settings, especially to Ethiopia, where animal-source food consumption is much lower due to economic, cultural and religious factors and carbohydrate intake is very high(Reference Zellelew46). The study populations in these Western studies were consuming 9–11 % of total energy from animal protein(Reference Lin, Bolca and Vandevijvere19,Reference Shang, Scott and Hodge21) , compared with 3 % of energy from animal protein in our Ethiopian study population. In addition, much of the protein intake in our study was from cereals (such as teff, wheat and rice) which may have less favourable effects on body weight than non-cereal plant proteins (such as pulses and nuts). A randomised controlled trial of protein intake and weight loss maintenance found that while substituting overall plant for animal protein was not associated with effects on body weight, a higher intake of cereal plant protein at the cost of non-cereal plant protein was associated with a larger increase in body weight(Reference van Baak, Larsen and Jebb47). Non-cereal plant protein sources such as pulses tend to have a lower glycaemic load and may increase satiety, compared with cereal protein sources (some of which may be refined carbohydrates)(Reference van Baak, Larsen and Jebb47). It is possible that in low-income settings with high intakes of carbohydrate and low consumption of animal-source foods, protein from animal sources may be beneficial in preventing overweight or obesity, given that much of the protein and energy in the diet comes from cereals.
A strength of this study is that we were able to approximate usual dietary intake through a non-quantitative FFQ adapted for use in Ethiopia. Most previous dietary surveys in Ethiopia have relied on a single 24-hour recall, which are not considered representative of usual dietary patterns. We also conducted a detailed assessment of confounders related to protein intake, overweight/obesity and underweight. Consistent with previous studies, we found a number of socio-demographic factors that were strongly related to protein intake including age, sex, religion, marital status, wealth and household food insecurity(Reference Steyn and McHiza48). It is well-documented that these confounders are also strongly associated with obesity risk in sub-Saharan Africa(Reference Tebekaw, Teller and Colón-Ramos5,Reference Steyn and McHiza48) .
FFQs are designed primarily to generate estimates of relative intake, which gives us the ability to distinguish between high and low consumers of a given nutrient, but they can only approximate absolute intake and are subject to additional limitations including recall and self-report biases. We were able to use a non-quantitative FFQ in this study because our analyses did not rely on absolute intakes. While it would have been ideal to use an FFQ validated in an Ethiopian population, no such instrument exists to our knowledge; therefore, we adapted this FFQ from one validated in Tanzania. Additionally, we were unable to capture seasonal variations in food intake because this FFQ was administered at a single time point. During the time period of the survey (Jan–Feb), urban energy intake tends to be higher than during the lean season in June and July(Reference Hirvonen, Taffesse and Worku Hassen49). Despite this, we found that a large proportion of our sample were below the estimated average requirement for protein intake (21 % of men, 23 % of non-lactating women and 62 % of lactating women). However, these findings should be interpreted with caution because the FFQ was not validated in this population and the low levels of protein intake we found are in part due to the low total energy intake in this sample, which may have been under-reported. That said, levels of intake reported by our study participants were consistent with previous estimates from Ethiopia(7).
An additional limitation of this study is possible residual confounding from measured and unmeasured factors, such as physical activity. We also lack some key covariates for the men in the sample (including age), but were able to use the female’s covariates as the male and female participants were living in the same household and in most cases were married. However, it is possible that age, education and employment of the female participant could be different than that of the male participant, which could lead to bias in the results of the adjusted models for the male participants. There is also potential for selection bias, since the study only included households who consented to participate in the survey and contained young children, which may not be representative of the overall population of Addis Ababa. However, this survey had a very high response rate (97 %) and the characteristics of the study population are similar to those of the greater urban Ethiopian population in terms of age, educational attainment and household size(2). In addition, Ethiopia has a high fertility rate and a long median duration of breastfeeding (2 years) which may also partially explain the high proportion of lactating women in this sample(2). To further address these limitations, we adjusted for socio-demographic and household characteristics in all analyses and adjusted for lactation status in secondary analyses. Lastly, this study focused on Addis Ababa, a major city and results may not be generalisable to rural areas of Ethiopia or other SSA contexts, where dietary patterns and prevalence of overweight differ. In Addis Ababa, consumption of animal source foods is twice as high as in rural areas of Ethiopia(Reference Abegaz, Hassen and Minten50) and overweight/obesity prevalence is six times higher among women and thirteen times higher among men in urban compared with rural areas of Ethiopia(2).
The cross-sectional design of this study is an additional limitation as it restricts our ability to make causal inferences and poses unique challenges when studying bodyweight as an outcome because participants may change their diets due to awareness of their weight. Repeated assessments of changes over time in dietary factors and body weight are needed to understand long-term effects of protein intake on overweight/obesity and underweight(Reference Mozaffarian, Hao and Rimm35). The cross-sectional design could explain why we did not observe any associations between underweight and animal or plant protein intake.
Animal protein has traditionally been thought of as beneficial in the context of undernutrition, as it provides higher-quality protein (defined by the amino acid content) than most plant sources, along with other essential nutrients(Reference Schönfeldt and Gibson Hall8,Reference Willett, Rockström and Loken15) . Although animal-source foods may be valuable in preventing and treating undernutrition, promotion of animal-source foods, especially meat, is controversial as they have greater environmental impact than plant protein sources(Reference Springmann, Godfray and Rayner51). One study in Addis Ababa found that meeting adult and child recommended daily protein intakes with plant-based foods would have a lower environmental impact than meeting this gap with animal-based foods, but would still result in an estimated 65 % increase in land and water use and a 2 % increase in greenhouse gas emissions. Meeting this protein gap with a milk and red meat strategy would result in an estimated 190 % increase in land and water use and 257 % increase in greenhouse gas emissions(Reference Blakstad, Danaei and Tadesse10).
Despite their environmental impacts, animal-source foods are a valuable source of protein and micronutrients and providing some animal-source foods alongside plant-based foods may offer additional health benefits in the Ethiopian context(Reference Blakstad, Danaei and Tadesse10). In Africa, animal-source food production can be sustainable through using grassland which cannot be used for crop production and by converting inedible crop by-products into edible food. It can also aid in reducing fertiliser use, cycling nutrients within the ecosystem and is an important mechanism to diversify farmers’ income in the case crop production is reduced(Reference Kim, Santo and Scatterday52–54). The Ethiopian Government has demonstrated their commitment to transitioning traditional livestock practices to climate-smart practices through policy interventions as outlined in their Livestock Master Plan and Climate Resilience and Green Economy Strategy(54).
For optimal human and planetary health, the EAT-Lancet Commission has recommended adopting a global reference diet, consisting of mainly plant-source foods (vegetables, fruits, whole grains, and legumes and nuts), and limited or modest amounts of animal-source foods (such as red or processed meat)(Reference Willett, Rockström and Loken15). However, in our sample, animal-source protein intake was only 11 g per day, which is less than half of the world average of animal-source protein availability in 2011 and much less than the amounts in the EAT-Lancet healthy reference diet(Reference Willett, Rockström and Loken15). In Ethiopia, animal protein consumption tends to be limited due to religious and social norms and the availability, accessibility and affordability of animal-source foods(Reference Bliznashka, Passarelli and Canavan55). It is estimated that to achieve global dietary recommendations in SSA, increases in food consumption-related greenhouse gas emissions would be necessary, but that these small increases would be far outweighed by higher-income countries adopting diets lower in meat consumption(Reference Springmann, Godfray and Rayner51).
Conclusion
Innovative strategies are needed to combat the complex, coexisting problems of over and undernutrition in SSA. Our results suggest that in the context of diets that are very high in carbohydrate such as those in Ethiopia, a greater proportion of energy intake from protein may be associated with lower prevalence of overweight and obesity. However, due to the cross sectional design of this study and limited generalisability to other contexts, our findings should be interpreted with caution. Community-based longitudinal data is scarce in Ethiopia, but such studies are needed to further clarify associations between dietary intake with underweight and overweight/obesity, track changes over time and account for seasonality and urban/rural differences in consumption. The majority of the literature on dietary protein and bodyweight is from Western settings and is likely not generalisable to SSA, where dietary patterns are very different. Longitudinal dietary surveys in Ethiopia are needed to inform local guidelines, policies and programmes to combat the double burden of malnutrition. Given the constraints limiting animal protein consumption in Ethiopia, future research should also assess the feasibility of incorporating modest amounts of animal products into Ethiopian diets.
Acknowledgements
Acknowledgements: The authors would like to thank the study participants and data collectors. Financial support: This work was supported by a United Kingdom Economic and Social Research Council (ESRC) and (DFID) Joint Fund for Poverty Alleviation Research (grant number: ES/R002118/1). Sabri Bromage received support from the National Institutes of Health (grant number: D43 TW010543). Authorship: E.C.H. and W.W. conceptualised the analyses presented in this paper. Y.B., A.W.T., C.R.C. and W.F. designed and implemented the Addis Ababa study. A.W.T. and Y.B. led the data collection activities. E.C.H. analysed the data and wrote the manuscript. S.B. advised on data analysis and critically revised the manuscript for important intellectual content. All authors have read and approved the final manuscript and are accountable for all aspects of the work. Ethics of human subject participation: This study was conducted according to the guidelines laid down in the Declaration of Helsinki and all procedures involving research study participants were approved by the Institutional Review Board of the Harvard T.H. Chan School of Public Health (reference number: IRB17-1825) and by Addis Continental Institute of Public Health (reference number: IRB17-1825). Written informed consent was obtained from all subjects/patients.
Conflict of interest:
The authors declare no conflicts of interest.
Supplementary material
For supplementary material/s referred to in this article, please visit https://doi.org/10.1017/S1368980022001100