Skip to main content
Download PDF

Analysis of the differentially expressed genes in the combs and testes of Qingyuan partridge roosters at different developmental stages

BMC Genomics volume 25, Article number: 33 (2024) Cite this article

Abstract

Background

The sexual maturity of chickens is an important economic trait, and the breeding of precocious and delayed puberty roosters is an important selection strategy for broilers. The comb serves as an important secondary sexual characteristic of roosters and determines their sexual precocity. Moreover, comb development is closely associated with gonad development in roosters. However, the underlying molecular mechanism regulating the sexual maturity of roosters has not yet been fully explored.

Results

In order to identify the genes related to precocious puberty in Qingyuan partridge roosters, and based on the synchrony of testis and combs development, combined with histological observation and RNA-seq method, the developmental status and gene expression profile of combs and testis were obtained. The results showed that during the early growth and development period (77 days of age), the development of combs and testis was significant in the high comb (H) group versus the low comb (L) group (p < 0.05); however, the morphological characteristic of the comb and testicular tissues converged during the late growth and development period (112 days of age) in the H and L groups. Based on these results, RNA-sequencing analysis was performed on the comb and testis tissues of the 77 and 112 days old Qingyuan Partridge roosters with different comb height traits. GO and KEGG analysis enrichment analysis showed that the differentially expressed genes were primarily enriched in MAPK signaling, VEGF signaling, and retinol metabolism pathways. Moreover, weighted correlation network analysis and module co-expression network analysis identified WNT6, AMH, IHH, STT3A, PEX16, KPNA7, CATHL2, ROR2, PAMR1, WISP2, IL17REL, NDRG4, CYP26B1, and CRHBP as the key genes associated with the regulation of precocity and delayed puberty in Qingyuan Partridge roosters.

Conclusions

In summary, we identified the key regulatory genes of sexual precocity in roosters, which provide a theoretical basis for understanding the developmental differences between precocious and delayed puberty in roosters.

Peer Review reports

Background

In animals, sexual maturation is accompanied by aging, changes in tissue morphology and an increase in reproductive capacity [1,2,3]. Sexual maturity of the chickens is an important economic trait, and the breeding of precocious or delayed puberty individuals is an important selection strategy for broilers. Economically, the reproductive ability of roosters is more important than that of hens, in a breeding flock, as the roosters can fertilize many eggs [4]. In general, the first crow of the rooster is considered an indicator of its sexual maturity [5]. The primary reason for the occurrence of precocious puberty among roosters is the premature activation of the hypothalamic-pituitary-gonadal axis and the secretion of multiple sex hormones, which lead to the development of secondary sexual characteristics [6]. With multiple generations of breeding, the age of mating advanced and the comb size increased in the precocious puberty roosters, leading to many economic benefits [7, 8]. Qingyuan Partridge chicken is a renowned high-quality chicken breed that is well-known for its delicious meat in China and neighboring countries [9]. To meet the market’s demand, it is essential to breed with exceptional growth performance and early maturity of Qingyuan Partridge roosters. Therefore, it is necessary to analyze the underlying molecular mechanism leading to precocious puberty in roosters.

The comb is an important secondary sexual characteristic of the roosters and serves as an indicator of the reproductive capacity and sexual precocity of chickens [10], thus it can be used to determine the sexual maturity of roosters for broiler production.Testes are the primary reproductive organ of male animals, and studies have shown an association between testis weight and sexual maturity [4]. Moreover, studies have shown that comb development is closely associated with gonad development in roosters and that testes weight is correlated with comb development in roosters aged < 98 days [11]. A study found that roosters with larger combs typically have higher testicular weights, semen production, and androgen levels [12]. Previous studies have shown that the genes Hydroxyacid Oxidase 1 (HAO1), bone morphogenetic protein-2 (BMP2), SRY-box transcription factor 5 (SOX5), and androgen receptor (AR) affect comb morphology and development [13,14,15,16]; however, their roles in the regulation of sexual precocity of Qingyuan Partridge roosters remain to be studied.

In chickens, sexual maturation is regulated by a complex network of genes [17]; however, studies on precocious puberty in roosters are limited compared to those in hens. A study found that an allele in the comb trait locus that leads to an increase in comb weight delays the initiation of sexual maturity in chickens [18]. In addition, studies have found that GTPase-activating Rap/Ran-GAP domain-like 1(GARNL1) [19], Dehydrogenase/Reductase (SDR family)Member 12 (DHRS12) [20], and G-protein coupled receptor 39 (GPR39) [21] are positively correlated, while go-nadotropin-releasing hormone-I (GnRH-I) [22] is negatively correlated with the egg laying age and ovarian development in hens. Moreover, a study found that Transcription factor 21 (TCF21) [23] and gamma-aminobutyric acid A receptor alpha 1 (GABRA1) [21] can stimulate testicular growth and development in roosters. However, the underlying molecular mechanism associated with the onset of sexual precocity in roosters remains unclear, necessitating further studies on the network of genetic factors that lead to this trait.

In this study, we used RNA-sequencing (RNA-seq) analysis to investigate the differential gene expression in the combs and testis of roosters with different comb heights. Thereafter, we identified key regulatory genes associated with precocious puberty and delayed puberty in roosters, which may provide a theoretical basis to improve the breeding efficiency of high-quality chickens and reduce the cost of breeding.

Result

Morphological differences in the comb and testes of Qingyuan partridge roosters at different developmental stages

This study demonstrated that the H group of Qingyuan Partridge roosters exhibited significantly higher testicular weight and testicular index during the rapid growth period (56 days of age) and early growth and development period (77 days of age) compared to the L group. However, no statistically significant difference was observed between these two parameters in the later stages of development (112 days of age) (Table 1). HE staining of the comb tissues (Fig. 1A) showed that the capillaries in the comb expanded, the collagen fibers became sparse, and the number of fibroblasts decreased with an increase in age. In addition, significant changes were observed in the comb dermis during 56 and 77 days of age. Comparison the comb tissues at the same developmental stage revealed that the capillaries in the comb tissues expanded earlier, and the collagen fibers decreased much earlier in the H group compared to the L group. However, the morphological characteristic changes in the comb tissues gradually converged between the H group and the L group in 112 days of age.

Table 1 Significance analysis of the correlated characteristics between Qingyuan partridge roosters at different age stages
Full size table
Fig. 1
figure 1

Morphological differences in the comb and testis tissues of the Qingyuan partridge roosters at different developmental stages. (A and B) Morphological analysis of the comb (A) and testis (B) tissues of the high-and low-comb groups at different developmental stages using hematoxylin and eosin staining. Original magnification: a-f in Figure A and B: 20x, scale bar = 100 μm and g-l in Figure A and B: 40x, scale bar = 50 μm. In Figure B, ▲ indicates spermatogonium and → indicates sperm cells

Full size image

The diameter of the seminiferous tubules and the number of spermatogenic epithelium cell layers gradually increased with age in testis tissues by HE staining (Fig. 1B).

The testicular development was faster in the H group compared with the L group at 56 and 77 days of age, with the presence of sperm cells observed in the H group at 77 days. Additionally, we have observed that the H group exhibited an earlier onset and greater diameter of the complete testicular seminiferous tubule structure in comparison to the L group. However, the morphological changes in the testis tissues of the H and L groups gradually converged at 112 days.

Altogether, these results suggest that the rapid development of the comb and the testis tissues is synchronous during the early growth and development periods in Qingyuan Partridge roosters. Furthermore, the results reveal that the morphology of the comb and testicular tissues is synchronized between the H and L groups during the late growth and development period. RNA-Seq analysis of the comb and testis tissues in the H and L groups was performed at 77 and 112 days of age, to explore the relationship between the spatiotemporal regulation of the key genes associated with comb and testis development and sexual maturity of roosters.

DEGs in the comb and testis tissues of the H and L groups at the same developmental stage

A principal component analysis (PCA) was conducted on the fragments per kilobase per million mapped fragments of each comb (Fig. 2A) and testis (Fig. 2C) sample, and it was found that all the subgroups had excellent within-group repeatability and between-group variability. The Wayne diagrams in Fig. 2B and D show the number of DEGs in the comb and testis tissues, respectively. RT-qPCR analysis of the DEGs (n = 6 each) from the comb and testis tissues of the 77 and 112-day-old roosters (Figure S1A-S1D) revealed a consistency between the RNA-seq and gene expression results, thus proving the reliability of the results.

Fig. 2
figure 2

Differentially expressed genes (DEGs) in the comb and testis tissues of the high-and low-comb groups at the same developmental stage. (A, C) PCA analysis of RNA-Seq data in the comb (A) and testis (C) tissues; (B, D) The Wayne map of differentially expressed genes in the comb (B) and testis (D) tissues; (E and F) The volcano maps of the DEGs in the comb tissues of the 77 (E) and 112 (F) day-old high- and low-comb groups and (G and H) the volcano maps of DEGs in the testis tissues of the 77 (G) and 112 (H) day-old high- and low-comb groups

Full size image

In this study, we identified a total of 222 and 44 DEGs (FDR < 0.05 and|log2FoldChange| >1) in the comb (126 upregulated and 96 downregulated genes; Fig. 2E and S2A) and testis (16 upregulated and 28 downregulated genes; Fig. 2G and S2C) tissues of the 77-days-old H and L group roosters, respectively. In contrast, we identified a total of 185 and 65 DEGs in the comb (93 upregulated and 92 downregulated genes; Fig. 2F and S2B) and testis (39 upregulated and 26 downregulated genes; Fig. 2H and S2D) tissues of the 112-days-old H and L group roosters, respectively.

These DEGs were subjected to GO and KEGG enrichment analyses to determine their biological functions. GO enrichment analysis showed that the DEGs were significantly enriched in signal-organism processes, cellular processes, response to stimulus, reproduction and developmental process (Fig. 3A and B). KEGG enrichment analysis showed that the DEGs in the comb tissues of the H and L groups at 77 days of age were significantly enriched in the adipocytokine signaling pathway, cytokine-cytokine receptor interaction, and Wnt signaling pathway (Fig. 3C) and that the DEGs in the testis tissues of the H and L groups at 77 days of age were significantly enriched in cytokine-cytokine receptor interaction and the MAPK signaling pathway (Fig. 3E). Furthermore, KEGG enrichment analysis showed that the DEGs in the comb tissues of the H and L groups at 112 days of age were significantly enriched in ECM-receptor interaction, glycosaminoglycan binding proteins, nuclear receptors, and PPAR signaling pathway (Fig. 3D) and that the DEGs in the testis tissues of the H and L groups at 112 days of age were significantly enriched in Drug metabolism-cytochrome P450 and retinol metabolism pathway (Fig. 3F).

Fig. 3
figure 3

Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of DEGs in the comb and testes tissues of the high- and low-comb groups at the same developmental stage. (A and B) GO enrichment analysis of the DEGs in the comb (A) and testes (B) tissues in the molecular function, cellular composition, and biological processes categories; (C and D) KEGG enrichment analysis of the DEGs in the comb tissues of the 77 (C) and 112 (D) day-old high-and low-comb groups; (E and F) KEGG enrichment analysis of the DEGs in the testes tissues of the 77 (E) and 112 (F) day-old high-and low-comb groups

Full size image

Comparative analysis of the DEGs in the comb and testis tissues of the H and L groups at 77 and 112 days of age revealed that secreted frizzled related protein 1 (SFRP1), cathelicidin 2 (CATHL2), cytochrome P450 family 26 subfamily B member 1 (CYP26B1), and WNT6 were significantly differentially expressed in the comb tissues, while Indian hedgehog signaling molecule (IHH), WNT6, and receptor tyrosine kinase like orphan receptor 2 (ROR2) were significantly differentially expressed in the testis tissues of the H and L groups at 77 days of age, suggesting that these genes may be responsible for the differences in the precocity and delayed puberty in Qingyuan Partridge roosters. Additionally, we found that corticotropin releasing hormone binding protein (CRHBP), AMH, and HtrA serine peptidase 1 (HTRA1) were significantly differentially expressed in the comb tissues, while CYP26B1 and AMH were significantly differentially expressed in testis tissues of the H and L groups at 112 days of age, suggesting that these downregulated genes may be responsible for the compensatory growth during the late developmental stages of the sex late-maturing Qingyuan Partridge roosters.

DEGs in the comb and testis tissues of the H or L group at different developmental stages

We conducted DEGs analysis of the comb and testis tissues of the H or L group roosters at 77 and 112 days old. We identified a total of 869 and 1060 DEGs in the comb (332 upregulated and 537 downregulated genes; 4 A and 4B) and testis (492 upregulated and 568 downregulated genes; Fig. 4E and F) tissues of the H group at 77 and 112 days of age, respectively. In contrast, we identified 199 and 745 DEGs in the comb (53 upregulated and 146 downregulated genes; 4 C and 4D) and testis (167 upregulated and 578 downregulated genes; Fig. 4G and H) tissues of the L group at 77 and 112 days of age, respectively.

Fig. 4
figure 4

DEGs in the comb and testis tissues of the high-or low-comb height groups at different developmental stages. (A and B) The volcano map (A) and heat map (B) of the DEGs in the comb tissues of the H group at 77 and 112 days of age; (C and D) the volcano map (C) and heat map (D) of the DEGs in the comb tissues of the L group at 77 and 112 days of age; (E and F) the volcano map (E) and heat map (F) of DEGs in the testes tissues of the H group at 77 and 112 days of age; and (G and H) the volcano maps (G) and heat map (H) of the DEGs in the testes tissues of the L group at 77 and 112 days of age

Full size image

These DEGs were subjected to GO and KEGG enrichment analyses to determine their biological functions. GO enrichment analysis showed that these DEGs were also significantly enriched in signal-organism processes, cellular processes, response to stimulus, and developmental processes (Fig. 5A and B). KEGG enrichment analysis showed that the DEGs in the comb tissues of the H group at 77 and 112 days of age were significantly enriched in glycosaminoglycan binding proteins, calcium signaling pathway, and VEGF signaling pathway (Fig. 5C) and that the DEGs in the testis tissues of the H group at 77 and 112 days of age were significantly enriched in ECM-receptor interaction, focal adhesion, and biosynthesis of unsaturated fatty acid (Fig. 5E). Additionally, KEGG enrichment analysis showed that the DEGs in the comb tissues of the L group at 77 and 112 days of age were significantly enriched in steroid biosynthesis, cytokine-cytokine receptor interaction, MAPK signaling pathway and cell adhesion molecules pathway (Fig. 5D) and that the DEGs in the testis tissues of the L group at 77 and 112 days of age were significantly enriched in focal adhesion and the VEGF signaling pathway (Fig. 5F).

Fig. 5
figure 5

GO and KEGG analysis of DEGs in the comb and testis tissues of the high- or low-comb groups at different developmental stages. (A and B) GO enrichment analysis of the DEGs in the comb (A) and testes (B) tissues in the molecular function, cellular composition, and biological processes categories; (C and D) KEGG enrichment analysis of the DEGs in the comb tissues of the high- (C) and low-comb groups (D) at 77 and 112 days of age; (E and F) KEGG enrichment analysis of the DEGs in the testes tissues of the high- (E) and low-comb groups (F) at 77 and 112 days of age

Full size image

We conducted a comparative analysis of the DEGs in the comb and testis tissues of the H or L group at 77 and 112 days of age and identified 33 DEGs in the comb and testis tissues of the H group. Among these, WNT1 inducible signaling pathway protein 2 (WISP2), peptidase domain containing associated with muscle regeneration 1 (PAMR1), interleukin 17 receptor E like (IL17REL), interleukin-13 receptor subunit alpha-2 precursor (IL13RA2), gliomedin (GLDN), myelin proteolipid protein (PLP1), galactose mutarotase (GALM), fibromodulin precursor (FMOD), sarcoglycan alpha (SGCA), actin gamma-enteric smooth muscle (ACTG2), and desmin (DES) were significantly downregulated and may be associated with the regulation of the precocious phenotype of the comb and testis tissues in the Qingyuan Partridge roosters. Additionally, we found 15 DEGs in the comb and testis tissues of the L groups, among which NDRG family member 4 (NDRG4), WNT6, and AR may be responsible for the regulation of the compensatory growth during the late developmental stages of the sex late-maturing in Qingyuan Partridge roosters.

WGCNA of the DEGs in the comb and testis tissues

To further explore the genes related to the difference of sexual maturity in different growth stages of Qingyuan Partridge roosters, we conducted a WGCNA of the DEGs in the comb and testis tissues of the H and L groups at different developmental stages. When the scale-free topology model fit (R2) reached 0.8, the soft thresholding power (β) was 4 (Fig. 6A). Hierarchical clustering among all the tissues was presented as an array tree (Fig. 6B). A total of 34,700 DEGs were clustered into 21 modules (Fig. 6C), among which, the brown module contained the highest number of genes (12,123 genes), while the orangered4 module contained the lowest number of genes (55 genes; Fig. 6D). A sample expression model heat map was used to compare the modules with similar expression coefficients for different traits at different ages, and the blue, darkolivegreen, sienna3, and plum1 modules were used to represent the 77-day-old H group, 77-day-old L group, 112-day-old H group, and 112-day-old L group, respectively (Fig. 6E).

Fig. 6
figure 6

Weighted correlation network analysis of the DEGs in the comb and testes tissues. (A) Analysis of the scale-free fitting index of various soft-thresholding powers (left) and the average connectivity of various soft-thresholding powers (right); (B) sample tree obtained from hierarchical clustering; (C) gene dynamic shearing clustering tree, each color represents a module; (D) classification of the DEGs into 21 modules according to the gene co-expression patterns; (E) the heat map of gene expression patterns of different phenotypic tissues at different developmental stages

Full size image

GO enrichment analysis showed that the DEGs in the blue module were significantly enriched in cellular processes, cell part, and signal-organism processes (Fig. 7A); DEGs in the darkolivegreen module were significantly enriched in cell and cell part (Fig. 7B); DEGs in the sienna3 module were significantly enriched in organelle and cellular component organization or biogenesis (Fig. 7C); and DEGs in the plum1 module were significantly enriched in organelle and structural molecule activity (Fig. 7D). Furthermore, we found that 36 GO pathways were enriched at both the developmental stages in the combs of the H group, while 14 GO pathways were enriched only at 77 days of age. Moreover, we found that 22 GO pathways were enriched at both the developmental stages in the combs of the L group, while 7 GO pathways were enriched only at 77 days of age and 2 GO pathways were enriched only at 112 days of age (Table S2). Further analysis of the DEGs revealed that SFRP1, WNT6, KPNA7, and CCN3, associated with gonadal, epidermal, and vascular development, may play a role in regulating precocious and delayed puberty in Qingyuan Partridge roosters.

Fig. 7
figure 7

GO and KEGG analysis of the high-expression modules. (A-D) GO enrichment analysis of the blue (A), darkolivegreen (B), sienna3 (C), and plum1 (D) modules in the molecular function, cellular composition, and biological processes categories; (E-H) KEGG enrichment analysis of the blue (E), darkolivegreen (F), sienna3 (G), and plum1 (H) modules

Full size image

KEGG enrichment analysis showed that the DEGs in the blue module were significantly enriched in ribosome, steroid hormone biosynthesis, and focal adhesion pathway (Fig. 7E); DEGs in the darkolivegreen module were significantly enriched in arachidonic acid metabolism and cell adhesion molecules (Fig. 7F); DEGs in the Sienna3 module were significantly enriched in cAMP signaling pathway, longevity regulating pathway-multiple species, and long-term potentiation (Fig. 7G); and genes in the plum1 module were significantly enriched in steroid biosynthesis (Fig. 7H). We further found that 13 KEGG pathways were enriched only at 77 days of age in the H group, while 3 KEGG pathways were enriched only at 77 days of age and 1 KEGG pathway was enriched only at 112 days of age in the L group (Table S3). These results suggest that the pathways enriched only in the H group at 77 days of age may be responsible for the differences in the comb height phenotypes between precocious and delayed puberty in Qingyuan Partridge roosters, while the pathways enriched only in the L group at 112 days of age may be responsible for regulating the compensatory development in delayed puberty Qingyuan Partridge roosters.

Cytoscape was used to construct gene interaction networks of the four modules, and Cytohubba analysis found that STT3 oligosaccharyltransferase complex catalytic subunit A (STT3A) and peroxisomal biogenesis factor 1(PEX1) were the hub genes in the blue module, while karyopherin subunit alpha 7 (KPNA7), ATHL, and ENSGALG00000012072 were the hub genes in the darkolivegreen, sienna3, and plum1 modules, respectively (Fig. 8A and D). In addition, previously discovered WNT6, AMH, IHH, ROR2, SFRP1, CYP26B1, PAMR1, IL17REL, WISP2, NDRG4, and HTRA1 were identified as the hub genes in the blue module, while CRHBP was identified as the hub gene in the sienna3 module. Altogether, these results reveal that these hub genes may be responsible for regulating the sexual maturity of Qingyuan Partridge roosters, via pathways associated with blood vessels, skin connective tissue development, energy storage, and cell growth.

Fig. 8
figure 8

Construction of the co-expression subnetwork and screening of the hub genes. (A-D) Gene co-expression network of the blue (A), darkolivegreen (B), sienna3 (C), and plum1 (D) modules

Full size image

Discussion

Comb height is a heritable trait [16], and it is usually used as the main reference trait in the selection of sexual precocity in roosters. Moreover, comb height can be used as an important indicator of sexual maturity and testis development [10]. In this study, we found that capillary dilation and loosening of collagen fiber arrangement occurred earlier in the H group compared to the L group at 56 and 77 days of age, consistent with the results observed by Yoshioka [24]. The density difference of collagen fibers may be the primary reason for the high and low phenotypes of the combs. However, the morphological changes in the comb tissues of the H and L groups converged gradually at 112 days of age, although the underlying mechanism is unknown and needs to be further investigated. Additionally, we found that the sexual maturation of the L group was delayed compared to that of the H group during early developmental stages; however, the L group showed compensatory growth during the late developmental stage. In our previous study, we found that the body weight of the low-comb roosters was higher than that of the high-comb roosters, during the later stages of growth and development. Moreover, studies have shown that slow-growing chickens contain higher levels of inosinic acid, suggesting that their flavor may be better than that of fast-growing chickens [25]. Therefore, further studies can be conducted on the breeding of low-comb roosters into a heavier and more flavorful late-maturing strain.

In poultry industry, the reproductive performance of roosters is crucially important [26]. In recent years, a multitude of scholars have conducted extensive research on the genetic mechanism underlying the process of sexual maturation in chickens. Previous studies have reported WNT6 and AMH as key genes associated with gonadal development in chickens [16, 27,28,29]. Similarly, the present study reveals that WNT6 exhibits the potential to impede the development of comb and testis in the H group during the initial stages of development, whereas AMH demonstrates the capacity to foster the advancement of comb and testis development in the L group during the later stages of development. The DEGs analysis of the RNA-seq data further demonstrated that the expression levels of genes WISP2, PAMR1, and IL17REL were significantly reduced in the H group, suggesting their potential role in regulating early maturation traits in roosters. Conversely, genes CRHBP, CYP26B1and NDRG4 exhibited differential expression exclusively in the L group, possibly indicating their involvement in regulating compensatory testicular and comb growth during the late developmental stages in the L group. WISP2 exhibits significant expression levels within the reproductive tissues and has been demonstrated to play a role in the modulation of avian development, cellular adhesion and proliferation, as well as the remodeling of the extracellular matrix [30]. Although the role of PAMR1 and IL17REL in chicken growth and development remains to be investigated, it is speculated that the downregulation of WISP2, PAMR1, and IL17REL genes observed in H group may indicate their inhibitory effects on comb and testicular development.

CYP26B1 is a retinoic acid-degrading enzyme that can regulate the initiation of meiosis in chicken germ cells [31]. NDRG4 is a downstream regulatory gene of N-myc and is involved in the regulation of cell proliferation, differentiation, and tissue development in chickens [32]. Therefore, it is postulated that the differentially expressed CYP26B1 and NDRG4 genes, exclusively observed in the L group, potentially participate in governing the compensatory growth of the L group’s testes and combs during the later phases of development by modulating cellular proliferation and tissue development.

GO and KEGG analysis revealed that the DEGs were primarily enriched in MAPK signaling, VEGF signaling, and retinol metabolism pathways. MAPK signaling is a common reproductive-related signaling pathway and is associated with angiogenesis and chondrocyte development [33]. Moreover, some studies have shown that long non-coding RNAs regulate the formation of primordial germ cells by activating the MAPK signaling pathway in chickens [34]. The VEGF signaling pathway is crucial for angiogenesis, development, and maintenance of physiological homeostasis [35]. Lastly, the retinol metabolic pathway plays an important role in the differentiation of male germ cells in chickens [36]. The presence of enriched DEGs within these pathways may account for the observed early onset of capillary anastomosis, loosening of the collagen fiber network, and intact structure of testicular seminiferous tubules in H-group roosters as compared to L-group roosters. Therefore, these DEGs can be considered potential candidate genes for the regulation of precocious and delayed puberty in Qingyuan Partridge roosters.

WGCNA has been shown to play an important role in screening key genes and signaling pathways [37, 38]. To our knowledge, this is the first study to use WGCNA to screen candidate genes that affect the development of comb and testis tissues of Qingyuan Partridge roosters. WGCNA of the DEGs identified four hub genes, including STT3A, PEX16, KPNA7, and CATHL2, associated with comb and testis development in the roosters. STT3A is involved in the regulation of glycosylation, which may affect the development and differentiation of cells and tissues [39, 40]. PEX16 plays a key role in the development of white adipocytes and lipid homeostasis, which may regulate energy storage during the development of sexual maturity in roosters [41]. KPNA7 is associated with cell mitosis, and its deletion can significantly reduce cell growth [42]. It is postulated that these core genes may exert an impact on sexual maturation via their regulatory roles in protein post-translational modification, adipogenesis, or cellular proliferation. CATHL2 has been reported to be associated with innate immunity in chickens, with the downregulation of CATHL2 leading to a decrease in innate immune function in chickens [43]. However, additional research is required to investigate its involvement in the regulation of sexual maturity. Moreover, we found additional hub genes (including IHH and ROR2) that were associated with comb and testis development in roosters. IHH plays a crucial role in regulating the proliferation and hypertrophy of chicken chondrocytes [44]. ROR2 is a signaling receptor for Wnt ligands and plays an important role in limb development [45]; additionally, it has been reported as a potential candidate gene for testicular growth and development in chickens [46]. The candidate genes associated with chicken chondrogenesis have the potential to impact the formation of dense connective tissue, consequently resulting in the accelerated onset of early sexual maturity in Qingyuan Partridge roosters.

Conclusion

In summary, these results provide a comprehensive transcriptomic analysis of the combs and testes of Qingyuan Partridge roosters at the age of 77 and 112 days. WGCNA was used to identify the candidate genes (WNT6, AMH, IHH, STT3A, PEX16, KPNA7, CATHL2, ROR2, PAMR1, WISP2, IL17REL, NDRG4, CYP26B1, and CRHBP) that affect precocious and delayed puberty in Qingyuan Partridge roosters. Therefore, our findings provide the basis for understanding the underlying molecular mechanism associated with sexual maturity and its related traits in chickens.

Materials and methods

Animals and tissues

Qingyuan Partridge roosters were obtained from the Original Breeding Farm of Guangdong Tiannong Food Co., Ltd. (Guangdong, China). The management procedures were conducted according to the standard requirements for the chicken breeding farm. The lighting conditions have equal periods of bright and dark. The roosters (42 days old) were divided into high-comb (H group; 22–29 mm, n = 260) and low-comb (L group; 10–17 mm, n = 210) groups, according to the comb height. The comb height difference between the two groups is 9.97 mm, which indicates a significant distinction.15 roosters were slaughtered from each group on days 56 (the rapid growth period of Qingyuan Partridge roosters), 77 (the early growth and development period), and 112 (the later stages of development), and comb and testicular tissues were collected. The tissues were stored in a 4% paraformaldehyde solution for histological analysis or frozen in liquid nitrogen and stored at − 80 ℃ for RNA extraction.

Determination of body weight and coefficients of testis

The pertinent characteristics of roosters at the ages of 56, 77, and 112 days were assessed utilizing an electronic scale, and subsequently computed in the following manner. Body weight was measured after roosters fasted for ~ 12 h. The comb height is determined by employing a vernier caliper to precisely measure the vertical extent between the uppermost comb tip and the base of the comb. Testicular weight was measured by extracting both left and right testicles (excluding the epididymis) in their entirety after dissecting the carcass. The long and short diameters of each testicle were measured using a vernier caliper. The long diameter was determined as the linear distance between the upper and lower extremities of the testicle, whereas the short diameter was determined as the linear distance between the anterior and posterior edges of the testicle. Subsequently, the testicular area was calculated as (short diameter×long diameter). The calculation for the testis index is as follows: Testis index (%) = Testis weight(mg)/body weight(g) × 100. All measured values are expressed as mean ± SD.

Histological analysis

Tissues for histological analysis were prepared as previously described [47]. The comb and testis tissues were fixed in 4% paraformaldehyde solution for 10 d. Thereafter, the tissues were paraffin-embedded, sectioned (5–7 μm), hematoxylin-eosin (HE) stained, and observed under a 400X optical microscope (Leica, Heidelberg, Germany).

RNA extraction and sequencing

The total RNA was extracted from the comb and testes tissues using the Total RNA Extraction Kit (Sangong, Shanghai, China) according to the manufacturer’s instructions. The concentration and quality of all the RNA tissues were determined using a micro nucleic acid/protein concentration tester (Eppendorf, Hamburg, Germany). Subsequently, equal amounts of total RNA of each group were pooled for comb and testes tissues and used to construct the libraries for transcriptome analysis. RNA-seq of these tissues was conducted by Genedenovo Biotech (Guangzhou, China). FASTP was used for quality control processing to filter the low-quality data, and HISAT2 was used to align clean data to the chicken reference genome. DESeq2 was used to analyze differential expression among the tissues at FDR < 0.05 and|log2FC| >1. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed to determine the biological functions of the differentially expressed genes (DEGs).

Construction of co-expression network

Two tissues from the same developmental stage were pooled for the weighted correlation network analysis (WGCNA). The co-expression network was constructed for the mRNA expression profile data using the R package, WGCNA. The sample cluster tree was established to filter the outliers. The adjacency matrix was converted to a topological overlap matrix to reduce the effects of spurious associations and noise. Thereafter, the dynamic tree cut algorithm was used to identify the gene co-expression modules, with each module containing at least 50 genes. Finally, highly similar modules were merged and represented with different colors (merge cut height = 0.25), and module membership and gene significance were calculated to identify hub genes.

Real-time quantitative PCR validation of candidate genes

Real-time quantitative PCR (RT-qPCR) was used to detect the mRNA expression of the target genes. Total RNA from comb and testis tissues was reverse transcribed using the Evo M-MLV Reverse Transcriptase Kit (Ruizhen, Guangzhou, China) according to the manufacturer’s instructions. The qPCR reactions were performed using the AceQ qPCR SYBR Green Master Mix (Vazyme, Nanjing, China) on the CFX96 Touch RT-qPCR System (Bio-Rad, California, USA). All the primers (Table S1) used in this study were designed using the Primer Premier 5.0 software (Premier Biosoft, Palo Alto, CA, USA). The GAPDH mRNA was used as an internal control, and the relative mRNA expression was determined using the 2-∆∆Ct method.

Statistical analysis

Statistical data analysis in this study was performed using the SPSS 20.0 software package (IBM). The t-test and one-way analysis of variance (ANOVA) were used appropriately to test the result of the comparison. Statistical significance was declared at p < 0.05.

Data availability

The raw data of mRNA sequencing has been submitted to the NCBI Sequence Read Archive (SRA) with the accession number PRJNA1013206.

References

  1. Lu J, Li YF, Qu L, Ma M, Yang XD, Shen MM, Wang XG, Guo J, Hu YP, Dou TC, et al. Effects of energy-restricted feeding during rearing on sexual maturation and reproductive performance of Rugao layer breeders. Poult Sci. 2021;100(8):101225.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Wang Y, Yuan J, Sun Y, Li Y, Wang P, Shi L, Ni A, Zong Y, Zhao J, Bian S, et al. Genetic basis of sexual maturation heterosis: insights from ovary lncRNA and mRNA repertoire in Chicken. Front Endocrinol (Lausanne). 2022;13:951534.

    Article  PubMed  Google Scholar 

  3. Guerra-Carvalho B, Carrageta DF, Crisostomo L, Carvalho RA, Alves MG, Oliveira PF. Molecular mechanisms regulating spermatogenesis in vertebrates: Environmental, metabolic, and epigenetic factor effects. Anim Reprod Sci. 2022;246:106896.

    Article  CAS  PubMed  Google Scholar 

  4. Mohammadi V, Sharifi SD, Sharafi M, Mohammadi-Sangcheshmeh A. Effects of dietary L-carnitine on puberty indices in the young breeder rooster. Heliyon. 2021;7(4):e6753.

    Article  Google Scholar 

  5. Joseph NS, Robinson FE, Renema RA, Thorsteinson KA. Comb growth during sexual maturation in female broiler breeders. J Appl Poultry Res. 2003;12(1):7–13.

    Article  Google Scholar 

  6. Zhu S, Long L, Hu Y, Tuo Y, Li Y, Yu Z. GnRHa/Stanozolol Combined Therapy maintains normal bone growth in central precocious puberty. Front Endocrinol (Lausanne). 2021;12:678797.

    Article  PubMed  Google Scholar 

  7. Wei FY, Liao YY, Liu MJ. Significance of canopy area selection for improving performance of local chickens. Guangxi J Anim Husb Veterinary Med. 2003;2:78–9.

    Google Scholar 

  8. Cong L, Hua LI, Zheng-Fen Z, Guo-Hong H, Jie-Ling L, Hui YU. Comparative analysis of laying performance between early and late feathering strains in Qingyuan Patridge Chickens. J South Agric. 2019;50(7):1625–31.

    Google Scholar 

  9. Li Y, Mei H, Liu Y, Li Z, Qamar H, Yu M, Ma X. Dietary supplementation with Rutin alters meat quality, fatty acid Profile, antioxidant capacity, and Expression Levels of Genes Associated with lipid metabolism in breast muscle of Qingyuan Partridge Chickens. Foods. 2023;12(12):2302.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Dong X, Li J, Zhang Y, Han D, Hua G, Wang J, Deng X, Wu C. Genomic analysis reveals pleiotropic alleles at EDN3 and BMP7 involved in chicken comb Color and Egg Production. Front Genet. 2019;10:612.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Mukhtar N, Khan SH. Comb: an important reliable visual ornamental trait for selection in chickens. World Poult Sci J. 2012;68:425–34.

    Article  Google Scholar 

  12. Mcgary S, Estevez I, Russek-Cohen E. Reproductive and aggressive behavior in male broiler breeders with varying fertility levels. Appl Anim Behav Sci. 2003;82(1):29–44.

    Article  Google Scholar 

  13. Johnsson M, Gustafson I, Rubin CJ, Sahlqvist AS, Jonsson KB, Kerje S, Ekwall O, Kampe O, Andersson L, Jensen P, et al. A sexual ornament in chickens is affected by pleiotropic alleles at HAO1 and BMP2, selected during domestication. Plos Genet. 2012;8(8):e1002914.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Tillgren V, Ho JC, Onnerfjord P, Kalamajski S. The novel small leucine-rich protein chondroadherin-like (CHADL) is expressed in cartilage and modulates chondrocyte differentiation. J Biol Chem. 2015;290(2):918–25.

    Article  CAS  PubMed  Google Scholar 

  15. Boije H, Harun-Or-Rashid M, Lee YJ, Imsland F, Bruneau N, Vieaud A, Gourichon D, Tixier-Boichard M, Bed’Hom B, Andersson L, et al. Sonic hedgehog-signalling patterns the developing chicken comb as revealed by exploration of the pea-comb mutation. PLoS ONE. 2012;7(12):e50890.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Shen M, Qu L, Ma M, Dou T, Lu J, Guo J, Hu Y, Yi G, Yuan J, Sun C, et al. Genome-Wide Association Studies for comb traits in chickens. PLoS ONE. 2016;11(7):e159081.

    Article  Google Scholar 

  17. Tan YG, Xu XL, Cao HY, Zhou W, Yin ZZ. Effect of age at first egg on reproduction performance and characterization of the hypothalamo-pituitary-gonadal axis in chickens. Poult Sci. 2021;100(9):101325.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Wright D, Rubin C, Schutz K, Kerje S, Kindmark A, Brandstrom H, Andersson L, Pizzari T, Jensen P. Onset of sexual maturity in female chickens is genetically linked to loci associated with fecundity and a sexual ornament. Reprod Domest Anim. 2012;47(Suppl 1):31–6.

    Article  PubMed  Google Scholar 

  19. Shen X, Zeng H, Xie L, He J, Li J, Xie X, Luo C, Xu H, Zhou M, Nie Q, et al. The GTPase activating Rap/RanGAP domain-like 1 gene is associated with chicken reproductive traits. PLoS ONE. 2012;7(4):e33851.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Sun C, Lu J, Yi G, Yuan J, Duan Z, Qu L, Xu G, Wang K, Yang N. Promising loci and genes for yolk and ovary weight in chickens revealed by a genome-wide Association study. PLoS ONE. 2015;10(9):e137145.

    Article  Google Scholar 

  21. Sun Y, Liu R, Zhao G, Zheng M, Li P, Liu L, Wen J. Genome-wide linkage analysis identifies loci for testicle and ovary traits in chickens. Anim Biotechnol. 2018;29(4):309–15.

    Article  CAS  PubMed  Google Scholar 

  22. Xu H, Zeng H, Luo C, Zhang D, Wang Q, Sun L, Yang L, Zhou M, Nie Q, Zhang X. Genetic effects of polymorphisms in candidate genes and the QTL region on chicken age at first egg. BMC Genet. 2011;12:33.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Zhang H, Na W, Zhang HL, Wang N, Du ZQ, Wang SZ, Wang ZP, Zhang Z, Li H. TCF21 is related to testis growth and development in broiler chickens. Genet Sel Evol. 2017;49(1):25.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Yoshioka K, Watahiki Y, Kanie A, Tsujio M, Ikadai H, Kashimoto T, Mutoh K. Morphology of the cockerel’s comb after androgen administration. Br Poult Sci. 2010;51(2):185–94.

    Article  CAS  PubMed  Google Scholar 

  25. Huang Z, Zhang J, Gu Y, Cai Z, Feng X, Yang C, Xin G. Research progress on inosine monophosphate deposition mechanism in chicken muscle. Crit Rev Food Sci Nutr. 2022;62(4):1062–78.

    Article  CAS  PubMed  Google Scholar 

  26. Lotfi S, Fakhraei J, Mansoori YH. Dietary supplementation of pumpkin seed oil and sunflower oil along with vitamin E improves sperm characteristics and reproductive hormones in roosters. Poult Sci. 2021;100(9):101289.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Takada S, Mano H, Koopman P. Regulation of amh during sex determination in chickens: Sox gene expression in male and female gonads. Cell Mol Life Sci. 2005;62(18):2140–6.

    Article  CAS  PubMed  Google Scholar 

  28. Schubert FR, Mootoosamy RC, Walters EH, Graham A, Tumiotto L, Munsterberg AE, Lumsden A, Dietrich S. Wnt6 marks sites of epithelial transformations in the chick embryo. Mech Dev. 2002;114(1–2):143–8.

    Article  CAS  PubMed  Google Scholar 

  29. Lambeth LS, Morris K, Ayers KL, Wise TG, O’Neil T, Wilson S, Cao Y, Sinclair AH, Cutting AD, Doran TJ, et al. Overexpression of anti-mullerian hormone disrupts gonadal sex differentiation, blocks sex hormone synthesis, and supports Cell Autonomous Sex Development in the chicken. Endocrinology. 2016;157(3):1258–75.

    Article  CAS  PubMed  Google Scholar 

  30. Alami T, Liu JL. Metabolic effects of CCN5/WISP2 Gene Deficiency and transgenic overexpression in mice. Int J Mol Sci. 2021;22(24):13418.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Yousefi TS, Mahdavi SN, Lari R, Bahrami AR, Dehghani H. Meiotic initiation in chicken germ cells is regulated by Cyp26b1 and mesonephros. J Exp Zool B Mol Dev Evol. 2019;332(7):269–78.

    Article  Google Scholar 

  32. Zhu M, Zheng R, Guo Y, Zhang Y, Zuo B. NDRG4 promotes myogenesis via Akt/CREB activation. Oncotarget. 2017;8(60):101720–34.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Jahejo AR, Niu S, Zhang D, Ning GB, Khan A, Mangi RA, Qadir MF, Khan A, Li JH, Tian WX. Transcriptome analysis of MAPK signaling pathway and associated genes to angiogenesis in chicken erythrocytes on response to thiram-induced tibial lesions. Res Vet Sci. 2019;127:65–75.

    Article  CAS  PubMed  Google Scholar 

  34. Zuo Q, Jin J, Jin K, Zhou J, Sun C, Song J, Chen G, Zhang Y, Li B. P53 and H3K4me2 activate N6-methylated LncPGCAT-1 to regulate primordial germ cell formation via MAPK signaling. J Cell Physiol. 2020;235(12):9895–909.

    Article  CAS  PubMed  Google Scholar 

  35. Apte RS, Chen DS, Ferrara N. VEGF in Signaling and Disease: Beyond Discovery and Development. Cell. 2019;176(6):1248–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Qisheng Z, Dong, Li L, Zhang. Ahmed, Kamel, Elsayed, Chao. Study on the regulatory mechanism of the lipid metabolism pathways during chicken male germ cell differentiation based on RNA-seq. PLoS ONE. 2015;10(2):e0109469.

    Article  Google Scholar 

  37. Zhang FL, Li WD, Zhang G, Zhang M, Liu ZJ, Zhu KX, Liu QC, Zhang SE, Shen W, Zhang XF. Identification of unique transcriptomic signatures through integrated multispecies comparative analysis and WGCNA in bovine oocyte development. BMC Genomics. 2023;24(1):265.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Wei J, Zhang J, Wei J, Hu M, Chen X, Qin X, Chen J, Lei F, Qin Y. Identification of AGXT2, SHMT1, and ACO2 as important biomarkers of acute kidney injury by WGCNA. PLoS ONE. 2023;18(2):e281439.

    Article  Google Scholar 

  39. Wen P, Chen J, Zuo C, Gao X, Fujita M, Yang G. Proteome and glycoproteome analyses reveal the protein N-Linked glycosylation specificity of STT3A and STT3B. Cells-Basel. 2022;11(18):2775.

    Article  CAS  Google Scholar 

  40. Zhang Y, Sun L, Lei C, Li W, Han J, Zhang J, Zhang Y. A Sweet warning: mucin-type O-Glycans in Cancer. Cells-Basel. 2022;11(22):3666.

    Article  CAS  Google Scholar 

  41. Hofer DC, Pessentheiner AR, Pelzmann HJ, Schlager S, Madreiter-Sokolowski CT, Kolb D, Eichmann TO, Rechberger G, Bilban M, Graier WF, et al. Critical role of the peroxisomal protein PEX16 in white adipocyte development and lipid homeostasis. Biochim Biophys Acta Mol Cell Biol Lipids. 2017;1862(3):358–68.

    Article  CAS  PubMed  Google Scholar 

  42. Vuorinen EM, Rajala NK, Ihalainen TO, Kallioniemi A. Depletion of nuclear import protein karyopherin alpha 7 (KPNA7) induces mitotic defects and deformation of nuclei in cancer cells. BMC Cancer. 2018;18(1):325.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Ma D, Liu Q, Zhang M, Feng J, Li X, Zhou Y, Wang X. iTRAQ-based quantitative proteomics analysis of the spleen reveals innate immunity and cell death pathways associated with heat stress in broilers (Gallus gallus). J Proteom. 2019;196:11–21.

    Article  CAS  Google Scholar 

  44. Ma RS, Zhou ZL, Luo JW, Zhang H, Hou JF. The Ihh signal is essential for regulating proliferation and hypertrophy of cultured chicken chondrocytes. Comp Biochem Physiol B Biochem Mol Biol. 2013;166(2):117–22.

    Article  CAS  PubMed  Google Scholar 

  45. Debebe Z, Rathmell WK. Ror2 as a therapeutic target in cancer. Pharmacol Ther. 2015;150:143–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Zhang H, Yu JQ, Yang LL, Kramer LM, Zhang XY, Na W, Reecy JM, Li H. Identification of genome-wide SNP-SNP interactions associated with important traits in chicken. BMC Genomics. 2017;18(1):892.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Guo HY, Sun J, Su Y, Zhang Q, Hua R, Xia X, Fang D. Thrombotic microangiopathy led to acute kidney injury in an infant with hemophagocytic lymphohistiocytosis: a case report. Exp Ther Med. 2021;21(4):396.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We thank the professors and students of Guangdong Provincial Key Laboratory of Animal Molecular Design and Precise Breeding for their assistance and support in this study.

Funding

This work was supported by the Science and Technology Project of Qingcheng District, Qingyuan City (No.: 2020A01) and the Science and Technology Program of Qingyuan (No. 2021SJXM011).

Author information

Author notes

  1. Hao Qi and Zhidan Deng contributed equally to this work.

Authors and Affiliations

  1. Guangdong Provincial Key Laboratory of Animal Molecular Design and Precise Breeding, College of Life Science and Engineering, Foshan University, 528225, Foshan, Guangdong, China

    Hao Qi, Zhidan Deng, Fei Ye, Junwei Gou, Miaoxin Huang, Hai Xiang & Hua Li

Authors
  1. Hao Qi
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Zhidan Deng
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Fei Ye
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Junwei Gou
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Miaoxin Huang
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Hai Xiang
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Hua Li
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

Conceived and designed the experiments and wrote the paper: H.Q. and H. L., Performed the experiments: H.Q., D. Z. and F. Y., Analyzed the data: J. G., M. H. and H. X. All authors have read and approved the final version of the manuscript.

Corresponding author

Correspondence to Hua Li.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the Animal Care Committee of Foshan University (Foshan, People’s Republic of China; Approval ID: FOSU#056). All animal experiments were performed according to the National and International Guidelines for Animal Welfare.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Qi, H., Deng, Z., Ye, F. et al. Analysis of the differentially expressed genes in the combs and testes of Qingyuan partridge roosters at different developmental stages. BMC Genomics 25, 33 (2024). https://doi.org/10.1186/s12864-024-09960-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-024-09960-2

Keywords

BMC Genomics

ISSN: 1471-2164

Contact us