Short communication: Confirmation of candidate causative variants on milk composition and cheesemaking properties in Montbéliarde cows

In a previous study, we identified candidate causative variants located in 24 functional candidate genes for milk protein and fatty acid composition in Montbéliarde, Normande, and Holstein cows. We designed these variants on the custom part of the EuroG10K BeadChip (Illumina Inc., San Diego, CA), which is routinely used for genomic selection analyses in French dairy cattle. To validate the effects of these candidate variants on milk composition and to estimate their effects on cheesemaking properties, a genome-wide association study was performed on milk protein, fatty acid and mineral composition, as well as on 9 cheesemaking traits (3 laboratory cheese yields, 5 coagulation traits, and milk pH). All the traits were predicted from midinfrared spectra in the Montbéliarde cow population of the Franche-Comté region. A total of 194 candidate variants located in 24 genes and 17 genomic regions were imputed on 19,862 cows with phenotypes and genotyped with either the BovineSNP50 (Illumina Inc.) or the EuroG10K BeadChip. We then tested the effect of each SNP in a mixed linear model including random polygenic effects estimated with a genomic relationship matrix. We confirm here the effects of candidate causative variants located in 17 functional candidate genes on both cheesemaking properties and milk composition traits. In each candidate gene, we identified the most plausible causative variant: 4 are missense in the ALPL, SLC26A4, CSN3, and SCD genes, 7 are located in 5′UTR (AGPAT6), 3′ untranslated region (GPT), or upstream (CSN1S1, CSN1S2, PAEP, DGAT1, and PICALM) regions, and 6 are located in introns of the SLC37A1, MGST1, CSN2, BRI3BP, FASN, and ANKH genes.     

Genomic analyses and biological validation of candidate genes for rectal temperature as an indicator of heat stress in Holstein cattle

Heat stress is a major cause of welfare issues and economic losses to the worldwide dairy cattle industry. Genetic selection for heat tolerance has a great potential to positively affect the dairy industry, as the gains are permanent and cumulative over generations. Rectal temperature (RT) is hypothesized to be a good indicator trait of heat tolerance. Therefore, this study investigated the genetic architecture of RT by estimating genetic parameters, performing genome-wide association studies, and biologically validating potential candidate genes identified to be related to RT in Holstein cattle. A total of 33,013 RT records from 7,598 cows were used in this study. In addition, 1,114 cows were genotyped using the Illumina 150K Bovine BeadChip (Illumina, San Diego, CA). Rectal temperature measurements taken in the morning (AMRT) and in the afternoon (PMRT) are moderately heritable traits, with estimates of 0.09 ± 0.02 and 0.04 ± 0.01, respectively. These 2 traits are also highly genetically correlated (r = 0.90 ± 0.08). A total of 10 SNPs (located on BTA3, BTA4, BTA8, BTA13, BTA14, and BTA29) were found to be significantly associated with AMRT and PMRT. Subsequently, gene expression analyses were performed to validate the key functional genes identified (SPAG17, FAM107B, TSNARE1, RALYL, and PHRF1). This was done through in vitro exposure of peripheral blood mononuclear cells (PBMC) to different temperatures (37°C, 39°C, and 42°C). The relative mRNA expression of 2 genes, FAM107B and PHRF1, significantly changed between the control and heat stressed PBMC. In summary, RT is heritable, and enough genetic variability exists to enable genetic improvement of heat tolerance in Holstein cattle. Important genomic regions were identified and biologically validated; FAM107B and PHRF1 are the main candidate genes identified to influence heat stress response in dairy cattle.     

Non-culture-based identification of mastitis-causing bacteria by MALDI-TOF mass spectrometry

The purpose of this study was to evaluate the detection limit of matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) for direct identification, without previous microbiological culture, of bovine mastitis-causing bacteria from milk samples. Milk samples (n = 15) were experimentally contaminated with Staphylococcus aureus, Streptococcus uberis, Streptococcus agalactiae, Streptococcus dysgalactiae, and Escherichia coli to have bacterial counts ranging from 10³ to 10⁹ cfu/mL. These contaminated milk samples were subjected to a preparation protocol for bacterial ribosomal protein extraction using the MALDI Sepsityper kit (Bruker Daltonik, Bremen, Germany), which allowed MALDI-TOF MS coupled with Biotyper software (Bruker Daltonik) to identify bacterial fingerprints based on intact ribosomal proteins. The ability of MALDI-TOF MS to correctly identify bacterial strains from experimentally contaminated milk (without previous microbiological culture) depended on the bacterial count of the samples and on the species of the bacteria evaluated. Adequate identification at the bacterial species level (score ≥2.0) directly from milk samples required bacterial counts in the following ranges: ≥10⁶ cfu/mL of Staph. aureus, ≥10⁷ cfu/mL of E. coli, and ≥10⁸ cfu/mL of Strep. agalactiae, Strep. dysgalactiae, and Strep. uberis. We concluded that direct identification of mastitis-causing pathogens is possible for Staph. aureus, E. coli, Strep. agalactiae, Strep. dysgalactiae, and Strep. uberis, but correct identification depended on the bacterial count in the milk samples.     

Elimination of selected mastitis pathogens during the dry period

We aimed to evaluate the elimination of 4 different mastitis pathogens, Streptococcus agalactiae, Mycoplasma bovis, Staphylococcus aureus, and Streptococcus uberis, from infected udder quarters during the dry period using quantitative PCR. The second purpose of this study was to evaluate the association between milk haptoglobin (Hp) concentration and the presence of udder pathogens (Strep. agalactiae, Staph. aureus, M. bovis, and Strep. uberis) in udder quarter milk samples before and after dry period. Aseptic udder quarter milk samples (n = 1,001) were collected from 133 dairy cows at dry off and at the first milking after calving from 1 large dairy herd. Bacterial DNA of Strep. agalactiae, Staph. aureus, Strep. uberis, and M. bovis in the udder quarter milk samples was identified with commercial quantitative PCR analysis Mastitis 4B (DNA Diagnostic A/S, Risskov, Denmark). Milk Hp concentration (mg/L) was measured from udder quarter milk samples. The elimination rates during the dry period for M. bovis, Staph. aureus, Strep. agalactiae, and Strep. uberis were 86.7, 93.6, 96.2, and 100.0%, respectively. The new IMI rate was 3.0% for M. bovis, 2.9% for Staph. aureus, 2.4% for Strep. agalactiae, and 3.1% for Strep. uberis. The milk Hp concentration was significantly higher in udder quarter milk samples with blood and in samples positive for Strep. agalactiae at dry off and for Staph. aureus postcalving. Elevated milk Hp concentration was not associated with the presence of M. bovis in the udder quarter milk samples. In conclusion, elimination of Staph. aureus, Strep. agalactiae, and Strep. uberis during the dry period was high; the elimination of M. bovis from infected udder quarters was lower, but probably spontaneous. Additionally, milk Hp concentration may be used as a marker for udder inflammation when combined with the bacteriological results at dry off and postpartum.     

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.     

Role of CSN2, CSN3, and BLG genes and the polygenic background in the cattle milk protein profile

To devise better selection strategies in dairy cattle breeding programs, a deeper knowledge of the role of the major genes encoding for milk protein fractions is required. The aim of the present study was to assess the effect of the CSN2, CSN3, and BLG genotypes on individual protein fractions (α_S1-CN, α_S2-CN, β-CN, κ-CN, β-LG, α-LA) expressed qualitatively as percentages of total nitrogen content (% N), quantitatively as contents in milk (g/L), and as daily production levels (g/d). Individual milk samples were collected from 1,264 Brown Swiss cows reared in 85 commercial herds in Trento Province (northeast Italy). A total of 989 cows were successfully genotyped using the Illumina Bovine SNP50 v.2 BeadChip (Illumina Inc.), and a genomic relationship matrix was constructed using the 37,519 SNP markers obtained. Milk protein fractions were quantified and the β-CN, κ-CN, and β-LG genetic variants were identified by reversed-phase HPLC (RP-HPLC). All protein fractions were analyzed through a Bayesian multitrait animal model implemented via Gibbs sampling. The effects of days in milk, parity order, and the CSN2, CSN3, and BLG genotypes were assigned flat priors in this model, whereas the effects of herd and animal additive genetic were assigned Gaussian prior distributions, and inverse Wishart distributions were assumed for the respective co-variance matrices. Marginal posterior distributions of the parameters of interest were compared before and after the inclusion of the effects of the 3 major genes in the model. The results showed that a high portion of the genetic variance was controlled by the major genes. This was particularly apparent in the qualitative protein profile, which was found to have a higher heritability than the protein fraction contents in milk and their daily yields. When the genes were included individually in the model, CSN2 was the major gene controlling all the casein fractions except for κ-CN, which was controlled directly by the CSN3 gene. The BLG gene had the most influence on the 2 whey proteins. The genetic correlations showed the major genes had only a small effect on the relationships between the protein fractions, but through comparison of the correlation coefficients of the proteins expressed in different ways they revealed potential mechanisms of regulation and competitive synthesis in the mammary gland. The estimates for the effects of the CSN2 and CSN3 genes on protein profiles showed overexpression of protein synthesis in the presence of the B allele in the genotype. Conversely, the β-LG B variant was associated with a lower concentration of β-LG compared with the β-LG A variant, independently of how the protein fractions were expressed, and it was followed by downregulation (or upregulation in the case of the β-LG B) of all other protein fractions. These results should be borne in mind when seeking to design more efficient selection programs aimed at improving milk quality for the efficiency of dairy industry and the effect of dairy products on human health.     

Exploring genome-wide differentiation and signatures of selection in Italian and North American Holstein populations

Among Italian dairy cattle, the Holstein is the most reared breed for the production of Parmigiano Reggiano protected designation of origin cheese, which represents one of the most renowned products in the entire Italian dairy industry. In this work, we used a medium-density genome-wide data set consisting of 79,464 imputed SNPs to study the genetic structure of Italian Holstein breed, including the population reared in the area of Parmigiano Reggiano cheese production, and assessing its distinctiveness from the North American population. Multidimensional scaling and ADMIXTURE approaches were used to explore the genetic structure among populations. We also investigated putative genomic regions under selection among these 3 populations by combining 4 different statistical methods based either on allele frequencies (single marker and window-based) or extended haplotype homozygosity (EHH; standardized log-ratio of integrated EHH and cross-population EHH). The genetic structure results allowed us to clearly distinguish the 3 Holstein populations; however, the most remarkable difference was observed between Italian and North American stock. Selection signature analyses identified several significant SNPs falling within or closer to genes with known roles in several traits such as milk quality, resistance to disease, and fertility. In particular, a total of 22 genes related to milk production have been identified using the 2 allele frequency approaches. Among these, a convergent signal has been found in the VPS8 gene which resulted to be involved in milk traits, whereas other genes (CYP7B1, KSR2, C4A, LIPE, DCDC1, GPR20, and ST3GAL1) resulted to be associated with quantitative trait loci related to milk yield and composition in terms of fat and protein percentage. In contrast, a total of 7 genomic regions were identified combining the results of standardized log-ratio of integrated EHH and cross-population EHH. In these regions candidate genes for milk traits were also identified. Moreover, this was also confirmed by the enrichment analyses in which we found that the majority of the significantly enriched quantitative trait loci were linked to milk traits, whereas the gene ontology and pathway enrichment analysis pointed to molecular functions and biological processes involved in AA transmembrane transport and methane metabolism pathway. This study provides information on the genetic structure of the examined populations, showing that they are distinguishable from each other. Furthermore, the selection signature analyses can be considered as a starting point for future studies in the identification of causal mutations and consequent implementation of more practical application.     

Strain dependent expression of stress response and virulence genes of Listeria monocytogenes in meat juices as determined by microarray

A subgenomic array, encompassing 54 probes targeting genes responsible for virulence, adhesion and stress response in Listeria monocytogenes, was used in order to study their expression in food systems. RNA extracted from L. monocytogenes inoculated in BHI and in situ (i.e. in minced meat and fermented sausage juices) and incubated at 4 °C, was hybridized on the array and the results obtained were compared in order to understand the effect that the food juice has on the expression. Three different strains of L. monocytogenes were tested, in order to determine the effect of the strain provenience. As determined by cluster analysis, each strain behaved in a different way when inoculated in food juices. The goal was to respond to acidic and osmotic stresses encountered in the food, particularly in the fermented sausage juice. No differences in the expression profile between the three strains were observed, when they were inoculated in BHI. On the other hand, in the meat and sausage juices, the iap, gadC and gadE genes, together with different internalin encoding genes, were significantly differentially expressed in the three strains.     

Short communication: Comparing real-time PCR and bacteriological cultures for Streptococcus agalactiae and Staphylococcus aureus in bulk-tank milk samples

For more than 30 yr, a control plan for Streptococcus agalactiae and Staphylococcus aureus has been carried out in more than 1,500 dairy herds of the province of Brescia (northern Italy). From 2010 to 2011, the apparent prevalence of Strep. agalactiae has been relatively stable around 10%, but the apparent prevalence of Staph. aureus has been greater than 40% with an increasing trend. The aim of this paper was to estimate and compare the diagnostic accuracy of 3 assays for the detection of Strep. agalactiae and Staph. aureus in bulk-tank milk samples (BTMS) in field conditions. The assays were a qualitative and a quantitative bacteriological culture (BC) for each pathogen and a homemade multiplex real-time PCR (rt-PCR). Because a gold standard was not available, the sensitivities (Se) and specificities (Sp) were evaluated using a Bayesian latent class approach. In 2012 we collected one BTMS from 165 dairy herds that were found positive for Strep. agalactiae in the previous 2-yr campaigns of eradication plan. In most cases, BTMS collected in these herds were positive for Staph. aureus as well, confirming the wide spread of this pathogen. At the same time we also collected composite milk samples from all the 8,624 lactating cows to evaluate the within-herd prevalence of Strep. agalactiae. Streptococcus agalactiae samples were cultured using a selective medium Tallium Kristalviolette Tossin, whereas for Staph. aureus, we used Baird Parker modified medium with added Rabbit Plasma Fibrinogen ISO-Formulation. In parallel, BTMS were tested using the rt-PCR. Regarding Strep. agalactiae, the posterior median of Se and Sp of the 2 BC was similar [qualitative BC: Se = 98%, posterior credible interval (95%PCI): 94–100%, and Sp = 99%, 95%PCI: 96–100%; quantitative BC: Se = 99%, 95%PCI: 96–100%, and Sp = 99%, 95%PCI: 95–100%] and higher than those of the rt-PCR (at 40 cycle threshold, Se = 92%, 95%PCI: 85–97%; Sp = 94%, 95%PCI: 88–98%). Also in case of Staph. aureus, the posterior medians of BC were generally higher than those of rt-PCR. In fact, although the Se of BC was slightly lower (rt-PCR at 40 cycle threshold, median Se = 99%, 95%PCI: 97–100%, and qualitative BC, median Se = 94%, 95%PCI: 87–99%), the Sp was much higher (rt-PCR at 40 cycle threshold, median Sp = 67%, 95%PCI: 38–97%; qualitative BC, median Sp = 95%; 95%PCI: 76–100%). Our study confirms that BC and rt-PCR are reliable diagnostic tools to detect Strep. agalactiae and Staph. aureus, and rt-PCR results should be confirmed by BC carried out on BTMS and possibly on composite milk samples.     

Molecular epidemiology,antimicrobial activity,and virulence gene clustering of Streptococcus agalactiae isolated from dairy cattle with mastitis in China

Streptococcus agalactiae is a contagious pathogen that causes bovine mastitis worldwide, resulting in considerable economic losses. In this study, we isolated 42 S. agalactiae strains in 379 milk samples from cows with subclinical mastitis on 15 dairy farms in 12 Chinese provinces. Analysis based on capsular typing and multilocus sequence typing, combined with patterns of virulence gene scanning and antimicrobial resistance, identified the lineages and populations of the isolates. We grouped the 42 isolates into 7 sequence types belonging to 6 clonal complexes, mainly CC103 (31/42 isolates; 73.8%). We identified an ST-23 strain named Sa 129 for the first time on Chinese dairy farms—this strain is usually associated with human isolates. Capsular types Ia and II were predominant in capsular typing. The prevalence of virulence profile 1 (bibA, cfb, cspA, cylE, fbsA, fbsB, hylB, and pavA) was 64.3%, and represented the main trend in China. With respect to antimicrobial resistance, most isolates were susceptible to β-lactams, rifamycin, glycopeptides, and oxazolidone; resistance to several antimicrobial agents, including lincomycin, clindamycin, and doxycycline, varied in 4 different regions. Our research provides a profile for the molecular epidemiology, multilocus sequence typing, antimicrobial resistance, and virulence gene clustering of S. agalactiae, and may be beneficial for the clinical monitoring, prevention, and control of mastitis in dairy cattle.     

Single gene and gene interaction effects on fertilization and embryonic survival rates in cattle

Decrease in fertility and conception rates is a major cause of economic loss and cow culling in dairy herds. Conception rate is the product of fertilization rate and embryonic survival rate. Identification of genetic factors that cause the death of embryos is the first step in eliminating this problem from the population and thereby increasing reproductive efficiency. A candidate pathway approach was used to identify candidate genes affecting fertilization and embryo survival rates using an in vitro fertilization experimental system. A total of 7,413 in vitro fertilizations were performed using oocytes from 504 ovaries and semen samples from 10 different bulls. Fertilization rate was calculated as the number of cleaved embryos 48 h postfertilization out of the total number of oocytes exposed to sperm. Survival rate of embryos was calculated as the number of blastocysts on d 7 of development out of the number of total embryos cultured. All ovaries were genotyped for 8 genes in the POU1F1 signaling pathway. Single-gene analysis revealed significant associations of GHR, PRLR, STAT5A, and UTMP with survival rate and of POU1F1, GHR, STAT5A, and OPN with fertilization rate. To further characterize the contribution of the entire integrated POU1F1 pathway to fertilization and early embryonic survival, a model selection procedure was applied. Comparisons among the different models showed that interactions between adjacent genes in the pathway revealed a significant contribution to the variation in fertility traits compared with other models that analyzed only bull information or only genes without interactions. Moreover, some genes that were not significant in the single-gene analysis showed significant effects in the interaction analysis. Thus, we propose that single genes as well as an entire pathway can be used in selection programs to improve reproduction performance in dairy cattle.     

Short communication: Global transcriptome analysis of Lactococcus lactis ssp. lactis in response to gradient freezing

Lactic acid bacteria are often preserved as starter cultures by freezing to extend shelf stability as well as maintain cell viability and acidification activity. Previous studies showed that the endocyte extracted from gradient-freezing pretreated cells could act as lyoprotectant in the lyophilization process of Lactococcus lactis ssp. lactis. In this study, the molecular mechanisms of L. lactis in response to gradient freezing exposure are described using high-throughput sequencing. Nineteen of 56 genes were upregulated after gradient freezing, whereas 37 genes were downregulated. Further validation results of quantitative real-time PCR experiments were consistent with the RNA sequencing. Gene Ontology (http://www.geneontology.org/) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG; https://www.genome.jp/kegg/) pathway were used to analyze the differentially expressed genes. Several pathways, such as glutathione metabolism, ATP-binding cassette transport, metabolism of cell wall and cell membrane components, and stress response-related pathways, were affected by gradient freezing. Six genes relevant to freezing stress response were selected for quantitative real-time PCR, including 3 upregulated genes (hisK, eutD, dukA) and 3 downregulated genes (als, yedF, pepN). The Gene Ontology enrichment and KEGG pathway analyses showed these genes may influence stress response-related pathways, improving the survival of the L. lactis under freezing stress. The identification of these genes deepened an understanding about their response under freezing stress, helping us find potential genes or pathways related to gradient freezing for further research on lyoprotectants.     

Bioactivities and sensory evaluation of Pu-erh teas made from three tea leaves in an improved pile fermentation process

Two commercial Pu-erh teas, 15-year-old Ta-Huang-In and 25-year-old Ta-Hon-In, were used for screening some species of fungi, yeasts, and bacteria, and six of them were isolated and identified as Actinoplanes aurantiacus, Actinoplanes pallidoaurantiacus, Actinoplanes purpeobrunneus, Streptomyces bacillaris, Streptomyces cavourensis subsp. cavourensis, and Streptomyces cinereus. They were selected for inoculated into thee tea leaves (Yun Nan from China, TTES-12 and C. S. Oolong from Taiwan) and fermented for 180 days. The total polyphenols and GABA content, and DPPH radical scavenging effects were determined to examine the tea infusion quality. The samples inoculated with S. cinereus had the highest total polyphenols content and maximum capacity to scavenge DPPH radicals; the highest GABA content was obtained while the sample inoculated with S. bacillaris. Further comparison of these samples with two commercial Pu-erh teas of different ages (Ta-Huang-In, 15-year storage and Ta-Hon-In, 25-year storage) showed that DPPH radical scavenging activity and GABA content of S. bacillaris and S. cinereus fermented tea leaf were higher than these two commercial teas. Sensory evaluation also demonstrated that the taste, flavor, and overall acceptance did not had significant differences between the experimental tea leaves and commercial samples. The present studies revealed that the fresh tea leaves inoculated with the suitable microbes in short period of time will possess a high-quality tea infusion as long-term storage Pu-erh tea.     

Detection of quantitative trait loci and putative causal variants affecting somatic cell score in dairy sheep by using a 50K SNP chip and whole-genome sequencing

This study presents a scan of the ovine genome to identify quantitative trait loci (QTL) influencing the somatic cell score (SCS), a classical indicator of subclinical mastitis in sheep, and a subsequent high-resolution analysis of one of the identified QTL regions based on the analysis of whole-genome sequence data sets. A half-sib commercial population of Churra sheep genotyped with a 50K SNP chip was analyzed using linkage analysis (LA) and combined linkage and linkage disequilibrium analysis (LDLA). By LA, 2 5% chromosome-wide significant QTL on OAR5 and OAR25 and one 5% genome-wide significant QTL on ovine chromosome 20 (OAR20) were detected, whereas 22 significant associations were identified by LDLA. Two of the associations detected by LDLA replicated LA-detected effects (OAR20, OAR25). We compared the detected associations with previously reported QTL in sheep and cattle, and functional candidate genes were identified within the estimated confidence intervals. We then performed a high-resolution analysis of the OAR20 QTL region, the most significant QTL region identified by LA that replicated a QTL previously described in Churra sheep for SCS using microsatellite markers. For that, 2 segregating trios of 2 segregating families for the OAR20 QTL (each including the Qq sire and 2 daughters, QQ and qq) were selected for whole-genome sequencing. The bioinformatic analysis of the 6 sequenced samples performed across the genomic interval considered (14.2–41.7 Mb) identified a total of 227,030 variants commonly identified by 2 independent software packages. For the 3 different concordance tests considered, due to discrepancies regarding the QTL peak in the segregating families, the list of mutations concordant with the QTL segregating pattern was processed to identify the variants identified in immune-related genes that show a moderate/high impact on the encoded protein function. Among a list of 85 missense variants concordant with the QTL segregation pattern that were within candidate immune-related genes, 13 variants distributed across 7 genes [PKHD1, NOTCH4, AGER, ENSOARG00000009395 (HLA-C, Homo sapiens), ENSOARG00000015002 (HLA-B, H. sapiens), MOG, and ENSOARG00000018075 (BoLA, Bos taurus, orthologous to human HLA-A] were predicted to cause deleterious effects on protein function. Future studies should assess the possible associations of the candidate variants identified herein in commercial populations with indicator traits of udder inflammation (SCS, clinical mastitis).     

Genome-wide association study of milk fatty acid composition in Italian Simmental and Italian Holstein cows using single nucleotide polymorphism arrays

Bovine milk is important for human nutrition, but its fat content is often criticized as a risk factor in cardiovascular disease. Selective breeding programs could be used to alter the fatty acid (FA) composition of bovine milk to improve the healthiness of dairy products for human consumption. Here, we performed a genome-wide association study (GWAS) on bovine milk to identify genomic regions or specific genes associated with FA profile and to investigate genetic differences between the Italian Simmental (IS) and Italian Holstein (IH) breeds. To achieve this, we first characterized milk samples from 416 IS cows and 436 IH cows for their fat profile by gas chromatography. Subjects were genotyped with single nucleotide polymorphism array and a single-marker regression model for GWAS was performed. Our findings confirm previously reported quantitative trait loci strongly associated with bovine milk fat composition. More specifically, our GWAS results revealed significant signals on chromosomes Bos taurus autosome 19 and 26 for milk FA. Further analysis using a gene-centric approach and pathway meta-analysis identified not only some well-known genes underlying quantitative trait loci for milk FA components, such as FASN, SCD, and DGAT1, but also other significant candidate genes, including some with functional roles in pathways related to “Lipid metabolism.” Highlighted genes related to FA profile include ECI2, PCYT2, DCXR, G6PC3, PYCR1, and ALG12 in IS, and CYP17A1, ACO2, PI4K2A, GOT1, GPT, NT5C2, PDE6G, POLR3H, and COX15 in IH. Overall, the breed-specific association outcomes reflect differences in the genetic backgrounds of the IS and IH breeds and their selective breeding histories.     

Identification of functional candidate variants and genes for feed efficiency in Holstein and Jersey cattle breeds using RNA-sequencing

The identification of functional genetic variants and associated candidate genes linked to feed efficiency may help improve selection for feed efficiency in dairy cattle, providing economic and environmental benefits for the dairy industry. This study used RNA-sequencing data obtained from liver tissue from 9 Holstein cows [n = 5 low residual feed intake (RFI), n = 4 high RFI] and 10 Jersey cows (n = 5 low RFI, n = 5 high RFI), which were selected from a single population of 200 animals. Using RNA-sequencing, 3 analyses were performed to identify: (1) variants within low or high RFI Holstein cattle; (2) variants within low or high RFI Jersey cattle; and (3) variants within low or high RFI groups, which are common across both Holstein and Jersey cattle breeds. From each analysis, all variants were filtered for moderate, modifier, or high functional effect, and co-localized quantitative trait loci (QTL) classes, enriched biological processes, and co-localized genes related to these variants, were identified. The overlapping of the resulting genes co-localized with functional SNP from each analysis in both breeds for low or high RFI groups were compared. For the first two analyses, the total number of candidate genes associated with moderate, modifier, or high functional effect variants fixed within low or high RFI groups were 2,810 and 3,390 for Holstein and Jersey breeds, respectively. The major QTL classes co-localized with these variants included milk and reproduction QTL for the Holstein breed, and milk, production, and reproduction QTL for the Jersey breed. For the third analysis, the common variants across both Holstein and Jersey breeds, uniquely fixed within low or high RFI groups were identified, revealing a total of 86,209 and 111,126 functional variants in low and high RFI groups, respectively. Across all 3 analyses for low and high RFI cattle, 12 and 31 co-localized genes were overlapping, respectively. Among the overlapping genes across breeds, 9 were commonly detected in both the low and high RFI groups (INSRR, CSK, DYNC1H1, GAB1, KAT2B, RXRA, SHC1, TRRAP, PIK3CB), which are known to play a key role in the regulation of biological processes that have high metabolic demand and are related to cell growth and regeneration, metabolism, and immune function. The genes identified and their associated functional variants may serve as candidate genetic markers and can be implemented into breeding programs to help improve the selection for feed efficiency in dairy cattle.     

Genome-wide association study for milking speed in French Holstein cows

Using a combination of data from the BovineSNP50 BeadChip SNP array (Illumina, San Diego, CA) and a EuroGenomics (Amsterdam, the Netherlands) custom single nucleotide polymorphism (SNP) chip with SNP pre-selected from whole genome sequence data, we carried out an association study of milking speed in 32,491 French Holstein dairy cows. Milking speed was measured by a score given by the farmer. Phenotypes were yield deviations as obtained from the French evaluation system. They were analyzed with a linear mixed model for association studies. We identified SNP on 22 chromosomes significantly associated with milking speed. As clinical mastitis and somatic cell score have an unfavorable genetic correlation with milking speed, we tested whether the most significant SNP on these 22 chromosomes associated with milking speed were also associated with clinical mastitis or somatic cell score. Nine hundred seventy-one genome-wide significant SNP were associated with milking speed. Of these, 86 were associated with clinical mastitis and 198 with somatic cell score. The most significant association signals for milking speed were observed on chromosomes 7, 8, 10, 14, and 18. The most significant signal was located on chromosome 14 (ZFAT gene). Eleven novel milking speed quantitative trait loci (QTL) were observed on chromosomes 7, 10, 11, 14, 18, 25, and 26. Twelve candidate SNP for milking speed mapped directly within genes. Of these, 10 were QTL lead SNP, which mapped within the genes HMHA1, POLR2E, GNB5, KLHL29, ZFAT, KCNB2, CEACAM18, CCL24, and LHPP. Limited pleiotropy was observed between milking speed QTL and clinical mastitis.