Comparative population genomics reveals genetic basis underlying body size of domestic chickens

bmp10 maintains cardiac function by regulating iron homeostasis

Heart disease remains the leading cause of death worldwide. Iron imbalance, whether deficiency or overload, contributes to heart failure. However, the molecular mechanisms governing iron homeostasis in the heart are poorly understood. Here, we demonstrate that mutation of bmp10, a heart-born morphogen crucial for embryonic heart development, results in severe anemia and cardiac hypertrophy in zebrafish. Initially, bmp10 deficiency causes cardiac iron deficiency, which later progresses to iron overload due to the dysregulated hepcidin/ferroportin axis in cardiac cells, leading to ferroptosis and heart failure. Early iron supplementation in bmp10-/- mutants rescues erythropoiesis, while iron chelation in juvenile fishes significantly alleviates cardiac hypertrophy. We further demonstrate that the interplay between HIF1α-driven hypoxic signaling and the IL6/p-STAT3 inflammatory pathways is critical for regulating cardiac iron metabolism. Our findings reveal BMP10 as a key regulator of iron homeostasis in the vertebrate heart and highlight the potential of targeting the BMP10-hepcidin-iron axis as a therapeutic strategy for iron-related cardiomyopathy.

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.

Genome-wide association studies for economically important traits in mink using copy number variation

Copy number variations (CNVs) are structural variants consisting of duplications and deletions of DNA segments, which are known to play important roles in the genetics of complex traits in livestock species. However, CNV-based genome-wide association studies (GWAS) have remained unexplored in American mink. Therefore, the purpose of the current study was to investigate the association between CNVs and complex traits in American mink. A CNV-based GWAS was performed with the ParseCNV2 software program using deregressed estimated breeding values of 27 traits as pseudophenotypes, categorized into traits of growth and feed efficiency, reproduction, pelt quality, and Aleutian disease tests. The study identified a total of 10,137 CNVs (6968 duplications and 3169 deletions) using the Affymetrix Mink 70K single nucleotide polymorphism (SNP) array in 2986 American mink. The association analyses identified 250 CNV regions (CNVRs) associated with at least one of the studied traits. These CNVRs overlapped with a total of 320 potential candidate genes, and among them, several genes have been known to be related to the traits such as ARID1B, APPL1, TOX, and GPC5 (growth and feed efficiency traits); GRM1, RNASE10, WNT3, WNT3A, and WNT9B (reproduction traits); MYO10, and LIMS1 (pelt quality traits); and IFNGR2, APEX1, UBE3A, and STX11 (Aleutian disease tests). Overall, the results of the study provide potential candidate genes that may regulate economically important traits and therefore may be used as genetic markers in mink genomic breeding programs.

Large-scale genomic and transcriptomic analyses elucidate the genetic basis of high meat yield in chickens

INTRODUCTION

Investigating the genetic markers and genomic signatures related to chicken meat production by combing multi-omics methods could provide new insights into modern chicken breeding technology systems.

OBJECT

Chicken is one of the most efficient and environmentally friendly livestock, especially the fast-growing white-feathered chicken (broiler), which is well known for high meat yield, but the underlying genetic basis is poorly understood.

METHOD

We generated whole-genome resequencing of three purebred broilers (n = 748) and six local breeds/lines (n = 114), and sequencing data of twelve chicken breeds (n = 199) were obtained from the NCBI database. Additionally, transcriptome sequencing of six tissues from two chicken breeds (n = 129) at two developmental stages was performed. A genome-wide association study combined with cis-eQTL mapping and the Mendelian randomization was applied.

RESULT

We identified > 17 million high-quality SNPs, of which 21.74% were newly identified, based on 21 chicken breeds/lines. A total of 163 protein-coding genes underwent positive selection in purebred broilers, and 83 genes were differentially expressed between purebred broilers and local chickens. Notably, muscle development was proven to be the major difference between purebred broilers and local chickens, or ancestors, based on genomic and transcriptomic evidence from multiple tissues and stages. The MYH1 gene family showed the top selection signatures and muscle-specific expression in purebred broilers. Furthermore, we found that the causal gene SOX6 influenced breast muscle yield and also related to myopathy occurrences. A refined haplotype was provided, which had a significant effect on SOX6 expression and phenotypic changes.

CONCLUSION

Our study provides a comprehensive atlas comprising the typical genomic variants and transcriptional characteristics for muscle development and suggests a new regulatory target (SOX6-MYH1s axis) for breast muscle yield and myopathy, which could aid in the development of genome-scale selective breeding aimed at high meat yield in broiler chickens.

Expression of bone morphogenetic protein 10 and its role in biomineralization in Hyriopsis cumingii

Shells and pearls are the products of biomineralization of shellfish after ingesting external mineral ions. Bone morphogenetic proteins (BMPs) play a role in a variety of biological function, and the genes that encode them, are considered important shell-forming genes in mollusks and are associated with shell and pearl formation, embryonic development, and other functions, but bone morphogenetic protein 10 (BMP10) is poorly understood in Hyriopsis cumingii. In this study, we cloned Hc-BMP10 and obtained a 2477 bp full-length sequence encoding 460 amino acids with a conserved TGF-β structural domain. During the embryonic developmental stages, the cleavage stage had the highest expression of Hc-BMP10, followed by juvenile clams; the expression in the mantle gradually decreased with increasing mussel age. A strong signal was detected on epidermal cells on the mantle edge by in situ hybridization. In both the shell notching and inserting operations of the pearl fragment assay, we found that the expression of Hc-BMP10 increased after the above treatments. RNA interference assays showed that the silencing of Hc-BMP10 resulted in a change in the morphology of the prismatic layer and nacreous layer, with the prismatic layer less closely aligned and the disordered aragonite flakes in the nacreous layer. These findings indicate that Hc-BMP10 is involved in the growth and development of H. cumingii, as well as the formation of shells and pearls. Therefore, this study provides some reference for selecting superior species for growth and pearl breeding of H. cumingii at a molecular level and further investigation of the molecular mechanism for biomineralization of Hc-BMP10.

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.

Whole-genome sequencing revealed genetic diversity, structure and patterns of selection in Guizhou indigenous chickens

BACKGROUND

The eight phenotypically distinguishable indigenous chicken breeds in Guizhou province of China are great resources for high-quality development of the poultry industry in China. However, their full value and potential have yet to be understood in depth. To illustrate the genetic diversity, the relationship and population structure, and the genetic variation patterns shaped by selection in Guizhou indigenous chickens, we performed a genome-wide analysis of 240 chickens from 8 phenotypically and geographically representative Guizhou chicken breeds and 60 chickens from 2 commercial chicken breeds (one broiler and one layer), together with 10 red jungle fowls (RJF) genomes available from previous studies.

RESULTS

The results obtained in this present study showed that Guizhou chicken breed populations harbored higher genetic diversity as compared to commercial chicken breeds, however unequal polymorphisms were present within Guizhou indigenous chicken breeds. The results from the population structure analysis markedly reflected the breeding history and the geographical distribution of Guizhou indigenous chickens, whereas, some breeds with complex genetic structure were ungrouped into one cluster. In addition, we confirmed mutual introgression within Guizhou indigenous chicken breeds and from commercial chicken breeds. Furthermore, selective sweep analysis revealed candidate genes which were associated with specific and common phenotypic characteristics evolved rapidly after domestication of Guizhou local chicken breeds and economic traits such as egg production performance, growth performance, and body size.

CONCLUSION

Taken together, the results obtained from the comprehensive analysis of the genetic diversity, genetic relationships and population structures in this study showed that Guizhou indigenous chicken breeds harbor great potential for commercial utilization, however effective conservation measures are currently needed. Additionally, the present study drew a genome-wide selection signature draft for eight Guizhou indigenous chicken breeds and two commercial breeds, as well as established a resource that can be exploited in chicken breeding programs to manipulate the genes associated with desired phenotypes. Therefore, this study will provide an essential genetic basis for further research, conservation, and breeding of Guizhou indigenous chickens.

The idiosyncratic genome of Korean long-tailed chicken as a valuable genetic resource

Today, breeds with ornamental traits such as exceptionally long tail feathers are economically valuable. However, the genetic basis of long-tail feathers is yet to be understood. To provide better understanding of long tail feathers, we sequenced Korean long-tailed chicken (KLC) genomes and compared them with genomes of other chicken breeds. We first analyzed the genome structure of KLC and its genomic relationship with other chickens and observed unique characteristics. Subsequently, we searched for genomic regions under selection. Feather keratin 1-like enriched region and several genes were found to have novel putative functions and effects on the long tail trait in KLC. Our findings support the value of KLC as a unique genetic resource and cast light on the genetic basis of long tail traits in avian species. We expect this novel knowledge to provide new genomic evidence and options for designing and implementing genetic improvements of ornamental chicken productivity through precision crossbreeding aids.

Origins, timing and introgression of domestic geese revealed by whole genome data

BACKGROUND

Geese are among the most important poultry species in the world. The current generally accepted hypothesis is that the European domestic geese originated from greylag geese (Anser anser), and Chinese domestic geese have two origins, most of which originated from swan geese (Anser cygnoides), and the Yili goose originated from greylag geese. To explain the origin and demographic history of geese, we selected 14 goose breeds from Europe and China and wild populations of swan and greylag geese, and whole genome sequencing data were obtained for 74 samples.

RESULTS

Population structure analysis and phylogenetic trees showed that the wild ancestor of Chinese domestic geese, except for Yili, is the swan geese, and the wild ancestor of Chinese Yili and European domestic geese is greylag geese. Analysis of the demographic history suggests that the domestication of Chinese geese occurred ~ 3499 years ago and that of the European geese occurred ~ 7552 years ago. Furthermore, gene flow was observed between domestic geese and their wild ancestors. Analysis of introgression showed that Yili geese had been introgressed by Chinese domestic geese, and the body size of Yili geese may be influenced by introgression events of some growth-related genes, including IGF-1.

CONCLUSIONS

Our study provides evidence for the origin of geese at the genome-wide level and advances the understanding of the history of goose domestication and the traits affected by introgression events.

High-density genotyping reveals candidate genomic regions for chicken body size in breeds of Asian origin

Body size is one of the main selection indices in chicken breeding. Although often investigated, knowledge of the underlying genetic mechanisms is incomplete. The aim of the current study was to identify genomic regions associated with body size differences between Asian Game and Asian Bantam type chickens. In this study, 94 and 107 chickens from 4 Asian Game and 5 Asian Bantam type breeds, respectively, were genotyped using the chicken 580K single nucleotide polymorphism (SNP) array. A genome-wide association study (GWAS) and principal component analyses (PCA) were performed to identify genomic regions associated with body size related-traits such as wing length, shank length, shank thickness, keel length, and body weight. Hierarchical clustering of genotype data showed a clear genetic difference between the investigated Asian Game and Asian Bantam chicken types. GWAS identified 16 genomic regions associated with wing length (2, FDR ≤ 0.018), shank thickness (6, FDR ≤ 0.008), keel length (5, FDR ≤ 0.023), and body weight (3, FDR ≤ 0.041). PCA showed that the first principal component (PC1) separated the 2 chicken types and significantly correlated with the measured body size related-traits (P ≤ 2.24e-40). SNPs contributing significantly to PC1 were subjected to a more detailed investigation. This analysis identified 11 regions potentially associated with differences in body size related-traits. A region on chromosome 4 (GGA4) (17.3-21.3 Mb) was detected in both analyses GWAS and PCA. This region harbors 60 genes. Among them are myotubularin 1 (MTM1) and secreted frizzled-related protein 2 (SFPR2) which can be considered as potential candidate genes for body size related-traits. Our results clearly show that the investigated Asian Game type chicken breeds are genetically different from the Asian Bantam breeds. A region on GGA4 between 17.3 and 21.3 Mb was identified which contributes to the phenotypic difference, though further validation of candidate genes is necessary.

Galbase: a comprehensive repository for integrating chicken multi-omics data

Background

Multi-omics data can provide a stereoscopic view to explore potential causal variations and genes, as well as underlying genetic mechanisms of complex traits. However, for many non-mammalian species, including chickens, these resources are poorly integrated and reused, greatly limiting genetic research and breeding processes of the species.

Results

Here, we constructed Galbase, an easily accessible repository that integrates public chicken multi-omics data from 928 re-sequenced genomes, 429 transcriptomes, 379 epigenomes, 15,275 QTL entries, and 7,526 associations. A total of 21.67 million SNPs, 2.71 million InDels, and 488,583 cis-regulatory elements were included. Galbase allows users to retrieve genomic variations in geographical maps, gene expression profiling in heatmaps, and epigenomic signals in peak patterns. It also provides modules for batch annotation of genes, regions, and loci based on multi-layered omics data. Additionally, a series of convenient tools, including the UCSC Genome Browser, WashU Epigenome Browser, BLAT, BLAST, and LiftOver, were also integrated to facilitate search, visualization, and analysis of sequence features.

Conclusion

Galbase grants new opportunities to research communities to undertake in-depth functional genomic studies on chicken. All features of Galbase make it a useful resource to identify genetic variations responsible for chicken complex traits. Galbase is publicly available at http://animal.nwsuaf.edu.cn/ChickenVar.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08598-2.

Large-scale genomic analysis reveals the genetic cost of chicken domestication

BACKGROUND

Species domestication is generally characterized by the exploitation of high-impact mutations through processes that involve complex shifting demographics of domesticated species. These include not only inbreeding and artificial selection that may lead to the emergence of evolutionary bottlenecks, but also post-divergence gene flow and introgression. Although domestication potentially affects the occurrence of both desired and undesired mutations, the way wild relatives of domesticated species evolve and how expensive the genetic cost underlying domestication is remain poorly understood. Here, we investigated the demographic history and genetic load of chicken domestication.

RESULTS

We analyzed a dataset comprising over 800 whole genomes from both indigenous chickens and wild jungle fowls. We show that despite having a higher genetic diversity than their wild counterparts (average π, 0.00326 vs. 0.00316), the red jungle fowls, the present-day domestic chickens experienced a dramatic population size decline during their early domestication. Our analyses suggest that the concomitant bottleneck induced 2.95% more deleterious mutations across chicken genomes compared with red jungle fowls, supporting the "cost of domestication" hypothesis. Particularly, we find that 62.4% of deleterious SNPs in domestic chickens are maintained in heterozygous states and masked as recessive alleles, challenging the power of modern breeding programs to effectively eliminate these genetic loads. Finally, we suggest that positive selection decreases the incidence but increases the frequency of deleterious SNPs in domestic chicken genomes.

CONCLUSION

This study reveals a new landscape of demographic history and genomic changes associated with chicken domestication and provides insight into the evolutionary genomic profiles of domesticated animals managed under modern human selection.

Integrating Genomic and Transcriptomic Data to Reveal Genetic Mechanisms Underlying Piao Chicken Rumpless Trait

Piao chicken, a rare Chinese native poultry breed, lacks primary tail structures, such as pygostyle, caudal vertebra, uropygial gland, and tail feathers. So far, the molecular mechanisms underlying tail absence in this breed remain unclear. In this study, we comprehensively employed comparative transcriptomic and genomic analyses to unravel potential genetic underpinnings of rumplessness in Piao chicken. Our results reveal many biological factors involved in tail development and several genomic regions under strong positive selection in this breed. These regions contain candidate genes associated with rumplessness, including Irx4, Il18, Hspb2, and Cryab. Retrieval of quantitative trait loci (QTL) and gene functions implies that rumplessness might be consciously or unconsciously selected along with the high-yield traits in Piao chicken. We hypothesize that strong selection pressures on regulatory elements might lead to changes in gene activity in mesenchymal stem cells of the tail bud. The ectopic activity could eventually result in tail truncation by impeding differentiation and proliferation of the stem cells. Our study provides fundamental insights into early initiation and genetic basis of the rumpless phenotype in Piao chicken.

Population analysis of the Korean native duck using whole-genome sequencing data

BACKGROUND

Advances in next-generation sequencing technologies have provided an opportunity to perform population-level comparative genomic analysis to discover unique genomic characteristics of domesticated animals. Duck is one of the most popular domesticated waterfowls, which is economically important as a source of meat, eggs, and feathers. The objective of this study is to perform population and functional analyses of Korean native duck, which has a distinct meat flavor and texture phenotype, using whole-genome sequencing data. To study the distinct genomic features of Korean native duck, we conducted population-level genomic analysis of 20 Korean native ducks together with 15 other duck breeds.

RESULTS

A total of 15.56 million single nucleotide polymorphisms were detected in Korean native duck. Based on the unique existence of non-synonymous single nucleotide polymorphisms in Korean native duck, a total of 103 genes related to the unique genomic characteristics of Korean native duck were identified in comparison with 15 other duck breeds, and their functions were investigated. The nucleotide diversity and population structures among the used duck breeds were then compared, and their phylogenetic relationship was analyzed. Finally, highly differentiated genomic regions among Korean native duck and other duck breeds were identified, and functions of genes in those regions were examined.

CONCLUSIONS

This is the first study to compare the population of Korean native duck with those of other duck breeds by using whole-genome sequencing data. Our findings can be used to expand our knowledge of genomic characteristics of Korean native duck, and broaden our understanding of duck breeds.

863 genomes reveal the origin and domestication of chicken

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

Comparative population genomic analysis uncovers novel genomic footprints and genes associated with small body size in Chinese pony

Background

Body size is considered as one of the most fundamental properties of an organism. Due to intensive breeding and artificial selection throughout the domestication history, horses exhibit striking variations for heights at withers and body sizes. Debao pony (DBP), a famous Chinese horse, is known for its small body size and lives in Guangxi mountains of southern China. In this study, we employed comparative population genomics to study the genetic basis underlying the small body size of DBP breed based on the whole genome sequencing data. To detect genomic signatures of positive selection, we applied three methods based on population comparison, fixation index (F_ST), cross population composite likelihood ratio (XP-CLR) and nucleotide diversity (θπ), and further analyzed the results to find genomic regions under selection for body size-related traits.

Results

A number of protein-coding genes in windows with the top 1% values of F_ST (367 genes), XP-CLR (681 genes), and log₂ (θπ ratio) (332 genes) were identified. The most significant signal of positive selection was mapped to the NELL1 gene, probably underlies the body size and development traits, and may also have been selected for short stature in the DBP population. In addition, some other loci on different chromosomes were identified to be potentially involved in the development of body size.

Conclusions

Results of our study identified some positively selected genes across the horse genome, which are possibly involved in body size traits. These novel candidate genes may be useful targets for clarifying our understanding of the molecular basis of body size and as such they should be of great interest for future research into the genetic architecture of relevant traits in horse breeding program.

A 51 bp indel polymorphism within the PTH1R gene is significantly associated with chicken growth and carcass traits

Parathyroid hormone (PTH) is a crucial regulator of calcium homeostasis and bone remodeling, and the parathyroid hormone 1 receptor (PTH1R) belongs to a class II G-protein-coupled receptor. PTH activates PTH1R, which mediates catabolic and anabolic processes in the skeleton. However, the functional mechanism of PTH1R has not been thoroughly elucidated in organisms. This study identified a 51 bp indel mutation in the first intron of the PTH1R gene and elucidated the effect of this gene mutation on the growth and carcass traits in chickens. The results indicated that the 51 bp indel was significantly associated with subcutaneous fat thickness, abdominal fat weight, body weight and daily gain over 4-8 weeks. Furthermore, we found that PTH1R gene expression was highest in the kidney and liver tissues, and it showed a trend of decreasing in leg and breast muscle tissues at different embryonic stages. In addition, we examined the expression of the three genotypes of the PTH1R gene in the liver, breast muscle and abdominal fat and found that the II genotype was significantly higher than the DD and ID genotypes. In summary, these findings suggest that the PTH1R gene can serve as a potential molecular marker for chicken breeding.

Genome-wide genetic structure and selection signatures for color in 10 traditional Chinese yellow-feathered chicken breeds

BACKGROUND

Yellow-feathered chickens (YFCs) have a long history in China. They are well-known for the nutritional and commercial importance attributable to their yellow color phenotype. Currently, there is a huge paucity in knowledge of the genetic determinants responsible for phenotypic and biochemical properties of these iconic chickens. This study aimed to uncover the genetic structure and the molecular underpinnings of the YFCs trademark coloration.

RESULTS

The whole-genomes of 100 YFCs from 10 major traditional breeds and 10 Huaibei partridge chickens from China were re-sequenced. Comparative population genomics based on autosomal single nucleotide polymorphisms (SNPs) revealed three geographically based clusters among the YFCs. Compared to other Chinese indigenous chicken genomes incorporated from previous studies, a closer genetic proximity within YFC breeds than between YFC breeds and other chicken populations is evident. Through genome-wide scans for selective sweeps, we identified RALY heterogeneous nuclear ribonucleoprotein (RALY), leucine rich repeat containing G protein-coupled receptor 4 (LGR4), solute carrier family 23 member 2 (SLC23A2), and solute carrier family 2 member 14 (SLC2A14), besides the classical beta-carotene dioxygenase 2 (BCDO2), as major candidates pigment determining genes in the YFCs.

CONCLUSION

We provide the first comprehensive genomic data of the YFCs. Our analyses show phylogeographical patterns among the YFCs and potential candidate genes giving rise to the yellow color trait of the YFCs. This study lays the foundation for further research on the genome-phenotype cross-talks that define important poultry traits and for formulating genetic breeding and conservation strategies for the YFCs.

An Evolutionary Genomic Perspective on the Breeding of Dwarf Chickens

The evolutionary history for dwarfism in chickens remains an enigma. Herein, we explore the evolution of the Serama, the smallest breed of chicken. Leveraging comparative population genomics, analyses identify several genes that are potentially associated with the growth and development of bones and muscles. These genes, and in particular both POU1F1 and IGF1, are under strong positive selection. Three allopatric dwarf bantams (Serama, Yuanbao, and Daweishan) with different breeding-histories, form distinct clusters and exhibit unique population structures. Parallel genetic mechanisms underlay their variation in body size. These findings provide insights into the multiple and complex pathways, depending on genomic variation, that chicken can take in response to aviculture selection for dwarfism.

Polygenic Basis and Variable Genetic Architectures Contribute to the Complex Nature of Body Weight -A Genome-Wide Study in Four Chinese Indigenous Chicken Breeds

Body weight (BW) is one of the most important economic traits for animal production and breeding, and it has been studied extensively for its phenotype–genotype associations. While mapping studies have mostly aimed at finding as many loci as possible that contributed to the variation in BW, the role of other factors in its genetic architecture, including their frequencies in the population and their interactions, have been largely overlooked. To comprehensively characterized the genetic architecture of BW, we performed a genome-wide association study (GWAS) both at the single-marker and haplotype level on birds from four indigenous Chinese chicken breeds (Chahua, Silkie, Langshan, and Beard), rather than studying crosses between two founder lines. Additionally, samples from two more breeds (Red Junglefowl and Recessive White) were included to better reflect variable genetic characteristics across populations. Six loci were mapped in this study, revealing the polygenic basis underlying BW. Moreover, by further examining the frequencies of the significantly associated haplotypes in each subpopulation and their effect sizes, most of the loci were found to affect BW in the Beard chicken breed alone. Two loci in GGA9 and GGA27, however, had a common effect on BW across subpopulations, showing that different underlying genetic mechanisms contribute to the phenotypic variability. These findings, particularly the variable genetic architectures found in different loci, improve our understanding of the overall genetic contributions to the large variability in BW among Chinese indigenous chicken breeds. These findings thus will have important implications for future chicken breeding.

Annotating long intergenic non-coding RNAs under artificial selection during chicken domestication

Background

Numerous biological functions of long intergenic non-coding RNAs (lincRNAs) have been identified. However, the contribution of lincRNAs to the domestication process has remained elusive. Following domestication from their wild ancestors, animals display substantial changes in many phenotypic traits. Therefore, it is possible that diverse molecular drivers play important roles in this process.

Results

We analyzed 821 transcriptomes in this study and annotated 4754 lincRNA genes in the chicken genome. Our population genomic analysis indicates that 419 lincRNAs potentially evolved during artificial selection related to the domestication of chicken, while a comparative transcriptomic analysis identified 68 lincRNAs that were differentially expressed under different conditions. We also found 47 lincRNAs linked to special phenotypes.

Conclusions

Our study provides a comprehensive view of the genome-wide landscape of lincRNAs in chicken. This will promote a better understanding of the roles of lincRNAs in domestication, and the genetic mechanisms associated with the artificial selection of domestic animals.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-1036-6) contains supplementary material, which is available to authorized users.