Genome-wide association study of temperament and tenderness using different Bayesian approaches in a Nellore–Angus crossbred population

Genetic evaluation of crossbred Bos indicus cow temperament at parturition

Cow temperament at parturition may be mostly a measure of aggressiveness. The heritability of cow temperament at parturition in Bos taurus cows has been reported to be low. The objectives of this study were to estimate the heritability of cow temperament at parturition, conduct a genome-wide association analysis of cow temperament at the time of parturition, and estimate the correspondence of cow temperament at the time of parturition with cow productive performance and early-life temperament traits in Bos indicus crossbreds. Cow temperament was assessed from 1 to 5 indicating increasing levels of aggressiveness of cows (937 cows and 4,337 parturitions) from 2005 to 2022. Estimates of heritability and repeatability were 0.12 ± 0.024 and 0.24 ± 0.018. The estimates of proportion of phenotypic variance were 0.13 ± 0.019 and 0.02 ± 0.011 for permanent and maternal permanent environmental components, respectively. Estimates of heritability for maximum lifetime temperament score and proportions of temperament scores >1 were 0.18 ± 0.07 and 0.13 ± 0.072. Within cycles (generations), 2-yr-old cows had lower temperament score means than cows in most other age categories. There were low to moderate positive estimates of unadjusted correlation coefficients (r = 0.22 to 0.29; P < 0.05) of unadjusted temperament score with temperament measured on the same females when they were 8 mo old. There were low to moderate positive estimates of correlation coefficients (r = 0.09 to 0.37; P < 0.05) of unadjusted temperament score with calving rate, weaning rate, weaning weight per cow exposed, and weaning weight per 454 kg cow weight at weaning. Cows with the lowest temperament score had lower (P < 0.05) calving and weaning rate than cows in other temperament categories. Within 3 of 5 cycles, cows with the lowest temperament score (totally docile) had lower (P < 0.05) weaning weight per cow exposed than cows in other temperament categories. There were 2 SNP on BTA 4 associated with maximum lifetime temperament score (FDR < 0.05). The non-genetic influence of a cow's mother was documented in her own temperament measured at the time of calving; this may be a consequence of learned behavior. Less aggressiveness displayed by cows at the time of calving may be accompanied by lower reproductive and maternal performance.

Genetic Variants and Their Putative Effects on microRNA-Seed Sites: Characterization of the 3' Untranslated Region of Genes Associated with Temperament

The 3' untranslated region has an important role in gene regulation through microRNAs, and it has been estimated that microRNAs regulate up to 50% of coding genes in mammals. With the aim of allelic variant identification of 3' untranslated region microRNA seed sites, the 3' untranslated region was searched for seed sites of four temperament-associated genes (CACNG4, EXOC4, NRXN3, and SLC9A4). The microRNA seed sites were predicted in the four genes, and the CACNG4 gene had the greatest number with 12 predictions. To search for variants affecting the predicted microRNA seed sites, the four 3' untranslated regions were re-sequenced in a Brahman cattle population. Eleven single nucleotide polymorphisms were identified in the CACNG4, and eleven in the SLC9A4. Rs522648682:T>G of the CACNG4 gene was located at the predicted seed site for bta-miR-191. Rs522648682:T>G evidenced an association with both exit velocity (p = 0.0054) and temperament score (p = 0.0097). The genotype TT had a lower mean exit velocity (2.93 ± 0.4 m/s) compared with the TG and GG genotypes (3.91 ± 0.46 m/s and 3.67 ± 0.46 m/s, respectively). The allele associated with the temperamental phenotype antagonizes the seed site, disrupting the bta-miR-191 recognition. The G allele of CACNG4-rs522648682 has the potential to influence bovine temperament through a mechanism associated with unspecific recognition of bta-miR-191.

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.

Identification of Genomic Regions Associated with Concentrations of Milk Fat, Protein, Urea and Efficiency of Crude Protein Utilization in Grazing Dairy Cows

The objective of this study was to identify genomic regions associated with milk fat percentage (FP), crude protein percentage (CPP), urea concentration (MU) and efficiency of crude protein utilization (ECPU: ratio between crude protein yield in milk and dietary crude protein intake) using grazing, mixed-breed, dairy cows in New Zealand. Phenotypes from 634 Holstein Friesian, Jersey or crossbred cows were obtained from two herds at Massey University. A subset of 490 of these cows was genotyped using Bovine Illumina 50K SNP-chips. Two genome-wise association approaches were used, a single-locus model fitted to data from 490 cows and a single-step Bayes C model fitted to data from all 634 cows. The single-locus analysis was performed with the Efficient Mixed-Model Association eXpedited model as implemented in the SVS package. Single nucleotide polymorphisms (SNPs) with genome-wide association p-values ≤ 1.11 × 10⁻⁶ were considered as putative quantitative trait loci (QTL). The Bayes C analysis was performed with the JWAS package and 1-Mb genomic windows containing SNPs that explained > 0.37% of the genetic variance were considered as putative QTL. Candidate genes within 100 kb from the identified SNPs in single-locus GWAS or the 1-Mb windows were identified using gene ontology, as implemented in the Ensembl Genome Browser. The genes detected in association with FP (MGST1, DGAT1, CEBPD, SLC52A2, GPAT4, and ACOX3) and CPP (DGAT1, CSN1S1, GOSR2, HERC6, and IGF1R) were identified as candidates. Gene ontology revealed six novel candidate genes (GMDS, E2F7, SIAH1, SLC24A4, LGMN, and ASS1) significantly associated with MU whose functions were in protein catabolism, urea cycle, ion transportation and N excretion. One novel candidate gene was identified in association with ECPU (MAP3K1) that is involved in post-transcriptional modification of proteins. The findings should be validated using a larger population of New Zealand grazing dairy cows.

Adrenergic Tone as an Intermediary in the Temperament Syndrome Associated With Flight Speed in Beef Cattle

The temperament of farm animals can influence their resilience to everyday variations within the managed production environment and has been under strong direct and indirect selection during the course of domestication. A prominent objective measure used for assessing temperament in beef cattle is the behavioral flight response to release from confinement in a crush or chute. This behavioral measure, termed flight speed (also known as escape velocity) is associated with physiological processes including body temperature, feeding behavior, growth rate, carcass composition, immune function, and health outcomes. This review examines the functional links between this suite of traits and adrenergic activity of the sympathetic nervous system and the adrenomedullary hormonal system. It is suggested that flight speed is the behavioral aspect of an underlying “flightiness” temperament syndrome, and that elevated adrenergic tone in animals with a high level of flightiness (i.e., flighty animals) tunes physiological activities toward a sustained “fight or flight” defense profile that reduces productivity and the capacity to flourish within the production environment. Nonetheless, despite a common influence of adrenergic tone on this suite of traits, variation in each trait is also influenced by other regulatory pathways and by the capacity of tissues to respond to a range of modulators in addition to adrenergic stimuli. It is suggested that tuning by adrenergic tone is an example of homeorhetic regulation that can help account for the persistent expression of behavioral and somatic traits associated with the flight speed temperament syndrome across the life of the animal. At a population level, temperament may modulate ecological fit within and across generations in the face of environmental variability and change. Associations of flight speed with the psychological affective state of the animal, and implications for welfare are also considered. The review will help advance understanding of the developmental biology and physiological regulation of temperament syndromes.

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

Simple Summary

This study is a comprehensive review of genomic regions associated with animal behavior in farmed mammals (beef and dairy cattle, pigs, and sheep) which contributes to a better understanding of the biological mechanisms influencing the target indicator trait and to gene expression studies by suggesting genes likely controlling the trait, and it will be useful in optimizing genomic predictions of breeding values incorporating biological information. Behavioral mechanisms are complex traits, genetically controlled by multiple genes spread across the whole genome. The majority of the genes identified in cattle, pigs, and sheep in association with a plethora of behavioral measurements (e.g., temperament, terrain use, milking speed, tail biting, and sucking reflex) are likely controlling stimuli reception (e.g., olfactory), internal recognition of stimuli (e.g., neuroactive ligand–receptor interaction), and body response to a stimulus (e.g., blood pressure, fatty acidy metabolism, hormone signaling, and inflammatory pathways). Six genes were commonly identified between cattle and pigs. About half of the genes for behavior identified in farmed mammals were also identified in humans for behavioral, mental, and neuronal disorders. Our findings indicate that the majority of the genes identified are likely controlling animal behavioral outcomes because their biological functions as well as potentially differing allele frequencies between two breed groups (subjectively) clustered based on their temperament characteristics.

Abstract

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.

Breed Ancestry, Divergence, Admixture, and Selection Patterns of the Simbra Crossbreed

In this study, we evaluated an admixed South African Simbra crossbred population, as well as the Brahman (Indicine) and Simmental (Taurine) ancestor populations to understand their genetic architecture and detect genomic regions showing signatures of selection. Animals were genotyped using the Illumina BovineLD v2 BeadChip (7K). Genomic structure analysis confirmed that the South African Simbra cattle have an admixed genome, composed of 5/8 Taurine and 3/8 Indicine, ensuring that the Simbra genome maintains favorable traits from both breeds. Genomic regions that have been targeted by selection were detected using the linkage disequilibrium-based methods iHS and Rsb. These analyses identified 10 candidate regions that are potentially under strong positive selection, containing genes implicated in cattle health and production (e.g., TRIM63, KCNA10, NCAM1, SMIM5, MIER3, and SLC24A4). These adaptive alleles likely contribute to the biological and cellular functions determining phenotype in the Simbra hybrid cattle breed. Our data suggested that these alleles were introgressed from the breed's original indicine and taurine ancestors. The Simbra breed thus possesses derived parental alleles that combine the superior traits of the founder Brahman and Simmental breeds. These regions and genes might represent good targets for ad-hoc physiological studies, selection of breeding material and eventually even gene editing, for improved traits in modern cattle breeds. This study represents an important step toward developing and improving strategies for selection and population breeding to ultimately contribute meaningfully to the beef production industry.

Novel genes involved in the genetic architecture of temperament in Brahman cattle

Cattle temperament is a complex and economically relevant trait. The objective of this study was to identify genomic regions and genes associated with cattle temperament. From a Brahman cattle population of 1,370 animals evaluated for temperament traits (Exit velocity-EV, Pen Score-PS, Temperament Score-TS), two groups of temperament-contrasting animals were identified based on their EV-average values ±1/2 standard deviation (SD). To be considered in the calm group, the EV of females ranged between 0.16–1.82 m/s (n = 50) and the EV of males ranged between 0.4–1.56 m/s (n = 48). Females were classified as temperamental if their EV ranged between 3.13–7.66 m/s (n = 46) and males were classified as temperamental if their EV ranged between 3.05–10.83 m/s (n = 45). Selected animals were genotyped using a total of 139,376 SNPs (GGP-HD-150K), evaluated for their association with EV. The Genome-Wide Association analysis (GWAS) identified fourteen SNPs: rs135340276, rs134895560, rs110190635, rs42949831, rs135982573, rs109393235, rs109531929, rs135087545, rs41839733, rs42486577, rs136661522, rs110882543, rs110864071, rs109722627, (P<8.1E-05), nine of them were located on intergenic regions, harboring seventeen genes, of which only ACER3, VRK2, FANCL and SLCO3A1 were considered candidate associated with bovine temperament due to their reported biological functions. Five SNPs were located at introns of the NRXN3, EXOC4, CACNG4 and SLC9A4 genes. The indicated candidate genes are implicated in a wide range of behavioural phenotypes and complex cognitive functions. The association of the fourteen SNPs on bovine temperament traits (EV, PS and TS) was evaluated; all these SNPs were significant for EV; only some were associated with PS and TS. Fourteen SNPs were associated with EV which allowed the identification of twenty-one candidate genes for Brahman temperament. From a functional point of view, the five intronic SNPs identified in this study, are candidates to address control of bovine temperament, further investigation will probe their role in expression of this trait.

Genome-wide association study for stayability measures in Nellore-Angus crossbred cows

Beef cow stayability is a complex, economically important trait often used as an indicator of a cow's potential lifetime productivity. Stayability is defined as capability of a cow to maintain a perfect record up to 6 yr of age. This age is commonly cited as a financial break-even point, where initial costs of cow development and maintenance are recovered by her cumulative net income from yearly calf receipts. Later-maturing Bos indicus-Bos taurus crossbred cows may experience reproductive difficulty early in life but have a high potential for a long reproductive life span. It was the objective of this study to identify genetic variants associated with measures of beef cow stayability. A population of B. indicus-B. taurus crossbred cows (n = 305) from central Texas was used. Phenotypes for various measures of stayability to 6 yr of age were produced by artificially imposing five different culling criteria on data from the population. Cows were scored either as a 1 (indicating a perfect record through 6 yr) or a 0 (indicating failure at or before 6 yr), under each criterion. Genome-wide association studies (GWAS) were conducted for each criterion using univariate procedures and prefitting the fixed effect of cow contemporary group. SNP associations for two criteria surpassed the false discovery threshold of 0.15, when a cow was scored as 0 upon her first failure to wean a calf, regardless of reason, through 6 yr (criterion 2), and when a cow was scored as 0 upon her first failure to give birth to a calf, through 6 yr (criterion 3). Associated SNP were found on bovine chromosomes (BTA) 1, 2, 5, 9, 18, and 21 for criterion 2 and on BTA 1, 5, 11, 15, and 24 for criterion 3. A critical region on BTA 5: 43-50 Mb was identified for each criterion. Due to the similarities to prior work, the tendency for B. indicus cattle to experience reproductive difficulties early in life, and due to the large proportion of cows that left the herd at an early age under these criteria, these results suggest that the associations are likely driven by an early life trait such as age at puberty or rate of heifer development.

Network Analysis Reveals Putative Genes Affecting Meat Quality in Angus Cattle

Improvements in eating satisfaction will benefit consumers and should increase beef demand which is of interest to the beef industry. Tenderness, juiciness, and flavor are major determinants of the palatability of beef and are often used to reflect eating satisfaction. Carcass qualities are used as indicator traits for meat quality, with higher quality grade carcasses expected to relate to more tender and palatable meat. However, meat quality is a complex concept determined by many component traits making interpretation of genome-wide association studies (GWAS) on any one component challenging to interpret. Recent approaches combining traditional GWAS with gene network interactions theory could be more efficient in dissecting the genetic architecture of complex traits. Phenotypic measures of 23 traits reflecting carcass characteristics, components of meat quality, along with mineral and peptide concentrations were used along with Illumina 54k bovine SNP genotypes to derive an annotated gene network associated with meat quality in 2,110 Angus beef cattle. The efficient mixed model association (EMMAX) approach in combination with a genomic relationship matrix was used to directly estimate the associations between 54k SNP genotypes and each of the 23 component traits. Genomic correlated regions were identified by partial correlations which were further used along with an information theory algorithm to derive gene network clusters. Correlated SNP across 23 component traits were subjected to network scoring and visualization software to identify significant SNP. Significant pathways implicated in the meat quality complex through GO term enrichment analysis included angiogenesis, inflammation, transmembrane transporter activity, and receptor activity. These results suggest that network analysis using partial correlations and annotation of significant SNP can reveal the genetic architecture of complex traits and provide novel information regarding biological mechanisms and genes that lead to complex phenotypes, like meat quality, and the nutritional and healthfulness value of beef. Improvements in genome annotation and knowledge of gene function will contribute to more comprehensive analyses that will advance our ability to dissect the complex architecture of complex traits.

Accuracy of genomic breeding values for meat tenderness in Polled Nellore cattle

Zebu () cattle, mostly of the Nellore breed, comprise more than 80% of the beef cattle in Brazil, given their tolerance of the tropical climate and high resistance to ectoparasites. Despite their advantages for production in tropical environments, zebu cattle tend to produce tougher meat than Bos taurus breeds. Traditional genetic selection to improve meat tenderness is constrained by the difficulty and cost of phenotypic evaluation for meat quality. Therefore, genomic selection may be the best strategy to improve meat quality traits. This study was performed to compare the accuracies of different Bayesian regression models in predicting molecular breeding values for meat tenderness in Polled Nellore cattle. The data set was composed of Warner-Bratzler shear force (WBSF) of longissimus muscle from 205, 141, and 81 animals slaughtered in 2005, 2010, and 2012, respectively, which were selected and mated so as to create extreme segregation for WBSF. The animals were genotyped with either the Illumina BovineHD (HD; 777,000 from 90 samples) chip or the GeneSeek Genomic Profiler (GGP Indicus HD; 77,000 from 337 samples). The quality controls of SNP were Hard-Weinberg Proportion -value ≥ 0.1%, minor allele frequency > 1%, and call rate > 90%. The FImpute program was used for imputation from the GGP Indicus HD chip to the HD chip. The effect of each SNP was estimated using ridge regression, least absolute shrinkage and selection operator (LASSO), Bayes A, Bayes B, and Bayes Cπ methods. Different numbers of SNP were used, with 1, 2, 3, 4, 5, 7, 10, 20, 40, 60, 80, or 100% of the markers preselected based on their significance test (-value from genomewide association studies [GWAS]) or randomly sampled. The prediction accuracy was assessed by the correlation between genomic breeding value and the observed WBSF phenotype, using a leave-one-out cross-validation methodology. The prediction accuracies using all markers were all very similar for all models, ranging from 0.22 (Bayes Cπ) to 0.25 (Bayes B). When preselecting SNP based on GWAS results, the highest correlation (0.27) between WBSF and the genomic breeding value was achieved using the Bayesian LASSO model with 15,030 (3%) markers. Although this study used relatively few animals, the design of the segregating population ensured wide genetic variability for meat tenderness, which was important to achieve acceptable accuracy of genomic prediction. Although all models showed similar levels of prediction accuracy, some small advantages were observed with the Bayes B approach when higher numbers of markers were preselected based on their -values resulting from a GWAS analysis.

Genome-wide association study for behavior, type traits, and muscular development in Charolais beef cattle

Behavior, type traits, and muscular development are of interest for beef cattle breeding. Genome-wide association studies (GWAS) enable the identification of candidate genes, which enables gene-based selection and provides insight in the genetic architecture of these traits. The objective of the current study was to perform a GWAS for 3 behavior traits, 12 type traits, and muscular development in Charolais cattle. Behavior traits, including aggressiveness at parturition, aggressiveness during gestation period, and maternal care, were scored by farmers. Type traits, including udder conformation, teat, feet and legs, and locomotion, were scored by trained classifiers. Data used in the GWAS consisted of 3,274 cows with phenotypic records and genotyping information for 44,930 SNP. When SNP had a false discovery rate (FDR) smaller than 0.05, they were referred to as significant. When SNP had a FDR between 0.05 and 0.20, they were referred to as suggestive. Four significant and 12 suggestive regions were detected for aggressiveness during gestation, maternal care, udder balance, teat thinness, teat length, foot angle, foot depth, and locomotion. These 4 significant and 12 suggestive regions were not supported by other significant SNP in close proximity. No SNP with major effects were detected for behavior and type traits, and SNP associations for these traits were spread across the genome, suggesting that behavior and type traits were influenced by many genes, each explaining a small part of genetic variance. The GWAS identified 1 region on chromosome 2 significantly associated with muscular development, which included the myostatin gene (), which is known to affect muscularity. No other regions associated with muscular development were found. Results showed that the myostatin region associated with muscular development had pleiotropic effects on udder volume, teat thinness, rear leg, and leg angle.

Crossbred steer temperament as yearlings and whole genome association of steer temperament as yearlings and calf temperament post-weaning

cattle often have the reputation for a poor or dangerous temperament. Identification of genomic regions that associate with temperament of such cattle may be useful for genetic improvement strategies. The objectives of this study were to evaluate subjective temperament scores (1 to 9; higher scores indicated more unfavorable temperament) for aggressiveness, nervousness, flightiness, gregariousness, and overall temperament of one-half steers in feedlot conditions at 1 yr of age and compare those scores of those steers when evaluated approximately 1 mo postweaning, and conduct whole genome association analyses using SNP markers and the temperament traits of those steers at 1 yr of age and for temperament traits of all calves at weaning. Contemporary groups ( < 0.001) were steers born in the same year and season, and fed in the same feedlot pen. Aggressiveness of steers at 1 yr of age was not associated with aggressiveness at weaning (linear regression coefficient did not differ from 0; = 0.96), but regressions of all other yearling scores of steers on the scores at weaning were positive (coefficients ranged from 0.26 ± 0.04 to 0.32 ± 0.04; < 0.001). Estimates of Pearson correlation coefficients (using unadjusted values and residual values) of the different traits measured at 1 yr of age were large ( > 0.63; < 0.008) except for aggressiveness with nervousness, flightiness, or gregariousness, which did not differ from 0 ( > 0.1). Five SNP on BTA 1, 24, and 29 had suggestive associations (0.17 < [adjusted for FDR] < 0.24) with aggressiveness, nervousness, or flightiness at evaluation postweaning and 13 SNP on 11 chromosomes had suggestive associations (0.07 < [adjusted for FDR] < 0.24) with aggressiveness, nervousness, flightiness, or overall temperament score of steers at 1 yr of age. Genes close to these loci with roles in neural systems of various organisms included synaptotagmin 4 (BTA 24), FAT atypical cadhedrin 3 (BTA 29), tubulin tyrosine ligase-like 1 (BTA 5), spermatogenesis associated 17 (BTA 16), stanniocalcin 2 (BTA 20), and GABA receptor γ 3 (BTA 21).

Association of SNPs in dopamine and serotonin pathway genes and their interacting genes with temperament traits in Charolais cows

Cattle temperament is a complex trait, and molecular studies aimed at defining this trait are scarce. We used an interaction networks approach to identify new genes (interacting genes) and to estimate their effects and those of 19 dopamine- and serotonin-related genes on the temperament traits of Charolais cattle. The genes proopiomelanocortin (POMC), neuropeptide Y (NPY), solute carrier family 18, member 2 (SLC18A2) and FBJ murine osteosarcoma viral oncogene homologue (FOSFBJ) were identified as new candidates. Their potential to be associated with temperament was estimated according to their reported biological activities, which included interactions with neural activity, receptor function, targeting or synthesis of neurotransmitters and association with behaviour. Pen score (PS) and exit velocity (EV) measures were determined from 412 Charolais cows to calculate their temperament score (TS). Based on the TS, calm (n = 55; TS, 1.09 ± 0.33) and temperamental (n = 58; TS, 2.27 ± 0.639) cows were selected and genotyped using a 248 single-nucleotide variation (SNV) panel. Of the 248 variations in the panel, only 151 were confirmed to be polymorphic (single-nucleotide polymorphisms; SNPs) in the tested population. Single-marker association analyses between genotypes and temperament measures (EV, PS and/or TS) indicated significant associations of six SNPs from four candidate genes. The markers rs109576799 and rs43696138, located in the DRD3 and HTR2A genes, respectively, were significantly associated with both EV and TS traits. Four markers, rs110365063 and rs137756569 from the POMC gene and rs110365063 and rs135155082 located in SLC18A2 and DRD2, respectively, were associated with PS. The variant rs110365063 located in bovine SLC18A2 causes a change in the amino acid sequence from Ala to Thr. Further studies are needed to confirm the association of genetic profile with cattle temperament; however, our study represents important progress in understanding the regulation of cattle temperament by different genes with divergent functions.

Genome-Wide Association Study between Single Nucleotide Polymorphisms and Flight Speed in Nellore Cattle

Introduction

Cattle temperament is an important factor that affects the profitability of beef cattle enterprises, due to its relationship with productivity traits, animal welfare and labor safety. Temperament is a complex phenotype often assessed by measuring a series of behavioral traits, which result from the effects of multiple environmental and genetic factors, and their interactions. The aims of this study were to perform a genome-wide association study and detect genomic regions, potential candidate genes and their biological mechanisms underlying temperament, measured by flight speed (FS) test in Nellore cattle.

Materials and Methods

The genome-wide association study (GWAS) was performed using a single-step procedure (ssGBLUP) which combined simultaneously all 16,600 phenotypes from genotyped and non-genotyped animals, full pedigree information of 162,645 animals and 1,384 genotyped animals in one step. The animals were genotyped with High Density Bovine SNP BeadChip which contains 777,962 SNP markers. After quality control (QC) a total of 455,374 SNPs remained.

Results

Heritability estimated for FS was 0.21 ± 0.02. Consecutive SNPs explaining 1% or more of the total additive genetic variance were considered as windows associated with FS. Nine candidate regions located on eight different Bos taurus chromosomes (BTA) (1 at 73 Mb, 2 at 65 Mb, 5 at 22 Mb and 119 Mb, 9 at 98 Mb, 11 at 67 Mb, 15 at 16 Mb, 17 at 63 Kb, and 26 at 47 Mb) were identified. The candidate genes identified in these regions were NCKAP5 (BTA2), PARK2 (BTA9), ANTXR1 (BTA11), GUCY1A2 (BTA15), CPE (BTA17) and DOCK1 (BTA26). Among these genes PARK2, GUCY1A2, CPE and DOCK1 are related to dopaminergic system, memory formation, biosynthesis of peptide hormone and neurotransmitter and brain development, respectively.

Conclusions

Our findings allowed us to identify nine genomic regions (SNP windows) associated with beef cattle temperament, measured by FS test. Within these windows, six promising candidate genes and their biological functions were identified. These results may contribute to a better comprehension into the genetic control of temperament expression in Nellore cattle.

Multi-breed and multi-trait co-association analysis of meat tenderness and other meat quality traits in three French beef cattle breeds

Background

Studies to identify markers associated with beef tenderness have focused on Warner–Bratzler shear force (WBSF) but the interplay between the genes associated with WBSF has not been explored. We used the association weight matrix (AWM), a systems biology approach, to identify a set of interacting genes that are co-associated with tenderness and other meat quality traits, and shared across the Charolaise, Limousine and Blonde d’Aquitaine beef cattle breeds.

Results

Genome-wide association studies were performed using ~500K single nucleotide polymorphisms (SNPs) and 17 phenotypes measured on more than 1000 animals for each breed. First, this multi-trait approach was applied separately for each breed across 17 phenotypes and second, between- and across-breed comparisons at the AWM and functional levels were performed. Genetic heterogeneity was observed, and most of the variants that were associated with WBSF segregated within rather than across breeds. We identified 206 common candidate genes associated with WBSF across the three breeds. SNPs in these common genes explained between 28 and 30 % of the phenotypic variance for WBSF. A reduced number of common SNPs mapping to the 206 common genes were identified, suggesting that different mutations may target the same genes in a breed-specific manner. Therefore, it is likely that, depending on allele frequencies and linkage disequilibrium patterns, a SNP that is identified for one breed may not be informative for another unrelated breed. Well-known candidate genes affecting beef tenderness were identified. In addition, some of the 206 common genes are located within previously reported quantitative trait loci for WBSF in several cattle breeds. Moreover, the multi-breed co-association analysis detected new candidate genes, regulators and metabolic pathways that are likely involved in the determination of meat tenderness and other meat quality traits in beef cattle.

Conclusions

Our results suggest that systems biology approaches that explore associations of correlated traits increase statistical power to identify candidate genes beyond the one-dimensional approach. Further studies on the 206 common genes, their pathways, regulators and interactions will expand our knowledge on the molecular basis of meat tenderness and could lead to the discovery of functional mutations useful for genomic selection in a multi-breed beef cattle context.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-016-0216-y) contains supplementary material, which is available to authorized users.

Detection of genetic variants affecting cattle behaviour and their impact on milk production: a genome-wide association study

Behaviour traits of cattle have been reported to affect important production traits, such as meat quality and milk performance as well as reproduction and health. Genetic predisposition is, together with environmental stimuli, undoubtedly involved in the development of behaviour phenotypes. Underlying molecular mechanisms affecting behaviour in general and behaviour and productions traits in particular still have to be studied in detail. Therefore, we performed a genome-wide association study in an F2 Charolais × German Holstein cross-breed population to identify genetic variants that affect behaviour-related traits assessed in an open-field and novel-object test and analysed their putative impact on milk performance. Of 37,201 tested single nucleotide polymorphism (SNPs), four showed a genome-wide and 37 a chromosome-wide significant association with behaviour traits assessed in both tests. Nine of the SNPs that were associated with behaviour traits likewise showed a nominal significant association with milk performance traits. On chromosomes 14 and 29, six SNPs were identified to be associated with exploratory behaviour and inactivity during the novel-object test as well as with milk yield traits. Least squares means for behaviour and milk performance traits for these SNPs revealed that genotypes associated with higher inactivity and less exploratory behaviour promote higher milk yields. Whether these results are due to molecular mechanisms simultaneously affecting behaviour and milk performance or due to a behaviour predisposition, which causes indirect effects on milk performance by influencing individual reactivity, needs further investigation.

Temperament type specific metabolite profiles of the prefrontal cortex and serum in cattle

In the past decade the number of studies investigating temperament in farm animals has increased greatly because temperament has been shown not only to affect handling but also reproduction, health and economically important production traits. However, molecular pathways underlying temperament and molecular pathways linking temperament to production traits, health and reproduction have yet to be studied in full detail. Here we report the results of metabolite profiling of the prefrontal cortex and serum of cattle with distinct temperament types that were performed to further explore their molecular divergence in the response to the slaughter procedure and to identify new targets for further research of cattle temperament. By performing an untargeted comprehensive metabolite profiling, 627 and 1097 metabolite features comprising 235 and 328 metabolites could be detected in the prefrontal cortex and serum, respectively. In total, 54 prefrontal cortex and 51 serum metabolite features were indicated to have a high relevance in the classification of temperament types by a sparse partial least square discriminant analysis. A clear discrimination between fearful/neophobic-alert, interested-stressed, subdued/uninterested-calm and outgoing/neophilic-alert temperament types could be observed based on the abundance of the identified relevant prefrontal cortex and serum metabolites. Metabolites with high relevance in the classification of temperament types revealed that the main differences between temperament types in the response to the slaughter procedure were related to the abundance of glycerophospholipids, fatty acyls and sterol lipids. Differences in the abundance of metabolites related to C21 steroid metabolism and oxidative stress indicated that the differences in the metabolite profiles of the four extreme temperament types could be the result of a temperament type specific regulation of molecular pathways that are known to be involved in the stress and fear response.

Comparison of breeding value prediction for two traits in a Nellore-Angus crossbred population using different Bayesian modeling methodologies

The objectives of this study were to 1) compare four models for breeding value prediction using genomic or pedigree information and 2) evaluate the impact of fixed effects that account for family structure. Comparisons were made in a Nellore-Angus population comprising F₂, F₃ and half-siblings to embryo transfer F₂ calves with records for overall temperament at weaning (TEMP; n = 769) and Warner-Bratzler shear force (WBSF; n = 387). After quality control, there were 34,913 whole genome SNP markers remaining. Bayesian methods employed were BayesB (π̃ = 0.995 or 0.997 for WBSF or TEMP, respectively) and BayesC (π = 0 and π̃), where π̃ is the ideal proportion of markers not included. Direct genomic values (DGV) from single trait Bayesian analyses were compared to conventional pedigree-based animal model breeding values. Numerically, BayesC procedures (using π̃) had the highest accuracy of all models for WBSF and TEMP (ρ̂_gĝ = 0.843 and 0.923, respectively), but BayesB had the least bias (regression of performance on prediction closest to 1, β̂_y,x = 2.886 and 1.755, respectively). Accounting for family structure decreased accuracy and increased bias in prediction of DGV indicating a detrimental impact when used in these prediction methods that simultaneously fit many markers.

Genetic selection for temperament traits in dairy and beef cattle

Animal temperament can be defined as a response to environmental or social stimuli. There are a number of temperament traits in cattle that contribute to their welfare, including their response to handling or milking, response to challenge such as human approach or intervention at calving, and response to conspecifics. In a number of these areas, the genetic basis of the trait has been studied. Heritabilities have been estimated and in some cases quantitative trait loci (QTL) have been identified. The variation is sometimes considerable and moderate heritabilities have been found for the major handling temperament traits, making them amenable to selection. Studies have also investigated the correlations between temperament and other traits, such as productivity and meat quality. Despite this, there are relatively few examples of temperament traits being used in selection programmes. Most often, animals are screened for aggression or excessive fear during handling or milking, with extreme animals being culled, or EBVs for temperament are estimated, but these traits are not commonly included routinely in selection indices, despite there being economic, welfare and human safety drivers for their. There may be a number of constraints and barriers. For some traits and breeds, there may be difficulties in collecting behavioral data on sufficiently large populations of animals to estimate genetic parameters. Most selection indices require estimates of economic values, and it is often difficult to assign an economic value to a temperament trait. The effects of selection primarily for productivity traits on temperament and welfare are discussed. Future opportunities include automated data collection methods and the wider use of genomic information in selection.