A genome-wide association study for clinical mastitis in first parity US Holstein cows using single-step approach and genomic matrix re-weighting procedure

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/or 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 known to affect the udder synthesis ability too, 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 also subclinical mastitis 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.

Genome-Wide Association Study for Milk Somatic Cell Score in Holstein Friesian Cows in Slovenia

Mastitis is a serious challenge for the dairy industry, leading to economic losses and affecting milk quality. The aim of this study is to identify genetic factors associated with mastitis resistance by conducting a genome-wide association study (GWAS) for the somatic cell score (SCS). Phenotypic records of 350 Holstein Friesian cows were obtained from the Slovenian Cattle Recording Scheme Database and consisted of around 1500 lactation data from 2012 to 2023 collected on a single farm in Slovenia. Corresponding genotypic data were also retrieved from the same database and genotyped using the Illumina BovineSNP50 BeadChip (Illumina, Inc., San Diego, CA, USA). For the association study, three SCS parameters were considered, including lactation mean somatic cell score (LM_SCS), maximum SCS value (SCSMAX), and top three mean value of SCS (TOP3). After performing a GWAS using FarmCPU and BLINK models, five significant SNPs associated with the TOP3 trait were found on BTA 14, 15, 22, and 29. The identified SNP markers were closely linked to six known candidate genes (DNASE1L3, SLC36A4, ARMC1, PDE7A, MMP13, CD44). These results indicate potential genetic markers associated with SCS in the Slovenian Holstein Friesian population.

Identification of candidate genes associated with milk production and mastitis based on transcriptome-wide association study

Genetic research for the assessment of mastitis and milk production traits simultaneously has a long history. The main issue that arises in this context is the known existence of a positive correlation between the risk of mastitis and lactation performance due to selection. The transcriptome-wide association study (TWAS) approach endeavors to combine the expression quantitative trait loci and genome-wide association study summary statistics to decode complex traits or diseases. Accordingly, we used the farmgtex project results as a complete bovine database for mastitis and milk production. The results of colocalization and TWAS approaches were used for the detection of functional associated candidate genes with milk production and mastitis traits on multiple tissue-based transcriptome records. Also, we used the david database for gene ontology to identify significant terms and associated genes. For the identification of interaction networks, the genemania and string databases were used. Also, the available z-scores in TWAS results were used for the calculation of the correlation between tissues. Therefore, the present results confirm that LYNX1, DGAT1, C14H8orf33, and LY6E were identified as significant genes associated with milk production in eight, six, five, and five tissues, respectively. Also, FBXL6 was detected as a significant gene associated with mastitis trait. CLN3 and ZNF34 genes emerged via both the colocalization and TWAS approaches as significant genes for milk production trait. It is expected that TWAS and colocalization can improve our perception of the potential health status control mechanism in high-yielding dairy cows.

Genetic factors underlying host resistance to Rhipicephalus microplus tick infestation in Braford cattle: a systems biology perspective

Approximately 80% of the world's cattle are raised in regions with a high risk of tick-borne diseases, resulting in significant economic losses due to parasitism by Rhipicephalus (Boophilus) microplus. However, the lack of a systemic biology approach hampers a comprehensive understanding of tick-host interactions that mediate tick resistance phenotypes. Here, we conducted a genome-wide association study (GWAS) of 2933 Braford cattle and found 340 single-nucleotide polymorphisms (SNPs) associated with tick counts. Gene expression analyses were performed on skin samples obtained from previously tick-exposed heifers with extremely high or low estimated breeding values for R. microplus counts. Evaluations were performed both before and after artificial infestation with ticks. Differentially expressed genes were found within 1-Mb windows centered at significant SNPs from GWAS. A total of 330 genes were related to the breakdown of homeostasis that was induced by larval attachment to bovine skin. Enrichment analysis pointed to a key role of proteolysis and signal transduction via JAK/STAT, NFKB and WNT/beta catenin signaling pathways. Integrative analysis on matrixEQTL revealed two cis-eQTLs and four significant SNPs in the genes peptidyl arginine deiminase type IV (PADI4) and LOC11449251. The integration of genomic data from QTL maps and transcriptome analyses has identified a set of twelve key genes that show significant associations with tick loads. These genes could be key candidates to improve the accuracy of genomic predictions for tick resistance in Braford cattle.

Inheritance of genomic regions and genes associated with number of oocytes and embryos in Gir cattle through daughter design

Over the past decades, daughter designs, including genotyped sires and their genotyped daughters, have been used as an approach to identify QTL related to economic traits. The aim of this study was to identify genomic regions inherited by Gir sire families and genes associated with number of viable oocytes (VO), total number of oocytes (TO), and number of embryos (EMBR) based on a daughter design approach. In total, 15 Gir sire families were selected. The number of daughters per family ranged from 26 to 395, which were genotyped with different SNP panels and imputed to the Illumina BovineHD BeadChip (777K) and had phenotypes for oocyte and embryo production. Daughters had phenotypic data for VO, TO, and EMBR. The search for QTL was performed through GWAS based on GBLUP. The QTL were found for each trait among and within families based on the top 10 genomic windows with the greatest genetic variance. For EMBR, genomic windows identified among families were located on BTA4, BTA5, BTA6, BTA7, BTA8, BTA13, BTA16, and BTA17, and they were most frequent on BTA7 within families. For VO, genomic windows were located on BTA2, BTA4, BTA5, BTA7, BTA17, BTA21, BTA22, BTA23, and BTA27 among families, being most frequent on BTA8 within families. For TO, the top 10 genomic windows were identified on BTA2, BTA4, BTA5, BTA7, BTA17, BTA21, BTA22, BTA26, and BTA27, being most frequent on BTA7 and BTA8 within families. Considering all results, the greatest number of genomic windows was found on BTA7, where the VCAN, XRCC4, TRNAC-ACA, HAPLN1, and EDIL3 genes were identified in the common regions. In conclusion, 15 Gir sire families with 26 to 395 daughters per family with phenotypes for oocyte and embryo production helped to identify the inheritance of several genomic regions, especially on BTA7, where the EDIL3, HAPLN1, and VCAN candidate genes were associated with number of oocytes and embryos in Gir cattle families.

Copy number variant scan in more than four thousand Holstein cows bred in Lombardy, Italy

Copy Number Variants (CNV) are modifications affecting the genome sequence of DNA, for instance, they can be duplications or deletions of a considerable number of base pairs (i.e., greater than 1000 bp and up to millions of bp). Their impact on the variation of the phenotypic traits has been widely demonstrated. In addition, CNVs are a class of markers useful to identify the genetic biodiversity among populations related to adaptation to the environment. The aim of this study was to detect CNVs in more than four thousand Holstein cows, using information derived by a genotyping done with the GGP (GeneSeek Genomic Profiler) bovine 100K SNP chip. To detect CNV the SVS 8.9 software was used, then CNV regions (CNVRs) were detected. A total of 123,814 CNVs (4,150 non redundant) were called and aggregated into 1,397 CNVRs. The PCA results obtained using the CNVs information, showed that there is some variability among animals. For many genes annotated within the CNVRs, the role in immune response is well known, as well as their association with important and economic traits object of selection in Holstein, such as milk production and quality, udder conformation and body morphology. Comparison with reference revealed unique CNVRs of the Holstein breed, and others in common with Jersey and Brown. The information regarding CNVs represents a valuable resource to understand how this class of markers may improve the accuracy in prediction of genomic value, nowadays solely based on SNPs markers.

Single-step genome-wide association analyses for milk urea concentration in Walloon Holstein cows

The aims of this study were to conduct a single-step genome-wide association to identify genomic regions associated with milk urea (MU) and to estimate genetic correlations between MU and milk yield (MY), milk composition (calcium content [CC], fat percentage [FP], protein percentage [PP], and casein percentage [CNP]), and cheese-making properties (CMP; coagulation time [CT], curd firmness after 30 min from rennet addition [a30], and titratable acidity [TA]). The used data have been collected from 2015 to 2020 on 78,073 first-parity (485,218 test-day records) and 48,766 second-parity (284,942 test-day records) Holstein cows distributed in 671 herds in the Walloon Region of Belgium. Data of 565,533 SNP located on 29 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 the top-3 genomic regions explaining the largest rate of the genetic variance were considered promising regions and used to identify potential candidate genes. Mean (SD) MU was 25.38 (8.02) mg/dL and 25.03 (8.06) mg/dL in the first and second lactation, respectively. Mean heritability estimates for daily MU were 0.21 and 0.23 for the first and second lactation, respectively. The genetic correlations estimated between MU and MY, milk composition, and CMP were quite low (ranged from -0.10 [CC] to 0.10 [TA] and -0.05 [CT] to 0.13 [TA] for the first and second lactations, respectively). The top-3 regions associated with MU were located from 80.61 to 80.74 Mb on BTA6, 103.26 to 103.41 Mb on BTA11, and 1.59 to 2.15 Mb on BTA14. Genes including KCNT1, MROH1, SHARPIN, TSSK5, CPSF1, HSF1, TONSL, DGAT1, SCX, and MAF1 were identified as positional candidate genes for MU. The findings of this study can be used for a better understanding of the genomic architecture underlying MU in Holstein cattle.

Multi-omics integration identifies regulatory factors underlying bovine subclinical mastitis

BACKGROUND

Mastitis caused by multiple factors remains one of the most common and costly disease of the dairy industry. Multi-omics approaches enable the comprehensive investigation of the complex interactions between multiple layers of information to provide a more holistic view of disease pathogenesis. Therefore, this study investigated the genomic and epigenomic signatures and the possible regulatory mechanisms underlying subclinical mastitis by integrating RNA sequencing data (mRNA and lncRNA), small RNA sequencing data (miRNA) and DNA methylation sequencing data of milk somatic cells from 10 healthy cows and 20 cows with naturally occurring subclinical mastitis caused by Staphylococcus aureus or Staphylococcus chromogenes.

RESULTS

Functional investigation of the data sets through gene set analysis uncovered 3458 biological process GO terms and 170 KEGG pathways with altered activities during subclinical mastitis, provided further insights into subclinical mastitis and revealed the involvement of multi-omics signatures in the altered immune responses and impaired mammary gland productivity during subclinical mastitis. The abundant genomic and epigenomic signatures with significant alterations related to subclinical mastitis were observed, including 30,846, 2552, 1276 and 57 differential methylation haplotype blocks (dMHBs), differentially expressed genes (DEGs), lncRNAs (DELs) and miRNAs (DEMs), respectively. Next, 5 factors presenting the principal variation of differential multi-omics signatures were identified. The important roles of Factor 1 (DEG, DEM and DEL) and Factor 2 (dMHB and DEM), in the regulation of immune defense and impaired mammary gland functions during subclinical mastitis were revealed. Each of the omics within Factors 1 and 2 explained about 20% of the source of variation in subclinical mastitis. Also, networks of important functional gene sets with the involvement of multi-omics signatures were demonstrated, which contributed to a comprehensive view of the possible regulatory mechanisms underlying subclinical mastitis. Furthermore, multi-omics integration enabled the association of the epigenomic regulatory factors (dMHBs, DELs and DEMs) of altered genes in important pathways, such as 'Staphylococcus aureus infection pathway' and 'natural killer cell mediated cytotoxicity pathway', etc., which provides further insights into mastitis regulatory mechanisms. Moreover, few multi-omics signatures (14 dMHBs, 25 DEGs, 18 DELs and 5 DEMs) were identified as candidate discriminant signatures with capacity of distinguishing subclinical mastitis cows from healthy cows.

CONCLUSION

The integration of genomic and epigenomic data by multi-omics approaches in this study provided a better understanding of the molecular mechanisms underlying subclinical mastitis and identified multi-omics candidate discriminant signatures for subclinical mastitis, which may ultimately lead to the development of more effective mastitis control and management strategies.

Genome-Wide Analysis of Milk Production Traits and Selection Signatures in Serbian Holstein-Friesian Cattle

To improve the genomic evaluation of milk-related traits in Holstein-Friesian (HF) cattle it is essential to identify the associated candidate genes. Novel SNP-based analyses, such as the genetic mapping of inherited diseases, GWAS, and genomic selection, have led to a new era of research. The aim of this study was to analyze the association of each individual SNP in Serbian HF cattle with milk production traits and inbreeding levels. The SNP 60 K chip Axiom Bovine BovMDv3 was deployed for the genotyping of 334 HF cows. The obtained genomic results, together with the collected phenotypic data, were used for a GWAS. Moreover, the identification of ROH segments was performed and served for inbreeding coefficient evaluation and ROH island detection. Using a GWAS, a polymorphism, rs110619097 (located in the intron of the CTNNA3 gene), was detected to be significantly (p < 0.01) associated with the milk protein concentration in the first lactation (adjusted to 305 days). The average genomic inbreeding value (FROH) was 0.079. ROH islands were discovered in proximity to genes associated with milk production traits and genomic regions under selection pressure for other economically important traits of dairy cattle. The findings of this pilot study provide useful information for a better understanding of the genetic architecture of milk production traits in Serbian HF dairy cows and can be used to improve lactation performances in Serbian HF cattle breeding programs.

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.

Single-step genome-wide association for selected milk fatty acids in Dual-Purpose Belgian Blue cows

The aim of this study was to estimate genetic parameters and identify genomic regions associated with selected individual and groups of milk fatty acids (FA) predicted by milk mid-infrared spectrometry in Dual-Purpose Belgian Blue cows. The used data were 69,349 test-day records of milk yield, fat percentage, and protein percentage along with selected individual and groups FA of milk (g/dL milk) collected from 2007 to 2020 on 7,392 first-parity (40,903 test-day records), and 5,185 second-parity (28,446 test-day records) cows distributed in 104 herds in the Walloon Region of Belgium. Data of 28,466 SNPs, located on 29 Bos taurus autosomes (BTA), of 1,699 animals (639 males and 1,060 females) 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 best linear unbiased prediction approach. The proportion of genetic variance explained by each 25-SNP sliding window (with an average size of ~2 Mb) was calculated, and regions accounting for at least 1.0% of the total additive genetic variance were used to search for candidate genes. Average daily heritability estimated for the included milk FA traits ranged from 0.01 (C4:0) to 0.48 (C12:0) and 0.01 (C4:0) to 0.42 (C12:0) in the first and second parities, respectively. Genetic correlations found between milk yield and the studied individual milk FA, except for C18:0, C18:1 trans, C18:1 cis-9, were positive. The results showed that fat percentage and protein percentage were positively genetically correlated with all studied individual milk FA. Genome-wide association analyses identified 11 genomic regions distributed over 8 chromosomes [BTA1, BTA4, BTA10, BTA14 (4 regions), BTA19, BTA22, BTA24, and BTA26] associated with the studied FA traits, though those found on BTA14 partly overlapped. The genomic regions identified differed between parities and lactation stages. Although these differences in genomic regions detected may be due to the power of quantitative trait locus detection, it also suggests that candidate genes underlie the phenotypic expression of the studied traits may vary between parities and lactation stages. These findings increase our understanding about the genetic background of milk FA and can be used for the future implementation of genomic evaluation to improve milk FA profile in Dual-Purpose Belgian Blue cows.

Single-step genome-wide association studies and post-GWAS analyses for the number of oocytes and embryos in Gir cattle

Genome-Wide Association Studies (GWAS) are used for identification of quantitate trait loci (QTL) and genes associated with several traits. We aimed to identify genomic regions, genes, and biological processes associated with number of total and viable oocytes, and number of embryos in Gir dairy cattle. A dataset with 17,526 follicular aspirations, including the following traits: number of viable oocytes (VO), number of total oocytes (TO), and number of embryos (EMBR) from 1641 Gir donors was provided by five different stock farms. A genotype file with 2093 animals and 395,524 SNP markers was used to perform a single-step GWAS analysis for each trait. The top 10 windows with the highest percentage of additive genetic variance explained by 100 adjacent SNPs were selected. The genomic regions identified in our work were overlapped with QTLs from QTL database on chromosomes 1, 2, 5, 6, 7, 8, 9, 13, 17, 18, 20, 21, 22, 24, and 29. These QTLs were classified as External, Health, Meat and carcass, Production or Reproduction traits, and about 38% were related to Reproduction. In total, 117 genes were identified, of which 111 were protein-coding genes. Exclusively associations were observed for 42 genes with EMBR, and 1 with TO. Also, 42 genes were in common between VO and TO, 28 between VO and EMBR, and four genes were in common among all traits. In conclusion, great part of the identified genes plays a functional role in initial embryo development or general cell functions. The protein-coding genes ARNT, EGR1, HIF1A, AHR, and PAX2 are good markers for the production of oocytes and embryos in Gir cattle.

Invited review: Good practices in genome-wide association studies to identify candidate sequence variants in dairy cattle

Genotype data from dairy cattle selection programs have greatly facilitated GWAS to identify variants related to economic traits. Results can enhance the accuracy of genomic prediction, analyze more complex models that go beyond additive effects, elucidate the genetic architecture of a trait, and finally, decipher the underlying biology of traits. The entire process, comprising data generation, quality control, statistical analyses, interpretation of association results, and linking results to biology should be designed and executed to minimize the generation of false-positive and false-negative associations and misleading links to biological processes. This review aims to provide general guidelines for data analysis that address data quality control, association tests, adjustment for population stratification, and significance evaluation to improve the reliability of conclusions. We also provide guidance on post-GWAS strategy and the interpretation of results. These guidelines are tailored to dairy cattle, which are characterized by long-range linkage disequilibrium, large half-sib families, and routinely collected phenotypes, requiring different approaches than those applied in human GWAS. We discuss common limitations and challenges that have been overlooked in the analysis and interpretation of GWAS to identify candidate sequence variants in dairy cattle.

Assessment of Genetic and Health Management of Tunisian Holstein Dairy Herds with a Focus on Longevity

In Tunisia, the recognition of the possibility of including longevity and disease resistance in dairy cattle selection objectives has been hypothesized as a useful strategy by both researchers and producers. However, in this paper, the state of the art, with a focus on health and longevity, is reviewed. Along the same lines, the heritability for the milk traits, fertility traits, and longevity of Tunisian Holstein dairy cows complies with the literature. Therefore, the influence of genetics on some diseases of the dairy cow was investigated. In addition, a decreasing efficiency in cow fertility has been observed over the last few years. The results showed that the risk of culling increased with common diseases. When analyzed with the Weibull model, functional lifespan was strongly influenced by milk yield; therefore, the risk increased with a reduced milk yield. In her first three lactations, the relative risk of selection increased gradually with lactation. Thus, the risk of thinning is highest at the beginning and end of the first feeding and the end of her second feeding. In conclusion, the risk of culling was reduced in parity. The factors that influence the life of the herd, such as health, husbandry, environmental conditions, and management, are often ignored when evaluating longevity.

Candidate genes for mastitis resistance in dairy cattle: a data integration approach

BACKGROUND

Inflammation of the mammary tissue (mastitis) is one of the most detrimental health conditions in dairy ruminants and is considered the most economically important infectious disease of the dairy sector. Improving mastitis resistance is becoming an important goal in dairy ruminant breeding programmes. However, mastitis resistance is a complex trait and identification of mastitis-associated alleles in livestock is difficult. Currently, the only applicable approach to identify candidate loci for complex traits in large farm animals is to combine different information that supports the functionality of the identified genomic regions with respect to a complex trait.

METHODS

To identify the most promising candidate loci for mastitis resistance we integrated heterogeneous data from multiple sources and compiled the information into a comprehensive database of mastitis-associated candidate loci. Mastitis-associated candidate genes reported in association, expression, and mouse model studies were collected by searching the relevant literature and databases. The collected data were integrated into a single database, screened for overlaps, and used for gene set enrichment analysis.

RESULTS

The database contains candidate genes from association and expression studies and relevant transgenic mouse models. The 2448 collected candidate loci are evenly distributed across bovine chromosomes. Data integration and analysis revealed overlaps between different studies and/or with mastitis-associated QTL, revealing promising candidate genes for mastitis resistance.

CONCLUSION

Mastitis resistance is a complex trait influenced by numerous alleles. Based on the number of independent studies, we were able to prioritise candidate genes and propose a list of the 22 most promising. To our knowledge this is the most comprehensive database of mastitis associated candidate genes and could be helpful in selecting genes for functional validation studies.

Identification of New Candidate Genes Related to Semen Traits in Duroc Pigs through Weighted Single-Step GWAS

Semen traits play a key role in the pig industry because boar semen is widely used in purebred and crossbred pigs. The production of high-quality semen is crucial to ensuring a good result in artificial insemination. With the wide application of artificial insemination in the pig industry, more and more attention has been paid to the improvement of semen traits by genetic selection. The purpose of this study was to identify the genetic regions and candidate genes associated with semen traits of Duroc boars. We used weighted single-step GWAS to identify candidate genes associated with sperm motility, sperm progressive motility, sperm abnormality rate and total sperm count in Duroc pigs. In Duroc pigs, the three most important windows for sperm motility-sperm progressive motility, sperm abnormality rate, and total sperm count-explained 12.45%, 9.77%, 15.80%, and 12.15% of the genetic variance, respectively. Some genes that are reported to be associated with spermatogenesis, testicular function and male fertility in mammals have been detected previously. The candidate genes CATSPER1, STRA8, ZSWIM7, TEKT3, UBB, PTBP2, EIF2B2, MLH3, and CCDC70 were associated with semen traits in Duroc pigs. We found a common candidate gene, STRA8, in sperm motility and sperm progressive motility, and common candidate genes ZSWIM7, TEKT3 and UBB in sperm motility and sperm abnormality rate, which confirms the hypothesis of gene pleiotropy. Gene network enrichment analysis showed that STRA8, UBB and CATSPER1 were enriched in the common biological process and participated in male meiosis and spermatogenesis. The SNPs of candidate genes can be given more weight in genome selection to improve the ability of genome prediction. This study provides further insight into the understanding the genetic structure of semen traits in Duroc boars.

New Insights on Nucleotide Sequence Variants and mRNA Levels of Candidate Genes Assessing Resistance/Susceptibility to Mastitis in Holstein and Montbéliarde Dairy Cows

A major factor in the propagation of an infectious disease is host genetics. In this study, 180 dairy cows (90 of each breed: Holstein and Montbéliarde) were used. Each breed's tested dairy cows were divided into two groups of comparable size (45 cows each), mastitis-free and mastitis-affected groups. Each cow's jugular vein was punctured to obtain blood samples for DNA and RNA extraction. In the examined Holstein and Montbéliarde dairy cows, single nucleotide polymorphisms (SNPs) related with mastitis resistance/susceptibility were found in the RASGRP1, NFkB, CHL1, MARCH3, PDGFD, MAST3, EPS15L1, C1QTNF3, CD46, COX18, NEURL1, PPIE, and PTX3 genes. Chi-square analysis of identified SNPs revealed a significant difference in gene frequency between mastitic and healthy cows. Except for CHL1, mastitic dairy cows of two breeds had considerably higher mRNA levels of the examined genes than did healthy ones. Marker-assisted selection and monitoring of dairy cows' susceptibility to mastitis may be accomplished through the use of discovered SNPs and changes in the gene expression profile of the studied genes. These findings also point to a possible method for reducing mastitis in dairy cows through selective breeding of animals using genetic markers linked to an animal's ability to resist infection.

Genome-wide scan for selective sweeps identifies novel loci associated with resistance to mastitis in German Holstein cattle

Domestication and selection significantly changed phenotypic and behavioural traits in modern domestic animals. In this study, to identify the genomic regions associated with mastitis, genomic data of German Holstein dairy cattle were analysed. The samples were genotyped using the Bovine 50 K SNP chip. For each defined healthy and sick group, 133 samples from 13,276 genotyped dairy cows were selected based on mastitis random residual effects. Grouping was done to infer selection signatures based on XP-EHH statistic. The results revealed that for the top 0.01 percentile of the obtained XP-EHH values, five genomic regions on chromosomes 8, 11, 12, 14 and 26 of the control group, and four regions on chromosomes 3, 4 (two regions) and 22 of the case group, have been under selection. Also, consideration of the top 0.1 percentile of the XP-EHH values, clarified 21 and 15 selective sweeps in the control and case group, respectively. This study identified some genomic regions containing potential candidate genes associated with resistance and susceptibility to mastitis, immune system and inflammation, milk traits, udder morphology and different types of cancers. In addition, these regions overlap with some quantitative trait loci linked to clinical mastitis, immunoglobulin levels, somatic cell score, udder traits, milk fat and protein, milk yield, milking speed and veterinary treatments. It is noteworthy that we found two regions in the healthy group (on chromosomes 12 and 14) with strong signals, which were not described previously. It is likely that future research could link these identified genomic regions to mastitis. The results of the current study contribute to the identification of causal mutations, genomic regions and genes affecting mastitis incidence in dairy cows.

Can breeders prevent pathogen adaptation when selecting for increased resistance to infectious diseases?

Background

Recent research shows that genetic selection has high potential to reduce the prevalence of infectious diseases in livestock. However, like all interventions that target infectious diseases, genetic selection of livestock can exert selection pressure on pathogen populations. Such selection on the pathogen may lead to escape strategies and reduce the effect of selection of livestock for disease resistance. Thus, to successfully breed livestock for lower disease prevalence, it is essential to develop strategies that prevent the invasion of pathogen mutants that escape host resistance. Here we investigate the conditions under which such “escape mutants” can replace wild-type pathogens in a closed livestock population using a mathematical model of disease transmission.

Results

Assuming a single gene that confers sufficient resistance, results show that genetic selection for resistance in livestock typically leads to an “invasion window” within which an escape mutant of the pathogen can invade. The bounds of the invasion window are determined by the frequency of resistant hosts in the population. The lower bound occurs when the escape mutant has an advantage over the wild-type pathogen in the population. The upper bound occurs when local eradication of the pathogen is expected. The invasion window is smallest when host resistance is strong and when infection with the wild-type pathogen provides cross immunity to infection with the escape mutant.

Conclusions

To minimise opportunities for pathogens to adapt, under the assumptions of our model, the aim of disease control through genetic selection should be to achieve herd-level eradication of the infection faster than the rate of emergence of escape mutants of the pathogen. Especially for microparasitic infections, this could be achieved by placing animals into herds according to their genetic resistance, such that these herds stay completely out of the invasion window. In contrast to classical breeding theory, our model suggests that multi-trait selection with gradual improvement of each trait of the breeding goal might not be the best strategy when resistance to infectious disease is part of the breeding goal. Temporally, combining genetic selection with other interventions helps to make the invasion window smaller, and thereby reduces the risk of invasion of escape mutants.

Immune mechanisms, resistance genes, and their roles in the prevention of mastitis in dairy cows

Mastitis is one of the most important diseases of the mammary gland. The increased incidence of this disease in cows is due to the breeding of dairy cattle for higher yields, which is accompanied by an increased susceptibility to mastitis. Therefore, the difficulty involved with preventing this disease has increased. An integral part of current research is the elimination of mastitis in order to reduce the consumption of antibiotic drugs, thereby reducing the resistance of microorganisms and decreasing companies' economic losses due to mastitis (i.e. decreased milk yield, increased drug costs, and reduced milk supply). Susceptibility to mastitis is based on dairy cows' immunity, health, nutrition, and welfare. Thus, it is important to understand the immune processes in the body in order to increase the resistance of animals. Recently, various studies have focused on the selection of mastitis resistance genes. An important point is also the prevention of mastitis. This publication aims to describe the physiology of the mammary gland along with its immune mechanisms and to approximate their connection with potential mastitis resistance genes. This work describes various options for mastitis elimination and focuses on genetic selection and a closer specification of resistance genes to mastitis. Among the most promising resistance genes for mastitis, we consider CD14, CXCR1, lactoferrin, and lactoglobulin.

Identifying genetic variants and pathways influencing daughter averages for twinning in North American Holstein cattle and evaluating the potential for genomic selection

Multiple birth in dairy cattle is a detrimental trait both economically for producers and for animal health. Genetics of twinning is complex and has led to several quantitative trait loci regions being associated with increased twinning. To identify variants associated with this trait, calving records from 2 time periods were used to estimate daughter averages for twinning for Holstein bulls. Multiple analyses were conducted and compared including GWAS, genomic prediction, and gene set enrichment analysis for pathway detection. Although pathway analysis did not yield many congruent pathways of interest between data sets, it did indicate two of interest. Both pathways have ties to the strong candidate region on BTA11 from the genome-wide association analysis across data sets. This region does not overlap with previously identified quantitative trait loci regions for twinning or ovulation rate in cattle. The strongest associated SNPs were upstream from 2 candidate genes LHCGR and FSHR, which are involved in folliculogenesis. Genomic prediction showed a moderate correlation accuracy (0.43) when predicting genomic breeding values for bulls with estimates from calving records from 2010 to 2016. Future analysis of the region on BTA11 and the relation of the candidate genes could improve this accuracy.

Genome-wide associations for heat stress response suggest potential candidate genes underlying milk fatty acid composition in dairy cattle

Contents of milk fatty acids (FA) display remarkable alterations along climatic gradients. Detecting candidate genes underlying such alterations might be beneficial for the exploration of climate sensitivity in dairy cattle. Consequently, we aimed on the definition of FA heat stress indicators, considering FA breeding values in response to temperature-humidity index (THI) alterations. Indicators were used in GWAS, in ongoing gene annotations and for the estimation of chromosome-wide variance components. The phenotypic data set consisted of 39,600 test-day milk FA records from 5,757 first-lactation Holstein dairy cows kept in 16 large-scale German cooperator herds. The FA traits were C18:0, polyunsaturated fatty acids (PUFA), saturated fatty acids (SFA), and unsaturated fatty acids (UFA). After genotype quality control, 40,523 SNP markers from 3,266 cows and 930 sires were considered. Meteorological data from the weather station in closest herd distance were used for the calculation of maximum hourly daily THI, which were allocated to 10 different THI classes. The same FA from 3 stages of lactation were considered as different, but genetically correlated traits. Consequently, a 3-trait reaction norm model was used to estimate genetic parameters and breeding values for FA along THI classes, considering either pedigree (A) or genomic (G) relationship matrices. De-regressed proofs and genomic estimated breeding values at the intermediate THI class 5 and at the extreme THI class 10 were used as pseudophenotypes in ongoing genomic analyses for thermoneutral (TNC) and heat stress conditions (HSC), respectively. The differences in de-regressed proofs and in genomic estimated breeding values from both THI classes were pseudophenotypes for heat stress response (HSR). Genetic correlations between the same FA under TNC and HSC were smallest in the first lactation stage and ranged from 0.20 for PUFA to 0.87 for SFA when modeling with the A matrix, and from 0.35 for UFA to 0.86 for SFA when modeling with the G matrix. In the first lactation stage, larger additive genetic variances under HSC compared with TNC indicate climate sensitivity for C18:0, PUFA, and UFA. Climate sensitivity was also reflected by pronounced chromosome-wide genetic variances for HSR of PUFA and UFA in the first stage of lactation. For all FA under TNC, HSC, and HSR, quite large genetic variance proportions were explained by BTA14. In GWAS, 30 SNP (within or close to 38 potential candidate genes) overlapped for HSR of the different FA. One unique potential candidate gene (AMFR) was detected for HSR of PUFA, 15 for HSR of SFA (ADGRB1, DENND3, DUSP16, EFR3A, EMP1, ENSBTAG00000003838, EPS8, MGP, PIK3C2G, STYK1, TMEM71, GSG1, SMARCE1, CCDC57, and FASN) and 3 for HSR of UFA (ENSBTAG00000048091, PAEP, and EPPK1). The identified unique genes play key roles in milk FA synthesis and are associated with disease resistance in dairy cattle. The results suggest consideration of FA in combination with climatic responses when inferring genetic mechanisms of heat stress in dairy cows.

Genetic Variability of Methane Production and Concentration Measured in the Breath of Polish Holstein-Friesian Cattle

The genetic architecture of methane (CH₄) production remains largely unknown. We aimed to estimate its heritability and to perform genome-wide association studies (GWAS) for the identification of candidate genes associated with two phenotypes: CH₄ in parts per million/day (CH₄ ppm/d) and CH₄ in grams/day (CH₄ g/d). We studied 483 Polish Holstein-Friesian cows kept on two commercial farms in Poland. Measurements of CH₄ and carbon dioxide (CO₂) concentrations exhaled by cows during milking were obtained using gas analyzers installed in the automated milking system on the farms. Genomic analyses were performed using a single-step BLUP approach. The percentage of genetic variance explained by SNPs was calculated for each SNP separately and then for the windows of neighbouring SNPs. The heritability of CH₄ ppm/d ranged from 0 to 0.14, with an average of 0.085. The heritability of CH₄ g/d ranged from 0.13 to 0.26, with an average of 0.22. The GWAS detected potential candidate SNPs on BTA 14 which explained ~0.9% of genetic variance for CH₄ ppm/d and ~1% of genetic variance for CH₄ g/d. All identified SNPs were located in the TRPS1 gene. We showed that methane traits are partially controlled by genes; however, the detected SNPs explained only a small part of genetic variation-implying that both CH₄ ppm/d and CH₄ g/d are highly polygenic traits.

Understanding the genomic architecture of clinical mastitis in Bos indicus

This study elucidated potential genetic variants and QTLs associated with clinical mastitis incidence traits in Bos indicus breed, Sahiwal. Estimated breeding values for the traits (calculated using Bayesian inference) were used as pseudo-phenotypes for association with genome-wide SNPs and further QTL regions underlying the traits were identified. In all, 25 SNPs were found to be associated with the traits at the genome-wide suggestive threshold (p ≤ 5 × 10^-4) and these SNPs were used to define QTL boundaries based on the linkage disequilibrium structure. A total of 16 QTLs were associated with the trait EBVs including seven each for clinical mastitis incidence (CMI) in first and second lactations and two for CMI in third lactation. Nine out of sixteen QTLs overlapped with the already reported QTLs for mastitis traits, whereas seven were adjudged as novel ones. Important candidates for clinical mastitis in the identified QTL regions included DNAJB9, ELMO1, ARHGAP26, NR3C1, CACNB2, RAB4A, GRB2, NUP85, SUMO2, RBPJ, and RAB33B genes. These findings shed light on the genetic architecture of the disease in Bos indicus, and present potential regions for fine mapping and downstream analysis in future.

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.

Identification of loci associated with susceptibility to Mycobacterium avium subsp. paratuberculosis infection in Holstein cattle using combinations of diagnostic tests and imputed whole-genome sequence data

Bovine paratuberculosis (PTB) is a chronic inflammatory disease caused by Mycobacterium avium susbp. paratuberculosis (MAP). Genome-wide association studies (GWAS) have identified single nucleotide polymorphisms (SNPs) significantly associated with susceptibility to bovine PTB. The main objective of this study was to identify quantitative trait loci (QTLs) associated with MAP infection in Spanish Holstein cows (N = 983) using combinations of diagnostic tests and imputed whole-genome sequence (WGS) data. The infection status of these animals was defined by three diagnostic methods including ELISA for MAP-antibodies detection, and tissue culture and PCR for MAP detection. The 983 cows included in this study were genotyped with the Bovine MD SNP50 Bead Chip, and the corresponding genotypes were imputed to WGS using the 1,000 Bull genomes reference population. In total, 33.77 million SNP variants per animal were identified across the genome. Linear mixed models were used to calculate the heritability (h²) estimates for each diagnostic test and test combinations. Next, we performed a case-control GWAS using the imputed WGS datasets and the phenotypes and combinations of phenotypes with h² estimates > 0.080. After performing the GWAS, the test combinations that showed SNPs with a significant association (P_FDR ≤ 0.05), were the ELISA-tissue PCR-tissue culture, ELISA-tissue culture, and ELISA-tissue PCR. A total of twelve quantitative trait loci (QTLs) highly associated with MAP infection status were identified on the Bos taurus autosomes (BTA) 4, BTA5, BTA11, BTA12, BTA14, BTA23, BTA24, and BTA28, and some of these QTLs were linked to immune-modulating genes. The identified QTLs on BTA23 spanning from 18.81 to 22.95 Mb of the Bos taurus genome overlapped with several QTLs previously found to be associated with PTB susceptibility, bovine tuberculosis susceptibility, and clinical mastitis. The results from this study provide more clues regarding the molecular mechanisms underlying susceptibility to PTB infection in cattle and might be used to develop national genetic evaluations for PTB in Spain.

Genome-Wide Association Study Identifies Candidate Genes Associated with Feet and Leg Conformation Traits in Chinese Holstein Cattle

Simple Summary

Feet and leg problems are among the major reasons for dairy cows leaving the herd, as well as having direct association with production and reproduction efficiency, health (e.g., claw disorders and lameness) and welfare. Hence, understanding the genetic architecture underlying feet and conformation traits in dairy cattle offers new opportunities toward the genetic improvement and long-term selection. Through a genome-wide association study on Chinese Holstein cattle, we identified several candidate genes associated with feet and leg conformation traits. These results could provide useful information about the molecular breeding basis of feet and leg traits, thus improving the longevity and productivity of dairy cattle.

Abstract

Feet and leg conformation traits are considered one of the most important economical traits in dairy cattle and have a great impact on the profitability of milk production. Therefore, identifying the single nucleotide polymorphisms (SNPs), genes and pathways analysis associated with these traits might contribute to the genomic selection and long-term plan selection for dairy cattle. We conducted genome-wide association studies (GWASs) using the fixed and random model circulating probability unification (FarmCPU) method to identify SNPs associated with bone quality, heel depth, rear leg side view and rear leg rear view of Chinese Holstein cows. Phenotypic measurements were collected from 1000 individuals of Chinese Holstein cattle and the GeneSeek Genomic Profiler Bovine 100 K SNP chip was utilized for individual genotyping. After quality control, 984 individual cows and 84,906 SNPs remained for GWAS work; as a result, we identified 20 significant SNPs after Bonferroni correction. Several candidate genes were identified within distances of 200 kb upstream or downstream to the significant SNPs, including ADIPOR2, INPP4A, DNMT3A, ALDH1A2, PCDH7, XKR4 and CADPS. Further bioinformatics analyses showed 34 gene ontology terms and two signaling pathways were significantly enriched (p ≤ 0.05). Many terms and pathways are related to biological quality, metabolism and development processes; these identified SNPs and genes could provide useful information about the genetic architecture of feet and leg traits, thus improving the longevity and productivity of Chinese Holstein dairy cattle.

Genes and Pathways Affecting Sheep Productivity Traits: Genetic Parameters, Genome-Wide Association Mapping, and Pathway Enrichment Analysis

Ewe productivity is a composite and maternal trait that is considered the most important economic trait in sheep meat production. The objective of this study was the application of alternative genome-wide association study (GWAS) approaches followed by gene set enrichment analysis (GSEA) on the ewes’ genome to identify genes affecting pregnancy outcomes and lamb growth after parturition in Iranian Baluchi sheep. Three maternal composite traits at birth and weaning were considered. The traits were progeny birth weight, litter mean weight at birth, total litter weight at birth, progeny weaning weight, litter mean weight at weaning, and total litter weight at weaning. GWASs were performed on original phenotypes as well as on estimated breeding values. The significant SNPs associated with composite traits at birth were located within or near genes RDX, FDX1, ARHGAP20, ZC3H12C, THBS1, and EPG5. Identified genes and pathways have functions related to pregnancy, such as autophagy in the placenta, progesterone production by the placenta, placental formation, calcium ion transport, and maternal immune response. For composite traits at weaning, genes (NR2C1, VEZT, HSD17B4, RSU1, CUBN, VIM, PRLR, and FTH1) and pathways affecting feed intake and food conservation, development of mammary glands cytoskeleton structure, and production of milk components like fatty acids, proteins, and vitamin B-12, were identified. The results show that calcium ion transport during pregnancy and feeding lambs by milk after parturition can have the greatest impact on weight gain as compared to other effects of maternal origin.

Genome-Wide Associations for Microscopic Differential Somatic Cell Count and Specific Mastitis Pathogens in Holstein Cows in Compost-Bedded Pack and Cubicle Farming Systems

Simple Summary

New free walk housing systems such as compost-bedded pack barns might positively influence animal welfare. However, udder health can be negatively affected due to the microbial environment in the pack. Udder health depends on many factors, such as the environment, the feed, the pathogen species, and the genetic mechanisms of the cow’s immune system. For a more precise evaluation of udder health, we examined novel traits including specific mastitis pathogens and differential somatic cell fractions in milk. In order to identify possible candidate genes for udder health, a genome-wide association study, including single-nucleotide polymorphisms (SNP) by housing system interactions (compost-bedded pack barn and conventional cubicle barn), was performed. We identified two potential candidate genes for the interaction effect in relation to udder health. The identified potential candidate gene HEMK1 (HemK methyltransferase family member 1) is involved in immune system development, and CHL1 (cell adhesion molecule L1 like) has an immunosuppressive effect during stress conditions. The results suggest housing system-specific breeding strategies in order to improve udder health in compost-bedded pack and conventional cubicle barns.

Abstract

The aim of the present study was to detect significant SNP (single-nucleotide polymorphism) effects and to annotate potential candidate genes for novel udder health traits in two different farming systems. We focused on specific mastitis pathogens and differential somatic cell fractions from 2198 udder quarters of 537 genotyped Holstein Friesian cows. The farming systems comprised compost-bedded pack and conventional cubicle barns. We developed a computer algorithm for genome-wide association studies allowing the estimation of main SNP effects plus consideration of SNPs by farming system interactions. With regard to the main effect, 35 significant SNPs were detected on 14 different chromosomes for the cell fractions and the pathogens. Six SNPs were significant for the interaction effect with the farming system for most of the udder health traits. We inferred two possible candidate genes based on significant SNP interactions. HEMK1 plays a role in the development of the immune system, depending on environmental stressors. CHL1 is regulated in relation to stress level and influences immune system mechanisms. The significant interactions indicate that gene activity can fluctuate depending on environmental stressors. Phenotypically, the prevalence of mastitis indicators differed between systems, with a notably lower prevalence of minor bacterial indicators in compost systems.

Genome-Wide Association Studies of Somatic Cell Count in the Assaf Breed

Simple Summary

Mastitis causes economic loss due to discarded milk and reduced milk production and quality, increased medical care costs and somatic cell count (SCC) penalties. The use of genetic markers associated with the variability of this trait through marker-assisted selection (MAS) could help traditional methods. Our objectives were to identify new single nucleotide polymorphisms (SNPs) and genes associated with mastitis resistance in Assaf sheep by using the Illumina Ovine Infinium® HD SNP BeadChip (680K). Firstly, corrected phenotype estimates for somatic cell score (SCS) were calculated using 6173 records from 1894 multiparous Assaf ewes, and were used to select 192 extreme animals (low SCS group: n = 96; and high SCS group: n = 96) for the genome-wide association study (GWAS). Four SNPs (rs419096188, rs415580501, rs410336647, and rs424642424), three of them totally linked, were found to be significant at the chromosome level (FDR 10%) in two different regions of OAR19 close to genes related to the immune system response. Validation studies of two SNPs (rs419096188 and rs424642424) by Kompetitive Allele-Specific PCR (KASP) genotyping in the total population (n = 1894) confirmed previous GWAS association results for the SCS trait. Finally, the SNP rs419096188 was also associated with lactose content trait.

Abstract

A genome-wide association study (GWAS) was performed to identify new single nucleotide polymorphisms (SNPs) and genes associated with mastitis resistance in Assaf sheep by using the Illumina Ovine Infinium^® HD SNP BeadChip (680K). In total, 6173 records from 1894 multiparous Assaf ewes with at least three test day records and aged between 2 and 7 years old were used to estimate a corrected phenotype for somatic cell score (SCS). Then, 192 ewes were selected from the top (n = 96) and bottom (n = 96) tails of the corrected SCS phenotype distribution to be used in a GWAS. Although no significant SNPs were found at the genome level, four SNPs (rs419096188, rs415580501, rs410336647, and rs424642424) were significant at the chromosome level (FDR 10%) in two different regions of OAR19. The SNP rs419096188 was located in intron 1 of the NUP210 and close to the HDAC11 genes (61 kb apart), while the other three SNPs were totally linked and located 171 kb apart from the ARPP21 gene. These three genes were related to the immune system response. These results were validated in two SNPs (rs419096188 and rs424642424) in the total population (n = 1894) by Kompetitive Allele-Specific PCR (KASP) genotyping. Furthermore, rs419096188 was also associated with lactose content.

Assessing the genetic background and genomic relatedness of red cattle populations originating from Northern Europe

Background

Local cattle breeds need special attention, as they are valuable reservoirs of genetic diversity. Appropriate breeding decisions and adequate genomic management of numerically smaller populations are required for their conservation. At this point, the analysis of dense genome-wide marker arrays provides encompassing insights into the genomic constitution of livestock populations. We have analyzed the genetic characterization of ten cattle breeds originating from Germany, The Netherlands and Denmark belonging to the group of red dairy breeds in Northern Europe. The results are intended to provide initial evidence on whether joint genomic breeding strategies of these populations will be successful.

Results

Traditional Danish Red and Groningen White-Headed were the most genetically differentiated breeds and their populations showed the highest levels of inbreeding. In contrast, close genetic relationships and shared ancestry were observed for the populations of German Red and White Dual-Purpose, Dutch Meuse-Rhine-Yssel, and Dutch Deep Red breeds, reflecting their common histories. A considerable amount of gene flow from Red Holstein to German Angler and to German Red and White Dual-Purpose was revealed, which is consistent with frequent crossbreeding to improve productivity of these local breeds. In Red Holstein, marked genomic signatures of selection were reported on chromosome 18, suggesting directed selection for important breeding goal traits. Furthermore, tests for signatures of selection between Red Holstein, Red and White Dual-Purpose, and Meuse-Rhine-Yssel uncovered signals for all investigated pairs of populations. The corresponding genomic regions, which were putatively under different selection pressures, harboured various genes which are associated with traits such as milk and beef production, mastitis and female fertility.

Conclusions

This study provides comprehensive knowledge on the genetic constitution and genomic connectedness of divergent red cattle populations in Northern Europe. The results will help to design and optimize breeding strategies. A joint genomic evaluation including some of the breeds studied here seems feasible.

Why genetic selection to reduce the prevalence of infectious diseases is way more promising than currently believed

Genetic selection for improved disease resistance is an important part of strategies to combat infectious diseases in agriculture. Quantitative genetic analyses of binary disease status, however, indicate low heritability for most diseases, which restricts the rate of genetic reduction in disease prevalence. Moreover, the common liability threshold model suggests that eradication of an infectious disease via genetic selection is impossible because the observed-scale heritability goes to zero when the prevalence approaches zero. From infectious disease epidemiology, however, we know that eradication of infectious diseases is possible, both in theory and practice, because of positive feedback mechanisms leading to the phenomenon known as herd immunity. The common quantitative genetic models, however, ignore these feedback mechanisms. Here, we integrate quantitative genetic analysis of binary disease status with epidemiological models of transmission, aiming to identify the potential response to selection for reducing the prevalence of endemic infectious diseases. The results show that typical heritability values of binary disease status correspond to a very substantial genetic variation in disease susceptibility among individuals. Moreover, our results show that eradication of infectious diseases by genetic selection is possible in principle. These findings strongly disagree with predictions based on common quantitative genetic models, which ignore the positive feedback effects that occur when reducing the transmission of infectious diseases. Those feedback effects are a specific kind of Indirect Genetic Effects; they contribute substantially to the response to selection and the development of herd immunity (i.e., an effective reproduction ratio less than one).

RNA-Seq Whole Transcriptome Analysis of Bovine Mammary Epithelial Cells in Response to Intracellular Staphylococcus aureus

Staphylococcus aureus (S. aureus), a common mastitis pathogen widespread in the natural environment of dairy farms, is capable of invading mammary epithelial cells making treatment difficult. However, the mechanism of the response of bovine mammary epithelial cell to S. aureus invasion remains elusive. In this study, transcriptomic analysis and bioinformatics tools were applied to explore the differentially expressed RNAs in bovine mammary epithelial cells (bMECs) between the control and S. aureus-treated group. A total of 259 differentially expressed mRNAs (DEmRNAs), 27 differentially expressed microRNAs (DEmiRNAs), and 21 differentially expressed long non-coding RNAs (DElncRNAs) were found. These RNAs mainly enrich the inflammatory response, immune response, endocytosis, and cytokine-cytokine receptor interaction. qRT-PCR was used to analyze the quality of the RNA-seq results. In particular, to the defense mechanism of bovine mammary epithelial cells against intracellular S. aureus, the PPAR signaling pathway and the genes (ACOX2, CROT, and NUDT12) were found to be up-regulated to promote the production of peroxisomes and ROS, DRAM1 expression was also up-regulated to facilitate the activation of autophagy, indicating that the above mechanisms were involved in the elimination of intracellular S. aureus in bovine mammary epithelial cells.

Role of the JAK-STAT Pathway in Bovine Mastitis and Milk Production

Simple Summary

The cytokine-activated Janus kinase (JAK)—signal transducer and activator of transcription (STAT) pathway has an important role in the regulation of immunity and inflammation. In addition, the signaling of this pathway has been reported to be associated with mammary gland development and milk production. Because of such important functions, the JAK-STAT pathway has been widely targeted in both human and animal diseases as a therapeutic agent. Recently, the JAK2, STATs, and inhibitors of the JAK-STAT pathway, especially cytokine signaling suppressors (SOCSs), have been reported to be associated with milk production and mastitis-resistance phenotypic traits in dairy cattle. Thus, in the current review, we attempt to overview the development of the JAK-STAT pathway role in bovine mastitis and milk production.

Abstract

The cytokine-activated Janus kinase (JAK)—signal transducer and activator of transcription (STAT) pathway is a sequence of communications between proteins in a cell, and it is associated with various processes such as cell division, apoptosis, mammary gland development, lactation, anti-inflammation, and immunity. The pathway is involved in transferring information from receptors on the cell surface to the cell nucleus, resulting in the regulation of genes through transcription. The Janus kinase 2 (JAK2), signal transducer and activator of transcription A and B (STAT5 A & B), STAT1, and cytokine signaling suppressor 3 (SOCS3) are the key members of the JAK-STAT pathway. Interestingly, prolactin (Prl) also uses the JAK-STAT pathway to regulate milk production traits in dairy cattle. The activation of JAK2 and STATs genes has a critical role in milk production and mastitis resistance. The upregulation of SOCS3 in bovine mammary epithelial cells inhibits the activation of JAK2 and STATs genes, which promotes mastitis development and reduces the lactational performance of dairy cattle. In the current review, we highlight the recent development in the knowledge of JAK-STAT, which will enhance our ability to devise therapeutic strategies for bovine mastitis control. Furthermore, the review also explores the role of the JAK-STAT pathway in the regulation of milk production in dairy cattle.

Lipopolysaccharide triggers different transcriptional signatures in taurine and indicine cattle macrophages: Reactive oxygen species and potential outcomes to the development of immune response to infections

Macrophages are classified upon activation as classical activated M1 and M2 anti-inflammatory regulatory populations. This macrophage polarization is well characterized in humans and mice, but M1/M2 profile in cattle has been far less explored. Bos primigenius taurus (taurine) and Bos primigenius indicus (indicine) cattle display contrasting levels of resistance to infection and parasitic diseases such as C57BL/6J and Balb/c murine experimental models of parasite infection outcomes based on genetic background. Thus, we investigated the differential gene expression profile of unstimulated and LPS stimulated monocyte-derived macrophages (MDMs) from Holstein (taurine) and Gir (indicine) breeds using RNA sequencing methodology. For unstimulated MDMs, the contrast between Holstein and Gir breeds identified 163 Differentially Expressed Genes (DEGs) highlighting the higher expression of C-C chemokine receptor type five (CCR5) and BOLA-DQ genes in Gir animals. LPS-stimulated MDMs from Gir and Holstein animals displayed 1,257 DEGs enriched for cell adhesion and inflammatory responses. Gir MDMs cells displayed a higher expression of M1 related genes like Nitric Oxide Synthase 2 (NOS2), Toll like receptor 4 (TLR4), Nuclear factor NF-kappa-B 2 (NFKB2) in addition to higher levels of transcripts for proinflammatory cytokines, chemokines, complement factors and the acute phase protein Serum Amyloid A (SAA). We also showed that gene expression of inflammatory M1 population markers, complement and SAA genes was higher in Gir in buffy coat peripheral cells in addition to nitric oxide concentration in MDMs supernatant and animal serum. Co-expression analyses revealed that Holstein and Gir animals showed different transcriptional signatures in the MDMs response to LPS that impact on cell cycle regulation, leukocyte migration and extracellular matrix organization biological processes. Overall, the results suggest that Gir animals show a natural propensity to generate a more pronounced M1 inflammatory response than Holstein, which might account for a faster immune response favouring resistance to many infection diseases.

Integrating Single-Step GWAS and Bipartite Networks Reconstruction Provides Novel Insights into Yearling Weight and Carcass Traits in Hanwoo Beef Cattle

Simple Summary

Hanwoo is an indigenous cattle breed in Korea and popular for meat production owing to its rapid growth and high-quality meat. Its yearling weight and carcass traits (backfat thickness, carcass weight, eye muscle area, and marbling score) are economically important for the selection of young and proven bulls. In recent decades, the advent of high throughput genotyping technologies has made it possible to perform genome-wide association studies (GWAS) for the detection of genomic regions associated with traits of economic interest in different species. In this study, we conducted a weighted single-step genome-wide association study which combines all genotypes, phenotypes and pedigree data in one step (ssGBLUP). It allows for the use of all SNPs simultaneously along with all phenotypes from genotyped and ungenotyped animals. Our results revealed 33 relevant genomic regions related to the traits of interest. Gene set enrichment analysis indicated that the identified candidate genes were related to biological processes and functional terms that were involved in growth and lipid metabolism. In conclusion, these results suggest that the incorporation of GWAS results and network analysis can help us to better understand the genetic bases underlying growth and carcass traits.

Abstract

In recent years, studies on the biological mechanisms underlying complex traits have been facilitated by innovations in high-throughput genotyping technology. We conducted a weighted single-step genome-wide association study (WssGWAS) to evaluate backfat thickness, carcass weight, eye muscle area, marbling score, and yearling weight in a cohort of 1540 Hanwoo beef cattle using BovineSNP50 BeadChip. The WssGWAS uncovered thirty-three genomic regions that explained more than 1% of the additive genetic variance, mostly located on chromosomes 6 and 14. Among the identified window regions, seven quantitative trait loci (QTL) had pleiotropic effects and twenty-six QTL were trait-specific. Significant pathways implicated in the measured traits through Gene Ontology (GO) term enrichment analysis included the following: lipid biosynthetic process, regulation of lipid metabolic process, transport or localization of lipid, regulation of growth, developmental growth, and multicellular organism growth. Integration of GWAS results of the studied traits with pathway and network analyses facilitated the exploration of the respective candidate genes involved in several biological functions, particularly lipid and growth metabolism. This study provides novel insight into the genetic bases underlying complex traits and could be useful in developing breeding schemes aimed at improving growth and carcass traits in Hanwoo beef cattle.

Single-step genome-wide association studies (GWAS) and post-GWAS analyses to identify genomic regions and candidate genes for milk yield in Brazilian Girolando cattle

Milk production is economically important to the Brazilian agribusiness, and the majority of the country's milk production derives from Girolando (Gir × Holstein) cows. This study aimed to identify quantitative trait loci (QTL) and candidate genes associated with 305-d milk yield (305MY) in Girolando cattle. In addition, we investigated the SNP-specific variances for Holstein and Gir breeds of origin within the sequence of candidate genes. A single-step genomic BLUP procedure was used to identify QTL associated with 305MY, and the most likely candidate genes were identified through follow-up analyses. Genomic breeding values specific for Holstein and Gir were estimated in the Girolando animals using a model that uses breed-specific partial relationship matrices, which were converted to breed of origin SNP effects. Differences between breed of origin were evaluated by comparing estimated SNP variances between breeds. From 10 genome regions explaining most additive genetic variance for 305MY in Girolando cattle, 7 candidate genes were identified on chromosomes 1, 4, 6, and 26. Within the sequence of these 7 candidate genes, Gir breed of origin SNP alleles showed the highest genetic variance. These results indicated QTL regions that could be further explored in genomic selection panels and which may also help in understanding the gene mechanisms involved in milk production in the Girolando breed.

Genetics of Omega-3 Long-Chain Polyunsaturated Fatty Acid Metabolism and Meat Eating Quality in Tattykeel Australian White Lambs

Meat eating quality with a healthy composition hinges on intramuscular fat (IMF), fat melting point (FMP), tenderness, juiciness, flavour and omega-3 long-chain polyunsaturated fatty acids (n-3 LC-PUFA) content. These health-beneficial n-3 LC-PUFA play significant roles in optimal cardiovascular, retinal, maternal and childhood brain functions, and include alpha linolenic (ALA), eicosapentaenoic (EPA), docosahexaenoic (DHA) and docosapentaenoic (DPA) acids. The primary objective of this review was to access, retrieve, synthesise and critically appraise the published literature on the synthesis, metabolism and genetics of n-3 LC-PUFA and meat eating quality. Studies on IMF content, FMP and fatty acid composition were reviewed to identify knowledge gaps that can inform future research with Tattykeel Australian White (TAW) lambs. The TAW is a new sheep breed exclusive to MARGRA brand of lamb with an outstanding low fat melting point (28-39°C), high n-3 LC-PUFA EPA+DHA content (33-69mg/100g), marbling (3.4-8.2%), tenderness (20.0-38.5N) and overall consumer liking (7.9-8.5). However, correlations between n-3 LC-PUFA profile, stearoyl-CoA desaturase (SCD), fatty acid binding protein 4 (FABP4), fatty acid synthase (FASN), other lipogenic genes and meat quality traits present major knowledge gaps. The review also identified research opportunities in nutrition-genetics interactions aimed at a greater understanding of the genetics of n-3 LC-PUFA, feedlot finishing performance, carcass traits and eating quality in the TAW sheep. It was concluded that studies on IMF, FMP and n-3 LC-PUFA profiles in parental and progeny generations of TAW sheep will be foundational for the genetic selection of healthy lamb eating qualities and provide useful insights into their correlations with SCD, FASN and FABP4 genes.

Short communication: Time-dependent genetic parameters and single-step genome-wide association analyses for predicted milk fatty acid composition in Ayrshire and Jersey dairy cattle

Milk fat content and fatty acid (FA) composition have great economic value to the dairy industry as they are directly associated with taste and chemical-physical characteristics of milk and dairy products. In addition, consumers' choices are not only based on the nutritional aspects of food, but also on products known to promote better health. Milk FA composition is also related to the metabolic status and physiological stages of cows and thus can also be used as indicator for other novel traits of interest (e.g., metabolic diseases and methane yield). Genetic selection is a promising alternative to manipulate milk FA composition. In this study, we aimed to (1) estimate time-dependent genetic parameters for 5 milk FA groups (i.e., short-chain, medium-chain, long-chain, saturated, and unsaturated) predicted based on milk mid-infrared spectroscopy, for Canadian Ayrshire and Jersey breeds, and (2) conduct a time-dependent, single-step genome-wide association study to identify genomic regions, candidate genes, and metabolic pathways associated with milk FA. We analyzed 31,709 test-day records of 9,648 Ayrshire cows from 268 herds, and 34,341 records of 11,479 Jersey cows from 883 herds. The genomic database contained a total of 2,330 Ayrshire and 1,019 Jersey animals. The average daily heritability ranged from 0.18 (long-chain FA) to 0.34 (medium-chain FA) in Ayrshire, and from 0.25 (long-chain and unsaturated FA) to 0.52 (medium-chain and saturated FA) in Jersey. Important genomic regions were identified in Bos taurus autosomes BTA3, BTA5, BTA12, BTA13, BTA14, BTA16, BTA18, BTA20, and BTA21. The proportion of the variance explained by 20 adjacent SNP ranged from 0.71% (saturated FA) to 1.11% (long-chain FA) in Ayrshire, and from 0.70% (unsaturated FA) to 3.09% (medium-chain FA) in Jersey cattle. Important candidate genes and pathways were also identified, such as the PTK2 and TRAPPC9 genes, associated with milk fat percentage, and HMGCS, FGF10, and C6 genes, associated with fertility traits and immune response. Our findings on the genetic parameters and candidate genes contribute to a better understanding of the genetic architecture of milk FA composition in Ayrshire and Jersey dairy cattle.

Genome-wide association study for carcass quality traits and growth in purebred and crossbred pigs1

Carcass quality traits such as back fat (BF), loin depth (LD), and ADG are of extreme economic importance for the swine industry. This study aimed to (i) estimate the genetic parameters for such traits and (ii) conduct a single-step genome-wide association study (ssGWAS) to identify genomic regions that affect carcass quality and growth traits in purebred (PB) and three-way crossbred (CB) pigs. A total of 28,497 PBs and 135,768 CBs pigs were phenotyped for BF, LD, and ADG. Of these, 4,857 and 3,532 were genotyped using the Illumina PorcineSNP60K Beadchip. After quality control, 36,328 SNPs were available and were used to perform an ssGWAS. A bootstrap analysis (n = 1,000) and a signal enrichment analysis were performed to declare SNP significance. Genome regions were based on the variance explained by significant 10-SNP sliding windows. Estimates of PB heritability (SE) were 0.42 (0.019) for BF, 0.39 (0.020) for LD, and 0.35 (0.021) for ADG. Estimates of CB heritability were 0.49 (0.042) for BF, 0.27 (0.029) for LD, and 0.12 (0.021) for ADG. Genetic correlations (SE) across the two populations were 0.81 (0.02), 0.79 (0.04), and 0.56 (0.05), for BF, LD, and ADG, respectively. The variance explained by significant regions for each trait in PBs ranged from 1.51% to 1.35% for BF, from 4.02% to 3.18% for LD, and from 2.26% to 1.45% for ADG. In CBs, the variance explained by significant regions ranged from 1.88% to 1.37% for BF, from 1.29% to 1.23% for LD, and from 1.54% to 1.32% for ADG. In this study, we have described regions of the genome that determine carcass quality and growth traits of PB and CB pigs. These results provide evidence that there are overlapping and nonoverlapping regions in the genome influencing carcass quality and growth traits in PBs and three-way CB pigs.

Genetic and genomic analyses underpin the feasibility of concomitant genetic improvement of milk yield and mastitis resistance in dairy sheep

Milk yield is the most important dairy sheep trait and constitutes the key genetic improvement goal via selective breeding. Mastitis is one of the most prevalent diseases, significantly impacting on animal welfare, milk yield and quality, while incurring substantial costs. Our objectives were to determine the feasibility of a concomitant genetic improvement programme for enhanced milk production and resistance to mastitis. Individual records for milk yield, and four mastitis-related traits (milk somatic cell count, California Mastitis Test score, total viable bacterial count in milk and clinical mastitis presence) were collected monthly throughout lactation for 609 ewes of the Chios breed. All ewes were genotyped with a mastitis specific custom-made 960 single nucleotide polymorphism (SNP) array. We performed targeted genomic association studies, (co)variance component estimation and pathway enrichment analysis, and characterised gene expression levels and the extent of allelic expression imbalance. Presence of heritable variation for milk yield was confirmed. There was no significant genetic correlation between milk yield and mastitis traits. Environmental factors appeared to favour both milk production and udder health. There were no overlapping of SNPs associated with mastitis resistance and milk yield in Chios sheep. Furthermore, four distinct Quantitative Trait Loci (QTLs) affecting milk yield were detected on chromosomes 2, 12, 16 and 19, in locations other than those previously identified to affect mastitis resistance. Five genes (DNAJA1, GHR, LYPLA1, NUP35 and OXCT1) located within the QTL regions were highly expressed in both the mammary gland and milk transcriptome, suggesting involvement in milk synthesis and production. Furthermore, the expression of two of these genes (NUP35 and OXCT1) was enriched in immune tissues implying a potentially pleiotropic effect or likely role in milk production during udder infection, which needs to be further elucidated in future studies. In conclusion, the absence of genetic antagonism between milk yield and mastitis resistance suggests that simultaneous genetic improvement of both traits be achievable.

Genome-wide association study applied to type traits related to milk yield in water buffaloes (Bubalus bubalis)

This research aimed to estimate genetic parameters for milk yield and type traits [withers height (WH), croup height (CH), body length (BL), croup length (CL), iliac width (ILW), ischial width (ISW), and thoracic circumference] in Murrah buffaloes and to identify genomic regions related to type traits by applying a single-step genome-wide association study. Data used to estimate the genetic parameters consisted of 601 records of milk yield in the first lactation and the aforementioned type traits. For the single-step genome-wide association study, 322 samples genotyped with a 90K Axiom Buffalo Genotyping array (Thermo Fisher Scientific, Santa Clara, CA) were used. Bivariate analysis revealed that heritability for milk yield (kg) at 305 d was 0.31 ± 0.11, whereas it ranged from 0.22 ± 0.07 to 0.34 ± 0.09 for the studied conformation traits. Based on the percentages of genetic variance explained by windows of 10 markers, there were 16 genomic regions explaining more than 0.5% of the variance for WH, CH, BL, CL, ILW, ISW, and thoracic circumference. Between those regions, 4 were associated with more than 1 trait, suggesting pleiotropic roles for some genes of Bos taurus autosome (BTA) 12 on CL and WH, BTA13 on ISW and ILW, BTA23 on CH and BL, and BTA28 on ISW and BL. Most of these regions coincide with known quantitative trait loci for milk traits. Thus, further studies based on sequence data will help to validate the association of this region with type traits and likely identify the causal mutations.

On the genomic regions associated with milk lactose in Fleckvieh cattle

Lactose is a sugar uniquely found in mammals' milk and it is the major milk solid in bovines. Lactose yield (LY, kg/d) is responsible for milk volume, whereas lactose percentage (LP) is thought to be more related to epithelial integrity and thus to udder health. There is a paucity of studies that have investigated lactose at the genomic level in dairy cows. This paper aimed to improve our knowledge on LP and LY, providing new insights into the significant genomic regions affecting these traits. A genome-wide association study for LP and LY was carried out in Fleckvieh cattle by using bulls' deregressed estimated breeding values of first lactation as pseudo-phenotypes. Heritabilities of first-lactation test-day LP and LY estimated using linear animal models were 0.38 and 0.25, respectively. A total of 2,854 bulls genotyped with a 54K SNP chip were available for the genome-wide association study; a linear mixed model approach was adopted for the analysis. The significant SNP of LP were scattered across the whole genome, with signals on chromosomes 1, 2, 3, 7, 12, 16, 18, 19, 20, 28, and 29; the top 4 significant SNP explained 4.90% of the LP genetic variance. The signals were mostly in regions or genes with involvement in molecular intra- or extracellular transport; for example, CDH5, RASGEF1C, ABCA6, and SLC35F3. A significant region within chromosome 20 was previously shown to affect mastitis or somatic cell score in cattle. As regards LY, the significant SNP were concentrated in fewer regions (chromosomes 6 and 14), related to mastitis/somatic cell score, immune response, and transport mechanisms. The 5 most significant SNP for LY explained 8.45% of genetic variance and more than one-quarter of this value has to be attributed to the variant within ADGRB1. Significant peaks in target regions remained even after adjustment for the 2 most significant variants previously detected on BTA6 and BTA14. The present study is a prelude for deeper investigations into the biological role of lactose for milk secretion and volume determination, stressing the connection with genes regulating intra- or extracellular trafficking and immune and inflammatory responses in dairy cows. Also, these results improve the knowledge on the relationship between lactose and udder health; they support the idea that LP and its derived traits are potential candidates as indicators of udder health in breeding programs aimed to enhance cows' resistance to mastitis.

Three Novel Players: PTK2B, SYK, and TNFRSF21 Were Identified to Be Involved in the Regulation of Bovine Mastitis Susceptibility via GWAS and Post-transcriptional Analysis

Bovine mastitis is a common inflammatory disease caused by multiple factors in early lactation or dry period. Genome wide association studies (GWAS) can provide a convenient and effective strategy for understanding the biological basis of mastitis and better prevention. 2b-RADseq is a high-throughput sequencing technique that offers a powerful method for genome-wide genetic marker development and genotyping. In this study, single nucleotide polymorphisms (SNPs) of the immune-regulated gene correlative with mastitis were screened and identified by two stage association analysis via GWAS-2b-RADseq in Chinese Holstein cows. We have screened 10,058 high quality SNPs from 7,957,920 tags and calculated their allele frequencies. Twenty-seven significant SNPs were co-labeled in two GWAS analysis models [Bayesian (P < 0.001) and Logistic regression (P < 0.01)], and only three SNPs (rs75762330, C > T, PIC = 0.2999; rs88640083, A > G, PIC = 0.1676; rs20438858, G > A, PIC = 0.3366) were annotated to immune-regulated genes (PTK2B, SYK, and TNFRSF21). Identified three SNPs are located in non-coding regions with low or moderate genetic polymorphisms. However, independent sample population validation (Case-control study) data showed that three important SNPs (rs75762330, P < 0.025, OR > 1; rs88640083, P < 0.005, OR > 1; rs20438858, P < 0.001, OR < 1) were significantly associated with clinical mastitis trait. Importantly, PTK2B and SYK expression was down-regulated in both peripheral blood leukocytes (PBLs) of clinical mastitis cows and in vitro LPS (E. coli)-stimulated bovine mammary epithelial cells, while TNFRSF21 was up-regulated. Under the same conditions, expression of Toll-like receptor 4 (TLR4), AKT1, and pro-inflammatory factors (IL-1β and IL-8) were also up-regulated. Interestingly, network analysis indicated that PTK2B and SYK are co-expressed in innate immune signaling pathway of Chinese Holstein. Taken together, these results provided strong evidence for the study of SNPs in bovine mastitis, and revealed the role of SYK, PTK2B, and TNFRSF21 in bovine mastitis susceptibility/tolerance.

Paraoxonase 1 in bovine milk and blood as marker of subclinical mastitis caused by Staphylococcus aureus

The aim of this study was to determine serum lipid values and parameters of oxidative stress in blood and milk of cows with subclinical mastitis (SCM) caused by Staphylococcus aureus and to establish association between these parameters. The study was performed on total of 104 cows assigned into control group of healthy animals (n = 12) and two groups of cows with SCM, either SCM₁ group (n = 37) with ˂ 1000 CFU/mL of S. aureus or SCM₂ group (n = 55) with ≥1000 CFU/mL of causative agent in milk. Significantly lower serum concentrations of high density lipoprotein-cholesterol (HDLC) in SCM₂ group and higher low density lipoprotein-cholesterol (LDL-C) in both SCM groups were recorded. Significantly lower paraoxonase-1 (PON1) activity and higher lipid hydroperoxides (LOOH) concentration in blood and milk were recorded in both groups of SCM cows vs control. In blood serum of cows from SCM₂ group significantly higher concentrations of hydrogen peroxide (H₂O₂) and malondialdehyde (MDA), but significantly lower total antioxidative capacity (TAC) were obtained. In milk serum of SCM cows were detected significantly higher concentrations of H₂O₂ and MDA, but significantly lower TAC vs control. There was a significant positive correlation between PON1 in serum with HDLC, but negative correlation with LDL-C and LOOH. In milk, PON1 negatively correlated with LOOH and somatic cell counts (SCC). Strong positive correlation was obtained between PON1 in blood and milk. Oxidative stress and inflammatory reaction induced by SCM significantly lowered PON1 activity in blood and milk of affected cows.

Genome-wide association analysis for β-hydroxybutyrate concentration in Milk in Holstein dairy cattle

BACKGROUND

Ketosis in dairy cattle has been shown to cause a high morbidity in the farm and substantial financial losses to dairy farmers. Ketosis symptoms, however, are difficult to identify, therefore, the amount of ketone bodies (mainly β-hydroxybutyric acid, BHB) is used as an indicator of subclinical ketosis in cows. It has also been shown that milk BHB concentrations have a strong correlation with ketosis in dairy cattle. Mid-infrared spectroscopy (MIR) has recently became a fast, cheap and high-throughput method for analyzing milk components. The aim of this study was to perform a genome-wide association study (GWAS) on the MIR-predicted milk BHB to identify genomic regions, genes and pathways potentially affecting subclinical ketosis in North American Holstein dairy cattle.

RESULTS

Several significant regions were identified associated with MIR-predicted milk BHB concentrations (indicator of subclinical ketosis) in the first lactation (SCK1) and second and later lactations (SCK2) in Holstein dairy cows. The strongest association was located on BTA6 for SCK1 and BTA14 on SCK2. Several SNPs on BTA6 were identified in regions and variants reported previously to be associated with susceptibility to ketosis and clinical mastitis in Jersey and Holstein dairy cattle, respectively. One highly significant SNP on BTA14 was found within the DGAT1 gene with known functions on fat metabolism and inflammatory response in dairy cattle. A region on BTA6 and three SNPs on BTA20 were found to overlap between SCK1 and SCK2. However, a novel region on BTA20 (55-63 Mb) for SCK2 was also identified, which was not reported in previous association studies. Enrichment analysis of the list of candidate genes within the identified regions for MIR-predicted milk BHB concentrations yielded molecular functions and biological processes that may be involved in the inflammatory response and lipid metabolism in dairy cattle.

CONCLUSIONS

The results of this study confirmed several SNPs and genes identified in previous studies as associated with ketosis susceptibility and immune response, and also found a novel region that can be used for further analysis to identify causal variations and key regulatory genes that affect clinical/ subclinical ketosis.

Frequentist p-values for large-scale-single step genome-wide association, with an application to birth weight in American Angus cattle

Background

Single-step genomic best linear unbiased prediction (SSGBLUP) is a comprehensive method for genomic prediction. Point estimates of marker effects from SSGBLUP are often used for genome-wide association studies (GWAS) without a formal framework of hypothesis testing. Our objective was to implement p-values for single-marker GWAS studies within the single-step GWAS (SSGWAS) framework by deriving computational algorithms and procedures, and by applying these to a large beef cattle population.

Methods

P-values were obtained based on the prediction error (co)variances for single nucleotide polymorphisms (SNPs), which were obtained from the prediction error (co)variances of genomic predictions based on the inverse of the coefficient matrix and formulas to estimate SNP effects.

Results

Computation of p-values took a negligible time for a dataset with almost 2 million animals in the pedigree and 1424 genotyped sires, and no inflation of statistics was observed. The SNPs that passed the Bonferroni threshold of 10^−5.9 were the same as those that explained the highest proportion of additive genetic variance, but even at the same significance levels and effects, some of them explained less genetic variance due to lower allele frequency.

Conclusions

The use of a p-value for SSGWAS is a very general and efficient strategy to identify quantitative trait loci (QTL). It can be used for complex datasets such as those used in animal breeding, where only a proportion of the pedigreed animals are genotyped.

Genome-wide association study identifies QTLs for displacement of abomasum in Chinese Holstein cattle1

Displacement of abomasum (DA) is one of the most common and important disorders in dairy cattle. The objective of the present study was to detect the quantitative trait loci (QTL) for DA in Chinese Holstein using single-step genomic BLUP methodology. A total of 60,556 producer-recorded DA event records from 32,190 cows, together with 2,336 genotyped animals with 40,054 SNP markers, were used for the analysis. Genomic data were incorporated into a threshold model for variance component estimation, and the estimated heritability of DA was 0.108 (SE = 0.086). Results of genome-wide association studies were reported as the proportion of genetic variance explained 20-SNP windows. Eight QTLs covering 129 genes on Bos taurus autosomes 2, 4, 7, 10, 14, 17, 20 showed associations with DA. Ten genes, namely BMP4, SOCS4, GCH1, DDHD1, ATG14, ACBP/DBI, SMO, AHCYL2, CYP7A1, and CACNA1A, involved in insulin metabolism and lipid metabolism pathways may be considered as candidate genes of DA in dairy. The identified QTLs, biological pathways, and associated genes underlying DA identified from the present study will contribute to the understanding of the genetic architecture of this complex disease.