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

Genetic diversity, demographic history, and selective signatures of Silkie chicken

BACKGROUND

Silkie is a traditional Chinese chicken breed characterized by its unique combination of specialized morphological traits. While previous studies have focused on the genetic basis of these traits, the overall genomic characteristics of the Silkie breed remain largely unexplored. In this study, we employed whole genome resequencing data to examine the genetic diversity, selective signals and demographic history of the Silkie breed through comparative analyses with seven other Chinese indigenous breeds (IDGBs), a commercial breed, and the wild ancestor Red Jungle Fowl.

RESULTS

In total, 20.8 million high-quality single nucleotide polymorphisms and 86 large structural variations were obtained. We discovered that Silkie exhibits a relatively high level of inbreeding and is genetically distinct from other IDGBs. Furthermore, our analysis indicated that Silkie has experienced a stronger historical population bottleneck and has a smaller effective population size compared with other IDGBs. We identified 45 putatively selected genes that are enriched in the melanogenesis pathway, which probably is related to the feather color. Among these genes, LMBR1 and PDSS2 have been previously associated with the extra toe and the hookless feathers, respectively. Six of the selected genes (KITLG, GSK3B, SOBP, CTBP1, ELMO2, SNRPN) are known to be associated with neurodevelopment and mental diseases in human, and are possibly related to the distinct behavior of Silkie. We further identified structural variants in Silkie and found previously reported variants linked to hyperpigmentation (END3), muff and beard (HOXB8), and Rose-comb phenotype (MNR2). Additionally, we found a 0.61 Mb inversion overlapping with the GMDS gene, which was previously linked to neurodevelopmental defects in zebrafish and humans. This may also be related to the behavior distinctiveness of Silkie.

CONCLUSIONS

Our study revealed that Silkie is genetically distinct and relatively highly inbred compared to other IDGB chicken populations, possibly attributed to more prolong population bottlenecks and selective breeding practice. These results enhance our understanding of how domestication and selective breeding have shaped the genome of Silkie. These findings contribute to the broader field of domestication and avian genomics, and have implications for the future conservation and breeding efforts.

Exploring Evolutionary Adaptations and Genomic Advancements to Improve Heat Tolerance in Chickens

Climate change poses a significant threat to the poultry industry, especially in hot climates that adversely affect chicken growth, development, and productivity through heat stress. This literature review evaluates the evolutionary background of chickens with the specific genetic characteristics that can help chickens to cope with hot conditions. Both natural selection and human interventions have influenced the genetic characteristics of the breeds used in the current poultry production system. The domestication of chickens from the Red junglefowl (Gallus gallus) has resulted in the development of various breeds with distinct genetic differences. Over the past few years, deliberate breeding for desirable traits (such as meat production and egg quality) in chickens has resulted in the emergence of various economically valuable breeds. However, this selective breeding has also caused a decrease in the genetic diversity of chickens, making them more susceptible to environmental stressors like heat stress. Consequently, the chicken breeds currently in use may possess a limited ability to adapt to challenging conditions, such as extreme heat. This review focuses on evaluating potential genes and pathways responsible for heat tolerance, including heat shock response, antioxidant defense systems, immune function, and cellular homeostasis. This article will also discuss the physiological and behavioral responses of chicken varieties that exhibit genetic resistance to heat, such as the naked neck and dwarf traits in different indigenous chickens. This article intends to review the current genomic findings related to heat tolerance in chickens that used methods such as the genome-wide association study (GWAS) and quantitative trait loci (QTL) mapping, offering valuable insights for the sustainability of poultry in the face of global warming.

Breeding practice of indigenous village chickens, and traits and breed preferences of smallholder farmers

BACKGROUND

Indigenous village chickens, or (IVCs), significantly contribute to rural livelihoods. Although natural selection has a disproportionate impact on the genetic structure of IVCs, farmers have developed locally tailored breeding practices to get the most out of their flocks.

OBJECTIVES

Small-scale farmers' insights on trait breed preferences and family flock breeding practices are presented in this cross-sectional study.

METHODS

A cross-sectional study was conducted in two agroecological zones using face-to-face individual interviews with 119 general informants.

RESULTS

Farmers prefer IVCs because they are multipurpose birds. Unlike policymakers, who usually underestimate the importance of IVCs, small-scale farmers acknowledge the coexistence of local and commercial chickens. Only 15.7% of farmers recruited homegrown cocks, whereas 47.9% outsourced breeding cocks from local markets and 36.4% from neighbours (χ-squared = 15.976, df = 2, p = 0.0003395). About 49.2% of small-scale farmers believed that consanguineous mating-induced inbreeding has only trivial effects. High flock turnover significantly reduces inbreeding. Regardless of the low production capacity, small-scale farmers prefer local (rank = 1.47) chickens to commercial (rank = 1.61). For cocks, fertility and growth traits were highly sought after, whereas for hens, maternal instincts and laying performance were prioritized. Compared to the highlands, the lowlands had a longer egg storage period (t = 2.677, df = 117, p = 0.009, 95% CI: -3.7607, -0.5622).

CONCLUSIONS

This study documented the wisdom of small-scale farmers and encouraged the incorporation of their insights into a sustainable genetic improvement program.

Genome-wide detections for runs of homozygosity and selective signatures reveal novel candidate genes under domestication in chickens

BACKGROUND

Indigenous chickens were developed through a combination of natural and artificial selection; essentially, changes in genomes led to the formation of these modern breeds via admixture events. However, their confusing genetic backgrounds include a genomic footprint regulating complex traits, which is not conducive to modern animal breeding.

RESULTS

To better evaluate the candidate regions under domestication in indigenous chickens, we considered both runs of homozygosity (ROHs) and selective signatures in 13 indigenous chickens. The genomes of Silkie feather chickens presented the highest heterozygosity, whereas the highest inbreeding status and ROH number were found in Luhua chickens. Short ROH (< 1 Mb), were the principal type in all chickens. A total of 291 ROH islands were detected, and QTLdb mapping results indicated that body weight and carcass traits were the most important traits. An ROH on chromosome 2 covering VSTM2A gene was detected in 12 populations. Combined analysis with the Tajima's D index revealed that 18 genes (e.g., VSTM2A, BBOX1, and RYR2) were under selection and covered by ROH islands. Transcriptional analysis results showed that RYR2 and BBOX1 were specifically expressed in the heart and muscle tissue, respectively.

CONCLUSION

Based on genome-wide scanning for ROH and selective signatures, we evaluated the genomic characteristics and detected significant candidate genes covered by ROH islands and selective signatures. The findings in this study facilitated the understanding of genetic diversity and provided valuable insights for chicken breeding and conservation strategies.

Detection of Runs of Homozygosity and Identification of Candidate Genes in the Whole Genome of Tunchang Pigs

Tunchang pigs are an indigenous pig population in China known for their high tolerance to roughage, delicious meat, and fecundity. However, the number of Tunchang pigs has been declining due to the influence of commercial breeds and African swine fever, which could potentially lead to inbreeding. To assess the inbreeding level and the genetic basis of important traits in Tunchang pigs, our research investigated the patterns in "runs of homozygosity" (ROHs) using whole genome resequencing data from 32 Tunchang pigs. The study aimed to determine the length, number, coverage, and distribution model of ROHs in Tunchang pigs, as well as genomic regions with high ROH frequencies. The results of the study revealed that a total of 20,499,374 single-nucleotide polymorphisms (SNPs) and 1953 ROH fragments were recognized in 32 individuals. The ROH fragments in Tunchang pigs were predominantly short, ranging from 0.5 to 1 megabases (Mb) in length. Furthermore, the coverage of ROHs varied across chromosomes, with chromosome 3 having the highest coverage and chromosome 11 having the lowest coverage. The genetic diversity of Tunchang pigs was found to be relatively high based on the values of HE (expected heterozygosity), HO (observed heterozygosity), pi (nucleotide diversity), Ne (effective population size), and MAF (minor allele frequency). The average inbreeding coefficients of Tunchang pigs, as determined by three different methods (FHOM, FGRM, and FROH), were 0.019, 0.0138, and 0.0304, respectively. These values indicate that the level of inbreeding in Tunchang pigs is currently low. Additionally, the study identified a total of 13 ROH islands on all chromosomes, which in total contained 38,913 SNPs and 120 genes. These ROH islands included genes associated with economically important traits, including meat quality (GYS1, PHLPP1, SLC27A5, and CRTC1), growth and development (ANKS1A, TAF11, SPDEF, LHB, and PACSIN1), and environmental adaptation (SLC26A7). The findings of this research offer valuable perspectives on the present status of Tunchang pig resources and offer a reference for breeding conservation plans and the efficient utilization of Tunchang pigs in the future. By understanding the inbreeding level and genetic basis of important traits in Tunchang pigs, conservation efforts can be targeted towards maintaining genetic diversity and promoting the sustainable development of this indigenous pig population.

Genetic admixture and diversity in Thai domestic chickens revealed through analysis of Lao Pa Koi fighting cocks

Lao Pa Koi (LPK) chicken is a popular fighting breed in Thailand, prized for (its unique characteristics acquired by selective breeding), and a valuable model for exploring the genetic diversity and admixture of red junglefowls and domestic chickens. In this study, genetic structure and diversity of LPK chicken were assessed using 28 microsatellite markers and mitochondrial DNA (mtDNA) D-loop sequences, and the findings were compared to a gene pool library from "The Siam Chicken Bioresource Project". High genetic variability was observed in LPK chickens using mtDNA D-loop haplotype analysis, and six haplotypes were identified. Microsatellite data revealed 182 alleles, with an average of 6.5 alleles per locus. These results confirmed the occurrence of genetic admixture of red junglefowl and Thai domestic chickens in LPK chicken breed. A maximum entropy modeling approach was used to analyze the spatial suitability and to assess the adaptive evolution of LPK chickens in diverse local environments. The model identified 82.52% of the area studied as unsuitable, and 9.34%, 7.11%, and 2.02% of the area indicated moderate, low, and high suitability, respectively. The highest contribution rate to land suitability for LPK chickens was found at an elevation of 100-250 m, suggesting the importance of elevation for their potential distribution. The results of this study provide valuable insights into the genetic origin of LPK chicken breed and identify resources for future genetic improvement.

Identification of Runs of Homozygosity Islands and Functional Variants in Wenchang Chicken

Wenchang chickens, a native breed in the Hainan province of China, are famous for their meat quality and adaptability to tropical conditions. For effective management and conservation, in the present study, we systematically investigated the characteristics of genetic variations and runs of homozygosity (ROH) along the genome using re-sequenced whole-genome sequencing data from 235 Wenchang chickens. A total of 16,511,769 single nucleotide polymorphisms (SNPs) and 53,506 ROH segments were identified in all individuals, and the ROH of Wenchang chicken were mainly composed of short segments (0-1 megabases (Mb)). On average, 5.664% of the genome was located in ROH segments across the Wenchang chicken samples. According to several parameters, the genetic diversity of the Wenchang chicken was relatively high. The average inbreeding coefficient of Wenchang chickens based on F_HOM, F_GRM, and F_ROH was 0.060 ± 0.014, 0.561 ± 0.020, and 0.0566 ± 0.01, respectively. A total of 19 ROH islands containing 393 genes were detected on 9 different autosomes. Some of these genes were putatively associated with growth performance (AMY1a), stress resistance (THEMIS2, PIK3C2B), meat traits (MBTPS1, DLK1, and EPS8L2), and fat deposition (LANCL2, PPARγ). These findings provide a better understanding of the degree of inbreeding in Wenchang chickens and the hereditary basis of the characteristics shaped under selection. These results are valuable for the future breeding, conservation, and utilization of Wenchang and other chicken breeds.

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.

Population Structure and Genetic Diversity Analysis of “Yufen 1” H Line Chickens Using Whole-Genome Resequencing

The effective protection and utilization of poultry resources depend on an accurate understanding of the genetic diversity and population structure. The breeding of the specialized poultry lineage “Yufen 1”, with its defined characteristics, was approved by the China Poultry Genetic Resource Committee in 2015. Thus, to investigate the relationship between the progenitor H line and other poultry breeds, the genetic diversity and population structure of “Yufen 1” H line (YF) were investigated and compared with those of 2 commercial chicken breeds, the ancestor breed Red Jungle Fowls, and 11 Chinese indigenous chicken breeds based on a whole-genome resequencing approach using 8,112,424 SNPs. The genetic diversity of YF was low, and the rate of linkage disequilibrium decay was significantly slower than that of the other Chinese indigenous breeds. In addition, it was shown that the YF population was strongly selected during intensive breeding and that genetic resources have been seriously threatened, which highlights the need to establish a systematic conservation strategy as well as utilization techniques to maintain genetic diversity within YF. Moreover, a principal component analysis, a neighbor-joining tree analysis, a structure analysis, and genetic differentiation indices indicated that YF harbors a distinctive genetic resource with a unique genetic structure separate from that of Chinese indigenous breeds at the genome level. The findings provide a valuable resource and the theoretical basis for the further conservation and utilization of YF.

Genome-Wide Genetic Structure of Henan Indigenous Chicken Breeds

There are five indigenous chicken breeds in Henan Province, China. These breeds have their own unique phenotypic characteristics in terms of morphology, behavior, skin and feather color, and productive performance, but their genetic basis is not well understood. Therefore, we analyzed the genetic structure, genomic diversity, and migration history of Henan indigenous chicken populations and the selection signals and genes responsible for Henan gamecock unique phenotypes using whole genome resequencing. The results indicate that Henan native chickens clustered most closely with the chicken populations in neighboring provinces. Compared to other breeds, Henan gamecock's inbreeding and selection intensity were more stringent. TreeMix analysis revealed the gene flow from southern chicken breeds into the Zhengyang sanhuang chicken and from the Xichuan black-bone chicken into the Gushi chicken. Selective sweep analysis identified several genes and biological processes/pathways that were related to body size, head control, muscle development, reproduction, and aggression control. Additionally, we confirmed the association between genotypes of SNPs in the strong selective gene LCORL and body size and muscle development in the Gushi-Anka F2 resource population. These findings made it easier to understand the traits of the germplasm and the potential for using the Henan indigenous chicken.

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.

Growth performance, survivability and profitability of improved smallholder chicken genetics in Nigeria: A COVID-19 intervention study

The impact of COVID-19 pandemic on smallholder farming households (SFH) includes increased poverty, and loss of livelihoods. Provision of livestock to SFH is a helpful intervention to mitigate this impact. This study provided a total of 150 smallholder poultry farmers, randomly selected from three states (Kebbi, Nasarawa, and Imo) in Nigeria, with ten 5-week-old chickens (mixed sexes) each, of either FUNAAB Alpha or Noiler chicken genetics. The improved, dual-purpose chickens were evaluated for growth performance (GP), survivability and profitability. The birds were managed under semi-scavenging production system. Body weight, mortality, and cost of production (COP) were recorded every 4 weeks until 21 weeks of age. Profitability was a function of the COP, and the selling price for live-birds (cocks). Body weight of Noiler (1,927 g) birds was not significantly (p > .05) higher than FUNAAB Alpha (1,792 g) at 21 weeks. Agroecology and genetics had significant (p < .05) effects on GP and survivability. Survivability of FUNAAB Alpha was higher (p < .05) than Noiler, with Nasarawa (81%-96%), having the highest (p < .0001) survival rate compared to Imo (62%-81%), and Kebbi (58%-75%). At 21 weeks, the number of cocks and hens differed significantly (p < .05) within the states (Imo: 2.4 ± .2 and 5.4 ± .3; Kebbi: 2.6 ± .2 and 5.5 ± .3; and Nasarawa: 2.9 ± .2 and 5.8 ± .3). Nasarawa (NGN 7,808; USD 19) ranked best for profitability, followed by Kebbi (NGN 6,545; USD 16) and Imo (NGN 5,875; USD 14). Overall, this study demonstrates that provision of improved chickens to vulnerable SFH in Nigeria holds great potential for economic growth, and resilience during emergencies, such as the COVID-19 pandemic.

Estimation of genetic variability and identification of regions under selection based on runs of homozygosity in Beijing-You Chickens

The genetic composition of populations is the result of a long-term process of selection and adaptation to specific environments and ecosystems. Runs of homozygosity (ROHs) are homozygous segments of the genome where the 2 haplotypes inherited from the parents are identical. The detection of ROH can be used to describe the genetic variability and quantify the level of inbreeding in an individual. Here, we investigated the occurrence and distribution of ROHs in 40 Beijing-You Chickens from the random breeding population (BJY_C) and 40 Beijing-You Chickens from the intramuscular fat (IMF) selection population (BJY_S). Principal component analysis (PCA) and maximum likelihood (ML) analyses showed that BJY_C was completely separated from the BJY_S. The nucleotide diversity of BJY_C was higher than that of BJY_S, and the decay rate of LD of BJY_C was faster. The ROHs were identified for a total of 7,101 in BJY_C and 9,273 in BJY_S, respectively. The ROH-based inbreeding estimate (F_ROH) of BJY_C was 0.079, which was significantly lower than that of BJY_S (F_ROH = 0.114). The results were the same as the estimates of the inbreeding coefficients calculated based on homozygosity (F_HOM), the correlation between uniting gametes (F_UNI), and the genomic relationship matrix (F_GRM). Additionally, the distribution and number of ROH islands in chromosomes of BJY_C and BJY_S were significantly different. The ROH islands of BJY_S that included genes associated with lipid metabolism and fat deposition, such as CIDEA and S1PR1, were absent in BJY_C. However, GPR161 was detected in both populations, which is a candidate gene for the formation of the unique five-finger trait in Beijing-You chickens. Our findings contributed to the understanding of the genetic diversity of random or artificially selected populations, and allowed the accurate monitoring of population inbreeding using genomic information, as well as the detection of genomic regions that affect traits under selection.

Comprehensive analysis of structural variants in chickens using PacBio sequencing

Structural variants (SVs) are one of the main sources of genetic variants and have a greater impact on phenotype evolution, disease susceptibility, and environmental adaptations than single nucleotide polymorphisms (SNPs). However, SVs remain challenging to accurately type, with several detection methods showing different limitations. Here, we explored SVs from 10 different chickens using PacBio technology and detected 49,501 high-confidence SVs. The results showed that the PacBio long-read detected more SVs than Illumina short-read technology genomes owing to some SV sites on chromosomes, which are related to chicken growth and development. During chicken domestication, some SVs beneficial to the breed or without any effect on the genomic function of the breed were retained, whereas deleterious SVs were generally eliminated. This study could facilitate the analysis of the genetic characteristics of different chickens and provide a better understanding of their phenotypic characteristics at the SV level, based on the long-read sequencing method. This study enriches our knowledge of SVs in chickens and improves our understanding of chicken genomic diversity.

Genomic insights into the conservation and population genetics of two Chinese native goat breeds

Chinese goats are an important group of goats worldwide. However, there are few studies on the conservation priority, genetic relationship, and potential gene flow between Chinese and global goat breeds. Here, we genotyped 239 goats from conservation populations of the Chinese Guangfeng and Ganxi breeds using the GoatSNP50 BeadChip. The conservation priority, population structure, selection signatures and introgression of these goats were analyzed in the context of 36 global goat breeds. First, we showed that Guangfeng and Ganxi goats had the largest effective population sizes across the global breeds 13 generations ago. Nevertheless, Ganxi goats have recently experienced a high degree of inbreeding, resulting in their conservation priority based on total gene and allelic diversities being lower than that of most other Chinese breeds (including Guangfeng goats). Population structure and admixture analyses showed that an average of 18% of Guangfeng genomic components were introgressed from Boer goats approximately 18-yr ago. Next, we reconstructed the subfamily structure of the core populations of Guangfeng and Ganxi goats, and proposed reasonable conservation strategies for inbreeding management. Moreover, a list of candidate genes under selection for fertility, immunity, growth, and meat quality were detected in Guangfeng and Ganxi goats. Finally, we identified some genes related to body development and reproduction, which were introgressed from Boer goats and may be beneficial for improving performance and productivity of Guangfeng goats. In conclusion, this study not only provides new insights into the conservation and utilization of Guangfeng and Ganxi goats but also enriches our understanding of artificial introgression from exotic goats into Chinese local goats.

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.

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.

Genetic Differentiation among Livestock Breeds-Values for Fst

(1) Background: The Fst statistic is widely used to characterize between-breed relationships. Fst = 0.1 has frequently been taken as indicating genetic distinctiveness between breeds. This study investigates whether this is justified. (2) Methods: A database was created of 35,080 breed pairs and their corresponding Fst values, deduced from microsatellite and SNP studies covering cattle, sheep, goats, pigs, horses, and chickens. Overall, 6560 (19%) of breed pairs were between breeds located in the same country, 7395 (21%) between breeds of different countries within the same region, 20,563 (59%) between breeds located far apart, and 562 (1%) between a breed and the supposed wild ancestor of the species. (3) Results: General values for between-breed Fst were as follows, cattle: microsatellite 0.06-0.12, SNP 0.08-0.15; sheep: microsatellite 0.06-0.10, SNP 0.06-0.17; horses: microsatellite 0.04-0.11, SNP 0.08-0.12; goats: microsatellite 0.04-0.14, SNP 0.08-0.16; pigs: microsatellite 0.06-0.27, SNP 0.15-0.22; chickens: microsatellite 0.05-0.28, SNP 0.08-0.26. (4) Conclusions: (1) Large amounts of Fst data are available for a substantial proportion of the world's livestock breeds, (2) the value for between-breed Fst of 0.1 is not appropriate owing to its considerable variability, and (3) accumulated Fst data may have value for interdisciplinary research.

Sequencing Reveals Population Structure and Selection Signatures for Reproductive Traits in Yunnan Semi-Fine Wool Sheep (Ovis aries)

Yunnan semi-fine wool sheep are among the most important cultivated sheep breeds in China. However, their population structure, genetic characteristics and traits of interest are poorly studied. In this study, we systematically studied the population characteristics and selection signatures of 40 Yunnan semi-fine wool sheep using SNPs obtained from whole-genome resequencing data. A total of 1393 Gb of clean data were acquired. The mapping rate against the reference genome was 91.23% on average (86.01%–92.26%), and the average sequence depth was 9.51X. After filtering, 28,593,198 SNPs and 4,725,259 indels with high quality were obtained. The heterozygosity rate, inbreeding coefficient and effective population size of the sheep were calculated to preliminarily explore their genetic characteristics. The average heterozygosity rate was 0.264, the average inbreeding coefficient was 0.0099, and the effective population size estimated from the heterozygote excess (HE) was 242.9. Based on the Tajima’s D and integrated haplotype score (iHS) approaches, 562 windows and 11,356 core SNPs showed selection signatures in the Yunnan semi-fine wool sheep population. After genome annotation and gene enrichment analysis, we found traces of early domestication in sensory organs, behavioural activity and the nervous system as well as adaptive changes in reproductive and wool traits under selection in this population. Some selected genes related to litter size, including FSHR, BMPR1B and OXT, were identified as being under selection. Specific missense mutations of the FSHR gene that differed from the reference genome were also identified in the population, and we found some SNP variations that may affect litter size. Our findings provide a theoretical basis for the conservation and utilization of Yunnan semi-fine wool sheep. Furthermore, our results reveal some changes common to sheep after domestication and provide a new opportunity to investigate the genetic variation influencing fecundity within a population evolving under artificial selection.

Genetic diversity and signatures of selection for heat tolerance and immune response in Iranian native chickens

Background

Understanding how evolutionary forces relating to climate have shaped the patterns of genetic variation within and between species is a fundamental pursuit in biology. Iranian indigenous chickens have evolved genetic adaptations to their local environmental conditions, such as hot and arid regions. In the present study, we provide a population genome landscape of genetic variations in 72 chickens representing nine Iranian indigenous ecotypes (Creeper, Isfahan, Lari, Marand, Mashhad, Naked neck, Sari, Shiraz and Yazd) and two commercial lines (White Leghorn and Arian). We further performed comparative population genomics to evaluate the genetic basis underlying variation in the adaptation to hot climate and immune response in indigenous chicken ecotypes. To detect genomic signatures of adaptation, we applied nucleotide diversity (θπ) and F_ST statistical measurements, and further analyzed the results to find genomic regions under selection for hot adaptation and immune response-related traits.

Results

By generating whole-genome data, we assessed the relationship between the genetic diversity of indigenous chicken ecotypes and their genetic distances to two different commercial lines. The results of genetic structure analysis revealed clustering of indigenous chickens in agreement with their geographic origin. Among all studied chicken groups, the highest level of linkage disequilibrium (LD) (~ 0.70) was observed in White Leghorn group at marker pairs distance of 1 Kb. The results from admixture analysis demonstrated evidence of shared ancestry between Arian individuals and indigenous chickens, especially those from the north of the country. Our search for potential genomic regions under selection in indigenous chicken ecotypes revealed several immune response and heat shock protein-related genes, such as HSP70, HSPA9, HSPH1, HSP90AB1 and PLCB4 that have been previously unknown to be involved in environmental-adaptive traits. In addition, we found some other candidate loci on different chromosomes probably related with hot adaptation and immune response-related traits.

Conclusions

The work provides crucial insights into the structural variation in the genome of Iranian indigenous chicken ecotypes, which up to now has not been genetically investigated. Several genes were identified as candidates for drought, heat tolerance, immune response and other phenotypic traits. These candidate genes may be helpful targets for understanding of the molecular basis of adaptation to hot environmental climate and as such they should be used in chicken breeding programs to select more efficient breeds for desert climate.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08434-7.

Runs of Homozygosity Uncover Potential Functional-Altering Mutation Associated With Body Weight and Length in Two Duroc Pig Lines

Runs of homozygosity (ROH) are widely used to investigate genetic diversity, demographic history, and positive selection signatures of livestock. Commercial breeds provide excellent materials to reveal the landscape of ROH shaped during the intense selection process. Here, we used the GeneSeek Porcine 50K single-nucleotide polymorphism (SNP) Chip data of 3,770 American Duroc (AD) and 2,096 Canadian Duroc (CD) pigs to analyze the genome-wide ROH. First, we showed that AD had a moderate genetic differentiation with CD pigs, and AD had more abundant genetic diversity and significantly lower level of inbreeding than CD pigs. In addition, sows had larger levels of homozygosity than boars in AD pigs. These differences may be caused by differences in the selective intensity. Next, ROH hotspots revealed that many candidate genes are putatively under selection for growth, sperm, and muscle development in two lines. Population-specific ROHs inferred that AD pigs may have a special selection for female reproduction, while CD pigs may have a special selection for immunity. Moreover, in the overlapping ROH hotspots of two Duroc populations, we observed a missense mutation (rs81216249) located in the growth and fat deposition-related supergene (ARSB-DMGDH-BHMT) region. The derived allele of this variant originated from European pigs and was nearly fixed in Duroc pigs. Further selective sweep and association analyses indicated that this supergene was subjected to strong selection and probably contributed to the improvement of body weight and length in Duroc pigs. These findings will enhance our understanding of ROH patterns in different Duroc lines and provide promising trait-related genes and a functional-altering marker that can be used for genetic improvement of pigs.

A bioinformatic search for correspondence between differentially expressed genes of domestic versus wild animals and orthologous human genes altering reproductive potential

One of the greatest achievements of genetics in the 20th century is D.K. Belyaev’s discovery of destabilizing selection during the domestication of animals and that this selection affects only gene expression regulation (not gene structure) and inf luences systems of neuroendocrine control of ontogenesis in a stressful environment. Among the experimental data generalized by Belyaev’s discovery, there are also f indings about accelerated extinc tion of testes’ hormonal function and disrupted seasonality of reproduction of domesticated foxes in comparison with their wild congeners. To date, Belyaev’s discovery has already been repeatedly conf irmed, for example, by independent observations during deer domestication, during the use of rats as laboratory animals, after the reintroduction of endangered species such as Przewalski’s horse, and during the creation of a Siberian reserve population of the Siberian grouse when it had reached an endangered status in natural habitats. A genome-wide comparison among humans, several domestic animals, and some of their wild congeners has given rise to the concept of self-domestication syndrome, which includes autism spectrum disorders. In our previous study, we created a bioinformatic model of human self-domestication syndrome using differentially expressed genes (DEGs; of domestic animals versus their wild congeners) orthologous to the human genes (mainly, nervous-system genes) whose changes in expression affect reproductive potential, i. e., growth of the number of humans in the absence of restrictions caused by limiting factors. Here, we applied this model to 68 human genes whose changes in expression alter the reproductive health of women and men and to 3080 DEGs of domestic versus wild animals. As a result, in domestic animals, we identif ied 16 and 4 DEGs, the expression changes of which are codirected with changes in the expression of the human orthologous genes decreasing and increasing human reproductive potential, respectively. The wild animals had 9 and 11 such DEGs, respectively. This difference between domestic and wild animals was signif icant according to Pearson’s χ2 test (p < 0.05) and Fisher’s exact test (p < 0.05). We discuss the results from the standpoint of restoration of endangered animal species whose natural habitats are subject to an anthropogenic impact.

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.

A genome-wide scan to identify signatures of selection in two Iranian indigenous chicken ecotypes

Background

Various regions of the chicken genome have been under natural and artificial selection for thousands of years. The substantial diversity that exits among chickens from different geographic regions provides an excellent opportunity to investigate the genomic regions under selection which, in turn, will increase our knowledge about the mechanisms that underlie chicken diversity and adaptation. Several statistics have been developed to detect genomic regions that are under selection. In this study, we applied approaches based on differences in allele or haplotype frequencies (F_ST and hapFLK, respectively) between populations, differences in long stretches of consecutive homozygous sequences (ROH), and differences in allele frequencies within populations (composite likelihood ratio (CLR)) to identify inter- and intra-populations traces of selection in two Iranian indigenous chicken ecotypes, the Lari fighting chicken and the Khazak or creeper (short-leg) chicken.

Results

Using whole-genome resequencing data of 32 individuals from the two chicken ecotypes, approximately 11.9 million single nucleotide polymorphisms (SNPs) were detected and used in genomic analyses after quality processing. Examination of the distribution of ROH in the two populations indicated short to long ROH, ranging from 0.3 to 5.4 Mb. We found 90 genes that were detected by at least two of the four applied methods. Gene annotation of the detected putative regions under selection revealed candidate genes associated with growth (DCN, MEOX2 and CACNB1), reproduction (ESR1 and CALCR), disease resistance (S1PR1, ALPK1 and MHC-B), behavior pattern (AGMO, GNAO1 and PSEN1), and morphological traits (IHH and NHEJ1).

Conclusions

Our findings show that these two phenotypically different indigenous chicken populations have been under selection for reproduction, immune, behavioral, and morphology traits. The results illustrate that selection can play an important role in shaping signatures of differentiation across the genomic landscape of two chicken populations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00664-9.

An Overview of the Use of Genotyping Techniques for Assessing Genetic Diversity in Local Farm Animal Breeds

Simple Summary

The number of local farm animal breeds is declining worldwide. However, these breeds have different degrees of genetic diversity. Measuring genetic diversity is important for the development of conservation strategies and, therefore, various genomic analysis techniques are available. The aim of the present work was to shed light on the use of these techniques in diversity studies of local breeds. In summary, a total of 133 worldwide studies that examined genetic diversity in local cattle, sheep, goat, chicken and pig breeds were reviewed. The results show that over time, almost all available genomic techniques were used and various diversity parameters were calculated. Therefore, the present results provide a comprehensive overview of the application of these techniques in the field of local breeds. This can provide helpful insights into the advancement of the conservation of breeds with high genetic diversity.

Abstract

Globally, many local farm animal breeds are threatened with extinction. However, these breeds contribute to the high amount of genetic diversity required to combat unforeseen future challenges of livestock production systems. To assess genetic diversity, various genotyping techniques have been developed. Based on the respective genomic information, different parameters, e.g., heterozygosity, allele frequencies and inbreeding coefficient, can be measured in order to reveal genetic diversity between and within breeds. The aim of the present work was to shed light on the use of genotyping techniques in the field of local farm animal breeds. Therefore, a total of 133 studies across the world that examined genetic diversity in local cattle, sheep, goat, chicken and pig breeds were reviewed. The results show that diversity of cattle was most often investigated with microsatellite use as the main technique. Furthermore, a large variety of diversity parameters that were calculated with different programs were identified. For 15% of the included studies, the used genotypes are publicly available, and, in 6%, phenotypes were recorded. In conclusion, the present results provide a comprehensive overview of the application of genotyping techniques in the field of local breeds. This can provide helpful insights to advance the conservation of breeds.

Genetic Divergence and Population Structure in Weedy and Cultivated Broomcorn Millets (Panicum miliaceum L.) Revealed by Specific-Locus Amplified Fragment Sequencing (SLAF-Seq)

Broomcorn millet (Panicum miliaceum L.) is one of the earliest domesticated crops in the world. Weedy broomcorn millet [Panicum ruderale (Kitag.) Chang or Panicum miliaceum subsp. ruderale (Kitag.) Tzvel] is thought to be the descendant of the wild ancestor or the feral type of this cereal. The genealogical relationships and genetic divergence among these taxa have not been clarified. In this study, the genetic diversity and population structure of weedy and cultivated broomcorn millets were investigated by using the high-throughput sequencing technology, i.e., the specific-locus amplified fragment sequencing (SLAF-seq). Our analyses consistently revealed both the wild and the feral genotypes in the weedy broomcorn millets. The single nucleotide polymorphisms (SNPs) at the genomic level provided useful evidence to distinguish the wild and the endoferal/exoferal types of weedy broomcorn millets. The genetic divergence revealed between the cultivated broomcorn millet from eastern Eurasia and those from central-western Eurasia was probably derived from either the genetic introgression from weedy broomcorn millets along the spread routes or the founder effect, while the limited gene flow of broomcorn millets from eastern and central-western Eurasia was probably due to the different uses of broomcorn millets and eating habits of the local people.

Genetic diversity among two native Indian chicken populations using cytochrome c oxidase subunit I and cytochrome b DNA barcodes

BACKGROUND AND AIM

India has large varieties (recognized, unrecognized) of native chickens (Desi) scattered throughout the country, managed under scavenging system different from commercial chicken breeds. However, they are less investigated for genetic diversity they harbor. The present study was planned to evaluate genetic diversity among two native chicken populations of North Gujarat (proposed Aravali breed) and South Gujarat (Ankleshwar breed). Aravali chicken, a distinct population with unique characters different from the registered chicken breeds of India is under process to be registered as a new chicken breed of Gujarat, India.

MATERIALS AND METHODS

Two mitochondrial markers, namely, cytochrome oxidase c subunit I (COX I) and cytochrome b (Cyt b) genes were studied across 10 birds from each population. Methodology included sample collection (blood), DNA isolation (manual), polymerase chain reaction amplification of mitochondrial genes, Sanger sequencing, and purification followed by data analysis using various softwares.

RESULTS

Haplotype analysis of the COX I gene unveiled a total eight and three haplotypes from the Aravali and Ankleshwar populations, respectively, with haplotype diversity (Hd) of 92.70 % for the Aravali and 34.50% for the Ankleshwar breed. Haplotype analysis of the Cyt b gene revealed a total of four haplotypes from the Aravali population with 60% Hd and no polymorphism in Ankleshwar breed. The phylogenetic analysis uncovered Red Jungle Fowl and Gray Jungle Fowl as prime roots for both populations and all domestic chicken breeds.

CONCLUSION

Study findings indicated high genetic variability in Aravali chicken populations with COX I mitochondrial marker being more informative for evaluating genetic diversity in chickens.

Population structure, inbreeding and admixture in local cattle populations managed by community-based breeding programs in Burkina Faso

High-throughput genomic markers provide an opportunity to assess important indicators of genetic diversity for populations managed in livestock breeding programs. While well-structured breeding programs are common in developed countries, in developing country situations, especially in West Africa, on-farm performance and pedigree recordings are rare, and thus, genomic markers provide insights to the levels of genetic diversity, inbreeding and introgression by other breeds. In this study, we analysed key population parameters such as population structure, admixture and levels of inbreeding in three neighbouring populations of African taurine and taurine × Zebu crosses managed by community-based breeding programs in the South-West of Burkina Faso. The three populations were pure Baoulé (called Lobi locally) in sedentary production systems, Baoulé x Zebu crossbreds in sedentary systems and Zebu × Baoulé crossbreds in transhumant production systems, respectively. The total sample analysed included 631 animals and 38,207 single nucleotide polymorphisms after quality control. Results of principal component and admixture analyses confirmed the genetic background of two distinct ancestral populations (taurine and zebuine) and levels of admixture in all three breeding populations, including the presumably pure Baoulé group of animals. Inbreeding levels were moderate, compared to European dairy and beef cattle populations and higher than those of Brazilian Nellore cattle. Very few animals with inbreeding levels indicating parent-offspring or full sib mating were observed, and inbreeding levels indicating half sib mating were also rare. For the management of breeding populations, farmers were advised to exchange best young bulls. The crossbreeding levels of presumably pure Baoulé animals are of concern to the breeding program due to the high level of endangerment of pure African taurine cattle populations across West Africa. Future rounds of bull selection in the community-based breeding program will make use of genomic information about admixture levels.

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.

Runs of Homozygosity in Modern Chicken Revealed by Sequence Data

Runs of homozygosity (ROH) are chromosomal stretches that in a diploid genome appear in a homozygous state and display identical alleles at multiple contiguous loci. This study aimed to systematically compare the genomic distribution of the ROH islands among five populations of wild vs commercial chickens of both layer and broiler type. To this end, we analyzed whole genome sequences of 115 birds including white layer (WL, n = 25), brown layer (BL, n = 25), broiler line A (BRA, n = 20), broiler line B (BRB, n = 20) and Red Junglefowl (RJF, n = 25). The ROH segments varied in size markedly among populations, ranging from 0.3 to 21.83 Mb reflecting their past genealogy. White layers contained the largest portion of the genome in homozygous state with an average ROH length of 432.1 Mb (±18.7) per bird, despite carrying it in short segments (0.3-1 Mb). Population-wise inbreeding measures based on Wright's (F_is) and genomic (F_ROH) metrics revealed highly inbred genome of layer lines relative to the broilers and Red Junglefowl. We further revealed the ROH islands, among commercial lines overlapped with QTL related to limb development (GREM1, MEOX2), body weight (Meis2a.1, uc_338), eggshell color (GLCCI1, ICA1, UMAD1), antibody response to Newcastle virus (ROBO2), and feather pecking. Comparison of ROH landscape in sequencing resolution demonstrated that a sizable portion of genome of commercial lines segregates in homozygote state, reflecting many generations of assortative mating and intensive selection in their recent history. In contrary, wild birds carry shorter ROH segments, likely suggestive of older evolutionary events.