Identification of candidate genes for chicken early- and late-feathering

The prolactin receptor: A cross-species comparison of gene structure, transcriptional regulation, tissue-specificity, and genetic variation

The conserved and multifaceted functions of prolactin (PRL) are coordinated through varied distribution and expression of its cell-surface receptor (PRLR) across a range of tissues and physiological states. The resultant heterogeneous expression of PRLR mRNA and protein across different organs and cell types supports a wide range of PRL-regulated processes including reproduction, lactation, development, and homeostasis. Genetic variation within the PRLR gene also accounts for several phenotypes impacting agricultural production and human pathology. The goal of this review is to highlight the many elements that control differential expression of the PRLR across tissues, and the various phenotypes that exist across species due to variation in the PRLR gene.

Artificial selection footprints in indigenous and commercial chicken genomes

BACKGROUND

Although many studies have been done to reveal artificial selection signatures in commercial and indigenous chickens, a limited number of genes have been linked to specific traits. To identify more trait-related artificial selection signatures and genes, we re-sequenced a total of 85 individuals of five indigenous chicken breeds with distinct traits from Yunnan Province, China.

RESULTS

We found 30 million non-redundant single nucleotide variants and small indels (< 50 bp) in the indigenous chickens, of which 10 million were not seen in 60 broilers, 56 layers and 35 red jungle fowls (RJFs) that we compared with. The variants in each breed are enriched in non-coding regions, while those in coding regions are largely tolerant, suggesting that most variants might affect cis-regulatory sequences. Based on 27 million bi-allelic single nucleotide polymorphisms identified in the chickens, we found numerous selective sweeps and affected genes in each indigenous chicken breed and substantially larger numbers of selective sweeps and affected genes in the broilers and layers than previously reported using a rigorous statistical model. Consistent with the locations of the variants, the vast majority (~ 98.3%) of the identified selective sweeps overlap known quantitative trait loci (QTLs). Meanwhile, 74.2% known QTLs overlap our identified selective sweeps. We confirmed most of previously identified trait-related genes and identified many novel ones, some of which might be related to body size and high egg production traits. Using RT-qPCR, we validated differential expression of eight genes (GHR, GHRHR, IGF2BP1, OVALX, ELF2, MGARP, NOCT, SLC25A15) that might be related to body size and high egg production traits in relevant tissues of relevant breeds.

CONCLUSION

We identify 30 million single nucleotide variants and small indels in the five indigenous chicken breeds, 10 million of which are novel. We predict substantially more selective sweeps and affected genes than previously reported in both indigenous and commercial breeds. These variants and affected genes are good candidates for further experimental investigations of genotype-phenotype relationships and practical applications in chicken breeding programs.

Identification of Duplication Genotypes of the Feathering Rate Gene in Chicken by a Multiplex PCR Following Electrophoresis and/or Sanger Sequencing

Sex-linked phenotypes of late feathering (LF) and early feathering (EF) are controlled by a pair of alleles K and k⁺. Autosexing based on the feathering rate is widely used in poultry production. It is reported that a tandem duplication of 176,324 base pairs linked to the K locus is responsible for LF expression and could be used as a molecular marker to detect LF chicken. So far, there is no genotyping method that can accurately and stably identify the LF homozygote and heterozygote in all chicken breeds. In the present study, a multiplex PCR test was developed to identify EF, LF homozygote, and heterozygote according to electrophoretic bands and the relative height of the peaks by Sanger sequencing. We tested 413 chickens of six native Chinese breeds with this method. The identification was consistent with the sex and phenotype records of the chickens. Band density analysis was performed, and the results supported our genotyping using the new assay. In order to further verify the accuracy of this test in distinguishing homozygote and heterozygote males, 152 LF males were mated with EF females, and the results of the offspring's phenotypes were consistent with our expectations. Our results support tandem duplication as molecular markers of LF, and this new test is applicable to all LF chickens associated with tandem duplication.

Genome-Wide Association Study and Transcriptome Differential Expression Analysis of the Feather Rate in Shouguang Chickens

The feather rate phenotype in chicks, including early-feathering and late-feathering phenotypes, are widely used as a sexing system in the poultry industry. The objective of this study was to obtain candidate genes associated with the feather rate in Shouguang chickens. In the present study, we collected 56 blood samples and 12 hair follicle samples of flight feathers from female Shouguang chickens. Then we identified the chromosome region associated with the feather rate by genome-wide association analysis (GWAS). We also performed RNA sequencing and analyzed differentially expressed genes between the early-feathering and late-feathering phenotypes using HISAT2, StringTie, and DESeq2. We identified a genomic region of 10.0–13.0 Mb of chromosome Z, which is statistically associated with the feather rate of Shouguang chickens at one-day old. After RNA sequencing analysis, 342 differentially expressed known genes between the early-feathering (EF) and late-feathering (LF) phenotypes were screened out, which were involved in epithelial cell differentiation, intermediate filament organization, protein serine kinase activity, peptidyl-serine phosphorylation, retinoic acid binding, and so on. The sperm flagellar 2 gene (SPEF2) and prolactin receptor (PRLR) gene were the only two overlapping genes between the results of GWAS and differential expression analysis, which implies that SPEF2 and PRLR are possible candidate genes for the formation of the chicken feathering phenotype in the present study. Our findings help to elucidate the molecular mechanism of the feather rate in chicks.

Whole-genome resequencing reveals aberrant autosomal SNPs affect chicken feathering rate

Previous studies have shown that the feather growth rate of chicks is determined by two alleles located on the sex chromosome Z; however, in chicken production, feathering is usually not consistently controlled by the sex chromosome. To identify whether the feathering rate is related to autosomal inheritance, whole-genome resequencing was performed in eight chickens with slow- and fast-feathering rate. A total of 54,984 autosomal single nucleotide polymorphisms (SNPs) were identified, including 393 and 376 exonic SNPs in slow-feathering and fast-feathering chickens, respectively. Mutated genes were mainly involved in response to stimuli and growth and reproduction processes. Mutated genes related to slow-feathering rate were mainly involved in wingless-type MMTV integration site signaling pathway and mitogen-activated protein kinase signaling pathway, whereas mutated genes associated with fast-feathering rate were primarily enriched in autophagy, calcium signaling pathway, extracellular matrix-receptor interaction, and Focal adhesion processes. Importantly, two SNPs, involved in feather development, were found in the exonic regions of Wnt signaling genes. These results shed new light on the relationship between genetic mutation and feather growth rate from the perspective of autosomal inheritance and may have economic significance in chicken breeding.

Changes in prolactin receptor homodimer availability may cause late feathering in chickens

Chicken early (EF) and late feathering (LF) are sex-linked phenotypes conferred by wild-type k⁺ and dominant K alleles on chromosome Z, respectively. Besides prolactin (PRL) receptor (PRLR) and sperm flagellar 2 (SPEF2) genes, the K allele contains a fusion gene in which partially duplicated PRLR (dPRLR) and SPEF2 (dSPEF2) genes are linked in a tail-to-tail manner. The causative dPRLR gene encodes a C-terminal truncated receptor. LF chickens have short or no primaries at hatching; however, their feather growth rate is higher than that of EF chickens. This study aimed to elucidate the molecular basis of the K allele's biphasic effect on feather development. By 3'RACE and RT-PCR analyses, we demonstrated that dSPEF2 gene transcription occurred beyond all coding exons of the dPRLR gene on the opposite strand and that dPRLR mRNA was less abundant than PRLR mRNA. In addition, a 5'UTR splice variant (SPV) of PRL receptor mRNAs was increased in LF chickens. In vitro expression analysis of 5'UTR linked to the luciferase reporter gene revealed higher translation efficiency of SPV. RT-qPCR showed that the dPRLR mRNA level was higher in embryos; conversely, SPV was higher in hatched chickens, as was dSPEF2 mRNA. These findings suggest that the K allele inhibits feather development at the fetal stage by expressing dPRLR to attenuate PRLR function and promotes feather growth after hatching by increasing PRLR through dSPEF2 mRNA expression. Increased SPV may cause greater feather growth than that in EF chickens by increasing the availability of PRLR homodimers and enhancing PRL signaling.

Analysis of a genetic factors contributing to feathering phenotype in chickens

In the current study, we sought to determine whether or not the endogenous retroviral ev21 influences feathering type of chickens, and if one mutation locus in the unoccupied repeat (UR) region can be used to predict the corresponding feathering type and genotype. The distribution of ev21 as well as the mutation locus in UR and occupied site (OR) regions was detected in HY-line gray progenitor (HYGP) 4 lines, HY-line brown (HYB) and Taihang chickens (TH). Furthermore, a detection method for the genotype resulting in late feathering (LF) phenotype was developed by double PCR using C line of HYGP, C line of Dawu progenitor, commercial line of HY-line gray (HYG) males, LF males of TH and Bashang long-tail chickens (BS). Results indicated that a product of 7590 bp from the long fragment amplification was observed to be a partial segment of ev21, and was linked with the LF phenotype in HYGP but not in HYB and TH chickens. A total of 2 of 35 males and 10 of 29 females of TH LF chickens were found to be ev21 negative. HaeIII RFLP mutations of 1450 bp of UR, 1440 bp of OR, and 538 bp in the UR and OR common region were analyzed, and genotypic features at the locus correlated with the feathering type phenotype in HYGP, but exhibited no significant effects in HYB and TH chickens. The cut-off of relative intensity of 857 and 1305 bp from the double PCR for distinction between homozygous and heterozygous LF males was 1.37. In conclusion, ev21 and the HaeIII RFLP patterns within the locus in UR cannot be used for prediction of feathering type phenotypes in Chinese heritage chickens. However, the partial duplication of PRLR and SPEF2 were able to predict the LF phenotype. Therefore, the double PCR detecting products of 857 and 1305 bp described herein could be used for the accurate identification of genotypes influencing feathering type.

Endogenous viral gene ev21 is not responsible for the expression of late feathering in chickens

The late-feathering (LF) gene K on the Z chromosome is an important gene in the chicken industry, which is frequently utilized for the feather sexing, a type of autosexing, of neonatal chicks. The K gene is closely associated with the endogenous ev21 gene from an avian leukosis virus and the incomplete duplication (ID) of prolactin receptor (PRLR) and sperm flagellar protein 2 (SPEF2) genes, and ev21 has been used as a molecular marker to detect LF birds. In the present study, a comprehensive survey for the presence or absence of ev21 and ID across 1,994 birds from 52 chicken breeds, three commercial hybrid groups, and the Red Jungle Fowl revealed that almost all LF breeds have both ev21 and ID. However, only one LF breed (Ingie) has only ID and no ev21. Moreover, this study revealed that almost all early (normal)-feathering (EF) breeds lack both ev21 and ID, but only one breed (White Plymouth Rock) included EF birds with ev21 but no ID. Therefore, regarding LF expression, the results indicated that ID is responsible, but ev21 is not required. Henceforth, ID should be used as a molecular marker to detect LF birds instead of ev21. Because ev21 contains the full genome of an avian leukosis virus, there is a risk of disease development in breeds with this gene. Therefore, the Ingie breed, which has no ev21 at the K locus, represents excellent material for the establishment of new LF stocks.

Characterization of microRNA and mRNA expression profiles in skin tissue between early-feathering and late-feathering chickens

Background

Early feathering and late feathering in chickens are sex-linked phenotypes, which have commercial application in the poultry industry for sexing chicks at hatch and have important impacts on performance traits. However, the genetic mechanism controlling feather development and feathering patterns is unclear. Here, miRNA and mRNA expression profiles in chicken wing skin tissues were analysed through high-throughput transcriptomic sequencing, aiming to understand the biological process of follicle development and the formation of different feathering phenotypes.

Results

Compared to the N1 group with no primary feathers extending out, 2893 genes and 31 miRNAs displayed significantly different expression in the F1 group with primary feathers longer than primary-covert feathers, and 1802 genes and 11 miRNAs in the L2 group displayed primary feathers shorter than primary-covert feathers. Only 201 altered genes and 3 altered miRNAs were identified between the N1 and L2 groups (fold change > 2, q value < 0.01). Both sequencing and qPCR tests revealed that PRLR was significantly decreased in the F1 and L2 groups compared to the N1 group, whereas SPEF2 was significantly decreased in the F1 group compared to the N1 or L2 group. Functional analysis revealed that the altered genes or targets of altered miRNAs were involved in multiple biological processes and pathways related to feather growth and development, such as the Wnt signalling pathway, the TGF-beta signalling pathway, the MAPK signalling pathway, epithelial cell differentiation, and limb development. Integrated analysis of miRNA and mRNA showed that 14 pairs of miRNA-mRNA negatively interacted in the process of feather formation.

Conclusions

Transcriptomic sequencing of wing skin tissues revealed large changes in F1 vs. N1 and L2 vs. N1, but few changes in F1 vs. L2 for both miRNA and mRNA expression. PRLR might only contribute to follicle development, while SPEF2 was highly related to the growth rate of primary feathers or primary-covert feathers and could be responsible for early and late feather formation. Interactions between miR-1574-5p/NR2F, miR-365-5p/JAK3 and miR-365-5p/CDK6 played important roles in hair or feather formation. In all, our results provide novel evidence to understand the molecular regulation of follicle development and feathering phenotype.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4773-z) contains supplementary material, which is available to authorized users.

Early and late feathering in turkey and chicken: same gene but different mutations

BACKGROUND

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.

RESULTS

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.

CONCLUSIONS

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.