Gene Regulatory Network Analysis Research Articles

MR-GGI: accurate inference of gene–gene interactions using Mendelian randomization

BackgroundResearchers have long studied the regulatory processes of genes to uncover their functions. Gene regulatory network analysis is one of the popular approaches for understanding these processes, requiring accurate identification of interactions among the genes to establish the gene regulatory network. Advances in genome-wide association studies and expression quantitative trait loci studies have led to a wealth of genomic data, facilitating more accurate inference of gene–gene interactions. However, unknown confounding factors may influence these interactions, making their interpretation complicated. Mendelian randomization (MR) has emerged as a valuable tool for causal inference in genetics, addressing confounding effects by estimating causal relationships using instrumental variables. In this paper, we propose a new statistical method, MR-GGI, for accurately inferring gene–gene interactions using Mendelian randomization.ResultsMR-GGI applies one gene as the exposure and another as the outcome, using causal cis-single-nucleotide polymorphisms as instrumental variables in the inverse-variance weighted MR model. Through simulations, we have demonstrated MR-GGI's ability to control type 1 error and maintain statistical power despite confounding effects. MR-GGI performed the best when compared to other methods using the F1 score on the DREAM5 dataset. Additionally, when applied to yeast genomic data, MR-GGI successfully identified six clusters. Through gene ontology analysis, we have confirmed that each cluster in our study performs distinct functional roles by gathering genes with specific functions.ConclusionThese findings demonstrate that MR-GGI accurately inferences gene–gene interactions despite the confounding effects in real biological environments.

Gene regulatory network analysis identifies MYL1, MDH2, GLS, and TRIM28 as the principal proteins in the response of mesenchymal stem cells to Mg2+ ions

Magnesium (Mg)-based implants have emerged as a promising alternative for orthopedic applications, owing to their bioactive properties and biodegradability. As the implants degrade, Mg2+ ions are released, influencing all surrounding cell types, especially mesenchymal stem cells (MSCs). MSCs are vital for bone tissue regeneration, therefore, it is essential to understand their molecular response to Mg2+ ions in order to maximize the potential of Mg-based biomaterials. In this study, we conducted a gene regulatory network (GRN) analysis to examine the molecular responses of MSCs to Mg2+ ions. We used time-series proteomics data collected at 11 time points across a 21-day period for the GRN construction. We studied the impact of Mg2+ ions on the resulting networks and identified the key proteins and protein interactions affected by the application of Mg2+ ions. Our analysis highlights MYL1, MDH2, GLS, and TRIM28 as the primary targets of Mg2+ ions in the response of MSCs during 1–21 days phase. Our results also identify MDH2-MYL1, MDH2-RPS26, TRIM28-AK1, TRIM28-SOD2, and GLS-AK1 as the critical protein relationships affected by Mg2+ ions. By offering a comprehensive understanding of the regulatory role of Mg2+ ions on MSCs, our study contributes valuable insights into the molecular response of MSCs to Mg-based materials, thereby facilitating the development of innovative therapeutic strategies for orthopedic applications.

The SoxE factor Sox9 is selectively expressed in indirect pathway striatal projection neurons and regulates synaptogenesis

Striatum, as the largest structure of the basal ganglia, serves as a center for information transmission and is critical for motor function and reward perception. However, the genetic mechanisms underlying its development require further exploration. Here, we found that Sox9,traditionally recognized as a glial marker, is uniquely expressed in striatal medium spiny neurons (MSNs), especially in Drd2-expressing indirect pathway MSNs (D2-MSNs). Intriguingly, Sox9 expression in the striatum, which is conserved in humans, is a dynamic process. It maintains a high level during the perinatal stage, and exhibits low expression levels or vanishes at the embryonic and postnatal stages, respectively. The peak period of Sox9 expression coincides with the transition from neurogenesis to synaptogenesis. Importantly, gene regulatory network analysis and gain-of-function experiments confirmed Sox9 is strongly correlated with synaptogenesis. Moreover, we identified that Sox9 regulates synaptogenesis by repressing Foxp2, a well-known synapse regulator. Furthermore, we demonstrated that the biased expression pattern of Sox9 in D2-MSNs is, at least in part, regulated by another SoxE family member Sox8, which is specifically expressed in Drd1-expressing direct pathway MSNs (D1-MSNs). Taken together, our findings reveal a new marker of D2-MSNs and identify its distinctive function in striatal development.

Adaptation in Unstable Environments and Global Gene Losses: Small but Stable Gene Networks by the May-Wigner Theory.

Although gene loss is common in evolution, it remains unclear whether it is an adaptive process. In a survey of seven major mangrove clades that are woody plants in the intertidal zones of daily environmental perturbations, we noticed that they generally evolved reduced gene numbers. We then focused on the largest clade of Rhizophoreae and observed the continual gene set reduction in each of the eight species. A great majority of gene losses are concentrated on environmental interaction processes, presumably to cope with the constant fluctuations in the tidal environments. Genes of the general processes for woody plants are largely retained. In particular, fewer gene losses are found in physiological traits such as viviparous seeds, high salinity, and high tannin content. Given the broad and continual genome reductions, we propose the May-Wigner theory (MWT) of system stability as a possible mechanism. In MWT, the most effective solution for buffering continual perturbations is to reduce the size of the system (or to weaken the total genic interactions). Mangroves are unique as immovable inhabitants of the compound environments in the land-sea interface, where environmental gradients (such as salinity) fluctuate constantly, often drastically. Extending MWT to gene regulatory network (GRN), computer simulations and transcriptome analyses support the stabilizing effects of smaller gene sets in mangroves vis-à-vis inland plants. In summary, we show the adaptive significance of gene losses in mangrove plants, including the specific role of promoting phenotype innovation and a general role in stabilizing GRN in unstable environments as predicted by MWT.

Genome-wide identification and in silico characterization of major RNAi gene families in date palm (Phoenix dactylifera)

BackgroundDates contain various minerals that are essential for good health. The major RNA interference (RNAi) gene families play a vital role in plant growth and development by controlling the expression of protein-coding genes against different biotic and abiotic stresses. However, these gene families for date palm are not yet studied. Therefore, this study has explored major RNAi genes and their characteristics in date palm.ResultsWe have identified 4 PdDCLs, 7 PdAGOs, and 3 PdRDRs as RNAi proteins from the date palm genome by using AtRNAi genes as query sequences in BLASTp search. Domain analysis of predicted RNAi genes has revealed the Helicase_C, Dicer_dimer, PAZ, RNase III, and Piwi domains that are associated with the gene silencing mechanisms. Most PdRNAi proteins have been found in the nucleus and cytosol associated with the gene silencing actions. The gene ontology (GO) enrichment analysis has revealed some important GO terms including RNA interference, dsRNA fragmentation, and ribonuclease_III activity that are related to the protein-coding gene silencing mechanisms. Gene regulatory network (GRN) analysis has identified PAZ and SNF2 as the transcriptional regulators of PdRNAi genes. Top-ranked 10 microRNAs including Pda-miR156b, Pda-miR396a, Pda-miR166a, Pda-miR167d, and Pda-miR529a have been identified as the key post-transcriptional regulators of PdRNAi genes that are associated with different biotic/abiotic stresses. The cis-acting regulatory element analysis of PdRNAi genes has detected some vital cis-acting elements including ABRE, MBS, MYB, MYC, Box-4, G-box, I-box, and STRE that are linked with different abiotic stresses.ConclusionThe results of this study might be valuable resources for the improvement of different characteristics in date palm by further studies in wet-lab.

Differential Gene Regulatory Network Analysis between Azacitidine-Sensitive and -Resistant Cell Lines.

Azacitidine, a DNA methylation inhibitor, is employed for the treatment of acute myeloid leukemia (AML). However, drug resistance remains a major challenge for effective azacitidine chemotherapy, though several studies have attempted to uncover the mechanisms of azacitidine resistance. With the aim to identify the mechanisms underlying acquired azacitidine resistance in cancer cell lines, we developed a computational strategy that can identify differentially regulated gene networks between drug-sensitive and -resistant cell lines by extending the existing method, differentially coexpressed gene sets (DiffCoEx). The technique specifically focuses on cell line-specific gene network analysis. We applied our method to gene networks specific to azacitidine sensitivity and identified differentially regulated gene networks between azacitidine-sensitive and -resistant cell lines. The molecular interplay between the metallothionein gene family, C19orf33, ELF3, GRB7, IL18, NRN1, and RBM47 were identified as differentially regulated gene network in drug resistant cell lines. The biological mechanisms associated with azacitidine and AML for the markers in the identified networks were verified through the literature. Our results suggest that controlling the identified genes (e.g., the metallothionein gene family) and "cellular response"-related pathways ("cellular response to zinc ion", "cellular response to copper ion", and "cellular response to cadmium ion", where the enriched functional-related genes are MT2A, MT1F, MT1G, and MT1E) may provide crucial clues to address azacitidine resistance in patients with AML. We expect that our strategy will be a useful tool to uncover patient-specific molecular interplay that provides crucial clues for precision medicine in not only gastric cancer but also complex diseases.

Enabling comparative gene regulatory network analysis on single-cell data with SCORPION.

Enabling comparative gene regulatory network analysis on single-cell data with SCORPION.

A Sorghum F-Box Protein Induces an Oxidative Burst in the Defense Against Colletotrichum sublineola.

The hemibiotrophic fungal pathogen Colletotrichum sublineola is the causal agent of anthracnose in sorghum (Sorghum bicolor), resulting in leaf blight, stalk rot, and head blight in susceptible genotypes, with yield losses of up to 50%. The development of anthracnose-resistant cultivars can reduce reliance on fungicides and provide a more sustainable and economical means for disease management. A previous genome-wide association study of the sorghum association panel identified the candidate resistance gene Sobic.005G172300 encoding an F-box protein. To better understand the role of this gene in the defense against C. sublineola, gene expression following infection with C. sublineola was monitored by RNA sequencing in seedlings of sorghum accession SC110, which harbored the resistance allele, and three accessions that harbored a susceptible allele. Only in SC110 did the expression of Sobic.005G172300 increase during the biotrophic phase of infection. Subsequent transcriptome analysis, gene co-expression networks, and gene regulatory networks of inoculated and mock-inoculated seedlings of resistant and susceptible accessions suggest that the increase in expression of Sobic.005G172300 induces an oxidative burst by lowering the concentration of ascorbic acid during the biotrophic phase of infection. Based on gene regulatory network analysis, the protein encoded by Sobic.005G172300 is proposed to target proteins involved in the biosynthesis of ascorbic acid for polyubiquitination through the SCF E3 ubiquitin ligase, causing their degradation via the proteasome.

Immune characteristics of kidney transplant recipients with acute respiratory distress syndrome induced by COVID-19 at single-cell resolution

BackgroundCOVID-19-induced acute respiratory distress syndrome (ARDS) can result in tissue damage and multiple organ dysfunction, especially in kidney transplant recipients (KTRs) receiving immunosuppressive drugs. Presently, single-cell research on COVID-19-induced ARDS is considerably advanced, yet knowledge about ARDS in KTRs is still constrained.MethodsSingle-cell RNA sequencing (scRNA-seq) analysis was performed to construct a comprehensive single-cell immune landscape of the peripheral blood mononuclear cells (PBMCs) of eight patients with COVID-19-induced ARDS, five KTRs with COVID-19-induced ARDS, and five healthy individuals. Subsequently, we conducted a comprehensive bioinformatics analysis, including cell clustering, enrichment analysis, trajectory analysis, gene regulatory network analysis, and cell–cell interaction analysis, to investigate the heterogeneity of the immune microenvironment in KTRs with ARDS.ResultOur study revealed that KTRs exhibit significant heterogeneity with COVID-19-induced ARDS compared with those of other individuals, with significant reductions in T cells, as well as an abnormal proliferation of B cells and monocytes. In the context of dual influences from immunosuppression and viral infection, KTRs exhibited more specific plasma cells, along with significant enrichment of dysfunctional GZMB and XAF1 double-positive effector T cells and IFI27-positive monocytes. Additionally, robust communication existed among T cells and monocytes in cytokine signaling. These effects impede the process of immune reconstitution in KTR patients.ConclusionOur findings suggest that KTRs with COVID-19-induced ARDS show elevated antibody levels, impaired T cell differentiation, and dysregulation of innate immunity. In summary, this study provides a theoretical foundation for a comprehensive understanding of COVID-19-induced ARDS in KTRs.

Key Genes in Olfactory Disorder in Experimental Autoimmune Encephalomyelitis Identified by Transcriptomic Analysis of the Olfactory Bulbs.

Experimental autoimmune encephalomyelitis (EAE) is an animal model of multiple sclerosis that shows demyelination in the central nervous system and functional deficits, including olfactory impairment. However, the genes related to olfactory impairment in EAE are unknown. We evaluated hub genes of the olfactory bulb in EAE mice. Differentially expressed genes (cut-offs, fold change > 2 and adjusted p < 0.05) and their related pathways in olfactory bulbs were subjected to gene ontology (GO) pathway analysis, gene set enrichment analysis (GSEA). Protein-protein interactions with selected genes were evaluated using the Search Tool for the Retrieval of Interacting Genes/Proteins. Gene regulatory networks (GRNs) which were constructed at the post-transcriptional level, including the genes-transcription factors (TFs) and gene-microRNAs (miRNAs) interaction networks. Twelve hub genes were found, three of which (Ctss, Itgb2, and Tlr2) were validated by RT-qPCR to be related to GO pathways such as immune response and regulation of immune response. GSEA showed that neuron-related genes-including Atp6v1g2, Egr1, and Gap43-and their pathways were significantly downregulated. GRNs analysis of six genes (Ctss, Itgb2, Tlr2, Atp6v1g2, Egr1, and Gap43) revealed 37 TFs and 84 miRNAs were identified as potential regulators of six genes, indicating significant interaction among six genes, TFs, and miRNAs. Collectively, these results suggest that transcriptomic analysis of the olfactory bulb of EAE mice can provide insight into olfactory dysfunction and reveal therapeutic targets for olfactory impairment.

The remission status of AML patients after allo-HCT is associated with a distinct single-cell bone marrow T-cell signature

AbstractAcute myeloid leukemia (AML) is a hematologic malignancy for which allogeneic hematopoietic cell transplantation (allo-HCT) often remains the only curative therapeutic approach. However, incapability of T cells to recognize and eliminate residual leukemia stem cells might lead to an insufficient graft-versus-leukemia (GvL) effect and relapse. Here, we performed single-cell RNA-sequencing (scRNA-seq) on bone marrow (BM) T lymphocytes and CD34+ cells of 6 patients with AML 100 days after allo-HCT to identify T-cell signatures associated with either imminent relapse (REL) or durable complete remission (CR). We observed a higher frequency of cytotoxic CD8+ effector and gamma delta (γδ) T cells in CR vs REL samples. Pseudotime and gene regulatory network analyses revealed that CR CD8+ T cells were more advanced in maturation and had a stronger cytotoxicity signature, whereas REL samples were characterized by inflammatory tumor necrosis factor/NF-κB signaling and an immunosuppressive milieu. We identified ADGRG1/GPR56 as a surface marker enriched in CR CD8+ T cells and confirmed in a CD33-directed chimeric antigen receptor T cell/AML coculture model that GPR56 becomes upregulated on T cells upon antigen encounter and elimination of AML cells. We show that GPR56 continuously increases at the protein level on CD8+ T cells after allo-HCT and confirm faster interferon gamma (IFN-γ) secretion upon re-exposure to matched, but not unmatched, recipient AML cells in the GPR56+ vs GPR56– CD8+ T-cell fraction. Together, our data provide a single-cell reference map of BM–derived T cells after allo-HCT and propose GPR56 expression dynamics as a surrogate for antigen encounter after allo-HCT.

Use of gene regulatory network analysis to repurpose drugs to treat bipolar disorder

BackgroundBipolar disorder (BD) presents significant challenges in drug discovery, necessitating alternative approaches. Drug repurposing, leveraging computational techniques and expanding biomedical data, holds promise for identifying novel treatment strategies. MethodsThis study utilized gene regulatory networks (GRNs) to identify significant regulatory changes in BD, using network-based signatures for drug repurposing. Employing the PANDA algorithm, we investigated the variations in transcription factor-GRNs between individuals with BD and unaffected individuals, incorporating binding motifs, protein interactions, and gene co-expression data. The differences in edge weights between BD and controls were then used as differential network signatures to identify drugs potentially targeting the disease-associated gene signature, employing the CLUEreg tool in the GRAND database. ResultsUsing a large RNA-seq dataset of 216 post-mortem brain samples from the CommonMind consortium, we constructed GRNs based on co-expression for individuals with BD and unaffected controls, involving 15,271 genes and 405 TFs. Our analysis highlighted significant influences of these TFs on immune response, energy metabolism, cell signalling, and cell adhesion pathways in the disorder. By employing drug repurposing, we identified 10 promising candidates potentially repurposed as BD treatments. LimitationsNon-drug-naïve transcriptomics data, bulk analysis of BD samples, potential bias of GRNs towards well-studied genes. ConclusionsFurther investigation into repurposing candidates, especially those with preclinical evidence supporting their efficacy, like kaempferol and pramocaine, is warranted to understand their mechanisms of action and effectiveness in treating BD. Additionally, novel targets such as PARP1 and A2b offer opportunities for future research on their relevance to the disorder.

A conserved graft formation process in Norway spruce and Arabidopsis identifies the PAT gene family as central regulators of wound healing

The widespread use of plant grafting enables eudicots and gymnosperms to join with closely related species and grow as one. Gymnosperms have dominated forests for over 200 million years, and despite their economic and ecological relevance, we know little about how they graft. Here we developed a micrografting method in conifers using young tissues that allowed efficient grafting with closely related species and between distantly related genera. Conifer graft junctions rapidly connected vasculature and differentially expressed thousands of genes including auxin and cell-wall-related genes. By comparing these genes to those induced during Arabidopsis thaliana graft formation, we found a common activation of cambium, cell division, phloem and xylem-related genes. A gene regulatory network analysis in Norway spruce (Picea abies) predicted that PHYTOCHROME A SIGNAL TRANSDUCTION 1 (PAT1) acted as a core regulator of graft healing. This gene was strongly up-regulated during both spruce and Arabidopsis grafting, and Arabidopsis mutants lacking PAT genes failed to attach tissues or successfully graft. Complementing Arabidopsis PAT mutants with the spruce PAT1 homolog rescued tissue attachment and enhanced callus formation. Together, our data show an ability for young tissues to graft with distantly related species and identifies the PAT gene family as conserved regulators of graft healing and tissue regeneration.

Deciphering cell fates of the human cornea with a multi‐omics meta atlas

The human cornea is a transparent tissue composed of several layers that allow light to enter the eye. Disruption of the function of various cells in different layers can give rise to diseases that affect the transparency of the cornea. Understanding the precise cell fates of the cornea and the underlying molecular mechanisms is pivotal for dissecting the pathology of corneal diseases and exploring treatments and regeneration. Currently, there is a lack of a precise cell fate map for the human cornea. In this study, we constructed a multi‐omics meta‐atlas to define cell fates in the human adult and embryonic cornea, by integrating publicly available datasets. This multi‐omics meta‐atlas is comprised of RNA‐seq data representing genome‐wide gene expression profiles and ATAC‐seq data annotating regulatory genomic regions that control gene expression in single‐cells of the cornea. We identified the key transcription factors and combinations that determine cell fates in the human cornea, with gene regulatory network analyses using the meta‐atlas. Some of these transcription factors are functionally validated. We also applied this meta‐atlas to analyse single cells of the cornea derived from keratoconus patients and showed an expansion of corneal stromal stem cells as a potential disease mechanism. In conclusion, this study provides an extensive reference of cell fates of the human cornea that can be used to dissect disease mechanisms. This meta‐atlas also provides insights for therapeutic discovery and molecular blueprints for corneal regeneration.

Bioinformatics and system biology approaches for identifying potential therapeutic targets for prostate cancer

Prostate cancer (PCa) is the primary cause of cancer-related deaths among men, and its incidence increases with age. Despite the availability of treatments, the high costs and complications associated with current cancer therapies pose challenges. Therefore, there is an urgent need to explore globally accessible and effective alternative medications to support prostate cancer treatment. Bioinformatics and system biology approaches were utilized to investigate repurposed medicines for prostate cancer treatment guided by host genomic biomarkers. The investigation identified 120 common differentially expressed genes (DEGs) and 24 hub genes (HubGs) as potential therapeutic targets by analyzing gene expression patterns and protein-protein interaction (PPI) networks. The top 10 hub genes CDK1, CCNA2, RRM2, ASPM, MELK, HJURP, DLGAP5, NCAPG, TTK, and HMMR are associated with the development and progression of prostate cancer. The DEGs enrichment analysis revealed significant axon guidance, cell-cell adhesion, extracellular space, and calcium signaling pathways associated with prostate cancer. Gene regulatory network analysis revealed transcription factors (FOXC1, GATA2 and E2F1 etc.) and micro-RNAs (hsa-mir-16–5p, hsa-mir-24–3p and hsa-mir-34a-5p etc.) as regulators of hub genes. The potential drug candidates, includingindirubin-3′-Monoxime (−7.7 kcal/mol), benzene sulfonamide (−8.8 kcal/mol), and cladribine (−6.6 kcal/mol) have been identified as the top 3 candidate drugs against CDK1, CCNA2 and RRM2, respectively. Besides, these drugs exhibit strong binding affinity with the target proteins and suggest their potential for repurposing in the treatment of prostate cancer. Collectively, the findings present valuable insights for advancing prostate cancer research, potentially leading to a deeper understanding and novel treatment strategies for relevant diseases.

Ascorbate oxidases in bread wheat: gene regulatory network, transcripts profiling, and interaction analyses provide insight into their role in plant development and stress response

Ascorbate oxidases in bread wheat: gene regulatory network, transcripts profiling, and interaction analyses provide insight into their role in plant development and stress response

A telomere-to-telomere reference genome of ficus (Ficus hispida) provides new insights into sex determination.

A high-quality reference genome is indispensable for resolving biologically essential traits. Ficus hispida is a dioecious plant. A complete Ficus reference genome will be crucial for understanding their sex evolution and important biological characteristics, such as aerial roots, mutualistic symbiosis with ficus-wasps, and fruiting from old stems. Here, we generated a telomere-to-telomere (T2T) genome for F. hispida using PacBio HiFi and Oxford Nanopore Ultra-long sequencing technologies. The genome contiguity and completeness has shown improvement compared with the previously released genome, with the annotation of six centromeres and 28 telomeres. We have refined our previously reported 2-Mb male-specific region into a 7.2-Mb genomic region containing 51 newly predicted genes and candidate sex-determination genes AG2 and AG3. Many of these genes showed extremely low expression, likely attributed to hypermethylation in the gene body and promoter regions. Gene regulatory networks (GRNs) revealed that AG2 and AG3 are related to the regulation of stamen development in male flowers, while the AG1 gene is responsible for regulating female flowers' defense responses and secondary metabolite processes. Comparative analysis of GRNs showed that the NAC, WRKY, and MYB transcription factor families dominate the female GRN, whereas the MADS and MYB transcription factor families are prevalent in the male GRN.

Systematic single cell RNA sequencing analysis reveals unique transcriptional regulatory networks of Atoh1-mediated hair cell conversion in adult mouse cochleae

Regeneration of mammalian cochlear hair cells (HCs) by modulating molecular pathways or transcription factors is a promising approach to hearing restoration; however, immaturity of the regenerated HCs in vivo remains a major challenge. Here, we analyzed a single cell RNA sequencing (scRNA-seq) dataset during Atoh1-induced supporting cell (SC) to hair cell (HC) conversion in adult mouse cochleae (Yamashita et al. (2018)) using multiple high-throughput sequencing analytical tools (WGCNA, SCENIC, ARACNE, and VIPER). Instead of focusing on differentially expressed genes, we established independent expression modules and confirmed the existence of multiple conversion stages. Gene regulatory network (GRN) analysis uncovered previously unidentified key regulators, including Nhlh1, Lhx3, Barhl1 and Nfia, that guide converted HC differentiation. Comparison of the late-stage converted HCs with the scRNA-seq data from neonatal mouse cochleae (Kolla et al. (2020)) revealed that they closely resemble postnatal day 1 wild-type OHCs, in contrast to other developmental stages. Using ARACNE and VIPER, we discovered multiple key regulators likely to promote conversion to a more mature OHC-like state, including Zbtb20, Nfia, Zmiz1, Gm14418, Bhlhe40, Six2, Fosb and Klf9. Our findings provide insights into the regulation of HC regeneration in adult mammalian cochleae in vivo and demonstrate an approach for analyzing GRNs in large scRNA-seq datasets.

Viral integration promotes SV40T-induced immortalization by disturbing the expression of DNA/chromosome- and ECM-associated functional genes

Lentivirus containing simian virus 40 large T antigen (SV40T) is routinely used to induce cell immortalization. However, the roles of viral integration itself in this progress is still controversial. Here, we transformed primary mouse embryonic fibroblasts (MEFs) with SV40T lentivirus and studied the roles of viral integration in the immortalization using RNA sequencing (RNA-seq) and whole genome sequencing (WGS). During the immortalization, differentially expressed genes (DGEs) are enriched in viral infection and several diverse activities. However, DEGs between immortalized and aging cells are significantly enriched in DNA/chromosome- and extracellular matrix (ECM)-associated activities. Gene regulatory network (GRN) analysis shows that although p53 is a key regulatory factor, many other transcription factors also play critical roles in the process, like STAT1. Of these DEGs, 32 genes have viral integration in their coding and/or regulatory regions. Our findings suggest that viral integration may promote SV40T-mediated immortalization by disturbing the expression of DNA/chromosome- and ECM-associated genes.

Gene regulatory network analysis predicts cooperating transcription factor regulons required for FLT3-ITD+ AML growth

Acute myeloid leukemia (AML) is a heterogeneous disease caused by different mutations. Previously, we showed that each mutational subtype develops its specific gene regulatory network (GRN) with transcription factors interacting within multiple gene modules, many of which are transcription factor genes themselves. Here, we hypothesize that highly connected nodes within such networks comprise crucial regulators of AML maintenance. We test this hypothesis using FLT3-ITD-mutated AML as a model and conduct an shRNA drop-out screen informed by this analysis. We show that AML-specific GRNs predict crucial regulatory modules required for AML growth. Furthermore, our work shows that all modules are highly connected and regulate each other. The careful multi-omic analysis of the role of one (RUNX1) module by shRNA and chemical inhibition shows that this transcription factor and its target genes stabilize the GRN of FLT3-ITD+ AML and that its removal leads to GRN collapse and cell death.