Identifying causal genes for stroke via integrating the proteome and transcriptome from brain and blood

Differential Expression of mRNA in Peripheral Blood Mononuclear Cells May Predict Postoperative Atrial Fibrillation in Coronary Artery Bypass Surgery Patients

Postoperative Atrial Fibrillation (POAF) frequently follows Coronary Artery Bypass Grafting (CABG) surgery. This prospective study investigates genes linked to POAF in CABG patients, aiming to create a predictive model. Employing differential gene and methylation analyses, the study identified four genes (WARS2, CKAP2, CHI3L1, HSD17B6) associated with POAF. Preoperative plasma samples and clinical data were collected from 139 CABG patients, categorized into POAF (+) (43) and POAF (-) (96). Real-time quantitative PCR assessed gene expression, and a predictive model using the LASSO method demonstrated robust performance, with AUC values of 0.8895 in the training set and 0.7840 in the test set. This pioneering study integrates genomics and clinical data, suggesting WARS2, CKAP2, and CHI3L1 as potential indicators for POAF prediction.

Identifying novel potential drug targets for endometriosis via plasma proteome screening

Background

Endometriosis (EM) is a chronic painful condition that predominantly affects women of reproductive age. Currently, surgery or medication can only provide limited symptom relief. This study used a comprehensive genetic analytical approach to explore potential drug targets for EM in the plasma proteome.

Methods

In this study, 2,923 plasma proteins were selected as exposure and EM as outcome for two-sample Mendelian randomization (MR) analyses. The plasma proteomic data were derived from the UK Biobank Pharmaceutical Proteomics Project (UKB-PPP), while the EM dataset from the FinnGen consortium R10 release data. Several sensitivity analyses were performed, including summary-data-based MR (SMR) analyses, heterogeneity in dependent instruments (HEIDI) test, reverse MR analyses, steiger detection test, and bayesian co-localization analyses. Furthermore, proteome-wide association study (PWAS) and single-cell transcriptomic analyses were also conducted to validate the findings.

Results

Six significant (p < 3.06 × 10-5) plasma protein-EM pairs were identified by MR analyses. These included EPHB4 (OR = 1.40, 95% CI: 1.20 - 1.63), FSHB (OR = 3.91, 95% CI: 3.13 - 4.87), RSPO3 (OR = 1.60, 95% CI: 1.38 - 1.86), SEZ6L2 (OR = 1.44, 95% CI: 1.23 - 1.68) and WASHC3 (OR = 2.00, 95% CI: 1.54 - 2.59) were identified as risk factors, whereas KDR (OR = 0.80, 95% CI: 0.75 - 0.90) was found to be a protective factor. All six plasma proteins passed the SMR test (P < 8.33 × 10-3), but only four plasma proteins passed the HEIDI heterogeneity test (PHEIDI > 0.05), namely FSHB, RSPO3, SEZ6L2 and EPHB4. These four proteins showed strong evidence of co-localization (PPH4 > 0.7). In particular, RSPO3 and EPHB4 were replicated in the validated PWAS. Single-cell analyses revealed high expression of SEZ6L2 and EPHB4 in stromal and epithelial cells within EM lesions, while RSPO3 exhibited elevated expression in stromal cells and fibroblasts.

Conclusion

Our study identified FSHB, RSPO3, SEZ6L2, and EPHB4 as potential drug targets for EM and highlighted the critical role of stromal and epithelial cells in disease development. These findings provide new insights into the diagnosis and treatment of EM.

Non-O blood types are associated with a greater risk of large artery atherosclerosis stroke and dysregulation of cholesterol metabolism: an observational study

BACKGROUND

Previous research on ABO blood types and stroke has been controversial, predominantly suggesting heightened risk of stroke in non-O blood types. Nonetheless, investigations into the correlation and underlying mechanisms between ABO blood groups and stroke subtypes, especially within Chinese cohorts, remain limited.

METHODS

The ABO blood types of 9,542 ischaemic stroke (IS) patients were inferred using two ABO gene loci (c.261G > del; c.802G > A). The healthy population was derived from the 1000 Genomes Project. Patients were classified by the causative classification system (CCS). Volcano plot and gene ontology (GO) analysis were employed to explore protein differential expression among blood types. Additionally, HT29 and SW480 cell lines with downregulated ABO expression were generated to evaluate its impact on cholesterol uptake and efflux.

RESULTS

A greater proportion of stroke patients had non-O blood types (70.46%) than did healthy individuals (61.54%). Notable differences in blood type distributions were observed among stroke subtypes, with non-O blood type patients mainly classified as having large artery atherosclerosis (LAA). Clinical baseline characteristics, such as the low-density lipoprotein cholesterol level, activated partial thromboplastin time and thrombin time, varied significantly among blood types. A volcano plot revealed 17 upregulated and 42 downregulated proteins in the O blood type. GO term analysis indicated that downregulated proteins were primarily associated with lipid metabolism pathways. In vitro experiments revealed that reducing ABO gene expression decreased cholesterol uptake and increased cholesterol efflux.

CONCLUSIONS

This study revealed that the non-O blood type increased the risk of LAA stroke through cholesterol metabolism.

Identifying novel proteins for suicide attempt by integrating proteomes from brain and blood with genome-wide association data

Genome-wide association studies (GWASs) have identified risk loci for suicide attempt (SA), but deciphering how they confer risk for SA remains largely unknown. This study aims to identify the key proteins and gain insights into SA pathogenesis. We integrated data from the brain proteome (N = 376) and blood proteome (N = 35,559) and combined it with the largest SA GWAS summary statistics to date (N = 518,612). A comprehensive set of methods was employed, including Mendelian randomization (MR), Steiger filtering, Bayesian colocalization, proteome‑wide association studies (PWAS), transcript-levels, cell-type specificity, correlation, and protein-protein interaction (PPI) network analysis. Validation was performed using other protein datasets and the SA dataset from FinnGen study. We identified ten proteins (GLRX5, GMPPB, B3GALTL, FUCA2, TTLL12, ADCK1, MMAA, HIBADH, ACP1, DOC2A) associated with SA in brain proteomics. GLRX5, GMPPB, and FUCA2 showed strong colocalization evidence and were supported by PWAS and transcript-level analysis, and were predominantly expressed in glutamatergic neuronal cells. In blood proteomics, one significant protein (PEAR1) and three near-significant proteins (NDE1, EVA1C, B4GALT2) were identified, but lacked colocalization evidence. Moreover, despite the limited correlation between the same protein in brain and blood, the PPI network analysis provided new insights into the interaction between brain and blood in SA. Furthermore, GLRX5 was associated with the GSTP1, the target of Clozapine. The comprehensive analysis provides strong evidence supporting a causal association between three genetically determined brain proteins (GLRX5, GMPPB, and FUCA2) with SA. These findings offer valuable insights into SA's underlying mechanisms and potential therapeutic approaches.

Bayesian colocalization of GWAS and eQTL signals reveals cell type-specific genes and regulatory variants for susceptibility to subtypes of ischemic stroke

A colocalization analysis of genome-wide association study (GWAS) signals and expression quantitative trait loci (eQTL) was conducted to pinpoint target genes and their regulatory nucleotide variants for subtypes of ischemic stroke. We utilized GWAS data from prominent meta-analysis consortia (MEGASTROKE and GIGASTROKE) and single-cell eQTL data in brain and blood tissues to enhance accuracy and minimize noise inherent in bulk RNA-seq. Employing Bayesian colocalization methods, we identified ten shared loci between GWAS and eQTL signals, targeting five eGenes. Specifically, RAPH1 and ICA1L were discovered for small vessel stroke (SVS), whereas SCYL3, CAV1, and CAV2 were for cardioembolic stroke (CS). However, no findings have been made for large artery stroke. The exploration and subsequent functional analysis of causal variants within the colocalized regions revealed their regulatory roles, particularly as enhancer variants (e.g., rs144505847 and rs72932755 targeting ICA1L; rs629234 targeting SCYL3; rs3807989 targeting CAV1 and CAV2). Notably, our study unveiled that all eQTL for CS were identified in oligodendrocytes, while those for SVS were across excitatory neurons, astrocytes, and oligodendrocyte precursor cells. This underscores the heterogeneous tissue-specific genetic factors by subtypes of ischemic stroke. The study emphasizes the need for intensive research efforts to discover causative genes and variants, unravelling the cell type-specific genetic architecture of ischemic stroke subtypes. This knowledge is crucial for advancing our understanding of the underlying pathophysiology and paving the way for precision neurology applications.

Traditional Chinese Medicine formula, Sanwujiao granule, attenuates ischemic stroke by promoting angiogenesis through early administration

ETHNOPHARMACOLOGICAL RELEVANCE

Ischemic stroke (IS) is one of the most lethal diseases with the insufficient pharmacology therapeutic approach. Sanwujiao granule (SW) is widely used for IS in China with little known about its underlying mechanism.

AIM OF THE STUDY

To investigate the characteristics of therapeutic effects and potential mechanisms of SW against IS.

MATERIALS AND METHODS

The fingerprint of SW was applied by high-performance liquid chromatography-mass spectrometry (HPLC-MS). Three different drug treatment strategies, including prophylactic administration, early administration and delayed administration, were applied in rats' permanent middle cerebral occlusion (pMCAO) model. The Garcia neurological deficit test, adhesive removal test, rotarod test, TTC and TUNEL staining were performed to evaluate the pathological changes. The transcriptomic analysis was used to predict the potential mechanism of SW. The vascular deficiency model of Tg(kdrl:eGFP) zebrafish larvae and oxygen-glucose deprivation model on bEnd.3 cells were used to verify SW's pharmacological effect. qRT-PCR, immunofluorescent staining and Western Blot were applied to detect the expression of genes and proteins. The network pharmacology approach was applied to discover the potential bioactive compounds in SW that contribute to its pharmacological effect.

RESULTS

SW early and delayed administration attenuated cerebral infarction, neurological deficit and cell apoptosis. The transcriptomic analysis revealed that SW activated angiogenesis-associated biological processes specifically by early administration. CD31 immunofluorescent staining further confirmed the microvessel intensity in peri-infarct regions was significantly elevated after SW early treatment. Additionally, on the vascular deficiency model of zebrafish larvae, SW showed the angiogenesis effect. Next, the cell migration and tube formation were also observed in the bEnd.3 cells with the oxygen-glucose deprivation induced cell injury. It's worth noting that both mRNA and protein levels of angiogenesis factor, insulin-like growth factor 1, were significantly elevated in the pMCAO rats' brains treated with SW. The network pharmacology approach was applied and chasmanine, karacoline, talatisamine, etc. were probably the main active compounds of SW in IS treatment as they affected the angiogenesis-associated targets.

CONCLUSIONS

These results demonstrate that SW plays a critical role in anti-IS via promoting angiogenesis through early administration, indicating that SW is a candidate herbal complex for further investigation in treating IS in the clinical.

The post-stroke young adult brain has limited capacity to re-express the gene expression patterns seen during early postnatal brain development

The developmental origins of the brain's response to injury can play an important role in recovery after a brain lesion. In this study, we investigated whether the ischemic young adult brain can re-express brain plasticity genes that were active during early postnatal development. Differentially expressed genes in the cortex of juvenile post-natal day 3 and the peri-infarcted cortical areas of young, 3-month-old post-stroke rats were identified using fixed-effects modeling within an empirical Bayes framework through condition-specific comparison. To further analyze potential biological processes, upregulated and downregulated genes were assessed for enrichment using GSEA software. The genes showing the highest expression changes were subsequently verified through RT-PCR. Our findings indicate that the adult brain partially recapitulates the gene expression profile observed in the juvenile brain but fails to upregulate many genes and pathways necessary for brain plasticity. Of the upregulated genes in post-stroke brains, specific roles have not been assigned to Apobec1, Cenpf, Ect2, Folr2, Glipr1, Myo1f, and Pttg1. New genes that failed to upregulate in the adult post-stroke brain include Bex4, Cd24, Klhl1/Mrp2, Trim67, and St8sia2. Among the upregulated pathways, the largest change was observed in the KEGG pathway "One carbon pool of folate," which is necessary for cellular proliferation, followed by the KEGG pathway "Antifolate resistance," whose genes mainly encode the family of ABC transporters responsible for the efflux of drugs that have entered the brain. We also noted three less-described downregulated KEGG pathways in experimental models: glycolipid biosynthesis, oxytocin, and cortisol pathways, which could be relevant as therapeutic targets. The limited brain plasticity of the adult brain is illustrated through molecular and histological analysis of the axonal growth factor, KIF4. Collectively, these results strongly suggest that further research is needed to decipher the complex genetic mechanisms that prevent the re-expression of brain plasticity-associated genes in the adult brain.

Novel insight into the etiology of ischemic stroke gained by integrative multiome-wide association study

Stroke, characterized by sudden neurological deficits, is the second leading cause of death worldwide. Although genome-wide association studies (GWAS) have successfully identified many genomic regions associated with ischemic stroke (IS), the genes underlying risk and their regulatory mechanisms remain elusive. Here, we integrate a large-scale GWAS (N = 1 296 908) for IS together with molecular QTLs data, including mRNA, splicing, enhancer RNA (eRNA), and protein expression data from up to 50 tissues (total N = 11 588). We identify 136 genes/eRNA/proteins associated with IS risk across 60 independent genomic regions and find IS risk is most enriched for eQTLs in arterial and brain-related tissues. Focusing on IS-relevant tissues, we prioritize 9 genes/proteins using probabilistic fine-mapping TWAS analyses. In addition, we discover that blood cell traits, particularly reticulocyte cells, have shared genetic contributions with IS using TWAS-based pheWAS and genetic correlation analysis. Lastly, we integrate our findings with a large-scale pharmacological database and identify a secondary bile acid, deoxycholic acid, as a potential therapeutic component. Our work highlights IS risk genes/splicing-sites/enhancer activity/proteins with their phenotypic consequences using relevant tissues as well as identify potential therapeutic candidates for IS.

Proteome-wide mendelian randomization identifies causal plasma proteins in venous thromboembolism development

Genome-wide association studies (GWAS) have identified numerous risk loci for venous thromboembolism (VTE), but it is challenging to decipher the underlying mechanisms. We employed an integrative analytical pipeline to transform genetic associations to identify novel plasma proteins for VTE. Proteome-wide association studies (PWAS) were determined by functional summary-based imputation leveraging data from a genome-wide association analysis (14,429 VTE patients, 267,037 controls), blood proteomes (1348 cases), followed by Mendelian randomization, Bayesian colocalization, protein-protein interaction, and pathway enrichment analysis. Twenty genetically regulated circulating protein abundances (F2, F11, ABO, PLCG2, LRP4, PLEK, KLKB1, PROC, KNG1, THBS2, SERPINA1, RARRES2, CEL, GP6, SERPINE2, SERPINA10, OBP2B, EFEMP1, F5, and MSR1) were associated with VTE. Of these 13 proteins demonstrated Mendelian randomized correlations. Six proteins (F2, F11, PLEK, SERPINA1, RARRES2, and SERPINE2) had strong support in colocalization analysis. Utilizing multidimensional data, this study suggests PLEK, SERPINA1, and SERPINE2 as compelling proteins that may provide key hints for future research and possible diagnostic and therapeutic targets for VTE.

Potential drug targets for asthma identified in the plasma and brain through Mendelian randomization analysis

Background

Asthma is a heterogeneous disease, and the involvement of neurogenic inflammation is crucial in its development. The standardized treatments focus on alleviating symptoms. Despite the availability of medications for asthma, they have proven to be inadequate in controlling relapses and halting the progression of the disease. Therefore, there is a need for novel drug targets to prevent asthma.

Methods

We utilized Mendelian randomization to investigate potential drug targets for asthma. We analyzed summary statistics from the UK Biobank and then replicated our findings in GWAS data by Demenais et al. and the FinnGen cohort. We obtained genetic instruments for 734 plasma and 73 brain proteins from recently reported GWAS. Next, we utilized reverse causal relationship analysis, Bayesian co-localization, and phenotype scanning as part of our sensitivity analysis. Furthermore, we performed a comparison and protein-protein interaction analysis to identify causal proteins. We also analyzed the possible consequences of our discoveries by the given existing asthma drugs and their targets.

Results

Using Mendelian randomization analysis, we identified five protein-asthma pairs that were significant at the Bonferroni level (P < 6.35 × 10-5). Specifically, in plasma, we found that an increase of one standard deviation in IL1R1 and ECM1 was associated with an increased risk of asthma, while an increase in ADAM19 was found to be protective. The corresponding odds ratios were 1.03 (95% CI, 1.02-1.04), 1.00 (95% CI, 1.00-1.01), and 0.99 (95% CI, 0.98-0.99), respectively. In the brain, per 10-fold increase in ECM1 (OR, 1.05; 95% CI, 1.03-1.08) and PDLIM4 (OR, 1.05; 95% CI, 1.04-1.07) increased the risk of asthma. Bayesian co-localization found that ECM1 in the plasma (coloc.abf-PPH4 = 0.965) and in the brain (coloc.abf-PPH4 = 0.931) shared the same mutation with asthma. The target proteins of current asthma medications were found to interact with IL1R1. IL1R1 and PDLIM4 were validated in two replication cohorts.

Conclusion

Our integrative analysis revealed that asthma risk is causally affected by the levels of IL1R1, ECM1, and PDLIM4. The results suggest that these three proteins have the potential to be used as drug targets for asthma, and further investigation through clinical trials is needed.

Identifying Novel Drug Targets for Epilepsy Through a Brain Transcriptome-Wide Association Study and Protein-Wide Association Study with Chemical-Gene-Interaction Analysis

Epilepsy is a severe neurological condition affecting 50-65 million individuals worldwide that can lead to brain damage. Nevertheless, the etiology of epilepsy remains poorly understood. Meta-analyses of genome-wide association studies involving 15,212 epilepsy cases and 29,677 controls of the ILAE Consortium cohort were used to conduct transcriptome-wide association studies (TWAS) and protein-wide association studies (PWAS). Furthermore, a protein-protein interaction (PPI) network was generated using the STRING database, and significant epilepsy-susceptible genes were verified using chip data. Chemical-related gene set enrichment analysis (CGSEA) was performed to determine novel drug targets for epilepsy. TWAS analysis identified 21,170 genes, of which 58 were significant (TWASfdr < 0.05) in ten brain regions, and 16 differentially expressed genes were verified based on mRNA expression profiles. The PWAS identified 2249 genes, of which 2 were significant (PWASfdr < 0.05). Through chemical-gene set enrichment analysis, 287 environmental chemicals associated with epilepsy were identified. We identified five significant genes (WIPF1, IQSEC1, JAM2, ICAM3, and ZNF143) that had causal relationships with epilepsy. CGSEA identified 159 chemicals that were significantly correlated with epilepsy (Pcgsea < 0.05), such as pentobarbital, ketone bodies, and polychlorinated biphenyl. In summary, we performed TWAS, PWAS (for genetic factors), and CGSEA (for environmental factors) analyses and identified several epilepsy-associated genes and chemicals. The results of this study will contribute to our understanding of genetic and environmental factors for epilepsy and may predict novel drug targets.

Proteome-Wide Association Studies for Blood Lipids and Comparison with Transcriptome-Wide Association Studies

Blood lipid traits are treatable and heritable risk factors for heart disease, a leading cause of mortality worldwide. Although genome-wide association studies (GWAS) have discovered hundreds of variants associated with lipids in humans, most of the causal mechanisms of lipids remain unknown. To better understand the biological processes underlying lipid metabolism, we investigated the associations of plasma protein levels with total cholesterol (TC), triglycerides (TG), high-density lipoprotein cholesterol (HDL), and low-density lipoprotein cholesterol (LDL) in blood. We trained protein prediction models based on samples in the Multi-Ethnic Study of Atherosclerosis (MESA) and applied them to conduct proteome-wide association studies (PWAS) for lipids using the Global Lipids Genetics Consortium (GLGC) data. Of the 749 proteins tested, 42 were significantly associated with at least one lipid trait. Furthermore, we performed transcriptome-wide association studies (TWAS) for lipids using 9,714 gene expression prediction models trained on samples from peripheral blood mononuclear cells (PBMCs) in MESA and 49 tissues in the Genotype-Tissue Expression (GTEx) project. We found that although PWAS and TWAS can show different directions of associations in an individual gene, 40 out of 49 tissues showed a positive correlation between PWAS and TWAS signed p-values across all the genes, which suggests a high-level consistency between proteome-lipid associations and transcriptome-lipid associations.

Transcriptome-wide association study reveals novel susceptibility genes for coronary atherosclerosis

Background

Genetic risk factors substantially contributed to the development of coronary atherosclerosis. Genome-wide association study (GWAS) has identified many risk loci for coronary atherosclerosis, but the translation of these loci into therapeutic targets is limited for their location in non-coding regions. Here, we aimed to screen the potential coronary atherosclerosis pathogenic genes expressed though TWAS (transcriptome wide association study) and explore the underlying mechanism association.

Methods

Four TWAS approaches (PrediXcan, JTI, UTMOST, and FUSION) were used to screen genes associated with coronary atherosclerosis. Enrichment analysis of TWAS-identified genes was applied through the Metascape website. The summary-data-based Mendelian randomization (SMR) analysis was conducted to provide the evidence of causal relationship between the candidate genes and coronary atherosclerosis. At last, the cell type-specific expression of the intersection genes was examined by using human coronary artery single-cell RNA-seq, interrogating the immune microenvironment of human coronary atherosclerotic plaque at different stages of maturity.

Results

We identified 19 genes by at least three approaches and 1 gene (NBEAL1) by four approaches. Enrichment analysis enriching the genes identified at least by two TWAS approaches, suggesting that these genes were markedly enriched in asthma and leukocyte mediated immunity reaction. Further, the summary-data-based Mendelian randomization (SMR) analysis provided the evidence of causal relationship between NBEAL1 gene and coronary atherosclerosis, confirming the protecting effects of NBEAL1 gene and coronary atherosclerosis. At last, the single cell cluster analysis demonstrated that NBEAL1 gene has differential expressions in macrophages, plasma cells and endothelial cells.

Conclusion

Our study identified the novel genes associated with coronary atherosclerosis and suggested the potential biological function for these genes, providing insightful guidance for further biological investigation and therapeutic approaches development in atherosclerosis-related diseases.

Novel insight into the etiology of ischemic stroke gained by integrative transcriptome-wide association study

Stroke, characterized by sudden neurological deficits, is the second leading cause of death worldwide. Although genome-wide association studies (GWAS) have successfully identified many genomic regions associated with ischemic stroke (IS), the genes underlying risk and their regulatory mechanisms remain elusive. Here, we integrate a large-scale GWAS (N=1,296,908) for IS together with mRNA, splicing, enhancer RNA (eRNA) and protein expression data (N=11,588) from 50 tissues. We identify 136 genes/eRNA/proteins associated with IS risk across 54 independent genomic regions and find IS risk is most enriched for eQTLs in arterial and brain-related tissues. Focusing on IS-relevant tissues, we prioritize 9 genes/proteins using probabilistic fine-mapping TWAS analyses. In addition, we discover that blood cell traits, particularly reticulocyte cells, have shared genetic contributions with IS using TWAS-based pheWAS and genetic correlation analysis. Lastly, we integrate our findings with a large-scale pharmacological database and identify a secondary bile acid, deoxycholic acid, as a potential therapeutic component. Our work highlights IS risk genes/splicing-sites/enhancer activity/proteins with their phenotypic consequences using relevant tissues as well as identify potential therapeutic candidates for IS.