AltAnalyze and DomainGraph: analyzing and visualizing exon expression data

DDX41 haploinsufficiency causes inefficient hematopoiesis under stress and cooperates with p53 mutations to cause hematologic malignancy

Germline heterozygous mutations in DDX41 predispose individuals to hematologic malignancies in adulthood. Most of these DDX41 mutations result in a truncated protein, leading to loss of protein function. To investigate the impact of these mutations on hematopoiesis, we generated mice with hematopoietic-specific knockout of one Ddx41 allele. Under normal steady-state conditions, there was minimal effect on lifelong hematopoiesis, resulting in a mild yet persistent reduction in red blood cell counts. However, stress induced by transplantation of the Ddx41+/- BM resulted in hematopoietic stem/progenitor cell (HSPC) defects and onset of hematopoietic failure upon aging. Transcriptomic analysis of HSPC subsets from the transplanted BM revealed activation of cellular stress responses, including upregulation of p53 target genes in erythroid progenitors. To understand how the loss of p53 affects the phenotype of Ddx41+/- HSPCs, we generated mice with combined Ddx41 and Trp53 heterozygous deletions. The reduction in p53 expression rescued the fitness defects in HSPC caused by Ddx41 heterozygosity. However, the combined Ddx41 and Trp53 mutant mice were prone to developing hematologic malignancies that resemble human myelodysplastic syndrome and acute myeloid leukemia. In conclusion, DDX41 heterozygosity causes dysregulation of the response to hematopoietic stress, which increases the risk of transformation with a p53 mutation.

Proximity labeling proteomics reveals Kv1.3 potassium channel immune interactors in microglia

Microglia are resident immune cells of the brain and regulate its inflammatory state. In neurodegenerative diseases, microglia transition from a homeostatic state to a state referred to as disease associated microglia (DAM). DAM express higher levels of proinflammatory signaling molecules, like STAT1 and TLR2, and show transitions in mitochondrial activity toward a more glycolytic response. Inhibition of Kv1.3 decreases the proinflammatory signature of DAM, though how Kv1.3 influences the response is unknown. Our goal was to identify the potential proteins interacting with Kv1.3 during transition to DAM. We utilized TurboID, a biotin ligase, fused to Kv1.3 to evaluate potential interacting proteins with Kv1.3 via mass spectrometry in BV-2 microglia following TLR4-mediated activation. Electrophysiology, western blotting, and flow cytometry were used to evaluate Kv1.3 channel presence and TurboID biotinylation activity. We hypothesized that Kv1.3 contains domain-specific interactors that vary during a TLR4-induced inflammatory response, some of which are dependent on the PDZ-binding domain on the C-terminus. We determined that the N-terminus of Kv1.3 is responsible for trafficking Kv1.3 to the cell surface and mitochondria (e.g. NUDC, TIMM50). Whereas, the C-terminus interacts with immune signaling proteins in an LPS-induced inflammatory response (e.g. STAT1, TLR2, and C3). There are 70 proteins that rely on the C-terminal PDZ-binding domain to interact with Kv1.3 (e.g. ND3, Snx3, and Sun1). Furthermore, we used Kv1.3 blockade to verify functional coupling between Kv1.3 and interferon-mediated STAT1 activation. Overall, we highlight that the Kv1.3 potassium channel functions beyond conducting the outward flux of potassium ions in an inflammatory context and that Kv1.3 modulates the activity of key immune signaling proteins, such as STAT1 and C3.

Egr2 drives the differentiation of Ly6Chi monocytes into fibrosis-promoting macrophages in metabolic dysfunction-associated steatohepatitis in mice

Metabolic dysfunction-associated steatohepatitis (MASH), previously called non-alcoholic steatohepatitis (NASH), is a growing concern worldwide, with liver fibrosis being a critical determinant of its prognosis. Monocyte-derived macrophages have been implicated in MASH-associated liver fibrosis, yet their precise roles and the underlying differentiation mechanisms remain elusive. In this study, we unveil a key orchestrator of this process: long chain saturated fatty acid-Egr2 pathway. Our findings identify the transcription factor Egr2 as the driving force behind monocyte differentiation into hepatic lipid-associated macrophages (hLAMs) within MASH liver. Notably, Egr2-deficiency reroutes monocyte differentiation towards a macrophage subset resembling resident Kupffer cells, hampering hLAM formation. This shift has a profound impact, suppressing the transition from benign steatosis to liver fibrosis, demonstrating the critical pro-fibrotic role played by hLAMs in MASH pathogenesis. Long-chain saturated fatty acids that accumulate in MASH liver emerge as potent inducers of Egr2 expression in macrophages, a process counteracted by unsaturated fatty acids. Furthermore, oral oleic acid administration effectively reduces hLAMs in MASH mice. In conclusion, our work not only elucidates the intricate interplay between saturated fatty acids, Egr2, and monocyte-derived macrophages but also highlights the therapeutic promise of targeting the saturated fatty acid-Egr2 axis in monocytes for MASH management.

Presurgery health influences outcomes following vertical sleeve gastrectomy in adolescents

OBJECTIVE

Weight loss following vertical sleeve gastrectomy (VSG) in youth can range from 10% to 50%. We examined whether there are differences in demographic or metabolic parameters before VSG in youth who achieve above-average weight loss (AAWL) versus below-average weight loss (BAWL) at 1 year post VSG and if youth with BAWL still achieve metabolic health improvements at 1 year post VSG.

METHODS

Demographic, anthropometric, and clinical lab data were collected before VSG and at 1, 3, 6, and 12 months after VSG.

RESULTS

Forty-three youth with a mean age of 16.9 (SD 1.7) years before VSG were studied; 70% were female, 19% non-Hispanic Black, 58% non-Hispanic White, and 23% mixed/other race. Mean baseline BMI was 51.1 (SD 10.5) kg/m2. Average weight loss was 25.8%. The AAWL group lost 18.6 kg/m2 (35.3%) versus the BAWL group, who lost 8.8 kg/m2 (17.5%). BMI, age, race, sex, and socioeconomic status at baseline were similar between AAWL and BAWL groups; however, the BAWL group had a higher frequency of pre-VSG dysglycemia, steatotic liver disease, and dyslipidemia. At 1 year post VSG, fewer youth in the BAWL group achieved ideal health parameters, and they had less resolution of comorbidities.

CONCLUSIONS

The presence of comorbidities before VSG is associated with less weight loss and reduced resolution of metabolic conditions at 1 year post VSG.

Human cytomegalovirus infection coopts chromatin organization to diminish TEAD1 transcription factor activity

Human cytomegalovirus (HCMV) infects up to 80% of the world's population. Here, we show that HCMV infection leads to widespread changes in human chromatin accessibility and chromatin looping, with hundreds of thousands of genomic regions affected 48 hours after infection. Integrative analyses reveal HCMV-induced perturbation of Hippo signaling through drastic reduction of TEAD1 transcription factor activity. We confirm extensive concordant loss of TEAD1 binding, active H3K27ac histone marks, and chromatin looping interactions upon infection. Our data position TEAD1 at the top of a hierarchy involving multiple altered important developmental pathways. HCMV infection reduces TEAD1 activity through four distinct mechanisms: closing of TEAD1-bound chromatin, reduction of YAP1 and phosphorylated YAP1 levels, reduction of TEAD1 transcript and protein levels, and alteration of TEAD1 exon-6 usage. Altered TEAD1-based mechanisms are highly enriched at genetic risk loci associated with eye and ear development, providing mechanistic insight into HCMV's established roles in these processes.

ASpedia-R: a package to retrieve junction-incorporating features and knowledge-based functions of human alternative splicing events

Motivation

Alternative splicing (AS) is a key regulatory mechanism that confers genetic diversity and phenotypic plasticity of human. The exons and their flanking regions include comprehensive junction-incorporating sequence features like splicing factor-binding sites and protein domains. These elements involve in exon usage and finally contribute to isoform-specific biological functions. Splicing-associated sequence features are involved in the multilayered regulation encompassing DNA and proteins. However, most analysis applications have investigated limited sequence features, like protein domains. It is insufficient to explain the comprehensive cause and effect of exon-specific biological processes.

Results

With the advent of RNA-seq technology, global AS event analysis has deduced more precise results. As accumulating analysis results, it could be a challenge to identify multi-omics sequence features for AS events. Therefore, application to investigate multi-omics sequence features is useful to scan critical evidence. ASpedia-R is an R package to interrogate junction-incorporating sequence features for human genes. Our database collected the heterogeneous profile encompassed from DNA to protein. Additionally, knowledge-based splicing genes were collected using text-mining to test the association with specific pathway terms. Our package retrieves AS events for high-throughput data analysis results via AS event ID converter. Finally, result profile could be visualized and saved to multiple formats: sequence feature result table, genome track figure, protein-protein interaction network, and gene set enrichment test result table. Our package is a convenient tool to understand global regulation mechanisms by splicing.

Availability and implementation

The package source code is freely available to non-commercial users at https://github.com/ncc-bioinfo/ASpedia-R.

Identification and validation of putative biomarkers by in silico analysis, mRNA expression and oxidative stress indicators for negative energy balance in buffaloes during transition period

OBJECTIVE

Transition period is considered from 3 weeks prepartum to 3 weeks postpartum, characterized with dramatic events (endocrine, metabolic, and physiological) leading to occurrence of production diseases (negative energy balance/ketosis, milk fever etc). The objectives of our study were to analyze the periodic concentration of serum beta-hydroxy butyric acid (BHBA), glucose and oxidative markers along with identification, and validation of the putative markers of negative energy balance in buffaloes using in-silico and quantitative real time-polymerase chain reaction (qRT-PCR) assay.

METHODS

Out of 20 potential markers of ketosis identified by in-silico analysis, two were selected and analyzed by qRT-PCR technique (upregulated; acetyl serotonin o-methyl transferase like and down regulated; guanylate cyclase activator 1B). Additional two sets of genes (carnitine palmotyl transferase A; upregulated and Insulin growth factor; downregulated) that have a role of hepatic fatty acid oxidation to maintain energy demands via gluconeogenesis were also validated. Extracted cDNA (complementary deoxyribonucleic acid) from the blood of the buffaloes were used for validation of selected genes via qRTPCR. Concentrations of BHBA, glucose and oxidative stress markers were identified with their respective optimized protocols.

RESULTS

The analysis of qRT-PCR gave similar trends as shown by in-silico analysis throughout the transition period. Significant changes (p<0.05) in the levels of BHBA, glucose and oxidative stress markers throughout this period were observed. This study provides validation from in-silico and qRT-PCR assays for potential markers to be used for earliest diagnosis of negative energy balance in buffaloes.

CONCLUSION

Apart from conventional diagnostic methods, this study improves the understanding of putative biomarkers at the molecular level which helps to unfold their role in normal immune function, fat synthesis/metabolism and oxidative stress pathways. Therefore, provides an opportunity to discover more accurate and sensitive diagnostic aids.

Broad de-regulated U2AF1 splicing is prognostic and augments leukemic transformation via protein arginine methyltransferase activation

The role of splicing dysregulation in cancer is underscored by splicing factor mutations; however, its impact in the absence of such rare mutations is poorly understood. To reveal complex patient subtypes and putative regulators of pathogenic splicing in Acute Myeloid Leukemia (AML), we developed a new approach called OncoSplice. Among diverse new subtypes, OncoSplice identified a biphasic poor prognosis signature that partially phenocopies U2AF1-mutant splicing, impacting thousands of genes in over 40% of adult and pediatric AML cases. U2AF1-like splicing co-opted a healthy circadian splicing program, was stable over time and induced a leukemia stem cell (LSC) program. Pharmacological inhibition of the implicated U2AF1-like splicing regulator, PRMT5, rescued leukemia mis-splicing and inhibited leukemic cell growth. Genetic deletion of IRAK4, a common target of U2AF1-like and PRMT5 treated cells, blocked leukemia development in xenograft models and induced differentiation. These analyses reveal a new prognostic alternative-splicing mechanism in malignancy, independent of splicing-factor mutations.

Epitranscriptomic Reader YTHDF2 Regulates SEK1( MAP2K4 )-JNK-cJUN Inflammatory Signaling in Astrocytes during Neurotoxic Stress

As the most abundant glial cells in the CNS, astrocytes dynamically respond to neurotoxic stress, however, the key molecular regulators controlling the inflammatory status of these sentinels during neurotoxic stress have remained elusive. Herein, we demonstrate that the m6A epitranscriptomic mRNA modification tightly regulates the pro-inflammatory functions of astrocytes. Specifically, the astrocytic neurotoxic stresser, manganese (Mn), downregulated the m6A reader YTHDF2 in human and mouse astrocyte cultures and in the mouse brain. Functionally, YTHDF2 knockdown augmented, while its overexpression dampened, neurotoxic stress induced proinflammatory response, suggesting YTHDF2 serves as a key upstream regulator of inflammatory responses in astrocytes. Mechnistically, YTHDF2 RIP-sequencing identified MAP2K4 ( MKK4; SEK1) mRNA as a YTHDF2 target influencing inflammatory signaling. Our target validation revealed Mn-exposed astrocytes mediates proinflammatory response by activating the phosphorylation of SEK1, JNK, and cJUN signaling. Collectively, YTHDF2 serves a key upstream 'molecular switch' controlling SEK1( MAP2K4 )-JNK-cJUN proinflammatory signaling in astrocytes.

Splicing neoantigen discovery with SNAF reveals shared targets for cancer immunotherapy

Immunotherapy has emerged as a crucial strategy to combat cancer by "reprogramming" a patient's own immune system. Although immunotherapy is typically reserved for patients with a high mutational burden, neoantigens produced from posttranscriptional regulation may provide an untapped reservoir of common immunogenic targets for new targeted therapies. To comprehensively define tumor-specific and likely immunogenic neoantigens from patient RNA-Seq, we developed Splicing Neo Antigen Finder (SNAF), an easy-to-use and open-source computational workflow to predict splicing-derived immunogenic MHC-bound peptides (T cell antigen) and unannotated transmembrane proteins with altered extracellular epitopes (B cell antigen). This workflow uses a highly accurate deep learning strategy for immunogenicity prediction (DeepImmuno) in conjunction with new algorithms to rank the tumor specificity of neoantigens (BayesTS) and to predict regulators of mis-splicing (RNA-SPRINT). T cell antigens from SNAF were frequently evidenced as HLA-presented peptides from mass spectrometry (MS) and predict response to immunotherapy in melanoma. Splicing neoantigen burden was attributed to coordinated splicing factor dysregulation. Shared splicing neoantigens were found in up to 90% of patients with melanoma, correlated to overall survival in multiple cancer cohorts, induced T cell reactivity, and were characterized by distinct cells of origin and amino acid preferences. In addition to T cell neoantigens, our B cell focused pipeline (SNAF-B) identified a new class of tumor-specific extracellular neoepitopes, which we termed ExNeoEpitopes. ExNeoEpitope full-length mRNA predictions were tumor specific and were validated using long-read isoform sequencing and in vitro transmembrane localization assays. Therefore, our systematic identification of splicing neoantigens revealed potential shared targets for therapy in heterogeneous cancers.

Population-level single-cell genomics reveals conserved gene programs in systemic juvenile idiopathic arthritis

Systemic autoimmune and autoinflammatory diseases are characterized by genetic and cellular heterogeneity. While current single-cell genomics methods provide insights into known disease subtypes, these analysis methods do not readily reveal novel cell-type perturbation programs shared among distinct patient subsets. Here, we performed single-cell RNA-Seq of PBMCs of patients with systemic juvenile idiopathic arthritis (SJIA) with diverse clinical manifestations, including macrophage activation syndrome (MAS) and lung disease (LD). We introduced two new computational frameworks called UDON and SATAY-UDON, which define patient subtypes based on their underlying disrupted cellular programs as well as associated biomarkers or clinical features. Among twelve independently identified subtypes, this analysis uncovered what we believe to be a novel complement and interferon activation program identified in SJIA-LD monocytes. Extending these analyses to adult and pediatric lupus patients found new but also shared disease programs with SJIA, including interferon and complement activation. Finally, supervised comparison of these programs in a compiled single-cell pan-immune atlas of over 1,000 healthy donors found a handful of normal healthy donors with evidence of early inflammatory activation in subsets of monocytes and platelets, nominating possible biomarkers for early disease detection. Thus, integrative pan-immune single-cell analysis resolved what we believe to be new conserved gene programs underlying inflammatory disease pathogenesis and associated complications.

Identification of Potential Genes and Critical Pathways in Postoperative Recurrence of Crohn's Disease by Machine Learning And WGCNA Network Analysis

Background

Crohn's disease (CD) is a chronic idiopathic inflammatory bowel disease affecting the entire gastrointestinal tract from the mouth to the anus. These patients often experience a period of symptomatic relapse and remission. A 20 - 30% symptomatic recurrence rate is reported in the first year after surgery, with a 10% increase each subsequent year. Thus, surgery is done only to relieve symptoms and not for the complete cure of the disease. The determinants and the genetic factors of this disease recurrence are also not well-defined. Therefore, enhanced diagnostic efficiency and prognostic outcome are critical for confronting CD recurrence.

Methods

We analysed ileal mucosa samples collected from neo-terminal ileum six months after surgery (M6=121 samples) from Crohn's disease dataset (GSE186582). The primary aim of this study is to identify the potential genes and critical pathways in post-operative recurrence of Crohn's disease. We combined the differential gene expression analysis with Recursive feature elimination (RFE), a machine learning approach to get five critical genes for the postoperative recurrence of Crohn's disease. The features (genes) selected by different methods were validated using five binary classifiers for recurrence and remission samples: Logistic Regression (LR), Decision tree classifier (DT), Support Vector Machine (SVM), Random Forest classifier (RF), and K-nearest neighbor (KNN) with 10-fold cross-validation. We also performed weighted gene co-expression network analysis (WGCNA) to select specific modules and feature genes associated with Crohn's disease postoperative recurrence, smoking, and biological sex. Combined with other biological interpretations, including Gene Ontology (GO) analysis, pathway enrichment, and protein-protein interaction (PPI) network analysis, our current study sheds light on the in-depth research of CD diagnosis and prognosis in postoperative recurrence.

Results

PLOD2, ZNF165, BOK, CX3CR1, and ARMCX4, are the important genes identified from the machine learning approach. These genes are reported to be involved in the viral protein interaction with cytokine and cytokine receptors, lysine degradation, and apoptosis. They are also linked with various cellular and molecular functions such as Peptidyl-lysine hydroxylation, Central nervous system maturation, G protein-coupled chemoattractant receptor activity, BCL-2 homology (BH) domain binding, Gliogenesis and negative regulation of mitochondrial depolarization. WGCNA identified a gene co-expression module that was primarily involved in mitochondrial translational elongation, mitochondrial translational termination, mitochondrial translation, mitochondrial respiratory chain complex, mRNA splicing via spliceosome pathways, etc.; Both the analysis result emphasizes that the mitochondrial depolarization pathway is linked with CD recurrence leading to oxidative stress in promoting inflammation in CD patients.

Conclusion

These key genes serve as the novel diagnostic biomarker for the postoperative recurrence of Crohn's disease. Thus, among other treatment options present until now, these biomarkers would provide success in both diagnosis and prognosis, aiming for a long-lasting remission to prevent further complications in CD.

Therapeutic targeting of the TPX2/TTK network in colorectal cancer

BACKGROUND

While the increased screening, changes in lifestyle, and recent advances in treatment regimen have decreased colorectal cancer (CRC) mortality, metastatic disease and recurrence remains a major clinical challenge. In the era of precision medicine, the identification of actionable novel therapeutic targets could ultimately offer an alternative treatment strategy for CRC.

METHODS

RNA-Seq was conducted using the illumina platform, while bioinformatics analyses were conducted using CLC genomics workbench and iDEP.951. Colony forming unit, flow cytometry, and fluorescent microscopy were used to assess cell proliferation, cell cycle distribution, and cell death, respectively. The growth potential of CRC cells under 3-dimensional (3D) conditions was assessed using Matrigel. STRING database (v11.5) and Ingenuity Pathway Analysis (IPA) tool were used for network and pathway analyses. CRISPR-Cas9 perturbational effects database was used to identify potential therapeutic targets for CRC, through integration with gene-drug interaction database. Structural modeling and molecular docking were used to assess the interaction between candidate drugs and their targets.

RESULTS

In the current study, we investigated the therapeutic potential of targeting TPX2, TTK, DDX39A, and LRP8, commonly upregulated genes in CRC identified through differential expression analysis in CRC and adjacent non-cancerous tissue. Targeted depletion of TPX2 and TTK impaired CRC proliferation, cell cycle progression, and organoid formation under 3D culture conditions, while suppression of DDX39A and LRP8 had modest effects on CRC colony formation. Differential expression analysis and bioinformatics on TPX2 and TTK-deficient cells identified cell cycle regulation as the hallmark associated with loss of TPX2 and TTK. Elevated expression of TPX2 and TTK correlated with an oncogenic state in tumor tissue from patients with colon adenocarcinoma, thus corroborating an oncogenic role for the TPX2/TTK network in the pathogenesis of CRC. Gene set enrichment and pathway analysis of TPX2high/TTKhigh CRC identified numerous additional gene targets as integral components of the TPX2/TTK network. Integration of TPX2/TTK enriched network with CRISPR-Cas9 functional screen data identified numerous novel dependencies for CRC. Additionally, gene-drug interaction analysis identified several druggable gene targets enriched in the TPX2/TTK network, including AURKA, TOP2A, CDK1, BIRC5, and many others.

CONCLUSIONS

Our data has implicated an essential role for TPX2 and TTK in CRC pathogenesis and identified numerous potential therapeutic targets and their drug interactions, suggesting their potential clinical use as a novel therapeutic strategy for patients with CRC. Video Abstract.

Advances in alternative splicing identification: deep learning and pantranscriptome

In plants, alternative splicing is a crucial mechanism for regulating gene expression at the post-transcriptional level, which leads to diverse proteins by generating multiple mature mRNA isoforms and diversify the gene regulation. Due to the complexity and variability of this process, accurate identification of splicing events is a vital step in studying alternative splicing. This article presents the application of alternative splicing algorithms with or without reference genomes in plants, as well as the integration of advanced deep learning techniques for improved detection accuracy. In addition, we also discuss alternative splicing studies in the pan-genomic background and the usefulness of integrated strategies for fully profiling alternative splicing.

Repression of TRIM13 by chromatin assembly factor CHAF1B is critical for AML development

Acute myeloid leukemia (AML) is an aggressive blood cancer that stems from the rapid expansion of immature leukemic blasts in the bone marrow. Mutations in epigenetic factors represent the largest category of genetic drivers of AML. The chromatin assembly factor CHAF1B is a master epigenetic regulator of transcription associated with self-renewal and the undifferentiated state of AML blasts. Upregulation of CHAF1B, as observed in almost all AML samples, promotes leukemic progression by repressing the transcription of differentiation factors and tumor suppressors. However, the specific factors regulated by CHAF1B and their contributions to leukemogenesis are unstudied. We analyzed RNA sequencing data from mouse MLL-AF9 leukemic cells and bone marrow aspirates, representing a diverse collection of pediatric AML samples and identified the E3 ubiquitin ligase TRIM13 as a target of CHAF1B-mediated transcriptional repression associated with leukemogenesis. We found that CHAF1B binds the promoter of TRIM13, resulting in its transcriptional repression. In turn, TRIM13 suppresses self-renewal of leukemic cells by promoting pernicious entry into the cell cycle through its nuclear localization and catalytic ubiquitination of cell cycle-promoting protein, CCNA1. Overexpression of TRIM13 initially prompted a proliferative burst in AML cells, which was followed by exhaustion, whereas loss of total TRIM13 or deletion of its catalytic domain enhanced leukemogenesis in AML cell lines and patient-derived xenografts. These data suggest that CHAF1B promotes leukemic development, in part, by repressing TRIM13 expression and that this relationship is necessary for leukemic progression.

Long non-coding RNA AC099850.4 correlates with advanced disease state and predicts worse prognosis in triple-negative breast cancer

Our understanding of the function of long non-coding RNAs (lncRNAs) in health and disease states has evolved over the past decades due to the many advances in genome research. In the current study, we characterized the lncRNA transcriptome enriched in triple-negative breast cancer (TNBC, n = 42) and estrogen receptor (ER+, n = 42) breast cancer compared to normal breast tissue (n = 56). Given the aggressive nature of TNBC, our data revealed selective enrichment of 57 lncRNAs in TNBC. Among those, AC099850.4 lncRNA was chosen for further investigation where it exhibited elevated expression, which was further confirmed in a second TNBC cohort (n = 360) where its expression correlated with a worse prognosis. Network analysis of AC099850.4high TNBC highlighted enrichment in functional categories indicative of cell cycle activation and mitosis. Ingenuity pathway analysis on the differentially expressed genes in AC099850.4high TNBC revealed the activation of the canonical kinetochore metaphase signaling pathway, pyridoxal 5'-phosphate salvage pathway, and salvage pathways of pyrimidine ribonucleotides. Additionally, upstream regulator analysis predicted the activation of several upstream regulator networks including CKAP2L, FOXM1, RABL6, PCLAF, and MITF, while upstream regulator networks of TP53, NUPR1, TRPS1, and CDKN1A were suppressed. Interestingly, elevated expression of AC099850.4 correlated with worse short-term relapse-free survival (log-rank p = 0.01). Taken together, our data are the first to reveal AC099850.4 as an unfavorable prognostic marker in TNBC, associated with more aggressive clinicopathological features, and suggest its potential utilization as a prognostic biomarker and therapeutic target in TNBC.

Uncovering Novel Roles of miR-122 in the Pathophysiology of the Liver: Potential Interaction with NRF1 and E2F4 Signaling

MicroRNA miR-122 plays a pivotal role in liver function. Despite numerous studies investigating this miRNA, the global network of genes regulated by miR-122 and its contribution to the underlying pathophysiological mechanisms remain largely unknown. To gain a deeper understanding of miR-122 activity, we employed two complementary approaches. Firstly, through transcriptome analysis of polyribosome-bound RNAs, we discovered that miR-122 exhibits potential antagonistic effects on specific transcription factors known to be dysregulated in liver disease, including nuclear respiratory factor-1 (NRF1) and the E2F transcription factor 4 (E2F4). Secondly, through proteome analysis of hepatoma cells transfected with either miR-122 mimic or antagomir, we discovered changes in several proteins associated with increased malignancy. Interestingly, many of these proteins were reported to be transcriptionally regulated by NRF1 and E2F4, six of which we validated as miR-122 targets. Among these, a negative correlation was observed between miR-122 and glucose-6-phosphate dehydrogenase levels in the livers of patients with hepatitis B virus-associated hepatocellular carcinoma. This study provides novel insights into potential alterations of molecular pathway occurring at the early stages of liver disease, driven by the dysregulation of miR-122 and its associated genes.

Differential Modulation of miR-122 Transcription by TGFβ1/BMP6: Implications for Nonresolving Inflammation and Hepatocarcinogenesis

Chronic inflammation is widely recognized as a significant factor that promotes and worsens the development of malignancies, including hepatocellular carcinoma. This study aimed to explore the potential role of microRNAs in inflammation-associated nonresolving hepatocarcinogenesis. By conducting a comprehensive analysis of altered microRNAs in animal models with liver cancer of various etiologies, we identified miR-122 as the most significantly downregulated microRNA in the liver of animals with inflammation-associated liver cancer. Although previous research has indicated the importance of miR-122 in maintaining hepatocyte function, its specific role as either the trigger or the consequence of underlying diseases remains unclear. Through extensive analysis of animals and in vitro models, we have successfully demonstrated that miR-122 transcription is differentially regulated by the immunoregulatory cytokines, by the transforming growth factor-beta 1 (TGFβ1), and the bone morphogenetic protein-6 (BMP6). Furthermore, we presented convincing evidence directly linking reduced miR-122 transcription to inflammation and in chronic liver diseases. The results of this study strongly suggest that prolonged activation of pro-inflammatory signaling pathways, leading to disruption of cytokine-mediated regulation of miR-122, may significantly contribute to the onset and exacerbation of chronic liver disease.

Unraveling the genetic mechanisms governing the host response to bovine anaplasmosis

Bovine anaplasmosis caused by Anaplasma marginale is a tick-borne disease of livestock with widespread prevalence and huge economic implications. In order to get new insights into modulation of host gene expression in response to natural infections of anaplasmosis, this study is the first attempt that compared the transcriptome profiles of peripheral blood mononuclear cells (PBMCs) of A. marginale infected and healthy crossbred cattle. Transcriptome analysis identified shared as well as unique functional pathways in the two groups. Translation and structural constituent of ribosome were the important terms for the genes abundantly expressed in the infected as well as healthy animals. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of the differentially expressed genes revealed that immunity and signal transduction related terms were enriched for the up-regulated genes in the infected animals. The over-represented pathways were cytokine-cytokine receptor interaction and signaling pathways involving chemokines, Interleukin 17 (IL17), Tumour Necrosis Factor (TNF), Nuclear Factor Kappa B (NFKB) etc. Interestingly, many genes previously associated with parasite-borne diseases such as amoebiasis, trypanosomiasis, toxoplasmosis, and leishmaniasis were profusely expressed in the dataset of the diseased animals. High expression was also evident for the genes for acute phase response proteins, anti-microbial peptides and many inflammatory cytokines. Role of cytokines in mediating communication between immune cells was the most conspicuous gene network identified through the Ingenuity Pathway Analysis. This study provides comprehensive information about the crosstalk of genes involved in host defense as well as parasite persistence in the host upon infection with A. marginale.

BAFF and APRIL counterregulate susceptibility to inflammation-induced preterm birth

Clinical evidence points to a function for B cell-activating factor (BAFF) in pregnancy. However, direct roles for BAFF-axis members in pregnancy have not been examined. Here, via utility of genetically modified mice, we report that BAFF promotes inflammatory responsiveness and increases susceptibility to inflammation-induced preterm birth (PTB). In contrast, we show that the closely related A proliferation-inducing ligand (APRIL) decreases inflammatory responsiveness and susceptibility to PTB. Known BAFF-axis receptors serve a redundant function in signaling BAFF/APRIL presence in pregnancy. Treatment with anti-BAFF/APRIL monoclonal antibodies or BAFF/APRIL recombinant proteins is sufficient to manipulate susceptibility to PTB. Notably, macrophages at the maternal-fetal interface produce BAFF, while BAFF and APRIL presence divergently shape macrophage gene expression and inflammatory function. Overall, our findings demonstrate that BAFF and APRIL play divergent inflammatory roles in pregnancy and provide therapeutic targets for mitigating risk of inflammation-induced PTB.

Lung Epithelial Cells from Obese Patients Have Impaired Control of SARS-CoV-2 Infection

Obesity is known to increase the complications of the COVID-19 coronavirus disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). However, the exact mechanisms of SARS-CoV-2 infection in obese patients have not been clearly elucidated. This study aims to better understand the effect of obesity on the course of SARS-CoV-2 infection and identify candidate molecular pathways involved in the progression of the disease, using an in vitro live infection model and RNA sequencing. Results from this study revealed the enhancement of viral load and replication in bronchial epithelial cells (NHBE) from obese subjects at 24 h of infection (MOI = 0.5) as compared to non-obese subjects. Transcriptomic profiling via RNA-Seq highlighted the enrichment of lipid metabolism-related pathways along with LPIN2, an inflammasome regulator, as a unique differentially expressed gene (DEG) in infected bronchial epithelial cells from obese subjects. Such findings correlated with altered cytokine and angiotensin-converting enzyme-2 (ACE2) expression during infection of bronchial cells. These findings provide a novel insight on the molecular interplay between obesity and SARS-CoV-2 infection. In conclusion, this study demonstrates the increased SARS-CoV-2 infection of bronchial epithelial cells from obese subjects and highlights the impaired immunity which may explain the increased severity among obese COVID-19 patients.

Quantifying tumor specificity using Bayesian probabilistic modeling for drug target discovery and prioritization

In diseases such as cancer, the design of new therapeutic strategies requires extensive, costly, and unfortunately sometimes deadly testing to reveal life threatening "off target" effects. A crucial first step in predicting toxicity are analyses of normal RNA and protein tissue expression, which are now possible using comprehensive molecular tissue atlases. However, no standardized approaches exist for target prioritization, which instead rely on ad-hoc thresholds and manual inspection. Such issues are compounded, given that genomic and proteomic data detection sensitivity and accuracy are often problematic. Thus, quantifiable probabilistic scores for tumor specificity that address these challenges could enable the creation of new predictive models for combinatorial drug design and correlative analyses. Here, we propose a Bayesian Tumor Specificity (BayesTS) score that can naturally account for multiple independent forms of molecular evidence derving from both RNA-Seq and protein expression while preserving the uncertainty of the inference. We applied BayesTS to 24,905 human protein-coding genes across 3,644 normal samples (GTEx and TCGA) spanning 63 tissues. These analyses demonstrate the ability of BayesTS to accurately incorporate protein, RNA and tissue distribution evidence, while effectively capturing the uncertainty of these inferences. This approach prioritized well-established drug targets, while deemphasizing those which were later found to induce toxicity. BayesTS allows for the adjustment of tissue importance weights for tissues of interest, such as reproductive and physiologically dispensable tissues (e.g., tonsil, appendix), enabling clinically translatable prioritizations. Our results show that BayesTS can facilitate novel drug target discovery and can be easily generalized to unconventional molecular targets, such as splicing neoantigens. We provide the code and inferred tumor specificity predictions as a database available online (https://github.com/frankligy/BayesTS). We envision that the widespread adoption of BayesTS will facilitate improved target prioritization for oncology drug development, ultimately leading to the discovery of more effective and safer drugs.

The early neutrophil-committed progenitors aberrantly differentiate into immunoregulatory monocytes during emergency myelopoiesis

Inflammatory stimuli cause a state of emergency myelopoiesis leading to neutrophil-like monocyte expansion. However, their function, the committed precursors, or growth factors remain elusive. In this study we find that Ym1⁺Ly6C^hi monocytes, an immunoregulatory entity of neutrophil-like monocytes, arise from progenitors of neutrophil 1 (proNeu1). Granulocyte-colony stimulating factor (G-CSF) favors the production of neutrophil-like monocytes through previously unknown CD81⁺CX3CR1^lo monocyte precursors. GFI1 promotes the differentiation of proNeu2 from proNeu1 at the cost of producing neutrophil-like monocytes. The human counterpart of neutrophil-like monocytes that also expands in response to G-CSF is found in CD14⁺CD16^- monocyte fraction. The human neutrophil-like monocytes are discriminated from CD14⁺CD16^- classical monocytes by CXCR1 expression and the capacity to suppress T cell proliferation. Collectively, our findings suggest that the aberrant expansion of neutrophil-like monocytes under inflammatory conditions is a process conserved between mouse and human, which may be beneficial for the resolution of inflammation.

Improving anti-tumor efficacy of low-dose Vincristine in rhabdomyosarcoma via the combination therapy with FOXM1 inhibitor RCM1

Rhabdomyosarcoma (RMS) is a highly metastatic soft-tissue sarcoma that often develops resistance to current therapies, including vincristine. Since the existing treatments have not significantly improved survival, there is a critical need for new therapeutic approaches for RMS patients. FOXM1, a known oncogene, is highly expressed in RMS, and is associated with the worst prognosis in RMS patients. In the present study, we found that the combination treatment with specific FOXM1 inhibitor RCM1 and low doses of vincristine is more effective in increasing apoptosis and decreasing RMS cell proliferation in vitro compared to single drugs alone. Since RCM1 is highly hydrophobic, we developed innovative nanoparticle delivery system containing poly-beta-amino-esters and folic acid (NP^FA), which efficiently delivers RCM1 to mouse RMS tumors in vivo. The combination of low doses of vincristine together with intravenous administration of NP^FA nanoparticles containing RCM1 effectively reduced RMS tumor volumes, increased tumor cell death and decreased tumor cell proliferation in RMS tumors compared to RCM1 or vincristine alone. The combination therapy was non-toxic as demonstrated by liver metabolic panels using peripheral blood serum. Using RNA-seq of dissected RMS tumors, we identified Chac1 as a uniquely downregulated gene after the combination treatment. Knockdown of Chac1 in RMS cells in vitro recapitulated the effects of the combination therapy. Altogether, combination treatment with low doses of vincristine and nanoparticle delivery of FOXM1 inhibitor RCM1 in a pre-clinical model of RMS has superior anti-tumor effects and decreases CHAC1 while reducing vincristine toxicity.

Tip60/KAT5 Histone Acetyltransferase Is Required for Maintenance and Neurogenesis of Embryonic Neural Stem Cells

Epigenetic regulation via epigenetic factors in collaboration with tissue-specific transcription factors is curtail for establishing functional organ systems during development. Brain development is tightly regulated by epigenetic factors, which are coordinately activated or inactivated during processes, and their dysregulation is linked to brain abnormalities and intellectual disability. However, the precise mechanism of epigenetic regulation in brain development and neurogenesis remains largely unknown. Here, we show that Tip60/KAT5 deletion in neural stem/progenitor cells (NSCs) in mice results in multiple abnormalities of brain development. Tip60-deficient embryonic brain led to microcephaly, and proliferating cells in the developing brain were reduced by Tip60 deficiency. In addition, neural differentiation and neuronal migration were severely affected in Tip60-deficient brains. Following neurogenesis in developing brains, gliogenesis started from the earlier stage of development in Tip60-deficient brains, indicating that Tip60 is involved in switching from neurogenesis to gliogenesis during brain development. It was also confirmed in vitro that poor neurosphere formation, proliferation defects, neural differentiation defects, and accelerated astrocytic differentiation in mutant NSCs are derived from Tip60-deficient embryonic brains. This study uncovers the critical role of Tip60 in brain development and NSC maintenance and function in vivo and in vitro.

miR-196b-TLR7/8 Signaling Axis Regulates Innate Immune Signaling and Myeloid Maturation in DNMT3A-Mutant AML

PURPOSE

DNMT3A mutations confer a poor prognosis in acute myeloid leukemia (AML), but the molecular mechanisms downstream of DNMT3A mutations in disease pathogenesis are not completely understood, limiting targeted therapeutic options. The role of miRNA in DNMT3A-mutant AML pathogenesis is understudied.

EXPERIMENTAL DESIGN

DNA methylation and miRNA expression was evaluated in human AML patient samples and in Dnmt3a/Flt3-mutant AML mice. The treatment efficacy and molecular mechanisms of TLR7/8-directed therapies on DNMT3A-mutant AML were evaluated in vitro on human AML patient samples and in Dnmt3a/Flt3-mutant AML mice.

RESULTS

miR-196b is hypomethylated and overexpressed in DNMT3A-mutant AML and is associated with poor patient outcome. miR-196b overexpression in DNMT3A-mutant AML is important to maintain an immature state and leukemic cell survival through repression of TLR signaling. The TLR7/8 agonist resiquimod induces dendritic cell-like differentiation with costimulatory molecule expression in DNMT3A-mutant AML cells and provides a survival benefit to Dnmt3a/Flt3-mutant AML mice. The small molecule bryostatin-1 augments resiquimod-mediated AML growth inhibition and differentiation.

CONCLUSIONS

DNMT3A loss-of-function mutations cause miRNA locus-specific hypomethylation and overexpression important for mutant DNMT3A-mediated pathogenesis and clinical outcomes. Specifically, the overexpression of miR-196b in DNMT3A-mutant AML creates a novel therapeutic vulnerability by controlling sensitivity to TLR7/8-directed therapies.

Fetal maturation revealed by amniotic fluid cell-free transcriptome in rhesus macaques

Accurate estimate of fetal maturity could provide individualized guidance for delivery of complicated pregnancies. However, current methods are invasive, have low accuracy, and are limited to fetal lung maturation. To identify diagnostic gestational biomarkers, we performed transcriptomic profiling of lung and brain, as well as cell-free RNA from amniotic fluid of preterm and term rhesus macaque fetuses. These data identify potentially new and prior-associated gestational age differences in distinct lung and neuronal cell populations when compared with existing single-cell and bulk RNA-Seq data. Comparative analyses found hundreds of genes coincidently induced in lung and amniotic fluid, along with dozens in brain and amniotic fluid. These data enable creation of computational models that accurately predict lung compliance from amniotic fluid and lung transcriptome of preterm fetuses treated with antenatal corticosteroids. Importantly, antenatal steroids induced off-target gene expression changes in the brain, impinging upon synaptic transmission and neuronal and glial maturation, as this could have long-term consequences on brain development. Cell-free RNA in amniotic fluid may provide a substrate of global fetal maturation markers for personalized management of at-risk pregnancies.

Prediction of potential molecular markers of bovine mastitis by meta-analysis of differentially expressed genes using combined p value and robust rank aggregation

Bovine mastitis causes significant economic loss to the dairy industry by affecting milk quality and quantity. Escherichia coli and Staphylococcus aureus are the two common mastitis-causing bacteria among the consortia of mastitis pathogens, wherein E. coli is an opportunistic environmental pathogen, and S. aureus is a contagious pathogen. This study was designed to predict molecular markers of bovine mastitis by meta-analysis of differentially expressed genes (DEG) in E. coli- or S. aureus-infected mammary epithelial cells (MECs) using p value combination and robust rank aggregation (RRA) methods. High-throughput transcriptome of bovine MECs, infected with E. coli or S. aureus, were analyzed, and correlation of z-scores were computed for the expression datasets to identify the lineage profile and functional ontology of DEGs. Key pathways enriched in infected MECs were deciphered by Gene Set Enrichment Analysis (GSEA), following which combined p value and RRA were used to perform DEG meta-analysis to limit type I error in the analysis. The miRNA-gene networks were then built to uncover potential molecular markers of mastitis. Lineage profiling of MECs showed that the gene expression levels were associated with mammary tissue lineage. The up-regulated genes were enriched in immune-related pathways, whereas down-regulated genes influenced the cellular processes. GSEA analysis of DEGs deciphered the involvement of Toll-like receptor (TLR), and NF-kappa B signaling pathway during infection. Comparison after meta-analysis yielded with genes ZC3H12A, RND1, and MAP3K8 having significant expression levels in both E. coli and S. aureus dataset, and on evaluating miRNA-gene network, 7 pairs were common to both sets identifying them as potential molecular markers.

First report on delineation of differentially expressed genes and pathways in milk somatic cells of mastitic and healthy Murrah buffaloes

Despite immense contribution of buffaloes as dairy species, limited studies have addressed the bubaline mastitis as compared to cattle. This was the first differential transcriptomic study investigating the alterations induced by clinical mastitis in buffalo milk relative to healthy controls. Comparative gene expression profiling of three biological replicates of each group identified 1014 up-regulated and 999 down-regulated genes in the diseased buffaloes (Fold change > 2, FDR < 0.05). Activation of immune and inflammatory responses were the most enriched GO terms in the mastitic animals, with higher transcript abundance of many genes coding for anti-microbial proteins such as β-defensins, perforin, granzymes, granulysin, cathelicidins etc. Analysis of the gene regulatory interactions of the up-regulated DEGs identified many hub genes that govern the cellular and macromolecular metabolic processes (E2F4, E2F1, RBL2, FOXM1, IRF1 and MYB). This study contributes to an insightful understanding of molecular mechanisms governing immune response of buffaloes to mastitis.

Recursive Feature Elimination-based Biomarker Identification for Open Neural Tube Defects

Background: Open spina bifida (myelomeningocele) is the result of the failure of spinal cord closing completely and is the second most common and severe birth defect. Open neural tube defects are multifactorial, and the exact molecular mechanism of the pathogenesis is not clear due to disease complexity for which prenatal treatment options remain limited worldwide. Artificial intelligence techniques like machine learning tools have been increasingly used in precision diagnosis. Objective: The primary objective of this study is to identify key genes for open neural tube defects using a machine learning approach that provides additional information about myelomeningocele in order to obtain a more accurate diagnosis. Materials and Methods: Our study reports differential gene expression analysis from multiple datasets (GSE4182 and GSE101141) of amniotic fluid samples with open neural tube defects. The sample outliers in the datasets were detected using principal component analysis (PCA). We report a combination of the differential gene expression analysis with recursive feature elimination (RFE), a machine learning approach to get 4 key genes for open neural tube defects. The features selected were validated using five binary classifiers for diseased and healthy samples: Logistic Regression (LR), Decision tree classifier (DT), Support Vector Machine (SVM), Random Forest classifier (RF), and K-nearest neighbour (KNN) with 5-fold cross-validation. Results: Growth Associated Protein 43 (GAP43), Glial fibrillary acidic protein (GFAP), Repetin (RPTN), and CD44 are the important genes identified in the study. These genes are known to be involved in axon growth, astrocyte differentiation in the central nervous system, post-traumatic brain repair, neuroinflammation, and inflammation-linked neuronal injuries. These key genes represent a promising tool for further studies in the diagnosis and early detection of open neural tube defects. Conclusion: These key biomarkers help in the diagnosis and early detection of open neural tube defects, thus evaluating the progress and seriousness in diseases condition. This study strengthens previous literature sources of confirming these biomarkers linked with open NTD's. Thus, among other prenatal treatment options present until now, these biomarkers help in the early detection of open neural tube defects, which provides success in both treatment and prevention of these defects in the advanced stage.

Premature Neural Progenitor Cell Differentiation Into Astrocytes in Retinoic Acid-Induced Spina Bifida Rat Model

During embryonic spinal cord development, neural progenitor cells (NPCs) generate three major cell lines: neurons, oligodendrocytes, and astrocytes at precise times and locations within the spinal cord. Recent studies demonstrate early astrogenesis in animal models of spina bifida, which may play a role in neuronal dysfunction associated with this condition. However, to date, the pathophysiological mechanisms related to this early astrocytic response in spina bifida are poorly understood. This study aimed to characterize the development of early astrogliosis over time from Pax6+, Olig2+, or Nkx2.2+ NPCs using a retinoic acid-induced spina bifida rat model. At three gestational ages (E15, E17, and E20), spinal cords from fetuses with retinoic acid-induced spina bifida, their healthy sibling controls, or fetuses treated with the vehicle control were analyzed. Results indicated that premature astrogliosis and astrocytic activation were associated with an altered presence of Pax6+, Olig2+, and Nkx2.2+ NPCs in the lesion compared to the controls. Finally, this response correlated with an elevation in genes involved in the Notch-BMP signaling pathway. Taken together, changes in NPC patterning factor expression with Notch-BMP signaling upregulation may be responsible for the altered astrogenesis patterns observed in the spinal cord in a retinoic acid-induced spina bifida model.

ASpediaFI: Functional Interaction Analysis of Alternative Splicing Events

Alternative splicing (AS) regulates biological processes governing phenotypes and diseases. Differential AS (DAS) gene test methods have been developed to investigate important exonic expression from high-throughput datasets. However, the DAS events extracted using statistical tests are insufficient to delineate relevant biological processes. In this study, we developed a novel application, Alternative Splicing Encyclopedia: Functional Interaction (ASpediaFI), to systemically identify DAS events and co-regulated genes and pathways. ASpediaFI establishes a heterogeneous interaction network of genes and their feature nodes (i.e., AS events and pathways) connected by co-expression or pathway gene set knowledge. Next, ASpediaFI explores the interaction network using the random walk with restart algorithm and interrogates the proximity from a query gene set. Finally, ASpediaFI extracts significant AS events, genes, and pathways. To evaluate the performance of our method, we simulated RNA sequencing (RNA-seq) datasets to consider various conditions of sequencing depth and sample size. The performance was compared with that of other methods. Additionally, we analyzed three public datasets of cancer patients or cell lines to evaluate how well ASpediaFI detects biologically relevant candidates. ASpediaFI exhibits strong performance in both simulated and public datasets. Our integrative approach reveals that DAS events that recognize a global co-expression network and relevant pathways determine the functional importance of spliced genes in the subnetwork. ASpediaFI is publicly available at https://bioconductor.org/packages/ASpediaFI.

Nuclear Vav3 is required for polycomb repression complex-1 activity in B-cell lymphoblastic leukemogenesis

Acute B-cell lymphoblastic leukemia (B-ALL) results from oligo-clonal evolution of B-cell progenitors endowed with initiating and propagating leukemia properties. The activation of both the Rac guanine nucleotide exchange factor (Rac GEF) Vav3 and Rac GTPases is required for leukemogenesis mediated by the oncogenic fusion protein BCR-ABL. Vav3 expression becomes predominantly nuclear upon expression of BCR-ABL signature. In the nucleus, Vav3 interacts with BCR-ABL, Rac, and the polycomb repression complex (PRC) proteins Bmi1, Ring1b and Ezh2. The GEF activity of Vav3 is required for the proliferation, Bmi1-dependent B-cell progenitor self-renewal, nuclear Rac activation, protein interaction with Bmi1, mono-ubiquitination of H2A(K119) (H2AK119Ub) and repression of PRC-1 (PRC1) downstream target loci, of leukemic B-cell progenitors. Vav3 deficiency results in de-repression of negative regulators of cell proliferation and repression of oncogenic transcriptional factors. Mechanistically, we show that Vav3 prevents the Phlpp2-sensitive and Akt (S473)-dependent phosphorylation of Bmi1 on the regulatory residue S314 that, in turn, promotes the transcriptional factor reprogramming of leukemic B-cell progenitors. These results highlight the importance of non-canonical nuclear Rho GTPase signaling in leukemogenesis.

Ph⁺ and Ph-like B-ALL remain poor prognosis leukemias. VAV3, a guanine nucleotide exchange factor, is activated and overexpressed in these leukemias. Here the authors reveal that leukemic VAV3 is predominantly nuclear. Nuclear VAV3, through its guanine nucleotide exchange factor and its effector nuclear RAC2, controls the repressive transcriptional activity of the polycomb repression complex-1 via nuclear AKT/PHLPP2 regulated BMI1.

BIN1 is a key regulator of proinflammatory and neurodegeneration-related activation in microglia

BACKGROUND

The BIN1 locus contains the second-most significant genetic risk factor for late-onset Alzheimer's disease. BIN1 undergoes alternate splicing to generate tissue- and cell-type-specific BIN1 isoforms, which regulate membrane dynamics in a range of crucial cellular processes. Whilst the expression of BIN1 in the brain has been characterized in neurons and oligodendrocytes in detail, information regarding microglial BIN1 expression is mainly limited to large-scale transcriptomic and proteomic data. Notably, BIN1 protein expression and its functional roles in microglia, a cell type most relevant to Alzheimer's disease, have not been examined in depth.

METHODS

Microglial BIN1 expression was analyzed by immunostaining mouse and human brain, as well as by immunoblot and RT-PCR assays of isolated microglia or human iPSC-derived microglial cells. Bin1 expression was ablated by siRNA knockdown in primary microglial cultures in vitro and Cre-lox mediated conditional deletion in adult mouse brain microglia in vivo. Regulation of neuroinflammatory microglial signatures by BIN1 in vitro and in vivo was characterized using NanoString gene panels and flow cytometry methods. The transcriptome data was explored by in silico pathway analysis and validated by complementary molecular approaches.

RESULTS

Here, we characterized microglial BIN1 expression in vitro and in vivo and ascertained microglia expressed BIN1 isoforms. By silencing Bin1 expression in primary microglial cultures, we demonstrate that BIN1 regulates the activation of proinflammatory and disease-associated responses in microglia as measured by gene expression and cytokine production. Our transcriptomic profiling revealed key homeostatic and lipopolysaccharide (LPS)-induced inflammatory response pathways, as well as transcription factors PU.1 and IRF1 that are regulated by BIN1. Microglia-specific Bin1 conditional knockout in vivo revealed novel roles of BIN1 in regulating the expression of disease-associated genes while counteracting CX3CR1 signaling. The consensus from in vitro and in vivo findings showed that loss of Bin1 impaired the ability of microglia to mount type 1 interferon responses to proinflammatory challenge, particularly the upregulation of a critical type 1 immune response gene, Ifitm3.

CONCLUSIONS

Our convergent findings provide novel insights into microglial BIN1 function and demonstrate an essential role of microglial BIN1 in regulating brain inflammatory response and microglial phenotypic changes. Moreover, for the first time, our study shows a regulatory relationship between Bin1 and Ifitm3, two Alzheimer's disease-related genes in microglia. The requirement for BIN1 to regulate Ifitm3 upregulation during inflammation has important implications for inflammatory responses during the pathogenesis and progression of many neurodegenerative diseases.

The Interplay of Microtubules with Mitochondria–ER Contact Sites (MERCs) in Glioblastoma

Gliomas are heterogeneous neoplasms, classified into grade I to IV according to their malignancy and the presence of specific histological/molecular hallmarks. The higher grade of glioma is known as glioblastoma (GB). Although progress has been made in surgical and radiation treatments, its clinical outcome is still unfavorable. The invasive properties of GB cells and glioma aggressiveness are linked to the reshaping of the cytoskeleton. Recent works suggest that the different susceptibility of GB cells to antitumor immune response is also associated with the extent and function of mitochondria–ER contact sites (MERCs). The presence of MERCs alterations could also explain the mitochondrial defects observed in GB models, including abnormalities of energy metabolism and disruption of apoptotic and calcium signaling. Based on this evidence, the question arises as to whether a MERCs–cytoskeleton crosstalk exists, and whether GB progression is linked to an altered cytoskeleton–MERCs interaction. To address this possibility, in this review we performed a meta-analysis to compare grade I and grade IV GB patients. From this preliminary analysis, we found that GB samples (grade IV) are characterized by altered expression of cytoskeletal and MERCs related genes. Among them, the cytoskeleton-associated protein 4 (CKAP4 or CLIMP-63) appears particularly interesting as it encodes a MERCs protein controlling the ER anchoring to microtubules (MTs). Although further in-depth analyses remain necessary, this perspective review may provide new hints to better understand GB molecular etiopathogenesis, by suggesting that cytoskeletal and MERCs alterations cooperate to exacerbate the cellular phenotype of high-grade GB and that MERCs players can be exploited as novel biomarkers/targets to enhance the current therapy for GB.

The Tumor Microenvironment Drives Intrahepatic Cholangiocarcinoma Progression

Intrahepatic cholangiocarcinoma (iCCA) is a highly aggressive cancer with limited therapeutic options and short overall survival. iCCA is characterized by a strong desmoplastic reaction in the surrounding ecosystem that likely affects tumoral progression. Overexpression of the Notch pathway is implicated in iCCA development and progression. Our aim was to investigate the effectiveness of Crenigacestat, a selective inhibitor of NOTCH1 signaling, against the cross-talk between cancer cells and the surrounding ecosystem in an in vivo HuCCT1-xenograft model. In the present study, a transcriptomic analysis approach, validated by Western blotting and qRT-PCR on iCCA tumor masses treated with Crenigacestat, was used to study the molecular pathways responsive to drug treatment. Our results indicate that Crenigacestat significantly inhibited NOTCH1 and HES1, whereas tumor progression was not affected. In addition, the drug triggered a strong immune response and blocked neovascularization in the tumor ecosystem of the HuCCT1-xenograft model without affecting the occurrence of fibrotic reactions. Therefore, although these data need further investigation, our observations confirm that Crenigacestat selectively targets NOTCH1 and that the desmoplastic response in iCCA likely plays a key role in both drug effectiveness and tumor progression.

A Pathogenic Th17/CD38+ Macrophage Feedback Loop Drives Inflammatory Arthritis through TNF-α

The pathobiology of rheumatoid inflammatory diseases, including rheumatoid arthritis (RA) and psoriatic arthritis, involves the interplay between innate and adaptive immune components and resident synoviocytes. Single-cell analyses of patient samples and relevant mouse models have characterized many cellular subsets in RA. However, the impact of interactions between cell types is not fully understood. In this study, we temporally profiled murine arthritic synovial isolates at the single-cell level to identify perturbations similar to those found in human RA. Notably, murine macrophage subtypes like those found in RA patients were expanded in arthritis and linked to promoting the function of Th17 cells in the joint. In vitro experiments identified a capacity for murine macrophages to maintain the functionality and expansion of Th17 cells. Reciprocally, murine Th17 cell-derived TNF-α induced CD38⁺ macrophages that enhanced Th17 functionality. Murine synovial CD38⁺ macrophages were expanded during arthritis, and their depletion or blockade via TNF-α neutralization alleviated disease while reducing IL-17A-producing cells. These findings identify a cellular feedback loop that promotes Th17 cell pathogenicity through TNF-α to drive inflammatory arthritis.

DJExpress: An Integrated Application for Differential Splicing Analysis and Visualization

RNA-seq analysis of alternative pre-mRNA splicing has facilitated an unprecedented understanding of transcriptome complexity in health and disease. However, despite the availability of countless bioinformatic pipelines for transcriptome-wide splicing analysis, the use of these tools is often limited to expert bioinformaticians. The need for high computational power, combined with computational outputs that are complicated to visualize and interpret present obstacles to the broader research community. Here we introduce DJExpress, an R package for differential expression analysis of transcriptomic features and expression-trait associations. To determine gene-level differential junction usage as well as associations between junction expression and molecular/clinical features, DJExpress uses raw splice junction counts as input data. Importantly, DJExpress runs on an average laptop computer and provides a set of interactive and intuitive visualization formats. In contrast to most existing pipelines, DJExpress can handle both annotated and de novo identified splice junctions, thereby allowing the quantification of novel splice events. Moreover, DJExpress offers a web-compatible graphical interface allowing the analysis of user-provided data as well as the visualization of splice events within our custom database of differential junction expression in cancer (DJEC DB). DJEC DB includes not only healthy and tumor tissue junction expression data from TCGA and GTEx repositories but also cancer cell line data from the DepMap project. The integration of DepMap functional genomics data sets allows association of junction expression with molecular features such as gene dependencies and drug response profiles. This facilitates identification of cancer cell models for specific splicing alterations that can then be used for functional characterization in the lab. Thus, DJExpress represents a powerful and user-friendly tool for exploration of alternative splicing alterations in RNA-seq data, including multi-level data integration of alternative splicing signatures in healthy tissue, tumors and cancer cell lines.

Time Course Transcriptome Analysis of Spina Bifida Progression in Fetal Rats

A better understanding of the transcriptomic modifications that occur in spina bifida may lead to identify mechanisms involved in the progression of spina bifida in utero and the development of new therapeutic strategies that aid in spinal cord regeneration after surgical interventions. In this study, RNA-sequencing was used to identify differentially expressed genes in fetal spinal cords from rats with retinoic acid-induced spina bifida at E15, E17, and E20. Gene ontology, KEGG, and protein–protein interaction analysis were conducted to predict pathways involved in the evolution of the disease. Approximately 3000, 1000 and 300 genes were differentially expressed compared to the control groups at E15, E17 and E20, respectively. Overall, the results suggest common alterations in certain pathways between gestational time points, such as upregulation in p53 and sonic hedgehog signaling at E15 and E17 and downregulation in the myelin sheath at E17 and E20. However, there were other modifications specific to gestational time points, including skeletal muscle development at E15, downregulated glucose metabolism at E17, and upregulated inflammation at E20. In conclusion, this work provides evidence that gestational age during spina bifida repair may be a significant variable to consider during the development of new regenerative therapeutics approaches.

Comprehensive Transcriptome and Pathway Analyses Revealed Central Role for Fascin in Promoting Triple-Negative Breast Cancer Progression

Recent years have witnessed major progress in development of novel therapeutic agents such as chemotherapy, targeted therapy and immune checkpoint inhibitors for breast cancer. However, cancer-related death remains high especially in triple-negative breast cancer (TNBC) due limited therapeutic options. Development of targeted therapies for TNBC requires better understanding of biology and signaling networks that promote disease progression. Fascin, an actin bundling protein, was identified as a key regulator of many signaling pathways that contribute to breast cancer progression. Herein, fascin ShRNA was used to generate stable fascin knockdown (FSCN1^KD) in the MDA-MB-231 TNBC cell line and then were subjected to comprehensive mRNA and miRNA transcriptome analysis. We identified 129 upregulated and 114 downregulated mRNA transcripts, while 14 miRNAs were differentially expressed in FSCN1^KD. Ingenuity pathway analysis (IPA) was used to predict the impact of differentially expressed transcripts on signaling pathways and functional categories and to construct miRNA-mRNA regulatory networks in the context of FSCN1 knockdown. Compared to FSCN1^KD, fascin-positive (FSCN1^CON) breast cancer cells showed enrichment in genes promoting cellular proliferation, migration, survival, DNA replication and repair. Expression of FSCN1^high (identified in BRCA dataset from TCGA) in conjunction with elevated expression of the top 10 upregulated or decreased expression of the top 10 downregulated genes (identified in our FSCN1^CON vs. FSCN1^KD) correlates with worst survival outcome. Taken together, these data confirmed fascin's role in promoting TNBC progression, and identified a novel opportunity for therapeutic interventions via targeting those FSCN1-related transcripts.

GeneCloudOmics: A Data Analytic Cloud Platform for High-Throughput Gene Expression Analysis

Gene expression profiling techniques, such as DNA microarray and RNA-Sequencing, have provided significant impact on our understanding of biological systems. They contribute to almost all aspects of biomedical research, including studying developmental biology, host-parasite relationships, disease progression and drug effects. However, the high-throughput data generations present challenges for many wet experimentalists to analyze and take full advantage of such rich and complex data. Here we present GeneCloudOmics, an easy-to-use web server for high-throughput gene expression analysis that extends the functionality of our previous ABioTrans with several new tools, including protein datasets analysis, and a web interface. GeneCloudOmics allows both microarray and RNA-Seq data analysis with a comprehensive range of data analytics tools in one package that no other current standalone software or web-based tool can do. In total, GeneCloudOmics provides the user access to 23 different data analytical and bioinformatics tasks including reads normalization, scatter plots, linear/non-linear correlations, PCA, clustering (hierarchical, k-means, t-SNE, SOM), differential expression analyses, pathway enrichments, evolutionary analyses, pathological analyses, and protein-protein interaction (PPI) identifications. Furthermore, GeneCloudOmics allows the direct import of gene expression data from the NCBI Gene Expression Omnibus database. The user can perform all tasks rapidly through an intuitive graphical user interface that overcomes the hassle of coding, installing tools/packages/libraries and dealing with operating systems compatibility and version issues, complications that make data analysis tasks challenging for biologists. Thus, GeneCloudOmics is a one-stop open-source tool for gene expression data analysis and visualization. It is freely available at http://combio-sifbi.org/GeneCloudOmics.

Germline DDX41 mutations cause ineffective hematopoiesis and myelodysplasia

DDX41 mutations are the most common germline alterations in adult myelodysplastic syndromes (MDSs). The majority of affected individuals harbor germline monoallelic frameshift DDX41 mutations and subsequently acquire somatic mutations in their other DDX41 allele, typically missense R525H. Hematopoietic progenitor cells (HPCs) with biallelic frameshift and R525H mutations undergo cell cycle arrest and apoptosis, causing bone marrow failure in mice. Mechanistically, DDX41 is essential for small nucleolar RNA (snoRNA) processing, ribosome assembly, and protein synthesis. Although monoallelic DDX41 mutations do not affect hematopoiesis in young mice, a subset of aged mice develops features of MDS. Biallelic mutations in DDX41 are observed at a low frequency in non-dominant hematopoietic stem cell clones in bone marrow (BM) from individuals with MDS. Mice chimeric for monoallelic DDX41 mutant BM cells and a minor population of biallelic mutant BM cells develop hematopoietic defects at a younger age, suggesting that biallelic DDX41 mutant cells are disease modifying in the context of monoallelic DDX41 mutant BM.

Gain-of-function cardiomyopathic mutations in RBM20 rewire splicing regulation and re-distribute ribonucleoprotein granules within processing bodies

Mutations in the cardiac splicing factor RBM20 lead to malignant dilated cardiomyopathy (DCM). To understand the mechanism of RBM20-associated DCM, we engineered isogenic iPSCs with DCM-associated missense mutations in RBM20 as well as RBM20 knockout (KO) iPSCs. iPSC-derived engineered heart tissues made from these cell lines recapitulate contractile dysfunction of RBM20-associated DCM and reveal greater dysfunction with missense mutations than KO. Analysis of RBM20 RNA binding by eCLIP reveals a gain-of-function preference of mutant RBM20 for 3' UTR sequences that are shared with amyotrophic lateral sclerosis (ALS) and processing-body associated RNA binding proteins (FUS, DDX6). Deep RNA sequencing reveals that the RBM20 R636S mutant has unique gene, splicing, polyadenylation and circular RNA defects that differ from RBM20 KO. Super-resolution microscopy verifies that mutant RBM20 maintains very limited nuclear localization potential; rather, the mutant protein associates with cytoplasmic processing bodies (DDX6) under basal conditions, and with stress granules (G3BP1) following acute stress. Taken together, our results highlight a pathogenic mechanism in cardiac disease through splicing-dependent and -independent pathways.

Ontogeny and function of the circadian clock in intestinal organoids

Circadian rhythms regulate diverse aspects of gastrointestinal physiology ranging from the composition of microbiota to motility. However, development of the intestinal circadian clock and detailed mechanisms regulating circadian physiology of the intestine remain largely unknown. In this report, we show that both pluripotent stem cell-derived human intestinal organoids engrafted into mice and patient-derived human intestinal enteroids possess circadian rhythms and demonstrate circadian phase-dependent necrotic cell death responses to Clostridium difficile toxin B (TcdB). Intriguingly, mouse and human enteroids demonstrate anti-phasic necrotic cell death responses to TcdB. RNA-Seq analysis shows that ˜3-10% of the detectable transcripts are rhythmically expressed in mouse and human enteroids. Remarkably, we observe anti-phasic gene expression of Rac1, a small GTPase directly inactivated by TcdB, between mouse and human enteroids, and disruption of Rac1 abolishes clock-dependent necrotic cell death responses. Our findings uncover robust functions of circadian rhythms regulating clock-controlled genes in both mouse and human enteroids governing organism-specific, circadian phase-dependent necrotic cell death responses, and lay a foundation for human organ- and disease-specific investigation of clock functions using human organoids for translational applications.

Alternative RNA splicing in stem cells and cancer stem cells: Importance of transcript-based expression analysis

Alternative ribonucleic acid (RNA) splicing can lead to the assembly of different protein isoforms with distinctive functions. The outcome of alternative splicing (AS) can result in a complete loss of function or the acquisition of new functions. There is a gap in knowledge of abnormal RNA splice variants promoting cancer stem cells (CSCs), and their prospective contribution in cancer progression. AS directly regulates the self-renewal features of stem cells (SCs) and stem-like cancer cells. Notably, octamer-binding transcription factor 4A spliced variant of octamer-binding transcription factor 4 contributes to maintaining stemness properties in both SCs and CSCs. The epithelial to mesenchymal transition pathway regulates the AS events in CSCs to maintain stemness. The alternative spliced variants of CSCs markers, including cluster of differentiation 44, aldehyde dehydrogenase, and doublecortin-like kinase, α6β1 integrin, have pivotal roles in increasing self-renewal properties and maintaining the pluripotency of CSCs. Various splicing analysis tools are considered in this study. LeafCutter software can be considered as the best tool for differential splicing analysis and identification of the type of splicing events. Additionally, LeafCutter can be used for efficient mapping splicing quantitative trait loci. Altogether, the accumulating evidence re-enforces the fact that gene and protein expression need to be investigated in parallel with alternative splice variants.

Antihypertensive drug treatment and susceptibility to SARS-CoV-2 infection in human PSC-derived cardiomyocytes and primary endothelial cells

The pathogenicity of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has been attributed to its ability to enter through the membrane-bound angiotensin-converting enzyme 2 (ACE2) receptor. Therefore, it has been heavily speculated that angiotensin-converting enzyme inhibitor (ACEI) or angiotensin receptor blocker (ARB) therapy may modulate SARS-CoV-2 infection. In this study, exposure of human pluripotent stem cell-derived cardiomyocytes (hPSC-CMs) and human endothelial cells (hECs) to SARS-CoV-2 identified significant differences in protein coding genes involved in immunity, viral response, and cardiomyocyte/endothelial structure. Specifically, transcriptome changes were identified in the tumor necrosis factor (TNF), interferon α/β, and mitogen-activated protein kinase (MAPK) (hPSC-CMs) as well as nuclear factor kappa-B (NF-κB) (hECs) signaling pathways. However, pre-treatment of hPSC-CMs or hECs with two widely prescribed antihypertensive medications, losartan and lisinopril, did not affect the susceptibility of either cell type to SARS-CoV-2 infection. These findings demonstrate the toxic effects of SARS-CoV-2 in hPSC-CMs/hECs and, taken together with newly emerging multicenter trials, suggest that antihypertensive drug treatment alone does not alter SARS-CoV-2 infection.

Transcriptional alterations of protein coding and noncoding RNAs in triple negative breast cancer in response to DNA methyltransferases inhibition

Background

DNA methylation plays a crucial role in multiple cellular processes such as gene regulation, chromatin stability, and genetic imprinting. In mammals, DNA methylation is achieved by DNA methyltransferases (DNMTs). A number of studies have associated alterations in DNMT activity to tumorigenesis; however, the exact role of DNMTs in shaping the genome in triple negative breast cancer (TNBC) is still being unraveled.

Methods

In the current study, we employed two DNMT inhibitors (Decitabine and 5-Azacytidine), two TNBC models (MDA-MB-231 and BT-549) and whole transcriptome RNA-Seq and characterized the transcriptional alterations associated with DNMT inhibition. Colony forming unit (CFU), flow cytometry, and fluorescent microscopy were used to assess cell proliferation, cell cycle distribution, and cell death, respectively. Ingenuity pathway analysis (IPA) was used for network and pathway analyses.

Results

Remarkably, DNMT inhibition induced the expression of genes involved in endoplasmic reticulum response to stress, response to unfolder protein, as well as cobalamin metabolic processes. In contrast, suppression of cellular processes related to cell cycle and mitosis were hallmarks of DNMT inhibition. Concordantly, DNMT inhibition led to significant inhibition of TNBC cell proliferation, G2-M cell cycle arrest and induction of cell death. Mechanistically, DNMT inhibition activated TP53, NUPR1, and NFkB (complex) networks, while RARA, RABL6, ESR1, FOXM1, and ERBB2 networks were suppressed. Our data also identified the long noncoding RNA (lncRNA) transcriptional portrait associated with DNMT inhibition and identified 25 commonly upregulated and 60 commonly downregulated lncRNAs in response to Decitabine and 5-Azacytidinec treatment in both TNBC models. TPT1-AS1 was the most highly induced (6.3 FC), while MALAT1 was the most highly suppressed (− 7.0 FC) lncRNA in response to DNMT inhibition.

Conclusions

Taken together, our data provides a comprehensive view of transcriptome alterations in the coding and noncoding transcriptome in TNBC in response to DNMT inhibition.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12935-021-02213-2.

OncoSplicing: an updated database for clinically relevant alternative splicing in 33 human cancers

Alternative splicing (AS) represents a crucial method in mRNA level to regulate gene expression and contributes to the protein complexity. Abnormal splicing has been reported to play roles in several diseases, including cancers. We developed the OncoSplicing database for visualization of survival-associated and differential alternative splicing in 2019. Here, we provide an updated version of OncoSplicing for an integrative view of clinically relevant alternative splicing based on 122 423 AS events across 33 cancers in the TCGA SpliceSeq project and 238 558 AS events across 32 cancers in the TCGA SplAdder project. The new version of the database contains several useful features, such as annotation of alternative splicing-associated transcripts, survival analysis based on median and optimal cut-offs, differential analysis between TCGA tumour samples and adjacent normal samples or GTEx normal samples, pan-cancer views of alternative splicing, splicing differences and results of Cox’PH regression, identification of clinical indicator-relevant and cancer-specific splicing events, and downloadable splicing data in the SplAdder project. Overall, the substantially updated version of OncoSplicing (www.oncosplicing.com) is a user-friendly and registration-free database for browsing and searching clinically relevant alternative splicing in human cancers.

Unique molecular characteristics and microglial origin of Kv1.3 channel-positive brain myeloid cells in Alzheimer's disease

Kv1.3 potassium channels, expressed by proinflammatory central nervous system mononuclear phagocytes (CNS-MPs), are promising therapeutic targets for modulating neuroinflammation in Alzheimer's disease (AD). The molecular characteristics of Kv1.3-high CNS-MPs and their cellular origin from microglia or CNS-infiltrating monocytes are unclear. While Kv1.3 blockade reduces amyloid beta (Aβ) burden in mouse models, the downstream immune effects on molecular profiles of CNS-MPs remain unknown. We show that functional Kv1.3 channels are selectively expressed by a subset of CD11b⁺CD45⁺ CNS-MPs acutely isolated from an Aβ mouse model (5xFAD) as well as fresh postmortem human AD brain. Transcriptomic profiling of purified CD11b⁺Kv1.3⁺ CNS-MPs, CD11b⁺CD45^int Kv1.3^neg microglia, and peripheral monocytes from 5xFAD mice revealed that Kv1.3-high CNS-MPs highly express canonical microglial markers (Tmem119, P2ry12) and are distinct from peripheral Ly6c^high/Ly6c^low monocytes. Unlike homeostatic microglia, Kv1.3-high CNS-MPs express relatively lower levels of homeostatic genes, higher levels of CD11c, and increased levels of glutamatergic transcripts, potentially representing phagocytic uptake of neuronal elements. Using irradiation bone marrow CD45.1/CD45.2 chimerism in 5xFAD mice, we show that Kv1.3⁺ CNS-MPs originate from microglia and not blood-derived monocytes. We show that Kv1.3 channels regulate membrane potential and early signaling events in microglia. Finally, in vivo blockade of Kv1.3 channels in 5xFAD mice by ShK-223 reduced Aβ burden, increased CD11c⁺ CNS-MPs, and expression of phagocytic genes while suppressing proinflammatory genes (IL1b). Our results confirm the microglial origin and identify unique molecular features of Kv1.3-expressing CNS-MPs. In addition, we provide evidence for CNS immunomodulation by Kv1.3 blockers in AD mouse models resulting in a prophagocytic phenotype.

Detecting Rewiring Events in Protein-Protein Interaction Networks Based on Transcriptomic Data

Proteins rarely carry out their cellular functions in isolation. Instead, eukaryotic proteins engage in about six interactions with other proteins on average. The aggregated protein interactome of an organism forms a “hairy ball”-type protein-protein interaction (PPI) network. Yet, in a typical human cell, only about half of all proteins are expressed at a particular time. Hence, it has become common practice to prune the full PPI network to the subset of expressed proteins. If RNAseq data is available, one can further resolve the specific protein isoforms present in a cell or tissue. Here, we review various approaches, software tools and webservices that enable users to construct context-specific or tissue-specific PPI networks and how these are rewired between two cellular conditions. We illustrate their different functionalities on the example of the interactions involving the human TNR6 protein. In an outlook, we describe how PPI networks may be integrated with epigenetic data or with data on the activity of splicing factors.