1
|
Feng X, Tong L, Ma L, Mu T, Yu B, Ma R, Li J, Wang C, Zhang J, Gu Y. Mining key circRNA-associated-ceRNA networks for milk fat metabolism in cows with varying milk fat percentages. BMC Genomics 2024; 25:323. [PMID: 38561663 PMCID: PMC10983688 DOI: 10.1186/s12864-024-10252-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2023] [Accepted: 03/26/2024] [Indexed: 04/04/2024] Open
Abstract
BACKGROUND Cow milk fat is an essential indicator for evaluating and measuring milk quality and cow performance. Growing research has identified the molecular functions of circular RNAs (circRNAs) necessary for mammary gland development and lactation in mammals. METHOD The present study analyzed circRNA expression profiling data in mammary epithelial cells (MECs) from cows with highly variable milk fat percentage (MFP) using differential expression analysis and weighted gene co-expression network analysis (WGCNA). RESULTS A total of 309 differentially expressed circRNAs (DE-circRNAs) were identified in the high and low MFP groups. WGCNA analysis revealed that the pink module was significantly associated with MFP (r = - 0.85, P = 0.007). Parental genes of circRNAs in this module were enriched mainly in lipid metabolism-related signaling pathways, such as focal adhesion, ECM-receptor interaction, adherens junction and AMPK. Finally, six DE-circRNAs were screened from the pink module: circ_0010571, circ_0007797, circ_0002746, circ_0003052, circ_0004319, and circ_0012840. Among them, circ_0002746, circ_0003052, circ_0004319, and circ_0012840 had circular structures and were highly expressed in mammary tissues. Subcellular localization revealed that these four DE-circRNAs may play a regulatory role in the mammary glands of dairy cows, mainly as competitive endogenous RNAs (ceRNAs). Seven hub target genes (GNB1, GNG2, PLCB1, PLCG1, ATP6V0C, NDUFS4, and PIGH) were obtained by constructing the regulatory network of their ceRNAs and then analyzed by CytoHubba and MCODE plugins in Cytoscape. Functional enrichment analysis revealed that these genes are crucial and most probable ceRNA regulators in milk fat metabolism. CONCLUSIONS Our study identified several vital circRNAs and ceRNAs affecting milk fat synthesis, providing new research ideas and a theoretical basis for cow lactation, milk quality, and breed improvement.
Collapse
Affiliation(s)
- Xiaofang Feng
- Key Laboratory of Ruminant Molecular and Cellular Breeding, School of Agriculture, Ningxia University, 750021, Yinchuan, China
| | - Lijia Tong
- Key Laboratory of Ruminant Molecular and Cellular Breeding, School of Agriculture, Ningxia University, 750021, Yinchuan, China
| | - Lina Ma
- NingXia Academy of Agriculture and Forestry Sciences, 750002, Yinchuan, China
| | - Tong Mu
- School of Life Science, Yan'an University, 716000, Yanan, China
| | - Baojun Yu
- Key Laboratory of Ruminant Molecular and Cellular Breeding, School of Agriculture, Ningxia University, 750021, Yinchuan, China
| | - Ruoshuang Ma
- Key Laboratory of Ruminant Molecular and Cellular Breeding, School of Agriculture, Ningxia University, 750021, Yinchuan, China
| | - Jiwei Li
- Key Laboratory of Ruminant Molecular and Cellular Breeding, School of Agriculture, Ningxia University, 750021, Yinchuan, China
| | - Chuanchuan Wang
- Key Laboratory of Ruminant Molecular and Cellular Breeding, School of Agriculture, Ningxia University, 750021, Yinchuan, China
| | - Juan Zhang
- Key Laboratory of Ruminant Molecular and Cellular Breeding, School of Agriculture, Ningxia University, 750021, Yinchuan, China.
| | - Yaling Gu
- Key Laboratory of Ruminant Molecular and Cellular Breeding, School of Agriculture, Ningxia University, 750021, Yinchuan, China
| |
Collapse
|
2
|
Iyer MS, Pal A, Venkatesh KV. A Systems Biology Approach To Disentangle the Direct and Indirect Effects of Global Transcription Factors on Gene Expression in Escherichia coli. Microbiol Spectr 2023; 11:e0210122. [PMID: 36749045 PMCID: PMC10100776 DOI: 10.1128/spectrum.02101-22] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2022] [Accepted: 01/19/2023] [Indexed: 02/08/2023] Open
Abstract
Delineating the pleiotropic effects of global transcriptional factors (TFs) is critical for understanding the system-wide regulatory response in a particular environment. Currently, with the availability of genome-wide TF binding and gene expression data for Escherichia coli, several gene targets can be assigned to the global TFs, albeit inconsistently. Here, using a systematic integrated approach with emphasis on metabolism, we characterized and quantified the direct effects as well as the growth rate-mediated indirect effects of global TFs using deletion mutants of FNR, ArcA, and IHF regulators (focal TFs) under glucose fermentative conditions. This categorization enabled us to disentangle the dense connections seen within the transcriptional regulatory network (TRN) and determine the exact nature of focal TF-driven epistatic interactions with other global and pathway-specific local regulators (iTFs). We extended our analysis to combinatorial deletions of these focal TFs to determine their cross talk effects as well as conserved patterns of regulatory interactions. Moreover, we predicted with high confidence several novel metabolite-iTF interactions using inferred iTF activity changes arising from the allosteric effects of the intracellular metabolites perturbed as a result of the absence of focal TFs. Further, using compendium level computational analyses, we revealed not only the coexpressed genes regulated by these focal TFs but also the coordination of the direct and indirect target expression in the context of the economy of intracellular metabolites. Overall, this study leverages the fundamentals of TF-driven regulation, which could serve as a better template for deciphering mechanisms underlying complex phenotypes. IMPORTANCE Understanding the pleiotropic effects of global TFs on gene expression and their relevance underlying a specific response in a particular environment has been challenging. Here, we distinguish the TF-driven direct effects and growth rate-mediated indirect effects on gene expression using single- and double-deletion mutants of FNR, ArcA, and IHF regulators under anaerobic glucose fermentation. Such dissection assists us in unraveling the precise nature of interactions existing between the focal TF(s) and several other TFs, including those altered by allosteric effects of intracellular metabolites. We were able to recapitulate the previously known metabolite-TF interactions and predict novel interactions with high confidence. Furthermore, we determined that the direct and indirect gene expression have a strong connection with each other when analyzed using the coexpressed- or coregulated-gene approach. Deciphering such regulatory patterns explicitly from the metabolism point of view would be valuable in understanding other unpredicted complex regulation existing in nature.
Collapse
Affiliation(s)
- Mahesh S. Iyer
- Department of Chemical Engineering, Indian Institute of Technology Bombay, Mumbai, India
| | - Ankita Pal
- Department of Chemical Engineering, Indian Institute of Technology Bombay, Mumbai, India
| | - K. V. Venkatesh
- Department of Chemical Engineering, Indian Institute of Technology Bombay, Mumbai, India
| |
Collapse
|
3
|
González-Morelo KJ, Galán-Vásquez E, Melis F, Pérez-Rueda E, Garrido D. Structure of co-expression networks of Bifidobacterium species in response to human milk oligosaccharides. Front Mol Biosci 2023; 10:1040721. [PMID: 36776740 PMCID: PMC9908966 DOI: 10.3389/fmolb.2023.1040721] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/09/2022] [Accepted: 01/12/2023] [Indexed: 01/27/2023] Open
Abstract
Biological systems respond to environmental perturbations and a large diversity of compounds through gene interactions, and these genetic factors comprise complex networks. Experimental information from transcriptomic studies has allowed the identification of gene networks that contribute to our understanding of microbial adaptations. In this study, we analyzed the gene co-expression networks of three Bifidobacterium species in response to different types of human milk oligosaccharides (HMO) using weighted gene co-expression analysis (WGCNA). RNA-seq data obtained from Geo Datasets were obtained for Bifidobacterium longum subsp. Infantis, Bifidobacterium bifidum and Bifidobacterium longum subsp. Longum. Between 10 and 20 co-expressing modules were obtained for each dataset. HMO-associated genes appeared in the modules with more genes for B. infantis and B. bifidum, in contrast with B. longum. Hub genes were identified in each module, and in general they participated in conserved essential processes. Certain modules were differentially enriched with LacI-like transcription factors, and others with certain metabolic pathways such as the biosynthesis of secondary metabolites. The three Bifidobacterium transcriptomes showed distinct regulation patterns for HMO utilization. HMO-associated genes in B. infantis co-expressed in two modules according to their participation in galactose or N-Acetylglucosamine utilization. Instead, B. bifidum showed a less structured co-expression of genes participating in HMO utilization. Finally, this category of genes in B. longum clustered in a small module, indicating a lack of co-expression with main cell processes and suggesting a recent acquisition. This study highlights distinct co-expression architectures in these bifidobacterial genomes during HMO consumption, and contributes to understanding gene regulation and co-expression in these species of the gut microbiome.
Collapse
Affiliation(s)
- Kevin J. González-Morelo
- Department of Chemical and Bioprocess Engineering, School of Engineering, Pontificia Universidad Católica de Chile, Santiago, Chile
| | - Edgardo Galán-Vásquez
- Departamento de Ingeniería de Sistemas Computacionales y Automatización, Instituto de Investigación en Matemáticas Aplicadas y en Sistemas. Universidad Nacional Autónoma de México, Ciudad Universitaria, México City, México
| | - Felipe Melis
- Department of Chemical and Bioprocess Engineering, School of Engineering, Pontificia Universidad Católica de Chile, Santiago, Chile
| | - Ernesto Pérez-Rueda
- Instituto de Investigaciones en Matemáticas Aplicadas y en Sistemas, Universidad Nacional Autónoma de México, Unidad Académica Yucatán, Mérida, Mexico
| | - Daniel Garrido
- Department of Chemical and Bioprocess Engineering, School of Engineering, Pontificia Universidad Católica de Chile, Santiago, Chile,*Correspondence: Daniel Garrido,
| |
Collapse
|
4
|
Dovrolis N, Filidou E, Tarapatzi G, Kokkotis G, Spathakis M, Kandilogiannakis L, Drygiannakis I, Valatas V, Arvanitidis K, Karakasiliotis I, Vradelis S, Manolopoulos VG, Paspaliaris V, Bamias G, Kolios G. Co-expression of fibrotic genes in inflammatory bowel disease; A localized event? Front Immunol 2022; 13:1058237. [PMID: 36632136 PMCID: PMC9826764 DOI: 10.3389/fimmu.2022.1058237] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2022] [Accepted: 12/08/2022] [Indexed: 12/27/2022] Open
Abstract
Introduction Extracellular matrix turnover, a ubiquitous dynamic biological process, can be diverted to fibrosis. The latter can affect the intestine as a serious complication of Inflammatory Bowel Diseases (IBD) and is resistant to current pharmacological interventions. It embosses the need for out-of-the-box approaches to identify and target molecular mechanisms of fibrosis. Methods and results In this study, a novel mRNA sequencing dataset of 22 pairs of intestinal biopsies from the terminal ileum (TI) and the sigmoid of 7 patients with Crohn's disease, 6 with ulcerative colitis and 9 control individuals (CI) served as a validation cohort of a core fibrotic transcriptomic signature (FIBSig), This signature, which was identified in publicly available data (839 samples from patients and healthy individuals) of 5 fibrotic disorders affecting different organs (GI tract, lung, skin, liver, kidney), encompasses 241 genes and the functional pathways which derive from their interactome. These genes were used in further bioinformatics co-expression analyses to elucidate the site-specific molecular background of intestinal fibrosis highlighting their involvement, particularly in the terminal ileum. We also confirmed different transcriptomic profiles of the sigmoid and terminal ileum in our validation cohort. Combining the results of these analyses we highlight 21 core hub genes within a larger single co-expression module, highly enriched in the terminal ileum of CD patients. Further pathway analysis revealed known and novel inflammation-regulated, fibrogenic pathways operating in the TI, such as IL-13 signaling and pyroptosis, respectively. Discussion These findings provide a rationale for the increased incidence of fibrosis at the terminal ileum of CD patients and highlight operating pathways in intestinal fibrosis for future evaluation with mechanistic and translational studies.
Collapse
Affiliation(s)
- Nikolas Dovrolis
- Laboratory of Pharmacology, Department of Medicine, Democritus University of Thrace, Alexandroupolis, Greece,Laboratory of Biology, Department of Medicine, Democritus University of Thrace, Alexandroupolis, Greece,Individualised Medicine & Pharmacological Research Solutions Center (IMPReS), Alexandroupolis, Greece,*Correspondence: George Kolios, ; Nikolas Dovrolis,
| | - Eirini Filidou
- Laboratory of Pharmacology, Department of Medicine, Democritus University of Thrace, Alexandroupolis, Greece,Individualised Medicine & Pharmacological Research Solutions Center (IMPReS), Alexandroupolis, Greece
| | - Gesthimani Tarapatzi
- Laboratory of Pharmacology, Department of Medicine, Democritus University of Thrace, Alexandroupolis, Greece,Individualised Medicine & Pharmacological Research Solutions Center (IMPReS), Alexandroupolis, Greece
| | - Georgios Kokkotis
- Gastrointestinal (GI) Unit, 3 Department of Internal Medicine, Sotiria Hospital, National and Kapodistrian University of Athens, Athens, Greece
| | - Michail Spathakis
- Laboratory of Pharmacology, Department of Medicine, Democritus University of Thrace, Alexandroupolis, Greece,Individualised Medicine & Pharmacological Research Solutions Center (IMPReS), Alexandroupolis, Greece
| | - Leonidas Kandilogiannakis
- Laboratory of Pharmacology, Department of Medicine, Democritus University of Thrace, Alexandroupolis, Greece,Individualised Medicine & Pharmacological Research Solutions Center (IMPReS), Alexandroupolis, Greece
| | - Ioannis Drygiannakis
- Gastroenterology and Hepatology Research Laboratory, Medical School, University of Crete, Heraklion, Greece
| | - Vassilis Valatas
- Laboratory of Pharmacology, Department of Medicine, Democritus University of Thrace, Alexandroupolis, Greece,Gastroenterology and Hepatology Research Laboratory, Medical School, University of Crete, Heraklion, Greece
| | - Konstantinos Arvanitidis
- Laboratory of Pharmacology, Department of Medicine, Democritus University of Thrace, Alexandroupolis, Greece,Individualised Medicine & Pharmacological Research Solutions Center (IMPReS), Alexandroupolis, Greece
| | - Ioannis Karakasiliotis
- Laboratory of Biology, Department of Medicine, Democritus University of Thrace, Alexandroupolis, Greece
| | - Stergios Vradelis
- Second Department of Internal Medicine, University Hospital of Alexandroupolis, Democritus University of Thrace, Alexandroupolis, Greece
| | - Vangelis G. Manolopoulos
- Laboratory of Pharmacology, Department of Medicine, Democritus University of Thrace, Alexandroupolis, Greece,Individualised Medicine & Pharmacological Research Solutions Center (IMPReS), Alexandroupolis, Greece
| | | | - Giorgos Bamias
- Gastrointestinal (GI) Unit, 3 Department of Internal Medicine, Sotiria Hospital, National and Kapodistrian University of Athens, Athens, Greece
| | - George Kolios
- Laboratory of Pharmacology, Department of Medicine, Democritus University of Thrace, Alexandroupolis, Greece,Individualised Medicine & Pharmacological Research Solutions Center (IMPReS), Alexandroupolis, Greece,*Correspondence: George Kolios, ; Nikolas Dovrolis,
| |
Collapse
|
5
|
Auger S, Mournetas V, Chiapello H, Loux V, Langella P, Chatel JM. Gene co-expression network analysis of the human gut commensal bacterium Faecalibacterium prausnitzii in R-Shiny. PLoS One 2022; 17:e0271847. [PMID: 36399439 PMCID: PMC9674144 DOI: 10.1371/journal.pone.0271847] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2022] [Accepted: 11/04/2022] [Indexed: 11/19/2022] Open
Abstract
Faecalibacterium prausnitzii is abundant in the healthy human intestinal microbiota, and the absence or scarcity of this bacterium has been linked with inflammatory diseases and metabolic disorders. F. prausnitzii thus shows promise as a next-generation probiotic for use in restoring the balance of the gut microbial flora and, due to its strong anti-inflammatory properties, for the treatment of certain pathological conditions. However, very little information is available about gene function and regulation in this species. Here, we utilized a systems biology approach—weighted gene co-expression network analysis (WGCNA)–to analyze gene expression in three publicly available RNAseq datasets from F. prausnitzii strain A2-165, all obtained in different laboratory conditions. The co-expression network was then subdivided into 24 co-expression gene modules. A subsequent enrichment analysis revealed that these modules are associated with different kinds of biological processes, such as arginine, histidine, cobalamin, or fatty acid metabolism as well as bacteriophage function, molecular chaperones, stress response, or SOS response. Some genes appeared to be associated with mechanisms of protection against oxidative stress and could be essential for F. prausnitzii’s adaptation and survival under anaerobic laboratory conditions. Hub and bottleneck genes were identified by analyses of intramodular connectivity and betweenness, respectively; this highlighted the high connectivity of genes located on mobile genetic elements, which could promote the genetic evolution of F. prausnitzii within its ecological niche. This study provides the first exploration of the complex regulatory networks in F. prausnitzii, and all of the “omics” data are available online for exploration through a graphical interface at https://shiny.migale.inrae.fr/app/faeprau.
Collapse
Affiliation(s)
- Sandrine Auger
- Université Paris-Saclay, INRAE, AgroParisTech, Micalis Institute, Jouy-en-Josas, France
- * E-mail: (SA); (VM)
| | - Virginie Mournetas
- ADLIN Science, Pépinière « Genopole Entreprises », Evry, France
- * E-mail: (SA); (VM)
| | | | - Valentin Loux
- Université Paris-Saclay, INRAE, MaIAGE, Jouy-en-Josas, France
- Université Paris-Saclay, INRAE, BioinfOmics, MIGALE Bioinformatics Facility, Jouy-en-Josas, France
| | - Philippe Langella
- Université Paris-Saclay, INRAE, AgroParisTech, Micalis Institute, Jouy-en-Josas, France
| | - Jean-Marc Chatel
- Université Paris-Saclay, INRAE, AgroParisTech, Micalis Institute, Jouy-en-Josas, France
| |
Collapse
|
6
|
Taboada-Castro H, Gil J, Gómez-Caudillo L, Escorcia-Rodríguez JM, Freyre-González JA, Encarnación-Guevara S. Rhizobium etli CFN42 proteomes showed isoenzymes in free-living and symbiosis with a different transcriptional regulation inferred from a transcriptional regulatory network. Front Microbiol 2022; 13:947678. [PMID: 36312930 PMCID: PMC9611204 DOI: 10.3389/fmicb.2022.947678] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2022] [Accepted: 09/05/2022] [Indexed: 11/13/2022] Open
Abstract
A comparative proteomic study at 6 h of growth in minimal medium (MM) and bacteroids at 18 days of symbiosis of Rhizobium etli CFN42 with the Phaseolus vulgaris leguminous plant was performed. A gene ontology classification of proteins in MM and bacteroid, showed 31 and 10 pathways with higher or equal than 30 and 20% of proteins with respect to genome content per pathway, respectively. These pathways were for energy and environmental compound metabolism, contributing to understand how Rhizobium is adapted to the different conditions. Metabolic maps based on orthology of the protein profiles, showed 101 and 74 functional homologous proteins in the MM and bacteroid profiles, respectively, which were grouped in 34 different isoenzymes showing a great impact in metabolism by covering 60 metabolic pathways in MM and symbiosis. Taking advantage of co-expression of transcriptional regulators (TF’s) in the profiles, by selection of genes whose matrices were clustered with matrices of TF’s, Transcriptional Regulatory networks (TRN´s) were deduced by the first time for these metabolic stages. In these clustered TF-MM and clustered TF-bacteroid networks, containing 654 and 246 proteins, including 93 and 46 TFs, respectively, showing valuable information of the TF’s and their regulated genes with high stringency. Isoenzymes were specific for adaptation to the different conditions and a different transcriptional regulation for MM and bacteroid was deduced. The parameters of the TRNs of these expected biological networks and biological networks of E. coli and B. subtilis segregate from the random theoretical networks. These are useful data to design experiments on TF gene–target relationships for bases to construct a TRN.
Collapse
Affiliation(s)
- Hermenegildo Taboada-Castro
- Proteomics Laboratory, Program of Functional Genomics of Prokaryotes, Center for Genomic Sciences, National Autonomous University of Mexico, Cuernavaca, Morelos, Mexico
| | - Jeovanis Gil
- Division of Oncology, Section for Clinical Chemistry, Department of Translational Medicine, Lund University, Lund, Sweden
| | - Leopoldo Gómez-Caudillo
- Proteomics Laboratory, Program of Functional Genomics of Prokaryotes, Center for Genomic Sciences, National Autonomous University of Mexico, Cuernavaca, Morelos, Mexico
| | - Juan Miguel Escorcia-Rodríguez
- Regulatory Systems Biology Research Group, Program of Systems Biology, Center for Genomic Sciences, National Autonomous University of Mexico, Mexico City, Mexico
| | - Julio Augusto Freyre-González
- Regulatory Systems Biology Research Group, Program of Systems Biology, Center for Genomic Sciences, National Autonomous University of Mexico, Mexico City, Mexico
| | - Sergio Encarnación-Guevara
- Proteomics Laboratory, Program of Functional Genomics of Prokaryotes, Center for Genomic Sciences, National Autonomous University of Mexico, Cuernavaca, Morelos, Mexico
- *Correspondence: Sergio Encarnacion Guevara,
| |
Collapse
|
7
|
Mu T, Hu H, Ma Y, Yang C, Feng X, Wang Y, Liu J, Yu B, Zhang J, Gu Y. Identification of critical lncRNAs for milk fat metabolism in dairy cows using WGCNA and the construction of a ceRNAs network. Anim Genet 2022; 53:740-760. [PMID: 36193627 DOI: 10.1111/age.13249] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2022] [Revised: 07/26/2022] [Accepted: 08/02/2022] [Indexed: 11/27/2022]
Abstract
As key regulators, long non-coding RNAs (lncRNAs) play a crucial role in the ruminant mammary gland. However, the function of lncRNAs in milk fat synthesis from dairy cows is largely unknown. In this study, we used the weighted gene co-expression network analysis (WGCNA) to comprehensive analyze the expression profile data of lncRNAs from the group's previous Illumina PE150 sequencing results based on bovine mammary epithelial cells from high- and low-milk-fat-percentage (MFP) cows, and identify core_lncRNAs significantly associated with MFP by module membership (MM) and gene significance (GS). Functional enrichment analysis (Gene Ontology, Kyoto Encyclopedia of Genes and Genomes) of core_lncRNA target genes (co-localization and co-expression) was performed to screen potential lncRNAs regulating milk fat metabolism and further construct an interactive regulatory network of lipid metabolism-related competing endogenous RNAs (ceRNAs). A total of 4876 lncRNAs were used to construct the WGCNA. The MEdarkturquoise module among the 19 modules obtained was significantly associated with MFP (r = 0.78, p-value <0.05) and contained 64 core_lncRNAs (MM > 0.8, GS > 0.4). Twenty-four lipid metabolism-related lncRNAs were identified by core_lncRNA target gene enrichment analysis. TCONS_00054233, TCONS_00152292, TCONS_00048619, TCONS_00033839, TCONS_00153791 and TCONS_00074642 were key candidate lncRNAs for regulating milk fat synthesis. The 22 ceRNAs most likely to be involved in milk fat metabolism were constructed by interaction network analysis, and TCONS_00133813 and bta-miR-2454-5p were located at the network's core. TCONS_00133813_bta-miR-2454-5p_TNFAIP3, TCONS_00133813_bta-miR-2454-5p_ARRB1 and TCONS_00133813_bta-miR-2454-5p_PIK3R1 are key candidate ceRNAs associated with milk fat metabolism. This study provides a framework for the co-expression module of MFP-related lncRNAs in ruminants, identifies several major lncRNAs and ceRNAs that influence milk fat synthesis, and provides a new understanding of the complex biology of milk fat synthesis in dairy cows.
Collapse
Affiliation(s)
- Tong Mu
- School of Agriculture, Ningxia University, Yinchuan, China
| | - Honghong Hu
- School of Agriculture, Ningxia University, Yinchuan, China
| | - Yanfen Ma
- School of Agriculture, Ningxia University, Yinchuan, China.,Key Laboratory of Ruminant Molecular and Cellular Breeding, Ningxia University, Yinchuan, China
| | - Chaoyun Yang
- School of Agriculture, Ningxia University, Yinchuan, China
| | - Xiaofang Feng
- School of Agriculture, Ningxia University, Yinchuan, China
| | - Ying Wang
- School of Agriculture, Ningxia University, Yinchuan, China
| | - Jiamin Liu
- School of Agriculture, Ningxia University, Yinchuan, China
| | - Baojun Yu
- School of Agriculture, Ningxia University, Yinchuan, China
| | - Juan Zhang
- School of Agriculture, Ningxia University, Yinchuan, China
| | - Yaling Gu
- School of Agriculture, Ningxia University, Yinchuan, China
| |
Collapse
|
8
|
Kuang J, Buchon N, Michel K, Scoglio C. A global [Formula: see text] gene co-expression network constructed from hundreds of experimental conditions with missing values. BMC Bioinformatics 2022; 23:170. [PMID: 35534830 PMCID: PMC9082846 DOI: 10.1186/s12859-022-04697-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2022] [Accepted: 04/25/2022] [Indexed: 11/21/2022] Open
Abstract
BACKGROUND Gene co-expression networks (GCNs) can be used to determine gene regulation and attribute gene function to biological processes. Different high throughput technologies, including one and two-channel microarrays and RNA-sequencing, allow evaluating thousands of gene expression data simultaneously, but these methodologies provide results that cannot be directly compared. Thus, it is complex to analyze co-expression relations between genes, especially when there are missing values arising for experimental reasons. Networks are a helpful tool for studying gene co-expression, where nodes represent genes and edges represent co-expression of pairs of genes. RESULTS In this paper, we establish a method for constructing a gene co-expression network for the Anopheles gambiae transcriptome from 257 unique studies obtained with different methodologies and experimental designs. We introduce the sliding threshold approach to select node pairs with high Pearson correlation coefficients. The resulting network, which we name AgGCN1.0, is robust to random removal of conditions and has similar characteristics to small-world and scale-free networks. Analysis of network sub-graphs revealed that the core is largely comprised of genes that encode components of the mitochondrial respiratory chain and the ribosome, while different communities are enriched for genes involved in distinct biological processes. CONCLUSION Analysis of the network reveals that both the architecture of the core sub-network and the network communities are based on gene function, supporting the power of the proposed method for GCN construction. Application of network science methodology reveals that the overall network structure is driven to maximize the integration of essential cellular functions, possibly allowing the flexibility to add novel functions.
Collapse
Affiliation(s)
- Junyao Kuang
- Department of Electrical and Computer Engineering, Kansas State University, Manhattan, KS 66506 USA
| | - Nicolas Buchon
- Department of Entomology, Cornell Institute of Host-Microbe Interactions and Disease, Cornell University, Ithaca, NY 14853 USA
| | - Kristin Michel
- Division of Biology, Kansas State University, Manhattan, KS 66506 USA
| | - Caterina Scoglio
- Department of Electrical and Computer Engineering, Kansas State University, Manhattan, KS 66506 USA
| |
Collapse
|
9
|
Mu T, Hu H, Ma Y, Wen H, Yang C, Feng X, Wen W, Zhang J, Gu Y. Identifying key genes in milk fat metabolism by weighted gene co-expression network analysis. Sci Rep 2022; 12:6836. [PMID: 35477736 PMCID: PMC9046402 DOI: 10.1038/s41598-022-10435-1] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2021] [Accepted: 03/21/2022] [Indexed: 12/13/2022] Open
Abstract
Milk fat is the most important and energy-rich substance in milk, and its content and composition are important reference elements in the evaluation of milk quality. However, the current identification of valuable candidate genes affecting milk fat is limited. IlluminaPE150 was used to sequence bovine mammary epithelial cells (BMECs) with high and low milk fat rates (MFP), the weighted gene co-expression network (WGCNA) was used to analyze mRNA expression profile data in this study. As a result, a total of 10,310 genes were used to construct WGCNA, and the genes were classified into 18 modules. Among them, violet (r = 0.74), yellow (r = 0.75) and darkolivegreen (r = − 0.79) modules were significantly associated with MFP, and 39, 181, 75 hub genes were identified, respectively. Combining enrichment analysis and differential genes (DEs), we screened five key candidate DEs related to lipid metabolism, namely PI4K2A, SLC16A1, ATP8A2, VEGFD and ID1, respectively. Relative to the small intestine, liver, kidney, heart, ovary and uterus, the gene expression of PI4K2A is the highest in mammary gland, and is significantly enriched in GO terms and pathways related to milk fat metabolism, such as monocarboxylic acid transport, phospholipid transport, phosphatidylinositol signaling system, inositol phosphate metabolism and MAPK signaling pathway. This study uses WGCNA to form an overall view of MFP, providing a theoretical basis for identifying potential pathways and hub genes that may be involved in milk fat synthesis.
Collapse
Affiliation(s)
- Tong Mu
- School of Agriculture, Ningxia University, Yinchuan, 750021, China
| | - Honghong Hu
- School of Agriculture, Ningxia University, Yinchuan, 750021, China
| | - Yanfen Ma
- School of Agriculture, Ningxia University, Yinchuan, 750021, China.,Key Laboratory of Ruminant Molecular and Cellular Breeding, Ningxia Hui Autonomous Region, Ningxia University, Yinchuan, 750021, China
| | - Huiyu Wen
- Maosheng Pasture of He Lanshan in Ningxia State Farm, Yinchuan, 750001, China
| | - Chaoyun Yang
- School of Agriculture, Ningxia University, Yinchuan, 750021, China
| | - Xiaofang Feng
- School of Agriculture, Ningxia University, Yinchuan, 750021, China
| | - Wan Wen
- Animal Husbandry Extension Station, Yinchuan, 750001, China
| | - Juan Zhang
- School of Agriculture, Ningxia University, Yinchuan, 750021, China
| | - Yaling Gu
- School of Agriculture, Ningxia University, Yinchuan, 750021, China.
| |
Collapse
|
10
|
Qi F, Zhang W, Huang J, Fu L, Zhao J. Single-Cell RNA Sequencing Analysis of the Immunometabolic Rewiring and Immunopathogenesis of Coronavirus Disease 2019. Front Immunol 2021; 12:651656. [PMID: 33936072 PMCID: PMC8079812 DOI: 10.3389/fimmu.2021.651656] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2021] [Accepted: 03/29/2021] [Indexed: 12/16/2022] Open
Abstract
Although immune dysfunction is a key feature of coronavirus disease 2019 (COVID-19), the metabolism-related mechanisms remain elusive. Here, by reanalyzing single-cell RNA sequencing data, we delineated metabolic remodeling in peripheral blood mononuclear cells (PBMCs) to elucidate the metabolic mechanisms that may lead to the progression of severe COVID-19. After scoring the metabolism-related biological processes and signaling pathways, we found that mono-CD14+ cells expressed higher levels of glycolysis-related genes (PKM, LDHA and PKM) and PPP-related genes (PGD and TKT) in severe patients than in mild patients. These genes may contribute to the hyperinflammation in mono-CD14+ cells of patients with severe COVID-19. The mono-CD16+ cell population in COVID-19 patients showed reduced transcription levels of genes related to lysine degradation (NSD1, KMT2E, and SETD2) and elevated transcription levels of genes involved in OXPHOS (ATP6V1B2, ATP5A1, ATP5E, and ATP5B), which may inhibit M2-like polarization. Plasma cells also expressed higher levels of the OXPHOS gene ATP13A3 in COVID-19 patients, which was positively associated with antibody secretion and survival of PCs. Moreover, enhanced glycolysis or OXPHOS was positively associated with the differentiation of memory B cells into plasmablasts or plasma cells. This study comprehensively investigated the metabolic features of peripheral immune cells and revealed that metabolic changes exacerbated inflammation in monocytes and promoted antibody secretion and cell survival in PCs in COVID-19 patients, especially those with severe disease.
Collapse
Affiliation(s)
- Furong Qi
- Institute for Hepatology, National Clinical Research Center for Infectious Disease, Shenzhen Third People's Hospital, The Second Affiliated Hospital, School of Medicine, Southern University of Science and Technology, Shenzhen, China.,Shenzhen Research Center for Communicable Disease Diagnosis and Treatment of Chinese Academy of Medical Science, Shenzhen, China
| | - Wenbo Zhang
- Trinity School of Durham and Chapel Hill, Durham, NC, United States
| | - Jialu Huang
- Electronic and Computer Engineering, China North Vehicle Research Institute, Beijing, China
| | - Lili Fu
- Center for Life Sciences, Tsinghua University, Beijing, China
| | - Jinfang Zhao
- Center for Life Sciences, Tsinghua University, Beijing, China
| |
Collapse
|
11
|
Zhao W, Wang M, Wang C, Liu Y, Liu H, Luo S. RACGAP1 is transcriptionally regulated by E2F3, and its depletion leads to mitotic catastrophe in esophageal squamous cell carcinoma. ANNALS OF TRANSLATIONAL MEDICINE 2020; 8:950. [PMID: 32953750 PMCID: PMC7475413 DOI: 10.21037/atm-20-2901] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Abstract
Background RACGAP1 has significant involvement in tumorigenesis of cancers, including liver cancer, stomach cancer, and colon cancer. However, the role and the exact mechanism of RACGAP1 in esophageal squamous cell carcinoma (ESCC) has not been explored. Methods QPCR and Western blots analysis was performed to analyze the expression of RACGAP1 in ESCC. MTT assays and colony formation assays were performed to explore the functional role of RACGAP1 in ESCC. Cell cycle analysis and immunofluorescence assays were used to investigate the function of RACGAP1 involvement in mitotic catastrophe. At last, we conducted the public datasets mining to explore the expression status and prognosis value of RACGAP1 as well as the correlation between RACGAP1 and E2F3 in various cancers. Results The high abnormal expression of RACGAP1 is observed in ESCC and associated with worse clinical outcomes of patients with ESCC. RACGAP1, a novel cell cycle associated gene regulated by E2F3, acts as an oncogenic driver in ESCC cell lines. Notably, for the first time, RACGAP1 depletion induced severe mitotic catastrophe, followed by massive cell death. Conclusions Our findings showed the essential role of RACGAP1 in ESCC cancer cell survival and the therapeutic potential of RACGAP1 as a molecular target for ESCC.
Collapse
Affiliation(s)
- Weifeng Zhao
- Department of Medical Oncology, the Affiliated Tumor Hospital of Zhengzhou University, Henan Cancer Hospital, Zhengzhou, China.,Department of Oncology, People's Hospital of Zhengzhou University, Henan Provincial People's Hospital, Zhengzhou, China
| | - Mengyao Wang
- Radiation Oncology Department, Affiliated Cancer Hospital and Institute of Guangzhou Medical University, Guangzhou, China
| | - Chaojie Wang
- Department of Oncology, People's Hospital of Zhengzhou University, Henan Provincial People's Hospital, Zhengzhou, China
| | - Yingjun Liu
- Department of General Surgery, the Affiliated Tumor Hospital of Zhengzhou University, Henan Cancer Hospital, Zhengzhou, China
| | - Huimin Liu
- Department of Oncology, People's Hospital of Zhengzhou University, Henan Provincial People's Hospital, Zhengzhou, China
| | - Suxia Luo
- Department of Medical Oncology, the Affiliated Tumor Hospital of Zhengzhou University, Henan Cancer Hospital, Zhengzhou, China
| |
Collapse
|