Across-country genetic and genomic analyses of foot score traits in American and Australian Angus cattle.

BACKGROUND: Hoof structure and health are essential for the welfare and productivity of beef cattle. Therefore, we assessed the genetic and genomic background of foot score traits in American (US) and Australian (AU) Angus cattle and investigated the feasibility of performing genomic evaluations combining data for foot score traits recorded in US and AU Angus cattle. The traits evaluated were foot angle (FA) and claw set (CS). In total, 109,294 and ~ 1.12 million animals had phenotypic and genomic information, respectively. Four sets of analyses were performed: (1) genomic connectedness between US and AU Angus cattle populations and population structure, (2) estimation of genetic parameters, (3) single-step genomic prediction of breeding values, and (4) single-step genome-wide association studies for FA and CS. RESULTS: There was no clear genetic differentiation between US and AU Angus populations. Similar heritability estimates (FA: 0.22-0.24 and CS: 0.22-0.27) and moderate-to-high genetic correlations between US and AU foot scores (FA: 0.61 and CS: 0.76) were obtained. A joint-genomic prediction using data from both populations outperformed within-country genomic evaluations. A genomic prediction model considering US and AU datasets as a single population performed similarly to the scenario accounting for genotype-by-environment interactions (i.e., multiple-trait model considering US and AU records as different traits), even though the genetic correlations between countries were lower than 0.80. Common significant genomic regions were observed between US and AU for FA and CS. Significant single nucleotide polymorphisms were identified on the Bos taurus (BTA) chromosomes BTA1, BTA5, BTA11, BTA13, BTA19, BTA20, and BTA23. The candidate genes identified were primarily from growth factor gene families, including FGF12 and GDF5, which were previously associated with bone structure and repair. CONCLUSIONS: This study presents comprehensive population structure and genetic and genomic analyses of foot scores in US and AU Angus cattle populations, which are essential for optimizing the implementation of genomic selection for improved foot scores in Angus cattle breeding programs. We have also identified candidate genes associated with foot scores in the largest Angus cattle populations in the world and made recommendations for genomic evaluations for improved foot score traits in the US and AU.

Family climate influences next-generation family business leader effectiveness and work engagement.

Effective next-generation leadership is central to the multi-generational survival of family businesses. This study of 100 next-generation family business leaders found that business-owning families that openly express their opinions, take time to listen to each other, and squarely address difficult issues positively influence the development of the emotional and social intelligence competencies in next-generation family leaders that drive their leadership effectiveness. That kind of open and transparent communication in the family also makes it more likely next-generation leaders will be held accountable for their leadership performance by others, which increases the degree to which they are positively engaged with their work in the family firm. On the other hand, the results suggest that senior-generation family leaders who lead autocratically, a leadership style often observed in entrepreneurs who found family firms, make it less likely that next-generation family leaders will learn the emotional and social intelligence competencies that predict their leadership effectiveness. The study also found that autocratic senior-generation leaders negatively affect next-generation leader self-efficacy and make it less likely that others will hold them accountable, which limits their engagement with work in the family business. One of the study's most important findings is that next-generation leader acceptance of personal responsibility for their leadership behaviors and results serves as a mediator through which the nature of the family climate influences their leadership effectiveness and work engagement. This suggests that while the nature of family relationships may make it easier or more difficult, next-generation family leaders have ultimate control over the development of their leadership talent and the inspiration, enthusiasm, energy, and pride they feel when working in the family business.

Unraveling the phenotypic and genomic background of behavioral plasticity and temperament in North American Angus cattle.

BACKGROUND: Longitudinal records of temperament can be used for assessing behavioral plasticity, such as aptness to learn, memorize, or change behavioral responses based on affective state. In this study, we evaluated the phenotypic and genomic background of North American Angus cow temperament measured throughout their lifetime around the weaning season, including the development of a new indicator trait termed docility-based learning and behavioral plasticity. The analyses included 273,695 and 153,898 records for yearling (YT) and cow at weaning (CT) temperament, respectively, 723,248 animals in the pedigree, and 8784 genotyped animals. Both YT and CT were measured when the animal was loading into/exiting the chute. Moreover, CT was measured around the time in which the cow was separated from her calf. A random regression model fitting a first-order Legendre orthogonal polynomial was used to model the covariance structure of temperament and to assess the learning and behavioral plasticity (i.e., slope of the regression) of individual cows. This study provides, for the first time, a longitudinal perspective of the genetic and genomic mechanisms underlying temperament, learning, and behavioral plasticity in beef cattle. RESULTS: CT measured across years is heritable (0.38-0.53). Positive and strong genetic correlations (0.91-1.00) were observed among all CT age-group pairs and between CT and YT (0.84). Over 90% of the candidate genes identified overlapped among CT age-groups and the estimated effect of genomic markers located within important candidate genes changed over time. A small but significant genetic component was observed for learning and behavioral plasticity (heritability = 0.02 ± 0.002). Various candidate genes were identified, revealing the polygenic nature of the traits evaluated. The pathways and candidate genes identified are associated with steroid and glucocorticoid hormones, development delay, cognitive development, and behavioral changes in cattle and other species. CONCLUSIONS: Cow temperament is highly heritable and repeatable. The changes in temperament can be genetically improved by selecting animals with favorable learning and behavioral plasticity (i.e., habituation). Furthermore, the environment explains a large part of the variation in learning and behavioral plasticity, leading to opportunities to also improve the overall temperament by refining management practices. Moreover, behavioral plasticity offers opportunities to improve the long-term animal and handler welfare through habituation.

Genetic Modeling and Genomic Analyses of Yearling Temperament in American Angus Cattle and Its Relationship With Productive Efficiency and Resilience Traits.

Cattle temperament has been considered by farmers as a key breeding goal due to its relevance for cattlemen's safety, animal welfare, resilience, and longevity and its association with many economically important traits (e.g., production and meat quality). The definition of proper statistical models, accurate variance component estimates, and knowledge on the genetic background of the indicator trait evaluated are of great importance for accurately predicting the genetic merit of breeding animals. Therefore, 266,029 American Angus cattle with yearling temperament records (1-6 score) were used to evaluate statistical models and estimate variance components; investigate the association of sex and farm management with temperament; assess the weighted correlation of estimated breeding values for temperament and productive, reproductive efficiency and resilience traits; and perform a weighted single-step genome-wide association analysis using 69,559 animals genotyped for 54,609 single-nucleotide polymorphisms. Sex and extrinsic factors were significantly associated with temperament, including conception type, age of dam, birth season, and additional animal-human interactions. Similar results were observed among models including only the direct additive genetic effect and when adding other maternal effects. Estimated heritability of temperament was equal to 0.39 on the liability scale. Favorable genetic correlations were observed between temperament and other relevant traits, including growth, feed efficiency, meat quality, and reproductive traits. The highest approximated genetic correlations were observed between temperament and growth traits (weaning weight, 0.28; yearling weight, 0.28). Altogether, we identified 11 genomic regions, located across nine chromosomes including BTAX, explaining 3.33% of the total additive genetic variance. The candidate genes identified were enriched in pathways related to vision, which could be associated with reception of stimulus and/or cognitive abilities. This study encompasses large and diverse phenotypic, genomic, and pedigree datasets of US Angus cattle. Yearling temperament is a highly heritable and polygenic trait that can be improved through genetic selection. Direct selection for temperament is not expected to result in unfavorable responses on other relevant traits due to the favorable or low genetic correlations observed. In summary, this study contributes to a better understanding of the impact of maternal effects, extrinsic factors, and various genomic regions associated with yearling temperament in North American Angus cattle.

Survival of the first rather than the fittest in a Shewanella electrode biofilm.

For natural selection to operate there must exist heritable variation among individuals that affects their survival and reproduction. Among free-living microbes, where differences in growth rates largely define selection intensities, competitive exclusion is common. However, among surface attached communities, these dynamics become less predictable. If extreme circumstances were to dictate that a surface population is immortal and all offspring must emigrate, the offspring would be unable to contribute to the composition of the population. Meanwhile, the immortals, regardless of reproductive capacity, would remain unchanged in relative abundance. The normal cycle of birth, death, and competitive exclusion would be broken. We tested whether conditions required to set up this idealized scenario can be approximated in a microbial biofilm. Using two differentially-reproducing strains of Shewanella oneidensis grown on an anode as the sole terminal electron acceptor - a system in which metabolism is obligately tied to surface attachment - we found that selection against a slow-growing competitor is drastically reduced. This work furthers understanding of natural selection dynamics in sessile microbial communities, and provides a framework for designing stable microbial communities for industrial and experimental applications.

Impact of Censored or Penalized Data in the Genetic Evaluation of Two Longevity Indicator Traits Using Random Regression Models in North American Angus Cattle.

This study aimed to evaluate the impact of different proportions (i.e., 20%, 40%, 60% and 80%) of censored (CEN) or penalized (PEN) data in the prediction of breeding values (EBVs), genetic parameters, and computational efficiency for two longevity indicators (i.e., traditional and functional longevity; TL and FL, respectively). In addition, three different criteria were proposed for PEN: (1) assuming that all cows with censored records were culled one year after their last reported calving; (2) assuming that cows with censored records older than nine years were culled one year after their last reported calving, while censored (missing) records were kept for cows younger than nine years; and (3) assuming that cows with censored records older than nine years were culled one year after their last reported calving, while cows younger than nine years were culled two years after their last reported calving. All analyses were performed using random regression models based on fourth order Legendre orthogonal polynomials. The proportion of commonly selected animals and EBV correlations were calculated between the complete dataset (i.e., without censored or penalized data; COM) and all simulated proportions of CEN or PEN. The computational efficiency was evaluated based on the total computing time taken by each scenario to complete 150,000 Bayesian iterations. In summary, increasing the CEN proportion significantly (p-value < 0.05 by paired t-tests) decreased the heritability estimates for both TL and FL. When compared to CEN, PEN tended to yield heritabilities closer to COM, especially for FL. Moreover, similar heritability patterns were observed for all three penalization criteria. High proportions of commonly selected animals and EBV correlations were found between COM and CEN with 20% censored data (for both TL and FL), and COM and all levels of PEN (for FL). The proportions of commonly selected animals and EBV correlations were lower for PEN than CEN for TL, which suggests that the criteria used for PEN are not adequate for TL. Analyses using COM and CEN took longer to finish than PEN analyses. In addition, increasing the amount of censored records also tended to increase the computational time. A high proportion (>20%) of censored data has a negative impact in the genetic evaluation of longevity. The penalization criteria proposed in this study are useful for genetic evaluations of FL, but they are not recommended when analyzing TL.

A Systematic Review of Genomic Regions and Candidate Genes Underlying Behavioral Traits in Farmed Mammals and Their Link with Human Disorders.

The main objectives of this study were to perform a systematic review of genomic regions associated with various behavioral traits in the main farmed mammals and identify key candidate genes and potential causal mutations by contrasting the frequency of polymorphisms in cattle breeds with divergent behavioral traits (based on a subjective clustering approach). A total of 687 (cattle), 1391 (pigs), and 148 (sheep) genomic regions associated with 37 (cattle), 55 (pigs), and 22 (sheep) behavioral traits were identified in the literature. In total, 383, 317, and 15 genes overlap with genomic regions identified for cattle, pigs, and sheep, respectively. Six common genes (e.g., NR3C2, PITPNM3, RERG, SPNS3, U6, and ZFAT) were found for cattle and pigs. A combined gene-set of 634 human genes was produced through identified homologous genes. A total of 313 out of 634 genes have previously been associated with behavioral, mental, and neurologic disorders (e.g., anxiety and schizophrenia) in humans. Additionally, a total of 491 candidate genes had at least one statistically significant polymorphism (p-value < 0.05). Out of those, 110 genes were defined as having polymorphic regions differing in greater than 50% of exon regions. Therefore, conserved genomic regions controlling behavior were found across farmed mammal species and humans.

Haplotype-Based Single-Step GWAS for Yearling Temperament in American Angus Cattle.

Behavior is a complex trait and, therefore, understanding its genetic architecture is paramount for the development of effective breeding strategies. The objective of this study was to perform traditional and weighted single-step genome-wide association studies (ssGWAS and WssGWAS, respectively) for yearling temperament (YT) in North American Angus cattle using haplotypes. Approximately 266 K YT records and 70 K animals genotyped using a 50 K single nucleotide polymorphisms (SNP) panel were used. Linkage disequilibrium thresholds (LD) of 0.15, 0.50, and 0.80 were used to create the haploblocks, and the inclusion of non-LD-clustered SNPs (NCSNP) with the haplotypes in the genomic models was also evaluated. WssGWAS did not perform better than ssGWAS. Cattle YT was found to be a highly polygenic trait, with genes and quantitative trait loci (QTL) broadly distributed across the whole genome. Association studies using LD-based haplotypes should include NCSNPs and different LD thresholds to increase the likelihood of finding the relevant genomic regions affecting the trait of interest. The main candidate genes identified, i.e., ATXN10, ADAM10, VAX2, ATP6V1B1, CRISPLD1, CAPRIN1, FA2H, SPEF2, PLXNA1, and CACNA2D3, are involved in important biological processes and metabolic pathways related to behavioral traits, social interactions, and aggressiveness in cattle. Future studies should further investigate the role of these candidate genes.

Using Random Regression Models to Genetically Evaluate Functional Longevity Traits in North American Angus Cattle.

This study aimed to propose novel longevity indicators by comparing genetic parameters for traditional (TL; i.e., the cow's lifespan after the first calving) and functional (FL; i.e., how long the cow stayed in the herd while also calving; assuming no missing (FLa) or missing (FLb) records for unknown calving) longevity, considering different culling reasons (natural death, structural problems, disease, fertility, performance, and miscellaneous). Longevity definitions were evaluated from 2 to 15 years of age, using single- and multiple-trait Bayesian random regression models (RRM). The RRM fitting heterogenous residual variance and fourth order Legendre polynomials were considered as the optimal models for the majority of longevity indicators. The average heritability estimates over ages for FLb (from 0.08 to 0.25) were always higher than those for FLa (from 0.07 to 0.19), and higher or equal to the ones estimated for TL (from 0.07 to 0.23), considering the different culling reasons. The average genetic correlations estimated between ages were low to moderate (~0.40), for all longevity definitions and culling reasons. However, removing the extreme ages (i.e., 2 and >12 years) increased the average correlation between ages (from ~0.40 to >0.70). The genetic correlations estimated between culling reasons were low (0.12 and 0.20 on average, considering all ages and ages between 3 and 12 years old, respectively), indicating that longevity based on different culling reasons should be considered as different traits in the genetic evaluations. Higher average genetic correlations (estimated from 3 to 12 years old) were observed between TL and FLb (0.73) in comparison to TL and FLa (0.64), or FLa and FLb (0.65). Consequently, a higher average proportion of commonly-selected sires, for the top 1% sires, was also observed between TL and FLb (91.74%), compared to TL and FLa (59.68%), or FLa and FLb (61.01%). Higher prediction accuracies for the expected daughter performances (calculated based on the pedigree information) were obtained for FLb in comparison to TL and FLa. Our findings indicate that FLb is preferred for the genetic evaluation of longevity. In addition, it is recommended including multiple longevity traits based on different groups of culling reasons in a selection sub-index, as they are genetically-different traits. Genetic selection based on breeding values at the age of four years is expected to result in greater selection responses for increased longevity in North American Angus cattle.

Development of a genetic evaluation for hair shedding in American Angus cattle to improve thermotolerance.

BACKGROUND: Heat stress and fescue toxicosis caused by ingesting tall fescue infected with the endophytic fungus Epichloë coenophiala represent two of the most prevalent stressors to beef cattle in the United States and cost the beef industry millions of dollars each year. The rate at which a beef cow sheds her winter coat early in the summer is an indicator of adaptation to heat and an economically relevant trait in temperate or subtropical parts of the world. Furthermore, research suggests that early-summer hair shedding may reflect tolerance to fescue toxicosis, since vasoconstriction induced by fescue toxicosis limits the ability of an animal to shed its winter coat. Both heat stress and fescue toxicosis reduce profitability partly via indirect maternal effects on calf weaning weight. Here, we developed parameters for routine genetic evaluation of hair shedding score in American Angus cattle, and identified genomic loci associated with variation in hair shedding score via genome-wide association analysis (GWAA). RESULTS: Hair shedding score was moderately heritable (h2 = 0.34 to 0.40), with different repeatability estimates between cattle grazing versus not grazing endophyte-infected tall fescue. Our results suggest modestly negative genetic and phenotypic correlations between a dam's hair shedding score (lower score is earlier shedding) and the weaning weight of her calf, which is one metric of performance. Together, these results indicate that economic gains can be made by using hair shedding score breeding values to select for heat-tolerant cattle. GWAA identified 176 variants significant at FDR < 0.05. Functional enrichment analyses using genes that were located within 50 kb of these variants identified pathways involved in keratin formation, prolactin signalling, host-virus interaction, and other biological processes. CONCLUSIONS: This work contributes to a continuing trend in the development of genetic evaluations for environmental adaptation. Our results will aid beef cattle producers in selecting more sustainable and climate-adapted cattle, as well as enable the development of similar routine genetic evaluations in other breeds.

Genome-wide association studies and genomic prediction of breeding values for calving performance and body conformation traits in Holstein cattle.

BACKGROUND: Our aim was to identify genomic regions via genome-wide association studies (GWAS) to improve the predictability of genetic merit in Holsteins for 10 calving and 28 body conformation traits. Animals were genotyped using the Illumina Bovine 50 K BeadChip and imputed to the Illumina BovineHD BeadChip (HD). GWAS were performed on 601,717 real and imputed single nucleotide polymorphism (SNP) genotypes using a single-SNP mixed linear model on 4841 Holstein bulls with breeding value predictions and followed by gene identification and in silico functional analyses. The association results were further validated using five scenarios with different numbers of SNPs. RESULTS: Seven hundred and eighty-two SNPs were significantly associated with calving performance at a genome-wise false discovery rate (FDR) of 5%. Most of these significant SNPs were on chromosomes 18 (71.9%), 17 (7.4%), 5 (6.8%) and 7 (2.4%) and mapped to 675 genes, among which 142 included at least one significant SNP and 532 were nearby one (100 kbp). For body conformation traits, 607 SNPs were significant at a genome-wise FDR of 5% and most of them were located on chromosomes 5 (30%), 18 (27%), 20 (13%), 6 (6%), 7 (5%), 14 (5%) and 13 (3%). SNP enrichment functional analyses for calving traits at a FDR of 1% suggested potential biological processes including musculoskeletal movement, meiotic cell cycle, oocyte maturation and skeletal muscle contraction. Furthermore, pathway analyses suggested potential pathways associated with calving performance traits including tight junction, oxytocin signaling, and MAPK signaling (P < 0.10). The prediction ability of the 1206 significant SNPs was between 78 and 83% of the prediction ability of the BovineSNP50 SNPs for calving performance traits and between 35 and 79% for body conformation traits. CONCLUSIONS: Various SNPs that are significantly associated with calving performance are located within or nearby genes with potential roles in tight junction, oxytocin signaling, and MAPK signaling. Combining the significant SNPs or SNPs within or nearby gene(s) from the HD panel with the BovineSNP50 panel yielded a marginal increase in the accuracy of prediction of genomic estimated breeding values for all traits compared to the use of the BovineSNP50 panel alone.

Evaluation of blood metabolites reflects presence or absence of liver abscesses in beef cattle.

Liver abscesses constitute a prominent concern regarding animal health and profitability of the beef industry. Our objective was to evaluate potential biliary and blood indicators of liver abscesses. Twenty-nine beef bulls (initially averaging 356±70.5 kg and 253±30 days of age) were fed a high-concentrate diet during a performance test of 112 days, during which blood was collected at nine time points spaced 0.5-13 days apart within 56 days before slaughter. At the abattoir, blood and bile were collected and livers were inspected for liver abscesses. Results indicated that liver abscesses are associated with elevated levels of plasma cortisol and aspartate aminotransferase, and decreased levels of albumin, cholesterol and testosterone over the period before slaughter. Based on the blood samples collected during exsanguination, the presence of liver abscesses was associated with lower concentrations of thyroxine, albumin, cholesterol and alkaline phosphatase, and is suggested to be associated with lower blood carbon dioxide (P=0.08) and lower biliary cortisol metabolites (P=0.07). Albumin and cholesterol are established indicators of hepatic function and are consistently related to the presence of liver abscesses. Identifying blood parameters that predict liver abscesses has practical implications for cattle husbandry and for ensuring food safety.

Maternal nutrient restriction in mid-to-late gestation influences fetal mRNA expression in muscle tissues in beef cattle.

BACKGROUND: Manipulating maternal nutrition during specific periods of gestation can result in re-programming of fetal and post-natal development. In this experiment we investigated how a feed restriction of 85% compared with 140% of total metabolizable energy requirements, fed to cows during mid-to-late gestation, influences phenotypic development of fetuses and mRNA expression of growth (Insulin-Like Growth Factor family and Insulin Receptor (INSR)), myogenic (Myogenic Differentiation 1 (MYOD1), Myogenin (MYOG), Myocyte Enhancer Factor 2A (MEF2A), Serum Response Factor (SRF)) and adipogenic (Peroxisome Proliferator Activated Receptor Gamma (PPARG)) genes in fetal longissimus dorsi (LD) and semitendinosus (ST) muscle. DNA methylation of imprinted genes, Insulin Like Growth Factor 2 (IGF2) and Insulin Like Growth Factor 2 Receptor (IGF2R), and micro RNA (miRNA) expression, were also examined as potential consequences of poor maternal nutrition, but also potential regulators of altered gene expression patterns. RESULTS: While the nutrient restriction impacted dam body weight, no differences were observed in phenotypic fetal measurements (weight, crown-rump length, or thorax circumference). Interestingly, LD and ST muscles responded differently to the differential pre-natal nutrient levels. While LD muscle of restricted fetal calves had greater mRNA abundances for Insulin Like Growth Factor 1 and its receptor (IGF1 and IGF1R), IGF2R, INSR, MYOD1, MYOG, and PPARG, no significant differences were observed for gene expression in ST muscle. Similarly, feed restriction had a greater impact on the methylation level of IGF2 Differentially Methylated Region 2 (DMR2) in LD muscle as compared to ST muscle between treatment groups. A negative correlation existed between IGF2 mRNA expression and IGF2 DMR2 methylation level in both LD and ST muscles. Differential expression of miRNAs 1 and 133a were also detected in LD muscle. CONCLUSIONS: Our data suggests that a nutrient restriction of 85% as compared to 140% of total metabolizable energy requirements during the 2nd half of gestation can alter the expression of growth, myogenic and adipogenic genes in fetal muscle without apparent differences in fetal phenotype. It also appears that the impact of feed restriction varies between muscles suggesting a priority for nutrient partitioning depending on muscle function and/or fiber composition. Differences in the methylation level in IGF2, a well-known imprinted gene, as well as differences in miRNA expression, may be functional mechanisms that precede the differences in gene expression observed, and could lead to trans-generational epigenetic programming.

Genetic diversity of a New Zealand multi-breed sheep population and composite breeds' history revealed by a high-density SNP chip.

BACKGROUND: Knowledge about the genetic diversity of a population is a crucial parameter for the implementation of successful genomic selection and conservation of genetic resources. The aim of this research was to establish the scientific basis for the implementation of genomic selection in a composite Terminal sheep breeding scheme by providing consolidated linkage disequilibrium (LD) measures across SNP markers, estimating consistency of gametic phase between breed-groups, and assessing genetic diversity measures, such as effective population size (Ne), and population structure parameters, using a large number of animals (n = 14,845) genotyped with a high density SNP chip (606,006 markers). Information generated in this research will be useful for optimizing molecular breeding values predictions and managing the available genetic resources. RESULTS: Overall, as expected, levels of pairwise LD decreased with increasing distance between SNP pairs. The mean LD r2 between adjacent SNP was 0.26 ± 0.10. The most recent effective population size for all animals (687) and separately per breed-groups: Primera (974), Lamb Supreme (380), Texel (227) and Dual-Purpose (125) was quite variable. The genotyped animals were outbred or had an average low level of inbreeding. Consistency of gametic phase was higher than 0.94 for all breed pairs at the average distance between SNP on the chip (~4.74 kb). Moreover, there was not a clear separation between the breed-groups based on principal component analysis, suggesting that a mixed-breed training population for calculation of molecular breeding values would be beneficial. CONCLUSIONS: This study reports, for the first time, estimates of linkage disequilibrium, genetic diversity and population structure parameters from a genome-wide perspective in New Zealand Terminal Sire composite sheep breeds. The levels of linkage disequilibrium indicate that genomic selection could be implemented with the high density SNP panel. The moderate to high consistency of gametic phase between breed-groups and overlapping population structure support the pooling of the animals in a mixed training population for genomic predictions. In addition, the moderate to high Ne highlights the need to genotype and phenotype a large training population in order to capture most of the haplotype diversity and increase accuracies of genomic predictions. The results reported herein are a first step toward understanding the genomic architecture of a Terminal Sire composite sheep population and for the optimal implementation of genomic selection and genome-wide association studies in this sheep population.

Prediction of genomic breeding values for growth, carcass and meat quality traits in a multi-breed sheep population using a HD SNP chip.

BACKGROUND: New Zealand has some unique Terminal Sire composite sheep breeds, which were developed in the last three decades to meet commercial needs. These composite breeds were developed based on crossing various Terminal Sire and Maternal breeds and, therefore, present high genetic diversity compared to other sheep breeds. Their breeding programs are focused on improving carcass and meat quality traits. There is an interest from the industry to implement genomic selection in this population to increase the rates of genetic gain. Therefore, the main objectives of this study were to determine the accuracy of predicted genomic breeding values for various growth, carcass and meat quality traits using a HD SNP chip and to evaluate alternative genomic relationship matrices, validation designs and genomic prediction scenarios. A large multi-breed population (n = 14,845) was genotyped with the HD SNP chip (600 K) and phenotypes were collected for a variety of traits. RESULTS: The average observed accuracies (± SD) for traits measured in the live animal, carcass, and, meat quality traits ranged from 0.18 ± 0.07 to 0.33 ± 0.10, 0.28 ± 0.09 to 0.55 ± 0.05 and 0.21 ± 0.07 to 0.36 ± 0.08, respectively, depending on the scenario/method used in the genomic predictions. When accounting for population stratification by adjusting for 2, 4 or 6 principal components (PCs) the observed accuracies of molecular breeding values (mBVs) decreased or kept constant for all traits. The mBVs observed accuracies when fitting both G and A matrices were similar to fitting only G matrix. The lowest accuracies were observed for k-means cross-validation and forward validation performed within each k-means cluster. CONCLUSIONS: The accuracies observed in this study support the feasibility of genomic selection for growth, carcass and meat quality traits in New Zealand Terminal Sire breeds using the Ovine HD SNP chip. There was a clear advantage on using a mixed training population instead of performing analyzes per genomic clusters. In order to perform genomic predictions per breed group, genotyping more animals is recommended to increase the size of the training population within each group and the genetic relationship between training and validation populations. The different scenarios evaluated in this study will help geneticists and breeders to make wiser decisions in their breeding programs.

Assessing accuracy of imputation using different SNP panel densities in a multi-breed sheep population.

BACKGROUND: Genotype imputation is a key element of the implementation of genomic selection within the New Zealand sheep industry, but many factors can influence imputation accuracy. Our objective was to provide practical directions on the implementation of imputation strategies in a multi-breed sheep population genotyped with three single nucleotide polymorphism (SNP) panels: 5K, 50K and HD (600K SNPs). RESULTS: Imputation from 5K to HD was slightly better (0.6 %) than imputation from 5K to 50K. Two-step imputation from 5K to 50K and then from 50K to HD outperformed direct imputation from 5K to HD. A slight loss in imputation accuracy was observed when a large fixed reference population was used compared to a smaller within-breed reference (including all 50K genotypes on animals from different breeds excluding those in the validation set i.e. to be imputed), but only for a few animals across all imputation scenarios from 5K to 50K. However, a major gain in imputation accuracy for a large proportion of animals (purebred and crossbred), justified the use of a fixed and large reference dataset for all situations. This study also investigated the loss in imputation accuracy specifically for SNPs located at the ends of each chromosome, and showed that only chromosome 26 had an overall imputation (5K to 50K) accuracy for 100 SNPs at each end higher than 60 % (r2). Most of the chromosomes displayed reduced imputation accuracy at least at one of their ends. Prediction of imputation accuracy based on the relatedness of low-density genotypes to those of the reference dataset, before imputation (without running an imputation software) was also investigated. FIMPUTE V2.2 outperformed BEAGLE 3.3.2 across all imputation scenarios. CONCLUSIONS: Imputation accuracy in sheep breeds can be improved by following a set of recommendations on SNP panels, software, strategies of imputation (one- or two-step imputation), and choice of the animals to be genotyped using both high- and low-density SNP panels. We present a method that predicts imputation accuracy for individual animals at the low-density level, before running imputation, which can be used to restrict genomic prediction only to the animals that can be imputed with sufficient accuracy.

Genome-wide association for milk production and female fertility traits in Canadian dairy Holstein cattle.

BACKGROUND: Genome-wide association studies (GWAS) are a powerful tool for detecting genomic regions explaining variation in phenotype. The objectives of the present study were to identify or refine the positions of genomic regions affecting milk production, milk components and fertility traits in Canadian Holstein cattle, and to use these positions to identify genes and pathways that may influence these traits. RESULT: Several QTL regions were detected for milk production (MILK), fat production (FAT), protein production (PROT) and fat and protein deviation (FATD, PROTD respectively). The identified QTL regions for production traits (including milk production) support previous findings and some overlap with genes with known relevant biological functions identified in earlier studies such as DGAT1 and CPSF1. A significant region on chromosome 21 overlapping with the gene FAM181A and not previous linked to fertility in dairy cattle was identified for the calving to first service interval and days open. A functional enrichment analysis of the GWAS results yielded GO terms consistent with the specific phenotypes tested, for example GO terms GO:0007595 (lactation) and GO:0043627 (response to estrogen) for milk production (MILK), GO:0051057 (positive regulation of small GTPase mediated signal transduction) for fat production (FAT), GO:0040019 (positive regulation of embryonic development) for first service to calving interval (CTFS) and GO:0043268 (positive regulation of potassium ion transport) for days open (DO). In other cases the connection between the enriched GO terms and the traits were less clear, for example GO:0003279 (cardiac septum development) for FAT and GO:0030903 (notochord development) for DO trait. CONCLUSION: The chromosomal regions and enriched pathways identified in this study confirm several previous findings and highlight new regions and pathways that may contribute to variation in production or fertility traits in dairy cattle.

Genome-wide association for heifer reproduction and calf performance traits in beef cattle.

The aim of this study was to identify SNP markers that associate with variation in beef heifer reproduction and performance of their calves. A genome-wide association study was performed by means of the generalized quasi-likelihood score (GQLS) method using heifer genotypes from the BovineSNP50 BeadChip and estimated breeding values for pre-breeding body weight (PBW), pregnancy rate (PR), calving difficulty (CD), age at first calving (AFC), calf birth weight (BWT), calf weaning weight (WWT), and calf pre-weaning average daily gain (ADG). Data consisted of 785 replacement heifers from three Canadian research herds, namely Brandon Research Centre, Brandon, Manitoba, University of Alberta Roy Berg Kinsella Ranch, Kinsella, Alberta, and Lacombe Research Centre, Lacombe, Alberta. After applying a false discovery rate correction at a 5% significance level, a total of 4, 3, 3, 9, 6, 2, and 1 SNPs were significantly associated with PBW, PR, CD, AFC, BWT, WWT, and ADG, respectively. These SNPs were located on chromosomes 1, 5-7, 9, 13-16, 19-21, 24, 25, and 27-29. Chromosomes 1, 5, and 24 had SNPs with pleiotropic effects. New significant SNPs that impact functional traits were detected, many of which have not been previously reported. The results of this study support quantitative genetic studies related to the inheritance of these traits, and provides new knowledge regarding beef cattle quantitative trait loci effects. The identification of these SNPs provides a starting point to identify genes affecting heifer reproduction traits and performance of their calves (BWT, WWT, and ADG). They also contribute to a better understanding of the biology underlying these traits and will be potentially useful in marker- and genome-assisted selection and management.