- Research Article
- Open Access
Early and late feathering in turkey and chicken: same gene but different mutations
Genetics Selection Evolution volume 50, Article number: 7 (2018)
Sex-linked slow (SF) and fast (FF) feathering rates at hatch have been widely used in poultry breeding for autosexing at hatch. In chicken, the sex-linked K (SF) and k+ (FF) alleles are responsible for the feathering rate phenotype. Allele K is dominant and a partial duplication of the prolactin receptor gene has been identified as the causal mutation. Interestingly, some domesticated turkey lines exhibit similar slow- and fast-feathering phenotypes, but the underlying genetic components and causal mutation have never been investigated. In this study, our aim was to investigate the molecular basis of feathering rate at hatch in domestic turkey.
We performed a sequence-based case–control association study and detected a genomic region on chromosome Z, which is statistically associated with rate of feathering at hatch in turkey. We identified a 5-bp frameshift deletion in the prolactin receptor (PRLR) gene that is responsible for slow feathering at hatch. All female cases (SF turkeys) were hemizygous for this deletion, while 188 controls (FF turkeys) were hemizygous or homozygous for the reference allele. This frameshift mutation introduces a premature stop codon and six novel amino acids (AA), which results in a truncated PRLR protein that lacks 98 C-terminal AA.
We present the causal mutation for feathering rate in turkey that causes a partial C-terminal loss of the prolactin receptor, and this truncated PRLR protein is strikingly similar to the protein encoded by the slow feathering K allele in chicken.
Sex identification is an important management factor within many commercial livestock operations. In poultry, sexing can be performed by examining feathering rate, a non-intrusive approach to separate males and females at hatch. In turkey and chicken layer breeds, sexing at hatch is crucial for production. However, for broiler breeds this method is mainly applied at the parent stock level. In chicken, the sex-linked dominant K locus, which is located on the Z-chromosome, is responsible for feather development and is associated with delayed emergence of primary and secondary flight feathers (SF), while the k + allele is associated with fast emergence (FF) of flight feathers . The status at this locus is widely used for autosexing at hatch . Elferink et al.  studied the molecular basis of the K allele and identified a 176-kb tandem duplication, which includes part of the genes prolactin receptor (PRLR) and sperm flagellar 2 (SPEF2) that are associated with the K allele. Moreover, a molecular test was developed to distinguish between homozygous and heterozygous late feathering males . The 176-kb duplication causes a 149-amino-acid (AA) C-terminal loss of the PRLR protein and is most likely the causal mutation for the SF phenotype . PRLR is a receptor of the anterior pituitary hormone prolactin that belongs to the type I cytokine receptor family  and is involved in various physiological processes including many reproductive and developmental processes, such as hair/coat morphology . The PRLR gene is widely expressed in all embryonic and somatic tissues and its expression is higher in SF than in FF chicks .
The domesticated turkey (Meleagris gallopavo), an important agricultural species and the second largest contributor to world poultry production , shows similar SF and FF phenotypes in some commercial lines , which are used for the same selection goal as in chicken, i.e. reliable and easy determination of sex at hatch. The SF phenotype differs between turkey and chicken with SF turkeys generally showing poor feathering even at a later age . Zakrzewska et al.  suggested that the dominant sex-linked inhibited feathering (IF) allele K is responsible for the genetic feathering defect in turkey. Interestingly, expression of this defect ranges from almost complete absence of feathers to full feather covering at a later age (> 4 weeks of age), although until 4 weeks of age no apparent differences between SF birds were observed. Moreover, SF turkeys show inferior reproductive efficiency compared with FF turkeys  and differences in body weight and carcass characteristics . The SF phenotype that is under study here differs from a late feathering phenotype that was described in turkey by Asmundson and Abbott , which consists in poor feathering at physical maturity (> 20 weeks of age). In chicken, the SF phenotype has been associated with the sex-linked allele K, whereas in turkey, the underlying genetic components and causal mutation have never been investigated. In this study, we used whole-genome sequence data that were obtained from either slow- or fast-feathering turkeys to perform a case–control genome-wide association study (GWAS) for feathering rate at hatch and to investigate its relation to the chicken allele K.
Dataset used for sequencing and mapping
We collected blood from 202 animals representing nine commercial turkey lines and that included 12 SF cases and 12 FF cases selected from the same line. For each sample, DNA was extracted and sequenced on the Illumina HiSeq 2000 sequencer, which generated paired-end 101 bp reads. We used the Sickle software to trim sequences , BWA-MEM (version 0.7.15) to map the whole-genome sequencing data to the turkey reference genome (Melgal5) , the Samtools dedup function to remove duplicate reads , the GATK IndelRealigner to perform local realignments of reads around indels  and Qualimap to obtain mapping statistics .
Variant detection and post-processing
We performed population-based variant calling using the Freebayes software with the following settings: (1) min-base-quality 10 (to exclude alleles with support base quality < 10), (2) min-alternate-fraction 0.2 (at least 20% of the reads should support the alternate allele in order to evaluate the position), (3) haplotype-length 0 (to avoid generating haplotypes in VCF), (4) ploidy 2 (assuming diploid organism), and (5) min-alternate-count 2 (to have at least two reads that support the alternate allele in order to evaluate position) . Post-processing was performed using bcftools , and variants that were located within 3 bp of an indel, or with a phred quality score and call rate lower than 20 and 0.7, respectively, were removed. The average call rate was about 0.985, and the average transition/transversion (TS/TV) ratio was 2.62, in line with previous findings in turkey .
PCA analysis was performed using PLINK  on the filtered vcf files and plotted using the default R plotting utilities.
Functional annotation of variants
SnpEffect  was used for variant annotation and the PROVEAN software for variant effect prediction in missense variants. The following variant classes were considered as potential candidate variants: missense, splice acceptor, splice donor, inframe indels, frameshift, stop lost, stop gained, and start lost variants.
Association study and identification of candidate variants
Single locus associations on the genotypes called by freebayes were tested for SNPs and indels in PLINK using permutations to generate uncorrected and corrected p values . p values were generated by applying the Fisher’s exact test and an adaptive Monte Carlo permutation test was performed with 5000 replications. Variants with a P lower than 1e-5 were considered significant. Manhattan plots were generated using qqman R package . We selected all significant protein-altering variants and evaluated their putative effect on the protein based on PROVEAN scores and SnpEffect annotations. Moreover, gene ontology (GO) annotations were obtained from the Uniprot database . Phenotype information on PRLR null-mutant mice was from Craven et al. . The ClustalO alignment software  was used to align chicken and turkey PRLR sequences.
CNV-seq was used to perform CNV analysis using a log2-threshold of 0.6 and a p value threshold of 0.001 . The optimum window size was automatically computed and ranged from 2.5 to 7.1 kb. The FF sample MG-WUR-121 and the SF sample MG-WUR-136 were used as control samples in CNV-seq analysis for analyses of SF and FF data, both exhibiting average to high coverage (see Additional file 1: Table S1). CNV-seq R utilities were used to plot the CNV events.
Case–control sequencing and variant detection
To study the molecular mechanisms that underlie feathering rate at hatch in turkey, we selected 12 animals within each group (SF and FF) from one commercial line for whole-genome re-sequencing (WGS) (All female, [see Additional file 1: Table S1]). Moreover, DNA from 178 FF turkeys from various commercial turkey lines was sequenced for additional control samples. The SF turkeys in the population analysed here have phenotypes that are similar to those described for the dominant sex-linked IF allele K by Zakrzewska et al. . Whole-genome DNA was sequenced and resulted in a total amount of 2.17 Tbp (tera base pairs) from 22.48 × 109 paired-end 101 bp reads. Mapping was performed with BWA—mem (version 0.7.15) to the Meleagris gallopavo build 5 (Melgal5: ) reference genome with an average mappability and coverage of 98.38%, and 10.5×, respectively. We performed population-based variant calling using Freebayes . Next, we filtered out variants with a low-quality (phred quality score < 20) or a call rate lower than 0.7, which resulted in 8,136,213 (post-filtering) variants including 6,595,059 SNPs, and 1,197,170 indels, with an average variant density of 8.4 variants per kb (see Additional file 1: Table S2). We performed PCA analysis on the 24 cases and control animals to assess population stratification; no distinct clustering was observed between the two groups (see Additional file 2: Figure S1).
Functional annotation of variants
We used SnpEff to assign a range of functional classes to the identified variants . The majority of the variants were located in intronic, ncRNA, or intergenic regions (see Additional file 1: Table S3). We identified 231,073 coding (90,370 protein-altering) variants with an overall missense/silent ratio of 0.545, which means that for every two silent mutations (synonymous) one missense mutation is found (see Additional file 1: Table S4).
Genome-wide association study for feathering rate at hatch
The GWAS revealed a significant signal for 134 SNPs on the Z chromosome. None of the detected variants is in perfect LD with the phenotype (see Additional file 1: Table S5). SNPs associated with the SF phenotype are all located on the short arm of the Z chromosome between 7.95 and 9.79 Mb (Fig. 1) and (see Additional file 2: Figure S2). This region contains 55 protein-coding genes including the PRLR and SPEF2 genes associated with the SF phenotype in chicken.
A 5-bp deletion in the PRLR gene is associated with slow feathering rate in turkey
We examined the putative effects of all significant variants associated with slow feathering. In total, we identified eight protein-altering variants (seven SNPs and one indel). The seven identified SNPs cause missense mutations in protein coding genes (Table 1). None of the missense variants were predicted to have a high impact on the corresponding protein (reaching a PROVEAN score < − 2.5). Moreover, none of the missense mutations were fixed within the group of SF turkeys (Table 1), and thus were unlikely to be the causal variant. The identified indel represents a hemizygous 5-bp deletion that is statistically associated with feathering rate and predicted to have a high impact on the PRLR gene by causing a frameshift (Fig. 2). This deletion, which is located within the terminal exon of the PRLR gene, produces a truncated PRLR protein by introducing a premature stop codon and adding six novel C-terminal amino acids (DSITET*, Fig. 2). All SF turkeys were hemizygous for the alternate allele, while ten FF turkey controls and all additional 178 FF turkeys were hemizygous or homozygous for the reference allele (Table 2). In addition, we performed a copy-number variation (CNV) analysis to test whether, as in chicken, a CNV event is associated with feathering rate at hatch. Although one region on chromosome Z between 7.9 and 8.1 Mb harboured copy number variants in various samples, none of them were associated with feathering rate at hatch (see Additional file 3).
Chicken and turkey slow feathering
Turkey and chicken PRLR proteins are very similar (90.24% sequence identity, (see Additional file 2: Figure S3) and both are 831 AA long. However, carriers of the 5-bp frameshift deletion in turkey lack the final 98 AA of the PRLR C-end tail (Fig. 2), whereas carriers of the K allele in chicken lack the terminal 149 AA of the PLRP C-end tail . The prolactin receptor forms a dimer to bind prolactin in the extracellular space on the N-terminal end of the protein . Moreover, PRLR contains two fibronectin type 3 domains (FN3), a WSXWS motif that is important for proper protein folding, and a Box 1 motif that is necessary for Janus kinase (JAK) interaction and activation . However, the function of the affected C-end tail, which is located in the cytoplasm, is mostly unknown, but likely shares similar functional relevance in chicken and turkey.
This study reveals the molecular mechanisms that underlie the rate of feathering at hatch in turkey. The use of NGS data provided a sufficient number of variants to describe the potential causal polymorphism, i.e. a 5-bp deletion within the last exon of the PRLR gene. This mutation is different from that of allele K in chicken, but impacts the same gene  and moreover, in a similar manner, i.e. by loss of a substantial part of the C-end tail (98 AA in turkey; 149 AA in chicken). Unlike SF chickens, SF turkeys are poorly feathered, even at physical maturity . Moreover, feathering of females can be so poor that carriers of this allele are not used commercially. Although strongly associated with SF, the PRLR 5-bp frameshift mutation is not in perfect LD with the phenotype since we observed two FF females that were hemizygous for this deletion. One possibility is that these two animals were mislabelled as FF turkeys, although they are SF turkeys; this is supported by the observation that none of the variants (including non-coding ones) is in complete LD with the phenotype.
The membrane-protein PRLR is a member of the cytokine receptor family that binds the prolactin hormone (PRL) within the extracellular space . This hormone is involved in a diverse range of biological activities including various reproductive and developmental processes, such as hair replacement and follicle development . Null mutant mice show different hair/coat morphologies and advanced hair replacement . Moreover, a frameshift variant, which introduces a premature stop codon in the bovine PRLR receptor and causes the loss of 120 C-terminal AA, is associated with abnormally short and sleek hair coat . Moreover, hair development and feather development are considered to have an evolutionary homologous origin. Thus, these findings support the PRLR gene as a likely candidate for feathering development within both commercial poultry species, chicken and turkey.
Other studies have suggested that feathering rate in chicken is caused by a higher expression of PRLR due to its partial duplication. Carriers of allele K show a 1.78-fold higher expression of PRLR in chicken . In contrast, Zhao et al.  found no difference in PRLR expression between SF and FF chicks, but that the expression of the other gene involved in the duplication, SPEF2, was significantly higher in SF than in FF chicks, which suggested that a mutation in this gene was responsible for the SF phenotype. We believe that the higher expression of SPEF2 in chicken is due to the large duplication that underlies the K allele. The duplication results in two partial PRLR genes (that lack both tails), while the SPEF2 gene remains complete . Incomplete PRLR mRNA could be subject to the nonsense-mediated decay mechanism resulting in a lower abundance of PRLR mRNA compared to SPEF2 mRNA. Thus, based on our findings, we believe that, rather than a higher expression of PRLR, it is the lack of the C-terminal end of the protein that is responsible for the slow feathering rate at hatch in both chicken and turkey. Interestingly, Nakamura et al.  reported that, in a late feathering chicken line, reversion to the fast feathering phenotype occurred in rare instances, but this was not observed in our population.
The PRLR mutations in chicken and turkey are clearly independent, but lead to similar phenotypes, which strongly suggests that they have been favoured by identical breeding goals being applied in these two species. Thus, the SF/FF phenotype shows a pattern that is similar to that observed for a small number of monogenic or oligogenic traits under domestication selection, which show independent mutations in the same genes in specific pathways . Coat colour is one of the most common domestication features, which is regulated by a small number of genes (e.g. KIT, MC1R, and TYR) in many domestic animals . Another example in poultry is comb morphology, which is a monogenic trait regulated by the same set of genes but with independent mutations in different breeds (e.g. EOMES, MNR2, and SOX5) . Thus, we hypothesize that the same independent selection applied for domestic feathering rate within and across species has resulted in independent mutations in PRLR.
We describe a case–control GWAS that detected a genomic region on the Z chromosome, which is statistically associated with rate of feathering at hatch in turkey. Within this genomic region, we identified a hemizygous 5-bp frameshift deletion in PRLR, which causes the loss of 98 C-terminal AA and is the causal polymorphism for low feathering phenotype in turkey. This is a clear example of similar selection pressures for the same trait (sexing at hatch) in two domestic poultry species that result in two distinct mutations but each affecting the C-terminal end of the same protein, i.e. PRLR. The function of the C-terminal end of this protein, located in the cytoplasm, remains mostly unknown, and further functional studies are necessary to gain more insight in the downstream molecular pathways affected by this mutation.
Siegel PB, Mueller CD, Craig JV. Some phenotypic differences among homozygous, heterozygous, and hemizygous late feathering chicks. Poult Sci. 1957;36:232–9.
Elferink MG, Vallée AAA, Jungerius AP, Crooijmans RPMA, Groenen MAM. Partial duplication of the PRLR and SPEF2 genes at the late feathering locus in chicken. BMC Genomics. 2008;9:391.
Bu GX, Huang G, Fu H, Li J, Huang SM, Wang YJ. Characterization of the novel duplicated PRLR gene at the late-feathering K locus in Lohmann chickens. J Mol Endocrinol. 2013;51:261–76.
Bu G, Ying Wang C, Cai G, Leung FC, Xu M, et al. Molecular characterization of prolactin receptor (cPRLR) gene in chickens: gene structure, tissue expression, promoter analysis, and its interaction with chicken prolactin (cPRL) and prolactin-like protein (cPRL-L). Mol Cell Endocrinol. 2013;370:149–62.
Luo CL, Shen X, Rao YS, Xu HP, Tang J, Sun L, et al. Differences of Z chromosome and genomic expression between early- and late-feathering chickens. Mol Biol Rep. 2012;39:6283–8.
Food and agriculture organization statistical division (FAOSTAT) of the United Nations. http://faostat.fao.org/. Accessed 11 Feb 2018.
Zakrzewska EI, Savage TF. Inhibited feathering: a new dominant sex-linked gene in the turkey. J Hered. 1997;88:238–47.
Renema RA, Sikur VR, Robinson FE, Korver DR, Zuidhof MJ. Effects of nutrient density and age at photostimulation on carcass traits and reproductive efficiency in fast- and slow-feathering turkey hens. Poult Sci. 2008;87:1897–908.
Sikur VR, Robinson FE, Korver DR, Renema RA, Zuidhof MJ. Effects of nutrient density on growth and carcass traits in fast- and slow-feathering female turkeys. Poult Sci. 2004;83:1507–17.
Asmundson VS, Abbott UK. Dominant sex-linked late-feathering in the turkey. J Hered. 1961;52:99–104.
Sickle: a sliding-window, adaptive, quality-based trimming tool for FastQ files https://github.com/najoshi/sickle. Accessed 10 Feb 2018.
Dalloul RA, Long JA, Zimin AV, Aslam L, Beal K, Blomberg L, et al. Multi-platform next-generation sequencing of the domestic turkey (Meleagris gallopavo): genome assembly and analysis. PLoS Biol. 2010;8:e1000475.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25:2078–9.
McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20:1297–303.
Okonechnikov K, Conesa A, Garcia-Alcalde F. Qualimap 2: advanced multi-sample quality control for high-throughput sequencing data. Bioinformatics. 2016;32:292–4.
Garrison E, Marth G. Haplotype-based variant detection from short-read sequencing. arXiv 2012:1207.3907 [q-bio.GN].
Aslam ML, Bastiaansen JW, Elferink MG, Megens HJ, Crooijmans RP, Blomberg A, et al. Whole genome SNP discovery and analysis of genetic diversity in Turkey (Meleagris gallopavo). BMC Genomics. 2012;13:391.
Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81:559–75.
Cingolani P, Platts A, Wang Ie L, Coon M, Nguyen T, Wang L, et al. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w(1118); iso-2; iso-3. Fly (Austin). 2012;6:80–92.
Turner SD. qqman: an R package for visualizing GWAS results using Q-Q and manhattan plots. bioRxiv 2014. https://doi.org/10.1101/005165.
The UniProt C. UniProt: the universal protein knowledgebase. Nucleic Acids Res. 2017;45:D158–69.
Craven AJ, Ormandy CJ, Robertson FG, Wilkins RJ, Kelly PA, Nixon AJ, et al. Prolactin signaling influences the timing mechanism of the hair follicle: analysis of hair growth cycles in prolactin receptor knockout mice. Endocrinology. 2001;142:2533–9.
Sievers F, Wilm A, Dineen D, Gibson TJ, Karplus K, Li W, et al. Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Mol Syst Biol. 2011;7:539.
Xie C, Tammi MT. CNV-seq, a new method to detect copy number variation using high-throughput sequencing. BMC Bioinformatics. 2009;10:80.
Ormandy CJ, Camus A, Barra J, Damotte D, Lucas B, Buteau H, et al. Null mutation of the prolactin receptor gene produces multiple reproductive defects in the mouse. Genes Dev. 1997;11:167–78.
Littlejohn MD, Henty KM, Tiplady K, Johnson T, Harland C, Lopdell T, et al. Functionally reciprocal mutations of the prolactin signalling pathway define hairy and slick cattle. Nat Commun. 2014;5:5861.
Zhao J, Yao J, Li F, Yang Z, Sun Z, Qu L, et al. Identification of candidate genes for chicken early- and late-feathering. Poult Sci. 2016;95:1498–503.
Nakamura A, Ishikawa A, Nagao K, Watanabe H, Uchida M, Kansaku N. Characteristics of reversion to early feathering phenotype in the late feathering line of Nagoya breed chickens. J Poult Sci. 2011;48:155–61.
Cieslak M, Reissmann M, Hofreiter M, Ludwig A. Colours of domestication. Biol Rev Camb Philos Soc. 2011;86:885–99.
Wright D. The genetic architecture of domestication in animals. Bioinform Biol Insights. 2015;9:11–20.
Imsland F, Feng C, Boije H, Bed’hom B, Fillon V, Dorshorst B, et al. The Rose-comb mutation in chickens constitutes a structural rearrangement causng both altered comb morphology and defective sperm motility. PLoS Genet. 2012;8:e1002775.
MAMG and AV conceived and designed the study. MFLD and JMH performed the data analysis. MFLD wrote the manuscript. JMH, RPMAC, AV, JAL, HJM, and MAMG provided useful comments and suggestions and helped to draft the manuscript. Data collection and sequencing were provided by RPMAC and JAL. All authors read and approved the final manuscript.
The authors would like to thank Chiara Bortoluzzi and Vinicius Da Silva for useful input on this work.
The authors declare that they have no competing interests.
Availability of data and materials
The data that support the findings of our study are available under BioProject: PRJNA429537. The sequence data is deposited on SRA, under Accession: SRS2833975.
Consent for publication
Ethics approval and consent to participate
The data used in this study were obtained as part of a data collection in a commercial breeding program. Sample collection and data recording were conducted strictly according to the Dutch law on animal protection and welfare (Gezondheids-en welzijnswet voor dieren).
This research was funded by the STW-Breed4Food Partnership, Project Number 14283: From sequence to phenotype: detecting deleterious variation by prediction of functionality. This study was financially supported by NWO-TTW and the Breed4Food partners Cobb Europe, CRV, Hendrix Genetics and Topigs-Norsvin. Additional funding was provided by the U.S. Department of Agriculture (USDA). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. The use of the HPC cluster was made possible by CAT-AgroFood (Shared Research Facilities Wageningen UR).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.