Erythrocyte n-6 polyunsaturated fatty acids, gut microbiota and incident type 2 diabetes: a prospective cohort study ==================================================================================================================== * Zelei Miao * Jie-sheng Lin * Yingying Mao * Geng-dong Chen * Fang-fang Zeng * Hong-li Dong * Zengliang Jiang * Jiali Wang * Congmei Xiao * Menglei Shuai * Wanglong Gou * Yuanqing Fu * Fumiaki Imamura * Yu-ming Chen * Ju-Sheng Zheng ## Abstract **Objective** To examine the association of erythrocyte n-6 polyunsaturated fatty acid (PUFA) biomarkers with incident type 2 diabetes (T2D) and explore the potential role of gut microbiota in the association. **Design** We evaluated 2,731 non-T2D participants recruited between 2008-2013 in the Guangzhou Nutrition and Health Study, China. T2D cases were identified with clinical and biochemical information collected at follow-up visits. Using stool samples collected during the follow-up in the subset (n=1,591), 16S rRNA profiling was conducted. Using multivariable-adjusted Poisson or linear regression, we examined associations of erythrocyte n-6 PUFA biomarkers with incident T2D, and diversity and composition of gut microbiota. **Results** Over 6.2 years of follow-up, 276 T2D cases were identified (risk=0.10). Higher levels of erythrocyte γ-linolenic acid (GLA), but not linoleic or arachidonic acid, were associated with higher T2D incidence. Comparing the top to the bottom quartile groups of GLA levels, relative risk was 1.72 (95% confidence intervals: 1.21, 2.44) adjusted for potential confounders. Baseline GLA was inversely associated with gut microbial richness and diversity (α-diversity, both *p*<0.05) during follow-up, and significantly associated with microbiota β-diversity (*p*=0.002). Seven genera (*Butyrivibrio, Blautia, Oscillospira, Odoribacter, S24-7 other, Rikenellaceae other*, and *Clostridiales other*) were enriched in quartile 1 of GLA, and in participants without T2D. **Conclusion** Relative concentrations of erythrocyte GLA were positively associated with incident T2D in a Chinese population and also with gut microbial profiles. These results highlight that gut microbiota may play an important role linking n-6 PUFA metabolism and T2D etiology. ## INTRODUCTION Type 2 diabetes (T2D) is one of the most prevalent metabolic conditions worldwide, and 463 million adults are living with diabetes. 1 To examine lifestyle-related risk factors for the development of T2D, the role of different fatty acids has been an important research agenda. During the past two decades, consumption of n-6 polyunsaturated fatty acids (PUFAs), rich in vegetable oil, has been increasing rapidly,2,3 while their relationship with human health remains controversial. The metabolic pathway of n-6 PUFAs is hypothesized to be closely involved in the T2D etiology,4 but has not been well characterized yet. A recent pooled analysis of 20 prospective cohort studies investigated the association of linoleic acid ((LA; 18:2n6) and arachidonic acid (AA; 20:4n6) biomarkers with incident T2D, reporting that LA, but not AA, was inversely associated with T2D.5 However, in this pooling project, only two individual n-6 PUFA biomarkers were evaluated while cohorts in Europe reported gamma-linolenic acid (GLA) and other n-6 PUFAs to be associated with higher T2D risk.5 Moreover, with only one small cohort from Taiwan in the above pooling project,5 evidence from Asian populations has been limited. As an Asian population has different metabolic and lifestyle characteristics in relation to diabetes compared to American or European populations,6 more investigation for n-6 PUFA biomarkers in Asian populations is highly warranted. The protective effect of n-6 PUFA or LA, the most abundant n-6 PUFA, on insulin homeostasis has been well characterized and discussed to link n-6 PUFA exposure to lower T2D risk.7,8 However, gut microbiota has been little investigated in research on PUFAs and any non-communicable diseases. On one hand, gut microbiota has been considered as an important risk factor for T2D and insulin resistance.9,10 On the other, rodent studies have suggested that high tissue levels of n-6 PUFAs are associated with differences in gut microbiota, including increased *Proteobacteria*, reduced *Actinobacteria* and *Bacteroidetes*, and altered microbial α-diversity.11–13 Yet, to our knowledge, there are no cohort studies investigating the potential link between n-6 PUFA exposure, T2D risk, and gut microbiota composition and diversity. Evaluating the prospective cohort in China, we aimed to examine the prospective associations 1) between n-6 PUFA biomarkers and T2D risk in a Chinese population, and 2) between n-6 PUFA biomarkers and gut microbiota composition and diversity, and subsequent health outcomes. ## METHODS ### Study design and population This study was based on the Guangzhou Nutrition and Health Study (GNHS), a community-based prospective cohort in southern China. Detailed information of the participants and study design was described previously.14 Briefly, 4,048 Chinese participants aged 40-75 years and lived in Guangzhou for at least 5 years were recruited between 2008-2010 (n=3,169) and 2012-2013 (n=879). All participants have been followed up every 3 years approximately. During each visit, demographic, dietary and lifestyle information was collected via questionnaires. The study was registered at [clinicaltrials.gov](http://clinicaltrials.gov) ([NCT03179657](http://medrxiv.org/lookup/external-ref?link_type=CLINTRIALGOV&access_num=NCT03179657&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F30%2F2020.03.29.20039693.atom)). According to our pre-specified inclusion and exclusion criteria of our participants for the present analysis (Figure S1), we excluded participants with missing information on age or sex (n=9), diet (n=7), self-reported/diagnosed T2D (n=400) or history of cancer (n=16) at baseline, or extreme levels of self-reported total energy intake (<800 kcal or >4,000 kcal per day for men, and <500 kcal or > 3,500 kcal per day for women) (n=53). We also excluded those without measurement of erythrocyte membrane fatty acid compositions at baseline (n=298) and further excluded those without any follow-up information on T2D status (n=534, follow-up rate 84%). The current analyses were censored by date of T2D ascertainment or Apr 30, 2019, which ever happened first. A total of 2,731 individuals were included in the present study, with a median follow-up time of 6.2 years. Among the above included 2,731 individuals, stool samples from 1,606 individuals were collected between 2014-2018 for the subsequent 16S profiling, with 1,591 individuals included in the microbial association analysis after excluding participants who used antibiotics within two weeks before stool collection (n=15). ### Measurement of erythrocyte membrane fatty acids and covariates We considered relative concentrations (% of total fatty acids) of erythrocyte LA, GLA, and AA as the main exposure variables. Blood samples were drawn after overnight fasting at the baseline visit. Erythrocytes were aliquoted within 2 hours of blood sampling and stored at −80°C. As previously described,15 fatty acid moieties of erythrocyte membranes were trans-methylated and measured as proportions (%) of total fatty acids by using gas chromatography (7890 GC, DB-23 capillary column 60m×0.25mm internal diameter×0.15μm film, Agilent, California, USA). Commercially available standards (Nu-Chek Prep, Minnesota, USA) were used to identify individual fatty acids (n=37) and quantify a relative peak strength of each. The intra-assay coefficients of variation for LA, GLA and AA were 6.4%, 12.8%, and 8.04%, respectively. Sociodemographic, lifestyle, and clinical covariates were collected with standardized questionnaires. Details of fatty acid and covariate measurements are available in Supplementary Methods. ### T2D case ascertainment At baseline (for exclusion) and follow-up visits, T2D cases were ascertained according to the criteria of American Diabetes Association for T2D diagnosis16 if a participant met one of the following criteria: a fasting blood glucose concentration ≥7.0 mmol/L (126 mg/dL), or HbA1c concentration ≥6.5%, or self-reported medical treatment for diabetes. ### Fecal sample collection and 16S rRNA profiling Fecal samples were collected on site during follow-up visits between 2014-2018. Microbial DNA extraction, polymerase chain reaction (PCR) and amplicon sequencing were performed as described (see Supplementary Methods). Fastq-files were demultiplexed, merge-paired, quality filtered by Quantitative Insights into Microbial Ecology software (version 1.9.0).17 Sequences were clustered into operational taxonomic units (OTUs) with 97% similarity and annotated based on the Greengenes Database (version 13.8).18 ### Statistical analysis Statistical analyses were performed using Stata 15 (StataCorp, Texas, USA) or R (version 3.6.3). We used multivariable Poisson regression models to estimate the risk ratio (RR) and 95% confidence interval (CI) of T2D comparing quartiles of the erythrocyte n-6 PUFA (LA, GLA and AA), adjusting for potential confounders. We fitted three statistical models: model 1 included age, sex, BMI, and waist-hip ratio (WHR); model 2, model 1 + education, household income, smoking, alcohol drinking status, physical activity, total energy intake, and family history of diabetes; and model 3, as model 2 + baseline erythrocyte total n-3 PUFAs and fasting glucose. *p*-value for trend was calculated based on per quartile difference in the corresponding n-6 PUFA variable. In sensitivity analyses, we repeated the above analysis based on model 3 to examine the impact of loss to follow up by a simple imputation (assuming that participants loss to follow up did not develop T2D), and 10 rounds of multiple imputations based on regression model using all covariates listed in model 3. We assessed influences of additional potential covariates on the models: dietary variables (dietary intake of dairy products, red and processed meat, fish, vegetable, and fruit, in quartiles), blood lipids (total triglycerides and low-density lipoprotein), prevalent coronary heart disease, treatment for hypertension or hyperlipidemia. We further examined the potential interaction of erythrocyte n-6 PUFA biomarkers with age, sex, BMI or n-3 PUFA biomarker concentrations on T2D risk. Post-hoc stratified analyses were performed if there was a significant interaction (*p**interaction*<0.05) for age, sex, BMI or n-3 PUFA biomarker concentrations with the n-6 PUFA biomarkers. For the gut microbiome analysis, we calculated the α-diversity metrics (observed OTUs, Shannon’s diversity index, Chao index and Simpson index) by using QIIME (version 1.9.0) based on rarefied OTU counts. We standardized α-diversity indicators (Z-score was calculated: the variable was subtracted by the mean and divided by the standard deviation) and then reduced assay-dependent variabilities by fitting a linear mixed model with an α-diversity indicator as an outcome variable and with assay-specific covariates including sequencing depth and Bristol scale as fixed effects, and sequencing batch as a random effect. The residuals from this modeling were then used in the following n-6 PUFA biomarker association analyses. We used a linear regression model to examine the association of individual n-6 PUFA biomarkers (by quartile) with the four α-diversity metrics, adjusted for age, sex, BMI, waist-hip ratio, education, household income, smoking status, alcohol drinking status, physical activity, total energy intake, and baseline erythrocyte total n-3 PUFAs. We also conducted an additional adjustment for dietary fiber intake to assess the potential marginal effect of dietary fiber, as fiber is a well-known factor influencing the gut microbiome.19 To investigate the relevance of significant findings in relation to T2D, we performed mediation analysis (R {*mediation*})20 with the residuals of α-diversity as a mediator of the relationship between n-6 PUFAs and T2D, adjusting for all covariates described in model 3. We examined the cross-sectional associations between the residuals of α-diversity metrics and T2D by using logistic regression models, adjusted for age, sex, BMI, education, household income, smoking status, alcohol drinking status, prevalent hypertension and dyslipidemia. We conducted a principal coordinate analysis (PCoA) and permutational multivariate analysis of variance (PERMANOVA) (R function adonis {*vegan*}, 999 permutations)21 based on Bray-Curtis distance to compare the whole gut composition (β-diversity) at genus level by quartiles of n-6 PUFA biomarkers. The potential confounders included in the PERMANOVA were sequencing depth, sequencing batch, Bristol scale, age, sex, BMI, WHR, education, household income, smoking status, alcohol drinking status, physical activity, total energy intake, and baseline erythrocyte total n-3 PUFAs. The above analysis was repeated for gut composition at OTU level. Scaled relative abundances (i.e., divided by the standard deviation) were used to lessen the influence of highly abundant genera on Bray-Curtis distance. The association of community structure with T2D was also assessed by PERMANOVA (R function adonis, 999 permutations) to obtain R2 indicating the proportion of the variability explained by the studied variables including n-6 PUFAs. We performed the post-hoc pairwise comparisons between different quartiles of GLA with Bonferroni correction (cut-off *p*<0.008). At genus level, we performed a linear discrimination analysis (LDA) as implemented using LEfSe.22 The default parameters were used (α value for Wilcoxon tests was 0.05 and the LDA score was 2.0) to identify biomarkers at genus level distinguishing different quartiles of n-6 PUFAs. Only taxon present in at least 10% of samples were included in the analyses. Then, we compared the relative abundances of genera by T2D status using the LEfSe method.22 Besides, we evaluated the potential correlation of the above identified genera with T2D-related traits, including fasting insulin, glucose, total cholesterol, total triglycerides, low-density lipoprotein, high-density lipoprotein cholesterol (HDL), non-HDL, HbA1c, homeostasis model assessment of insulin resistance (HOMA-IR) and β-cell function (HOMA-β). We first modelled T2D-related traits as outcome variables in linear regression with age, sex BMI and obtained the residuals of each T2D-related trait respectively, then we calculated the Spearman correlation coefficients between the genera and the residuals of T2D-related traits. Results displaying a *p*<0.05 after Bonferroni adjustment were considered as significant. ## RESULTS ### Characteristics of the study participants Table 1 and Table S1 summarize the baseline characteristics of the population by quartiles of erythrocyte GLA, LA and AA. The median levels (interquartile range) of baseline erythrocyte LA, GLA, AA and total n-6 PUFAs (% total fatty acids) were 9.84 (8.89-10.76), 0.037 (0.025-0.052), 11.51 (10.00-12.60) and 21.6 (19.6-23.1), respectively. Table S2 and S3 present baseline population characteristics of participants with and without follow-up information, 16S profiling respectively. Participants who dropped out had higher serum glucose levels, consumed fewer vegetables and energy, tended to be older, less educated and physically active, and were more likely to smoke. View this table: [Table 1.](http://medrxiv.org/content/early/2020/03/30/2020.03.29.20039693/T1) Table 1. Baseline population characteristics by quartiles of erythrocyte γ-linolenic acid* ### Association of erythrocyte n-6 PUFAs with incident T2D Higher levels of baseline GLA were associated with higher risk of T2D across the three statistical models (*p*-trend < 0.001 for all models, Table 2). For model 3, compared with the first quartile group (Q1), RRs (95% CIs) of T2D at Q2, Q3 and Q4 were 1.22 (0.85, 1.74), 1.43 (1.01, 2.03) and 1.72 (1.21, 2.44), respectively. No association was found for LA or AA. Similar results were obtained after using the simple or multiple imputation for the missing data in the sensitivity analysis (Table S4). Adjustment for additional dietary factors and disease histories showed similar results (Table S5). We did not observe significant interaction for n-6 PUFA biomarkers with age, sex, BMI or n-3 PUFAs on T2D risk (*p*interaction>0.05). View this table: [Table 2.](http://medrxiv.org/content/early/2020/03/30/2020.03.29.20039693/T2) Table 2. Association of erythrocyte n-6 fatty acids with incident type 2 diabetes* ### The prospective association of erythrocyte n-6 PUFA biomarkers with gut microbiota diversity Lower GLA was associated with higher gut microbiota α-diversity (Q4 *vs*. Q1: *p*=0.021 for observed OTUs, 0.028 for Shannon’s diversity index, 0.023 for Chao index, respectively) (Figure 1). However, baseline LA, AA or total n-6 PUFAs were not associated with the microbial α-diversity (Table S6). These results remained similar after further adjusting for dietary fiber intake (Table S7). Among 1,591 individuals included in the microbial association analysis, we also observed positive association between GLA and T2D (Table S8). α-diversity acted as a potential partial mediator of the association between GLA and T2D (Table 3). For β-diversity, GLA and AA significantly (*p*<0.01) contributed to dissimilarities in gut microbiota composition at genus level (Table S9). Specifically, for GLA, the gut microbiota composition at genus level was different across quartile groups of GLA concentrations (Table S10), with *p* = 0.006 for comparison between Q1 and Q4, for example, as visualized in Figure 2. View this table: [Table 3.](http://medrxiv.org/content/early/2020/03/30/2020.03.29.20039693/T3) Table 3. Summary statistics of the mediation analysis for α-diversity indicators* ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/30/2020.03.29.20039693/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2020/03/30/2020.03.29.20039693/F1) Figure 1. Gut microbial community richness and diversity between quartile 1 and quartile 4 of erythrocyte γ-linolenic acid. The box plots feature the median (centerline), upper and lower quartiles (box limits), 1.5× the interquartile range (whiskers). Gut microbial community richness was represented by Observed OTUs (A) and Chao index (B), and community diversity was represented by Shannon’s diversity index (C) and Simpson index (D). *p* values were calculated for Q4 (n=397) of the erythrocyte n-6 fatty acids using Q1 (n=398) as the reference group using a linear mixed model: sequencing batch as random effect; sequencing depth, Bristol scale, age, sex, BMI, waist-hip ratio, education, household income, smoking and alcohol drinking status, physical activity, total energy intake, family history of diabetes, baseline erythrocyte total n-3 PUFAs and fasting glucose as fixed effects. Abbreviations: BMI, body mass index, Q, quartile. ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/30/2020.03.29.20039693/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2020/03/30/2020.03.29.20039693/F2) Figure 2. Dissimilarities in gut microbiota composition (β-diversity) between quartile 1 and quartile 4 of erythrocyte γ-linolenic acid. Dissimilarities in gut microbiota composition between quartile 1 and quartile 4 of erythrocyte γ-linolenic acid (n=795) represented by unconstrained PCoA with the Bray-Curtis dissimilarity index calculated on scaled relative abundances at genus level. Baseline erythrocyte γ-linolenic acid level explained 0.3% of the dissimilarities in gut microbiota composition (PERMANOVA, *p*= 0.002). Abbreviation: Q, quartile. ### The cross-sectional association of gut microbiota diversity with T2D T2D patients had a lower α-diversity represented by observed OTUs, Shannon’s diversity index, Simpson index and Chao index compared with patients without T2D (Table S11). We also found that the overall β-diversity was different between participants with and without T2D (Table S12). ### Taxonomic profiles of gut microbial community in different levels of GLA We then focused on exploring gut genera potentially underlying the association between GLA and incident T2D. At the genus level, the gut microbiota of participants within Q1 and Q4 of GLA was both dominated by *Bacteroides*, while the relative abundance of *Bacteroides* was significantly higher in Q4 group (*p*=0.016). We also identified *Rothia, [Eubacterium], Coprococcus* enriched in participants with high GLA levels (Q4) using LDA method (Figure 3). We found several overlapping taxonomic biomarkers for non-T2D status and low GLA status, including *Butyrivibrio, Blautia, Oscillospira, Odoribacter, S24-7 other, Rikenellaceae other*, and *Clostridiales other* (Figure 3). In the cross-sectional analysis relating taxonomic factors to metabolic traits, the genus *[Eubacterium], S24-7 other, Blautia, Oscillospira, Odoribacter, Rikenellaceae other, Coriobacteriaceae Other, Faecalibacterium* and *Christensenellaceae Other* were negatively correlated with total triglycerides (Bonferroni corrected *p*<0.05) (Figure 4). No correlation was observed for other metabolic trait-microbiota pairs (Table S13). ![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/30/2020.03.29.20039693/F3.medium.gif) [Figure 3.](http://medrxiv.org/content/early/2020/03/30/2020.03.29.20039693/F3) Figure 3. Gut microbial biomarkers at genus level by γ-linolenic acid and type 2 diabetes. Gut microbial taxonomic biomarkers identified by the linear discriminant analysis (LDA) at quartile 1 and 4 of erythrocyte γ-linolenic acid (A), and by type 2 diabetes status (B). Color indicates the group in which a differentially abundant taxon is enriched ((A) dark blue: quartile 1 of GLA; red: quartile 4 of GLA; (B) gray: participants without T2D; red: participants with T2D). Abbreviation: Q, quartile; T2D, type 2 diabetes; Non-T2D, participants without type 2 diabetes. ![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/30/2020.03.29.20039693/F4.medium.gif) [Figure 4.](http://medrxiv.org/content/early/2020/03/30/2020.03.29.20039693/F4) Figure 4. Heatmap of the Spearman correlation coefficients between γ-linolenic acid-related microbes and T2D-related traits. The intensity of the colors represents the degree of association between γ-linolenic acid-related microbes and ten T2D-related traits as measured by the Spearman’s correlations. All significant correlations are marked with asterisk (Bonferroni corrected *p*<0.05). Abbreviation: TG, total triglycerides; TC, total cholesterol; LDL, low-density lipoprotein cholesterol; HDL, high-density lipoprotein cholesterol; HOMA-IR, homeostasis model assessment of insulin resistance; HOMA-β, homeostasis model assessment of β-cell function; T2D, type 2 diabetes. ## DISCUSSION In the present community-based prospective cohort of a Chinese population, we found that baseline erythrocyte GLA levels were positively associated with incident T2D, independent of BMI and other potential confounders. Proportions of LA, AA or total n-6 PUFAs were not associated with T2D incidence. Our further investigation integrating gut microbiota data revealed that GLA may be associated with T2D risk through a mechanism which varies diversity and composition of gut microbiota. Our findings suggest that n-6 PUFAs and gut microbiota co-vary in the development of T2D risk, highlighting the presence of a novel mechanism of how fatty acids or gut microbiota influence T2D risk. Our findings suggest that high circulating GLA might be a risk factor of T2D, which are consistent with several prospective studies in western countries.23,24 The null finding for AA in our present study is also consistent with those from the previous studies.5,23,25 The findings for LA were inconsistent between ours and previous ones.5,23 Circulating LA reflects habitual consumption of n-6 PUFA, typically from plant-derived oils. The inconsistency may reflect that European and Chinese populations use plant oils differently with a combination of different foods. For instance, Chinese adults may consume plant oils with meat by which positive and negative effects of those would lead to the null finding for LA. By contrast, an American or European population may use plant oils with vegetables or in a healthy diet pattern. Such dietary confounding might influence GLA finding little, as GLA is a product of tightly-regulated desaturation of LA, and result in the consistent findings between studies. Given the putative relationship between gut microbiome and T2D, we hypothesized and confirmed that gut microbiome diversity may play a role linking n-6 PUFAs and T2D and that n-6 PUFAs may prospectively affect gut microbiota composition. Our findings are in line with a randomized controlled trial in which α-diversity, *Blautia* was reduced, while *Bacteroides* were increased after higher-fat diet intervention (rich in n-6 PUFAs owing to exclusive use of soybean oil).26 Although we did not observe higher relative abundance of *Bacteroides* in T2D patients in our study, *Bacteroides* has been found to be more abundant in T2D in another Chinese population.9 We found that *Butyrivibrio, Blautia, Oscillospira*, and *Odoribacter*, genera known to contain a series of butyrate-producing bacteria, are all enriched in participants within Q1 of GLA and participants without T2D. *Butyrivibrio fibrisolvens* was reported to have the ability to utilize fructose polymers and it showed anticarcinogenic effect in mice.27,28 *Blautia* was reported to have a lower relative abundance in T2D patients and have inverse association with visceral fat.29,30 *Oscillospira* was positively associated with leanness and was reduced in inflammatory diseases.31 In addition, evidence suggested that some species of *Oscillospira* were able to utilize glucuronate, a common animal-derived sugar.32 *Odoribacter splanchnicus* was considered as a potential probiotics candidate for obesity prevention as a higher abundance in obesity-depleted group.33 Moreover, in our study, *Blautia, Oscillospira* and *Odoribacter* were inversely associated with circulating triglycerides, which was found to be positively associated with T2D and higher level of GLA.34,35 Thus, our study suggests that lower level of GLA may be correlated with factors increasing the abundance of these bacteria and preventing the development of T2D. The more detailed mechanism behind these associations remains unclear but may be related to chronic inflammation. Accumulating evidence supports the hypothesis that chronic low-grade inflammation is a risk factor for the development of T2D,36 and that n-6 PUFA may play a vital role in the regulation of inflammation.37 GLA is produced from LA by the enzyme delta-6-desaturase, and is further metabolized to dihomo-γ-linolenic acid (DGLA; 20:3n6) which undergoes oxidative metabolism by cyclooxygenases and lipoxygenases to produce anti-inflammatory eicosanoids such as prostaglandins and leukotrienes.38 Moreover, imbalance of the gut microbiota and some specific microbes have been considered as contributors for inflammation by gut microbiota-derived lipopolysaccharide(LPS)39,40. The LPS stimulates Toll-like receptors (TLRs) to influence innate immunity via its interaction with LPS receptors (i.e., CD14/TLR 4complex).41 It can be proposed that LPS level is regulated by fat and fatty acids, and associated with insulin resistance and T2D.40,42,43 Moreover, the enrichment of inflammation-related bacteria (i.e., *Blautia*44,45) in participants with lower levels of GLA supports an important role of inflammation in the relationship between GLA, gut microbiota and T2D. Strengths of the present study include the prospective cohort study design, which is rarely conducted in prior research on fatty acid-gut microbiota associations, and the collection and adjustment for a variety of information including socio-demographic, lifestyle and dietary factors, anthropometric parameters and circulating blood biomarkers. Several limitations merit attention. First, although we have adjusted for several potential confounders, due to the observational nature of the present study, we are not able to fully exclude the influence of residual confounding. Second, erythrocyte n-6 PUFAs were measured at baseline only and may not represent long-term status. Third, fecal samples were only collected during follow-up but not at baseline. Repeated measure of both circulating n-6 PUFAs and gut microbiota or an intervention of n-6 PUFAs would enable us to characterize bidirectional effects of n-6 PUFA intakes on gut microbiota and of gut microbiota on circulating n-6 PUFAs. Nevertheless, our study fills the gap of the prospective associations between n-6 PUFA biomarkers and gut microbiota. In conclusion, our study suggests that erythrocyte GLA biomarker is positively associated with incident T2D in the Chinese population. Gut microbiota may play an important role to explain the finding. The identification of GLA-related gut microbiota features may serve as potential intervention targets for further investigation. The detailed mechanism behind the association of GLA with T2D risk and gut microbiota needs further clarification in the future. ## Data Availability The raw data for 16â∈%S rRNA gene sequences are available in the CNSA ([https://db.cngb.org/cnsa/](https://db.cngb.org/cnsa/)) of CNGBdb at accession number CNP0000829. ## Contributors JSZ, YMC conceived the study concept and design. GDC, FFZ, HLD, JW, CX collected the data. ZM, JSL, YM did the data analysis. JSZ, ZM, JSL, YM wrote the first draft of the article. ZJ, MS, WG, YF, FI, YMC contributed to the discussion and critical revision of the manuscript. All authors read, revised and approved the final draft. ## Funding This study was funded by the National Natural Science Foundation of China (81903316, 81773416), Zhejiang Ten-thousand Talents Program (101396522001), Westlake University (101396021801) and the 5010 Program for Clinical Researches (2007032) of the Sun Yat-sen University (Guangzhou, China). The funder had no role in study design, data collection and analysis, decision to publish, or writing of the manuscript. ## Competing interests None declared ## Patient consent Obtained ## Ethics approval This study was approved by the Ethics Committee of the School of Public Health at Sun Yat-sen University (2018048) and Ethics Committee of Westlake University (20190114ZJS0003). ## Acknowledgments We thank all the participants and staff of the GNHS study. We thank Westlake University Supercomputer Center for computational resource and related assistance. * Received March 29, 2020. * Revision received March 29, 2020. * Accepted March 30, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## REFERENCES 1. 1.IDF. Diabetes Atlas. 9th Edition.; 2019. [https://diabetesatlas.org/en/](https://diabetesatlas.org/en/). Accessed March 18, 2020. 2. 2.Simopoulos AP. An increase in the Omega-6/Omega-3 fatty acid ratio increases the risk for obesity. Nutrients. 2016;8(3). doi:10.3390/nu8030128 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/nu8030128&link_type=DOI) 3. 3.Shen X, Fang A, He J, et al. Trends in dietary fat and fatty acid intakes and related food sources among Chinese adults: A longitudinal study from the China Health and Nutrition Survey (1997-2011). Public Health Nutr. 2017;20(16):2927–2936. doi:10.1017/S1368980017001781 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S1368980017001781&link_type=DOI) 4. 4.Patterson E, Wall R, Fitzgerald GF, Ross RP, Stanton C. Health implications of high dietary omega-6 polyunsaturated fatty acids. J Nutr Metab. 2012;2012. doi:10.1155/2012/539426 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1155/2012/539426&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22570770&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F30%2F2020.03.29.20039693.atom) 5. 5.Wu JHY, Marklund M, Imamura F, et al. Omega-6 fatty acid biomarkers and incident type 2 diabetes: pooled analysis of individual-level data for 39[740 adults from 20 prospective cohort studies. Lancet Diabetes Endocrinol. 2017;5(12):965–974. doi:10.1016/S2213-8587(17)30307-8 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S2213-8587(17)30307-8&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29032079&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F30%2F2020.03.29.20039693.atom) 6. 6.Chan JCN, Malik V, Jia W, et al. Diabetes in Asia: epidemiology, risk factors, and pathophysiology. JAMA - J Am Med Assoc. 2009;301(20):2129–2140. doi:10.1001/jama.2009.726 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.2009.726&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19470990&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F30%2F2020.03.29.20039693.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000266328800025&link_type=ISI) 7. 7.Varga T, Czimmerer Z, Nagy L. PPARs are a unique set of fatty acid regulated transcription factors controlling both lipid metabolism and inflammation. Biochim Biophys Acta - Mol Basis Dis. 2011;1812(8):1007–1022. doi:10.1016/j.bbadis.2011.02.014 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.bbadis.2011.02.014&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21382489&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F30%2F2020.03.29.20039693.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000292350400023&link_type=ISI) 8. 8.Imamura F, Sharp SJ, Koulman A, et al. A combination of plasma phospholipid fatty acids and its association with incidence of type 2 diabetes: The EPIC-InterAct case-cohort study. PLoS Med. 2017;14(10). doi:10.1371/journal.pmed.1002409 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pmed.1002409&link_type=DOI) 9. 9.Wang J, Qin J, Li Y, et al. A metagenome-wide association study of gut microbiota in type 2 diabetes. Nature. 2012;490(7418):55–60. doi:10.1038/nature11450 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature11450&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23023125&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F30%2F2020.03.29.20039693.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000309446800031&link_type=ISI) 10. 10.Caricilli AM, Saad MJA. The Role of Gut Microbiota on Insulin Resistance. Nutrients. 2013;5:829–851. doi:10.3390/nu5030829 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/nu5030829&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23482058&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F30%2F2020.03.29.20039693.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000316610400014&link_type=ISI) 11. 11.Lee HC, Yu SC, Lo YC, Lin IH, Tung TH, Huang SY. A high linoleic acid diet exacerbates metabolic responses and gut microbiota dysbiosis in obese rats with diabetes mellitus. Food Funct. 2019;10(2):786–798. doi:10.1039/c8fo02423e [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1039/c8fo02423e&link_type=DOI) 12. 12.Kaliannan K, Li XY, Wang B, et al. Multi-omic analysis in transgenic mice implicates omega-6/omega-3 fatty acid imbalance as a risk factor for chronic disease. Commun Biol. 2019;2(1). doi:10.1038/s42003-019-0521-4 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s42003-019-0521-4&link_type=DOI) 13. 13.Turnbaugh PJ, Bäckhed F, Fulton L, Gordon JI. Diet-Induced Obesity Is Linked to Marked but Reversible Alterations in the Mouse Distal Gut Microbiome. Cell Host Microbe. 2008;3(4):213–223. doi:10.1016/j.chom.2008.02.015 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.chom.2008.02.015&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18407065&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F30%2F2020.03.29.20039693.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000255326100005&link_type=ISI) 14. 14.Zhang ZQ, He LP, Liu YH, Liu J, Su YX, Chen YM. Association between dietary intake of flavonoid and bone mineral density in middle aged and elderly Chinese women and men. Osteoporos Int. 2014;25(10):2417–2425. doi:10.1007/s00198-014-2763-9 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00198-014-2763-9&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25062726&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F30%2F2020.03.29.20039693.atom) 15. 15.Zeng F, Sun L, Liu Y, et al. Higher Erythrocyte n–3 PUFAs Are Associated with Decreased Blood Pressure in Middle-Aged and Elderly Chinese Adults. J Nutr. 2014;144(8):1240–1246. doi:10.3945/jn.114.192286 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6OToibnV0cml0aW9uIjtzOjU6InJlc2lkIjtzOjEwOiIxNDQvOC8xMjQwIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDMvMzAvMjAyMC4wMy4yOS4yMDAzOTY5My5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 16. 16.Diagnosis and classification of diabetes mellitus. Diabetes Care. 2013;36(SUPPL.1). doi:10.2337/dc13-S067 [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiZGlhY2FyZSI7czo1OiJyZXNpZCI7czoxOToiMzYvU3VwcGxlbWVudF8xL1M2NyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzAzLzMwLzIwMjAuMDMuMjkuMjAwMzk2OTMuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 17. 17.Caporaso JG, Kuczynski J, Stombaugh J, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7(5):335–336. doi:10.1038/nmeth.f.303 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nmeth.f.303&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20383131&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F30%2F2020.03.29.20039693.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000277175100003&link_type=ISI) 18. 18.DeSantis TZ, Hugenholtz P, Larsen N, et al. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol. 2006;72(7):5069–5072. doi:10.1128/AEM.03006-05 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYWVtIjtzOjU6InJlc2lkIjtzOjk6IjcyLzcvNTA2OSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzAzLzMwLzIwMjAuMDMuMjkuMjAwMzk2OTMuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 19. 19.Zhao L, Zhang F, Ding X, et al. Gut bacteria selectively promoted by dietary fibers alleviate type 2 diabetes. Science (80-). 2018;359(6380):1151–1156. doi:10.1126/science.aao5774 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzNTkvNjM4MC8xMTUxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDMvMzAvMjAyMC4wMy4yOS4yMDAzOTY5My5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 20. 20.Tingley D, Yamamoto T, Hirose K, Keele L, Imai K. mediation: R Package for Causal Mediation Analysis. J Stat Softw. 2014;59(5). [http://www.jstatsoft.org/](http://www.jstatsoft.org/). Accessed March 17, 2020. 21. 21.Oksanen J, Blanchet FG, Kindt R, et al. Package “vegan” Title Community Ecology Package.; 2013. [http://vegan.r-forge.r-project.org/](http://vegan.r-forge.r-project.org/). Accessed March 17, 2020. 22. 22.Segata N, Izard J, Waldron L, et al. Metagenomic biomarker discovery and explanation. Genome Biol. 2011;12(6):R60. doi:10.1186/gb-2011-12-6-r60 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/gb-2011-12-6-r60&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21702898&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F30%2F2020.03.29.20039693.atom) 23. 23.Forouhi NG, Imamura F, Sharp SJ, et al. Association of Plasma Phospholipid n-3 and n-6 Polyunsaturated Fatty Acids with Type 2 Diabetes: The EPIC-InterAct Case-Cohort Study. PLoS Med. 2016;13(7). doi:10.1371/journal.pmed.1002094 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pmed.1002094&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27434045&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F30%2F2020.03.29.20039693.atom) 24. 24.Yary T, Voutilainen S, Tuomainen TP, Ruusunen A, Nurmi T, Virtanen JK. Serum n-6 polyunsaturated fatty acids, Δ5-and D6-desaturase activities, and risk of incident type 2 diabetes in men: The Kuopio Ischaemic Heart Disease Risk Factor Study. Am J Clin Nutr. 2016;103(5):1337–1343. doi:10.3945/ajcn.115.128629 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiYWpjbiI7czo1OiJyZXNpZCI7czoxMDoiMTAzLzUvMTMzNyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzAzLzMwLzIwMjAuMDMuMjkuMjAwMzk2OTMuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 25. 25.Akter S, Kurotani K, Sato M, et al. High Serum Phospholipid Dihomo-γ-Linoleic Acid Concentration and Low Δ5-Desaturase Activity Are Associated with Increased Risk of Type 2 Diabetes among Japanese Adults in the Hitachi Health Study. J Nutr. 2017;147(8):1558–1566. doi:10.3945/jn.117.248997 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6OToibnV0cml0aW9uIjtzOjU6InJlc2lkIjtzOjEwOiIxNDcvOC8xNTU4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDMvMzAvMjAyMC4wMy4yOS4yMDAzOTY5My5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 26. 26.Wan Y, Wang F, Yuan J, et al. Effects of dietary fat on gut microbiota and faecal metabolites, and their relationship with cardiometabolic risk factors: a 6-month randomised controlled-feeding trial. Gut. 2019:gutjnl-2018-317609. doi:10.1136/gutjnl-2018-317609 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiZ3V0am5sIjtzOjU6InJlc2lkIjtzOjk6IjY4LzgvMTQxNyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzAzLzMwLzIwMjAuMDMuMjkuMjAwMzk2OTMuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 27. 27.Ohkawara S, Furuya H, Nagashima K, Asanuma N, Hino T. Oral Administration of Butyrivibrio fibrisolvens, a Butyrate-Producing Bacterium, Decreases the Formation of Aberrant Crypt Foci in the Colon and Rectum of Mice. J Nutr. 2005;135(12):2878–2883. doi:10.1093/jn/135.12.2878 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6OToibnV0cml0aW9uIjtzOjU6InJlc2lkIjtzOjExOiIxMzUvMTIvMjg3OCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzAzLzMwLzIwMjAuMDMuMjkuMjAwMzk2OTMuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 28. 28.Cendrowska M, Kasperowicz A, Michałowski T. Preliminary assessment of the capability of the rumen bacterium, Butyrivibrio fibrisolvens, to utilize fructose polymers for growth. J Anim Feed Sci. 2006;15(Suppl. 1):11–14. doi:10.22358/jafs/70131/2006 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.22358/jafs/70131/2006&link_type=DOI) 29. 29.Inoue R, Ohue-Kitano R, Tsukahara T, et al. Prediction of functional profiles of gut microbiota from 16S rRNA metagenomic data provides a more robust evaluation of gut dysbiosis occurring in Japanese type 2 diabetic patients. J Clin Biochem Nutr. 2017;61(3):217–221. doi:10.3164/jcbn.17-44 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3164/jcbn.17-44&link_type=DOI) 30. 30.Ozato N, Saito S, Yamaguchi T, et al. Blautia genus associated with visceral fat accumulation in adults 20–76 years of age. npj Biofilms Microbiomes. 2019;5(1). doi:10.1038/s41522-019-0101-x [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41522-019-0101-x&link_type=DOI) 31. 31.Konikoff T, Gophna U. Oscillospira: a Central, Enigmatic Component of the Human Gut Microbiota. Trends Microbiol. 2016;24(7):523–524. doi:10.1016/j.tim.2016.02.015 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.tim.2016.02.015&link_type=DOI) 32. 32.Gophna U, Konikoff T, Nielsen HB. Oscillospira and related bacteria – From metagenomic species to metabolic features. Environ Microbiol. 2017;19(3):835–841. doi:10.1111/1462-2920.13658 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/1462-2920.13658&link_type=DOI) 33. 33.Luo G, Li B, Yang C, et al. Major Traditional Probiotics: Comparative Genomic Analyses and Roles in Gut Microbiome of Eight Cohorts. Front Microbiol. 2019;10. doi:10.3389/fmicb.2019.00712 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fmicb.2019.00712&link_type=DOI) 34. 34.Salma MU, Nessa A, Saki SA. Relationship of HDL-C and Triglyceride with Type 2 Diabetes Mellitus in Adult Female. Mymensingh Med J. 2019;28(2):306–310. [http://www.ncbi.nlm.nih.gov/pubmed/31086143](http://www.ncbi.nlm.nih.gov/pubmed/31086143). Accessed December 30, 2019. 35. 35.Lyte JM, Gabler NK, Hollis JH. Postprandial serum endotoxin in healthy humans is modulated by dietary fat in a randomized, controlled, cross-over study. Lipids Health Dis. 2016;15(1):186. doi:10.1186/s12944-016-0357-6 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12944-016-0357-6&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27816052&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F30%2F2020.03.29.20039693.atom) 36. 36.Calle MC, Fernandez ML. Inflammation and type 2 diabetes. Diabetes Metab. 2012;38(3):183–191. doi:10.1016/j.diabet.2011.11.006 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.diabet.2011.11.006&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22252015&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F30%2F2020.03.29.20039693.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000307431200001&link_type=ISI) 37. 37.Ertunc ME, Hotamisligil GS. Lipid signaling and lipotoxicity in metaflammation: Indications for metabolic disease pathogenesis and treatment. J Lipid Res. 2016;57(12):2099–2114. doi:10.1194/jlr.R066514 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamxyIjtzOjU6InJlc2lkIjtzOjEwOiI1Ny8xMi8yMDk5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDMvMzAvMjAyMC4wMy4yOS4yMDAzOTY5My5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 38. 38.Lee JM, Lee H, Kang SB, Park WJ. Fatty acid desaturases, polyunsaturated fatty acid regulation, and biotechnological advances. Nutrients. 2016;8(1). doi:10.3390/nu8010023 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/nu8010023&link_type=DOI) 39. 39.Cani PD, Osto M, Geurts L, Everard A. Involvement of gut microbiota in the development of low-grade inflammation and type 2 diabetes associated with obesity. Gut Microbes. 2012;3(4). doi:10.4161/gmic.19625 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.4161/gmic.19625&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22572877&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F30%2F2020.03.29.20039693.atom) 40. 40.Cani PD, Amar J, Iglesias MA, et al. Metabolic endotoxemia initiates obesity and insulin resistance. Diabetes. 2007;56(7):1761–1772. doi:10.2337/db06-1491 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZGlhYmV0ZXMiO3M6NToicmVzaWQiO3M6OToiNTYvNy8xNzYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDMvMzAvMjAyMC4wMy4yOS4yMDAzOTY5My5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 41. 41.Bäckhed F, Normark S, Schweda EKH, Oscarson S, Richter-Dahlfors A. Structural requirements for TLR4-mediated LPS signalling: a biological role for LPS modifications. Microbes Infect. 2003;5(12):1057–1063. doi:10.1016/s1286-4579(03)00207-7 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1286-4579(03)00207-7&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=14554246&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F30%2F2020.03.29.20039693.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000186306200001&link_type=ISI) 42. 42.Suganami T, Tanimoto-Koyama K, Nishida J, et al. Role of the Toll-like receptor 4/NF-κB pathway in saturated fatty acid-induced inflammatory changes in the interaction between adipocytes and macrophages. Arterioscler Thromb Vasc Biol. 2007;27(1):84–91. doi:10.1161/01.ATV.0000251608.09329.9a [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYXR2YmFoYSI7czo1OiJyZXNpZCI7czo3OiIyNy8xLzg0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDMvMzAvMjAyMC4wMy4yOS4yMDAzOTY5My5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 43. 43.Pendyala S, Walker JM, Holt PR. A high-fat diet is associated with endotoxemia that originates from the gut. Gastroenterology. 2012;142(5). doi:10.1053/j.gastro.2012.01.034 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1053/j.gastro.2012.01.034&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22326433&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F30%2F2020.03.29.20039693.atom) 44. 44.Jenq RR, Taur Y, Devlin SM, et al. Intestinal Blautia Is Associated with Reduced Death from Graft-versus-Host Disease. Biol Blood Marrow Transplant. 2015;21(8):1373–1383. doi:10.1016/j.bbmt.2015.04.016 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.bbmt.2015.04.016&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25977230&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F30%2F2020.03.29.20039693.atom) 45. 45.Chen W, Liu F, Ling Z, Tong X, Xiang C. Human intestinal lumen and mucosa-associated microbiota in patients with colorectal cancer. PLoS One. 2012;7(6):e39743. doi:10.1371/journal.pone.0039743 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0039743&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22761885&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F30%2F2020.03.29.20039693.atom)