Genetic variability in swine leukocyte antigen class II and Toll-like receptors affects immune responses to vaccination for bacterial infections in pigs

Single-Cell RNA Sequencing Reveals Microevolution of the Stickleback Immune System

The risk and severity of pathogen infections in humans, livestock, or wild organisms depend on host immune function, which can vary between closely related host populations or even among individuals. This immune variation can entail between-population differences in immune gene coding sequences, copy number, or expression. In recent years, many studies have focused on population divergence in immunity using whole-tissue transcriptomics. But, whole-tissue transcriptomics cannot distinguish between evolved differences in gene regulation within cells, versus changes in cell composition within the focal tissue. Here, we leverage single-cell transcriptomic approaches to document signatures of microevolution of immune system structure in a natural system, the three-spined stickleback (Gasterosteus aculeatus). We sampled nine adult fish from three populations with variability in resistance to a cestode parasite, Schistocephalus solidus, to create the first comprehensive immune cell atlas for G. aculeatus. Eight broad immune cell types, corresponding to major vertebrate immune cells, were identified. We were also able to document significant variation in both abundance and expression profiles of the individual immune cell types among the three populations of fish. Furthermore, we demonstrate that identified cell type markers can be used to reinterpret traditional transcriptomic data: we reevaluate previously published whole-tissue transcriptome data from a quantitative genetic experimental infection study to gain better resolution relating infection outcomes to inferred cell type variation. Our combined study demonstrates the power of single-cell sequencing to not only document evolutionary phenomena (i.e., microevolution of immune cells) but also increase the power of traditional transcriptomic data sets.

Vaccination Failures in Pigs-The Impact of Chosen Factors on the Immunisation Efficacy

Infectious diseases that often lead to economic losses still pose a severe problem in the pig production sector. Because of increasing restrictions on antibiotic usage, vaccines may become one of the major approaches to controlling infectious diseases; much research has proved that they could be very efficient. Nevertheless, during their life, pigs are exposed to various factors that can interfere with vaccination efficacy. Therefore, in the present paper, we reviewed the influence of chosen factors on the pig immunisation process, such as stress, faecal microbiota, host genetics, the presence of MDAs, infections with immunosuppressive pathogens, and treatment with antibiotics and mycotoxins. Many of them turned out to have an adverse impact on vaccine efficacy.

Polymorphisms in Pattern Recognition Receptor Genes Are Associated with Respiratory Disease Severity in Pig Farms

Reduced productivity caused by infections, particularly respiratory diseases, is a serious problem in pig farming. We have previously reported polymorphisms in porcine pattern recognition receptor genes affecting molecular functions and demonstrated that the 2197A/C polymorphism in the nucleotide-binding oligomerization domain containing 2 (NOD2) gene influences porcine circovirus 2-induced mortality. Here, we investigated how these polymorphisms affect respiratory disease-induced lesions, using samples from a slaughterhouse dealing with pigs from two farms. Lung lesions were evaluated using two scoring systems, Goodwin (GW) and slaughterhouse pleuritis evaluation system (SPES), to determine the influence of Mycoplasma hyopneumoniae (Mhp) and Actinobacillus pleuropneumoniae (App), respectively. SPES scores were significantly higher when the 1205T allele of Toll-like receptor 5 (TLR5-1205T), rather than TLR5-1205C, was present. On the farm with more severe Mhp invasion, lower GW lesion scores were significantly associated with the presence of the NOD-like receptor family pyrin domain containing 3 (NLRP3)-2906G allele; where App invasion was worse, lower SPES scores were significantly associated with the presence of the NOD2-2197C allele. Combinations of polymorphisms in pattern recognition receptor genes can therefore be utilized for breeding for resistance against respiratory diseases in pigs. DNA markers of these polymorphisms can thus be used to improve productivity by reducing respiratory diseases due to bacterial pathogens in pig livestock.

Regulation and Molecular Mechanism of TLR5 on Resistance to Escherichia coli F18 in Weaned Piglets

Toll-like receptor 5 (TLR5) plays an important role in immune system. In this study, we performed transcriptome analysis of the duodenum in E. coli F18-resistant and -sensitive Sutai weaned piglets and analyzed the differential expression of TLR5. The cellular localization of TLR5 was investigated, and the effect of TLR5 expression on E. coli invasion was evaluated after pig small intestinal epithelial cell lines (IPEC-J2) were stimulated by E. coli. The results showed that TLR5 expression level in duodenum and jejunum were significantly higher in E. coli F18-sensitive than in E. coli F18-resistant piglets. TLR5 protein was mainly expressed in the cytoplasm and cell membrane. The expression of genes associated with the TLR5 signaling pathway were significantly higher in TLR5-overexpressed cells than in control cells. Bacterial adhesion was higher in TLR5-overexpressed cells than in blank cells and lower in TLR5 interference than in blank cells. The core promoter region of TLR5 included two CpG islands and 16 acting elements. The methylation of the mC-6 site in the second CpG island of the promoter region had a regulatory effect on TLR5 expression. Therefore, TLR5 plays an important regulatory role on E. coli invasion. Low expression of TLR5 inhibited the immune response and decreased cell damage, which was conducive to the resistance to E. coli stimulation. In conclusion, this study preliminarily revealed the molecular mechanism of TLR5 gene regulating the resistance of piglets to Escherichia coli, and provided a new candidate gene for screening Escherichia coli resistance markers in pigs.

Influence of swine leukocyte antigen haplotype on serum antibody titers against swine erysipelas vaccine and reproductive and meat production traits of SLA-defined selectively bred Duroc pigs

We investigated possible associations of SLA class II haplotypes with serum antibody titers against a swine erysipelas vaccine, reproductive and meat production traits using a population of selective breeding Duroc pigs. In the selective breeding Duroc pigs, four SLA class II-DRB1 and -DQB1 alleles were assigned by using PCR-sequence specific primer technique. Low-resolution haplotype (Lr)-0.30 and/or Lr-0.13 were deduced from the SLA class II alleles in the population of SLA-defined Duroc pigs. SLA-homozygous piglets with the Lr-0.30 haplotype had relatively lower serum antibody titers against the vaccine compared to those with Lr-0.13. In contrast, there were no statistically significant differences in reproductive performance between the SLA-defined pigs with two SLA class II haplotypes. Weaning and rearing rates until the body weight of 105 kg was reached in homozygous piglets with Lr-0.30 were significantly lower than those in homozygous piglets with Lr-0.13. The SLA-defined pigs had lower birth and weaning weights, body weights at 60 days of age, and daily weight gains than non-selective breeding Duroc pigs. Furthermore, the SLA-defined pigs had slightly lower back fat thickness compared to the non-selective breeding pigs. The rib eye areas of homozygous or heterozygous pigs with Lr-0.13 were larger than those of homozygous pigs with Lr-0.30 and non-selective breeding pigs. These data suggested that SLA haplotypes had the potential as useful genetic markers for selective breeding in the population of SLA-defined Duroc pigs.

Effect of swine leukocyte antigen-DQA gene variation on diarrhea in Large White, Landrace, and Duroc piglets

Piglet diarrhea is one of the most common factors that affects the benefits of the swine industry. Although recent studies have shown that exon 2 of SLA-DQA is associated with piglet resistance to diarrhea, contributions of genetic variation in the additional exon coding regions of this gene remain unclear. Here, we investigated variation in exons 1, 3 and 4 of the SLA-DQA gene and evaluated their effects on diarrheal infection in 425 suckling piglets. No variation was identified in exon 1. In exon 3, there were eight alleles detected, generated by 14 single nucleotide polymorphisms (SNPs) and three nucleotide deletions, eight SNPs being newly identified. Four allele sequences and three SNPs were identified in exon 4, only one SNP being newly identified. Statistical analysis showed that the genotypes of exon 3 are significantly associated with piglet diarrhea; indeed, genotypes DQA*wb01/wb02 and wb04/wb05 are clearly associated with resistance to piglet diarrhea, as they have the lowest probabilities of infection (P < 0.05). However, no significant association was found between the genotypes of exon 4 and diarrhea (P > 0.05). These results provide important new information concerning the level of genetic diversity at the SLA-DQA locus and suggest that further genetic association studies of piglet diarrhea resistance should include analyses of both exons 2 and 3 of this locus.

Identification of the Q969R gain-of-function polymorphism in the gene encoding porcine NLRP3 and its distribution in pigs of Asian and European origin

The nucleotide-binding domain, leucine-rich-containing family, pyrin-domain containing-3 (NLRP3) inflammasome comprises the major components caspase-1, apoptosis-associated speck-like protein containing a caspase recruitment domain (ASC), and NLRP3. NLRP3 plays important roles in maintaining immune homeostasis mediated by intestinal microorganisms and in the immunostimulatory properties of vaccine adjuvants used to induce an immune response. In the present study, we first cloned a complementary DNA (cDNA) encoding porcine ASC because its genomic sequence was not completely determined. The availability of the ASC cDNA enabled us to reconstitute porcine NLRP3 inflammasomes using an in vitro system that led to the identification of the immune functions of porcine NLRP3 and ASC based on the production of interleukin-1β (IL-1β). Further, we identified six synonymous and six nonsynonymous single-nucleotide polymorphisms (SNPs) in the coding sequence of NLRP3 of six breeds of pigs, including major commercial breeds. Among the nonsynonymous SNPs, the Q969R polymorphism is associated with an increased release of IL-1β compared with other porcine NLRP3 variants, indicating that this polymorphism represents a gain-of-function mutation. This allele was detected in 100 % of the analyzed Chinese Jinhua and Japanese wild boars, suggesting that the allele is maintained in the major commercial native European breeds Landrace, Large White, and Berkshire. These findings represent an important contribution to our knowledge of the diversity of NLRP3 nucleotide sequences among various pig populations. Moreover, efforts to exploit the gain of function induced by the Q969R polymorphism promise to improve pig breeding and husbandry by conferring enhanced resistance to pathogens as well as contributing to vaccine efficacy.

Identification of genetic variation in the swine toll-like receptors and development of a porcine TLR genotyping array

BACKGROUND

Toll-like receptors (TLR) are crucial in innate immunity for the recognition of a broad range of microbial pathogens and are expressed in multiple cell types. There are 10 TLR genes described in the pig genome.

RESULTS

With a twofold objective i.e. to catalogue genetic variants in porcine TLR genes and develop a genotyping array for genetic association studies on immune-related traits, we combined targeted sub-genome enrichment and high-throughput sequencing to sequence the 10 porcine TLR genes in 266 pigs from 10 breeds and wild boars using a DNA-pooling strategy. We identified 306 single nucleotide variants across the 10 TLR and 11 populations, 87 of which were novel. One hundred and forty-seven positions i.e. six stop-gains and 141 non-synonymous substitutions were predicted to alter the protein sequence. Three positions were unique to a single breed with alternative allele frequencies equal to or higher than 0.5. We designed a genotyping array for future applications in genetic association studies, with a selection of 126 variants based on their predicted impact on protein sequence. Since TLR4, TLR7 and TLR9 were underrepresented in this selection, we also included three variants that were located in the 3'UTR of these genes. We tested the array by genotyping 214 of the 266 sequenced pigs. We found that 93 variants that involved the 10 TLR genes were polymorphic in these animals. Twelve of these variants were novel. Furthermore, seven known variants that are associated with immune-related phenotypes are present on the array and can thus be used to test such associations in additional populations.

CONCLUSIONS

We identified genetic variations that potentially have an impact on the protein sequence of porcine TLR. A genotyping array with 80 non-synonymous, 10 synonymous and three 3'UTR polymorphisms in the 10 TLR genes is now available for association studies in swine populations with measures on immune-related traits.

Genome-wide identification of allele-specific expression in response to Streptococcus suis 2 infection in two differentially susceptible pig breeds

Although allele expression imbalance has been recognized in many species, and strongly linked to diseases, no whole transcriptome allele imbalance has been detected in pigs during pathogen infections. The pathogen Streptococcus suis 2 (SS2) causes serious zoonotic disease. Different pig breeds show differential susceptibility/resistance to pathogen infection, but the biological insight is little known. Here we analyzed allele-specific expression (ASE) using the spleen transcriptome of four pigs belonging to two phenotypically different breeds after SS2 infection. The comparative analysis of allele specific SNPs between control and infected animals revealed 882 and 1096 statistically significant differentially expressed allele SNPs (criteria: ratio ≧ 2 or ≦ 0.5) in Landrace and Enshi black pig, respectively. Twenty nine allelically imbalanced SNPs were further verified by Sanger sequencing, and later six SNPs were quantified by pyrosequencing assay. The pyrosequencing results are in agreement with the RNA-seq results, except two SNPs. Looking at the role of ASE in predisposition to diseases, the discovery of causative variants by ASE analysis might help the pig industry in long term to design breeding programs for improving SS2 resistance.

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.

Swine Leukocyte Antigen-DQA Gene Variation and Its Association with Piglet Diarrhea in Large White, Landrace and Duroc

The swine leukocyte antigen class II molecules are possibly associated with the induction of protective immunity. The study described here was to investigate the relationship between polymorphisms in exon 2 of the swine DQA gene and piglet diarrhea. This study was carried out on 425 suckling piglets from three purebred pig strains (Large White, Landrace and Duroc). The genetic diversity of exon 2 in swine DQA was detected by PCR-SSCP and sequencing analysis, eight unique SSCP patterns (AB, BB, BC, CC, CD, BD, BE and DD) representing five specific allele (A to E) sequences were detected. Sequence analysis revealed 21 nucleotide variable sites and resulting in 12 amino acid substitutions in the populations. A moderate level polymorphism and significant deviations from Hardy-Weinberg equilibrium of the genotypes distribution were observed in the populations (p<0.01). The association analysis indicated that there was a statistically significant difference in the score of piglet diarrhea between different genotypes, individuals with genotype CC showed a lower diarrhea score than genotypes AB (0.98±0.09), BB (0.85±0.77) and BC (1.25±0.23) (p<0.05), and significantly low than genotype BE (1.19±0.19) (p<0.01), CC genotype may be a most resistance genotype for piglet diarrhea.