Expression Genome-wide Association Study Identifies Key Regulatory Variants Enriched with Metabolic and Immune Functions in Four Porcine Tissues

Overview

Journal BMC Genomics

Publisher Biomed Central

Specialty Genetics

Date 2024 Jul 11

PMID 38992576

Authors

Samin Farhangi

Marta Godia

Martijn F L Derks

Barbara Harlizius

Bert Dibbits

Rayner Gonzalez-Prendes

Richard P M A Crooijmans

Ole Madsen

Martien A M Groenen

Affiliations

Soon will be listed here.

Abstract

Background: Integration of high throughput DNA genotyping and RNA-sequencing data enables the discovery of genomic regions that regulate gene expression, known as expression quantitative trait loci (eQTL). In pigs, efforts to date have been mainly focused on purebred lines for traits with commercial relevance as such growth and meat quality. However, little is known on genetic variants and mechanisms associated with the robustness of an animal, thus its overall health status. Here, the liver, lung, spleen, and muscle transcriptomes of 100 three-way crossbred female finishers were studied, with the aim of identifying novel eQTL regulatory regions and transcription factors (TFs) associated with regulation of porcine metabolism and health-related traits.

Results: An expression genome-wide association study with 535,896 genotypes and the expression of 12,680 genes in liver, 13,310 genes in lung, 12,650 genes in spleen, and 12,595 genes in muscle resulted in 4,293, 10,630, 4,533, and 6,871 eQTL regions for each of these tissues, respectively. Although only a small fraction of the eQTLs were annotated as cis-eQTLs, these presented a higher number of polymorphisms per region and significantly stronger associations with their target gene compared to trans-eQTLs. Between 20 and 115 eQTL hotspots were identified across the four tissues. Interestingly, these were all enriched for immune-related biological processes. In spleen, two TFs were identified: ERF and ZNF45, with key roles in regulation of gene expression.

Conclusions: This study provides a comprehensive analysis with more than 26,000 eQTL regions identified that are now publicly available. The genomic regions and their variants were mostly associated with tissue-specific regulatory roles. However, some shared regions provide new insights into the complex regulation of genes and their interactions that are involved with important traits related to metabolism and immunity.

References

Goddard M, Hayes B . Mapping genes for complex traits in domestic animals and their use in breeding programmes. Nat Rev Genet. 2009; 10(6):381-91. DOI: 10.1038/nrg2575. View

Kern C, Wang Y, Xu X, Pan Z, Halstead M, Chanthavixay G . Functional annotations of three domestic animal genomes provide vital resources for comparative and agricultural research. Nat Commun. 2021; 12(1):1821. PMC: 7988148. DOI: 10.1038/s41467-021-22100-8. View

Ernst C, Steibel J . Molecular advances in QTL discovery and application in pig breeding. Trends Genet. 2013; 29(4):215-24. DOI: 10.1016/j.tig.2013.02.002. View

Ma J, Yang J, Zhou L, Ren J, Liu X, Zhang H . A splice mutation in the PHKG1 gene causes high glycogen content and low meat quality in pig skeletal muscle. PLoS Genet. 2014; 10(10):e1004710. PMC: 4207639. DOI: 10.1371/journal.pgen.1004710. View

Ponsuksili S, Murani E, Trakooljul N, Schwerin M, Wimmers K . Discovery of candidate genes for muscle traits based on GWAS supported by eQTL-analysis. Int J Biol Sci. 2014; 10(3):327-37. PMC: 3957088. DOI: 10.7150/ijbs.8134. View

Puig-Oliveras A, Revilla M, Castello A, Fernandez A, Folch J, Ballester M . Expression-based GWAS identifies variants, gene interactions and key regulators affecting intramuscular fatty acid content and composition in porcine meat. Sci Rep. 2016; 6:31803. PMC: 4989154. DOI: 10.1038/srep31803. View

Ibragimov E, Pedersen A, Xiao L, Cirera S, Fredholm M, Karlskov-Mortensen P . Analysis of merged transcriptomic and genomic datasets to identify genes and pathways underlying residual feed intake in growing pigs. Sci Rep. 2022; 12(1):21946. PMC: 9763391. DOI: 10.1038/s41598-022-26496-1. View

Ballester M, Ramayo-Caldas Y, Revilla M, Corominas J, Castello A, Estelle J . Integration of liver gene co-expression networks and eGWAs analyses highlighted candidate regulators implicated in lipid metabolism in pigs. Sci Rep. 2017; 7:46539. PMC: 5396199. DOI: 10.1038/srep46539. View

Ponsuksili S, Murani E, Brand B, Schwerin M, Wimmers K . Integrating expression profiling and whole-genome association for dissection of fat traits in a porcine model. J Lipid Res. 2011; 52(4):668-78. PMC: 3284160. DOI: 10.1194/jlr.M013342. View

10.

Revilla M, Puig-Oliveras A, Crespo-Piazuelo D, Criado-Mesas L, Castello A, Fernandez A . Expression analysis of candidate genes for fatty acid composition in adipose tissue and identification of regulatory regions. Sci Rep. 2018; 8(1):2045. PMC: 5794915. DOI: 10.1038/s41598-018-20473-3. View

11.

Clark E, Archibald A, Daetwyler H, Groenen M, Harrison P, Houston R . From FAANG to fork: application of highly annotated genomes to improve farmed animal production. Genome Biol. 2020; 21(1):285. PMC: 7686664. DOI: 10.1186/s13059-020-02197-8. View

12.

Teng J, Gao Y, Yin H, Bai Z, Liu S, Zeng H . A compendium of genetic regulatory effects across pig tissues. Nat Genet. 2024; 56(1):112-123. PMC: 10786720. DOI: 10.1038/s41588-023-01585-7. View

13.

Luo L, Van Der Zande L, van Marwijk M, Knol E, Rodenburg T, Bolhuis J . Impact of Enrichment and Repeated Mixing on Resilience in Pigs. Front Vet Sci. 2022; 9:829060. PMC: 8988148. DOI: 10.3389/fvets.2022.829060. View

14.

Warr A, Affara N, Aken B, Beiki H, Bickhart D, Billis K . An improved pig reference genome sequence to enable pig genetics and genomics research. Gigascience. 2020; 9(6). PMC: 7448572. DOI: 10.1093/gigascience/giaa051. View

15.

Love M, Huber W, Anders S . Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014; 15(12):550. PMC: 4302049. DOI: 10.1186/s13059-014-0550-8. View

16.

Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira M, Bender D . PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007; 81(3):559-75. PMC: 1950838. DOI: 10.1086/519795. View

17.

Derks M, Gross C, Lopes M, Reinders M, Bosse M, Gjuvsland A . Accelerated discovery of functional genomic variation in pigs. Genomics. 2021; 113(4):2229-2239. DOI: 10.1016/j.ygeno.2021.05.017. View

18.

Yang J, Lee S, Goddard M, Visscher P . GCTA: a tool for genome-wide complex trait analysis. Am J Hum Genet. 2010; 88(1):76-82. PMC: 3014363. DOI: 10.1016/j.ajhg.2010.11.011. View

19.

Pan Z, Yao Y, Yin H, Cai Z, Wang Y, Bai L . Pig genome functional annotation enhances the biological interpretation of complex traits and human disease. Nat Commun. 2021; 12(1):5848. PMC: 8494738. DOI: 10.1038/s41467-021-26153-7. View

20.

Bailey T, Johnson J, Grant C, Noble W . The MEME Suite. Nucleic Acids Res. 2015; 43(W1):W39-49. PMC: 4489269. DOI: 10.1093/nar/gkv416. View