Causal role of immune cell traits in stroke: A Mendelian randomization study.

OBJECTIVES: Immune mechanisms play a crucial role in the development of stroke. However, immune-related phenotypes are diverse and their associations with stroke are largely unknown. Here, we aimed to systematically explore the causal role of immune cell traits in stroke and its subtypes by leveraging data from genome-wide association studies (GWASs). MATERIALS AND METHODS: Exposure data were obtained from a recent GWAS on 731 immune cell traits profiled by flow cytometry involving 3757 individuals. By conducting two-sample univariable Mendelian randomization (MR) analyses, each immune cell trait was assessed for causal relationships with stroke outcomes from the MEGASTROKE Consortium (40,585 cases and 406,111 controls). The robustness of the MR results was verified by a series of sensitivity analyses. RESULTS: We identified three significant associations after Bonferroni correction (P < 1.37E-05). Increased CD27 expression on memory B cell (OR = 1.23, 95% CI = 1.14-1.33, P = 2.78E-08), IgD-CD38dim B cell (OR = 1.16, 95% CI = 1.09-1.23, P = 5.98E-06) and unswitched memory B cell (OR = 1.18, 95% CI = 1.10-1.27, P = 1.09E-05) were associated with a higher risk of large-artery atherosclerotic stroke (LAS). Furthermore, expression quantitative trait loci data also indicated elevated blood CD27 mRNA level was a risk factor for LAS (OR = 1.37, 95% CI = 1.02-1.84, P = 0.037). CONCLUSIONS: This study provided genetic evidence of the causal relationship between immune cell traits and stroke, highlighting the role of CD27 on memory B cell as a novel factor for LAS risk.

Genetic and Epigenetic Analysis Revealing Variants in the NCAM1-TTC12-ANKK1-DRD2 Cluster Associated Significantly With Nicotine Dependence in Chinese Han Smokers.

BACKGROUNDS: Although studies have demonstrated that the NCAM1-TTC12-ANKK1-DRD2 gene cluster plays essential roles in addictions in subjects of European and African origin, study of Chinese Han subjects is limited. Further, the underlying biological mechanisms of detected associations are largely unknown. METHODS: Sixty-four single-nucleotide polymorphisms (SNPs) in this cluster were analyzed for association with Fagerstrom Test for Nicotine Dependence score (FTND) and cigarettes per day (CPD) in male Chinese Han smokers (N = 2616). Next-generation bisulfite sequencing was used to discover smoking-associated differentially methylated regions (DMRs). Both cis-eQTL and cis-mQTL analyses were applied to assess the cis-regulatory effects of these risk SNPs. RESULTS: Association analysis revealed that rs4648317 was significantly associated with FTND and CPD (p = .00018; p = .00072). Moreover, 14 additional SNPs were marginally significantly associated with FTND or CPD (p = .05-.01). Haplotype-based association analysis showed that one haplotype in DRD2, C-T-A-G, formed by rs4245148, rs4581480, rs4648317, and rs11214613, was significantly associated with CPD (p = .0005) and marginally associated with FTND (p = .003). Further, we identified four significant smoking-associated DMRs, three of which are located in the DRD2/ANKK1 region (p = .0012-.00005). Finally, we found five significant CpG-SNP pairs (p = 7.9 × 10-9-6.6 × 10-6) formed by risk SNPs rs4648317, rs11604671, and rs2734849 and three methylation loci. CONCLUSIONS: We found two missense variants (rs11604671; rs2734849) and an intronic variant (rs4648317) with significant effects on ND and further explored their mechanisms of action through expression and methylation analysis. We found the majority of smoking-related DMRs are located in the ANKK1/DRD2 region, indicating a likely causative relation between non-synonymous SNPs and DMRs. IMPLICATIONS: This study shows that there exist significant association of variants and haplotypes in ANKK1/DRD2 region with ND in Chinese male smokers. Further, this study also shows that DNA methylation plays an important role in mediating such associations.

Single-Cell Transcriptome Analysis Identified Core Genes and Transcription Factors in Mesenchymal Cell Differentiation during Liver Cirrhosis.

BACKGROUND: Mesenchymal cells, including hepatic stellate cells (HSCs), fibroblasts (FBs), myofibroblasts (MFBs), and vascular smooth muscle cells (VSMCs), are the main cells that affect liver fibrosis and play crucial roles in maintaining tissue homeostasis. The dynamic evolution of mesenchymal cells is very important but remains to be explored for researching the reversible mechanism of hepatic fibrosis and its evolution mechanism of hepatic fibrosis to cirrhosis. METHODS: Here, we analysed the transcriptomes of more than 50,000 human single cells from three cirrhotic and three healthy liver tissue samples and the mouse hepatic mesenchymal cells of two healthy and two fibrotic livers to reconstruct the evolutionary trajectory of hepatic mesenchymal cells from a healthy to a cirrhotic state, and a subsequent integrative analysis of bulk RNA sequencing (RNA-seq) data of HSCs from quiescent to active (using transforming growth factor ß1 (TGF-ß1) to stimulate LX-2) to inactive states. RESULTS: We identified core genes and transcription factors (TFs) involved in mesenchymal cell differentiation. In healthy human and mouse livers, the expression of NR1H4 and members of the ZEB families (ZEB1 and ZEB2) changed significantly with the differentiation of FB into HSC and VSMC. In cirrhotic human livers, VSMCs transformed into HSCs with downregulation of MYH11, ACTA2, and JUNB and upregulation of PDGFRB, RGS5, IGFBP5, CD36, A2M, SOX5, and MEF2C. Following HSCs differentiation into MFBs with the upregulation of COL1A1, TIMP1, and NR1H4, a small number of MFBs reverted to inactivated HSCs (iHSCs). The differentiation trajectory of mouse hepatic mesenchymal cells was similar to that in humans; however, the evolution trajectory and proportion of cell subpopulations that reverted from MFBs to iHSCs suggest that the mouse model may not accurately reflect disease progression and outcome in humans. CONCLUSIONS: Our analysis elucidates primary genes and TFs involved in mesenchymal cell differentiation during liver fibrosis using scRNA-seq data, and demonstrated the core genes and TFs in process of HSC activation to MFB and MFB reversal to iHSC using bulk RNA-seq data of human fibrosis induced by TGF-ß1. Furthermore, our findings suggest promising targets for the treatment of liver fibrosis and provide valuable insights into the molecular mechanisms underlying its onset and progression.

Hsa_circ_0007401 regulates gemcitabine resistance of pancreatic cancer through the hsa-miR-6509-3p/fli1 axis.

Pancreatic cancer is a highly malignant cancer with a poor prognosis. Owing to the strong drug resistance of pancreatic cancer, adjuvant chemotherapy has failed to achieve good results in clinical practice. The expression profile data of circular RNA (circRNA) (GSE110580), microRNA (miRNA) (GSE79234), and messenger RNA (mRNA) (GSE140077, GES35141) were obtained from the gene expression omnibus database. The Cancer-Specific circRNA Database identified the structural pattern of circRNA, and the starBase and circBank databases together predicted the miRNA of circRNA. The mirDIP database predicts the target mRNAs of miRNAs and identifies the ceRNA network of circRNA-miRNA-mRNA via negative regulatory mechanisms. The final validation was performed using clinical data from the cancer treatment response gene signature database of patients treated with gemcitabine for pancreatic cancer of the cancer genome atlas. By differential expression analysis, 22 differential circRNAs (8 upregulated and 14 downregulated), 70 differential microRNAs (37 upregulated and 33 downregulated), and 256 differential messenger RNA (DEmRNA) (161 upregulated and 95 downregulated) were obtained. Gene ontology and Kyoto encyclopedia of genes and genomes enrichment analyses showed that DEmRNAs were associated with drug response, exogenous cellular stimulation, and the tumor necrosis factor signaling pathway. The screened downregulated differential circular RNA (hsa_circ_0007401), upregulated differential microRNA (hsa-miR-6509-3p), and downregulated DEmRNA (FLI1) were consistent with the negative regulation mechanism of the ceRNA network, and FLI1 was significantly downregulated in the data of gemcitabine-resistant pancreatic cancer patients in the cancer genome atlas (n = 26).

Integrated analysis of circRNA-associated ceRNA network in ischemic stroke.

Introduction: Stroke, of which ischemic stroke (IS) is the major type, is the second leading cause of disability and death worldwide. Circular RNAs (circRNAs) are reported to play important role in the physiology and pathology of IS. CircRNAs often act as competing endogenous RNA (ceRNA) to regulate gene expression by acting as miRNA sponges. However, whole transcriptome-wide screenings of circRNA-mediated ceRNA networks associated with IS are still lacking. In the present study, we constructed a circRNA-miRNA-mRNA ceRNA network by whole transcriptome-wide analysis. Methods: CircRNAs, miRNAs and mRNAs expression profiles were downloaded from the Gene Expression Omnibus (GEO) datasets. We identified differentially expressed (DE) circRNAs, miRNAs, and mRNAs in IS patients. StarBase and CircBank databases were used to predict the miRNA targets of DEcircRNAs, and mirDIP database was used to predict the mRNA targets of DEmiRNAs. CircRNA-miRNA pairs and miRNA-mRNA pairs were established. Then, we identified hub genes via protein-protein interaction analysis and constructed a core ceRNA sub-network. Results: In total, 276 DEcircRNAs, 43 DEmiRNAs, and 1926 DEmRNAs were explored. The ceRNA network included 69 circRNAs, 24 miRNAs, and 92 mRNAs. The core ceRNA subnetwork included hsa_circ_0011474, hsa_circ_0023110, CDKN1A, FHL2, RPS2, CDK19, KAT6A, CBX1, BRD4, and ZFHX3. Discussion: In conclusion, we established a novel hsa_circ_0011474 - hsa-miR-20a-5p/hsa-miR-17-5p - CDKN1A ceRNA regulatory axis associated with IS. Our findings provide new insights into the pathogenesis of IS and offer promising diagnostic and predictive biomarkers.

Risk factors for stroke outcomes in adults: Stroke in China.

This study aimed to determine the possible risk factors for stroke outcomes based on prospective cohort study in China. A total of 146 stroke patients were recruited and divided into 2 groups, which assessed using the modified Rankin Scale (mRS), good outcomes (mRS <= 2) and poor outcomes group (mRS > 2). Demographic, clinical, and laboratory characteristics of participants were obtained from the medical record. The multivariable logistic regression analysis was employed to assess the risk factors for stroke outcomes. Of 146 participants, 28 (19.18%) were presented with poor outcomes at day 90. As a result of multivariable logistic regression analysis, a significantly increased risk of stroke outcomes was found in patients with Barthel Index (BI) score (stroke (OR 1.50, 95% CI 1.21 ~ 1.85, P < .001) and IS (OR 1.48, 95% CI 1.20 ~ 1.83, P < .001)).

A-GSTCN: An Augmented Graph Structural-Temporal Convolution Network for Medication Recommendation Based on Electronic Health Records.

Medication recommendation based on electronic health records (EHRs) is a significant research direction in the biomedical field, which aims to provide a reasonable prescription for patients according to their historical and current health conditions. However, the existing recommended methods have many limitations in dealing with the structural and temporal characteristics of EHRs. These methods either only consider the current state while ignoring the historical situation, or fail to adequately assess the structural correlations among various medical events. These factors result in poor recommendation quality. To solve this problem, we propose an augmented graph structural-temporal convolutional network (A-GSTCN). Firstly, an augmented graph attention network is used to model the structural features among medical events of patients' EHRs. Next, the dilated convolution combined with residual connection is applied in the proposed model, which can improve the temporal prediction capability and further reduce the complexity. Moreover, the cache memory module further enhances the model's learning of the history of EHRs. Finally, the A-GSTCN model is compared with the baselines through experiments, and the efficiency of the A-GSTCN model is verified by Jaccard, F1 and PRAUC. Not only that, the proposed model also reduces the training parameters by an order of magnitude.

Genome-wide methylation and expression analyses reveal the epigenetic landscape of immune-related diseases for tobacco smoking.

BACKGROUND: Smoking is a major causal risk factor for lung cancer, chronic obstructive pulmonary disease (COPD), cardiovascular disease (CVD), and is the main preventable cause of deaths in the world. The components of cigarette smoke are involved in immune and inflammatory processes, which may increase the prevalence of cigarette smoke-related diseases. However, the underlying molecular mechanisms linking smoking and diseases have not been well explored. This study was aimed to depict a global map of DNA methylation and gene expression changes induced by tobacco smoking and to explore the molecular mechanisms between smoking and human diseases through whole-genome bisulfite sequencing (WGBS) and RNA-sequencing (RNA-seq). RESULTS: We performed WGBS on 72 samples (36 smokers and 36 nonsmokers) and RNA-seq on 75 samples (38 smokers and 37 nonsmokers), and cytokine immunoassay on plasma from 22 males (9 smokers and 13 nonsmokers) who were recruited from the city of Jincheng in China. By comparing the data of the two groups, we discovered a genome-wide methylation landscape of differentially methylated regions (DMRs) associated with smoking. Functional enrichment analyses revealed that both smoking-related hyper-DMR genes (DMGs) and hypo-DMGs were related to synapse-related pathways, whereas the hypo-DMGs were specifically related to cancer and addiction. The differentially expressed genes (DEGs) revealed by RNA-seq analysis were significantly enriched in the "immunosuppression" pathway. Correlation analysis of DMRs with their corresponding gene expression showed that genes affected by tobacco smoking were mostly related to immune system diseases. Finally, by comparing cytokine concentrations between smokers and nonsmokers, we found that vascular endothelial growth factor (VEGF) was significantly upregulated in smokers. CONCLUSIONS: In sum, we found that smoking-induced DMRs have different distribution patterns in hypermethylated and hypomethylated areas between smokers and nonsmokers. We further identified and verified smoking-related DMGs and DEGs through multi-omics integration analysis of DNA methylome and transcriptome data. These findings provide us a comprehensive genomic map of the molecular changes induced by smoking which would enhance our understanding of the harms of smoking and its relationship with diseases.

Genome-wide DNA methylation analysis reveals significant impact of long-term ambient air pollution exposure on biological functions related to mitochondria and immune response.

Exposure to long-term ambient air pollution is believed to have adverse effects on human health. However, the mechanisms underlying these impacts are poorly understood. DNA methylation, a crucial epigenetic modification, is susceptible to environmental factors and likely involved in these processes. We conducted a whole-genome bisulfite sequencing study on 120 participants from a highly polluted region (HPR) and a less polluted region (LPR) in China, where the HPR had much higher concentrations of five air pollutants (PM2.5, PM10, SO2, NO2, and CO) (fold difference 1.6 to 6.6 times; P value 1.80E-07 to 3.19E-23). Genome-wide methylation analysis revealed 371 DMRs in subjects from the two areas and these DMRs were located primarily in gene regulatory elements such as promoters and enhancers. Gene enrichment analysis showed that DMR-related genes were significantly enriched in diseases related to pulmonary disorders and cancers and in biological processes related to mitochondrial assembly and cytokine production. Further, HPR participants showed a higher mtDNA copy number. Of those identified DMRs, 15 were significantly correlated with mtDNA copy number. Finally, cytokine assay indicated that an increased plasma interleukin-5 level was associated with greater air pollution. Taken together, our findings suggest that exposure to long-term ambient air pollution can lead to alterations in DNA methylation whose functions relate to mitochondria and immune responses.

Identification of 34 genes conferring genetic and pharmacological risk for the comorbidity of schizophrenia and smoking behaviors.

The prevalence of smoking is significantly higher in persons with schizophrenia (SCZ) than in the general population. However, the biological mechanisms of the comorbidity of smoking and SCZ are largely unknown. This study aimed to reveal shared biological pathways for the two diseases by analyzing data from two genome-wide association studies with a total sample size of 153,898. With pathway-based analysis, we first discovered 18 significantly enriched pathways shared by SCZ and smoking, which were classified into five groups: postsynaptic density, cadherin binding, dendritic spine, long-term depression, and axon guidance. Then, by using an integrative analysis of genetic, epigenetic, and expression data, we found not only 34 critical genes (e.g., PRKCZ, ARHGEF3, and CDKN1A) but also various risk-associated SNPs in these genes, which convey susceptibility to the comorbidity of the two disorders. Finally, using both in vivo and in vitro data, we demonstrated that the expression profiles of the 34 genes were significantly altered by multiple psychotropic drugs. Together, this multi-omics study not only reveals target genes for new drugs to treat SCZ but also reveals new insights into the shared genetic vulnerabilities of SCZ and smoking behaviors.

Association and cis-mQTL analysis of variants in CHRNA3-A5, CHRNA7, CHRNB2, and CHRNB4 in relation to nicotine dependence in a Chinese Han population.

Nicotine dependence (ND) is a worldwide health problem. Numerous genetic studies have demonstrated a significant association of variants in nicotinic acetylcholine receptors (nAChRs) with smoking behaviors. However, most of these studies enrolled only subjects of European or African ancestry. In addition, although an increasing body of evidence implies a causal connection of single-nucleotide polymorphisms (SNPs) and epigenetic regulation of gene expression, few studies of this issue have been reported. In this study, we performed both association and interaction analysis for 67 SNPs in CHRNA3-A5, CHRNA7, CHRNB2, and CHRNB4 with ND in a Chinese Han population (N = 5055). We further analyzed cis-mQTL for the three most significant SNPs and 5580 potential methylation loci within these target gene regions. Our results indicated that the SNPs rs1948 and rs7178270 in CHRNB4 and rs3743075 in CHRNA3 were significantly associated with the Fagerström Test for Nicotine Dependence (FTND) score (p = 6.6 × 10-5; p = 2.0 × 10-4, and p = 7.0 × 10-4, respectively). Haplotype-based association analysis revealed that two major haplotypes, T-G and C-A, formed by rs3743075-rs3743074 in CHRNA3, and other two major haplotypes, A-G-C and G-C-C, formed by rs1948-rs7178270-rs17487223 in CHRNB4, were significantly associated with the FTND score (p ≤ 8.0 × 10-4). Further, we found evidence for the presence of significant interaction among variants within CHRNA3/B4/A5, CHRNA4/B2/A5, and CHRNA7 in affecting ND, with corresponding p values of 5.8 × 10-6, 8.0 × 10-5, and 0.012, respectively. Finally, we identified two CpG sites (CpG_2975 and CpG_3007) in CHRNA3 that are significantly associated with three cis-mQTL SNPs (rs1948, rs7178270, rs3743075) in the CHRNA5/A3/B4 cluster (p ≤ 1.9 × 10-6), which formed four significant CpG-SNP pairs in our sample. Together, we revealed at least three novel SNPs in CHRNA3 and CHRNB4 to be significantly associated with the FTND score. Further, we showed that these significant variants contribute to ND via two methylated sites, and we demonstrated significant interaction affecting ND among variants in CHRNA5/A3/B4, CHRNA7, and CHRNA4/B2/A5. In sum, these findings provide robust evidence that SNPs in nAChR genes convey a risk of ND in the Chinese Han population.

Meta-Analysis of the COMT Val158Met Polymorphism in Major Depressive Disorder: Effect of Ethnicity.

The COMT (catechol-O-methyltransferase) Val158Met polymorphism (rs4680) is a potential susceptibility variant for major depressive disorder (MDD). Although many genetic studies have examined the association between MDD and this polymorphism, the results were inconclusive. In the present study, we conducted a series of meta-analyses of samples consisting of 2905 MDD cases and 2403 controls with the goal of determining whether this variant indeed has any effect on MDD. We revealed a significant association in the comparison of Val/Val + Val/Met vs. Met/Met (OR =1.180; 95 % CI = 1.019, 1.367; P = 0.027), Val/Met vs. Val/Val (OR =1.18; 95 % CI = 1.038, 1.361; P = 0.013), and Val/Met vs. Met/Met (OR =1.229; 95 % CI = 1.053, 1.435; P = 0.009). Further meta-analyses of samples with European ancestry demonstrated a significant association of this SNP with MDD susceptibility in Val/Val + Val/Met vs. Met/Met (OR =1.231, 95 % CI = 1.046, 1.449; P = 0.013) and Val/Met vs. Met/Met (OR =1.284, 95 % CI = 1.050, 1.484; P = 0.012). For the samples with East Asian ancestry, we found a significant association in both allelic (Val vs. Met: OR =0.835; 95 % CI = 0.714, 0.975; P = 0.023) and genotypic (Met/Met + Val/Met vs. Val/Val: OR =1.431, 95 % CI = 1.143, 1.791; P = 0.002; Val/Met vs. Val/Val: OR =1.482, 95 % CI = 1.171, 1.871; P = 0.001) analyses. No evidence of heterogeneity among studies or publication bias was observed. Together, our results indicate that the COMT Val158Met polymorphism is a vulnerability factor for MDD with distinct effects in different ethnic populations.