Analysis of nonsynonymous SNPs in candidate genes that influence bovine temperament and evaluation of their effect in Brahman cattle

BACKGROUND

Temperament is an important production trait in cattle and multiple strategies had been developed to generate molecular markers to assist animal selection. As nonsynonymous single nucleotide polymorphisms are markers with the potential to affect gene functions, they could be useful to predict phenotypic effects. Genetic selection of less stress-responsive, temperamental animals is desirable from an economic and welfare point of view.

METHODS AND RESULTS

Two nonsynonymous single nucleotide polymorphisms identified in HTR1B and SLC18A2 candidate genes for temperament were analyzed in silico to determine their effects on protein structure. Those nsSNPs allowing changes in proteins were selected for a temperament association analysis in a Brahman population. Transversion effects on protein structure were evaluated in silico for each amino acid change model, revealing structural changes in the proteins of the HTR1B and SLC18A2 genes. The selected nsSNPs were genotyped in a Brahman population (n = 138), and their genotypic effects on three temperament traits were analyzed: exit velocity, pen score, and temperament score. Only the SNP rs209984404-HTR1B (C/A) showed a significant association (P = 0.0144) with pen score. The heterozygous genotype showed a pen score value 1.17 points lower than that of the homozygous CC genotype.

CONCLUSION

The results showed that in silico analysis could direct the selection of nsSNPs with the potential to change the protein. Non-synonymous single nucleotide polymorphisms causing structural changes and reduced protein stability were identified. Only rs209984404-HTR1B shows that the allele affecting protein stability was associated with the genotype linked to docility in cattle.

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.

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.

Large-Scale Phenotyping of Livestock Welfare in Commercial Production Systems: A New Frontier in Animal Breeding

Genomic breeding programs have been paramount in improving the rates of genetic progress of productive efficiency traits in livestock. Such improvement has been accompanied by the intensification of production systems, use of a wider range of precision technologies in routine management practices, and high-throughput phenotyping. Simultaneously, a greater public awareness of animal welfare has influenced livestock producers to place more emphasis on welfare relative to production traits. Therefore, management practices and breeding technologies in livestock have been developed in recent years to enhance animal welfare. In particular, genomic selection can be used to improve livestock social behavior, resilience to disease and other stress factors, and ease habituation to production system changes. The main requirements for including novel behavioral and welfare traits in genomic breeding schemes are: (1) to identify traits that represent the biological mechanisms of the industry breeding goals; (2) the availability of individual phenotypic records measured on a large number of animals (ideally with genomic information); (3) the derived traits are heritable, biologically meaningful, repeatable, and (ideally) not highly correlated with other traits already included in the selection indexes; and (4) genomic information is available for a large number of individuals (or genetically close individuals) with phenotypic records. In this review, we (1) describe a potential route for development of novel welfare indicator traits (using ideal phenotypes) for both genetic and genomic selection schemes; (2) summarize key indicator variables of livestock behavior and welfare, including a detailed assessment of thermal stress in livestock; (3) describe the primary statistical and bioinformatic methods available for large-scale data analyses of animal welfare; and (4) identify major advancements, challenges, and opportunities to generate high-throughput and large-scale datasets to enable genetic and genomic selection for improved welfare in livestock. A wide variety of novel welfare indicator traits can be derived from information captured by modern technology such as sensors, automatic feeding systems, milking robots, activity monitors, video cameras, and indirect biomarkers at the cellular and physiological levels. The development of novel traits coupled with genomic selection schemes for improved welfare in livestock can be feasible and optimized based on recently developed (or developing) technologies. Efficient implementation of genetic and genomic selection for improved animal welfare also requires the integration of a multitude of scientific fields such as cell and molecular biology, neuroscience, immunology, stress physiology, computer science, engineering, quantitative genomics, and bioinformatics.

Use of random regression to estimate genetic parameters of temperament across an age continuum in a crossbred cattle population

The objective was to estimate genetic parameters of temperament in beef cattle across an age continuum. The population consisted predominantly of Brahman-British crossbred cattle. Temperament was quantified by 1) pen score (PS), the reaction of a calf to a single experienced evaluator on a scale of 1 to 5 (1 = calm, 5 = excitable); 2) exit velocity (EV), the rate (m/s) at which a calf traveled 1.83 m upon exiting a squeeze chute; and 3) temperament score (TS), the numerical average of PS and EV. Covariates included days of age and proportion of Bos indicus in the calf and dam. Random regression models included the fixed effects determined from the repeated measures models, except for calf age. Likelihood ratio tests were used to determine the most appropriate random structures. In repeated measures models, the proportion of B. indicus in the calf was positively related with each calf temperament trait (0.41 ± 0.20, 0.85 ± 0.21, and 0.57 ± 0.18 for PS, EV, and TS, respectively; P < 0.01). There was an effect of contemporary group (combinations of season, year of birth, and management group) and dam age (P < 0.001) in all models. From repeated records analyses, estimates of heritability (h2) were 0.34 ± 0.04, 0.31 ± 0.04, and 0.39 ± 0.04, while estimates of permanent environmental variance as a proportion of the phenotypic variance (c2) were 0.30 ± 0.04, 0.31 ± 0.03, and 0.34 ± 0.04 for PS, EV, and TS, respectively. Quadratic additive genetic random regressions on Legendre polynomials of age were significant for all traits. Quadratic permanent environmental random regressions were significant for PS and TS, but linear permanent environmental random regressions were significant for EV. Random regression results suggested that these components change across the age dimension of these data. There appeared to be an increasing influence of permanent environmental effects and decreasing influence of additive genetic effects corresponding to increasing calf age for EV, and to a lesser extent for TS. Inherited temperament may be overcome by accumulating environmental stimuli with increases in age, especially after weaning.

Genetic relationships between temperament of calves at auction and carcass traits in Japanese Black cattle

Correlations of calves' temperament with carcass traits were estimated to clarify the genetic relationships between them in Japanese Black cattle. The temperament records for 3128 calves during auction at a calf market were scored on a scale of 1 (calm) to 5 (nervous) as temperament score (TS), and the TS were divided into two groups (TSG): TS 1 and 2 comprised TSG 1, and 3 to 5 constituted TSG 2. Carcass data were obtained from 33 552 fattened cattle. A threshold animal model was used for analyzing the underlying liability for TSG, whereas a linear one was used for TS and carcass traits. The heritability estimates for TS and TSG were 0.12 and 0.11, respectively. On the other hand, moderate to high heritability estimates were obtained for carcass traits (0.40 to 0.68). The temperament scores were negatively correlated with carcass weight, rib thickness and subcutaneous fat thickness (-0.13 to -0.59). In contrast, weak to moderate positive correlations were found between the temperament scores and rib eye area or yield estimate (0.16 to 0.45). The temperament scores and beef marbling score had no correlation. These results showed that it is possible to improve temperament and carcass traits simultaneously.

Prenatal transportation stress alters temperament and serum cortisol concentrations in suckling Brahman calves

This experiment examined the relationship between prenatal stress and subsequent calf temperament through weaning. The prenatal stressor used was repeated transportation of pregnant Brahman cows for 2 h at 60 ± 5, 80 ± 5, 100 ± 5, 120 ± 5, and 140 ± 5 d of gestation. Prenatally stressed calves ( = 41) were compared with controls ( = 44; dams did not undergo transportation during pregnancy) from 2 wk of age until weaning (average age at weaning = 174.8 ± 1.3 d). Temperament was defined by pen score (PS; 1 = calm and 5 = excitable), exit velocity (EV; m/sec), and temperament score (TS; (PS + EV)/2) and was recorded for each calf on d -168, -140, -112, -84, -56, -28, and 0 relative to weaning (d 0 = weaning). Cortisol concentrations were determined in serum samples obtained on d -168, -140, -28, and 0 relative to weaning. Birth weight and weaning weight were not different between treatment groups ( > 0.1). Pen score was greater ( = 0.03) in prenatally stressed calves (2.84 ± 0.21) relative to controls (2.31 ± 0.21). Exit velocity was greater ( < 0.01) in prenatally stressed calves (2.1 ± 0.14 m/sec) than in controls (1.61 ± 0.14 m/sec). Exit velocity was affected by a treatment × calf sex interaction ( = 0.04) and was greater in prenatally stressed females. Exit velocity was also affected by day ( < 0.0001). Temperament score was greater ( = 0.01) in prenatally stressed calves (2.45 ± 0.16) than in controls (1.95 ± 0.16). Temperament score was affected by day ( < 0.01). Basal cortisol concentrations were greater ( = 0.04) in prenatally stressed calves (15.87 ± 1.04 ng/mL) than in controls (13.42 ± 1.03 ng/mL). Basal cortisol concentrations were greater ( < 0.01) in females (16.61 ± 1.06 ng/mL) than in males (12.68 ± 1.02 ng/mL). Cortisol concentrations were positively correlated ( < 0.01) with PS ( = 0.55, < 0.01), EV ( = 0.4, < 0.01), and TS ( = 0.55, < 0.01). Overall, suckling Brahman calves that were prenatally stressed were more temperamental and had greater circulating serum concentrations of cortisol than control calves.

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.

Genetic parameters for large-scale behavior traits and type traits in Charolais beef cows

In the last decades, beef cattle breeding mainly focused on improving production and reproduction traits. Nowadays, there is a growing interest to include behavior and type traits in the breeding goal. There is an interest in behavior as it is associated with human safety and workability and in type traits as they might be associated with longevity of cows. The objective of the current study was to estimate the heritability for behavior and type traits in Charolais and to estimate the genetic correlations among these traits. Behavior traits, including aggressiveness at parturition, aggressiveness during gestation period, and maternal care, were scored by farmers using an on-farm recording system to enable large-scale collection of phenotypes. Type traits, including udder traits ( = 3), teat traits (3), feet and leg traits (5), and locomotion (1), were scored by 10 trained classifiers. Data was available on 6,649 cows in parity 1 to 12 and located in 380 herds. Results showed that differences between herds explained up to 23% of the total phenotypic variance in behavior traits. This might be due to differences in management or to consistent differences in scoring between farmers. Aggressiveness at parturition had higher heritability (0.19) and higher genetic coefficient of variation (CV = 11%) than aggressiveness during gestation ( = 0.06 and CV = 4%) and maternal care ( = 0.02 and CV = 2%). Heritabilities for udder traits (0.14 to 0.20) and teat traits (0.17 to 0.35) were higher than for feet and leg traits (0.02 to 0.19). Genetic coefficients of variation for udder and teat traits were also higher (up to 21%) than for feet and leg traits (up to 11%). Strong genetic correlations were found between behavior traits (with absolute values from 0.71 to 0.98). The genetic correlations indicate that it is difficulty to simultaneously improve maternal care and reduce aggressiveness. We concluded that there are good opportunities to implement selection for improved udder and teat traits and against aggressiveness at parturition using a simple on-farm recording system of behavior.

Behavioural linear standardized scoring system of the Lidia cattle breed by testing in herd: estimation of genetic parameters

Docility is very important for cattle production, and many behavioural tests to measure this trait have been developed. However, very few objective behavioural tests to measure the opposite approach 'aggressive behaviour' have been described. Therefore, the aim of this work was to validate in the Lidia cattle breed a behavioural linear standardized scoring system that measure the aggressiveness and enable genetic analysis of behavioural traits expressing fearless and fighting ability. Reproducibility and repeatability measures were calculated for the 12 linear traits of this scoring system to assess its accuracy, and ranged from 85.3 and 94.2%, and from 66.7 to 97.9%, respectively. Genetic parameters were estimated using an animal model with a Bayesian approach. A total of 1202 behavioural records were used. The pedigree matrix contained 5001 individuals. Heritability values (with standard deviations) ranged between 0.13 (0.04) (Falls of the bull) and 0.41 (0.08) (Speed of approach to horse). Genetic correlations varied from 0.01 (0.07) to 0.90 (0.13). Finally, an exploratory factor analysis using the genetic correlation matrix was calculated. Three main factors were retained to describe the traditional genetic indexes aggressiveness, strength and mobility.

Genetic evaluation of aspects of temperament in Nellore-Angus calves

The objective of this work was to estimate heritability of each of 5 subjectively measured aspects of temperament of cattle and the genetic correlations of pairs of those traits. From 2003 to 2013, Nellore-Angus F2 and F3 calves (n = 1,816) were evaluated for aspects of temperament at an average 259 d of age, which was approximately 2 mo after weaning. Calves were separated from a group and subjectively scored from 1 (calm, good temperament) to 9 (wild, poor temperament) for aggressiveness (willingness to hit an evaluator), nervousness, flightiness, gregariousness (willingness to separate from the group), and a distinct overall score by 4 evaluators. Data were analyzed using threshold and linear models with additive genetic random effects. Two-trait animal models (nonthreshold) included the additive genetic covariance for pairs of traits and were used to estimate additive genetic correlations. Contemporary groups (n = 104) represented calves penned together for evaluation on given evaluation days. Heifers had greater (worse) means for all traits than steers (P < 0.05). The regression of score on age in days was included in final models for flightiness (P = 0.05; -0.006 ± 0.003) and gregariousness (P = 0.025; -0.007 ± 0.003). Estimates of heritability were large (0.51, 0.4, 0.45, 0.49, and 0.47 for aggressiveness, nervousness, flightiness, gregariousness, and overall temperament, respectively; SE = 0.07 for each). The ability to use this methodology to distinctly separate different aspects of calf temperament appeared to be limited, as estimates of additive genetic correlations were near unity for all pairs of traits; estimates of phenotypic correlation ranged from 0.88 ± 0.01 to 0.99 ± 0.002 for pairs of traits. Distinct subsequent analyses indicated a significant negative relationship of 4 of the various temperament scores with weight at weaning (regression coefficients ranged from -0.008 ± 0.002 for nervousness, flightiness, and gregariousness to -0.003 ± 0.002 for aggressiveness). In subsequent analyses, the regression of temperament trait on sequence of evaluation within a pen was highly significant and solutions ranged from 0.05 ± 0.007 for aggressiveness to 0.08 ± 0.007 for all other traits. The apparent large additive genetic variance for any one of these traits may be useful in identification of genes responsible for differences in cattle temperament.