Genome wide association study of thyroid hormone levels following challenge with porcine reproductive and respiratory syndrome virus

Introduction: Porcine reproductive and respiratory syndrome virus (PRRSV) causes respiratory disease in piglets and reproductive disease in sows. Piglet and fetal serum thyroid hormone (i.e., T3 and T4) levels decrease rapidly in response to Porcine reproductive and respiratory syndrome virus infection. However, the genetic control of T3 and T4 levels during infection is not completely understood. Our objective was to estimate genetic parameters and identify quantitative trait loci (QTL) for absolute T3 and/or T4 levels of piglets and fetuses challenged with Porcine reproductive and respiratory syndrome virus. Methods: Sera from 5-week-old pigs (N = 1792) at 11 days post inoculation (DPI) with Porcine reproductive and respiratory syndrome virus were assayed for T3 levels (piglet_T3). Sera from fetuses (N = 1,267) at 12 or 21 days post maternal inoculation (DPMI) with Porcine reproductive and respiratory syndrome virus of sows (N = 145) in late gestation were assayed for T3 (fetal_T3) and T4 (fetal_T4) levels. Animals were genotyped using 60 K Illumina or 650 K Affymetrix single nucleotide polymorphism (SNP) panels. Heritabilities, phenotypic correlations, and genetic correlations were estimated using ASREML; genome wide association studies were performed for each trait separately using Julia for Whole-genome Analysis Software (JWAS). Results: All three traits were low to moderately heritable (10%-16%). Phenotypic and genetic correlations of piglet_T3 levels with weight gain (0-42 DPI) were 0.26 ± 0.03 and 0.67 ± 0.14, respectively. Nine significant quantitative trait loci were identified for piglet_T3, on Sus scrofa chromosomes (SSC) 3, 4, 5, 6, 7, 14, 15, and 17, and collectively explaining 30% of the genetic variation (GV), with the largest quantitative trait loci identified on SSC5, explaining 15% of the genetic variation. Three significant quantitative trait loci were identified for fetal_T3 on SSC1 and SSC4, which collectively explained 10% of the genetic variation. Five significant quantitative trait loci were identified for fetal_T4 on SSC1, 6, 10, 13, and 15, which collectively explained 14% of the genetic variation. Several putative immune-related candidate genes were identified, including CD247, IRF8, and MAPK8. Discussion: Thyroid hormone levels following Porcine reproductive and respiratory syndrome virus infection were heritable and had positive genetic correlations with growth rate. Multiple quantitative trait loci with moderate effects were identified for T3 and T4 levels during challenge with Porcine reproductive and respiratory syndrome virus and candidate genes were identified, including several immune-related genes. These results advance our understanding of growth effects of both piglet and fetal response to Porcine reproductive and respiratory syndrome virus infection, revealing factors associated with genomic control of host resilience.

Host response to successive challenges with lentogenic and velogenic Newcastle disease virus in local chickens of Ghana

Newcastle disease (ND) is a highly contagious viral disease that constantly threatens poultry production. The velogenic (highly virulent) form of ND inflicts the most damage and can lead to 100% mortality in unvaccinated village chicken flocks. This study sought to characterize responses of local chickens in Ghana after challenging them with lentogenic and velogenic Newcastle disease virus (NDV) strains. At 4 wk of age, chicks were challenged with lentogenic NDV. Traits measured were pre- and post-lentogenic infection growth rates (GR), viral load at 2 and 6 d post-lentogenic infection (DPI), viral clearance rate and antibody levels at 10 DPI. Subsequently, the chickens were naturally exposed to velogenic NDV (vNDV) after anti-NDV antibody titers had waned to levels ≤1:1,700. Body weights and blood samples were again collected for analysis. Finally, chickens were euthanized and lesion scores (LS) across tissues were recorded. Post-velogenic exposure GR; antibody levels at 21 and 34 days post-velogenic exposure (DPE); LS for trachea, proventriculus, intestines, and cecal tonsils; and average LS across tissues were measured. Variance components and heritabilities were estimated for all traits using univariate animal models. Mean pre- and post-lentogenic NDV infection GRs were 6.26 g/day and 7.93 g/day, respectively, but mean post-velogenic NDV exposure GR was −1.96 g/day. Mean lesion scores ranged from 0.52 (trachea) to 1.33 (intestine), with males having significantly higher (P < 0.05) lesion scores compared to females. Heritability estimates for the lentogenic NDV trial traits ranged from moderate (0.23) to high (0.55) whereas those for the vNDV natural exposure trial were very low (≤ 0.08). Therefore, in contrast to the vNDV exposure trial, differences in the traits measured in the lentogenic challenge were more affected by genetics and thus selection for these traits may be more feasible compared to those following vNDV exposure. Our results can form the basis for identifying local chickens with improved resilience in the face of NDV infection for selective breeding to improve productivity.

Whole Genome Analyses Reveal Novel Genes Associated with Chicken Adaptation to Tropical and Frigid Environments

INTRODUCTION

Investigating the genetic footprints of historical temperature selection can get insights to the local adaptation and feasible influences of climate change on long-term population dynamics.

OBJECT

Chicken is a significative species to study genetic adaptation on account of its similar domestication track related to human activity with the most diversified varieties. Yet, few studies have demonstrated the genetic signatures of its adaptation to naturally tropical and frigid environments.

METHOD

Here, we generated whole genome resequencing of 119 domesticated chickens in China including the following breeds which are in order of breeding environmental temperature from more tropical to more frigid: Wenchang chicken (WCC), green-shell chicken (GSC), Tibetan chicken (TBC), and Lindian chicken (LDC).

RESULTS

Our results showed WCC branched off earlier than LDC with an evident genetic admixture between WCC and LDC, suggesting their closer genetic relationship. Further comparative genomic analyses solute carrier family 33 member 1 (SLC33A1) and thyroid stimulating hormone receptor (TSHR) genes exhibited stronger signatures for positive selection in the genome of the more tropical WCC. Furthermore, genotype data from about 3,000 African local ecotypes confirmed that allele frequencies of single nucleotide polymorphisms (SNPs) in these 2 genes appeared strongly associated with tropical environment adaptation. In addition, the NADH:ubiquinone oxidoreductase subunit S4 (NDUFS4) gene exhibited a strong signature for positive selection in the LDC genome, and SNPs with marked allele frequency differences indicated a significant relationship with frigid environment adaptation.

CONCLUSION

Our findings partially clarify how selection footprints from environmental temperature stress can lead to advantageous genomic adaptions to tropical and frigid environments in poultry and provide a valuable resource for selective breeding of chickens.

44 Development and Application of Novel Phenotyping Technologies to Swine Production and Breeding

Advancement within numerous categories of technology are beginning to see application for improving production and genetic advancement in pigs. Digital technologies, particularly video, are creating a step change in improved animal environments, biosecurity, vaccine management and disease intervention. For breeders, video technology creates opportunity to begin assessing animal behavior on a useful scale and to expand the quantity and type of animals measured for growth rate, body composition, structural integrity and feed efficiency. A seemingly endless stream of ‘omics’ data is becoming less costly to obtain and includes such datasets as common SNP arrays in addition to DNA sequence, transcriptome, proteome, metabolome, gut microbiome and more. Leveraging the capabilities of these data to predict a phenotype is key to unlocking the full potential of ‘omics’ data. Novel, at least to animal scientists, analytical techniques are required to understand gene networks and their associated expression profiles. These techniques apply equally to video data, which in its purest form, is an interpretation of complex patterns. Scaling and applying these techniques will be important over the coming decade. There exists today a very applied approach to using these technologies in the swine industry. One of the most complex and well-trained pattern recognition technologies lies within our own brain. Experienced managers reviewing video of pigs within facilities, team members entering and performing critical activities impacting biosecurity and coordinating these observations with simple interventions can make large, positive improvements in productivity. Investing in hardware and time will produce results that reduce stressors, improve processes that protect the farm from pathogens and result in better overall animal well-being. This talk will focus on examples of current applications that are of benefit today. In addition, we hope to promote discussion around priorities for genetic improvement, pig management and the practical needs to move a variety of technologies from research to development in the field.

Complex immunogenomic control of serum thyroid hormone levels in fetuses and piglets challenged with Porcine Reproductive and Respiratory Syndrome Virus

Porcine Reproductive and Respiratory Syndrome Virus (PRRSV) causes respiratory disease in piglets and reproductive disease in sows. However, the interaction between host immunity and hormone homeostasis is not completely understood but piglet and fetal serum thyroid hormone (i.e., T3 and T4) levels decrease rapidly in response to PRRSV infection. Our objective was to estimate genetic parameters, identify QTL, and uncover nearby candidate genes for T3 and/or T4 levels of piglets and fetuses challenged with PRRSV2. Serum from 5-week old pigs (N=1792) at 11 days post infection (dpi) with PRRSV2 were assayed for T3 levels (piglet_T3). While serum from fetuses (N=1267) at 12 or 21 dpi from sows (N=150) challenged with PRRSV2 in late gestation were assayed for T3 (fetal_T3) and T4 (fetal_T4) levels. Animals were genotyped for 60K or 650K SNPs, heritabilities were estimated and genome wide association studies were performed for each trait separately. We found all three traits to be moderately heritable (18–28%) and the significant QTL (N=37) across several chromosomes collectively explained 24–34% of the genetic variation. A shared QTL identified for all traits on chromosome 6 contained the IL12RB2 and IFI44 genes. The largest QTL identified for both piglet_T3 and fetal_T3 on chromosome 5 contained IL17RA. A large QTL identified for piglet_T3 on chromosome 7 contained 9 MHC genes (e.g., SLA-DMA, SLA-DOA, SLA-DQA, SLA-DRA, and MHC-IA7). In conclusion, our data support the concept that there are complex interactions between thyroid hormone levels and multiple immunogenetic pathways during PRRSV infection, and that there is potential to select for pleiotropic QTL to simultaneously improve host immunity for reproductive and respiratory PRRS.

A genome-wide analysis of cardiac lesions of pigs that die during transport: Is heart failure of in-transit-loss pigs associated with a heritable cardiomyopathy?

While heart failure is a primary cause of death for many in-transit-loss (ITL) pigs, the underlying cause of these deaths is not known. Cardiomyopathies are considered a common cause of heart failure in humans and often have a genetic component. The objective of this study was to determine if genes associated with cardiomyopathies could be identified in ITL pigs. Samples from the hearts of pigs that died during transport to an abattoir in Ontario, Canada were collected and genotyped along with samples from pigs that did not die during transport (ILT hearts: n = 149; non-ITL/control hearts: n = 387). Genome-wide analyses were carried out on each of the determined phenotypes (gross cardiac lesions) using a medium density single nucleotide polymorphism (SNP) chip and 500 kb windows/regions for analysis, with 250 kb regions of overlap. The distribution derived by a multidimensional scaling (MDS) analysis of all phenotypes demonstrated a lack of complete separation between phenotypes of affected and unaffected animals, which made diagnosis difficult. Although genetic differences were small, a few genes associated with dilated cardiomyopathy (DCM) and arrhythmogenic right ventricular cardiomyopathy (ARVM) were identified. In addition, multiple genes associated with cardiac arrhythmias and ventricular hypertrophy were identified that can possibly result in heart failure. The results of this preliminary study did not provide convincing evidence that a single, heritable cardiomyopathy is the cause of heart failure in ITL pigs.

An initial genome-wide investigation of protein-losing enteropathy in Gordon setters: Exploratory observations

The objective of this preliminary study was to identify genomic regions that may predispose Gordon setters from the United Kingdom to familial protein-losing enteropathy (PLE) at a young age. A total of 106 related Gordon setters was used, including 6 affected dogs from an affected litter, 6 case controls from the same litter, 10 related/affected dogs, and 84 related/unaffected dogs. Genomic DNA was collected from each Gordon setter and extracted from buccal mucosal swabs. Genotyping of affected and unaffected dogs was carried out using the Canine Illumina HD SNP array and data generated were analyzed with PLINK software, using fixation index (Fst) and runs of homozygosity (ROH) methods. Pairwise Fst analyses between the affected and unaffected Gordon setter dogs identified various regions of differentiation on chromosomes 10, 18, 21, and 23 that contained several important genes. These regions revealed 5 candidate genes, including RARB, TTC7A, SOCS5, PIGF, and RHOD, that are associated with human inflammatory bowel disease (IBD) and could potentially be associated with PLE in Gordon setters. Run of homozygosity (ROH) analyses revealed additional unique regions on chromosomes 15 and 17. These regions contained genes SYT1, UCN, and FNDC that could also be potential candidates for PLE in Gordon setters. The biological functions of the identified genes provided initial insights into the pathophysiology of PLE. Further large-scale studies are warranted to investigate the possible causality of these genomic regions and any possible genetic markers that could be used in predicting susceptibility to PLE syndrome.

Phenotypic variability and population structure analysis of Tanzanian free-range local chickens

BACKGROUND

Free-range local chickens (FRLC) farming is an important activity in Tanzania, however, they have not been well-characterized. This study aimed to phenotypically characterize three Tanzanian FRLCs and to determine their population structure. A total of 389 mature breeder chickens (324 females and 65 males) from three popular Tanzanian FRLC ecotypes (Kuchi, Morogoro-medium and Ching'wekwe) were used for the phenotypic characterization. Progenies of these chickens were utilized to assess population structure. The ecotypes were collected from four geographical zones across Tanzania: Lake, Central, Northern and Coastal zones. Body weights and linear measurements were obtained from the mature breeders, including body, neck, shanks, wingspan, chest girth, and shank girth. Descriptive statistics were utilized to characterize the chickens. Correlations between the linear measurements and differences among the means of measured linear traits between ecotypes and between sexes were assessed. A total of 1399 progeny chicks were genotyped using a chicken 600 K high density single nucleotide polymorphism (SNP) panel for determination of population structure.

RESULTS

The means for most traits were significantly higher in Kuchi relative to Ching'wekwe and Morogoro-medium. However, shank length and shank girth were similar between Kuchi and Morogoro-medium females. All traits were correlated with the exception of shank girth in Morogoro-medium. Admixture analyses revealed that Morogoro-medium and Ching'wekwe clustered together as one population, separate from Kuchi.

CONCLUSIONS

Phenotypic traits could be used to characterize FRLCs, however, there were variations in traits among individuals within ecotypes; therefore, complementary genomic methods should be considered to improve the characterization for selective breeding.

Genetic Basis of Response of Ghanaian Local Chickens to Infection With a Lentogenic Newcastle Disease Virus

Newcastle disease (ND) is a global threat to domestic poultry, especially in rural areas of Africa and Asia, where the loss of entire backyard local chicken flocks often threatens household food security and income. To investigate the genetics of Ghanaian local chicken ecotypes to Newcastle disease virus (NDV), in this study, three popular Ghanaian chicken ecotypes (regional populations) were challenged with a lentogenic NDV strain at 28 days of age. This study was conducted in parallel with a similar study that used three popular Tanzanian local chicken ecotypes and after two companion studies in the United States, using Hy-line Brown commercial laying birds. In addition to growth rate, NDV response traits were measured following infection, including anti-NDV antibody levels [pre-infection and 10 days post-infection (dpi)], and viral load (2 and 6 dpi). Genetic parameters were estimated, and two genome-wide association study analysis methods were used on data from 1,440 Ghanaian chickens that were genotyped on a chicken 600K Single Nucleotide Polymorphism (SNP) chip. Both Ghana and Tanzania NDV challenge studies revealed moderate to high (0.18 – 0.55) estimates of heritability for all traits, except viral clearance where the heritability estimate was not different from zero for the Tanzanian ecotypes. For the Ghana study, 12 quantitative trait loci (QTL) for growth and/or response to NDV from single-SNP analyses and 20 genomic regions that explained more than 1% of genetic variance using the Bayes B method were identified. Seven of these windows were also identified as having at least one significant SNP in the single SNP analyses for growth rate, anti-NDV antibody levels, and viral load at 2 and 6 dpi. An important gene for growth during stress, CHORDC1 associated with post-infection growth rate was identified as a positional candidate gene, as well as other immune related genes, including VAV2, IL12B, DUSP1, and IL17B. The QTL identified in the Ghana study did not overlap with those identified in the Tanzania study. However, both studies revealed QTL with genes vital for growth and immune response during NDV challenge. The Tanzania parallel study revealed an overlapping QTL on chromosome 24 for viral load at 6 dpi with the US NDV study in which birds were challenged with NDV under heat stress. This QTL region includes genes related to immune response, including TIRAP, ETS1, and KIRREL3. The moderate to high estimates of heritability and the identified QTL suggest that host response to NDV of local African chicken ecotypes can be improved through selective breeding to enhance increased NDV resistance and vaccine efficacy.

Genetic Analyses of Tanzanian Local Chicken Ecotypes Challenged with Newcastle Disease Virus

Newcastle Disease (ND) is a continuing global threat to domestic poultry, especially in developing countries, where severe outbreaks of velogenic ND virus (NDV) often cause major economic losses to households. Local chickens are of great importance to rural family livelihoods through provision of high-quality protein. To investigate the genetic basis of host response to NDV, three popular Tanzanian chicken ecotypes (regional populations) were challenged with a lentogenic (vaccine) strain of NDV at 28 days of age. Various host response phenotypes, including anti-NDV antibody levels (pre-infection and 10 days post-infection, dpi), and viral load (2 and 6 dpi) were measured, in addition to growth rate. We estimated genetic parameters and conducted genome-wide association study analyses by genotyping 1399 chickens using the Affymetrix 600K chicken SNP chip. Estimates of heritability of the evaluated traits were moderate (0.18–0.35). Five quantitative trait loci (QTL) associated with growth and/or response to NDV were identified by single-SNP analyses, with some regions explaining ≥1% of genetic variance based on the Bayes-B method. Immune related genes, such as ETS1, TIRAP, and KIRREL3, were located in regions associated with viral load at 6 dpi. The moderate estimates of heritability and identified QTL indicate that NDV response traits may be improved through selective breeding of chickens to enhance increased NDV resistance and vaccine efficacy in Tanzanian local ecotypes.

Detection of Selection Signatures Among Brazilian, Sri Lankan, and Egyptian Chicken Populations Under Different Environmental Conditions

Extreme environmental conditions are a major challenge in livestock production. Changes in climate, particularly those that contribute to weather extremes like drought or excessive humidity, may result in reduced performance and reproduction and could compromise the animal's immune function. Animal survival within extreme environmental conditions could be in response to natural selection and to artificial selection for production traits that over time together may leave selection signatures in the genome. The aim of this study was to identify selection signatures that may be involved in the adaptation of indigenous chickens from two different climatic regions (Sri Lanka = Tropical; Egypt = Arid) and in non-indigenous chickens that derived from human migration events to the generally tropical State of São Paulo, Brazil. To do so, analyses were conducted using fixation index (Fst) and hapFLK analyses. Chickens from Brazil (n = 156), Sri Lanka (n = 92), and Egypt (n = 96) were genotyped using the Affymetrix Axiom^®600k Chicken Genotyping Array. Pairwise Fst analyses among countries did not detect major regions of divergence between chickens from Sri Lanka and Brazil, with ecotypes/breeds from Brazil appearing to be genetically related to Asian-Indian (Sri Lanka) ecotypes. However, several differences were detected in comparisons of Egyptian with either Sri Lankan or Brazilian populations, and common regions of difference on chromosomes 2, 3 and 8 were detected. The hapFLK analyses for the three separate countries suggested unique regions that are potentially under selection on chromosome 1 for all three countries, on chromosome 4 for Sri Lankan, and on chromosomes 3, 5, and 11 for the Egyptian populations. Some of identified regions under selection with hapFLK analyses contained genes such as TLR3, SOCS2, EOMES, and NFAT5 whose biological functions could provide insights in understanding adaptation mechanisms in response to arid and tropical environments.

Evaluation of growth, deposition of back fat, and loin muscle for purebred Berkshire pigs housed in bedded hoop buildings

This study was conducted to evaluate the accretion of BW, back fat, and loin muscle from purebred Berkshire pigs raised in bedded hoop barns in Iowa. The growth of a total of 144 purebred Berkshire pigs (18 barrows and 18 gilts per trial) was evaluated from 4 trials (2 winter and 2 summer trials). Pigs were fed ad libitum utilizing a 5-phase standard corn-soybean meal feeding program that met or exceeded NRC nutrient requirements. Pigs were housed in bedded hoop barns (unheated) to approximate common niche market requirements. At 21-d intervals, pigs were serially weighed, and ultrasonic back fat depth and loin muscle area (LMA) measurements were taken. Live BW measurements began at the trial initiation at approximately 18 to 32 kg, but ultrasonic scans for 10th-rib back fat depth and LMA began at between 36 and 45 kg until market weight of about 122 ± 2.5 kg. The rate (µ) of live body growth (weight) and ultrasonic back fat depth were influenced ( < 0.01) by trial and sex, with no significant interactions between trial and sex. Both live BW and back fat deposition were significantly greater in trial 1 than all other trials (2, 3, and 4). The rate of accretion and maximum growth of LMA depth were not affected ( > 0.05) by trial and sex. Overall, barrows averaged 31 mm of back fat at 125 kg, whereas gilts had an average of about 23 mm at 121 kg of market weight. Results suggest that because of the sex differences in growth and back fat deposition between Berkshire barrows and gilts, it may be important to formulate their diets differently in commercial pork production systems.