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

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.

Physiological Stress Responses in Cattle Used in the Spanish Rodeo

Certain events can cause distress in cattle. In Spain, there is a sport similar to rodeo called persecution and takedown, in which calves are harassed and knocked down by riders. In this study, the physiological stress response of calves (n = 260) is assessed by measuring hormonal physiological parameters. Salivary samples were collected from Salers (n = 110) and Lidia (n = 150) calves before, during, and after the persecution and takedown event. The hormones epinephrine, cortisol, serotonin, and dopamine were determined in saliva samples using enzyme-immunoassay techniques. The results obtained revealed that epinephrine and cortisol levels increased during the event in Salers calves, with a significant increase (p < 0.05) in the case of epinephrine, although after the event, these values returned to their initial state. Therefore, this sport supposes an assumable punctual stressor stimulus for the animal. In contrast, in Lidia calves, cortisol and epinephrine levels decreased, with a significant decrease (p < 0.05) in the case of cortisol, which may be related to the temperament of this breed and facing a stressful situation in a different manner. This is confirmed by serotonin and dopamine levels that were altered in Lidia calves with respect to the other group studied. In conclusion, the sport of persecution and takedown produces a physiological response of adaptive stress assumable for the animals.

Behavioural Traits in Bos taurus Cattle, Their Heritability, Potential Genetic Markers, and Associations with Production Traits

Simple Summary

Cattle have the potential to seriously injure humans and cause damage to property. The risk of cattle reacting in a dangerous manner can be reduced through genetic selection for cattle which have a better temperament. A literature search was undertaken which returned papers which met the criteria of “Bovine”, “Genetics” and “Behaviour” or terms therein. Behavioural traits were grouped and their heritability, genomic associations and correlations with production traits examined. It was found that heritability estimates were more accurate in studies with large populations (n > 1000). Gene associations with behavioural traits were found on all chromosomes except for chromosome 13, with associated SNPs reported on all chromosomes except 5, 13, 17, 18 and 23. Generally, it was found that correlations between behaviour and production traits were low or negligible, suggesting that genetic improvement can be undertaken without negatively affecting production. There was variation between the results of the studies examined, and this underlines that any genetic study is population specific. Thus, to assess the heritability, genetic associations with production and genomic areas of interest for behavioural traits, a large-scale study of the population of interest would be required.

Abstract

People who work with cattle are at severe risk of serious injury due to the size and strength of the cattle. This risk can be minimised by breeding less dangerous cattle, which have a more favourable reaction to humans. This study provides a systematic review of literature pertaining to cattle genetics relating to behaviour. The review protocol was developed using the Preferred Reporting Items for Systematic reviews and Meta-Analyses (PRISMA) framework, with Population, Exposure and Outcome components identified as Bovine, Genetics and Behaviour respectively. Forty-nine studies were identified in the sifting and assigned non-exclusively to groups of heritability (22), genomic associations (13) and production traits related to behaviour (24). Behavioural traits were clustered into the following groups: “temperament, disposition and/ or docility”, “aggression”, “chute score”, “flight speed”, “milking temperament”, “non-restrained methods” and “restrained methods”. Fourteen papers reported high accuracy (Standard Error ≤ 0.05) estimates of heritability, the majority (n = 12) of these studies measured over 1000 animals. The heritability estimates were found to vary between studies. Gene associations with behavioural traits were found on all chromosomes except for chromosome 13, with associated SNPs reported on all chromosomes except 5, 13, 17, 18 and 23. Generally, it was found that correlations between behaviour and production traits were low or negligible. These studies suggest that additive improvement of behavioural traits in cattle is possible and would not negatively impact performance. However, the variation between studies demonstrates that the genetic relationships are population specific. Thus, to assess the heritability, genetic associations with production and genomic areas of interest for behavioural traits, a large-scale study of the population of interest would be required.

Exploration of Potential Ewing Sarcoma Drugs from FDA-Approved Pharmaceuticals through Computational Drug Repositioning, Pharmacogenomics, Molecular Docking, and MD Simulation Studies

Novel drug development is a time-consuming process with relatively high debilitating costs. To overcome this problem, computational drug repositioning approaches are being used to predict the possible therapeutic scaffolds against different diseases. In the current study, computational drug repositioning approaches were employed to fetch the promising drugs from the pool of FDA-approved drugs against Ewing sarcoma. The binding interaction patterns and conformational behaviors of screened drugs within the active region of Ewing sarcoma protein (EWS) were confirmed through molecular docking profiles. Furthermore, pharmacogenomics analysis was employed to check the possible associations of selected drugs with Ewing sarcoma genes. Moreover, the stability behavior of selected docked complexes (drugs-EWS) was checked by molecular dynamics simulations. Taken together, astemizole, sulfinpyrazone, and pranlukast exhibited a result comparable to pazopanib and can be used as a possible therapeutic agent in the treatment of Ewing sarcoma.

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.

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.

Association between temperament related traits and single nucleotide polymorphisms in the serotonin and oxytocin systems in Merino sheep

Animal temperament is defined as the consistent behavioral and physiological differences that are seen between individuals in response to the same stressor. Neurotransmitter systems, like serotonin and oxytocin in the central nervous system, underlie variation in behavioral traits in humans and other animals. Variations like single nucleotide polymorphisms (SNPs) in the genes for tryptophan 5-hydroxylase (TPH2), the serotonin transporter (SLC6A4), the serotonin receptor (HTR2A), and the oxytocin receptor (OXTR) are associated with behavioral phenotype in humans. Thus, the objective of this study was to identify SNPs in those genes and to test if those variations are associated with the temperament in Merino sheep. Using ewes from the University of Western Australia temperament flock, which has been selected on emotional reactivity for more than 20 generations, eight SNPs (rs107856757, rs107856818, rs107856856 and rs107857156 in TPH2, rs20917091 in SLC6A4, rs17196799 and rs17193181 in HTR2A, and rs17664565 in OXTR) were found to be distributed differently between calm and nervous sheep. These eight SNPs were then genotyped in 260 sheep from a flock that has never been selected on emotional reactivity, followed by the estimation of the behavioral traits of those 260 sheep using an arena test and an isolation box test. We found that several SNPs in TPH2 (rs107856757, rs107856818, rs107856856 and rs107857156) were in strong linkage disequilibrium, and all were associated with behavioral phenotype in the nonselected sheep. Similarly, rs17196799 in HTR2A was also associated with the behavioral phenotype.

Using imputed whole-genome sequence variants to uncover candidate mutations and genes affecting milking speed and temperament in Holstein cattle

Milking speed (MS) and temperament (MT) are 2 workability traits of great importance in dairy cattle production and breeding. This is mainly due to an increased intensification of the worldwide production systems and greater adoption of precision technologies with less human-cattle interaction. Both MS and MT are heritable traits and thus, genomic selection is a promising tool to expedite their genetic progress. However, the genetic architecture and biological mechanisms underlying the phenotypic expression of these traits remain underexplored. In this study, we investigated the association of >5.7 million imputed whole-genome sequence variants with MT and MS in 4,381 and 4,219 North American Holstein cattle, respectively. The statistical analyses were performed using a mixed linear model fitting a polygenic effect. We detected 40 and 35 significant SNPs independently associated with MT and MS, respectively, which were distributed across 26 chromosomes. Eight candidate genes (GRIN3A, KCNJ3, BOSTAUV1R417, BOSTAUV1R419, MAP2K5, KCTD3, GAP43, and LSAMP) were suggested to play an important role in MT as they are involved in biologically relevant pathways, such as glutamatergic synapse, vomeronasal receptor and oxytocin signaling. Within their coding and upstream sequences, we used an independent data set to further detect or validate significantly differentiated SNP between cattle breeds with known differences in MT. There were fewer candidate genes potentially implicated in MS, but immunity-related genes (e.g., BOLA-NC1 and LOC512672), also identified in other populations, were validated in this study. The significant SNP and novel candidate genes identified contribute to a better understanding of the biological mechanisms underlying both traits in dairy cattle. This information will also be useful for the optimization of prediction of genomic breeding values by giving greater weights to SNP located in the genomic regions identified.

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.

Influence of prenatal transportation stress-induced differential DNA methylation on the physiological control of behavior and stress response in suckling Brahman bull calves

The objective of this experiment was to examine potential differential methylation of DNA as a mechanism for altered behavioral and stress responses in prenatally stressed (PNS) compared with nonprenatally stressed (Control) young bull calves. Mature Brahman cows (n = 48) were transported for 2-h periods at 60 ± 5, 80 ± 5, 100 ± 5, 120 ± 5, and 140 ± 5 d of gestation (Transported group) or maintained as nontransported Controls (n = 48). From the offspring born to Transported and Control cows, a subset of 28-d-old intact bulls (n = 7 PNS; n = 7 Control) were evaluated for methylation of DNA of behavior and stress response-associated genes. Methylation of DNA from white blood cells was assessed via reduced representation bisulfite sequencing methods. Because increased methylation of DNA within gene promoter regions has been associated with decreased transcriptional activity of the corresponding gene, differentially methylated (P ≤ 0.05) CG sites (cytosine followed by a guanine nucleotide) located within promoter regions (n = 1,205) were used to predict (using Ingenuity Pathway Analysis software) alterations to canonical pathways in PNS compared with Control bull calves. Among differentially methylated genes (P ≤ 0.05) related to behavior and the stress response were OPRK1, OPRM1, PENK, POMC, NR3C2, TH, DRD1, DRD5, COMT, HTR6, HTR5A, GABRA4, GABRQ, and GAD2. Among altered (P < 0.05) signaling pathways related to behavior and the stress response were Opioid Signaling, Corticotropin-Releasing Hormone Signaling, Dopamine Receptor Signaling, Dopamine-DARPP32 Feedback in cAMP Signaling, Serotonin Receptor Signaling, and GABA Receptor Signaling. Alterations to behavior and stress response-related genes and canonical pathways supported previously observed elevations in temperament score and serum cortisol through weaning in the larger population of PNS calves from which bulls in this study were derived. Differential methylation of DNA and predicted alterations to behavior and stress response-related pathways in PNS compared with Control bull calves suggest epigenetic programming of behavior and the stress response in utero.

Aggressive behavior in cattle is associated with a polymorphism in the MAOA gene promoter

Molecular mechanisms underlying aggressive behavior are primitive and similar among the subphylum Vertebrata. In humans, a primary goal in the study of aggression is to determine the neurobehavioral molecular factors triggering violence. Although several species have been used to study agonistic responses, researchers are limited by the difficulty of artificially inducing aggression in animals not selected for it. Conversely, the Lidia cattle breed has been selected since the eighteenth century to display agonistic responses based on traits such as aggressiveness, ferocity and mobility, all of them showing significant heritability values. This intensive selection may have driven shifts in specific allele frequencies. In a previous analysis across the autosomes, we revealed long-term selection regions including genes involved in behavioral development. In the present study, we focus on mapping recent signatures of selection associated with aggressiveness at chromosome X, by comparing Lidia cattle samples with two non-specialized Spanish breeds showing tamed behavior. The most significant markers peaked around the monoamine oxidase A (MAOA) gene, and thus the associations of three functionally important regions located near the promoter of this gene were further investigated. A polymorphism consisting of a variable number of tandem repeats of the nucleotide 'C' (BTX:105,462,494) and displaying lower number of repetitions in the Lidia breed when compared with the tamed breeds was detected. In silico analyses predicted that the g.105,462,494delsinsC variant may code for the Sp1 binding motif, one of the major transcription factors controlling the core promoter and expression of the MAOA gene in humans.

Demography and rapid local adaptation shape Creole cattle genome diversity in the tropics

The introduction of Iberian cattle in the Americas after Columbus' arrival imposed high selection pressures on a limited number of animals over a brief period of time. Knowledge of the genomic regions selected during this process may help in enhancing climatic resilience and sustainable animal production. We first determined taurine and indicine contributions to the genomic structure of modern Creole cattle. Second, we inferred their demographic history using approximate Bayesian computation (ABC), linkage disequilibrium (LD) and N e Slope (NeS) analysis. Third, we performed whole genome scans for selection signatures based on cross-population extended haplotype homozygosity (XP-EHH) and population differentiation (F ST) to disentangle the genetic mechanisms involved in adaptation and phenotypic change by a rapid and major environmental transition. To tackle these questions, we combined SNP array data (~54,000 SNPs) in Creole breeds with their modern putative Iberian ancestors. Reconstruction of the population history of Creoles from the end of the 15th century indicated a major demographic expansion until the introduction of zebu and commercial breeds into the Americas ~180 years ago, coinciding with a drastic N e contraction. NeS analysis provided insights into short-term complexity in population change and depicted a decrease/expansion episode at the end of the ABC-inferred expansion, as well as several additional fluctuations in N e with the attainment of the current small N e only towards the end of the 20th century. Selection signatures for tropical adaptation pinpointed the thermoregulatory slick hair coat region, identifying a new candidate gene (GDNF), as well as novel candidate regions involved in immune function, behavioural processes, iron metabolism and adaptation to new feeding conditions. The outcomes from this study will help in future-proofing farm animal genetic resources (FAnGR) by providing molecular tools that allow selection for improved cattle performance, resilience and welfare under climate change.

Prenatal transportation stress alters genome-wide DNA methylation in suckling Brahman bull calves

The objective of this experiment was to identify genome-wide differential methylation of DNA in young prenatally stressed (PNS) bull calves. Mature Brahman cows (n = 48) were transported for 2-h periods at 60 ± 5, 80 ± 5, 100 ± 5, 120 ± 5, and 140 ± 5 d of gestation or maintained as nontransported Controls (n = 48). Methylation of DNA from white blood cells from a subset of 28-d-old intact male offspring (n = 7 PNS; n = 7 Control) was assessed via reduced representation bisulfite sequencing. Samples from PNS bulls contained 16,128 CG, 226 CHG, and 391 CHH (C = cytosine; G = guanine; H = either adenine, thymine, or cytosine) sites that were differentially methylated compared to samples from Controls. Of the CG sites, 7,407 were hypermethylated (at least 10% more methylated than Controls; P ≤ 0.05) and 8,721 were hypomethylated (at least 10% less methylated than Controls; P ≤ 0.05). Increased DNA methylation in gene promoter regions typically results in decreased transcriptional activity of the region. Therefore, differentially methylated CG sites located within promoter regions (n = 1,205) were used to predict (using Ingenuity Pathway Analysis software) alterations to canonical pathways in PNS compared with Control bull calves. In PNS bull calves, 113 pathways were altered (P ≤ 0.05) compared to Controls. Among these were pathways related to behavior, stress response, metabolism, immune function, and cell signaling. Genome-wide differential DNA methylation and predicted alterations to pathways in PNS compared with Control bull calves suggest epigenetic programming of biological systems in utero.

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.