In the last decade, different omic technologies have experiences an exponential technological advancement. However, metabolomics has not followed a similarly vertiginous improving improvement and are far from genomics or transcriptomics in terms of throughput, cost and even accuracy. Therefore, genome-scale in-silico methodologies to estimate metabolic activities from genomic data constitute an active field. The solutions available fall into two extremes: those with few assumptions about the relationships among proteins and metabolites, which are easy-to-use but less accurate and those that account for the complex relationships among molecules and proteins defined in the metabolic pathways, which are more accurate but require mathematical skills. Here, we introduce Metabolica, an algorithm that considers the complex functional relationships among all the molecules and proteins involved in the metabolic pathway analyzed but keeping an easy use that do not require of advanced mathematical skills. Metabolica has been implemented in a freely available software R package. The software inputs transcriptomic data and infers the activities of the reactions that produce the different metabolites in the pathway analyzed. An example shows how detected dysregulated metabolites in several cancers are related to patient survival.

The methods used to introduce CRISPR (Clustered Regularly Interspaced Short Palindromic Repeats)/Cas-mediated genome editing into fruit species, as well as the impacts of the application of this technology to activate and knock out of target genes in different fruit trees species including tree development, yield, fruit quality, and tolerance to biotic and abiotic stresses have been firstly described in this review. The application of this gene editing technology could allow the development of new generations of fruit crops with improved traits by targeting different genetic segments or even could facilitate the introduction of traits in elite cultivars without changing other traits. However, at this moment, the scarcity of efficient regeneration and transformation protocols in some species, the fact that many of those procedures are the genotype-dependent or the convenience of segregating the transgenic parts of the CRISPR system represent the main handicaps limiting the potential of genetic editing techniques for fruit trees. Finally, latest news on the legislation and regulations about the use of plants modified through CRISPR/Cas systems has been also discussed.

Biofouling is the growth of organisms on wetted surfaces. Biofouling includes micro- (bacteria and unicellular algae) and macrofouling (mussels, barnacles, tube worms, bryozoans, etc.) and is a major problem for industries. However, the settlement and growth of some biofouling species, like oysters and corals, can be desirable. Thus, it is important to understand the process of biofouling in detail. Modern “omic” techniques, metabolomics, metagenomics, transcriptomics, and proteomics provide unique opportunities to study biofouling organisms and communities and investigate their metabolites and interactions with the environment. In this review, we analyze the recent publications that employ metagenomic, metabolomic, and proteomic techniques for the investigation of biofouling and biofouling organisms. Specific emphasis is given to metagenomics and proteomics as well as those publications that use a combination of different “omics” techniques. Finally, this review presents the future outlook for the use of “omics” techniques in marine biofouling studies. Like all trans-disciplinary research, environmental “omics” is in its infancy and will advance rapidly as researchers develop the necessary expertise, theory, and technology.

Glioma accounts for 80% of all malignant brain tumors and is the most common adult primary brain tumor. Age is an important factor affecting the development of cancer as somatic mutations accumulate with age. In this study, we aimed to analyze the significance of age-related non-silent somatic mutations in glioma prognosis. Histological tumor grade depends on age at diagnosis in patients with IDH1, TP53, ATRX, and EGFR mutations. The hierarchical clustering of patients was dominantly separated by IDH1 and EGFR mutations. Furthermore, patients with IDH1 mutation were dominantly separated by TP53 and ATRX double mutation and its double wildtype counterpart. Patients with the double mutation showed poorer prognosis than those with the double wild type genotype. In conclusion, among the many somatic mutations, those in IDH1, TP53, ATRX, and EGFR are important for glioma classification based on histological grade. Patients with EGFR mutation had the poorest prognosis, whereas those with only IDH1 mutation showed the best prognosis.

Vibrio splendidus is the main pathogen infect a wide range of hosts, especially sea cucumber Apostichopus japonicus. Our previous studies had shown that the level of L-glutamic acid (L-Glu) was significantly increased under heat stress and it was one of the best carbon sources that could be used by V. splendidus AJ01. In this study, the effects of exogenous L-Glu on the coelomocyte viability, tissue status and individual mortality of sea cucumbers were analyzed. The results showed that 10 mM L-Glu decreased coelomocyte viability and increased mortality, with tissue rupture and pyknosis, while 0.1 mM L-Glu slightly affected the survival of sea cucumbers without obvious damages at both cellular and tissue levels. Transcriptomic analysis showed that there were 343 upregulated and 206 downregulated genes regulated by exogenous L-Glu. GO analysis showed that differentially expressed genes (DEGs) were mainly enriched in signaling and membrane formation, while KEGG analysis showed that DEGs were significantly enriched in the upregulated endocytosis and downregulated lysosomal pathways. The coelomocyte viability further decreased 20% in the simultaneous presence of exogenous L-Glu and V. splendidus AJ01, compared to that of V. splendidus AJ01 infection alone. Consequently, a higher mortality of sea cucumber was also observed in the presence of exogenous L-Glu challenged by V. splendidus AJ01. Real time reverse transcriptase PCR showed that L-Glu specifically upregulated the expression of the fliC gene, conding the subunit protein of flagellar filament, with promoted swimming motility activity of V. splendidus. Our result indicated that L-Glu should be kept in a state of equilibrium, and the excess L-Glu in the host-pathogen interface would prompt the virulence of V. splendidus via increasing bacterial motility.

Adipose tissue composition contributes greatly to the quality and nutritional value of meat. Transcriptomic and lipidomic techniques were used to investigate the molecular mechanism of fat deposition difference among Ningxiang, Berkshire and F1 pigs. Transcriptomic analysis identified 680, 592 and 380 differentially expressed genes (DEGs) by comparing the groups of Ningxiang pigs vs Berkshire pigs, Berkshire pigs vs F1 pigs, Ningxiang pigs vs F1 pigs. Lipidomic analysis screened 423, 252, and 50 significantly changed lipids (SCLs) by comparing the groups of Ningxiang pigs vs Berkshire pigs, Berkshire pigs vs F1 pigs, Ningxiang pigs vs F1 pigs. Lycine, serine and threonine metabolism pathway, fatty acid biosynthesis and metabolism related pathways were siginificant enriched in the groups of Berkshire pigs vs Ningxiang pigs and Berkshire pigs vs F1 pigs. The DEGs (PHGDH, LOC110256000) and the SCLs (Phosphatidylserines) may have a great impact on lycine, serine and threonine metabolism pathway. Moreover, the DEGs (FASN, ACACA, CBR4, SCD, ELOV6, HACD, CYP3A46, CYP2B22, GPX1, GPX3) and the SCLs (palmitoleic acid, linoleic acid, arachidonic acid, icosadienoic acid) play important role in fatty acid biosynthesis and metabolism. Thus, the fat deposition difference among Ningxiang, Berkshire and F1 pigs may be caused by the difference of the expression pattern of key genes in multiple enrichment KEGG pathways. The research presented multiple lipids that were potentially available biological indicator and screened key genes that were potentially gene targets for molecular design breeding. The research also explored the molecular mechanisms of fat deposition difference among Ningxiang, Berkshire and F1 pigs and provides insight into the selection of backfat thickness and fat composition of adipose tissue for future breeding strategies.

This study analyzes oral supplement of molecular iodine (I2) alone and in combination with the neoadjuvant therapy 5-fluorouracil/epirubicin/cyclophosphamide or taxotere/epirubicin (FEC/TE) in women with Early (stage II) and Advanced (stage III) breast cancer. In the Early group, 30 women were treated with I2 (5 mg/day) or placebo (colored water) for 7–35 days before surgery. For the Advanced group, 30 patients received I2 or placebo along with FEC/TE treatment. After surgery, all patients received FEC/TE + I2 for 170 days. I2 supplementation showed a significant attenuation of the side effects and absence of tumor chemoresistance. Control, I2, FEC/TE, and FEC/TE+I2 groups exhibit response rate of 0, 33%, 73%, and 100%, respectively, and a pathologic complete response of 18%, and 36% in the last two groups. Five-year disease-free survival rate was significantly higher in patients with I2 supplement before and after surgery compared to those receiving the supplement only after surgery (82% vs. 46%). I2-treated tumors exhibit less invasive potential, and significant increases in apoptosis, estrogen receptor expression, and immune cell infiltration. Transcriptomic analysis indicated activation of the antitumoral immune response. The results led us to register a phase III clinical trial analyzing chemotherapy + I2 treatment for advanced breast cancer.

Lentiviral vectors (LVs) are a powerful tool for gene and cell therapy and human embryonic kidney cells (HEK293) have been extensively used as a platform for production of these vectors. Like most cells and cellular tissues, HEK293 cells release extracellular vesicles (EVs). EVs released by cells share similar size, biophysical characteristics and even a biogenesis pathway with cell-produced enveloped viruses, making it a challenge to efficiently separate EVs from LVs. Thus, EVs co-purify with LVs during downstream processing, becoming “impurities” in the context of cell therapy. To characterize EVs from an inducible lentivirus producing cell line, two conditions were studied: non-induced and induced. EVs’ identity was confirmed by transmission electron microscopy and western blot. Seven proteins were identified by mass spectrometry as potential EV markers. Lipid composition of EVs and LVs showed similar enrichment in phosphatidylserine. RNA cargos in EVs showed enrichment in genes involved in viral processes and binding functions. Flow virometry, GTA and ddPCR results also confirmed the heterogenic nature of EVs and LVs populations. These findings provide insights on the product profile of lentiviral preparation and could help develop separation strategies of co-produced EVs.

The Loropetalum chinense and Loropetalum chinense var. rubrum are typical and traditional ornamental and Chinese herbal medicine in Asia. However, more information is needed on the mechanisms underlying its flower coloring. Here, we profiled the flavonoid metabolome, full-length sequencing, and transcriptome analysis to investigate the flavonoid biosynthesis and global transcriptome changes among different petal coloring cultivars of L. chinense and L. chinense var. rubrum. The total anthocyanins and phenotypic of the petal were highly consistent with the petal color. Moreover, a total of 207 flavonoid components were identified. Of these, 12 flavonoid components were considered significantly different expression compounds among the four samples. Meanwhile, the first reference full-length transcriptome of L. chinense var. rubrum was being built, which had 171,783 high-quality non-redundant transcripts with correcting with next-generation sequencing (NGS). Among them, 52,851 transcripts were annotated in the seven database of NR, KOG, GO, NT, Pfam, Swiss-Port and KEGG. Combined with NGS analysis, the DETs involved in flavonoids and anthocyanins contributed greatest to the flower-coloring. Additionally, the different expressed of eight LcDFRs and four LcANS genes were positively correlated with flavonoid biosynthesis, and the four LcBZ1 and one Lc3Mat1 were positively correlated with the content of seven anthocyanins revealed by coupling with metabolomics and transcriptomics analysis. Together these results were used to mine candidate genes by analyzing flower coloring changes in a comprehensive metabolic and transcriptomic level in L. chinense and L. chinense var. rubrum.

Prefrontal cortex (PFC) is one of the brain regions with more prominent changes in human aging. The molecular processes related to the aging cognitive decline and mood changes are not completely understood. In order to improve our knowledge, we integrated transcriptomic data of four studies of human PFC from old people -58-80 years old- compared with young people -20-40 years old- using a meta-analytic approximation combined with molecular signature analysis. We identified 1816 differentially expressed genes -561 up-regulated and 1256 down-regulated. Pathway analysis revealed down-regulation of synaptic genes with conservation of gene expression of other neuronal regions. Additionally, we identified up-regulation of markers of astrogliosis with transcriptomic signature compatible with A1 neurotoxic astrocytes and A2 neuroprotective astrocytes. Response to interferon is related to A1 astrocytes and the A2 phenotype is mediated in aging by activation of SHH pathway and up-regulation of metallothioneins I and genes of the family EZR -ezrin, radixin, and moesin-. The main conclusions of our study are the confirmation that in aged PFC there is a global dysfunction of the synapses and we reported for the first time opposite phenotypes of astrogliosis because of brain aging.

Fragile X syndrome (FXS) is a neurodevelopmental disorder associated with a wide range of cognitive, behavioral and medical problems. It arises from the silencing of the fragile X mental retardation 1 (FMR1) gene, and consequently, in the absence of its encoded protein, FMRP (Fragile X Mental Retardation Protein). FMRP is a ubiquitously expressed and multifunctional RNA-binding protein, primarily considered as a translational regulator. Pre-clinical studies of the past two decades have therefore focus on this function to relate FMRP’s absence to the molecular mechanisms underlying FXS physiopathology. Based on these data, successful pharmacological strategies were developed to rescue fragile X phenotype in animal models. Unfortunately, these results did not translate into human, as clinical trials using same therapeutic approaches did not reach the expected outcomes. These failures highlight the need to put into perspectives the different functions of FMRP in order to get a more comprehensive understanding of FXS pathophysiology. In this review, FMRP’s involvement on noteworthy molecular mechanisms are pointed out; ultimately contributing to various biochemicals alterations composing the fragile X phenotype.

An outbreak of a virus-like disease has caused severe damage to noni plants (Morinda citrifolia L.) in Xishuangbanna area of China's southwestern Yunnan province since 2015. The diseased plants displayed typical mosaic symptom with light and dark green patches on leaves. Flexuous filamentous virus particles of about 800 nm in length were observed from the leaf saps by transmission electron microscope. Illumina transcriptomic sequencing further revealed the presence of a potyvirus and its near complete genome was obtained from de novo assembly. The complete genome of 9,659 nts was obtained by Sanger sequencing of eight amplicons generate by RT-PCR and 5’ and 3’ RACE. BLASTp analysis of the polyprotein sequence showed that the virus was most closely related to Tobacco vein banding mosaic virus (TVBMV), but these two viruses only shared 50.7% amino acid sequence similarity. Both phylogenetic analyses of the polyprotein and CP amino acid sequences indicated that this virus is a member of genus Potyvirus. However, the low sequence homology with all known potyviruses established this virus as a new species in the genus, tentatively named as Noni mosaic virus (NoMV). Our field surveys showed that 100% of the symptomatic samples and 28.57% of the asymptomatic samples were infected with this novel potyvirus. Aphids collected from diseased leaves were also detected carrying the virus. In summary, our data indicated that a novel species of potyvirus, NoMV, is prevalent in Yunnan, China and is associated with an emerging mosaic disease on M. citrifolia.

Background: The demand in the development of cancer nanomedicine has increased due to various limitations in conventional cancer therapy. This study assessed the mRNA transcriptomic profiling of human HepG2 cells exposed to C. roseus-AgNPs. Methods: The proliferative activity of hepatocellular carcinoma (HepG2) and normal human liver (THLE3) cells treated with C. roseus‑AgNPs were measured using MTT assay. The RNA samples were extracted and sequenced using BGIseq500 platform. This is followed by data filtering, mapping, gene expression analysis, DEG analysis, GO analysis, and pathway analysis. Results: The mean IC₅₀ values of C. roseus‑AgNPs on HepG2 was 4.38±1.59 µg/mL while on THLE3 cells was 800±1.55 µg/mL. Transciptomic profiling revealed an alteration of 296 genes. C. roseus‑AgNPs induced the expression of stress-associated genes such as MT, HSP and HMOX-1. Cellular signaling pathways were potentially activated through MAPK, TNF and TGF pathways that responsible for apoptosis and cell cycle arrest. The alteration of ARF6, EHD2, FGFR3, RhoA, EEA1, VPS28, VPS25, TSG101 indicated the uptake of C. roseus-AgNPs via both clathrin-dependent and clathrin-independent endocytosis. Conclusions: This study provides the new insights on gene expression study of biosynthesized AgNPs on cancer cells. The cytotoxicity effect is mediated by the aberrant gene alteration, and more interestingly the unique selective antiproliferative properties indicates the C. roseus‑AgNPs as an ideal anticancer candidate.

Decades of research identified numerous gene biomarkers of cardiac diseases whose restored sequence or/and expression level was hoped to recover the normal cardiac function. However, each human has unique and dynamic pathophysiological characteristics resulting from the unrepeatable combination of favoring factors such are: race, sex, age, medical history, diet, stress, exposure to toxins, habits etc. As such, no treatment fits everybody and finding personalized solutions is a top priority for medicine of 21st century. The Genomic Fabric Paradigm (GFP) provides the most theoretically possible comprehensive characterization of the transcriptome, its alterations in disease and recovery following a treatment. By attaching to each gene the independent average expression level, expression variation and expression coordination with each other gene, GFP delivers thousands times more information than the traditional analysis. This report presents the theoretical bases of the GFP and some applications to our microarray data from mouse models of post ischemic, and constant and intermittent hypoxia-induced heart failure. The GFP analyses revealed novel transcriptomic aspects of the gene expression control and networking under ischemic conditions. Through all-inclusive characterization of the transcriptome and the unrepeatable gene hierarchy in each condition, GFP is an essential avenue towards development of a truly personalized cardiogenomic therapy.

Prostate cancer is a leading cause of death among men but its genomic characterization and best therapeutic strategy are still under debate. The Genomic Fabric Paradigm (GFP) considers the transcriptome as a multi-dimensional mathematical object subjected to a dynamic set of expression correlations among the genes. Here, GFP is applied to gene expression profiles of three (one primary, and two secondary) cancer nodules and the surrounding normal tissue from a surgically removed prostate tumor. GFP was used to determine the regulation and rewiring of the P53 signaling, apoptosis, prostate cancer and several other pathways involved in survival and proliferation of the cancer cells. Genes responsible for the block of differentiation, evading apoptosis, immortality, insensitivity to anti-growth signals, proliferation, resistance to chemotherapy and sustained angiogenesis were found as differently regulated in the three cancer nodules with respect to the normal tissue. The analysis indicates that even histo-pathologically equally graded cancer nodules from the same tumor have substantially different transcriptomic organizations, raising legitimate questions about the validity of meta-analyses comparing large populations of healthy and cancer humans. The study suggests that GFP may provide a personalized alternative to the biomarkers’ approach of cancer genomics.

The plant graft healing process is an intricate development influenced by numerous endogenous and environmental factors. This process involves the histological changes, physiological and biochemical reactions, signal transduction, and hormone exchanges in the grafting junction. Studies have shown that applying exogenous plant growth regulators can effectively promote the graft healing process and improve the quality of grafted plantlets. However, the physiological and molecular mechanism of graft healing formation remains unclear. In our present study, transcriptome changes in the melon and cucurbita genomes were analyzed between control and NAA treatment, and we provided the first view of complex networks to regulate graft healing under exogenous NAA application. The results showed that the exogenous NAA application could accelerate the graft healing process of oriental melon scion grafted onto squash rootstock through histological observation, increase the SOD, POD, PAL, and PPO activities during graft union development and enhance the contents of IAA, GA3, and ZR except for the IL stage. The DEGs were identified in the plant hormone signal-transduction, phenylpropanoid biosynthesis, and phenylalanine metabolism through transcriptome analysis of CK vs. NAA at the IL, CA, and VB stage by KEGG pathway enrichment analysis. Moreover, the exogenous NAA application significantly promoted the expression of genes involved in the hormone signal-transduction pathway, ROS scavenging system, and vascular bundle formation.

Idiopathic lung fibrosis (IPF) is a progressive and fatal degenerative lung disease of unknown etiology. Although in its final stages it implicates in a reactive manner all lung cell types, the initial damage involves the alveolar epithelial compartment, in particular the alveolar epithelial type 2 cells (AEC2s). AEC2s serve dual progenitor and surfactant secreting functions, both of which are deeply impacted in IPF. Thus, we hypothesize that the size of the surfactant processing compartment, as measured by Lysotracker incorporation, allows the identification of different epithelial states in the IPF lung. Flow cytometry analysis of epithelial Lysotracker incorporation delineates two populations (Lysohigh and Lysolow) of AEC2s which behave in a compensatory manner during bleomycin injury and in the donor/IPF lung. Employing flow cytometry and transcriptomic analysis of cells isolated from donor and IPF lungs, we demonstrate that the Lysohigh population expresses all classical AEC2 markers and is drastically diminished in IPF. The Lysolow population, which is increased in proportion in IPF, co-expresses AEC2s and basal cell markers resembling the phenotype of the previously identified intermediate AEC2 population in the IPF lung. In that regard, we provide an in-depth flow-cytometry characterization of Lysotracker uptake, HTII-280, proSP-C, mature SP-B, NGFR, KRT5 and CD24 expression in human lung epithelial cells. Combining functional analysis with extra- and intra- cellular marker expression and transcriptomic analysis, we advance the current understanding of epithelial cell behavior and fate in lung fibrosis.

For the majority of lipophilic compounds adipose tissue is traditionally considered as storage depot and only rarely as a target organ. Meanwhile, abnormalities in adipose tissue physiology induced by chemical exposures may contribute to the current epidemic of obesity and metabolic diseases. Polybrominated diphenyl ethers (PBDEs) is a group of lipophilic flame retardants found in majority of human samples in North America. Their ability to alter physiology of adipose tissue is unknown. We exposed pregnant mice to 0.2 mg/kg body weight/day of BDE-47 perinatally. Transcriptomic changes in gonadal adipose tissue were analyzed in male offspring using RNA-seq approach with subsequent bioinformatic analysis. Genes of coagulation and complement cascade, de novo lipogenesis, and xenobiotic metabolism were altered in expression in response to BDE-47 exposure. The affected molecular network included the following hubs: PPARα, HNF1A and HNF4. These findings suggest that adipose tissue should be considered a target tissue for BDE-47, in addition to its role as a storage depot. This study also builds a background for a targeted search of sensitive phenotypic endpoints of BDE-47 exposure, including lipid profile parameters and coagulation factors in circulation. Additional studies are needed to investigate the role of PBDEs as an obesogen.

Pituitary adenylate cyclase-activating polypeptide receptor 1 (PAC1-R) has a higher affinity for pituitary adenylate cyclase- activating polypeptide (PACAP) than for vasoactive intestinal peptide (VIP). PACAP exerts its neuroprotective effects through PAC1-R. PAC1-R has become a target for drug development, and we targeted PAC1-R to screen small positive allosteric modulator 1 (SPAM1). SPAM1 has a neuroprotective effects, and in the D-galactose (D-gal) mouse model, SPAM1 improved the structure of the hippocampal dentate gyrus and restored the number of neurons. Compared with D-gal model mice, SPAM1 treatment upregulated the expression of Sirtuin 6 (SIRT6) and Lamin B1, and downregulated the expression of YinYang 1 (YY1) and p16. Consistent with this, the same results were obtained in senescent RGC-5 cells. Meanwhile, SPAM1 was also found to up-regulate SIRT6 and down-regulate YY1 expression in normal RGC-5 cells. Then, we further explored the mechanism of SPAM1 neuroprotection by whole transcriptome sequencing and proteomic analysis, and found that SPAM is involved in the longevity regulating pathway, and confirmed the up-regulation of neurofilament light polypeptide and neurofilament medium polypeptide by western blot, thus laying a pharmacological foundation for SPAM1 as a therapeutic agent for neurodegenerative diseases. polypeptide and neurofilament medium polypeptide. Thus, it lays a pharmacological foundation for SPAM1 as a therapeutic drug for neurodegenerative diseases.

De novo metastatic hormone-sensitive PC (mHSPC) accounts for 5-10% of all prostate cancer (PC) diagnoses but it is responsible for nearly 50% of PC-related deaths. Since 2015, the prognosis of mHSPC has slightly improved thanks to the introduction of new hormonal agents and chemotherapy combined with androgen deprivation therapy from the first-line setting. This review describes the current therapeutic opportunities in de novo mHSPC, focusing on potential molecular biomarkers identified in the main clinical trials that have changed the standard of care, the genomic features of de novo mHSPC, and the principal ongoing trials that are investigating new therapeutic approaches and the efficacy of a biomarker-guided treatment in this setting. The road towards a personalized treatment for de novo mHSPC is still long considering that the randomized clinical trials, which have furnished the basis of the current therapeutic options, stratified patients according to clinical criteria that not necessarily reflect the biological rationale of the chosen therapy. The role of transcriptomic profiling of mHSPC as predictive biomarker requires further validation, as well as it remains to ascertain how the genomic alterations detected in mHSPC, that are considered predictive in the castration-resistant disease, can be exploited in the mHSPC setting.

The survival analysis of the Cancer Genome Atlas (TCGA) dataset is a well-known method to discover the gene expression-based prognostic biomarkers of head and neck squamous cell carcinoma (HNSCC). A cutoff point is usually used in survival analysis for the patients' dichotomization in the continuous gene expression. There is some optimization software for cutoff determination. However, the software's predetermined cutoffs are usually set at the median or quantiles of gene expression value to perform the analyses. There are also few clinicopathological features available on their pre-processed data sets. We applied an in-house workflow, including data retrieving and pre-processing, feature selection, sliding-window cutoff selection, Kaplan-Meier survival analysis, and Cox proportional hazard modeling for biomarker discovery. In our approach for the TCGA HNSCC cohort, we scanned human protein-coding genes to find optimal cutoff values. After adjustment with confounders, the clinical tumor stage and the surgical margin involvement are independent risk factors for patients' prognosis. According to the resulting tables with Bonferroni-adjusted P value under the optimal cutoff and the hazard ratio, three biomarker candidates, CAMK2N1, CALML5, and FCGBP, are significantly associated with the patients' overall survival. We validated this discovery by using the other independent HNSCC dataset (GSE65858). Thus, we suggest the transcriptomic analysis could help for biomarker discovery.