Genetic analysis of reproductive performance in sows during porcine reproductive and respiratory syndrome (PRRS) and porcine epidemic diarrhea (PED) outbreaks

Elucidation of novel SNPs affecting immune response to classical swine fever vaccination in pigs using immunogenomics approach

The host genetic makeup plays a significant role in causing the within-breed variation among individuals after vaccination. The present study was undertaken to elucidate the genetic basis of differential immune response between high and low responder Landlly (Landrace X Ghurrah) piglets vis-à-vis CSF vaccination. For the purpose, E2 antibody response against CSF vaccination was estimated in sampled animals on the day of vaccination and 21-day post-vaccination as a measure of humoral immune response. Double-digestion restriction associated DNA (ddRAD) sequencing was undertaken on 96 randomly chosen Landlly piglets using Illumina HiSeq platform. SNP markers were called using standard methodology. Genome-wide association study (GWAS) was undertaken in PLINK program to identify the informative SNP markers significantly associated with differential immune response. The results revealed significant SNPs associated with E2 antibody response against CSF vaccination. The genome-wide informative SNPs for the humoral immune response against CSF vaccination were located on SSC10, SSC17, SSC9, SSC2, SSC3 and SSC6. The overlapping and flanking genes (500Kb upstream and downstream) of significant SNPs were CYB5R1, PCMTD2, WT1, IL9R, CD101, TMEM64, TLR6, PIGG, ADIPOR1, PRSS37, EIF3M, and DNAJC24. Functional enrichment and annotation analysis were undertaken for these genes in order to gain maximum insights into the association of these genes with immune system functionality in pigs. The genetic makeup was associated with differential immune response against CSF vaccination in Landlly piglets while the identified informative SNPs may be used as suitable markers for determining variation in host immune response against CSF vaccination in pigs.

Gut bacterial species in late trimester of pregnant sows influence the occurrence of stillborn piglet through pro-inflammation response

Maternal gut microbiota is an important regulator for the metabolism and immunity of the fetus during pregnancy. Recent studies have indicated that maternal intestinal microbiota is closely linked to the development of fetus and infant health. Some bacterial metabolites are considered to be directly involved in immunoregulation of fetus during pregnancy. However, the detailed mechanisms are largely unknown. In this study, we exploited the potential correlation between the gut microbiota of pregnant sows and the occurrence of stillborn piglets by combining the 16S rRNA gene and metagenomic sequencing data, and fecal metabolome in different cohorts. The results showed that several bacterial species from Bacteroides, potential pathogens, and LPS-producing bacteria exhibited significantly higher abundances in the gut of sows giving birth to stillborn piglets. Especially, Bacteroides fragilis stood out as the key driver in both tested cohorts and showed the most significant association with the occurrence of stillborn piglets in the DN1 cohort. However, several species producing short-chain fatty acids (SCFAs), such as Prevotella copri, Clostridium butyricum and Faecalibacterium prausnitzii were enriched in the gut of normal sows. Functional capacity analysis of gut microbiome revealed that the pathways associated with infectious diseases and immune diseases were enriched in sows giving birth to stillborn piglets. However, energy metabolism had higher abundance in normal sows. Fecal metabolome profiling analysis found that Lysophosphatidylethanolamine and phosphatidylethanolamine which are the main components of cell membrane of Gram-negative bacteria showed significantly higher concentration in stillbirth sows, while SCFAs had higher concentration in normal sows. These metabolites were significantly associated with the stillborn-associated bacterial species including Bacteroides fragilis. Lipopolysaccharide (LPS), IL-1β, IL-6, FABP2, and zonulin had higher concentration in the serum of stillbirth sows, indicating increased intestinal permeability and pro-inflammatory response. The results from this study suggested that certain sow gut bacterial species in late trimester of pregnancy, e.g., an excess abundance of Bacteroides fragilis, produced high concentration of LPS which induced sow pro-inflammatory response and might cause the death of the relatively weak piglets in a farrow. This study provided novel evidences about the effect of maternal gut microbiota on the fetus development and health.

Host Genetics of Response to Porcine Reproductive and Respiratory Syndrome in Sows: Reproductive Performance

Porcine Reproductive and Respiratory Syndrome (PRRS) is historically the most economically important swine disease worldwide that severely affects the reproductive performance of sows. However, little is still known about the genetic basis of reproductive performance in purebred herds during a PRRS outbreak through the comparison of maternal and terminal breeds. Thus, the objective of this work was to explore the host genetics of response to PRRS in purebred sows from two breeds. Reproductive data included 2546 Duroc and 2522 Landrace litters from 894 and 813 purebred sows, respectively, which had high-density genotype data available (29,799 single nucleotide polymorphisms; SNPs). The data were split into pre-PRRS, PRRS, and post-PRRS phases based on standardized farrow-year-week estimates. Heritability estimates for reproductive traits were low to moderate (≤0.20) for Duroc and Landrace across PRRS phases. On the other hand, genetic correlations of reproductive traits between PRRS phases were overall moderate to high for both breeds. Several associations between MARC0034894, a candidate SNP for response to PRRS, with reproductive performance were identified (P-value < 0.05). Genomic analyses detected few QTL for reproductive performance across all phases, most explaining a small percentage of the additive genetic variance (≤8.2%, averaging 2.1%), indicating that these traits are highly polygenic. None of the identified QTL within a breed and trait overlapped between PRRS phases. Overall, our results indicate that Duroc sows are phenotypically more resilient to PRRS than Landrace sows, with a similar return to PRRS-free performance between breeds for most reproductive traits. Genomic prediction results indicate that genomic selection for improved reproductive performance under a PRRS outbreak is possible, especially in Landrace sows, by training markers using data from PRRS-challenged sows. On the other hand, the high genetic correlations with reproductive traits between PRRS phases suggest that selection for improved reproductive performance in a clean environment could improve performance during PRRS, but with limited efficiency due to their low heritability estimates. Thus, we hypothesize that an indicator trait that could be indirectly selected to increase the response to selection for these traits would be desirable and would also improve the reproductive performance of sows during a PRRS outbreak.

Host Genetics of Response to Porcine Reproductive and Respiratory Syndrome in Sows: Antibody Response as an Indicator Trait for Improved Reproductive Performance

Antibody response to porcine reproductive and respiratory syndrome (PRRS) virus (PRRSV) infection, measured as sample-to-positive (S/P) ratio, has been proposed as an indicator trait for improved reproductive performance during a PRRS outbreak in Landrace sows. However, this result has not yet been validated in Landrace sows or evaluated in terminal sire lines. The main objectives of this work were to validate the use of S/P ratio as an indicator trait to select pigs during a PRRS outbreak and to explore the genetic basis of antibody response to PRRSV. Farrowing data included 2,546 and 2,522 litters from 894 Duroc and 813 Landrace sows, respectively, split into pre-PRRS, PRRS, and post-PRRS phases. Blood samples were taken from 1,231 purebred sows (541 Landrace and 690 Duroc) following a PRRS outbreak for subsequent PRRSV ELISA analysis for S/P ratio measurement. All animals had high-density genotype data available (29,799 single nucleotide polymorphisms; SNPs). Genetic parameters and genome-wide association studies (GWAS) for S/P ratio were performed for each breed separately. Heritability estimates (± standard error) of S/P ratio during the PRRS outbreak were moderate, with 0.35 ± 0.08 for Duroc and 0.34 ± 0.09 for Landrace. During the PRRS outbreak, favorable genetic correlations of S/P ratio with the number of piglets born alive (0.61 ± 0.34), number of piglets born dead (-0.33 ± 0.32), and number of stillborn piglets (-0.27 ± 0.31) were observed for Landrace sows. For Duroc, the GWAS identified a major quantitative trait locus (QTL) on chromosome (Chr) 7 (24-15 megabases; Mb) explaining 15% of the total genetic variance accounted for by markers (TGVM), and another one on Chr 8 (25 Mb) explaining 2.4% of TGVM. For Landrace, QTL on Chr 7 (24-25 Mb) and Chr 7 (108-109 Mb), explaining 31% and 2.2% of TGVM, respectively, were identified. Some of the SNPs identified in these regions for S/P ratio were associated with reproductive performance but not during the PRRS outbreak. Genomic prediction accuracies for S/P ratio were moderate to high for the within-breed analysis. For the between-breed analysis, these were overall low. These results further support the use of S/P ratio as an indicator trait for improved reproductive performance during a PRRS outbreak in Landrace sows.

Genomics of response to porcine reproductive and respiratory syndrome virus in purebred and crossbred sows: antibody response and performance following natural infection vs. vaccination

Antibody response, measured as sample-to-positive (S/P) ratio, to porcine reproductive and respiratory syndrome virus (PRRSV) following a PRRSV-outbreak (S/POutbreak) in a purebred nucleus and following a PRRSV-vaccination (S/PVx) in commercial crossbred herds have been proposed as genetic indicator traits for improved reproductive performance in PRRSV-infected purebred and PRRSV-vaccinated crossbred sows, respectively. In this study, we investigated the genetic relationships of S/POutbreak and S/PVx with performance at the commercial (vaccinated crossbred sows) and nucleus level (non-infected and PRRSV-infected purebred sows), respectively, and tested the effect of previously identified SNP for these indicator traits. Antibody response was measured on 541 Landrace sows ~54 d after the start of a PRRSV outbreak, and on 906 F1 (Landrace × Large White) gilts ~50 d after vaccination with a commercial PRRSV vaccine. Reproductive performance was recorded for 711 and 428 Landrace sows before and during the PRRSV outbreak, respectively, and for 811 vaccinated F1 animals. The estimate of the genetic correlation (rg) of S/POutbreak with S/PVx was 0.72 ± 0.18. The estimates of rg of S/POutbreak with reproductive performance in vaccinated crossbred sows were low to moderate, ranging from 0.05 ± 0.23 to 0.30 ± 0.20. The estimate of rg of S/PVx with reproductive performance in non-infected purebred sows was moderate and favorable with number born alive (0.50 ± 0.23) but low (0 ± 0.23 to -0.11 ± 0.23) with piglet mortality traits. The estimates of rg of S/PVx were moderate and negative (-0.38 ± 0.21) with number of mummies in PRRSV-infected purebred sows and low with other traits (-0.30 ± 0.18 to 0.05 ± 0.18). Several significant associations (P0 > 0.90) of previously reported SNP for S/P ratio (ASGA0032063 and H3GA0020505) were identified for S/P ratio and performance in non-infected purebred and PRRSV-exposed purebred and crossbred sows. Genomic regions harboring the major histocompatibility complex class II region significantly contributed to the genetic correlation of antibody response to PRRSV with most of the traits analyzed. These results indicate that selection for antibody response in purebred sows following a PRRSV outbreak in the nucleus and for antibody response to PRRSV vaccination measured in commercial crossbred sows are expected to increase litter size in purebred and commercial sows.

Assessing the statistical training in animal science graduate programs in the United States: survey on statistical training

Statistical analysis of data and understanding of experimental design are critical skills needed by animal science graduate students (ASGS). These skills are even more valuable with the increased development of high-throughput technologies. The objective of this study was to evaluate the perceived statistical training of U.S. ASGS. A survey with 38 questions was shared across U.S. universities, and 416 eligible ASGS from 43 universities participated in this study. The survey included questions on the demographics and overall training, graduate education on statistics, and self-assessment on statistics and career path of ASGS. Several analyses were performed: relationship between perceived received education (PRE; i.e., how ASGS evaluated their graduate education in statistics) and perceived knowledge (PK; i.e., how ASGS evaluated their knowledge in statistics from their education); ranking of statistical topics based on PRE, PK, and confidence in performing statistical analyses (CPSA); cluster analysis of statistical topics for PRE, PK, and CPSA; and factors (demographic, overall training, interest in statistics, and field of study) associated with the overall scores (OS) for PRE, PK, and CPSA. Students had greater (P < 0.05) PRE than PK for most of the statistical topics included in this study. The moderate to high repeatability of answers within statistical topics indicates substantial correlations in ASGS answers between PRE and PK. The cluster analysis resulted in distinct groups of "Traditional" and "Nontraditional" statistical topics. ASGS showed lower (P < 0.05) scores of PRE, PK, and CPSA in "Nontraditional" compared with "Traditional" statistical methods. Several factors were associated (P < 0.05) with the OS of PRE, PK, and CSPA. In general, factors related to greater training and interest in statistics of ASGS were associated with greater OS, such as taking more credits in statistics courses, having additional training in statistics outside the classroom, knowing more than one statistics software, and more. This study provided comprehensive information on the perceived level of education, knowledge, and confidence in statistics in ASGS in the United States. Although objective measurements of their training in statistics are needed, the current study suggests that ASGS have limited statistical training on topics of major importance for the current and future trends of data-driven research in animal sciences.

A descriptive survey of porcine epidemic diarrhea in pig populations in northern Vietnam

Porcine epidemic diarrhea (PED) virus (PEDV) is a globally emerging and re-emerging epizootic swine virus that causes massive economic losses in the swine industry, with high mortality in piglets. In Vietnam, PED first emerged in 2009 and has now developed to an endemic stage. This is the first cross-sectional survey performed to evaluate the proportion of PEDV-positive swine farms in Vietnam from January 2018 to February 2019. Fecal samples from 327 pig farms in northern Vietnam were collected and tested for PEDV infection by reverse transcription-loop-mediated isothermal amplification (RT-LAMP) method. The proportion of PEDV-positive farms was 30.9% and PEDV-positive farms were distributed throughout the study area. The highest proportion of PEDV-positive farms was 70% (7/10) among nucleus production type farms (P < 0.05). Higher proportions of PEDV-positive farms were found in the Northeast and Red River Delta areas, which are the major areas of pig production (P < 0.05). The proportion of PEDV-positive farms was higher among larger farms (P < 0.05). Our findings illustrate the high proportion of PEDV-positive farms in the Vietnamese pig population and will help to better understand the epidemiological dynamics of PED infection, to estimate impact, and establish and improve prevention and control measures.

Evaluation of a commercial genetic test for fescue toxicosis in pregnant Angus beef cattle

Most tall fescue [Lolium arundinaceum (Schreb.) Darbysh] in the Southeastern United States contains an endophyte that causes fescue toxicosis (FT) in grazing animals, a serious disease that causes approximately $1 billion in economic losses to the animal industries in the United States. Recently, a genetic test called T-Snip (AgBotanica, LCC, Columbia, MO), was developed with the objective of identifying animals with genetic variation for FT tolerance. The aim of this study was to validate the use of this genetic test in mature, pregnant cows. Over 13 wk, weekly phenotypic data, including body weight, rectal temperatures, hair coat scores, hair shedding scores, and body condition scores, were collected on 148 pregnant purebred Angus cows at 2 locations in NC where infected fescue was the primary source of feed. Birth weights (cBW) and 205-d adjusted weaning weights (adjWW) from these cow’s calves were recorded. All cows were genotyped for T-Snip. At the end of the trial, each phenotypic trait was calculated as the slope of the linear regression of performance on weeks. The effect of T-Snip rating genotypes (4 levels) on slope traits was tested using a linear model also including the fixed-effects of location, parity, and the initial measurement for each trait (covariate). For cBW and adjWW, the model also included the sex of the calf and the month of birth as categorical effects. Associations of T-Snip genotypes were observed for body weight gain (aBWd) of pregnant cows (P = 0.15; interaction with location), change in body condition score (aBCSd; P = 0.13), and adjWW (P = 0.06; interaction with location). For aBWd and adjWW, associations were found just within one location (P = 0.017 and 0.047, respectively), which was the location with higher endophyte infection rate. For all associations, the direction of the T-Snip genotypes was the same and as expected: the greater the genotype score, the better performance. No associations were found for the other traits (P > 0.10). These results indicate that the T-Snip test may be predictive of cow performance (aBWd, aBCSd, and adjWW) in an endophyte-infected tall fescue environment.

Integrative analysis of microRNA and mRNA expression profiles in MARC-145 cells infected with PRRSV

Porcine reproductive and respiratory syndrome virus (PRRSV) causes tremendous economic losses to the swine industry worldwide. miRNAs are crucial regulators of gene expression and a wide range of complex interactions of miRNAs-mRNAs is possible during virus infection. However, there is no comprehensive integrated study of miRNA and mRNA networks in MARC-145 cells after infection with PRRSV. We analyzed the differential expressions, co-relations, annotations, and putative functions of miRNA and mRNA networks in PRRSV-infected MARC-145 cells. Based on the filtering criterion, 22 differentially expressed miRNAs (DEmiRs) (15 up- and 7 downregulated) were filtered out. miRNA-mRNA interaction networks were constructed. For the 18 selected miRNAs, 390 potential target genes were predicted from the differentially expressed mRNAs (DEmRs). GO and KEGG pathway annotations predicted 34 KEGG pathways, 12 of which are known to be involved in virus infection. Real-time PCR validated the RNA-seq results. Our analysis showed that miR-27a-5p and miR-21-3p downregulate the expression of two of their potential target genes-SPARC, CLIC1, and cofilin-1, COX7A2, respectively. Further experiments proved that miR-21-3p and miR-27a-5p can promote PRRSV replication significantly. It is the first report that these two miRNAs participate in the interaction of host cells with PRRSV. Our results provide insights into the role of miRNAs in response to PRRSV infection, which will aid the research for developing novel therapies against PRRSV.

Investigating the relationship between vaginal microbiota and host genetics and their impact on immune response and farrowing traits in commercial gilts

Our objectives were to evaluate the interaction between host genetics and vaginal microbiota and their relationships with antibody (Ab) response to porcine reproductive and respiratory syndrome virus (PRRSV) vaccination and farrowing performance in commercial gilts. The farrowing performance traits were number born alive, number weaning (NW), total number born, number born dead, stillborn, mummies and preweaning mortality (PWM). The vaginal microbiota was collected on days 4 (D4) and 52 (D52) after vaccination for PRRSV. Blood samples were collected on D52 for Ab measurement. Actinobacteria, Bacterioidetes, Firmicutes, Proteobacteria and Tenericutes were the most abundant Phyla identified in the vaginal microbiota. Heritability ranged from ~0 to 0.60 (Fusobacterium) on D4 and from ~0 to 0.63 (Terrisporobacter) on D52, with 43 operational taxonomic units (OTUs) presenting moderate to high heritability. One major QTL on chromosome 12 was identified for 5 OTUs (Clostridiales, Acinetobacter, Ruminococcaceae, Campylobacter and Anaerococcus), among other 19 QTL. The microbiability for Ab response to PRRSV vaccination was low for both days (<0.07). For farrowing performance, microbiability varied from <0.001 to 0.15 (NW on D4). For NW and PWM, the microbiability was greater than the heritability estimates. Actinobacillus, Streptococcus, Campylobacter, Anaerococcus, Mollicutes, Peptostreptococcus, Treponema and Fusobacterium showed different abundance between low and high Ab responders. Finally, canonical discriminant analyses revealed that vaginal microbiota was able to classify gilts in high and low Ab responders to PRRSV vaccination with a misclassification rate of <0.02. Although the microbiota explained limited variation in Ab response and farrowing performance traits, there is still potential to explore the use of vaginal microbiota to explain variation in traits such as NW and PWM. In addition, these results revealed that there is a partial control of host genetic over vaginal microbiota, suggesting a possibility for genetic selection on the vaginal microbiota.

Genetic Markers Associated with Field PRRSV-Induced Abortion Rates

In gilts and sows, the more severe clinical manifestation of porcine reproductive and respiratory syndrome virus (PRRSV) occurs in late gestation and can result in up to a 40% abortion incidence. Despite the known genetic component in resilience to PRRSV, there is scarce information regarding the abortive outcome of this disease. We tested the relationship between eight molecular markers (six from published studies and two identified in the present study in the HDAC6 gene) and the probability of abortion during a PRRSV outbreak, using data from two commercial Landrace x Large White sow farms with an incidence of abortion of 35% and 17%. From the markers tested, USP18_-1533G>A did not segregate in these populations, and CD163_c.3534C>T and HDAC6_g.2360C>T did not affect the abortion rate. In contrast, the minor allele of two markers in SSC4 (WUR1000125 in GBP1 and rs340943904 in GBP5), which lower viremia in growing pigs, and the major alleles of CD163_rs1107556229 and HDAC6_rs325981825 were associated with a lower probability of abortion during PRRSV outbreaks. The more striking result was for the MX1 gene, where the odds ratio of aborting versus not aborting was nine times lower in the sows homozygous for a 275-bp insertion than in the other genotypes. Interactions between markers were not relevant. All together, we bring here the first evidence that mutations in the host genome can predispose or protect from complete reproductive failure in sows infected with PRRSV.