Identification of key modules and hub genes involved in regulating the feather follicle development of Wannan chickens using WGCNA

Carcass appearance is important economic trait, which affects customers in making purchase decisions. Both density and diameter of feather follicles are two important indicators of carcass appearance. However, the regulatory network and key genes be involved in feather follicle development remain poorly understood. To identify key genes and modules that involved in feather follicle development in chickens, 16 transcriptome datasets of Wannan chickens skin tissue (3 birds at the E9, E11, and E14, respectively, and 7 birds at the 12W) were used for weighted gene co-expression network analysis (WGCNA) analysis, and 12 skin tissue samples (3 birds for each stage) were selected for DEGs analysis. A total of 5,025, 2,337, and 10,623 DEGs were identified in 3 comparison groups, including the E9 vs. E11, the E11 vs. E14, and the E14 vs. 12W. Additionally, 31 co-expression gene modules were identified by WGCNA and the dark-orange, cyan, and blue module were found to be significantly associated with feather follicle development (p < 0.01). In total, 92,898 and 8,448 hub genes were obtained in the dark-orange, cyan, and blue modules, respectively. We focused on the cyan and blue modules, as 6 and 336 hub genes of these modules were identified to overlap with the DEGs of the three comparison groups, respectively. The 6 overlapped genes such as LAMC2, COL6A3, and COL6A2 etc., were over-represented in 12 categories such as focal adhesion and ECM-receptor interaction signaling pathway. Among the 336 genes that overlapped between the blue module and different DEGs comparison groups several genes including WNT7A and WNT9B were enriched in Wnt and ECM-receptor interaction signaling pathway. These results suggested that the LAMC2, COL6A3, COL6A2, WNT7A, and WNT9B genes may play a crucial role in the regulation of feather follicle development in Wannan chickens. Our results provided a reference for the molecular regulatory network and key genes in the development of feather follicles and contribute to molecular breeding for carcass appearance traits in chickens.

Molecular genetic foundation of a sex-linked tailless trait in Hongshan chicken by whole genome data analysis

As a Chinese local chicken breed, Hongshan chickens have 2 kinds of tail feather phenotypes, normal and taillessness. Our previous studies showed that taillessness was a sex-linked dominant trait. Abnormal development of the tail vertebrae could be explained this phenomenon in some chicken breeds. However, the number of caudal vertebrae in rumpless Hongshan chickens was normal, so rumplessness in Hongshan chicken was not related to the development of the caudal vertebrae. Afterwards, we found that rumplessness in Hongshan was due to abnormal development of tail feather rather than abnormal development of caudal vertebrae. In order to understand the genetic foundation of the rumplessness of Hongshan chickens, we compared and reanalyzed 2 sets of data in normal and rumpless Hongshan chickens from our previous studies. By joint analysis of genome-wide selection signature analysis and genome-wide association approach, we found that 1 overlapping gene (EDIL3) and 16 peak genes (ENSGALG00000051843, ENSGALG00000053498, ENSGALG00000054800, KIF27, PTPRD, ENSGALG00000047579, ENSGALG00000041052, ARHGEF28, CAMK4, SERINC5, ENSGALG00000050776, ERCC8, MCC, ADAMTS19, ENSGALG00000053322, CHRNA8) located on the Z chromosome was associated with the rumpless trait. The results of this study furtherly revealed the molecular mechanism of the rumpless trait in Hongshan chickens, and identified the candidate genes associated with this trait. Our results will help to improve the shape of chicken tail feathers and to rise individual economic value in some specific market in China.

In ovo injection of AZD6244 suppresses feather follicle development by the inhibition of ERK and Wnt/β-catenin pathways in goose embryos (Anser cygnoides)

1. Feathers are an important product from poultry, and the state of feather growth and development plays an important role in their economic value.2. In total, 120 eggs were selected for immunoblotting and immunolocalisation experiments of ERK and β-catenin proteins in different developmental stages of goose embryos. The ERK protein was highly expressed in the early stage of goose embryo development, while β-catenin protein was highly expressed in the middle stage of embryo development.3. The 120 eggs were divided into four treatment groups, including an uninjected group (BLANK), a group injected with 100 µl of cosolvent (CK), a group injected with 100 µl of AZD6244 containing cosolvent in a dose of 5 mg/kg AZD6244 containing cosolvent (AZD5) and a group injected with 100 µl of AZD6244 containing cosolvent in a dose of 15 mg/kg AZD6244 containing cosolvent (AZD15). The eggs were injected on the ninth day of embryonic development (E9). Samples were collected at E21.5 to observe feather width, feather follicle diameter, ERK and Wnt/β-catenin pathway protein expression.4. The AZD5 and AZD15 doses were within the embryonic safety range compared to the BLANK and CK groups and had no significant effect on the survival rate and weight at the inflection point, but significantly reduced the feather width and feather follicle diameter (p < 0.05). The AZD6244 treatment inhibited ERK protein phosphorylation levels and blocked the Wnt/β-catenin pathway, which in turn significantly down-regulated the expression levels of FZD4, β-catenin, TCF4 and LEF1 (p < 0.05), with an inhibitory effect in the AZD15 group being more significant. The immunohistochemical results of β-catenin and p-ERK were consistent with Western blot results.5. The small molecule inhibitor AZD6244 regulated the growth and development of feather follicles in goose embryos by the ERK and Wnt/β-catenin pathways.

Effects of Sex on Growth Performance, Carcass Traits, Blood Biochemical Parameters, and Meat Quality of XueShan Chickens

The demand for high-quality chilled chicken has continued to increase in China. Chickens are sexually dimorphic, and to better understand the specific differences in chicken production based on sex, we examined how sex affects growth performance, carcass traits, and meat quality of yellow-feathered chickens. Male and female Xueshan chickens were used as the experimental model. Although males exhibited better growth performance, including body weight (BW), body slope, keel, shank length, and shank girth (p < 0.05), as well as carcass traits, such as dressed weight, leg muscle, and lean meat, females had higher carcass and breast muscle yields (p < 0.05). Males had higher follicle density and yellowness (b*) of the skin and better skin than females (p < 0.05). Among blood biochemical parameters, the serum content of corticosterone (CORT) was higher in males, while those of superoxide dismutase (SOD), glutathione peroxidase (GSH-PX), total antioxidant capacity (T-AOC), and catalase (CAT) were lower in males than in females (p < 0.05). The pH levels, shear force, and moisture content quality were better in male breast meat, while the intramuscular fat content (IMF) was lower in males than in females (p < 0.05). The redness (a*) and moisture content were higher in male leg meat, while the pH, water-loss rate (WLR), lightness (L*), and IMF were lower (p < 0.05). The muscle fiber diameter and cross-sectional area were also higher in males (p < 0.05). Consumers felt that soup of male chicken was better than female (p < 0.05), while mouthfeel and tenderness acceptance of breast meat were different between the sexes. These results indicate that female chickens can be marketed as a whole carcass, while males are more suitable for processed carcass products. This study provides significant insights into the production and processing methodologies of yellow-feathered chickens.

Transcriptome Sequencing and Mass Spectrometry Reveal Genes Involved in the Non-mendelian Inheritance-Mediated Feather Growth Rate in Chicken

The feather growth rate in chickens included early and late feathering. We attempted to characterize the genes and pathways associated with the feather growth rate in chickens that are not in agreement with Mendelian inheritance. Gene expression profiles in the hair follicle tissues of late-feathering cocks (LC), early-feathering cocks (EC), late-feathering hens (LH), and early-feathering hens (EH) were acquired using RNA sequencing (RNA-seq), mass spectrometry (MS), and quantitative reverse transcription PCR (qRT‑PCR). A total of 188 differentially expressed genes (DEGs) were ascertained in EC vs. LC and 538 DEGs were identified in EH vs. LH. We observed that 14 up-regulated genes and 9 down-regulated genes were screened both in EC vs. LC and EH vs. LH. MS revealed that 41 and 138 differentially expressed proteins (DEPs) were screened out in EC vs. LC and EH vs. LH, respectively. Moreover, these DEGs and DEPs were enriched in multiple feather-related pathways, including JAK-STAT, MAPK, WNT, TGF-β, and calcium signaling pathways. qRT-PCR assay showed that the expression of WNT8A was decreased in LC compared with EC, while ALK and GRM4 expression were significantly up-regulated in EH relative to LH. This study helps to elucidate the potential mechanism of the feather growth rate in chickens that do not conform to genetic law.

Transcriptome and Weighted Gene Co-Expression Network Analysis for Feather Follicle Density in a Chinese Indigenous Breed

Feather follicle density plays an important role in appealing to consumers' first impressions when making purchasing decisions. However, the molecular network that contributes to this trait remains largely unknown. The aim of this study was to perform transcriptome and weighted gene co-expression network analyses to determine the candidate genes relating to feather follicle density in Wannan male chickens. In total, five hundred one-day-old Wannan male chickens were kept in a conventional cage system. Feather follicle density was recorded for each bird at 12 weeks of age. At 12 weeks, fifteen skin tissue samples were selected for weighted gene co-expression network analysis, of which six skin tissue samples (three birds in the H group and three birds in the L group) were selected for transcriptome analysis. The results showed that, in total, 95 DEGs were identified, and 56 genes were upregulated and 39 genes were downregulated in the high-feather-follicle-density group when compared with the low-feather-follicle-density group. Thirteen co-expression gene modules were identified. The red module was highly significantly negatively correlated with feather follicle density (p < 0.01), with a significant negative correlation coefficient of -0.72. In total, 103 hub genes from the red module were screened. Upon comparing the 103 hub genes with differentially expressed genes (DEGs), it was observed that 13 genes were common to both sets, including MELK, GTSE1, CDK1, HMMR, and CENPE. From the red module, FOXM1, GTSE1, MELK, CDK1, ECT2, and NEK2 were selected as the most important genes. These genes were enriched in the DNA binding pathway, the heterocyclic compound binding pathway, the cell cycle pathway, and the oocyte meiosis pathway. This study suggests that FOXM1, GTSE1, MELK, CDK1, ECT2, and NEK2 may be involved in regulating the development of feather follicle density in Wannan male chickens. The results of this study reveal the genetic structure and molecular regulatory network of feather follicle density in Wannan male chickens, and provide a basis for further elucidating the genetic regulatory mechanism and identifying molecular markers with breeding value.

Multi-Stage Transcriptome Analysis Revealed the Growth Mechanism of Feathers and Hair Follicles during Induction Molting by Fasting in the Late Stage of Egg Laying

Induced molting is a common method to obtain a new life in laying hens, in which periodic changes in feathers are the prominent feature. Nevertheless, its precise molecular mechanism remains unclear. In this study, feather and hair follicle samples were collected during fasting-induced physiological remodeling for hematoxylin-eosin staining, hormone changes and follicle traits, and transcriptome sequencing. Feather shedding was observed in F13 to R25, while newborns were observed in R3 to R32. Triiodothyronine and tetraiodothyronine were significantly elevated during feather shedding. The calcium content was significantly higher, and the ash content was significantly lower after the changeover. The determination of hair follicle traits revealed an increasing trend in pore density and a decrease in pore diameter after the resumption of feeding. According to RNA-seq results, several core genes were identified, including DSP, CDH1, PKP1, and PPCKB, which may have an impact on hair follicle growth. The focus was to discover that starvation may trigger changes in thyroid hormones, which in turn regulate feather molting through thyroid hormone synthesis, calcium signaling, and thyroid hormone signaling pathways. These data provide a valuable resource for the analysis of the molecular mechanisms underlying the cyclical growth of hair follicles in the skin during induced molting.

Molecular Regulatory Mechanisms in Chicken Feather Follicle Morphogenesis

In China, the sale of freshly slaughtered chickens is becoming increasingly popular in comparison with that of live chickens, and due to this emerging trend, the skin and feather follicle traits of yellow-feathered broilers have attracted a great deal of research attention. The feather follicle originates from the interaction between the epidermis and dermis in the early embryonic stage. Feather follicle morphogenesis is regulated by the Wnt, ectodysplasin (Eda), epidermal growth factor (EGF), fibroblast growth factor (FGF), bone morphogenetic protein (BMP), sonic hedgehog (Shh), Notch, and other signaling pathways that exist in epithelial and mesenchymal cells. The Wnt pathway is essential for feather follicle and feather morphogenesis. Eda interacts with Wnt to induce FGF expression, which attracts mesenchymal cell movement and aggregates to form feather follicle primordia. BMP acts as an inhibitor of the above signaling pathways to limit the size of the feather tract and distance between neighboring feather primordia in a dose-dependent manner. The Notch/Delta pathway can interact with the FGF pathway to promote feather bud formation. While not a part of the early morphogenesis of feather follicles, Shh and BMP signaling are involved in late feather branching. This review summarizes the roles of miRNAs/lncRNA in the regulation of feather follicle and feather growth and development and suggests topics that need to be solved in a future study. This review focuses on the regulatory mechanisms involved in feather follicle morphogenesis and analyzes the impact of SNP sites on feather follicle traits in poultry. This work may help us to understand the molecular regulatory networks influencing feather follicle growth and provide basic data for poultry carcass quality.

Genome-wide association analysis of the primary feather growth traits of duck: identification of potential Loci for growth regulation

The feather is an important epidermal appendage, plays an important role in the life activities of avian specie, and has important economic value. Revealing the molecular regulation mechanism of feather growth has a significant meaning in studying adaptive evolution, physiology, and mating of avian species and also provides a theoretical reference for poultry breeding. In this study, the genome-wide association analysis (GWAS) of 358 ducks was based on primary feather length phenotypic data (28-60 d), length growth rates (LGRs), and maturity scores (60 d) to explore the genetic basis affecting feather growth and maturation. The results showed that, among the primary feather 1 to 5 in ducks, the mean LGR of primary feather 2 was the fastest, with the longest length. The primary feathers in males grew and matured slightly faster than in females. The mean maturity scores of primary feather 10∼7 were higher than primary feather 1 to 3 in ducks. GWAS further showed 116 SNPs associated with feather length traits. In addition, 2 candidate regions (Chr1: 127,407,230-127,524,879 bp and Chr21: 182,061,707-183,616,298 bp) were associated with LGR, which contain total 13 candidate genes (The extremely significant SNPs were mainly located in 2 genes: Chr1: REPS2 and Chr21: PTPRT). Four candidate regions (Chr1: 29,113,036-28,675,018 bp, Chr2: 18,253,612-149,111,290 bp, Chr15: 6,489,774 to 12,138,221 bp and Chr21: 6,578,021-8,472,904 bp) were associated with feather maturity, which contain total 24 candidate genes (The extremely significant SNPs were mainly located in 4 genes: Chr1: IMMP2L, DOCK4 and DDX10, Chr2: LDLRAD4). In conclusion, sex factors influence feather growth and maturity, and the genetic basis of the growth /maturity trait between different feathers is similar. REPS2, PTPRT genes, and IMMP2L, DOCK4, DDX10, and LDLRAD4 are important candidate genes that influence feather growth and maturity, respectively.

Indigenous, Yellow-Feathered Chickens Body Measurements, Carcass Traits, and Meat Quality Depending on Marketable Age

Given an increasing trend in slaughter and chilling for the sale of chickens in China, it is important to determine the marketable age of chickens for chilled sales. This study determined the effects of two marketable ages on the body measurements, carcass traits, and meat quality of yellow-feathered chickens. A total of 360 healthy one-day-old male Xueshan chickens were raised in six pens (straw-covered floor, numbered 1 to 6) and treated in the same manner (free access to food and water) until day 100. Sixty chickens from pens numbered 1 to 3 and 4 to 6 were selected to determine the body measurements, carcass traits, and meat quality at two slaughter ages (90 and 100 days), respectively. One hundred-day-old chickens had a higher body slope, cockscomb, keel, shank lengths, and higher live and dressed weights (p < 0.05). The abdomen skin follicle density, a*(redness) and b*(yellowness) values were higher in 100-day-old chickens (p < 0.05), whereas the 90-day-old chickens were characterized by better spotted skin. For the breast muscle, pH, shear force, a*, moisture, and protein and intramuscular fat contents were lower; moreover, L*(lightness) and b* were higher in 90-day-old chickens. In leg muscles, the pH, shear force, L*, b* and collagen content were lower; furthermore, the a* and moisture contents were higher in 90-day-old chickens (p < 0.05). These findings indicate that two marketable ages both have pros and cons, but 90 days chickens perform better on carcass appearance, and producers can adjust the marketable age to meet needs of different consumers. This study provides a unique idea and theoretical reference for breeding and marketing yellow-feathered chickens.

Transcriptional Characteristics Showed That miR-144-y/FOXO3 Participates in Embryonic Skin and Feather Follicle Development in Zhedong White Goose

Simple Summary

Feather is one of the most valuable and economical products in goose farming and plays a crucial physiological role in birds. For avian biology and the poultry industry, it is essential to comprehend and regulate how skin and feather follicles develop during embryogenesis. This study showed that several key regulatory genes (FOXO3, CTGF, and PTCH1, among others) and miRNAs (miR-144-y) participated in the developmental process of the skin and feather follicles in Zhedong white goose. Our findings are particularly important because they will serve as a valuable resource for upcoming studies on down feathers in agricultural economic growth regarding complex molecular mechanisms and breeding techniques.

Abstract

Skin and feather follicle development are essential processes for goose embryonic growth. Transcriptome and next-generation sequencing (NGS) network analyses were performed to improve the genome of Zhedong White goose and discover the critical genes, miRNAs, and pathways involved in goose skin and feather follicle morphogenesis. Sequencing output generated 6,002,591,668 to 8,675,720,319 clean reads from fifteen libraries. There were 1234, 3024, 4416, and 5326 different genes showing differential expression in four stages, E10 vs. E13, E10 vs. E18, E10 vs. E23, and E10 vs. E28, respectively. The differentially expressed genes (DEGs) were found to be implicated in multiple biological processes and pathways associated with feather growth and development, such as the Wnt signaling pathway, cell adhesion molecules, ECM–receptor interaction signaling pathways, and cell cycle and DNA replication pathways, according to functional analysis. In total, 8276 DEGs were assembled into twenty gene profiles with diverse expression patterns. The reliability of transcriptome results was verified by real-time quantitative PCR by selecting seven DEGs and five miRNAs. The localization of forkhead box O3 (FOXO3), connective tissue growth factor (CTGF), protein parched homolog1 (PTCH1), and miR-144-y by in situ hybridization showed spatial-temporal expression patterns and that FOXO3 and miR-144-y have an antagonistic targeting relationship. The correlation coefficient of FOXO3 and miR-144-y was -0.948, showing a strong negative correlation. Dual-luciferase reporter assay results demonstrated that miR-144-y could bind to the expected location to suppress the expression of FOXO3, which supports that there is a targeting relationship between them. The detections in this report will provide critical insight into the complex molecular mechanisms and breeding practices underlying the developmental characteristics of skin and feather follicles in Zhedong white geese.

Genome assembly and evolutionary analysis of the mandarin duck Aix galericulata reveal strong genome conservation among ducks

The mandarin duck, Aix galericulata, is popular in East Asian cultures and displays exaggerated sexual dimorphism, especially in feather traits during breeding seasons. We generated and annotated the first mandarin duck de novo assembly, which was 1.08 Gb in size and encoded 16,615 proteins. Using a phylogenomic approach calibrated with fossils and molecular divergences, we inferred that the last common ancestor of ducks occurred 13.3–26.7 Ma. The majority of the mandarin duck genome repetitive sequences belonged to the chicken repeat 1 (CR1) retroposon CR1-J2_Pass, which underwent a duck lineage-specific burst. Synteny analyses among ducks revealed infrequent chromosomal rearrangements in which breaks were enriched in LINE retrotransposons and DNA transposons. The calculation of the d_N/d_S ratio revealed that the majority of duck genes were under strong purifying selection. The expanded gene families in the mandarin duck are primarily involved in olfactory perception as well as the development and morphogenesis of feather and branching structures. This new reference genome will improve our understanding of the morphological and physiological characteristics of ducks and provide a valuable resource for functional genomics studies to investigate the feather traits of the mandarin duck.

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.