Whole-genome resequencing of Xishuangbanna fighting chicken to identify signatures of selection

Genetic characteristics and selection signatures between Southern Chinese local and commercial chickens

The introduction of exotic breeds and the cultivation of new lines by breeding companies have posed challenges to native chickens in South China, including loss of breed characteristics, decreased genetic diversity, and declining purity. Understanding the population genetic structure and genetic diversity of native chickens in South China is crucial for further advancements in breeding efforts. In this study, we analyzed the population genetic structure and genetic diversity of 321 individuals from 10 different breeds in South China. By comparing commercial chickens with native ones, we identified selection signatures occurring between local chickens and commercial breeds. The analysis of population genetic structure revealed that the native chicken populations in South China exhibited a considerable level of genetic diversity. Moreover, the commercial lines of Xiaobai chicken and Huangma chicken displayed even higher levels of genetic diversity, which distinguished them from other native varieties at the clustering level. However, certain individuals within these commercial varieties showed a discernible genetic relationship with the native populations. Notably, both commercial varieties also retained a significant degree of genetic similarity to their respective native counterparts. In order to investigate the genomic changes occurring during the commercialization of native chickens, we employed 4 methods (Fst, ROD, XPCLR, and XPEHH) to identify potential candidate regions displaying selective signatures in Southern Chinese native chicken population. A total of 168 (identified by Fst and ROD) and 86 (identified by XPCLR and XPEHH) overlapping genes were discovered. Functional annotation analysis revealed that these genes may be associated with reproduction and growth (SAMSN1, HYLS1, ROBO3, FGF14, PRSS23), musculoskeletal development (DNER, MYBPC1, DGKB, ORC1, KLF10), disease resistance and environmental adaptability (PUS3, CRB2, CALD1, USP15, SGCD, LTBP1), as well as egg production (ADGRB3, ACSF3). Overall, native chickens in South China harbor numerous selective sweep regions compared to commercial chickens, enriching valuable genomic resources for future genetic research and breeding conservation.

Integrating genomics and transcriptomics to identify candidate genes for high-altitude adaptation and egg production in Nixi chicken

1. This study combined genome-wide selection signal analysis with RNA-sequencing to identify candidate genes associated with high altitude adaptation and egg production performance in Nixi chickens (NXC).2. Based on the whole-genome data from 20 NXC (♂:10; ♀:10), the population selection signal was analysed by sliding window analysis. The selected genes were screened by combination with the population differentiation statistic (FST). The sequence diversity statistic (θπ). RNA-seq was performed on the ovarian tissues of NXC (n = 6) and Lohmann laying hens (n = 6) to analyse the differentially expressed genes (DEGs) between the two groups. The functional enrichment analysis of the selected genes and differentially expressed genes was performed.3. There were 742 genes under strong positive selection and 509 differentially expressed genes screened in NXC. Integrated analysis of the genome and transcriptome revealing 26 overlapping genes. The candidate genes for adaptation to a high-altitude environment, as well as for egg production, disease resistance, vision and pigmentation in NXC were preliminarily screened.4. The results provided theoretical guidance for further research on the genetic resource protection and utilisation of NXC.

Artificial selection footprints in indigenous and commercial chicken genomes

BACKGROUND

Although many studies have been done to reveal artificial selection signatures in commercial and indigenous chickens, a limited number of genes have been linked to specific traits. To identify more trait-related artificial selection signatures and genes, we re-sequenced a total of 85 individuals of five indigenous chicken breeds with distinct traits from Yunnan Province, China.

RESULTS

We found 30 million non-redundant single nucleotide variants and small indels (< 50 bp) in the indigenous chickens, of which 10 million were not seen in 60 broilers, 56 layers and 35 red jungle fowls (RJFs) that we compared with. The variants in each breed are enriched in non-coding regions, while those in coding regions are largely tolerant, suggesting that most variants might affect cis-regulatory sequences. Based on 27 million bi-allelic single nucleotide polymorphisms identified in the chickens, we found numerous selective sweeps and affected genes in each indigenous chicken breed and substantially larger numbers of selective sweeps and affected genes in the broilers and layers than previously reported using a rigorous statistical model. Consistent with the locations of the variants, the vast majority (~ 98.3%) of the identified selective sweeps overlap known quantitative trait loci (QTLs). Meanwhile, 74.2% known QTLs overlap our identified selective sweeps. We confirmed most of previously identified trait-related genes and identified many novel ones, some of which might be related to body size and high egg production traits. Using RT-qPCR, we validated differential expression of eight genes (GHR, GHRHR, IGF2BP1, OVALX, ELF2, MGARP, NOCT, SLC25A15) that might be related to body size and high egg production traits in relevant tissues of relevant breeds.

CONCLUSION

We identify 30 million single nucleotide variants and small indels in the five indigenous chicken breeds, 10 million of which are novel. We predict substantially more selective sweeps and affected genes than previously reported in both indigenous and commercial breeds. These variants and affected genes are good candidates for further experimental investigations of genotype-phenotype relationships and practical applications in chicken breeding programs.

Analysis of the Selection Signal of the Tibetan Black Chicken Genome Based on Whole-Genome Sequencing

BACKGROUND

The Tibetan chicken has adapted well to high altitudes genetically after its long-term habitation in the plateau. In this study, we analyzed the selection signal of Tibetan black chickens (TBCs) and discovered genes associated with the characteristics of germplasm.

METHODS

Whole-genome sequencing (WGS) was used to identify the single-nucleotide polymorphism (SNP) markers and genetic structures in the genome of Tibetan black chickens. Further, we performed a comparative population genomics analysis between the genomic data obtained in this present study and the genomic data for five wild red jungle fowls (RJFs) accessed from the NCBI database (GenBank accession number PRJNA241474). Thereafter, the Fst and Pi selections were used to identify genes under positive selection in the Tibetan black chicken genome.

RESULTS

A total of 9,490,690 SNPs were identified in the Tibetan black chickens. In addition, the results from the gene ontology (GO) analysis showed that 732 genes of TBCs were enriched in a total of 210 GO terms with specific molecular functions such as regulation of cellular catabolic process, the MAPK signaling pathway, regulation of ion transport, growth, morphogenesis and lung alveolus development which may provide a better mechanism to facilitate oxygen transport and utilization in TBCs. Moreover, the results from the KEGG analysis showed that 732 genes of the TBCs were significantly enriched in the calcium signaling pathway, circadian entrainment (ADCY1, GNG7 and PER3), oxytocin signaling pathway and pathways of multiple neurodegeneration diseases. In addition, the CD86 antigen (CD86) was identified as a gene associated with the immune response in chickens. It was also revealed that genes such as TRIT1, HPCAL4, NT5C1A and HEYL were discovered under selection in Tibetan black chickens on chromosome 23. These genes may be related to the local adaptive characteristics of Tibetan black chickens, for instance, NT5C1A and HEYL may be involved in the high-altitude adaption of oxygen delivery in Tibetan black chickens.

CONCLUSIONS

In summary, we found that selection mainly affects the disease resistance and cold acclimatization of Tibetan black chickens. Hence, these results may provide important genetic information for the evolution and breeding of Tibetan black chickens.

Systematic Selection Signature Analysis of Chinese Gamecocks Based on Genomic and Transcriptomic Data

Selection pressures driven by natural causes or human interference are key factors causing genome variants and signatures of selection in specific regions of the genome. Gamecocks were bred for cockfighting, presenting pea-combs, larger body sizes, stronger limbs, and higher levels of aggression than other chickens. In this study, we aimed to explore the genomic differences between Chinese gamecocks and commercial, indigenous, foreign, and cultivated breeds by detecting the regions or sites under natural or artificial selection using genome-wide association studies (GWAS), genome-wide selective sweeps based on the genetic differentiation index (F_ST), and transcriptome analyses. Ten genes were identified using GWAS and F_ST: gga-mir-6608-1, SOX5, DGKB, ISPD, IGF2BP1, AGMO, MEOX2, GIP, DLG5, and KCNMA1. The ten candidate genes were mainly associated with muscle and skeletal development, glucose metabolism, and the pea-comb phenotype. Enrichment analysis results showed that the differentially expressed genes between the Luxi (LX) gamecock and Rhode Island Red (RIR) chicken were mainly related to muscle development and neuroactive-related pathways. This study will help to understand the genetic basis and evolution of Chinese gamecocks and support the further use of gamecocks as an excellent breeding material from a genetic perspective.

The study of selection signature and its applications on identification of candidate genes using whole genome sequencing data in chicken - a review

Chicken is a major source of protein for the increasing human population and is useful for research purposes. There are almost 1,600 distinct regional breeds of chicken across the globe, among which a large body of genetic and phenotypic variations has been accumulated due to extensive natural and artificial selection. Moreover, natural selection is a crucial force for animal domestication. Several approaches have been adopted to detect selection signatures in different breeds of chicken using whole genome sequencing (WGS) data including integrated haplotype score (iHS), cross-populated extend haplotype homozygosity test (XP-EHH), fixation index (F_ST), cross-population composite likelihood ratio (XP-CLR), nucleotide diversity (Pi), and others. In addition, gene enrichment analyses are utilized to determine KEGG pathways and gene ontology (GO) terms related to traits of interest in chicken. Herein, we review different studies that have adopted diverse approaches to detect selection signatures in different breeds of chicken. This review systematically summarizes different findings on selection signatures and related candidate genes in chickens. Future studies could combine different selection signatures approaches to strengthen the quality of the results thereby providing more affirmative inference. This would further aid in deciphering the importance of selection in chicken conservation for the increasing human population.

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.

Genome-wide scan for selection signatures and genes related to heat tolerance in domestic chickens in the tropical and temperate regions in Asia

Heat stress is one of the major environmental stressors challenging the global poultry industry. Identifying the genes responsible for heat tolerance is fundamentally important for direct breeding programs. To uncover the genetic basis underlying the ambient temperature adaptation of chickens, we analyzed a total of 59 whole genomes from indigenous chickens that inhabit South Asian tropical regions and temperate regions from Northern China. We applied F_ST and π-ratio to scan selective sweeps and identified 34 genes with a signature of positive selection in chickens from tropical regions. Several of these genes are functionally implicated in metabolism (FABP2, RAMP3, SUGCT, and TSHR) and vascular smooth muscle contractility (CAMK2), and they may be associated with adaptation to tropical regions. In particular, we found a missense mutation in thyroid-stimulating hormone receptor (41020238:G>A) that shows significant differences in allele frequency between the chicken populations of the two regions. To evaluate whether the missense mutation in TSHR could enhance the heat tolerance of chickens, we constructed segregated chicken populations and conducted heat stress experiments using homozygous mutations (AA) and wild-type (GG) chickens. We found that GG chickens exhibited significantly higher concentrations of alanine aminotransferase, lactate dehydrogenase, and creatine kinase than AA chickens under heat stress (35 ± 1°C) conditions (P < 0.05). These results suggest that TSHR (41020238:G>A) can facilitate heat tolerance and adaptation to higher ambient temperature conditions in tropical climates. Overall, our results provide potential candidate genes for molecular breeding of heat-tolerant chickens.

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.

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.

Analysis of selection signatures on the Z chromosome of bidirectional selection broiler lines for the assessment of abdominal fat content

Background

The discovery of selection signatures has enabled the identification of genomics regions under selective pressure, enhancing knowledge of evolutionary genotype-phenotypes. Sex chromosomes play an important role in species formation and evolution. Therefore, the exploration of selection signatures on sex chromosomes has important biological significance.

Results

In this study, we used the Cross Population Extend Haplotype Homozygosity Test (XPEHH), F-statistics (F_ST) and EigenGWAS to assess selection signatures on the Z chromosome in 474 broiler chickens via Illumina chicken 60 K SNP chips. SNP genotype data were downloaded from publicly available resources. We identified 17 selection regions, amongst which 1, 11 and 12 were identified by XPEHH, F_ST, and EigenGWAS, respectively. Each end of the Z chromosome appeared to undergo the highest levels of selection pressure. A total of 215 candidate genes were located in 17 selection regions, some of which mediated lipogenesis, fatty acid production, fat metabolism, and fat decomposition, including FGF10, ELOVL7, and IL6ST. Using abdominal adipose tissue expression data of the chickens, 187 candidate genes were expressed with 15 differentially expressed genes (DEGs) in fat vs. lean lines identified. Amongst the DEGs, VCAN was related to fat metabolism. GO pathway enrichment analysis and QTL annotations were performed to fully characterize the selection mechanism(s) of chicken abdominal fat content.

Conclusions

We have found some selection regions and candidate genes involving in fat metabolism on the Z chromosome. These findings enhance our understanding of sex chromosome selection signatures.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-021-00971-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.

Revisiting the evolutionary history of domestic and wild ducks based on genomic analyses

Although domestic ducks have been important poultry species throughout human history, their origin remains enigmatic, with mallards and/or Chinese spot-billed ducks being proposed as the direct wild ancestor(s) of domestic ducks. Here, we analyzed 118 whole genomes from mallard, Chinese spot-billed, and domestic ducks to reconstruct their evolutionary history. We found pervasive introgression patterns among these duck populations. Furthermore, we showed that domestic ducks separated from mallard and Chinese spot-billed ducks nearly 38 thousand years ago (kya) and 54 kya, respectively, which is considerably outside the time period of presumed duck domestication. Thus, our results suggest that domestic ducks may have originated from another wild duck population that is currently undefined or unsampled, rather than from present-day mallard and/or Chinese spot-billed ducks, as previously thought. Overall, this study provides new insight into the complex evolution of ducks.

Identifying Signatures of Selection Related to Comb Development

The aim of this study was to identify genes involved in comb development to provide insights into the molecular mechanism of chickens' comb formation. Fixation index (FST) and average number of base differences (π) of males with large and small combs were calculated based on whole-genome resequencing data. Chromosome regions with larger FST values and smaller π were considered candidate selection regions. Through further annotation of gene functions and pathways, we sought to screen possible selected genes associated with comb development. By screening whole genome resequencing data, FST and π were calculated using a 40 Kb sliding window strategy and eight regions were identified. Quantitative trait loci (QTL; FOX1 gene) related to comb length were found on chromosome 1. QTL (GLP1R, BTBD9, MIR6633, and MDGA1 genes) related to comb weight were found on chromosome 3. QTL (ALDH1A1, TMC1, and ANXA1 genes) associated with comb area were found on the Z chromosome. Nineteen genes, Wnt signaling pathway and neuroactive ligand-receptor interaction signaling pathway directly or indirectly related to comb growth and development were found through functional annotation and GO analysis. Among the selected genes LYN, GLP1R, FOX1, TBK1, STRAP, ST6GALNAC, and Wnt signaling pathways were related to immunity. MDGA1, BTBD9, MTSS1, SrGAPs, and neuroactive ligand receptor interaction signaling pathways related to neural function were screened. ALDH1A1, ANXAl, THBS, HIF-1α, and ACTN1 genes were related to heat dissipation. Among the selected genes FOX1, MDGAl, and ANXAl associated with immunity, neurological function, and heat dissipation function coincided with genes affecting the length, weight, and area of the comb. Comprehensive analysis suggested that comb development was due to multiple genes and signaling pathways.

Chickens productivity selection affects immune system genes

The quantitative trait loci associated with the immune properties of chickens are of interest from the point of view of obtaining animals resistant to infectious agents using marker-assisted selection. In the process of selecting markers for genomic selection in broiler-type chickens, a non-standard genotype frequency of the RACK1 gene allele (SNP Gga_rs15788101) in the B5 line of broiler-type chicken cross Smena 8 was identified and it was suggested that this gene was involved in selection. Therefore, it was decided to investigate the available polymorphisms in the three genes responsible for the IgY titer (DMA, RACK1 and CD1B). Molecular typing of single nucleotide polymorphisms of three loci revealed an approach to fixation of the unfavorable allele of the DMA gene (SNP Gga_rs15788237), an approach to fixation of the unfavorable allele of the RACK1 gene and the prevalence of the favorable CD1B gene allele (SNP Gga_rs16057130). Analysis of the haplotypes revealed a strong linkage disequilibrium of these genes. This suggests that these genes experience selection pressure. Analysis of the protein-coding sequences of the CD1B and DMA genes of various breeds of chickens revealed a negative selection of these genes. In order to understand whether the fixation of the studied alleles is the result of artificial selection of the B5 line of the cross Smena 8, an analysis of similar loci in layer chickens Hisex White was carried out. The frequencies of the alleles at the loci of the CD1B gene (Gga_rs16057130) and the RACK1 gene (Gga_rs15788101) in the Hisex White chicken genome differ from the frequencies of the alleles obtained for chickens of the B5 line of the cross Smena 8. It can be assumed that the fixation of the allele in the DMA gene (SNP Gga_rs15723) is associated with artificial or natural selection, consistent in broilers and layers. Changes in the loci Gga_rs16057130 and Gga_rs15788101 in the B5 line of the Smena 8 chickens are most likely associated with artificial selection of broiler productivity traits, which can subsequently lead to fixation of alleles at these loci. Artificial breeding of chickens leads to degradation of the variability of genes encoding elements of the immune system, which can cause a decrease in resistance to various diseases. The study of the negative impact of selection of economic traits on immunity should provide means to mitigate negative consequences and help find ways to obtain disease-resistant animals.

Analyses of mitogenomic markers shed light on the divergence, population dynamics, and demographic history of Pakistani chickens

Pakistan is one of a few sites, associated with the earliest known independent domestication event in the evolutionary history of chicken, which is socio-economically and historically the most important poultry bird in the country. However, the divergence, past population dynamics, and demographic history of Pakistani chickens have not been addressed so far. Therefore, we herein investigated the indigenous Pakistani chickens using mitogenomic markers. We first prepared individual DNA samples from the chicken feathers, and generated nucleotide sequence data, which was then subjected to various population genetics analyses. In molecular phylogenetic analysis, the Pakistani chickens were clustered under nine different clades. Among the wild fowls, the Indian red jungle fowl (IRJF) shared very close affinities to Pakistani chickens. The Bayesian skyline plot showed an increase in the effective population size of Pakistani chickens during the last 50 years. Finally, a time-calibrated phylogeny inferred molecular divergence of the Pakistani chickens. A molecular rate of 3.6 × 10^-6 mutations/site/year (95% HPD interval: 2.28 × 10^-8 to 9.32 × 10^-6) was estimated for the data set. In a rooted tree with root-age of 12058 years (95% HPD interval: 1161-38411), the Pakistani chicken haplotypes showed divergence from IRJF haplotypes around 6987 years (95% HPD interval: 1132-20746) ago, and they shared their most recent common ancestor with Gallus gallus spadiceus, and G. g. jabouillei at the root of the tree. Overall, these results suggest that Pakistani chicken haplotypes share their ancestral gene pool with the IRJF as compared to other red jungle fowl subspecies.

Genomic variations and signatures of selection in Wuhua yellow chicken

Wuhua yellow chicken (WHYC) is an important traditional yellow-feathered chicken from China, which is characterized by its white tail feathers, white flight feathers, and strong disease resistance. However, the genomic basis of these unique traits associated with WHYC is poorly understood. In this study, whole-genome resequencing was performed with an average coverage of 20.77-fold to investigate heritable variation and identify selection signals in WHYC. Reads were mapped onto the chicken reference genome (Galgal5) with a coverage of 85.95%. After quality control, 11,953,471 single nucleotide polymorphisms and 1,069,574 insertion/deletions were obtained. In addition, 41,408 structural variants and 33,278 copy number variants were found. Comparative genomic analysis of WHYC and other yellow-feathered chicken breeds showed that selected regions were enriched in genes involved in transport and catabolism, immune system, infectious diseases, signal transduction, and signaling molecules and interactions. Several genes associated with disease resistance were also identified, including IFNA, IFNB, CD86, IL18, IL11RA, VEGFC, and ATG10. Furthermore, our results suggest that PMEL and TYRP1 may contribute to the white feather coloring in WHYC. These findings can improve our understanding of the genetic characteristics of WHYC and may contribute to future breed improvement.

Whole-genome resequencing of Dulong Chicken reveal signatures of selection

1. Dulong Chickens (DLCs) live at high altitude (~3000 m) and humidity (~90%), are endemic to the Yunnan province, and have gradually developed unique physiological characteristics, but their genetic basis is still unclear. Using the fixation index (F_ST ) approach, based on whole-genome resequencing, DLCs were analysed to uncover the genomic architecture of the population and candidate genes involved in selection during domestication. 2. A total of 469 candidate genes were obtained to be putatively under selection in DLCs. Further investigations revealed the genic footprint for local adaptation (high-altitude and high-humidity) as the genic signatures that are involved in economic traits (related to egg production). 3. Candidate genes were identified that may be associated with disease resistance, aggressiveness, small body size and positive selection of vision in DLCs. 4. These data revealed loci of selective signals that operate during selection for production at high altitude and humidity.

Genome diversity of Chinese indigenous chicken and the selective signatures in Chinese gamecock chicken

Gamecock chickens are one of the earliest recorded birds in China, and have accumulated some unique morphological and behavioral signatures such as large body size, muscularity and aggressive behavior, whereby being excellent breeding materials and a good model for studying bird muscular development and behavior. In this study, we sequenced 126 chicken genomes from 19 populations, including four commercial chicken breeds that are commonly farmed in China, 13 nationwide Chinese typical indigenous chicken breeds (including two Chinese gamecock breeds), one red jungle fowl from Guangxi Province of China and three gamecock chickens from Laos. Combined with 31 published chicken genomes from three populations, a comparative genomics analysis was performed across 157 chickens. We found a severe confounding effect on potential cold adaptation exerted by introgression from commercial chickens into Chinese indigenous chickens, and argued that the genetic introgression from commercial chickens into indigenous chickens should be seriously considered for identifying selection footprint in indigenous chickens. LX gamecock chickens might have played a core role in recent breeding and conservation of other Chinese gamecock chickens. Importantly, AGMO (Alkylglycerol monooxygenase) and CPZ (Carboxypeptidase Z) might be crucial for determining the behavioral pattern of gamecock chickens, while ISPD (Isoprenoid synthase domain containing) might be essential for the muscularity of gamecock chickens. Our results can further the understanding of the evolution of Chinese gamecock chickens, especially the genetic basis of gamecock chickens revealed here was valuable for us to better understand the mechanisms underlying the behavioral pattern and the muscular development in chicken.

863 genomes reveal the origin and domestication of chicken

Despite the substantial role that chickens have played in human societies across the world, both the geographic and temporal origins of their domestication remain controversial. To address this issue, we analyzed 863 genomes from a worldwide sampling of chickens and representatives of all four species of wild jungle fowl and each of the five subspecies of red jungle fowl (RJF). Our study suggests that domestic chickens were initially derived from the RJF subspecies Gallus gallus spadiceus whose present-day distribution is predominantly in southwestern China, northern Thailand and Myanmar. Following their domestication, chickens were translocated across Southeast and South Asia where they interbred locally with both RJF subspecies and other jungle fowl species. In addition, our results show that the White Leghorn chicken breed possesses a mosaic of divergent ancestries inherited from other subspecies of RJF. Despite the strong episodic gene flow from geographically divergent lineages of jungle fowls, our analyses show that domestic chickens undergo genetic adaptations that underlie their unique behavioral, morphological and reproductive traits. Our study provides novel insights into the evolutionary history of domestic chickens and a valuable resource to facilitate ongoing genetic and functional investigations of the world's most numerous domestic animal.

Breeding history and candidate genes responsible for black skin of Xichuan black-bone chicken

Background

Domesticated chickens have a wide variety of phenotypes, in contrast with their wild progenitors. Unlike other chicken breeds, Xichuan black-bone chickens have blue-shelled eggs, and black meat, beaks, skin, bones, and legs. The breeding history and the economically important traits of this breed have not yet been explored at the genomic level. We therefore used whole genome resequencing to analyze the breeding history of the Xichuan black-bone chickens and to identify genes responsible for its unique phenotype.

Results

Principal component and population structure analysis showed that Xichuan black-bone chicken is in a distinct clade apart from eight other breeds. Linkage disequilibrium analysis showed that the selection intensity of Xichuan black-bone chickens is higher than for other chicken breeds. The estimated time of divergence between the Xichuan black-bone chickens and other breeds is 2.89 ka years ago. Fst analysis identified a selective sweep that contains genes related to melanogenesis. This region is probably associated with the black skin of the Xichuan black-bone chickens and may be the product of long-term artificial selection. A combined analysis of genomic and transcriptomic data suggests that the candidate gene related to the black-bone trait, EDN3, might interact with the upstream ncRNA LOC101747896 to generate black skin color during melanogenesis.

Conclusions

These findings help explain the unique genetic and phenotypic characteristics of Xichuan black-bone chickens, and provide basic research data for studying melanin deposition in animals.

Identifying Candidate Genes for Hypoxia Adaptation of Tibet Chicken Embryos by Selection Signature Analyses and RNA Sequencing

The Tibet chicken (Gallus gallus) lives on the Qinghai-Tibet Plateau and adapts to the hypoxic environment very well. The objectives of this study was to obtain candidate genes associated with hypoxia adaptation in the Tibet chicken embryos. In the present study, we used the fixation index (Fst) and cross population extended haplotype homozygosity (XPEHH) statistical methods to detect signatures of positive selection of the Tibet chicken, and analyzed the RNA sequencing data from the embryonic liver and heart with HISAT, StringTie and Ballgown for differentially expressed genes between the Tibet chicken and White leghorn (Gallus gallus, a kind of lowland chicken) embryos hatched under hypoxia condition. Genes which were screened out by both selection signature analysis and RNA sequencing analysis could be regarded as candidate genes for hypoxia adaptation of chicken embryos. We screened out 1772 genes by XPEHH and 601 genes by Fst, and obtained 384 and 353 differentially expressed genes in embryonic liver and heart, respectively. Among these genes, 89 genes were considered as candidate genes for hypoxia adaptation in chicken embryos. ARNT, AHR, GSTK1 and FGFR1 could be considered the most important candidate genes. Our findings provide references to elucidate the molecular mechanism of hypoxia adaptation in Tibet chicken embryos.

Efficiency of using SNP markers in the <i>MSTN</i> gene in the selection of the Pushkin breed chickens

In the poultry industry, indicators reflecting the growth rate of young stock and the exterior characteristics of chickens are important benchmarks for breeding. Traditional selection based on phenotypic evaluation is characterized by low efficiency with a low character inheritance ratio and is difficult to apply in small groups of animals and birds bred in bioresource collections. The use of molecular genetic markers associated with economically important traits makes it possible to carry out early selection of birds. This entails an increase in the profitability of the poultry industry. Recently, single nucleotide polymorphisms (SNPs) have served as convenient markers for selection purposes. For five generations (P1–P5), an experimental selection of hens of the Pushkin breed was carried out for live weight. It was based on selection for single nucleotide polymorphism rs313744840 in the MSTN gene. As a result, a significant increase in the frequency of allele A in this gene, from 0.11 to 0.50, took place. The association of SNP markers with meat qualities in the experimental group led to changes in the exterior profile of an adult bird at 330 days of age. The individuals with the AA and AG genotypes had the greatest live weight and longest body. As a result of selection, the bird on average became larger due to an increase in the number of heterozygous individuals with long bodies and large chest girths. The depth of the chest and the width of the pelvis increased due to an increase in the frequency of allele A in the experimental population. A tendency towards an increase in these indicators with the substitution of G with A in the genotype was found. Saturation of the population with desirable alleles led to an increase in the average live weight of the chickens. Analysis of the exterior parameters of adult birds showed that this growth is achieved by increasing the depth and volume of the bird body, and not by increasing the length of the limbs. Thus, marker selection carried out for five generations in the experimental population of Pushkin breed chickens to increase body weight has reliably (p < 0.001) changed the exterior profile of adult birds.

Genome-wide re-sequencing and transcriptome analysis reveal candidate genes associated with the pendulous comb phenotype in domestic chickens

To determine the causative variations associated with two chicken comb phenotypes, pendulous comb (PC) or upright comb (UC), two pooled genomic DNA samples from PC and UC chickens were re-sequenced by Next-Generation Sequencer, and genome-wide Single nucleotide polymorphisms (SNPs) were detected. Using three selective sweep approaches, F_ST , θπ, and Tajima's D, with top 5% window values serving as the threshold, a total of 84 positively selective genes (PSGs) were identified. There were no SNPs in exons of the PSGs with significant differences in allele frequencies between the two comb phenotype groups. Then, 515 differentially expressed genes (DEGs) between the PC and UC were identified by RNA-seq. Three genes including CD36 (CD36 molecule), ADAMTSL3 (ADAMTS-like 3), and AOX1 (aldehyde oxidases 1) are overlapped between PSGs and DEGs. After genotyping seven candidate SNPs in the regulatory regions of the three overlapping genes in 120 chickens from two other breeds, two variants (rs14607046 and rs731818051) in the regulatory regions of AOX1 and ADAMTSL3 were found to have significant differences in allele frequency between the PC and UC, suggesting that the two variants may be causative mutations for PC. Overall, our study shed light on the genetic basis underlying the PC phenotype in chickens.

Feather colour affects the aggressive behaviour of chickens with the same genotype on the dominant white (I) locus

Aggression in chickens is a serious economic and animal welfare issue in poultry farming. Pigmentation traits have been documented to be associated with animal behaviour. Chicken pecking behaviour has been found to be related to feather colour, with premelanosome protein 17 (PMEL17) being one of the candidate genes. In the present study, we performed a genotypic and phenotypic association analysis between chicken plumage colour (red and white) and aggressive behaviour in an F1 hybrid group generated by crossing the autosomal dominant white-feathered breed White Leghorn (WL) and the red-feathered breed Rhode Island Red (RIR). In genetic theory, all the progeny should have white feathers because WL is homozygous autosomal dominant for white feathers. However, we found a few red-feathered female chickens. We compared the aggressiveness between the red and white females to determine whether the feather colour alone affected the behaviour, given that the genetic background should be the same except for feather colour. The aggressiveness was recorded 5 days after sexual maturity at 26 weeks. Generally, white plumage hens showed significantly higher aggressiveness compared to the red ones in chasing, attacking, pecking, and threatening behaviour traits. We assessed three candidate feather colour genes—PMEL17, solute carrier family 45 member 2 (SLC45A2), and SRY-box 10 (SOX10)—to determine the genetic basis for the red and white feather colour in our hybrid population; there was no association between the three loci and feather colour. The distinct behavioural findings observed herein provide clues to the mechanisms underlying the association between phenotype and behaviour in chickens. We suggest that mixing red and white chickens together might reduce the occurrence of aggressive behaviours.

Population genomics identifies patterns of genetic diversity and selection in chicken

BACKGROUND

There are hundreds of phenotypically distinguishable domestic chicken breeds or lines with highly specialized traits worldwide, which provide a unique opportunity to illustrate how selection shapes patterns of genetic variation. There are many local chicken breeds in China.

RESULTS

Here, we provide a population genome landscape of genetic variations in 86 domestic chickens representing 10 phenotypically diverse breeds. Genome-wide analysis indicated that sex chromosomes have less genetic diversity and are under stronger selection than autosomes during domestication and local adaptation. We found an evidence of admixture between Tibetan chickens and other domestic population. We further identified strong signatures of selection affecting genomic regions that harbor genes underlying economic traits (typically related to feathers, skin color, growth, reproduction and aggressiveness) and local adaptation (to high altitude). By comparing the genomes of the Tibetan and lowland fowls, we identified genes associated with high-altitude adaptation in Tibetan chickens were mainly involved in energy metabolism, body size maintenance and available food sources.

CONCLUSIONS

The work provides crucial insights into the distinct evolutionary scenarios occurring under artificial selection for agricultural production and under natural selection for success at high altitudes in chicken. Several genes were identified as candidates for chicken economic traits and other phenotypic traits.

Population genomic data reveal genes related to important traits of quail

Background

Japanese quail (Coturnix japonica), a recently domesticated poultry species, is important not only as an agricultural product, but also as a model bird species for genetic research. However, most of the biological questions concerning genomics, phylogenetics, and genetics of some important economic traits have not been answered. It is thus necessary to complete a high-quality genome sequence as well as a series of comparative genomics, evolution, and functional studies.

Results

Here, we present a quail genome assembly spanning 1.04 Gb with 86.63% of sequences anchored to 30 chromosomes (28 autosomes and 2 sex chromosomes Z/W). Our genomic data have resolved the long-term debate of phylogeny among Perdicinae (Japanese quail), Meleagridinae (turkey), and Phasianinae (chicken). Comparative genomics and functional genomic data found that four candidate genes involved in early maturation had experienced positive selection, and one of them encodes follicle stimulating hormone beta (FSHβ), which is correlated with different FSHβ levels in quail and chicken. We re-sequenced 31 quails (10 wild, 11 egg-type, and 10 meat-type) and identified 18 and 26 candidate selective sweep regions in the egg-type and meat-type lines, respectively. That only one of them is shared between egg-type and meat-type lines suggests that they were subject to an independent selection. We also detected a haplotype on chromosome Z, which was closely linked with maroon/yellow plumage in quail using population resequencing and a genome-wide association study. This haplotype block will be useful for quail breeding programs.

Conclusions

This study provided a high-quality quail reference genome, identified quail-specific genes, and resolved quail phylogeny. We have identified genes related to quail early maturation and a marker for plumage color, which is significant for quail breeding. These results will facilitate biological discovery in quails and help us elucidate the evolutionary processes within the Phasianidae family.

Animal Morality: What It Means and Why It Matters

It has been argued that some animals are moral subjects, that is, beings who are capable of behaving on the basis of moral motivations (Rowlands 2011, 2012, 2017). In this paper, we do not challenge this claim. Instead, we presuppose its plausibility in order to explore what ethical consequences follow from it. Using the capabilities approach (Nussbaum 2004, 2007), we argue that beings who are moral subjects are entitled to enjoy positive opportunities for the flourishing of their moral capabilities, and that the thwarting of these capabilities entails a harm that cannot be fully explained in terms of hedonistic welfare. We explore the implications of this idea for the assessment of current practices involving animals.

Population structure, genetic diversity and phylogenetic analysis of different rural and commercial chickens of Pakistan using complete sequence of mtDNA D-loop

Chicken is the most important poultry bird in Pakistan that not only provides nutrition but also contributes to country's economy. The Pakistani chicken and its germplasm resources are not genetically characterized and explored. Here, we focused at rural and commercial chickens of Pakistan and explored genetic diversity, population structure and phylogeny. We first collected feather samples from Rural and Broiler populations of Pakistani chickens, isolated DNA and sequenced complete D-loop of mtDNA. The length of complete D-loop ranged from 1231 to 1234 bp in Pakistani chickens. The GC content was 39%. Hotspots of mutations were three hypervariable sites (HVS). Most of the variations (77%) were in HVS1. In a total, 26 polymorphic sites defined 12 haplotypes and all major haplogroups (A-I) in genetic structure of Pakistani chickens. Genetic diversity remained relatively very low in Broiler (Pi = 0.00212 ± 0.00136). There was a low sharing of matrilineages between the two populations (Fst = 0.170). With high Hd value (0.825 ± 0.051) and presence of all nine major haplogroups the rural chicken population showed relatively rich genepool. Finally we did molecular phylogenetic analysis and inferred phylogeny. Presence of subcontinent specific haplogroups E3 and I and clustering of Indian red junglefowl closely with Pakistani chickens in Bayesian inference tree, provide further evidence for an independent domestication event of chicken in subcontinent.

Exploring evidence of positive selection signatures in cattle breeds selected for different traits

Since domestication, the genome landscape of cattle has been changing due to natural and artificial selection forces resulting in several general and specialized cattle breeds of the world. Identifying genomic regions affected due to these forces in livestock gives an insight into the history of selection for economically important traits and genetic adaptation to specific environments of the populations under consideration. This study explores the genes/genomic regions under selection in relation to the phenotypes of Holstein, Hanwoo, and N'Dama cattle breeds using Tajima's D, XP-CLR, and XP-EHH population statistical methods. The whole genomes of 10 Holstein (South Korea), 11 Hanwoo (South Korea), and 10 N'Dama (West Africa-Guinea) cattle breeds re-sequenced to ~11x coverage and retained 37 million SNPs were used for the study. Selection signature analysis revealed 441, 512, and 461 genes under selection from Holstein, Hanwoo, and N'Dama cattle breeds, respectively. Among all these, seven genes including ARFGAP3, SNORA70, and other RNA genes were common between the breeds. From each of the gene lists, significant functional annotation cluster terms including milk protein and thyroid hormone signaling pathway (Holstein), histone acetyltransferase activity (Hanwoo), and renin secretion (N'Dama) were enriched. Genes that are related to the phenotypes of the respective breeds were also identified. Moreover, significant breed-specific missense variants were identified in CSN3, PAPPA2 (Holstein), C1orf116 (Hanwoo), and COMMD1 (N'Dama) genes. The genes identified from this study provide an insight into the biological mechanisms and pathways that are important in cattle breeds selected for different traits of economic significance.

Limited awareness of animal influenza prevention and control among Dai Lue smallholder farmers in Southwest China

Awareness of animal influenza and its prevention and control is important for ensuring livestock health, production and welfare. In China, a country stereotyped as a major source of emerging zoonotic infectious diseases, research on the public understanding of animal influenza is limited to the Han, the main ethnic group. The present qualitative study in Southwest China investigated awareness of animal influenza among the Dai, an ethnic minority. The participants (15 men and 10 women, ages 18–83) were smallholder farmers of pigs and poultry in rural areas of Jinghong, Xishuangbanna, Yunnan Province. A mixture of interviews and group discussions took place in homes and villages. The participants were asked about their knowledge of avian influenza (H7N9), swine influenza (H1N1), precautions taken to protect against influenza, procedures when animals were sick and perceived risk of animal influenza. The data were analysed following coding and thematic analysis. The findings demonstrated a limited understanding of animal health and welfare among participants. Specifically, they were largely unaware of animal influenza (H7N9, H1N1) including its causes, symptoms, prevention and treatment. The farmers were also uninformed of the risks they faced and unknowingly engaged in behaviours which increased direct or indirect exposure to infected animals, a risk factor for human infection. They also reported poor usage of veterinary services. In order to guarantee the health, welfare and production of their livestock, immediate action is needed to enable Dai smallholder farmers to prevent and respond to animal influenza effectively and timely.

Genomic analysis reveals selection in Chinese native black pig

Identification of genomic signatures that help reveal mechanisms underlying desirable traits in domesticated pigs is of significant biological, agricultural and medical importance. To identify the genomic footprints left by selection during domestication of the Enshi black pig, a typical native and meat-lard breed in China, we generated about 72-fold coverage of the pig genome using pools of genomic DNA representing three different populations of Enshi black pigs from three different locations. Combining this data with the available whole genomes of 13 Chinese wild boars, we identified 417 protein-coding genes embedded in the selected regions of Enshi black pigs. These genes are mainly involved in developmental and metabolic processes, response to stimulus, and other biological processes. Signatures of selection were detected in genes involved in body size and immunity (RPS10 and VASN), lipid metabolism (GSK3), male fertility (INSL6) and developmental processes (TBX19). These findings provide a window into the potential genetic mechanism underlying development of desirable phenotypes in Enshi black pigs during domestication and subsequent artificial selection. Thus, our results illustrate how domestication has shaped patterns of genetic variation in Enshi black pigs and provide valuable genetic resources that enable effective use of pigs in agricultural production.