Gelsolin expression in sheep milk somatic cells during lactation

The identification of genes involved in phenotypes related to milk quality is important for both economic and health aspects in livestock production. The aim of this study was to assess the level of gelsolin gene expression in two breeds of dairy sheep – Sarda and Gentile – with pronounced differences in quantitative and qualitative milk traits. Gelsolin, a type of actin-modulating proteins is involved in the processes of actin remodeling during cell growth and apoptosis; therefore a role of this protein in mammary changes during lactation was here hypothesized. Individual milk samples were collected three times during lactation from 26 ewes of the two breeds. The differential gene expression of gelsolin in the two breeds and the three lactation times was estimated by quantitative PCR on RNA extracted from milk somatic cells. Correlations of gelsolin gene expression with milk yield and quality and days of lactation were also estimated. The results showed that gelsolin gene expression was significantly higher in the Sarda compared to the Gentile at each lactation stage, in agreement with the longer lactation duration and the higher daily milk yield of the first breed. Significant correlations of gelsolin gene expression were found with milk fat content in Sarda breed (−0.46, P<0.05). Gelsolin expression analysis confirmed the link between gelsolin gene function and milk fat content of sheep.     

Novel SNPs and INDEL polymorphisms in the 3′UTR of DGAT1 gene: in silico analyses and a possible association

Diacylglycerol–O–acyltransferase (DGAT1) gene encodes the rate-limiting enzyme of triglyceride synthesis. A polymorphism in this gene, DGAT1 K232A, has been associated with milk production and composition in taurine breeds. However, this polymorphism is not a good tool for ascertaining the effects of this QTL in Bos indicus (Zebu), since the frequency of the DGAT1 232A allele is too low in these breeds. We sequenced the 3′-untranslated region of DGAT1 gene in a sample of bulls of the breeds Guzerá (Bos indicus) and Holstein (Bos taurus) and, using in silico analysis, we searched for genetic variation, evolutionary conservation, regulatory elements, and possible substitution effects. Six single nucleotide (SNPs) and one insertion-deletion (INDEL) polymorphisms were found in the Guzerá bulls. Additionally, we developed a preliminary association study, using this INDEL polymorphism as a genetic marker. A significant association was detected (P ≤ 0.05) between the INDEL (DGAT1 3′UTR INDEL) and the breeding values (BV) for protein, fat, and milk yields over a 305-day lactation period. The DGAT1 3′ UTR INDEL genotype I/I (I, for insertion) was associated with lower BVs (?38.77 kg for milk, ?1.86 kg for fat, and ?1.48 kg for protein yields), when compared to the genotype I/D (D, for deletion). I/D genotype was lower D/D genotype (?34.98 kg milk, ?1.73 kg fat, and ?1.09 kg protein yields). This study reports the first polymorphism of DGAT1 3′UTR in the Guzerá breed, as well as its association with BV for milk protein, fat, and milk yields.     

Genome‐wide analysis of Italian sheep diversity reveals a strong geographic pattern and cryptic relationships between breeds

Italy counts several sheep breeds, arisen over centuries as a consequence of ancient and recent genetic and demographic events. To finely reconstruct genetic structure and relationships between Italian sheep, 496 subjects from 19 breeds were typed at 50K single nucleotide polymorphism loci. A subset of foreign breeds from the Sheep HapMap dataset was also included in the analyses. Genetic distances (as visualized either in a network or in a multidimensional scaling analysis of identical by state distances) closely reflected geographic proximity between breeds, with a clear north–south gradient, likely because of high levels of past gene flow and admixture all along the peninsula. Sardinian breeds diverged more from other breeds, a probable consequence of the combined effect of ancient sporadic introgression of feral mouflon and long‐lasting genetic isolation from continental sheep populations. The study allowed the detection of previously undocumented episodes of recent introgression (Delle Langhe into the endangered Altamurana breed) as well as signatures of known, or claimed, historical introgression (Merino into Sopravissana and Gentile di Puglia; Bergamasca into Fabrianese, Appenninica and, to a lesser extent, Leccese). Arguments that would question, from a genomic point of view, the current breed classification of Bergamasca and Biellese into two separate breeds are presented. Finally, a role for traditional transhumance practices in shaping the genetic makeup of Alpine sheep breeds is proposed. The study represents the first exhaustive analysis of Italian sheep diversity in an European context, and it bridges the gap in the previous HapMap panel between Western Mediterranean and Swiss breeds.     

Analysis of polymorphism within <Emphasis Type="Italic">POU1F1</Emphasis> gene in relation to milk production traits in dairy Sarda sheep breed

The ovine POU1F1 gene is localized on chromosome 1 and it contains five introns and six exons. In different mammalian species some mutations in different exons are associated with different production traits. The aim of our research was to study the POU1F1 gene nucleotide sequence to detect possible polymorphisms and their relationships with milk productive traits in Sarda breed sheep. The study had been conducted on 140 ewes, 4 or 5 years old coming from a farm located in Sardinia. All the animals were multiparous, lactating and in their third to fifth lactation. Individual milk yield had been recorded monthly and for each sample fat, protein, casein, lactose, and somatic cell count values were analysed. A jugular blood sample was collected from each ewe to perform genomic DNA extraction. PCR, SSCP and sequencing analysis were carried out to examine the six exons to highlight possible SNPs. One-way ANOVA was used to analyse association of variants with milk yield and/or its composition. Two novel SNP were found: 121 C>T in the 5′UTR of the fourth intron fragment and 249 G>A in the 3′UTR of the sixth exon fragment. The statistical analysis did not shown association between milk productive traits and the found polymorphisms. However, further investigations about the promoter region or the prophet genes, like the PROP-1, could clarify its exact role in regulating the productive traits in sheep.     

Characterization and Genetic Study of the Ovine α S2 -Casein (CSN1S2) Allele B

An interesting and quite complex protein pattern has been described at ovine milk proteins but the genetic control of the variation observed was assessed only in few cases. The aim of this work was to characterize the ovine α _s2 -casein (CSN1S2) B variant, first observed in the Italian Gentile di Puglia, a fine-wooled ovine breed, and to investigate its occurrence in two further breeds, the Sarda and Camosciata, which are the most widespread dairy breeds in Italy. The B variant differs from the most common form A with two amino acid exchanges: Asp₇₅ → Tyr₇₅ and Ile₁₀₅ → Val₁₀₅. The first substitution, resulting in a loss of a negative charge, is responsible for the higher isoelectric point of the B protein variant, which allows its detection by isoelectric focusing electrophoresis (IEF). The occurrence of CSN1S2*B in Sarda and Comisana was demonstrated. Since the Asp₇₅ → Tyr₇₅ substitution modifies the protein electric charge, milk properties may result affected to some extent.     

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.     

Milk Fat Content and DGAT1 Genotype Determine Lipid Composition of the Milk Fat Globule Membrane

During secretion of milk fat globules, triacylglycerol (TAG) droplets are enveloped by a phospholipid (PL) trilayer. Globule size has been found to be related to polar lipid composition and fat content, and milk fat content and fatty acid composition have been associated with the diacylglycerol acyltransferase 1 (DGAT1) K232A polymorphism; however, the association between the DGAT1 polymorphism and fat globule size and polar lipid composition has not been studied. The ratio between polar and neutral lipids as well as the composition of the polar lipids in milk has industrial as well as nutritional and health implications. Understanding phenotypic and genotypic factors influencing these parameters could contribute to improving milk lipid composition for dairy products. The focus of the present study was to determine the effect of both fat content and DGAT1 polymorphism on PL/TAG ratio, as a marker for milk fat globule size, and detailed PL composition. Milk samples were selected from 200 cows such that there were equal numbers of samples for the different fat contents as well as per DGAT1 genotype. Samples were analyzed for neutral and polar lipid concentration and composition. PL/TAG ratio was significantly associated with both fat content and DGAT1 genotype. Phosphatidylinositol and phosphatidylserine concentrations were associated with fat content*DGAT1 genotype with a stronger association for the AA than the KK genotype. Sphingomyelin concentration tended to interact with fat content*DGAT1 genotype. Phosphatidylethanolamine (PE) concentration showed a biphasic response to fat content, suggesting that multiple biological processes influence its concentration. These results provide a new direction for controlling polar lipid concentration and composition in milk through selective breeding of cows.     

Novel Polymorphisms in Porcine 3′UTR of the Leptin Gene,Including a Rare Variant within Target Sequence for MIR-9 Gene in Duroc Breed,not Associated with Production Traits

The leptin gene (LEP) is considered as a functional candidate for production traits in livestock due to its crucial role in energy homeostasis. Because polymorphisms in regulatory sequences may affect gene expression, we searched for them in the 3′UTR of LEP and analyzed their association with production traits. Four breeds and a composite line were studied. In the Polish Landrace and Polish Large White breeds, 8 SNPs and 1 indel were observed; whereas, in the Duroc breed, 9 specific SNPs were found. Pietrain and Line 990 were monomorphic. One SNP (g.+168C>T), observed in the Duroc breed only, was located within a target site for microRNA (miR-9). Association studies showed a weak association between one SNP (c.+846C>T) and abdominal fat weight in the Polish Landrace only. Thus, we concluded that contribution of polymorphisms in the 3′UTR to phenotypic variability of pig production traits is marginal. Moreover, we presented an overview of known polymorphisms (128) in the pig leptin gene.     

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.     

11个绵羊品种MSTN基因非翻译区的变异

利用PCR-RFLP技术对特克塞尔羊、夏洛莱羊、小尾寒羊、蒙古羊、乌珠穆沁羊、阿勒泰羊、呼伦贝尔羊、塔什库尔干羊、多浪羊、湖羊和岗巴羊11个品种的345个个体的肌肉生长抑制素(Myostatin, MSTN)基因非翻译区(UTR)的变异进行了多态性分析。结果表明大小为271 bp和1 003 bp的扩增片段经限制性内切酶MboⅡ和BsaⅠ酶切后表现多态, 经卡方检验所有品种在该基因座位均处于平衡状态(P>0.05), 3种基因型在11个绵羊品种中的分布差异极显著(P<0.01)。通过限制性内切酶HpyCH4Ⅳ 酶切实验, 证明我国9个地方绵羊品种不存在特克塞尔绵羊中发现的导致肌肉发达的SNP位点, 并在3′UTR区发现了个别碱基突变位点能够形成miRNA作用的靶基序, 测序表明3′UTR区的突变频率较高。     

Application of Selection Mapping to Identify Genomic Regions Associated with Dairy Production in Sheep

In Europe, especially in Mediterranean areas, the sheep has been traditionally exploited as a dual purpose species, with income from both meat and milk. Modernization of husbandry methods and the establishment of breeding schemes focused on milk production have led to the development of “dairy breeds.” This study investigated selective sweeps specifically related to dairy production in sheep by searching for regions commonly identified in different European dairy breeds. With this aim, genotypes from 44,545 SNP markers covering the sheep autosomes were analysed in both European dairy and non-dairy sheep breeds using two approaches: (i) identification of genomic regions showing extreme genetic differentiation between each dairy breed and a closely related non-dairy breed, and (ii) identification of regions with reduced variation (heterozygosity) in the dairy breeds using two methods. Regions detected in at least two breeds (breed pairs) by the two approaches (genetic differentiation and at least one of the heterozygosity-based analyses) were labeled as core candidate convergence regions and further investigated for candidate genes. Following this approach six regions were detected. For some of them, strong candidate genes have been proposed (e.g. ABCG2, SPP1), whereas some other genes designated as candidates based on their association with sheep and cattle dairy traits (e.g. LALBA, DGAT1A) were not associated with a detectable sweep signal. Few of the identified regions were coincident with QTL previously reported in sheep, although many of them corresponded to orthologous regions in cattle where QTL for dairy traits have been identified. Due to the limited number of QTL studies reported in sheep compared with cattle, the results illustrate the potential value of selection mapping to identify genomic regions associated with dairy traits in sheep.     

Identification of novel SNPs in the Sarda breed goats POU1F1 gene and their association with milk productive performance

The aim of the study was to detect polymorphism in the POU1F1 gene in Sarda breed goat, as well as to establish if SNPs could be associated with milk productive traits. The research was conducted on 129 Sarda breed goats from 4 to 5 years old, multiparous, lactating and in their third to fifth lactation. We report nine exonic and seven non-coding regions SNPs within the Sarda goat POU1F1 gene, namely, Ex 1 61 G>C; Ex 1 108 G>A; Ex 3 C>T; Ex 3 92 C>T; Ex 4 110 A>G; Ex 5 34 G>A resulting in Arg213Lys change; IVS4 641 G>A, IVS4 643 A>C, IVS4 659 G>A, IVS4 677 A>C, IVS4 G699Del, IVS4 709 C>G, Ex 6 17 G>A resulting in Arg228Ser change, Ex 6 58 G>T, Ex 6 172 T>C, 3′UTR 110 T>C. A statistically significant association was found between genotype TT, in position 17 of the exon 6 (3.1 % of frequency), and increased milk yield (P < 0.01) while genotype GT (25.6 % of frequency) was associated with a higher fat content. Genotype TT in position 58 of the exon 6 (3.9 % of frequency) was found to be associated with a higher fat (P < 0.01) and protein content (P < 0.05). Twenty-eight haplotypes were detected, but no significant association between the haplotypes and the milk production traits have been found. Our data, as well as providing new SNPs extending the POU1F1 gene characterization, evidence a relationship between polymorphism and milk production traits in Sarda goat breed.     

Study of polymorphisms in the promoter region of ovine β-lactoglobulin gene and phylogenetic analysis among the Valle del Belice breed and other sheep breeds considered as ancestors

The aim of this work was to sequence the promoter region of β-lactoglobulin (BLG) gene in four sheep breeds, in order to identify polymorphisms, infer and analyze haplotypes, and phylogenetic relationship among the Valle del Belice breed and the other three breeds considered as ancestors. Sequencing analysis and alignment of the obtained sequences showed the presence of 36 single nucleotide polymorphisms (SNPs) and one deletion. A total of 22 haplotypes found in “best” reconstruction were inferred considering the 37 polymorphic sites identified. Haplotypes were used for the reconstruction of a phylogenetic tree using the Neighbor-Joining algorithm. The number of polymorphisms identified showed high variability within breeds. Analysis of genetic diversity indexes showed that the Sarda breed presented the lowest nucleotide diversity, whereas the Comisana breed presented the highest one. Comparing the nucleotide diversity among breeds, the highest value was obtained between Valle del Belice and Pinzirita breeds, whereas the lowest one was between Valle del Belice and Sarda breeds. Considering that polymorphisms in the promoter region of BLG gene could have a functional role associated with milk composition, the lowest value of nucleotide diversity between Valle del Belice and Sarda breeds may be related to a higher similarity of milk composition of these two breeds compared to the others. Further analyses will be conducted in order to evaluate the possible correlation between the genetic diversity indexes and the BLG content in milk of our breeds.     

Influence of DGAT1 K232A Polymorphism on Milk Fat Percentage and Fatty Acid Profiles in Romanian Holstein Cattle

Milk and dairy products are considered the main sources of saturated fatty acids, which are a valuable source of nutrients in the human diet. Fat composition can be adjusted through guided nutrition of dairy animals but also through selective breeding. Recently, a dinucleotide substitution located in the exon 8 of the gene coding for acyl CoA: diacylglycerol acyltransferase 1 (DGAT1), that alters the amino acid sequence from a lysine to an alanine (p.Lys232Ala) in the mature protein, was shown to have a strong effect on milk fat content in some cattle breeds. Therefore, the objectives of this work were to study the occurrence of the DGAT1 p.Lys232Ala polymorphism in Romanian Holstein cattle and Romanian Buffalo breeds and to further investigate its possible influence on fat percentage and fatty acid profiles. The results obtained in this study show that in Romanian Holstein cattle the K allele is associated with increased fat percentage and higher levels of C16:0 and C18:0 fatty acids. The ratio of saturated fatty acids versus unsaturated fatty acids (SFA/UFA) was also higher in KK homozygous individuals, whereas the fractions of C14:0, unsaturated C18 decreased. The DGAT1 p.Lys232Ala polymorphism revealed a high genetic variance for fat percentage, unsaturated C18, C16:0, and SFA/UFA. Although the effect of this polymorphism was not so evident for short chain fatty acids such as C4:0–C8:0, it was significant for C14:0 fatty acids. We concluded that selective breeding of carriers of the A allele in Romanian Holsteins can contribute to improvement in unsaturated fatty acids content of milk. However, in buffalo, the lack of the A allele makes selection inapplicable because only the K allele, associated with higher saturated fatty acids contents in milk, was identified.     

Genome Wide Analysis of Fertility and Production Traits in Italian Holstein Cattle

A genome wide scan was performed on a total of 2093 Italian Holstein proven bulls genotyped with 50K single nucleotide polymorphisms (SNPs), with the objective of identifying loci associated with fertility related traits and to test their effects on milk production traits. The analysis was carried out using estimated breeding values for the aggregate fertility index and for each trait contributing to the index: angularity, calving interval, non-return rate at 56 days, days to first service, and 305 day first parity lactation. In addition, two production traits not included in the aggregate fertility index were analysed: fat yield and protein yield. Analyses were carried out using all SNPs treated separately, further the most significant marker on BTA14 associated to milk quality located in the DGAT1 region was treated as fixed effect. Genome wide association analysis identified 61 significant SNPs and 75 significant marker-trait associations. Eight additional SNP associations were detected when SNP located near DGAT1 was included as a fixed effect. As there were no obvious common SNPs between the traits analyzed independently in this study, a network analysis was carried out to identify unforeseen relationships that may link production and fertility traits.     

Identification of novel SNPs of ovine <Emphasis Type="Italic">PRL</Emphasis> gene and their association with milk production traits

Prolactin (PRL) is a lactogenic hormone that plays a significant role in milk production; its depletion in sheep provokes a severe reduction of milk secretion. Thus, PRL also could be used as a positional marker gene associated with milk production and composition traits. Therefore, the purpose of the study was to identify genotype frequencies of single nucleotide polymorphisms the intron 2 in ovine PRL gene and its possible association genotypes with milk traits in dairy sheep breeds. The genetic structures of ovine PRL gene were examined by PCR-RFLP and DNA sequencing methods in three sheep populations. Four hundred and fifty blood and milk samples were used in the study, which included 150 samples from each of Sakiz, Akkaraman and Awassi ewes respectively. As a result, PRL genotype AA showed a strong association with milk yields content, whereas the animals carrying BB genotype had a higher fat percentage value in the three sheep breeds. Haplotype analysis of the obtained sequences showed the presence of 12 haplotypes in the PRL intron 2 region. In the present study, we have reported for the first time 48 SNPs of the PRL gene for intron 2 in dairy sheep breeds. These preliminary results indicate that the identified SNPs lend themselves readily for further research regarding physiological impacts such as milk production and reproductive traits in other dairy sheep populations.     

Genomic organization of the DGAT2/MOGAT gene family in cattle (Bos taurus) and other mammals

We report the cloning and initial characterization of the genes encoding DGAT2 (diacylglycerol transferase 2), MOGAT1 and MOGAT2 (monoacylglycerol transferases 1 and 2) in domestic cattle (Bos taurus). The three closely related genes belong to a gene family with at least eight members in mammals and are candidate genes for quantitative traits related to dietary fat uptake, lipid synthesis and storage. MOGAT2 and DGAT2 form a tandem and were mapped to bovine chromosome (BTA) 15q25-->q26 by fluorescence in situ hybridization. MOGAT1 was localized to BTA 2q43-->q44. The three genes were investigated for polymorphisms that might be associated with breeding values for milk fat percentage in the dairy breeds German Holstein, German Simmental and German Brown. All the detected polymorphisms were located outside exons or, with one exception, were silent. In MOGAT1, a missense mutation in exon 4 was found that causes a non-conservative substitution of cysteine170 (uncharged, hydrophobic) by lysine (positively charged, hydrophilic). However, allele frequency estimates from pooled DNA samples revealed no significant association of the observed polymorphisms with breeding values for milk fat percentage. A comparative analysis of chromosomal locations and exon-intron structure of the known members of the DGAT2/MOGAT gene family in humans, rodents and cattle indicates an ancient tandem duplication of the ancestor gene combined with an intron gain (or loss) in one copy. Further members of the family may have arisen by duplications of this gene tandem via two rounds of interchromosomal or genome duplications as well as further local (single) gene duplication and loss events.     

DGAT1, a new positional and functional candidate gene for intramuscular fat deposition in cattle

Intramuscular fat content, also assessed as marbling of meat, represents an important beef quality trait. Recent work has mapped a quantitative trait locus (QTL) with an effect on marbling to the centromeric region of bovine chromosome 14, with the gene encoding thyroglobulin (TG) being proposed as a positional and functional candidate gene for this QTL. Recently, the gene encoding diacylglycerol O-acyltransferase (DGAT1), which also has been mapped within the region of the marbling QTL, has been demonstrated to affect the fat content of milk. In the present study, the effects of a 5'-polymorphism of TG and of a lysine/alanine polymorphism of DGAT1 on the fat content of musculus (m.) semitendinosus and m. longissimus dorsi in 55 bovine animals (28 German Holstein and 27 Charolais) has been investigated. Significant effects were found for both candidate genes in both the breeds. These effects seem to be independent of one another because the alleles of the two polymorphisms showed no statistically significant disequilibrium. The DGAT1 effect is mainly on the m. semitendinosus. The TG polymorphism only affects m. longissimus dorsi. However, both intramuscular fat enhancing effects seem to be recessive. The possibility of two linked loci, acting recessively on intramuscular fat content, will require special strategies when selecting for higher marbling scores.