The eQTL colocalization and transcriptome-wide association study identify potentially causal genes responsible for economic traits in Simmental beef cattle

Wentao Cai; Yapeng Zhang; Tianpeng Chang; Zezhao Wang; Bo Zhu; Yan Chen; Xue Gao; Lingyang Xu; Lupei Zhang; Huijiang Gao; Jiuzhou Song; Junya Li

doi:10.1186/s40104-023-00876-7

The eQTL colocalization and transcriptome-wide association study identify potentially causal genes responsible for economic traits in Simmental beef cattle

J Anim Sci Biotechnol. 2023 May 11;14(1):78. doi: 10.1186/s40104-023-00876-7.

Authors

Wentao Cai¹, Yapeng Zhang¹, Tianpeng Chang¹, Zezhao Wang¹, Bo Zhu¹, Yan Chen¹, Xue Gao¹, Lingyang Xu¹, Lupei Zhang¹, Huijiang Gao¹, Jiuzhou Song², Junya Li³

Affiliations

¹ Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, 100193, China.
² Department of Animal and Avian Science, University of Maryland, College Park, MD, 20742, USA. songj88@umd.edu.
³ Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, 100193, China. lijunya@caas.cn.

Abstract

Background: A detailed understanding of genetic variants that affect beef merit helps maximize the efficiency of breeding for improved production merit in beef cattle. To prioritize the putative variants and genes, we ran a comprehensive genome-wide association studies (GWAS) analysis for 21 agronomic traits using imputed whole-genome variants in Simmental beef cattle. Then, we applied expression quantitative trait loci (eQTL) mapping between the genotype variants and transcriptome of three tissues (longissimus dorsi muscle, backfat, and liver) in 120 cattle.

Results: We identified 1,580 association signals for 21 beef agronomic traits using GWAS. We then illuminated 854,498 cis-eQTLs for 6,017 genes and 46,970 trans-eQTLs for 1,903 genes in three tissues and built a synergistic network by integrating transcriptomics with agronomic traits. These cis-eQTLs were preferentially close to the transcription start site and enriched in functional regulatory regions. We observed an average of 43.5% improvement in cis-eQTL discovery using multi-tissue eQTL mapping. Fine-mapping analysis revealed that 111, 192, and 194 variants were most likely to be causative to regulate gene expression in backfat, liver, and muscle, respectively. The transcriptome-wide association studies identified 722 genes significantly associated with 11 agronomic traits. Via the colocalization and Mendelian randomization analyses, we found that eQTLs of several genes were associated with the GWAS signals of agronomic traits in three tissues, which included genes, such as NADSYN1, NDUFS3, LTF and KIFC2 in liver, GRAMD1C, TMTC2 and ZNF613 in backfat, as well as TIGAR, NDUFS3 and L3HYPDH in muscle that could serve as the candidate genes for economic traits.

Conclusions: The extensive atlas of GWAS, eQTL, fine-mapping, and transcriptome-wide association studies aid in the suggestion of potentially functional variants and genes in cattle agronomic traits and will be an invaluable source for genomics and breeding in beef cattle.

Keywords: Cattle; Colocalization; GWAS; TWAS; eQTL mapping.

Abstract

Grants and funding