A systematic review with meta-analysis of heritability estimates for temperament-related traits in beef and dairy cattle populations

Temperament (docility) is a key breeding goal in the cattle industry due to its direct relationship with animal welfare, cattle handler's safety and animal productivity. Over the past six decades, numerous studies have reported heritability estimates for temperament-related traits in cattle populations ranging from low to high values. Therefore, the primary objective of this study was to perform a comprehensive systematic review with meta-analysis to obtain weighted estimates of heritability for temperament-related traits in worldwide cattle populations. After data editing and quality control, 106 studies were included in the systematic review, of which 29.2% and 70.8% reported estimates of heritability for temperament-related traits in dairy and beef cattle populations, respectively. Meta-analyses were performed for 95 heritability estimates using a random model approach. The weighted heritability estimates were as follow: (a) flight score at weaning = 0.23 (95% CI: 0.15-0.32); (b) flight speed at weaning = 0.30 (95% CI: 0.26-0.33); (c) joint analysis of flight speed and flight score at weaning = 0.27 (95% CI: 0.22-0.31); (d) flight speed at yearling = 0.26 (95% CI: 0.21-0.30); (e) joint analysis of flight speed at weaning and yearling = 0.27 (95% CI: 0.24-0.30); (f) movement score = 0.12 (95% CI: 0.08-0.15); (g) crush score at weaning = 0.21 (95% CI: 0.17-0.25); (h) pen score at weaning = 0.27 (95% CI: 0.19-0.34); (i) pen score at yearling = 0.20 (95% CI: 0.17-0.23); (j) joint analysis of pen score at weaning and yearling = 0.22 (95% CI: 0.18-0.26); (k) cow's aggressiveness at calving = 0.10 (95% CI: 0.01-0.19); (l) general temperament = 0.13 (95% CI: 0.06-0.19); (m) milking temperament = 0.16 (95% CI: 0.11-0.21); and (n) joint analysis of general and milking temperament = 0.14 (95% CI: 0.11-0.18). The heterogeneity index ranged from 0% to 77%, and the Q-test was significant (p < 0.05) for four single-trait meta-analyses. In conclusion, temperament is moderately heritable in beef cattle populations, and flight speed at weaning had the highest weighted heritability estimate. Moreover, between-study heterogeneity was low or moderate in beef cattle traits, suggesting reasonable standardization across studies. On the other hand, low-weighted heritability and high between-study heterogeneity were estimated for temperament-related traits in dairy cattle, suggesting that more studies are needed to better understand the genetic inheritance of temperament in dairy cattle populations.

Genetic and phenotypic parameters for sexual precocity and parasite resistance traits in Nellore cattle

Indicator traits of sexual precocity are widely used as selection criteria for the genetic improvement of beef cattle; however, the impact of selection for these traits on resistance to endoparasites and ectoparasites is unknown. Therefore, this study aimed to estimate the genetic and phenotypic parameters for indicator traits of sexual precocity and parasite resistance in Nellore cattle. The sexual precocity traits evaluated were probability of first calving (PFC) and scrotal circumference at 12 and 18 months of age (SC12 and SC18). The resistance-related traits included tick (TC), gastrointestinal nematode egg (NEC), and Eimeria spp. oocyst (EOC) counts. (Co)variance components were estimated by Bayesian inference using multitrait animal models. The mean heritabilities for PFC, SC12, SC18, TC, NEC, and EOC were 0.23, 0.38, 0.42, 0.14, 0.16, and 0.06, respectively, and suggest that selection will change the mean values of these traits over time. The genetic and phenotypic correlations for most pairs formed by a precocity and a resistance trait were not different from zero, suggesting that selection for sexual precocity traits will not result in changes in resistance traits. Thus, selection for indicator traits of sexual precocity does not elicit unfavorable correlated responses in resistance to endoparasites and ectoparasites, and joint selection aimed at improving these traits can be performed using multitrait selection methods, when necessary.

Genetic parameters for behavioral and growth traits of Nellore cattle

The growing concern of consumers with the welfare of production animals searches welfare in a production system extremely important; thus, the study of animal temperament is necessary to select less excitable temperament animals resulting in healthy development and fewer accidents. The objective of this study was to estimate genetic parameters for traits related to animal temperament and growth traits of Nellore cattle. In addition to exploring the genetic pattern of these traits through cluster and principal component analysis (PCA), to reveal possible groups of individuals that express less excitable temperament and greater growth. A total of 2,332 measurements from 1,245 male and female Nellore cattle born between 2008 and 2016 were utilized in the study. The (co)variance components were estimated by Bayesian inference using a two-trait animal model. The heritability for temperament score (TS), flight speed (FS), body condition score (BCS), live weight (LW), and hip height (HH) were 0.08, 0.12, 0.06, 0.13, and 0.48, respectively. The genetic correlation between the temperament indicator traits was strong and positive (0.78 ± 0.24). The TS and FS showed a favorable or null genetic correlation with LW, BCS, and HH. The third cluster included animals with low EBV for TS and FS and with high EBV for BCS, LW, and HH. In the PCA, the PC1 was what best evidenced the aim of this study; thus, our findings suggest that we could explore select animals based on cluster 3 and PC1 in breeding programs to select Nellore cattle with less excitable temperament and greater growth.

Prevalence and tick loads in Nguni cattle reared in different environmental conditions across four provinces of South Africa

Background and Aim:

In tropical and subtropical countries, ixodid ticks are among livestock’s most economically important ectoparasites. Although Nguni cattle from South Africa have adapted to harsh environments, it is unknown whether they will be resistant to ticks, and the diseases carried by ticks under various climatic conditions. Therefore, this study aimed to compare tick load and estimate the prevalence of different tick species among Nguni cattle under different environmental conditions.

Materials and Methods:

Tick counts were conducted monthly under natural challenges over 2 years on 586 Nguni cattle located at ARC-Roodeplaat and Loskop farms (warmer climate), Mukhuthali Nguni Community and the University of Fort Hare farms (cooler climate). The generalized linear model procedure of the Statistical Analysis System was used to analyze the data. It fitted the location (farm), sex, year, month or season, and animal age as covariates.

Results:

The tick species (relative prevalence) observed were as follows: Amblyomma hebraeum (42%), Rhipicephalus evertsi (22%), Rhipicephalus (Boophilus) spp. (16%), Rhipicephalus appendiculatus (11%), Hyalomma marginatum (5%), and Rhipicephalus simus (4%). Tick infestation was significantly affected by location, season, year, month of the tick counting and age of the animal. Loskop farm had the highest tick count (m = 30.69) and showed the largest variation in tick count. Compared to the other seasons, higher tick counts were seen during the hot-dry (September–November) and hot-wet (December–February) seasons. A. hebraeum was the dominant tick species across all four farms, followed by R. evertsi. The perianal region (under the tail head), the perineum and the belly body locations were the most preferred tick attachment sites.

Conclusion:

These results provide useful information for developing appropriate control strategies for ticks and tick-borne diseases in these provinces of South Africa. Further work must investigate the feasibility of genetic improvement for tick resistance.

Genome variation in tick infestation and cryptic divergence in Tunisian indigenous sheep

Background

Ticks are obligate haematophagous ectoparasites considered second to mosquitos as vectors and reservoirs of multiple pathogens of global concern. Individual variation in tick infestation has been reported in indigenous sheep, but its genetic control remains unknown.

Results

Here, we report 397 genome-wide signatures of selection overlapping 991 genes from the analysis, using ROH, LR-GWAS, XP-EHH, and F_ST, of 600 K SNP genotype data from 165 Tunisian sheep showing high and low levels of tick infestations and piroplasm infections. We consider 45 signatures that are detected by consensus results of at least two methods as high-confidence selection regions. These spanned 104 genes which included immune system function genes, solute carriers and chemokine receptor. One region spanned STX5, that has been associated with tick resistance in cattle, implicating it as a prime candidate in sheep. We also observed RAB6B and TF in a high confidence candidate region that has been associated with growth traits suggesting natural selection is enhancing growth and developmental stability under tick challenge. The analysis also revealed fine-scale genome structure indicative of cryptic divergence in Tunisian sheep.

Conclusions

Our findings provide a genomic reference that can enhance the understanding of the genetic architecture of tick resistance and cryptic divergence in indigenous African sheep.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08321-1.

The Genetics of Thermoregulation in Pigs: A Review

Heat stress (HS) affects pig performance, health and welfare, resulting in a financial burden to the pig industry. Pigs have a limited number of functional sweat glands and their thermoregulatory mechanisms used to maintain body temperature, are challenged by HS to maintain body temperature. The genetic selection of genotypes tolerant to HS is a promising long-term (adaptation) option that could be combined with other measures at the production system level. This review summarizes the current knowledge on the genetics of thermoregulation in pigs. It also discusses the different phenotypes that can be used in genetic studies, as well as the variability in thermoregulation between pig breeds and the inheritance of traits related to thermoregulation. This review also considers on-going challenges to face for improving heat tolerance in pigs.

Suitability of GWAS as a Tool to Discover SNPs Associated with Tick Resistance in Cattle: A Review

Understanding the biological mechanisms underlying tick resistance in cattle holds the potential to facilitate genetic improvement through selective breeding. Genome wide association studies (GWAS) are popular in research on unraveling genetic determinants underlying complex traits such as tick resistance. To date, various studies have been published on single nucleotide polymorphisms (SNPs) associated with tick resistance in cattle. The discovery of SNPs related to tick resistance has led to the mapping of associated candidate genes. Despite the success of these studies, information on genetic determinants associated with tick resistance in cattle is still limited. This warrants the need for more studies to be conducted. In Africa, the cost of genotyping is still relatively expensive; thus, conducting GWAS is a challenge, as the minimum number of animals recommended cannot be genotyped. These population size and genotype cost challenges may be overcome through the establishment of collaborations. Thus, the current review discusses GWAS as a tool to uncover SNPs associated with tick resistance, by focusing on the study design, association analysis, factors influencing the success of GWAS, and the progress on cattle tick resistance studies.

Multiple Country and Breed Genomic Prediction of Tick Resistance in Beef Cattle

Ticks cause substantial production losses for beef and dairy cattle. Cattle resistance to ticks is one of the most important factors affecting tick control, but largely neglected due to the challenge of phenotyping. In this study, we evaluate the pooling of tick resistance phenotyped reference populations from multi-country beef cattle breeds to assess the possibility of improving host resistance through multi-trait genomic selection. Data consisted of tick counts or scores assessing the number of female ticks at least 4.5 mm length and derived from seven populations, with breed, country, number of records and genotyped/phenotyped animals being respectively: Angus (AN), Brazil, 2,263, 921/1,156, Hereford (HH), Brazil, 6,615, 1,910/2,802, Brangus (BN), Brazil, 2,441, 851/851, Braford (BO), Brazil, 9,523, 3,062/4,095, Tropical Composite (TC), Australia, 229, 229/229, Brahman (BR), Australia, 675, 675/675, and Nguni (NG), South Africa, 490, 490/490. All populations were genotyped using medium density Illumina SNP BeadChips and imputed to a common high-density panel of 332,468 markers. The mean linkage disequilibrium (LD) between adjacent SNPs varied from 0.24 to 0.37 across populations and so was sufficient to allow genomic breeding values (GEBV) prediction. Correlations of LD phase between breeds were higher between composites and their founder breeds (0.81 to 0.95) and lower between NG and the other breeds (0.27 and 0.35). There was wide range of estimated heritability (0.05 and 0.42) and genetic correlation (-0.01 and 0.87) for tick resistance across the studied populations, with the largest genetic correlation observed between BN and BO. Predictive ability was improved under the old-young validation for three of the seven populations using a multi-trait approach compared to a single trait within-population prediction, while whole and partial data GEBV correlations increased in all cases, with relative improvements ranging from 3% for BO to 64% for TC. Moreover, the multi-trait analysis was useful to correct typical over-dispersion of the GEBV. Results from this study indicate that a joint genomic evaluation of AN, HH, BN, BO and BR can be readily implemented to improve tick resistance of these populations using selection on GEBV. For NG and TC additional phenotyping will be required to obtain accurate GEBV.

Tropical Beef: Is There an Axiomatic Basis to Define the Concept?

Cattle production in tropical regions has been estimated to account for just over half of cattle worldwide, yet it has not been demonstrated that sufficient similarities in the cattle exist to describe tropical cattle and, even less so, to characterize the meat from these animals. The aim of this review is to investigate the quality and nutrient composition of meat from cattle raised in the Tropics to determine if there is an axiomatic basis that would allow the definition of a concept of “tropical beef”. Tropical beef is the meat obtained from cattle raised in tropical environments, the population of which remains largely uncharacterized. Production systems in the Tropics are highly diverse but converge on the use of indigenous and Bos indicus breeds or Bos indicus-influenced crossbreeds under pasture feeding regimes. While some systems allow cattle to be slaughtered at ≤2 years of age, most often animals are ≥3 years. These production systems generally produce lean, low-yielding carcasses and tough (>46 N), lean (≤3.6% intramuscular fat) meat with a macronutrient composition otherwise similar to beef from animals raised elsewhere (72–74% moisture and 20–24% protein). Fatty acid profiles depend on the breed and production systems, while mineral content is influenced by the environment. Although lean and tough, tropical beef is highly acceptable to the consumers it serves, is culturally and traditionally relevant and, in many countries, contributes to food security. Consolidating the findings from animal and meat science studies in the Tropics has allowed the demonstration of an axiomatic basis defining “tropical beef” as a concept.

Serum proteomes of Santa Gertrudis cattle before and after infestation with Rhipicephalus australis ticks

Previous studies have applied genomics and transcriptomics to identify immune and genetic markers as key indicator traits for cattle tick susceptibility/resistance; however, results differed between breeds, and there is lack of information on the use of host proteomics. Serum samples from Santa Gertrudis cattle (naïve and phenotyped over 105 days as tick-resistant [TR] or tick-susceptible [TS]) were used to conduct differential abundance analyses of protein profiles. Serum proteins were digested into peptides followed by identification and quantification using sequential window acquisition of all instances of theoretical fragment ion mass spectrometry. Before tick infestation, abundance of 28 proteins differed significantly (adjusted P < 10^-5 ) between TR and TS. These differences were also observed following tick infestation (TR vs TS) with a further eight differentially abundant proteins in TR cattle, suggesting possible roles in adaptive responses. The intragroup comparisons (TS-0 vs TS and TR-0 vs TR) showed that tick infestation elicited quite similar responses in both groups of cattle, but with relatively stronger responses in TR cattle. Many of the significantly differentially abundant proteins in TR Santa Gertrudis cattle (before and after tick infestation) were associated with immune responses including complement factors, chemotaxis for immune cells and acute-phase responses.

Gene regulation could be attributed to TCF3 and other key transcription factors in the muscle of pubertal heifers

Puberty is a whole-body event, driven by the hypothalamic integration of peripheral signals such as leptin or IGF-1. In the process of puberty, reproductive development is simultaneous to growth, including muscle growth. To enhance our understanding of muscle function related to puberty, we performed transcriptome analyses of muscle samples from six pre- and six post-pubertal Brahman heifers (Bos indicus). Our aims were to perform differential expression analyses and co-expression analyses to derive a regulatory gene network associate with puberty. As a result, we identified 431 differentially expressed (DEx) transcripts (genes and non-coding RNAs) when comparing pre- to post-pubertal average gene expression. The DEx transcripts were compared with all expressed transcripts in our samples (over 14,000 transcripts) for functional enrichment analyses. The DEx transcripts were associated with "extracellular region," "inflammatory response" and "hormone activity" (adjusted p < .05). Inflammatory response for muscle regeneration is a necessary aspect of muscle growth, which is accelerated during puberty. The term "hormone activity" may signal genes that respond to progesterone signalling in the muscle, as the presence of this hormone is an important difference between pre- and post-pubertal heifers in our experimental design. The DEx transcript with the highest average expression difference was a mitochondrial gene, ENSBTAG00000043574 that might be another important link between energy metabolism and puberty. In the derived co-expression gene network, we identified six hub genes: CDC5L, MYC, TCF3, RUNX2, ATF2 and CREB1. In the same network, 48 key regulators of DEx transcripts were identified, using a regulatory impact factor metric. The hub gene TCF3 was also a key regulator. The majority of the key regulators (22 genes) are members of the zinc finger family, which has been implicated in bovine puberty in other tissues. In conclusion, we described how puberty may affect muscle gene expression in cattle.

Cattle adapted to tropical and subtropical environments: genetic and reproductive considerations

Efforts to understand biological functions and develop management schemes specific to Bos indicus-influenced cattle raised in tropical and subtropical environments are critical to meet the increasing global demand for protein. In the United States, B. indicus breeds are mostly used to generate B. indicus × B. taurus crosses with increased thermal and parasite tolerance, while retaining some productive characteristics of B. taurus cattle. Although crossbreeding represents a proven strategy to improve cattle adaptation almost immediately, research has also attempted to identify B. taurus genetics that can withstand subtropical and tropical climates. Reduced milk production and delayed reproductive maturation appear to be related with tropical adaptation of B. taurus breeds, as a means to conserve energy under stressful conditions and limited nutrition. Moreover, longevity may be the ultimate adaptation response to unfavorable environments, and retention of bulls and heifers from proven cows is the recommended strategy to improve longevity in B. indicus-influenced herds. Besides selection for longevity, other aspects should be considered when planning reproductive management in tropical and subtropical regions. Bos indicus and B. taurus breeds have multiple differences pertaining to reproductive function, including age at puberty, ovarian dynamics, and pregnancy development. Nutritional strategies such as the stair-step regimen, and use of exogenous progesterone (P4) inserts are options to hasten puberty attainment of late-maturing B. indicus-influenced heifers. Yet, limited pharmacological alternatives are available for reproductive management of B. indicus-influenced females in the United States, which rely on GnRH-based protocols not specifically designed to the reproductive function of B. indicus breeds. In contrast, hormonal protocols based on exogenous P4, estradiol esters, and equine chorionic gonadotropin are available for use in B. indicus females in South America. These include protocols tailored to prepubertal heifers, anestrous cows, and cycling nulliparous or parous females, which often yield pregnancy rates of 50% to fixed-time artificial insemination. The global dairy industry also faces similar challenges in increasing demand and production as the beef industry. Selection of cows capable of sustaining optimal milk yield, reproductive success, and health status in hot and humid conditions is essential for optimal dairy production in subtropical and tropical regions.

Shedding and characterization of gastrointestinal nematodes of growing beef heifers in Central Texas

•

Bos indicus crossbred females had higher fecal eggs as calves than as yearlings.

•

More Bos indicus background had higher eggs per gram of feces through 1 year of age.

•

Cooperia and Haemonchus species predominated through 1 year of heifer age.

•

Ostertagia and Trichostrongylus species were detected mostly after 1 year of age.

•

Eggs per gram of feces correlated with few production traits of heifers.

Gastrointestinal nematodes (GIN) do not detrimentally affect cattle to the extent of small ruminants. However, they are developing resistance to drugs used to treat them. Genetic strategies to control the nematodes and/or their detrimental effects could be a sustainable alternative to treatment with drugs. An essential first step in development of such a strategy is characterization of nematode populations in commonly used breed types of cattle in local conditions. Fecal egg counts (FEC) were obtained every two months on a cohort of 53 crossbred Nellore-Angus heifers grazing Central Texas pastures from an average heifer age of 3 months to approximately 2 years of age. For 10 of those 12 sets of samples, coprocultures were set up to characterize gastrointestinal nematode species present. Heifers were ½ Nellore ½ Angus (n = 18) or ¾ Angus ¼ Nellore (n = 35). They were born in the spring of 2014 to cows that were from 3–5 years old. They were maintained as a group throughout weaning, postweaning, exposure to bulls as yearlings, and as pregnant cows through the birth of their first calves. An interaction of breed group with sampling time (P < 0.0001) highlighted favorable FEC of ¼ Nellore heifers as compared to ½ Nellore in all but two sampling times. Fecal egg count means were, in general, higher for heifers of both groups in sampling times up to one year of age. Season effects on FEC may be important, but the effect of age may have obscured their detection. There were few significant correlation coefficients for FEC traits with a variety of production traits of these females. Average FEC residuals were positively correlated (r = 0.28 and 0.41; P < 0.05) with winter coat shedding score evaluated at approximately 17 and 24 months of age. Residual correlations of average FEC with calf weaning weight and incidence of shedding with calf age at weaning (r> 0.3) may be indicative of the increased susceptibility of females that lactate heavily or longer to internal parasite infection. Proportions of GIN genera by sampling day differed from χ2 expectation (P < 0.0001). Cooperia and Haemonchus species were detected in large proportions in sampling dates that corresponded to heifers less than one year of age. Ostertagia and Trichostrongylus species predominated in sampling dates after heifers reached one year of age.

Flight speed as an indicator of Zebu cattle growth performance and suitability to the feedlot environment

Context Studies assessing the relationship between cattle reactivity and growth performance, as well as suitability to the feedlot, in Zebu cattle are limited. Based on previous studies using Bos Taurus cattle, it is expected that reactivity measured at weaning could be used as a predictor of feedlot performance, and that cattle with greater flight speed (FS) would have reduced average daily gain (ADG). Aims The aims of this study were to (1) determine if FS (m/s) measured before entry into the feedlot could be used as a predictor of cattle performance, (2) evaluate the consistency of FS over time, and (3) determine when FS should be measured to best predict animal suitability to the feedlot. Methods A total of 54 Nellore and 40 Guzerat calves were assessed by using FS, which was collected at the time of (FS0) and 1 week after weaning (FS1), as well as on Day 28 (FS28), 56 (FS56) and 168 (FS168) after entry into the feedlot, and a FS average was calculated using FS28, FS56 and FS168 (FS28–168). Performance using ADG (kg/calf.day) was calculated for each animal over three periods: ADG1–28, ADG28–168 and ADG1–168. Bunk attendance (min/calf.day) was recorded over a 173-day feeding period. Key results The ADG28–168 was reduced by 0.14 kg/calf.day for each unit increase in FS0 (P &lt; 0.05). Both ADG28–168 and ADG1–168 were reduced by 0.02 kg/calf.day for each unit increase in FS1 (P ≤ 0.01), but neither FS0 nor FS1 was associated with ADG1–28 (P &gt; 0.10). A significant increase was observed between FS1 and FS28 (P &lt; 0.01), followed by a reduction in FS56 (P &lt; 0.01) and FS168 (P &lt; 0.01). Conclusion FS measured at or within 1 week of weaning has the potential to be used as an indicator of feedlot growth performance in Nellore and Guzerat cattle. Implications Although our results do not indicate the efficacy of FS as a predictor of early suitability of cattle to the feedlot, it does have the potential to support its recommendation as an indicator of overall feedlot growth performance in Nellore and Guzerat breeds.

Towards a new phenotype for tick resistance in beef and dairy cattle: a review

About 80% of the world’s cattle are affected by ticks and tick-borne diseases, both of which cause significant production losses. Cattle host resistance to ticks is the most important factor affecting the economics of tick control, but it is largely neglected in tick-control programs due to technical difficulties and costs associated with identifying individual-animal variation in resistance. The present paper reviews the scientific literature to identify factors affecting resistance of cattle to ticks and the biological mechanisms of host tick resistance, to develop alternative phenotype(s) for tick resistance. If new cost-effective phenotype(s) can be developed and validated, then tick resistance of cattle could be genetically improved using genomic selection, and incorporated into breeding objectives to simultaneously improve cattle productive attributes and tick resistance. The phenotype(s) could also be used to improve tick control by using cattle management. On the basis of the present review, it is recommended that three possible phenotypes (haemolytic analysis; measures of skin hypersensitivity reactions; simplified artificial tick infestations) be further developed to determine their practical feasibility for consistently, cost-effectively and reliably measuring cattle tick resistance in thousands of individual animals in commercial and smallholder farmer herds in tropical and subtropical areas globally. During evaluation of these potential new phenotypes, additional measurements should be included to determine the possibility of developing a volatile-based resistance phenotype, to simultaneously improve cattle resistance to both ticks and biting flies. Because the current measurements of volatile chemistry do not satisfy the requirements of a simple, cost-effective phenotype for use in commercial cattle herds, consideration should also be given to inclusion of potentially simpler measures to enable indirect genetic selection for volatile-based resistance to ticks.

Genome-wide association studies for tick resistance in Bos taurus × Bos indicus crossbred cattle: A deeper look into this intricate mechanism

Rhipicephalus (Boophilus) microplus is the main cattle ectoparasite in tropical areas. Gir × Holstein crossbred cows are well adapted to different production systems in Brazil. In this context, we performed genome-wide association study (GWAS) and post-GWAS analyses for R. microplus resistance in an experimental Gir × Holstein F₂ population. Single nucleotide polymorphisms (SNP) identified in GWAS were used to build gene networks and to investigate the breed of origin for its alleles. Tick artificial infestations were performed during the dry and rainy seasons. Illumina BovineSNP50 BeadChip (Illumina Inc., San Diego, CA) and single-step BLUP procedure was used for GWAS. Post-GWAS analyses were performed by gene ontology terms enrichment and gene transcription factors networks, generated from enriched transcription factors, identified from the promoter sequences of selected gene sets. The genetic origin of marker alleles in the F₂ population was assigned using the breed of origin of alleles approach. Heritability estimates for tick counts were 0.40 ± 0.11 in the rainy season and 0.54 ± 0.11 in the dry season. The top ten 0.5-Mbp windows with the highest percentage of genetic variance explained by SNP markers were found in chromosomes 10 and 23 for both the dry and rainy seasons. Gene network analyses allowed the identification of genes involved with biological processes relevant to immune system functions (TREM1, TREM2, and CD83). Gene-transcription factors network allowed the identification of genes involved with immune functions (MYO5A, TREML1, and PRSS16). In resistant animals, the average proportion of animals showing significant SNPs with paternal and maternal alleles originated from Gir breed was 44.8% whereas the proportion of animals with both paternal and maternal alleles originated from Holstein breed was 11.3%. Susceptible animals showing both paternal and maternal alleles originated from Holstein breed represented 44.6% on average, whereas both paternal and maternal alleles originated from Gir breed animals represented 9.3%. This study allowed us to identify candidate genes for tick resistance in Gir × Holstein crossbreds in both rainy and dry seasons. According to the origin of alleles analysis, we found that most animals classified as resistant showed 2 alleles from Gir breed, while the susceptible ones showed alleles from Holstein. Based on these results, the identified genes may be thoroughly investigated in additional experiments aiming to validate their effects on tick resistance phenotype in cattle.

Cattle Tick Rhipicephalus microplus-Host Interface: A Review of Resistant and Susceptible Host Responses

Ticks are able to transmit tick-borne infectious agents to vertebrate hosts which cause major constraints to public and livestock health. The costs associated with mortality, relapse, treatments, and decreased production yields are economically significant. Ticks adapted to a hematophagous existence after the vertebrate hemostatic system evolved into a multi-layered defense system against foreign invasion (pathogens and ectoparasites), blood loss, and immune responses. Subsequently, ticks evolved by developing an ability to suppress the vertebrate host immune system with a devastating impact particularly for exotic and crossbred cattle. Host genetics defines the immune responsiveness against ticks and tick-borne pathogens. To gain an insight into the naturally acquired resistant and susceptible cattle breed against ticks, studies have been conducted comparing the incidence of tick infestation on bovine hosts from divergent genetic backgrounds. It is well-documented that purebred and crossbred Bos taurus indicus cattle are more resistant to ticks and tick-borne pathogens compared to purebred European Bos taurus taurus cattle. Genetic studies identifying Quantitative Trait Loci markers using microsatellites and SNPs have been inconsistent with very low percentages relating phenotypic variation with tick infestation. Several skin gene expression and immunological studies have been undertaken using different breeds, different samples (peripheral blood, skin with tick feeding), infestation protocols and geographic environments. Susceptible breeds were commonly found to be associated with the increased expression of toll like receptors, MHC Class II, calcium binding proteins, and complement factors with an increased presence of neutrophils in the skin following tick feeding. Resistant breeds had higher levels of T cells present in the skin prior to tick infestation and thus seem to respond to ticks more efficiently. The skin of resistant breeds also contained higher numbers of eosinophils, mast cells and basophils with up-regulated proteases, cathepsins, keratins, collagens and extracellular matrix proteins in response to feeding ticks. Here we review immunological and molecular determinants that explore the cattle tick Rhipicephalus microplus-host resistance phenomenon as well as contemplating new insights and future directions to study tick resistance and susceptibility, in order to facilitate interventions for tick control.

Climate Adaptation of Tropical Cattle

There is sustained growth in the number of tropical cattle, which represent more than half of all cattle worldwide. By and large, most research in tropical areas is still focused on breeds of cattle, their particular advantages or disadvantages in tropical areas, and the tropical forages or feeds that could be usefully fed to them. A consistent issue for adaptation to climate is the heat of tropical environments. Changing the external characteristics of the animal, such as color and coat characteristics, is one way to adapt, and there are several major genes for these traits. However, further improvement in heat tolerance and other adaptation traits will need to use the entire genome and all physical and physiological systems. Apart from the response to heat, climate forcing through methane emission identifies dry season weight loss as an important if somewhat neglected trait in climate adaptation of cattle. The use of genome-estimated breeding values in tropical areas is in its infancy and will be difficult to implement, but will be essential for rapid, coordinated genetic improvement. The difficulty of implementation cannot be exaggerated and may require major improvements in methodology.

BREEDING AND GENETICS SYMPOSIUM: Breeding for resilience to heat stress effects in dairy ruminants. A comprehensive review

Selection for heat tolerant (HT) animals in dairy production has been so far linked to estimation of declines in production using milk recording and meteorological information on the day of control using reaction norms. Results from these models show that there is a reasonable amount of genetic variability in the individual response to high heat loads, which makes feasible selection of HT animals at low costs. However, the antagonistic relationship between level of production and response to heat stress (HS) implies that selection for HT animals under this approach must be done with caution so that productivity is not damaged. Decomposition of the genetic variability in principal components (PC) can provide selection criteria independent of milk production level although biological interpretation of PC is difficult. Moreover, given that response to heat stress for each animal is estimated with very sparse information collected under different physiological and management circumstances, biased (normally underestimation) and lack of accuracy may be expected. Alternative phenotypic characterization of HT can come from the use of physiological traits, which have also shown moderate heritability. However, costs of a large scale implementation based on physiological characteristics has precluded its use. Another alternative is the use of biomarkers that define heat tolerance. A review of biomarkers of HS from more recent studies is provided. Of particular interest are milk biomarkers, which together with infrared spectra prediction equations can provide useful tools for genetic selection. In the 'omics' era, genomics, transcriptomics, proteomics and metabolomics have been already used to detect genes affecting HT. A review of findings in these areas is also provided. Except for the slick hair gene, there are no other genes for which variants have been clearly associated with HT. However, integration of omics information could help in pointing at knots of the HS control network and, in the end, to a panel of markers to be used in the selection of HT animals. Overall, HT is a complex phenomenon that requires integration of fine phenotypes and omics information to provide accurate tools for selection without damaging productivity. Technological developments to make on-farm implementation feasible and with greater insight into the key biomarkers and genes involved in HT are needed.

Genetic parameters for tick counts across months for different tick species and anatomical locations in South African Nguni cattle

The objective of the study was to characterise genetic parameters across months for different tick species and anatomical locations in South African Nguni cattle. Tick counts were conducted monthly, over a 2-year period, on 586 Nguni cattle under natural infestation, from four herds located in different provinces of South Africa. The counts were recorded for six species of ticks (Amblyomma hebraeum, Rhipicephalus evertsi evertsi, Rhipicephalus decoleratus and microplus (Boofilids), Rhipicephalus appendiculatus, Rhipicephalus simus and Hyalomma marginatum) attached on eight anatomical locations on the animals and were summed by species and anatomical location. Heritability estimates, phenotypic and genetic correlations were estimated on a monthly basis using mixed linear models, fitting univariate and bivariate sire models. Fixed effects considered were location, sex, year and age as a covariate. Tick counts were higher in the hot months, and A. hebraeum was the most dominant tick species. Heritability estimates for tick count varied by month and trait and ranged from 0 to 0.89. Genetic correlations were mostly positive, and low to high, with some negative correlations with high standard error. Phenotypic correlations were low to moderate. In general, high genetic correlations were observed between whole body count and the anatomical location counts, suggesting that it may not be necessary to conduct whole body counts. Counts from the belly and perineum appeared to be the most suitable surrogate traits for whole body count. These findings provide useful information for developing strategies for the practical implementation of genetic selection, as a supplement to the traditional tick control measures.

Cattle phenotypes can disguise their maternal ancestry

Background

Cattle are bred for, amongst other factors, specific traits, including parasite resistance and adaptation to climate. However, the influence and inheritance of mitochondrial DNA (mtDNA) are not usually considered in breeding programmes. In this study, we analysed the mtDNA profiles of cattle from Victoria (VIC), southern Australia, which is a temperate climate, and the Northern Territory (NT), the northern part of Australia, which has a tropical climate, to determine if the mtDNA profiles of these cattle are indicative of breed and phenotype, and whether these profiles are appropriate for their environments.

Results

A phylogenetic tree of the full mtDNA sequences of different breeds of cattle, which were obtained from the NCBI database, showed that the mtDNA profiles of cattle do not always reflect their phenotype as some cattle with Bos taurus phenotypes had Bos indicus mtDNA, whilst some cattle with Bos indicus phenotypes had Bos taurus mtDNA. Using D-loop sequencing, we were able to contrast the phenotypes and mtDNA profiles from different species of cattle from the 2 distinct cattle breeding regions of Australia. We found that 67 of the 121 cattle with Bos indicus phenotypes from NT (55.4%) had Bos taurus mtDNA. In VIC, 92 of the 225 cattle with Bos taurus phenotypes (40.9%) possessed Bos indicus mtDNA. When focusing on oocytes from cattle with the Bos taurus phenotype in VIC, their respective oocytes with Bos indicus mtDNA had significantly lower levels of mtDNA copy number compared with oocytes possessing Bos taurus mtDNA (P < 0.01). However, embryos derived from oocytes with Bos indicus mtDNA had the same ability to develop to the blastocyst stage and the levels of mtDNA copy number in their blastocysts were similar to blastocysts derived from oocytes harbouring Bos taurus mtDNA. Nevertheless, oocytes originating from the Bos indicus phenotype exhibited lower developmental potential due to low mtDNA copy number when compared with oocytes from cattle with a Bos taurus phenotype.

Conclusions

The phenotype of cattle is not always related to their mtDNA profiles. MtDNA profiles should be considered for breeding programmes as they also influence phenotypic traits and reproductive capacity in terms of oocyte quality.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-017-0523-5) contains supplementary material, which is available to authorized users.

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

Introduction

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

Materials and Methods

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

Results

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

Conclusions

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

Genome-wide association study of tick resistance in South African Nguni cattle

Ticks and tick-borne diseases are among the main causes of economic loss in the South African cattle industry through high morbidity and mortality rates. Concerns of the general public regarding chemical residues may tarnish their perceptions of food safety and environmental health when the husbandry of cattle includes frequent use of acaricides to manage ticks. The primary objective of this study was to identify single nucleotide polymorphism (SNP) markers associated with host resistance to ticks in South African Nguni cattle. Tick count data were collected monthly from 586 Nguni cattle reared in four herds under natural grazing conditions over a period of two years. The counts were recorded for six species of ticks attached in eight anatomical locations on the animals and were summed by species and anatomical location. This gave rise to 63 measured phenotypes or traits, with results for 12 of these traits being reported here. Tick count (x) data were transformed using log10(x+1) and the resulting values were examined for normality. DNA was extracted from hair and blood samples and was genotyped using the Illumina BovineSNP50 assay. After quality control (call rate >90%, minor allele frequency >0.02), 40,436 SNPs were retained for analysis. Genetic parameters were estimated and association analysis for tick resistance was carried out using two approaches: a genome-wide association (GWA) analysis using the GenABEL package and a regional heritability mapping (RHM) analysis. The Bonferroni genome-wide (P<0.05) corrected significance threshold was 1.24×10(-6), with 2.47×10(-5) as the suggestive significance threshold (P<0.10) (i.e., one false positive per genome scan) in the GWA analysis. Likelihood ratio test (LRT) thresholds for genome-wide and suggestive significance were 13.5 and 9.15 for the RHM analysis. Six ixodid tick species were identified, with Amblyomma hebraeum (the vector for Heartwater disease) being the dominant species. Heritability estimates (h(2)) from the fitted animal and sire models ranged from 0.02±0.00 to 0.17±0.04 for the transformed tick count data. Several genomic regions harbouring quantitative trait loci (QTL) were identified for different tick count traits by both the GWA and RHM approaches. Three genome-wide significant regions on chromosomes 7, 10 and 19 were identified for total tick count on the head, total body A. hebraeum tick count and total A. hebraeum on the perineum region, respectively. Additional regions significant at the suggestive level were identified on chromosomes 1, 3, 6, 7, 8, 10, 11, 12, 14, 15, 17, 19 and 26 for several of the traits. The GWA approach identified more genomic regions than did the RHM approach. The chromosomal regions identified here as harbouring QTL underlying variation in tick burden form the basis for further analyses to identify specific candidate genes and polymorphisms related to cattle tick resistance and provide the potential for marker-assisted selection in Nguni cattle.

Accuracy of genomic selection for age at puberty in a multi-breed population of tropically adapted beef cattle

Genomic selection is becoming a standard tool in livestock breeding programs, particularly for traits that are hard to measure. Accuracy of genomic selection can be improved by increasing the quantity and quality of data and potentially by improving analytical methods. Adding genotypes and phenotypes from additional breeds or crosses often improves the accuracy of genomic predictions but requires specific methodology. A model was developed to incorporate breed composition estimated from genotypes into genomic selection models. This method was applied to age at puberty data in female beef cattle (as estimated from age at first observation of a corpus luteum) from a mix of Brahman and Tropical Composite beef cattle. In this dataset, the new model incorporating breed composition did not increase the accuracy of genomic selection. However, the breeding values exhibited slightly less bias (as assessed by deviation of regression of phenotype on genomic breeding values from the expected value of 1). Adding additional Brahman animals to the Tropical Composite analysis increased the accuracy of genomic predictions and did not affect the accuracy of the Brahman predictions.

Host resistance in cattle to infestation with the cattle tick Rhipicephalus microplus

Resistance to Rhipicephalus microplus infestation in cattle has many effector mechanisms, each of which is likely to be modulated by complex, interacting factors. Some of the mechanisms of host resistance and their modulating factors have been identified and quantified, although much remains to be explained. The variation in resistance to tick infestation is most marked between Bos taurus and Bos indicus cattle, taurine cattle given the same exposure carrying between five and 10 times as many ticks as indicine cattle. Tick resistance is mostly manifest against attaching larvae, which attempt to feed often and without success, death occurring mostly within 24 h of finding a host. There is evidence of innate and adaptive immune response to tick infestation, and it appears that the relative importance of each differs between indicine and taurine cattle. There is conflicting information regarding the role of humoral immunity in tick resistance, and recent studies indicate that strong IgG responses to tick antigens are not protective. A strong T-cell-mediated response directed against larval stages, as mounted by indicine cattle, seems to be protective. Variation in the extracellular matrix of skin (epidermal growth factors, collagens and other matrix components such as lumican) also contributes to variation in host resistance.

Genetic association between temperament and sexual precocity indicator traits in Nellore cattle

The aim of this study was to estimate the genetic associations between temperament and sexual precocity indicators in Nellore cattle. Temperament was evaluated at approximately 550 days of age, using three traits: i) movement score (MOV), recording the animals' movement inside the crush; ii) temperament score (TS) for assessing the reactions of animals in a corral pen; and iii) flight speed test (FS), recording the speed (in m/s) at which the animals exit the crush after being weighed. Three sexual precocity indicators were used: i) age at first calving (AFC, in days); ii) occurrence of precocious pregnancy (OPP, as a binary trait); and iii) yearling scrotal circumference (SC, in cm). The (co)variance and genetic parameters were estimated by Bayesian Inference via Gibbs sampling, assuming a linear animal model for FS, AFC and SC and a threshold animal model for MOV, TS and OPP in multitrait analyses. The mean posterior heritability estimates for MOV, FS, TS, AFC, OPP and SC were 0.11 ± 0.03, 0.27 ± 0.07, 0.16 ± 0.09, 0.09 ± 0.01, 0.44 ± 0.06 and 0.43 ± 0.02, respectively. The mean posterior correlation estimates were: 0.13 ± 0.14 (MOV-AFC), -0.03 ± 0.16 (MOV-OPP), 0.08 ± 0.10 (MOV-SC), 0.14 ± 0.11 (FS-AFC), -0.19 ± 0.15 (FS-OPP), -0.07 ± 0.08 (FS-SC), 0.09 ± 0.09 (TS-AFC), -0.03 ± 0.09 (TS-OPP) and -0.28 ± 0.08 (TS-SC). These results indicate that all the traits had sufficient genetic variability to respond to selection. The low magnitude of the genetic correlations indicated that selection for improvements in sexual precocity would not change the temperament of Nellore cattle. In order to obtain genetic progress in cattle temperament, direct selection for these traits is recommended.

The genetic architecture of climatic adaptation of tropical cattle

Adaptation of global food systems to climate change is essential to feed the world. Tropical cattle production, a mainstay of profitability for farmers in the developing world, is dominated by heat, lack of water, poor quality feedstuffs, parasites, and tropical diseases. In these systems European cattle suffer significant stock loss, and the cross breeding of taurine x indicine cattle is unpredictable due to the dilution of adaptation to heat and tropical diseases. We explored the genetic architecture of ten traits of tropical cattle production using genome wide association studies of 4,662 animals varying from 0% to 100% indicine. We show that nine of the ten have genetic architectures that include genes of major effect, and in one case, a single location that accounted for more than 71% of the genetic variation. One genetic region in particular had effects on parasite resistance, yearling weight, body condition score, coat colour and penile sheath score. This region, extending 20 Mb on BTA5, appeared to be under genetic selection possibly through maintenance of haplotypes by breeders. We found that the amount of genetic variation and the genetic correlations between traits did not depend upon the degree of indicine content in the animals. Climate change is expected to expand some conditions of the tropics to more temperate environments, which may impact negatively on global livestock health and production. Our results point to several important genes that have large effects on adaptation that could be introduced into more temperate cattle without detrimental effects on productivity.

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.

Genetics of early and lifetime annual reproductive performance in cows of two tropical beef genotypes in northern Australia

Reproduction records from 2137 cows first mated at 2 years of age and recorded through to 8.5 years of age were used to study the genetics of early and lifetime reproductive performance from two genotypes (1020 Brahman and 1117 Tropical Composite) in tropical Australian production systems. Regular ultrasound scanning of the reproductive tract, coupled with full recording of mating, calving and weaning histories, allowed a comprehensive evaluation of a range of reproductive traits. Results showed components traits of early reproductive performance had moderate to high heritabilities, especially in Brahmans. The heritability of lactation anoestrous interval in 3-year-old cows was 0.51 ± 0.18 and 0.26 ± 0.11 for Brahman and Tropical Composite, respectively. Heritabilities of binary reproductive output traits (conception rate, pregnancy rate, calving rate and weaning rate) from first and second matings were generally moderate to high on the underlying scale. Estimates ranged from 0.15 to 0.69 in Brahman and 0.15 to 0.34 in Tropical Composite, but were considerably lower when expressed on the observed scale, particularly for those traits with high mean levels. Heritabilities of lifetime reproduction traits were low, with estimates of 0.11 ± 0.06 and 0.07 ± 0.06 for lifetime annual weaning rate in Brahman and Tropical Composite, respectively. Significant differences in mean reproductive performance were observed between the two genotypes, especially for traits associated with anoestrus in first-lactation cows. Genetic correlations between early-in-life reproductive measures and lifetime reproduction traits were moderate to high. Genetic correlations between lactation anoestrous interval and lifetime annual weaning rate were –0.62 ± 0.24 in Brahman and –0.87 ± 0.32 in Tropical Composite. The results emphasise the substantial opportunity that exists to genetically improve weaning rates in tropical beef cattle breeds by focusing recording and selection on early-in-life female reproduction traits, particularly in Brahman for traits associated with lactation anoestrus.

Multi-trait assessment of early-in-life female, male and genomic measures for use in genetic selection to improve female reproductive performance of Brahman cattle

Early-in-life female and male measures with potential to be practical genetic indicators were chosen from earlier analyses and examined together with genomic measures for multi-trait use to improve female reproduction of Brahman cattle. Combinations of measures were evaluated on the genetic gains expected from selection of sires and dams for each of age at puberty (AGECL, i.e. first observation of a corpus luteum), lactation anoestrous interval in 3-year-old cows (LAI), and lifetime annual weaning rate (LAWR, i.e. the weaning rate of cows based on the number of annual matings they experienced over six possible matings). Selection was on an index of comparable records for each combination. Selection intensities were less than theoretically possible but assumed a concerted selection effort was able to be made across the Brahman breed. The results suggested that substantial genetic gains could be possible but need to be confirmed in other data. The estimated increase in LAWR in 10 years, for combinations without or with genomic measures, ranged from 8 to 12 calves weaned per 100 cows from selection of sires, and from 12 to 15 calves weaned per 100 cows from selection of sires and dams. Corresponding reductions in LAI were 60–103 days or 94–136 days, and those for AGECL were 95–125 or 141–176 days, respectively. Coat score (a measure of the sleekness or wooliness of the coat) and hip height in females, and preputial eversion and liveweight in males, were measures that may warrant wider recording for Brahman female reproduction genetic evaluation. Pregnancy-test outcomes from Matings 1 and 2 also should be recorded. Percentage normal sperm may be important to record for reducing LAI and scrotal size and serum insulin-like growth factor-I concentration in heifers at 18 months for reducing AGECL. Use of a genomic estimated breeding value (EBV) in combination with other measures added to genetic gains, especially at genomic EBV accuracies of 40%. Accuracies of genomic EBVs needed to approach 60% for the genomic EBV to be the most important contributor to gains in the combinations of measures studied.

Genetic relationships of female reproduction with growth, body composition, maternal weaning weight and tropical adaptation in two tropical beef genotypes

The genetic relationships of female reproduction with growth and body composition, tropical adaptation traits and maternal weaning weight (descriptive of genetic potential milk production) were estimated in 1027 Brahman (BRAH) and 1132 Tropical Composite (TCOMP) females. Female reproduction was evaluated at puberty, as outcomes of the first and second annual mating periods (Mating 1 and Mating 2, which commenced when females averaged 27 and 39 months of age, respectively), as well as annual averages over up to six matings. Traits evaluated included age at puberty, Mating 1 and 2 pregnancy rate, weaning rate and days to calving, and lifetime annual calving and weaning rate. Traits describing growth and body composition (liveweight, hip height, ultrasound-scanned P8 fat depth and eye muscle area, subjective body condition score and blood IGF-I concentration) were measured in the animals as heifers (at ~18 months of age), and again at the start of Mating 2. Traits describing tropical adaptation included coat-length scores in both genotypes and, in BRAH, buffalo fly lesion scores. Previously reported analyses of these data identified heifer IGF-I and coat and buffalo-fly-lesion scores as potential genetic indicators for age at puberty in BRAH. The results of the present study found that exploiting these relationships would have no unfavourable genetic consequences for later female reproduction and, in some cases, may be indicators of female reproduction, when evaluated as outcomes of Matings 1 or 2, or as lifetime annual calving or weaning rates. For BRAH, heifer liveweight was a genetic indicator for Mating 1 weaning rate (rg = 0.70), and, while standard errors were high, there were also positive genetic correlations of heifer hip height, eye muscle area and blood IGF-I concentration with Mating 1 weaning rate (rg = 0.61, 0.58 and 0.43, respectively). For TCOMP, significant genetic relationships of heifer growth, body composition and tropical adaptation traits with female reproduction were virtually absent, suggesting that there is less opportunity to identify earlier in life measures as genetic indicators of reproduction for this genotype. Higher maternal weaning weight was significantly genetically related to lower lifetime annual weaning rate (rg = –0.50) in BRAH, and with lower Mating 2 calving and weaning rate (rg = –0.72 and –0.59, respectively) in TCOMP, which will need to be considered when making selection decisions that affect genetic milk in these genotypes. Importantly, the results presented revealed no strong genetic antagonisms of heifer growth and body composition traits with female reproduction, suggesting that selection could be undertaken to improve these simultaneously.

Genome-wide association study for inhibin, luteinizing hormone, insulin-like growth factor 1, testicular size and semen traits in bovine species

The fertility of young bulls impacts on reproduction rates, farm profit and the rate of genetic progress in beef herds. Cattle researchers and industry therefore routinely collect data on the reproductive performance of bulls. Genome-wide association studies were carried out to identify genomic regions and genes associated with reproductive traits measured during the pubertal development of Tropical Composite bulls, from 4 to 24 months of age. Data from 1 085 bulls were collected for seven traits: blood hormone levels of inhibin at 4 months (IN), luteinizing hormone following a gonadotropin releasing hormone challenge at 4 months (LH), insulin-like growth factor 1 at 6 months (IGF1), scrotal circumference at 12 months (SC), sperm motility at 18 months (MOT), percentage of normal spermatozoa at 24 months (PNS) and age at a scrotal circumference of 26 cm (AGE26, or pubertal age). Data from 729 068 single-nucleotide polymorphisms were used in the association analysis. Significant polymorphism associations were discovered for IN, IGF1, SC, AGE26 and PNS. Based on these associations, INHBE, INHBC and HELB are proposed as candidate genes for IN regulation. Polymorphisms associated with IGF1 mapped to the PLAG1 gene region, validating a reported quantitative trait locus on chromosome 14 for IGF1. The X chromosome contained most of the significant associations found for SC, AGE26 and PNS. These findings will contribute to the identification of diagnostic genetic markers and informed genomic selection strategies to assist breeding of cattle with improved fertility. Furthermore, this work provides evidence contributing to gene function annotation in the context of male fertility.

Breeding strategies for tick resistance in tropical cattle: a sustainable approach for tick control

About 80 % of world cattle population is under the risk of ticks and tick borne diseases (TTBDs). Losses caused by bovine tick burdens in tropical countries have a tremendous economic impact on production systems. Chemical control of disease has been found to be ineffective and also involving large cost. To reduce our reliance on these chemical products, it is necessary to embark on programs that include habitat management, genetic selection of hosts, and development of a strain capable of inducing host resistance to ticks. Selection for disease resistance provide alternate method for sustainable control of TTBDs. Domestic livestock manifests tick-resistance by skin thickness, coat type, coat color, hair density and skin secretions etc. Zebu cattle have, on average, greater tick resistance than either European cattle or African cattle. Heritability for tick burden in cattle has been shown to range about 0.30, which is sufficient to result in the success of some programs of selection for tick resistance in cattle. To select animals at younger age, to reduce generation interval and to increase genetic gain, marker assisted selection is an important tool. There are also various MHC molecules which are associated with resistance to TTBDs.

Male traits and herd reproductive capability in tropical beef cattle. 2. Genetic parameters of bull traits

A total of 4063 young bulls of two tropical genotypes (1639 Brahman and 2424 Tropical Composite) raised in northern Australia were evaluated for a comprehensive range of production and reproduction traits up to 24 months of age. Prior to weaning, peripheral blood concentrations of luteinising hormone (LH) and inhibin were measured at 4 months of age. At weaning (6 months) blood insulin-like growth factor-1 (IGF-I) and flight time were recorded. Body composition traits of fat depth and eye-muscle area were determined by ultrasonography at 15 months of age when additional measurements of liveweight, hip height and body condition score were recorded. Bull breeding soundness was evaluated at ~12, 18 and 24 months of age when measurements of scrotal circumference, sheath score, semen mass activity, progressive motility of individual sperm and percent morphologically normal sperm were recorded. Magnitude of heritability and genetic correlations changed across time for some traits. Heritability of LH, inhibin, IGF-I and of 18-month scrotal circumference, mass activity, progressive motility and percent normal sperm was 0.31, 0.74, 0.44, 0.75, 0.24, 0.15 and 0.25, respectively, for Brahmans and 0.48, 0.72, 0.36, 0.43, 0.13, 0.15 and 0.20, respectively, for Tropical Composites. Inhibin and IGF-I had moderate genetic association with percent normal sperm at 24 months in Brahmans but low to negligible associations in Tropical Composites. Body condition score in Brahmans and sperm motility (mass and individual) traits in both genotypes had moderate to strong genetic correlation with percent normal sperm and may prove useful candidates for indirect selection. There is scope to increase scrotal circumference by selection and this will be associated with favourable correlated responses of improved semen quality in both genotypes. The lack of genetic antagonism among bull traits indicates that selection for improved semen quality will not adversely affect other production traits.

Male traits and herd reproductive capability in tropical beef cattle. 1. Experimental design and animal measures

Research into the genetics of whole herd profitability has been a focus of the Beef Cooperative Research Centre for Beef Genetic Technologies over the past decade and it has been identified that measures of male reproduction may offer a potential indirect means of selecting for improved female reproduction. This paper describes the experimental design and provides a descriptive analysis of an array of male traits in Brahman and Tropical Composite genotypes managed under the medium to high stress, semi-extensive to extensive production systems of northern Australia. A total of 1639 Brahman and 2424 Tropical Composite bulls with known pedigrees, bred and raised in northern Australia, were evaluated for a comprehensive range of productive and reproductive traits. These included blood hormonal traits (luteinising hormone, inhibin and insulin-like growth factor-I); growth and carcass traits (liveweight, body condition score, ultrasound scanned 12–13th rib fat, rump P8 fat, eye muscle area and hip height); adaptation traits (flight time and rectal temperature); and a bull breeding soundness evaluation (leg and hoof conformation, sheath score, length of everted prepuce, penile anatomy, scrotal circumference, semen mass activity, sperm motility and sperm morphology). Large phenotypic variation was evident for most traits, with complete overlap between genotypes, indicating that there is likely to be a significant opportunity to improve bull fertility traits through management and bull selection.

Genetic parameters of infectious bovine keratoconjunctivitis and its relationship with weight and parasite infestations in Australian tropical Bos taurus cattle

Background

Infectious bovine keratoconjunctivitis (IBK) or ‘pinkeye’ is an economically important ocular disease that significantly impacts animal performance. Genetic parameters for IBK infection and its genetic and phenotypic correlations with cattle tick counts, number of helminth (unspecified species) eggs per gram of faeces and growth traits in Australian tropically adapted Bos taurus cattle were estimated.

Methods

Animals were clinically examined for the presence of IBK infection before and after weaning when the calves were 3 to 6 months and 15 to 18 months old, respectively and were also recorded for tick counts, helminth eggs counts as an indicator of intestinal parasites and live weights at several ages including 18 months.

Results

Negative genetic correlations were estimated between IBK incidence and weight traits for animals in pre-weaning and post-weaning datasets. Genetic correlations among weight measurements were positive, with moderate to high values. Genetic correlations of IBK incidence with tick counts were positive for the pre-weaning and negative for the post-weaning datasets but negative with helminth eggs counts for the pre-weaning dataset and slightly positive for the post-weaning dataset. Genetic correlations between tick and helminth eggs counts were moderate and positive for both datasets. Phenotypic correlations of IBK incidence with helminth eggs per gram of faeces were moderate and positive for both datasets, but were close to zero for both datasets with tick counts.

Conclusions

Our results suggest that genetic selection against IBK incidence in tropical cattle is feasible and that calves genetically prone to acquire IBK infection could also be genetically prone to have a slower growth. The positive genetic correlations among weight traits and between tick and helminth eggs counts suggest that they are controlled by common genes (with pleiotropic effects). Genetic correlations between IBK incidence and tick and helminth egg counts were moderate and opposite between pre-weaning and post-weaning datasets, suggesting that the environmental and (or) maternal effects differ between these two growth phases. This preliminary study provides estimated genetic parameters for IBK incidence, which could be used to design selection and breeding programs for tropical adaptation in beef cattle.

The RIPK2 gene: a positional candidate for tick burden supported by genetic associations in cattle and immunological response of knockout mouse

Ticks and tick-borne diseases have a detrimental impact on livestock production causing estimated losses of around $200 million per year in Australia alone. Host resistance to ticks is heritable, within-breed heritability estimates being around 0.35, and with large differences between breeds. Previously a QTL for tick burden was detected on BTA14 at ~72 Mb distal to the centromere, near the gene receptor-interacting serine-threonine kinase 2 (RIPK2). To identify polymorphisms in this region, we sequenced all exons of the RIPK2 gene, identifying 46 single nucleotide polymorphism (SNP). Using SNP from RIPK2 as well as SNP from the bovine genome sequence, we genotyped two samples, one of 1,122 taurine dairy cattle and one of 761 zebu and zebu composite beef cattle. We confirmed that SNP and haplotypes from this region, including from RIPK2, were associated with tick burden in both dairy and beef cattle. To determine whether RIPK2 influences response to tick salivary gland extract (SGE), an immunisation experiment with tick SGE in a RIPK2 knockout (RIPK2 −/−) mouse strain was conducted. There was a significant (P < 0.05) reduction in IgG production in the RIPK2 −/− mouse in response to the SGE compared to its background strain C57BL/ 6 as well as the outbred CD1 mouse strain. In addition, antibodies generated by RIPK2 −/− mice recognised a different set of antigens within SGE when compared to parental-derived antibodies. In summary, the SNP association with tick burden at BTA14 was confirmed and quantitative and qualitative differences in antibody production were observed between RIPK2 −/− and wild-type mice.

Genetic associations between flight speed and growth traits in Nellore cattle

The aim of the present study was to estimate genetic parameters for flight speed and its association with growth traits in Nellore beef cattle. The flight speed (FS) of 7,402 yearling animals was measured, using a device composed of a pair of photoelectric cells. Time interval data (s) were converted to speed (m/s) and faster animals were regarded as more reactive. The growth traits analyzed were weaning weight (WW), ADG from weaning to yearling age, and yearling scrotal circumference (SC). The (co)variance components were estimated using REML in a multitrait analysis applying an animal model. The model included random direct additive genetic and residual effects, fixed effects of contemporary groups, age of dam (classes), and age of animal as covariable. For WW, the model also included maternal genetic and permanent environmental random effects. The direct heritability estimate for FS was 0.26 ± 0.05 and direct heritability estimates for WW, SC, and ADG were 0.30 ± 0.01, 0.48 ± 0.02, and 0.19 ± 0.01, respectively. Estimates of the genetic correlation between FS and the growth traits were -0.12 ± 0.07 (WW), -0.13 ± 0.08 (ADG), and -0.11 ± 0.07 (SC). Although the values were low, these correlations showed that animals with better temperaments (slower FS) tended to present better performance. It is possible to infer that longterm selection for weight and scrotal circumference can promote a positive genetic response in the temperament of animals. Nevertheless, to obtain faster genetic progress in temperament, it would be necessary to perform direct selection for such trait. Flight speed is an easily measured indicator of temperament and can be included as a selection criterion in breeding programs for Nellore cattle.

The IGF1 pathway genes and their association with age of puberty in cattle

Insulin-like growth factor I (somatomedin C) (IGF1) influences gonadotrophin-releasing hormone (GnRH) neurons during puberty, and GnRH release guides pubertal development. Therefore, genes of the IGF1 pathway are biological candidates for the identification of single-nucleotide polymorphisms (SNPs) affecting age of puberty. In a genome-wide association study, genotyped heifers were Tropical Composite (TCOMP, n = 866) or Brahman (BRAH, n = 843), with observation of age at first corpus luteum defining puberty. We examined SNPs in or near genes of the IGF1 pathway and report seven genes associated with age at puberty in cattle: IGF1R, IGFBP2, IGFBP4, PERK (HUGO symbol EIF2AK3), PIK3R1, GSK3B and IRS1. SNPs in the IGF1 receptor (IGF1R) showed the most promising associations: two SNPs were associated with puberty in TCOMP (P < 0.05) and one in BRAH (P = 0.00009). This last SNP explained 2% of the genetic variation (R(2) = 2.04%) for age of puberty in BRAH. Hence, IGF1R was examined further. Additional SNPs were genotyped, and haplotypes were analysed. To test more SNPs in this gene, four new SNPs from dbSNP were selected and genotyped. Single SNP and haploytpe analysis revealed associations with age of puberty in both breeds. There were two haplotypes of 12 IGF1R SNPs associated with puberty in BRAH (P < 0.05) and one in TCOMP (P < 0.05). One haplotype of two SNPs was associated (P < 0.01) with puberty in BRAH, but not in TCOMP. In conclusion, the IGF1 pathway appeared more relevant for age of puberty in Brahman cattle, and IGF1R showed higher significance when compared with other genes from the pathway.