Most human transcription factor (TF) genes encode multiple protein isoforms differing in DNA-binding domains, effector domains, or other protein regions. The global extent to which this results in functional differences between isoforms remains unknown. Here, we systematically compared 693 isoforms of 246 TF genes, assessing DNA binding, protein binding, transcriptional activation, subcellular localization, and condensate formation. Relative to reference isoforms, two-thirds of alternative TF isoforms exhibit differences in one or more molecular activities, which often could not be predicted from sequence. We observed two primary categories of alternative TF isoforms: "rewirers" and "negative regulators," both of which were associated with differentiation and cancer. Our results support a model wherein the relative expression levels of, and interactions involving, TF isoforms add an understudied layer of complexity to gene regulatory networks, demonstrating the importance of isoform-aware characterization of TF functions and providing a rich resource for further studies.

In hepatocellular carcinoma (HCC), successful translation of experimental targets identified in mouse models to human patients has proven challenging. In this study, we used a comprehensive transcriptomic approach in mice to identify novel potential targets for therapeutic intervention in humans.

We analyzed combined genome-wide miRNA and mRNA expression data in three pathogenically distinct mouse models of liver cancer. Effects of target genes on hepatoma cell fitness were evaluated by proliferation, survival and motility assays. TCGA and GEO databases, in combination with tissue microarrays, were used to validate the mouse targets and their impact on human HCC prognosis. Finally, the functional effects of the identified targets on tumorigenesis and tumor therapy were tested in hydrodynamic tail vein injection-based preclinical HCC models in vivo.

The expression of miR-107 was found to be significantly reduced in mouse models of liver tumors of various etiologies and in cohorts of humans with HCC. Overexpression of miR-107 or inhibition of its novel target kinesin family member 23 (Kif23) significantly reduced proliferation by interfering with cytokinesis, thereby controlling survival and motility of mouse and human hepatoma cells. In humans, KIF23 expression was found to be a prognostic marker in liver cancer, with high expression associated with poor prognosis. Hydrodynamic tail vein injection of vectors carrying either pre-miR-107 or anti-Kif23 shRNA inhibited the development of highly aggressive c-Myc-NRAS-induced liver cancers in mice.

Disruption of the miR-107/Kif23 axis inhibited hepatoma cell proliferation in vitro and prevented oncogene-induced liver cancer development in vivo, offering a novel potential avenue for the treatment of HCC in humans.

Our study revealed the central role of the miR-107/KIF23 axis in controlling tumor cell fitness and hepatocellular carcinoma progression. The results demonstrate that the overexpression of miR-107 or silencing of its target, KIF23, markedly suppresses the proliferation, survival, and motility of human and mouse hepatoma cells. In this work, we demonstrate that the disruption of miR-107/Kif23 signaling effectively protects mice from an aggressive form of oncogene-induced liver cancer in vivo, implying that targeting miR-107/KIF23 might be a novel therapeutic approach for hepatocellular carcinoma in humans.

The discovery of tumor-derived neoantigens which elicit an immune response through major histocompatibility complex (MHC-I/II) binding has led to significant advancements in immunotherapy. While many neoantigens have been discovered through the identification of non-synonymous mutations, the rate of these is low in some cancers, including head and neck squamous cell carcinoma. Therefore, the identification of neoantigens through additional means, such as aberrant splicing, is necessary. To achieve this, we developed the splice isoform neoantigen evaluator (SINE) pipeline. Our tool documents peptides present on spliced or inserted genomic regions of interest using Patient Harmonic-mean Best Rank scores, calculating the MHC-I/II binding affinity across the complete human leukocyte antigen landscape. Here, we found 125 potentially immunogenic events and 9 principal binders in a cohort of head and neck cancer patients where the corresponding wild-type peptides display no MHC-I/II affinity. Further, in a melanoma cohort of patients treated with anti-PD1 therapy, the expression of immunogenic splicing events identified by SINE predicted response, potentially indicating the existence of immune editing in these tumors. Overall, we demonstrate SINE's ability to identify clinically relevant immunogenic neojunctions, thus acting as a useful tool for researchers seeking to understand the neoantigen landscape from aberrant splicing in cancer.

Background. Very low-calorie ketogenic diets (VLCKD) are an effective weight-loss strategy for obese individuals, reducing risks of liver conditions such as non-alcoholic steatohepatitis and fibrosis. Small extracellular vesicles (sEVs) are implicated in liver fibrosis by influencing hepatic cell phenotypes and contributing to liver damage. This study investigates sEVs derived from serum of 60 obese adults categorized into low fibrosis risk (LR) and intermediate/high fibrosis risk (IHR) groups based on FibroScan elastography (FIB E scores, limit value 8 kPa) and all participants underwent an 8-week VLCKD intervention. Methods. The study examines the impact of these sEVs on fibrosis markers, inflammation, and autophagy in a hepatocyte cell line (HEPA-RG) using bioinformatics, RNA sequencing, lipidomics, RT-PCR, and Western blotting before (T0) and after (T1) VLCKD. Results. sEVs from LR patients post-VLCKD reduced fibrosis related gene expression (e.g., ACTA2) and enhanced proteins associated with regeneration and inflammation (e.g., HDAC6). Conversely, sEVs from IHR patients increased fibrosis and inflammation related gene expression (PIK3CB, AKT1, ACTA2) in hepatocytes, raising concerns about VLCKD suitability for IHR patients. IHR sEVs also decreased expression of HDAC10, HDAC6, HDAC3, MMP19, and MMP2, while increasing modulation of p-AKT, α-SMA, and VIM. Conclusion. These findings underscore the critical role of sEVs in regulating inflammation, remodeling, and hepatic stress responses, particularly in IHR patients, and suggest sEVs could complement instrumental evaluations like FibroScan in fibrosis assessment.

Cop1 encodes a ubiquitin E3 ligase that has been well preserved during evolution in both plants and metazoans. In metazoans, the C/EBP family transcription factors are targets for degradation by Cop1, and this process is regulated by the Tribbles pseudokinase family. Over-expression of Tribbles homolog 1 (Trib1) induces acute myeloid leukemia (AML) via Cop1-dependent degradation of the C/EBPα p42 isoform. Here, we induced rapid growth arrest and granulocytic differentiation of Trib1-expressing AML cells using a Cop1 conditional knockout (KO), which is associated with a transient increase in the C/EBPα p42 isoform. The growth-suppressive effect of Cop1 KO was canceled by silencing of Cebpa and reinforced by exogenous expression of the p42 isoform. Moreover, Cop1 KO improved the survival of recipients transplanted with Trib1-expressing AML cells. We further identified a marked increase in Trib1 protein expression in Cop1 KO, indicating that Trib1 is self-degraded by the Cop1 degradosome. COP1 downregulation also inhibits the proliferation of human AML cells in a TRIB1-dependent manner. Taken together, our results provide new insights into the role of Trib1/Cop1 machinery in the C/EBPα p42-dependent leukemogenic activity, and a novel idea to develop new therapeutics.

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. We hypothesized that the disease specificity of targets can be systematically learned for all genes by jointly evaluating complementary molecular measurements of healthy tissues using a hierarchical Bayesian modeling approach. Our method, BayesTS, integrates protein and gene expression evidence and includes tunable parameters to moderate tissue essentiality. Applied to all protein coding genes, BayesTS outperforms alternative strategies to define therapeutic targets and nominates previously unknown targets while allowing for incorporation of new types of modalities. To expand target repertoires, we show that extension of BayesTS to splicing antigens and combinatorial target pairs results in more specific targets for therapy. We expect that BayesTS will facilitate improved target prioritization for oncology drug development, ultimately leading to the discovery of more effective and safer treatments.

Acute myeloid leukemia (AML) is a multi-clonal disease, existing as a milieu of clones with unique but related genotypes as initiating clones acquire subsequent mutations. However, bulk sequencing cannot fully capture AML clonal architecture or the clonal evolution that occurs as patients undergo therapy. To interrogate clonal evolution, we performed simultaneous single cell molecular profiling and immunophenotyping on 43 samples from 32 NPM1-mutant AML patients at different stages of disease. Here we show that diagnosis and relapsed AML samples display similar clonal architecture patterns, but signaling mutations can drive increased clonal diversity specifically at relapse. We uncovered unique genotype-immunophenotype relationships regardless of disease state, suggesting leukemic lineage trajectories can be hard-wired by the mutations present. Analysis of longitudinal samples from patients on therapy identified dynamic clonal, transcriptomic, and immunophenotypic changes. Our studies provide resolved understanding of leukemic clonal evolution and the relationships between genotype and cell state in leukemia biology.

Cancer-associated fibroblasts (CAFs), mainly responsible for the desmoplastic reaction hallmark of intrahepatic Cholangiocarcinoma (iCCA), likely have a role in tumor aggressiveness and resistance to therapy, although the molecular mechanisms involved are unknown. Aim of the study is to investigate how targeting hCAF/iCCA cross-talk with a Notch1 inhibitor, namely Crenigacestat, may affect cancer progression.

We used different in vitro models in 2D and established new 3D hetero-spheroids with iCCA cells and human (h)CAFs. The results were confirmed in a xenograft model, and explanted tumoral tissues underwent transcriptomic and bioinformatic analysis.

hCAFs/iCCA cross-talk sustains increased migration of both KKU-M213 and KKU-M156 cells, while Crenigacestat significantly inhibits only the cross-talk stimulated migration. Hetero-spheroids grew larger than homo-spheroids, formed by only iCCA cells. Crenigacestat significantly reduced the invasion and growth of hetero- but not of homo-spheroids. In xenograft models, hCAFs/KKU-M213 tumors grew significantly larger than KKU-M213 tumors, but were significantly reduced in volume by Crenigacestat treatment, which also significantly decreased the fibrotic reaction. Ingenuity pathway analysis revealed that genes of hCAFs/KKU-M213 but not of KKU-M213 tumors increased tumor lesions, and that Crenigacestat treatment inhibited the modulated canonical pathways. Cell cycle checkpoints were the most notably modulated pathway and Crenigacestat reduced CCNE2 gene expression, consequently inducing cell cycle arrest. In hetero-spheroids, the number of cells increased in the G2/M cell cycle phase, while Crenigacestat significantly decreased cell numbers in the G2/M phase in hetero but not in homo-spheroids.

The hCAFs/iCCA cross-talk is a new target for reducing cancer progression with drugs such as Crenigacestat.

As the most abundant glial cells in the central nervous system (CNS), astrocytes dynamically respond to neurotoxic stress, however, the key molecular regulators controlling the inflammatory status of these sentinels during neurotoxic stress are many and complex. Herein, we demonstrate that the m6A epitranscriptomic mRNA modification tightly regulates the pro-inflammatory functions of astrocytes. Specifically, the astrocytic neurotoxic stressor, 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, the neurotoxic stress-induced proinflammatory response, suggesting YTHDF2 serves as a key upstream regulator of inflammatory responses in astrocytes. Mechanistically, YTHDF2 RIP-sequencing identified MAP2K4 (MKK4; SEK1) mRNA as a YTHDF2 target influencing inflammatory signaling. Our target validation revealed that Mn-exposed astrocytes mediate proinflammatory responses by activating the phosphorylation of SEK1, JNK, and cJUN signaling. Collectively, YTHDF2 serves as a key upstream 'molecular switch' controlling SEK1(MAP2K4)-JNK-cJUN proinflammatory signaling in astrocytes.

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 a lipopolysaccharide-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.

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.

The severity of COVID-19 is linked to excessive inflammation. Neutrophils represent a critical arm of the innate immune response and are major mediators of inflammation, but their role in COVID-19 pathophysiology remains poorly understood. We conducted transcriptomic profiling of neutrophils obtained from patients with mild and severe COVID-19, as well as from SARS-CoV-2 infected mice, in comparison to non-infected healthy controls. In addition, we investigated the inflammasome formation potential in neutrophils from patients and mice upon SARS-CoV-2 infection. Transcriptomic analysis of polymorphonuclear cells (PMNs), consisting mainly of mature neutrophils, revealed a striking type I interferon (IFN-I) gene signature in severe COVID-19 patients, contrasting with mild COVID-19 and healthy controls. Notably, low-density granulocytes (LDGs) from severe COVID-19 patients exhibited an immature neutrophil phenotype and lacked this IFN-I signature. Moreover, PMNs from severe COVID-19 patients showed heightened nigericin-induced caspase1 activation, but reduced responsiveness to exogenous inflammasome priming. Furthermore, IFN-I emerged as a priming stimulus for neutrophil inflammasomes. These findings suggest a potential role for neutrophil inflammasomes in driving inflammation during severe COVID-19. Altogether, these findings open promising avenues for targeted therapeutic interventions to mitigate the pathological processes associated with the disease.

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.

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.

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

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.

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 (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.

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.

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.

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

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.