the online version of this article ( doi:10.1007/s13353 - 014 - 0199 - 8 ) contains supplementary material , which is available to authorized users . salmonella is a widespread foodborne pathogen with the ability to adapt to different environments , consequently creating significant challenges to food - producing industries in controlling this pathogen in food chain products . swine ( sus scrofa ) are an important reservoir of salmonella because colonization and shedding of this bacterium occurs within asymptomatic pigs , imposing elevated risks to public and animal health . thus , diverse intervention strategies are needed to control the transmission of salmonella from pig products to humans and to the environment . in bacterial infections , the severity of infection is impacted by the pathogenicity of the microorganism and its interaction with the host immune defense system ( zanella et al . toll - like receptor 4 ( tlr4 ) is a well - characterized gram - negative bacterial lipopolysaccharide ( lps ) recognition receptor and a host inflammatory response activator well conserved among animal species ( noreen et al . 2012 ; yang et al . 2012 ) . schrder and schumann ( schrder and schumann 2005 ) suggested that mutations in the tlr4 regions involved with pathogen recognition and transduction signaling may affect host susceptibility to infection . polymorphisms in the tlr4 gene have been associated with different infectious diseases in humans , such as meningitis and tuberculosis , as well as some types of cancers ( noreen et al . 2012 ) and with infection and disease in cattle , chicken and pigs ( yang et al . tlr4 is located on sus scrofa 1 ( ssc1 ) v10.2 ( 289,776,058  bp to 289,785,087  bp ) . ( 2006 ) identified the genomic structure of porcine tlr4 , and shinkai et al . specifically for tlr4 , 13 snps were widely distributed in 11 pig breeds , and of those , seven were non - synonymous . thirty four snps were identified in tlr4 using pigs representing european commercial breeds and some traditional breeds ( n  =  259 ) , and of these , 17 snps were located in the non - coding region and 17 snps were found in the coding region ( palermo et al . furthermore , polymorphisms in the tlr4 gene have been identified as potential genetic markers for disease susceptibility in pigs ( uenishi and shinkay 2009 ) . our collaborative group has reported up - regulation of tlr4 and its target genes in pigs challenged with salmonella enterica serovar typhimurium ( huang et al . 2011 ) . therefore , to determine if tlr4 is a possible candidate gene associated with salmonella shedding , we first , identified snps in the tlr4 gene of our previously described low and persistent shedder pig populations ( huang et al . 2011 ; uthe et al . second , we investigated associations of the tlr4 snps with salmonella shedding status . selecting for pigs with reduced salmonella fecal shedding would decrease environmental contamination and lower pathogen transmission to other animals and humans ; thus , identification of loci in tlr4 associated with salmonella fecal shedding is the focus of this study . all procedures involving animals in the nadc-40 and nadc-77 populations were approved by the usda , ars , nadc animal care and use committee . briefly , all the pigs used in this study were intranasally challenged at 7  weeks of age with 1    10  cfu of s. typhimurium 4232 as previously described ( huang et al . 2011 ; uthe et al . 2009 ) . at days 2 , 7 , 14 , and 20/21 post - inoculation ( pi ) , salmonella fecal shedding was quantified using a standard bacteriological test previously described ( uthe et al . 2009 ) . of the initial 117 animals , 40 ( n  =  40 ) pigs were chosen based on their extreme fecal culture status ; quantitative classification of the phenotype was scored based on cumulative salmonella fecal shedding on days 2 , 7 , 14 , and 20/21 pi ( huang et al . 2011 ) . genomic dna was extracted from blood samples and purified as previously described ( uthe et al . nine sets of primers were designed based on the ensembl gene sequence for enssscg00000005503 using beacon designer ( table  1 ) . primers were selected to cover all exons ( n  =  3 ) of tlr4 including a 713  bp upstream region and a 1225  bp downstream region ( ssc1 : 289,775,345  bp289,786,312  bp ensembl genome build 10.2 ) . the pcr mix for each reaction contained 16.75  l dh2o , 1.25  l each primer ( 10  m ) , 2  l of dntps , 2.5  l 10x buffer with mgcl2 , 0.25  l platinum taq dna polymerase high fidelity ( invitrogen grand island , ny , usa ) and 1  l of 10  ng/l pig dna . dna samples were amplified using the mj research ptc-200 pcr thermal cycler ( biorad laboratories , hercules , ca ) . the pcr reaction was performed as follows : 94  c for 2  min , 30 cycles of 94  c for 30  s , 58  c for 30  s , 72  c for 1  min and a final step of 72  c for 7  min . pcr products were visualized by agarose gel electrophoresis to confirm a single correct size product and purified using minelute 96 uf pcr purification kit ( qiagen ) prior to dna sequencing using an ab 3730xl dna analyzer ( applied biosystems ) at iowa state university , ames , ia.table 1identified snps and position in the tlr4 gene of salmonella low and persistent shedder pigssingle marker association ( p - value)genbank accession number for the snplocation in sus scrofa genome ( bp)amino acidnlocationprimerssnp designationqualitativequantitativesiteamino acid15upstream5gaaccatgcagtagaacagg3ctggaagtctgtagtcaagg 
 5u : a-1082g#0.0330.064nossc1:289,774,9832 
 5u : t-1019c#0.0330.064nossc1:289,775,0463 
 5u : c-984t#0.0330.064nossc1:289,775,08145cacaagaaggaagagatagc3 caccaagggaagctctagg 
 5u : c-522t#0.1330.244nossc1:289,775,5435 
 5u : g-400a0.3630.550rs80830544ssc1:289,775,6656 
 5u : g-75c#0.0250.056nossc1:289,775,9797intron 25acagaagattggatggaagga3 gagataagaaagctgagacc 
 2i : a232c0.0040.029rs80881287ssc1:289,780,2268 
 2i : c298t**0.0020.013rs80787918ssc1:289,780,2929intron 25cctcacttgatatgtttgcc3gttcctccaggacagatttg 
 2i : c2567t#0.0010.025nossc1:289,782,76110exon 3 
 c318a0.0030.037rs80923358ssc1:289,782,83411 
 g417a0.0030.037rs80951861ssc1:289,782,93312 
 t611a*0.0070.054rs80811682ssc1:289,783,127204l / h13exon 35attcaaggtctggctggttc3 tgaagacatcaggaagcaag 
 g826a*0.2850.514shinkai et al . ( 2006)ssc1:289,783,342276v / i14 
 g960a0.0640.105rs80981701ssc1:289,783,47615 
 g962a*0.0340.046rs80955017ssc1:289,783,478321r / h16 
 c1027a*0.1760.231rs80894552ssc1:289,783,543343q / kexon 35acatccacgttgtcttccg3cagttcattcctcacccag17exon 35cttcctcctggtatctgtgg3ggcagtcctgtgtatctcg 
 g2397a0.0250.056rs80834103ssc1:289,784,913183downstream5actcccaacgtgtcccttg3ccaagaagtgccactttcaac 
 3d : c208t**0.0020.011rs80907449ssc1:289,785,250 
 to first codon of exon 1 ; position in intron 2 ; position in coding region ; position in 3utr downstream of last codon ; * non - synonymous snps ; * * haplotype components ; # novel snps identified snps and position in the tlr4 gene of salmonella low and persistent shedder pigs 
 to first codon of exon 1 ; position in intron 2 ; position in coding region ; position in 3utr downstream of last codon ; * non - synonymous snps ; * * haplotype components ; # novel snps sequences were analyzed and polymorphisms were identified using phred / phrap / consed / polyphred software ( nickerson et al . snps were removed if the minor allele frequency ( maf ) was less than 10  % , if the snps failed to genotype in more than 10  % of the samples , or if the snps failed the hardy  weinberg equilibrium ( p  <  0.001 ) . statistical analyses were conducted within plink and r statistical environment ( version 1.07 , ( purcell et al . 2007 ) ) . a chi - squared test (  ) was used to test associations of snps located in tlr4 and the qualitative measurement of salmonella shedding ( persistent versus low ) . the wald test was used to verify associations with salmonella shedding as a quantitative trait . a significance threshold for the association analysis was set to p    0.05 . following the single marker association test , a haplotype test was conducted within plink to identify if a haplotype was more informative than a single snp . first , an omnibus association test was performed to identify the overall association of the haplotype with the qualitative measurement of salmonella shedding . if an association was identified , a haplotype - specific test was performed to identify which combination of the alleles provided the strongest evidence for an association with salmonella shedding in swine . following haplotype construction , a stepwise regression using a backward - elimination process was performed to identify the effect of each associated snp in relationship to the haplotype ; in this test , all associated snps were included and excluded individually from the analysis , and the association of the haplotype was tested each time using plink . huang et al . ( 2011 ) identified the tlr4-dependent set of genes ( tlr4 regulon ) as a major inducer of the transcriptional response in salmonella persistently shedding pigs , and this tlr4 regulon was not significantly affected in the low shedding pigs . thus , tlr4 is considered a potential candidate to analyze the association of genetic polymorphisms with the diverse phenotypic patterns of salmonella shedding in swine . in this study , two swine populations were investigated , nadc-40 and nadc-77 ( uthe et al . 2011 ) , each population with 10 low and 10 persistent salmonella shedding animals ( fig .  1 ) . for the quantitative measurement of salmonella shedding per pig , a cumulative measurement taken within days 2 , 7 , 14 , and 20/21 pi was calculated ( huang et al . sequencing analysis of those 40 ( n  =  40 ) animals identified 18 snps ; 12 were previously described in the literature and/or annotated in genbank and six are novel snps ( table  1 ) . five of the novel snps are located within the 5untranslated region ( utr ) and one is within intron 2.fig . 1area under the log curve illustrating the log of cumulative colony forming units ( cfu ) . quantitative bacteriology of salmonella shedding in swine fecal samples was performed at day 2 , 7 , 14 , and 20/21  days post - challenge with salmonella enterica serovar typhimurium , and cfu were determined area under the log curve illustrating the log of cumulative colony forming units ( cfu ) . quantitative bacteriology of salmonella shedding in swine fecal samples was performed at day 2 , 7 , 14 , and 20/21  days post - challenge with salmonella enterica serovar typhimurium , and cfu were determined the swine tlr4 gene ( 9030  bp / sgsc sscrofa10.2/susscr3 ) is composed of three exons ( 93 , 167 , and 2266  bp ) . taking together these results and the literature , 50 snps have been identified in tlr4 , with 22 snps located in the coding regions ( thomas et al . our investigation identified four of those nine non - synonymous snps in exon 3 segregating in the nadc-40 and nadc-77 pig populations . of the 18 snps identified in the two pig populations , 13 ( n  =  13 ) snps were associated ( p    0.05 ) with salmonella shedding as a qualitative phenotype using a chi - squared test ; of those 13 snps , seven were also associated with salmonella shedding as a quantitative phenotype using a wald statistical test ( table  1 ) . using a haplotype construction and the backward - elimination process , the most significant haplotype for both measurements of salmonella shedding , qualitative ( p    7.9    10 ) and quantitative ( p    4.0    10 ) ( table  2 ) comprised a region of 4.9  kb composed of snps , rs80787918 ( snp8 ) and rs80907449 ( snp18 ) ( r  =  0.902 ) located at ssc1:289,780,292  bp and ssc1:289,785,250  bp , respectively ( table  1).table 2haplotypes frequency ( snps rs80787918 and rs80907449 ) and associations with qualitative and quantitative phenotypes of salmonella sheddinghaplotype frequencyqualitative 
 p valuequantitative 
 p valuehaplotypepersistent shedderslow shedderscc0.6750.30.000790.004201tc00.0250.31430.1054ct00.0250.31430.1445tt0.3250.650.003340.02912 haplotypes frequency ( snps rs80787918 and rs80907449 ) and associations with qualitative and quantitative phenotypes of salmonella shedding four snps , rs80811682 ( snp12 ) , snp13 , rs80955017 ( snp15 ) , and rs80894552 ( snp16 ) , located on exon 3 of tlr4 gene are non - synonymous mutations and they are positioned between markers rs80787918 ( snp8 ) and rs80907449 ( snp18 ) . when the additive effect of those markers was tested within the haplotype constructed with markers rs80787918 ( snp8 ) and rs80907449 ( snp18 ) , we did not observe any improvement in the association test . however , analyzing together the markers rs80787918 ( snp8 ) , rs80811682 ( snp12 ) , snp13 , rs80955017 ( snp15 ) , rs80894552 ( snp16 ) and rs80907449 ( snp18 ) , the haplotype composed of alleles ( ctggcc ) was found in higher frequency ( 65  % ) in persistent shedders than 31  % in low shedders pigs ( p  <  0.003 ) . possibly , the addition of more markers in the haplotype is being penalized by the increased number of degrees of freedom and reduced number of samples per each class affecting the significance of our association results . when markers snp13 and rs80894552 ( snp16 ) , which were not significant in the single marker association test , were removed from the haplotype , the significance improved to p  <  0.001 . the haplotype composed of alleles ctgc was found in 65  % of the high shedders and 30  % of low shedders and the haplotype ( taat ) was found in 20  % of the high shedders and 41  % of the low shedders ( p  <  0.04 ) . finally , a specific haplotype ( ctc ) constructed with markers rs80787918 ( snp8 ) , rs80811682 ( snp12 ) and rs80907449 ( snp18 ) , was observed in 65  % of the persistent shedders and 30  % on the low shedders pigs ( p  <  0.001 ) . the opposite haplotype ( tat ) was observed in 27.5  % of the persistent shedders and in 59.4  % of the low shedders pigs ( p  <  0.003 ) . a trend was observed between haplotypes constructed with markers : rs80787918 ( snp8 ) , snp13 and rs80907449 ( snp18 ) ( cgc ) ; rs80787918 ( snp8 ) , rs80955017 ( snp15 ) and rs80907449 ( snp18 ) ( cgc ) ; rs80787918 ( snp8 ) , rs80894552 ( snp16 ) and rs80907449 ( snp18 ) ( ccc ) , where they were observed in 67.5  % of the persistent shedders and 30  % of low shedder pigs . haplotype cc of snps rs80787918 ( snp8 ) and rs80907449 ( snp18 ) was identified in higher frequency in persistent shedding pigs ( 67.5  % : n  =  14 ) compared to low shedding pigs ( 30  % ; n  =  6 ) ; furthermore , the frequency of haplotype tt in low shedding pigs ( 65  % ; n  =  13 ) was greater when compared to persistent shedding pigs ( 32.5  % ; n  =  6 ) . no animals from the persistent shedding group were identified with the haplotype tc or ct , while it was observed in low frequency in the low shedding group ( 2.5  % ) . together , these results suggest that the region located between markers rs80787918 and rs80907449 , more specifically on exon 3 , is possibly harboring the causative mutation for salmonella colonization and shedding variation in swine . the results from this study support the concept that tlr4 is an important modulator associated with the porcine response to salmonella infection in swine . particularly interesting is that the haplotype with the highest significant association to the shedding phenotypes was found most often ( 65  % ) in the persistent shedder pigs than in low shedder pigs . genetic variation in molecular functional regions , such as a ligand recognition site , can alter host resistance / susceptibility to specific pathogens ( uenishi et al . furthermore , synonymous mutations in a gene can play a significant role in transcriptional regulation ( sauna and kimchi - sarfaty 2011 ; sato et al . ( 2011 ) who demonstrated polymorphisms in tlr5 and tlr2 alter the cellular response to s. choleraesuis , our results highlight the importance of linking genetic variations that may influence the molecular function of a key transcriptional regulator ( tlr4 ) with salmonella shedding in swine . mention of trade names or commercial products in this article is solely for the purpose of providing specific information and does not imply recommendations or endorsement by the u.s .
toll - like receptor 4 ( tlr4 ) is a key factor in the innate immune recognition of lipopolysaccharide ( lps ) from gram - negative bacteria . 
 previous studies from our group identified differences in the expression profile of tlr4 and genes affected by the tlr4 signaling pathway among pigs that shed varying levels of salmonella , a gram - negative bacterium . therefore , genetic variation in this gene may be involved with the host s immune response to bacterial infections . the current study screened for single nucleotide polymorphisms ( snps ) in the tlr4 gene and tested their association with salmonella fecal shedding . 
 pigs ( n  =  117 ) were intranasally challenged at 7  weeks of age with 1    109  cfu of s. typhimurium 4232 and were classified as low or persistent salmonella shedders based on the levels of salmonella being excreted in fecal material . 
 salmonella fecal shedding was determined by quantitative bacteriology on days 2 , 7 , 14 , and 20/21 post exposure , and the cumulative levels of salmonella were calculated to identify the low ( n  =  20 ) and persistent ( n  =  20 ) salmonella shedder pigs . from those 40 animals , 
 the tlr4 region was sequenced , and 18 single nucleotide polymorphisms ( snps ) in tlr4 were identified . 
 twelve snps have been previously described and six are novel snps of which five are in the 5 untranslated region and one is in intron 2 . 
 single marker association test identified 13 snps associated with the qualitative trait of salmonella fecal shedding , and seven of those snps were also associated with a quantitative measurement of fecal shedding ( p  <  0.05 ) . using a stepwise regression process , a haplotype composed of snps rs80787918 and rs80907449 ( p    4.0    103 ) spanning a region of 4.9  kb was identified , thereby providing additional information of the influence of those snps on salmonella fecal shedding in pigs.electronic supplementary materialthe online version of this article ( doi:10.1007/s13353 - 014 - 0199 - 8 ) contains supplementary material , which is available to authorized users .