Identification of key genes in atrial fibrillation using bioinformatics analysis
BMC Cardiovascular Disorders volume 20, Article number: 363 (2020)
Atrial fibrillation (AF) is one of the most common arrhythmia, which brings huge burden to the individual and the society. However, the mechanism of AF is not clear. This paper aims at screening the key differentially expressed genes (DEGs) of atrial fibrillation and to construct enrichment analysis and protein-protein interaction (PPI) network analysis for these DEGs.
The datasets were collected from the Gene Expression Omnibus database to extract data of left atrial appendage (LAA) RNA of patients with or without AF in GSE79768, GSE31821, GSE115574, GSE14975 and GSE41177. Batch normalization, screening of the differential genes and gene ontology analysis were finished by R software. Reactome analysis was used for pathway analysis. STRING platform was utilized for PPI network analysis. At last, we performed reverse transcription-quantitative polymerase chain reaction (RT-qPCR) to validate the expression of key genes in 20 sinus rhythm (SR) LAA tissues and 20 AF LAA tissues.
A total of 106 DEGs were screened in the merged dataset. Among these DEGs, 74 genes were up-regulated and 32 genes down-regulated. DEGs were mostly enriched in extracellular matrix organization, protein activation cascade and extracellular structure organization. In PPI network, we identified SPP1, COL5A1 and VCAN as key genes which were associated with extracellular matrix. RT-qPCR showed the same expression trend of the three key genes as in our bioinformatics analysis. The expression levels of SPP1, COL5A1 and VCAN were increased in AF tissues compared to SR tissues (P < 0.05).
According to the analyses which were conducted by bioinformatics tools, genes related to extracellular matrix were involved in pathology of AF and may become the possible targets for the diagnosis and treatment of AF.
Atrial fibrillation (AF) is one of the most common arrhythmia in the world and the incidence of AF rises as the age increases [1,2,3,4]. AF could induce palpitation, heart failure and thrombus formation [5, 6] which may impair the quality of life  and increase the risk of mortality [8, 9]. Since the mechanism of AF is not clear, it is important to reveal the mechanism of pathogenesis of AF.
The data of gene expression profiles have been increased rapidly in recent years and bioinformatics method, as a powerful tool, has been used to deeply explore the pathophysiological processes. Differentially expressed genes (DEGs) represented the genes whose expression level changed in pathologic condition compared to healthy condition. Alterations of gene expression were involved in the pathogenesis of AF, thus studying the DEGs between AF tissues and sinus rhythm (SR) tissues was necessary. In this study, we used the publicly available gene microarray data to perform de novo analysis and aimed to construct an updated prediction of new genes as biomarkers of AF to supplement currently available data.
Acquisition of data of gene expression profiles
The data of gene expression profiles were acquired by logging in the Gene Expression Omnibus (GEO) database affiliated to National Center for Biotechnology Information (NCBI) (https://www.ncbi.nlm.nih.gov) and editing key word “atrial fibrillation”. Three hundred eight results were shown including 42 series, 3 datasets, 2 platforms and 261 samples with detailed information. Among the 42 series, 29 series belonged to Homo sapiens, 9 series belonged to Mus musculus, 2 series belonged to Rattus norvegicus, 1 series belonged to Canis lupus familiaris and 1 series belonged to Sus scrofa. We divided the 29 series which belonged to Homo sapiens into 2 groups: expression profiling by high throughput sequencing and expression profiling by array and chose the latter group for further research. Microarray data which included RNA content of left atrial appendage (LAA) were selected. According to the criteria mentioned above, five microarray data (GSE79768, GSE31821, GSE115574, GSE14975 and GSE41177) were selected and downloaded. In GSE79768, there were 7 cases of samples from patients with AF and 6 cases of samples from patients without AF. In GSE31821, the number of cases from AF patients and SR patients was 2 and 2, respectively. In GSE115574, the number of patients with AF was 14 and the number of patients without AF was 15. GSE14975 included 5 atrial tissue samples from patients with AF and 5 from patients without AF. Sixteen samples from AF patients and 3 samples from SR patients were included in GSE41177. More detailed information is shown in Table 1.
Preprocessing of raw data, screening of differential genes and drawing of volcano plot and heat map
First, GEO Sample (GSM) which stands for LAA was selected. Then the downloaded platform files and series matrix files were transformed by Perl language software 5.28.1. The probe ID in the matrix files was converted into gene symbol by Perl language and saved as a TXT file. After that, five datasets were merged by Perl language and batch normalization was conducted using the “sva” package of the R software. We deleted duplicated genes and values lacking specific gene symbols. DEGs analysis was conducted using the “limma” package of R software. DEGs were selected by using cut-off values of P < 0.05 and |log2FC| > 0.5, where FC = fold change. Finally, volcano plot was drawn according to the data above and the “pheatmap” package of R software was applied to construct heat map.
DEGs functional enrichment analysis
On the basis of the DEGs from the merged five datasets, gene ontology (GO) enrichment analysis and pathway enrichment analysis were performed using the “clusterProfiler” package of R software and reactome [10, 11] (https://reactome.org/), respectively. GO analysis is a widely used bioinformatics tool to investigate the annotation of genes and proteins. It can be used to integrate annotation data and provide tools access to all the data provided by the project. Reactome is an open-source, open access, manually curated and peer-reviewed pathway database. The goal is to provide intuitive bioinformatics tools for the visualization, interpretation and analysis of pathway knowledge to support basic and clinical research, genome analysis, modeling, systems biology and education. A q value< 0.05 was identified as significant difference in GO analysis; false discovery rate (FDR) < 0.05 was identified as significant difference in pathway analysis.
Protein-protein interaction network analysis
We applied Search Tool for the Retrieval of Interacting Genes (STRING)  11.0 (https://string-db.org/) to construct a protein-protein interaction (PPI) network based on the DEGs detected above. Interactions of proteins in the STRING database were screened by using an interaction score criterion. PPI pairs with an interaction score > 0.4 were considered significant in this analysis.
Reverse transcription-quantitative PCR
From 2017.6 to 2018.2, LAA tissues were collected from 20 AF patients who underwent mitral valve repair and maze procedure and 20 SR patients with left atrial thrombus who underwent mitral valve repair, LAA resection and thrombectomy. The study was approved by the Ethics Committee of the First Affiliated Hospital of China Medical University (NO.2017-69-2) and was carried out in accordance with Declaration of Helsinki. Signed informed consents were obtained from all patients prior to tissue collection.
We selected three key genes for validation using reverse transcription-quantitative PCR (RT-qPCR). Total RNA was extracted with TRIzol Reagent (Invitrogen, USA) according to the manufacturer’s instructions. RNA was converted to cDNA using the PrimeScript™ RT reagent Kit (TaKaRa, Japan). RT-qPCR was performed using TB Green™ Premix Ex Taq™ II (TaKaRa, Japan) and the data were analyzed with LightCycler® 480 Instrument (Roche, Switzerland) and the relative changes were calculated using the 2-△△Ct method. SPP1 primers were: forward 5′- CTGTGTTG GTGGAGGATGTCTGC − 3′, and reverse 5′- GTCGGCGTTTGGCTGAGAAGG -3′. COL5A1 primers were: forward 5′- CGTATGATGACCTCACCTATGG − 3′, and reverse 5′- CGTAGTAGTTCTCGTCAAGGTT-3′. VCAN primers were: forward 5′- ACTGAAACTTCCTACGTATGCA − 3′, and reverse 5′- CTCACAAAGTG CACCAACATAA − 3′. β-actin primers were: forward 5′-CCTGGCACCCAGCA CAAT-3′, and reverse 5′- GGGCCGGACTCGTCATAC-3′.
All data analysis was performed with the SPSS system for statistics (IBM SPSS Statistics 22.0). Data were expressed as mean ± standard deviation. The Student’s t-test was used for comparison between two groups. Differences were considered to be significant at a level of P < 0.05.
Identification of DEGs between AF and SR LAA tissues
A total of 106 DEGs were obtained in the merged dataset. Among these DEGs, 74 genes were up-regulated and 32 genes down-regulated. The top 10 DEGs are shown in Table 2. The detailed information about the 106 DEGs is shown in Additional file 1. Volcano plot is shown in Fig. 1 and heat map shown in Fig. 2.
GO analysis and reactome pathway enrichment
GO analysis of genes includes biological process, cell composition and molecular function. The results are shown in Fig. 3. GO biological process analysis found that DEGs were mainly enriched in extracellular matrix organization, protein activation cascade, synapse pruning, extracellular structure organization, complement activation and collagen fibril organization. In the cell composition part, the DEGs were involved in collagen-containing extracellular matrix, extracellular matrix, collagen trimer, endoplasmic reticulum lumen, basement membrane and extracellular matrix component. In the molecular function section, the genes participated mainly in extracellular matrix structure constituent, extracellular matrix structural constituent conferring tensile strength, extracellular matrix structural constituent conferring compression resistance and collagen binding. The results suggested that DEGs were mostly involved in extracellular matrix structure. We also performed our reactome pathway analysis, which showed DEGs were mainly enriched in Extracellular matrix organization and Collagen formation. All the pathways involved are shown in Table 3.
PPI network integration
We used the STRING database to investigate PPI network, including 74 up-regulated genes and 32 down-regulated genes. Four major proteins were discovered in the PPI network analysis, with secreted phosphoprotein 1(SPP1) connecting 12 nodes, with collagen type V alpha 1 chain (COL5A1) connecting 11 nodes, with versican (VCAN) connecting 11 nodes and TYRO protein tyrosine kinase binding protein (TYROBP) connecting 11 nodes. A PPI network of DEGs and the key proteins were performed as shown in Figs. 4 and 5, respectively. Disconnected nodes in the network were hidden.
The validation of key genes expression
We performed RT-qPCR in LAA tissues to further verify the expression of three key genes. As illustrated in Fig. 6, compared with SR group, the transcription levels of SPP1, COL5A1 and VCAN in AF group were significantly increased (P < 0.05), which is consistent with our bioinformatics analysis results.
AF is a major medical problem bringing huge burden to the individual and the society . But the treatment of AF is far from satisfactory, so research on the pathogenesis and biomarkers of AF is still necessary and exploring the molecular level dysfunction can provide effective treatment and more diagnostic biomarkers. For example, Wang et al.  identified some miRNAs and genes as the target molecules in AF development from two GEO datasets by DEGs screening, construction of regulatory network and GO enrichment analysis. Barth et al.  identified key genes in permanent AF from one GEO dataset by integrated method including DEGs screening, hierarchical clustering and GO analysis.
In the present study, we used publicly available microarray data and bioinformatic approaches to predict the potential key genes associated with the pathogenesis of AF. The gene chip data of 5 datasets from GEO were selected and integrated to find a breakthrough point. Using multiple samples and a large amount of microarray data made the experimental results more reliable, providing some valuable references for prevention and treatment of AF. All the samples included were LAA tissue, which avoided the diversity of tissue from different parts. A total of 106 DEGs were obtained. The vast majority of them were up-regulated and the top 10 were CHGB, IGFBP2, LBH, FHL2, ANGPTL2, SNAI2, C1orf105, COLQ, C1orf168, OTOGL.
After that, GO enrichment analysis, pathway analysis and PPI network analysis were applied to analyze the differences between AF and SR patients and investigate the molecular pathogenesis of AF. GO analysis showed that the DEGs mainly participated in extracellular matrix organization, protein activation cascade, extracellular structure organization and collagen fibril organization. These biological processes are important processes involved in the pathophysiological mechanism of AF. Recent research showed that there is a positive correlation between fibrosis and AF, and down-regulation of EZH2 could inhibit Ang-II-induced atrial enlargement and fibrosis and reduce AF vulnerability . Reactome pathway analysis revealed that the DEGs were associated with the extracellular matrix (ECM) organization and collagen formation. More and more studies showed structural remodeling induced by atrial fibrosis and ECM alteration was one of the causes of initiation and perpetuation of AF.
Then, we mapped a PPI network in which the top 3 most influential proteins were SPP1, COL5A1 and VCAN, and they were all enriched in extracellular matrix organization.
SPP1 is involved in the formation and degradation of matrix. Kramerova  showed that Spp1 acts upstream of TGFβ to promote fibrosis in x-linked muscular dystrophy (mdx) muscles. Postnatal pharmacological inhibition of Spp1 reduces fibrosis and improves muscle function in mdx mice. Morse  demonstrated that in idiopathic pulmonary fibrosis (IPF), macrophages highly expressing SPP1 contribute importantly to lung fibrosis. Ruberti  found that SPP1 plasma levels are significantly higher in Primary myelofibrosis compared with essential thrombocytemia and polycytemia vera patients.
COL5A1 is an alpha chain for one of the low abundance fibrillar collagens. Fibrillar collagen molecules are trimers that can be composed of one or more types of alpha chains. Mutations in this gene are associated with Ehlers-Danlos syndrome [20, 21]. Alternative splicing of this gene results in multiple transcript variants. In Sun’s study , the collagen-V-knockout stroma demonstrated severe dysfunctional regulation of fibrillogenesis, which indicate central regulatory roles for collagen V in fibril and matrix assembly during tissue development.
VCAN is a member of the aggrecan/versican proteoglycan family. The protein encoded is a large chondroitin sulfate proteoglycan and is a major component of the extracellular matrix. This protein is involved in cell adhesion, proliferation, migration  and angiogenesis  and plays a central role in tissue morphogenesis and maintenance. Mutations in this gene are the cause of Wagner syndrome [25,26,27]. Burns  found loss of exon 7 in the Vcan gene in mice resulted in ventricular septal defects and altered pulmonary and aortic outflow tracts.
To further validate the expression of the key genes, we tested the LAA tissue by qPCR. We found that the relative expression levels of SPP1, COL5A1 and VCAN showed the same expression trend as in our bioinformatics analysis. There have been some researches on their effects on extracellular matrix. Extracellular matrix alteration, especially fibrosis, is one of the main causes for the development of AF . The DEGs detected in this article were mainly related to the extracellular matrix, which indicated that it is valuable to further study the DEGs and the key genes may become new diagnostic and therapeutic targets for AF.
There are some advantages in our research. First, to the best of our knowledge, our research is the first study to integrate all published microarray datasets related to human AF and LAA mRNA content in the GEO database for bioinformatics analysis. Second, the present study truly reflects the genetic-level changes in the pathophysiological process of AF and can provide possible target molecules for further research. Third, different bioinformatics methods are applied in this research, which can provide more ideas for data analysis. Finally, the chips in this article are published in recent years and all the samples are LAA tissue, which makes the analysis results more accurate and believable. However, the cut-off value we selected in the present study are relatively loose. Further experimental studies are needed to confirm the identified genes and pathways.
In the merged dataset, we found that SPP1, COL5A1 and VCAN, which are related to extracellular matrix, function together in AF. The present study provided useful information for further exploration of the pathogenesis of AF. Strengthening the research on DEGs of AF will have an important impact on the diagnosis and treatment of the disease.
Availability of data and materials
The datasets generated and/or analyzed during the current study are available in the Gene Expression Omnibus repository and the web links are shown below: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE79768 to GSE79768, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE31821 to GSE31821, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE115574 to GSE115574, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE14975 to GSE14975 and https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE41177 to GSE41177.
Differentially expressed genes
Left atrial appendage
Reverse transcription-quantitative polymerase chain reaction
Gene expression omnibus
National Center for biotechnology information
Secreted phosphoprotein 1
Collagen type V alpha 1 chain
TYRO protein tyrosine kinase binding protein
Friberg L, Bergfeldt L. Atrial fibrillation prevalence revisited. J Intern Med. 2013;274(5):461–8.
Piccini JP, Hammill BG, Sinner MF, Jensen PN, Hernandez AF, Heckbert SR, Benjamin EJ, Curtis LH. Incidence and prevalence of atrial fibrillation and associated mortality among Medicare beneficiaries, 1993-2007. Circ Cardiovasc Qual. 2012;5(1):85–U162.
Rahman F, Kwan GF, Benjamin EJ. Global epidemiology of atrial fibrillation. Nat Rev Cardiol. 2014;11(11):639–54.
Staerk L, Wang B, Preis SR, Larson MG, Lubitz SA, Ellinor PT, McManus DD, Ko D, Weng LC, Lunetta KL, et al. Lifetime risk of atrial fibrillation according to optimal, borderline, or elevated levels of risk factors: cohort study based on longitudinal data from the Framingham heart study. BMJ. 2018;361:k1453.
Freedman B, Potpara TS, Lip GYH. Stroke prevention in atrial fibrillation. Lancet. 2016;388(10046):806–17.
Stewart S, Hart CL, Hole DJ, McMurray JJ. A population-based study of the long-term risks associated with atrial fibrillation: 20-year follow-up of the Renfrew/Paisley study. Am J Med. 2002;113(5):359–64.
Thrall G, Lane D, Carroll D, Lip GY. Quality of life in patients with atrial fibrillation: a systematic review. Am J Med. 2006;119(5):448 e441–19.
Benjamin EJ, Wolf PA, D'Agostino RB, Silbershatz H, Kannel WB, Levy D. Impact of atrial fibrillation on the risk of death: the Framingham heart study. Circulation. 1998;98(10):946–52.
Friberg L, Hammar N, Pettersson H, Rosenqvist M. Increased mortality in paroxysmal atrial fibrillation: report from the Stockholm cohort-study of atrial fibrillation (SCAF). Eur Heart J. 2007;28(19):2346–53.
Fabregat A, Jupe S, Matthews L, Sidiropoulos K, Gillespie M, Garapati P, Haw R, Jassal B, Korninger F, May B, et al. The Reactome pathway knowledgebase. Nucleic Acids Res. 2018;46(D1):D649–55.
Fabregat A, Sidiropoulos K, Viteri G, Forner O, Marin-Garcia P, Arnau V, D'Eustachio P, Stein L, Hermjakob H. Reactome pathway analysis: a high-performance in-memory approach. BMC Bioinformatics. 2017;18(1):142.
Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, Simonovic M, Doncheva NT, Morris JH, Bork P, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47(D1):D607–13.
Chugh SS, Roth GA, Gillum RF, Mensah GA. Global burden of atrial fibrillation in developed and developing nations. Glob Heart. 2014;9(1):113–9.
Wang T, Wang B. Identification of microRNA-mRNA interactions in atrial fibrillation using microarray expression profiles and bioinformatics analysis. Mol Med Rep. 2016;13(6):4535–40.
Barth AS, Merk S, Arnoldi E, Zwermann L, Kloos P, Gebauer M, Steinmeyer K, Bleich M, Kaab S, Hinterseer M, et al. Reprogramming of the human atrial transcriptome in permanent atrial fibrillation: expression of a ventricular-like genomic signature. Circ Res. 2005;96(9):1022–9.
Song S, Zhang R, Mo B, Chen L, Liu L, Yu Y, Cao W, Fang G, Wan Y, Gu Y, et al. EZH2 as a novel therapeutic target for atrial fibrosis and atrial fibrillation. J Mol Cell Cardiol. 2019;135:119–33.
Kramerova I, Kumagai-Cresse C, Ermolova N, Mokhonova E, Marinov M, Capote J, Becerra D, Quattrocelli M, Crosbie RH, Welch E, et al. Spp1 (osteopontin) promotes TGFbeta processing in fibroblasts of dystrophin deficient muscles through matrix metalloproteinases. Hum Mol Genet. 2019;28(20):3431–42.
Morse C, Tabib T, Sembrat J, Buschur KL, Bittar HT, Valenzi E, Jiang Y, Kass DJ, Gibson K, Chen W, et al. Proliferating SPP1/MERTK-expressing macrophages in idiopathic pulmonary fibrosis. Eur Respir J. 2019;54(2):1802441.
Ruberti S, Bianchi E, Guglielmelli P, Rontauroli S, Barbieri G, Tavernari L, Fanelli T, Norfo R, Pennucci V, Fattori GC, et al. Involvement of MAF/SPP1 axis in the development of bone marrow fibrosis in PMF patients. Leukemia. 2018;32(2):438–49.
Junkiert-Czarnecka A, Pilarska-Deltow M, Bak A, Heise M, Haus O. New variants in COL5A1 gene among polish patients with Ehlers-Danlos syndrome: analysis of nine cases. Postepy Dermatol Alergol. 2019;36(1):29–33.
Ritelli M, Dordoni C, Venturini M, Chiarelli N, Quinzani S, Traversa M, Zoppi N, Vascellaro A, Wischmeijer A, Manfredini E, et al. Clinical and molecular characterization of 40 patients with classic Ehlers-Danlos syndrome: identification of 18 COL5A1 and 2 COL5A2 novel mutations. Orphanet J Rare Dis. 2013;8(1):58.
Sun M, Chen S, Adams SM, Florer JB, Liu H, Kao WW, Wenstrup RJ, Birk DE. Collagen V is a dominant regulator of collagen fibrillogenesis: dysfunctional regulation of structure and function in a corneal-stroma-specific Col5a1-null mouse model. J Cell Sci. 2011;124(Pt 23):4096–105.
Zhao X, Sun Z, Li H, Jiang F, Zhou J, Zhang L. MiR-135a-5p modulates biological functions of thyroid carcinoma cells via targeting VCAN 3′-UTR. Cancer Biomark. 2017;20(2):207–16.
Sathyan S, Koshy LV, Balan S, Easwer HV, Premkumar S, Nair S, Bhattacharya RN, Alapatt JP, Banerjee M. Association of Versican (VCAN) gene polymorphisms rs251124 and rs2287926 (G428D), with intracranial aneurysm. Meta Gene. 2014;2:651–60.
Ankala A, Jain N, Hubbard B, Alexander JJ, Shankar SP. Is exon 8 the most critical or the only dispensable exon of the VCAN gene? Insights into VCAN variants and clinical spectrum of Wagner syndrome. Am J Med Genet A. 2018;176(8):1778–83.
Burin-des-Roziers C, Rothschild PR, Layet V, Chen JM, Ghiotti T, Leroux C, Cremers FP, Brezin AP, Valleix S. Deletions overlapping VCAN exon 8 are new molecular defects for Wagner disease. Hum Mutat. 2017;38(1):43–7.
Kloeckener-Gruissem B, Neidhardt J, Magyar I, Plauchu H, Zech JC, Morle L, Palmer-Smith SM, Macdonald MJ, Nas V, Fry AE, et al. Novel VCAN mutations and evidence for unbalanced alternative splicing in the pathogenesis of Wagner syndrome. Eur J Hum Genet. 2013;21(3):352–6.
Burns TA, Dours-Zimmermann MT, Zimmermann DR, Krug EL, Comte-Walters S, Reyes L, Davis MA, Schey KL, Schwacke JH, Kern CB, et al. Imbalanced expression of Vcan mRNA splice form proteins alters heart morphology and cellular protein profiles. PLoS One. 2014;9(2):e89133.
Dzeshka MS, Lip GY, Snezhitskiy V, Shantsila E. Cardiac fibrosis in patients with atrial fibrillation: mechanisms and clinical implications. J Am Coll Cardiol. 2015;66(8):943–59.
We thank Dr. Lei Yu for assistance with the tissue collecting and colleagues at the laboratory for their support.
This work was supported by the Department of Education of Liaoning Province, China [grant number XLYC1802066].
Ethics approval and consent to participate
The study was approved by the Ethics Committee of the First Affiliated Hospital of China Medical University (NO.2017-69-2) and was carried out in accordance with Declaration of Helsinki. Signed informed consents were obtained from all patients prior to tissue collection.
Consent for publication
The authors have no conflicts of interest to declare.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Liu, Y., Tang, R., Zhao, Y. et al. Identification of key genes in atrial fibrillation using bioinformatics analysis. BMC Cardiovasc Disord 20, 363 (2020). https://doi.org/10.1186/s12872-020-01653-4
- Atrial fibrillation