Lina Tortorici

Genome wide linkage disequilibrium and genetic structure in Sicilian dairy sheep breeds

BackgroundThe recent availability of sheep genome-wide SNP panels allows providing background information concerning genome structure in domestic animals. The aim of this work was to investigate the patterns of linkage disequilibrium (LD), the genetic diversity and population structure in Valle del Belice, Comisana, and Pinzirita dairy sheep breeds using the Illumina Ovine SNP50K Genotyping array.ResultsAverage r2 between adjacent SNPs across all chromosomes was 0.155 ± 0.204 for Valle del Belice, 0.156 ± 0.208 for Comisana, and 0.128 ± 0.188 for Pinzirita breeds, and some variations in LD value across chromosomes were observed, in particular for Valle del Belice and Comisana breeds. Average values of r2 estimated for all pairwise combinations of SNPs pooled over all autosomes were 0.058 ± 0.023 for Valle del Belice, 0.056 ± 0.021 for Comisana, and 0.037 ± 0.017 for Pinzirita breeds. The LD declined as a function of distance and average r2 was lower than the values observed in other sheep breeds. Consistency of results among the several used approaches (Principal component analysis, Bayesian clustering, FST, Neighbor networks) showed that while Valle del Belice and Pinzirita breeds formed a unique cluster, Comisana breed showed the presence of substructure. In Valle del Belice breed, the high level of genetic differentiation within breed, the heterogeneous cluster in Admixture analysis, but at the same time the highest inbreeding coefficient, suggested that the breed had a wide genetic base with inbred individuals belonging to the same flock. The Sicilian breeds were characterized by low genetic differentiation and high level of admixture. Pinzirita breed displayed the highest genetic diversity (He, Ne) whereas the lowest value was found in Valle del Belice breed.ConclusionsThis study has reported for the first time estimates of LD and genetic diversity from a genome-wide perspective in Sicilian dairy sheep breeds. Our results indicate that breeds formed non-overlapping clusters and are clearly separated populations and that Comisana sheep breed does not constitute a homogenous population. The information generated from this study has important implications for the design and applications of association studies as well as for development of conservation and/or selection breeding programs.

Genetic Characterisation of CSN2 Gene in Girgentana Goat Breed

Among calcium sensitive caseins, β-casein is the most abundant in goat milk, representing up to 50% of total casein content. The goat β-casein locus has been widely investigated and at least ten alleles have been identified in different goat breeds. The aim of this work was to investigate the polymorphisms of β-casein gene in Girgentana dairy goat breed in order to assess the genotype distribution and evaluate how frequencies have changed during the last 10 years, as genotype is known to influence technological and nutritional milk properties. Sequencing analysis and alignment of the obtained sequences of β-casein exon 7, showed the presence of C, C1, and A strong alleles, and 0′ null allele, with frequencies of 0.597, 0.326, 0.023, and 0.054, respectively. Seven genotypic classes were found in Girgentana goat breed and the most frequent genotype was CC1 (0.423) followed by CC (0.326), C1C1 (0.110), and C0′ (0.096). No AA nor 0′0′ homozygous individuals were found. The presence of strong alleles at CSN2 gene in Girgentana goat breed could be useful for the production of milk with high protein content and good cheese-making properties. Moreover, food business operators should consider the possibility of reviving interest in Girgentana goat milk using weak and null genotypes at CSN2 locus to make peculiar food products, such as drinking milk.

Molecular Characterisation of κ–Casein Gene in Girgentana Dairy Goat Breed and Identification of Two New Alleles

The κ-casein fraction plays an important role in the formation, stabilisation and aggregation on casein micelles and thus affects technological and nutritional properties of milk. In this study, exon 4 of κ-casein (CSN3) gene was sequenced and analysed in Girgentana goat breed. Analyses of the obtained sequences showed the presence of A, B, D, and G known alleles and two new genetic variants, named D’ and N. The new D’ allele differs from D in one transition, G284→A284, which did not cause amino acid change. The new N allele differs from A in five single nucleotide polymorphisms (SNPs): T245/C245, G284/A284, G309/A309, G471/A471 and T591/C591, while it differs from C in one transition, i.e. T583→C583. Comparing the amino acid sequences of N and A alleles, the first two SNPs caused no amino acid change, whereas the other SNPs produced changes (Val65/Ile65, Val119/Ile119, and Ser159/Pro159, respectively). Comparison of N allele with C revealed the amino acid change Val156→Ala156. The most frequent allele was A (0.480) followed by B (0.363), D (0.112), and N (0.034). The D’ and G alleles were identified only in two animals and in heterozygous conditions with a very low frequency (0.005). The most common genotype was AB (39.5%) followed by AA (19.5%), AD (12.7%), and BB (11.7%). Homozygous D’D’, GG, and NN individuals were not found. Further analysis will be performed in order to establish associations among genotypes and quantitative and qualitative milk traits.

Population genetic structure and milk production traits in Girgentana goat breed

The aim of this work was to evaluate the genetic status of the Girgentana goat, an endangered breed from Sicily (Italy), using microsatellite markers. Furthermore, as the main purpose of the Girgentana breed is milk production, quantitative milk traits were investigated, including fatty acid profile. Molecular data from CSN1S1, CSN2, CSN1S2, and CSN3 casein genes were also used to infer haplotypes. A total of 264 individuals were collected. Samples of Maltese (n = 41) and Derivata di Siria (n = 33) goat breeds were also used to understand the genetic relationship among breeds. Test-day records for milk production were collected to determine daily milk yield, fat, protein, casein, lactose, and somatic cell count. Individual milk samples were also collected for fatty acid extraction. Wright’s statistics, gene flow, Nei genetic distance, factorial correspondence analysis, and Bayesian assignment test showed the existence of genetic variability and differentiation among breeds. The AMOVA results indicated that 89.96% of the total variance was partitioned within populations. The Girgentana breed appears to have a subdivided population, and has not experienced a recent bottleneck. A high variability in milk yield was observed. Mean morning milk yield was 1448 ± 404 g, with 4.30 ± 0.87% and 3.72 ± 0.44% of fat and protein percentages, respectively. The average somatic cell count found in Girgentana goat milk was higher than the threshold of 1 500 000 cells/mL advised in Europe for fresh milk. Gross milk and fatty acid composition were similar to that reported in the literature for other local goat breeds.

Full-length sequencing and identification of novel polymorphisms in the ACACA gene of Valle del Belice sheep breed

The essential role of the acetyl-CoA carboxylase (ACACA) enzyme in milk fatty acid (FA) synthesis suggests that it may be responsible for the phenotypic variability observed in milk. Before attempting association analyses between this gene and/or enzyme and phenotypic traits, a study on the genetic variability within this locus is required. The aim of this work was to sequence the entire coding region of ACACA gene in Valle del Belice sheep breed to identify polymorphic sites. A total of 51 coding exons of ACACA gene were sequenced in 32 individuals of Valle del Belice sheep breed. Sequencing analysis and alignment of obtained sequences showed the presence of 23 polymorphic sites. The most polymorphic was exon 53 which showed presence of 12 single-nucleotide polymorphisms (SNPs), of which eight were missense mutations, caused amino acid changes and therefore may affect protein function or stability causing variation in phenotype. The identified polymorphisms showed high variability of the ACACA gene. Sequences analysis allowed to find six new SNPs in exon 53 (6832C>T; 6835C>A; 6840G>A; 6847G>T; 6852C>T and 6860G>C). A total of 31 haplotypes were inferred. Although this study could not provide association study with production traits, it shows finding of novel SNPs that might be important in future studies and laid the basis for further association analyses needed to evaluate the potential use of these SNPs as genetic markers for fat content and FAs composition in milk of Valle del Belice sheep breed.