Search

ABSTRACT: At the individual cow level, suboptimum fertility, mastitis, negative energy balance, and ketosis are major issues in dairy farming. These problems are widespread on dairy farms and have an important economic impact. The objectives of this study were (1) to assess the potential of milk mid-infrared (MIR) spectra to predict key biomarkers of energy deficit (citrate, isocitrate, glucose-6 phosphate [glucose-6P], free glucose), ketosis (β-hydroxybutyrate [BHB] and acetone), mastitis (N-acetyl-β-d-glucosaminidase activity [NAGase] and lactate dehydrogenase), and fertility (progesterone); (2) to test alternative methodologies to partial least squares (PLS) regression to better account for the specific asymmetric distribution of the biomarkers; and (3) to create robust models by merging large datasets from 5 international or national projects. Benefiting from this international collaboration, the dataset comprised a total of 9,143 milk samples from 3,758 cows located in 589 herds across 10 countries and represented 7 breeds. The samples were analyzed by reference chemistry for biomarker contents, whereas the MIR analyses were performed on 30 instruments from different models and brands, with spectra harmonized into a common format. Four quantitative methodologies were evaluated to address the strongly skewed distribution of some biomarkers. Partial least squares regression was used as the reference basis, and compared with a random modification of distribution associated with PLS (random-downsampling-PLS), an optimized modification of distribution associated with PLS (KennardStone-downsampling-PLS), and support vector machine (SVM). When the ability of MIR to predict biomarkers was too low for quantification, different qualitative methodologies were tested to discriminate low versus high values of biomarkers. For each biomarker, 20% of the herds were randomly removed within all countries to be used as the validation dataset. The remaining 80% of herds were used as the calibration dataset. In calibration, the 3 alternative methodologies outperform the PLS performances for the majority of biomarkers. However, in the external herd validation, PLS provided the best results for isocitrate, glucose-6P, free glucose, and lactate dehydrogenase (coefficient of determination in external herd validation [R2v] = 0.48, 0.58, 0.28, and 0.24, respectively). For other molecules, PLS-random-downsampling and PLS-KennardStone-downsampling outperformed PLS in the majority of cases, but the best results were provided by SVM for citrate, BHB, acetone, NAGase, and progesterone (R2v = 0.94, 0.58, 0.76, 0.68, and 0.15, respectively). Hence, PLS and SVM based on the entire dataset provided the best results for normal and skewed distributions, respectively. Complementary to the quantitative methods, the qualitative discriminant models enabled the discrimination of high and low values for BHB, acetone, and NAGase with a global accuracy around 90%, and glucose-6P with an accuracy of 83%. In conclusion, MIR spectra of milk can enable quantitative screening of citrate as a biomarker of energy deficit and discrimination of low and high values of BHB, acetone, and NAGase, as biomarkers of ketosis and mastitis. Finally, progesterone could not be predicted with sufficient accuracy from milk MIR spectra to be further considered. Consequently, MIR spectrometry can bring valuable information regarding the occurrence of energy deficit, ketosis, and mastitis in dairy cows, which in turn have major influences on their fertility and survival.

Crossbreeding of indigenous tropical and improved western dairy cattle breeds as tool to improve dairy cattle performance on smallholder farms has been widely advocated, criticised and yet applied. The government of Ethiopia supported this technology for decades but adoption rate is low. Constraints are documented but there is little information about farm level introduction and development of crossbreeding. A total 122 smallholders with mixed crop livestock farms and at least 8 years of successful crossbreeding were interviewed using a pre-tested questionnaire in two contexts in Amhara Regional state in north-western Ethiopia. Crossbreeding initiator was either uncoordinated government extension or a coordinated development project, also implemented with governmental support. Qualitative and quantitative data on farmers’ motivations, crossbreeding introduction, initiator support, breeding adaptation and impacts at farm level were analysed. Results show that even though motives vary between contexts the underlying reason to introduce crossbreeding was economic profit. To be able to introduce crossbreeding support of initiators (e.g. extension) and other farmers was essential. The crossbreeding introduction context had some influence. Governmental actors were the main source of support and supplier of exotic genetics but the farmer network acted as safety net filling gaps of government support. Breeding strategies focused on performance increase. A lack of basic understanding of crossbreeding has been identified. A surprising, probably biased, result was general satisfaction with initiator support and with breeding services. It was challenged by the high proportion of farmers unable to follow a breeding strategy due to insufficient bull and/or semen supply. Crossbreeding changed the smallholder production system to a high input – high output system. Except for crossbred adaptation problems, challenges were ranked context specific and influenced by the initiator. Farmers perceived crossbreeding as success and recommended it. We conclude that farmers can realize income increase with crossbreeding. The complexity of this technology, high initial investment and the need for support services and external production inputs are probable reasons why crossbreeding uptake is low. Improving the availability of semen and/or bulls must be the top priority for breeding service providers to enable farmers to follow a breeding strategy and reach a suitable and sustainable herd performance. Access to investment capital, input supply, strong technical support and market linkages are crucial for successful crossbreeding.

Breeding values for length of productive life in Slovak Pinzgau cattle were estimated using survival analysis. As the results were corrected for milk production, the final breeding values represented the ability of cows to avoid culling from reasons other than milk production. In addition to the relative milk yield, the risk of culling was also studied in connection with the herd and year of calving, parity and stage of lactation, herd size change and age at the first calving. Among the fixed effects, the low milk production, high age at the first calving, and decreasing herd size were associated with increased risk of culling. The risk was non-linear for parity × stage of lactation classes, decreasing within the first parity and increasing during later parities. Two genetic random effects were considered in separate models: the sire of the cow and the animal itself, both with the corresponding pedigree records up to the third generation. The genetic effects were estimated in separate runs, but the rest of the model remained unchanged. Heritability was h2 = 0.08 for sire model, confirming the results of an earlier study in the Slovak Pinzgau population. The computational feasibility of the animal model for estimation of breeding values for cows was confirmed in this study. The new breeding values could be computed for each cow, accounting for all relationships within the population. The estimated heritability for the animal model was h2 = 0.11, which should be used later on for animal breeding purposes.

Identification of breeding objective traits pertinent to specific production environments with the involvement of target beneficiaries is crucial to the success of a breed improvement program. A choice experiment was conducted in four locations representing different production systems and agro-ecologies that are habitat to four indigenous sheep breeds (Afar, Bonga, Horro and Menz) of Ethiopia with the objective of identifying farmers’/pastoralists’ preferences for sheep breeding traits. Following a synthesis of secondary information and diagnostic surveys, two communities per location consisting of 60 households each having at least four breeding ewes were identified. Producers’ priority attributes used in the choice sets were identified through in-depth production system studies conducted from December 2007 to March 2008. On the basis of prior information, four to seven attributes were used to design choice sets with different profiles in order to capture results that mimic real life of the different communities. The attributes and levels chosen for the sheep profile were as follows: body size (large/small), coat color (brown/white/black), tail type (good/bad) for both rams and ewes; horn (polled/horned) and libido (active/poor) for rams; and lambing interval (three lambings in 2 years/two lambings in 2 years time), mothering ability (good mother/bad mother), twinning rate (twin bearer/single bearer) and milk yield (two cups per milking/one cup per milking) for ewes. A fractional factorial design was implemented to construct the alternatives included in the choice sets. The design resulted in a randomized selection of 48 sheep profiles (24 sets) for both sexes, which were grouped into four blocks with six choice sets each. An individual respondent was presented with one of the four blocks to make his/her choices. Results indicate that producers’ trait preferences were heterogeneous except for body size in rams and mothering ability in ewes where nearly homogeneous preferences were investigated. In the pastoral production system, attention was given to coat color of both breeding rams and ewes, favoring brown and white colors over black. Ram libido influenced producers’ decisions in Bonga, Horro and Menz areas. The influence of milk yield and twinning on respondents’ decision making was high in Afar and Horro, respectively. Breeders in all areas attempt to combine production and reproduction traits as well as they can in order to maximize benefits from their sheep. The elicited measurable objective traits were used to design alternative community-based sheep breeding plans for the four indigenous sheep breeds in their production environments that have been implemented since.

The objective of this study was to compare linear models and survival analysis for genetic evaluation of ovulatory disorders, which included veterinary treatments of silent heat/anestrus and cystic ovaries. Data of 23 450 daughters of 274 Austrian Fleckvieh sires were analyzed. For linear model analyses, ovulatory disorders were defined as a binary response (presence or absence) in the time periods from calving to 150 days after calving and from calving to 300 days after calving. For survival analysis, ovulatory disorders were defined either as the number of days from calving to the day of the first treatment for an ovulatory disorder (uncensored record) or from calving to the day of culling, or the last day of the period under investigation (until 150 or 300 days after calving; censored record). Estimates of heritability were very similar (0.016 to 0.020) across methods and periods. Correlations between sire estimated breeding value from linear model and survival analysis were 0.98, whereas correlations between different time periods were somewhat lower (0.95 and 0.96). The results showed that the length of time period had a larger effect on genetic evaluation than methodology.

Survival Kit V3.12 was used to analyse the length of productive life of cattle in theCzechRepublic. The data set consisted of 230 028 registeredHolstein cows. The model included the time-dependent effects parity × stage of lactation interaction, herd × year × season interaction, class of milk production within herd and year, breed within years and the time-independent effect of age at first calving and the random effect of sire. The highest risk of culling was found for cows at the beginning and at the end of the first lactation and at the end of any other lactation. The risk of culling decreased with parity. The risk of culling of cows assigned to the lowest milk production class was five times higher than that of cows assigned to the average milk production class. Risk of culling diminished with a decreasing percentage ofHolstein breed. Cows younger at first calving showed a lower risk of culling. Breeding values for sires expressed as a risk ratio of their daughters were between 0.7 and 1.45. Estimated heritability of functional longevity was 0.025 on the log scale and 0.041 on the original scale.

The analysis was conducted of the national database of station tests carried out between May 1996- February 2001, using the Hungarian Large White breed. Days of test, total amount of consumed feed and valuable cuts were taken into the analysis. Using the method of cross validation, small subsets of the data were excluded and then predicted using the remaining part of the data treating herd × year effects either as fixed or as random. The size of the data excluded was 50 or 10 records at a time and the process was repeated 100 or 500 times, respectively. Mean squared error, bias and correlation between the excluded and predicted observations were calculated for all the excluded subsets. There was no significant difference between the fixed and random models but in the case of valuable cuts the random models showed a lower mean squared error and higher correlation between the excluded and predicted observations than the fixed models.

Numerous studies have found a large number of QTL associated with productive and functional traits in the cattle genome. Several Countries have already established research programs aiming at identification and exploitation of QTLs in dairy cattle (Bovenhuis and Schrooten 2002). However, the published results are not directly exploitable for MAS, as commercial interests limit the information in the literature, and because marker-QTL phase is specific to each family.

Milking speed is an important trait influencing udder health of dairy cows as well as labor efficiency. Yet, it has received little attention in genomic association studies. The main objective of this study was to determine regions and genes on the genome with a potential effect on milking speed in Fleckvieh (dual purpose Simmental) cattle. Genome-wide association studies were conducted using de-regressed breeding values of bulls as phenotypes. Six SNP on 4 autosomes were significantly associated with milking speed for additive effects. Significant regions on BTA4 and BTA19 correspond with findings for other dairy cattle breeds. Based on the observation of Fleckvieh breed managers, variation of milking speed in batches of daughters of some bulls is much higher than in daughter groups of other bulls. This difference in within family variation may be caused by transmission of alternative alleles of bulls being heterozygous for a gene affecting milking speed. To check on this, we considered standard deviation of yield deviations in milking speed of half-sib daughters as a new trait and performed GWAS for dominance effects. One signal on BTA5 passed the genome wide Bonferroni threshold that corresponded to the significant signal from standard GWAS on de-regressed breeding values. The key conclusion of this study is that several strong genomic signals were found for milking speed in Fleckvieh cattle and that the strongest of them are supported by similar findings in Brown Swiss and Holstein Friesian cattle. Milking speed is a complex trait whose sub-processes have not yet been elucidated in detail. Hence, it remains a challenge to link the associated regions on the genome with causal genes and their functions., (The Authors. Published by Elsevier Inc. on behalf of the American Dairy Science Association®. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).)

While dairy goat production, characterized by traditional production on small farms, is an important source of income in the Czech Republic and Slovakia, locally adapted breeds have not been fully consolidated over the last 100 years due to large fluctuations in population size and inconsistent breeding programs that allowed for different crossbreeding strategies. Our main objective in this study was therefore to assess the conservation status of 4 Czech (Alpine Goat, White Shorthair, Brown Shorthair and Czech Landrace) and one Slovak (Slovak White Shorthair) local goat breeds, to analyze their population structure and admixture, and to estimate their relatedness to several neighboring breeds. Our analyses included 142 goats belonging to 5 local breeds genotyped with the Illumina 50K BeadChip and 618 previously genotyped animals representing 15 goat breeds from Austria and Switzerland (all analyses based on 46,862 autosomal SNPs and 760 animals). In general, the conservation status of the Czech and Slovak local goat breeds was satisfactory, with the exception of the Brown Shorthair goat, as the analyzed parameters (heterozygosity, haplotype richness, ROH-based inbreeding and effective population size) were mostly above the median of 20 breeds. However, for all 5 Czech and Slovakian breeds, an examination of historical effective population size indicated a substantial decline about 8 to 22 generations ago. In addition, our study revealed that the Czech and Slovakian breeds are not fully consolidated; for instance, White Shorthair and Brown Shorthair were not clearly distinguishable. Considerable admixture, especially in Czech Landrace (effective number of parental clusters equal to 4.2), and low but numerous migration rates from other Austrian and Swiss breeds were found. These results provide valuable insights for future breeding programs and genetic diversity management of local Czech and Slovak goat breeds., (The Authors. Published by Elsevier Inc. on behalf of the American Dairy Science Association®. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).)

The pig industry is usually considered an intensive livestock industry, mainly supported by hybrid breeding between commercial pig breeds. However, people's pursuit of a more natural environment and higher meat quality has led to an increasing demand for eco-friendly and diverse pig feeding systems. Therefore, the importance of rearing and conserving local pig breeds is increasing. The Livni pig is a local breed with good adaptability to the environmental and fodder conditions in central Russia. In this study, we aimed to analyze the genetic diversity and population structure of Livni pigs using whole-genome single nucleotide polymorphism (SNP) data. We utilized the Porcine GGP HD BeadChip on genotype samples from old ( n = 32, 2004) and modern ( n = 32, 2019) populations of Livni pigs. For the museum samples of Livni pigs ( n = 3), we extracted DNA from their teeth, performed genomic sequencing, and obtained SNP genotypes from the whole-genome sequences. SNP genotypes of Landrace ( n = 32) and Large White ( n = 32) pigs were included for comparative analysis. We observed that the allelic richness of Livni pigs was higher than those of Landrace and Large White pigs ( A R = 1.775-1.798 vs. 1.703 and 1.668, respectively). The effective population size estimates ( N E 5 = 108 for Livni pigs, N E 5 = 59 for Landrace and Large White pigs) confirmed their genetic diversity tendency. This was further supported by the length and number of runs of homozygosity, as well as the genomic inbreeding coefficient (almost twofold lower in Livni pigs compared to Landrace and Large White pigs). These findings suggest that the Livni pig population exhibits higher genetic diversity and experiences lower selection pressure compared to commercial pig populations. Furthermore, both principal component and network tree analyses demonstrated a clear differentiation between Livni pigs and transboundary commercial pigs. The TreeMix results indicated gene flow from Landrace ancestors to Livni pigs (2019) and from Large White ancestors to Livni pigs (2004), which was consistent with their respective historical breeding backgrounds. The comparative analysis of museum, old, and modern Livni pigs indicated that the modern Livni pig populations have preserved their historical genomic components, suggesting their potential suitability for future design selection programs.

Genetic improvement of udder health in dairy cows is of high relevance as mastitis is one of the most prevalent diseases. Since it is known that the heritability of mastitis is low and direct data on mastitis cases are often not available in large numbers, auxiliary traits, such as somatic cell count (SCC) are used for the genetic evaluation of udder health. In previous studies, models to predict clinical mastitis based on mid-infrared (MIR) spectral data and a somatic cell count-derived score (SCS) were developed. Those models can provide a probability of mastitis for each cow at every test-day, which is potentially useful as an additional auxiliary trait for the genetic evaluation of udder health. Furthermore, MIR spectral data were used to estimate contents of lactoferrin, a glycoprotein positively associated with immune response. The present study aimed to estimate heritabilities (h 2 ) and genetic correlations (r a ) for clinical mastitis diagnosis (CM), SCS, MIR-predicted mastitis probability (MIRprob), MIR + SCS-predicted mastitis probability (MIRSCSprob) and lactoferrin estimates (LF). Data for this study were collected within the routine milk recording and health monitoring system of Austria from 2014 to 2021 and included records of approximately 54,000 Fleckvieh cows. Analyses were performed in two datasets, including test-day records from 5 to 150 or 5 to 305 days in milk. Prediction models were applied to obtain MIR- and SCS-based phenotypes (MIRprob, MIRSCSprob, LF). To estimate heritabilities and genetic correlations bivariate linear animal models were applied for all traits. A lactation model was used for CM, defined as a binary trait, and a test-day model for all other continuous traits. In addition to the random animal genetic effect, the fixed effects year-season of calving and parity-age at calving and the random permanent environmental effect were considered in all models. For CM the random herd-year effect, for continuous traits the random herd-test day effect and the covariate days in milk (linear and quadratic) were additionally fitted. The obtained genetic parameters were similar in both datasets. The heritability found for CM was expectedly low (h 2  = 0.02). For SCS and MIRSCSprob, heritability estimates ranged from 0.23 to 0.25, and for MIRprob and LF from 0.15 to 0.17. CM was highly correlated with SCS and MIRSCSprob (r a  = 0.85 to 0.88). Genetic correlations of CM were moderate with MIRprob (r a  = 0.26 and 0.37) during 150 and 305 days in milk, respectively and low with LF (h 2  = 0.10 and 0.11). However, basic selection index calculations indicate that the added value of the new MIR-predicted phenotypes is limited for genetic evaluation of udder health., (© 2024 The Authors. Journal of Animal Breeding and Genetics published by John Wiley & Sons Ltd.)

At the individual cow level, suboptimum fertility, mastitis, negative energy balance, and ketosis are major issues in dairy farming. These problems are widespread on dairy farms and have an important economic impact. The objectives of this study were (1) to assess the potential of milk mid-infrared (MIR) spectra to predict key biomarkers of energy deficit (citrate, isocitrate, glucose-6 phosphate [glucose-6P], free glucose), ketosis (β-hydroxybutyrate [BHB] and acetone), mastitis (N-acetyl-β-d-glucosaminidase activity [NAGase] and lactate dehydrogenase), and fertility (progesterone); (2) to test alternative methodologies to partial least squares (PLS) regression to better account for the specific asymmetric distribution of the biomarkers; and (3) to create robust models by merging large datasets from 5 international or national projects. Benefiting from this international collaboration, the dataset comprised a total of 9,143 milk samples from 3,758 cows located in 589 herds across 10 countries and represented 7 breeds. The samples were analyzed by reference chemistry for biomarker contents, whereas the MIR analyses were performed on 30 instruments from different models and brands, with spectra harmonized into a common format. Four quantitative methodologies were evaluated to address the strongly skewed distribution of some biomarkers. Partial least squares regression was used as the reference basis, and compared with a random modification of distribution associated with PLS (random-downsampling-PLS), an optimized modification of distribution associated with PLS (KennardStone-downsampling-PLS), and support vector machine (SVM). When the ability of MIR to predict biomarkers was too low for quantification, different qualitative methodologies were tested to discriminate low versus high values of biomarkers. For each biomarker, 20% of the herds were randomly removed within all countries to be used as the validation dataset. The remaining 80% of herds were used as the calibration dataset. In calibration, the 3 alternative methodologies outperform the PLS performances for the majority of biomarkers. However, in the external herd validation, PLS provided the best results for isocitrate, glucose-6P, free glucose, and lactate dehydrogenase (coefficient of determination in external herd validation [R 2 v] = 0.48, 0.58, 0.28, and 0.24, respectively). For other molecules, PLS-random-downsampling and PLS-KennardStone-downsampling outperformed PLS in the majority of cases, but the best results were provided by SVM for citrate, BHB, acetone, NAGase, and progesterone (R 2 v = 0.94, 0.58, 0.76, 0.68, and 0.15, respectively). Hence, PLS and SVM based on the entire dataset provided the best results for normal and skewed distributions, respectively. Complementary to the quantitative methods, the qualitative discriminant models enabled the discrimination of high and low values for BHB, acetone, and NAGase with a global accuracy around 90%, and glucose-6P with an accuracy of 83%. In conclusion, MIR spectra of milk can enable quantitative screening of citrate as a biomarker of energy deficit and discrimination of low and high values of BHB, acetone, and NAGase, as biomarkers of ketosis and mastitis. Finally, progesterone could not be predicted with sufficient accuracy from milk MIR spectra to be further considered. Consequently, MIR spectrometry can bring valuable information regarding the occurrence of energy deficit, ketosis, and mastitis in dairy cows, which in turn have major influences on their fertility and survival., (The Authors. Published by Elsevier Inc. and Fass Inc. on behalf of the American Dairy Science Association®. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).)

The search for SNPs and candidate genes that determine the manifestation of major selected traits is one crucial objective for genomic selection aimed at increasing poultry production efficiency. Here, we report a genome-wide association study (GWAS) for traits characterizing meat performance in the domestic quail. A total of 146 males from an F 2 reference population resulting from crossing a fast (Japanese) and a slow (Texas White) growing breed were examined. Using the genotyping-by-sequencing technique, genomic data were obtained for 115,743 SNPs (92,618 SNPs after quality control) that were employed in this GWAS. The results identified significant SNPs associated with the following traits at 8 weeks of age: body weight (nine SNPs), daily body weight gain (eight SNPs), dressed weight (33 SNPs), and weights of breast (18 SNPs), thigh (eight SNPs), and drumstick (three SNPs). Also, 12 SNPs and five candidate genes ( GNAL , DNAJC6 , LEPR , SPAG9 , and SLC27A4 ) shared associations with three or more traits. These findings are consistent with the understanding of the genetic complexity of body weight-related traits in quail. The identified SNPs and genes can be used in effective quail breeding as molecular genetic markers for growth and meat characteristics for the purpose of genetic improvement.

The community-based breeding program (CBBP) is an innovative approach recommended for genetic improvement and sustainable use of animal genetic resources in extensive farming systems. Successful implementation of this approach requires an understanding of the characteristics of production systems, breeding objectives, and farmers' trait preference. This study aimed to identify the selection criteria of goat farmers in rural areas of Burkina Faso and their potential implications in establishing CBBP. Following focus group discussions, a well-structured questionnaire was designed and administered to 372 randomly selected goat farmers in two different agro-ecological zones. A list of traits obtained during focus group discussions was provided to farmers individually, and they were asked to rank the ones they preferentially use to select breeding animals. Statistical tests were conducted to compare data between the two agro-ecological zones. The results showed that the average goat flock per household was higher (P < 0.05) in the Sudanian (15.68 ± 13.76), compared to the Sudano-Sahelian area (12.93 ± 13.3). Adult females were the dominant age-sex group in both areas. Reasons for culling, keeping breeding bucks, and castration practice were significantly different (P < 0.05) among agro-ecological zones. The most important common criterion for selection in the two zones was body size, coat color, and growth rate for the bucks and does, while fertility (0.06) parameters including twining ability (0.18), kidding frequency (0.11), and mothering ability (0.15) were furthermore considered for breeding does selection. These findings provide valuable insights for developing CBBPs tailored to goat production in the study areas., (© 2024. The Author(s), under exclusive licence to Springer Nature B.V.)

After the domestication of goats around 10,000 years before the present (BP), humans transported goats far beyond the range of their wild ancestor, the bezoar goat. This brought domestic goats into contact with many wild goat species such as ibex and markhor, enabling introgression between domestic and wild goats. To investigate this, while shedding light on the taxonomic status of wild and domestic goats, we analysed genome-wide SNP data of 613 specimens from 14 taxonomic units, including Capra hircus, C. pyrenaica, C. ibex (from Switzerland, Austria, Germany and Slovenia), C. aegagrus aegagrus, C. a. cretica, C. h. dorcas, C. caucasica caucasica, C. c. severtzovi, C. c. cylindricornis, C. falconeri, C. sibirica sibirica, C. s. alaiana and C. nubiana, as well as Oreamnos americanus (mountain goat) as an outgroup. To trace gene flow between domestic and wild goats, we integrated genotype data of local goat breeds from the Alps as well as from countries such as Spain, Greece, Türkiye, Egypt, Sudan, Iran, Russia (Caucasus and Altai) and Pakistan. Our phylogenetic analyses displayed a clear separation between bezoar-type and ibex-type clades with wild goats from the Greek islands of Crete and Youra clustered within domestic goats, confirming their feral origin. Our analyses also revealed gene flow between the lineages of Caucasian tur and domestic goats that most likely occurred before or during early domestication. Within the clade of domestic goats, analyses inferred gene flow between African and Iberian goats. The detected events of introgression were consistent with previous reports and offered interesting insights into the historical relationships among domestic and wild goats., (© 2023 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

The goal of our study was to identify signatures of selection in the Turopolje pigs and other commercial pig breeds. We conducted a comprehensive analysis of five datasets, including one local pig breed (Turopolje) and four commercial pig breeds (Large White, Landrace, Pietrain, and Duroc), using strict quality control measures. Our final dataset consisted of 485 individuals and 54,075 single nucleotide polymorphisms (SNPs). To detect selection signatures within these pig breeds, we utilized the XP-EHH and XP-nSL methodologies, which allowed us to identify candidate genes that have been subject to positive selection. Our analysis consistently highlighted the PTBP2 and DPYD genes as commonly targeted by selection in the Turopolje breed. DPYD is associated with muscular development in pigs and other species and PTBP2 emerges as one of the potential genes linked to seminal characteristics. Furthermore, in the Large White breed, a number of genes were detected with the two methods, such as ATP1A1 , CASQ2 , CD2 , IGSF3 , MAB21L3 , NHLH2 , SLC22A15 , VANGL1 . In the Duroc breed, a different set of genes was detected, such as ARSB , BHMT , BHMT2 , DMGDH , JMY . The function of these genes was related to body weight, production efficiency and meat quality, average daily gain, and other similar traits. Overall, our results have identified a number of genomic regions that are under selective pressure between local and commercial pig breeds. This information can help to improve our understanding of the mechanisms underlying pig breeding, and ultimately contribute to the development of more efficient and sustainable pig production practices. Our study highlights the power of using multiple genomic methodologies to detect genetic signatures of selection, and provides important insights into the genetic diversity of pig breeds.

Traces of long-term artificial selection can be detected in genomes of domesticated birds via whole-genome screening using single-nucleotide polymorphism (SNP) markers. This study thus examined putative genomic regions under selection that are relevant to the development history, divergence and phylogeny among Japanese quails of various breeds and utility types. We sampled 99 birds from eight breeds (11% of the global gene pool) of egg (Japanese, English White, English Black, Tuxedo and Manchurian Golden), meat (Texas White and Pharaoh) and dual-purpose (Estonian) types. The genotyping-by-sequencing analysis was performed for the first time in domestic quails, providing 62,935 SNPs. Using principal component analysis, Neighbor-Net and Admixture algorithms, the studied breeds were characterized according to their genomic architecture, ancestry and direction of selective breeding. Japanese and Pharaoh breeds had the smallest number and length of homozygous segments indicating a lower selective pressure. Tuxedo and Texas White breeds showed the highest values of these indicators and genomic inbreeding suggesting a greater homozygosity. We revealed evidence for the integration of genomic and performance data, and our findings are applicable for elucidating the history of creation and genomic variability in quail breeds that, in turn, will be useful for future breeding improvement strategies., Competing Interests: The authors declare no conflicts of interest.

In most cases, inbreeding is expected to have unfavourable effects on traits in livestock. The consequences of inbreeding depression could be substantial, primarily in reproductive and sperm quality traits, and thus lead to decreased fertility. Therefore, the objectives of this study were (i) to compute inbreeding coefficients using pedigree (F PED ) and genomic data based on runs of homozygosity (ROH) in the genome (F ROH ) of Austrian Pietrain pigs, and (ii) to assess inbreeding depression on four sperm quality traits. In total, 74,734 ejaculate records from 1034 Pietrain boars were used for inbreeding depression analyses. Traits were regressed on inbreeding coefficients using repeatability animal models. Pedigree-based inbreeding coefficients were lower than ROH-based inbreeding values. The correlations between pedigree and ROH-based inbreeding coefficients ranged from 0.186 to 0.357. Pedigree-based inbreeding affected only sperm motility while ROH-based inbreeding affected semen volume, number of spermatozoa, and motility. For example, a 1% increase in pedigree inbreeding considering 10 ancestor generations (F PED10 ) was significantly (p < 0.05) associated with a 0.231% decrease in sperm motility. Almost all estimated effects of inbreeding on the traits studied were unfavourable. It is advisable to properly manage the level of inbreeding to avoid high inbreeding depression in the future. Further, analysis of effects of inbreeding depression for other traits, including growth and litter size for the Austrian Pietrain population is strongly advised., (© 2023 John Wiley & Sons Ltd. Published by John Wiley & Sons Ltd.)

Ethiopia is a developing nation that could highly benefit from securing food via improved smallholder poultry farming. To support farmer and breeding decisions regarding which chicken strain to use in which Ethiopian environment, G*E analyses for body weight (BW) of growing male and female chickens were conducted. Research questions were (1) if a G*E is present for BW and (2) which strain performs best in which environment in terms of predicted BW. Analyses were performed using predicted BW at four different ages (90, 120, 150, and 180 days) of five strains (Horro, Koekoek, Kuroiler, Sasso-Rhode Island Red (S-RIR), and Sasso) tested in five Ethiopian regions (Addis Ababa, Amhara, Oromia, South Region, and Tigray) that are part of three Agro-Ecological Zones (AEZ) (cool humid, cool sub-humid, and warm semi-arid). The indigenous Horro strain was used as a control group to compare four other introduced tropically adapted strains. The dataset consisted of 999 female and 989 male farm-average BW measurements. G*E was strongly present ( p < 0.001) for all combinations of strain and region analyzed. In line with previous research, Sasso was shown to have the highest predicted BW, especially at an early age, followed by Kuroiler. Horro had the lowest predicted BW at most ages and in most regions, potentially due to its young breeding program. The highest predicted BW were observed in Tigray, Oromia, and Amhara regions, which are in the main part of the cool sub-humid AEZ.

Several factors, including breed, lead to divergent performance of pigs for production and reproduction traits in different environments. A recent genomics study showed that Modern European (ME) pig breeds contribute to the ancestry of smallholder pigs in the Hoima and Kamuli districts, Uganda. These pigs were also involved in a longitudinal study with several traits recorded, including 540 body weights (WT) of 374 growing pigs, 195 records of total number of piglets born alive (TBA) of 157 sows, and 110 total number weaned (TNW) records of 94 sows. Linear mixed-effects models were used to test for the significance of environmental effects, including housing system, geographic location, and the season when the events occurred as well as animal-specific effects like age, sex, parity, and farrow-to-weaning interval. Stepwise model reduction starting from models with all main effects and pairwise interactions was applied. The final models were then expanded to include proportions of Modern European (ME) ancestry for the subset of animals genotyped, following genomic ancestry analysis based on a Porcine 50K SNP Chip. ME ancestry proportions ranged from 0.02 to 0.50 and were categorized into three classes (low/medium/high ME) based on 33.3% quantiles. The effects of ME classes on WT and TBA were not significant. ME showed a significant effect on TNW. Sows with a high proportion of ME weaned 2.4 piglets more than the low group, the medium ME group being intermediate. This study used genomic data to investigate the effects of genetic ancestry on the performance of smallholder pigs in Uganda. The proportion of Modern European ancestry did not exceed 0.50, therefore not allowing for the comparison of local versus pure "exotic" types of pigs. For the range of ancestries observed, which is the relevant one for current smallholder systems in Uganda, differences were small for the body weight of growing pigs and the number of piglets born alive, while higher proportions of ME ancestry resulted in significantly more piglets weaned. The availability of genotypes of a higher number of growing pigs would have been beneficial for drawing conclusions on the effect of ME ancestry on the growth rates of smallholder pigs in Uganda., Competing Interests: Author CP was employed by PIG Austria GmbH. Author CL was employed by Pig Improvement Company. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Babigumira, Sölkner, Mészáros, Wurzinger, Pfeiffer, Lewis, Lukuyu, Ouma and Marshall.)

The African Goat Improvement Network (AGIN) is a collaborative group of scientists focused on genetic improvement of goats in small holder communities across the African continent. The group emerged from a series of workshops focused on enhancing goat productivity and sustainability. Discussions began in 2011 at the inaugural workshop held in Nairobi, Kenya. The goals of this diverse group were to: improve indigenous goat production in Africa; characterize existing goat populations and to facilitate germplasm preservation where appropriate; and to genomic approaches to better understand adaptation. The long-term goal was to develop cost-effective strategies to apply genomics to improve productivity of small holder farmers without sacrificing adaptation. Genome-wide information on genetic variation enabled genetic diversity studies, facilitated improved germplasm preservation decisions, and provided information necessary to initiate large scale genetic improvement programs. These improvements were partially implemented through a series of community-based breeding programs that engaged and empowered local small farmers, especially women, to promote sustainability of the production system. As with many international collaborative efforts, the AGIN work serves as a platform for human capacity development. This paper chronicles the evolution of the collaborative approach leading to the current AGIN organization and describes how it builds capacity for sustained research and development long after the initial program funds are gone. It is unique in its effectiveness for simultaneous, multi-level capacity building for researchers, students, farmers and communities, and local and regional government officials. The positive impact of AGIN capacity building has been felt by participants from developing, as well as developed country partners., Competing Interests: Author TSS was employed by Acceligen Inc. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Van Tassell, Rosen, Woodward-Greene, Silverstein, Huson, Sölkner, Boettcher, Rothschild, Mészáros, Nakimbugwe, Gondwe, Muchadeyi, Nandolo, Mulindwa, Banda, Kaumbata, Getachew, Haile, Soudre, Ouédraogo, Rischkowsky, Mwai, Dzomba, Nash, Abegaz, Masiga, Wurzinger, Sayre, Stella, Tosser-Klopp and Sonstegard.)

This study aimed to estimate the heterosis for productive traits in a two-way crossbreeding scheme. Four guinea pig lines were originally selected for the following traits: line P1 for the growth rate, P2 for the partial feed conversion rate, M1 for the growth rate of the litter at 10 days of age, and M2 for the litter size at birth. The comparison included 176 purebreds (P1: 46, P2: 43, M1: 54 and M2: 33) and 150 crosses (P1P2: 42, P2P1: 38, M1M2: 11 and M2M1: 59); body weights at birth, 10 days, weaning and 60 days of age were analyzed. A linear fixed-effect model was used, and heterosis was estimated as the difference between the average performance of the crossbred and pure-line animals. The pure line comparisons showed that P2 was lower than P1 for weight at 10 days and weaning weight, while all other comparisons between the paternal and maternal pure lines were not significant. The results indicated significant positive heterosis effects for both types of crosses, but only for birth weight: 3.7% for paternal crosses and 12.7% for maternal crosses. The heterosis estimates were mostly positive but not significant for all other traits. A reason for the low levels of heterosis could be that the lines are not very genetically differentiated. These results suggest that applying a two-way crossbreeding scheme within paternal and maternal guinea pig lines for meat production is not recommended due to the absence of heterosis for growth traits.

In Burkina Faso, goats are the second most numerous ruminant livestock population, with almost exclusively indigenous breeds being reared in extensive production systems in various agroecological zones. This study was carried out to understand the morphological variation of local goat breeds in the Sudano-Sahelian and Sudanian agroecological zones. A total of 511 adult female animals belonging to two presumed populations (Mossi breed in Sudano-Sahelian zone and Djallonké breed in Sudanian zone) were sampled and body weight as well as a range of linear body measurements, following FAO guidelines, were recorded. The least squares means of body measurements of indicated that Sudano-Sahelian goats have significantly ( p < 0.001) larger body measurements than Sudanian goats. Furthermore, relative high variability of the two populations in morphometric traits was observed. Principal Component Analysis (PCA) suggested structure between Mossi breed on one side and Djallonké on the other side, but no strict separation was observed, suggesting that gene flow is occurring among the different populations. A dispersion map with four clusters was built based on the first two factors. The least square means of body measurements ranked the four groups from small to large body size, namely Djallonké, Mossi × Djallonké, Mossi, and Sahelian × Mossi. Gene flow from Sahelian goat into other populations of the country, based on migration of the Fulani ethnic group from the Sahel into areas with Mossi and Djallonké breeds, could explain this configuration and confirms the continuous erosion of genetic identity of these two local breeds. The sustainable use of these adapted local goat genetic resources calls for the promotion of sustainable genetic improvement using participatory breeding approaches.

We used genome-wide SNP data from 18 local cattle breeds from six countries of the Alpine region to characterize population structure and identify genomic regions underlying positive selection. The geographically close breeds Evolèner, Eringer, Valdostana Pezzata Nera, and Valdostana Castana were found to differ from all other Alpine breeds. In addition, three breeds, Simmental, and Original Braunvieh from Switzerland and Pinzgauer from Austria built three separate clusters. Of the 18 breeds studied, the intra-alpine Swiss breed Evolèner had the highest average inbreeding based on runs of homozygosity (F ROH ) and the highest average genomic relationship within the breed. In contrast, Slovenian Cika cattle had the lowest average genomic inbreeding and the lowest average genomic relationship within the breed. We found selection signatures on chromosome 6 near known genes such as KIT and LCORL explaining variation in coat color and body size in cattle. The most prominent selection signatures were similar regardless of marker density and the breeds in the data set. In addition, using available high-density SNP data from 14 of the breeds we identified 47 genome regions as ROH islands. The proportion of homozygous animals was higher in all studied animals of local breeds than in Holstein and Brown Swiss cattle, the two most important commercial breeds in the Alpine region. We report ROH islands near genes related to thermoregulation, coat color, production, and stature. The results of this study serve as a basis for the search for causal variants underlying adaptation to the alpine environment and other specific characteristics selected during the evolution of local Alpine cattle breeds., (© 2023 The Authors. Animal Genetics published by John Wiley & Sons Ltd on behalf of Stichting International Foundation for Animal Genetics.)

Background: The unambiguous diagnosis of asphyxiation is still a major challenge for the forensic pathologist, especially in terms of highly advanced decomposed corps., Methodology: In order to demonstrate asphyxiation particularly in profoundly putrid bodies we hypothesized that hypoxic stress is basically responsible for generalized fatty degeneration of visceral organs which can be detected by histological examination using a special staining technique referred to as Oil-Red-O Stain (Sudan III-red-B-stain). To test this hypothesis we examined different tissues (myocardium, liver, lung and kidney) of 107 people divided into 5 groups. These are: (i) 71 case-victims who were found in a truck and died most likely due to asphyxiation, whereby any other violent or natural cause of death was ruled out by postmortem examination; (ii) 10 barely decomposed positive-control-victims; (iii) 6 non-decomposed positive-control-victims; iv) 10 drowning non-decomposed positive-control victims, and v) 10 negative-control-victims. Apart from general histological special staining methods, an immunohistochemically approach as a case-control-study on lung tissues of same individuals was carried out by means of using two polyclonal rabbit-antibodies against (i) HIF-1-α (Hypoxia Inducing Factor-1 alpha) and (ii) SP-A (pulmonary surfactant-associated protein A) to detect both the transcription factor and pulmonary surfactants. The positive proof of already either of them gives evidence of death caused by hypoxia., Results: Histological examination of myocardium, liver and kidney of the 71 case-victims and the 10 positive-control-victims using Oil-Red-O Stain showed a fatty degeneration of small droplet type; there was no evidence for fatty degeneration in tissues of the 10 negative-control-victims. These findings strongly indicate a causal association between oxygen deficiency and generalized fatty degeneration of viscera due to insufficient oxygen supply. In terms of methodology, this special staining technique seems to be very informative, even applicable on decomposed corps. Results of immunohistochemistry indicate that on the one hand the detection of HIF-1α is not possible to perform on (advanced) putrid bodies, whereas the verification of SP-A is still feasible on the other., Conclusion: Positive Oil-Red-O staining and the immunohistochemical detection of SP-A can serve as a serious hint for the diagnosis of asphyxia on putrid corpses, considering other circumstances of death that have been determined., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 Elsevier B.V. All rights reserved.)

Mid-infrared (MIR) spectroscopy is routinely applied to determine major milk components, such as fat and protein. Moreover, it is used to predict fine milk composition and various traits pertinent to animal health. MIR spectra indicate an absorbance value of infrared light at 1060 specific wavenumbers from 926 to 5010 cm -1 . According to research, certain parts of the spectrum do not contain sufficient information on traits of dairy cows. Hence, the objective of the present study was to identify specific regions of the MIR spectra of particular importance for the prediction of mastitis and ketosis, performing variable selection analysis. Partial least squares discriminant analysis (PLS-DA) along with three other statistical methods, support vector machine (SVM), least absolute shrinkage and selection operator (LASSO), and random forest (RF), were compared. Data originated from the Austrian milk recording and associated health monitoring system (GMON). Test-day data and corresponding MIR spectra were linked to respective clinical mastitis and ketosis diagnoses. Certain wavenumbers were identified as particularly relevant for the prediction models of clinical mastitis (23) and ketosis (61). Wavenumbers varied across four distinct statistical methods as well as concerning different traits. The results indicate that variable selection analysis could potentially be beneficial in the process of modeling.

Computer vision is a tool that could provide livestock producers with digital body measures and records that are important for animal health and production, namely body height and length, and chest girth. However, to build these tools, the scarcity of labeled training data sets with uniform images (pose, lighting) that also represent real-world livestock can be a challenge. Collecting images in a standard way, with manual image labeling is the gold standard to create such training data, but the time and cost can be prohibitive. We introduce the PreciseEdge image segmentation algorithm to address these issues by employing a standard image collection protocol with a semi-automated image labeling method, and a highly precise image segmentation for automated body measurement extraction directly from each image. These elements, from image collection to extraction are designed to work together to yield values highly correlated to real-world body measurements. PreciseEdge adds a brief preprocessing step inspired by chromakey to a modified GrabCut procedure to generate image masks for data extraction (body measurements) directly from the images. Three hundred RGB (red, green, blue) image samples were collected uniformly per the African Goat Improvement Network Image Collection Protocol (AGIN-ICP), which prescribes camera distance, poses, a blue backdrop, and a custom AGIN-ICP calibration sign. Images were taken in natural settings outdoors and in barns under high and low light, using a Ricoh digital camera producing JPG images (converted to PNG prior to processing). The rear and side AGIN-ICP poses were used for this study. PreciseEdge and GrabCut image segmentation methods were compared for differences in user input required to segment the images. The initial bounding box image output was captured for visual comparison. Automated digital body measurements extracted were compared to manual measures for each method. Both methods allow additional optional refinement (mouse strokes) to aid the segmentation algorithm. These optional mouse strokes were captured automatically and compared. Stroke count distributions for both methods were not normally distributed per Kolmogorov-Smirnov tests. Non-parametric Wilcoxon tests showed the distributions were different (p< 0.001) and the GrabCut stroke count was significantly higher (p = 5.115 e-49), with a mean of 577.08 (std 248.45) versus 221.57 (std 149.45) with PreciseEdge. Digital body measures were highly correlated to manual height, length, and girth measures, (0.931, 0.943, 0.893) for PreciseEdge and (0.936, 0. 944, 0.869) for GrabCut (Pearson correlation coefficient). PreciseEdge image segmentation allowed for masks yielding accurate digital body measurements highly correlated to manual, real-world measurements with over 38% less user input for an efficient, reliable, non-invasive alternative to livestock hand-held direct measuring tools., Competing Interests: The authors have declared that no competing interests exist.

Cattle have been essential for the development of human civilization since their first domestication few thousand years ago. Since then, they have spread across vast geographic areas following human activities. Throughout generations, the cattle genome has been shaped with detectable signals induced by various evolutionary processes, such as natural and human selection processes and demographic events. Identifying such signals, called selection signatures, is one of the primary goals of population genetics. Previous studies used various selection signature methods and normalized the outputs score using specific windows, in kbp or based on the number of SNPs, to identify the candidate regions. The recent method of iSAFE claimed for high accuracy in pinpointing the candidate SNPs. In this study, we analyzed whole-genome resequencing (WGS) data of ten individuals from Austrian Fleckvieh ( Bos taurus ) and fifty individuals from 14 Chinese indigenous breeds ( Bos taurus, Bos taurus indicus, and admixed). Individual WGS reads were aligned to the cattle reference genome of ARS. UCD1.2 and subsequently undergone single nucleotide variants (SNVs) calling pipeline using GATK. Using these SNVs, we examined the population structure using principal component and admixture analysis. Then we refined selection signature candidates using the iSAFE program and compared it with the classical iHS approach. Additionally, we run Fst population differentiation from these two cattle groups. We found gradual changes of taurine in north China to admixed and indicine to the south. Based on the population structure and the number of individuals, we grouped samples to Fleckvieh, three Chinese taurines (Kazakh, Mongolian, Yanbian), admixed individuals (CHBI&#95;Med), indicine individuals (CHBI&#95;Low), and a combination of admixed and indicine (CHBI) for performing iSAFE and iHS tests. There were more significant SNVs identified using iSAFE than the iHS for the candidate of positive selection and more detectable signals in taurine than in indicine individuals. However, combining admixed and indicine individuals decreased the iSAFE signals. From both within-population tests, significant SNVs are linked to the olfactory receptors, production, reproduction, and temperament traits in taurine cattle, while heat and parasites tolerance in the admixed individuals. Fst test suggests similar patterns of population differentiation between Fleckvieh and three Chinese taurine breeds against CHBI. Nevertheless, there are genes shared only among the Chinese taurine, such as PAX5, affecting coat color, which might drive the differences between these yellowish coated breeds, and those in the greater Far East region., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Naji, Jiang, Utsunomiya, Rosen, Sölkner, Wang, Jiang, Zhang, Zhang, Ding and Mészáros.)

By their paternal transmission, Y-chromosomal haplotypes are sensitive markers of population history and male-mediated introgression. Previous studies identified biallelic single-nucleotide variants in the SRY, ZFY and DDX3Y genes, which in domestic goats identified four major Y-chromosomal haplotypes, Y1A, Y1B, Y2A and Y2B, with a marked geographical partitioning. Here, we extracted goat Y-chromosomal variants from whole-genome sequences of 386 domestic goats (75 breeds) and seven wild goat species, which were generated by the VarGoats goat genome project. Phylogenetic analyses indicated domestic haplogroups corresponding to Y1B, Y2A and Y2B, respectively, whereas Y1A is split into Y1AA and Y1AB. All five haplogroups were detected in 26 ancient DNA samples from southeast Europe or Asia. Haplotypes from present-day bezoars are not shared with domestic goats and are attached to deep nodes of the trees and networks. Haplogroup distributions for 186 domestic breeds indicate ancient paternal population bottlenecks and expansions during migrations into northern Europe, eastern and southern Asia, and Africa south of the Sahara. In addition, sharing of haplogroups indicates male-mediated introgressions, most notably an early gene flow from Asian goats into Madagascar and the crossbreeding that in the 19th century resulted in the popular Boer and Anglo-Nubian breeds. More recent introgressions are those from European goats into the native Korean goat population and from Boer goat into Uganda, Kenya, Tanzania, Malawi and Zimbabwe. This study illustrates the power of the Y-chromosomal variants for reconstructing the history of domestic species with a wide geographical range., (© 2022 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

Monitoring for mastitis on dairy farms is of particular importance, as it is one of the most prevalent bovine diseases. A commonly used indicator for mastitis monitoring is somatic cell count. A supplementary tool to predict mastitis risk may be mid-infrared (MIR) spectroscopy of milk. Because bovine health status can affect milk composition, this technique is already routinely used to determine standard milk components. The aim of the present study was to compare the performance of models to predict clinical mastitis based on MIR spectral data and/or somatic cell count score (SCS), and to explore differences of prediction accuracies for acute and chronic clinical mastitis diagnoses. Test-day data of the routine Austrian milk recording system and diagnosis data of its health monitoring, from 59,002 cows of the breeds Fleckvieh (dual purpose Simmental), Holstein Friesian and Brown Swiss, were used. Test-day records within 21 days before and 21 days after a mastitis diagnosis were defined as mastitis cases. Three different models (MIR, SCS, MIR + SCS) were compared, applying Partial Least Squares Discriminant Analysis. Results of external validation in the overall time window (-/+21 days) showed area under receiver operating characteristic curves (AUC) of 0.70 when based only on MIR, 0.72 when based only on SCS, and 0.76 when based on both. Considering as mastitis cases only the test-day records within 7 days after mastitis diagnosis, the corresponding areas under the curve were 0.77, 0.83 and 0.85. Hence, the model combining MIR spectral data and SCS was performing best. Mastitis probabilities derived from the prediction models are potentially valuable for routine mastitis monitoring for farmers, as well as for the genetic evaluation of the trait udder health.

Nguni cattle are a Sanga type breed with mixed B. taurus and B. indicus ancestry and proven resistance to ticks, diseases and other harsh conditions of the African geographical landscape. The multi-coloured Nguni coats have found a niche market in the leather industry leading to breeding objectives towards the promotion of such diversity. However, there is limited studies on the genomic architecture underlying the coat colour and patterns hampering any potential breeding and improvement of such trait. This study investigated the genetics of base coat colour, colour-sidedness and the white forehead stripe in Nguni cattle using coat colour phenotyped Nguni cattle and Illumina Bovine HD (770K) genotypes. Base coat colour phenotypes were categorised into eumelanin ( n = 45) and pheomelanin ( n = 19). Animals were categorised into either colour-sided ( n = 46) or non-colour-sided ( n = 94) and similarly into presence ( n = 15) or absence ( n = 67) of white forehead stripe. Genome-wide association tests were conducted using 622,103 quality controlled SNPs and the Efficient Mixed Model Association eXpedited method (EMMAX) implemented in Golden Helix SNP Variation Suite. The genome-wide association studies for base coat colour (eumelanin vs. pheomelanin) resulted into four indicative SNPs on BTA18 and a well-known gene, MC1R, was observed within 1 MB from the indicative SNPs (p < 0.00001) and found to play a role in the melanogenesis (core pathway for melanin production) and the MAPK signalling pathway. GWAS for colour-sidedness resulted in four indicative SNPs, none of which were in close proximity to the KIT candidate gene known for colour-sidedness. GWAS for the white forehead stripe resulted in 17 indicative SNPs on BTA6. Four genes MAPK10, EFNA5, PPP2R3C and PAK1 were found to be associated with the white forehead stripe and were part of the MAPK, adrenergic and Wnt signalling pathways that are synergistically associated with the synthesis of melanin. Overall, our results prove prior knowledge of the role of MC1R in base coat colours in cattle and suggested a different genetic mechanism for forehead stripe phenotypes in Nguni cattle., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Kunene, Muchadeyi, Hadebe, Mészáros, Sölkner, Dugmore and Dzomba.)