A common regulatory haplotype doubles lactoferrin concentration in milk

BACKGROUND

Bovine lactoferrin (Lf) is an iron absorbing whey protein with antibacterial, antiviral, and antifungal activity. Lactoferrin is economically valuable and has an extremely variable concentration in milk, partly driven by environmental influences such as milking frequency, involution, or mastitis. A significant genetic influence has also been previously observed to regulate lactoferrin content in milk. Here, we conducted genetic mapping of lactoferrin protein concentration in conjunction with RNA-seq, ChIP-seq, and ATAC-seq data to pinpoint candidate causative variants that regulate lactoferrin concentrations in milk.

RESULTS

We identified a highly-significant lactoferrin protein quantitative trait locus (pQTL), as well as a cis lactotransferrin (LTF) expression QTL (cis-eQTL) mapping to the LTF locus. Using ChIP-seq and ATAC-seq datasets representing lactating mammary tissue samples, we also report a number of regions where the openness of chromatin is under genetic influence. Several of these also show highly significant QTL with genetic signatures similar to those highlighted through pQTL and eQTL analysis. By performing correlation analysis between these QTL, we revealed an ATAC-seq peak in the putative promotor region of LTF, that highlights a set of 115 high-frequency variants that are potentially responsible for these effects. One of the 115 variants (rs110000337), which maps within the ATAC-seq peak, was predicted to alter binding sites of transcription factors known to be involved in lactation-related pathways.

CONCLUSIONS

Here, we report a regulatory haplotype of 115 variants with conspicuously large impacts on milk lactoferrin concentration. These findings could enable the selection of animals for high-producing specialist herds.

RNA-Seq Analysis of Peripheral Whole Blood from Dairy Bulls with High and Low Antibody-Mediated Immune Responses-A Preliminary Study

Enhancing the immune response through breeding is regarded as an effective strategy for improving animal health, as dairy cattle identified as high immune responders are reported to have a decreased prevalence of economically significant diseases. The identification of differentially expressed genes (DEGs) associated with immune responses might be an effective tool for breeding healthy dairy cattle. In this study, antibody-mediated immune responses (AMIRs) were induced by the immunization of hen egg white lysozyme (HEWL) in six Chinese Holstein dairy bulls divided into high- and low-AMIR groups based on their HEWL antibody level. Then, RNA-seq was applied to explore the transcriptome of peripheral whole blood between the two comparison groups. As a result, several major upregulated and downregulated genes were identified and attributed to the regulation of locomotion, tissue development, immune response, and detoxification. In addition, the result of the KEGG pathway analysis revealed that most DEGs were enriched in pathways related to disease, inflammation, and immune response, including antigen processing and presentation, Staphylococcus aureus infection, intestinal immune network for IgA production, cytokine-cytokine receptor interaction, and complement and coagulation cascades. Moreover, six genes (BOLA-DQA5, C5, CXCL2, HBA, LTF, and COL1A1) were validated using RT-qPCR, which may provide information for genomic selection in breeding programs. These results broaden the knowledge of the immune response mechanism in dairy bulls, which has strong implications for breeding cattle with an enhanced AMIR.

Kompetitive Allele Specific PCR Genotyping of 89 SNPs in Romanian Spotted and Romanian Brown Cattle Breeds and Their Association with Clinical Mastitis

Mastitis is the most common production disease in the dairy sector worldwide, its incidence being associated with both cows' exposure to bacteria and the cows' genetic make-up for resistance to pathogens. The objective of our study was to analyse 89 missense SNPs belonging to six genes (CXCR2, CXCL8, TLR4, BRCA1, LTF, BOLA-DRB3), which were found to be associated with genetic resistance or susceptibility to mastitis. A total of 298 cattle (250 Romanian Spotted and 48 Romanian Brown) were genotyped by Kompetitive Allele Specific PCR (KASP) and a chi-squared test was used for genetic association studies with clinical mastitis. A total of 35 SNPs (39.3%) among the selected 89 SNPs were successfully genotyped, of which 31 markers were monomorphic. The polymorphic markers were found in two genes: TLR4 (rs460053411) and BOLA-DRB3 (rs42309897, rs208816121, rs110124025). The polymorphic SNPs with MAF > 5% and call rates > 95% were used for the association study. The results showed that rs110124025 in the BOLA-DRB3 gene was significantly associated with mastitis prevalence (p ≤ 0.05) in both investigated breeds. Current results show that the SNP rs110124025 in the BOLA-DRB3 gene can be used as a candidate genetic marker in selection for mastitis resistance in Romanian dairy cattle.

Gene expression adjustment of inflammatory mechanisms in dairy cow mammary gland parenchyma during host defense against staphylococci

The aim of the study was to identify differences in the expression of splice variants of the PRMT2, LTF and C4A genes in the mammary glands of healthy dairy cows and those infected with staphylococci. An expression study was conducted on 38 Polish Holstein-Friesian dairy cows who were removed from the herd owing to subclinical and chronic mastitic or reproductive issues. Two days before slaughter, milk samples were taken for microbiological analysis and examined for the presence of bacteria. The mammary gland parenchyma samples with a predominance of secretory tissue were taken; these were divided into three groups according to the health status of the mammary gland: H (without pathogenic bacteria in milk), CoNS (with coagulase-negative staphylococci in milk), and CoPS (with coagulase-positive staphylococci in milk). Two of the investigated genes, LTF and C4A, demonstrated variants unequivocally expressed in infected tissue. Two LTF gene variants were found to be associated with cow health status, and with the type of bacteria causing mastitis (CoPS or CoNS). In addition, the expression of C4A isoforms differed with regard to mastitis etiology groups. The comprehensive evaluation of PRMT2 transcript suggested that the gene may also be involved in course of mastitis: two of four PRMT2 transcripts showed increased expression in the mammary gland of the CoPS group compared to controls. The obtained results are important for the knowledge on the etiology of bovine mastitis. The effects of the identified mastitis-relevant splice variants need to be further explored on the protein level to verify the suitability of splice variants and recognize their contribution towards the disease phenotypes and course.

Acute phase protein expressions in secretory and cistern lining epithelium tissues of the dairy cattle mammary gland during chronic mastitis caused by staphylococci

Background

Mastitis is the most common disease in dairy cattle and the costliest for the dairy farming industry, as it lowers milk yield and quality. Mastitis occurs as a result of interactions between microorganisms and the individual genetic predispositions of each animal. Thus, it is important to fully understand the mechanisms underlying these interactions. Elucidating the immune response mechanisms can determine which genetic background makes an animal highly resistant to mastitis. We analyzed the innate immune responses of dairy cows naturally infected with coagulase-positive staphylococci (CoPS; N = 8) or coagulase-negative staphylococci (CoNS; N = 7), causing persistent mastitis (after several failed treatments) vs. infection-free (i.e., healthy [H]; N = 8) dairy cows. The expressions of the acute phase protein genes serum amyloid A3 (SAA3), haptoglobin (HP), ceruloplasmin (CP) genes in the tissues most exposed to pathogens— mammary gland cistern lining epithelial cells (CLECs) and mammary epithelial cells (MECs)—were analyzed.

Results

We found constitutive and extrahepatic expressions of the studied genes in both tissue types. HP expression in the MECs of the CoPS-infected group was higher than in the H group (p ≤ 0.05). Moreover, higher SAA3 expression in the CoPS and CoNS groups than in the H group (p = 0.06 and 0.08, respectively) was found. No differences between SAA3 and HP in CLECs were revealed, regardless of the pathogen type. However, higher expression of CP (p ≤ 0.05) in the CoPS group than in the H group was noted.

Conclusions

The expressions of selected acute phase proteins were similar between CLECs and MECs, which means that CLECs are not only a mechanical barrier but are also responsible for the biological immune response. Our findings agree with the results of other authors describing the immunological response of MECs during chronic mastitis, but the results for CLECs are novel.

Lactoferrin concentration and expression in New Zealand cows milked once or twice a day

This study evaluated the concentration and expression of lactoferrin (LF) in cows selected for once a day (OAD) milking compared to twice a day (TAD) milking. Milk samples were collected from the Massey University TAD and OAD herds. Milk traits and expression of LF and insulin-like growth factor 1 (IGF-1) were analyzed with a general linear model that included the fixed effects of milking frequency, lactation number, interaction between milking frequency and lactation number, and as covariates proportion of F, heterosis F × J and deviation from the herd median calving date. Cows milked OAD produced milk with higher (p < .01) concentrations of protein and lactose than TAD milked cows. Compared to TAD cows, cows milked OAD had higher expression of the LF gene (1.40 vs. 1.29 folds, p = .03) and the IGF-1 gene (1.69 vs. 1.48 folds, p = .007). The correlation between the expression of LF gene and the concentration of LF in milk was strong (r = .66 p < .001), but the correlation between the expression of the IGF-1 gene and LF concentration was stronger (r = .94, p < .001). These results suggest that milking frequency affects the milk composition and expression of milk composition genes at early lactation.

Structural and functional analysis of the signaling lymphocytic activation molecule family 7 (SLAMF7) gene in response to infection with coagulase-negative and coagulase-positive staphylococci

Splice variants of the signaling lymphocytic activation molecule family 7 (SLAMF7) gene have been identified, and differences in the expression of this gene have been demonstrated at the mRNA level in the mammary glands of healthy and mastitis-infected dairy cows. At the same time, significant associations have been found between a deletion in the SLAM7 gene exon, the occurrence of different splice variants, and the occurrence of mastitis in one group of dairy cows. An expression study was conducted on 40 Polish Holstein-Friesian dairy cows of the Black and White variety (group I). Milk samples were taken for microbiological analysis 2 d before slaughter and examined for the presence of bacteria. Immediately after slaughter, mammary tissue samples were taken and divided into 3 groups according to the health status of the mammary gland: healthy (without pathogenic bacteria in milk), coagulase-negative staphylococci (CNS), and coagulase-positive staphylococci (CPS). Based on different SLAMF7 gene DNA fragments, 2 alternative variants of this gene (V1 and V2) and complete gene expression were identified. Separate analyses performed for each isoform showed that the health status of the cow was strongly associated with the expression level of individual variants. The highest expression was detected for the SLAMF7 complete amplicon in healthy cows, and in the CNS and CPS cows the expression of this variant was also higher than V1 and V2. Sanger sequencing was applied to detect the polymorphism/indel variant in the second exon of the SLAMF7 gene probably having the greatest effect on the protein structure and function of SLAMF7. Two genotypes were detected: AA (wild-type) and AB (insertion A). In healthy cows, the frequency of homozygotes AA was higher than the heterozygotes, whereas in the infected animals, the genotypic distribution was the opposite. An association analysis between the identified polymorphism and production traits-including somatic cell count, as well as lactose, protein, and casein content and yield as indicators of subclinical mastitis occurrence-was performed on the group II cows (166 Polish Holstein-Friesian dairy cows). Unfortunately, due to the low number of AB animals, no relationship was demonstrated between genotype in the second exon and the health status of cows. Additionally, the difference in the percentage of SLAMF7-targeted DNA methylation between the groups of animals was not significant, with an average of ∼66 to 68%.