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

Genomic insights into growth traits in German Black Pied cattle: a dual-purpose breed at risk

The German Black Pied cattle (DSN) is an endangered dual-purpose breed valued for its genetic diversity and high milk fat and protein content. However, due to competition with higher-yielding dairy breeds, the DSN population has declined, leading to its designation as an endangered breed. While previous research has focused on the milk production traits of DSN, this study aims to address meat traits to further understand the genetic determination of the dual-purpose characteristics of the breed. We conducted genome-wide association studies on 669 DSN bulls to identify genetic loci associated with birth weight, BW, and BW gain at different growth stages. Using imputed whole-genome sequencing data, we identified 14 quantitative trait loci across ten chromosomes. Significant associations were found for birth weight on chromosomes 5 and 18, for body weight at 3 weeks (BW3w) on chromosomes 3 and 16, for body weight at 7 months (BW7m) on chromosomes 3 and 10, and for body weight gain from birth or 3 weeks to 18 months (BWG0d-18m, BWG3w-18m) on chromosomes 4 and 7. Key positional candidate genes influencing muscle and fat tissue development included RERGL and LMO3 (identified for birth weight), MET and CAPZA2 (identified for BWG0d-18m) which are essential for skeletal muscle development and actin filament regulation, respectively, TLN2 (identified for BW7m), MYO1F and ADAMTS10 (identified for BWG3w-18m) which are critical for actin filament assembly, cytoskeletal function, and skeletal development, respectively. Candidate genes such as CPT2 (identified or BW3w) and VPS13C (identified or BW7m) are involved in lipid metabolism and mitochondrial function. Additionally, candidate genes such as IGSF3 (identified for BW7m), KLRC1 and members of the C-type lectin family (identified for birth weight) are associated with immune regulation, and thus, suggest a potential interplay between metabolism, immune function, and growth efficiency. These findings highlight the distinct genetic mechanisms underlying growth at various developmental stages, underscoring the importance of breed-specific genetic evaluations. The identified loci also overlap with previously reported loci for meat and production traits in other cattle breeds, underscoring their relevance and potential utility in DSN breeding strategies. This study provides a foundation for conservation and genomic breeding strategies to maintain the dual-purpose characteristics of DSN through optimising both meat and milk production.

Genome-wide association studies for milk production traits in two autochthonous Aosta cattle breeds

Genome-wide association studies (GWASs) are used to identify quantitative trait loci for phenotypic traits of interest. The use of multilocus mixed models allows to correct for population stratification and account for long-range linkage disequilibrium. In this study, GWASs were conducted to identify the genetic bases of milk production (milk yield, protein and fat composition, and yield) in two autochthonous dual-purpose cattle breeds from the Aosta Valley. Using either the breeding values or the deregressed proofs, common significative single nucleotide polymorphisms have been identified for milk yield, protein percentage, and fat percentage. Two major quantitative trait loci regions have been identified on the chromosomes 5 and 14 for the fat percentage, harbouring the MGST1, CYHR1, VPS28, and CPSF1 genes. For the protein percentage, a candidate region has been identified on BTA 6; in this region, the CSN1S1, CSN2, HSTN, CSN3, and RUFY3 genes are annotated. Most of the identified genes have already been associated with milk composition in other studies on cosmopolitan and local cattle. These results show that the genes involved in milk composition quantitative traits in the Aosta cattle are common also in other cattle breeds and they can be further investigated with the use of whole genome sequencing data.

Improving the accuracy of genomic prediction in dairy cattle using the biologically annotated neural networks framework

BACKGROUND

Biologically annotated neural networks (BANNs) are feedforward Bayesian neural network models that utilize partially connected architectures based on SNP-set annotations. As an interpretable neural network, BANNs model SNP and SNP-set effects in their input and hidden layers, respectively. Furthermore, the weights and connections of the network are regarded as random variables with prior distributions reflecting the manifestation of genetic effects at various genomic scales. However, its application in genomic prediction has yet to be explored.

RESULTS

This study extended the BANNs framework to the area of genomic selection and explored the optimal SNP-set partitioning strategies by using dairy cattle datasets. The SNP-sets were partitioned based on two strategies-gene annotations and 100 kb windows, denoted as BANN_gene and BANN_100kb, respectively. The BANNs model was compared with GBLUP, random forest (RF), BayesB and BayesCπ through five replicates of five-fold cross-validation using genotypic and phenotypic data on milk production traits, type traits, and one health trait of 6,558, 6,210 and 5,962 Chinese Holsteins, respectively. Results showed that the BANNs framework achieves higher genomic prediction accuracy compared to GBLUP, RF and Bayesian methods. Specifically, the BANN_100kb demonstrated superior accuracy and the BANN_gene exhibited generally suboptimal accuracy compared to GBLUP, RF, BayesB and BayesCπ across all traits. The average accuracy improvements of BANN_100kb over GBLUP, RF, BayesB and BayesCπ were 4.86%, 3.95%, 3.84% and 1.92%, and the accuracy of BANN_gene was improved by 3.75%, 2.86%, 2.73% and 0.85% compared to GBLUP, RF, BayesB and BayesCπ, respectively across all seven traits. Meanwhile, both BANN_100kb and BANN_gene yielded lower overall mean square error values than GBLUP, RF and Bayesian methods.

CONCLUSION

Our findings demonstrated that the BANNs framework performed better than traditional genomic prediction methods in our tested scenarios, and might serve as a promising alternative approach for genomic prediction in dairy cattle.

Unraveling genomic diversity and positive selection signatures of Qaidam cattle through whole-genome re-sequencing

Qaidam cattle are a typical Chinese native breed inhabiting northwest China. They bear the characteristics of high cold and roughage tolerance, low-oxygen adaptability and good meat quality. To analyze the genetic diversity of Qaidam cattle, 60 samples were sequenced using whole-genome resequencing technology, along with 192 published sets of whole-genome sequencing data of Indian indicine cattle, Chinese indicine cattle, North Chinese cattle breeds, East Asian taurine cattle, Eurasian taurine cattle and European taurine cattle as controls. It was found that Qaidam cattle have rich genetic diversity in Bos taurus, but the degree of inbreeding is also high, which needs further protection. The phylogenetic analysis, principal component analysis and ancestral component analysis showed that Qaidam cattle mainly originated from East Asian taurine cattle. Qaidam cattle had a closer genetic relationship with the North Chinese cattle breeds and the least differentiation from Mongolian cattle. Annotating the selection signals obtained by composite likelihood ratio, nucleotide diversity analysis, integrated haplotype score, genetic differentiation index, genetic diversity ratio and cross-population extended haplotype homozygosity methods, several genes associated with immunity, reproduction, meat, milk, growth and adaptation showed strong selection signals. In general, this study provides genetic evidence for understanding the germplasm characteristics of Qaidam cattle. At the same time, it lays a foundation for the scientific and reasonable protection and utilization of genetic resources of Chinese local cattle breeds, which has great theoretical and practical significance.

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.

Genome-Wide Association Study of Lactation Traits in Chinese Holstein Cows in Southern China

Lactation traits are economically important for dairy cows. Southern China has a high-temperature and high-humidity climate, and environmental and genetic interactions greatly impact dairy cattle performance. The aim of this study was to identify novel single-nucleotide polymorphism sites and novel candidate genes associated with lactation traits in Chinese Holstein cows under high-temperature and humidity conditions in southern China. A genome-wide association study was performed for the lactation traits of 392 Chinese Holstein cows, using GGP Bovine 100 K SNP gene chips. Some 23 single nucleotide polymorphic loci significantly associated with lactation traits were screened. Among them, 16 were associated with milk fat rate, 7 with milk protein rate, and 3 with heat stress. A quantitative trait locus that significantly affects milk fat percentage in Chinese Holstein cows was identified within a window of approximately 0.5 Mb in the region of 0.4-0.9 Mb on Bos taurus autosome 14. According to Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses, ten genes (DGAT1, IDH2, CYP11B1, GFUS, CYC1, GPT, PYCR3, OPLAH, ALDH1A3, and NAPRT) associated with lactation fat percentage, milk yield, antioxidant activity, stress resistance, and inflammation and immune response were identified as key candidates for lactation traits. The results of this study will help in the development of an effective selection and breeding program for Chinese Holstein cows in high-temperature and humidity regions.

Assessing Genomic Diversity and Signatures of Selection in Chinese Red Steppe Cattle Using High-Density SNP Array

Chinese Red Steppe Cattle (CRS), a composite cattle breed, is well known for its milk production, high slaughter rate, carcass traits, and meat quality. Nowadays, it is widely bred in Jilin and Hebei Province and the Inner Mongolia Autonomous region. However, the population structure and the genetic basis of prominent characteristics of CRS are still unknown. In this study, we systematically describe their population structure, genetic diversity, and selection signature based on genotyping data from 61 CRS individuals with GGP Bovine 100 K chip. The results showed that CRS cattle had low inbreeding levels and had formed a unique genetic structure feature. Using two complementary methods (including comprehensive haplotype score and complex likelihood ratio), we identified 1291 and 1285 potentially selected genes, respectively. There were 141 genes annotated in common 106 overlapping genomic regions covered 5.62 Mb, including PLAG1, PRKG2, DGAT1, PARP10, TONSL, ADCK5, and BMP3, most of which were enriched in pathways related to muscle growth and differentiation, milk production, and lipid metabolism. This study will contribute to understanding the genetic mechanism behind artificial selection and give an extensive reference for subsequent breeding.

Longitudinal genome-wide association studies of milk production traits in Holstein cattle using whole-genome sequence data imputed from medium-density chip data

Longitudinal traits, such as milk production traits in dairy cattle, are featured by having phenotypic values at multiple time points, which change dynamically over time. In this study, we first imputed SNP chip (50-100K) data to whole-genome sequence (WGS) data in a Chinese Holstein population consisting of 6,470 cows. The imputation accuracies were 0.88 to 0.97 on average after quality control. We then performed longitudinal GWAS in this population based on a random regression test-day model using the imputed WGS data. The longitudinal GWAS revealed 16, 39, and 75 quantitative trait locus regions associated with milk yield, fat percentage, and protein percentage, respectively. We estimated the 95% confidence intervals (CI) for these quantitative trait locus regions using the logP drop method and identified 581 genes involved in these CI. Further, we focused on the CI that covered or overlapped with only 1 gene or the CI that contained an extremely significant top SNP. Twenty-eight candidate genes were identified in these CI. Most of them have been reported in the literature to be associated with milk production traits, such as DGAT1, HSF1, MGST1, GHR, ABCG2, ADCK5, and CSN1S1. Among the unreported novel genes, some also showed good potential as candidate genes, such as CCSER1, CUX2, SNTB1, RGS7, OSR2, and STK3, and are worth being further investigated. Our study provided not only new insights into the candidate genes for milk production traits, but also a general framework for longitudinal GWAS based on random regression test-day model using WGS data.

Genome-wide association analysis of milk production, somatic cell score, and body conformation traits in Holstein cows

Milk production and body conformation traits are critical economic traits for dairy cows. To understand the basic genetic structure for those traits, a genome wide association study was performed on milk yield, milk fat yield, milk fat percentage, milk protein yield, milk protein percentage, somatic cell score, body form composite index, daily capacity composite index, feed, and leg conformation traits, based on the Illumina Bovine HD100k BeadChip. A total of 57, 12 and 26 SNPs were found to be related to the milk production, somatic cell score and body conformation traits in the Holstein cattle. Genes with pleiotropic effect were also found in this study. Seven significant SNPs were associated with multi-traits and were located on the PLEC, PLEKHA5, TONSL, PTGER4, and LCORL genes. In addition, some important candidate genes, like GPAT3, CEBPB, AGO2, SLC37A1, and FNDC3B, were found to participate in fat metabolism or mammary gland development. These results can be used as candidate genes for milk production, somatic cell score, and body conformation traits of Holstein cows, and are helpful for further gene function analysis to improve milk production and quality.

Diversity of Mitochondrial DNA Haplogroups and Their Association with Bovine Antral Follicle Count

Maternal origins based on the bovine mitochondrial D-loop region are proven to have two main origins: Bos taurus and Bos indicus. To examine the association between the maternal origins of bovine and reproductive traits, the complete mitochondrial D-loop region sequences from 501 Chinese Holstein cows and 94 individuals of other breeds were analyzed. Based on the results obtained from the haplotype analysis, 260 SNPs (single nucleotide polymorphism), 32 indels (insertion/deletion), and 219 haplotypes were identified. Moreover, the nucleotide diversity (π) and haplotype diversity (Hd) were 0.024 ± 0.001 and 0.9794 ± 0.003, respectively, indicating the abundance of genetic resources in Chinese Holstein cows. The results of the median-joining network analysis showed two haplogroups (HG, including HG1 and HG2) that diverged in genetic distance. Furthermore, the two haplogroups were significantly (p < 0.05) correlated with the antral follicle (diameter ≥ 8 mm) count, and HG1 individuals had more antral follicles than HG2 individuals, suggesting that these different genetic variants between HG1 and HG2 correlate with reproductive traits. The construction of a neighbor-joining phylogenetic tree and principal component analysis also revealed two main clades (HG1 and HG2) with different maternal origins: Bos indicus and Bos taurus, respectively. Therefore, HG1 originating from the maternal ancestors of Bos indicus may have a greater reproductive performance, and potential genetic variants discovered may promote the breeding process in the cattle industry.

Genetic Parameter Estimation and Genome-Wide Association Study-Based Loci Identification of Milk-Related Traits in Chinese Holstein

Accurately estimating the genetic parameters and revealing more genetic variants underlying milk production and quality are conducive to the genetic improvement of dairy cows. In this study, we estimate the genetic parameters of five milk-related traits of cows-namely, milk yield (MY), milk fat percentage (MFP), milk fat yield (MFY), milk protein percentage (MPP), and milk protein yield (MPY)-based on a random regression test-day model. A total of 95,375 test-day records of 9,834 cows in the lower reaches of the Yangtze River were used for the estimation. In addition, genome-wide association studies (GWASs) for these traits were conducted, based on adjusted phenotypes. The heritability, as well as the standard errors, of MY, MFP, MFY, MPP, and MPY during lactation ranged from 0.22 ± 0.02 to 0.31 ± 0.04, 0.06 ± 0.02 to 0.15 ± 0.03, 0.09 ± 0.02 to 0.28 ± 0.04, 0.07 ± 0.01 to 0.16 ± 0.03, and 0.14 ± 0.02 to 0.27 ± 0.03, respectively, and the genetic correlations between different days in milk (DIM) within lactations decreased as the time interval increased. Two, six, four, six, and three single nucleotide polymorphisms (SNPs) were detected, which explained 5.44, 12.39, 8.89, 10.65, and 7.09% of the phenotypic variation in MY, MFP, MFY, MPP, and MPY, respectively. Ten Kyoto Encyclopedia of Genes and Genomes pathways and 25 Gene Ontology terms were enriched by analyzing the nearest genes and genes within 200 kb of the detected SNPs. Moreover, 17 genes in the enrichment results that may play roles in milk production and quality were selected as candidates, including CAMK2G, WNT3A, WNT9A, PLCB4, SMAD9, PLA2G4A, ARF1, OPLAH, MGST1, CLIP1, DGAT1, PRMT6, VPS28, HSF1, MAF1, TMEM98, and F7. We hope that this study will provide useful information for in-depth understanding of the genetic architecture of milk production and quality traits, as well as contribute to the genomic selection work of dairy cows in the lower reaches of the Yangtze River.

Genome-Wide Association Study Candidate Genes on Mammary System-Related Teat-Shape Conformation Traits in Chinese Holstein Cattle

In the dairy industry, mammary system traits are economically important for dairy animals, and it is important to explain their fundamental genetic architecture in Holstein cattle. Good and stable mammary system-related teat traits are essential for producer profitability in animal fitness and in the safety of dairy production. In this study, we conducted a genome-wide association study on three traits—anterior teat position (ATP), posterior teat position (PTP), and front teat length (FTL)—in which the FarmCPU method was used for association analyses. Phenotypic data were collected from 1000 Chinese Holstein cattle, and the GeneSeek Genomic Profiler Bovine 100K single-nucleotide polymorphisms (SNP) chip was used for cattle genotyping data. After the quality control process, 984 individual cattle and 84,406 SNPs remained for GWAS work analysis. Nine SNPs were detected significantly associated with mammary-system-related teat traits after a Bonferroni correction (p < 5.92 × 10⁻⁷), and genes within a region of 200 kb upstream or downstream of these SNPs were performed bioinformatics analysis. A total of 36 gene ontology (GO) terms and 3 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were significantly enriched (p < 0.05), and these terms and pathways are mainly related to metabolic processes, immune response, and cellular and amino acid catabolic processes. Eleven genes including MMS22L, E2F8, CSRP3, CDH11, PEX26, HAL, TAMM41, HIVEP3, SBF2, MYO16 and STXBP6 were selected as candidate genes that might play roles in the teat traits of cows. These results identify SNPs and candidate genes that give helpful biological information for the genetic architecture of these teat traits, thus contributing to the dairy production, health, and genetic selection of Chinese Holstein cattle.

Genomewide Association Analyses of Lactation Persistency and Milk Production Traits in Holstein Cattle Based on Imputed Whole-Genome Sequence Data

Lactation persistency and milk production are among the most economically important traits in the dairy industry. In this study, we explored the association of over 6.1 million imputed whole-genome sequence variants with lactation persistency (LP), milk yield (MILK), fat yield (FAT), fat percentage (FAT%), protein yield (PROT), and protein percentage (PROT%) in North American Holstein cattle. We identified 49, 3991, 2607, 4459, 805, and 5519 SNPs significantly associated with LP, MILK, FAT, FAT%, PROT, and PROT%, respectively. Various known associations were confirmed while several novel candidate genes were also revealed, including ARHGAP35, NPAS1, TMEM160, ZC3H4, SAE1, ZMIZ1, PPIF, LDB2, ABI3, SERPINB6, and SERPINB9 for LP; NIM1K, ZNF131, GABRG1, GABRA2, DCHS1, and SPIDR for MILK; NR6A1, OLFML2A, EXT2, POLD1, GOT1, and ETV6 for FAT; DPP6, LRRC26, and the KCN gene family for FAT%; CDC14A, RTCA, HSTN, and ODAM for PROT; and HERC3, HERC5, LALBA, CCL28, and NEURL1 for PROT%. Most of these genes are involved in relevant gene ontology (GO) terms such as fatty acid homeostasis, transporter regulator activity, response to progesterone and estradiol, response to steroid hormones, and lactation. The significant genomic regions found contribute to a better understanding of the molecular mechanisms related to LP and milk production in North American Holstein cattle.

Genome-Wide Association Study Based on Random Regression Model Reveals Candidate Genes Associated with Longitudinal Data in Chinese Simmental Beef Cattle

Simple Summary

Genome-wide association study (GWAS) has become the main approach for detecting functional genes that affects complex traits. For growth traits, the conventional GWAS method can only deal with the single-record traits observed at specific time points, rather than the longitudinal traits measured at multiple time points. Previous studies have reported the random regression model (RRM) for longitudinal data could overcome the limitation of the traditional GWAS model. Here, we present an association analysis based on RRM (GWAS-RRM) for 808 Chinese Simmental beef cattle at four stages of age. Ultimately, 37 significant single-nucleotide polymorphisms (SNPs) and several important candidate genes were screened to be associated with the body weight. Enrichment analysis showed these genes were significantly enriched in the signaling transduction pathway and lipid metabolism. This study not only offers a further understanding of the genetic basis for growth traits in beef cattle, but also provides a robust analytics tool for longitudinal traits in various species.

Abstract

Body weight (BW) is an important longitudinal trait that directly described the growth gain of bovine in production. However, previous genome-wide association study (GWAS) mainly focused on the single-record traits, with less attention paid to longitudinal traits. Compared with traditional GWAS models, the association studies based on the random regression model (GWAS-RRM) have better performance in the control of the false positive rate through considering time-stage effects. In this study, the BW trait data were collected from 808 Chinese Simmental beef cattle aged 0, 6, 12, and 18 months, then we performed a GWAS-RRM to fit the time-varied SNP effect. The results showed a total of 37 significant SNPs were associated with BW. Gene functional annotation and enrichment analysis indicated FGF4, ANGPT4, PLA2G4A, and ITGA5 were promising candidate genes for BW. Moreover, these genes were significantly enriched in the signaling transduction pathway and lipid metabolism. These findings will provide prior molecular information for bovine gene-based selection, as well as facilitate the extensive application of GWAS-RRM in domestic animals.

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.

Sequence-based genome-wide association study of individual milk mid-infrared wavenumbers in mixed-breed dairy cattle

BACKGROUND

Fourier-transform mid-infrared (FT-MIR) spectroscopy provides a high-throughput and inexpensive method for predicting milk composition and other novel traits from milk samples. While there have been many genome-wide association studies (GWAS) conducted on FT-MIR predicted traits, there have been few GWAS for individual FT-MIR wavenumbers. Using imputed whole-genome sequence for 38,085 mixed-breed New Zealand dairy cattle, we conducted GWAS on 895 individual FT-MIR wavenumber phenotypes, and assessed the value of these direct phenotypes for identifying candidate causal genes and variants, and improving our understanding of the physico-chemical properties of milk.

RESULTS

Separate GWAS conducted for each of 895 individual FT-MIR wavenumber phenotypes, identified 450 1-Mbp genomic regions with significant FT-MIR wavenumber QTL, compared to 246 1-Mbp genomic regions with QTL identified for FT-MIR predicted milk composition traits. Use of mammary RNA-seq data and gene annotation information identified 38 co-localized and co-segregating expression QTL (eQTL), and 31 protein-sequence mutations for FT-MIR wavenumber phenotypes, the latter including a null mutation in the ABO gene that has a potential role in changing milk oligosaccharide profiles. For the candidate causative genes implicated in these analyses, we examined the strength of association between relevant loci and each wavenumber across the mid-infrared spectrum. This revealed shared association patterns for groups of genomically-distant loci, highlighting clusters of loci linked through their biological roles in lactation and their presumed impacts on the chemical composition of milk.

CONCLUSIONS

This study demonstrates the utility of FT-MIR wavenumber phenotypes for improving our understanding of milk composition, presenting a larger number of QTL and putative causative genes and variants than found from FT-MIR predicted composition traits. Examining patterns of significance across the mid-infrared spectrum for loci of interest further highlighted commonalities of association, which likely reflects the physico-chemical properties of milk constituents.

Selection signatures in two oldest Russian native cattle breeds revealed using high-density single nucleotide polymorphism analysis

Native cattle breeds can carry specific signatures of selection reflecting their adaptation to the local environmental conditions and response to the breeding strategy used. In this study, we comprehensively analysed high-density single nucleotide polymorphism (SNP) genotypes to characterise the population structure and detect the selection signatures in Russian native Yaroslavl and Kholmogor dairy cattle breeds, which have been little influenced by introgression with transboundary breeds. Fifty-six samples of pedigree-recorded purebred animals, originating from different breeding farms and representing different sire lines, of the two studied breeds were genotyped using a genome-wide bovine genotyping array (Bovine HD BeadChip). Three statistical analyses—calculation of fixation index (F_ST) for each SNP for the comparison of the pairs of breeds, hapFLK analysis, and estimation of the runs of homozygosity (ROH) islands shared in more than 50% of animals—were combined for detecting the selection signatures in the genome of the studied cattle breeds. We confirmed nine and six known regions under putative selection in the genomes of Yaroslavl and Kholmogor cattle, respectively; the flanking positions of most of these regions were elucidated. Only two of the selected regions (localised on BTA 14 at 24.4–25.1 Mbp and on BTA 16 at 42.5–43.5 Mb) overlapped in Yaroslavl, Kholmogor and Holstein breeds. In addition, we detected three novel selection sweeps in the genome of Yaroslavl (BTA 4 at 4.74–5.36 Mbp, BTA 15 at 17.80–18.77 Mbp, and BTA 17 at 45.59–45.61 Mbp) and Kholmogor breeds (BTA 12 at 82.40–81.69 Mbp, BTA 15 at 16.04–16.62 Mbp, and BTA 18 at 0.19–1.46 Mbp) by using at least two of the above-mentioned methods. We expanded the list of candidate genes associated with the selected genomic regions and performed their functional annotation. We discussed the possible involvement of the identified candidate genes in artificial selection in connection with the origin and development of the breeds. Our findings on the Yaroslavl and Kholmogor breeds obtained using high-density SNP genotyping and three different statistical methods allowed the detection of novel putative genomic regions and candidate genes that might be under selection. These results might be useful for the sustainable development and conservation of these two oldest Russian native cattle breeds.

Genome-wide association for milk production and lactation curve parameters in Holstein dairy cows

The aim of this study was to identify genomic regions associated with 305‐day milk yield and lactation curve parameters on primiparous (n = 9,910) and multiparous (n = 11,158) Holstein cows. The SNP solutions were estimated using a weighted single‐step genomic BLUP approach and imputed high‐density panel (777k) genotypes. The proportion of genetic variance explained by windows of 50 consecutive SNP (with an average of 165 Kb) was calculated, and regions that accounted for more than 0.50% of the variance were used to search for candidate genes. Estimated heritabilities were 0.37, 0.34, 0.17, 0.12, 0.30 and 0.19, respectively, for 305‐day milk yield, peak yield, peak time, ramp, scale and decay for primiparous cows. Genetic correlations of 305‐day milk yield with peak yield, peak time, ramp, scale and decay in primiparous cows were 0.99, 0.63, 0.20, 0.97 and −0.52, respectively. The results identified three windows on BTA14 associated with 305‐day milk yield and the parameters of lactation curve in primi‐ and multiparous cows. Previously proposed candidate genes for milk yield supported by this work include GRINA, CYHR1, FOXH1, TONSL, PPP1R16A, ARHGAP39, MAF1, OPLAH and MROH1, whereas newly identified candidate genes are MIR2308, ZNF7, ZNF34, SLURP1, MAFA and KIFC2 (BTA14). The protein lipidation biological process term, which plays a key role in controlling protein localization and function, was identified as the most important term enriched by the identified genes.

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.

Identification of genetic markers associated with milk production traits in Chinese Holstein cattle based on post genome-wide association studies

With the rapid development of dairy industry, the breeding process of dairy cows has been accelerated. In previous genome-wide association studies (GWAS), a large number of genetic markers have been reported which may contribute to the selection of Holstein populations with superior milk-producing traits, but they remain to be further verified before practical application. In this study, 90 single nucleotide polymorphisms (SNPs) were selected, which were reported to be significantly associated with five milk production traits, including 305-day milk yield (305MY), 305-day milk fat percent (305FC), 305-day milk protein percent (305PC), 305-day milk fat yield (305FY) and 305-day milk protein yield (305PY). Effective 305-day data and fresh DNA samples were obtained from 295 healthy cows with gestational age of 1-4. Matrix-assisted laser desorption/ionization-time of flight (MALDI-TOF) was used to perform precise genotyping of these loci, followed by site association and haplotype analysis. Results showed that 36 out of 90 loci were supported to be used as genetic markers. In particular, several novel and effective haplotypes were also presented. Overall, our results verified tens of useful markers and provided a basis for further development of breeding strategies.

Genetic Effects of LPIN1 Polymorphisms on Milk Production Traits in Dairy Cattle

Our initial RNA sequencing work identified that lipin 1 (LPIN1) was differentially expressed during dry period, early lactation, and peak of lactation in dairy cows, and it was enriched into the fat metabolic Gene Ontology (GO) terms and pathways, thus we considered LPIN1 as the candidate gene for milk production traits. In this study, we detected the polymorphisms of LPIN1 and verified their genetic effects on milk yield and composition in a Chinese Holstein cow population. We found seven SNPs by re-sequencing the entire coding region and partial flanking region of LPIN1, including one in 5′ flanking region, four in exons, and two in 3′ flanking region. Of these, four SNPs, c.637T > C, c.708A > G, c.1521C > T, and c.1555A > C, in the exons were predicted to result in the amino acid replacements. With the Haploview 4.2, we found that seven SNPs in LPIN1 formed two haplotype blocks (D′ = 0.98–1.00). Single-SNP association analyses showed that SNPs were significantly associated with milk yield, fat yield, fat percentage, or protein yield in the first or second lactation (p = < 0.0001–0.0457), and only g.86049389C > T was strongly associated with protein percentage in both lactations (p = 0.0144 and 0.0237). The haplotype-based association analyses showed that the two haplotype blocks were significantly associated with milk yield, fat yield, protein yield, or protein percentage (p = < 0.0001–0.0383). By quantitative real-time PCR (qRT-PCR), we found that LPIN1 had relatively high expression in mammary gland and liver tissues. Furthermore, we predicted three SNPs, c.637T > C, c.708A > G, and c.1521C > T, using SOPMA software, changing the LPIN1 protein structure that might be potential functional mutations. In summary, we demonstrated the significant genetic effects of LPIN1 on milk production traits, and the identified SNPs could serve as genetic markers for dairy breeding.