The SNPs in the ACACA gene are effective on fatty acid composition in Holstein milk

Whole Genome Scan Uncovers Candidate Genes Related to Milk Production Traits in Barka Cattle

In this study, our primary aim was to explore the genomic landscape of Barka cattle, a breed recognized for high milk production in a semi-arid environment, by focusing on genes with known roles in milk production traits. We employed genome-wide analysis and three selective sweep detection methods (ZFST, θπ ratio, and ZHp) to identify candidate genes associated with milk production and composition traits. Notably, ACAA1, P4HTM, and SLC4A4 were consistently identified by all methods. Functional annotation highlighted their roles in crucial biological processes such as fatty acid metabolism, mammary gland development, and milk protein synthesis. These findings contribute to understanding the genetic basis of milk production in Barka cattle, presenting opportunities for enhancing dairy cattle production in tropical climates. Further validation through genome-wide association studies and transcriptomic analyses is essential to fully exploit these candidate genes for selective breeding and genetic improvement in tropical dairy cattle.

Genomic analysis uncovers novel candidate genes related to adaptation to tropical climates and milk production traits in native goats

BACKGROUND

Since domestication, both evolutionary forces and human selection have played crucial roles in producing adaptive and economic traits, resulting in animal breeds that have been selected for specific climates and different breeding goals. Pakistani goat breeds have acquired genomic adaptations to their native climate conditions, such as tropical and hot climates. In this study, using next-generation sequencing data, we aimed to assess the signatures of positive selection in three native Pakistani goats, known as milk production breeds, that have been well adapted to their local climate.

RESULTS

To explore the genomic relationship between studied goat populations and their population structure, whole genome sequence data from native goat populations in Pakistan (n = 26) was merged with available worldwide goat genomic data (n = 184), resulting in a total dataset of 210 individuals. The results showed a high genetic correlation between Pakistani goats and samples from North-East Asia. Across all populations analyzed, a higher linkage disequilibrium (LD) level (- 0.59) was found in the Pakistani goat group at a genomic distance of 1 Kb. Our findings from admixture analysis (K = 5 and K = 6) showed no evidence of shared genomic ancestry between Pakistani goats and other goat populations from Asia. The results from genomic selection analysis revealed several candidate genes related to adaptation to tropical/hot climates (such as; KITLG, HSPB9, HSP70, HSPA12B, and HSPA12B) and milk production related-traits (such as IGFBP3, LPL, LEPR, TSHR, and ACACA) in Pakistani native goat breeds.

CONCLUSIONS

The results from this study shed light on the structural variation in the DNA of the three native Pakistani goat breeds. Several candidate genes were discovered for adaptation to tropical/hot climates, immune responses, and milk production traits. The identified genes could be exploited in goat breeding programs to select efficient breeds for tropical/hot climate regions.

Identification of Copy Number Variations and Selection Signatures in Wannan Spotted Pigs by Whole Genome Sequencing Data: A Preliminary Study

Copy number variation (CNV) is an important structural variation used to elucidate complex economic traits. In this study, we sequenced 25 Wannan spotted pigs (WSPs) to detect their CNVs and identify their selection signatures compared with those of 10 Asian wild boars. A total of 14,161 CNVs were detected in the WSPs, accounting for 0.72% of the porcine genome. The fixation index (Fst) was used to identify the selection signatures, and 195 CNVs with the top 1% of the Fst value were selected. Eighty genes were identified in the selected CNV regions. Functional GO and KEGG analyses revealed that the genes within these selected CNVs are associated with key traits such as reproduction (GAL3ST1 and SETD2), fatty acid composition (PRKG1, ACACA, ACSL3, UGT8), immune system (LYZ), ear size (WIF1), and feed efficiency (VIPR2). The findings of this study contribute novel insights into the genetic CNVs underlying WSP characteristics and provide essential information for the protection and utilization of WSP populations.

Integrated transcriptomic and proteomic analyses reveal the mechanism of easy acceptance of artificial pelleted diets during food habit domestication in Largemouth bass (Micropterus salmoides)

Acceptance of artificial pelleted diets contributes to increasing the cultured areas and output of carnivorous fish. However, the mechanism of acceptance of artificial pelleted diets remains largely unknown. In this study, the easy acceptance of artificial pelleted diets (EAD) group and the not easy acceptance of artificial pelleted diets (NAD) group of Largemouth bass (Micropterus salmoides) were divided based on the ratios of stomach weight/body weight (SB) after 0.5 h feeding, which was bigger than 18% in the EAD group and ranged from 8 to 12% in the NAD group. Through transcriptome and proteome sequencing, a total of 2463 differentially expressed genes (DEGs) and 230 differentially expressed proteins (DEPs) were identified, respectively. Integrated analyses of transcriptome and proteome data revealed that 152 DEPs were matched with the corresponding DEGs (named co-DEGs-DEPs), and 54 co-DEGs-DEPs were enriched in 16 KEGG pathways, including the metabolic pathways, steroid biosynthesis, fatty acid biosynthesis, etc. Furthermore, 3 terpenoid backbone biosynthesis-related genes (Hmgcr, Hmgcs, and Fdps) in metabolic pathways, 10 steroid biosynthesis-related genes (Fdft1, Sqle, Lss, Cyp51a1, Tm7sf2, Nsdhl, Hsd17b7, Dhcr24, Sc5d, and Dhcr7), and 3 fatty acid biosynthesis-related genes (Acaca, Fasn, and Ascl) were all up-regulated in the EAD group, suggesting that the lipid metabolism pathway and steroid biosynthesis pathway play important roles in early food habit domestication in Largemouth bass. In addition, the detection results of randomly selected 15 DEGs and 15 DEPs indicated that both transcriptome and proteome results in the study were reliable. Our study provides useful information for further research on the mechanisms of food habit domestication in fish.

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.

Identification of Genetic Polymorphisms of PI, PIII, and Exon 53 in the Acetyl-CoA Carboxylase-α (ACACα) Gene and Their Association with Milk Composition Traits of Najdi Sheep

Recently, increasing attention has been paid to sheep milk products, which are high in saturated fatty acids (SFA), and the extent of their impact on human health. This study aimed to identify SNPs for PI, PIII, and Exon 53 in the ACACα gene and their association with the MC and FA profiles in Najdi sheep milk. A total of 76 multiparous Najdi ewes were used, and they were maintained using the same feeding system. Milk and blood samples were collected during the first lactation. A genetic polymorphism analysis identified 20 SNPs: 4 SNPs on PI, 6 SNPs on PIII, and 10 SNPs on Exon 53. In PI, the SNP g.4412G > A was associated (p < 0.05) with palmitic acid (C16:0), palmitoleic acid (16:1 n-7) and linoleic acid (LA), while SNP g.4485C > G was associated with CLA and vaccenic acid (VA) (p < 0.05). Furthermore, in PIII, two SNPs (g.1168A > G and g.1331G > T) were associated with milk protein (p < 0.05), while the SNP g.6860G > C in Exon 53 was associated with milk fat (p < 0.05). SNPs in the Najdi breed have been shown to be strongly related to milk fat and EFA contents. This could support a genetic selection program and the control of milk traits in the Najdi breed of high-quality dairy sheep.

Using QTL to Identify Genes and Pathways Underlying the Regulation and Production of Milk Components in Cattle

Milk is a complex liquid, and the concentrations of many of its components are under genetic control. Many genes and pathways are known to regulate milk composition, and the purpose of this review is to highlight how the discoveries of quantitative trait loci (QTL) for milk phenotypes can elucidate these pathways. The main body of this review focuses primarily on QTL discovered in cattle (Bos taurus) as a model species for the biology of lactation, and there are occasional references to sheep genetics. The following section describes a range of techniques that can be used to help identify the causative genes underlying QTL when the underlying mechanism involves the regulation of gene expression. As genotype and phenotype databases continue to grow and diversify, new QTL will continue to be discovered, and although proving the causality of underlying genes and variants remains difficult, these new data sets will further enhance our understanding of the biology of lactation.

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.

Integrating genomics and transcriptomics to identify candidate genes for subcutaneous fat deposition in beef cattle

Fat deposition is a complex economic trait regulated by polygenic genetic basis and environmental factors. Therefore, integrating multi-omics data to uncover its internal regulatory mechanism has attracted extensive attention. Here, we performed genomics and transcriptomics analysis to detect candidates affecting subcutaneous fat (SCF) deposition in beef cattle. The association of 770K SNPs with the backfat thickness captured nine significant SNPs within or near 11 genes. Additionally, 13 overlapping genes regarding fat deposition were determined via the analysis of differentially expressed genes and weighted gene co-expression network analysis (WGCNA). We then calculated the correlations of these genes with BFT and constructed their interaction network. Finally, seven biomarkers including ACACA, SCD, FASN, ACOX1, ELOVL5, HACD2, and HSD17B12 were screened. Notably, ACACA, identified by the integration of genomics and transcriptomics, was more likely to exert profound effects on SCF deposition. These findings provided novel insights into the regulation mechanism underlying bovine fat accumulation.

Restoration of the adipogenic gene expression by naringenin and naringin in 3T3-L1 adipocytes

Background

Naringenin and its glycoside naringin are well known citrus flavonoids with several therapeutic benefits. Although the anti-adipogenic effects of naringenin and naringin have been reported previously, the detailed mechanism underlying their anti-adipogenesis effects is poorly understood.

Objectives

This study examined the anti-adipogenic effects of naringenin and naringin by determining differential gene expression patterns in these flavonoids-treated 3T3-L1 adipocytes.

Methods

Lipid accumulation and triglyceride (TG) content were determined by Oil red O staining and TG assay. Glucose uptake was measured using a 2-[N-(7-Nitrobenz-2-oxa-1,3-diazol-4-yl)amino]-2-deoxy-d-glucose fluorescent d-glucose analog. The phosphorylation levels of AMP-activated protein kinase (AMPK) and acetyl Co-A carboxylase (ACC) were observed via Western blot analysis. Differential gene expressions in 3T3-L1 adipocytes were evaluated via RNA sequencing analysis.

Results

Naringenin and naringin inhibited both lipid accumulation and TG content, increased phosphorylation levels of both AMPK and ACC and decreased the expression level of 3-hydroxy-3-methylglutaryl CoA reductase (HMGCR) in 3T3-L1 adipocytes. RNA sequencing analysis revealed that 32 up-regulated (> 2-fold) and 17 down-regulated (< 0.6-fold) genes related to lipid metabolism, including Acaca, Fasn, Scd1, Mogat1, Dgat, Lipin1, Cpt1a, and Lepr, were normalized to the control level in naringenin-treated adipocytes. In addition, 25 up-regulated (> 2-fold) and 25 down-regulated (< 0.6-fold) genes related to lipid metabolism, including Acaca, Fasn, Fabp5, Scd1, Srebf1, Hmgcs1, Cpt1c, Lepr, and Lrp1, were normalized to the control level by naringin.

Conclusions

The results indicate that naringenin and naringin have anti-adipogenic potentials that are achieved by normalizing the expression levels of lipid metabolism-related genes that were perturbed in differentiated 3T3-L1 cells.

RNA-Seq based genetic variant discovery provides new insights into controlling fat deposition in the tail of sheep

Genetic basis of fat deposition in sheep tail have not been completely elucidated yet. Understanding the genetic mechanisms controlling fat-tail size can improve breeding strategies to modulate fat deposition. RNA sequencing has made it possible to discover genetic variants that may underlie various phenotypic differences. Hence, to identify genetic variants that are important for describing different fat-tail phenotypes in sheep, RNA sequencing was used for single nucleotide polymorphism (SNP) calling in two Iranian sheep breeds (Lori-Bakhtiari, fat-tailed; n = 4, vs Zel, thin-tailed; n = 4). Using a stringent pipeline, a total of 112,344 known SNPs were genotyped, of which 30,550 and 42,906 SNPs were shared by at least two Lori-Bakhtiari and Zel, respectively. Comparing these SNPs showed 2,774 (including 209 missense and 25 deleterious SNPs) and 10,470 (including 1,054 missense and 116 deleterious SNPs) breed-specific SNPs in Lori-Bakhtiari and Zel sheep, respectively. Potential breed-specific SNPs were detected by considering those located in QTL regions associated with fatness or reported as important candidates in previous similar studies. Of the breed-specific SNPs, 724 and 2,905 were located in the QTL regions. Functional enrichment analysis of the affected genes revealed several enriched gene ontologies and KEGG pathways related to fat metabolism. Based on the results, several affected genes were proposed to be strongly linked with fat deposition such as DGAT2, ACSL1, ACACA, ADIPOQ, ACLY, FASN, CPT2, SCD, ADCY6, PER3, CSF1R, SLC22A4, GFPT1, CDS2, BMP6, ACSS2, ELOVL6, HOXA10 and FABP4. Moreover, several SNPs were found in the candidate genes related to fatty acid oxidation introducing them as promising candidates responsible for lower fat content in tail of Zel. Our findings provided new insights into the genetic mechanisms of fat deposition in sheep, which can serve to designing appropriate breeding programs.

Current Evidences and Future Perspectives for AMPK in the Regulation of Milk Production and Mammary Gland Biology

Activated protein kinase (AMP)-activated protein kinase (AMPK) senses the cellular energy status and coordinates catabolic and anabolic processes. Extensive studies have proposed that AMPK regulates energy homeostasis, cell growth, autophagy, mitochondrial biology and inflammation. The biological functions of AMPK vary in different tissues or organs. As a unique organ that produces milk, the mammary gland has recently attracted substantial research attention. This review discusses how AMPK in the mammary gland is activated by energy deprivation and heat stress via the activation of canonical and non-canonical pathways. In addition, the important downstream targets of AMPK and their functions in the mammary gland, especially during milk synthesis, are summarized in the review.

Impact of SNPs in ACACA, SCD1, and DGAT1 Genes on Fatty Acid Profile in Bovine Milk with Regard to Lactation Phases

Simple Summary

Fatty acids are an important component of milk fat. Because of their wide spectrum of effects on human health, it is important to better understand the regulation of their profile in milk. This study aims to analyzing the relation between selected genes with milk fatty acid content. As increased concentration of unhealthy fatty acids and lower concentration of healthy ones in less frequent homozygotes and a strong influence of the genes on fatty acids with 18 carbon atoms were observed, these findings could be useful in future dairy cattle selection aiming production of more healthy milk.

Abstract

Milk fat is a dietary source of fatty acids (FA), which can be health promoting or can increase risks of some diseases. FA profile composition depends on many factors, among them gene polymorphism. This study analyzed the relation between polymorphism of acetyl-CoA carboxylase α (ACACA), stearoyl-CoA desaturase 1 (SCD1), diacylglycerol acyltransferase 1 (DGAT1) genes with FA profile in milk from Polish Holstein-Friesian cattle and determined changes of FA percentage during lactation with regard to polymorphism. Milk samples were collected twice: during the first phase of lactation (<90 Days in milk; DIM) and at the end of lactation (>210 DIM). During the first milk collection, blood samples were taken to analyze three chosen single nucleotide polymorphisms (SNPs): AJ312201.1g.1488C > G SNP in ACACA gene, A293V SNP in SCD1 gene, and K232A SNP in DGAT1 gene. Increased concentration of FA that are less beneficial for human health and have lower concentration of healthy FA in homozygotes: GG in ACACA, VV in SCD1, and KK in DGAT1 were observed, as well as a strong influence of the analyzed genes on FA with 18C atoms was also found. Moreover, it was demonstrated that lactation phase significantly affected FA percentage in milk depending on the phenotype. These results may contribute their part to knowledge toward obtaining more beneficial milk composition.

Comparative analysis of the skim milk and milk fat globule membrane proteomes produced by Jersey cows grazing pastures with different plant species diversity

The objective of this experiment was to identify and characterize the bovine milk proteome within the skim milk fraction and milk fat globule membrane (MFGM)-associated fraction from 16 organically certified lactating Jersey cows after a short term of grazing pastures with or without annual forage crops (AFC). Cows were offered a partial mixed ration (~60% of dry matter intake) and approximately 40% of their total dry matter intake as herbage. Eight cows were offered a cool-season grass–legume herbage (GLH), which included orchardgrass (Dactylis glomerata), timothy (Phleum pratense), Kentucky bluegrass (Poa pratensis), and white clover (Trifolium repens). The other 8 cows were offered the same GLH strip-tilled with the AFC, including oat (Avena sativa), millet (Pennisetum glaucum), teff (Eragrostis tef), buckwheat (Fagopyrum esculentum), and chickling vetch (Lathyrus sativus). Milk samples were collected from each cow during a.m. and p.m. milkings on d 19 to 21 of grazing, and composite milk samples per cow were analyzed for (1) the high-abundance milk protein profile, (2) the skim milk low-abundance protein-enriched proteome, and (3) the MFGM proteome. Of the 443 proteins identified in the skim and MFGM proteomes, 433 were included in statistical analysis, including 68 proteins identified in the skim milk fraction and 365 in the MFGM-associated fraction. Analysis of the skim and MFGM proteomes encompassed unique gene ontology profiles and proportions of functional classifications. In response to diet, α_S1-casein as well as 8 low-abundance proteins were present in higher concentration or abundance in milk from cows grazing the GLH strip-tilled with the AFC compared with milk from cows grazing GLH, suggesting that even short-term grazing of pastures including some AFC may affect the milk proteome.

Identification of key pathways and genes in different types of chronic kidney disease based on WGCNA

Chronic kidney disease (CKD) is a highly heterogeneous nephrosis that occurs when the structure and function of the kidney is damaged. Gene expression studies have been widely used to elucidate various biological processes; however, the gene expression profile of CKD is currently unclear. The present study aimed to identify diagnostic biomarkers and therapeutic targets using renal biopsy sample data from patients with CKD. Gene expression data from 30 patients with CKD and 21 living donors were analyzed by weighted gene co-expression network analysis (WGCNA), in order to identify gene networks and profiles for CKD, as well as its specific characteristics, and to potentially uncover diagnostic biomarkers and therapeutic targets for patients with CKD. In addition, functional enrichment analysis was performed on co-expressed genes to determine modules of interest. Four co-expression modules were constructed from the WGCNA. The number of genes in the constructed modules ranged from 269 genes in the Turquoise module to 60 genes in the Yellow module. All four co-expression modules were correlated with CKD clinical traits (P<0.05). For example, the Turquoise module, which mostly contained genes that were upregulated in CKD, was positively correlated with CKD clinical traits, whereas the Blue, Brown and Yellow modules were negatively correlated with clinical traits. Functional enrichment analysis revealed that the Turquoise module was mainly enriched in genes associated with the ‘defense response’, ‘mitotic cell cycle’ and ‘collagen catabolic process’ Gene Ontology (GO) terms, implying that genes involved in cell cycle arrest and fibrogenesis were upregulated in CKD. Conversely, the Yellow module was mainly enriched in genes associated with ‘glomerulus development’ and ‘kidney development’ GO terms, indicating that genes associated with renal development and damage repair were downregulated in CKD. The hub genes in the modules were acetyl-CoA carboxylase α, cyclin-dependent kinase 1, Wilm's tumour 1, NPHS2 stomatin family member, podocin, JunB proto-oncogene, AP-1 transcription factor subunit, activating transcription factor 3, forkhead box O1 and v-abl Abelson murine leukemia viral oncogene homolog 1, which were confirmed to be significantly differentially expressed in CKD biopsies. Combining the eight hub genes enabled a high capacity for discrimination between patients with CKD and healthy subjects, with an area under the receiver operating characteristic curve of 1.00. In conclusion, this study provided a framework for co-expression modules of renal biopsy samples from patients with CKD and living donors, and identified several potential diagnostic biomarkers and therapeutic targets for CKD.

Single Nucleotide Polymorphisms of NUCB2 and their Genetic Associations with Milk Production Traits in Dairy Cows

We previously used the RNA sequencing technique to detect the hepatic transcriptome of Chinese Holstein cows among the dry period, early lactation, and peak of lactation, and implied that the nucleobindin 2 (NUCB2) gene might be associated with milk production traits due to its expression being significantly increased in early lactation or peak of lactation as compared to dry period (q value < 0.05). Hence, in this study, we detected the single nucleotide polymorphisms (SNPs) of NUCB2 and analyzed their genetic associations with milk yield, fat yield, fat percentage, protein yield, and protein percentage. We re-sequenced the entire coding and 2000 bp of 5' and 3' flanking regions of NUCB2 by pooled sequencing, and identified ten SNPs, including one in 5' flanking region, two in 3' untranslated region (UTR), and seven in 3' flanking region. The single-SNP association analysis results showed that the ten SNPs were significantly associated with milk yield, fat yield, fat percentage, protein yield, or protein percentage in the first or second lactation (p values <= 1 × 10-4 and 0.05). In addition, we estimated the linkage disequilibrium (LD) of the ten SNPs by Haploview 4.2, and found that the SNPs were highly linked in one haplotype block (D' = 0.98-1.00), and the block was also significantly associated with at least one milk traits in the two lactations (p values: 0.0002-0.047). Further, we predicted the changes of transcription factor binding sites (TFBSs) that are caused by the SNPs in the 5' flanking region of NUCB2, and considered that g.35735477C>T might affect the expression of NUCB2 by changing the TFBSs for ETS transcription factor 3 (ELF3), caudal type homeobox 2 (CDX2), mammalian C-type LTR TATA box (VTATA), nuclear factor of activated T-cells (NFAT), and v-ets erythroblastosis virus E26 oncogene homolog (ERG) (matrix similarity threshold, MST > 0.85). However, the further study should be performed to verify the regulatory mechanisms of NUCB2 and its polymorphisms on milk traits. Our findings first revealed the genetic effects of NUCB2 on the milk traits in dairy cows, and suggested that the significant SNPs could be used in genomic selection to improve the accuracy of selection for dairy cattle breeding.

Identification and Expression Analysis of Long Noncoding RNAs in Fat-Tail of Sheep Breeds

Emerging evidence suggests that long non-coding RNAs (lncRNAs) participate in the regulation of a diverse range of biological processes. However, most studies have been focused on a few established model organisms and little is known about lncRNAs in fat-tail development in sheep. Here, the first profile of lncRNA in sheep fat-tail along with their possible roles in fat deposition were investigated, based on a comparative transcriptome analysis between fat-tailed (Lori-Bakhtiari) and thin-tailed (Zel) Iranian sheep breeds. Among all identified lncRNAs candidates, 358 and 66 transcripts were considered novel intergenic (lincRNAs) and novel intronic (ilncRNAs) corresponding to 302 and 58 gene loci, respectively. Our results indicated that a low percentage of the novel lncRNAs were conserved. Also, synteny analysis identified 168 novel lincRNAs with the same syntenic region in human, bovine and chicken. Only seven lncRNAs were identified as differentially expressed genes between fat and thin tailed breeds. Q-RT-PCR results were consistent with the RNA-Seq data and validated the findings. Target prediction analysis revealed that the novel lncRNAs may act in cis or trans and regulate the expression of genes that are involved in the lipid metabolism. A gene regulatory network including lncRNA-mRNA interactions were constructed and three significant modules were found, with genes relevant to lipid metabolism, insulin and calcium signaling pathway. Moreover, integrated analysis with AnimalQTLdb database further suggested six lincRNAs and one ilncRNAs as candidates of sheep fat-tail development. Our results highlighted the putative contributions of lncRNAs in regulating expression of genes associated with fat-tail development in sheep.

Genetic Effects of LPIN1 Polymorphisms on Milk Production Traits in Dairy Cattle

Our initial RNA sequencing work identified that lipin 1 (LPIN1) was differentially expressed during dry period, early lactation, and peak of lactation in dairy cows, and it was enriched into the fat metabolic Gene Ontology (GO) terms and pathways, thus we considered LPIN1 as the candidate gene for milk production traits. In this study, we detected the polymorphisms of LPIN1 and verified their genetic effects on milk yield and composition in a Chinese Holstein cow population. We found seven SNPs by re-sequencing the entire coding region and partial flanking region of LPIN1, including one in 5′ flanking region, four in exons, and two in 3′ flanking region. Of these, four SNPs, c.637T > C, c.708A > G, c.1521C > T, and c.1555A > C, in the exons were predicted to result in the amino acid replacements. With the Haploview 4.2, we found that seven SNPs in LPIN1 formed two haplotype blocks (D′ = 0.98–1.00). Single-SNP association analyses showed that SNPs were significantly associated with milk yield, fat yield, fat percentage, or protein yield in the first or second lactation (p = < 0.0001–0.0457), and only g.86049389C > T was strongly associated with protein percentage in both lactations (p = 0.0144 and 0.0237). The haplotype-based association analyses showed that the two haplotype blocks were significantly associated with milk yield, fat yield, protein yield, or protein percentage (p = < 0.0001–0.0383). By quantitative real-time PCR (qRT-PCR), we found that LPIN1 had relatively high expression in mammary gland and liver tissues. Furthermore, we predicted three SNPs, c.637T > C, c.708A > G, and c.1521C > T, using SOPMA software, changing the LPIN1 protein structure that might be potential functional mutations. In summary, we demonstrated the significant genetic effects of LPIN1 on milk production traits, and the identified SNPs could serve as genetic markers for dairy breeding.

Transcriptome comparison between prenatal and postnatal Large White livers identifies differences in the expression level of genes related to metabolism and postnatal growth

The current study examined the liver transcriptomic profiles of the Large White different in developmental periods. It was performed on pigs of two developmental stages: 70-day fetus (P70) and 70-day piglets (D70). The objective of the study was to identify genes associated with Large White liver lipid metabolism, growth and development. We sequenced eight sRNA libraries of liver samples from four Large White at P70 and D70 respectively. We totally obtained 19,202 genes. 4916 of them were found to be differentially expressed (DEGs) (p < 0.05, fold change ≥ 1), of which 2502 were up-regulated and 2414 were down-regulated. GO enrichment and KEGG pathway analysis indicated that ACACA, ACADM, ACAA2 and HADH were simultaneously enriched in diverse pathways related to lipid metabolism, and so they were considered to be the promising candidate genes which could affect the porcine liver lipid metabolism. Notably, the gene insulin-like growth factor 1 (IGF1) which participated in somatotropic axis signaling was found to be up-regulated in D70 compared with P70. miRWalk and TargetScan softwares were used to screen the miRNAs which bound to the 3' untranslated region (3'UTR) of IGF1. After integration analysis with miRNAs sequencing data, miR-18b and miR-130b-3p were selected for further study. MiR-18b and miR-130b-3p were down-regulated in D70 compared with P70. Dual luciferase assays indicated that miR-18b and miR-130b-3p could obviously decrease (p < 0.05) the fluorescence activity of the group transfected with the wild-type vector of IGF1 3'UTR, while the relative luciferase activity of the group transfected with the mutant vector of IGF1 3'UTR did not change significantly. Taken together, it indicated that miR-18b and miR-130b-3p could target IGF1 directly.

The Japanese Wagyu beef industry: current situation and future prospects - A review

In Japan, Wagyu cattle include four Japanese breeds; Black, Brown, Shorthorn, and Polled. Today, the renowned brand name Wagyu includes not only cattle produced in Japan, but also cattle produced in countries such as Australia and the United States. In recent years, the intramuscular fat percentage in beef (longissimus muscle) from Japanese Black cattle has increased to be greater than 30%. The Japanese Black breed is genetically predisposed to producing carcass lipids containing higher concentrations of monounsaturated fatty acids than other breeds. However, there are numerous problems with the management of this breed including high production costs, disposal of untreated excrement, the requirement for imported feed, and food security risks resulting from various viral diseases introduced by imported feed. The feeding system needs to shift to one that is more efficient, and improves management for farmers, food security for consumers, and the health environment for residents of Japan. Currently, we are developing a metabolic programming and an information and communications technology (ICT, or Interne of Things) management system for Wagyu beef production as future systems. If successful, we will produce safe, high-quality Wagyu beef using domestic pasture resources while solving the problems of how to utilize increasing areas of abandoned agricultural land and to make use of the plant-based feed resources in Japan’s mountainous areas.

Longissimus dorsi muscle transcriptomic analysis of Yunling and Chinese simmental cattle differing in intramuscular fat content and fatty acid composition

Intramuscular fat (IMF) content and fatty acid (FA) composition vary significantly across beef cattle breeds, which play an important role in taste and nutritional value. However, the molecular mechanisms underlying these phenotypic differences remain unknown. The present study compared meat quality traits between Yunling cattle and Chinese Simmental cattle. Yunling cattle showed a lower IMF content and proportion of monounsaturated fatty acids (MUFA), as well as higher proportions of saturated fatty acids (SFA), polyunsaturated fatty acids (PUFA), and short-chain fatty acids (sc-FA) in the longissimus dorsi (LD) muscle than Chinese Simmental cattle. To further identify the candidate genes and pathways responsible for these phenotypic differences, the transcriptome of LD muscle from the two breeds were measured using RNA-seq. A total of 1347 differentially expressed genes were identified. The major metabolic pathways that were differentially modulated were lipolysis and glycometabolism. Yunling cattle showed a higher expression of lipolysis genes (ALDH9A1, ACSL5, ACADM, ACAT2, ACOT2) and a lower expression of genes related to glycometabolism (PGM1, GALM, PGM1, GPI, LDHA). This research identified candidate genes and pathways for IMF content and FA composition in the LD muscle of beef cattle, which may facilitate the design of new selection strategies to improve meat quality.

Full-length sequencing and identification of novel polymorphisms in the ACACA gene of Valle del Belice sheep breed

The essential role of the acetyl-CoA carboxylase (ACACA) enzyme in milk fatty acid (FA) synthesis suggests that it may be responsible for the phenotypic variability observed inmilk.Before attempting association analyses between this gene and/or enzyme and phenotypic traits, a study on the genetic variability within this locus is required. The aim of this work was to sequence the entire coding region of ACACA gene in Valle del Belice sheep breed to identify polymorphic sites. A total of 51 coding exons of ACACA gene were sequenced in 32 individuals of Valle del Belice sheep breed. Sequencing analysis and alignment of obtained sequences showed the presence of 23 polymorphic sites. The most polymorphic was exon 53 which showed presence of 12 single-nucleotide polymorphisms (SNPs), ofwhich eightweremissensemutations, caused amino acid changes and therefore may affect protein function or stability causing variation in phenotype. The identified polymorphisms showed high variability of the ACACA gene. Sequences analysis allowed to find six new SNPs in exon 53 (6832C>T; 6835C>A; 6840G>A; 6847G>T; 6852C>T and 6860G>C). A total of 31 haplotypes were inferred. Although this study could not provide association study with production traits, it shows finding of novel SNPs that might be important in future studies and laid the basis for further association analyses needed to evaluate the potential use of these SNPs as genetic markers for fat content and FAs composition in milk of Valle del Belice sheep breed.

Genetic association of the ACACB gene with milk yield and composition traits in dairy cattle

Previously, we re-sequenced the whole genomes of eight Holstein bulls with high or low milk protein and fat percentage, and we detected two indels in the ACACB (acetyl-CoA carboxylase beta) gene that were polymorphic between the two groups. Thus, we considered ACACB as a promising candidate gene potentially affecting milk composition traits. Herein, we verified the genetic effects of ACACB on five milk traits in a Chinese Holstein population. We identified six SNPs in the 5'-promoter region, five in the 5'- untranslated region (UTR), 11 in exons, four in the 3'-UTR and three in the 3'-flanking region by re-sequencing the entire coding and regulatory regions of ACACB. One of these SNPs (ss1987461005) is reported here for the first time, and three of the SNPs (rs109482081, rs110819816 and rs109281947) were predicted to result in amino acid replacements. Genotype-phenotype association analyses showed that all the identified SNPs, except for ss1987461005, rs208919019 and rs134447911, were significantly associated with milk yield, fat yield, fat percentage, protein yield or protein percentage (P < 0.0001 to 0.0484). Linkage disequilibrium analyses were conducted among the identified SNPs to confirm the genetic associations. Two SNPs-rs135874354 (g.66218726T>C) and rs210928430 (g.66218117G>A)-were predicted to alter transcription factor binding sites in the 5'-promoter region of ACACB. A luciferase activity assay showed that the promoter activity of haplotype TG was significantly higher than that of CG (P = 0.0002) and that the promoter activity of haplotype TA was remarkably higher than that of CA (P = 7.4285E-09), showing that the T allele of rs135874354 increased promoter activity. Thus, rs135874354 was considered to be a potentially functional mutation. Our findings have, for the first time, profiled the genetic effect of ACACB on milk production traits in dairy cattle and revealed a potentially causal mutation that requires further the in-depth validation.

SNP co-association and network analyses identify E2F3, KDM5A and BACH2 as key regulators of the bovine milk fatty acid profile

The fatty acid (FA) profile has a considerable impact on the nutritional and technological quality of milk and dairy products. The molecular mechanism underlying the regulation of fat metabolism in bovine mammary gland have been not completely elucidated. We conducted genome-wide association studies (GWAS) across 65 milk FAs and fat percentage in 1,152 Brown Swiss cows. In total, we identified 175 significant single nucleotide polymorphism (SNPs) spanning all chromosomes. Pathway analyses revealed that 12:0 was associated with the greatest number of overrepresented categories/pathways (e.g. mitogen-activated protein kinase (MAPK) activity and protein phosphorylation), suggesting that it might play an important biological role in controlling milk fat composition. An Associated Weight Matrix approach based on SNP co-associations predicted a network of 791 genes related to the milk FA profile, which were involved in several connected molecular pathways (e.g., MAPK, lipid metabolism and hormone signalling) and undetectable through standard GWAS. Analysis of transcription factors and their putative target genes within the network identified BACH2, E2F3 and KDM5A as key regulators of milk FA metabolism. These findings contribute to increasing knowledge of FA metabolism and mammary gland functionality in dairy cows and may be useful in developing targeted breeding practices to improve milk quality.

The response of gene expression associated with lipid metabolism, fat deposition and fatty acid profile in the longissimus dorsi muscle of Gannan yaks to different energy levels of diets

The energy available from the diet, which affects fat deposition in vivo, is a major factor in the expression of genes regulating fat deposition in the longissimus dorsi muscle. Providing high-energy diets to yaks might increase intramuscular fat deposition and fatty acid concentrations under a traditional grazing system in cold seasons. A total of fifteen adult castrated male yaks with an initial body weight 274.3 ± 3.14 kg were analyzed for intramuscular adipose deposition and fatty acid composition. The animals were divided into three groups and fed low-energy (LE: 5.5 MJ/kg), medium-energy (ME: 6.2 MJ/kg) and high-energy (HE: 6.9 MJ/kg) diets, respectively. All animals were fed ad libitum twice daily at 08:00–09:00 am and 17:00–18:00 pm and with free access to water for 74 days, including a 14-d period to adapt to the diets and the environment. Intramuscular fat (IMF) content, fatty acid profile and mRNA levels of genes involved in fatty acid synthesis were determined. The energy levels of the diets significantly (P<0.05) affected the content of IMF, total SFA, total MUFA and total PUFA. C16:0, C18:0 and C18:1n9c account for a large proportion of total fatty acids. Relative expression of acetyl-CoA carboxylase (ACACA), fatty acid synthase (FASN), stearoyl-CoA desaturase (SCD), sterol regulatory element-binding protein-1c (SREBP-1c), peroxisome proliferator-activated receptor γ (PPARγ) and fatty acid-binding protein 4 (FABP4) was greater in HE than in LE yaks (P<0.05). Moreover, ME yaks had higher (P<0.05) mRNA expression levels of PPARγ, ACACA, FASN, SCD and FABP4 than did the LE yaks. The results demonstrate that the higher energy level of the diets increased IMF deposition and fatty acid content as well as increased intramuscular lipogenic gene expression during the experimental period.

Functional validation of GPIHBP1 and identification of a functional mutation in GPIHBP1 for milk fat traits in dairy cattle

In a previous genome-wide association study (GWAS) on milk production traits in a Chinese Holstein population, we revealed that GPIHBP1 is a novel promising candidate gene for milk fat content traits. In this study, we performed over-expression and RNAi experiments on GPIHBP1 in bovine primary mammary epithelial cells. The results showed that the expression of several important milk fat-related genes (LPL, CD36, VLDLR, ACACA and FASN) increased or decreased when the expression of GPIHBP1 was up- or down-regulated. To identify the potential functional SNP involved, we explored the genetic variants of GPIHBP1 and found that a G/A mutation (chr14:2553998) in the promoter region of GPIHBP1 significantly reduced promoter activity and had an effect on transcription factor binding sites. This finding was consistent with the lower expression of GPIHBP1 observed in the mammary gland tissue of cows harboring the homozygous AA mutation compared with wild-type homozygous GG or heterozygous AG. Furthermore, association analysis showed that cows with the AA genotype outperformed those with the GG and AG genotypes in terms of the milk fat percentage. Our study demonstrates that GPIHBP1 could be a strong candidate gene for milk fat content traits and, in particular, the G to A mutation at chr14:2553998 within GPIHBP1 could be a functional mutation related to its effects.

Differential expression of genes in milk of dairy cattle during lactation

The milk fat globule (MFG) is one of the most representative of mammary gland tissues and can be utilized to study gene expression of lactating cows during lactation. In this study, RNA-seq technology was employed to detect differential expression of genes in MFGs at day 10 and day 70 after calving between two groups of cows with extremely high (H group) and low (L group) 305-day milk yield, milk fat yield and milk protein yield. In total, 1232, 81, 429 and 178 significantly differentially expressed genes (false discovery rate q < 0.05) were detected between H10 and L10, H70 and L70, H10 and H70, and L10 and L70 respectively. Gene Ontology enrichment and pathway analysis revealed that these differentially expressed genes were enriched in biological processes involved in mammary gland development, protein and lipid metabolism process, signal transduction, cellular process, differentiation and immune function. Among these differentially expressed genes, 178 (H10 vs. L10), 4 (H70 vs. L70), 68 (H10 vs. H70) and 22 (L10 vs. L70) were found to be located within previously reported QTL regions for milk production traits. Based on these results, some promising candidate genes for milk production traits in dairy cattle were suggested.