Association between Genetic Polymorphism in the Swine Leukocyte Antigen-DRA Gene and Piglet Diarrhea in Three Chinese Pig Breeds

Article information

Asian-Australas J Anim Sci. 2014;27(9):1228-1235
1College of Animal Science and Technology, Gansu Agricultural University, Lanzhou 730070, China
2Animal Husbandry and Veterinary Institute of Gansu Province, Pingliang 744000, China.
3Gansu Research Center for Swine Production Engineering and Technology, Lanzhou 730070, China.
*Corresponding Author: S. B. Gun. Tel: +86-931-7631804, Fax: +86-931-7631239, E-mail: gunsbao056@126.com
Received 2013 September 12; Revised 2013 December 30; Accepted 2014 March 27.

Abstract

The swine leukocyte antigen (SLA)-DRA locus is noteworthy among other SLA class II loci for its limited variation and has not been investigated in depth. This study was investigated to detect polymorphisms of four exons of SLA-DRA gene and its association with piglet diarrhea in Landrace, Large White and Duroc pigs. No polymorphisms were detected in exon 3, while 2 SNPs (c.178G>A and c.211T>C), 2 SNPs (c.3093A>C and c.3104C>T) and 5 SNPs (c.4167A>G, c.4184A>G, c.4194A>G, c.4246A>G and c.4293G>A) were detected in exon 1, exon 2 and exon 4 respectively, and 1 SNP (c.4081T>C) in intron 3. Statistical results showed that genotype had significant effect on piglet diarrhea, individuals with genotype BC had a higher diarrhea score when compared with the genotypes AA, AB, AC and CC. Futhermore, genotype AC had a higher diarrhea score than the genotype CC in exon 1 (p<0.05); diarrhea scores of genotype AA and BB were higher than those of genotypes AC and CC in exon 2 (p<0.05); individuals with genotype AA had a higher diarrhea score than individuals with genotype AB and BB in exon 4 (p<0.05). Fourteen common haplotypes were founded by haplotype constructing of all SNPs in the three exons, its association with piglet diarrhea appeared that Hap2, 5, 8, 10, and 14 may be the susceptible haplotypes and Hap9 may be the resistant haplotype to piglet diarrhea. The genetic variations identified of the SLA-DRA gene may potentially be functional mutations related to piglet diarrhea.

INTRODUCTION

The major histocompatibility complex (MHC) class II molecules are tightly linked and highly polymorphic cell surface glycoproteins, which preferentially recognize foreign antigen to T cells and subsequently, initiate specific immune responses against infectious disease (Kedrick and Kim, 2000; Patch et al., 2011). In swine, swine leukocyte antigen (SLA) class II antigens are heterodimeric proteins that are comprised of non-covalently linked α (34 kDa) and β (29 kDa) chains, with the α chain coded by the SLA-DRA and SLA-DQA genes and β chain coded by the SLA-DRB1 and SLA-DQB1 genes (Lunney et al., 2009).

It is evidenced that some SLA alleles and haplotypes are associated with some immune-related diseases such as antibody responses to sheep erythrocytes and cell-mediated immune responses to tuberculin-purified protein (Mallard et al., 1989; Shinkai et al., 2012), PRRS (Lewis et al., 2007; Molina et al., 2008), and other swine immune response to various infectious agents (Renard and Vaiman, 1989; Chen et al., 2005; Ho et al., 2010a). To investigate the role of SLA genes in immune responsiveness, it is essential to have a better understanding of SLA polymorphism.

The SLA-DRA gene, which encodes the α chain of DR antigen is closely linked to DRB. The full-length genome sequence of 4.6 kb, encodes for a polypeptide of 252 amino acid residues, which includes four exons and is considered to have limited allelic variation (Ho et al., 2009; Lunney et al., 2009). Therefore previous studies on genetic polymorphism of SLA class II genes have primarily focused on polymorphism in DQB, DRB, and DQA genes (Luetkemeier et al., 2009; Ho et al., 2010a; Essler et al., 2013), but rarely targeted the swine DRA gene for comprehensive analysis. Consequently, the extent of variation at the locus is not yet known.

Due to large scale intensive production systems, pigs are exposed to various pathogens that cause infectious diseases and pathological symptoms (Lunney et al., 2009). Diarrhea in newborns and pre-weaning piglets, which causes increased mortality rates, poor growth and additional medical costs, is one of the most serious threats in the swine industry. It is estimated that 50% of piglet mortality is due to diarrhea deaths (Morris et al., 2002). Our preliminary studies indicated that piglet diarrhea during pre-weaning was linked to the exon 2 of SLA-DRA and SLA-DQA gene polymorphism (Yang et al., 2012; Yang et al., 2013).

The objective of this study was to detect the genetic variation of the four exon regions of the SLA-DRA gene, and investigate the association with piglet diarrhea from Landrace, Large White and Duroc pigs using polymerase chain reaction-single strand conformation polymorphism (PCR-SSCP) and sequencing, to provide a theoretical foundation for further pig anti-disease breeding.

MATERIALS AND METHODS

Animals and genomic DNA extraction

A total of 216 pure strain piglets from 24 sires were selected from Linze Xinghua Pig Breeding Farm (Zhangye, Gansu, China), including Landrace (n = 72), Large White (n = 72) and Duroc (n = 72).

Experimental piglets were born in the same time period (within 20 days) and the experiment conducted under forage and feeding management conditions and disease control in accordance with the Guide for Chinese Feeding Standard of Swine and approved by the institute of Ministry of Agriculture Feed Industry Center, China Agricultural University (Beijing, China). All piglets were monitored two times a day during the entire suckling period (from birth to weaning, 0 to 28 days) and assigned a daily score based on visual analysis of symptoms traits as follows: 0: normal, solid feces; 1: slight diarrhea, soft and loose feces; 2: moderate diarrhea, semi-liquid feces; 3: severe diarrhea, liquid and unformed feces (Kelly et al., 1990).

Genomic DNA was isolated from swine ear tissues samples by standard phenol-chloroform extraction meth odology (Sambrook and Russell, 2001) and quantified by using a Nanodrop ND-1000 spectrophotometer (Ther mo Fisher Scientific, Wilmington, DE, USA).

Polymerase chain reaction amplification and polymerase chain reaction-single strand conformation polymorphism analyses

Four primer pairs (Table 1), were designed to amplify exons 1, 2, 3 and 4 of SLA-DRA gene according to the published SLA-DRA gene sequence in GenBank (accession No. AY303990). PCR analyses were performed in 25 μL reactions containing of 100 ng genomic DNA, 2.5 μL of 10×PCR buffer (containing Mg2+, 15 mM), 2.5 mM dNTP mixture, 10 pmol of each primer, and 5 U of Taq polymerase (TaKaRa, Dalian, China). The thermal profile consisted of 2 min at 94°C, followed by 35 cycles of 20 s at 94°C, 30 s at annealing temperatures (Table 1) and 30 s at 72°C, with a final extension of 10 min at 72°C.

Primer sequences, annealing temperature (Ta) and size of amplification products

For SSCP analysis, 2 μL PCR products were mixed with 7 μL of the denaturing solution (98% formamide, 10 mM ethylene diamine tetraacetic acid (EDTA), 0.0 25% bromophenol blue, and 0.025% xylene-cyanole), af ter denaturation at 98°C for 10 min, samples were coo led on ice immediately and loaded on polyacrylamide gel (170×170×1 mm) in 1×TBE under optimal conditio ns. Gels were silver-stained according to the method of Byun et al. (2009).

Cloning and sequencing

PCR products representative of different SSCP patterns were purified using DNA Fragment Quick Recover Kit (TaKaRa, Dalian, China). The purified PCR products were ligated into the pMD19-T vectors (TaKaRa, Dalian, China) according to the manufacturer’s recommendations. The mixture was incubated at 4°C overnight, and transformed into Escherichia coli DH5α strains. Recombinant clones were detected by blue/white screening, 10 to 15 insert-positive clones for each transformation were incubated overnight at 37°C, in a vibrating rotary incubator (200 rpm). Verification of targeted inserts were screened using the PCR-SSCP approach as described by Zhou et al. (2008), and only those clones for which the SSCP patterns matched those of the corresponding genomic DNA were selected and subsequently sequenced by Sangon Biotech (Shanghai, China) Co., Ltd.

Statistical analysis

Nucleotide and amino acid sequences were analyzed using MEGA5 (Tamura et al., 2011) and new allele sequences were submitted to GenBank database (http://www.ncbi.nlm.nih.gov/nuccore/); gene heterozygosity (He) was calculated using POPGENE software (version 3.2); polymorphism information content (PIC) was calculated by Botstein’s methods (1980). Haplotype construction and linkage disequilibrium were analyzed by PHASE 2.0 (Stephens and Donnelly, 2003) and DNAsp 5.0 (Hudson, 2001) software respectively.

The effects of genotype on piglet diarrhea were analyzed using the general linear model procedure of SAS software (SAS, 1999) according to the following statistical model:

Yijklm=μ+Bi+Sj+Gk+eijkl

Where Yijkl is the observation value of piglet diarrhea score, μ is the overall population means, Bi is the fixed effect of breed, Sj is fixed effect of gender, Gk is the effect of genotype, and e is the random residual effect. Preliminary analysis also included the fixed interaction effects of breed, sex and genotype, however, these interaction effects were subsequently removed because they did not have a significant effect.

RESULTS

Genetic polymorphism of the SLA-DRA gene

There was no polymorphism in exon 3 of the DRA locus. Six genotypes were detected in exon 1 and exon 2 (Figure 1A and Figure 1B), respectively. They were defined as AA, AB, BB, AC, CC, and BC. For exon 1, genotypes AA, BC, AC, and CC were detected in all three pig populations, genotype BB was only detected in Duroc. Three alleles were described as: A (A-T), B (G-C), and C (G-T) in this locus (Table 2). Allele C was the most common in Large White with a frequency of 0.55, whereas allele A was predominate in Landrace and Duroc with frequencies of 0.51 and 0.44, respectively.

Figure 1

PCR-SSCP analysis of PCR amplified products for exon 1, exon 2, and exon 4 of the SLA-DRA gene, respectively. Exon 3 with no polymorphism detected in the amplification region is not listed in this figure. (A) PCR-SSCP band patterns for exon 1 on 12% acrylamide:bisacrylamide gel (39:1) after 10 min (250 V) and then 20 h (200 V) of the run at 4°C, (B) PCR-SSCP band patterns for exon 2 on 10% acrylamide:bisacrylamide gel (29:1) after 30 min (250 V) and then 20 h (200 V) of the run at 4°C, (C) PCR-SSCP band patterns exon 4 on 10% acrylamide:bisacrylamide gel (29:1) after 30 min (250 V) and then 20 h (200 V) of the run at 12°C. PCR-SSCP, polymerase chain reaction-single strand conformation polymorphism.

Genotype and allele frequency, He and PIC of the SLA-DRA exon 1 in the three pig strains

For exon 2, genotypes AA, AB, BB, and BC were detected in all three pig breeds, but no AC or CC genotype were observed in the Duroc. Three alleles were described as: A (A-C), B (C-C), and C (C-T) at this locus (Table 3). Allele B was predominant in Large White and Landrace with frequencies of 0.51 and 0.61, respectively. The frequency of allele A (0.49) was identical with allele B in Duroc.

Genotype and allele frequency, He and PIC of the SLA-DRA exon 2 in the three pig strains

Seven different genotypes (AA, AB, BB, BC, BC, BD, CE, and CF) were identified in exon 4 of the DRA locus (Figure 1C). Seven genotypes were evident in the Large White and Landrace, but only genotypes AA, AB, and CE were detected in Duroc. Seven alleles were described as: A (T-A-A-A-G-G), B (T-G-A-A-G-G), C (C-A-A-A-G-G), D (T-A-A-A-A-G), E (T-A-A-G-G-A), and F (T-G-G-A-G-G) (Table 4). A was the most common allele in the three breeds with frequencies of 0.52, 0.54, and 0.90, respectively.

Genotype and allele frequency, He and PIC of the SLA-DRA exon 4 in the three pig strains

Results of genetic indexes (He and PIC) indicated that all PIC values ranged from 0.4 to 0.59 (Tables 2, 3, and 4) with the exception of low polymorphism in Duroc in exon 4 site (PIC = 0.19) (Table 4) indicating all of these loci exhibited abundant genetic diversity within the different pig strains according to the protocol of Vaiman et al. (1994). The PIC and He for Duroc pigs were lower than Large White and Landrace in exon 4 locus, which is not unexpected. However, it may be explained by the differences in genetic backgrounds as well as selection pressures between Duroc and Large White, Landrace.

Nucleotide and amino acid variation

The nucleotide characterizations were investigated by cloning and sequencing. According to references sequence AY303990 of SLA-DRA gene, two SNPs (c.178G>A and c.211T>C) and two SNPs (c.3093A>C and c.3104C>T) were detected in exon 1 and exon 2, respectively (Tables 2 and 3). Additionally, in exon 4 locus, one SNP (c.4081T>C) in intron 3 and five SNPs (c.4167A>G, c.4184A>G, c.4194A>G, c.4246A>G, and c.4293G>A) in exon 4 were detected (Table 4). When compared with the existing SLA-DRA alleles in the Immuno Polymorphism Database (IPD)-MHC SLA sequence database (http://www.ebi.ac.uk/ipd/mhc/sla/), mutations c.178G>A, c.211T>C, c.4081T>C c.4184A>G and c.4293G>A were found for the first time in swine, which formed four new allele sequences and have been submitted to the GenBank database (http://www.ncbi.nlm.nih.gov/nuccore/) and assigned the accession number JX135565, JX135566, JX135567, and JX135568.

Effects of genotypes and haplotypes on piglet diarrhea

Results shows differences at SLA-DRA exon 1 and exon 2 loci for genotype and piglet diarrhea (p<0.05) (Table 5). The least square means along with standard errors for piglet diarrhea score among different genotypes are given in Table 6.

Effects of different factors on piglet diarrhea score

Effects of different genotypes on piglet diarrhea score

In exon 1, the individuals with genotype BC showed a higher diarrhea score (1.77±0.28) than those with genotypes AA (0.69±0.18), AB (0.85±0.21), AC (0.97±0.11) and CC (0.52±0.19), and genotype AC showed a higher diarrhea score (0.97±0.11) than the genotype CC (0.52±0.19) (p<0.05). In exon 2, diarrhea score of genotypes AA (1.12±0.19) and BB (1.07±0.14) were significantly higher than those of AC (0.22±0.35) and CC (−0.40±0.53) (p<0.05). In exon 4, although the observed significance level for the F test was 0.1552, the genotypic linear contrasts suggested that diarrhea score was higher in the AA genotype (1.06±0.10) than AB (0.59±0.16) and BB (0.29±0.35) genotypes (p<0.05). The animals with genotype BD were excluded when association analysis was carried out because the sample size was less than three.

Haplotype construction of all SNPs in the three exons and its association with piglet diarrhea are analyzed in Table 7. A total of 54 haplotypes were found, 14 of which had frequencies greater than 2% and were defined as common haplotypes while the others were combined in this study (Breast Cancer Association Consortium, 2006). Diarrhea individuals were mainly distributed in Hap2, 5, 8, 10, and 14, which takes 47.92% of the total diarrhea individuals (46/96), therefore, those 5 haplotypes may be susceptible to piglet diarrhea; no diarrhea piglets were in Hap9, it may be the resistant haplotype.

SLA-DRA gene haplotype and its association with piglet diarrhea

DISCUSSION

Genetic polymorphism of the SLA-DRA gene

Early studies have described limited variation in the SLA-DRA locus for the NIH miniature pig (Sachs et al., 1988; Chardon et al., 1999), but subsequent researches revealed that SLA-DRA locus exhibits limited polymorphism with 13 alleles representing 3 allele groups, compared with other SLA class II genes (Ho et al., 2009; Lunney et al., 2009; Ho et al., 2010b). In this study, we detected polymorphisms throughout the entire four exon regions of SLA-DRA gene firstly, 3, 2, and 6 nucleotide mutations were detected in exon 1, exon 2, and exon 4 of DRA gene, respectively, which confirmed the result of Ho et al. (2009) and Lunney et al. (2009) who showed that the SLA-DRA gene is polymorphic.

Moreover, data on SLA class II genes now predicts that the determinants of antigen binding reside in exon 2 (Chardon et al., 2000). In the SLA-DRA locus, to date only two polymorphic sites have been found in the α1 domain (exon 2). One is an A, C, or G (the G has only been found in one allele) mutation, which results in a substitution of a methionine, leucine, or valine, respectively. The second is a C or T which both code for an arginine. The two SNPs (c.3093A>C and c.3104C>T) detected in our study are accord with the A and C polymorphisms and the C and T polymorphism. The only two nucleotide mutations in this region may be attributed to the conservative property of SLA-DRA gene, these conservative substitutions result from the survival of only mutations that retain the the crucial function of SLA-DRA gene in regulation of immunity, its conservative property is a result of evolutionary selection.

Only two homozygous genotypes (AA and BB) were present in exon 4 of the SLA-DRA gene, alleles C, D, E, and F exist in heterozygous forms. There are two factors contributing to this: i) there could be a selection effect as these alleles have negative effects on individual performance, heterozygous individuals have a greater advantage over homozygous against a variety of pathogens (Sommer et al., 2005), consequently the missing homozygous genotypes may have been eliminated in the breeding process (Beuzen et al., 2000); ii) there could be a sampling effect because the sample sizes were relatively small and few homozygous animals with these alleles were expected. Thus, this polymorphism mechanism may be a necessity for resisting environment stress and maintaining normal function of SLA-DRA gene.

Effects of genotypes on piglet diarrhea

Previously studies have established that piglet mortality rate during pre-weaning caused by diarrhea was linked to SLA haplotype (Renard and Vaiman, 1989). Our current study reported that the diarrhea score of the DQA exon 2 genotype CC (0.36±0.19) was significantly lower than AB (0.98±0.09), BB (0.85±0.77), BC (1.25±0.23) and BE (1.19±0.19) (Yang et al., 2013). In the present study, our analysis also revealed significant differences in the score of piglet diarrhea among different genotypes of SLA-DRA gene, this supports the previous observations by Renard and Vaiman (1989) and Yang et al. (2013).

The DRA gene is located at the most centromeric of SLA II gene cluster, encodes the α chain with DQA gene. The presence of a small number of haplotypes shared by multiple individuals is indicative of linkage disequilibrium (LD). Our research group had observed a significant association between DQA gene and piglet diarrhea relating to nucleotide mutations c.4086A>G and c.4087C>T, c.4089A>G and c.4090C>T as linked codons variation in allele C of Large White, Landrace and Duroc pigs (Yang et al., 2013). According to the results of Yang et al. (2013), linkage disequilibrium analysis between the exon 2 loci of SLA-DRA-DQA gene revealed five haplotypes of DRA*A-DQA*B, DRA*B-DQA*B, DRA*B-DQA*C, DRA*B-DQA*D, DRA*C-DQA*B. Absolute value of D′ and r2 for all SNP loci between DQA and DRA gene showed that of the 210 pairwise comparisons, 95 showed significant LD at p<0.01 level, indicating the comprising haplotypes of SLA-DQA and DRA gene are in complete linkage disequilibrium (|D′| = 1.0000, r2 = 0.7952). In this study, we obtained SNPs in SLA-DRA gene and significant genetic correlation with piglet diarrhea were found in these 3 observed pig breeds, which together may explain the linkage of those two SLA II α chain markers reflecting the important function in disease resistance.

Amino acid mutations in coding regions have significant impact on the stability of proteins, which largely depend on the chemical and physical properties of mutation amino acid (molecular weight, volume and hydrophobic), mutation position (whether in transmembrane domain), and charge properties (Wang and Moult, 2001). The mutations in this study were associated with three synonymous mutations and seven non-synonymous mutations, in which c.178G>A and c.4167A>G both caused changes of uncharged amino acids to positively charged amino acid (p.10Gly>Arg, p.206Gln>Arg), and c.211T>C and c.4184A>G which resulted in a hydrophilic amino acid change to hydrophobic amino acid (p.21Ser>Pro, p.212Thr>Ala). Especially, the function of transmembrane protein is influenced by single peptide, c.178G>A that occurred in highly hydrophobic areas of the SLA-DRA single peptide region, and results in a change of hydrophobic residual base to hydrophilic and therefore may lead to the functional loss of this peptide. These are notable variations as they not only occurred in an important functional region but also caused changes in amino acid properties. The constructed haplotypes of all SNPs in the three exons and their association with piglet diarrhea revealed that haplotypes Hap2 (ATACTGAAGG), Hap5 (GTACTGAAGG), Hap8 (GTCCTGAAGG), Hap10 (GTCTTAAAGG) and Hap14 (GCCCTAAAGG) may be the susceptible haplotypes to piglet diarrhea and Hap9 (GTCCCAAGGA) may be the resistant haplotype. The nucleotide substitutions may be linked to variations of each other. Thus, we ought to focus on optimizing individuals of Hap9 and firstly eliminating individuals of Hap2, 5, 8, 10, and 14 in molecular breeding.

In conclusion, the present study firstly reports five novel SNPs within the exon regions of the SLA-DRA gene, which will provide more genetic information on the SLA-DRA gene variation. Our data revealed that these detected SNPs significantly influenced piglet diarrhea, and could be used as potential genetic markers for application in pig disease resistance breeding.

ACKNOWLEDGMENTS

This study is supported by Bio-Technology Special Program of Agriculture and Animal Husbandry of Gansu Province, China (Grant No. GNSW-2008-04, No.092NKDA036 and No.0804NKCA065) and Natural Science Foundation of China (Grant No.31101682). We sincerely thank Prof. Michael Adrian Brown (Grassland Research Laboratory, USDA-ARS, EL Reno, OK, USA) for his assistance in manuscript writing.

References

Beuzen ND, Stear MJ, Chang KC. 2000;Molecular markers and their use in animal breeding. Vet J 160:42–52.
Breast Cancer Association Consortium. 2006;Commonly studied single-nucleotide polymorphisms and breast cancer: Results from the breast cancer association consortium. J Natl Cancer Inst 98:1382–1396.
Botstein D, White RL, Skolnick M, Davis RW. 1980;Construction of a genetic linkage map in man using restriction fragment length polymorphisms. Am J Hum Genet 32:314–331.
Byun SO, Fang Q, Zhou H, Hickford JGH. 2009;An effective method for silver-staining DNA in large numbers of polyacrylamide gels. Anal Biochem 385:174–175.
Chardon P, Renard C, Vaiman M. 1999;The major histocompatibility complex in swine. Immunol Rev 167:179–192.
Chardon P, Renard C, Gaillard CR, Vaiman M. 2000;The porcine major histocompatibility complex and related paralogous regions: A review. Genet Sel Evol 32:109–128.
Chen FX, Xie J, Zhou Y, Li NL, Chou KY. 2004;Novel SLA-DR alleles of three Chinese pig strains and the related function in human T cell response. Cell Mol Immunol 1:212–218.
Essler SE, Ertl W, Deutsch J, Ruetgen BC, Groiss S, Stadler M, Wysoudil B, Gerner W, Ho CS, Saalmueller A. 2013;Molecular characterization of swine leukocyte antigen gene diversity in purebred Pietrain pigs. Anim Genet 44:202–205.
Hedrick PW, Kim TJ. 2000. Genetics of complex polymorphisms: Parasites and maintenance of MHC variation. Evolutionary Genetics from Molecules to Morphology In : Singh RS, Krimbas CB, eds. Cambridge, NY, USA: p. 204–234.
Ho CS, Lunney JK, Ando A, Rogel-Gaillard C, Lee JH, Schook LB, Smith DM. 2009;Nomenclature for factors of the SLA system, update 2008. Tissue Antigens 73:307–315.
Ho CS, Martens GW, Amoss MS Jr, Gomez-Raya LN, Beattie CW, Smith DM. 2010a;The SLA alleles and haplotypes present in Sinclair and Hanford swine. Dev Comp Immunol 34:250–257.
Ho CS, Lunney JK, Lee JH, Franzo-Romain MH, Martens GW, Rowland RR, Smith DM. 2010b;Molecular characterization of swine leucocyte antigen class II genes in outbred pig populations. Anim Genet 41:428–432.
Hudson RR. 2001. Linkage Disequilibrium and Recombination. Handbook of Statistical Genetics In: Handbook of Statistical Genetics. In : Balding DJ, Bishop M, Cannings C, eds. John Wiley and Sons, Ltd. Chichester, UK: p. 309–324.
Kelly DJ, O’Brien JJ, McCracken KJ. 1990;Effect of creep feeding on the incidence, duration and severity of post-weaning diarrhea in pigs. Res Vet Sci 49:223–228.
Lewis CR, Ait-Ali T, Clapperton M, Archibald AL, Bishop S. 2007;Genetic perspectives on host responses to porcine reproductive and respiratory syndrome (PRRS). Viral Immunol 20:343–358.
Luetkemeier ES, Malhi RS, Beever JE, Schook LB. 2009;Diversification of porcine MHC class II genes: evidence for selective advantage. Immunogenetics 61:119–129.
Lunney JK, Ho CS, Wysocki M, Smith DM. 2009;Molecular genetics of the swine major histocompatibility complex, the SLA complex. Dev Comp Immunol 33:362–374.
Mallard BA, Wilkie BN, Kennedy BW. 1989;Genetic and other effects on antibody and cell mediated immune response in swine leucocyte antigen (SLA)-defined miniature pigs. Anim Genet 20:167–178.
Molina RM, Cha SH, Chittick W, Lawson S, Murtaugh MP, Nelson EA, Christopher-Hennings J, Yoon KJ, Evans R, Rowland RRR, Zimmerman JJ. 2008;Immune response against porcine reproductive and respiratory syndrome virus during acute and chronic infection. Vet Immunol Immunopathol 126:283–292.
Morris RS, Davies PR, Lawton DE. 2002. Evolution of diseases in the world’s pig industry. In : 17th International Pig Veterinary Society Congress Proceedings. Ames, IA, USA. p. 1–10.
Patch JR, Pedersen LE, Toka FN, Moraes M, Grubman MJ, Nielsen M, Jungersen G, Soren B, Golde WT. 2011;Induction of foot-and-mouth disease virus-specific cytotoxic T cell killing by vaccination. Clin Vaccine Immunol 18:280–288.
Renard C, Vaiman M. 1989;Possible relationships between SLA and porcine reproduction. Reprod Nutr Dev 29:569–576.
Sachs DH, Germana S, Ei-Gamil M, Gustafsson K, Hirsch F, Pratt K. 1988;Class II genes of miniature swine. Immunogeneties 28:22–29.
Sambrook J, Russell DW. 2001. Molecular Cloning: A Laboratory Manual 3rd edth ed. Cold Spring Harbor: Cold Spring Harbor Laboratory Press.
SAS. 1999. SAS User’s Guide: Statistics (Version 8.01 Ed) SAS Institute Inc. Cary, NC, USA:
Shinkai H, Arakawa A, Tanaka-Matsuda M, Ide-Okumura H, Terada K, Chikyu M, Kawarasaki T, Ando A, Uenishi H. 2012;Genetic variability in swine leukocyte antigen class II and Toll-like receptors affects immune responses to vaccination for bacterial infections in pigs. Comp Immunol Microbiol Infect Dis 35:523–532.
Sommer S. 2005;The importance of immune gene variability (MHC) in evolutionary ecology and conservation. Front Zool 2:16.
Stephens M, Donnelly P. 2003;A comparison of Bayesian methods for haplotype reconstruction from population genotype data. Am J Hum Genet 73:1162–1169.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. 2011;MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol 28:2731–2739.
Vaiman D, Mecier D, Moazami-Goudarzi K, Eggen A, Ciampolini R, Lepingle A, Velmala R, Kaukinen J, Varvio SL, Martin P. 1994;A set of 99 cattle microsatellites: characterization, synteny mapping, and polymorphism. Mamm Genome 5:288–297.
Wang Z, Moult J. 2001;SNPs, proteins structer, and disease. Hum Mutat 17:263–270.
Yang QL, Kong JJ, Zhao SG, Liu LX, Wang DW, Li FS, Jiang TT, Feng Y. 2012;Polymorphism of SLA-DRA gene exon 2 and its correlation with piglet diarrhea. Acta Vet Zootech Sin 43:1020–1027. [in Chinese, English abstract].
Yang QL, Kong JJ, Wang DW, Zhao SG, Gun SB. 2013;Swine leukocyte antigen-DQA gene variation and its association with piglet diarrhea in Large White, Landrace and Duroc. Asian Australas J Anim Sci 26:1065–1071.
Zhou H, Hickford JGH. 2008;Clonal polymerase chain reaction-single-strand conformational polymorphism analysis: An effective approach for identifying cloned sequences. Anal Biochem 378:111–112.

Article information Continued

Figure 1

PCR-SSCP analysis of PCR amplified products for exon 1, exon 2, and exon 4 of the SLA-DRA gene, respectively. Exon 3 with no polymorphism detected in the amplification region is not listed in this figure. (A) PCR-SSCP band patterns for exon 1 on 12% acrylamide:bisacrylamide gel (39:1) after 10 min (250 V) and then 20 h (200 V) of the run at 4°C, (B) PCR-SSCP band patterns for exon 2 on 10% acrylamide:bisacrylamide gel (29:1) after 30 min (250 V) and then 20 h (200 V) of the run at 4°C, (C) PCR-SSCP band patterns exon 4 on 10% acrylamide:bisacrylamide gel (29:1) after 30 min (250 V) and then 20 h (200 V) of the run at 12°C. PCR-SSCP, polymerase chain reaction-single strand conformation polymorphism.

Table 1

Primer sequences, annealing temperature (Ta) and size of amplification products

Locus Primer sequence (5′→3′) Ta (°C) Position Size (bp)
Exon 1 F: CTTTGCTTGTATTGC
R: ACCTAACTACCCCTC
56.8 79–264 186
Exon 2 F: CAGAGAATCACGTGATCAT
R: ACAGGTACCATTGGTGTT
55.0 288–3137 257
Exon 3 F: TGCTAAACAGGGAAGGCT
R: ACAAAGGAGACTGAGGGATG
56.8 3531–3882 352
Exon 4 F: TCCCGTAATACATCGTTC
R: TTCCTTTCCTTGGCTCAT
55.6 4018–4374 357

Primer positions refer to SLA-DRA gene sequence in GenBank (Accession No.AY303990), exon 1, 2, 3, and 4 are located at 151–226, 2884–3129, 3573–3854, 4155–4309, respectively.

Table 2

Genotype and allele frequency, He and PIC of the SLA-DRA exon 1 in the three pig strains

Exon1 GenBank (No.AY303990) 5′→3′ Frequency


178 211 Large White Landrace Duroc
72 72 72
Genotype AA A T 0.17 (12) 0.13 (9) 0.18 (13)
AB A/G T/C 0.01 (1) 0.06 (4) 0.35 (25)
BB G C 0 (0) 0 (0) 0.10 (7)
BC G T/C 0.04 (3) 0.04 (3) 0.11 (8)
AC A/G T 0.50 (36) 0.71 (51) 0.17 (12)
CC G T 0.28 (20) 0.07 (5) 0.10 (7)
Allele A A T 0.42 0.51 0.44
B G C 0.03 0.05 0.33
C G T 0.55 0.44 0.24
He 0.51 0.54 0.67
PIC 0.40 0.44 0.57

He, observed heterozygosity; PIC, polymorphism information content; SLA, swine leukocyte antigen. The following marker means are same. Numbers in brackets are population size tested that belong to the respective genotypes. Bases in the boxes are mutation bases in amplified regions.

Table 3

Genotype and allele frequency, He and PIC of the SLA-DRA exon 2 in the three pig strains

Exon2 GenBank (No. AY303990) 5′→3′ Frequency


3093 3104 Large White Landrace Duroc
72 72 72
Genotype AA A C 0.13 (9) 0.08 (6) 0.21 (15)
AB A/C C 0.33 (24) 0.33 (24) 0.57 (41)
BB C C 0.26 (19) 0.38 (27) 0.18 (13)
BC C C/T 0.17 (12) 0.14 (10) 0.04 (3)
AC A/C C/T 0.07 (5) 0.06 (4) 0 (0)
CC C T 0.04 (3) 0.01 (1) 0 (0)
Allele A A C 0.33 0.28 0.49
B C C 0.51 0.61 0.49
C C T 0.16 0.11 0.02
He 0.61 0.54 0.52
PIC 0.53 0.47 0.40

He, observed heterozygosity; PIC, polymorphism information content; SLA, swine leukocyte antigen.

Table 4

Genotype and allele frequency, He and PIC of the SLA-DRA exon 4 in the three pig strains

Exon4 GenBank (No. AY303990) 5′→3′ Frequency


4081 4167 4184 4194 4246 4293 Large White Landrace Duroc
T 72 72 72
Genotype AA T A A A G G 0.35 (25) 0.43 (31) 0.85 (61)
AB T A/G A A G G 0.35 (25) 0.22 (16) 0.97 (7)
BB T G A A G G 0.07 (5) 0.06 (4) 0 (0)
BC T/C A/G A A G G 0.10 (7) 0.06 (4) 0 (0)
BD T A/G A A A/G G 0.01 (1) 0.01 (1) 0 (0)
CE T/C A A A/G G G/A 0.08 (6) 0.18 (13) 0.06 (4)
CF T/C A/G A/G A G G 0.04 (3) 0.04 (3) 0 (0)
Allele A T A A A G G 0.52 0.54 0.90
B T G A A G G 0.30 0.20 0.05
C C A A A G G 0.11 0.14 0.03
D T A A A A G 0.01 0.01 0
E T A A G G A 0.04 0.09 0.03
F T G G A G G 0.02 0.02 0
He 0.63 0.64 0.19
PIC 0.57 0.59 0.19

He, observed heterozygosity; PIC, polymorphism information content; SLA, swine leukocyte antigen.

SNP c.4081T>C was found in intron 3.

Table 5

Effects of different factors on piglet diarrhea score

Factors Exon1 Exon2 Exon4



Breed Sex Genotype Breed Sex Genotype Breed Sex Genotype
F value 0.65 0.01 3.21 0.34 0.02 2.12 0.65 0.03 1.58
p 0.52 0.93 0.008* 0.71 0.90 0.04* 0.52 0.86 0.16
*

Significant at 0.05 level.

Table 6

Effects of different genotypes on piglet diarrhea score

Exon1 Exon2 Exon4



Genotype Number LSM±SE Genotype Number LSM±SE Genotype Number LSM±SE
AA 34 0.69±0.18bc AA 30 1.12±0.19a AA 117 1.06±0.10a
BB 7 0.83±0.41abc BB 59 1.07±0.14a BB 9 0.29±0.35b
AB 30 0.85±0.21bc AB 89 0.83±0.11ab AB 48 0.59±0.16b
BC 14 1.77±0.28a BC 25 0.79±0.21ab BC 11 0.76±0.32ab
AC 99 0.97±0.11b AC 9 0.22±0.35b BD 2 0.97±0.74ab
CC 32 0.52±0.19c CC 4 −0.40±0.53b CE 23 1.05±0.22ab
CF 6 0.63±0.43ab

LSM±SE, Diarrhea score was estimated for each genotype class and shown as means least square means±standard error.

a,b

Values with different superscripts within the same line differ significantly at p<0.05.

Table 7

SLA-DRA gene haplotype and its association with piglet diarrhea

Name Haplotypes Frequency (%) Diarrhea piglets

Number Frequency (%)
Hap1 ATACTAAAGG 8.63 4 4.17
Hap2 ATACTGAAGG 4.27 7 7.29
Hap3 ATCCTGAAGG 2.60 1 1.04
Hap4 GTACTAAAGG 8.89 3 3.13
Hap5 GTACTGAAGG 4.13 8 8.33
Hap6 GTACCAAGGA 2.08 4 4.17
Hap7 GTCCTAAAGG 9.67 4 4.17
Hap8 GTCCTGAAGG 8.13 13 13.54
Hap9 GTCCCAAGGA 2.64 0 0
Hap10 GTCTTAAAGG 6.25 10 10.42
Hap11 GTCTCAAGGA 3.28 3 3.12
Hap12 GTCTCGGAGG 2.64 3 3.12
Hap13 GCACTGAAGG 6.64 5 5.21
Hap14 GCCCTAAAGG 6.54 8 8.33
Others 23.61 23 23.96

SLA, swine leukocyte antigen.