Comparative analysis of immune related genes between domestic pig and germ-free minipig
Laboratory Animal Research volume 36, Article number: 44 (2020)
Recently, minipig has been considered as an animal model that is appropriate for human disease model to study toxicology, pharmacology, and xenotransplantation. Nevertheless, minipigs are bred in various environment according to their use. Here, we suggest that minipigs used for research should be bred in well-controlled facility, comparing immune status of pigs raised in different breeding environment. DNA microarray was performed with ear skin and placenta of Landrace domestic pigs (DPs) and Minnesota germ-free minipigs (GPs). Their immune transcriptome was analyzed by gene ontology (GO) annotation database, based on criteria of |log2 fold change| ≥1 with P ≤ 0.05. As a result, we found that immune related genes in the ear skin of DPs were highly activated, compared to GPs. On the other hand, no significant s were found in the placenta. Quantitative real-time PCR (qRT-PCR) was performed in five candidate immune genes. Their fold changes were consistent with the results from DNA microarray (P ≤ 0.05). In conclusion, we experimentally proved that porcine immune system was affected by different breeding environment, suggesting the importance of controlling microbes in animal room for the qualified research.
Experimental animal models have been considered as important research tool to conduct preclinical study and identify human disease mechanism. The species of animal model have been also diversified, such as rodents, non-rodents, and non-human primates . In general toxicity test, one species of non-rodents is essentially required to meet Good Laboratory Practice (GLP) guideline . However, arguments to utilize those animals for research have been constantly raised due to ethical issues and lacks of experiment reproducibility . To overcome these problems, the minipig has been raised as an alternative animal due to similarity of pathological features with human, and being relatively free from animal ethic issues. The advantages of utilizing minipig are as follows . First, the minipig is relatively free from animal ethic issue than other animals because of their industrial use as the livestock in human society. Second, the minipig has been improved for proper size to be managed and bred easier than traditional farm pigs, which enables to reduce economic costs incurred by securing breeding space for installment of cages and their feed consumption. In addition, these features produce reliable data for researcher by increasing the number of animals in experiment. Finally, the minipig shares similar histological and physiological characteristics with human, which makes them as standard resource for xenotransplantation. Despite these advantages, it is highly necessary to breed minipigs in the well-controlled facility such as Specific Pathogen Free or Germ-free grade . Immunogen, such as bacteria, in the conventional facility could bring out infection and induce activation of immune response, which leads to affect experiment data. Moreover, a recent report warns of transmission of zoonotic disease from porcine organ to human patient through xenotransplantation .
Thus, it should be experimentally identified that immunological change of pigs in different breeding environment for researchers to perform safe and qualified experiments. We chose ear skin, primary immunological barrier against external environment, as target tissue . Placenta known for immune privilege region was chosen as control tissue to confirm endogenous variation caused by immunological feature of tissues . Transcriptome profiling and comparative analysis between domestic pigs (DPs) and germ-free minipigs (GPs) were performed using DNA microarray and quantitative real-time PCR (qRT-PCR).
A pregnant Minnesota minipig raised in sterile condition (Laboratory Animal Research Center in Konkuk University, Seoul, Korea) was randomly selected from pregnant sows. The facility continuously maintained positive pressure with HEPA-filtered air, and temperature (22 ± 2 °C) and humidity (50 ± 5%). Newborn germ-free piglets (3 males and 1 females; n = 4) were produced by hysterectomy of pregnant sow under aseptic conditions. Average body weights of the piglets ranged from 450 to 550 g. They were transferred to the aseptic isolator (1200 W × 900D × 950H m/m, SK-ISO1700HBP600; Three-shine Inc., Daejeon, Korea), and fed sterilized soy bean milk (Fig. 1a). The procedures that contained animal welfare and ethical problem were approved by the Konkuk University Institutional Animal Care and Use Committee (IACUC; approval number: KU01410). A pregnant Landrace pig raised in conventional condition (Dongsan Farm, Younchun, Gyunggido, Korea) was also selected from pregnant sows. The average body weights of domestic piglets (3 males and 1 females; n = 4) ranged from 1.05 to 1.25 kg. All breeding process for the piglets were performed under same condition with that for Germ-free piglets except controlling microbial environment.
RNA sample preparation
Chorion parts of placenta (n = 4) were specifically collected and frozen with liquid nitrogen immediately at birth. Ear skins of germ-free piglets (n = 4) were obtained in the aseptic isolator on 3 weeks after birth. To prevent RNA degradation, frozen placenta samples are emerged in 1 ml of a RNAlater-ICE (Ambion, Invitrogen Carlsbad, CA, USA) solution and incubated in 4 °C for over-night. Total RNA was extracted from the ear skins and the placentas by using Phenol/Chloroform extraction method. All RNA samples were pooled for DNA microarray. Reverse transcription reaction for qRT-PCR was performed with a QuantiNova Reverse Transcription Kit (QIAGEN Science Inc., Germantown, MD, USA).
DNA microarray and gene ontology analysis
DNA microarray was performed using Porcine (v2) Gene expression 4 × 44 K microarray (Cat.no G2519F-025440; Agilent technology, Santa Clara, CA, USA) covered with 43,803 probes in capable of detecting every porcine gene. The procedure was followed by technical manual of SurePrint Gene Expression Array (Agilent technology). Quality check of RNA samples was performed using Agilent 2100 Bioanlayzer (Agilent technology). Gene Ontology (GO) was analyzed by using DAVID Bioinformatics Resources 6.87 .
Quantitative real-time PCR
qRT-PCR was performed in 20 μl PCR reaction mixtures that included 2x iQ SYBR Green Supermix (Bio-rad Laboratories. Hercules, CA, USA, each primer at a concentration of 0.5 μM, and template cDNA according to the manufacturer’s protocol. it was run in a QuantStudio 5 Real-Time PCR System (Applied Biosystmes, Foster city, CA, USA) using the amplification parameters: the initial denaturation step at 95 °C for 3 min was followed by 40 cycles of denaturation at 95 °C for 20 s, annealing and elongation at 60 °C for 1 min without the final elongation step. Gene expression level of Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) was used for an endogenous control to normalize each RNA samples. All primer sequences are shown in Table 1 [10,11,12,13,14,15]. 2-ΔΔCt method was used for the relative quantification of target genes.
All samples for qRT-PCR was run in triplicate and result graph was expressed as means ± S.D. The P value was determined using Student’s t test. Pearson’s correlation coefficient was used to measure statistical relationship between qRT-PCR and DNA microarray.
Summary of DNA microarray data
A total of 13,673 and 3599 differentially expressed genes (DEGs) were found in each ear skin and placenta (P < 0.05, |log2 fold| ≥ 1). Of total DEGs, we classified immune-related genes to confirm immunological difference between DPs and GPs. As a result, 90 and 44 genes were found in each ear skin and placenta (Fig. 1b). Specific information was provided in Table S1.
Comparison of immune related genes in the ear skin and placenta of DP compared to GP
First, we found that the number of up-regulated genes in the ear skin (67/73) was about 7 times greater than that in the placenta (10/73), compared DPs with GPs. However, no significant differences were found in comparison result of down-regulated genes between ear skin (23/53) and placenta (34/53) (Fig. 1c). In addition, we also sorted out 20 immune related genes that were commonly expressed in both tissues for accurate comparison. As a result, the number of up-regulated genes in ear skin (14/20) about 3 times greater than that in placenta (4/20), showing adverse pattern between two tissues (Fig. 1d). Thus, it could be expected that the DPs was infected with pathogens under conventional condition, resulting in immune activation.
Systemic immune activation in DP
Furthermore, we performed clustering analysis using DNA microarray data of the ear skin to confirm systemic immune activation. Immune related genes were grouped by biological functions that associated with a process of systemic immune activation. Thus, 22 genes were divided into 4 groups (Complement, Innate immunity, T cell signaling, B cell signaling; (Fig. 1e). This result indicated that immune activation was systemically occurred in the body of DPs.
Validation of DEG data by using quantitative real-time PCR
qRT-PCR was performed with five candidate genes including Interleukin 18 receptor alpha (IL-18R1), Cluster of differentiation 40 ligand (CD40LG) for ear skin, Interleukin 1 beta (IL-1B) for placenta, Interleukin 6 (IL-6), and 2′-5′-Oligoadenylate synthetase 1 (OAS1) to validate DNA microarray result. As a result, the fold changes of all candidate genes were highly consistent with DNA microarray results (Fig. 2a). Moreover, we calculate correlation coefficient to measure significance of a linear relationship between qRT-PCR and DNA microarray. In conclusion, expression pattern of candidate genes in the ear skin and placenta were highly correlated (Fig. 2b). Specific gene expression result was provided in Table S2.
Our analysis results demonstrated that porcine immune system was affected by microbes in breeding environment. The number of immune-related genes in the ear skin was greater than that in the placenta, which meant that intense immunological changes occur in the marginal tissue (Fig. 1b). We also could verify different expression patterns in placenta, implicating that unidentified mechanism suppress aggressive immune response to prevent its fetus from being killed [16, 17]. According to recent study, it was identified that specific mechanism of immune tolerance between maternal endothelium and fetal tissues during pregnancy .
Contrary to expectations, the expression patterns of some genes were not consistent with our conclusion. These results might be explained by the fact that differences of expression value were too small to be statistically significant (Ras-related protein Rab-27A (RAB27A), Interleukin 1 beta (IL-1B), KIT proto-oncogene, receptor tyrosine kinase (KIT), and MHC class II histocompatibility antigen SLA-DQA (SLA-DQA1)). In addition, homologous expression pattern of Colony stimulating factor 2 (CSF2), Vitronectin (VTN), and OAS1 in both ear skin and placenta or increased expression of immune genes in the placenta of germ-free minipigs have not been clearly explained. However, some studies reported that expression of some immune genes in placenta could be unexpected because of variable reasons, such as epigenetic modification and its unique microenvironment [18, 19]. In addition, pregnant mice raised in conventional and germ-free environment represented different immunological adaptations . These studies suggested that un identified immune modulation were occurred during pregnancy.
In spite of some experimental limitations, we identified that uncontrolled breeding facility could spoil immune system of experimental animals. Clinically, porcine skin have been well known for good model in toxicology studies due to its similarity of histological structure and skin immune system with human . This study could support the idea that researchers should breed minipigs in the bacteria-controlled facility to perform qualified animal experiment.
Availability of data and materials
All experiment data during this study are included in this manuscript and supporting files.
Germ free minipig
Good laboratory practice
Quantitative real-time PCR
Differentially expressed gene
Interleukin 18 receptor alpha
Cluster of differentiation 40 ligand
2′-5-Oligoadenylate synthetase 1
Ras-related protein Rab-27A
KIT proto-oncogene, receptor tyrosine kinase
MHC class II histocompatibility antigen SLA-DQA
Colony stimulating factor 2
Interleukin 1 beta
Ganderup NC, Harvey W, Mortensen JT, Harrouk W. The minipig as nonrodent species in toxicology--where are we now? Int J Toxicol. 2012;31(6):507–28.
Helen P, Paul B, de Lolke H, Downes N, Jones K, Cassen EM, Kimber I. Reviewing the utility of two species in general toxicology related to drug development. Int J Toxicol. 2018;37(2):121–4.
Hartung T. Thoughts on limitations of animal models. Parkinsonism Related Disord. 2008;14(Suppl 2):S81–3.
Vodicka P, Smetana K, Dvoránková B, Emerick T, Xu YZ, Ourednik J, Ourednik V, Motlík J. The miniature pig as an animal model in biomedical research. Ann N Y Acad Sci. 2005;1049:161–71.
Safron J, Gonder JC. The SPF pig in research. ILAR J. 1997;38(1):28–31.
Niu D, Wei HJ, Lin L, Geroge H, Wang T, Lee IH, Zhao HY, Wang Y, Kan Y, Shrock E, Lesha E, Wang G, Luo Y, Qing Y, Jiao D, Zhao H, Zhhou X, Wnag S, Wei H, Güell M, Church GM, Yang L. Inactivation of porcine endogenous retrovirus in pigs using CRISPR-Cas9. Science. 2017;357(6357):1303–7.
Prescott SL, Larcombe DL, Logan AC, West C, Burks W, Caraballo L, Levin M, Etten VE, Horwitz P, Kozyrskyj A, Campbell DE. The skin microbiome: impact of modern environments on skin ecology, barrier integrity, and systemic immune programming. World Allergy Organ J. 2017;10(1):29.
Gobert M, Lafaille JJ. Maternal-fetal immune tolerance, block by block. Cell. 2012;150(1):7–9.
Huang DW, Sherman BT, Lempicki RA. Bioinformatics enrichment tools:paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37(1):1–13.
de Goffau MC, Lager S, Sovio U, Gaccioli F, Cook E, Peacock SJ, Parkhill J, Charnock-Jones DS, Smith GCS. Human placenta has no microbiome but can contain potential pathogens. Nature. 2019;572(7769):329–34.
Geraldes P, Gagnon S, Hadjadj S, Merhi Y, Sirois MG, Cloutier I, Tanguay JF. Estradiol blocks the induction of CD40 and CD40L expression on endothelial cells and prevents neutrophil adhesion: an ERalpha-mediated pathway. Cardiovasc Res. 2006;71(3):566–73.
Meijer RI, Hoevenaars FPM, Serné EH, Yudkin JS, Kokhuis TJA, Weijers EM, Hinsbergh VWMV, Smulders YM, Eringa EC. JNK2 in myeloid cells impairs insulin's vasodilator effects in muscle during early obesity development through perivascular adipose tissue dysfunction. Am J Physiol Heart Circ Physiol. 2019;317(2):H364–74.
Nagai Y, Watanabe K, Aso H, Ohwada S, Muneta Y, Yamaguchi T. Cellular localization of IL-18 and IL-18 receptor in pig anterior pituitary gland. Domest Anim Endocrinol. 2006;30(2):144–54.
Whyte JJ, Meyer AE, Spate LD, Benne JA, Cecil R, Samuel MS, Murphy CN, Prather RS, Geisert RD. Inactivation of porcine interleukin-1beta results in failure of rapid conceptus elongation. Proc Natl Acad Sci U S A. 2018;115(2):307–12.
Zheng S, Zhu D, Lian X, Liu W, Cao R, Chen P. Porcine 2′, 5′-oligoadenylate synthetases inhibit Japanese encephalitis virus replication in vitro. J Med Virol. 2016;88(5):760–8.
Svensson-Arvelund J, Ernerudh J, Buse E, Cline JM, Haeger JDDD, Markert UR, Pfarrer C, Vos PD, Faas MM. The placenta in toxicology. Part II: systemic and local immune adaptations in pregnancy. Toxicol Pathol. 2014;42(2):327–38.
Alijotas-Reig J, Llurba E, Gris JM. Potentiating maternal immune tolerance in pregnancy: a new challenging role for regulatory T cells. Placenta. 2014;35(4):241–8.
Tomlinson MS, Bommarito PA, Martin EM, Smeester L, Fichorova RN, Onderdonk AB, Kuban KCK, O’Shea TM, Fry RC. Microorganisms in the human placenta are associated with altered CpG methylation of immune and inflammation-related genes. PLoS One. 2017;12(12):e0188664.
Parnell LA, Briggs CM, Cao B, Delannoy-Bruno O, Schrieffer AE, Mysorekar IU. Microbial communities in placentas from term normal pregnancy exhibit spatially variable profiles. Sci Rep. 2017;7:11200.
Faas M, Liu Y, Borghuis T, Loo-Bouwman C, Harmsen H, de Vos P. Microbiota induced changes in the immune response in pregnant mice. Front immunol. 2020;10:2976.
Summerfield A, Meurens F, Ricklin ME. The immunology of the porcine skin and its value as a model for human skin. Mol Immunol. 2015;66(1):14–21.
This work was supported by the Korea Institute of Toxicology (KIT, Korea) grant funded by the Ministry of science and ICT (MIST, Korea). [Project number: KK-1911].
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Lee, J.Y., Kim, S.E., Lee, H.T. et al. Comparative analysis of immune related genes between domestic pig and germ-free minipig. Lab Anim Res 36, 44 (2020). https://doi.org/10.1186/s42826-020-00077-7
- Germ-free facility
- Differentially expressed gene
- Systemic immune activation