Skip to main content
Download PDF
  • Short communication
  • Open access
  • Published:

A novel deletion in KRT75L4 mediates the frizzle trait in a Chinese indigenous chicken

Genetics Selection Evolution volume 50, Article number: 68 (2018) Cite this article

Abstract

Background

Highly diversified in morphology and structure, feathers have evolved into various forms. Frizzle feathers, which result from a developmental defect of the feather, are observed in several domestic chicken breeds. The frizzle phenotype is consistent with incomplete dominance of a major gene, but the molecular mechanisms that underlie this phenotype remain obscure. Kirin, a Chinese indigenous chicken breed that originated in the Guangdong province, is famous for its frizzle feathers. The KRT75 gene is considered as the dominant gene responsible for the frizzle trait in several chicken breeds, but this is not the case in the Kirin breed. Thus, the objective of our study was to investigate the genomic region and mutation responsible for this phenotype in this particular breed.

Results

A resource population was produced by crossing Kirin and Huaixiang chickens to produce F1 and F2 generations. DNA samples from 75 frizzle feather and normal feather individuals were sequenced with double-digest genotyping by sequencing (dd-GBS). After the detection of 525,561 high-quality variants, a genome-wide association analysis was carried out and the gene responsible for the frizzle phenotype was localized within the type II α-keratin cluster on chromosome 33. Sanger sequencing was used to screen for mutations in the exons of five genes of this type II α-keratin cluster. A 15-bp deletion in exon 3 of KRT75L4 that showed complete segregation with the frizzle phenotype was detected within the F2 population. Transcriptome sequencing demonstrated that KRT75L4 was expressed but that the transcript was shorter in Kirin than in Huaixiang chickens. In addition, by using Sanger sequencing, we were able to confirm that the deletion was in complete linkage with frizzle feathers.

Conclusions

A deletion in the KRT75L4 gene is responsible for the frizzle feather phenotype in the Kirin chicken. The identification of this mutation, which causes a developmental defect of avian integument appendages, will improve our understanding of the mechanisms that are involved in feather formation.

Background

Birds display a high degree of diversity in feather structure. In chickens, the frizzle trait is a varietal characteristic, and birds that carry this trait have all their contour feathers curling outward and upward [1, 2]. The rachises of frizzle and normal feathers are different i.e., frizzle feathers turn outward relative to the skin, whereas normal body and flight feathers have a dorsal orientation. In frizzle chickens, rectrices and remiges are less affected but have an irregular appearance, and display other modifications, such as thickening of the barbs and barbules, alteration of the hooklets and in some cases, other structural abnormalities [3, 4]. These modifications of the feather structure reduce the insulating effect of feathers [5,6,7].

To date, the frizzle mutation has been reported to occur in a single autosomal gene, denoted as the F gene (see [5]), and shows incomplete dominance inheritance [4, 5]. Although the curling of feathers of the frizzle phenotype is visible in heterozygous individuals, it is even more pronounced in homozygous individuals, but frizzle feathers are more ornamental (the curve is less severe and feathers look like daisy petals) in heterozygous than in homozygous chickens. Because of their enhanced heat dissipation under tropical conditions, both homozygous and heterozygous frizzle chickens have a relatively higher meat and egg production than wild type individuals [8, 9].

Kirin chicken (KRC) is a Chinese indigenous breed of chicken that is noted for its frizzle feathers. The KRC breed originates from the Guangdong province, and is raised exclusively in this region with a warm and humid subtropical climate. Adult KRC have a distinct disorientation of feathers (Fig. 1a–d). In general, the frizzle phenotype is not observed in the down of frizzle dwarf chicken at hatch but becomes more obvious in the second-generation bilateral feathers [10]. In contrast, in KRC, Tao et al. [11] observed the presence of frizzle downy feathers in chicks at hatch (Fig. 1e, f), which suggests a different mechanism for the frizzle trait in this breed. Ng et al. [10] demonstrated that a deletion at the junction between exon 5 and intron 5 of the KRT75 gene causes the frizzle feather phenotype in frizzle dwarf chicken, but this deletion is absent in KRC, which suggests that another gene is involved in the formation of feathers [10,11,12] and that the molecular mechanism responsible for frizzle feathers in KRC is different.

Fig. 1
figure 1

Comparison between frizzle and normal phenotypes. a–d Adult Kirin (frizzle feathers) and Huaixiang (normal feathers) chickens. Adult frizzle feathers curve away from the body, which gives a highly divergent appearance from normal feathers. e One-day-old Kirin and f Huaixiang chickens. The downy feathers appear frizzled, which contrasts with feathers from previously reported breeds

Full size image

To dissect the genetic basis of the frizzle phenotype in KRC, we crossed homozygous frizzle KRC with normal feathered Huaixiang chickens (HX). A double-digest genotyping by sequencing (dd-GBS) based [13] genome-wide association analysis on individuals from the F2 generation, allowed us to identify a candidate region for the F gene, in KRC. Then, we used Sanger sequencing to ascertain the causative mutation in several candidate genes. We also compared the transcriptome profiles of KRC and HX. Finally, we validated our finding in other breeds, e.g. Houdan (H), Luhua (L) and HX.

Methods

Animals and dd-GBS sequencing

A three-generation resource population was produced by crossing purebred KRC with frizzle feathers and HX with normal feathers. The F1 generation was derived by mating one KRC male with 10 HX females (KRC × HX), and by mating one HX male with 10 KRC females (HX × KRC). Three male and eight female F1 individuals from KRC × HX and three male and eight female F1 individuals from HX×KRC were mated, respectively, to produce the F2 generation. We conducted a χ2-test to study inheritance patterns of the frizzle trait by using the R package software [14].

Eighteen F1 and 57 F2 homozygous frizzle feather and normal feather individuals were sampled from 12 families [see Additional file 1: Table S1]. Genomic DNA was extracted following the instructions of the AxyPrep™ DNA Gel Extraction kit, double-digested with different restriction enzyme pairs and sequenced. We found that digestion with the enzyme combination, Pst I and Msp I, was appropriate and highly consistent with in silico digestion-site prediction. Thus, these two restriction enzymes were used to digest total DNA. Adapters marked with specific barcodes were ligated to the restriction-digested overhanging sequences by T4 ligase, following the dd-GBS protocol [13]. DNA samples from the above 75 individuals, marked by different barcode adapters, were mixed and purified according to the manufacturer’s instructions. DNA fragments that carried ligated adapters were amplified with primers that had complementary sequences for each adapter. The constructed 2 × 150 bp libraries were sequenced by the Illumina HiSeq 2000 sequencing system in Shanghai Personal Biotechnology Co., Ltd.

All raw sequences are deposited at NCBI under the BioProject registration number PRJNA445355. Removal of adapter sequences from the reads was done by using AdapterRemoval v2 [15] and the trimmed reads were mapped to the reference genome of Gallus gallus (galGal5) by the Burrows–Wheeler Alignment (BWA) tool v0.7.12 [16] in order to identify variants. High-confidence single nucleotide polymorphisms (SNPs) were called by using the GATK software [17], and the outputs were transformed into PED and MAP files format, which were the standard input files of PLINK-1.9 [18].

Genome-wide association analysis

The dd-GBS based SNPs were filtered out based on the following criteria: (1) minor allele frequency lower than 0.05; (2) SNP missing rate across all individuals higher than 0.2; (3) P value for Hardy–Weinberg equilibrium less than 1 × 10e−40; and (4) Mendelian errors based on parental genotypes. GWAS on 57 F2 samples were performed by using the case-control model in PLINK.

Sanger sequencing

In order to detect all possible causative mutations in five keratin genes (KRT5, KRT6A, KRT75, KRT75L2 and KRT75L4), we used exon sequencing. Primers for these candidate genes were designed by using primer premier 5.0 and are provided in Additional file 2: Table S2. DNA obtained from blood samples of five homozygous normal and five frizzle feather individuals was pooled, PCR-amplified and sequenced, respectively, with three replicates for each pool.

To validate that the mutation was in exon 3 of the KRT75L4 gene, blood was sampled from individuals that were randomly selected from KRC, HX, H and L chicken. Heterozygous frizzle samples (F) were obtained from crosses between KRC and other breeds with normal feather breeds. DNA was extracted and pooled from 49 H, 34 HX, 12 L, 86 F and 104 KRC chickens, resulting in five pools. DNA samples of 10 H, 10 HX, 10 L, 40 F and 102 KRC chickens were PCR-amplified individually. PCR conditions were as follows: 5 min at 94 °C for one cycle, and 30 s at 94 °C, 30 s at 60 °C, and 30 s at 72 °C for 35 cycles. The PCR products were sequenced on an ABI 3730 sequencer.

Transcriptome analysis of frizzle and normal Kirin chicken

We performed the transcriptome analysis of follicle tissues from three KRC and three HX individuals at the embryonic age of E13. Under an RNase-free laminar flow cabinet, feather follicles with skin were removed surgically from the back of each individual, sprayed with RNAsafety, and immediately frozen in liquid nitrogen. Total RNA was extracted with Trizol reagent, and mRNA sequencing was carried out using poly (A)-enriched mRNA. Each of these mRNA samples was converted first into barcoded RNA-seq libraries, and then sequenced on an Illumina HiSeq 2000 sequencer.

The FASTQ sequence of each library was submitted to NCBI (PRJNA445349). Using the Tophat2 [19] program with default parameters, stock-specific transcripts were mapped to the chicken reference genome. The normalized gene expression levels were measured in fragments per kb of exon per million fragments mapped (FPKM) using Cufflinks. Differentially expressed genes (DEG) were identified using the DESeq package [20] based on a P value less than 0.05 and |log2(K/H)| higher or equal to  1, where K and H denote the expression level of each gene in the KRC and HX groups, respectively.

Results

Resource population from crossbreeding between KRC and HX

Reciprocal matings between homozygous KRC and HX resulted in a heterozygous F1, and, as previously reported [3,4,5], the frizzle phenotype was incompletely penetrant in this population. Penetrance was similar regardless of whether the F gene was inherited from a cockerel or a hen, thus ruling out sex-linked inheritance. In order to determine the gene underlying the frizzle trait, six males and 16 females were selected from the F1 generation to produce 486 F2 offspring. Of these, 113 individuals were homozygous frizzle, 241 were heterozygous frizzle and 132 were homozygous normal feathered. The χ2-test (df = 2) showed that the frizzle phenotype is consistent with Mendelian inheritance (P = 0.22) and incompletely dominant in KRC.

QTL mapping by genome-wide association analysis

dd-GBS sequencing applied to 75 F1 and F2 individuals generated 1464 million reads for an average of 3× sequencing depth. After filtering, these reads were aligned to the reference genome, generating about 100,000 effective fragments, which is equivalent to 30 Mb and represents 2.8% of the entire genome with an average depth of 54×. A total of 1,948,658 raw variants were identified, loaded into PLINK 1.9 [18], and filtered according to the criteria presented in the Methods section, which yielded 525,561 variants. The distribution of these variants across the chromosomes was relatively even [see Additional file 3: Figure S1]. Eighteen of the 75 F1 individuals were used as references to remove variants with Mendelian errors.

We applied a case-control model to calculate the association of variants with the frizzle phenotype across the whole genome (see Fig. 2). SNPs that showed a significant association and their nearest gene are listed in Table 1. A majority of the most significant SNPs were located within a region of 130 kb (P = 9.46 × 10−26), which corresponds to intron 5 of the keratin 75-like 4 gene (KRT75L4). The other most significant SNPs were located in the vicinity of the KRT75L2, KRT5, KRT7 and KRT6A genes, which encompass a region between 1.19 and 1.52 Mb on chromosome 33 and are putative candidates for the F locus. In summary, the candidate QTL overlaps (1.19–1.39 Mb) with the cluster of α-keratin genes on chromosome 33 [21]. To our knowledge, the other detected genes, such as FAIM2, AQP5 and GALNT6 are not involved in feather formation.

Fig. 2
figure 2

Genome-wide significance of associations. a Manhattan plot showing the association of all SNPs with the frizzle feather phenotype. Chromosomes 1–33, W and Z are in different colors. The red solid line indicates a significant association with a P value of 1e−10. b Manhattan plot of chromosome 33. The red solid line indicates a significant association with a P value of 1e−10

Full size image
Table 1 Markers that are associated with frizzle feather
Full size table

Identification of a mutation in KRT75L4 in the KRC

The candidate region identified by GWAS lies within a genomic interval on chromosome 33, which covers the cluster of type II α-keratin genes. Although in frizzle dwarf chicken, the causative mutation for the frizzle phenotype is localized in the KRT75 gene [10], previous results indicate that this is not the case in KRC [11]. To detect the putative causative mutations in KRC, we PCR-amplified and sequenced all the exons of five keratin genes (KRT5, KRT6A, KRT75, KRT75L2, and KRT75L4) located within the region detected by GWAS. The KRT75L4 gene has three isoforms (ENSGALT00000086819, ENSGALT00000090244 and ENSGALT00000047324) and the transcript is associated with two gene entries (ENSGALG00000044875 and ENSGALG00000043689). We analyzed the DNA extracted from 10 homozygous F2 individuals, i.e. five with normal and five with frizzle feathers and detected 90 synonymous substitutions, five missense substitutions, two deletions and three insertions (indels). Among these five indels, only one was located within the coding sequence of KRT75L4, i.e. a 15-bp deletion (Fig. 3), which can alter the sequence of the corresponding protein. This deletion showed complete segregation with the frizzle phenotype in the F2 samples that were sequenced; and thus, was named F deletion. All the other variants were excluded because at least one wild type chicken had the same genotype as that in the chickens with frizzle feathers.

Fig. 3
figure 3

Fluorescent peaks of the 15-bp deletion and its flanking sequences. The sequence on chromosome 33 starts at position 1,294,968 bp of the reference genome (Gallus_gallus-5.0). The upper and lower sequences are from HX and KRC individuals, respectively, with a deletion of 15 bp in KRC. The upstream and downstream sequences of the deletion have blue and red backgrounds, respectively

Full size image

Whole-transcriptome sequencing

Whole-transcriptome sequencing produced an average of 49 million reads for each sample. After the trimming and filtering steps, 37.4 to 53.6 million high-quality sequence reads from each of the individuals remained and were aligned to the reference genome. On average, 88% of the reads were successfully aligned.

Of the 94 significant differentially expressed genes (DEG) based on a |log2(K/H)| higher than  1 and a P less than 0.05 [see Additional file 4: Table S4], 50 DEG were up-regulated and 44 were down-regulated in KRC. However, none of these DEG were located in the 1.19–1.52 Mb region detected by GWAS on chromosome 33 or in the type II α-keratin cluster.

The level of expression of KRT75L4 was similar in the KRC and HX groups. However, the sequence of the KRT75L4 transcript (displayed by the integrative genomics viewer (IGV) [22, 23]) differed between the KRC and HX groups (Fig. 4). Due to the F deletion in one exon (exon 3 of ENSGALT00000090244.1 or exon 4 of ENSGALT00000086819.1), the KRT75L4 transcript was comparatively shorter in KRC than in HX, leading to the loss of Val-Leu-Pro-Ala-Ser in the translation product, which thus may have a different function compared to the normal protein. The transcripts of the other keratin genes on chromosome 33 were similar for both KRC and HX, which is concordant with the results from the analysis of the KRC’s genome sequence. Thus, our findings indicate that KRT75L4 is probably the F gene.

Fig. 4
figure 4

Integrative genomics viewer (IGV) display of KRT75L4 transcripts in KRC and HXC breeds. a Part of the alignment tracks of RNA-seq reads in the region of KRT75L4 is shown. The read coverage tracks of HX and KRC groups are in orange and green, respectively. b Zoomed-in image of the second exon in the upper graph. Missing coverage in the middle-front region is observed in the three samples of the KRC group, which is due to the F-deletion that is shown in c

Full size image

Validation in other breeds

To validate furthermore the association detected between KRT75L4 and the frizzle feather phenotype, we PCR-amplified and sequenced the F deletion and its flanking sequences located in KRT75L4 exon 3 from DNA obtained from other breeds of chickens. The F deletion was not found in HX, H and L chickens, which are breeds with normal feathers (Table 2). These results suggest that even if the F deletion is not the causative variant, it is strongly linked to the phenotype of frizzle feathers.

Table 2 Genotype at the F-deletion in the validation groups
Full size table

Discussion

Heat stress can be a major concern for chickens that are raised in warm climates [24, 25]. A reduction in feather coverage increases heat dissipation and body heat irradiation, and such adaptation to high environmental temperatures has been well-documented [26,27,28]. Adomako et al. [8] reported favorable effects of the F locus on production traits, such as egg weight, egg mass, body weight and productivity index. Thus, identification of the F gene is important for breeding programs that are designed to enhance the productive performance of chickens maintained in hot humid climates.

Feathers consist mainly of two types of keratin proteins: α- and β-keratins, which are encoded by multigene families. β-keratins are found only in reptiles and birds, whereas α-keratins exist in all vertebrates [29]. Prum and Williamson [30] attributed the alteration of feather shape to mutations in keratin genes. Cellular and biochemical evidence indicates that α-keratin may have a key role in the early formation of rachides, barbs, and barbules [31]. In addition, the human and mouse homologs of KRT75, KRT6A and KRT6B are expressed in distinct parts of the hair follicle [32,33,34]. The morphological and structural diversity of feathers probably results from different combinations of α- and β-keratin genes in different intra-feather tissues [21].

Analyses of mouse and human cells showed that some α-keratins are involved in intracellular signaling pathways, and that mutations in the corresponding genes may interfere in the construction of the cell cytoskeleton [35, 36]. Although type II keratin genes are located within a single cluster on human chromosome 12 and mouse chromosome 15 [37], studies on avian keratin genes have been challenging because of their conservation and duplications. In spite of considerable efforts to establish a correspondence between individual α-keratin chicken genes and their homologs in mammals [21], they are not included in the most recent chicken genome assembly and annotation release. Ng et al. [10] reported that the loss of the authentic splice site at the exon5/intron 5 junction of KRT75 resulted in a 69-bp in-frame deletion within its coding region and played a key role in frizzle feather development. The gene referred to as KRT75 [10] is actually KRT6A according to RefSeq. Furthermore, KRT75L4 lies within the noncoding region of KRT6A, along with KRT75L1, KRT75L2, KRT5 and other predicted loci. The official full name of KRT75L4 is keratin, type II cytoskeletal 75-like 4, and encodes an uncharacterized protein, with an annotation score of 1 [38]. Our results of transcriptome analysis, combined with other evidence, show that the gene is expressed in chicken tissues, although its function is unknown.

Conclusions

We have demonstrated that a deletion allele in KRT75L4 is the major determinant of frizzle feathers in Kirin and, not KRT6A, as in other breeds of frizzle chickens. This feather defect seems to be amplified in chickens due to the elaborate morphogenesis of their feathers. Thus, identification of another gene that has a role in the morphology of feathers contributes to a better understanding of the genetics and development of feathers. Finally, our work provides information on a key gene involved in the formation of frizzle feathers and should help to explore furthermore the role of the keratin genes in normal and aberrant vertebrate development.

References

  1. Somes RG. Mutations and major variants of plumage and skin in chickens. In: Crawford RD, editor. In poultry breeding and genetics. Amsterdam: Elsevier; 1990. p. 169–208.

    Google Scholar 

  2. Somes RG. Guinea fowl plumage color inheritance, with particular attention on the Dun color. J Hered. 1996;87:138–42.

    Article  Google Scholar 

  3. Landauer W, Dunn LC. The “Frizzle” character of fowls: its expression and inheritance. J Hered. 1930;21:291–305.

    Article  Google Scholar 

  4. Hutt FB. The genetics of the fowl. J Genet. 1930;22:109–27.

    Article  Google Scholar 

  5. Haaren-Kiso AV, Horst P, Valle-Zarate A. The effect of the frizzle gene (F) for the productive adaptability of laying hens under warm and temperate environmental conditions. In: Proceedings of the 18th world’s poultry congress: 4–9 Sept 1988; Nagoya; 1988. p. 381–88.

  6. Mérat P. Pleiotropic and associated effects of major genes. In: Crawford RD, editor. In poultry breeding and genetics. Amsterdam: Elsevier; 1990. p. 429–67.

    Google Scholar 

  7. Yunis R, Cahaner A. The effects of naked neck (Na) and frizzle (F) genes on growth and meat yield of broilers and their interactions with ambient temperatures and potential growth rate. Poult Sci. 1999;78:1347–52.

    Article  CAS  Google Scholar 

  8. Adomako K, Olympio OS, Hagan JK, Hamidu JA. Effect of the frizzle gene (F) on egg production and egg quality of laying hens kept in tropical villages. Br Poult Sci. 2014;55:709–14.

    Article  CAS  Google Scholar 

  9. Mathur PK. Genotype-environment interactions: problems associated with selection for increased production. In: Muir WM, Aggrey SE, editors. In poultry genetics, breeding and biotechnology. Wallingford: CABI Pub; 2003. p. 546.

    Google Scholar 

  10. Ng CS, Wu P, Foley J, Foley A, McDonald ML, Juan WT, et al. The chicken frizzle feather is due to an α-keratin (KRT75) mutation that causes a defective rachis. PLoS Genet. 2012;8:e1002748.

    Article  CAS  Google Scholar 

  11. Tao L, Du BW, Zhang L. The development of frizzled follicle and genetic characteristics of candidate gene KRT75 in frizzled feather chicken. Sci Agric Sin. 2015;48:821–30.

    CAS  Google Scholar 

  12. Dong J, Zhang Q, Zhang L, et al. Genetic characteristics of frizzled and incomplete frizzled feather and SNP detection of KRT75 gene. China Anim Husb Vet Med. 2016;43:3003–10.

    Google Scholar 

  13. Peterson BK, Weber JN, Kay EH, Fisher HS, Hoekstra HE. Double digest RADseq: an inexpensive method for de novo SNP discovery and genotyping in model and non-model species. PLoS One. 2012;7:e37135.

    Article  CAS  Google Scholar 

  14. The R Foundation. The R Project for Statistical Computing. https://www.r-project.org/foundation/. Accessed 30 Nov 2017.

  15. Schubert M, Lindgreen S, Orlando L. AdapterRemoval v2: rapid adapter trimming, identification, and read merging. BMC Res Notes. 2016;9:88.

    Article  Google Scholar 

  16. Li H, Durbin R. Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics. 2009;25:1754–60.

    Article  CAS  Google Scholar 

  17. 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.

    Article  CAS  Google Scholar 

  18. Chang CC, Chow CC, Tellier LCAM, Vattikuti S, Purcell SM, Lee JJ. Second-generation PLINK: rising to the challenge of larger and richer datasets. GigaScience. 2015;4:7.

    Article  Google Scholar 

  19. Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013;14:R36.

    Article  Google Scholar 

  20. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.

    Article  Google Scholar 

  21. Ng CS, Wu P, Fan W, Yan J, Chen CK, Lai YT, et al. Genomic organization, transcriptomic analysis, and functional characterization of avian α- and β-keratins in diverse feather forms. Genome Biol Evol. 2014;6:2258–73.

    Article  Google Scholar 

  22. Robinson JT, Thorvaldsdóttir H, Winckler W, Guttman M, Lander ES, Getz G, Mesirov JP. Integrative genomics viewer. Nat Biotechnol. 2011;29:24–6.

    Article  CAS  Google Scholar 

  23. Thorvaldsdóttir H, Robinson JT, Mesirov JP. Integrative genomics viewer (IGV): high-performance genomics data visualization and exploration. Brief Bioinform. 2013;14:178–92.

    Article  Google Scholar 

  24. Cahaner A, Ajuh JA, Siegmund-Schultze M, Azoulay Y, Druyan S, Zárate AV. Effects of the genetically reduced feather coverage in naked neck and featherless broilers on their performance under hot conditions. Poult Sci. 2008;87:2517–27.

    Article  CAS  Google Scholar 

  25. Chen CF, Gourichon D, Huang NZ, Lee YP, Bordas A, Tixier-Boichard M. Performance comparison of dwarf laying hens segregating for the naked neck gene in temperate and subtropical environments. Genet Sel Evol. 2009;41:13.

    Article  Google Scholar 

  26. Eberhart DE, Washburn KW. Variation in body temperature response of naked neck and normally feathered chickens to heat stress. Poult Sci. 1993;72:1385–90.

    Article  CAS  Google Scholar 

  27. Deeb N, Cahaner A. The effects of naked neck genotypes, ambient temperature, and feeding status and their interactions on body temperature and performance of broilers. Poult Sci. 1999;78:1341–6.

    Article  CAS  Google Scholar 

  28. N’dri AL, Mignon-Grasteau S, Sellier N, Beaumont C, Tixier-Boichard M. Interactions between the naked neck gene, sex, and fluctuating ambient temperature on heat tolerance, growth, body composition, meat quality, and sensory analysis of slow growing meat-type broilers. Livest Sci. 2007;110:33–45.

    Article  Google Scholar 

  29. Maderson PF, Hillenius WJ, Hiller U, Dove CC. Towards a comprehensive model of feather regeneration. J Morphol. 2009;270:1166–208.

    Article  Google Scholar 

  30. Prum RO, Williamson S. Theory of the growth and evolution of feather shape. J Exp Zool. 2001;291:30–57.

    Article  CAS  Google Scholar 

  31. Alibardi L, Toni M. Cytochemical and molecular characteristics of the process of cornification during feather morphogenesis. Prog Histochem Cytochem. 2008;43:1–69.

    Article  CAS  Google Scholar 

  32. Winter H, Langbein L, Praetzel S, Jacobs M, Rogers MA, Leigh IM, et al. A novel human type II cytokeratin, K6hf, specifically expressed in the companion layer of the hair follicle. J Invest Dermatol. 1998;111:955–92.

    Article  CAS  Google Scholar 

  33. Wojcik SM, Longley MA, Roop DR. Discovery of a novel murine keratin 6 (K6) isoform explains the absence of hair and nail defects in mice deficient for K6a and K6b. J Cell Biol. 2001;154:619–30.

    Article  CAS  Google Scholar 

  34. Wang ZL, Wong P, Langbein L, Schweizer J, Coulombe PA. Type II epithelial keratin 6hf (K6hf) is expressed in the companion layer, matrix, and medulla in anagen-stage hair follicles. J Invest Dermatol. 2003;121:1276–82.

    Article  CAS  Google Scholar 

  35. Fuchs E, Cleveland DW. A structural scaffolding of intermediate filaments in health and disease. Science. 1998;279:514–9.

    Article  CAS  Google Scholar 

  36. Pan X, Hobbs RP, Coulombe PA. The expanding significance of keratin intermediate filaments in normal and diseased epithelia. Curr Opin Cell Biol. 2013;25:47–56.

    Article  CAS  Google Scholar 

  37. Hesse M, Zimek A, Weber K, Magin TM. Comprehensive analysis of keratin gene clusters in humans and rodents. Eur J Cell Biol. 2004;83:19–26.

    Article  CAS  Google Scholar 

  38. Apweiler R, Bairoch A, Wu CH, Barker WC, Boeckmann B, Ferro S, et al. UniProt: the universal protein knowledgebase. Nucleic Acids Res. 2004;32:D115–9.

    Article  CAS  Google Scholar 

Download references

Authors’ contributions

HM and BD conceived and designed the experiments. JD, ZW, and DL performed the experiments. JD, CH, ZW and LT analyzed the data. YL, SL, JC and FY contributed reagents, materials and tools. NL, QZ, LZ and GW collected the samples. CH, JD and HM wrote the manuscript. CH and FA revised the manuscript. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

Raw reads (base-called FASTQ) used for dd-GBS project and RNA-seq reads are available from the European Nucleotide Archive, project PRJNA445355 and project PRJNA445349, respectively.

Consent for publication

Not applicable.

Ethics approval and consent to participate

All animal experiments were conducted according to the guidelines established by the ethics committee for the Care and Use of Laboratory Animals in Guangdong Ocean University.

Funding

This study was supported by program grant funding from a key project of Zhanjiang modern agriculture (2016A03010): “Strain breeding and molecular mechanism research on resistance to high temperature characteristics of Kirin chicken (frizzle fowl)”. The funding bodies had no role neither in the design of the study, nor in the collection, analysis, interpretation of data, and writing the manuscript.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Author notes

  1. Jing Dong, Chuan He and Zhibing Wang Co-first authors

Authors and Affiliations

  1. Animal Science Department of Agricultural College, Guangdong Ocean University, Huguangyan East, Zhanjiang, 524088, Guangdong, China

    Jing Dong, Zhibing Wang, Yanqing Li, Shanshan Li, Lin Tao, Jiebo Chen, Donghua Li, Fenxia Yang, Naibin Li, Quan Zhang, Li Zhang & Bingwang Du

  2. Department of Animal Science, School of Agriculture and Biology, Shanghai Jiao Tong University, Shanghai Key Laboratory of Veterinary Biotechnology, Shanghai, 200240, China

    Chuan He, Fisayo Akinyemi & He Meng

  3. Zhanjiang Jinsheng Animal Husbandry Science and Technology Ltd., Zhanjiang, 524025, Guangdong, China

    Guangqin Wang

Authors
  1. Jing Dong
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Chuan He
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Zhibing Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Yanqing Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Shanshan Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Lin Tao
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Jiebo Chen
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Donghua Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Fenxia Yang
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Naibin Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  11. Quan Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  12. Li Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  13. Guangqin Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  14. Fisayo Akinyemi
    View author publications

    You can also search for this author in PubMed Google Scholar

  15. He Meng
    View author publications

    You can also search for this author in PubMed Google Scholar

  16. Bingwang Du
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding authors

Correspondence to He Meng or Bingwang Du.

Additional files

Additional file 1: Table S1.

Pedigree and phenotypes of the GWAS population. Description: Columns A and B = sire and dam IDs, respectively; column C = sample IDs; column D = phenotype of the samples.

Additional file 2: Table S2.

Sequence of the primers used to screen all the exons of the five candidate genes.

Additional file 3: Figure S1.

Distribution of high-quality SNPs across the whole genome generated by dd-GBS. Each bar represents a chromosome, with a length proportional to the physical length (Mbp) in the Gallus–gallus 5.0 assembly. Yellow and green bars indicate SNP density in 20,000-bp regions.

Additional file 4: Table S4.

Summary of the differentially expressed genes detected between KRC and HX. Column A = Ensembl IDs of differentially expressed genes; columns B and C = average RPKM of HX and KRC groups, respectively; column D = fold-change both groups; column E = significance P-value.; columns F, G and H = chromosome number and positions of the DEG; column I = gene names.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Dong, J., He, C., Wang, Z. et al. A novel deletion in KRT75L4 mediates the frizzle trait in a Chinese indigenous chicken. Genet Sel Evol 50, 68 (2018). https://doi.org/10.1186/s12711-018-0441-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12711-018-0441-7

Genetics Selection Evolution

ISSN: 1297-9686