Fatty acid desaturase 1 (FADS1) is a cancer marker for patient survival and a potential novel target for precision cancer treatment

Disruption of polyunsaturated fatty acid biosynthesis drives STING-dependent acute myeloid leukemia cell maturation and death

The role of polyunsaturated fatty acid (PUFA) biosynthesis in acute myeloid leukemia (AML) remains largely undefined. A comparative expression analysis of 35 genes encoding fatty acid biosynthesis enzymes showed that fatty acid desaturase 1 (FADS1) was highly expressed across multiple AML subtypes relative to healthy controls and that elevated FADS1 expression correlates with worse overall AML patient survival. Functionally, shRNA-mediated inhibition of FADS1 reduced AML cell growth in vitro and significantly delayed leukemia onset in an AML mouse model. AML cell lines depleted of FADS1 arrested in the G1/S-phase of the cell cycle, acquired characteristics of myeloid maturation and subsequently died. To understand the molecular consequences of FADS1 inhibition, a combination of mass spectrometry-based analysis of complex lipids and gene expression analysis (RNA-seq) was performed. FADS1 inhibition caused AML cells to exhibit significant lipidomic remodeling, including depletion of PUFAs from the phospholipids, phosphatidylserine, and phosphatidylethanolamine. These lipidomic alterations were accompanied by an increase induction of inflammatory and stimulator of interferon genes (STING)-mediated type-1 interferon signaling. Remarkably, genetic deletion of STING largely prevented the AML cell maturation and death phenotypes mediated by FADS1 inhibition. Highlighting the therapeutic implications of these findings, pharmacological blockade of PUFA biosynthesis reduced patient-derived AML cell numbers ex vivo but not that of healthy donor cells. Similarly, STING agonism attenuated patient-derived-AML survival; however, STING activation also reduced healthy granulocyte numbers. Collectively, these data unveil a previously unrecognized importance of PUFA biosynthesis in leukemogenesis and that imbalances in PUFA metabolism can drive STING-mediated AML maturation and death.

Targeting Fatty Acid Desaturase I Inhibits Renal Cancer Growth Via ATF3-mediated ER Stress Response

Monounsaturated fatty acids (MUFAs) play a pivotal role in maintaining endoplasmic reticulum (ER) homeostasis, an emerging hallmark of cancer. However, the role of polyunsaturated fatty acid (PUFAs) desaturation in persistent ER stress driven by oncogenic abnormalities remains elusive. Fatty Acid Desaturase 1 (FADS1) is a rate-limiting enzyme controlling the bioproduction of long-chain PUFAs. Our previous research has demonstrated the significant role of FADS1 in cancer survival, especially in kidney cancers. We explored the underlying mechanism in this study. We found that pharmacological inhibition or knockdown of the expression of FADS1 effectively inhibits renal cancer cell proliferation and induces cell cycle arrest. The stable knockdown of FADS1 also significantly inhibits tumor formation in vivo. Mechanistically, we show that while FADS1 inhibition induces ER stress, its expression is also augmented by ER-stress inducers. Notably, FADS1-inhibition sensitized cellular response to ER stress inducers, providing evidence of FADS1's role in modulating the ER stress response in cancer cells. We show that, while FADS1 inhibition-induced ER stress leads to activation of ATF3, ATF3-knockdown rescues the FADS1 inhibition-induced ER stress and cell growth suppression. In addition, FADS1 inhibition results in the impaired biosynthesis of nucleotides and decreases the level of UPD-N-Acetylglucosamine, a critical mediator of the unfolded protein response. Our findings suggest that PUFA desaturation is crucial for rescuing cancer cells from persistent ER stress, supporting FADS1 as a new therapeutic target.

PancanQTLv2.0: a comprehensive resource for expression quantitative trait loci across human cancers

Expression quantitative trait locus (eQTL) analysis is a powerful tool used to investigate genetic variations in complex diseases, including cancer. We previously developed a comprehensive database, PancanQTL, to characterize cancer eQTLs using The Cancer Genome Atlas (TCGA) dataset, and linked eQTLs with patient survival and GWAS risk variants. Here, we present an updated version, PancanQTLv2.0 (https://hanlaboratory.com/PancanQTLv2/), with advancements in fine-mapping causal variants for eQTLs, updating eQTLs overlapping with GWAS linkage disequilibrium regions and identifying eQTLs associated with drug response and immune infiltration. Through fine-mapping analysis, we identified 58 747 fine-mapped eQTLs credible sets, providing mechanic insights of gene regulation in cancer. We further integrated the latest GWAS Catalog and identified a total of 84 592 135 linkage associations between eQTLs and the existing GWAS loci, which represents a remarkable ∼50-fold increase compared to the previous version. Additionally, PancanQTLv2.0 uncovered 659516 associations between eQTLs and drug response and identified 146948 associations between eQTLs and immune cell abundance, providing potentially clinical utility of eQTLs in cancer therapy. PancanQTLv2.0 expanded the resources available for investigating gene expression regulation in human cancers, leading to advancements in cancer research and precision oncology.

FADS2 function at the major cancer hotspot 11q13 locus alters fatty acid metabolism in cancer

Dysregulation of fatty acid metabolism and de novo lipogenesis is a key driver of several cancer types through highly unsaturated fatty acid (HUFA) signaling precursors such as arachidonic acid. The human chromosome 11q13 locus has long been established as the most frequently amplified in a variety of human cancers. The fatty acid desaturase genes (FADS1, FADS2 and FADS3) responsible for HUFA biosynthesis localize to the 11q12-13.1 region. FADS2 activity is promiscuous, catalyzing biosynthesis of several unsaturated fatty acids by Δ6, Δ8, and Δ4 desaturation. Our main aim here is to review known and putative consequences of FADS2 dysregulation due to effects on the 11q13 locus potentially driving various cancer types. FADS2 silencing causes synthesis of sciadonic acid (5Z,11Z,14Z-20:3) in MCF7 cells and breast cancer in vivo. 5Z,11Z,14Z-20:3 is structurally identical to arachidonic acid (5Z,8Z,11Z,14Z-20:4) except it lacks the internal Δ8 double bond required for prostaglandin and leukotriene synthesis, among other eicosanoids. Palmitic acid has substrate specificity for both SCD and FADS2. Melanoma, prostate, liver and lung cancer cells insensitive to SCD inhibition show increased FADS2 activity and sapienic acid biosynthesis. Elevated serum mead acid levels found in hepatocellular carcinoma patients suggest an unsatisfied demand for arachidonic acid. FADS2 circular RNAs are at high levels in colorectal and lung cancer tissues. FADS2 circular RNAs are associated with shorter overall survival in colorectal cancer patients. The evidence thusfar supports an effort for future research on the role of FADS2 as a tumor suppressor in a range of neoplastic disorders.

Melanoma cells induce dedifferentiation and metabolic changes in adipocytes present in the tumor niche

BACKGROUND

One of the factors that affect the progression of melanoma is the tumor microenvironment, which consists of cellular elements, extracellular matrix, acidification, and a hypoxic state. Adipocytes are one of the types of cell present in the niche and are localized in the deepest layer of the skin. However, the relationship between fat cells and melanoma remains unclear.

METHODS

We assessed the influence of melanoma cells on adipocytes using an indirect coculture system. We estimated the level of cancer-associated adipocyte (CAA) markers through quantitative PCR analysis. The fibroblastic phenotype of CAAs was confirmed by cell staining and western blotting analysis. The lipid content was estimated by lipid detection in CAAs using LipidSpot and by quantitative analysis using Oil Red O. The expression of proteins involved in lipid synthesis, delipidation, and metabolic processes were assessed through quantitative PCR or western blotting analysis. Lactate secretion was established using a Lactate-Glo™ assay. Proteins secreted by CAAs were identified in cytokine and angiogenesis arrays. The proliferation of melanoma cells cocultured with CAAs was assessed using an XTT proliferation assay. Statistical analysis was performed using a one-way ANOVA followed by Tukey's test in GraphPad Prism 7 software.

RESULTS

Obtained CAAs were identified by decreased levels of leptin, adiponectin, resistin, and FABP4. Adipocytes cocultured with melanoma presented fibroblastic features, such as a similar proteolytic pattern to that of 3T3L1 fibroblasts and increased levels of vimentin and TGFβRIII. Melanoma cells led to a reduction of lipid content in CAAs, possibly by downregulation of lipid synthesis pathways (lower FADS, SC4MOL, FASN) or enhancement of lipolysis (higher level of phosphorylation of ERK and STAT3). Adipocytes cocultured with melanoma cells secreted higher IL6 and SerpinE1 levels and produced less CCL2, CXCL1, and angiogenic molecules. CAAs also showed metabolic changes comprising the increased secretion of lactate and enhanced production of glucose, lactate, and ion transporters. In addition, changes in adipocytes observed following melanoma coculture resulted in a higher proliferation rate of cancer cells.

CONCLUSIONS

Melanoma cells led to decreased lipid content in adipocytes, which might be related to enhanced delipidation or reduction of lipid synthesis. Fibroblast-like CAAs showed metabolic changes that may be the reason for accelerated proliferation of melanoma cells.

Delineation of fatty acid metabolism in gastric cancer: Therapeutic implications

BACKGROUND

The prognosis of gastric cancer is extremely poor. Metabolic reprogramming involving lipids has been associated with cancer occurrence and progression.

AIM

To illustrate fatty acid metabolic mechanisms in gastric cancer, detect core genes, develop a prognostic model, and provide treatment options.

METHODS

Raw data from The Cancer Genome Atlas and Gene Expression Omnibus databases were collected and analyzed. Differentially expressed fatty acid metabolism genes were identified and incorporated into a risk model based on least absolute shrinkage and selection operator regression analysis. Then, patients from The Cancer Genome Atlas were assigned to high- and low-risk cohorts according to the mean value of the risk score as the threshold, which was verified in the Gene Expression Omnibus database. Relationships between chemotherapeutic sensitivity and tumor microenvironment features were assessed.

RESULTS

An integrated evaluation was performed in this study. Fatty acid metabolism-related genes were used to construct the risk model. Patients classified into the high-risk cohort were considered to be resistant to chemotherapy based on results of the "pRRophetic" R package. Patients in the high-risk cohort were associated with type I/II interferon activation, increased inflammation level, immune cell infiltration, and tumor immune dysfunction based on the exclusion algorithm, indicating the potential benefit of immunotherapy in these patients.

CONCLUSION

We constructed a fatty acid-related risk score model to assess the comprehensive fatty acid features in gastric cancer and validated its vital role in prognosis, chemotherapy sensitivity, and immunotherapy.

Delineation of fatty acid metabolism in gastric cancer: Therapeutic implications

BACKGROUND

The prognosis of gastric cancer is extremely poor. Metabolic reprogramming involving lipids has been associated with cancer occurrence and progression.

AIM

To illustrate fatty acid metabolic mechanisms in gastric cancer, detect core genes, develop a prognostic model, and provide treatment options.

METHODS

Raw data from The Cancer Genome Atlas and Gene Expression Omnibus databases were collected and analyzed. Differentially expressed fatty acid metabolism genes were identified and incorporated into a risk model based on least absolute shrinkage and selection operator regression analysis. Then, patients from The Cancer Genome Atlas were assigned to high- and low-risk cohorts according to the mean value of the risk score as the threshold, which was verified in the Gene Expression Omnibus database. Relationships between chemotherapeutic sensitivity and tumor microenvironment features were assessed.

RESULTS

An integrated evaluation was performed in this study. Fatty acid metabolism-related genes were used to construct the risk model. Patients classified into the high-risk cohort were considered to be resistant to chemotherapy based on results of the “pRRophetic” R package. Patients in the high-risk cohort were associated with type I/II interferon activation, increased inflammation level, immune cell infiltration, and tumor immune dysfunction based on the exclusion algorithm, indicating the potential benefit of immunotherapy in these patients.

CONCLUSION

We constructed a fatty acid-related risk score model to assess the comprehensive fatty acid features in gastric cancer and validated its vital role in prognosis, chemotherapy sensitivity, and immunotherapy.

Assessment of Colorectal Cancer Risk Factors through the Application of Network-Based Approaches in a Racially Diverse Cohort of Colon Organoid Stem Cells

Numerous demographic factors have been associated with colorectal cancer (CRC) risk. To better define biological mechanisms underlying these associations, we performed RNA sequencing of stem-cell-enriched organoids derived from the healthy colons of seven European Americans and eight African Americans. A weighted gene co-expression network analysis was performed following RNA sequencing. Module-trait relationships were determined through the association testing of each module and five CRC risk factors (age, body mass index, sex, smoking history, and race). Only modules that displayed a significantly positive correlation for gene significance and module membership were considered for further investigation. In total, 16 modules were associated with known CRC risk factors (p < 0.05). To contextualize the role of risk modules in CRC, publicly available RNA-sequencing data from TCGA-COAD were downloaded and re-analyzed. Differentially expressed genes identified between tumors and matched normal-adjacent tissue were overlaid across each module. Loci derived from CRC genome-wide association studies were additionally overlaid across modules to identify robust putative targets of risk. Among them, MYBL2 and RXRA represented strong plausible drivers through which cigarette smoking and BMI potentially modulated CRC risk, respectively. In summary, our findings highlight the potential of the colon organoid system in identifying novel CRC risk mechanisms in an ancestrally diverse and cellularly relevant population.

LncRNA LINC01569 promotes M2 macrophage polarization to accelerate hypopharyngeal carcinoma progression through the miR-193a-5p/FADS1 signaling axis

Background: Long non-coding RNA (lncRNA) LINC01569 plays an important role in regulating the tumor microenvironment (TME) and macrophage polarization. However, whether it participates in the progression of hypopharyngeal carcinoma by regulating the TME remains unclear. Methods: An online database was used to analyze clinical data. Macrophage polarization was detected using qRT-PCR and flow cytometry. In vivo experiments were performed using tumor-bearing nude mice. A co-culture system of hypopharyngeal carcinoma cells and macrophages was used to explore the interactions between the two cell types. Results: LINC01569 enhancement was observed in hypopharyngeal carcinoma tumor-associated macrophages (TAMs). In IL4-induced M2 macrophages, the expression of LINC01569 increased, while LINC01569 expression declined significantly in LPS-induced M1 macrophages. SiRNA-mediated downregulation of LINC01569 inhibits IL4-induced M2 macrophage polarization. Using online databases and a dual-luciferase reporter, miR-193a-5p was confirmed as a potential downstream sponge of LINC01569. MiR-193a-5p expression decreased in IL4-mediated M2 macrophages, which was restored by LINC01569 downregulation. Additionally, LINC01569 inhibition-mediated blocking of M2 macrophage polarization was moderately abolished by transfection with the miR-193a-5p inhibitor. Fatty acid desaturase 1 (FADS1) was verified as a downstream target of miR-193a-5p, and LINC01569 downregulation-mediated inhibition of FADS1 was blocked by miR-193a-5p mimics. Importantly, LINC01569 downregulation-mediated decline in M2 macrophage polarization was abolished by miR-193a-5p mimics, which was further reversed by FADS1 knockdown. Implantation of a mixture of FaDu cells and IL4-induced macrophages promoted tumor growth and proliferation, which were abrogated by the knockdown of LINC01569 in macrophages. Using an in co-culture system of FaDu cells and macrophages in vitro, M2 macrophage-regulated cell growth and apoptosis of FaDu cells were found to be mediated by the LINC01569/miR-193a-5p signaling axis. Conclusion: LINC01569 is highly expressed in the TAMs of hypopharyngeal carcinoma. LINC01569 downregulation restrains macrophages from polarizing toward M2 through the miR-193a-5p/FADS1 signaling axis, thereby helping tumor cells escape inherent immune surveillance and promoting the occurrence and development of hypopharyngeal carcinoma.

Postoperative Metabolic Phenoreversion in Clear Cell Renal Cell Carcinoma

The ultimate goal of surgical treatment in cancer is to remove the tumor mass for restoring a healthy state. A 16-lipid panel that discriminated healthy controls from clear cell renal cell carcinoma (ccRCC) patients in a prior study was evaluated in the present work in paired-serum samples collected from patients (n = 41) before and after nephrectomy. Changes in the lipid and metabolite fingerprints from ccRCC patients were investigated and compared with fingerprints from healthy individuals obtained by means of ultra-performance liquid chromatography-high-resolution mass spectrometry. The lipid panel differentiated phenotypes associated with metabolic restoration after surgery, representing a serum signature of phenoreversion to a healthy metabolic state. In particular, PC 16:0/0:0, PC 18:2/18:2, and linoleic acid allowed discriminating serum samples from ccRCC patients with poor prognosis from those with an improved outcome during the follow-up period. Ratios of PC 16:0/0:0 and PC 18:2/18:2 with linoleic acid levels may contribute as prognostic tools to support decision-making during the patient follow-up care. The preliminary character of these results should be validated with larger cohorts, including subjects with different ethnicities, life style, and diets. MetaboLights study references: MTBLS1839, MTBLS3838, and MTBLS4629.

Investigating the role of FADS family members in breast cancer based on bioinformatic analysis and experimental validation

Breast cancer (BC) is the most common malignant tumor in women worldwide. Emerging evidence indicates the significance of fatty acid metabolism in BC. Fatty acid desaturase (FADS) is closely associated with cancer occurrence and development. Here, bioinformatic analysis and experimental validation were applied to investigate the potential functions of FADS in BC. Several public databases, including TCGA, GEO, HPA, Kaplan-Meier plotter, STRING, DAVID, cBioPortal, TIMER, TRRUST, and LinkedOmics were used to determine mRNA/protein expression levels, prognostic significance, functional enrichment, genetic alterations, association with tumor-infiltrating immune cells, and related transcription factors and kinases. BC tissues showed higher and lower mRNA expression of FADS2/6/8 and FADS3/4/5, respectively. FADS1/2/6 and FADS3/4/5 showed higher and lower protein expression levels, respectively, in BC tissues. Moreover, FADS1/7 up- and FADS3/8 down-regulation predicted poor overall and recurrence-free survival, while FADS2/5 up- and FADS4 down-regulation were associated with poor recurrence-free survival. Receiver operating characteristic curves revealed that FADS2/3/4/8 were indicative diagnostic markers. FADS family members showing differential expression levels were associated with various clinical subtypes, clinical stages, lymph node metastasis status, copy number variants, DNA methylation, and miRNA regulation in BC. The mRNA expression level of FADS1/2/3/4/5/7/8 was observed to be significantly negatively correlated with DNA methylation. FADS1/2 upregulation was significantly correlated with clinical stages. FADS1/4 expression was obviously lower in BC patients with higher lymph node metastasis than lower lymph node metastasis, while FADS7/8 expression was obviously higher in BC patients with higher lymph node metastasis than lower lymph node metastasis. FADS family members showed varying degrees of genetic alterations, and Gene Ontology and KEGG pathway enrichment analyses suggested their involvement in lipid metabolism. Their expression level was correlated with immune cell infiltration levels. FADS2 was chosen for further validation analyses. We found FADS2 to be significantly over-expressed in clinical BC tissue samples. The proliferation, migration, and invasion abilities of MDA-MB-231 and BT474 cells were significantly reduced after FADS2 knockdown. Furthermore, FADS2 may promote the occurrence and development of BC cells via regulating the epithelial-mesenchymal transition (EMT) pathway. Altogether, our results suggest that FADS1/2/3/4 can serve as potential therapeutic targets, prognostic indicators, and diagnostic markers in patients with BC.