Salvatore Mastrangelo

Genomic inbreeding estimation in small populations: evaluation of runs of homozygosity in three local dairy cattle breeds

In the local breeds with small population size, one of the most important problems is the increase of inbreeding coefficient (F). High levels of inbreeding lead to reduced genetic diversity and inbreeding depression. The availability of high-density single nucleotide polymorphism (SNP) arrays has facilitated the quantification of F by genomic markers in farm animals. Runs of homozygosity (ROH) are contiguous lengths of homozygous genotypes and represent an estimate of the degree of autozygosity at genome-wide level. The current study aims to quantify the genomic F derived from ROH (F ROH) in three local dairy cattle breeds. F ROH values were compared with F estimated from the genomic relationship matrix (F GRM), based on the difference between observed v. expected number of homozygous genotypes (F HOM) and the genomic homozygosity of individual i (F MOL i ). The molecular coancestry coefficient (f MOL ij ) between individuals i and j was also estimated. Individuals of Cinisara (71), Modicana (72) and Reggiana (168) were genotyped with the 50K v2 Illumina BeadChip. Genotypes from 96 animals of Italian Holstein cattle breed were also included in the analysis. We used a definition of ROH as tracts of homozygous genotypes that were >4 Mb. Among breeds, 3661 ROH were identified. Modicana showed the highest mean number of ROH per individual and the highest value of F ROH, whereas Reggiana showed the lowest ones. Differences among breeds existed for the ROH lengths. The individuals of Italian Holstein showed high number of short ROH segments, related to ancient consanguinity. Similar results showed the Reggiana with some extreme animals with segments covering 400 Mb and more of genome. Modicana and Cinisara showed similar results between them with the total length of ROH characterized by the presence of large segments. High correlation was found between F HOM and F ROH ranged from 0.83 in Reggiana to 0.95 in Cinisara and Modicana. The correlations among F ROH and other estimated F coefficients were generally lower ranged from 0.45 (F MOL i -F ROH) in Cinisara to 0.17 (F GRM-F ROH) in Modicana. On the basis of our results, recent inbreeding was observed in local breeds, considering that 16 Mb segments are expected to present inbreeding up to three generations ago. Our results showed the necessity of implementing conservation programs to control the rise of inbreeding and coancestry in the three Italian local dairy cattle breeds.

The genome-wide structure of two economically important indigenous Sicilian cattle breeds.

Genomic technologies, such as high-throughput genotyping based on SNP arrays, provided background information concerning genome structure in domestic animals. The aim of this work was to investigate the genetic structure, the genome-wide estimates of inbreeding, coancestry, effective population size (Ne), and the patterns of linkage disequilibrium (LD) in 2 economically important Sicilian local cattle breeds, Cinisara (CIN) and Modicana (MOD), using the Illumina Bovine SNP50K v2 BeadChip. To understand the genetic relationship and to place both Sicilian breeds in a global context, genotypes from 134 other domesticated bovid breeds were used. Principal component analysis showed that the Sicilian cattle breeds were closer to individuals of Bos taurus taurus from Eurasia and formed nonoverlapping clusters with other breeds. Between the Sicilian cattle breeds, MOD was the most differentiated, whereas the animals belonging to the CIN breed showed a lower value of assignment, the presence of substructure, and genetic links with the MOD breed. The average molecular inbreeding and coancestry coefficients were moderately high, and the current estimates of Ne were low in both breeds. These values indicated a low genetic variability. Considering levels of LD between adjacent markers, the average r(2) in the MOD breed was comparable to those reported for others cattle breeds, whereas CIN showed a lower value. Therefore, these results support the need of more dense SNP arrays for a high-power association mapping and genomic selection efficiency, particularly for the CIN cattle breed. Controlling molecular inbreeding and coancestry would restrict inbreeding depression, the probability of losing beneficial rare alleles, and therefore the risk of extinction. The results generated from this study have important implications for the development of conservation and/or selection breeding programs in these 2 local cattle breeds.

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.

Genome-wide analysis in endangered populations: A case study in Barbaresca sheep

Analysis of genomic data is becoming increasingly common in the livestock industry and the findings have been an invaluable resource for effective management of breeding programs in small and endangered populations. In this paper, with the goal of highlighting the potential of genomic analysis for small and endangered populations, genome-wide levels of linkage disequilibrium, measured as the squared correlation coefficient of allele frequencies at a pair of loci, effective population size, runs of homozygosity (ROH) and genetic diversity parameters, were estimated in Barbaresca sheep using Illumina OvineSNP50K array data. Moreover, the breeds genetic structure and its relationship with other breeds were investigated. Levels of pairwise linkage disequilibrium decreased with increasing distance between single nucleotide polymorphisms. An average correlation coefficient <0.25 was found for markers located up to 50 kb apart. Therefore, these results support the need to use denser single nucleotide polymorphism panels for high power association mapping and genomic selection efficiency in future breeding programs. The estimate of past effective population size ranged from 747 animals 250 generations ago to 28 animals five generations ago, whereas the contemporary effective population size was 25 animals. A total of 637 ROH were identified, most of which were short (67%) and ranged from 1 to 10 Mb. The genetic analyses revealed that the Barbaresca breed tended to display lower variability than other Sicilian breeds. Recent inbreeding was evident, according to the ROH analysis. All the investigated parameters showed a comparatively narrow genetic base and indicated an endangered status for Barbaresca. Multidimensional scaling, model-based clustering, measurement of population differentiation, neighbor networks and haplotype sharing distinguished Barbaresca from other breeds, showed a low level of admixture with the other breeds considered in this study, and indicated clear genetic differences compared with other breeds. Attention should be given to the conservation of Barbaresca due to its critical conservation status. In this context, genomic information may have a crucial role in management of small and endangered populations.

Genetic polymorphism at the CSN1S1 gene in Girgentana dairy goat breed

The aim of this work was to evaluate the variability of the αs1-casein locus in the endangered Girgentana dairy goat breed in order to define genetic improvement and a conservation program for this breed. The study was performed on 200 dairy goats by means of different PCR protocols. The most frequent alleles were A (0.590) and F (0.290) followed by B (0.065) and N (0.047). CSN1S1 E allele was identified with a very low frequency (0.008). The most common genotype was AF (0.365) followed by AA (0.340). The high frequency of the strong genotypes is associated with the production of milk with high fat and protein content and with optimal technological properties. In Girgentana goat breed, the CSN1S1 genotype information could be utilised in selection strategies for milk protein content and milk yield, in order to select genetic lines for the production of ‘drinking milk’ using weak and null genotypes, and for niche products using strong genotypes.

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.

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.

Genetic variability at αs2-casein gene in Girgentana dairy goat breed

Casein genes are highly polymorphic and the high degree of variability has qualitative and quantitative effects on milk composition thereby affecting chemical, physical and technological properties of goat milk. The aim of this work was to evaluate the genetic polymorphisms of the as2-casein (CSN1S2) gene in the endangered Girgentana dairy goat breed in order to assess the genotypes distribution, as it is known genotype influences technological and nutritional milk properties. The study was performed on 207 sample of Girgentana goat breed, analysed with different PCR protocols. The most frequent alleles was A (0.722), followed by F (0.225), C (0.051) and E (0.002) while B, D and 0 alleles were not found. Genotypes detected were AA (0.512), AF (0.338), AC (0.082), FF (0.043), CF (0.020) and EF (0.005). Our results suggested that Girgentana goat breed could be used for the production of milk with high fat and protein content and with optimal technological ability, suitable for cheese making.

Quantitative determination of casein genetic variants in goat milk: Application in Girgentana dairy goat breed.

The study was conducted to develop a high-performance liquid chromatographic (HPLC) method to quantify casein genetic variants (αs2-, β-, and κ-casein) in milk of homozygous individuals of Girgentana goat breed. For calibration experiments, pure genetic variants were extracted from individual milk samples of animals with known genotypes. The described HPLC approach was precise, accurate and highly suitable for quantification of goat casein genetic variants of homozygous individuals. The amount of each casein per allele was: αs2-casein A = 2.9 ± 0.8 g/L and F = 1.8 ± 0.4 g/L; β-casein C = 3.0 ± 0.8 g/L and C1 = 2.0 ± 0.7 g/L and κ-casein A = 1.6 ± 0.3 g/L and B = 1.1 ± 0.2 g/L. A good correlation was found between the quantities of αs2-casein genetic variants A and F, and β-casein C and C1 with other previously described method. The main important result was obtained for κ-casein because, till now, no data were available on quantification of single genetic variants for this protein.