Genome-wide analysis reveals the patterns of genetic diversity and population structure of 8 Italian local chicken breeds.

The aim of this study was to conduct a genome-wide comparative analysis of 8 local Italian chicken breeds (Ermellinata di Rovigo, Millefiori di Lonigo [PML], Polverara Bianca, Polverara Nera, Padovana, Pepoi [PPP], Robusta Lionata, and Robusta Maculata), all under a conservation plan, to understand their genetic diversity and population structure. A total of 152 animals were analyzed using the Affymetrix Axiom 600 K Chicken Genotyping Array. The levels of genetic diversity were highest and lowest in PML and PPP, respectively. The results of genomic inbreeding based on runs of homozygosity (ROH; FROH) showed marked differences among breeds and ranged from 0.161 (PML) to 0.478 (PPP). Furthermore, in all breeds, short ROH (<4 Mb in length) were more frequent than long segments. Patterns of genetic differentiation, model-based clustering, and neighbor networks showed that most breeds formed nonoverlapping clusters and were clearly separate populations. The 2 Polverara breeds shared a similar genetic background and showed the lowest genetic differentiation in comparison with purebred lines; the local populations showed separated groups. PPP and PML were closer to the group of the purebred broiler lines (BRSA, BRSB, BRDA, and BRDB). Six genomic regions are presented as hotspots of autozygosity among the Italian chicken breeds, with candidate genes involved in multiple morphological phenotypes as breast muscle, muscle dry matter content, and body weight. This study is the first exhaustive genome-wide analysis of the diversity of these Italian local chickens from Veneto region. We conclude that breeds have conserved authentic genetic patterns. The results are of significant importance because they will help design and implement conservation strategies. In fact, the conservation of these breeds may also have positive impacts on the local economy, niche traditional markets, and offering a source of high-quality products to consumers. In this context, genomic information may play a crucial role in the management of local breeds.

Genome-wide association studies for milk production traits in Valle del Belice sheep using repeated measures.

Genome-wide association studies (GWASes) have become a powerful tool for identifying genomic regions associated with important traits in livestock. Milk production traits in dairy sheep are measured at different time points during their life span. Using phenotypic data generated from longitudinal traits could improve the power of association studies but until now have received less attention in GWASes as a methodology and has not been implemented. The aim of this study was to carry out a GWAS for milk production traits in Valle del Belice sheep using repeated measures. After quality control, 469 ewes and 37 228 SNPs were retained for the analysis, and phenotypic data included 5586 test-day records for five milk production traits (milk yield, MY; fat yield and percentage, FY and F%; protein yield and percentage, PY and P%). Nine SNPs located within or close to known genes were found to be associated with milk production traits. In particular, rs398340969, associated with both milk yield and protein yield, is located within the DCPS gene. In addition, rs425417915 and rs417079368, both associated with both fat percentage and protein percentage, are located within the TTC7B gene and at 0.37 Mb within the SUCNR1 gene respectively. In summary, the use of repeated records was beneficial for mapping genomic regions affecting milk production traits in the Valle del Belice sheep.

Genome-wide association study reveals the locus responsible for microtia in Valle del Belice sheep breed.

Microtia is a congenital deformity of the outer ear with phenotypes varying from a small auricle to total absence (anotia). The genetic basis is still poorly understood, and very few studies have been performed in sheep. Valle del Belice sheep is a breed showing microtia. The aim of this study was to identify the potential genomic regions involved in microtia in sheep. A total of 40 individuals, 20 with microtia and 20 normal, were genotyped with the Illumina OvineSNP50 BeadChip. The comparison among the results from a genome-wide association study, Fisher's exact test and FST analysis revealed a single strong association signal: rs419889303 on chromosome 1, located within intron 3 of the CLRN1 gene. Our study suggests for the first time that this novel candidate gene is responsible for microtia in sheep. Additional analysis based on the sequencing would help confirm our findings and allow for the proposal of a precise genetic basis for microtia in sheep.

Genome-wide identification of runs of homozygosity islands and associated genes in local dairy cattle breeds.

Runs of homozygosity (ROH) are widely used as predictors of whole-genome inbreeding levels in cattle. They identify regions that have an unfavorable effect on a phenotype when homozygous, but also identify the genes associated with traits of economic interest present in these regions. Here, the distribution of ROH islands and enriched genes within these regions in four dairy cattle breeds were investigated. Cinisara (71), Modicana (72), Reggiana (168) and Italian Holstein (96) individuals were genotyped using the 50K v2 Illumina BeadChip. The genomic regions most commonly associated with ROHs were identified by selecting the top 1% of the single nucleotide polymorphisms (SNPs) most commonly observed in the ROH of each breed. In total, 11 genomic regions were identified in Cinisara and Italian Holstein, and eight in Modicana and Reggiana, indicating an increased ROH frequency level. Generally, ROH islands differed between breeds. The most homozygous region (>45% of individuals with ROH) was found in Modicana on chromosome 6 within a quantitative trail locus affecting milk fat and protein concentrations. We identified between 126 and 347 genes within ROH islands, which are involved in multiple signaling and signal transduction pathways in a wide variety of biological processes. The gene ontology enrichment provided information on possible molecular functions, biological processes and cellular components under selection related to milk production, reproduction, immune response and resistance/susceptibility to infection and diseases. Thus, scanning the genome for ROH could be an alternative strategy to detect genomic regions and genes related to important economic traits.

Penalized classification for optimal statistical selection of markers from high-throughput genotyping: application in sheep breeds.

The identification of individuals' breed of origin has several practical applications in livestock and is useful in different biological contexts such as conservation genetics, breeding and authentication of animal products. In this paper, penalized multinomial regression was applied to identify the minimum number of single nucleotide polymorphisms (SNPs) from high-throughput genotyping data for individual assignment to dairy sheep breeds reared in Sicily. The combined use of penalized multinomial regression and stability selection reduced the number of SNPs required to 48. A final validation step on an independent population was carried out obtaining 100% correctly classified individuals. The results using independent analysis, such as admixture, F st, principal component analysis and random forest, confirmed the ability of these methods in selecting distinctive markers. The identified SNPs may constitute a starting point for the development of a SNP based identification test as a tool for breed assignment and traceability of animal products.

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 breed's 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.

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.