WebCSEA: web-based cell-type-specific enrichment analysis of genes

Comparison of stromal vascular fraction cell composition between Coleman fat and extracellular matrix/stromal vascular fraction gel

As a mechanically condensed product of Coleman fat, extracellular matrix/stromal vascular fraction gel (ECM/SVF-gel) eliminates adipocytes, concentrates SVF cells, and improves fat graft retention. This study aims to compare SVF cell composition between Coleman fat and ECM/SVF-gel. Matched Coleman fat and ECM/SVF-gel of 28 healthy women were subjected to RNA-seq, followed by functional enrichment and cell-type-specific enrichment analyses, and deconvolution of SVF cell subsets, reconstructing SVF cell composition in the transcriptome level. ECM/SVF-gels had 9 upregulated and 73 downregulated differentially expressed genes (DEGs). Downregulated DEGs were mainly associated with inflammatory and immune responses, and enriched in fat macrophages. M2 macrophages, resting CD4+ memory T cells, M1 macrophages, resting mast cells, and M0 macrophages ranked in the top five most prevalent immune cells in the two groups. The proportions of the principal non-immune cells (e.g., adipose-derived stem cells, pericytes, preadipocytes, microvascular endothelial cells) had no statistical differences between the two groups. Our findings reveal ECM/SVF-gels share the same dominant immune cells beneficial to fat graft survival with Coleman fat, but exhibiting obvious losses of immune cells (especially macrophages), while non-immune cells necessary for adipose regeneration might have no significant loss in ECM/SVF-gels and their biological effects could be markedly enhanced by the ECM/SVF-gel's condensed nature.

High-throughput identification of repurposable neuroactive drugs with potent anti-glioblastoma activity

Glioblastoma, the most aggressive primary brain cancer, has a dismal prognosis, yet systemic treatment is limited to DNA-alkylating chemotherapies. New therapeutic strategies may emerge from exploring neurodevelopmental and neurophysiological vulnerabilities of glioblastoma. To this end, we systematically screened repurposable neuroactive drugs in glioblastoma patient surgery material using a clinically concordant and single-cell resolved platform. Profiling more than 2,500 ex vivo drug responses across 27 patients and 132 drugs identified class-diverse neuroactive drugs with potent anti-glioblastoma efficacy that were validated across model systems. Interpretable molecular machine learning of drug-target networks revealed neuroactive convergence on AP-1/BTG-driven glioblastoma suppression, enabling expanded in silico screening of more than 1 million compounds with high patient validation accuracy. Deep multimodal profiling confirmed Ca2+-driven AP-1/BTG-pathway induction as a neuro-oncological glioblastoma vulnerability, epitomized by the anti-depressant vortioxetine synergizing with current standard-of-care chemotherapies in vivo. These findings establish an actionable framework for glioblastoma treatment rooted in its neural etiology.

Shared Genetic Architecture Among Gastrointestinal Diseases, Schizophrenia, and Brain Subcortical Volumes

BACKGROUND AND HYPOTHESIS

The gut-brain axis plays important roles in both gastrointestinal diseases (GI diseases) and schizophrenia (SCZ). Moreover, both GI diseases and SCZ exhibit notable abnormalities in brain subcortical volumes. However, the genetic mechanisms underlying the comorbidity of these diseases and the shared alterations in brain subcortical volumes remain unclear.

STUDY DESIGN

Using the genome-wide association studies data of SCZ, 14 brain subcortical volumes, and 8 GI diseases, the global polygenic overlap and local genetic correlations were identified, as well as the shared genetic variants among those phenotypes. Furthermore, we conducted multi-trait colocalization analyses to bolster our findings. Functional annotations, cell-type enrichment, and protein-protein interaction (PPI) analyses were carried out to reveal the critical etiology and pathology mechanisms.

STUDY RESULTS

The global polygenic overlap and local genetic correlations informed the close relationships between SCZ and both GI diseases and brain subcortical volumes. Moreover, 84 unique lead-shared variants were identified. The associated genes were linked to vital biological processes within the immune system. Additionally, significant correlations were observed with key immune cells and the PPI analysis identified several histone-associated hub genes. These findings highlighted the pivotal roles played by the immune system for both SCZ and GI diseases, along with the shared alterations in brain subcortical volumes.

CONCLUSIONS

These findings revealed the shared genetic architecture contributing to SCZ and GI diseases, as well as their shared alterations in brain subcortical volumes. These insights have substantial implications for the concurrent development of intervention and therapy targets for these diseases.

TransTEx: novel tissue-specificity scoring method for grouping human transcriptome into different expression groups

MOTIVATION

Although human tissues carry out common molecular processes, gene expression patterns can distinguish different tissues. Traditional informatics methods, primarily at the gene level, overlook the complexity of alternative transcript variants and protein isoforms produced by most genes, changes in which are linked to disease prognosis and drug resistance.

RESULTS

We developed TransTEx (Transcript-level Tissue Expression), a novel tissue-specificity scoring method, for grouping transcripts into four expression groups. TransTEx applies sequential cut-offs to tissue-wise transcript probability estimates, subsampling-based P-values and fold-change estimates. Application of TransTEx on GTEx mRNA-seq data divided 199 166 human transcripts into different groups as 17 999 tissue-specific (TSp), 7436 tissue-enhanced, 36 783 widely expressed (Wide), 79 191 lowly expressed (Low), and 57 757 no expression (Null) transcripts. Testis has the most (13 466) TSp isoforms followed by liver (890), brain (701), pituitary (435), and muscle (420). We found that the tissue specificity of alternative transcripts of a gene is predominantly influenced by alternate promoter usage. By overlapping brain-specific transcripts with the cell-type gene-markers in scBrainMap database, we found that 63% of the brain-specific transcripts were enriched in nonneuronal cell types, predominantly astrocytes followed by endothelial cells and oligodendrocytes. In addition, we found 61 brain cell-type marker genes encoding a total of 176 alternative transcripts as brain-specific and 22 alternative transcripts as testis-specific, highlighting the complex TSp and cell-type specific gene regulation and expression at isoform-level. TransTEx can be adopted to the analysis of bulk RNA-seq or scRNA-seq datasets to find tissue- and/or cell-type specific isoform-level gene markers.

AVAILABILITY AND IMPLEMENTATION

TransTEx database: https://bmi.cewit.stonybrook.edu/transtexdb/ and the R package is available via GitHub: https://github.com/pallavisurana1/TransTEx.

Genetic overlap between Alzheimer's disease and immune-mediated diseases: an atlas of shared genetic determinants and biological convergence

The occurrence of immune disease comorbidities in Alzheimer's disease (AD) has been observed in both epidemiological and molecular studies, suggesting a neuroinflammatory basis in AD. However, their shared genetic components have not been systematically studied. Here, we composed an atlas of the shared genetic associations between 11 immune-mediated diseases and AD by analyzing genome-wide association studies (GWAS) summary statistics. Our results unveiled a significant genetic overlap between AD and 11 individual immune-mediated diseases despite negligible genetic correlations, suggesting a complex shared genetic architecture distributed across the genome. The shared loci between AD and immune-mediated diseases implicated several genes, including GRAMD1B, FUT2, ADAMTS4, HBEGF, WNT3, TSPAN14, DHODH, ABCB9, and TNIP1, all of which are protein-coding genes and thus potential drug targets. Top biological pathways enriched with these identified shared genes were related to the immune system and cell adhesion. In addition, in silico single-cell analyses showed enrichment of immune and brain cells, including neurons and microglia. In summary, our results suggest a genetic relationship between AD and the 11 immune-mediated diseases, pinpointing the existence of a shared however non-causal genetic basis. These identified protein-coding genes have the potential to serve as a novel path to therapeutic interventions for both AD and immune-mediated diseases and their comorbidities.

DeepFace: Deep-learning-based framework to contextualize orofacial-cleft-related variants during human embryonic craniofacial development

Orofacial clefts (OFCs) are among the most common human congenital birth defects. Previous multiethnic studies have identified dozens of associated loci for both cleft lip with or without cleft palate (CL/P) and cleft palate alone (CP). Although several nearby genes have been highlighted, the "casual" variants are largely unknown. Here, we developed DeepFace, a convolutional neural network model, to assess the functional impact of variants by SNP activity difference (SAD) scores. The DeepFace model is trained with 204 epigenomic assays from crucial human embryonic craniofacial developmental stages of post-conception week (pcw) 4 to pcw 10. The Pearson correlation coefficient between the predicted and actual values for 12 epigenetic features achieved a median range of 0.50-0.83. Specifically, our model revealed that SNPs significantly associated with OFCs tended to exhibit higher SAD scores across various variant categories compared to less related groups, indicating a context-specific impact of OFC-related SNPs. Notably, we identified six SNPs with a significant linear relationship to SAD scores throughout developmental progression, suggesting that these SNPs could play a temporal regulatory role. Furthermore, our cell-type specificity analysis pinpointed the trophoblast cell as having the highest enrichment of risk signals associated with OFCs. Overall, DeepFace can harness distal regulatory signals from extensive epigenomic assays, offering new perspectives for prioritizing OFC variants using contextualized functional genomic features. We expect DeepFace to be instrumental in accessing and predicting the regulatory roles of variants associated with OFCs, and the model can be extended to study other complex diseases or traits.

Investigating gene functions and single-cell expression profiles of de novo variants in orofacial clefts

Orofacial clefts (OFCs) are common congenital birth defects with various etiologies, including genetic variants. Online Mendelian Inheritance in Man (OMIM) annotated several hundred genes involving OFCs. Furthermore, several hundreds of de novo variants (DNVs) have been identified from individuals with OFCs. Some DNVs are related to known OFC genes or pathways, but there are still many DNVs whose relevance to OFC development is unknown. To explore novel gene functions and their cellular expression profiles, we focused on DNVs in genes that were not listed in OMIM. We collected 960 DNVs in 853 genes from published studies and curated these genes, based on the DNVs' deleteriousness, into 230 and 23 genes related to cleft lip with or without cleft palate (CL/P) and cleft palate only (CPO), respectively. For comparison, we curated 178 CL/P and 277 CPO genes from OMIM. In CL/P, the pathways enriched in DNV and OMIM genes were significantly overlapped (p = 0.002). Single-cell RNA sequencing (scRNA-seq) analysis of mouse lip development revealed that both gene sets had abundant expression in the ectoderm (DNV genes: adjusted p = 0.032, OMIM genes: adjusted p < 0.0002), while only DNV genes were enriched in the endothelium (adjusted p = 0.032). Although we did not achieve significant findings using CPO gene sets, which was mainly due to the limited number of DNV genes, scRNA-seq analysis implicated various expression patterns among DNV and OMIM genes. Our results suggest that combinatory pathway and scRNA-seq data analyses are helpful for contextualizing genes in OFC development.

Abdominal aortic aneurysm and cardiometabolic traits share strong genetic susceptibility to lipid metabolism and inflammation

Abdominal aortic aneurysm has a high heritability and often co-occurs with other cardiometabolic disorders, suggesting shared genetic susceptibility. We investigate this commonality leveraging recent GWAS studies of abdominal aortic aneurysm and 32 cardiometabolic traits. We find significant genetic correlations between abdominal aortic aneurysm and 21 of the cardiometabolic traits investigated, including causal relationships with coronary artery disease, hypertension, lipid traits, and blood pressure. For each trait pair, we identify shared causal variants, genes, and pathways, revealing that cholesterol metabolism and inflammation are shared most prominently. Additionally, we show the tissue and cell type specificity in the shared signals, with strong enrichment across traits in the liver, arteries, adipose tissues, macrophages, adipocytes, and fibroblasts. Finally, we leverage drug-gene databases to identify several lipid-lowering drugs and antioxidants with high potential to treat abdominal aortic aneurysm with comorbidities. Our study provides insight into the shared genetic mechanism between abdominal aortic aneurysm and cardiometabolic traits, and identifies potential targets for pharmacological intervention.

Genome-wide DNA methylation and transcriptomic analysis of liver tissues subjected to early ischemia/reperfusion injury upon human liver transplantation

INTRODUCTION AND OBJECTIVES

Epigenetic changes represent a mechanism connecting external stresses with long-term modifications of gene expression programs. In solid organ transplantation, ischemia-reperfusion injury (IRI) appears to induce epigenomic changes in the graft, although the currently available data are extremely limited. The present study aimed to characterize variations in DNA methylation and their effects on the transcriptome in liver transplantation from brain-dead donors.

PATIENTS AND METHODS

12 liver grafts were evaluated through serial biopsies at different timings in the procurement-transplantation process: T0 (warm procurement, in donor), T1 (bench surgery), and T2 (after reperfusion, in recipient). DNA methylation (DNAm) and transcriptome profiles of biopsies were analyzed using microarrays and RNAseq.

RESULTS

Significant variations in DNAm were identified, particularly between T2 and T0. Functional enrichment of the best 1000 ranked differentially methylated promoters demonstrated that 387 hypermethylated and 613 hypomethylated promoters were involved in spliceosomal assembly and response to biotic stimuli, and inflammatory immune responses, respectively. At the transcriptome level, T2 vs. T0 showed an upregulation of 337 and downregulation of 61 genes, collectively involved in TNF-α, NFKB, and interleukin signaling. Cell enrichment analysis individuates macrophages, monocytes, and neutrophils as the most significant tissue-cell type in the response.

CONCLUSIONS

In the process of liver graft procurement-transplantation, IRI induces significant epigenetic changes that primarily act on the signaling pathways of inflammatory responses dependent on TNF-α, NFKB, and interleukins. Our DNAm datasets are the early IRI methylome literature and will serve as a launch point for studying the impact of epigenetic modification in IRI.

Transcriptional pathobiology and multi-omics predictors for Parkinson's disease

Early diagnosis and biomarker discovery to bolster the therapeutic pipeline for Parkinson's disease (PD) are urgently needed. In this study, we leverage the large-scale whole-blood total RNA-seq dataset from the Accelerating Medicine Partnership in Parkinson's Disease (AMP PD) program to identify PD-associated RNAs, including both known genes and novel circular RNAs (circRNA) and enhancer RNAs (eRNAs). There were 1,111 significant marker RNAs, including 491 genes, 599 eRNAs, and 21 circRNAs, that were first discovered in the PPMI cohort (FDR < 0.05) and confirmed in the PDBP/BioFIND cohorts (nominal p < 0.05). Functional enrichment analysis showed that the PD-associated genes are involved in neutrophil activation and degranulation, as well as the TNF-alpha signaling pathway. We further compare the PD-associated genes in blood with those in post-mortem brain dopamine neurons in our BRAINcode cohort. 44 genes show significant changes with the same direction in both PD brain neurons and PD blood, including neuroinflammation-associated genes IKBIP, CXCR2, and NFKBIB. Finally, we built a novel multi-omics machine learning model to predict PD diagnosis with high performance (AUC = 0.89), which was superior to previous studies and might aid the decision-making for PD diagnosis in clinical practice. In summary, this study delineates a wide spectrum of the known and novel RNAs linked to PD and are detectable in circulating blood cells in a harmonized, large-scale dataset. It provides a generally useful computational framework for further biomarker development and early disease prediction.

Elucidating immune-related gene transcriptional programs via factorization of large-scale RNA-profiles

Recent developments in immunotherapy, including immune checkpoint blockade (ICB) and adoptive cell therapy (ACT), have encountered challenges such as immune-related adverse events and resistance, especially in solid tumors. To advance the field, a deeper understanding of the molecular mechanisms behind treatment responses and resistance is essential. However, the lack of functionally characterized immune-related gene sets has limited data-driven immunological research. To address this gap, we adopted non-negative matrix factorization on 83 human bulk RNA sequencing (RNA-seq) datasets and constructed 28 immune-specific gene sets. After rigorous immunologist-led manual annotations and orthogonal validations across immunological contexts and functional omics data, we demonstrated that these gene sets can be applied to refine pan-cancer immune subtypes, improve ICB response prediction and functionally annotate spatial transcriptomic data. These functional gene sets, informing diverse immune states, will advance our understanding of immunology and cancer research.

Susceptibility gene identification and risk evaluation model construction by transcriptome-wide association analysis for salt sensitivity of blood pressure

BACKGROUND

Salt sensitivity of blood pressure (SSBP) is an intermediate phenotype of hypertension and is a predictor of long-term cardiovascular events and death. However, the genetic structures of SSBP are uncertain, and it is difficult to precisely diagnose SSBP in population. So, we aimed to identify genes related to susceptibility to the SSBP, construct a risk evaluation model, and explore the potential functions of these genes.

METHODS AND RESULTS

A genome-wide association study of the systemic epidemiology of salt sensitivity (EpiSS) cohort was performed to obtain summary statistics for SSBP. Then, we conducted a transcriptome-wide association study (TWAS) of 12 tissues using FUSION software to predict the genes associated with SSBP and verified the genes with an mRNA microarray. The potential roles of the genes were explored. Risk evaluation models of SSBP were constructed based on the serial P value thresholds of polygenetic risk scores (PRSs), polygenic transcriptome risk scores (PTRSs) and their combinations of the identified genes and genetic variants from the TWAS. The TWAS revealed that 2605 genes were significantly associated with SSBP. Among these genes, 69 were differentially expressed according to the microarray analysis. The functional analysis showed that the genes identified in the TWAS were enriched in metabolic process pathways. The PRSs were correlated with PTRSs in the heart atrial appendage, adrenal gland, EBV-transformed lymphocytes, pituitary, artery coronary, artery tibial and whole blood. Multiple logistic regression models revealed that a PRS of P < 0.05 had the best predictive ability compared with other PRSs and PTRSs. The combinations of PRSs and PTRSs did not significantly increase the prediction accuracy of SSBP in the training and validation datasets.

CONCLUSIONS

Several known and novel susceptibility genes for SSBP were identified via multitissue TWAS analysis. The risk evaluation model constructed with the PRS of susceptibility genes showed better diagnostic performance than the transcript levels, which could be applied to screen for SSBP high-risk individuals.

Human lineage mutations regulate RNA-protein binding of conserved genes NTRK2 and ITPR1 involved in human evolution

Background

The role of human lineage mutations (HLMs) in human evolution through post-transcriptional modification is unclear.

Aims

To investigate the contribution of HLMs to human evolution through post-transcriptional modification.

Methods

We applied a deep learning model Seqweaver to predict how HLMs impact RNA-binding protein affinity.

Results

We found that only 0.27% of HLMs had significant impacts on RNA-binding proteins at the threshold of the top 1% of human common variations. These HLMs enriched in a set of conserved genes highly expressed in adult excitatory neurons and prenatal Purkinje neurons, and were involved in synapse organisation and the GTPase pathway. These genes also carried excess damaging coding mutations that caused neurodevelopmental disorders, ataxia and schizophrenia. Among these genes, NTRK2 and ITPR1 had the most aggregated evidence of functional importance, suggesting their essential roles in cognition and bipedalism.

Conclusions

Our findings suggest that a small subset of human-specific mutations have contributed to human speciation through impacts on post-transcriptional modification of critical brain-related genes.

Elucidating immune-related gene transcriptional programs via factorization of large-scale RNA-profiles

Recent developments in immunotherapy, including immune checkpoint blockade (ICB) and adoptive cell therapy, have encountered challenges such as immune-related adverse events and resistance, especially in solid tumors. To advance the field, a deeper understanding of the molecular mechanisms behind treatment responses and resistance is essential. However, the lack of functionally characterized immune-related gene sets has limited data-driven immunological research. To address this gap, we adopted non-negative matrix factorization on 83 human bulk RNA-seq datasets and constructed 28 immune-specific gene sets. After rigorous immunologist-led manual annotations and orthogonal validations across immunological contexts and functional omics data, we demonstrated that these gene sets can be applied to refine pan-cancer immune subtypes, improve ICB response prediction and functionally annotate spatial transcriptomic data. These functional gene sets, informing diverse immune states, will advance our understanding of immunology and cancer research.

Unraveling transcriptomic signatures and dysregulated pathways in systemic lupus erythematosus across disease states

OBJECTIVES

This study aims to elucidate the transcriptomic signatures and dysregulated pathways in patients with Systemic Lupus Erythematosus (SLE), with a particular focus on those persisting during disease remission.

METHODS

We conducted bulk RNA-sequencing of peripheral blood mononuclear cells (PBMCs) from a well-defined cohort comprising 26 remission patients meeting the Low Lupus Disease Activity State (LLDAS) criteria, 76 patients experiencing disease flares, and 15 healthy controls. To elucidate immune signature changes associated with varying disease states, we performed extensive analyses, including the identification of differentially expressed genes and pathways, as well as the construction of protein-protein interaction networks.

RESULTS

Several transcriptomic features recovered during remission compared to the active disease state, including down-regulation of plasma and cell cycle signatures, as well as up-regulation of lymphocytes. However, specific innate immune response signatures, such as the interferon (IFN) signature, and gene modules involved in chromatin structure modification, persisted across different disease states. Drug repurposing analysis revealed certain drug classes that can target these persistent signatures, potentially preventing disease relapse.

CONCLUSION

Our comprehensive transcriptomic study revealed gene expression signatures for SLE in both active and remission states. The discovery of gene expression modules persisting in the remission stage may shed light on the underlying mechanisms of vulnerability to relapse in these patients, providing valuable insights for their treatment.

Unveiling Gene Interactions in Alzheimer's Disease by Integrating Genetic and Epigenetic Data with a Network-Based Approach

Alzheimer's Disease (AD) is a complex disease and the leading cause of dementia in older people. We aimed to uncover aspects of AD's pathogenesis that may contribute to drug repurposing efforts by integrating DNA methylation and genetic data. Implementing the network-based tool, a dense module search of genome-wide association studies (dmGWAS), we integrated a large-scale GWAS dataset with DNA methylation data to identify gene network modules associated with AD. Our analysis yielded 286 significant gene network modules. Notably, the foremost module included the BIN1 gene, showing the largest GWAS signal, and the GNAS gene, the most significantly hypermethylated. We conducted Web-based Cell-type-Specific Enrichment Analysis (WebCSEA) on genes within the top 10% of dmGWAS modules, highlighting monocyte as the most significant cell type (p < 5 × 10-12). Functional enrichment analysis revealed Gene Ontology Biological Process terms relevant to AD pathology (adjusted p < 0.05). Additionally, drug target enrichment identified five FDA-approved targets (p-value = 0.03) for further research. In summary, dmGWAS integration of genetic and epigenetic signals unveiled new gene interactions related to AD, offering promising avenues for future studies.

Drug Repositioning of Inflammatory Bowel Disease Based on Co-Target Gene Expression Signature of Glucocorticoid Receptor and TET2

The glucocorticoid receptor (GR) and ten-eleven translocation 2 (TET2), respectively, play a crucial role in regulating immunity and inflammation, and GR interacts with TET2. However, their synergetic roles in inflammatory bowel disease (IBD), including ulcerative colitis (UC) and Crohn's disease (CD), remain unclear. This study aimed to investigate the co-target gene signatures of GR and TET2 in IBD and provide potential therapeutic interventions for IBD. By integrating public data, we identified 179 GR- and TET2-targeted differentially expressed genes (DEGs) in CD and 401 in UC. These genes were found to be closely associated with immunometabolism, inflammatory responses, and cell stress pathways. In vitro inflammatory cellular models were constructed using LPS-treated HT29 and HCT116 cells, respectively. Drug repositioning based on the co-target gene signatures of GR and TET2 derived from transcriptomic data of UC, CD, and the in vitro model was performed using the Connectivity Map (CMap). BMS-536924 emerged as a top therapeutic candidate, and its validation experiment within the in vitro inflammatory model confirmed its efficacy in mitigating the LPS-induced inflammatory response. This study sheds light on the pathogenesis of IBD from a new perspective and may accelerate the development of novel therapeutic agents for inflammatory diseases including IBD.

COUP-TFII regulates early bipotential gonad signaling and commitment to ovarian progenitors

BACKGROUND

The absence of expression of the Y-chromosome linked testis-determining gene SRY in early supporting gonadal cells (ESGC) leads bipotential gonads into ovarian development. However, genetic variants in NR2F2, encoding three isoforms of the transcription factor COUP-TFII, represent a novel cause of SRY-negative 46,XX testicular/ovotesticular differences of sex development (T/OT-DSD). Thus, we hypothesized that COUP-TFII is part of the ovarian developmental network. COUP-TFII is known to be expressed in interstitial/mesenchymal cells giving rise to steroidogenic cells in fetal gonads, however its expression and function in ESGCs have yet to be explored.

RESULTS

By differentiating induced pluripotent stem cells into bipotential gonad-like cells in vitro and by analyzing single cell RNA-sequencing datasets of human fetal gonads, we identified that NR2F2 expression is highly upregulated during bipotential gonad development along with markers of bipotential state. NR2F2 expression was detected in early cell populations that precede the steroidogenic cell emergence and that retain a multipotent state in the undifferentiated gonad. The ESGCs differentiating into fetal Sertoli cells lost NR2F2 expression, whereas pre-granulosa cells remained NR2F2-positive. When examining the NR2F2 transcript variants individually, we demonstrated that the canonical isoform A, disrupted by frameshift variants previously reported in 46,XX T/OT-DSD patients, is nearly 1000-fold more highly expressed than other isoforms in bipotential gonad-like cells. To investigate the genetic network under COUP-TFII regulation in human gonadal cell context, we generated a NR2F2 knockout (KO) in the human granulosa-like cell line COV434 and studied NR2F2-KO COV434 cell transcriptome. NR2F2 ablation downregulated markers of ESGC and pre-granulosa cells. NR2F2-KO COV434 cells lost the enrichment for female-supporting gonadal progenitor and acquired gene signatures more similar to gonadal interstitial cells.

CONCLUSIONS

Our findings suggest that COUP-TFII has a role in maintaining a multipotent state necessary for commitment to the ovarian development. We propose that COUP-TFII regulates cell fate during gonad development and impairment of its function may disrupt the transcriptional plasticity of ESGCs. During early gonad development, disruption of ESGC plasticity may drive them into commitment to the testicular pathway, as observed in 46,XX OT-DSD patients with NR2F2 haploinsufficiency.

Monocyte-derived transcriptomes explain the ineffectiveness of abatacept in rheumatoid arthritis

BACKGROUND

The biological mechanisms underlying the differential response to abatacept in patients with rheumatoid arthritis (RA) are unknown. Here, we aimed to identify cellular, transcriptomic, and proteomic features that predict resistance to abatacept in patients with RA.

METHODS

Blood samples were collected from 22 RA patients treated with abatacept at baseline and after 3 months of treatment. Response to treatment was defined by the European League Against Rheumatism (EULAR) response criteria at 3 months, and seven patients were classified as responders and the others as non-responders. We quantified gene expression levels by RNA sequencing, 67 plasma protein levels, and the expression of surface molecules (CD3, 19, and 56) by flow cytometry. In addition, three gene expression data sets, comprising a total of 27 responders and 50 non-responders, were used to replicate the results.

RESULTS

Among the clinical characteristics, the number of monocytes was significantly higher in the non-responders before treatment. Cell type enrichment analysis showed that differentially expressed genes (DEGs) between responders and non-responders were enriched in monocytes. Gene set enrichment analysis, together with single-cell analysis and deconvolution analysis, identified that Toll-like receptor 5 (TLR5) and interleukin-17 receptor A (IL17RA) pathway in monocytes was upregulated in non-responders. Hepatocyte growth factor (HGF) correlated with this signature showed higher concentrations in non-responders before treatment. The DEGs in the replication set were also enriched for the genes expressed in monocytes, not for the TLR5 and IL17RA pathway but for the oxidative phosphorylation (OXPHOS) pathway.

CONCLUSIONS

Monocyte-derived transcriptomic features before treatment underlie the differences in abatacept efficacy in patients with RA. The pathway activated in monocytes was the TLR5 and IL17RA-HGF signature in the current study, while it was the OXPHOS pathway in the replication set. Elevated levels of HGF before treatment may serve as a potential biomarker for predicting poor responses to abatacept. These findings provide insights into the biological mechanisms of abatacept resistance, contributing valuable evidence for stratifying patients with RA.

Disentangling Accelerated Cognitive Decline from the Normal Aging Process and Unraveling Its Genetic Components: A Neuroimaging-Based Deep Learning Approach

Background

The progressive cognitive decline, an integral component of Alzheimer's disease (AD), unfolds in tandem with the natural aging process. Neuroimaging features have demonstrated the capacity to distinguish cognitive decline changes stemming from typical brain aging and AD between different chronological points.

Objective

To disentangle the normal aging effect from the AD-related accelerated cognitive decline and unravel its genetic components using a neuroimaging-based deep learning approach.

Methods

We developed a deep-learning framework based on a dual-loss Siamese ResNet network to extract fine-grained information from the longitudinal structural magnetic resonance imaging (MRI) data from the Alzheimer's Disease Neuroimaging Initiative (ADNI) study. We then conducted genome-wide association studies (GWAS) and post-GWAS analyses to reveal the genetic basis of AD-related accelerated cognitive decline.

Results

We used our model to process data from 1,313 individuals, training it on 414 cognitively normal people and predicting cognitive assessment for all participants. In our analysis of accelerated cognitive decline GWAS, we identified two genome-wide significant loci: APOE locus (chromosome 19 p13.32) and rs144614292 (chromosome 11 p15.1). Variant rs144614292 (G > T) has not been reported in previous AD GWA studies. It is within the intronic region of NELL1, which is expressed in neurons and plays a role in controlling cell growth and differentiation. The cell-type-specific enrichment analysis and functional enrichment of GWAS signals highlighted the microglia and immune-response pathways.

Conclusions

Our deep learning model effectively extracted relevant neuroimaging features and predicted individual cognitive decline. We reported a novel variant (rs144614292) within the NELL1 gene.

Transcriptomic Landscape of Circulating Extracellular Vesicles in Heart Transplant Ischemia-Reperfusion

Ischemia-reperfusion injury (IRI) is an inevitable event during heart transplantation, which is known to exacerbate damage to the allograft. However, the precise mechanisms underlying IRI remain incompletely understood. Here, we profiled the whole transcriptome of plasma extracellular vesicles (EVs) by RNA sequencing from 41 heart transplant recipients immediately before and at 12 h after transplant reperfusion. We found that the expression of 1317 protein-coding genes in plasma EVs was changed at 12 h after reperfusion. Upregulated genes of plasma EVs were related to metabolism and immune activation, while downregulated genes were related to cell survival and extracellular matrix organization. In addition, we performed correlation analyses between EV transcriptome and intensity of graft IRI (i.e., cardiomyocyte injury), as well as EV transcriptome and primary graft dysfunction, as well as any biopsy-proven acute rejection after heart transplantation. We ultimately revealed that at 12 h after reperfusion, 4 plasma EV genes (ITPKA, DDIT4L, CD19, and CYP4A11) correlated with both cardiomyocyte injury and primary graft dysfunction, suggesting that EVs are sensitive indicators of reperfusion injury reflecting lipid metabolism-induced stress and imbalance in calcium homeostasis. In conclusion, we show that profiling plasma EV gene expression may enlighten the mechanisms of heart transplant IRI.

Investigation into Molecular Brain Aging in Senescence-Accelerated Mouse (SAM) Model Employing Whole Transcriptomic Analysis in Search of Potential Molecular Targets for Therapeutic Interventions

With the progression of an aging society, cognitive aging has emerged as a pressing concern necessitating attention. The senescence-accelerated mouse-prone 8 (SAMP8) model has proven instrumental in investigating the early stages of cognitive aging. Through an extensive examination of molecular changes in the brain cortex, utilizing integrated whole-genome transcriptomics, our principal aim was to uncover potential molecular targets with therapeutic applications and relevance to drug screening. Our investigation encompassed four distinct conditions, comparing the same strain at different time points (1 year vs. 16 weeks) and the same time point across different strains (SAMP8 vs. SAMR1), namely: physiological aging, accelerated aging, early events in accelerated aging, and late events in accelerated aging. Focusing on key functional alterations associated with aging in the brain, including neurogenesis, synapse dynamics, neurometabolism, and neuroinflammation, we identified candidate genes linked to these processes. Furthermore, employing protein-protein interaction (PPI) analysis, we identified pivotal hub genes involved in interactions within these functional domains. Additionally, gene-set perturbation analysis allowed us to uncover potential upstream genes or transcription factors that exhibited activation or inhibition across the four conditions. In summary, our comprehensive analysis of the SAMP8 mouse brain through whole-genome transcriptomics not only deepens our understanding of age-related changes but also lays the groundwork for a predictive model to facilitate drug screening for cognitive aging.

Disentangling accelerated cognitive decline from the normal aging process and unraveling its genetic components: A neuroimaging-based deep learning approach

Background

The progressive cognitive decline that is an integral component of AD unfolds in tandem with the natural aging process. Neuroimaging features have demonstrated the capacity to distinguish cognitive decline changes stemming from typical brain aging and Alzheimer's disease between different chronological points.

Methods

We developed a deep-learning framework based on dual-loss Siamese ResNet network to extract fine-grained information from the longitudinal structural magnetic resonance imaging (MRI) data from the Alzheimer's Disease Neuroimaging Initiative (ADNI) study. We then conducted genome-wide association studies (GWAS) and post-GWAS analyses to reveal the genetic basis of AD-related accelerated cognitive decline.

Results

We used our model to process data from 1,313 individuals, training it on 414 cognitively normal people and predicting cognitive assessment for all participants. In our analysis of accelerated cognitive decline GWAS, we identified two genome-wide significant loci: APOE locus (chromosome 19 p13.32) and rs144614292 (chromosome 11 p15.1). Variant rs144614292 (G>T) has not been reported in previous AD GWA studies. It is within the intronic region of NELL1, which is expressed in neuron and plays a role in controlling cell growth and differentiation. In addition, MUC7 and PROL1/OPRPNon chromosome 4 were significant at the gene level. The cell-type-specific enrichment analysis and functional enrichment of GWAS signals highlighted the microglia and immune-response pathways. Furthermore, we found that the cognitive decline slope GWAS was positively correlated with previous AD GWAS.

Conclusion

Our deep learning model was demonstrated effective on extracting relevant neuroimaging features and predicting individual cognitive decline. We reported a novel variant (rs144614292) within the NELL1 gene. Our approach has the potential to disentangle accelerated cognitive decline from the normal aging process and to determine its related genetic factors, leveraging opportunities for early intervention.

Longitudinal APOE4- and amyloid-dependent changes in the blood transcriptome in cognitively intact older adults

BACKGROUND

Gene expression is dysregulated in Alzheimer's disease (AD) patients, both in peripheral blood and post mortem brain. We investigated peripheral whole-blood gene (co)expression to determine molecular changes prior to symptom onset.

METHODS

RNA was extracted and sequenced for 65 cognitively healthy F-PACK participants (65 (56-80) years, 34 APOE4 non-carriers, 31 APOE4 carriers), at baseline and follow-up (interval: 5.0 (3.4-8.6) years). Participants received amyloid PET at both time points and amyloid rate of change derived. Accumulators were defined with rate of change ≥ 2.19 Centiloids. We performed differential gene expression and weighted gene co-expression network analysis to identify differentially expressed genes and networks of co-expressed genes, respectively, with respect to traits of interest (APOE4 status, amyloid accumulation (binary/continuous)), and amyloid positivity status, followed by Gene Ontology annotation.

RESULTS

There were 166 significant differentially expressed genes at follow-up compared to baseline in APOE4 carriers only, whereas 12 significant differentially expressed genes were found only in APOE4 non-carriers, over time. Among the significant genes in APOE4 carriers, several had strong evidence for a pathogenic role in AD based on direct association scores generated from the DISQOVER platform: NGRN, IGF2, GMPR, CLDN5, SMIM24. Top enrichment terms showed upregulated mitochondrial and metabolic pathways, and an exacerbated upregulation of ribosomal pathways in APOE4 carriers compared to non-carriers. Similarly, there were 33 unique significant differentially expressed genes at follow-up compared to baseline in individuals classified as amyloid negative at baseline and positive at follow-up or amyloid positive at both time points and 32 unique significant differentially expressed genes over time in individuals amyloid negative at both time points. Among the significant genes in the first group, the top five with the highest direct association scores were as follows: RPL17-C18orf32, HSP90AA1, MBP, SIRPB1, and GRINA. Top enrichment terms included upregulated metabolism and focal adhesion pathways. Baseline and follow-up gene co-expression networks were separately built. Seventeen baseline co-expression modules were derived, with one significantly negatively associated with amyloid accumulator status (r² =  - 0.25, p = 0.046). This was enriched for proteasomal protein catabolic process and myeloid cell development. Thirty-two follow-up modules were derived, with two significantly associated with APOE4 status: one downregulated (r² =  - 0.27, p = 0.035) and one upregulated (r² = 0.26, p = 0.039) module. Top enrichment processes for the downregulated module included proteasomal protein catabolic process and myeloid cell homeostasis. Top enrichment processes for the upregulated module included cytoplasmic translation and rRNA processing.

CONCLUSIONS

We show that there are longitudinal gene expression changes that implicate a disrupted immune system, protein removal, and metabolism in cognitively intact individuals who carry APOE4 or who accumulate in cortical amyloid. This provides insight into the pathophysiology of AD, whilst providing novel targets for drug and therapeutic development.

Gene-based association study reveals a distinct female genetic signal in primary hypertension

Hypertension is a polygenic disease that affects over 1.2 billion adults aged 30-79 worldwide. It is a major risk factor for renal, cerebrovascular, and cardiovascular diseases. The heritability of hypertension is estimated to be high; nevertheless, our understanding of its underlying mechanisms remains scarce and incomplete. This study covered the entries from European ancestry from the UK-Biobank (UKB), with 74,090 cases diagnosed with essential (primary) hypertension and 200,734 controls. We compared the findings from large-scale genome-wide association studies (GWAS) to the gene-based method of proteome-wide association studies (PWAS). We focused on 70 statistically significant associated genes, most of which failed to reach significance in variant-based GWAS. A total of 30% of the PWAS-associated genes were validated against independent cohorts, including the Finnish Biobank. Furthermore, gene-based analyses that were performed on both sexes revealed sex-dependent genetics with a stronger genetic component associated with females. Analysis of systolic and diastolic blood pressure measurements confirms a strong genetic effect associated with females. We demonstrated that gene-based approaches provide insight into the underlying biology of hypertension. Specifically, the expression profiles of the identified genes exposed the enrichment of endothelial cells from multiple organs. Furthermore, females' top-ranked significant genes are involved in cellular immunity. We conclude that studying hypertension and blood pressure via gene-based association methods improves interpretability and exposes sex-dependent genetic effects, which enhances clinical utility.

De novo mutations disturb early brain development more frequently than common variants in schizophrenia

Investigating functional, temporal, and cell-type expression features of mutations is important for understanding a complex disease. Here, we collected and analyzed common variants and de novo mutations (DNMs) in schizophrenia (SCZ). We collected 2,636 missense and loss-of-function (LoF) DNMs in 2,263 genes across 3,477 SCZ patients (SCZ-DNMs). We curated three gene lists: (a) SCZ-neuroGenes (159 genes), which are intolerant to LoF and missense DNMs and are neurologically important, (b) SCZ-moduleGenes (52 genes), which were derived from network analyses of SCZ-DNMs, and (c) SCZ-commonGenes (120 genes) from a recent GWAS as reference. To compare temporal gene expression, we used the BrainSpan dataset. We defined a fetal effect score (FES) to quantify the involvement of each gene in prenatal brain development. We further employed the specificity indexes (SIs) to evaluate cell-type expression specificity from single-cell expression data in cerebral cortices of humans and mice. Compared with SCZ-commonGenes, SCZ-neuroGenes and SCZ-moduleGenes were highly expressed in the prenatal stage, had higher FESs, and had higher SIs in fetal replicating cells and undifferentiated cell types. Our results suggested that gene expression patterns in specific cell types in early fetal stages might have impacts on the risk of SCZ during adulthood.

Genome-wide analysis of genetic pleiotropy and causal genes across three age-related ocular disorders

Age-related macular degeneration (AMD), cataract, and glaucoma are leading causes of blindness worldwide. Previous genome-wide association studies (GWASs) have revealed a variety of susceptible loci associated with age-related ocular disorders, yet the genetic pleiotropy and causal genes across these diseases remain poorly understood. By leveraging large-scale genetic and observational data from ocular disease GWASs and UK Biobank (UKBB), we found significant pairwise genetic correlations and consistent epidemiological associations among these ocular disorders. Cross-disease meta-analysis uncovered seven pleiotropic loci, three of which were replicated in an additional cohort. Integration of variants in pleiotropic loci and multiple single-cell omics data identified that Müller cells and astrocytes were likely trait-related cell types underlying ocular comorbidity. In addition, we comprehensively integrated eye-specific gene expression quantitative loci (eQTLs), epigenomic profiling, and 3D genome data to prioritize causal pleiotropic genes. We found that pleiotropic genes were essential in nerve development and eye pigmentation, and targetable by aflibercept and pilocarpine for the treatment of AMD and glaucoma. These findings will not only facilitate the mechanistic research of ocular comorbidities but also benefit the therapeutic optimization of age-related ocular diseases.

A gene regulatory network approach harmonizes genetic and epigenetic signals and reveals repurposable drug candidates for multiple sclerosis

Multiple sclerosis (MS) is a complex dysimmune disorder of the central nervous system. Genome-wide association studies (GWAS) have identified 233 genetic variations associated with MS at the genome-wide significant level. Epigenetic studies have pinpointed differentially methylated CpG sites in MS patients. However, the interplay between genetic risk factors and epigenetic regulation remains elusive. Here, we employed a network model to integrate GWAS summary statistics of 14 802 MS cases and 26 703 controls with DNA methylation profiles from 140 MS cases and 139 controls and the human interactome. We identified differentially methylated genes by aggregating additive effects of differentially methylated CpG sites within promoter regions. We reconstructed a gene regulatory network (GRN) using literature-curated transcription factor knowledge. Colocalization of the MS GWAS and methylation quantitative trait loci (mQTL) was performed to assess the GRN. The resultant MS-associated GRN highlighted several single nucleotide polymorphisms with GWAS-mQTL colocalization: rs6032663, rs6065926 and rs2024568 of CD40 locus, rs9913597 of STAT3 locus, and rs887864 and rs741175 of CIITA locus. Moreover, synergistic mQTL and expression QTL signals were identified in CD40, suggesting gene expression alteration was likely induced by epigenetic changes. Web-based Cell-type Specific Enrichment Analysis of Genes (WebCSEA) indicated that the GRN was enriched in T follicular helper cells (P-value = 0.0016). Drug target enrichment analysis of annotations from the Therapeutic Target Database revealed the GRN was also enriched with drug target genes (P-value = 3.89 × 10-4), revealing repurposable candidates for MS treatment. These candidates included vorinostat (HDAC1 inhibitor) and sivelestat (ELANE inhibitor), which warrant further investigation.

New Insight Into Mechanisms of Hepatic Encephalopathy: An Integrative Analysis Approach to Identify Molecular Markers and Therapeutic Targets

Hepatic encephalopathy (HE) is a set of complex neurological complications that arise from advanced liver disease. The precise molecular and cellular mechanism of HE is not fully understood. Differentially expressed genes (DEGs) from microarray technologies are powerful approaches to obtain new insight into the pathophysiology of HE. We analyzed microarray data sets of cirrhotic patients with HE from Gene Expression Omnibus to identify DEGs in postmortem cerebral tissues. Consequently, we uploaded significant DEGs into the STRING to specify protein-protein interactions. Cytoscape was used to reconstruct the genetic network and identify hub genes. Target genes were uploaded to different databases to perform comprehensive enrichment analysis and repurpose new therapeutic options for HE. A total of 457 DEGs were identified in 2 data sets totally from 12 cirrhotic patients with HE compared with 12 healthy subjects. We found that 274 genes were upregulated and 183 genes were downregulated. Network analyses on significant DEGs indicated 12 hub genes associated with HE. Enrichment analysis identified fatty acid beta-oxidation, cerebral organic acidurias, and regulation of actin cytoskeleton as main involved pathways associated with upregulated genes; serotonin receptor 2 and ELK-SRF/GATA4 signaling, GPCRs, class A rhodopsin-like, and p38 MAPK signaling pathway were related to downregulated genes. Finally, we predicted 39 probable effective drugs/agents for HE. This study not only confirms main important involved mechanisms of HE but also reveals some yet unknown activated molecular and cellular pathways in human HE. In addition, new targets were identified that could be of value in the future study of HE.

Transcriptional and cell type profiles of cortical brain regions showing ultradian cortisol rhythm dependent responses to emotional face stimulation

The characteristic endogenous circadian rhythm of plasma glucocorticoid concentrations is made up from an underlying ultradian pulsatile secretory pattern. Recent evidence has indicated that this ultradian cortisol pulsatility is crucial for normal emotional response in man. In this study, we investigate the anatomical transcriptional and cell type signature of brain regions sensitive to a loss of ultradian rhythmicity in the context of emotional processing. We combine human cell type and transcriptomic atlas data of high spatial resolution with functional magnetic resonance imaging (fMRI) data. We show that the loss of cortisol ultradian rhythm alters emotional processing response in cortical brain areas that are characterized by transcriptional and cellular profiles of GABAergic function. We find that two previously identified key components of rapid non-genomic GC signaling - the ANXA1 gene and retrograde endocannabinoid signaling - show most significant differential expression (q = 3.99e-10) and enrichment (fold enrichment = 5.56, q = 9.09e-4). Our results further indicate that specific cell types, including a specific NPY-expressing GABAergic neuronal cell type, and specific G protein signaling cascades underly the cerebral effects of a loss of ultradian cortisol rhythm. Our results provide a biological mechanistic underpinning of our fMRI findings, indicating specific cell types and cascades as a target for manipulation in future experimental studies.

Drug-Target Network Study Reveals the Core Target-Protein Interactions of Various COVID-19 Treatments

The coronavirus disease 2019 (COVID-19) pandemic has caused a dramatic loss of human life and devastated the worldwide economy. Numerous efforts have been made to mitigate COVID-19 symptoms and reduce the death rate. We conducted literature mining of more than 250 thousand published works and curated the 174 most widely used COVID-19 medications. Overlaid with the human protein-protein interaction (PPI) network, we used Steiner tree analysis to extract a core subnetwork that grew from the pharmacological targets of ten credible drugs ascertained by the CTD database. The resultant core subnetwork consisted of 34 interconnected genes, which were associated with 36 drugs. Immune cell membrane receptors, the downstream cellular signaling cascade, and severe COVID-19 symptom risk were significantly enriched for the core subnetwork genes. The lung mast cell was most enriched for the target genes among 1355 human tissue-cell types. Human bronchoalveolar lavage fluid COVID-19 single-cell RNA-Seq data highlighted the fact that T cells and macrophages have the most overlapping genes from the core subnetwork. Overall, we constructed an actionable human target-protein module that mainly involved anti-inflammatory/antiviral entry functions and highly overlapped with COVID-19-severity-related genes. Our findings could serve as a knowledge base for guiding drug discovery or drug repurposing to confront the fast-evolving SARS-CoV-2 virus and other severe infectious diseases.