Meta- and Orthogonal Integration of Influenza "OMICs" Data Defines a Role for UBR4 in Virus Budding

Deficiency of ATF2 retards senescence induced by replication stress and pamidronate in mouse jaw bone marrow stem cells

The aging process is associated with a loss of bone mass and an accumulation of senescent cells, which is under epigenetic control. Morphological and molecular analysis revealed a notable reduction in bone mass and alveolar crest height in aged mice, accompanied by increased levels of senescent mouse jaw bone marrow stem cells (mJBMSCs). To investigate whether specific transcription factors are involved, assay for transposase-accessible chromatin with sequencing (ATAC-seq) was performed on mJBMSCs isolated from 2-, 4-, 8-, and 20-month-old mice. In 20-month-old mJBMSCs, increased chromatin accessibility was observed alongside elevated expression of activating transcription factor 2 (ATF2) in both cells and alveolar bone. Silencing Atf2 in mJBMSCs failed to reverse physiological aging, but delayed replication stress and pamidronate (PAM) induced senescence. The analysis of ATAC-seq and RNA sequencing indicated that the differentially expressed genes upregulated by PAM but downregulated by ATF2 deficiency were related to some key biological processes, including negative regulation of cell proliferation, inflammatory response, adipogenesis, and cellular senescence. The dual-luciferase assay was conducted to demonstrate that ATF2 enhances Cdkn2a transcription by binding to its promoter region. Our findings suggest significant chromatin alterations in aged mJBMSCs, positioning ATF2 as a potential target for combating externally induced senescence.

Inhibition of SARS-CoV-2 growth in the lungs of mice by a peptide-conjugated morpholino oligomer targeting viral RNA

Further development of direct-acting antiviral agents against human SARS-CoV-2 infections remains a public health priority. Here, we report that an antisense peptide-conjugated morpholino oligomer (PPMO) named 5'END-2, targeting a highly conserved sequence in the 5' UTR of SARS-CoV-2 genomic RNA, potently suppressed SARS-CoV-2 growth in vitro and in vivo. In HeLa-ACE 2 cells, 5'END-2 produced IC50 values of between 40 nM and 1.15 μM in challenges using six genetically disparate strains of SARS-CoV-2, including JN.1. In vivo, using K18-hACE2 mice and the WA-1/2020 virus isolate, two doses of 5'END-2 at 10 mg/kg, administered intranasally on the day before and the day after infection, produced approximately 1.4 log10 virus titer reduction in lung tissue at 3 days post-infection. Under a similar dosing schedule, intratracheal administration of 1.0-2.0 mg/kg 5'END-2 produced over 3.5 log10 virus growth suppression in mouse lungs. Electrophoretic mobility shift assays characterized specific binding of 5'END-2 to its complementary target RNA. Furthermore, using reporter constructs containing SARS-CoV-2 5' UTR leader sequence, in an in-cell system, we observed that 5'END-2 could interfere with translation in a sequence-specific manner. The results demonstrate that direct pulmonary delivery of 5'END-2 PPMO is a promising antiviral strategy against SARS-CoV-2 infections and warrants further development.

Hyper-Interferon Sensitive Influenza Induces Adaptive Immune Responses and Overcomes Resistance to Anti-PD-1 in Murine Non-Small Cell Lung Cancer

Despite recent advances in immunotherapy with immune checkpoint inhibitors, many patients with non-small cell lung cancer (NSCLC) fail to respond or develop resistance after an initial response. In situ vaccination (ISV) with engineered viruses has emerged as a promising antigen-agnostic strategy that can both condition the tumor microenvironment and augment antitumor T-cell responses to overcome immune resistance. We engineered a live attenuated viral vaccine, hyper-IFN-sensitive (HIS) virus, by conducting a genome-wide functional screening and introducing eight IFN-sensitive mutations in the influenza genome to enhance host IFN response. Compared with wild-type influenza, HIS replication was attenuated in immunocompetent hosts, enhancing its potential as a safe option for cancer therapy. HIS ISV elicited robust yet transient type I IFN responses in murine NSCLCs, leading to an enrichment of polyfunctional effector Th1 CD4+ T cells and cytotoxic CD8+ T cells into the tumor. HIS ISV demonstrated enhanced antitumor efficacy compared with wild-type in multiple syngeneic murine models of NSCLC with distinct driver mutations and varying mutational burden. This efficacy was dependent on host type 1 IFN responses and T lymphocytes. HIS ISV overcame resistance to anti-PD-1 in LKB1-deficient murine NSCLC, resulting in improved overall survival and systemic tumor-specific immunity. These studies provide compelling evidence to support further clinical evaluation of HIS as an "off-the-shelf" ISV strategy for patients with NSCLC refractory to immune checkpoint inhibitors.

CaMKII-dependent non-canonical RIG-I pathway promotes influenza virus propagation in the acute-phase of infection

Ca2+/calmodulin-dependent protein kinase II (CaMKII) is one of hundreds of host-cell factors involved in the propagation of type A influenza virus (IAV), although its mechanism of action is unknown. Here, we identified CaMKII inhibitory peptide M3 by targeting its kinase domain using affinity-based screening of a tailored random peptide library. M3 inhibited IAV cytopathicity and propagation in cells by specifically inhibiting the acute-phase activation of retinoic acid-inducible gene I (RIG-I), which is uniquely regulated by CaMKII. Downstream of the RIG-I pathway activated TBK1 and then IRF3, which induced small but sufficient amounts of transcripts of the genes for IFN α/β to provide the capped 5'-ends that were used preferentially as primers to synthesize viral mRNAs by the cap-snatching mechanism. Importantly, knockout of RIG-I in cells almost completely inhibited the expression of IFN mRNAs and subsequent viral NP mRNA early in infection (up to 6 h after infection), which then protected cells from cytopathicity 24 h after infection. Thus, CaMKII-dependent acute-phase activation of RIG-I promoted IAV propagation, whereas the canonical RIG-I pathway stimulated antiviral activity by inducing large amounts of mRNA for IFNs and then for antiviral proteins later in infection. Co-administration of M3 with IAV infection rescued mice from the lethality and greatly reduced proinflammatory cytokine mRNA expression in the lung, indicating that M3 is highly effective against IAV in vivo. Thus, regulation of the CaMKII-dependent non-canonical RIG-I pathway may provide a novel host-factor-directed antiviral therapy.IMPORTANCEThe recent emergence of IAV strains resistant to commonly used therapeutic agents that target viral proteins has exacerbated the need for innovative strategies. Here, we originally identified CaMKII-inhibitory peptide M3, which efficiently inhibits IAV-lethality in vitro and in vivo. M3 specifically inhibited the acute-phase activation of RIG-I, which is a novel pathway to promote IAV propagation. Thus, this pathway acts in an opposite manner compared with the canonical RIG-I pathway, which plays essential roles in antiviral innate immune response later in infection. The CaMKII-dependent non-canonical RIG-I pathway can be a promising and novel drug target for the treatment of infections.

Neutrophil-centric analysis of gastric cancer: prognostic modeling and molecular insights

Gastric cancer remains a significant global health concern with poor prognosis. This study investigates the role of neutrophils in gastric cancer progression and their potential as prognostic indicators. Using multi-omics approaches, including Weighted Gene Co-expression Network Analysis (WGCNA), machine learning, and single-cell analysis, we identified neutrophil-associated gene signatures and developed a robust prognostic model. Our findings reveal distinct gastric cancer subtypes based on neutrophil-associated genes, with one subtype showing increased neutrophil infiltration and poorer prognosis. Single-cell analysis uncovered neutrophil-associated alterations in cell composition, gene expression profiles, and intercellular communication within the tumor microenvironment. Additionally, we explored the relationship between neutrophil-associated genes, microbiota composition, and alternative splicing events in gastric cancer. Furthermore, we identified QKI as a key regulator of alternative splicing and demonstrated its role in promoting malignant phenotypes and enhancing TGF-beta signaling and epithelial-mesenchymal transition in gastric cancer cells by wet experiment. Lastly, the role of QKI in the association with drug resistance and the identification of specific agents for treating QKI-associated drug resistance were also explored. This comprehensive study provides novel insights into the complex interplay between neutrophils, the tumor microenvironment, microbiota, alternative splicing and gastric cancer progression, offering potential new targets for therapeutic intervention.

Rapid and large-scale glycopeptide enrichment strategy based on chemical ligation

Protein glycosylation, the most universal post-translational modification, is thought to play a crucial role in regulating multiple essential cellular processes. However, the low abundance of glycoproteins and the heterogeneity of glycans complicate their comprehensive analysis. Here, we develop a rapid and large-scale glycopeptide enrichment strategy via bioorthogonal ligation and trypsin cleavage. The enrichment process is performed in one tube to minimize sample loss and time costs. This method combines convenience and practicality, identifying over 900 O-GlcNAc sites from a 500 μg sample. Surprisingly, it allows simultaneous identification of N-glycosites, O-GlcNAc sites, O-GalNAc sites and N-glycans via a two-step enzymatic release strategy. Combined with quantitative analysis, it reveals the distinct O-GlcNAcylation patterns in different compartments during oxidative stress. In summary, our study offers a convenient and robust tool for glycoproteome and glycome profiling, facilitating in-depth analysis to elucidate the biological functions of glycosylation.

Protein turnover regulation is critical for influenza A virus infection

The abundance of a protein is defined by its continuous synthesis and degradation, a process known as protein turnover. Here, we systematically profiled the turnover of proteins in influenza A virus (IAV)-infected cells using a pulse-chase stable isotope labeling by amino acids in cell culture (SILAC)-based approach combined with downstream statistical modeling. We identified 1,798 virus-affected proteins with turnover changes (tVAPs) out of 7,739 detected proteins (data available at pulsechase.innatelab.org). In particular, the affected proteins were involved in RNA transcription, splicing and nuclear transport, protein translation and stability, and energy metabolism. Many tVAPs appeared to be known IAV-interacting proteins that regulate virus propagation, such as KPNA6, PPP6C, and POLR2A. Notably, our analysis identified additional IAV host and restriction factors, such as the splicing factor GPKOW, that exhibit significant turnover rate changes while their total abundance is minimally affected. Overall, we show that protein turnover is a critical factor both for virus replication and antiviral defense.

Hederagenin regulates the migration and invasion of hepatocellular carcinoma cells through FOXO signaling pathway

OBJECTIVE

This study aimed to elucidate the effects of Hederagenin (HG) on hepatocellular carcinoma (HCC) and explore its potential molecular mechanisms.

MATERIALS AND METHODS

Virtual screening was employed to identify potential targets within core pathways of liver cancer and to analyze the possible mechanisms of HG. CCK-8 assays were used to assess the viability of HCC cells, while Hoechst 33342/PI staining was utilized to evaluate apoptosis. The migration and invasion abilities of HCC cells were examined using Transwell and scratch assays, and single-cell cloning ability was assessed via colony formation assays. Subsequent qRT-PCR was conducted to determine the mRNA expression levels of FOXO1 and FOXO6 following HG treatment. Western blot (WB) analysis was employed to measure the protein expression levels of IGF1R, FOXO1, FOXO6, MMP2, MMP9, and VEGFA, as well as the phosphorylation status of FOXO1 Ser249.

RESULTS

Virtual screening indicated that HG might exert antitumor effects through the FOXO signaling pathway. Experimental results demonstrated that HG induces apoptosis in a dose-dependent manner and inhibits the proliferation, migration, invasion, and single-cell cloning ability of HCC cells. After HG treatment, FOXO1 expression was upregulated, while the expression levels of IGF1R, phosphorylated FOXO1 Ser249, MMP2, MMP9, and VEGFA were downregulated.

CONCLUSION

In summary, our study is the first to demonstrate that HG regulates the phosphorylation of FOXO1, affecting the proliferation, migration, and invasion of HCC cells. The findings suggest that HG can inhibit the migration of HCC cells in vitro. The data indicate that HG-mediated targeting of the FOXO1/FOXO6 pathway holds promise as a novel therapeutic approach.

Inhibition of novel human-HPV hybrid ecDNA enhancers reduces oncogene expression and tumor growth in oropharyngeal cancer

Extrachromosomal circular DNA (ecDNA) have been found in most types of human cancers, and ecDNA incorporating viral genomes has recently been described, specifically in human papillomavirus (HPV)-mediated oropharyngeal cancer (OPC). However, the molecular mechanisms of human-viral hybrid ecDNA (hybrid ecDNA) for carcinogenesis remains elusive. We characterized the epigenetic status of hybrid ecDNA using HPVOPC cell lines and patient-derived tumor xenografts, identifying HPV oncogenes E6/E7 in hybrid ecDNA were flanked by novel somatic DNA enhancers and HPV L1 enhancers, with strong cis-interaction. Targeting of these enhancers by clustered regularly interspaced short palindromic repeats interference or hybrid ecDNA by bromodomain and extra-terminal inhibitor reduced E6/E7 expression, and significantly inhibited in vitro and/or in vivo growth only in ecDNA(+) models. HPV DNA in hybrid ecDNA structures are associated with novel somatic and HPV enhancers in hybrid ecDNA that drive HPV ongogene expression and carcinogenesis, and can be targeted with ecDNA disrupting therapeutics.

Stabilization of β-Catenin Directs HEB to Limit Thymic Selection

Activation of β-catenin in CD4+CD8+ double-positive (DP) thymocytes halts development before the thymic selection stage and predisposes to transformation. Leukemogenesis, but not the developmental block, depends on TCF-1, β-catenin's DNA-binding partner. In this study, we show that β-catenin activation directs the DNA-binding protein HEB to block DP thymocyte development. Conditional loss of HEB in DP thymocytes with stabilized β-catenin restores the frequencies of postselection TCRβhi/CCR7+ and TCRβhi/CD69+ DPs and their cell-cycle profile. This recovery is associated with significant reversal of β-catenin-induced expression changes, particularly those related to the CD69+ DP cell signature and to cell-cycle pathways. Stabilizing β-catenin in DP thymocytes directs HEB binding to ≈11,000 novel DNA sites throughout the genome. Novel HEB sites mark most CD69+ DP cell signature genes that change expression upon activation of β-catenin and then revert after loss of HEB. Moreover, many of the novel HEB sites occupy promoter regions of genes enriched in mitotic cell cycle pathways. HEB binding to those regions correlates with downregulation of the associated genes, and HEB inactivation restores expression to physiologic levels. These findings highlight a molecular interplay between HEB and β-catenin that can impair thymic development.

Single-cell EpiChem jointly measures drug-chromatin binding and multimodal epigenome

Studies of molecular and cellular functions of small-molecule inhibitors in cancer treatment, eliciting effects by targeting genome and epigenome associated proteins, requires measurement of drug-target engagement in single-cell resolution. Here we present EpiChem for in situ single-cell joint mapping of small molecules and multimodal epigenomic landscape. We demonstrate single-cell co-assays of three small molecules together with histone modifications, chromatin accessibility or target proteins in human colorectal cancer (CRC) organoids. Integrated multimodal analysis reveals diverse drug interactions in the context of chromatin states within heterogeneous CRC organoids. We further reveal drug genomic binding dynamics and adaptive epigenome across cell types after small-molecule drug treatment in CRC organoids. This method provides a unique tool to exploit the mechanisms of cell type-specific drug actions.

The M2 Protein of the Influenza A Virus Interacts with PEX19 to Facilitate Virus Replication by Disrupting the Function of Peroxisome

The peroxisomal biogenesis factor 19 (PEX19) is necessary for early peroxisomal biogenesis. PEX19 has been implicated in the replication of a variety of viruses, but the details pertaining to the mechanisms of how PEX19 engages in the life cycle of these viruses still need to be elucidated. Here, we demonstrated that the C terminus of PEX19 interacted with the cytoplasmic tail region of the M2 protein of the influenza A virus (IAV) and inhibited the viral growth titers. IAV infection or PEX19 knockdown triggered a reduction in the peroxisome pool and led to the accumulation of ROS and cell damage, thereby creating favorable conditions for IAV replication. Moreover, a reduction in the peroxisome pool led to the attenuation of early antiviral response mediated by peroxisome MAVS and downstream type III interferons. This study also showed that the interaction between IAV M2 and PEX19 affected the binding of PEX19 to the peroxisome-associated protein PEX14 and peroxisome membrane protein 24 (PMP24). Collectively, our data demonstrate that host factor PEX19 suppresses the replication of the IAV, and the IAV employs its M2 protein to mitigate the restricting role of PEX19.

Integrated genomic/epigenomic analysis stratifies subtypes of clear cell ovarian carcinoma, highlighting their cellular origin

The cellular origin of clear cell ovarian carcinoma (CCOC), a major histological subtype of ovarian carcinoma remains elusive. Here, we explored the candidate cellular origin and identify molecular subtypes using integrated genomic/epigenomic analysis. We performed whole exome-sequencing, microarray, and DNA methylation array in 78 CCOC samples according to the original diagnosis. The findings revealed that ARID1A and/or PIK3CA mutations were mutually exclusive with DNA repair related genes, including TP53, BRCA1, and ATM. Clustering of CCOC and other ovarian carcinomas (n = 270) with normal tissues from the fallopian tube, ovarian surface epithelium, endometrial epithelium, and pelvic peritoneum mesothelium (PPM) in a methylation array showed that major CCOC subtypes (with ARID1A and/or PIK3CA mutations) were associated with the PPM-lile cluster (n = 64). This cluster was sub-divided into three clusters: (1) mismatch repair (MMR) deficient with tumor mutational burden-high (n = 2), (2) alteration of ARID1A (n = 51), and (3) ARID1A wild-type (n = 11). The remaining samples (n = 14) were subdivided into (4) ovarian surface epithelium-like (n = 11) and (5) fallopian tube-like (considered as high-grade serous histotype; n = 3). Among these, subtypes (1-3) and others (4 and 5) were found to be associated with immunoreactive signatures and epithelial-mesenchymal transition, respectively. These results contribute to the stratification of CCOC into biological subtypes.

OmicScope unravels systems-level insights from quantitative proteomics data

Shotgun proteomics analysis presents multifaceted challenges, demanding diverse tool integration for insights. Addressing this complexity, OmicScope emerges as an innovative solution for quantitative proteomics data analysis. Engineered to handle various data formats, it performs data pre-processing - including joining replicates, normalization, data imputation - and conducts differential proteomics analysis for both static and longitudinal experimental designs. Empowered by Enrichr with over 224 databases, OmicScope performs Over Representation Analysis (ORA) and Gene Set Enrichment Analysis (GSEA). Additionally, its Nebula module facilitates meta-analysis from independent datasets, providing a systems biology approach for enriched insights. Complete with a data visualization toolkit and accessible as Python package and a web application, OmicScope democratizes proteomics analysis, offering an efficient and high-quality pipeline for researchers.

Tropism for ciliated cells is the dominant driver of influenza viral burst size in the human airway

Influenza viruses pose a significant burden on global human health. Influenza has a broad cellular tropism in the airway, but how infection of different epithelial cell types impacts replication kinetics and burden in the airways is not fully understood. Using primary human airway cultures, which recapitulate the diverse epithelial cell landscape of the human airways, we investigated the impact of cell type composition on virus tropism and replication kinetics. Cultures were highly diverse across multiple donors and 30 independent differentiation conditions and supported a range of influenza replication. Although many cell types were susceptible to influenza, ciliated and secretory cells were predominantly infected. Despite the strong tropism preference for secretory and ciliated cells, which consistently make up 75% or more of infected cells, only ciliated cells were associated with increased virus production. Surprisingly, infected secretory cells were associated with overall reduced virus output. The disparate response and contribution to influenza virus production could be due to different pro- and antiviral interferon-stimulated gene signatures between ciliated and secretory populations, which were interrogated with single-cell RNA sequencing. These data highlight the heterogeneous outcomes of influenza virus infections in the complex cellular environment of the human airway and the disparate impacts of infected cell identity on multiround burst size, even among preferentially infected cell types.

Possible association of G6PC2 and MUC6 induced by low‑dose‑rate irradiation in mouse intestine with inflammatory bowel disease

Although there are several types of radiation exposure, it is debated whether low‑dose‑rate (LDR) irradiation (IR) affects the body. Since the small intestine is a radiation‑sensitive organ, the present study aimed to evaluate how it changes when exposed to LDR IR and identify the genes sensitive to these doses. After undergoing LDR (6.0 mGy/h) γ radiation exposure, intestinal RNA from BALB/c mice was extracted 1 and 24 h later. Mouse whole genome microarrays were used to explore radiation‑induced transcriptional alterations. Reverse transcription‑quantitative (RT‑q) PCR was used to examine time‑ and dose‑dependent radiation responses. The histopathological status of the jejunum in the radiated mouse was not changed by 10 mGy of LDR IR; however, 23 genes were upregulated in response to LDR IR of the jejunum in mice after 1 and 24 h of exposure. Upregulated genes were selected to validate the results of the RNA sequencing analysis for RT‑qPCR detection and results showed that only Na+/K+ transporting subunit α4, glucose‑6‑phosphatase catalytic subunit 2 (G6PC2), mucin 6 (MUC6) and transient receptor potential cation channel subfamily V member 6 levels significantly increased after 24 h of LDR IR. Furthermore, G6PC2 and MUC6 were notable genes induced by LDR IR exposure according to protein expression via western blot analysis. The mRNA levels of G6PC2 and MUC6 were significantly elevated within 24 h under three conditions: i) Exposure to LDR IR, ii) repeated exposure to LDR IR and iii) exposure to LDR IR in the presence of inflammatory bowel disease. These results could contribute to an improved understanding of immediate radiation reactions and biomarker development to identify radiation‑susceptible individuals before histopathological changes become noticeable. However, further investigation into the specific mechanisms involving G6PC2 and MUC6 is required to accomplish this.

Nuclear receptor Nur77 and Yin-Yang 1 synergistically increase mitochondrial abundance and activity in macrophages

Mitochondrial biogenesis requires precise regulation of both mitochondrial-encoded and nuclear-encoded genes. Nuclear receptor Nur77 is known to regulate mitochondrial metabolism in macrophages and skeletal muscle. Here, we compared genome-wide Nur77 binding site and target gene expression in these two cell types, which revealed conserved regulation of mitochondrial genes and enrichment of motifs for the transcription factor Yin-Yang 1 (YY1). We show that Nur77 and YY1 interact, that YY1 increases Nur77 activity, and that their binding sites are co-enriched at mitochondrial ribosomal protein gene loci in macrophages. Nur77 and YY1 co-expression synergistically increases Mrpl1 expression as well as mitochondrial abundance and activity in macrophages but not skeletal muscle. As such, we identify a macrophage-specific Nur77-YY1 interaction that enhances mitochondrial metabolism.

A Carbon 21 Steroidal Glycoside with Pregnane Skeleton from Cynanchum atratum Bunge Promotes Megakaryocytic and Erythroid Differentiation in Erythroleukemia HEL Cells through Regulating Platelet-Derived Growth Factor Receptor Beta and JAK2/STAT3 Pathway

Erythroleukemia is a rare form of acute myeloid leukemia (AML). Its molecular pathogenesis remains vague, and this disease has no specific therapeutic treatments. Previously, our group isolated a series of Carbon 21 (C-21) steroidal glycosides with pregnane skeleton from the root of Cynanchum atratum Bunge. Among them, we found that a compound, named BW18, can induce S-phase cell cycle arrest and apoptosis via the mitogen-activated protein kinase (MAPK) pathway in human chronic myeloid leukemia K562 cells. However, its anti-tumor activity against erythroleukemia remains largely unknown. In this study, we aimed to investigate the anti-erythroleukemia activity of BW18 and the underlying molecular mechanisms. Our results demonstrated that BW18 exhibited a good anti-erythroleukemia activity in the human erythroleukemia cell line HEL and an in vivo xenograft mouse model. In addition, BW18 induced cell cycle arrest at the G2/M phase and promoted megakaryocytic and erythroid differentiation in HEL cells. Furthermore, RNA sequencing (RNA-seq) and rescue assay demonstrated that overexpression of platelet-derived growth factor receptor beta (PDGFRB) reversed BW18-induced megakaryocytic differentiation in HEL cells, but not erythroid differentiation. In addition, the network pharmacology analysis, the molecular docking and cellular thermal shift assay (CETSA) revealed that BW18 could inactivate Janus tyrosine kinase 2 (JAK2)/signal transducer and activator of transcription 3 (STAT3) pathway, which might mediate BW18-induced erythroid differentiation. Taken together, our findings elucidated a novel role of PDGFRB in regulating erythroleukemia differentiation and highlighted BW18 as an attractive lead compound for erythroleukemia treatment.

Systematic review and meta-analysis of genome-wide pooled CRISPR screens to identify host factors involved in influenza A virus infection

The host-virus interactome is increasingly recognized as an important research field to discover new therapeutic targets to treat influenza. Multiple pooled genome-wide CRISPR-Cas screens have been reported to identify new pro- and antiviral host factors of the influenza A virus. However, at present, a comprehensive summary of the results is lacking. We performed a systematic review of all reported CRISPR studies in this field in combination with a meta-analysis using the algorithm of meta-analysis by information content (MAIC). Two ranked gene lists were generated based on evidence in 15 proviral and 4 antiviral screens. Enriched pathways in the proviral MAIC results were compared to those of a prior array-based RNA interference (RNAi) meta-analysis. The top 50 proviral MAIC list contained genes whose role requires further elucidation, such as the endosomal ion channel TPCN1 and the kinase WEE1. Moreover, MAIC indicated that ALYREF, a component of the transcription export complex, has antiviral properties, whereas former knockdown experiments attributed a proviral role to this host factor. CRISPR-Cas-pooled screens displayed a bias toward early-replication events, whereas the prior RNAi meta-analysis covered early and late-stage events. RNAi screens led to the identification of a larger fraction of essential genes than CRISPR screens. In summary, the MAIC algorithm points toward the importance of several less well-known pathways in host-influenza virus interactions that merit further investigation. The results from this meta-analysis of CRISPR screens in influenza A virus infection may help guide future research efforts to develop host-directed anti-influenza drugs.

IMPORTANCE

Viruses rely on host factors for their replication, whereas the host cell has evolved virus restriction factors. These factors represent potential targets for host-oriented antiviral therapies. Multiple pooled genome-wide CRISPR-Cas screens have been reported to identify pro- and antiviral host factors in the context of influenza virus infection. We performed a comprehensive analysis of the outcome of these screens based on the publicly available gene lists, using the recently developed algorithm meta-analysis by information content (MAIC). MAIC allows the systematic integration of ranked and unranked gene lists into a final ranked gene list. This approach highlighted poorly characterized host factors and pathways with evidence from multiple screens, such as the vesicle docking and lipid metabolism pathways, which merit further exploration.

E4BP4 in macrophages induces an anti-inflammatory phenotype that ameliorates the severity of colitis

Macrophages are versatile cells of the innate immune system that work by altering their pro- or anti-inflammatory features. Their dysregulation leads to inflammatory disorders such as inflammatory bowel disease. We show that macrophage-specific upregulation of the clock output gene and transcription factor E4BP4 reduces the severity of colitis in mice. RNA-sequencing and single-cell analyses of macrophages revealed that increased expression of E4BP4 leads to an overall increase in expression of anti-inflammatory genes including Il4ra with a concomitant reduction in pro-inflammatory gene expression. In contrast, knockout of E4BP4 in macrophages leads to increased proinflammatory gene expression and decreased expression of anti-inflammatory genes. ChIP-seq and ATAC-seq analyses further identified Il4ra as a target of E4BP4, which drives anti-inflammatory polarization in macrophages. Together, these results reveal a critical role for E4BP4 in regulating macrophage inflammatory phenotypes and resolving inflammatory bowel diseases.

Defining the relationship of salivary gland malignancies to novel cell subpopulations in human salivary glands using single nucleus RNA-sequencing

Salivary glands have essential roles in maintaining oral health, mastication, taste and speech, by secreting saliva. Salivary glands are composed of several types of cells, and each cell type is predicted to be involved in the carcinogenesis of different types of cancers including adenoid cystic carcinoma (ACC), acinic cell carcinoma (AciCC), salivary duct carcinoma (SDC), myoepithelial carcinoma (MECA) and other histology. In our study, we performed single nucleus RNA-seq on three human salivary gland samples to clarify the gene expression profile of each complex cellular component of the salivary glands and related these expression patterns to expression found in salivary gland cancers (SGC) to infer cell of origin. By single nucleus RNA-seq, salivary gland cells were stratified into four clusters: acinar cells, ductal cells 1, ductal cells 2 and myoepithelial cells/stromal cells. The localization of each cell group was verified by IHC of each cluster marker gene, and one group of ductal cells was found to represent intercalated ductal cells labeled with HES1. Furthermore, in comparison with SGC RNA-seq data, acinar cell markers were upregulated in AciCC, but downregulated in ACC and ductal cell markers were upregulated in SDC but downregulated in MECA, suggesting that markers of origin are highly expressed in some SGC. Cell type expressions in specific SGC histology are similar to those found in normal salivary gland populations, indicating a potential etiologic relationship.

Xbp1 promotes odontoblastic differentiation through modulating mitochondrial homeostasis

Odontoblast differentiation depends on the orderly recruitment of transcriptional factors (TFs) in the transcriptional regulatory network. The depletion of crucial TFs disturbs dynamic alteration of the chromatin landscape and gene expression profile, leading to developmental defects. Our previous studies have revealed that the basic leucine zipper (bZIP) TF family is crucial in odontoblastic differentiation, but the function of bZIP TF family member XBP1 is still unknown. Here, we showed the stage-specific expression patterns of the spliced form Xbp1s during tooth development. Elevated Xbp1 expression and nuclear translocation of XBP1S in mesenchymal stem cells (MSCs) were induced by differentiation medium in vitro. Diminution of Xbp1 expression impaired the odontogenic differentiation potential of MSCs. The further integration of ATAC-seq and RNA-seq identified Hspa9 as a direct downstream target, an essential mitochondrial chaperonin gene that modulated mitochondrial homeostasis. The amelioration of mitochondrial dysfunction rescued the impaired odontogenic differentiation potential of MSCs caused by the diminution of Xbp1. Furthermore, the overexpression of Hspa9 rescued Xbp1-deficient defects in odontoblastic differentiation. Our study illustrates the crucial role of Xbp1 in odontoblastic differentiation via modulating mitochondrial homeostasis and brings evidence to the therapy of mitochondrial diseases caused by genetic defects.

Full-spectral genome analysis of natural killer/T cell lymphoma highlights impacts of genome instability in driving its progression

BACKGROUND

Natural killer/T cell lymphoma (NKTCL) is a clinically and genetically heterogeneous disease with poor prognosis. Genome sequencing and mutation characterization provides a powerful approach for patient stratification, treatment target discovery, and etiology identification. However, previous studies mostly concentrated on base-level mutations in primary NKTCL, whereas the large-scale genomic alterations in NKTCL and the mutational landscapes in relapsed/refractory NKTCL remain largely unexplored.

METHODS

Here, we assembled whole-genome sequencing and whole-exome sequencing data from 163 patients with primary or relapsed/refractory NKTCL and compared their somatic mutational landscapes at both nucleotide and structure levels.

RESULTS

Our study not only confirmed previously reported common NKTCL mutational targets like STAT3, TP53, and DDX3X but also unveiled several novel high-frequency mutational targets such as PRDM9, DST, and RBMX. In terms of the overall mutational landscape, we observed striking differences between primary and relapsed/refractory NKTCL patient groups, with the latter exhibits higher levels of tumor mutation burden, copy number variants (CNVs), and structural variants (SVs), indicating a strong signal of genomic instability. Complex structural rearrangements such as chromothripsis and focal amplification are also significantly enriched in relapsed/refractory NKTCL patients, exerting a substantial impact on prognosis. Accordingly, we devised a novel molecular subtyping system (i.e., C0-C4) with distinct prognosis by integrating potential driver mutations at both nucleotide and structural levels, which further provides an informative guidance for novel treatments that target these specific driver mutations and genome instability as a whole.

CONCLUSIONS

The striking differences underlying the mutational landscapes between the primary and relapsed/refractory NKTCL patients highlight the importance of genomic instability in driving the progression of NKTCL. Our newly proposed molecular subtyping system is valuable in assisting patient stratification and novel treatment design towards a better prognosis in the age of precision medicine.

Network-based identification of key proteins and repositioning of drugs for non-small cell lung cancer

BACKGROUND

NSCLC is a lethal cancer that is highly prevalent and accounts for 85% of cases of lung cancer. Conventional cancer treatments, such as chemotherapy and radiation, frequently exhibit limited efficacy and notable adverse reactions. Therefore, a drug repurposing method is proposed for effective NSCLC treatment.

AIMS

This study aims to evaluate candidate drugs that are effective for NSCLC at the clinical level using a systems biology and network analysis approach.

METHODS

Differentially expressed genes in transcriptomics data were identified using the systems biology and network analysis approaches. A network of gene co-expression was developed with the aim of detecting two modules of gene co-expression. Following that, the Drug-Gene Interaction Database was used to find possible drugs that target important genes within two gene co-expression modules linked to non-small cell lung cancer (NSCLC). The use of Cytoscape facilitated the creation of a drug-gene interaction network. Finally, gene set enrichment analysis was done to validate candidate drugs.

RESULTS

Unlike previous research on repositioning drugs for NSCLC, which uses a gene co-expression network, this project is the first to research both gene co-expression and co-occurrence networks. And the co-occurrence network also accounts for differentially expressed genes in cancer cells and their adjacent normal cells. For effective management of non-small cell lung cancer (NSCLC), drugs that show higher gene regulation and gene affinity within the drug-gene interaction network are thought to be important. According to the discourse, NSCLC genes have a lot of control over medicines like vincristine, fluorouracil, methotrexate, clotrimazole, etoposide, tamoxifen, sorafenib, doxorubicin, and pazopanib.

CONCLUSION

Hence, there is a possibility of repurposing these drugs for the treatment of non-small-cell lung cancer.

The function role of HIGD1A in nonalcoholic steatohepatitis from chronic hepatitis B

BACKGROUND

Accompanied by the growing prevalence of nonalcoholic fatty liver disease (NAFLD), the coexistence of chronic hepatitis B (CHB) and NAFLD has increased. In the context of CHB, there is limited understanding of the factors that influence the development of NASH.

METHODS

We enrolled CHB combined NAFLD patients who had liver biopsy and divided them to NASH vs. non-NASH groups. A whole transcriptome chip was used to examine the expression profiles of long noncoding RNAs (lncRNAs) and mRNA in biopsied liver tissues. The function analysis of HIGD1A were performed. We knocked down or overexpressed HIGD1A in HepG2.2.15 cells by transient transfection of siRNA-HIGD1A or pcDNA-HIGD1A. In vivo investigations were conducted using hepatitis B virus (HBV) transgenic mice.

RESULTS

In 65 patients with CHB and NAFLD, 28 were patients with NASH, and 37 were those without NASH. After screening 582 differentially expressed mRNAs, GO analysis revealed differentially expressed mRNAs acting on nicotinamide adenine dinucleotide phosphate (NADPH), which influenced redox enzyme activity. KEGG analysis also shown that they were involved in the NAFLD signaling pathway. The function analysis revealed that HIGD1A was associated with the mitochondrion. Then, both in vivo and in vitro CHB model, HIGD1A was significantly higher in the NASH group than in the non-NASH group. HIGD1A knockdown impaired mitochondrial transmembrane potential and induced cell apoptosis in HepG2.2.15 cells added oleic acid and palmitate. On the contrary, hepatic HIGD1A overexpression ameliorated free fatty acids-induced apoptosis and oxidative stress. Furthermore, HIGD1A reduced reactive oxygen species (ROS) level by increasing glutathione (GSH) expression, but Adenosine 5'-monophosphate (AMP)-activated protein kinase (AMPK)/Acetyl-CoA carboxylase (ACC) pathway was not involved.

CONCLUSION

Both in vivo and in vitro CHB model, an upward trend of HIGD1A was observed in the NASH-related inflammatory response. HIGDIA played a protective role in cells against oxidative stress. Our data suggested that HIGD1A may be a positive regulator of NASH within the CHB context.

Transcriptional Changes in Radiation-Induced Lung Injury: A Comparative Analysis of Two Radiation Doses for Preclinical Research

In a recent stereotactic body radiation therapy animal model, radiation pneumonitis and radiation pulmonary fibrosis were observed at around 2 and 6 weeks, respectively. However, the molecular signature of this model remains unclear. This study aimed to examine the molecular characteristics at these two stages using RNA-seq analysis. Transcriptomic profiling revealed distinct transcriptional patterns for each stage. Inflammatory response and immune cell activation were involved in both stages. Cell cycle processes and response to type II interferons were observed during the inflammation stage. Extracellular matrix organization and immunoglobulin production were noted during the fibrosis stage. To investigate the impact of a 10 Gy difference on fibrosis progression, doses of 45, 55, and 65 Gy were tested. A dose of 65 Gy was selected and compared with 75 Gy. The 65 Gy dose induced inflammation and fibrosis as well as the 75 Gy dose, but with reduced lung damage, fewer inflammatory cells, and decreased collagen deposition, particularly during the inflammation stage. Transcriptomic analysis revealed significant overlap, but differences were observed and clarified in Gene Ontology and KEGG pathway analysis, potentially influenced by changes in interferon-gamma-mediated lipid metabolism. This suggests the suitability of 65 Gy for future preclinical basic and pharmaceutical research connected with radiation-induced lung injury.

Bioinformatics and machine learning driven key genes screening for hepatocellular carcinoma

Liver cancer, a global menace, ranked as the sixth most prevalent and third deadliest cancer in 2020. The challenge of early diagnosis and treatment, especially for hepatocellular carcinoma (HCC), persists due to late-stage detections. Understanding HCC's complex pathogenesis is vital for advancing diagnostics and therapies. This study combines bioinformatics and machine learning, examining HCC comprehensively. Three datasets underwent meticulous scrutiny, employing various analytical tools such as Gene Ontology (GO) function and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis, protein interaction assessment, and survival analysis. These rigorous investigations uncovered twelve pivotal genes intricately linked with HCC's pathophysiological intricacies. Among them, CYP2C8, CYP2C9, EPHX2, and ESR1 were significantly positively correlated with overall patient survival, while AKR1B10 and NQO1 displayed a negative correlation. Moreover, the Adaboost prediction model yielded an 86.8 % accuracy, showcasing machine learning's potential in deciphering complex dataset patterns for clinically relevant predictions. These findings promise to contribute valuable insights into the elusive mechanisms driving liver cancer (HCC). They hold the potential to guide the development of more precise diagnostic methods and treatment strategies in the future. In the fight against this global health challenge, unraveling HCC's intricacies is of paramount importance.

Drug target prediction through deep learning functional representation of gene signatures

Many machine learning applications in bioinformatics currently rely on matching gene identities when analyzing input gene signatures and fail to take advantage of preexisting knowledge about gene functions. To further enable comparative analysis of OMICS datasets, including target deconvolution and mechanism of action studies, we develop an approach that represents gene signatures projected onto their biological functions, instead of their identities, similar to how the word2vec technique works in natural language processing. We develop the Functional Representation of Gene Signatures (FRoGS) approach by training a deep learning model and demonstrate that its application to the Broad Institute's L1000 datasets results in more effective compound-target predictions than models based on gene identities alone. By integrating additional pharmacological activity data sources, FRoGS significantly increases the number of high-quality compound-target predictions relative to existing approaches, many of which are supported by in silico and/or experimental evidence. These results underscore the general utility of FRoGS in machine learning-based bioinformatics applications. Prediction networks pre-equipped with the knowledge of gene functions may help uncover new relationships among gene signatures acquired by large-scale OMICs studies on compounds, cell types, disease models, and patient cohorts.

Remodeling of the endothelial cell transcriptional program via paracrine and DNA-binding activities of MPO

Myeloperoxidase (MPO) is an enzyme that functions in host defense. MPO is released into the vascular lumen by neutrophils during inflammation and may adhere and subsequently penetrate endothelial cells (ECs) coating vascular walls. We show that MPO enters the nucleus of ECs and binds chromatin independently of its enzymatic activity. MPO drives chromatin decondensation at its binding sites and enhances condensation at neighboring regions. It binds loci relevant for endothelial-to-mesenchymal transition (EndMT) and affects the migratory potential of ECs. Finally, MPO interacts with the RNA-binding factor ILF3 thereby affecting its relative abundance between cytoplasm and nucleus. This interaction leads to change in stability of ILF3-bound transcripts. MPO-knockout mice exhibit reduced number of ECs at scar sites following myocardial infarction, indicating reduced neovascularization. In summary, we describe a non-enzymatic role for MPO in coordinating EndMT and controlling the fate of endothelial cells through direct chromatin binding and association with co-factors.

Transcriptional programming of translation by BCL6 controls skeletal muscle proteostasis

Skeletal muscle is dynamically controlled by the balance of protein synthesis and degradation. Here we discover an unexpected function for the transcriptional repressor B cell lymphoma 6 (BCL6) in muscle proteostasis and strength in mice. Skeletal muscle-specific Bcl6 ablation in utero or in adult mice results in over 30% decreased muscle mass and force production due to reduced protein synthesis and increased autophagy, while it promotes a shift to a slower myosin heavy chain fibre profile. Ribosome profiling reveals reduced overall translation efficiency in Bcl6-ablated muscles. Mechanistically, tandem chromatin immunoprecipitation, transcriptomic and translational analyses identify direct BCL6 repression of eukaryotic translation initiation factor 4E-binding protein 1 (Eif4ebp1) and activation of insulin-like growth factor 1 (Igf1) and androgen receptor (Ar). Together, these results uncover a bifunctional role for BCL6 in the transcriptional and translational control of muscle proteostasis.

Time-of-day-dependent variation of the human liver transcriptome and metabolome is disrupted in MASLD

Background & Aims

Liver homeostasis is ensured in part by time-of-day-dependent processes, many of them being paced by the molecular circadian clock. Liver functions are compromised in metabolic dysfunction-associated steatotic liver disease (MASLD) and metabolic dysfunction-associated steatohepatitis (MASH), and clock disruption increases susceptibility to MASLD progression in rodent models. We therefore investigated whether the time-of-day-dependent transcriptome and metabolome are significantly altered in human steatotic and MASH livers.

Methods

Liver biopsies, collected within an 8 h-window from a carefully phenotyped cohort of 290 patients and histologically diagnosed to be either normal, steatotic or MASH hepatic tissues, were analyzed by RNA sequencing and unbiased metabolomic approaches. Time-of-day-dependent gene expression patterns and metabolomes were identified and compared between histologically normal, steatotic and MASH livers.

Results

Herein, we provide a first-of-its-kind report of a daytime-resolved human liver transcriptome-metabolome and associated alterations in MASLD. Transcriptomic analysis showed a robustness of core molecular clock components in steatotic and MASH livers. It also revealed stage-specific, time-of-day-dependent alterations of hundreds of transcripts involved in cell-to-cell communication, intracellular signaling and metabolism. Similarly, rhythmic amino acid and lipid metabolomes were affected in pathological livers. Both TNFα and PPARγ signaling were predicted as important contributors to altered rhythmicity.

Conclusion

MASLD progression to MASH perturbs time-of-day-dependent processes in human livers, while the differential expression of core molecular clock components is maintained.

Impact and implications

This work characterizes the rhythmic patterns of the transcriptome and metabolome in the human liver. Using a cohort of well-phenotyped patients (n = 290) for whom the time-of-day at biopsy collection was known, we show that time-of-day variations observed in histologically normal livers are gradually perturbed in liver steatosis and metabolic dysfunction-associated steatohepatitis. Importantly, these observations, albeit obtained across a restricted time window, provide further support for preclinical studies demonstrating alterations of rhythmic patterns in diseased livers. On a practical note, this study indicates the importance of considering time-of-day as a critical biological variable which may significantly affect data interpretation in animal and human studies of liver diseases.

Plasma cell-free DNA 5-hydroxymethylcytosine and whole-genome sequencing signatures for early detection of esophageal cancer

Esophageal cancer is a highly incidence and deadly disease with a poor prognosis, especially in developing countries. Owing to the lack of specific symptoms and early diagnostic biomarkers, most patients are diagnosed with advanced disease, leading to a 5-year survival rate of less than 15%. Early (n = 50) and middle-advanced (n = 50) esophageal squamous cell carcinoma (ESCC) patients, as well as 71 healthy individuals, underwent 5-hydroxymethylcytosine (5hmC) sequencing on their plasma cell-free DNA (cfDNA). A Northern Chinese cohort of cfDNA 5hmC dataset of 150 ESCC patients and 183 healthy individuals were downloaded for validation. A diagnostic model was developed using cfDNA 5hmC signatures and then improved by low-pass whole genome sequencing (WGS) features of cfDNA. Conserved cfDNA 5hmC modification motifs were observed in the two independent ESCC cohorts. The diagnostic model with 5hmC features achieved an AUC of 0.810 and 0.862 in the Southern and Northern cohorts, respectively, with sensitivities of 69.3-74.3% and specificities of 82.4-90.7%. The performance was well maintained in Stage I to Stage IV, with accuracy of 70-100%, but low in Stage 0, 33.3%. Low-pass WGS of cfDNA improved the AUC to 0.934 with a sensitivity of 82.4%, a specificity of 88.2%, and an accuracy of 84.3%, particularly significantly in Stage 0, with an accuracy up to 80%. 5hmC and WGS could efficiently differentiate very early ESCC from healthy individuals. These findings imply a non-invasive and convenient method for ESCC detection when clinical treatments are available and may eventually prolong survival.

GCLiPP: global crosslinking and protein purification method for constructing high-resolution occupancy maps for RNA binding proteins

GCLiPP is a global RNA interactome capture method that detects RNA-binding protein (RBP) occupancy transcriptome-wide. GCLiPP maps RBP-occupied sites at a higher resolution than phase separation-based techniques. GCLiPP sequence tags correspond with known RBP binding sites and are enriched for sites detected by RBP-specific crosslinking immunoprecipitation (CLIP) for abundant cytosolic RBPs. Comparison of human Jurkat T cells and mouse primary T cells uncovers shared peaks of GCLiPP signal across homologous regions of human and mouse 3' UTRs, including a conserved mRNA-destabilizing cis-regulatory element. GCLiPP signal overlapping with immune-related SNPs uncovers stabilizing cis-regulatory regions in CD5, STAT6, and IKZF1.

Integrated environmental, lifestyle, and epigenetic risk prediction of primary gastric neoplasia using the longitudinally monitored cohorts

BACKGROUND

DNA methylation accumulates in non-malignant gastric mucosa after exposure to pathogens. To elucidate how environmental, methylation, and lifestyle factors interplay to influence primary gastric neoplasia (GN) risk, we analyzed longitudinally monitored cohorts in Japan and Singapore.

METHODS

Asymptomatic subjects who underwent a gastric mucosal biopsy on the health check-up were enrolled. We analyzed the association between clinical factors and GN development using Cox hazard models. We further conducted comprehensive methylation analysis on selected tissues, including (i) mucosae from subjects developing GN later, (ii) mucosae from subjects not developing GN later, and (iii) GN tissues and surrounding mucosae. We also use the methylation data of mucosa collected in Singapore. The association between methylation and GN risk, as well as lifestyle and methylation, were analyzed.

FINDINGS

Among 4234 subjects, GN was developed in 77 subjects. GN incidence was correlated with age, drinking, smoking, and Helicobacter pylori (HP) status. Accumulation of methylation in biopsied gastric mucosae was predictive of higher future GN risk and shorter duration to GN incidence. Whereas methylation levels were associated with HP positivity, lifestyle, and morphological alterations, DNA methylation remained an independent GN risk factor through multivariable analyses. Pro-carcinogenic epigenetic alterations initiated by HP exposure were amplified by unfavorable but modifiable lifestyle choices. Adding DNA methylation to the model with clinical factors improved the predictive ability for the GN risk.

INTERPRETATION

The integration of environmental, lifestyle, and epigenetic information can provide increased resolution in the stratification of primary GN risk.

FUNDING

The funds are listed in Acknowledgements section.

Single-cell epigenomics and spatiotemporal transcriptomics reveal human cerebellar development

Human cerebellar development is orchestrated by molecular regulatory networks to achieve cytoarchitecture and coordinate motor and cognitive functions. Here, we combined single-cell transcriptomics, spatial transcriptomics and single cell chromatin accessibility states to systematically depict an integrative spatiotemporal landscape of human fetal cerebellar development. We revealed that combinations of transcription factors and cis-regulatory elements (CREs) play roles in governing progenitor differentiation and cell fate determination along trajectories in a hierarchical manner, providing a gene expression regulatory map of cell fate and spatial information for these cells. We also illustrated that granule cells located in different regions of the cerebellar cortex showed distinct molecular signatures regulated by different signals during development. Finally, we mapped single-nucleotide polymorphisms (SNPs) of disorders related to cerebellar dysfunction and discovered that several disorder-associated genes showed spatiotemporal and cell type-specific expression patterns only in humans, indicating the cellular basis and possible mechanisms of the pathogenesis of neuropsychiatric disorders.

Kaempferol prevents aseptic loosening via enhance the Wnt/β-catenin signaling pathway in vitro and in vivo

Kaempferol has demonstrated notable positive effects on the osteogenic differentiation of mesenchymal stem cells (MSC) and osteoblasts. A substantial body of research has emphasized the role of dislodged titanium particles in aseptic loosening following joint replacement surgery. This study predominantly investigates the suppressive influence of Kaempferol on osteolysis induced by titanium (Ti) alloy particles. In vitro investigations disclosed that Kaempferol effectively enhanced mineralization and alkaline phosphatase (ALP) activity in bone-marrow mesenchymal stem cells exposed to Ti particles. In addition, we conducted a comprehensive analysis of osteogenic differentiation microarray data_sets (GSE37676, GSE79814, and GSE114474) to identify differentially expressed genes. Significantly, Kaempferol upregulated the expression of critical osteogenic markers, including Runt-related transcription factor 2 (Runx2), osteocalcin (OCN), osterix/Sp-7, and β-catenin. In vivo experiments, including H&E staining and Immunohistochemistry, provided compelling evidence that Kaempferol exerted a robust inhibitory effect on periprosthetic osteolysis in mice, with particularly pronounced results at higher doses. Moreover, it elevated the expression levels of osteogenic factors and Wnt/β-catenin signaling components. These findings collectively indicate that Kaempferol mitigates the hindrance to osteogenesis posed by titanium particles by activating the Runx2 and Wnt/β-catenin signaling pathways. This research lays a solid foundation for the prospective utilization of Kaempferol in the management of aseptic loosening following arthroplasty, offering promising therapeutic potential.

A PPARγ/long noncoding RNA axis regulates adipose thermoneutral remodeling in mice

Interplay between energy-storing white adipose cells and thermogenic beige adipocytes contributes to obesity and insulin resistance. Irrespective of specialized niche, adipocytes require the activity of the nuclear receptor PPARγ for proper function. Exposure to cold or adrenergic signaling enriches thermogenic cells though multiple pathways that act synergistically with PPARγ; however, the molecular mechanisms by which PPARγ licenses white adipose tissue to preferentially adopt a thermogenic or white adipose fate in response to dietary cues or thermoneutral conditions are not fully elucidated. Here, we show that a PPARγ/long noncoding RNA (lncRNA) axis integrates canonical and noncanonical thermogenesis to restrain white adipose tissue heat dissipation during thermoneutrality and diet-induced obesity. Pharmacologic inhibition or genetic deletion of the lncRNA Lexis enhances uncoupling protein 1-dependent (UCP1-dependent) and -independent thermogenesis. Adipose-specific deletion of Lexis counteracted diet-induced obesity, improved insulin sensitivity, and enhanced energy expenditure. Single-nuclei transcriptomics revealed that Lexis regulates a distinct population of thermogenic adipocytes. We systematically map Lexis motif preferences and show that it regulates the thermogenic program through the activity of the metabolic GWAS gene and WNT modulator TCF7L2. Collectively, our studies uncover a new mode of crosstalk between PPARγ and WNT that preserves white adipose tissue plasticity.

Functional diversification and dynamics of CAR-T cells in patients with B-ALL

Understanding of cellular evolution and molecular programs of chimeric antigen receptor-engineered (CAR)-T cells post-infusion is pivotal for developing better treatment strategies. Here, we construct a longitudinal high-precision single-cell transcriptomic landscape of 7,578 CAR-T cells from 26 patients with B cell acute lymphoblastic leukemia (B-ALL) post-infusion. We molecularly identify eight CAR-T cell subtypes, including three cytotoxic subtypes with distinct kinetics and three dual-identity subtypes with non-T cell characteristics. Remarkably, long-term remission is coincident with the dominance of cytotoxic subtypes, while leukemia progression is correlated with the emergence of subtypes with B cell transcriptional profiles, which have dysfunctional features and might predict relapse. We further validate in vitro that the generation of B-featured CAR-T cells is induced by excessive tumor antigen stimulation or suppressed TCR signaling, while it is relieved by exogenous IL-12. Moreover, we define transcriptional hallmarks of CAR-T cell subtypes and reveal their molecular changes along computationally inferred cellular evolution in vivo. Collectively, these results decipher functional diversification and dynamics of peripheral CAR-T cells post-infusion.

An Inflammatory Checkpoint Generated by IL1RN Splicing Offers Therapeutic Opportunity for KRAS-Mutant Intrahepatic Cholangiocarcinoma

KRAS mutations are causally linked to protumor inflammation and are identified as driving factors in tumorigenesis. Here, using multiomics data gathered from a large set of patients, we showed that KRAS mutation was associated with a specific landscape of alternative mRNA splicing that connected to myeloid inflammation in intrahepatic cholangiocarcinoma (iCCA). Then, we identified a negative feedback mechanism in which the upregulation of interleukin 1 receptor antagonist (IL1RN)-201/203 due to alternative splicing confers vital anti-inflammatory effects in KRAS-mutant iCCA. In KRAS-mutant iCCA mice, both IL1RN-201/203 upregulation and anakinra treatment ignited a significant antitumor immune response by altering neutrophil recruitment and phenotypes. Furthermore, anakinra treatment synergistically enhanced anti-PD-1 therapy to activate intratumoral GZMB+ CD8+ T cells in KRAS-mutant iCCA mice. Clinically, we found that high IL1RN-201/203 levels in patients with KRAS-mutant iCCA were significantly associated with superior response to anti-PD-1 immunotherapy.

SIGNIFICANCE

This work describes a novel inflammatory checkpoint mediated by IL1RN alternative splicing variants that may serve as a promising basis to develop therapeutic options for KRAS-mutant iCCA and other cancers. This article is featured in Selected Articles from This Issue, p. 2109.

RANK ligand converts the NCoR/HDAC3 co-repressor to a PGC1β- and RNA-dependent co-activator of osteoclast gene expression

The nuclear receptor co-repressor (NCoR) complex mediates transcriptional repression dependent on histone deacetylation by histone deacetylase 3 (HDAC3) as a component of the complex. Unexpectedly, we found that signaling by the receptor activator of nuclear factor κB (RANK) converts the NCoR/HDAC3 co-repressor complex to a co-activator of AP-1 and NF-κB target genes that are required for mouse osteoclast differentiation. Accordingly, the dominant function of NCoR/HDAC3 complexes in response to RANK signaling is to activate, rather than repress, gene expression. Mechanistically, RANK signaling promotes RNA-dependent interaction of the transcriptional co-activator PGC1β with the NCoR/HDAC3 complex, resulting in the activation of PGC1β and inhibition of HDAC3 activity for acetylated histone H3. Non-coding RNAs Dancr and Rnu12, which are associated with altered human bone homeostasis, promote NCoR/HDAC3 complex assembly and are necessary for RANKL-induced osteoclast differentiation in vitro. These findings may be prototypic for signal-dependent functions of NCoR in other biological contexts.

Proteomic and genetic analyses of influenza A viruses identify pan-viral host targets

Influenza A Virus (IAV) is a recurring respiratory virus with limited availability of antiviral therapies. Understanding host proteins essential for IAV infection can identify targets for alternative host-directed therapies (HDTs). Using affinity purification-mass spectrometry and global phosphoproteomic and protein abundance analyses using three IAV strains (pH1N1, H3N2, H5N1) in three human cell types (A549, NHBE, THP-1), we map 332 IAV-human protein-protein interactions and identify 13 IAV-modulated kinases. Whole exome sequencing of patients who experienced severe influenza reveals several genes, including scaffold protein AHNAK, with predicted loss-of-function variants that are also identified in our proteomic analyses. Of our identified host factors, 54 significantly alter IAV infection upon siRNA knockdown, and two factors, AHNAK and coatomer subunit COPB1, are also essential for productive infection by SARS-CoV-2. Finally, 16 compounds targeting our identified host factors suppress IAV replication, with two targeting CDK2 and FLT3 showing pan-antiviral activity across influenza and coronavirus families. This study provides a comprehensive network model of IAV infection in human cells, identifying functional host targets for pan-viral HDT.

The Rho GTPase exchange factor Vav2 promotes extensive age-dependent rewiring of the hair follicle stem cell transcriptome

Both the number and regenerative activity of hair follicle stem cells (HFSCs) are regulated by Vav2, a GDP/GTP exchange factor involved in the catalytic stimulation of the GTPases Rac1 and RhoA. However, whether Vav2 signaling changes in HFSCs over the mouse lifespan is not yet known. Using a mouse knock-in mouse model, we now show that the expression of a catalytically active version of Vav2 (Vav2Onc) promotes an extensive rewiring of the overall transcriptome of HFSCs, the generation of new transcription factor hubs, and the synchronization of many transcriptional programs associated with specific HFSC states and well-defined signaling pathways. Interestingly, this transcriptome rewiring is not fixed in time, as it involves the induction of 15 gene expression waves with diverse distribution patterns during the life of the animals. These expression waves are consistent with the promotion by Vav2Onc of several functional HFSC states that differ from those normally observed in wild-type HFSCs. These results further underscore the role of Vav2 in the regulation of the functional state of HFSCs. They also indicate that, unlike other Vav2-dependent biological processes, the signaling output of this exchange factor is highly contingent on age-dependent intrinsic and/or extrinsic HFSC factors that shape the final biological readouts triggered in this cell type.

Potential selection for lipid kinase activity and spermatogenesis in Henan native pig breeds and growth shaping by introgression of European genes

BACKGROUND

China has one third of the worldwide indigenous pig breeds. The Henan province is one of the earliest pig domestication centers of China (about 8000 years ago). However, the precise genetic characteristics of the Henan local pig breeds are still obscure. To understand the origin and the effects of selection on these breeds, we performed various analyses on lineage composition, genetic structure, and detection of selection sweeps and introgression in three of these breeds (Queshan, Nanyang and Huainan) using genotyping data on 125 Queshan, 75 Nanyang, 16 Huainan pigs and 878 individuals from 43 Eurasian pig breeds.

RESULTS

We found no clear evidence of ancestral domestic pig DNA lineage in the Henan local breeds, which have an extremely complicated genetic background. Not only do they share genes with some northern Chinese pig breeds, such as Erhualian, Hetaodaer, and Laiwu, but they also have a high admixture of genes from foreign pig breeds (33-40%). Two striking selection sweeps in small regions of chromosomes 2 and 14 common to the Queshan and Nanyang breeds were identified. The most significant enrichment was for lipid kinase activity (GO:0043550) with the genes FII, AMBRA1, and PIK3IP1. Another interesting 636.35-kb region on chromosome 14 contained a cluster of spermatogenesis genes (OSBP2, GAL3ST1, PLA2G3, LIMK2, and PATZ1), a bisexual sterility gene MORC2, and a fat deposition gene SELENOM. Reproduction and growth genes LRP4, FII, and ARHGAP1 were present in a 238.05-kb region on SSC2 under selection. We also identified five loci associated with body length (P = 0.004) on chromosomes 1 and 12 that were introgressed from foreign pig breeds into the Henan breeds. In addition, the Chinese indigenous pig breeds fell into four main types instead of the previously reported six, among which the Eastern type could be divided into two subgroups.

CONCLUSIONS

Admixture of North China, East China and foreign pigs contributed to high genetic diversity of Henan local pigs. Ontology terms associated with lipid kinase activity and spermatogenesis and growth shaping by introgression of European genes in Henan pigs were identified through selective sweep analyses.

Exploring the dynamics and influencing factors of CD4 T cell activation using single-cell RNA-seq

T cell activation is a key event in adaptive immunity. However, the dynamics and influencing factors of T cell activation remain unclear. Here, we analyzed CD4 T cells that were stimulated with anti-CD3/CD28 under several conditions to explore the factors affecting T cell activation. We found a stimulated T subset (HSPhi T) highly expressing heat shock proteins, which was derived from stimulated naive T. We identified and characterized inert T, a stimulated T cell subset in transitional state from resting T to activated T. Interestingly, resting CXCR4low T responded to stimulation more efficiently than resting CXCR4hi T. Furthermore, stimulation of CD4 T in the presence of CD8 T resulted in more effector T and more homogeneous expressions of CD25, supporting that presence of CD8 T reduces the extreme response of T cells, which can be explained by regulation of CD4 T activation through CD8 T-initiated cytokine signaling and FAS/FASLG signaling.

Dynamic Glycoprotein Hyposialylation Promotes Chemotherapy Evasion and Metastatic Seeding of Quiescent Circulating Tumor Cell Clusters in Breast Cancer

Most circulating tumor cells (CTC) are detected as single cells, whereas a small proportion of CTCs in multicellular clusters with stemness properties possess 20- to 100-times higher metastatic propensity than the single cells. Here we report that CTC dynamics in both singles and clusters in response to therapies predict overall survival for breast cancer. Chemotherapy-evasive CTC clusters are relatively quiescent with a specific loss of ST6GAL1-catalyzed α2,6-sialylation in glycoproteins. Dynamic hyposialylation in CTCs or deficiency of ST6GAL1 promotes cluster formation for metastatic seeding and enables cellular quiescence to evade paclitaxel treatment in breast cancer. Glycoproteomic analysis reveals newly identified protein substrates of ST6GAL1, such as adhesion or stemness markers PODXL, ICAM1, ECE1, ALCAM1, CD97, and CD44, contributing to CTC clustering (aggregation) and metastatic seeding. As a proof of concept, neutralizing antibodies against one newly identified contributor, PODXL, inhibit CTC cluster formation and lung metastasis associated with paclitaxel treatment for triple-negative breast cancer.

SIGNIFICANCE

This study discovers that dynamic loss of terminal sialylation in glycoproteins of CTC clusters contributes to the fate of cellular dormancy, advantageous evasion to chemotherapy, and enhanced metastatic seeding. It identifies PODXL as a glycoprotein substrate of ST6GAL1 and a candidate target to counter chemoevasion-associated metastasis of quiescent tumor cells. This article is featured in Selected Articles from This Issue, p. 1949.

NETosis promotes chronic inflammation and fibrosis in systemic lupus erythematosus and COVID-19

Pulmonary fibrosis, a serious complication of systemic lupus erythematosus (SLE) and coronavirus disease 2019 (COVID-19), leads to irreversible lung damage. However, the underlying mechanism of this condition remains unclear. In this study, we revealed the landscape of transcriptional changes in lung biopsies from individuals with SLE, COVID-19-induced pulmonary fibrosis, and idiopathic pulmonary fibrosis (IPF) using histopathology and RNA sequencing, respectively. Despite the diverse etiologies of these diseases, lung expression of matrix metalloproteinase genes in these diseases showed similar patterns. Particularly, the differentially expressed genes were significantly enriched in the pathway of neutrophil extracellular trap formation, showing similar enrichment signature between SLE and COVID-19. The abundance of Neutrophil extracellular traps (NETs) was much higher in the lungs of individuals with SLE and COVID-19 compared to those with IPF. In-depth transcriptome analyses revealed that NETs formation pathway promotes epithelial-mesenchymal transition (EMT). Furthermore, stimulation with NETs significantly up-regulated α-SMA, Twist, Snail protein expression, while decreasing the expression of E-cadherin protein in vitro. This indicates that NETosis promotes EMT in lung epithelial cells. Given drugs that are efficacious in degrading damaged NETs or inhibiting NETs production, we identified a few drug targets that were aberrantly expressed in both SLE and COVID-19. Among these targets, the JAK2 inhibitor Tofacitinib could effectively disrupted the process of NETs and reversed NET-induced EMT in lung epithelial cells. These findings support that the NETs/EMT axis, activated by SLE and COVID-19, contributes to the progression of pulmonary fibrosis. Our study also highlights that JAK2 as a potential target for the treatment of fibrosis in these diseases.

Picolinic acid is a broad-spectrum inhibitor of enveloped virus entry that restricts SARS-CoV-2 and influenza A virus in vivo

The COVID-19 pandemic highlights an urgent need for effective antivirals. Targeting host processes co-opted by viruses is an attractive antiviral strategy with a high resistance barrier. Picolinic acid (PA) is a tryptophan metabolite endogenously produced in mammals. Here, we report the broad-spectrum antiviral activity of PA against enveloped viruses, including severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), influenza A virus (IAV), flaviviruses, herpes simplex virus, and parainfluenza virus. Mechanistic studies reveal that PA inhibits enveloped virus entry by compromising viral membrane integrity, inhibiting virus-cellular membrane fusion, and interfering with cellular endocytosis. More importantly, in pre-clinical animal models, PA exhibits promising antiviral efficacy against SARS-CoV-2 and IAV. Overall, our data establish PA as a broad-spectrum antiviral with promising pre-clinical efficacy against pandemic viruses SARS-CoV-2 and IAV.

Selective activation of PPARα maintains thermogenic capacity of beige adipocytes

Beige adipocytes are inducible thermogenic adipocytes used for anti-obesity treatment. Beige adipocytes rapidly lose their thermogenic capacity once external cues are removed. However, long-term administration of stimulants, such as PPARγ and β-adrenergic receptor agonists, is unsuitable due to various side effects. Here, we reported that PPARα pharmacological activation was the preferred target for maintaining induced beige adipocytes. Pemafibrate used in clinical practice for dyslipidemia was developed as a selective PPARα modulator (SPPARMα). Pemafibrate administration regulated the thermogenic capacity of induced beige adipocytes, repressed body weight gain, and ameliorated impaired glucose tolerance in diet-induced obese mouse models. The transcriptome analysis revealed that the E-twenty-six transcription factor ELK1 acted as a cofactor of PPARα. ELK1 was mobilized to the Ucp1 transcription regulatory region with PPARα and modulated its expression by pemafibrate. These results suggest that selective activation of PPARα by pemafibrate is advantageous to maintain the function of beige adipocytes.

Proteomics and phosphoproteomics profiling in glutamatergic neurons and microglia in an iPSC model of Jansen de Vries Syndrome

Background

Jansen de Vries Syndrome (JdVS) is a rare neurodevelopmental disorder (NDD) caused by gain-of-function (GOF) truncating mutations in PPM1D exons 5 or 6. PPM1D is a serine/threonine phosphatase that plays an important role in the DNA damage response (DDR) by negatively regulating TP53 (P53). JdVS-associated mutations lead to the formation of a truncated PPM1D protein that retains catalytic activity and has a GOF effect because of reduced degradation. Somatic PPM1D exons 5 and 6 truncating mutations are well-established factors in a number of cancers, due to excessive dephosphorylation and reduced function of P53 and other substrates involved in DDR. Children with JdVS have a variety of neurodevelopmental, psychiatric, and physical problems. In addition, a small fraction has acute neuropsychiatric decompensation apparently triggered by infection or severe non-infectious environmental stress factors.

Methods

To understand the molecular basis of JdVS, we developed an induced pluripotent stem cell (iPSC) model system. iPSCs heterozygous for the truncating variant (PPM1D+/tr), were made from a patient, and control lines engineered using CRISPR-Cas9 gene editing. Proteomics and phosphoprotemics analyses were carried out on iPSC-derived glutamatergic neurons and microglia from three control and three PPM1D+/tr iPSC lines. We also analyzed the effect of the TLR4 agonist, lipopolysaccharide, to understand how activation of the innate immune system in microglia could account for acute behavioral decompensation.

Results

One of the major findings was the downregulation of POGZ in unstimulated microglia. Since loss-of-function variants in the POGZ gene are well-known causes of autism spectrum disorder, the decrease in PPM1D+/tr microglia suggests this plays a role in the neurodevelopmental aspects of JdVS. In addition, neurons, baseline, and LPS-stimulated microglia show marked alterations in the expression of several E3 ubiquitin ligases, most notably UBR4, and regulators of innate immunity, chromatin structure, ErbB signaling, and splicing. In addition, pathway analysis points to overlap with neurodegenerative disorders.

Limitations

Owing to the cost and labor-intensive nature of iPSC research, the sample size was small.

Conclusions

Our findings provide insight into the molecular basis of JdVS and can be extrapolated to understand neuropsychiatric decompensation that occurs in subgroups of patients with ASD and other NDDs.

The Polyamine-Hypusine Circuit Controls an Oncogenic Translational Program Essential for Malignant Conversion in MYC-Driven Lymphoma

The MYC oncoprotein is activated in a broad spectrum of human malignancies and transcriptionally reprograms the genome to drive cancer cell growth. Given this, it is unclear if targeting a single effector of MYC will have therapeutic benefit. MYC activates the polyamine-hypusine circuit, which posttranslationally modifies the eukaryotic translation factor eIF5A. The roles of this circuit in cancer are unclear. Here we report essential intrinsic roles for hypusinated eIF5A in the development and maintenance of MYC-driven lymphoma, where the loss of eIF5A hypusination abolishes malignant transformation of MYC-overexpressing B cells. Mechanistically, integrating RNA sequencing, ribosome sequencing, and proteomic analyses revealed that efficient translation of select targets is dependent upon eIF5A hypusination, including regulators of G1-S phase cell-cycle progression and DNA replication. This circuit thus controls MYC's proliferative response, and it is also activated across multiple malignancies. These findings suggest the hypusine circuit as a therapeutic target for several human tumor types.

SIGNIFICANCE

Elevated EIF5A and the polyamine-hypusine circuit are manifest in many malignancies, including MYC-driven tumors, and eIF5A hypusination is necessary for MYC proliferative signaling. Not-ably, this circuit controls an oncogenic translational program essential for the development and maintenance of MYC-driven lymphoma, supporting this axis as a target for cancer prevention and treatment. See related commentary by Wilson and Klein, p. 248. This article is highlighted in the In This Issue feature, p. 247.