- Open Access
FAM201A, a long noncoding RNA potentially associated with atrial fibrillation identified by ceRNA network analyses and WGCNA
BMC Medical Genomics volume 15, Article number: 80 (2022)
Being the most common arrhythmia in clinic, atrial fibrillation (AF) causes various comorbidities to patients such as heart failure and stroke. LncRNAs were reported involved in pathogenesis of AF, yet, little is known about AF-associated lncRNAs. The present study aims to explore lncRNAs associated with AF susceptibility based on competing endogenous RNA (ceRNA) network analysis and weighted gene co-expression network analysis (WGCNA).
GSE41177 and GSE79768 datasets were obtained from the Gene Expression Omnibus (GEO) database. Competing endogenous RNA (ceRNA) network analysis was performed using GSE41177. Differentially expressed lncRNAs (DElncRNAs), mRNAs (DEmRNAs) between AF patients and patients with sinus rhythm (SR) were identified from GSE41177 using R software. Then, the ceRNA network was constructed based on DElncRNAs, the predicted target miRNAs and DEmRNAs. Weighted gene co-expression network analysis (WGCNA) was performed using GSE79768 to validate the AF-related lncRNAs identified from GSE41177. LncRNA modules and crucial lncRNAs relevant to AF and were identified.
In summary, 18 DElncRNAs and 350 DEmRNAs were found between AF patients and SR patients. A total of 5 lncRNAs, 10 miRNAs, and 21 mRNAs were contained in the final ceRNA network. Taking into consideration both the ceRNA theory and inference scores from the comparative toxicogenomics database (CTD) database, the ceRNA axis FAM201A-miR-33a-3p-RAC3 was identified as mostly relevant to AF susceptibility. FAM201A (Gene significance, GS = − 0.62; Module membership, MM = 0.75) was also proved in the blue module, which was identified most highly relevant with AF by WGCNA.
These results demonstrated that decreased expression of FAM201A might be associated with susceptibility of AF. Working as the ceRNA to regulate RAC3 might be one function of FAM201A in AF susceptibility, which requires further exploration in future research.
Atrial fibrillation (AF), the most common cardiac arrhythmia, affects approximately 34 million people worldwide, and the number increases with aging [1, 2]. Being a major contributor to stroke, heart failure, sudden death and myocardial infarction, AF poses a significant burden to patients and society . The clinical risk factors of AF include gender, alcohol consumption, smoking, obesity and some clinical comorbidities . Besides, a number of protein-coding genes have been identified significantly associated with risk of AF. Though, only a small proportion of heritability for AF has been uncovered [4, 5]. Until now, the pathogenesis of AF still remains poorly understood, limiting the discovery of novel therapeutic targets for AF.
Long noncoding RNAs (lncRNAs), with more than 200 nucleotides in length, are non-protein-coding RNAs, taking an important part in transcriptional and epigenetic gene regulation . Recently, a large number of studies have demonstrated that lncRNAs are involved in pathogenesis of various diseases, such as cancers, diabetes, heart failure and myocardial infarction [7,8,9,10]. More importantly, some lncRNAs were proved to be new biomarkers or therapeutic targets for AF, such as lncRNA HNRNPU-AS1, lncRNA PVT1, lncRNA GAS5 [11,12,13]. Still, little is known about lncRNAs associated with AF susceptibility, which deserves further discoveries.
One pivotal regulatory function of lncRNAs is working as competing endogenous RNA (ceRNA) to regulate mRNA transcription. It was stated in this ceRNA hypothesis that lncRNA could compete for miRNA via shared microRNA response elements, and subsequently sponge miRNA to indirectly regulate mRNA expression . In the pathogenesis of AF, for example, lncRNA PVT1, which was increased in atrial samples of AF patients, acted as a sponge for miR-128-3p and up-regulated Sp1 expression to facilitate atrial fibrosis . LncRNA TCONS-00106987 acted as a sponge for miR-26 and up-regulated KCNJ2 expression to trigger atrial electrical remodeling . Thus, analyses of ceRNA networks turned out to be an efficient method to discover AF-related lncRNAs.
Besides, another biology algorithm, weighted gene co-expression network analysis (WGCNA) , can be used to analyze crucial genes or novel biomarkers of various diseases including AF . Briefly, gene modules were constructed and associated with clinical traits to identify disease related modules and the key genes in each module. It is worth noting that not only protein-coding RNAs, but also non-protein-coding RNAs could be used to construct the co-expression networks associated with diseases [18, 19].
For short of AF-related lncRNA microarray datasets, combination of diverse methods or biology algorithms are deadly required to deeply discover lncRNAs involved in pathogenesis of AF. Therefore, in the present study, we aim to identify lncRNAs associated with AF susceptibility based on ceRNA network analyses, as well as WGCNA.
Materials and methods
Microarray data collection
The lncRNA and mRNA expression profiles of two datasets were retrieved from the Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo/). The dataset GSE41177 was used to identify lncRNAs associated with AF susceptibility based on ceRNA network construction, and GSE79768 was used to validate the potential role of lncRNAs for AF based on WGCNA. GSE41177 (platform: GPL570) contains left atrial appendages and paired pulmonary vein and the surrounding left atrial junctions from 16 AF patients and 3 patients with sinus rhythm (SR). GSE79768 (platform: GPL570) contains left atrial appendages and right atrial appendages from 7 AF patients and 6 SR patients. Data of GSE41177 were obtained from AF and SR patients undergoing valvular heart disease surgery. LA-PV junction specimens were taken from the area between the right-superior and right-inferior PVs adjacent to the atriotomy site, as described elsewhere . Data of GSE79768 were obtained from patients receiving surgery for mitral valve or coronary artery disease, as described elsewhere . In the present study, to achieve the sample consistency of the two datasets, only data of left atrial appendages in GSE41177 and GSE79768 were used for analyses. Thus, a total of 23 left arial appendages of AF patients and 9 left arial appendages of SR were included in the study. Additionally, there is evidence that gene expression files of left atrial and right atrial are not usually identical and left atrial is proved more pivotal in AF initiation and maintenance [22, 23].
Differential lncRNA and mRNA expression analyses
The lncRNA and mRNA expression profiles were retrieved from dataset GSE41177 for differential expression analyses and ceRNA network construction. Identification of differentially expressed lncRNAs (DElncRNAs), mRNAs (DEmRNAs) between AF patients and SR patients were performed using the “limma” package of R software (Version 3.6.3) . The adjusted P value < 0.05 and |log2-fold change (FC)|> 1 was taken as the threshold to select DElncRNAs and DEmRNAs. Hierarchical cluster heatmaps were generated to represent expression intensity and direction of DElncRNAs and DEmRNAs using the “pheatmap” package of R software based on Euclidean distance.
Construction of the ceRNA network
DElncRNAs associated with target miRNAs were predicted using the miRcode database (http://www.mircode.org/), which contains the putative interactions between lncRNAs and miRNAs . Then, these miRNAs associated target mRNAs were predicted based on the miRDB (http://mirdb.org), miRTarBase (http://miRTarBase.cuhk.edu.cn/), and TargetScan (http://www.targetscan.org) databases [26,27,28]. Only those identified in all three databases were selected as target mRNAs. The overlap of target mRNAs and DEmRNAs, together with DElnRNAs and predicted miRNAs, were obtained to construct the final ceRNA network using cytoscape software (Version 3.7.2) .
Gene ontology enrichment analyses for mRNAs in the ceRNA network
To reveal the potential biological functions of mRNAs in the ceRNA network, Gene Ontology (GO) enrichment analyses was performed using the “clusterProfiler” package of R software . The bar graph was generated to display the enrichment results. The threshold was set as the adjusted P value < 0.05.
Potential AF susceptibility lncRNAs and mRNAs prediction
The Comparative Toxicogenomics Database (CTD) (http://ctd.mdibl.org) database contains data of associations between chemicals, gene products, phenotypes, diseases, and environmental exposures . The CTD database was used to predict the potential AF-related lncRNAs and mRNAs in the ceRNA network, with the inference score reflecting the association between AF and lncRNAs or mRNAs.
Validation for potential role of lncRNAs through WGCNA
The lncRNA expression profiles were retrieved from dataset GSE79768 for WGCNA. All 1210 lncRNAs were selected to construct a co-expression network using the “WGCNA” package in R software. A proper soft threshold of 9 was chosen to satisfy the degree of independence of 0.85 with the minimum value. An adjacency matrix was constructed and converted into a topological overlap matrix to reflect the correlation strength in the co-expression network. Hierarchical clustering was performed using DynamicTreeCut algorithm to construct the network modules with the minimum module size of 30 and height cutoff of 0.25. The module eigengene (ME) is the first principal component of the module, representing the overall expression profiles of each module. To identify modules relevant to AF susceptibility, the correlation between MEs and atrial rhythm phenotypes were analyzed with Pearson’s correlation and visualized by the heatmap. A threshold of P < 0.05 was used to screen the modules significantly associated with AF susceptibility. Gene significance (GS) is defined as the correlation between each gene and the clinical trait. Module membership (MM) is defined as the correlation between each gene and each module eigengene. The lncRNAs satisfying the criteria of |GS|> 0.6 and |MM|> 0.5 were identified as crucial lncRNAs relevant to AF.
Differential expression analyses of lncRNAs and mRNAs in atrial fibrillation
Microarray data of GSE41177 were downloaded and used for analyses. A total of 18 DElncRNAs (Additional file 1) and 350 DEmRNAs (Additional file 2) were identified between AF patients and SR patients, using the threshold of an adjusted P value < 0.05 and |log2 FC|> 1 as cut-off. The volcano plots showed all RNA expression levels in AF samples compared to SR samples (Additional file 3). The heatmap showed the distinguished expression levels of lncRNAs (Fig. 1) and mRNAs (Additional file 4) between AF and SR groups.
Construction of the ceRNA network
We predict the miRNAs that interact with DElncRNAs based on miRcode database. A total of 222 interactions between 6 DElncRNAs and 157miRNAs were identified. Then the target mRNAs of abovementioned 157 miRNAs were predicted based on the miRDB, miRTarBase, and TargetScan databases. Overall, 992 mRNAs identified in all three databases were selected as target mRNAs. The overlap of these 992 mRNAs and 350 DEmRNAs were used to construct the ceRNA network, which finally contained a total of 5 lncRNAs, 10 miRNAs, and 21 mRNAs (Fig. 2, Tables 1 and 2, Additional file 5). Except for one up-regulated lncRNA CTD-3080P12, the other four lncRNAs (FAM201A, LINC00326, LINC00029, LINC00355) were all down-regulated in AF patients compared with SR patients.
Gene ontology enrichment analyses for mRNAs in the ceRNA network
The results of GO enrichment analyses with the screening criteria of adjusted P value < 0.05 were presented in Additional file 6 and Additional file 7. The mRNAs in the ceRNA network were mainly enriched in transmembrane receptor protein kinase activity, protein phosphatase binding, phosphatase binding, transmembrane receptor protein tyrosine kinase activity. Neither specific function nor specific pathway was identified in the enrichment analyses. The reason for that was speculated to be the small number of mRNAs enrolled in the enrichment analyses.
Potential AF associated lncRNAs and mRNAs prediction in the ceRNA netwrok
LncRNAs are supposed to positively regulate mRNAs by working as miRNA sponges based on ceRNA theory. In our ceRNA network, expression direction of seven down-regulated lncRNA and mRNA pairs were in accord with the ceRNA theory, including two DElncRNAs (FAM201A, LINC00355) and seven DEmRNAs (RAC3, MATR3, NR3C2, KAT7, CHL1, SUCO, STAM). Then, the CTD database was used to predict the potential role of above DElncRNAs and DEmRNAs in AF. Inference scores were used to reflect the associations between AF and above RNAs. Finally, one lncRNA and seven mRNAs including FAM201A (Inference score: 3.08), NR3C2 (Inference score: 24.69), RAC3 (Inference score: 22.57), KAT7 (Inference score: 16.61), STAM (Inference score: 9.62), MATR3(Inference score: 5.27), CHL1 (Inference score: 3.77), and SUCO (Inference score: 1.79) turned out to have potential associations with AF. Regarding mRNAs, NR3C2 and RAC3 have relatively higher inference scores for AF. Taking inference scores of both lncRNAs and mRNAs into consideration, the ceRNA axis FAM201A-miR-33a-3p-RAC3 was predicted to have potential role for AF susceptibility.
Validation for potential role of lncRNAs through WGCNA
WGCNA was performed to further validate the potential role of lncRNAs in AF. All 1210 lncRNA expression matrices from GSE7976 were screened for construction of the co-expression network. The sample clustering tree and trait heatmap was illustrated in Fig. 3a. The soft‐threshold of 9 was set to construct a scale-free network, with the scale‐free topology fit index > 0.85 (Fig. 3b). The final eight modules identified based on average hierarchical clustering and dynamic tree cutting were shown in Fig. 3c.
The module‐trait association analysis was performed by calculating the correlation between MEs and the AF phenotype. Blue module (r = − 0.82; P = 5 × 10–4) and black module (r = − 0.61; P = 0.03) were highly relevant with AF (Fig. 4a). The significant correlation was observed between MM and GS for AF in the blue module, as shown in Fig. 4b. A total of 95 lncRNAs in blue module were identified as crucial lncRNAs associated with AF, according to the criteria of |GS|> 0.6 and |MM|> 0.5 (Additional file 8). We mainly focused on the MM and GS for AF of lncRNAs, namely, FAM201A, LINC00355, LINC00326, LINC00029, LINC00355, CTD-3080P12, which were identified in the ceRNA network. Finally, FAM201A (GS: − 0.62, P.GS: 0.02; MM: 0.75, P.MM: 3.35 × 10−3) and LINC00326 (GS: − 0.60, P.GS: 0.03; MM: 0.78, P.MM: 1.58 × 10−3) were in the blue module and proved to be highly related to AF (Table 3).
Taken together, these results demonstrated that FAM201A might have great potential for susceptibility of AF based on the ceRNA network, CTD database and WGCNA. FAM201A may function, at least partly, as ceRNA to regulate RAC3 in AF susceptibility.
In the present study, 23 left arial appendages of AF patients and 9 left atrial appendages of SR patients were enrolled from two datasets, one for screening AF-related lnRNAs by ceRNA network analyses and another for validation by WGCNA. By construction of a ceRNA network, combined with a CTD database, the ceRNA axis FAM201A-miR-33a-3p-RAC3 was identified associated with AF susceptibility. Subsequently, by WGCNA for lncRNAs, two co-expression lncRNA modules were proved associated with AF and FAM201A was finally validated to be highly negatively related to AF. Collectively, lncRNA FAM201A was speculated to function, at least partly, as ceRNA to regulate RAC3 in AF susceptibility.
FAM201A refers to the lncRNA family with sequence similarity 201-member A, located in genomic 9p13.1 with 2.9 Kbp long . FAM201A was demonstrated to be involved in various diseases previously, especially cancers. In patients with lung squamous cell cancer, FAM201A was up-regulated. By regulating ATP-binding cassette transporter E1, FAM201A participated in cell proliferation, migration, invasion and influenced the survival of these patients . In tissues from non-small cell lung cancer patients, elevated expression level of FAM201A was detected related to radioresistance. For these patients, FAM201A-miR-370-EGFR was suggested as a key axis in regulation of radiotherapy sensitivity . Another in vivo experiment showed that high level expression of FAM201A involved in development of lung adenocarcinoma and down-regulation of FAM201A exerted the opposite effect . In addition to cancers, down-regulation of FAM201A was reported to be associated with Osteonecrosis of the femoral head through bioinformatics analysis and quantitative real-time polymerase chain reaction experiments . For the first time, FAM201A was revealed to be associated with AF susceptibility in the present study. In order to confirm the potential function of FAM201A on AF, both the ceRNA network analyses and WGCNA were applied. Identical results concerning the function of FAM201A on AF were achieved based on the two different methods and two different microarray datasets. Thus, these reliable results indicated that down-regulation of FAM201A may serve as potential prediction of AF susceptibility.
The underlying mechanism of FAM201A on AF could be elucidated, at least partly, through ceRNA network analyses. In our constructed ceRNA network, FAM201A regulated RAC3 by sponging miR-33a-3p. Down-regulation of FAM201A and the consequent down-regulation of RAC3 were detected to correlate with AF susceptibility. RAC3 refers to the Rac family of small guanosine triphosphatases . This is the first time that alteration of RAC3 expression from left atrial appendages of AF patients compared with SR patients has been reported. However, expression level of RAC3 was previously reported to be elevated in leukocytes from AF patients compared with controls . Actually, there is a contradiction between their study and the present study. It was believed that the different samples obtained might be responsible for the discrepancy between the findings from this previous study and our study. As stated above, different samples obtained for research could contribute to distinct gene expression .
Although we have not confirmed through experiments that FAM201A plays the role of ceRNA and regulates RAC3 in AF susceptibility. There was evidence that the ceRNA axis FAM201A-miR-33a-3p-RAC3 might increase AF susceptibility through autophagy. The fact is that FAM201A, miR-33a-3p and RAC3 all correlate with autophagy, directly or indirectly. He et al. reported that in RAC3 knockdown human umbilical vein endothelial cells, level of autophagy was detected to be much higher, which was related to inhibition of endothelial dysfunction caused by oxidized low‐density lipoprotein. The results indicated the role of RAC3 on endothelial dysfunction by down-regulating autophagy . Rubio et al. demonstrate that the expression level of RAC3 was negatively associated with autophagy, influencing chemoresistance of colorectal cancer . In a word, RAC3 was verified to have a capacity for down-regulating autophagy. Han et al. found that in primary hepatocellular cancer, the down-regulation of miR-33a-3p expression is related to metastases and poor prognosis, and confirmed that increased miR-33a-3p expression could inhibit cell migration and invasion . A large number of studies have confirmed that decreased level of cell migration and invasion are often significantly related to increased level of autophagy [42, 43]. It can be seen that the increased expression of miR-33a-3p is indirectly related to the increase of autophagy. One of the main functions of FAM201A is to regulate cell proliferation and invasion. Specifically, low expression of FAM201A can inhibit cell proliferation and invasion . There is evidence that decreased cell proliferation and invasion are related to autophagy activation [44, 45]. Therefore, the decrease of lncRNA FAM201A expression may be indirectly related to the increase of autophagy. Interestingly, high levels of autophagy was observed in both AF patients and rabbit models of atrial rapid pacing . Promotion of autophagy by overexpression of autophagy-related gene7 could lead to decreased L-type calcium channel, diminished L-type calcium current, abbreviated action potential duration, and higher AF incidence. It was also testified that up-regulated autophagy could aggravate the autophagic degradation of L-type calcium channels and the related electrical remodeling in AF . In addition to these findings, the high level of autophagy or autophagy-related gene expression were also observed elsewhere in experimental and clinical AF [47, 48]. Thus, we speculated based on the findings of our study that decreased expression of RAC3 resulting from decreased FAM201A might potentially promote autophagy and the consequent vulnerability to AF. We would indeed elucidate this mechanism in the future experiments.
There were still some limitations in the present study. First, the sample size in our study was not large enough. Although we testified the role of lncRNAs through different methods, more validations were still needed to confirm the role of some key lncRNAs in AF. Second, the expression levels of lncRNAs in the dataset used for WGCNA were relatively much lower than those of protein-coding genes. Much data about lnRNAs identified in the ceRNA network would be missed if we include all protein-coding genes in the analyses. Therefore, in order to obtain as much information about lncRNAs as possible, we performed WGCNA only for lnRNAs, which would inevitably lead to the missing information about lncRNA-mRNA interactions. Third, the GO enrichment analyses failed to provide us significant clues on the underlying mechanisms of AF related to genes identified in the present study, which could be explained by the small number of mRNAs included in the final ceRNA network used for this enrichment analyses. Most importantly, the underlying mechanism of FAM201A related to AF was speculated on previous research and our bioinformatics analysis, and further experiments were badly required to verify this theory, in order to better explain the mechanisms of AF.
In conclusion, we identified an important ceRNA axis FAM201A-miR-33a-3p-RAC3 associated with AF susceptibility through analyses of the ceRNA network. The pivotal role of FAM201A on AF was then validated by WGCNA. These findings indicated that decreased expression of FAM201A exerted an important role on AF susceptibility through down-regulating RAC3 and gave us a novel clue on further experiments about the underlying mechanisms of AF.
Availability of data and materials
The datasets (accession number: GSE41177; accession number: GSE79768) used in the present study are available in the GEO database. The dataset GSE41177 can be accessed in the website: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE41177. The dataset GSE79768 can be accessed in the website: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE79768.
Competing endogenous RNA
Weighted gene co-expression network analysis
Gene expression omnibus
Differentially expressed lncRNAs
Differentially expressed mRNAs
Comparative toxicogenomics database
Long noncoding RNA
Family with sequence similarity 201-member A
Rac family of small guanosine triphosphatases
Chugh SS, Havmoeller R, Narayanan K, Singh D, Rienstra M, Benjamin EJ, Gillum RF, Kim YH, McAnulty JH Jr, Zheng ZJ, et al. Worldwide epidemiology of atrial fibrillation: a Global Burden of Disease 2010 Study. Circulation. 2014;129(8):837–47.
Hindricks G, Potpara T, Dagres N, Arbelo E, Bax JJ, Blomstrom-Lundqvist C, Boriani G, Castella M, Dan GA, Dilaveris PE, et al. 2020 ESC guidelines for the diagnosis and management of atrial fibrillation developed in collaboration with the European Association for Cardio-Thoracic Surgery (EACTS): the task force for the diagnosis and management of atrial fibrillation of the European Society of Cardiology (ESC) Developed with the special contribution of the European Heart Rhythm Association (EHRA) of the ESC. Eur Heart J. 2021;42(5):373–498.
Chung MK, Refaat M, Shen WK, Kutyifa V, Cha YM, Di Biase L, Baranchuk A, Lampert R, Natale A, Fisher J, et al. Atrial fibrillation: JACC council perspectives. J Am Coll Cardiol. 2020;75(14):1689–713.
Roselli C, Chaffin MD, Weng LC, Aeschbacher S, Ahlberg G, Albert CM, Almgren P, Alonso A, Anderson CD, Aragam KG, et al. Multi-ethnic genome-wide association study for atrial fibrillation. Nat Genet. 2018;50(9):1225–33.
Zhang J, Huang X, Wang X, Gao Y, Liu L, Li Z, Chen X, Zeng J, Ye Z, Li G. Identification of potential crucial genes in atrial fibrillation: a bioinformatic analysis. BMC Med Genom. 2020;13(1):104.
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, Li X, Lu Y. 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, Zheng P, Pan P, Feng T, Jin L, 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, Zhai M, Huang Y, Yan KW, Dong YH, et al. LncRNA CAIF inhibits autophagy and attenuates myocardial infarction by blocking p53-mediated myocardin transcription. Nat Commun. 2018;9(1):29.
Ruan ZB, Wang F, Gongben BD, Chen GC, Zhu L. Identification of circulating lncRNA Expression profiles in patients with atrial fibrillation. Dis Markers. 2020;2020:8872142.
Cao F, Li Z, Ding WM, Yan L, Zhao QY. LncRNA PVT1 regulates atrial fibrosis via miR-128-3p-SP1-TGF-beta1-Smad axis in atrial fibrillation. Mol Med. 2019;25(1):7.
Shi J, Chen L, Chen S, Wu B, Yang K, Hu X. Circulating long noncoding RNA, GAS5, as a novel biomarker for patients with atrial fibrillation. J Clin Lab Anal. 2021;35(1):e23572.
Salmena L, Poliseno L, Tay Y, Kats L, Pandolfi PP. A ceRNA hypothesis: the Rosetta Stone of a hidden RNA language? Cell. 2011;146(3):353–8.
Du J, Li Z, Wang X, Li J, Liu D, Wang X, Wei J, Ma S, Zhang Y, Hou Y. Long noncoding RNA TCONS-00106987 promotes atrial electrical remodelling during atrial fibrillation by sponging miR-26 to regulate KCNJ2. J Cell Mol Med. 2020;24(21):12777–88.
Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinform. 2008;9:559.
Liu Y, Liang X, Wang J, Wang Y, Liu X, Shangguan W, Miao S, Liu Y, Jiao Z. Identification of atrial fibrillation-associated lncRNAs and exploration of their functions based on WGCNA and ceRNA network analyses. Gen Physiol Biophys. 2021;40(4):289–305.
Hu Z, Yang D, Tang Y, Zhang X, Wei Z, Fu H, Xu J, Zhu Z, Cai Q. 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, Shuai X, Wu C, Tao K. Identification of a 5microRNA signature and hub miRNAmRNA interactions associated with pancreatic cancer. Oncol Rep. 2019;41(1):292–300.
Yeh YH, Kuo CT, Lee YS, Lin YM, Nattel S, Tsai FC, Chen WJ. Region-specific gene expression profiles in the left atria of patients with valvular atrial fibrillation. Heart Rhythm. 2013;10(3):383–91.
Tsai FC, Lin YC, Chang SH, Chang GJ, Hsu YJ, Lin YM, Lee YS, Wang CL, Yeh YH. Differential left-to-right atria gene expression ratio in human sinus rhythm and atrial fibrillation: implications for arrhythmogenesis and thrombogenesis. Int J Cardiol. 2016;222:104–12.
Sinner MF, Tucker NR, Lunetta KL, Ozaki K, Smith JG, Trompet S, Bis JC, Lin H, Chung MK, Nielsen JB, et al. Integrating genetic, transcriptional, and functional analyses to identify 5 novel genes for atrial fibrillation. Circulation. 2014;130(15):1225–35.
Thomas AM, Cabrera CP, Finlay M, Lall K, Nobles M, Schilling RJ, Wood K, Mein CA, Barnes MR, Munroe PB, et al. Differentially expressed genes for atrial fibrillation identified by RNA sequencing from paired human left and right atrial appendages. Physiol Genom. 2019;51(8):323–32.
Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, Smyth GK. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47.
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(Database issue):D146–52.
Chou CH, Shrestha S, Yang CD, Chang NW, Lin YL, Liao KW, Huang WC, Sun TH, Tu SJ, Lee WH, 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.
Kohl M, Wiese S, Warscheid B. Cytoscape: software for visualization and analysis of biological networks. Methods Mol Biol. 2011;696:291–303.
Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–7.
Davis AP, Grondin CJ, Johnson RJ, Sciaky D, McMorran R, Wiegers J, Wiegers TC, Mattingly CJ. The comparative toxicogenomics database: update 2019. Nucleic Acids Res. 2019;47(D1):D948–54.
Humphray SJ, Oliver K, Hunt AR, Plumb RW, Loveland JE, Howe KL, Andrews TD, Searle S, Hunt SE, Scott CE, et al. DNA sequence and analysis of human chromosome 9. Nature. 2004;429(6990):369–74.
He W, Qiao ZX, Ma B. Long noncoding RNA FAM201A mediates the metastasis of lung squamous cell cancer via regulating ABCE1 expression. Eur Rev Med Pharmacol Sci. 2019;23(23):10343–53.
Liu AM, Zhu Y, Huang ZW, Lei L, Fu SZ, Chen Y. Long noncoding RNA FAM201A involves in radioresistance of non-small-cell lung cancer by enhancing EGFR expression via miR-370. Eur Rev Med Pharmacol Sci. 2019;23(13):5802–14.
Huang J, Yu Q, Zhou Y, Chu Y, Jiang F, Wang Q. FAM201A knockdown inhibits proliferation and invasion of lung adenocarcinoma cells by regulating miR-7515/GLO1 axis. J Cell Physiol. 2021;236(8):5620–32.
Huang G, Zhao G, Xia J, Wei Y, Chen F, Chen J, Shi J. FGF2 and FAM201A affect the development of osteonecrosis of the femoral head after femoral neck fracture. Gene. 2018;652:39–47.
Haataja L, Groffen J, Heisterkamp N. Characterization of RAC3, a novel member of the Rho family. J Biol Chem. 1997;272(33):20384–8.
Duzen IV, Yavuz F, Vuruskan E, Saracoglu E, Poyraz F, Cekici Y, Alici H, Goksuluk H, Candemir B, Sucu M, et al. Investigation of leukocyte RHO/ROCK gene expressions in patients with non-valvular atrial fibrillation. Exp Ther Med. 2019;18(4):2777–82.
He D, Xu L, Wu Y, Yuan Y, Wang Y, Liu Z, Zhang C, Xie W, Zhang L, Geng Z, et al. Rac3, but not Rac1, promotes ox-LDL induced endothelial dysfunction by downregulating autophagy. J Cell Physiol. 2020;235(2):1531–42.
Rubio MF, Lira MC, Rosa FD, Sambresqui AD, Salazar Guemes MC, Costas MA. RAC3 influences the chemoresistance of colon cancer cells through autophagy and apoptosis inhibition. Cancer Cell Int. 2017;17:111.
Han SY, Han HB, Tian XY, Sun H, Xue D, Zhao C, Jiang ST, He XR, Zheng WX, Wang J, et al. MicroRNA-33a-3p suppresses cell migration and invasion by directly targeting PBX3 in human hepatocellular carcinoma. Oncotarget. 2016;7(27):42461–73.
Liu L, Wang HJ, Meng T, Lei C, Yang XH, Wang QS, Jin B, Zhu JF. lncRNA GAS5 inhibits cell migration and invasion and promotes autophagy by targeting miR-222-3p via the GAS5/PTEN-signaling pathway in CRC. Mol Ther Nucleic Acids. 2019;17:644–56.
Song P, Li Y, Dong Y, Liang Y, Qu H, Qi D, Lu Y, Jin X, Guo Y, Jia Y, et al. Estrogen receptor beta inhibits breast cancer cells migration and invasion through CLDN6-mediated autophagy. J Exp Clin Cancer Res. 2019;38(1):354.
Sun R, Zhai R, Ma C, Miao W. Combination of aloin and metformin enhances the antitumor effect by inhibiting the growth and invasion and inducing apoptosis and autophagy in hepatocellular carcinoma through PI3K/AKT/mTOR pathway. Cancer Med. 2020;9(3):1141–51.
Zhu H, Diao S, Lim V, Hu L, Hu J. FAM83D inhibits autophagy and promotes proliferation and invasion of ovarian cancer cells via PI3K/AKT/mTOR pathway. Acta Biochim Biophys Sin (Shanghai). 2019;51(5):509–16.
Yuan Y, Zhao J, Gong Y, Wang D, Wang X, Yun F, Liu Z, Zhang S, Li W, Zhao X, et al. Autophagy exacerbates electrical remodeling in atrial fibrillation by ubiquitin-dependent degradation of L-type calcium channel. Cell Death Dis. 2018;9(9):873.
Shingu Y, Takada S, Yokota T, Shirakawa R, Yamada A, Ooka T, Katoh H, Kubota S, Matsui Y. Correlation between increased atrial expression of genes related to fatty acid metabolism and autophagy in patients with chronic atrial fibrillation. PLoS ONE. 2020;15(4):e0224713.
Wiersma M, Meijering RAM, Qi XY, Zhang D, Liu T, Hoogstra-Berends F, Sibon OCM, Henning RH, Nattel S, Brundel B. Endoplasmic reticulum stress is associated with autophagy and cardiomyocyte remodeling in experimental and human atrial fibrillation. J Am Heart Assoc. 2017;6(10):e006458.
This work was supported by grants obtained from the National Natural Science Foundation of China (No.82070328; 81870249). The funding body played no role in the design of the study and collection, analyses, and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. Differentially expressed lncRNAs in AF samples compared with SR samples.
. GO enrichment of DEmRNAs.
. Volcano plot of all RNA expression levels in AF samples compared to SR samples.
. Heatmap of differentially expressed mRNAs in AF and SR samples.
. Bar plot of GO enrichment of DEmRNAs.
. Differentially expressed mRNAs in AF samples compared with SR samples.
. The association between lncRNAs and miRNAs, miRNA and mRNAs in the ceRNA network.
. Crucial lncRNAs in blue module associated with AF based on WGCNA.
About this article
Cite this article
Chen, X., He, XY., Dan, Q. et al. FAM201A, a long noncoding RNA potentially associated with atrial fibrillation identified by ceRNA network analyses and WGCNA. BMC Med Genomics 15, 80 (2022). https://doi.org/10.1186/s12920-022-01232-w
- Atrial fibrillation
- Competing endogenous RNA
- Weighted gene co-expression network analysis
- Long noncoding RNA family with sequence similarity 201-member A