Comprehensive RNA-Seq Data Analysis Identifies Key mRNAs and lncRNAs in Atrial Fibrillation.
Ontology highlight
ABSTRACT:
SUBMITTER: Wu DM
PROVIDER: S-EPMC6783610 | biostudies-literature | 2019
REPOSITORIES: biostudies-literature
Publications
Comprehensive RNA-Seq Data Analysis Identifies Key mRNAs and lncRNAs in Atrial Fibrillation.
Wu Dong-Mei DM Zhou Zheng-Kun ZK Fan Shao-Hua SH Zheng Zi-Hui ZH Wen Xin X Han Xin-Rui XR Wang Shan S Wang Yong-Jian YJ Zhang Zi-Feng ZF Shan Qun Q Li Meng-Qiu MQ Hu Bin B Lu Jun J Chen Gui-Quan GQ Hong Xiao-Wu XW Zheng Yuan-Lin YL
Frontiers in genetics 20191002
Long non-coding RNAs (lncRNAs) are an emerging class of RNA species that may play a critical regulatory role in gene expression. However, the association between lncRNAs and atrial fibrillation (AF) is still not fully understood. In this study, we used RNA sequencing data to identify and quantify the both protein coding genes (PCGs) and lncRNAs. The high enrichment of these up-regulated genes in biological functions concerning response to virus and inflammatory response suggested that chronic vi ...[more]
Similar Datasets
Project description:BackgroundThe GENCODE project has collected over 10,000 human long non-coding RNA (lncRNA) genes. However, the vast majority of them remain to be functionally characterized. Computational investigation of potential functions of human lncRNA genes is helpful to guide further experimental studies on lncRNAs.ResultsIn this study, based on expression correlation between lncRNAs and protein-coding genes across 19 human normal tissues, we used the hypergeometric test to functionally annotate a single lncRNA or a set of lncRNAs with significantly enriched functional terms among the protein-coding genes that are significantly co-expressed with the lncRNA(s). The functional terms include all nodes in the Gene Ontology (GO) and 4,380 human biological pathways collected from 12 pathway databases. We successfully mapped 9,625 human lncRNA genes to GO terms and biological pathways, and then developed the first ontology-driven user-friendly web interface named lncRNA2Function, which enables researchers to browse the lncRNAs associated with a specific functional term, the functional terms associated with a specific lncRNA, or to assign functional terms to a set of human lncRNA genes, such as a cluster of co-expressed lncRNAs. The lncRNA2Function is freely available at http://mlg.hit.edu.cn/lncrna2function.ConclusionsThe LncRNA2Function is an important resource for further investigating the functions of a single human lncRNA, or functionally annotating a set of human lncRNAs of interest.
| S-EPMC4331805 | biostudies-literature
Project description:BackgroundAtrial 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.MethodsNinety-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.ResultsA 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.ConclusionsThis 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.
| S-EPMC10626701 | biostudies-literature
Project description:Circular RNAs (circRNAs) are an emerging class of RNA species that may play a critical regulatory role in gene expression control, which can serve as diagnostic biomarkers for many diseases due to their abundant, stable, and cell- or tissue-specific expression. However, the association between circRNAs and atrial fibrillation (AF) is still not clear. In this study, we used RNA sequencing data to identify and quantify the circRNAs. Differential expression analysis of the circRNAs identified 250 up- and 126 down-regulated circRNAs in AF subjects compared with healthy donors, respectively. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of the parental genes of the dysregulated circRNAs indicated that the up-regulated parental genes may participate in the process of DNA damage under oxidative stress. Furthermore, to annotate the dysregulated circRNAs, we constructed and merged the competing endogenous RNA (ceRNA) network and protein-protein interaction (PPI) network, respectively. In the merged network, 130 of 246 dysregulated circRNAs were successfully characterized by more than one pathway. Notably, the five circRNAs, including chr9:15474007-15490122, chr16:75445723-75448593, hsa_circ_0007256, chr12:56563313-56563992, and hsa_circ_0003533, showed the highest significance by the enrichment analysis, and four of them were enriched in cytokine-cytokine receptor interaction. These dysregulated circRNAs may mainly participate in biological processes of inflammatory response. In conclusion, the present study identified a set of dysregulated circRNAs, and characterized their potential functions, which may be associated with inflammatory responses in AF. To our knowledge, this is the first study to uncover the association between circRNAs and AF, which not only improves our understanding of the roles of circRNAs in AF, but also provides candidates of potentially functional circRNAs for AF researchers.
| S-EPMC6389718 | biostudies-literature
Project description:To systematically evaluate the whole-transcriptome sequencing data of cholangiocarcinoma (CHOL) to gain more insights into the transcriptomic landscape and molecular mechanism of this cancer, we performed whole-transcriptome sequencing based on the tumorous (C) and their corresponding non-tumorous adjacent to the tumors (CP) from eight CHOL patients. Subsequently, differential expression analysis was performed on the C and CP groups, followed by functional interaction prediction analysis to investigate gene-regulatory circuits in CHOL. In addition, The Cancer Genome Atlas (TCGA) for CHOL data was used to validate the results. In total, 2,895 differentially expressed messenger RNAs (dif-mRNAs), 56 differentially expressed microRNAs (dif-miRNAs), 151 differentially expressed long non-coding RNAs (dif-lncRNAs), and 110 differentially expressed circular RNAs (dif-circRNAs) were found in CHOL samples compared with controls. Enrichment analysis on those differentially expressed genes (DEGs) related to miRNA, lncRNA, and circRNA also identified the function of spliceosome. The downregulated hsa-miR-144-3p were significantly enriched in the competing endogenous RNA (ceRNA) complex network, which also included 7 upregulated and 13 downregulated circRNAs, 7 upregulated lncRNAs, and 90 upregulated and 40 downregulated mRNAs. Moreover, most of the DEGs and a few of the miRNAs (such as hsa-miR-144-3p) were successfully validated by TCGA data. The genes involved in RNA splicing and protein degradation processes and miR-144-3p may play fundamental roles in the pathogenesis of CHOL.
| S-EPMC7390861 | biostudies-literature
Project description:Recently, long non-coding RNAs (lncRNAs) have received increased research interest owing to their participation via distinct mechanisms in the biological processes of nonfunctional pituitary adenomas. However, changes in the expression of lncRNAs in gonadotrophin adenoma, which is the most common nonfunctional pituitary adenomas, have not yet been reported. In this study, we performed a genome-wide analysis of lncRNAs and mRNAs obtained from gonadotrophin adenoma patients' samples and normal pituitary tissues using RNA-seq. The differentially expressed lncRNAs and mRNAs were identified using fold-change filtering. We identified 839 lncRNAs and 1015 mRNAs as differentially expressed. Gene Ontology analysis indicated that the biological functions of differentially expressed mRNAs were related to transcription regulator activity and basic metabolic processes. Ingenuity Pathway Analysis was performed to identify 64 canonical pathways that were significantly enriched in the tumor samples. Furthermore, to investigate the potential regulatory roles of the differentially expressed lncRNAs on the mRNAs, we constructed general co-expression networks for 100 coding and 577 non-coding genes that showed significantly correlated expression patterns in tumor cohort. In particular, we built a special sub-network of co-expression involving 186 lncRNAs interacting with 15 key coding genes of the mTOR pathway, which might promote the pathogenesis of gonadotrophin tumor. This is the first study to explore the patterns of genome-wide lncRNAs expression and co-expression with mRNAs, which might contribute to the molecular pathogenesis of gonadotrophin adenoma.
| S-EPMC5354857 | biostudies-literature
Project description:BackgroundType 1 diabetes (T1D, named insulin-dependent diabetes) has a relatively rapid onset and significantly decreases life expectancy. This study is conducted to reveal the long non-coding RNA (lncRNA)-microRNA (miRNA)-mRNA regulatory axises implicated in T1D.MethodsThe gene expression profile under GSE55100 (GPL570 and GPL8786 datasets; including 12 T1D samples and 10 normal samples for each dataset) was extracted from Gene Expression Omnibus database. Using limma package, the differentially expressed mRNAs (DE-mRNAs), miRNAs (DE-miRNAs), and lncRNAs (DE-lncRNAs) between T1D and normal samples were analyzed. For the DE-mRNAs, the functional terms were enriched by DAVID tool, and the significant pathways were enriched using gene set enrichment analysis. The interactions among DE-lncRNAs, DE-miRNAs and DE-mRNAs were predicted using mirwalk and starbase. The lncRNA-miRNA-mRNA interaction network analysis was visualized by Cytoscape. The key genes in the interaction network were verified by quantitatively real-time PCR.ResultsIn comparison to normal samples, 236 DE-mRNAs, 184 DE-lncRNAs, and 45 DE-miRNAs in T1D samples were identified. For the 236 DE-mRNAs, 16 Gene Ontology (GO)_biological process (BP) terms, four GO_cellular component (CC) terms, and 57 significant pathways were enriched. A network involving 36 DE-mRNAs, 8 DE- lncRNAs, and 15 DE-miRNAs was built, such as TRG-AS1-miR-23b/miR-423-PPM1L and GAS5-miR-320a/miR-23b/miR-423-SERPINA1 regulatory axises. Quantitatively real-time PCR successfully validated the expression levels of TRG-AS1- miR-23b -PPM1L and GAS5-miR-320a- SERPINA1.ConclusionTRG-AS1-miR-23b-PPM1L and GAS5-miR-320a-SERPINA1 regulatory axises might impact the pathogenesis of T1D.
| S-EPMC7852109 | biostudies-literature
Project description:Osteoarthritis (OA) is a chronic disease common in the elderly population and imposes significant health and economic burden. Total joint replacement is the only currently available treatment but does not prevent cartilage degeneration. The molecular mechanism of OA, especially the role of inflammation in disease progression, is incompletely understood. We collected knee joint synovial tissue samples of eight OA patients and two patients with popliteal cysts (controls), measured the expression levels of lncRNAs, miRNAs, and mRNAs in these tissues by RNA-seq, and identified differentially expressed genes (DEGs) and key pathways. In the OA group, 343 mRNAs, 270 lncRNAs, and 247 miRNAs were significantly upregulated, and 232 mRNAs, 109 lncRNAs, and 157 miRNAs were significantly downregulated. mRNAs potentially targeted by lncRNAs were predicted. Nineteen overlapped miRNAs were screened based on our sample data and GSE 143514 data. Pathway enrichment and functional annotation analyses showed that the inflammation-related transcripts CHST11, ALDH1A2, TREM1, IL-1β, IL-8, CCL5, LIF, miR-146a-5p, miR-335-5p, lncRNA GAS5, LINC02288, and LOC101928134 were differentially expressed. In this study, inflammation-related DEGs and non-coding RNAs were identified in synovial samples, suggesting that competing endogenous RNAs have a role in OA. TREM1, LIF, miR146-5a, and GAS5 were identified to be OA-related genes and potential regulatory pathways. This research helps elucidate the pathogenesis of OA and identify novel therapeutic targets for this disorder.
| S-EPMC9968173 | biostudies-literature
Project description:BackgroundAdult degenerative scoliosis (ADS) is forecast to be a prevalent disabling condition in an aging society. Universally, its pathogenesis is perceived as intervertebral disc degeneration (IDD), however, a thought-provoking issue is why precisely a subset of patients with disc degeneration develop ADS. Exploring the diversities between common IDD and ADS would contribute to unraveling the etiological mechanisms of ADS. Therefore, we aimed to integrate the circRNA, lncRNA, miRNA, and mRNA expression profiles from normal adults (Normal), patients with lumbar disc herniation (LDH), and ADS by whole transcriptome sequencing, which identifies critical functional ncRNA and ceRNA networks and crosstalk between the various transcripts.MethodsThe fresh whole blood samples (n = 3/group) were collected from ADS patients, LDH patients, and healthy volunteers (Normal group), which were examined for mRNA, miRNA, lncRNA, and circRNA expression and screened for differentially expressed (DE) ncRNAs. Then, Gene Ontology (GO) and KEGG analyses were performed for gene annotation and enrichment pathways on the DE RNAs, which were constructed as a lncRNA-miRNA-mRNA network. Eventually, DE RNAs were validated by qRT-PCR targeting disc nucleus pulposus (NP) tissue in ADS and LDH group (n = 10/group).ResultsCompared to the LDH group, we identified 3322 DE mRNAs, 221 DE lncRNAs, 20 DE miRNAs, and 15 DE circRNAs in the ADS. In contrast to Normal, 21 miRNAs and 19 circRNAs were differentially expressed in the ADS. The expression of multiple differentially expressed ncRNAs was confirmed by qRT-PCR analysis to be consistent with the sequencing results. In addition, GO, and KEGG analysis demonstrated that most DE mRNAs and ncRNAs target genes are involved in various biological processes, including Endocytosis, Apoptosis, Rap1 signaling pathway, Notch signaling pathway, and others. The constructed lncRNA-miRNA-mRNA co-expression network was primarily related to angiogenesis and regulation.ConclusionBy focusing on comparing asymmetric and symmetric disc degeneration, whole-transcriptome sequencing and bioinformatics analysis systematically screened for key ncRNAs in the development of ADS, which provided an abundance of valuable candidates for the elucidation of regulatory mechanisms. The DE ncRNAs and the lncRNA-miRNA-mRNA network are intrinsically involved in the regulation of mediator and angiogenesis, which may contribute to the insight into the pathogenesis of ADS.
| S-EPMC10098162 | biostudies-literature
Project description:BackgroundRNA-Seq has become one of the most widely used applications based on next-generation sequencing technology. However, raw RNA-Seq data may have quality issues, which can significantly distort analytical results and lead to erroneous conclusions. Therefore, the raw data must be subjected to vigorous quality control (QC) procedures before downstream analysis. Currently, an accurate and complete QC of RNA-Seq data requires of a suite of different QC tools used consecutively, which is inefficient in terms of usability, running time, file usage, and interpretability of the results.ResultsWe developed a comprehensive, fast and easy-to-use QC pipeline for RNA-Seq data, RNA-QC-Chain, which involves three steps: (1) sequencing-quality assessment and trimming; (2) internal (ribosomal RNAs) and external (reads from foreign species) contamination filtering; (3) alignment statistics reporting (such as read number, alignment coverage, sequencing depth and pair-end read mapping information). This package was developed based on our previously reported tool for general QC of next-generation sequencing (NGS) data called QC-Chain, with extensions specifically designed for RNA-Seq data. It has several features that are not available yet in other QC tools for RNA-Seq data, such as RNA sequence trimming, automatic rRNA detection and automatic contaminating species identification. The three QC steps can run either sequentially or independently, enabling RNA-QC-Chain as a comprehensive package with high flexibility and usability. Moreover, parallel computing and optimizations are embedded in most of the QC procedures, providing a superior efficiency. The performance of RNA-QC-Chain has been evaluated with different types of datasets, including an in-house sequencing data, a semi-simulated data, and two real datasets downloaded from public database. Comparisons of RNA-QC-Chain with other QC tools have manifested its superiorities in both function versatility and processing speed.ConclusionsWe present here a tool, RNA-QC-Chain, which can be used to comprehensively resolve the quality control processes of RNA-Seq data effectively and efficiently.
| S-EPMC5813327 | biostudies-literature
Project description:BackgroundUsing human humoral metabolomic profiling, we can discover the diagnostic biomarkers and pathogenesis of disease. The specific characterization of atrial fibrillation (AF) subtypes with metabolomics may facilitate effective and targeted treatment, especially in early stages.ObjectivesBy investigating disturbed metabolic pathways, we could evaluate the diagnostic value of biomarkers based on metabolomics for different types of AF.MethodsA cohort of 363 patients was enrolled and divided into a discovery and validation set. Patients underwent an electrocardiogram (ECG) for suspected AF. Groups were divided as follows: healthy individuals (Control), suspected AF (Sus-AF), first diagnosed AF (Fir-AF), paroxysmal AF (Par-AF), persistent AF (Per-AF), and AF causing a cardiogenic ischemic stroke (Car-AF). Serum metabolomic profiles were determined by gas chromatography-mass spectrometry (GC-MS) and liquid chromatography-quadrupole time-of-flight mass spectrometry (LC-QTOF-MS). Metabolomic variables were analyzed with clinical information to identify relevant diagnostic biomarkers.ResultsThe metabolic disorders were characterized by 16 cross-comparisons. We focused on comparing all of the types of AF (All-AFs) plus Car-AF vs. Control, All-AFs vs. Car-AF, Par-AF vs. Control, and Par-AF vs. Per-AF. Then, 117 and 94 metabolites were identified by GC/MS and LC-QTOF-MS, respectively. The essential altered metabolic pathways during AF progression included D-glutamine and D-glutamate metabolism, glycerophospholipid metabolism, etc. For differential diagnosis, the area under the curve (AUC) of specific metabolomic biomarkers ranged from 0.8237 to 0.9890 during the discovery phase, and the predictive values in the validation cohort were 78.8-90.2%.ConclusionsSerum metabolomics is a powerful way to identify metabolic disturbances. Differences in small-molecule metabolites may serve as biomarkers for AF onset, progression, and differential diagnosis.
| S-EPMC9393302 | biostudies-literature