Set Of Target Genes Research Articles

Profile of miRNAs induced during sheep fat tail development and roles of four key miRNAs in proliferation and differentiation of sheep preadipocytes

BackgroundThe fat tail of sheep is an adaptive trait that facilitates their adaptation to harsh natural environments. MicroRNAs (miRNAs) have been demonstrated to play crucial roles in the regulation of tail fat deposition.MethodsIn this study, miRNA-Seq was employed to investigate the expression profiles of miRNAs during different developmental stages of sheep fat tails and elucidate the functions of differentially expressed miRNAs (DE miRNAs).ResultsA total of 350 DE miRNAs were identified, among which 191, 60, 26, and 21 were significantly upregulated in tail fat tissues of fetal, lamb, hogget Altay sheep, and adult Xinjiang fine wool (XFW) sheep but downregulated in other stages. Furthermore, we predicted a set of candidate target genes (4,476) for the top 20 DE miRNAs. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that they involve in several adipogenesis-related pathways. Subsequent investigations indicated that four DE miRNAs, miR-433-3p, miR-485-3p, miR-409-3p, and miR-495-3p, could suppress the expression of peroxisome proliferator-activated receptor gamma (PPARγ) and phosphoinositide-3-kinase regulatory subunit 3 (PIK3R3) and regulate the preadipocyte development in sheep. Meanwhile, the lipid metabolism-related genes, fatty acid-binding protein (FABP3), perilipin 1 (PLIN1), adiponectin C1Q and collagen domain containing (ADIPOQ), and lipoprotein lipase (LPL), were significantly downregulated (p < 0.01).ConclusionThe expression patterns of miRNAs exhibited significant fluctuations during different development periods of the fat tail, and some of them may participate in the regulation of tail fat deposition by modulating the proliferation and differentiation of preadipocytes.

Sainsc: A Computational Tool for Segmentation-Free Analysis of In Situ Capture Data.

Spatially resolved transcriptomics (SRT) has become the method of choice for characterising the complexity of biomedical tissue samples. Until recently, scientists were restricted to SRT methods that can profile a limited set of target genes at high spatial resolution or transcriptome-wide but at a low spatial resolution. Through recent developments, there are now methods that offer both subcellular spatial resolution and full transcriptome coverage. However, utilising these new methods' high spatial resolution and gene resolution remains elusive due to several factors, including low detection efficiency and high computational costs. Here, we present Sainsc (Segmentation-free analysis of in situ capture data), which combines a cell-segmentation-free approach with efficient data processing of transcriptome-wide nanometre-resolution spatial data. Sainsc can generate cell-type maps with accurate cell-type assignment at the nanometre scale, together with corresponding maps of the assignment scores that facilitate interpretation of the local confidence of cell-type assignment. We demonstrate its utility and accuracy for different tissues and technologies. Compared to other methods, Sainsc requires lower computational resources and has scalable performance, enabling interactive data exploration. Sainsc is compatible with common data analysis frameworks and is available as open-source software in multiple programming languages.

Identification of miR-451 target genes as prognostic markers in diffuse large B-cell lymphoma

ABSTRACT Background B-cell lymphoma, a diverse malignancy, is intricately regulated by multiple factors. MicroRNAs (miRNAs) have been demonstrated to be important regulators of the initiation and progression of human B-cell lymphoma, but their functions need to be further explored. Research design and methods Two B-cell lymphoma cell lines, Romas and HBL-1, were engineered to overexpress miR-451, with cell proliferation assessed via cell counting assays. Flow cytometry was used to study the effects of miR-451 on cell cycle and apoptosis. Bioinformatics analyses were performed using GEO datasets (GSE181063, GSE32918, GSE56315) to identify DEGs, and potential miR-451 target genes were predicted using tools like ENCORI, TargetScan, and R packages. A risk model for DLBCL prognosis was developed using Cox and LASSO regression. qRT-PCR validated the expression of these target genes. Results This study revealed that miR-451 inhibited cell proliferation, arrested the cell cycle, and induced apoptosis in human DLBCL cell lines. Bioinformatics analysis identified 9 target genes (MMP9, AQP9, RIN2, EOMES, LCP2, SELPLG, MAL, SOCS5, S1PR3) significantly associated with DLBCL prognosis, suggesting a potential mechanism by which miR-451 suppresses DLBCL development. Conclusions Our study indicates that a specific set of miR-451 target genes may significantly influence DLBCL patient outcomes.

Uncoupling FRUITFULL's functions through modification of a protein motif identified by co-ortholog analysis.

Many plant transcription factors (TFs) are multifunctional and regulate growth and development in more than one tissue. These TFs can generally associate with different protein partners depending on the tissue type, thereby regulating tissue-specific target gene sets. However, how interaction specificity is ensured is still largely unclear. Here, we examine protein-protein interaction specificity using subfunctionalized co-orthologs of the FRUITFULL (FUL) subfamily of MADS-domain TFs. In Arabidopsis, FUL is multifunctional, playing important roles in flowering and fruiting, whereas these functions have partially been divided in the tomato co-orthologs FUL1 and FUL2. By linking protein sequence and function, we discovered a key amino acid motif that determines interaction specificity of MADS-domain TFs, which in Arabidopsis FUL determines the interaction with AGAMOUS and SEPALLATA proteins, linked to the regulation of a subset of targets. This insight offers great opportunities to dissect the biological functions of multifunctional MADS TFs.

Network pharmacology to unveil the mechanism of Astragali Radix in the treatment of lupus nephritis via PI3K/AKT/mTOR pathway

We used network pharmacology, molecular docking, and in vitro experiments to explore the mechanisms of Astragali Radix in the treatment of lupus nephritis. We screened compounds and targets of Astragali Radix, as well as related genes of lupus nephritis from databases. We identified 211 common genes and 44 compounds between the herb and the disease, and constructed global, narrowed, hierarchical Compound-Target Interaction networks to illustrate the possible mechanism. We found that the PI3K/AKT/mTOR pathway is a core target gene set identified through enrichment analysis, PPI analysis and MCODE analysis. In vitro experiments showed that freeze-dried Astragali powder inhibits activation of PI3K, AKT1 and mTOR in TGF-β1 stimulated HK-2 cells. Molecular docking demonstrated that (R)-isomucronulatol, 3,9,10-trimethoxypterocarpan and astrapterocarpan exhibited promising binding affinity to PI3K, AKT, and mTOR proteins.

The transcription elongation factors Spt4 and Spt5 control neural progenitor proliferation and are implicated in neuronal remodeling during Drosophila mushroom body development.

Spt4 and Spt5 form the DRB sensitivity inducing factor (DSIF) complex that regulates transcription elongation at multiple steps including promotor-proximal pausing, processivity and termination. Although this implicated a general role in transcription, several studies pointed to smaller sets of target genes and indicated a more specific requirement in certain cellular contexts. To unravel common or distinct functions of Spt4 and Spt5 in vivo, we generated knock-out alleles for both genes in Drosophila melanogaster. Using the development of the mushroom bodies as a model, we provided evidence for two common functions of Spt4 and Spt5 during mushroom body development, namely control of cell proliferation of neural progenitor cells and remodeling of axonal projections of certain mushroom body neurons. This latter function is not due to a general requirement of Spt4 and Spt5 for axon pathfinding of mushroom body neurons, but due to distinct effects on the expression of genes controlling remodeling.

P53motifDB: integration of genomic information and tumor suppressor p53 binding motifs.

The tumor suppressor gene TP53 encodes the DNA binding transcription factor p53 and is one of the most commonly mutated genes in human cancer. Tumor suppressor activity requires binding of p53 to its DNA response elements and subsequent transcriptional activation of a diverse set of target genes. Despite decades of close study, the logic underlying p53 interactions with its numerous potential genomic binding sites and target genes is not yet fully understood. Here, we present a database of DNA and chromatin-based information focused on putative p53 binding sites in the human genome to allow users to generate and test new hypotheses related to p53 activity in the genome. Users can query genomic locations based on experimentally observed p53 binding, regulatory element activity, genetic variation, evolutionary conservation, chromatin modification state, and chromatin structure. We present multiple use cases demonstrating the utility of this database for generating novel biological hypotheses, such as chromatin-based determinants of p53 binding and potential cell type-specific p53 activity. All database information is also available as a precompiled sqlite database for use in local analysis or as a Shiny web application.

Pharmacological restriction of genomic binding sites redirects PU.1 pioneer transcription factor activity

Transcription factor (TF) DNA-binding dynamics govern cell fate and identity. However, our ability to pharmacologically control TF localization is limited. Here we leverage chemically driven binding site restriction leading to robust and DNA-sequence-specific redistribution of PU.1, a pioneer TF pertinent to many hematopoietic malignancies. Through an innovative technique, ‘CLICK-on-CUT&Tag’, we characterize the hierarchy of de novo PU.1 motifs, predicting occupancy in the PU.1 cistrome under binding site restriction. Temporal and single-molecule studies of binding site restriction uncover the pioneering dynamics of native PU.1 and identify the paradoxical activation of an alternate target gene set driven by PU.1 localization to second-tier binding sites. These transcriptional changes were corroborated by genetic blockade and site-specific reporter assays. Binding site restriction and subsequent PU.1 network rewiring causes primary human leukemia cells to differentiate. In summary, pharmacologically induced TF redistribution can be harnessed to govern TF localization, actuate alternate gene networks and direct cell fate.

Direct targets of MEF2C are enriched for genes associated with schizophrenia and cognitive function and are involved in neuron development and mitochondrial function.

Myocyte Enhancer Factor 2C (MEF2C) is a transcription factor that plays a crucial role in neurogenesis and synapse development. Genetic studies have identified MEF2C as a gene that influences cognition and risk for neuropsychiatric disorders, including autism spectrum disorder (ASD) and schizophrenia (SCZ). Here, we investigated the involvement of MEF2C in these phenotypes using human-derived neural stem cells (NSCs) and glutamatergic induced neurons (iNs), which represented early and late neurodevelopmental stages. For these cellular models, MEF2C function had previously been disrupted, either by direct or indirect mutation, and gene expression assayed using RNA-seq. We integrated these RNA-seq data with MEF2C ChIP-seq data to identify dysregulated direct target genes of MEF2C in the NSCs and iNs models. Several MEF2C direct target gene-sets were enriched for SNP-based heritability for intelligence, educational attainment and SCZ, as well as being enriched for genes containing rare de novo mutations reported in ASD and/or developmental disorders. These gene-sets are enriched in both excitatory and inhibitory neurons in the prenatal and adult brain and are involved in a wide range of biological processes including neuron generation, differentiation and development, as well as mitochondrial function and energy production. We observed a trans expression quantitative trait locus (eQTL) effect of a single SNP at MEF2C (rs6893807, which is associated with IQ) on the expression of a target gene, BNIP3L. BNIP3L is a prioritized risk gene from the largest genome-wide association study of SCZ and has a function in mitophagy in mitochondria. Overall, our analysis reveals that either direct or indirect disruption of MEF2C dysregulates sets of genes that contain multiple alleles associated with SCZ risk and cognitive function and implicates neuron development and mitochondrial function in the etiology of these phenotypes.

Genome binding properties of Zic transcription factors underlie their changing functions during neuronal maturation

BackgroundThe Zic family of transcription factors (TFs) promote both proliferation and maturation of cerebellar granule neurons (CGNs), raising the question of how a single, constitutively expressed TF family can support distinct developmental processes. Here we use an integrative experimental and bioinformatic approach to discover the regulatory relationship between Zic TF binding and changing programs of gene transcription during postnatal CGN differentiation.ResultsWe first established a bioinformatic pipeline to integrate Zic ChIP-seq data from the developing mouse cerebellum with other genomic datasets from the same tissue. In newborn CGNs, Zic TF binding predominates at active enhancers that are co-bound by developmentally regulated TFs including Atoh1, whereas in mature CGNs, Zic TF binding consolidates toward promoters where it co-localizes with activity-regulated TFs. We then performed CUT&RUN-seq in differentiating CGNs to define both the time course of developmental shifts in Zic TF binding and their relationship to gene expression. Mapping Zic TF binding sites to genes using chromatin looping, we identified the set of Zic target genes that have altered expression in RNA-seq from Zic1 or Zic2 knockdown CGNs.ConclusionsOur data show that Zic TFs are required for both induction and repression of distinct, developmentally regulated target genes through a mechanism that is largely independent of changes in Zic TF binding. We suggest that the differential collaboration of Zic TFs with other TF families underlies the shift in their biological functions across CGN development.

The JNK and Hippo pathways control epithelial integrity and prevent tumor initiation by regulating an overlapping transcriptome

Epithelial organs maintain their integrity and prevent tumor initiation by actively removing defective cells, such as those that have lost apicobasal polarity. Here, we identify how transcription factors of two key signaling pathways-Jun-N-terminal kinase (JNK) and Hippo-regulate epithelial integrity by controlling transcription of an overlapping set of target genes. Targeted DamID experiments reveal that, in proliferating cells of the Drosophila melanogaster eye, the AP-1 transcription factor Jun and the Hippo pathway transcription regulators Yorkie and Scalloped bind to a common suite of target genes that promote organ growth. In defective neoplastic cells, AP-1 transcription factors repress transcription of growth genes together with the C-terminal binding protein (CtBP) co-repressor. If gene repression by AP-1/CtBP fails, neoplastic tumor growth ensues, driven by Yorkie/Scalloped. Thus, AP-1/CtBP eliminates defective cells and prevents tumor initiation by acting in parallel to Yorkie/Scalloped to repress expression of a shared transcriptome. These findings shed new light on the maintenance of epithelial integrity and tumor suppression.

Angiogenesis related genes based prognostic model of glioma patients developed by multi-omics approach

IntroductionGlioma, particularly glioblastoma (GBM), is a highly malignant brain tumor with poor prognosis despite current therapeutic approaches. The tumor microenvironment (TME), plays a crucial role in glioma progression by promoting invasion and drug resistance. Angiogenesis, the formation of new blood vessels, is a tightly regulated process involving endothelial cell activation, proliferation, and migration. In cancer, angiogenesis becomes dysregulated, leading to excessive blood vessel formation.MethodsWe enrolled bulk data of TCGA-LGG/GBM, CGGA-693, and CGGA-325 cohorts, scRNA data of GSE162631, GSE84465, and GSE138794 cohorts. Identification of malignant cells was conducted by “copycat” R package. The “AUCell” R package scored the activity of target gene set of each single cell. Consensus clustering was applied using the “ConsensusClusterPlus” R package, while tumor-infiltrating immune cells were determined using “IOBR” R package. To construct a prognostic model, we used LASSO and multiCOX algorithms based on the expression levels of the 15 hub genes, the efficacy of which was verified by KM and ROC analysis.ResultsWe identified 4 different malignant cell subclusters in glioma and disclosed their distinct gene expression patterns and interactions within TME. We identified differentially expressed immune-related genes (DE-ARGs) in glioma and found 15 genes that were specifically expressed in the malignant glioma cell populations. Glioma cells with higher expression of these DE-ARGs were associated with gliogenesis, glial cell development, and vasculature development. We found that tumor-infiltrating monocytes were the main interacting cell type within glioma TME. Using the expression patterns of the 15 screened DE-ARGs, we categorized glioma samples into 2 molecular clusters with distinct immune features, suggesting a possible relationship between angiogenesis and immune activation and recruitment. We constructed a prognostic model based on the expression levels of the 15 DE-ARGs and evaluated its predictive ability for glioma patient outcomes, which displayed exceedingly high efficacy.ConclusionWe characterized different malignant cell subclusters in glioma and investigate their gene expression patterns and interactions within TME. We constructed a prognostic model based on the expression levels of the 15 DE-ARGs and evaluated its predictive ability for glioma patient outcomes, which displayed exceedingly high efficacy.

Body-Wide Inactivation of the Myc-Like Mlx Transcription Factor Network Accelerates Aging and Increases the Lifetime Cancer Incidence.

The "Mlx" and "Myc" transcription factor networks cross-communicate and share many common gene targets. Myc's activity depends upon its heterodimerization with Max, whereas the Mlx Network requires that the Max-like factor Mlx associate with the Myc-like factors MondoA or ChREBP. The current work demonstrates that body-wide Mlx inactivation, like that of Myc, accelerates numerous aging-related phenotypes pertaining to body habitus and metabolism. The deregulation of numerous aging-related Myc target gene sets is also accelerated. Among other functions, these gene sets often regulate ribosomal and mitochondrial structure and function, genomic stability, and aging. Whereas "MycKO" mice have an extended lifespan because of a lower cancer incidence, "MlxKO" mice have normal lifespans and a higher cancer incidence. Like Myc, the expression of Mlx, MondoA, and ChREBP and their control over their target genes deteriorate with age in both mice and humans. Collectively, these findings underscore the importance of lifelong and balanced cross-talk between the two networks to maintain proper function and regulation of the many factors that can affect normal aging.

Methynissolin confers protection against gastric carcinoma via targeting RIPK2

Gastric carcinoma (GC) is one of the most common malignant tumors worldwide, and ranks fifth in incidence and top 5 in tumor-associated mortality. Unfortunately, the molecular mechanisms contributing to GC development are still unclear, urgently requiring more studies. It has been proposed that the multifunctional regulator receptor-interacting protein 2 (RIPK2) is both a regulator that induces stomach cancer and its complications such as cachexia. However, the potential molecular function of RIPK2 on regulation of GC is largely unknown. Moreover, 3,9-di-O-methylnissolin, a flavonoid derived from the roots of Astragalus membranaceus, exhibits protective effects against cancer, inflammation, and bacterial infections. Unfortunately the antineoplastic mechanism of oral exposure of such flavonoid are poorly understood. It remains unknown whether 3,9-di-O-methylnissolin is capable of mitigating GC progression, and potential molecular mechanism. This study aims to provide mechanistic evidence that the 3,9-di-O-methylnissolin-triggered mitigation of GC by repressing RIPK2 activity and its downstream ASK1 signaling. The TCGA database, GEO database, subjects with GC phenotype and GC cell lines were used in this work. 3,9-di-O-methylnissolin intervention reduces stomach cancer cell migration, proliferation and invasion in vitro. Also, 3,9-di-O-methylnissolin confers protection against tumor growth in xenograft model in vivo. Mechanistically, on the other hand, in the para-carcinoma tissue, overexpression of RIPK2 reduces lipogenesis and disturbs lipid homeostasis. Overexpression of RIPK2 activates downstream signaling of ASK1, which in turn downregulates PPARα, a master regulator of lipid metabolism. This process results in lower lipogenesis and downregulation of a set of PPARα target genes. Astragalus membranaceus naturally produces 3,9-di-O-methylnissolin, a bioactive substance that when suppression of ASK1 by methylnissolin reduces RIPK2-induced GC and its related cachexia and restores lipid homeostasis. These findings offer important new understandings of stomach cancer and could aid in the development of treatment plans for gastric carcinoma and cachexia related to cancer.

A chronic signaling TGFb zebrafish reporter identifies immune response in melanoma.

Developmental signaling pathways associated with growth factors such as TGFb are commonly dysregulated in melanoma. Here we identified a human TGFb enhancer specifically activated in melanoma cells treated with TGFB1 ligand. We generated stable transgenic zebrafish with this TGFb Induced Enhancer driving green fluorescent protein (TIE:EGFP). TIE:EGFP was not expressed in normal melanocytes or early melanomas but was expressed in spatially distinct regions of advanced melanomas. Single-cell RNA-sequencing revealed that TIE:EGFP+ melanoma cells down-regulated interferon response while up-regulating a novel set of chronic TGFb target genes. ChIP-sequencing demonstrated that AP-1 factor binding is required for activation of chronic TGFb response. Overexpression of SATB2, a chromatin remodeler associated with tumor spreading, showed activation of TGFb signaling in early melanomas. Confocal imaging and flow cytometric analysis showed that macrophages localize to TIE:EGFP+ regions and preferentially phagocytose TIE:EGFP+ melanoma cells compared to TIE:EGFP- melanoma cells. This work identifies a TGFb induced immune response and demonstrates the need for the development of chronic TGFb biomarkers to predict patient response to TGFb inhibitors.

Exploring the multicomponent synergy mechanism of Zuogui Wan on postmenopausal osteoporosis by a systems pharmacology strategy.

To explore the multi-component synergistic mechanism of Zuogui Wan (, ZGW) in treating postmenopausal osteoporosis (PMOP). The main components and target genes of ZGW were screened via the Traditional Chinese Medicine Systems Pharmacology (TCMSP). In addition, the target gene sets of PMOP were derived from the GeneCards and Online Mendelian Inheritance in Man databases. The search tool for recurring instances of neighbouring genes (STRING) 11.0 software was used to analyze the interaction among intersecting genes. Cytoscape 3.6.1 software and the Matthews correlation coefficient (MCC) algorithm were used to screen the core genes. Fifty Sprague-Dawley female rats were randomly divided into the sham-operated (Sham) group and the four ovariectomized (OVX) subgroups. Rats subjected to Sham or OVX were administered with the vehicle (OVX, 1 mL water/100 g weight), 17β-estradiol (E2, 50 μg·kg-1·d-1), and lyophilized powder of ZGW at a low dose of 2.3 (ZGW-L) and high dose of 4.6 (ZGW-H) g·kg-1·d-1 for three months. The bone density and bone strength were assessed using dual-energy X-ray and three-point bending tests, respectively. Furthermore, enzyme-linked immun-osorbent assay, Hematoxylin-eosin staining, and western blot analysis were used to determine the potential pharmacological mechanisms of action of ZGW in PMOP. A total of 117 active compounds of ZGW were screened from the TCMSP. Furthermore, 108 intersecting genes of drugs and diseases were identified. Using STRING software and the MCC algorithm, ten core genes, including C-X-C chemokine living 8 (CXCL8), C-C chemokine receptor type 2 (CCR2), alpha-2a active receptor (ADRA2A), melatonin receptor type 1B (MTNR1B), and amyloid-beta A4 protein (APP), were identified. The anti-osteoporosis regulation network of ZGW was constructed using the Cytoscape software. The animal experiments demonstrated that ZGW groups significantly reduced the serum levels of β-C-terminal telopeptide of type I collagen (β-CTX) and increased serum levels of bone-specific alkaline phosphatase (BALP) (P < 0.05, P < 0.01). The OVX group exhibited a significant decrease in bone mineral density and bone strength compared with the Sham group (P < 0.01). Moreover, treatment with ZGW resulted in increased trabecular thickness, improved arrangement of trabecular structure, and reduced empty bone lacunae. Furthermore, treatment with ZGW significantly increased the protein expression of CXCL8, ADRA2A, and CCR2 (P < 0.05, P < 0.01), and significantly decreased the protein expression of Runx2 (P < 0.01). Furthermore, the ZGW and E2 groups demonstrated significantly increased BMD (P < 0.05, P < 0.01), improved bone strength (P < 0.05, P < 0.01), reduced expression of CXCL8, ADRA2A, and CCR2, and increased runt-related transcription factor 2 levels in bone tissue (P < 0.05, P < 0.01) compared with the OVX group. However, there were no significant differences in MTNR1B and APP expression among the groups. ZGW shows synergistic mechanisms in PMOP through multiple components, targets, and pathways.

HIF2α-dependent Dock4/Rac1-signaling regulates formation of adherens junctions and cell polarity in normoxia.

Hypoxia-inducible factors (HIF) 1 and 2 regulate similar but distinct sets of target genes. Although HIFs are best known for their roles in mediating the hypoxia response accumulating evidence suggests that under certain conditions HIFs, particularly HIF2, may function also under normoxic conditions. Here we report that HIF2α functions under normoxic conditions in kidney epithelial cells to regulate formation of adherens junctions. HIF2α expression was required to induce Dock4/Rac1/Pak1-signaling mediating stability and compaction of E-cadherin at nascent adherens junctions. Impaired adherens junction formation in HIF2α- or Dock4-deficient cells led to aberrant cyst morphogenesis in 3D kidney epithelial cell cultures. Taken together, we show that HIF2α functions in normoxia to regulate epithelial morphogenesis.

IrGSEA: the integration of single-cell rank-based gene set enrichment analysis.

irGSEA is an R package designed to assess the outcomes of various gene set scoring methods when applied to single-cell RNA sequencing data. This package incorporates six distinct scoring methods that rely on the expression ranks of genes, emphasizing relative expression levels over absolute values. The implemented methods include AUCell, UCell, singscore, ssGSEA, JASMINE and Viper. Previous studies have demonstrated the robustness of these methods to variations in dataset size and composition, generating enrichment scores based solely on the relative gene expression of individual cells. By employing the robust rank aggregation algorithm, irGSEA amalgamates results from all six methods to ascertain the statistical significance of target gene sets across diverse scoring methods. The package prioritizes user-friendliness, allowing direct input of expression matrices or seamless interaction with Seurat objects. Furthermore, it facilitates a comprehensive visualization of results. The irGSEA package and its accompanying documentation are accessible on GitHub (https://github.com/chuiqin/irGSEA).

A transcription factor-mediated regulatory network controls fungal pathogen colonization of insect body cavities.

Successful host tissue colonization is crucial for fungal pathogens to cause mycosis and complete the infection cycle, in which fungal cells undergo a series of morphological transition-included cellular events to combat with hosts. However, many transcription factors (TFs) and their mediated networks regulating fungal pathogen colonization of host tissue are not well characterized. Here, a TF (BbHCR1)-mediated regulatory network was identified in an insect pathogenic fungus, Beauveria bassiana, that controlled insect hemocoel colonization. BbHCR1 was highly expressed in fungal cells after reaching insect hemocoel and controlled the yeast (in vivo blastospores)-to-hyphal morphological switch, evasion of immune defense response, and fungal virulence. Comparative analysis of RNA sequencing and chromatin immunoprecipitation sequencing identified a core set of BbHCR1 target genes during hemocoel colonization, in which abaA and brlA were targeted to limit the rapid switch from blastospores to hyphae and fungal virulence. Two targets encoding hypothetical proteins, HP1 and HP2, were activated and repressed by BbHCR1, respectively, which acted as a virulence factor and repressor, respectively, suggesting that BbHCR1 activated virulence factors but repressed virulence repressors during the colonization of insect hemocoel. BbHCR1 tuned the expression of two dominant hemocoel colonization-involved metabolite biosynthetic gene clusters, which linked its regulatory role in evasion of immune response. Those functions of BbHCR1 were found to be collaboratively regulated by Fus3- and Hog1-MAP kinases via phosphorylation. These findings have drawn a regulatory network in which Fus3- and Hog1-MAP kinases phosphorylate BbHCR1, which in turn controls the colonization of insect body cavities by regulating fungal morphological transition and virulence-implicated genes.IMPORTANCEFungal pathogens adopt a series of tactics for successful colonization in host tissues, which include morphological transition and the generation of toxic and immunosuppressive molecules. However, many transcription factors (TFs) and their linked pathways that regulate tissue colonization are not well characterized. Here, we identified a TF (BbHCR1)-mediated regulatory network that controls the insect fungal pathogen, Beauveria bassiana, colonization of insect hemocoel. During these processes, BbHCR1 targeted the fungal central development pathway for the control of yeast (blastospores)-to-hyphae morphological transition, activated virulence factors, repressed virulence repressors, and tuned the expression of two dominant hemocoel colonization-involved immunosuppressive and immunostimulatory metabolite biosynthetic gene clusters. The BbHCR1 regulatory function was governed by Fus3- and Hog1-MAP kinases. These findings led to a new regulatory network composed of Fus3- and Hog1-MAP kinases and BbHCR1 that control insect body cavity colonization by regulating fungal morphological transition and virulence-implicated genes.

P53 rapidly restructures 3D chromatin organization to trigger a transcriptional response

Activation of the p53 tumor suppressor triggers a transcriptional program to control cellular response to stress. However, the molecular mechanisms by which p53 controls gene transcription are not completely understood. Here, we uncover the critical role of spatio-temporal genome architecture in this process. We demonstrate that p53 drives direct and indirect changes in genome compartments, topologically associating domains, and DNA loops prior to one hour of its activation, which escort the p53 transcriptional program. Focusing on p53-bound enhancers, we report 340 genes directly regulated by p53 over a median distance of 116 kb, with 74% of these genes not previously identified. Finally, we showcase that p53 controls transcription of distal genes through newly formed and pre-existing enhancer-promoter loops in a cohesin dependent manner. Collectively, our findings demonstrate a previously unappreciated architectural role of p53 as regulator at distinct topological layers and provide a reliable set of new p53 direct target genes that may help designs of cancer therapies.