Identification of novel genes associated with atherosclerosis in Bama miniature pig

Abstract

Background: Atherosclerosis is a chronic cardiovascular disease of great concern. However, it is difficult to establish a direct connection between conventional small animal models and clinical practice. The pig’s genome, physiology, and anatomy reflect human biology better than other laboratory animals, which is crucial for studying the pathogenesis of atherosclerosis.

Methods: We used whole-genome sequencing data from nine Bama minipigs to perform a genome-wide linkage analysis, and further used bioinformatic tools to filter and identify underlying candidate genes. Candidate gene function prediction was performed using the online prediction tool STRING 12.0. Immunohistochemistry and immunofluorescence were used to detect the expression of proteins encoded by candidate genes.

Results: We mapped differential single nucleotide polymorphisms (SNPs) to genes and obtained a total of 102 differential genes, then we used GO and KEGG pathway enrichment analysis to identify four candidate genes, including SLA-1, SLA-2, SLA-3, and TAP2. nsSNPs cause changes in the primary and tertiary structures of SLA-I and TAP2 proteins, the primary structures of these two proteins have undergone amino acid changes, and the tertiary structures also show slight changes. In addition, immunohistochemistry and immunofluorescence results showed that the expression changes of TAP2 protein in coronary arteries showed a trend of increasing from the middle layer to the inner layer.

Conclusions: We have identified SLA-I and TAP2 as potential susceptibility genes of atherosclerosis, highlighting the importance of antigen processing and immune response in atherogenesis.

Keywords

atherosclerosis / candidate genes / genome-wide linkage analysis / major histocompatibility complex / whole genome sequencing

References

[1] Kobiyama K, Ley K. Atherosclerosis. Circ Res. 2018;123:1118-1120.

[2] Soliman GA. Dietary fiber, atherosclerosis, and cardiovascular disease. Nutrients. 2019;11:1155.

[3] Chen Z, Schunkert H. Genetics of coronary artery disease in the post-GWAS era. J Intern Med. 2021;290:980-992.

[4] Nikpay M, Goel A, Won H-H. et al. A comprehensive 1, 000 genomes-based genome-wide association meta-analysis of coronary artery disease. Nat Genet. 2015;47:1121-1130.

[5] Schunkert H, König IR, Kathiresan S, et al. Large-scale association analysis identifies 13 new susceptibility loci for coronary artery disease. Nat Genet. 2011;43:333-338.

[6] Ott J, Wang J, Leal SM. Genetic linkage analysis in the age of whole-genome sequencing. Nat Rev Genet. 2015;16:275-284.

[7] An T, Zhang J, Liu YF, et al. Combined analysis of whole-exon sequencing and lncRNA sequencing in type 2 diabetes mellitus patients with obesity. J Cell Mol Med. 2020;24:2451-2463.

[8] van der Harst P, Verweij N. Identification of 64 novel genetic loci provides an expanded view on the genetic architecture of coronary artery disease. Circ Res. 2018;122:433-443.

[9] Erdmann J, Kessler T, Munoz Venegas L, Schunkert H. A decade of genome-wide association studies for coronary artery disease: the challenges ahead. Cardiovasc Res. 2018;114:1241-1257.

[10] Schunkert H. Family or SNPs: what counts for hereditary risk of coronary artery disease? Eur Heart J. 2016;37:568-571.

[11] Getz GS, Reardon CA. Animal models of atherosclerosis. Arterioscler Thromb Vasc Biol. 2012;32:1104-1115.

[12] Zhang Y, Fatima M, Hou S, Bai L, Zhao S, Liu E. Research methods for animal models of atherosclerosis (review). Mol Med Rep. 2021;24:871.

[13] Horodinschi R-N, Stanescu AMA, Bratu OG, et al. Treatment with statins in elderly patients. Medicina (Kaunas). 2019;55:721.

[14] Bentzon JF, Falk E. Atherosclerotic lesions in mouse and man: is it the same disease? Curr Opin Lipidol. 2010;21:434-440.

[15] Zhao Y, Xiang L, Liu Y, Niu M, Yuan J, Chen H. Atherosclerosis induced by a high-cholesterol and high-fat diet in the inbred strain of the Wuzhishan miniature pig. Anim Biotechnol. 2018;29:110-118.

[16] Hamamdzic D, Wilensky RL. Porcine models of accelerated coronary atherosclerosis: role of diabetes mellitus and hypercholesterolemia. J Diabetes Res. 2013;2013:761415.

[17] Skold BH, Getty R, Ramsey FK. Spontaneous atherosclerosis in the arterial system of aging swine. Am J Vet Res. 1966;27:257-273.

[18] Stary HC, Chandler AB, Glagov S, et al. A definition of initial, fatty streak, and intermediate lesions of atherosclerosis. A report from the committee on vascular lesions of the council on arteriosclerosis, American Heart Association. Circulation. 1994;89:2462-2478.

[19] Stary HC, Chandler AB, Glagov S, et al. A definition of initial, fatty streak, and intermediate lesions of atherosclerosis. A report from the committee on vascular lesions of the council on arteriosclerosis, American Heart Association. Arterioscler Thromb. 1994;14:840-856.

[20] Kawaguchi H, Yamada T, Miura N, et al. Rapid development of atherosclerosis in the world’s smallest microminipig fed a high-fat/high-cholesterol diet. J Atheroscler Thromb. 2014;21:186-203.

[21] Zhao Y, Niu M, Jia Y, et al. Establishment of type 2 diabetes mellitus models using streptozotocin after 3 months high-fat diet in Bama minipigs. Anim Biotechnol. 2023;34:2295-2312.

[22] Kovacic S, Bakran M. Genetic susceptibility to atherosclerosis. Stroke Res Treat. 2012;2012:362941.

[23] Mirhadi S, Tam S, Li Q, et al. Integrative analysis of non-small cell lung cancer patient-derived xenografts identifies distinct proteotypes associated with patient outcomes. Nat Commun. 2022;13:1811.

[24] Li H, Durbin R. Fast and accurate short read alignment with burrows-wheeler transform. Bioinformatics. 2009;25:1754-1760.

[25] Li H, Handsaker B, Wysoker A, et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25:2078-2079.

[26] uinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26:841-842.

[27] McKenna A, Hanna M, Banks E, et al. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20:1297-1303.

[28] Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010;38:e164.

[29] Yu G, Wang L-G, Han Y, He Y. clusterProfiler: an R package for comparing biological themes among gene clusters. Omics. 2012;16:284-287.

[30] Jumper J, Evans R, Pritzel A, et al. Highly accurate protein structure prediction with AlphaFold. Nature. 2021;596:583-589.

[31] Sonzogni O, Haynes J, Seifried LA, et al. Reporters to mark and eliminate basal or luminal epithelial cells in culture and in vivo. PLoS Biol. 2018;16:e2004049.

[32] Nunes Dos Santos K, Florentino RM, França A, et al. Polymorphism in the promoter region of NFE2L2 gene is a genetic marker of susceptibility to cirrhosis associated with alcohol abuse. Int J Mol Sci. 2019;20:3859.

[33] Wang Y, Zhao J, Deng X, et al. Genome-wide analysis of serine carboxypeptidase-like protein (SCPL) family and functional validation of Gh_SCPL42 unchromosome conferring cotton verticillium der verticillium wilt stress in Gossypium hirsutum. BMC Plant Biol. 2022;22:421.

[34] Pan X, Zhang N, Wei X, et al. Illumination of PRRSV cytotoxic T lymphocyte epitopes by the three-dimensional structure and peptidome of swine lymphocyte antigen class I (SLA-I). Front Immunol. 2019;10:2995.

[35] Liu X, Yin L, Xue M, et al. Coronavirus porcine Deltacoronavirus upregulates MHC class I expression through RIG-I/IRF1-mediated NLRC5 induction. J Virol. 2022;96:e0015822.

[36] Fan J, Watanabe T. Atherosclerosis: known and unknown. Pathol Int. 2022;72:151-160.

[37] Hemmat N, Ebadi A, Badalzadeh R, Memar MY, Baghi HB. Viral infection and atherosclerosis. Eur J Clin Microbiol Infect Dis. 2018;37:2225-2233.

[38] Joly JH, Lowry WE, Graham NA. Differential gene set enrichment analysis: a statistical approach to quantify the relative enrichment of two gene sets. Bioinformatics. 2021;36:5247-5254.

[39] Burger F, Baptista D, Roth A, et al. Single-cell RNA-seq reveals a crosstalk between hyaluronan receptor LYVE-1-expressing macrophages and vascular smooth muscle cells. Cells. 2022;11:411.

[40] Randolph TC, Broderick S, Shaw LK, et al. Race and sex differences in RS interval and associated outcome among patients with left ventricular systolic dysfunction. J Am Heart Assoc. 2017;6:e004381.

[41] Newby PR, Pickles OJ, Mazumdar S, et al. Follow-up of potential novel Graves’disease susceptibility loci, identified in the UK WTCCC genome-wide nonsynonymous SNP study. Eur J Hum Genet. 2010;18:1021-1026.

[42] Fan S, Wu Y, Wang S, et al. Structural and biochemical analyses of swine major histocompatibility complex class I complexes and prediction of the epitope map of important influenza a virus strains. J Virol. 2016;90:6625-6641.

[43] Zhang N, i J, Feng S, et al. Crystal structure of swine major histocompatibility complex class I SLA-1 0401 and identification of 2009 pandemic swine-origin influenza A H1N1 virus cytotoxic T lymphocyte epitope peptides. J Virol. 2011;85:11709-11724.

[44] Lunney JK, Ho C-S, Wysocki M, Smith DM. Molecular genetics of the swine major histocompatibility complex, the SLA complex. Dev Comp Immunol. 2009;33:362-374.

[45] Grainger AT, Jones MB, Li J, Chen MH, Manichaikul A, Shi W. Genetic analysis of atherosclerosis identifies a major susceptibility locus in the major histocompatibility complex of mice. Atherosclerosis. 2016;254:124-132.

[46] Kobiyama K, Vassallo M, Mitzi J, et al. A clinically applicable adjuvant for an atherosclerosis vaccine in mice. Eur J Immunol. 2018;48:1580-1587.

[47] Abele R, Tampé R. The ABCs of immunology: structure and function of TAP, the transporter associated with antigen processing. Physiology (Bethesda). 2004;19:216-224.

[48] Ozbas-Gerceker F, Bozman N, Gezici S, et al. Association of TAP1 and TAP2 gene polymorphisms with hematological malignancies. Asian Pac J Cancer Prev. 2013;14:5213-5217.

[49] Moins-Teisserenc H, Semana G, Alizadeh M, et al. TAP2 gene polymorphism contributes to genetic susceptibility to multiple sclerosis. Hum Immunol. 1995;42:195-202.

[50] Roh EY, Yoon JH, Shin S, Song EY, Park MH. Association of TAP1 and TAP2 genes with susceptibility to pulmonary tuberculosis in Koreans. APMIS. 2015;123:457-464.

[51] Cresswell P, Ackerman AL, Giodini A, Peaper DR, Wearsch PA. Mechanisms of MHC class I-restricted antigen processing and cross-presentation. Immunol Rev. 2005;207:145-157.

[52] Dimayuga PC, Zhao X, Yano J, et al. Identification of apoB-100 peptide-specific CD8+T cells in atherosclerosis. J Am Heart Assoc. 2017;6:6.

[53] Seijkens TTP, Poels K, Meiler S, et al. Deficiency of the T cell regulator casitas B-cell lymphoma-B aggravates atherosclerosis by inducing CD8+T cell-mediated macrophage death. Eur Heart J. 2019;40:372-382.

[54] Huang L, Zheng Y, Yuan X, et al. Decreased frequencies and impaired functions of the CD31+subpopulation in Treg cells associated with decreased FoxP3 expression and enhanced Treg cell defects in patients with coronary heart disease. Clin Exp Immunol. 2017;187:441-454.

[55] Jukema RA, Ahmed TAN, Tardif J-C. Does low-density lipoprotein cholesterol induce inflammation? If so, does it matter? Current insights and future perspectives for novel therapies. BMC Med. 2019;17:197.

[56] Vitale M, Rezzani R, Rodella L, et al. HLA class I antigen and transporter associated with antigen processing (TAP1 and TAP2) down-regulation in high-grade primary breast carcinoma lesions. Cancer Res. 1998;58:737-742.

RIGHTS & PERMISSIONS

2024 The Author(s). Animal Models and Experimental Medicine published by John Wiley & Sons Australia, Ltd on behalf of The Chinese Association for Laboratory Animal Sciences.