Bacterial culture and susceptibility test results for clinical mastitis samples from Australia's subtropical dairy region.

This study aimed to identify the pathogens isolated from the milk of cows with clinical mastitis in the subtropical region of Australia and to determine the antimicrobial susceptibility of these bacteria. Thirty dairy herds in the subtropical dairy region were asked to submit milk samples for the first 5 cases of clinical mastitis each month for 12 mo. Samples underwent aerobic culture, and isolates were identified via MALDI-TOF mass spectrometry. Antimicrobial susceptibility was determined for Escherichia coli, Enterococcus spp., Streptococcus agalactiae, Streptococcus uberis, Streptococcus dysgalactiae, Staphylococcus aureus, and non-aureus staphylococci and mammaliicocci (NASM). Between March 2021 and July 2022, 1,230 milk samples were collected. A positive culture result was recorded for 812 (66%) of the milk samples; from these samples, 909 isolates were obtained, including 49 isolates where no identification was possible. The remaining samples were classified as having no growth (16.8%) or as being contaminated (17.2%). The most common isolates with a MALDI-TOF diagnosis (n = 909) were Strep. uberis (23.6%), followed by the NASM group (15.0%). Farms enrolled in the study were in 3 distinct locations within the subtropical dairy region: North Queensland, Southeast Queensland, and Northern New South Wales. Some variation in isolate prevalence occurred between these 3 locations. We found lower odds of a sample being positive for E. coli in North Queensland (odds ratio [OR]: 0.25; 95% confidence interval [CI]: 0.07-0.87) and higher odds in Southeast Queensland (OR: 4.01; 95% CI: 1.96-8.20) compared with the reference, Northern New South Wales. We further found higher odds of Strep. dysgalactiae in North Queensland (OR: 5.69; 95% CI: 1.85-17.54) and Southeast Queensland compared with Northern New South Wales (OR: 3.99; 95% CI: 1.73-9.22). Although some seasonal patterns were observed, season was not significant for any of the analyzed isolates. Farm-level differences in pathogen profiles were obvious. Overall, clinical mastitis pathogens had low levels of resistance to the antimicrobials tested. This research demonstrates that Strep. uberis and the NASM bacterial group are the most common pathogens causing clinical mastitis in the subtropical dairy region. It highlights the importance of understanding pathogenic causes of mastitis at the farm and regional level for targeted control and therapy.

Concurrently mapping quantitative trait loci associations from multiple subspecies within hybrid populations.

Many of the world's agriculturally important plant and animal populations consist of hybrids of subspecies. Cattle in tropical and sub-tropical regions for example, originate from two subspecies, Bos taurus indicus (Bos indicus) and Bos taurus taurus (Bos taurus). Methods to derive the underlying genetic architecture for these two subspecies are essential to develop accurate genomic predictions in these hybrid populations. We propose a novel method to achieve this. First, we use haplotypes to assign SNP alleles to ancestral subspecies of origin in a multi-breed and multi-subspecies population. Then we use a BayesR framework to allow SNP alleles originating from the different subspecies differing effects. Applying this method in a composite population of B. indicus and B. taurus hybrids, our results show that there are underlying genomic differences between the two subspecies, and these effects are not identified in multi-breed genomic evaluations that do not account for subspecies of origin effects. The method slightly improved the accuracy of genomic prediction. More significantly, by allocating SNP alleles to ancestral subspecies of origin, we were able to identify four SNP with high posterior probabilities of inclusion that have not been previously associated with cattle fertility and were close to genes associated with fertility in other species. These results show that haplotypes can be used to trace subspecies of origin through the genome of this hybrid population and, in conjunction with our novel Bayesian analysis, subspecies SNP allele allocation can be used to increase the accuracy of QTL association mapping in genetically diverse populations.

Bacterial culture and antimicrobial susceptibility results from bovine milk samples submitted to four veterinary diagnostic laboratories in Australia from 2015 to 2019.

A 5-year retrospective study was conducted to describe the mastitis-causing organisms isolated from bovine milk samples submitted to four veterinary diagnostic laboratories in Australia. The aim of this study was to identify temporal, geographical, and seasonal patterns of occurrence for the organisms and report the in vitro susceptibility of the most common mastitis-causing pathogens. In total, 22,102 milk samples were submitted between 2015 and 2019. The results were reported as positive growth for at least one significant organism (n = 11,407; 51.6%), no growth (n = 5,782; 26.2%), and mixed/contaminated growth (n = 4,913; 22.2%). Culture results for no growth, gram-negative bacteria, and eukaryotic organisms were combined for each region, and they were accounted for between 23 and 46% of submissions. These results represent a subset of mastitis cases for which the antibiotic treatment may not be warranted. A total of 11,907 isolates were cultured from 11,407 milk samples. The most common isolated organisms were Streptococcus uberis [41.3%; 95% confidence interval (CI): 40.4-42.1%] and Staphylococcus aureus (23.6%; 95% CI: 22.8-24.3%). For S. uberis and S. aureus, there was an association between a positive culture result and the dairy region. All regions except for the Sub-tropical Dairy region were more likely to culture S. uberis compared to the reference, Dairy NSW (P < 0.001). Similarly, for S. aureus, a positive culture result was more likely in all other dairy regions compared to Dairy NSW (P < 0.001). The LISA cluster analysis identified differences between High-High (hotspot) postcodes for S. aureus and S. uberis throughout all the analyzed dairy regions. These results highlight the need for further investigations into specific risk factors, such as environmental factors and herd-level predictors, which may have influenced the observed regional variations. Common mastitis-causing pathogens showed overall good susceptibility to a range of antimicrobials used in the treatment of mastitis. On-going surveillance of mastitis-causing pathogens and their antimicrobial susceptibilities will facilitate targeted mastitis control and treatment programs.

Investigation and Management of an Outbreak of Lead Intoxication in an Extensively Managed Beef Herd.

Fifteen hundred 12−15-month-old tropically adapted heifers inadvertently grazed a paddock which had a refuse dump in it containing burnt out vehicle batteries. The cattle grazed this paddock for approximately seven days. Subsequently these cattle were managed as two cohorts (cull and potential replacement breeding animals). Deaths commenced in the cull heifer group approximately 18 days after initial exposure to the refuse dump during relocation to a feedlot. Mortalities continued for 12 days, with other heifers showing clinical signs of marked central nervous system dysfunction requiring euthanasia. Necropsy of several clinically affected cattle plus blood sampling for lead analysis confirmed a diagnosis of lead intoxication. The crude mortality rate in the cull heifers was 6.6% (n = 685). Following confirmation of the diagnosis most of the potential replacement heifers (second cohort) were also relocated to the feedlot. The estimated crude mortality rate in this cohort was 5.8% (n = 815). All possible lead intoxication deaths occurred within 34 days of initial exposure, and apparently after day 16 at the feedlot no further heifers showed any clinical signs which could be attributed to lead intoxication. Longitudinal monitoring of blood lead concentrations was used to identify cattle suitable for slaughter. Overall, 70% of heifers initially blood sampled (n = 1408) had no detectable lead in their blood, however 16% had markedly elevated blood lead concentrations (> 0.7µmol/L) which persisted, and 2% had above the maximum normal threshold 1.5 years later. These latter cattle were subsequently euthanized, and necropsy revealed that visible pieces of lead were still present in the reticulum of several animals. At no time did any of these heifers with persistently high blood lead concentrations show clinical signs of lead intoxication.

Multi-breed genomic predictions and functional variants for fertility of tropical bulls.

Worldwide, most beef breeding herds are naturally mated. As such, the ability to identify and select fertile bulls is critically important for both productivity and genetic improvement. Here, we collected ten fertility-related phenotypes for 6,063 bulls from six tropically adapted breeds. Phenotypes were comprised of four bull conformation traits and six traits directly related to the quality of the bull's semen. We also generated high-density DNA genotypes for all the animals. In total, 680,758 single nucleotide polymorphism (SNP) genotypes were analyzed. The genomic correlation of the same trait observed in different breeds was positive for scrotal circumference and sheath score on most breed comparisons, but close to zero for the percentage of normal sperm, suggesting a divergent genetic background for this trait. We confirmed the importance of a breed being present in the reference population to the generation of accurate genomic estimated breeding values (GEBV) in an across-breed validation scenario. Average GEBV accuracies varied from 0.19 to 0.44 when the breed was not included in the reference population. The range improved to 0.28 to 0.59 when the breed was in the reference population. Variants associated with the gene HDAC4, six genes from the spermatogenesis-associated (SPATA) family of proteins, and 29 transcription factors were identified as candidate genes. Collectively these results enable very early in-life selection for bull fertility traits, supporting genetic improvement strategies currently taking place within tropical beef production systems. This study also improves our understanding of the molecular basis of male fertility in mammals.

Pregnancy rate per cycle is heritable and reduces with cycle in naturally mated tropically adapted beef cows.

Pregnancy in cattle is the outcome of the complex process of initiation of cycling, fertilization, maternal recognition of pregnancy and foeto-placental development. Though much is known about initiation of cycling and associated risk factors, there are virtually no data on pregnancy rate per cycle for naturally mated cattle, especially for extensively managed, tropically adapted genotypes, which this study aimed to determine. Tropical composite (Bos indicus and African Sanga crosses with Bos taurus) and Brahman cattle (n = 2,181) of known pedigree in four-year groups at four sites were mated annually for 84 days. Body condition, ovarian function, pregnancies, calving and lactation were monitored through six full reproductive cycles using 4-8 weekly ultrasound of the reproductive tract outside the calving period and daily monitoring during calving. From this, dates of commencement of cycling and conception in each year were estimated for each animal, enabling calculation of established pregnancy for consecutive 21-day periods while cycling and of pregnancies within four months of calving while lactating (P4M). Pregnancy per 21-day period (cycle) during mating for cycling animals averaged 63%, 71%, 41% and 28% in four consecutive cycles. Pregnant per cycle was 2%-11% higher in tropical composites than in Brahmans. The only other consistently significant risk to becoming pregnant was if cycling commenced later than three weeks before mating commenced. P4M averaged 62% and was lower for cows in sub-optimal body condition and in first-parity and later-calving cows. Pregnant per cycle was moderately heritable (~20%), while heritability was moderate to high (33%) for P4M. Selection for pregnant per cycle could be achieved indirectly by selection for P4M, a trait that is readily measured.

Comparison of Genetic Merit for Weight and Meat Traits between the Polled and Horned Cattle in Multiple Beef Breeds.

Breeding for polled animals is deemed the most practical solution to eradicate horns naturally and circumvent management costs and risks on health and welfare. However, there has been a historical reluctance by some farmers to select polled animals due to perceived lower productivity of their calves. This study has compared estimated breeding values (EBVs) between horned and polled animals (N = 2,466,785) for 12 production and carcass traits to assess historical (before 2000) and recent (2000-2018) genetic implications of poll breeding. Older generations of the polled animals in most breeds had significantly lower (Bonferroni-corrected p = 0.05) genetic merits for live (birth to maturity) and carcass weights, milk, meat quality, and fat content traits. Substantial gains of genetic potential were achieved during 2000 to 2018 in each breed, such that polled animals have significantly improved for the majority of traits studied. Generally, polled cohorts showed advantageous EBVs for live and carcass weights irrespective of the lower birth weights in some breeds. While Polled Brahman showed inferior production parameters, the poll genetics' effect size (d) and correlation (r) were very small on recent birth weight (d = -0.30, r = -0.08), 200 days (-0.19, -0.05), 400 days (-0.06, -0.02), 600 days (-0.05, -0.01), mature cow live weight (-0.08, -0.02), and carcass weight (-0.19, -0.05). In conclusion, although there is some evidence that historical selection for polled breeding animals may have reduced productivity, there is strong evidence that more recent selection for polled genotypes in the breeds studied has not resulted in any adverse effects on genetic merit.

Use of whole-genome sequence data and novel genomic selection strategies to improve selection for age at puberty in tropically-adapted beef heifers.

BACKGROUND: In tropically-adapted beef heifers, application of genomic prediction for age at puberty has been limited due to low prediction accuracies. Our aim was to investigate novel methods of pre-selecting whole-genome sequence (WGS) variants and alternative analysis methodologies; including genomic best linear unbiased prediction (GBLUP) with multiple genomic relationship matrices (MGRM) and Bayesian (BayesR) analyses, to determine if prediction accuracy for age at puberty can be improved. METHODS: Genotypes and phenotypes were obtained from two research herds. In total, 868 Brahman and 960 Tropical Composite heifers were recorded in the first population and 3695 Brahman, Santa Gertrudis and Droughtmaster heifers were recorded in the second population. Genotypes were imputed to 23 million whole-genome sequence variants. Eight strategies were used to pre-select variants from genome-wide association study (GWAS) results using conditional or joint (COJO) analyses. Pre-selected variants were included in three models, GBLUP with a single genomic relationship matrix (SGRM), GBLUP MGRM and BayesR. Five-way cross-validation was used to test the effect of marker panel density (6 K, 50 K and 800 K), analysis model, and inclusion of pre-selected WGS variants on prediction accuracy. RESULTS: In all tested scenarios, prediction accuracies for age at puberty were highest in BayesR analyses. The addition of pre-selected WGS variants had little effect on the accuracy of prediction when BayesR was used. The inclusion of WGS variants that were pre-selected using a meta-analysis with COJO analyses by chromosome, fitted in a MGRM model, had the highest prediction accuracies in the GBLUP analyses, regardless of marker density. When the low-density (6 K) panel was used, the prediction accuracy of GBLUP was equal (0.42) to that with the high-density panel when only six additional sequence variants (identified using meta-analysis COJO by chromosome) were included. CONCLUSIONS: While BayesR consistently outperforms other methods in terms of prediction accuracies, reasonable improvements in accuracy can be achieved when using GBLUP and low-density panels with the inclusion of a relatively small number of highly relevant WGS variants.

Optimized Genetic Testing for Polledness in Multiple Breeds of Cattle.

Many breeds of modern cattle are naturally horned, and for sound husbandry management reasons the calves frequently undergo procedures to physically remove the horns by disbudding or dehorning. These procedures are however a welfare concern. Selective breeding for polledness - absence of horns - has been effective in some cattle breeds but not in others (Bos indicus genotypes) due in part to the complex genetics of horn phenotype. To address this problem different approaches to genetic testing which provide accurate early-in-life prediction of horn phenotype have been evaluated, initially using microsatellites (MSAT) and more recently single nucleotide polymorphism (SNP). A direct gene test is not effective given the genetic heterogeneity and large-sized sequence variants associated with polledness in different breeds. The current study investigated 39,943 animals of multiple breeds to assess the accuracy of available poll testing assays. While the standard SNP-based test was an improvement on the earlier MSAT haplotyping method, 1999 (9.69%) out of 20,636 animals tested with this SNP-based assay did not predict a genotype, most commonly associated with the Indicus-influenced breeds. The current study has developed an optimized poll gene test that resolved the vast majority of these 1999 unresolved animals, while the predicted genotypes of those previously resolved remained unchanged. Hence the optimized poll test successfully predicted a genotype in 99.96% of samples assessed. We demonstrated that a robust set of 5 SNPs can effectively determine PC and PF alleles and eliminate the ambiguous and undetermined results of poll gene testing previously identified as an issue in cattle.

Multivariate genomic predictions for age at puberty in tropically adapted beef heifers.

Heifers that have an earlier age at puberty often have greater lifetime productivity. Age at puberty is moderately heritable so selection should effectively reduce the number of days to puberty, and improve heifer productivity and profitability as a result. However, recording age at puberty is intensive, requiring repeat ovarian scanning to determine age at first corpus luteum (AGECL). Genomic selection has been proposed as a strategy to select for earlier age at puberty; however, large reference populations of cows with AGECL records and genotypes would be required to generate accurate GEBV for this trait. Reproductive maturity score (RMS) is a proxy trait for age at puberty for implementation in northern Australia beef herds, where large scale recording of AGECL is not feasible. RMS assigns a score of 0 to 5 from a single ovarian scan to describe ovarian maturity at ~600 d. Here we use multivariate genomic prediction to evaluate the value of a large RMS data set to improve accuracy of GEBV for age at puberty (AGECL). There were 882 Brahman and 990 Tropical Composite heifers with AGECL phenotypes, and an independent set of 974 Brahman, 1,798 Santa Gertrudis, and 910 Droughtmaster heifers with RMS phenotypes. All animals had 728,785 real or imputed SNP genotypes. The correlation of AGECL and RMS (h2 = 0.23) was estimated as -0.83 using the genomic information. This result also demonstrates that using genomic information it is possible to estimate genetic correlations between traits collected on different animals in different herds, with minimal or unknown pedigree linkage between them. Inclusion of heifers with RMS in the multi-trait model improved the accuracy of genomic evaluations for AGECL. Accuracy of RMS GEBV generally did not improve by adding heifers with AGECL phenotypes into the reference population. These results suggest that RMS and AGECL may be used together in a multi-trait prediction model to increase the accuracy of prediction for age at puberty in tropically adapted beef cattle.

Towards multi-breed genomic evaluations for female fertility of tropical beef cattle.

Developing accurate genomic evaluations of fertility for tropical beef cattle must deal with at least two major challenges (i) recording cow fertility traits in extensive production systems on large numbers of cows and (ii) the genomic evaluations should work across the breeds, crossbreds, and composites used in tropical beef production. Here, we assess accuracy of genomic evaluations for a trait which can be collected on a large scale in extensive conditions, corpus luteum score (CLscore), which is 1 if ovarian scanning indicates a heifer has cycled by 600 d and 0 if not, in a multi-breed population. A total of 3,696 heifers, including 979 Brahmans, 914 Droughtmasters, and 1,803 Santa Gertrudis in seven herds across 3-yr cohorts with CLscores, were genotyped for 24,211 SNPs. Genotypes were imputed to 728,785 SNPs. GBLUP and BayesR were used to predict GEBV. Accuracy of GEBV was evaluated with two validation strategies. In the first strategy, the last year cohort of heifers from each herd was used for validation, such that every herd had heifers in both reference and validation populations. In the second validation strategy, each herd in turn was removed in its entirety from the reference population, and was used for validation. For both validation strategies, accuracy of GEBV for single breed and multi-breed reference populations was assessed. For the first validation strategy, accuracy of GEBV ranged from 0.2 for Brahmans to 0.4 for Droughtmasters. Increasing marker density from 24K SNPs to 728K SNPs resulted in a small increase in accuracy, and including multiple-breeds in the reference did not help improve accuracy. These results suggest that provided a herd has animals in the reference population, the accuracy of the GEBV is largely determined by within herd (linkage) information. The situation was very different when entire herds were predicted in the second validation. In this case accuracy of GEBV using only 24K SNPs and only a within breed reference was close to zero for all breeds. Accuracy increased substantially when 728K SNPs, BayesR, and a multi-breed reference were used, from 0.15 for Brahmans to 0.35 for Santa Gertrudis. Given the second validation strategy is more likely to reflect the situation for many herds in tropical beef production (no animals in the reference), genomic evaluations for fertility in tropical beef cattle should be based on high-density markers (728K SNPs) and should be multi-breed.

Scrotal circumference of Australian beef bulls.

Normal range for scrotal circumference in Australian beef bulls was established using more than 300,000 measurements of breed, management group, age, liveweight, and scrotal circumference. The data used were derived from Australian bull breeders and two large research projects in northern Australia. Most bulls were within 250 to 750 kg liveweight and 300 to 750 days of age. The differences between breeds and variances within breeds were higher when scrotal circumference was predicted from age rather than liveweight, because of variance in growth rates. The average standard deviation for predicted scrotal circumference from liveweight and age was 25 and 30 mm, respectively. Scrotal circumference by liveweight relationships have a similar pattern across all breeds, except in Waygu, with a 50 to 70 mm range in average scrotal circumference at liveweights between 250 and 750 kg. Temperate breed bulls tended to have higher scrotal circumference at the same liveweight than tropically adapted breeds. Five groupings of common beef breeds in Australian were identified, within which there were similar predictions of scrotal circumference from liveweight. It was concluded that liveweight and breed are required to identify whether scrotal circumference is within normal range for Australian beef bulls that experience a wide range of nutritional conditions.

The expression of prostaglandin endoperoxide synthase 2 messenger RNA and the proportion of smooth muscle and collagen in the sheep cervix during the estrous cycle.

The use of transcervical artificial insemination in sheep is limited because of the anatomy of the cervix, which restricts the passage of an inseminating pipette into the uterine lumen. There is a degree of natural cervical relaxation at estrus that enables greater penetration with an inseminating pipette. We hypothesize that this relaxation may be regulated by cervical prostaglandin synthesis and remodeling of the cervical extracellular matrix. The present study investigated the changes in prostaglandin endoperoxide synthase 2 (PTGS2) mRNA expression and the proportion of smooth muscle and collagen in the sheep cervix during the estrous cycle. Sheep cervices were collected at four stages of the estrous cycle: prior to the LH surge, during the LH surge, after the LH surge, and during the luteal phase. The expression of cervical PTGS2 mRNA was determined by in situ hybridization, and the proportion of smooth muscle and collagen in the cervix was investigated by Masson trichrome staining. The expression of PTGS2 mRNA in the sheep cervix was greatest prior to the LH surge, when estradiol concentrations were also greatest. The increase in PTGS2 mRNA expression was associated with an increase in the proportion of collagen in the sheep cervix. We propose that prior to the LH surge, estradiol may stimulate PTGS2 mRNA expression and hence prostaglandin E2 synthesis in the sheep cervix to regulate cervical relaxation, most likely through the rearrangement of collagen bundles within the cervical extracellular matrix.

The anatomy of the sheep cervix and its influence on the transcervical passage of an inseminating pipette into the uterine lumen.

The anatomy of the sheep cervix is highly variable between animals and may explain the differing success of transcervical AI between individuals. This study aims to quantify the variation in cervical morphology between ewes and investigate the relationship between cervical anatomy and cervical penetration. Two series of reproductive tracts were collected. Series A: 132 adult anoestrous ewes, and series B: 165 cycling adult ewes and ewe lambs which were identified as luteal or non-luteal based on the presence of a corpus luteum. The morphology of the cervical external os was classified as slit, papilla, duckbill, flap or rose. An inseminating pipette was inserted into the lumen and the depth of penetration recorded. The cervix was opened longitudinally, its length recorded, the number of cervical rings counted and the arrangement of those rings graded. The maximum depth of cervical penetration was affected by cervical grade (series A: P=0.021; series B: P=0.037) and the stage of the oestrous cycle (P=0.008). Grade 1 cervices were more penetrable than grade 2, with grade 3 the least penetrable and non-luteal cervices could be penetrated further than luteal cervices. The distribution of os types differed with age, with rose types more common in adult ewes, and papilla os types more common in ewe lambs. These results indicate that the depth of cervical penetration is affected by the anatomy of the cervical lumen. Cervices with a less convoluted lumen (grade 1) were more penetrable. Non-luteal cervices are likely to have higher oestradiol concentrations than luteal, stimulating cervical relaxation and enabling deeper penetration. The difference in os types with age may be contributable to a morphological alteration at parturition.