Identification of functional lncRNAs in atrial fibrillation based on RNA sequencing

Abstract

Abstract Background Atrial fibrillation (AF) is one of the most common arrhythmia contributing to serious conditions such as stroke and heart failure. Recent studies demonstrated that long noncoding RNAs (lncRNAs) were related to cardiovascular disease. However, the molecular mechanisms of AF are not fully clear. This study intended to discover lncRNAs that are differentially expressed in AF compared with controls and evaluate the potential functions of these lncRNAs. Methods Ninety-seven patients (49 patients with AF and 48 patients without AF) were included in this study. Among these patients, leucocyte suspensions of 3 AF patients and 3 controls were sent for RNA-seq analysis to select differentially expressed lncRNA and mRNA. Different lncRNA expressions were validated in another samples (46 AF patients and 45 controls). Gene ontology (GO) enrichment analysis was conducted to annotate the function of selected mRNAs. Alternative splicing (AS) analysis was performed and a lncRNA-mRNA network was also constructed. The receiver operating characteristics (ROC) curve was used to evaluate diagnostic values. Logistic regression analysis was utilized to assess the risk or protective factor of AF. Results A total of 223 mRNAs and 105 lncRNAs were detected in AF patients compared with controls. Total 4 lncRNAs (LINC01781, AC009509.2, AL662844.3, AL662844.4) associated with AF were picked out for validation in another samples by quantitative real-time PCR (qRT-PCR), detecting that upregulated AC009509.2 and downregulated LINC01781 in AF patients. Multivariate logistic regression analysis illustrated that left atrial diameter (OR 1.201; 95% CI 1.093-1.320; P=0.000) and AC009509.2 (OR 1.732; 95% CI 1.092-2.747; P=0.020) were related to AF respectively. ROC curve showed that AC009509.2, LINC01781 and left atrial diameter (LAD) were predictors of AF. For LINC01781, the area under the curve (AUC) was 0.654 (95% CI 0.541-0.767, P=0.0113). For AC009509.2, the AUC was 0.710 (95% CI 0.599-0.822, P=0.0005). Bioinformatic methods (GO enrichment, AS analysis and lncRNA-mRNA network construction) were performed to reveal the role of lncRNAs. Conclusions This study discussed differentially expressed lncRNA and their potential interaction with mRNA in AF. LncRNA AC009509.2 could be a new potential biomarker for AF prediction.