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

Comparative Analysis of Runs of Homozygosity Islands in Indigenous and Commercial Chickens Revealed Candidate Loci for Disease Resistance and Production Traits

Runs of homozygosity (ROH) are contiguous stretches of identical genomic regions inherited from both parents. Assessment of ROH in livestock species contributes significantly to our understanding of genetic health, population genetic structure, selective pressure and conservation efforts. In this study, whole genome re-sequencing data from 140 birds of 10 Iranian indigenous chicken ecotypes, 3 commercial chicken breeds and 1 red junglefowl (RJF) population were used to investigate their population genetic structure, ROH characteristics (length and frequency) and genomic inbreeding coefficients (FROH). Additionally, we examined ROH islands for selection footprints in the indigenous chicken populations. Our results revealed distinct genetic backgrounds, among which the White Leghorn breed exhibited the greatest genetic distance from other populations, while the gamecock populations formed a separate cluster. We observed significant differences in ROH characteristics, in which the commercial breeds showed a higher number of ROH compared to indigenous chickens and red junglefowls. Short ROH ranging from 0.1 to 1 Mb were dominant among the populations. The Arian line had the highest mean length of ROH, while the White Leghorn breed showed the highest number of ROH. Among indigenous chickens, the Lari-Afghani ecotype exhibited the highest FROH, but the Sari inherited the richest genetic diversity. Interestingly, GGA16 carried no ROH in the red junglefowls, whereas GGA22 had the highest FROH across all populations, except in the Isfahan ecotype. We also identified ROH islands associated with genetic adaptations in indigenous ecotypes. These islands harboured immune-related genes contributing to disease resistance (TLR2, TICAM1, IL22RA1, NOS2, CCL20 and IFNLR1), heat tolerance and oxidative stress response (NFKB1, HSF4, OSGIN1 and BDNF), and muscle development, lipid metabolism and reproduction (MEOX2, CEBPB, CDS2 and GnRH-I). Overall, this study highlights the genetic potential of indigenous chickens to survive and adapt to their respective environments.

Genome-Wide Scan for Selective Sweeps Reveals Novel Loci Associated with Prolificacy in Iranian Sheep Breeds in Comparison with Highly Prolific Exotic Breed

Domestication and selection significantly changed phenotypic traits in modern domestic animals. To identify the genomic regions associated with prolificacy in this study, 837 ewes from three Iranian indigenous sheep breeds, consisting of Baluchi, Lori-Bakhtiari, and Zandi uniparous breeds, and one Greek highly prolific dairy sheep, namely Chios, were genotyped using OvineSNP50K arrays. Statistical tests were then performed using different and complementary methods based on either site frequency (FST) and haplotype (hapFLK) between populations, followed by a pathway analysis of the genes contained in the selected regions. The results revealed that for the top 0.01 percentile of the obtained FST values, 16 genomic regions on chromosomes 2, 3, 4, 7, 8, 9, 13, 14, 16, 18, 19, and 20, and for hapFLK values, 3 regions located on chromosomes 3, 7, and 13, were under selection. A bioinformatic analysis of these genomic regions showed that these loci overlapped with potential candidate genes associated with prolificacy in sheep including GNAQ, COL5A2, COL3A1, HECW1, FBN1, COMMD3, RYR1, CCL28, SERPINA14, and HSPA2. These regions also overlapped with some quantitative trait loci (QTLs) linked to prolificacy traits, milk yield, and body weight. These findings suggest that future research could further link these genomic regions to prolificacy traits in sheep.

Danzhou chicken: a unique genetic resource revealed by genome-wide resequencing data

Danzhou chicken (DZ) is a local breed in China noted for its strong adaptability, roughage resistance, strong wildness, and delicious taste, thus containing important genetic resources. In this study, genome re-sequencing data was generated from 200 DZ chickens. Combined with previously generated data from 72 additional chickens across six other exotic and local breeds, these data were used to systematically evaluate the germplasm characteristics of DZ chickens from a genomic perspective. Unlike exotic breeds, both DZ and southern local chicken varieties exhibited high genetic diversity, and the genetic distance between DZ and southern local chickens was smaller than the genetic distance between DZ and exotic chickens. A reconstructed Neighbor-Joining phylogenetic tree indicated that all sampled populations clustered into single independent populations, with DZ chickens showing clear evidence of intra-population differentiation, forming 2 subpopulations. Principal component analysis and ADMIXTURE analysis showed that DZ was significantly different from other breeds. These results indicate that DZ is a unique genetic resource that is different from other southern native and exotic chickens. The results of the study will improve our understanding of the genetic structure and current status of the DZ breed, which is of great significance in promoting the conservation of genetic resources of DZ chickens and fostering breed innovations and genetic improvement.

Integrated transcriptomic analysis on chicken ovary reveals CYP21A1 affects follicle granulosa cell development and steroid hormone synthesis

Egg production is an economically important trait in poultry breeding and production. Follicular development was regulated by several hormones released and genes expressed in the granulosa cells, impacting the egg production and fecundity of hens. However, the molecular functions of these candidate genes that modulate these processes remain largely unknown. In the present study, bioinformatics analyses were performed to identify the candidate genes related to egg production in the ovarian tissue of White Leghorns with high egg production and Beijing You chicken with low egg production during sexual maturity and peak laying periods. The ovarian granulosa cells were used to assess the function of CYP21A1 by transfecting with CYP21A1-specific small interfering RNAs (siRNAs) and overexpression plasmids. We identified 514 differentially expressed genes (|Log2(fold change) | >1, P <0.05) between the 2 chicken breeds in both laying periods. Among these genes, CYP21A1, which is involved in the steroid hormone biosynthesis pathway was consistently upregulated in White Leghorns. Weighted gene co-expression network analysis (WGCNA) further suggested that CYP21A1 was a hub gene, which could positively respond to treatment with follicle stimulation hormone (FSH), affecting egg production. The interference of CYP21A1 significantly inhibited cell proliferation and promoted cell apoptosis. Overexpression of CYP21A1 promotes cell proliferation and inhibits cell apoptosis. Furthermore, the interference with CYP21A1 significantly downregulated the expression of STAR, CYP11A1, HSD3B1, and FSHR and also decreased the synthesis of progesterone (P4) and estradiol (E2) in granulosa cells. Overexpression of CYP21A1 increased the synthesis of P4 and estradiol E2 and the expression of steroid hormone synthesis-related genes in granulosa cells. Our findings provide new evidence for the biological role of CYP21A1 on granulosa cell proliferation, apoptosis, and steroid hormone synthesis, which lays the theoretical basis for improving egg production.

Genome-wide scans for selection signatures in indigenous chickens reveal candidate genes associated with local adaptation

Population growth and climate change pose challenges to the sustainability of poultry farming. The emphasis on high-yield traits in commercialized breeds has led to a decline in their adaptability. Chicken varieties adapted to the local environment, possessing traits that facilitate adaptation to climate change, such as disease resistance and tolerance to extreme weather conditions, can improve hybridization outcomes. In this study, we conducted an analysis of the population structure and genetic diversity of 110 chickens representing indigenous breeds from southern China and two different commercial breeds. Further, we performed comparative population genomics, utilizing nucleotide diversity and fixation statistics, to characterize genomic features of natural selection and to identify unique genetic traits and potential selection markers developed by indigenous breeds after adapting to the local environment. Results based on genetic diversity and population structure analyses showed that indigenous varieties exhibited high levels of genetic diversity. Commercial breeds that have been indigenously bred demonstrated higher levels of genetic diversity than those that have not, and breeds with different selection practices displayed significant differences in genetic structure. Additionally, we further searched for potential genomic regions in native chicken ecotypes, uncovering several candidate genes related to ecological adaptations affecting local breeds, such as IKBKB, S1PR1, TSHR, IL1RAPL1 and AMY2A, which are involved in disease resistance, heat tolerance, immune regulation and behavioral traits. This work provides important insights into the genomic characterization of ecotypes of native chicken in southern China. The identification of candidate genes associated with traits such as disease resistance, heat tolerance, immunomodulation, and behavioral traits is a significant outcome. These candidate genes may contribute to the understanding of the molecular basis of the adaptation of the southern native chicken to the local environment. It is recommended that these genes be integrated into chicken breeding programs to enhance sustainable agriculture and promote effective conservation and utilization strategies.

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.

Genomic Inbreeding and Runs of Homozygosity Analysis of Cashmere Goat

Cashmere goats are valuable genetic resources which are famous worldwide for their high-quality fiber. Runs of homozygosity (ROHs) have been identified as an efficient tool to assess inbreeding level and identify related genes under selection. However, there is limited research on ROHs in cashmere goats. Therefore, we investigated the ROH pattern, assessed genomic inbreeding levels and examined the candidate genes associated with the cashmere trait using whole-genome resequencing data from 123 goats. Herein, the Inner Mongolia cashmere goat presented the lowest inbreeding coefficient of 0.0263. In total, we identified 57,224 ROHs. Seventy-four ROH islands containing 50 genes were detected. Certain identified genes were related to meat, fiber and milk production (FGF1, PTPRM, RERE, GRID2, RARA); fertility (BIRC6, ECE2, CDH23, PAK1); disease or cold resistance and adaptability (PDCD1LG2, SVIL, PRDM16, RFX4, SH3BP2); and body size and growth (TMEM63C, SYN3, SDC1, STRBP, SMG6). 135 consensus ROHs were identified, and we found candidate genes (FGF5, DVL3, NRAS, KIT) were associated with fiber length or color. These findings enhance our comprehension of inbreeding levels in cashmere goats and the genetic foundations of traits influenced by selective breeding. This research contributes significantly to the future breeding, reservation and use of cashmere goats and other goat breeds.

Genotyping-by-Sequencing Strategy for Integrating Genomic Structure, Diversity and Performance of Various Japanese Quail (Coturnix japonica) Breeds

Traces of long-term artificial selection can be detected in genomes of domesticated birds via whole-genome screening using single-nucleotide polymorphism (SNP) markers. This study thus examined putative genomic regions under selection that are relevant to the development history, divergence and phylogeny among Japanese quails of various breeds and utility types. We sampled 99 birds from eight breeds (11% of the global gene pool) of egg (Japanese, English White, English Black, Tuxedo and Manchurian Golden), meat (Texas White and Pharaoh) and dual-purpose (Estonian) types. The genotyping-by-sequencing analysis was performed for the first time in domestic quails, providing 62,935 SNPs. Using principal component analysis, Neighbor-Net and Admixture algorithms, the studied breeds were characterized according to their genomic architecture, ancestry and direction of selective breeding. Japanese and Pharaoh breeds had the smallest number and length of homozygous segments indicating a lower selective pressure. Tuxedo and Texas White breeds showed the highest values of these indicators and genomic inbreeding suggesting a greater homozygosity. We revealed evidence for the integration of genomic and performance data, and our findings are applicable for elucidating the history of creation and genomic variability in quail breeds that, in turn, will be useful for future breeding improvement strategies.

Whole Genome Resequencing Revealed the Genetic Relationship and Selected Regions among Baicheng-You, Beijing-You, and European-Origin Broilers

As the only two You-chicken breeds in China, Baicheng-You (BCY) and Beijing-You (BJY) chickens are famous for their good meat quality. However, so far, the molecular basis of germplasm of the two You-chicken breeds is not yet clear. The genetic relationship among BCY, BJY, and European-origin broilers (BRs) was analyzed using whole genome resequencing data to contribute to this issue. A total of 18,852,372 single nucleotide polymorphisms (SNPs) were obtained in this study. After quality control, 8,207,242 SNPs were applied to subsequent analysis. The data indicated that BJY chickens possessed distant distance with BRs (genetic differentiation coefficient (FST) = 0.1681) and BCY (FST = 0.1231), respectively, while BCY and BRs had a closer relationship (FST = 0.0946). In addition, by using FST, cross-population extended haplotype homozygosity (XP-EHH), and cross-population composite likelihood ratio (XP-CLR) methods, we found 374 selected genes between BJY and BRs chickens and 279 selected genes between BCY and BJY chickens, respectively, which contained a number of important candidates or genetic variations associated with feather growth and fat deposition of BJY chickens and potential disease resistance of BCY chickens. Our study demonstrates a genome-wide view of genetic diversity and differentiation among BCY, BJY, and BRs. These results may provide useful information on a molecular basis related to the special characteristics of these broiler breeds, thus enabling us to better understand the formation mechanism of Chinese-You chickens.

Comparative population genomics analysis uncovers genomic footprints and genes influencing body weight trait in Chinese indigenous chicken

Body weight of chicken is a typical quantitative trait, which shows phenotypic variations due to selective breeding. Despite some QTL loci have been obtained, the body weight of native chicken breeds in different geographic regions varies greatly, its genetic basis remains unresolved questions. To address this issue, we analyzed 117 Chinese indigenous chickens from 10 breeds (Huiyang Bearded, Xinhua, Hotan Black, Baicheng You, Liyang, Yunyang Da, Jining Bairi, Lindian, Beijing You, Tibetan). We applied fixation index (FST) analysis to find selected genomic regions and genes associated with body weight traits. Our study suggests that NELL1, XYLT1, and NCAPG/LCORL genes are strongly selected in the body weight trait of Chinese indigenous chicken breeds. In addition, the IL1RAPL1 gene was strongly selected in large body weight chickens, while the PCDH17 and CADM2 genes were strongly selected in small body weight chickens. This result suggests that the patterns of genetic variation of native chicken and commercial chicken, and/or distinct local chicken breeds may follow different evolutionary mechanisms.

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.

Genome-wide mapping of signatures of selection using a high-density array identified candidate genes for growth traits and local adaptation in chickens

BACKGROUND

Availability of single nucleotide polymorphism (SNP) genotyping arrays and progress in statistical analyses have allowed the identification of genomic regions and genes under selection in chicken. In this study, SNP data from the 600 K Affymetrix chicken array were used to detect signatures of selection in 23 local Italian chicken populations. The populations were categorized into four groups for comparative analysis based on live weight (heavy vs light) and geographical area (Northern vs Southern Italy). Putative signatures of selection were investigated by combining three extended haplotype homozygosity (EHH) statistical approaches to quantify excess of haplotype homozygosity within (iHS) and between (Rsb and XP-EHH) groups. Presence of runs of homozygosity (ROH) islands was also analysed for each group.

RESULTS

After editing, 541 animals and 313,508 SNPs were available for statistical analyses. In total, 15 candidate genomic regions that are potentially under selection were detected among the four groups: eight within a group by iHS and seven by combining the results of Rsb and XP-EHH, which revealed divergent selection between the groups. The largest overlap between genomic regions identified to be under selection by the three approaches was on chicken chromosome 8. Twenty-one genomic regions were identified with the ROH approach but none of these overlapped with regions identified with the three EHH-derived statistics. Some of the identified regions under selection contained candidate genes with biological functions related to environmental stress, immune responses, and disease resistance, which indicate local adaptation of these chicken populations.

CONCLUSIONS

Compared to commercial lines, local populations are predominantly reared as backyard chickens, and thus, may have developed stronger resistance to environmental challenges. Our results indicate that selection can play an important role in shaping signatures of selection in local chicken populations and can be a starting point to identify gene mutations that could have a useful role with respect to climate change.

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.

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.

A genome-wide scan for signatures of selection in Kurdish horse breed

The genetic diversity and genomic regions being under putative natural selection in Kurdish horse population were studied. The samples from 72 horses were genotyped by using GGP Equine 70K SNP arrays. The Ne Slope (NeS) analyses revealed that a sharp decline in Ne has probably occurred around four generations ago, and high frequency of ROH with 2-4 Mbp in length suggested that the inbreeding has probably occurred around 20 generations ago. The effective population size (Ne) was 104 horses up to three generations ago and the average inbreeding (F_ROH) was 0.047(±0.045). Using de-correlated composite of multiple selection signals (DCMS) and runs of homozygosity (ROH) analyses the genomic regions being under putative selection were detected. By using DCMS, a total of 148 significant SNP (FDR < 0.05) were identified, 40% of which were located on ECA9, where the greatest peak was observed. This genomic region harbors several known QTL which are associated with withers height (body size). Also significant genomic regions (FDR<0.05), harboring QTL associated with insect bite hypersensitivity (IBH), hair density and coat texture, alternate gaits, guttural pouch tympany and temperament were identified. By using outputs of ROH analyses, two hotspot regions (i.e. 30% of individuals was considered as threshold), were identified on ECA7 (50.11-54.36 Mbp) and ECA11 (26.10-29.07 Mbp) harboring QTL associated with withers height, alternate gait and IBH. In summary, the genomic regions being under putative natural selection which harbors known QTL associated with body size and IBH, among others, were introduced. Nevertheless, additional functional and comparative studies are necessary to corroborate their effect on the observed genetic and phenotypic diversity of the Kurdish horses.

The SLC27A1 Gene and Its Enriched PPAR Pathway Are Involved in the Regulation of Flavor Compound Hexanal Content in Chinese Native Chickens

The role of hexanal in flavor as an indicator of the degree of oxidation of meat products is undeniable. However, the genes and pathways of hexanal formation have not been characterized in detail. In this study, we performed differential gene expression analysis and weighted gene co-expression network analysis (WGCNA) on groups of Tiannong partridge chickens with different relative hexanal content in order to find the genes involved in the formation of hexanal and the specific pathways of hexanal formation. Then we confirmed the relationship of these candidate genes with hexanal using Jingxing Yellow chicken and Wenchang chicken. In this study, WGCNA revealed a module of co-expressed genes that were highly associated with the volatile organic compound hexanal. We also compared transcriptome gene expression data of samples from chicken groups with high and low relative contents of hexanal and identified a total of 651 differentially expressed genes (DEGs). Among them, 356 genes were up regulated, and 295 genes were downregulated. The different biological functions associated with the DEGs, hub genes and hexanal were identified by functional analysis using the Kyoto Encyclopedia of Genes and Genomes (KEGG) annotations. Among all the hub genes in the significant module identified by WGCNA, more were enriched in the PPAR signaling pathway, the proteasome pathway, etc. Additionally, we found that DEGs and hub genes, including SLC27A1, ACOX3, NR4A1, VEGFA, JUN, EGR1, CACNB1, GADD45A and DUSP1, were co-enriched in the peroxisome proliferator-activated receptor (PPAR) signaling pathway, p53 signaling pathway and mitogen-activated protein kinases (MAPK) signaling pathway, etc. Transcriptome results of the Jingxing Yellow chicken population showed that the SLC27A1 gene was significantly associated with hexanal and enriched in the PPAR pathway. Our study provides a comprehensive insight into the key genes related to hexanal content, and can be further explored by functional and molecular studies.