Reduced Rank Regression-Derived Dietary Patterns Related to the Fatty Liver Index and Associations with Type 2 Diabetes Mellitus among Ghanaian Populations under Transition: The RODAM Study
by
, , , , , , ,
Tracy Bonsu Osei
1
,
Anne-Marieke van Dijk
2,
Sjoerd Dingerink
2,
Felix Patience Chilunga
3,
Erik Beune
3,
Karlijn Anna Catharina Meeks
3,4,
Silver Bahendeka
5,
Matthias Bernd Schulze
6,7,
Charles Agyemang
3
,
Mary Nicolaou
3,
Adriaan Georgius Holleboom
2 and
Ina Danquah
1,6,*
1
Heidelberg Institute of Global Health (HIGH), Faculty of Medicine and University Hospital, Heidelberg University, 69120 Heidelberg, Germany
2
Department of Vascular Medicine, Amsterdam UMC, University of Amsterdam, 1105 AZ Amsterdam, The Netherlands
3
Department of Public Health, Amsterdam UMC, University of Amsterdam, 1105 AZ Amsterdam, The Netherlands
4
Center for Research on Genomics and Global Health, National Human Genome Research Institute, Bethesda, MD 20892-2152, USA
5
Department of Internal Medicine, St. Francis Hospital Nsambya, MKPGMS-Uganda Martyrs University, Kampala 5498, Uganda
6
Department of Molecular Epidemiology, German Institute of Human Nutrition Potsdam-Rehbruecke, 14558 Nuthetal, Germany
7
Institute of Nutritional Science, University of Potsdam, 14558 Nuthetal, Germany
*
Author to whom correspondence should be addressed.
Nutrients 2021, 13(11), 3679; https://doi.org/10.3390/nu13113679
Submission received: 27 August 2021
/
Revised: 8 October 2021
/
Accepted: 16 October 2021
/
Published: 20 October 2021
(This article belongs to the Section Nutrition and Public Health)
Abstract
:The Fatty Liver Index (FLI) is a proxy for the steatotic component of non-alcoholic fatty liver disease (NAFLD). For sub-Saharan African populations, the contribution of dietary factors to the development of NAFLD in the etiology of type 2 diabetes mellitus (T2DM) remains to be clarified. We identified sex-specific dietary patterns (DPs) related to the FLI using reduced ranked regression (RRR) and evaluated the associations of these DPs with T2DM. This analysis used data from the RODAM, a multi-center cross-sectional study of Ghanaian populations living in Ghana and Europe. The daily intake frequencies of 30 food groups served as the predictor variables, while the FLI was the response variable. The odds ratios and 95% confidence intervals for T2DM were calculated per one standard deviation increase in the DP score using logistic regression. In males, the DP score explained 9.9% of the variation in their food intake and 16.0% of the variation in the FLI. This DP was characterized by high intakes of poultry, whole-grain cereals, coffee and tea, condiments, and potatoes, and the chance of T2DM was 45% higher per 1 DP score-SD (Model 2). Our results indicate that the intake of modernized foods was associated with proxies of NAFLD, possibly underlying the metabolic pathways to developing T2DM.
1. Introduction
In sub-Saharan Africa (SSA) and in African migrants living in Europe, type 2 diabetes mellitus (T2DM) constitutes a major health problem [1,2], and the relationships with non-alcoholic fatty liver disease are not fully understood [3]. NAFLD is a spectrum of liver disease primarily characterized by the accumulation of fat in the liver in the absence of excessive alcohol consumption and other chronic liver diseases, such as viral infection [4,5,6]. It is a common liver disorder with an estimated global prevalence of 25% [5,7]. NAFLD strongly co-occurs in patients with components of metabolic syndrome, especially central obesity and T2DM. Indeed, up to 70% of patients with T2DM have NAFLD, and they stand a higher risk of developing non-alcoholic steatohepatitis and advanced fibrosis [6,8,9]. In fact, it is estimated that about 68–74% of individuals with T2DM have evidence of steatosis, and approximately 15% may have advanced fibrosis [9].
To date, the rapid economic transition and urbanization in SSA as well as the sudden lifestyle changes upon migration from Africa to Europe are among the proposed reasons for the upsurge of T2DM among these populations [10]. The parallel nutrition transition comprises shifts in dietary practices from plant-based, fiber-rich and low-fat diets to modernized diets rich in saturated fats and simple carbohydrates. These dietary changes may play a key role in the development of T2DM and NAFLD among SSA populations [11]. Diets high in refined carbohydrates and saturated fats are associated with increased intrahepatic triglycerides, obesity and being overweight, insulin resistance and NAFLD [12,13]. Conversely, the Mediterranean diet, which is characterized by high intakes of fiber-rich foods such as vegetables and whole grains, improves liver function [14,15]. Although dietary interventions are essential in the prevention and primary management of NAFLD, the contributions of dietary factors to the development of NAFLD and the etiology of T2DM remain to be clarified, particularly in SSA populations with a high prevalence of obesity and T2DM [1,2].
Here, we present two methodological advancements to study the etiologic pathways from a diet through NAFLD up to T2DM in SSA populations under various transitions. First, the method of reduced rank regression (RRR) is applied to derive dietary patterns (DPs) related to the proxy markers of NAFLD. RRR determines linear combinations of food groups as the predictor variables by maximizing the explained variation in these proxy markers as the response variables [16]. This approach advances the exploratory nature of principal component analysis and the hypothesis-driven techniques of constructing a priori indices [17,18]. Second, the Fatty Liver Index (FLI) is used as a non-invasive proxy for hepatic steatosis in NAFLD, incorporating both anthropometric and biochemical markers [19,20], and it seems to perform well at any degree of obesity [21].
Considering that NAFLD and T2DM frequently co-occur within subjects and share similar risk factors, we hypothesized that RRR-derived DPs related to NAFLD and characterized by modernized foods are associated with T2DM. Therefore, this study aimed to establish sex-specific DPs related to FLI by RRR and determine the associations between these DPs and T2DM among Ghanaian populations under transition.
2. Materials and Methods
2.1. Study Design and Population
The Research on Obesity and Diabetes among African Migrants (RODAM) study is a multi-center cross-sectional study that aims to investigate the relative contributions of lifestyle and other risk factors to obesity and T2DM among African migrants in Europe and their counterparts living in rural and urban Ghana. A detailed description of the study design and procedures has been presented elsewhere [22]. In brief, participant recruitment was carried out between July 2012 and September 2015, and it included 6385 Ghanaian adults between the ages of 18 and 96 years living in rural Ghana (n = 946), urban Ghana (n = 1619), Amsterdam (n = 1900), London (n = 1258) and Berlin (n = 662). For recruitment in Ghana, the census data of 2010 were used to draw urban (Kumasi and Obuasi) and rural participants from the Ashanti Region. The response rates in rural and urban Ghana were 76% and 74%, respectively. In Amsterdam, random selection was performed by using the municipal register to select Ghanaian migrants, who were invited by postal mail and home visits. In Amsterdam, 67% replied by response card or after a home visit, and of these, 53% of the respondents consented to participate in the study. In London, of those individuals who were invited based on their registration in Ghanaian organizations, 75% agreed and participated in the study. Berlin had a participation rate of 68%.
The study protocol was reviewed and approved by the local ethics committees in Ghana (School of Medical Sciences, Komfo Anokye Teaching Hospital, Committee on Human Research, Publication and Ethical Review Board), the Netherlands (Institutional Review Board of AMC, University of Amsterdam), the UK (London School of Hygiene and Tropical Medicine Research Ethics Committee) and Germany (Ethics Committee of Charite-Universitatsmedizin Berlin). All participants who took part in the study were thoroughly informed about the purpose and procedures of the study and gave written informed consent.
2.2. Assessment of NAFLD Proxy Markers
Fasting venous blood samples were collected by trained personnel according to standard operating procedures. Serum and plasma samples were obtained through centrifugation and stored at −80 °C until analysis. The enzymatic method (hexokinase) was used to measure the fasting blood glucose in the plasma. All biochemical analyses were performed using an ABX Pentra 400 chemistry analyzer (ABX Pentra; Horiba ABX, Germany). T2DM was defined according to the WHO guidelines as fasting plasma glucose ≥7.0 mmol/L or documented use of glucose-lowering medication or self-reported diabetes. The serum lipid profile (total cholesterol, triglyceride, high-density lipoprotein cholesterol and low-density lipoprotein cholesterol) was measured using immunoturbidimetric assays. The biomarkers of inflammation and liver function comprised C-reactive protein (CRP), aspartate aminotransferase (AST), alanine aminotransferase (ALT) and γ-glutamyl transferase (GGT). We used two sets of proxy markers to define NAFLD. These were the FLI as a single marker and a selection of biomarkers as the second set of proxy markers. The FLI was calculated for each participant according to the algorithm by Bedogni et al. (2006) [23]. The formula considers the BMI, waist circumference, triglycerides and GGT. In the absence of excessive alcohol intake (>21 units (168 g) per week for men and >14 units (112 g) per week for women), the FLI supports the definition of NAFLD in our epidemiological study.
2.3. Dietary Assessment
A semi-quantitative Ghana-specific food propensity questionnaire (Ghana-FPQ) was used to determine the usual dietary intake of 134 food items. The usual weekly intake frequencies of food groups were documented in pre-defined portion sizes for the past 12 months. We calculated the amounts of food intake (g/d) using common household measures and translated these into energy intake and nutrient consumption by applying the West African Food Composition Table (2012) and the German Nutrient Database (2010). The food items were grouped into 30 food groups according to their culinary use and nutrient profiles [10]. The food groups subjected to RRR included whole-grain cereals, refined cereals, sweet spreads, dairy products, fruits, nuts and seeds, roots, tubers and plantains, potatoes, fermented maize products, vegetables, legumes, vegetable soups, stews, sauce, rice and pasta, eggs, red meat dishes, cakes and sweets, coffee and tea, alcoholic beverages, sodas and juices, palm oil, olive oil, other oils, margarine, cooking fats and condiments.
2.4. Assessment of Covariates
All information on the socio-demographics, education levels, medical histories and lifestyle factors was collected using a structured questionnaire, either self-administered or by an interviewer in the preferred language. The educational status of the participants was categorized into four levels: never or elementary, lower, intermediate and higher or tertiary. Physical activity was assessed using the WHO surveillance questionnaire, and the answers were subsequently classified based on the guidelines of the IPAQ group into three levels of total physical activity: low, moderate and high [24]. Smoking status was classified as current, former or non-smokers. At all study sites, physical examinations were conducted with validated devices according to standard operating procedures and by trained personnel. For all the study participants, anthropometric measurements were taken in light clothes and without shoes. These included weight (kg) by a person scale, height (cm) by a stadiometer, and waist circumference (cm) using a measuring tape (all devices SECA, Germany). We calculated the Body Mass Index (BMI) as the weight over the height squared (kg/m2).
3. Statistical Analysis
3.1. Descriptive Statistics
Data processing and analysis was performed using the Statistical Analysis Software (SAS) (version 9.4). Participants with missing data for the variables of interest were excluded from the present analysis. Details of the exclusions are presented in Supplementary Figure S1. The final analytical sample comprised 3687 participants. Normally distributed continuous variables are presented as means ± standard deviations (SDs) with the skewed variables as medians (interquartile ranges). All categorical variables are presented as percentages.
3.2. Reduced Rank Regression
Owing to the differences in dietary intakes, adipose tissue distribution, and observed T2DM prevalence between males and females, all subsequent analyses were conducted separately for males and females. First, we used RRR to derive dietary pattern (DP) scores related to the FLI as a proxy marker of NAFLD. The RRR method determined the linear combinations of predictor variables (= 30 food groups) that explained as much variation as possible in the response variable (log transformed FLI). This method has been explained in detail by Hoffmann et al. (2004) [25]. Second, we created RRR-derived DP scores using established biomarkers of NAFLD as the response variables (i.e., AST, ALT, GGT, CRP, total cholesterol, HDL-cholesterol, LDL-cholesterol and triglycerides). Hence, four different DPs were derived that explained a maximum of variation in the FLI (for men and women) and in the biomarkers of NAFLD (for men and women).
3.3. Association Analyses
The associations between the RRR-derived DPs and T2DM were calculated by logistic regression analyses. The odds ratios and their 95% confidence intervals (CIs) for T2DM were calculated across the quintiles of the RRR-derived DP scores by using the first quintile as the reference category and per 1 score-SD. To assess the linearity assumption, we calculated the odds ratios and their 95% confidence intervals (Cls) across the quintiles of the RRR-derived DP scores, using the first quintile as the reference category. Linear trends across the quintiles were tested by modeling the median of each quintile as a continuous variable. For the associations with the FLI-related DP scores, we sequentially adjusted for the age (years) and study site (categorical) (Model 1) as well as the education level (4 categories), energy intake (kcal/d), smoking status (yes or no), physical activity (METs-h/week), and alcohol consumption (g/day) (Model 2). For the associations with the DPs related to the biomarkers, we additionally adjusted for the BMI and waist circumference (Model 3).
3.4. Sensitivity Analysis
In sensitivity analyses, we excluded participants with self-reported T2DM to address a potential source of reverse causation. Specifically, individuals with known T2DM who received medical treatment might have presented altered biomarkers and could have changed their diets. Consequently, we repeated all RRR analyses using the FLI or NAFLD biomarkers as response variables for both males and females. In addition, we repeated the final regression analysis for the association between RRR-derived DPs and T2DM, applying an HbA1c-based definition of the outcome, namely HbA1c ≥ 48 mmol/mol, self-reported diabetes or the use of glucose-lowering medication.
4. Results
4.1. Study Population
Table 1 shows the socio-demographic, anthropometric, lifestyle, and clinical characteristics of the RODAM study population by study site and sex. The majority of the participants were females (63.0%) and middle-aged (46.1 ± 11.1 years). The males were older than the females, but the mean age was similar across the study sites (Table 1). Males had higher attained educations than females, and the proportion of participants with higher education was highest in London. Males were more physically active, more likely to be former or current smokers and reported higher alcohol consumption than females. The mean BMI (24.8 ± 4.4 kg/m2 vs. 27.8 ± 5.7 kg/m2) and waist circumference (87 ± 12.1 cm vs. 91 ± 12.5 cm) were lower in males than in females. The participants in London showed the lowest level of physical activity (median: 28; IQR: 5–112 min/day) and the highest BMI (29.4 ± 4.8 kg/m2) and waist circumference (95.4 ± 11.3cm), while rural Ghana showed the lowest BMI (22.7 ± 4.3 kg/m2) and waist circumference (81.2 ± 10.9 cm) and the highest physical activity level (median: 90; IQR: 36–161 min/day).
4.2. Intakes of Energy, Nutrients and Food Groups
The estimated mean energy intake was 2533 ± 837 kcal/d, and this was higher among males than females and highest in Berlin, followed by London, rural Ghana, Amsterdam and urban Ghana (Table 1). The mean contributions to the daily energy intake were 53.0 ± 9.1 energy% for carbohydrates, 32.2 ± 8.2 energy% for total fats and 13.8 ± 2.9% for protein. This was similar between males and females. The contribution of macronutrients to energy intake varied across study sites. In rural and urban Ghana, carbohydrates were the major contribution to the energy percentage, while in London and Amsterdam, fats and proteins contributed the most to energy intake. The consumption of certain foods like coffee and tea, sodas and juice and alcohol was higher in males compared with females (Figure 1a). This was also seen in red meat, poultry, processed meat and eggs, while women consumed more vegetarian mixed dishes, cakes and sweets than men (Figure 1b).
4.3. RRR-Derived Dietary Patterns Related to the Fatty Liver Index
Table 2 shows the explained variation and the factor loadings of those food groups related to the FLI as derived by RRR for males and females. The DP scores explained 9.9% and 6.5% of the total variation in food group intakes in males and females, respectively. Furthermore, the RRR-derived pattern scores were positively related with the FLI and explained 16.0% of the FLI variation among males (beta = 0.4) and 8.8% of the FLI variation among females (beta = 0.3). The DPs scores in males were characterized by high intakes of poultry, whole-grain cereals, coffee and tea, condiments and potatoes, and they were inversely related with palm oil, roots, tubers and plantains, refined cereals and fermented maize products. Among females, the DPs scores were characterized by frequent consumption of coffee and tea, poultry, whole-grain cereals, margarine, fish and alcoholic beverages and by low intakes of roots, tubers and plantains, fermented maize products, palm oil and refined cereals.
In the sensitivity analysis (Table S1), similar DP scores were obtained when self-reported T2DM cases were excluded. Among males, the explained variation in food group intake was 9.6%, and the explained variation in the FLI was 15.5%. In females, these figures were 6.3% and 8.9%, respectively.
4.4. Associations of FLI-Related Pattern Scores with T2DM
The associations with T2DM are shown per quintile and per one SD of the RRR-derived DP scores in Table 3. Among males, a higher adherence to the DP score was associated with increased odds of T2DM. As a trend, this was seen across the DP score quintiles, and these associations were significant per one score-SD increase. In the crude model, the odds of T2DM were 55% higher per 1 score-SD (95% CI: 1.30–1.86). The association attenuated after adjustment for the age and study site (Model 1: 1.34; 95% CI: 1.04–1.73) and remained stable after further adjustment for socioeconomic and lifestyle factors (Model 2: 1.45; 95% CI: 1.10–1.93). Among females, the associations were weaker. The RRR-derived DP score was positively associated with T2DM per one score-SD increase (crude model: 1.24; 95% CI: 1.07–1.44) but attenuated after adjustment for demographic, socioeconomic and lifestyle factors (Model 2: 1.16; 95% CI: 0.95–1.42).
In the sensitivity analysis (Table S2), we applied an HbA1c-based definition of T2DM as an outcome variable. This produced a stronger association between the RRR-derived DP scores and T2DM for both men and women. In the final model, the associations were OR: 1.54 (95% CI: 1.19–1.99) for men and OR: 1.27 (95% CI: 1.07–1.52) for women (Table S2).
4.5. RRR-Derived Dietary Patterns Related to the NAFLD Biomarkers
In addition to the FLI-related DP scores, we also extracted two DP scores related to the biomarkers of NAFLD (i.e., liver enzymes, blood lipids and CRP). Table 4 shows these DP scores and the corresponding factor loadings of the food groups. Among males, the DP score explained 11.2% of the total variation in the food group intake, and it was characterized by high intakes of whole-grain cereals, poultry, dairy products, coffee and tea, condiments, potatoes, margarine and olive oil. Again, this DP correlated with rare intakes of palm oil, roots, tubers and plantains and fermented maize products. In women, the RRR-derived DP score explained 12.1% of the total variation of the food intake, and it was correlated with frequent intakes of palm oil, vegetarian mixed dishes and fish. Interestingly, the DP score among females correlated inversely with food groups that loaded positively among females, namely whole-grain cereals, poultry, dairy products, coffee and tea, condiments, potatoes, margarine and olive oil.
The explained variations in the NAFLD biomarkers and the response weights for the RRR-derived DP scores are presented in Table 5. Among males, the DP score explained 2.8% of the total biomarker variation, with the major contribution to AST variation (7.4%) and a response weight of −0.58. This was followed by total cholesterol (4.5%) and LDL-cholesterol (4.2%), both showing positive response weights with the RRR-derived DP score. The relationships with AST were mainly driven by whole-grain cereals and coffee and tea (both r = −0.20). In addition, the positive response weights for the total cholesterol and LDL-cholesterol were mainly attributed to whole-grain cereals, potatoes and condiments (Table S3). In females, the explained total variation in the NAFLD biomarkers was 4.5%, with major contributions to AST at 12.3% (response weight: 0.58), triglycerides at 11.1% (response weight: 0.55) and HDL-cholesterol at 8.6% (response weight: −0.49) (Table 5). These relationships for AST and triglycerides were mainly driven by vegetarian mixed dishes and palm oil (both r = 0.20), and that for HDL-cholesterol was driven by coffee and tea and potatoes (Table S4). In our sensitivity analysis, after exclusion of participants with self-reported T2DM, we observed similar DP scores related to the NAFLD proxy markers as in the full dataset. Among males, the DP scores explained 11.0% of the total variation in the food group intake and 3.3% of the total variation in the biomarkers of NAFLD (Table S5). Among females, the DP scores explained 12.0% of the variation in food groups and 4.5% of the total variation in the biomarkers of NAFLD (Table S4).
4.6. Associations of NAFLD Biomarker-Related Pattern Scores with T2DM
Table 6 shows the odds of T2DM across quintiles and per one SD of the RRR-derived DP scores related to the NAFLD biomarkers, stratified by sex. Among males, the crude model showed a positive association with T2DM status (p for trend across quintiles = 0.0002; OR per 1 score-SD increase: 1.42; 95% CI: 1.20–1.68). However, this association attenuated to the null when adjusting for demographics, socioeconomics, lifestyle and anthropometry. Contrary to the males, there was no clear linear trend for the associations across quintiles with T2DM status (p for trend >0.05). Positive associations ranged from 8% increased odds of T2DM in the crude model (95% CI: 0.93–1.25) to 30% increased odds of T2DM in Model 3 (95% CI: 0.99–1.71) per one score-SD increase.
Again, we conducted sensitivity analysis using an HbA1c definition as the outcome variable (Table S6). The corresponding associations of the biomarker-related DP scores were weaker than those in the primary analysis. In the final Model 2 results, the associations were OR: 1.01 (95% Cl: 0.78–1.32) among men and OR: 1.27 (95% CI: 0.99–1.62) among women (Table S6).
5. Discussion
The relationships between NAFLD and T2DM with respect to their dietary risk factors informed the derivation of DP scores by applying the RRR method. We used two sets of proxy markers for NAFLD (1: FLI; 2: serum liver enzymes, blood lipids and CRP) as response variables and the habitual consumption of 30 food groups as predictor variables. The FLI-related DPs were characterized by frequent intakes of modernized foods and were positively associated with the odds of T2DM. This association was stronger in males than in females. The RRR-derived DP scores using CRP, liver enzymes and blood lipids as the response variables showed inconsistent relationships with food intake between males and females and implausible associations with T2DM.
5.1. Dietary Patterns and Proxy Markers of NALFD and T2DM
In our first approach, the FLI was used as a proxy marker for NALFD and as the response variable in the RRR analysis. The respective FLI-related DPs were similar between males and females and were characterized by frequent intakes of poultry, whole-grain cereals, coffee and tea, condiments, potatoes, alcoholic beverages, margarine and fish, while fermented maize products, refines cereals, roots, tubers and plantains and palm oil showed inverse correlations with the FLI. These DPs correspond to the nutrition transition observed for SSA populations when experiencing rapid economic growth, accelerated urbanization and migration to Europe [10,26].
Our findings partly accord with recent evidence from a meta-analysis on dietary risk factors of NAFLD [27]. Red meat and sugar-sweetened beverages are robustly associated with a higher risk of NAFLD, while nuts are consistently associated with lower NAFLD risk. This meta-analysis demonstrates that many foods have neutral relationships with NAFLD, such as dairy, fish, rice, noodles and whole-grain products [27]. In our study, we also identified positive factor loadings for condiments, meat (poultry) and coffee and tea, which showed positive correlations with the FLI. For coffee and tea, this appears to contradict the observed lower risk of elevated liver enzymes associated with frequent coffee and tea consumption [28]. The antioxidant properties of these beverages may exhibit a suppressive effect of hyperglycemia by enhancing insulin sensitivity and subsequently improving a fatty liver [29]. However, coffee and tea are often sweetened with sugar and condensed milk, which happens to be underreported during dietary assessments. Thus, the high glycemic load of sweetened beverages and condiments may explain the positive correlations with the FLI in our study. Reassuringly, in a large Chinese population, an RRR-derived DP related to inflammatory markers is also characterized by high intakes of sugar-containing foods and positively associated with NAFLD [30]. In fact, low glycemic index diets are currently promoted for the amelioration of NAFLD [31].
While many studies confirm the benefits of poultry and fish as part of a healthy diet [32,33], this may not apply uniformly. Review data and recent findings from the multi-ethnic cohort show that poultry intake can be associated with higher risk of NAFLD [27,34]. The underlying mechanisms remain to be investigated, as most studies did not examine poultry separately but focused on overall meat consumption. Additionally, dietary patterns rich in fish intake are known to prevent NAFLD, mainly due to the beneficial effects of omega-3 poly-unsaturated fatty acids on health [33]. Yet, the increased contamination of sea fish with environmental toxins might cause harmful effects to the human endocrine system and promote NAFLD [35], which could contribute to our findings. In light of this, our study reaffirms the current dietary recommendations for the prevention of NAFLD and T2DM, including foods with low glycemic indexes that are rich in poly-unsaturated fatty acids [27].
Interestingly, we observed stronger relationships for the FLI-related DPs among males compared with females, despite the fact that the characteristic food groups and the explained variation in the FLI were similar between sexes. This may reflect the sex-specific contributions of NAFLD to the etiology of T2DM. Indeed, recent findings from the RODAM study confirmed the increased odds of T2DM among participants with an FLI ≥60 compared with an FLI <60. This association is stronger in men (OR: 2.43; 95% CI: 1.73–3.41) than in women (OR: 2.02; 95% CI: 1.52–2.69) [36]. Further studies found a positive association between FLI and metabolic-associated fatty liver diseases (MAFLD) [36]. This observation was attributed to several factors, including diet-related changes in the composition of the gut microbiota, which in turn may influence the development of MAFLD. Male adults from Ghana and Kenya present a lower prevalence of an increased FLI and higher prevalence of T2DM than females [37,38]. The explanation for the observed sex-specific relationships between the FLI—as a proxy marker of NAFLD—and T2DM may lay in the FLI components. It is well-established that measures of central adiposity, such as the waist circumference and waist-to-hip ratio, predict the risk of cardiometabolic diseases among SSA populations better than the measure of general obesity (BMI), particularly among females [39,40]. The fact that females of African ancestry tend to have more subcutaneous and less visceral fat at any given level of the BMI than men supports the notion that female individuals can metabolically tolerate a certain degree of obesity as opposed to males [41]. In addition, the biomarker components of the FLI may contribute to the observed sex differences, including serum triglycerides and HDL-cholesterol [37].
Additionally, in an urban Ghanaian case control study, Frank et al. derived a DP using RRR that has been positively related with serum triglycerides and inversely related with adiponectin. This pattern is associated with higher odds of T2DM and shares similar features with the RRR-derived DPs in the present analysis, being characterized by high consumption of starchy foods and low intake of fruits and vegetables [16]. The combination of condiments, margarine and higher glycemic index foods such potatoes may contribute to the observed associations with T2DM [42,43]. These foods characterize modernized diets, which are known to confer an unhealthy body composition and increased risk of T2DM [44]. In addition, our findings are partly corroborated by recent evidence from a meta-analysis of 16 cohort studies. In this study, the RRR-derived DPs were characterized by high intakes of refined grains, sugar-sweetened soft drinks and processed meat and were significantly associated with T2DM risk [45]. In our study, the DPs were also characterized by sugary foods such as condiments and sweetened hot beverages and were positively associated with T2DM. Notably, previous exploratory approaches in this Ghanaian study population identified three DPs using PCA that had either no or a negative association with T2DM [46]. A rice, pasta, meat and fish DP showed increased odds of T2DM by one SD-score increase but overlapped with the current RRR-derived DPs only regarding condiments and fish. The strongest similarities of the RRR-derived DPs were seen for the previously identified mixed DP, which was not associated with T2DM. At the same time, the strongest dissimilarities were observed for the roots, tubers and plantains DP, which tended to be positively associated with T2DM. These findings may reflect the multiple pathways through which dietary factors contribute to the etiology of T2DM, with NAFLD being only one possible link.
In our second approach, we used a selection of proxy markers for NAFLD as response variables for the RRR analysis, comprising liver enzymes, blood lipids and CRP. In this analysis, we observed inconsistent DP scores between men and women and biologically implausible associations with T2DM. Presumably, this selection of biomarkers may not be specific enough to operationalize NAFLD, as they also reflect other metabolic pathways to T2DM, including chronic inflammation (CRP) and dyslipidemia (blood lipids) [47], which may explain the unexpected results in this analysis. In addition, unmeasured confounding might have distorted our findings with regard to biomarker-related DPs.
5.2. Strengths and Limitations
When interpreting the findings of our study, some limitations need to be acknowledged. First, the definition of NAFLD by the use of the FLI as a proxy marker without liver biopsy or imaging may shield detailed information on the severity classification of a fatty liver [48]. However, liver biopsy is invasive, associated with discomfort and often technically unfeasible in large-scale epidemiological studies. In contrast, the FLI has proven to be a useful predictor for a fatty liver, since it shows high concordance with the imaging and histological criteria for NAFLD [49]. At the same time, we acknowledge that validation of the FLI among sub-Saharan African populations is still pending [36]. As data on viral hepatitis status was not available in our study, misclassification of NAFLD cannot not fully be excluded. However, the study provides unique insights into diet-disease pathways for T2DM in a large sample of adults from sub-Saharan Africa. The method for DP construction accounts for the potential links between NAFLD and T2DM through projecting the response variable of the FLI on the food groups as predictor variables. We are cognizant that this approach has been criticized for its choice of proxy measures, which could simply reflect the early disease stages of T2DM [50]. In fact, previous data from Ghanaian adults suggest that increased serum triglycerides serve as an important risk factor for T2DM [16,51], while the BMI does not [1,39]. For the associations of FLI-related DP scores with T2DM, we did not adjust for anthropometric components of the FLI in our final regression model, which might have caused some residual confounding and overestimation of the observed DP–T2DM relationships, as the FLI formula includes the BMI and waist circumference. Like in all cross-sectional studies, reverse causation cannot be ruled out in our study either, particularly for the bidirectional relationship between NAFLD and T2DM, given the fact that some study participants with T2DM had received treatment, including lifestyle counseling. Still, the DP scores remained when individuals with self-reported T2DM were excluded. Finally, we have addressed potential confounding by our comprehensive sets of adjustment variables. However, unmeasured confounding, such as a family history of diabetes, cannot be ruled out.
6. Conclusions
We identified DPs that showed adherence to modernized foods and a positive relation with the FLI, a proxy marker for fatty liver disease. These DPs were positively associated with the odds of T2DM, particularly in males. These findings give credence to the concept that modernized dietary practices among Ghanaian populations under transition may underly metabolic pathways to NAFLD and T2DM.
Supplementary Materials
The following are available online at https://www.mdpi.com/article/10.3390/nu13113679/s1, Figure S1. present the flow diagram of excluded participants because of missing or implausible data resulting in a simple size of 3687 participants for the characterization of dietary pattern. Table S1: Sensitivity analysis excluding study participants with self-reported diabetes. Percentage of explained variation in food intake and factor loadings of RRR-derived dietary pattern scores related with Fatty Liver Index among men and women. Table S2. Associations of the FLI-related and RRR-derived dietary pattern scores with type 2 diabetes among males and females 1. Table S3. Pearson correlations between the biomarker-related dietary pattern score, food intake frequencies, and NAFLD biomarkers (log-transformed) among men (n = 1366). Table S4. Pearson correlations between the biomarker-related dietary pattern score, food intake frequencies, and NAFLD biomarkers (log-transformed) among women (n = 2321). Table S5. Sensitivity analysis excluding study participants with self-reported diabetes. Percentage of explained variation in food intake and factor loadings of RRR-derived dietary pattern scores, and explained variation in NAFLD biomarkers and response weights among men and women. Table S6. Associations of the biomarker-related and RRR-derived dietary pattern scores with type 2 diabetes among males and females 1.
Author Contributions
T.B.O. performed the statistical analysis and data interpretation and wrote the manuscript draft. I.D. and T.B.O. conceived and designed the study. The RODAM consortium members were responsible for on-site recruitment and data collection. A.-M.v.D., S.D., F.P.C., E.B., K.A.C.M., S.B., M.B.S., C.A., M.N. and A.G.H. read, edited, and contributed to the manuscript revision. I.D. supervised the study. All authors have read and agreed to the published version of the manuscript.
Funding
The RODAM study was supported by the European Commission under the Framework Programme (Grant Number: 278901). K.A.C.M. is supported by the Intramural Research Program of the National Institutes of Health in the Center for Research on Genomics and Global Health (CRGGH). A. G. Holleboom is supported by the Amsterdam UMC Fellowship, Health Holland TKI public-private partnership grants and research grants from Gilead and Novo Nordisk. Tracy Bonsu Osei is supported by the German Academic Exchange Service (DAAD), and Ina Danquah is supported by the Robert Bosch Foundation (RBS) (0100035-002).
Institutional Review Board Statement
The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Ethics committee of Ghana (School of Medical Sciences, Komfo Anokye Teaching Hospital, Committee on Human Research, Publication and Ethical Review Board- CHRPE/AP/200/12), the Netherlands (Institutional Review Board of AMC, University of Amsterdam- E2-236), the UK (London School of Hygiene and Tropical Medicine Research Ethics Committee- 6208) and Germany (Ethics Committee of Charite-Universitatsmedizin Berlin-EA1/307/12).
Informed Consent Statement
Informed consent was obtained from all subjects involved in the study.
Data Availability Statement
Data will be available upon request from the corresponding author.
Acknowledgments
The authors are very grateful to the advisory board members for their valuable support in shaping the methods, to the research assistants, interviewers and other staff of the five research locations who took part in gathering the data and, most of all, to the Ghanaian volunteers who participated in this project. We gratefully acknowledge J. van Straalen from the Department of Clinical Chemistry of the Academic Medical Centre (Amsterdam, The Netherlands) for his valuable support with standardization of the laboratory procedures and the Academic Medical Center (AMC) Biobank for support in biobank management and storage of the collected samples.
Conflicts of Interest
The authors declare no conflict of interest.
References
- Agyemang, C.; Meeks, K.; Beune, E.; Owusu-Dabo, E.; Mockenhaupt, F.P.; Addo, J.; de Graft Aikins, A.; Bahendeka, S.; Danquah, I.; Schulze, M.B.; et al. Obesity and type 2 diabetes in sub-Saharan Africans–Is the burden in today’s Africa similar to African mi-grants in Europe? The RODAM study. BMC Med. 2016, 14, 166. [Google Scholar] [CrossRef] [Green Version]
- Onyekwere, C.A.; Ogbera, A.O.; Balogun, B.O. Non-alcoholic fatty liver disease and the metabolic syndrome in an urban hos-pital serving an African community. Ann. Hepatol. 2016, 10, 119–124. [Google Scholar] [CrossRef]
- Ahmed, M.H.; Noor, S.K.; Bushara, S.O.; Husain, N.E.; Elmadhoun, W.M.; Ginawi, I.A.; Osman, M.M.; Mahmoud, A.O.; Almobarak, A.O. Non-alcoholic fatty liver disease in Africa and Middle East: An attempt to predict the present and future impli-cations on the healthcare system. Gastroenterol. Res. 2017, 10, 271. [Google Scholar] [CrossRef] [Green Version]
- Campos-Murguía, A.; Ruiz-Margáin, A.; González-Regueiro, J.; Macías-Rodríguez, R.U. Clinical assessment and management of liver fibrosis in non-alcoholic fatty liver disease. World J. Gastroenterol. 2020, 26, 5919–5943. [Google Scholar] [CrossRef]
- Lallukka, S.; Yki-Järvinen, H. Non-alcoholic fatty liver disease and risk of type 2 diabetes. Best Pract. Res. Clin. Endocrinol. Metab. 2016, 30, 385–395. [Google Scholar] [CrossRef] [Green Version]
- Machado, M.V.; Cortez-Pinto, H. Non-alcoholic fatty liver disease: What the clinician needs to know. World J. Gas-Troenterol. 2014, 20, 12956. [Google Scholar] [CrossRef]
- Idalsoaga, F.; Kulkarni, A.V.; Mousa, O.Y.; Arrese, M.; Arab, J.P. Non-alcoholic Fatty Liver Disease and Alcohol-Related Liver Disease: Two Intertwined Entities. Front. Med. 2020, 7, 7. [Google Scholar] [CrossRef]
- Lee, B.-W.; Lee, Y.H.; Park, C.Y.; Rhee, E.J.; Lee, W.Y.; Kim, N.H.; Choi, M.K.; Park, K.-G.; Choi, Y.-K.; Cha, B.S.; et al. Non-alcoholic fatty liver disease in patients with type 2 diabetes mellitus: A position statement of the fatty liver re-search group of the korean diabetes association. Diabetes Metab. J. 2020, 44, 382–401. [Google Scholar] [CrossRef]
- Ciardullo, S.; Monti, T.; Perseghin, G. High Prevalence of Advanced Liver Fibrosis Assessed by Transient Elastography among U.S. Adults with Type 2 Diabetes. Diabetes Care 2021, 44, 519–525. [Google Scholar] [CrossRef]
- Galbete, C.; Nicolaou, M.; Meeks, K.; Aikins, A.D.-G.; Addo, J.; Amoah, S.K.; Smeeth, L.; Owusu-Dabo, E.; Klipstein-Grobusch, K.; Bahendeka, S.; et al. Food consumption, nutrient intake, and dietary patterns in Ghanaian migrants in Europe and their compatriots in Ghana. Food Nutr. Res. 2017, 61, 1341809. [Google Scholar] [CrossRef] [Green Version]
- Erridge, C. Diet, commensals and the intestine as sources of pathogen-associated molecular patterns in atherosclerosis, type 2 diabetes and non-alcoholic fatty liver disease. Atherosclerosis 2011, 216, 1–6. [Google Scholar] [CrossRef]
- Jegatheesan, P.; De Bandt, J. Fructose and NAFLD: The Multifaceted Aspects of Fructose Metabolism. Nutrients 2017, 9, 230. [Google Scholar] [CrossRef] [Green Version]
- Yki-Järvinen, H.; Luukkonen, P.K.; Hodson, L.; Moore, J.B. Dietary carbohydrates and fats in nonalcoholic fatty liver disease. Nat. Rev. Gastroenterol. Hepatol. 2021, 1–17. [Google Scholar] [CrossRef]
- Della Pepa, G.; Vetrani, C.; Lombardi, G.; Bozzetto, L.; Annuzzi, G.; Rivellese, A.A. Isocaloric Dietary Changes and Non-Alcoholic Fatty Liver Disease in High Cardiometabolic Risk Individuals. Nutrients 2017, 9, 1065. [Google Scholar] [CrossRef] [Green Version]
- Tana, C.; Ballestri, S.; Ricci, F.; Di Vincenzo, A.; Ticinesi, A.; Gallina, S.; Giamberardino, M.A.; Cipollone, F.; Sutton, R.; Vettor, R.; et al. Cardiovascular Risk in Non-Alcoholic Fatty Liver Disease: Mechanisms and Therapeutic Implications. Int. J. Environ. Res. Public Health 2019, 16, 3104. [Google Scholar] [CrossRef] [Green Version]
- Frank, L.K.; Jannasch, F.; Kröger, J.; Bedu-Addo, G.; Mockenhaupt, F.P.; Schulze, M.B.; Danquah, I. A Dietary Pattern Derived by Reduced Rank Regression is Associated with Type 2 Diabetes in An Urban Ghanaian Population. Nutrients 2015, 7, 5497–5514. [Google Scholar] [CrossRef] [Green Version]
- Livingstone, K.M.; McNaughton, S.A. Association between diet quality, dietary patterns and cardiometabolic health in Aus-tralian adults: A cross-sectional study. Nutr. J. 2018, 17, 19. [Google Scholar] [CrossRef] [Green Version]
- Vermeulen, E.; Stronks, K.; Visser, M.; Brouwer, I.A.; Schene, A.H.; Mocking, R.J.T.; Colpo, M.; Bandinelli, S.; Ferrucci, L.; Nicolaou, M. The association between dietary patterns derived by reduced rank regression and depressive symptoms over time: The Invecchiare in Chianti (InCHIANTI) study. Br. J. Nutr. 2016, 115, 2145–2153. [Google Scholar] [CrossRef] [Green Version]
- Byrne, C.D.; Targher, G. NAFLD: A multisystem disease. J. Hepatol. 2015, 62, S47–S64. [Google Scholar] [CrossRef] [Green Version]
- Drescher, H.K.; Weiskirchen, S.; Weiskirchen, R. Current Status in Testing for Nonalcoholic Fatty Liver Disease (NAFLD) and Nonalcoholic Steatohepatitis (NASH). Cells 2019, 8, 845. [Google Scholar] [CrossRef] [Green Version]
- Younossi, Z.M.; Loomba, R.; Anstee, Q.M.; Rinella, M.E.; Bugianesi, E.; Marchesini, G.; Neuschwander-Tetri, B.A.; Serfaty, L.; Negro, F.; Caldwell, S.H.; et al. Diagnostic modalities for nonalcoholic fatty liver disease, nonalcoholic steatohepatitis, and associated fibrosis. Hepatology 2018, 68, 349–360. [Google Scholar] [CrossRef] [Green Version]
- Agyemang, C.; Beune, E.; Meeks, K.; Owusu-Dabo, E.; Agyei-Baffour, P.; Aikins, A.D.-G.; Dodoo, F.; Smeeth, L.; Addo, J.; Mockenhaupt, F.; et al. Rationale and cross-sectional study design of the Research on Obesity and type 2 Diabetes among African Migrants: The RODAM study. BMJ Open 2014, 4, e004877. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Bedogni, G.; Bellentani, S.; Miglioli, L.; Masutti, F.; Passalacqua, M.; Castiglione, A.; Tiribelli, C. The Fatty Liver Index: A simple and accurate predictor of hepatic steatosis in the general population. BMC Gastroenterol. 2006, 6, 33. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Haskell, W.L.; Lee, I.M.; Pate, R.R.; Powell, K.E.; Blair, S.N.; Franklin, B.A.; Macera, C.A.; Heath, G.W.; Thompson, P.T.; Bauman, A. Physical activity and public health: Updated recommendation for adults from the American College of Sports Medicine and the American Heart Association. Circulation 2007, 116, 1081. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Hoffmann, K.; Schulze, M.B.; Schienkiewitz, A.; Nöthlings, U.; Boeing, H. Application of a New Statistical Method to Derive Dietary Patterns in Nutritional Epidemiology. Am. J. Epidemiol. 2004, 159, 935–944. [Google Scholar] [CrossRef]
- Addo, J.; Agyemang, C.; de-Graft Aikins, A.; Beune, E.; Schulze, M.B.; Danquah, I.; Galbete, C.; Nicolaou, M.; Meeks, K.; Klipstein-Grobusch, K.; et al. Association between socioeconomic position and the prevalence of type 2 diabetes in Ghanaians in different geo-graphic locations: The RODAM study. J. Epidemiol. Community Health 2017, 71, 633–639. [Google Scholar] [CrossRef] [Green Version]
- He, K.; Li, Y.; Guo, X.; Zhong, L.; Tang, S. Food groups and the likelihood of non-alcoholic fatty liver disease: A systematic review and meta-analysis. Br. J. Nutr. 2020, 124, 1–13. [Google Scholar] [CrossRef]
- Birerdinc, A.; Stepanova, M.; Pawloski, L.; Younossi, Z.M. Caffeine is protective in patients with non-alcoholic fatty liver disease. Aliment. Pharmacol. Ther. 2012, 35, 76–82. [Google Scholar] [CrossRef]
- Molloy, J.W.; Calcagno, C.J.; Williams, C.D.; Jones, F.J.; Torres, D.M.; Harrison, S.A. Association of coffee and caffeine consumption with fatty liver disease, nonalcoholic steatohepatitis, and degree of hepatic fibrosis. Hepatology 2012, 55, 429–436. [Google Scholar] [CrossRef]
- Xia, Y.; Zhang, Q.; Liu, L.; Meng, G.; Wu, H.; Bao, X.; Gu, Y.; Sun, S.; Wang, X.; Zhou, M.; et al. Intermediary effect of inflammation on the association between dietary patterns and non-alcoholic fatty liver disease. Nutrients 2020, 71, 110562. [Google Scholar] [CrossRef]
- Franco, I.; Bianco, A.; Mirizzi, A.; Campanella, A.; Bonfiglio, C.; Sorino, P.; Notarnicola, M.; Tutino, V.; Cozzolongo, R.; Giannuzzi, V.; et al. Physical Activity and Low Glycemic Index Mediterranean Diet: Main and Modification Effects on NAFLD Score. Results from a Randomized Clinical Trial. Nutrients 2020, 13, 66. [Google Scholar] [CrossRef] [PubMed]
- Abenavoli, L.; Boccuto, L.; Federico, A.; Dallio, M.; Loguercio, C.; Di Renzo, L.; De Lorenzo, A. Diet and Non-Alcoholic Fatty Liver Disease: The Mediterranean Way. Int. J. Environ. Res. Public Health 2019, 16, 3011. [Google Scholar] [CrossRef] [Green Version]
- Giraldi, L.; Miele, L.; Aleksovska, K.; Manca, F.; Leoncini, E.; Biolato, M.; Pirro, M.A.; Marrone, G.; Cefalo, C.; Racco, S.; et al. Mediterranean diet and the prevention of non-alcoholic fatty liver disease: Results from a case-control study. Eur. Rev. Med. Pharmacol. Sci. 2020, 24, 7391–7398. [Google Scholar] [PubMed]
- Noureddin, M.; Jones, C.; Alkhouri, N.; Gomez, E.V.; Dieterich, D.T.; Rinella, M.E.; Therapondos, G.; Girgrah, N.; Mantry, P.S.; Sussman, N.L.; et al. Screening for Nonalcoholic Fatty Liver Disease in Persons with Type 2 Diabetes in the United States Is Cost-effective: A Comprehensive Cost-Utility Analysis. Gastroenterology 2020, 159, 1985–1987.e4. [Google Scholar] [CrossRef] [PubMed]
- Treviño, L.S.; Katz, T.A. Endocrine Disruptors and Developmental Origins of Nonalcoholic Fatty Liver Disease. Endocrinology 2017, 159, 20–31. [Google Scholar] [CrossRef] [PubMed]
- Van Dijk, A.-M.; Dingerink, S.; Chilunga, F.P.; Meeks, K.A.C.; Bahendeka, S.; Schulze, M.B.; Danquah, I.; Osei, T.B.; Serné, E.; Agyemang, C.; et al. Metabolic-associated Fatty Liver Disease as Assessed by the Fatty Liver Index Among Migrant and Non-migrant Ghanaian Populations. J. Clin. Transl. Hepatol. 2021, 9, 494–502. [Google Scholar] [PubMed]
- Christensen, D.L.; Faurholt-Jepsen, D.; Birkegaard, L.; Mwaniki, D.L.; Boit, M.K.; Kilonzo, B.; Brage, S.; Friis, H.; Tetens, I.; Borch-Johnsen, K.; et al. Cardiovascular risk factors in rural Kenyans are associated with differential age gradients, but not modified by sex or ethnicity. Ann. Hum. Biol. 2015, 43, 42–49. [Google Scholar] [CrossRef] [Green Version]
- Paquette, M.; Gauthier, D.; Chamberland, A.; Prat, A.; Rolfe, E.D.L.; Rasmussen, J.J.; Kaduka, L.; Seidah, N.G.; Bernard, S.; Christensen, D.L.; et al. Circulating PCSK9 is associated with liver biomarkers and hepatic steatosis. Clin. Biochem. 2020, 77, 20–25. [Google Scholar] [CrossRef]
- Danquah, I.; Othmer, T.; Frank, L.K.; Bedu-Addo, G.; Schulze, M.B.; Mockenhaupt, F.P. The TCF7L2 rs7903146 (T) allele is associated with type 2 diabetes in urban Ghana: A hospital-based case–control study. BMC Med. Genet. 2013, 14, 96. [Google Scholar] [CrossRef] [Green Version]
- Darko, S.N.; Meeks, K.A.; Owiredu, W.K.; Laing, E.F.; Boateng, D.; Beune, E.; Addo, J.; Aikins, A.D.-G.; Bahendeka, S.; Mockenhaupt, F.; et al. Anthropometric indices and their cut-off points in relation to type 2 diabetes among Ghanaian migrants and non-migrants: The RODAM study. Diabetes Res. Clin. Pr. 2021, 173, 108687. [Google Scholar] [CrossRef] [PubMed]
- Goedecke, J.H.; Olsson, T. Pathogenesis of type 2 diabetes risk in black Africans: A South African perspective. J. Intern. Med. 2020, 288, 284–294. [Google Scholar] [CrossRef] [Green Version]
- Aljefree, N.; Ahmed, F. Association between dietary pattern and risk of cardiovascular disease among adults in the Middle East and North Africa region: A systematic review. Food Nutr. Res. 2015, 59, 27486. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Cempaka, A.R.; Tseng, S.-H.; Yuan, K.-C.; Bai, C.-H.; Tinkov, A.A.; Skalny, A.V.; Chang, J.-S. Dysregulated Iron Metabolism-Associated Dietary Pattern Predicts an Altered Body Composition and Metabolic Syndrome. Nutrients 2019, 11, 2733. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Fung, T.T.; Schulze, M.B.; Manson, J.E.; Willett, W.C.; Hu, F.B. Dietary Patterns, Meat Intake, and the Risk of Type 2 Diabetes in Women. Arch. Intern. Med. 2004, 164, 2235–2240. [Google Scholar] [CrossRef]
- Jannasch, F.; Kröger, J.; Schulze, M.B. Dietary Patterns and Type 2 Diabetes: A Systematic Literature Review and Me-ta-Analysis of Prospective Studies. J. Nutr. 2017, 147, 1174–1182. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Galbete, C.; Nicolaou, M.; Meeks, K.; Klipstein-Grobusch, K.; Aikins, A.D.-G.; Addo, J.; Amoah, S.K.; Smeeth, L.; Owusu-Dabo, E.; Spranger, J.; et al. Dietary patterns and type 2 diabetes among Ghanaian migrants in Europe and their compatriots in Ghana: The RODAM study. Nutr. Diabetes 2018, 8, 25. [Google Scholar] [CrossRef] [PubMed]
- Feinglos, M.N.; Bethel, M.A. Type 2 Diabetes Mellitus: An Evidence-Based Approach to Practical Management; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2008. [Google Scholar]
- Von Eckardstein, A. Laboratory diagnostics of non-alcoholic fatty liver disease. Laboratoriums Medizin 2015, 38, s1. [Google Scholar] [CrossRef] [Green Version]
- Olubamwo, O.O.; Virtanen, J.K.; Pihlajamaki, J.; Mantyselka, P.; Tuomainen, T.P. Fatty liver index as a predictor of increased risk of cardiometabolic disease: Finding from the Kuopio Is-chaemic Heart Disease Risk Factor Study Cohort. BMJ Open 2019, 9, e031420. [Google Scholar] [CrossRef] [Green Version]
- Weikert, C.; Schulze, M.B. Evaluating dietary patterns: The role of reduced rank regression. Curr. Opin. Clin. Nutr. Metab. Care 2016, 19, 341–346. [Google Scholar] [CrossRef] [PubMed]
- Danquah, I.; Bedu-Addo, G.; Terpe, K.-J.; Micah, F.; Amoako, Y.A.; Awuku, Y.A.; Dietz, E.; Van Der Giet, M.; Spranger, J.; Mockenhaupt, F.P. Diabetes mellitus type 2 in urban Ghana: Characteristics and associated factors. BMC Public Health 2012, 12, 210. [Google Scholar] [CrossRef] [Green Version]
Figure 1.
Mean intakes and standard deviations (g/day) of (a) 16 food groups with a mean intake of ˃50 g/day and (b) 14 food groups with a mean intake of ≤50 g/day.
Figure 1.
Mean intakes and standard deviations (g/day) of (a) 16 food groups with a mean intake of ˃50 g/day and (b) 14 food groups with a mean intake of ≤50 g/day.
Table 1.
Characteristics and biomarkers of the RODAM study population by sex and site 1.
| Characteristics | Total (n = 3687) | Men (n = 1366) | Women (n = 2321) | Rural Ghana (n = 820) | Urban Ghana (n = 1358) | Amsterdam (n = 707) | Berlin (n = 451) | London (n = 351) |
|---|---|---|---|---|---|---|---|---|
| Sex (female%) | 63.0 | - | - | 61.6 | 72.2 | 60.1 | 45.0 | 59.0 |
| Age (years) | 46.1 ± 11.1 | 46.9 ± 11.3 | 45.6 ± 10.9 | 46.7 ± 12.6 | 45.3 ± 11.4 | 46.6 ± 8.5 | 45.2 ± 10.4 | 47.9 ± 10.9 |
| Education | ||||||||
| Never or elementary% | 37.6 | 21.9 | 46.9 | 56.8 | 43.9 | 35.5 | 9.3 | 9.1 |
| Low | 37.7 | 41.6 | 35.5 | 31.6 | 38.9 | 37.8 | 50.1 | 31.6 |
| Intermediate | 16.2 | 22.5 | 12.5 | 7.9 | 12.5 | 21.8 | 26.6 | 24.8 |
| Higher vocational | 8.5 | 14.1 | 5.2 | 3.7 | 4.71 | 5.0 | 14.0 | 34.5 |
| Length of stay (years) | - | - | - | - | - | 16.4 ± 8.1 | 17.0 ± 10.9 | 17.2 ± 11.0 |
| Body mass index (kg/m2) | 26.7 ± 5.5 | 24.8 ± 4.4 | 27.8 ± 5.7 | 22.7 ± 4.3 | 26.9 ± 5.4 | 28.9 ± 5 | 27.6 ± 4.8 | 29.4 ± 4.8 |
| Waist circumference (cm) | 89.5 ± 12.5 | 87 ± 12.1 | 91 ± 12.5 | 81.2 ± 10.9 | 89.4 ± 11.8 | 94.6 ± 11.6 | 92.2 ± 11.5 | 95.4 ± 11.3 |
| Smoking (current or former%) | 9.3 | 19.6 | 3.2 | 7.9 | 6.9 | 11.6 | 18.4 | 5.4 |
| Physical activities (MET-min/day) | 72 (14–168) | 96 (28–196) | 56 (10–156) | 90 (36–161) | 60 (6–156) | 88.7 (26–258) | 72 (12–198) | 28 (5–112) |
| Total energy intake (kcal/day) | 2533 ± 837 | 2628 ± 827 | 2477 ± 817 | 2594 ± 828 | 2298 ± 661 | 2478 ± 854 | 2929 ± 944 | 2898 ± 953 |
| Carbohydrate intake (energy%) | 53.0 ± 9.1 | 52.2 ± 9.5 | 53.5 ± 8.9 | 56.4.5 ± 8.3 | 54.5 ± 8.1 | 50.5 ± 8.3 | 48.5 ± 10.9 | 50.2 ± 9.6 |
| Fat intake (energy%) | 32.2 ± 8.2 | 32 ± 8.6 | 32.3 ± 8 | 31.4 ± 7.3 | 31.6 ± 7.2 | 32.1 ± 8.3 | 33.7 ± 10.6 | 34.1 ± 9.6 |
| Protein intake (energy%) | 13.8 ± 2.9 | 13.9 ± 3.1 | 13.8 ± 2.9 | 11.5 ± 2.2 | 13.6 ± 2.4 | 15.8 ± 2.7 | 14.8 ± 3.1 | 15.1 ± 2.9 |
| Alcohol intake (g/day) | 0 (0–0.1) | 0 (0–0.3) | 0 (0–0.1) | 0 (0–0.1) | 0 (0–0.1) | 0.1 (0–0.4) | 0.1 (0–0.6) | 0 (0–0.1) |
| AST U/L | 32.3 (26.6–39.8) | 35.1 (29.1–43.1) | 30.6 (25.3–37.7) | 36.1 (30.4–43.1) | 34.4 (28.7–41.5) | 26.1 (22.4–30.8) | 28.9 (24.7–34.9) | 34.1 (28.0–43.1) |
| ALT U/L | 19.2 (14.9–25.7) | 23.0 (17.4–31.2) | 17.6 (13.9–22.7) | 19.2 (15–24.9) | 19.3 (15–25.8) | 17.4 (13.7–23.0) | 19.9 (14.8–26.9) | 22.5 (18.3–30.3) |
| GGT U/L | 30.8 (23.2–43.1) | 37.4 (27.3–52.6) | 27.9 (21.7–37.3) | 29.5 (22.3–42.2) | 31.4 (23.9–42.9) | 30.2 (22.8–42.0) | 32.9 (24.7–46.1) | 30.6 (22.7–43.9) |
| CRP mg/L | 0.7 (0.2–2.5) | 0.5 (0.1–1.5) | 0.9 (0.2- 3.2) | 0.7 (0.2–2.6) | 0.9 (0.2–3.1) | 0.8 (0.2–2.3) | 0.5 (0.2–1.9) | 0.8 (0.2–2.3) |
| Total cholesterol (mmol/L) | 5.0 ± 1.1 | 4.9 ± 1.1 | 5.1 ± 1.1 | 4.6 ± 1.1 | 5.2 ± 1.2 | 5.0 ± 1.1 | 5.1 ± 1.1 | 5.0 ± 1.0 |
| LDL-cholesterol (mmol/L) | 3.2 ± 1.0 | 3.1 ± 1.0 | 3.2 ± 1.0 | 2.8 ± 1.0 | 3.4 ± 1.0 | 3.2 ± 0.9 | 3.2 ± 1.0 | 3.3 ± 0.9 |
| HDL-cholesterol (mmol/L) | 1.3 ± 0.4 | 1.3 ± 0.4 | 1.3 ± 0.4 | 1.2 ± 0.4 | 1.3 ± 0.3 | 1.4 ± 0.3 | 1.5 ± 0.4 | 1.3 ± 0.3 |
| Triglycerides (mmol/L) | 0.9 (0.7–1.2) | 1.0 (0.7–1.3) | 0.9 (0.7–1.2) | 1.0 (0.8- 1.3) | 1.0 (0.8–1.3) | 0.8 (0.6–1.0) | 0.9 (0.6–1.1) | 0.8 (0.6–1.1) |
| Fatty Liver Index | 2.6 ± 6.3 | 2.0 ± 5.4 | 2.9 ± 6.8 | 1.0 ± 3.2 | 2.8 ± 6.7 | 3.3 ± 7.2 | 3.1 ± 7.3 | 3.1 ± 6.0 |
1 Continuous variables are expressed as means ± SDs or medians (IQR). Categorical variables are expressed as percentages (%). AST: aspartate aminotransferase; ALT: alanine aminotransferase; CRP: C-reactive protein; GGT: gamma-glutamyl transferase; HDL: high-density lipoprotein; LDL: low-density lipoprotein.
Table 2.
Percentage of explained variation in food intake and factor loadings of 30 food groups of the RRR-derived dietary pattern scores related to the FLI for males and females 1.
Table 2.
Percentage of explained variation in food intake and factor loadings of 30 food groups of the RRR-derived dietary pattern scores related to the FLI for males and females 1.
| Food Group | Men (n = 1366) | Women (n = 2321) | ||
|---|---|---|---|---|
| Explained Variation (%) | Factor Loading | Explained Variation (%) | Factor Loading | |
| Poultry | 30.8 | 0.32 | 16.6 | 0.29 |
| Whole-grain cereals | 21.5 | 0.27 | 15.8 | 0.29 |
| Coffee and tea | 19.9 | 0.26 | 21.0 | 0.33 |
| Condiments | 18.3 | 0.25 | 15.7 | 0.28 |
| Potatoes | 17.8 | 0.25 | 2.1 | 0.10 |
| Alcoholic beverages | 10.3 | 0.19 | 7.5 | 0.20 |
| Margarine | 10.2 | 0.19 | 10.6 | 0.23 |
| Olive oil | 7.7 | 0.16 | 0.4 | 0.04 |
| Processed meat | 7.2 | 0.16 | 0.3 | 0.04 |
| Other oils | 5.5 | 0.14 | 2.4 | 0.11 |
| Dairy products | 4.9 | 0.13 | 0.4 | 0.04 |
| Sodas and juices | 3.8 | 0.11 | 1.8 | 0.10 |
| Cakes and sweets | 3.4 | 0.11 | 0.6 | −0.05 |
| Red meat | 2.6 | 0.09 | 0.0 | 0.01 |
| Vegetables | 2.5 | 0.09 | 5.2 | 0.16 |
| Sweet spreads | 1.7 | 0.08 | 0.0 | −0.01 |
| Cooking fats | 1.5 | 0.07 | 0.2 | −0.03 |
| Eggs | 1.3 | 0.07 | 2.1 | −0.10 |
| Rice and pasta | 1.1 | 0.06 | 1.7 | 0.09 |
| Vegetable soups, stews and sauces | 0.9 | 0.06 | 0.0 | −0.01 |
| Nuts and seeds | 0.5 | 0.04 | 1.5 | 0.09 |
| Fish | 0.1 | 0.02 | 9.5 | 0.22 |
| Meat mixed dishes | 2.1 | −0.08 | 1.3 | −0.08 |
| Fruits | 4.1 | −0.12 | 1.6 | −0.09 |
| Legumes | 4.8 | −0.13 | 2.4 | −0.11 |
| Vegetarian mixed dishes | 7.2 | −0.16 | 0.9 | −0.07 |
| Fermented maize products | 18.4 | −0.25 | 19.8 | −0.32 |
| Refined cereal | 22.4 | −0.27 | 13.3 | −0.26 |
| Roots, tubers and plantains | 28.7 | −0.31 | 22.0 | −0.34 |
| Palm oil | 35.1 | −0.34 | 17.3 | −0.30 |
| Total | 9.9 | 6.5 | ||
1 Factor loadings are correlations between food groups and the dietary pattern scores. Figures in bold represent food items with relevant contributions to the dietary pattern score (≥0.20% explained variation in the factor loadings for either males or females).
Table 3.
Associations of the FLI-related and RRR-derived dietary pattern scores with type 2 diabetes among males and females 1.
Table 3.
Associations of the FLI-related and RRR-derived dietary pattern scores with type 2 diabetes among males and females 1.
| Model | Odds Ratio (95% Confidence Interval) | ||||||
|---|---|---|---|---|---|---|---|
| Q1 | Q2 | Q3 | Q4 | Q5 | p for Trend | Per 1 Score-SD | |
| Men | |||||||
| Diabetes/Control | 17/256 | 22/251 | 31/243 | 46/227 | 43/230 | ||
| Crude | 1 (reference) | 1.32 (0.69, 2.54) | 1.92 (1.06, 3.56) | 3.05 (1.70, 5.47) | 2.82 (1.56, 5.07) | <0.0001 | 1.55 (1.30, 1.86) |
| Model 1 | 1 (reference) | 1.41 (0.71, 2.79) | 1.79 (0.90, 3.56) | 2.30 (1.07, 4.97) | 1.97 (0.89, 4.35) | 0.11 | 1.34 (1.04, 1.73) |
| Model 2 | 1(reference) | 1.25 (0.62, 2.49) | 1.58 (0.79, 3.16) | 2.14 (0.98, 4.68) | 2.03 (0.90, 4.60) | 0.07 | 1.45 (1.10, 1.93) |
| Women | |||||||
| Diabetes/Control | 25/439 | 38/426 | 32/433 | 47/417 | 47/417 | ||
| Crude | 1(reference) | 1.57 (0.93, 2.64) | 1.30 (0.76, 2.23) | 1.98 (1.20, 3.27) | 1.98 (1.20, 3.27) | <0.005 | 1.24 (1.07, 1.44) |
| Model 1 | 1(reference) | 1.59 (0.93, 2.71) | 1.31 (0.75, 2.30) | 2.05 (1.19, 3.54) | 1.98 (1.09, 3.59) | <0.02 | 1.23 (1.03, 1.48) |
| Model 2 | 1(reference) | 1.28 (0.74, 2.21) | 1.03 (0.58, 1.83) | 1.64 (0.94, 2.84) | 1.65 (0.90, 3.02) | 0.07 | 1.16 (0.95, 1.42) |
1 Odds ratios (ORs) at 95% confidence intervals (CIs) were calculated by logistic regression, and the p-values for the trend were calculated by modeling the median of the dietary pattern scores as the independent variable. Model 1: adjusted for age (years) and study site (categorical); Model 2: Model 1 + education (4 categories), energy intake (kcal/day), smoking (yes or no), physical activity (METs-h/week) and alcohol (alcohol/day).
Table 4.
Percentage of explained variation in food intake and factor loadings of 30 food groups of the RRR-derived dietary pattern scores related to NAFLD biomarkers among males and females 1.
Table 4.
Percentage of explained variation in food intake and factor loadings of 30 food groups of the RRR-derived dietary pattern scores related to NAFLD biomarkers among males and females 1.
| Food Group | Men (n = 1366) | Women (n = 2321) | ||
|---|---|---|---|---|
| Explained Variation (%) | Factor Loading | Explained Variation (%) | Factor Loading | |
| Whole-grain cereals | 36.5 | 0.33 | 25.3 | −0.26 |
| Poultry | 26.6 | 0.28 | 35.1 | −0.31 |
| Dairy products | 25.8 | 0.28 | 13.4 | −0.19 |
| Coffee and tea | 23.8 | 0.27 | 44.9 | −0.35 |
| Condiments | 21.4 | 0.25 | 33.0 | −0.30 |
| Potatoes | 19.0 | 0.24 | 28.7 | −0.28 |
| Margarine | 12.8 | 0.20 | 14.6 | −0.20 |
| Olive oil | 13.6 | 0.20 | 18.3 | −0.22 |
| Sodas and juices | 7.9 | 0.15 | 6.1 | −0.13 |
| Sweet spreads | 7.8 | 0.15 | 6.6 | −0.13 |
| Rice and pasta | 7.5 | 0.15 | 0.7 | −0.05 |
| Processed meat | 5.5 | 0.13 | 5.7 | −0.13 |
| Palm oil | 32.7 | −0.31 | 27.1 | 0.27 |
| Roots, tubers and plantains | 30.6 | −0.30 | 12.4 | 0.18 |
| Fermented maize products | 23.5 | −0.26 | 6.1 | 0.13 |
| Vegetarian mixed dishes | 10.5 | −0.18 | 27.3 | 0.27 |
| Refined cereals | 5.6 | −0.13 | 7.4 | 0.14 |
| Cakes and sweets | 4.0 | 0.11 | 7.5 | −0.14 |
| Vegetables | 3.9 | 0.11 | 9.2 | −0.16 |
| Meaty mixed dishes | 4.1 | −0.11 | 0.9 | 0.05 |
| Legumes | 3.9 | −0.11 | 0.5 | −0.04 |
| Other oils | 2.8 | 0.09 | 4.2 | −0.11 |
| Cooking fats | 1.5 | 0.07 | 0.1 | −0.01 |
| Fish | 1.6 | −0.07 | 13.1 | 0.19 |
| Fruits | 1.4 | −0.06 | 1.7 | −0.07 |
| Eggs | 0.9 | 0.05 | 5.9 | −0.13 |
| Vegetable soups, stews and sauces | 0.3 | 0.03 | 0.1 | 0.01 |
| Red meat | 0.3 | −0.03 | 1.1 | −0.05 |
| Nuts and seeds | 0.4 | −0.03 | 0.0 | −0.01 |
| Alcoholic beverages | 0.1 | −0.02 | 6.3 | −0.13 |
| Total | 11.2 | 12.1 | ||
1 Factor loadings are correlations between food groups and the dietary pattern score. Figures in bold represent food groups with relevant contributions to the dietary pattern score (≥0.20% explained variation in the factor loading for either males or females).
Table 5.
Percentage of explained variation in the NALFD biomarkers and response weights of the RRR-derived dietary pattern scores among males and females 1.
Table 5.
Percentage of explained variation in the NALFD biomarkers and response weights of the RRR-derived dietary pattern scores among males and females 1.
| Biomarker | Men (n = 1366) | Women (n = 2321) | ||
|---|---|---|---|---|
| Explained Variation (%) | Response Weight | Explained Variation (%) | Response Weight | |
| Cholesterol | 4.5 | 0.45 | 0.3 | 0.09 |
| LDL-cholesterol | 4.2 | 0.44 | 1.0 | 0.16 |
| HDL-cholesterol | 2.9 | 0.36 | 8.6 | −0.49 |
| AST | 7.4 | −0.58 | 12.3 | 0.58 |
| GGT | 1.5 | −0.26 | 0.3 | 0.09 |
| Triglycerides | 1.3 | −0.24 | 11.1 | 0.55 |
| C-reactive protein | 0.3 | −0.12 | 0.4 | 0.11 |
| ALT | 0.2 | 0.08 | 2.2 | 0.25 |
| Total | 2.80 | 4.50 | ||
1 Figures in bold represent biomarkers with relevant relationships with the dietary pattern scores (response weight >|0.35|).
Table 6.
Associations of the biomarker-related and RRR-derived dietary pattern scores with type 2 diabetes among males and females 1.
Table 6.
Associations of the biomarker-related and RRR-derived dietary pattern scores with type 2 diabetes among males and females 1.
| Odds Ratio (95% Confidence Interval) | |||||||
|---|---|---|---|---|---|---|---|
| Q1 | Q2 | Q3 | Q4 | Q5 | p for Trend | Per 1 Score-SD | |
| Men | |||||||
| Diabetes/Control | 15/258 | 27/246 | 33/241 | 47/226 | 37/236 | ||
| Crude | 1 (reference) | 1.89 (0.98, 3.63) | 2.36 (1.25, 4.44) | 3.58 (1.95, 6.57) | 2.7 0(1.44, 5.04) | <0.0002 | 1.42 (1.20, 1.68) |
| Model 1 | 1 (reference) | 1.91 (0.97, 3.76) | 2.06 (1.00, 4.22) | 2.42 (1.07, 5.50) | 1.64 (0.72, 3.75) | 0.578 | 1.15 (0.89, 1.48) |
| Model 2 | 1 (reference) | 1.76 (0.89, 3.50) | 1.8 (0.86, 3.75) | 2.07 (0.90, 4.78) | 1.44 (0.61, 3.40) | 0.743 | 1.13 (0.86, 1.49) |
| Model 3 | 1 (reference) | 1.62 (0.81–3.22) | 1.46 (0.69–3.10) | 1.70 (0.72–4.01) | 1.21 (0.50, 2.89) | 0.962 | 1.07 (0.81, 1.42) |
| Women | |||||||
| Diabetes/Control | 37/427 | 38/426 | 36/429 | 36/430 | 44/420 | ||
| Crude | 1 (reference) | 1.03 (0.64, 1.65) | 0.97 (0.60, 1.56) | 0.91 (0.56, 1.48) | 1.21 (0.77, 1.91) | 0.592 | 1.08 (0.93, 1.25) |
| Model 1 | 1 (reference) | 1.24 (0.76, 2.05) | 1.63 (0.84, 3.14) | 1.46 (0.73, 2.92) | 1.69 (0.84, 3.39) | 0.229 | 1.35 (1.05, 1.73) |
| Model 2 | 1 (reference) | 1.16 (0.70, 1.94) | 1.39 (0.69, 2.79) | 1.20 (0.57, 2.52) | 1.40 (0.67, 2.93) | 0.522 | 1.29 (0.99, 1.68) |
| Model 3 | 1 (reference) | 1.11 (0.66, 1.86) | 1.32 (0.65, 2.67) | 1.20 (0.57, 2.57) | 1.42 (0.67, 3.02) | 0.422 | 1.30 (0.99, 1.71) |
1 Odds ratios (ORs) at 95% confidence intervals (CIs) were calculated by logistic regression, and the p-values for the trend were calculated by modeling the median of the dietary pattern scores as the independent variable, Q1 as quintile 1, etc. Model 1: adjusted for age (years) and study site (categorical); Model 2: Model 1 + education (4 categories), energy intake (kcal/day), smoking (yes or no), physical activity (METs-h/week); Model 3: Model 2+ BMI and waist circumference.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. |
© 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
MDPI and ACS Style
Osei, T.B.; van Dijk, A.-M.; Dingerink, S.; Chilunga, F.P.; Beune, E.; Meeks, K.A.C.; Bahendeka, S.; Schulze, M.B.; Agyemang, C.; Nicolaou, M.; et al. Reduced Rank Regression-Derived Dietary Patterns Related to the Fatty Liver Index and Associations with Type 2 Diabetes Mellitus among Ghanaian Populations under Transition: The RODAM Study. Nutrients 2021, 13, 3679. https://doi.org/10.3390/nu13113679
AMA Style
Osei TB, van Dijk A-M, Dingerink S, Chilunga FP, Beune E, Meeks KAC, Bahendeka S, Schulze MB, Agyemang C, Nicolaou M, et al. Reduced Rank Regression-Derived Dietary Patterns Related to the Fatty Liver Index and Associations with Type 2 Diabetes Mellitus among Ghanaian Populations under Transition: The RODAM Study. Nutrients. 2021; 13(11):3679. https://doi.org/10.3390/nu13113679
Chicago/Turabian StyleOsei, Tracy Bonsu, Anne-Marieke van Dijk, Sjoerd Dingerink, Felix Patience Chilunga, Erik Beune, Karlijn Anna Catharina Meeks, Silver Bahendeka, Matthias Bernd Schulze, Charles Agyemang, Mary Nicolaou, and et al. 2021. "Reduced Rank Regression-Derived Dietary Patterns Related to the Fatty Liver Index and Associations with Type 2 Diabetes Mellitus among Ghanaian Populations under Transition: The RODAM Study" Nutrients 13, no. 11: 3679. https://doi.org/10.3390/nu13113679
APA StyleOsei, T. B., van Dijk, A. -M., Dingerink, S., Chilunga, F. P., Beune, E., Meeks, K. A. C., Bahendeka, S., Schulze, M. B., Agyemang, C., Nicolaou, M., Holleboom, A. G., & Danquah, I. (2021). Reduced Rank Regression-Derived Dietary Patterns Related to the Fatty Liver Index and Associations with Type 2 Diabetes Mellitus among Ghanaian Populations under Transition: The RODAM Study. Nutrients, 13(11), 3679. https://doi.org/10.3390/nu13113679
Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.
