Somatic mutations in arachidonic acid metabolism pathway genes enhance oral cancer post-treatment disease-free survival

Selenomethionine in gelatin methacryloyl hydrogels: Modulating ferroptosis to attenuate skin aging

During skin aging, the degeneration of epidermal stem cells (EpiSCs) leads to diminished wound healing capabilities and epidermal disintegration. This study tackles this issue through a comprehensive analysis combining transcriptomics and untargeted metabolomics, revealing age-dependent alterations in the Gpx gene family and arachidonic acid (AA) metabolic networks, resulting in enhanced ferroptosis. Selenomethionine (Se-Met) could enhance GPX4 expression, thereby assisting EpiSCs in countering AA-induced mitochondrial damage and ferroptosis. Additionally, Se-Met demonstrates antioxidative characteristics and extensive ultraviolet absorption. For the sustained and controllable release of Se-Met, it was covalently grafted to UV-responsive GelMA hydrogels via AC-PEG-NHS tethers. The Se-Met@GelMA hydrogel effectively accelerated wound healing in a chronological aging mice model, by inhibiting lipid peroxidation and ferroptosis with augmented GPX4 expression. Moreover, in a photoaging model, this hydrogel significantly mitigated inflammatory responses, extracellular matrix remodeling, and ferroptosis in UV-exposed mice. These characteristics render Se-Met@GelMA hydrogel valuable in practical clinical applications.

Pathological and pharmacological functions of the metabolites of polyunsaturated fatty acids mediated by cyclooxygenases, lipoxygenases, and cytochrome P450s in cancers

Oxylipins have garnered increasing attention because they were consistently shown to play pathological and/or pharmacological roles in the development of multiple cancers. Oxylipins are the metabolites of polyunsaturated fatty acids via both enzymatic and nonenzymatic pathways. The enzymes mediating the metabolism of PUFAs include but not limited to lipoxygenases (LOXs), cyclooxygenases (COXs), and cytochrome P450s (CYPs) pathways, as well as the down-stream enzymes. Here, we systematically summarized the pleiotropic effects of oxylipins in different cancers through pathological and pharmacological aspects, with specific reference to the enzyme-mediated oxylipins. We discussed the specific roles of oxylipins on cancer onset, growth, invasion, and metastasis, as well as the expression changes in the associated metabolic enzymes and the associated underlying mechanisms. In addition, we also discussed the clinical application and potential of oxylipins and related metabolic enzymes as the targets for cancer prevention and treatment. We found the specific function of most oxylipins in cancers, especially the underlying mechanisms and clinic applications, deserves and needs further investigation. We believe that research on oxylipins will provide not only more therapeutic targets for various cancers but also dietary guidance for both cancer patients and healthy humans.

Multiomics characterization of fatty acid metabolism for the clinical management of hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is a prevalent malignancy and there is a lack of effective biomarkers for HCC diagnosis. Living organisms are complex, and different omics molecules interact with each other to implement various biological functions. Genomics and metabolomics, which are the top and bottom of systems biology, play an important role in HCC clinical management. Fatty acid metabolism is associated with malignancy, prognosis, and immune phenotype in cancer, which is a potential hallmark in malignant tumors. In this study, the genes and metabolites related to fatty acid metabolism were thoroughly investigated by a dynamic network construction algorithm named EWS-DDA for the early diagnosis and prognosis of HCC. Three gene ratios and eight metabolite ratios were identified by EWS-DDA as potential biomarkers for HCC clinical management. Further analysis using biological analysis, statistical analysis and document validation in the discovery and validation sets suggested that the selected potential biomarkers had great clinical prognostic value and helped to achieve effective early diagnosis of HCC. Experimental results suggested that in-depth evaluation of fatty acid metabolism from different omics viewpoints can facilitate the further understanding of pathological alterations associated with HCC characteristics, improving the performance of early diagnosis and clinical prognosis.

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.

Progression of prostate cancer reprograms MYC-mediated lipid metabolism via lysine methyltransferase 2A

BACKGROUND

The activities of MYC, the androgen receptor, and its associated pioneer factors demonstrate substantial reprogramming between early and advanced prostate cancer. Although previous studies have shown a shift in cellular metabolic requirements associated with prostate cancer progression, the epigenetic regulation of these processes is incompletely described. Here, we have integrated chromatin immunoprecipitation sequencing (ChIP-seq) and whole-transcriptome sequencing to identify novel regulators of metabolism in advanced prostate tumors characterized by elevated MYC activity.

RESULTS

Using ChIP-seq against MYC, HOXB13, and AR in LNCaP cells, we observed redistribution of co-bound sites suggestive of differential KMT2A activity as a function of MYC expression. In a cohort of 177 laser-capture microdissected foci of prostate tumors, KMT2A expression was positively correlated with MYC activity, AR activity, and HOXB13 expression, but decreased with tumor grade severity. However, KMT2A expression was negatively correlated with these factors in 25 LuCaP patient-derived xenograft models of advanced prostate cancer and 99 laser-capture microdissected foci of metastatic castration-resistant prostate cancer. Stratified by KMT2A expression, ChIP-seq against AR and HOXB13 in 15 LuCaP patient-derived xenografts showed an inverse association with sites involving genes implicated in lipid metabolism, including the arachidonic acid metabolic enzyme PLA2G4F. LuCaP patient-derived xenograft models grown as organoids recapitulated the inverse association between KMT2A expression and fluorine-18 labeled arachidonic acid uptake in vitro.

CONCLUSIONS

Our study demonstrates that the epigenetic activity of transcription factor oncogenes exhibits a shift during prostate cancer progression with distinctive phenotypic effects on metabolism. These epigenetically driven changes in lipid metabolism may serve as novel targets for the development of novel imaging agents and therapeutics.

Aspirin Inhibition of Group VI Phospholipase A2 Induces Synthetic Lethality in AAM Pathway Down-Regulated Gingivobuccal Squamous Carcinoma

Background: To elucidate the role of iPLA2/PLA2G6 in gingivobuccal squamous cell carcinoma (GB-SCC) and to ascertain the synthetic lethality-based chemoprevention role of aspirin in arachidonic acid metabolism (AAM) pathway down-regulated GB-SCC. Methods: The in vitro efficacy of aspirin on GB-SCC cells (ITOC-03 and ITOC-04) was assessed by cell proliferation, colony formation, apoptosis, cell migration, cell cycle assay and RNA-seq, while inhibition of PLA2G6 and AAM pathway components was affirmed by qPCR, Western blot and immunofluorescence staining. The in vivo effect of aspirin was evaluated using NOD-SCID mice xenografts and immunohistochemical analysis. Results: We found that aspirin, which has been reported to act through the COX pathway, is inhibiting PLA2G6, and thereby the COX and LOX components of the AAM pathway. The findings were validated using PLA2G6 siRNA and immunohistochemical marker panel. Moreover, a pronounced effect in ITOC-04 cells and xenografts implied aspirin-induced synthetic lethality in the AAM pathway down-regulated GB-SCC. Conclusions: This study reveals that aspirin induces the anti-tumor effect by a previously unrecognized mechanism of PLA2G6 inhibition. In addition, the effect of aspirin is influenced by the baseline AAM pathway status and could guide precision prevention clinical trials of AAM pathway inhibitors.

Inhibition of PLA2G4E/cPLA2 promotes survival of random skin flaps by alleviating Lysosomal membrane permeabilization-Induced necroptosis

Necrosis that appears at the ischemic distal end of random-pattern skin flaps increases the pain and economic burden of patients. Necroptosis is thought to contribute to flap necrosis. Lysosomal membrane permeabilization (LMP) plays an indispensable role in the regulation of necroptosis. Nonetheless, the mechanisms by which lysosomal membranes become leaky and the relationship between necroptosis and lysosomes are still unclear in ischemic flaps. Based on Western blotting, immunofluorescence, enzyme-linked immunosorbent assay, and liquid chromatography-mass spectrometry (LC-MS) analysis results, we found that LMP was presented in the ischemic distal portion of random-pattern skin flaps, which leads to disruption of lysosomal function and macroautophagic/autophagic flux, increased necroptosis, and aggravated necrosis of the ischemic flaps. Moreover, bioinformatics analysis of the LC-MS results enabled us to focus on the role of PLA2G4E/cPLA2 (phospholipase A2, group IVE) in LMP of the ischemic flaps. In vivo inhibition of PLA2G4E with an adeno-associated virus vector attenuated LMP and necroptosis, and promoted flap survival. In addition, microRNA-seq helped us determine that Mir504-5p was differentially expressed in ischemic flaps. A string of in vitro and in vivo tests was employed to verify the inhibitory effect of Mir504-5p on PLA2G4E, LMP and necroptosis. Finally, we concluded that the inhibition of PLA2G4E by Mir504-5p reduced LMP-induced necroptosis, thereby promoting the survival of random-pattern skin flaps.

The Impact of m1A Methylation Modification Patterns on Tumor Immune Microenvironment and Prognosis in Oral Squamous Cell Carcinoma

N1-methyladenosine (m1A) modification widely participates in the occurrence and progression of numerous diseases. Nevertheless, the potential roles of m1A in the tumor immune microenvironment (TIME) are still not fully understood. Based on 10 m1A methylation regulators, we comprehensively explored the m1A modification patterns in 502 patients with oral squamous cell carcinoma (OSCC). The m1A modification patterns were correlated with TIME characteristics and the m1A score was established to evaluate the effect of the m1A modification patterns on individual OSCC patients. The TIME characteristics and survival outcomes under the three m1A modification patterns were significantly distinct. OSCC patients in the high m1A score group were characterized by poorer prognosis, lower immune infiltration, lower ssGSEA score, lower expression levels of immune checkpoint molecules, and higher tumor mutation loads. The present study revealed that m1A modification might be associated with the TIME in OSCC, and has potential predictive ability for the prognosis of OSCC.

Molecular Subtypes of Oral Squamous Cell Carcinoma Based on Immunosuppression Genes Using a Deep Learning Approach

Background: The mechanisms through which immunosuppressed patients bear increased risk and worse survival in oral squamous cell carcinoma (OSCC) are unclear. Here, we used deep learning to investigate the genetic mechanisms underlying immunosuppression in the survival of OSCC patients, especially from the aspect of various survival-related subtypes. Materials and methods: OSCC samples data were obtained from The Cancer Genome Atlas (TCGA), International Cancer Genome Consortium (ICGC), and OSCC-related genetic datasets with survival data in the National Center for Biotechnology Information (NCBI). Immunosuppression genes (ISGs) were obtained from the HisgAtlas and DisGeNET databases. Survival analyses were performed to identify the ISGs with significant prognostic values in OSCC. A deep learning (DL)-based model was established for robustly differentiating the survival subpopulations of OSCC samples. In order to understand the characteristics of the different survival-risk subtypes of OSCC samples, differential expression analysis and functional enrichment analysis were performed. Results: A total of 317 OSCC samples were divided into one inferring cohort (TCGA) and four confirmation cohorts (ICGC set, GSE41613, GSE42743, and GSE75538). Eleven ISGs (i.e., BGLAP, CALCA, CTLA4, CXCL8, FGFR3, HPRT1, IL22, ORMDL3, TLR3, SPHK1, and INHBB) showed prognostic value in OSCC. The DL-based model provided two optimal subgroups of TCGA-OSCC samples with significant differences (p = 4.91E-22) and good model fitness [concordance index (C-index) = 0.77]. The DL model was validated by using four external confirmation cohorts: ICGC cohort (n = 40, C-index = 0.39), GSE41613 dataset (n = 97, C-index = 0.86), GSE42743 dataset (n = 71, C-index = 0.87), and GSE75538 dataset (n = 14, C-index = 0.48). Importantly, subtype Sub1 demonstrated a lower probability of survival and thus a more aggressive nature compared with subtype Sub2. ISGs in subtype Sub1 were enriched in the tumor-infiltrating immune cells-related pathways and cancer progression-related pathways, while those in subtype Sub2 were enriched in the metabolism-related pathways. Conclusion: The two survival subtypes of OSCC identified by deep learning can benefit clinical practitioners to divide immunocompromised patients with oral cancer into two subpopulations and give them target drugs and thus might be helpful for improving the survival of these patients and providing novel therapeutic strategies in the precision medicine area.

NMFNA: A Non-negative Matrix Factorization Network Analysis Method for Identifying Modules and Characteristic Genes of Pancreatic Cancer

Pancreatic cancer (PC) is a highly fatal disease, yet its causes remain unclear. Comprehensive analysis of different types of PC genetic data plays a crucial role in understanding its pathogenic mechanisms. Currently, non-negative matrix factorization (NMF)-based methods are widely used for genetic data analysis. Nevertheless, it is a challenge for them to integrate and decompose different types of genetic data simultaneously. In this paper, a non-NMF network analysis method, NMFNA, is proposed, which introduces a graph-regularized constraint to the NMF, for identifying modules and characteristic genes from two-type PC data of methylation (ME) and copy number variation (CNV). Firstly, three PC networks, i.e., ME network, CNV network, and ME-CNV network, are constructed using the Pearson correlation coefficient (PCC). Then, modules are detected from these three PC networks effectively due to the introduced graph-regularized constraint, which is the highlight of the NMFNA. Finally, both gene ontology (GO) and pathway enrichment analyses are performed, and characteristic genes are detected by the multimeasure score, to deeply understand biological functions of PC core modules. Experimental results demonstrated that the NMFNA facilitates the integration and decomposition of two types of PC data simultaneously and can further serve as an alternative method for detecting modules and characteristic genes from multiple genetic data of complex diseases.

Worldwide prevalence of PI3K-AKT-mTOR pathway mutations in head and neck cancer: A systematic review and meta-analysis

A systematic review (SR) and meta-analysis were conducted to determine the prevalence of PI3K-AKT-mTOR signaling pathway mutations in patients with head and neck cancer (HNC). Overall, 105 studies comprising 8630 patients and 1306 mutations were selected. The estimated mutations prevalence was 13 % for PIK3CA (95 % confidence interval [CI] = 11-14; I² = 82 %; p < 0.0001), 4% for PTEN (95 % CI = 3-5; I² = 55 %; p < 0.0001), 3% for MTOR (95 % CI = 2-4; I² = 5%; p = 0.40), and 2% for AKT (95 % CI = 1-2; I² = 50 %; p = 0.0001). We further stratified the available data of the participants according to risk factors and tumor characteristics, including HPV infection, tobacco use, alcohol exposure, TNM stage, and histological tumor differentiation, and performed subgroup analysis. We identified significant associations between PI3K-AKT-mTOR pathway-associated mutations and advanced TNM stage (odds ratio [OR] = 0.20; 95 % CI = 0.09-0.44; I² = 71 %; p = 0.0001) and oropharyngeal HPV-positive tumors and PIK3CA mutations (OR = 17.48; 95 % CI = 4.20-72.76; I² = 69 %; p < 0.0002). No associations were found between alcohol and tobacco exposure, and tumor differentiation grade. This SR demonstrated that the PI3K-AKT-mTOR pathway emerges as a potential prognostic factor and could offer a molecular basis for future studies on therapeutic targeting in HNC patients.

Biomolecules as Model Indicators of In Vitro and In Vivo Cold Plasma Safety

The potential applications for cold plasma in medicine are extensive, from microbial inactivation and induction of apoptosis in cancer cells to stimulating wound healing and enhancing the blood coagulation cascade. The safe bio-medical application of cold plasma and subsequent effect on complex biological pathways requires precision and a distinct understanding of how physiological redox chemistry is manipulated. Chemical modification of biomolecules such as carbohydrates, proteins, and lipids treated with cold plasma have been characterized, however, the context of how alterations of these molecules affect cell behavior or in vivo functionality has not been determined. Thus, this study examines the cytotoxic and mutagenic effects of plasma-treated molecules in vitro using CHO-K1 cells and in vivo in Galleria mellonella larvae. Specifically, albumin, glucose, cholesterol, and arachidonic acid were chosen as representative biomolecules, with established involvement in diverse bioprocesses including; cellular respiration, intracellular transport, cell signaling or membrane structure. Long- and short-term effects depended strongly on the molecule type and the treatment milieu indicating the impact of chemical and physical modifications on downstream biological pathways. Importantly, absence of short-term toxicity did not always correlate with absence of longer-term effects, indicating the need to comprehensively assess ongoing effects for diverse biological applications.

Identification of core miRNA prognostic markers in patients with laryngeal cancer using bioinformatics analysis

PURPOSE

Lots of studies indicated that many microRNAs (miRNAs) are associated with the prognosis of patients with laryngeal cancer (LC). The objective of our study is to identify potential core miRNAs associated with the pathogenesis and prognosis of LC.

METHODS

Using the Cancer Genome Atlast database, we identified 70 differentially expressed miRNAs between LC tumor specimens and non-tumor specimens. Then Cox regression analyses and the least absolute shrinkage and selection operator regression signature were performed to detect miRNA prognostic markers. A nomogram integrating miRNA prognostic markers was constructed to predict overall survival (OS) for LC patients. The potential target genes of the key miRNA were predicted by miRTarBase and miRDB databases. Subsequently, their potential functions were revealed by gene ontology annotation and kyoto encyclopedia of genes and genomes pathway enrichment analysis. Related biological pathways of the key target gene involved in LC were detected through gene set enrichment analysis (GSEA).

RESULTS

A prognostic miRNA signature was constructed. The up-regulated miR-105-1 was related to a worse OS (p = 0.043), which suggested that miR-105-1 may likely be the key miRNA prognostic marker. Survival analyses and paired expression analyses of target genes indicated that ENDOU may be the key target gene. Finally, we conducted GSEA to elucidate the pathways enriched between low- and high-ENDOU expression datasets.

CONCLUSION

Our findings might bring some new light on the pathogenesis of LC. Then, it might facilitate doctors to predict the prognosis and improve treatment outcomes for LC patients. However, the behaviors of LC are relatively heterogeneous, and the TCGA database cannot provide detailed information about the subsites and treatment modalities of LC. Further molecular biological experiments and clinical investigations would be required to confirm this conclusion.

Cytochrome P450 2U1 Is a Novel Independent Prognostic Biomarker in Breast Cancer Patients

Background: The susceptibility of breast cancer is largely affected by the metabolic capacity of breast tissue. This ability depends in part on the expression profile of cytochrome P450 (CYPs). CYPs are a superfamily of enzymes with related catalysis to endogenous and exogenous bioactive substances, including xenobiotic metabolism, drugs, and some endogenous substances metabolism which activate cells and stimulate cell signaling pathways, such as arachidonic acid metabolism, steroid metabolism, fatty acid metabolism. Interestingly, CYP was electively expressed in different tumors, and mediated the metabolic activation of multiple carcinogens and participated in the activation and deactivation of tumor therapeutic drugs. However, the biological action of cytochrome P450 2U1 (CYP2U1) in breast carcinoma is little understood so far.

Methods: To investigate the biological value of CYP2U1 in breast carcinoma, we performed immunohistochemical (IHC) analysis and survival analysis based on clinico-pathological data of breast cancer.

Results: IHC analysis showed that the abundance of CYP2U1 protein was inversely proportional to the state of estrogen receptor(ER) (P < 0.05), and the lower the degree of tumor differentiation, the higher the protein abundance (P < 0.001). Additionally, compared with luminal tumors, the CYP2U1 protein content was more abundant in triple negative breast cancer (P < 0.05). Importantly, survival analysis showed that higher CYP2U1 protein levels predicted poor 5-year overall survival rate (P < 0.01), 5-year disease-free survival rate (P < 0.05), and 5-year metastatic-free survival rate (P < 0.01) for the entire enrolled breast cancer patients.

Conclusions: CYP2U1 is generally closely related to the clinicopathological characteristics and is also an adverse prognostic factor for breast carcinoma patients, indicating that CYP2U1 is engaged in the malignant progression of breast carcinoma.

Identification of cancer stem cell characteristics in liver hepatocellular carcinoma by WGCNA analysis of transcriptome stemness index

Cancer stem cells (CSCs) are characterized by self-renewal and -differential potential as compared to common cancer cells and play an important role in the development and therapeutic resistance of liver hepatocellular carcinoma (LIHC). However, the specific pathogenesis of LIHC stem cells is still unclear, and the genes involved in the stemness of LIHC stem cells are currently unknown. In this study, we investigated novel biomarkers associated with LIHC and explored the expression characteristics of stem cell-related genes in LIHC. We found that mRNA expression-based stemness index (mRNAsi) was significantly overexpressed in liver cancer tissues. Further, mRNAsi expression in LIHC increased with the tumor pathological grade, with grade 4 tumors harboring the greatest stem cell features. Upon establishing mRNAsi scores based on mRNA expression of every gene, we found an association with poor overall survival in LIHC. Moreover, modules of interest were determined based on weighted gene co-expression network analysis (WGCNA) inclusion criteria, and three significant modules (red, green, and brown) and 21 key genes (DCN, ECM1, HAND2, PTGIS, SFRP1, SRPX, COLEC10, GRP182, ADAMTS7, CD200, CDH11, COL8A1, FAP, LZTS1, MAP1B, NAV1, NOTCH3, OLFML2A, PRR16, TMEM119, and VCAN) were identified. Functional analysis of these 21 genes demonstrated their enrichment in pathways involved in angiogenesis, negative regulation of DNA-binding transcription factor activity, apoptosis, and autophagy. Causal relationship with proteins indicated that the Wnt, Notch, and Hypoxia pathways are closely related to LIHC tumorigenesis. To our knowledge, this is the first report of a novel CSC biomarker, mRNAsi, to predict the prognosis of LIHC. Further, we identified 21 key genes through mRNA expression network analysis, which could be potential therapeutic targets to inhibit the stemness of cancer cells in LIHC.

Epigenomic dysregulation-mediated alterations of key biological pathways and tumor immune evasion are hallmarks of gingivo-buccal oral cancer

BACKGROUND

Gingivo-buccal oral squamous cell carcinoma (OSCC-GB) is the most common cancer among men in India and is associated with high mortality. Although OSCC-GB is known to be quite different from tongue cancer in its genomic presentation and its clinical behavior, it is treated identically as tongue cancer. Predictive markers of prognosis and therapy that are specific to OSCC-GB are, therefore, required. Although genomic drivers of OSCC-GB have been identified by whole exome and whole genome sequencing, no epigenome-wide study has been conducted in OSCC-GB; our study has filled this gap, and has discovered and validated epigenomic hallmarks of gingivobuccal oral cancer.

METHODS

We have carried out integrative analysis of epigenomic (n = 87) and transcriptomic (n = 72) profiles of paired tumor-normal tissues collected from OSCC-GB patients from India. Genome-wide DNA methylation assays and RNA-sequencing were performed on high-throughput platforms (Illumina) using a half-sample of randomly selected patients to discover significantly differentially methylated probes (DMPs), which were validated on the remaining half-sample of patients.

RESULTS

About 200 genes showed significant inverse correlation between promoter methylation and expression, of which the most significant genes included genes that act as transcription factors and genes associated with other cancer types. Novel findings of this study include identification of (a) potential immunosuppressive effect in OSCC-GB due to significant promoter hypomethylation driven upregulation of CD274 and CD80, (b) significant dysregulation by epigenetic modification of DNMT3B (upregulation) and TET1 (downregulation); and (c) known drugs that can reverse the direction of dysregulation of gene expression caused by promoter methylation.

CONCLUSIONS

In OSCC-GB patients, there are significant alterations in expression of key genes that (a) regulate normal cell division by maintenance of balanced DNA methylation and transcription process, (b) maintain normal physiological signaling (PPAR, B cell receptor) and metabolism (arachidonic acid) pathways, and (c) provide immune protection against antigens, including tumor cells. These findings indicate novel therapeutic targets, including immunotherapeutic, for treatment of OSCC-GB.

Fatty acids distribution and content in oral squamous cell carcinoma tissue and its adjacent microenvironment

Squamous cell carcinoma of the oral cavity mucosa grows under conditions of poor oxygenation and nutrient scarcity. Reprogramming of lipid biosynthesis accompanies tumor growth, but the conditions under which it occurs are not fully understood. The fatty acid content of the serum, tumor tissue and adjacent tumor microenvironment was measured by gas chromatography in 30 patients with squamous cell carcinoma grade 1-3. Twenty-five fatty acids were identified; their frequencies and percentages in each of the environments were assessed. Nineteen of the twenty-five fatty acids were found in tumor tissue, tumor adjacent tissue and blood serum. Of them, 8 were found in all thirty patients. Percentages of C16:0 and C18:1n9 were highest in the tumor, C18:1n9 and C16:0 were highest in tumor adjacent tissue, and C16:0 and C18:0 were highest in blood serum. The frequencies and amounts of C22:1n13, C22:4n6, C22:5n3 and C24:1 in tumor adjacent tissues were higher than those in blood serum, independent of the tumor grade. The correlations between the amount of fatty acid and tumor grade were the strongest in tumor adjacent tissues. The correlations between particular fatty acids were most prevalent for grade 1+2 tumors and were strongest for grade 3 tumors. In the adjacent tumor microenvironment, lipogenesis was controlled by C22:6w3. In blood serum, C18:1trans11 limited the synthesis of long-chain fatty acids. Our research reveals intensive lipid changes in oral cavity SCC adjacent to the tumor microenvironment and blood serum of the patients. Increase in percentage of some of the FAs in the path: blood serum-tumor adjacent microenvironment-tumor, and it is dependent on tumor grade. This dependency is the most visible in the tumor adjacent environment.

Establishment and genomic characterization of gingivobuccal carcinoma cell lines with smokeless tobacco associated genetic alterations and oncogenic PIK3CA mutation

Smokeless tobacco associated Gingivobuccal squamous cell carcinoma (GB-SCC) is a major public health problem but available oral cancer cell lines are mostly from smoking associated tongue SCC raising the need for pertinent GB-SCC cell line models. As part of the International Cancer Genome Consortium (ICGC) Project, 4 novel cell lines, namely, Indian Tata Memorial Centre Oral Cancer (ITOC) -01 to -04 were established and characterized with conventional methods, karyotyping, ultrastructure, in vivo tumourigenicity, Whole exome sequencing (WES) and RNA sequencing. These hyperploid cell lines form xenografts in mice and show metabolically active and necrotic areas on fluorodeoxyglucose-positron emission tomography (FDG-PET) imaging. WES of ITOC cell lines recapitulate the genomic tumor profile of ICGC GB-SCC database. We further identified smokeless tobacco associated genetic alterations (PCLO, FAT3 and SYNE2) and oncogenic PIK3CA mutation in GB-SCC cell lines. Transcriptome profiling identified deregulation of pathways commonly altered in cancer and down-regulation of arachidonic acid metabolism pathway, implying its possible role in GB-SCC. Clinical application of high throughput sequencing data depends on relevant cell line models to validate potential targets. Extensively characterized, these oral SCC cell lines are particularly suited for mechanistic studies and pre-clinical drug development for smokeless tobacco associated oral cancer.

Harnessing genomics to improve outcomes for women with cancer in India: key priorities for research

Cumulatively, breast, cervical, ovarian, and uterine cancer account for more than 70% of cancers in women in India. Distinct differences in the clinical presentation of women with cancer suggest underlying differences in cancer biology and genetics. The peak age of onset of breast and ovarian cancer appears to be a decade earlier in India (age 45-50 years) than in high-income countries (age >60 years). Understanding these differences through research to develop diagnosis, screening, prevention, and treatment frameworks that ar e specific to the Indian population are critical and essential to improving women's health in India. Since the sequencing of the human genome in 2001, applications of advanced technologies, such as massively parallel sequencing, have transformed the understanding of the genetic and environmental drivers of cancer. How can advanced technologies be harnessed to provide health-care solutions at a scale and to a budget suitable for a country of 1·2 billion people? What research programmes are necessary to answer questions specific to India, and to build capacity for innovative solutions using these technologies? In order to answer these questions, we convened a workshop with key stakeholders to address these issues. In this Series paper, we highlight challenges in tackling the growing cancer burden in India, discuss ongoing genomics research and developments in infrastructure, and suggest key priorities for future research in cancer in India.

Identification of differentially expressed genes in the development of osteosarcoma using RNA-seq

Objective

Osteosarcoma (OS) is a malignant bone tumor with high morbidity in young adults and adolescents. This study aimed to discover potential early diagnosis biomarkers in OS.

Results

In total, 111 differentially expressed genes (DEGs) were identified in primary OS compared with normal controls and 235 DEGs were identified in metastatic OS compared with primary OS. AURKB and PPP2R2B were the significantly up-regulated and down-regulated hub proteins, respectively, in the PPI protein-protein network (PPI) network of primary OS. ISG15 and BTRC were the significantly up-regulated and down-regulated hub proteins, respectively, in the network of metastatic OS. The DEGs in metastatic OS compared with primary OS were significantly enriched in the arachidonic acid metabolism, malaria, and chemokine signaling pathways. Finally, we employed quantitative real-time polymerase chain reaction (qRT-PCR) to validate the expression levels of candidate DEGs and the results indicated that our bioinformatics approach was acceptable.

Materials and Methods

The mRNA expression profiling of 20 subjects was obtained through high-throughput RNA-sequencing. DEGs were identified between primary OS and normal Control, and between primary OS and metastatic OS, respectively. Functional annotation and PPI networks were used to obtain insights into the functions of DEGs. qRT-PCR was performed to detect the expression levels of dysregulated genes in OS.

Conclusions

Our work might provide groundwork for the further exploration of tumorigenesis and metastasis mechanisms of OS.

Genomic characterization of tobacco/nut chewing HPV-negative early stage tongue tumors identify MMP10 asa candidate to predict metastases

OBJECTIVES

Nodal metastases status among early stage tongue squamous cell cancer patients plays a decisive role in the choice of treatment, wherein about 70% patients can be spared from surgery with an accurate prediction of negative pathological lymph node status. This underscores an unmet need for prognostic biomarkers to stratify the patients who are likely to develop metastases.

MATERIALS AND METHODS

We performed high throughput sequencing of fifty four samples derived from HPV negative early stage tongue cancer patients habitual of chewing betel nuts, areca nuts, lime or tobacco using whole exome (n=47) and transcriptome (n=17) sequencing that were analyzed using in-house computational tools. Additionally, gene expression meta-analyses were carried out for 253 tongue cancer samples. The candidate genes were validated using qPCR and immuno-histochemical analysis in an extended set of 50 early primary tongue cancer samples.

RESULTS AND CONCLUSION

Somatic analysis revealed a classical tobacco mutational signature C:G>A:T transversion in 53% patients that were mutated in TP53, NOTCH1, CDKN2A, HRAS, USP6, PIK3CA, CASP8, FAT1, APC, and JAK1. Similarly, significant gains at genomic locus 11q13.3 (CCND1, FGF19, ORAOV1, FADD), 5p15.33 (SHANK2, MMP16, TERT), and 8q24.3 (BOP1); and, losses at 5q22.2 (APC), 6q25.3 (GTF2H2) and 5q13.2 (SMN1) were observed in these samples. Furthermore, an integrated gene-expression analysis of 253 tongue tumors suggested an upregulation of metastases-related pathways and over-expression of MMP10 in 48% tumors that may be crucial to predict nodal metastases in early tongue cancer patients. In overall, we present the first descriptive portrait of somatic alterations underlying the genome of tobacco/nut chewing HPV-negative early tongue cancer, and identify MMP10 asa potential prognostic biomarker to stratify those likely to develop metastases.

Analysis of the whole transcriptome from gingivo-buccal squamous cell carcinoma reveals deregulated immune landscape and suggests targets for immunotherapy

Background

Gingivo-buccal squamous cell carcinoma (GBSCC) is one of the most common oral cavity cancers in India with less than 50% patients surviving past 5 years. Here, we report a whole transcriptome profile on a batch of GBSCC tumours with diverse tobacco usage habits. The study provides an entire landscape of altered expression with an emphasis on searching for targets with therapeutic potential.

Methods

Whole transcriptomes of 12 GBSCC tumours and adjacent normal tissues were sequenced and analysed to explore differential expression of genes. Expression changes were further compared with those in TCGA head and neck cohort (n = 263) data base and validated in an independent set of 10GBSCC samples.

Results

Differentially expressed genes (n = 2176) were used to cluster the patients based on their tobacco habits, resulting in 3 subgroups. Immune response was observed to be significantly aberrant, along with cell adhesion and lipid metabolism processes. Different modes of immune evasion were seen across 12 tumours with up-regulation or consistent expression of CD47, unlike other immune evasion genes such as PDL1, FUT4, CTLA4 and BTLA which were downregulated in a few samples. Variation in infiltrating immune cell signatures across tumours also indicates heterogeneity in immune evasion strategies. A few actionable genes such as ITGA4, TGFB1 and PTGS1/COX1 were over expressed in most samples.

Conclusion

This study found expression deregulation of key immune evasion genes, such as CD47 and PDL1, and reasserts their potential as effective immunotherapeutic targets for GBSCC, which requires further clinical studies. Present findings reiterate the idea of using transcriptome profiling to guide precision therapeutic strategies.

Cytochrome P450 2U1, a very peculiar member of the human P450s family

Cytochrome P450 2U1 (CYP2U1) exhibits several distinctive characteristics among the 57 human CYPs, such as its presence in almost all living organisms with a highly conserved sequence, its particular gene organization with only five exons, its major location in thymus and brain, and its protein sequence involving an unusually long N-terminal region containing 8 proline residues and an insert of about 20 amino acids containing 5 arginine residues after the transmembrane helix. Few substrates, including fatty acids, N-arachidonoylserotonin (AS), and some drugs, have been reported so far. However, its biological roles remain largely unknown, even though CYP2U1 mutations have been involved in some pathological situations, such as complicated forms of hereditary spastic paraplegia. These data together with its ability to hydroxylate some fatty acids and AS suggest its possible role in lipid metabolism.

Transcriptomic characterization of differential gene expression in oral squamous cell carcinoma: a meta-analysis of publicly available microarray data sets

Oral squamous cell carcinoma (OSCC) is a highly prevalent cancer worldwide, and OSCC often goes undiagnosed until advanced disease is present, which contributes to a low survival rate for OSCC patients. The identification of biomarkers for the early detection OSCC and novel therapeutic targets for OSCC treatment is an important research objective. We performed bioinformatics analyses of the gene expression profile of OSCC using microarray data to identify genes that contribute to the development of OSCC. We also predicted the transcription factors involved in the regulation of differential gene expression in OSCC. Our results showed that PI3K, EGFR, STAT1, and CPBP are important contributors to the changes in cellular physiology that occur during the development of OSCC. Therefore, these genes represent potential diagnostic biomarkers and therapeutic targets for OSCC.