Identification of Candidate Genes for Reactivity in Guzerat (Bos indicus) Cattle: A Genome-Wide Association Study

Behavioural and physiological responses to stressors in sheep with temperament classified by genotype or phenotype

The single nucleotide polymorphism (SNP) rs107856856, located in the tryptophan hydroxylase-2 gene, is associated with the behavioural phenotype for sheep temperament measured at weaning. Here, we tested the association between that SNP and physiological and behavioural responses to stressors in adult sheep. Two groups of adult sheep, one with genotype A/A (calm genotype) and the other with G/G (nervous genotype) in rs107856856, were selected from 160 sheep and were exposed, twice, to an open-field arena and an isolation box test (IBT). During each repeat, the behaviour and physiological responses (cortisol, prolactin, dehydroepiandrosterone [DHEA], brain derived neurotrophic factor [BDNF], characteristics of the response of body temperature, and oxidative stress) were measured. The behavioural and physiological responses of the sheep were compared between genotypes and also between groups classified on their phenotype as assessed by their initial isolation box score ("low responders" and "high responders"). The SNP rs107856856 had some effects on the behavioural phenotype (IBT score) but no effects on the physiological response to stress (cortisol, prolactin, DHEA, BDNF, oxidative stress or changes in body temperature) in the adult sheep, probably because the sheep were exposed, and therefore had adapted, to human contact during their life.

Assessment of runs of homozygosity, heterozygosity-rich regions and genomic inbreeding estimates in a subpopulation of Guzerá (Bos indicus) dual-purpose cattle

For decades, inbreeding in cattle has been evaluated using pedigree information. Nowadays, inbreeding coefficients can be obtained using genomic information such as runs of homozygosity (ROH). The aims of this study were to quantify ROH and heterozygosity-rich regions (HRR) in a subpopulation of Guzerá dual-purpose cattle, to examine ROH and HRR islands, and to compare inbreeding coefficients obtained by ROH with alternative genomic inbreeding coefficients. A subpopulation of 1733 Guzerá animals genotyped for 50k SNPs was used to obtain the ROH and HRR segments. Inbreeding coefficients by ROH (FROH ), by genomic relationship matrix based on VanRaden's method 1 using reference allele frequency in the population (FGRM ), by genomic relationship matrix based on VanRaden's method 1 using allele frequency fixed in 0.5 (FGRM_0.5 ), and by the proportion of homozygous loci (FHOM ) were calculated. A total of 15,660 ROH were identified, and the chromosome with the highest number of ROH was BTA6. A total of 4843 HRRs were identified, and the chromosome with the highest number of HRRs was BTA23. No ROH and HRR islands were identified according to established criteria, but the regions closest to the definition of an island were examined from 64 to 67 Mb of BTA6, from 36 to 37 Mb of BTA2 and from 0.50 to 1.25 Mb of BTA23. The genes identified in ROH islands have previously been associated with dairy and beef traits, while genes identified on HRR islands have previously been associated with reproductive traits and disease resistance. FROH was equal to 0.095 ± 0.084, and its Spearman correlation with FGRM was low (0.44) and moderate-high with FHOM (0.79) and with FGRM_0.5 (0.80). The inbreeding coefficients determined by ROH were higher than other cattle breeds' and higher than pedigree-based inbreeding in the Guzerá breed obtained in previous studies. It is recommended that future studies investigate the effects of inbreeding determined by ROH on the traits under selection in the subpopulation studied.

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.

Genomic Selection for Dairy Cattle Behaviour Considering Novel Traits in a Changing Technical Production Environment

Cow behaviour is a major factor influencing dairy herd profitability and is an indicator of animal welfare and disease. Behaviour is a complex network of behavioural patterns in response to environmental and social stimuli and human handling. Advances in agricultural technology have led to changes in dairy cow husbandry systems worldwide. Increasing herd sizes, less time availability to take care of the animals and modern technology such as automatic milking systems (AMSs) imply limited human-cow interactions. On the other hand, cow behaviour responses to the technical environment (cow-AMS interactions) simultaneously improve production efficiency and welfare and contribute to simplified "cow handling" and reduced labour time. Automatic milking systems generate objective behaviour traits linked to workability, milkability and health, which can be implemented into genomic selection tools. However, there is insufficient understanding of the genetic mechanisms influencing cow learning and social behaviour, in turn affecting herd management, productivity and welfare. Moreover, physiological and molecular biomarkers such as heart rate, neurotransmitters and hormones might be useful indicators and predictors of cow behaviour. This review gives an overview of published behaviour studies in dairy cows in the context of genetics and genomics and discusses possibilities for breeding approaches to achieve desired behaviour in a technical production environment.

Genomic prediction in Brazilian Guzerá cattle: application of a single-step approach to productive and reproductive traits

This study aimed to investigate the feasibility of genomic prediction for productive and reproductive traits in Guzerá cattle using single-step genomic best linear unbiased prediction (ssGBLUP). Evaluations included the 305-day cumulative yields (first lactation, in kg) of milk, lactose, protein, fat, and total solids; adjusted body weight (kg) at the ages of 450, 365, and 210 days; and age at first calving (in days), from a database containing 197,283 measurements from Guzerá males and females born between 1954 and 2018. The pedigree included 433,823 animals spanning up to 14 overlapping generations. A total of 1618 animals were genotyped. The analyses were performed using ssGBLUP and traditional BLUP methods. Predictive ability and bias were accessed using cross-validation: predictive ability was similar between the methods and ranged from 0.27 to 0.47 for the genomic-based model and from 0.30 to 0.45 for the pedigree-based model; the bias was also similar between the methods, ranging from 0.88 to 1.35 in the genomic-based model and from 0.96 to 1.41 in the pedigree-based model. The individual accuracies of breeding values were evidently increased in the genomic evaluation, with values ranging from 0.41 to 0.56 in the genomic-based model and from 0.26 to 0.54 in the pedigree-based model. Even based on a small number of genotyped animals and a small database for some traits, the results suggest that ssGBLUP is feasible and may be applied to national genetic evaluation of the breed to increase the accuracy of breeding values without greatly impacting predictive ability and bias.

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.

Cross-population selection signatures in Canchim composite beef cattle

Analyses of livestock genomes have been used to detect selection signatures, which are genomic regions associated with traits under selection leading to a change in allele frequency. The objective of the present study was to characterize selection signatures in Canchim composite beef cattle using cross-population analyses with the founder Nelore and Charolais breeds. High-density single nucleotide polymorphism genotypes were available on 395 Canchim representing the target population, along with genotypes from 809 Nelore and 897 Charolais animals representing the reference populations. Most of the selection signatures were co-located with genes whose functions agree with the expectations of the breeding programs; these genes have previously been reported to associate with meat quality, as well as reproductive traits. Identified genes were related to immunity, adaptation, morphology, as well as behavior, could give new perspectives for understanding the genetic architecture of Canchim. Some selection signatures identified genes that were recently introduced in Canchim, such as the loci related to the polled trait.

Identification of important genomic footprints using eight different selection signature statistics in domestic cattle breeds

In the present study, the population genomic data of different cattle breeds were explored to decipher the genomic regions affected due to selective events and reflected in the productive, reproductive, thermo-tolerance, and health-related traits. To find out these genomic deviations due to selective sweeps, we used eight different statistical tools (Tajima's D, Fu & Li's D*, CLR, ROH, iHS, F_ST, FLK, and hapFLK) on seven indigenous and five exotic cattle breeds. We further performed composite analysis by comparing their covariance matrix. Several candidate genes were found to be related to milk production (ADARB, WDR70, and CA8), reproductive (PARN, FAM134B2, and ZBTB20), and health-related traits (SP110, CXCL2, CLXCL3, CXCL5, IRF8, and MYOM1). The outcome of this investigation provides a basis for detecting selective sweeps that explain the genetic variation of traits. They may possess functional importance for multiple cattle breeds in different subcontinents. However, further studies are required to improve the findings using high-density arrays or whole-genome sequencing with higher resolution and greater sample sizes.

Genomic Signatures of Domestication Selection in the Australasian Snapper (Chrysophrys auratus)

Domestication of teleost fish is a recent development, and in most cases started less than 50 years ago. Shedding light on the genomic changes in key economic traits during the domestication process can provide crucial insights into the evolutionary processes involved and help inform selective breeding programmes. Here we report on the recent domestication of a native marine teleost species in New Zealand, the Australasian snapper (Chrysophrys auratus). Specifically, we use genome-wide data from a three-generation pedigree of this species to uncover genetic signatures of domestication selection for growth. Genotyping-By-Sequencing (GBS) was used to generate genome-wide SNP data from a three-generation pedigree to calculate generation-wide averages of FST between every generation pair. The level of differentiation between generations was further investigated using ADMIXTURE analysis and Principal Component Analysis (PCA). After that, genome scans using Bayescan, LFMM and XP-EHH were applied to identify SNP variants under putative selection following selection for growth. Finally, genes near candidate SNP variants were annotated to gain functional insights. Analysis showed that between generations FST values slightly increased as generational time increased. The extent of these changes was small, and both ADMIXTURE analysis and PCA were unable to form clear clusters. Genome scans revealed a number of SNP outliers, indicative of selection, of which a small number overlapped across analyses methods and populations. Genes of interest within proximity of putative selective SNPs were related to biological functions, and revealed an association with growth, immunity, neural development and behaviour, and tumour repression. Even though few genes overlapped between outlier SNP methods, gene functionalities showed greater overlap between methods. While the genetic changes observed were small in most cases, a number of outlier SNPs could be identified, of which some were found by more than one method. Multiple outlier SNPs appeared to be predominately linked to gene functionalities that modulate growth and survival. Ultimately, the results help to shed light on the genomic changes occurring during the early stages of domestication selection in teleost fish species such as snapper, and will provide useful candidates for the ongoing selective breeding in the future of this and related species.

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.

A genome-wide association study for prolificacy in three Polish sheep breeds

Reproductive traits (especially litter size) are usually characterized by low heritability, and thus, phenotypic selection is often ineffective and slow. In order to improve fertility characteristics such as ovulation rate and litter size, it seems more effective to select breeding animals based on their genotype. The aim of the study was to use genome-wide association study (GWAS) in three sheep breeds to identify the genetic variants affecting the litter size in sheep. The study allowed us to identify one genome-wide significant SNP (rs402032081—located in ephrin type-A receptor 6, EPHA6) showing an association with litter size in Polish Mountain Sheep. We suggest that the EPHA6 gene can be a candidate gene for prolificacy trait in selected breeds of sheep; however, it needs further functional data for validation.

Identification of QTL regions and candidate genes for growth and feed efficiency in broilers

BACKGROUND

Feed accounts for about 70% of the total cost of poultry meat production. Residual feed intake (RFI) has become the preferred measure of feed efficiency because it is phenotypically independent of growth rate and body weight. In this study, our aim was to estimate genetic parameters and identify quantitative trait loci (QTL) for feed efficiency in 3314 purebred broilers using a genome-wide association study. Broilers were genotyped using a custom 55 K single nucleotide polymorphism (SNP) array.

RESULTS

Estimates of genomic heritability for seven growth and feed efficiency traits, including body weight at 28 days of age (BW28), BW42, average daily feed intake (ADFI), RFI, and RFI adjusted for weight of abdominal fat (RFIa), ranged from 0.12 to 0.26. Eleven genome-wide significant SNPs and 15 suggestively significant SNPs were detected, of which 19 clustered around two genomic regions. A region on chromosome 16 (2.34-2.66 Mb) was associated with both BW28 and BW42, and the most significant SNP in this region, AX_101003762, accounted for 7.6% of the genetic variance of BW28. The other region, on chromosome 1 (91.27-92.43 Mb) was associated with RFI and ADFI, and contains the NSUN3 and EPHA6 as candidate genes. The most significant SNP in this region, AX_172588157, accounted for 4.4% of the genetic variance of RFI. In addition, a genomic region containing the gene AGK on chromosome 1 was found to be associated with RFIa. The NSUN3 and AGK genes were found to be differentially expressed in breast muscle, thigh muscle, and abdominal fat between male broilers with high and low RFI.

CONCLUSIONS

We identified QTL regions for BW28 and BW42 (spanning 0.32 Mb) and RFI (spanning 1.16 Mb). The NSUN3, EPHA6, and AGK were identified as the most likely candidate genes for these QTL. These genes are involved in mitochondrial function and behavioral regulation. These results contribute to the identification of candidate genes and variants for growth and feed efficiency in poultry.

Genome-wide detection of signatures of selection in indicine and Brazilian locally adapted taurine cattle breeds using whole-genome re-sequencing data

BACKGROUND

The cattle introduced by European conquerors during the Brazilian colonization period were exposed to a process of natural selection in different types of biomes throughout the country, leading to the development of locally adapted cattle breeds. In this study, whole-genome re-sequencing data from indicine and Brazilian locally adapted taurine cattle breeds were used to detect genomic regions under selective pressure. Within-population and cross-population statistics were combined separately in a single score using the de-correlated composite of multiple signals (DCMS) method. Putative sweep regions were revealed by assessing the top 1% of the empirical distribution generated by the DCMS statistics.

RESULTS

A total of 33,328,447 biallelic SNPs with an average read depth of 12.4X passed the hard filtering process and were used to access putative sweep regions. Admixture has occurred in some locally adapted taurine populations due to the introgression of exotic breeds. The genomic inbreeding coefficient based on runs of homozygosity (ROH) concurred with the populations' historical background. Signatures of selection retrieved from the DCMS statistics provided a comprehensive set of putative candidate genes and revealed QTLs disclosing cattle production traits and adaptation to the challenging environments. Additionally, several candidate regions overlapped with previous regions under selection described in the literature for other cattle breeds.

CONCLUSION

The current study reported putative sweep regions that can provide important insights to better understand the selective forces shaping the genome of the indicine and Brazilian locally adapted taurine cattle breeds. Such regions likely harbor traces of natural selection pressures by which these populations have been exposed and may elucidate footprints for adaptation to the challenging climatic conditions.

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 Identifies Genomic Loci Associated With Neurotransmitter Concentration in Cattle

Abnormal neurotransmitter concentration is one of the factors that affect the health status, behavioral personality, and welfare level of animals, but the genetic basis of the abnormality is still largely unknown. The objective of this study is to identify putative genomic loci associated with neurotransmitter concentration in cattle. We measured serotonin (5HT), dopamine (DA), cortisol, glutamate (Glu), and ACTH concentrations in blood serum using double-antibody sandwich ELISA in 30 Brahman cattle and 127 Yunling cattle. Interestingly, we found that ACTH concentration was positively correlated with body weight, cannon circumference, and hip width (P < 0.05). Genome-wide association study (GWAS) was performed with mixed linear models using autosomal SNPs derived from the whole-genome sequence. We identified five, five, two, three, and five suggestive loci associated with 5HT, DA, cortisol, Glu, and ACTH concentration, respectively. These 20 associated loci implicated 18 candidate genes. For Glu concentration, the most significant association locus was assigned to MCHR1, a G-coupled receptor that could modulate glutamate release. For dopamine concentration, a very strong association locus was located in the intron of SLC18A2, which is a critical mediator of dopamine dynamics. However, for ACTH concentration, a very strong association locus was assigned to HTR1F, a G protein-coupled receptor that can influence the release of ACTH. Other candidate genes of interest identified for neurotransmitter concentration were PRMT6, GADD45A, PCCA, ANGPT1, ACCS, LOC100336971, TNR, GSDMA, CNTN3, CARMIL1, CDKAL1, RBFOX1, PCDH15, and LGALS12. Our findings will provide targets for the genetic improvement of neurotransmitter-related traits in domestic cattle and basic materials for studying the mechanism of neurotransmitter synthesis, release, and transport in human and animals.

Genome-wide association study identifies genomic loci associated with flight reaction in cattle

Lower flight reaction is closely related to higher production in cattle, but the genetic basis of lower flight reaction is not clearly understood. Here, we sampled a total of 45 Brahman cattles and 166 Yunling cattles with flight distance (FD), and 73 Brahman cattles and 288 Yunling cattles with crush score (CS) and flight speed (FS), whereas there were 45 Brahman cattles and 161 Yunling cattles with all three traits. The FD, CS and FS in Brahman cattle were significantly lower than those in Yunling cattle. The flight reaction traits had negative correlation with conformational traits (e.g., body weight, withers height and body length). Based on SNPs derived from a subset of 162 whole genomes (25 Brahman genomes and 100 Yunling genomes with FD, 30 Brahman and 131 Yunling genomes with CS and FS), genome-wide association study with mixed linear model was performed to test potential associations between flight reaction traits and genomic variants. We identified five, two and two genomic loci suggestively associated with FD, CS and FS, respectively. Five out of five candidate genes for FD (LOC789753, LRP6, CTIF, SLC9A9 and ZEB1) were reported to be related to Alzheimer's disease representing cognitive impairment in human, which was consistent with the finding that cognitive-behavioural intervention decreased the FD of cows to human. In CS, a very strong association locus was assigned to CDH8, a cadherin involved in synaptic adhesion, axon outgrowth and guidance, whose deletion was associated with autism spectrum disorder. In FS, a very strong association locus was assigned to GABRG2, a gamma-aminobutyric acid (a major inhibitory neurotransmitter in brain) receptor, whose polymorphisms were associated with suicidal behaviour in schizophrenia patients. Our findings will provide targets for molecular-marker selection and genetic manipulation of cattle improvement to meeting the growing demand for lower flight reaction to human.

Genome-wide detection of copy-number variations in local cattle breeds

The aim of the present study was to identify copy-number variations (CNVs) in Cinisara (CIN) and Modicana (MOD) cattle breeds on the basis of signal intensity (logR ratio) and B allele frequency of each marker, using Illumina’s BovineSNP50K Genotyping BeadChip. The CNVs were detected with the PennCNV and SVS 8.7.0 software and were aggregated into CNV regions (CNVRs). PennCNV identified 487 CNVs in CIN that aggregated into 86 CNVRs, and 424 CNVs in MOD that aggregated into 81 CNVRs. SVS identified a total of 207 CNVs in CIN that aggregated into 39 CNVRs, and 181 CNVs in MOD that aggregated into 41 CNVRs. The CNVRs identified with the two softwares contained 29 common CNVRs in CIN and 17 common CNVRs in MOD. Only a small number of CNVRs identified in the present study have been identified elsewhere, probably because of the limitations of the array used. In total, 178 and 208 genes were found within the CNVRs of CIN and MOD respectively. Gene Ontology and KEGG pathway analyses showed that several of these genes are involved in milk production, reproduction and behaviour, the immune response, and resistance/susceptibility to infectious diseases. Our results have provided significant information for the construction of more-complete CNV maps of the bovine genome and offer an important resource for the investigation of genomic changes and traits of interest in the CIN and MOD cattle breeds. Our results will also be valuable for future studies and constitute a preliminary report of the CNV distribution resources in local cattle genomes.

Autozygosity islands and ROH patterns in Nellore lineages: evidence of selection for functionally important traits

Background

The aim of this study was to assess genome-wide autozygosity in a Nellore cattle population and to characterize ROH patterns and autozygosity islands that may have occurred due to selection within its lineages. It attempts also to compare estimates of inbreeding calculated from ROH (F_ROH), genomic relationship matrix (F_GRM), and pedigree-based coefficient (F_PED).

Results

The average number of ROH per animal was 55.15 ± 13.01 with an average size of 3.24 Mb. The Nellore genome is composed mostly by a high number of shorter segments accounting for 78% of all ROH, although the proportion of the genome covered by them was relatively small. The genome autozygosity proportion indicates moderate to high inbreeding levels for classical standards, with an average value of 7.15% (178.70 Mb). The average of F_PED and F_ROH, and their correlations (− 0.05 to 0.26) were low. Estimates of correlation between F_GRM-F_PED was zero, while the correlation (− 0.01 to − 0.07) between F_GRM-F_ROH decreased as a function of ROH length, except for F_{ROH > 8Mb} (− 0.03). Overall, inbreeding coefficients were not high for the genotyped animals. Autozygosity islands were evident across the genome (n = 62) and their genomic location did not largely differ within lineages. Enriched terms (p < 0.01) associated with defense response to bacteria (GO:0042742), immune complex reaction (GO:0045647), pregnancy-associated glycoproteins genes (GO:0030163), and organism growth (GO:0040014) were described within the autozygotic islands.

Conclusions

Low F_PED-F_ROH correlation estimates indicate that F_PED is not the most suitable method for capturing ancient inbreeding when the pedigree does not extend back many generations and F_ROH should be used instead. Enriched terms (p < 0.01) suggest a strong selection for immune response. Non-overlapping islands within the lineages greatly explain the mechanism underlying selection for functionally important traits in Nellore cattle.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5060-8) contains supplementary material, which is available to authorized users.