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.

Successful editing and maintenance of lactogenic gene expression in primary bovine mammary epithelial cells.

In vitro investigation of bovine lactation processes is limited by a lack of physiologically representative cell models. This deficiency is most evident through the minimal or absent expression of lactation-specific genes in cultured bovine mammary tissues. Primary bovine mammary epithelial cells (pbMECs) extracted from lactating mammary tissue and grown in culture initially express milk protein transcripts at relatively representative levels. However, expression drops dramatically after only three or four passages, which greatly reduces the utility of primary cells to model and further examine lactogenesis. To investigate the effects of alternate alleles in pbMECs including effects on transcription, we have developed methods to deliver CRISPR-Cas9 gene editing reagents to primary mammary cells, resulting in very high editing efficiencies. We have also found that culturing the cells on an imitation basement membrane composed of Matrigel, results in the restoration of a more representative lactogenic gene expression profile and the cells forming three-dimensional structures in vitro. Here, we present data from four pbMEC lines recovered from pregnant cows and detail the expression profile of five key milk synthesis genes in these MECs grown on Matrigel. Additionally, we describe an optimised method for preferentially selecting CRISPR-Cas9-edited cells conferring a knock-out of DGAT1, using fluorescence-activated cell sorting (FACS). The combination of these techniques facilitates the use of pbMECs as a model to investigate the effects of gene introgressions and genetic variation in lactating mammary tissue.