Polymorphisms and genetic effects of PRLR, MOGAT1, MINPP1 and CHUK genes on milk fatty acid traits in Chinese Holstein

Signatures of positive selection after the introduction of genomic selection in the Finnish Ayrshire population

The Finnish Ayrshire (FAY) belongs to the Nordic Red breeds and is characterized by high milk yield, high milk components, good fertility, and functional conformation. The FAY breeding program is based on genomic selection. Despite the benefits of selection on breeding values, autozygosity in the genome may increase due to selection, and increased autozygosity may cause inbreeding depression in selected traits. However, there is lack of studies concerning selection signatures in the FAY after genomic selection introduction. The aim of this study was to identify signatures of selection in FAY after the introduction of genomic selection. Genomic data included 45,834 SNPs. The genotyped animals were divided into 2 groups: animals born before genomic selection introduction (6,108 cows) and animals born after genomic selection introduction (47,361 cows). We identified the selection signatures using 3 complementary methods: 2 based on identification of selection signatures from runs of homozygosity (ROH) islands and one based on the decay of site-specific extended haplotype between populations at SNP sites (Rsb). In total, we identified 34 ROH islands on chromosomes 1, 3, 6, 8, 12-15, 17, 19, 22, and 26 in FAY animals born before genomic selection (between 1980 and 2011) and 30 ROH islands on chromosomes 1-3, 13-17, 22, and 25-26 in FAY animals born after genomic selection introduction (between 2015 and 2020). We additionally detected 22 ΔROH islands on chromosomes 2-3, 11, 13, 14, 16, 18, 20, and 25-26. Finally, a total of 31 Rsb regions on chromosomes 2, 3, 14, 18, 20, and 25 were identified. Based on the results, genomic selection has favored certain alleles and haplotypes on genomic regions related to traits relevant in the FAY breeding program: milk production, fertility, growth, beef production traits, and feed efficiency. Several genes related to these traits (e.g., PLA2G4A, MECR, CHUK, COX15, RICTOR, SHISA9, and SEMA4G) overlapped or partially overlapped the observed selection signature regions. The association of genotypes within these regions and their effects on traits relevant in the FAY breeding program should be studied and genetic regions undergoing selection monitored in the FAY population.

Abnormal expression of serum miR-4746-5p in liver cancer patients after interventional chemotherapy and its possible mechanism

Interventional chemotherapy is a common operation in the clinical treatment of liver cancer. The aim of this study was to investigate the expression and molecular mechanism of serum miR-4746-5p in liver cancer patients before and after interventional chemotherapy. The levels of miR-4746-5p and CDKN1C in serum samples from liver cancer patients were detected using real-time fluorescence quantitative polymerase chain reaction. Receiver operating characteristic curves revealed the diagnostic value of miR-4746-5p in tumors. Differences in clinical indicators between liver cancer patients and healthy controls were assessed using Pearson correlation analysis. Luciferase reporter gene assays confirmed the targeted interaction between miR-4746-5p and CDKN1C. In vitro cellular assays were validated by Cell Counting Kit-8, Transwell assay, and chemoresistance assay. Serum miR-4746-5p levels were increased in liver cancer patients but were downregulated after chemotherapy intervention. CDKN1C expression showed the opposite trend. Low levels of miR-4746-5p mediated cell growth and metastasis by targeting and negatively regulating CDKN1C expression, while silencing CDKN1C restored cell activity. Inhibition of miR-4746-5p reduced chemoresistance, while downregulation of CDKN1C affected cell sensitivity. miR-4746-5p may be a potential therapeutic factor for liver cancer diagnosis and interventional chemotherapy.

Single-step genome-wide association for selected milk fatty acids in Dual-Purpose Belgian Blue cows

The aim of this study was to estimate genetic parameters and identify genomic regions associated with selected individual and groups of milk fatty acids (FA) predicted by milk mid-infrared spectrometry in Dual-Purpose Belgian Blue cows. The used data were 69,349 test-day records of milk yield, fat percentage, and protein percentage along with selected individual and groups FA of milk (g/dL milk) collected from 2007 to 2020 on 7,392 first-parity (40,903 test-day records), and 5,185 second-parity (28,446 test-day records) cows distributed in 104 herds in the Walloon Region of Belgium. Data of 28,466 SNPs, located on 29 Bos taurus autosomes (BTA), of 1,699 animals (639 males and 1,060 females) were used. Random regression test-day models were used to estimate genetic parameters through the Bayesian Gibbs sampling method. The SNP solutions were estimated using a single-step genomic best linear unbiased prediction approach. The proportion of genetic variance explained by each 25-SNP sliding window (with an average size of ~2 Mb) was calculated, and regions accounting for at least 1.0% of the total additive genetic variance were used to search for candidate genes. Average daily heritability estimated for the included milk FA traits ranged from 0.01 (C4:0) to 0.48 (C12:0) and 0.01 (C4:0) to 0.42 (C12:0) in the first and second parities, respectively. Genetic correlations found between milk yield and the studied individual milk FA, except for C18:0, C18:1 trans, C18:1 cis-9, were positive. The results showed that fat percentage and protein percentage were positively genetically correlated with all studied individual milk FA. Genome-wide association analyses identified 11 genomic regions distributed over 8 chromosomes [BTA1, BTA4, BTA10, BTA14 (4 regions), BTA19, BTA22, BTA24, and BTA26] associated with the studied FA traits, though those found on BTA14 partly overlapped. The genomic regions identified differed between parities and lactation stages. Although these differences in genomic regions detected may be due to the power of quantitative trait locus detection, it also suggests that candidate genes underlie the phenotypic expression of the studied traits may vary between parities and lactation stages. These findings increase our understanding about the genetic background of milk FA and can be used for the future implementation of genomic evaluation to improve milk FA profile in Dual-Purpose Belgian Blue cows.

Genetic polymorphisms of PKLR gene and their associations with milk production traits in Chinese Holstein cows

Our previous work had confirmed that pyruvate kinase L/R (PKLR) gene was expressed differently in different lactation periods of dairy cattle, and participated in lipid metabolism through insulin, PI3K-Akt, MAPK, AMPK, mTOR, and PPAR signaling pathways, suggesting that PKLR is a candidate gene to affect milk production traits in dairy cattle. Here, we verified whether this gene has significant genetic association with milk yield and composition traits in a Chinese Holstein cow population. In total, we identified 21 single nucleotide polymorphisms (SNPs) by resequencing the entire coding region and partial flanking region of PKLR gene, in which, two SNPs were located in 5′ promoter region, two in 5′ untranslated region (UTR), three in introns, five in exons, six in 3′ UTR and three in 3′ flanking region. The single marker association analysis displayed that all SNPs were significantly associated with milk yield, fat and protein yields or protein percentage (p ≤ 0.0497). The haplotype block containing all the SNPs, predicted by Haploview, had a significant association with fat yield and protein percentage (p ≤ 0.0145). Further, four SNPs in 5′ regulatory region and eight SNPs in UTR and exon regions were predicted to change the transcription factor binding sites (TFBSs) and mRNA secondary structure, respectively, thus affecting the expression of PKLR, leading to changes in milk production phenotypes, suggesting that these SNPs might be the potential functional mutations for milk production traits in dairy cattle. In conclusion, we demonstrated that PKLR had significant genetic effects on milk production traits, and the SNPs with significant genetic effects could be used as candidate genetic markers for genomic selection (GS) in dairy cattle.

Assessment of the association of the MOGAT1 and MOGAT3 gene with growth traits in different growth stages in Holstein calves

The members of the monoacylglycerol acyltransferase (MOGAT) family are essential candidate genes that influence economic traits associated with triglyceride synthesis, dietary fat absorption, and storage in livestock. In addition, the MOGAT gene family may also play an essential function in human polygenic diseases, like type 2 diabetes and obesity. The present study was conducted on Holstein calves to find the association between MOGAT1, MOGAT3/g.A229G, and MOGAT3/g.G1627A and growth traits. The polymerase chain reaction–restriction fragment length polymorphism (PCR-RFLP) method was performed for genotyping the MOGAT1, MOGAT3/g.A229G, and MOGAT3/g.G1627A genes' locus using the TaqI, MspI, and BsuRI restriction enzyme. The allele frequency of A and G of the MOGAT1 locus was 0.79 and 0.21, respectively, while the genotype frequency was 0.65, 0.28, and 0.07 for AA, AG, and GG, respectively. While the allele and genotype frequencies of the MOGAT3/g.A229G locus were 00.57(A1), 0.43(G1), 0.35(A1A1), 0.45(A1G1), and 0.20(G1G1), the allele and genotype frequencies of the MOGAT3/g.G1627A locus were 0.49(A2), 0.51(G2), 0.25(A2A2), 0.49(A2G2), and 0.26(G2G2). Chi-square analysis showed that MOGAT3/g.G1627A distribution was at the Hardy–Weinberg disequilibrium (p < 0.05), and MOGAT1 and MOGAT3/g.A229G distribution was at the Hardy–Weinberg equilibrium (p > 0.05). In total, two statistical methods (general linear model (GLM) and PROC MIXED) were used to identify an association between gene locus and growth traits. An association analysis showed a statistically significant difference between the MOGAT1 and body weight, body length, and chest circumference, MOGAT3/g.A229G with average daily gain (ADG) and withers height, and MOGAT3/g.G1627A with body weight and body length (p < 0.05). The results confirmed that the MOGAT1, MOGAT3/g.A229G, and MOGAT3/g.G1627A locus are strong candidate genes that could be considered molecular markers for growth traits in cattle breeding.

Identification and functional characterization of multiple inositol polyphosphate phosphatase1 (Minpp1) isoform-2 in exosomes with potential to modulate tumor microenvironment

Inositol polyphosphates (InsPs) play key signaling roles in diverse cellular functions, including calcium homeostasis, cell survival and death. Multiple inositol polyphosphate phosphatase 1 (Minpp1) affects the cellular levels of InsPs and cell functions. The Minpp1 is an endoplasmic reticulum (ER) resident but localizes away from its cytosolic InsPs substrates. The current study examines the heterogeneity of Minpp1 and the potential physiologic impact of Minpp1 isoforms, distinct motifs, subcellular distribution, and enzymatic potential. The NCBI database was used to analyze the proteome diversity of Minpp1 using bioinformatics tools. The analysis revealed that translation of three different Minpp1 variants resulted in three isoforms of Minpp1 of varying molecular weights. A link between the minpp1 variant-2 gene and ER-stress, using real-time PCR, suggests a functional similarity between minpp1 variant-1 and variant-2. A detailed study on motifs revealed Minpp1 isoform-2 is the only other isoform, besides isoform-1, that carries a phosphatase motif for InsPs hydrolysis but no ER-retention signal. The confocal microscopy revealed that the Minpp1 isoform-1 predominantly localized near the nucleus with a GRP-78 ER marker, while Minpp1 isoform-2 was scattered more towards the cell periphery where it co-localizes with the plasma membrane-destined multivesicular bodies biomarker CD63. MCF-7 cells were used to establish that Minpp1 isoform-2 is secreted into exosomes. Brefeldin A treatment resulted in overexpression of the exosome-associated Minpp1 isoform-2, suggesting its secretion via an unconventional route involving endocytic-generated vesicles and a link to ER stress. Results further demonstrated that the exosome-associated Minpp1 isoform-2 was enzymatically active. Overall, the data support the possibility that an extracellular form of enzymatically active Minpp1 isoform-2 mitigates any anti-proliferative actions of extracellular InsPs, thereby also impacting the makeup of the tumor microenvironment.

In silico identification of variations in microRNAs with a potential impact on dairy traits using whole ruminant genome SNP datasets

MicroRNAs are small noncoding RNAs that have important roles in the lactation process and milk biosynthesis. Some polymorphisms have been studied in various livestock species from the perspective of pathology or production traits. To target variants that could be the causal variants of dairy traits, genetic variants of microRNAs expressed in the mammary gland or present in milk and localized in dairy quantitative trait loci (QTLs) were investigated in bovine, caprine, and ovine species. In this study, a total of 59,124 (out of 28 millions), 13,427 (out of 87 millions), and 4761 (out of 38 millions) genetic variants in microRNAs expressed in the mammary gland or present in milk were identified in bovine, caprine, and ovine species, respectively. A total of 4679 of these detected bovine genetic variants are located in dairy QTLs. In caprine species, 127 genetic variants are localized in dairy QTLs. In ovine species, no genetic variant was identified in dairy QTLs. This study leads to the detection of microRNA genetic variants of interest in the context of dairy production, taking advantage of whole genome data to identify microRNA genetic variants expressed in the mammary gland and localized in dairy QTLs.

Molecular mechanisms of fat deposition: IL-6 is a hub gene in fat lipolysis, comparing thin-tailed with fat-tailed sheep breeds

Tail fat content affects meat quality and varies significantly among different breeds of sheep. Ghezel (fat-tailed) and Zel (thin-tailed) are two important Iranian local sheep breeds with different patterns of fat storage. The current study presents the transcriptome characterization of tail fat using RNA sequencing in order to get a better comprehension of the molecular mechanism of lipid storage in the two mentioned sheep breeds. Seven (Zel  =  4 and Ghezel  =  3) 7-month-old male lambs were used for this experiment. The results of sequencing were analyzed with bioinformatics methods, including differentially expressed genes (DEGs) identification, functional enrichment analysis, structural classification of proteins, protein-protein interaction (PPI) and network and module analyses. Some of the DEGs, such as LIPG, SAA1, SOCS3, HIF-1 α , and especially IL-6, had a close association with lipid metabolism. Furthermore, functional enrichment analysis revealed pathways associated with fat deposition, including "fatty acid metabolism", "fatty acid biosynthesis" and "HIF-1 signaling pathway". The structural classification of proteins showed that major down-regulated DEGs in the Zel (thin-tailed) breed were classified under transporter class and that most of them belonged to the solute carrier transporter (SLC) families. In addition, DEGs under the transcription factor class with an important role in lipolysis were up-regulated in the Zel (thin-tailed) breed. Also, network analysis revealed that IL-6 and JUNB were hub genes for up-regulated PPI networks, and HMGCS1, VPS35 and VPS26A were hub genes for down-regulated PPI networks. Among the up-regulated DEGs, the IL-6 gene seems to play an important role in lipolysis of tail fat in thin-tailed sheep breeds via various pathways such as tumor necrosis factor (TNF) signaling and mitogen-activated protein kinase (MAPK) signaling pathways. Due to the probable role of the IL-6 gene in fat lipolysis and also due to the strong interaction of IL-6 with the other up-regulated DEGs, it seems that IL-6 accelerates the degradation of lipids in tail fat cells.

Genome-wide identification of Diacylglycerol Acyltransferases (DGAT) family genes influencing Milk production in Buffalo

Background

The diacylglycerol acyltransferases (DGAT) are a vital group of enzymes in catalyzing triacylglycerol biosynthesis. DGAT genes like DGAT1 and DGAT2, have been identified as two functional candidate genes affecting milk production traits, especially for fat content in milk. Buffalo milk is famous for its excellent quality, which is rich in fat and protein content. Therefore, this study aimed to characterize DGAT family genes in buffalo and to find candidate markers or DGAT genes influencing lactation performance.

Results

We performed a genome-wide study and identified eight DGAT genes in buffalo. All the DGAT genes classified into two distinct clades (DGAT1 and DGAT2 subfamily) based on their phylogenetic relationships and structural features. Chromosome localization displayed eight buffalo DGAT genes distributed on five chromosomes. Collinearity analysis revealed that the DGAT family genes were extensive homologous between buffalo and cattle. Afterward, we discovered genetic variants loci within the genomic regions that DGAT genes located in buffalo. Seven haplotype blocks were constructed and were associated with buffalo milk production traits. Single marker association analyses revealed four most significant single nucleotide polymorphisms (SNPs) mainly affecting milk protein percentage or milk fat yield in buffalo. Genes functional analysis indicated that these DGAT family genes could influence lactation performance in the mammal through regulating lipid metabolism.

Conclusion

In the present study, we performed a comprehensive analysis for the DGAT family genes in buffalo, which including identification, structural characterization, phylogenetic classification, chromosomal distribution, collinearity analysis, association analysis, and functional analysis. These findings provide useful information for an in-depth study to determine the role of DGAT family gens play in the regulation of milk production and milk quality improvement in buffalo.