The Genetic Architecture of Bovine Telomere Length in Early Life and Association With Animal Fitness.

Health and survival are key goals for selective breeding in farm animals. Progress, however, is often limited by the low heritability of these animal fitness traits in addition to measurement difficulties. In this respect, relevant early-life biomarkers may be useful for breeding purposes. Telomere length (TL), measured in leukocytes, is a good candidate biomarker since TL has been associated with health, ageing, and stress in humans and other species. However, telomere studies are very limited in farm animals. Here, we examined the genetic background, genomic architecture, and factors affecting bovine TL measurements in early life, and the association of the latter with animal fitness traits expressed later in life associated with survival, longevity, health, and reproduction. We studied two TL measurements, one at birth (TLB) and another during the first lactation (TLFL) of a cow. We performed a genome-wide association study of dairy cattle TL, the first in a non-human species, and found that TLB and TLFL are complex, polygenic, moderately heritable, and highly correlated traits. However, genomic associations with distinct chromosomal regions were identified for the two traits suggesting that their genomic architecture is not identical. This is reflected in changes in TL throughout an individual's life. TLB had a significant association with survival, length of productive life and future health status of the animal, and could be potentially used as an early-life biomarker for disease predisposition and longevity in dairy cattle.

Characterization of copy number variants in a large multibreed population of beef and dairy cattle using high-density single nucleotide polymorphism genotype data.

Copy number variants (CNVs) are a form of genomic variation that changes the structure of the genome through deletion or duplication of stretches of DNA. The objective of the present study was to characterize CNVs in a large multibreed population of beef and dairy bulls. The CNVs were called on the autosomes of 5,551 cattle from 22 different beef and dairy breeds, using 2 freely available software suites, QuantiSNP and PennCNV. All CNVs were classified into either deletions or duplications. The median concordance between PennCNV and QuantiSNP, per animal, was 18.5% for deletions and 0% for duplications. The low concordance rate between PennCNV and QuantiSNP indicated that neither algorithm, by itself, could identify all CNVs in the population. In total, PennCNV and QuantiSNP collectively identified 747,129 deletions and 432,523 duplications; 80.2% of all duplications and 69.1% of all deletions were present only once in the population. Only 0.154% of all CNVs identified were present in more than 50 animals in the population. The distribution of the percentage of the autosomes that were composed of deletions, per animal, was positively skewed, as was the distribution for the percentage of the autosomes that were composed of duplications, per animal. The first quartile, median, and third quartile of the distribution of the percentage of the autosomes that were composed of deletions were 0.019%, 0.037%, and 0.201%, respectively. The first quartile, median, and third quartile of the distribution of the percentage of the autosomes that were composed of duplications were 0.013%, 0.028%, and 0.076%, respectively. The distributions of the number of deletions and duplications per animal were both positively skewed. The interquartile range for the number of deletions per animal in the population was between 16 and 117, whereas for duplications it was between 8 and 23. Per animal, there tended to be twice as many deletions as duplications. The distribution of the length of deletions was positively skewed, as was the distribution of the length of duplications. The interquartile range for the length of deletions in the population was between 25 and 101 kb, and for duplications the interquartile range was between 46 and 235 kb. Per animal, duplications tended to be twice as long as deletions. This study provides a description of the characteristics and distribution of CNVs in a large multibreed population of beef and dairy cattle.

Immune-associated traits measured in milk of Holstein-Friesian cows as proxies for blood serum measurements.

Previous work has highlighted that immune-associated (IA) traits measurable in blood are associated with health, productivity, and reproduction in dairy cows. The aim of the present study was to determine relationships between IA traits measured in blood serum and those simultaneously measured in milk as well as their association with disease phenotypes. All animals were Holstein-Friesian cows from the Langhill research herd (n = 546) housed at the SRUC Dairy Research Centre in Scotland. Milk and serum samples were collected on 20 separate occasions between July 2010 and March 2015 and analyzed by ELISA for haptoglobin (Hp), tumor necrosis factor-α (TNF-α), and natural antibodies binding keyhole limpet hemocyanin (NAbKLH) and lipopolysaccharide (NAbLPS). Data were analyzed using mixed linear models that included pedigree information. Analyses revealed positive phenotypic correlations between milk and serum NAb (0.59 ≤ r ≤ 0.77), Hp (r = 0.37), and TNF-α (r = 0.12). Milk and serum NAb were also found to have a strong genetic correlation (0.81 ≤ r ≤ 0.94) and were genetically correlated with cow lameness (0.66 and 0.79 for milk NAbKLH and serum NAbLPS, respectively). Clinical mastitis was found to be phenotypically correlated with both milk and serum Hp (0.09 ≤ r ≤ 0.23). Serum Hp was also strongly genetically correlated with other cellular IA traits such as percent NKp46+ (a natural killer cell marker; 0.35) and percent peripheral blood mononuclear cells (PBMC; -0.90). Similarly, genetic correlations were found to exist between serum TNF-α and percent NKp46+ (0.22), percent PBMC (0.41), and percent lymphocytes (0.47). Excluding serum Hp, all milk and serum IA traits were repeatable, ranging from 0.11 (milk Hp) to 0.43 (serum NAbLPS). Between-animal variation was highest in milk and serum NAb (0.34-0.43) and significant estimates of heritability were also observed in milk and serum NAb (0.17-0.37). Our findings show that certain IA traits, such as NAbKLH and NAbLPS, found in milk and serum are strongly correlated and highlight the potential of using routinely collected milk samples as a less invasive and cost-effective source of informative data for predictive modeling of animal IA traits.

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.

Herd-specific random regression carcass profiles for beef cattle after adjustment for animal genetic merit.

Abattoir data are an important source of information for the genetic evaluation of carcass traits, but also for on-farm management purposes. The present study aimed to quantify the contribution of herd environment to beef carcass characteristics (weight, conformation score and fat score) with particular emphasis on generating finishing herd-specific profiles for these traits across different ages at slaughter. Abattoir records from 46,115 heifers and 78,790 steers aged between 360 and 900days, and from 22,971 young bulls aged between 360 and 720days, were analysed. Finishing herd-year and animal genetic (co)variance components for each trait were estimated using random regression models. Across slaughter age and gender, the ratio of finishing herd-year to total phenotypic variance ranged from 0.31 to 0.72 for carcass weight, 0.21 to 0.57 for carcass conformation and 0.11 to 0.44 for carcass fat score. These parameters indicate that the finishing herd environment is an important contributor to carcass trait variability and amenable to improvement with management practices.

Estimating genetic and phenotypic parameters of cellular immune-associated traits in dairy cows.

Data collected from an experimental Holstein-Friesian research herd were used to determine genetic and phenotypic parameters of innate and adaptive cellular immune-associated traits. Relationships between immune-associated traits and production, health, and fertility traits were also investigated. Repeated blood leukocyte records were analyzed in 546 cows for 9 cellular immune-associated traits, including percent T cell subsets, B cells, NK cells, and granulocytes. Variance components were estimated by univariate analysis. Heritability estimates were obtained for all 9 traits, the highest of which were observed in the T cell subsets percent CD4+, percent CD8+, CD4+:CD8+ ratio, and percent NKp46+ cells (0.46, 0.41, 0.43 and 0.42, respectively), with between-individual variation accounting for 59 to 81% of total phenotypic variance. Associations between immune-associated traits and production, health, and fertility traits were investigated with bivariate analyses. Strong genetic correlations were observed between percent NKp46+ and stillbirth rate (0.61), and lameness episodes and percent CD8+ (-0.51). Regarding production traits, the strongest relationships were between CD4+:CD8+ ratio and weight phenotypes (-0.52 for live weight; -0.51 for empty body weight). Associations between feed conversion traits and immune-associated traits were also observed. Our results provide evidence that cellular immune-associated traits are heritable and repeatable, and the noticeable variation between animals would permit selection for altered trait values, particularly in the case of the T cell subsets. The associations we observed between immune-associated, health, fertility, and production traits suggest that genetic selection for cellular immune-associated traits could provide a useful tool in improving animal health, fitness, and fertility.

Phenotypic effects of subclinical paratuberculosis (Johne's disease) in dairy cattle.

The effect of subclinical paratuberculosis (or Johne's disease) risk status on performance, health, and fertility was studied in 58,096 UK Holstein-Friesian cows with 156,837 lactations across lactations 1 to 3. Low-, medium-, and high-risk group categories were allocated to cows determined by a minimum of 4 ELISA milk tests taken at any time during their lactating life. Lactation curves of daily milk, protein, and fat yields and protein and fat percentage, together with loge-transformed somatic cell count, were estimated using a random regression model to quantify differences between risk groups. The effect of subclinical paratuberculosis risk groups on fertility, lactation-average somatic cell count, and mastitis were analyzed using linear regression fitting risk group as a fixed effect. Milk yield losses associated with high-risk cows compared with low-risk cows in lactations 1, 2, and 3 for mean daily yield were 0.34, 1.05, and 1.61kg; likewise, accumulated 305-d yields were 103, 316, and 485kg, respectively. The total loss was 904kg over the first 3 lactations. Protein and fat yield losses associated with high-risk cows were significant, but primarily a feature of decreasing milk yield. Similar trends were observed for both test-day and lactation-average somatic cell count measures with higher somatic cell counts from medium- and high-risk cows compared with low-risk cows, and differences were in almost all cases significant. Likewise, mastitis incidence was significantly higher in high-risk cows compared with low-risk cows in lactations 2 and 3. Whereas the few significant differences between risk groups among fertility traits were inconsistent with no clear trend. These results are expected to be conservative, as some animals that were considered negative may become positive after the timeframe of this study, particularly if the animal was tested when relatively young. However, the magnitude of milk yield losses together with higher somatic cell counts and an increase in mastitis incidence should motivate farmers to implement the appropriate control measures to reduce the spread of the disease.

Evaluating alternate models to estimate genetic parameters of calving traits in United Kingdom Holstein-Friesian dairy cattle.

BACKGROUND: The focus in dairy cattle breeding is gradually shifting from production to functional traits and genetic parameters of calving traits are estimated more frequently. However, across countries, various statistical models are used to estimate these parameters. This study evaluates different models for calving ease and stillbirth in United Kingdom Holstein-Friesian cattle. METHODS: Data from first and later parity records were used. Genetic parameters for calving ease, stillbirth and gestation length were estimated using the restricted maximum likelihood method, considering different models i.e. sire (-maternal grandsire), animal, univariate and bivariate models. Gestation length was fitted as a correlated indicator trait and, for all three traits, genetic correlations between first and later parities were estimated. Potential bias in estimates was avoided by acknowledging a possible environmental direct-maternal covariance. The total heritable variance was estimated for each trait to discuss its theoretical importance and practical value. Prediction error variances and accuracies were calculated to compare the models. RESULTS AND DISCUSSION: On average, direct and maternal heritabilities for calving traits were low, except for direct gestation length. Calving ease in first parity had a significant and negative direct-maternal genetic correlation. Gestation length was maternally correlated to stillbirth in first parity and directly correlated to calving ease in later parities. Multi-trait models had a slightly greater predictive ability than univariate models, especially for the lowly heritable traits. The computation time needed for sire (-maternal grandsire) models was much smaller than for animal models with only small differences in accuracy. The sire (-maternal grandsire) model was robust when additional genetic components were estimated, while the equivalent animal model had difficulties reaching convergence. CONCLUSIONS: For the evaluation of calving traits, multi-trait models show a slight advantage over univariate models. Extended sire models (-maternal grandsire) are more practical and robust than animal models. Estimated genetic parameters for calving traits of UK Holstein cattle are consistent with literature. Calculating an aggregate estimated breeding value including direct and maternal values should encourage breeders to consider both direct and maternal effects in selection decisions.