mBAT-combo: A more powerful test to detect gene-trait associations from GWAS data.

Gene-based association tests aggregate multiple SNP-trait associations into sets defined by gene boundaries and are widely used in post-GWAS analysis. A common approach for gene-based tests is to combine SNPs associations by computing the sum of χ2 statistics. However, this strategy ignores the directions of SNP effects, which could result in a loss of power for SNPs with masking effects, e.g., when the product of two SNP effects and the linkage disequilibrium (LD) correlation is negative. Here, we introduce "mBAT-combo," a set-based test that is better powered than other methods to detect multi-SNP associations in the context of masking effects. We validate the method through simulations and applications to real data. We find that of 35 blood and urine biomarker traits in the UK Biobank, 34 traits show evidence for masking effects in a total of 4,273 gene-trait pairs, indicating that masking effects is common in complex traits. We further validate the improved power of our method in height, body mass index, and schizophrenia with different GWAS sample sizes and show that on average 95.7% of the genes detected only by mBAT-combo with smaller sample sizes can be identified by the single-SNP approach with a 1.7-fold increase in sample sizes. Eleven genes significant only in mBAT-combo for schizophrenia are confirmed by functionally informed fine-mapping or Mendelian randomization integrating gene expression data. The framework of mBAT-combo can be applied to any set of SNPs to refine trait-association signals hidden in genomic regions with complex LD structures.

Machine learning enables pan-cancer identification of mutational hotspots at persistent CTCF binding sites.

CCCTC-binding factor (CTCF) is an insulator protein that binds to a highly conserved DNA motif and facilitates regulation of three-dimensional (3D) nuclear architecture and transcription. CTCF binding sites (CTCF-BSs) reside in non-coding DNA and are frequently mutated in cancer. Our previous study identified a small subclass of CTCF-BSs that are resistant to CTCF knock down, termed persistent CTCF binding sites (P-CTCF-BSs). P-CTCF-BSs show high binding conservation and potentially regulate cell-type constitutive 3D chromatin architecture. Here, using ICGC sequencing data we made the striking observation that P-CTCF-BSs display a highly elevated mutation rate in breast and prostate cancer when compared to all CTCF-BSs. To address whether P-CTCF-BS mutations are also enriched in other cell-types, we developed CTCF-INSITE-a tool utilising machine learning to predict persistence based on genetic and epigenetic features of experimentally-determined P-CTCF-BSs. Notably, predicted P-CTCF-BSs also show a significantly elevated mutational burden in all 12 cancer-types tested. Enrichment was even stronger for P-CTCF-BS mutations with predicted functional impact to CTCF binding and chromatin looping. Using in vitro binding assays we validated that P-CTCF-BS cancer mutations, predicted to be disruptive, indeed reduced CTCF binding. Together this study reveals a new subclass of cancer specific CTCF-BS DNA mutations and provides insights into their importance in genome organization in a pan-cancer setting.

Diffusible signal factor (DSF)-mediated quorum sensing modulates swarming in Xanthomonas albilineans.

Xanthomonas spp. are plant pathogens known for significantly impacting crop yields. Among them, Xanthomonas albilineans (Xal) is notable for colonizing the xylem and causing sugarcane leaf scald disease. This study employed homologous recombination to mutate quorum sensing (QS) regulatory genes (rpf) to investigate their role in Xal pathogenicity. Deletions of rpfF (ΔrpfF), rpfC (ΔrpfC), and rpfG (ΔrpfG) led to reduced swarming, growth, and virulence. However, DSF supplementation restored swarming and growth in the ΔrpfF mutant. Deleting rpfC, rpfG, and rpfF also reduced twitching motility and affected Type IV Pilus (T4P) expression. Transcriptomic analysis revealed that ΔrpfF positively regulates flagellar genes. DSF supplementation in ΔrpfF (ΔrpfF-DSF) modulated the expression of flagellar, chemotaxis, and T4P genes. These findings elucidate the DSF-mediated swarming pathway in Xal and provide valuable insights into its regulatory mechanisms.

Interface Engineering in 2D/2D Heterogeneous Photocatalysts.

Assembling different 2D nanomaterials into heterostructures with strong interfacial interactions presents a promising approach for novel artificial photocatalytic materials. Chemically implementing the 2D nanomaterials' construction/stacking modes to regulate different interfaces can extend their functionalities and achieve good performance. Herein, based on different fundamental principles and photochemical processes, multiple construction modes (e.g., face-to-face, edge-to-face, interface-to-face, edge-to-edge) are overviewed systematically with emphasis on the relationships between their interfacial characteristics (e.g., point, linear, planar), synthetic strategies (e.g., in situ growth, ex situ assembly), and enhanced applications to achieve precise regulation. Meanwhile, recent efforts for enhancing photocatalytic performances of 2D/2D heterostructures are summarized from the critical factors of enhancing visible light absorption, accelerating charge transfer/separation, and introducing novel active sites. Notably, the crucial roles of surface defects, cocatalysts, and surface modification for photocatalytic performance optimization of 2D/2D heterostructures are also discussed based on the synergistic effect of optimization engineering and heterogeneous interfaces. Finally, perspectives and challenges are proposed to emphasize future opportunities for expanding 2D/2D heterostructures for photocatalysis.

Low dose PFDA induces DNA damage and DNA repair inhibition by promoting nuclear cGAS accumulation in ovarian epithelial cells.

Per- and polyfluoroalkyl substances (PFAS), the versatile anthropogenic chemicals, are popular with the markets and manufactured in large quantities yearly. Accumulation of PFAS has various adverse health effects on human. Albeit certain members of PFAS were found to have genotoxicity in previous studies, the mechanisms underlying their effects on DNA damage repair remain unclear. Here, we investigated the effects of Perfluorodecanoic acid (PFDA) on DNA damage and DNA damage repair in ovarian epithelial cells through a series of in vivo and in vitro experiments. At environmentally relevant concentration, we firstly found that PFDA can cause DNA damage in primary mouse ovarian epithelial cells and IOSE-80 cells. Moreover, nuclear cGAS increased in PFDA-treated cells, which leaded to the efficiency of DNA homologous recombination (HR) decreased and DNA double-strand breaks perpetuated. In vivo experiments also verified that PFDA can induce more DNA double-strand breaks lesions and nuclear cGAS in ovarian tissue. Taken together, our results unveiled that low dose PFDA can cause deleterious effects on DNA and DNA damage repair (DDR) in ovarian epithelial cells and induce genomic instability.

Identification of a ZC3H12D-regulated competing endogenous RNA network for prognosis of lung adenocarcinoma at single-cell level.

BACKGROUND: To identify hub genes from the competing endogenous RNA (ceRNA) network of lung adenocarcinoma (LUAD) and to explore their potential functions on prognosis of patients from a single-cell perspective. METHODS: We performed RNA-sequencing of LUAD to construct ceRNA regulatory network, integrating with public databases to identify the vital pathways related to patients' prognosis and to reveal the expression level of hub genes under different conditions, the functional enrichment of co-expressed genes and their potential immune-related mechanisms. RESULTS: ZC3H12D-hsa-miR-4443-ENST00000630242 axis was found to be related with LUAD. Lower ZC3H12D expression was significantly associated with shorter overall survival (OS) of patients (HR = 2.007, P < 0.05), and its expression was higher in early-stage patients, including T1 (P < 0.05) and N0 (P < 0.05). Additionally, ZC3H12D expression was higher in immune cells displayed by single-cell RNA-sequencing data, especially in Treg cells of lung cancer and CD8 T cells, B cells and CD4 T cells of LUAD. The functional enrichment analysis showed that the co-expressed genes mainly played a role in lymphocyte activation and cytokine-cytokine receptor interaction. In addition, ZC3H12D was associated with multiple immune cells and immune molecules, including immune checkpoints CTLA4, CD96 and TIGIT. CONCLUSION: ZC3H12D-hsa-miR-4443-ENST00000630242 ceRNA network was identified in LUAD. ZC3H12D could affect prognosis of patients by regulating mRNA, miRNA, lncRNA, immune cells and immune molecules. Therefore, it may serve as a vital predictive marker and could be regarded as a potential therapeutic target for LUAD in the future.

Whole genome deep sequencing analysis of cell-free DNA in samples with low tumour content.

BACKGROUND: Circulating cell-free DNA (cfDNA) in the plasma of cancer patients contains cell-free tumour DNA (ctDNA) derived from tumour cells and it has been widely recognized as a non-invasive source of tumour DNA for diagnosis and prognosis of cancer. Molecular profiling of ctDNA is often performed using targeted sequencing or low-coverage whole genome sequencing (WGS) to identify tumour specific somatic mutations or somatic copy number aberrations (sCNAs). However, these approaches cannot efficiently detect all tumour-derived genomic changes in ctDNA. METHODS: We performed WGS analysis of cfDNA from 4 breast cancer patients and 2 patients with benign tumours. We sequenced matched germline DNA for all 6 patients and tumour samples from the breast cancer patients. All samples were sequenced on Illumina HiSeqXTen sequencing platform and achieved approximately 30x, 60x and 100x coverage on germline, tumour and plasma DNA samples, respectively. RESULTS: The mutational burden of the plasma samples (1.44 somatic mutations/Mb of genome) was higher than the matched tumour samples. However, 90% of high confidence somatic cfDNA variants were not detected in matched tumour samples and were found to comprise two background plasma mutational signatures. In contrast, cfDNA from the di-nucleosome fraction (300 bp-350 bp) had much higher proportion (30%) of variants shared with tumour. Despite high coverage sequencing we were unable to detect sCNAs in plasma samples. CONCLUSIONS: Deep sequencing analysis of plasma samples revealed higher fraction of unique somatic mutations in plasma samples, which were not detected in matched tumour samples. Sequencing of di-nucleosome bound cfDNA fragments may increase recovery of tumour mutations from plasma.

C1QTNF6 regulated by miR-29a-3p promotes proliferation and migration in stage I lung adenocarcinoma.

OBJECTIVE: C1QTNF6 has been implicated as an essential component in multiple cellular and molecular preliminary event, including inflammation, glucose metabolism, endothelial cell modulation and carcinogenesis. However, the biological process and potential mechanism of C1QTNF6 in lung adenocarcinoma (LUAD) are indefinite and remain to be elucidated. Therefore, we investigated the interaction among the traits of C1QTNF6 and LUAD pathologic process. METHODS: RT-qPCR and western blot were conducted to determine the expression levels of C1QTNF6. RNA interference and overexpression of C1QTNF6 were constructed to identify the biological function of C1QTNF6 in cellular proliferative, migratory and invasive potentials in vitro. Dual-luciferase reporter assay was applied to identify the possible interaction between C1QTNF6 and miR-29a-3p. Moreover, RNA sequencing analysis of C1QTNF6 knockdown was performed to identify the potential regulatory pathways. RESULTS: C1QTNF6 was upregulated in stage I LUAD tissues compared with adjacent non-cancerous tissues. Concurrently, C1QTNF6 knockdown could remarkably inhibit cell proliferation, migratory and invasive abilities, while overexpression of C1QTNF6 presented opposite results. Additionally, miR-29a-3p may serve as an upstream regulator of C1QTNF6 and reduce the expression of C1QTNF6. Subsequent experiments showed that miR-29a-3p could decrease the cell mobility and proliferation positive cell rates, as well as reduce the migratory and invasive possibilities in LUAD cells via downregulating C1QTNF6. Moreover, RNA sequencing analysis demonstrated that the cytokine-cytokine receptor interaction pathway may participate in the process of C1QTNF6 regulating tumor progression. CONCLUSION: Our study first demonstrated that downregulation of C1QTNF6 could inhibit tumorigenesis and progression in LUAD cells negatively regulated by miR-29a-3p. These consequences could reinforce our awareness and understanding of the underlying mechanism and provide a promising therapeutic target for LUAD.

Bent crystal Laue analyser combined with total reflection fluorescence X-ray absorption fine structure (BCLA + TRF-XAFS) and its application to surface studies.

A bent crystal Laue analyser (BCLA) is an X-ray energy analyser used for fluorescence X-ray absorption fine-structure (XAFS) spectroscopy to separate the fluorescence X-ray emission line of a target atom from the elastic scattering X-rays and other fluorescence emission lines. Here, the feasibility of the BCLA for total reflection fluorescence XAFS (TRF-XAFS), which has a long X-ray footprint on the substrate surface owing to grazing incidence, was tested. The focal line of the BCLA was adjusted on the X-ray footprint and the XAFS signal for one monolayer of Pt deposited on a 60 nm Au film with high sensitivity was obtained. Although range-extended XAFS was expected by the rejection of Au fluorescence arising from the Au substrate, a small glitch was found in the Au L3 edge because of the sudden change of the complex refraction index of the Au substrate at the Au edge. This abnormal spectrum feature can be removed by reflectivity correction using Au foil absorption data. BCLA combined with TRF-XAFS spectroscopy (BCLA + TRF-XAFS) is a new technique for the in situ surface analysis of highly dispersed systems even in the presence of a liquid overlayer.

sCNAphase: using haplotype resolved read depth to genotype somatic copy number alterations from low cellularity aneuploid tumors.

Accurate identification of copy number alterations is an essential step in understanding the events driving tumor progression. While a variety of algorithms have been developed to use high-throughput sequencing data to profile copy number changes, no tool is able to reliably characterize ploidy and genotype absolute copy number from tumor samples that contain less than 40% tumor cells. To increase our power to resolve the copy number profile from low-cellularity tumor samples, we developed a novel approach that pre-phases heterozygote germline single nucleotide polymorphisms (SNPs) in order to replace the commonly used 'B-allele frequency' with a more powerful 'parental-haplotype frequency'. We apply our tool-sCNAphase-to characterize the copy number and loss-of-heterozygosity profiles of four publicly available breast cancer cell-lines. Comparisons to previous spectral karyotyping and microarray studies revealed that sCNAphase reliably identified overall ploidy as well as the individual copy number mutations from each cell-line. Analysis of artificial cell-line mixtures demonstrated the capacity of this method to determine the level of tumor cellularity, consistently identify sCNAs and characterize ploidy in samples with as little as 10% tumor cells. This novel methodology has the potential to bring sCNA profiling to low-cellularity tumors, a form of cancer unable to be accurately studied by current methods.

GtTR: Bayesian estimation of absolute tandem repeat copy number using sequence capture and high throughput sequencing.

BACKGROUND: Tandem repeats comprise significant proportion of the human genome including coding and regulatory regions. They are highly prone to repeat number variation and nucleotide mutation due to their repetitive and unstable nature, making them a major source of genomic variation between individuals. Despite recent advances in high throughput sequencing, analysis of tandem repeats in the context of complex diseases is still hindered by technical limitations. We report a novel targeted sequencing approach, which allows simultaneous analysis of hundreds of repeats. We developed a Bayesian algorithm, namely - GtTR - which combines information from a reference long-read dataset with a short read counting approach to genotype tandem repeats at population scale. PCR sizing analysis was used for validation. RESULTS: We used a PacBio long-read sequenced sample to generate a reference tandem repeat genotype dataset with on average 13% absolute deviation from PCR sizing results. Using this reference dataset GtTR generated estimates of VNTR copy number with accuracy within 95% high posterior density (HPD) intervals of 68 and 83% for capture sequence data and 200X WGS data respectively, improving to 87 and 94% with use of a PCR reference. We show that the genotype resolution increases as a function of depth, such that the median 95% HPD interval lies within 25, 14, 12 and 8% of the its midpoint copy number value for 30X, 200X WGS, 395X and 800X capture sequence data respectively. We validated nine targets by PCR sizing analysis and genotype estimates from sequencing results correlated well with PCR results. CONCLUSIONS: The novel genotyping approach described here presents a new cost-effective method to explore previously unrecognized class of repeat variation in GWAS studies of complex diseases at the population level. Further improvements in accuracy can be obtained by improving accuracy of the reference dataset.

Enhancing Upconversion from NaYF4:Yb,Er@NaYF4 Core-Shell Nanoparticles Assembled on Metallic Nanostructures.

We report a simple method for the fabrication of a three-layered plasmonic structure of silicon substrate-Au nanospheres-upconversion particles (UCNPs) that displays up to 101-fold fluorescence enhancement. Monodispersed pure hexagonal-phase NaYF4:Yb,Er core and NaYF4:Yb,Er@NaYF4 core@shell nanocrystals were prepared by a solvothermal method. Two dimensional (2D) assembled Au spheres were prepared on a Si substrate, and then, 2D arrays of UCNPs were deposited on the grown 2D monolayered Au spheres by a self-organizing process. The distance between plasmonic Au NPs and rare-earth (RE) core was finely adjusted by changing the undoped NaYF4 shell thickness. The UC emission enhancement shows a pronounced shell thickness dependence. For the non-plasmonic structured samples, a single peak in upconversion luminescence (UCL) enhancement was observed as the undoped NaYF4 shell thickness increases from 0 nm to 23.0 nm. In contrast, for the plasmonic structured samples, multi-oscillations in UCL enhancement were observed in the undoped NaYF4 shell thickness range of 0-23.0 nm, where the UCL enhancement factors of three bands (521 nm, 540 nm and 654 nm) are high up to 65, 101 and 61, respectively, at 19.6 nm-thick NaYF4 shell. The multi-oscillations in UCL enhancement in the plasmonic samples can be associated with plasmonic coupling between arrays of core-shell UCNPs with various sizes and the underlying 2D Au spheres. The related mechanisms of the UCL enhancements are discussed.

VirB11, a traffic ATPase, mediated flagella assembly and type IV pilus morphogenesis to control the motility and virulence of Xanthomonas albilineans.

Xanthomonas albilineans (Xal) is a gram-negative bacterial pathogen responsible for developing sugarcane leaf scald disease, which engenders significant economic losses within the sugarcane industry. In the current study, homologous recombination exchange was carried out to induce mutations within the virB/D4-like type IV secretion system (T4SS) genes of Xal. The results revealed that the virB11-deletion mutant (ΔvirB11) exhibited a loss in swimming and twitching motility. Application of transmission electron microscopy analysis further demonstrated that the ΔvirB11 failed to develop flagella formation and type IV pilus morphology and exhibited reduced swarming behaviour and virulence. However, these alterations had no discernible impact on bacterial growth. Comparative transcriptome analysis between the wild-type Xal JG43 and the deletion-mutant ΔvirB11 revealed 123 differentially expressed genes (DEGs), of which 28 and 10 DEGs were notably associated with flagellar assembly and chemotaxis, respectively. In light of these findings, we postulate that virB11 plays an indispensable role in regulating the processes related to motility and chemotaxis in Xal.

Anaerobic fermentation for hydrogen production and tetracycline degradation: Biodegradation mechanism and microbial community succession.

The misuse and continues discharge of antibiotics can cause serious pollution, which is urgent to take steps to remit the environment pollution. In this study, anaerobic bacteria isolated from the aeration tank of a local sewage treatment plant were employed to investigate hydrogen production and tetracycline (TC) degradation during anaerobic fermentation. Results indicate that low concentrations of TC enhanced hydrogen production, increasing from 366 mL to a maximum of 480 mL. This increase is attributed to stimulated hydrolysis and acidogenesis, coupled with significant inhibition of homoacetogenesis. Furthermore, the removal of TC, facilitated by adsorption and biodegradation, exceeded 90 %. During the fermentation process, twenty-one by-products were identified, leading to the proposal of four potential degradation pathways. Analysis of the microbial community revealed shifts in diversity and a decrease in the abundance of hydrogen-producing bacteria, whereas bacteria harboring tetracycline resistance genes became more prevalent. This study provides a possibility to treat tetracycline-contaminated wastewater and to produce clean energy simultaneously by anaerobic fermentation.

Simultaneous CO2 and biogas slurry treatment using a newly isolated microalga with high CO2 tolerance.

The fixation of carbon dioxide (CO2) using microalgae is a promising CO2 capture and utilization technology. Microalgae have also been suggested as a potential treatment for biogas slurry (BS). This study screened microalgae capable of tolerating both high CO2 concentrations and BS, assessed their CO2 fixation and pollutant removal capabilities, and evaluated the potential use of the resulting algal biomass. Chlamydopodium sp. HS01, which showed the highest tolerance to 15% CO2 and BS, was selected due to its strong growth, CO2 fixation, and ammonia nitrogen removal abilities. The generated biomass also demonstrated significant potential for bioenergy production. Metabolomics analysis revealed that the lipid composition of HS01 underwent substantial changes under 15% CO2 alone and in combination with BS, likely as a stress adaptation strategy. Overall, HS01 presents high potential for resource utilization of CO2 coupled with actual BS.

Overcoming tyrosine kinase inhibitor resistance in lung cancer brain metastasis with CTLA4 blockade.

Lung cancer brain metastasis (LCBM) poses a significant clinical challenge due to acquired resistance to tyrosine kinase inhibitor (TKI) treatment. To elucidate its underlying mechanisms, we employed single-cell RNA sequencing analysis on surgically obtained LCBM samples with diverse genetic backgrounds and TKI treatment histories. Our study uncovers that TKI treatment elevates the immune checkpoint CTLA4 expression in T cells, promoting an immune-suppressive microenvironment. This immunomodulation is initiated by tumor-derived HMGB1 in response to TKIs. In LCBM syngeneic murine models with TKI-sensitive or TKI-resistant EGFR mutations, combining CTLA4 blockade with TKIs demonstrates enhanced efficacy over TKI monotherapy or TKIs with PD1 blockade. These findings provide insights into the TKI resistance mechanisms and highlight the potential of CTLA4 blockade in effectively overcoming TKI resistance in LCBM.

The potential of epigenetic therapy to target the 3D epigenome in endocrine-resistant breast cancer.

Three-dimensional (3D) epigenome remodeling is an important mechanism of gene deregulation in cancer. However, its potential as a target to counteract therapy resistance remains largely unaddressed. Here, we show that epigenetic therapy with decitabine (5-Aza-mC) suppresses tumor growth in xenograft models of pre-clinical metastatic estrogen receptor positive (ER+) breast tumor. Decitabine-induced genome-wide DNA hypomethylation results in large-scale 3D epigenome deregulation, including de-compaction of higher-order chromatin structure and loss of boundary insulation of topologically associated domains. Significant DNA hypomethylation associates with ectopic activation of ER-enhancers, gain in ER binding, creation of new 3D enhancer-promoter interactions and concordant up-regulation of ER-mediated transcription pathways. Importantly, long-term withdrawal of epigenetic therapy partially restores methylation at ER-enhancer elements, resulting in a loss of ectopic 3D enhancer-promoter interactions and associated gene repression. Our study illustrates the potential of epigenetic therapy to target ER+ endocrine-resistant breast cancer by DNA methylation-dependent rewiring of 3D chromatin interactions, which are associated with the suppression of tumor growth.

E-Spin: A Stochastic Ising Spin Based on Electrically-Controlled MTJ for Constructing Large-Scale Ising Annealing Systems.

With its unique computer paradigm, the Ising annealing machine has become an emerging research direction. The Ising annealing system is highly effective at addressing combinatorial optimization (CO) problems that are difficult for conventional computers to tackle. However, Ising spins, which comprise the Ising system, are difficult to implement in high-performance physical circuits. We propose a novel type of Ising spin based on an electrically-controlled magnetic tunnel junction (MTJ). Electrical operation imparts true randomness, great stability, precise control, compact size, and easy integration to the MTJ-based spin. In addition, simulations demonstrate that the frequency of electrically-controlled stochastic Ising spin (E-spin) is 50 times that of the thermal disturbance MTJ-based spin (p-bit). To develop a large-scale Ising annealing system, up to 64 E-spins are implemented. Our Ising annealing system demonstrates factorization of integers up to 264 with a temporal complexity of around O(n). The proposed E-spin shows superiority in constructing large-scale Ising annealing systems and solving CO problems.

Identification of key pathways and genes in nasopharyngeal carcinoma based on WGCNA.

OBJECTIVE: We aim to identify the potential genes and signaling pathways associated with the nasopharyngeal carcinoma (NPC) prognosis using Weighted Gene Co-Expression Network Analysis (WGCNA). METHODS: Gene Expression Omnibus (GEO) query was utilized to download two NPC mRNA microarray data. WGCNA was conducted on differentially expressed genes (DEGs) to obtain tumor-associated gene modules. Genes in core modules were intersected with DEGs for gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analysis. GSE102349 dataset was devoted to identifying prognostic hub genes by survival analysis and the results were confirmed by quantitative polymerase chain reaction (qPCR). RESULTS: Co-expression networks were built, and we detected 12 gene modules. The Brown module and Magenta module were extremely associated with NPC samples. GO functional analysis and KEGG pathway analysis was carried out to the genes in the Brown and Magenta modules. Our data indicated that DEGs in Brown module and Magenta module were correlated with the biological regulation, metabolic process, reproduction, and cellular proliferation. Twenty-six hub genes were obtained and were considered to be closely related to NPC. GSE102349 dataset was devoted to identifying prognostic hub genes by survival analysis. The expression of IL33, MPP3 and SLC16A7 in GSE102349 dataset was significantly correlated with the progression-free survival (PFS). The results of qPCR indicated a strong correlation between SLC16A7 expression and the overall survival (OS). CONCLUSIONS: WGCNA contributed to the detection of gene modules and identification of hub genes and crucial genes. These crucial genes might be potential targets for pharmaceutic therapies with potential clinical significance.

Multiomics Integrative Analysis Identifying EPC1 as a Prognostic Biomarker in Head and Neck Squamous Cell Carcinoma.

Background: Biomarker research in head and neck squamous cell carcinoma (HNSCC) is constantly revealing promising findings. An enhancer of polycomb homolog 1 (EPC1) was found to play a procancer role in nasopharyngeal carcinoma (NPC), but its role in HNSCC with strong heterogeneity is still unclear. Herein, we investigated the prognostic significance and related mechanisms of EPC1 in HNSCC. Methods: The Kaplan-Meier plotter was used to evaluate the prognostic significance of EPC1. Based on a range of published public databases, the multiomics expression of EPC1 in HNSCC was explored to investigate the mechanisms affecting prognosis. Results: According to the clinical data, high EPC1 expression in HNSCC was a predictor of patient prognosis (hazard ratio (HR) = 0.64; 95% confidence interval (CI) 0.49-0.83; P < 0.01). EPC1 expression varied among clinical subtypes and was related to key factors, such as TP53 and human papillomavirus (HPV) (P < 0.05). At the genetic level, EPC1 expression level may be associated with protein phosphorylation, cell adhesion, cancer-related pathways, etc. For the noncoding region, a competing endogenous RNA network was constructed, and 6 microRNAs and 12 long noncoding RNAs were identified. At the protein level, a protein-protein interaction (PPI) network related to EPC1 expression was constructed and found to be involved in HPV infection, endocrine resistance, and multiple cancer pathways. At the immune level, EPC1 expression was correlated with a variety of immune cells and immune molecules, which together constituted the immune microenvironments of tumors. Conclusion: High EPC1 expression may predict a better prognosis in HNSCC, as it is more frequently found in HNSCC with HPV infection. EPC1 may participate in the genomics, transcriptomics, proteomics, and immunomics of HNSCC, and the results can provide a reference for the development of targeted drugs and evaluation of patient prognosis.