Analysis of ddRAD-seq data provides new insights into the genomic structure and patterns of diversity in Italian donkey populations.

With more than 150 recognized breeds, donkeys assume relevant economic importance, especially in developing countries. Even if the estimated number of heads worldwide is 53M, this species received less attention than other livestock species. Italy has traditionally been considered one of the cradles of European donkey breeding, and despite a considerable loss of biodiversity, today still counts nine autochthonous populations. A total of 220 animals belonging to nine different populations were genotyped using the double-digest restriction site associated DNA (ddRAD) sequencing to investigate the pattern of diversity using a multi-technique approach. A total of 418,602,730 reads were generated and successfully demultiplexed to obtain a medium-density SNP genotypes panel with about 27K markers. The diversity indices showed moderate levels of variability. The genetic distances and relationships, largely agree with the breeding history of the donkey populations under investigation. The results highlighted the separation of populations based on their genetic origin or geographical proximity between breeding areas, showed low to moderate levels of admixture, and indicated a clear genetic difference in some cases. For some breeds, the results also validate the success of proper management conservation plans. Identified runs of homozygosity islands, mapped within genomic regions related to immune response and local adaptation, are consistent with the characteristics of the species known for its rusticity and adaptability. This study is the first exhaustive genome-wide analysis of the diversity of Italian donkey populations. The results emphasized the high informativeness of genome-wide markers retrieved through the ddRAD approach. The findings take on great significance in designing and implementing conservation strategies. Standardized genotype arrays for donkey species would make it possible to combine worldwide datasets to provide further insights into the evolution of the genomic structure and origin of this important genetic resource.

Donkeys assume relevant economic importance in several countries worldwide. However, the genetic structure of these populations is less investigated compared to other species. The aim of this study was to investigate the genetic background of nine different Italian donkey populations. A total of 220 animals were genotyped with about 27K markers extracted by the double-digest restriction site associated DNA sequencing. The consistency of the results across different approaches agreed with the demographic history, the origin, and previous results on the nine donkey populations, suggesting that our conclusions are robust. Moreover, the results of the present study highlighted low to moderate levels of admixture and, for some breeds, confirmed the success of proper management conservation plans.

Effect of genotyping density on the detection of runs of homozygosity and heterozygosity in cattle.

Runs of homozygosity (ROHom) are contiguous stretches of homozygous regions of the genome. In contrast, runs of heterozygosity (ROHet) are heterozygosity-rich regions. The detection of these two types of genomic regions (ROHom and ROHet) is influenced by the parameters involved in their identification and the number of available single-nucleotide polymorphisms (SNPs). The present study aimed to test the effect of chip density in detecting ROHom and ROHet in the Italian Simmental cattle breed. A sample of 897 animals were genotyped at low density (50k SNP; 397 individuals), medium density (140k SNP; 348 individuals), or high density (800k SNP; 152 individuals). The number of ROHom and ROHet per animal (nROHom and nROHet, respectively) and their average length were calculated. ROHom or ROHet shared by more than one animal and the number of times a particular SNP was inside a run were also computed (SNPROHom and SNPROHet). As the chip density increased, the nROHom increased, whereas their average length decreased. In contrast, the nROHet decreased and the average length increased as the chip density increased. The most repeated ROHom harbored no genes, whereas in the most repeated ROHet four genes (SNRPN, SNURF, UBE3A, and ATP10A) previously associated with reproductive traits were found. Across the 3 datasets, 31 SNP, located on Bos taurus autosome (BTA) 6, and 37 SNP (located on BTA21) exceeded the 99th percentile in the distribution of the SNPROHom and SNPROHet, respectively. The genomic region on BTA6 mapped the SLIT2, PACRGL, and KCNIP4 genes, whereas 19 and 18 genes were mapped on BTA16 and BTA21, respectively. Interestingly, most of genes found through the ROHet analysis were previously reported to be related to health, reproduction, and fitness traits. The results of the present study confirm that the detection of ROHom is more reliable when the chip density increases, whereas the ROHet trend seems to be the opposite. Genes and quantitative trait loci (QTL) mapped in the highlighted regions confirm that ROHet can be due to balancing selection, thus related to fitness traits, health, and reproduction, whereas ROHom are mainly involved in production traits. The results of the present study strengthened the usefulness of these parameters in analyzing the genomes of livestock and their biological meaning.

Runs of homozygosity (ROHom), continuous stretches of homozygous loci, and runs of heterozygosity (ROHet), continuous stretches of heterozygous loci, may be due to directional (ROHom) or balancing selection (ROHet) and are interesting to analyze those shared among animals within a population and the genes they harbor. The detection of both types of genomic regions is influenced by genotyping density and involved parameters. Thus, this work aimed to study the impact of the BeadChip density on the ROHom and ROHet detection in the Italian Simmental cattle breed. Results showed that the ROHom detection is more reliable as the density increases, whereas a more cryptic pattern was observed for ROHet. Interestingly, the hypothesis on how these two types of runs arise was supplied by the results of this study. The genes mapped on the highlighted ROHet were mainly associated with fitness traits, health, and reproduction, whereas those found in the ROHom were associated with production traits.

Characterization of heterozygosity-rich regions in Italian and worldwide goat breeds.

Heterozygosity-rich regions (HRR) are genomic regions of high heterozygosity, which may harbor loci related to key functional traits such as immune response, survival rate, fertility, and other fitness traits. This study considered 30 Italian and 19 worldwide goat breeds genotyped with the Illumina GoatSNP50k BeadChip. The aim of the work was to study inter-breed relationships and HRR patterns using Sliding Window (SW) and Consecutive Runs (CR) detection methods. Genetic relationships highlighted a clear separation between non-European and European breeds, as well as the north-south geographic cline within the latter. The Pearson correlation coefficients between the descriptive HRR parameters obtained with the SW and CR methods were higher than 0.9. A total of 166 HRR islands were detected. CHI1, CHI11, CHI12 and CHI18 were the chromosomes harboring the highest number of HRR islands. The genes annotated in the islands were linked to various factors such as productive, reproductive, immune, and environmental adaptation mechanisms. Notably, the Montecristo feral goat showed the highest number of HRR islands despite the high level of inbreeding, underlining potential balancing selection events characterizing its evolutionary history. Identifying a species-specific HRR pattern could provide a clearer view of the mechanisms regulating the genome modelling following anthropogenic selection combined with environmental interaction.

Genome-wide population structure, homozygosity, and heterozygosity patterns of Nero Siciliano pig in the framework of Italian and cosmopolitan breeds.

Analysis of genomic data is becoming more and more common for the effective management of livestock breeding programmes, even in the case of local populations. In this work, the genome-wide data of Nero Siciliano pig breed were compared to that of wild boar, Italian local and cosmopolitan breeds to investigate its genetic structure, and runs of homozygosity (ROH) and heterozygosity patterns. The Nero Siciliano has been reported to have the highest rate of genetic diversity among the Italian breeds, and a genetic variability comparable to that of the cosmopolitan breeds. Analyses of genomic structure and relationships underlined its proximity to wild boar, and an internal substructure probably linked to different family lines. The breed showed a low value of inbreeding estimated from ROH, and the highest diversity index among the Italian breeds, even if lower than that of the cosmopolitans. Four ROH islands in three chromosomes (SSC8, SSC11, and SSC14) and one heterozygosity-rich region (SSC1) were identified in Nero Siciliano, highlighting genomic regions related to productive QTL. Across breeds, SSC8 and SSC14 were the chromosomes with most ROH islands, with Mora Romagnola and wild boar showing the highest level of autozygosity. Chromosomes SSC2, SSC6, SSC8 and SSC13 showed the majority of runs of heterozygosity regions, mainly found in the cosmopolitan pig breeds, which reported several genes associated with health-related QTL. The outlined results can help to better identify the genomic profile of this local breed in order to plan matings, maintain adequate internal diversity and exploit the production system.

A genome-wide comparison between selected and unselected Valle del Belice sheep reveals differences in population structure and footprints of recent selection.

About three decades of breeding and selection in the Valle del Belìce sheep are expected to have left several genomic footprints related to milk production traits. In this study, we have assembled a dataset with 451 individuals of the Valle del Belìce sheep breed: 184 animals that underwent directional selection for milk production and 267 unselected animals, genotyped for 40,660 single-nucleotide polymorphisms (SNPs). Three different statistical approaches, both within (iHS and ROH) and between (Rsb) groups, were used to identify genomic regions potentially under selection. Population structure analyses separated all individuals according to their belonging to the two groups. A total of four genomic regions on two chromosomes were jointly identified by at least two statistical approaches. Several candidate genes for milk production were identified, corroborating the polygenic nature of this trait and which may provide clues to potential new selection targets. We also found candidate genes for growth and reproductive traits. Overall, the identified genes may explain the effect of selection to improve the performances related to milk production traits in the breed. Further studies using high-density array data, would be particularly relevant to refine and validate these results.

High-density SNP markers elucidate the genetic divergence and population structure of Noticiana sheep breed in the Mediterranean context.

Among livestock species, sheep have played an early major role in the Mediterranean area. Italy has a long history of sheep breeding and, despite a dramatic contraction in numbers, still raise several local populations that may represent a unique source of genetic diversity. The Noticiana is a breed of the south-eastern part of Sicily appreciated both for its dairy products and for its resistance to harsh environment. In this study, the high-density Illumina Ovine SNP600K BeadChip array was used for the first genome-wide characterization of 48 individuals of Noticiana sheep to investigate its diversity, the genome structure and the relationship within the context of worldwide and Italian breeds. Moreover, the runs of homozygosity (ROH) pattern and the pairwise FST-outliers were examined. Noticiana reported moderate levels of genetic diversity. The high percentage of short and medium length ROH segments (93% under 4 Mb) is indicative of a within breed relatedness dating back to ancient times, despite the absence of management for the mating plans and the reduced population size. In the worldwide context, the Southern Italian, Spanish and Albanian breeds overlapped in a macro cluster which also included the Noticiana sheep. The results highlighted ancestral genetic components of Noticiana shared with Comisana breed, and showed the clear separation from the other Italian sheep. This is likely the consequence of the combined effects of genetic drift, small population size and reproductive isolation. ROH islands and FST-outliers approaches in Noticiana identified genes and QTLs involved in milk and meat production, as well as related to the local adaptation, and therefore are consistent with the phenotypic traits of the studied breed. Although a wider sampling could be useful to deepen the genomic survey on Noticiana, these results represent a crucial starting point for the characterization of an important local genetic resource, with a view of supporting the local economy and preserving the biodiversity of the sheep species.

High-density single nucleotide polymorphism markers reveal the population structure of 2 local chicken genetic resources.

Italy counts a large number of local chicken populations, some without a recognized genetic structure, such as Val Platani (VPL) and Cornuta (COS), which represent noteworthy local genetic resources. In this study, the genotype data of 34 COS and 42 VPL, obtained with the Affymetrix Axiom600KChicken Genotyping Array, were used with the aim to investigate the genetic diversity, the runs of homozygosity (ROH) pattern, as well as the population structure and relationship within the framework of other local Italian and commercial chickens. The genetic diversity indices, estimated using different approaches, displayed moderate levels of genetic diversity in both populations. The identified ROH hotspots harbored genes related to immune response and adaptation to local hot temperatures. The results on genetic relationship and population structure reported a clear clustering of the populations according to their geographic origin. The COS formed a nonoverlapping genomic cluster and clearly separated from the other populations, but showed evident proximity to the Siciliana breed (SIC). The VPL highlighted intermediate relationships between the COS-SIC group and the rest of the sample, but closer to the other Italian local chickens. Moreover, VPL showed a complex genomic structure, highlighting the presence of 2 subpopulations that match with the different source of the samples. The results obtained from the survey on genetic differentiation underline the hypothesis that Cornuta is a population with a defined genetic structure. The substructure that characterizes the Val Platani chicken is probably the consequence of the combined effects of genetic drift, small population size, reproductive isolation, and inbreeding. These findings contribute to the understanding of genetic diversity and population structure, and represent a starting point for designing programs to monitor and safeguard these local genetic resources, in order to define a possible official recognition program as breeds.

Genomic Tools for the Characterization of Local Animal Genetic Resources: Application in Mascaruna Goat.

Italy contains a large number of local goat populations, some of which do not have a recognized genetic structure. The "Mascaruna" is a goat population reared for milk production in Sicily. In this study, a total of 72 individuals were genotyped with the Illumina Goat_IGGC_65K_v2 BeadChip with the aim to characterize the genetic diversity, population structure and relatedness with another 31 Italian goat populations. The results displayed a moderate level of genetic variability for Mascaruna, in concordance with the estimated values for Italian goats. Runs of homozygosity islands are linked to genes involved in milk production, immune response and local adaptation. Population structure analyses separated Mascaruna from the other goat populations, indicating a clear genetic differentiation. Although they are not conclusive, our current results represent a starting point for the creation of monitoring and conservation plans. Additional analyses and a wider sampling would contribute to refine and validate these results. Finally, our study describing the diversity and structure of Mascaruna confirms the usefulness of applied genomic analyses as valid tools for the study of the local uncharacterized genetic resources.

How Geography and Climate Shaped the Genomic Diversity of Italian Local Cattle and Sheep Breeds.

Understanding the relationships among geography, climate, and genetics is increasingly important for animal farming and breeding. In this study, we examine these inter-relationships in the context of local cattle and sheep breeds distributed along the Italian territory. To this aim, we used redundancy analysis on genomic data from previous projects combined with geographical coordinates and corresponding climatic data. The effect of geographic factors (latitude and longitude) was more important in sheep (26.4%) than that in cattle (13.8%). Once geography had been partialled out of analysis, 10.1% of cattle genomic diversity and 13.3% of that of sheep could be ascribed to climatic effects. Stronger geographic effects in sheep can be related to a combination of higher pre-domestication genetic variability together with biological and productive specificities. Climate alone seems to have had less impact on current genetic diversity in both species, even if climate and geography are greatly confounded. Results confirm that both species are the result of complex evolutionary histories triggered by interactions between human needs and environmental conditions.

A 20-SNP Panel as a Tool for Genetic Authentication and Traceability of Pig Breeds.

Food authentication in local breeds has important implications from both an economic and a qualitative point of view. Meat products from autochthonous breeds are of premium value, but can easily incur fraudulent or accidental substitution or mislabeling. The aim of this study was to identify a small number of SNPs using the Illumina PorcineSNP60 BeadChip for breed traceability, in particular of the Italian Nero Siciliano pig and its derived products. A panel of 12 SNPs was sufficient to discriminate Nero Siciliano pig from cosmopolitan breeds and wild boars. After adding 8 SNPs, the final panel of 20 SNPs allowed us to discriminate all the breeds involved in the study, to correctly assign each individual to its breed, and, moreover, to discriminate Nero Siciliano from first-generation hybrids. Almost all livestock breeds are being genotyped with medium- or high-density SNP panels, providing a large amount of information for many applications. Here, we proposed a method to select a reduced SNP panel to be used for the traceability of pig breeds.

Genome-wide survey on three local horse populations with a focus on runs of homozygosity pattern.

Purosangue Orientale Siciliano, Sanfratellano and Siciliano represent the Sicilian equine genetic resource. This study aimed to investigate the genetic diversity, population structure and the pattern of autozygosity of Sicilian horse populations using genome-wide single-nucleotide polymorphism (SNP) data generated with the Illumina Equine SNP70 array. The genotyping data of 17 European and Middle East populations were also included in the study. The patterns of genetic differentiation, model-based clustering and Neighbour-Net showed the expected positioning of Sicilian populations within the wide analysed framework and the close connections between the Purosangue Orientale Siciliano and the Arab as well as between Sanfratellano, Siciliano and Maremmano. The highest expected heterozygosity (He ) and contemporary effective population size (cNe) were reported in Siciliano (He  = 0.323, cNe = 397), and the lowest were reported in Purosangue Orientale Siciliano (He  = 0.277, cNe = 10). The analysis of the runs of homozygosity and the relative derived inbreeding revealed high internal homogeneity in Purosangue Orientale Siciliano and Arab horses, intermediate values in Maremmano and Sanfratellano and high heterogeneity in the Siciliano population. The genome-wide SNP analysis showed the selective pressure on Purosangue Orientale Siciliano towards traits related to endurance performance. Our results underline the importance of planning adequate conservation and exploitation programmes to reduce the level of inbreeding and, therefore, the loss of genetic diversity.

Feeding System Resizes the Effects of DGAT1 Polymorphism on Milk Traits and Fatty Acids Composition in Modicana Cows.

The interaction between genetic polymorphism and feeding system on milk traits and fatty acid composition was investigated in Modicana cows. Two DGAT1 K232A genotypes (AK and AA) and two feeding regimes, extensive system (EX) with 8 h of grazing without concentrate (EX) and semi-intensive systems (SI) with 2 h of grazing with concentrate, were investigated. DGAT1 genotype did not influence milk yield and composition. The feeding system affected milk composition: protein was significantly higher in SI and lactose in the EX system. A significant genotype × feeding system interaction was observed: the protein and casein levels of AK cows were higher in the SI compared to the EX system. Milk fatty acids profile, total saturated to total unsaturated fatty acids, n-6 to n-3 ratios, and atherogenic index were affected by the feeding system, improving the healthy properties of milk from animals reared in the extensive system. DGAT1 genotype influenced the fatty acid composition: milk from AA cows had a more favorable fatty acid composition due to lower total saturated fatty acids, saturated to unsaturated ratio, atherogenic index, and higher levels of oleic acid and total unsaturated fatty acids. Furthermore, an interaction genotype x feeding system was observed: the AK milk was richer in short-chain FAs (C4:0-C8:0) and C10:0 only in the EX but not in the SI system. Our data suggest that a high amount of green forage in the diet of Modicana cows can resize the effect of the DGAT1 genotype on milk traits and fatty acids composition.

Genome-Wide SNP Analysis Reveals the Population Structure and the Conservation Status of 23 Italian Chicken Breeds.

The genomic variability of local Italian chicken breeds, which were monitored under a conservation plan, was studied using single nucleotide polymorphisms (SNPs) to understand their genetic diversity and population structure. A total of 582 samples from 23 local breeds and four commercial stocks were genotyped using the Affymetrix 600 K Chicken SNP Array. In general, the levels of genetic diversity, investigated through different approaches, were lowest in the local chicken breeds compared to those in the commercial stocks. The level of genomic inbreeding, based on runs of homozygosity (FROH), was markedly different among the breeds and ranged from 0.121 (Valdarnese) to 0.607 (Siciliana). In all breeds, short runs of homozygosity (ROH) (<4 Mb in length) were more frequent than long segments. The patterns of genetic differentiation, model-based clustering, and neighbor networks showed that most breeds formed non-overlapping clusters and were clearly separate populations, which indicated the presence of gene flow, especially among breeds that originated from the same geographical area. Four genomic regions were identified as hotspots of autozygosity (islands) among the breeds, where the candidate genes are involved in morphological traits, such as body weight and feed conversion ratio. We conclude that the investigated breeds have conserved authentic genetic patterns, and these results can improve conservation strategies; moreover, the conservation of local breeds may play an important role in the local economy as a source of high-quality products for consumers.

New Insights into Chemical and Mineral Composition of Donkey Milk throughout Nine Months of Lactation.

Donkey milk is increasingly being proposed as a natural alternative milk for various categories of consumers, especially infants and the elderly population. However, its potential production, gross, and mineral composition have not been deeply investigated yet. Sixty-two individual milk samples were collected monthly from nine Ragusano donkeys reared in a specialized dairy farm. Milk yield as well as chemical and mineral composition, including macro and micro elements, were investigated over an entire lactation, from the second to the ninth month of milking. Milk yield averaged 1.64 kg/day, which highlights good aptitude of the Ragusano breed for the production of milk. Gross composition was characterized by low content of dry matter (8.19%), a high amount of lactose (6.07%), low protein (1.34%), and very low-fat content (0.16%). Whey proteins represented 58% of the total protein, and proteose-peptones accounted for 0.35 ± 0.07 g per 100 g. Total ash content was 0.36 g per 100 g and represented 4.40% of the dry matter. The most abundant element was K, which was followed by Ca, Na, and P. As expected, the micro elements Fe, Zn, and Cu were found in low amounts or in traces. Dry matter, fat, whey proteins. The total ash, Ca, P, Mg, and mineral ratios were significantly affected by the lactation stage.

Casein haplotype diversity in seven dairy goat breeds.

Selection, drift, gene flow and breeding have extensively shaped the genomic variability of domestic animals. In goat species, several mutations identified within the casein genes have been shown to affect the level of gene expression of milk production traits. The four casein genes - CSN1S1, CSN2, CSN1S2 and CSN3 - are organized in a cluster of 250 kb located in chromosome 6, and due to tight linkage, their genetic variability is well depicted by haplotypes which are transmitted to the progeny. Thirty single nucleotide polymorphisms (SNPs) located within the casein gene cluster were used to characterize the haplotype variability of six southern Italian goat breeds (Girgentana, Maltese, Rossa Mediterranea, Argentata dell'Etna, Messinese, Capra dell'Aspromonte). A representative sample of the Norwegian dairy goat breed (Norsk melkegeit) has been used as an out-group to obtain a weighted measure of genetic diversity in the metapopulation. A total of 54 haplotypes were detected among the seven breeds: 26, 9, 8 and 11 haplotypes were found at CSN1S1, CSN2, CSN1S2 and CSN3 respectively. The number of haplotypes per breed was 14 (Norwegian), 26 (Messinese), 27 (Rossa Mediterranea and Girgentana) and 31 (Maltese, Argentata dell'Etna and Capra dell'Aspromonte). The Maltese breed showed the highest number of private haplotypes, whereas the Norwegian goat recorded the highest number of shared haplotypes. The linkage disequilibrium analysis showed higher levels of association for the SNP pairs within casein loci than SNP pairs between casein loci, likely reflecting low levels of intra-genic recombination. The highest linkage disequilibrium values were found in CSN1S1 and CSN2 genes in all the breeds, except for Argentata dell'Etna and Rossa Mediterranea. The resolution of the haplotype diversity at the casein cluster can be exploited both for selective and conservative plans.

Leptin Gene Polymorphism in Goats Fed with Diet at Different Energy Level: Effects on Feed Intake, Milk Traits, Milk Fatty Acids Composition, and Metabolic State.

The study investigated the effects of a polymorphism at the LEP gene intron 1 microsatellite region and its interaction with diet energy level on feed intake, milk traits, milk fatty acid composition, and metabolic state in goats. Sixteen Girgentana lactating goats at mid-lactation, selected on the basis of their genotype (8 goats homozygous 266 bp/266 bp, L genotype; 8 goats heterozygous 266 bp/264 bp, H genotype), were fed ad libitum according to a change-over design, with two diets at different energy levels reached with different hay inclusion: low energy diet (LE)-100% of hay; and high energy diet (HE)-65% of hay. No differences in milk yield and composition or in dry matter intake were found between leptin genotypes or between diets. Leptin genotype had no effect on plasma metabolite concentrations. The differences between diets were recorded for plasma ß-hydroxybutyric acid (BHBA) concentrations with higher (p = 0.01) values for the HE compared to the LE diet (0.44 vs. 0.24 mmol/L, respectively). Nonesterified fatty acid (NEFA) values seem to indicate a positive energy balance in goats. No interaction genotype per diet was evident for most of the studied parameters. Fatty acid composition was strongly influenced by LEP genotype: L goats, compared to H goats, showed higher levels of monounsaturated fatty acids (MUFA), polyunsaturated fatty acids (PUFA), and 14:1/14:0 desaturation index; lower levels of saturated fatty acids (SFA); and a more favorable atherogenic index. These results seem to suggest an improvement of health characteristics of milk with the L genotype.

Polymorphism at donkey ß-lactoglobulin II locus: identification and characterization of a new genetic variant with a very low expression.

In the last years, donkey milk had evidenced a renewed interest as a potential functional food and a breast milk substitute. In this light, the study of the protein composition assumes an important role. In particular, ß-lactoglobulin (ß-LG), which is considered as one of the main allergenic milk protein, in donkey species consists of two molecular forms, namely ß-LG I and ß-LG II. In the present research, a genetic analysis coupled with a proteomic approach showed the presence of a new allele, here named F, which is apparently associated with a null or a severely reduced expression of ß-LG II protein. The new ß-LG II F genetic variant shows a theoretical average mass (Mav) of 18,310.64 Da, a value practically corresponding with that of the variant D (∆mass < 0.07 Da), but differs from ß-LG II D for two amino acid substitutions: Thr100 (variant F) → Ala100 (variant D) and Thr118 (variant F) → Met118 (variant D). Proteomic investigation of the whey protein fraction of an individual milk sample, homozygous FF at ß-LG II locus, allowed to identify, as very minor component, the new ß-LG II F genetic variant. By MS/MS analysis of enzymatic digests, the sequence of the ß-LG II F was characterized, and the predicted genomic data confirmed.

Genetic diversity and variability in endangered Pantesco and two other Sicilian donkey breeds assessed by microsatellite markers.

The genetic variability of Pantesco and other two Sicilian autochthonous donkey breeds (Ragusano and Grigio Siciliano) was assessed using a set of 14 microsatellites. The main goals were to describe the current differentiation among the breeds and to provide genetic information useful to safeguard the Pantesco breed as well as to manage Ragusano and Grigio Siciliano. In the whole sample, that included 108 donkeys representative of the three populations, a total of 85 alleles were detected. The mean number of alleles was lower in Pantesco (3.7), than in Grigio Siciliano and Ragusano (4.4 and 5.9, resp.). The three breeds showed a quite low level of gene diversity (He) ranging from 0.471 in Pantesco to 0.589 in Grigio. The overall genetic differentiation index (Fst) was quite high; more than 10% of the diversity was found among breeds. Reynolds' (D(R)) genetic distances, correspondence, and population structure analysis reproduced the same picture, revealing that, (a) Pantesco breed is the most differentiated in the context of the Sicilian indigenous breeds, (b) within Ragusano breed, two well-defined subgroups were observed. This information is worth of further investigation in order to provide suitable data for conservation strategies.

Genetic analysis of Sicilian autochthonous horse breeds using nuclear and mitochondrial DNA markers.

Genetic diversity and relationship among 3 Sicilian horse breeds were investigated using 16 microsatellite markers and a 397-bp length mitochondrial D-loop sequence. The analysis of autosomal DNA was performed on 191 horses (80 Siciliano [SIC], 61 Sanfratellano [SAN], and 50 Sicilian Oriental Purebred [SOP]). SIC and SAN breeds were notably higher in genetic variability than the SOP. Genetic distances and cluster analysis showed a close relationship between SIC and SAN breeds, as expected according to the breeds' history. Sequencing of hypervariable mitochondrial DNA region was performed on a subset of 60 mares (20 for each breed). Overall, 20 haplotypes with 31 polymorphic sites were identified: A higher haplotype diversity was detected in SIC and SAN breeds, with 13 and 11 haplotypes respectively, whereas only one haplotype was found in SOP. These were compared with 118 sequences from GenBank. BLAST showed that 17 of the 20 haplotypes had been reported previously in other breeds. One haplotype, found in SIC, traces back to a Bronze Age archaeological site (Inner Mongolia). The 3 Sicilian breeds are now rare, and 2 of them are officially endangered. Our results represent a valuable tool for management strategies as well as for conservation purposes.