A genome-wide association analysis for body weight at 35 days measured on 137,343 broiler chickens.

BACKGROUND: Body weight (BW) is an economically important trait in the broiler (meat-type chickens) industry. Under the assumption of polygenicity, a "large" number of genes with "small" effects is expected to control BW. To detect such effects, a large sample size is required in genome-wide association studies (GWAS). Our objective was to conduct a GWAS for BW measured at 35 days of age with a large sample size. METHODS: The GWAS included 137,343 broilers spanning 15 pedigree generations and 392,295 imputed single nucleotide polymorphisms (SNPs). A false discovery rate of 1% was adopted to account for multiple testing when declaring significant SNPs. A Bayesian ridge regression model was implemented, using AlphaBayes, to estimate the contribution to the total genetic variance of each region harbouring significant SNPs (1 Mb up/downstream) and the combined regions harbouring non-significant SNPs. RESULTS: GWAS revealed 25 genomic regions harbouring 96 significant SNPs on 13 Gallus gallus autosomes (GGA1 to 4, 8, 10 to 15, 19 and 27), with the strongest associations on GGA4 at 65.67-66.31 Mb (Galgal4 assembly). The association of these regions points to several strong candidate genes including: (i) growth factors (GGA1, 4, 8, 13 and 14); (ii) leptin receptor overlapping transcript (LEPROT)/leptin receptor (LEPR) locus (GGA8), and the STAT3/STAT5B locus (GGA27), in connection with the JAK/STAT signalling pathway; (iii) T-box gene (TBX3/TBX5) on GGA15 and CHST11 (GGA1), which are both related to heart/skeleton development); and (iv) PLAG1 (GGA2). Combined together, these 25 genomic regions explained ~ 30% of the total genetic variance. The region harbouring significant SNPs that explained the largest portion of the total genetic variance (4.37%) was on GGA4 (~ 65.67-66.31 Mb). CONCLUSIONS: To the best of our knowledge, this is the largest GWAS that has been conducted for BW in chicken to date. In spite of the identified regions, which showed a strong association with BW, the high proportion of genetic variance attributed to regions harbouring non-significant SNPs supports the hypothesis that the genetic architecture of BW35 is polygenic and complex. Our results also suggest that a large sample size will be required for future GWAS of BW35.

Telomere attrition rates are associated with weather conditions and predict productive lifespan in dairy cattle.

Telomere length is predictive of adult health and survival across vertebrate species. However, we currently do not know whether such associations result from among-individual differences in telomere length determined genetically or by early-life environmental conditions, or from differences in the rate of telomere attrition over the course of life that might be affected by environmental conditions. Here, we measured relative leukocyte telomere length (RLTL) multiple times across the entire lifespan of dairy cattle in a research population that is closely monitored for health and milk production and where individuals are predominantly culled in response to health issues. Animals varied in their change in RLTL between subsequent measurements and RLTL shortened more during early life and following hotter summers which are known to cause heat stress in dairy cows. The average amount of telomere attrition calculated over multiple repeat samples of individuals predicted a shorter productive lifespan, suggesting a link between telomere loss and health. TL attrition was a better predictor of when an animal was culled than their average TL or the previously for this population reported significant TL at the age of 1 year. Our present results support the hypothesis that TL is a flexible trait that is affected by environmental factors and that telomere attrition is linked to animal health and survival traits. Change in telomere length may represent a useful biomarker in animal welfare studies.

The usage of Mate Select, a web-based selection tool for pedigree dogs for promoting sustainable breeding.

BACKGROUND: Inbreeding is a phenomenon that accumulates through the mating of relatives within closed populations, such as pedigree dog breeds, and results in reduced genetic variation within breeds, and may lead to poorer health and fertility from inbreeding depression. The impact of inbreeding is driven by the selection and mating of parents, but information on choices to reduce inbreeding is difficult to assess for individual breeders. Tools to inform dog breeders on the current state of the inbreeding and the relationships among possible parents are potentially useful for providing guidance towards choices that are more beneficial to the breed. However, their utility depends on their usage and this study examines the usage of Mate Select, a web-based tool offered by The Kennel Club, covering 222 breeds for a period of 7 years following its launch in 2011. RESULTS: The average usage was 2830 searches/week in 2012 with a slight fall of 2.2% per year (P < 0.001) to 2480 searches/week in 2018. Of these, 4% originated from outside the UK, across all continents except Antarctica, with the majority coming from English speaking countries. Searches/week showed a cyclical pattern with two cycles of 26.0 and 50.1 weeks. Since Mate Select's launch there has been a steady increase in searches from mobile devices, from 11% in 2012 to 43% in 2018. For the 197 breeds with at least 10 dams registered with the Kennel Club during the study period, there was a relationship between usage and registrations, with the average number of searches as a multiple of the number of dams increasing from 2 to 10 for breeds with up to 70 dams and declining towards 2 again for the largest breeds with approximately 20,000 registered dams. However, there remained substantial variation among breeds of similar size, and breeds for which EBVs had become available during the study period had a 2.46 fold greater frequency of searches per registered bitch (P < 0.001), but this was not linked directly to the publication of EBVs. CONCLUSIONS: Mate Select has sustained and substantial usage, although there is also substantial variation in usage among breeds, which offers an opportunity to develop further guidance.

Genetic analysis of novel phenotypes for farm animal resilience to weather variability.

BACKGROUND: Climate change is expected to have a negative impact on food availability. While most efforts have been directed to reducing greenhouse gas emissions, complementary strategies are necessary to control the detrimental effects of climate change on farm animal performance. The objective of this study was to develop novel animal resilience phenotypes using reaction norm slopes, and examine their genetic and genomic parameters. A closely monitored dairy goat population was used for this purpose. RESULTS: Individual animals differed in their response to changing atmospheric temperature and a temperature-humidity index. Significant genetic variance and heritability estimates were derived for these animal resilience phenotypes. Furthermore, some resilience traits had a significant unfavourable genetic correlation with animal performance. Genome-wide association analyses identified several candidate genes related to animal resilience to environment change. CONCLUSIONS: Heritable variation exists among dairy goats in their production response to fluctuating weather variables. Results may inform future breeding programmes aimed to ensure efficient animal performance under changing climatic conditions.

Bovine telomere dynamics and the association between telomere length and productive lifespan.

Average telomere length (TL) in blood cells has been shown to decline with age in a range of vertebrate species, and there is evidence that TL is a heritable trait associated with late-life health and mortality in humans. In non-human mammals, few studies to date have examined lifelong telomere dynamics and no study has estimated the heritability of TL, despite these being important steps towards assessing the potential of TL as a biomarker of productive lifespan and health in livestock species. Here we measured relative leukocyte TL (RLTL) in 1,328 samples from 308 Holstein Friesian dairy cows and in 284 samples from 38 female calves. We found that RLTL declines after birth but remains relatively stable in adult life. We also calculated the first heritability estimates of RLTL in a livestock species which were 0.38 (SE = 0.03) and 0.32 (SE = 0.08) for the cow and the calf dataset, respectively. RLTL measured at the ages of one and five years were positively correlated with productive lifespan (p < 0.05). We conclude that bovine RLTL is a heritable trait, and its association with productive lifespan may be used in breeding programmes aiming to enhance cow longevity.

Longitudinal changes in telomere length and associated genetic parameters in dairy cattle analysed using random regression models.

Telomeres cap the ends of linear chromosomes and shorten with age in many organisms. In humans short telomeres have been linked to morbidity and mortality. With the accumulation of longitudinal datasets the focus shifts from investigating telomere length (TL) to exploring TL change within individuals over time. Some studies indicate that the speed of telomere attrition is predictive of future disease. The objectives of the present study were to 1) characterize the change in bovine relative leukocyte TL (RLTL) across the lifetime in Holstein Friesian dairy cattle, 2) estimate genetic parameters of RLTL over time and 3) investigate the association of differences in individual RLTL profiles with productive lifespan. RLTL measurements were analysed using Legendre polynomials in a random regression model to describe TL profiles and genetic variance over age. The analyses were based on 1,328 repeated RLTL measurements of 308 female Holstein Friesian dairy cattle. A quadratic Legendre polynomial was fitted to the fixed effect of age in months and to the random effect of the animal identity. Changes in RLTL, heritability and within-trait genetic correlation along the age trajectory were calculated and illustrated. At a population level, the relationship between RLTL and age was described by a positive quadratic function. Individuals varied significantly regarding the direction and amount of RLTL change over life. The heritability of RLTL ranged from 0.36 to 0.47 (SE = 0.05-0.08) and remained statistically unchanged over time. The genetic correlation of RLTL at birth with measurements later in life decreased with the time interval between samplings from near unity to 0.69, indicating that TL later in life might be regulated by different genes than TL early in life. Even though animals differed in their RLTL profiles significantly, those differences were not correlated with productive lifespan (p = 0.954).

Use and optimization of different sources of information for genomic prediction.

BACKGROUND: Molecular data is now commonly used to predict breeding values (BV). Various methods to calculate genomic relationship matrices (GRM) have been developed, with some studies proposing regression of coefficients back to the reference matrix of pedigree-based relationship coefficients (A). The objective was to compare the utility of two GRM: a matrix based on linkage analysis (LA) and anchored to the pedigree, i.e. [Formula: see text] and a matrix based on linkage disequilibrium (LD), i.e. [Formula: see text], using genomic and phenotypic data collected on 5416 broiler chickens. Furthermore, the effects of regressing the coefficients of [Formula: see text] back to A (LDA) and to [Formula: see text] (LDLA) were evaluated, using a range of weighting factors. The performance of the matrices and their composite products was assessed by the fit of the models to the data, and the empirical accuracy and bias of the BV that they predicted. The sensitivity to marker choice was examined by using two chips of equal density but including different single nucleotide polymorphisms (SNPs). RESULTS: The likelihood of models using GRM and composite matrices exceeded the likelihood of models based on pedigree alone and was highest with intermediate weighting factors for both the LDA and LDLA approaches. For these data, empirical accuracies were not strongly affected by the weighting factors, although they were highest when different sources of information were combined. The optimum weighting factors depended on the type of matrices used, as well as on the choice of SNPs from which the GRM were constructed. Prediction bias was strongly affected by the chip used and less by the form of the GRM. CONCLUSIONS: Our findings provide an empirical comparison of the efficacy of pedigree and genomic predictions in broiler chickens and examine the effects of fitting GRM with coefficients regressed back to a reference anchored to the pedigree, either A or [Formula: see text]. For the analysed dataset, the best results were obtained when [Formula: see text] was combined with relationships in A or [Formula: see text], with optimum weighting factors that depended on the choice of SNPs used. The optimum weighting factor for broiler body weight differed from weighting factors that were based on the density of SNPs and theoretically derived using generalised assumptions.

Method Specific Calibration Corrects for DNA Extraction Method Effects on Relative Telomere Length Measurements by Quantitative PCR.

Telomere length (TL) is increasingly being used as a biomarker in epidemiological, biomedical and ecological studies. A wide range of DNA extraction techniques have been used in telomere experiments and recent quantitative PCR (qPCR) based studies suggest that the choice of DNA extraction method may influence average relative TL (RTL) measurements. Such extraction method effects may limit the use of historically collected DNA samples extracted with different methods. However, if extraction method effects are systematic an extraction method specific (MS) calibrator might be able to correct for them, because systematic effects would influence the calibrator sample in the same way as all other samples. In the present study we tested whether leukocyte RTL in blood samples from Holstein Friesian cattle and Soay sheep measured by qPCR was influenced by DNA extraction method and whether MS calibration could account for any observed differences. We compared two silica membrane-based DNA extraction kits and a salting out method. All extraction methods were optimized to yield enough high quality DNA for TL measurement. In both species we found that silica membrane-based DNA extraction methods produced shorter RTL measurements than the non-membrane-based method when calibrated against an identical calibrator. However, these differences were not statistically detectable when a MS calibrator was used to calculate RTL. This approach produced RTL measurements that were highly correlated across extraction methods (r > 0.76) and had coefficients of variation lower than 10% across plates of identical samples extracted by different methods. Our results are consistent with previous findings that popular membrane-based DNA extraction methods may lead to shorter RTL measurements than non-membrane-based methods. However, we also demonstrate that these differences can be accounted for by using an extraction method-specific calibrator, offering researchers a simple means of accounting for differences in RTL measurements from samples extracted by different DNA extraction methods within a study.