1
|
Fonseca PADS, Suarez-Vega A, Esteban-Blanco C, Marina H, Pelayo R, Gutiérrez-Gil B, Arranz JJ. Integration of epigenomic and genomic data to predict residual feed intake and the feed conversion ratio in dairy sheep via machine learning algorithms. BMC Genomics 2025; 26:313. [PMID: 40165084 PMCID: PMC11956460 DOI: 10.1186/s12864-025-11520-1] [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: 11/18/2024] [Accepted: 03/24/2025] [Indexed: 04/02/2025] Open
Abstract
BACKGROUND Feed efficiency (FE) is an essential trait in livestock species because of the constant demand to increase the productivity and sustainability of livestock production systems. A better understanding of the biological mechanisms associated with FEs might help improve the estimation and selection of superior animals. In this work, differentially methylated regions (DMRs) were identified via genome-wide bisulfite sequencing (GWBS) by comparing the DNA methylation profiles of milk somatic cells from dairy ewes that were divergent in terms of residual feed intake. The DMRs were identified by comparing divergent groups for residual feed intake (RFI), the feed conversion ratio (FCR), and the consensus between both metrics (Cons). Additionally, the predictive performance of these DMRs and genetic variants mapped within these regions was evaluated via three machine learning (ML) models (xgboost, random forest (RF), and multilayer feedforward artificial neural network (deeplearning)). The average performance of each model was based on the root mean squared error (RMSE) and squared Spearman correlation (rho2). Finally, the best model for each scenario was selected on the basis of the highest ratio between rho2 and RMSE. RESULTS In total, 12,257, 9,328, and 6,723 genes were annotated for DMRs detected in the RFI, FCR, and Cons groups, respectively. These genes are associated with important pathways for regulating FE in dairy sheep, such as protein digestion and absorption, hormone synthesis and secretion, control of energy availability, cellular signaling, and feed behavior pathways. With respect to the ML predictions, the smallest mean RMSE (0.17) was obtained using RF, which was used to predict RFI. The highest mean rho2 (0.20) was obtained when the RFI was predicted via the mean methylation within the DMRs identified, the consensus groups were compared, and the genetic variants mapped within these DMRs were included. The best overall models were obtained for the prediction of RFI using the DMRs obtained in the comparison of RFI groups (RMSE = 0.10, rho2 = 0.86) using xgboost and the DMRs plus the genetic variants identified via the Cons groups (RMSE = 0.07, rho2 = 0.62) using RF. CONCLUSIONS The results provide new insights into the biological mechanisms associated with FE and the control of these processes through epigenetic mechanisms. Additionally, the potential use of epigenetic information as a biomarker for the prediction of FE can be suggested based on the obtained results.
Collapse
Affiliation(s)
- Pablo Augusto de Souza Fonseca
- Departamento de ProducciĂłn Animal, Facultad de Veterinaria, Universidad de LeĂłn, Campus de Vegazana, Leon, 24007, Spain
| | - Aroa Suarez-Vega
- Departamento de ProducciĂłn Animal, Facultad de Veterinaria, Universidad de LeĂłn, Campus de Vegazana, Leon, 24007, Spain
| | - Cristina Esteban-Blanco
- Departamento de ProducciĂłn Animal, Facultad de Veterinaria, Universidad de LeĂłn, Campus de Vegazana, Leon, 24007, Spain
| | - Héctor Marina
- Departamento de ProducciĂłn Animal, Facultad de Veterinaria, Universidad de LeĂłn, Campus de Vegazana, Leon, 24007, Spain
| | - RocĂo Pelayo
- Departamento de ProducciĂłn Animal, Facultad de Veterinaria, Universidad de LeĂłn, Campus de Vegazana, Leon, 24007, Spain
| | - Beatriz Gutiérrez-Gil
- Departamento de ProducciĂłn Animal, Facultad de Veterinaria, Universidad de LeĂłn, Campus de Vegazana, Leon, 24007, Spain
| | - Juan-José Arranz
- Departamento de ProducciĂłn Animal, Facultad de Veterinaria, Universidad de LeĂłn, Campus de Vegazana, Leon, 24007, Spain.
| |
Collapse
|
2
|
Ibrahim H, Balboa D, Saarimäki-Vire J, Montaser H, Dyachok O, Lund PE, Omar-Hmeadi M, Kvist J, Dwivedi OP, Lithovius V, Barsby T, Chandra V, Eurola S, Ustinov J, Tuomi T, Miettinen PJ, Barg S, Tengholm A, Otonkoski T. RFX6 haploinsufficiency predisposes to diabetes through impaired beta cell function. Diabetologia 2024; 67:1642-1662. [PMID: 38743124 PMCID: PMC11343796 DOI: 10.1007/s00125-024-06163-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] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 12/08/2023] [Accepted: 03/21/2024] [Indexed: 05/16/2024]
Abstract
AIMS/HYPOTHESIS Regulatory factor X 6 (RFX6) is crucial for pancreatic endocrine development and differentiation. The RFX6 variant p.His293LeufsTer7 is significantly enriched in the Finnish population, with almost 1:250 individuals as a carrier. Importantly, the FinnGen study indicates a high predisposition for heterozygous carriers to develop type 2 and gestational diabetes. However, the precise mechanism of this predisposition remains unknown. METHODS To understand the role of this variant in beta cell development and function, we used CRISPR technology to generate allelic series of pluripotent stem cells. We created two isogenic stem cell models: a human embryonic stem cell model; and a patient-derived stem cell model. Both were differentiated into pancreatic islet lineages (stem-cell-derived islets, SC-islets), followed by implantation in immunocompromised NOD-SCID-Gamma mice. RESULTS Stem cell models of the homozygous variant RFX6-/- predictably failed to generate insulin-secreting pancreatic beta cells, mirroring the phenotype observed in Mitchell-Riley syndrome. Notably, at the pancreatic endocrine stage, there was an upregulation of precursor markers NEUROG3 and SOX9, accompanied by increased apoptosis. Intriguingly, heterozygous RFX6+/- SC-islets exhibited RFX6 haploinsufficiency (54.2% reduction in protein expression), associated with reduced beta cell maturation markers, altered calcium signalling and impaired insulin secretion (62% and 54% reduction in basal and high glucose conditions, respectively). However, RFX6 haploinsufficiency did not have an impact on beta cell number or insulin content. The reduced insulin secretion persisted after in vivo implantation in mice, aligning with the increased risk of variant carriers to develop diabetes. CONCLUSIONS/INTERPRETATION Our allelic series isogenic SC-islet models represent a powerful tool to elucidate specific aetiologies of diabetes in humans, enabling the sensitive detection of aberrations in both beta cell development and function. We highlight the critical role of RFX6 in augmenting and maintaining the pancreatic progenitor pool, with an endocrine roadblock and increased cell death upon its loss. We demonstrate that RFX6 haploinsufficiency does not affect beta cell number or insulin content but does impair function, predisposing heterozygous carriers of loss-of-function variants to diabetes. DATA AVAILABILITY Ultra-deep bulk RNA-seq data for pancreatic differentiation stages 3, 5 and 7 of H1 RFX6 genotypes are deposited in the Gene Expression Omnibus database with accession code GSE234289. Original western blot images are deposited at Mendeley ( https://data.mendeley.com/datasets/g75drr3mgw/2 ).
Collapse
Affiliation(s)
- Hazem Ibrahim
- Stem Cells and Metabolism Research Program, Faculty of Medicine, University of Helsinki, Helsinki, Finland.
| | - Diego Balboa
- Stem Cells and Metabolism Research Program, Faculty of Medicine, University of Helsinki, Helsinki, Finland
| | - Jonna Saarimäki-Vire
- Stem Cells and Metabolism Research Program, Faculty of Medicine, University of Helsinki, Helsinki, Finland
| | - Hossam Montaser
- Stem Cells and Metabolism Research Program, Faculty of Medicine, University of Helsinki, Helsinki, Finland
| | - Oleg Dyachok
- Department of Medical Cell Biology, Uppsala University, Uppsala, Sweden
| | - Per-Eric Lund
- Department of Medical Cell Biology, Uppsala University, Uppsala, Sweden
| | | | - Jouni Kvist
- Stem Cells and Metabolism Research Program, Faculty of Medicine, University of Helsinki, Helsinki, Finland
| | - Om P Dwivedi
- Institute for Molecular Medicine Finland, FIMM, HiLIFE, Helsinki, Finland
- Research Program of Clinical and Molecular Metabolism, University of Helsinki, Helsinki, Finland
| | - Väinö Lithovius
- Stem Cells and Metabolism Research Program, Faculty of Medicine, University of Helsinki, Helsinki, Finland
| | - Tom Barsby
- Stem Cells and Metabolism Research Program, Faculty of Medicine, University of Helsinki, Helsinki, Finland
| | - Vikash Chandra
- Stem Cells and Metabolism Research Program, Faculty of Medicine, University of Helsinki, Helsinki, Finland
| | - Solja Eurola
- Stem Cells and Metabolism Research Program, Faculty of Medicine, University of Helsinki, Helsinki, Finland
| | - Jarkko Ustinov
- Stem Cells and Metabolism Research Program, Faculty of Medicine, University of Helsinki, Helsinki, Finland
| | - Tiinamaija Tuomi
- Institute for Molecular Medicine Finland, FIMM, HiLIFE, Helsinki, Finland
- Research Program of Clinical and Molecular Metabolism, University of Helsinki, Helsinki, Finland
- Folkhälsan Institute of Genetics, Folkhälsan Research Center, Biomedicum Helsinki, Finland
- Abdominal Center, Endocrinology, University of Helsinki and Helsinki University Hospital, Helsinki, Finland
- Lund University Diabetes Centre, Department of Clinical Sciences, Lund University, Lund, Sweden
| | - Päivi J Miettinen
- Department of Pediatrics, Helsinki University Hospital, Helsinki, Finland
| | - Sebastian Barg
- Department of Medical Cell Biology, Uppsala University, Uppsala, Sweden
| | - Anders Tengholm
- Department of Medical Cell Biology, Uppsala University, Uppsala, Sweden
| | - Timo Otonkoski
- Stem Cells and Metabolism Research Program, Faculty of Medicine, University of Helsinki, Helsinki, Finland.
- Department of Pediatrics, Helsinki University Hospital, Helsinki, Finland.
| |
Collapse
|
3
|
Wicik Z, Nowak A, Jarosz-Popek J, Wolska M, Eyileten C, Siller-Matula JM, von Lewinski D, Sourij H, Filipiak‬ KJ, Postuła M. Characterization of the SGLT2 Interaction Network and Its Regulation by SGLT2 Inhibitors: A Bioinformatic Analysis. Front Pharmacol 2022; 13:901340. [PMID: 36046822 PMCID: PMC9421436 DOI: 10.3389/fphar.2022.901340] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2022] [Accepted: 06/22/2022] [Indexed: 11/29/2022] Open
Abstract
Background: Sodium–glucose cotransporter 2 (SGLT2), also known as solute carrier family 5 member 2 (SLC5A2), is a promising target for a new class of drugs primarily established as kidney-targeting, effective glucose-lowering agents used in diabetes mellitus (DM) patients. Increasing evidence indicates that besides renal effects, SGLT2 inhibitors (SGLT2i) have also a systemic impact via indirectly targeting the heart and other tissues. Our hypothesis states that the pleiotropic effects of SGLT2i are associated with their binding force, location of targets in the SGLT2 networks, targets involvement in signaling pathways, and their tissue-specific expression. Methods: Thus, to investigate differences in SGLT2i impact on human organisms, we re-created the SGLT2 interaction network incorporating its inhibitors and metformin and analyzed its tissue-specific expression using publicly available datasets. We analyzed it in the context of the so-called key terms ( autophagy, oxidative stress, aging, senescence, inflammation, AMPK pathways, and mTOR pathways) which seem to be crucial to elucidating the SGLT2 role in a variety of clinical manifestations. Results: Analysis of SGLT2 and its network components’ expression confidence identified selected organs in the following order: kidney, liver, adipose tissue, blood, heart, muscle, intestine, brain, and artery according to the TISSUES database. Drug repurposing analysis of known SGLT2i pointed out the influence of SGLT1 regulators on the heart and intestine tissue. Additionally, dapagliflozin seems to also have a stronger impact on brain tissue through the regulation of SGLT3 and SLC5A11. The shortest path analysis identified interaction SIRT1-SGLT2 among the top five interactions across six from seven analyzed networks associated with the key terms. Other top first-level SGLT2 interactors associated with key terms were not only ADIPOQ, INS, GLUT4, ACE, and GLUT1 but also less recognized ILK and ADCY7. Among other interactors which appeared in multiple shortest-path analyses were GPT, COG2, and MGAM. Enrichment analysis of SGLT2 network components showed the highest overrepresentation of hypertensive disease, DM-related diseases for both levels of SGLT2 interactors. Additionally, for the extended SGLT2 network, we observed enrichment in obesity (including SGLT1), cancer-related terms, neuroactive ligand–receptor interaction, and neutrophil-mediated immunity. Conclusion: This study provides comprehensive and ranked information about the SGLT2 interaction network in the context of tissue expression and can help to predict the clinical effects of the SGLT2i.
Collapse
Affiliation(s)
- Zofia Wicik
- Center for Preclinical Research and Technology CEPT, Department of Experimental and Clinical Pharmacology, Medical University of Warsaw, Warsaw, Poland
| | - Anna Nowak
- Center for Preclinical Research and Technology CEPT, Department of Experimental and Clinical Pharmacology, Medical University of Warsaw, Warsaw, Poland
- Doctoral School, Medical University of Warsaw, Warsaw, Poland
| | - Joanna Jarosz-Popek
- Center for Preclinical Research and Technology CEPT, Department of Experimental and Clinical Pharmacology, Medical University of Warsaw, Warsaw, Poland
- Doctoral School, Medical University of Warsaw, Warsaw, Poland
| | - Marta Wolska
- Center for Preclinical Research and Technology CEPT, Department of Experimental and Clinical Pharmacology, Medical University of Warsaw, Warsaw, Poland
- Doctoral School, Medical University of Warsaw, Warsaw, Poland
| | - Ceren Eyileten
- Center for Preclinical Research and Technology CEPT, Department of Experimental and Clinical Pharmacology, Medical University of Warsaw, Warsaw, Poland
- Genomics Core Facility, Centre of New Technologies, University of Warsaw, Warsaw, Poland
| | - Jolanta M. Siller-Matula
- Center for Preclinical Research and Technology CEPT, Department of Experimental and Clinical Pharmacology, Medical University of Warsaw, Warsaw, Poland
- Department of Cardiology, Medical University of Vienna, Vienna, Austria
| | - Dirk von Lewinski
- Department of Internal Medicine, Division of Cardiology, Medical University of Graz, Graz, Austria
| | - Harald Sourij
- Division of Endocrinology and Diabetology, Interdisciplinary Metabolic Medicine Trials Unit, Medical University of Graz, Graz, Austria
| | | | - Marek Postuła
- Center for Preclinical Research and Technology CEPT, Department of Experimental and Clinical Pharmacology, Medical University of Warsaw, Warsaw, Poland
- *Correspondence: Marek Postuła,
| |
Collapse
|