Identification of biomarkers associated with metabolic cardiovascular disease using mRNA-SNP-miRNA regulatory network analysis

Author(s):  
Zhiyuan Fan ◽  
Wenjuan Peng ◽  
Zhiwen Wang ◽  
ling Zhang ◽  
Kuo Liu

Abstract Background: CVD is the leading cause of death in T2DM patients. However, few biomarkers have been identified to detect and diagnose CVD in the early stage of T2DM. The aim of our study was to identify the important mRNAs, micro (mi)RNAs and SNPs (single nucleotide polymorphisms) that are associated with metabolic cardiovascular disease. Materials and methods: Expression profiles and GWAS data were obtained from Gene Expression Omnibus (GEO) database. MiRNA-sequencing was conducted by Illumina HiSeq 2000 platform in T2DM patients and T2DM with CVD patients. EQTL analysis and gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were conducted. MRNA-miRNA co-expression network and mRNA-SNP-miRNA interaction network were established and visualized by Cytoscape 3.7.2.Results: In our study, we identified 56 genes and 16 miRNAs that were significantly differentially expressed. GO and KEGG analyses results indicated that B cell receptor signaling pathway and hematopoietic cell lineage were included in the biological functions of differentially expressed genes. MRNA-miRNA co-expression network and mRNA-SNP-miRNA interaction network illustrated that let-7i-5p, RASGRP3, KRT1 and CEP41 may be potential biomarkers for the early detection and diagnosis of CVD in T2DM patients.Conclusion: Our results suggested that downregulated let-7i-5p, and upregulated RASGRP3, KRT1 and CEP41 may play crucial roles in molecular mechanisms underlying the initiation and development of CVD in T2DM patients.

2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Zhiyuan Fan ◽  
Wenjuan Peng ◽  
Zhiwen Wang ◽  
Ling Zhang ◽  
Kuo Liu

Abstract Background CVD is the leading cause of death in T2DM patients. However, few biomarkers have been identified to detect and diagnose CVD in the early stage of T2DM. The aim of our study was to identify the important mRNAs, micro (mi)RNAs and SNPs (single nucleotide polymorphisms) that are associated with metabolic cardiovascular disease. Materials and methods Expression profiles and GWAS data were obtained from Gene Expression Omnibus (GEO) database. MiRNA-sequencing was conducted by Illumina HiSeq 2000 platform in T2DM patients and T2DM with CVD patients. EQTL analysis and gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were conducted. MRNA-miRNA co-expression network and mRNA-SNP-miRNA interaction network were established and visualized by Cytoscape 3.7.2. Results In our study, we identified 56 genes and 16 miRNAs that were significantly differentially expressed. KEGG analyses results indicated that B cell receptor signaling pathway and hematopoietic cell lineage were included in the biological functions of differentially expressed genes. MRNA-miRNA co-expression network and mRNA-SNP-miRNA interaction network illustrated that let-7i-5p, RASGRP3, KRT1 and CEP41 may be potential biomarkers for the early detection and diagnosis of CVD in T2DM patients. Conclusion Our results suggested that downregulated let-7i-5p, and upregulated RASGRP3, KRT1 and CEP41 may play crucial roles in molecular mechanisms underlying the initiation and development of CVD in T2DM patients.


2018 ◽  
Author(s):  
yuanshuai Fu ◽  
Zhe Xu ◽  
Zaizhong Chen ◽  
Bin Wen ◽  
Jianzhong Gao

The discus fish (Symphysodon aequifasciatus) is an ornamental fish that is well-known around the world. Phenotype investigation indicated that there are no significant differences in appearance between males and females of the discus fish. To better understand the sexual development mechanisms and obtain a high efficiency sex identification method in the artificial reproduction process of the discus fish, we constructed six cDNA libraries from three adult testes and three adult ovaries, and perform RNA-sequencing for identifying sex-biased candidate genes, microRNA (miRNA), and metabolic pathway using the Illumina Hiseq 4000. A total of 50,082 non-redundant genes (unigenes) were identified, of which 18,570 unigenes were significantly overexpressed in testes, and 11,182 unigenes were significantly overexpressed in ovaries, and 8 differentially expressed unigenes were validated by quantitative Real-Time PCR (qPCR). A total of 551 miRNAs were identified, of which 47 miRNAs were differentially expressed between testes and ovaries, and 7 differentially expressed miRNAs and one non-differential miRNA were validated by qPCR. Twenty-four of these differentially expressed miRNAs and their 15 predicted target genes constituted 41 important miRNA-mRNA interaction pairs, which may be important candidates for sex-related miRNAs and sex-related genes in the discus fish. Some of vital sex-related metabolic pathways were also identified that may play key roles in regulating gonad development of the discus fish. These results can provide important insights to better understand molecular mechanisms for sexual dimorphism in gonads development.


Vaccines ◽  
2021 ◽  
Vol 9 (12) ◽  
pp. 1427
Author(s):  
Mumdooh J. Sabir ◽  
Ross Low ◽  
Neil Hall ◽  
Majid Rasool Kamli ◽  
Md. Zubbair Malik

Cryptosporidium parvum (C. parvum) is a protozoan parasite known for cryptosporidiosis in pre-weaned calves. Animals and patients with immunosuppression are at risk of developing the disease, which can cause potentially fatal diarrhoea. The present study aimed to construct a network biology framework based on the differentially expressed genes (DEGs) of C. parvum infected subjects. In this way, the gene expression profiling analysis of C. parvum infected individuals can give us a snapshot of actively expressed genes and transcripts under infection conditions. In the present study, we have analyzed microarray data sets and compared the gene expression profiles of the patients with the different data sets of the healthy control. Using a network medicine approach to identify the most influential genes in the gene interaction network, we uncovered essential genes and pathways related to C. parvum infection. We identified 164 differentially expressed genes (109 up- and 54 down-regulated DEGs) and allocated them to pathway and gene set enrichment analysis. The results underpin the identification of seven significant hub genes with high centrality values: ISG15, MX1, IFI44L, STAT1, IFIT1, OAS1, IFIT3, RSAD2, IFITM1, and IFI44. These genes are associated with diverse biological processes not limited to host interaction, type 1 interferon production, or response to IL-gamma. Furthermore, four genes (IFI44, IFIT3, IFITM1, and MX1) were also discovered to be involved in innate immunity, inflammation, apoptosis, phosphorylation, cell proliferation, and cell signaling. In conclusion, these results reinforce the development and implementation of tools based on gene profiles to identify and treat Cryptosporidium parvum-related diseases at an early stage.


2020 ◽  
Vol 15 ◽  
Author(s):  
Yeqing Sun ◽  
Lei Chen ◽  
Yingqi Zhang ◽  
Jincheng Zhang ◽  
Shashi Ranjan Tiwari

Background: Osteoarthritis (OA), one of the most important causes leading to joint disability, was considered as an untreatable disease. A series of genes were reported to regulate the pathogenesis of OA, including microRNAs, Long non-coding RNAs and Circular RNA. So far, the expression profiles and functions of lncRNAs, mRNAs, and circRNAs in OA are not fully understood. Objective: The present study aimed to identify differently expressed genes in OA. Methods: The present study conducted RNA-seq to identify differently expressed genes in OA. Ontology (GO) analysis was used to analysis the Molecular Function and Biological Process. KEGG pathway analysis was used to perform the differentially expressed lncRNAs in biological pathways. Results: Hierarchical clustering revealed a total of 943 mRNAs, 518 lncRNAs, and 300 circRNAs were dysregulated in OA compared to normal samples. Furthermore, we constructed differentially expressed mRNAs mediated proteinprotein interaction network, differentially expressed lncRNAs mediated trans regulatory networks, and competitive endogenous RNA (ceRNA) to reveal the interaction among these genes in OA. Bioinformatics analysis revealed these dysregulated genes were involved in regulating multiple biological processes, such as wound healing, negative regulation of ossification, sister chromatid cohesion, positive regulation of interleukin-1 alpha production, sodium ion transmembrane transport, positive regulation of cell migration, and negative regulation of inflammatory response. To the best of our knowledge, this study for the first time revealed the expression pattern of mRNAs, lncRNAs and circRNAs in OA. Conclusion: This study provided novel information to validate these differentially expressed RNAs may be as possible biomarkers and targets in OA.


2021 ◽  
Vol 9 (2) ◽  
pp. 385 ◽  
Author(s):  
Zongmin Liu ◽  
Lingzhi Li ◽  
Qianwen Wang ◽  
Faizan Ahmed Sadiq ◽  
Yuankun Lee ◽  
...  

Biofilm formation has evolved as an adaptive strategy for bacteria to cope with harsh environmental conditions. Currently, little is known about the molecular mechanisms of biofilm formation in bifidobacteria. A time series transcriptome sequencing analysis of both biofilm and planktonic cells of Bifidobacterium longum FGSZY16M3 was performed to identify candidate genes involved in biofilm formation. Protein–protein interaction network analysis of 1296 differentially expressed genes during biofilm formation yielded 15 clusters of highly interconnected nodes, indicating that genes related to the SOS response (dnaK, groS, guaB, ruvA, recA, radA, recN, recF, pstA, and sufD) associated with the early stage of biofilm formation. Genes involved in extracellular polymeric substances were upregulated (epsH, epsK, efp, frr, pheT, rfbA, rfbJ, rfbP, rpmF, secY and yidC) in the stage of biofilm maturation. To further investigate the genes related to biofilm formation, weighted gene co-expression network analysis (WGCNA) was performed with 2032 transcript genes, leading to the identification of nine WGCNA modules and 133 genes associated with response to stress, regulation of gene expression, quorum sensing, and two-component system. These results indicate that biofilm formation in B. longum is a multifactorial process, involving stress response, structural development, and regulatory processes.


Hereditas ◽  
2021 ◽  
Vol 158 (1) ◽  
Author(s):  
Haoming Li ◽  
Linqing Zou ◽  
Jinhong Shi ◽  
Xiao Han

Abstract Background Alzheimer’s disease (AD) is a fatal neurodegenerative disorder, and the lesions originate in the entorhinal cortex (EC) and hippocampus (HIP) at the early stage of AD progression. Gaining insight into the molecular mechanisms underlying AD is critical for the diagnosis and treatment of this disorder. Recent discoveries have uncovered the essential roles of microRNAs (miRNAs) in aging and have identified the potential of miRNAs serving as biomarkers in AD diagnosis. Methods We sought to apply bioinformatics tools to investigate microarray profiles and characterize differentially expressed genes (DEGs) in both EC and HIP and identify specific candidate genes and pathways that might be implicated in AD for further analysis. Furthermore, we considered that DEGs might be dysregulated by miRNAs. Therefore, we investigated patients with AD and healthy controls by studying the gene profiling of their brain and blood samples to identify AD-related DEGs, differentially expressed miRNAs (DEmiRNAs), along with gene ontology (GO) analysis, KEGG pathway analysis, and construction of an AD-specific miRNA–mRNA interaction network. Results Our analysis identified 10 key hub genes in the EC and HIP of patients with AD, and these hub genes were focused on energy metabolism, suggesting that metabolic dyshomeostasis contributed to the progression of the early AD pathology. Moreover, after the construction of an miRNA–mRNA network, we identified 9 blood-related DEmiRNAs, which regulated 10 target genes in the KEGG pathway. Conclusions Our findings indicated these DEmiRNAs having the potential to act as diagnostic biomarkers at an early stage of AD.


PeerJ ◽  
2019 ◽  
Vol 7 ◽  
pp. e7821 ◽  
Author(s):  
Xiaoming Zhang ◽  
Jing Zhuang ◽  
Lijuan Liu ◽  
Zhengguo He ◽  
Cun Liu ◽  
...  

Background Cumulative evidence suggests that long non-coding RNAs (lncRNAs) play an important role in tumorigenesis. This study aims to identify lncRNAs that can serve as new biomarkers for breast cancer diagnosis or screening. Methods First, the linear fitting method was used to identify differentially expressed genes from the breast cancer RNA expression profiles in The Cancer Genome Atlas (TCGA). Next, the diagnostic value of all differentially expressed lncRNAs was evaluated using a receiver operating characteristic (ROC) curve. Then, the top ten lncRNAs with the highest diagnostic value were selected as core genes for clinical characteristics and prognosis analysis. Furthermore, core lncRNA-mRNA co-expression networks based on weighted gene co-expression network analysis (WGCNA) were constructed, and functional enrichment analysis was performed using the Database for Annotation, Visualization and Integrated Discovery (DAVID). The differential expression level and diagnostic value of core lncRNAs were further evaluated by using independent data set from Gene Expression Omnibus (GEO). Finally, the expression status and prognostic value of core lncRNAs in various tumors were analyzed based on Gene Expression Profiling Interactive Analysis (GEPIA). Results Seven core lncRNAs (LINC00478, PGM5-AS1, AL035610.1, MIR143HG, RP11-175K6.1, AC005550.4, and MIR497HG) have good single-factor diagnostic value for breast cancer. AC093850.2 has a prognostic value for breast cancer. AC005550.4 and MIR497HG can better distinguish breast cancer patients in early-stage from the advanced-stage. Low expression of MAGI2-AS3, LINC00478, AL035610.1, MIR143HG, and MIR145 may be associated with lymph node metastasis in breast cancer. Conclusion Our study provides candidate biomarkers for the diagnosis and prognosis of breast cancer, as well as a bioinformatics basis for the further elucidation of the molecular pathological mechanism of breast cancer.


2021 ◽  
Author(s):  
Nana Yang ◽  
Qianghua Wang ◽  
Biao Ding ◽  
Yinging Gong ◽  
Yue Wu ◽  
...  

Abstract Background: The accumulation of ROS resulting from upregulated levels of oxidative stress is commonly implicated in preeclampsia (PE). Ferroptosis is a novel form of iron-dependent cell death instigated by lipid peroxidation likely plays important role in PE pathogenesis. This study aims to investigate expression profiles and functions of the ferroptosis-related genes (FRGs) in early- and late-onset preeclampsia.Methods: The gene expression data and clinical information were downloaded from GEO database. The “limma” R package was used for screening differentially expressed genes. GO(Gene Ontology), Kyoto Encyclopedia of Genes and Genomes(KEGG) and protein protein interaction (PPI) network analyses were conducted to investigate the bioinformatics functions and molecular interactions of significantly different FRGs. Quantitative real-time reverse transcriptase PCR was used to verify the expression of hub FRGs in PE.Results: A total number of 4,215 DEGs were identified between EOPE and preterm cases and 3,356 DEGs were found between EOPE and LOPE subtypes. 20 significantly different FRGs were identified in EOPE, while only 3 in LOPE. Functional enrichment analysis revealed that the differentially expressed FRGs was mainly involved in EOPE and enriched in hypoxia- and iron-related pathways, such as response to hypoxia, iron homeostasis and iron ion binding process. The PPI network analysis and verification by RT-qPCR resulted in the identification of the following six interesting FRGs: FTH1, HIF1A, FTL, IREB2, MAPK8 and PLIN2. Conclusions: EOPE and LOPE owned distinct underlying molecular mechanisms and ferroptosis may be mainly implicated in pathogenesis of EOPE. Further studies are necessary for deeper inquiry into placental ferroptosis and its role in the pathogenesis of EOPE.


2020 ◽  
Author(s):  
Na Li ◽  
Ru-feng Bai ◽  
Chun Li ◽  
Li-hong Dang ◽  
Qiu-xiang Du ◽  
...  

Abstract Background: Muscle trauma frequently occurs in daily life. However, the molecular mechanisms of muscle healing, which partly depend on the extent of the damage, are not well understood. This study aimed to investigate gene expression profiles following mild and severe muscle contusion, and to provide more information about the molecular mechanisms underlying the repair process.Methods: A total of 33 rats were divided randomly into control (n = 3), mild contusion (n = 15), and severe contusion (n = 15) groups; the contusion groups were further divided into five subgroups (1, 3, 24, 48, and 168 h post-injury; n = 3 per subgroup). Then full genome microarray of RNA isolated from muscle tissue was performed to access the gene expression changes during healing process.Results: A total of 2,844 and 2,298 differentially expressed genes were identified in the mild and severe contusion groups, respectively. The analysis of the overlapping differentially expressed genes showed that there are common mechanisms of transcriptomic repair of mild and severe contusion within 48 h post-contusion. This was supported by the results of principal component analysis, hierarchical clustering, and weighted gene co‐expression network analysis of the 1,620 coexpressed genes in mildly and severely contused muscle. From these analyses, we discovered that the gene profiles in functional modules and temporal clusters were similar between the mild and severe contusion groups; moreover, the genes showed time-dependent patterns of expression, which allowed us to identify useful markers of wound age. We then performed an analysis of the functions of genes (including Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway annotation, and protein–protein interaction network analysis) in the functional modules and temporal clusters, and the hub genes in each module–cluster pair were identified. Interestingly, we found that genes downregulated within 24−48 h of the healing process were largely associated with metabolic processes, especially oxidative phosphorylation of reduced nicotinamide adenine dinucleotide phosphate, which has been rarely reported. Conclusions: These results improve our understanding of the molecular mechanisms underlying muscle repair, and provide a basis for further studies of wound age estimation.


2021 ◽  
Author(s):  
Weihao Chen ◽  
Zhifeng Li ◽  
Wei Sun ◽  
Mingxing Chu

Abstract Background:In sheep, FecB is the essential biomarker of the fertility, previous researches have provided a detailed insight on the regulation involved estrus phase and FecB in the reproductive-related tissues including hypothalamus, pituitary, and ovary. However, as the host of embryo development and connection between the ovary and the uterus, little is known about the interaction between mRNAs and lncRNAs in sheep oviduct. In the present study, RNA-Seq was performed to identify the transcriptomic profiles of mRNAs and lncRNAs in oviduct during estrus phase of sheep with FecBBB/++ genotypes.Results:In total, 21,863 lncRNAs and 43,674 mRNAs were identified, 57 DE lncRNAs and 637 DE mRNAs were revealed in the comparisons between follicular phase and luteal phase, 26 DE lncRNAs and 421 DE lncRNAs were revealed in the comparisons between FecB BB genotype and FecB ++ genotype. Functional enrichment analysis suggested that GO and KEGG terms related to reproduction such as SAGA complex, ATP-binding cassette (ABC), Nestin, and Hippo signalling pathway. DE-interaction network suggested that LNC_018420 maybe the key regulators related to embryo development in sheep oviduct.Conclusion:This was the first study to reveal the transcriptomic profiles of mRNAs and lncRNAs in the oviduct of FecB BB/++ sheep at estrus phase using RNA-Seq. Our findings can provide new understanding on the molecular mechanisms of mRNAs and lncRNAs underlying sheep embryo development and also opening new lines of investigation in sheep reproduction.


Sign in / Sign up

Export Citation Format

Share Document