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

Dermal FOXO3 activity in response to Wnt/β-catenin signaling is required for feather follicle development of goose embryos (Anser cygnoides)

Feather is an important economic trait of poultry, and growth and development state of feathers plays an important role in the economic value of poultry. Dermal fibroblasts are required for structural integrity of the skin and for feather follicle development. How FOXO3 affects feather follicle development as skin tissues change during goose embryo (Anser cygnoides) development and growth is not well understood. Here, we demonstrate that in vitro culture of single feathers and skin tissue results in changes in feather morphological structure by adding drugs to the culture medium that affect FOXO3 expression. We used feather follicles to show that during growth, the root location of feathers, the dermis layer, affects cell proliferation and apoptosis and regulates the expression of major genes in the Wingless-types/beta-catenin (Wnt/β-catenin) signaling pathway through the activity of FOXO3 in dermal fibroblasts. Feathers and dorsal skin tissues develop the correct structure, but feather length and width and feather follicle diameter change significantly (p < 0.05) without significant changes in feather follicle density (p > 0.05). Transfected dermal fibroblasts also showed that FOXO3 affected the formation and development of feather follicles in the embryonic stage by regulating the Wnt/β-catenin signaling pathway. Therefore, this study reveals the critical role of dermal fibroblast-FOXO3-induced Wnt/β-catenin signaling in promoting the formation and development of embryonic feather follicles.

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.

Expression of FOXO3 in the skin follicles of goose embryos during embryonic development

1. The Forkhead box O3 (FOXO3) transcription factor is a crucial regulator in controlling cell metabolism, proliferation, apoptosis, migration and response to oxidative stress. However, FOXO3 has not previously been studied much in the embryonic skin follicles of geese.2. This study used Zhedong white geese (Anser cygnoides), Jilin white geese (Anser cygnoides) and Hungarian white geese (Anser anser). The feather follicle structure in the dorsal skin during embryonic stages was examined with haematoxylin and eosin (HE) and Pollak staining. The FOXO3 protein content in the embryonic dorsal skin from feather follicles was detected using western blotting and quantitative real-time PCR.3. The mRNA expression level of FOXO3 in the dorsal skin of Jilin white geese was highly expressed on embryonic day 23 (E23; P < 0.01), while mRNA expression of FOXO3 was highly expressed in the feather follicle of Hungarian white geese at E28 (P < 0.01). The expression of FOXO3 protein mainly concentrated in the early embryonic phase among these goose breeds (P < 0.05). This suggested that FOXO3 plays a crucial role in the development and growth of embryonic dorsal skin of feather follicles. The location of the FOXO3 protein was determined using the IHC technique, which further verified the effect of FOXO3 in the dorsal skin for feather follicles during embryogenesis.4. The study demonstrated the differential expression and localisation of the FOXO3 gene among different goose species. It was speculated that the gene could potentially improve goose feather follicle development and feather-related traits and provide a basis for further understanding of FOXO3 function in the dorsal tissue of goose embryos.

Selection signatures for fiber production in commercial species: A review

Natural fibers derived from diverse animal species have gained increased attention in recent years due to their favorable environmental effects, long-term sustainability benefits, and remarkable physical and mechanical properties that make them valuable raw materials used for textile and non-textile production. Domestication and selective breeding for the economically significant fiber traits play an imperative role in shaping the genomes and, thus, positively impact the overall productivity of the various fiber-producing species. These selection pressures leave unique footprints on the genome due to alteration in the allelic frequencies at specific loci, characterizing selective sweeps. Recent advances in genomics have enabled the discovery of selection signatures across the genome using a variety of methods. The increased demand for 'green products' manufactured from natural fibers necessitates a detailed investigation of the genomes of the various fiber-producing plant and animal species to identify the candidate genes associated with important fiber attributes such as fiber diameter/fineness, color, length, and strength, among others. The objective of this review is to present a comprehensive overview of the concept of selection signature and selective sweeps, discuss the main methods used for its detection, and address the selection signature studies conducted so far in the diverse fiber-producing animal species.

Nonlinear model fitting analysis of feather growth and development curves in the embryonic stages of Jilin white geese (Anser cygnoides)

Poultry is subject to varying degrees of feather loss and feather pecking during production, which seriously affects the live appearance and carcass appearance of their commercial traits and greatly reduces the production profitability of the farming enterprise. It also has an impact on down production and quality in the case of geese. In this study, mathematical models (Logistic, Gompertz, and Von Bertalanffy) were used to assess feather growth and development during the embryonic period in Jilin white geese (Anser cygnoides) predicting the weight and length of feathers from the back, chest, and belly tracts at different embryonic ages, to determine which growth model more accurately described feather growth patterns. The result first showed that the primary feather follicles of the Jilin white goose developed at E14 and secondary feather follicles at E18; primary feather follicle density increased and then decreased, whereas secondary feather follicle density increased continuously and the primary and secondary feather follicles developed independently. Secondly, the embryonic feather growth followed a slow-fast-slow pattern, with feathers growing slowly from E12 to E18, quickly from E18 to E24, and then decreasing after E24 until just before emergence (E30). In addition, before E14, feathers were concentrated in the back tracts, and no feathers were found on the head, neck, chest, abdomen, or wings. By E22, the whole body of the embryo was covered with feathers, and the back feathers were the earliest and fastest to develop. Compared to the Gompertz and von Bertalanffy models, the logistic model fit (R2 = 0.997) was the highest, while the sum of residual squares (RSS = 25661.67), Akaike's information criterion (AIC = 77.600), Bayesian information criterion (BIC = 78.191), and mean square error (MSE = 2851.296) were the lowest. Therefore, the logistic model was more suitable for describing the changes in whole-body feather growth during the embryonic period in Jilin white geese. In conclusion, using the growth curve model to explain the relationship between feather growth and embryonic age in geese will potentially speed up the process of genetic improvement in Jilin white geese (A. cygnoides) and thus provide scientific support for molecular genetic breeding.

Regulation of feather follicle development and Msx2 gene SNP degradation in Hungarian white goose

BACKGROUND

Hungarian white goose has excellent down production performance and was introduced to China in 2010. The growth and development of feather follicles has an important impact on down production. Goose feather follicles can be divided into primary and secondary feather follicles, both of which originate in the embryonic stage. Msx2 (Msh Homeobox 2) plays a regulatory role in tissues and organs such as eyes, teeth, bones and skin. However, its regulatory mechanism on goose feather follicles development remains unclear.

RESULTS

Msx2 gene first increased, then decreased and increased at the end (E13, E18, E23, E28) during embryonic feather follicle development, and the expression level was the highest at E18. The pEGFP-N1-Msx2 overexpression vector and si-Msx2 siRNA vector were constructed to transfect goose embryo dermal fibroblasts. The results showed that the cell viability of ov-Msx2 group was significantly increased, and the gene expression levels of FGF5 and TGF-β1 genes were significantly down-regulated (P < 0.05), the expressions of PCNA, Bcl2, CDK1, FOXN1 and KGF genes were significantly up-regulated (P < 0.05). After transfection of siRNA vector, the cell viability of the si-Msx2 group was significantly decreased (P < 0.01) compared with the si-NC group. TGF-β1 expression was significantly up-regulated (P < 0.05), FGF5 expression was extremely significantly up-regulated (P < 0.01), while PCNA, Bcl2, CDK1, FOXN1 and KGF gene expression was significantly down-regulated (P < 0.05). High-throughput sequencing technology was used to mine the exon SNPs of Msx2. A total of 11 SNP loci were screened, four of the SNPs located in exon 1 were missense mutations. The feather follicle diameter of the GC genotype at the G78C site is significantly larger than that of the other two genotypes.

CONCLUSIONS

Msx2 maybe inhibit the apoptosis of goose dermal fibroblasts and promotes their proliferation. G78C can be used as a potential molecular marker for downy Variety.

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.

Using comparative genomics to detect mutations regulating plumage variations in graylag (A. anser) and swan geese (A. cygnoides)

Although graylag geese (A. anser) showed similar plumages of white, grey, and white with grey patches compared to those in swan geese (A. cygnoides), it was believed the substantial molecular mechanism for plumage variations were different. To date, studies on genes responsible for diverse plumages among graylag geese were limited and causal mutations remain unknown. In this study, genomes from 57 individuals belonging to six breeds showing different plumages were sequenced at ∼10X depth. Firstly, the allele frequency differences (AFD) of variants on the scaffold394 (NW_013185915.1) between grey and white goose breeds (A. anser) was calculated and a genomic region between 768,290-779,889 bp was detected to carry candidate variants associated with plumages, including one SNP (g. 775,151G > T, ∼18.6 kb upstream of EDNRB2) found to be fixed in white geese. This region was overlapped with the one detected by the haplotype-based sweep analysis, in which significant signals defined a candidate region of 736,610-820,622 bp on the same scaffold. Results from the transcriptomic data showed that expression levels of EDNRB2 and many other melanogenesis-related genes were significantly decreased among white geese compared to that in grey geese, especially at late embryonic stages (>E15). Modifications at transcriptional levels might result in abnormal melanocyte developments and thus the white plumages when they grow up. In addition, a frameshift mutation (C > -) in exon4 of MLANA gene on scaffold176 (NW_013185876.1) was suggested as the causal mutation for sex-linked dilution phenotype in graylag geese although this requires more demonstration experiments. Together with observed white plumages caused by EDNRB2 mutations in coding regions among swan geese and chicken, our study provided new examples to study the parallel evolution.

In Ovo Injection of CHIR-99021 Promotes Feather Follicle Development via Modulating the Wnt Signaling Pathway and Transcriptome in Goose Embryos (Anser cygnoides)

Feather performs important physiological functions in birds, and it is also one of the economic productions in goose farming. Understanding and modulating feather follicle development during embryogenesis are essential for bird biology and the poultry industry. CHIR-99021 is a potent Wnt/β-catenin signaling pathway activator associated with feather follicle development. In this study, goose embryos (Anser cygnoides) received an in ovo injection of CHIR-9902, which was conducted at the beginning of feather follicle development (E9). The results showed that feather growth and feather follicle development were promoted. The Wnt signaling pathway was activated by the inhibition of GSK-3β. Transcriptomic analyses showed that the transcription changes were related to translation, metabolism, energy transport, and stress in dorsal tissue of embryos that received CHIR-99021, which might be to adapt and coordinate the promoting effects of CHIR-99021 on feather follicle development. This study suggests that in ovo injection of CHIR-99021 is a potential strategy to improve feather follicle development and feather-related traits for goose farming and provides profiling of the Wnt signaling pathway and transcriptome in dorsal tissue of goose embryos for further understanding of feather follicle development.

Identifying molecular pathways and candidate genes associated with knob traits by transcriptome analysis in the goose (Anser cygnoides)

Anser cygnoides has a spherical crest on the beak roof, which is described as knob. However, the mechanisms affecting knob morphology are unclear. Here, we investigated the phenotypic characteristics and molecular basis of knob-size differences in Yangzhou geese. Anatomically, the knob was identified as frontal hump in the frontal area of the skull, rather than hump of upper beak. Although the frontal hump length, and height varied greatly in geese with different knob phenotypes, little was changed in the width. Histologically, knob skin in large-size knobs geese have a greater length in the stratum corneum, stratum spinosum, and stratum reticular than that in small-size knobs geese. Moveover, the 415 differentially expressed genes were found between the large knobs and small ones through transcriptome profiling. In addition, GO enrichment and KEGG pathway analysis revealed 455 significant GO terms and 210 KEGG pathways were enriched, respectively. Among these, TGF-β signaling and thyroid hormone synthesis-signaling pathways were identified to determine knob-size phenotype. Furthermore, BMP5, DCN, TSHR and ADCY3 were recognized to involve in the growth and development of knob. Our data provide comprehensive molecular determinants of knob size phenotype, which can potentially promote the genetic improvement of goose knobs.

Genome-wide association studies detects candidate genes for wool traits by re-sequencing in Chinese fine-wool sheep

BACKGROUND

The quality and yield of wool determine the economic value of the fine-wool sheep. Therefore, discovering markers or genes relevant to wool traits is the cornerstone for the breeding of fine-wool sheep. In this study, we used the Illumina HiSeq X Ten platform to re-sequence 460 sheep belonging to four different fine-wool sheep breeds, namely, Alpine Merino sheep (AMS), Chinese Merino sheep (CMS), Aohan fine-wool sheep (AHS) and Qinghai fine-wool sheep (QHS). Eight wool traits, including fiber diameter (FD), fiber diameter coefficient of variance (FDCV), fiber diameter standard deviation (FDSD), staple length (SL), greasy fleece weight (GFW), clean wool rate (CWR), staple strength (SS) and staple elongation (SE) were examined. A genome-wide association study (GWAS) was performed to detect the candidate genes for the eight wool traits.

RESULTS

A total of 8.222 Tb of raw data was generated, with an average of approximately 8.59X sequencing depth. After quality control, 12,561,225 SNPs were available for analysis. And a total of 57 genome-wide significant SNPs and 30 candidate genes were detected for the desired wool traits. Among them, 7 SNPs and 6 genes are related to wool fineness indicators (FD, FDCV and FDSD), 10 SNPs and 7 genes are related to staple length, 13 SNPs and 7 genes are related to wool production indicators (GFW and CWR), 27 SNPs and 10 genes associated with staple elongation. Among these candidate genes, UBE2E3 and RHPN2 associated with fiber diameter, were found to play an important role in keratinocyte differentiation and cell proliferation. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment results, revealed that multitude significant pathways are related to keratin and cell proliferation and differentiation, such as positive regulation of canonical Wnt signaling pathway (GO:0090263).

CONCLUSION

This is the first GWAS on the wool traits by using re-sequencing data in Chinese fine-wool sheep. The newly detected significant SNPs in this study can be used in genome-selective breeding for the fine-wool sheep. And the new candidate genes would provide a good theoretical basis for the fine-wool sheep breeding.

Identification of key genes and pathways associated with duck (Anas platyrhynchos) embryonic skin development using weighted gene co-expression network analysis

Skin and feather follicle morphogenesis are important processes for duck development; however, the mechanisms underlying morphogenesis at the embryonic stage remain unclear. To improve the understanding of these processes, we used transcriptome and weighted gene co-expression network analyses to identify the critical genes and pathways involved in duck skin development. Five modules were found to be the most related to five key stages in skin development that span from embryonic day 8 (E8) to postnatal day 7 (D7). Using STEM software, 6519 genes from five modules were clustered into 10 profiles to reveal key genes. Above all, we obtained several key module genes including WNT3A, NOTCH1, SHH, BMP2, NOG, SMAD3, and TGFβ2. Furthermore, we revealed that several pathways play critical roles throughout the skin development process, including the Wnt pathway and cytoskeletal rearrangement-related pathways, whereas others are involved in specific stages of skin development, such as the Notch, Hedgehog, and TGF-beta signaling pathways. Overall, this study identified the pathways and genes that play critical roles in skin development, which may provide a basis for high-quality down-type meat duck breeding.

The role of CTNNB1 and LEF1 in feather follicles development of Anser cygnoides and Anser anser

BACKGROUND

Wingless-types/beta-catenin (Wnt/β-catenin) signaling pathway is one of the most extensively studied transcriptional cascades involved in various types of organogenesis including embryonic and postnatal development. Downy feather quantity is primarily affected by follicular development and gene regulations.

OBJECTIVE

This research was aimed to investigate the role of catenin beta-1(CTNNB1) and lymphoid enhancerbinding factor-1 (LEF1) on feather follicles development at different developmental stages.

METHODS

Fluorescence quantitative PCR, Western-blot and immunohistochemical methods were used in Anser cygnoides and Anser anser embryos (E12, E13 E18, and E28) and after birth gosling stages (G18, G48, G88) for gene expression analysis.

RESULTS

CTNNB1 and LEF1 genes were expressed in Anser cygnoides and Anser anser at different embryonic and after-birth gosling developmental stages and the expression levels were significantly different in different stages (p < 0.05). The mRNA expression of CTNNB1 and LEF1 genes reached the highest level at D88 in Anser cygnoides, while the highest expression levels were at D18 and D88 in Anser anser, and the expression levels of CTNNB1 genes at D88 in all embryonic stages were significantly lower than after-birth stages. CTNNB1 and LEF1 protein expression were the highest at E12 and E28 for Anser cygnoides feather follicles development. While at a similar stage for Anser anser, the expression of CTNNB1 and LEF1 protein was the highest at D48 and D18. Protein expression at embryonic stages was in the epidermis (E) and the hair basal plate (P), the expression site for after-birth stages was in the dermal papilla (DP).

CONCLUSION

Our study illustrated that CTNNB1 and LEF1 has an impact on Anser cygnoides and Anser anser feather follicles growth and development.

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.

Exploration of key regulators driving primary feather follicle induction in goose skin

The primary feather follicles are universal skin appendages widely distributed in the skin of feathered birds. The morphogenesis and development of the primary feather follicles in goose skin remain largely unknown. Here, the induction of primary feather follicles in goose embryonic skin (pre-induction vs induction) was investigated by de novo transcriptome analyses to reveal 409 differentially expressed genes (DEGs). The DEGs were characterized to potentially regulate the de novo formation of feather follicle primordia consisting of placode (4 genes) and dermal condensate (12 genes), and the thickening of epidermis (5 genes) and dermal fibroblasts (17 genes), respectively. Further analyses enriched DEGs into GO terms represented as cell adhesion and KEGG pathways including Wnt and Hedgehog signaling pathways that are highly correlated with cell communication and molecular regulation. Six selected Wnt pathway genes were detected by qPCR with up-regulation in goose skin during the induction of primary feather follicles. The localization of WNT16, SFRP1 and FRZB by in situ hybridization showed weak expression in the primary feather primordia, whereas FZD1, LEF1 and DKK1 were expressed initially in the inter-follicular skin and feather follicle primordia, then mainly restricted in the feather primordia. The spatial-temporal expression patterns indicate that Wnt pathway genes DKK1, FZD1 and LEF1 are the important regulators functioned in the induction of primary feather follicle in goose skin. The dynamic molecular changes and specific gene expression patterns revealed in this report provide the general knowledge of primary feather follicle and skin development in waterfowl, and contribute to further understand the diversity of hair and feather development beyond the mouse and chicken models.

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.

Transcriptome analysis of leg muscles in fast and slow growth Bian chickens

Chicken is popular among consumers in the market, but the mechanism for regulating its growth is still unclear. In this experiment, two groups of Bian chickens of different body weights at 16 weeks of age were studied. The leg muscles were taken for transcriptome sequencing after slaughter. In the differential gene screening, all the genes obtained by sequencing the fast and slow growth groups were screened by Fold Change ≥2 and False Discovery Rate (FDR) <0.05, and 108 differentially expressed genes were obtained. The slow growth group has 17 up-regulated genes and 91 down-regulated genes compared with the fast growing group. Significance analysis of differentially expressed genes in gene ontology (GO) enrichment indicates that there are 65, 16 and 6 significantly enriched entries in the three main categories of biological processes, cellular components and molecular functions (P-value <0.05), respectively. Pathway enrichment analysis yielded three significantly enriched signal pathways: Adrenergic signaling in cardiomyocytes, Cardiac muscle contraction and Tight junction. The experiment would contribute to reveal the molecular mechanism of chicken growth and provide a theoretical basis for improving the performance of Bian chicken.