A Convergent Functional Genomics Analysis to Identify Biological Regulators Mediating Effects of Creatine Supplementation

Creatine (Cr) and phosphocreatine (PCr) are physiologically essential molecules for life, given they serve as rapid and localized support of energy- and mechanical-dependent processes. This evolutionary advantage is based on the action of creatine kinase (CK) isozymes that connect places of ATP synthesis with sites of ATP consumption (the CK/PCr system). Supplementation with creatine monohydrate (CrM) can enhance this system, resulting in well-known ergogenic effects and potential health or therapeutic benefits. In spite of our vast knowledge about these molecules, no integrative analysis of molecular mechanisms under a systems biology approach has been performed to date; thus, we aimed to perform for the first time a convergent functional genomics analysis to identify biological regulators mediating the effects of Cr supplementation in health and disease. A total of 35 differentially expressed genes were analyzed. We identified top-ranked pathways and biological processes mediating the effects of Cr supplementation. The impact of CrM on miRNAs merits more research. We also cautiously suggest two dose–response functional pathways (kinase- and ubiquitin-driven) for the regulation of the Cr uptake. Our functional enrichment analysis, the knowledge-based pathway reconstruction, and the identification of hub nodes provide meaningful information for future studies. This work contributes to a better understanding of the well-reported benefits of Cr in sports and its potential in health and disease conditions, although further clinical research is needed to validate the proposed mechanisms.

Download Full-text

ALGAEFUN with MARACAS, microALGAE FUNctional enrichment tool for MicroAlgae RnA-seq and Chip-seq AnalysiS

10.1101/2021.08.14.456338 ◽

2021 ◽

Author(s):

Ana B Romero-Losada ◽

Christina Arvanitidou ◽

Pedro de los Reyes ◽

Mercedes Garcia-Gonzalez ◽

Francisco J Romero-Campero

Keyword(s):

Molecular Mechanisms ◽

Enrichment Analysis ◽

Functional Enrichment ◽

Soil Bioremediation ◽

Full Potential ◽

Rna Seq ◽

Raw Data ◽

Cell Factories ◽

Gene Sets ◽

The Impact

Background: Microalgae are emerging as promising sustainable sources for biofuels, biostimulants in agriculture, soil bioremediation, feed and human nutrients. Nonetheless, the molecular mechanisms underpinning microalgae physiology and the biosynthesis of compounds of biotechnological interest are largely uncharacterized. This hinders the development of microalgae full potential as cell-factories. The recent application of omics technologies into microalgae research aims at unraveling these systems. Nevertheless, the lack of specific tools for analysing omics raw data generated from microalgae to provide biological meaningful information are hampering the impact of these technologies. The purpose of ALGAEFUN with MARACAS consists in providing researchers in microalgae with an enabling tool that will allow them to exploit transcriptomic and cistromic high-throughput sequencing data. Results: ALGAEFUN with MARACAS consists of two different tools. First, MARACAS (MicroAlgae RnA-seq and Chip-seq AnalysiS) implements a fully automatic computational pipeline receiving as input RNA-seq (RNA sequencing) or ChIP-seq (chromatin immunoprecipitation sequencing) raw data from microalgae studies. MARACAS generates sets of differentially expressed genes or lists of genomic loci for RNA-seq and ChIP-seq analysis respectively. Second, ALGAEFUN (microALGAE FUNctional enrichment tool) is a web-based application where gene sets generated from RNA-seq analysis as well as lists of genomic loci from ChIP-seq analysis can be used as input. On the one hand, it can be used to perform Gene Ontology and biological pathways enrichment analysis over gene sets. On the other hand, using the results of ChIP-seq data analysis, it identifies a set of potential target genes and analyses the distribution of the loci over gene features. Graphical representation of the results as well as tables with gene annotations are generated and can be downloaded for further analysis. Conclusions: ALGAEFUN with MARACAS provides an integrated environment for the microalgae research community that facilitates the process of obtaining relevant biological information from raw RNA-seq and ChIP-seq data. These applications are designed to assist researchers in the interpretation of gene lists and genomic loci based on functional enrichment analysis. ALGAEFUN with MARACAS is publicly available on https://greennetwork.us.es/AlgaeFUN/ .

Download Full-text

Systems Biology Approaches Reveal a Multi-stress Responsive WRKY Transcription Factor and Stress Associated Gene Co-expression Networks in Chickpea

Current Bioinformatics ◽

10.2174/1574893614666190204152500 ◽

2019 ◽

Vol 14 (7) ◽

pp. 591-601 ◽

Cited By ~ 1

Author(s):

Aravind K. Konda ◽

Parasappa R. Sabale ◽

Khela R. Soren ◽

Shanmugavadivel P. Subramaniam ◽

Pallavi Singh ◽

...

Keyword(s):

Expression Pattern ◽

Gene Networks ◽

Molecular Mechanisms ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Wrky Transcription Factor ◽

Functional Enrichment ◽

Differentially Expressed ◽

Callose Deposition ◽

Significant Probability

Background: Chickpea is a nutritional rich premier pulse crop but its production encounters setbacks due to various stresses and understanding of molecular mechanisms can be ascribed foremost importance. Objective: The investigation was carried out to identify the differentially expressed WRKY TFs in chickpea in response to herbicide stress and decipher their interacting partners. Methods: For this purpose, transcriptome wide identification of WRKY TFs in chickpea was done. Behavior of the differentially expressed TFs was compared between other stress conditions. Orthology based cofunctional gene networks were derived from Arabidopsis. Gene ontology and functional enrichment analysis was performed using Blast2GO and STRING software. Gene Coexpression Network (GCN) was constructed in chickpea using publicly available transcriptome data. Expression pattern of the identified gene network was studied in chickpea-Fusarium interactions. Results: A unique WRKY TF (Ca_08086) was found to be significantly (q value = 0.02) upregulated not only under herbicide stress but also in other stresses. Co-functional network of 14 genes, namely Ca_08086, Ca_19657, Ca_01317, Ca_20172, Ca_12226, Ca_15326, Ca_04218, Ca_07256, Ca_14620, Ca_12474, Ca_11595, Ca_15291, Ca_11762 and Ca_03543 were identified. GCN revealed 95 hub genes based on the significant probability scores. Functional annotation indicated role in callose deposition and response to chitin. Interestingly, contrasting expression pattern of the 14 network genes was observed in wilt resistant and susceptible chickpea genotypes, infected with Fusarium. Conclusion: This is the first report of identification of a multi-stress responsive WRKY TF and its associated GCN in chickpea.

Download Full-text

Structural variations in papaya genomes

BMC Genomics ◽

10.1186/s12864-021-07665-4 ◽

2021 ◽

Vol 22 (1) ◽

Author(s):

Zhenyang Liao ◽

Xunxiao Zhang ◽

Shengcheng Zhang ◽

Zhicong Lin ◽

Xingtan Zhang ◽

...

Keyword(s):

Agronomic Traits ◽

Enrichment Analysis ◽

Copy Number Variant ◽

Functional Enrichment Analysis ◽

Functional Enrichment ◽

Structural Variations ◽

Reproduction Traits ◽

Depth Study ◽

Environmental Adaptability ◽

The Impact

Abstract Background Structural variations (SVs) are a type of mutations that have not been widely detected in plant genomes and studies in animals have shown their role in the process of domestication. An in-depth study of SVs will help us to further understand the impact of SVs on the phenotype and environmental adaptability during papaya domestication and provide genomic resources for the development of molecular markers. Results We detected a total of 8083 SVs, including 5260 deletions, 552 tandem duplications and 2271 insertions with deletion being the predominant, indicating the universality of deletion in the evolution of papaya genome. The distribution of these SVs is non-random in each chromosome. A total of 1794 genes overlaps with SV, of which 1350 genes are expressed in at least one tissue. The weighted correlation network analysis (WGCNA) of these expressed genes reveals co-expression relationship between SVs-genes and different tissues, and functional enrichment analysis shows their role in biological growth and environmental responses. We also identified some domesticated SVs genes related to environmental adaptability, sexual reproduction, and important agronomic traits during the domestication of papaya. Analysis of artificially selected copy number variant genes (CNV-genes) also revealed genes associated with plant growth and environmental stress. Conclusions SVs played an indispensable role in the process of papaya domestication, especially in the reproduction traits of hermaphrodite plants. The detection of genome-wide SVs and CNV-genes between cultivated gynodioecious populations and wild dioecious populations provides a reference for further understanding of the evolution process from male to hermaphrodite in papaya.

Download Full-text

Developmental cannabidiol exposure increases anxiety and modifies genome-wide brain DNA methylation in adult female mice

Clinical Epigenetics ◽

10.1186/s13148-020-00993-4 ◽

2021 ◽

Vol 13 (1) ◽

Author(s):

Nicole M. Wanner ◽

Mathia Colwell ◽

Chelsea Drown ◽

Christopher Faulk

Keyword(s):

Dna Methylation ◽

Molecular Mechanisms ◽

Coat Color ◽

Brain Regions ◽

Functional Enrichment ◽

Positive Effects ◽

Genome Wide ◽

Nulliparous Female ◽

The Impact ◽

The Brain

Abstract Background Use of cannabidiol (CBD), the primary non-psychoactive compound found in cannabis, has recently risen dramatically, while relatively little is known about the underlying molecular mechanisms of its effects. Previous work indicates that direct CBD exposure strongly impacts the brain, with anxiolytic, antidepressant, antipsychotic, and other effects being observed in animal and human studies. The epigenome, particularly DNA methylation, is responsive to environmental input and can direct persistent patterns of gene regulation impacting phenotype. Epigenetic perturbation is particularly impactful during embryogenesis, when exogenous exposures can disrupt critical resetting of epigenetic marks and impart phenotypic effects lasting into adulthood. The impact of prenatal CBD exposure has not been evaluated; however, studies using the psychomimetic cannabinoid Δ9-tetrahydrocannabinol (THC) have identified detrimental effects on psychological outcomes in developmentally exposed adult offspring. We hypothesized that developmental CBD exposure would have similar negative effects on behavior mediated in part by the epigenome. Nulliparous female wild-type Agouti viable yellow (Avy) mice were exposed to 20 mg/kg CBD or vehicle daily from two weeks prior to mating through gestation and lactation. Coat color shifts, a readout of DNA methylation at the Agouti locus in this strain, were measured in F1 Avy/a offspring. Young adult F1 a/a offspring were then subjected to tests of working spatial memory and anxiety/compulsive behavior. Reduced-representation bisulfite sequencing was performed on both F0 and F1 cerebral cortex and F1 hippocampus to identify genome-wide changes in DNA methylation for direct and developmental exposure, respectively. Results F1 offspring exposed to CBD during development exhibited increased anxiety and improved memory behavior in a sex-specific manner. Further, while no significant coat color shift was observed in Avy/a offspring, thousands of differentially methylated loci (DMLs) were identified in both brain regions with functional enrichment for neurogenesis, substance use phenotypes, and other psychologically relevant terms. Conclusions These findings demonstrate for the first time that despite positive effects of direct exposure, developmental CBD is associated with mixed behavioral outcomes and perturbation of the brain epigenome.

Download Full-text

In Silico Identification of Contradictory Role of ADAMTS5 in Hepatocellular Carcinoma

Technology in Cancer Research & Treatment ◽

10.1177/1533033820986826 ◽

2021 ◽

Vol 20 ◽

pp. 153303382098682

Author(s):

Zhipeng Zhu ◽

Jiuhua Xu ◽

Xiaofang Wu ◽

Sihao Lin ◽

Lulu Li ◽

...

Keyword(s):

Hepatocellular Carcinoma ◽

Subgroup Analysis ◽

Prognostic Value ◽

Molecular Mechanisms ◽

Enrichment Analysis ◽

Gene Set Enrichment Analysis ◽

Clinicopathological Parameters ◽

Free Survival ◽

Cox Regression Analysis ◽

The Impact

Background: ADAMTS5 has different roles in multiple types of cancers and participates in various molecular mechanisms. However, the prognostic value of ADAMTS5 in patients with hepatocellular carcinoma (HCC) still remains unclear. We carried the study to evaluate the prognostic value and identified underlying molecular mechanisms in HCC. Methods: Firstly, the association of ADAMTS5 expression and clinicopathological parameters was evaluated by in GSE14520. Next, ADAMTS5 expression in HCC was performed using GSE14520, GSE36376, GSE76427 and The Cancer Genome Atlas (TCGA) profile. Furthermore, Kaplan-Meier analysis, Univariate and Multivariate Cox regression analysis, subgroup analysis was performed to evaluate the prognostic value of ADAMTS5 in HCC. Finally, GO enrichment analysis, gene set enrichment analysis (GSEA) and weighted gene co-expression network analysis (WGCNA) were performed to revealed underlying molecular mechanisms. Result: The expression of ADAMTS5 was positively correlated with the development of HCC. Next, high ADAMTS5 expression was significantly associated with poorer survival (all P < 0.05) and the impact of ADAMTS5 on all overall survival (OS), disease-free survival (DFS), relapse-free survival (RFS), disease specific survival (DSS) and progression free interval (PFI) was specific for HCC among other 29 cancer types. Subgroup analysis showed that ADAMTS5 overexpression was significantly associated with poorer OS in patients with HCC. Finally, ADAMTS5 might participate in the status conversion from metabolic-dominant to extracellular matrix-dominant, and the activation of ECM-related biological process might contribute to high higher mortality risk for patients with HCC. Conclusion: ADAMTS5 may play an important role in the progression of HCC, and may be considered as a novel and effective biomarker for predicting prognosis for patients with HCC.

Download Full-text

Creatine in Health and Disease

Nutrients ◽

10.3390/nu13020447 ◽

2021 ◽

Vol 13 (2) ◽

pp. 447

Author(s):

Richard B. Kreider ◽

Jeffery R. Stout

Keyword(s):

Scientific Evidence ◽

Critical Role ◽

Creatine Supplementation ◽

Cellular Metabolism ◽

Medical Evidence ◽

Therapeutic Benefits ◽

Nutritional Strategy ◽

Severity Of Injury ◽

Health And Disease

Although creatine has been mostly studied as an ergogenic aid for exercise, training, and sport, several health and potential therapeutic benefits have been reported. This is because creatine plays a critical role in cellular metabolism, particularly during metabolically stressed states, and limitations in the ability to transport and/or store creatine can impair metabolism. Moreover, increasing availability of creatine in tissue may enhance cellular metabolism and thereby lessen the severity of injury and/or disease conditions, particularly when oxygen availability is compromised. This systematic review assesses the peer-reviewed scientific and medical evidence related to creatine’s role in promoting general health as we age and how creatine supplementation has been used as a nutritional strategy to help individuals recover from injury and/or manage chronic disease. Additionally, it provides reasonable conclusions about the role of creatine on health and disease based on current scientific evidence. Based on this analysis, it can be concluded that creatine supplementation has several health and therapeutic benefits throughout the lifespan.

Download Full-text

A Systematic Analysis of Candidate Genes Associated with Nicotine Addiction

BioMed Research International ◽

10.1155/2015/313709 ◽

2015 ◽

Vol 2015 ◽

pp. 1-9 ◽

Cited By ~ 2

Author(s):

Meng Liu ◽

Xia Li ◽

Rui Fan ◽

Xinhua Liu ◽

Ju Wang

Keyword(s):

Molecular Mechanisms ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Nicotine Addiction ◽

Clustering Coefficient ◽

Functional Enrichment ◽

Moderate Degree ◽

Cancer Genes ◽

Systematic Analysis ◽

The Central Nervous System

Nicotine, as the major psychoactive component of tobacco, has broad physiological effects within the central nervous system, but our understanding of the molecular mechanism underlying its neuronal effects remains incomplete. In this study, we performed a systematic analysis on a set of nicotine addiction-related genes to explore their characteristics at network levels. We found that NAGenes tended to have a more moderate degree and weaker clustering coefficient and to be less central in the network compared to alcohol addiction-related genes or cancer genes. Further, clustering of these genes resulted in six clusters with themes in synaptic transmission, signal transduction, metabolic process, and apoptosis, which provided an intuitional view on the major molecular functions of the genes. Moreover, functional enrichment analysis revealed that neurodevelopment, neurotransmission activity, and metabolism related biological processes were involved in nicotine addiction. In summary, by analyzing the overall characteristics of the nicotine addiction related genes, this study provided valuable information for understanding the molecular mechanisms underlying nicotine addiction.

Download Full-text

Identification of Key Genes and Pathways Associated with Differences of Subchondral Bone in Osteoarthritis

10.21203/rs.3.rs-51799/v1 ◽

2020 ◽

Author(s):

Yiyuan Zhang ◽

Rongguo Yu ◽

Jiayu Zhang ◽

Eryou Feng ◽

Haiyang Wang ◽

...

Keyword(s):

Subchondral Bone ◽

Molecular Mechanisms ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Functional Enrichment ◽

Hub Genes ◽

Common Chronic Disease ◽

Drug Candidates ◽

Therapeutic Modalities ◽

Pathogenic Genes

Abstract BackgroundOsteoarthritis (OA) is a common chronic disease worldwide. Subchondral bone is an important pathological change in OA and responds more rapidly to adverse loading and events compared to cartilage. However, the pathogenic genes and pathways of subchondral bone are largely unclear.ObjectiveThis study aimed to identify signature differences in genes involved in knee lateral tibial (LT) and medial tibial (MT) plateaus of subchondral bone tissue while exploring their potential molecular mechanisms via bioinformatics analysis.MethodsFirst, the gene expression data of GSE51588 was downloaded from the GEO database. Differentially expressed genes (DEGs) between knee LT and MT were identified, and functional enrichment analyses were performed. Then, a protein-protein interactive network was constructed in order to acquire the hub genes, and modules analysis was conducted using STRING and Cytoscape for further analysis. The enriched hub genes were queried in DGIdb database to find suitable drug candidates in OA.ResultsA total of 202 DEGs (112 upregulated genes and 84 downregulated genes) were determined. In the PPI network, ten hub genes were identified. Five significant modules were identified using the MCODE plugin unit. Functional enrichment analysis revealed the most important signaling pathways. Six of the ten hub genes were targetable by a total of 35 drugs, suggesting their possible therapeutic use for OA .ConclusionsThe identified hub genes and functional enrichment pathways were implicated in the development and progression of subchondral bone in OA, thus improving our understanding of OA and offering molecular targets for future therapeutic modalities.

Download Full-text

Integrated Analysis of lncRNAs, mRNAs, and TFs to Identify Regulatory Networks Underlying MAP Infection in Cattle

Frontiers in Genetics ◽

10.3389/fgene.2021.668448 ◽

2021 ◽

Vol 12 ◽

Author(s):

Maryam Heidari ◽

Abbas Pakdel ◽

Mohammad Reza Bakhtiarizadeh ◽

Fariba Dehghanian

Keyword(s):

Regulatory Networks ◽

Molecular Mechanisms ◽

Enrichment Analysis ◽

Johne's Disease ◽

Johne’S Disease ◽

Functional Enrichment ◽

Integrated Analysis ◽

Complex Nature ◽

Hub Genes ◽

Cis And Trans

Johne’s disease is a chronic infection of ruminants that burdens dairy herds with a significant economic loss. The pathogenesis of the disease has not been revealed clearly due to its complex nature. In order to achieve deeper biological insights into molecular mechanisms involved in MAP infection resulting in Johne’s disease, a system biology approach was used. As far as is known, this is the first study that considers lncRNAs, TFs, and mRNAs, simultaneously, to construct an integrated gene regulatory network involved in MAP infection. Weighted gene coexpression network analysis (WGCNA) and functional enrichment analysis were conducted to explore coexpression modules from which nonpreserved modules had altered connectivity patterns. After identification of hub and hub-hub genes as well as TFs and lncRNAs in the nonpreserved modules, integrated networks of lncRNA-mRNA-TF were constructed, and cis and trans targets of lncRNAs were identified. Both cis and trans targets of lncRNAs were found in eight nonpreserved modules. Twenty-one of 47 nonpreserved modules showed significant biological processes related to the immune system and MAP infection. Some of the MAP infection’s related pathways in the most important nonpreserved modules comprise “positive regulation of cytokine-mediated signaling pathway,” “negative regulation of leukocyte migration,” “T-cell differentiation,” “neutrophil activation,” and “defense response.” Furthermore, several genes were identified in these modules, including SLC11A1, MAPK8IP1, HMGCR, IFNGR1, CMPK2, CORO1A, IRF1, LDLR, BOLA-DMB, and BOLA-DMA, which are potentially associated with MAP pathogenesis. This study not only enhanced our knowledge of molecular mechanisms behind MAP infection but also highlighted several promising hub and hub-hub genes involved in macrophage-pathogen interaction.

Download Full-text

Identifcation of The Ferroptosis-Related Gene Signature In Placenta of Patients With Early-Onset Preeclampsia

10.21203/rs.3.rs-618453/v1 ◽

2021 ◽

Author(s):

Nana Yang ◽

Qianghua Wang ◽

Biao Ding ◽

Yinging Gong ◽

Yue Wu ◽

...

Keyword(s):

Iron Homeostasis ◽

Molecular Mechanisms ◽

Late Onset ◽

Expression Profiles ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Functional Enrichment ◽

Differentially Expressed ◽

Ppi Network ◽

Network Analyses

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.

Download Full-text