Skip to main content
Download PDF

Transcriptome profiling analysis of uterus during chicken laying periods

BMC Genomics volume 24, Article number: 433 (2023) Cite this article

Abstract

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.

Peer Review reports

Introduction

In the uterus, the egg absorbs fluid into albumin, which then calcifies on the eggshell membrane to form the eggshell. The formation of eggshells is a highly complex process that is precisely controlled by genetic and biological pathways in the uterus [1,2,3]. Eggshell formation requires the involvement of multiple ion transporter genes to supply the required ions and minerals [4]. ATP binding cassette subfamily C member 9 (ABCC9), inositol 1,4,5-trisphosphate receptor type 2 (ITPR2), potassium inwardly rectifying channel subfamily J member 8 (KCNJ8) and WNK lysine deficient protein kinase 1 (WNK1) modulate eggshell thickness by participating in ion transport [5]. Eggshells contain a variety of matrix proteins [6,7,8]. Three matrix protein genes ovocleidin-11, ovalbumin and ovocalyxin-32 were shown to affect eggshell thickness [9]. Milk fat globule EGF and factor V/VIII domain containing (MFGE8) and EGF like repeats and discoidin domains 3 (EDIL3) are two glycoproteins involved in the regulation of vesicle-mediated eggshell calcification [10]. These proteins and other organic components interact with minerals to form the eggshell in the uterus [11,12,13].

Eggshell quality decreases as hens age in the laying period [14,15,16,17]. Aged hens have a disrupted uterine structure with reduced glandular density, fibrosis and atrophy of the endometrium, probably due to continued egg-laying behavior, which leads to reduced protein synthesis, ion transport and immune defense of the uterus [17, 18]. The expression of matrix proteins and ion transporters (BPI fold containing family B member 3, ATPase sarcoplasmic/endoplasmic reticulum Ca2+ transporting 2, ovalbumin, carbonic anhydrase 2, transferrin, sodium channel epithelial 1 gamma subunit, ovocleidin 17) in the uterus of laying hens change with age [19,20,21]. Age-related dysregulation of gene expression in uterus may lead to deterioration of eggshell quality.

The rapid development of RNA sequencing (RNA-seq) has provided a powerful tool for the study of the transcriptome [22, 23]. RNA-seq has been widely used in the study of various species and traits in recent years [24,25,26,27,28,29]. There have been many studies on poultry uterus over the years [3, 30,31,32]. However, the avian uterine functional differences in different laying periods at the transcriptome level remains obscure. In this study, we screened differentially expressed genes and pathways in chicken uterus at different laying stages by RNA-seq and bioinformatics analysis, which laid a foundation for further studies on uterine function maintenance in chicken.

Materials and methods

Animals and sample collection

Twelve Nandan-Yao hens (Gallus gallus) used in this study, laying continuously 3 eggs, were sourced from Guangxi Guigang Gangfeng Agriculture and Husbandry Co., Ltd. Euthanasia was performed by cervical dislocation, with all effort made to minimize suffering. All chickens were caged and reared individually according to standard feeding management protocols. The uterine samples were obtained from three periods (four chickens per group), including early laying period (22 weeks old), peak laying period (31 weeks old), and late laying period (51 weeks old).

Total RNA extraction

The total RNA was extracted from the uterus with TRIzol reagent (Invitrogen Life Technologies, USA) following the manufacturer's instruction. Agarose gel electrophoresis was used to assess RNA concentration. RNA integrity was determined using UV–Vis Spectrophotometer Q5000 (Quawell, USA).

RNA sequencing and quality control

The cDNA libraries were constructed and sequenced following the manufacturer’s standard procedures on an Illumina HiSeq 2500 (Illumina, San Diego, CA, USA) in Novogene Bioinformatics Technology Co., Ltd., Beijing, China. Raw reads of FASTQ format were processed with the Trim Galore [33]. Low quality sequences including rate of N base > 10%, quality score < 20 and adaptor sequences were removed to generate clean data for downstream analysis. FastQC software was used to calculate sequence duplication levels, GC content and Q20 scores of clean data [34].

Bioinformatics analysis

Reference genome and gene model annotation files were downloaded from the ENSEMBL (http://ftp.ensembl.org/pub/release-102/gtf/gallus_gallus/, http://ftp.ensembl.org/pub/release-102/fasta/gallus_gallus/dna/). The clean reads were.mapped to the chicken reference genome by Hisat2v2.1.0 [35, 36]. The stringtiev2.1.1 was then used to annotate the transcripts [37]. Differential expression analysis was performed with the DESeq2 [38], a package in R software. Differentially expressed genes (DEGs) were screened by P-value < 0.05 and |log2foldchange|≥ 1. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) [39,40,41] analysis were conducted by the R package clusterProfiler 3.14.3 to identify critical pathways [42]. P-adjust < 0.05 denoted statistical significance. Interactions between DEGs were explored by STRING [42] database, with confidence score > 0.9 being valid [43]. We applied Cytoscape [44] (http://cytoscape.org/) to visualize the STRING analysis results and CytoHubba application to find key hub genes through the MCC algorithm [45]. Module detection (MCODE) in Cytoscape was used to screen the key modules in the protein–protein interaction network [46].

Validation of RNA-Seq

Six genes retinoic acid receptor responder 1 (RARRES1), dickkopf WNT signaling pathway inhibitor 3(DKK3), R-spondin 3(RSPO3), relaxin 3(RLN3), glypican 4(GPC4) and potassium inwardly rectifying channel subfamily J member 2(KCNJ2) related to uterine functional maintenance were selected for qRT-PCR analysis to validate RNA-seq data [9, 47,48,49,50]. Primer sequences of target and reference genes were shown in Supplemental table 1. RNA was reverse transcribed into cDNA using RT Reagent Kit (Takara, Dalian, China). qRT-PCR was performed by SYBR Green Supermix kit (Takara, Dalian, China) in Bio-RAD CFX96 Real Time Detection system. The qRT-PCR was performed according to the following program: the volume of the reaction mixture was 20 μl, with 2 μl of cDNA, 0.4 μl of primers, 10 μl of SYBR (Takara, Dalian, China), and 7.2 μl of RNA-free water. Then the QRT-PCR was run at 95 ◦C for 30 s, 95 ◦C for 5 s, 60 ◦C for 30 s, 35 cycles. The expression levels of genes were calculated relative to the expression of the β-actin using the 2−ΔΔCT method.

Results

Transcriptome data

As shown in Supplementary table 2, a total of 270,440,113 clean reads were acquired, ranging from 19,122,152 to 26,346,160. The clean reads were characterized by the average GC content of 51.95% and more than 97.12% of Q20, while more than 90.59% of clean reads were mapped to the reference genome.

Analysis of differentially expressed genes

There were 520, 706 and 736 DEGs identified by comparison of the uterus in W31 vs W22 group, W51 vs W31 group, W51 vs W22 group, respectively (Fig. 1A). Of the DEGs, 197 up-regulated and 323 down-regulated genes in the W31 vs W22 group (Fig. 1B), 479 up-regulated and 227 down-regulated genes in the W51 vs W31(Fig. 1C), 380 up-regulated and 356 down-regulated genes in the W51 vs W22(Fig. 1D). Furthermore, the intersection of the three comparison groups showed a total of 18 DEGs including sodium voltage-gated channel beta subunit 3 (SCN3B), avidin (AVD), family with sequence similarity 178, member B (FAM178B), solute carrier family 17 member 9 (SLC17A9), unc-13 homolog C (UNC13C), multiple EGF-like-domains 6-like(LOC424998), proopiomelanocortin (POMC), serum/glucocorticoid regulated kinase 1 (SGK1), ADAM metallopeptidase with thrombospondin type 1 motif 18(ADAMTS18), HEAT repeat-containing protein 8-like (HRC8L), transglutaminase 4 (TGM4), transmembrane protein 2-like (LOC415478), amphiregulin (AREG), ligand of Numb protein X 2-like (LOC422320), protein kinase cGMP-dependent type II (PRKG2) and so on (Fig. 1A).

Fig. 1
figure 1

Venn diagram and Volcano plot of DEGs. A represents the Venn diagram of intersection DEGs of W31 vs W22, W51 vs W31 and W51 vs W22. B, C and D represent the volcano plot of DEGs in W31 vs W22, W51 vs W31 and W51 vs W22, respectively

Full size image

GO and KEGG analysis for DEGs

In the W31 vs W22 group, GO analysis indicated that DEGs were enriched in 9 terms, including extracellular matrix, extracellular region part, extracellular region, collagen-containing extracellular matrix, extracellular space, glycosaminoglycan binding, heparin binding, sulfur compound binding and extracellular matrix structural constituent (Fig. 2A, Supplementary table 3).

Fig. 2
figure 2

GO enrichment analysis of DEGs. A, B and C represent GO enrichment analysis of DEGs in W31 vs W22, W51 vs W31 and W51 vs W22, respectively

Full size image

The results of functional enrichment analysis of the W51 vs W31 group showed that the extracellular region and extracellular region part were significantly enriched in GO analysis (Fig. 2B, Supplementary table 3), and neuroactive ligand-receptor interaction were significantly enriched in KEGG analysis (Fig. 3A, Supplementary table 5).

Fig. 3
figure 3

KEGG enrichment analysis of DEGs. A and B represent KEGG analysis of DEGs in W51 vs W31, W51 vs W22, respectively

Full size image

In the W51 vs W22 group, GO analysis of DEGs were significantly enriched in 9 terms, containing collagen-containing extracellular matrix, extracellular matrix, collagen trimer, extracellular region, extracellular region part, supramolecular complex, supramolecular polymer, supramolecular fiber and extracellular space (Fig. 2C, Supplementary table 6). And ECM-receptor interaction was significantly enriched in KEGG analysis (Fig. 3B, Supplementary table 7).

Integration of PPI network

In the W31 vs W22 group, the top ten genes were cell division cycle 45(CDC45), cell division cycle 6(CDC6), ribonucleotide reductase regulatory subunit M2(RRM2), minichromosome maintenance complex component 5(MCM5), structural maintenance of chromosomes 4(SMC4), structural maintenance of chromosomes 2 (SMC2), WD repeat and HMG-box DNA binding protein 1(WDHD1), maternal embryonic leucine zipper kinase(MELK), non-SMC condensin II complex subunit G2(NCAPG2), and thymidine kinase 1(TK1) (Fig. 4A). Module analysis showed modules 1, 2 and 3 included 15, 10 and 7 genes with module scores of 11.091, 6.333 and 6.25, respectively(Fig. 4B-D).

Fig. 4
figure 4

Top ten genes and three protein–protein interaction (PPI) hub network modules of W31 vs W22 group. A represents the top ten genes of W31 vs W22. The darker the color of the node represents the higher ranking of the gene. B, C and D represent the top three modules, respectively. Red and blue nodes indicate up- and down-regulation of genes, respectively. Larger nodes represent more interaction relationships

Full size image

The top ten genes were erb-b2 receptor tyrosine kinase 4(ERBB4), amphiregulin (AREGB), epiregulin (EREG), fibroblast growth factor 1(FGF1), betacellulin (BTC), phosphoinositide-3-kinase adaptor protein 1(PIK3AP1), neuregulin 2(NRG2), kinesin family member 11(KIF11), marker of proliferation Ki-67(MKI67) and anillin actin binding protein (ANLN) in the W51 vs W31 group (Fig. 5A). Module analysis indicated modules 1, 2 and 3 contained 7, 6 and 5 genes with module scores of 7.143, 6.444 and 6, respectively (Fig. 5B-D).

Fig. 5
figure 5

Top ten genes and three protein–protein interaction (PPI) hub network modules of W51 vs W31 group. A represents the top ten genes of W51 vs W31. The darker the color of the node represents the higher ranking of the gene. B, C and D represent the top three modules, respectively. Red and blue nodes indicate up- and down-regulation of genes, respectively. Larger nodes represent more interaction relationships

Full size image

In the W51 vs W22 group, the top ten genes contained collagen type I alpha 1 chain(COL1A1), collagen type I alpha 2 chain(COL1A2), collagen type VI alpha 1 chain(COL6A1), collagen type VI alpha 3 chain(COL6A3), periostin (POSTN), collagen type V alpha 1 chain(COL5A1), lumican (LUM), collagen type VI alpha 2 chain(COL6A2), secreted protein acidic and cysteine rich (SPARC), and fibrillin 1(FBN1) (Fig. 6A). Module analysis revealed modules 1, 2 and 3 consisted of 12, 7 and 9 genes with module scores of 6, 4 and 3, respectively (Fig. 6B-D).

Fig. 6
figure 6

Top ten genes and three protein–protein interaction (PPI) hub network modules of W51 vs W22 group. A represents the top ten genes of W51 vs W22. The darker the color of the node represents the higher ranking of the gene. B, C and D represent the top three modules, respectively. Red and blue nodes indicate up- and down-regulation of genes, respectively. Larger nodes represent more interaction relationships

Full size image

Dynamic transcriptional profiling of the uterus

To study the dynamic transcriptome profile of chicken uterus during the laying period (22W, 31W, 51W), eight major expression patterns were determined using the fuzzy c-means algorithm in the Mfuzz package, and the results are shown in Fig. 7.

Fig. 7
figure 7

Dynamic transcriptional profiling of the uterus. Note: Left: Time series expression spectrum of eight clusters. The red line indicates the average gene expression level. Middle plot: Expression heatmap plotted with TPM values for each gene in eight clusters. Right: Top five significant gene GO enrichment in each cluster

Full size image

Cluster 1 (3057 genes), cluster 3 (2412 genes) and cluster 8 (2039 genes) were specifically highly expressed at 22W. Cluster 1 significantly enriched in 147 GO terms. which mainly includes collagen-containing extracellular matrix, regulation of cell size, collagen trimer, and external encapsulating structure, extracellular matrix, etc. Cluster 3 significantly enriched in 34 GO terms, which mainly includes supramolecular polymer, supramolecular fiber, microtubule-based transport, supramolecular complex, contractile fiber, etc. Cluster 8 significantly enriched in 76 GO terms, which mainly includes epidermis development, epidermal cell differentiation, epithelial cell differentiation and carbohydrate derivative catabolic process, carbohydrate derivative metabolic process, etc. KEGG enrichment analysis showed that these genes were significantly enriched in cytokine-cytokine receptor interaction signaling pathway.

Cluster 2 (2288 genes) and cluster 4 (2931 genes) were specifically highly expressed at 31W. Cluster 2 significantly enriched in 89 GO terms, which mainly includes transferase complex, catalytic complex, glycolytic process, ATP generation from ADP, chromosome, etc. KEGG enrichment analysis showed that these genes were significantly enriched in the proteasome signaling pathway. Cluster 4 significantly enriched in 14 GO terms, which mainly includes endomembrane system, organelle membrane, peptidase complex, bounding membrane of organelle, basal part of cell and so on. KEGG enrichment analysis showed that these genes in cluster 4 were significantly enriched in 2 signaling pathways, including oxidative phosphorylation and proteasome. Cluster 6 (2047 genes) was highly expressed at 31W and 51W. KEGG enrichment analysis showed that these genes were significantly enriched in the ribosome signaling pathway.

Cluster 5 (2152 genes) and cluster 7 (3380 genes) were specifically highly expressed at 51W. Cluster 5 significantly enriched in 9 GO terms, which mainly includes endoplasmic reticulum lumen, centrosome, microtubule organizing center, protein folding chaperone, centriole, etc. Cluster 7 significantly enriched in 175 GO terms, which mainly includes double-stranded DNA binding, regulation of gene expression, DNA-binding transcription factor activity, RNA polymerase II-specific, mesonephros development, mesonephric epithelium development, etc. KEGG enrichment analysis showed that these genes in cluster 7 were significantly enriched in 2 signaling pathways, including MAPK signaling pathway and neuroactive ligand-receptor interaction.

Validation of RNA-seq

We chose 6 genes (RARRES1, DKK3, RSPO3, RLN3, GPC4, KCNJ2) related to uterine functional maintenance for qRT-PCR analysis to validate RNA-seq data. The analysis data showed that the trends of six genes were corroborated with RNA-seq analysis results(Fig. 8). The results indicated that the results of RNA-seq analysis were accurate and reliable.

Fig. 8
figure 8

RNA-Seq validation using QRT-PCR. Six DEGs (RARRES1, DKK3, RSPO3, RLN3, GPC4, KCNJ2) were used to validate the accuracy of RNA sequencing. N = 4

Full size image

Discussion

The structure and function of the uterus have a critical impact on eggshell quality. Ion transport and eggshell matrix protein secretion in the uterus are regulated by numerous genes [2, 51, 52]. In this study, high-throughput transcriptome analysis was used to investigate the differential gene expression profiles in three different periods of uterus. The identified differentially expressed genes and pathways generated by this study are a rich resource that can be used to deepen our understanding of functional maintenance of the chicken uterus.

DEGs were significantly enriched in the extracellular matrix, extracellular region, extracellular region part, ECM receptor interaction, extracellular matrix structural constituent, collagen-containing extracellular matrix and collagen trimer. The extracellular matrix (ECM) is involved in uterine function [53]. The ECM consists of collagen, proteoglycans and glycoproteins such as fibronectin and laminin, which is a complex mixture of structural proteins that act as a support structure to maintain three-dimensional tissue architecture. The ECM also regulates cellular responses through ligand-integrin interactions, participating in cell adhesion, survival, proliferation, differentiation and migration [54,55,56].

Protein network interaction analysis of DEGs identified 15 genes related to uterine functional maintenance including collagen type V alpha 2 chain(COL5A2), COL5A1, COL1A1, COL1A2, fibronectin 1(FN1), lysyl oxidase (LOX), thrombospondin 2(THBS2), periostin (POSTN), matrix metallopeptidase 13(MMP13), VANGL planar cell polarity protein 2(VANGL2), RAD54 homolog B (RAD54B), secreted phosphoprotein 1(SPP1), Syndecan-1 (SDC1), Betacellulin (BTC), and angiopoietin like 3(ANGPTL3).

In our results, FN1, LOX, THBS2, COL1A1, COL1A2, COL5A1, COL5A2 and POSTN were specifically highly expressed at 22W.Studies have shown that collagens are major structural protein components of the extracellular matrix. Collagen type V (COLV) acts in conjunction with collagen type I (COL1A1, COL1A2) and collagen type III to regulate fiber formation and deposition in the ECM [57]. There are three different isomers of COLV: COL5A1, COL5A2 and COL5A3 [58]. FN1 is a glycoprotein in the extracellular matrix in plasma and cell surface. It mediates cellular interactions with the ECM and takes part in cell adhesion, migration, growth, differentiation and metastasis [59, 60]. FN1 plays an important role in the endometrium, influencing placentation, blastocyst adhesion and implantation [61,62,63]. LOX is a copper-containing amine oxidase that catalyzes cross-linking of collagen and elastin in the ECM. The inhibition of LOX activity decreases endometrial stromal cell migration and embryo adhesion [64]. THBS2 is a homologous glycoprotein that is linked to the ECM disulfide bond. It is involved in adhesion, migration and proliferation of cell, interactions of cell-to-cell and cell–matrix, and angiogenesis [65, 66]. Our results showed that COL1A1, COL1A2, THBS2, LOX, FN1 were significantly enriched in terms such as extracellular matrix, collagen-containing extracellular matrix, external encapsulating structure, extracellular matrix, extracellular region, extracellular space, etc. POSTN is a matricellular protein that can combine with structural components of the extracellular matrix which is abundantly expressed in collagen-rich fibrous connective tissue [67].

In our results, MMP13, VANGL2 and RAD54B were specifically highly expressed at 31W. Matrix metallopeptidases (MMPs) and their inhibitors (TIMPs) regulate COL turnover and ECM remodeling [68]. Matrix metalloproteinases belong to the zinc proteases family that hydrolyze components of the extracellular matrix [69]. MMP13 participates in the degradation of the collagen network of the ECM [69, 70]. VANGL2 is the vertebrate ortholog of core Planar cell polarity (PCP) components. PCP refers to the capacity of a tissue, typically, but not exclusively, an epithelium, to transmit directional information across the tissue plane such that its cellular constituents can differentiate, divide or move in a coordinated manner and along a common axis, generally orthogonal to the apical-basal axis [71]. Research has shown that VANGL2 gene play a crucial role in regulating planar cell polarity and convergent extension movements [72]. RAD54B involved in the initial stages of homologous recombination. Vertebrate cells lacking RAD54B exhibit reduced levels of gene conversion [73, 74]. RAD54B is involved in cell cycle regulation after DNA damage and participates in homologous recombinational repair, which ensures the precise repair of the most deleterious DNA lesions, double-stranded breaks [75]. The Peak laying period is accompanied by degradation of the collagen network of the ECM, homologous recombination, and cellular information transfer.

In our results, SPP1, SDC1, BTC and ANGPTL3 were specifically highly expressed at 51W. SPP1 is a secretory protein present in the extracellular matrix. SPP1 protein is localized to the endometrium and participates in the formation and calcification of the eggshell [76,77,78]. Considering the reported inhibitory function of osteopontin in other mineralized systems, together with its main occurrence in the non-mineralized parts of the eggshell and at the outermost part of the shell, suggests that this molecule could be part of the mechanism regulating the eggshell calcification [79]. SDC1 play an important role in the placenta, in trophoblast development, and in pregnancy complications, highlighting one of the most important members of this family. Proteoglycan macromolecules play key roles in several physiological processes (e.g., adhesion, proliferation, migration, invasion, angiogenesis, and apoptosis), all of which are important for placentation and healthy pregnancy [80]. BTC was shown to be expressed exclusively at the sites of blastocyst apposition in the mouse uterus, suggesting an involvement in embryo implantation [81]. Retarded (but not impaired) embryo implantation was observed in BTC-overexpressing transgenic mice [82]. ANGPTL3 was an endometrial gene associated with angiogenesis, increasing uterine blood flow during early pregnancy to enhance the availability of micronutrients for transport into the uterine lumen [83]. These genes may interact to regulate the functional maintenance of the uterus.

Conclusions

The results of the present study suggest show that uterine tissues of hens at different laying stages differ in their transcriptomes. A series of key genes and pathways involved in the maintenance of uterine function were identified using RNA-seq and bioinformatics analysis. These key genes (FN1, LOX, THBS2, COL1A1, COL1A2, COL5A1, COL5A2, POSTN, MMP13, VANGL2, RAD54B, SPP1, SDC1, BTC, ANGPTL3) may take part in uterine functional maintenance through the extracellular matrix, extracellular region, extracellular region part, ECM receptor interaction, extracellular matrix structural constituent, collagen-containing extracellular matrix and collagen trimer. We have generated a rich resource, which will provide insight into chicken uterine functional maintenance.

Availability of data and materials

The RNA sequencing data used during the current study are available from the NCBI (accession number: PRJNA849016 and the link of the website: http://www.ncbi.nlm.nih.gov/bioproject/849016).

Abbreviations

RNA-Seq:

RNA sequencing

DEGs:

Differentially expressed genes

GO:

Gene Ontology

KEGG:

Kyoto Encyclopedia of Genes and Genomes

qRT-PCR:

Quantitative Real-time PCR

References

  1. Jeong W, Lim W, Kim J, Ahn SE, Lee HC, Jeong JW, Han JY, Song G, Bazer FW. Cell-specific and temporal aspects of gene expression in the chicken oviduct at different stages of the laying cycle. Biol Reprod. 2012;86(6):172.

    PubMed  Google Scholar 

  2. Khan S, Wu SB, Roberts J. RNA-sequencing analysis of shell gland shows differences in gene expression profile at two time-points of eggshell formation in laying chickens. BMC Genomics. 2019;20(1):89.

    PubMed  PubMed Central  Google Scholar 

  3. Yin Z, Lian L, Zhu F, Zhang ZH, Hincke M, Yang N, Hou ZC. The transcriptome landscapes of ovary and three oviduct segments during chicken (Gallus gallus) egg formation. Genomics. 2020;112(1):243–51.

    CAS  PubMed  Google Scholar 

  4. Brionne A, Nys Y, Hennequet-Antier C, Gautron J. Hen uterine gene expression profiling during eggshell formation reveals putative proteins involved in the supply of minerals or in the shell mineralization process. BMC Genomics. 2014;15:220.

    PubMed  PubMed Central  Google Scholar 

  5. Duan Z, Sun C, Shen M, Wang K, Yang N, Zheng J, Xu G. Genetic architecture dissection by genome-wide association analysis reveals avian eggshell ultrastructure traits. Sci Rep. 2016;6:28836.

    CAS  PubMed  PubMed Central  Google Scholar 

  6. Gautron J, Guyot N, Brionne A, R´ehault-Godbert S. Bioactive minor egg components; 2019.

  7. Marie P, Labas V, Brionne A, Harichaux G, Hennequet-Antier C, Nys Y, Gautron J. Quantitative proteomics and bioinformatic analysis provide new insight into protein function during avian eggshell biomineralization. J Proteomics. 2015;113:178–93.

    CAS  PubMed  Google Scholar 

  8. Rose-Martel M, Smiley S, Hincke MT. Novel identification of matrix proteins involved in calcitic biomineralization. J Proteomics. 2015;116:81–96.

    CAS  PubMed  Google Scholar 

  9. Dunn IC, Joseph NT, Bain M, Edmond A, Wilson PW, Milona P, Nys Y, Gautron J, Schmutz M, Preisinger R, et al. Polymorphisms in eggshell organic matrix genes are associated with eggshell quality measurements in pedigree Rhode Island Red hens. Anim Genet. 2009;40(1):110–4.

    CAS  PubMed  Google Scholar 

  10. Stapane L, Le Roy N, Hincke MT, Gautron J. The glycoproteins EDIL3 and MFGE8 regulate vesicle-mediated eggshell calcification in a new model for avian biomineralization. J Biol Chem. 2019;294(40):14526–45.

    CAS  PubMed  PubMed Central  Google Scholar 

  11. Gautron J, Hincke MT, Nys Y. Precursor matrix proteins in the uterine fluid change with stages of eggshell formation in hens. Connect Tissue Res. 1997;36(3):195–210.

    CAS  PubMed  Google Scholar 

  12. Weiner S. Biomineralization: a structural perspective. J Struct Biol. 2008;163(3):229–34.

    CAS  PubMed  Google Scholar 

  13. Nys Y, Guyot N. Improving the Safety and Quality of Eggs and Egg Products. In: Egg formation and chemistry. Edited by Van Immerseel F, Nys Y, Bain M. Sawston, Cambridge: Woodhead Publishing; 2011. p. 400-408.

  14. Bain MM, Nys Y, Dunn IC. Increasing persistency in lay and stabilising egg quality in longer laying cycles. What are the challenges. Br Poult Sci. 2016;57:330–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  15. Sah N, Kuehu DL, Khadka VS, Deng Y, Peplowska K, Jha R, Mishra B. RNA sequencing-based analysis of the laying hen uterus revealed the novel genes and biological pathways involved in the eggshell biomineralization. Sci Rep. 2018;8(1):16853.

    PubMed  PubMed Central  Google Scholar 

  16. Moln’ar A, Maertens L, Ampe B, Buyse J, Kempen I, Zoons J, Delezie E. Changes in egg quality traits during the last phase of production: is there potential for an extended laying cycle. Br Poult Sci. 2016;57:842–7.

    CAS  Google Scholar 

  17. Park JA, Sohn SH. The Influence of Hen Aging on Eggshell Ultrastructure and Shell Mineral Components. Korean J Food Sci Anim Resour. 2018;38(5):1080–91.

    PubMed  PubMed Central  Google Scholar 

  18. Anna W, Yvonne R, Helena W, Lena H. Age-related changes in the shell gland and duodenum in relation to shell quality and bone strength in commercial laying hen hybrids. Acta Vet Scand. 2019;61(1):14.

    Google Scholar 

  19. Park JA, Cho EJ, Park JY, Sohn SH. Histological change of uterus endometrium and expression of the eggshell-related genes according to hen age. Korean J Poult Sci. 2017;44(1):19–28.

    Google Scholar 

  20. Veronika G, Monika P, Ivan M, Eva T. Concentration of egg white antimicrobial and immunomodulatory proteins is related to eggshell pigmentation across traditional chicken breeds. Poult Sci. 2019;98(12):6931–41.

    Google Scholar 

  21. Panheleux M, Nys Y, Williams J, Gautron J, Boldicke T, Hincke MT. Extraction and quantification by ELISA of eggshell organic matrix proteins (ovocleidin-17, ovalbumin, ovotransferrin) in shell from young and old hens. Poult Sci. 2000;79(4):580–8.

    CAS  PubMed  Google Scholar 

  22. Emrich SJ, Barbazuk WB, Li L, Schnable PS. Gene discovery and annotation using LCM-454 transcriptome sequencing. Genome Res. 2007;17(1):69–73.

    CAS  PubMed  PubMed Central  Google Scholar 

  23. Lister R, O’Malley RC, Tonti-Filippini J, Gregory BD, Berry CC, Millar AH, Ecker JR. Highly integrated single-base resolution maps of the epigenome in Arabidopsis. Cell. 2008;133(3):523–36.

    CAS  PubMed  PubMed Central  Google Scholar 

  24. Jin L, Zhao L, Hu S, Long K, Liu P, Liu R, Zhou X, Wang Y, Huang Z, Lin X, et al. Transcriptional Differences of Coding and Non-Coding Genes Related to the Absence of Melanocyte in Skins of Bama Pig. Genes. 2019;11(1):47.

    PubMed  PubMed Central  Google Scholar 

  25. Zhou S, Ma Y, Zhao D, Mi Y, Zhang C. Transcriptome profiling analysis of underlying regulation of growing follicle development in the chicken. Poult Sci. 2020;99(6):2861–72.

    CAS  PubMed  PubMed Central  Google Scholar 

  26. Zou X, Wang J, Qu H, Lv XH, Shu DM, Wang Y, Ji J, He YH, Luo CL, Liu DW. Comprehensive analysis of miRNAs, lncRNAs, and mRNAs reveals potential players of sexually dimorphic and left-right asymmetry in chicken gonad during gonadal differentiation. Poult Sci. 2020;99(5):2696–707.

    CAS  PubMed  PubMed Central  Google Scholar 

  27. Jia B, Liu Y, Li Q, Zhang J, Ge C, Wang G, Chen G, Liu D, Yang F. Altered miRNA and mRNA Expression in Sika Deer Skeletal Muscle with Age. Genes. 2020;11(2):172.

    CAS  PubMed  PubMed Central  Google Scholar 

  28. Li B, Zhang K, Ye Y, Xing J, Wu Y, Ma Y, Li Y. Effects of Castration on miRNA, lncRNA, and mRNA Profiles in Mice Thymus. Genes. 2020;11(2):147.

    PubMed  PubMed Central  Google Scholar 

  29. Li TT, Wang SY, Wu RM, Zhou XY, Zhu DH, Zhang Y. Identification of long non-protein coding RNAs in chicken skeletal muscle using next generation sequencing. Genomics. 2012;99(5):292–8.

    CAS  PubMed  Google Scholar 

  30. Khogali MK, Wen K, Jauregui D, Liu L, Zhao M, Gong D, Geng T. Uterine structure and function contributes to the formation of the sandpaper-shelled eggs in laying hens. Anim Reprod Sci. 2021;232: 106826.

    CAS  PubMed  Google Scholar 

  31. Cui Z, Zhang Z, Amevor FK, Du X, Li L, Tian Y, Kang X, Shu G, Zhu Q, Wang Y, et al. Circadian miR-449c-5p regulates uterine Ca(2+) transport during eggshell calcification in chickens. BMC Genomics. 2021;22(1):764.

    CAS  PubMed  PubMed Central  Google Scholar 

  32. Zhang F, Yin ZT, Zhang JF, Zhu F, Hincke M, Yang N, Hou ZC. Integrating transcriptome, proteome and QTL data to discover functionally important genes for duck eggshell and albumen formation. Genomics. 2020;112(5):3687–95.

    CAS  PubMed  Google Scholar 

  33. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20.

    CAS  PubMed  PubMed Central  Google Scholar 

  34. Schmieder R, Edwards R. Quality control and preprocessing of metagenomic datasets. Bioinformatics. 2011;27(6):863–4.

    CAS  PubMed  PubMed Central  Google Scholar 

  35. Pertea M, Kim D, Pertea GM, Leek JT, Salzberg SL. Transcript-level expression analysis of RNA-seq experiments with HISAT. StringTie and Ballgown Nature protocols. 2016;11(9):1650–67.

    CAS  PubMed  Google Scholar 

  36. Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12(4):357–60.

    CAS  PubMed  PubMed Central  Google Scholar 

  37. Pertea M, Pertea GM, Antonescu CM, Chang TC, Mendell JT, Salzberg SL. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat Biotechnol. 2015;33(3):290–5.

    CAS  PubMed  PubMed Central  Google Scholar 

  38. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.

    PubMed  PubMed Central  Google Scholar 

  39. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.

    CAS  PubMed  PubMed Central  Google Scholar 

  40. Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 2019;28(11):1947–51.

    CAS  PubMed  PubMed Central  Google Scholar 

  41. Kanehisa M, Furumichi M, Sato Y, Kawashima M, Ishiguro-Watanabe M. KEGG for taxonomy-based analysis of pathways and genomes. Nucleic Acids Res. 2023;51(D1):D587-d592.

    CAS  PubMed  Google Scholar 

  42. Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;5(16):284–7.

    Google Scholar 

  43. Franceschini A, Szklarczyk D, Frankild S, Kuhn M, Simonovic M, Roth A, Lin J, Minguez P, Bork P, von Mering C. STRINGv9.1: protein-protein interaction networks, with increased coverage andintegration. Nucleic Acids Res. 2013;41(D1):D808–15.

    CAS  PubMed  Google Scholar 

  44. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.

    CAS  PubMed  PubMed Central  Google Scholar 

  45. Liu Z, Meng J, Li X, Zhu F, Liu T, Wu G, Zhang L. Identification of Hub Genes and Key Pathways Associated with Two Subtypes of Diffuse Large B-Cell Lymphoma Based on Gene Expression Profiling via Integrated Bioinformatics. Biomed Res Int. 2018;2018:3574534.

    PubMed  PubMed Central  Google Scholar 

  46. Bader GD, Hogue CW. An automated method for finding molecular complexes in large protein interaction networks. BMC Bioinformatics. 2003;4(4):2.

    PubMed  PubMed Central  Google Scholar 

  47. Sah N, Kuehu DL, Khadka VS, Deng Y, Jha R, Wasti S, Mishra B. RNA sequencing-based analysis of the magnum tissues revealed the novel genes and biological pathways involved in the egg-white formation in the laying hen. BMC Genomics. 2021;22(1):318.

    CAS  PubMed  PubMed Central  Google Scholar 

  48. Lavelin I, Meiri N, Einat M, Genina O, Pines M. Mechanical strain regulation of the chicken glypican-4 gene expression in the avian eggshell gland. Am J Physiol Regul Integr Comp Physiol. 2002;283(4):R853–61.

    PubMed  Google Scholar 

  49. Marie P, Labas V, Brionne A, Harichaux G, Hennequet-Antier C, Rodriguez-Navarro AB, Nys Y, Gautron J. Quantitative proteomics provides new insights into chicken eggshell matrix protein functions during the primary events of mineralisation and the active calcification phase. J Proteomics. 2015;126:140–54.

    CAS  PubMed  Google Scholar 

  50. Jonchère V, Brionne A, Gautron J, Nys Y. Identification of uterine ion transporters for mineralisation precursors of the avian eggshell. BMC Physiol. 2012;12:10.

    PubMed  PubMed Central  Google Scholar 

  51. Jonch’ere V, R’ehault-Godbert S, Hennequet-Antier C, Cabau C, Sibut V, Cogburn LA, Nys Y, Gautron J. Gene expression profiling to identify eggshell proteins involved in physical defense of the chicken egg. BMC Genomic. 2010;11:57.

    Google Scholar 

  52. Jonch’ere V, Brionne A, Gautron J, Nys Y. Identification of uterine ion transporters for mineralisation precursors of the avian eggshell. BMC Physiol. 2012;12:10.

    CAS  Google Scholar 

  53. O’Connor BB, Pope BD, Peters MM, Ris-Stalpers C, Parker KK. The role of extracellular matrix in normal and pathological pregnancy: Future applications of microphysiological systems in reproductive medicine. Exp Biol Med (Maywood). 2020;245(13):1163–74.

    CAS  PubMed  Google Scholar 

  54. Järveläinen H, Sainio A, Koulu M, Wight TN, Penttinen R. Extracellular matrix molecules: potential targets in pharmacotherapy. Pharmacol Rev. 2009;61(2):198–223.

    PubMed  PubMed Central  Google Scholar 

  55. Malik M, Norian J, McCarthy-Keith D, Britten J, Catherino WH. Why leiomyomas are called fibroids: the central role of extracellular matrix in symptomatic women. Semin Reprod Med. 2010;28(3):169–79.

    PubMed  Google Scholar 

  56. Leppert PC, Jayes FL, Segars JH. The extracellular matrix contributes to mechanotransduction in uterine fibroids. Obstet Gynecol Int. 2014;2014: 783289.

    PubMed  PubMed Central  Google Scholar 

  57. Mak KM, Png CY, Lee DJ. Type V Collagen in Health, Disease, and Fibrosis. Anat Rec (Hoboken). 2016;299(5):613–29.

    CAS  PubMed  Google Scholar 

  58. Roulet M, Välkkilä M, Chanut-Delalande H, Hämäläinen ER, Kessler E, Ala-Kokko L, Männikkö M, Bonod-Bidaud C, Ruggiero F. The collagen V homotrimer [alpha1(V)](3) production is unexpectedly favored over the heterotrimer [alpha1(V)](2)alpha2(V) in recombinant expression systems. J biomed biotechnol. 2010;2010:376927.

    PubMed  PubMed Central  Google Scholar 

  59. Pankov R, Yamada KM. Fibronectin at a glance. J Cell Sci. 2002;115(20):3861–3.

    CAS  PubMed  Google Scholar 

  60. Soikkeli J, Podlasz P, Yin M, Nummela P, Jahkola T, Virolainen S, Krogerus L, Heikkilä P, von Smitten K, Saksela O, et al. Metastatic outgrowth encompasses COL-I, FN1, and POSTN up-regulation and assembly to fibrillar networks regulating cell adhesion, migration, and growth. Am J Pathol. 2010;177(1):387–403.

    CAS  PubMed  PubMed Central  Google Scholar 

  61. Bronson RA, Fusi FM. Integrins and human reproduction. Mol Hum Reprod. 1996;2(3):153–68.

    CAS  PubMed  Google Scholar 

  62. Wang J, Armant DR. Integrin-mediated adhesion and signaling during blastocyst implantation. Cells Tissues Organs. 2002;172(3):190–201.

    CAS  PubMed  Google Scholar 

  63. Wang JML, Armant DR. Trophoblast adhesion of the peri-implantation mouse blastocyst is regulated by integrin signaling that targets phospholipase C. Dev Biol. 2007;302:143–53.

    CAS  PubMed  Google Scholar 

  64. Li SY, Yan JQ, Song Z, Liu YF, Song MJ, Qin JW, Yang ZM, Liang XH. Molecular characterization of lysyl oxidase-mediated extracellular matrix remodeling during mouse decidualization. FEBS Lett. 2017;591(10):1394–407.

    CAS  PubMed  Google Scholar 

  65. Oguri M, Kato K, Yokoi K, Watanabe S, Metoki N, Yoshida H, Satoh K, Aoyagi Y, Nishigaki Y, Yoshida H, et al. Association of polymorphisms of THBS2 and HSPA8 with hypertension in Japanese individuals with chronic kidney disease. Mol Med Rep. 2009;2(2):205–11.

    CAS  PubMed  Google Scholar 

  66. Laherty CD, O’Rourke K, Wolf FW, Katz R, Seldin MF, Dixit VM. Characterization of mouse thrombospondin 2 sequence and expression during cell growth and development. J Biol Chem. 1992;267(5):3274–81.

    CAS  PubMed  Google Scholar 

  67. Norris RA, Damon B, Mironov V, Kasyanov V, Ramamurthi A, Moreno-Rodriguez R, Trusk T, Potts JD, Goodwin RL, Davis J, et al. Periostin regulates collagen fibrillogenesis and the biomechanical properties of connective tissues. J Cell Biochem. 2007;101(3):695–711.

    CAS  PubMed  PubMed Central  Google Scholar 

  68. Szóstek-Mioduchowska AZ, Baclawska A, Okuda K, Skarzynski DJ. Effect of proinflammatory cytokines on endometrial collagen and metallopeptidase expression during the course of equine endometrosis. Cytokine. 2019;123: 154767.

    PubMed  Google Scholar 

  69. Dollery CM, McEwan JR, Henney AM. Matrix metalloproteinases and cardiovascular disease. Circ Res. 1995;77(5):863–8.

    CAS  PubMed  Google Scholar 

  70. T SS. Regulation of proteolytic activity in tissues. Crit Rev Biochem Mol Biol. 1994;29:315–83.

    Google Scholar 

  71. Bailly E, Walton A, Borg JP. The planar cell polarity Vangl2 protein: From genetics to cellular and molecular functions. Semin Cell Dev Biol. 2018;81:62–70.

    CAS  PubMed  Google Scholar 

  72. Torban E, Kor C, Gros P. Van Gogh-like2 (Strabismus) and its role in planar cell polarity and convergent extension in vertebrates. Trends Genet. 2004;20(11):570–7.

    CAS  PubMed  Google Scholar 

  73. Miyagawa K, Tsuruga T, Kinomura A, Usui K, Katsura M, Tashiro S, Mishima H, Tanaka K. A role for RAD54B in homologous recombination in human cells. Embo j. 2002;21(1–2):175–80.

    CAS  PubMed  PubMed Central  Google Scholar 

  74. Sonoda E, Takata M, Yamashita YM, Morrison C, Takeda S. Homologous DNA recombination in vertebrate cells. Proc Natl Acad Sci U S A. 2001;98(15):8388–94.

    CAS  PubMed  PubMed Central  Google Scholar 

  75. McAndrew EN, McManus KJ. The enigmatic oncogene and tumor suppressor-like properties of RAD54B: Insights into genome instability and cancer. Genes Chromosomes Cancer. 2017;56(7):513–23.

    CAS  PubMed  Google Scholar 

  76. Chaen T, Konno T, Egashira M, Bai R, Nomura N, Nomura S, Hirota Y, Sakurai T, Imakawa K. Estrogen-dependent uterine secretion of osteopontin activates blastocyst adhesion competence. PLoS ONE. 2012;7(11): e48933.

    CAS  PubMed  PubMed Central  Google Scholar 

  77. Qi QR, Xie QZ, Liu XL, Zhou Y. Osteopontin is expressed in the mouse uterus during early pregnancy and promotes mouse blastocyst attachment and invasion in vitro. PLoS ONE. 2014;9(8): e104955.

    PubMed  PubMed Central  Google Scholar 

  78. Pines M, Bar VK A. Involvement of osteopontin in egg shell formation in the laying chicken. Matrix Biol. 1995;14:765–71.

    CAS  PubMed  Google Scholar 

  79. Fernandez MS, Escobar C, Lavelin I, Pines M, Arias JL. Localization of osteopontin in oviduct tissue and eggshell during different stages of the avian egg laying cycle. J Struct Biol. 2003;143(3):171–80.

    CAS  PubMed  Google Scholar 

  80. Oravecz O, Balogh A, Romero R, Xu Y, Juhasz K, Gelencser Z, Xu Z, Bhatti G, Pique-Regi R, Peterfia B, et al. Proteoglycans. Systems-Level Insight into Their Expression in Healthy and Diseased Placentas. Int J Mol Sci. 2022;23(10):5798.

  81. Das SK, Das N, Wang J, Lim H, Schryver B, Plowman GD, Dey SK. Expression of betacellulin and epiregulin genes in the mouse uterus temporally by the blastocyst solely at the site of its apposition is coincident with the “window” of implantation. Dev Biol. 1997;190(2):178–90.

    CAS  PubMed  Google Scholar 

  82. Gratao AA, Dahlhoff M, Sinowatz F, Wolf E, Schneider MR. Betacellulin overexpression in the mouse ovary leads to MAPK3/MAPK1 hyperactivation and reduces litter size by impairing fertilization. Biol Reprod. 2008;78(1):43–52.

    CAS  PubMed  Google Scholar 

  83. Dorniak P, Bazer FW, Wu G, Spencer TE. Conceptus-derived prostaglandins regulate endometrial function in sheep. Biol Reprod. 2012;87(1):9 1–7.

    PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was supported in part by the Science and Technology Key Project of Guangxi (GK AB21075010) and the Science and Technology Major Project of Guangxi (GK AA17204027).

Author information

Author notes

  1. Tiantian Sun, and Cong Xiao contributed equally to this work.

Authors and Affiliations

  1. College of Animal Science and Technology, Guangxi University, Nanning, 530004, China

    Tiantian Sun, Cong Xiao, Zhuliang Yang, Jixian Deng & Xiurong Yang

  2. Guangxi Key Laboratory of Animal Breeding, Disease Control and Prevention, Guangxi University, Nanning, 530004, China

    Xiurong Yang

Authors
  1. Tiantian Sun
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Cong Xiao
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Zhuliang Yang
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Jixian Deng
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Xiurong Yang
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

Xiurong Yang and Tiantian Sun designed the study. Tiantian Sun, Cong Xiao, Zhuliang Yang, Jixian Deng, Xiurong Yang conducted the sample preparation. Cong Xiao did the bioinformatic analysis. Tiantian Sun performed the experiment and wrote the main manuscript. Tiantian Sun and Cong Xiao prepared figures. All authors reviewed the manuscript.

Corresponding author

Correspondence to Xiurong Yang.

Ethics declarations

Ethics approval and consent to participate

All experimental and sample collection procedures were approved by the Institutional Animal Care and Use Committee (IACUC) of the College of Animal Science and Technology of Guangxi University (Guangxi, China), with approval number GXU2018-058. All methods were carried out under relevant guidelines and regulations. All methods are reported in accordance with ARRIVE guidelines for the reporting of animal experiments.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Supplementary Information

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

Sun, T., Xiao, C., Yang, Z. et al. Transcriptome profiling analysis of uterus during chicken laying periods. BMC Genomics 24, 433 (2023). https://doi.org/10.1186/s12864-023-09521-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-023-09521-z

Keywords

BMC Genomics

ISSN: 1471-2164

Contact us