<figure><img src="image_1.png" /></figure>

HAL
open science

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

Jing Dong, Chuan He, Zhibing Wang, Yanqing Li, Shanshan Li, Lin Tao, Jiebo Chen, Donghua Li, Fenxia Yang, Naibin Li, et al.

<figure><img src="image_2.png" /></figure>

## ► To cite this version:

Jing Dong, Chuan He, Zhibing Wang, Yanqing Li, Shanshan Li, et al. A novel deletion in KRT75L4 mediates the frizzle trait in a Chinese indigenous chicken. Genetics Selection Evolution, 2018, 50 (1), pp.68. 10.1186/s12711-018-0441-7. hal-02405296

HAL Id: hal-02405296

https://hal.science/hal-02405296

Submitted on 11 Dec 2019

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_3.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_4.png" /></figure>

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

Jing Dong<fnref n="1" />†, Chuan He<fnref n="2" />†, Zhibing Wang<fnref n="1" />†, Yanqing Li<fnref n="1" />, Shanshan Li<fnref n="1" />, Lin Tao<fnref n="1" />, Jiebo Chen<fnref n="1" />, Donghua Li<fnref n="1" />, Fenxia Yang<fnref n="1" />, Naibin Li<fnref n="1" />, Quan Zhang<fnref n="1" />, Li Zhang<fnref n="1" />, Guangqin Wang³, Fisayo Akinyemi<fnref n="2" />, He Meng<fnref n="2" />* and Bingwang Du<fnref n="1" />*

**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 F₁ and F₂ 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 F₂ 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-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

<footnote>*Correspondence: menghe@sjtu.edu.cn; dubw@gdou.edu.cn</footnote>

<footnote>†Jing Dong, Chuan He and Zhibing Wang Co-first authors</footnote>

<footnote n="1">¹Animal Science Department of Agricultural College, Guangdong Ocean University, Huguangyan East, Zhanjiang 524088, Guangdong, China</footnote>

<footnote n="2">²Department of Animal Science, School of Agriculture and Biology, Shanghai Jiao Tong University, Shanghai Key Laboratory of Veterinary Biotechnology, Shanghai 200240, China</footnote>

<footnote>Full list of author information is available at the end of the article</footnote>
[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–12] and that the molecular mechanism responsible for frizzle feathers in KRC is different.

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 $F_2$ 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.

<figure><img src="image_5.png" /><figcaption>**Fig. 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</figcaption></figure>
# 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 $F_1$ 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 $F_1$ individuals from KRC × HX and three male and eight female $F_1$ individuals from HX × KRC were mated, respectively, to produce the $F_2$ generation. We conducted a $\chi^2$-test to study inheritance patterns of the frizzle trait by using the R package software [14].

Eighteen $F_1$ and 57 $F_2$ 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* (gal-Gal5) 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 \times 10^{-40}$; and (4) Mendelian errors based on parental genotypes. GWAS on 57 $F_2$ 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 $F_1$, and, as previously
reported [3–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 F₁ generation to produce 486 F₂ offspring. Of these, 113 individuals were homozygous frizzle, 241 were heterozygous frizzle and 132 were homozygous

normal feathered. The $\chi^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 F₁ and F₂ individuals generated 1464 million reads for an average of 3× sequencing depth. After filtering, these reads were aligned to the reference genome, generating about

<figure><img src="image_6.png" /></figure>
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 F₁ 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⁻²⁶), 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.

## 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 F₂ 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 F₂ samples that were sequenced; and thus, was named F deletion.

<table><caption>**Table 1 Markers that are associated with frizzle feather**</caption><thead><tr><th>Chromosome</th><th>Position</th><th>P value</th><th>Nearest gene</th><th>Distance</th></tr></thead><tbody><tr><td>33</td><td>1206216</td><td>1.79E-22</td><td>KRT75L1</td><td>5183</td></tr><tr><td>33</td><td>1235337</td><td>9.46E-26</td><td>KRT75L2</td><td>2919</td></tr><tr><td>33</td><td>1250655</td><td>9.46E-26</td><td>LOC431300</td><td>0</td></tr><tr><td>33</td><td>1280282</td><td>4.09E-24</td><td>LOC768978</td><td>0</td></tr><tr><td>33</td><td>1287383</td><td>2.60E-25</td><td>KRT75L4</td><td>5598</td></tr><tr><td>33</td><td>1287932</td><td>3.75E-21</td><td>KRT75L4</td><td>5049</td></tr><tr><td>33</td><td>1287940</td><td>3.75E-21</td><td>KRT75L4</td><td>5041</td></tr><tr><td>33</td><td>1288000</td><td>1.38E-22</td><td>KRT75L4</td><td>4981</td></tr><tr><td>33</td><td>1288087</td><td>6.19E-18</td><td>KRT75L4</td><td>4894</td></tr><tr><td>33</td><td>1292993</td><td>9.80E-26</td><td>KRT75L4</td><td>0</td></tr><tr><td>33</td><td>1307081</td><td>9.46E-26</td><td>KRT75L4</td><td>0</td></tr><tr><td>33</td><td>1307123</td><td>9.46E-26</td><td>KRT75L4</td><td>0</td></tr><tr><td>33</td><td>1307957</td><td>9.80E-26</td><td>KRT75L4</td><td>0</td></tr><tr><td>33</td><td>1316459</td><td>6.42E-25</td><td>KRT5</td><td>0</td></tr><tr><td>33</td><td>1317239</td><td>1.54E-18</td><td>KRT5</td><td>0</td></tr><tr><td>33</td><td>1318080</td><td>8.31E-20</td><td>KRT5</td><td>0</td></tr><tr><td>33</td><td>1318882</td><td>1.76E-24</td><td>KRT5</td><td>0</td></tr><tr><td>33</td><td>1318889</td><td>1.76E-24</td><td>KRT5</td><td>0</td></tr><tr><td>33</td><td>1318898</td><td>1.76E-24</td><td>KRT5</td><td>0</td></tr><tr><td>33</td><td>1328936</td><td>6.19E-18</td><td>KRT6A</td><td>0</td></tr><tr><td>33</td><td>1356055</td><td>2.78E-18</td><td>KRT7</td><td>0</td></tr><tr><td>33</td><td>1356078</td><td>2.78E-18</td><td>KRT7</td><td>0</td></tr><tr><td>33</td><td>1356154</td><td>2.78E-18</td><td>KRT7</td><td>0</td></tr><tr><td>33</td><td>1357751</td><td>3.65E-19</td><td>KRT7</td><td>0</td></tr><tr><td>33</td><td>1432037</td><td>1.18E-21</td><td>FAIM2</td><td>0</td></tr><tr><td>33</td><td>1432581</td><td>1.79E-22</td><td>FAIM2</td><td>0</td></tr><tr><td>33</td><td>1432677</td><td>1.79E-22</td><td>FAIM2</td><td>0</td></tr><tr><td>33</td><td>1439835</td><td>1.33E-18</td><td>FAIM2</td><td>4398</td></tr><tr><td>33</td><td>1467224</td><td>1.25E-18</td><td>AQP5</td><td>3314</td></tr><tr><td>33</td><td>1510761</td><td>6.15E-18</td><td>LOC100858973</td><td>0</td></tr><tr><td>33</td><td>1522910</td><td>1.75E-18</td><td>GALNT6</td><td>0</td></tr><tr><td>33</td><td>1523316</td><td>2.78E-18</td><td>GALNT6</td><td>0</td></tr></tbody></table>

GWAS identified multiple significant markers. For each marker, position information (chromosome and position), significances (P value) for each genotype at the position, the nearest gene and distances from corresponding markers are provided

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.

## 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 |log₂(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.

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

## 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–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.

<figure><img src="image_7.png" /></figure>
<figure><img src="image_8.png" /></figure>

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–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.
<table><caption>**Table 2** Genotype at the *F*-deletion in the validation groups</caption><thead><tr><th></th><th>-/-</th><th>-/w</th><th>w/w</th></tr></thead><tbody><tr><td>Houdan</td><td>0</td><td>0</td><td>59</td></tr><tr><td>Huaixiang</td><td>0</td><td>0</td><td>44</td></tr><tr><td>Luhua</td><td>0</td><td>0</td><td>22</td></tr><tr><td>Heterozygous F</td><td>0</td><td>126</td><td>0</td></tr><tr><td>Kirin</td><td>206</td><td>0</td><td>0</td></tr></tbody></table>

Genotypes: -/- homozygous *F*-deletion, -/w heterozygous *F*-deletion, w/w homozygous wild type

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

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

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

## Author details

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

Received: 16 April 2018 Accepted: 16 December 2018

Published online: 20 December 2018

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

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

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

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

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.

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.

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.

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.

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.

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.

<footnote n="1">¹ Animal Science Department of Agricultural College, Guangdong Ocean University, Huguangyan East, Zhanjiang 524088, Guangdong, China.</footnote>

<footnote n="2">² Department of Animal Science, School of Agriculture and Biology, Shanghai Jiao Tong University, Shanghai Key Laboratory of Veterinary Biotechnology, Shanghai 200240, China. ³ Zhanjiang Jinsheng Animal Husbandry Science and Technology Ltd., Zhanjiang 524025, Guangdong, China.</footnote>
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.

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.

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.

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

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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

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

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

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.

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.

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.

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

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.

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.

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.

Ready to submit your research? Choose BMC and benefit from:

* fast, convenient online submission

* thorough peer review by experienced researchers in your field

* rapid publication on acceptance

* support for research data, including large and complex data types

* gold Open Access which fosters wider collaboration and increased citations

* maximum visibility for your research: over 100M website views per year

At BMC, research is always in progress.

Learn more biomedcentral.com/submissions

<figure><img src="image_9.png" /></figure>