Sequence analysis and genetic variability of stearoyl CoA desaturase (SCD) gene in the Italian Mediterranean river buffalo

CSN1S1, CSN3 and LPL: Three Validated Gene Polymorphisms Useful for More Sustainable Dairy Production in the Mediterranean River Buffalo

The search for DNA polymorphisms useful for the genetic improvement of dairy farm animals has spanned more than 40 years, yielding relevant findings in cattle for milk traits, where the best combination of alleles for dairy processing has been found in casein genes and in DGAT1. Nowadays, similar results have not yet been reached in river buffaloes, despite the availability of advanced genomic technologies and accurate phenotype records. The aim of the present study was to investigate and validate the effect of four single nucleotide polymorphisms (SNP) in the CSN1S1, CSN3, SCD and LPL genes on seven milk traits in a larger buffalo population. These SNPs have previously been reported to be associated with, or affect, dairy traits in smaller populations often belonging to one farm. A total of 800 buffaloes were genotyped. The following traits were individually recorded, monthly, throughout each whole lactation period from 2010 to 2021: daily milk yield (dMY, kg), protein yield (dPY, kg) and fat yield (dFY, kg), fat and protein contents (dFP, % and dPP, %), somatic cell count (SCC, 103 cell/mL) and urea (mg/dL). A total of 15,742 individual milk test day records (2496 lactations) were available for 680 buffalo cows, with 3.6 ± 1.7 parities (from 1 to 13) and an average of 6.1 ± 1.2 test day records per lactation. Three out four SNPs in the CSN1S1, CSN3 and LPL genes were associated with at least one of analyzed traits. In particular, the CSN1S1 (AJ005430:c.578C>T) gave favorable associations with all yield traits (dMY, p = 0.022; dPY, p = 0.014; dFY, p = 0.029) and somatic cell score (SCS, p = 0.032). The CSN3 (HQ677596: c.536C>T) was positively associated with SCS (p = 0.005) and milk urea (p = 0.04). Favorable effects on daily milk yield (dMY, p = 0.028), fat (dFP, p = 0.027) and protein (dPP, p = 0.050) percentages were observed for the LPL. Conversely, the SCD did not show any association with milk traits. This is the first example of a confirmation study carried out in the Mediterranean river buffalo for genes of economic interest in the dairy field, and it represents a very important indication for the preselection of young bulls destined for breeding programs aimed at more sustainable dairy production.

Genetic variability among and within domestic Old and New World camels at the α-lactalbumin gene (LALBA) reveals new alleles and polymorphisms responsible for differential expression

α-Lactalbumin (α-LA), which is encoded by the LALBA gene, is a major whey protein that binds to Ca2+ and facilitates lactose synthesis as a regulatory subunit of the synthase enzyme complex. In addition, it has been shown to play central roles in immune modulation, cell-growth regulation, and antimicrobial activity. In this study, a multitechnical approach was used to fully characterize the LALBA gene and its variants in both coding and regulatory regions for domestic camelids (dromedary, Bactrian camel, alpaca, and llama). The gene analysis revealed a conserved structure among the camelids, but a slight difference in size (2,012 bp on average) due to intronic variations. Promoters were characterized for the transcription factor binding sites (11 found in total). Intraspecies sequence comparison showed 36 SNPs in total (2 in the dromedary, none in the Bactrian camel, 22 in the alpaca, and 12 in the llama), whereas interspecies comparison showed 86 additional polymorphic sites. Eight SNPs were identified as trans-specific polymorphisms, and 2 of them (g.112A>G and g.1229A>G) were particularly interesting in the New World camels. The first creates a new binding site for transcription factor SP1. An enhancing effect of the g.112G variant on the expression was demonstrated by 3 independent pGL3 gene reporter assays. The latter is responsible for the p.78Ile>Val AA replacement and represents novel allelic variants (named LALBA A and B). A link to protein variants has been established by isoelectric focusing (IEF), and bioinformatics analysis revealed that carriers of valine (g.1229G) have a higher glycosylation rate. Genotyping methods based on restriction fragment length polymorphism (PCR-RFLP) were set up for both SNPs. Overall, adenine was more frequent (0.54 and 0.76) at both loci. Four haplotypes were found, and the AA and GA were the most common with a frequency of 0.403 and 0.365, respectively. Conversely, a putative biological gain characterizes the haplotype GG. Therefore, opportunities for rapid directional selection can be realized if this haplotype is associated with favorable milk protein properties. This study adds knowledge at the gene and protein level for α-LA (LALBA) in camelids and importantly contributes to a relatively unexplored research area in these species.

Complete CSN1S2 Characterization, Novel Allele Identification and Association With Milk Fatty Acid Composition in River Buffalo

The αs2-casein is one of the phosphoproteins secreted in all ruminants' milk, and it is the most hydrophilic of all caseins. However, this important gene (CSN1S2) has not been characterized in detail in buffaloes with only two alleles detected (reported as alleles A and B), and no association studies with milk traits have been carried out unlike what has been achieved for other species of ruminants. In this study, we sequenced the whole gene of two Mediterranean river buffalo homozygotes for the presence/absence of the nucleotide C (g.7539G>C) realized at the donor splice site of exon 7 and, therefore, responsible for the skipping of the same exon at mRNA level (allele B). A high genetic variability was found all over the two sequenced CSN1S2 alleles. In particular, 74 polymorphic sites were found in introns, six in the promoter, and three SNPs in the coding region (g.11072C>T, g.12803A>T, and g.14067A>G) with two of them responsible for amino acid replacements. Considering this genetic diversity, those found in the database and the SNP at the donor splice site of exon 7, it is possible to deduce at least eight different alleles (CSN1S2 A, B, B1, B2, C, D, E, and F) responsible for seven different possible translations of the buffalo αs2-casein. Haplotype data analysis suggests an evolutionary pathway of buffalo CSN1S2 gene consistent with our proposal that the published allele CSN1S2 A is the ancestral αs2-CN form, and the B2 probably arises from interallelic recombination (single crossing) between the alleles D and B (or B1). The allele CSN1S2 C is of new identification, while CSN1S2 B, B1, and B2 are deleted alleles because all are characterized by the mutation g.7539G>C. Two SNPs (g.7539G>C and g.14067A>G) were genotyped in 747 Italian buffaloes, and major alleles had a relative frequency of 0.83 and 0.51, respectively. An association study between these SNPs and milk traits including fatty acid composition was carried out. The SNP g.14067A>G showed a significant association (P < 0.05) on the content of palmitic acid in buffalo milk, thus suggesting its use in marker-assisted selection programs aiming for the improvement of buffalo milk fatty acid composition.

Fatty acid biosynthesis and transcriptional regulation of Stearoyl-CoA Desaturase 1 (SCD1) in buffalo milk

Background

Buffalo milk is considered as a highly nutritious food owing to its higher contents of fatty acids (FA) and rich nutrient profile. Higher fat contents of buffalo milk make it suitable for processing to develop various healthy and nutritious products. Moreover, buffalo milk contains more unsaturated FAs (UFA) such as oleic and linolenic acid, which are important from the human health point of view owing to their desirable physiological effects. However, inadequate information is available about the chemical composition and mechanism of FA synthesis in buffalo milk. In this study, we hypothesized that expression of SCD1 gene could alter the biosynthesis of FA in epithelial cells of mammary gland and subsequently affect the FA contents in buffalo milk. We investigated the transcriptional and biological role of Stearoyl-CoA Desaturase 1 (SCD1) in the buffalo mammary epithelial cells (BMECs) during FA and triacylglycerol (TAG) synthesis.

Results

Results revealed that unsaturated fatty acid contents were much higher in concentration in buffalo milk as compared to Holstein cow. Significant increase in the expression level of FAS, ACACA, SREBP1, PPARG, GPAT, and AGPAT genes was observed in response to altered expression of SCD1 in buffalo milk. Moreover, change in SCD1 gene in BMECs also mediated the expression of genes related to FA biosynthesis subsequently leading to alter the FA composition. Overexpression of SCD1 significantly increased the expression of genes associated with FA and TAG synthesis leading to enhance FA and unsaturated FA contents in BMECs. However, down-regulation of SCD1 exhibited opposite consequences.

Conclusion

Our study provides mechanistic insights on transcriptional regulation of SCD1 to alter FA and TAG synthesis through directly or indirectly mediating biosynthesis and metabolic pathways in BMECs. We provide preliminary findings regarding engineering of FA contents in buffalo milk through SCD1 signaling.

The single nucleotide polymorphism g.133A>C in the stearoyl CoA desaturase gene (SCD) promoter affects gene expression and quali-quantitative properties of river buffalo milk

The stearoyl-CoA desaturase (SCD) gene has been investigated in depth in ruminants because of its effect on milk fat composition. In river buffalo, the single nucleotide polymorphism (SNP) g.133A>C in the gene promoter has been associated with milk quality and yield. However, the biological reason for such effects remains unexplored. In this study, we combined mRNA profile analysis, an electromobility shift assay, and quantitative PCR to elucidate the role of this SNP on gene transcription and its effects on milk fat traits. A preliminary genotyping of g.133A>C was carried out on a group of 303 river buffaloes to choose individuals for the downstream applications. Analysis of allele frequencies showed an increase in the minor allele C (0.25) compared with previous findings (0.16). Six animals (2 for each genotype) were chosen for cloning and 216 positive cDNA recombinant clones for SCD (72 per genotype) were analyzed by PCR. All clones showed the same length on agarose gel; therefore, random clones were chosen for sequencing. No qualitative differences were found and all gene transcripts assembled correctly. An electrophoretic mobility shift assay was performed to evaluate the binding of the transcription factor Sp1 to DNA sequences including g.133A>C. Genotype CC showed a higher binding (mean ± standard error of the mean) than genotype AA in 2 different conditions [Enzo buffer (EB), Enzo Life Science Inc., Farmingdale, NY: 201.77 ± 4.06 vs. 141.65 ± 3.77 band intensity values and Poletto buffer (PB): 95.90 ± 1.15 vs. 67.30 ± 2.14 band intensity values]. The subsequent quantitative PCR confirmed the upregulation of the CC genotype compared with the AA and AC genotypes. The association study with milk fat traits revealed a favorable effect of allele C. The heterozygous genotype had the highest values for monounsaturated fatty acids, oleic acid (C18:1 cis-9), polyunsaturated fatty acids, and odd- and branched-chain fatty acids, and the lowest values for saturated fatty acids and atherogenic and thrombogenic indices; the heterozygous genotype differed significantly from the AA genotype. The AC genotype has previously been associated with higher milk yield. Therefore, the g.133A>C SNP is a marker with dual effects and is an interesting candidate for assisted selection programs in river buffalo. These data clarified the biological role of the SNP g.133A>C in the SCD promoter and how it affects gene function, providing important knowledge on the genetic background of lipid metabolism, including the future possibility of selecting alleles with quantitatively or qualitatively favorable effects.

Extruded linseed alone or in combination with fish oil modifies mammary gene expression profiles in lactating goats

Nutrition is a major factor that regulates ruminant milk components, particularly its fatty acid (FA) composition, which is an important determinant of milk nutritional quality. In the mammary gland, milk component biosynthesis involves a large number of genes under nutritional regulation that are not well understood. Thus, the objective of the present study was to evaluate the effects of extruded linseeds (EL) alone or in combination with fish oil (ELFO) on goat mammary gene expression. In total, 14 goats were fed one of the following three diets: a natural grassland hay basal diet (CTRL) alone, CTRL supplemented with 530 g/day of EL, or 340 g/day of EL plus 39 g/day of fish oil (ELFO). Mammary secretory tissues were collected after slaughter on day 28, to determine the expression of 14 lipogenic genes and five lipogenic enzyme activities and transcriptomic profiles. The mRNA abundance decreased for SCD1 (P<0.1) with ELFO v. CTRL, and for AZGP1 (P<0.1) and ACSBG1 (P<0.05) decreased with EL v. ELFO and the CTRL diets (only for ACSBG1), respectively. Transcriptomic analyses performed using a bovine microarray revealed 344 and 314 differentially expressed genes (DEG) in the EL and ELFO diets, respectively, compared with the CTRL diet, with 76 common DEGs. In total, 21 and 27 DEGs were involved in lipid metabolism and transport class in the EL and ELFO v. the CTRL diets, respectively, with eight common genes (ALDH3B1, ALDH18A1, DGKD1, ENPP1, IL7, NSMAF, PI4KA and SERINC5) down-regulated by these two treatments. In EL v. CTRL diets, a gene network related to lipid metabolism and transport was detected. Although this network was not detected in the ELFO v. CTRL analysis, five genes known to be involved in lipid metabolism and transport were up-regulated (SREBF1, PPARG and GPX4) or down-regulated (FABP1 and ENPP6) by ELFO. The protein metabolism and transport biological processes were largely altered by both EL and ELFO v. CTRL diets without changes in major milk protein secretion. Amino acid metabolism was highlighted as an enriched network by Ingenuity Pathway Analysis and was similar to cellular growth and proliferation function. Two regulation networks centered on the estrogen receptor (ESR1) and a transcriptional factor (SP1) were identified in EL and ELFO v. CTRL diets. In conclusion, these results show that these two supplemented diets, which largely changed milk FA composition, had more effects on mRNA linked to protein metabolism and transport pathways than to lipid metabolism, and could affect mammary remodeling.

Quality of buffalo milk as affected by dietary protein level and flaxseed supplementation

The aim of the present research was to evaluate the effects of protein level and flaxseed supplementation on the yield and quality of buffalo milk. In particular, the fatty acid profile of milk from buffalo cows subjected to different diets has been investigated. A 2×3 factorial design was tested with buffalo cows receiving 2 dietary crude protein (CP) and 3 flaxseed (FS) supplementation levels. Treatments were (1) low dietary CP level [12% of dry matter (DM)] and no flaxseed supplementation (LP); (2) low dietary CP level (12% of DM) and low flaxseed supplementation (500g/d) (LPFS500); (3) low dietary CP level (12% of DM) and moderate flaxseed supplementation (1,000g/d) (LPFS1000); (4) moderate dietary CP level (15% of DM) and no flaxseed supplementation (MP); (5) moderate dietary CP level (15% of DM) and low flaxseed supplementation (500g/d) (MPFS500); and (6) moderate dietary CP level (15% of DM) and moderate flaxseed supplementation (1,000g/d) (MPFS1000). Milk protein and casein were affected by flaxseed supplementation being higher in MP, intermediate in LP, and lower in flaxseed-supplemented diets. However, the results from the present study highlighted that low protein diets sustained milk yield, protein, and casein synthesis in milk when whole flaxseed was administered. Short-chain fatty acids, in particular C8:0 and C10:0, were the lowest in milk from buffalo cows fed the highest level of flaxseed supplementation. Medium-chain fatty acids were the lowest in FS1000, intermediate in FS500, and the highest in the HP and LP groups. Long-chain fatty acids were the highest in FS1000, intermediate in FS500 groups, and the lowest in milk from buffalo receiving no flaxseed supplementation. Protein level of the diet influenced the percentage of C18:0, which was higher in MP than LP groups. Total conjugated linoleic acid content evidenced the same trend of long-chain fatty acids, with an increase of about 7% in FL500 and of 22% in FL1000 than the control. Apart from protein level of the diet, atherogenic index, thrombogenic index, and n-6/n-3 were the lowest in FS1000 groups; thrombogenic index and n-6/n-3 were intermediate in milk from animals receiving FS500. Nutritional value of the acidic profile in buffalo milk is influenced by flaxseed supplementation, and its improvement reflects the level of dietary flaxseed supplementation.

Genetic variability detected at the lactoferrin locus (LTF) in the Italian Mediterranean river buffalo

Lactoferrin (LTF) is a multi-functional protein belonging to the whey protein fractions of milk. The gene LTF encoding for such protein is considered a potential candidate for body measurement, milk composition and yield. This study reports on the genetic variability at the LTF locus in the Italian Mediterranean river buffalo and its possible association with milk yield. Eleven polymorphic sites were found in the DNA fragment spanning exons 15–16. In particular, intron 15 was extremely polymorphic with nine single nucleotide polymorphisms (SNPs) detected, whereas the remaining two SNPs were exonic mutations (g.88G > A at exon 15 and g.1351G > A at exon 16) and both synonymous. The genotyping of the informative samples evidenced three haplotypes, whose frequencies were 0.6, 0.3 and 0.1, respectively, whereas the analysis of the exonic SNPs showed a perfect condition of linkage disequilibrium (g.88A/g.1351G and g.88G/g.1351A). The association study carried out by using the SNP g.88G > A showed that buffalo LTF gene has no statistically significant influence on daily milk yield. This study adds knowledge to the genetic variability of a species less investigated than the other ruminant species. These findings may serve as a useful tool for large-scale screening of buffalo populations.

Mediterranean river buffalo CSN1S1 gene: search for polymorphisms and association studies

The aim of the present study was to investigate the variability at CSN1S1 locus of the Italian Mediterranean river buffalo and to study possible allele effects on milk yield and its composition. Effects of parity, calving season and month of production were also evaluated. Three single-nucleotide polymorphisms were detected. The first mutation, located at position 89 of the 17th exon (c.628C>T), is responsible for the amino acid change p.Ser178 (B allele)/Leu178 (A allele). The other two polymorphisms, detected at the positions 144 (c.882G>A) and 239 (c.977A>G) of 19th exon, respectively, are silent (3ʹ UTR, untranslated region). Associations between the CSN1S1 genotypes and milk production traits were investigated using 4122 test day records of 503 lactations from 175 buffalo cows. Milk yield, fat and protein percentages were analysed using a mixed linear model. A significant association between the c.628C>T SNP and the protein percentage was found. In particular, the CC genotype showed an average value ~0.04% higher than the CT and TT genotypes. The allele substitution effect of cytosine into thymine was –0.014, with a quite low (0.3%) protein percentage contribution to total phenotypic variance. A large dominance effect was detected. Characterisation of the CSN1S1 transcripts and a method based on MboI amplification created restriction site PCR for a rapid genotyping of c.628C>T are provided.

A single nucleotide polymorphism in the promoter region of river buffalo stearoyl CoA desaturase gene (SCD) is associated with milk yield

An association study between the milk yield trait and the stearoyl-CoA desaturase (SCD) polymorphism (g.133A > C) in Italian Mediterranean river buffalo was carried out. A full characterization of the river buffalo SCD promoter region was presented. Genotyping information was provided and a quick method for allelic discrimination was developed. The frequency of the C allele was 0·16. Test-day (TD) records (43 510) of milk production belonging to 226 lactations of 169 buffalo cows were analysed with a mixed linear model in order to estimate the effect of g.133A > C genotype, as well as the effect of parity and calving season. The SCD genotype was significantly associated with milk yield (P = 0·02). The genotype AC showed an over-dominance effect with an average daily milk yield approximately 2 kg/d higher than CC buffaloes. Such a difference represents about 28% more milk/d. The effect of the genotype was constant across lactation stages. The contribution of SCD genotype (r2SCD) to the total phenotypic variance in milk yield was equal to 0·12. This report is among the first indications of genetic association between a trait of economic importance in river buffalo. Although such results need to be confirmed with large-scale studies in the same and other buffalo populations, they might offer useful indications for the application of MAS programmes in river buffalo and in the future they might be of great economic interest for the river buffalo dairy industry.