The SYNBREED chicken diversity panel: a global resource to assess chicken diversity at high genomic resolution

Selection of Male Donors in Local Chicken Breeds to Implement the Italian Semen Cryobank: Variability in Semen Quality, Freezability and Fertility

Native breed conservation is an important component of poultry biodiversity. The aim of this work is to describe different steps that lead to donor selection for the implementation of the Italian Semen Cryobank of Autochthonous Chicken and Turkey Breeds. The variability within and between breeds was evaluated, and the stored semen reproductive capacity was in vivo tested using artificial insemination. Semen from Bionda Piemontese, Bianca di Saluzzo and Pepoi roosters was collected and processed. Concentration, volume, sperm membrane integrity, total motile sperm, progressive motile sperm and kinetic parameters were analyzed; sperm parameters accounting for bird variability were used to select male donors. Fresh semen quality parameters measured in donor ejaculates showed significant differences between breeds; no differences were found after cryopreservation. Variability in the fertilizing ability of cryopreserved semen was found within a breed (5-16%) and between birds within a breed (BP = 3-7%; BS = 7-31%; PP = 6-22%); only sperm quality parameters measured in fresh ejaculates, not frozen/thawed, may be associated with in vivo fertility results. In conclusion, sperm concentration and progressive motility were successfully used as selection parameters to identify chicken male donors with improved sperm quality for sperm cryobanking. However, new reliable sperm markers to predict cryopreserved semen's fertilizing ability are required.

Raising laying hens: housing complexity and genetic strain affect startle reflex amplitude and behavioural response to fear-inducing stimuli

Individual variation in fearfulness can be modified during ontogeny, and high levels of fear can affect animal welfare. We asked whether early-life environmental complexity and genetic strain affect fear behaviour in young laying hens (pullets). Four replicates of brown (B) and white (W) genetic strains (breeds) of layers were each raised in four environmental treatments (housing): conventional cages (Conv) and different rearing aviaries with increasing space and complexity (Low < Mid < High). We used a startle reflex test (weeks 4 and 14) to measure startle amplitude and autonomic response (i.e. comb temperature). A combination of novel arena (NA) and novel object (NO) tests was used (week 14) to assess NA exploration and alertness, latency to approach the centre and initial NO avoidance and investigation. Housing × strain affected startle amplitude (B-Conv, B-High < B-Low, B-Mid; B > W; no housing effect in W) but not autonomic response. Fear behaviour was affected by housing (NA exploration, investigation: Conv < Low, Mid, High; NO avoidance: Conv, High < Low, Mid), strain (NA alertness: B > W, NO avoidance: W > B) and their interaction (NA centre approach: B-Conv < all other groups). We present evidence for strain-specific fear responses depending on early experience.

Comparative peculiarities of genomic diversity in Gallus gallus domesticus chickens with decorative plumage: the muffs and beard phenotype

Throughout history, humans have been attempting to develop the ornamental features of domestic animals in addition to their productive qualities. Many chicken breeds have developed tufts of elongated feathers that jut out from the sides and bottom of the beak, leading to the phenotype known as muffs and beard. It is an incomplete autosomal dominant phenotype determined by the Mb locus localised on chromosome GGA27. This project aimed to analyse the genetic diversity of chicken breeds using full genomic genotyping with the Chicken 60K BeadChip. A total of 53,313 Single Nucleotide Polymorphisms were analysed. DNA was obtained from breeds with the muffs and beard as a marker phenotype: Faverolles (n = 20), Ukrainian Muffed (n = 18), Orloff (n = 20), Novopavlov White (n = 20), and Novopavlov Coloured (n = 15). The Russian White (n = 20) was selected as an alternative breed without the muffs and beard phenotype. The chickens are owned by the Centre of Collective Use "Genetic Collection of Rare and Endangered Breeds of Chickens" (St. Petersburg region, Pushkin), and are also included in the Core Shared Research Facility (CSRF) and/or Large-Scale Research Facility (LSRF). Multidimensional scaling revealed that the Novopavlov White and the Novopavlov Coloured populations formed a separate group. The Ukrainian Muffed and the Orloff have also been combined into a separate group. Based on cluster analysis, with the cross-validation error and the most probable number of clusters K = 4 taken into account, the Orloff was singled out as a separate group. The Ukrainian Muffed exhibited a notable similarity with the Orloff under the same conditions. At K = 5, the populations of the Novopavlov White and the Novopavlov Coloured diverged. Only at K = 6, a distinct and separate cluster was formed by the Ukrainian Muffed. The Russian White had the greatest number of short (1-2 Mb) homozygous regions. If the HOXB8 gene is located between 3.402 and 3.404 Mb on chromosome GGA27, homozygous regions are rarely found in the chickens with the muffs and beard phenotype. Scanning the chicken genome with the Chicken 60K BeadChip provided enough information about the genetic diversity of the chicken breeds for the peculiarities of the development of the ornamental muffs and beard phenotypes in them to be understood. For example, Phoenix bantams, whose tail feathers grow throughout their lives, require greater consideration of husbandry conditions.

The novel IMAGE001 genotyping array as a valuable alternative for genetic diversity screening in chicken: a demonstration in a local chicken breed in Belgium

Screening for genetic diversity in livestock species breeds is of utmost importance, especially for local, small populations that are at the risk of extinction. Luckily, recent developments in technology increase access to genotyping, also for numerically small breeds. One of these new technologies is the IMAGE001 single nucleotide polymorphism genotyping array that includes markers for 6 different species (cow, pig, sheep, chicken, horse and goat). For our current study, we studied the Turkey-headed Malines chicken, a local chicken breed in Belgium, for the first time. A total of 110 animals were genotyped, together with 29 samples from 4 supposedly related breeds. The genotypes were used to assess the genetic diversity of this local breed. Our analysis revealed an average inbreeding coefficient of 0.20 through runs of homozygosity analysis, and effective population size estimation based on linkage disequilibrium indicated a low genetic diversity (Ne = 34). Moreover, a principal component analysis and a genetic differentiation study (FST) were performed using these marker data to position the Turkey-headed Malines relative to the 4 other indigenous Belgian chicken breeds. Finally, we discussed the practical implications of the overlap between the IMAGE001 array and other existing chicken genotyping arrays. This study is the first use of the novel IMAGE001 array to evaluate a local chicken breed, and demonstrates it as a viable option for genomic characterization a breed. Moreover, with this research, we are able to provide a good basis for further evaluation of the Belgian chicken heritage.

A pangenome graph reference of 30 chicken genomes allows genotyping of large and complex structural variants

BACKGROUND

The red junglefowl, the wild outgroup of domestic chickens, has historically served as a reference for genomic studies of domestic chickens. These studies have provided insight into the etiology of traits of commercial importance. However, the use of a single reference genome does not capture diversity present among modern breeds, many of which have accumulated molecular changes due to drift and selection. While reference-based resequencing is well-suited to cataloging simple variants such as single-nucleotide changes and short insertions and deletions, it is mostly inadequate to discover more complex structural variation in the genome.

METHODS

We present a pangenome for the domestic chicken consisting of thirty assemblies of chickens from different breeds and research lines.

RESULTS

We demonstrate how this pangenome can be used to catalog structural variants present in modern breeds and untangle complex nested variation. We show that alignment of short reads from 100 diverse wild and domestic chickens to this pangenome reduces reference bias by 38%, which affects downstream genotyping results. This approach also allows for the accurate genotyping of a large and complex pair of structural variants at the K feathering locus using short reads, which would not be possible using a linear reference.

CONCLUSIONS

We expect that this new paradigm of genomic reference will allow better pinpointing of exact mutations responsible for specific phenotypes, which will in turn be necessary for breeding chickens that meet new sustainability criteria and are resilient to quickly evolving pathogen threats.

A Study on Egg Production and Quality According to the Age of Four Italian Chicken Dual-Purpose Purebred Hens Reared Outdoors

The month of hatching and the rearing management, especially temperature and photoperiod, are important factors for pullets and hens reared outdoors. The yield performance and egg quality of dual-purpose chicken breeds from the Veneto region (Italy), Pepoi (PP), Ermellinata di Rovigo (ER), Robusta Maculata (RM) and Robusta Lionata (RL), with different adult body weights (ABW, kg, PP = 1.3; ER = 2.3, RM and RL = 3.1), were studied, using a factorial model (4 × 2), considering breed and age (26-33 weeks, first age, summer-autumn, under decreasing natural photoperiod-on average, 12L:12D, and 42-53 weeks, second age, winter, under implemented photoperiod-14L:10D) as the main effects and interaction. The chicks hatched in spring, and they started laying at the end of summer/beginning of autumn. Significant (p < 0.05) results were shown for many traits. ER showed higher hen-day egg production than that of PP, and RM and RL were the lowest; ER, RM and RL showed medium-size eggs and PP showed small-size eggs. RM produced the most spherical eggs and ER the most ovoid, and they showed the highest and the lowest eggshell thickness, respectively. RM showed the highest yolk to albumen ratio, and RL showed the lowest. The age increased the laying rate and the egg weight in all the groups. At 26-33 weeks, ER and PP showed higher hen-day egg production (on average 24%) than RM and RL (on average, less than 10%). The onset of laying (at least 10% laying rate) was shown, at different ages, according to the % ABW the breeds had reached: PP was the first, followed by ER, then RM, and RL was the last. At 42-53 weeks, the hen-day egg production ranged, on average, from 38 to 52%, according to the breeds; orthogonal contrasts on two-weekly data showed, at first age, increasing linear (ER) and quadratic (other groups) trends, and at second age, positive linear (ER, RM) and cubic (PP, RL) trends. Age (32 vs. 53 weeks) affected almost all the eggshell traits in PP and ER, whereas in RL, and especially RM, fewer traits changed. The age increased the yolk to albumen ratio (unchanged in PP). These results may be useful for the effective management of local purebred chickens, with the purpose to ensure the wellbeing of the hens and for supplying eggs of different quality throughout the year.

Discovery of novel MHC-B haplotypes in Chantecler chickens

The major histocompatibility complex (MHC) is a highly polymorphic cluster of genes which contribute to immune response. Located on chromosome 16, the chicken MHC has great influence over disease resistance and susceptibility. Through the use of a high-density SNP panel which encompasses the MHC-B region, haplotypes can be easily identified. This study aims to use an MHC-B SNP panel to evaluate the MHC-B variability in the Chantecler breed. This breed is native to Quebec, Canada, and is a dual-purpose breed known for its strong resistance to extreme cold temperatures. The Chantecler breed faced a near extinction event in the 1970s, which most likely resulted in a genetic bottleneck and loss of diversity. Despite this, SNP haplotype diversity was observed among 4 Chantecler populations. A total of 8 haplotypes were observed. Of these haplotypes, 6 were previously defined in other breeds, and the other 2 were unique to the Chantecler. Within the populations, the number of haplotypes ranged from 4 to 7, with 3 haplotypes, including the novel BSNP-Chant01, being present in all the groups. This study shows existence of reasonable diversity in the MHC-B region of the Chantecler breed and our results further contribute to understanding the variability of this region in chickens.

A gene bank's collection of genetic diversity among minor chicken breeds

Genetic differences among heritage or fancier breeds of chickens have not been quantified in the United States. Gene banks collecting germplasm for conserving these breeds need this information as do breeders and companies raising them. Our goal was to evaluate genetic diversity of 10 heritage/fancier chicken breeds that are a component of the national collection and to use this information to establish a baseline of their genetic diversity and future conservation efforts. Breeds could be broadly classified as European, Asian, Mediterranean, and United States (US) in origin. The US breeds were composite breeds developed between the 1849 and 1935. Animals (n = 24-31 per breed) were sampled for DNA analysis from 2 or 3 hatcheries per breed and a total of 8 hatcheries. The hatcheries were assumed to maintain and breed their own populations of the studied breeds. Effective population sizes ranged from 47 to 145 and used to estimate probabilities of extinction for a 50-generation timeline. It was determined that Crevecoeur and Aseel had a probability of extinction that exceeded 40%, the remaining 8 breeds had probabilities of <28%. ADMIXTURE analysis indicated the minimal CV corresponded to 9 populations. In that analysis New Hampshire and Rhode Island Red were classified as the same population, which was not unusual given that New Hampshire was developed as a subpopulation of Rhode Island Red. Crevecoeur and Buttercup were the 2 most genetically divergent breeds based on pairwise Fst among the breeds and principal component analysis, which was supported by the ADMIXTURE results. Inbreeding coefficients computed from genomic information was lowest for Crevecoeur, Rhode Island Red, Buttercup, and Andalusian (0.8-2.6%), while New Hampshire, Buckeye, and Aseel were highest (12.8-14.3%). Within breed Fst among hatcheries supplying animals for sampling generally indicated a genetic structure was present on a breed-by-breed basis. Genetic relationships within hatchery were also computed for each breed. Several of the hatcheries had sent samples that suggested genetic relationships as high as half-sibs while several others had genetic relationships closer to first cousins. We conclude that the chicken breeds evaluated have substantial genetic variability within the in situ populations and the gene bank has captured this diversity for future use.

The development of laying hen locomotion in 3D space is affected by early environmental complexity and genetic strain

Adult laying hens are increasingly housed in spatially complex systems, e.g., non-cage aviaries, where locomotion between elevated structures can be challenging for these gallinaceous birds. This study assessed the effect of early environmental complexity on spatial skills in two genetic strains. Brown (B) and white (W) feathered birds were raised in: Conventional cages with minimal complexity (Conv) or rearing aviaries with low (Low), intermediate (Mid), or high complexity (High). Birds from each housing treatment were challenged at three different time points in three different, age-appropriate vertical spatial tasks. Whites performed better than brown birds in all tests regardless of rearing environment. In chicks, test performance was predominantly explained by variation between replicates and differences in motivation for test participation. Treatment effects were seen in pubertal birds (pullets), with pullets from aviaries performing better than those from Conv. White High pullets performed better than white Mid or Low, an effect that was not found in browns. Pullets preferred to use a ramp to move downwards, but only when ramps had previously been experienced and when the ramp was not too steep. Overall, early environmental complexity affected spatial skills of laying hen pullets with stronger effects in white than brown feathered birds.

Assessment of long-term trends in genetic mean and variance after the introduction of genomic selection in layers: a simulation study

Nucleus-based breeding programs are characterized by intense selection that results in high genetic gain, which inevitably means reduction of genetic variation in the breeding population. Therefore, genetic variation in such breeding systems is typically managed systematically, for example, by avoiding mating the closest relatives to limit progeny inbreeding. However, intense selection requires maximum effort to make such breeding programs sustainable in the long-term. The objective of this study was to use simulation to evaluate the long-term impact of genomic selection on genetic mean and variance in an intense layer chicken breeding program. We developed a large-scale stochastic simulation of an intense layer chicken breeding program to compare conventional truncation selection to genomic truncation selection optimized with either minimization of progeny inbreeding or full-scale optimal contribution selection. We compared the programs in terms of genetic mean, genic variance, conversion efficiency, rate of inbreeding, effective population size, and accuracy of selection. Our results confirmed that genomic truncation selection has immediate benefits compared to conventional truncation selection in all specified metrics. A simple minimization of progeny inbreeding after genomic truncation selection did not provide any significant improvements. Optimal contribution selection was successful in having better conversion efficiency and effective population size compared to genomic truncation selection, but it must be fine-tuned for balance between loss of genetic variance and genetic gain. In our simulation, we measured this balance using trigonometric penalty degrees between truncation selection and a balanced solution and concluded that the best results were between 45° and 65°. This balance is specific to the breeding program and depends on how much immediate genetic gain a breeding program may risk vs. save for the future. Furthermore, our results show that the persistence of accuracy is better with optimal contribution selection compared to truncation selection. In general, our results show that optimal contribution selection can ensure long-term success in intensive breeding programs using genomic selection.

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.

Differences among domestic chicken breeds in tonic immobility responses as a measure of fearfulness

Background

One priority for animal welfare is for animals to experience less fear, especially during human contact. For domestic animals, breeds that are less fearful may provide genetic resources to develop strains with improved welfare due to lower susceptibility to fear. Genetic predispositions inherited in these breeds might reflect the large diversity of chicken breeds. The goal of the present study was to systematically test a diverse group of chicken breeds to search for breeds that experience less fear.

Methods

Nineteen chicken breeds from commercial hybrid lines, native layer-type, meat-type and dual-purpose breeds, ornamental breeds as well as bantam breeds were tested in a standardized tonic immobility (TI) test. Chickens were manually restrained on their back, and the time to first head movement and first leg movement, the duration of TI, as well as the number of attempts needed to induce TI were measured.

Results

The TI response differed among chicken breeds (p ≤ 0.001) for naïve, mature hens. The median number of attempts required to induce TI ranged from 1 to 2 and did not differ significantly among breeds. Median durations were much more variable, with Lohmann Brown showing shortest durations (6 s, 12 s, 58 s for time to first head movement, first leg movement and total duration of TI, respectively). In contrast, medians reached the maximum of 600 s for all three measures in German Creepers. Repeated tests on the same individuals did not affect attempts needed to induce TI nor TI durations. Breeds clustered into two main groups, with layer-type native breeds and ornamental breeds having longer TI durations, and bantam, dual-purpose and meat-type native breeds having shorter TI durations.

Conclusions

Our findings provide evidence for substantial variation of fearfulness among breeds. This variation could be linked to the intended use during the breed's specific history. Knowledge and quantitative measurement of these behavioural responses provide the opportunity to improve welfare through selection and future breeding.

Fearfulness in Commercial Laying Hens: A Meta-Analysis Comparing Brown and White Egg Layers

High fearfulness in commercial laying hens can negatively affect production parameters and animal welfare. Brown and white egg layers differ in several behavioral characteristics, though reported differences in fearfulness are inconsistent. A meta-analysis was conducted to determine whether there are systematic differences in measures of fearfulness between brown and white layers. Twenty-three studies that examined either 1 or both of 2 behavioral tests were included: tonic immobility (TI) (longer duration = higher fearfulness, 16 studies) and novel object (NO) test (lower approach rate = higher fearfulness, 11 studies). The 2 tests were analyzed separately. TI analyses: A generalized linear mixed effect model (GLMM) with a lognormal distribution was fitted to describe the data with experiment nested in study as a random effect. Explanatory (X) variables were considered through backward selection, where potential X-variables included color (brown vs. white layers), decade (1980s, 2000s, 2020s), age (prelay vs. in lay), genetic stock (hybrid vs. grand-/parent stock), and methodology (back vs. side position). NO test analyses: univariable GLMMs with a beta distribution were fitted with approach rate as the Y-variable and color, decade, age, stock, or 2 methodological factors (test duration, single vs. group testing) as X-variables. Models were evaluated by assessing information criteria, residuals/random effect normality, significance of X-variables and model evaluation statistics (mean square prediction error, concordance correlation coefficient). TI duration was best explained by a color-by-decade interaction (P = 0.0006). Whites in the 1980s had longer TI durations (709.43 ± 143.88 s) than browns in the 1980s (282.90 ± 59.70 s), as well as in comparison to browns (208.80 ± 50.82 s) or whites (204.85 ± 49.60 s) in the 2020s. The NO approach rate was best explained by color (P ≤ 0.05 in 3 models), age (P < 0.05 in 3 models), and decade (P = 0.04). Whites had a higher approach rate (0.7 ± 0.07) than browns (0.5 ± 0.11), birds in lay a higher rate (0.8 ± 0.07) than birds prelay (0.4 ± 0.12), and approach rate for papers published in the 2000s (0.8 ± 0.09) was higher than in the 2020s (0.2 ± 0.12). The phylogenetic difference in the 1980s was no longer detectable after enforcing an upper limit on TI durations (10 min), as became common practice in later studies. Our findings suggest that phylogenetic differences in fearfulness and changes over time are test dependent, and this raises important questions and potential consequences for assessing hen welfare in commercial egg production.

Analysis of Conservation Priorities and Runs of Homozygosity Patterns for Chinese Indigenous Chicken Breeds

To achieve sustainable development of the poultry industry, the effective conservation of genetic resources has become increasingly important. In the present study, we systematically elucidated the population structure, conservation priority, and runs of homozygosity (ROH) patterns of Chinese native chicken breeds. We used a high-density genotyping dataset of 157 native chickens from eight breeds. The population structure showed different degrees of population stratification among the breeds. Chahua chicken was the most differentiated breed from the other breeds (Nei = 0.0813), and the Wannan three-yellow chicken (WanTy) showed the lowest degree of differentiation (Nei = 0.0438). On the basis of contribution priority, Xiaoshan chicken had the highest contribution to the total gene diversity (1.41%) and the maximum gene diversity of the synthetic population (31.1%). WanTy chicken showed the highest contribution to the total allelic diversity (1.31%) and the maximum allelic diversity of the syntenic population (17.0%). A total of 5242 ROH fragments and 5 ROH island regions were detected. The longest ROH fragment was 41.51 Mb. A comparison of the overlapping genomic regions between the ROH islands and QTLs in the quantitative trait loci (QTL) database showed that the annotated candidate genes were involved in crucial economic traits such as immunity, carcass weight, drumstick and leg muscle development, egg quality and egg production, abdominal fat precipitation, body weight, and feed intake. In conclusion, our findings revealed that Chahua, Xiaoshan, and WanTy should be the priority conservation breeds, which will help optimize the conservation and breeding programs for Chinese indigenous chicken breeds.

Temporal transcriptome profiling of floating apical out chicken enteroids suggest stability and reproducibility

Enteroids are miniature self-organising three-dimensional (3D) tissue cultures which replicate much of the complexity of the intestinal epithelium. We recently developed an apical-out leukocyte-containing chicken enteroid model providing a novel physiologically relevant in vitro tool to explore host-pathogen interactions in the avian gut. However, the replicate consistency and culture stability have not yet been fully explored at the transcript level. In addition, causes for the inability to passage apical-out enteroids were not determined. Here we report the transcriptional profiling of chicken embryonic intestinal villi and chicken enteroid cultures using bulk RNA-seq. Comparison of the transcriptomes of biological and technical replicate enteroid cultures confirmed their high level of reproducibility. Detailed analysis of cell subpopulation and function markers revealed that the mature enteroids differentiate from late embryonic intestinal villi to recapitulate many digestive, immune and gut-barrier functions present in the avian intestine. These transcriptomic results demonstrate that the chicken enteroid cultures are highly reproducible, and within the first week of culture they morphologically mature to appear similar to the in vivo intestine, therefore representing a physiologically-relevant in vitro model of the chicken intestine.

Temporal expression of avian β defensin 10 and cathelicidins in the yolk sac tissue of broiler and layer embryos

The yolk sac is a multifunctional organ, which not only participates in nutrient absorption, but also plays an important role in immune function. The objective of this study was to compare the mRNA abundance of avian β-defensin 10 (AvBD10) and 3 cathelicidins (CATH1, CATH2, and CATH3) in the yolk sac tissue (YST) of commercial broilers and white egg and brown egg commercial layers. AvBD10 and CATH mRNA abundance was analyzed using two-way ANOVA and Tukey's test, with P < 0.05 being considered significant. AvBD10 and CATH mRNA showed similar temporal expression patterns in the YST of both broiler and layers, with an increase from embryonic day (E) 7 to E9 through E13 followed by a decrease to day of hatch. AvBD10 mRNA showed a breed × age interaction with greater expression in the YST of both layers compared to broilers at E9 and E11. CATH1 mRNA was greater in the YST of brown egg layers than broilers. CATH2 mRNA showed a breed × age interaction, with greater expression in the YST of brown egg layers than broilers at E11. CATH3 mRNA showed no difference in the YST between layers and broilers. Because broilers and brown egg layers are genetically related, these results show that selection for production parameters (broiler vs. layer) and not genetic relatedness (white egg layer vs. brown egg layer and broilers) is the basis for the differences in AvBD10, CATH1, and CATH2 mRNA in the YST of broilers and layers. The yolk-free body weights of broiler embryos were greater than that of both brown and white egg layers from E9 to 17. One possible explanation is that the reduced expression of AvBD10, CATH1 and CATH2 mRNA in the YST of broilers compared to layers at E9 and 11 may be due to faster embryonic growth at the expense of host defense peptide expression in broilers compared to layers.

High-density genotyping reveals candidate genomic regions for chicken body size in breeds of Asian origin

Body size is one of the main selection indices in chicken breeding. Although often investigated, knowledge of the underlying genetic mechanisms is incomplete. The aim of the current study was to identify genomic regions associated with body size differences between Asian Game and Asian Bantam type chickens. In this study, 94 and 107 chickens from 4 Asian Game and 5 Asian Bantam type breeds, respectively, were genotyped using the chicken 580K single nucleotide polymorphism (SNP) array. A genome-wide association study (GWAS) and principal component analyses (PCA) were performed to identify genomic regions associated with body size related-traits such as wing length, shank length, shank thickness, keel length, and body weight. Hierarchical clustering of genotype data showed a clear genetic difference between the investigated Asian Game and Asian Bantam chicken types. GWAS identified 16 genomic regions associated with wing length (2, FDR ≤ 0.018), shank thickness (6, FDR ≤ 0.008), keel length (5, FDR ≤ 0.023), and body weight (3, FDR ≤ 0.041). PCA showed that the first principal component (PC1) separated the 2 chicken types and significantly correlated with the measured body size related-traits (P ≤ 2.24e-40). SNPs contributing significantly to PC1 were subjected to a more detailed investigation. This analysis identified 11 regions potentially associated with differences in body size related-traits. A region on chromosome 4 (GGA4) (17.3-21.3 Mb) was detected in both analyses GWAS and PCA. This region harbors 60 genes. Among them are myotubularin 1 (MTM1) and secreted frizzled-related protein 2 (SFPR2) which can be considered as potential candidate genes for body size related-traits. Our results clearly show that the investigated Asian Game type chicken breeds are genetically different from the Asian Bantam breeds. A region on GGA4 between 17.3 and 21.3 Mb was identified which contributes to the phenotypic difference, though further validation of candidate genes is necessary.

Comparative Analysis of Molecular RFLP and SNP Markers in Assessing and Understanding the Genetic Diversity of Various Chicken Breeds

Monitoring the genetic diversity of small populations is important with respect to conserving rare and valuable chicken breeds, as well as discovery and innovation in germplasm research and application. Restriction fragment length polymorphisms (RFLPs), the molecular markers that underlie multilocus DNA fingerprinting (MLDF), have historically been employed for this purpose, but over the past two decades, there has been an irreversible shift toward high-throughput single-nucleotide polymorphisms (SNPs). In this study, we conducted a comparative analysis of archived MLDF results and new data from whole-genome SNP genotyping (SNPg) among 18 divergently selected breeds representing a large sample of the world gene pool. As a result, we obtained data that fit the general concept of the phylogenetic distribution of the studied breeds and compared them with RFLP and SNP markers. RFLPs were found to be useful markers for retrospective assessment of changes in the genetic architecture and variability underlying the phenotypic variation in chicken populations, especially when samples from previous generations used for MLDF are unavailable for SNPg. These results can facilitate further research necessary to assess the possibility of extrapolating previous MLDF results to study the long-term dynamics of genetic diversity in various small chicken germplasm populations over time. In general, the whole-genome characterization of populations and breeds by multiple SNP loci will further form the basis for the development and implementation of genomic selection with the aim of effective use of the genetic potential of the domestic gene pool in the poultry industry.

Single nucleotide polymorphism marker combinations for classifying Yeonsan Ogye chicken using a machine learning approach

Genetic analysis has great potential as a tool to differentiate between different species and breeds of livestock. In this study, the optimal combinations of single nucleotide polymorphism (SNP) markers for discriminating the Yeonsan Ogye chicken (Gallus gallus domesticus) breed were identified using high-density 600K SNP array data. In 3,904 individuals from 198 chicken breeds, SNP markers specific to the target population were discovered through a case-control genome-wide association study (GWAS) and filtered out based on the linkage disequilibrium blocks. Significant SNP markers were selected by feature selection applying two machine learning algorithms: Random Forest (RF) and AdaBoost (AB). Using a machine learning approach, the 38 (RF) and 43 (AB) optimal SNP marker combinations for the Yeonsan Ogye chicken population demonstrated 100% accuracy. Hence, the GWAS and machine learning models used in this study can be efficiently utilized to identify the optimal combination of markers for discriminating target populations using multiple SNP markers.

Unveiling new perspective of phylogeography, genetic diversity, and population dynamics of Southeast Asian and Pacific chickens

The complex geographic and temporal origins of chicken domestication have attracted wide interest in molecular phylogeny and phylogeographic studies as they continue to be debated up to this day. In particular, the population dynamics and lineage-specific divergence time estimates of chickens in Southeast Asia (SEA) and the Pacific region are not well studied. Here, we analyzed 519 complete mitochondrial DNA control region sequences and identified 133 haplotypes with 70 variable sites. We documented 82.7% geographically unique haplotypes distributed across major haplogroups except for haplogroup C, suggesting high polymorphism among studied individuals. Mainland SEA (MSEA) chickens have higher overall genetic diversity than island SEA (ISEA) chickens. Phylogenetic trees and median-joining network revealed evidence of a new divergent matrilineage (i.e., haplogroup V) as a sister-clade of haplogroup C. The maximum clade credibility tree estimated the earlier coalescence age of ancestral D-lineage (i.e., sub-haplogroup D2) of continental chickens (3.7 kya; 95% HPD 1985–4835 years) while island populations diverged later at 2.1 kya (95% HPD 1467–2815 years). This evidence of earlier coalescence age of haplogroup D ancestral matriline exemplified dispersal patterns to the ISEA, and thereafter the island clade diversified as a distinct group.

Identifying cancer tissue-of-origin by a novel machine learning method based on expression quantitative trait loci

Cancer of unknown primary (CUP) refers to cancer with primary lesion unidentifiable by regular pathological and clinical diagnostic methods. This kind of cancer is extremely difficult to treat, and patients with CUP usually have a very short survival time. Recent studies have suggested that cancer treatment targeting primary lesion will significantly improve the survival of CUP patients. Thus, it is critical to develop accurate yet fast methods to infer the tissue-of-origin (TOO) of CUP. In the past years, there are a few computational methods to infer TOO based on single omics data like gene expression, methylation, somatic mutation, and so on. However, the metastasis of tumor involves the interaction of multiple levels of biological molecules. In this study, we developed a novel computational method to predict TOO of CUP patients by explicitly integrating expression quantitative trait loci (eQTL) into an XGBoost classification model. We trained our model with The Cancer Genome Atlas (TCGA) data involving over 7,000 samples across 20 types of solid tumors. In the 10-fold cross-validation, the prediction accuracy of the model with eQTL was over 0.96, better than that without eQTL. In addition, we also tested our model in an independent data downloaded from Gene Expression Omnibus (GEO) consisting of 87 samples across 4 cancer types. The model also achieved an f1-score of 0.7-1 depending on different cancer types. In summary, eQTL was an important information in inferring cancer TOO and the model might be applied in clinical routine test for CUP patients in the future.

Managing genetic diversity in breeding programs of small populations: the case of French local chicken breeds

Background

On-going climate change will drastically modify agriculture in the future, with a need for more sustainable systems, in particular regarding animal production. In this context, genetic diversity is a key factor for adaptation to new conditions: local breeds likely harbor unique adaptive features and represent a key component of diversity to reach resilience. However, local breeds often suffer from small population sizes, which puts these valuable resources at risk of extinction. In chickens, population management programs were initiated a few decades ago in France, relying on a particular niche market that aims at promoting and protecting local breeds. We conducted a unique comprehensive study of 22 French local breeds, along with four commercial lines, to evaluate their genetic conservation status and the efficiency of the population management programs.

Results

Using a 57K single nucleotide polymorphism (SNP) chip, we demonstrated that both the between- and within-breed genetic diversity levels are high in the French local chicken populations. Diversity is mainly structured according to the breeds’ selection and history. Nevertheless, we observed a prominent sub-structuring of breeds according to farmers’ practices in terms of exchange, leading to more or less isolated flocks. By analysing demographic parameters and molecular information, we showed that consistent management programs are efficient in conserving genetic diversity, since breeds that integrated such programs earlier had older inbreeding.

Conclusions

Management programs of French local chicken breeds have maintained their genetic diversity at a good level. We recommend that future programs sample as many individuals as possible, with emphasis on both males and females from the start, and focus on a quick and strong increase of population size while conserving as many families as possible. We also stress the usefulness of molecular tools to monitor small populations for which pedigrees are not always available. Finally, the breed appears to be an appropriate operational unit for the conservation of genetic diversity, even for local breeds, for which varieties, if present, could also be taken into account.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-022-00746-2.

Whole Genome Analyses Reveal Novel Genes Associated with Chicken Adaptation to Tropical and Frigid Environments

INTRODUCTION

Investigating the genetic footprints of historical temperature selection can get insights to the local adaptation and feasible influences of climate change on long-term population dynamics.

OBJECT

Chicken is a significative species to study genetic adaptation on account of its similar domestication track related to human activity with the most diversified varieties. Yet, few studies have demonstrated the genetic signatures of its adaptation to naturally tropical and frigid environments.

METHOD

Here, we generated whole genome resequencing of 119 domesticated chickens in China including the following breeds which are in order of breeding environmental temperature from more tropical to more frigid: Wenchang chicken (WCC), green-shell chicken (GSC), Tibetan chicken (TBC), and Lindian chicken (LDC).

RESULTS

Our results showed WCC branched off earlier than LDC with an evident genetic admixture between WCC and LDC, suggesting their closer genetic relationship. Further comparative genomic analyses solute carrier family 33 member 1 (SLC33A1) and thyroid stimulating hormone receptor (TSHR) genes exhibited stronger signatures for positive selection in the genome of the more tropical WCC. Furthermore, genotype data from about 3,000 African local ecotypes confirmed that allele frequencies of single nucleotide polymorphisms (SNPs) in these 2 genes appeared strongly associated with tropical environment adaptation. In addition, the NADH:ubiquinone oxidoreductase subunit S4 (NDUFS4) gene exhibited a strong signature for positive selection in the LDC genome, and SNPs with marked allele frequency differences indicated a significant relationship with frigid environment adaptation.

CONCLUSION

Our findings partially clarify how selection footprints from environmental temperature stress can lead to advantageous genomic adaptions to tropical and frigid environments in poultry and provide a valuable resource for selective breeding of chickens.

Assessing the genomic diversity and relatedness in 10 Canadian heritage chicken lines using whole-genome sequence data

In the past 50 years, there has been a steep increase in the demand for poultry products, met by increasing production along with genetic selection for improved growth, efficiency, health and reproduction. The selection tends to reduce the number and type of genetic resources contributing to the majority of production. The University of Alberta maintains 10 heritage chicken lines (Brown Leghorn (BL), Light Sussex (LS), New Hampshire (NH), Saskatchewan Barred Rock (SaskBR), Shaver Barred Rock (ShaverBR), Shaver Rhode Island Red (RIR), White Leghorn (WL) and three commercial crosses called Alberta Meat Control strains 1957 (AMC-1957), 1978 sire line (AMC-1978-20S) and 1978 dam line (AMC-1978-30D), that played a large role in the evolution of the poultry industry in Canada. Since these lines have not been subjected to the same intensive selection pressures as commercial counterparts, they may contain unique genetic variants lost in commercial lines. Thus, for conservation management of these lines, the first step is to assess their genetic diversity. 71 male samples from across 10 lines were analysed using whole-genome sequencing and patterns of genetic diversity and relatedness among these lines were explored. AMC-1978-30D showed the highest genetic diversity as reflected in observed and expected heterozygosity (0.327 and 0.250), percentage of polymorphic markers (~ 65%) and average recent inbreeding coefficient (-0.039), followed by AMC-1978-20S and AMC-1957. BL showed the lowest genetic diversity as reflected in observed and expected heterozygosity (0.130 and 0.116), percentage of polymorphic markers (~31%) and average recent inbreeding coefficient (0.577), followed by LS, WL and NH. Our findings highlight the need for special attention for the populations of BL, WL, LS and NH, with the largest levels of inbreeding. Our results can be used to develop a breeding strategy to optimize and conserve the genetic variation present in heritage lines.

Recovering High-Quality Host Genomes from Gut Metagenomic Data through Genotype Imputation

Metagenomic datasets of host-associated microbial communities often contain host DNA that is usually discarded because the amount of data is too low for accurate host genetic analyses. However, genotype imputation can be employed to reconstruct host genotypes if a reference panel is available. Here, the performance of a two-step strategy is tested to impute genotypes from four types of reference panels built using different strategies to low-depth host genome data (≈2× coverage) recovered from intestinal samples of two chicken genetic lines. First, imputation accuracy is evaluated in 12 samples for which both low- and high-depth sequencing data are available, obtaining high imputation accuracies for all tested panels (>0.90). Second, the impact of reference panel choice in population genetics statistics on 100 chickens is assessed, all four panels yielding comparable results. In light of the observations, the feasibility and application of the applied imputation strategy are discussed for different species with regard to the host DNA proportion, genomic diversity, and availability of a reference panel. This method enables leveraging insofar discarded host DNA to get insights into the genetic structure of host populations, and in doing so, facilitates the implementation of hologenomic approaches that jointly analyze host and microbial genomic data.

Avian leukosis virus (ALV) is highly prevalent in fancy-chicken flocks in Saxony

The current prevalence of avian leukosis virus (ALV) in fancy chickens in Germany is unknown. Therefore, 537 cloacal swabs from 50 purebred fancy-chicken flocks in Saxony were tested for the presence of the ALV p27 protein using a commercial antigen-capture ELISA. The detection rate was 28.7% at the individual-animal level and 56.0% at the flock level. Phylogenetic analysis of PCR products obtained from 22 different flocks revealed the highest similarity to ALV subtype K. When classifying breeds by their origin, ALV detection rates differed significantly. Evaluation of questionnaire data revealed no significant differences between ALV-positive and negative flocks regarding mortality.

Assessment of linkage disequilibrium patterns between structural variants and single nucleotide polymorphisms in three commercial chicken populations

BACKGROUND

Structural variants (SV) are causative for some prominent phenotypic traits of livestock as different comb types in chickens or color patterns in pigs. Their effects on production traits are also increasingly studied. Nevertheless, accurately calling SV remains challenging. It is therefore of interest, whether close-by single nucleotide polymorphisms (SNPs) are in strong linkage disequilibrium (LD) with SVs and can serve as markers. Literature comes to different conclusions on whether SVs are in LD to SNPs on the same level as SNPs to other SNPs. The present study aimed to generate a precise SV callset from whole-genome short-read sequencing (WGS) data for three commercial chicken populations and to evaluate LD patterns between the called SVs and surrounding SNPs. It is thereby the first study that assessed LD between SVs and SNPs in chickens.

RESULTS

The final callset consisted of 12,294,329 bivariate SNPs, 4,301 deletions (DEL), 224 duplications (DUP), 218 inversions (INV) and 117 translocation breakpoints (BND). While average LD between DELs and SNPs was at the same level as between SNPs and SNPs, LD between other SVs and SNPs was strongly reduced (DUP: 40%, INV: 27%, BND: 19% of between-SNP LD). A main factor for the reduced LD was the presence of local minor allele frequency differences, which accounted for 50% of the difference between SNP - SNP and DUP - SNP LD. This was potentially accompanied by lower genotyping accuracies for DUP, INV and BND compared with SNPs and DELs. An evaluation of the presence of tag SNPs (SNP in highest LD to the variant of interest) further revealed DELs to be slightly less tagged by WGS SNPs than WGS SNPs by other SNPs. This difference, however, was no longer present when reducing the pool of potential tag SNPs to SNPs located on four different chicken genotyping arrays.

CONCLUSIONS

The results implied that genomic variance due to DELs in the chicken populations studied can be captured by different SNP marker sets as good as variance from WGS SNPs, whereas separate SV calling might be advisable for DUP, INV, and BND effects.

Genome-wide run of homozygosity analysis reveals candidate genomic regions associated with environmental adaptations of Tibetan native chickens

Background

In Tibet, the two most important breeds are Tibetan chicken and Lhasa white chicken, and the duo exhibit specific adaptations to the high altitude thereby supplying proteins for humans living in the plateau. These breeds are partly included in the conservation plans because they represent important chicken genetic resources. However, the genetic diversity of these chickens is rarely investigated. Based on whole-genome sequencing data of 113 chickens from 4 populations of Tibetan chicken including Shigatse (SH), Nyemo (NM), Dagze (DZ) and Nyingchi (LZ), as well as Lhasa white (LW) chicken breed, we investigated the genetic diversity of these chicken breeds by genetic differentiation, run of homozygosity (ROH), genomic inbreeding and selection signature analyses.

Results

Our results revealed high genetic diversity across the five chicken populations. The linkage disequilibrium decay was highest in LZ, while subtle genetic differentiation was found between LZ and other populations (Fst ranging from 0.05 to 0.10). Furthermore, the highest ROH-based inbreeding estimate (F_ROH) of 0.11 was observed in LZ. In other populations, the F_ROH ranged from 0.04 to 0.06. In total, 74, 111, 62, 42 and 54 ROH islands containing SNPs ranked top 1% for concurrency were identified in SH, NM, DZ, LZ and LW, respectively. Genes common to the ROH islands in the five populations included BDNF, CCDC34, LGR4, LIN7C, GLS, LOC101747789, MYO1B, STAT1 and STAT4. This suggested their essential roles in adaptation of the chickens. We also identified a common candidate genomic region harboring AMY2A, NTNG1 and VAV3 genes in all populations. These genes had been implicated in digestion, neurite growth and high-altitude adaptation.

Conclusions

High genetic diversity is observed in Tibetan native chickens. Inbreeding is more intense in the Nyingchi population which is also genetically distant from other chicken populations. Candidate genes in ROH islands are likely to be the drivers of adaptation to high altitude exhibited by the five Tibetan native chicken populations. Our findings contribute to the understanding of genetic diversity offer valuable insights for the genetic mechanism of adaptation, and provide veritable tools that can help in the design and implementation of breeding and conservation strategies for Tibetan native chickens.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08280-z.

Increasing calling accuracy, coverage, and read-depth in sequence data by the use of haplotype blocks

High-throughput genotyping of large numbers of lines remains a key challenge in plant genetics, requiring geneticists and breeders to find a balance between data quality and the number of genotyped lines under a variety of different existing genotyping technologies when resources are limited. In this work, we are proposing a new imputation pipeline (“HBimpute”) that can be used to generate high-quality genomic data from low read-depth whole-genome-sequence data. The key idea of the pipeline is the use of haplotype blocks from the software HaploBlocker to identify locally similar lines and subsequently use the reads of all locally similar lines in the variant calling for a specific line. The effectiveness of the pipeline is showcased on a dataset of 321 doubled haploid lines of a European maize landrace, which were sequenced at 0.5X read-depth. The overall imputing error rates are cut in half compared to state-of-the-art software like BEAGLE and STITCH, while the average read-depth is increased to 83X, thus enabling the calling of copy number variation. The usefulness of the obtained imputed data panel is further evaluated by comparing the performance of sequence data in common breeding applications to that of genomic data generated with a genotyping array. For both genome-wide association studies and genomic prediction, results are on par or even slightly better than results obtained with high-density array data (600k). In particular for genomic prediction, we observe slightly higher data quality for the sequence data compared to the 600k array in the form of higher prediction accuracies. This occurred specifically when reducing the data panel to the set of overlapping markers between sequence and array, indicating that sequencing data can benefit from the same marker ascertainment as used in the array process to increase the quality and usability of genomic data.

High-throughput genotyping of large numbers of lines remains a key challenge in plant genetics and breeding. Cost, precision, and throughput must be balanced to achieve optimal efficiency given available technologies and finite resources. Although genotyping arrays are still considered the gold standard in high-throughput quantitative genetics, recent advances in sequencing provide new opportunities. Both the quality and cost of genomic data generated based on sequencing are highly dependent on the used read-depth. In this work, we propose a new imputation pipeline (“HBimpute”) that uses haplotype blocks to detect individuals of the same genetic origin and subsequently uses all reads of those individuals in the variant calling. Thus, the obtained virtual read-depth is artificially increased, leading to higher calling accuracy, coverage, and the ability to call copy number variation based on low read-depth sequencing data. To conclude, our approach makes sequencing a cost-competitive alternative to genotyping arrays with the added benefit of allowing the calling of structural variation.

Genetic Diversity of 17 Autochthonous Italian Chicken Breeds and Their Extinction Risk Status

The preservation of genetic variability of autochthonous poultry breeds is crucial in global biodiversity. A recent report revealed small breed size and potential risk of extinction of all native Italian poultry breeds; therefore, a correct assessment of their genetic diversity is necessary for a suitable management of their preservation. In this work, we provided an overview of the contribution to poultry biodiversity of some Italian autochthonous breeds reared in conservation centers devoted to local biodiversity preservation. The level of genetic diversity, molecular kinship, inbreeding, contribution to overall genetic diversity, and rate of extinction of each breed were analyzed with a set of 14 microsatellite loci in 17 autochthonous chicken breeds. To evaluate genetic variability, total number (Na), and effective number (Ne) of alleles, observed (Ho) and expected (He) heterozygosity, and F (Wright’s inbreeding coefficient) index were surveyed. The contribution of each analyzed breed to genetic diversity of the whole dataset was assessed using MolKin3.0; global genetic diversity and allelic richness contributions were evaluated. All the investigated loci were polymorphic; 209 alleles were identified (94 of which private alleles). The average number of alleles per locus was 3.62, and the effective number of alleles was 2.27. The Ne resulted lower in all breeds due to the presence of low-frequency alleles that can be easily lost by genetic drift, thus reducing the genetic variability of the breeds, and increasing their risk of extinction. The global molecular kinship was 27%, the average breed molecular kinship was 53%, and the mean inbreeding rate 43%, with a self-coancestry of 78%. Wright’s statistical analysis showed a 41% excess of homozygous due to breed genetic differences (34%) and to inbreeding within the breed (9%). Genetic variability analysis showed that 11 breeds were in endangered status. The contribution to Italian poultry genetic diversity, estimated as global genetic diversity, and ranged from 30.2 to 98.5%. In conclusion, the investigated breeds maintain a unique genetic pattern and play an important role in global Italian poultry biodiversity, providing a remarkable contribution to genetic variability.

Comparison of major histocompatibility complex-B variability in Sri Lankan indigenous chickens with five global chicken populations using MHC-B SNP panel

In the present study, we investigated the major histocompatibility complex (MHC)-B haplotypes diversity of Sri Lankan indigenous chickens from three different geographical sites consisting of highly mixed populations using 90 SNPs in the MHC-B region. A total of 48 haplotypes were identified. Those included 37 novel haplotypes and 11 previously identified 'standard' haplotypes. The MHC-linked marker, LEI0258, had 23 alleles showing less diversity than defined by MHC-B SNP haplotypes. Among those identified haplotypes, five standard haplotypes-BSNP-O02, BSNP-M01, BSNP-A04, BSNP-K03, BSNP-T04-were most commonly observed, suggesting past introgression of imported breeds. Comparison of the MHC-B haplotypes of Sri Lankan and four other global populations with previously defined haplotypes indicated the sharing of 23 standard haplotypes with common origins. Novel haplotypes are population-specific and not shared among the geographical boundaries. Backyard indigenous chickens are unselected, highly crossbred, and generally thrive under dynamic environmental conditions. Hence free-range production systems may be responsible for maintaining high diversity in the MHC-B region with novel haplotypes.

Effects of the domestic thyroid stimulating hormone receptor (TSHR) variant on the hypothalamic-pituitary-thyroid axis and behavior in chicken

Domestic chickens are less fearful, have a faster sexual development, grow bigger, and lay more eggs than their primary ancestor, the red junglefowl. Several candidate genetic variants selected during domestication have been identified, but only a few studies have directly linked them with distinct phenotypic traits. Notably, a variant of the thyroid stimulating hormone receptor (TSHR) gene has been under strong positive selection over the past millennium, but it's function and mechanisms of action are still largely unresolved. We therefore assessed the abundance of the domestic TSHR variant and possible genomic selection signatures in an extensive data set comprising multiple commercial and village chicken populations as well as wild-living extant members of the genus Gallus. Furthermore, by mean of extensive backcrossing we introgressed the wild-type TSHR variant from red junglefowl into domestic White Leghorn chickens and investigated gene expression, hormone levels, cold adaptation, and behavior in chickens possessing either the wild-type or domestic TSHR variant. While the domestic TSHR was the most common variant in all studied domestic populations and in one of two red junglefowl population, it was not detected in the other Gallus species. Functionally, the individuals with the domestic TSHR variant had a lower expression of the TSHR in the hypothalamus and marginally higher in the thyroid gland than wild-type TSHR individuals. Expression of TSHB and DIO2, two regulators of sexual maturity and reproduction in birds, was higher in the pituitary gland of the domestic-variant chickens. Furthermore, the domestic variant was associated with higher activity in the open field test. Our findings confirm that the spread of the domestic TSHR variant is limited to domesticated chickens, and to a lesser extent, their wild counterpart, the red junglefowl. Furthermore, we showed that effects of genetic variability in TSHR mirror key differences in gene expression and behavior previously described between the red junglefowl and domestic chicken.

The link between broiler flock heterogeneity and cecal microbiome composition

BACKGROUND

Despite low genetic variation of broilers and deployment of considerate management practices, there still exists considerable body weight (BW) heterogeneity within broiler flocks which adversely affects the commercial value. The purpose of this study was to investigate the role of the cecal microbiome in weight differences between animals. Understanding how the gut microbiome may contribute to flock heterogeneity helps to pave the road for identifying methods to improve flock uniformity and performance.

RESULTS

Two hundred eighteen male broiler chicks were housed in the same pen, reared for 37 days, and at study end the 25 birds with highest BW (Big) and the 25 birds with lowest BW (Small) were selected for microbiome analysis. Cecal contents were analyzed by a hybrid metagenomic sequencing approach combining long and short read sequencing. We found that Big birds displayed higher microbial alpha diversity, higher microbiome uniformity (i.e. lower beta diversity within the group of Big birds), higher levels of SCFA-producing and health-associated bacterial taxa such as Lachnospiraceae, Faecalibacterium, Butyricicoccus and Christensenellales, and lower levels of Akkermansia muciniphila and Escherichia coli as compared to Small birds.

CONCLUSION

Cecal microbiome characteristics could be linked to the size of broiler chickens. Differences in alpha diversity, beta diversity and taxa abundances all seem to be directly associated with growth differences observed in an otherwise similar broiler flock.

The chicken pan-genome reveals gene content variation and a promoter region deletion in IGF2BP1 affecting body size

Domestication and breeding have reshaped the genomic architecture of chicken, but the retention and loss of genomic elements during these evolutionary processes remain unclear. We present the first chicken pan-genome constructed using 664 individuals, which identified an additional ∼66.5 Mb sequences that are absent from the reference genome (GRCg6a). The constructed pan-genome encoded 20,491 predicated protein-coding genes, of which higher expression level are observed in conserved genes relative to dispensable genes. Presence/absence variation (PAV) analyses demonstrated that gene PAV in chicken was shaped by selection, genetic drift, and hybridization. PAV-based GWAS identified numerous candidate mutations related to growth, carcass composition, meat quality, or physiological traits. Among them, a deletion in the promoter region of IGF2BP1 affecting chicken body size is reported, which is supported by functional studies and extra samples. This is the first time to report the causal variant of chicken body size QTL located at chromosome 27 which was repeatedly reported. Therefore, the chicken pan-genome is a useful resource for biological discovery and breeding. It improves our understanding of chicken genome diversity and provides materials to unveil the evolution history of chicken domestication.

Domestic chicken diversity: Origin, distribution, and adaptation

Chicken is the most numerous among the domesticated livestock species. Across cultures, religions, and societies, chicken is widely accepted with little or no taboo compared to other domestic animals. Its adaptability to diverse environmental conditions and demonstrated potential for breeding improvement provide a unique genetic resource for addressing the challenges of food security in a world impacted by climatic change and human population growth. Recent studies, shedding new knowledge on the chicken genomes, have helped reconstruct its past evolutionary history. Here, we review the literature concerning the origin, dispersion, and adaptation of domestic chicken. We highlight the role of human and natural selection in shaping the diversity of the species and provide a few examples of knowledge gaps that may be the focus of future research.

The Effects of Domestication on the Brain and Behavior of the Chicken in the Light of Evolution

The avian class is characterized by particularly strong variability in their domesticated species. With more than 250 breeds and highly efficient commercial lines, domestic chickens represent the outcome of a really long period of artificial selection. One characteristic of domestication is the alterations in brain size and brain composition. The influence of domestication on brain morphology has been reviewed in the past, but mostly with a focus on mammals. Studies on avian species have seldom been taken into account. In this review, we would like to give an overview about the changes and variations in (brain) morphology and behavior in the domestic chicken, taking into consideration the constraints of evolutionary theory and the sense or nonsense of excessive artificial selection.

Keel Bone Damage in Laying Hens-Its Relation to Bone Mineral Density, Body Growth Rate and Laying Performance

Simple Summary

Keel bone damage is an important welfare issue for laying hens. Four lines of laying hens, differing in phylogenetic origin and laying rate kept in single cages or a floor housing system were weighed and deformities of the keel bone were evaluated regularly between 15 and 69 weeks of age. Deformities, fractures and the bone mineral density of the keels were assessed after hens were euthanized. We analyzed the relationship between bone mineral density and total egg number as well as body growth. Hens kept in cages showed more deformities, but fewer fractures and a lower bone mineral density of the keel bone than did floor-housed hens. Keel bones of white-egg layers had a lower mineral density and were more often deformed compared with brown-egg layers. Keel bones were more often broken in hens of the layer lines with a high laying rate compared to the lines with a moderate laying rate. Laying rate and adult body weight had an effect on the keel bone mineral density. The study contributes to the understanding of factors causing keel bone damage in laying hens. It showed that the bone mineral density greatly affects keel bone deformities.

Abstract

Keel bone damage is an important animal welfare problem in laying hens. Two generations of four layer lines, differing in phylogenetic background and performance level and kept in single cages or floor pens were weighed and scored for keel bone deformities (KBD) during the laying period. KBD, keel bone fractures (KBF) and the bone mineral density (BMD) of the keels were assessed post mortem. For BMD, relationships to laying performance and body growth were estimated. Caged hens showed more deformities, but fewer fractures and a lower BMD of the keel bone than floor-housed hens. White-egg layers had a lower BMD (0.140–0.165 g/cm²) and more KBD than brown-egg layers (0.179–0.184 g/cm²). KBF occurred more often in the high-performing lines than the moderate-performing ones. However, in the high-performing lines, BMD was positively related to total egg number from 18 to 29 weeks of age. The adult body weight derived from fitted growth curves (Gompertz function) had a significant effect (p < 0.001) on keels’ BMD. The study contributes to the understanding of predisposing factors for keel bone damage in laying hens. It showed that the growth rate has a rather subordinate effect on keels’ BMD, while the BMD itself greatly affects KBD.

How imputation can mitigate SNP ascertainment Bias

Background

Population genetic studies based on genotyped single nucleotide polymorphisms (SNPs) are influenced by a non-random selection of the SNPs included in the used genotyping arrays. The resulting bias in the estimation of allele frequency spectra and population genetics parameters like heterozygosity and genetic distances relative to whole genome sequencing (WGS) data is known as SNP ascertainment bias. Full correction for this bias requires detailed knowledge of the array design process, which is often not available in practice. This study suggests an alternative approach to mitigate ascertainment bias of a large set of genotyped individuals by using information of a small set of sequenced individuals via imputation without the need for prior knowledge on the array design.

Results

The strategy was first tested by simulating additional ascertainment bias with a set of 1566 chickens from 74 populations that were genotyped for the positions of the Affymetrix Axiom™ 580 k Genome-Wide Chicken Array. Imputation accuracy was shown to be consistently higher for populations used for SNP discovery during the simulated array design process. Reference sets of at least one individual per population in the study set led to a strong correction of ascertainment bias for estimates of expected and observed heterozygosity, Wright’s Fixation Index and Nei’s Standard Genetic Distance. In contrast, unbalanced reference sets (overrepresentation of populations compared to the study set) introduced a new bias towards the reference populations. Finally, the array genotypes were imputed to WGS by utilization of reference sets of 74 individuals (one per population) to 98 individuals (additional commercial chickens) and compared with a mixture of individually and pooled sequenced populations. The imputation reduced the slope between heterozygosity estimates of array data and WGS data from 1.94 to 1.26 when using the smaller balanced reference panel and to 1.44 when using the larger but unbalanced reference panel. This generally supported the results from simulation but was less favorable, advocating for a larger reference panel when imputing to WGS.

Conclusions

The results highlight the potential of using imputation for mitigation of SNP ascertainment bias but also underline the need for unbiased reference sets.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07663-6.

Researching on the fine structure and admixture of the worldwide chicken population reveal connections between populations and important events in breeding history

Here, we have evaluated the general genomic structure and diversity and studied the divergence resulting from selection and historical admixture events for a collection of worldwide chicken breeds. In total, 636 genomes (43 populations) were sequenced from chickens of American, Chinese, Indonesian, and European origin. Evaluated populations included wild junglefowl, rural indigenous chickens, breeds that have been widely used to improve modern western poultry populations and current commercial stocks bred for efficient meat and egg production. In‐depth characterizations of the genome structure and genomic relationships among these populations were performed, and population admixture events were investigated. In addition, the genomic architectures of several domestication traits and central documented events in the recent breeding history were explored. Our results provide detailed insights into the contributions from population admixture events described in the historical literature to the genomic variation in the domestic chicken. In particular, we find that the genomes of modern chicken stocks used for meat production both in eastern (Asia) and western (Europe/US) agriculture are dominated by contributions from heavy Asian breeds. Further, by exploring the link between genomic selective divergence and pigmentation, connections to functional genes feather coloring were confirmed.

Influence of Age and Phylogenetic Background on Blood Parameters Associated With Bone Metabolism in Laying Hens

The high laying performance of today’s laying hens places enormous demands on their mineral metabolism. While up-to-date data are rare, the present study aimed to describe blood parameters associated with egg laying and bone metabolism during the pre-laying period, in the course of the laying period and the daily egg laying cycle. Ten to 15 laying hens of two high-performing, phylogenetically divergent lines (BLA: brown-egg layer; WLA: white-egg layer), kept in single cages were blood sampled at 17, 25, 29, 49, and 69 weeks of age. Sampling was made at 6 a.m., 10 a.m., 2 p.m. and, with the exception of week 17, 6 p.m. Blood samples were analyzed for concentrations of total and ionized calcium, inorganic phosphate (PO4), markers of bone formation (osteocalcin) and resorption [carboxyterminal crosslinked telopeptide of type I collagen (CTX-I)], 25-hydroxycholecalciferol (25(OH)D₃) and estradiol-17β. In the pre-laying period (17 week), the estradiol-17β level calculated for WLA was more than twice as high as the level calculated for BLA, while no significant difference could be observed in the laying period (25 to 69 weeks). BLA hens had significantly higher total calcium concentrations at 49 weeks of age as well as up to twice as high levels of osteocalcin and 25(OH)D₃ than WLA at any time of the day from 25 to 69 weeks of age. While osteocalcin, CTX-I and 25(OH)D₃ concentrations were significantly higher before the onset of lay, total calcium and estradiol-17β levels significantly increased from 17 to 69 weeks of age. In contrast, PO4 values varied only slightly during the experimental period and ionized calcium was highest at 17 and 49 weeks of age and lowest around peak production (29 week). In the course of the daily egg laying cycle blood concentrations clearly reflected the stage of egg formation. Our results provide up-to-date data of bone- and egg laying-associated blood parameters of two contemporary purebred layer lines over the course of the pre- and egg-laying period and the daily egg laying cycle. Differences between brown- and white-egg layers raise questions, whether phylogenetic background determines their efficiency to cope with high calcium demands relating to egg production.

Genetic diversity in global chicken breeds in relation to their genetic distances to wild populations

Background

Migration of a population from its founder population is expected to cause a reduction of its genetic diversity and facilitates differentiation between the population and its founder population, as predicted by the theory of genetic isolation by distance. Consistent with that theory, a model of expansion from a single founder predicts that patterns of genetic diversity in populations can be explained well by their geographic expansion from their founders, which is correlated with genetic differentiation.

Methods

To investigate this in chicken, we estimated the relationship between the genetic diversity of 160 domesticated chicken populations and their genetic distances to wild chicken populations.

Results

Our results show a strong inverse relationship, i.e. 88.6% of the variation in the overall genetic diversity of domesticated chicken populations was explained by their genetic distance to the wild populations. We also investigated whether the patterns of genetic diversity of different types of single nucleotide polymorphisms (SNPs) and genes are similar to that of the overall genome. Among the SNP classes, the non-synonymous SNPs deviated most from the overall genome. However, genetic distance to the wild chicken still explained more variation in domesticated chicken diversity across all SNP classes, which ranged from 83.0 to 89.3%.

Conclusions

Genetic distance between domesticated chicken populations and their wild relatives can predict the genetic diversity of the domesticated populations. On the one hand, genes with little genetic variation across populations, regardless of the genetic distance to the wild population, are associated with major functions such as brain development. Changes in such genes may be detrimental to the species. On the other hand, genetic diversity seems to change at a faster rate within genes that are associated with e.g. protein transport and protein and lipid metabolic processes. In general, such genes may be flexible to changes according to the populations’ needs. These results contribute to the knowledge of the evolutionary patterns of different functional genomic regions in the chicken.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00628-z.

How array design creates SNP ascertainment bias

Single nucleotide polymorphisms (SNPs), genotyped with arrays, have become a widely used marker type in population genetic analyses over the last 10 years. However, compared to whole genome re-sequencing data, arrays are known to lack a substantial proportion of globally rare variants and tend to be biased towards variants present in populations involved in the development process of the respective array. This affects population genetic estimators and is known as SNP ascertainment bias. We investigated factors contributing to ascertainment bias in array development by redesigning the Axiom^™ Genome-Wide Chicken Array in silico and evaluating changes in allele frequency spectra and heterozygosity estimates in a stepwise manner. A sequential reduction of rare alleles during the development process was shown. This was mainly caused by the identification of SNPs in a limited set of populations and a within-population selection of common SNPs when aiming for equidistant spacing. These effects were shown to be less severe with a larger discovery panel. Additionally, a generally massive overestimation of expected heterozygosity for the ascertained SNP sets was shown. This overestimation was 24% higher for populations involved in the discovery process than not involved populations in case of the original array. The same was observed after the SNP discovery step in the redesign. However, an unequal contribution of populations during the SNP selection can mask this effect but also adds uncertainty. Finally, we make suggestions for the design of specialized arrays for large scale projects where whole genome re-sequencing techniques are still too expensive.

High Levels of Genetic Variation in MHC-Linked Microsatellite Markers from Native Chicken Breeds

The major histocompatibility complex (MHC) is a highly polymorphic gene region that regulates cellular communication in all specific immune responses. In this study, we investigated 11 microsatellite (MS) markers in the MHC-B region of chicken populations from four countries: Sri Lanka, Bangladesh, South Korea, and Nigeria. The MS markers were divided into two sets. Set 1 included five novel MS markers, which we assessed using 192 samples from 21 populations. Set 2 included six previously reported markers, which we assessed using 881 samples from 29 populations. The Set 1 MS markers had lower polymorphism (polymorphic information content (PIC) < 0.5) than the Set 2 markers (PIC = 0.4-0.9). In all populations, the LEI0258 marker was the most polymorphic, with a total of 38 alleles (PIC = 0.912, expected heterozygosity (He) = 0.918). Local populations from Sri Lanka, Bangladesh, and Nigeria had higher allele diversity and more haplotypes for Set 2 MS markers than Korean and commercial populations. The Sri Lankan Karuwalagaswewa village population had the highest MHC diversity (mean allele number = 8.17, He = 0.657), whereas the white leghorn population had the lowest (mean allele number = 2.33, He = 0.342). A total of 409 haplotypes (89 shared and 320 unique), with a range of 4 (Rhode Island red) to 46 (Karuwalagaswewa village (TA)), were identified. Among the shared haplotypes, the B21-like haplotype was identified in 15 populations. The genetic relationship observed in a neighbour-joining tree based on the DA distance agreed with the breeding histories and geographic separations. The results indicated high MHC diversity in the local chicken populations. The difference in the allelic pattern among populations presumably reflects the effects of different genotypes, environments, geographic variation, and breeding policies in each country. The selection of MHC allele in domestic poultry can vary due to intensification of poultry production. Preserved MHC diversity in local chicken provides a great opportunity for future studies that address the relationships between MHC polymorphisms and differential immune responses.

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.

Unraveling the history of the genus Gallus through whole genome sequencing

The genus Gallus is distributed across a large part of Southeast Asia and has received special interest because the domestic chicken, Gallus gallus domesticus, has spread all over the world and is a major protein source for humans. There are four species: the red junglefowl (G. gallus), the green junglefowl (G. varius), the Lafayette's junglefowl (G. lafayettii) and the grey junglefowl (G. sonneratii). The aim of this study is to reconstruct the history of these species by a whole genome sequencing approach and resolve inconsistencies between well supported topologies inferred using different data and methods. Using deep sequencing, we identified over 35 million SNPs and reconstructed the phylogeny of the Gallus genus using both distance (BioNJ) and maximum likelihood (ML) methods. We observed discrepancies according to reconstruction methods and genomic components. The two most supported topologies were previously reported and were discriminated by using phylogenetic and gene flow analyses, based on ABBA statistics. Terminology fix requested by the deputy editor led to support a scenario with G. gallus as the earliest branching lineage of the Gallus genus, instead of G. varius. We discuss the probable causes for the discrepancy. A likely one is that G. sonneratii samples from parks or private collections are all recent hybrids, with roughly 10% of their autosomal genome originating from G. gallus. The removal of those regions is needed to provide reliable data, which was not done in previous studies. We took care of this and additionally included two wild G. sonneratii samples from India, showing no trace of introgression. This reinforces the importance of carefully selecting and validating samples and genomic components in phylogenomics.

Genome-Wide Population Genetic Analysis of Commercial, Indigenous, Game, and Wild Chickens Using 600K SNP Microarray Data

Following chicken domestication, diversified chicken breeds were developed by both natural and artificial selection, which led to the accumulation of abundant genetic and phenotypic variations, making chickens an ideal genetic research model. To better understand the genetic structure of chicken breeds under different selection pressures, we genotyped various chicken populations with specific selection targets, including indigenous, commercial, gamecock, and wild ancestral chickens, using the 600K SNP array. We analyzed the population structure, genetic relationships, run of homozygosity (ROH), effective population number (Ne), and other genetic parameters. The wild ancestral population, red junglefowl (RJF), possessed the highest diversity, in comparison with all other domesticated populations, which was supported by linkage disequilibrium decay (LD), effective population number, and ROH analyses. The gamecock breeds, which were subjected to stronger male-biased selection for fighting-related traits, also presented higher variation than the commercial and indigenous breeds. Admixture analysis also indicated that game breed is a relatively independent branch of Chinese local breeds. Following intense selection for reproductive and productive traits, the commercial lines showed the least diversity. We also observed that the European local chickens had lower genetic variation than the Chinese local breeds, which could be attributed to the shorter history of the European breed. ROH were present in a breed specific manner and 191 ROH island were detected on four groups (commercial, local, game and wild chickens). These ROH islands were involved in egg production, growth and silky feathers and other traits. Moreover, we estimated the effective sex ratio of these breeds to demonstrate the change in the ratio of the two sexes. We found that commercial chickens had a greater sex imbalance between females and males. The commercial lines showed the highest female-to-male ratios. Interestingly, RJF comprised a greater proportion of males than females. Our results show the population genetics of chickens under selection pressures, and can aid in the development of better conservation strategies for different chicken breeds.

Bone Health or Performance? Adaptation Response of Genetically Divergent Chicken Layer Lines to a Nutritive Calcium Depletion

Simple Summary

Selection for high egg production in laying hens has led to an increased calcium requirement and consequently to an intensified calcium mobilization from the bones. However, excessive demineralization can lead to osteoporosis, which is manifested by a high incidence of bone-damaged hens. In this study, we characterized the adaptation response of laying hens to a repeated dietary calcium restriction (1.09% instead of 4.26% calcium) by means of egg production, eggshell quality, body weight and bone stability. The animal model included four layer lines differing in performance level (high vs. moderately performing lines) and phylogenetic origin (white-egg vs. brown-egg layers). We assumed that the high performing lines would respond by maintaining egg production level at the expense of eggshell quality and bone health. Egg production and eggshell quality declined considerably and bone demineralization occurred in all lines. Contrary to our hypothesis, there was evidence that phylogeny rather than performance level influenced the hens′ response. The brown-egg lines appeared to be more tolerant to the calcium depletion, while the white-egg lines were more sensitive. Our findings demonstrate the influence of genetics on the adaptive capacity of chickens and underline the importance of preserving genetic variability to cope with potential future environmental challenges.

Abstract

In modern laying hybrids, calcium (Ca) homeostasis is immensely challenged by daily eggshell calcification. However, excessive mobilization of Ca from bones may lead to osteoporosis, which then manifests in a high incidence of poor bone quality. The aim of this study was to characterize the hens’ adaptation response to an alternating dietary Ca restriction. The animal model consisted of four purebred layer lines, differing in laying performance (high vs. moderately performing lines) and phylogenetic origin (white- vs. brown-egg lines). According to the resource allocation theory, hens selected for high egg production were assumed to show a different response pattern to cope with this nutritive challenge compared to moderately performing lines. Data collected included egg number, egg quality traits, body weight and bone characteristics. The Ca depletion led to a temporary drop in egg production and shell quality and a loss of bone stability due to Ca mobilization. The white-egg lines response was more pronounced, whereas the brown-egg lines were less sensitive towards reduced Ca supply. Our study shows that the hens’ responsiveness to coping with a nutritive Ca depletion is not ultimately linked to genetic selection for increased egg production but rather to phylogenetic origin.

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.