Enrichr: interactive and collaborative HTML5 gene list enrichment analysis tool

Transcriptomic Evaluation of a Stress Vulnerability Network Using Single-Cell RNA Sequencing in Mouse Prefrontal Cortex

BACKGROUND

Increased vulnerability to stress is a major risk factor for several mood disorders, including major depressive disorder. Although cellular and molecular mechanisms associated with depressive behaviors following stress have been identified, little is known about the mechanisms that confer the vulnerability that predisposes individuals to future damage from chronic stress.

METHODS

We used multisite in vivo neurophysiology in freely behaving male and female C57BL/6 mice (n = 12) to measure electrical brain network activity previously identified as indicating a latent stress vulnerability brain state. We combined this neurophysiological approach with single-cell RNA sequencing of the prefrontal cortex to identify distinct transcriptomic differences between groups of mice with inherent high and low stress vulnerability.

RESULTS

We identified hundreds of differentially expressed genes (padjusted < .05) across 5 major cell types in animals with high and low stress vulnerability brain network activity. This unique analysis revealed that GABAergic (gamma-aminobutyric acidergic) neuron gene expression contributed most to the network activity of the stress vulnerability brain state. Upregulation of mitochondrial and metabolic pathways also distinguished high and low vulnerability brain states, especially in inhibitory neurons. Importantly, genes that were differentially regulated with vulnerability network activity significantly overlapped (above chance) with those identified by genome-wide association studies as having single nucleotide polymorphisms significantly associated with depression as well as genes more highly expressed in postmortem prefrontal cortex of patients with major depressive disorder.

CONCLUSIONS

This is the first study to identify cell types and genes involved in a latent stress vulnerability state in the brain.

Pax3/7 gene function in Oikopleura dioica supports a neuroepithelial-like origin for its house-making Fol territory

Larvacean tunicates feature a spectacular innovation not seen in other animals - the trunk oikoplastic epithelium (OE). This epithelium produces a house, a large and complex extracellular structure used for filtering and concentrating food particles. Previously we identified several homeobox transcription factor genes expressed during early OE patterning. Among these are two Pax3/7 copies that we named pax37A and pax37B. The vertebrate homologs, PAX3 and PAX7 are involved in developmental processes related to neural crest and muscles. In the ascidian tunicate Ciona intestinalis, Pax3/7 plays a role in the development of cells deriving from the neural plate border, including trunk epidermal sensory neurons and tail nerve cord neurons, as well as in the neural tube closure. Here we have investigated the roles of Oikopleura dioica pax37A and pax37B in the development of the OE, by using CRISPR-Cas9 mutant lines and analyzing scRNA-seq data from wild-type animals. We found that pax37B but not pax37A is essential for the differentiation of cell fields that produce the food concentrating filter of the house: the anterior Fol, giant Fol and Nasse cells. Trajectory analysis supported a neuroepithelial-like or a preplacodal ectoderm transcriptional signature in these cells. We propose that the highly specialized secretory epithelial cells of the Fol region either maintained or evolved neuroepithelial features. This is supported by a fragmented gene regulatory network involved in their development that also operates in ascidian epidermal neurons.

Uncovering the role of aquaporin and chromobox family members as potential biomarkers in head and neck squamous cell carcinoma via integrative multiomics and in silico approach

Head and neck squamous cell carcinoma (HNSC) is a diverse group of tumors arising from oral cavity, oropharynx, larynx, and hypopharynx squamous epithelium, posing significant morbidity. Aquaporins (AQPs) are membrane proteins forming water channels, some associated with carcinomas. Chromobox (CBX) family is known to modulate physiological and oncological processes. In our study, we analyzed AQPs and CBXs having significant expression followed by their prognostic and mutational assessment. Next, we performed enrichment and tumor infiltration analysis followed by HPA validation. Lastly, we established a 3-node miRNA-TF-mRNA regulatory network and performed protein-protein docking of the highest-degree subnetwork motif between TF and mRNA. Significant upregulation of CBX3/2 and downregulation of AQP3/5/7 correlated with poor overall survival (OS) in HNSC patients. The most significant pathway, GO-BP, GO-MF, and GO-CC terms associated with AQP3 and CBX3 were passive transport by aquaporins, response to vitamin, glycerol channel activity, and condensed chromosome, centromeric region. AQP3 negatively correlated withCD 4 + T cells, positively withCD 8 + T cells and B cells, and negatively with tumor purity, whereas CBX3 positively correlated withCD 4 + T cells, negatively withCD 8 + T cells and B cells, and positively with tumor purity. Three-node miRNA-TF-mRNA regulatory network revealed a highest-degree subnetwork motif comprising one TF (SMAD3), one miRNA (miR-423-5p), and one mRNA (AQP3). Protein-protein interaction studies suggested a direct interaction between AQP3 and Smad3 proteins. We concluded that AQP3 and CBX3 hold potential as treatment strategies and individual prognostic biomarkers, while further protein-protein interaction studies of AQP3 could offer insights into its interactions with Smad3 proteins.

Integrative Bioinformatics Analysis: Unraveling Variant Signatures and Single-Nucleotide Polymorphism Markers Associated with 5-FU-Based Chemotherapy Resistance in Colorectal Cancer Patients

BACKGROUND

Drug resistance in colorectal cancer (CRC) is modulated by multiple molecular factors, which can be ascertained through genetic investigation. Single nucleotide polymorphisms (SNPs) within key genes have the potential to impair the efficacy of chemotherapeutic agents such as 5-fluorouracil (5-FU). Therefore, the identification of SNPs linked to drug resistance can significantly contribute to the advancement of tailored therapeutic approaches and the enhancement of treatment outcomes in patients with CRC.

MATERIAL AND METHOD

To identify dysregulated genes in 5-FU-based chemotherapy responder or non-responder CRC patients, a meta-analysis was performed. Next, the protein-protein interaction (PPI) network of the identified genes was analyzed using the STRING database. The most significant module was chosen for further analysis. In addition, a literature review was conducted to identify drug resistance-related genes. Enrichment analysis was conducted to validate the main module genes and the genes identified from the literature review. The associations between SNPs and drug resistance were investigated, and the consequences of missense variants were assessed using in silico tools.

RESULT

The meta-analysis identified 796 dysregulated genes. Then, to conduct PPI analysis and enrichment analysis, we were able to discover 23 genes that are intricately involved in the cell cycle pathway. Consequently, these 23 genes were chosen for SNP analysis. By using the dbSNP database and ANNOVAR, we successfully detected and labeled SNPs in these specific genes. Additionally, after careful exclusion of SNPs with allele frequencies below 0.01, we evaluated 6 SNPs from the HDAC1, MCM2, CDK1, BUB1B, CDC14B, and CCNE1 genes using 8 bioinformatics tools. Therefore, these SNPs were identified as potentially harmful by multiple computational tools. Specifically, rs199958833 in CDK1 (Val124Gly) was predicted to be damaging by all tools used. Our analysis strongly indicates that this specific SNP could negatively affect the stability and functionality of the CDK1 protein.

CONCLUSION

Based on our current understanding, the evaluation of CDK1 polymorphisms in the context of drug resistance in CRC has yet to be undertaken. In this investigation, we showed that rs199958833 variant in the CDK1 gene may favor resistance to 5-FU-based chemotherapy. However, these findings need validation in an independent cohort of patients.

Resource: A curated database of brain-related functional gene sets (Brain.GMT)

Transcriptional profiling has become a common tool for investigating the nervous system. During analysis, differential expression results are often compared to functional ontology databases, which contain curated gene sets representing well-studied pathways. This dependence can cause neuroscience studies to be interpreted in terms of functional pathways documented in better studied tissues (e.g., liver) and topics (e.g., cancer), and systematically emphasizes well-studied genes, leaving other findings in the obscurity of the brain "ignorome". To address this issue, we compiled a curated database of 918 gene sets related to nervous system function, tissue, and cell types ("Brain.GMT") that can be used within common analysis pipelines (GSEA, limma, edgeR) to interpret results from three species (rat, mouse, human). Brain.GMT includes brain-related gene sets curated from the Molecular Signatures Database (MSigDB) and extracted from public databases (GeneWeaver, Gemma, DropViz, BrainInABlender, HippoSeq) and published studies containing differential expression results. Although Brain.GMT is still undergoing development and currently only represents a fraction of available brain gene sets, "brain ignorome" genes are already better represented than in traditional Gene Ontology databases. Moreover, Brain.GMT substantially improves the quantity and quality of gene sets identified as enriched with differential expression in neuroscience studies, enhancing interpretation. •We compiled a curated database of 918 gene sets related to nervous system function, tissue, and cell types ("Brain.GMT").•Brain.GMT can be used within common analysis pipelines (GSEA, limma, edgeR) to interpret neuroscience transcriptional profiling results from three species (rat, mouse, human).•Although Brain.GMT is still undergoing development, it substantially improved the interpretation of differential expression results within our initial use cases.

An epigenetic human cytomegalovirus infection score predicts viremia risk in seropositive lung transplant recipients

Cytomegalovirus (CMV) infection and reactivation in solid organ transplant (SOT) recipients increases the risk of viremia, graft failure and death. Clinical studies of CMV serostatus indicate that donor positive recipient negative (D+/R-) patients have greater viremia risk than D-/R-. The majority of patients are R+ having intermediate serologic risk. To characterize the long-term impact of CMV infection and assess viremia risk, we sought to measure the effects of CMV on the recipient immune epigenome. Specifically, we profiled DNA methylation in 156 individuals before lung or kidney transplant. We found that the methylome of CMV positive SOT recipients is hyper-methylated at loci associated with neural development and Polycomb group (PcG) protein binding, and hypo-methylated at regions critical for the maturation of lymphocytes. In addition, we developed a machine learning-based model to predict the recipient CMV serostatus after correcting for cell type composition and ancestry. This CMV episcore measured at baseline in R+ individual stratifies viremia risk accurately in the lung transplant cohort, and along with serostatus the CMV episcore could be a potential biomarker for identifying R+ patients at high viremia risk.

Control of nuclear localization of the nucleocapsid protein of SARS-CoV-2

The nucleocapsid (N) protein of coronaviruses is a structural protein that binds viral RNA for assembly into the mature virion, a process that occurs in the cytoplasm. Several coronavirus N proteins also localize to the nucleus. Herein, we identify that two sequences (NLSs) are required for nuclear localization of the SARS-CoV-2 N protein. Deletion or mutation of these two sequences creates an N protein that does not localize to the nucleus in HEK293T cells. Overexpression of both wild-type and NLS-mutated N proteins dysregulate a largely overlapping set of mRNAs in HEK293T cells, suggesting that these N proteins do not have direct nuclear effects on transcription. Consistent with that hypothesis, both N proteins induce nuclear localization of NF-κB p65 and dysregulate a set of previously identified NF-κB-dependent genes. The effects of N on nuclear properties are proposed to alter host cell functions that contribute to viral pathogenesis or replication.

Correlating tissue and plasma‑specific piRNA changes to predict their possible role in pancreatic malignancy and chronic inflammation

The aggressiveness of pancreatic ductal adenocarcinoma is primarily due to lack of effective early detection biomarkers. Circulating non-coding RNAs serve as diagnostic or prognostic biomarkers in multiple types of cancer. Comparison of their expression between diseased tissue and relevant body fluids such as saliva, urine, bile, pancreatic juice, blood etc. may reveal mechanistic involvement of common non-coding RNAs. piwi-interacting RNAs (piRNAs) are a class of non-coding RNAs. The aim of the present study was to investigate plasma and tumour tissue piRNA changes in patients with pancreatic cancer (PC) and explore the possible role in tumorigenesis and pancreatic inflammation. Sequencing of circulating plasma small RNAs from patients with PC and chronic pancreatitis (CP) was performed and differentially expressed piRNAs were compared with those in tissues. Subsequent search for target genes for those piRNAs was performed followed by pathway and cluster analysis. A total of 36 piRNAs were shown to be deregulated in pancreatic tumour tissue and alteration of 11 piRNAs was detected in plasma of patients with PC. piRNAs hsa-piR-23246, hsa-piR-32858 and hsa-piR-9137 may serve a key role in PC development as their expression was correlated in both plasma and tumour tissue. Key piRNA-target interactions interfering with key biological pathways were also characterized. A total of 19 deregulated piRNAs in plasma samples of patients with CP was identified; these targeted genes responsible for chronic inflammation. Therefore, the present study provides a comprehensive description of piRNA alteration in pancreatic malignancy and inflammation; these may be explored for biomarker potential in future.

A Copernican revolution of multigenic analysis: A retrospective study on clinical exome sequencing in unclear genetic disorders

Despite the inevitable shift in medical practice towards a deeper understanding of disease etiology and progression through multigenic analysis, the profound historical impact of Mendelian diseases cannot be overlooked. These diseases, such as cystic fibrosis and thalassemia, are characterized by a single variant in a single gene leading to clinical conditions, and have significantly shaped our medical knowledge and treatments. In this respect, the monogenic approach inevitably results in the underutilization of Next-Generation Sequencing (NGS) data. Herein, a retrospective study was performed to assess the diagnostic value of the clinical exome in 32 probands with specific phenotypic characteristics (patients with autoinflammation and immunological dysregulation, N = 20; patients diagnosed with Hemolytic uremic syndrome N = 9; and patients with Waldenström macroglobulinemia, N = 3). A gene enrichment analysis was performed using the *. VCF file generated by SOPHiA-DDM-v4. This analysis selected a subset of genes containing pathogenic or likely pathogenic variants with autosomal dominant (AD) inheritance. In addition, all variants of uncertain significance (VUS) were included, filtered by AD inheritance mode, the presence of compound heterozygotes, and a minor allele frequency (MAF) cutoff of 0.05 %. The aim of the pipeline described here is based on a perspective shift that focuses on analyzing patients' gene assets, offering new light on the complex interplay between genetics and disease presentation. Integrating this approach into clinical practices could significantly enhance the management of patients with rare genetic disorders.

Comprehensive analysis of m6A methylome alterations after azacytidine plus venetoclax treatment for acute myeloid leukemia by nanopore sequencing

N6 adenosine methylation (m6A), one of the most prevalent internal modifications on mammalian RNAs, regulates RNA transcription, stabilization, and splicing. Growing evidence has focused on the functional role of m6A regulators on acute myeloid leukemia (AML). However, the global m6A levels after azacytidine (AZA) plus venetoclax (VEN) treatment in AML patients remain unclear. In our present study, bone marrow (BM) sample pairs (including pre-treatment [AML] and post-treatment [complete remission (CR)] samples) were harvested from three AML patients who had achieved CR after AZA plus VEN treatment for Nanopore direct RNA sequencing. Notably, the amount of m6A sites and the m6A levels in CR BMs was significantly lower than those in the AML BMs. Such a significant reduction in the m6A levels was also detected in AZA-treated HL-60 cells. Thirteen genes with decreased m6A and expression levels were identified, among which three genes (HPRT1, SNRPC, and ANP32B) were closely related to the prognosis of AML. Finally, we speculated the mechanism via which m6A modifications affected the mRNA stability of these three genes. In conclusion, we illustrated for the first time the global landscape of m6A levels in AZA plus VEN treated AML (CR) patients and revealed that AZA had a significant demethylation effect at the RNA level in AML patients. In addition, we identified new biomarkers for AZA plus VEN-treated AML via Nanopore sequencing technology in RNA epigenetics.

Maternal glucose levels and late pregnancy circulating extracellular vesicle and particle miRNAs in the MADRES pregnancy cohort

Maternal hyperglycemia during pregnancy adversely affects maternal and child outcomes. While mechanisms are not fully understood, maternal circulating miRNAs may play a role. We examined whether continuous glucose levels and hyperglycemia subtypes (gestational diabetes, type 2 diabetes, and glucose intolerance) were associated with circulating miRNAs during late pregnancy. Seven miRNAs (hsa-miR-107, hsa-let-7b-5p, hsa-miR-126-3p, hsa-miR-181a-5p, hsa-miR-374a-5p, hsa-miR-382-5p, and hsa-miR-337-5p) were associated (p < 0.05) with either hyperglycemia or continuous glucose levels prior to multiple testing correction. These miRNAs target genes involved in pathways relevant to maternal and child health, including insulin signaling, placental development, energy balance, and appetite regulation.

Chronic pain-induced methylation in the prefrontal cortex targets gene networks associated with cognition and Alzheimer's disease

Chronic pain is prevalent among aging adults. Epidemiologic evidence has demonstrated that individuals with chronic pain have accelerated memory decline and increased probability of dementia. Neurophysiologic, molecular, and pharmacologic hypotheses have been proposed to explain the relationship between chronic pain and cognitive decline, but there remains currently limited evidence supporting any of these. Here, we integrate multi-omic data across human cohorts and rodent species and demonstrate that methylation in the prefrontal cortex induced by chronic pain specifically targets transcriptional networks associated with cognitive ability, memory, and Alzheimer's disease in humans. We validate this with multiple independent data sets and identify cortical microglia as a likely mechanism by which chronic pain can increase dementia risk. Our analyses support the molecular hypothesis for the role of chronic pain in cognitive decline and identifies several potential therapeutic targets.

Impact of HDAC6-mediated progesterone receptor expression on the response of breast cancer cells to hormonal therapy

Modulation of estrogen receptor (ER) and progesterone receptor (PR) expression, as well as their emerging functional crosstalk, remains a potential approach for enhancing the response to hormonal therapy in breast cancer. Aberrant epigenetic alterations induced by histone deacetylases (HDACs) were massively implicated in dysregulating the function of hormone receptors in breast cancer. Although much is known about the regulation of ER signaling by HDAC, the precise role of HDAC in modulating the expression of PR and its impact on the outcomes of hormonal therapy is poorly defined. Here, we demonstrate the involvement of HDAC6 in regulating PR expression in breast cancer cells. The correlation between HDAC6 and hormone receptors was investigated in patients' tissues by immunohistochemistry (n = 80) and publicly available data (n = 3260) from breast cancer patients. We explored the effect of modulating the expression of HDAC6 as well as its catalytic inhibition on the level of hormone receptors by a variety of molecular analyses, including Western blot, immunofluorescence, Real-time PCR, RNA-seq analysis and chromatin immunoprecipitation. Based on our in-silico and immunohistochemistry analyses, HDAC6 levels were negatively correlated with PR status in breast cancer tissues. The downregulation of HDAC6 enhanced the expression of PR-B in hormone receptor-positive and triple-negative breast cancer (TNBC) cells. The selective targeting of HDAC6 by tubacin resulted in the enrichment of the H3K9 acetylation mark at the PGR-B gene promoter region and enhanced the expression of PR-B. Additionally, transcriptomic analysis of tubacin-treated cells revealed enhanced activity of acetyltransferase and growth factor signaling pathways, along with the enrichment of transcription factors involved in the transcriptional activity of ER, underscoring the crucial role of HDAC6 in regulating hormone receptors. Notably, the addition of HDAC6 inhibitor potentiated the effects of anti-ER and anti-PR drugs mainly in TNBC cells. Together, these data highlight the role of HDAC6 in regulating PR expression and provide a promising therapeutic approach for boosting breast cancer sensitivity to hormonal therapy.

Integrated multi-omics analysis of PBX1 in mouse adult neural stem- and progenitor cells identifies a transcriptional module that functionally links PBX1 to TCF3/4

Developmental transcription factors act in networks, but how these networks achieve cell- and tissue specificity is still poorly understood. Here, we explored pre-B cell leukemia homeobox 1 (PBX1) in adult neurogenesis combining genomic, transcriptomic, and proteomic approaches. ChIP-seq analysis uncovered PBX1 binding to numerous genomic sites. Integration of PBX1 ChIP-seq with ATAC-seq data predicted interaction partners, which were subsequently validated by mass spectrometry. Whole transcriptome spatial RNA analysis revealed shared expression dynamics of Pbx1 and interacting factors. Among these were class I bHLH proteins TCF3 and TCF4. RNA-seq following Pbx1, Tcf3 or Tcf4 knockdown identified proliferation- and differentiation associated genes as shared targets, while sphere formation assays following knockdown argued for functional cooperativity of PBX1 and TCF3 in progenitor cell proliferation. Notably, while physiological PBX1-TCF interaction has not yet been described, chromosomal translocation resulting in genomic TCF3::PBX1 fusion characterizes a subtype of acute lymphoblastic leukemia. Introducing Pbx1 into Nalm6 cells, a pre-B cell line expressing TCF3 but lacking PBX1, upregulated the leukemogenic genes BLK and NOTCH3, arguing that functional PBX1-TCF cooperativity likely extends to hematopoiesis. Our study hence uncovers a transcriptional module orchestrating the balance between progenitor cell proliferation and differentiation in adult neurogenesis with potential implications for leukemia etiology.

Extracellular vesicle transfer of miR-1 to adipose tissue modifies lipolytic pathways following resistance exercise

Extracellular vesicles (EVs) have emerged as important mediators of intertissue signaling and exercise adaptations. In this human study, we provide evidence that muscle-specific microRNA-1 (miR-1) was transferred to adipose tissue via EVs following an acute bout of resistance exercise. Using a multimodel machine learning automation tool, we discovered muscle primary miR-1 transcript and CD63+ EV count in circulation as top explanatory features for changes in adipose miR-1 levels in response to resistance exercise. RNA-Seq and in-silico prediction of miR-1 target genes identified caveolin 2 (CAV2) and tripartite motif containing 6 (TRIM6) as miR-1 target genes downregulated in the adipose tissue of a subset of participants with the highest increases in miR-1 levels following resistance exercise. Overexpression of miR-1 in differentiated human adipocyte-derived stem cells downregulated these miR-1 targets and enhanced catecholamine-induced lipolysis. These data identify a potential EV-mediated mechanism by which skeletal muscle communicates with adipose tissue and modulates lipolysis via miR-1.

Viscoelastic High-Molecular-Weight Hyaluronic Acid Hydrogels Support Rapid Glioblastoma Cell Invasion with Leader-Follower Dynamics

Hyaluronic acid (HA), the primary component of brain extracellular matrix, is increasingly used to model neuropathological processes, including glioblastoma (GBM) tumor invasion. While elastic hydrogels based on crosslinked low-molecular-weight (LMW) HA are widely exploited for this purpose and have proven valuable for discovery and screening, brain tissue is both viscoelastic and rich in high-MW (HMW) HA, and it remains unclear how these differences influence invasion. To address this question, hydrogels comprised of either HMW (1.5 MDa) or LMW (60 kDa) HA are introduced, characterized, and applied in GBM invasion studies. Unlike LMW HA hydrogels, HMW HA hydrogels relax stresses quickly, to a similar extent as brain tissue, and to a greater extent than many conventional HA-based scaffolds. GBM cells implanted within HMW HA hydrogels invade much more rapidly than in their LMW HA counterparts and exhibit distinct leader-follower dynamics. Leader cells adopt dendritic morphologies similar to invasive GBM cells observed in vivo. Transcriptomic, pharmacologic, and imaging studies suggest that leader cells exploit hyaluronidase, an enzyme strongly enriched in human GBMs, to prime a path for followers. This study offers new insight into how HA viscoelastic properties drive invasion and argues for the use of highly stress-relaxing materials to model GBM.

Polygenic and transcriptional risk scores identify chronic obstructive pulmonary disease subtypes in the COPDGene and ECLIPSE cohort studies

BACKGROUND

Genetic variants and gene expression predict risk of chronic obstructive pulmonary disease (COPD), but their effect on COPD heterogeneity is unclear. We aimed to define high-risk COPD subtypes using genetics (polygenic risk score, PRS) and blood gene expression (transcriptional risk score, TRS) and assess differences in clinical and molecular characteristics.

METHODS

We defined high-risk groups based on PRS and TRS quantiles by maximising differences in protein biomarkers in a COPDGene training set and identified these groups in COPDGene and ECLIPSE test sets. We tested multivariable associations of subgroups with clinical outcomes and compared protein-protein interaction networks and drug repurposing analyses between high-risk groups.

FINDINGS

We examined two high-risk omics-defined groups in non-overlapping test sets (n = 1133 NHW COPDGene, n = 299 African American (AA) COPDGene, n = 468 ECLIPSE). We defined "high activity" (low PRS, high TRS) and "severe risk" (high PRS, high TRS) subgroups. Participants in both subgroups had lower body-mass index (BMI), lower lung function, and alterations in metabolic, growth, and immune signalling processes compared to a low-risk (low PRS, low TRS) subgroup. "High activity" but not "severe risk" participants had greater prospective FEV1 decline (COPDGene: -51 mL/year; ECLIPSE: -40 mL/year) and proteomic profiles were enriched in gene sets perturbed by treatment with 5-lipoxygenase inhibitors and angiotensin-converting enzyme (ACE) inhibitors.

INTERPRETATION

Concomitant use of polygenic and transcriptional risk scores identified clinical and molecular heterogeneity amongst high-risk individuals. Proteomic and drug repurposing analysis identified subtype-specific enrichment for therapies and suggest prior drug repurposing failures may be explained by patient selection.

FUNDING

National Institutes of Health.

Therapeutic CRISPR epigenome editing of inflammatory receptors in the intervertebral disc

Low back pain (LBP) ranks among the leading causes of disability worldwide and generates a tremendous socioeconomic cost. Disc degeneration, a leading contributor to LBP, can be characterized by the breakdown of the extracellular matrix of the intervertebral disc (IVD), disc height loss, and inflammation. The inflammatory cytokine tumor necrosis factor α (TNF-α) has multiple signaling pathways, including proinflammatory signaling through tumor necrosis factor receptor 1 superfamily, member 1a (TNFR1 or TNFRSF1A), and has been implicated as a primary mediator of disc degeneration. We tested our ability to regulate the TNFR1 signaling pathway in vivo, utilizing CRISPR epigenome editing to slow the progression of disc degeneration in rats. Sprague-Dawley rats were treated with TNF-α and CRISPR interference (CRISPRi)-based epigenome-editing therapeutics targeting TNFR1, showing decreased behavioral pain in a disc degeneration model. Surprisingly, while treatment with the vectors alone was therapeutic, the TNF-α injection became therapeutic after TNFR1 modulation. These results suggest direct inflammatory receptor modulation as a potent strategy for treating disc degeneration.

How the liver transcriptome and lipid composition influence the progression of nonalcoholic fatty liver disease to hepatocellular carcinoma in a murine model

The incidence of nonalcoholic fatty liver disease (NAFLD) has been steadily increasing in Western society in recent years and has been recognized as a risk factor for the development of hepatocellular carcinoma (HCC). However, the molecular mechanisms underlying the progression from NAFLD to HCC are still unclear, despite the use of suitable mouse models. To identify the transcriptional and lipid profiles of livers from mice with NAFLD-HCC, we induced both NAFLD and NAFLD-HCC pathologies in C57BL/6J mice and performed RNA-sequencing (RNA-seq) and targeted lipidomic analysis. Our RNA-seq analysis revealed that the transcriptional signature of NAFLD in mice is characterized by changes in inflammatory response and fatty acid metabolism. Moreover, the signature of NAFLD-HCC is characterized by processes typically observed in cancer, such as epithelial to mesenchymal transition, angiogenesis and inflammatory responses. Furthermore, we found that the diet used in this study inhibited cholesterol synthesis in both models. The analysis of lipid composition also showed a significant impact of the provided diet. Therefore, our study supports the idea that a Western diet (WD) affects metabolic processes and hepatic lipid composition. Additionally, the combination of a WD with the administration of a carcinogen drives the progression from NAFLD to HCC.

Network medicine informed multiomics integration identifies drug targets and repurposable medicines for Amyotrophic Lateral Sclerosis

Amyotrophic Lateral Sclerosis (ALS) is a devastating, immensely complex neurodegenerative disease by lack of effective treatments. We developed a network medicine methodology via integrating human brain multi-omics data to prioritize drug targets and repurposable treatments for ALS. We leveraged non-coding ALS loci effects from genome-wide associated studies (GWAS) on human brain expression quantitative trait loci (QTL) (eQTL), protein QTL (pQTL), splicing QTL (sQTL), methylation QTL (meQTL), and histone acetylation QTL (haQTL). Using a network-based deep learning framework, we identified 105 putative ALS-associated genes enriched in known ALS pathobiological pathways. Applying network proximity analysis of predicted ALS-associated genes and drug-target networks under the human protein-protein interactome (PPI) model, we identified potential repurposable drugs (i.e., Diazoxide and Gefitinib) for ALS. Subsequent validation established preclinical evidence for top-prioritized drugs. In summary, we presented a network-based multi-omics framework to identify drug targets and repurposable treatments for ALS and other neurodegenerative disease if broadly applied.

Olfr2-positive macrophages originate from monocytes proliferate in situ and present a pro-inflammatory foamy-like phenotype

AIMS

Olfactory receptor 2 (Olfr2) has been identified in a minimum of 30% of vascular macrophages, and its depletion was shown to reduce atherosclerosis progression. Mononuclear phagocytes, including monocytes and macrophages within the vessel wall, are major players in atherosclerosis. Single-cell RNA sequencing studies revealed that atherosclerotic artery walls encompass several monocytes and vascular macrophages, defining at least nine distinct subsets potentially serving diverse functions in disease progression. This study investigates the functional phenotype and ontogeny of Olfr2-expressing vascular macrophages in atherosclerosis.

METHODS AND RESULTS

Olfr2+ macrophages rapidly increase in Apoe-/- mice's aorta when fed a Western diet (WD). Mass cytometry showed that Olfr2+ cells are clustered within the CD64 high population and enriched for CD11c and Ccr2 markers. Olfr2+ macrophages express many pro-inflammatory cytokines, including Il1b, Il6, Il12, and Il23, and chemokines, including Ccl5, Cx3cl1, Cxcl9, and Ccl22. By extracting differentially expressed genes from bulk RNA sequencing (RNA-seq) of Olfr2+ vs. Olfr2- macrophages, we defined a signature that significantly mapped to single-cell data of plaque myeloid cells, including monocytes, subendothelial MacAir, and Trem2Gpnmb foamy macrophages. By adoptive transfer experiments, we identified that Olfr2 competent monocytes from CD45.1Apoe-/-Olfr2+/+ mice transferred into CD45.2Apoe-/-Olfr2-/- recipient mice fed WD for 12 weeks, accumulate in the atherosclerotic aorta wall already at 72 h, and differentiate in macrophages. Olfr2+ macrophages showed significantly increased BrdU incorporation compared to Olfr2- macrophages. Flow cytometry confirmed that at least 50% of aortic Olfr2+ macrophages are positive for BODIPY staining and have increased expression of both tumour necrosis factor and interleukin 6 compared to Olfr2- macrophages. Gene set enrichment analysis of the Olfr2+ macrophage signature revealed a similar enrichment pattern in human atherosclerotic plaques, particularly within foamy/TREM2hi-Mφ and monocytes.

CONCLUSIONS

In summary, we conclude that Olfr2+ macrophages in the aorta originate from monocytes and can accumulate at the early stages of disease progression. These cells can undergo differentiation into MacAir and Trem2Gpnmb foamy macrophages, exhibiting proliferative and pro-inflammatory potentials. This dynamic behaviour positions them as key influencers in shaping the myeloid landscape within the atherosclerotic plaque.

Comparative single-cell analyses identify shared and divergent features of human and mouse kidney development

The mammalian kidney maintains fluid homeostasis through diverse epithelial cell types generated from nephron and ureteric progenitor cells. To extend a developmental understanding of the kidney's epithelial networks, we compared chromatin organization (single-nuclear assay for transposase-accessible chromatin sequencing [ATAC-seq]; 112,864 nuclei) and gene expression (single-cell/nuclear RNA sequencing [RNA-seq]; 109,477 cells/nuclei) in the developing human (10.6-17.6 weeks; n = 10) and mouse (post-natal day [P]0; n = 10) kidney, supplementing analysis with published mouse datasets from earlier stages. Single-cell/nuclear datasets were analyzed at a species level, and then nephron and ureteric cellular lineages were extracted and integrated into a common, cross-species, multimodal dataset. Comparative computational analyses identified conserved and divergent features of chromatin organization and linked gene activity, identifying species-specific and cell-type-specific regulatory programs. In situ validation of human-enriched gene activity points to human-specific signaling interactions in kidney development. Further, human-specific enhancer regions were linked to kidney diseases through genome-wide association studies (GWASs), highlighting the potential for clinical insight from developmental modeling.

Breast cancer genomic analyses reveal genes, mutations, and signaling networks

Breast cancer (BC) is the most commonly diagnosed cancer and the predominant cause of death in women. BC is a complex disorder, and the exploration of several types of BC omic data, highlighting genes, perturbations, signaling and cellular mechanisms, is needed. We collected mutational data from 9,555 BC samples using cBioPortal. We classified 1174 BC genes (mutated ≥ 40 samples) into five tiers (BCtier_I-V) and subjected them to pathway and protein‒protein network analyses using EnrichR and STRING 11, respectively. BCtier_I possesses 12 BC genes with mutational frequencies > 5%, with only 5 genes possessing > 10% frequencies, namely, PIK3CA (35.7%), TP53 (34.3%), GATA3 (11.5%), CDH1 (11.4%) and MUC16 (11%), and the next seven BC genes are KMT2C (8.8%), TTN (8%), MAP3K1 (8%), SYNE1 (7.2%), AHNAK2 (7%), USH2A (5.5%), and RYR2 (5.4%). Our pathway analyses revealed that the five top BC pathways were the PI3K-AKT, TP53, NOTCH, HIPPO, and RAS pathways. We found that BC panels share only seven genes. These findings show that BC arises from genetic disruptions evident in BC signaling and protein networks.

Hox gene activity directs physical forces to differentially shape chick small and large intestinal epithelia

Hox transcription factors play crucial roles in organizing developmental patterning across metazoa, but how these factors trigger regional morphogenesis has largely remained a mystery. In the developing gut, Hox genes help demarcate identities of intestinal subregions early in embryogenesis, which ultimately leads to their specialization in both form and function. Although the midgut forms villi, the hindgut develops sulci that resolve into heterogeneous outgrowths. Combining mechanical measurements of the embryonic chick intestine and mathematical modeling, we demonstrate that the posterior Hox gene HOXD13 regulates biophysical phenomena that shape the hindgut lumen. We further show that HOXD13 acts through the transforming growth factor β (TGF-β) pathway to thicken, stiffen, and promote isotropic growth of the subepithelial mesenchyme-together, these features lead to hindgut-specific surface buckling. TGF-β, in turn, promotes collagen deposition to affect mesenchymal geometry and growth. We thus identify a cascade of events downstream of positional identity that direct posterior intestinal morphogenesis.

Cancer-associated fibroblasts shape early myeloid cell response to chemotherapy-induced immunogenic signals in next generation tumor organoid cultures

BACKGROUND

Patient-derived colorectal cancer (CRC) organoids (PDOs) solely consisting of malignant cells led to major advances in the understanding of cancer treatments. Yet, a major limitation is the absence of cells from the tumor microenvironment, thereby prohibiting potential investigation of treatment responses on immune and structural cells. Currently there are sparse reports describing the interaction of PDOs, cancer-associated fibroblasts (CAFs) and tumor-associated macrophages (TAMs) in complex primary co-culture assay systems.

METHODS

Primary PDOs and patient matched CAF cultures were generated from surgical resections. Co-culture systems of PDOs, CAFs and monocytic myeloid cells were set up to recapitulate features seen in patient tumors. Single-cell transcriptomics and flow cytometry was used to show effects of culture systems on TAM populations in the co-culture assays under chemotherapeutic and oncolytic viral treatment.

RESULTS

In contrast to co-cultures of tumor cells and monocytes, CAF/monocyte co-cultures and CAF/monocyte/tumor cell triple cultures resulted in a partial differentiation into macrophages and a phenotypic switch, characterized by the expression of major immunosuppressive markers comparable to TAMs in CRC. Oxaliplatin and 5-fluorouracil, the standard-of-care chemotherapy for CRC, induced polarization of macrophages to a pro-inflammatory phenotype comparable to the immunogenic effects of treatment with an oncolytic virus. Monitoring phagocytosis as a functional proxy to macrophage activation and subsequent onset of an immune response, revealed that chemotherapy-induced cell death, but not virus-mediated cell death, is necessary to induce phagocytosis of CRC cells. Moreover, CAFs enhanced the phagocytic activity in chemotherapy treated CRC triple cultures.

CONCLUSIONS

Primary CAF-containing triple cultures successfully model TAM-like phenotypes ex vivo and allow the assessment of their functional and phenotypic changes in response to treatments following a precision medicine approach.

Dermal TRPV1 innervations engage a macrophage- and fibroblast-containing pathway to activate hair growth in mice

Pain, detected by nociceptors, is an integral part of injury, yet whether and how it can impact tissue physiology and recovery remain understudied. Here, we applied chemogenetics in mice to locally activate dermal TRPV1 innervations in naive skin and found that it triggered new regenerative cycling by dormant hair follicles (HFs). This was preceded by rapid apoptosis of dermal macrophages, mediated by the neuropeptide calcitonin gene-related peptide (CGRP). TRPV1 activation also triggered a macrophage-dependent induction of osteopontin (Spp1)-expressing dermal fibroblasts. The neuropeptide CGRP and the extracellular matrix protein Spp1 were required for the nociceptor-triggered hair growth. Finally, we showed that epidermal abrasion injury induced Spp1-expressing dermal fibroblasts and hair growth via a TRPV1 neuron and CGRP-dependent mechanism. Collectively, these data demonstrated a role for TRPV1 nociceptors in orchestrating a macrophage and fibroblast-supported mechanism to promote hair growth and enabling the efficient restoration of this mechano- and thermo-protective barrier after wounding.

MiR-107 and Its Association With House Dust Mite Sensitisation: Implications for Asthma

INTRODUCTION

MicroRNAs (miRNAs) have been linked to allergic diseases but their effects on sensitisation to allergens in individuals with asthma are unknown. We aimed to identify miRNAs associated with house dust mite (HDM) sensitisation in childhood asthma.

METHODS

Serum samples from 1126 children with asthma who participated in the Genetics of Asthma in Costa Rica Study (GACRS) were profiled for 304 miRNAs. We first divided according to HDM sensitisation and then tested whether miRNAs were differentially expressed (DE) between the two groups. Gene enrichment analysis for target genes of the DE miRNAs was then performed to identify potential causal pathways. Replication analysis was performed in the Childhood Asthma Management Program (CAMP), in which expression data of 258 miRNAs in 491 children were available. A mediation analysis was conducted to discern relationships between miRNA and phenotype differences according to HDM sensitisation in GACRS cohort.

RESULTS

There were 906 (80.5%) and 220 (19.5%) subjects in the GACRS HDM+ and HDM- groups. Compared with HDM- participants, those in the HDM+ group were more likely to be severe in variables including pulmonary function, oral corticosteroid use and blood tests. A total of 17 miRNAs were DE (p < 0.05) between the two groups, with miR-642a-3p, let-7c-5p and miR-107 most significantly associated with HDM sensitisation. In CAMP, there were 39 DE miRNAs, and increased expression of miR-107 in HDM+ children was replicated in this cohort. In both GACRS and CAMP, the cadherin-binding pathway was enriched in an analysis of target genes for DE miRNA. In a mediation analysis, miR-107 showed significant indirect effects on eosinophil count and total IgE that were mediated by HDM sensitisation.

CONCLUSION

In children with asthma, miR-107 is associated with HDM sensitisation. Furthermore, miR-107 was indirectly associated with total IgE and eosinophil count through HDM sensitisation.

TFEB controls syncytiotrophoblast formation and hormone production in placenta

TFEB, a bHLH-leucine zipper transcription factor belonging to the MiT/TFE family, globally modulates cell metabolism by regulating autophagy and lysosomal functions. Remarkably, loss of TFEB in mice causes embryonic lethality due to severe defects in placentation associated with aberrant vascularization and resulting hypoxia. However, the molecular mechanism underlying this phenotype has remained elusive. By integrating in vivo analyses with multi-omics approaches and functional assays, we have uncovered an unprecedented function for TFEB in promoting the formation of a functional syncytiotrophoblast in the placenta. Our findings demonstrate that constitutive loss of TFEB in knock-out mice is associated with defective formation of the syncytiotrophoblast layer. Indeed, using in vitro models of syncytialization, we demonstrated that TFEB translocates into the nucleus during syncytiotrophoblast formation and binds to the promoters of crucial placental genes, including genes encoding fusogenic proteins (Syncytin-1 and Syncytin-2) and enzymes involved in steroidogenic pathways, such as CYP19A1, the rate-limiting enzyme for the synthesis of 17β-Estradiol (E2). Conversely, TFEB depletion impairs both syncytial fusion and endocrine properties of syncytiotrophoblast, as demonstrated by a significant decrease in the secretion of placental hormones and E2 production. Notably, restoration of TFEB expression resets syncytiotrophoblast identity. Our findings identify that TFEB controls placental development and function by orchestrating both the transcriptional program underlying trophoblast fusion and the acquisition of endocrine function, which are crucial for the bioenergetic requirements of embryonic development.

Genomic Analysis Identifies Risk Factors in Restless Legs Syndrome

OBJECTIVE

Restless legs syndrome (RLS) is a neurological condition that causes uncomfortable sensations in the legs and an irresistible urge to move them, typically during periods of rest. The genetic basis and pathophysiology of RLS are incompletely understood. We sought to identify additional novel genetic risk factors associated with RLS susceptibility.

METHODS

We performed a whole-genome sequencing and genome-wide association meta-analysis of RLS cases (n = 9,851) and controls (n = 38,957) in 3 population-based biobanks (All of Us, Canadian Longitudinal Study on Aging, and CARTaGENE).

RESULTS

Genome-wide association analysis identified 9 independent risk loci, of which 8 had been previously reported, and 1 was a novel risk locus (LMX1B, rs35196838, OR 1.14, 95% CI 1.09-1.19, p value = 2.2 × 10-9). Furthermore, a transcriptome-wide association study also identified GLO1 and a previously unreported gene, ELFN1. A genetic correlation analysis revealed significant common variant overlaps between RLS and neuroticism (rg = 0.40, se = 0.08, p value = 5.4 × 10-7), depression (rg = 0.35, se = 0.06, p value = 2.17 × 10-8), and intelligence (rg = -0.20, se = 0.06, p value = 4.0 × 10-4).

INTERPRETATION

Our study expands the understanding of the genetic architecture of RLS, and highlights the contributions of common variants to this prevalent neurological disorder. ANN NEUROL 2024;96:994-1005.

The autophagy protein RUBCNL/PACER represses RIPK1 kinase-dependent apoptosis and necroptosis

Mesenchymal stem cells (MSCs) are used in cell therapy; nonetheless, their application is limited by their poor survival after transplantation in a proinflammatory microenvironment. Macroautophagy/autophagy activation in MSCs constitutes a stress adaptation pathway, promoting cellular homeostasis. Our proteomics data indicate that RUBCNL/PACER (RUN and cysteine rich domain containing beclin 1 interacting protein like), a positive regulator of autophagy, is also involved in cell death. Hence, we screened MSC survival upon various cell death stimuli under loss or gain of function of RUBCNL. MSCs were protected from TNF (tumor necrosis factor)-induced regulated cell death when RUBCNL was expressed. TNF promotes inflammation by inducing RIPK1 kinase-dependent apoptosis or necroptosis. We determine that MSCs succumb to RIPK1 kinase-dependent apoptosis upon TNF sensing and necroptosis when caspases are inactivated. We show that RUBCNL is a negative regulator of both RIPK1-dependent apoptosis and necroptosis. Furthermore, RUBCNL mutants that lose the ability to regulate autophagy, retain their function in negatively regulating cell death. We also found that RUBCNL forms a complex with RIPK1, which disassembles in response to TNF. In line with this finding, RUBCNL expression limits assembly of RIPK1-TNFRSF1A/TNFR1 complex I, suggesting that complex formation between RUBCNL and RIPK1 represses TNF signaling. These results provide new insights into the crosstalk between the RIPK1-mediated cell death and autophagy machineries and suggest that RUBCNL, due to its functional duality in autophagy and apoptosis/necroptosis, could be targeted to improve the therapeutic efficacy of MSCs. Abbreviations: BAF: bafilomycin A1; CASP3: caspase 3; Caspases: cysteine-aspartic proteases; cCASP3: cleaved CASP3; CQ: chloroquine; CHX: cycloheximide; cPARP: cleaved poly (ADP-ribose) polymerase; DEPs: differential expressed proteins; ETO: etoposide; MEF: mouse embryonic fibroblast; MLKL: mixed lineage kinase domain-like; MSC: mesenchymal stem cell; MTORC1: mechanistic target of rapamycin kinase complex 1; Nec1s: necrostatin 1s; NFKB/NF-kB: nuclear factor of kappa light polypeptide gene enhancer in B cells; PLA: proximity ligation assay; RCD: regulated cell death; RIPK1: receptor (TNFRSF)-interacting serine-threonine kinase 1; RIPK3: receptor-interacting serine-threonine kinase 3; RUBCNL/PACER: RUN and cysteine rich domain containing beclin 1 interacting protein like; siCtrl: small interfering RNA nonsense; siRNA: small interfering RNA; TdT: terminal deoxynucleotidyl transferase; Tm: tunicamycin; TNF: tumor necrosis factor; TNFRSF1A/TNFR1: tumor necrosis factor receptor superfamily, member 1a.

Colony-stimulating factor 2 (CSF2) as a gut microbiome dependent immune factor that alters molecular and behavioral responses to cocaine in male mice

Cocaine use disorder is a condition that leads to tremendous morbidity and mortality for which there are currently no FDA-approved pharmacotherapies. Previous research has demonstrated an important role for the resident population of bacteria of the large intestine, collectively dubbed the gut microbiome, in modulating brain and behavior in models of cocaine and other substance use disorders. Importantly, previous work has repeatedly shown that depletion of the gut microbiome leads to increased cocaine taking and seeking behaviors in multiple models. While the precise mechanism of these gut-brain signaling pathways in models of cocaine use is not fully clear, and intriguing possibility is through gut microbiome influences on innate immune system function. In this manuscript we identify the cytokine colony stimulating factor 2 (CSF2) as an immune factor that is increased by cocaine in a gut microbiome dependent manner. Peripherally injected CSF2 crosses the blood-brain barrier into the nucleus accumbens, a brain region central to behavioral responses to cocaine. Treatment with peripheral CSF2 reduces acute and sensitized locomotor responses to cocaine as well as reducing cocaine place preference at high doses. On a molecular level, we find that peripheral injections of CSF2 alter the transcriptional response to both acute and repeated cocaine in the nucleus accumbens. Finally, treatment of microbiome depleted mice with CSF2 reverses the behavioral effects of microbiome depletion on the conditioned place preference assay. Taken together, this work identifies an innate immune factor that represents a novel gut-brain signaling cascade in models of cocaine use and lays the foundations for further translational work targeting this pathway.

A multi-ancestry cerebral cortex transcriptome-wide association study identifies genes associated with smoking behaviors

Transcriptome-wide association studies (TWAS) have provided valuable insight in identifying genes that may impact cigarette smoking. Most of previous studies, however, mainly focused on European ancestry. Limited TWAS studies have been conducted across multiple ancestries to explore genes that may impact smoking behaviors. In this study, we used cis-eQTL data of cerebral cortex from multiple ancestries in MetaBrain, including European, East Asian, and African samples, as reference panels to perform multi-ancestry TWAS analyses on ancestry-matched GWASs of four smoking behaviors including smoking initiation, smoking cessation, age of smoking initiation, and number of cigarettes per day in GWAS & Sequencing Consortium of Alcohol and Nicotine use (GSCAN). Multiple-ancestry fine-mapping approach was conducted to identify credible gene sets associated with these four traits. Enrichment and module network analyses were further performed to explore the potential roles of these identified gene sets. A total of 719 unique genes were identified to be associated with at least one of the four smoking traits across ancestries. Among those, 249 genes were further prioritized as putative causal genes in multiple ancestry-based fine-mapping approach. Several well-known smoking-related genes, including PSMA4, IREB2, and CHRNA3, showed high confidence across ancestries. Some novel genes, e.g., TSPAN3 and ANK2, were also identified in the credible sets. The enrichment analysis identified a series of critical pathways related to smoking such as synaptic transmission and glutamate receptor activity. Leveraging the power of the latest multi-ancestry GWAS and eQTL data sources, this study revealed hundreds of genes and relevant biological processes related to smoking behaviors. These findings provide new insights for future functional studies on smoking behaviors.

Neuroprotective effects of G9a inhibition through modulation of peroxisome-proliferator activator receptor gamma-dependent pathways by miR-128

JOURNAL/nrgr/04.03/01300535-202419110-00033/figure1/v/2024-03-08T184507Z/r/image-tiff Dysregulation of G9a, a histone-lysine N-methyltransferase, has been observed in Alzheimer's disease and has been correlated with increased levels of chronic inflammation and oxidative stress. Likewise, microRNAs are involved in many biological processes and diseases playing a key role in pathogenesis, especially in multifactorial diseases such as Alzheimer's disease. Therefore, our aim has been to provide partial insights into the interconnection between G9a, microRNAs, oxidative stress, and neuroinflammation. To better understand the biology of G9a, we compared the global microRNA expression between senescence-accelerated mouse-prone 8 (SAMP8) control mice and SAMP8 treated with G9a inhibitor UNC0642. We found a downregulation of miR-128 after a G9a inhibition treatment, which interestingly binds to the 3' untranslated region (3'-UTR) of peroxisome-proliferator activator receptor γ (PPARG) mRNA. Accordingly, Pparg gene expression levels were higher in the SAMP8 group treated with G9a inhibitor than in the SAMP8 control group. We also observed modulation of oxidative stress responses might be mainly driven Pparg after G9a inhibitor. To confirm these antioxidant effects, we treated primary neuron cell cultures with hydrogen peroxide as an oxidative insult. In this setting, treatment with G9a inhibitor increases both cell survival and antioxidant enzymes. Moreover, up-regulation of PPARγ by G9a inhibitor could also increase the expression of genes involved in DNA damage responses and apoptosis. In addition, we also described that the PPARγ/AMPK axis partially explains the regulation of autophagy markers expression. Finally, PPARγ/GADD45α potentially contributes to enhancing synaptic plasticity and neurogenesis after G9a inhibition. Altogether, we propose that pharmacological inhibition of G9a leads to a neuroprotective effect that could be due, at least in part, by the modulation of PPARγ-dependent pathways by miR-128.

Evaluation of cerebrospinal fluid (CSF) and interstitial fluid (ISF) mouse proteomes for the validation and description of Alzheimer's disease biomarkers

BACKGROUND

Mass spectrometry (MS)-based cerebrospinal fluid (CSF) proteomics is an important method for discovering biomarkers of neurodegenerative diseases. CSF serves as a reservoir for interstitial fluid (ISF), and extensive communication between the two fluid compartments helps to remove waste products from the brain.

NEW METHOD

We performed proteomic analyses of both CSF and ISF fluid compartments using intracerebral microdialysis to validate and detect novel biomarkers of Alzheimer's disease (AD) in APPtg and C57Bl/6J control mice.

RESULTS

We identified up to 625 proteins in ISF and 4483 proteins in CSF samples. By comparing the biofluid profiles of APPtg and C57Bl/6J mice, we detected 37 and 108 significantly up- and downregulated candidates, respectively. In ISF, 7 highly regulated proteins, such as Gfap, Aldh1l1, Gstm1, and Txn, have already been implicated in AD progression, whereas in CSF, 9 out of 14 highly regulated proteins, such as Apba2, Syt12, Pgs1 and Vsnl1, have also been validated to be involved in AD pathogenesis. In addition, we also detected new interesting regulated proteins related to the control of synapses and neurotransmission (Kcna2, Cacng3, and Clcn6) whose roles as AD biomarkers should be further investigated.

COMPARISON WITH EXISTING METHODS

This newly established combined protocol provides better insight into the mutual communication between ISF and CSF as an analysis of tissue or CSF compartments alone.

CONCLUSIONS

The use of multiple fluid compartments, ISF and CSF, for the detection of their biological communication enables better detection of new promising AD biomarkers.

A dual role for PSIP1/LEDGF in T cell acute lymphoblastic leukemia

T cell acute lymphoblastic leukemia (T-ALL) is an aggressive hematological malignancy. Current intensified therapeutic protocols coincide with severe side effects, and no salvage therapy is available for primary therapy-resistant or relapsed patients. This highlights the need to identify new therapeutic targets in T-ALL. PSIP1, dispensable for normal hematopoiesis, is a dependency factor in KMT2A-rearranged myeloid leukemia. Nonetheless, loss-of-function mutations suggest a tumor suppressor role for PSIP1 in T-ALL. Here, we demonstrate that the loss of Psip1 accelerates T-ALL initiation in mice which we correlated with reduced H3K27me3 binding. Contrastingly, loss of PSIP1 impaired cell proliferation in several T-ALL cell lines. In cell lines, PSIP1 down-regulation leads to a reduction of COX20, an assembly factor of the cytochrome c oxidase in the mitochondria, and to a reduction in mitochondrial respiration. This indicates that PSIP1 can exert a dual role in the context of T-ALL, either as a tumor suppressor gene during tumor initiation or as a dependency factor in tumor maintenance.

Neuropeptide signalling orchestrates T cell differentiation

The balance between T helper type 1 (TH1) cells and other TH cells is critical for antiviral and anti-tumour responses1-3, but how this balance is achieved remains poorly understood. Here we dissected the dynamic regulation of TH1 cell differentiation during in vitro polarization, and during in vivo differentiation after acute viral infection. We identified regulators modulating T helper cell differentiation using a unique TH1-TH2 cell dichotomous culture system and systematically validated their regulatory functions through multiple in vitro and in vivo CRISPR screens. We found that RAMP3, a component of the receptor for the neuropeptide CGRP (calcitonin gene-related peptide), has a cell-intrinsic role in TH1 cell fate determination. Extracellular CGRP signalling through the receptor RAMP3-CALCRL restricted the differentiation of TH2 cells, but promoted TH1 cell differentiation through the activation of downstream cAMP response element-binding protein (CREB) and activating transcription factor 3 (ATF3). ATF3 promoted TH1 cell differentiation by inducing the expression of Stat1, a key regulator of TH1 cell differentiation. After viral infection, an interaction between CGRP produced by neurons and RAMP3 expressed on T cells enhanced the anti-viral IFNγ-producing TH1 and CD8+ T cell response, and timely control of acute viral infection. Our research identifies a neuroimmune circuit in which neurons participate in T cell fate determination by producing the neuropeptide CGRP during acute viral infection, which acts on RAMP3-expressing T cells to induce an effective anti-viral TH1 cell response.

The TβRI promotes migration and metastasis through thrombospondin 1 and ITGAV in prostate cancer cells

TGFβ potently modifies the extracellular matrix (ECM), which is thought to favor tumor cell invasion. However, the mechanism whereby the cancer cells employ the ECM proteins to facilitate their motility is largely unknown. In this study we used RNA-seq and proteomic analysis to examine the proteins secreted by castration-resistant prostate cancer (CRPC) cells upon TGFβ treatment and found that thrombospondin 1 (THBS1) was observed to be one of the predominant proteins. The CRISPR Cas9, or siRNA techniques was used to downregulate TGFβ type I receptor (TβRI) to interfere with TGFβ signaling in various cancer cells in vitro. The interaction of ECM proteins with the TβRI in the migratory prostate cancer cells in response to TGFβ1 was demonstrated by several different techniques to reveal that THBS1 mediates cell migration by interacting with integrin subunit alpha V (ITGAV) and TβRI. Deletion of TβRI or THBS1 in cancer cells prevented their migration and invasion. THBS1 belongs to a group of tumorigenic ECM proteins induced via TGFβ signaling in CRPC cells, and high expression of THBS1 in human prostate cancer tissues correlated with the degree of malignancy. TGFβ-induced production of THBS1 through TβRI facilitates the invasion and metastasis of CRPC cells as shown in vivo xenograft animal experiments.

The interplay of mutagenesis and ecDNA shapes urothelial cancer evolution

Advanced urothelial cancer is a frequently lethal disease characterized by marked genetic heterogeneity1. In this study, we investigated the evolution of genomic signatures caused by endogenous and external mutagenic processes and their interplay with complex structural variants (SVs). We superimposed mutational signatures and phylogenetic analyses of matched serial tumours from patients with urothelial cancer to define the evolutionary dynamics of these processes. We show that APOBEC3-induced mutations are clonal and early, whereas chemotherapy induces mutational bursts of hundreds of late subclonal mutations. Using a genome graph computational tool2, we observed frequent high copy-number circular amplicons characteristic of extrachromosomal DNA (ecDNA)-forming SVs. We characterized the distinct temporal patterns of APOBEC3-induced and chemotherapy-induced mutations within ecDNA-forming SVs, gaining new insights into the timing of these mutagenic processes relative to ecDNA biogenesis. We discovered that most CCND1 amplifications in urothelial cancer arise within circular ecDNA-forming SVs. ecDNA-forming SVs persisted and increased in complexity, incorporating additional DNA segments and contributing to the evolution of treatment resistance. Oxford Nanopore Technologies long-read whole-genome sequencing followed by de novo assembly mapped out CCND1 ecDNA structure. Experimental modelling of CCND1 ecDNA confirmed its role as a driver of treatment resistance. Our findings define fundamental mechanisms that drive urothelial cancer evolution and have important therapeutic implications.

Advances in omics data for eosinophilic esophagitis: moving towards multi-omics analyses

Eosinophilic esophagitis (EoE) is a chronic, allergic inflammatory disease of the esophagus characterized by eosinophil accumulation and has a growing global prevalence. EoE significantly impairs quality of life and poses a substantial burden on healthcare resources. Currently, only two FDA-approved medications exist for EoE, highlighting the need for broader research into its management and prevention. Recent advancements in omics technologies, such as genomics, epigenetics, transcriptomics, proteomics, and others, offer new insights into the genetic and immunologic mechanisms underlying EoE. Genomic studies have identified genetic loci and mutations associated with EoE, revealing predispositions that vary by ancestry and indicating EoE's complex genetic basis. Epigenetic studies have uncovered changes in DNA methylation and chromatin structure that affect gene expression, influencing EoE pathology. Transcriptomic analyses have revealed a distinct gene expression profile in EoE, dominated by genes involved in activated type 2 immunity and epithelial barrier function. Proteomic approaches have furthered the understanding of EoE mechanisms, identifying potential new biomarkers and therapeutic targets. However, challenges in integrating diverse omics data persist, largely due to their complexity and the need for advanced computational methods. Machine learning is emerging as a valuable tool for analyzing extensive and intricate datasets, potentially revealing new aspects of EoE pathogenesis. The integration of multi-omics data through sophisticated computational approaches promises significant advancements in our understanding of EoE, improving diagnostics, and enhancing treatment effectiveness. This review synthesizes current omics research and explores future directions for comprehensively understanding the disease mechanisms in EoE.

RNA sequencing analysis reveals distinct gene expression patterns in infrapatellar fat pads of patients with end-stage osteoarthritis or rheumatoid arthritis

Osteoarthritis (OA) and autoimmune-driven rheumatoid arthritis (RA) are inflammatory joint diseases that share partly similar symptoms but have different, inadequately understood pathogeneses. Adipose tissues, including intra-articular infrapatellar fat pad (IFP), may contribute to their development. Analysis of differentially expressed genes (DEGs) in IFPs could improve the diagnostics of these conditions and help to develop novel treatment strategies. The aim was to identify potentially crucial genes and pathways discriminating OA and RA IFPs using RNA sequencing analysis. We aimed to distinguish genetically distinct patient groups as a starting point for further translational studies with the eventual goal of personalized medicine. Samples were collected from arthritic knees during total knee arthroplasty of sex- and age-matched OA and seropositive RA patients (n = 5-6/group). Metabolic pathways of interest were investigated by whole transcriptome sequencing, and DEGs were analyzed with univariate tests, hierarchical clustering (HC), and pathway analyses. There was significant interindividual variation in mRNA expression patterns, but distinct subgroups of OA and RA patients emerged that reacted similarly to their disease states based on HC. Compared to OA, RA samples showed 703 genes to be upregulated and 691 genes to be downregulated. Signaling pathway analyses indicated that these DEGs had common pathways in lipid metabolism, fatty acid biosynthesis and degradation, adipocytokine and insulin signaling, inflammatory response, and extracellular matrix organization. The divergent mRNA expression profiles in RA and OA suggest contribution of IFP to the regulation of synovial inflammatory processes and articular cartilage degradation and could provide novel diagnostic and therapeutic targets.

Alzheimer's disease-associated protective variant Plcg2-P522R modulates peripheral macrophage function in a sex-dimorphic manner

Genome-wide association studies have identified a protective mutation in the phospholipase C gamma 2 (PLCG2) gene which confers protection against Alzheimer's disease (AD)-associated cognitive decline. Therefore, PLCG2, which is primarily expressed in immune cells, has become a target of interest for potential therapeutic intervention. The protective allele, known as P522R, has been shown to be hyper-morphic in microglia, increasing phagocytosis of amyloid-beta (Aβ), and increasing the release of inflammatory cytokines. However, the effect of this protective mutation on peripheral tissue-resident macrophages, and the extent to which sex modifies this effect, has yet to be assessed. Herein, we show that peripheral macrophages carrying the P522R mutation do indeed show functional differences compared to their wild-type (WT) counterparts, however, these alterations occur in a sex-dependent manner. In macrophages from females, the P522R mutation increases lysosomal protease activity, cytokine secretion, and gene expression associated with cytokine secretion and apoptosis. In contrast, in macrophages from males, the mutation causes decreased phagocytosis and lysosomal protease activity, modest increases in cytokine secretion, and induction of gene expression associated with negative regulation of the immune response. Taken together, these results suggest that the mutation may be conferring different effects dependent on sex and cell type, and highlight the importance of considering sex as a biological variable when assessing the effects of genetic variants and implications for potential immune system-targeted therapies.

Investigating the role of prognostic mitophagy-related genes in non-small cell cancer pathogenesis via multiomics and network-based approach

As one of the most prevalent malignancies, lung cancer displays considerable biological variability in both molecular and clinical characteristics. Lung cancer is broadly categorized into small cell lung cancer (SCLC) and non-small cell lung cancer (NSCLC) with the latter being most prevalent. The primary histological subtypes of NSCLC are lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC). In the present work, we primarily extracted mRNA count data from a publicly accessible database followed by differentially expressed genes (DEGs) and differentially expressed mitophagy-related genes (DEMRGs) identification in case of both LUAD and LUSC cohorts. Next, we identified important DEMRGs via clustering approach followed by enrichment, survival, and mutational analyses. Lastly, the finalized prognostic biomarker was validated using wet-lab experimentations. Primarily, we obtained 986 and 1714 DEGs across LUAD and LUSC cohorts. Only 7 DEMRGs from both cohorts had significant membership values as indicated by the clustering analysis. Most significant pathway, Gene Ontology (GO)-biological process (BP), GO-molecular function (MF), GO-cellular compartment (CC) terms were macroautophagy, GTP metabolic process, magnesium ion binding, mitochondrial outer membrane. Among all, only TDRKH reported significant overall survival (OS) and 14% amplification across LUAD patients. Lastly, we validated TDRKH via immunohistochemistry (IHC) and semi-quantitative polymerase chain reaction (PCR). In conclusion, our findings advocate for the exploration of TDRKH and their genetic alterations in precision oncology therapeutic approaches for LUAD, emphasizing the potential for target-driven therapy and early diagnostics.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-04127-y.

PPARγ antagonism as a new tool for preventing or overcoming endocrine resistance in luminal A breast cancers

PURPOSE

This research investigates the role of PPARγ in the complex molecular events underlying the acquisition of resistance to tamoxifen (Tam) in luminal A breast cancer (BC) cells. Furthermore, it focuses on evaluating the possibility of repurposing Imatinib mesylate, an FDA-approved anticancer agent recently recognized also as a PPARγ antagonist, for the personalized therapy of endocrine-resistant BC with increased PPARγ expression.

METHODS

Differential gene expression between parental and Tam-resistant MCF7 cells was assessed by RNA-seq followed by bioinformatics analysis and validation by RT-qPCR. PPARγ was downregulated by esiRNAs or inhibited by the antagonist GW9662. Cell viability and proliferation were measured by MTT and colony formation assays. Spheroids were prepared from parental and Tam-resistant MCF7 cells. Other luminal A BC cell lines resistant to Tam were generated.

RESULTS

In MCF7-TamR cells, PPARγ and several of its target genes were significantly upregulated. Increased PPARγ expression was due to the modulation of its positive/negative transcriptional regulators. Downregulating PPARγ with esiRNAs or GW9662 effectively killed parental and Tam-resistant cells and spheroids. Imatinib revealed to be as effective as GW9662 in restoring Tam susceptibility of these cells. PPARγ overexpression was also observed in the newly-selected Tam-resistant luminal A BC cells, in which GW9662 and Imatinib restored their susceptibility to Tam.

CONCLUSION

Our findings demonstrate that the overexpression of PPARγ is a frequent occurrence during acquisition of Tam resistance in luminal A BC cells, and that PPARγ antagonism represents an alternative therapeutic approach for the personalized treatment of BC showing dysregulation of this nuclear receptor.

Chronic Alcohol Intake Compromises Lung Immunity by Altering Immunometabolism in Humans and Mouse Models

Chronic alcohol consumption disrupts lung immunity and host defense mechanisms, rendering individuals with alcohol use disorder more susceptible to developing inflammatory lung conditions with poor prognoses. Here, we focused on investigating the molecular and cellular effects of alcohol ingestion on lung immunity in male and female subjects using population-based human lung transcriptomics analysis and an experimental mouse model of chronic alcohol drinking using the National Institute on Alcohol Abuse and Alcoholism alcohol feeding model. Flow cytometry and transcriptomics analyses in lungs revealed a sexually dimorphic effect of chronic alcohol drinking on lung immunity in both human and mouse. Male lungs were more sensitive to chronic alcohol drinking-induced dysregulation of lung immunity compared with female lungs. Furthermore, comparative transcriptomics analysis using lungs and liver samples from matched human and mouse subjects demonstrated that lungs were more sensitive than liver to the effects of alcohol in downregulating immune-related genes and pathways. Furthermore, the transcriptomics analysis provided evidence that immunometabolic change is a central driver in lung alteration by downregulating the immune pathways and upregulating metabolic pathways. Chronic alcohol consumption resulted in reduced mTOR signaling and decreased immune cell populations. The mTOR signaling axis may serve as an upstream regulator of alcohol-induced dysregulation in lung immunity.

Regular exercise suppresses steatosis-associated liver cancer development by degrading E2F1 and c-Myc via circadian gene upregulation

Regular exercise is believed to suppress cancer progression. However, the precise molecular mechanisms by which exercise prevents cancer development remain unclear. In this study, using a steatosis-associated liver cancer mouse model, we found that regular exercise at a speed of 18 m/min for 20 min daily suppressed liver cancer development. To explore the underlying mechanisms, we examined the gene expression profiles in the livers of the exercise and non-exercise groups. The expressions of circadian genes, such as Per1 and Cry2, were upregulated in the exercise group. As circadian rhythm disruption is known to cause various diseases, including cancer, improving circadian rhythm through exercise could contribute to cancer prevention. We further found that the expression of a series of E2F1 and c-Myc target genes that directly affect the proliferation of cancer cells was downregulated in the exercise group. However, the expression of E2F1 and c-Myc was transcriptionally unchanged but degraded at the post-translational level by exercise. Cry2, which is regulated by the Skp1-Cul1-FBXL3 (SCFFBXL3) ubiquitin ligase complex by binding to FBXL3, can form a complex with E2F1 and c-Myc, which we think is the mechanism to degrade them. Our study revealed a previously unknown mechanism by which exercise prevents cancer development.

Sex dimorphism and tissue specificity of gene expression changes in aging mice

BACKGROUND

Aging is a complex process that involves all tissues in an organism and shows sex dimorphism. While transcriptional changes in aging have been well characterized, the majority of studies have focused on a single sex and sex differences in gene expression in aging are poorly understood. In this study, we explore sex dimorphism in gene expression in aging mice across three tissues.

METHODS

We collected gastrocnemius muscle, liver and white adipose tissue from young (6 months, n = 14) and old (24 months, n = 14) female and male C57BL/6NIA mice and performed RNA-seq. To investigate sex dimorphism in aging, we considered two levels of comparisons: (a) differentially expressed genes between females and males in the old age group and (b) comparisons between females and males across the aging process. We utilized differential expression analysis and gene feature selection to investigate candidate genes. Gene set enrichment analysis was performed to identify candidate molecular pathways. Furthermore, we performed a co-expression network analysis and chose the gene module(s) associated with aging independent of sex or tissue-type.

RESULTS

We identified both tissue-specific and tissue-independent genes associated with sex dimorphism in aged mice. Unique differentially expressed genes between old males and females across tissues were mainly enriched for pathways related to specific tissue function. We found similar results when exploring sex differences in the aging process, with the exception that in the liver genes enriched for lipid metabolism and digestive system were identified in both females and males. Combining enriched pathways across analyses, we identified amino acid metabolism, digestive system, and lipid metabolism as the core mechanisms of sex dimorphism in aging. Although the vast majority of age-related genes were sex and tissue specific, we identified 127 hub genes contributing to aging independent of sex and tissue that were enriched for the immune system and signal transduction.

CONCLUSIONS

There are clear sex differences in gene expression in aging across liver, muscle and white adipose. Core pathways, including amino acid metabolism, digestive system and lipid metabolism, contribute to sex differences in aging.

Heterogeneous enhancer states orchestrate β cell responses to metabolic stress

Obesity-induced β cell dysfunction contributes to the onset of type 2 diabetes. Nevertheless, elucidating epigenetic mechanisms underlying islet dysfunction at single cell level remains challenging. Here we profile single-nuclei RNA along with enhancer marks H3K4me1 or H3K27ac in islets from lean or obese mice. Our study identifies distinct gene signatures and enhancer states correlating with β cell dysfunction trajectory. Intriguingly, while many metabolic stress-induced genes exhibit concordant changes in both H3K4me1 and H3K27ac at their enhancers, expression changes of specific subsets are solely attributable to either H3K4me1 or H3K27ac dynamics. Remarkably, a subset of H3K4me1+H3K27ac- primed enhancers prevalent in lean β cells and occupied by FoxA2 are largely absent after metabolic stress. Lastly, cell-cell communication analysis identified the nerve growth factor (NGF) as protective paracrine signaling for β cells through repressing ER stress. In summary, our findings define the heterogeneous enhancer responses to metabolic challenges in individual β cells.

LCoRL Regulates Growth and Metabolism

Genome-wide association studies (GWAS) in humans and livestock have identified genes associated with metabolic traits. However, the causality of many of these genes on metabolic homeostasis is largely unclear due to a lack of detailed functional analyses. Here we report ligand-dependent corepressor-like (LCoRL) as a metabolic regulator for body weight and glucose homeostasis. Although GWAS data show that LCoRL is strongly associated with body size, glucose homeostasis, and other metabolic traits in humans and livestock, functional investigations had not been performed. We generated Lcorl knockout mice (Lcorl-/-) and characterized the metabolic traits. We found that Lcorl-/- pups are born smaller than the wild-type (WT) littermates before reaching normal weight by 7 to 9 weeks of age. While aging, Lcorl-/- mice remain lean compared to WT mice, which is associated with a decrease in daily food intake. Glucose tolerance and insulin sensitivity are improved in Lcorl-/- mice. Mechanistically, this stunted growth is linked to a reduction of circulating levels of IGF-1. The expression of the genes downstream of GH signaling and the genes involved in glucose and lipid metabolism are altered in the liver of Lcorl-/- mice. Furthermore, Lcorl-/- mice are protected against a high-fat diet challenge and show reduced exercise capacity in an exercise stress test. Collectively, our results are congruent with many of the metabolic parameters linked to the Lcorl locus as reported in GWAS in humans and livestock.

SARS-CoV-2 variants mediated tissue-specific metabolic reprogramming determines the disease pathophysiology in a hamster model

Despite significant effort, a clear understanding of host tissue-specific responses and their implications for immunopathogenicity against the severe acute respiratory syndrome coronavirus2 (SARS-CoV-2) variant infection has remained poorly defined. To shed light on the interaction between tissues and SARS-CoV-2 variants, we sought to characterize the complex relationship among acute multisystem manifestations, dysbiosis of the gut microbiota, and the resulting implications for SARS-CoV-2 variant-specific immunopathogenesis in the Golden Syrian Hamster (GSH) model using multi-omics approaches. Our investigation revealed the presence of increased SARS-CoV-2 genomic RNA in diverse tissues of delta-infected GSH compared to the omicron variant. Multi-omics analyses uncovered distinctive metabolic responses between the delta and omicron variants, with the former demonstrating dysregulation in synaptic transmission proteins associated with neurocognitive disorders. Additionally, delta-infected GSH exhibited an altered fecal microbiota composition, marked by increased inflammation-associated taxa and reduced commensal bacteria compared to the omicron variant. These findings underscore the SARS-CoV-2-mediated tissue insult, characterized by modified host metabolites, neurological protein dysregulation, and gut dysbiosis, highlighting the compromised gut-lung-brain axis during acute infection.

ZnO Nanoparticles as Potent Inducers of Dermal Immunosuppression in Contact Hypersensitivity in Mice

Nanosized zinc oxide (nZnO) metal particles are used in skin creams and sunscreens to enhance their texture and optical properties as UV filters. Despite their common use, little is known about the molecular mechanisms of nZnO exposure on damaged skin. We studied the effects of topically applied nZnO particles on allergic skin inflammation in an oxazolone (OXA)-induced contact hypersensitivity (CHS) mouse model. We investigated whether exposure to nZnO during the sensitization or challenge phase would induce immunological changes and modulate transcriptional responses. We followed skin thickness, cellular infiltration, and changes in the local transcriptome up to 28 days after the challenge. The responses peaked at 24 h and were fully resolved by 28 days. Co-exposure to nZnO and hapten did not interfere with the formation of the sensitization process. Conversely, during the hapten challenge, the application of nZnO fully suppressed the development of the CHS response by the inhibition of pro-inflammatory pathways, secretion of pro-inflammatory cytokines, and proliferation of immune cells. In differentiated and stimulated THP-1 cells and the CHS mouse model, we found that nZnO particles and Zn ions contributed to anti-inflammatory responses. The immunosuppressive properties of nZnO in inflamed skin are mediated by impaired IL-1R-, CXCR2-, and LTB4-mediated pathways. nZnO-induced dermal immunosuppression may be beneficial for individuals with contact allergies who use nZnO-containing cosmetic products. Our findings also provide a deeper understanding of the mechanisms of nZnO, which could be considered when developing nanoparticle-containing skin products.