On the genomic regions associated with milk lactose in Fleckvieh cattle

Sequence-based GWAS in 180,000 German Holstein cattle reveals new candidate variants for milk production traits

BACKGROUND

Milk production traits are complex and influenced by many genetic and environmental factors. Although extensive research has been performed for these traits, with many associations unveiled thus far, due to their crucial economic importance, complex genetic architecture, and the fact that causal variants in cattle are still scarce, there is a need for a better understanding of their genetic background. In this study, we aimed to identify new candidate loci associated with milk production traits in German Holstein cattle, the most important dairy breed in Germany and worldwide. For that purpose, 180,217 cattle were imputed to the sequence level and large-scale genome-wide association study (GWAS) followed by fine-mapping and evolutionary and functional annotation were carried out to identify and prioritize new association signals.

RESULTS

Using the imputed sequence data of a large cattle dataset, we identified 50,876 significant variants, confirming many known and identifying previously unreported candidate variants for milk (MY), fat (FY), and protein yield (PY). Genome-wide significant signals were fine-mapped with the Bayesian approach that determines the credible variant sets and generates the probability of causality for each signal. The variants with the highest probabilities of being causal were further classified using external information about the function and evolution, making the prioritization for subsequent validation experiments easier. The top potential causal variants determined with fine-mapping explained a large percentage of genetic variance compared to random ones; 178 variants explained 11.5%, 104 explained 7.7%, and 68 variants explained 3.9% of the variance for MY, FY, and PY, respectively, demonstrating the potential for causality.

CONCLUSIONS

Our findings proved the power of large samples and sequence-based GWAS in detecting new association signals. In order to fully exploit the power of GWAS, one should aim at very large samples combined with whole-genome sequence data. These can also come with both computational and time burdens, as presented in our study. Although milk production traits in cattle are comprehensively investigated, the genetic background of these traits is still not fully understood, with the potential for many new associations to be revealed, as shown. With constantly growing sample sizes, we expect more insights into the genetic architecture of milk production traits in the future.

Mastitis has a cumulative and lasting effect on milk yield and lactose content in dairy cows

Milk lactose content (LC) physiologically decreases with parity order in dairy cows, but also after udder health inflammation(s) and in presence of elevated milk SCC in subclinical cases. Therefore, the progressive decrease in milk LC observed along cows' productive life can be attributed to a combination of factors that altogether impair the epithelial integrity, resulting in weaker tight junctions, e.g., physiological aging of epithelium, mechanical epithelial stress due to milking, and experienced clinical or subclinical mastitis. Mastitis is also known to affect the udder synthesis ability, so our intention through this study was to evaluate if there is a cumulative and lasting effect of mammary gland inflammation(s) on milk yield (MY) and LC. For this purpose, we used diagnoses of clinical mastitis and milk data of Austrian Fleckvieh cows to evaluate the effect of cumulative mastitis events on LC and MY. Only mastitis diagnoses recorded by trained veterinarians were used. Finally, we investigated if cumulative mastitis is a heritable trait and whether it is genetically correlated with either LC or MY. Estimates were obtained using univariate and bivariate linear animal models. A significant reduction in LC and MY was observed in cows that suffered from mastitis compared with those that did not experience udder inflammation. The h2 of cumulative mastitis is promising and much greater (0.09) than the h2 of the binary event itself (≤0.03). The genetic correlations between cumulative mastitis with LC and MY were negative, suggesting that cows with a great genetic merit for MY and LC are expected to be more resistant to repeated inflammations and less recidivist. When we used number of lifetime SCC peaks (≥200,000 or 400,000 cells/mL) to calculate cumulative inflammation events, h2 was even higher (up to 0.38), implying that subclinical mastitis also has a relevant negative impact on both LC and MY. Finally, the present study demonstrated how repeated mastitis events can permanently affect the mammary gland epithelial integrity and synthesis ability, and that the number of cumulative mastitis is a promising phenotype to be used in selection index in combination with other indicator traits toward more resistant and resilient mammary glands.

A Pilot Detection and Associate Study of Gene Presence-Absence Variation in Holstein Cattle

Presence-absence variations (PAVs) are important structural variations, wherein a genomic segment containing one or more genes is present in some individuals but absent in others. While PAVs have been extensively studied in plants, research in cattle remains limited. This study identified PAVs in 173 Holstein bulls using whole-genome sequencing data and assessed their associations with 46 economically important traits. Out of 28,772 cattle genes (from the longest transcripts), a total of 26,979 (93.77%) core genes were identified (present in all individuals), while variable genes included 928 softcore (present in 95-99% of individuals), 494 shell (present in 5-94%), and 371 cloud genes (present in <5%). Cloud genes were enriched in functions associated with hormonal and antimicrobial activities, while shell genes were enriched in immune functions. PAV-based genome-wide association studies identified associations between gene PAVs and 16 traits including milk, fat, and protein yields, as well as traits related to health and reproduction. Associations were found on multiple chromosomes, illustrating important associations on cattle chromosomes 7 and 15, involving olfactory receptor and immune-related genes, respectively. By examining the PAVs at the population level, the results of this research provided crucial insights into the genetic structures underlying the complex traits of Holstein cattle.

The effect of first-lactation calving season, milk production, and morphology on the survival of Simmental cows

Longevity in dairy and dual-purpose cattle is a complex trait which depends on many individual and managerial factors. The purpose of the present study was to investigate the survival (SURV) rate of Italian Simmental dual-purpose cows across different parities. Data of this study referred to 2 173 primiparous cows under official milk recording that calved between 2002 and 2020. Only cows linearly classified for type traits, including muscularity (MU) and body condition score (BCS) were kept. Survival analysis was carried out, through the Cox regression model, for different pairwise combinations of classes of milk productivity MU, BCS, and calving season. Herd-year of first calving was also considered in the model. SURV (0 = culled; 1 = survived) at each lactation up to the 6th were the dependent variables, so that, for example, SURV2 equal to 1 was attributed to cows that entered the 2nd lactation. Survival rates were 98, 71, 63, 56, and 53% for 2nd, 3rd, 4th, 5th, and 6th lactation, respectively. Results revealed that SURV2 was not dependent on milk yield, while in subsequent parities, low-producing cows were characterized by higher SURV compared to high-producing ones. Additionally, cows starting the lactation in autumn survived less (47.38%) than those starting in spring (53.49%), suggesting that facing the late gestation phase in summer could increase the culling risk. The present study indicates that SURV in Italian Simmental cows is influenced by various factors in addition to milk productivity. However, it is important to consider that in this study all first-calving cows culled before the linear evaluation - carried out between mid- and late lactation in this breed - were not accounted for. Finding can be transferred to other dual-purpose breeds, where the cows' body conformation and muscle development - i.e. meat-related features - are often considered as important as milk performance by farmers undertaking culling decisions.

Single and longitudinal genome-wide association studies for dairy traits available in goats with three recorded lactations

Milk yield and composition phenotypes are systematically recorded across several lactations in goats, but the majority of genome-wide association studies (GWAS) performed so far have rather ignored the longitudinal nature of such data. Here, we have used two different GWAS approaches to analyse data from three lactations recorded in Murciano-Granadina goats. In Analysis 1, independent GWAS have been carried out for each trait and lactation, while a single longitudinal GWAS, jointly considering all data, has been performed in Analysis 2. In both analyses, genome-wide significant QTL for lactose percentage on chromosome 2 (129.77-131.01 Mb) and for milk protein percentage on the chromosome 6 (74.8-94.6 Mb) casein gene cluster region were detected. In Analysis 1, several QTL were not replicated in all three lactations, possibly due to the existence of lactation-specific genetic determinants. In Analysis 2, we identified several genome-wide significant QTL related to milk yield and protein content that were not uncovered in Analysis 1. The increased number of QTL identified in Analysis 2 suggests that the longitudinal GWAS is particularly well suited for the genetic analysis of dairy traits. Moreover, our data confirm that variability within or close to the casein complex is the main genetic determinant of milk protein percentage in Murciano-Granadina goats.

Enhancing Genomic Prediction Accuracy for Body Conformation Traits in Korean Holstein Cattle

The Holstein breed is the mainstay of dairy production in Korea. In this study, we evaluated the genomic prediction accuracy for body conformation traits in Korean Holstein cattle, using a range of π levels (0.75, 0.90, 0.99, and 0.995) in Bayesian methods (BayesB and BayesC). Focusing on 24 traits, we analyzed the impact of different π levels on prediction accuracy. We observed a general increase in accuracy at higher levels for specific traits, with variations depending on the Bayesian method applied. Notably, the highest accuracy was achieved for rear teat angle when using deregressed estimated breeding values including parent average as a response variable. We further demonstrated that incorporating parent average into deregressed estimated breeding values enhances genomic prediction accuracy, showcasing the effectiveness of the model in integrating both offspring and parental genetic information. Additionally, we identified 18 significant window regions through genome-wide association studies, which are crucial for future fine mapping and discovery of causal mutations. These findings provide valuable insights into the efficiency of genomic selection for body conformation traits in Korean Holstein cattle and highlight the potential for advancements in the prediction accuracy using larger datasets and more sophisticated genomic models.

Single-step genome-wide association analyses for selected infrared-predicted cheese-making traits in Walloon Holstein cows

This study aimed to perform genome-wide association study to identify genomic regions associated with milk production and cheese-making properties (CMP) in Walloon Holstein cows. The studied traits were milk yield, fat percentage, protein percentage, casein percentage (CNP), calcium content, somatic cell score (SCS), coagulation time, curd firmness after 30 min from rennet addition, and titratable acidity. The used data have been collected from 2014 to 2020 on 78,073 first-parity (485,218 test-day records), 48,766 second-parity (284,942 test-day records), and 21,948 third-parity (105,112 test-day records) Holstein cows distributed in 671 herds in the Walloon Region of Belgium. Data of 565,533 single nucleotide polymorphisms (SNP), located on 29 Bos taurus autosomes (BTA) of 6,617 animals (1,712 males), were used. Random regression test-day models were used to estimate genetic parameters through the Bayesian Gibbs sampling method. The SNP solutions were estimated using a single-step genomic BLUP approach. The proportion of the total additive genetic variance explained by windows of 50 consecutive SNPs (with an average size of ∼216 KB) was calculated, and regions accounting for at least 1.0% of the total additive genetic variance were used to search for positional candidate genes. Heritability estimates for the studied traits ranged from 0.10 (SCS) to 0.53 (CNP), 0.10 (SCS) to 0.50 (CNP), and 0.12 (SCS) to 0.49 (CNP) in the first, second, and third parity, respectively. Genome-wide association analyses identified 6 genomic regions (BTA1, BTA14 [4 regions], and BTA20) associated with the considered traits. Genes including the SLC37A1 (BTA1), SHARPIN, MROH1, DGAT1, FAM83H, TIGD5, MROH6, NAPRT, ADGRB1, GML, LYPD2, JRK (BTA14), and TRIO (BTA20) were identified as positional candidate genes for the studied CMP. The findings of this study help to unravel the genomic background of a cow's ability for cheese production and can be used for the future implementation and use of genomic evaluation to improve the cheese-making traits in Walloon Holstein cows.

Understanding the phenotypic and genetic background of the lactose content in Sarda dairy sheep

Lactose, the principal carbohydrate found in milk, plays an important role in the physiological processes of milk production because it is related to milk volume, and it is responsible for the osmotic equilibrium between blood and milk in the mammary gland. In this study, factors affecting lactose content (LC) in sheep milk are investigated. For this purpose, 2,358 test-day records were sampled from 509 ewes (3-7 records per animal). The LC and other main milk traits were analyzed using a mixed linear model that included days in milk (DIM) class, parity, lambing month, and type of lambing as fixed effects and animal, permanent environment, and flock test day as random effects. The pedigree-based approach was used to estimate the heritability and repeatability of LC. Moreover, the genomic background of LC was investigated through a GWAS. The LC was affected by all tested factors (i.e., DIM class, parity, lambing month, and type of lambing). Low heritability (0.10 ± 0.05) and moderate repeatability (0.42 ± 0.02) were estimated for LC. High negative genetic correlations were estimated between LC and NaCl (-0.99 ± 0.01) and between LC and somatic cell count (-0.94 ± 0.05). Only 2 markers passed the chromosome-wide Bonferroni threshold. Results of the present study, although obtained on a relatively small sample, suggest the possibility to include LC in the breeding programs, particularly because of its strong relationship with NaCl and somatic cell count.

Analysis of runs of homozygosity of cattle living in different climate zones

Aim of this study was to analyze the distribution and characteristics of runs of homozygosity in Bos taurus taurus and Bos taurus indicus breeds, as well as their crosses, farmed all around the world. With this aim in view, we used single-nucleotide polymorphisms (SNP) genotypes for 3,263 cattle belonging to 204 different breeds. After quality control, 23,311 SNPs were retained for the analysis. Animals were divided into seven different groups: 1) continental taurus, 2) temperate taurus, 3) temperate indicus, 4) temperate composite, 5) tropical taurus, 6) tropical indicus, and 7) tropical composite. The climatic zones were created according to the latitude of the breeds' country of origin: i) continental, latitude ≥ 45°; ii) temperate, 45°< Latitude >23.26°; iii) tropics, latitude ≤ 23.26°. Runs of homozygosity were computed as 15 SNPs spanning in at least 2 Mb; number of ROH per animal (nROH), average ROH length (meanMb), and ROH-based inbreeding coefficients (FROH) were also computed. Temperate indicus showed the largest nROH, whereas Temperate taurus the lowest value. Moreover, the largest meanMb was observed for Temperate taurus, whereas the lowest value for Tropics indicus. Temperate indicus breeds showed the largest FROH values. Genes mapped in the identified ROH were reported to be associated with the environmental adaptation, disease resistance, coat color determinism, and production traits. Results of the present study confirmed that runs of homozygosity could be used to identify genomic signatures due to both artificial and natural selection.

Genomic characterization of autozygosity and recent inbreeding trends in all major breeds of US dairy cattle

Maintaining a genetically diverse dairy cattle population is critical to preserving adaptability to future breeding goals and avoiding declines in fitness. This study characterized the genomic landscape of autozygosity and assessed trends in genetic diversity in 5 breeds of US dairy cattle. We analyzed a sizable genomic data set containing 4,173,679 pedigreed and genotyped animals of the Ayrshire, Brown Swiss, Guernsey, Holstein, and Jersey breeds. Runs of homozygosity (ROH) of 2 Mb or longer in length were identified in each animal. The within-breed means for number and the combined length of ROH were highest in Jerseys (62.66 ± 8.29 ROH and 426.24 ± 83.40 Mb, respectively; mean ± SD) and lowest in Ayrshires (37.24 ± 8.27 ROH and 265.05 ± 85.00 Mb, respectively). Short ROH were the most abundant, but moderate to large ROH made up the largest proportion of genome autozygosity in all breeds. In addition, we identified ROH islands in each breed. This revealed selection patterns for milk production, productive life, health, and reproduction in most breeds and evidence for parallel selective pressure for loci on chromosome 6 between Ayrshire and Brown Swiss and for loci on chromosome 20 between Holstein and Jersey. We calculated inbreeding coefficients using 3 different approaches, pedigree-based (F_PED), marker-based using a genomic relationship matrix (F_GRM), and segment-based using ROH (F_ROH). The average inbreeding coefficient ranged from 0.06 in Ayrshires and Brown Swiss to 0.08 in Jerseys and Holsteins using F_PED, from 0.22 in Holsteins to 0.29 in Guernsey and Jerseys using F_GRM, and from 0.11 in Ayrshires to 0.17 in Jerseys using F_ROH. In addition, the effective population size at past generations (5-100 generations ago), the yearly rate of inbreeding, and the effective population size in 3 recent periods (2000-2009, 2010-2014, and 2015-2018) were determined in each breed to ascertain current and historical trends of genetic diversity. We found a historical trend of decreasing effective population size in the last 100 generations in all breeds and breed differences in the effect of the recent implementation of genomic selection on inbreeding accumulation.

Genetic parameters and genome-wide association for milk production traits and somatic cell score in different lactation stages of Shanghai Holstein population

The aim of this study was to investigate the genetic parameters and genetic architectures of six milk production traits in the Shanghai Holstein population. The data used to estimate the genetic parameters consisted of 1,968,589 test-day records for 305,031 primiparous cows. Among the cows with phenotypes, 3,016 cows were genotyped with Illumina Bovine SNP50K BeadChip, GeneSeek Bovine 50K BeadChip, GeneSeek Bovine LD BeadChip v4, GeneSeek Bovine 150K BeadChip, or low-depth whole-genome sequencing. A genome-wide association study was performed to identify quantitative trait loci and genes associated with milk production traits in the Shanghai Holstein population using genotypes imputed to whole-genome sequences and both fixed and random model circulating probability unification and a mixed linear model with rMVP software. Estimated heritabilities (h2) varied from 0.04 to 0.14 for somatic cell score (SCS), 0.07 to 0.22 for fat percentage (FP), 0.09 to 0.27 for milk yield (MY), 0.06 to 0.23 for fat yield (FY), 0.09 to 0.26 for protein yield (PY), and 0.07 to 0.35 for protein percentage (PP), respectively. Within lactation, genetic correlations for SCS, FP, MY, FY, PY, and PP at different stages of lactation estimated in random regression model were ranged from -0.02 to 0.99, 0.18 to 0.99, 0.04 to 0.99, 0.04 to 0.99, 0.01 to 0.99, and 0.33 to 0.99, respectively. The genetic correlations were highest between adjacent DIM but decreased as DIM got further apart. Candidate genes included those related to production traits (DGAT1, MGST1, PTK2, and SCRIB), disease-related (LY6K, COL22A1, TECPR2, and PLCB1), heat stress-related (ITGA9, NDST4, TECPR2, and HSF1), and reproduction-related (7SK and DOCK2) genes. This study has shown that there are differences in the genetic mechanisms of milk production traits at different stages of lactation. Therefore, it is necessary to conduct research on milk production traits at different stages of lactation as different traits. Our results can also provide a theoretical basis for subsequent molecular breeding, especially for the novel genetic loci.

Genome-wide association study of beef bull semen attributes

Background

Cattle production is dependent upon fertility because it results in producing offspring to offset production costs. A number of semen attributes are believed to affect fertility and are frequently measured as part of routine breeding soundness exams or semen collection procedures. The objective of this study was to perform a single-step genome-wide association study (ssGWAS) for beef bull semen attributes. Beef bull fertility phenotypes including volume (VOL), concentration (CONC), number of spermatozoa (NSP), initial motility (IMot), post-thaw motility (PTMot), three-hour post-thaw motility (3HRPTMot), percentage of normal spermatozoa (%NORM), primary abnormalities (PRIM), and secondary abnormalities (SEC) were obtained from two artificial insemination (AI) centers. A total of 1819 Angus bulls with 50,624 collection records were used for ssGWAS. A five-generation pedigree was obtained from the American Angus Association and consisted of 6521 sires and 17,136 dams. Genotypes on 1163 bulls were also obtained from the American Angus Association and utilized in ssGWAS.

Results

A multi-trait animal model was used for the estimation of single nucleotide polymorphism (SNP) effects. Significant SNP were those with a -log10P-value threshold greater than 4.0. Volume, CONC, NSP, IMot, PTMot, 3HRPTMot, %NORM, PRIM, and SEC have five, three, six, seven, two, six, six, and two genome-wide significant SNP, respectively.

Conclusions

Several significant SNP were determined to be near or within quantitative trait loci (QTL) associated with beef bull semen attributes. In addition, genes associated with fertility were found to contain or be near the significant SNP found in the study. The results indicate there are regions of the genome that impact fertility, proving inclusion of genomic information into genetic evaluation should be advantageous for genetic improvement of male fertility traits.

Genome-wide association study for meat tenderness in beef cattle identifies patterns of the genetic contribution in different post-mortem stages

The beef tenderization process during the post-mortem period is one of the most important sensorial attributes and it is well-established. The aim of this study was to identify the genetic contribution pattern to meat tenderness at 7-(LMD7), 14-(LMD14), and 21-(LMD21) days post-mortem. The heritabilities for LMD7 (0.194), LMD14 (0.142) and LMD21 (0.048) are well established in the population evaluated here. However, its genetic contribution in terms of genomic candidate regions is still poorly understood. Tenderness was measured in the Longissiums thoracis using Warner-Bratzler shear force in the three post-mortem periods. A total of 4323 crossbred beef cattle were phenotyped and genotyped using the Illumina BovineSNP50K. The percentage of the total genetic variance was estimated using the weighted single-step genomic best linear unbiased prediction method. The main candidate windows for LMD7 were associated with proteolysis of myofibrillar structures and the weakening endomysium and perimysium. Candidate windows for LMD14 and LMD21 were mapped in bovine QTLs for body composition, height and growth. Results presented herein highlight, the largest contribution of proteolysis related processes before 14-days post-mortem and body composition characteristics in later stages for meat tenderness.

Investigation of β-hydroxybutyrate in early lactation of Simmental cows: Genetic parameters and genomic predictions

Genomic information allows for a more accurate calculation of relationships among animals than the pedigree information, leading to an increase in accuracy of breeding values. Here, we used pedigree-based and single-step genomic approaches to estimate variance components and breeding values for β-hydroxybutyrate milk content (BHB). Additionally, we performed a genome-wide association study (GWAS) to depict its genetic architecture. BHB concentrations within the first 90 days of lactation, estimated from milk medium infrared spectra, were available for 30,461 cows (70,984 records). Genotypes at 42,152 loci were available for 9,123 animals. Low heritabilities were found for BHB using pedigree-based (0.09 ± 0.01) and genomic (0.10 ± 0.01) approaches. Genetic correlation between BHB and milk traits ranged from -0.27 ± 0.06 (BHB and protein percentage) to 0.13 ± 0.07 (BHB and fat-to-protein ratio) using pedigree and from -0.26 ± 0.05 (BHB and protein percentage) to 0.13 ± 0.06 (BHB and fat-to-protein ratio) using genomics. Breeding values were validated for 344 genotyped cows using linear regression method. The genomic EBV (GEBV) had greater accuracy (0.51 vs. 0.45) and regression coefficient (0.98 vs. 0.95) compared to EBV. The correlation between two subsequent evaluations, without and with phenotypes for validation cows, was 0.85 for GEBV and 0.82 for EBV. Predictive ability (correlation between (G)EBV and adjusted phenotypes) was greater when genomic information was used (0.38) than in the pedigree-based approach (0.31). Validation statistics in the pairwise two-trait models (milk yield, fat and protein percentage, urea, fat/protein ratio, lactose and logarithmic transformation of somatic cells count) were very similar to the ones highlighted for the single-trait model. The GWAS allowed discovering four significant markers located on BTA20 (57.5-58.2 Mb), where the ANKH gene is mapped. This gene has been associated with lactose, alpha-lactalbumin and BHB. Results of this study confirmed the usefulness of genomic information to provide more accurate variance components and breeding values, and important insights about the genomic determination of BHB milk content.

Genome-wide association study on milk production and somatic cell score for Thai dairy cattle using weighted single-step approach with random regression test-day model

Genome-wide association studies are a powerful tool to identify genomic regions and variants associated with phenotypes. However, only limited mutual confirmation from different studies is available. The objectives of this study were to identify genomic regions as well as genes and pathways associated with the first-lactation milk, fat, protein, and total solid yields; fat, protein, and total solid percentage; and somatic cell score (SCS) in a Thai dairy cattle population. Effects of SNPs were estimated by a weighted single-step GWAS, which back-solved the genomic breeding values predicted using single-step genomic BLUP (ssGBLUP) fitting a single-trait random regression test-day model. Genomic regions that explained at least 0.5% of the total genetic variance were selected for further analyses of candidate genes. Despite the small number of genotyped animals, genomic predictions led to an improvement in the accuracy over the traditional BLUP. Genomic predictions using weighted ssGBLUP were slightly better than the ssGBLUP. The genomic regions associated with milk production traits contained 210 candidate genes on 19 chromosomes [Bos taurus autosome (BTA) 1 to 7, 9, 11 to 16, 20 to 21, 26 to 27 and 29], whereas 21 candidate genes on 3 chromosomes (BTA 11, 16, and 21) were associated with SCS. Many genomic regions explained a small fraction of the genetic variance, indicating polygenic inheritance of the studied traits. Several candidate genes coincided with previous reports for milk production traits in Holstein cattle, especially a large region of genes on BTA14. We identified 141 and 5 novel genes related to milk production and SCS, respectively. These novel genes were also found to be functionally related to heat tolerance (e.g., SLC45A2, IRAG1, and LOC101902172), longevity (e.g., SYT10 and LOC101903327), and fertility (e.g., PAG1). These findings may be attributed to indirect selection in our population. Identified biological networks including intracellular cell transportation and protein catabolism implicate milk production, whereas the immunological pathways such as lymphocyte activation are closely related to SCS. Further studies are required to validate our findings before exploiting them in genomic selection.

Rule Discovery in Milk Content towards Mastitis Diagnosis: Dealing with Farm Heterogeneity over Multiple Years through Classification Based on Associations

Simple Summary

Invisible (subclinical) mastitis decreases milk quality and production. Invisible mastitis is linked to an increased use of antimicrobials. The risk of the emergence of antimicrobial-resistant bacteria is a major public health concern worldwide. Therefore, early detection of infected cows is of great importance. Machine learning has opened a new avenue for early mastitis prediction based on simple and accessible milking parameters, such as milk volume, fat, protein, lactose, electrical conductivity (EC), milking time, and milking peak flow. However, farm heterogeneity is a major challenge where multiple patterns can predict mastitis. Here, we employed a classification based on associations and scaling approach for multiple pattern discovery over multiple years. The approach we have developed helps to address farm heterogeneity and generalise machine learning-based diagnosis of mastitis worldwide.

Abstract

Subclinical mastitis, an economically challenging disease of dairy cattle, is associated with an increased use of antimicrobials which reduces milk quantity and quality. It is more common than clinical mastitis and far more difficult to detect. Recently, much attention has been paid to the development of machine-learning expert systems for early detection of subclinical mastitis from milking features. However, differences between animals within a farm as well as between farms, particularly across multiple years, are major obstacles to the generalisation of machine learning models. Here, for the first time, we integrated scaling by quartiling with classification based on associations in a multi-year study to deal with farm heterogeneity by discovery of multiple patterns towards mastitis. The data were obtained from one farm comprising Holstein Friesian cows in Ongaonga, New Zealand, using an electronic automated monitoring system. The data collection was repeated annually over 3 consecutive years. Some discovered rules, such as when the milking peak flow is low, electrical conductivity (EC) of milk is low, milk lactose is low, milk fat is high, and milk volume is low, the cow has subclinical mastitis, reached high confidence (>70%) in multiple years. On averages, over 3 years, low level of milk lactose and high value of milk EC were part of 93% and 83.8% of all subclinical mastitis detecting rules, offering a reproducible pattern of subclinical mastitis detection. The scaled year-independent combinational rules provide an easy-to-apply and cost-effective machine-learning expert system for early detection of hidden mastitis using milking parameters.

Environmental Factors and Genetic Parameters of Beef Traits in Fleckvieh Cattle Using Field and Station Testing

Simple Summary

Meat production plays an important role in the efficiency of rearing dual-purpose breeds of cattle in Europe. It is generally known and accepted that most of the beef produced in the EU is produced from dairy herds. This fact provided the motivation for analyzing the influence of genetic and environmental factors on the indicators of meat production from Fleckvieh (known also as dairy Simmental) cattle in the Czech Republic and the interrelationships between these factors. The analysis included evaluating the statistical relationship between the degree of fatness, the heritability and genetic correlation with other traits and also the relationship with exterior muscularity. Together with this analysis, the normal indicators, which are usually evaluated in the population, were also analyzed. The results of our study indicate that meat yield traits are positively related to the live weight of the animal. These findings create possibilities for application in breeding strategies of dual-purpose cattle.

Abstract

The goal of this study was to analyze the genetic and environmental factors of selected meat yield indicators in Fleckvieh cattle in the Czech Republic, through the application of station (S) and field (F) testing methods. Data collected from fattened bulls were analyzed for F (n = 9378) and for S (n = 6346). In the F method and the S method, the values of the main meat yield indicators were as follows: carcass weight 402.91 kg (F), 339.37 kg (S); carcass daily gain 626.05 g/day (F), 609.74 g/day (S); SEUROP carcass classification 2.73 (F), 2.19 (S). Environmental factors were found to have a significant impact on the selected meat yield indicators; their heritability ranged from 0.14 (SEUROP classification) to 0.33 (dressing percentage). The genetic trend was significantly positive only in relation to those meat yield traits, which had a positive link to the size or weight of the animal. The genetic correlation between observations obtained in the S and F methods of testing was very high in relation to the carcass daily gain (0.8351) and carcass weight (0.8244), while slightly lower correlations were calculated for the SEUROP classification. A genetic evaluation of the degree of fatness is not routinely performed in Fleckvieh populations, and the newly established heritability for this trait ranges between 0.17–0.20. The genetic correlation between beef yield indicators and the exterior trait of muscularity was also established, and shows a strong link to the net daily gain, the SEUROP classification and body weight (0.79–0.97). The aim of the study was to evaluate the genetic and environmental effects on meat yield and also estimate genetic parameters for new traits. We can also state, based on the results, that a strong positive genetic trend is confirmed, especially in traits related to the size or weight of animals. This result can be used in breeding programs of dual-purpose cattle, where we can genetically improve the meat and milk yield through the body size.

Genome-wide association study for genotype by lactation stage interaction of milk production traits in dairy cattle

Substantial evidence demonstrates that the genetic background of milk production traits changes during lactation. However, most GWAS for milk production traits assume that genetic effects are constant during lactation and therefore might miss those quantitative trait loci (QTL) whose effects change during lactation. The GWAS for genotype by lactation stage interaction are aimed at explicitly detecting the QTL whose effects change during lactation. The purpose of this study was to perform GWAS for genotype by lactation stage interaction for milk yield, lactose yield, lactose content, fat yield, fat content, protein yield, and somatic cell score to detect QTL with changing effects during lactation. For this study, 19,286 test-day records of 1,800 first-parity Dutch Holstein cows were available and cows were genotyped using a 50K SNP panel. A total of 7 genomic regions with effects that change during lactation were detected in the GWAS for genotype by lactation stage interaction. Two regions on Bos taurus autosome (BTA)14 and BTA19 were also significant based on a GWAS that assumed constant genetic effects during lactation. Five regions on BTA4, BTA10, BTA11, BTA16, and BTA23 were only significant in the GWAS for genotype by lactation stage interaction. The biological mechanisms that cause these changes in genetic effects are still unknown, but negative energy balance and effects of pregnancy may play a role. These findings increase our understanding of the genetic background of lactation and may contribute to the development of better management indicators based on milk composition.