The harmonizome: a collection of processed datasets gathered to serve and mine knowledge about genes and proteins

UBR5 promotes antiviral immunity by disengaging the transcriptional brake on RIG-I like receptors

The Retinoic acid-Inducible Gene I (RIG-I) like receptors (RLRs) are the major viral RNA sensors essential for the initiation of antiviral immune responses. RLRs are subjected to stringent transcriptional and posttranslational regulations, of which ubiquitination is one of the most important. However, the role of ubiquitination in RLR transcription is unknown. Here, we screen 375 definite ubiquitin ligase knockout cell lines and identify Ubiquitin Protein Ligase E3 Component N-Recognin 5 (UBR5) as a positive regulator of RLR transcription. UBR5 deficiency reduces antiviral immune responses to RNA viruses, while increases viral replication in primary cells and mice. Ubr5 knockout mice are more susceptible to lethal RNA virus infection than wild type littermates. Mechanistically, UBR5 mediates the Lysine 63-linked ubiquitination of Tripartite Motif Protein 28 (TRIM28), an epigenetic repressor of RLRs. This modification prevents intramolecular SUMOylation of TRIM28, thus disengages the TRIM28-imposed brake on RLR transcription. In sum, UBR5 enables rapid upregulation of RLR expression to boost antiviral immune responses by ubiquitinating and de-SUMOylating TRIM28.

The glucocorticoid receptor as a master regulator of the Müller cell response to diabetic conditions in mice

Diabetic retinopathy (DR) is considered a primarily microvascular complication of diabetes. Müller glia cells are at the centre of the retinal neurovascular unit and play a critical role in DR. We therefore investigated Müller cell-specific signalling pathways that are altered in DR to identify novel targets for gene therapy. Using a multi-omics approach on purified Müller cells from diabetic db/db mice, we found the mRNA and protein expression of the glucocorticoid receptor (GR) to be significantly decreased, while its target gene cluster was down-regulated. Further, oPOSSUM TF analysis and ATAC- sequencing identified the GR as a master regulator of Müller cell response to diabetic conditions. Cortisol not only increased GR phosphorylation. It also induced changes in the expression of known GR target genes in retinal explants. Finally, retinal functionality was improved by AAV-mediated overexpression of GR in Müller cells. Our study demonstrates an important role of the glial GR in DR and implies that therapeutic approaches targeting this signalling pathway should be aimed at increasing GR expression rather than the addition of more ligand.

Long Noncoding RNA VLDLR-AS1 Levels in Serum Correlate with Combat-Related Chronic Mild Traumatic Brain Injury and Depression Symptoms in US Veterans

More than 75% of traumatic brain injuries (TBIs) are mild (mTBI) and military service members often experience repeated combat-related mTBI. The chronic comorbidities concomitant with repetitive mTBI (rmTBI) include depression, post-traumatic stress disorder or neurological dysfunction. This study sought to determine a long noncoding RNA (lncRNA) expression signature in serum samples that correlated with rmTBI years after the incidences. Serum samples were obtained from Long-Term Impact of Military-Relevant Brain-Injury Consortium Chronic Effects of Neurotrauma Consortium (LIMBIC CENC) repository, from participants unexposed to TBI or who had rmTBI. Four lncRNAs were identified as consistently present in all samples, as detected via droplet digital PCR and packaged in exosomes enriched for CNS origin. The results, using qPCR, demonstrated that the lncRNA VLDLR-AS1 levels were significantly lower among individuals with rmTBI compared to those with no lifetime TBI. ROC analysis determined an AUC of 0.74 (95% CI: 0.6124 to 0.8741; p = 0.0012). The optimal cutoff for VLDLR-AS1 was ≤153.8 ng. A secondary analysis of clinical data from LIMBIC CENC was conducted to evaluate the psychological symptom burden, and the results show that lncRNAs VLDLR-AS1 and MALAT1 are correlated with symptoms of depression. In conclusion, lncRNA VLDLR-AS1 may serve as a blood biomarker for identifying chronic rmTBI and depression in patients.

Uveitic glaucoma-like features in Yap conditional knockout mice

Glaucoma is a multifactorial neurodegenerative disease characterized by the progressive and irreversible degeneration of the optic nerve and retinal ganglion cells. Despite medical advances aiming at slowing degeneration, around 40% of treated glaucomatous patients will undergo vision loss. It is thus of utmost importance to have a better understanding of the disease and to investigate more deeply its early causes. The transcriptional coactivator YAP, an important regulator of eye homeostasis, has recently drawn attention in the glaucoma research field. Here we show that Yap conditional knockout mice (Yap cKO), in which the deletion of Yap is induced in both Müller glia (i.e. the only retinal YAP-expressing cells) and the non-pigmented epithelial cells of the ciliary body, exhibit a breakdown of the aqueous-blood barrier, accompanied by a progressive collapse of the ciliary body. A similar phenotype is observed in human samples that we obtained from patients presenting with uveitis. In addition, aged Yap cKO mice harbor glaucoma-like features, including deregulation of key homeostatic Müller-derived proteins, retinal vascular defects, optic nerve degeneration and retinal ganglion cell death. Finally, transcriptomic analysis of Yap cKO retinas pointed to early-deregulated genes involved in extracellular matrix organization potentially underlying the onset and/or progression of the observed phenotype. Together, our findings reveal the essential role of YAP in preserving the integrity of the ciliary body and retinal ganglion cells, thereby preventing the onset of uveitic glaucoma-like features.

Hepatitis C virus E1 recruits high-density lipoprotein to support infectivity and evade antibody recognition

Hepatitis C virus (HCV) is a member of the Flaviviridae family; however, unlike other family members, the HCV virion has an unusually high lipid content. HCV has two envelope glycoproteins, E1 and E2. E2 contributes to receptor binding, cell membrane attachment, and immune evasion. In contrast, the functions of E1 are poorly characterized due, in part, to challenges in producing the protein. This manuscript describes the expression and purification of a soluble E1 ectodomain (eE1) that is recognized by conformational, human monoclonal antibodies. eE1 forms a complex with apolipoproteins AI and AII, cholesterol, and phospholipids by recruiting high-density lipoprotein (HDL) from the extracellular media. We show that HDL binding is a function specific to eE1 and HDL hinders recognition of E1 by a neutralizing monoclonal antibody. Either low-density lipoprotein or HDL increases the production and infectivity of cell culture-produced HCV, but E1 preferentially selects HDL, influencing both viral life cycle and antibody evasion.IMPORTANCEHepatitis C virus (HCV) infection is a significant burden on human health, but vaccine candidates have yet to provide broad protection against this infection. We have developed a method to produce high quantities of soluble E1 or E2, the viral proteins located on the surface of HCV. HCV has an unusually high lipid content due to the recruitment of apolipoproteins. We found that E1 (and not E2) preferentially recruits host high-density lipoprotein (HDL) extracellularly. This recruitment of HDL by E1 prevents binding of E1 by a neutralizing antibody and furthermore prevents antibody-mediated neutralization of the virus. By comparison, low-density lipoprotein does not protect the virus from antibody-mediated neutralization. Our findings provide mechanistic insight into apolipoprotein recruitment, which may be critical for vaccine development.

VISTA promotes the metabolism and differentiation of myeloid-derived suppressor cells by STAT3 and polyamine-dependent mechanisms

Myeloid-derived suppressor cells (MDSCs) impair antitumor immune responses. Identifying regulatory circuits during MDSC development may bring new opportunities for therapeutic interventions. We report that the V-domain suppressor of T cell activation (VISTA) functions as a key enabler of MDSC differentiation. VISTA deficiency reduced STAT3 activation and STAT3-dependent production of polyamines, which causally impaired mitochondrial respiration and MDSC expansion. In both mixed bone marrow (BM) chimera mice and myeloid-specific VISTA conditional knockout mice, VISTA deficiency significantly reduced tumor-associated MDSCs but expanded monocyte-derived dendritic cells (DCs) and enhanced T cell-mediated tumor control. Correlated expression of VISTA and arginase-1 (ARG1), a key enzyme supporting polyamine biosynthesis, was observed in multiple human cancer types. In human endometrial cancer, co-expression of VISTA and ARG1 on tumor-associated myeloid cells is associated with poor survival. Taken together, these findings unveil the VISTA/polyamine axis as a central regulator of MDSC differentiation and warrant therapeutically targeting this axis for cancer immunotherapy.

GMFGRN: a matrix factorization and graph neural network approach for gene regulatory network inference

The recent advances of single-cell RNA sequencing (scRNA-seq) have enabled reliable profiling of gene expression at the single-cell level, providing opportunities for accurate inference of gene regulatory networks (GRNs) on scRNA-seq data. Most methods for inferring GRNs suffer from the inability to eliminate transitive interactions or necessitate expensive computational resources. To address these, we present a novel method, termed GMFGRN, for accurate graph neural network (GNN)-based GRN inference from scRNA-seq data. GMFGRN employs GNN for matrix factorization and learns representative embeddings for genes. For transcription factor-gene pairs, it utilizes the learned embeddings to determine whether they interact with each other. The extensive suite of benchmarking experiments encompassing eight static scRNA-seq datasets alongside several state-of-the-art methods demonstrated mean improvements of 1.9 and 2.5% over the runner-up in area under the receiver operating characteristic curve (AUROC) and area under the precision-recall curve (AUPRC). In addition, across four time-series datasets, maximum enhancements of 2.4 and 1.3% in AUROC and AUPRC were observed in comparison to the runner-up. Moreover, GMFGRN requires significantly less training time and memory consumption, with time and memory consumed <10% compared to the second-best method. These findings underscore the substantial potential of GMFGRN in the inference of GRNs. It is publicly available at https://github.com/Lishuoyy/GMFGRN.

Standigm ASK™: knowledge graph and artificial intelligence platform applied to target discovery in idiopathic pulmonary fibrosis

Standigm ASK™ revolutionizes healthcare by addressing the critical challenge of identifying pivotal target genes in disease mechanisms-a fundamental aspect of drug development success. Standigm ASK™ integrates a unique combination of a heterogeneous knowledge graph (KG) database and an attention-based neural network model, providing interpretable subgraph evidence. Empowering users through an interactive interface, Standigm ASK™ facilitates the exploration of predicted results. Applying Standigm ASK™ to idiopathic pulmonary fibrosis (IPF), a complex lung disease, we focused on genes (AMFR, MDFIC and NR5A2) identified through KG evidence. In vitro experiments demonstrated their relevance, as TGFβ treatment induced gene expression changes associated with epithelial-mesenchymal transition characteristics. Gene knockdown reversed these changes, identifying AMFR, MDFIC and NR5A2 as potential therapeutic targets for IPF. In summary, Standigm ASK™ emerges as an innovative KG and artificial intelligence platform driving insights in drug target discovery, exemplified by the identification and validation of therapeutic targets for IPF.

Rare loss-of-function variants in matrisome genes are enriched in Ebstein's anomaly

Ebstein's anomaly, a rare congenital heart disease, is distinguished by the failure of embryological delamination of the tricuspid valve leaflets from the underlying primitive right ventricle myocardium. Gaining insight into the genetic basis of Ebstein's anomaly allows a more precise definition of its pathogenesis. In this study, two distinct cohorts from the Chinese Han population were included: a case-control cohort consisting of 82 unrelated cases and 125 controls without cardiac phenotypes and a trio cohort comprising 36 parent-offspring trios. Whole-exome sequencing data from all 315 participants were utilized to identify qualifying variants, encompassing rare (minor allele frequency < 0.1% from East Asians in the gnomAD database) functional variants and high-confidence (HC) loss-of-function (LoF) variants. Various statistical models, including burden tests and variance-component models, were employed to identify rare variants, genes, and biological pathways associated with Ebstein's anomaly. Significant associations were noted between Ebstein's anomaly and rare HC LoF variants found in genes related to the matrisome, a collection of extracellular matrix (ECM) components. Specifically, 47 genes with HC LoF variants were exclusively or predominantly identified in cases, while nine genes showed such variants in the probands. Over half of unrelated cases (n = 42) and approximately one-third of probands (n = 12) were found to carry one or two LoF variants in these prioritized genes. These results highlight the role of the matrisome in the pathogenesis of Ebstein's anomaly, contributing to a better understanding of the genetic architecture underlying this condition. Our findings hold the potential to impact the genetic diagnosis and treatment approaches for Ebstein's anomaly.

PhpC modulates G-quadruplex-RNA landscapes in human cells

Stabilizing DNA/RNA G-quadruplexes (G4s) using small molecules (ligands) has proven an efficient strategy to decipher G4 biology. Quite paradoxically, this search has also highlighted the need for finding molecules able to disrupt G4s to tackle G4-associated cellular dysfunctions. We report here on both qualitative and quantitative investigations that validate the G4-RNA-destabilizing properties of the leading compound PhpC in human cells.

The gene expression profile and cell of origin of canine peripheral T-cell lymphoma

BACKGROUND

Peripheral T-cell lymphoma (PTCL) refers to a heterogenous group of T-cell neoplasms with poor treatment responses and survival times. Canine PTCL clinically and immunophenotypically resembles the most common human subtype, PTCL-not otherwise specified (PTCL-NOS), leading to interest in this canine disease as a naturally occurring model for human PTCL. Gene expression profiling in human PTCL-NOS has helped characterize this ambiguous diagnosis into distinct subtypes, but similar gene expression profiling in canine PTCL is lacking.

METHODS

Bulk RNA-sequencing was performed on tumor samples from 33 dogs with either CD4+ (26/33), CD8+ (4/33), or CD4-CD8- (3/33) PTCL as diagnosed by flow cytometry, and sorted CD4+ and CD8+ lymphocytes from healthy control dogs. Following normalization of RNA-seq data, we performed differential gene expression and unsupervised clustering methods. Gene set enrichment analysis was performed to determine the enrichment of canine CD4+ PTCL for human PTCL-NOS, oncogenic pathways, and various stages of T-cell development gene signatures. We utilized gene set variation analysis to evaluate individual canine CD4+ PTCLs for various human and murine T-cell and thymocyte gene signatures. Cultured canine PTCL cells were treated with a pan-PI3K inhibitor, and cell survival and proliferation were compared to DMSO-treated controls. Expression of GATA3 and phosphorylated AKT was validated by immunohistochemistry.

RESULTS

While the canine CD4+ PTCL phenotype exhibited a consistent gene expression profile, the expression profiles of CD8+ and CD4-CD8- canine PTCLs were more heterogeneous. Canine CD4+ PTCL had increased expression of GATA3, upregulation of its target genes, enrichment for PI3K/AKT/mTOR signaling, and downregulation of PTEN, features consistent with the more aggressive GATA3-PTCL subtype of human PTCL-NOS. In vitro assays validated the reliance of canine CD4+ PTCL cells on PI3K/AKT/mTOR signaling for survival and proliferation. Canine CD4+ PTCL was enriched for thymic precursor gene signatures, exhibited increased expression of markers of immaturity (CD34, KIT, DNTT, and CCR9), and downregulated genes associated with the T-cell receptor, MHC class II associated genes (DLA-DQA1, DLA-DRA, HLA-DQB1, and HLA-DQB2), and CD25.

CONCLUSIONS

Canine CD4+ PTCL most closely resembled the GATA3-PTCL subtype of PTCL-NOS and may originate from an earlier stage of T-cell development than the more conventionally posited mature T-helper cell origin.

Carcinogenic roles of MAFG-AS1 in human cancers

The MAF bZIP transcription factor G-antisense RNA 1 (MAFG-AS1) is located on chromosome 17. MAFG-AS1 was upregulated in 15 human cancers. MAFG-AS1 not only suppresses 16 miRNAs but also directly impacts 22 protein-coding genes' expression. Notably, abnormal MAFG-AS1 expression is connected to clinicopathological characteristics and a worse prognosis in a variety of cancers. Moreover, MAFG-AS1 takes its part in the tumorigenesis and progression of various human malignancies by suppressing apoptosis and promoting proliferation, migration, invasion, aerobic glycolysis, ferroptosis, angiogenesis, EMT, and metastasis. Besides, it can predict treatment effectiveness in ER + breast cancer, urothelial bladder carcinoma, and liver cancer by functioning as a trigger of resistance to tamoxifen, sorafenib, and cisplatin. This study systematically presents the functions of MAFG-AS1 in various cancers, as well as the findings of bioinformatics analyses of the MAFG-AS1, which should give clear advice for future research.

Aggressive high-grade NF2 mutant meningiomas downregulate oncogenic YAP signaling via the upregulation of VGLL4 and FAT3/4

Background

Meningiomas are the most common primary central nervous system tumors in adults. Although generally benign, a subset is of higher grade and ultimately fatal. Around half of all meningiomas harbor inactivating mutations in NF2, leading to deregulation of oncogenic YAP1 activity. While benign NF2 mutant meningiomas exhibit few genetic events in addition to NF2 inactivation, aggressive high-grade NF2 mutant meningiomas frequently harbor a highly aberrant genome. It is unclear if NF2 mutant meningiomas of different grades are equally reliant on YAP activity.

Methods

We analyzed bulk and single-cell RNA-Seq data from a large cohort of human meningiomas for the expression of YAP1 target genes and Hippo effectors as well as in vitro cell line experiments.

Results

Aggressive NF2 mutant meningiomas harbor decreased expression levels of YAP1 target genes and increased expression levels of the YAP1 antagonist VGLL4 and the upstream regulators FAT3/4 compared to their benign counterparts. Decreased expression of YAP1 target genes as well as high expression of VGLL4 and FAT3/4 is significantly associated with an increased risk of recurrence. In vitro, overexpression of VGLL4 resulted in the downregulation of YAP activity in benign NF2 mutant meningioma cells, confirming the direct link between VGLL4 expression and decreased levels of YAP activity observed in aggressive NF2 mutant meningiomas.

Conclusions

Our results shed new insight into the biology of benign and aggressive NF2 mutant meningiomas and may have important implications for the efficacy of therapies targeting oncogenic YAP1 activity in NF2 mutant meningiomas.

Phenelzine protects against acetaminophen induced apoptosis in HepG2 cells

Acetaminophen (APAP) overdosing is the most common cause of drug-induced liver failure. Despite extensive study, N-acetylcysteine is currently the only antidote utilized for treatment. The purpose of this study was to evaluate the effect and mechanisms of phenelzine, an FDA-approved antidepressant, on APAP-induced toxicity in HepG2 cells. The human liver hepatocellular cell line HepG2 was used to investigate APAP-induced cytotoxicity. The protective effects of phenelzine were determined by examining the cell viability, combination index calculation, Caspase 3/7 activation, Cytochrome c release, H2O2 levels, NO levels, GSH activity, PERK protein levels, and pathway enrichment analysis. Elevated H2O2 production and decreased glutathione (GSH) levels were indicators of APAP-induced oxidative stress. The combination index of 2.04 indicated that phenelzine had an antagonistic effect on APAP-induced toxicity. When compared to APAP alone, phenelzine treatment considerably reduced caspase 3/7 activation, cytochrome c release, and H2O2 generation. However, phenelzine had minimal effect on NO and GSH levels and did not alleviate ER stress. Pathway enrichment analysis revealed a potential connection between APAP toxicity and phenelzine metabolism. These findings suggested that phenelzine's protective effect against APAP-induced cytotoxicity could be attributed to the drug's capacity to reduce APAP-mediated apoptotic signaling.

Insulin-like growth factor 2 mRNA-binding protein 3 enhanced melanoma migration through regulation of AKT1 and RELA expression

IMP-3 expression is a poor prognostic factor of melanomas and it promotes melanoma cell migration and invasion by a pathway modulating HMGA2 mRNA expression. We tried to identify other putative targets of IMP-3. We identified putative IMP-3-binding RNAs, including AKT1, MAPK3, RB1 and RELA, by RNA immunoprecipitation coupled with next-generation sequencing. IMP-3 overexpression increased AKT and RELA levels in MeWo cells. siRNAs against AKT1 and RELA inhibited MeWo/Full-length IMP-3 cell migration. IMP-3 knockdown of A2058 cells decreased AKT1 and RELA expression and lowered migration ability. Co-transfection of A2058 cells with AKT1- or RELA-expressing plasmids with IMP-3 siRNA restored the inhibitory effects of IMP-3 knockdown on migration. HMGA2 did not influence AKT1 and RELA expression in melanoma cells. Human melanoma samples with high IMP-3 levels also showed high HMGA2, AKT1 and RELA expression. Our results show that IMP-3 enhances melanoma cell migration through the regulation of the AKT1 and RELA axis.

Anoikis-related CTNND1 is associated with immuno-suppressive tumor microenvironment and predicts unfavorable immunotherapeutic outcome in non-small cell lung cancer

Background: Immunotherapy has greatly changed the treatment of advanced non-small cell lung cancer (NSCLC). Anoikis is a programmed cell death process associated with cancer. However, the correlation between anoikis-related genes and the tumor microenvironment (TME) features and immunotherapeutic outcome in NSCLC has not been fully explored. Methods: The bulk and single-cell transcriptome data of NSCLC were downloaded from TCGA and GEO databases. The distribution of anoikis-related genes on different cell types at the single-cell level was analyzed, and these genes specifically expressed by tumor cells and immunotherapy-related were further extracted. Next, the candidate gene CTNND1 was identified and its correlations with the TME features and immunotherapeutic outcome in NSCLC were explored in multiple public cohorts. Finally, an in-house cohort was used to determine the CTNND1 expression and immuno-correlation in NSCLC. Results: At single-cell atlas, we found that anoikis-related genes expressed specifically in tumor cells of NSCLC. By intersecting anoikis-related genes, immunotherapy-associated genes, and the genes expressed in tumor cells, we obtained a special biomarker CTNND1. In addition, cell-cell communication analysis revealed that CTNND1+ tumor cells communicated with immune subpopulations frequently. Moreover, we found that high expression of CTNND1 was related to immuno-suppressive status of NSCLC. The expression of CTNND1 and its immuno-correlation were also validated, and the results showed that CTNND1 was highly expressed in NSCLC tissues and tumors with high CTNND1 expression accompanied with low CD8+ T cells infiltration. Conclusions: Overall, our study reported that CTNND1 can be considered as a novel biomarker for the predication of immunotherapeutic responses and a potential target for NSCLC therapy.

Endothelial YAP Mediates Hyperglycemia-Induced Platelet Hyperactivity and Arterial Thrombosis

BACKGROUND

Hyperglycemia-a symptom that characterizes diabetes-is highly associated with atherothrombotic complications. However, the underlying mechanism by which hyperglycemia fuels platelet activation and arterial thrombus formation is still not fully understood.

METHODS

The profiles of polyunsaturated fatty acid metabolites in the plasma of patients with diabetes and healthy controls were determined with targeted metabolomics. FeCl3-induced carotid injury model was used to assess arterial thrombus formation in mice with endothelial cell (EC)-specific YAP (yes-associated protein) deletion or overexpression. Flow cytometry and clot retraction assay were used to evaluate platelet activation. RNA sequencing and multiple biochemical analyses were conducted to unravel the underlying mechanism.

RESULTS

The plasma PGE2 (prostaglandin E2) concentration was elevated in patients with diabetes with thrombotic complications and positively correlated with platelet activation. The PGE2 synthetases COX-2 (cyclooxygenase-2) and mPGES-1 (microsomal prostaglandin E synthase-1) were found to be highly expressed in ECs but not in other type of vessel cells in arteries from both patients with diabetes and hyperglycemic mice, compared with nondiabetic individuals and control mice, respectively. A combination of RNA sequencing and ingenuity pathway analyses indicated the involvement of YAP signaling. EC-specific deletion of YAP limited platelet activation and arterial thrombosis in hyperglycemic mice, whereas EC-specific overexpression of YAP in mice mimicked the prothrombotic state of diabetes, without affecting hemostasis. Mechanistically, we found that hyperglycemia/high glucose-induced endothelial YAP nuclear translocation and subsequently transcriptional expression of COX-2 and mPGES-1 contributed to the elevation of PGE2 and platelet activation. Blockade of EP3 (prostaglandin E receptor 3) activation by oral administration of DG-041 reversed the hyperactivity of platelets and delayed thrombus formation in both EC-specific YAP-overexpressing and hyperglycemic mice.

CONCLUSIONS

Collectively, our data suggest that hyperglycemia-induced endothelial YAP activation aggravates platelet activation and arterial thrombus formation via PGE2/EP3 signaling. Targeting EP3 with DG-041 might be therapeutic for diabetes-related thrombosis.

Incorporating Sequence-Dependent DNA Shape and Dynamics into Transcriptome Data Analysis

Differentially expressed genes in a cellular context may be co-regulated by the same transcription factor. However, in the absence of a concurrent transcription factor binding data, such interactions are difficult to detect, especially at the single cell expression level. Motif enrichments in such genes can be used to gain insight into differential expressions caused by the shared upstream TFs. However, it is now established that many genes are co-regulated by the same TF due to a shared DNA shape or sequence-dependent conformational dynamics instead of sequence motif. In this work, we demonstrate how, starting from a gene expression data, such DNA shape and dynamics signatures can be potentially detected using publicly available tools, including DynaSeq, developed in our group for predicting the sequence-dependent components of these DNA shape features.

Cerebrospinal fluid proteomics in patients with Alzheimer's disease reveals five molecular subtypes with distinct genetic risk profiles

Alzheimer's disease (AD) is heterogenous at the molecular level. Understanding this heterogeneity is critical for AD drug development. Here we define AD molecular subtypes using mass spectrometry proteomics in cerebrospinal fluid, based on 1,058 proteins, with different levels in individuals with AD (n = 419) compared to controls (n = 187). These AD subtypes had alterations in protein levels that were associated with distinct molecular processes: subtype 1 was characterized by proteins related to neuronal hyperplasticity; subtype 2 by innate immune activation; subtype 3 by RNA dysregulation; subtype 4 by choroid plexus dysfunction; and subtype 5 by blood-brain barrier impairment. Each subtype was related to specific AD genetic risk variants, for example, subtype 1 was enriched with TREM2 R47H. Subtypes also differed in clinical outcomes, survival times and anatomical patterns of brain atrophy. These results indicate molecular heterogeneity in AD and highlight the need for personalized medicine.

Mapping of the gene network that regulates glycan clock of ageing

Glycans are an essential structural component of immunoglobulin G (IgG) that modulate its structure and function. However, regulatory mechanisms behind this complex posttranslational modification are not well known. Previous genome-wide association studies (GWAS) identified 29 genomic regions involved in regulation of IgG glycosylation, but only a few were functionally validated. One of the key functional features of IgG glycosylation is the addition of galactose (galactosylation), a trait which was shown to be associated with ageing. We performed GWAS of IgG galactosylation (N=13,705) and identified 16 significantly associated loci, indicating that IgG galactosylation is regulated by a complex network of genes that extends beyond the galactosyltransferase enzyme that adds galactose to IgG glycans. Gene prioritization identified 37 candidate genes. Using a recently developed CRISPR/dCas9 system we manipulated gene expression of candidate genes in the in vitro IgG expression system. Upregulation of three genes, EEF1A1, MANBA and TNFRSF13B, changed the IgG glycome composition, which confirmed that these three genes are involved in IgG galactosylation in this in vitro expression system.

Histone H2A variant H2A.B is enriched in transcriptionally active HSV-1 lytic chromatin

Herpes simplex virus 1 (HSV-1) transcription is restricted in latently infected neurons and the genomes are in mostly silenced chromatin, whereas all viral genes are transcribed in lytically infected cells, in which the genomes are dynamically chromatinized. Epigenetic regulation modulates HSV-1 transcription during lytic, latent, and reactivating infections, but the precise mechanisms are not fully defined. Nucleosomes are dynamic; they slide, breathe, assemble and disassemble. We and others have proposed that the most dynamic HSV-1 chromatin is transcriptionally competent whereas the least dynamic is silenced. However, the mechanisms yielding the unusually dynamic viral chromatin remain unknown. Histone variants affect nucleosome dynamics. The dynamics of H2A, H2A.X and macroH2A were enhanced in infected cells, whereas those of H2A.B uniquely decreased. We constructed stably transduced cells expressing tagged histone H2A, H2A.B, macroH2A, or H2B, which assembles the H2A/H2B nucleosome dimers with all H2A variants. All H2A variants, ectopic, and endogenous H2B, were assembled into HSV-1 chromatin evenly throughout the genome, but canonical H2A was relatively depleted from the viral chromatin whereas H2A.B was enriched in the most dynamic viral chromatin. When viral transcription was restricted, H2A.B became as depleted from the viral chromatin through the entire genome as H2A. We propose that lytic HSV-1 nucleosomes are enriched in the dynamic variant H2A.B/H2B dimers to promote HSV-1 chromatin dynamics and transcriptional competency, and conclude that the dynamics of HSV-1 chromatin are determined in part by the H2A variants.

Importance

HSV-1 transcription is epigenetically regulated during latent and lytic infections, and epigenetic inhibitors have been proposed as potential antiviral drugs to modulate latency and reactivation. However, the detailed mechanisms of regulation of HSV-1 transcription by epigenetics have not been fully characterized and may differ from those regulating cellular transcription. In particular, the lytic HSV-1 chromatin is unusually dynamic, whereas the latent silenced one is not, but the mechanisms resulting in the unique dynamics of the lytic chromatin remain unknown. Here we identify the enrichment on the highly dynamic histone 2A variant H2A in the most dynamic viral chromatin, which provides a mechanistic understanding for its unique dynamics. Future work to identify the mechanisms of enrichment in H2A.B on the viral chromatin may identify novel druggable epigenetic regulators that modulate HSV-1 latency and reactivation.

Deciphering driver regulators of cell fate decisions from single-cell transcriptomics data with CEFCON

Single-cell technologies enable the dynamic analyses of cell fate mapping. However, capturing the gene regulatory relationships and identifying the driver factors that control cell fate decisions are still challenging. We present CEFCON, a network-based framework that first uses a graph neural network with attention mechanism to infer a cell-lineage-specific gene regulatory network (GRN) from single-cell RNA-sequencing data, and then models cell fate dynamics through network control theory to identify driver regulators and the associated gene modules, revealing their critical biological processes related to cell states. Extensive benchmarking tests consistently demonstrated the superiority of CEFCON in GRN construction, driver regulator identification, and gene module identification over baseline methods. When applied to the mouse hematopoietic stem cell differentiation data, CEFCON successfully identified driver regulators for three developmental lineages, which offered useful insights into their differentiation from a network control perspective. Overall, CEFCON provides a valuable tool for studying the underlying mechanisms of cell fate decisions from single-cell RNA-seq data.

Club cell CREB regulates the goblet cell transcriptional network and pro-mucin effects of IL-1B

Introduction: Club cells are precursors for mucus-producing goblet cells. Interleukin 1β (IL-1B) is an inflammatory mediator with pro-mucin activities that increases the number of mucus-producing goblet cells. IL-1B-mediated mucin production in alveolar adenocarcinoma cells requires activation of the cAMP response element-binding protein (CREB). Whether the pro-mucin activities of IL-1B require club cell CREB is unknown. Methods: We challenged male mice with conditional loss of club cell Creb1 and wild type littermates with intra-airway IL-1B or vehicle. Secondarily, we studied human "club cell-like" H322 cells. Results: IL-1B increased whole lung mRNA of secreted (Mucin 5ac, Mucin 5b) and tethered (Mucin 1, Mucin 4) mucins independent of genotype. However, loss of club cell Creb1 increased whole lung mRNA of member RAS oncogene family (Rab3D), decreased mRNA of the muscarinic receptor 3 (M3R) and prevented IL-1B mediated increases in purinergic receptor P2Y, (P2ry2) mRNA. IL-1B increased the density of goblet cells containing neutral mucins in wildtype mice but not in mice with loss of club cell Creb1. These findings suggested that club cell Creb1 regulated mucin secretion. Loss of club cell Creb1 also prevented IL-1B-mediated impairments in airway mechanics. Four days of pharmacologic CREB inhibition in H322 cells increased mRNA abundance of forkhead box A2 (FOXA2), a repressor of goblet cell expansion, and decreased mRNA expression of SAM pointed domain containing ETS transcription factor (SPDEF), a driver of goblet cell expansion. Chromatin immunoprecipitation demonstrated that CREB directly bound to the promoter region of FOXA2, but not to the promoter region of SPDEF. Treatment of H322 cells with IL-1B increased cAMP levels, providing a direct link between IL-1B and CREB signaling. Conclusion: Our findings suggest that club cell Creb1 regulates the pro-mucin properties of IL-1B through pathways likely involving FOXA2.

The Antiviral Activity of the Lectin Griffithsin against SARS-CoV-2 Is Enhanced by the Presence of Structural Proteins

Although COVID-19 transmission has been reduced by the advent of vaccinations and a variety of rapid monitoring techniques, the SARS-CoV-2 virus itself has shown a remarkable ability to mutate and persist. With this long track record of immune escape, researchers are still exploring prophylactic treatments to curtail future SARS-CoV-2 variants. Specifically, much focus has been placed on the antiviral lectin Griffithsin in preventing spike protein-mediated infection via the hACE2 receptor (direct infection). However, an oft-overlooked aspect of SARS-CoV-2 infection is viral capture by attachment receptors such as DC-SIGN, which is thought to facilitate the initial stages of COVID-19 infection in the lung tissue (called trans-infection). In addition, while immune escape is dictated by mutations in the spike protein, coronaviral virions also incorporate M, N, and E structural proteins within the particle. In this paper, we explored how several structural facets of both the SARS-CoV-2 virion and the antiviral lectin Griffithsin can affect and attenuate the infectivity of SARS-CoV-2 pseudovirus. We found that Griffithsin was a better inhibitor of hACE2-mediated direct infection when the coronaviral M protein is present compared to when it is absent (possibly providing an explanation regarding why Griffithsin shows better inhibition against authentic SARS-CoV-2 as opposed to pseudotyped viruses, which generally do not contain M) and that Griffithsin was not an effective inhibitor of DC-SIGN-mediated trans-infection. Furthermore, we found that DC-SIGN appeared to mediate trans-infection exclusively via binding to the SARS-CoV-2 spike protein, with no significant effect observed when other viral proteins (M, N, and/or E) were present. These results provide etiological data that may help to direct the development of novel antiviral treatments, either by leveraging Griffithsin binding to the M protein as a novel strategy to prevent SARS-CoV-2 infection or by narrowing efforts to inhibit trans-infection to focus on DC-SIGN binding to SARS-CoV-2 spike protein.

Genetic analysis of cognitive preservation in the midwestern Amish reveals a novel locus on chromosome 2

INTRODUCTION

Alzheimer disease (AD) remains a debilitating condition with limited treatments and additional therapeutic targets needed. Identifying AD protective genetic loci may identify new targets and accelerate identification of therapeutic treatments. We examined a founder population to identify loci associated with cognitive preservation into advanced age.

METHODS

Genome-wide association and linkage analyses were performed on 946 examined and sampled Amish individuals, aged 76-95, who were either cognitively unimpaired (CU) or impaired (CI).

RESULTS

12 SNPs demonstrated suggestive association (P≤5×10-4) with cognitive preservation. Genetic linkage analyses identified >100 significant (LOD≥3.3) SNPs, some which overlapped with the association results. Only one locus on chromosome 2 retained significance across multiple analyses.

DISCUSSION

A novel significant result for cognitive preservation on chromosome 2 includes the genes LRRTM4 and CTNNA2. Additionally, the lead SNP, rs1402906, impacts the POU3F2 transcription factor binding affinity, which regulates LRRTM4 and CTNNA2.

Spatiotemporal expression patterns of anxiety disorder-associated genes

Anxiety disorders (ADs) are the most common form of mental disorder that affects millions of individuals worldwide. Although physiological studies have revealed the neural circuits related to AD symptoms, how AD-associated genes are spatiotemporally expressed in the human brain still remains unclear. In this study, we integrated genome-wide association studies of four human AD subtypes-generalized anxiety disorder, social anxiety disorder, panic disorder, and obsessive-compulsive disorder-with spatial gene expression patterns. Our investigation uncovered a novel division among AD-associated genes, marked by significant and distinct expression enrichments in the cerebral nuclei, limbic, and midbrain regions. Each gene cluster was associated with specific anxiety-related behaviors, signaling pathways, region-specific gene networks, and cell types. Notably, we observed a significant negative correlation in the temporal expression patterns of these gene clusters during various developmental stages. Moreover, the specific brain regions enriched in each gene group aligned with neural circuits previously associated with negative decision-making and anxious temperament. These results suggest that the two distinct gene clusters may underlie separate neural systems involved in anxiety. As a result, our findings bridge the gap between genes and neural circuitry, shedding light on the mechanisms underlying AD-associated behaviors.

Genes with epigenetic alterations in human pancreatic islets impact mitochondrial function, insulin secretion, and type 2 diabetes

Epigenetic dysregulation may influence disease progression. Here we explore whether epigenetic alterations in human pancreatic islets impact insulin secretion and type 2 diabetes (T2D). In islets, 5,584 DNA methylation sites exhibit alterations in T2D cases versus controls and are associated with HbA1c in individuals not diagnosed with T2D. T2D-associated methylation changes are found in enhancers and regions bound by β-cell-specific transcription factors and associated with reduced expression of e.g. CABLES1, FOXP1, GABRA2, GLR1A, RHOT1, and TBC1D4. We find RHOT1 (MIRO1) to be a key regulator of insulin secretion in human islets. Rhot1-deficiency in β-cells leads to reduced insulin secretion, ATP/ADP ratio, mitochondrial mass, Ca2+, and respiration. Regulators of mitochondrial dynamics and metabolites, including L-proline, glycine, GABA, and carnitines, are altered in Rhot1-deficient β-cells. Islets from diabetic GK rats present Rhot1-deficiency. Finally, RHOT1methylation in blood is associated with future T2D. Together, individuals with T2D exhibit epigenetic alterations linked to mitochondrial dysfunction in pancreatic islets.

Identification of novel genes including NAV2 associated with isolated tall stature

Very tall people attract much attention and represent a clinically and genetically heterogenous group of individuals. Identifying the genetic etiology can provide important insights into the molecular mechanisms regulating linear growth. We studied a three-generation pedigree with five isolated (non-syndromic) tall members and one individual with normal stature by whole exome sequencing; the tallest man had a height of 211 cm. Six heterozygous gene variants predicted as damaging were shared among the four genetically related tall individuals and not present in a family member with normal height. To gain insight into the putative role of these candidate genes in bone growth, we assessed the transcriptome of murine growth plate by microarray and RNA Seq. Two (Ift140, Nav2) of the six genes were well-expressed in the growth plate. Nav2 (p-value 1.91E-62) as well as Ift140 (p-value of 2.98E-06) showed significant downregulation of gene expression between the proliferative and hypertrophic zone, suggesting that these genes may be involved in the regulation of chondrocyte proliferation and/or hypertrophic differentiation. IFT140, NAV2 and SCAF11 have also significantly associated with height in GWAS studies. Pathway and network analysis indicated functional connections between IFT140, NAV2 and SCAF11 and previously associated (tall) stature genes. Knockout of the all-trans retinoic acid responsive gene, neuron navigator 2 NAV2, in Xenopus supports its functional role as a growth promotor. Collectively, our data expand the spectrum of genes with a putative role in tall stature phenotypes and, among other genes, highlight NAV2 as an interesting gene to this phenotype.

Stat5 opposes the transcription factor Tox and rewires exhausted CD8+ T cells toward durable effector-like states during chronic antigen exposure

Rewiring exhausted CD8+ T (Tex) cells toward functional states remains a therapeutic challenge. Tex cells are epigenetically programmed by the transcription factor Tox. However, epigenetic remodeling occurs as Tex cells transition from progenitor (Texprog) to intermediate (Texint) and terminal (Texterm) subsets, suggesting development flexibility. We examined epigenetic transitions between Tex cell subsets and revealed a reciprocally antagonistic circuit between Stat5a and Tox. Stat5 directed Texint cell formation and re-instigated partial effector biology during this Texprog-to-Texint cell transition. Constitutive Stat5a activity antagonized Tox and rewired CD8+ T cells from exhaustion to a durable effector and/or natural killer (NK)-like state with superior anti-tumor potential. Temporal induction of Stat5 activity in Tex cells using an orthogonal IL-2:IL2Rβ-pair fostered Texint cell accumulation, particularly upon PD-L1 blockade. Re-engaging Stat5 also partially reprogrammed the epigenetic landscape of exhaustion and restored polyfunctionality. These data highlight therapeutic opportunities of manipulating the IL-2-Stat5 axis to rewire Tex cells toward more durably protective states.

Single cell multi-omics reveal intra-cell-line heterogeneity across human cancer cell lines

Human cancer cell lines have long served as tools for cancer research and drug discovery, but the presence and the source of intra-cell-line heterogeneity remain elusive. Here, we perform single-cell RNA-sequencing and ATAC-sequencing on 42 and 39 human cell lines, respectively, to illustrate both transcriptomic and epigenetic heterogeneity within individual cell lines. Our data reveal that transcriptomic heterogeneity is frequently observed in cancer cell lines of different tissue origins, often driven by multiple common transcriptional programs. Copy number variation, as well as epigenetic variation and extrachromosomal DNA distribution all contribute to the detected intra-cell-line heterogeneity. Using hypoxia treatment as an example, we demonstrate that transcriptomic heterogeneity could be reshaped by environmental stress. Overall, our study performs single-cell multi-omics of commonly used human cancer cell lines and offers mechanistic insights into the intra-cell-line heterogeneity and its dynamics, which would serve as an important resource for future cancer cell line-based studies.

ACE2 knockout hinders SARS-CoV-2 propagation in iPS cell-derived airway and alveolar epithelial cells

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), the causative agent of COVID-19, continues to spread around the world with serious cases and deaths. It has also been suggested that different genetic variants in the human genome affect both the susceptibility to infection and severity of disease in COVID-19 patients. Angiotensin-converting enzyme 2 (ACE2) has been identified as a cell surface receptor for SARS-CoV and SARS-CoV-2 entry into cells. The construction of an experimental model system using human iPS cells would enable further studies of the association between viral characteristics and genetic variants. Airway and alveolar epithelial cells are cell types of the lung that express high levels of ACE2 and are suitable for in vitro infection experiments. Here, we show that human iPS cell-derived airway and alveolar epithelial cells are highly susceptible to viral infection of SARS-CoV-2. Using gene knockout with CRISPR-Cas9 in human iPS cells we demonstrate that ACE2 plays an essential role in the airway and alveolar epithelial cell entry of SARS-CoV-2 in vitro. Replication of SARS-CoV-2 was strongly suppressed in ACE2 knockout (KO) lung cells. Our model system based on human iPS cell-derived lung cells may be applied to understand the molecular biology regulating viral respiratory infection leading to potential therapeutic developments for COVID-19 and the prevention of future pandemics.

RL2 Enhances the Elimination of Breast Cancer Cells by Doxorubicin

RL2 (recombinant lactaptin 2), a recombinant analogon of the human milk protein Κ-Casein, induces mitophagy and cell death in breast carcinoma cells. Furthermore, RL2 was shown to enhance extrinsic apoptosis upon long-term treatment while inhibiting it upon short-term stimulation. However, the effects of RL2 on the action of chemotherapeutic drugs that induce the intrinsic apoptotic pathway have not been investigated to date. Here, we examined the effects of RL2 on the doxorubicin (DXR)-induced cell death in breast cancer cells with three different backgrounds. In particular, we used BT549 and MDA-MB-231 triple-negative breast cancer (TNBC) cells, T47D estrogen receptor alpha (ERα) positive cells, and SKBR3 human epidermal growth factor receptor 2 (HER2) positive cells. BT549, MDA-MB-231, and T47D cells showed a severe loss of cell viability upon RL2 treatment, accompanied by the induction of mitophagy. Furthermore, BT549, MDA-MB-231, and T47D cells could be sensitized towards DXR treatment with RL2, as evidenced by loss of cell viability. In contrast, SKBR3 cells showed almost no RL2-induced loss of cell viability when treated with RL2 alone, and RL2 did not sensitize SKBR3 cells towards DXR-mediated loss of cell viability. Bioinformatic analysis of gene expression showed an enrichment of genes controlling metabolism in SKBR3 cells compared to the other cell lines. This suggests that the metabolic status of the cells is important for their sensitivity to RL2. Taken together, we have shown that RL2 can enhance the intrinsic apoptotic pathway in TNBC and ERα-positive breast cancer cells, paving the way for the development of novel therapeutic strategies.

The genomic landscape of mammal domestication might be orchestrated by selected transcription factors regulating brain and craniofacial development

Domestication transforms once wild animals into tamed animals that can be then exploited by humans. The process entails modifications in the body, cognition, and behavior that are essentially driven by differences in gene expression patterns. Although genetic and epigenetic mechanisms were shown to underlie such differences, less is known about the role exerted by trans-regulatory molecules, notably transcription factors (TFs) in domestication. In this paper, we conducted extensive in silico analyses aimed to clarify the TF landscape of mammal domestication. We first searched the literature, so as to establish a large list of genes selected with domestication in mammals. From this list, we selected genes experimentally demonstrated to exhibit TF functions. We also considered TFs displaying a statistically significant number of targets among the entire list of (domestication) selected genes. This workflow allowed us to identify 5 candidate TFs (SOX2, KLF4, MITF, NR3C1, NR3C2) that were further assessed in terms of biochemical and functional properties. We found that such TFs-of-interest related to mammal domestication are all significantly involved in the development of the brain and the craniofacial region, as well as the immune response and lipid metabolism. A ranking strategy, essentially based on a survey of protein-protein interactions datasets, allowed us to identify SOX2 as the main candidate TF involved in domestication-associated evolutionary changes. These findings should help to clarify the molecular mechanics of domestication and are of interest for future studies aimed to understand the behavioral and cognitive changes associated to domestication.

Enzymatic vitamin A2 production enables red-shifted optogenetics

Optogenetics is a technology using light-sensitive proteins to control signaling pathways and physiological processes in cells and organs and has been applied in neuroscience, cardiovascular sciences, and many other research fields. Most commonly used optogenetic actuators are sensitive to blue and green light, but red-light activation would allow better tissue penetration and less phototoxicity. Cyp27c1 is a recently deorphanized cytochrome P450 enzyme that converts vitamin A1 to vitamin A2, thereby red-shifting the spectral sensitivity of visual pigments and enabling near-infrared vision in some aquatic species.Here, we investigated the ability of Cyp27c1-generated vitamin A2 to induce a shift in spectral sensitivity of the light-gated ion channel Channelrhodopsin-2 (ChR2) and its red-shifted homolog ReaChR. We used patch clamp to measure photocurrents at specific wavelengths in HEK 293 cells expressing ChR2 or ReaChR. Vitamin A2 incubation red-shifted the wavelength for half-maximal currents (λ50%) by 6.8 nm for ChR2 and 12.4 nm for ReaChR. Overexpression of Cyp27c1 in HEK 293 cells showed mitochondrial localization, and HPLC analysis showed conversion of vitamin A1 to vitamin A2. Notably, the λ50% of ChR2 photocurrents was red-shifted by 10.5 nm, and normalized photocurrents at 550 nm were about twofold larger with Cyp27c1 expression. Similarly, Cyp27c1 shifted the λ50% of ReaChR photocurrents by 14.3 nm and increased normalized photocurrents at 650 nm almost threefold.Since vitamin A2 incubation is not a realistic option for in vivo applications and expression of Cyp27c1 leads to a greater red-shift in spectral sensitivity, we propose co-expression of this enzyme as a novel strategy for red-shifted optogenetics.

Integrative Analysis of PAIP2B to Identify a Novel Biomarker for Pancreatic Ductal Adenocarcinoma

The aim of the study was to evaluate the potential diagnostic and prognostic value of gene, Poly A-Binding Protein Interacting Protein 2B ( PAIP2B ) in pancreatic cancer. We used the gene expression data and clinical information of pancreatic adenocarcinoma patients from The Cancer Genome Atlas database and Gene Expression Omnibus database to analyze the expression of PAIP2B in pancreatic cancer samples, and validated the expression of PAIP2B in tumor tissue, using bioinformatics technology to explore the prognostic value of PAIP2B and its possible biological function. A significantly lower level of PAIP2B was observed in pancreatic cancer patients than in controls, and validated by immunohistochemistry. PAIP2B reduced the proliferation and invasion of cancer cells and had a significantly high expression in early stage. Patients with lower levels of PAIP2B had a significantly shorter median survival time than those with higher levels. DNA demethylation played an important role in PAIP2B expression. In addition, PAIP2B expression was significantly associated with the tumor-infiltrating immune cells, especially T cells CD8, T cells CD4 memory resting, macrophages M0, and dendritic cells resting. Our study also found that PAIP2B regulated miRNA function leading to disease progression in pancreatic cancer patients. Our study explored the potential value of PAIP2B as a biological link between prognosis and pancreatic cancer, and provided reference for the follow-up study on the role of PAIP2B in pancreatic cancer.

Hepatocyte nuclear factor 4α mediated quinolinate phosphoribosylltransferase (QPRT) expression in the kidney facilitates resilience against acute kidney injury

Nicotinamide adenine dinucleotide (NAD+) levels decline in experimental models of acute kidney injury (AKI). Attenuated enzymatic conversion of tryptophan to NAD+ in tubular epithelium may contribute to adverse cellular and physiological outcomes. Mechanisms underlying defense of tryptophan-dependent NAD+ production are incompletely understood. Here we show that regulation of a bottleneck enzyme in this pathway, quinolinate phosphoribosyltransferase (QPRT) may contribute to kidney resilience. Expression of QPRT declined in two unrelated models of AKI. Haploinsufficient mice developed worse outcomes compared to littermate controls whereas novel, conditional gain-of-function mice were protected from injury. Applying these findings, we then identified hepatocyte nuclear factor 4 alpha (HNF4α) as a candidate transcription factor regulating QPRT expression downstream of the mitochondrial biogenesis regulator and NAD+ biosynthesis inducer PPARgamma coactivator-1-alpha (PGC1α). This was verified by chromatin immunoprecipitation. A PGC1α - HNF4α -QPRT axis controlled NAD+ levels across cellular compartments and modulated cellular ATP. These results propose that tryptophan-dependent NAD+ biosynthesis via QPRT and induced by HNF4α may be a critical determinant of kidney resilience to noxious stressors.

Comparative analysis identifies genetic and molecular factors associated with prognostic clusters of PANoptosis in glioma, kidney and melanoma cancer

The importance of inflammatory cell death, PANoptosis, in cancer is increasingly being recognized. PANoptosis can promote or inhibit tumorigenesis in context-dependent manners, and a computational approach leveraging transcriptomic profiling of genes involved in PANoptosis has shown that patients can be stratified into PANoptosis High and PANoptosis Low clusters that have significant differences in overall survival for low grade glioma (LGG), kidney renal cell carcinoma (KIRC) and skin cutaneous melanoma (SKCM). However, the molecular mechanisms that contribute to differential prognosis between PANoptosis clusters require further elucidation. Therefore, we performed a comprehensive comparison of genetic, genomic, tumor microenvironment, and pathway characteristics between the PANoptosis High and PANoptosis Low clusters to determine the relevance of each component in driving the differential associations with prognosis for LGG, KIRC and SKCM. Across these cancer types, we found that activation of the proliferation pathway was significantly different between PANoptosis High and Low clusters. In LGG and SKCM, we also found that aneuploidy and immune cell densities and activations contributed to differences in PANoptosis clusters. In individual cancers, we identified important roles for barrier gene pathway activation (in SKCM) and the somatic mutation profiles of driver oncogenes as well as hedgehog signaling pathway activation (in LGG). By identifying these genetic and molecular factors, we can possibly improve the prognosis for at risk-stratified patient populations based on the PANoptosis phenotype in LGG, KIRC and SKCM. This not only advances our mechanistic understanding of cancer but will allow for the selection of optimal treatment strategies.

Mapping cardiac remodeling in chronic kidney disease

Patients with advanced chronic kidney disease (CKD) mostly die from sudden cardiac death and recurrent heart failure. The mechanisms of cardiac remodeling are largely unclear. To dissect molecular and cellular mechanisms of cardiac remodeling in CKD in an unbiased fashion, we performed left ventricular single-nuclear RNA sequencing in two mouse models of CKD. Our data showed a hypertrophic response trajectory of cardiomyocytes with stress signaling and metabolic changes driven by soluble uremia-related factors. We mapped fibroblast to myofibroblast differentiation in this process and identified notable changes in the cardiac vasculature, suggesting inflammation and dysfunction. An integrated analysis of cardiac cellular responses to uremic toxins pointed toward endothelin-1 and methylglyoxal being involved in capillary dysfunction and TNFα driving cardiomyocyte hypertrophy in CKD, which was validated in vitro and in vivo. TNFα inhibition in vivo ameliorated the cardiac phenotype in CKD. Thus, interventional approaches directed against uremic toxins, such as TNFα, hold promise to ameliorate cardiac remodeling in CKD.

Mouse models of pediatric high-grade gliomas with MYCN amplification reveal intratumoral heterogeneity and lineage signatures

Pediatric high-grade gliomas of the subclass MYCN (HGG-MYCN) are highly aggressive tumors frequently carrying MYCN amplifications, TP53 mutations, or both alterations. Due to their rarity, such tumors have only recently been identified as a distinct entity, and biological as well as clinical characteristics have not been addressed specifically. To gain insights into tumorigenesis and molecular profiles of these tumors, and to ultimately suggest alternative treatment options, we generated a genetically engineered mouse model by breeding hGFAP-cre::Trp53Fl/Fl::lsl-MYCN mice. All mice developed aggressive forebrain tumors early in their lifetime that mimic human HGG-MYCN regarding histology, DNA methylation, and gene expression. Single-cell RNA sequencing revealed a high intratumoral heterogeneity with neuronal and oligodendroglial lineage signatures. High-throughput drug screening using both mouse and human tumor cells finally indicated high efficacy of Doxorubicin, Irinotecan, and Etoposide as possible therapy options that children with HGG-MYCN might benefit from.

Construction of an anoikis-associated lncRNA-miRNA-mRNA network reveals the prognostic role of β-elemene in non-small cell lung cancer

β-Elemene is the main active ingredient in Curcumae Rhizoma that exerts antitumour effects. Anoikis affects tumour development through various biological pathways in non-small cell lung cancer (NSCLC), but the regulation between β-elemene and anoikis remains to be explored. First, we explored the molecular expression patterns of anoikis-associated genes (AAGs) using consensus clustering and characterized the impact of AAGs on patient prognosis, clinical characteristics, and genomic instability. In addition, we revealed that AAG regulatory genes have rich interactions with β-elemene targets, and established a lncRNA-miRNA-mRNA network to explain the effect of β-elemene on anoikis. Finally, to reveal the prognostic effect of their correlation, the prognostic scoring model and clinical nomogram of β-elemene and anoikis were successfully established by least absolute shrinkage and selection operator (LASSO) and random forest algorithms. This prognostic scoring model containing noncoding RNA (ncRNA) can indicate the immunotherapy and mutational landscape, providing a novel theoretical basis and direction for the study of the antitumour mechanism of β-elemene in NSCLC patients.

Cross-talk between PARN and EGFR-STAT3 Signaling Facilitates Self-Renewal and Proliferation of Glioblastoma Stem Cells

Glioblastoma is the most common type of malignant primary brain tumor and displays highly aggressive and heterogeneous phenotypes. The transcription factor STAT3 has been reported to play a key role in glioblastoma malignancy. Thus, discovering targets and functional downstream networks regulated by STAT3 that govern glioblastoma pathogenesis may lead to improved treatment strategies. In this study, we identified that poly(A)-specific ribonuclease (PARN), a key modulator of RNA metabolism, activates EGFR-STAT3 signaling to support glioblastoma stem cells (GSC). Functional integrative analysis of STAT3 found PARN as the top-scoring transcriptional target involved in RNA processing in patients with glioblastoma, and PARN expression was strongly correlated with poor patient survival and elevated malignancy. PARN positively regulated self-renewal and proliferation of GSCs through its 3'-5' exoribonuclease activity. EGFR was identified as a clinically relevant target of PARN in GSCs. PARN positively modulated EGFR by negatively regulating the EGFR-targeting miRNA miR-7, and increased EGFR expression created a positive feedback loop to increase STAT3 activation. PARN depletion in GSCs reduced infiltration and prolonged survival in orthotopic brain tumor xenografts; similar results were observed using siRNA nanocapsule-mediated PARN targeting. Pharmacological targeting of STAT3 also confirmed PARN regulation by STAT3 signaling. In sum, these results suggest that a STAT3-PARN regulatory network plays a pivotal role in tumor progression and thus may represent a target for glioblastoma therapeutics.

SIGNIFICANCE

A positive feedback loop comprising PARN and EGFR-STAT3 signaling supports self-renewal and proliferation of glioblastoma stem cells to drive tumor progression and can be targeted in glioblastoma therapeutics.

STAMarker: determining spatial domain-specific variable genes with saliency maps in deep learning

Spatial transcriptomics characterizes gene expression profiles while retaining the information of the spatial context, providing an unprecedented opportunity to understand cellular systems. One of the essential tasks in such data analysis is to determine spatially variable genes (SVGs), which demonstrate spatial expression patterns. Existing methods only consider genes individually and fail to model the inter-dependence of genes. To this end, we present an analytic tool STAMarker for robustly determining spatial domain-specific SVGs with saliency maps in deep learning. STAMarker is a three-stage ensemble framework consisting of graph-attention autoencoders, multilayer perceptron (MLP) classifiers, and saliency map computation by the backpropagated gradient. We illustrate the effectiveness of STAMarker and compare it with serveral commonly used competing methods on various spatial transcriptomic data generated by different platforms. STAMarker considers all genes at once and is more robust when the dataset is very sparse. STAMarker could identify spatial domain-specific SVGs for characterizing spatial domains and enable in-depth analysis of the region of interest in the tissue section.

Identification of a Novel Oxidative Stress- and Anoikis-Related Prognostic Signature and Its Immune Landscape Analysis in Non-Small Cell Lung Cancer

The objective of this study was to identify a kind of prognostic signature based on oxidative stress- and anoikis-related genes (OARGs) for predicting the prognosis and immune landscape of NSCLC. Initially, We identified 47 differentially expressed OARGs that primarily regulate oxidative stress and epithelial cell infiltration through the PI3K-Akt pathway. Subsequently, 10 OARGs related to prognosis determined two potential clusters. A cluster was associated with a shorter survival level, lower immune infiltration, higher stemness index and tumor mutation burden. Next, The best risk score model constructed by prognostic OARGs was the Random Survival Forest model, and it included SLC2A1, LDHA and PLAU. The high-risk group was associated with cluster A and poor prognosis, with a higher tumor mutation burden, stemness index and proportion of M0-type macrophages, and a lower immune checkpoint expression level, immune function score and IPS score. The calibration curve and decision-making curve showed that the risk score combined with clinical pathological characteristics could be used to construct a nomogram for guiding the clinical treatment strategies. Finally, We found that all three hub genes were highly expressed in tumor tissues, and LDHA expression was mainly regulated by has-miR-338-3p, has-miR-330-5p and has-miR-34c-5p. Altogether, We constructed an OARG-related prognostic signature to reveal potential relationships between the signature and clinical characteristics, TME, stemness, tumor mutational burden, drug sensitivity and immune landscape in NSCLC patients.

DNA methylation sites in early adulthood characterised by pubertal timing and development: a twin study

BACKGROUND

Puberty is a highly heritable and variable trait, with environmental factors having a role in its eventual timing and development. Early and late pubertal onset are both associated with various diseases developing later in life, and epigenetic characterisation of pubertal timing and development could lead to important insights. Blood DNA methylation, reacting to both genotype and environment, has been associated with puberty; however, such studies are relatively scarce. We investigated peripheral blood DNA methylation profiles (using Illumina 450 K and EPIC platforms) of 1539 young adult Finnish twins associated with pubertal development scale (PDS) at ages 12 and 14 as well as pubertal age (PA).

RESULTS

Fixed effect meta-analysis of the two platforms on 347,521 CpGs in common identified 58 CpG sites associated (p < 1 × 10-5) with either PDS or PA. All four CpGs associated with PA and 45 CpGs associated with PDS were sex-specific. Thirteen CpGs had a high heritability (h2: 0.51-0.98), while one CpG site (mapped to GET4) had a high shared environmental component accounting for 68% of the overall variance in methylation at the site. Utilising twin discordance analysis, we found 6 CpG sites (5 associated with PDS and 1 with PA) that had an environmentally driven association with puberty. Furthermore, genes with PDS- or PA-associated CpGs were consistently linked to various developmental processes and diseases such as breast, prostate and ovarian cancer, while methylation quantitative trait loci of associated CpG sites were enriched in immune pathways developing during puberty.

CONCLUSIONS

By identifying puberty-associated DNA methylation sites and examining the effects of sex, environment and genetics, we shed light on the intricate interplay between environment and genetics in the context of puberty. Through our comprehensive analysis, we not only deepen the understanding of the significance of both genetic and environmental factors in the complex processes of puberty and its timing, but also gain insights into potential links with disease risks.

M2 macrophage-derived cathepsin S promotes peripheral nerve regeneration via fibroblast-Schwann cell-signaling relay

BACKGROUND

Although peripheral nerves have an intrinsic self-repair capacity following damage, functional recovery is limited in patients. It is a well-established fact that macrophages accumulate at the site of injury. Numerous studies indicate that the phenotypic shift from M1 macrophage to M2 macrophage plays a crucial role in the process of axon regeneration. This polarity change is observed exclusively in peripheral macrophages but not in microglia and CNS macrophages. However, the molecular basis of axonal regeneration by M2 macrophage is not yet fully understood. Herein, we aimed to identify the M2 macrophage-derived axon regeneration factor.

METHODS

We established a peripheral nerve injury model by transection of the inferior alveolar nerve (IANX) in Sprague-Dawley rats. Transcriptome analysis was performed on the injured nerve. Recovery from sensory deficits in the mandibular region and histological reconnection of IAN after IANX were assessed in rats with macrophage depletion by clodronate. We investigated the effects of adoptive transfer of M2 macrophages or M2-derived cathepsin S (CTSS) on the sensory deficit. CTSS initiating signaling was explored by western blot analysis in IANX rats and immunohistochemistry in co-culture of primary fibroblasts and Schwann cells (SCs).

RESULTS

Transcriptome analysis revealed that CTSS, a macrophage-selective lysosomal protease, was upregulated in the IAN after its injury. Spontaneous but partial recovery from a sensory deficit in the mandibular region after IANX was abrogated by macrophage ablation at the injured site. In addition, a robust induction of c-Jun, a marker of the repair-supportive phenotype of SCs, after IANX was abolished by macrophage ablation. As in transcriptome analysis, CTSS was upregulated at the injured IAN than in the intact IAN. Endogenous recovery from hypoesthesia was facilitated by supplementation of CTSS but delayed by pharmacological inhibition or genetic silencing of CTSS at the injured site. Adoptive transfer of M2-polarized macrophages at this site facilitated sensory recovery dependent on CTSS in macrophages. Post-IANX, CTSS caused the cleavage of Ephrin-B2 in fibroblasts, which, in turn, bound EphB2 in SCs. CTSS-induced Ephrin-B2 cleavage was also observed in human sensory nerves. Inhibition of CTSS-induced Ephrin-B2 signaling suppressed c-Jun induction in SCs and sensory recovery.

CONCLUSIONS

These results suggest that M2 macrophage-derived CTSS contributes to axon regeneration by activating SCs via Ephrin-B2 shedding from fibroblasts.

The RioK1 network determines p53 activity at multiple levels

By responding to a host of adverse conditions, ranging from DNA damage to viral infection, transcription factor p53 supports genomic stability, cellular health, and survival. Not surprisingly, tumours across the cancer spectrum carry mutations in p53, misexpress the protein, or dysregulate its activity. Several signalling pathways, many of which comprise oncogenic proteins, converge upon p53 to control its stability and activity. We here present the conserved kinase/ATPase RioK1 as an upstream factor that determines p53 activity at the DNA, RNA, and protein levels. It achieves this task by integrating the regulatory events that act on p53 into a coherent response circuit. We will also discuss how RIOK1 overexpression represents an alternative mechanism for cancers to inactivate p53, and how targeting RioK1 could eradicate malignancies that are driven by a dysregulated RioK1-p53 network.

Genomic inbreeding estimation, runs of homozygosity, and heterozygosity-enriched regions uncover signals of selection in the Quarter Horse racing line

With the advent of genomics, significant progress has been made in the genetic improvement of livestock species, particularly through increased accuracy in predicting breeding values for selecting superior animals and the possibility of performing a high-resolution genetic scan throughout the genome of an individual. The main objectives of this study were to estimate the individual genomic inbreeding coefficient based on runs of homozygosity (FROH ), to identify and characterize runs of homozygosity and heterozygosity (ROH and ROHet, respectively; length and distribution) throughout the genome, and to map selection signatures in relevant chromosomal regions in the Quarter Horse racing line. A total of 336 animals registered with the Brazilian Association of Quarter Horse Breeders (ABQM) were genotyped. One hundred and twelve animals were genotyped using the Equine SNP50 BeadChip (Illumina, USA), with 54,602 single nucleotide polymorphisms (SNPs; 54K). The remaining 224 samples were genotyped using the Equine SNP70 BeadChip (Illumina, USA) with 65,157 SNPs (65K). To ensure data quality, we excluded animals with a call rate below 0.9. We also excluded SNPs located on non-autosomal chromosomes, as well as those with a call rate below 0.9 or a p-value below 1 × 10-5 for Hardy-Weinberg equilibrium. The results indicate moderate to high genomic inbreeding, with 46,594 ROH and 16,101 ROHet detected. In total, 30 and 14 candidate genes overlap with ROH and ROHet regions, respectively. The ROH islands showed genes linked to crucial biological processes, such as cell differentiation (CTBP1, WNT5B, and TMEM120B), regulation of glucose metabolic process (MAEA and NKX1-1), heme transport (PGRMC2), and negative regulation of calcium ion import (VDAC1). In ROHet, the islands showed genes related to respiratory capacity (OR7D19, OR7D4G, OR7D4E, and OR7D4J) and muscle repair (EGFR and BCL9). These findings could aid in selecting animals with greater regenerative capacity and developing treatments for muscle disorders in the QH breed. This study serves as a foundation for future research on equine breeds. It can contribute to developing reproductive strategies in animal breeding programs to improve and preserve the Quarter Horse breed.

Intracellular fraction of zona pellucida protein 3 is required for the oocyte-to-embryo transition in mice

In oocyte biology, the zona pellucida has long been known to operate three extracellular functions downstream of the secretory pathway, namely, encasing the oocytes in ovarian follicles, mediating sperm-oocyte interaction, and preventing premature embryo contact with oviductal epithelium. The present study uncovers a fourth function that is fundamentally distinct from the other three, being critical for embryonic cell survival in mice. Intriguingly, the three proteins of the mouse zona pellucida (ZP1, ZP2, ZP3) were found abundantly present also inside the embryo 4 days after fertilization, as shown by mass spectrometry, immunoblotting, and immunofluorescence. Contrary to current understanding of the roles of ZP proteins, ZP3 was associated more with the cytoskeleton than with secretory vesicles in the subcortical region of metaphase II oocytes and zygotes, and was excluded from regions of cell-cell contact in cleavage-stage embryos. Trim-away-mediated knockdown of ZP3 in fertilized oocytes hampered the first zygotic cleavage, while ZP3 overexpression supported blastocyst formation. Transcriptome analysis of ZP3-knockdown embryos pointed at defects of cytoplasmic translation in the context of embryonic genome activation. This conclusion was supported by reduced protein synthesis in the ZP3-knockdown and by the lack of cleavage arrest when Trim-away was postponed from the one-cell to the late two-cell stage. These data place constraints on the notion that zona proteins only operate in the extracellular space, revealing also a role during the oocyte-to-embryo transition. Ultimately, these data recruit ZP3 into the family of maternal factors that contribute to developmental competence of mouse oocytes.

angptl4 gene expression as a marker of adaptive homeostatic response to social isolation across the lifespan in zebrafish

Social isolation has detrimental health effects, but the underlying mechanisms are unclear. Here, we investigated the impact of 2 weeks of isolation on behavior and gene expression in the central nervous system at different life stages of zebrafish. Results showed that socially deprived young adult zebrafish experienced increased anxiety, accompanied by changes in gene expression. Most gene expression patterns returned to normal within 24 hours of reintroduction to a social environment, except angptl4, which was upregulated after reintroduction, suggesting an adaptive mechanism. Similarly, aging zebrafish displayed heightened anxiety and increased central nervous system expression of angptl4 during isolation, but effects were reversed upon reintroduction to a social group. The findings imply that angptl4 plays a homeostatic role in response to social isolation, which varies across the lifespan. The study emphasizes the importance of social interactions for psychological well-being and highlights the negative consequences of isolation, especially in older individuals. Further research may unravel how social isolation affects angptl4 expression and its developmental and aging effects.

scapGNN: A graph neural network-based framework for active pathway and gene module inference from single-cell multi-omics data

Although advances in single-cell technologies have enabled the characterization of multiple omics profiles in individual cells, extracting functional and mechanistic insights from such information remains a major challenge. Here, we present scapGNN, a graph neural network (GNN)-based framework that creatively transforms sparse single-cell profile data into the stable gene-cell association network for inferring single-cell pathway activity scores and identifying cell phenotype-associated gene modules from single-cell multi-omics data. Systematic benchmarking demonstrated that scapGNN was more accurate, robust, and scalable than state-of-the-art methods in various downstream single-cell analyses such as cell denoising, batch effect removal, cell clustering, cell trajectory inference, and pathway or gene module identification. scapGNN was developed as a systematic R package that can be flexibly extended and enhanced for existing analysis processes. It provides a new analytical platform for studying single cells at the pathway and network levels.