The harmonizome: a collection of processed datasets gathered to serve and mine knowledge about genes and proteins

Crosstalk between degradation and bioenergetics: how autophagy and endolysosomal processes regulate energy production

Cells undergo metabolic reprogramming to adapt to changes in nutrient availability, cellular activity, and transitions in cell states. The balance between glycolysis and mitochondrial respiration is crucial for energy production, and metabolic reprogramming stipulates a shift in such balance to optimize both bioenergetic efficiency and anabolic requirements. Failure in switching bioenergetic dependence can lead to maladaptation and pathogenesis. While cellular degradation is known to recycle precursor molecules for anabolism, its potential role in regulating energy production remains less explored. The bioenergetic switch between glycolysis and mitochondrial respiration involves transcription factors and organelle homeostasis, which are both regulated by the cellular degradation pathways. A growing body of studies has demonstrated that both stem cells and differentiated cells exhibit bioenergetic switch upon perturbations of autophagic activity or endolysosomal processes. Here, we highlighted the current understanding of the interplay between degradation processes, specifically autophagy and endolysosomes, transcription factors, endolysosomal signaling, and mitochondrial homeostasis in shaping cellular bioenergetics. This review aims to summarize the relationship between degradation processes and bioenergetics, providing a foundation for future research to unveil deeper mechanistic insights into bioenergetic regulation.

ANXA5 predicts prognosis and immune response and mediates proliferation and migration in head and neck squamous cell carcinoma

BACKGROUND

Head and neck squamous cell carcinoma (HNSCC) is a common malignancy that often develops unnoticed. Typically, these tumors are identified at advanced stages, resulting in a relatively low chance of successful treatment. Anoikis serves as a natural defense against the spread of tumor cells, meaning circumventing anoikis can effectively inhibit tumor metastasis. Nonetheless, studies focusing on anoikis in the context of HNSCC remain scarce.

METHODS

Anoikis-related genes (ARGs) were identified by using the GeneCards and Harmonizome databases. Expression data of these genes and relevant clinical features were downloaded from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. A LASSO regression and a prognostic risk score model were developed to determine their prognostic significance. The analysis included the use of the CIBERSORT algorithm to quantify immune and stromal cell presence. Furthermore, in vitro and in vivo, we confirmed the expression and functional roles of proteins and mRNA of genes independently predictive of prognosis.

RESULTS

The study identified eight genes linked to prognosis (ANXA5, BAK1, CDKN2A, PPARG, CCR7, MAPK11, CRYAB, CRYBA1) and developed a prognostic model that effectively forecasts the survival outcomes for patients with HNSCC. A higher survival likelihood is associated with lower risk scores. In addition, a significant relationship was found between immune and risk score, and ANXA5 deletion promoted the killing of HNSCC cells by activated CD8+ T cells. During the screening process, 65 different chemotherapeutic drugs were found to have significant differences in IC50 values when comparing high- and low-risk categories. ANXA5 emerged as a gene with independent prognostic significance, exhibiting notably elevated protein and mRNA levels in HNSCC tissue compared to non-tumorous tissue. The suppression of ANXA5 gene activity resulted in a substantial decrease in both the growth and mobility of HNSCC cells. Animal model experiments demonstrated that inhibiting ANXA5 suppressed HNSCC growth and migration in vivo.

CONCLUSION

Through bioinformatics, a prognostic risk model of high precision was developed, offering valuable insights into the survival rates and immune responses in patients with HNSCC. ANXA5 is highlighted as a significant prognostic factor among the identified genes, indicating its promise as a potential therapeutic target for those with HNSCC.

Molecular classification of geriatric breast cancer displays distinct senescent subgroups of prognostic significance

Breast cancer in the elderly presents distinct biological characteristics and clinical treatment responses compared with cancer in younger patients. Comprehensive Geriatric Assessment is recommended for evaluating treatment efficacy in elderly cancer patients based on physiological classification. However, research on molecular classification in older cancer patients remains insufficient. In this study, we identified two subgroups with distinct senescent clusters among geriatric breast cancer patients through multi-omics analysis. Using various machine learning algorithms, we developed a comprehensive scoring model called "Sene_Signature," which more accurately distinguished elderly breast cancer patients compared with existing methods and better predicted their prognosis. The Sene_Signature was correlated with tumor immune cell infiltration, as supported by single-cell transcriptomics, RNA sequencing, and pathological data. Furthermore, we observed increased drug responsiveness in patients with a high Sene_Signature to treatments targeting the epidermal growth factor receptor and cell-cycle pathways. We also established a user-friendly web platform to assist investigators in assessing Sene_Signature scores and predicting treatment responses for elderly breast cancer patients. In conclusion, we developed a novel model for evaluating prognosis and therapeutic responses, providing a potential molecular classification that assists in the pre-treatment assessment of geriatric breast cancer.

NADPH oxidases in healthy white adipose tissue and in obesity: function, regulation, and clinical implications

Reactive oxygen species, when produced in a controlled manner, are physiological modulators of healthy white adipose tissue (WAT) expansion and metabolic function. By contrast, unbridled production of oxidants is associated with pathological WAT expansion and the establishment of metabolic dysfunctions, most notably insulin resistance and type 2 diabetes mellitus. NADPH oxidases (NOXs) produce oxidants in an orderly fashion and are present in adipocytes and in other diverse WAT-constituent cell types. Recent studies have established several links between aberrant NOX-derived oxidant production, adiposity, and metabolic homeostasis. The objective of this review is to highlight the physiological roles attributed to diverse NOX isoforms in healthy WAT and summarize current knowledge of the metabolic consequences related to perturbations in their adequate oxidant production. We detail WAT-related alterations in preclinical investigations conducted in NOX-deficient murine models. In addition, we review clinical studies that have employed NOX inhibitors and currently available data related to human NOX mutations in metabolic disturbances. Future investigations aimed at understanding the integration of NOX-derived oxidants in the regulation of the WAT cellular redox network are essential for designing successful redox-related precision therapies to curb obesity and attenuate obesity-associated metabolic pathologies.

Actin filament dynamics at barbed ends: New structures, new insights

The dynamic actin cytoskeleton contributes to many critical biological processes by providing the structural support underlying the morphology of most cells, facilitating intracellular transport, and generating forces required for cell motility and division. To execute many of these functions, actin monomers polymerize into polarized filaments that display different structural and biochemical properties at each end. Filament dynamics are regulated by diverse regulatory proteins which collaborate to dictate rates of elongation and disassembly, particularly at the fast-growing barbed (plus) end. This review highlights the biochemical mechanisms of six barbed end regulatory proteins: formin, profilin, capping protein, IQGAP1, cyclase-associated protein, and twinfilin. We discuss how individual proteins influence actin dynamics and how several intriguing complex assemblies influence the polymerization fate of actin filaments. Understanding these mechanisms offers insights into how actin is regulated in essential cell processes and dysregulated in disease.

Characterization of the AGR2-NPM3 axis uncovers the AGR2 involvement in PD-L1 regulation in colorectal cancer

Despite extensive research, the molecular role of AGR2 in the progression and metastasis of colorectal cancer (CRC) has not been fully characterized. We used quantitative mass spectrometry (SWATH MS) to identify differentially expressed proteins in paired CRC cell models of the SW480 and SW620 cell lines in response to AGR2 protein level manipulation. Relying on the results from SWATH MS and subsequent immunochemical validation, we selected NMP3 as the top candidate protein associated with AGR2 in CRC tumour cells in our screen. RT‒qPCR and immunochemical analysis confirmed the involvement of AGR2-mediated regulation of NPM3 at the transcriptional and posttranscriptional levels. Since PD-L1 is a constituent of the NPM3 regulatory axis, we aimed to correlate the changes in PD-L1 to the differential expression of AGR2 in our cell models. We found that AGR2 positively regulates PD-L1 levels in both SW480 and SW620 cell lines; additionally, several different CRC patient transcriptome cohorts confirmed the association of AGR2 with PD-L1. Our work reveals a new AGR2-NPM3 regulatory axis and the involvement of AGR2 in the regulation of PD-L1, which paves the way for the association of AGR2 with immune evasion in CRC cells.

Non-Negative Matrix Tri-Factorization for Representation Learning in Multi-Omics Datasets with Applications to Drug Repurposing and Selection

The vast corpus of heterogeneous biomedical data stored in databases, ontologies, and terminologies presents a unique opportunity for drug design. Integrating and fusing these sources is essential to develop data representations that can be analyzed using artificial intelligence methods to generate novel drug candidates or hypotheses. Here, we propose Non-Negative Matrix Tri-Factorization as an invaluable tool for integrating and fusing data, as well as for representation learning. Additionally, we demonstrate how representations learned by Non-Negative Matrix Tri-Factorization can effectively be utilized by traditional artificial intelligence methods. While this approach is domain-agnostic and applicable to any field with vast amounts of structured and semi-structured data, we apply it specifically to computational pharmacology and drug repurposing. This field is poised to benefit significantly from artificial intelligence, particularly in personalized medicine. We conducted extensive experiments to evaluate the performance of the proposed method, yielding exciting results, particularly compared to traditional methods. Novel drug-target predictions have also been validated in the literature, further confirming their validity. Additionally, we tested our method to predict drug synergism, where constructing a classical matrix dataset is challenging. The method demonstrated great flexibility, suggesting its applicability to a wide range of tasks in drug design and discovery.

HLA-E and NKG2A Mediate Resistance to M. bovis BCG Immunotherapy in Non-Muscle-Invasive Bladder Cancer

Mycobacterium bovis Bacillus Calmette-Guerin (BCG) is the primary treatment for non-muscle-invasive bladder cancer (NMIBC), known to stimulate inflammatory cytokines, notably interferon (IFN)-γ. We observed that prolonged IFN-γ exposure fosters adaptive resistance in recurrent tumors, aiding immune evasion and tumor proliferation. We identify HLA-E and NKG2A, part of a novel NK and T cell checkpoint pathway, as key mediators of resistance in BCG-unresponsive NMIBC. IFN-γ enhances HLA-E and PD-L1 expression in recurrent tumors, with an enrichment of intra-tumoral NKG2A-expressing NK and CD8 T cells. CXCL9+ macrophages and dendritic cells and CXCL12-expressing stromal cells likely recruit CXCR3/CXCR4-expressing NK and T cells and CXCR7+ HLA-EHIGH tumor cells. NK and CD8 T cells remain functional within BCG-unresponsive tumors but are inhibited by HLA-E and PD-L1, providing a framework for combined NKG2A and PD-L1 blockade strategy for bladder-sparing treatment of BCG-unresponsive NMIBC.

Design of linked-domain protein inhibitors of UBE2D as tools to study cellular ubiquitination

Ubiquitin (Ub) is a post-translational modification that largely controls proteostasis through mechanisms spanning transcription, translation, and notably, protein degradation. Ub conjugation occurs through a hierarchical cascade of three enzyme classes (E1, E2, and E3s) involving >1000 proteins that regulate the ubiquitination of proteins. The E2 Ub-conjugating enzymes are the midpoint, yet their cellular roles remain under-characterized, partly due to a lack of inhibitors. For example, the cellular roles of the promiscuous E2 UBE2D/UBCH5 are not well described. Here, we develop a highly selective, multivalent, engineered protein inhibitor for the UBE2D family that simultaneously targets the RING- and backside-binding sites. In HeLa cells, these inhibitors phenocopy knockdown of UBE2D by reducing the IC 50 to cisplatin and whole-cell proteomics reveal an increased abundance of ∼20% of the identified proteins, consistent with reduced Ub degradation and proteotoxic stress. These precision tools will enable new studies probing UBE2D's central role in proteome management.

A cochlear progenitor pool influences patterning of the mammalian sensory epithelium via MYBL2

During embryonic development, Wnt signaling influences both proliferation and sensory formation in the cochlea. How this dual nature of Wnt signaling is coordinated is unknown. In this study, we define a novel role for a Wnt-regulated gene, Mybl2, which was already known to be important for proliferation, in determining the size and patterning of the sensory epithelium in the murine cochlea. Using a quantitative spatial analysis approach and analyzing Mybl2 loss-of-function, we show that Mybl2 promoted proliferation in the inner sulcus domain but limited the size of the sensory domain by influencing their adjoining boundary position via Jag1 regulation during development. Mybl2 loss-of-function simultaneously decreased proliferation in the inner sulcus and increased the size of the sensory domain, resulting in a wider sensory epithelium with ectopic inner hair cell formation during late embryonic stages. These data suggest that progenitor cells in the inner sulcus determine boundary formation and pattern the sensory epithelium via MYBL2.

Discovery of genomic loci for liver health and steatosis reveals overlap with glutathione redox genetics

Non-alcoholic fatty liver disease (NAFLD) is the most common chronic liver condition in the United States, encompassing a wide spectrum of liver pathologies including steatosis, steatohepatitis, fibrosis, and cirrhosis. Despite its high prevalence, there are no medications currently approved by the Food and Drug Administration for the treatment of NAFLD. Recent work has suggested that NAFLD has a strong genetic component and identifying causative genes will improve our understanding of the molecular mechanisms contributing to NAFLD and yield targets for future therapeutic investigations. Oxidative stress is known to play an important role in NAFLD pathogenesis, yet the underlying mechanisms accounting for disturbances in redox status are not entirely understood. To better understand the relationship between the glutathione redox system and signs of NAFLD in a genetically-diverse population, we measured liver weight, serum biomarkers aspartate aminotransferase (AST) and alanine aminotransferase (ALT), and graded liver pathology in a large cohort of Diversity Outbred mice. We compared hepatic endpoints to those of the glutathione redox system previously measured in the livers and kidneys of the same mice, and we screened for statistical and genetic associations using the R/qtl2 software. We discovered several novel genetic loci associated with markers of liver health, including loci that were associated with both liver steatosis and glutathione redox status. Candidate genes within each locus point to possible new mechanisms underlying the complex relationship between NAFLD and the glutathione redox system, which could have translational implications for future studies targeting NAFLD pathology.

Comparative analysis of single-cell pathway scoring methods and a novel approach

Single-cell gene set analysis (scGSA) provides a useful approach for quantifying molecular functions and pathways in high-throughput transcriptomic data, facilitating the biological interpretation of complex human datasets. However, various factors such as gene set size, quality of the gene sets and the dropouts impact the performance of scGSA. To address these limitations, we present a single-cell Pathway Score (scPS) method to measure gene set activity at single-cell resolution. Furthermore, we benchmark our method with six other methods: AUCell, AddModuleScore, JASMINE, UCell, SCSE and ssGSEA. The comparison across all the methods using two different simulation approaches highlights the effect of cell count, gene set size, noise, condition-specific genes and zero imputation on their performance. The results of our study indicate that the scPS is comparable with other single-cell scoring methods and detects fewer false positives. Importantly, this work reveals critical variables in the scGSA.

Single-cell transcriptomic analysis of immune cell dynamics in the healthy human endometrium

The microenvironment of the endometrial immune system is crucial to the success of placental implantation and healthy pregnancy. However, the functionalities of immune cells across various stages of the reproductive cycle have yet to be fully comprehended. To address this, we conducted advanced bioinformatic analysis on 230,049 high-quality single-cell transcriptomes from healthy endometrial samples obtained during the proliferative, secretory, early pregnancy, and late pregnancy stages. Our investigation has unveiled that proliferative natural killer (NK) cells, a potential source of endometrial NK cells, exhibit the most robust proliferative and differentiation potential during non-pregnant stages. We have also identified similar differentiation trajectories of NK cells originating from proliferative NK cells across four stages. Notably, during early pregnancy, NK cells demonstrate the highest oxidative phosphorylation metabolism activity, and, in conjunction with macrophages and T cells, exhibit the strongest type II interferon response. With spatial transcriptome data, we have discerned that the most robust immune-non-immune interactions are associated with the promotion and inhibition of cell proliferation, differentiation and migration during four stages. Furthermore, we have compiled lists of stage-specific risk genes implicated in reproductive diseases, which hold promise as potential disease biomarkers. Our study provides insights into the dynamics of the endometrial immune microenvironment during different reproductive cycle stages, thus serving as a reference for detecting pathological changes during pregnancy.

Dimerization of hub protein DYNLL1 and bZIP transcription factor CREB3L1 enhances transcriptional activation of CREB3L1 target genes like arginine vasopressin

bZIP transcription factors can function as homodimers or heterodimers through interactions with their disordered coiled-coil domain. Such dimer assemblies are known to influence DNA-binding specificity and/or the recruitment of binding partners, which can cause a functional switch of a transcription factor from being an activator to a repressor. We recently identified the genomic targets of a bZIP transcription factor called CREB3L1 in rat hypothalamic supraoptic nucleus by ChIP-seq. The objective of this study was to investigate the CREB3L1 protein-to-protein interactome of which little is known. For this approach, we created and screened a rat supraoptic nucleus yeast two-hybrid prey library with the bZIP region of rat CREB3L1 as the bait. Our yeast two-hybrid approach captured five putative CREB3L1 interacting prey proteins in the supraoptic nucleus. One interactor was selected by bioinformatic analyses for more detailed investigation by co-immunoprecipitation, immunofluorescent cellular localisation, and reporter assays in vitro. Here we identify dimerisation hub protein Dynein Light Chain LC8-Type 1 as a CREB3L1 interacting protein that in vitro enhances CREB3L1 activation of target genes.

Phenotype and energy metabolism differ between osteoarthritic chondrocytes from male compared to female patients: Implications for sexual dimorphism in osteoarthritis development?

OBJECTIVES

The prevalence and severity of knee osteoarthritis (OA) are greater in females than males. The purpose of this study was to determine whether there is an underlying difference in the biology of OA chondrocytes between males and females.

METHODS

Chondrocytes were obtained following knee arthroplasty from male and female patients with primary OA. Phenotype marker expression, glucose and fat consumption, and rates of glycolysis and oxidative phosphorylation were compared between females and males. RNAi was used to determine the consequences of differential expression of Sry-box transcription factor 9 (SOX9) and PGC1α between males and females.

RESULTS

OA chondrocytes from male donors showed elevated ribonucleic acid (RNA) and protein levels of SOX9, elevated COL2A1 protein synthesis, higher glucose consumption, and higher usage of glycolysis compared to females. OA chondrocytes from females had higher PGC1α protein levels, higher fat consumption, and higher oxidative energy metabolism than males. Knockdown of SOX9 reduced expression of COL2A1 to a greater extent in male OA chondrocytes than females whereas knockdown of PGC1α reduced COL2A1 expression in females but not males. Expression of ACAN and the glycolytic enzyme PGK1 was also reduced in males but not females following SOX9 knockdown.

CONCLUSIONS

OA chondrocyte phenotype and energy metabolism differ between males and females. Our results indicate transcriptional control of COL2A1 differs between the two. Differences in chondrocyte biology between males and females imply the underlying mechanisms involved in OA may also differ, highlighting the need to consider sex and gender when investigating pathogenesis and potential treatments for OA.

Anoikis-related subtype and prognosis analyses based on bioinformatics, and an expression verification of ANGPTL4 based on experiments of lung adenocarcinoma

Background

Lung adenocarcinoma (LUAD) is one of the most common malignant tumors with high mortality. Anoikis resistance is an important mechanism of tumor cell proliferation and migration. Our research is devoted to exploring the role of anoikis in the diagnosis, classification, and prognosis of LUAD.

Methods

We downloaded the expression profile, mutation, and clinical data of LUAD from The Cancer Genome Atlas (TCGA) database. The "ConsensusClusterPlus" package was then used for the cluster analysis, and least absolute shrinkage and selection operator (LASSO) and multivariate Cox regression analyses were used to establish the prognostic model. We verified the reliability of the model using a Gene Expression Omnibus (GEO) data set. A gene set variation analysis (GSVA) was conducted to investigate the functional enrichment differences in the different clusters and risk groups. The CIBERSORT algorithm and a single-sample gene set enrichment analysis (ssGSEA) were used to analyze immune cell infiltration. The tumor mutation burden (TMB) and Tumor Immune Dysfunction and Exclusion (TIDE) scores were used to evaluate the patients' sensitivity to immunotherapy. Immunohistochemical staining of tissue microarrays was used to verify the correlation between ANGPTL4 expression and the clinicopathological characteristics and prognosis of LUAD patients.

Results

First, we screened 135 differentially expressed anoikis-related genes (ARGs) and 23 prognosis-related ARGs from TCGA-LUAD data set. Next, 494 LUAD samples were allocated to cluster A and cluster B based on the 23 prognosis-related ARGs. The Kaplan-Meier (K-M) analysis showed the overall survival (OS) of cluster B was better than that of cluster A. The clinicopathological characteristics and functional enrichment analyses revealed significant differences between clusters A and B. The tumor microenvironment (TME) analysis showed that cluster B had more immune cell infiltration and a higher TME score than cluster A. Subsequently, a LASSO Cox regression model of LUAD was constructed with ten ARGs. The K-M analysis showed that the low-risk patients had longer OS than the high-risk patients. The receiver operating characteristic curve, nomogram, and GEO data set verification results showed that the model had high accuracy and reliability. The level of immune cell infiltration and TME score were higher in the low-risk group than the high-risk group. The high-risk group had stronger sensitivity to immune checkpoint block therapy and weaker sensitivity to chemotherapy drugs than the low-risk group. ANGPTL4 expression was correlated with stage, tumor differentiation, tumor size, lymph node metastasis, and OS.

Conclusions

We discovered novel molecular subtypes and constructed a novel prognostic model of LUAD. Our findings provide important insights into subtype classification and the accurate survival prediction of LUAD. We also identified ANGPTL4 as a prognostic indicator of LUAD.

Deep learning of multimodal networks with topological regularization for drug repositioning

MOTIVATION

Computational techniques for drug-disease prediction are essential in enhancing drug discovery and repositioning. While many methods utilize multimodal networks from various biological databases, few integrate comprehensive multi-omics data, including transcriptomes, proteomes, and metabolomes. We introduce STRGNN, a novel graph deep learning approach that predicts drug-disease relationships using extensive multimodal networks comprising proteins, RNAs, metabolites, and compounds. We have constructed a detailed dataset incorporating multi-omics data and developed a learning algorithm with topological regularization. This algorithm selectively leverages informative modalities while filtering out redundancies.

RESULTS

STRGNN demonstrates superior accuracy compared to existing methods and has identified several novel drug effects, corroborating existing literature. STRGNN emerges as a powerful tool for drug prediction and discovery. The source code for STRGNN, along with the dataset for performance evaluation, is available at https://github.com/yuto-ohnuki/STRGNN.git .

Decoding the Role of NEIL1 Gene in DNA Repair and Lifespan: A Literature Review with Bioinformatics Analysis

Longevity, the length of an organism's lifespan, is impacted by environmental factors, metabolic processes, and genetic determinants. The base excision repair (BER) pathway is crucial for maintaining genomic integrity by repairing oxidatively modified base lesions. Nei-like DNA Glycosylase 1 (NEIL1), part of the BER pathway, is vital in repairing oxidative bases in G-rich DNA regions, such as telomeres and promoters. Hence, in this comprehensive review, it have undertaken a meticulous investigation of the intricate association between NEIL1 and longevity. The analysis delves into the multifaceted aspects of the NEIL1 gene, its various RNA transcripts, and the diverse protein isoforms. In addition, a combination of bioinformatic analysis is conducted to identify NEIL1 mutations, transcription factors, and epigenetic modifications, as well as its lncRNA/pseudogene/circRNA-miRNA-mRNA regulatory network. The findings suggest that the normal function of NEIL1 is a significant factor in human health and longevity, with defects in NEIL1 potentially leading to various cancers and related syndromes, Alzheimer's disease, obesity, and diabetes.

DNA Methylation Patterns Associated with Tinnitus in Young Adults-A Pilot Study

PURPOSE

Tinnitus, the perception of sound without any external sound source, is a prevalent hearing health concern. Mounting evidence suggests that a confluence of genetic, environmental, and lifestyle factors can influence the pathogenesis of tinnitus. We hypothesized that alteration in DNA methylation, an epigenetic modification that occurs at cytosines of cytosine-phosphate-guanine (CpG) dinucleotide sites, where a methyl group from S-adenyl methionine gets transferred to the fifth carbon of the cytosine, could contribute to tinnitus. DNA methylation patterns are tissue-specific, but the tissues involved in tinnitus are not easily accessible in humans. This pilot study used saliva as a surrogate tissue to identify differentially methylated CpG regions (DMRs) associated with tinnitus. The study was conducted on healthy young adults reporting bilateral continuous chronic tinnitus to limit the influence of age-related confounding factors and health-related comorbidities.

METHODS

The present study evaluated the genome-wide methylation levels from saliva-derived DNA samples from 24 healthy young adults with bilateral continuous chronic tinnitus (> 1 year) and 24 age, sex, and ethnicity-matched controls with no tinnitus. Genome-wide DNA methylation was evaluated for > 850,000 CpG sites using the Infinium Human Methylation EPIC BeadChip. The association analysis used the Bumphunter algorithm on 23 cases and 20 controls meeting the quality control standards. The methylation level was expressed as the area under the curve of CpG sites within DMRs.The FDR-adjusted p-value threshold of 0.05 was used to identify statistically significant DMRs associated with tinnitus.

RESULTS

We obtained 25 differentially methylated regions (DMRs) associated with tinnitus. Genes within or in the proximity of the hypermethylated DMRs related to tinnitus included LCLAT1, RUNX1, RUFY1, NUDT12, TTC23, SLC43A2, C4orf27 (STPG2), and EFCAB4B. Genes within or in the proximity of hypomethylated DMRs associated with tinnitus included HLA-DPB2, PM20D1, TMEM18, SNTG2, MUC4, MIR886, MIR596, TXNRD1, EID3, SDHAP3, HLA-DPB2, LASS3 (CERS3), C10orf11 (LRMDA), HLA-DQB1, NADK, SZRD1, MFAP2, NUP210L, TPM3, INTS9, and SLC2A14. The burden of genetic variation could explain the differences in the methylation levels for DMRs involving HLA-DPB2, HLA-DQB1, and MUC4, indicating the need for replication in large independent cohorts.

CONCLUSION

Consistent with the literature on comorbidities associated with tinnitus, we identified genes within or close to DMRs involved in auditory functions, chemical dependency, cardiovascular diseases, psychiatric conditions, immune disorders, and metabolic syndromes. These results indicate that epigenetic mechanisms could influence tinnitus, and saliva can be a good surrogate for identifying the epigenetic underpinnings of tinnitus in humans. Further research with a larger sample size is needed to identify epigenetic biomarkers and investigate their influence on the phenotypic expression of tinnitus.

CRISPRi-based screens in iAssembloids to elucidate neuron-glia interactions

The sheer complexity of the brain has complicated our ability to understand the cellular and molecular mechanisms underlying its function in health and disease. Genome-wide association studies have uncovered genetic variants associated with specific neurological phenotypes and diseases. In addition, single-cell transcriptomics have provided molecular descriptions of specific brain cell types and the changes they undergo during disease. Although these approaches provide a giant leap forward towards understanding how genetic variation can lead to functional changes in the brain, they do not establish molecular mechanisms. To address this need, we developed a 3D co-culture system termed iAssembloids (induced multi-lineage assembloids) that enables the rapid generation of homogenous neuron-glia spheroids. We characterize these iAssembloids with immunohistochemistry and single-cell transcriptomics and combine them with large-scale CRISPRi-based screens. In our first application, we ask how glial and neuronal cells interact to control neuronal death and survival. Our CRISPRi-based screens identified that GSK3β inhibits the protective NRF2-mediated oxidative stress response in the presence of reactive oxygen species elicited by high neuronal activity, which was not previously found in 2D monoculture neuron screens. We also apply the platform to investigate the role of APOE-ɛ4, a risk variant for Alzheimer's Disease, in its effect on neuronal survival. We find that APOE-ɛ4-expressing astrocytes may promote more neuronal activity as compared to APOE-ɛ3-expressing astrocytes. This platform expands the toolbox for the unbiased identification of mechanisms of cell-cell interactions in brain health and disease.

Disentangling cell-intrinsic and extrinsic factors underlying evolution

A key goal of developmental biology is to determine the extent to which cells and organs develop autonomously, as opposed to requiring interactions with other cells or environmental factors. Chimeras have played a foundational role in this by enabling qualitative classification of cell-intrinsically vs. extrinsically driven processes. Here, we extend this framework to precisely decompose evolutionary divergence in any quantitative trait into cell-intrinsic, extrinsic, and intrinsic-extrinsic interaction components. Applying this framework to thousands of gene expression levels in reciprocal rat-mouse chimeras, we found that the majority of their divergence is attributable to cell-intrinsic factors, though extrinsic factors also play an integral role. For example, a rat-like extracellular environment extrinsically up-regulates the expression of a key transcriptional regulator of the endoplasmic reticulum (ER) stress response in some but not all cell types, which in turn strongly predicts extrinsic up-regulation of its target genes and of the ER stress response pathway as a whole. This effect is also seen at the protein level, suggesting propagation through multiple regulatory levels. Applying our framework to a cellular trait, neuronal differentiation, revealed a complex interaction of intrinsic and extrinsic factors. Finally, we show that imprinted genes are dramatically mis-expressed in species-mismatched environments, suggesting that mismatch between rapidly evolving intrinsic and extrinsic mechanisms controlling gene imprinting may contribute to barriers to interspecies chimerism. Overall, our conceptual framework opens new avenues to investigate the mechanistic basis of developmental processes and evolutionary divergence across myriad quantitative traits in any multicellular organism.

A multiomic characterization of the leukemia cell line REH using short- and long-read sequencing

The B-cell acute lymphoblastic leukemia (ALL) cell line REH, with the t(12;21) ETV6::RUNX1 translocation, is known to have a complex karyotype defined by a series of large-scale chromosomal rearrangements. Taken from a 15-yr-old at relapse, the cell line offers a practical model for the study of pediatric B-ALL. In recent years, short- and long-read DNA and RNA sequencing have emerged as a complement to karyotyping techniques in the resolution of structural variants in an oncological context. Here, we explore the integration of long-read PacBio and Oxford Nanopore whole-genome sequencing, IsoSeq RNA sequencing, and short-read Illumina sequencing to create a detailed genomic and transcriptomic characterization of the REH cell line. Whole-genome sequencing clarified the molecular traits of disrupted ALL-associated genes including CDKN2A, PAX5, BTG1, VPREB1, and TBL1XR1, as well as the glucocorticoid receptor NR3C1 Meanwhile, transcriptome sequencing identified seven fusion genes within the genomic breakpoints. Together, our extensive whole-genome investigation makes high-quality open-source data available to the leukemia genomics community.

Epigenomic landscapes during prefrontal cortex development and aging in rhesus

The prefrontal cortex (PFC) is essential for higher-level cognitive functions. How epigenetic dynamics participates in PFC development and aging is largely unknown. Here, we profiled epigenomic landscapes of rhesus monkey PFCs from prenatal to aging stages. The dynamics of chromatin states, including higher-order chromatin structure, chromatin interaction and histone modifications are coordinated to regulate stage-specific gene transcription, participating in distinct processes of neurodevelopment. Dramatic changes of epigenetic signals occur around the birth stage. Notably, genes involved in neuronal cell differentiation and layer specification are pre-configured by bivalent promoters. We identified a cis-regulatory module and the transcription factors (TFs) associated with basal radial glia development, which was associated with large brain size in primates. These TFs include GLI3, CREB5 and SOX9. Interestingly, the genes associated with the basal radial glia (bRG)-associated cis-element module, such as SRY and SOX9, are enriched in sex differentiation. Schizophrenia-associated single nucleotide polymorphisms are more enriched in super enhancers (SEs) than typical enhancers, suggesting that SEs play an important role in neural network wiring. A cis-regulatory element of DBN1 is identified, which is critical for neuronal cell proliferation and synaptic neuron differentiation. Notably, the loss of distal chromatin interaction and H3K27me3 signal are hallmarks of PFC aging, which are associated with abnormal expression of aging-related genes and transposon activation, respectively. Collectively, our findings shed light on epigenetic mechanisms underlying primate brain development and aging.

Multiome Perturb-seq unlocks scalable discovery of integrated perturbation effects on the transcriptome and epigenome

Single-cell CRISPR screens link genetic perturbations to transcriptional states, but high-throughput methods connecting these induced changes to their regulatory foundations are limited. Here we introduce Multiome Perturb-seq, extending single-cell CRISPR screens to simultaneously measure perturbation-induced changes in gene expression and chromatin accessibility. We apply Multiome Perturb-seq in a CRISPRi screen of 13 chromatin remodelers in human RPE-1 cells, achieving efficient assignment of sgRNA identities to single nuclei via an improved method for capturing barcode transcripts from nuclear RNA. We organize expression and accessibility measurements into coherent programs describing the integrated effects of perturbations on cell state, finding that ARID1A and SUZ12 knockdowns induce programs enriched for developmental features. Pseudotime analysis of perturbations connects accessibility changes to changes in gene expression, highlighting the value of multimodal profiling. Overall, our method provides a scalable and simply implemented system to dissect the regulatory logic underpinning cell state.

Identifying cell type-specific transcription factor-mediated activity immune modules reveal implications for immunotherapy and molecular classification of pan-cancer

Systematic investigation of tumor-infiltrating immune (TII) cells is important to the development of immunotherapies, and the clinical response prediction in cancers. There exists complex transcriptional regulation within TII cells, and different immune cell types display specific regulation patterns. To dissect transcriptional regulation in TII cells, we first integrated the gene expression profiles from single-cell datasets, and proposed a computational pipeline to identify TII cell type-specific transcription factor (TF) mediated activity immune modules (TF-AIMs). Our analysis revealed key TFs, such as BACH2 and NFKB1 play important roles in B and NK cells, respectively. We also found some of these TF-AIMs may contribute to tumor pathogenesis. Based on TII cell type-specific TF-AIMs, we identified eight CD8+ T cell subtypes. In particular, we found the PD1 + CD8+ T cell subset and its specific TF-AIMs associated with immunotherapy response. Furthermore, the TII cell type-specific TF-AIMs displayed the potential to be used as predictive markers for immunotherapy response of cancer patients. At the pan-cancer level, we also identified and characterized six molecular subtypes across 9680 samples based on the activation status of TII cell type-specific TF-AIMs. Finally, we constructed a user-friendly web interface CellTF-AIMs (http://bio-bigdata.hrbmu.edu.cn/CellTF-AIMs/) for exploring transcriptional regulatory pattern in various TII cell types. Our study provides valuable implications and a rich resource for understanding the mechanisms involved in cancer microenvironment and immunotherapy.

Enhancer-driven gene regulatory networks inference from single-cell RNA-seq and ATAC-seq data

Deciphering the intricate relationships between transcription factors (TFs), enhancers, and genes through the inference of enhancer-driven gene regulatory networks (eGRNs) is crucial in understanding gene regulatory programs in a complex biological system. This study introduces STREAM, a novel method that leverages a Steiner forest problem model, a hybrid biclustering pipeline, and submodular optimization to infer eGRNs from jointly profiled single-cell transcriptome and chromatin accessibility data. Compared to existing methods, STREAM demonstrates enhanced performance in terms of TF recovery, TF-enhancer linkage prediction, and enhancer-gene relation discovery. Application of STREAM to an Alzheimer's disease dataset and a diffuse small lymphocytic lymphoma dataset reveals its ability to identify TF-enhancer-gene relations associated with pseudotime, as well as key TF-enhancer-gene relations and TF cooperation underlying tumor cells.

Transcription factor 7 like 2 promotes metastasis in hepatocellular carcinoma via NEDD9-mediated activation of AKT/mTOR signaling pathway

BACKGROUND

Hepatocellular carcinoma (HCC) is one of the most common malignant tumors of the digestive system, and the exact mechanism of HCC is still unclear. Transcription factor 7 like 2 (TCF7L2) plays a pivotal role in cell proliferation and stemness maintenance. However, the exact mechanism of TCF7L2 in HCC remains unclear.

METHODS

Clinical samples and public databases were used to analyze the expression and prognosis of TCF7L2 in HCC. The function of TCF7L2 in HCC was studied in vitro and in vivo. ChIP and luciferase assays were used to explore the molecular mechanism of TCF7L2. The relationship between TCF7L2 and NEDD9 was verified in HCC clinical samples by tissue microarrays.

RESULTS

The expression of TCF7L2 was upregulated in HCC, and high expression of TCF7L2 was associated with poor prognosis of HCC patients. Overexpression of TCF7L2 promoted the metastasis of HCC in vitro and in vivo, while Knockdown of TCF7L2 showed the opposite effect. Mechanically, TCF7L2 activated neural precursor cell expressed developmentally downregulated protein 9 (NEDD9) transcription by binding to the -1522/-1509 site of the NEDD9 promoter region, thereby increasing the phosphorylation levels of AKT and mTOR. The combination of TCF7L2 and NEDD9 could distinguish the survival of HCC patients.

CONCLUSIONS

This study demonstrated that TCF7L2 promotes HCC metastasis by activating AKT/mTOR pathway in a NEDD9-dependent manner, suggesting that potential of TCF7L2 and NEDD9 as prognostic markers and therapeutic targets for HCC.

The Xenorhabdus nematophila LrhA transcriptional regulator modulates production of γ-keto-N-acyl amides with inhibitory activity against mutualistic host nematode egg hatching

Xenorhabdus nematophila is a symbiotic Gammaproteobacterium that produces diverse natural products that facilitate mutualistic and pathogenic interactions in their nematode and insect hosts, respectively. The interplay between X. nematophila secondary metabolism and symbiosis stage is tuned by various global regulators. An example of such a regulator is the LysR-type protein transcription factor LrhA, which regulates amino acid metabolism and is necessary for virulence in insects and normal nematode progeny production. Here, we utilized comparative metabolomics and molecular networking to identify small molecule factors regulated by LrhA and characterized a rare γ-ketoacid (GKA) and two new N-acyl amides, GKA-Arg (1) and GKA-Pro (2) which harbor a γ-keto acyl appendage. A lrhA null mutant produced elevated levels of compound 1 and reduced levels of compound 2 relative to wild type. N-acyl amides 1 and 2 were shown to be selective agonists for the human G-protein-coupled receptors (GPCRs) C3AR1 and CHRM2, respectively. The CHRM2 agonist 2 deleteriously affected the hatch rate and length of Steinernema nematodes. This work further highlights the utility of exploiting regulators of host-bacteria interactions for the identification of the bioactive small molecule signals that they control.

IMPORTANCE

Xenorhabdus bacteria are of interest due to their symbiotic relationship with Steinernema nematodes and their ability to produce a variety of natural bioactive compounds. Despite their importance, the regulatory hierarchy connecting specific natural products and their regulators is poorly understood. In this study, comparative metabolomic profiling was utilized to identify the secondary metabolites modulated by the X. nematophila global regulator LrhA. This analysis led to the discovery of three metabolites, including an N-acyl amide that inhibited the egg hatching rate and length of Steinernema carpocapsae nematodes. These findings support the notion that X. nematophila LrhA influences the symbiosis between X. nematophila and S. carpocapsae through N-acyl amide signaling. A deeper understanding of the regulatory hierarchy of these natural products could contribute to a better comprehension of the symbiotic relationship between X. nematophila and S. carpocapsae.

Elevating PLK1 overcomes BETi resistance in prostate cancer via triggering BRD4 phosphorylation-dependent degradation in mitosis

Bromodomain-containing protein 4 (BRD4) has emerged as a promising therapeutic target in prostate cancer (PCa). Understanding the mechanisms of BRD4 stability could enhance the clinical response to BRD4-targeted therapy. In this study, we report that BRD4 protein levels are significantly decreased during mitosis in a PLK1-dependent manner. Mechanistically, we show that BRD4 is primarily phosphorylated at T1186 by the CDK1/cyclin B complex, recruiting PLK1 to phosphorylate BRD4 at S24/S1100, which are recognized by the APC/CCdh1 complex for proteasome pathway degradation. We find that PLK1 overexpression lowers SPOP mutation-stabilized BRD4, consequently rendering PCa cells re-sensitized to BRD4 inhibitors. Intriguingly, we report that sequential treatment of docetaxel and JQ1 resulted in significant inhibition of PCa. Collectively, the results support that PLK1-phosphorylated BRD4 triggers its degradation at M phase. Sequential treatment of docetaxel and JQ1 overcomes BRD4 accumulation-associated bromodomain and extra-terminal inhibitor (BETi) resistance, which may shed light on the development of strategies to treat PCa.

A de novo PRPF8 Pathogenic Variant in Transient Severe Hypophosphatemia with Delayed Puberty and Growth Failure

INTRODUCTION

Childhood hypophosphatemia is a rare condition and may be caused by malabsorption, malignancies, or genetic factors. Prolonged hypophosphatemia leads to impaired growth and radiographic signs of rickets.

METHODS

We performed a detailed clinical and genetic evaluation of an adolescent boy with repeatedly low plasma phosphate concentrations (below 0.60 mmol/L) and growth failure.

RESULTS

At 14 years, the patient presented with decelerating growth and delayed puberty. Biochemistry showed hypophosphatemia due to increased urinary phosphate loss; kidney function and vitamin D status were normal. Radiographs showed mild metaphyseal changes. A gene panel for known genetic hypophosphatemia was negative. Trio exome analysis followed by Sanger sequencing identified a pathogenic heterozygous de novo stop-gain variant in PRPF8 gene, c.5548C>T p.(Arg1850*), in the conserved RNase H homology domain. PRPF8 encodes the pre-RNA protein 8, which has a role in RNA processing. Heterozygous PRPF8 variants have been associated with retinitis pigmentosa and neurodevelopmental disorders but not with phosphate metabolism. The patient underwent growth hormone (GH) stimulation tests which confirmed GH deficiency. Head MRI indicated partially empty sella. GH treatment was started at 15 years. Surprisingly, phosphate metabolism normalized during GH treatment, suggesting that hypophosphatemia was at least partly secondary to GH deficiency.

CONCLUSION

The evaluation of an adolescent with profound long-term hypophosphatemia revealed a pituitary developmental defect associated with a stop-gain variant in PRPF8. Hypophosphatemia alleviated with GH treatment. The pathological PRPF8 variant may contribute to abnormal pituitary development; however, its role in phosphate metabolism remains uncertain.

Involvement of the choroid plexus in Alzheimer's disease pathophysiology: findings from mouse and human proteomic studies

BACKGROUND

Structural and functional changes of the choroid plexus (ChP) have been reported in Alzheimer's disease (AD). Nonetheless, the role of the ChP in the pathogenesis of AD remains largely unknown. We aim to unravel the relation between ChP functioning and core AD pathogenesis using a unique proteomic approach in mice and humans.

METHODS

We used an APP knock-in mouse model, APPNL-G-F, exhibiting amyloid pathology, to study the association between AD brain pathology and protein changes in mouse ChP tissue and CSF using liquid chromatography mass spectrometry. Mouse proteomes were investigated at the age of 7 weeks (n = 5) and 40 weeks (n = 5). Results were compared with previously published human AD CSF proteomic data (n = 496) to identify key proteins and pathways associated with ChP changes in AD.

RESULTS

ChP tissue proteome was dysregulated in APPNL-G-F mice relative to wild-type mice at both 7 and 40 weeks. At both ages, ChP tissue proteomic changes were associated with epithelial cells, mitochondria, protein modification, extracellular matrix and lipids. Nonetheless, some ChP tissue proteomic changes were different across the disease trajectory; pathways related to lysosomal function, endocytosis, protein formation, actin and complement were uniquely dysregulated at 7 weeks, while pathways associated with nervous system, immune system, protein degradation and vascular system were uniquely dysregulated at 40 weeks. CSF proteomics in both mice and humans showed similar ChP-related dysregulated pathways.

CONCLUSIONS

Together, our findings support the hypothesis of ChP dysfunction in AD. These ChP changes were related to amyloid pathology. Therefore, the ChP could become a novel promising therapeutic target for AD.

Parkinson's disease-associated shifts between DNA methylation and DNA hydroxymethylation in human brain in PD-related genes, including PARK19 (DNAJC6) and PTPRN2 (IA-2β)

Background

The majority of Parkinson's disease (PD) cases are due to a complex interaction between aging, genetics, and environmental factors; epigenetic mechanisms are thought to act as important mediators of these risk factors. While multiple studies to date have explored the role of DNA modifications in PD, few focus on 5-hydroxymethylcytosine (5hmC). Because 5hmC occurs at its highest levels in the brain and is thought to be particularly important in the central nervous system, particularly in the response to neurotoxicants, it is important to explore the potential role of 5hmC in PD. This study expands on our previously published epigenome-wide association study (EWAS) performed on DNA isolated from neuron-enriched nuclei from human postmortem parietal cortex from the Banner Sun Health Research Institute Brain Bank. The study aimed to identify paired changes in 5hmC and 5mC in PD in enriched neuronal nuclei isolated from PD post-mortem parietal cortex and age- and sex-matched controls. We performed oxidative bisulfite (oxBS) conversion and paired it with our previously published bisulfite (BS)-based EWAS on the same samples to identify cytosines with significant shifts between these two related epigenetic marks. Interaction differentially modified cytosines (iDMCs) were identified using our recently published mixed-effects model for co-analyzing βmC and βhmC data.

Results

We identified 1,030 iDMCs with paired changes in 5mC and 5hmC (FDR < 0.05) that map to 695 genes, including PARK19 (DNAJC6), a familial PD gene, and PTPRN2 (IA-2), which has been previously implicated in PD in both epigenetic and mechanistic studies. The majority of iDMC-containing genes have not previously been implicated in PD and were not identified in our previous BS-based EWAS.

Conclusions

These data potentially link epigenetic regulation of the PARK19 and PTPRN2 loci in the pathogenesis of idiopathic PD. In addition, iDMC-containing genes have known functions in synaptic formation and function, cell cycle and senescence, neuroinflammation, and epigenetic regulation. These data suggest that there are significant shifts between 5mC and 5hmC associated with PD in genes relevant to PD pathogenesis that are not captured by analyzing BS-based data alone or by analyzing each mark as a distinct dataset.

Analysis of the Actions of RARγ Agonists on Growing Osteochondromas in a Mouse Model

The actions of the retinoic acid nuclear receptor gamma (RARγ) agonist, palovarotene, on pre-existing osteochondromas were investigated using a mouse multiple osteochondroma model. This approach was based on the knowledge that patients often present to the clinic after realizing the existence of osteochondroma masses, and the findings from preclinical investigations are the effects of drugs on the initial formation of osteochondromas. Systemic administration of palovarotene, with increased doses (from 1.76 to 4.0 mg/kg) over time, fully inhibited tumor growth, keeping the tumor size (0.31 ± 0.049 mm3) similar to the initial size (0.27 ± 0.031 mm3, p = 0.66) while the control group tumor grew (1.03 ± 0.23 mm3, p = 0.023 to the drug-treated group). Nanoparticle (NP)-based local delivery of the RARγ agonist also inhibited the growth of osteochondromas at an early stage (Control: 0.52 ± 0.11 mm3; NP: 0.26 ± 0.10, p = 0.008). Transcriptome analysis revealed that the osteoarthritis pathway was activated in cultured chondrocytes treated with palovarotene (Z-score = 2.29), with the upregulation of matrix catabolic genes and the downregulation of matrix anabolic genes, consistent with the histology of palovarotene-treated osteochondromas. A reporter assay performed in cultured chondrocytes demonstrated that the Stat3 pathway, but not the Stat1/2 pathway, was stimulated by RARγ agonists. The activation of Stat3 by palovarotene was confirmed using immunoblotting and immunohistochemistry. These findings suggest that palovarotene treatment is effective against pre-existing osteochondromas and that the Stat3 pathway is involved in the antitumor actions of palovarotene.

CSF proteomic profiles of neurodegeneration biomarkers in Alzheimer's disease

INTRODUCTION

We aimed to unravel the underlying pathophysiology of the neurodegeneration (N) markers neurogranin (Ng), neurofilament light (NfL), and hippocampal volume (HCV), in Alzheimer's disease (AD) using cerebrospinal fluid (CSF) proteomics.

METHODS

Individuals without dementia were classified as A+ (CSF amyloid beta [Aβ]42), T+ (CSF phosphorylated tau181), and N+ or N- based on Ng, NfL, or HCV separately. CSF proteomics were generated and compared between groups using analysis of covariance.

RESULTS

Only a few individuals were A+T+Ng-. A+T+Ng+ and A+T+NfL+ showed different proteomic profiles compared to A+T+Ng- and A+T+NfL-, respectively. Both Ng+ and NfL+ were associated with neuroplasticity, though in opposite directions. Compared to A+T+HCV-, A+T+HCV+ showed few proteomic changes, associated with oxidative stress.

DISCUSSION

Different N markers are associated with distinct neurodegenerative processes and should not be equated. N markers may differentially complement disease staging beyond amyloid and tau. Our findings suggest that Ng may not be an optimal N marker, given its low incongruency with tau pathophysiology.

HIGHLIGHTS

In Alzheimer's disease, neurogranin (Ng)+, neurofilament light (NfL)+, and hippocampal volume (HCV)+ showed differential protein expression in cerebrospinal fluid. Ng+ and NfL+ were associated with neuroplasticity, although in opposite directions. HCV+ showed few proteomic changes, related to oxidative stress. Neurodegeneration (N) markers may differentially refine disease staging beyond amyloid and tau. Ng might not be an optimal N marker, as it relates more closely to tau.

A liver immune rheostat regulates CD8 T cell immunity in chronic HBV infection

Chronic hepatitis B virus (HBV) infection affects 300 million patients worldwide1,2, in whom virus-specific CD8 T cells by still ill-defined mechanisms lose their function and cannot eliminate HBV-infected hepatocytes3-7. Here we demonstrate that a liver immune rheostat renders virus-specific CD8 T cells refractory to activation and leads to their loss of effector functions. In preclinical models of persistent infection with hepatotropic viruses such as HBV, dysfunctional virus-specific CXCR6+ CD8 T cells accumulated in the liver and, as a characteristic hallmark, showed enhanced transcriptional activity of cAMP-responsive element modulator (CREM) distinct from T cell exhaustion. In patients with chronic hepatitis B, circulating and intrahepatic HBV-specific CXCR6+ CD8 T cells with enhanced CREM expression and transcriptional activity were detected at a frequency of 12-22% of HBV-specific CD8 T cells. Knocking out the inhibitory CREM/ICER isoform in T cells, however, failed to rescue T cell immunity. This indicates that CREM activity was a consequence, rather than the cause, of loss in T cell function, further supported by the observation of enhanced phosphorylation of protein kinase A (PKA) which is upstream of CREM. Indeed, we found that enhanced cAMP-PKA-signalling from increased T cell adenylyl cyclase activity augmented CREM activity and curbed T cell activation and effector function in persistent hepatic infection. Mechanistically, CD8 T cells recognizing their antigen on hepatocytes established close and extensive contact with liver sinusoidal endothelial cells, thereby enhancing adenylyl cyclase-cAMP-PKA signalling in T cells. In these hepatic CD8 T cells, which recognize their antigen on hepatocytes, phosphorylation of key signalling kinases of the T cell receptor signalling pathway was impaired, which rendered them refractory to activation. Thus, close contact with liver sinusoidal endothelial cells curbs the activation and effector function of HBV-specific CD8 T cells that target hepatocytes expressing viral antigens by means of the adenylyl cyclase-cAMP-PKA axis in an immune rheostat-like fashion.

SFINN: inferring gene regulatory network from single-cell and spatial transcriptomic data with shared factor neighborhood and integrated neural network

MOTIVATION

The rise of single-cell RNA sequencing (scRNA-seq) technology presents new opportunities for constructing detailed cell type-specific gene regulatory networks (GRNs) to study cell heterogeneity. However, challenges caused by noises, technical errors, and dropout phenomena in scRNA-seq data pose significant obstacles to GRN inference, making the design of accurate GRN inference algorithms still essential. The recent growth of both single-cell and spatial transcriptomic sequencing data enables the development of supervised deep learning methods to infer GRNs on these diverse single-cell datasets.

RESULTS

In this study, we introduce a novel deep learning framework based on shared factor neighborhood and integrated neural network (SFINN) for inferring potential interactions and causalities between transcription factors and target genes from single-cell and spatial transcriptomic data. SFINN utilizes shared factor neighborhood to construct cellular neighborhood network based on gene expression data and additionally integrates cellular network generated from spatial location information. Subsequently, the cell adjacency matrix and gene pair expression are fed into an integrated neural network framework consisting of a graph convolutional neural network and a fully-connected neural network to determine whether the genes interact. Performance evaluation in the tasks of gene interaction and causality prediction against the existing GRN reconstruction algorithms demonstrates the usability and competitiveness of SFINN across different kinds of data. SFINN can be applied to infer GRNs from conventional single-cell sequencing data and spatial transcriptomic data.

AVAILABILITY AND IMPLEMENTATION

SFINN can be accessed at GitHub: https://github.com/JGuan-lab/SFINN.

Identifying new cancer genes based on the integration of annotated gene sets via hypergraph neural networks

MOTIVATION

Identifying cancer genes remains a significant challenge in cancer genomics research. Annotated gene sets encode functional associations among multiple genes, and cancer genes have been shown to cluster in hallmark signaling pathways and biological processes. The knowledge of annotated gene sets is critical for discovering cancer genes but remains to be fully exploited.

RESULTS

Here, we present the DIsease-Specific Hypergraph neural network (DISHyper), a hypergraph-based computational method that integrates the knowledge from multiple types of annotated gene sets to predict cancer genes. First, our benchmark results demonstrate that DISHyper outperforms the existing state-of-the-art methods and highlight the advantages of employing hypergraphs for representing annotated gene sets. Second, we validate the accuracy of DISHyper-predicted cancer genes using functional validation results and multiple independent functional genomics data. Third, our model predicts 44 novel cancer genes, and subsequent analysis shows their significant associations with multiple types of cancers. Overall, our study provides a new perspective for discovering cancer genes and reveals previously undiscovered cancer genes.

AVAILABILITY AND IMPLEMENTATION

DISHyper is freely available for download at https://github.com/genemine/DISHyper.

Elucidating prognosis in cervical squamous cell carcinoma and endocervical adenocarcinoma: a novel anoikis-related gene signature model

Background

Cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC) are among the most prevalent gynecologic malignancies globally. The prognosis is abysmal once cervical cancer progresses to lymphatic metastasis. Anoikis, a specialized form of apoptosis induced by loss of cell adhesion to the extracellular matrix, plays a critical role. The prediction model based on anoikis-related genes (ARGs) expression and clinical data could greatly aid clinical decision-making. However, the relationship between ARGs and CESC remains unclear.

Methods

ARGs curated from the GeneCards and Harmonizome portals were instrumental in delineating CESC subtypes and in developing a prognostic framework for patients afflicted with this condition. We further delved into the intricacies of the immune microenvironment and pathway enrichment across the identified subtypes. Finally, our efforts culminated in the creation of an innovative nomogram that integrates ARGs. The utility of this prognostic tool was underscored by Decision Curve Analysis (DCA), which illuminate its prospective benefits in guiding clinical interventions.

Results

In our study, We discerned a set of 17 survival-pertinent, anoikis-related differentially expressed genes (DEGs) in CESC, from which nine were meticulously selected for the construction of prognostic models. The derived prognostic risk score was subsequently validated as an autonomous prognostic determinant. Through comprehensive functional analyses, we observed distinct immune profiles and drug response patterns among divergent prognostic stratifications. Further, we integrated the risk scores with the clinicopathological characteristics of CESC to develop a robust nomogram. DCA corroborated the utility of our model, demonstrating its potential to enhance patient outcomes through tailored clinical treatment strategies.

Conclusion

The predictive signature, encompassing nine pivotal genes, alongside the meticulously constructed nomogram developed in this research, furnishes clinicians with a sophisticated tool for tailoring treatment strategies to individual patients diagnosed with CESC.

Endonucleolytic RNA cleavage drives changes in gene expression during the innate immune response

Viral infection triggers several double-stranded RNA (dsRNA) sensors that lead to changes in gene expression in the cell. One of these sensors activates an endonuclease, ribonuclease L (RNase L), that cleaves single-stranded RNA. However, how the resultant widespread RNA fragmentation affects gene expression is not fully understood. Here, we show that this fragmentation induces the ribotoxic stress response via ZAKα, potentially through stalled ribosomes and/or ribosome collisions. The p38 and JNK pathways that are activated as part of this response promote outcomes that inhibit the virus, such as programmed cell death. We also show that RNase L limits the translation of stress-responsive genes. Intriguingly, we found that the activity of the generic endonuclease, RNase A, recapitulates many of the same molecular phenotypes as activated RNase L, demonstrating how widespread RNA cleavage can evoke an antiviral program.

The non-mitotic role of HMMR in regulating the localization of TPX2 and the dynamics of microtubules in neurons

A functional nervous system is built upon the proper morphogenesis of neurons to establish the intricate connection between them. The microtubule cytoskeleton is known to play various essential roles in this morphogenetic process. While many microtubule-associated proteins (MAPs) have been demonstrated to participate in neuronal morphogenesis, the function of many more remains to be determined. This study focuses on a MAP called HMMR in mice, which was originally identified as a hyaluronan binding protein and later found to possess microtubule and centrosome binding capacity. HMMR exhibits high abundance on neuronal microtubules and altering the level of HMMR significantly affects the morphology of neurons. Instead of confining to the centrosome(s) like cells in mitosis, HMMR localizes to microtubules along axons and dendrites. Furthermore, transiently expressing HMMR enhances the stability of neuronal microtubules and increases the formation frequency of growing microtubules along the neurites. HMMR regulates the microtubule localization of a non-centrosomal microtubule nucleator TPX2 along the neurite, offering an explanation for how HMMR contributes to the promotion of growing microtubules. This study sheds light on how cells utilize proteins involved in mitosis for non-mitotic functions.

Systematic Analysis of Biological Processes Reveals Gene Co-expression Modules Driving Pathway Dysregulation in Alzheimer's Disease

Alzheimer's disease (AD) manifests as a complex systems pathology with intricate interplay among various genes and biological processes. Traditional differential gene expression (DEG) analysis, while commonly employed to characterize AD-driven perturbations, does not sufficiently capture the full spectrum of underlying biological processes. Utilizing single-nucleus RNA-sequencing data from postmortem brain samples across key regions-middle temporal gyrus, superior frontal gyrus, and entorhinal cortex-we provide a comprehensive systematic analysis of disrupted processes in AD. We go beyond the DEG-centric analysis by integrating pathway activity analysis with weighted gene co-expression patterns to comprehensively map gene interconnectivity, identifying region- and cell-type-specific drivers of biological processes associated with AD. Our analysis reveals profound modular heterogeneity in neurons and glia as well as extensive AD-related functional disruptions. Co-expression networks highlighted the extended involvement of astrocytes and microglia in biological processes beyond neuroinflammation, such as calcium homeostasis, glutamate regulation, lipid metabolism, vesicle-mediated transport, and TOR signaling. We find limited representation of DEGs within dysregulated pathways across neurons and glial cells, suggesting that differential gene expression alone may not adequately represent the disease complexity. Further dissection of inferred gene modules revealed distinct dynamics of hub DEGs in neurons versus glia, suggesting that DEGs exert more impact on neurons compared to glial cells in driving modular dysregulations underlying perturbed biological processes. Interestingly, we observe an overall downregulation of astrocyte and microglia modules across all brain regions in AD, indicating a prevailing trend of functional repression in glial cells across these regions. Notable genes from the CALM and HSP90 families emerged as hub genes across neuronal modules in all brain regions, suggesting conserved roles as drivers of synaptic dysfunction in AD. Our findings demonstrate the importance of an integrated, systems-oriented approach combining pathway and network analysis to comprehensively understand the cell-type-specific roles of genes in AD-related biological processes.

Inhibiting AGTR1 reduces AML burden and protects the heart from cardiotoxicity in mouse models

Clinical treatment of acute myeloid leukemia (AML) largely relies on intensive chemotherapy. However, the application of chemotherapy is often hindered by cardiotoxicity. Patient sequence data revealed that angiotensin II receptor type 1 (AGTR1) is a shared target between AML and cardiovascular disease (CVD). We found that inhibiting AGTR1 sensitized AML to chemotherapy and protected the heart against chemotherapy-induced cardiotoxicity in a human AML cell-transplanted mouse model. These effects were regulated by the AGTR1-Notch1 axis in AML cells and cardiomyocytes from mice. In mouse cardiomyocytes, AGTR1 was hyperactivated by AML and chemotherapy. AML leukemogenesis increased the expression of the angiotensin-converting enzyme and led to increased production of angiotensin II, the ligand of AGTR1, in an MLL-AF9-driven AML mouse model. In this model, the AGTR1-Notch1 axis regulated a variety of genes involved with cell stemness and chemotherapy resistance. AML cell stemness was reduced after Agtr1a deletion in the mouse AML cell transplant model. Mechanistically, Agtr1a deletion decreased γ-secretase formation, which is required for transmembrane Notch1 cleavage and release of the Notch1 intracellular domain into the nucleus. Using multiomics, we identified AGTR1-Notch1 signaling downstream genes and found decreased binding between these gene sequences with Notch1 and chromatin enhancers, as well as increased binding with silencers. These findings describe an AML/CVD association that may be used to improve AML treatment.

Plasma Extracellular MicroRNAs Associated With Cardiovascular Disease Risk Factors in Middle-Aged and Older Adults

BACKGROUND

Extracellular microRNAs (miRNAs) are a class of noncoding RNAs that remain stable in the extracellular milieu, where they contribute to various physiological and pathological processes by facilitating intercellular signaling. Previous studies have reported associations between miRNAs and cardiovascular diseases (CVDs); however, the plasma miRNA signatures of CVD and its risk factors have not been fully elucidated at the population level.

METHODS AND RESULTS

Plasma miRNA levels were measured in 4440 FHS (Framingham Heart Study) participants. Linear regression analyses were conducted to test the cross-sectional associations of each miRNA with 8 CVD risk factors. Prospective analyses of the associations of miRNAs with new-onset obesity, hypertension, type 2 diabetes, CVD, and all-cause mortality were conducted using proportional hazards regression. Replication was carried out in 1999 RS (Rotterdam Study) participants. Pathway enrichment analyses were conducted and target genes were predicted for miRNAs associated with ≥5 risk factors in the FHS. In the FHS, 6 miRNAs (miR-193b-3p, miR-122-5p, miR-365a-3p, miR-194-5p, miR-192-5p, and miR-193a-5p) were associated with ≥5 risk factors. This miRNA signature was enriched for pathways associated with CVD and several genes annotated to these pathways were predicted targets of the identified miRNAs. Furthermore, miR-193b-3p, miR-194-5p, and miR-193a-5p were each associated with ≥2 risk factors in the RS. Prospective analysis revealed 8 miRNAs associated with all-cause mortality in the FHS.

CONCLUSIONS

These findings highlight associations between miRNAs and CVD risk factors that may provide valuable insights into the underlying pathogenesis of CVD.

Single-cell signatures identify microenvironment factors in tumors associated with patient outcomes

The cellular components of tumors and their microenvironment play pivotal roles in tumor progression, patient survival, and the response to cancer treatments. Unveiling a comprehensive cellular profile within bulk tumors via single-cell RNA sequencing (scRNA-seq) data is crucial, as it unveils intrinsic tumor cellular traits that elude identification through conventional cancer subtyping methods. Our contribution, scBeacon, is a tool that derives cell-type signatures by integrating and clustering multiple scRNA-seq datasets to extract signatures for deconvolving unrelated tumor datasets on bulk samples. Through the employment of scBeacon on the The Cancer Genome Atlas (TCGA) cohort, we find cellular and molecular attributes within specific tumor categories, many with patient outcome relevance. We developed a tumor cell-type map to visually depict the relationships among TCGA samples based on the cell-type inferences.

Progressive senescence programs induce intrinsic vulnerability to aging-related female breast cancer

Cancer incidence escalates exponentially with advancing age; however, the underlying mechanism remains unclear. In this study, we build a chronological molecular clock at single-cell transcription level with a mammary stem cell-enriched population to depict physiological aging dynamics in female mice. We find that the mammary aging process is asynchronous and progressive, initiated by an early senescence program, succeeded by an entropic late senescence program with elevated cancer associated pathways, vulnerable to cancer predisposition. The transition towards senescence program is governed by a stem cell factor Bcl11b, loss of which accelerates mammary ageing with enhanced DMBA-induced tumor formation. We have identified a drug TPCA-1 that can rejuvenate mammary cells and significantly reduce aging-related cancer incidence. Our findings establish a molecular portrait of progressive mammary cell aging and elucidate the transcriptional regulatory network bridging mammary aging and cancer predisposition, which has potential implications for the management of cancer prevalence in the aged.

Repurposing Castanea sativa Spiny Burr By-Products Extract as a Potentially Effective Anti-Inflammatory Agent for Novel Future Biotechnological Applications

The concept of a "circular bioeconomy" holds great promise for the health, cosmetic, and nutrition sectors by re-using Castanea sativa (Mill.) by-products. This sustainable resource is rich in bioactive secondary metabolites with antioxidant and anti-inflammatory properties. By transforming these by-products into high-value products for human health, we can promote sustainable economic growth and reduce the environmental impact of traditional waste disposal, adding value to previously underutilized resources. In the present study, we investigated the antioxidant capacity, phytochemical composition, and in vitro antioxidant and anti-inflammatory activity of C. sativa burr (CSB) aqueous extract. The spectrophotometric study revealed high total phenolic content (TPC) values with significant antioxidant and anti-radical properties. Using UPLC-MS/MS techniques, the phytochemical investigation identified 56 metabolites, confirming the presence of phenolic compounds in CSBs. In addition, CSBs significantly downregulated pro-inflammatory mediators in LPS-stimulated RAW 264.7 macrophage cells without significant cell toxicity. Lastly, in silico studies pinpointed three kinases from RAW 264.7 cells as binding partners with ellagic acid, the predominant compound found in our extract. These findings strongly advocate for the recycling and valorization of C. sativa by-products, challenging their conventional classification as mere "waste".

Context-Dependent Regulation of Peripheral Nerve Abundance by the PI3K Pathway in the Tumor Microenvironment of Head and Neck Squamous Cell Carcinoma

Recent studies have highlighted neurons and their associated Schwann cells (SCs) as key regulators of cancer development. However, the mode of their interaction with tumor cells or other components of the tumor microenvironment (TME) remains elusive. We established an SC-related 43-gene set as a surrogate for peripheral nerves in the TME. Head and neck squamous cell carcinoma (HNSCC) from The Cancer Genome Atlas (TCGA) were classified into low, intermediate and high SC score groups based on the expression of this gene set. Perineural invasion (PNI) and TGF-β signaling were hallmarks of SChigh tumors, whereas SClow tumors were enriched for HPV16-positive OPSCC and higher PI3K-MTOR activity. The latter activity was partially explained by a higher frequency of PTEN mutation and PIK3CA copy number gain. The inverse association between PI3K-MTOR activity and peripheral nerve abundance was context-dependent and influenced by the TP53 mutation status. An in silico drug screening approach highlighted the potential vulnerabilities of HNSCC with variable SC scores and predicted a higher sensitivity of SClow tumors to DNA topoisomerase inhibitors. In conclusion, we have established a tool for assessing peripheral nerve abundance in the TME and provided new clinical and biological insights into their regulation. This knowledge may pave the way for new therapeutic strategies and impart proof of concept in appropriate preclinical models.

System analysis based on Anoikis-related genes identifies MAPK1 as a novel therapy target for osteosarcoma with neoadjuvant chemotherapy

BACKGROUND

Osteosarcoma (OS) is the most common bone malignant tumor in children, and its prognosis is often poor. Anoikis is a unique mode of cell death.However, the effects of Anoikis in OS remain unexplored.

METHOD

Differential analysis of Anoikis-related genes was performed based on the metastatic and non-metastatic groups. Then LASSO logistic regression and SVM-RFE algorithms were applied to screen out the characteristic genes. Later, Univariate and multivariate Cox regression was conducted to identify prognostic genes and further develop the Anoikis-based risk score. In addition, correlation analysis was performed to analyze the relationship between tumor microenvironment, drug sensitivity, and prognostic models.

RESULTS

We established novel Anoikis-related subgroups and developed a prognostic model based on three Anoikis-related genes (MAPK1, MYC, and EDIL3). The survival and ROC analysis results showed that the prognostic model was reliable. Besides, the results of single-cell sequencing analysis suggested that the three prognostic genes were closely related to immune cell infiltration. Subsequently, aberrant expression of two prognostic genes was identified in osteosarcoma cells. Nilotinib can promote the apoptosis of osteosarcoma cells and down-regulate the expression of MAPK1.

CONCLUSIONS

We developed a novel Anoikis-related risk score model, which can assist clinicians in evaluating the prognosis of osteosarcoma patients in clinical practice. Analysis of the tumor immune microenvironment and chemotherapeutic drug sensitivity can provide necessary insights into subsequent mechanisms. MAPK1 may be a valuable therapeutic target for neoadjuvant chemotherapy in osteosarcoma.

Chronic TNF exposure induces glucocorticoid-like immunosuppression in the alveolar macrophages of aged mice that enhances their susceptibility to pneumonia

Chronic low-grade inflammation, particularly elevated tumor necrosis factor (TNF) levels, occurs due to advanced age and is associated with greater susceptibility to infection. One reason for this is age-dependent macrophage dysfunction (ADMD). Herein, we use the adoptive transfer of alveolar macrophages (AM) from aged mice into the airway of young mice to show that inherent age-related defects in AM were sufficient to increase the susceptibility to Streptococcus pneumoniae, a Gram-positive bacterium and the leading cause of community-acquired pneumonia. MAPK phosphorylation arrays using AM lysates from young and aged wild-type (WT) and TNF knockout (KO) mice revealed multilevel TNF-mediated suppression of kinase activity in aged mice. RNAseq analyses of AM validated the suppression of MAPK signaling as a consequence of TNF during aging. Two regulatory phosphatases that suppress MAPK signaling, Dusp1 and Ptprs, were confirmed to be upregulated with age and as a result of TNF exposure both ex vivo and in vitro. Dusp1 is known to be responsible for glucocorticoid-mediated immune suppression, and dexamethasone treatment increased Dusp1 and Ptprs expression in cells and recapitulated the ADMD phenotype. In young mice, treatment with dexamethasone increased the levels of Dusp1 and Ptprs and their susceptibility to infection. TNF-neutralizing antibody reduced Dusp1 and Ptprs levels in AM from aged mice and reduced pneumonia severity following bacterial challenge. We conclude that chronic exposure to TNF increases the expression of the glucocorticoid-associated MAPK signaling suppressors, Dusp1 and Ptprs, which inhibits AM activation and increases susceptibility to bacterial pneumonia in older adults.

Chemo-small extracellular vesicles released in cisplatin-resistance ovarian cancer cells are regulated by the lysosomal function

Chemoresistance is a common problem in ovarian cancer (OvCa) treatment, where resistant cells, in response to chemotherapy, secrete small extracellular vesicles (sEVs), known as chemo-sEVs, that transfer resistance to recipient cells. sEVs are formed as intraluminal vesicles (ILVs) within multivesicular endosomes (MVEs), whose trafficking is regulated by Ras-associated binding (RAB) GTPases that mediate sEVs secretion or lysosomal degradation. A decrease in lysosomal function can promote sEVs secretion, but the relationship between MVEs trafficking pathways and sEVs secretion in OvCa chemoresistance is unclear. Here, we show that A2780cis cisplatin (CCDP) resistant OvCa cells had an increased number of MVEs and ILVs structures, higher levels of Endosomal Sorting Complex Required for Transport (ESCRTs) machinery components, and RAB27A compared to A2780 CDDP-sensitive OvCa cells. CDDP promoted the secretion of chemo-sEVs in A2780cis cells, enriched in DNA damage response proteins. A2780cis cells exhibited poor lysosomal function with reduced levels of RAB7, essential in MVEs-Lysosomal trafficking. The silencing of RAB27A in A2780cis cells prevents the Chemo-EVs secretion, reduces its chemoresistance and restores lysosomal function and levels of RAB7, switching them into an A2780-like cellular phenotype. Enhancing lysosomal function with rapamycin reduced chemo-sEVs secretion. Our results suggest that adjusting the balance between secretory MVEs and lysosomal MVEs trafficking could be a promising strategy for overcoming CDDP chemoresistance in OvCa.