Polymorphisms of CYP7A1 and HADHB Genes and Their Effects on Milk Production Traits in Chinese Holstein Cows

Our preliminary research proposed the cytochrome P450 family 7 subfamily A member 1 (CYP7A1) and hydroxyacyl-coenzyme A dehydrogenase trifunctional multienzyme complex beta subunit (HADHB) genes as candidates for association with milk-production traits in dairy cattle because of their differential expression across different lactation stages in the liver tissues of Chinese Holstein cows and their potential roles in lipid metabolism. Hence, we identified single-nucleotide polymorphisms (SNPs) of the CYP7A1 and HADHB genes and validated their genetic effects on milk-production traits in a Chinese Holstein population with the goal of providing valuable genetic markers for genomic selection (GS) in dairy cattle, This study identified five SNPs, 14:g.24676921A>G, 14:g.24676224G>A, 14:g.24675708G>T, 14:g.24665961C>T, and 14:g.24664026A>G, in the CYP7A1 gene and three SNPs, 11:g.73256269T>C, 11:g.73256227A>C, and 11:g.73242290C>T, in HADHB. The single-SNP association analysis revealed significant associations (p value ≤ 0.0461) between the eight SNPs of CYP7A1 and HADHB genes and 305-day milk, fat and protein yields. Additionally, using Haploview 4.2, we found that the five SNPs of CYP7A1 formed two haplotype blocks and that the two SNPs of HADHB formed one haplotype block; notably, all three haplotype blocks were also significantly associated with milk, fat and protein yields (p value ≤ 0.0315). Further prediction of transcription factor binding sites (TFBSs) based on Jaspar software (version 2023) showed that the 14:g.24676921A>G, 14:g.24675708G>T, 11:g.73256269T>C, and 11:g.73256227A>C SNPs could alter the 5' terminal TFBS of the CYP7A1 and HADHB genes. The 14:g.24665961C>T SNP caused changes in the structural stability of the mRNA for the CYP7A1 gene. These alterations have the potential to influence gene expression and, consequently, the phenotype associated with milk-production traits. In summary, we have confirmed the genetic effects of CYP7A1 and HADHB genes on milk-production traits in dairy cattle and identified potential functional mutations that we suggest could be used for GS of dairy cattle and in-depth mechanistic studies of animals.

Genetic Background of Blood β-Hydroxybutyrate Acid Concentrations in Early-Lactating Holstein Dairy Cows Based on Genome-Wide Association Analyses

Ketosis is a common metabolic disorder in the early lactation of dairy cows. It is typically diagnosed by measuring the concentration of β-hydroxybutyrate (BHB) in the blood. This study aimed to estimate the genetic parameters of blood BHB and conducted a genome-wide association study (GWAS) based on the estimated breeding value. Phenotypic data were collected from December 2019 to August 2023, comprising blood BHB concentrations in 45,617 Holstein cows during the three weeks post-calving across seven dairy farms. Genotypic data were obtained using the Neogen Geneseek Genomic Profiler (GGP) Bovine 100 K SNP Chip and GGP Bovine SNP50 v3 (Illumina Inc., San Diego, CA, USA) for genotyping. The estimated heritability and repeatability values for blood BHB levels were 0.167 and 0.175, respectively. The GWAS result detected a total of ten genome-wide significant associations with blood BHB. Significant SNPs were distributed in Bos taurus autosomes (BTA) 2, 6, 9, 11, 13, and 23, with 48 annotated candidate genes. These potential genes included those associated with insulin regulation, such as INSIG2, and those linked to fatty acid metabolism, such as HADHB, HADHA, and PANK2. Enrichment analysis of the candidate genes for blood BHB revealed the molecular functions and biological processes involved in fatty acid and lipid metabolism in dairy cattle. The identification of novel genomic regions in this study contributes to the characterization of key genes and pathways that elucidate susceptibility to ketosis in dairy cattle.

Understanding the underlying genetic mechanisms for age at first calving, inter-calving period and scrotal circumference in Bonsmara cattle

BACKGROUND

Reproduction is a key feature of the sustainability of a species and thus represents an important component in livestock genetic improvement programs. Most reproductive traits are lowly heritable. In order to gain a better understanding of the underlying genetic basis of these traits, a genome-wide association was conducted for age at first calving (AFC), first inter-calving period (ICP) and scrotal circumference (SC) within the South African Bonsmara breed. Phenotypes and genotypes (120,692 single nucleotide polymorphisms (SNPs) post editing) were available on 7,128 South African Bonsmara cattle; the association analyses were undertaken using linear mixed models.

RESULTS

Genomic restricted maximum likelihood analysis of the 7,128 SA Bonsmara cattle yielded genomic heritability's of 0.183 (SE = 0.021) for AFC, 0.207 (SE = 0.022) for ICP and 0.209 (SE = 0.019) for SC. A total of 16, 23 and 51 suggestive (P ≤ 4 × 10-6) SNPs were associated with AFC, ICP and SC, while 11, 11 and 44 significant (P ≤ 4 × 10-7) SNPs were associated with AFC, ICP and SC respectively. A total of 11 quantitative trait loci (QTL) and 11 candidate genes were co-located with these associated SNPs for AFC, with 10 QTL harbouring 11 candidate genes for ICP and 41 QTL containing 40 candidate genes for SC. The QTL identified were close to genes previously associated with carcass, fertility, growth and milk-related traits. The biological pathways influenced by these genes include carbohydrate catabolic processes, cellular development, iron homeostasis, lipid metabolism and storage, immune response, ovarian follicle development and the regulation of DNA transcription and RNA translation.

CONCLUSIONS

This was the first attempt to study the underlying polymorphisms associated with reproduction in South African beef cattle. Genes previously reported in cattle breeds for numerous traits bar AFC, ICP or SC were detected in this study. Over 20 different genes have not been previously reported in beef cattle populations and may have been associated due to the unique genetic composite background of the SA Bonsmara breed.

Genetic marker identification of SEC13 gene for milk production traits in Chinese holstein

SEC13 homolog, nuclear pore and COPII coat complex component (SEC13) is the core component of the cytoplasmic COPII complex, which mediates material transport from the endoplasmic reticulum to the Golgi complex. Our preliminary work found that SEC13 gene was differentially expressed in dairy cows during different stages of lactation, and involved in metabolic pathways of milk synthesis such as citric acid cycle, fatty acid, starch and sucrose metabolisms, so we considered that the SEC13 might be a candidate gene affecting milk production traits. In this study, we detected the polymorphisms of SEC13 gene and verified their genetic effects on milk yield and composition traits in a Chinese Holstein cow population. By sequencing the whole coding and partial flanking regions of SEC13, we found four single nucleotide polymorphisms (SNPs). Subsequent association analysis showed that these four SNPs were significantly associated with milk yield, fat yield, protein yield or protein percentage in the first and second lactations (p ≤.0351). We also found that two SNPs in SEC13 formed one haplotype block by Haploview4.2, and the block was significantly associated with milk yield, fat yield, fat percentage, protein yield or protein percentage (p ≤ .0373). In addition, we predicted the effect of SNP on 5'region on transcription factor binding sites (TFBSs), and found that the allele A of 22:g.54362761A>G could bind transcription factors (TFs) GATA5, GATA3, HOXD9, HOXA10, CDX1 and Hoxd13; and further dual-luciferase reporter assay verified that the allele A of this SNP inhibited the fluorescence activity. We speculate that the A allele of 22:g.54362761A>G might inhibit the transcriptional activity of SEC13 gene by binding the TFs, which may be a cause mutation affecting the formation of milk production traits in dairy cows. In summary, we proved that SEC13 has a significant genetic effect on milk production traits and the identified significant SNPs could be used as candidate genetic markers for GS SNP chips development; on the other hand, we verified the transcriptional regulation of 22:g.54362761A>G on SEC13 gene, providing research direction for further function validation tests.

Analysis of runs of homozygosity of cattle living in different climate zones

Aim of this study was to analyze the distribution and characteristics of runs of homozygosity in Bos taurus taurus and Bos taurus indicus breeds, as well as their crosses, farmed all around the world. With this aim in view, we used single-nucleotide polymorphisms (SNP) genotypes for 3,263 cattle belonging to 204 different breeds. After quality control, 23,311 SNPs were retained for the analysis. Animals were divided into seven different groups: 1) continental taurus, 2) temperate taurus, 3) temperate indicus, 4) temperate composite, 5) tropical taurus, 6) tropical indicus, and 7) tropical composite. The climatic zones were created according to the latitude of the breeds' country of origin: i) continental, latitude ≥ 45°; ii) temperate, 45°< Latitude >23.26°; iii) tropics, latitude ≤ 23.26°. Runs of homozygosity were computed as 15 SNPs spanning in at least 2 Mb; number of ROH per animal (nROH), average ROH length (meanMb), and ROH-based inbreeding coefficients (FROH) were also computed. Temperate indicus showed the largest nROH, whereas Temperate taurus the lowest value. Moreover, the largest meanMb was observed for Temperate taurus, whereas the lowest value for Tropics indicus. Temperate indicus breeds showed the largest FROH values. Genes mapped in the identified ROH were reported to be associated with the environmental adaptation, disease resistance, coat color determinism, and production traits. Results of the present study confirmed that runs of homozygosity could be used to identify genomic signatures due to both artificial and natural selection.

Assessing Population Structure and Signatures of Selection in Wanbei Pigs Using Whole Genome Resequencing Data

Wanbei pig (WBP) is one of the indigenous pig resources in China and has many germplasm characteristics. However, research on its genome is lacking. To assess the genomic variation, population structure, and selection signatures, we resequenced 18 WBP for the first time and performed a comprehensive analysis with resequenced data of 10 Asian wild boars. In total, 590.03 Gb of data and approximately 41 million variants were obtained. Polymorphism level (θπ) ratio and genetic differentiation (fixation index)-based cross approaches were applied, and 539 regions, which harbored 176 genes, were selected. Functional analysis of the selected genes revealed that they were associated with lipid metabolism (SCP2, APOA1, APOA4, APOC3, CD36, BCL6, ADCY8), backfat thickness (PLAG1, CACNA2D1), muscle (MYOG), and reproduction (CABS1). Overall, our results provide a valuable resource for characterizing the uniqueness of WBP and a basis for future breeding.

Bta-miR-106b Regulates Bovine Mammary Epithelial Cell Proliferation, Cell Cycle, and Milk Protein Synthesis by Targeting the CDKN1A Gene

Our previous studies found that bta-miR-106b and its corresponding target gene, CDKN1A, were differentially expressed between the mammary epithelium of lactating Holstein cows with extremely high and low milk protein and fat percentage, implying the potential role of bta-miR-106b in milk composition synthesis. In this study, with luciferase assay experiment, bta-miR-106b was validated to target the 3'-untranslated region (UTR) of bovine CDKN1A, thereby regulating its expression. Moreover, in bovine mammary epithelial cells (BMECs), over-expression of bta-miR-106b significantly down-regulated the CDKN1A expression at both mRNA and protein levels, and inhibitors of bta-miR-106b increased CDKN1A expression. Of note, we observed that bta-miR-106b accelerated cell proliferation and cell cycle, and changed the expressions of protein synthesis related pathways such as JAK-STAT and PI3K/AKT/mTOR through regulating CDKN1A expression. Our findings highlight the important regulatory role of bta-miR-106b in milk protein synthesis by targeting CDKN1A in dairy cattle.

Single Nucleotide Polymorphisms of ALDH18A1 and MAT2A Genes and Their Genetic Associations with Milk Production Traits of Chinese Holstein Cows

Our preliminary work had suggested two genes, aldehyde dehydrogenase 18 family member A1 (ALDH18A1) and methionine adenosyltransferase 2A (MAT2A), related to amino acid synthesis and metabolism as candidates affecting milk traits by analyzing the liver transcriptome and proteome of dairy cows at different lactation stages. In this study, the single nucleotide polymorphisms (SNPs) of ALDH18A1 and MAT2A genes were identified and their genetic effects and underlying causative mechanisms on milk production traits in dairy cattle were analyzed, with the aim of providing effective genetic information for the molecular breeding of dairy cows. By resequencing the entire coding and partial flanking regions of ALDH18A1 and MAT2A, we found eight SNPs located in ALDH18A1 and two in MAT2A. Single-SNP association analysis showed that most of the 10 SNPs of these two genes were significantly associated with the milk yield traits, 305-day milk yield, fat yield, and protein yield in the first and second lactations (corrected p ≤ 0.0488). Using Haploview 4.2, we found that the seven SNPs of ALDH18A1 formed two haplotype blocks; subsequently, the haplotype-based association analysis showed that both haplotypes were significantly associated with 305-day milk yield, fat yield, and protein yield (corrected p ≤ 0.014). Furthermore, by Jaspar and Genomatix software, we found that 26:g.17130318 C>A and 11:g.49472723G>C, respectively, in the 5′ flanking region of ALDH18A1 and MAT2A genes changed the transcription factor binding sites (TFBSs), which might regulate the expression of corresponding genes to affect the phenotypes of milk production traits. Therefore, these two SNPs were considered as potential functional mutations, but they also require further verification. In summary, ALDH18A1 and MAT2A were proved to probably have genetic effects on milk production traits, and their valuable SNPs might be used as candidate genetic markers for dairy cattle’s genomic selection (GS).

Whole Transcriptome Profiling of the Effects of Cadmium on the Liver of the Xiangxi Yellow Heifer

Cadmium (Cd) is a major heavy metal toxicant found in industrial zones. Humans and animals are exposed to it through their diet, which results in various physiological problems. In the current study, the toxic effects of Cd on the liver were investigated by whole-transcriptome sequencing (RNA-seq) of the livers of Xiangxi heifers fed a diet with excess Cd. We randomly divided six healthy heifers into two groups. The first group received a control diet, whereas the second group received Cd-exceeding diets for 100 days. After 100 days, the livers were collected. A total of 551 differentially expressed mRNAs, 24 differentially expressed miRNAs, and 169 differentially expressed lncRNAs were identified (p < 0.05, |log2FC| >1). Differentially expressed genes (DEGs) were analyzed by gene ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses. We found that under Cd exposure, DEGs were enriched in the adenosine 5'-monophosphate-activated protein kinase pathway, which is involved in autophagy regulation, and the peroxisome proliferator-activated receptor pathway, which is involved in lipid metabolism. In addition, the apolipoprotein A4 gene, which has anti-inflammatory and antioxidant effects, the anti-apoptotic gene ATPase H+/K+ transporting the nongastric alpha2 subunit, and the cholesterol metabolism-associated gene endothelial lipase gene were significantly downregulated. C-X-C motif chemokine ligand 3, cholesterol 7α-hydroxylase, and stearoyl-CoA desaturase, which are involved in the development of fatty liver, were significantly upregulated. These genes revealed the main effects of Cd on the liver of Xiangxi yellow heifers. The current study provides insightful information regarding the DEGs involved in autophagy regulation, apoptosis, lipid metabolism, anti-inflammation, and antioxidant enzyme activity. These may serve as useful biomarkers for predicting and treating Cd-related diseases in the future.

Longitudinal liver proteome profiling in dairy cows during the transition from gestation to lactation: Investigating metabolic adaptations and their interactions with fatty acids supplementation via repeated measurements ANOVA-simultaneous component analysis

Repeated measurements analysis of variance - simultaneous component analysis (ASCA) has been developed to handle complex longitudinal omics datasets and combine novel information with existing data. Herein, we aimed at applying ASCA to 64 liver proteomes collected at 4-time points (day -21, +1, +28, and + 63 relative to parturition) from 16 Holstein cows treated from 9 wk. antepartum to 9 wk. postpartum (PP) with coconut oil (CTRL) or a mixture of essential fatty acids (EFA) and conjugated linoleic acid (CLA) (EFA + CLA). The ASCA modeled 116, 43, and 97 differentially abundant proteins (DAP) during the transition to lactation, between CTRL and EFA + CLA, and their interaction, respectively. Time-dependent DAP were annotated to pathways related to the metabolism of carbohydrates, FA, and amino acid in the PP period. The DAP between FA and the interaction effect were annotated to the metabolism of xenobiotics by cytochrome P450, drug metabolism - cytochrome P450, retinol metabolism, and steroid hormone biosynthesis. Collectively, ASCA provided novel information on molecular markers of metabolic adaptations and their interactions with EFA + CLA supplementation. Bioinformatics analysis suggested that supplemental EFA + CLA amplified hepatic FA oxidation; cytochrome P450 was enriched to maintain metabolic homeostasis by oxidation/detoxification of endogenous compounds and xenobiotics. SIGNIFICANCE: This report is among the first ones applying repeated measurement analysis of variance-simultaneous component analysis (ASCA) to deal with longitudinal proteomics results. ASCA separately identified differentially abundant proteins (DAP) in 'transition time', 'between fatty acid treatments', and 'their interaction'. We first identified the molecular signature of hepatic metabolic adaptations during postpartum negative energy balance; the enriched pathways were well-known pathways related to mobilizing fatty acids (FA) and amino acids to support continuous energy production through fatty acid oxidation, TCA cycle, and gluconeogenesis. Some of the DAP were not previously reported in transition dairy cows. Secondly, we provide novel information on the mechanisms by which supplemented essential FA and conjugated linoleic acids interact with hepatic metabolism. In this regard, FA amplified hepatic detoxifying and oxidation capacity through ligand activation of nuclear receptors. Finally, we briefly compared the strengths and weaknesses of the ASCA model with PLS-DA and outlined why these methods are complementary.

Isobaric Tags for Relative and Absolute Quantification-Based Proteomics Reveals Candidate Proteins of Fat Deposition in Chinese Indigenous Sheep With Morphologically Different Tails

Background: Chinese indigenous sheep can be classified into two types according to their tail morphology: fat-rumped and thin-tailed sheep, of which the typical breeds are Altay sheep and Tibetan sheep, respectively.

Methods: To identify the differentially expressed proteins (DEPs) underlying the phenotypic differences between tail types, we used isobaric tags for relative and absolute quantification (iTRAQ) combined with multi-dimensional liquid chromatography tandem-mass spectrometry (LC-MS/MS) technology to detect candidate proteins. We then subjected these to a database search and identified the DEPs. Finally, bioinformatics technology was used to carry out Gene Ontology (GO) functional and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses.

Results: A total of 3,248 proteins were identified, of which 44 were up-regulated and 40 were down-regulated DEPs. Analyzing their GO function terms and KEGG pathways revealed that the functions of these DEPs are mainly binding, catalytic activity, structural molecule activity, molecular function regulator, and transporter activity. Among the genes encoding the DEPs, APOA2, GALK1, ADIPOQ, and NDUFS4 are associated with fat formation and metabolism.

Conclusion: The APOA2, GALK1, ADIPOQ, and NDUFS4 genes may be involved in the deposition of fat in the tail of sheep. This study provides a scientific basis for the breeding of thin-tailed sheep.

When Two plus Two Is More than Four: Evidence for a Synergistic Effect of Fatty Acids on Peroxisome Proliferator-Activated Receptor Activity in a Bovine Hepatic Model

The inclusion of fat in livestock diets represents a valuable and cost-effective way to increase the animal’s caloric intake. Beyond their caloric value, fatty acids can be understood in terms of their bioactivity, via the modulation of the ligand-dependent nuclear peroxisome proliferator-activated receptors (PPAR). Isotypes of PPAR regulate important metabolic processes in both monogastric and ruminant animals, including the metabolism of fatty acids (FA), the production of milk fat, and the immune response; however, information on the modulation of bovine PPAR by fatty acids is limited. The objective of this study was to expand our understanding on modulation of bovine PPAR by FA, both when used individually and in combination, in an immortalized cell culture model of bovine liver. Of the 10 FA included in the study, the greatest activation of the PPAR reporter was detected with saturated FA C12:0, C16:0, and C18:0, as well as phytanic acid, and the unsaturated FA C16:1 and C18:1. When supplemented in mixtures of 2 FA, the most effective combination was C12:0 + C16:0, while in mixtures of 3 FA, the greatest activation was caused by combinations of C12:0 with C16:0 and either C18:0, C16:1, or C18:1. Some mixtures display a synergistic effect that leads to PPAR activation greater than the sum of their parts, which may be explained by structural dynamics within the PPAR ligand-binding pocket. Our results provide fundamental information for the development of tailored dietary plans that focus on the use of FA mixtures for nutrigenomic purposes.

Identification and analysis of lipid droplet-related proteome in the adipose tissue of grass carp (Ctenopharyngodon idella) under fed and starved conditions

Fat accumulation in the mesenteric adipose tissue is a serious problem in grass carp (Ctenopharyngodon idella) culture. Lipid droplet-related proteins (LDRPs) are involved in the formation, degradation, and biological functions of lipid droplets. In this study, we aimed to provide reference proteomics data to study lipid droplet regulation in fish. We isolated LDRPs from the mesenteric adipose tissue of grass carp (1-year-old) after normal feeding and 7 days of starvation, and identified and analysed them using isobaric tags for relative and absolute quantitation (iTRAQ) technology. Short-term starvation had no significant effect on the body weight, condition factor, visceral index, hepatopancreas index, intraperitoneal fat index, adipose tissue triglyceride content, and adipocyte size of grass carp. Nine hundred and fifty proteins were identified and annotated using the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases; they are involved in a variety of metabolic and signalling pathways, including amino acid, lipid, and carbohydrate metabolism, and the PI3K-Akt signalling pathway. There were 296 differentially expressed proteins (DEPs), with 143 up-regulated and 153 down-regulated proteins. Three proteins involved in triglyceride and fatty acid syntheses and two proteins involved in autophagy were up-regulated, and six proteins involved in lipid catabolism were down-regulated. These results indicate that under short-term starvation, lipid droplets in the adipose tissue of grass carp may maintain their shape by promoting fat production and inhibiting lipolysis, and autophagy may be one of the main strategies for coping with short-term energy deprivation.

MiR-16a Regulates Milk Fat Metabolism by Targeting Large Tumor Suppressor Kinase 1 (LATS1) in Bovine Mammary Epithelial Cells

Milk contains a number of beneficial fatty acids including short and medium chain and unsaturated conjugated and nonconjugated fatty acids. In this study, microRNA sequencing of mammary tissue collected in early-, peak-, mid-, and late-lactation periods was performed to determine the miRNA expression profiles. miR-16a was one of the differentially expressed miRNA and was selected for in-depth functional studies pertaining to fatty acid metabolism. The mimic of miR-16a impaired fat metabolism [triacylglycerol (TAG) and cholesterol] while knock-down of miR-16a promoted fat metabolism in vitro in bovine mammary epithelial cells (BMECs). In addition, the in vitro work with BMECs also revealed that miR-16a had a negative effect on the cellular concentration of cis 9-C18:1, total C18:1, C20:1, and C22:1 and long-chain polyunsaturated fatty acids. Therefore, these data suggesting a negative effect on fatty acid metabolism extend the discovery of the key role of miR-16a in mediating adipocyte differentiation. Through a combination of bioinformatics analysis, target gene 3' UTR luciferase reporter assays, and western blotting, we identified large tumor suppressor kinase 1 (LATS1) as a target of miR-16a. Transfection of siRNA-LATS1 into BMECs led to increases in TAG, cholesterol, and cellular fatty acid concentrations, suggesting a positive role of LATS1 in mammary cell fatty acid metabolism. In summary, data suggest that miR-16a regulates biological processes associated with intracellular TAG, cholesterol, and unsaturated fatty acid synthesis through LATS1. These data provide a theoretical and experimental framework for further clarifying the regulation of lipid metabolism in mammary cells of dairy cows.