Whole-genome resequencing reveals genetic indels of feathered-leg traits in domestic chickens

Whole genome resequencing reveals the adaptability of native chickens to drought, tropical and frigid environments in Xinjiang

Chickens exhibit extensive genetic diversity and are distributed worldwide. Different chicken breeds have evolved to thrive in diverse environmental conditions. However, research on the genetic mechanisms underlying chicken adaptation to extreme environments, such as tropical, frigid and drought-prone regions, remains limited. In this study, we conducted whole-genome sequencing of 240 individuals from six native chicken breeds in Xinjiang, China, as well as 4 publicly available chicken breeds inhabiting regions with varying annual precipitations, temperatures, and altitudes. Our analysis revealed several genetic variants among the examined breeds. Furthermore, we investigated the genetic diversity and population structure of breeds residing in extreme drought and temperature environments by comparing them. Notably, native chicken breeds exhibited different genetic diversity and population structures. Moreover, we identified candidate genes associated with chicken adaptability to the environment, such as CORO2A, CTNNA3, AGMO, GRID2, BBOX1, COL3A1, INSR, SOX5, MAP2 and PLPPR1. Additionally, pathways such as lysosome, cysteine and methionine metabolism, glycosaminoglycan degradation, and Wnt signaling may be play crucial roles in regulating chicken adaptation to drought environments. Overall, these findings contribute to our understanding of the genetic mechanisms governing chicken adaptation to extreme environments, and also offer insights for enhancing the resilience of chicken breeds to different climatic conditions.

A chromosome-level genome assembly of the Asian house martin implies potential genes associated with the feathered-foot trait

The presence of feathers is a vital characteristic among birds, yet most modern birds had no feather on their feet. The discoveries of feathers on the hind limbs of basal birds and dinosaurs have sparked an interest in the evolutionary origin and genetic mechanism of feathered feet. However, the majority of studies investigating the genes associated with this trait focused on domestic populations. Understanding the genetic mechanism underpinned feathered-foot development in wild birds is still in its infancy. Here, we assembled a chromosome-level genome of the Asian house martin (Delichon dasypus) using the long-read High Fidelity sequencing approach to initiate the search for genes associated with its feathered feet. We employed the whole-genome alignment of D. dasypus with other swallow species to identify high-SNP regions and chromosomal inversions in the D. dasypus genome. After filtering out variations unrelated to D. dasypus evolution, we found six genes related to feather development near the high-SNP regions. We also detected three feather development genes in chromosomal inversions between the Asian house martin and the barn swallow genomes. We discussed their association with the wingless/integrated (WNT), bone morphogenetic protein, and fibroblast growth factor pathways and their potential roles in feathered-foot development. Future studies are encouraged to utilize the D. dasypus genome to explore the evolutionary process of the feathered-foot trait in avian species. This endeavor will shed light on the evolutionary path of feathers in birds.

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.

Integration of Whole-Genome Resequencing and Transcriptome Sequencing Reveals Candidate Genes in High Glossiness of Eggshell

Eggshell gloss is an important characteristic for the manifestation of eggshell appearance. However, no study has yet identified potential candidate genes for eggshell gloss between high-gloss (HG) and low-gloss (LG) chickens. The aim of this study was to perform a preliminary investigation into the formation mechanism of eggshell gloss and to identify potential genes. The eggshell gloss of 300-day-old Rhode Island Red hens was measured from three aspects. Uterine tissues of the selected HG and LG (n = 5) hens were collected for RNA-seq. Blood samples were also collected for whole-genome resequencing (WGRS). RNA-seq analysis showed that 150 differentially expressed genes (DEGs) were identified in the uterine tissues of HG and LG hens. These DEGs were mainly enriched in the calcium signaling pathway and the neuroactive ligand-receptor interaction pathway. Importantly, these two pathways were also significantly enriched in the WGRS analysis results. Further joint analysis of WGRS and RNA-seq data revealed that 5-hydroxytryptamine receptor 1F (HTR1F), zinc finger protein 536 (ZNF536), NEDD8 ubiquitin-like modifier (NEDD8), nerve growth factor (NGF) and calmodulin 1 (CALM1) are potential candidate genes for eggshell gloss. In summary, our research provides a reference for the study of eggshell gloss and lays a foundation for improving egg glossiness in layer breeding.

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.

Whole Genome Resequencing Helps Study Important Traits in Chickens

The emergence of high-throughput sequencing technology promotes life science development, provides technical support to analyze many life mechanisms, and presents new solutions to previously unsolved problems in genomic research. Resequencing technology has been widely used for genome selection and research on chicken population structure, genetic diversity, evolutionary mechanisms, and important economic traits caused by genome sequence differences since the release of chicken genome sequence information. This article elaborates on the factors influencing whole genome resequencing and the differences between these factors and whole genome sequencing. It reviews the important research progress in chicken qualitative traits (e.g., frizzle feather and comb), quantitative traits (e.g., meat quality and growth traits), adaptability, and disease resistance, and provides a theoretical basis to study whole genome resequencing in chickens.

Selective footprints and genes relevant to cold adaptation and other phenotypic traits are unscrambled in the genomes of divergently selected chicken breeds

BACKGROUND

The genomes of worldwide poultry breeds divergently selected for performance and other phenotypic traits may also be affected by, and formed due to, past and current admixture events. Adaptation to diverse environments, including acclimation to harsh climatic conditions, has also left selection footprints in breed genomes.

RESULTS

Using the Chicken 50K_CobbCons SNP chip, we genotyped four divergently selected breeds: two aboriginal, cold tolerant Ushanka and Orloff Mille Fleur, one egg-type Russian White subjected to artificial selection for cold tolerance, and one meat-type White Cornish. Signals of selective sweeps were determined in the studied breeds using three methods: (1) assessment of runs of homozygosity islands, (2) F_ST based population differential analysis, and (3) haplotype differentiation analysis. Genomic regions of true selection signatures were identified by two or more methods or in two or more breeds. In these regions, we detected 540 prioritized candidate genes supplemented them with those that occurred in one breed using one statistic and were suggested in other studies. Amongst them, SOX5, ME3, ZNF536, WWP1, RIPK2, OSGIN2, DECR1, TPO, PPARGC1A, BDNF, MSTN, and beta-keratin genes can be especially mentioned as candidates for cold adaptation. Epigenetic factors may be involved in regulating some of these important genes (e.g., TPO and BDNF).

CONCLUSION

Based on a genome-wide scan, our findings can help dissect the genetic architecture underlying various phenotypic traits in chicken breeds. These include genes representing the sine qua non for adaptation to harsh environments. Cold tolerance in acclimated chicken breeds may be developed following one of few specific gene expression mechanisms or more than one overlapping response known in cold-exposed individuals, and this warrants further investigation.

Identification of Key Candidate Genes for Beak Length Phenotype by Whole-Genome Resequencing in Geese

The domestic goose is an important economic animal in agriculture and its beak, a trait with high heritability, plays an important role in promoting food intake and defending against attacks. In this study, we sequenced 772 420-day-old Xingguo gray geese (XGG) using a low-depth (~1 ×) whole-genome resequencing strategy. We detected 12,490,912 single nucleotide polymorphisms (SNPs) using the standard GATK and imputed with STITCH. We then performed a genome-wide association study on the beak length trait in XGG. The results indicated that 57 SNPs reached genome-wide significance levels for the beak length trait and were assigned to seven genes, including TAPT1, DHX15, CCDC149, LGI2, SEPSECS, ANAPC4, and Slc34a2. The different genotypes of the most significant SNP (top SNP), which was located upstream of LGI2 and explained 7.24% of the phenotypic variation in beak length, showed significant differences in beak length. Priority-based significance analysis concluded that CCDC149, LGI2, and SEPSECS genes in the most significant quantitative trait locus interval were the most plausible positional and functional candidate genes for beak length development in the XGG population. These findings not only enhance our understanding of the genetic mechanism of the beak length phenotype in geese, but also lay the foundation for further studies to facilitate the genetic selection of traits in geese.

Unveiling Comparative Genomic Trajectories of Selection and Key Candidate Genes in Egg-Type Russian White and Meat-Type White Cornish Chickens

Comparison of genomic footprints in chicken breeds with different selection history is a powerful tool in elucidating genomic regions that have been targeted by recent and more ancient selection. In the present work, we aimed at examining and comparing the trajectories of artificial selection in the genomes of the native egg-type Russian White (RW) and meat-type White Cornish (WC) breeds. Combining three different statistics (top 0.1% SNP by F_ST value at pairwise breed comparison, hapFLK analysis, and identification of ROH island shared by more than 50% of individuals), we detected 45 genomic regions under putative selection including 11 selective sweep regions, which were detected by at least two different methods. Four of such regions were breed-specific for each of RW breed (on GGA1, GGA5, GGA8, and GGA9) and WC breed (on GGA1, GGA5, GGA8, and GGA28), while three remaining regions on GGA2 (two sweeps) and GGA3 were common for both breeds. Most of identified genomic regions overlapped with known QTLs and/or candidate genes including those for body temperatures, egg productivity, and feed intake in RW chickens and those for growth, meat and carcass traits, and feed efficiency in WC chickens. These findings were concordant with the breed origin and history of their artificial selection. We determined a set of 188 prioritized candidate genes retrieved from the 11 overlapped regions of putative selection and reviewed their functions relative to phenotypic traits of interest in the two breeds. One of the RW-specific sweep regions harbored the known domestication gene, TSHR. Gene ontology and functional annotation analysis provided additional insight into a functional coherence of genes in the sweep regions. We also showed a greater candidate gene richness on microchromosomes relative to macrochromosomes in these genomic areas. Our results on the selection history of RW and WC chickens and their key candidate genes under selection serve as a profound information for further conservation of their genomic diversity and efficient breeding.

Whole genome re-sequencing identifies unique adaption of single nucleotide polymorphism, insertion/deletion and structure variation related to hypoxia in Tibetan chickens

BACKGROUND

The adaptation to hypoxia in high altitude areas has great research value in the field of biological sciences. Tibetan chicken has unique adaptability to high-altitude, low pressure and anoxic conditions, and served as a biological model to search for genetic diversity of hypoxia adaption.

METHODS

The whole genome re-sequencing technology was conducted to investigate the genetic diversity.

RESULTS

In this study, we obtained quantity genetic resource, contained 5164926 single nucleotide polymorphisms (SNPs), 237504 Insertion/Deletion (InDel), 55606 structural variation types in all chromosomes of Tibetan chicken. Moreover, 17154 non-synonymous mutations, 45763 synonymous mutations, 258 InDel mutations and 9468 structural mutations were detected in coding sequencing (CDS) region. Furthermore, SNPs occur in 591 genes, including HIF1A, VEGF, MAPK 8/9/10/11, PPARA/D/G, NOTCH2, and ABCs, which were involved in 14 hypoxia-related pathways, such as VEGF signaling pathway, MAPK signaling pathway, PPAR signaling pathway and Notch signaling pathway. Among them, 19 genes with non-synonymous SNP variation in CDS were identified. Moreover, structure variation in CDS also occurred in the mentioned above genes with SNPs.

CONCLUSIONS

This study provides useful targets for clarifying the hypoxia adaptability of the domestication of chickens in Tibetan and may help breeding efforts to develop improved breeds for the highlands.

A Novel 65-bp Indel in the GOLGB1 Gene Is Associated with Chicken Growth and Carcass Traits

Simple Summary

Many Chinese-local chickens show slow-growing and low-producing performance, which is not conductive to the development of the poultry industry. The identification of thousands of indels in the last twenty years has helped us to make progress in animal genetics and breeding. Golgin subfamily B member 1 (GOLGB1) is located on chromosome 1 in chickens. Previous study showed that a large number of QTLs on the chicken chromosome 1 were related to the important economic traits. However, the biological function of GOLGB1 gene in chickens is still unclear. In this study, we detected a novel 65-bp indel in the fifth intron of the chicken GOLGB1 gene. Correlation analysis between the 65-bp indel and chicken growth and carcass traits was performed through a yellow chicken population, which is commercial. Results revealed that this 65-bp indel was significantly associated with chicken body weight, highly significantly associated with neck weight, abdominal fat weight, abdominal fat percentage, and the yellow index b of breast. These findings hinted that the 65-bp indel in GOLGB1 could be assigned to a molecular marker in chicken breeding and enhance production in the chicken industry.

Abstract

Golgin subfamily B member 1 (GOLGB1) gene encodes the coat protein 1 vesicle inhibiting factor, giantin. Previous study showed that mutations of the GOLGB1 gene are associated with dozens of human developmental disorders and diseases. However, the biological function of GOLGB1 gene in chicken is still unclear. In this study, we detected a novel 65-bp insertion/deletion (indel) polymorphism in the chicken GOLGB1 intron 5. Association of this indel with chicken growth and carcass traits was analyzed in a yellow chicken population. Results showed that this 65-bp indel was significantly associated with chicken body weight (p < 0.05), highly significantly associated with neck weight, abdominal fat weight, abdominal fat percentage and the yellow index b of breast (p < 0.01). Analysis of genetic parameters indicated that “I” was the predominant allele. Except for the yellow index b of breast, II genotype individuals had the best growth characteristics, by comparison with the ID genotype and DD genotype individuals. Moreover, the mRNA expression of GOLGB1 was detected in the liver tissue of chicken with different GOLGB1 genotypes, where the DD genotype displayed high expression levels. These findings hinted that the 65-bp indel in GOLGB1 could be assigned to a molecular marker in chicken breeding and enhance production in the chicken industry.

Exploration of key regulators driving primary feather follicle induction in goose skin

The primary feather follicles are universal skin appendages widely distributed in the skin of feathered birds. The morphogenesis and development of the primary feather follicles in goose skin remain largely unknown. Here, the induction of primary feather follicles in goose embryonic skin (pre-induction vs induction) was investigated by de novo transcriptome analyses to reveal 409 differentially expressed genes (DEGs). The DEGs were characterized to potentially regulate the de novo formation of feather follicle primordia consisting of placode (4 genes) and dermal condensate (12 genes), and the thickening of epidermis (5 genes) and dermal fibroblasts (17 genes), respectively. Further analyses enriched DEGs into GO terms represented as cell adhesion and KEGG pathways including Wnt and Hedgehog signaling pathways that are highly correlated with cell communication and molecular regulation. Six selected Wnt pathway genes were detected by qPCR with up-regulation in goose skin during the induction of primary feather follicles. The localization of WNT16, SFRP1 and FRZB by in situ hybridization showed weak expression in the primary feather primordia, whereas FZD1, LEF1 and DKK1 were expressed initially in the inter-follicular skin and feather follicle primordia, then mainly restricted in the feather primordia. The spatial-temporal expression patterns indicate that Wnt pathway genes DKK1, FZD1 and LEF1 are the important regulators functioned in the induction of primary feather follicle in goose skin. The dynamic molecular changes and specific gene expression patterns revealed in this report provide the general knowledge of primary feather follicle and skin development in waterfowl, and contribute to further understand the diversity of hair and feather development beyond the mouse and chicken models.