Mutations Upstream of the TBX5 and PITX1 Transcription Factor Genes Are Associated with Feathered Legs in the Domestic Chicken

Genetic diversity, demographic history, and selective signatures of Silkie chicken

BACKGROUND

Silkie is a traditional Chinese chicken breed characterized by its unique combination of specialized morphological traits. While previous studies have focused on the genetic basis of these traits, the overall genomic characteristics of the Silkie breed remain largely unexplored. In this study, we employed whole genome resequencing data to examine the genetic diversity, selective signals and demographic history of the Silkie breed through comparative analyses with seven other Chinese indigenous breeds (IDGBs), a commercial breed, and the wild ancestor Red Jungle Fowl.

RESULTS

In total, 20.8 million high-quality single nucleotide polymorphisms and 86 large structural variations were obtained. We discovered that Silkie exhibits a relatively high level of inbreeding and is genetically distinct from other IDGBs. Furthermore, our analysis indicated that Silkie has experienced a stronger historical population bottleneck and has a smaller effective population size compared with other IDGBs. We identified 45 putatively selected genes that are enriched in the melanogenesis pathway, which probably is related to the feather color. Among these genes, LMBR1 and PDSS2 have been previously associated with the extra toe and the hookless feathers, respectively. Six of the selected genes (KITLG, GSK3B, SOBP, CTBP1, ELMO2, SNRPN) are known to be associated with neurodevelopment and mental diseases in human, and are possibly related to the distinct behavior of Silkie. We further identified structural variants in Silkie and found previously reported variants linked to hyperpigmentation (END3), muff and beard (HOXB8), and Rose-comb phenotype (MNR2). Additionally, we found a 0.61 Mb inversion overlapping with the GMDS gene, which was previously linked to neurodevelopmental defects in zebrafish and humans. This may also be related to the behavior distinctiveness of Silkie.

CONCLUSIONS

Our study revealed that Silkie is genetically distinct and relatively highly inbred compared to other IDGB chicken populations, possibly attributed to more prolong population bottlenecks and selective breeding practice. These results enhance our understanding of how domestication and selective breeding have shaped the genome of Silkie. These findings contribute to the broader field of domestication and avian genomics, and have implications for the future conservation and breeding efforts.

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.

How structural variants shape avian phenotypes: Lessons from model systems

Despite receiving significant recent attention, the relevance of structural variation (SV) in driving phenotypic diversity remains understudied, although recent advances in long-read sequencing, bioinformatics and pangenomic approaches have enhanced SV detection. We review the role of SVs in shaping phenotypes in avian model systems, and identify some general patterns in SV type, length and their associated traits. We found that most of the avian SVs so far identified are short indels in chickens, which are frequently associated with changes in body weight and plumage colouration. Overall, we found that relatively short SVs are more frequently detected, likely due to a combination of their prevalence compared to large SVs, and a detection bias, stemming primarily from the widespread use of short-read sequencing and associated analytical methods. SVs most commonly involve non-coding regions, especially introns, and when patterns of inheritance were reported, SVs associated primarily with dominant discrete traits. We summarise several examples of phenotypic convergence across different species, mediated by different SVs in the same or different genes and different types of changes in the same gene that can lead to various phenotypes. Complex rearrangements and supergenes, which can simultaneously affect and link several genes, tend to have pleiotropic phenotypic effects. Additionally, SVs commonly co-occur with single-nucleotide polymorphisms, highlighting the need to consider all types of genetic changes to understand the basis of phenotypic traits. We end by summarising expectations for when long-read technologies become commonly implemented in non-model birds, likely leading to an increase in SV discovery and characterisation. The growing interest in this subject suggests an increase in our understanding of the phenotypic effects of SVs in upcoming years.

Morphological variation of tail bone among two chicken breeds and their F1 progeny

Fancy breeds of Japanese indigenous chicken display extensive morphological diversity, particularly in tail feathers. Although marked differences in tail and bone traits have been reported between Tosa-jidori (wild type) and Minohikichabo (rich type) breeds, little is known about the pattern of genetic inheritance in cross experiments. Therefore, this study aimed to investigate the strain and sex effects, and inheritance patterns, in the morphometric variation of pygostyle bones among Tosa-jidori, Minohikichabo, and their F1 hybrids. Five morphological traits, angle of the apex of the pygostyle, pygostyle length, margo cranialis length, tail feather number, and body weight, were evaluated at the adult stage. A significant strain difference was detected in all traits, whereas significant sex differences were observed in only three traits, but not in the angle of the apex of the pygostyle and tail feather number. In F1 hybrids, the angle of the apex of the pygostyle was significantly different to that of Tosa-jidori but not that of Minohikichabo, whereas the pygostyle length and tail number of F1 hybrids were significantly different from those of Minohikichabo but not those of Tosa-jidori. A significant heterosis effect was found in the margo cranialis length and body weight. All five traits showed nonadditive inheritance patterns but varied in each trait between partial dominance (angle of the apex of pygostyle), full dominance (pygostyle length and tail feather number), and over-dominance (margo cranialis length and body weight). Interestingly, different patterns of genetic inheritance in the F1 hybrid were observed at different locations, even within the same pygostyle bone. Using the Japanese indigenous chicken model, these results provide a substantial step toward understanding the genetic architecture of morphology in chickens.

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.

Shaping gene expression and its evolution by chromatin architecture and enhancer activity

Transcriptional regulation plays a pivotal role in orchestrating the intricate genetic programs governing embryonic development. The expression of developmental genes relies on the combined activity of several cis-regulatory elements (CREs), such as enhancers and silencers, which can be located at long linear distances from the genes that they regulate and that interact with them through establishment of chromatin loops. Mutations affecting their activity or interaction with their target genes can lead to developmental disorders and are thought to have importantly contributed to the evolution of the animal body plan. The income of next-generation-sequencing approaches has allowed identifying over a million of sequences with putative regulatory potential in the human genome. Characterizing their function and establishing gene-CREs maps is essential to decode the logic governing developmental gene expression and is one of the major challenges of the post-genomic era. Chromatin 3D organization plays an essential role in determining how CREs specifically contact their target genes while avoiding deleterious off-target interactions. Our understanding of these aspects has greatly advanced with the income of chromatin conformation capture techniques and fluorescence microscopy approaches to visualize the organization of DNA elements in the nucleus. Here we will summarize relevant aspects of how the interplay between CRE activity and chromatin 3D organization regulates developmental gene expression and how it relates to pathological conditions and the evolution of animal body plan.

Duck pan-genome reveals two transposon insertions caused bodyweight enlarging and white plumage phenotype formation during evolution

Structural variations (SVs) are a major source of domestication and improvement traits. We present the first duck pan-genome constructed using five genome assemblies capturing ∼40.98 Mb new sequences. This pan-genome together with high-depth sequencing data (∼46.5×) identified 101,041 SVs, of which substantial proportions were derived from transposable element (TE) activity. Many TE-derived SVs anchoring in a gene body or regulatory region are linked to duck's domestication and improvement. By combining quantitative genetics with molecular experiments, we, for the first time, unraveled a 6945 bp Gypsy insertion as a functional mutation of the major gene IGF2BP1 associated with duck bodyweight. This Gypsy insertion, to our knowledge, explains the largest effect on bodyweight among avian species (27.61% of phenotypic variation). In addition, we also examined another 6634 bp Gypsy insertion in MITF intron, which triggers a novel transcript of MITF, thereby contributing to the development of white plumage. Our findings highlight the importance of using a pan-genome as a reference in genomics studies and illuminate the impact of transposons in trait formation and livestock 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.

HOXB8 overexpression induces morphological changes in chicken mandibular skin: an RNA-seq analysis

The Huiyang beard chicken is a well-known Chinese local breed known for its elongated feathers gathered from both sides of the face (muffs) and below the beak (beard), as well as short wattles (SW). The muff and beard (Mb) mutation is caused by ectopic upregulation of the homeobox B8 (HOXB8) gene in the mandibular skin; and the chi-square test showed a significant correlation between SW and Mb genotypes. However, the underlying molecular mechanisms that regulate Mb and SW variations remain unclear. In this study, we investigated the transcriptomes of the mandibular skin and wattles of chickens with and without the Mb genotype to elucidate the molecular basis of these traits. Our results show that HOXB8 is expressed at significantly higher levels in both the mandibular skin and wattles of Mb chickens than in those of wild-type chickens, indicating that HOXB8 regulates both the Mb and SW phenotypes. Key genes for keratin synthesis were highly expressed in the mandibular skin of Mb chickens, suggesting that HOXB8 may play a role in feather development. In wattles, changes in the expression of extracellular matrix synthesis genes may contribute to SW traits. DNA-binding motif analyses revealed that differentially expressed genes were likely to be directly regulated by HOXB8 binding, indicating that HOXB8 may directly or indirectly regulate feather follicle development and wattle growth. Our study identified both known and novel targets, including several genes not previously implicated in feather development and mesenchymal formation. These findings provide insights into the molecular mechanisms of skin appendage variation in birds and offer potential applications in breeding poultry breeds with unique phenotypes.

Identification of genes involved in regulating the development of feathered feet in chicken embryo

The genetic and developmental factors driving the diverse distribution and morphogenesis of feathers and scales on bird feet are yet unclear. Within a single species, Guangxi domestic chickens exhibit dramatic variety in feathered feet, making them an accessible model for research into the molecular basis of variations in skin appendages. In this study, we used H&E staining to observe the morphogenesis of feathered feet, scaled feet and wings skin at different embryonic stages in Longsheng-Feng chickens and Guangxi Partridge chickens. We selected 4 periods (E6, E7, E8, and E12) that play an important role in feather development and performed transcriptome sequencing to screen for candidate genes associated with feathered feet. Through comparison and analysis of transcriptome data, we identified a set of differently expressed genes (DGEs), which were enriched in appendage organ development, hindlimb morphogenesis, activation of transcription factor binding, and binding of sequence-specific DNA in the cis-regulatory region. In addition, we identified some feathered feet-related genes by analyzing the classical signaling pathways that regulate feather development. Finally, we identified candidate genes that regulate feathered feet formation, which include TBX5, PITX1, ZIC1, FGF20, WNT11, WNT7A, WNT16, and SHH. Interestingly, we found that TBX5 was significantly overexpressed in the skin of the feathered feet and had the highest expression at E7 (P < 0.01), whereas PITX1 expression was significantly reduced at E7(P < 0.01). It is hypothesized that TBX5 and PITX1 regulate the development of hair follicles through the Wnt/β-catenin signaling pathway at E7. Our results provide a theoretical basis for investigating the molecular regulatory mechanisms underlying the formation of chicken feathered feet.

An atlas of regulatory elements in chicken: A resource for chicken genetics and genomics

A comprehensive characterization of regulatory elements in the chicken genome across tissues will have substantial impacts on both fundamental and applied research. Here, we systematically identified and characterized regulatory elements in the chicken genome by integrating 377 genome-wide sequencing datasets from 23 adult tissues. In total, we annotated 1.57 million regulatory elements, representing 15 distinct chromatin states, and predicted about 1.2 million enhancer-gene pairs and 7662 super-enhancers. This functional annotation of the chicken genome should have wide utility on identifying regulatory elements accounting for gene regulation underlying domestication, selection, and complex trait regulation, which we explored. In short, this comprehensive atlas of regulatory elements provides the scientific community with a valuable resource for chicken genetics and genomics.

Effect of polymorphisms in the 5’-flanking sequence of MC1R on feather color in Taihang chickens

MC1R plays an important role in the regulation of the formation, transfer, and deposition of melanin in animals and is important for determining coat color. Many studies have reported on single nucleotide polymorphisms (SNPs) in the coding sequence of MC1R. However, few studies have investigated the polymorphisms in the 5’-flanking sequence of MC1R. In this study, we sequenced 2000 bp of the 5’-flanking sequence of MC1R in 300 Taihang chickens with brown feathers (MTH) and 300 Taihang chickens with black feathers (HTH). The sequencing results showed that 4 SNPs (MC1R g.18838722 G > C, g.18838624 T > C, g.18838694 G > A, and g.18838624 C > T) were located in the 5’-flanking sequence of MC1R between the MTH and HTH groups. Association analysis showed that there was a significant correlation between the 4 SNPs and feather color in Taihang chickens. The correlation between MC1R g.18838624 T >C and feather color of Taihang chicken was 100%, of which the CC (E1) genotype is MTH and the TT (E2) genotype is HTH. Furthermore, there was a significant correlation between MC1R g.18838624 T > C and egg production at 302 d. E1 (184.14 ± 0.674) was significantly higher than that in E2 (181.75 ± 0.577) (P < 0.05). Luciferase reporter assays were used to detect the transcriptional activity of MC1R with different SNP genotypes. The results showed that the luciferase activity of E2 was significantly higher than that of E1 (P < 0.05). In addition, transcription factor-binding site predictions showed that E2 creates a new binding site for ZEB1. RT‒qPCR results revealed that the expression of MC1R in E2 was significantly lower than that in E1 (P < 0.05), and the expression of ZEB1 in E2 was significantly higher than that in E1 (P < 0.05). Overexpression and shRNA experiments demonstrated that ZEB1 regulates the expression of MC1R in DF1 cells. ZEB1 has a negative regulatory effect on the transcriptional activity of MC1R; it inhibits the expression of MC1R and affects the feather color of Taihang chickens. This study provides new insight into the molecular mechanism of feather color formation and the transcriptional regulation of MC1R in Taihang chickens.

Using comparative genomics to detect mutations regulating plumage variations in graylag (A. anser) and swan geese (A. cygnoides)

Although graylag geese (A. anser) showed similar plumages of white, grey, and white with grey patches compared to those in swan geese (A. cygnoides), it was believed the substantial molecular mechanism for plumage variations were different. To date, studies on genes responsible for diverse plumages among graylag geese were limited and causal mutations remain unknown. In this study, genomes from 57 individuals belonging to six breeds showing different plumages were sequenced at ∼10X depth. Firstly, the allele frequency differences (AFD) of variants on the scaffold394 (NW_013185915.1) between grey and white goose breeds (A. anser) was calculated and a genomic region between 768,290-779,889 bp was detected to carry candidate variants associated with plumages, including one SNP (g. 775,151G > T, ∼18.6 kb upstream of EDNRB2) found to be fixed in white geese. This region was overlapped with the one detected by the haplotype-based sweep analysis, in which significant signals defined a candidate region of 736,610-820,622 bp on the same scaffold. Results from the transcriptomic data showed that expression levels of EDNRB2 and many other melanogenesis-related genes were significantly decreased among white geese compared to that in grey geese, especially at late embryonic stages (>E15). Modifications at transcriptional levels might result in abnormal melanocyte developments and thus the white plumages when they grow up. In addition, a frameshift mutation (C > -) in exon4 of MLANA gene on scaffold176 (NW_013185876.1) was suggested as the causal mutation for sex-linked dilution phenotype in graylag geese although this requires more demonstration experiments. Together with observed white plumages caused by EDNRB2 mutations in coding regions among swan geese and chicken, our study provided new examples to study the parallel evolution.

Transcriptome Analysis Reveals the Differentially Expressed Genes Associated with Growth in Guangxi Partridge Chickens

The Guangxi Partridge chicken is a well-known chicken breed in southern China with good meat quality, which has been bred as a meat breed to satisfy the increased demand of consumers. Compared with line D whose body weight is maintained at the average of the unselected group, the growth rate and weight of the selected chicken group (line S) increased significantly after breeding for four generations. Herein, transcriptome analysis was performed to identify pivotal genes and signal pathways of selective breeding that contributed to potential mechanisms of growth and development under artificial selection pressure. The average body weight of line S chickens was 1.724 kg at 90 d of age, which showed a significant increase at 90 d of age than line D chickens (1.509 kg), although only the internal organ ratios of lung and kidney changed after standardizing by body weight. The myofiber area and myofiber density of thigh muscles were affected by selection to a greater extent than that of breast muscle. We identified 51, 210, 31, 388, and 100 differentially expressed genes (DEGs) in the hypothalamus, pituitary, breast muscle, thigh muscle, and liver between the two lines, respectively. Several key genes were identified in the hypothalamus-pituitary-muscle axis, such as FST, THSB, PTPRJ, CD36, PITX1, PITX2, AMPD1, PRKAB1, PRKAB2, and related genes for muscle development, which were attached to the cytokine–cytokine receptor interaction signaling pathway, the PPAR signaling pathway, and lipid metabolism. However, signaling molecular pathways and the cell community showed that elevated activity in the liver of line S fowl was mainly involved in focal adhesion, ECM-receptor interaction, cell adhesion molecules, and signal transduction. Collectively, muscle development, lipid metabolism, and several signaling pathways played crucial roles in the improving growth performance of Guangxi Partridge chickens under artificial selection for growth rate. These results support further study of the adaptation of birds under selective pressure.

Who are you? A framework to identify and report genetic sample mix-ups

Sample mix-ups occur when samples have accidentally been duplicated, mislabelled or swapped. When samples are subsequently genotyped or sequenced, this can lead to individual IDs being incorrectly linked to genetic data, resulting in incorrect or biased research results, or reduced power to detect true biological patterns. We surveyed the community and found that almost 80% of responding researchers have encountered sample mix-ups. However, many recent studies in the field of molecular ecology do not appear to systematically report individual assignment checks as part of their publications. Although checks may be done, lack of consistent reporting means that it is difficult to assess whether sample mix-ups have occurred or been detected. Here, we present an easy-to-follow sample verification framework that can utilise existing metadata, including species, population structure, sex and pedigree information. We demonstrate its application to a dataset representing individuals of a threatened Aotearoa New Zealand bird species, the hihi, genotyped on a 50K SNP array. We detected numerous incorrect genotype-ID associations when comparing observed and genetic sex or comparing to relationships in a verified microsatellite pedigree. The framework proposed here helped to confirm 488 individuals (39%), correct another 20 bird-genotype links, and detect hundreds of incorrect sample IDs, emphasizing the value of routinely checking genetic and genomic datasets for their accuracy. We therefore promote the implementation and reporting of this simple yet effective sample verification framework as a standardized quality control step for studies in the field of molecular ecology.

A ROR2 coding variant is associated with craniofacial variation in domestic pigeons

Vertebrate craniofacial morphogenesis is a highly orchestrated process that is directed by evolutionarily conserved developmental pathways.1,2 Within species, canalized development typically produces modest morphological variation. However, as a result of millennia of artificial selection, the domestic pigeon displays radical craniofacial variation within a single species. One of the most striking cases of pigeon craniofacial variation is the short-beak phenotype, which has been selected in numerous breeds. Classical genetic experiments suggest that pigeon beak length is regulated by a small number of genetic factors, one of which is sex linked (Ku2 locus).3-5 However, the genetic underpinnings of pigeon craniofacial variation remain unknown. Using geometric morphometrics and quantitative trait locus (QTL) mapping on an F2 intercross between a short-beaked Old German Owl (OGO) and a medium-beaked Racing Homer (RH), we identified a single Z chromosome locus that explains a majority of the variation in beak morphology in the F2 population. Complementary comparative genomic analyses revealed that the same locus is strongly differentiated between breeds with short and medium beaks. Within the Ku2 locus, we identified an amino acid substitution in the non-canonical Wnt receptor ROR2 as a putative regulator of pigeon beak length. The non-canonical Wnt pathway serves critical roles in vertebrate neural crest cell migration and craniofacial morphogenesis.6,7 In humans, ROR2 mutations cause Robinow syndrome, a congenital disorder characterized by skeletal abnormalities, including a widened and shortened facial skeleton.8,9 Our results illustrate how the extraordinary craniofacial variation among pigeons can reveal genetic regulators of vertebrate craniofacial diversity.

Cis-acting mutation affecting GJA5 transcription is underlying the Melanotic within-feather pigmentation pattern in chickens

Melanotic (Ml) is a mutation in chickens that extends black (eumelanin) pigmentation in normally brown or red (pheomelanin) areas, thus affecting multiple within-feather patterns [J. W. Moore, J. R. Smyth Jr, J. Hered. 62, 215-219 (1971)]. In the present study, linkage mapping using a back-cross between Dark Cornish (Ml/Ml) and Partridge Plymouth Rock (ml+/ml+ ) chickens assigned Ml to an 820-kb region on chromosome 1. Identity-by-descent mapping, via whole-genome sequencing and diagnostic tests using a diverse set of chickens, refined the localization to the genomic region harboring GJA5 encoding gap-junction protein 5 (alias connexin 40) previously associated with pigmentation patterns in zebrafish. An insertion/deletion polymorphism located in the vicinity of the GJA5 promoter region was identified as the candidate causal mutation. Four different GJA5 transcripts were found to be expressed in feather follicles and at least two showed differential expression between genotypes. The results showed that Melanotic constitutes a cis-acting regulatory mutation affecting GJA5 expression. A recent study established the melanocortin-1 receptor (MC1R) locus and the interaction between the MC1R receptor and its antagonist agouti-signaling protein as the primary mechanism underlying variation in within-feather pigmentation patterns in chickens. The present study advances understanding the mechanisms underlying variation in plumage color in birds because it demonstrates that the activity of connexin 40/GJA5 can modulate the periodic pigmentation patterns within individual feathers.

A Missense Mutation in the MYBPH Gene Is Associated With Abdominal Fat Traits in Meat-Type Chickens

Chicken is an important source of protein for human nutrition and a model system for growth and developmental biology. Although the genetic architecture of quantitative traits in meat-type chickens has been the subject of ongoing investigation, the identification of mutations associated with carcass traits of economic interest remains challenging. Therefore, our aim was to identify predicted deleterious mutation, which potentially affects protein function, and test if they were associated with carcass traits in chickens. For that, we performed a genome-wide association analysis (GWAS) for breast, thigh and drumstick traits in meat-type chickens and detected 19 unique quantitative trait loci (QTL). We then used: (1) the identified windows; (2) QTL for abdominal fat detected in a previous study with the same population and (3) previously obtained whole genome sequence data, to identify 18 predicted deleterious single nucleotide polymorphisms (SNPs) in those QTL for further association with breast, thigh, drumstick and abdominal fat traits. Using the additive model, a predicted deleterious SNP c.482C > T (SIFT score of 0.4) was associated (p-value < 0.05) with abdominal fat weight and percentage. This SNP is in the second exon of the MYBPH gene, and its allele frequency deviates from Hardy–Weinberg equilibrium. In conclusion, our study provides evidence that the c.482C > T SNP in the MYBPH gene is a putative causal mutation for fat deposition in meat-type chickens.

The chicken pan-genome reveals gene content variation and a promoter region deletion in IGF2BP1 affecting body size

Domestication and breeding have reshaped the genomic architecture of chicken, but the retention and loss of genomic elements during these evolutionary processes remain unclear. We present the first chicken pan-genome constructed using 664 individuals, which identified an additional ∼66.5 Mb sequences that are absent from the reference genome (GRCg6a). The constructed pan-genome encoded 20,491 predicated protein-coding genes, of which higher expression level are observed in conserved genes relative to dispensable genes. Presence/absence variation (PAV) analyses demonstrated that gene PAV in chicken was shaped by selection, genetic drift, and hybridization. PAV-based GWAS identified numerous candidate mutations related to growth, carcass composition, meat quality, or physiological traits. Among them, a deletion in the promoter region of IGF2BP1 affecting chicken body size is reported, which is supported by functional studies and extra samples. This is the first time to report the causal variant of chicken body size QTL located at chromosome 27 which was repeatedly reported. Therefore, the chicken pan-genome is a useful resource for biological discovery and breeding. It improves our understanding of chicken genome diversity and provides materials to unveil the evolution history of chicken domestication.

Parallel Genetic Origin of Foot Feathering in Birds

Understanding the genetic basis of similar phenotypes shared between lineages is a long-lasting research interest. Even though animal evolution offers many examples of parallelism, for many phenotypes little is known about the underlying genes and mutations. We here use a combination of whole-genome sequencing, expression analyses, and comparative genomics to study the parallel genetic origin of ptilopody (Pti) in chicken. Ptilopody (or foot feathering) is a polygenic trait that can be observed in domesticated and wild avian species and is characterized by the partial or complete development of feathers on the ankle and feet. In domesticated birds, ptilopody is easily selected to fixation, though extensive variation in the type and level of feather development is often observed. By means of a genome-wide association analysis, we identified two genomic regions associated with ptilopody. At one of the loci, we identified a 17-kb deletion affecting PITX1 expression, a gene known to encode a transcription regulator of hindlimb identity and development. Similarly to pigeon, at the second loci, we observed ectopic expression of TBX5, a gene involved in forelimb identity and a key determinant of foot feather development. We also observed that the trait evolved only once as foot-feathered birds share the same haplotype upstream TBX5. Our findings indicate that in chicken and pigeon ptilopody is determined by the same set of genes that affect similar molecular pathways. Our study confirms that ptilopody has evolved through parallel evolution in chicken and pigeon.