Unveiling regulatory variants in the blood transcriptome and their association with immunity traits in pigs

Teodor Jové Juncà¹Daniel Crespo Piazuelo^1,2Carles Hernández Banqué¹Olga Gonzalez¹Lingzhao Fang³Raquel Quintanilla^1*Maria Ballester^1*

• ¹Institute of Agrifood Research and Technology (IRTA), Barcelona, Spain
• ²Cuarte SL, Zaragoza, Spain
• ³Aarhus University, Aarhus, Central Denmark Region, Denmark

Genome wide association studies (GWAS) have been widely used to investigate the association of genetic markers with complex traits in both humans and livestock species. A particular trait of interest, when studying animal robustness and general immunocompetence, is the transcriptomic profile of blood. To identify genetic variants affecting gene expression in pig blood, we performed expression GWAS (eGWAS) in 255 animals from a commercial Duroc population between 8,499,177 imputed single nucleotide polymorphisms (SNPs) and the expression levels of 14,642 genes obtained from RNA sequencing. Out of the nearly 125 associations tested, 23 million were found to be significant, grouped in 9,930 expression quantitative trait loci (eQTLs) associated to the expression levels of 6,051 genes. Over 36% of detected eQTLs mapped in close proximity to the genomic location of their associated gene and were classified as cis-eQTLs. Moreover, 430,694 variants were found to be associated with the expression of 10 or more different genes and were annotated as transcriptional hotspots. Among genes regulated by these hotspots, we identified transcription factors and co-factors regulating immune responses, such as ARNT, or co-expressed genes related to immunity (CSF3R, JAK2, SOCS3, STAT5B and UBE2D1) and associated with health traits, such as phagocytic activity or haptoglobin concentration. In addition, several of the cisregulating variants for immunity candidate genes overlapped with previously described immunity QTLs. Colocalization studies revealed putative common causal variants between the proportion of memory and helper T cells and the candidate genes CLEC12B, IGKV2D, KLRC1, KLRD1 and ZAP70. In conclusion, the associations identified in this study enable the characterization of transcriptional regulators of the pig blood transcriptome. Moreover, the colocalization between immunity QTLs and eQTLs has revealed potential causative mutations regulating immunocompetence in pigs. All these data and results contribute to shedding light on the regulatory mechanisms of blood gene expression and porcine immune regulation.