Exploring genomic inbreeding and selection signatures in a commercial Brangus herd through functional annotation

Composite breeds, including Brangus, are widely utilized in subtropical and tropical regions to harness the advantages of both Bos t. taurus and Bos t. indicus breeds. The formation and subsequent selection of composite breeds may result in discernible signatures of selection and shifts in genomic population structure. The objectives of this study were to 1) assess genomic inbreeding, 2) identify signatures of selection, 3) assign functional roles to these signatures in a commercial Brangus herd, and 4) contrast signatures of selection between selected and non-selected cattle from the same year. A total of 4035 commercial Brangus cattle were genotyped using the GGP-F250K array. Runs of Homozygosity (ROH) were used to identify signatures of selection and calculate genomic inbreeding. Quantitative trait loci (QTL) enrichment analysis and literature search identified phenotypic traits linked to ROH islands. Genomic inbreeding averaged 5%, primarily stemming from ancestors five or more generations back. A total of nine ROH islands were identified, QTL enrichment analysis revealed traits related to growth, milk composition, carcass, reproductive, and meat quality traits. Notably, the ROH island on BTA14 encompasses the pleiomorphic adenoma (PLAG1) gene, which has been linked to growth, carcass, and reproductive traits. Moreover, ROH islands associated with milk yield and composition were more pronounced in selected replacement heifers of the population, underscoring the importance of milk traits in cow-calf production. In summary, our research sheds light on the changing genetic landscape of the Brangus breed due to selection pressures and reveals key genomic regions impacting production traits.

Transcription factors and miRNA act as contrary regulators of gene expression in the testis and epididymis of Bos indicus animals

Spermatogenesis is highly conserved among mammalians, but its gene expression and regulatory profile are not entirely known. As transcription factors (TFs) and miRNAs are crucial for gene expression regulation, identifying genes negatively regulated by miRNAs and positively regulated by TFs in the testis and epididymis can provide a deeper understanding of gene expression and regulatory patterns. To do this, we used expression data coming from RNA-Seq and miRNA-Seq experiments made with biopsies from testicular parenchyma, head of the epididymis, and tail of the epididymis of four Brahman bulls. We identified miRNA differentially expressed (DE) by comparing the three distinct tissues. A co-expression analysis combined with a regulatory impact factor approach identified miRNAs and TFs with regulatory impact over gene expression regulation in the Bos indicus tissues studied. We identified 116 DE miRNAs, 206 miRNAs and 237 TFs with a significant regulatory impact on mRNA patterns in the tissues' comparisons. bta-miR-196b was the only DE miRNA for all tissue comparisons and it may be a regulator of spermatogenesis through its links with adipogenesis and insulin biosynthesis. DE genes and TFs involved in contrary regulations between the epididymis head and testis parenchyma were associated with spermatogenesis, sexual reproduction, and sperm motility. Our results provide possible mechanisms, governed by the contrary effect of miRNA and TF, leading to the differential expression between the studied tissues. We have demonstrated that our predictions of miRNAs and TFs co-regulations over target DE genes can retrieve known regulatory mechanisms and predict new ones that merit further validation.

MicroRNA breed and parent-of-origin effects provide insights into biological pathways differentiating cattle subspecies in fetal liver

Introduction: MicroRNAs (miRNAs) play a crucial role in regulating gene expression during key developmental processes, including fetal development. Brahman (Bos taurus indicus) and Angus (Bos taurus taurus) cattle breeds represent two major cattle subspecies with strikingly different phenotypes. Methods: We analyzed miRNA expression in liver samples of purebred and reciprocal crosses of Angus and Brahman to investigate breed and parent-of-origin effects at the onset of accelerated fetal growth. Results: We identified eight novel miRNAs in fetal liver samples and 14 differentially expressed miRNAs (DEMs) between purebred samples. Correlation of gene expression modules and miRNAs by breed and parent-of-origin effects revealed an enrichment of genes associated with breed-specific differences in traits such as heat tolerance (Brahman) and fat deposition (Angus). We demonstrate that genes predicted to be targets of DEMs were more likely to be differentially expressed than non-targets (p-value < 0.05). We identified several miRNAs (bta-miR-187, bta-miR-216b, bta-miR-2284c, bta-miR-2285c, bta-miR-2285cp, bta-miR-2419-3p, bta-miR-2419-5p, and bta-miR-11984) that showed similar correlation patterns as bta-miR-2355-3p, which has been associated with the glutamatergic synapse pathway, a key facilitator of heat tolerance. Furthermore, we report Angus-breed-specific miRNAs (bta-miR-2313-5p, btamiR-490, bta-miR-2316, and bta-miR-11990) that may be involved in fat deposition. Finally, we showed that the DEMs identified in fetal liver are involved in Rap1, MAPK, and Ras signalling pathways, which are important for fetal development, muscle development and metabolic traits such as fat metabolism. Conclusion: Our work sheds light on the miRNA expression patterns that contribute to gene expression differences driving phenotypic differences in indicine and taurine cattle.

Viperin (RSAD2) gene expression in peripheral blood mononuclear cells of pregnant crossbred beef cows is altered by Bos indicus genetics

The expression of interferon (IFN) stimulated genes (ISGs) in lymphocytes has been used for pregnancy diagnosis in cattle. However, among-cow variability has yielded sub-optimal predictive accuracy. We hypothesized that the expression of ISGs (ISG15, OAS1, RSAD2, CLEC3B, and AKR1B1) in early pregnancy varies according to the proportion of Bos indicus (B. indicus) genetics on females. Multiparous cows were classified in three genetic groups, High Angus (HA; n = 45 [0-33% Brahman influence]), Angus-Brahman (AB; n = 30 [34-67%]), and High Brahman (HB; n = 19 [68-100%]) and submitted to a Select-Synch + CIDR protocol. Cows that displayed estrus (n = 94) were artificially inseminated (Day0; D0). On D19, blood samples were collected to obtain peripheral blood mononuclear cells (PBMC) and measure progesterone (P4) concentrations. On D30, pregnancy diagnosis was performed. The expression of RSAD2 in PBMC of pregnant cows was positively related to the proportion of B. indicus genetics of the groups, but not the expression of ISG15 and OAS1. In pregnant cows, the proportion of B. indicus genetics was negatively associated to circulating levels of P4 concentrations. The P4 concentrations were related positively with RSAD2 expression. ROC curve results determined that for cattle with B. indicus genetics lower than 67%, the CLEC3B and AKR1B1 combination was the most accurate option to predict the outcome of pregnancy. In cows with more than 68% of B. indicus genetics, RSAD2 provided the best accuracy. In conclusion, there is a relationship between the proportion of B. indicus genetics and the ISGs gene expression in PBMC during pregnancy.

Comparison of Telomere Length in Age-Matched Primiparous and Multiparous Brahman Cows

Physiological and psychological stressors have been associated with the attrition of telomeres, which are the protective caps of chromosomes. This study compares the telomere length (TL) in 4-year-old Brahman cows grouped by the first parity (n = 8) and the second parity (n = 11). The cows were bled via jugular venipuncture, weighed, and had their body condition scores recorded at Day -28 prior to calving and at Day + 7 and Day + 28 post-calving. The duration of labor (Dlabor) and parturition ease were recorded. The peripheral leukocytes were isolated, the leukocyte blood count with differential was recorded, and the genomic DNA was extracted. The relative quantity of telomere products, which is proportional to the average TL, was determined via multiplex quantitative PCR using the ratio (T/S ratio) of bovine telomere and β-globulin DNA. Standards of the bovine telomere (10¹²-10⁷ dilution series) and β-globulin (10⁹-10⁴ dilution series) genes were utilized to produce relative copy numbers. The samples were assayed in triplicate and were included if the triplicate C_q difference was less than 0.25 cycles. The parity was the fixed effect, and the random effects included the sire and day repeated with the cow as the subject. Statistical significance was not observed in the leukocyte number or type (p > 0.1). A reduction in the TL of approximately 9225 telomeric copies was found between Parity 1 and Parity 2 (p = 0.02). A trend was found between the TL and Dlabor (p = 0.06). The stress of parturition and raising the first calf of a cow's life may be responsible for TL attenuation. Parity may be considered a stressor of cow longevity.

Impact of Brahman genetics on skin histology characteristics with implications for heat tolerance in cattle

Cattle lose heat predominantly through cutaneous evaporation at the skin-hair coat interface when experiencing heat stress. Sweating ability, sweat gland properties, and hair coat properties are a few of the many variables determining the efficacy of evaporative cooling. Sweating is a significant heat dissipation mechanism responsible for 85% of body heat loss when temperatures rise above 86⁰F. The purpose of this study was to characterize skin morphological parameters in Angus, Brahman, and their crossbred cattle. Skin samples were taken during the summer of 2017 and 2018 from a total of 319 heifers from six breed groups ranging from 100% Angus to 100% Brahman. Epidermis thickness decreased as the percentage of Brahman genetics increased where the 100% Angus group had a significantly thicker epidermis compared to the 100% Brahman animals. A more extended epidermis layer was identified in Brahman animals due to more pronounced undulations in this skin layer. Breed groups with 75% and 100% Brahman genes were similar and had the largest sweat gland area, indicative of superior resilience to heat stress, compared to breed groups with 50% or lower Brahman genetics. There was a significant linear breed group effect on sweat gland area indicating an increase of 862.0 µm² for every 25% increase in Brahman genetics. Sweat gland length increased as the Brahman percentage increased, while the sweat gland depth showed an opposite trend, decreasing from 100% Angus to 100% Brahman. The number of sebaceous glands was highest in 100% Brahman animals which had about 1.77 more sebaceous glands (p < 0.05) per 4.6 mm²area. Conversely, the sebaceous gland area was greatest in the 100% Angus group. This study identified significant differences in skin properties related to heat exchange ability between Brahman and Angus cattle. Equally important, these differences are also accompanied by significant levels of variation within each breed, which is indicative that selection for these skin traits would improve the heat exchange ability in beef cattle. Further, selecting beef cattle for these skin traits would lead to increased resilience to heat stress without disrupting production traits.

Inter-Individual Variation in DNA Methylation Patterns across Two Tissues and Leukocytes in Mature Brahman Cattle

Quantifying the natural inter-individual variation in DNA methylation patterns is important for identifying its contribution to phenotypic variation, but also for understanding how the environment affects variability, and for incorporation into statistical analyses. The inter-individual variation in DNA methylation patterns in female cattle and the effect that a prenatal stressor has on such variability have yet to be quantified. Thus, the objective of this study was to utilize methylation data from mature Brahman females to quantify the inter-individual variation in DNA methylation. Pregnant Brahman cows were transported for 2 h durations at days 60 ± 5; 80 ± 5; 100 ± 5; 120 ± 5; and 140 ± 5 of gestation. A non-transport group was maintained as a control. Leukocytes, amygdala, and anterior pituitary glands were harvested from eight cows born from the non-transport group (Control) and six from the transport group (PNS) at 5 years of age. The DNA harvested from the anterior pituitary contained the greatest variability in DNA methylation of cytosine-phosphate-guanine (mCpG) sites from both the PNS and Control groups, and the amygdala had the least. Numerous variable mCpG sites were associated with retrotransposable elements and highly repetitive regions of the genome. Some of the genomic features that had high variation in DNA methylation are involved in immune responses, signaling, responses to stimuli, and metabolic processes. The small overlap of highly variable CpG sites and features between tissues and leukocytes supports the role of variable DNA methylation in regulating tissue-specific gene expression. Many of the CpG sites that exhibited high variability in DNA methylation were common between the PNS and Control groups within a tissue, but there was little overlap in genomic features with high variability. The interaction between the prenatal environment and the genome could be responsible for the differences in location of the variable DNA methylation.

Comparisons of Hematological and Biochemical Profiles in Brahman and Yunling Cattle

Brahman cattle are tolerant to parasite challenges and heat stress. Yunling cattle are three-way hybrids that are half Brahman cattle, a fourth Murray Grey cattle and a fourth Yunnan Yellow cattle, with good beef performance. The hematological and biochemical parameters can reflect the physiology and metabolic conditions of cattle, and there are valuable indicators of production performance and adaptability that can be found by studying the cattle. To assess the health status and differences, we compared 55 hematological and biochemical parameters of 28 Brahman cattle and 65 Yunling cattle using an automatic biochemical analyzer. Our results showed that 27 hematological and biochemical indices of Brahman cattle were lower than those of Yunling cattle, whereas the other parameters were higher. There are 20 indices with significant differences that were detected between Brahman and Yunling cattle (with p ≤ 0.01 or 0.01 ≤ p ≤ 0.05, respectively), and no significant differences were found for other indices (p > 0.05). Based on these results, Yunling cattle may have a better physical condition than Brahman cattle, may be better at adapting to local environments, and can maintain a good production and reproduction performance. As a new breed that is half Brahman, the abilities of Yunling cattle, including adaptability, stress resistance and tolerance to crude feed, were better than Brahman cattle under the same management conditions.

Genetic variation in PLAG1 is associated with early fertility in Australian Brahman cattle

Variation in the genome region coding for PLAG1 has well-documented associations with skeletal growth and age at puberty in cattle. However, the influence of PLAG1 on other economically important traits such as cow stayability has not yet been explored. Here we investigate the effect of PLAG1 variation on early and later in life female fertility, as well as size and growth, in a well phenotyped Australian Brahman herd. Yearly pregnancy and productivity records were collected from 2,839 genotyped Brahman cows and used to generate fertility, growth, and weight phenotypes. A variant on chromosome 14 in PLAG1 (NC_037341.1:g.23338890G>T, rs109815800) was previously determined to be a putative causative mutation associated with variation in cattle stature. The imputed PLAG1 genotype at this variant was isolated for each animal and the effect of PLAG1 genotype on each trait was estimated using linear modelling. Regardless of how heifer fertility was measured, there was a significant (P < 0.05) and desirable relationship between the additive effects of PLAG1 genotype and successful heifer fertility. Heifers with two copies of the alternate allele (TT) conceived earlier and had higher pregnancy and calving rates. However, the effects of PLAG1 genotype on fertility began to diminish as cows aged and did not significantly influence stayability at later ages. While there was no effect of genotype on growth, PLAG1 had a negative effect on mature cow weight (P < 0.01), where females with two copies of the alternate allele (TT) were significantly smaller than those with either one or none. Selection emphasis on improved Brahman heifer fertility will likely increase the frequency of the T allele of rs109815800, which may also increase herd profitability and long-term sustainability through improved reproductive efficiency and reduced mature cow size.

Backfat Thickness Does Affect the Restoration of Ovarian Activity Postpartum in Different Breeding Programs in Zebu Cattle

The backfat thickness (BFT) was used to forecast the onset of ovarian activity and predict the calf growth. Eighty Brahman cows with their calves were allocated in two groups of 40 distributed in 4 months according to the month of calving, starting in March and finishing in June. One was synchronized and inseminated at fixed time following by natural mating (TAI+NM); whilst the other only by natural mating (NM). The programs started at 60 ± 5 days postpartum and ended 60 days later. From day 30 postpartum, serial ultrasound examinations and progesterone samples were used to monitor the onset of ovarian activity. The BFT in the rump area was measured by ultrasound from 30 days postpartum and every 15 days thereafter. The weight of the calves was recorded at birth and at weaning on 160 days. The adjusted effect of BFT on ovarian activity and the calves' development was assessed by binomial logistic regression at 30, 60, 75, and 120 days postpartum. The cycling cows averaged higher BFT irrespective of breeding program (P < 0.001). Also, slower changes in BFT were recorded during the follow-up at each time for all cows. However, the former had the higher BFT values from calving to the end of the study (P < 0.001). At 60, 75, and 120 days, the BFT measured, at the preceding time, was the only factor predicting the commencement of cyclicity (P < 0.001). The accumulative pregnancy through time was higher in TAI+NM (P = 0.003). Daily weight gain and weaning weight of the calves born in March was significantly heavier (P < 0.001) than peers born in April, May, or June. The most critical element to forecast the onset of ovarian activity is the monitoring of BFT around calving regardless of the breeding program. BFT to estimate the development of the calves until weaning was unpredictable.

Comparison of Gene Editing Versus Conventional Breeding to Introgress the POLLED Allele Into the Tropically Adapted Australian Beef Cattle Population

Dehorning is the process of physically removing horns to protect animals and humans from injury, but the process is costly, unpleasant, and faces increasing public scrutiny. Genetic selection for polled (hornless), which is genetically dominant to horned, is a long-term solution to eliminate the need for dehorning. However, due to the limited number of polled Australian Brahman bulls, the northern Australian beef cattle population remains predominantly horned. The potential to use gene editing to produce high-genetic-merit polled cattle was recently demonstrated. To further explore the concept, this study simulated introgression of the POLLED allele into a tropically adapted Australian beef cattle population via conventional breeding or gene editing (top 1% or 10% of seedstock bulls/year) for 3 polled mating schemes and compared results to baseline selection on genetic merit (Japan Ox selection index, $JapOx) alone, over the course of 20 years. The baseline scenario did not significantly decrease the 20-year HORNED allele frequency (80%), but resulted in one of the fastest rates of genetic gain ($8.00/year). Compared to the baseline, the conventional breeding scenarios where polled bulls were preferentially used for breeding, regardless of their genetic merit, significantly decreased the 20-year HORNED allele frequency (30%), but resulted in a significantly slower rate of genetic gain ($6.70/year, P ≤ 0.05). The mating scheme that required the exclusive use of homozygous polled bulls, resulted in the lowest 20-year HORNED allele frequency (8%), but this conventional breeding scenario resulted in the slowest rate of genetic gain ($5.50/year). The addition of gene editing the top 1% or 10% of seedstock bull calves/year to each conventional breeding scenario resulted in significantly faster rates of genetic gain (up to $8.10/year, P ≤ 0.05). Overall, our study demonstrates that, due to the limited number of polled Australian Brahman bulls, strong selection pressure on polled will be necessary to meaningfully increase the number of polled animals in this population. Moreover, these scenarios illustrate how gene editing could be a tool for accelerating the development of high-genetic-merit homozygous polled sires to mitigate the current trade-off of slower genetic gain associated with decreasing HORNED allele frequency in the Australian Brahman population.

The effect of Brahman genes on body temperature plasticity of heifers on pasture under heat stress

Bos taurus indicus cattle have the superior ability for the regulation of body temperature during heat stress due to a number of physiological and cellular level adaptive traits. The objectives of this study were to quantify the change in body temperature in heifers with various proportions of Brahman genes per unit increase in heat stress as measured by temperature-humidity index (THI) and to assess how different breed groups responded to varying intensity and duration of heat stress. A total of 299 two-yr-old heifers from six breed groups ranging from 100% Angus to 100% Brahman were evaluated under hot and humid conditions during 2017 and 2018 summer days. Two strategies were used to estimate the plasticity in body temperature of breed groups in response to environmental challenges: 1) a random regression mixed model was used to estimate reaction norm parameters for each breed group in response to a specified environmental heat stress and 2) a repeated measures mixed model was used to evaluate the response to different environmental heat loads. The reaction norm model estimated an intercept and slope measuring the change in body temperature per unit increase in THI environmental heat stress for different breed groups of animals and allowed the identification of genotypes which are robust, with low slope values indicative of animals that are able to maintain normal body temperature across a range of environments. The repeated measures mixed model showed that Brahman cattle have an advantage under moderate or high heat stress conditions but both Angus and Brahman breed groups are greatly affected when heat stress is severe. A critical factor appears to be the opportunity to cool down during the night hours more than the number of hours with extreme THI. With heat stress conditions predicted to intensify and expand into currently temperate zones, developing effective strategies to ensure sustainable beef production systems are imperative. Effective strategies will require the identification of the genes conferring the superior thermotolerance in Brahman cattle.

Genome-wide DNA methylation alteration in prenatally stressed Brahman heifer calves with the advancement of age

Possible phenotypic impairments associated with maternal stress during gestation in beef cattle may be explained by epigenetic effects. This study examined the impact of prenatal transportation stress on DNA methylation of lymphocytes of Brahman cows over the first 5 years of life. Methylation analysis through reduced representation bisulphite sequencing was conducted on DNA from lymphocytes from 28 paired samples from 6 prenatally stressed (PNS) and 8 control (Control) females obtained initially when they were 28 days of age and 5 years of age. There were 14,386 CpG (C = cytosine; p = phosphate; G = guanine) sites differentially methylated (P < 0.01) in 5-yr-old Control cows compared to their lymphocyte DNA at 28 days of age, this number was slightly decreased in 5-yr-old PNS with 13,378 CpG sites. Only 2,749 age-related differentially methylated CpG sites were seen within PNS females. There were 2,637 CpG sites differentially methylated (P < 0.01) in PNS cows relative to Controls at 5 years of age. There were differentially methylated genes in 5-yr-old cows that contributed similarly to altered gene pathways in both treatment groups. Canonical pathways altered in PNS compared to Control cows at 5 years of age were mostly related to development and growth, nervous system development and function, and immune response. Prenatal stress appeared to alter the epigenome in Brahman cows compared to Control at 5 years of age, which implies a persistent intervention in DNA methylation in lymphocytes, and may confer long-lasting effects on gene expression, and consequently relevant phenotypic changes.

Mitochondrial oxygen consumption in early postmortem permeabilized skeletal muscle fibers is influenced by cattle breed

Functional properties and integrity of skeletal muscle mitochondria (mt) during the early postmortem period may influence energy metabolism and pH decline, thereby impacting meat quality development. Angus typically produce more tender beef than Brahman, a Bos indicus breed known for heat tolerance. Thus, our objectives were to compare mt respiratory function in muscle collected early postmortem (1 h) from Angus and Brahman steers (n = 26); and to evaluate the effect of normal and elevated temperature on mt function ex vivo. We measured mt oxygen consumption rate (OCR) in fresh-permeabilized muscle fibers from Longissimus lumborum (LL) at 2 temperatures (38.5 and 40.0 °C) and determined citrate synthase (CS) activity and expression of several mt proteins. The main effects of breed, temperature, and their interaction were tested for mt respiration, and breed effect was tested for CS activity and protein expression. Breed, but not temperature (P > 0.40), influenced mt OCR (per tissue weight), with Brahman exhibiting greater complex I+II-mediated oxidative phosphorylation capacity (P = 0.05). Complex I- and complex II-mediated OCR also tended to be greater in Brahman (P = 0.07 and P = 0.09, respectively). Activity of CS was higher in LL from Brahman compared to Angus (P = 0.05). Expression of specific mt proteins did not differ between breeds, except for higher expression of adenosine triphosphate (ATP) synthase subunit 5 alpha in Brahman muscle (P = 0.04). Coupling control ratio differed between breeds (P = 0.05), revealing greater coupling between oxygen consumption and phosphorylation in Brahman. Our data demonstrate that both Angus and Brahman mt retained functional capacity and integrity 1-h postmortem; greater oxidative phosphorylation capacity and coupling in Brahman mt could be related to heat tolerance and impact early postmortem metabolism.

Genetic diversity of BoLA-DRB3 in South American Zebu cattle populations

Background

Bovine leukocyte antigens (BoLAs) are used extensively as markers of disease and immunological traits in cattle. However, until now, characterization of BoLA gene polymorphisms in Zebu breeds using high resolution typing methods has been poor. Here, we used a polymerase chain reaction sequence-based typing (PCR-SBT) method to sequence exon 2 of the BoLA class II DRB3 gene from 421 cattle (116 Bolivian Nellore, 110 Bolivian Gir, and 195 Peruvian Nellore-Brahman). Data from 1416 Taurine and Zebu samples were also included in the analysis.

Results

We identified 46 previously reported alleles and no novel variants. Of note, 1/3 of the alleles were detected only in Zebu cattle. Comparison of the degree of genetic variability at the population and sequence levels with genetic distance in the three above mentioned breeds and nine previously reported breeds revealed that Zebu breeds had a gene diversity score higher than 0.86, a nucleotide diversity score higher than 0.06, and a mean number of pairwise differences greater than 16, being similar to those estimated for other cattle breeds. A neutrality test revealed that only Nellore-Brahman cattle showed the even gene frequency distribution expected under a balanced selection scenario. The F_ST index and the exact G test showed significant differences across all cattle populations (F_ST = 0.057; p <  0.001). Neighbor-joining trees and principal component analysis identified two major clusters: one comprising mainly European Taurine breeds and a second comprising Zebu breeds. This is consistent with the historical and geographical origin of these breeds. Some of these differences may be explained by variation of amino acid motifs at antigen-binding sites.

Conclusions

The results presented herein show that the historical divergence between Taurine and Zebu cattle breeds is a result of origin, selection, and adaptation events, which would explain the observed differences in BoLA-DRB3 gene diversity between the two major bovine types. This allelic information will be important for investigating the relationship between the major histocompatibility complex and disease, and contribute to an ongoing effort to catalog bovine MHC allele frequencies according to breed and location.

Electronic supplementary material

The online version of this article (10.1186/s12863-018-0618-7) contains supplementary material, which is available to authorized users.

Population Structure and Genomic Breed Composition in an Angus-Brahman Crossbred Cattle Population

Crossbreeding is a common strategy used in tropical and subtropical regions to enhance beef production, and having accurate knowledge of breed composition is essential for the success of a crossbreeding program. Although pedigree records have been traditionally used to obtain the breed composition of crossbred cattle, the accuracy of pedigree-based breed composition can be reduced by inaccurate and/or incomplete records and Mendelian sampling. Breed composition estimation from genomic data has multiple advantages including higher accuracy without being affected by missing, incomplete, or inaccurate records and the ability to be used as independent authentication of breed in breed-labeled beef products. The present study was conducted with 676 Angus–Brahman crossbred cattle with genotype and pedigree information to evaluate the feasibility and accuracy of using genomic data to determine breed composition. We used genomic data in parametric and non-parametric methods to detect population structure due to differences in breed composition while accounting for the confounding effect of close familial relationships. By applying principal component analysis (PCA) and the maximum likelihood method of ADMIXTURE to genomic data, it was possible to successfully characterize population structure resulting from heterogeneous breed ancestry, while accounting for close familial relationships. PCA results offered additional insight into the different hierarchies of genetic variation structuring. The first principal component was strongly correlated with Angus–Brahman proportions, and the second represented variation within animals that have a relatively more extended Brangus lineage—indicating the presence of a distinct pattern of genetic variation in these cattle. Although there was strong agreement between breed proportions estimated from pedigree and genetic information, there were significant discrepancies between these two methods for certain animals. This was most likely due to inaccuracies in the pedigree-based estimation of breed composition, which supported the case for using genomic information to complement and/or replace pedigree information when estimating breed composition. Comparison with a supervised analysis where purebreds are used as the training set suggest that accurate predictions can be achieved even in the absence of purebred population information.

Objective Measures for the Assessment of Post-Operative Pain in Bos indicus Bull Calves Following Castration

Simple Summary

Surgical castration of cattle is a common husbandry procedure, and although this procedure is known to cause pain in cattle and other species, in some countries it is often performed without anaesthesia or analgesia. Society is increasingly aware of this animal welfare issue and it is creating pressure to drive research into animal welfare science with the aim of identifying practical and economical approaches to pain management in livestock. To effectively manage pain, a pain assessment must be performed. Pain assessment methods are often subjective and therefore influenced by the observer. Ideally, objective assessments that generate consistent and repeatable results between observers should be identified. Bos indicus bull calves were divided into four groups: no castration (NC, n = 6); castration with pre-operative local anaesthetic (CL n = 12); castration with pre-operative anti-inflammatory medication (CM, n = 12); and, castration without pain relief (C, n = 12). A range of objective assessments was performed: bodyweight measurements, activity, and rest levels, and four different compounds in the blood. The results of this study suggest that animals rest for longer periods after the pre-operative administration of anti-inflammatory medication. The other objective assessments measured in this study were not able to consistently differentiate between treatment groups. These findings emphasise the need for alternative quantifiable and objective indicators of pain in Bos indicus bull calves.

Abstract

The aim of the study was to assess pain in Bos indicus bull calves following surgical castration. Forty-two animals were randomised to four groups: no castration (NC, n = 6); castration with pre-operative lidocaine (CL, n = 12); castration with pre-operative meloxicam (CM, n = 12); and, castration alone (C, n = 12). Bodyweight was measured regularly and pedometers provided data on activity and rest from day −7 (7 days prior to surgery) to 13. Blood was collected for the measurement of serum amyloid A (SAA), haptoglobin, fibrinogen, and iron on days 0, 3 and 6. Bodyweight and pedometry data were analysed with a mixed effect model. The blood results were analysed with repeated measure one-way analysis of variance (ANOVA). There was no treatment effect on bodyweight or activity. The duration of rest was greatest in the CM group and lowest in the C group. There was a significant increase in the concentrations of SAA, haptoglobin, and fibrinogen in all of the groups from day 0 to 3. Iron concentrations were not different at the time points it was measured. The results of this study suggest that animals rest for longer periods after the pre-operative administration of meloxicam. The other objective assessments measured in this study were not able to consistently differentiate between treatment groups.

A field investigation of a modified intravaginal progesterone releasing device and oestradiol benzoate based ovulation synchronisation protocol designed for fixed-time artificial insemination of Brahman heifers

Pregnancy rates (PR) to fixed-time AI (FTAI) in Brahman heifers were compared after treatment with a traditional oestradiol-based protocol (OPO-8) or a modified protocol (OPO-6) where the duration of intravaginal progesterone releasing device (IPRD) was reduced from 8 to 6 days, and the interval from IPRD removal to oestradiol benzoate (ODB) was increased from 24 to 36 h. Rising 2 yo heifers on Farm A: (n = 238 and n = 215; two consecutive days AI); B (n = 271); and C (n = 393) were allocated to OPO-8 or OPO-6. An IPRD was inserted and 1mg ODB i.m. on Day 0 for OPO-8 heifers and Day 2 for OPO-6 heifers. On Day 8, the IPRD was removed and 500 μg cloprostenol i.m. At 24h, for OPO-8 heifers, and 36 h, for OPO-6 heifers, post IPRD removal all heifers received 1mg ODB i.m. FTAI was conducted at 54 and 72 h post IPRD removal for OPO-8 and OPO-6 heifers. At Farm A, OPO-6 heifers, AI on the second day, the PR was 52.4% to FTAI (P = 0.024) compared to 36.8% for OPO-8 heifers. However, no differences were found between OPO-8 and OPO-6 protocols at Farm A (first day of AI) (39.9 vs. 35.7%), or Farms B (26.2 vs. 35.4%) and C (43.2% vs. 40.3%). Presence of a corpus luteum at IPRD insertion affected PR to FTAI (43.9% vs. 28.8%; P < 0.001). This study has shown that the modified ovulation synchronisation protocol OPO-6 may be a viable alternative to the OPO-8 protocol for FTAI in B. indicus heifers.