Comprehensive Analysis of Differential Immunocyte Infiltration and Potential ceRNA Networks Involved in the Development of Atrial Fibrillation

Identification of atrial fibrillation-related genes through transcriptome data analysis and Mendelian randomization

Background

Atrial fibrillation (AF) is a common persistent arrhythmia characterized by rapid and chaotic atrial electrical activity, potentially leading to severe complications such as thromboembolism, heart failure, and stroke, significantly affecting patient quality of life and safety. As the global population ages, the prevalence of AF is on the rise, placing considerable strains on individuals and healthcare systems. This study utilizes bioinformatics and Mendelian Randomization (MR) to analyze transcriptome data and genome-wide association study (GWAS) summary statistics, aiming to identify biomarkers causally associated with AF and explore their potential pathogenic pathways.

Methods

We obtained AF microarray datasets GSE41177 and GSE79768 from the Gene Expression Omnibus (GEO) database, merged them, and corrected for batch effects to pinpoint differentially expressed genes (DEGs). We gathered exposure data from expression quantitative trait loci (eQTL) and outcome data from AF GWAS through the IEU Open GWAS database. We employed inverse variance weighting (IVW), MR-Egger, weighted median, and weighted model approaches for MR analysis to assess exposure-outcome causality. IVW was the primary method, supplemented by other techniques. The robustness of our results was evaluated using Cochran's Q test, MR-Egger intercept, MR-PRESSO, and leave-one-out sensitivity analysis. A "Veen" diagram visualized the overlap of DEGs with significant eQTL genes from MR analysis, referred to as common genes (CGs). Additional analyses, including Gene Ontology (GO) enrichment, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, and immune cell infiltration studies, were conducted on these intersecting genes to reveal their roles in AF pathogenesis.

Results

The combined dataset revealed 355 differentially expressed genes (DEGs), with 228 showing significant upregulation and 127 downregulated. Mendelian randomization (MR) analysis identified that the autocrine motility factor receptor (AMFR) [IVW: OR = 0.977; 95% CI, 0.956-0.998; P = 0.030], leucine aminopeptidase 3 (LAP3) [IVW: OR = 0.967; 95% CI, 0.934-0.997; P = 0.048], Rab acceptor 1 (RABAC1) [IVW: OR = 0.928; 95% CI, 0.875-0.985; P = 0.015], and tryptase beta 2 (TPSB2) [IVW: OR = 0.971; 95% CI, 0.943-0.999; P = 0.049] are associated with a reduced risk of atrial fibrillation (AF). Conversely, GTPase-activating SH3 domain-binding protein 2 (G3BP2) [IVW: OR = 1.030; 95% CI, 1.004-1.056; P = 0.024], integrin subunit beta 2 (ITGB2) [IVW: OR = 1.050; 95% CI, 1.017-1.084; P = 0.003], glutaminyl-peptide cyclotransferase (QPCT) [IVW: OR = 1.080; 95% CI, 1.010-0.997; P = 1.154], and tripartite motif containing 22 (TRIM22) [IVW: OR = 1.048; 95% CI, 1.003-1.095; P = 0.035] are positively associated with AF risk. Sensitivity analyses indicated a lack of heterogeneity or horizontal pleiotropy (P > 0.05), and leave-one-out analysis did not reveal any single nucleotide polymorphisms (SNPs) impacting the MR results significantly. GO and KEGG analyses showed that CG is involved in processes such as protein polyubiquitination, neutrophil degranulation, specific and tertiary granule formation, protein-macromolecule adaptor activity, molecular adaptor activity, and the SREBP signaling pathway, all significantly enriched. The analysis of immune cell infiltration demonstrated associations of CG with various immune cells, including plasma cells, CD8T cells, resting memory CD4T cells, regulatory T cells (Tregs), gamma delta T cells, activated NK cells, activated mast cells, and neutrophils.

Conclusion

By integrating bioinformatics and MR approaches, genes such as AMFR, G3BP2, ITGB2, LAP3, QPCT, RABAC1, TPSB2, and TRIM22 are identified as causally linked to AF, enhancing our understanding of its molecular foundations. This strategy may facilitate the development of more precise biomarkers and therapeutic targets for AF diagnosis and treatment.

CDK5R1, GSE1, HSPG2 and WDFY3 as indirect epigenetic-sensitive genes in atrial fibrillation

BACKGROUND

Although mounting evidence supports that aberrant DNA methylation occurs in the hearts of patients with atrial fibrillation (AF), noninvasive epigenetic characterization of AF has not yet been defined.

METHODS

We investigated DNA methylome changes in peripheral blood CD4+ T cells isolated from 10 patients with AF relative to 11 healthy subjects (HS) who were enrolled in the DIANA clinical trial (NCT04371809) via reduced-representation bisulfite sequencing (RRBS).

RESULTS

An atrial-specific PPI network revealed 18 hub differentially methylated genes (DMGs), wherein ROC curve analysis revealed reasonable diagnostic performance of DNA methylation levels found within CDK5R1 (AUC = 0.76; p = 0.049), HSPG2 (AUC = 0.77; p = 0.038), WDFY3 (AUC = 0.78; p = 0.029), USP49 (AUC = 0.76; p = 0.049), GSE1 (AUC = 0.76; p = 0.049), AIFM1 (AUC = 0.76; p = 0.041), CDK5RAP2 (AUC = 0.81; p = 0.017), COL4A1 (AUC = 0.86; p < 0.001), SEPT8 (AUC = 0.90; p < 0.001), PFDN1 (AUC = 0.90; p < 0.01) and ACOT7 (AUC = 0.78; p = 0.032). Transcriptional profiling of the hub DMGs provided a significant overexpression of PSDM6 (p = 0.004), TFRC (p = 0.01), CDK5R1 (p < 0.001), HSPG2 (p = 0.01), WDFY3 (p < 0.001), USP49 (p = 0.004) and GSE1 (p = 0.021) in AF patients vs HS.

CONCLUSIONS

CDK5R1, GSE1, HSPG2 and WDFY3 resulted the best discriminatory genes both at methylation and gene expression level. Our results provide several candidate diagnostic biomarkers with the potential to advance precision medicine in AF.

Integrative identification of immune-related key genes in atrial fibrillation using weighted gene coexpression network analysis and machine learning

Background

The immune system significantly participates in the pathologic process of atrial fibrillation (AF). However, the molecular mechanisms underlying this participation are not completely explained. The current research aimed to identify critical genes and immune cells that participate in the pathologic process of AF.

Methods

CIBERSORT was utilized to reveal the immune cell infiltration pattern in AF patients. Meanwhile, weighted gene coexpression network analysis (WGCNA) was utilized to identify meaningful modules that were significantly correlated with AF. The characteristic genes correlated with AF were identified by the least absolute shrinkage and selection operator (LASSO) logistic regression and support vector machine recursive feature elimination (SVM-RFE) algorithm.

Results

In comparison to sinus rhythm (SR) individuals, we observed that fewer activated mast cells and regulatory T cells (Tregs), as well as more gamma delta T cells, resting mast cells, and M2 macrophages, were infiltrated in AF patients. Three significant modules (pink, red, and magenta) were identified to be significantly associated with AF. Gene enrichment analysis showed that all 717 genes were associated with immunity- or inflammation-related pathways and biological processes. Four hub genes (GALNT16, HTR2B, BEX2, and RAB8A) were revealed to be significantly correlated with AF by the SVM-RFE algorithm and LASSO logistic regression. qRT–PCR results suggested that compared to the SR subjects, AF patients exhibited significantly reduced BEX2 and GALNT16 expression, as well as dramatically elevated HTR2B expression. The AUC measurement showed that the diagnostic efficiency of BEX2, HTR2B, and GALNT16 in the training set was 0.836, 0.883, and 0.893, respectively, and 0.858, 0.861, and 0.915, respectively, in the validation set.

Conclusions

Three novel genes, BEX2, HTR2B, and GALNT16, were identified by WGCNA combined with machine learning, which provides potential new therapeutic targets for the early diagnosis and prevention of AF.

Identification of Potential Key Biomarkers of Atrial Fibrillation and Their Correlation with Immune Infiltration in Atrial Tissue

Objective

To identify potential key biomarkers and characterize immune infiltration in atrial tissue of patients with atrial fibrillation (AF) through bioinformatics analysis.

Methods

Differentially expressed genes (DEGs) were identified by the LIMMA package in Bioconductor, and functional and pathway enrichment analyses were undertaken using GO and KEGG. The LASSO logistic regression and BORUTA algorithm were employed to screen for potential novel key markers of AF from all DEGs. Gene set variation analysis was also performed. Single-sample gene set enrichment analysis was employed to quantify the infiltration levels for each immune cell type, and the correlation between hub genes and infiltrating immune cells was analyzed.

Results

A total of 52 DEGs were identified, including of 26 downregulated DEGs and 26 upregulated DEGs. DEGs were primarily enriched in the Major Histocompatibility Complex class II protein complex, glucose homeostasis, protein tetramerization, regulation of synapse organization, cytokine activity, heart morphogenesis, and blood circulation. Three downregulated genes and three upregulated genes were screened by LASSO logistic regression and the BORUTA algorithm. Finally, immune infiltration analysis indicated that the atrial tissue of AF patients contained significant infiltration of APC_co_inhibition, Mast_cell, neutrophils, pDCs, T_cell_costimulation, and Th1_cells compared with paired sinus rhythm (SR) atrial tissue, and the three downregulated genes were negatively correlated with the six kinds of immune cells mentioned above.

Conclusion

The hub genes identified in this study and the differences in immune infiltration of atrial tissue observed between AF and SR tissue might help to characterize the occurrence and progression of AF.

Construction and Analysis of the lncRNA-miRNA-mRNA Network Based on Competing Endogenous RNA in Atrial Fibrillation

Background

Accumulated studies have revealed that long non-coding RNAs (lncRNAs) play critical roles in human diseases by acting as competing endogenous RNAs (ceRNAs). However, functional roles and regulatory mechanisms of lncRNA-mediated ceRNA in atrial fibrillation (AF) remain unknown. In the present study, we aimed to construct the lncRNA-miRNA-mRNA network based on ceRNA theory in AF by using bioinformatic analyses of public datasets.

Methods

Microarray data sets of GSE115574 and GSE79768 from the Gene Expression Omnibus database were downloaded. Twenty-one AF right atrial appendage (RAA) samples and 22 sinus rhythm (SR) subjects RAA samples were selected for subsequent analyses. After merging all microarray data and adjusting for batch effect, differentially expressed genes were identified. Gene Ontology (GO) categories and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were carried out. A ceRNA network was constructed.

Result

A total of 8 lncRNAs and 43 mRNAs were significantly differentially expressed with fold change >1.5 (p < 0.05) in RAA samples of AF patients when compared with SR. GO and KEGG pathway analysis showed that cardiac muscle contraction pathway were involved in AF development. The ceRNA was predicted by co-expressing LOC101928304/ LRRC2 from the constructional network analysis, which was competitively combined with miR-490-3p. The expression of LOC101928304 and LRRC were up-regulated in myocardial tissue of patients with AF, while miR-490-3p was down-regulated.

Conclusion

We constructed the LOC101928304/miR-490-3p/LRRC2 network based on ceRNA theory in AF in the bioinformatic analyses of public datasets. The ceRNA network found from this study may help improve our understanding of lncRNA-mediated ceRNA regulatory mechanisms in the pathogenesis of AF.

Molecular Insights in Atrial Fibrillation Pathogenesis and Therapeutics: A Narrative Review

The prevalence of atrial fibrillation (AF) is bound to increase globally in the following years, affecting the quality of life of millions of people, increasing mortality and morbidity, and beleaguering health care systems. Increasingly effective therapeutic options against AF are the constantly evolving electroanatomic substrate mapping systems of the left atrium (LA) and ablation catheter technologies. Yet, a prerequisite for better long-term success rates is the understanding of AF pathogenesis and maintenance. LA electrical and anatomical remodeling remains in the epicenter of current research for novel diagnostic and treatment modalities. On a molecular level, electrical remodeling lies on impaired calcium handling, enhanced inwardly rectifying potassium currents, and gap junction perturbations. In addition, a wide array of profibrotic stimuli activates fibroblast to an increased extracellular matrix turnover via various intermediaries. Concomitant dysregulation of the autonomic nervous system and the humoral function of increased epicardial adipose tissue (EAT) are established mediators in the pathophysiology of AF. Local atrial lymphomononuclear cells infiltrate and increased inflammasome activity accelerate and perpetuate arrhythmia substrate. Finally, impaired intracellular protein metabolism, excessive oxidative stress, and mitochondrial dysfunction deplete atrial cardiomyocyte ATP and promote arrhythmogenesis. These overlapping cellular and molecular alterations hinder us from distinguishing the cause from the effect in AF pathogenesis. Yet, a plethora of therapeutic modalities target these molecular perturbations and hold promise in combating the AF burden. Namely, atrial selective ion channel inhibitors, AF gene therapy, anti-fibrotic agents, AF drug repurposing, immunomodulators, and indirect cardiac neuromodulation are discussed here.

In Silico Identification of miRNA-lncRNA Interactions in Male Reproductive Disorder Associated with COVID-19 Infection

Coronavirus disease 2019 (COVID-19), a global pandemic, is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Angiotensin-converting enzyme 2 (ACE2) is the receptor for SARS-CoV-2 and transmembrane serine protease 2 (TMPRSS2) facilitates ACE2-mediated virus entry. Moreover, the expression of ACE2 in the testes of infertile men is higher than normal, which indicates that infertile men may be susceptible to be infected and SARS-CoV-2 may cause reproductive disorder through the pathway induced by ACE2 and TMPRSS2. Little is known about the pathway regulation of ACE2 and TMPRSS2 expression in male reproductive disorder. Since the regulation of gene expression is mediated by microRNAs (miRNAs) and long non-coding RNAs (lncRNAs) at the post-transcriptional level, the aim of this study was to analyze the dysregulated miRNA-lncRNA interactions of ACE2 and TMPRSS2 in male reproductive disorder. Using bioinformatics analysis, we speculate that the predicted miRNAs including miR-125a-5p, miR-125b-5p, miR-574-5p, and miR-936 as regulators of ACE2 and miR-204-5p as a modulator of TMPRSS2 are associated with male infertility. The lncRNAs with a tissue-specific expression for testis including GRM7-AS3, ARHGAP26-AS1, BSN-AS1, KRBOX1-AS1, CACNA1C-IT3, AC012361.1, FGF14-IT1, AC012494.1, and GS1-24F4.2 were predicted. The identified miRNAs and lncRNAs are proposed as potential biomarkers to study the possible association between COVID-19 and male infertility. This study encourages further studies of miRNA-lncRNA interactions to explain the molecular mechanisms of male infertility in COVID-19 patients.