The NCATS BioPlanet - An Integrated Platform for Exploring the Universe of Cellular Signaling Pathways for Toxicology, Systems Biology, and Chemical Genomics

NKX2‑1 copy number alterations are associated with oncogenic, immunological and prognostic remodeling in non‑small cell lung cancer

NK2 homeobox 1 (NKX2-1) copy number alterations (CNAs) are frequently observed in lung cancer. However, little is known about the complete landscape of focal alterations in NKX2-1 copy number (CN), their clinical significance and their therapeutic implications in non-small cell lung cancer (NSCLC). The correlations between NKX2-1 expression and EGFR driver mutations and programmed death ligand 1 (PD-L1) co-expression were studied using immunohistochemistry and PCR from the tumors of recruited Filipino patients (n=45). Clinical features of NSCLC with NKX2-1 CNAs were resolved at the tumor and clonal levels using the molecular profiles of patients with lung adenocarcinoma and lung squamous cell carcinoma from The Cancer Genome Atlas (n=1,130), and deconvoluted single-cell RNA-seq data from the Bivona project (n=1,654), respectively. Despite a significant and positive correlation between expression and CN (r=0.264; P<0.001), NKX2-1 CNAs exerted a stronger influence on the combined EGFR and PD-L1 status of NSCLC tumors than expression. NKX2-1 CN gain was prognostic of favorable survival (P=0.018) and a better response to targeted therapy. NKX2-1 CN loss predicted a worse survival (P=0.041). Mutational architecture in the Y-chromosome differentiated the two prognostic groups. There were 19,941 synonymous mutations and 1,408 genome-wide CN perturbations associated with NKX2-1 CNAs. Tumors with NKX2-1 CN gain expressed lymphocyte markers more heterogeneously than those with CN loss. Higher expression of tumor-infiltrating lymphocyte gene signatures in CN gain was prognostic of longer disease-free survival (P=0.005). Tumors with NKX2-1 CN gain had higher B-cell (P<0.001) and total T-cell estimates (P=0.003). NKX2-1 CN loss was associated with immunologically colder tumors due to higher M2 macrophage infiltrates (P=0.011) and higher expression of immune checkpoint proteins, CD274 (P=0.025), VTCN1 (P<0.001) and LGALS9 (P=0.002). In conclusion, NKX2-1 CNAs are associated with tumors that exhibit clinically diverse characteristics, and with unique oncogenic, immunological and prognostic signatures.

Knocking Out Sodium Glucose-Linked Transporter 5 Prevents Fructose-Induced Renal Oxidative Stress and Salt-Sensitive Hypertension

BACKGROUND

A fructose high-salt (FHS) diet increases systolic blood pressure and Ang II (angiotensin II)-stimulated proximal tubule (PT) superoxide (O2-) production. These increases are prevented by scavenging O2- or an Ang II type 1 receptor antagonist. SGLT4 (sodium glucose-linked cotransporters 4) and SGLT5 are implicated in PT fructose reabsorption, but their roles in fructose-induced hypertension are unclear. We hypothesized that PT fructose reabsorption by SGLT5 initiates a genetic program enhancing Ang II-stimulated oxidative stress in males and females, thereby causing fructose-induced salt-sensitive hypertension.

METHODS

We measured systolic blood pressure in male and female Sprague-Dawley (wild type [WT]), SGLT4 knockout (-/-), and SGLT5-/- rats. Then, we measured basal and Ang II-stimulated (37 nmol/L) O2- production by PTs and conducted gene coexpression network analysis.

RESULTS

In male WT and female WT rats, FHS increased systolic blood pressure by 15±3 (n=7; P<0.0027) and 17±4 mm Hg (n=9; P<0.0037), respectively. Male and female SGLT4-/- had similar increases. Systolic blood pressure was unchanged by FHS in male and female SGLT5-/-. In male WT and female WT fed FHS, Ang II stimulated O2- production by 14±5 (n=6; P<0.0493) and 8±3 relative light units/µg protein/s (n=7; P<0.0218), respectively. The responses of SGTL4-/- were similar. Ang II did not stimulate O2- production in tubules from SGLT5-/-. Five gene coexpression modules were correlated with FHS. These correlations were completely blunted in SGLT5-/- and partially blunted by chronically scavenging O2- with tempol.

CONCLUSIONS

SGLT5-mediated PT fructose reabsorption is required for FHS to augment Ang II-stimulated proximal nephron O2- production, and increases in PT oxidative stress likely contribute to FHS-induced hypertension.

Bone morphogenetic protein (BMP) signaling determines neuroblastoma cell fate and sensitivity to retinoic acid

Retinoic acid (RA) is a standard-of-care neuroblastoma drug thought to be effective by inducing differentiation. Curiously, RA has little effect on primary human tumors during upfront treatment but can eliminate neuroblastoma cells from the bone marrow during post-chemo consolidation therapy-a discrepancy that has never been explained. To investigate this, we treated a large cohort of neuroblastoma cell lines with RA and observed that the most RA-sensitive cells predominantly undergo apoptosis or senescence, rather than differentiation. We conducted genome-wide CRISPR knockout screens under RA treatment, which identified BMP signaling as controlling the apoptosis/senescence vs differentiation cell fate decision and determining RA's overall potency. We then discovered that BMP signaling activity is markedly higher in neuroblastoma patient samples at bone marrow metastatic sites, providing a plausible explanation for RA's ability to clear neuroblastoma cells specifically from the bone marrow, seemingly mimicking interactions between BMP and RA during normal development.

Transcriptomic point of departure determination: a comparison of distribution-based and gene set-based approaches

A key step in assessing the potential human and environmental health risks of industrial and agricultural chemicals is to determine the toxicity point of departure (POD), which is the highest dose level that causes no adverse effect. Transcriptomic POD (tPOD) values have been suggested to accurately estimate toxicity POD values. One step in the most common approach for tPOD determination involves mapping genes to annotated gene sets, a process that might lead to substantial information loss particularly in species with poor gene annotation. Alternatively, methods that calculate tPOD values directly from the distribution of individual gene POD values omit this mapping step. Using rat transcriptome data for 79 molecules obtained from Open TG-GATEs (Toxicogenomics Project Genomics Assisted Toxicity Evaluation System), the hypothesis was tested that methods based on the distribution of all individual gene POD values will give a similar tPOD value to that obtained via the gene set-based method. Gene set-based tPOD values using four different gene set structures were compared to tPOD values from five different individual gene distribution methods. Results revealed a high tPOD concordance for all methods tested, especially for molecules with at least 300 dose-responsive probesets: for 90% of those molecules, the tPOD values from all methods were within 4-fold of each other. In addition, random gene sets based upon the structure of biological knowledge-derived gene sets produced tPOD values with a median absolute fold change of 1.3-1.4 when compared to the original biological knowledge-derived gene set counterparts, suggesting that little biological information is used in the gene set-based tPOD generation approach. These findings indicate using individual gene distributions to calculate a tPOD is a viable and parsimonious alternative to using gene sets. Importantly, individual gene distribution-based tPOD methods do not require knowledge of biological organization and can be applied to any species including those with poorly annotated gene sets.

The genetic regulation of the gastric transcriptome is associated with metabolic and obesity-related traits and diseases

Tissue-specific gene expression and gene regulation lead to a better understanding of tissue-specific physiology and pathophysiology. We analyzed the transcriptome and genetic regulatory profiles of two distinct gastric sites, corpus and antrum, to identify tissue-specific gene expression and its regulation. Gastric corpus and antrum mucosa biopsies were collected during routine gastroscopies from up to 431 healthy individuals. We obtained genotype and transcriptome data and performed transcriptome profiling and expression quantitative trait locus (eQTL) studies. We further used data from genome-wide association studies (GWAS) of various diseases and traits to partition their heritability and to perform transcriptome-wide association studies (TWAS). The transcriptome data from corpus and antral mucosa highlights the heterogeneity of gene expression in the stomach. We identified enriched pathways revealing distinct and common physiological processes in gastric corpus and antrum. Furthermore, we found an enrichment of the single nucleotide polymorphism (SNP)-based heritability of metabolic, obesity-related, and cardiovascular traits and diseases by considering corpus- and antrum-specifically expressed genes. Particularly, we could prioritize gastric-specific candidate genes for multiple metabolic traits, like NQO1 which is involved in glucose metabolism, MUC1 which contributes to purine and protein metabolism or RAB27B being a regulator of weight and body composition. Our findings show that gastric corpus and antrum vary in their transcriptome and genetic regulatory profiles indicating physiological differences which are mostly related to digestion and epithelial protection. Moreover, our findings demonstrate that the genetic regulation of the gastric transcriptome is linked to biological mechanisms associated with metabolic, obesity-related, and cardiovascular traits and diseases. NEW & NOTEWORTHY We analyzed the transcriptomes and genetic regulatory profiles of gastric corpus and for the first time also of antrum mucosa in 431 healthy individuals. Through tissue-specific gene expression and eQTL analyses, we uncovered unique and common physiological processes across both primary gastric sites. Notably, our findings reveal that stomach-specific eQTLs are enriched in loci associated with metabolic traits and diseases, highlighting the pivotal role of gene expression regulation in gastric physiology and potential pathophysiology.

The landscape of immune dysregulation in pediatric sepsis at a single-cell resolution

Recognizing immune dysregulation as a hallmark of sepsis pathophysiology, leukocytes have attracted major attention of investigation. While adult and pediatric sepsis are clinically distinct, their immunological delineation remains limited. Single cell technologies facilitated the characterization of immune signatures. We tackled to delineate immunological profiles of pediatric sepsis at a single-cell level by analyzing blood samples from six septic children, at both acute and recovery phases, and four healthy children. 16 single-cell transcriptomic datasets were analyzed and compared to adult sepsis dataset. We showed a unique shift in neutrophil subpopulations and functions between acute and recovery phases, along with the regulatory role of resistin. Neutrophil signatures were comparable between adult and pediatric sepsis. Innate-like CD4 T cells were predominantly and uniquely observed in acute phase of pediatric sepsis. Our study serves as a rich source of information about the phenotypic diversity and trajectory of circulating immune cells during pediatric sepsis.

Hypoxia-inducible factor 1 recruits FACT and RNF20/40 to mediate histone ubiquitination and transcriptional activation of target genes

Hypoxia-inducible factor 1 (HIF-1) is a transcriptional activator that mediates cellular adaptation to decreased oxygen availability. HIF-1 recruits chromatin-modifying enzymes leading to changes in histone acetylation, citrullination, and methylation at target genes. Here, we demonstrate that hypoxia-inducible gene expression in estrogen receptor (ER)-positive MCF7 and ER-negative SUM159 human breast cancer cells requires the histone H2A/H2B chaperone facilitates chromatin transcription (FACT) and the H2B ubiquitin ligase RING finger protein 20/40 (RNF20/40). Knockdown of FACT or RNF20/40 expression leads to decreased transcription initiation and elongation at HIF-1 target genes. Mechanistically, FACT and RNF20/40 are recruited to hypoxia response elements (HREs) by HIF-1 and stabilize binding of HIF-1 (and each other) at HREs. Hypoxia induces the monoubiquitination of histone H2B at lysine 120 at HIF-1 target genes in an HIF-1-dependent manner. Together, these findings delineate a cooperative molecular mechanism by which FACT and RNF20/40 stabilize multiprotein complex formation at HREs and mediate histone ubiquitination to facilitate HIF-1 transcriptional activity.

phuEGO: A Network-Based Method to Reconstruct Active Signaling Pathways From Phosphoproteomics Datasets

Signaling networks are critical for virtually all cell functions. Our current knowledge of cell signaling has been summarized in signaling pathway databases, which, while useful, are highly biased toward well-studied processes, and do not capture context specific network wiring or pathway cross-talk. Mass spectrometry-based phosphoproteomics data can provide a more unbiased view of active cell signaling processes in a given context, however, it suffers from low signal-to-noise ratio and poor reproducibility across experiments. While progress in methods to extract active signaling signatures from such data has been made, there are still limitations with respect to balancing bias and interpretability. Here we present phuEGO, which combines up-to-three-layer network propagation with ego network decomposition to provide small networks comprising active functional signaling modules. PhuEGO boosts the signal-to-noise ratio from global phosphoproteomics datasets, enriches the resulting networks for functional phosphosites and allows the improved comparison and integration across datasets. We applied phuEGO to five phosphoproteomics data sets from cell lines collected upon infection with SARS CoV2. PhuEGO was better able to identify common active functions across datasets and to point to a subnetwork enriched for known COVID-19 targets. Overall, phuEGO provides a flexible tool to the community for the improved functional interpretation of global phosphoproteomics datasets.

Craniofacial chondrogenesis in organoids from human stem cell-derived neural crest cells

Knowledge of cell signaling pathways that drive human neural crest differentiation into craniofacial chondrocytes is incomplete, yet essential for using stem cells to regenerate craniomaxillofacial structures. To accelerate translational progress, we developed a differentiation protocol that generated self-organizing craniofacial cartilage organoids from human embryonic stem cell-derived neural crest stem cells. Histological staining of cartilage organoids revealed tissue architecture and staining typical of elastic cartilage. Protein and post-translational modification (PTM) mass spectrometry and snRNA-seq data showed that chondrocyte organoids expressed robust levels of cartilage extracellular matrix (ECM) components: many collagens, aggrecan, perlecan, proteoglycans, and elastic fibers. We identified two populations of chondroprogenitor cells, mesenchyme cells and nascent chondrocytes, and the growth factors involved in paracrine signaling between them. We show that ECM components secreted by chondrocytes not only create a structurally resilient matrix that defines cartilage, but also play a pivotal autocrine cell signaling role in determining chondrocyte fate.

Neuronal-specific methylome and hydroxymethylome analysis reveal significant loci associated with alcohol use disorder

Background: Alcohol use disorder (AUD) is a complex condition associated with adverse health consequences that affect millions of individuals worldwide. Epigenetic modifications, including DNA methylation (5 mC), have been associated with AUD and other alcohol-related traits. Epigenome-wide association studies (EWAS) have identified differentially methylated genes associated with AUD in human peripheral and brain tissue. More recently, epigenetic studies of AUD have also evaluated DNA hydroxymethylation (5 hmC) in the human brain. However, most of the epigenetic work in postmortem brain tissue has examined bulk tissue. In this study, we investigated neuronal-specific 5 mC and 5 hmC alterations at CpG sites associated with AUD in the human orbitofrontal cortex (OFC). Methods: Neuronal nuclei from the OFC were evaluated in 34 human postmortem brain samples (10 AUD, 24 non-AUD). Reduced representation oxidative bisulfite sequencing was used to assess 5 mC and 5 hmC at the genome-wide level. Differential 5 mC and 5 hmC were evaluated using the methylKit R package and significance was set at false discovery rate < 0.05 and differential methylation > 2. Functional enrichment analyses were performed, and gene-level convergence was evaluated in an independent dataset that assessed 5 mC and 5 hmC of AUD in bulk cortical tissue. Results: We identified 417 5 mC and 363 5hmC significant differential CpG sites associated with AUD, with 59% in gene promoters. Some of the identified genes have been previously implicated in alcohol consumption, including SYK, DNMT3A for 5 mC, GAD1, DLX1, DLX2, for 5 hmC and GATA4 in both. Convergence with a previous AUD 5 mC and 5 hmC study was observed for 28 genes. We also identified 5 and 35 differential regions for 5 mC and 5 hmC, respectively. Lastly, GWAS enrichment analysis showed an association with AUD for differential 5 mC genes. Discussion: This study reveals neuronal-specific methylome and hydroxymethylome dysregulation associated with AUD, identifying both previously reported and potentially novel gene associations with AUD. Our findings provide new insights into the epigenomic dysregulation of AUD in the human brain.

Loss-of-function mutation in PRMT9 causes abnormal synapse development by dysregulation of RNA alternative splicing

Protein arginine methyltransferase 9 (PRMT9) is a recently identified member of the PRMT family, yet its biological function remains largely unknown. Here, by characterizing an intellectual disability associated PRMT9 mutation (G189R) and establishing a Prmt9 conditional knockout (cKO) mouse model, we uncover an important function of PRMT9 in neuronal development. The G189R mutation abolishes PRMT9 methyltransferase activity and reduces its protein stability. Knockout of Prmt9 in hippocampal neurons causes alternative splicing of ~1900 genes, which likely accounts for the aberrant synapse development and impaired learning and memory in the Prmt9 cKO mice. Mechanistically, we discover a methylation-sensitive protein-RNA interaction between the arginine 508 (R508) of the splicing factor 3B subunit 2 (SF3B2), the site that is exclusively methylated by PRMT9, and the pre-mRNA anchoring site, a cis-regulatory element that is critical for RNA splicing. Additionally, using human and mouse cell lines, as well as an SF3B2 arginine methylation-deficient mouse model, we provide strong evidence that SF3B2 is the primary methylation substrate of PRMT9, thus highlighting the conserved function of the PRMT9/SF3B2 axis in regulating pre-mRNA splicing.

The Immunometabolic Gene N-Acetylglucosamine Kinase Is Uniquely Involved in the Heritability of Multiple Sclerosis Severity

The clinical severity of multiple sclerosis (MS), an autoimmune disorder of the central nervous system, is thought to be determined by environmental and genetic factors that have not yet been identified. In a recent genome-wide association study (GWAS), a single nucleotide polymorphism (SNP), rs10191329, has been associated with MS severity in two large independent cohorts of patients. Different approaches were followed by the authors to prioritize the genes that are transcriptionally regulated by such an SNP. It was concluded that the identified SNP regulates a group of proximal genes involved in brain resilience and cognitive abilities rather than immunity. Here, by conducting an alternative strategy for gene prioritization, we reached the opposite conclusion. According to our re-analysis, the main target of rs10191329 is N-Acetylglucosamine Kinase (NAGK), a metabolic gene recently shown to exert major immune functions via the regulation of the nucleotide-binding oligomerization domain-containing protein 2 (NOD2) pathway. To gain more insights into the immunometabolic functions of NAGK, we analyzed the currently known list of NAGK protein partners. We observed that NAGK integrates a dense network of human proteins that are involved in glucose metabolism and are highly expressed by classical monocytes. Our findings hold potentially major implications for the understanding of MS pathophysiology.

Functional and Multi-Omics Effects of an Optimized CRISPR-Mediated FURIN Depletion in U937 Monocytes

The pro-protein convertase FURIN (PCSK3) is implicated in a wide range of normal and pathological biological processes such as infectious diseases, cancer and cardiovascular diseases. Previously, we performed a systemic inhibition of FURIN in a mouse model of atherosclerosis and demonstrated significant plaque reduction and alterations in macrophage function. To understand the cellular mechanisms affected by FURIN inhibition in myeloid cells, we optimized a CRISPR-mediated gene deletion protocol for successfully deriving hemizygous (HZ) and nullizygous (NZ) FURIN knockout clones in U937 monocytic cells using lipotransfection-based procedures and a dual guide RNA delivery strategy. We observed differences in monocyte and macrophage functions involving phagocytosis, lipid accumulation, cell migration, inflammatory gene expression, cytokine release patterns, secreted proteomics (cytokines) and whole-genome transcriptomics between wild-type, HZ and NZ FURIN clones. These studies provide a mechanistic basis on the possible roles of myeloid cell FURIN in cardiovascular disorders.

Serine/arginine-rich splicing factor 7 promotes the type I interferon response by activating Irf7 transcription

Tight regulation of macrophage immune gene expression is required to fight infection without risking harmful inflammation. The contribution of RNA-binding proteins (RBPs) to shaping the macrophage response to pathogens remains poorly understood. Transcriptomic analysis reveals that a member of the serine/arginine-rich (SR) family of mRNA processing factors, SRSF7, is required for optimal expression of a cohort of interferon-stimulated genes in macrophages. Using genetic and biochemical assays, we discover that in addition to its canonical role in regulating alternative splicing, SRSF7 drives transcription of interferon regulatory transcription factor 7 (IRF7) to promote antiviral immunity. At the Irf7 promoter, SRSF7 maximizes STAT1 transcription factor binding and RNA polymerase II elongation via cooperation with the H4K20me1 histone methyltransferase KMT5a (SET8). These studies define a role for an SR protein in activating transcription and reveal an RBP-chromatin network that orchestrates macrophage antiviral gene expression.

Novel lncRNAs LINC01221, RP11-472G21.2 and CRNDE are markers of differential expression in pediatric patients with T cell acute lymphoblastic leukemia

INTRODUCTION

Pediatric T-cell acute lymphoblastic leukemia (T-ALL) poses significant challenges due to its aggressive nature and resistance to standard treatments. Long non-coding RNAs (lncRNAs) have emerged as potential biomarkers and therapeutic targets in leukemia. This study aims to characterize the lncRNA landscape in pediatric T-ALL, identify specific lncRNAs signatures, and assess their clinical relevance.

METHODS

RNA sequencing was performed on T-ALL patient and control samples. Differential expression analysis identified dysregulated lncRNAs and mRNAs. Functional enrichment analysis revealed potential roles of these lncRNAs in cancer pathogenesis. Validation of candidate lncRNAs was conducted using real-time PCR. Clinical correlations were assessed, including associations with patients' clinical characteristics and survival outcomes.

RESULTS

Analysis identified 674 dysregulated lncRNAs in pediatric T-ALL, with LINC01221 and CRNDE showing the most interactions in cancer progression pathways. Functional enrichment indicated involvement in apoptosis, survival, proliferation, and metastasis. Top 10 lncRNAs based on adjusted p value < 0.05 and Fold Change > 2 were selected for validation. Seven lncRNAs LINC01221, PCAT18, LINC00977, RP11-620J15.3, RP11-472G21.2, CTD-2291D10.4, and CRNDE showed correlation with RNA sequencing data. RP11-472G21.2 and CTD-2291D10.4 were highly expressed in T-ALL patients, with RP11-620J15.3 correlating significantly with better overall survival (p = 0.0007) at a median follow up of 32 months. The identified lncRNAs were further analysed in B-ALL patients. Distinct lncRNAs signatures were noted, distinguishing T-ALL from B-ALL and healthy controls, with lineage-specific overexpression of LINC01221 (p < 0.0001), RP11-472G21.2 (p < 0.001) and CRNDE (p = 0.04) in T-ALL.

CONCLUSION

This study provides insights into the lncRNA landscape of pediatric T-ALL, offering potential diagnostic and prognostic markers. RP11-620J15.3 emerges as a promising prognostic marker, and distinct lncRNAs signatures may aid in the differentiation of T-ALL subtypes. Further research with larger cohorts is warranted to validate these findings and advance personalized treatment strategies for pediatric T-ALL patients.

Correlation of Myeloid-Derived Suppressor Cell Expansion with Upregulated Transposable Elements in Severe COVID-19 Unveiled in Single-Cell RNA Sequencing Reanalysis

Some studies have investigated the potential role of transposable elements (TEs) in COVID-19 pathogenesis and complications. However, to the best of our knowledge, there is no study to examine the possible association of TE expression in cell functions and its potential role in COVID-19 immune response at the single-cell level. In this study, we reanalyzed single-cell RNA seq data of bronchoalveolar lavage (BAL) samples obtained from six severe COVID-19 patients and three healthy donors to assess the probable correlation of TE expression with the immune responses induced by the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) in COVID-19 patients. Our findings indicate that the expansion of myeloid-derived suppressor cells (MDSCs) may be a characteristic feature of COVID-19. Additionally, a significant increase in TE expression in MDSCs was observed. This upregulation of TEs in COVID-19 may be linked to the adaptability of these cells in response to their microenvironments. Furthermore, it appears that the identification of overexpressed TEs by pattern recognition receptors (PRRs) in MDSCs may enhance the suppressive capacity of these cells. Thus, this study emphasizes the crucial role of TEs in the functionality of MDSCs during COVID-19.

The landscape of immune dysregulation in pediatric sepsis at a single-cell resolution

Recognizing immune dysregulation as a hallmark of sepsis pathophysiology, leukocytes have attracted major attention of investigation. While adult and pediatric sepsis are clinically distinct, their immunological delineation remains limited. Breakthrough of single cell technologies facilitated the characterization of immune signatures. We tackled to delineate immunological profiles of pediatric sepsis at a single-cell level by analyzing blood samples from six septic children, at both acute and recovery phases, and four healthy children. 16 single-cell transcriptomic datasets (96,156 cells) were analyzed and compared to adult sepsis dataset. We showed a unique shift in neutrophil subpopulations and functions between acute and recovery phases, along with examining the regulatory role of resistin. Neutrophil signatures were comparable between adult and pediatric sepsis. Innate-like CD4 T cells were predominantly and uniquely observed in acute phase of pediatric sepsis. Our study provides a thorough and comprehensive understanding of immune dysregulation in pediatric sepsis.

Discovery of circulating miRNAs as biomarkers of chronic Chagas heart disease via a small RNA-Seq approach

Chagas disease affects approximately 7 million people worldwide in Latin America and is a neglected tropical disease. Twenty to thirty percent of chronically infected patients develop chronic Chagas cardiomyopathy decades after acute infection. Identifying biomarkers of Chagas disease progression is necessary to develop better therapeutic and preventive strategies. Circulating microRNAs are increasingly reliable biomarkers of disease and therapeutic targets. To identify new circulating microRNAs for Chagas disease, we performed exploratory small RNA sequencing from the plasma of patients and performed de novo miRNA prediction, identifying potential new microRNAs. The levels of the new microRNAs temporarily named miR-Contig-1519 and miR-Contig-3244 and microRNAs that are biomarkers for nonchagasic cardiomyopathies, such as miR-148a-3p and miR-224-5p, were validated by quantitative reverse transcription. We found a specific circulating microRNA signature defined by low miR-Contig-3244, miR-Contig-1519, and miR-148a-3 levels but high miR-224-5p levels for patients with chronic Chagas disease. Finally, we predicted in silico that these altered circulating microRNAs could affect the expression of target genes involved in different cellular pathways and biological processes, which we will explore in the future.

Integrative single-cell characterization of a frugivorous and an insectivorous bat kidney and pancreas

Frugivory evolved multiple times in mammals, including bats. However, the cellular and molecular components driving it remain largely unknown. Here, we use integrative single-cell sequencing (scRNA-seq and scATAC-seq) on insectivorous (Eptesicus fuscus; big brown bat) and frugivorous (Artibeus jamaicensis; Jamaican fruit bat) bat kidneys and pancreases and identify key cell population, gene expression and regulatory differences associated with the Jamaican fruit bat that also relate to human disease, particularly diabetes. We find a decrease in loop of Henle and an increase in collecting duct cells, and differentially active genes and regulatory elements involved in fluid and electrolyte balance in the Jamaican fruit bat kidney. The Jamaican fruit bat pancreas shows an increase in endocrine and a decrease in exocrine cells, and differences in genes and regulatory elements involved in insulin regulation. We also find that these frugivorous bats share several molecular characteristics with human diabetes. Combined, our work provides insights from a frugivorous mammal that could be leveraged for therapeutic purposes.

PathDIP 5: improving coverage and making enrichment analysis more biologically meaningful

Pathway Data Integration Portal (PathDIP) is an integrated pathway database that was developed to increase functional gene annotation coverage and reduce bias in pathway enrichment analysis. PathDIP 5 provides multiple improvements to enable more interpretable analysis: users can perform enrichment analysis using all sources, separate sources or by combining specific pathway subsets; they can select the types of sources to use or the types of pathways for the analysis, reducing the number of resulting generic pathways or pathways not related to users' research question; users can use API. All pathways have been mapped to seven representative types. The results of pathway enrichment can be summarized through knowledge-based pathway consolidation. All curated pathways were mapped to 53 pathway ontology-based categories. In addition to genes, pathDIP 5 now includes metabolites. We updated existing databases, included two new sources, PathBank and MetabolicAtlas, and removed outdated databases. We enable users to analyse their results using Drugst.One, where a drug-gene network is created using only the user's genes in a specific pathway. Interpreting the results of any analysis is now improved by multiple charts on all the results pages. PathDIP 5 is freely available at https://ophid.utoronto.ca/pathDIP.

Phenelzine protects against acetaminophen induced apoptosis in HepG2 cells

Acetaminophen (APAP) overdosing is the most common cause of drug-induced liver failure. Despite extensive study, N-acetylcysteine is currently the only antidote utilized for treatment. The purpose of this study was to evaluate the effect and mechanisms of phenelzine, an FDA-approved antidepressant, on APAP-induced toxicity in HepG2 cells. The human liver hepatocellular cell line HepG2 was used to investigate APAP-induced cytotoxicity. The protective effects of phenelzine were determined by examining the cell viability, combination index calculation, Caspase 3/7 activation, Cytochrome c release, H2O2 levels, NO levels, GSH activity, PERK protein levels, and pathway enrichment analysis. Elevated H2O2 production and decreased glutathione (GSH) levels were indicators of APAP-induced oxidative stress. The combination index of 2.04 indicated that phenelzine had an antagonistic effect on APAP-induced toxicity. When compared to APAP alone, phenelzine treatment considerably reduced caspase 3/7 activation, cytochrome c release, and H2O2 generation. However, phenelzine had minimal effect on NO and GSH levels and did not alleviate ER stress. Pathway enrichment analysis revealed a potential connection between APAP toxicity and phenelzine metabolism. These findings suggested that phenelzine's protective effect against APAP-induced cytotoxicity could be attributed to the drug's capacity to reduce APAP-mediated apoptotic signaling.

Transcriptome changes in the nucleus of the solitary tract induced by repeated stress, alcohol dependence, or stress-induced drinking in dependent mice

Stress increases alcohol consumption in dependent animals and contributes to the development of alcohol use disorder. The nucleus of the solitary tract (NTS) is a critical brainstem region for integrating and relaying central and peripheral signals to regulate stress responses, but it is not known if it plays a role in alcohol dependence- or in stress-induced escalations in alcohol drinking in dependent mice. Here, we used RNA-sequencing and bioinformatics analyses to study molecular adaptations in the NTS of C57BL/6J male mice that underwent an ethanol drinking procedure that uses exposure to chronic intermittent ethanol (CIE) vapor, forced swim stress (FSS), or both conditions (CIE + FSS). Transcriptome profiling was performed at three different times after the last vapor cycle (0-hr, 72-hr, and 168-hr) to identify changes in gene expression associated with different stages of ethanol intoxication and withdrawal. In the CIE and CIE + FSS groups at 0-hr, there was upregulation of genes enriched for cellular response to type I interferon (IFN) and type I IFN- and cytokine-mediated signaling pathways, while the FSS group showed upregulation of neuronal genes. IFN signaling was the top gene network positively correlated with ethanol consumption levels in the CIE and CIE + FSS groups. Results from different analyses (differential gene expression, weighted gene coexpression network analysis, and rank-rank hypergeometric overlap) indicated that activation of type I IFN signaling would be expected to increase ethanol consumption. The CIE and CIE + FSS groups also shared an immune signature in the NTS as has been demonstrated in other brain regions after chronic ethanol exposure. A temporal-based clustering analysis revealed a unique expression pattern in the CIE + FSS group that suggests the interaction of these two stressors produces adaptations in synaptic and glial functions that may drive stress-induced drinking.

Exploring the effects of experimental parameters and data modeling approaches on in vitro transcriptomic point-of-departure estimates

Multiple new approach methods (NAMs) are being developed to rapidly screen large numbers of chemicals to aid in hazard evaluation and risk assessments. High-throughput transcriptomics (HTTr) in human cell lines has been proposed as a first-tier screening approach for determining the types of bioactivity a chemical can cause (activation of specific targets vs. generalized cell stress) and for calculating transcriptional points of departure (tPODs) based on changes in gene expression. In the present study, we examine a range of computational methods to calculate tPODs from HTTr data, using six data sets in which MCF7 cells cultured in two different media formulations were treated with a panel of 44 chemicals for 3 different exposure durations (6, 12, 24 hr). The tPOD calculation methods use data at the level of individual genes and gene set signatures, and compare data processed using the ToxCast Pipeline 2 (tcplfit2), BMDExpress and PLIER (Pathway Level Information ExtractoR). Methods were evaluated by comparing to in vitro PODs from a validated set of high-throughput screening (HTS) assays for a set of estrogenic compounds. Key findings include: (1) for a given chemical and set of experimental conditions, tPODs calculated by different methods can vary by several orders of magnitude; (2) tPODs are at least as sensitive to computational methods as to experimental conditions; (3) in comparison to an external reference set of PODs, some methods give generally higher values, principally PLIER and BMDExpress; and (4) the tPODs from HTTr in this one cell type are mostly higher than the overall PODs from a broad battery of targeted in vitro ToxCast assays, reflecting the need to test chemicals in multiple cell types and readout technologies for in vitro hazard screening.

The genomic landscape of mammal domestication might be orchestrated by selected transcription factors regulating brain and craniofacial development

Domestication transforms once wild animals into tamed animals that can be then exploited by humans. The process entails modifications in the body, cognition, and behavior that are essentially driven by differences in gene expression patterns. Although genetic and epigenetic mechanisms were shown to underlie such differences, less is known about the role exerted by trans-regulatory molecules, notably transcription factors (TFs) in domestication. In this paper, we conducted extensive in silico analyses aimed to clarify the TF landscape of mammal domestication. We first searched the literature, so as to establish a large list of genes selected with domestication in mammals. From this list, we selected genes experimentally demonstrated to exhibit TF functions. We also considered TFs displaying a statistically significant number of targets among the entire list of (domestication) selected genes. This workflow allowed us to identify 5 candidate TFs (SOX2, KLF4, MITF, NR3C1, NR3C2) that were further assessed in terms of biochemical and functional properties. We found that such TFs-of-interest related to mammal domestication are all significantly involved in the development of the brain and the craniofacial region, as well as the immune response and lipid metabolism. A ranking strategy, essentially based on a survey of protein-protein interactions datasets, allowed us to identify SOX2 as the main candidate TF involved in domestication-associated evolutionary changes. These findings should help to clarify the molecular mechanics of domestication and are of interest for future studies aimed to understand the behavioral and cognitive changes associated to domestication.

AMPK Associates with Chromatin and Phosphorylates the TAF-1 Subunit of the Transcription Initiation Complex to Regulate Histone Gene Expression in ALL Cells

The survival rates for relapsed/refractory acute lymphoblastic leukemia (ALL) remain poor. We and others have reported that ALL cells are vulnerable to conditions inducing energy/ER-stress mediated by AMP-activated protein kinase (AMPK). To identify the target genes directly regulated by AMPKα2, we performed genome-wide RNA-seq and ChIP-seq in CCRF-CEM (T-ALL) cells expressing HA-AMPKα2 (CN2) under normal and energy/metabolic stress conditions. CN2 cells show significantly altered AMPKα2 genomic binding and transcriptomic profile under metabolic stress conditions, including reduced histone gene expression. Proteomic analysis and in vitro kinase assays identified the TATA-Box-Binding Protein-Associated Factor 1 (TAF1) as a novel AMPKα2 substrate that downregulates histone gene transcription in response to energy/metabolic stress. Knockdown and knockout studies demonstrated that both AMPKα2 and TAF1 are required for histone gene expression. Mechanistically, upon activation, AMPKα2 phosphorylates TAF1 at Ser-1353 which impairs TAF1 interaction with RNA polymerase II (Pol II), leading to a compromised state of p-AMPKα2/p-TAF1/Pol II chromatin association and suppression of transcription. This mechanism was also observed in primary ALL cells and in vivo in NSG mice. Consequently, we uncovered a non-canonical function of AMPK that phosphorylates TAF1, both members of a putative chromatin-associated transcription complex that regulate histone gene expression, among others, in response to energy/metabolic stress.

IMPLICATIONS

Fully delineating the protein interactome by which AMPK regulates adaptive survival responses to energy/metabolic stress, either via epigenetic gene regulation or other mechanisms, will allow the rational development of strategies to overcome de novo or acquired resistance in ALL and other cancers.

Different PfEMP1-expressing Plasmodium falciparum variants induce divergent endothelial transcriptional responses during co-culture

The human malaria parasite Plasmodium falciparum is responsible for the majority of mortality and morbidity caused by malaria infection and differs from other human malaria species in the degree of accumulation of parasite-infected red blood cells in the microvasculature, known as cytoadherence or sequestration. In P. falciparum, cytoadherence is mediated by a protein called PfEMP1 which, due to its exposure to the host immune system, undergoes antigenic variation resulting in the expression of different PfEMP1 variants on the infected erythrocyte membrane. These PfEMP1s contain various combinations of adhesive domains, which allow for the differential engagement of a repertoire of endothelial receptors on the host microvasculature, with specific receptor usage associated with severe disease. We used a co-culture model of cytoadherence incubating human brain microvascular endothelial cells with erythrocytes infected with two parasite lines expressing different PfEMP1s that demonstrate different binding profiles to vascular endothelium. We determined the transcriptional profile of human brain microvascular endothelial cells (HBMEC) following different incubation periods with infected erythrocytes, identifying different transcriptional profiles of pathways previously found to be involved in the pathology of severe malaria, such as inflammation, apoptosis and barrier integrity, induced by the two PfEMP1 variants.

Neuronal-specific methylome and hydroxymethylome analysis reveal replicated and novel loci associated with alcohol use disorder

Alcohol use disorder (AUD) is a complex condition associated with adverse health consequences that affect millions of individuals worldwide. Epigenetic modifications, including DNA methylation (5mC), have been associated with AUD and other alcohol-related traits. Epigenome-wide association studies (EWAS) have identified differentially methylated genes associated with AUD in human peripheral and brain tissue. More recently, epigenetic studies of AUD have also evaluated DNA hydroxymethylation (5hmC) in the human brain. However, most of the epigenetic work in postmortem brain tissue has examined bulk tissue. In this study, we investigated neuronal-specific 5mC and 5hmC alterations at CpG sites associated with AUD in the human orbitofrontal cortex (OFC). Neuronal nuclei from the OFC were evaluated in 34 human postmortem brain samples (10 AUD, 24 non-AUD). Reduced representation oxidative bisulfite sequencing was used to assess 5mC and 5hmC at the genome-wide level. Differential 5mC and 5hmC were evaluated using the methylKit R package and significance was set at false discovery rate <0.05 and differential methylation >2. Functional enrichment analyses were performed and replication was evaluated replication in an independent dataset that assessed 5mC and 5hmC of AUD in bulk cortical tissue. We identified 417 5mC and 363 5hmC genome-wide significant differential CpG sites associated with AUD, with 59% in gene promoters. We also identified genes previously implicated in alcohol consumption, such as SYK, CHRM2, DNMT3A, and GATA4, for 5mC and GATA4, and GAD1, GATA4, DLX1 for 5hmC. Replication was observed for 28 CpG sites from a previous AUD 5mC and 5hmC study, including FOXP1. Lastly, GWAS enrichment analysis showed an association with AUD for differential 5mC genes. This study reveals neuronal-specific methylome and hydroxymethylome dysregulation associated with AUD. We replicated previous findings and identified novel associations with AUD for both 5mC and 5hmC marks within the OFC. Our findings provide new insights into the epigenomic dysregulation of AUD in the human brain.

PDE3A Is a Highly Expressed Therapy Target in Myxoid Liposarcoma

Liposarcomas (LPSs) are a heterogeneous group of malignancies that arise from adipose tissue. Although LPSs are among the most common soft-tissue sarcoma subtypes, precision medicine treatments are not currently available. To discover LPS-subtype-specific therapy targets, we investigated RNA sequenced transcriptomes of 131 clinical LPS tissue samples and compared the data with a transcriptome database that contained 20,218 samples from 95 healthy tissues and 106 cancerous tissue types. The identified genes were referred to the NCATS BioPlanet library with Enrichr to analyze upregulated signaling pathways. PDE3A protein expression was investigated with immunohistochemistry in 181 LPS samples, and PDE3A and SLFN12 mRNA expression with RT-qPCR were investigated in 63 LPS samples. Immunoblotting and cell viability assays were used to study LPS cell lines and their sensitivity to PDE3A modulators. We identified 97, 247, and 37 subtype-specific, highly expressed genes in dedifferentiated, myxoid, and pleomorphic LPS subtypes, respectively. Signaling pathway analysis revealed a highly activated hedgehog signaling pathway in dedifferentiated LPS, phospholipase c mediated cascade and insulin signaling in myxoid LPS, and pathways associated with cell proliferation in pleomorphic LPS. We discovered a strong association between high PDE3A expression and myxoid LPS, particularly in high-grade tumors. Moreover, myxoid LPS samples showed elevated expression levels of SLFN12 mRNA. In addition, PDE3A- and SLFN12-coexpressing LPS cell lines SA4 and GOT3 were sensitive to PDE3A modulators. Our results indicate that PDE3A modulators are promising drugs to treat myxoid LPS. Further studies are required to develop these drugs for clinical use.

Identification of new candidate drugs for primary Sjögren's syndrome using a drug repurposing transcriptomic approach

OBJECTIVES

To date, no immunomodulatory drug has demonstrated its efficacy in primary SS (pSS). We sought to analyse potential commonalities between pSS transcriptomic signatures and signatures of various drugs or specific knock-in or knock-down genes.

METHODS

Gene expression from peripheral blood samples of patients with pSS was compared with that of healthy controls in two cohorts and three public databases. In each of the five datasets, we analysed the 150 most up- and downregulated genes between pSS patients and controls with regard to the differentially expressed genes resulting from the biological action on nine cell lines of 2837 drugs, 2160 knock-in and 3799 knock-down genes in the Connectivity Map database.

RESULTS

We analysed 1008 peripheral blood transcriptomes from five independent studies (868 patients with pSS and 140 healthy controls). Eleven drugs could represent potential candidate drugs, with histone deacetylases and PI3K inhibitors among the most significantly associated. Twelve knock-in genes were associated with a pSS-like profile and 23 knock-down genes were associated with a pSS-revert profile. Most of those genes (28/35, 80%) were interferon-regulated.

CONCLUSION

This first drug repositioning transcriptomic approach in SS confirms the interest of targeting interferons and identifies histone deacetylases and PI3K inhibitors as potential therapeutic targets.

Deletion of Endogenous Neuregulin-4 Limits Adaptive Immunity During Interleukin-10 Receptor-Neutralizing Colitis

BACKGROUND

Growth factors are essential for maintenance of intestinal health. We previously showed that exogenous neuregulin-4 (NRG4) promotes colonocyte survival during cytokine challenge and is protective against acute models of intestinal inflammation. However, the function(s) of endogenous NRG4 are not well understood. Using NRG4-/- mice, we tested the role of endogenous NRG4 in models of colitis skewed toward either adaptive (interleukin-10 receptor [IL-10R] neutralization) or innate (dextran sulfate sodium [DSS]) immune responses.

METHODS

NRG4-/- and wild-type cage mate mice were subjected to chronic IL-10R neutralization colitis and acute DSS colitis. Disease was assessed by histological examination, inflammatory cytokine levels, fecal lipocalin-2 levels, and single cell mass cytometry immune cell profiling. Homeostatic gene alterations were evaluated by RNA sequencing analysis from colonic homogenates, with real-time quantitative polymerase chain reaction confirmation in both tissue and isolated epithelium.

RESULTS

During IL-10R neutralization colitis, NRG4-/- mice had reduced colonic inflammatory cytokine expression, histological damage, and colonic CD8+ T cell numbers vs wild-type cage mates. Conversely, in DSS colitis, NRG4-/- mice had elevated cytokine expression, fecal lipocalin-2 levels, and impaired weight recovery. RNA sequencing showed a loss of St3gal4, a sialyltransferase involved in immune cell trafficking, in NRG4-null colons, which was verified in both tissue and isolated epithelium. The regulation of St3gal4 by NRG4 was confirmed with ex vivo epithelial colon organoid cultures from NRG4-/- mice and by induction of St3gal4 in vivo following NRG4 treatment.

CONCLUSIONS

NRG4 regulates colonic epithelial ST3GAL4 and thus may allow for robust recruitment of CD8+ T cells during adaptive immune responses in colitis. On the other hand, NRG4 loss exacerbates injury driven by innate immune responses.

Integrated gene expression profiling and functional enrichment analyses to discover biomarkers and pathways associated with Guillain-Barré syndrome and autism spectrum disorder to identify new therapeutic targets

Guillain-Barré syndrome (GBS) is one of the most prominent and acute immune-mediated peripheral neuropathy, while autism spectrum disorders (ASD) are a group of heterogeneous neurodevelopmental disorders. The complete mechanism regarding the neuropathophysiology of these disorders is still ambiguous. Even after recent breakthroughs in molecular biology, the link between GBS and ASD remains a mystery. Therefore, we have implemented well-established bioinformatic techniques to identify potential biomarkers and drug candidates for GBS and ASD. 17 common differentially expressed genes (DEGs) were identified for these two disorders, which later guided the rest of the research. Common genes identified the protein-protein interaction (PPI) network and pathways associated with both disorders. Based on the PPI network, the constructed hub gene and module analysis network determined two common DEGs, namely CXCL9 and CXCL10, which are vital in predicting the top drug candidates. Furthermore, coregulatory networks of TF-gene and TF-miRNA were built to detect the regulatory biomolecules. Among drug candidates, imatinib had the highest docking and MM-GBSA score with the well-known chemokine receptor CXCR3 and remained stable during the 100 ns molecular dynamics simulation validated by the principal component analysis and the dynamic cross-correlation map. This study predicted the gene-based disease network for GBS and ASD and suggested prospective drug candidates. However, more in-depth research is required for clinical validation.Communicated by Ramaswamy H. Sarma.

Proteogenomic insights suggest druggable pathways in endometrial carcinoma

We characterized a prospective endometrial carcinoma (EC) cohort containing 138 tumors and 20 enriched normal tissues using 10 different omics platforms. Targeted quantitation of two peptides can predict antigen processing and presentation machinery activity, and may inform patient selection for immunotherapy. Association analysis between MYC activity and metformin treatment in both patients and cell lines suggests a potential role for metformin treatment in non-diabetic patients with elevated MYC activity. PIK3R1 in-frame indels are associated with elevated AKT phosphorylation and increased sensitivity to AKT inhibitors. CTNNB1 hotspot mutations are concentrated near phosphorylation sites mediating pS45-induced degradation of β-catenin, which may render Wnt-FZD antagonists ineffective. Deep learning accurately predicts EC subtypes and mutations from histopathology images, which may be useful for rapid diagnosis. Overall, this study identified molecular and imaging markers that can be further investigated to guide patient stratification for more precise treatment of EC.

SARS-CoV-2 niches in human placenta revealed by spatial transcriptomics

BACKGROUND

Functional placental niches are presumed to spatially separate maternal-fetal antigens and restrict the vertical transmission of pathogens. We hypothesized a high-resolution map of placental transcription could provide direct evidence for niche microenvironments with unique functions and transcription profiles.

METHODS

We utilized Visium Spatial Transcriptomics paired with H&E staining to generate 17,927 spatial transcriptomes. By integrating these spatial transcriptomes with 273,944 placental single-cell and single-nuclei transcriptomes, we generated an atlas composed of at least 22 subpopulations in the maternal decidua, fetal chorionic villi, and chorioamniotic membranes.

FINDINGS

Comparisons of placentae from uninfected healthy controls (n = 4) with COVID-19 asymptomatic (n = 4) and symptomatic (n = 5) infected participants demonstrated that severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) detection in syncytiotrophoblasts occurred in both the presence and the absence of maternal clinical disease. With spatial transcriptomics, we found that the limit of detection for SARS-CoV-2 was 1/7,000 cells, and placental niches without detectable viral transcripts were unperturbed. In contrast, niches with high SARS-CoV-2 transcript levels were associated with significant upregulation in pro-inflammatory cytokines and interferon-stimulated genes, altered metallopeptidase signaling (TIMP1), with coordinated shifts in macrophage polarization, histiocytic intervillositis, and perivillous fibrin deposition. Fetal sex differences in gene expression responses to SARS-CoV-2 were limited, with confirmed mapping limited to the maternal decidua in males.

CONCLUSIONS

High-resolution placental transcriptomics with spatial resolution revealed dynamic responses to SARS-CoV-2 in coordinate microenvironments in the absence and presence of clinically evident disease.

FUNDING

This work was supported by the NIH (R01HD091731 and T32-HD098069), NSF (2208903), the Burroughs Welcome Fund and the March of Dimes Preterm Birth Research Initiatives, and a Career Development Award from the American Society of Gene and Cell Therapy.

Exploring molecular mechanisms underlying the pathophysiological association between knee osteoarthritis and sarcopenia

Objectives

Accumulating evidence indicates a strong link between knee osteoarthritis (KOA) and sarcopenia. However, the mechanisms involved have not yet been elucidated. This study primarily aims to explore the molecular mechanisms that explain the connection between these 2 disorders.

Methods

The gene expression profiles for KOA and sarcopenia were obtained from the Gene Expression Omnibus database, specifically from GSE55235, GSE169077, and GSE1408. Various bioinformatics techniques were employed to identify and analyze common differentially expressed genes (DEGs) across the 3 datasets. The techniques involved the analysis of Gene Ontology and pathways to enhance understanding, examining protein-protein interaction (PPI) networks, and identifying hub genes. In addition, we constructed the network of interactions between transcription factors (TFs) and genes, the co-regulatory network of TFs and miRNAs for hub genes, and predicted potential drugs.

Results

In total, 14 common DEGs were found between KOA and sarcopenia. Detailed information on biological processes and signaling pathways of common DEGs was obtained through enrichment analysis. After performing PPI network analysis, we discovered 4 hub genes (FOXO3, BCL6, CDKN1A, and CEBPB). Subsequently, we developed coregulatory networks for these hub genes involving TF-gene and TF-miRNA interactions. Finally, we identified 10 potential chemical compounds.

Conclusions

By conducting bioinformatics analysis, our study has successfully identified common gene interaction networks between KOA and sarcopenia. The potential of these findings to offer revolutionary understanding into the common development of these 2 conditions could lead to the identification of valuable targets for therapy.

A co-formulation of interferons alpha2b and gamma distinctively targets cell cycle in the glioblastoma-derived cell line U-87MG

BACKGROUND

HeberFERON is a co-formulation of α2b and γ interferons, based on their synergism, which has shown its clinical superiority over individual interferons in basal cell carcinomas. In glioblastoma (GBM), HeberFERON has displayed promising preclinical and clinical results. This led us to design a microarray experiment aimed at identifying the molecular mechanisms involved in the distinctive effect of HeberFERON compared to the individual interferons in U-87MG model.

METHODS

Transcriptional expression profiling including a control (untreated) and three groups receiving α2b-interferon, γ-interferon and HeberFERON was performed using an Illumina HT-12 microarray platform. Unsupervised methods for gene and sample grouping, identification of differentially expressed genes, functional enrichment and network analysis computational biology methods were applied to identify distinctive transcription patterns of HeberFERON. Validation of most representative genes was performed by qPCR. For the cell cycle analysis of cells treated with HeberFERON for 24 h, 48 and 72 h we used flow cytometry.

RESULTS

The three treatments show different behavior based on the gene expression profiles. The enrichment analysis identified several mitotic cell cycle related events, in particular from prometaphase to anaphase, which are exclusively targeted by HeberFERON. The FOXM1 transcription factor network that is involved in several cell cycle phases and is highly expressed in GBMs, is significantly down regulated. Flow cytometry experiments corroborated the action of HeberFERON on the cell cycle in a dose and time dependent manner with a clear cellular arrest as of 24 h post-treatment. Despite the fact that p53 was not down-regulated, several genes involved in its regulatory activity were functionally enriched. Network analysis also revealed a strong relationship of p53 with genes targeted by HeberFERON. We propose a mechanistic model to explain this distinctive action, based on the simultaneous activation of PKR and ATF3, p53 phosphorylation changes, as well as its reduced MDM2 mediated ubiquitination and export from the nucleus to the cytoplasm. PLK1, AURKB, BIRC5 and CCNB1 genes, all regulated by FOXM1, also play central roles in this model. These and other interactions could explain a G2/M arrest and the effect of HeberFERON on the proliferation of U-87MG.

CONCLUSIONS

We proposed molecular mechanisms underlying the distinctive behavior of HeberFERON compared to the treatments with the individual interferons in U-87MG model, where cell cycle related events were highly relevant.

Thymic lymphoma detection in RORγ knockout mice using 5-hydroxymethylcytosine profiling of circulating cell-free DNA

A high incidence of thymic lymphoma has been noted in mice deficient of retinoid-related orphan receptor γ2 (RORγ2), which is required for differentiation of naïve CD4+ T cells into TH17 cells. Using a RORγ homozygous knockout (KO) mouse model of thymic lymphoma, we characterized this tumor progression and investigated the utility of 5-hydroxymethylcytosine (5hmC) signatures as a non-invasive circulating biomarker for early prediction of malignancy. No evidence for malignancy was noted in the wild-type mice, while primary thymic lymphoma with multi-organ metastasis was observed microscopically in 97% of the homozygous RORγ KO mice. The severity of thymic lymphoma was not age-dependent in the KO mice of 2 to 4 months old. Differential enrichment of 5hmC in thymic DNA and plasma cell-free DNA (cfDNA) was compared across different stages of tumor progression. Random forest modeling of plasma cfDNA achieved good predictivity (AUC = 0.74) in distinguishing early non-metastatic thymic lymphoma compared to cancer-free controls, while perfect predictivity was achieved with advanced multi-organ metastatic disease (AUC = 1.00). Lymphoid-specific genes involved in thymocyte selection during T cell development (Themis, Tox) were differentially enriched in both plasma and thymic tissue. This could help in differentiating thymic lymphoma from other tumors commonly detected in rodent carcinogenicity studies used in pharmaceutical drug development to inform human malignancy risk. Overall, these results provide a proof-of-concept for using circulating cfDNA profiles in rodent carcinogenicity studies for early risk assessment of novel pharmaceutical targets.

Gene expression profiling in whole blood stimulated ex vivo with lipopolysaccharide as a tool to predict post-stroke depressive symptoms: Proof-of-concept study

Prediction of post-stroke depressive symptoms (DSs) is challenging in patients without a history of depression. Gene expression profiling in blood cells may facilitate the search for biomarkers. The use of an ex vivo stimulus to the blood helps to reveal differences in gene profiles by reducing variation in gene expression. We conducted a proof-of-concept study to determine the usefulness of gene expression profiling in lipopolysaccharide (LPS)-stimulated blood for predicting post-stroke DS. Out of 262 enrolled patients with ischemic stroke, we included 96 patients without a pre-stroke history of depression and not taking any anti-depressive medication before or during the first 3 months after stroke. We assessed DS at 3 months after stroke using the Patient Health Questionnaire-9. We used RNA sequencing to determine the gene expression profile in LPS-stimulated blood samples taken on day 3 after stroke. We constructed a risk prediction model using a principal component analysis combined with logistic regression. We diagnosed post-stroke DS in 17.7% of patients. Expression of 510 genes differed between patients with and without DS. A model containing 6 genes (PKM, PRRC2C, NUP188, CHMP3, H2AC8, NOP10) displayed very good discriminatory properties (area under the curve: 0.95) with the sensitivity of 0.94 and specificity of 0.85. Our results suggest the potential utility of gene expression profiling in whole blood stimulated with LPS for predicting post-stroke DS. This method could be useful for searching biomarkers of post-stroke depression.

Programmed cell death 1 genetic variant and liver damage in nonalcoholic fatty liver disease

BACKGROUND AND AIMS

Programmed cell death 1/programmed cell death-ligand 1 (PD-1/PDL-1) axis has been reported to modulate liver inflammation and progression to hepatocellular carcinoma (HCC) in patients with nonalcoholic fatty liver disease (NAFLD). Here, we examined whether the PDCD1 variation is associated with NAFLD severity in individuals with liver biopsy.

METHODS

We examined the impact of PDCD1 gene variants on HCC, as robust severe liver disease phenotype in UK Biobank participants. The strongest genetic association with the rs13023138 G>C variation was subsequently tested for association with liver damage in 2889 individuals who underwent liver biopsy for suspected nonalcoholic steatohepatitis (NASH). Hepatic transcriptome was examined by RNA-Seq in a subset of NAFLD individuals (n = 121). Transcriptomic and deconvolution analyses were performed to identify biological pathways modulated by the risk allele.

RESULTS

The rs13023138 C>G showed the most robust association with HCC in UK Biobank (p = 5.28E-4, OR = 1.32, 95% CI [1.1, 1.5]). In the liver biopsy cohort, rs13023138 G allele was independently associated with severe steatosis (OR 1.17, 95% CI 1.02-1.34; p = .01), NASH (OR 1.22, 95% CI 1.09-1.37; p < .001) and advanced fibrosis (OR 1.26, 95% CI 1.06-1.50; p = .007). At deconvolution analysis, rs13023138 G>C allele was linked to higher hepatic representation of M1 macrophages, paralleled by upregulation of pathways related to inflammation and higher expression of CXCR6.

CONCLUSIONS

The PDCD1 rs13023138 G allele was associated with HCC development in the general population and with liver disease severity in patients at high risk of NASH.

Drug approval prediction based on the discrepancy in gene perturbation effects between cells and humans

BACKGROUND

Poor translation between in vitro and clinical studies due to the cells/humans discrepancy in drug target perturbation effects leads to safety failures in clinical trials, thus increasing drug development costs and reducing patients' life quality. Therefore, developing a predictive model for drug approval considering the cells/humans discrepancy is needed to reduce drug attrition rates in clinical trials.

METHODS

Our machine learning framework predicts drug approval in clinical trials based on the cells/humans discrepancy in drug target perturbation effects. To evaluate the discrepancy to predict drug approval (1404 approved and 1070 unapproved drugs), we analysed CRISPR-Cas9 knockout and loss-of-function mutation rate-based gene perturbation effects on cells and humans, respectively. To validate the risk of drug targets with the cells/humans discrepancy, we examined the targets of failed and withdrawn drugs due to safety problems.

FINDINGS

Drug approvals in clinical trials were correlated with the cells/humans discrepancy in gene perturbation effects. Genes tolerant to perturbation effects on cells but intolerant to those on humans were associated with failed drug targets. Furthermore, genes with the cells/humans discrepancy were related to drugs withdrawn due to severe side effects. Motivated by previous studies assessing drug safety through chemical properties, we improved drug approval prediction by integrating chemical information with the cells/humans discrepancy.

INTERPRETATION

The cells/humans discrepancy in gene perturbation effects facilitates drug approval prediction and explains drug safety failures in clinical trials.

FUNDING

S.K. received grants from the Korean National Research Foundation (2021R1A2B5B01001903 and 2020R1A6A1A03047902) and IITP (2019-0-01906, Artificial Intelligence Graduate School Program, POSTECH).

Serine arginine-rich splicing factor (SRSF7) cooperates with the histone methyltransferase KMT5a to promote the type I interferon response via transcriptional activation of IRF7

Tight regulation of macrophage immune gene expression is required to fight infection without risking harmful inflammation. The contribution of RNA binding proteins (RBPs) to shaping the macrophage response to pathogens remains poorly understood. Transcriptomic analysis revealed that a member of the serine/arginine-rich (SR) family of mRNA processing factors, SRSF7, is required for optimal expression of a cohort of interferon stimulated genes (ISGs) in macrophages. Using genetic and biochemical assays, we discovered that in addition to its canonical role in regulating alternative splicing, SRSF7 drives transcription of interferon regulatory transcription factor 7 (IRF7) to promote antiviral immunity. At the Irf7 promoter, SRSF7 maximizes STAT1 transcription factor binding and RNA polymerase II elongation via cooperation with the H4K20me1 histone methyltransferase KMT5a (SET8). These studies define an unorthodox role for an SR protein in activating transcription and reveal an unappreciated RNA binding protein-chromatin network that orchestrates macrophage antiviral gene expression.

Profiling of miRNAs and their interfering targets in peripheral blood mononuclear cells from patients with chronic myeloid leukaemia

Introduction

MicroRNAs may be implicated in the acquisition of drug resistance in chronic myeloid leukemia as they regulate the expression of not only BCR-ABL1 but also genes associated with the activation of drug transfer proteins or essential signaling pathways.

Methods

To understand the impact of specifically expressed miRNAs in chronic myeloid leukemia and their target genes, we collected peripheral blood mononuclear cells (PBMC) from patients diagnosed with chronic myeloid leukemia (CML) and healthy donors to determine whole miRNA expression by small RNA sequencing and screened out 31 differentially expressed microRNAs (DE-miRNAs) with high expression. With the utilization of miRNA set enrichment analysis tools, we present here a comprehensive analysis of the relevance of DE-miRNAs to disease and biological function. Furthermore, the literature-based miRNA-target gene database was used to analyze the overall target genes of the DE-miRNAs and to define their associated biological responses. We further integrated DE-miRNA target genes to identify CML miRNA targeted gene signature singscore (CMTGSS) and used gene-set enrichment analysis (GSEA) to analyze the correlation between CMTGSS and Hallmark gene-sets in PBMC samples from clinical CML patients. Finally, the association of CMTGSS stratification with multiple CML cell lineage gene sets was validated in PBMC samples from CML patients using GSEA.

Results

Although individual miRNAs have been reported to have varying degrees of impact on CML, overall, our results show that abnormally upregulated miRNAs are associated with apoptosis and aberrantly downregulated miRNAs are associated with cell cycle. The clinical database shows that our defined DE-miRNAs are associated with the prognosis of CML patients. CMTGSS-based stratification analysis presented a tendency for miRNAs to affect cell differentiation in the blood microenvironment.

Conclusion

Collectively, this study defined differentially expressed miRNAs by miRNA sequencing from clinical samples and comprehensively analyzed the biological functions of the differential miRNAs in association with the target genes. The analysis of the enrichment of specific myeloid differentiated cells and immune cells also suggests the magnitude and potential targets of differentially expressed miRNAs in the clinical setting. It helps us to make links between the different results obtained from the multi-faceted studies to provide more potential research directions.

Gene expression profiling unveils the temporal dynamics of CIGB-300-regulated transcriptome in AML cell lines

BACKGROUND

Protein kinase CK2 activity is implicated in the pathogenesis of various hematological malignancies like Acute Myeloid Leukemia (AML) that remains challenging concerning treatment. This kinase has emerged as an attractive molecular target in therapeutic. Antitumoral peptide CIGB-300 blocks CK2 phospho-acceptor sites on their substrates but it also binds to CK2α catalytic subunit. Previous proteomic and phosphoproteomic experiments showed molecular and cellular processes with relevance for the peptide action in diverse AML backgrounds but earlier transcriptional level events might also support the CIGB-300 anti-leukemic effect. Here we used a Clariom S HT assay for gene expression profiling to study the molecular events supporting the anti-leukemic effect of CIGB-300 peptide on HL-60 and OCI-AML3 cell lines.

RESULTS

We found 183 and 802 genes appeared significantly modulated in HL-60 cells at 30 min and 3 h of incubation with CIGB-300 for p < 0.01 and FC >  = │1.5│, respectively; while 221 and 332 genes appeared modulated in OCI-AML3 cells. Importantly, functional enrichment analysis evidenced that genes and transcription factors related to apoptosis, cell cycle, leukocyte differentiation, signaling by cytokines/interleukins, and NF-kB, TNF signaling pathways were significantly represented in AML cells transcriptomic profiles. The influence of CIGB-300 on these biological processes and pathways is dependent on the cellular background, in the first place, and treatment duration. Of note, the impact of the peptide on NF-kB signaling was corroborated by the quantification of selected NF-kB target genes, as well as the measurement of p50 binding activity and soluble TNF-α induction. Quantification of CSF1/M-CSF and CDKN1A/P21 by qPCR supports peptide effects on differentiation and cell cycle.

CONCLUSIONS

We explored for the first time the temporal dynamics of the gene expression profile regulated by CIGB-300 which, along with the antiproliferative mechanism, can stimulate immune responses by increasing immunomodulatory cytokines. We provided fresh molecular clues concerning the antiproliferative effect of CIGB-300 in two relevant AML backgrounds.

Functional specialization of short-lived and long-lived macrophage subsets in human tonsils

Macrophages play a central role in tissue homeostasis and host defense. However, the properties of human macrophages in non-diseased tissues remain poorly understood. Here, we characterized human tonsil macrophages and identified three subsets with distinct phenotype, transcriptome, life cycle, and function. CD36hi macrophages were related to monocytes, while CD36lo macrophages showed features of embryonic origin and CD36int macrophages had a mixed profile. scRNA-seq on non-human primate tonsils showed that monocyte recruitment did not pre-exist an immune challenge. Functionally, CD36hi macrophages were specialized for stimulating T follicular helper cells, by producing Activin A. Combining reconstruction of ligand-receptor interactions and functional assays, we identified stromal cell-derived TNF-α as an inducer of Activin A secretion. However, only CD36hi macrophages were primed for Activin A expression, via the activity of IRF1. Our results provide insight into the heterogeneity of human lymphoid organ macrophages and show that tonsil CD36hi macrophage specialization is the result of both intrinsic features and interaction with stromal cells.

The effect of inhibition of receptor tyrosine kinase AXL on DNA damage response in ovarian cancer

AXL is a receptor tyrosine kinase that is often overexpressed in cancers. It contributes to pathophysiology in cancer progression and therapeutic resistance, making it an emerging therapeutic target. The first-in-class AXL inhibitor bemcentinib (R428/BGB324) has been granted fast track designation by the U.S. Food and Drug Administration (FDA) in STK11-mutated advanced metastatic non-small cell lung cancer and was also reported to show selective sensitivity towards ovarian cancers (OC) with a Mesenchymal molecular subtype. In this study, we further explored AXL's role in mediating DNA damage responses by using OC as a disease model. AXL inhibition using R428 resulted in the increase of DNA damage with the concurrent upregulation of DNA damage response signalling molecules. Furthermore, AXL inhibition rendered cells more sensitive to the inhibition of ATR, a crucial mediator for replication stress. Combinatory use of AXL and ATR inhibitors showed additive effects in OC. Through SILAC co-immunoprecipitation mass spectrometry, we identified a novel binding partner of AXL, SAM68, whose loss in OC cells harboured phenotypes in DNA damage responses similar to AXL inhibition. In addition, AXL- and SAM68-deficiency or R428 treatment induced elevated levels of cholesterol and upregulated genes in the cholesterol biosynthesis pathway. There might be a protective role of cholesterol in shielding cancer cells against DNA damage induced by AXL inhibition or SMA68 deficiency.

Myeloid cell influx into the colonic epithelium is associated with disease severity and non-response to anti-Tumor Necrosis Factor Therapy in patients with Ulcerative Colitis

Ulcerative colitis (UC) is an idiopathic chronic inflammatory disease of the colon with sharply rising global prevalence. Dysfunctional epithelial compartment (EC) dynamics are implicated in UC pathogenesis although EC-specific studies are sparse. Applying orthogonal high-dimensional EC profiling to a Primary Cohort (PC; n=222), we detail major epithelial and immune cell perturbations in active UC. Prominently, reduced frequencies of mature BEST4+OTOP2+ absorptive and BEST2+WFDC2+ secretory epithelial enterocytes were associated with the replacement of homeostatic, resident TRDC+KLRD1+HOPX+ γδ+ T cells with RORA+CCL20+S100A4+ TH17 cells and the influx of inflammatory myeloid cells. The EC transcriptome (exemplified by S100A8, HIF1A, TREM1, CXCR1) correlated with clinical, endoscopic, and histological severity of UC in an independent validation cohort (n=649). Furthermore, therapeutic relevance of the observed cellular and transcriptomic changes was investigated in 3 additional published UC cohorts (n=23, 48 and 204 respectively) to reveal that non-response to anti-Tumor Necrosis Factor (anti-TNF) therapy was associated with EC related myeloid cell perturbations. Altogether, these data provide high resolution mapping of the EC to facilitate therapeutic decision-making and personalization of therapy in patients with UC.

Enhanced chondrogenic differentiation of iPS cell-derived mesenchymal stem/stromal cells via neural crest cell induction for hyaline cartilage repair

Background: To date, there is no effective long-lasting treatment for cartilage tissue repair. Primary chondrocytes and mesenchymal stem/stromal cells are the most commonly used cell sources in regenerative medicine. However, both cell types have limitations, such as dedifferentiation, donor morbidity, and limited expansion. Here, we report a stepwise differentiation method to generate matrix-rich cartilage spheroids from induced pluripotent stem cell-derived mesenchymal stem/stromal cells (iMSCs) via the induction of neural crest cells under xeno-free conditions. Methods: The genes and signaling pathways regulating the chondrogenic susceptibility of iMSCs generated under different conditions were studied. Enhanced chondrogenic differentiation was achieved using a combination of growth factors and small-molecule inducers. Results: We demonstrated that the use of a thienoindazole derivative, TD-198946, synergistically improves chondrogenesis in iMSCs. The proposed strategy produced controlled-size spheroids and increased cartilage extracellular matrix production with no signs of dedifferentiation, fibrotic cartilage formation, or hypertrophy in vivo. Conclusion: These findings provide a novel cell source for stem cell-based cartilage repair. Furthermore, since chondrogenic spheroids have the potential to fuse within a few days, they can be used as building blocks for biofabrication of larger cartilage tissues using technologies such as the Kenzan Bioprinting method.

A Descriptive Whole-Genome Transcriptomics Study in a Stem Cell-Based Tool Predicts Multiple Tissue-Specific Beneficial Potential and Molecular Targets of Carnosic Acid

Carnosic acid (CA) is a phenolic diterpene widely distributed in herbal plants, rosemary and sage. Although its medicinal properties, such as antioxidant, antimicrobial, and neuroprotective effects, have been well-documented, its relevant biochemical processes and molecular targets have not been fully explored yet. In the present study, we conducted an untargeted whole-genome transcriptomics analysis to investigate CA-induced early biological and molecular events in human amniotic epithelial stem cells (hAESCs) with the aim of exploring its multiple tissue-specific functionalities and potential molecular targets. We found that seven days of CA treatment in hAESCs could induce mesoderm-lineage-specific differentiation. Tissue enrichment analysis revealed that CA significantly enriched lateral plate mesoderm-originated cardiovascular and adipose tissues. Further tissue-specific PPI analysis and kinase and transcription factor enrichment analyses identified potential upstream regulators and molecular targets of CA in a tissue-specific manner. Gene ontology enrichment analyses revealed the metabolic, antioxidant, and antifibrotic activities of CA. Altogether, our comprehensive whole-genome transcriptomics analyses offer a thorough understanding of the possible underlying molecular mechanism of CA.

CD103+ regulatory T cells underlie resistance to radio-immunotherapy and impair CD8+ T cell activation in glioblastoma

Glioblastomas are aggressive primary brain tumors with an inherent resistance to T cell-centric immunotherapy due to their low mutational burden and immunosuppressive tumor microenvironment. Here we report that fractionated radiotherapy of preclinical glioblastoma models induce a tenfold increase in T cell content. Orthogonally, spatial imaging mass cytometry shows T cell enrichment in human recurrent tumors compared with matched primary glioblastoma. In glioblastoma-bearing mice, α-PD-1 treatment applied at the peak of T cell infiltration post-radiotherapy results in a modest survival benefit compared with concurrent α-PD-1 administration. Following α-PD-1 therapy, CD103⁺ regulatory T cells (Tregs) with upregulated lipid metabolism accumulate in the tumor microenvironment, and restrain immune checkpoint blockade response by repressing CD8⁺ T cell activation. Treg targeting elicits tertiary lymphoid structure formation, enhances CD4⁺ and CD8⁺ T cell frequency and function and unleashes radio-immunotherapeutic efficacy. These results support the rational design of therapeutic regimens limiting the induction of immunosuppressive feedback pathways in the context of T cell immunotherapy in glioblastoma.

Enhancer hijacking at the ARHGAP36 locus is associated with connective tissue to bone transformation

Heterotopic ossification is a disorder caused by abnormal mineralization of soft tissues in which signaling pathways such as BMP, TGFβ and WNT are known key players in driving ectopic bone formation. Identifying novel genes and pathways related to the mineralization process are important steps for future gene therapy in bone disorders. In this study, we detect an inter-chromosomal insertional duplication in a female proband disrupting a topologically associating domain and causing an ultra-rare progressive form of heterotopic ossification. This structural variant lead to enhancer hijacking and misexpression of ARHGAP36 in fibroblasts, validated here by orthogonal in vitro studies. In addition, ARHGAP36 overexpression inhibits TGFβ, and activates hedgehog signaling and genes/proteins related to extracellular matrix production. Our work on the genetic cause of this heterotopic ossification case has revealed that ARHGAP36 plays a role in bone formation and metabolism, outlining first details of this gene contributing to bone-formation and -disease.

qHTSWaterfall: 3-dimensional visualization software for quantitative high-throughput screening (qHTS) data

High throughput screening (HTS) is widely used in drug discovery and chemical biology to identify and characterize agents having pharmacologic properties often by evaluation of large chemical libraries. Standard HTS data can be simply plotted as an x-y graph usually represented as % activity of a compound tested at a single concentration vs compound ID, whereas quantitative HTS (qHTS) data incorporates a third axis represented by concentration. By virtue of the additional data points arising from the compound titration and the incorporation of logistic fit parameters that define the concentration-response curve, such as EC50 and Hill slope, qHTS data has been challenging to display on a single graph. Here we provide a flexible solution to the rapid plotting of complete qHTS data sets to produce a 3-axis plot we call qHTS Waterfall Plots. The software described here can be generally applied to any 3-axis dataset and is available as both an R package and an R shiny application.