Genome-wide association study of body weight in Wenshang Barred chicken based on the SLAF-seq technology

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.

Investigating the impact of pigmentation variation of breast muscle on growth traits, melanin deposition, and gene expression in Xuefeng black-bone chickens

The blackness traits, considered an important economic factor in the black-bone chicken industry, still exhibits a common phenomenon of significant difference in blackness of breast muscle. To improve this phenomenon, this study compared growth traits, blackness traits, and transcriptome of breast muscles between the High Blackness Group (H group) and Low Blackness Group (L group) in the Xuefeng black-bone chickens. The results are as follows: 1) There was no significant difference in growth traits between the H group and the L group (P > 0.05). 2) The skin/breast muscle L values in the H group were significantly lower than those in the L group, while the breast muscle melanin content exhibited the opposite trend (P < 0.05). 3) A significant negative correlation was observed between breast muscle melanin content and skin/breast muscle L value (P < 0.05), and skin L value exhibiting a significant positive correlation with breast muscle L value (P < 0.05). 4) The breast muscle transcriptome comparison between the H group and L group revealed 831 and 405 DEGs in female and male chickens, respectively. This included 37 shared DEGs significantly enriched in melanosome, pigment granule, and the melanogenesis pathway. Seven candidate genes (DCT, PMEL, MLANA, TYRP1, OCA2, EDNRB2, and CALML4) may play a crucial role in the melanin production of breast muscle in Xuefeng black-bone chicken. The findings could accelerate the breeding process for achieving desired levels of breast muscle blackness and contribute to the exploration of the mechanisms underlying melanin production in black-bone chickens.

Whole-Genome Resequencing Revealed Selective Signatures for Growth Traits in Hu and Gangba Sheep

A genomic study was conducted to uncover the selection signatures in sheep that show extremely significant differences in growth traits under the same breed, age in months, nutrition level, and management practices. Hu sheep from Gansu Province and Gangba sheep from the Tibet Autonomous Region in China were selected. We collected whole-genome data from 40 sheep individuals (24 Hu sheep and 16 Gangba sheep), through whole-genome sequencing. Selection signals were analyzed using parameters such as FST, π ratio, and Tajima's D. We have identified several candidate genes that have undergone strong selection, particularly those associated with growth traits. Specifically, five growth-related genes were identified in both the Hu sheep group (HDAC1, MYH7B, LCK, ACVR1, GNAI2) and the Gangba sheep group (RBBP8, ACSL3, FBXW11, PLAT, CRB1). Additionally, in a genomic region strongly selected in both the Hu and Gangba sheep groups (Chr 22: 51,425,001-51,500,000), the growth-associated gene CYP2E1 was identified, further highlighting the genetic factors influencing growth characteristics in these breeds. This study analyzes the genetic basis for significant differences in sheep phenotypes, identifies candidate genes related to sheep growth traits, lays the foundation for molecular genetic breeding in sheep, and accelerates the genetic improvement in livestock.

Mitogenome-wide association study on body measurement traits of Wenshang Barred chickens

Mitochondria are best known for synthesizing ATP through the tricarboxylic acid cycle and oxidative phosphorylation. The cytoplasmic mitochondrial DNA (mtDNA) is important for maintaining the function. This study was designed to reveal the effect of mtDNA on chicken body measurement traits (BMTs). A population of 605 Wenshang Barred chickens were recorded BMTs, including body slope length, keel length, chest width, etc. The single-nucleotide polymorphisms (SNPs) of their mitogenomes were detected by PCR amplification and DNA sequencing. Totally 69 mutations in mitogenome were discovered, including 18 in noncoding region and 51 in coding region. By multi-sequence alignment and haplotype construction, the chickens were clustered into eight haplotypes and further three haplogroups. The association between BMTs and mtDNA SNPs, haplotypes and haplogroups were analyzed in the linear model by ASReml, respectively. Among them, the SNP mt11086 T/C in ND3 was found to significantly affect chest dept (p < .05) and was highly conservative by phylogenetic conservation analyses, which reflected the genetic effect on body size and growth of chickens. No significant association between the mitochondrial haplotypes or haplogroups and BMTs was found. The polymorphic site reflecting body size could be put into chicken breeding programs as the genetic marker.

Integrated Analysis of the Effects of Cecal Microbiota and Serum Metabolome on Market Weights of Chinese Native Chickens

The gut microbiota plays an important role in the physiological activities of the host and affects the formation of important economic traits in livestock farming. The effects of cecal microbiota on chicken weights were investigated using the Guizhou yellow chicken as a model. Experimental cohorts from chickens with high- (HC, n = 16) and low-market-weights (LC, n = 16) were collected. Microbial 16S rRNA gene sequencing and non-targeted serum metabolome data were integrated to explore the effect and metabolic mechanism of cecal microbiota on market weight. The genera Lachnoclostridium, Alistipes, Negativibacillus, Sellimonas, and Ruminococcus torques were enriched in the HC group, while Phascolarctobacterium was enriched in the LC group (p < 0.05). Metabolomic analysis determined that pantothenic acid (vitamin B5), luvangetin (2H-1-benzopyran-6-acrylic acid), and menadione (vitamin K3) were significantly higher in HC serum, while beclomethasone dipropionate (a glucocorticoid) and chlorophene (2-benzyl-4-chlorophenol) were present at higher levels in the LC group. The microbes enriched in HC were significantly positively correlated with metabolites, including pantothenic acid and menadione, and negatively correlated with beclomethasone dipropionate and chlorophene. These results indicated that specific cecal bacteria in Guizhou yellow chickens alter the host metabolism and growth performance. This study provides a reference for revealing the mechanism of cecal microbe actions that affect chicken body weight.

Identification and Validation of Marketing Weight-Related SNP Markers Using SLAF Sequencing in Male Yangzhou Geese

Growth performance is a complex economic trait for avian production. The swan goose (Anser cygnoides) has never been exploited genetically like chickens or other waterfowl species such as ducks. Traditional phenotypic selection is still the main method for genetic improvement of geese body weight. In this study, specific locus amplified fragment sequencing (SLAF-seq) with bulked segregant analysis (BSA) was conducted for discovering and genotyping single nucleotide polymorphisms (SNPs) associated with marketing weight trait in male geese. A total of 149,045 SNPs were obtained from 427,093 SLAF tags with an average sequencing depth of 44.97-fold and a Q30 value of 93.26%. After SNPs' filtering, a total of 12,917 SNPs were included in the study. The 31 highest significant SNPs-which had different allelic frequencies-were further validated by individual-based AS-PCR genotyping in two populations. The association between 10 novel SNPs and the marketing weight of male geese was confirmed. The 10 significant SNPs were involved in linear regression model analysis, which confirmed single-SNP associations and revealed three types of SNP networks for marketing weight. The 10 significant SNPs were located within or close to 10 novel genes, which were identified. The qPCR analysis showed significant difference between genotypes of each SNP in seven genes. Developed SLAF-seq and identified genes will enrich growth performance studies, promoting molecular breeding applications to boost the marketing weight of Chinese geese.

Genome-Wide Association Study Identifies 12 Loci Associated with Body Weight at Age 8 Weeks in Korean Native Chickens

Meat from Korean native chickens (KNCs) has high consumer demand; however, slow growth performance and high variation in body weight (BW) of KNCs remain an issue. Genome-wide association study (GWAS) is a powerful method to identify quantitative trait-associated genomic loci. A GWAS, based on a large-scale KNC population, is needed to identify underlying genetic mechanisms related to its growth traits. To identify BW-associated genomic regions, we performed a GWAS using the chicken 60K single nucleotide polymorphism (SNP) panel for 1328 KNCs. BW was measured at 8 weeks of age, from 2018 to 2020. Twelve SNPs were associated with BW at the suggestive significance level (p < 2.95 × 10⁻⁵) and located near or within 11 candidate genes, including WDR37, KCNIP4, SLIT2, PPARGC1A, MYOCD and ADGRA3. Gene set enrichment analysis based on the GWAS results at p < 0.05 (1680 SNPs) showed that 32 Gene Ontology terms and two Kyoto Encyclopedia of Genes and Genomes pathways, including regulation of transcription, motor activity, the mitogen-activated protein kinase signaling pathway, and tight junction, were significantly enriched (p < 0.05) for BW-associated genes. These pathways are involved in cell growth and development, related to BW gain. The identified SNPs are potential biomarkers in KNC breeding.

Genome-wide association analysis reveals that EDNRB2 causes a dose-dependent loss of pigmentation in ducks

Background

Birds have various plumage color patterns, and spot is a common phenotype. Herein, we conducted genome-wide association studies (GWAS) in a population of 225 ducks with different sized black spots to reveal the genetic basis of this phenomenon.

Results

First, we quantified the black spot phenotype within the duck population. The results showed that the uncolored area of the body surface first appeared on the ventral side. With increasing duck age, the area of the black spots was highly conserved across the whole body surface. The GWAS results identified a 198 kb (Chr4: 10,149,651 bp to 10,348,068 bp) genetic region that was significantly associated with the black spot phenotype. The conditional GWAS and linkage disequilibrium (LD) analysis further narrowed the ultimate candidate region to 167 kb (Chr4: 10,180,939 bp to 10,348,068 bp). A key gene regulating melanoblast migration and differentiation, EDNRB2 (Endothelin B receptor-like), was found in the candidate region and having significant mRNA expression level changes in embryonic duck skin tissue with different spot sizes. The significant SNPs (single nucleotide polymorphisms) associated with the EDNRB2 gene were annotated, and two mutations (Chr4: 10,180,939 T > C and Chr4: 10,190,671 A > T) were found to result in the loss of binding sites for two trans-factors, XBP1 and cMYB. The phenotypic effect of these two mutations suggested that they can regulate the size of black spots in a dose-dependent manner, and Chr4: 10,180,939 T > C was the major allele locus.

Conclusions

Our results revealed that EDNRB2 was the gene responsible for the variation in duck body surface spot size. Chr4: 10,180,939 T > C was the major allele that explained 49.5 % (dorsal side) and 32.9 % (ventral side) of the variation in duck body surface spot size, while 32.1 % (dorsal side) and 19.1 % (ventral side) of the variation could be explained by Chr4: 10,190,671 A > T. The trans-factor prediction also suggested that XBP1 and cMYB have the potential to interact with EDNRB2, providing new insights into the mechanism of action of these genes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07719-7.

Identification of candidate genes associated with slaughter traits in F2 chicken population using genome-wide association study

Slaughter traits are crucial economic traits of chickens. We performed a GWAS to discover critical loci and candidate genes for 21 slaughter traits in an F2 chicken population resulting from crossing Luxi gamecocks and recessive white feather broilers. We found some SNPs and genes which were significantly associated with keel length, head length, body slope length, bilateral leg weight without shin, bilateral foot weight, subcutaneous fat thickness, heart weight, muscular stomach weight and glandular stomach weight. This study provides references for further investigation of slaughter traits and molecular breeding in chicken.

Genome-wide association study reveals novel loci associated with feeding behavior in Pekin ducks

Background

Feeding behavior traits are an essential part of livestock production. However, the genetic base of feeding behavior traits remains unclear in Pekin ducks. This study aimed to determine novel loci related to feeding behavior in Pekin ducks.

Results

In this study, the feeding information of 540 Pekin ducks was recorded, and individual genotype was evaluated using genotyping-by-sequencing methods. Genome-wide association analysis (GWAS) was conducted for feeding behavior traits. Overall, thirty significant (P-value < 4.74E-06) SNPs for feeding behavior traits were discovered, and four of them reached the genome-wide significance level (P-value < 2.37E-07). One genome-wide significance locus associated with daily meal times was located in a 122.25 Mb region on chromosome 2, which was within the intron of gene ubiquitin-conjugating enzyme E2 E2 (UBE2E2), and could explain 2.64% of the phenotypic variation. This locus was also significantly associated with meal feed intake, and explained 2.72% of this phenotypic variation.

Conclusions

This study is the first GWAS for feeding behavior traits in ducks. Our results provide a list of candidate genes associated with feeding behavior, and also help to better understand the genetic mechanisms of feeding behavior patterns in ducks.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07668-1.

Genome-wide association study of body size traits in Wenshang Barred chickens based on the specific-locus amplified fragment sequencing technology

Chicken body size (BS) is an economically important trait, which has been assessed in many studies for genetic selection. However, previous reports detected functional chromosome mutations or regions using gene chips. The present study used the specific-locus amplified fragment sequencing (SLAF-seq) technology to perform a genome-wide association study (GWAS) of purebred Wenshang Barred chickens. A total of 250 one-day-old male chickens were assessed in this study. Body size in individual birds was measured at 56 days. SLAF-seq was used to genotype and GWAS analysis was carried out using the general linear model (GLM) of the TASSEL program. A total of 1,286,715 single-nucleotide polymorphisms (SNPs) were detected, of which 175,211 were tested as candidate SNPs for genome-wide association analysis using the TASSEL general linear model. Three SNPs markers reached genome-wide significance. Of these, chrZ:81729634, chrZ:81841715, and chrZ:81954149 at 81,729,634, 81,841,715, and 81,954,149 bp of GGA Z were significantly associated with body diagonal length at 56 days (BDL56); and tibia length at 56 days (TL56). These SNPs were close to three genes, including ZCCHC7, PAX5, and MELK. These results open new horizons for studies on BS and should promote the use of Chinese chickens, especially Wenshang Barred chickens.

Genomic analysis of GBS data reveals genes associated with facial pigmentation in Xinyang blue-shelled layers

Facial pigmentation is an important economic trait of chickens, especially for laying hens, which will affect the carcass appearance of eliminated layers. Therefore, identifying the genomic regions and exploring the function of this region that contributes to understanding the variation of skin color traits is significant for breeding. In the study, 291 pure-line Xinyang blue-shelled laying hens were selected, of which 75 were dark-faced chickens and 216 were white-faced chickens. The population was sequenced and typed by GBS genotyping technology. The obtained high-quality SNPs and pigmentation phenotypes were analyzed by a genome-wide association study (GWAS) and a F ST scan. Based on the two analytical methods, we identified a same genomic region (10.70-11.60 Mb) on chromosome 20 with 68 significant SNPs ( - log⁡ 10 ( P ) > 6 ), mapped to 10 known genes, including NPEPL1, EDN3, GNAS, C20orf85, VAPB, BMP7, TUBB1, ELMO2, DDX27, and NCOA5, which are associated with dermal hyperpigmentation.

Genome-wide association study reveals the genetic determinism of growth traits in a Gushi-Anka F2 chicken population

Chicken growth traits are economically important, but the relevant genetic mechanisms have not yet been elucidated. Herein, we performed a genome-wide association study to identify the variants associated with growth traits. In total, 860 chickens from a Gushi-Anka F₂ resource population were phenotyped for 68 growth and carcass traits, and 768 samples were genotyped based on the genotyping-by-sequencing (GBS) method. Finally, 734 chickens and 321,314 SNPs remained after quality control and removal of the sex chromosomes, and these data were used to carry out a GWAS analysis. A total of 470 significant single-nucleotide polymorphisms (SNPs) for 43 of the 68 traits were detected and mapped on chromosomes (Chr) 1-6, -9, -10, -16, -18, -23, and -27. Of these, the significant SNPs in Chr1, -4, and -27 were found to be associated with more than 10 traits. Multiple traits shared significant SNPs, indicating that the same mutation in the region might have a large effect on multiple growth or carcass traits. Haplotype analysis revealed that SNPs within the candidate region of Chr1 presented a mosaic pattern. The significant SNPs and pathway enrichment analysis revealed that the MLNR, MED4, CAB39L, LDB2, and IGF2BP1 genes could be putative candidate genes for growth and carcass traits. The findings of this study improve our understanding of the genetic mechanisms regulating chicken growth and carcass traits and provide a theoretical basis for chicken breeding programs.

A Genome-Wide Association Study Identifying Genetic Variants Associated with Growth, Carcass and Meat Quality Traits in Rabbits

Simple Summary

Rabbit meat has been widely consumed in China and is considered as an ideal food source due to its high protein, low fat, low cholesterol and low sodium contents. The growth rate, carcass characteristics and meat quality are considered economically important traits in the rabbit industry. Genomic selection (GS) could facilitate genetic selection for important economic traits, however, the lack of molecular markers for these traits limits the application of GS in rabbits. Genome-wide association study (GWAS) has the potential to comprehensively identify trait-associated molecular markers and has been applied in animal and plant research. In this study, GWAS was used to examine growth, carcass and meat quality traits of meat rabbits based on specific-locus amplified fragment sequencing (SLAF-seq) technology to identify significantly associated SNPs and functional genes, to be used as a basis for prompting the application of GS in rabbits.

Abstract

Growth, carcass characteristics and meat quality are the most important traits used in the rabbit industry. Identification of the candidate markers and genes significantly associated with these traits will be beneficial in rabbit breeding. In this study, we enrolled 465 rabbits, including 16 male Californian rabbits and 17 female Kangda5 line rabbits as the parental generation, along with their offspring (232 male and 200 female), in a genome-wide association study (GWAS) based on SLAF-seq technology. Bodyweight at 35, 42, 49, 56, 63 and 70 d was recorded for growth traits; and slaughter liveweight (84 d) and dressing out percentage were measured as carcass traits; and cooking loss and drip loss were measured as meat quality traits. A total of 5,223,720 SLAF markers were obtained by digesting the rabbit genome using RsaI + EcoRV-HF^® restriction enzymes. After quality control, a subset of 317,503 annotated single-nucleotide polymorphisms (SNPs) was retained for subsequent analysis. A total of 28, 81 and 10 SNPs for growth, carcass and meat quality traits, respectively, were identified based on genome-wide significance (p < 3.16 × 10⁻⁷). Additionally, 16, 71 and 9 candidate genes were identified within 100 kb upstream or downstream of these SNPs. Further analysis is required to determine the biological roles of these candidate genes in determining rabbit growth, carcass traits and meat quality.