Arachidonic acid-derived hydroxyeicosatetraenoic acids are positively associated with colon polyps in adult males: a cross-sectional study

Oxylipids are potent lipid mediators associated with inflammation-induced colon carcinomas and colon tumor survival. Therefore, oxylipid profiles may be useful as novel biomarkers of colon polyp presence. The aim of this study was to investigate the relationship between plasma non-esterified oxylipids and the presence of colon polyps. A total of 123 Caucasian men, ages 48 to 65, were categorized into three groups: those with no polyps, those with one or more hyperplastic polyps, and those with one or more adenomas. Plasma non-esterified oxylipids were analyzed using solid phase extraction and quantified using a targeted HPLC tandem mass spectrometric analysis. Statistical analyses included Kruskal-Wallis one-way ANOVA with Dunn’s test for multiple comparison and generalized linear models to adjust for confounding factors such as age, anthropometrics, and smoking status. In general, monohydroxy omega-6-derived oxylipids were significantly increased in those with polyps. Concentrations of 5-hydroxyeicosatetraenoic acid (HETE) and 11-HETE were significantly higher in those with hyperplastic polyps and adenomas compared to those with no polyps. Arachidonic acid-derived HETEs were significantly associated with colon polyp types, even after adjusting for age, smoking, and body mass index or waist circumference in regression models. Since many of these oxylipids are formed through oxygenation by lipoxygenases (i.e., 5-, 12-, and 15-HETE, and 15- hydroxyeicosatrienoic acid [HETrE]) or auto-oxidative reactions (i.e., 11-HETE), this may indicate that lipoxygenase activity and lipid peroxidation are increased in those with colon polyps. In addition, since oxylipids such as 5-, 12-, and 15-HETE are signaling molecules involved in inflammation regulation, these oxylipids may have important functions in inflammation-associated polyp presence. Future studies should be performed in a larger cohorts to investigate if these oxylipids are useful as potential biomarkers of colon polyps.

Pathology associated with obesity contributes to chronic low-grade inflammation which plays a role in the development and pathogenesis of colon polyps 1 . In obesity, chronic low-grade inflammation is associated with excess lipid accumulation in white adipose tissue (WAT) which triggers monocyte recruitment and activation of M1 macrophages 7,8 . In fact, macrophages are major sources of oxidized lipid mediators deemed oxylipids 9 , which are potent bioactive lipid mediators endogenously synthesized. Oxylipids can be synthesized enzymatically by cytochrome P450 (CYP) enzymes 10 , cyclooxygenases (COX) 11 , and lipoxygenase (LOX) 12 . However, oxylipids are also produced through non-enzymatic autoxidation mechanisms (i.e., free radicals) 9 . Omega-3 (n-3) and omega-6 (n-6) polyunsaturated fatty acids (PUFAs) are structurally similar oxylipid substrates, and accumulation of one PUFA family can lead to enzymatic competition and alter downstream lipid metabolism. Most n-3 PUFA metabolites are considered to be anti-inflammatory while those from n-6 PUFA are considered to be pro-inflammatory in cancer models 9 . Oxylipid function is dependent on the oxidation pathway and the precursor PUFA substrate, and oxylipid signaling mechanisms include angiogenesis, cell growth, and either activating/ enhancing or resolving inflammatory responses 13 .
Western diets are typically high in n-6 PUFAs such as linoleic acid (LA), which can alter the ratio of PUFAs leading to a shift towards n-6 fatty acid elongation and desaturation, resulting in accumulation of arachidonic acid (ARA) and ARA-derived oxylipids. It is suspected that altered lipid metabolism can play a role in the transformation of colorectal polyps to colorectal cancer [14][15][16] . Previous studies investigated PUFA dietary intake effects on adenomas and colon cancer risk, however, the results from these studies are mixed [17][18][19] . Dietary recalls and nutritional assessments are subject to bias from participants, thus, investigating PUFA biochemical markers and metabolites may identify biomarkers associated with colon polyp risk. In our study, 63 plasma oxylipids were quantified to determine oxylipid classes and concentrations significantly associated with colon polyp presence and type in 123 males. After false discovery correction and across all statistical analyses, the plasma ARA-derived oxylipids 5-and 11-hydroxyeciosatetreanoic acid (HETEs) were significantly elevated and associated with polyp type. The effects remained significant even after normalizing these HETEs to ARA concentrations and adjusting for confounding factors such as age, smoking, body mass index (BMI), and waist circumference (WC). In addition, concentration ranges of these HETEs were associated with an increased likelihood of polyp presence.

Methods and Materials ethics statement. The study was approved by the Biomedical and Health Institutional Review Board (IRB)
of Michigan State University (IRB# 08-786). All methods were performed in accordance with the relevant guidelines and regulations of the Michigan State University IRB. At the time of enrollment, immediately prior to routine colonoscopy, written informed consent was obtained from each participant.

Study population and polyp identification.
A group of asymptomatic and >96% Caucasian males ranging from 48-65 years of age were recruited between August 2009 and February 2011 in a cross-sectional study as previously described [6][7][8][20][21][22][23][24] . Exclusion criteria included but not limited to: 1) cancer in the past two years, 2) diabetes, 3) autoimmune diseases, etc. A detailed list of exclusion criteria can be found in the previous published paper 21 . Clinical metadata on participants' co-morbidities, medication status, and family history were obtained at the time of enrollment immediately prior to full colonoscopy. During the colonoscopy, a gastroenterologist (MSU-affiliated clinics, MI) categorized each polyp based on the location and board-certified pathologists assigned a polyp type to each specimen collected. The participants were grouped based on their polyp types as: 1) those with no polyp, 2) those with one or more hyperplastic polyps, and 3) those with one or more adenomas. Anthropometric measures (i.e. height, weight, and waist circumference) of study subjects were taken by trained medical professionals and these anthropometrics measures were used to calculate BMI. Venous blood of participants was drawn and plasma fraction was isolated by centrifugation and stored at −80 °C until analysis. Smoking status was assigned to subjects as "ever smoked" or "never smoked". A complete description of the study design can be found elsewhere 6 .
Non-esterified plasma PUFA and oxylipid extraction, isolation and analysis. Non-esterified plasma PUFA and oxylipid extraction with internal standards and isolation were carried out as previously reported by our group and are described in detail 20,25 . In brief, 500 μL of each patient's plasma was used for lipid extraction. Deuterium labeled internal standards were added. Phenomenex Strata-X (60 mg/3 mL, Phenomenex, Torrane, CA) was used to perform solid-phase extraction for each sample as previously described 7,20 . The prepared samples were stored under high-purity argon at −20 °C for no longer than a week. Liquid chromatography tandem mass spectrometry (LC/MS/MS) analysis was conducted as previously reported 20 . In brief, liquid chromatography separations were carried out on a Waters ACQUITY UPLC system (Waters, Milford, MA) using an Ascentis Express C18 column (10 cm × 2.1 mm; 2.7 μm particles, Sigma-Aldrich, St. Louis, MO) which was connected to a Waters Xevo TQ-S triple quadrupole mass spectrometer (Waters, Milford, MA). Peak detection, integration, and quantification were performed using TargetLynx (Waters) software.
The participants were grouped based on the polyp types as previously mentioned 20,25 . Statistical differences in age, BMI, WC, and non-esterified PUFA and oxylipids concentrations among three groups of participants were determined by Kruskal Wallis one-way analysis of variance (ANOVA) and Dunn's test for multiple comparison. P-values were corrected for false discoveries according to Benjamini-Hochberg (BH) and Bonferroni (Bon) 26,27 . Polyp types (i.e., hyperplastic and adenoma) were regressed on the concentrations of oxylipids and clinical covariates, including age and smoking status of participants, using generalized linear model assuming a Poisson distribution. Model 1 included log-transformed explanatory variables (PUFA, oxylipids, or oxylipid product-to precursor ratios) and was adjusted for age and smoking status. Since BMI and WC are highly correlated, two additional models were run adjusting for age and smoking along with either BMI or WC (Table 2). Next, each individual oxylipid concentration was normalized to its precursor PUFA to determine whether differences in oxylipid concentrations were due to altered PUFA availability or metabolism (Table 3). Finally, simple logistic regression models, adjusted for age, smoking, and BMI, were used to calculate polyp presence likelihood calculate across oxylipids concentrations categorized into tertiles (Table 4). Test for exposure was conducted to determine if increases in oxylipid tertiles compared to the lowest tertile were associated with having polyps. The test for trend was conducted to determine if increase in oxylipid tertiles were associated with polyps compared to individuals without polyps. All the statistical analyses, except logistic regressions, were performed using R v3.2.2 28 . Logistic regressions were performed using SAS version 9.4 (SAS, Cary, NC).

Results
characteristics and oxylipid concentrations of study populations. Age, BMI, WC, smoking status, and non-esterified PUFA and plasma oxylipid concentrations for the overall population (n = 123) are separated into three groups based on the polyp type (Table 1). Of the 123 participants, 67 (54.5%) had no polyps, 20 (16.3%) had hyperplastic polyps, and 36 (29.3%) had adenomas. Non-parametric ANOVAs were carried out for each variable across these three groups. Age did not differ across the three groups. BMI was significantly higher in participants with adenomas. Participants with hyperplastic polyps or adenomas had increased WCs compared to those with no polyps. LA-derived oxylipid concentrations did not differ across the different polyp types. The levels of 15-HETrE, which is a DGLA-derived oxylipid, were higher in individuals with adenomas compared to no polyps. Concentrations of ARA-derived 5-, 11-, and 15-HETE were all increased significantly in participants with hyperplastic polyps and adenomas compared to those with no polyps.
Non-esterified plasma oxylipids are associated with polyp type. Generalized linear models were used to estimate the effects of each PUFA and oxylipid associated with polyp type. Poisson distribution was assumed for polyp type due to the distribution of patients in each group (Table 1), and models were adjusted for age and smoking status. In all regression models, the response variable polyp type was regressed on concentrations of oxylipids and their precursor PUFAs (Table 2 and Fig. 1 top), or PUFA-normalized oxylipid concentrations (Table 3 and Fig. 1 bottom). The −log10(p-values) from these single lipid regressions are presented as a Manhattan plot (Fig. 1), and p-values that were significant after BH FDR adjustment were circled, and p-values significant below the Bonferroni cutoff are displayed above the dashed line. 15-HETrE derived from DGLA, 5-, 11-, 15-HETEs and AEA derived from ARA, and 17,18-DiHETE derived from EPA were significantly associated with polyp type. 5-HETE and 11-HETE remained significant even after FDR p-value adjustment.
PUFAs serve as substrates for oxylipids and are obtained through dietary intake. In order to assess whether differences in these oxylipids are due to altered levels of PUFA availability in the body or altered lipid metabolism, each oxylipid was normalized to its parent PUFA. Polyp type was then regressed on PUFA-normalized oxylipids with models adjusting for age and smoking status ( Fig. 1 bottom and Table 3). ARA-derived 5-, 11-, 12-, and 15-HETE were positively associated with polyp type. Only 5-and 11-HETE remained significant after adjusting p-values for FDR. Results from non-parametric ANOVAs indicated BMI and WC were significantly different across polyp type. Therefore, two additional models were analyzed to adjust for BMI and WC separately to see if any lipids remained associated with polyp type. Model 1 adjusted for age and smoking, model 2 adjusted for age, smoking, and BMI, and model 3 adjusted for age, smoking, and WC (Tables 2 and 3). Across all three models, only ARA-derived 5-and 11-HETE remained significant after adjusting for age, smoking, and either BMI and WC (Table 2). After normalization to parent PUFA concentrations, 5-, and 11-HETE remained significant across all models even after FDR correction (   www.nature.com/scientificreports www.nature.com/scientificreports/ Presence of colon polyps is associated with specific concentrations of oxylipids. Next, odds ratios (OR) were determined using polytomous logistic regressions to determine concentration ranges of 5-and 11-HETE that were associated with presence of 1 or more colon polyps compared to individuals with no colon polyps. The concentrations 5-and 11-HETE were separated into tertiles to provide the specific ranges that were most likely to be associated with the presence of polyps. In all models, polyp presence was regressed on oxylipid tertiles and all data is referenced against the no polyps' group. Models were adjusted for age, smoking, and BMI, and odds ratios and confidence intervals were calculated. Individuals with >1.86 nM 5-HETE were five to six times more likely to have polyps than individuals with ≤ 1.86 nM 5-HETE. Results for 11-HETE were similar to

Discussion
The purpose of this study was to identify PUFA substrates and oxylipids associated with presence of colon polyps. We employed LC-MS/MS to analyze plasma PUFAs and oxylipids of participants. We report associations between higher concentrations of pro-inflammatory oxylipids and the presence of colon polyps. Specifically, we report elevated concentrations of ARA-derived 5-HETE and 11-HETE were associated with an increased likelihood of polyp presence. Associations between 5-and 11-HETE concentrations with polyp type remained significant after normalization to parent PUFA (i.e., ARA) concentrations, accounting for age, smoking, and either BMI or WC, and adjusting p-values for false discoveries. However, non-esterified plasma PUFAs were not associated with polyp in our study. Our results suggest that presence of colon polyps is associated with oxylipids that are metabolized through certain pathways, not entirely due to altered dietary PUFA intake 20 . Cyclooxygenases and lipoxygenases are PUFA oxygenating enzymes. ARA-derived oxylipids produced through COX or LOX pathways are often pro-inflammatory mediators involved in promoting inflammation 29 . Early epidemiology research focused on the COX pathways and investigated the role of ARA metabolism in cancer 30 , and there are fewer reports on the role of LOX pathway in colon polyps and colon cancer. Among the LOX pathways, 5-LOX is closely associated with inflammation and carcinogenesis 29 , and 5-LOX expression is associated with colon adenomas and increased polyp size 31 . Overexpression of 5-LOX has been found in human cancers at different sites including prostate, pancreatic, colon, bladder, esophageal, and testicular cancer [32][33][34][35] . 5-LOX products formed by these cancer cells or other tumor associated cells can upregulate and contribute to tumor growth 36 .
Here we report elevated plasma 5-HETE in patients with at least one or more colon polyps, and in patients with hyperplastic polyps or adenomas. These results are in agreement with the findings from a murine and human studies which reported 5-LOX is up-regulated in colon polyps compared to normal colonic mucosa 31,37 . In our study, the association between elevated 5-HETE and colon polyp type remained significant even after adjustment of BMI and WC. Individuals with elevated 5-HETE concentrations were at least six times more likely to have colon polyps compared to patients with no polyps. Moreover, other ARA-derived oxylipids synthesized through LOX pathways were not significantly associated with any responses from our study. Taken together, our observation may indicate that plasma 5-HETE concentrations could be an early indicator of malignant progression. 11-HETE is an ARA-derived oxylipid synthesized de novo through non-enzymatic oxidation and has been described as a marker of lipid peroxidation 38 . 11-HETE can also be produced enzymatically as a by-product of prostaglandin biosynthesis via the COX pathway 39 . Previous studies indicate increased serum 11-HETE is associated with diseases ranging from coronary events to cancers 20 . To our knowledge, we are the first group to report an association between elevated 11-HETE concentrations and colon polyps. The elevated levels of 11-HETE were observed in patients with hyperplastic polyps and adenomas. One possible explanation is that 11-HETE, generated from enzymatic oxidations through COX pathways, could have some function involved in tissue microenvironment that favors tumor cell growth, proliferation, and metastases 29 . On the other hand, since 11-HETE has also been described as a marker of lipid peroxidation, increased plasma 11-HETE could also be indicative of elevated oxidative stress and increased reactive oxygen species (ROS) 38 . Unfortunately our method did not employ a chiral HPLC column so we were unable to further speculate the potential source of 11-HETE production, however, it is possible that the source 11-HETE (i.e., COX-2 or ROS) could be determined by assessing concentration of 11-HETE enantiomers 40 . Future studies should verify that 11-HETE is associated with colon polyp presence and type in larger more diverse population, and determine the stereochemistry of 11-HETE to determine if www.nature.com/scientificreports www.nature.com/scientificreports/ elevated levels are due to biosynthesis through enzymatic, non-enzymatic, or a combination of both reactions in order better understand the association between 11-HETE and colon polyps.
Interestingly, our group previously reported concentrations of non-esterified 5-HETE and 11-HETE were positively associated with increased BMI and WC, in the same study population as our current study 20 . Our group also previously reported in these study participants that increased BMI was associated with an increased likelihood of adenoma presence compared to those with no polyps 6 . In our current study, we performed two statistical models that were adjusted for BMI and WC, and observed 5-and 11-HETE remained significantly associated with polyp type. Therefore, we speculate it is possible that increased levels of 5-and 11-HETE may serve as indicators of colon polyps independent of obesity in Caucasian males. One limitation of our study is that patients' dietary patterns were not collected at the time of their colonoscopy. Previous epidemiologic studies suggest that high-fat diets are associated with increased incidence of tumors at different organ sites including colon 41,42 . Therefore, it is not possible to calculate associations between PUFA dietary intake and colon polyp types in this study. However, our results were consistent even after normalizing oxylipid concentrations to their parent PUFA, to account for the potential effects of altered levels of PUFA substrates (e.g., through dietary intake). Recently Butler et al. reported that PUFA percentage of total free plasma fatty acids were associated with colon cancer risk 43 . Since there is no longer plasma remaining from our study subjects for a similar follow-up lipid analysis, we Figure 1. Manhattan plots of polyp type regressed on non-esterified plasma polyunsaturated fatty acids, oxylipids, and PUFA-normalized oxylipids. Manhattan plots of the −log10(p-value) for polyp type regressed on log transformed PUFAs, oxylipids, and oxylipid parent-to-precursor ratios individually (i.e., one at a time). Single non-esterified plasma fatty acid, oxylipid, and product to precursor ratio regression models are defined as: polyp type = age + smoking + log(oxylipid) and polyp type = age + smoking + log(oxylipid/PUFA). Each y-axis represents −log10(p-value) for each respective model, and x-axis represent the plasma non-esterified fatty acid, oxylipid, and oxylipid product-to-precursor ratios. Plasma lipids with Benjamini-Hochberg false discovery rate p-values ≤ 0.05 are circled. Oxylipids totals calculated as follows: Total HODE = ∑ 9-HODE + 13-HODE; Total KODE = ∑ 9-KODE + 13-KODE; Total EpOME = ∑ 9,10-EpOME + 12,13-EpOME; Total HETE = ∑ 5-HETE + 11-HETE + 12-HETE + 15-HETE + 20-HETE; Total DHET = ∑ 8,9-DHET + 11,12-DHET + 14,15-DHET; Total DiHETE = ∑ 5,6-DiHETE + 17,18-DiHETE; 9-KO/ HODE = 9-KODE/9-HODE; 13-KO/HODE = 13-KODE/13-HODE; 9,10-Di/EpOME = 9,10-DiHOME/9,10-EpOME; 12,13-Di/EpOME = 12,13-DiHOME/12,13-EpOME; 14,15-DH/EET = 14,15-DHET/14,15-EET. (2019) 9:12033 | https://doi.org/10.1038/s41598-019-48381-0 www.nature.com/scientificreports www.nature.com/scientificreports/ were unable to normalize oxylipids to precursor PUFA percentage of total free plasma fatty acids. Recently, Haid et al. reported that plasma metabolite levels were altered after storage at −80 °C for 5 years 44 , and the oxylipid analysis of our patient samples were conducted 5 years after collection. While degradation across our patient samples would have been consistent since fresh unthawed plasma was used in our analysis, we recognize reported concentrations of oxylipids in our study may be lower than concentrations reported in plasma stored for a shorter period of time. This is an important consideration for follow-up studies targeting the tertile concentration cutoffs reported in our study.

conclusion
In conclusion, we report specific ARA-derived oxylipids produced through distinct oxygenating pathways are highly associated with colon polyp type. The oxylipids 5-and 11-HETE were significantly associated with polyps in all statistical models, even after adjusting for false discoveries and accounting for known risk factors such as age, smoking, and either BMI and WC. Future studies could be carried out in a large cohort to see if these oxylipids levels could be used as potential diagnostic measures for colon polyps, or if they can serve as independent predictors for colon polyp risks.