- Open Access
Knockout of butyrophilin subfamily 1 member A1 (BTN1A1) alters lipid droplet formation and phospholipid composition in bovine mammary epithelial cells
Journal of Animal Science and Biotechnology volume 11, Article number: 72 (2020)
Milk lipids originate from cytoplasmic lipid droplets (LD) that are synthesized and secreted from mammary epithelial cells by a unique membrane-envelopment process. Butyrophilin 1A1 (BTN1A1) is one of the membrane proteins that surrounds LD, but its role in bovine mammary lipid droplet synthesis and secretion is not well known.
The objective was to knockout BTN1A1 in bovine mammary epithelial cells (BMEC) via the CRISPR/Cas9 system and evaluate LD formation, abundance of lipogenic enzymes, and content of cell membrane phospholipid (PL) species. Average LD diameter was determined via Oil Red O staining, and profiling of cell membrane phospholipid species via liquid chromatography-tandem mass spectrometry (LC-MS/MS).
Lentivirus-mediated infection of the Cas9/sgRNA expression vector into BMEC resulted in production of a homozygous clone BTN1A1(−/−). The LD size and content decreased following BTN1A1 gene knockout. The mRNA abundance of fatty acid synthase (FASN) and peroxisome proliferator-activated receptor-gamma (PPARG) was downregulated in the BTN1A1(−/−) clone. Subcellular analyses indicated that BTN1A1 and LD were co-localized in the cytoplasm. BTN1A1 gene knockout increased the percentage of phosphatidylethanolamine (PE) and decreased phosphatidylcholine (PC), which resulted in a lower PC/PE ratio.
Results suggest that BTN1A1 plays an important role in regulating LD synthesis via a mechanism involving membrane phospholipid composition.
The process of milk fat secretion begins in the endoplasmic reticulum of mammary epithelial cells (MEC) from which lipid droplets (LD) begin to form . During secretion from cells into milk, LD are progressively enveloped by the apical plasma membrane forming milk fat globules (MFG) with sizes ranging from 0.2–15 μm [2,3,4]. Inside MFG, there is a lipid core surrounded by a complex phospholipid- and protein-coated membrane (MFGM) . The size and features of MFG have received particular interest due to their influence on the manufacturing properties and nutritional quality of dairy products [6,7,8]. It has been suggested that intracellular LD formed by mammary epithelial cells are related to the size of secreted MFG [9, 10] and, thus, the mechanism regulating LD size in the cytoplasm is relevant to the properties of MFG.
The abundance of some proteins in the MFGM of MEC affects formation of LD [11,12,13,14]. Proteomic analysis of the composition of the MFGM revealed that the most-abundant proteins include butyrophilin (BTN1A1), perilipin 2 (PLIN2; formerly adipophilin) and xanthine dehydrogenase (XDH) . The BTN1A1 protein accounts for more than 20% of total protein in the MFGM [3, 16]. Abundance of these proteins increases during lactation in the bovine mammary gland . McManaman et al.  proposed that BTN may function as an integral receptor for cytoplasmic LD, and budding of the droplets at the cell surface initiates formation of complexes between BTN and XDH. Robenek et al.  suggested that a crucial feature of MFG secretion is the establishment of a network of adhesion sites containing BTN in the secretory granule monolayer. Clearly, BTN1A1 proteins are involved in MFG formation, but little is known about the mechanism whereby BTN1A1 aids in LD formation in bovine MEC.
The clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated 9 (Cas9) system can be used for genome editing in mammalian cells via single-guide RNA (sgRNA) [19,20,21]. Using sgRNA leads to Cas9-mediated cleavage of the target DNA and the protospacer adjacent motif, resulting in double-stranded breaks . In mammalian cells, double-stranded breaks (DSB) are mostly repaired by the non-homologous end-joining pathway that enables efficient construction of knockout alleles by introducing small insertions or deletions that result in loss of target protein expression . The aim of this study was to investigate the effect of knockout of BTN1A1 on LD formation using the CRISPR/Cas9 system. The LD size was measured by Oil Red O staining, and profiling of cell membrane phospholipid species was analysed via liquid chromatography-tandem mass spectrometry.
Materials and methods
Bovine mammary epithelial cell isolation and culture
Bovine mammary epithelial cell (BMEC) were isolated and purified as described previously . Briefly, mammary tissue (150 mg) was harvested from 4 to 5 years old late-lactation dairy cows and immediately transported to the laboratory. Samples were washed with DPBS (D8662, Sigma) and cut into three 1 mm pieces. Tissue was transferred with tweezers onto clean cell culture dishes at 37 °C in an atmosphere of 5% CO2 and 95% air. After 4 h, 5 mL of basal medium was added into the culture. The basal medium was replaced with fresh medium every 24 h until cells were thoroughly distributed across the bottom of the dish. Basal medium was composed of DMEM/F12 (12400–024, Gibco, New York, USA) and 10% (v/v) fetal bovine serum (16000–044, Gibco, New York, USA). Subsequently, BMEC were enriched by selective detachment with 0.25% trypsin (Gibco, Grand Island, NY, USA). After 3 min of enzyme digestion, detached fibroblasts were removed by washing with DPBS and epithelial cells attached to the dish were allowed to grow by addition of fresh medium. The BMEC were continuously purified using the same method. The purified BMEC (1 × 106 cells/mL) were suspended in freezing medium for cryopreservation in liquid nitrogen.
After thawing, the BMEC were cultured in basal medium similar to that of Peterson et al.  with modifications. The basal medium was composed of DMEM/F12 with 10% fetal bovine serum, insulin (5 μg/mL, I6634, Sigma, St. Louis, MO, USA), hydrocortisone (1 μg/mL, H0888, Sigma), transferrin (5 μg/mL, T1428, Sigma), progesterone (1 μg/mL, P8783, Sigma), and epidermal growth factor (100 ng/mL, SRP3027, Sigma). At approximately 24 h before applying treatments (approximately 70–80% confluence), cells were cultured in lactogenic medium as reported by Kadegowda et al.  with modifications. The lactogenic medium was DMEM/F12 without serum, containing bovine Serum Albumin (1 mg/mL, A1933, Sigma), insulin (5 μg/mL, I6634, Sigma), hydrocortisone (1 μg/mL, H0888, Sigma) and prolactin (1 μg/mL, 682-PL, R&D Systems) to induce protein synthesis.
We designed primers and measured mRNA abundance of β-casein in BMEC by RT-PCR to confirm the lactogenic state of BMEC. In order to determine the suitable puromycin tolerance concentration in BMEC, cells were seeded in 35-mm cell culture dish and when they approached 70–80% confluence, puromycin (P8833, Sigma) was added to the culture medium at a final concentration of 0, 2.5, 5.0, 7.5 or 10.0 μg/mL. After culturing for 72 h, the minimum lethal dose was selected based on cell survival determined via microscopy.
BTN1A1-GFP plasmid construction and overexpression
To construct the BTN1A1 overexpression plasmid, full-length open reading frame sequences were amplified from cDNA using PCR and subcloned into the restriction sites of the pEGFP-N1 vector. The following PCR primer pair was used to clone the bovine BTN1A1 open reading frame (Gene ID: NM_174508.2): Forward 5′-CTCGAGATGGCAGTCTTTCCAAACT-3′, encoding a XhoI restriction site; Reverse 5′-AAGCTTAGGCACCCCTTGGCTG-3′, encoding a HindIII restriction site. Total RNA was extracted from bovine mammary tissue samples using TRIzol reagent (15,596, Invitrogen) and reverse-transcribed to cDNA as a template for PCR. The PCR products were transformed into TOP10 competent cells. Positive clones were selected, and plasmids extracted. After digestion and identification, the BTN1A1-GFP expression vector was sequenced and verified.
The BMEC were counted and seeded at a density of 5 × 105 in 6-well plates and then cultured with lactogenic medium until 70–80% confluence before transient transfection with BTN1A1-GFP expression vector using Lipofectamine 3000 transfection reagent (L3000008, Invitrogen) according to the manufacturer’s protocol. Briefly, the vector DNA and lipofectamine 3000 reagent mixture was prepared by diluting in OptiMEM medium followed by thorough mixing. The mixture was incubated for 15 min at room temperature before addition of the DNA-lipid complex to the cells. After culture for 48 h, washing with PBS, the LD were stained with Nile Red. Abundance of the BTN1A1-GFP fusion protein and LD were determined under an inverted fluorescence microscope to analyze the co-localization of the BTN1A1 and lipid droplet.
Construction of Cas9/sgRNA expression vectors
Exons 2 and 3 of the bovine BTN1A1 gene (GenBank ID 282157) were selected as target loci for editing, and sgRNA were designed as reported previously  . All potential 20 bp primer sequences followed by 5′-CACC-3′ were scored and analysed based on several factors, including mismatches and the number of off-target sites. Three sgRNA (Table 1) were selected based on their predicted score, synthesised by Invitrogen (Beijing, China), and subcloned into the lentiCRISPR v2 vector (52961, Addgene), resulting in BTN1A1-sgRNA1, BTN1A1-sgRNA2 and BTN1A1-sgRNA3.
The HEK293T/17 cells were seeded in six-well plates at 1 × 106 cells/well. The following day, HEK293T/17 cells were co-transfected with 2 μg of lentiCRISPR v2 -BTN sgRNA 1–3 plasmids, 1.5 μg of psPAX2 plasmid, and 0.5 μg of pMD2.G plasmid. All transfections were carried out with TurboFect reagent (Thermo Scientific) as recommended by the manufacturer. After 8 h (37 °C, 5% CO2), the medium was replaced with growth medium supplemented with 10% FBS, and after a further 24 h the medium was changed. Lentiviruses were harvested at 48–72 h, virus-containing media were pooled, centrifuged at 800×g for 5 min, and supernatants were used to infect BMEC.
Cell infection and T7E1 assays
The BMEC were transfected with lentivirus supernatant until 60–70% confluence. After transfecting for 48 h, the medium was replaced with one containing puromycin (5 μg/mL). After puromycin selection of BMEC for 72 h, genomic DNA was isolated using a MiniBEST Universal Genomic DNA Extraction Kit (Takara, China) according to the manufacturer’s instructions and stored − 20 °C until use. PCR amplification of the target gene using primers BTN1A1-Test1 and BTN1A1-Test2 primer (Table 1) was performed using Q5 High-Fidelity DNA Polymerase (MD491S, New England Biolabs) with 35 cycles at 98 °C for 10 s, 65 °C for 20 s, and 72 °C for 20 s.
For T7E1 assays, 200 ng of PCR product was incubated at 95 °C for 5 min in 1× NEBuffer 2, then slowly cooled from 95 °C to 85 °C at a rate of − 2 °C/s, and to 25 °C at a rate of − 0.1 °C/s. After annealing, 5 U of T7 endonuclease I (M0302, New England Biolabs) was added to each sample and reactions were incubated at 37 °C for 15 min. T7E1-digested products were separated on a 2.5% agarose gel, stained with ethidium bromide, and visualised using a gel documentation system (BIO-RAD).
Generation of individual cell clones and sequencing verification
Individual cell clones were isolated from the positive cell pool by the limiting dilution method as described previously . After transfection with Cas9/sgRNA- BTN1A1–3 expression plasmids as described above, positive cells were diluted to one cell per 100 μL of medium, inoculated into 96-well cell culture plates, and cultured for 10–14 days to obtain single clone colonies. Genomic DNA was isolated from individual clones, and PCR was performed using the BTN1A1-Test primer. PCR products were purified and subjected to sequencing and TA cloning sequencing.
Real-time quantitative PCR
The WT (normal BMEC) and BTN1A1(−/−) BMEC were plated at 50,000 cells per well in 6-well plates for cellular RNA extraction (n = 3). After culturing with lactogenic medium for 48 h, total RNA was extracted from BMEC using a TRIzol Plus RNA Purification Kit (12183555, Thermo Fisher Scientific) according to the manufacturer’s protocol. Briefly, first-strand cDNA was synthesised using a Prime Script RT Reagent kit with gDNA Eraser (Takara, Japan), and qPCR analysis was performed using SYBR fast qPCR mix (Takara) and a 7500 fast Real-Time PCR System (ABI). Primers for amplification of target genes by real-time PCR are listed in Table 1. The overall percentage relative mRNA abundance for each gene among all those measured was calculated using the inverse of PCR efficiency raised to ∆Ct. The efficiency of PCR amplification for each gene was calculated using the standard curve method [E 10(− 1/standard curve slope)]. The ubiquitously expressed transcript (UXT), ribosomal protein S15 (RPS15) and ribosomal protein S9 (RPS9) were selected as internal control genes for normalization .
Lipid droplet staining
The WT and BTN1A1(−/−) BMEC were plated at 30,000 cells per well in 12-well plates for Nile red staining (n = 3). After culturing with lactogenic medium for 48 h, LD staining was performed as described previously . Briefly, BMEC were washed twice with phosphate-buffered saline (PBS) for 3 min, fixed with 4% paraformaldehyde for 20 min, and Nile Red (10 μg/mL) added to each well and incubated for 15 min at room temperature. After washing with PBS three times, plates were placed in a Varioskan LUX full-wavelength microplate reader (Thermo Scientific) and fluorescence intensity was measured at an excitation wavelength of 480 nm and an emission wavelength of 575 nm. All experiments were conducted in triplicate.
Oil red O staining and lipid droplet size determination
The WT and BTN1A1(−/−) BMEC were plated at 60,000 cells per 35-mm Glass Bottom Cell Culture Dish (NEST Biotechnology, China) for oil red staining (n = 3). After culturing with lactogenic medium for 48 h, Oil-Red-O staining was performed using Oil-Red-O working solution (G1262, Solarbio) with modifications. Briefly, cells were washed thrice with PBS and fixed with 4% paraformaldehyde for 20 min at room temperature. After three washes in PBS, cells were washed with isopropanol and stained with Oil-Red-O working solution (G1260, Solarbio) for 25 min, then washed with PBS, and nuclei re-stained with hematoxylin for 1–2 min. Cells were covered with distilled water and images were visualised using an Olympus IX73 fluorescence microscope equipped with an Olympus DP80 digital camera. Forty LD were randomly selected from each cell culture dish and diameter measured using Cell Sens standard software (version 1.13; Olympus). Average LD diameter values represent the mean size of all 120 LD per group. Lipid droplets were divided into four size groups (0 < size < 2.0 μm, 2.0 μm < size < 2.5 μm, 2.5 μm < size < 3.0 μm and size > 3.0 μm). The number of LD were counted within each of these four sizes.
Lipid extraction and liquid chromatography-tandem mass spectrometry (LC-MS/MS)
The WT and BTN1A1(−/−) BMEC were plated at 150,000 cells per 60-mm in a Cell Culture Dish for lipid extraction (n = 6). After culturing with lactogenic medium for 48 h, cells were harvested using trypsin, homogenized in 200 μL of ultra-pure water, mixed with 200 μL of pre-cooled methanol, and 800 μL methyl tertiary-butyl ether. After vortexing and sonication in a low temperature water bath for 20 min, and room temperature for 30 min, samples were centrifuged at 14,000×g for 15 min at 4 °C, and the upper organic phase was carefully removed and dried using nitrogen gas. Samples were dissolved in 200 μL of isopropanol/methanol (1/1, v/v) prior to LC-MS/MS detection.
Separation of lipids was performed on an ACQUITY UPLC CSH C18 column (1.7 μm internal diameter, 2.1 mm × 100 mm; Waters Corp.) using an ultra-high-performance liquid chromatography (UHPLC) Nexera LC-30A instrument (Shimadzu Technology, Japan). Samples separated by UHPLC were analysed by a Q-Exactive Plus Mass Spectrometer (Thermo Fisher Scientific) in positive ion mode with a heater temperature of 300 °C, a sheath gas flow rate of 45 arb, an aux gas flow rate of 15 arb, a sweep gas flow rate of 1 arb, a spray voltage of 3.0 kV, a capillary temp of 350 °C, an S-lens RF level of 50%, and an MS1 scan range of 200–1800. Spectra were also obtained in negative ion mode with a heater temperature of 300 °C, a sheath gas flow rate of 45 arb, an aux gas flow rate of 15 arb, a sweep gas flow rate of 1 arb, a spray voltage of 2.5 kV, a capillary temp of 350 °C, an S-lens RF level of 60%, and an MS1 scan range 250–1800.
Peak recognition, lipid identification, peak extraction, peak alignment, and quantitative analysis were performed with Lipid Search software version 4.1 (Thermo Fisher Scientific). Lipid Search software is composed of 8 categories, 300 subclasses, and a database of approximately 1.7 million lipid species. Through the identification algorithm of sub ion, parent ion and neutral loss scanning, the software performs systematic and reliable qualitative analysis of lipids. For data extracted by Lipid Search, lipid molecules with RSD > 30% were deleted. Total peak area was normalised and used for further quantification and analysis. The PLS-DA model scores R2 = 0.996 and Q2 = 0.975 were obtained by partial least square discriminant analysis of all experimental and QC samples. According to the variable importance in projection (VIP) obtained with the PLS-DA model, influence intensity of each lipid expression pattern was measured, and different lipids with biological significance were obtained. Total phospholipid content was the sum of all peak intensities of phospholipids species identified. Peak intensities of lipid species were summed into their respective classes. The relative percentage of lipid classes was calculated by dividing peak intensities of each classes by total phospholipid peak intensities.
Statistical analysis of all data was performed via SPSS 20.0 (IBM, Armonk, NY, USA). All experiments were conducted in triplicate. Data are reported as Means ± SE. Significant differences were determined using t-tests. All P values < 0.05 were considered statistically significant. The most differentially expressed phospholipid species in the lipidomics data were identified at a cutoff value P < 0.05, VIP > 1, Fold change> 2.5 or < 0.5.
The β-casein gene was highly-expressed in cultured BMEC (Additional file 1: Fig. S1). Most cells survived at the dose of 2.5 μg/mL puromycin. When dose of puromycin was greater than 5.0 μg/mL, most cells died. Thus, the minimum lethal dose of puromycin against BMEC was 5.0 μg/mL (Additional file 1: Fig. S2). This concentration was used for subsequent BMEC selection.
BTN1A1 protein and lipid droplet co-localization
CRISPR/Cas9-mediated knockout of the BTN1A1 gene
We designed three different sgRNA targeting exons 2 and 3 of the BTN1A1 gene (Fig. 2a). T7E1 cleavage assays yielded two DNA bands for sgRNA3 (Fig. 2b), indicating that sgRNA3 successfully edited the genomic DNA of selected cells. Based on the T7E1 cleavage results, BTN1A1-sgRNA3 was used in subsequent experiments and infected into BMEC to generate single clones. We obtained 4 single clones following gene editing. To assess knockout efficiency, the BTN1A1 protein levels in different mutants were measured by western blotting (Additional file 1: Fig. S3). Clone 8 was deemed to have sufficient BTN1A1 knockout, and sequencing revealed this clone included a 10-bp deletion in the target sequence in both alleles (Fig. 2c and d). Clone 8 was, thus, used as the homozygous BTN1A1(−/−) clone in further experiments.
BTN1A1 knockout affects abundance of lipogenic genes
Knockout of BTN1A1downregulated (P < 0.05) mRNA abundance of FASN and PPARG (Fig. 3), but had no significant effect (P > 0.05) on PLIN2 abundance.
BTN1A1 knockout decreases lipid droplet content and diameter
After BTN1A1 knockout, Nile Red staining revealed that compared with the WT group, LD formation was decreased. With increasing culture time, LD content of BTN1A1(−/−) cells decreased from 1065 (12 h) to 644 (48 h), demonstrating a time-dependent effect (Fig. 4a). After knockout of BTN1A1, average LD diameter decreased from 2.31 μm in WT cells to 2.16 μm in BTN1A1(−/−) cells (P < 0.01; Fig. 4b Fig. 4c). Number of LD in a specific diameter range was determined for WT and BTN1A1(−/−) groups (Fig.4d). Compared with WT, the number of small LD (diameter < 2.0 μm) was greater, and large LD (diameter > 3.0 μm) lower in BTN1A1(−/−) cells (P < 0.05).
Knockout of BTN1A1 affects the cell membrane phospholipid composition
A total of 1030 phospholipid species were identified in WT and BTN1A1(−/−) cells by LCMS/MS analysis (Fig. 5a, Additional file 1 Table S1). PC and PE were the most predominant lipids classes (285 and 212, respectively). Peak intensities of lipid species were summed into their respective classes. The percentage of the polar lipid class out of total phospholipids in both WT and BTN1A1(−/−) cell is shown in Table 2. The PC, PE, PI, PS and SM accounted for more than 95% of total phospholipids. Compared with WT cells, there was greater relative percentage of PE (22.55% vs. 16.90%) and lower of PC (49.24% vs. 55.35%) in the BTN1A1(−/−) cell membrane (P < 0.01, Table 2) resulting in a significantly lower ratio of PC/PE (P < 0.001, Fig. 5b). The most different phospholipid species are shown in Table 3, with LPE(18:1) increasing to 4.13-fold and PC(44:1) decreasing to 0.12-fold.
Members of the butyrophilin (BTN) gene family are attracting increased attention because they may play multifunctional roles in a variety of physiological processes including lactation, regulation of T cells in the immune system, and autoimmune diseases . The BTN1A1 protein is associated with the secretion of milk fat. It is highly-expressed in bovine  and caprine  mammary gland during lactation. In the present study, the high mRNA abundance of β-casein confirmed that BMEC maintained lactogenic capacity expected of mammary cells. Subcellular co-localization analyses revealed that BTN1A1 was bound to cytoplasmic LD . Knockout of BTN1A1 in BMEC demonstrated not only an important role of this protein in determining LD size, but also a role in cell phospholipid composition.
The CRISPR/Cas9 system has advantages compared with other genome editing technologies such as transcription activator-like effector nuclease and zinc finger nuclease approaches, including ease of use, high efficiency, and adaptability to different organisms [33, 34]. In the present study, use of the lentiCRISPR v2 vector, which is constructed around a 3rd generation lentiviral backbone, allowed simultaneous infection/transfection of Cas9 and sgRNA followed by selection via puromycin resistance . Editing via CRISP/Cas9 can leave a variety of mutations because when the target region is amplified, and PCR products are denatured and renatured, there are mismatches at the target site. As shown in Fig. 2b, the presence of two shorter bands in agarose gel electrophoresis following product cleavage suggests that CRISPR/Cas9-sgRNA3 successfully introduced insertion or deletion mutations in the genomic target DNA. After selection of a single clone, DNA sequencing confirmed a 10 bp deletion mutation in the BTN1A1 exon (Fig. 2c and d). Furthermore, BTN1A protein abundance was completely blocked in the BTN1A1(−/−) clone. These results indicate that the CRISPR/Cas9 system successfully knocked out the BTN1A1 gene, yielding a homozygous clone.
In the present study, knockout of BTN1A1 gene in BMEC decreased content and diameter of LD. In liver [36, 37] and adipose  cells, the two main factors controlling LD size are phospholipids and proteins on their surface. BTN1A1 is a type I membrane protein that is incorporated into the surface membrane coat surrounding LD . Thus, the fact that BTN1A1 deficiency decreased the size of LD in the present study underscores its importance in the context of LD fusion.
The nuclear receptor PPARG is a member of the PPAR family, which mediates lipid accumulation through the induction of lipogenic genes such as PLIN2 and FASN in mammary cells [26, 39, 40]. Thus, the lower abundance of PPARG and FASN suggests a regulatory role for these genes in LD formation. Some reports also link PPAR activation via phospholipids. For instance, sphingomyelin treatment regulated the expression of PPARG [41, 42]. Shi et al.  reported that activation of PPARD affected abundance of genes related to fatty acid metabolism (e.g. FASN, PLIN2) in goat MEC. In the present study, knockout of BTN1A1 affected PPARG abundance and phospholipid classes, which we speculate provides some evidence of a role for phospholipids on PPAR activity in BMEC. Clearly, additional experiments are needed to better define the mechanisms linking phospholipids and PPAR in the overall process of LD formation.
It is well-known that cellular membranes are made primarily with phospholipids, and alterations in their availability may enable better dispersion of triacylglycerol in the form of smaller MFG. In turn, such effect may greatly affect the size of MFG . It has been reported that the size of LD in BMEC is determined by the phospholipid composition of the cell membrane [9, 44]. The ratio of PC/PE has been suggested as a predictor for LD size. We evaluated the effect of BTN1A1 knockout on membrane phospholipid composition and found that the percentage of PC decreased and the PE increased (Table 2), resulting in a lower PC/PE ratio (Fig. 5b) in the BTN1A1(−/−) cells. Consistent with these results, SCD-deficient C. elegans displayed decreased LD size with lower percentage of PC, higher PE, and lower PC/PE ratio .
Phosphatidylethanolamine can be transformed into PC by the action of transferases, and PE could be produced by decarboxylation of PS . Thus, factors affecting phospholipid synthesis may be involved in the regulation of LD size. The chain length, unsaturation, and polarity of phospholipids affect the diffusion and transport rate of phospholipids in the membrane . We also found that the most different phospholipid species contained monounsaturated fatty acids in the sn-1 position primarily (C16:1 or C18:1, Table 3). These results suggest that a variety of phospholipids may play roles in regulating the size of LD.
The abundance of BTN1A1 gene in ruminant mammary cells is critical for formation of lipid droplets, i.e. regulate the content and size of LD in BMEC. Besides its role in lipid droplet formation, BTN1A1 exerts some degree of control on transcription of lipogenic genes and the profile of cell membrane phospholipids. Although the present study cannot explain the exact mechanisms of BTN1A1 action, the determination of LD size within ruminant mammary cells seems determined by complex regulatory mechanisms including milk fat-globule membrane protein abundance and phospholipid composition of the membranes.
Availability of data and materials
The data during and /or analysed during the current study are available from the corresponding author on reasonable request.
Butyrophilin subfamily 1 member A1
Mammary epithelial cells
Milk fat globule
CRISPR and CRISPR-associated 9 (Cas9)
Fatty acid synthase
Peroxisome proliferator-activated receptor-gamma
McManaman JL, Russell TD, Schaack J, Orlicky DJ, Robenek H. Molecular determinants of milk lipid secretion. J Mammary Gland Biol Neoplasia. 2007;12(4):259–68.
Russell TD, Palmer CA, Orlicky DJ, Fischer A, Rudolph MC, Neville MC, et al. Cytoplasmic lipid droplet accumulation in developing mammary epithelial cells: roles of adipophilin and lipid metabolism. J Lipid Res. 2007;48(7):1463–75.
Chong BM, Reigan P, Mayle-Combs KD, Orlicky DJ, McManaman JL. Determinants of adipophilin function in milk lipid formation and secretion. Trends Endocrinol Metab. 2011;22(6):211–7.
Henry C, Saadaoui B, Bouvier F, Cebo C. Phosphoproteomics of the goat milk fat globule membrane: new insights into lipid droplet secretion from the mammary epithelial cell. Proteomics. 2015;15(13):2307–17.
Nakatani H, Yasueda T, Oshima K, Nadano D, Matsuda T. Re-evaluation of Milk-fat globule EGF-factor VIII (MFG-E8) as an intrinsic component of the mouse Milk-fat globule membrane. Biosci Biotechnol Biochem. 2012;76(11):2055–60.
Liang L, Qi C, Wang X, Jin Q, McClements DJ. Influence of homogenization and thermal processing on the gastrointestinal fate of bovine Milk fat: in vitro digestion study. J Agric Food Chem. 2017;65(50):11109–17.
Minegishi Y, Ota N, Soga S, Shimotoyodome A. Effects of Nutritional Supplementation with Milk Fat Globule Membrane on Physical and Muscle Function in Healthy Adults Aged 60 and Over with Semiweekly Light Exercise: A Randomized Double-Blind, Placebo-Controlled Pilot Trial. J Nutr Sci Vitaminol (Tokyo). 2016;62(6):409–15.
Couvreur S, Hurtaud C. Relationships between milks differentiated on native milk fat globule characteristics and fat, protein and calcium compositions. Animal. 2017;11(3):507–18.
Cohen BC, Shamay A, Argov-Argaman N. Regulation of Lipid Droplet Size in Mammary Epithelial Cells by Remodeling of Membrane Lipid Composition-A Potential Mechanism. PLoS One. 2015;10(3).
Smoczynski M. Role of phospholipid flux during Milk secretion in the mammary gland. J Mammary Gland Biol Neoplasia. 2017;22(2):117–29.
Monks J, Dzieciatkowska M, Bales ES, Orlicky DJ, Wright RM, et al. Xanthine oxidoreductase mediates membrane docking of milk-fat droplets but is not essential for apocrine lipid secretion. J Physiol. 2016;594(20):5899–921.
Kang Y, Hengbo S, Jun L, Jun L, Wangsheng Z, Huibin T, et al. PPARG modulated lipid accumulation in dairy GMEC via regulation of ADRP gene. J Cell Biochem. 2015;116(1):192–201.
Ren CF, Wang LZ, Fan YX, Jia RX, Zhang GM, Deng MT, et al. Scd1 Contributes to Lipid Droplets Formation in GMEC via Transcriptional Regulation of Tip47 and Adrp. Eur J Lipid Sci Technol. 2018;120(2).
Tan R, Wang WJ, Wang SC, Wang Z, Sun LX, He W, et al. Small GTPase Rab40c Associates with Lipid Droplets and Modulates the Biogenesis of Lipid Droplets. PLoS One. 2013;8(4).
Nguyen HTH, Ong L, Hoque A, Kentish SE, Williamson N, Ang CS, et al. A proteomic characterization shows differences in the milk fat globule membrane of buffalo and bovine milk. Food Biosci. 2017;19:7–16.
Mather IH, Keenan TW. Origin and secretion of milk lipids. J Mammary Gland Biol Neoplasia. 1998;3(3):259–73.
Bionaz M, Loor JJ. Gene networks driving bovine milk fat synthesis during the lactation cycle. BMC Genomics. 2008;9:366.
Robenek H, Hofnagel O, Buers I, Lorkowski S, Schnoor M, Robenek MJ, et al. Butyrophilin controls milk fat globule secretion. Proc Natl Acad Sci U S A. 2006;103(27):10385–90.
Mali P, Yang LH, Esvelt KM, Aach J, Guell M, DiCarlo JE, Norville JE. Church GM.RNA-guided human genome engineering via Cas9. Science. 2013;339(6121):823–6.
Zhou Z, Tan H, Li Q, Chen J, Gao S, Wang Y, et al. CRISPR/Cas9-mediated efficient targeted mutagenesis of RAS in Salvia miltiorrhiza. Phytochemistry. 2018;148:63–70.
Cong L, Ran FA, Cox D, Lin SL, Barretto R, Habib N, et al. Multiplex genome engineering using CRISPR/Cas systems. Science. 2013;339(6121):819–23.
Sander JD. Joung JK.CRISPR-Cas systems for editing, regulating and targeting genomes. Nat Biotechnol. 2014;32(4):347–55.
Shalem O, Sanjana NE, Hartenian E, Shi X, Scott DA, Mikkelson T, et al. Genome-scale CRISPR-Cas9 knockout screening in human cells. Science. 2014;343(6166):84–7.
Ma YF, Wu ZH, Gao M, Loor JJ. Nuclear factor erythroid 2-related factor 2-antioxidant activation through the action of ataxia telangiectasia-mutated serine/threonine kinase is essential to counteract oxidative stress in bovine mammary epithelial cells. J Dairy Sci. 2018;101(6):5317–28.
Peterson DG, Matitashvili EA, Bauman DE. The inhibitory effect of trans-10, cis-12 CLA on lipid synthesis in bovine mammary epithelial cells involves reduced proteolytic activation of the transcription factor SREBP-1. J Nutr. 2004;134(10):2523–7.
Kadegowda AK, Bionaz M, Piperova LS, Erdman RA, Loor JJ. Peroxisome proliferator-activated receptor-gamma activation and long-chain fatty acids alter lipogenic gene networks in bovine mammary epithelial cells to various extents. J Dairy Sci. 2009;92(9):4276–89.
Moyer TC, Holland AJ. Generation of a conditional analog-sensitive kinase in human cells using CRISPR/Cas9-mediated genome engineering. Centrosome Centriole. 2015;129:19–36.
Tian H, Luo J, Zhang Z, Wu J, Zhang T, Busato S, et al. CRISPR/Cas9-mediated Stearoyl-CoA Desaturase 1 (SCD1) deficiency affects fatty acid metabolism in goat mammary epithelial cells. J Agric Food Chem. 2018;66(38):10041–52.
Bionaz M, Loor JJ. Identification of reference genes for quantitative real-time PCR in the bovine mammary gland during the lactation cycle. Physiol Genomics. 2007;29(3):312–9.
Han LQ, Gao TY, Yang GY, Loor JJ. Overexpression of SREBF chaperone (SCAP) enhances nuclear SREBP1 translocation to upregulate fatty acid synthase (FASN) gene expression in bovine mammary epithelial cells. J Dairy Sci. 2018;101(7):6523–31.
Smith IA, Knezevic BR, Ammann JU, Rhodes DA, Aw D, Palmer DB, et al. BTN1A1, the mammary gland Butyrophilin, and BTN2A2 are both inhibitors of T cell activation. J Immunol. 2010;184(7):3514–25.
Shi HB, Yu K, Luo J, Li J, Tian HB, Zhu JJ, et al. Loor JJ.Adipocyte differentiation-related protein promotes lipid accumulation in goat mammary epithelial cells. J Dairy Sci. 2015;98(10):6954–64.
Cobb RE, Wang YJ, Zhao HM. High-efficiency multiplex genome editing of Streptomyces species using an engineered CRISPR/Cas system. ACS Synth Biol. 2015;4(6):723–8.
Liu H, Liu C, Zhao YH, Han XJ, Zhou ZW, Wang C, et al. Comparing successful gene knock-in efficiencies of CRISPR/Cas9 with ZFNs and TALENs gene editing systems in bovine and dairy goat fetal fibroblasts. J Integr Agric. 2018;17(2):406–14.
Sanjana NE, Shalem O, Zhang F. Improved vectors and genome-wide libraries for CRISPR screening. Nat Methods. 2014;11(8):783–4.
Suzuki M, Shinohara Y, Ohsaki Y, Fujimoto T. Lipid droplets: size matters. J Electron Microsc. 2011;60:S101–16.
Yu JH, Li P. The size matters: regulation of lipid storage by lipid droplet dynamics. Sci China Life Sci. 2017;60(1):46–56.
Horl G, Wagner A, Cole LK, Malli R, Reicher H, Kotzbeck P, et al. Sequential synthesis and methylation of Phosphatidylethanolamine promote lipid droplet biosynthesis and stability in tissue culture and in vivo. J Biol Chem. 2011;286(19):17338–50.
Shi HB, Luo J, Yao DW, Zhu JJ, Xu HF, Shi HP, et al. Peroxisome proliferator-activated receptor-gamma stimulates the synthesis of monounsaturated fatty acids in dairy goat mammary epithelial cells via the control of stearoyl-coenzyme a desaturase. J Dairy Sci. 2013;96(12):7844–53.
Schadinger SE, Bucher NLR, Schreiber BM. Farmer SR.PPAR gamma 2 regulates lipogenesis and lipid accumulation in steatotic hepatocytes. Am J Physiol Endocrinol Metab. 2005;288(6):E1195–205.
Klingler C, Zhao XJ, Adhikary T, Li J, Xu GW, Haring HU, et al. Lysophosphatidylcholines activate PPAR delta and protect human skeletal muscle cells from lipotoxicity. Biochim Biophys Acta. 2016;1861(12):1980–92.
Norris GH, Porter CM, Jiang C, Millar CL, Blesso CN. Dietary sphingomyelin attenuates hepatic steatosis and adipose tissue inflammation in high-fat-diet-induced obese mice. J Nutr Biochem. 2017;40:36–43.
Shi HB, Zhang CH, Zhao W, Luo J, Loor JJ. Peroxisome proliferator-activated receptor delta facilitates lipid secretion and catabolism of fatty acids in dairy goat mammary epithelial cells. J Dairy Sci. 2017;100(1):797–806.
Cohen BC, Raz C, Shamay A, Argov-Argaman N. Lipid droplet fusion in mammary epithelial cells is regulated by Phosphatidylethanolamine metabolism. J Mammary Gland Biol Neoplasia. 2017;22(4):235–49.
Shi X, Li J, Zou XJ, Greggain J, Rodkaer SV, Faergeman NJ, et al. Regulation of lipid droplet size and phospholipid composition by stearoyl-CoA desaturase. J Lipid Res. 2013;54(9):2504–14.
Pavlovic Z, Bakovic M. Regulation of Phosphatidylethanolamine homeostasis-the critical role of CTP: Phosphoethanolamine Cytidylyltransferase (Pcyt2). Int J Mol Sci. 2013;14(2):2529–50.
We thank Shanghai Applied Protein Technology Co., Ltd., for technical support during quantitative phospholipid analysis by LC-MS/MS.
This research was supported by the National Natural Science Foundation of China (U1904116), Special Funds for Modern Agricultural Industry Technology System (CARS-37), National Key Research and Development Program of China (Beijing, China; 2016YFD0500503).
All of the animal studies were conducted in accordance with the experimental practices and standards approved by the Animal Welfare and Research Ethics Committee at Henan Agricultural University.
Consent for publication
The authors declare that they have no competing interests.
mRNA abundance of β-casein gene in BMEC. Figure S2. Optimum lethal dose of puromycin against BMEC. Figure S3. BTN1A1 protein expression levels in different mutants. Table S1. Phospholipid lipidomics data by LC-MS/MS
About this article
Cite this article
Han, L., Zhang, M., Xing, Z. et al. Knockout of butyrophilin subfamily 1 member A1 (BTN1A1) alters lipid droplet formation and phospholipid composition in bovine mammary epithelial cells. J Animal Sci Biotechnol 11, 72 (2020). https://doi.org/10.1186/s40104-020-00479-6
- Lipid droplet
- Mammary epithelial cell
- Milk fat globule