Identification of genes associated with sperm storage capacity in hens at different times after insemination by RNA-seq and Ribo-seq.

BACKGROUND: Sperm storage capacity (SSC) determines the duration of fertility in hens and is an important reproduction trait that cannot be ignored in production. Currently, the genetic mechanism of SSC is still unclear in hens. Therefore, to explore the genetic basis of SSC, we analyzed the uterus-vagina junction (UVJ) of hens with different SSC at different times after insemination by RNA-seq and Ribo-seq. RESULTS: Our results showed that 589, 596, and 527 differentially expressed genes (DEGs), 730, 783, and 324 differentially translated genes (DTGs), and 804, 625, and 467 differential translation efficiency genes (DTEGs) were detected on the 5th, 10th, and 15th days after insemination, respectively. In transcription levels, we found that the differences of SSC at different times after insemination were mainly reflected in the transmission of information between cells, the composition of intercellular adhesion complexes, the regulation of ion channels, the regulation of cellular physiological activities, the composition of cells, and the composition of cell membranes. In translation efficiency (TE) levels, the differences of SSC were mainly related to the physiological and metabolic activities in the cell, the composition of the organelle membrane, the physiological activities of oxidation, cell components, and cell growth processes. According to pathway analysis, SSC was related to neuroactive ligand-receptor interaction, histidine metabolism, and PPAR signaling pathway at the transcriptional level and glutathione metabolism, oxidative phosphorylation, calcium signaling pathway, cell adhesion molecules, galactose metabolism, and Wnt signaling pathway at the TE level. We screened candidate genes affecting SSC at transcriptional levels (COL4A4, MUC6, MCHR2, TACR1, AVPR1A, COL1A1, HK2, RB1, VIPR2, HMGCS2) and TE levels(COL4A4, MUC6, CYCS, NDUFA13, CYTB, RRM2, CAMK4, HRH2, LCT, GCK, GALT). Among them, COL4A4 and MUC6 were the key candidate genes differing in transcription, translation, and translation efficiency. CONCLUSIONS: Our study used the combined analysis of RNA-seq and Ribo-seq for the first time to investigate the SSC and reveal the physiological processes associated with SSC. The key candidate genes affecting SSC were screened, and the theoretical basis was provided for the analysis of the molecular regulation mechanism of SSC.

Molecular phylogeny and morphology of Sporodiniella sinensis sp. nov. (Syzygitaceae, Mucorales), an invertebrate-associated species from Yunnan, China.

During investigations of invertebrate-associated fungi in Yunnan Province of China, a new species, Sporodiniella sinensis sp. nov., was collected. Morphologically, S. sinensis is similar to Sporodiniella umbellata; however, it is distinguished from S. umbellata by its greater number of sporangiophore branches, longer sporangiophores, larger sporangiospores, and columellae. The novel species exhibits similarities of 91.62 % for internal transcribed spacer (ITS), 98.66-99.10 % for ribosomal small subunit (nrSSU), and 96.36-98.22 % for ribosomal large subunit (nrLSU) sequences, respectively, compared to S. umbellata. Furthermore, phylogenetic analyses based on combined sequences of ITS, nrLSU and nrSSU show that it forms a separate clade in Sporodiniella, and clusters closely with S. umbellata with high statistical support. The phylogenetic and morphological evidence support S. sinensis as a distinct species. Here, it is formally described and illustrated, and compared with other relatives.

Transcriptome profiling analysis of uterus during chicken laying periods.

The avian eggshell is formed in the uterus. Changes in uterine function may have a significant effect on eggshell quality. To identify the vital genes impacting uterine functional maintenance in the chicken, uteri in three different periods (22W, 31W, 51W) were selected for RNA sequencing and bioinformatics analysis. In our study, 520, 706 and 736 differentially expressed genes (DEGs) were respectively detected in the W31 vs W22 group, W51 vs W31 group and W51 vs W22 group. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated DEGs were enriched in the extracellular matrix, extracellular region part, extracellular region, extracellular matrix structural constituent, ECM receptor interaction, collagen-containing extracellular matrix and collagen trimer in the uterus (P < 0.05). Protein-protein interaction analysis revealed that FN1, LOX, THBS2, COL1A1, COL1A2, COL5A1, COL5A2, POSTN, MMP13, VANGL2, RAD54B, SPP1, SDC1, BTC, ANGPTL3 might be key candidate genes for uterine functional maintenance in chicken. This study discovered dominant genes and pathways which enhanced our knowledge of chicken uterine functional maintenance.

Grade follicles transcriptional profiling analysis in different laying stages in chicken.

During follicular development, a series of key events such as follicular recruitment and selection are crucially governed by strict complex regulation. However, its molecular mechanisms remain obscure. To identify the dominant genes controlling chicken follicular development, the small white follicle (SWF), the small yellow follicle (SYF), and the large yellow follicle (LYF) in different laying stages (W22, W31, W51) were collected for RNA sequencing and bioinformatics analysis. There were 1866, 1211, and 1515 differentially expressed genes (DEGs) between SWF and SYF in W22, W31, and W51, respectively. 4021, 2295, and 2902 DEGs were respectively identified between SYF and LYF in W22, W31, and W51. 5618, 4016, and 4809 DEGs were respectively identified between SWF and LYF in W22, W31, and W51. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated that extracellular matrix, extracellular region, extracellular region part, ECM-receptor interaction, collagen extracellular matrix, and collagen trimer were significantly enriched (P < 0.05). Protein-protein interaction analysis revealed that COL4A2, COL1A2, COL4A1, COL5A2, COL12A1, ELN, ALB, and MMP10 might be key candidate genes for follicular development in chicken. The current study identified dominant genes and pathways contributing to our understanding of chicken follicular development.

Evolutionary innovations through gain and loss of genes in the ectomycorrhizal Boletales.

We aimed to identify genomic traits of transitions to ectomycorrhizal ecology within the Boletales by comparing the genomes of 21 symbiotrophic species with their saprotrophic brown-rot relatives. Gene duplication rate is constant along the backbone of Boletales phylogeny with large loss events in several lineages, while gene family expansion sharply increased in the late Miocene, mostly in the Boletaceae. Ectomycorrhizal Boletales have a reduced set of plant cell-wall-degrading enzymes (PCWDEs) compared with their brown-rot relatives. However, the various lineages retain distinct sets of PCWDEs, suggesting that, over their evolutionary history, symbiotic Boletales have become functionally diverse. A smaller PCWDE repertoire was found in Sclerodermatineae. The gene repertoire of several lignocellulose oxidoreductases (e.g. laccases) is similar in brown-rot and ectomycorrhizal species, suggesting that symbiotic Boletales are capable of mild lignocellulose decomposition. Transposable element (TE) proliferation contributed to the higher evolutionary rate of genes encoding effector-like small secreted proteins, proteases, and lipases. On the other hand, we showed that the loss of secreted CAZymes was not related to TE activity but to DNA decay. This study provides novel insights on our understanding of the mechanisms influencing the evolutionary diversification of symbiotic boletes.

Transcriptome landscapes of differentially expressed genes related to fat deposits in Nandan-Yao chicken.

Nandan-Yao chicken is a Chinese native chicken with lower fat deposition and better meat quality. Fat deposition is a quite complex and important economic trait. However, its molecular mechanism is still unknown in chickens. In the current study, Nandan-Yao chicken was divided into two groups based on the rate of abdominal fat at 120 days old, namely the high-fat group and low-fat group. The total RNAs were isolated and sequenced by RNA sequencing (RNA-seq). After quality control, we gained 1222, 902, 784, 624, and 736 differentially expressed genes (DEGs) in abdominal fat, back skin, liver, pectoral muscle, and leg muscle, respectively. Analysis of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) showed that significantly enriched GO term and KEGG signaling pathway mainly involved cytosolic ribosome, growth development, PPAR signaling pathway, Wnt signaling pathway, and linoleic acid metabolism in abdominal fat, back skin, and liver. While in pectoral muscle and leg muscle, it is mainly enriched in phosphatidylinositol signaling system, adrenergic signaling in cardiomyocytes, cytosolic ribosome, and cytosolic part. Sixteen genes were differentially expressed in all five tissues. Among them, PLA2G4A and RPS4Y1 might be the key regulators for fat deposition in Nandan-Yao chicken. The protein-protein interaction (PPI) network analysis of DEGs showed that PCK1 was the most notable genes. The findings in the current study will help to understand the regulation mechanism of abdominal fat and intramuscular fat in Nandan-Yao chicken and provide a theoretical basis for Chinese local chicken breeding.

Genome-wide analysis and prediction of genes involved in the biosynthesis of polysaccharides and bioactive secondary metabolites in high-temperature-tolerant wild Flammulina filiformis.

BACKGROUND: Flammulina filiformis (previously known as Asian F. velutipes) is a popular commercial edible mushroom. Many bioactive compounds with medicinal effects, such as polysaccharides and sesquiterpenoids, have been isolated and identified from F. filiformis, but their biosynthesis and regulation at the molecular level remains unclear. In this study, we sequenced the genome of the wild strain F. filiformis Liu355, predicted its biosynthetic gene clusters (BGCs) and profiled the expression of these genes in wild and cultivar strains and in different developmental stages of the wild F. filiformis strain by a comparative transcriptomic analysis. RESULTS: We found that the genome of the F. filiformis was 35.01 Mb in length and harbored 10,396 gene models. Thirteen putative terpenoid gene clusters were predicted and 12 sesquiterpene synthase genes belonging to four different groups and two type I polyketide synthase gene clusters were identified in the F. filiformis genome. The number of genes related to terpenoid biosynthesis was higher in the wild strain (119 genes) than in the cultivar strain (81 genes). Most terpenoid biosynthesis genes were upregulated in the primordium and fruiting body of the wild strain, while the polyketide synthase genes were generally upregulated in the mycelium of the wild strain. Moreover, genes encoding UDP-glucose pyrophosphorylase and UDP-glucose dehydrogenase, which are involved in polysaccharide biosynthesis, had relatively high transcript levels both in the mycelium and fruiting body of the wild F. filiformis strain. CONCLUSIONS: F. filiformis is enriched in a number of gene clusters involved in the biosynthesis of polysaccharides and terpenoid bioactive compounds and these genes usually display differential expression between wild and cultivar strains, even in different developmental stages. This study expands our knowledge of the biology of F. filiformis and provides valuable data for elucidating the regulation of secondary metabolites in this unique F. filiformis strain.

Spiroseoflosterol, a Rearranged Ergostane-Steroid from the Fruiting Bodies of Butyriboletus roseoflavus.

Spiroseoflosterol (1), a 7(8→9)-abeo-ergostane steroid with a unique spiro[4.5]decan-6-one system, was isolated from the fruiting bodies of Butyriboletus roseoflavus. Its structure was determined by interpretation of comprehensive spectroscopic, X-ray, and computational data. A plausible biogenetic pathway for spiroseoflosterol (1) was postulated based on a key semipinacol rearrangement. Compound 1 was cytotoxic to HepG2 and Huh7/S (sorafenib-resistant Huh7) with IC50 values of 9.1 and 6.2 µM, respectively.

Bioactive Sesquiterpenes from the Edible Mushroom Flammulina velutipes and Their Biosynthetic Pathway Confirmed by Genome Analysis and Chemical Evidence.

Twelve putative sesquiterpene synthases genes were found in clades along with enzymes with 1,6-, 1,10-, and 1,11-cyclase activities in the genome of Flammulina velutipes. Chemistry investigation of F. velutipes led to the identification of two seco-cuparane sesquiterpenes, flammufuranone A (1) and B (2); 13 new sesquiterpenes with nor-eudesmane, spiroaxane, cadinane, and cuparane skeletons (3-14, 16); as well as two new ergosterol derivatives (17 and 18). Sesquiterpenes (3-14) derived from 1,10-cyclizing enzyme were first reported from this mushroom. The absolute configurations in 1 (3R,7S) and 2 (3R,7R) were assigned by electronic circular dichroism (ECD) calculation. The absolute configuration in 3 was confirmed by X-ray diffraction analysis. The absolute configurations in the 1,2-diol moiety of 13, and in the 1,3-diol moiety of 17 and 18 were determined using Snatzke's method. Among these compounds, 3, 5, 13, and 14 were found to inhibit the HMG-CoA reductase with IC50 of 114.7, 77.6, 55.5, and 87.1 µM, respectively. Compounds 5, 6, 7, 10, 13, and 14 showed DPP-4 inhibitory activity with IC50 of 75.9, 83.7, 70.9, 79.7, 80.5, and 74.8 µM, respectively. The biosynthesis for sesquiterpenes in F. velutipes was also discussed.

Species clarification of the culinary Bachu mushroom in western China.

The Bachu mushroom, previously identified as Helvella leucopus, is characterized by a saddle-shaped, to irregularly lobed pileus, with a gray, brown to blackish hymenium and a whitish to pale receptacle surface and white, terete stipe with enlarged basal grooves. It has high economic value, mostly as a dietary supplement in western China, and its medicinal functions have raised broad interest. In the present paper species of the Bachu mushroom in Xinjiang Autonomous Region, western China were investigated with morphology and DNA sequence data. Phylogenetic analyses inferred from ITS, 28S and TEF1 sequence data strongly supported lineages corresponding to morphological features. The Bachu mushroom, which differs from the European Helvella leucopus, comprises two distinct new species, namely Helvella bachu and Helvella subspadicea. In this paper we introduce the new species with descriptions and figures and compare them with similar taxa. The European Helvella spadicea is also re-examined, described and illustrated.

Antiviral effects of two Ganoderma lucidum triterpenoids against enterovirus 71 infection.

Enterovirus 71 (EV71) is a major causative agent for hand, foot and mouth disease (HFMD), and fatal neurological and systemic complications in children. However, there is currently no clinical approved antiviral drug available for the prevention and treatment of the viral infection. Here, we evaluated the antiviral activities of two Ganoderma lucidum triterpenoids (GLTs), Lanosta-7,9(11),24-trien-3-one,15;26-dihydroxy (GLTA) and Ganoderic acid Y (GLTB), against EV71 infection. The results showed that the two natural compounds display significant anti-EV71 activities without cytotoxicity in human rhabdomyosarcoma (RD) cells as evaluated by 3-(4,5-Dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) cell proliferation assay. The mechanisms by which the two compounds affect EV71 infection were further elucidated by three action modes using Ribavirin, a common antiviral drug, as a positive control. The results suggested that GLTA and GLTB prevent EV71 infection through interacting with the viral particle to block the adsorption of virus to the cells. In addition, the interactions between EV71 virion and the compounds were predicated by computer molecular docking, which illustrated that GLTA and GLTB may bind to the viral capsid protein at a hydrophobic pocket (F site), and thus may block uncoating of EV71. Moreover, we demonstrated that GLTA and GLTB significantly inhibit the replication of the viral RNA (vRNA) of EV71 replication through blocking EV71 uncoating. Thus, GLTA and GLTB may represent two potential therapeutic agents to control and treat EV71 infection.

Molecular phylogenetics of the Ophiocordyceps sinensis-species complex lineage (Ascomycota, Hypocreales), with the discovery of new species and predictions of species distribution.

Ophiocordyceps sinensis is a famous traditional Chinese medicine adapted to the alpine environment of the Qinghai-Tibet Plateau and adjacent regions. Clarification of the species diversity of Ophiocordyceps sinensis and its relatives could expand the traditional medicinal resources and provide insights into the speciation and adaptation. The study is prompted by the discovery of a new species, O. megala, described here from a biodiversity hotspot in the Hengduan Mountains, China. Combined morphological, ecological, and phylogenetic evidence supports its distinctiveness from O. sinensis, O. xuefengensis, and O. macroacicularis. Additionally, based on the phylogenetic construction of Ophiocordyceps, a special clade was focused phylogenetically on the more closely related O. sinensis complex, which was defined as the O. sinensis- species complex lineage. A total of 10 species were currently confirmed in this lineage. We made a comprehensive comparison of the sexual/asexual morphological structures among this species complex, distinguishing their common and distinctive features. Furthermore, using the method of species distribution modelling, we studied the species ocurrences in relation to climatic, edaphic, and altitudinal variables for the eight species in the O. sinensis-species complex, and determined that their potential distribution could extend from the southeastern Qinghai-Tibet Plateau to the Xuefeng Mountains without isolating barrier. Thus, the biodiversity corridor hypothesis was proposed around the O. sinensis-species complex. Our study highlights the phylogeny, species diversity, and suitable distribution of the O. sinensis-species complex lineage, which should have a positive implication for the resource discovery and adaptive evolution of this unique and valuable group.

The sperm storage capacity in hens was correlated with the morphological differences of the oviduct and uterus-vagina junction.

The fertilization rate is an important index to evaluate the reproductive capacity of hens, which is mainly affected by semen quality, timing of artificial insemination (AI), and the ability to store sperm. A high sperm storage (SS) capacity can extend the interval, reduce the frequency, and decrease the labor costs of AI. However, relatively few studies have investigated the SS capacity of hens. Therefore, the aims of the present study were to identify factors influencing the SS capacity of Guangxi partridge chickens and to explore the impact of the sperm count in different sections of the oviduct and sperm storage tubules (SSTs), in addition to the number and surface area of SSTs on SS capacity at different time points after AI. We found that SS capacity was positively correlated to the egg production rate (P < 0.01) and body length (P < 0.05). On post-AI days 5, 10, and 15, the sperm count was higher in uterus-vagina junction (UVJ) than the magnum, isthmus, and infundibulum (P < 0.01), but gradually decreased over time. Also, the duration of SS and the sperm count of the UVJ was greater in the high SS group than the low SS group (P < 0.05). Histopathological analysis of the UVJ showed that the number and surface area of the SSTs (P < 0.01), as well as the proportion of SSTs containing sperm, were greater in the high SS group at all time points post AI (P < 0.01), while the proportion of SSTs containing sperm gradually decreased over time. Collectively, these results highlight the potential for selective breeding of SS capacity and show that SS capacity is related to laying performance and body length of Guangxi partridge hens. In addition, SS capacity was positively correlated to the surface area of SSTs and the proportion containing sperm. A greater sperm count stored in the UVJ was correlated to more sperm transported to the infundibulum and subsequent greater SS capacity of hens.

Phylogenetic and taxonomic updates of Agaricales, with an emphasis on Tricholomopsis.

The order Agaricales was divided into eight suborders. However, the phylogenetic relationships among some suborders are largely unresolved, and the phylogenetic positions and delimitations of some taxa, such as Sarcomyxaceae and Tricholomopsis, remain unsettled. In this study, sequence data of 38 genomes were generated through genome skimming on an Illumina sequencing system. To anchor the systematic position of Sarcomyxaceae and Tricholomopsis, a phylogenetic analysis based on 555 single-copy orthologous genes from the aforementioned genomes and 126 publicly accessible genomes was performed. The results fully supported the clustering of Tricholomopsis with Phyllotopsis and Pleurocybella within Phyllotopsidaceae, which formed a divergent monophyletic major lineage together with Pterulaceae, Radulomycetaceae, and Macrotyphula in Agaricales. The analysis also revealed that Sarcomyxaceae formed a unique major clade. Therefore, two new suborders, Phyllotopsidineae and Sarcomyxineae, are proposed for the two major lineages. Analyses of 450 single-copy orthologous genes and four loci suggested that Tricholomopsis consisted of at least four clades. Tricholomopsis is subsequently subdivided into four distinct sections. Seventeen Tricholomopsis species in China, including six new species, are reported. Conoloma is established to accommodate T. mucronata. The substrate preference of Tricholomopsis species and the transitions of the pileate ornamentations among the species within the genus are discussed.

Four New Species and One New Record of Thelephora from China.

Species of the genus Thelephora (Thelephorales, Thelephoraceae) are ectomycorrhizal symbionts of coniferous and broad-leaved plants, and some of them are well-known edible mushrooms, making it an exceptionally important group ecologically and economically. However, the diversity of the species from China has not been fully elucidated. In this study, we conducted a phylogenetic analysis based on the internal transcribed spacer (ITS) regions, using Maximum Likelihood and Bayesian analyses, along with morphological observations of this genus. Four new species from China are proposed, viz., T. dactyliophora, T. lacunosa, T. petaloides, and T. pinnatifida. In addition, T. sikkimensis originally described from India is reported for the first time from China. Thelephora dactyliophora, T. pinnatifida, and T. sikkimensis are distributed in subtropical forests and mainly associated with plants of the families Fagaceae and Pinaceae. Thelephora lacunosa and T. petaloides are distributed in tropical to subtropical forests. Thelephora lacunosa is mainly associated with plants of the families Fagaceae and Pinaceae, while T. petaloides is mainly associated with plants of the family Fagaceae. Line drawings of microstructures, color pictures of fresh basidiomes, and detailed descriptions of these five species are provided.

Molecular phylogeny and morphology reveal two new entomopathogenic species of Ophiocordyceps (Ophiocordycipitaceae, Hypocreales) parasitic on termites from China.

Two new termite-pathogenic species, Ophiocordycepsglobiperitheciata and O.longistipes, are described from Yunnan Province, China. Six-locus (ITS, nrSSU, nrLSU, tef-1α, rpb1 and rpb2) phylogenetic analyses in combination with morphological observations were employed to characterize these two species. Phylogenetically, O.globiperitheciata is most closely related to Hirsutellacryptosclerotium and O.communis, whereas O.longistipes shares a sister relationship with O.fusiformis. However, O.globiperitheciata differs from H.cryptosclerotium by parasitizing Blattodea and producing clavate, unbifurcated stromata. Ophiocordycepsglobiperitheciata is distinguished from O.communis by multiple stromata, shorter asci and ascospores. Ophiocordycepslongistipes differs from O.fusiformis in producing larger stromata, perithecia, asci and ascospores, as well as smaller citriform or oval conidia. Morphological descriptions of the two new species and a dichotomous key to the 19 termite-pathogenic Ophiocordyceps species are presented.

High Muscle Expression of IGF2BP1 Gene Promotes Proliferation and Differentiation of Chicken Primary Myoblasts: Results of Transcriptome Analysis.

Muscle development is a multifaceted process influenced by numerous genes and regulatory networks. Currently, the regulatory network of chicken muscle development remains incompletely elucidated, and its molecular genetic mechanisms require further investigation. The Longsheng-Feng chicken, one of the elite local breeds in Guangxi, serves as an excellent resource for the selection and breeding of high-quality broiler chickens. In this study, we conducted transcriptome sequencing of the pectoral muscles of Longsheng-Feng chickens and AA broiler chickens with different growth rates. Through comprehensive bioinformatics analysis, we identified differentially expressed genes that affect muscle growth and showed that IGF2BP1 is a key participant in chicken muscle development. Subsequently, we employed QRT-PCR, EdU staining, and flow cytometry to further investigate the role of IGF2BP1 in the proliferation and differentiation of chicken myogenic cells. We identified 1143 differentially expressed genes, among which IGF2BP1 is intimately related to the muscle development process and is highly expressed in muscle tissues. Overexpression of IGF2BP1 significantly promotes the proliferation and differentiation of chicken primary myoblasts, while knockdown of IGF2BP1 significantly inhibits these processes. In summary, these results provide valuable preliminary insights into the regulatory roles of IGF2BP1 in chicken growth and development.

Phylogeny, Taxonomy and Evolutionary Trade-Offs in Reproductive Traits of Gomphoid Fungi (Gomphaceae, Gomphales).

Although functional ecology is a well-established field, our understanding of the evolutionary and ecological significance of the reproductive traits in macrofungi is still limited. Here, we reconstructed a phylogeny tree of gomphoid fungi in the narrower sense, including the species of the genera Gomphus and Turbinellus and used it to uncover the evolution of reproductive traits. Our analyses indicated that fungal fruit bodies and spores did not enlarge at a steady rate over time. Early gomphoid fungi essentially maintained their fruit body size, spore size and spore shape through the Mesozoic. In the Cenozoic, gomphoid fungi acquired significantly larger and more spherical spores by simultaneously expanding in length and width, with the fruit body size first decreasing and then enlarging. We argue that these trade-offs were driven by the effect of biological extinction and the dramatic climate changes of the Cenozoic. Gomphoid fungi initially increased in spore size and fruit body number as extinction survivors filled vacant niches. Both fruit bodies and spores eventually became larger as ecosystems saturated and competition intensified. One new species of Gomphus and nine new species of Turbinellus are described.

Identification of genes involved in regulating the development of feathered feet in chicken embryo.

The genetic and developmental factors driving the diverse distribution and morphogenesis of feathers and scales on bird feet are yet unclear. Within a single species, Guangxi domestic chickens exhibit dramatic variety in feathered feet, making them an accessible model for research into the molecular basis of variations in skin appendages. In this study, we used H&E staining to observe the morphogenesis of feathered feet, scaled feet and wings skin at different embryonic stages in Longsheng-Feng chickens and Guangxi Partridge chickens. We selected 4 periods (E6, E7, E8, and E12) that play an important role in feather development and performed transcriptome sequencing to screen for candidate genes associated with feathered feet. Through comparison and analysis of transcriptome data, we identified a set of differently expressed genes (DGEs), which were enriched in appendage organ development, hindlimb morphogenesis, activation of transcription factor binding, and binding of sequence-specific DNA in the cis-regulatory region. In addition, we identified some feathered feet-related genes by analyzing the classical signaling pathways that regulate feather development. Finally, we identified candidate genes that regulate feathered feet formation, which include TBX5, PITX1, ZIC1, FGF20, WNT11, WNT7A, WNT16, and SHH. Interestingly, we found that TBX5 was significantly overexpressed in the skin of the feathered feet and had the highest expression at E7 (P < 0.01), whereas PITX1 expression was significantly reduced at E7(P < 0.01). It is hypothesized that TBX5 and PITX1 regulate the development of hair follicles through the Wnt/ß-catenin signaling pathway at E7. Our results provide a theoretical basis for investigating the molecular regulatory mechanisms underlying the formation of chicken feathered feet.