A combination of genome-wide association study and selection signature analysis dissects the genetic architecture underlying bone traits in chickens

Research Note: SOCS2 contributes to reduction of the third digit during development of the chicken forelimb

The development of the avian wing pattern has been the subject of heated debate due to its special shape. The Suppressor of cytokine signaling 2 (SOCS2) gene encodes a negative regulator of growth hormone (GH) signaling and bone growth and is known to be strongly expressed in the third digit of chicken forelimbs. These observations suggest that SOCS2 might regulate the morphology of the avian wing, however, the function of SOCS2 in avian limb development remains unknown. Here, we reexamined SOCS2 expression in successive developmental stages of chicken limb development by in situ hybridization (ISH) and describe extended expression from the posterior of the stypolod to the third digit of the forelimbs. We used the RCAS avian retrovirus to overexpress SOCS2 in the developing chicken limb buds, which resulted in reduced or malformed chicken wings while hindlimbs developed normally. Transcriptome sequencing (mRNA-Seq) revealed changes in expression of genes known to be associated with growth and development in forelimbs with overexpressed SOCS2. This study highlights a pivotal role for SOCS2 during the development of the wing in the chicken and provides new insight into molecular mechanisms regulating avian limb development.

Genome-wide association study for bone quality of ducks during the laying period

The cage-rearing model of the modern poultry industry makes the bones of birds, especially egg-laying birds, more vulnerable to fracture, which poses serious damage to the health of birds. Research confirms that genetic material plays an important role in regulating bone growth, development, and remodeling. However, the genetic architecture underlying bone traits is not well understood. The objectives of this study are to identify valuable genes and genetic markers through a genome-wide association study (GWAS) for breeding to improve the duck bone quality. First, we quantified the tibia and femur quality traits of 260 laying ducks. Based on GWAS, a total of 75 SNP loci significantly associated with bone quality traits were identified, and 67 potential candidate genes were annotated. According to gene function analysis, genes P4HA2, WNT3A, and BST1 et al may influence bone quality by regulating bone cell activity, calcium and phosphate metabolism, or bone collagen maturation and cross-linking. Meanwhile, combined with the transcriptome results, we found that HOXB cluster genes are also important in bone growth and development. Therefore, our findings were helpful in further understanding the genetic architecture of the duck bone quality and provided a worthy theoretical basis and technological support to improve duck bone quality by breeding.

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.

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

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

Genome-wide association study of growth curve parameters reveals novel genomic regions and candidate genes associated with metatarsal bone traits in chickens

The growth and development of chicken bones have an enormous impact on the health and production performance of chickens. However, the development pattern and genetic regulation of the chicken skeleton are poorly understood. This study aimed to evaluate metatarsal bone growth and development patterns in chickens via non-linear models, and to identify the genetic determinants of metatarsal bone traits using a genome-wide association study (GWAS) based on growth curve parameters. Data on metatarsal length (MeL) and metatarsal circumference (MeC) were obtained from 471 F2 chickens (generated by crossing broiler sires, derived from a line selected for high abdominal fat, with Baier layer dams) at 4, 6, 8, 10, and 12 weeks of age. Four non-linear models (Gompertz, Logistic, von Bertalanffy, and Brody) were used to fit the MeL and MeC growth curves. Subsequently, the estimated growth curve parameters of the mature MeL or MeC (A), time-scale parameter (b), and maturity rate (K) from the non-linear models were utilized as substitutes for the original bone data in GWAS. The Logistic and Brody models displayed the best goodness-of-fit for MeL and MeC, respectively. Single-trait and multi-trait GWASs based on the growth curve parameters of the Logistic and Brody models revealed 4 618 significant single nucleotide polymorphisms (SNPs), annotated to 332 genes, associated with metatarsal bone traits. The majority of these significant SNPs were located on Gallus gallus chromosome (GGA) 1 (167.433-176.318 Mb), GGA2 (96.791-103.543 Mb), GGA4 (65.003-83.104 Mb) and GGA6 (64.685-95.285 Mb). Notably, we identified 12 novel GWAS loci associated with chicken metatarsal bone traits, encompassing 35 candidate genes. In summary, the combination of single-trait and multi-trait GWASs based on growth curve parameters uncovered numerous genomic regions and candidate genes associated with chicken bone traits. The findings benefit an in-depth understanding of the genetic architecture underlying metatarsal growth and development in chickens.

Weighted single-step GWAS identified candidate genes associated with carcass traits in a Chinese yellow-feathered chicken population

Carcass traits in broiler chickens are complex traits that are influenced by multiple genes. To gain deeper insights into the genetic mechanisms underlying carcass traits, here we conducted a weighted single-step genome-wide association study (wssGWAS) in a population of Chinese yellow-feathered chicken. The objective was to identify genomic regions and candidate genes associated with carcass weight (CW), eviscerated weight with giblets (EWG), eviscerated weight (EW), breast muscle weight (BMW), drumstick weight (DW), abdominal fat weight (AFW), abdominal fat percentage (AFP), gizzard weight (GW), and intestine length (IL). A total of 1,338 broiler chickens with phenotypic and pedigree information were included in this study. Of these, 435 chickens were genotyped using a 600K single nucleotide polymorphism chip for association analysis. The results indicate that the most significant regions for 9 traits explained 2.38% to 5.09% of the phenotypic variation, from which the region of 194.53 to 194.63Mb on chromosome 1 with the gene RELT and FAM168A identified on it was significantly associated with CW, EWG, EW, BMW, and DW. Meanwhile, the 5 traits have a strong genetic correlation, indicating that the region and the genes can be used for further research. In addition, some candidate genes associated with skeletal muscle development, fat deposition regulation, intestinal repair, and protection were identified. Gene ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses suggested that the genes are involved in processes such as vascular development (CD34, FGF7, FGFR3, ITGB1BP1, SEMA5A, LOXL2), bone formation (FGFR3, MATN1, MEF2D, DHRS3, SKI, STC1, HOXB1, HOXB3, TIPARP), and anatomical size regulation (ADD2, AKT1, CFTR, EDN3, FLII, HCLS1, ITGB1BP1, SEMA5A, SHC1, ULK1, DSTN, GSK3B, BORCS8, GRIP2). In conclusion, the integration of phenotype, genotype, and pedigree information without creating pseudo-phenotype will facilitate the genetic improvement of carcass traits in chickens, providing valuable insights into the genetic architecture and potential candidate genes underlying carcass traits, enriching our understanding and contributing to the breeding of high-quality broiler chickens.

Integrated analysis of whole genome and transcriptome sequencing uncovers genetic differences between Zi goose and Xianghai flying goose

Zi goose is a famous indigenous breed originating from northeast China with high annual egg production. Xianghai flying goose is a composite breed and is bred by crosses of the wild swan goose and the Zi goose. Our previous study revealed significant differences in muscle fiber characteristics between the two populations. Here, we aimed to reveal the underlying genetic basis of the above phenotype differences through whole-genome and transcriptome analysis. A total of 20 blood samples (10 Zi geese and 10 Xianghai flying geese) were used for whole genome sequencing, and eight breast muscle tissue samples (four Zi geese and four Xianghai flying geese) were used for RNA sequencing. Using the FST and XP-EHH analysis, some highly differentiated genome regions annotated with egg production (RORB, WNT4, BMPR1B) and breast muscle development (WNT7B) between the two populations were detected. RNA-sequencing analysis revealed differentially expressed genes related to muscle development (IGF1, PAX7). Moreover, several genes were detected by both genome and transcriptome analysis, and some of them were reported to be associated with muscle growth (SLIT2, PREX1) and intramuscular fat (COL6A1). These findings will help researchers better understand the genetic basis related to egg production and muscle development in geese.

Study of genes polymorphisms in RANK/RANKL/OPG and WNT signaling pathways and their associations with bone parameters in broiler chicken

Limb problems are one of the most common problems with fast-growing meat-type chickens. Various bone abnormalities, which can lead to limping, bone weakness, or even fractures, bring overall discomfort to birds and a loss of production. Genetic aspects are often associated with these side effects on bone stability and are also cited as the dominant cause. These points to a close negative relationship of genetic selection for rapid growth with traits involved in bone integrity. Due to the assumption of an additive genetic background, improvements through genetic tools can be used. Our study is focused on selected genes of important signaling pathways for bone metabolism. We tried to detect polymorphisms that would show associations with selected bone parameters in a total of 48 broilers. Those were fast-growing Ross 308 hybrids and slow-growing Hubbard M22BxJA87A hybrids. The TNFRSF11A and WISP1 genes were tested. A total of fourteen polymorphisms were found, three of them were synonymous and five in the intron. In the case of four polymorphisms found in exons of the TNFRSF11A gene (c.11G > T, c.31G > A, c.37C > G, c.514G > A), associations with the observed bone parameters (bone strength, bone dimensions and bone mass) were demonstrated. The genetic architecture of bone traits is not fully understood, therefore the present study and the knowledge gained can help to increase the potential in poultry breeding processes and thus reduce the death of individuals.

RNA-seq analysis of the active chick embryo chorioallantoic membrane reveals genes that encode proteins assigned to ion transport and innate immunity

The chicken chorioallantoic membrane (CAM) is an extraembryonic membrane that is vital for the embryo. It undergoes profound cell differentiation between 11 and 15 days of embryonic incubation (EID), which corresponds to the acquisition of its physiological functions. To gain insight into the functional genes that accompany these biological changes, RNA sequencing of the CAM at EID11 and EID15 was performed. Results showed that CAM maturation coincides with the overexpression of 4225 genes, including many genes encoding proteins involved in mineral metabolism, innate immunity, homeostasis, angiogenesis, reproduction, and regulation of hypoxia. Of these genes, some exhibit variability in expression depending on the chicken breed (broiler versus layer breeds). Besides the interest of these results for the poultry sector, the identification of new functional gene candidates opens additional research avenues in the field of developmental biology.

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 Association Study Revealed the Effect of rs312715211 in ZNF652 Gene on Abdominal Fat Percentage of Chickens

Abdominal fat percentage (AFP) is an important economic trait in chickens. Intensive growth selection has led to the over-deposition of abdominal fat in chickens, but the genetic basis of AFP is not yet clear. Using 520 female individuals from selection and control lines of Jingxing yellow chicken, we investigated the genetic basis of AFP using a genome-wide association study (GWAS) and fixation indices (F_ST). A 0.15 MB region associated with AFP was located on chromosome 27 and included nine significant single nucleotide polymorphisms (SNPs), which could account for 3.34-5.58% of the phenotypic variation. In addition, the π value, genotype frequency, and dual-luciferase results identified SNP rs312715211 in the intron region of ZNF652 as the key variant. The wild genotype was associated with lower AFP and abdominal fat weight (AFW), but higher body weight (BW). Finally, annotated genes based on the top 1% SNPs were used to investigate the physiological function of ZNF652. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis suggested that ZNF652 may reduce AFW and BW in broilers through the TGF-β1/SMad2/3 and MAPK/FoxO pathways via EGFR and TGFB1. Our findings elucidated the genetic basis of chicken AFP, rs312715211 on the ZNF652 gene, which can affect BW and AFW and was the key variant associated with AFP. These data provide new insight into the genetic mechanism underlying AF deposition in chickens and could be beneficial in breeding chickens for AF.

Genomic signatures reveal selection in Lingxian white goose

Lingxian white goose (LXW) is a goose breed indigenous to China that is famous for its meat quality and fast growth. However, the genomic evidence underlying such excellent breeding characteristics remains poorly understood. Therefore, we performed whole-genome resequencing of 141 geese from 3 indigenous breeds to scan for selection signatures and detect genomic regions related to breed features of LXW. We identified 5 reproduction-related genes (SYNE1, ESR1, NRIP1, CCDC170, and ARMT1) in highly differentiated regions and 11 notable genes in 26 overlapping windows, some of which are responsible for meat quality (DHX15), growth traits (LDB2, SLIT2, and RBPJ), reproduction (KCNIP4), and unique immunity traits (DHX15 and SLIT2). These findings provide insights into the genetic characteristics of LXW and identify genes affecting important traits in LXW, which extends the genetic resources and basis for facilitating genetic improvement in domestic geese breeds.

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.

Genetic Parameter Estimation and Whole Sequencing Analysis of the Genetic Architecture of Chicken Keel Bending

Bone health is particularly important for high-yielding commercial layer chickens. The keel of poultry is an extension of the abdomen side of the sternum along the sagittal plane and is one of the most important bones. In this study, the keel phenotype of White Leghorns laying hen flocks showed significant individual differences. To clarify its genetic mechanism, we first estimated the heritability of keel bend (KB) in White Leghorn, recorded the production performance of the chicken flock, examined the blood biochemical indexes and bone quality in KB and keel normal (KN) chickens, and performed whole-genome pooled sequencing in KB and KN chickens. We then performed selection elimination analysis to determine the genomic regions that may affect the keel phenotypes. The results show that KB is a medium heritability trait. We found that cage height had a significant effect on the KB (p < 0.01). At 48 weeks, there were significant differences in the number of eggs, the number of normal eggs, and eggshell strength (p < 0.05). The content of parathyroid hormone was lower (p < 0.01) and that of calcitonin was higher (p < 0.01) in KB chickens than in KN chickens. The differences in bone mineral density, bone strength, and bone cortical thickness of the humerus and femur were extremely significant (p < 0.01), with all being lower in KB chickens than in KN chickens. In addition, the bones of KB chickens contained more fat organization. A total of 128 genes were identified in selective sweep regions. We identified 10 important candidate genes: ACP5, WNT1, NFIX, CNN1, CALR, FKBP11, TRAPPC5, MAP2K7, RELA, and ENSGALG00000047166. Among the significantly enriched Kyoto Encyclopedia of Genes and Genomes pathways found, we identifed two bone-related pathways, one involving “osteoclast differentiation” and the other the “MAPK signaling pathway.” These results may help us better understand the molecular mechanism of bone traits in chickens and other birds and provide new insights for the genetic breeding of chickens.

Genome-Wide Association Studies Provide Insight Into the Genetic Determination for Hyperpigmentation of the Visceral Peritoneum in Broilers

Hyperpigmentation of the visceral peritoneum (HVP) has been becoming one of the most challenging problems in yellow-feathered chicken production, which seriously affected chicken carcass quality traits. Detecting which genes dominantly impact pigmentation in the peritoneum tissues is of great benefit to the genetic improvement of HVP. To investigate the genetic mechanism of HVP in yellow-feathered broilers, genome-wide association studies (GWASs) were conducted in the F₂ generation of a cross broiler population with 395 birds. A total of 115,706 single-nucleotide polymorphisms (SNPs) of 122,415 were retained to identify quantitative trait loci (QTL) associated to HVP in chicken. The GWAS results based on the logistic mixed model (LMM) revealed that a narrow genomic location on chromosomes 1 (49.2–51.3 Mb) was significantly associated (p ≤ 4.32 × 10⁻⁷) with HVP, which contained 23 SNP makers related to 14 functional genes (MFNG, POLDIP3, POLR2F, PICK1, PDXP, SGSM3, RANGAP1, MYH9, RPL3, GALP3, LGALS1, MICALL1, ATF4, and CYP2D6). Four highly associated (p < 10⁻⁵) haplotype blocks of 0.80 kb (two SNPs), 0.06 kb (two SNPs), 0.95 kb (two SNPs), and 0.03 kb (two SNPs) were identified with two, two, four, and four distinct haplotypes, respectively. As a melanoma-associated gene, CYP2D6 were also possibly involved in the development of HVP occurring in chicken with two significant variations (rs314284996 and rs317955795) in the promoter regions. Further tests revealed that the expression of CYP2D6 was obviously higher in the visceral peritoneum tissue of chicken with HVP than that in the normal group (p < 0.05). Our results provide a novel clue to understand the genetic mechanism of HVP generation in chicken, and the mapped QTL or candidate genes might serve for genomic selection to improve carcass quality in the yellow-feathered chicken industry.