Effects of polymorphism of the GPAM gene on milk quality traits and its relation to triglyceride metabolism in bovine mammary epithelial cells of dairy cattle

Analysis of miRNAs in milk of four livestock species

BACKGROUND

Milk is essential for mammalian nutrition because it provides vital nutrients for growth and development. Milk composition, which is influenced by genetic and environmental factors, supports lactation, a complex process crucial for milk production and quality. Recent research has focused on noncoding RNAs, particularly microRNAs (miRNAs), which are present in body fluids and regulate gene expression post-transcriptionally. This study comprehensively characterizes miRNAs in milk of four livestock species, namely Bubalus bubalis, Capra hircus, Equus asinus, and Ovis aries and identifies potential target genes.

RESULTS

High-throughput sequencing of milk RNA resulted in distinct read counts across species: B. bubalis (8,790,441 reads), C. hircus (12,976,275 reads), E. asinus (9,385,067 reads), and O. aries (7,295,297 reads). E. asinus had the highest RNA mapping rate (94.6%) and O. aries the lowest (84.8%). A substantially greater proportion of miRNAs over other small RNAs was observed for the donkey milk sample (7.74%) compared to buffalo (0.87%), goat (1.57%), and sheep (1.12%). Shared miRNAs, which included miR-200a, miR-200b, miR-200c, and miR-23a among others, showed varying expression levels across species, confirmed by qPCR analysis. Functional annotation of predicted miRNA target genes highlighted diverse roles, with an enrichment in functions linked to metabolism and immunity. Pathway analysis identified immune response pathways as significant, with several miRNAs targeting specific genes across species, suggesting their regulatory function in milk.

CONCLUSIONS

Both conserved and species-specific miRNAs were detected in milk of the investigated species. The identified target genes of these miRNAs have important roles in neonatal development, adaptation, growth, and immune response. Furthermore, they influence milk and meat production traits in livestock.

Comparison of miR-106b, miR-191, and miR-30d expression dynamics in milk with regard to its composition in Holstein and Ayrshire cows

OBJECTIVE

Milk composition varies considerably and depends on paratypical, genetic, and epigenetic factors. MiRNAs belong to the class of small non-coding RNAs; they are one of the key tools of epigenetic control because of their ability to regulate gene expression at the post-transcriptional level. We compared the relative expression levels of miR-106b, miR-191, and miR-30d in milk to demonstrate the relationship between the content of these miRNAs with protein and fat components of milk in Holstein and Ayrshire cattle.

METHODS

Milk fat, protein, and casein contents were determined in the obtained samples, as well as the content of the main fatty acids (g/100 g milk), including: saturated acids, such as myristic (C14:0), palmitic (C16:0), and stearic (C18:0) acids; monounsaturated acids, including oleic (C18:1) acid; as well as long-, medium- and short-chain, polyunsaturated, and trans fatty acids. Real-time stem-loop one-tube reverse transcription polymerase chain reaction with TaqMan probes was used to measure the miRNA expression levels.

RESULTS

The miRNA expression levels in milk samples were found to be decreased in the first two months in Holstein breed, and in the first four months in Ayrshire breed. Correlation analysis did not reveal any dependence between changes in the expression level of miRNA and milk fat content, but showed a multidirectional relationship with individual milk fatty acids. Positive associations between the expression levels of miR-106b and miR-30d and protein and casein content were found in the Ayrshire breed. Receiver operating characteristic curve analysis showed that miR-106b and miR-30d expression levels can cause changes in fatty acid and protein composition of milk in Ayrshire cows, whereas miR-106b expression level determines the fatty acid composition in Holsteins.

CONCLUSION

The data obtained in this study showed that miR-106b, miR-191, and miR-30d expression levels in milk samples have peculiarities associated with breed affiliation and the lactation period.

Polymorphisms of ITGA9 Gene and Their Correlation with Milk Quality Traits in Yak (Bos grunniens)

A single-nucleotide polymorphism (SNP) is a genome-level trait that arises from a variation in a single nucleotide, leading to diversity in DNA sequences. SNP screening is commonly used to provide candidate genes for yak breeding efforts. Integrin Subunit Alpha 9 (ITGA9) is an integrin protein. It plays an important role in cell adhesion, signalling, and other processes. The aim of this study was to discuss the association between genetic polymorphisms in the ITGA9 gene and milk quality traits and to identify potential molecular marker loci for yak breeding quality. We genotyped 162 yaks using an Illumina Yak cGPS 7K liquid chip and identified the presence of polymorphisms at nine SNP loci in the ITGA9 gene of yaks. The results showed that the mutant genotypes in the loci g.285,808T>A, g.306,600T>C, and g.315,413C>T were positively correlated with the contents of casein, protein, total solids (TS), and solid nonfat (SNF) in yak milk. In other loci, heterozygous genotypes had a positive correlation with nutrient content in yak milk. Then, two ITGA9 haplotype blocks were constructed based on linkage disequilibrium, which facilitated a more accurate screening of ITGA9 as a candidate gene for yak milk quality improvement. In conclusion, we identified SNPs and haplotype blocks related to yak milk quality traits and provided genetic resources for marker-assisted selection in yak breeding.

Potential Role of Lauric Acid in Milk Fat Synthesis in Chinese Holstein Cows Based on Integrated Analysis of Ruminal Microbiome and Metabolome

The composition and metabolic profile of the ruminal microbiome have an impact on milk composition. To unravel the ruminal microbiome and metabolome affecting milk fat synthesis in dairy cows, 16S rRNA and internal transcribed spacer (ITS) gene sequencing, as well as ultra-high performance liquid chromatography-mass spectrometry (UHPLC-MS/MS) methods were used to investigate the significant differences in ruminal bacterial and fungal communities as well as metabolome among Chinese Holstein cows with contrasting milk fat contents under the same diet (H-MF 5.82 ± 0.41% vs. L-MF 3.60 ± 0.12%). Another objective was to culture bovine mammary epithelial cells (BMECs) to assess the effect of metabolites on lipid metabolism. Results showed that the acetate-to-propionate ratio and xylanase activity in ruminal fluid were both higher in H-MF. Microbiome sequencing identified 10 types of bacteria and four types of fungi differently abundant at the genus level. Metabolomics analysis indicated 11 different ruminal metabolites between the two groups, the majority of which were lipids and organic acids. Among these, lauric acid (LA) was enriched in fatty acid biosynthesis with its concentration in milk fat of H-MF cows being greater (217 vs. 156 mg per 100 g milk), thus, it was selected for an in vitro study with BMECs. Exogenous LA led to a marked increase in intracellular triglyceride (TG) content and lipid droplet formation, and it upregulated the mRNA abundance of fatty acid uptake and activation (CD36 and ACSL1), TG synthesis (DGAT1, DGAT2 and GPAM), and transcriptional regulation (SREBP1) genes. Taken together, the greater relative abundance of xylan-fermenting bacteria and fungi, and lower abundance of bacteria suppressing short-chain fatty acid-producing bacteria or participating in fatty acid hydrogenation altered lipids and organic acids in the rumen of dairy cows. In BMECs, LA altered the expression of genes involved in lipid metabolism in mammary cells, ultimately promoting milk fat synthesis. Thus, it appears that this fatty acid plays a key role in milk fat synthesis.

Diacylglycerol acyltransferase 2 promotes the adipogenesis of intramuscular preadipocytes in goat

Diacylglycerol acyltransferase 2 (DGAT2) is the key enzyme that catalyzes the last step of triglyceride synthesis. However, its role in intramuscular fat (IMF) deposition in goat remains unclear. The purpose of this study was to explore the role of DGAT2 in regulating goat IMF deposition. In the present study, the expression of DGAT2 was highest in goat triceps brachii, and highest on the first day after oleic acid induction in goat intramuscular preadipocytes. The overexpression of DGAT2 promoted the accumulation of lipid droplets and triglyceride synthesis, accompanied by the expression upregulation of DGAT1, TIP47, ACC and ACOX1 significantly, and expression downregulation of AGPAT6, LPIN1, LPL, HSL, ATGL and ADRP significantly. In contrast, the silencing of DGAT2 decreased the accumulation of lipid droplets, inhibited the expression of DGAT1, GPAM, ADRP, AGPAT6, LPL, HSL, ATGL, ACC, FASN, ACOX1 significantly, and enhanced that of TIP47 significantly. Overall, these data underscore DGAT2 may play a potentially important role in lipid droplets formation and triglyceride accumulation, so as to maintain intramuscular fat deposition, beyond triglyceride synthesis in goat.

Development of TaqMan PCR assay for genotyping SNP rs211250281 of the bovine agpat6 gene

Milk fat percentage is an important production trait of dairy cattle and is one of the goals of breeding programs. Over 95% of the milk fat accounts for triacylglycerols. AGPAT6 (1-acylglycerol-3-phosphate O-acyltransferase 6) catalyzes an intermediary step of triglyceride synthesis in the mammary cells. Genome-wide association studies identified SNP rs211250281 (g27: 36520069 G/T) in the agpat6 gene associated with milk fat content and fat-to-protein ratio in dairy cattle. The article presents data on the development of TaqMan PCR assay for genotyping SNP rs211250281 of the bovine agpat6 gene. In this method, a primer pair, initiating amplification of 75-bp fragments of the agpat6 gene, and two allele-specific TaqMan probes are used. Identification of the G and T alleles is based on a comparison of the final fluorescence intensity of FAM and VIC dyes, respectively. The developed TaqMan PCR assay was validated by Sanger sequencing method. The results of both methods fully coincided, that confirmed high accuracy of the developed TaqMan PCR assay. This reliable, simple, rapid, and high-throughput method could be a suitable tool for studying the distribution of the SNP rs211250281 among different cattle breeds and its association with milk fat content.

Review: genetic background of milk fatty acid synthesis in bovines

Milk fat composition is an important trait for the dairy industry as it directly influences the nutritional and technological properties of milk and other dairy products. The synthesis of milk fat is a complex process regulated by a network of genes. Thus, understanding the genetic variation and molecular mechanisms regulating milk fat synthesis will help to improve the nutritional quality of dairy products. In this review, we provide an overview of milk fat synthesis in bovines along with the candidate genes involved in the pathway. We also discuss de novo synthesis of fatty acids (ACSS, ACACA, FASN), uptake of FAs (FATP, FAT, LPL), intracellular activation and channelling of FAs (ACSL, FABP), elongation (EVOLV6), desaturation (SCD, FADS), formation of triglycerides (GPAM, AGPAT, LIPIN, DGAT), and milk lipid secretion (BTN1A1, XDH, PLIN2). The genetic variability of individual fatty acids will help to develop selection strategies for obtaining a healthier milk fat profile in bovines. Thus, this review will offer a potential understanding of the molecular mechanisms that regulate milk fat synthesis in bovines.

ranchSATdb: A Genome-Wide Simple Sequence Repeat (SSR) Markers Database of Livestock Species for Mutant Germplasm Characterization and Improving Farm Animal Health

Microsatellites, also known as simple sequence repeats (SSRs), are polymorphic loci that play an important role in genome research, animal breeding, and disease control. Ranch animals are important components of agricultural landscape. The ranch animal SSR database, ranchSATdb, is a web resource which contains 15,520,263 putative SSR markers. This database provides a comprehensive tool for performing end-to-end marker selection, from SSRs prediction to generating marker primers and their cross-species feasibility, visualization of the resulting markers, and finding similarities between the genomic repeat sequences all in one place without the need to switch between other resources. The user-friendly online interface allows users to browse SSRs by genomic coordinates, repeat motif sequence, chromosome, motif type, motif frequency, and functional annotation. Users may enter their preferred flanking area around the repeat to retrieve the nucleotide sequence, they can investigate SSRs present in the genic or the genes between SSRs, they can generate custom primers, and they can also execute in silico validation of primers using electronic PCR. For customized sequences, an SSR prediction pipeline called miSATminer is also built. New species will be added to this website's database on a regular basis throughout time. To improve animal health via genomic selection, we hope that ranchSATdb will be a useful tool for mapping quantitative trait loci (QTLs) and marker-assisted selection. The web-resource is freely accessible at https://bioinfo.usu.edu/ranchSATdb/.

miR-497 Regulates LATS1 through the PPARG Pathway to Participate in Fatty Acid Synthesis in Bovine Mammary Epithelial Cells

Nutrient metabolism is required to maintain energy balance in animal organisms, and fatty acids play an irreplaceable role in fat metabolism. In this study, microRNA sequencing was performed on mammary gland tissues collected from cows during early, peak, and late lactation to determine miRNA expression profiles. Differentially expressed miRNA (miR-497) was selected for functional studies of fatty acid substitution. Simulants of miR-497 impaired fat metabolism [triacylglycerol (TAG) and cholesterol], whereas knockdown of miR-497 promoted fat metabolism in bovine mammary epithelial cells (BMECs) in vitro. In addition, in vitro experiments on BMECs showed that miR-497 could down-regulate C16:1, C17:1, C18:1, and C20:1 as well as long-chain polyunsaturated fats. Thus, these data expand the discovery of a critical role for miR-497 in mediating adipocyte differentiation. Through bioinformatics analysis and further validation, we identified large tumor suppressor kinase 1 (LATS1) as a target of miR-497. siRNA-LATS1 increased concentrations of fatty acids, TAG, and cholesterol in cells, indicating an active role of LATS1 in milk fat metabolism. In summary, miR-497/LATS1 can regulate the biological processes associated with TAG, cholesterol, and unsaturated fatty acid synthesis in cells, providing an experimental basis for further elucidating the mechanistic regulation of lipid metabolism in BMECs.

Yi-Shen-Hua-Shi granule ameliorates diabetic kidney disease by the "gut-kidney axis"

ETHNOPHARMACOLOGICAL RELEVANCE

Yi-Shen-Hua-Shi (YSHS) granule is an effective prescription widely used in traditional Chinese medicine to treat diabetic kidney disease (DKD), its exact efficacy in treating DKD has been confirmed but the underlying regulatory mechanism has not been fully elucidated.

AIM OF THE STUDY

To explore the mechanism by which YSHS granule regulates intestinal flora and serum metabolites and then regulates renal mRNA expression through the "gut-kidney axis", so as to improve DKD.

MATERIALS AND METHODS

40 rats were divided into five groups: Normal group (N) (normal saline), model group (M) (STZ + normal saline), YSHS granule low-dose group (YL) (STZ + 2.27 g kg^-1 d^-1), YSHS granule high-dose group (YH) (STZ + 5.54g kg^-1 d^-1) and valsartan group (V) (STZ + 7.38mg kg^-1 d^-1). After 6 weeks, changes in blood glucose, blood lipids, and renal function related indexes were observed, as well as pathological changes in the kidney and colon. Intestinal microbiota was sequenced by 16S rDNA, serum differential metabolites were identified by LC-MS/MS, and renal differences in mRNA expression were observed by RNA-seq. Further, through the association analysis of intestinal differential microbiota, serum differential metabolites and kidney differential mRNAs, the target flora, target metabolites and target genes of YSHS granule were screened and verified, and the "gut-metabolism-transcription" co-expression network was constructed.

RESULTS

In group M, blood glucose, blood lipid and proteinuria were increased, inflammation, oxidative stress and renal function were aggravated, with the proliferation of mesangial matrix, vacuolar degeneration of renal tubules, accumulation of collagen and lipid, and increased intestinal permeability, and YSHS granule and valsartan improved these disorders to varying degrees. High dose of YSHS granule improved the diversity and abundance of flora, decreased the F/B value, greatly increased the abundance of Lactobacillus and Lactobacillus_murinus, and decreased the abundance of Prevoella UCG_001. 14 target metabolites of YSHS granule were identified, which were mainly enriched in 20 KEGG pathways, such as Glycerophospholipid metabolism, Sphingolipid metabolism and Phenylalanine, tyrosine and tryptophan biosynthesis. 96 target mRNAs of YSHS granule were also identified. The enriched top 20 pathways were closely related to glucose and lipid metabolism, of which a total of 21 differential mRNAs were expressed. Further correlation analysis revealed that Lactobacillus, Lactobacillus_murinus and Prevotella UCG_001 were highly correlated with Glycerophospholipid metabolism, Sphingolipid metabolism and Phenylalanine, tyrosine and tryptophan biosynthesis pathways. At the same time, 6 pathways including Glycerophospholipid metabolism, Arachidonic acid metabolism, Purine metabolism, Primary bile acid biosynthesis, Ascorbate and aldarate metabolism and Galactose metabolism were co-enriched by the target metabolites and the target mRNAs of YSHS granule, including 7 differential metabolites such as phosphatidylethanolamine and 7 differential genes such as Adcy3. The 7 differential metabolites had high predictive value of AUC, and the validation of 7 differential genes were highly consistent with the sequencing results.

CONCLUSION

YSHS granule could improve DKD through the "gut-kidney axis". Lactobacillus and Lactobacillus_murinus were the main driving forces. 6 pathways related to glucose and lipid metabolism, especially Glycerophospholipid metabolism, may be an important follow-up response and regulatory mechanism.

Feed efficiency in dairy sheep: An insight from the milk transcriptome

Introduction

As higher feed efficiency in dairy ruminants means a higher capability to transform feed nutrients into milk and milk components, differences in feed efficiency are expected to be partly linked to changes in the physiology of the mammary glands. Therefore, this study aimed to determine the biological functions and key regulatory genes associated with feed efficiency in dairy sheep using the milk somatic cell transcriptome.

Material and methods

RNA-Seq data from high (H-FE, n = 8) and low (L-FE, n = 8) feed efficiency ewes were compared through differential expression analysis (DEA) and sparse Partial Least Square-Discriminant analysis (sPLS-DA).

Results

In the DEA, 79 genes were identified as differentially expressed between both conditions, while the sPLS-DA identified 261 predictive genes [variable importance in projection (VIP) > 2] that discriminated H-FE and L-FE sheep.

Discussion

The DEA between sheep with divergent feed efficiency allowed the identification of genes associated with the immune system and stress in L-FE animals. In addition, the sPLS-DA approach revealed the importance of genes involved in cell division (e.g., KIF4A and PRC1) and cellular lipid metabolic process (e.g., LPL, SCD, GPAM, and ACOX3) for the H-FE sheep in the lactating mammary gland transcriptome. A set of discriminant genes, commonly identified by the two statistical approaches, was also detected, including some involved in cell proliferation (e.g., SESN2, KIF20A, or TOP2A) or encoding heat-shock proteins (HSPB1). These results provide novel insights into the biological basis of feed efficiency in dairy sheep, highlighting the informative potential of the mammary gland transcriptome as a target tissue and revealing the usefulness of combining univariate and multivariate analysis approaches to elucidate the molecular mechanisms controlling complex traits.

Deciphering the Key Regulatory Roles of KLF6 and Bta-miR-148a on Milk Fat Metabolism in Bovine Mammary Epithelial Cells

MicroRNAs (miRNAs) are non-coding RNAs that regulate the expression of their target genes involved in many cellular functions at the post-transcriptional level. Previously, bta-miR-148a showed significantly high expression in bovine mammary epithelial cells (BMECs) of Chinese Holstein cows producing high milk fat compared to those with low milk fat content. Here, we investigated the role of bta-miR-148a through targeting Krüppel-like factor 6 (KLF6) and further analyzed the role of KLF6 in regulating fat metabolism through targeting PPARA, AMPK/mTOR/PPARG, and other fat marker genes in BMECs of Chinese Holstein. The bioinformatics analysis showed that the 3’ UTR of KLF6 mRNA possesses the binding sites for bta-miR-148a, which was further verified through dual-luciferase reporter assay. The BMECs were transfected with bta-miR-148a-mimic, inhibitor, and shNC, and the expression of KLF6 was found to be negatively regulated by bta-miR-148a. Moreover, the contents of triglyceride (TG), and cholesterol (CHO) in BMECs transfected with bta-miR-148a-mimic were significantly lower than the contents in BMECs transfected with bta-miR-148a-shNC. Meanwhile, the TG and CHO contents were significantly increased in BMECs transfected with bta-miR-148a-inhibitor than in BMECs transfected with bta-miR-148a-shNC. In addition, the TG and CHO contents were significantly decreased in BMECs upon the down-regulation of KLF6 through transfection with pb7sk-KLF6-siRNA1 compared to the control group. Contrarily, when KLF6 was overexpressed in BMECs through transfection with pBI-CMV3-KLF6, the TG and CHO contents were significantly increased compared to the control group. Whereas, the qPCR and Western blot evaluation of PPARA, AMPK/mTOR/PPARG, and other fat marker genes revealed that all of the genes were considerably down-regulated in the KLF6-KO-BMECs compared to the normal BMECs. Taking advantage of deploying new molecular markers and regulators for increasing the production of better-quality milk with tailored fat contents would be the hallmark in dairy sector. Hence, bta-miR-148a and KLF6 are potential candidates for increased milk synthesis and the production of valuable milk components in dairy cattle through marker-assisted selection in molecular breeding. Furthermore, this study hints at the extrapolation of a myriad of functions of other KLF family members in milk fat synthesis.

Acetate-Induced Milk Fat Synthesis Is Associated with Activation of the mTOR Signaling Pathway in Bovine Mammary Epithelial Cells

Acetate is a precursor substance for fatty acid synthesis in bovine mammary epithelial cells (BMECs), and the mTOR signaling pathway plays an important role in milk fat synthesis. However, the mechanism of the regulatory effects of acetate on lipogenic genes via the mTOR signaling pathway in BMEC remains unknown. We hypothesized that acetate can enhance the expression of lipogenic genes and triglyceride (TG) production by activating the mTOR signaling pathway in BMECs. Therefore, the aim of this study was to investigate the network of acetate-regulated lipid metabolism by the mTOR signaling pathway in BMECs. These results showed that TG synthesis was elevated (p < 0.01) in BMECs with acetate treatment. The lipid droplets were increased in the acetate-treated groups compared with those in the control group through the Bodipy staining of the lipids. In addition, the fatty acid profile in BMECs treated with acetate was affected, with an elevation in the proportions of C14:0, C16:0, and C18:0. The mRNA levels of the sterol-response-element-binding protein 1 (SREBP1), stearoyl-CoA desaturase 1 (SCD1), and fatty acid synthase (FAS) genes involved in the lipogenesis and transcriptional factors were upregulated (p < 0.05) in BMECs with acetate treatment. Remarkably, the expression of acetyl-CoA carboxylase α (ACCα) and FAS rate-limiting enzymes involved in lipogenesis was upregulated in BMECs with acetate treatment. Moreover, the addition of acetate enhanced the key protein expression of S6K1, which is related to the mTOR signaling pathway. Taken together, our data suggest that TG accumulation and expression of lipogenic genes induced by acetate are associated with the activation of the mTOR signaling pathway, which provides new insights into the understanding of the molecular mechanism in the expression of mTOR-signaling-pathway-regulated lipogenic genes.

Fatty acid metabolism-related genes are associated with flavor-presenting aldehydes in Chinese local chicken

Aldehydes are primary volatile organic compounds (VOCs) in local Chinese chicken meat and contribute green grass, fatty, citrus, and bitter almond aromas to chicken meat. To understand the genetic basis of these aldehyde VOC aromas, we used approximately 500 Chinese Jingxing Yellow (JXY) chickens to conduct genome-wide association studies (GWAS) on the flavor traits with the data of single nucleotide polymorphisms (SNPs) and insertions and deletions (INDELs). In total, 501 association variants (253 SNPs and 248 INDELs) were found to be suggestively (SNPs: p-value < 2.77e-06 and INDELs: p-value < 3.78e-05) associated with total aldehydes (the sum of nine aldehydes), hexanal, heptanal, benzaldehyde, (E,E)-2,4-nonadienal, octanal, (E)-2-decenal, nonanal, decanal, and octadecanal. Of them, six SNPs and 23 INDELs reached a genome-wide significance level (SNPs: p-value < 1.38e-07 and INDELs: p-value < 1.89e-06). Potential candidate aldehyde genes were functionally annotated for lipid metabolism, especially fatty acid-related pathways and phospholipid-related gene ontology (GO) terms. Moreover, the GWAS analysis of total aldehydes, hexanal, and nonanal generated the most significant signals, and phenotypic content differed between different genotypes at candidate gene-related loci. For total aldehydes and hexanal traits, candidate genes were annotated based on the significant and suggestive variants on chromosomes 3 and 8 with highly polymorphic linkage blocks. The following candidate genes were also identified: GALM, MAP4K3, GPCPD1, RPS6KA2, CRLS1, ASAP1, TRMT6, SDC1, PUM2, ALDH9A1, MGST3, GMEB1, MECR, LDLRAP1, GPAM and ACSL5. We also found that polyunsaturated fatty acids (PUFAs) (C18:2n6c linoleic acid and C18:3n3 linolenic acid) were significantly correlated with total aldehydes and hexanal contents. PUFAs are important aldehyde precursors, and consistently, our results suggested that candidate genes involved in fatty acid pathways and phospholipid GO terms were identified in association loci. This work provides an understanding of the genetic basis of aldehyde formation, which is a key flavor-forming compound.

Transcriptome analysis of CRISPR/Cas9-mediated GPAM-/- in bovine mammary epithelial cell-line unravelled the effects of GPAM gene on lipid metabolism

Glycerol-3-phosphate acyltransferase mitochondrial (GPAM) is an enzyme in animal lipid metabolism pathways that catalyzes the initial and most committed step of glycerolipid biosynthesis. The present study mainly focused on exploring the relationship between the GPAM gene and the lipid metabolism of mammary epithelial cells and the effect of GPAM on the related pathways of lipid metabolism. The GPAM gene was knocked out entirely in bovine mammary epithelial cells(BMECs) using CRISPR/Cas9 technology, and the mechanism by which the GPAM gene regulates lipid metabolism in BMECs was confirmed. Furthermore, after the complete loss of GPAM, BMECs' triglycerides (TGs) and cholesterol (CHOL) levels were significantly decreased (p < 0.05). Concurrently, the content of octanoic acid, a medium-chain saturated fatty acid, increased substantially in BMECs. RNA-seq of GPAM^-/- BMECs revealed that GPAM could affect the expression of genes related to lipid metabolism, downregulated the expression of Acyl-CoA synthetase long-chain family member 5 (ACSL5), Fatty Acid Binding Protein 3 (FABP3), Hormone-sensitive lipase (HSL), Protease, serine-2 (PRSS2), 1-Acylglycerol-3-Phosphate O Acyltransferase 4 (AGPAT4), and regulated the milk synthesis metabolism pathway.The findings revealed that a number of genes were expressed, a number of genes were differentially expressed genes (DEGs), and a number of GO terms were enriched, with a number of GO terms considerably increased. Further, the differentially expressed genes (DEGs) were significantly enriched in Fat digestion and absorption pathway, Fatty acid metabolic pathway, Biosynthesis of unsaturated fatty acids, Biosynthesis of unsaturated fatty acids and steroids, NF-kappa B signalling pathway, MAPK signalling pathway. In conclusion, the current research results show that GPAM is a crucial regulator of BMEC lipid metabolism. GPAM^-/- BMEC may also become useful genetic materials and tools for future research on gene functions related to lipid and fatty acid metabolism. This study will contribute to the discovery of gene regulation and molecular mechanisms in milk fat synthesis.