A DIO2 missense mutation and its impact on fetal response to PRRSV infection.

BACKGROUND: Porcine reproductive and respiratory syndrome virus 2 (PRRSV-2) infection during late gestation substantially lowers fetal viability and survival. In a previous genome-wide association study, a single nucleotide polymorphism on chromosome 7 was significantly associated with probability of fetuses being viable in response to maternal PRRSV-2 infection at 21 days post maternal inoculation. The iodothyronine deiodinase 2 (DIO2) gene, located ~ 14 Kilobase downstream of this SNP, was selected as a priority candidate related to fetal susceptibility following maternal PRRSV-2 infection. Our objectives were to identify mutation(s) within the porcine DIO2 gene and to determine if they were associated with fetal outcomes after PRRSV-2 challenge. Sequencing of the DIO2, genotyping identified variants, and association of DIO2 genotypes with fetal phenotypes including DIO2 mRNA levels, viability, survival, viral loads, cortisol and thyroid hormone levels, and growth measurements were conducted. RESULTS: A missense variant (p.Asn91Ser) was identified in the parental populations from two independent PRRSV-2 challenge trials. This variant was further genotyped to determine association with fetal PRRS outcomes. DIO2 mRNA levels in fetal heart and kidney differed by the genotypes of Asn91Ser substitution with significantly greater DIO2 mRNA expression in heterozygotes compared with wild-type homozygotes (P < 0.001 for heart, P = 0.002 for kidney). While Asn91Ser did not significantly alter fetal viability and growth measurements, interaction effects of the variant with fetal sex or trial were identified for fetal viability or crown rump length, respectively. However, this mutation was not related to dysregulation of the hypothalamic-pituitary-adrenal and thyroid axis, indicated by no differences in circulating cortisol, T4, and T3 levels in fetuses of the opposing genotypes following PRRSV-2 infection. CONCLUSIONS: The present study suggests that a complex relationship among DIO2 genotype, DIO2 expression, fetal sex, and fetal viability may exist during the course of fetal PRRSV infection. Our study also proposes the increase in cortisol levels, indicative of fetal stress response, may lead to fetal complications, such as fetal compromise, fetal death, or premature farrowing, during PRRSV infection.

Genome-wide association studies for additive and dominance effects for body composition traits in commercial crossbred Piétrain pigs.

Fat depth (FD) and muscle depth (MD) are economically important traits and used to estimate carcass lean content (LMP), which is one of the main breeding objectives in pig breeding programmes. We assessed the genetic architectures of body composition traits for additive and dominance effects in commercial crossbred Piétrain pigs using both 50 K array and sequence genotypes. We first performed a genome-wide association study (GWAS) using single-marker association analysis with a false discovery rate of 0.1. Then, we estimated the additive and dominance effects of the most significant variant in the quantitative trait loci (QTL) regions. It was investigated whether the use of whole-genome sequence (WGS) will improve the QTL detection (both additive and dominance) with a higher power compared with lower density SNP arrays. Our results showed that more QTL regions were detected by WGS compared with 50 K array (n = 54 vs. n = 17). Of the novel associated regions associated with FD and LMP and detected by WGS, the most pronounced peak was on SSC13, situated at ~116-118, 121-127 and 129-134 Mbp. Additionally, we found that only additive effects contributed to the genetic architecture of the analysed traits and no significant dominance effects were found for the tested SNPs at QTL regions, regardless of panel density. The associated SNPs are located in or near several relevant candidate genes. Of these genes, GABRR2, GALR1, RNGTT, CDH20 and MC4R have been previously reported as being associated with fat deposition traits. However, the genes on SSC1 (ZNF292, ORC3, CNR1, SRSF12, MDN1, TSHZ1, RELCH and RNF152) and SSC18 (TTC26 and KIAA1549) have not been reported previously to our best knowledge. Our current findings provide insights into the genomic regions influencing composition traits in Piétrain pigs.

Genetic analysis of disease resilience of wean-to-finish pigs under a natural disease challenge model using reaction norms.

BACKGROUND: Disease resilience is the ability to maintain performance across environments with different disease challenge loads (CL). A reaction norm describes the phenotypes that a genotype can produce across a range of environments and can be implemented using random regression models. The objectives of this study were to: (1) develop measures of CL using growth rate and clinical disease data recorded under a natural polymicrobial disease challenge model; and (2) quantify genetic variation in disease resilience using reaction norm models. METHODS: Different CL were derived from contemporary group effect estimates for average daily gain (ADG) and clinical disease phenotypes, including medical treatment rate (TRT), mortality rate, and subjective health scores. Resulting CL were then used as environmental covariates in reaction norm analyses of ADG and TRT in the challenge nursery and finisher, and compared using model loglikelihoods and estimates of genetic variance associated with CL. Linear and cubic spline reaction norm models were compared based on goodness-of-fit and with multi-variate analyses, for which phenotypes were separated into three traits based on low, medium, or high CL. RESULTS: Based on model likelihoods and estimates of genetic variance explained by the reaction norm, the best CL for ADG in the nursery was based on early ADG in the finisher, while the CL derived from clinical disease traits across the nursery and finisher was best for ADG in the finisher and for TRT in the nursery and across the nursery and finisher. With increasing CL, estimates of heritability for nursery and finisher ADG initially decreased, then increased, while estimates for TRT generally increased with CL. Genetic correlations for ADG and TRT were low between high versus low CL, but high for close CL. Linear reaction norm models fitted the data significantly better than the standard genetic model without genetic slopes, while the cubic spline model fitted the data significantly better than the linear reaction norm model for most traits. Reaction norm models also fitted the data better than multi-variate models. CONCLUSIONS: Reaction norm models identified genotype-by-environment interactions related to disease CL. Results can be used to select more resilient animals across different levels of CL, high-performance animals at a given CL, or a combination of these.

Identification of candidate genes and enriched biological functions for feed efficiency traits by integrating plasma metabolites and imputed whole genome sequence variants in beef cattle.

BACKGROUND: Feed efficiency is one of the key determinants of beef industry profitability and sustainability. However, the cellular and molecular background behind feed efficiency is largely unknown. This study combines imputed whole genome DNA variants and 31 plasma metabolites to dissect genes and biological functions/processes that are associated with residual feed intake (RFI) and its component traits including daily dry matter intake (DMI), average daily gain (ADG), and metabolic body weight (MWT) in beef cattle. RESULTS: Regression analyses between feed efficiency traits and plasma metabolites in a population of 493 crossbred beef cattle identified 5 (L-valine, lysine, L-tyrosine, L-isoleucine, and L-leucine), 4 (lysine, L-lactic acid, L-tyrosine, and choline), 1 (citric acid), and 4 (L-glutamine, glycine, citric acid, and dimethyl sulfone) plasma metabolites associated with RFI, DMI, ADG, and MWT (P-value < 0.1), respectively. Combining the results of metabolome-genome wide association studies using 10,488,742 imputed SNPs, 40, 66, 15, and 40 unique candidate genes were identified as associated with RFI, DMI, ADG, and MWT (P-value < 1 × 10-5), respectively. These candidate genes were found to be involved in some key metabolic processes including metabolism of lipids, molecular transportation, cellular function and maintenance, cell morphology and biochemistry of small molecules. CONCLUSIONS: This study identified metabolites, candidate genes and enriched biological functions/processes associated with RFI and its component traits through the integrative analyses of metabolites with phenotypic traits and DNA variants. Our findings could enhance the understanding of biochemical mechanisms of feed efficiency traits and could lead to improvement of genomic prediction accuracy via incorporating metabolite data.

Investigating the genetic architecture of disease resilience in pigs by genome-wide association studies of complete blood count traits collected from a natural disease challenge model.

BACKGROUND: Genetic improvement for disease resilience is anticipated to be a practical method to improve efficiency and profitability of the pig industry, as resilient pigs maintain a relatively undepressed level of performance in the face of infection. However, multiple biological functions are known to be involved in disease resilience and this complexity means that the genetic architecture of disease resilience remains largely unknown. Here, we conducted genome-wide association studies (GWAS) of 465,910 autosomal SNPs for complete blood count (CBC) traits that are important in an animal's disease response. The aim was to identify the genetic control of disease resilience. RESULTS: Univariate and multivariate single-step GWAS were performed on 15 CBC traits measured from the blood samples of 2743 crossbred (Landrace × Yorkshire) barrows drawn at 2-weeks before, and at 2 and 6-weeks after exposure to a polymicrobial infectious challenge. Overall, at a genome-wise false discovery rate of 0.05, five genomic regions located on Sus scrofa chromosome (SSC) 2, SSC4, SSC9, SSC10, and SSC12, were significantly associated with white blood cell traits in response to the polymicrobial challenge, and nine genomic regions on multiple chromosomes (SSC1, SSC4, SSC5, SSC6, SSC8, SSC9, SSC11, SSC12, SSC17) were significantly associated with red blood cell and platelet traits collected before and after exposure to the challenge. By functional enrichment analyses using Ingenuity Pathway Analysis (IPA) and literature review of previous CBC studies, candidate genes located nearby significant single-nucleotide polymorphisms were found to be involved in immune response, hematopoiesis, red blood cell morphology, and platelet aggregation. CONCLUSIONS: This study helps to improve our understanding of the genetic basis of CBC traits collected before and after exposure to a polymicrobial infectious challenge and provides a step forward to improve disease resilience.

Quantitative analysis of the blood transcriptome of young healthy pigs and its relationship with subsequent disease resilience.

BACKGROUND: Disease resilience, which is the ability of an animal to maintain performance under disease, is important for pigs in commercial herds, where they are exposed to various pathogens. Our objective was to investigate population-level gene expression profiles in the blood of 912 healthy F1 barrows at ~ 27 days of age for associations with performance and health before and after their exposure to a natural polymicrobial disease challenge at ~ 43 days of age. RESULTS: Most significant (q < 0.20) associations of the level of expression of individual genes in blood of young healthy pigs were identified for concurrent growth rate and subjective health scores prior to the challenge, and for mortality, a combined mortality-treatment trait, and feed conversion rate after the challenge. Gene set enrichment analyses revealed three groups of gene ontology biological process terms that were related to disease resilience: 1) immune and stress response-related terms were enriched among genes whose increased expression was unfavorably associated with both pre- and post-challenge traits, 2) heme-related terms were enriched among genes that had favorable associations with both pre- and post-challenge traits, and 3) terms related to protein localization and viral gene expression were enriched among genes that were associated with reduced performance and health traits after but not before the challenge. CONCLUSIONS: Gene expression profiles in blood from young healthy piglets provide insight into their performance when exposed to disease and other stressors. The expression of genes involved in stress response, heme metabolism, and baseline expression of host genes related to virus propagation were found to be associated with host response to disease.

Synaptogyrin-2 influences replication of Porcine circovirus 2.

Porcine circovirus 2 (PCV2) is a circular single-stranded DNA virus responsible for a group of diseases collectively known as PCV2 Associated Diseases (PCVAD). Variation in the incidence and severity of PCVAD exists between pigs suggesting a host genetic component involved in pathogenesis. A large-scale genome-wide association study of experimentally infected pigs (n = 974), provided evidence of a host genetic role in PCV2 viremia, immune response and growth during challenge. Host genotype explained 64% of the phenotypic variation for overall viral load, with two major Quantitative Trait Loci (QTL) identified on chromosome 7 (SSC7) near the swine leukocyte antigen complex class II locus and on the proximal end of chromosome 12 (SSC12). The SNP having the strongest association, ALGA0110477 (SSC12), explained 9.3% of the genetic and 6.2% of the phenotypic variance for viral load. Dissection of the SSC12 QTL based on gene annotation, genomic and RNA-sequencing, suggested that a missense mutation in the SYNGR2 (SYNGR2 p.Arg63Cys) gene is potentially responsible for the variation in viremia. This polymorphism, located within a protein domain conserved across mammals, results in an amino acid variant SYNGR2 p.63Cys only observed in swine. PCV2 titer in PK15 cells decreased when the expression of SYNGR2 was silenced by specific-siRNA, indicating a role of SYNGR2 in viral replication. Additionally, a PK15 edited clone generated by CRISPR-Cas9, carrying a partial deletion of the second exon that harbors a key domain and the SYNGR2 p.Arg63Cys, was associated with a lower viral titer compared to wildtype PK15 cells (>24 hpi) and supernatant (>48hpi)(P < 0.05). Identification of a non-conservative substitution in this key domain of SYNGR2 suggests that the SYNGR2 p.Arg63Cys variant may underlie the observed genetic effect on viral load.

The genetic basis of natural antibody titers of young healthy pigs and relationships with disease resilience.

BACKGROUND: Disease resilience is the ability to maintain performance under pathogen exposure but is difficult to select for because breeding populations are raised under high health. Selection for resilience requires a trait that is heritable, easy to measure on healthy animals, and genetically correlated with resilience. Natural antibodies (NAb) are important parts of the innate immune system and are found to be heritable and associated with disease susceptibility in dairy cattle and poultry. Our objective was to investigate NAb and total IgG in blood of healthy, young pigs as potential indicator traits for disease resilience. RESULTS: Data were from Yorkshire x Landrace pigs, with IgG and IgM NAb (four antigens) and total IgG measured by ELISA in blood plasma collected ~ 1 week after weaning, prior to their exposure to a natural polymicrobial challenge. Heritability estimates were lower for IgG NAb (0.12 to 0.24, + 0.05) and for total IgG (0.19 + 0.05) than for IgM NAb (0.33 to 0.53, + 0.07) but maternal effects were larger for IgG NAb (0.41 to 0.52, + 0.03) and for total IgG (0.19 + 0.05) than for IgM NAb (0.00 to 0.10, + 0.04). Phenotypically, IgM NAb titers were moderately correlated with each other (average 0.60), as were IgG NAb titers (average 0.42), but correlations between IgM and IgG NAb titers were weak (average 0.09). Phenotypic correlations of total IgG were moderate with NAb IgG (average 0.46) but weak with NAb IgM (average 0.01). Estimates of genetic correlations among NAb showed similar patterns but with small SE, with estimates averaging 0.76 among IgG NAb, 0.63 among IgM NAb, 0.17 between IgG and IgM NAb, 0.64 between total IgG and IgG NAb, and 0.13 between total IgG and IgM NAb. Phenotypically, pigs that survived had slightly higher levels of NAb and total IgG than pigs that died. Genetically, higher levels of NAb tended to be associated with greater disease resilience based on lower mortality and fewer parenteral antibiotic treatments. Genome-wide association analyses for NAb titers identified several genomic regions, with several candidate genes for immune response. CONCLUSIONS: Levels of NAb in blood of healthy young piglets are heritable and potential genetic indicators of resilience to polymicrobial disease.

The effects of a globin blocker on the resolution of 3'mRNA sequencing data in porcine blood.

BACKGROUND: Gene expression profiling in blood is a potential source of biomarkers to evaluate or predict phenotypic differences between pigs but is expensive and inefficient because of the high abundance of globin mRNA in porcine blood. These limitations can be overcome by the use of QuantSeq 3'mRNA sequencing (QuantSeq) combined with a method to deplete or block the processing of globin mRNA prior to or during library construction. Here, we validated the effectiveness of QuantSeq using a novel specific globin blocker (GB) that is included in the library preparation step of QuantSeq. RESULTS: In data set 1, four concentrations of the GB were applied to RNA samples from two pigs. The GB significantly reduced the proportion of globin reads compared to non-GB (NGB) samples (P = 0.005) and increased the number of detectable non-globin genes. The highest evaluated concentration (C1) of the GB resulted in the largest reduction of globin reads compared to the NGB (from 56.4 to 10.1%). The second highest concentration C2, which showed very similar globin depletion rates (12%) as C1 but a better correlation of the expression of non-globin genes between NGB and GB (r = 0.98), allowed the expression of an additional 1295 non-globin genes to be detected, although 40 genes that were detected in the NGB sample (at a low level) were not present in the GB library. Concentration C2 was applied in the rest of the study. In data set 2, the distribution of the percentage of globin reads for NGB (n = 184) and GB (n = 189) samples clearly showed the effects of the GB on reducing globin reads, in particular for HBB, similar to results from data set 1. Data set 3 (n = 84) revealed that the proportion of globin reads that remained in GB samples was significantly and positively correlated with the reticulocyte count in the original blood sample (P < 0.001). CONCLUSIONS: The effect of the GB on reducing the proportion of globin reads in porcine blood QuantSeq was demonstrated in three data sets. In addition to increasing the efficiency of sequencing non-globin mRNA, the GB for QuantSeq has an advantage that it does not require an additional step prior to or during library creation. Therefore, the GB is a useful tool in the quantification of whole gene expression profiles in porcine blood.

Genome-wide association scan for heterotic quantitative trait loci in multi-breed and crossbred beef cattle.

BACKGROUND: Heterosis has been suggested to be caused by dominance effects. We performed a joint genome-wide association analysis (GWAS) using data from multi-breed and crossbred beef cattle to identify single nucleotide polymorphisms (SNPs) with significant dominance effects associated with variation in growth and carcass traits and to understand the mode of action of these associations. METHODS: Illumina BovineSNP50 genotypes and phenotypes for 11 growth and carcass traits were available for 6796 multi-breed and crossbred beef cattle. After performing quality control, 42,610 SNPs and 6794 animals were used for further analyses. A single-SNP GWAS for the joint association of additive and dominance effects was conducted in purebred, crossbred, and combined datasets using the ASReml software. Genomic breed composition predicted from admixture analyses was included in the mixed effect model to account for possible population stratification and breed effects. A threshold of 10% genome-wide false discovery rate was applied to declare associations as significant. The significant SNPs with dominance association were mapped to their corresponding genes at 100 kb. RESULTS: Seven SNPs with significant dominance associations were detected for birth weight, weaning weight, pre-weaning daily gain, yearling weight and marbling score across the three datasets at a false discovery rate of 10%. These SNPs were located on bovine chromosomes 1, 3, 4, 6 and 21 and mapped to six putative candidate genes: U6atac, AGBL4, bta-mir-2888-1, REPIN1, ICA1 and NXPH1. These genes have interesting biological functions related to the regulation of gene expression, glucose and lipid metabolism and body fat mass. For most of the identified loci, we observed over-dominance association with the studied traits, such that the heterozygous individuals at any of these loci had greater genotypic values for the trait than either of the homozygous individuals. CONCLUSIONS: Our results revealed very few regions with significant dominance genetic effects across all the traits studied in the three datasets used. Regarding the SNPs that were detected with dominance associations, further investigation is needed to determine their relevance in crossbreeding programs assuming that dominance effects are the main cause of (or contribute usefully to) heterosis.

Functional Annotation of All Salmonid Genomes (FAASG): an international initiative supporting future salmonid research, conservation and aquaculture.

We describe an emerging initiative - the 'Functional Annotation of All Salmonid Genomes' (FAASG), which will leverage the extensive trait diversity that has evolved since a whole genome duplication event in the salmonid ancestor, to develop an integrative understanding of the functional genomic basis of phenotypic variation. The outcomes of FAASG will have diverse applications, ranging from improved understanding of genome evolution, to improving the efficiency and sustainability of aquaculture production, supporting the future of fundamental and applied research in an iconic fish lineage of major societal importance.

Genome-wide analysis of the transcriptional response to porcine reproductive and respiratory syndrome virus infection at the maternal/fetal interface and in the fetus.

BACKGROUND: Porcine Reproductive and Respiratory Syndrome Virus (PRRSV) infection of pregnant pigs can result in congenital infection and ultimately fetal death. Little is known about immune responses to infection at the maternal-fetal interface and in the fetus itself, or the molecular events behind virus transmission and disease progression in the fetus. To investigate these processes, RNA-sequencing of two sites, uterine endothelium with adherent placental tissue and fetal thymus, was performed 21 days post-challenge on four groups of fetuses selected from a large PRRSV challenge experiment of pregnant gilts: control (CON), uninfected (UNINF), infected (INF), and meconium-stained (MEC) (n = 12/group). Transcriptional analyses consisted of multiple contrasts between groups using two approaches: differential gene expression analysis and weighted gene co-expression network analysis (WGCNA). Biological functions, pathways, and regulators enriched for differentially expressed genes or module members were identified through functional annotation analyses. Expression data were validated by reverse transcription quantitative polymerase chain reaction (RTqPCR) carried out for 16 genes of interest. RESULTS: The immune response to infection in endometrium was mainly adaptive in nature, with the most upregulated genes functioning in either humoral or cell-mediated immunity. In contrast, the expression profile of infected fetal thymus revealed a predominantly innate immune response to infection, featuring the upregulation of genes regulated by type I interferon and pro-inflammatory cytokines. Fetal infection was associated with an increase in viral load coupled with a reduction in T cell signaling in the endometrium that could be due to PRRSV-controlled apoptosis of uninfected bystander cells. There was also evidence for a reduction in TWIST1 activity, a transcription factor involved in placental implantation and maturation, which could facilitate virus transmission or fetal pathology through dysregulation of placental function. Finally, results suggested that events within the fetus could also drive fetal pathology. Thymus samples of meconium-stained fetuses exhibited an increase in the expression of pro-inflammatory cytokine and granulocyte genes previously implicated in swine infectious disease pathology. CONCLUSIONS: This study identified major differences in the response to PRRSV infection in the uterine endometrium and fetus at the gene expression level, and provides insight into the molecular basis of virus transmission and disease progression.

Bioinformatic analyses in early host response to Porcine Reproductive and Respiratory Syndrome virus (PRRSV) reveals pathway differences between pigs with alternate genotypes for a major host response QTL.

BACKGROUND: A region on Sus scrofa chromosome 4 (SSC4) surrounding single nucleotide polymorphism (SNP) marker WUR10000125 (WUR) has been reported to be strongly associated with both weight gain and serum viremia in pigs after infection with PRRS virus (PRRSV). A proposed causal mutation in the guanylate binding protein 5 gene (GBP5) is predicted to truncate the encoded protein. To investigate transcriptional differences between WUR genotypes in early host response to PRRSV infection, an RNA-seq experiment was performed on globin depleted whole blood RNA collected on 0, 4, 7, 10 and 14 days post-infection (dpi) from eight littermate pairs with one AB (favorable) and one AA (unfavorable) WUR genotype animal per litter. RESULTS: Gene Ontology (GO) enrichment analysis of transcripts that were differentially expressed (DE) between dpi across both genotypes revealed an inflammatory response for all dpi when compared to day 0. However, at the early time points of 4 and 7dpi, several GO terms had higher enrichment scores compared to later dpi, including inflammatory response (p < 10(-7)), specifically regulation of NFkappaB (p < 0.01), cytokine, and chemokine activity (p < 0.01). At 10 and 14dpi, GO term enrichment indicated a switch to DNA damage response, cell cycle checkpoints, and DNA replication. Few transcripts were DE between WUR genotypes on individual dpi or averaged over all dpi, and little enrichment of any GO term was found. However, there were differences in expression patterns over time between AA and AB animals, which was confirmed by genotype-specific expression patterns of several modules that were identified in weighted gene co-expression network analyses (WGCNA). Minor differences between AA and AB animals were observed in immune response and DNA damage response (p = 0.64 and p = 0.11, respectively), but a significant effect between genotypes pointed to a difference in ion transport/homeostasis and the participation of G-coupled protein receptors (p = 8e-4), which was reinforced by results from regulatory and phenotypic impact factor analyses between genotypes. CONCLUSION: We propose these pathway differences between WUR genotypes are the result of the inability of the truncated GBP5 of the AA genotyped pigs to inhibit viral entry and replication as quickly as the intact GBP5 protein of the AB genotyped pigs.

Genetic and genomic basis of antibody response to porcine reproductive and respiratory syndrome (PRRS) in gilts and sows.

BACKGROUND: Our recent research showed that antibody response to porcine reproductive and respiratory syndrome (PRRS), measured as sample-to-positive (S/P) ratio, is highly heritable and has a high genetic correlation with reproductive performance during a PRRS outbreak. Two major quantitative trait loci (QTL) on Sus scrofa chromosome 7 (SSC7; QTLMHC and QTL130) accounted for ~40 % of the genetic variance for S/P. Objectives of this study were to estimate genetic parameters for PRRS S/P in gilts during acclimation, identify regions associated with S/P, and evaluate the accuracy of genomic prediction of S/P across populations with different prevalences of PRRS and using different single nucleotide polymorphism (SNP) sets. METHODS: Phenotypes and high-density SNP genotypes of female pigs from two datasets were used. The outbreak dataset included 607 animals from one multiplier herd, whereas the gilt acclimation (GA) dataset included data on 2364 replacement gilts from seven breeding companies placed on health-challenged farms. Genomic prediction was evaluated using GA for training and validation, and using GA for training and outbreak for validation. Predictions were based on SNPs across the genome (SNPAll), SNPs in one (SNPMHC and SNP130) or both (SNPSSC7) QTL, or SNPs outside the QTL (SNPRest). RESULTS: Heritability of S/P in the GA dataset increased with the proportion of PRRS-positive animals in the herd (from 0.28 to 0.47). Genomic prediction accuracies ranged from low to moderate. Average accuracies were highest when using only the 269 SNPs in both QTL regions (SNPSSC7, with accuracies of 0.39 and 0.31 for outbreak and GA validation datasets, respectively. Average accuracies for SNPALL, SNPMHC, SNP130, and SNPRest were, respectively, 0.26, 0.39, 0.21, and 0.05 for the outbreak, and 0.28, 0.25, 0.22, and 0.12, for the GA validation datasets. CONCLUSIONS: Moderate genomic prediction accuracies can be obtained for PRRS antibody response using SNPs located within two major QTL on SSC7, while the rest of the genome showed limited predictive ability. Results were obtained using data from multiple genetic sources and farms, which further strengthens these findings. Further research is needed to validate the use of S/P ratio as an indicator trait for reproductive performance during PRRS outbreaks.

Comparison of host genetic factors influencing pig response to infection with two North American isolates of porcine reproductive and respiratory syndrome virus.

BACKGROUND: Porcine reproductive and respiratory syndrome (PRRS) is one of the most important swine diseases in the world and genetic selection of pigs for increased resistance to PRRS is an attractive method to improve the health status of the swine herd. This study compared phenotypic and genetic responses to infection with one of two genetically distinct type 2 PRRS virus (PRRSV) isolates: NVSL-97-7895 (NVSL) and KS-2006-72109 (KS06), and evaluated whether the single nucleotide polymorphism (SNP) WUR10000125 (WUR) on chromosome 4 that was associated with viral load and weight gain under infection with NVSL also has an effect on response to infection across North American PRRSV isolates. Wood's lactation curve was fitted to repeated viremia measurements to derive five curve characteristics that were evaluated. RESULTS: Infection with NVSL was characterized by reaching a 14 ± 2 % higher peak viremia (PV) 2.5 ± 0.6 days earlier (time to peak; TP) than KS06, followed by 36 ± 1 % faster virus clearance, which occurred 3.9 ± 0.7 days sooner. Weight gain from 0 to 42 days post-infection (WG) tended to be higher under infection with KS06 than NVSL (3.7 ± 1.5 kg). Estimates of heritability were moderate for both PRRSV isolates for viral load from 0 to 21 days post-infection (VL) (NVSL: 0.31 ± 0.06; KS06: 0.51 ± 0.09) and WG (NVSL: 0.33 ± 0.06; KS06: 0.31 ± 0.09). Strong negative genetic correlations were observed between VL and WG for both NVSL (-0.74 ± 0.10) and KS06 (-0.52 ± 0.17) infected pigs. Pigs with genotype AB at the WUR SNP had a more desirable phenotype than AA pigs for all traits under infection with NVSL, but only for VL and PV with KS06; effects on other traits were smaller and not significantly different from zero (P > 0.05). Genetic correlations of host response between isolates were strong for VL, WG and PV. Accounting for WUR genotype had little impact on these correlations, suggesting that response to PRRSV infection has a substantial polygenic component that is common between these two isolates. CONCLUSIONS: These results suggest that the KS06 PRRSV isolate is less virulent than NVSL but that genetic selection for increased resistance to either of these genetically distinct isolates is expected to increase resistance to the other isolate.

GO-FAANG meeting: a Gathering On Functional Annotation of Animal Genomes.

The Functional Annotation of Animal Genomes (FAANG) Consortium recently held a Gathering On FAANG (GO-FAANG) Workshop in Washington, DC on October 7-8, 2015. This consortium is a grass-roots organization formed to advance the annotation of newly assembled genomes of domesticated and non-model organisms (www.faang.org). The workshop gathered together from around the world a group of 100+ genome scientists, administrators, representatives of funding agencies and commodity groups to discuss the latest advancements of the consortium, new perspectives, next steps and implementation plans. The workshop was streamed live and recorded, and all talks, along with speaker slide presentations, are available at www.faang.org. In this report, we describe the major activities and outcomes of this meeting. We also provide updates on ongoing efforts to implement discussions and decisions taken at GO-FAANG to guide future FAANG activities. In summary, reference datasets are being established under pilot projects; plans for tissue sets, morphological classification and methods of sample collection for different tissues were organized; and core assays and data and meta-data analysis standards were established.

Identification of a putative quantitative trait nucleotide in guanylate binding protein 5 for host response to PRRS virus infection.

BACKGROUND: Previously, we identified a major quantitative trait locus (QTL) for host response to Porcine Respiratory and Reproductive Syndrome virus (PRRSV) infection in high linkage disequilibrium (LD) with SNP rs80800372 on Sus scrofa chromosome 4 (SSC4). RESULTS: Within this QTL, guanylate binding protein 5 (GBP5) was differentially expressed (DE) (p < 0.05) in blood from AA versus AB rs80800372 genotyped pigs at 7,11, and 14 days post PRRSV infection. All variants within the GBP5 transcript in LD with rs80800372 exhibited allele specific expression (ASE) in AB individuals (p < 0.0001). A transcript re-assembly revealed three alternatively spliced transcripts for GBP5. An intronic SNP in GBP5, rs340943904, introduces a splice acceptor site that inserts five nucleotides into the transcript. Individuals homozygous for the unfavorable AA genotype predominantly produced this transcript, with a shifted reading frame and early stop codon that truncates the 88 C-terminal amino acids of the protein. RNA-seq analysis confirmed this SNP was associated with differential splicing by QTL genotype (p < 0.0001) and this was validated by quantitative capillary electrophoresis (p < 0.0001). The wild-type transcript was expressed at a higher level in AB versus AA individuals, whereas the five-nucleotide insertion transcript was the dominant form in AA individuals. Splicing and ASE results are consistent with the observed dominant nature of the favorable QTL allele. The rs340943904 SNP was also 100 % concordant with rs80800372 in a validation population that possessed an alternate form of the favorable B QTL haplotype. CONCLUSIONS: GBP5 is known to play a role in inflammasome assembly during immune response. However, the role of GBP5 host genetic variation in viral immunity is novel. These findings demonstrate that rs340943904 is a strong candidate causal mutation for the SSC4 QTL that controls variation in host response to PRRSV.

Gene co-expression network analysis identifies porcine genes associated with variation in Salmonella shedding.

BACKGROUND: Salmonella enterica serovar Typhimurium is a gram-negative bacterium that can colonise the gut of humans and several species of food producing farm animals to cause enteric or septicaemic salmonellosis. While many studies have looked into the host genetic response to Salmonella infection, relatively few have used correlation of shedding traits with gene expression patterns to identify genes whose variable expression among different individuals may be associated with differences in Salmonella clearance and resistance. Here, we aimed to identify porcine genes and gene co-expression networks that differentiate distinct responses to Salmonella challenge with respect to faecal Salmonella shedding. RESULTS: Peripheral blood transcriptome profiles from 16 pigs belonging to extremes of the trait of faecal Salmonella shedding counts recorded up to 20 days post-inoculation (low shedders (LS), n = 8; persistent shedders (PS), n = 8) were generated using RNA-sequencing from samples collected just before (day 0) and two days after (day 2) Salmonella inoculation. Weighted gene co-expression network analysis (WGCNA) of day 0 samples identified four modules of co-expressed genes significantly correlated with Salmonella shedding counts upon future challenge. Two of those modules consisted largely of innate immunity related genes, many of which were significantly up-regulated at day 2 post-inoculation. The connectivity at both days and the mean gene-wise expression levels at day 0 of the genes within these modules were higher in networks constructed using LS samples alone than those using PS alone. Genes within these modules include those previously reported to be involved in Salmonella resistance such as SLC11A1 (formerly NRAMP1), TLR4, CD14 and CCR1 and those for which an association with Salmonella is novel, for example, SIGLEC5, IGSF6 and TNFSF13B. CONCLUSIONS: Our analysis integrates gene co-expression network analysis, gene-trait correlations and differential expression to provide new candidate regulators of Salmonella shedding in pigs. The comparatively higher expression (also confirmed in an independent dataset) and the significantly higher connectivity of genes within the Salmonella shedding associated modules in LS compared to PS even before Salmonella challenge may be factors that contribute to the decreased faecal Salmonella shedding observed in LS following challenge.

Increasing gene discovery and coverage using RNA-seq of globin RNA reduced porcine blood samples.

BACKGROUND: Transcriptome analysis of porcine whole blood has several applications, which include deciphering genetic mechanisms for host responses to viral infection and vaccination. The abundance of alpha- and beta-globin transcripts in blood, however, impedes the ability to cost-effectively detect transcripts of low abundance. Although protocols exist for reduction of globin transcripts from human and mouse/rat blood, preliminary work demonstrated these are not useful for porcine blood Globin Reduction (GR). Our objectives were to develop a porcine specific GR protocol and to evaluate the GR effects on gene discovery and sequence read coverage in RNA-sequencing (RNA-seq) experiments. RESULTS: A GR protocol for porcine blood samples was developed using RNase H with antisense oligonucleotides specifically targeting porcine hemoglobin alpha (HBA) and beta (HBB) mRNAs. Whole blood samples (n = 12) collected in Tempus tubes were used for evaluating the efficacy and effects of GR on RNA-seq. The HBA and HBB mRNA transcripts comprised an average of 46.1% of the mapped reads in pre-GR samples, but those reads reduced to an average of 8.9% in post-GR samples. Differential gene expression analysis showed that the expression level of 11,046 genes were increased, whereas 34 genes, excluding HBA and HBB, showed decreased expression after GR (FDR <0.05). An additional 815 genes were detected only in post-GR samples. CONCLUSIONS: Our porcine specific GR primers and protocol minimize the number of reads of globin transcripts in whole blood samples and provides increased coverage as well as accuracy and reproducibility of transcriptome analysis. Increased detection of low abundance mRNAs will ensure that studies relying on transcriptome analyses do not miss information that may be vital to the success of the study.

Variation in time and magnitude of immune response and viremia in experimental challenges with Porcine circovirus 2b.

BACKGROUND: Porcine circovirus 2 is the primary agent responsible for inducing a group of associated diseases known as Porcine Circovirus Associated Diseases (PCVAD), which can have detrimental effects on production efficiency as well as causing significant mortality. The objective of this study was to evaluate variation in viral replication, immune response and growth across pigs (n = 974) from different crossbred lines. The approach used in this study was experimental infection with a PCV2b strain of pigs at an average of 43 days of age. RESULTS: The sequence of the PCV2b isolate used in the challenge was similar with a cluster of PCV2b isolates known to induce PCVAD and increased mortality rates. The swine leukocyte antigen class II (SLAII) profile of the population was diverse, with nine DQB1 haplotypes being present. Individual viremia and antibody profiles during challenge demonstrate variation in magnitude and time of viral surge and immune response. The correlations between PCV2 specific antibodies and average daily gain (ADG) were relatively low and varied between - 0.14 to 0.08 for IgM and -0.02 and 0.11 for IgG. In contrast, PCV2 viremia was an important driver of ADG decline following infection; a moderate negative correlation was observed between viral load and overall ADG (r = - 0.35, P < 0.001). The pigs with the lowest 10% level of viral load maintained a steady increase in weekly ADG (P < 0.0001) compared to the pigs that had the 10% greatest viral load (P < 0.55). In addition, the highly viremic group expressed higher IgM and IgG starting with d 14 and d 21 respectively, and higher tumor necrosis factor - alpha (TNF-α) at d 21 (P < 0.005), compared to low viremic group. CONCLUSIONS: Molecular sources of the observed differences in viremia and immune response could provide a better understanding of the host factors that influence the development of PCVAD and lead to improved knowledge of swine immunity.