New genetic polymorphisms within ovine β- and αS2-caseins

Single and combined effects of CSN1S1 and CSN2-casein genes on Awassi sheep milk quantity and quality

Background and Aim: Milk produced from Awassi sheep is of high nutritive value; its production is relatively low in Awassi sheep, so the genetic improvement programs targeted milk production and its components are of high importance, especially when using genes that have an important signal to milk traits. This study was aimed at assessing the influence of alpha S1 (CSN1S1) and beta-casein (CSN2) genes genotypes interaction on Awassi ewes milk productivity. Materials and Methods: A total number of 391 milk yield and its composition records (taken through five consecutive years, 2007-2011) of 167 ewes were utilized for this study. DNA samples were extracted from the ewe's blood samples, then the polymerase chain reaction products of alpha S1 (CSN1S1) and beta-casein (CSN2) genes were sequenced. The obtained sequences were analyzed; thereafter, the detected variants were tested for their possible association with milk traits. Results: The CSN1S1 and CSN2 variants allelic frequencies were 0.85 and 0.15, and 0.95 and 0.05, respectively. Lactose and solid not fat (SNF) % were associated with TC CSN1S1 genotypes. No association was found among CSN1S1 polymorphic genotypes with milk production, lactose, and SNF % were associated with TC CSN1S1 genotypes. Ewes of CSN2 AC genotype showed higher milk production traits, while no association was found between milk composition traits and CNS2 genotypes. Nevertheless, CSN1S1∗CSN2 interaction showed the highest SNF, fat percentages, and milk production. Conclusion: The substantial interaction effects between CSN1S1×CSN2 genes were significantly affected the amount of milk, fat, and SNF% produced. The detected variants should be included in the breeding programs of Awassi sheep that are designed for improving their milk quantity and quality.

DNA polymorphisms of FGFBP1, leptin, κ-casein, and αs1-casein genes and their association with reproductive performance in dromedary she-camels

In livestock, selection programs on reproductive traits are time consuming due to long generation interval and low heritability; however understanding the genetic factors that influence reproductive performance facilitates an accurate selection and genetic improvement of these quantitative traits. Therefore, the aim of this study was to identify the effect of lactation order and polymorphisms of FGFBP1, leptin, κ-casein, and αs1-casein genes on reproductive performance in dromedary she-camels. Blood samples were collected from one hundred she-camels using vacutainer tube containing EDTA as an anticoagulant for DNA extraction. PCR-DNA sequencing approach was performed to investigate the polymorphisms in fragments for coding sequence (CDS) parts of FGFBP1 (539-bp) and leptin (359-bp), 5'-flanking region of κ-casein (488-bp) and intervening non-coding fragment spanning from exon 4 to exon 6 (924-bp) of αs1-casein genes. Results revealed that there were single nucleotide polymorphisms (SNPs) in the investigated genes between enrolled she-camels except for leptin gene that elicited a monomorphic pattern. The identified SNPs allocated camels into different groups according to which animal harbored the explored SNPs. Multiple linear regression analysis (MLR) reveled that identified SNPs and parity order had significant effects on age at first calving (AFC), days open (DO), calving interval (CI), number of services/conception (NS/C) and milk yield traits. Determination coefficients (R²) value indicated also that variability in the phenotypic measurements of studied traits might be correlated with the identified SNPs in the reproduction related genes. This study highlights that parity order has a remarkable effect on reproductive performance. Single nucleotide polymorphisms in the investigated genes could be also used as proxy predictors and genetic markers that permit a marker assisted selection for reproductive performance traits in dromedary she-camels.

Variant discovery in the sheep milk transcriptome using RNA sequencing

BACKGROUND

The identification of genetic variation underlying desired phenotypes is one of the main challenges of current livestock genetic research. High-throughput transcriptome sequencing (RNA-Seq) offers new opportunities for the detection of transcriptome variants (SNPs and short indels) in different tissues and species. In this study, we used RNA-Seq on Milk Sheep Somatic Cells (MSCs) with the goal of characterizing the genetic variation within the coding regions of the milk transcriptome in Churra and Assaf sheep, two common dairy sheep breeds farmed in Spain.

RESULTS

A total of 216,637 variants were detected in the MSCs transcriptome of the eight ewes analyzed. Among them, a total of 57,795 variants were detected in the regions harboring Quantitative Trait Loci (QTL) for milk yield, protein percentage and fat percentage, of which 21.44% were novel variants. Among the total variants detected, 561 (2.52%) and 1,649 (7.42%) were predicted to produce high or moderate impact changes in the corresponding transcriptional unit, respectively. In the functional enrichment analysis of the genes positioned within selected QTL regions harboring novel relevant functional variants (high and moderate impact), the KEGG pathway with the highest enrichment was "protein processing in endoplasmic reticulum". Additionally, a total of 504 and 1,063 variants were identified in the genes encoding principal milk proteins and molecules involved in the lipid metabolism, respectively. Of these variants, 20 mutations were found to have putative relevant effects on the encoded proteins.

CONCLUSIONS

We present herein the first transcriptomic approach aimed at identifying genetic variants of the genes expressed in the lactating mammary gland of sheep. Through the transcriptome analysis of variability within regions harboring QTL for milk yield, protein percentage and fat percentage, we have found several pathways and genes that harbor mutations that could affect dairy production traits. Moreover, remarkable variants were also found in candidate genes coding for major milk proteins and proteins related to milk fat metabolism. Several of the SNPs found in this study could be included as suitable markers in genotyping platforms or custom SNP arrays to perform association analyses in commercial populations and apply genomic selection protocols in the dairy production industry.

Variations at regulatory regions of the milk protein genes are associated with milk traits and coagulation properties in the Sarda sheep

Regulatory variation at the ovine casein genes could have important effects on the composition and coagulation properties of milk. Herewith, we have partially resequenced the promoters and the 3'-UTR of the four casein genes in 25 Sarda sheep. Alignment of these sequences allowed us to identify a total of 29 SNPs. This level of polymorphism (one SNP every 250 bp) is remarkably high if compared with SNP densities estimated in human genic regions (approximately one SNP per bp). The 29 SNPs identified in our resequencing experiment, plus three previously reported SNPs mapping to the lactalbumin, alpha (LALBA) and β-lactoglobulin (BLG, also known as progestagen-associated endometrial protein, PAEP) genes, were genotyped with a multiplex TaqMan Open Array Real-Time PCR assay in 760 Sarda sheep with records for milk composition and coagulation properties. Association analysis revealed the existence of significant associations of CSN1S2 and CSN3 genotypes with milk protein and casein contents. Moreover, genotypes at CSN1S1 were significantly associated with rennet coagulation time, curd firming time and curd firmness, whereas CSN2 was associated with curd firming time. These results suggest that SNPs mapping to the promoters and 3'-UTRs of ovine casein genes may exert regulatory effects on gene expression and that they could be used for improving sheep milk quality and technological traits at the population level through marker assisted selection.

Investigating the genetic polymorphism of sheep milk proteins: a useful tool for dairy production

Sheep is the second most important dairy species after cow worldwide, and especially in the Mediterranean and Middle East regions. In some countries, the difficult environmental conditions require a peculiar adaptation and, in these contexts, sheep are able to provide higher quality protein than cattle. In the least-developed countries, the amount of dairy sheep and ovine milk production is progressively increasing. In order to improve dairy productions, in particular those with local connotations, it is necessary to obtain in-depth information regarding milk quality and rheological properties. The genetic polymorphisms of milk proteins are often associated with quantitative and qualitative parameters in milk and are potential candidate markers that should be included in breeding strategies similar to those already available for cattle. Due to the current and growing interest in this topic and considering the large amount of new information, the aim of this study was to review the literature on sheep milk protein polymorphisms with a particular emphasis on recent findings in order to give scientists useful support. Moreover, the effects of different protein variants on milk yield and composition are discussed.

The β-casein in camels: molecular characterization of the CSN2 gene, promoter analysis and genetic variability

The β-casein is the most abundant protein in camel milk and its encoding gene (CSN2) is considered in other species a 'major' gene for the presence of alleles associated to different level of expression. In the present paper, we report for the first time the characterization of the nucleotide sequence of the whole β-casein-encoding gene (CSN2) plus 2,141 bp at the 5'-flanking region in Camelus dromedarius. The promoter region and the complete cDNA are also provided for the first time in Camelus bactrianus. The gene is spread over 7.8 kb and consists of 9 exons varying in length from 24 bp (exon 5) to 519 bp (exon 7), and 8 introns from 95 bp (intron 5) to 1,950 bp (intron 1). The composite response element (CoRE) region was identified in the promoter, whereas the presence of mature microRNA sequences improves the knowledge on the factors putatively involved in the gene regulation. A total of 46 polymorphic sites have been detected. The transition g.2126A>G falls within the TATA-box of dromedary CSN2 promoter with a putative influence on the transcription factor binding activity. The frequency of the G allele is 0.35 in a population of 180 she-camels belonging to 4 different ecotypes. In the same population, a conservative SNP (g.4175C>A) was found at the codon 7 of the signal peptide, whereas a comparative analysis with a cDNA sequence available in the database evidenced a missense SNP (g.4180T(Leu)>G(Arg)) at exon 2. Four SNPs were found in the bactrian camel. The SNP c.666G>A is responsible for the amino acid change Met(201)→Ile and it represents the first missense allele at the β-casein in camels. Finally, five interspersed repeated elements were identified at intronic level, whereas the presence of putative bio-functional peptides belonging to ACE-inhibitor and anti-oxidative families confirms the potential protective role of the camel milk for the human nutrition.

Characterization of β-casein gene in Indian riverine buffalo

The study aimed at characterization of buffalo β-casein gene and its promoter by PCR-SSCP analysis. Complete β-casein exon VII region analysis revealed two SSCP band patterns, with pattern-I representing predominant allele B (85%) present in homozygous (genotype BB) condition and pattern-II representing a rare allele A1 present in heterozygous condition (genotype A1B). Sequencing of two patterns revealed three nucleotide substitutions at codon 68, 151 and 193 of exon VII. The cDNA sequence of buffalo β-casein gene indicated three further nucleotide substitutions between allele A1 and B at codon 10, 39, and 41. Analysis of β-casein proximal promoter region (-350 upstream to +32) revealed four SSCP band patterns. These SSCP patterns corresponded to nucleotide substitutions at seven locations within 382 bp 5' UTR region of β-casein gene. Haplotype analysis suggested pattern-I of exon VII (wild type) was associated with three types of promoters and pattern-II of exon VII (rare type) corresponded to one exclusive type of promoter. The study suggested two haplotypes of exon VII and four haplotypes of promoter for buffalo β-casein.

Identification of novel single nucleotide polymorphisms in the DGAT1 gene of buffaloes by PCR-SSCP

Diacylglycerol O-acyltransferase 1 (DGAT1) is a microsomal enzyme that catalyzes the final step of triglyceride synthesis. The DGAT1 gene is a strong functional candidate for determining milk fat content in cattle. In this work, we used PCR-SSCP (polymerase chain reaction-single-strand conformation polymorphism) and DNA sequencing to examine polymorphism in the region spanning exon 7 to exon 9 of the DGAT1 gene in Murrah and Pandharpuri buffaloes. Three alleles (A, B and C) and four novel single-nucleotide polymorphisms were identified in the buffalo DGAT1 gene. The frequencies of the alleles differed between the two buffalo breeds, with allele C being present in Murrah but not in Pandharpuri buffalo. The allele variation detected in this work may influence DGAT1 expression and function. The results described here could be useful in examining the association between the DGAT1 gene and milk traits in buffalo.

Molecular genetic characterization of ovine CSN1S2 variants C and D reveal further important variability within CSN1S2

Within this study, the recently identified ovine CSN1S2 variants C and D were characterized at the molecular genetic level. Sequencing of the cDNA and of parts of the DNA identified several sequence differences within CSN1S2*C and D in comparison to CSN1S2*A and B. CSN1S2*C is characterized by two non-synonymous single nucleotide polymorphisms (SNPs) within exon 7 (c.178A>G, c.187G>T) leading to the amino acid substitutions p.Val45Ile and p.Ala48Ser. CSN1S2*D is caused by the SNP c.183G>C, leading to an amino acid replacement at position 46 (p.Arg46Ser). A very common c.527G>A-SNP within exon 15, resulting in the amino acid substitution p.Arg161His and producing the new variant CSN1S2*G, not detectable by isoelectric focusing and previously misidentified as CSN1S2*A, was also identified. On the basis of the identified sequence differences, a new nomenclature is proposed and a possible phylogenetic pathway shown for ovine CSN1S2 variants.