Genome-Wide Linkage Analysis Identifies Loci for Physical Appearance Traits in Chickens

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.

Genetic markers associated with bone composition in Rhode Island Red laying hens

BACKGROUND

Bone damage has welfare and economic impacts on modern commercial poultry and is known as one of the major challenges in the poultry industry. Bone damage is particularly common in laying hens and is probably due to the physiological link between bone and the egg laying process. Previous studies identified and validated quantitative trait loci (QTL) for bone strength in White Leghorn laying hens based on several measurements, including bone composition measurements on the cortex and medulla of the tibia bone. In a previous pedigree-based analysis, bone composition measurements showed heritabilities ranging from 0.18 to 0.41 and moderate to strong genetic correlations with tibia strength and density. Bone composition was measured using infrared spectroscopy and thermogravimetry. The aim of this study was to combine these bone composition measurements with genotyping data via a genome-wide association study (GWAS) to investigate genetic markers that contribute to genetic variance in bone composition in Rhode Island Red laying hens. In addition, we investigated the genetic correlations between bone composition and bone strength.

RESULTS

We found novel genetic markers that are significantly associated with cortical lipid, cortical mineral scattering, medullary organic matter, and medullary mineralization. Composition of the bone organic matter showed more significant associations than bone mineral composition. We also found interesting overlaps between the GWAS results for tibia composition traits, particularly for cortical lipid and tibia strength. Bone composition measurements by infrared spectroscopy showed more significant associations than thermogravimetry measurements. Based on the results of infrared spectroscopy, cortical lipid showed the highest genetic correlations with tibia density, which was negative (- 0.20 ± 0.04), followed by cortical CO3/PO4 (0.18 ± 0.04). Based on the results of thermogravimetry, medullary organic matter% and mineral% showed the highest genetic correlations with tibia density (- 0.25 ± 0.04 and 0.25 ± 0.04, respectively).

CONCLUSIONS

This study detected novel genetic associations for bone composition traits, particularly those involving organic matter, that could be used as a basis for further molecular genetic investigations. Tibia cortical lipids displayed the strongest genetic associations of all the composition measurements, including a significantly high genetic correlation with tibia density and strength. Our results also highlighted that cortical lipid may be a key measurement for further avian bone studies.

Genomic and gene expression associations to morphology of a sexual ornament in the chicken

How sexual selection affects the genome ultimately relies on the strength and type of selection, and the genetic architecture of the involved traits. While associating genotype with phenotype often utilizes standard trait morphology, trait representations in morphospace using geometric morphometric approaches receive less focus in this regard. Here, we identify genetic associations to a sexual ornament, the comb, in the chicken system (Gallus gallus). Our approach combined genome-wide genotype and gene expression data (>30k genes) with different aspects of comb morphology in an advanced intercross line (F8) generated by crossing a wild-type Red Junglefowl with a domestic breed of chicken (White Leghorn). In total, 10 quantitative trait loci were found associated to various aspects of comb shape and size, while 1,184 expression QTL were found associated to gene expression patterns, among which 98 had overlapping confidence intervals with those of quantitative trait loci. Our results highlight both known genomic regions confirming previous records of a large effect quantitative trait loci associated to comb size, and novel quantitative trait loci associated to comb shape. Genes were considered candidates affecting comb morphology if they were found within both confidence intervals of the underlying quantitative trait loci and eQTL. Overlaps between quantitative trait loci and genome-wide selective sweeps identified in a previous study revealed that only loci associated to comb size may be experiencing on-going selection under domestication.

Genetic basis of chicken plumage color in artificial population of complex epistasis

Chicken plumage color, the genetic basis of which is often affected by epistasis, has long interested scientists. In the current study, a population of complex epistasis was constructed by crossing dominant White Leghorn chickens with recessive white feather chickens. Through a genome-wide association study, we identified single nucleotide polymorphisms and genes significantly associated with white and colored plumage in hens at different developmental stages. Interestingly, white plumage in adulthood was associated with the recessive white feather gene (TYR), whereas white feathers at birth stage were associated with the dominant white feather gene (PMEL), indicating age-related roles for these genes. TYR was shown to exert an epistatic effect on PMEL in adult hens. Additionally, TYR had an epistatic effect on barred plumage, while barred plumage had an epistatic effect on black plumage. TYR had no epistatic effect on the yellow plumage. We confirmed that the barred plumage gene is CDKN2A, as reported in previous studies. Golgb1 and REEP3, which play important roles in the Golgi network and affect the formation of feather pigments, are important candidate genes for yellow plumage. The candidate genes for black plumage are CAMKK1 and IFT22. Further research is warranted to elucidate the molecular mechanisms underlying these traits.

RNA-sequence reveals differentially expressed genes affecting the crested trait of Wumeng crested chicken

Wumeng crested chicken has a cluster of slender feathers on its head, and the underlying skull region exhibits an obvious tumor-like protrusion. This is the typical skull structure of crested chickens. The associated regulatory genes are located on autosomes and are incompletely dominant. This trait is related to brain herniation, but the genetic mechanisms of its formation and development are unclear. In this study, RNA sequencing (RNA-Seq) analysis was conducted on 6 skull tissue samples from 3 Wumeng crested chickens with prominent skull protrusions and 3 without a prominent skull protrusion phenotype. A total of 46,376,934 to 43,729,046 clean reads were obtained, the percentage of uniquely mapped reads compared with the reference genome was between 89.73%-91.00%, and 39,795,458-41,836,502 unique reads were obtained. Among different genomic regions, the highest frequency of sequencing reads occurred in exon regions (85.44-88.28%). Additionally, a total of 423 new transcripts and 26,999 alternative splicings (AS) events were discovered in this sequencing analysis. This study identified 1,089 differentially expressed genes (DEGs), among which 485 were upregulated and 604 were downregulated. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses indicated that the DEGs were enriched in terms related to signal transduction, cell development, cell differentiation, the lysosome, serine, and threonine metabolism, and the interaction of cytokines with cytokine receptors. Based on the comprehensive analysis of DEGs combined with reported quantitative trait loci (QTLs), the expression of BMP2, EPHA3, EPHB1, HOXC6, SCN2B, BMP7, and HOXC10 was verified by real-time quantitative polymerase chain reaction (qRT-PCR). The qRT-PCR results were consistent with the RNA-Seq results, indicating that these 7 genes may be candidates genes regulating the crested trait.

Genome-wide detection of CNVs and their association with performance traits in broilers

Background

Copy number variations (CNVs) are a major type of structural genomic variants that underlie genetic architecture and phenotypic variation of complex traits, not only in humans, but also in livestock animals. We identified CNVs along the chicken genome and analyzed their association with performance traits. Genome-wide CNVs were inferred from Affymetrix® high density SNP-chip data for a broiler population. CNVs were concatenated into segments and association analyses were performed with linear mixed models considering a genomic relationship matrix, for birth weight, body weight at 21, 35, 41 and 42 days, feed intake from 35 to 41 days, feed conversion ratio from 35 to 41 days and, body weight gain from 35 to 41 days of age.

Results

We identified 23,214 autosomal CNVs, merged into 5042 distinct CNV regions (CNVRs), covering 12.84% of the chicken autosomal genome. One significant CNV segment was associated with BWG on GGA3 (q-value = 0.00443); one significant CNV segment was associated with BW35 (q-value = 0.00571), BW41 (q-value = 0.00180) and BW42 (q-value = 0.00130) on GGA3, and one significant CNV segment was associated with BW on GGA5 (q-value = 0.00432). All significant CNV segments were verified by qPCR, and a validation rate of 92.59% was observed. These CNV segments are located nearby genes, such as KCNJ11, MyoD1 and SOX6, known to underlie growth and development. Moreover, gene-set analyses revealed terms linked with muscle physiology, cellular processes regulation and potassium channels.

Conclusions

Overall, this CNV-based GWAS study unravels potential candidate genes that may regulate performance traits in chickens. Our findings provide a foundation for future functional studies on the role of specific genes in regulating performance in chickens.

Supplementary Information

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

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

Feathered leg is a trait in domestic chickens that has undergone intense selection by fancy breeders. Previous studies have shown that two major loci controlling feathered leg are located on chromosomes 13 and 15. Here, we present genetic evidence for the identification of candidate causal mutations at these loci. This was accomplished by combining classical linkage mapping using an experimental cross segregating for feathered leg and high-resolution identical-by-descent mapping using whole-genome sequence data from 167 samples of chicken with or without feathered legs. The first predicted causal mutation is a single-base change located 25 kb upstream of the gene for the forelimb-specific transcription factor TBX5 on chromosome 15. The second is a 17.7-kb deletion located ∼200 kb upstream of the gene for the hindlimb-specific transcription factor PITX1 on chromosome 13. These mutations are predicted to activate TBX5 and repress PITX1 expression, respectively. The study reveals a remarkable convergence in the evolution of the feathered-leg phenotype in domestic chickens and domestic pigeons, as this phenotype is caused by noncoding mutations upstream of the same two genes. Furthermore, the PITX1 causal variants are large overlapping deletions, 17.7 kb in chicken and 44 kb in pigeons. The results of the present study are consistent with the previously proposed model for pigeon that feathered leg is caused by reduced PITX1 expression and ectopic expression of TBX5 in hindlimb buds resulting in a shift of limb identity from hindlimb to more forelimb-like identity.

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

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

Copy number variation in HOXB7 and HOXB8 involves in the formation of beard trait in chickens

The derived feathering phenotype beard in domestic birds is an ideal resource to investigate the genetic mechanisms controlling feather development and differentiation. In the present study, we performed a GWAS and QTL linkage analysis on the trait of beard in Beijing fatty chicken. One major QTL (1.2-1.9 Mb) was identified that could explain 34% of the phenotypic variation. The copy number variation that was copied from the region (GGA27:3 578 409-3 592 890 bp) containing homebox B7 (HOXB7) and homebox B8 (HOXB8) was validated to be only exhibited in the genome of bearded chickens. Protein-protein interaction analysis indicated that HOXB7 and HOXB8 proteins could highly interact with the HOXB family members, including HOXB4, HOXB5 and HOXB6, whose genomic locations near HOXB7 and HOXB8 suggested that they may regulate their family members to involve in the formation of the beard trait in chickens. Overall, our work provides basic data for understanding the mechanisms regulating beard development and differentiation.

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.

Study on CHADL as a candidate gene for comb growth traits in Partridge Shank roosters

The comb is an important secondary sexual characteristic and comb growth traits, such as size and color of the comb, are widely used as indicators in chicken breeding programs. However, the genetic basis for these traits remains mostly unknown. It was found that the chondroadherin-like (CHADL) gene was up-regulated in large combs and was located in reported comb growth quantitative trait loci. In this study, tissue-specific expressions, expression patterns in combs of different ages, and CHADL polymorphisms were analyzed to investigate the relationship between this gene and comb growth traits of Partridge Shank roosters. The results showed that CHADL was more highly expressed in combs than in 10 other tissues, and its expressions in combs tended to gradually increase from the 5-wk-old mark to the 26-wk-old mark. The single-nucleotide polymorphism rs316423539 in the CHADL gene was significantly associated with the comb area and height, whereas rs14822286 was highly correlated with the comb color. Moreover, H1H5, H1H6, and H3H6 were the most advantageous genotype combinations for comb growth traits. Our results might help understand the molecular mechanism of comb growth traits and improve these traits directly by marker assistant selections.

Efficiency of using SNP markers in the <i>MSTN</i> gene in the selection of the Pushkin breed chickens

In the poultry industry, indicators reflecting the growth rate of young stock and the exterior characteristics of chickens are important benchmarks for breeding. Traditional selection based on phenotypic evaluation is characterized by low efficiency with a low character inheritance ratio and is difficult to apply in small groups of animals and birds bred in bioresource collections. The use of molecular genetic markers associated with economically important traits makes it possible to carry out early selection of birds. This entails an increase in the profitability of the poultry industry. Recently, single nucleotide polymorphisms (SNPs) have served as convenient markers for selection purposes. For five generations (P1–P5), an experimental selection of hens of the Pushkin breed was carried out for live weight. It was based on selection for single nucleotide polymorphism rs313744840 in the MSTN gene. As a result, a significant increase in the frequency of allele A in this gene, from 0.11 to 0.50, took place. The association of SNP markers with meat qualities in the experimental group led to changes in the exterior profile of an adult bird at 330 days of age. The individuals with the AA and AG genotypes had the greatest live weight and longest body. As a result of selection, the bird on average became larger due to an increase in the number of heterozygous individuals with long bodies and large chest girths. The depth of the chest and the width of the pelvis increased due to an increase in the frequency of allele A in the experimental population. A tendency towards an increase in these indicators with the substitution of G with A in the genotype was found. Saturation of the population with desirable alleles led to an increase in the average live weight of the chickens. Analysis of the exterior parameters of adult birds showed that this growth is achieved by increasing the depth and volume of the bird body, and not by increasing the length of the limbs. Thus, marker selection carried out for five generations in the experimental population of Pushkin breed chickens to increase body weight has reliably (p < 0.001) changed the exterior profile of adult birds.

Genome-wide association studies targeting the yield of extraembryonic fluid and production traits in Russian White chickens

Background

The Russian White is a gene pool breed, registered in 1953 after crossing White Leghorns with local populations and, for 50 years, selected for cold tolerance and high egg production (EL). The breed has great potential in meeting demands of local food producers, commercial farmers and biotechnology sector of specific pathogen-free (SPF) eggs, the former valuing the breed for its egg weight (EW), EL, age at first egg (AFE), body weight (BW), and the latter for its yield of extraembryonic fluid (YEF) in 12.5-day embryos, ratio of extraembryonic fluid to egg weight, and embryo mass. Moreover, its cold tolerance has been presumably associated with day-old chick down colour (DOCDC) – white rather than yellow, the genetic basis of these traits being however poorly understood.

Results

We undertook genome-wide association studies (GWASs) for eight performance traits using single nucleotide polymorphism (SNP) genotyping of 146 birds and an Illumina 60KBeadChip. Several suggestive associations (p < 5.16*10^− 5) were found for YEF, AFE, BW and EW. Moreover, on chromosome 2, an association with the white DOCDC was found where there is an linkage disequilibrium block of SNPs including genes that are responsible not for colour, but for immune resistance.

Conclusions

The obtained GWAS data can be used to explore the genetics of immunity and carry out selection for increasing YEF for SPF eggs production.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5605-5) contains supplementary material, which is available to authorized users.

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

Chicken body weight (BW) is an economically important trait, and many studies have been conducted on genetic selection for BW. However, previous studies have detected functional chromosome mutations or regions using gene chips. The present study used the specific-locus amplified fragment sequencing (SLAF-seq) technology to perform a genome-wide association study (GWAS) on purebred Wengshang Barred chicken. A total of 1,286,715 single-nucleotide polymorphisms (SNPs) were detected, and 175,211 SNPs were selected as candidate SNPs for genome-wide association analysis using TASSEL general linear models. Six SNP markers reached genome-wide significance. Of these, rs732048524, rs735522839, rs738991545, and rs15837818 were significantly associated with body weight at 28 days (BW28), while rs314086457 and rs315694878 were significantly associated with BW120. These SNPs are close to seven genes (PRSS23, ME3, FAM181B, NABP1, SDPR, TSSK6L2, and RBBP8). Moreover, 24 BW-associated SNPs reached "suggestive" genome-wide significance. Of these, 6, 13, 1, and 4 SNPs were associated with BW28, BW56, BW80, and BW120, respectively. These results would enrich the studies on BW and promote the use of Chinese chicken, especially the Wenshang Barred chicken.

Genomic determinants of epidermal appendage patterning and structure in domestic birds

Variation in regional identity, patterning, and structure of epidermal appendages contributes to skin diversity among many vertebrate groups, and is perhaps most striking in birds. In pioneering work on epidermal appendage patterning, John Saunders and his contemporaries took advantage of epidermal appendage diversity within and among domestic chicken breeds to establish the importance of mesoderm-ectoderm signaling in determining skin patterning. Diversity in chickens and other domestic birds, including pigeons, is driving a new wave of research to dissect the molecular genetic basis of epidermal appendage patterning. Domestic birds are not only outstanding models for embryonic manipulations, as Saunders recognized, but they are also ideal genetic models for discovering the specific genes that control normal development and the mutations that contribute to skin diversity. Here, we review recent genetic and genomic approaches to uncover the basis of epidermal macropatterning, micropatterning, and structural variation. We also present new results that confirm expression changes in two limb identity genes in feather-footed pigeons, a case of variation in appendage structure and identity.

A Complex Structural Variation on Chromosome 27 Leads to the Ectopic Expression of HOXB8 and the Muffs and Beard Phenotype in Chickens

Muffs and beard (Mb) is a phenotype in chickens where groups of elongated feathers gather from both sides of the face (muffs) and below the beak (beard). It is an autosomal, incomplete dominant phenotype encoded by the Muffs and beard (Mb) locus. Here we use genome-wide association (GWA) analysis, linkage analysis, Identity-by-Descent (IBD) mapping, array-CGH, genome re-sequencing and expression analysis to show that the Mb allele causing the Mb phenotype is a derived allele where a complex structural variation (SV) on GGA27 leads to an altered expression of the gene HOXB8. This Mb allele was shown to be completely associated with the Mb phenotype in nine other independent Mb chicken breeds. The Mb allele differs from the wild-type mb allele by three duplications, one in tandem and two that are translocated to that of the tandem repeat around 1.70 Mb on GGA27. The duplications contain total seven annotated genes and their expression was tested during distinct stages of Mb morphogenesis. A continuous high ectopic expression of HOXB8 was found in the facial skin of Mb chickens, strongly suggesting that HOXB8 directs this regional feather-development. In conclusion, our results provide an interesting example of how genomic structural rearrangements alter the regulation of genes leading to novel phenotypes. Further, it again illustrates the value of utilizing derived phenotypes in domestic animals to dissect the genetic basis of developmental traits, herein providing novel insights into the likely role of HOXB8 in feather development and differentiation.

Genetic variation is a key part for the study of evolution, development and differentiation. In domestic animals, many breeds display striking phenotypes that differentiate them from their wild ancestors. Several of these have been related to structural variations, including Fibromelanosis and Rose-comb in chickens, Double-muscled and Osteopetrosis in cattle, Cone degeneration in dogs, and White coat color in pigs. The feather is a type of skin appendages that exists in multiple variants on different body parts, and the derived feathering phenotypes in domestic birds are perfect resources to decipher the mechanisms regulating feather development and differentiation. Here we study the genetics of the Muffs and beard trait, a variant that alters the feather development in the facial area of chickens. We show that this phenotype is associated with a genomic structural variant that leads to an ectopic expression of HOXB8 in the facial skin during feather development. This is thus another example of how structural variants in the genome lead to novel, derived phenotypic changes in domestic animals and suggests an important role for HOXB8 in feather development.