Transcriptional Network Analysis Research Articles

Viral and host network analysis of the human cytomegalovirus transcriptome in latency.

HCMV genes UL135 and UL138 play opposing roles regulating latency and reactivation in CD34+ human progenitor cells (HPCs). Using the THP-1 cell line model for latency and reactivation, we designed an RNA sequencing study to compare the transcriptional profile of HCMV infection in the presence and absence of these genes. The loss of UL138 results in elevated levels of viral gene expression and increased differentiation of cell populations that support HCMV gene expression and genome synthesis. The loss of UL135 results in diminished viral gene expression during an initial burst that occurs as latency is established and no expression of eleven viral genes from the ULb' region even following stimulation for differentiation and reactivation. Transcriptional network analysis revealed host transcription factors with potential to regulate the ULb' genes in coordination with pUL135. These results reveal roles for UL135 and UL138 in regulation of viral gene expression and potentially hematopoietic differentiation.

Integrated multi-omics investigation revealed the importance of phenylpropanoid metabolism in the defense response of Lilium regale Wilson to Fusarium wilt

Abstract Lilies (genus: Lilium) play a significant role in the global cut-flower industry, but they are highly susceptible to Fusarium wilt caused by Fusarium oxysporum. However, Lilium regale, a wild lily species, exhibits remarkable resistance to F. oxysporum. To investigate the quantitative resistance of L. regale to Fusarium wilt, a comprehensive multi-omics analysis was conducted. Upon inoculation with F. oxysporum, L. regale roots showed a significant accumulation of phenylpropane metabolites, including lignin precursors, flavonoids, and hydroxycinnamic acids. These findings were consistent with the up-regulated expression of phenylpropanoid biosynthesis-related genes encoding various enzymes, as revealed by transcriptomics and proteomics analyses. Furthermore, metabolomics and proteomics data demonstrated differential activation of monoterpenoid and isoquinoline alkaloid biosynthesis. Colorimetry and high-performance liquid chromatography analyses revealed significantly higher levels of total flavonoids, lignin, ferulic acid, phlorizin, and quercetin contents in L. regale scales compared to susceptible lily ‘Siberia’ scales during F. oxysporum infection. These phenylpropanes exhibited inhibitory effects on F. oxysporum growth and suppressed the expression of pathogenicity-related genes. Transcriptional regulatory network analysis suggested that ethylene-responsive transcription factors (ERFs) may positively regulate phenylpropanoid biosynthesis. Therefore, LrERF4 was cloned and transiently overexpressed in the Fusarium wilt-susceptible Oriental hybrid lily, ‘Siberia’. The overexpression of LrERF4 resulted in increased levels of total flavonoids, lignin, ferulic acid, phlorizin, and quercetin, while the silencing of LrERF4 in L. regale through RNAi had opposite effect. In conclusion, phenylpropanoid metabolism plays a crucial role in the defense response of L. regale against Fusarium wilt, with LrERF4 acting as a positive regulator of phenylpropane biosynthesis.

Identification of master regulator genes controlling pathogenic CD4+ T cell fate in inflammatory bowel disease through transcriptional network analysis

Inflammatory bowel diseases (IBD) are a group of chronic inflammatory conditions of the gastrointestinal tract associated with multiple pathogenic factors, including dysregulation of the immune response. Effector CD4+ T cells and regulatory CD4+ T cells (Treg) are central players in maintaining the balance between tolerance and inflammation. Interestingly, genetic modifications in these cells have been implicated in regulating the commitment of specific phenotypes and immune functions. However, the transcriptional program controlling the pathogenic behavior of T helper cells in IBD progression is still unknown. In this study, we aimed to find master transcription regulators controlling the pathogenic behavior of effector CD4+ T cells upon gut inflammation. To achieve this goal, we used an animal model of IBD induced by the transfer of naïve CD4+ T cells into recombination-activating gene 1 (Rag1) deficient mice, which are devoid of lymphocytes. As a control, a group of Rag1−/− mice received the transfer of the whole CD4+ T cells population, which includes both effector T cells and Treg. When gut inflammation progressed, we isolated CD4+ T cells from the colonic lamina propria and spleen tissue, and performed bulk RNA-seq. We identified differentially up- and down-regulated genes by comparing samples from both experimental groups. We found 532 differentially expressed genes (DEGs) in the colon and 30 DEGs in the spleen, mostly related to Th1 response, leukocyte migration, and response to cytokines in lamina propria T-cells. We integrated these data into Gene Regulatory Networks to identify Master Regulators, identifying four up-regulated master gene regulators (Lef1, Dnmt1, Mybl2, and Jup) and only one down-regulated master regulator (Foxo3). The altered expression of master regulators observed in the transcriptomic analysis was confirmed by qRT-PCR analysis and found an up-regulation of Lef1 and Mybl2, but without differences on Dnmt1, Jup, and Foxo3. These two master regulators have been involved in T cells function and cell cycle progression, respectively. We identified two master regulator genes associated with the pathogenic behavior of effector CD4+ T cells in an animal model of IBD. These findings provide two new potential molecular targets for treating IBD.

Abstract 6883: Investigation of tumor-immune interactions in pediatric sarcomas using single-cell multi-omics

Abstract Background: Pediatric sarcoma represents the second most common solid tumor in children. Despite advances in multimodal therapies, one-third of sarcoma patients do not survive the disease. Over the past few decades, the development of novel therapeutic strategies for sarcoma has been impeded by an insufficient understanding of the tumor microenvironment (TME). In particular, the potential of immunotherapy remains constrained by the gap of knowledge about the molecular and functional phenotypes of the immune cells in the sarcoma TME. Methods: To address this knowledge gap, we collected treatment-naïve tumor samples from 32 patients with five pediatric sarcoma subtypes. We employed a single-nuclei RNA+ATAC-Seq multi-omics assay to comprehensively characterize the transcriptomic and epigenomic heterogeneity of tumor and immune cells within the sarcoma TME. Results: Our analysis revealed that tumor-associated macrophages (TAMs) were the most abundant immune cell type in the sarcoma TME. We identified seven distinct TAM subpopulations, each exhibiting unique gene expression and chromatin accessibility patterns. Among these, TAM subsets with higher expression of SPP1 and C1QC genes displayed the strongest immunosuppressive potential. Transcriptional regulatory network (TRN) analysis identified transcription factors that were most likely to drive the phenotypes of SPP1+ and C1QC+ TAMs, including ELF2, JUND and RUNX1. Through cell-cell communication analysis, we discovered that SPP1+ TAMs interact with T cells via SPP1-CD44 signaling, while C1QC+ TAMs interact with NK cells through HLA-E-KLRC1 signaling, suggesting the potential involvement of these signaling pathways in TAM-mediated suppression of these cytotoxic lymphocytes. We prioritized ligand-receptor gene pairs, such as SEMA3D-NRP1 and LTBP3-ITGB5, which likely mediate the interaction between malignant sarcoma cells and SPP1+/C1QC+ TAMs based on their dominance across sarcoma subtypes and association with poorer clinical outcomes. On the other hand, our analysis revealed that sarcoma subtype-specific interactions between TAMs and malignant cells are common, highlighting the heterogeneity of signaling events across sarcoma subtypes. Conclusion and Significance: Our study presents a comprehensive single-cell multi-omics characterization of pediatric sarcomas and sheds light on its tumor microenvironment. By identifying unique markers for immunosuppressive TAM subpopulations and elucidating intercellular interactions promoting tumor development, this study proposes potential therapeutic targets for sarcoma treatment based on the role of TAMs in sarcoma TME. Citation Format: Zhan Zhang, Kyung Jin Ahn, Changya Chen, Rumeysa Biyik Sit, Hesham Mohei, Brian Lockhart, Malay Haldar, Kathrin Bernt, Vinodh Pillai, Kai Tan. Investigation of tumor-immune interactions in pediatric sarcomas using single-cell multi-omics [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 6883.

The relationship between extreme inter-individual variation in macrophage gene expression and genetic susceptibility to inflammatory bowel disease

The differentiation of resident intestinal macrophages from blood monocytes depends upon signals from the macrophage colony-stimulating factor receptor (CSF1R). Analysis of genome-wide association studies (GWAS) indicates that dysregulation of macrophage differentiation and response to microorganisms contributes to susceptibility to chronic inflammatory bowel disease (IBD). Here, we analyzed transcriptomic variation in monocyte-derived macrophages (MDM) from affected and unaffected sib pairs/trios from 22 IBD families and 6 healthy controls. Transcriptional network analysis of the data revealed no overall or inter-sib distinction between affected and unaffected individuals in basal gene expression or the temporal response to lipopolysaccharide (LPS). However, the basal or LPS-inducible expression of individual genes varied independently by as much as 100-fold between subjects. Extreme independent variation in the expression of pairs of HLA-associated transcripts (HLA-B/C, HLA-A/F and HLA-DRB1/DRB5) in macrophages was associated with HLA genotype. Correlation analysis indicated the downstream impacts of variation in the immediate early response to LPS. For example, variation in early expression of IL1B was significantly associated with local SNV genotype and with subsequent peak expression of target genes including IL23A, CXCL1, CXCL3, CXCL8 and NLRP3. Similarly, variation in early IFNB1 expression was correlated with subsequent expression of IFN target genes. Our results support the view that gene-specific dysregulation in macrophage adaptation to the intestinal milieu is associated with genetic susceptibility to IBD.

Systematic identification of wheat spike developmental regulators by integrated multi-omics, transcriptional network, GWAS, and genetic analyses

The spike architecture of wheat plays a crucial role in determining grain number, making it a key trait for optimization in wheat breeding programs. In this study, we used a multi-omic approach to analyze the transcriptome and epigenome profiles of the young spike at eight developmental stages, revealing coordinated changes in chromatin accessibility and H3K27me3 abundance during the flowering transition. We constructed a core transcriptional regulatory network (TRN) that drives wheat spike formation and experimentally validated a multi-layer regulatory module involving TaSPL15, TaAGLG1, and TaFUL2. By integrating the TRN with genome-wide association studies, we identified 227 transcription factors, including 42 with known functions and 185 with unknown functions. Further investigation of 61 novel transcription factors using multiple homozygous mutant lines revealed 36 transcription factors that regulate spike architecture or flowering time, such as TaMYC2-A1, TaMYB30-A1, and TaWRKY37-A1. Of particular interest, TaMYB30-A1, downstream of and repressed by WFZP, was found to regulate fertile spikelet number. Notably, the excellent haplotype of TaMYB30-A1, which contains a C allele at the WFZP binding site, was enriched during wheat breeding improvement in China, leading to improved agronomic traits. Finally, we constructed a free and open access Wheat Spike Multi-Omic Database (http://39.98.48.156:8800/#/). Our study identifies novel and high-confidence regulators and offers an effective strategy for dissecting the genetic basis of wheat spike development, with practical value for wheat breeding.

Alveolar epithelial progenitor cells require Nkx2-1 to maintain progenitor-specific epigenomic state during lung homeostasis and regeneration

Lung epithelial regeneration after acute injury requires coordination cellular coordination to pattern the morphologically complex alveolar gas exchange surface. During adult lung regeneration, Wnt-responsive alveolar epithelial progenitor (AEP) cells, a subset of alveolar type 2 (AT2) cells, proliferate and transition to alveolar type 1 (AT1) cells. Here, we report a refined primary murine alveolar organoid, which recapitulates critical aspects of in vivo regeneration. Paired scRNAseq and scATACseq followed by transcriptional regulatory network (TRN) analysis identified two AT1 transition states driven by distinct regulatory networks controlled in part by differential activity of Nkx2-1. Genetic ablation of Nkx2-1 in AEP-derived organoids was sufficient to cause transition to a proliferative stressed Krt8+ state, and AEP-specific deletion of Nkx2-1 in adult mice led to rapid loss of progenitor state and uncontrolled growth of Krt8+ cells. Together, these data implicate dynamic epigenetic maintenance via Nkx2-1 as central to the control of facultative progenitor activity in AEPs.

The ferroptosis landscape in acute myeloid leukemia.

Ferroptosis induction through the suppression of glutathione peroxidase 4 (GPX4) and apoptosis-inducing factor mitochondria-associated 2 (AIFM2) has proven to be an effective approach in eliminating chemotherapy-resistant cells of various types. However, a comprehensive understanding of the roles of GPX4 and AIFM2 in acute myeloid leukemia (AML) has not yet been achieved. Using cBioPortal, DepMap, GEPIA, Metascape, and ONCOMINE, we compared the transcriptional expression, survival data, gene mutation, methylation, and network analyses of GPX4- and AIFM2-associated signaling pathways in AML. The results revealed that high expression levels of GPX4 and AIFM2 are associated with an adverse prognosis for AML patients. Overexpression of AIFM2 correlated with elevated mutation frequencies in NPM1 and DNMT3A. GPX4 upregulation modulated the following pathways: GO:0045333, cellular respiration; R-HSA-5389840, mitochondrial translation elongation; GO:0009060, aerobic respiration; R-HSA-9609507, protein localization; and R-HSA-8953854, metabolism of RNA. On the other hand, the overexpression of AIFM2 influenced the following processes: GO:0048704, embryonic skeletal system morphogenesis; GO:0021546, rhombomere development; GO:0009954, proximal/distal pattern formation; and GO:0048732, gland development. This study identifies the high expression of GPX4 and AIFM2 as novel biomarkers predicting a poor prognosis for AML patients. Furthermore, ferroptosis induction may improve the stratified treatment of AML.

Research on radiotherapy related genes and prognostic target identification of rectal cancer based on multi-omics

BackgroundRadiosensitivity of rectal cancer is related to the radiotherapy efficacy and prognosis of patients with rectal cancer, and the genes and molecular mechanisms related to radiosensitivity of rectal cancer have not been clarified. We explored the radiosensitivity related genes of rectal cancer at a multi omics level.MethodsmRNA expression data and rectum adenocarcinoma (READ) data were obtained from the Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus Database (GEO) (GSE150082, GSE60331, GSE46862, GSE46862). Differentially expressed genes between radiotherapy sensitive group and radiotherapy insensitive group were screened. GO analysis and KEGG pathway analysis were performed for differentially expressed genes. Among the differentially expressed genes, five core genes associated with rectal cancer prognosis were selected using random survival forest analysis. For these five core genes, drug sensitivity analysis, immune cell infiltration analysis, TISIDB database immune gene correlation analysis, GSEA enrichment analysis, construction of Nomogram prediction model, transcriptional regulatory network analysis, and qRT-PCR validation was performed on human rectal adenocarcinoma tissue.ResultsWe found that 600 up-regulated genes and 553 down-regulated genes were significantly different between radiotherapy sensitive group and radiotherapy insensitive group in rectal cancer. Five key genes, TOP2A, MATR3, APOL6, JOSD1, and HOXC6, were finally screened by random survival forest analysis. These five key genes were associated with different immune cell infiltration, immune-related genes, and chemosensitivity. A comprehensive transcriptional regulatory network was constructed based on these five core genes. qRT-PCR revealed that MATR3 expression was different in rectal cancer tissues and adjacent non-cancerous tissues, while APOL6, HOXC6, JOSD1, and TOP2A expression was not different.ConclusionFive radiosensitivity-related genes related to the prognosis of rectal cancer: TOP2A, MATR3, APOL6, JOSD1, HOXC6, are involved in multiple processes such as immune cell infiltration, immune-related genes, chemosensitivity, signaling pathways and transcriptional regulatory networks and may be potential biomarkers for radiotherapy of rectal cancer.

Multi-omics network analysis on samples from sequential biopsies reveals vital role of proliferation arrest for Macrosteatosis related graft failure in rats after liver transplantation

To investigate the molecular impact of graft MaS on post-transplant prognosis, based on multi-omics integrative analysis. Rats were fed by methionine-choline deficient diet (MCD) for MaS grafts. Samples were collected from grafts by sequential biopsies. Transcriptomic and metabolomic profilings were assayed. Post-transplant MaS status showed a close association with graft failure. Differentially expressed genes (DEGs) for in-vivo MaS were mainly enriched on pathways of cell cycle and DNA replication. Post-transplant MaS caused arrests of graft regeneration via inhibiting the E2F1 centered network, which was confirmed by an in vitro experiment. Data from metabolomics assays found insufficient serine/creatine which is located on one‑carbon metabolism was responsible for MaS-related GF. Pre-transplant MaS caused severe fibrosis in long-term survivors. DEGs for grafts from long-term survivors with pre-transplant MaS were mainly enriched in pathways of ECM-receptor interaction and focal adhesion. Transcriptional regulatory network analysis confirmed SOX9 as a key transcription factor (TF) for MaS-related fibrosis. Metabolomic assays found elevation of aromatic amino acid (AAA) was a major feature of fibrosis in long-term survivors. Graft MaS in vivo increased post-transplant GF via negative regulations on graft regeneration. Pre-transplant MaS induced severe fibrosis in long-term survivors via activations on ECM-receptor interaction and AAA metabolism.

GOAT: Gene-level biomarker discovery from multi-Omics data using graph ATtention neural network for eosinophilic asthma subtype.

Asthma is a heterogeneous disease where various subtypes are established and molecular biomarkers of the subtypes are yet to be discovered. Recent availability of multi-omics data paved a way to discover molecular biomarkers for the subtypes. However, multi-omics biomarker discovery is challenging because of the complex interplay between different omics layers. We propose a deep attention model named Gene-level biomarker discovery from multi-Omics data using graph ATtention neural network (GOAT) for identifying molecular biomarkers for eosinophilic asthma subtypes with multi-omics data. GOAT identifies genes that discriminate subtypes using a graph neural network by modeling complex interactions among genes as the attention mechanism in the deep learning model. In experiments with multi-omics profiles of the COREA (Cohort for Reality and Evolution of Adult Asthma in Korea) asthma cohort of 300 patients, GOAT outperforms existing models and suggests interpretable biological mechanisms underlying asthma subtypes. Importantly, GOAT identified genes that are distinct only in terms of relationship with other genes through attention. To better understand the role of biomarkers, we further investigated two transcription factors, CTNNB1 and JUN, captured by GOAT. We were successful in showing the role of the transcription factors in eosinophilic asthma pathophysiology in a network propagation and transcriptional network analysis, which were not distinct in terms of gene expression level differences. Source code is available https://github.com/DabinJeong/Multi-omics_biomarker. The preprocessed data underlying this article is accessible in data folder of the github repository. Raw data are available in Multi-Omics Platform at http://203.252.206.90:5566/, and it can be accessible when requested.

Comparative transcriptional and co-expression network analysis of two upland cotton accessions with extreme phenotypic differences reveals molecular mechanisms of fiber development.

Upland cotton (Gossypium hirsutum) is the main source of natural fiber in the global textile industry, and thus its fiber quality and yield are important parameters. In this study, comparative transcriptomics was used to analyze differentially expressed genes (DEGs) due to its ability to effectively screen candidate genes during the developmental stages of cotton fiber. However, research using this method is limited, particularly on fiber development. The aim of this study was to uncover the molecular mechanisms underlying the whole period of fiber development and the differences in transcriptional levels. Comparative transcriptomes are used to analyze transcriptome data and to screen for differentially expressed genes. STEM and WGCNA were used to screen for key genes involved in fiber development. qRT-PCR was performed to verify gene expression of selected DEGs and hub genes. Two accessions of upland cotton with extreme phenotypic differences, namely EZ60 and ZR014121, were used to carry out RNA sequencing (RNA-seq) on fiber samples from different fiber development stages. The results identified 704, 376, 141, 269, 761, and 586 genes that were upregulated, and 1,052, 476, 355, 259, 702, and 847 genes that were downregulated at 0, 5, 10, 15, 20, and 25 days post anthesis, respectively. Similar expression patterns of DEGs were monitored using short time-series expression miner (STEM) analysis, and associated pathways of DEGs within profiles were investigated. In addition, weighted gene co-expression network analysis (WGCNA) identified five key modules in fiber development and screened 20 hub genes involved in the development of fibers. Through the annotation of the genes, it was found that the excessive expression of resistance-related genes in the early fiber development stages affects the fiber yield, whereas the sustained expression of cell elongation-related genes is critical for long fibers. This study provides new information that can be used to improve fibers in newly developed upland cotton genotypes.

Canagliflozin mediates tumor suppression alone and in combination with radiotherapy in non-small cell lung cancer (NSCLC) through inhibition of HIF-1α.

Non-small cell lung cancer (NSCLC) has a poor prognosis and effective therapeutic strategies are lacking. The diabetes drug canagliflozin inhibits NSCLC cell proliferation and the mammalian target of rapamycin (mTOR) pathway, which mediates cell growth and survival, but it is unclear whether this drug can enhance response rates when combined with cytotoxic therapy. Here, we evaluated the effects of canagliflozin on human NSCLC response to cytotoxic therapy in tissue cultures and xenografts. Ribonucleic acid sequencing (RNA-seq), real-time quantitative PCR (RT-qPCR), metabolic function, small interfering ribonucleic acid (siRNA) knockdown and protein expression assays were used in mechanistic analyses. We found that canagliflozin inhibited proliferation and clonogenic survival of NSCLC cells, and augmented the efficacy of radiotherapy to mediate these effects and inhibit NSCLC xenograft growth. Canagliflozin treatment alone moderately inhibited mitochondrial oxidative phosphorylation and exhibited greater anti-proliferative capacity than specific mitochondrial complex-I inhibitors. The treament downregulated genes mediating hypoxia-inducible factor (HIF)-1α stability, metabolism and survival, activated adenosine monophosphate-activated protein kinase (AMPK) and inhibited mTOR, a critical activator of HIF-1α signaling. HIF-1α knockdown and stabilization experiments suggested that canagliflozin mediates anti-proliferative effects, in part, through suppression of HIF-1α. Transcriptional regulatory network analysis pinpointed histone deacetylase 2 (HDAC2), a gene suppressed by canagliflozin, as a key mediator of canagliflozin's transcriptional reprogramming. HDAC2 knockdown eliminated HIF-1α levels and enhanced the anti-proliferative effects of canagliflozin. HDAC2-regulated genes suppressed by canagliflozin are associated with poor prognosis in several clinical NSCLC datasets. In addition, we include evidence that canagliflozin also improves NSCLC response to chemotherapy. In summary, canagliflozin may be a promising therapy to develop in combination with cytotoxic therapy in NSCLC.

Identification of Genes Related to Resistance to Ichthyophthirius multifiliis Based on Co-expression Network Analysis in Grass Carp.

The ciliate protozoan Ichthyophthirius multifiliis is an essential parasite causing white spot disease in grass carp, leading to significant economic losses. Understanding the molecular basis of grass carp's response to I. multifiliis has important scientific and environmental values. The transcriptional network analysis offers a valuable strategy to decipher the changes in gene expression in grass carp infected with I. multifiliis. Our goal was to screen the genes and pathways involved in resistance to I. multifiliis in grass carp. The different traits exhibited by grass carp infected with I. multifiliis may be caused by the differences in gene expression among grass carp individuals. Herein, to reveal those resistance-associated genes against I. multifiliis infection, we performed RNA sequencing using weighted gene co-expression network analysis (WGCNA). The biological function analysis and hub gene annotation for highly relevant modules revealed that different pathogen recognition and clearance responses resulted in different resistance to I. multifiliis infection. Furthermore, gene enrichment analysis revealed that I. multifiliis invasion in the disease-resistant group mainly activated immune pathways, including scavenger receptor activity and kappa B kinase/NF-kappa B signaling. By the annotation of the highly correlated module of the hub gene, we revealed that the apoptosis and ribosome biogenesis-related genes were enriched in the disease-resistant grass carp. The results of the dark grey module showed that several genes were mainly enriched in the two-component system (ko02020) and steroid biosynthesis (ko00100), suggesting that they are resistance-associated and energy metabolism-associated genes. In the disease resistance group, hub genes mainly included Nlrc3, fos, AAP8, HAP2, HAX, cho2, and zgc:113,036. This study revealed the gene network associated with disease resistance after I. multifiliis infection. The disease resistance-related pathways and central genes identified in this study are candidate references for breeders breeding disease-resistant. The results of this study may also provide some references for the development of drugs to antagonize I. multifiliis infection.

Pien-Tze-Huang prevents hepatocellular carcinoma by inducing ferroptosis via inhibiting SLC7A11-GSH-GPX4 axis

BackgroundMalignant transformation from hepatic fibrosis to carcinogenesis may be a therapeutic target for hepatocellular carcinoma (HCC). The aim of this study was to evaluate anti-cancer efficacy of Pien-Tze-Huang (PZH), and to investigate the underlying mechanisms by integrating transcriptional regulatory network analysis and experimental validation.MethodsA diethylnitrosamine (DEN)-induced HCC model in rats was established and used to evaluate the anti-cancer efficacy of PZH. After detecting a transcriptomic profiling, the “disease-related gene–drug effective target” interaction network was constructed, and the candidate targets of PZH against malignant transformation from hepatic fibrosis to HCC were identified and verified in vitro.ResultsPZH effectively alleviated the pathological changes of hepatic fibrosis and cirrhosis, and inhibited tumor formation and growth in DEN-induced HCC rats. Additionally, the administration of PZH reduced the levels of various hepatic function-related serological indicators significantly. Mechanically, a ferroptosis-related SLC7A11-GSH-GPX4 axis might be one of potential targets of PZH against malignant transformation from hepatic fibrosis to HCC. Especially, high SLC7A11 expression may be associated with poor prognosis of HCC patients. Experimentally, the administration of PZH markedly increased the trivalent iron and ferrous ion, suppressed the expression levels of SLC7A11 and GPX4 proteins, and reduced the GSH/GSSG ratio in the liver tissues of DEN-induced HCC rats.ConclusionsOur data offer an evidence that PZH may effectively improve the hepatic fibrosis microenvironment and prevent the occurrence of HCC through promoting ferroptosis in tumor cells via inhibiting the SLC7A11-GSH-GPX4 axis, implying that PZH may be a potential candidate drug for prevention and treatment of HCC at an early stage.

Single-cell transcriptional profiling of hearts during cardiac hypertrophy reveals the role of MAMs in cardiomyocyte subtype switching

Pathological cardiac hypertrophy is the main predecessor of heart failure. Its pathology is sophisticated, and its progression is associated with multiple cellular processes. To explore new therapeutic approaches, more precise examination of cardiomyocyte subtypes and involved biological processes is required in response to hypertrophic stimuli. Mitochondria and the endoplasmic reticulum (ER) are two crucial organelles associated with the progression of cardiac hypertrophy and are connected through junctions known as mitochondria-associated endoplasmic reticulum membranes (MAMs). Although MAM genes are altered in cardiac hypertrophy, the importance of MAMs in cardiac hypertrophy and the expression pattern of MAMs in certain cardiac cell types require a comprehensive analysis. In this study, we analyzed the temporal expression of MAM proteins in the process of cardiac hypertrophy and observed that MAM-related proteins preferentially accumulated in cardiomyocytes at the initial stage of cardiac hypertrophy and underwent a gradual decline, which was synchronized with the proportion of two cardiomyocyte subtypes (CM2 and CM3). Meanwhile, these subtypes went through a functional switch during cardiac hypertrophy. Trajectory analysis suggested that there was a differentiation trajectory of cardiomyocyte subtypes from high to low MAM protein expression. Distinct regulon modules across different cardiomyocyte cell types were revealed by transcriptional regulatory network analysis. Furthermore, scWGCNA revealed that MAM-related genes were clustered into a module that correlated with diabetic cardiomyopathy. Altogether, we identified cardiomyocyte subtype transformation and the potential critical transcription factors involved, which may serve as therapeutic targets in combating cardiac hypertrophy.

Bromodomain and extraterminal (BET) domain proteins regulate immune checkpoints through AMPK-dependent and -independent pathways in peripheral T cell subsets

Abstract Immune checkpoint blockade unleashes anti-tumor immunity of tumor-infiltrating CD8+ T cells. Yet, many cancer patients are refractory: unknown mechanisms resolve responders from non-responders. Analysis of transcriptional networks regulating immune checkpoint genes may provide insight. We used multi-parameter FACS to explore expression of several checkpoint inhibition molecules, focusing on functional roles of 5′-Adenosine Monophosphate-activated Protein Kinase (AMPK) and BET transcriptional co-regulators (BRD2, BRD3 and BRD4), both known to be key signaling pathways disturbed in obesity and cancer, resulting in exhaustion of effector CD8+ T cells and tumor progression. Here, human primary PBMCs from normal, healthy, adult donors were activated ex vivowith anti-CD3 and anti-CD28. We then measured expression of inhibitory receptors PD1, CTLA4, TIM3 and TIGIT on T cell subsets. Small molecule inhibitors of BET protein expression, including pan-BET inhibitor JQ1 and BRD4-selective PROTAC degrader MZ1, identified BET protein-regulated targets. We found that BRD4 regulates PD1, CTLA4 and TIM3 (but not TIGIT) in both CD4+ and CD8+ T cells. Treatment with Compound C, a modulator of AMPK, defined a unique PD1 and CTLA4 network. We confirmed that T cells with elevated expression of inhibitory receptors were impaired in cytokine secretion and proliferation. We conclude that separate signaling pathways likely converge on BRD2 and/or BRD3 and/or BRD4, acting independently as effectors at each promoter. Patterns of AMPK-regulated metabolism and BET gene expression are personalized and must be considered when evaluating immune checkpoint blockade in cancer therapy. Supported by grants from the Cancer Systems Biology Consortium of the National Cancer Institute (U01CA182898,U01CA243004; PI Denis), NCI Cancer Moonshot (R01CA222170; PI Denis) and Boston University's Immunology Training Program (T32AI007309-32; PI Kepler).

Single-cell transcriptomics uncovers the differentiation of a subset of murine esophageal progenitors into taste buds in vivo

Mouse esophagus is lined with a stratified epithelium, which is maintained by the constant renewal of unipotent progenitors. In this study, we profiled mouse esophagus by single-cell RNA sequencing and found taste buds specifically in the cervical segment of the esophagus. These taste buds have the same cellular composition as the ones from the tongue but express fewer taste receptor types. State-of-the-art transcriptional regulatory network analysis allowed the identification of specific transcription factors associated to the differentiation of immature progenitors into the three different taste bud cell types. Lineage tracing experiments revealed that esophageal taste buds arise from squamous bipotent progenitor, thus demonstrating that all esophageal progenitors are not unipotent. Our cell resolution characterization of cervical esophagus epithelium will enable a better understanding of esophageal progenitor potency and insights into the mechanisms involved in the development of taste buds.

TEAD4 is a master regulator of high-risk nasopharyngeal carcinoma

The molecular basis underlying nasopharyngeal carcinoma (NPC) remains unclear. Recent progress in transcriptional regulatory network analysis helps identify the master regulator (MR) proteins that transcriptionally define malignant tumor phenotypes. Here, we investigated transcription factor-target interactions and identified TEA domain transcription factor 4 (TEAD4) as an MR of high-risk NPC. Precisely, TEAD4 promoted NPC migration, invasion and cisplatin resistance, depending on its autopalmitoylation. Mechanistically, YTHDF2 (YTH domain family 2) recognized WTAP (Wilms tumor 1-associating protein)-mediated TEAD4 m6A methylation to facilitate its stability and led to aberrant up-regulation of TEAD4. Up-regulated TEAD4 further drove NPC progression by transcriptionally activating BZW2 (basic leucine zipper and W2 domains 2) to induce the oncogenic AKT pathway. Moreover, the transcriptional activity of TEAD4 was independent of its canonical coactivators YAP/TAZ. Clinically, TEAD4 serves as an independent predictor of unfavorable prognosis and cisplatin response in NPC. Our data revealed the crucial role of TEAD4 in driving tumor malignancy, thus, may provide therapeutic vulnerability in NPC.

Nuclear Receptor Subfamily 4A Signaling as a Key Disease Pathway of CD1c+ Dendritic Cell Dysregulation in Systemic Sclerosis.

This study was undertaken to identify key disease pathways driving conventional dendritic cell (cDC) alterations in systemic sclerosis (SSc). Transcriptomic profiling was performed on peripheral blood CD1c+ cDCs (cDC2s) isolated from 12 healthy donors and 48 patients with SSc, including all major disease subtypes. We performed differential expression analysis for the different SSc subtypes and healthy donors to uncover genes dysregulated in SSc. To identify biologically relevant pathways, we built a gene coexpression network using weighted gene correlation network analysis. We validated the role of key transcriptional regulators using chromatin immunoprecipitation (ChIP) sequencing and in vitro functional assays. We identified 17 modules of coexpressed genes in cDCs that correlated with SSc subtypes and key clinical traits, including autoantibodies, skin score, and occurrence of interstitial lung disease. A module of immunoregulatory genes was markedly down-regulated in patients with the diffuse SSc subtype characterized by severe fibrosis. Transcriptional regulatory network analysis performed on this module predicted nuclear receptor 4A (NR4A) subfamily genes (NR4A1, NR4A2, NR4A3) as the key transcriptional regulators of inflammation. Indeed, ChIP-sequencing analysis indicated that these NR4A members target numerous differentially expressed genes in SSc cDC2s. Inclusion of NR4A receptor agonists in culture-based experiments provided functional proof that dysregulation of NR4As affects cytokine production by cDC2s and modulates downstream T cell activation. NR4A1, NR4A2, and NR4A3 are important regulators of immunosuppressive and fibrosis-associated pathways in SSc cDCs. Thus, the NR4A family represents novel potential targets to restore cDC homeostasis in SSc.