Impact Factor 6.244 | CiteScore 3.9
More on impact ›

ORIGINAL RESEARCH article

Front. Oncol., 11 May 2018 | https://doi.org/10.3389/fonc.2018.00134

Lipidomic Impacts of an Obesogenic Diet Upon Lewis Lung Carcinoma in Mice

imageSneha Sundaram1, imagePetr Žáček1†, imageMichael R. Bukowski1, imageAaron A. Mehus1, imageLin Yan1 and imageMatthew J. Picklo1,2*
  • 1USDA-ARS Grand Forks Human Nutrition Research Center, Grand Forks, ND, United States
  • 2Department of Chemistry, University of North Dakota, Grand Forks, ND, United States

Metabolic reprogramming of lipid metabolism is a hallmark of cancer. Consumption of a high-fat obesogenic diet enhances spontaneous metastasis using a Lewis lung carcinoma (LLC) model. In order to gain further insights into the mechanisms by which dietary fats impact cancer progression, we conducted a lipidomic analysis of primary tumors originated from LLC from mice fed with a standard AIN93G diet or a soybean oil-based high-fat diet (HFD). Hierarchical clustering heatmap analysis of phosphatidylcholine (PC) lipids and phosphatidylethanolamine (PE) lipids demonstrated an increase in polyunsaturated fatty acids (PUFA)-containing phospholipids and a decrease in monounsaturated fatty acids (MUFA)-containing lipids in tumors from mice fed the HFD. The quantities of 51 PC and 24 PE lipids differed in primary tumors of LLC from mice fed the control diet and the HFD. Analysis of triacylglycerol (TAG) lipids identified differences in 32 TAG (by brutto structure) between the two groups; TAG analysis by neutral loss identified 46 PUFA-containing TAG species that were higher in mice fed with the HFD than in the controls. Intake of the HFD did not alter the expression of the de novo lipogenesis enzymes (fatty acid synthase, acetyl-CoA carboxylase-1, and stearoyl-CoA desaturase-1). Our results demonstrate that the dietary fatty acid composition of the HFD is reflected in the higher order lipidomic composition of primary tumors. Subsequent studies are needed to investigate how these lipidomic changes may be used for targeted dietary intervention to reduce tumor growth and malignant progression.

Introduction

Obesity, resulting from an imbalance between energy intake and energy expenditure, is a leading risk factor for cancer. Excess body fat mass is positively associated with the risk of multiple cancers at different organ sites (1, 2). Furthermore, cancer patients who are overweight or obese are at a greater risk of recurrence and shorter disease-free intervals compared to patients with normal body weight (35). In experimental animals, obesity enhances both primary tumorigenesis (6) and metastasis (7, 8). Multiple studies demonstrate that obesity and type of dietary lipid intake modify lipidomic profiles in humans (9, 10) and laboratory rodents (11, 12). Understanding how the tumorigenic lipid profile responds to obesity and dietary lipid intake may provide insights into novel preventive and therapeutic strategies.

Clinical studies show altered lipid metabolism in cancer patients (1315). Phosphatidylcholines (PC), phosphatidylethanolamines (PE), and sphingomyelins are the most dominant phospholipids (PL) comprising 80% of the cell membranes (16). Patients with breast (17), lung (18), and ovarian cancer (19) exhibit elevated concentrations of PCs and PEs in cancerous tissues compared to adjacent normal tissues. On the other hand, significantly lower concentrations of total PEs in plasma are related to clinical stage and pathologic grade of prostate cancer (20). The levels of ceramides (CER), a dominant member of the sphingomyelin family, are elevated in head and neck squamous cell carcinoma compared to the normal tissue (21). The expression of ceramide synthase mRNAs is higher in cancerous breast tissue than in normal breast tissue (22). Triacylglycerols (TAG), in the form of lipid droplets, are the main constituents of body fat. Increased numbers of lipid droplets in colon adenocarcinoma has been suggested to be active in regulation of cancer pathogenesis (23). It is proposed that these altered cellular lipid profiles and distribution change cellular functions, which may affect the development and progression of cancer (24).

Lewis lung carcinoma (LLC) is a highly malignant carcinoma commonly used in murine studies of metastasis (7, 8), angiogenesis (25), and cachexia (26). When LLC cells are subcutaneously transplanted, they produce a rapidly growing primary tumor that metastasizes to the lungs. We found that obesity enhances the aggressiveness of pulmonary metastasis of LLC (7, 8). The purpose of this study was to determine changes in lipidomic signatures in primary tumors of LLC from mice caused by the intake of a high-fat, obesogenic diet.

Materials and Methods

Chemicals

Chloroform, potassium chloride, butylated hydroxytoluene (BHT), ammonium chloride, and ammonium acetate were purchased from Sigma-Aldrich (St. Louis, MO, USA). Methanol was from Avantor Performance Materials, Inc. (Center Valley, PA, USA). PL and CER standards, including LIPID MAPS mass spectrometry (MS) standards—Core H, were from Avanti Polar Lipids Co. (Alabaster, AL, USA). The compounds PC(17:0/14:1), PC(17:0/20:4), PC(21:0/22:6), LPC(17:1), PE(17:0/14:1), and PE(21:0/22:6) of certified concentrations were used as the internal standards for quantitation of PCs, lyso PCs (LPCs), CERs, PEs, and lyso PEs (LPEs). TAG species were from Nu-Chek Inc. (Elysian, MN, USA). Antibodies against fatty acid synthase (FASN), acetyl-CoA carboxylase-1 (ACC1), stearoyl-CoA desaturase-1 (SCD1) were from Cell Signaling Technology (Danvers, MA, USA), that against phosphorylated ACC1 (pACC) was from EMD Millipore (Burlington, MA, USA), and that against the system control protein was from ProteinSimple (San Jose, CA, USA).

Animals and Diets

Three-week-old male C57BL/6 mice (Harlan, Madison, WI, USA) were maintained in a pathogen-free room with a temperature of 22 ± 1°C and a 12:12 h light/dark cycle. The standard AIN93G diet (27) providing 16% of energy from soybean oil and a modified AIN93G diet providing 45% of energy from soybean oil [high-fat diet (HFD)] were used in this study (Table 1). Both diets were powder diets and were stored at −20°C until feeding.

TABLE 1
www.frontiersin.org

Table 1. Composition of experimental diets.

Lewis Lung Carcinoma

The LLC cell line, a variant that metastasizes to lungs (28), was obtained from Dr. Pnina Brodt, McGill University, Montreal, Quebec, Canada. The cells were cultured with RPMI-1640 medium containing 10% heat-inactivated fetal bovine serum and maintained in a humidified atmosphere of 5% CO2 in air at 37°C. Hoechst DNA staining and direct culture tests (performed by American Type Cell Collection, Manassas, VA, USA) showed that cells were free of mycoplasma.

Experimental Design

After acclimation with the AIN93G diet for 1 week, mice were randomly assigned into two groups (n = 6 per group) and fed the AIN93G diet or the HFD. Mice had free access to their diets and deionized water. Ten weeks after initiation of the experimental feeding, mice were subcutaneously injected with 2.5 × 105 viable LLC cells per mouse into the lower dorsal region. Intake of the HFD for 10 weeks was sufficient to demonstrate an elevation in adiposity as observed in our previous work (29). Mice were assessed for body composition 1 week before cell injection by using an Echo Whole-Body Composition Analyzer (Model 100, Echo Medical System, Houston, TX, USA). Ten days after cell injection when the resulting subcutaneous tumor was approximately 1 cm in diameter, mice were anesthetized with a mixture of ketamine and xylazine; the subcutaneous tumors were harvested and cleaned to remove skin, weighed, flash-frozen in liquid nitrogen, and stored at −80°C.

Lipidomic Analysis

The primary tumor (approximately 100 mg) was homogenized with 2 mL of extraction solution (hexane:isopropanol 1:1, 50 µM BHT) and centrifuged at 2,000 g for 10 min at 10°C. The liquid phase was collected. The process was repeated once. The pooled liquid phases were evaporated under an N2 stream at 30°C. The dried extract was reconstituted in 1 mL of chloroform:methanol (1:1) to yield a primary extract. The extract (10 μL) was mixed with 980 μL of mobile phase (chloroform: methanol 1:1, 10 mM ammonium acetate, 50 μM of BHT) and 10 μL of a CER(d18:1/17:0) to yield a secondary extract for MS analysis. This analysis determined isobaric and isomeric PCs using MS3 fragmentation (LPC fragments were employed), isobaric PEs using MS2 scan (fatty acid fragments were employed), and CER analysis using precursor ion scan (PIS) m/z +264. The secondary extract (25 µL) was further mixed with 975 µL of mobile phase and spiked with 15 µL of a mixture of PC and LPC internal standard and 8 µL of PE internal standard for analysis of PC, LPC, PE, and LPE using experiments PIS m/z +184 and neutral loss scan (NLS) m/z +141, respectively. The mass spectrometric analysis for PCs and PEs were performed using the methods previously described (30, 31). In addition to the positive mode, PEs were analyzed using the MS2 fragmentation of precursor [M-H] performed in negative mode (32); this combination enabled identification and quantitation of isobaric PE. Furthermore, analysis of the isobaric PE employing MS2 was corrected for possible interferences of PC species. LPEs were quantified only relatively with respect to the PE(17:0/14:1) due to the lack of commercially available LPE internal standards. The identified lipids were normalized by the fresh tumor weight (nmol/g tissue). PL nomenclature was used according to Liebisch et al. (33).

Ceramides were analyzed by using the method of PIS of m/z 264 in positive mode (PIS m/z +264) (34). The scanned fragment m/z 264 resulted from the loss of N-linked fatty acids and two molecules of water. The nature of this fragment that contained the sphingosine backbone enabled direct identification of the amide-linked fatty acids from the mass of the molecular ion [M + H]+ or [M + H–H2O]+. For the quantification purposes only the [M + H]+ ion was employed. Parameters for the experiment were set as follows: collision energy 33 V, ion spray voltage 5,200 V, declustering potential 80 V, entrance potential 10 V, and scanning rate 200 Da/s. CER(17:0) was used for quantification of CER. The following CER including the CER(17:0) were employed for an ionization efficiency correction: CER(10:0), CER(20:0), and CER(24:1). Data collection and processing was performed as previously described (31). Along with CER species, PIS m/z +264 detected hexosylceramides (HexCer), which were identified by precursor mass. These were quantified following the same method as other CER, resulting in semi-quantitative data that allows for detection of differences between treatment groups.

Triacylglycerol Analyses

Mass spectrometry was performed for analysis of TAG as previously described (31). The ion source was calibrated with an equimolar TAG standard as previously shown (35). Neutral mass losses were used to determine the relative content of selected fatty acids across brutto structures in multiplexed experiments, and correction factors for other brutto structures were imputed as detailed previously (31).

Triacylglycerol appeared in the mass spectrum as ammonium complexes between m/z 750 and m/z 1,000, with cluster of TAG species around m/z 825, 850, 875, 900, 925, and 950, corresponding to TAG species with fatty acyl chains containing 48, 50, 52, 54, 56, and 58 carbons, respectively. Within these clusters, signals for the monoisotopic species differed in m/z by two units, representing the desaturation level across the three acyl chains.

A numerical expression for the distribution of saturation within a given TAG group with C carbons as the summation of the product of the saturation number (N) and the mole fraction of the species with that saturation value for saturation values within the group. The resultant weighted desaturation index (DC) describes the average desaturation level of the group, and is higher for more desaturated species.

DC=i=0NmaxNicNi(cN0+cN1+cNmax)

Analysis of De Novo Lipogenesis Enzymes

The primary tumor (approximately 35 mg) was pulverized and lysed for protein extraction as previously described (36). Protein concentration was determined by using the Bradford protein assay (Bio-Rad Laboratories, Inc., Hercules, CA, USA).

Diluted protein lysates were combined with Simple Western™ sample dilution buffer (ProteinSimple, San Jose, CA, USA) which contains reducing agent dithiothreitol, fluorescent standards, and a system control protein [26 or 90 kDa-dependent on the molecular weight (MW) kit used]. A ratio of 4:1 (protein lysate-to-sample dilution buffer) was used. The samples were then analyzed on the Sally Sue instrument (ProteinSimple) as previously reported (37). The ProteinSimple system control protein served as an internal control and was also used to normalize protein expression. For the lower MW SCD1, the 26 kDa system control protein was used for normalization and for the higher MW proteins (FASN, ACC1, and pACC), the 90 kDa system control protein was used.

Statistical Analyses

Resulting lipidomic data were analyzed by using MetaboAnalyst software (version 3.0, McGill University, Sainte Anne de Bellevue, QC, Canada) (38, 39). A false discovery rate (FDR) of 0.05 was used and the FDR-corrected p-values are reported (39). Student’s t-test was performed to compare differences in weight gain and body composition between the two groups by using SAS software (version 9.4, SAS Institute, Cary, NC, USA). Data are presented as mean ± SD. Differences with a p ≤ 0.05 were considered significant.

Results

Morphometric results are presented in Table 2. Mice on the HFD weighed more than those on the AIN93G control diet. Weight gain was 35% higher in mice fed the HFD than in those fed the control diet (Table 2). The percent body fat mass was 44% higher in mice fed the HFD, but there was no difference in absolute lean mass weight between the two groups (Table 2).

TABLE 2
www.frontiersin.org

Table 2. Weight gain and body composition of mice fed the AIN93G or the high-fat diet.

Phospholipids

Analysis of PLs identified 169 PC lipids (discriminating isobaric and regioisomeric PC, LPC, PC ethers, and sphingomyelins) and 81 PE lipids (discriminating isobaric species). Values and statistical comparisons for all PC and PE lipids analyzed are provided in Tables S1 and S2 in Supplementary Material.

Of the 169 PC species identified, 51 were significant in quantity between the two groups (29 were higher and 22 were lower in mice fed the HFD). A volcano plot (a combination of fold change and p-values of comparisons) of PC lipids showed that the difference of eight of the identified PC were at least twofold (Figure 1A); six [PC(18:2/20:2), PC(20:2/18:2), PC(18:3/18:2), PC(18:2/18:3), PC(18:2/18:2), and PC(16:2/18:2)] showed an increase (Figures 1B,C,E,F,G,I) and two [PC(20:1_18:1) and PC(16:1/16:1)] showed a decrease (Figures 1D,H) by HFD. The identity of the 16:2 fatty acid detected in PC(16:2/18:2) was confirmed by GC-MS using authentic standards as the n6 species, 7,10-hexadecadienoic acid.

FIGURE 1
www.frontiersin.org

Figure 1. Volcano plot (A) and concentrations of selected phosphatidylcholines [PCs; (B–I)] in the primary tumor of Lewis lung carcinoma with a twofold or greater difference between the AIN93G diet and the high-fat diet. Values are mean ± SD (n = 6 per group). *p ≤ 0.05 for each comparison.

Hierarchical clustering heatmap analysis makes the identified lipids visualized and shows the variation of each lipid between the comparisons. The heatmap of PC (Figure 2) showed greater amounts of polyunsaturated fatty acids (PUFA) in tumors from mice fed the HFD than in those fed the control diet, whereas greater amounts of monounsaturated fatty acids (MUFA) were in mice fed the control diet than in those fed the HFD.

FIGURE 2
www.frontiersin.org

Figure 2. Hierarchical clustering heatmap analysis of phosphatidylcholine (PC) lipids in the primary tumor of Lewis lung carcinoma from mice fed the AIN93G diet or the high-fat diet. Each colored cell on the map corresponds to a concentration value. The top 25 lipid species from the two diet groups are presented and ranked by t-tests to retain the most contrasting patterns. Values are measured by Euclidean distance with a Ward clustering algorithm (n = 6 per group). *p ≤ 0.05 for each comparison.

We tested the hypothesis that intake of the HFD changed the regioisomeric (sn1 vs sn2) distribution of PC lipids in the tumors. Of the 47 regioisomers analyzed, the percent of the regioisomer in the sn1 vs sn2 form was lower in HFD vs the control for PC(16:0/18:3) (63.5 ± 3.9 vs 49.1 ± 2.7%), PC(16:0/18:2) (80.0 ± 1.7 vs 76.3 ± 1.0%), and PC(16:0/18:1) (59.1 ± 1.1 vs 54.5 ± 1.4%) (Figure 3). These results indicate that for these three PC species, the unsaturated fatty acid was observed in the sn1 position to a greater extent in tumors from mice fed the HFD.

FIGURE 3
www.frontiersin.org

Figure 3. Changes in phosphatidylcholine (PC) regioisomers in primary tumors of Lewis lung carcinoma in response to consumption of high-fat diet. Data are presented as the % of sn1 isomer for each lipid (n = 6 per group). *p ≤ 0.05 for each comparison.

Twenty-four PE are significantly different in quantity between the two groups [12 were higher and 12 were lower in mice fed the HFD (Table S2 in Supplementary Material)]; however, unlike for PC, regioisomers (i.e., sn1 vs sn2) for PEs could not be determined. The volcano plot of PE identified difference in PE(18:2_20:2) (Figure 4A); the amount of PE(18:2_20:2) was twofold greater in mice fed the HFD than in the control mice (Figure 4B). Similar to the PC lipids, there was a shift to an elevation in PUFA-containing PE lipids in tumors from mice fed the HFD as demonstrated by heatmapping with cluster analysis (Figure 5).

FIGURE 4
www.frontiersin.org

Figure 4. Volcano plot (A) and concentration of phosphatidylethanolamine [(PE); 18:2_20:2] (B) in the primary tumors of Lewis lung carcinoma with greater than twofold difference between the AIN93G diet and the high-fat diet. Values are mean ± SD (n = 6 per group). *p < 0.05 for each comparison.

FIGURE 5
www.frontiersin.org

Figure 5. Hierarchical clustering heatmap analysis of phosphatidylethanolamine (PE) lipids in the primary tumor of Lewis lung carcinoma from mice fed the AIN93G diet or high-fat diet. Each colored cell on the map corresponds to a concentration value. The top 25 lipid species from the two diet groups are presented and ranked by t-tests to retain the most contrasting patterns. Values are measured by Euclidean distance with a Ward clustering algorithm (n = 6 per group). *p ≤ 0.05 for each comparison.

Ceramides and HexCer

Changes in CER and HexCer metabolism have been studied with respect to cancer development (21, 22). Our analyses identified and quantified 15 separate CER species and 5 HexCer species. However, there were no significant differences in CER and HexCer species between the two dietary groups (data not shown). Values and statistical comparisons for CER and HexCer analyzed are provided in Table S3 in Supplementary Material.

Triacylglycerol

Analysis of TAG composition identified 62 separated TAG by brutto structure. The volcano plot analysis of TAGs showed that the difference of 15 of the identified TAGs were at least twofold (Figure 6A). Twelve TAGs (53:4, 54:4, 54:5, 54:6, 54:7, 56:5, 56:6, 56:7, 56:8, 58:5, 58:9, and 58:10) showed significant increases and three (46:1, 48:1, and 48:2) showed significant decreases by the HFD (Figure 6B). Of the 62 TAG structures identified, 32 TAG differed between treatment groups and were significant in quantity between the two diets. Values and statistical comparisons for the TAG concentration (for brutto structure) are provided in Table S4 in Supplementary Material. One of the samples from the control, AIN93G group was omitted from statistical analysis for TAG as its values were 2.5-fold greater than that of the others in the dietary group. Analysis of this sample was repeated with similar results.

FIGURE 6
www.frontiersin.org

Figure 6. Volcano plot with a twofold difference (A) and concentrations of triacylglycerols [(TAG); (B)] in the primary tumor of Lewis lung carcinoma from mice fed the AIN93G diet or the high-fat diet (HFD). Values are mean ± SD (n = 5 for the AIN93G group and n = 6 for the HFD group). *p ≤ 0.05 for each comparison.

Enhanced mass spectra for the four most intense TAG clusters (TAG48, TAG50, TAG52, and TAG54) are shown for the control and the HFD groups. The intensity of signals for TAG48 and TAG50 decreased in the HFD group with TAG48:1 and 48:2 having most significant contribution to this decrease (Figure 7A). Signals for TAG52 remained unchanged; however, TAG54 increased in the HFD group, with TAG54:4, 54:5, 54:6, and 54:7 providing significant contributions (Figure 7B). These results can be expressed in terms of the weighted saturation index (DC). For the TAG46, TAG48, and TAG50, DC was unchanged between the control and HFD groups (Table 3). For TAG species with 52 or more carbons, the DC increased between the control and HFD groups, indicating an increase in the concentration of TAG with higher levels of desaturation.

FIGURE 7
www.frontiersin.org

Figure 7. Enhanced mass spectrum for triacylglycerols (TAG) 48 and TAG50 (A) and TAG52 and TAG54 (B). Note the higher concentrations of C48 and polyunsaturated C54 containing TAG species in tumors from the AIN93G (A) and high-fat diet (B) fed mice, respectively. *denote significant differences in TAG species as determined by volcano plot analysis (Figure 6).

TABLE 3
www.frontiersin.org

Table 3. Weighted saturation index (DC) of triacylglycerols (TAG) in tumors of Lewis lung carcinoma from mice fed the AIN93G or the high-fat diet.

While the Dc for the TAG48 and TAG48 clusters do not change between the two groups, they do decrease in concentration. By total TAG concentration for each group, the tumors from HFD-fed mice contained less than 50% of C46 and C48 when compared to tumors from the controls, but had an increase of 159 and 310% of species containing 52 or more carbons thereby contributing to the overall increase in the TAG concentration (Table 3).

Using neutral loss experiments, it was possible to interpret the differences in TAG concentration between the two groups in terms of specific fatty acid moieties. Decreases in the concentration for TAG48 species corresponded to statistically significant decreases in the same species as observed in neutral loss experiments for both palmitic (16:0) or palmitoleic (16:1) acid (Figure 8A). The only 16:0-containing species to experience significant increases in concentration were TAG52:4 and TAG52:5. When the 16:0 contributions are subtracted from these brutto structures, the remaining acyl chain C and saturation values are 36:4 and 36:5, which correspond to fatty acid combinations of 18:2_18:2 and 18:2_18:3, respectively (Figure 8B). In neutral loss experiments for 18:1, 18:2, 18:3, 20:4, and 22:6, species with the 18:2_18:2 and 18:2_18:3 remainder motifs experienced statistically significant increases in the HFD-fed mice. Values and statistical comparisons for the TAG concentration (for NLS determinations) are provided in Table S5 in Supplementary Material. The increase in these linoleic and linolenic acid-containing species is a reflection of the increased concentration of these fatty acids in the HFD.

FIGURE 8
www.frontiersin.org

Figure 8. Relative concentrations of triacylglycerols (TAG) corresponding to the fatty acids that contribute to significant (p ≤ 0.05) concentration changes including16:0 and 16:1 (A) and subtracted from brutto structure including 36:4 and 36:5 (B) in the primary tumor of Lewis lung carcinoma from mice fed the AIN93G diet or the high-fat diet (HFD). Values are mean ± SD (n = 5 for the AIN93G group and n = 6 for the HFD group).

De Novo Lipogenesis

Available studies indicate that the de novo lipogenesis pathway proteins are involved with HFD-induced tumorigenicity (40). However, analysis of the de novo lipogenic proteins demonstrated no differences in expression levels of de novo lipogenic enzymes ACC1, FASN, and SCD1, nor the pACC/ACC1 ratio, between the two groups (Figure 9).

FIGURE 9
www.frontiersin.org

Figure 9. Expression of de novo lipogenesis enzymes fatty acid synthase (FASN), acetyl-CoA carboxylase-1 (ACC1), phosphorylated ACC1 (pACC), stearoyl-CoA desaturase-1 (SCD1), and the ratio of pACC to ACC1. The chemiluminescent signals of the target proteins were normalized to the signals of the system control protein. Values are mean ± SD (n = 6 per group).

Discussion

The alteration of lipid metabolism in cancer cells may affect cancer progression, as evidenced by results demonstrating that lipid profiles of cancer tissues are different from those of adjacent normal tissues (1315). The composition and concentrations of fatty acids of a diet may alter lipid profiles of cancer, and thus determine the fate of the diet on tumorigenesis (11). LLC is an aggressive carcinoma, which metastasizes spontaneously to the lungs from a subcutaneously transplanted primary tumor (28, 41); feeding mice an obesogenic, HFD enhances its pulmonary metastasis (7, 41). In the present study, consumption of the HFD resulted in quantitative changes in PLs and TAG in the primary tumor of LLC.

Significant alterations in cellular glycerophospholipids occurred in the tumors following intake of HFD. Malignant progression alters the cellular composition of glycerophospholipids and may result in enhanced malignancy as a result of changes in cellular structure and signaling. PC and PE lipids in cancer tissue are altered compared to normal tissue as a result of malignant progression in humans (lung cancer and breast cancer) (17, 18), in rodent models of cancer (P53 and Brca1 double-knockout), and in in vitro works (21MT-2 metastatic breast cancer cells vs 76N normal epithelial cells; hormone-resistant MB435 and MB231 breast cancer cells vs hormone-sensitive MCF-7 breast cancer cells) (42, 43).

Our results demonstrate that intake of HFD decreased the content of MUFA and saturated fatty acid (SFA) containing PC and PE lipids in tumors relative to those from normal fat, lean controls; where as PUFA-containing PE and PC lipids were elevated (see Figures 2 and 5). Soybean oil, rich in the n-6 and n-3 PUFAs α-linolenic acid (ALA; 18:3n-3) and linoleic acid (LA; 18:2n-6) respectively, was used in this study as the source of dietary fat as part of the AIN93 formulation (27). Thus, the total amount of fatty acids, but not the proportions of the fatty acids, differed between the dietary groups. The shift in PE and PC composition to a higher PUFA composition in tumors from HFD fed mice is reflective of higher total PUFA intake. The decrease in SFA and MUFA containing PE and PC lipids may be the result of (1) greater selectivity for PUFA vs SFA and MUFA in these PLs, (2) selective shunting of MUFA and SFA to other lipid pools, and/or (3) a decrease in de novo synthesis of SFA and MUFA either in tumors or in the mouse.

The infusion MS approach we used is able to determine sn1 vs sn2 position of fatty acids in PC esters. For most PC esters, the percentage of fatty acid in sn1 vs sn2 is relative stable. For example, PC18:0/18:2 is observed a 90% of the total PC18:0_18:2 species, a percentage that is not altered by intake of the HFD (Table S1 in Supplementary Material). On the other hand, the sn1 regioisomer percentage decreased for palmitic acid-containing PC(16:0/18:3), PC(16:0/18:1), and PC(16:0/18:2). Except for PC(16:0/18:1), the decrease in the sn1 percentage was due to an increase in the content of corresponding sn2 regioisomer. The biochemical impact of these changes is not clear. However, such changes in PC structures may impact membrane fluidity or interactions with proteins.

Cancer is metabolically active to provide energy for its rapid development and growth. TAG stored in lipid droplets of the cells, provide a reservoir of fatty acids that can be utilized for energy generation. TAG species are elevated in non-small cell lung cancer compared to the normal lung tissue (14). In this study, the HFD altered TAG composition and concentration in the tumor, evidenced by an increase in the overall concentrations of TAG in the tumor. These findings indicate that consumption of a diet high in PUFAs leads to the incorporation of dietary fatty acids into the TAG species and suggest that exogenous lipids may alter the energy metabolism of LLC cells. Furthermore, the fatty acid reservoir has been suggested to not only serve as an energy depot, but may also have active roles in cancer pathogenesis; cancer promoting cyclooxygenase-2 and prostaglandin synthase are found in lipid droplets of colon cancer cells (23). It remains to be investigated whether the elevation of TAG is purely an increase in energy storage, our findings suggest that this elevation may contribute to the enhanced aggressiveness of LLC.

Several reports indicate that CER metabolism is implicated in tumor development in part through extracellular matrix interactions (44, 45). However, we observed no differences in the concentrations of several CER or HexCer species in tumors as a result of intake of the HFD. While we recognize that this lack of change may be the result of the small number of tumors analyzed, the CER and HexCer concentrations were similar between both groups as opposed to the changes for PC, PE, and TAG. We cannot exclude that this lack of change results from the intake of a high PUFA, obesogenic diet vs a high SFA, obesogenic diet. Intake of a high SFA diet is associated with increases in CER (46, 47).

Several studies describe the ability of tumors to become more lipogenic, thereby giving them the ability to provide their own fatty acids essential for promoting tumorigenesis (40). The expression of de novo lipogenic enzymes in this study was not altered by the diets; it suggests that generation of fatty acids through lipogenesis is not likely responsible for the modified fatty acid levels observed. It has been shown that a HFD reduces hepatic de novo lipogenesis (48). It is within reason that the HFD may cause liver to reduce production of MUFA and SFA and that the tumor fatty acid profile may reflect hepatic production of MUFA and SFA. We cannot rule out that decreases of MUFA and SFA in TAG may be the result of lipolytic metabolism of MUFA and SFA containing TAG or perhaps modifications in fatty acid transport in the tumor (40, 49). Subsequent tracer-based studies are needed to define mechanisms underlying changes in lipid metabolic pathways.

Our work characterizes the lipidomic profile of LLC in control and obese mice and demonstrated that an obesogenic, HFD altered the lipidome. This alteration reflected the fatty acid composition of soybean oil used in the diet. We demonstrated previously that intake of a HFD enhances lung metastasis from a subcutaneous primary tumor of LLC (8, 50), indicating that HFD increases the aggressiveness of malignant cells. The changes in lipidomic profiles of LLC by the HFD suggest that this alteration may contribute, at least partly, to the enhanced metastasis. Lipidomic research provides us a useful tool in identifying and quantifying thousands of cellular lipid species and their interactions with other nutrients and metabolites. Further studies are warranted to define the biochemical mechanisms of obesity in cancer promotion, so as to the strategies for its prevention, through identifying alterations in cellular lipid metabolism, trafficking, and homeostasis.

Ethics Statement

This study was performed in accordance with the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health (51) and was approved by the Institutional Animal Care and Use Committee of Grand Forks Human Nutrition Research Center.

Author Contributions

MP, LY, and SS contributed to the conception and design of the study, data interpretation, and writing of the manuscript. PŽ and MB performed lipidomic analyses. AM performed de novo lipogenesis analyses, and they contributed to data interpretation and manuscript preparation. All authors contributed to review and revision of the manuscript and agreed to be accountable for the content of the work.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The authors acknowledge the assistant of the following staff of Grand Forks Human Nutrition Research Center: Lana DeMars and Kay Keehr for technical support, James Lindlauf for preparing diets, and vivarium staff for providing high quality animal care.

Funding

This work was supported by National Institute of Food and Agriculture, U.S. Department of Agriculture, under award number 2014-67017-21758 and through U.S. Department of Agriculture, Agricultural Research Service, research projects 3062-51000-053-00D and 3062-51000-050-00D.

Supplementary Material

The Supplementary Material for this article can be found online at https://www.frontiersin.org/articles/10.3389/fonc.2018.00134/full#supplementary-material.

Abbreviations

ACC1, acetyl-CoA carboxylase-1; BHT, butylated hydroxytoluene; CER, ceramides; FASN, fatty acid synthase; FDR, false discovery rate; HexCer, hexosylceramide; HFD, high-fat diet; LLC, Lewis lung carcinoma; LPC, lysophosphatidylcholine; LPE, lysophosphatidylethanolamine; MS, mass spectrometry; MUFA, monounsaturated fatty acid; NLS, neutral loss scan; pACC, phosphorylated ACC1; PC, phosphatidylcholine; PE, phosphatidylethanolamine; PIS, precursor ion scan; PL: phospholipids; PUFA, polyunsaturated fatty acid; SCD1, stearoyl-CoA desaturase-1; SFA, saturated fatty acid; sPLS-DA, sparse partial least squares-discriminant analysis; TAG, triacylglycerol.

References

1. Lauby-Secretan B, Scoccianti C, Loomis D, Grosse Y, Bianchini F, Straif K, et al. Body fatness and cancer – viewpoint of the IARC working group. N Engl J Med (2016) 375:794–8. doi:10.1056/NEJMsr1606602

CrossRef Full Text | Google Scholar

2. Picon-Ruiz M, Morata-Tarifa C, Valle-Goffin JJ, Friedman ER, Slingerland JM. Obesity and adverse breast cancer risk and outcome: mechanistic insights and strategies for intervention. CA Cancer J Clin (2017) 67(5):378–97. doi:10.3322/caac.21405

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Daniell HW. Increased lymph node metastases at mastectomy for breast cancer associated with host obesity, cigarette smoking, age, and large tumor size. Cancer (1988) 62:429–35. doi:10.1002/1097-0142(19880715)62:2<429::AID-CNCR2820620230>3.0.CO;2-4

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Bassett WW, Cooperberg MR, Sadetsky N, Silva S, Duchane J, Pasta DJ, et al. Impact of obesity on prostate cancer recurrence after radical prostatectomy: data from CaPSURE. Urology (2005) 66:1060–5. doi:10.1016/j.urology.2005.05.040

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Moon HG, Ju YT, Jeong CY, Jung EJ, Lee YJ, Hong SC, et al. Visceral obesity may affect oncologic outcome in patients with colorectal cancer. Ann Surg Oncol (2008) 15:1918–22. doi:10.1245/s10434-008-9891-4

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Sundaram S, Yan L. High-fat diet enhances mammary tumorigenesis and pulmonary metastasis and alters inflammatory and angiogenic profiles in MMTV-PyMT mice. Anticancer Res (2016) 36:6279–87. doi:10.21873/anticanres.11223

CrossRef Full Text | Google Scholar

7. Sundaram S, Yan L. Dietary energy restriction reduces high-fat diet-enhanced metastasis of Lewis lung carcinoma in mice. Oncotarget (2016) 7:65669–75. doi:10.18632/oncotarget.11598

CrossRef Full Text | Google Scholar

8. Yan L, Sundaram S. Monocyte chemotactic protein-1 deficiency reduces spontaneous metastasis of Lewis lung carcinoma in mice fed a high-fat diet. Oncotarget (2016) 7:24792–9. doi:10.18632/oncotarget.8364

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Pietilainen KH, Rog T, Seppanen-Laakso T, Virtue S, Gopalacharyulu P, Tang J, et al. Association of lipidome remodeling in the adipocyte membrane with acquired obesity in humans. PLoS Biol (2011) 9:e1000623. doi:10.1371/journal.pbio.1000623

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Tonks KT, Coster AC, Christopher MJ, Chaudhuri R, Xu A, Gagnon-Bartsch J, et al. Skeletal muscle and plasma lipidomic signatures of insulin resistance and overweight/obesity in humans. Obesity (Silver Spring) (2016) 24:908–16. doi:10.1002/oby.21448

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Djuric Z, Aslam MN, Simon BR, Sen A, Jiang Y, Ren J, et al. Effects of fish oil supplementation on prostaglandins in normal and tumor colon tissue: modulation by the lipogenic phenotype of colon tumors. J Nutr Biochem (2017) 46:90–9. doi:10.1016/j.jnutbio.2017.04.013

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Wang W, Yang J, Qi W, Yang H, Wang C, Tan B, et al. Lipidomic profiling of high-fat diet-induced obesity in mice: importance of cytochrome P450-derived fatty acid epoxides. Obesity (Silver Spring) (2017) 25:132–40. doi:10.1002/oby.21692

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Braicu EI, Darb-Esfahani S, Schmitt WD, Koistinen KM, Heiskanen L, Poho P, et al. High-grade ovarian serous carcinoma patients exhibit profound alterations in lipid metabolism. Oncotarget (2017) 8:102912–22. doi:10.18632/oncotarget.22076

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Eggers LF, Muller J, Marella C, Scholz V, Watz H, Kugler C, et al. Lipidomes of lung cancer and tumour-free lung tissues reveal distinct molecular signatures for cancer differentiation, age, inflammation, and pulmonary emphysema. Sci Rep (2017) 7:11087. doi:10.1038/s41598-017-11339-1

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Ros-Mazurczyk M, Jelonek K, Marczyk M, Binczyk F, Pietrowska M, Polanska J, et al. Serum lipid profile discriminates patients with early lung cancer from healthy controls. Lung Cancer (2017) 112:69–74. doi:10.1016/j.lungcan.2017.07.036

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Fajardo VA, Mcmeekin L, Leblanc PJ. Influence of phospholipid species on membrane fluidity: a meta-analysis for a novel phospholipid fluidity index. J Membr Biol (2011) 244:97–103. doi:10.1007/s00232-011-9401-7

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Hilvo M, Denkert C, Lehtinen L, Muller B, Brockmoller S, Seppanen-Laakso T, et al. Novel theranostic opportunities offered by characterization of altered membrane lipid metabolism in breast cancer progression. Cancer Res (2011) 71:3236–45. doi:10.1158/0008-5472.CAN-10-3894

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Lee GK, Lee HS, Park YS, Lee JH, Lee SC, Lee JH, et al. Lipid MALDI profile classifies non-small cell lung cancers according to the histologic type. Lung Cancer (2012) 76:197–203. doi:10.1016/j.lungcan.2011.10.016

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Kang S, Lee A, Park YS, Lee SC, Park SY, Han SY, et al. Alteration in lipid and protein profiles of ovarian cancer: similarity to breast cancer. Int J Gynecol Cancer (2011) 21:1566–72. doi:10.1097/IGC.0b013e318226c5f5

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Cvetkovic B, Vucic V, Cvetkovic Z, Popovic T, Glibetic M. Systemic alterations in concentrations and distribution of plasma phospholipids in prostate cancer patients. Med Oncol (2012) 29:809–14. doi:10.1007/s12032-011-9914-z

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Karahatay S, Thomas K, Koybasi S, Senkal CE, Elojeimy S, Liu X, et al. Clinical relevance of ceramide metabolism in the pathogenesis of human head and neck squamous cell carcinoma (HNSCC): attenuation of C(18)-ceramide in HNSCC tumors correlates with lymphovascular invasion and nodal metastasis. Cancer Lett (2007) 256:101–11. doi:10.1016/j.canlet.2007.06.003

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Erez-Roman R, Pienik R, Futerman AH. Increased ceramide synthase 2 and 6 mRNA levels in breast cancer tissues and correlation with sphingosine kinase expression. Biochem Biophys Res Commun (2010) 391:219–23. doi:10.1016/j.bbrc.2009.11.035

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Accioly MT, Pacheco P, Maya-Monteiro CM, Carrossini N, Robbs BK, Oliveira SS, et al. Lipid bodies are reservoirs of cyclooxygenase-2 and sites of prostaglandin-E2 synthesis in colon cancer cells. Cancer Res (2008) 68:1732–40. doi:10.1158/0008-5472.CAN-07-1999

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Fernandis AZ, Wenk MR. Lipid-based biomarkers for cancer. J Chromatogr B Analyt Technol Biomed Life Sci (2009) 877:2830–5. doi:10.1016/j.jchromb.2009.06.015

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Liu TG, Huang Y, Cui DD, Huang XB, Mao SH, Ji LL, et al. Inhibitory effect of ginsenoside Rg3 combined with gemcitabine on angiogenesis and growth of lung cancer in mice. BMC Cancer (2009) 9:250. doi:10.1186/1471-2407-9-250

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Choi E, Carruthers K, Zhang L, Thomas N, Battaglino RA, Morse LR, et al. Concurrent muscle and bone deterioration in a murine model of cancer cachexia. Physiol Rep (2013) 1:e00144. doi:10.1002/phy2.144

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Reeves PG, Nielsen FH, Fahey GCJ. AIN-93 purified diets for laboratory rodents: final report of the American Institute of Nutrition Ad Hoc writing committee on the reformulation of the AIN-76A rodent diet. J Nutr (1993) 123:1939–51. doi:10.1093/jn/123.11.1939

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Brodt P. Characterization of two highly metastatic variants of Lewis lung carcinoma with different organ specificities. Cancer Res (1986) 46:2442–8.

PubMed Abstract | Google Scholar

29. Sundaram S, Yan L. Time-restricted feeding reduces adiposity in mice fed a high-fat diet. Nutr Res (2016) 36:603–11. doi:10.1016/j.nutres.2016.02.005

CrossRef Full Text | Google Scholar

30. Zacek P, Bukowski M, Rosenberger TA, Picklo M. Quantitation of isobaric phosphatidylcholine species in human plasma using a hybrid quadrupole linear ion-trap mass spectrometer. J Lipid Res (2016) 57:2225–34. doi:10.1194/jlr.D070656

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Zacek P, Bukowski M, Johnson LK, Raatz SK, Picklo M. Selective enrichment of n-3 fatty acids in human plasma lipid motifs following intake of marine fish. J Nutr Biochem (2017) 54:57–65. doi:10.1016/j.jnutbio.2017.11.002

CrossRef Full Text | Google Scholar

32. Surma MA, Herzog R, Vasilj A, Klose C, Christinat N, Morin-Rivron D, et al. An automated shotgun lipidomics platform for high throughput, comprehensive, and quantitative analysis of blood plasma intact lipids. Eur J Lipid Sci Technol (2015) 117:1540–9. doi:10.1002/ejlt.201500145

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Liebisch G, Vizcaino JA, Kofeler H, Trotzmuller M, Griffiths WJ, Schmitz G, et al. Shorthand notation for lipid structures derived from mass spectrometry. J Lipid Res (2013) 54:1523–30. doi:10.1194/jlr.M033506

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Liebisch G, Drobnik W, Reil M, Trumbach B, Arnecke R, Olgemoller B, et al. Quantitative measurement of different ceramide species from crude cellular extracts by electrospray ionization tandem mass spectrometry (ESI-MS/MS). J Lipid Res (1999) 40:1539–46.

PubMed Abstract | Google Scholar

35. Han X, Gross RW. Quantitative analysis and molecular species fingerprinting of triacylglyceride molecular species directly from lipid extracts of biological samples by electrospray ionization tandem mass spectrometry. Anal Biochem (2001) 295:88–100. doi:10.1006/abio.2001.5178

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Picklo MJ, Thyfault JP. Vitamin E and vitamin C do not reduce insulin sensitivity but inhibit mitochondrial protein expression in exercising obese rats. Appl Physiol Nutr Metab (2015) 40:343–52. doi:10.1139/apnm-2014-0302

CrossRef Full Text | Google Scholar

37. Harris VM. Protein detection by simple western analysis. Methods Mol Biol (2015) 1312:465–8. doi:10.1007/978-1-4939-2694-7_47

CrossRef Full Text | Google Scholar

38. van den Berg RA, Hoefsloot HCJ, Westerhuis JA, Smilde AK, Van Der Werf MJ. Centering, scaling, and transformations: improving the biological information content of metabolomics data. BMC Genomics (2006) 7:142. doi:10.1186/1471-2164-7-142

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Xia J, Wishart DS. Using metaboanalyst 3.0 for comprehensive metabolomics data analysis. Curr Protoc Bioinformatics (2016) 55:14.10.1–91. doi:10.1002/cpbi.11

CrossRef Full Text | Google Scholar

40. Zaidi N, Lupien L, Kuemmerle NB, Kinlaw WB, Swinnen JV, Smans K. Lipogenesis and lipolysis: the pathways exploited by the cancer cells to acquire fatty acids. Prog Lipid Res (2013) 52:585–9. doi:10.1016/j.plipres.2013.08.005

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Yan L, Demars LC. Effects of dietary fat on spontaneous metastasis of Lewis lung carcinoma in mice. Clin Exp Metastasis (2010) 27:581–90. doi:10.1007/s10585-010-9347-7

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Singer S, Souza K, Thilly WG. Pyruvate utilization, phosphocholine and adenosine triphosphate (ATP) are markers of human breast tumor progression: a 31P- and 13C-nuclear magnetic resonance (NMR) spectroscopy study. Cancer Res (1995) 55:5140–5.

PubMed Abstract | Google Scholar

43. Sterin M, Cohen JS, Ringel I. Hormone sensitivity is reflected in the phospholipid profiles of breast cancer cell lines. Breast Cancer Res Treat (2004) 87:1–11. doi:10.1023/B:BREA.0000041572.07837.ec

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Zheng M, Fang H, Tsuruoka T, Tsuji T, Sasaki T, Hakomori S. Regulatory role of GM3 ganglioside in alpha 5 beta 1 integrin receptor for fibronectin-mediated adhesion of FUA169 cells. J Biol Chem (1993) 268:2217–22.

PubMed Abstract | Google Scholar

45. Kazarian T, Jabbar AA, Wen FQ, Patel DA, Valentino LA. Gangliosides regulate tumor cell adhesion to collagen. Clin Exp Metastasis (2003) 20:311–9. doi:10.1023/A:1024021707604

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Chavez JA, Summers SA. Characterizing the effects of saturated fatty acids on insulin signaling and ceramide and diacylglycerol accumulation in 3T3-L1 adipocytes and C2C12 myotubes. Arch Biochem Biophys (2003) 419:101–9. doi:10.1016/j.abb.2003.08.020

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Holland WL, Bikman BT, Wang LP, Yuguang G, Sargent KM, Bulchand S, et al. Lipid-induced insulin resistance mediated by the proinflammatory receptor TLR4 requires saturated fatty acid-induced ceramide biosynthesis in mice. J Clin Invest (2011) 121:1858–70. doi:10.1172/JCI43378

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Duarte JA, Carvalho F, Pearson M, Horton JD, Browning JD, Jones JG, et al. A high-fat diet suppresses de novo lipogenesis and desaturation but not elongation and triglyceride synthesis in mice. J Lipid Res (2014) 55:2541–53. doi:10.1194/jlr.M052308

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Nath A, Chan C. Genetic alterations in fatty acid transport and metabolism genes are associated with metastatic progression and poor prognosis of human cancers. Sci Rep (2016) 6:18669. doi:10.1038/srep18669

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Yan L, DeMars LC. Effects of a high-fat diet on spontaneous metastasis of Lewis lung carcinoma in plasminogen activator inhibitor-1 deficient and wild-type mice. PLoS One (2014) 9:e110869. doi:10.1371/journal.pone.0110869

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Institute for Laboratory Animal Research. Guide for the Care and Use of Laboratory Animals. Washington, DC: National Academies Press (2011).

Google Scholar

Keywords: lipidomics, Lewis lung carcinoma, primary tumor, high-fat diet, mice

Citation: Sundaram S, Žáček P, Bukowski MR, Mehus AA, Yan L and Picklo MJ (2018) Lipidomic Impacts of an Obesogenic Diet Upon Lewis Lung Carcinoma in Mice. Front. Oncol. 8:134. doi: 10.3389/fonc.2018.00134

Received: 28 February 2018; Accepted: 13 April 2018;
Published: 11 May 2018

Edited by:

Saverio Marchi, University of Ferrara, Italy

Reviewed by:

Eugenia Morselli, Pontificia Universidad Católica de Chile, Chile
Alessia Ligresti, National Research Council of Italy – ICB, Italy

Copyright: © 2018 Sundaram, Žáček, Bukowski, Mehus, Yan and Picklo. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Matthew J. Picklo, matthew.picklo@ars.usda.gov

Present address: Petr Žáček, Charles University, BIOCEV, Prague, Czechia