GAS6-AS1, a long noncoding RNA, functions as a key candidate gene in atrial fibrillation related stroke determined by ceRNA network analysis and WGCNA
BMC Medical Genomics volume 16, Article number: 51 (2023)
Stroke attributable to atrial fibrillation (AF related stroke, AFST) accounts for 13 ~ 26% of ischemic stroke. It has been found that AFST patients have a higher risk of disability and mortality than those without AF. Additionally, it’s still a great challenge to treat AFST patients because its exact mechanism at the molecular level remains unclear. Thus, it’s vital to investigate the mechanism of AFST and search for molecular targets of treatment. Long non-coding RNAs (lncRNAs) are related to the pathogenesis of various diseases. However, the role of lncRNAs in AFST remains unclear. In this study, AFST-related lncRNAs are explored using competing endogenous RNA (ceRNA) network analysis and weighted gene co-expression network analysis (WGCNA).
GSE66724 and GSE58294 datasets were downloaded from GEO database. After data preprocessing and probe reannotation, differentially expressed lncRNAs (DELs) and differentially expressed mRNAs (DEMs) between AFST and AF samples were explored. Then, functional enrichment analysis and protein-protein interaction (PPI) network analysis of the DEMs were performed. At the meantime, ceRNA network analysis and WGCNA were performed to identify hub lncRNAs. The hub lncRNAs identified both by ceRNA network analysis and WGCNA were further validated by Comparative Toxicogenomics Database (CTD).
In all, 19 DELs and 317 DEMs were identified between the AFST and AF samples. Functional enrichment analysis suggested that the DEMs associated with AFST were mainly enriched in the activation of the immune response. Two lncRNAs which overlapped between the three lncRNAs identified by the ceRNA network analysis and the 28 lncRNAs identified by the WGCNA were screened as hub lncRNAs for further validation. Finally, lncRNA GAS6-AS1 turned out to be associated with AFST by CTD validation.
These findings suggested that low expression of GAS6-AS1 might exert an essential role in AFST through downregulating its downstream target mRNAs GOLGA8A and BACH2, and GAS6-AS1 might be a potential target for AFST therapy.
Atrial fibrillation (AF), affecting 25% of adults worldwide, is the most common clinical tachyarrhythmia  and is independently associated with a two-fold risk of mortality [2, 3]. Stroke attributable to atrial fibrillation (AF related stroke, AFST) accounts for 13 ~ 26% of ischemic stroke , and this proportion increases with age . AFST is characterized by a high percentage of early recurrent ischemic stroke  and hemorrhagic transformation (HT) in the days immediately following the index stroke . AFST patients have a worse prognosis, including higher risk of disability and mortality, than those without AF . Nowadays, a growing number of studies focus on preventing and intervening stroke in AF patients, however, the molecular mechanism of AFST is still not clearly understood, making its treatment a big challenge. Therefore, investigating the mechanism of AFST, as well as searching for the molecular targets for treatment, are of great clinical importance.
Long non-coding RNAs (lncRNAs) are a new kind of non-coding RNAs that lack of functional protein-coding ability , and are found of pronounced lower amounts than protein-coding genes. The function of lncRNAs in human transcription and epigenetics has been widely demonstrated . Numerous research has shown that lncRNAs are related to various diseases, including cancer, heart failure, myocardial infarction and diabetes [11,12,13,14]. Despite these findings, the mechanism of lncRNAs in AFST remains unclear. According to the competing endogenous RNA (ceRNA) hypothesis, lncRNA can regulate messenger RNA (mRNA) expression as miRNA sponge . By constructing disease-associated lncRNA-miRNA-mRNA regulatory ceRNA network, it is possible to identify disease-associated hub lncRNAs.
The weighted gene co-expression network analysis (WGCNA) is a relatively recent method to investigate the complex association between genes and clinical characteristics . WGCNA can aggregate co-expressed genes into modules to identify disease-related hub genes. Co-expression modules associated with diseases can be constructed not only using mRNAs, but also miRNAs or lncRNAs [17, 18]. The method has been widely used to study plenty of diseases, including cancer , severe asthma , and proved to be an effective method to identify potential therapeutic molecular targets.
In this study, we aimed to identify potential hub lncRNAs associated with AFST using ceRNA network analysis and WGCNA.
Materials and methods
In the current study, we integrated two datasets from the Gene Expression Omnibus (GEO) database. To uncover lncRNAs involved in AFST pathogenesis, it was imperative to combine diverse methods or biology algorithms, thus we conducted a series of analyses including differential expression analysis, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses, protein–protein interaction (PPI) network of the differentially expressed mRNA (DEMs) and cluster analysis, WGCNA, ceRNA network analysis, Comparative Toxicogenomics Database (CTD) validation, prognostic analysis based on Receiver operating characteristics (ROC). The workflow was illustrated in Fig. 1.
GEO is a public genomic data repository containing array-based data . Following screening, two datasets of GSE66724  and GSE58294 , both of which were annotated using GPL570 [HG-U133_Plus_2] Affymetrix Human Genome U133 Plus 2.0 Array, were downloaded from GEO database. Since the two datasets shared the same platform, these two candidates were selected for the integrated analysis. In all, 16 blood samples were collected from 8 patients with AF but no stroke (AF group), and 8 patients with both AF and stroke (AFST group) in GSE66724. Blood samples of GSE58294 were collected from patients with AF and stroke (AFST group, n = 69) and patients with AF but no stroke (AF group, n = 23). In GSE58294, all blood samples were obtained during the acute phase of the stroke.
Data preprocessing and probe reannotation
R packages of “affy” and “limma” were applied to assess GSE66724 and GSE58294 RAW data. The data were preprocessed by Robust Multi-array Average (RMA) procedure, and then the data of these two datasets were integrated for the subsequent analysis. Then, we marked different datasets as different batches, and used the “Combat” function in the “sva” package of R software to adjust the batch effect between the two datasets, then the principal component analysis (PCA) cluster plot was drawn to illustrate the samples before and after the batch effect removal. Reannotation of Affymetrix microarray probes to lncRNAs was performed according to the literature . Only lncRNAs with mean expression values > 0.5 in each sample were selected, finally, 1347 lncRNAs were obtained. Before proceeding to the next step, the expression value was normalized using “normalizeBetweenArrays” function in the “limma” package. The repeatability of the data was also validated by the PCA . The PCA and PCA cluster plots were carried out by the “FactoMineR” and “Factoextra” packages.
Differentially expressed lncRNA (DELs) and differentially expressed mRNAs (DEMs) analyses
The “limma” package was used to explore DELs and DEMs between AFST and AF samples using the empirical bayes method . Benjamin multiple test calibration was used to calculate the false discovery rate (FDR). Finally, the FDR < 0.05 and Fold change (FC) > 1.5 was taken as the threshold to select DELs and DEMs. Thereafter, a volcano plot of the DELs and DEMs was plotted using the “ggplot2” package. A hierarchical cluster heatmap was plotted to represent DEL and DEM expression intensity using the “pheatmap” package.
Functional enrichment analysis of the DEMs
With GO enrichment analysis, genes could be annotated using dynamic, controlled terms, which were distributed into biological processes (BP), cellular components (CC), and molecular functions (MF). In KEGG analysis, genomic information was linked to higher-order functional information and specific pathways. We used the “clusterProfiler” package to analyze the enrichment of GO terms and KEGG pathways in DEMs. Adjusted p value < 0.05 as well as q value < 0.05 were applied as the detection threshold, and the enrichment results were displayed using a dot graph and GOcircle plot.
At the same time, GO enrichment analysis and KEGG enrichment analysis were also performed based on Metascape . The p-value < 0.01 was applied as the detection threshold. Then, a network representing the enriched GO terms and KEGG pathways was constructed. The network was visualized using Cytoscape software (V3.6.0) and nodes that represent the enriched terms and pathways were colored according to cluster ID and p-value . Based on the DEMs identified in our study, we performed the gene-pathway crosstalk analysis to investigate the interactions among significantly enriched genes and pathways using the ClueGO and Cluepedia plug-in of Cytoscape, and the enriched genes and pathways were mapped into a crosstalk network.
Identification of protein-protein interaction (PPI) networks of DEMs
PPI network analysis of DEMs was performed using Metascape. A network was constructed when proteins interacted with each other. Subsequently, Cytoscape software (V3.6.0) was applied to visualize and analyze the network, and the topological features including the degree, closeness, betweenness of the nodes in the PPI network were calculated using the CentiScaPe plug-in of Cytoscape. In order to search clusters, the Molecular Complex Detection (MCODE) plug-in of Cytoscape was used.
CeRNA network construction
The MiRcode database (http://www.mircode.org/), which included presumed interactions between lncRNAs and miRNAs, was used to predict DELs’ relevant target miRNAs . Then, according to the miRTarBase (http://miRTarBase.cuhk.edu.cn/), miRDB (http://mirdb.org), and TargetScan (http://www.targetscan.org) databases [29,30,31], the aforementioned miRNAs’ relevant target mRNAs were predicted. Only the mRNAs that were identified in all three databases were screened as target mRNAs. In summary, the final ceRNA network contained the DELs, the predicted miRNAs, and the intersection of the target mRNAs and DEMs.
Identification of Hub lncRNAs through WGCNA
To explore the association between genes and clinical traits, the lncRNA expression matrix was extracted from the merged dataset. All 1347 lncRNAs were chosen to construct the co-expression modules following the instruction of “WGCNA” package . First, we used the “picksoftthreshold” function in the “WGCNA” package to calculate the soft threshold power β for each module. Following the β being settled down, the adjacency matrix was constructed and transformed into a topological overlap matrix (TOM). Then, hierarchical clustering and dynamic tree cut were performed with a merging cut-off value of 0.25 to determine co-expression modules.
The module eigengene (ME) was a weighted average gene expression value and indicated the overall expression level of the module. Then, pearson's correlation analysis was performed on MEs and clinical traits, allowing the identification of the modules which were significantly associated with the external traits. To further verify the module-trait correlation, we also calculated the module significance (MS, defined as the average absolute GS of all genes in the module). In general, modules with high MS values were considered as key modules. For each module, gene significance (GS) represented the association between genes and clinical traits, and module membership (MM) represented the association between genes and MEs. In the key modules, lncRNAs with |GS|> 0.6 and |MM|> 0.5 were identified as AFST related hub lncRNAs.
Using a Venn diagram, the intersection between the hub lncRNAs identified by WGCNA and ceRNA network analysis was determined. Next, using Cytoscape software (V3.6.0), we built a sub-ceRNA regulatory network including the overlapped hub lncRNAs, its target miRNAs, and the downstream mRNAs.
Further validation of the lncRNAs and mRNAs in the sub-ceRNA network
The CTD (http://ctd.mdibl.org) provided information about the associations between gene products, phenotypes, and diseases . Using the CTD, we were able to identify the potential relationship between lncRNAs and mRNAs in our sub-ceRNA network and the diseases of AF and stroke, with the inference score indicating the strength of association. The genes with high inference scores were identified as having potential clinical implications. Then the expression profiles of the genes were shown and ROC curves were generated to evaluate their diagnostic accuracy, and sensitivity and specificity were assessed using the area under the curve (AUC).
Identification of DELs and DEMs in AFST
After data preprocessing, merging, and reannotation of GSE66724 and GSE58294 (Additional files 1 and 2), 54,674 probes corresponding to 18,084 genes, which contained 1347 lncRNAs and 16,737 protein-coding genes, were obtained. According to PCA, significant differences between AF and AFST samples were found (Fig. 2A). Using a threshold of FC > 1.5 and FDR < 0.05, a total of 19 DELs and 317 DEMs were identified between AFST samples and AF samples (Additional files 3 and 4). In the AFST samples, 6 DELs were upregulated, 13 were downregulated; while out of 317 DEMs, 168 were upregulated, 149 were downregulated. A volcano plot and a heatmap of the DELs or DEMs were shown in Fig. 2. In the heatmap, the top 100 DELs or DEMs according to the value of |logFC| were shown and the AFST samples and AF samples were clearly distinguishable from the heatmap.
Functional enrichment analysis of the DEMs
The enrichment analyses of GO and KEGG pathways with a cut-off value of adjusted p-value < 0.05 as well as q-value < 0.05 were presented in Additional files 5 and 6, where the top 20 GO terms and KEGG pathways were shown according to the adjusted p-value. As shown in Fig. 3A and Additional file 7, activation of immune response, immune response-regulating cell surface receptor signaling pathway, and antigen receptor-mediated signaling pathway were dominant enriched BP terms, meanwhile, the enriched MF term was immune receptor activity. The concentric circle diagram of the GO analysis was shown in Additional file 8. Moreover, the KEGG enrichment analysis showed that the complement and coagulation cascade was the most enriched pathway, followed by hematopoietic cell lineage, NF-kappa B signaling pathway, B cell receptor signaling pathway (Fig. 3B). The significantly enriched terms and pathways might contribute to a further understanding of the role played by DEMs in AFST.
Additionally, we used Metascape to analyze functional enrichment, and the enriched terms were integrated into the networks by cluster ID and p-value. Nodes with the same cluster ID were colored the same in Fig. 3C, and terms enriched with more genes tended to be more significant in Fig. 3D. At the same time, we performed the gene-pathway crosstalk analysis to investigate the interactions among significantly enriched genes and pathways using the ClueGO and Cluepedia plug-in of Cytoscape, a gene-pathway network was constructed to visualize the associations between the significantly enriched pathways and genes (Fig. 4).
As a result of the enrichment analysis described above, the DEMs associated with AFST were mainly enriched in the activation of immune response and complement and coagulation cascades. The results showed that AFST might be closely associated with the process of immune response and complement and coagulation cascades.
PPI network and cluster analysis
In order to better understand the DEM interactions, we used Metascape to analyze PPI network. The PPI network was composed of 216 nodes and 339 edges (Additional file 9), and the topological features including the degree, closeness, and betweenness of the nodes in the PPI network were showed in Additional file 10. Then we used the MCODE plug-in of Cytoscape to search for clusters in the network. Finally, according to k-core = 2, four clusters were identified (Additional files 9 and 11).
Construction of the ceRNA network
First, the miRcode database was applied to predict miRNAs interacting with DELs. In all, 165 interactions between 4 DELs and 109 unique miRNAs were determined (Additional file 12). Following that, the target mRNAs of the 109 miRNAs were predicted using the miRTarBase, miRDB, and TargetScan databases. In total, 688 interactions between 109 miRNAs and 599 distinct mRNAs were identified (Additional file 13). Based on the overlapped mRNAs of the 599 mRNAs and 317 DEMs, a ceRNA network consisting of 3 lncRNAs, 7 miRNAs, and 11 mRNAs was constructed (Table 1, Additional file 14). All the three lncRNAs (LINC00323, LINC00342, GAS6-AS1) were downregulated in AFST patients.
Identification hub lncRNAs through WGCNA
In order to further verify the hub lncRNAs, we performed WGCNA in which all 1347 lncRNAs were included to construct the co-expression modules. The samples were analyzed using hierarchical clustering, and four obvious outliers (GSM1406037, GSM1406065, GSM1630733, GSM1630739) were removed from the cohort before WGCNA (Fig. 5A). It was shown in Fig. 5B that a threshold power of 3 was sufficient for WGCNA. As illustrated in Fig. 5C, the final 7 modules were identified based on a hierarchical clustering and dynamic tree cutting algorithm (cut-off value was 0.25). The largest module (blue) contained 906 lncRNAs while the smallest one (pink) contained 21 lncRNAs. By WGCNA, genes without a distinct module assignment were grouped in a gray module and were dismissed in the following analysis. Furthermore, interactions between the seven modules were analyzed. Together with the eigengene adjacency heatmap, the dendrogram of the modules demonstrated a high level of co-expression module independence (Fig. 5D).
Using correlation analysis, we investigated the relationship between modules and external traits. The green module had the most negative correlation with AFST (r = − 0.74), while the brown module had the most positive correlation with AFST (r = 0.73). (Fig. 6A). Moreover, across all modules, the green module had the highest MS values, followed by the red module and the brown module (Fig. 6B). Therefore, taking together the results of correlation analysis and MS, the red module, green module, and brown module were identified as the core modules for AFST. In addition, the genes in the 3 modules were analyzed using GS and MM. The genes in the upper right section of Fig. 6C–E, which had high values of GS and MM, were significantly associated with AFST and were the most important elements of the three modules at the same time. Consequently, a total of 28 lncRNAs (Table 2) in the upper right section of Fig. 6C–E were considered for further analysis.
The overlapped lncRNAs of the three lncRNAs in the ceRNA network and the 28 lncRNAs identified through WGCNA, GAS6-AS1 and LINC00342, were identified as hub lncRNAs (Fig. 6F). These two lncRNAs, together with their target miRNAs and mRNAs, were applied to construct a sub-ceRNA network (Fig. 7). According to ceRNA theory, as miRNA sponges, lncRNAs were supposed to regulate mRNAs positively. In our sub-ceRNA network, two downregulated lncRNAs (GAS6-AS1, LINC00342) and four downregulated mRNAs (BCL7A, BACH2, GOLGA8A, EBF1) aligned with the ceRNA theory, and were considered for further investigation.
Further validation of the lncRNAs and mRNAs in the sub-ceRNA network
Then, using the CTD, we predicted the potential role of the aforementioned six genes in AF and stroke. The inference score for the RNAs targeted AF and stroke was shown in Table 3. Finally, one lncRNA, GAS6-AS1, and three mRNAs including BCL7A, BACH2, GOLGA8A turned out to be associated with AFST based on ceRNA network analysis and WGCNA, as well as CTD validation. GAS6-AS1 might function, at least in part, as a ceRNA to regulate BCL7A, BACH2, and GOLGA8A in AFST.
The expression levels of the four hub genes were shown in Additional file 15, which showed that GAS6-AS1, BCL7A, BACH2, and GOLGA8A expression were significantly lower in the AFST samples compared with the AF samples. Subsequently, ROC curves were performed to assess the diagnostic value of the hub genes for AFST, and it was shown that the AUC for GAS6-AS1 was 0.828. Similar results for BCL7A, BACH2, and GOLGA8A were presented in Additional file 16.
In the current study, 31 blood samples from AF patients and 77 blood samples from AFST patients were enrolled from two datasets. For the first time, we found that lncRNA GAS6-AS1 might be associated with AFST. Both ceRNA network analysis and WGCNA were performed to confirm the role of GAS6-AS1 in AFST. The two different methods yielded identical results regarding the function of GAS6-AS1 in AFST, which was further confirmed by CTD. The reliable results indicated that lncRNA GAS6-AS1 might be a potential predictor of AFST or a potential therapeutic target in treating AFST.
Several studies had assessed the biomarkers in AFST previously. It was suggested by Allende et al.  that Hsp70 protected AFST patients by preventing thrombosis without increasing bleeding risk and it would be a new target to treat AFST patients. Using the datasets of GSE79768 and GSE58294, Zou et al.  found that the expression of ZNF566, PDZK1IP1, ZFHX3, and PITX2 genes were related to AFST and may be potential therapeutic targets for it. Based on the datasets of GSE66724 and GSE58294, Zhang et al.  found that ten genes including SMURF2, CDC42, UBE3A, RBBP6, CDC5L, NEDD4L, UBE2D2, UBE2B, UBE2I, and MAPK1 were overexpressed in AFST patients. According to Li et al. , the factor of inflammation was supposed to be considered when treating AFST patients, and certain genes, including MEF2A, CAND1, PELI1, and PDCD4 were identified and might contribute to the pathogenesis of AFST. The inconsistency of the hub genes in different studies might be attributed to the different samples included and different analysis protocols. It was intriguing that all the aforementioned studies focused on the differentially expressed mRNAs, to our knowledge, no previous study had investigated the role of lncRNA in AFST.
In 1988, for the first time, Schneider and his colleagues identified six members of the growth-arrest-specific (GAS) family of genes . Located on chromosome 13q34, the GAS6 gene has been shown to contribute to cell proliferation. An antisense RNA of GAS6, named GAS6-AS1, which is transcribed from chromosome 13q34 too, also plays an important role in the pathogenesis of many kinds of cancers. In different cancers, the role of GAS6-AS1 on patients' prognosis is extensively inconsistent. GAS6-AS1 may play a tumor suppressor role in lung cancer . Similarly, a higher level of GAS6-AS1 expression is associated with a better survival in Non-Small-Cell Lung Cancer (NSCLC) patients . Nevertheless, GAS6-AS1 promotes the migration and proliferation of gastric cancer cells by enhancing their entry into S-phase . By sponging miR-370-3p, GAS6-AS1 contributes to the development of acute myeloid leukemia . The opposite results that both the oncogenic  effect and anti-oncogenic  effect are obtained in papillary renal cell carcinoma.
The role of GAS6-AS1 in stroke has rarely been investigated. It’s suggested that GAS6-AS1 may be related to an increased risk of HT after intravenous thrombolysis in acute ischemic stroke patients . In the current study, the association between GAS6-AS1 and AFST is reported for the first time. GOLGA8A, one of the target mRNAs of GAS6-AS1 in our ceRNA network, has been shown to be related to intracerebral hemorrhage too . So, the GAS6-AS1/hsa-miR-363-3p/GOLGA8A axis in our ceRNA network seems to be related to intracerebral hemorrhage. Meanwhile, AFST is characterized by a high percentage of HT in the days immediately after the stroke . Therefore, it is plausible to postulate an association between the GAS6-AS1/hsa-miR-363-3p/GOLGA8A axis and HT after AFST, which warrants further investigation.
Increasing evidence suggests that ischemic stroke is associated with profound immune responses in the blood and the activation of multiple immune cell subsets. However, there is still a debate over whether these immune responses are beneficial or detrimental . Therefore, it is crucial to identify specific molecular targets to develop a new immunomodulatory treatment to prevent the detrimental effect of immune responses after stroke . Functional enrichment analyses in our study reveal that the DEMs related to AFST are primarily enriched in the biological processes of activation of the immune response and complement and coagulation cascades. The result proposes that AFST may be correlated with the process of immune response. Therefore, the hub genes identified in our study may be the molecular targets that we are looking for to develop new immunomodulatory therapies.
Among the three target mRNAs of GAS6-AS1 in our ceRNA network, BACH2 has higher inference scores for AF and Stroke, at the same time, Bach2 has been suggested as an influential immune-regulating transcription factor in T helper 2 (Th2), Follicular T helper (Tfh), regulatory T cell (Treg), B cells and plays a key role in Th2 immune response previously . BCL7A tends to be related to cancer , but not stroke. Taking into account the inference scores and the biological function of the target mRNAs, it is possible that GAS6-AS1 downregulation may function in AFST patients by regulating BACH2 as a ceRNA through the immune response.
Collectively, we thus propose an association between the ceRNA axis GAS6-AS1/hsa-miR-363-3p/GOLGA8A and HT after AFST, and predict that the GAS6-AS1/hsa-miR-507/BACH2 axis has a potential role in AFST through inflammatory and immune responses. They may be potential targets for AFST therapy. The detailed mechanisms may need further investigation.
There are still limitations in our current study. First, although different approaches have been used to demonstrate the role of lncRNA GAS6-AS1, further validation is needed to confirm it. Second, due to the lower expression levels of lncRNAs compared to mRNAs, WGCNA is performed only for lncRNAs, and as a result, lncRNA-mRNA interactions may be missing. Most importantly, the potential mechanisms of the association between GAS6-AS1 and AFST was speculated on the basis of previous studies and bioinformatics analysis. Further experiments (both in vivo and in vitro) are desperately needed to verify our findings. In addition, gene expression differs in different stroke phase . All blood samples in GSE58294 are taken during the acute phase of the stroke, and we cannot rule out that samples taken at a different stroke phase may have yielded different results.
In conclusion, we identified a hub lncRNA of GAS6-AS1 associated with AFST by ceRNA network analysis and WGCNA. It was subsequently validated by CTD that GAS6-AS1 played a pivotal role in AFST. These findings suggested that low expression of GAS6-AS1 might exert an essential role in AFST through downregulating GOLGA8A and BACH2, by affecting post-AFST hemorrhagic transformation and post-AFST immune response, and pointed out the direction for further research. Altogether, these analyses suggested that GAS6-AS1 might represent a potential target for AFST therapy.
Availability of data and materials
The datasets generated for this study can be found in the GEO database (GSE66724 and GSE58294; https://www.ncbi.nlm.nih.gov/geo/).
Atrial fibrillation related stroke
The area under the curve
- BACH2 :
BTB domain and CNC homolog 2
- BCL7A :
B-cell CLL/lymphoma 7 protein family, member A
Competing endogenous RNA
Chronic obstructive pulmonary disease
Comparative toxicogenomics database
Differentially expressed mRNAs
Differentially expressed lncRNAs
False discovery rate
- GAS6-AS1 :
Growth-arrest-specific 6-antisense RNA1
Gene expression omnibus
Golgin A8 family, member A
Kyoto encyclopedia of genes and genomes
Long non-coding RNA
Molecular complex detection
Non-small-cell lung cancer
Principal component analysis
Robust multi-array average
Receiver operating characteristic
Follicular T helper
T helper 2
Topological overlap matrix
Regulatory T cell
Weighted gene co-expression network analysis
Zulkifly H, Lip G, Lane DA. Epidemiology of atrial fibrillation. Int J Clin Pract. 2018;72(3):e13070.
Chugh SS, Havmoeller R, Narayanan K, Singh D, Rienstra M, Benjamin EJ, et al. Worldwide epidemiology of atrial fibrillation: a Global Burden of Disease 2010 study. Circulation. 2014;129(8):837–47.
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.
Seiffge DJ, Werring DJ, Paciaroni M, Dawson J, Warach S, Milling TJ, et al. Timing of anticoagulation after recent ischaemic stroke in patients with atrial fibrillation. Lancet Neurol. 2019;18(1):117–26.
Wolf PA, Abbott RD, Kannel WB. Atrial fibrillation as an independent risk factor for stroke: the Framingham study. Stroke. 1991;22(8):983–8.
Hart RG, Coull BM, Hart D. Early recurrent embolism associated with nonvalvular atrial fibrillation: a retrospective study. Stroke. 1983;14(5):688–93.
D’Amelio M, Terruso V, Famoso G, Di Benedetto N, Realmuto S, Valentino F, et al. Early and late mortality of spontaneous hemorrhagic transformation of ischemic stroke. J Stroke Cerebrovasc Dis. 2014;23(4):649–54.
Lamassa M, Di Carlo A, Pracucci G, Basile AM, Trefoloni G, Vanni P, et al. Characteristics, outcome, and care of stroke associated with atrial fibrillation in Europe: data from a multicenter multinational hospital-based registry (The European Community Stroke Project). Stroke. 2001;32(2):392–8.
Ponting CP, Oliver PL, Reik W. Evolution and functions of long noncoding RNAs. Cell. 2009;136(4):629–41.
Kopp F, Mendell JT. Functional classification and experimental dissection of long noncoding RNAs. Cell. 2018;172(3):393–407.
Miao H, Lu J, Guo Y, Qiu H, Zhang Y, Yao X, et al. LncRNA TP73-AS1 enhances the malignant properties of pancreatic ductal adenocarcinoma by increasing MMP14 expression through miRNA -200a sponging. J Cell Mol Med. 2021;25(7):3654–64.
Yan L, Zhang Y, Zhang W, Deng SQ, Ge ZR. lncRNA-NRF is a potential biomarker of heart failure after acute myocardial infarction. J Cardiovasc Transl Res. 2020;13(6):1008–15.
Yang F, Chen Y, Xue Z, Lv Y, Shen L, Li K, et al. High-throughput sequencing and exploration of the lncRNA-circRNA-miRNA-mRNA network in Type 2 diabetes mellitus. Biomed Res Int. 2020;2020:8162524.
Liu CY, Zhang YH, Li RB, Zhou LY, An T, Zhang RC, et al. LncRNA CAIF inhibits autophagy and attenuates myocardial infarction by blocking p53-mediated myocardin transcription. Nat Commun. 2018;9(1):29.
Salmena L, Poliseno L, Tay Y, Kats L, Pandolfi PP. A ceRNA hypothesis: the Rosetta Stone of a hidden RNA language. Cell. 2011;46(3):353–8.
Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinform. 2008;9:559.
Hu Z, Yang D, Tang Y, Zhang X, Wei Z, Fu H, et al. Five-long non-coding RNA risk score system for the effective prediction of gastric cancer patient survival. Oncol Lett. 2019;17(5):4474–86.
Ma X, Tao R, Li L, Chen H, Liu Z, Bai J, et al. Identification of a 5-microRNA signature and hub miRNA-mRNA interactions associated with pancreatic cancer. Oncol Rep. 2019;41(1):292–300.
Wan Q, Tang J, Han Y, Wang D. Co-expression modules construction by WGCNA and identify potential prognostic markers of uveal melanoma. Exp Eye Res. 2018;166:13–20.
Modena BD, Bleecker ER, Busse WW, Erzurum SC, Gaston BM, Jarjour NN, et al. Gene expression correlated with severe asthma characteristics reveals heterogeneous mechanisms of severe disease. Am J Respir Crit Care Med. 2017;195(11):1449–63.
Wang Z, Monteiro CD, Jagodnik KM, Fernandez NF, Gundersen GW, Rouillard AD, et al. Extraction and analysis of signatures from the gene expression omnibus by the crowd. Nat Commun. 2016;7:12846.
Allende M, Molina E, Guruceaga E, Tamayo I, González-Porras JR, Gonzalez-López TJ, et al. Hsp70 protects from stroke in atrial fibrillation patients by preventing thrombosis without increased bleeding risk. Cardiovasc Res. 2016;110(3):309–18.
Stamova B, Jickling GC, Ander BP, Zhan X, Liu D, Turner R, et al. Gene expression in peripheral immune cells following cardioembolic stroke is sexually dimorphic. PLoS ONE. 2014;9(7):e102550.
Peng F, Wang R, Zhang Y, Zhao Z, Zhou W, Chang Z, et al. Differential expression analysis at the individual level reveals a lncRNA prognostic signature for lung adenocarcinoma. Mol Cancer. 2017;16(1):98.
Ringnér M. What is principal component analysis. Nat Biotechnol. 2008;26(3):303–4.
Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47.
Zhou Y, Zhou B, Pache L, Chang M, Khodabakhshi AH, Tanaseichuk O, et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. 2019;10(1):1523.
Jeggari A, Marks DS, Larsson E. miRcode: a map of putative microRNA target sites in the long non-coding transcriptome. Bioinformatics. 2012;28(15):2062–3.
Wong N, Wang X. miRDB: an online resource for microRNA target prediction and functional annotations. Nucleic Acids Res. 2015;43:D146-52.
Chou CH, Shrestha S, Yang CD, Chang NW, Lin YL, Liao KW, et al. miRTarBase update 2018: a resource for experimentally validated microRNA-target interactions. Nucleic Acids Res. 2018;46(D1):D296–302.
Agarwal V, Bell GW, Nam JW, Bartel DP. Predicting effective microRNA target sites in mammalian mRNAs. Elife. 2015;4:e05005.
Davis AP, Grondin CJ, Johnson RJ, Sciaky D, Wiegers J, Wiegers TC, et al. Comparative toxicogenomics database (CTD): update 2021. Nucleic Acids Res. 2021;49(D1):D1138–43.
Zou R, Zhang D, Lv L, Shi W, Song Z, Yi B, et al. Bioinformatic gene analysis for potential biomarkers and therapeutic targets of atrial fibrillation-related stroke. J Transl Med. 2019;17(1):45.
Zhang YF, Meng LB, Hao ML, Yang JF, Zou T. Identification of co-expressed genes between atrial fibrillation and stroke. Front Neurol. 2020;11:184.
Li Y, Tan W, Ye F, Wen S, Hu R, Cai X, et al. Inflammation as a risk factor for stroke in atrial fibrillation: data from a microarray data analysis. J Int Med Res. 2020;48(5):300060520921671.
Schneider C, King RM, Philipson L. Genes specifically expressed at growth arrest of mammalian cells. Cell. 1988;54(6):787–93.
Wang Y, Ma M, Li C, Yang Y, Wang M. GAS6-AS1 overexpression increases GIMAP6 expression and inhibits lung adenocarcinoma progression by sponging miR-24-3p. Front Oncol. 2021;11:645771.
Han L, Kong R, Yin DD, Zhang EB, Xu TP, De W, et al. Low expression of long noncoding RNA GAS6-AS1 predicts a poor prognosis in patients with NSCLC. Med Oncol. 2013;30(4):694.
Zhang P, Dong Q, Zhu H, Li S, Shi L, Chen X. Long non-coding antisense RNA GAS6-AS1 supports gastric cancer progression via increasing GAS6 expression. Gene. 2019;696:1–9.
Li Y, Ma J, Li B, Zhu X, Wang J. Cirrhosis of Wilson’s disease: High and low cutoff using acoustic radiation force impulse (ARFI) -comparison and combination with serum fibrosis index. Clin Hemorheol Microcirc. 2021;79(4):575–85.
Lan H, Zeng J, Chen G, Huang H. Survival prediction of kidney renal papillary cell carcinoma by comprehensive LncRNA characterization. Oncotarget. 2017;8(67):110811–29.
Yang F, Song Y, Ge L, Zhao G, Liu C, Ma L. Long non-coding RNAs as prognostic biomarkers in papillary renal cell carcinoma. Oncol Lett. 2019;18(4):3691–7.
Guo ZN, Liu J, Chang J, Zhang P, Jin H, Sun X, et al. GAS6/Axl signaling modulates blood-brain barrier function following intravenous thrombolysis in acute ischemic stroke. Front Immunol. 2021;12:742359.
Merino-Zamorano C, Delgado P, Fernández de Retana S, Fernández-Cadenas I, Rodríguez-Luna D, Montaner J, et al. Identification of plasma biomarkers of human intracerebral hemorrhage subtypes through microarray technology. J Stroke Cerebrovasc Dis. 2016;25(3):665–71.
Gelderblom M, Leypoldt F, Steinbach K, Behrens D, Choe CU, Siler DA, et al. Temporal and spatial dynamics of cerebral immune cell accumulation in stroke. Stroke. 2009;40(5):1849–57.
Rayasam A, Hsu M, Kijak JA, Kissel L, Hernandez G, Sandor M, et al. Immune responses in stroke: how the immune system contributes to damage and healing after stroke and how this knowledge could be translated to better cures. Immunology. 2018;154(3):363–76.
Liu G, Liu F. Bach2: a key regulator in Th2-related immune cells and Th2 immune response. J Immunol Res. 2022;2022:2814510.
Liu J, Gao L, Ji B, Geng R, Chen J, Tao X, et al. BCL7A as a novel prognostic biomarker for glioma patients. J Transl Med. 2021;19(1):335.
Oh SH, Kim OJ, Shin DA, Song J, Yoo H, Kim YK, et al. Alteration of immunologic responses on peripheral blood in the acute phase of ischemic stroke: blood genomic profiling study. J Neuroimmunol. 2012;249(1–2):60–5.
This study was supported by the department of Science and Technology of Hebei Province (grant number: 22377749D) and Medical Research Program of Hebei Province (grant number: 20230580).
Ethics approval and consent to participate
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. FigS1. Data distribution. (A) Data distribution of GSE66724 before normalization (B) Data distribution of GSE58294 before normalization (C) Data distribution of the merged dataset after data normalization.
FigS2. PCA plot of the data before and after the batch effect removal. (A) PCA results before the batch effect removal. (B) PCA results after the batch effect removal.
. Differentially expressed lncRNAs.
. Differentially expressed mRNAs.
. GO enrichment analysis of the DEMs.
. KEGG enrichment analysis of the DEMs.
. FigS3. GO terms plot of the DEMs, Colors in different plots indicate the level of significance. (A) Biological processes (B) Cellular components (C) Molecular functions.
. FigS4. Circle plot of the GO enrichment analysis. The left outer semicircle represents the logFC value of the genes, and the right semicircle corresponds to GO terms enriched.
. FigS5. Clusters of the PPI network based on the Metascape and MCODE analysis. Four colors of red, bule, yellow and green indicate four clusters identified by MCODE analysis.
. Topological features of the nodes in the PPI network.
. Genes in different cluster identified by MCODE.
. lncRNA target miRNA prediction.
. miRNA target mRNA prediction.
. FigS6. CeRNA regulatory network. Red rhombuses represent lncRNAs, green triangles represent miRNAs and blue circles represent mRNAs, respectively.
. FigS7. Boxplot of the expression level for four hub genes.
. FigS8. The receiver operator characteristic curves of GAS6-AS1, GOLGA8A, BACH2 and BCL7A for AFST.
About this article
Cite this article
Li, Rb., Yang, Xh., Zhang, Jd. et al. GAS6-AS1, a long noncoding RNA, functions as a key candidate gene in atrial fibrillation related stroke determined by ceRNA network analysis and WGCNA. BMC Med Genomics 16, 51 (2023). https://doi.org/10.1186/s12920-023-01478-y
- Atrial fibrillation
- Competing endogenous RNA
- Weighted gene co‑expression network analysis
- Long non-coding RNA