Genetic Variability in Local and Imported Germplasm Chicken Populations as Revealed by Analyzing Runs of Homozygosity

Genetic analysis of a Kaijiang duck conservation population through genome-wide scan

1. The Kaijiang duck is a native Chinese breed known for its excellent egg laying performance, killing-out percentage (88.57%), and disease resistance. The assessment of population genetic structure is the basis for understanding the genetics of indigenous breeds and for their protection and management.2. In this study, whole-genome sequencing was performed on 60 Kaijiang ducks to identify genetic variations and investigate the population structure. Homozygosity (ROH) analysis was conducted to assess inbreeding levels in the population.3. The study revealed a moderate level of inbreeding, indicated by an average inbreeding coefficient of 0.1043. This may impact the overall genetic diversity.4. Genomic Regions of Interest identified included 168 genomic regions exhibiting high levels of autozygosity. These regions were associated with processes including muscle growth, pigmentation, neuromodulation, and growth and reproduction.5. The significance of these pathways indicated their potential role in shaping the desirable traits of the Kaijiang duck. These findings provide insights into the genetic basis of the Kaijiang duck's desirable traits and can inform future breeding and conservation efforts.

Dissecting Selective Signatures and Candidate Genes in Grandparent Lines Subject to High Selection Pressure for Broiler Production and in a Local Russian Chicken Breed of Ushanka

Breeding improvements and quantitative trait genetics are essential to the advancement of broiler production. The impact of artificial selection on genomic architecture and the genetic markers sought remains a key area of research. Here, we used whole-genome resequencing data to analyze the genomic architecture, diversity, and selective sweeps in Cornish White (CRW) and Plymouth Rock White (PRW) transboundary breeds selected for meat production and, comparatively, in an aboriginal Russian breed of Ushanka (USH). Reads were aligned to the reference genome bGalGal1.mat.broiler.GRCg7b and filtered to remove PCR duplicates and low-quality reads using BWA-MEM2 and bcftools software; 12,563,892 SNPs were produced for subsequent analyses. Compared to CRW and PRW, USH had a lower diversity and a higher genetic distinctiveness. Selective sweep regions and corresponding candidate genes were examined based on ZFST, hapFLK, and ROH assessment procedures. Twenty-seven prioritized chicken genes and the functional projection from human homologs suggest their importance for selection signals in the studied breeds. These genes have a functional relationship with such trait categories as body weight, muscles, fat metabolism and deposition, reproduction, etc., mainly aligned with the QTLs in the sweep regions. This information is pivotal for further executing genomic selection to enhance phenotypic traits.

Large-scale genome-wide SNP analysis reveals the rugged (and ragged) landscape of global ancestry, phylogeny, and demographic history in chicken breeds

The worldwide chicken gene pool encompasses a remarkable, but shrinking, number of divergently selected breeds of diverse origin. This study was a large-scale genome-wide analysis of the landscape of the complex molecular architecture, genetic variability, and detailed structure among 49 populations. These populations represent a significant sample of the world's chicken breeds from Europe (Russia, Czech Republic, France, Spain, UK, etc.), Asia (China), North America (USA), and Oceania (Australia). Based on the results of breed genotyping using the Illumina 60K single nucleotide polymorphism (SNP) chip, a bioinformatic analysis was carried out. This included the calculation of heterozygosity/homozygosity statistics, inbreeding coefficients, and effective population size. It also included assessment of linkage disequilibrium and construction of phylogenetic trees. Using multidimensional scaling, principal component analysis, and ADMIXTURE-assisted global ancestry analysis, we explored the genetic structure of populations and subpopulations in each breed. An overall 49-population phylogeny analysis was also performed, and a refined evolutionary model of chicken breed formation was proposed, which included egg, meat, dual-purpose types, and ambiguous breeds. Such a large-scale survey of genetic resources in poultry farming using modern genomic methods is of great interest both from the viewpoint of a general understanding of the genetics of the domestic chicken and for the further development of genomic technologies and approaches in poultry breeding. In general, whole genome SNP genotyping of promising chicken breeds from the worldwide gene pool will promote the further development of modern genomic science as applied to poultry.

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.

Whole Genome Screening Procures a Holistic Hold of the Russian Chicken Gene Pool Heritage and Demographic History

A study for genomic variation that may reflect putative selective signaling and be associated with economically important traits is instrumental for obtaining information about demographic and selection history in domestic animal species and populations. A rich variety of the Russian chicken gene pool breeds warrants a further detailed study. Specifically, their genomic features can derive implications from their genome architecture and selective footprints for their subsequent breeding and practical efficient exploitation. In the present work, whole genome genotyping of 19 chicken breeds (20 populations with up to 71 samples each) was performed using the Chicken 50 K BeadChip DNA chip. The studied breed sample included six native Russian breeds of chickens developed in the 17th-19th centuries, as well as eight Russian chicken breeds, including the Russian White (RW), created in the 20th century on the basis of improving local chickens using breeds of foreign selection. Five specialized foreign breeds of chickens, including the White Leghorn (WL), were used along with other breeds representing the Russian gene pool. The characteristics of the genetic diversity and phylogenetic relationships of the native breeds of chickens were represented in comparison with foreign breeds. It was established that the studied native breeds demonstrate their own genetic structure that distinguishes them from foreign breeds, and from each other. For example, we previously made an assumption on what could cause the differences between two RW populations, RW1 and RW2. From the data obtained here, it was verified that WL was additionally crossed to RW2, unlike RW1. Thus, inherently, RW1 is a purer population of this improved Russian breed. A significant contribution of the gene pool of native breeds to the global genetic diversity of chickens was shown. In general, based on the results of a multilateral survey of this sample of breeds, it can be concluded that phylogenetic relationships based on their genetic structure and variability robustly reflect the known, previously postulated and newly discovered patterns of evolution of native chickens. The results herein presented will aid selection and breeding work using this gene pool.

Selective footprints and genes relevant to cold adaptation and other phenotypic traits are unscrambled in the genomes of divergently selected chicken breeds

BACKGROUND

The genomes of worldwide poultry breeds divergently selected for performance and other phenotypic traits may also be affected by, and formed due to, past and current admixture events. Adaptation to diverse environments, including acclimation to harsh climatic conditions, has also left selection footprints in breed genomes.

RESULTS

Using the Chicken 50K_CobbCons SNP chip, we genotyped four divergently selected breeds: two aboriginal, cold tolerant Ushanka and Orloff Mille Fleur, one egg-type Russian White subjected to artificial selection for cold tolerance, and one meat-type White Cornish. Signals of selective sweeps were determined in the studied breeds using three methods: (1) assessment of runs of homozygosity islands, (2) F_ST based population differential analysis, and (3) haplotype differentiation analysis. Genomic regions of true selection signatures were identified by two or more methods or in two or more breeds. In these regions, we detected 540 prioritized candidate genes supplemented them with those that occurred in one breed using one statistic and were suggested in other studies. Amongst them, SOX5, ME3, ZNF536, WWP1, RIPK2, OSGIN2, DECR1, TPO, PPARGC1A, BDNF, MSTN, and beta-keratin genes can be especially mentioned as candidates for cold adaptation. Epigenetic factors may be involved in regulating some of these important genes (e.g., TPO and BDNF).

CONCLUSION

Based on a genome-wide scan, our findings can help dissect the genetic architecture underlying various phenotypic traits in chicken breeds. These include genes representing the sine qua non for adaptation to harsh environments. Cold tolerance in acclimated chicken breeds may be developed following one of few specific gene expression mechanisms or more than one overlapping response known in cold-exposed individuals, and this warrants further investigation.

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.

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.

Identification of Key Candidate Genes in Runs of Homozygosity of the Genome of Two Chicken Breeds, Associated with Cold Adaptation

Simple Summary

The search for genomic regions related to adaptive abilities preserved in the chicken gene pool of two breeds, which have not been under intensive selection pressure, is of great importance for breeding in the future. This study aimed to identify key candidate genes associated with the adaptation of chickens to cold environments (using the example of the Russian White breed) by using molecular genetic methods. A total of 12 key genes on breed-specific ROH (runs of homozygosity) islands were identified, which may be potential candidate genes associated with the high level of adaptability of chickens to cold environments in the early postnatal period. These genes were associated with lipid metabolism, maintaining body temperature in cold environments, non-shivering thermogenesis and muscle development and are perspectives for further research.

Abstract

It is well known that the chicken gene pools have high adaptive abilities, including adaptation to cold environments. This research aimed to study the genomic distribution of runs of homozygosity (ROH) in a population of Russian White (RW) chickens as a result of selection for adaptation to cold environments in the early postnatal period, to perform a structural annotation of the discovered breed-specific regions of the genome (compared to chickens of the Amroks breed) and to suggest key candidate genes associated with the adaptation of RW chickens to cold environments. Genotyping of individual samples was performed using Illumina Chicken 60K SNP BeadChip^® chips. The search for homozygous regions by individual chromosomes was carried out using the PLINK 1.9 program and the detectRuns R package. Twelve key genes on breed-specific ROH islands were identified. They may be considered as potential candidate genes associated with the high adaptive ability of chickens in cold environments in the early postnatal period. Genes associated with lipid metabolism (SOCS3, NDUFA4, TXNRD2, IGFBP 1, IGFBP 3), maintaining body temperature in cold environments (ADIPOQ, GCGR, TRPM2), non-shivering thermogenesis (RYR2, CAMK2G, STK25) and muscle development (METTL21C) are perspectives for further research. This study contributes to our understanding of the mechanisms of adaptation to cold environments in chickens and provides a molecular basis for selection work.

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.

Analysis of Autozygosity Using Whole-Genome Sequence Data of Full-Sib Families in Pikeperch (Sander lucioperca)

Pikeperch (Sander lucioperca) has emerged as a high value species to the aquaculture industry. However, its farming techniques are at an early stage and its production is often performed without a selective breeding program, potentially leading to high levels of inbreeding. In this study, we identified and characterized autozygosity based on genome-wide runs of homozygosity (ROH) on a sample of parental and offspring individuals, determined effective population size (N_e), and assessed relatedness among parental individuals. A mean of 2,235 ± 526 and 1,841 ± 363 ROH segments per individual, resulting in a mean inbreeding coefficient of 0.33 ± 0.06 and 0.25 ± 0.06 were estimated for the progeny and parents, respectively. N_e was about 12 until four generations ago and at most 106 for 63 generations in the past, with varying genetic relatedness amongst the parents. This study shows the importance of genomic information when family relationships are unknown and the need of selective breeding programs for reproductive management decisions in the aquaculture industry.

Runs of homozygosity analysis reveals population characteristics of yellow-feathered chickens using re-sequencing data

1. To evaluate the inbreeding of yellow-feathered chickens (YFCs) and identify genes related to their unique characteristics, whole-genome re-sequencing data were applied to detect runs of homozygosity (ROH) in the genomes of ten YFC breeds. The number, length, distribution of ROH, and inbreeding coefficient in different YFC populations were calculated. Genomic regions with high frequency in ROH were annotated.2. In total, 25,547 ROH with an average length of 335 kb were detected, with most being <1 Mb. The domination of short ROH reflected the long breeding history of this chicken. The number, length, frequency, and distribution of ROH varied among chicken populations, and high genetic diversity was maintained.3. Numerous genes related to YFC characteristics were identified in the high-frequency ROH regions. Among these, IFNA, IFNB, IL11RA, IL22RA1, IFNLR1, and TRIF genes were involved in disease resistance. The AMY, G6PC, SDHB, GCNT4, and ACO genes were associated with energy material metabolism; and FABPL, AQP7, ACAA2, and RYR2 were related to meat quality and flavour. The KITLG, CREB3, RYR2, and LGR4 genes, related to pigmentation, were detected.4. This ROH-based inbreeding evaluation laid the foundation for breeding and conservation of YFC populations, and the candidate genes identified can be used for marker-assisted selection.

Unveiling Comparative Genomic Trajectories of Selection and Key Candidate Genes in Egg-Type Russian White and Meat-Type White Cornish Chickens

Comparison of genomic footprints in chicken breeds with different selection history is a powerful tool in elucidating genomic regions that have been targeted by recent and more ancient selection. In the present work, we aimed at examining and comparing the trajectories of artificial selection in the genomes of the native egg-type Russian White (RW) and meat-type White Cornish (WC) breeds. Combining three different statistics (top 0.1% SNP by F_ST value at pairwise breed comparison, hapFLK analysis, and identification of ROH island shared by more than 50% of individuals), we detected 45 genomic regions under putative selection including 11 selective sweep regions, which were detected by at least two different methods. Four of such regions were breed-specific for each of RW breed (on GGA1, GGA5, GGA8, and GGA9) and WC breed (on GGA1, GGA5, GGA8, and GGA28), while three remaining regions on GGA2 (two sweeps) and GGA3 were common for both breeds. Most of identified genomic regions overlapped with known QTLs and/or candidate genes including those for body temperatures, egg productivity, and feed intake in RW chickens and those for growth, meat and carcass traits, and feed efficiency in WC chickens. These findings were concordant with the breed origin and history of their artificial selection. We determined a set of 188 prioritized candidate genes retrieved from the 11 overlapped regions of putative selection and reviewed their functions relative to phenotypic traits of interest in the two breeds. One of the RW-specific sweep regions harbored the known domestication gene, TSHR. Gene ontology and functional annotation analysis provided additional insight into a functional coherence of genes in the sweep regions. We also showed a greater candidate gene richness on microchromosomes relative to macrochromosomes in these genomic areas. Our results on the selection history of RW and WC chickens and their key candidate genes under selection serve as a profound information for further conservation of their genomic diversity and efficient breeding.