The genetic architecture of milk ELISA scores as an indicator of Johne's disease (paratuberculosis) in dairy cattle

Johne's disease (or paratuberculosis), caused by Mycobacterium avium ssp. paratuberculosis (MAP) infection, is a globally prevalent disease with severe economic and welfare implications. With no effective treatment available, understanding the role of genetics influencing host infection status is essential to develop selection strategies to breed for increased resistance to MAP infection. The main objectives of this study were to estimate genetic parameters for the MAP-specific antibody response using milk ELISA scores in Canadian Holstein cattle as an indicator of resistance to Johne's disease, and to unravel genomic regions and candidate genes significantly associated with MAP infection. After data editing, 168,987 milk ELISA records from 2,306 herds, obtained from CanWest Dairy Herd Improvement, were used for further analyses. Variance and heritability estimates for MAP infection status were determined using univariate linear animal models under 3 scenarios: (a) SCEN1: the complete data set (all herds); (b) SCEN2: herds with at least one suspect or test-positive animal (ELISA optical density ≥0.07); and (c) SCEN3: herds with at least one test-positive animal (ELISA optical density ≥0.11). Heritability estimates were calculated as 0.066, 0.064, and 0.063 for SCEN1, SCEN2, and SCEN3, respectively. The correlations between estimated breeding values for resistance to MAP infection and other economically important traits, when significant, were favorable and of low magnitude. Genome-wide association analyses identified important genomic regions on Bos taurus autosome (BTA)1, BTA7, BTA9, BTA14, BTA15, BTA17, BTA19, and BTA25 showing significant association with MAP infection status. These regions included 2 single nucleotide polymorphisms located 2 kb upstream of positional candidate genes CD86 and WNT9B, which play key roles in host immune response and tissue homeostasis. This study revealed the genetic architecture of MAP infection in Canadian Holstein cattle as measured by milk ELISA scores by estimating genetic parameters along with the identification of genomic regions potentially influencing MAP infection status. These findings will be of significant value toward implementing genetic and genomic evaluations for resistance to MAP infection in Holstein cattle.     

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.     

Genome-wide association study of normal and atypical progesterone profiles in Holstein-Friesian dairy cows

Female fertility has a major role in dairy production and affects the profitability of dairy cattle. The genetic progress obtained by traditional selection can be slow because of the low heritability of classical fertility traits. Endocrine fertility traits based on progesterone concentration in milk have higher heritability and more directly reflect the cow's own reproductive physiology. The aim of our study was to identify genomic regions for 7 endocrine fertility traits in dairy cows by performing a genome-wide association study with 54,000 SNP. The next step was to fine-map targeted genomic regions with significant SNP using imputed sequences to identify potential candidate genes associated with the normal and atypical progesterone profiles. The association between a SNP and a phenotype was assessed by a single SNP analysis, using a linear mixed model that included a random polygenic effect. Phenotypes and genotypes were available for 1,126 primiparous and multiparous Holstein-Friesian cows from research herds in Ireland, the Netherlands, Sweden, and the United Kingdom. In total, 44 significant SNP associated with 7 endocrine fertility traits were identified on Bos taurus autosome (BTA) 1–4, 6, 8–9, 11–12, 14–17, 19, 21–24, and 29. Three chromosomes, BTA8, BTA17, and BTA23, were imputed from 54,000 SNP genotypes to the whole-genome sequence level with Beagle version 4.1. The fine-mapping identified several significant associations with delayed cyclicity, cessation of cyclicity, commencement of luteal activity, and inter-ovulatory interval. These associations may contribute to an index of markers for genetic improvement of fertility. Several potential candidate genes reported to affect reproduction were also identified in the targeted genomic regions. However, due to high linkage disequilibrium, it was not possible to identify putative causal genes or polymorphisms for any of the regions.     

Fine-mapping of BTA17 using imputed sequences for associations with de novo synthesized fatty acids in bovine milk

The aim of this study was to fine-map a genomic region associated with milk fatty acids (FA) on Bos taurus autosome (BTA) 17. This genomic region has been discovered with 50,000 (50k) single nucleotide polymorphisms (SNP) imputed to 777,000 (777k) SNP. In this study, high-density genotypes were imputed to whole-genome sequences level to identify candidate gene(s) associated with milk FA composition on BTA17. Phenotypes and genotypes were available for 1,640 cows sampled in winter, and for 1,581 cows sampled in summer. Phenotypes consisted of gas chromatography measurements in winter and in summer milk samples of 6 individual FA and the indicator of de novo synthesis, C6:0–C14:0. Genotypes consisted of imputed 777k SNP, and 89 sequenced ancestors of the population of genotyped cows. In addition, 450 whole-genome sequences from the 1,000 Bull Genome Consortium were available. Using 495 Holstein-Friesian sequences as a reference population, the 777k SNP genotypes of the cows were imputed to sequence level. We then applied single-variant analyses with an animal model, and identified thousands of significant associations with C6:0, C8:0, C10:0, C12:0, C14:0, and C6:0–C14:0. For C8:0 in summer milk samples, the genomic region located between 29 and 34 Mbp on BTA17 revealed a total of 646 significant associations. The most significant associations [–log₁₀(P-value) = 7.82] were 8 SNP in perfect linkage disequilibrium. After fitting one of these 8 SNP as a fixed effect in the model, and re-running the single-variant analyses, no further significant associations were found for any of the 6 FA or C6:0–C14:0. These findings suggest that one polymorphism underlying this QTL on BTA17 influences multiple de novo synthesized milk FA. Thirteen genes in the QTL region were identified and analyzed carefully. Six out of the 8 SNP that showed the strongest associations were located in the La ribonucleoprotein domain family, member 1B (LARP1B) gene, and we suggest LARP1B as a primary candidate gene. Another gene of interest for this QTL region might be PKL4. None of these suggested candidate genes have previously been associated with milk fat synthesis or milk FA composition.     

Genetic and nongenetic profiling of milk pregnancy-associated glycoproteins in Holstein cattle

Pregnancy-associated glycoproteins (PAG) are secreted by the trophoblast and are detectable in maternal circulation around the time of attachment of the fetal placenta, as well as in blood and milk throughout gestation. The understanding of the genetic mechanisms controlling PAG levels can confer advantages for livestock breeding programs given the precocity and the ease of obtaining this phenotype from routine pregnancy diagnosis. The aim of this study was to characterize PAG levels by estimating genetic parameters and correlations with other dairy traits, and to identify genomic regions and candidate genes associated with PAG levels in milk. The PAG data consisted of pregnancy diagnoses using commercial assays from 2012 to 2017, and genotype data consisted of 54,123 SNP markers for 2,352 individuals (embryos and dams). The model included contemporary group (herd, year, and season) and embryo age as fixed effects, and random embryonic (direct) and maternal (indirect) genetic effects. Using genomic data, the estimated heritability for direct and maternal genetic effects (± standard deviations) were 0.23 ± 0.05 and 0.11 ± 0.05, respectively. The genetic correlation between these effects was almost zero (0.001 ± 0.06). A preliminary analysis revealed low correlations between milk PAG levels and other dairy traits. The genome-wide association analysis was performed using 2 approaches: single-marker and single-step using all markers. Four genomic regions with direct genetic effects were detected on Bos taurus autosome (BTA) 6, BTA7, BTA19, and BTA29 of the embryonic genome. The BTA29 locus was within the bovine PAG gene cluster, suggesting a cis-regulatory quantitative trait locus on the PAG expression. However, other associations, without an obvious link to PAG expression, could be related to the transportation of PAG and their concentration in milk. Only 1 region from the maternal genome, on BTA4, had a significant indirect effect, where WNT2 is a candidate gene related to placenta vascularization and gestation health. Collectively, our results suggest a moderate genetic control of milk PAG levels from both maternal and fetal genomes, but larger studies are needed to fully evaluate the usefulness of milk PAG in the genetic evaluation of fetal growth and cow fertility.     

Whole-genome association study for milk protein composition in dairy cattle

Our objective was to perform a genome-wide association study for content in bovine milk of α_S1-casein (α_S1-CN), α_S2-casein (α_S2-CN), β-casein (β-CN), κ-casein (κ-CN), α-lactalbumin (α-LA), β-lactoglobulin (β-LG), casein index, protein percentage, and protein yield using a 50K single nucleotide polymorphism (SNP) chip. In total, 1,713 Dutch Holstein-Friesian cows were genotyped for 50,228 SNP and a 2-step association study was performed. The first step involved a general linear model and the second step used a mixed model accounting for all family relationships. Associations with milk protein content and composition were detected on 20 bovine autosomes. The main genomic regions associated with milk protein composition or protein percentage were found on chromosomes 5, 6, 11, and 14. The number of chromosomal regions showing significant (false discovery rate <0.01) effects ranged from 3 for β-CN and 3 for β-LG to 12 for α_S2-CN. A genomic region on Bos taurus autosome (BTA) 6 was significantly associated with all 6 major milk proteins, and a genomic region on BTA 11 was significantly associated with the 4 caseins and β-LG. In addition, regions were detected that only showed a significant effect on one of the milk protein fractions: regions on BTA 13 and 22 with effects on α_S1-CN; regions on BTA 1, 9, 10, 17, 19, and 28 with effects on α_S2-CN; a region on BTA 6 with an effect on β-CN; regions on BTA 13 and 21 with effects on κ-CN; regions on BTA 1, 5, 9, 16, 17, and 26 with effects on α-LA; and a region on BTA 24 with an effect on β-LG. The proportion of genetic variance explained by the SNP showing the strongest association in each of these genomic regions ranged from <1% for α_S1-CN on BTA 22 to almost 100% for casein index on BTA 11. Variation associated with regions on BTA 6, 11, and 14 could in large part but not completely be explained by known protein variants of β-CN (BTA 6), κ-CN (BTA 6), and β-LG (BTA 11) or DGAT1 variants (BTA 14). Our results indicate 3 regions with major effects on milk protein composition, in addition to several regions with smaller effects involved in the regulation of milk protein composition.     

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 association study for milk production traits in a Brazilian Holstein population

Advances in the molecular area of selection have expanded knowledge of the genetic architecture of complex traits through genome-wide association studies (GWAS). Several GWAS have been performed so far, but confirming these results is not always possible due to several factors, including environmental conditions. Thus, our objective was to identify genomic regions associated with traditional milk production traits, including milk yield, somatic cell score, fat, protein and lactose percentages, and fatty acid composition in a Holstein cattle population producing under tropical conditions. For this, 75,228 phenotypic records from 5,981 cows and genotypic data of 56,256 SNP from 1,067 cows were used in a weighted single-step GWAS. A total of 46 windows of 10 SNP explaining more than 1% of the genetic variance across 10 Bos taurus autosomes (BTA) harbored well-known and novel genes. The MGST1 (BTA5), ABCG2 (BTA6), DGAT1 (BTA14), and PAEP (BTA11) genes were confirmed within some of the regions identified in our study. Potential novel genes involved in tissue damage and repair of the mammary gland (COL18A1), immune response (LTTC19), glucose homeostasis (SLC37A1), synthesis of unsaturated fatty acids (LTBP1), and sugar transport (SLC37A1 and MFSD4A) were found for milk yield, somatic cell score, fat percentage, and fatty acid composition. Our findings may assist genomic selection by using these regions to design a customized SNP array to improve milk production traits on farms with similar environmental conditions.     

Genome-wide association studies to identify quantitative trait loci affecting milk production traits in water buffalo

Water buffalo is the second largest resource of milk supply around the world, and it is well known for its distinctive milk quality in terms of fat, protein, lactose, vitamin, and mineral contents. Understanding the genetic architecture of milk production traits is important for future improvement by the buffalo breeding industry. The advance of genome-wide association studies (GWAS) provides an opportunity to identify potential genetic variants affecting important economical traits. In the present study, GWAS was performed for 489 buffaloes with 1,424 lactation records using the 90K Affymetrix Buffalo SNP Array (Affymetrix/Thermo Fisher Scientific, Santa Clara, CA). Collectively, 4 candidate single nucleotide polymorphisms (SNP) in 2 genomic regions were found to associate with buffalo milk production traits. One region affecting milk fat and protein percentage was located on the equivalent of Bos taurus autosome (BTA)3, spanning 43.3 to 43.8 Mb, which harbored the most likely candidate genes MFSD14A, SLC35A3, and PALMD. The other region on the equivalent of BTA14 at 66.5 to 67.0 Mb contained candidate genes RGS22 and VPS13B and influenced buffalo total milk yield, fat yield, and protein yield. Interestingly, both of the regions were reported to have quantitative trait loci affecting milk performance in dairy cattle. Furthermore, we suggest that buffaloes with the C allele at AX-85148558 and AX-85073877 loci and the G allele at AX-85106096 locus can be selected to improve milk fat yield in this buffalo-breeding program. Meanwhile, the G allele at AX-85063131 locus can be used as the favorable allele for improving milk protein percentage. Genomic prediction showed that the reliability of genomic estimated breeding values (GEBV) of 6 milk production traits ranged from 0.06 to 0.22, and the correlation between estimated breeding values and GEBV ranged from 0.23 to 0.35. These findings provide useful information to understand the genetic basis of buffalo milk properties and may play a role in accelerating buffalo breeding programs using genomic approaches.     

Estimation of genetic parameters and single-step genome-wide association studies for milk urea nitrogen in Holstein cattle

The main objectives of this study were to estimate genetic parameters for milk urea nitrogen (MUN) in Holstein cattle and to conduct a single-step (ss)GWAS to identify candidate genes associated with MUN. Phenotypic measurements from 24,435 Holstein cows were collected from March 2013 to July 2019 in 9 dairy farms located in the Beijing area, China. A total of 2,029 cows were genotyped using the Illumina 150K Bovine Bead Chip, containing 121,188 SNP. A single-trait repeatability model was used to evaluate the genetic background of MUN. We found that MUN is a trait with low heritability (0.06 ± 0.004) and repeatability (0.12). Considering similar milk production levels, a lower MUN concentration indicates higher nitrogen digestibility. The genetic correlations between MUN and milk yield, net energy concentration, fat percentage, protein percentage, and lactose percentage were positive and ranged from 0.02 to 0.26. The genetic correlation between MUN and somatic cell score (SCS) was negative (?0.18), indicating that animals with higher MUN levels tend to have lower SCS. Both ssGWAS and pathway enrichment analyses were used to explore the genetic mechanisms underlying MUN. A total of 18 SNP (located on BTA11, BTA12, BTA14, BTA17, and BTA18) were found to be significantly associated with MUN. The genes CFAP77, CAMSAP1, CACNA1B, ADGRB1, FARP1, and INTU are considered to be candidate genes for MUN. These candidate genes are associated with important biological processes such as protein and lipid metabolism and binding to specific proteins. This set of candidate genes, metabolic pathways, and their functions provide a better understanding of the genomic architecture and physiological mechanisms underlying MUN in Holstein cattle.     

Genome-wide association analysis and pathways enrichment for lactation persistency in Canadian Holstein cattle

Lactation persistency (LP), defined as the rate of declining milk yield after milk peak, is an economically important trait for dairy cattle. Improving LP is considered a good alternative method for increasing overall milk production because it does not cause the negative energy balance and other health issues that cows experience during peak milk production. However, little is known about the biology of LP. A genome-wide association study (GWAS) and pathway enrichment were used to explore the genetic mechanisms underlying LP. The GWAS was performed using a univariate regression mixed linear model on LP data of 3,796 cows and 44,100 single nucleotide polymorphisms (SNP). Eight and 47 SNP were significantly and suggestively associated with LP, respectively. The 2 most important quantitative trait loci regions for LP were (1) a region from 106 to 108 Mb on Bos taurus autosome (BTA) 5, where the most significant SNP (ARS-BFGL-NGS-2399) was located and also formed a linkage disequilibrium block with 3 other SNP; and (2) a region from 29.3 to 31.3 Mb on BTA 20, which contained 3 significant SNP. Based on physical positions, MAN1C1, MAP3K5, HCN1, TSPAN9, MRPS30, TEX14, and CCL28 are potential candidate genes for LP because the significant SNP were located in their intronic regions. Enrichment analyses of a list of 536 genes in 0.5-Mb flanking regions of significant and suggestive SNP indicates that synthesis of milk components, regulation of cell apoptosis processes and insulin, and prolactin signaling pathways are important for LP. Upstream regulators relevant for LP positional candidate genes were prolactin (PRL), peroxisome proliferator-activated receptor gamma (PPARG), and Erb-B2 receptor tyrosine kinase 2 (ERBB2). Several networks related to cellular development, proliferation and death were significantly enriched for LP positional candidate genes. In conclusion, this study detected several SNP, genes, and interesting regions for fine mapping and validation of candidate genes and SNP for potential use in selection for improved LP. This study also provided further insights on the biology of LP which will help to prioritize selected candidate genes for functional validation and application.     

Genome-wide association study for selected cheese-making properties in Dual-Purpose Belgian Blue cows

This study aimed to estimate genetic parameters and identify genomic region(s) associated with selected cheese-making properties (CMP) in Dual-Purpose Belgian Blue (DPBB) cows. Edited data were 46,301 test-day records of milk yield, fat percentage, protein percentage, casein percentage, milk calcium content (CC), coagulation time (CT), curd firmness after 30 min from rennet addition (a30), and milk titratable acidity (MTA) collected from 2014 to 2020 on 4,077 first-parity (26,027 test-day records), and 3,258 second-parity DPBB cows (20,274 test-day records) distributed in 124 herds in the Walloon Region of Belgium. Data of 28,266 SNP, located on 29 Bos taurus autosomes (BTA) of 1,699 animals 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 25 consecutive SNPs (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. Heritability estimates for the included CMP ranged from 0.19 (CC) to 0.50 (MTA), and 0.24 (CC) to 0.41 (MTA) in the first and second parity, respectively. The genetic correlation estimated between CT and a30 varied from ?0.61 to ?0.41 and from ?0.55 to ?0.38 in the first and second lactations, respectively. Negative genetic correlations were found between CT and milk yield and composition, while those estimated between curd firmness and milk composition were positive. Genome-wide association analyses results identified 4 genomic regions (BTA1, BTA3, BTA7, and BTA11) associated with the considered CMP. The identified genomic regions showed contrasting results between parities and among the different stages of each parity. It suggests that different sets of candidate genes underlie the phenotypic expression of the considered CMP between parities and lactation stages of each parity. The findings of this study can be used for future implementation and use of genomic evaluation to improve the cheese-making traits in DPBB cows.     

Genetic variation of natural antibodies in milk of Dutch Holstein-Friesian cows

Defense mechanisms of dairy cows against diseases partly rest on their naturally present disease resistance capacity. Natural antibodies (NAb) form a soluble part of the innate immune system, being defined as antibodies circulating in animals without prior intentional antigenic stimulation. Genetic selection on NAb titers in milk, therefore, might improve disease resistance. We estimated genetic parameters of NAb titers binding lipopolysaccharide, lipoteichoic acid (LTA), peptidoglycan, and keyhole limpet hemocyanin, and titers of the NAb isotypes IgG1, IgM, and IgA binding LTA in milk of Dutch Holstein-Friesian heifers. Natural antibody titers were measured in 1 milk sample from each of 1,939 Holstein-Friesian heifers and used for estimating genetic parameters of NAb titers. The data show that phenotypic variation exists among heifers in NAb titers binding lipopolysaccharide, LTA, peptidoglycan, and keyhole limpet hemocyanin, and the NAb isotypes IgG1, IgM, and IgA binding LTA in milk. High genetic correlations among NAb (ranging from 0.45 to 0.99) indicated a common genetic basis for the levels of different NAb in bovine milk. Intra-herd heritability estimates for NAb ranged from 0.10 to 0.53. The results indicated that NAb levels have potential for genetic selection.     

Short communication: Replication of genome-wide association studies for milk production traits in Chinese Holstein by an efficient rotated linear mixed model

Milk is regarded as an important nutrient for humans, and Chinese Holstein cows provide high-quality milk for billions of Chinese people. Therefore, detecting quantitative trait nucleotides (QTN) or candidate genes for milk production traits in Chinese Holstein is important. In this study, we performed genome-wide association studies (GWAS) in a Chinese Holstein population of 6,675 cows and 71,633 SNP using deregressed proofs (DRP) as phenotypes to replicate our previous study in a population of 1,815 cows and 39,163 SNP using estimated breeding values (EBV) as phenotypes. The associations between 3 milk production traits—milk yield (MY), fat percentage (FP), and protein percentage (PP)—and the SNP were determined by using an efficient rotated linear mixed model, which benefits from linear transformations of genomic estimated values and Eigen decomposition of the genomic relationship matrix algorithm. In total, we detected 94 SNP that were significantly associated with one or more milk production traits, including 7 SNP for MY, 76 for FP, and 36 for PP; 87% of these SNP were distributed across Bos taurus autosomes 14 and 20. In total, 83 SNP were found to be located within the reported quantitative trait loci (QTL) regions, and one novel segment (between 1.41 and 1.49 Mb) on chromosome 14 was significantly associated with FP, which could be an important candidate QTL region. In addition, the detected intervals were narrowed down from the reported regions harboring causal variants. The top significant SNP for the 3 traits was ARS-BFGL-NGS-4939, which is located within the DGAT1 gene. Five detected genes (CYHR1, FOXH1, OPLAH, PLEC, VPS28) have effects on all 3 traits. Our study provides a suite of QTN, candidate genes, and a novel QTL associated with milk production traits, and thus forms a solid basis for genomic selection and molecular breeding for milk production traits in Chinese Holstein.     

Genome-wide association study for cheese yield and curd nutrient recovery in dairy cows

Cheese production and consumption are increasing in many countries worldwide. As a result, interest has increased in strategies for genetic selection of individuals for technological traits of milk related to cheese yield (CY) in dairy cattle breeding. However, little is known about the genetic background of a cow's ability to produce cheese. Recently, a relatively large panel (1,264 cows) of different measures of individual cow CY and milk nutrient and energy recoveries in the cheese (REC) became available. Genetic analyses showed considerable variation for CY and for aptitude to retain high proportions of fat, protein, and water in the coagulum. For the dairy industry, these characteristics are of major economic importance. Nevertheless, use of this knowledge in dairy breeding is hampered by high costs, intense labor requirement, and lack of appropriate technology. However, in the era of genomics, new possibilities are available for animal breeding and genetic improvement. For example, identification of genomic regions involved in cow CY might provide potential for marker-assisted selection. The objective of this study was to perform genome-wide association studies on different CY and REC measures. Milk and DNA samples from 1,152 Italian Brown Swiss cows were used. Three CY traits expressing the weight (wt) of fresh curd (%CY_CURD), curd solids (%CY_SOLIDS), and curd moisture (%CY_WATER) as a percentage of weight of milk processed, and 4 REC (REC_FAT, REC_PROTEIN, REC_SOLIDS, and REC_ENERGY, calculated as the % ratio between the nutrient in curd and the corresponding nutrient in processed milk) were analyzed. Animals were genotyped with the Illumina BovineSNP50 Bead Chip v.2. Single marker regressions were fitted using the GenABEL R package (genome-wide association using mixed model and regression–genomic control). In total, 103 significant associations (88 single nucleotide polymorphisms) were identified in 10 chromosomes (2, 6, 9, 11, 12, 14, 18, 19, 27, 28). For REC_FAT and REC_PROTEIN, high significance peaks were identified in Bos taurus autosome (BTA) 6 and BTA11, respectively. Marker ARS-BFGL-NGS-104610 (～104.3 Mbp) was highly associated with REC_PROTEIN and Hapmap52348-rs29024684 (～87.4 Mbp), closely located to the casein genes on BTA6, with REC_FAT. Genomic regions identified may enhance marker-assisted selection in bovine cheese breeding beyond the use of protein (casein) and fat contents, whereas new knowledge will help to unravel the genomic background of a cow's ability for cheese production.     

Inbreeding coefficients and runs of homozygosity islands in Brazilian water buffalo

Characterization of autozygosity is relevant to monitor genetic diversity and manage inbreeding levels in breeding programs. Identification of autozygosity hotspots can unravel genomic regions targeted by selection for economically important traits and can help identify candidate genes for selection. In this study, we estimated the inbreeding levels of a Brazilian population of Murrah buffalo undergoing selection for milk production traits, particularly milk yield. We also studied the distribution of runs of homozygosity (ROH) islands and identified putative genes and quantitative trait loci (QTL) under selection. We genotyped 422 Murrah buffalo for 51,611 SNP; 350 of these had ROH longer than 10 Mb, indicating the occurrence of inbreeding in the last 5 generations. The mean length of the ROH per animal was 4.28 ± 1.85 Mb. Inbreeding coefficients were calculated from the genomic relationship matrix, the pedigree, and the ROH, with estimates varying between 0.242 and 0.035. Inbreeding estimates from the pedigree had a low correlation with the genomic estimates, and estimates from the genomic relationship matrix were much higher than those from the pedigree or the ROH. Signatures of selection were identified in 6 genomic regions, located on chromosomes 1, 2, 3, 5, 16, and 18, encompassing a total of 190 genes and 174 QTL. Many of the genes (e.g., APRT and ACSF3) and QTL identified are related to milk production traits, such as milk yield, milk fat yield and percentage, and milk protein yield and percentage. Other genes are associated with reproduction and immune response traits as well as morphological aspects of the buffalo species. Inbreeding levels in this population are still low but are increasing due to selection and should be managed to avoid future losses due to inbreeding depression. The proximity of genes linked to milk production traits with genes associated with reproduction and immune system traits suggests the need to include these latter genes in the breeding program to avoid negatively affecting them due to selection for production traits.     

A quantitative trait locus on Bos taurus autosome 17 explains a large proportion of the genetic variation in de novo synthesized milk fatty acids

A genomic region associated with milk fatty acid (FA) composition has been detected on Bos taurus autosome (BTA)17 based on 50,000 (50K) single nucleotide polymorphism (SNP) genotypes. The aim of our study was to fine-map BTA17 with imputed 777,000 (777K) SNP genotypes to identify candidate genes associated with milk FA composition. Phenotypes consisted of gas chromatography measurements of 14 FA based on winter and summer milk samples. Phenotypes and genotypes were available on 1,640 animals in winter milk, and on 1,581 animals in summer milk samples. Single-SNP analyses showed that several SNP in a region located between 29.0 and 34.0 Mbp were in strong association with C6:0, C8:0, and C10:0. This region was further characterized based on haplotypes. In summer milk samples, for example, these haplotypes explained almost 10% of the genetic variance in C6:0, 9% in C8:0, 3.5% in C10:0, 1.8% in C12:0, and 0.9% in C14:0. Two groups of haplotypes with distinct predicted effects could be defined, suggesting the presence of one causal variant. Predicted haplotype effects tended to increase from C6:0 to C14:0; however, the proportion of genetic variance explained by the haplotypes tended to decrease from C6:0 to C14:0. This is an indication that the quantitative trait locus (QTL) region is involved either in the elongation process or in early termination of de novo synthesized FA. Although many genes are present in this QTL region, most of these genes on BTA17 have not been characterized yet. The strongest association was found close to the progesterone receptor membrane component 2 (PGRMC2) gene, which has not yet been associated with milk FA composition. Therefore, no clear candidate gene associated with milk FA composition could be identified for this QTL.     

High-density genome-wide association study for residual feed intake in Holstein dairy cattle

Improving feed efficiency (FE) of dairy cattle may boost farm profitability and reduce the environmental footprint of the dairy industry. Residual feed intake (RFI), a candidate FE trait in dairy cattle, can be defined to be genetically uncorrelated with major energy sink traits (e.g., milk production, body weight) by including genomic predicted transmitting ability of such traits in genetic analyses for RFI. We examined the genetic basis of RFI through genome-wide association (GWA) analyses and post-GWA enrichment analyses and identified candidate genes and biological pathways associated with RFI in dairy cattle. Data were collected from 4,823 lactations of 3,947 Holstein cows in 9 research herds in the United States. Of these cows, 3,555 were genotyped and were imputed to a high-density list of 312,614 SNP. We used a single-step GWA method to combine information from genotyped and nongenotyped animals with phenotypes as well as their ancestors' information. The estimated genomic breeding values from a single-step genomic BLUP were back-solved to obtain the individual SNP effects for RFI. The proportion of genetic variance explained by each 5-SNP sliding window was also calculated for RFI. Our GWA analyses suggested that RFI is a highly polygenic trait regulated by many genes with small effects. The closest genes to the top SNP and sliding windows were associated with dry matter intake (DMI), RFI, energy homeostasis and energy balance regulation, digestion and metabolism of carbohydrates and proteins, immune regulation, leptin signaling, mitochondrial ATP activities, rumen development, skeletal muscle development, and spermatogenesis. The region of 40.7 to 41.5 Mb on BTA25 (UMD3.1 reference genome) was the top associated region for RFI. The closest genes to this region, CARD11 and EIF3B, were previously shown to be related to RFI of dairy cattle and FE of broilers, respectively. Another candidate region, 57.7 to 58.2 Mb on BTA18, which is associated with DMI and leptin signaling, was also associated with RFI in this study. Post-GWA enrichment analyses used a sum-based marker-set test based on 4 public annotation databases: Gene Ontology, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, Reactome pathways, and medical subject heading (MeSH) terms. Results of these analyses were consistent with those from the top GWA signals. Across the 4 databases, GWA signals for RFI were highly enriched in the biosynthesis and metabolism of amino acids and proteins, digestion and metabolism of carbohydrates, skeletal development, mitochondrial electron transport, immunity, rumen bacteria activities, and sperm motility. Our findings offer novel insight into the genetic basis of RFI and identify candidate regions and biological pathways associated with RFI in dairy cattle.     

Distribution and location of genetic effects for dairy traits

Genetic effects for many dairy traits and for total economic merit are evenly distributed across all chromosomes. A high-density scan using 38,416 single nucleotide polymorphism markers for 5,285 bulls confirmed 2 previously known major genes on Bos taurus autosomes (BTA) 6 and 14 but revealed few other large effects. Markers on BTA18 had the largest effects on calving ease, several conformation traits, longevity, and total merit. Prediction accuracy was highest using a heavy-tailed prior assuming that each marker had an effect on each trait, rather than assuming a normal distribution of effects as in a linear model, or that only some loci have nonzero effects. A prior model combining heavy tails with finite alleles produced results that were intermediate compared with the individual models. Differences between models were small (1 to 2%) for traits with no major genes and larger for heavy tails with traits having known quantitative trait loci (QTL; 6 to 8%). Analysis of bull recessive codes suggested that marker effects from genomic selection may be used to identify regions of chromosomes to search in detail for candidate genes, but individual single nucleotide polymorphisms were not tracking causative mutations with the exception of diacylglycerol O-acyltransferase 1. Additive genetic merits were constructed for each chromosome, and the distribution of BTA14-specific estimated breeding value (EBV) showed that selection primarily for milk yield has not changed the distribution of EBV for fat percentage even in the presence of a known QTL. Such chromosomal EBV also may be useful for identifying complementary mates in breeding programs. The QTL affecting dystocia, conformation, and economic merit on BTA18 appear to be related to calf size or birth weight and may be the result of longer gestation lengths. Results validate quantitative genetic assumptions that most traits are due to the contributions of a large number of genes of small additive effect, rather than support the finite locus model.     

Gene mapping,gene-set analysis,and genomic prediction of postpartum blood calcium in Holstein cows

The onset of lactation results in a sudden irreversible loss of Ca for colostrum and milk synthesis. Some cows are unable to quickly adapt to this demand and succumb to clinical hypocalcemia, whereas a larger proportion of cows develop subclinical hypocalcemia that predisposes them to other peripartum diseases. The objective of this study was to perform a comprehensive genomic analysis of blood total Ca concentration in periparturient Holstein cows. We first performed a genomic scan and a subsequent gene-set analysis to identify candidate genes, biological pathways, and molecular mechanisms affecting postpartum Ca concentration. Then, we assessed the prediction of postpartum Ca concentration using genomic information. Data consisted of 7,691 records of plasma or serum concentrations of Ca measured in the first, second, and third day after parturition of 959 primiparous and 1,615 multiparous cows that calved between December 2015 and June 2020 in 2 dairy herds. All cows were genotyped with 80k SNPs. The statistical model included lactation (1 to 5+), calf category (male, females, twins), and day as fixed effects, and season-treatment-experiment, animal, and permanent environmental as random effects. Model predictive ability was evaluated using 10-fold cross-validation. Heritability and repeatability estimates were 0.083 (standard error = 0.017) and 0.444 (standard error = 0.028). The association mapping identified 2 major regions located on Bos taurus autosome (BTA)6 and BTA16 that explained 1.2% and 0.7% of additive genetic variance of Ca concentration, respectively. Interestingly, the region on BTA6 harbors the GC gene, which encodes the vitamin D binding protein, and the region on BTA16 harbors LRRC38, which is actively involved in K transport. Other sizable peaks were identified on BTA5, BTA2, BTA7, BTA14, and BTA9. These regions harbor genes associated with Ca channels (CACNA1S, CRACR2A), K channels (KCNK9), bone remodeling (LRP6), and milk production (SOCS2). The gene-set analysis revealed terms related to vitamin transport, calcium ion transport, calcium ion binding, and calcium signaling. Genomic predictions of phenotypic and genomic estimated breeding values of Ca concentration yielded predictive correlations up to 0.50 and 0.15, respectively. Overall, the present study contributes to a better understanding of the genetic basis of postpartum blood Ca concentration in Holstein cows. In addition, the findings may contribute to the development of novel selection and management strategies for reducing periparturient hypocalcemia in dairy cattle.