Transcriptomic analyses of regenerating adult feathers in chicken

The potential influence of genome-wide adaptive divergence on conservation translocation outcome in an isolated greater sage-grouse population

Conservation translocations are an important conservation tool commonly employed to augment declining or reestablish extirpated populations. One goal of augmentation is to increase genetic diversity and reduce the risk of inbreeding depression (i.e., genetic rescue). However, introducing individuals from significantly diverged populations risks disrupting coadapted traits and reducing local fitness (i.e., outbreeding depression). Genetic data are increasingly more accessible for wildlife species and can provide unique insight regarding the presence and retention of introduced genetic variation from augmentation as an indicator of effectiveness and adaptive similarity as an indicator of source and recipient population suitability. We used 2 genetic data sets to evaluate augmentation of isolated populations of greater sage-grouse (Centrocercus urophasianus) in the northwestern region of the species range (Washington, USA) and to retrospectively evaluate adaptive divergence among source and recipient populations. We developed 2 statistical models for microsatellite data to evaluate augmentation outcomes. We used one model to predict genetic diversity after augmentation and compared these predictions with observations of genetic change. We used the second model to quantify the amount of observed reproduction attributed to transplants (proof of population integration). We also characterized genome-wide adaptive divergence among source and recipient populations. Observed genetic diversity (HO = 0.65) was higher in the recipient population than predicted had no augmentation occurred (HO = 0.58) but less than what was predicted by our model (HO = 0.75). The amount of shared genetic variation between the 2 geographically isolated resident populations increased, which is evidence of periodic gene flow previously assumed to be rare. Among candidate adaptive genes associated with elevated fixation index (FST) (143 genes) or local environmental variables (97 and 157 genes for each genotype-environment association method, respectively), we found clusters of genes with related functions that may influence the ability of transplants to use local resources and navigate unfamiliar environments and their reproductive potential, all possible reasons for low genetic retention from augmentation.

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.

Transition from natal downs to juvenile feathers: conserved regulatory switches in Neoaves

The transition from natal downs for heat conservation to juvenile feathers for simple flight is a remarkable environmental adaptation process in avian evolution. However, the underlying epigenetic mechanism for this primary feather transition is mostly unknown. Here we conducted time-ordered gene co-expression network construction, epigenetic analysis, and functional perturbations in developing feather follicles to elucidate four downy-juvenile feather transition events. We discovered that LEF1 works as a key hub of Wnt signaling to build rachis and converts radial downy to bilateral symmetry. Extracellular matrix reorganization leads to peripheral pulp formation, which mediates epithelial -mesenchymal interactions for branching morphogenesis. ACTA2 compartments dermal papilla stem cells for feather cycling. Novel usage of scale keratins strengthens feather sheath with SOX14 as the epigenetic regulator. We found this primary feather transition largely conserved in chicken (precocious) and zebra finch (altricial) and discussed the possibility that this evolutionary adaptation process started in feathered dinosaurs.

The idiosyncratic genome of Korean long-tailed chicken as a valuable genetic resource

Today, breeds with ornamental traits such as exceptionally long tail feathers are economically valuable. However, the genetic basis of long-tail feathers is yet to be understood. To provide better understanding of long tail feathers, we sequenced Korean long-tailed chicken (KLC) genomes and compared them with genomes of other chicken breeds. We first analyzed the genome structure of KLC and its genomic relationship with other chickens and observed unique characteristics. Subsequently, we searched for genomic regions under selection. Feather keratin 1-like enriched region and several genes were found to have novel putative functions and effects on the long tail trait in KLC. Our findings support the value of KLC as a unique genetic resource and cast light on the genetic basis of long tail traits in avian species. We expect this novel knowledge to provide new genomic evidence and options for designing and implementing genetic improvements of ornamental chicken productivity through precision crossbreeding aids.

Embryonic transcriptome unravels mechanisms and pathways underlying embryonic development with respect to muscle growth, egg production, and plumage formation in native and broiler chickens

Background: Muscle development, egg production, and plumage colors are different between native and broiler chickens. The study was designed to investigate why improved Aseel (PD4) is colorful, stronger, and grew slowly compared with the control broiler (CB).

Methods: A microarray was conducted using the 7th-day embryo (7EB) and 18th-day thigh muscle (18TM) of improved Aseel and broiler, respectively. Also, we have selected 24 Gallus gallus candidate reference genes from NCBI, and total RNA was isolated from the broiler, improved Aseel embryo tissues, and their expression profiles were studied by real-time quantitative PCR (qPCR). Furthermore, microarray data were validated with qPCR using improved Aseel and broiler embryo tissues.

Results: In the differential transcripts screening, all the transcripts obtained by microarray of slow and fast growth groups were screened by fold change ≥ 1 and false discovery rate (FDR) ≤ 0.05. In total, 8,069 transcripts were differentially expressed between the 7EB and 18TM of PD4 compared to the CB. A further analysis showed that a high number of transcripts are differentially regulated in the 7EB of PD4 (6,896) and fewer transcripts are differentially regulated (1,173) in the 18TM of PD4 compared to the CB. On the 7th^- and 18th-day PD4 embryos, 3,890, 3,006, 745, and 428 transcripts were up- and downregulated, respectively. The commonly up- and downregulated transcripts are 91 and 44 between the 7th- and 18th-day of embryos. In addition, the best housekeeping gene was identified. Furthermore, we validated the differentially expressed genes (DEGs) related to muscle growth, myostatin signaling and development, and fatty acid metabolism genes in PD4 and CB embryo tissues by qPCR, and the results correlated with microarray expression data.

Conclusion: Our study identified DEGs that regulate the myostatin signaling and differentiation pathway; glycolysis and gluconeogenesis; fatty acid metabolism; Jak-STAT, mTOR, and TGF-β signaling pathways; tryptophan metabolism; and PI3K-Akt signaling pathways in PD4. The results revealed that the gene expression architecture is present in the improved Aseel exhibiting embryo growth that will help improve muscle development, differentiation, egg production, protein synthesis, and plumage formation in PD4 native chickens. Our findings may be used as a model for improving the growth in Aseel as well as optimizing the growth in the broiler.

Gut microbiota and transcriptome analysis reveals a genetic component to dropping moisture in chickens

High dropping moisture (DM) in poultry production has deleterious effects on the environment, feeding cost, and public health of people and animals. To explore the contributing genetic components, we classified DM of 67-wk-old Rhode Island Red (RIR) hens at 4 different levels and evaluated the underlying genetic heritability. We found the heritability of DM to be 0.219, indicating a moderately heritable trait. We then selected chickens with the highest and lowest DM levels. Using transcriptome, we only detected 12 differentially expressed genes (DEGs) between these 2 groups from the spleen, and 1,507 DEGs from intestinal tissues (jejunum and cecum). The low number of DEGs observed in the spleen suggests that differing moisture levels are not attributed to pathogenic infection. Fourteen of the intestinal high expressed genes are associated with water-salt metabolism (WSM). We also investigated the gut microbial composition by 16S rRNA gene amplicon sequencing. Six different microbial operational taxonomic units (OTUs) (Cetobacterium, Sterolibacterium, Elusimicrobium, Roseburia, Faecalicoccus, and Megamonas) between the 2 groups from jejunum and cecum are potentially biomarkers related to DM levels. Our results identify a genetic component to chicken DM, and can guide breeding strategies.

Establishment of a culture model for the prolonged maintenance of chicken feather follicles structure in vitro

Protocols allowing the in vitro culture of human hair follicles in a serum free-medium up to 9 days were developed 30 years ago. By using similar protocols, we achieved the prolonged maintenance in vitro of juvenile feather follicles (FF) microdissected from young chickens. Histology showed a preservation of the FF up to 7 days as well as feather morphology compatible with growth and/or differentiation. The integrity of the FF wall epithelium was confirmed by transmission electron microscopy at Day 5 and 7 of culture. A slight elongation of the feathers was detected up to 5 days for 75% of the examined feathers. By immunochemistry, we demonstrated the maintenance of expression and localization of two structural proteins: scaffoldin and fibronectin. Gene expression (assessed by qRT-PCR) of NCAM, LCAM, Wnt6, Notch1, and BMP4 was not altered. In contrast, Shh and HBS1 expression collapsed, DKK3 increased, and KRT14 transiently increased upon cultivation. This indicates that cultivation modifies the mRNA expression of a few genes, possibly due to reduced growth or cell differentiation in the feather, notably in the barb ridges. In conclusion, we have developed the first method that allows the culture and maintenance of chicken FF in vitro that preserves the structure and biology of the FF close to its in vivo state, despite transcriptional modifications of a few genes involved in feather development. This new culture model may serve to study feather interactions with pathogens or toxics and constitutes a way to reduce animal experimentation.

The evolutionary history and mechanistic basis of female ornamentation in a tropical songbird

Ornamentation, such as the showy plumage of birds, is widespread among female vertebrates, yet the evolutionary pressures shaping female ornamentation remain uncertain. In part this is due to a poor understanding of the mechanistic route to ornamentation in females. To address this issue, we evaluated the evolutionary history of ornament expression in a tropical passerine bird, the White-shouldered Fairywren, whose females, but not males, strongly vary between populations in occurrence of ornamented black-and-white plumage. We first use phylogenomic analysis to demonstrate that female ornamentation is derived and that female ornamentation evolves independently of changes in male plumage. We then use exogenous testosterone in a field experiment to induce partial ornamentation in naturally unornamented females. By sequencing the transcriptome of experimentally induced ornamented and natural feathers, we identify genes expressed during ornament production and evaluate the degree to which female ornamentation in this system is associated with elevated testosterone, as is common in males. We reveal that some ornamentation in females is linked to testosterone and that sexes differ in ornament-linked gene expression. Lastly, using genomic outlier analysis we identify a candidate melanogenesis gene that lies in a region of high genomic divergence among populations that is also differentially expressed in feather follicles of different female plumages. Taken together, these findings are consistent with sex-specific selection favoring the evolution of female ornaments and demonstrate a key role for testosterone in generating population divergence in female ornamentation through gene regulation. More broadly, our work highlights similarities and differences in how ornamentation evolves in the sexes.

Analysis and comparison of protein secondary structures in the rachis of avian flight feathers

Avians have evolved many different modes of flying as well as various types of feathers for adapting to varied environments. However, the protein content and ratio of protein secondary structures (PSSs) in mature flight feathers are less understood. Further research is needed to understand the proportions of PSSs in feather shafts adapted to various flight modes in different avian species. Flight feathers were analyzed in chicken, mallard, sacred ibis, crested goshawk, collared scops owl, budgie, and zebra finch to investigate the PSSs that have evolved in the feather cortex and medulla by using nondestructive attenuated total reflection Fourier transform infrared spectroscopy (ATR-FTIR). In addition, synchrotron radiation-based, Fourier transform infrared microspectroscopy (SR-FTIRM) was utilized to measure and analyze cross-sections of the feather shafts of seven bird species at a high lateral resolution to resolve the composition of proteins distributed within the sampled area of interest. In this study, significant amounts of α-keratin and collagen components were observed in flight feather shafts, suggesting that these proteins play significant roles in the mechanical strength of flight feathers. This investigation increases our understanding of adaptations to flight by elucidating the structural and mechanistic basis of the feather composition.

Vertebrate keratinization evolved into cornification mainly due to transglutaminase and sulfhydryl oxidase activities on epidermal proteins: An immunohistochemical survey

The epidermis of vertebrates forms an extended organ to protect and exchange gas, water, and organic molecules with aquatic and terrestrial environments. Herein, the processes of keratinization and cornification in aquatic and terrestrial vertebrates were compared using immunohistochemistry. Keratins with low cysteine and glycine contents form the main bulk of proteins in the anamniote epidermis, which undergoes keratinization. In contrast, specialized keratins rich in cysteine-glycine and keratin associated corneous proteins rich in cysteine, glycine, and tyrosine form the bulk of proteins of amniote soft cornification in the epidermis and hard cornification in scales, claws, beak, feathers, hairs, and horns. Transglutaminase (TGase) and sulfhydryl oxidase (SOXase) are the main enzymes involved in cornification. Their evolution was fundamental for the terrestrial adaptation of vertebrates. Immunohistochemistry results revealed that TGase and SOXase were low to absent in fish and amphibian epidermis, while they increased in the epidermis of amniotes with the evolution of the stratum corneum and skin appendages. TGase aids the formation of isopeptide bonds, while SOXase forms disulfide bonds that generate numerous cross-links between keratins and associated corneous proteins, likely increasing the mechanical resistance and durability of the amniote epidermis and its appendages. TGase is low to absent in the beta-corneous layers of sauropsids but is detected in the softer but pliable alpha-layers of sauropsids, mammalian epidermis, medulla, and inner root sheath of hairs. SOXase is present in hard and soft corneous appendages of reptiles, birds, and mammals, and determines cross-linking among corneous proteins of scales, claws, beaks, hairs, and feathers.

Immunolocalization of epidermal differentiation complex proteins reveals distinct molecular compositions of cells that control structure and mechanical properties of avian skin appendages

The epidermal differentiation complex (EDC) is a cluster of genes that encode structural proteins of skin derivatives with variable mechanical performances, from the scales of reptiles and birds to the hard claws and beaks, and to the flexible but resistant corneous material of feathers. Corneous proteins with or without extended beta-regions are produced from avian genomes, and include the largely prevalent corneous beta proteins (CβPs, formerly indicated as beta-keratins), and minor contribution from histidine-rich proteins, trichohyalin-like proteins (scaffoldin), loricrin, and other proteins rich in cysteine or other types of amino acids. The light-microscopic and ultrastructural immunolocalization of major and minor EDC-proteins in avian skin (feather CβPs, EDKM, EDWM, EDMTFH, EDDM, and scaffoldin) suggests that each specific appendage consists of a particular mix of these proteins in addition to the main proteins containing a peculiar beta-region of 34 amino acids, indicated as feather/scale/claw/beak CβPs (fCβPs, sCβPs, cCβPs, bCβPs). This indicates that numerous proteins of the EDC are added to the variable meshwork of intermediate filament keratins to produce avian epidermis with different mechanical and functional properties. Although the specific roles for these proteins are not known they likely make an important contribution to the final material properties of the different skin appendages of birds. The highest number of sauropsid CβPs is found in birds, suggesting a relation to the evolution of feathers, and additional epidermal differentiation proteins have contributed to the evolutionary adaptations of avian skin.

Genomic signatures of drift and selection driven by predation and human pressure in an insular lizard

Genomic divergence was studied in 10 small insular populations of the endangered Balearic Islands lizard (Podarcis lilfordi) using double digest restriction-site associated DNA sequencing. The objectives were to establish levels of divergence among populations, investigate the impact of population size on genetic variability and to evaluate the role of different environmental factors on local adaptation. Analyses of 72,846 SNPs supported a highly differentiated genetic structure, being the populations with the lowest population size (Porros, Foradada and Esclatasang islets) the most divergent, indicative of greater genetic drift. Outlier tests identified ~ 2% of loci as candidates for selection. Genomic divergence-Enviroment Association analyses were performed using redundancy analyses based on SNPs putatively under selection, detecting predation and human pressure as the environmental variables with the greatest explanatory power. Geographical distributions of populations and environmental factors appear to be fundamental drivers of divergence. These results support the combined role of genetic drift and divergent selection in shaping the genetic structure of these endemic island lizard populations.

Transcriptome analysis of feather follicles reveals candidate genes and pathways associated with pheomelanin pigmentation in chickens

Yellow plumage is common in chickens, especially in breeds such as the Huiyang Bearded chicken, which is indigenous to China. We evaluated plumage colour distribution in F1, F2, and F3 populations of an Huiyang Bearded chicken × White Leghorn chicken cross, the heredity of the yellow plumage trait was distinguished from that of the gold plumage and other known plumage colours. Microscopic analysis of the feather follicles indicated that pheomelanin particles were formed in yellow but not in white feathers. To screen genes related to formation of the pheomelanin particles, we generated transcriptome data from yellow and white feather follicles from 7- and 11-week-old F3 chickens using RNA-seq. We identified 27 differentially expressed genes (DEGs) when comparing the yellow and white feather follicles. These DEGs were enriched in the Gene Ontology classes ‘melanosome’ and ‘melanosome organization’ related to the pigmentation process. Down-regulation of TYRP1, DCT, PMEL, MLANA, and HPGDS, verified using quantitative reverse transcription PCR, may lead to reduced eumelanin and increased pheomelanin synthesis in yellow plumage. Owing to the presence of the Dominant white locus, both white and yellow plumage lack eumelanin, and white feathers showed no pigments. Our results provide an understanding of yellow plumage formation in chickens.

Transcriptome Comparison Reveals Key Components of Nuptial Plumage Coloration in Crested Ibis

Nuptial plumage coloration is critical in the mating choice of the crested ibis. This species has a characteristic nuptial plumage that develops from the application of a black sticky substance, secreted by a patch of skin in the throat and neck region. We aimed to identify the genes regulating its coloring, by comparing skin transcriptomes between ibises during the breeding and nonbreeding seasons. In breeding season skins, key eumelanin synthesis genes, TYR, DCT, and TYRP1 were upregulated. Tyrosine metabolism, which is closely related to melanin synthesis, was also upregulated, as were transporter proteins belonging to multiple SLC families, which might act during melanosome transportation to keratinocytes. These results indicate that eumelanin is likely an important component of the black substance. In addition, we observed upregulation in lipid metabolism in breeding season skins. We suggest that the lipids contribute to an oil base, which imbues the black substance with water insolubility and enhances its adhesion to feather surfaces. In nonbreeding season skins, we observed upregulation in cell adhesion molecules, which play critical roles in cell interactions. A number of molecules involved in innervation and angiogenesis were upregulated, indicating an ongoing expansion of nerves and blood vessels in sampled skins. Feather β keratin, a basic component of avian feather filament, was also upregulated. These results are consistent with feather regeneration in the black skin of nonbreeding season ibises. Our results provide the first molecular evidence indicating that eumelanin is the key component of ibis coloration.

Characterization of Embryonic Skin Transcriptome in Anser cygnoides at Three Feather Follicles Developmental Stages

In order to enrich the Anser cygnoides genome and identify the gene expression profiles of primary and secondary feather follicles development, de novo transcriptome assembly of skin tissues was established by analyzing three developmental stages at embryonic day 14, 18, and 28 (E14, E18, E28). Sequencing output generated 436,730,608 clean reads from nine libraries and de novo assembled into 56,301 unigenes. There were 2,298, 9,423 and 12,559 unigenes showing differential expression in three stages respectively. Furthermore, differentially expressed genes (DEGs) were functionally classified according to genes ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and series-cluster analysis. Relevant specific GO terms such as epithelium development, regulation of keratinocyte proliferation, morphogenesis of an epithelium were identified. In all, 15,144 DEGs were clustered into eight profiles with distinct expression patterns and 2,424 DEGs were assigned to 198 KEGG pathways. Skin development related pathways (mitogen-activated protein kinase signaling pathway, extra-cellular matrix -receptor interaction, Wingless-type signaling pathway) and genes (delta like canonical Notch ligand 1, fibroblast growth factor 2, Snail family transcriptional repressor 2, bone morphogenetic protein 6, polo like kinase 1) were identified, and eight DEGs were selected to verify the reliability of transcriptome results by real-time quantitative PCR. The findings of this study will provide the key insights into the complicated molecular mechanism and breeding techniques underlying the developmental characteristics of skin and feather follicles in Anser cygnoides.

Mid-Cretaceous amber inclusions reveal morphogenesis of extinct rachis-dominated feathers

We describe three-dimensionally preserved feathers in mid-Cretaceous Burmese amber that share macro-morphological similarities (e.g., proportionally wide rachis with a "medial stripe") with lithic, two-dimensionally preserved rachis-dominated feathers, first recognized in the Jehol Biota. These feathers in amber reveal a unique ventrally concave and dorsoventrally thin rachis, and a dorsal groove (sometimes pigmented) that we identify as the "medial stripe" visible in many rachis-dominated rectrices of Mesozoic birds. The distally pennaceous portion of these feathers shows differentiated proximal and distal barbules, the latter with hooklets forming interlocking barbs. Micro-CT scans and transverse sections demonstrate the absence of histodifferentiated cortex and medullary pith of the rachis and barb rami. The highly differentiated barbules combined with the lack of obvious histodifferentiation of the barb rami or rachis suggests that these feathers could have been formed without the full suite and developmental interplay of intermediate filament alpha keratins and corneous beta-proteins that is employed in the cornification process of modern feathers. This study thus highlights how the development of these feathers might have differed from that of their modern counterparts, namely in the morphogenesis of the ventral components of the rachis and barb rami. We suggest that the concave ventral surface of the rachis of these Cretaceous feathers is not homologous with the ventral groove of modern rachises. Our study of these Burmese feathers also confirms previous claims, based on two-dimensional fossils, that they correspond to an extinct morphotype and it cautions about the common practice of extrapolating developmental aspects (and mechanical attributes) of modern feathers to those of stem birds (and their dinosaurian outgroups) because the latter need not to have developed through identical pathways.

Combined transcriptomics and proteomics forecast analysis for potential genes regulating the Columbian plumage color in chickens

BACKGROUND

Coloration is one of the most recognizable characteristics in chickens, and clarifying the coloration mechanisms will help us understand feather color formation. "Yufen I" is a commercial egg-laying chicken breed in China that was developed by a three-line cross using lines H, N and D. Columbian plumage is a typical feather character of the "Yufen I" H line. To elucidate the molecular mechanism underlying the pigmentation of Columbian plumage, this study utilizes high-throughput sequencing technology to compare the transcriptome and proteome differences in the follicular tissue of different feathers, including the dorsal neck with black and white striped feather follicles (Group A) and the ventral neck with white feather follicles (Group B) in the "Yufen I" H line.

RESULTS

In this study, we identified a total of 21,306 genes and 5,203 proteins in chicken feather follicles. Among these, 209 genes and 382 proteins were differentially expressed in two locations, Group A and Group B, respectively. A total of 8 differentially expressed genes (DEGs) and 9 differentially expressed proteins (DEPs) were found to be involved in the melanogenesis pathway. Additionally, a specifically expressed MED23 gene and a differentially expressed GNAQ protein were involved in melanin synthesis. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis mapped 190 DEGs and 322 DEPs to 175 and 242 pathways, respectively, and there were 166 pathways correlated with both DEGs and DEPs. 49 DEPs/DEGs overlapped and were enriched for 12 pathways. Transcriptomic and proteomic analyses revealed that the following pathways were activated: melanogenesis, cardiomyocyte adrenergic, calcium and cGMP-PKG. The expression of DEGs was validated by real-time quantitative polymerase chain reaction (qRT-PCR) that produced results similar to those from RNA-seq. In addition, we found that the expression of the MED23, FZD10, WNT7B and WNT11 genes peaked at approximately 8 weeks in the "Yufen I" H line, which is consistent with the molting cycle. As both groups showed significant differences in terms of the expression of the studied genes, this work opens up avenues for research in the future to assess their exact function in determining plumage color.

CONCLUSION

Common DEGs and DEPs were enriched in the melanogenesis pathway. MED23 and GNAQ were also reported to play a crucial role in melanin synthesis. In addition, this study is the first to reveal gene and protein variations in in the "Yufen I" H line during Columbian feather color development and to discover principal genes and proteins that will aid in functional genomics studies in the future. The results of the present study provide a significant conceptual basis for the future breeding schemes with the "Yufen I" H line and provide a basis for research on the mechanisms of feather pigmentation.

Genetic architecture related to contour feathers density in an F2 resource population via a genome-wide association study

The density of contour feathers is an important trait as it is closely related to heat dissipation in birds. Thus, identification of the major genes that control this trait will be useful to improve heat tolerance in chicken. So far, no GWAS study for the density of contour feathers in birds has been previously published; therefore, this study was aimed to identify genomic regions controlling the density of contour feathers. A total of 1252 hens were genotyped, using the 600 K Affymetrix Axiom Chicken Genotyping Array. The association analyses were performed using the GenABEL package in the R program. In brief, 146 significant SNP markers were mainly located on chromosome 1 and were identified to associate with the density of contour feathers in the current GWAS analysis. Moreover, we identified several within/nearby candidate genes (SUCLA2, DNAJC15, DHRS12, MLNR, and RB1) that are either directly or indirectly involved in the genetic control of the density of contour feathers in chicken. This study laid the foundation for studying the mechanism that underlies the density of chicken feathers. Furthermore, it is feasible to shear the back feathers of live chickens and measure the density of the feathers to improve heat tolerance in breeding practice.

Stable reference genes for expression studies in breast muscle of normal and white striping-affected chickens

The normalization with proper reference genes is a crucial step to obtain accurate mRNA expression levels in quantitative PCR (qPCR) studies. Therefore, in this study, 10 reference candidate genes were evaluated to determine their stability in normal pectoralis major muscle of broilers and those counterparts affected with White Striping (WS) myopathy at 42 days age. Four different tools were used for ranking the most stable genes: GeNorm, NormFinder, BestKeeper and Comparative Ct (ΔCt), and a general ranking was performed using the RankAggreg tool to select the best reference genes among all tools. From the 10 genes evaluated in the breast muscle of broilers, 8 were amplified. Most of the algorithms/tools indicated the same two genes, RPL30 and RPL5, as the most stable in the broilers breast muscle. In addition, there was agreement among the tools for the least stable genes: MRPS27, GAPDH and RPLP1 in the broilers breast muscle. Therefore, it is interesting to note that even with different tools for evaluating gene expression, there was consensus on the most and least stable genes. These results indicate that the Ribosomal protein L30 (RPL30) and Ribosomal protein L5 (RPL5) can be recommended for accurate normalization in qPCR studies with chicken pectoralis major muscle affected with White Striping and other myopathies.

Feather Evolution from Precocial to Altricial Birds

Birds are the most abundant terrestrial vertebrates and their diversity is greatly shaped by the feathers. How avian evolution is linked to feather evolution has long been a fascinating question. Numerous excellent studies have shed light on this complex relationship by investigating feather diversity and its underlying molecular mechanisms. However, most have focused on adult domestic birds, and the contribution of feather diversity to environmental adaptation has not been well-studied. In this review, we described bird diversity using the traditional concept of the altricial-precocial spectrum in bird hatchlings. We combined the spectrum with a recently published avian phylogeny to profile the spectrum evolution. We then focused on the discrete diagnostic character of the spectrum, the natal down, and propose a hypothesis for the precocial-to-altricial evolution. For the underlying molecular mechanisms in feather diversity and bird evolution, we reviewed the literature and constructed the known mechanisms for feather tract definition and natal down development. Finally, we suggested some future directions for research on altricial-precocial divergence, which may expand our understanding of the relationship between natal down diversity and bird evolution.

Immunolocalization and phylogenetic profiling of the feather protein with the highest cysteine content

Feathers are the most complex skin appendages of vertebrates. Mature feathers consist of interconnected dead keratinocytes that are filled with heavily cross-linked proteins. Although the molecular architecture determines essential functions of feathers, only few feather proteins have been characterized with regard to their amino acid sequences and evolution. Here, we identify Epidermal Differentiation protein containing DPCC Motifs (EDDM) as a cysteine-rich protein that has co-evolved with other feather proteins. The EDDM gene is located within the avian epidermal differentiation complex (EDC), a cluster of genes that has originated and diversified in amniotes. EDDM shares the exon-intron organization with EDC genes of other amniotes, including humans, and a gene encoding an EDDM-like protein is present in crocodilians, suggesting that avian EDDM arose by sequence modification of an epidermal differentiation gene present in a common ancestor of archosaurs. The EDDM protein contains multiple sequence repeats and a higher number of cysteine residues than any other protein encoded in the EDC. Immunohistochemical analysis of chicken skin and skin appendages showed expression of EDDM in barb and barbules of feathers as well as in the subperiderm on embryonic scutate scales. These results suggest that the diversification and differential expression of EDDM, besides other EDC genes, was instrumental in facilitating the evolution of the most complex molecular architecture of feathers.

De Novo Assembly and Comparative Transcriptome Profiling of Anser anser and Anser cygnoides Geese Species' Embryonic Skin Feather Follicles

Geese feather production and the quality of downy feathers are additional economically important traits in the geese industry. However, little information is available about the molecular mechanisms fundamental to feather formation and the quality of feathers in geese. This study conducted de novo transcriptome sequencing analysis of two related geese species using the Illumina 4000 platform to determine the genes involved in embryonic skin feather follicle development. A total of 165,564,278 for Anser anser and 144,595,262 for Anser cygnoides clean reads were generated, which were further assembled into 77,134 unigenes with an average length of 906 base pairs in Anser anser and 66,041 unigenes with an average length of 922 base pairs in Anser cygnoides. To recognize the potential regulatory roles of differentially expressed genes (DEGs) during geese embryonic skin feather follicle development, the obtained unigenes were annotated to Gene Ontology (GO), Eukaryotic Orthologous Groups (KOG), and Kyoto Encyclopedia of Genes and Genomes (KEGG) for functional analysis. In both species, GO and KOG had shown similar distribution patterns during functional annotation except for KEGG, which showed significant variation in signaling enrichment. Anser asnser was significantly enriched in the calcium signaling pathway, whereas Anser cygnoides was significantly enriched with glycerolipid metabolism. Further analysis indicated that 14,227 gene families were conserved between the species, among which a total of 20,715 specific gene families were identified. Comparative RNA-Seq data analysis may reveal inclusive knowledge to assist in the identification of genetic regulators at a molecular level to improve feather quality production in geese and other poultry species.

De Novo Transcriptome Sequencing Analysis of Goose (Anser anser) Embryonic Skin and the Identification of Genes Related to Feather Follicle Morphogenesis at Three Stages of Development

The objective of this study was to evaluate the changes in the goose embryo transcriptome during feather development. RNA-Sequencing (RNA-Seq) was used to find the transcriptome profiles of feather follicles from three stages of embryonic dorsal skin at embryonic day 13, 18, and 28 (E13, E18, E28). The results showed that 3001, 6634, and 13,780 genes were differently expressed in three stages. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that differentially expressed genes (DEGs) in E13 vs. E18 were significantly mapped into the GO term of extracellular structure organization and the pathway of extracellular matrix (ECM)-receptor interaction. In E18 vs. E28, the top significantly mapped into GO term was the single-organism developmental process; the pathway was also the ECM-receptor interaction. DEGs in E13 vs. E28 were significantly mapped into the GO term of the multicellular organismal process and the pathway of cell adhesion molecules. Subsequently, the union of DEGs was categorized by succession cluster into eight profiles, which were then grouped into four ideal profiles. Lastly, the seven genes spatio-temporal expression pattern was confirmed by real-time PCR. Our findings advocate that interleukin 20 receptor subunit alpha (IL20RA), interleukin 6 receptor (IL6R), interleukin 1 receptor type 1 (IL-1R1), Wnt family member 3A (WNT3A), insulin-like growth factor binding protein 3 (IGFBP3), bone morphogenetic protein 7 (BMP7), and secreted-frizzled related protein 2 (SFRP2) might possibly play vital roles in skin and feather follicle development and growth processes.

Genetic and Molecular Basis of Feather Diversity in Birds

Feather diversity is striking in many aspects. Although the development of feather has been studied for decades, genetic and genomic studies of feather diversity have begun only recently. Many questions remain to be answered by multidisciplinary approaches. In this review, we discuss three levels of feather diversity: Feather morphotypes, intraspecific variations, and interspecific variations. We summarize recent studies of feather evolution in terms of genetics, genomics, and developmental biology and provide perspectives for future research. Specifically, this review includes the following topics: 1) Diversity of feather morphotype; 2) feather diversity among different breeds of domesticated birds, including variations in pigmentation pattern, in feather length or regional identity, in feather orientation, in feather distribution, and in feather structure; and 3) diversity of feathers among avian species, including plumage color and morph differences between species and the regulatory differences in downy feather development between altricial and precocial birds. Finally, we discussed future research directions.

Dynamic transcriptome profiling towards understanding the morphogenesis and development of diverse feather in domestic duck

Background

Feathers with complex and fine structure are hallmark avian integument appendages, which have contributed significantly to the survival and breeding for birds. Here, we aimed to explore the differentiation, morphogenesis and development of diverse feathers in the domestic duck.

Results

Transcriptome profiles of skin owing feather follicle from two body parts at three physiological stages were constructed to understand the molecular network and excavate the candidate genes associated with the development of plumulaceous and flight feather structures. The venn analysis of differentially expressed genes (DEGs) between abdomen and wing skin tissues at three developmental stages showed that 38 genes owing identical differentially expression pattern. Together, our data suggest that feather morphological and structural diversity can be possibly related to the homeobox proteins. The key series-clusters, many candidate biological processes and genes were identified for the morphogenesis, growth and development of two feather types. Through comparing the results of developmental transcriptomes from plumulaceous and flight feather, we found that DEGs belonging to the family of WNT, FGF and BMP have certain differences; even the consistent DEGs of skin and feather follicle transcriptomes from abdomen and wing have the different expression patterns.

Conclusions

Overall, this study detected many functional genes and showed differences in the molecular mechanisms of diverse feather developments. The findings in WNT, FGF and BMP, which were consistent with biological experiments, showed more possible complex modulations. A correlative role of HOX genes was also suggested but future biological verification experiments are required. This work provided valuable information for subsequent research on the morphogenesis of feathers.

Electronic supplementary material

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

Disulfide-bond-mediated cross-linking of corneous beta-proteins in lepidosaurian epidermis

Corneous beta-proteins (CBPs), formerly referred to as beta-keratins, are major protein components of the epidermis in lepidosaurian reptiles and are largely responsible for their material properties. These proteins have been suggested to form filaments of 3.4nm in thickness and to interact with themselves or with other proteins, including intermediate filament (IF) keratins. Here, we performed immunocytochemical labeling of CBPs in the epidermis of different lizards and snakes and investigated by immunoblotting analysis whether the reduction of disulfide bonds or protein oxidation affects the solubility and mobility of these CBPs. Immunogold labeling suggested that CBPs are partly co-localized with IF-keratins in differentiating and mature beta-cells. The chemical reduction of epidermal proteins from lizard and snake epidermis increased the abundance of CBP-immunoreactive bands in the size range of CBP monomers on Western blots. Conversely, in vitro oxidation of epidermal proteins reduced the abundance of putative CBP monomers. Some modifications in the IF-keratin range were also noted. These results strongly indicate that CBPs associate with IF-keratins and other proteins via disulfide bonds in the epidermis of lizards and snakes, which likely contributes to the resilience of the cornified beta- and alpha-layers of the lepidosaurian epidermis in live animals and after shedding.

Using scale and feather traits for module construction provides a functional approach to chicken epidermal development

Gene co-expression network analysis has been a research method widely used in systematically exploring gene function and interaction. Using the Weighted Gene Co-expression Network Analysis (WGCNA) approach to construct a gene co-expression network using data from a customized 44K microarray transcriptome of chicken epidermal embryogenesis, we have identified two distinct modules that are highly correlated with scale or feather development traits. Signaling pathways related to feather development were enriched in the traditional KEGG pathway analysis and functional terms relating specifically to embryonic epidermal development were also enriched in the Gene Ontology analysis. Significant enrichment annotations were discovered from customized enrichment tools such as Modular Single-Set Enrichment Test (MSET) and Medical Subject Headings (MeSH). Hub genes in both trait-correlated modules showed strong specific functional enrichment toward epidermal development. Also, regulatory elements, such as transcription factors and miRNAs, were targeted in the significant enrichment result. This work highlights the advantage of this methodology for functional prediction of genes not previously associated with scale- and feather trait-related modules.

Review: cornification, morphogenesis and evolution of feathers

Feathers are corneous microramifications of variable complexity derived from the morphogenesis of barb ridges. Histological and ultrastructural analyses on developing and regenerating feathers clarify the three-dimensional organization of cells in barb ridges. Feather cells derive from folds of the embryonic epithelium of feather germs from which barb/barbule cells and supportive cells organize in a branching structure. The following degeneration of supportive cells allows the separation of barbule cells which are made of corneous beta-proteins and of lower amounts of intermediate filament (IF)(alpha) keratins, histidine-rich proteins, and corneous proteins of the epidermal differentiation complex. The specific protein association gives rise to a corneous material with specific biomechanic properties in barbules, rami, rachis, or calamus. During the evolution of different feather types, a large expansion of the genome coding for corneous feather beta-proteins occurred and formed 3-4-nm-thick filaments through a different mechanism from that of 8-10 nm IF keratins. In the chick, over 130 genes mainly localized in chromosomes 27 and 25 encode feather corneous beta-proteins of 10-12 kDa containing 97-105 amino acids. About 35 genes localized in chromosome 25 code for scale proteins (14-16 kDa made of 122-146 amino acids), claws and beak proteins (14-17 kDa proteins of 134-164 amino acids). Feather morphogenesis is periodically re-activated to produce replacement feathers, and multiple feather types can result from the interactions of epidermal and dermal tissues. The review shows schematic models explaining the translation of the morphogenesis of barb ridges present in the follicle into the three-dimensional shape of the main types of branched or un-branched feathers such as plumulaceous, pennaceous, filoplumes, and bristles. The temporal pattern of formation of barb ridges in different feather types and the molecular control from the dermal papilla through signaling molecules are poorly known. The evolution and diversification of the process of morphogenesis of barb ridges and patterns of their formation within feathers follicle allowed the origin and diversification of numerous types of feathers, including the asymmetric planar feathers for flight.

Animal tracking meets migration genomics: transcriptomic analysis of a partially migratory bird species

Seasonal migration is a widespread phenomenon, which is found in many different lineages of animals. This spectacular behaviour allows animals to avoid seasonally adverse environmental conditions to exploit more favourable habitats. Migration has been intensively studied in birds, which display astonishing variation in migration strategies, thus providing a powerful system for studying the ecological and evolutionary processes that shape migratory behaviour. Despite intensive research, the genetic basis of migration remains largely unknown. Here, we used state-of-the-art radio-tracking technology to characterize the migratory behaviour of a partially migratory population of European blackbirds (Turdus merula) in southern Germany. We compared gene expression of resident and migrant individuals using high-throughput transcriptomics in blood samples. Analyses of sequence variation revealed a nonsignificant genetic structure between blackbirds differing by their migratory phenotype. We detected only four differentially expressed genes between migrants and residents, which might be associated with hyperphagia, moulting and enhanced DNA replication and transcription. The most pronounced changes in gene expression occurred between migratory birds depending on when, in relation to their date of departure, blood was collected. Overall, the differentially expressed genes detected in this analysis may play crucial roles in determining the decision to migrate, or in controlling the physiological processes required for the onset of migration. These results provide new insights into, and testable hypotheses for, the molecular mechanisms controlling the migratory phenotype and its underlying physiological mechanisms in blackbirds and other migratory bird species.

The Feather Model for Chemo- and Radiation Therapy-Induced Tissue Damage

Chemo- and radiation therapy are the main modalities for cancer treatment. A major limiting factor is their toxicity to normal tissue, thus reducing the dose and duration of the therapy. The hair follicle, gastrointestinal tract, and hematopoietic system are among the target organs that often show side effects in cancer therapy . Although these organs are highly mitotic in common, the molecular mechanism of the damage remains unclear. The feather follicle is a fast-growing mini-organ, which allows observation and manipulation on each follicle individually. As a model system, the feather follicle is advantageous because of the following reasons: (1) its complex structure is regulated by a set of evolutionarily conserved molecular pathways, thus facilitating the effort to dissect the specific signaling events involved; (2) its morphology allows the continuity of normal-perturbed-normal structure in a single feather, thus "recording" the damaging effect of chemo- and radiation therapy; (3) further histological and molecular analysis of the damage response can be performed on each plucked feather; thus, it is not necessary to sacrifice the experimental animal. Here, we describe methods of applying the feather model to study the molecular mechanism of chemo- and radiation therapy-induced tissue damage.

Immunolocalization of a Histidine-Rich Epidermal Differentiation Protein in the Chicken Supports the Hypothesis of an Evolutionary Developmental Link between the Embryonic Subperiderm and Feather Barbs and Barbules

The morphogenesis of feathers is a complex process that depends on a tight spatiotemporal regulation of gene expression and assembly of the protein components of mature feathers. Recent comparative genomics and gene transcription studies have indicated that genes within the epidermal differentiation complex (EDC) encode numerous structural proteins of cornifying skin cells in amniotes including birds. Here, we determined the localization of one of these proteins, termed EDMTFH (Epidermal Differentiation Protein starting with a MTF motif and rich in Histidine), which belongs to a group of EDC-encoded proteins rich in aromatic amino acid residues. We raised an antibody against an EDMTFH-specific epitope and performed immunohistochemical investigations by light microscopy and immunogold labeling by electron microscopy of chicken embryos at days 14-18 of development. EDMTFH was specifically present in the subperiderm, a transient layer of the embryonic epidermis, and in barbs and barbules of feathers. In the latter, it partially localized to bundles of so-called feather beta-keratins (corneous beta-proteins, CBPs). Cells of the embryonic periderm, the epidermis proper, and the feather sheath were immunonegative for EDMTFH. The results of this study indicate that EDMTFH may contribute to the unique mechanical properties of feathers and define EDMTFH as a common marker of the subperiderm and the feather barbules. This expression pattern of EDMTFH resembles that of epidermal differentiation cysteine-rich protein (EDCRP) and feather CBPs and is in accordance with the hypothesis that a major part of the cyclically regenerating feather follicle is topologically, developmentally and evolutionarily related to the embryonic subperiderm.