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

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.

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: CTTTGCTTGTATTGCR: ACCTAACTACCCCTC	56.8	79–264	186
Exon 2	F: CAGAGAATCACGTGATCATR: ACAGGTACCATTGGTGTT	55.0	288–3137	257
Exon 3	F: TGCTAAACAGGGAAGGCTR: ACAAAGGAGACTGAGGGATG	56.8	3531–3882	352
Exon 4	F: TCCCGTAATACATCGTTCR: 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.

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:

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.

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.

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.

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.

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.

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.

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.

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.

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

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.

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.

DISCUSSION

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.