Genome-wide association study identifies quantitative trait loci regions involved in muscle acidic profile in Large White heavy pigs.

The widespread use of genome-wide association studies resulted in the discovery of genomic regions associated with fatty acid (FA) composition in different porcine tissues, but little information exists about the genes involved in FA composition of meat obtained from heavy pigs selected for the production of Italian dry-cured hams. To this objective, we genotyped with a single nucleotide polymorphism (SNP) panel 795 Italian Large White heavy pigs to identify the markers and genomic regions associated with Semimembranosus muscle FA profile. Heritability estimates for intramuscular fat FA profile were of low-to-moderate magnitude, suggesting that these traits may be improved with genomic selection. On the whole, 45 SNPs were significantly associated with 14 FAs, and 4 of them (ALGA008109, ALGA0081097, CASI0010164 and SIRI0000267) were associated with more than 1 FA. The palmitoleic : palmitic and oleic : stearic ratios displayed the highest number of significant markers and the most significant associations (Bonferroni adjusted P < 5.00E-07). Of particular interest, the palmitoleic : palmitic ratio was strongly associated with markers located at 111 to 114 Mb on chromosome 14, in the same chromosomal region where Stearoyl-CoA desaturase Δ9 (SCD) gene is located. Several significant chromosomal regions were found; some of them harbour key genes playing pivotal roles in FA desaturation and elongation, such as SCD and some members of the Elongation of Very Long-Chain FA (ELOVL) gene family. The results suggest that the identification of causal mutations in these regions may provide a set of markers useful for selection schemes aimed at improving FA composition in pork products.

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.

Genetic parameters of backfat fatty acids and carcass traits in Large White pigs.

Subcutaneous fat thickness and fatty acid composition (FAC) play an important role on seasoning loss and organoleptic characteristics of seasoned hams. Dry-cured ham industry prefers meats with low contents of polyunsaturated fatty acids (PUFA) because these negatively affect fat firmness and ham quality, whereas consumers require higher contents in those fatty acids (FA) for their positive effect on human health. A population of 950 Italian Large White pigs from the Italian National Sib Test Selection Programme was investigated with the aim to estimate heritabilities, genetic and phenotypic correlations of backfat FAC, Semimembranosus muscle intramuscular fat (IMF) content and other carcass traits. The pigs were reared in controlled environmental condition at the same central testing station and were slaughtered at reaching 150 kg live weight. Backfat samples were collected to analyze FAC by gas chromatography. Carcass traits showed heritability levels from 0.087 for estimated carcass lean percentage to 0.361 for hot carcass weight. Heritability values of FA classes were low-to-moderate, all in the range 0.245 for n-3 PUFA to 0.264 for monounsaturated FA (MUFA). Polyunsaturated fatty acids showed a significant genetic correlation with loin thickness (0.128), backfat thickness (-0.124 for backfat measured by Fat-O-Meat'er and -0.175 for backfat measured by calibre) and IMF (-0.102). Obviously, C18:2(n-6) shows similar genetic correlations with the same traits (0.211 with loin thickness, -0.206 with backfat measured by Fat-O-Meat'er, -0.291 with backfat measured by calibre and -0.171 with IMF). Monounsaturated FA, except with the backfat measured by calibre (0.068; P<0.01), do not show genetic correlations with carcass characteristics, whereas a negative genetic correlation was found between MUFA and saturated FA (SFA; -0.339; P<0.001). These results suggest that MUFA/SFA ratio could be increased without interfering with carcass traits. The level of genetic correlations between FA and carcass traits should be taken into account in dealing with the development of selection schemes addressed to modify carcass composition and/or backfat FAC.

Genome-wide association study for longevity in the Holstein cattle population.

Longevity is one of the most important traits determining dairy cow profitability. In the last decades dairy cows suffered a lowering in the age at culling. With the aim to identify the genes involved in longevity, dates of birth, yields, dates of calving during lifespan and culling dates were collected for 946 culled cows which had been genotyped with the Bovine High Density panel. Using the GenABEL package in R, genome-wide association analysis was performed on three potential traits of longevity: (1) 'days in production,' (2) 'days in herd,' (3) number of calvings over lifespan.' Five genome-wide significant single nucleotide polymorphisms (SNPs) associated with all three longevity traits were detected. Several consecutive SNPs identified on chromosomes 16 and 30 indicated the presence of two suggestive quantitative trait loci (QTL). The genes comprised in the QTL regions had biological functions related to fertility, reproductive disorders, heat stress and welfare of cows. These findings might contribute to improving breeding strategies to improve longevity.

Genetic parameters and genome-wide associations of twinning rate in a local breed, the Maremmana cattle.

This study seeks to verify the feasibility of increasing twinning in a herd of the Italian autochtonous Maremmana breed. The data set included 1260 individuals born from 1963 to 2014, 527 males and 733 females, 402 of them calving at least once from 1983 through 2015. Breeding values for twinning were estimated by a single-trait linear animal model. However, since twinning is a dichotomous trait and the frequency of twins is far smaller than the frequency of single births, breeding values were also estimated by a single-trait animal threshold model. Heritability of twinning was 0.014±0.018 and 0.062±0.093 for the linear and the threshold models, respectively. Repeatability was 0.071±0.004 and 0.286± 0.012, respectively, for the two models. Genotyping with the Illumina BovineSNP54 BeadChip was performed for cows living on farm in 2012 (119 cows) and a genome-wide association analysis was performed on the corrected phenotype of all calving during the lifespan of each cow, using the GenABEL package in R and a three step GRAMMAR-GC approach. Genomic heritability, calculated from the genomic kinship matrix estimated through genomic marker data, was 0.29±0.021. The most significant detected single nucleotide polymorphisms (Hapmap22923-BTA-129564) was located in proximity of two genes, ARHGAP8 and TMEM200C, which might be potential functional candidates for twinning rate in cattle.

The ACACA gene is a potential candidate gene for fat content in sheep milk.

No major gene has yet been reported in sheep that explains the variation of milk fat content. The coding region of the acetyl-CoA carboxylase alpha (ACACA) gene, which plays an important role in de novo fatty acid synthesis, had been investigated, but no non-synonymous mutations have been reported. In this study, the genomic regions encoding the three promoters of the ACACA gene were directly sequenced in 264 sheep of three different breeds, and 10 SNPs were identified. Allele frequencies of most SNPs significantly differed (P = 0.05-0.0001) between breeds. The SNPs that potentially altered either gene regulatory elements or putative binding sites of transcription factors were made evident through in silico analysis. The association analysis with milk traits, performed for one SNP of PIII (GenBank AJ292286, g.1330G>T), showed a significant allelic substitution effect (+0.33%, P < 0.0001 and +0.35%, P < 0.01) in the Altamurana and Gentile breeds respectively. Because this SNP was located in the binding site of the paired box protein transcription factors, which was shown to function as an efficient promoter element, and because PIII transcripts are expressed in the mammary gland, the SNP in PIII of the ACACA gene might affect the variation of fat content in sheep milk.

Leptin gene haplotypes are associated with change in immunological and hematological variables in dairy cow during the peripartum period.

In this study, the effect of polymorphisms in the leptin gene on the hematological variables in periparturient dairy cows was investigated. The hematological profile of 67 Holstein cows was assessed for 6 wk around calving. The DNA of the cows was genotyped at 6 polymorphic loci within the leptin gene, and 7 haplotypes were reconstructed. Significant haplotype substitution effects were found, for haplotype 1, on total white blood cell count for 2 wk around calving (+0.70 10(3)/µL, P = 0.05; +1.38 10(3)/µL, P = 0.0001); on neutrophil cell count in the first week after calving (+0.94 10(3)/µL, P = 0.001); on lymphocyte count during the 3 wk before and the first week after calving (+0.32 10(3)/µL, P = 0.05; +0.27 10(3)/µL, P = 0.03; +0.26 10(3)/µL, P = 0.04; +0.34 10(3)/µL, P = 0.01); on red blood cell count during the last week before calving and wk 1 and 2 after calving (+0.21 10(6)/µL, P = 0.02; +0.23 10(6)/µL, P = 0.01; +0.20 10(6)/µL, P = 0.03); on mean corpuscular volume (-1.35 fL, P = 0.01; -1.29 fL, P = 0.002; -1.18 fL, P = 0.004; -1.09 fL, P = 0.008; -1.23 fL, P = 0.003; -1.31 fL, P = 0.003); and on mean corpuscular hemoglobin (-0.37 pg, P = 0.05; -0.38 pg, P = 0.02; -0.39 pg, P = 0.01; -0.34 pg, P = 0.03; -0.40 pg, P = 0.01; -0.40 pg, P = 0.01) during all 6 wk analyzed. Significant haplotype substitution effects, but opposite those of haplotype-1, were found for haplotype-2 on white blood cell count (-1.10 10(3)/µL, P = 0.01; -1.30 10(3)/µL, P = 0.002; -1.09 10(3)/µL, P = 0.01) and neutrophil count (-0.82 10(3)/µL, P = 0.02; -0.95 10(3)/µL, P = 0.005; -0.92 10(3)/µL, P = 0.01). Haplotype-3 influenced red blood cell count (-0.23 10(6)/µL, P = 0.03; -0.28 10(6)/µL, P = 0.01; -0.34 10(6)/µL, P = 0.002) during the last 2 wk before and the first week after calving, and also, with effects evident only in wk 3 and 2 before calving, mean corpuscular volume (+1.38 fL, P = 0.03; +0.97 fL, P = 0.05; +1.08 fL, P = 0.05), mean corpuscular hemoglobin (+0.58 pg, P = 0.02; +0.38 pg, P = 0.04; 0.51 pg, P = 0.01), and red blood cell distribution width (-0.56% P = 0.02; -0.47%, P = 0.05). The current study provided evidence that several polymorphisms in the leptin gene play a role in the variability of hematological variables during the peripartum period, and might be used as genetic markers for improving the immunological conditions of dairy cows in critical productive periods.

Exploring polymorphisms and effects of candidate genes on milk fat quality in dairy sheep.

The aim of the present study was to investigate the genetic control of the fatty acid (FA) composition in milk from 3 breeds of sheep: Altamurana, Gentile di Puglia, and Sarda. Single nucleotide polymorphisms within genes, encoding enzymes putatively involved in the synthesis and metabolism of milk fat, were selected for analysis, and the allele substitution effects were determined for 16 genes, which were polymorphic in the 3 sheep breeds, upon the milk fat composition. Four genes (alpha-1-antichymotrypsin-2; diacylglycerol O-acyltransferase homolog-2; propionyl Coenzyme A carboxylase, beta polypeptide; and insulin-like growth factor-I) play a role in the desaturation of stearic FA into polyunsaturated fatty acids. Furthermore, 2 genes (growth hormone receptor and zona pellucida glycoprotein-2) affect the variability of the total fat content in addition to the butyric and stearic FA profile, and the fatty acid synthetase gene has an influence on the medium-chain FA. Milk FA profiles play an important role in dairy sheep farming because they have a large effect on cheese characteristics and also because sheep milk may be marketed as a source of nutraceuticals because it contains higher levels of conjugated linoleic acid than milk from other ruminants. The current study evaluated the global effects of a large number of single nucleotide polymorphisms and haplotypes on traits that are not commonly investigated in sheep but that are potentially very useful for improving milk quality.

Ovine acyl CoA:diacylglycerol acyltransferase 1- molecular characterization, polymorphisms and association with milk traits.

The objective of this work was to characterize the complete coding region of the ovine acylCoA:diacylglycerol acyltransferase 1 (DGAT1) gene of three Italian sheep breeds: Sarda, Altamurana and Gentile di Puglia. Characterization was accomplished by direct sequencing of 8676 bp of the relevant DNA, including introns and partial 5' and 3' untranslated regions (UTRs). We detected five novel SNPs; one SNP (g.5553C>T) is located in intron 2, has similar frequencies in the three breeds and showed a negative association with milk fat content. More interesting is an SNP in the 5' UTR (g.127C>A), the occurrence of which is rare in the higher milk-fat breeds (Altamurana and Gentile di Puglia); it is located in the core sequence of Sp1, a putative binding site of transcription factors. This SNP showed a significant negative association with milk fat content in the Sarda sheep. Because DGAT1 plays a fundamental role in triacylglycerol synthesis, the novel detected SNP in the 5' UTR of the DGAT1 gene might explain, at least partially, the variation of fat content in the milk of Sarda sheep.

Short communication: Effect of stearoyl-coenzyme A desaturase polymorphism on fatty acid composition of milk.

The effect of the stearoyl-CoA desaturase (SCD) gene on milk fatty acid composition was tested. Cows of 3 breeds of northern Italy, Piedmontese, Valdostana, and Jersey, were genotyped at exon 5 of the SCD gene. This has been suggested as a primary candidate gene to change the proportion of saturated vs. unsaturated fatty acids in milk, wherein a single nucleotide polymorphism (C/T) gives rise to a different AA codon. It was possible to ascribe a reduced desaturase activity to the T allele only in the case of caproleic and myristoleic fatty acids. In contrast with the findings of SCD effects on carcass fat, it was not possible to confirm the higher desaturation activity of this single nucleotide polymorphism on long-chain fatty acids, due to the different pathways that originate milk fatty acids of different carbon length; long-chain fatty acids are highly influenced by the complex metabolic events that affect the ingested nutrients during their transfer to milk fat.

Meat molecular traceability: How to choose the best set of microsatellites?

To set up routine assays for molecular meat traceability along the food chain, the availability of a simple and low cost test for the identification of cattle carcasses is required. For this purpose, we evaluated 13 microsatellites for their ability in the identification of animals belonging to four Italian cattle breeds. Here we propose a criterion for a microsatellite-based test with the best reliability when reducing the number of loci to be analysed. The method is based on the observation that in the same loci breeds can show differences in frequencies and number of fixed alleles. This non-uniform distribution of alleles between breeds results in differences in the informative content of the same loci in different breeds. Taking into account these differences, it is possible to perform tests for the allocation of samples to specific animals utilizing a small number of microsatellites. The proposed approach allows cost reduction and ease in performing the analyses.

Single nucleotide polymorphism detection in promoter III of the acetyl-CoA carboxylase-alpha gene in sheep.

Three novel SNPs were identified in the locus OAR292286, encoding the DNA sequence of promoter III of the ovine acetyl-CoA carboxylase-alpha gene, in Italian sheep of four breeds: Gentile di Puglia (25 individuals) and Sopravissana (31) which are triple-purpose local endangered breeds, Comisana (25) which is a local non-selected, non-endangered dairy breed and Sarda (15) which is a popular selected high yielding dairy breed. Variant alleles are: G/T at 1330 bp, C/G at 1338 bp and C/T at 1430 bp. Frequencies of the variant alleles were calculated and chi-squared analysis of the differences in allele frequency between breed pairs was performed. Allele frequencies of the Sarda breed differ significantly from the other considered breeds.

Genetic diversity between Piedmontese, Maremmana, and Podolica cattle breeds.

One hundred twenty animals of three native Italian cattle breeds, Piedmontese, Maremmana, and Podolica, were genetically characterized at 21 microsatellite loci located on 13 chromosomes. Allele numbers ranged from 3 to 19; average gene diversity ranged from 0.206 to 0.878 (average 0.738). The breed that preserved the highest genetic variability was the Podolica, where the chosen markers show the highest gene diversity (0.741) and the highest heterozygosity (0.155). The lowest inbreeding rate (0.102) was registered by the Piedmontese. Genetic distances were 0.069 (Piedmontese versus Maremmana), 0.050 (Piedmontese versus Podolica), and 0.041 (Maremmana versus Podolica) and reflect the different phylogenetic origins of the breeds: Maremmana and Podolica originated from the Grey Steppe cattle group, while Piedmontese belongs to the Northern Italy Lowland cattle group. Observed heterozygosity was not significantly different from expected in any of the breeds, which is an indication that they maintain a random mating structure. The probabilistic assignment of all sampled individuals to three theoretical populations, on the basis of allele frequencies, indicated that 82% of Piedmontese, 66% of Maremmana, and 33% of Podolica can be assigned to the appropriate breed with a probability higher than 90%. This result very well reflects the impact of the selection activity on the breed genetic structure. The chosen microsatellites proved to be a good tool for describing the correct reality of the analyzed populations, but they are not sufficient to discriminate between breeds.

Effects of age and calving season on lactation curves of milk production traits in Italian water buffaloes.

Test day (TD) records of milk production traits (milk yield, fat, and protein percentages) of 534 Italian buffalo cows were analyzed with a mixed linear model in order to estimate lactation curves pertaining to different ages at calving and different seasons of calving. Milk yield lactation curves of younger animals were lower than those of older animals until 20 wk from parturition. No effect of age at calving could be observed for fat and protein percentages. Season of calving affected milk yield only in the first phase of lactation, with the lowest production levels for summer calvings; no effect could be observed on fat and protein contents. Average correlations among TD measures within lactation were 0.59, 0.31, and 0.36 for milk yield, fat, and protein percentages, respectively. Five standard linear functions of time were able to reconstruct the average lactation curves. Goodness of fit was satisfactory for all models considered, although only the five-parameter model was flexible enough to fit all the three traits considered with excellent results.