Multiplatform metabolomic interlaboratory study of a whole human stool candidate reference material from omnivore and vegan donors.

INTRODUCTION: Human metabolomics has made significant strides in understanding metabolic changes and their implications for human health, with promising applications in diagnostics and treatment, particularly regarding the gut microbiome. However, progress is hampered by issues with data comparability and reproducibility across studies, limiting the translation of these discoveries into practical applications. OBJECTIVES: This study aims to evaluate the fit-for-purpose of a suite of human stool samples as potential candidate reference materials (RMs) and assess the state of the field regarding harmonizing gut metabolomics measurements. METHODS: An interlaboratory study was conducted with 18 participating institutions. The study allowed for the use of preferred analytical techniques, including liquid chromatography-mass spectrometry (LC-MS), gas chromatography-mass spectrometry (GC-MS), and nuclear magnetic resonance (NMR). RESULTS: Different laboratories used various methods and analytical platforms to identify the metabolites present in human stool RM samples. The study found a 40% to 70% recurrence in the reported top 20 most abundant metabolites across the four materials. In the full annotation list, the percentage of metabolites reported multiple times after nomenclature standardization was 36% (LC-MS), 58% (GC-MS) and 76% (NMR). Out of 9,300 unique metabolites, only 37 were reported across all three measurement techniques. CONCLUSION: This collaborative exercise emphasized the broad chemical survey possible with multi-technique approaches. Community engagement is essential for the evaluation and characterization of common materials designed to facilitate comparability and ensure data quality underscoring the value of determining current practices, challenges, and progress of a field through interlaboratory studies.

Allelic effects on KLHL17 expression likely mediated by JunB/D underlie a PDAC GWAS signal at chr1p36.33.

Pancreatic Ductal Adenocarcinoma (PDAC) is the third leading cause of cancer-related deaths in the U.S. Both rare and common germline variants contribute to PDAC risk. Here, we fine-map and functionally characterize a common PDAC risk signal at 1p36.33 (tagged by rs13303010) identified through a genome wide association study (GWAS). One of the fine-mapped SNPs, rs13303160 (r2=0.93 in 1000G EUR samples, OR=1.23, P value=2.74x10-9) demonstrated allele-preferential gene regulatory activity in vitro and allele-preferential binding of JunB and JunD in vitro and in vivo. Expression Quantitative Trait Locus (eQTL) analysis identified KLHL17 as a likely target gene underlying the signal. Proteomic analysis identified KLHL17 as a member of the Cullin-E3 ubiquitin ligase complex in PDAC-derived cells. In silico differential gene expression analysis of the GTExv8 pancreas data suggested an association between lower KLHL17 (risk associated) and pro-inflammatory pathways. We hypothesize that KLHL17 may mitigate inflammation by recruiting pro-inflammatory proteins for ubiquitination and degradation thereby influencing PDAC risk.

Large-scale multi-omic analysis identifies noncoding somatic driver mutations and nominates ZFP36L2 as a driver gene for pancreatic ductal adenocarcinoma.

Identification of somatic driver mutations in the noncoding genome remains challenging. To comprehensively characterize noncoding driver mutations for pancreatic ductal adenocarcinoma (PDAC), we first created genome-scale maps of accessible chromatin regions (ACRs) and histone modification marks (HMMs) in pancreatic cell lines and purified pancreatic acinar and duct cells. Integration with whole-genome mutation calls from 506 PDACs revealed 314 ACRs/HMMs significantly enriched with 3,614 noncoding somatic mutations (NCSMs). Functional assessment using massively parallel reporter assays (MPRA) identified 178 NCSMs impacting reporter activity (19.45% of those tested). Focused luciferase validation confirmed negative effects on gene regulatory activity for NCSMs near CDKN2A and ZFP36L2. For the latter, CRISPR interference (CRISPRi) further identified ZFP36L2 as a target gene (16.0 - 24.0% reduced expression, P = 0.023-0.0047) with disrupted KLF9 binding likely mediating the effect. Our integrative approach provides a catalog of potentially functional noncoding driver mutations and nominates ZFP36L2 as a PDAC driver gene.

Targeted sulfur(VI) fluoride exchange-mediated covalent modification of a tyrosine residue in the catalytic pocket of tyrosyl-DNA phosphodiesterase 1.

Developing effective inhibitors of the DNA repair enzyme tyrosyl-DNA phosphodiesterase 1 (TDP1) has been challenging because of the enzyme shallow catalytic pocket and non-specific substrate binding interactions. Recently, we discovered a quinolone-binding hot spot in TDP1's active site proximal to the evolutionary conserved Y204 and F259 residues that position DNA. Sulfur (VI) fluoride exchange (SuFEx) is a biocompatible click chemistry reaction that enables acylation of protein residues, including tyrosine. Selective protein modifications can provide insights into the biological roles of proteins and inform ligand design. As we report herein, we used SuFEx chemistries to prepare covalent TDP1-bound binders showing site-specific covalent bonds with Y204. Our work presents the first application of SuFEx chemistries to TDP1 ligands. It validates the ability to covalently modify specific TDP1 residues by designed targeting and adds to the chemical biology resource toolbox for studying TDP1.

A patient-derived cell model for malignant transformation in IDH-mutant glioma.

Malignant transformation (MT) is commonly seen in IDH-mutant gliomas. There has been a growing research interest in revealing its underlying mechanisms and intervening prior to MT at the early stages of the transforming process. Here we established a unique pair of matched 3D cell models: 403L, derived from a low-grade glioma (LGG), and 403H, derived from a high-grade glioma (HGG), by utilizing IDH-mutant astrocytoma samples from the same patient when the tumor was diagnosed as WHO grade 2 (tumor mutational burden (TMB) of 3.96/Mb) and later as grade 4 (TMB of 70.07/Mb), respectively. Both cell models were authenticated to a patient's sample retaining endogenous expression of IDH1 R132H. DNA methylation profiles of the parental tumors referred to LGG and HGG IDH-mutant glioma clusters. The immunopositivity of SOX2, NESTIN, GFAP, OLIG2, and beta 3-Tubulin suggested the multilineage potential of both models. 403H was more prompt to cell invasion and developed infiltrative HGG in vivo. The differentially expressed genes (DEGs) from the RNA sequencing analysis revealed the tumor invasion and aggressiveness related genes exclusively upregulated in the 403H model. Pathway analysis showcased an enrichment of genes associated with epithelial-mesenchymal transition (EMT) and Notch signaling pathways in 403H and 403L, respectively. Mass spectrometry-based targeted metabolomics and hyperpolarized (HP) 1-13C pyruvate in-cell NMR analyses demonstrated significant alterations in the TCA cycle and fatty acid metabolism. Citrate, glutamine, and 2-HG levels were significantly higher in 403H. To our knowledge, this is the first report describing the development of a matched pair of 3D patient-derived cell models representative of MT and temozolomide (TMZ)-induced hypermutator phenotype (HMP) in IDH-mutant glioma, providing insights into genetic and metabolic changes during MT/HMP. This novel in vitro model allows further investigation of the mechanisms of MT at the cellular level.

Transfer RNA acetylation regulates in vivo mammalian stress signaling.

Transfer RNA (tRNA) modifications are crucial for protein synthesis, but their position-specific physiological roles remain poorly understood. Here we investigate the impact of N4-acetylcytidine (ac4C), a highly conserved tRNA modification, using a Thumpd1 knockout mouse model. We find that loss of Thumpd1-dependent tRNA acetylation leads to reduced levels of tRNALeu, increased ribosome stalling, and activation of eIF2α phosphorylation. Thumpd1 knockout mice exhibit growth defects and sterility. Remarkably, concurrent knockout of Thumpd1 and the stress-sensing kinase Gcn2 causes penetrant postnatal lethality, indicating a critical genetic interaction. Our findings demonstrate that a modification restricted to a single position within type II cytosolic tRNAs can regulate ribosome-mediated stress signaling in mammalian organisms, with implications for our understanding of translation control as well as therapeutic interventions.

Tyrosine phosphorylation of both STAT5A and STAT5B is necessary for maximal IL-2 signaling and T cell proliferation.

Cytokine-mediated STAT5 protein activation is vital for lymphocyte development and function. In vitro tyrosine phosphorylation of a C-terminal tyrosine is critical for activation of STAT5A and STAT5B; however, the importance of STAT5 tyrosine phosphorylation in vivo has not been assessed. Here we generate Stat5a and Stat5b tyrosine-to-phenylalanine mutant knockin mice and find they have greatly reduced CD8+ T-cell numbers and profoundly diminished IL-2-induced proliferation of these cells, and this correlates with reduced induction of Myc, pRB, a range of cyclins and CDKs, and a partial G1→S phase-transition block. These mutant CD8+ T cells also exhibit decreased IL-2-mediated activation of pERK and pAKT, which we attribute in part to diminished expression of IL-2Rß and IL-2Rγ. Our findings thus demonstrate that tyrosine phosphorylation of both STAT5A and STAT5B is essential for maximal IL-2 signaling. Moreover, our transcriptomic and proteomic analyses elucidate the molecular basis of the IL-2-induced proliferation of CD8+ T cells.

Identification of molecular determinants of gene-specific bursting patterns by high-throughput imaging screens.

Stochastic transcriptional bursting is a universal property of active genes. While different genes exhibit distinct bursting patterns, the molecular mechanisms for gene-specific stochastic bursting are largely unknown. We have developed and applied a high-throughput-imaging based screening strategy to identify cellular factors and molecular mechanisms that determine the bursting behavior of human genes. Focusing on epigenetic regulators, we find that protein acetylation is a strong acute modulator of burst frequency, burst size and heterogeneity of bursting. Acetylation globally affects the Off-time of genes but has gene-specific effects on the On-time. Yet, these effects are not strongly linked to promoter acetylation, which do not correlate with bursting properties, and forced promoter acetylation has variable effects on bursting. Instead, we demonstrate acetylation of the Integrator complex as a key determinant of gene bursting. Specifically, we find that elevated Integrator acetylation decreases bursting frequency. Taken together our results suggest a prominent role of non-histone proteins in determining gene bursting properties, and they identify histone-independent acetylation of a transcription cofactor as an allosteric modulator of bursting via a far-downstream bursting checkpoint.

Engineering Tumor Stroma Morphogenesis Using Dynamic Cell-Matrix Spheroid Assembly.

The tumor microenvironment consists of resident tumor cells organized within a compositionally diverse, three-dimensional (3D) extracellular matrix (ECM) network that cannot be replicated in vitro using bottom-up synthesis. We report a new self-assembly system to engineer ECM-rich 3D MatriSpheres wherein tumor cells actively organize and concentrate microgram quantities of decellularized ECM dispersions which modulate cell phenotype. 3D colorectal cancer (CRC) MatriSpheres were created using decellularized small intestine submucosa (SIS) as an orthotopic ECM source that had greater proteomic homology to CRC tumor ECM than traditional ECM formulations such as Matrigel. SIS ECM was rapidly concentrated from its environment and assembled into ECM-rich 3D stroma-like regions by mouse and human CRC cell lines within 4-5 days via a mechanism that was rheologically distinct from bulk hydrogel formation. Both ECM organization and transcriptional regulation by 3D ECM cues affected programs of malignancy, lipid metabolism, and immunoregulation that corresponded with an in vivo MC38 tumor cell subpopulation identified via single cell RNA sequencing. This 3D modeling approach stimulates tumor specific tissue morphogenesis that incorporates the complexities of both cancer cell and ECM compartments in a scalable, spontaneous assembly process that may further facilitate precision medicine.

IRF4 requires ARID1A to establish plasma cell identity in multiple myeloma.

Multiple myeloma (MM) is an incurable plasma cell malignancy that exploits transcriptional networks driven by IRF4. We employ a multi-omics approach to discover IRF4 vulnerabilities, integrating functional genomics screening, spatial proteomics, and global chromatin mapping. ARID1A, a member of the SWI/SNF chromatin remodeling complex, is required for IRF4 expression and functionally associates with IRF4 protein on chromatin. Deleting Arid1a in activated murine B cells disrupts IRF4-dependent transcriptional networks and blocks plasma cell differentiation. Targeting SWI/SNF activity leads to rapid loss of IRF4-target gene expression and quenches global amplification of oncogenic gene expression by MYC, resulting in profound toxicity to MM cells. Notably, MM patients with aggressive disease bear the signature of SWI/SNF activity, and SMARCA2/4 inhibitors remain effective in immunomodulatory drug (IMiD)-resistant MM cells. Moreover, combinations of SWI/SNF and MEK inhibitors demonstrate synergistic toxicity to MM cells, providing a promising strategy for relapsed/refractory disease.

Genome-scale exon perturbation screens uncover exons critical for cell fitness.

CRISPR-Cas technology has transformed functional genomics, yet understanding of how individual exons differentially shape cellular phenotypes remains limited. Here, we optimized and conducted massively parallel exon deletion and splice-site mutation screens in human cell lines to identify exons that regulate cellular fitness. Fitness-promoting exons are prevalent in essential and highly expressed genes and commonly overlap with protein domains and interaction interfaces. Conversely, fitness-suppressing exons are enriched in nonessential genes, exhibiting lower inclusion levels, and overlap with intrinsically disordered regions and disease-associated mutations. In-depth mechanistic investigation of the screen-hit TAF5 alternative exon-8 revealed that its inclusion is required for assembly of the TFIID general transcription initiation complex, thereby regulating global gene expression output. Collectively, our orthogonal exon perturbation screens established a comprehensive repository of phenotypically important exons and uncovered regulatory mechanisms governing cellular fitness and gene expression.

Chemical tools to define and manipulate interferon-inducible Ubl protease USP18.

Ubiquitin-specific protease 18 (USP18) is a multifunctional cysteine protease primarily responsible for deconjugating interferon-inducible ubiquitin-like (Ubl) modifier ISG15 from protein substrates. Here, we report the design and synthesis of activity-based probes (ABPs) capable of selectively detecting USP18 activity over other ISG15 cross-reactive deubiquitinases (DUBs) by incorporating unnatural amino acids into the C-terminal tail of ISG15. Combining with a ubiquitin-based DUB ABP, the selective USP18 ABP is employed in a chemoproteomic screening platform to identify and assess inhibitors of DUBs including USP18. We further demonstrate that USP18 ABPs can be utilized to profile differential activities of USP18 in lung cancer cell lines, providing a strategy that will help define the activity-related landscape of USP18 in different disease states and unravel important (de)ISGylation-dependent biological processes.

Discovery of a Tunable Heterocyclic Electrophile 4-Chloro-pyrazolopyridine That Defines a Unique Subset of Ligandable Cysteines.

Electrophilic small molecules with novel reactivity are powerful tools that enable activity-based protein profiling and covalent inhibitor discovery. Here, we report a reactive heterocyclic scaffold, 4-chloro-pyrazolopyridine (CPzP) for selective modification of proteins via a nucleophilic aromatic substitution (SNAr) mechanism. Chemoproteomic profiling reveals that CPzPs engage cysteines within functionally diverse protein sites including ribosomal protein S5 (RPS5), inosine monophosphate dehydrogenase 2 (IMPDH2), and heat shock protein 60 (HSP60). Through the optimization of appended recognition elements, we demonstrate the utility of CPzP for covalent inhibition of prolyl endopeptidase (PREP) by targeting a noncatalytic active-site cysteine. This study suggests that the proteome reactivity of CPzPs can be modulated by both electronic and steric features of the ring system, providing a new tunable electrophile for applications in chemoproteomics and covalent inhibitor design.

Sphingolipid biosynthesis is essential for metabolic rewiring during TH17 cell differentiation.

T helper 17 (TH17) cells are implicated in autoimmune diseases, and several metabolic processes are shown to be important for their development and function. In this study, we report an essential role for sphingolipids synthesized through the de novo pathway in TH17 cell development. Deficiency of SPTLC1, a major subunit of serine palmitoyl transferase enzyme complex that catalyzes the first and rate-limiting step of de novo sphingolipid synthesis, impaired glycolysis in differentiating TH17 cells by increasing intracellular reactive oxygen species (ROS) through enhancement of nicotinamide adenine dinucleotide phosphate oxidase 2 activity. Increased ROS leads to impaired activation of mammalian target of rapamycin C1 and reduced expression of hypoxia-inducible factor 1-alpha and c-Myc-induced glycolytic genes. SPTLCI deficiency protected mice from developing experimental autoimmune encephalomyelitis and experimental T cell transfer colitis. Our results thus show a critical role for de novo sphingolipid biosynthetic pathway in shaping adaptive immune responses with implications in autoimmune diseases.

A structure-based designed small molecule depletes hRpn13Pru and a select group of KEN box proteins.

Proteasome subunit hRpn13 is partially proteolyzed in certain cancer cell types to generate hRpn13Pru by degradation of its UCHL5/Uch37-binding DEUBAD domain and retention of an intact proteasome- and ubiquitin-binding Pru domain. By using structure-guided virtual screening, we identify an hRpn13 binder (XL44) and solve its structure ligated to hRpn13 Pru by integrated X-ray crystallography and NMR to reveal its targeting mechanism. Surprisingly, hRpn13Pru is depleted in myeloma cells following treatment with XL44. TMT-MS experiments reveal a select group of off-targets, including PCNA clamp-associated factor PCLAF and ribonucleoside-diphosphate reductase subunit M2 (RRM2), that are similarly depleted by XL44 treatment. XL44 induces hRpn13-dependent apoptosis and also restricts cell viability by a PCLAF-dependent mechanism. A KEN box, but not ubiquitination, is required for XL44-induced depletion of PCLAF. Here, we show that XL44 induces ubiquitin-dependent loss of hRpn13Pru and ubiquitin-independent loss of select KEN box containing proteins.

Rewiring of RNA methylation by the oncometabolite fumarate in renal cell carcinoma.

Metabolic reprogramming is a hallmark of cancer that facilitates changes in many adaptive biological processes. Mutations in the tricarboxylic acid cycle enzyme fumarate hydratase (FH) lead to fumarate accumulation and cause hereditary leiomyomatosis and renal cell cancer (HLRCC). HLRCC is a rare, inherited disease characterized by the development of non-cancerous smooth muscle tumors of the uterus and skin, and an increased risk of an aggressive form of kidney cancer. Fumarate has been shown to inhibit 2-oxoglutarate-dependent dioxygenases (2OGDDs) involved in the hydroxylation of HIF1α, as well as in DNA and histone demethylation. However, the link between fumarate accumulation and changes in RNA post-transcriptional modifications has not been defined. Here, we determine the consequences of fumarate accumulation on the activity of different members of the 2OGDD family targeting RNA modifications. By evaluating multiple RNA modifications in patient-derived HLRCC cell lines, we show that mutation of FH selectively affects the levels of N6-methyladenosine (m6A), while the levels of 5-formylcytosine (f5C) in mitochondrial tRNA are unaffected. This supports the hypothesis of a differential impact of fumarate accumulation on distinct RNA demethylases. The observation that metabolites modulate specific subsets of RNA-modifying enzymes offers new insights into the intersection between metabolism and the epitranscriptome.

Metabolism-focused CRISPR screen unveils mitochondrial pyruvate carrier 1 as a critical driver for PARP inhibitor resistance in lung cancer.

Homologous recombination (HR) and poly ADP-ribosylation are partially redundant pathways for the repair of DNA damage in normal and cancer cells. In cell lines that are deficient in HR, inhibition of poly (ADP-ribose) polymerase (poly (ADP-ribose) polymerase [PARP]1/2) is a proven target with several PARP inhibitors (PARPis) currently in clinical use. Resistance to PARPi often develops, usually involving genetic alterations in DNA repair signaling cascades, but also metabolic rewiring particularly in HR-proficient cells. We surmised that alterations in metabolic pathways by cancer drugs such as Olaparib might be involved in the development of resistance to drug therapy. To test this hypothesis, we conducted a metabolism-focused clustered regularly interspaced short palindromic repeats knockout screen to identify genes that undergo alterations during the treatment of tumor cells with PARPis. Of about 3000 genes in the screen, our data revealed that mitochondrial pyruvate carrier 1 (MPC1) is an essential factor in desensitizing nonsmall cell lung cancer (NSCLC) lung cancer lines to PARP inhibition. In contrast to NSCLC lung cancer cells, triple-negative breast cancer cells do not exhibit such desensitization following MPC1 loss and reprogram the tricarboxylic acid cycle and oxidative phosphorylation pathways to overcome PARPi treatment. Our findings unveil a previously unknown synergistic response between MPC1 loss and PARP inhibition in lung cancer cells.

Loss of NAT10 disrupts enhancer organization via p300 mislocalization and suppresses transcription of genes necessary for metastasis progression.

Acetylation of protein and RNA represent a critical event for development and cancer progression. NAT10 is the only known RNA acetylase that catalyzes the N4-actylcytidine (ac4C) modification of RNAs. Here, we show that the loss of NAT10 significantly decreases lung metastasis in allograft and genetically engineered mouse models of breast cancer. NAT10 interacts with a mechanosensitive, metastasis susceptibility protein complex at the nuclear pore. In addition to its canonical role in RNA acetylation, we find that NAT10 interacts with p300 at gene enhancers. NAT10 loss is associated with p300 mislocalization into heterochromatin regions. NAT10 depletion disrupts enhancer organization, leading to alteration of gene transcription necessary for metastatic progression, including reduced myeloid cell-recruiting chemokines that results in a less metastasis-prone tumor microenvironment. Our study uncovers a distinct role of NAT10 in enhancer organization of metastatic tumor cells and suggests its involvement in the tumor-immune crosstalk dictating metastatic outcomes.

Specialized replication of heterochromatin domains ensures self-templated chromatin assembly and epigenetic inheritance.

Heterochromatin, defined by histone H3 lysine 9 methylation (H3K9me), spreads across large domains and can be epigenetically inherited in a self-propagating manner. Heterochromatin propagation depends upon a read-write mechanism, where the Clr4/Suv39h methyltransferase binds to preexisting trimethylated H3K9 (H3K9me3) and further deposits H3K9me. How the parental methylated histone template is preserved during DNA replication is not well understood. Here, we demonstrate using Schizosaccharomyces pombe that heterochromatic regions are specialized replication domains demarcated by their surrounding boundary elements. DNA replication throughout these domains is distinguished by an abundance of replisome components and is coordinated by Swi6/HP1. Although mutations in the replicative helicase subunit Mcm2 that affect histone binding impede the maintenance of a heterochromatin domain at an artificially targeted ectopic site, they have only a modest impact on heterochromatin propagation via the read-write mechanism at an endogenous site. Instead, our findings suggest a crucial role for the replication factor Mcl1 in retaining parental histones and promoting heterochromatin propagation via a mechanism involving the histone chaperone FACT. Engagement of FACT with heterochromatin requires boundary elements, which position the heterochromatic domain at the nuclear peripheral subdomain enriched for heterochromatin factors. Our findings highlight the importance of replisome components and boundary elements in creating a specialized environment for the retention of parental methylated histones, which facilitates epigenetic inheritance of heterochromatin.

hRpn13 shapes the proteome and transcriptome through epigenetic factors HDAC8, PADI4, and transcription factor NF-κB p50.

The anti-cancer target hRpn13 is a proteasome substrate receptor. However, hRpn13-targeting molecules do not impair its interaction with proteasomes or ubiquitin, suggesting other critical cellular activities. We find that hRpn13 depletion causes correlated proteomic and transcriptomic changes, with pronounced effects in myeloma cells for cytoskeletal and immune response proteins and bone-marrow-specific arginine deiminase PADI4. Moreover, a PROTAC against hRpn13 co-depletes PADI4, histone deacetylase HDAC8, and DNA methyltransferase MGMT. PADI4 binds and citrullinates hRpn13 and proteasomes, and proteasomes from PADI4-inhibited myeloma cells exhibit reduced peptidase activity. When off proteasomes, hRpn13 can bind HDAC8, and this interaction inhibits HDAC8 activity. Further linking hRpn13 to transcription, its loss reduces nuclear factor κB (NF-κB) transcription factor p50, which proteasomes generate by cleaving its precursor protein. NF-κB inhibition depletes hRpn13 interactors PADI4 and HDAC8. Altogether, we find that hRpn13 acts dually in protein degradation and expression and that proteasome constituency and, in turn, regulation varies by cell type.