A novel allelic donkey ß-lactoglobulin I protein isoform generated by a non-AUG translation initiation codon is associated with a nonsynonymous SNP.

ß-lactoglobulin I (ß-LG I) is one of the most important whey proteins in donkey milk. However, to our knowledge, there has been no study focusing on the full nucleotide sequences of this gene (BLG I). Current investigation of donkey BLG I gene is very limited with only 2 variants (A and B) characterized so far at the protein level. Recently, a new ß-LG I variant, with a significantly higher mass (+1,915 Da) than known variants has been detected. In this study, we report the whole nucleotide sequence of the BLG I gene from 2 donkeys, whose milk samples are characterized by the ß-LG I SDS-PAGE band with a normal electrophoretic mobility (18,514.25 Da, ß-LG I B1 form) the first, and by the presence of a unique ß-LG I band with a higher electrophoretic mobility (20,428.5 Da, ß-LG I D form) the latter. A high genetic variability was found all over the 2 sequenced BLG I alleles. In particular, 16 polymorphic sites were found in introns, one in the 5' flanking region, 3 SNPs in the 5' untranslated region and one SNP in the coding region (g.1871G > A) located at the 40th nucleotide of exon 2 and responsible for the AA substitutions p.Asp28 > Asn in the mature protein. Two SNPs (g.920-922CAC > TGT and g.1871G/A) were genotyped in 93 donkeys of 2 Italian breeds (60 Ragusana and 33 Amiatina, respectively) and the overall frequencies of g.920-922CAC and g.1871A were 0.3065 and 0.043, respectively. Only the rare allele g.1871A was observed to be associated with the slower migrating ß-LG I. Considering this genetic diversity and those found in the database, it was possible to deduce at least 5 different alleles (BLG I A, B, B1, C, D) responsible for 4 potential ß-LG I translations. Among these alleles, B1 and D are those characterized in the present research, with the D allele of real novel identification. Haplotype data analysis suggests an evolutionary pathway of donkey BLG I gene and a possible phylogenetic map is proposed. Analyses of mRNA secondary structure showed relevant changes in the structures, as consequence of the g.1871G > A polymorphism, that might be responsible for the recognition of an alternative initiation site providing an additional signal peptide. The extension of 19 AA sequence to the mature protein, corresponding to the canonical signal peptide with an additional alanine residue, is sufficient to provide the observed molecular weight of the slower migrating ß-LG I encoded by the BLG I D allele.

Remarkable genetic diversity detected at river buffalo prolactin receptor (PRLR) gene and association studies with milk fatty acid composition.

Prolactin is an anterior pituitary peptide hormone involved in many different endocrine activities and is essential for reproductive performance. This action is mediated by its receptor, the prolactin receptor, encoded by the PRLR gene. In this study, we sequenced and characterized the Mediterranean river buffalo PRLR gene (from exon 3 to 10), and we found remarkable genetic diversity. In particular, we found 24 intronic polymorphisms and 13 exonic SNPs, seven of which were non-synonymous. Furthermore, the polymorphisms identified in the 3'-UTR were investigated to establish their possible influence on microRNA binding sites. Considering all the amino acid changes and the observed allelic combinations, it is possible to deduce at least six different translations of the buffalo prolactin receptor and, consequently, the presence at the PRLR gene of at least six alleles. Furthermore, we identified a deletion of a CACTACC heptamer between nucleotides 1102 and 1103 of exon 10 (3'-UTR), and we developed an allele-specific PCR to identify the carriers of this genetic marker. Finally, the SNP g.11188A>G, detected in exon 10 and responsible for the amino acid replacement p.His328Arg, was genotyped in 308 Italian Mediterranean river buffaloes, and an association study with milk fat traits was carried out. The statistical analysis showed a tendency that approached significance for the AA genotype with higher contents of odd branched-chain fatty acids. Thus, our results suggest that the PRLR gene is a good candidate for gene association studies with qualitative traits related to buffalo milk production.