Uterine Transcriptome Analysis Reveals Mechanism Underlying the Ultrastructure Differences of Eggshell in Young and Aged Laying Hens

doi:10.21203/rs.3.rs-26128/v3

Download PDF

Research article

Uterine Transcriptome Analysis Reveals Mechanism Underlying the Ultrastructure Differences of Eggshell in Young and Aged Laying Hens

https://doi.org/10.21203/rs.3.rs-26128/v3

This work is licensed under a CC BY 4.0 License

Journal Publication

published 09 Nov, 2020

Read the published version in BMC Genomics →

You are reading this older preprint version

Read the latest preprint version →

Background

Lower eggshell quality in the late laying period leads to economic loss. It is a major threat to the quality and safety of egg products. Age-related variations in ultrastructure were thought to induce this deterioration. Eggshell mineralization is a highly complex process under precise regulation of genes and biological pathways in uterus of laying hens. Herein, we evaluated the physical, mechanical and ultrastructure properties of eggshell and conducted RNA sequencing to learn the transcriptomic differences in uterus between laying hens in the peak (young hens) and late phase (aged hens) of production.

Results

The declined breaking strength and fracture toughness of eggshell were observed in aged hen group compared to those in young hen group, accompanied with ultrastructure variations including the increased thickness of mammillary layer and the decreased incidence of early fusion. During the initial stage of eggshell formation, a total of 183 differentially expressed genes (DEGs; 125 upregulated and 58 downregulated) were identified in uterus of laying hens in the late phase in relative to those at peak production. The DEGs annotated to Gene Ontology terms related to antigen processing and presentation were downregulated in aged hens compared to young hens. The contents of proinflammatory cytokine IL-1β in uterus were higher in aged hens relative to those in young hens. Besides, the genes of some matrix proteins potentially involved in eggshell mineralization, such as ovalbumin, versican and glypican 3, were also differentially expressed between two groups.

Conclusions

Altered gene expression of matrix proteins along with the compromised immune function in uterus of laying hens in the late phase of production may conduce to age-related impairments of eggshell ultrastructure and mechanical properties. The current study enhances our understanding of the age-related deteriorations in eggshell ultrastructure and provides potential targets for improvement of eggshell quality in the late laying period.

Epigenetics & Genomics

Animal Science

Laying hen

Late laying period

Eggshell quality

Ultrastructure

Uterus

Transcriptome

Great economic loss caused by lower eggshell quality is a major concern for egg industry and this problem can be more serious in the late laying period [1-3]. It has been estimated that the incidence of cracked and broken eggs during the late phase of production could reach as high as 12-20% [4], which is one of the key obstacles for extending the laying cycle of commercial flocks [5]. Eggshell is a highly ordered structure comprising membranes (inner and outer), mammillary layer, palisade layer, vertical crystal layer and cuticle. The key roles of ultrastructure have been increasingly recognized in determining shell mechanical properties and quality [6, 7]. Eggshell in the late phase of production possessed lower breaking strength, accompanied with a greater variability in structural properties such as thickness, grain morphology and crystal texture [7]. Besides, the percentages of abnormal structure such as type B mammillae and late fusion in eggshell from aged hens were higher than those from young hens [8], which was thought to contribute to the impairment of mechanical properties. Therefore, age-related abnormalities in ultrastructure could be one reasonable explanation for the compromised shell quality in the late laying period.

The ultrastructure results from the sequential precipitation of mineral carbonate and organic matrix during three stages of mineralization (the initiation, rapid growth and the termination) [9]. Organic matrix plays a crucial role in modulating ultrastructure and mechanical properties of eggshell, as confirmed by in vitro and in vivo observations and by analogies with other biominerals [10, 11]. Eggshell biomineralization is a highly complex process under precise control of genes and biological pathways in uterus [12, 13]. More than 37 ion transporter genes (such as calbindin, ATP2A3, TRPV6 and CA2), which participate in the supply of ions and minerals for biomineralization, have been identified to establish the uterine ion transport model [14]. A recent transcriptome analysis indicated some novel genes (such as ATP2C2, stanniocalcin 2 and calcitonin-related polypeptide B) and three canonical pathways including calcium transport I, cAMP-mediated signaling and cardiac β-adrenergic signaling potentially involved in regulating eggshell formation [15]. Furthermore, the genetic bases of ultrastructure traits and mechanical properties of eggshell have also been extensively explored. For example, genetic variations of ovocleidin-116, ovotransferrin and secreted phosphoprotein 1 were linked to the decreased eggshell quality [16]. Matrix protein genes such as ovalbumin, ovocleidin-116 and ovocalyxin-32 were associated with shell thickness and mammillary layer thickness [17], while ABCC9, ITPR2, KCNJ8 and WNK1, which are involved in ion transport, were suggested to be the key genes modulating eggshell thickness and effective thickness [18]. Hence, deteriorations in eggshell ultrastructure and quality in the late phase of production were speculated to be associated with age-related dysregulation of gene expression in uterus. It was supported by the findings that the expression of genes encoding ion transporters and matrix proteins (ATP2A2, SCNN1G, CA2 and ovocalyxin-36) changed with the age of laying hens [19]. Age-related changes in the expression of matrix proteins (ovalbumin, ovotransferrin and ovocleidin-17) in uterus of laying hens were suggested by their varying concentration in eggshell [20, 21]. Besides, a destructed uterine structure, the reduced gland density [22] and the fibrosis and atrophy of the endometrium [8] of aged hens probably caused by sustained laying behavior would contribute to the uterine hypofunction of protein synthesis, ion transport and immune defense. However, a comprehensive understanding of uterine functional differences in terms of eggshell biomineralization between young and aged hens at the transcriptome level remains obscure.

This study was performed to evaluate age-related changes in eggshell quality and ultrastructure and thus determine the crucial ultrastructure layer as well as its corresponding formation stage responsible for the declined eggshell quality in the late phase of production. Then, RNA sequence was employed to investigate transcriptomic differences in uterus between laying hens in different laying periods at a specific stage of eggshell formation. Our study may contribute to the identification of important genes associated with age-related deterioration in eggshell quality and providing potential targets to improve eggshell quality in the late phase of production.

The physical and mechanical properties of eggshell

The contents of major inorganic materials (calcium and phosphorus) and total matrix protein did not change significantly (P > 0.05) with hen age (Table 1). Egg and shell weight were higher (P < 0.05) in the aged hens relative to the young hens, whereas shell thickness and ratio were similar (P > 0.05) for these two groups. In terms of shell mechanical properties, breaking strength and fracture toughness of eggshell from aged hens were significantly lower (P < 0.05) than those of young ones.

Age-related deterioration in eggshell ultrastructure

There was a significant increase in the thickness of mammillary layer (72.87±6.47 vs. 81.98±9.12, young vs. aged hens, P = 0.037), but a significant decline (P < 0.05) in effective thickness (316.38±6.20 vs. 300.26±17.64, P = 0.029) of aged hens relative to those of young ones (Fig. 1). No differences (P > 0.05) were observed between groups in mammillary knob width, total shell thickness, the ratio of mammillary and effective layer. Most ultrastructure variants in mammillary layer were not affected by hen age (P > 0.05), while there was a reduction in the incidence of early fusion (P < 0.05) in mammillary layer of eggshell in aged hens as compared to young ones (Table 2 & Fig. 2). As a result, total score for mammillary variants was higher (P < 0.05) in aged hen group relative to young one.

Proinflammatory cytokine contents in uterus tissue

The content of proinflammatory cytokine IL-1β in uterus in the aged hens was higher (P < 0.05, Table 3) than that in the young hens, while IL-6 and IFN-γ contents were not affected (P > 0.05) by hen age.

Identification of differentially expressed genes (DEGs) in uterus between two groups

A total of 739,940,606 (220.78 Gb of clean data) clean reads were generated from the thirty-one libraries divided into two groups, with more than 20 million clean reads from each sample (Additional file 1). About 37-52 ´ 10⁶ reads (88-92% of the total raw reads) were uniquely mapped to Gallus gallus genome. More than 91% bases had a quality score of ≥ Q30 and the GC contents of the libraries were ranged from 47.95 to 51.15%, close to 50%, indicating a reliable quality of RNA sequence results. A total of 183 DEGs were identified in uterus between the young and aged hen groups (Fig. 3A & Additional file 2). Among these DEGs, there were 125 significantly up-regulated and 58 significantly down-regulated genes in the aged hens relative to the young ones and the differences in the gene expression profile between groups were visualized in volcano plot (Fig. 3B). We selected 8 genes for quantitative real-time (qRT) -PCR to verify the transcriptome data, and the relative RNA expression was consistent with the results from transcriptome analysis (Fig. 4).

Functional annotation and pathway enrichment analysis of DEGs between two groups

Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis were performed to obtain valuable information for functional prediction of DEGs. In total, the DEGs between groups were functionally distributed into 13 COG (Clusters of Orthologous Groups) categories (Additional file 3). Among them, the greatest number of DEGs were assigned to the category of general function prediction only (31.91%), followed by the category of amino acid transport and metabolism (12.77%), secondary metabolites biosynthesis, transport and catabolism (10.64%), posttranslational modification, protein turnover, chaperones (8.51%), cytoskeleton (8.51%), and inorganic ion transport and metabolism (6.38%). The DEGs were annotated into three major functional categories (Fig. 5) including biological process, cellular component and molecular function. The most enriched terms in the category of biological process were cellular process, single-organism process and biological regulation. Cell, cell part and organelle were most enriched in the category of cellular component and binding, catalytic activity and signal transducer activity were most enriched in the category of molecular function. The DEGs in the aged hen group relative to the young hen group were only enriched in the pathway of peroxisome proliferator-activated receptor (PPAR) signaling pathway (rich factor (RF) = 8.90, Fig. 6). In the PPAR signaling pathway, stearoyl-CoA desaturase (SCD, fold change (FC) =0.61), acyl-CoA synthetase bubblegum family member 1 (ACSBG1, FC=1.51), apolipoprotein C3 (APOC3, FC=0.64), and fatty acid binding protein 5 (FABP5, FC=0.66) were differentially expressed between two groups.

GO clustering analysis of DEGs between two groups

The DEGs were primarily enriched in the following terms within the category of biological process (Table 4): antigen processing and presentation of peptide antigen via major histocompatibility complex (MHC) class I, positive regulation of T cell mediated cytotoxicity, regulation of catabolic process and antigen processing and presentation. Within the category of cellular component, the clusters were enriched in MHC class I, membrane, integral component of membrane, mitochondrion and cytoplasmic part. In terms of the category of molecular function, the enriched terms were peptide antigen binding, catalytic activity, identical protein binding, transferase activity and antigen binding. Since DEGs were predominantly enriched in the GO terms associated with immune response, further analysis was conducted on the DEGs in relation to immune response and immune defense, as well as the DEGs involved in calcification (Table 5). Ten DEGs encoding matrix proteins or calcium transport proteins were suspected to play a role in eggshell biomineralization. Among the immune-related DEGs, MHC class I antigen BF2 (BF2), MHC class I A molecules 5, 6 (MHCIA5 and MHCIA6) and MHC class I α-chain (YF5) in association with antigen processing and presentation were all downregulated in aged hens relative to young ones. The expression of ABCA13 and LOC107056420 related to defense mechanism were all upregulated in the aged hen group as compared to young hen group.

The breaking strength of eggshell declines during the production period of laying hens [23, 24], which was also observed in this study. A recent study evaluated shell quality of eggs sampled every 5 wks, from 30 to 81 wk of hen age and the results showed that eggshell from the peak phase (30-53 wk of age) and the late phase (58-81 wk of age) of production could be distinguished in terms of egg physical and eggshell quality attributes [23]. The decreased shell quality in the late phase of production could be attributed to the increased egg size but no proportional increases in shell weight. This was consistent with our results that a reduction in the contents of calcium in eggshell concurrent with an increase in the weight of egg and eggshell were observed in aged hen group. Besides, it has been suggested that age-associated deteriorations in ultrastructure would contribute to the compromised mechanical properties of eggshell in the late laying period. The variability of ultrastructure quality attributes such as thickness, grain morphology and crystal texture [7], and the occurrence of abnormal structures including type A and B mammillae and mammillary alignment in eggshell [8] were reported to increase with the advance of hen age. In the current study, eggshell ultrastructure in aged hens was characterized by thicker mammillary layer, thinner effective layer and less frequency of early fusion. The thickness and proportion of effective layer were the major ultrastructural characteristics determining shell strength and a reduction in its thickness could compromise shell strength conducing to a higher incidence of breakage [25]. Stress concentrations would occur more frequently and rapidly in the site where the fusion of adjacent mammillary columns was delayed or where structural abnormalities arose. Less early fusion of mammillae in aged hen group suggested a weak binding connection between mammillary knobs, where a crack would propagate easily through shell and outwards from the force points, and therefore compromised its mechanical properties [25, 26]. These changes might further induce negative effects on the structure and quality of following layers [27]. Therefore, ultrastructural variations in eggshell especially in the mammillary layer, may trigger the impairment of shell mechanical properties in aged hen group. RNA sequencing analysis was employed to reveal differences between young and aged hen groups in transcriptome profile of uterus during the initial phase (8.5 h post-oviposition (PO)), corresponding to the stage of mammillary layer formation.

Structural formation is controlled by precise interactions between mineral and organic precursors, which are secreted by specialized oviduct cell populations under gene regulation [13, 28]. Calcium is the major component of eggshell and many calcium transport genes (such as calbindin, ATP2B1, 2 and TRPV6) have been identified to be involved in the supply of calcium for mineralization. In this study, there were three calcium transport-related DEGs including thyrotropin-releasing hormone receptor, calcium voltage-gated channel subunit alpha1 C and 1 D (CACNA1C and CACNA1D) in uterus between young and aged hen groups. CACNA1C and CACNA1D mediate the entry of calcium ions into excitable cells and are also involved in a variety of calcium-dependent processes and MAPK signaling pathways [29]. Besides, solute carriers were proposed to actively involve in ions and minerals transport for eggshell biomineralization in the uterus [15, 30]. The overexpression of SLC9A9 encoding protein would cause excess leak of protons from the recycling endosomes leading to more alkaline endosomes [31]. In our study, the relative higher expression levels of SLC9A9 in the aged hens as compared to the young hens may indicate an age-related change of proton exclusion and Na⁺influx in the uterus. However, their functions in biomineralization have never been reported and need further investigation.

Calcium-binding property is the prerequisite for matrix proteins to participate in biomineralization. Numerus matrix proteins have been identified to play a regulatory role in the initial stage of biomineralization [10, 14, 32]. Ovalbumin, one of the most dominant matrix proteins in this period, could serve as an effective stabilization agent for transient precursors and prevent undirected mineralization [33]. It would preferentially interact with the obtuse rather than the acute steps on the (104) calcite face and thus crystals from the nucleation site could only grow in specific directions. Besides, this effect is closely dependent on its concentration and it would be striking at a higher concentration [34]. In this study, the upregulated expression of ovalbumin with age, consistent with previous observation in shell extracts [21], suggested that calcite crystals would grow in specific and limited directions rather than in all directions, presumptively leading to the delay of mammillary fusion and a thicker mammillary layer. This was supported by the findings that the variations in ovalbumin gene were associated with mammillary layer thickness and breaking strength of eggshell through candidate gene association analysis [17]. Phospholipase C like 1 (PLCL1) protein can inhibit inositol 1,4,5-trisphosphate (IP3)-mediated calcium signaling regulating mechanical sensing of bone cells, supporting its important role in bone formation [35]. Versican (VCAN) encoding a large chondroitin sulfate proteoglycan, could interact with other matrix/cell surface molecules to facilitate establishment or maintenance of early joint interzone structure [36]. Cadherin-6 (CDH6), a single chain transmembrane glycoprotein, can interact selectively and non-covalently with calcium ions and the expression of CDH6 in mouse stromal cells was reported to be associated with the physiological control of bone remodeling [37]. Therefore, the increased expression of ovalbumin, PLCL1, VCAN and CDH6 with the resultant upregulations of GO terms of calcium ion binding might suggest a disordered regulation of matrix proteins, resulting in ultrastructural defects and poor quality of eggshell in the aged hen group.

Proteoglycans (keratan and dermatan sulfate proteoglycans) have an affinity for calcium and were thought to affect the nucleation and assemble of calcite crystals [38]. Glypicans (GPC), belonging to the heparan sulphate proteoglycans, were suspected to be involved in biomineralization. The expression of GPC4 in the shell gland of hens was only found during shell calcification and regulated by egg mechanical strain [39]. GPC3 was regarded as a critical molecule to regulate the activity of bone morphogenetic proteins and promote bone formation. The deletion of GPC3 gene was reported to cause a disorder of bone mineralization and malformation in ribs and skull sutures [40]. In addition, bony suture fusion was associated with the altered expression of GPC3. Accordingly, the downregulations of GPC3 in aged hen group might induce abnormal fusion of mammillary knobs, presumptively contributing to the compromised structure and quality of eggshell. Dental matrix protein-1 (DMP-1) present in eggshell, belongs to the family of the secretory calcium-binding phosphoproteins, some of which (like ovocleidin-116 and osteopontin) have been demonstrated to regulate eggshell calcification [9, 41]. Although DMP-1 has been extensively studied in association with tissue mineralization [42, 43], its effects on mineralization are still not characterized. Our previous data showed that the elevation of mammillary-knob density was associated with the increased expression of proteoglycan and glycoprotein genes in uterus, which contributed to the improvement of eggshell ultrastructure and mechanical properties in the late laying period [44]. However, we did not observe any modifications in mammillary-knob density and the concentration-dependent effects of glycoproteins on ultrastructure attributes need further investigation. Besides, protein modifications were reported to affect the interactions between calcium and matrix proteins, and affect the occlusions of proteins or organics inside the crystals [45]. The enriched term for protein modification process in the current study indicated that age-related co-translational or post-translational modifications of proteins (e.g., glycosylation and phosphorylation) might cause dysfunction of matrix proteins in shell formation. In summary, a disturbance of the expression of matrix protein could inevitably obstruct eggshell calcification, presumably contributing to the impaired shell ultrastructure and quality in aged hen group.

A decline of immune function is a common feature in aged animals. Immune hypofunction was also reported in reproductive tract of aged hens, evidenced by less frequency of immunocompetent cells (Bu-1⁺, CD3⁺, CD4⁺ and CD8⁺cells), which can respond to bacterial and viral stimuli by releasing classical immune mediators [46-48]. In this study, the expression of BF2, YF5, MHCIA5 and MHCIA6, together with the GO terms in association with antigen processing and presentation, were all downregulated in aged hen group. MHC class I and II molecules present antigen peptides to CD8⁺ and CD4⁺ T cells, respectively, which can define the specificity of adaptive immune responses and determine the susceptibility to infectious pathogens [49]. BF2 is the single gene in chickens for classical peptide antigen presentation and MHC class I α-chain (YF) genes encoded molecules are capable of presenting antigen in a non-classical manner but similar to the classical class I molecules [50]. MHCIA5 and MHCIA6 were reported to be involved in the presentation of peptides derived from the endoplasmic reticulum lumen [51]. Therefore, the decreased expression of BF2, YF5, MHCIA5 and MHCIA6 with relevant downregulations of GO clusters of antigen processing and presentation might suggest a hypofunction of antigen processing and presentation, resulting in a compromised immune function and an elevated susceptibility to viral and bacterial infection in uterus of aged hen group. Analogously, age-related disorder in human immune system was associated with the dysfunction in antigen processing and presentation [52, 53]. These findings may provide clues for understanding age-related immune hypofunction of the reproductive tract in other poultry species and in mammals.

The compromised ability to mediate effective immune response against pathogens would elevate the risk of infection and trigger inflammation [54]. Higher contents of IL-1β in this study and the emigration of leukocytes previously reported [55] indicated an inflammation state in uterine mucosa of aged hens. Age-related inflammation characterized by high levels of proinflammatory cytokines has also been reported in bovine endometrial cells [56] and mammal intestinal mucosa [57]. Infection or inflammation could further stimulate the upregulations of proteins possessing antimicrobial defense [58]. The current data revealed that proteins in relation to immune defense, such as LOC107056420 (intestinal mucin-like protein) and ABCA13, were upregulated in aged hens as compared to young hens. Similarly, it was reported that aging induced increased expression of tight junction proteins (occludin) and antibacterial proteins (Avian β-defensins-1, 2, 10 and 11) in the uterine mucosa of laying hens [59]. In support of this view, higher concentrations of immune-defense proteins ovotransferrin and ovocleidin-17 were observed in eggshell from aged hens [21]. However, proteomic study on eggshell and uterine fluid suggested a negative relationship between eggshell quality and proteins involving in defense against bacterial [60]. The upregulation of these proteins in this study may have an adverse impact on eggshell quality. Furthermore, IL-1β has been demonstrated to disrupt the expression of eggshell formation-related proteins such as calcium-binding protein D28K (CaBP-D28K), presumably due to the roles of IL-1β in promoting protein degradation [61]. However, no difference was observed in the expression of calcium transporter and the calcium contents of eggshell between groups. This may be because aging is characterized by a state of chronic, low-grade, systemic inflammation with a moderate increase in IL-1β, which was not enough to cause perturbations of Ca²⁺ transport. In contrast, an intense increase in IL-1β induced by virus infection [62, 63] or uterus culture with 100 ng/mL recombinant chicken IL-1β [61] resulted in a disturbance of CaBP-D28K expression. Recent findings have indicated the important roles of proteoglycans in both immune response and eggshell formation [38, 64]. Age-related alterations in immune functions inevitably disturb the expression of proteoglycans in uterus and consequently disrupt eggshell biomineralization. Therefore, age-related immune hypofunction might disrupt the expression of matrix proteins, probably contributing to the impaired eggshell quality in aged hen group.

Age-related increases in uterus lipid accumulation were regarded as an obstacle to uterine functions, since removal of lipid accumulation by induced molting could result in oviductal regression [65], presumptively favoring the improvement of ultrastructure and eggshell quality. In this study, DEGs including SCD, ACSBG1, APOC3 and FABP5, were mapped to PPAR signaling pathway, which has been considered as important sites for lipid metabolism [66]. It may suggest a disordered regulation of lipid metabolism in the uterus of aged hens, but the effects of these changes on uterine lipid deposition or eggshell biomineralization need further investigation.

This study suggested that decreased incidence of early fusion and increased thickness of mammillary layer are the crucial variations leading to the compromised mechanical properties of eggshell in the late laying period. A disturbed regulation of matrix protein expression and immune hypofunction in uterus of aged laying hens could contribute to age-induced ultrastructural deterioration of eggshell. These findings provide insights into the mechanism underlying age-related deterioration in eggshell ultrastructure, which may contribute to future studies on improving eggshell quality in the late laying period.

Animals and housing

The animal protocol for our study was approved by the Animal Care and Use Committee of the Feed Research Institute of the Chinese Academy of Agricultural Sciences. A total of 96 healthy young (42 wk of age) and 96 healthy aged (72 wk of age) Hy-Line Brown laying hens were individually divided into 8 replicates of 12 birds each in a randomized block design, and 3 birds were placed in one cage. The layer chicks were commercially obtained from Xiaoming Agriculture and Animal Husbandry Co. Ltd. (Ningxia, China) and were separately reared in two houses with similar configurations. Before this experiment, all the laying hens used in this study were raised in one hen house to acclimate the environment for 4 wks. All birds were fed with the same corn-soybean meal basal diet (Additional file 4) and provided with feed and water ad libitum with exposure to 16 h of light/d and the control temperature. All hens remained in good health during the feeding period. The egg production in young and aged hen groups were 95.05±1.90% and 86.94±0.49%, respectively. There were no culled birds, and medical intervention was not applied to any bird.

Sample collection

A total of 24 egg samples (8 eggs/replicate/d × 3 d) from each replicate were collected on three successive days and weighed. Another 12 egg samples (4 eggs/replicate/d × 3 d) from each replicate were collected only for the measurement of stiffness, elastic modules and fracture toughness. An automatic-monitoring control system (FRI, CAAS, Beijing, China) was used to record the daily oviposition time and the total time an egg spent in the oviduct could be calculated. A total of 16 birds (8 replicates with 2 birds each) from each group were sacrificed at 8.5 h PO corresponding to the initiation stage of eggshell calcification. One sample from the aged hen group was rejected for further analysis, because no egg with incomplete shell was present in hen oviduct. The mucosa of uterine tissues surrounding the eggs were collected and snapped frozen immediately in liquid nitrogen, then stored at -80°C until further analysis.

Eggshell physical and mechanical properties

Eggshell thickness was measured by the Egg Shell Gauge and breaking strength was determined by Egg Force Reader (Israel Orka Food Technology Ltd., Ramat Hashron, Israel). Shell stiffness was measured by quasi-static compression using a Texture Lab Pro (TMS-Pro, Food Technology Ltd., VA, USA) fitted with a 25-N load cell at a compression speed of 1 mm/min. Three consecutive measurements were carried out at different points around the equator of each egg. The elastic modules and fracture toughness of each egg were calculated according to the formulas as previously described [67]. The elastic modulus (E, N/mm²) reflects the contribution of eggshell material to the overall stiffness characteristics: E = C[(S_d·R)/T²], where S_d = stiffness (N/mm), R = radius of curvature (breadth/2 (mm)), T = thickness of eggshell (mm), and C = non-dimensionalised constant. The C is determined by the shape index (SI), radius of curvature (R), and the thickness (T) of the eggshell and was calculated as follows: C = A(0.408+(3.026T)/R), where A = (-0.666+(1.8666×(SI))-(0.907×(SI)²) + (0.153×(SI)³))/0.444. Fracture toughness, or resistance to fracture (KC, N/mm^3/2), reflects the nature and magnitude of inherent defects within a material and was calculated as follows: KC = K_nd (F/T^3/2), where K_nd = 0.777 (2.388 + (2.9934 (6/R))), F = breaking strength (N). After removing egg contents, shell was washed, air-dried at room temperature and weighed. The shell ratio was calculated as shell weight/egg weight ´ 100.

Eggshell ultrastructure

Three pieces of shell sample (~0.5 cm²) from the equatorial section of each egg were assessed for the ultrastructure of the external shell surface and the cross section by scanning electronic microscopy (SEM, FEI Quanta 600, Thermo Fisher Scientific Ltd., Portland, OR, USA). Shell specimens for the ultrastructure of mammillary layer were prepared according to previous methods [68, 69]. Briefly, the shell membrane was carefully removed after soaking in distilled water. Then, each sample was immersed overnight in a solution of 6% sodium hypochlorite, 4.12% sodium chloride, and 0.15% sodium hydroxide, washed with distilled water and air-dried at room temperature. Shell samples were glued on a copper block, coated with gold powder, and analyzed with SEM. Mammillary thickness, the effective thickness (total thickness of palisade, vertical crystal layer and cuticle) and the width of the mammillary knobs were determined and calculated with SEM ruler. Mammillary thickness is the length from the top of the membrane to the bottom of the palisade layer. The width of the mammillary knobs was calculated as follows: width = the length of mammillary knobs/ the number of mammillary knobs. Mammillary knob density was counted and expressed as the number of knobs per unit. Each ultrastructural variant in the mammillary layer was assigned a score depending on its incidence in the eggshell. The total score was calculated by the sum of all variant scores [68, 69].

Calcium and phosphorus contents in diets and eggshell

Diet and eggshell samples were dissolved in a solution of 3 mL HNO₃ and 3 mL H₂O₂ and then digested by the microwave dissolution instrument (MDS-10, Shanghai Xinyi Instrument Technology co., Ltd, Shanghai, China). Calcium and phosphorus contents were measured by flame atomic absorption spectrophotometry (Zeenit700P, Analytik Jena, Germany) and a spectrophotometer (UV-2700, Shimadzu, Japan).

Proinflammatory cytokine contents in uterus tissue

The contents of IL-1β, IL-6 and IFN-α in uterus tissues were measured with ELISA kits (Shanghai Enzyme-linked Biotechnology Co., Ltd., Shanghai, China) as recommended by the manufacturer’s instructions. The results were normalized against total protein content in each sample for comparison.

RNA extraction, library preparation and sequencing

Total RNA in the uterus tissues was extracted using TRIzol reagent according to the manufacturers’ instructions (Tiangen Biotech Co. Ltd., Beijing, China). Before library construction, the RNA integrity was assessed using the Agilent Bioanalyzer 2100 system (Agilent Technologies, CA, USA) with the RNA Nano 6000 Assay Kit. A total of 1 μg RNA per sample was used as input material for the RNA sample preparation. RNA library construction was performed using NEBNext^® Ultra^TM RNA Library Prep Kit for Illumina^® (NEB Inc., Ipswich, MA, USA) following the manufacturer’s instruction. Briefly, poly-T oligo-attached magnetic beads were used for the purification of mRNA from total RNA. Fragmentation was performed using divalent cations in NEBNext First Strand Synthesis Reaction Buffer (5×) under elevated temperature. First strand cDNA was synthesized using random hexamer primer and M-MuLV Reverse Transcriptase and second strand cDNA synthesis was subsequently carried out using DNA Polymerase I and RNase H. The rest of the overhangs were converted into blunt ends via exonuclease/polymerase activities. After adenylation of 3’ ends of DNA fragments, NEBNext Adaptor with hairpin loop structure were ligated to prepare for hybridization. AMPure XP system (Beckman Coulter, Beverly, USA) was used to purify the library fragments to select cDNA fragments of preferentially 240 bp in length. The size-selected, adaptor-ligated cDNA was mixed with 3 μl USER Enzyme (NEB Inc., Ipswich, MA, USA) and then incubated at 37°C for 15 min and 95°C for 5 min before PCR. The PCR was carried out with Phusion High-Fidelity DNA polymerase, Universal PCR primers and Index (X) Primer. The products were purified with AMPure XP system (Beckman Coulter, Beverly, USA) and library quality was assessed with Agilent Bioanalyzer 2100 system (Agilent Technologies, CA, USA). The clustering of the index-coded samples was performed with a cBot Cluster Generation System and TruSeq PE Cluster Kit v4-cBot-HS (Illumia) and the library preparations were sequenced on an Illumina platform and paired-end reads were generated.

Sequence quality control and functional annotation of DEGs

Raw reads with adapter, ploy-N and of low quality were removed to obtain clean reads and Q20, Q30, GC-content and sequence duplication level of the clean data were assessed. These clean reads were then mapped to the chicken reference genome (Gallus gallus 5.0) using Hisat2 tools soft. Gene function was annotated based on the following databases: Nt (NCBI non-redundant nucleotide sequences), COG and GO. Gene expression levels were estimated by fragments per kilobase of transcript per million fragments mapped (FPKM). Differentially expression analysis was performed using the DESeq2 and the resulting P values were adjusted using the Benjamini and Hochberg’s approach for controlling the false discovery rate (FDR). The genes at least 1.5-FC and FDR < 0.05 were defined as DEGs between groups. GO analysis of DEGs was performed using the GOseq R packages based Wallenius non-central hyper-geometric distribution [70] and the KEGG pathway enrichment analysis of DEGs was performed using KOBAS software [71].

qRT-PCR validation of RNA sequencing results

Eight genes were selected for qRT-PCR validation. The RNA samples were reverse transcribed with the FastQuant RT kit (KR106, Tiangen Biotech Co. Ltd., Beijing, China) to prepare cDNA. The mRNA expression of target genes was examined by qRT-PCR using CFX96 touch RT-PCR detection system (Bio-rad laboratories. Inc., CA) with a 20 μL PCR reaction mixture (primer concentration: 0.3 μM) according to instructions of the SuperReal PreMix Plus kit (SYBR Green, FP205, Tiangen Co., Beijing, China). Primers used in this study are shown in Additional file 5. The reaction conditions were as follows: 95ºC for 15 min; 40 cycles of 95 ºC for 10 s, 60ºC for 30 s. Each sample was measured in duplicate. The size of all amplified products was confirmed by electrophoresis on a 1.5% (w/v) agarose gel with gelred (SolarGelRed Nucleic Acid Gel Stain, Beijing Solarbio Science & Technology Co., Ltd., Beijing, China) and visualized in Gel Doc XR+ System (Bio-rad laboratories. Inc., CA, USA). The relative mRNA expression levels were normalized to avian β-actin by the 2^-ΔΔCt method [72].

Statistical analysis

Unpaired t-tests (two tailed) were used to analyze the significant differences between groups using SPSS (version 23.0 for Windows; SPSS Inc., Chicago, IL, USA). Data were presented as mean with standard deviation (SD) and statistical significance was defined as a P value < 0.05.

ACSBG1: acyl-CoA synthetase bubblegum family member 1; APOC3: apolipoprotein C3; BF2: MHC class I antigen BF2; CaBP-D28K: calcium-binding protein D28K; CACNA1C: calcium voltage-gated channel subunit alpha1 C; CACNA1D: calcium voltage-gated channel subunit alpha1 D; CDH6: cadherin-6; COG: Clusters of Orthologous Groups; DEGs: differentially expressed genes; DMP1: dental matrix protein-1; FABP5: fatty acid binding protein 5; FC: fold change; FDR: false discovery rate; FPKM: fragments per kilobase of transcript per million fragments mapped; GO: Gene Ontology; GPC: glypicans; KEGG: Kyoto Encyclopedia of Genes and Genomes; MHC: major histocompatibility complex; MHCIA: MHC class I A molecules; Nt: NCBI non-redundant nucleotide sequences; PLCL1: phospholipase C like 1; PO: post-oviposition; PPAR: peroxisome proliferator-activated receptor; RF: rich factor; SCD: stearoyl-CoA desaturase; SEM: scanning electronic microscopy; VCAN: versican; YF5: MHC class I α-chain.

Ethics approval and consent to participate

The animal protocol was approved by the Animal Care and Use Committee of the Feed Research Institute of Chinese Academy of Agricultural Sciences.

Consent for publication

Not applicable.

Availability of data and materials

The RNA-Seq datasets are available in the Sequence Read Archive (SRA) of the NCBI (SRA accession: PRJNA627801).

Competing interests

The authors declare that they have no competing interests.

Funding

This study was financially supported by the National Natural Science Foundation of China (3187131721), the earmarked fund for Modern Agro-industry Technology Research System (CARS-40-K12), China Agriculture Research System-Beijing Team for Poultry Industry (CARS-PSTP, Beijing, China), and the Agricultural Science and Technology Innovation Program (ASTIP) of the Chinese Academy of Agricultural Sciences. Apart from providing funds, they were not involved in the study design, data collection, analysis, interpretation, or manuscript writing.

Authors’ contributions

JF conducted the animal trial, performed the sample analysis and wrote the manuscript. JW, HJZ and SGW assisted with data analysis. JW and GHQ contributed to the experimental design and revised the manuscript. All authors have read and approved the manuscript.

Acknowledgements

We thank Yaoming Cui, Yanan Zhang and Haifeng Zhang for their support in animal feeding, sample collection and sample analysis, and thank Weiwei Wang for his assistance in the revision of manuscript.

Dhawale A. Abnormal eggs cause subnormal profits. World Poult. 2008;24(6):20-23.
Ren ZZ, Piepenburg AJ, Bütz DE, Claus JR, Cook ME. Vaccine to fibroblast growth factor 23 peptides increases eggshell strength. Poult Sci. 2018;97(3):882-889.
Fathi MM, Galal A, Ali UM, Abou-Emera OK. Physical and mechanical properties of eggshell as affected by chicken breed and flock age. Br Poult Sci. 2019;60(5):506-512.
Travel A, Nys Y, Bain M. Effect of hen age, moult, laying environment and egg storage on egg quality. In: Nys Y, Bain M, van Immerseel F, editors. Improving the safety and quality of eggs and egg products. Cambridge: Woodhead Publishing limited; 2011. P. 300-329.
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(3):330-338.
Ketta M, Tůmova E. Eggshell structure, measurements, and quality-affecting factors in laying hens: a review. Czech J Anim Sci. 2016;61(7):299-309.
Athanasiadou D, Jiang W, Goldbaum D, Saleem A, Basu K, Pacella MS, Böhm CF, Chromik RR, Hincke MT, Rodríguez-Navarro AB, Vali H, Wolf SE, Gray JJ, Bui KH, McKee MD. Nanostructure, osteopontin, and mechanical properties of calcitic avian eggshell. Sci Adv. 2018; 4:eaar3219.
Ji-Ae Park, Sea-Hwan Sohn. The influence of hen aging on eggshell ultrastructure and shell mineral components. Korean J Food Sci Anim Resour. 2018;38(5):1080-1091.
Hincke MT, Nys Y, Gautron J, Mann K, Rodriguez-Navarro AB, Mckee MD. The eggshell: structure, composition and mineralization. Front Biosci. 2012;17(1):1266-1280.
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-154.
Gautron J. Proteomics analysis of avian eggshell matrix proteins: toward new advances on biomineralisation. Proteomics. 2019;19(13):1900120.
Jeong W, Lim W, Kim J, Ahn SE, Lee HC, Jeong J-W, 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.
Khan S, Wu S-B, 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.
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.
Sah N, Kuehu D L, Khadka V S, Deng YP, 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[J]. Sci Rep, 2018;8(1):16853.
Zhang Q, Zhu F, Liu L, Zheng CW, Wang DH, Hou ZC, Ning ZH. Integrating transcriptome and genome re-sequencing data to identify key genes and mutations affecting chicken eggshell qualities. Plos One. 2015;10(5):e0125890.
Dunn IC, Joseph NT, Bain M, Edmond A, Wilson PW, Milona P, Nys Y, Gautron J, Schmutz M, Preisinger R, Waddington D. Polymorphisms in eggshell organic matrix genes are associated with eggshell quality measurements in pedigree Rhode Island Red hens. Anim Genet. 2009;40(1):110-114.
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.
Sea, Hwan, Sohn, Jung, Yeon, Park, Ji, Ae, Park, Eun. 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.
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-6941.
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-588.
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.
Sirri F, Zampiga M, Berardinelli A, Meluzzi A. Variability and interaction of some egg physical and eggshell quality attributes during the entire laying hen cycle. Poult Sci. 2018;97(5):1818-1823.
Kemps B, Govaerts T, De Ketelaere B, Mertens K, Bamelis F, Bain M, Decuypere E, De Baerdemaeker J. The influence of line and laying period on the relationship between different eggshell and membrane strength parameters. Poult Sci. 2006;85(7):1309-1317.
Radwan LM. Eggshell quality: a comparison between Fayoumi, Gimieizah and Brown Hy-Line strains for mechanical properties and ultrastructure of their eggshells. Anim Prod Sci. 2015;56(5):908-912.
Dunn IC, Rodríguez-Navarro AB, Mcdade K, Schmutz M, Preisinger R, Waddington D, Wilson PW, Bain MM. Genetic variation in eggshell crystal size and orientation is large and these traits are correlated with shell thickness and are associated with eggshell matrix protein markers. Anim Genet. 2012;43(4):410-418.
Arias JI, Jure C, Wiff JP, Fernández MS, Fuenzalida V, Arias JL. Effect of sulfate content of biomacromolecules on the crystallization of calcium carbonate. Mrs Proc. 2001;711:HH1.7.1.
Yin L, Yu L, Zhang L, Ran J, Li J, Yang C, Jiang X, Du H, Hu X, Liu Y. Transcriptome analysis reveals differentially expressed genes and pathways for oviduct development and defense in prelaying and laying hens. Am J Reprod Immunol. 2019;82(3):e13159.
Mortensen OV. MKP3 eliminates depolarization-dependent neurotransmitter release through downregulation of L-type calcium channel Cav1.2 expression. Cell Calcium. 2013;53(3):224-230.
Zhang J, Wang Y, Zhang C, Xiong M, Rajput SA, Liu Y, Qi D. The differences of gonadal hormones and uterine transcriptome during shell calcification of hens laying hard or weak-shelled eggs. BMC Genomics. 2019;20(1):707.
Kondapalli KC, Llongueras JP, Capilla-González V, Prasad H, Hack A, Smith C, Guerrero-Cázares H, Quiñones-Hinojosa A, Rao R. A leak pathway for luminal protons in endosomes drives oncogenic signalling in glioblastoma. Nat Commun. 2015;6:6289.
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-193.
Wolf SE, Leiterer J, Pipich V, Barrea R, Emmerling F, Tremel W. Strong stabilization of amorphous calcium carbonate emulsion by ovalbumin: gaining insight into the mechanism of 'polymer-induced liquid precursor' processes. J Am Chem Soc. 2011;133(32):12642-12649.
Wang X, Wu C, Tao K, Zhao K, Wang J, Xu H, Xia D, Shan H, Lu JR. Influence of ovalbumin on CaCO3 precipitation during in vitro biomineralization. J Phys Chem B. 2010;114(16):5301-5308.
Liu YZ, Wilson SG, Wang L, Liu XG, Guo YF, Li J, Yan H, Deloukas P, Soranzo N, Chinnapenhorsley U, Cervino A, Williams FM, Xiong DH, Zhang YP, Jin TB, Levy S, Papasian CJ, Drees BM, Hamilton JJ, Recker RR, Spector TM, Deng HW. Identification of PLCL1 gene for hip bone size variation in females in a genome-wide association study. Plos One. 2008;3(9):e3160.
Nagchowdhuri PS, Andrews KN, Robart S, Capehart AA. Versican knockdown reduces interzone area during early stages of chick synovial joint development. Anat Rec. 2012;295(3):397-409.
Marie PJ. Role of N‐cadherin in bone formation. J Cell Physiol. 2002;190(3):297-305.
Fernandez MS, Moya A, Lopez L, Arias JL. Secretion pattern, ultrastructural localization and function of extracellular matrix molecules involved in eggshell formation. Matrix Biol. 2001;19(8):793-803.
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-R861.
Dwivedi PP, Grose RH, Filmus J, Hii CST, Xian CJ, Anderson PJ, Powell BC. Regulation of bone morphogenetic protein signalling and cranial osteogenesis by Gpc1 and Gpc3. Bone. 2013;55(2):367-376.
Horvat-Gordon M, Yu F, Burns D, Leach RM. Ovocleidin (OC 116) is present in avian skeletal tissues. Poult Sci. 2008;87(8):1618-1623.
Fisher LW, Torchia DA, Fohr B, Young MF, Fedarko NS. Flexible structures of sibling proteins, bone sialoprotein, and osteopontin. Biochem Biophys Res Commun. 2001;280(2):460-465.
Horvat-Gordon M, Yu F, Burns D, Leach RM. Ovocleidin (oc 116) is present in avian skeletal tissues. Poult Sci. 2008;87(8):1618-1623.
Zhang YN, Zhang HJ, Wu SG, Wang J, Qi GH. Dietary manganese supplementation affects mammillary knobs of eggshell ultrastructure in laying hens. Poult Sci. 2018;97(4):1253-1262.
Du J, Xu G, Liu C, Zhang R. The role of phosphorylation and dephosphorylation of shell matrix proteins in shell formation: an in vivo and in vitro study. Crystengcomm. 2018;20(27):3905-3916.
Zheng WM, Yoshimura Y, Tamura T. Effects of age and gonadal steroids on the localization of antigen-presenting cells, and T and B cells in the chicken oviduct. J Reprod Fertil. 1998;114(1):45-54.
Barua A, Yoshimura Y. Immunolocalization of MHC-II+ cells in the ovary of immature, young laying and old laying hens Gallus domesticus. J Reprod Fertil 1999;116(2):385-389.
Barua A, Yoshimura Y. Effects of aging and sex steroids on the localization of T cell subsets in the ovary of chicken, Gallus domesticus. Gen Comp Endocrinol. 1999;114(1):28-35.
Koch M, Camp S, Collen T, Avila D, Salomonsen J, Wallny H-J, Hateren Av, Hunt L, Jacob JP, Johnston F, Marston DA, Shaw I, Dunbar PR, Cerundolo V, Jones EY, Kaufman J. Structures of an MHC Class I molecule from B21 chickens illustrate promiscuous peptide binding. Immunity. 2007;27(6):885-899.
Miller MM, Taylor RL. Brief review of the chicken Major Histocompatibility Complex: the genes, their distribution on chromosome 16, and their contributions to disease resistance. Poult Sci. 2016;95(2):375-392.
Gromme M, Uytdehaag FGCM, Janssen H, Calafat J, Binnendijk RSV, Kenter MJH, Tulp A, Verwoerd D, Neefjes J. Recycling MHC Class I molecules and endosomal peptide loading. Proc Natl Acad Sci U. S. A. 1999;96(18):10326-10331.
Clark HL, Banks R, Jones L, Hornick TR, Higgins PA, Burant CJ, Canaday DH. Characterization of MHC-II antigen presentation by B cells and monocytes from older individuals. Clin Immunol. 2012;144(2):172-177.
Wong C, Goldstein DR. Impact of aging on antigen presentation cell function of dendritic cells. Curr Opin Immunol. 2013;25(4):535-541.
Gavazzi G, Krause KH. Ageing and infection. Lancet Infect Dis. 2002;2(11):659-666.
Elhamouly M, Nii T, Isobe N, Yoshimura Y. Expression of pro- and anti-inflammatory cytokines and chemokines during the ovulatory cycle and effects of aging on their expression in the uterine mucosa of laying hens. Cytokine. 2018;111:303-308.
Tanikawa N, Ohtsu A, Kawahara-Miki R, Kimura K, Matsuyama S, Iwata H, Kuwayama T, Shirasuna K. Age-associated mRNA expression changes in bovine endometrial cells in vitro. Reprod Biol Endocrinol. 2017;15(1):63.
Man AL, Bertelli E, Rentini S, Regoli M, Briars G, Marini M, Watson AJM, Nicoletti C. Age-associated modifications of intestinal permeability and innate immunity in human small intestine. Clin Sci. 2015;129(7):515-527.
Yarbrough VL, Winkle S, Herbst-Kralovetz MM. Antimicrobial peptides in the female reproductive tract: a critical component of the mucosal immune barrier with physiological and clinical implications. Hum Reprod Update. 2015;21(3):353-377.
Elhamouly M, Nii T, Isobe N, Yoshimura Y. Age-related modulation of the isthmic and uterine mucosal innate immune defense system in laying hens. Poult Sci. 2019;98(7):3022-3028.
Sun C, Xu G, Yang N. Differential label-free quantitative proteomic analysis of avian eggshell matrix and uterine fluid proteins associated with eggshell mechanical property. Proteomics. 2013;13(23-24):3523-3536.
Takahiro N, Naoki I, Yukinori Y. Effects of interleukin-1 beta and-6 on the expression of ion transporters involved in eggshell mineralization in cultured hen uterine mucosal tissue. J Poult Sci. 2018;55(2):142-149.
Nii T, Isobe N, Yoshimura Y. Effects of avian infectious bronchitis virus antigen on eggshell formation and immunoreaction in hen oviduct. Theriogenology. 2014;81(8):1129-1138.
Qi X, Tan D, Wu C, Tang C, Li T, Han X, Wang J, Liu C, Li R, Wang J. Deterioration of eggshell quality in laying hens experimentally infected with H9N2 avian influenza virus. Vet Res. 2016;47(1):35.
Taylor KR, Gallo RL. Glycosaminoglycans and their proteoglycans: host-associated molecular patterns for initiation and modulation of inflammation. FASEB J. 2006;20(1):9-22.
Gongruttananun N. Induced molt using cassava meal. 2. Effects on eggshell quality, ultrastructure, and pore density in late-phase laying hens. Poult Sci. 2018;97(3):1050-1058.
Levy E, Beaulieu JF, Delvin E, Seidman E, Yotov W, Basque JR, Ménard D. Human crypt intestinal epithelial cells are capable of lipid production, apolipoprotein synthesis, and lipoprotein assembly. J Lipid Res. 2000;41(1):12–22.
Mabe I, Rapp C, Bain M, Nys Y. Supplementation of a corn-soybean meal diet with manganese, copper, and zinc from organic or inorganic sources improves eggshell quality in aged laying hens. Poult Sci. 2003;82(12):1903-1913.
Bain MM. Eggshell strength: a mechanical/ultrastructural evaluation. Ph.D. Thesis. 1990; University of Glasgow: Glasgow.
Bain MM. Eggshell strength: A relationship between the mechanism of failure and the ultrastructural organisation of the mammillary layer. Br Poult Sci. 1992;33(2):303-319.
Young MD, Wakefield MJ, Smyth GK, Oshlack A. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 2010;11(2):R14.
Mao X, Tao C, Olyarchuk JG, Wei L. Automated genome annotation and pathway identification using the KEGG Orthology (KO) as a controlled vocabulary. Bioinformatics. 2005;21(19):3787-3793.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2^−ΔΔCtmethod. Methods. 2001;25(4):402-408.

Additional file 1: Sequence quality and alignment information of uterus in young (42 wk of age; T1-16) and aged hen groups (72 wk of age; T17-31).

Additional file 2: The information of differentially expressed genes (DEGs) in uterus between young (42 wk of age) and aged hen groups (72 wk of age).

Additional file 3: Cluster of Orthologous Genes (COG) classification of differentially expressed genes of uterus in aged hen group (72 wk of age) relative to young hen group (42 wk of age).

Additional file 4: Ingredient and nutrient levels of the experimental diets (air-dried basis).

Additional file 5: Sequences for real-time PCR primers.

Table 1 Comparison of chemical composition, physical and mechanical properties of eggshell between young (42 wk of age) and aged hen groups (72 wk of age)

Items	Young hen	Aged hen	P-value
Calcium content (mg/g)	370.62±19.72	344.68±33.96	0.083
Phosphorus content (μg/g)	1.11±0.12	1.00±0.17	0.169
Matrix protein content (mg/g)	13.72±1.18	13.32±1.18	0.507
Shell thickness (mm)	0.37±0.02	0.39±0.02	0.234
Egg weight (g)	64.42±0.54	65.83±1.12	0.006
Shell weight (g)	6.37±0.05	6.49±0.10	0.008
Shell ratio (%)	9.89±0.05	9.86±0.04	0.244
Breaking strength (N)	46.32±1.64	36.85±2.74	<0.001
Stiffness (N/mm)	74.88±1.91	75.39±3.16	0.702
Elastic modules (N/mm²)	7,636.34±133.61	7,583.75±261.38	0.620
Fracture toughness (N/mm^3/2)	573.11±40.29	456.35±38.75	<0.001

Data are the mean of 8 replicates with 24 eggs each (another 12 eggs of each replicate for the measurement of stiffness, elastic modules and fracture toughness).

Table 2 Comparison of ultrastructural variations in mammillary layer of eggshell between young (42 wk of age) and aged hen groups (72 wk of age)

Ultrastructural variation	Young hen	Aged hen	P-value
Mammillae density	224.53±13.48	221.92±11.80	0.686
Confluence	2.59±1.04	2.97±0.91	0.456
Type B	1.72±0.76	1.94±0.40	0.483
Type A	1.13±0.13	1.16±0.19	0.705
Aragonite	1.03±0.09	1.16±0.19	0.117
Early fusion	1.72±0.62	2.59±0.92	0.042
Late fusion	3.13±1.13	4.22±1.16	0.077
Cuffing	2.97±1.39	3.28±0.75	0.583
Pitted	1.13±0.35	1.25±0.46	0.554
Caps	2.66±0.87	3.06±0.79	0.343
Total score	19.06±2.98	22.63±2.68	0.025

Data are the mean of 8 replicates with 24 eggs each.

Table 3 Comparison of the contents of proinflammatory cytokines in uterine mucosa between young (42 wk of age) and aged hen groups (72 wk of age)

Items	Young hen	Aged hen	P-value
IL-1β (pg/mg protein)	7.47±0.57	8.17±0.47	0.019
IL-6 (pg/mg protein)	0.33±0.04	0.36±0.05	0.155
IFN-γ (pg/mg protein)	0.78±0.07	0.82±0.09	0.348

Data are the mean of 8 replicates with 2 birds each (one sample from the aged hen group was rejected due to no egg with incomplete shell present in hen oviduct).

Table 4 Gene oncology (GO) clustering analysis of differentially expressed genes (DEGs; |fold change| > 1.5 at a false discovery rate < 0.05)

GO Terms	GO.ID	Number of DEGs	P-value
Biological Process
antigen processing and presentation of peptide antigen via MHC class I	GO:0002474	4	<0.001
positive regulation of T cell mediated cytotoxicity	GO:0001916	4	<0.001
regulation of catabolic process	GO:0009894	2	0.005
antigen processing and presentation	GO:0019882	5	0.006
protein ubiquitination involved in ubiquitin-dependent protein catabolic process	GO:0071947	2	0.011
protein modification by small protein conjugation or removal	GO:0070647	5	0.016
cellular protein modification process	GO:0006464	10	0.043
protein modification process	GO:0036211	10	0.043
macromolecule modification	GO:0043412	10	0.045
Cellular Component
MHC class I protein complex	GO:0042612	4	<0.001
membrane	GO:0016020	41	<0.001
integral component of membrane	GO:0016021	25	<0.001
mitochondrion	GO:0005739	4	0.001
cytoplasmic part	GO:0044444	16	0.004
endomembrane system	GO:0012505	8	0.016
membrane part	GO:0044425	34	0.026
endosome	GO:0005768	3	0.029
intrinsic component of membrane	GO:0031224	27	0.042
extracellular matrix	GO:0031012	3	0.043
Molecular Function
peptide antigen binding	GO:0042605	4	<0.001
catalytic activity	GO:0003824	29	<0.001
identical protein binding	GO:0042802	2	0.004
transferase activity	GO:0016740	11	0.005
antigen binding	GO:0003823	4	0.005
calcium ion binding	GO:0005509	5	0.008
adenyl ribonucleotide binding	GO:0032559	3	0.015
ATP binding	GO:0005524	3	0.015
ion binding	GO:0043167	31	0.018
adenyl nucleotide binding	GO:0030554	3	0.018
oxidoreductase activity	GO:0016491	3	0.020
oxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen	GO:0016705	2	0.026
voltage-gated calcium channel activity	GO:0005245	2	0.047

Table 5 Analysis of differentially expressed genes (|fold change| > 1.5 at a false discovery rate < 0.05) in association with matrix protein, calcium transport, immune response and defense mechanism

Potential function	Gene (Fold change)	Annotation
Matrix protein	VCAN (1.55)	Versican, chondroitin sulfate proteoglycan, calcium ion binding
	CDH6 (1.71)	Cadherin 6, calcium binding, in association with the physiological control of bone remodeling
	LOC107056462 (1.61)	Alpha-actinin-1-like, calcium binding
	PLCL1 (1.61)	Phospholipase C like 1, calcium binding
	OVAL (1.52)	Ovalbumin, regulate shell structure and present in eggshell
	GPC3 (0.64)	Glypican 3, regulate the activity of bone morphogenetic proteins and promote bone formation
	DMP1 (0.66)	Dental matrix protein-1, in association with tissue mineralization
Calcium transport	CACNA1C (1.79)	Voltage-dependent L-type calcium channel subunit alpha-1C, ion transport protein, voltage-gated calcium channel activity
	CACNA1D (1.68)	Voltage-dependent L-type calcium channel subunit alpha-1D, ion transport protein, voltage-gated calcium channel activity
	NMUR2 (1.52)	Thyrotropin-releasing hormone receptor, calcium ion transport
Immune response	BF2 (0.66)	MHC class I antigen BF2, antigen processing and presentation
	MHCIA5 (0.58)	Major histocompatibility complex, class I, A5, antigen processing and presentation
	MHCIA6 (0.62)	Major histocompatibility complex, class I, A6, antigen processing and presentation
	YF5 (0.64)	MHC class I antigen YF5, antigen processing and presentation
	LOC101748499 (1.64)	Class II histocompatibility antigen, B-L beta chain-like, immune response
	IL7 (1.63)	Interleukin 7, a cytokine important for B and T cell development
	C7 (1.61)	Complement C7, as part of the terminal complement pathway of the innate immune system
	KCTD7 (1.71)	Potassium channel tetramerization domain containing 7, negative regulation of inflammation
Defense mechanism	ABCA13 (1.65)	ATP-binding cassette, sub-family A, membrane 13, defense mechanism
	LOC107056420 (1.50)	Intestinal mucin-like protein, defense mechanism

Download PDF

Journal Publication

published 09 Nov, 2020

Read the published version in BMC Genomics →

Editorial decision: Major revision
01 Sep, 2020
Review #2 received at journal
31 Aug, 2020
Review #1 received at journal
27 Aug, 2020
Reviewers invited by journal
10 Aug, 2020
Reviewer #1 agreed at journal
10 Aug, 2020
Reviewer #2 agreed at journal
10 Aug, 2020
Editor assigned by journal
05 Aug, 2020
Submission checks completed at journal
04 Aug, 2020
Editor invited by journal
04 Aug, 2020

You are reading this older preprint version

Read the latest preprint version →

Uterine Transcriptome Analysis Reveals Mechanism Underlying the Ultrastructure Differences of Eggshell in Young and Aged Laying Hens

Status:

Journal Publication

Version 3

Abstract

Figures

Background

Results

Discussion

Conclusions

Methods

Abbreviations

Declarations

References

Additional Files

Tables

Supplementary Files

Status:

Journal Publication

Version 3