- Open Access
DNA methylome and transcriptome identified Key genes and pathways involved in Speckled Eggshell formation in aged laying hens
BMC Genomics volume 24, Article number: 31 (2023)
The quality of poultry eggshells is closely related to the profitability of egg production. Eggshell speckles reflect an important quality trait that influences egg appearance and customer preference. However, the mechanism of speckle formation remains poorly understood. In this study, we systematically compared serum immune and antioxidant indices of hens laying speckled and normal eggs. Transcriptome and methylome analyses were used to elucidate the mechanism of eggshell speckle formation.
The results showed that seven differentially expressed genes (DEGs) were identified between the normal and speckle groups. Gene set enrichment analysis (GSEA) revealed that the expressed genes were mainly enriched in the calcium signaling pathway, focal adhesion, and MAPK signaling pathway. Additionally, 282 differentially methylated genes (DMGs) were detected, of which 15 genes were associated with aging, including ARNTL, CAV1, and GCLC. Pathway analysis showed that the DMGs were associated with T cell-mediated immunity, response to oxidative stress, and cellular response to DNA damage stimulus. Integrative analysis of transcriptome and DNA methylation data identified BFSP2 as the only overlapping gene, which was expressed at low levels and hypomethylated in the speckle group.
Overall, these results indicate that aging- and immune-related genes and pathways play a crucial role in the formation of speckled eggshells, providing useful information for improving eggshell quality.
Poultry eggs are one of the most important protein sources, and their relatively low cost makes them popular among consumers. Over the last four decades, egg production has improved considerably due to the development of specialized egg breeds and genetic selection, with the target of “Feeding laying hens to 100 weeks to produce 500 eggs” [1, 2]. However, achieving this target has been limited by the gradual decline in eggshell quality and physiology associated with hen aging, resulting in increased eggshell weight, lighter eggshell color, and speckled eggshells .
The reddish-brown speckle, an important eggshell quality trait, often appears on the blunt end of the brown eggshell, considerably affecting egg appearance and customer preference. The degree of eggshell speckling is evaluated using a scoring method. Speckles can be scored according to pigment intensity, distribution, and speckle size . Moreover, the heritability of speckled eggshells ranges from 0.15–0.2, indicating genetic determination . Additionally, a previous study indicated that aged hens have a higher rate of speckled eggshells than younger hens, reaching 20% after 60 weeks of age .
The eggshell gland is an egg formation organ that plays a critical role in eggshell structure and color formation. During egg formation, the yolk travels through the infundibulum, magnum, and isthmus and reaches the eggshell gland, which secretes a large amount of calcium, pigment, cuticle, and other substances, forming a complete eggshell structure and an outer cuticle [7,8,9]. Any modification or damage to the eggshell gland affects eggshell structure formation and pigmentation [10,11,12].
DNA methylation is one of the earliest known modification pathways, and involves the transfer of methyl groups to the fifth carbon site of cytosine to form 5-methylcytosine . DNA methylation plays an important role in the aging process of animals, regulating age-related gene expression and the etiologies of neurological, immunological, and metabolic diseases [14,15,16,17]. Additionally, DNA methylation combined with environment factors can cause different phenotypes during aging . Several complex livestock phenotypes have been linked to DNA methylation [19, 20]. Recently, RNA sequencing (RNA-seq) has been useful in revealing the genes and pathways underlying traits at the transcriptional level [21,22,23,24], such as embryonic muscle development, feed efficiency, and litter size.
Our previous study showed that, although the occurrence of speckled eggshells does not affect the performance of laying hens (unpublished data), speckles can affect the appearance of eggs and reduce their economic value considerably. Moreover, the molecular mechanisms of speckled egg formation are poorly understood. Therefore, this study aimed to elucidate the mechanism of speckled eggshell formation using transcriptomic and DNA methylation techniques. Since speckled eggshells are heritable and age-related, we used a combination of transcriptomic and DNA methylation analyses to explore the key genes and pathways involved in speckled eggshell formation. It is anticipated that the findings of this study will improve the understanding of the molecular mechanism of eggshell trait formation, which would be useful for animal breeding.
Serum biochemical parameters
A typical egg and speckled egg are shown in Fig. 1. Serum antioxidant and immune indices were measured to determine the physiological status of the laying hens. Serum biochemical parameters are listed in Table 1. Serum levels of immunoglobulin G (IgG) and immunoglobulin A (IgA) are common indicators of humoral immune function. Birds in the normal group had a higher (p = 0.028) IgA content than those in the speckle group. Superoxide dismutase (SOD), catalase (CAT), total antioxidant capacity (T-AOC), glutathione (GSH), and glutathione peroxidase (GSH-PX) are important antioxidant enzymes in the body. MDA is one of the products formed by the reaction of lipids with oxygen radicals, and its content represents the degree of lipid peroxidation. These indices are important in evaluating the oxidative stress process. However, there were no differences (p > 0.05) in antioxidant parameters between the normal and speckle groups.
Transcriptome profile of the eggshell gland
Six cDNA libraries were constructed from the speckle and normal groups. After quality control, a total of 615,170,158 raw reads and 604,275,600 clean reads were obtained (98.22% of the raw reads). After alignment using HISAT2 software, the mapping rate was 90.75–93.18%, and the unique mapping rate in all samples was greater than 73.47% (Supplementary Table S3). Gene expression levels are illustrated using a cluster heatmap and principal component analysis (PCA). There were no significant differences in the gene expression profiles of samples from the speckle and normal groups, as the samples did not form distinct clusters (Fig. 2A, B). A total of seven differentially expressed genes (DEGs) were identified between the normal and speckle groups (p < 0.05, |log2 Fold Change|> 1), including two upregulated and five downregulated genes (Table 2).
Gene set enrichment analysis (GSEA) showed that four pathways were significantly enriched by the expressed genes (Supplementary Table S4). Negative normalized enrichment scores (NES) indicated lower expression levels for some pathways in the normal group compared with those of the speckle group, with the calcium signaling pathway, neuroactive ligand-receptor interaction, focal adhesion, and MAPK signaling pathway being the pathways with the lowest expression (Fig. 2C).
Three DEGs (BFSP2, IQSEC, TMOD4) identified by RNA-seq were verified using quantitative real-time PCR (qRT-PCR). A similar expression of the three genes was evident using RNA-Seq and qRT-PCR, and the coefficient of determination (R2) reached 0.93 (Fig. 2D), indicating that the RNA-seq data were reliable.
DNA methylation profile of the eggshell gland
A total of 564,415,302 and 581,414,308 clean reads were obtained from the speckle and normal groups, respectively, after quality control (Supplementary Table S5), of which 73–78% were uniquely mapped to the converted chicken reference genome (GRCg6a). The cytosine (C) methylation rate of the six eggshell gland samples was approximately 3.4%, and the cytosine site methylation of CpG ranged from 55.5–63.9% in the two groups. The cytosine site methylation of CHH and CHG (H represents A, C, or T) was detected at a low proportion (0.3–0.4%) (Supplementary Table S5). Pearson correlation analysis of the CpG bases suggested that all samples were highly correlated (r > 0.89) (Fig. 3A). PCA showed that samples from the two groups were not significantly different, as they did not form separate clusters (Fig. 3B). There were no significant differences in the methylation levels of CG, CHG, and CHH between the two groups (Fig. 3C). However, the speckle group showed a higher CG methylation level than the normal group. The repeat and exon regions exhibited the highest CG methylation levels, whereas the 5′ UTR region had the lowest CG methylation levels (Fig. 3D).
A total of 2788 differentially methylated regions (DMRs) were identified between the normal and speckle groups. The DMRs were mainly located in the introns (47.45%), followed by the intergenic region (36.05%), exon (8.29%), promoter (5.95%), 3′-UTR (1.18%), and 5′-UTR (0.97%) regions (Fig. 4). Additionally, 282 differentially methylated genes (DMGs) were identified, including 172 hypermethylated and 74 hypomethylated genes in the promoter region. Moreover, 36 DMGs were found in the gene body, including 30 hypermethylated genes and six hypomethylated genes. We converted DMGs into their human orthologs and obtained 158 gene symbols that were uploaded in Metascape for functional annotation, gene ontology (GO), and pathway analyses. The genes were enriched in 176 GO biological processes, including regulation of cellular response to growth factor stimulus and T cell-mediated immunity, response to oxidative stress, and cellular response to DNA damage stimulus. Pathway analysis showed that the genes were enriched in 19 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, including RNA degradation, inflammatory mediator regulation of TRP channels, and the TNF signaling pathway. Moreover, 25 gene sets were detected, including signaling by Rho GTPases, the RHO GTPase cycle, and the CDC42 GTPase cycle. The top 20 enriched ontology clusters are shown in Fig. 5A. Among the 158 homologous human genes, 15 were related to aging or longevity (Supplementary Table S6).
An integrative analysis of the whole-genome bisulfite sequencing (WGBS) and RNA-seq data was performed to determine the relationship between DNA methylation and gene expression levels (Fig. 5B and C). There was a negative correlation between DNA methylation and gene expression levels upstream of the transcription start site (TSS) and downstream of the transcription termination site (TTS); however, there was no correlation between DNA methylation and gene expression levels in the gene body. The Venn diagram showed that BFSP2 was the only overlapping gene between DMGs and DEGs (Fig. 5D).
Studies over several years have shown that the formation of speckles on eggshells is heritable, and aged laying hens produce a higher rate of speckled eggshells than younger hens . Aging of laying hens is often accompanied by chronic inflammation and oxidative stress [25, 26]. The blood parameters of animals reflect their physiological and nutritional status. SOD, GSH-PX, GSH, CAT, and T-AOC are components of the antioxidant defense system. The results of the present study indicate that there was no significant difference in the antioxidant capacity of aged hens between the speckle and normal groups, which is contrary to the findings of Moreno and Osorno . Moreno reported that birds laying eggs with speckled shells may suffer physiological stress because of the pro-oxidant function of the main speckle component and have a higher tolerance of oxidative stress . The differences in the results could be attributed to the different species used in the studies.
Immunoglobulins are proteins involved in anti-inflammatory processes and play regulatory roles in inflammatory reactions . A previous study showed that Blue Tit birds laying speckled eggs exhibit low total immunoglobulin levels . Similarly, the results of the present study showed that birds in the speckle group had substantially lower serum IgA levels compared with those in the normal group, indicating that hens in the speckle group may have a lower anti-inflammatory ability, which is consistent with the findings of Martínez and Merino . IgA, as one of the important indicators in evaluating the humoral immune function of poultry, has antiviral, antibacterial, and antitoxin functions . Studies have shown that dietary supplementation of yeast β-glucan can increase the serum IgA content. The increased IgA content can increase an individual's ability to maintain immune homeostasis, leading to possible health benefits . Studies have also shown that as laying hens age, their resistance to external pathogens decreases . When pathogens enter the body of aged laying hens, the mucosal defense system of oviduct tissues may be damaged due to the concomitant decrease in immunity . However, a higher level of IgA may increase the body's ability to maintain immune homeostasis, preventing the invasion of pathogens into submucosal tissue of the oviduct, thus maintaining good eggshell quality.
Furthermore, RNA-seq analysis identified seven DEGs between the speckle and normal groups, including IQSEC3, BFSP2, TMOD4, LOC112530987, GABRA2, TRIQK, and a pseudogene. IQSEC3, a member of the brefeldin A-resistant ARF guanine nucleotide exchange factors (GEFs) family [35, 36], promotes the development of inhibitory synapses by binding with gephyrin . IQSEC3 is functionally important for maintaining network activity in vivo. IQSEC3 knockdown in the hippocampal dentate gyrus in rodents decreases the density of GABAergic synapses and increases susceptibility to severe seizures . However, there are no studies on the function of IQSEC3 in poultry.
Gamma-aminobutyric acid (GABA) is the main inhibitory neurotransmitter in the central nervous system of vertebrates , and can change the structure of the receptor and the ion permeability of the corresponding receptor . GABA can act directly on tubular smooth muscle cells through GABAA-R or GABAB-R, which are located on the fallopian tube wall and participate by regulating tubular contractility in rabbits , humans , and rats . Studies have shown that GABRA1 plays an important role in egg production. High expression levels of GABRA1 can inhibit the proliferation of granulosa cells, enhance cell apoptosis, and inhibit the synthesis and secretion of progesterone, resulting in reduced egg production [43,44,45]. The egg-laying rhythm of laying hens is a neuromodulatory process [46, 47]. Based on these results, we speculate that there is a subtle relationship between eggshell speckle formation and the nervous system, although the precise link requires further investigation.
TMOD4, a member of a family of proteins that cap the pointed ends of actin filaments , is expressed in skeletal muscle and the heart [49, 50]. Studies have reported that TMOD4 is present in adult chicken lenses, erythrocytes, and fast twitch skeletal muscle fibers . TMOD4 has been found to participate in myofibril assembly, muscle contraction, and differentiation [52,53,54]. The contraction of oviduct muscle makes the egg rotate in the uterus, allowing pigment to be evenly deposited on the eggshell surface. We speculate that the formation of speckled eggs may be related to oviduct muscle contraction, but further research is needed.
GSEA was used to analyze the biological function of all expressed genes to avoid the loss of some interesting gene sets through the cutoff-free strategy. The expressed genes were mainly enriched in the calcium signaling pathway, focal adhesion, the MAPK signaling pathway, and neuroactive ligand–receptor interaction. It has been reported that in the process of eggshell mineralization, genes related to the calcium signaling pathway are involved in the absorption of calcium and carbonate ions from the blood and are transported to the uterine fluid via the epithelial cells of the oviduct to participate in eggshell mineralization [46, 55]. When the calcium and carbonate ions continue to mineralize in the eggshell, the egg rotates continuously in the uterus, and the eggshell pigment, protoporphyrin-IX, can be uniformly deposited on the eggshell surface . We speculate that the difference in the supply of carbonate and calcium ions between the two groups might have led to the uneven distribution of pigment in the deposition process on the eggshell surface, causing eggshell speckles.
Focal adhesions are macromolecular structures that form mechanical connections between the intracellular actin cytoskeleton and extracellular matrix components . It has been shown that focal adhesions play a vital role in maintaining the morphology and function of the oviduct in Chinese brown frogs , while focal adhesions have also been involved in the mechanism of egg production differences in Jinghai Yellow chickens and Nandan-Yao domestic chickens [59, 60]. In the current study, focal adhesions were involved in the formation of eggshell speckles, but the specific participatory mechanism still requires further investigation.
The MAPK pathway has three prominent members: p38 mitogen activated protein kinase (p38 MAPK), Jun N-terminal kinase (JNK), and extracellular response kinase (ERK), which together regulate cell growth, differentiation, apoptosis, inflammation, and other important physiological responses . The ERK 1/2 MAPK pathway has been demonstrated to play an important role in the growth, development, and differentiation of the oviduct and uterus . Wang et al. used Vanadium to induce the growth of oviduct epithelial cells, and found that MAPK family members were activated, leading to oxidative stress of the oviduct, reduced cell activity, and cell apoptosis . In the current study, we found that the MAPK signal pathway in the speckled group was substantially enriched, and therefore speculate that the hens laying speckled eggs might have experienced some stress, with some impact on the oviduct. Neuroactive ligand–receptor interactions are related to steroid hormone synthesis in the gonads. They play an essential role in the regulation of egg production and ovarian function in poultry [59, 64]. The differentially expressed gene, GABRA2, also belongs to the neuroactive ligand–receptor interaction pathway. Since eggshell speckles appear as a result of egg production, and egg production is a rhythmic process, the substantial enrichment of neuroactive ligand–receptor interactions in the speckled group leads us to speculate that there may be differences in egg production performance between hens laying speckled eggs and hens laying normal eggs, but this requires further exploration.
Environmental factors can affect gene expression via epigenetic modifications. The combination of genetic and epigenetic modifications could be helpful in explaining the formation mechanism of complex traits [19, 65]. The relationship between genome-wide DNA methylation and gene expression has been studied for years [66, 67]. Generally, DNA methylation represses gene expression . The results of the present study are consistent with previous findings in that high gene expression levels are associated with low DNA methylation in the promoter region [67, 68]. However, no obvious trend was observed in the gene body region, which may be because gene expression patterns are also regulated by other factors [69, 70]. Functional enrichment analysis revealed that DMGs were mainly enriched in T cell-mediated immunity, response to oxidative stress, and cellular response to DNA damage stimuli, and most of the pathways are associated with aging [71,72,73]. Other aging-related genes have also been identified, including GCLC , CAV1 , and LYN .
A combination of transcriptome and DNA methylation data showed that BFSP2 was the only overlapping gene and was significantly expressed at both DNA methylation and transcriptomic levels. BFSP2 was hypomethylated and its expression levels were low in the speckle group, suggesting gene expression is regulated by other transcription factors in addition to methylation modification. Moreover, BFSP2 has been identified as a candidate gene in autosomal-dominant congenital cataracts  and progressive cataract disease . Although it has been shown that BFSP2 is related to eye development in chickens , studies on other functions in chickens are limited. Therefore, further studies are necessary to elucidate its role in speckled eggshell formation.
In conclusion, the serum immune indices indicate that the IgA content in the speckled group was substantially lower than that in the normal group, and we speculate that a decreased immune function is closely correlated with eggshell speckle formation. Transcriptome analysis detected seven DEGs between the speckle and normal groups, of which IQSEC3, GABRA2, and BFSP2 were identified as potentially important genes associated with speckled egg formation. DNA methylation analysis identified DMGs associated with T cell-mediated immunity, response to oxidative stress, and cellular response to DNA damage stimulus. Of the 282 DMGs identified, 15 were associated with aging. Integrative analysis of transcriptome and DNA methylation revealed that the only overlapping gene was BFSP2, which has barely been studied in chickens. The data presented here suggest that the immune and aging pathways of laying hens may contribute to speckled eggshell formation, improving our understanding of the mechanism of generation of eggshell speckles.
All chicken were obtained from Rongde Breeding Company Ltd. (Hebei, China). Ethical and animal welfare issues were approved by the Ethics Committee of China Agricultural University (permit number: AW02602202-1–1). All experimental protocols were performed according to the Ministry of Science and Technology (Beijing, China). We declare that this study is reported in accordance with ARRIVE guidelines.
A total of 2000 Island Red hens (60-weeks old) were raised in the same environment in individual cages at Hebei Rongde Poultry Breeding Co. Ltd., Hebei, China. The appearance of the eggs was observed according to the method of Cheng  to identify individuals that consistently laid speckled and normal eggs, and those laying heavily speckled eggs and normal eggs; hens were divided into two groups (8 hens per group): speckle and normal groups. The hens had free access to drinking water and were fed at a fixed time every day in an enclosed chicken house under a standard lighting program of 16:8 h (light:dark). The temperature in the chicken house was maintained at 18 ± 1 °C and the relative humidity at 55 ± 5%. All hens were free of avian leukosis viruses and Salmonella pullorum.
Eggshell gland collection
Eggshell speckles are mainly located in the cuticle and vertical crystal layer (unpublished data). Eggshell mineralization stops and the cuticle is deposited in the final 1.5 h before egg expulsion into the eggshell gland [9, 56]. Considering the gene expression dynamics during different stages of egg production, each chicken was sampled at approximately 1.5 h before laying by observing and recording the laying time of the two groups for 7 d. On d 8, eggshell glandular mucosa was collected during the eggshell pigment formation period, placed in cryopreservation tubes, and quickly immersed in liquid nitrogen for preservation.
Serum biochemical parameters
Blood samples were collected from the wing vein of laying hens in the normal (n = 8) and speckle (n = 8) groups into coagulant collection tubes, placed at room temperature for 3 h, and then centrifuged for 10 min at a rate of 3,500 rpm. Serum and plasma were collected for the analysis of blood immune and antioxidant indices. IgG and IgA were used to evaluate the immune performance of an individual hen. SOD, GSH-PX, GSH, CAT, and T-AOC are antioxidant enzymes, while MDA is the product of oxidative stress and is usually used together with antioxidant enzymes as a key indicator to measure antioxidant capacity. Lysozyme (LZM) is used to measure an individual's ability to clear pathogens. IgG, IgA, SOD, MDA, GSH-PX, GSH, CAT, T-AOC, and LZM levels were measured using the colorimetric method with a spectrophotometer (UV-1750, Shimadsu, Japan). The commercial kits of IgG (E026-1–1), IgA (E027-1–1), SOD (A001-1–2), MDA (A003-1–2), GSH-PX (A005-1–2,), GSH (A006-1–1), CAT (A007-1–1), T-AOC (A015-1–2), and LZM (A050-1–1) were purchased from Nanjing Jiancheng Bioengineering Institute (Nanjing, China).
RNA and DNA extraction, library preparation, and sequencing
DNA samples (n = 3 per group, Supplementary Table S1) and total RNA samples (n = 6 per group, Supplementary Table S1) were extracted from the eggshell gland tissue using the animal genomic DNA kit (DP304, Tiangen, Beijing, China) and Trizol reagent (Cat. No. 155–96-026, Invitrogen Life Technologies, Carlsbad, USA), according to the manufacturer’s instructions. DNA and RNA quality and concentration were measured using a spectrophotometer (NanoDrop Technologies, Rockland, DE, USA). Qualified DNA and RNA samples were sent to Personalbio (Shanghai, China) and Majorbio (Shanghai, China), respectively, for library construction. Transcriptome and whole-genome bisulfite sequencing were performed using the Illumina HiSeq X Ten (Illumina, San Diego, CA, United States).
RNA-seq data analysis
Raw reads were quality assessed and filtered using fastp software (0.20.1) . Clean reads were aligned to the chicken reference genome (GRCg6a) using HISAT2 (2.2.1) . The sequencing alignment/mapping (SAM) file was converted to a binary alignment/mapping (BAM) file and sorted using SAMtools (1.11) . The number of reads mapped to each gene was calculated using featureCounts . DEGs were identified using DEseq2 based on the following criteria: P-adjust < 0.05 and |log2 Fold Change|> 1 . GSEA was used to identify significantly different genes in the two groups . Compared with the functional enrichment analysis of DEGs, GSEA uses a cut-free strategy and accumulates subtle expression changes in each gene to compare biological differences. All expressed genes were ranked according to the fold-change value between the two groups. The enrichment score and NES were calculated for each gene set, and significantly different gene sets were identified based on |NES|> 1 and a p-value < 0.05.
Quantitative real-time PCR
Total RNA was extracted from 12 eggshell gland tissues (six hens laying normal eggs and six hens laying speckled eggs) using an RNA extraction kit (Tiangen, Beijing, China) following the manufacturer’s specifications. The mRNA was reversely transcribed into cDNA using a PrimeScript™ RT reagent kit with gDNA Eraser (Takara, Dalian, China). Three DEGs (Due to insufficient samples, we only verified three DEGs, BFSP2, IQSEC, TMOD4) and GAPDH were selected to design primers based on the chicken coding sequence from the NCBI database (Supplementary Table S2). GAPDH was selected as the internal control. The qRT-PCR was conducted using the TB Green® Premix Ex Taq™ Kit (Takara, Kusatsu, Japan) with an ABI 7500 system (Applied Biosystems), using the following program: 95 °C for 30 s, 40 cycles of 95 °C for 5 s, 55 °C for 30 s, 72 °C for 30 s. The 2−ΔΔCt method was used to calculate the relative expression level. The R2 between qRT-PCR and RNA-seq was calculated using simple linear regression in R Studio.
WGBS data analysis
Sequence quality, adapter filtering, and low-quality reads were screened using Fastp (0.20.1) . The clean data were fully bisulfite-converted, aligned against the bisulfite-converted reference sequence, and de-duplicated using Bismark software (0.22.3) . A bismark methylation extractor was used to determine methylation status. The methylation levels of CpG, CHH, and CHG were calculated using the sliding-window approach . The MethylKit package was employed to identify DMRs and DMGs . DMRs were defined when the differential methylation level was > 20% at p ≤ 0.05. Genes that overlapped with the DMRs were defined as DMGs. The region from the TSS to the TTS was defined as the gene body region. The region 2 kb upstream of the TSS and 2 kb downstream of TTS was defined as the promoter region. Pathway enrichment analysis was conducted using Metascape (https://metascape.org/) with default parameters to annotate the functions of DMGs .
Integrative analysis of transcriptome and WGBS Data
To obtain the expression profile of DNA methylation and transcriptome levels, all genes were divided into four groups according to their expression levels, and the average methylation level was calculated based on each group. Overlapping genes between the DMGs and DEGs were analyzed using a Venn plot.
EXCEL and SPSS 25.0 (SPSS Inc., Chicago, IL, USA) were used for statistical analysis. Student’s t-test was used to compare the two groups at a significance level of 5% (p < 0.05). The results are presented as the mean ± standard error.
Availability of data and materials
The RNA sequencing raw data are available at: https://dataview.ncbi.nlm.nih.gov/object/PRJNA850950?reviewer=u9h3k9jracbqg9fkjeabjl2sm2.
The WGBS raw data are available at: https://dataview.ncbi.nlm.nih.gov/object/PRJNA851109?reviewer=k9miop7kfajh6l2mcoa6jha6pq.
Differentially expressed gene
Differentially methylated gene
Whole-genome bisulfite sequencing
Differentially methylated region
Transcriptional start site
Transcriptional terminal site
Principal component analysis
Kyoto Encyclopedia of Genes and Genomes
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–8.
Pottgüter R. Feeding laying hens to 100 weeks of age. Lohmann Inf. 2016;50:18–21.
Molnar 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(6):842–7.
Gosler AG, Higham JP, James Reynolds S. Why are birds’ eggs speckled? Ecol Lett. 2005;8(10):1105–13.
Arango J, Settar P, Arthur J, O’Sullivan N. Relationship between shell color and incidence of speckles in brown egg lines. In: Proc XIIth European Poultry Conference: 2006. 2006. p. 10–4.
Cheng X, Fan C, Ning Z. Quality of freckle eggs and its influencing factors. Chinese Poult Sci. 2019;41(19):6–9 (In Chinese).
Hincke MT, Nys Y, Gautron J, Mann K, Rodriguez-Navarro AB, McKee MD. The eggshell: structure, composition and mineralization. Front Biosci (Landmark Ed). 2012;17(4):1266–80.
Samiullah S, Roberts JR, Chousalkar K. Eggshell color in brown-egg laying hens - a review. Poult Sci. 2015;94(10):2566–75.
Wilson PW, Suther CS, Bain MM, Icken W, Jones A, Quinlan-Pluck F, Olori V, Gautron J, Dunn IC. Understanding avian egg cuticle formation in the oviduct: a study of its origin and deposition. Biol Reprod. 2017;97(1):39–49.
Zhu M, Li H, Miao L, Li L, Dong X, Zou X. Dietary cadmium chloride impairs shell biomineralization by disrupting the metabolism of the eggshell gland in laying hens. J Anim Sci. 2020;98(2):skaa025.
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:35.
Wang J, Yuan Z, Zhang K, Ding X, Bai S, Zeng Q, Peng H, Celi P. Epigallocatechin-3-gallate protected vanadium-induced eggshell depigmentation via P38MAPK-Nrf2/HO-1 signaling pathway in laying hens. Poult Sci. 2018;97(9):3109–18.
Avery OT, Macleod CM, McCarty M. Studies on the chemical nature of the substance inducing transformation of pneumococcal types : induction of transformation by a desoxyribonucleic acid fraction isolated from pneumococcus type Iii. J Exp Med. 1944;79(2):137–58.
Kochmanski J, Marchlewicz EH, Cavalcante RG, Sartor MA, Dolinoy DC. Age-related epigenome-wide DNA methylation and hydroxymethylation in longitudinal mouse blood. Epigenetics. 2018;13(7):779–92.
Zhang X, Hu M, Lyu X, Li C, Thannickal VJ, Sanders YY. DNA methylation regulated gene expression in organ fibrosis. Biochim Biophys Acta Mol Basis Dis. 2017;1863(9):2389–97.
Ehrlich M. DNA hypermethylation in disease: mechanisms and clinical relevance. Epigenetics. 2019;14(12):1141–63.
Ling C, Ronn T. Epigenetics in Human Obesity and Type 2 Diabetes. Cell Metab. 2019;29(5):1028–44.
Kaminen-Ahola N, Ahola A, Maga M, Mallitt KA, Fahey P, Cox TC, Whitelaw E, Chong S. Maternal ethanol consumption alters the epigenotype and the phenotype of offspring in a mouse model. PLoS Genet. 2010;6(1):e1000811.
Tan X, Liu R, Xing S, Zhang Y, Li Q, Zheng M, Zhao G, Wen J. Genome-Wide Detection of Key Genes and Epigenetic Markers for Chicken Fatty Liver. Int J Mol Sci. 2020;21(5):1800.
Hwang JH, An SM, Kwon S, Park DH, Kim TW, Kang DG, Yu GE, Kim IS, Park HC, Ha J, et al. DNA methylation patterns and gene expression associated with litter size in Berkshire pig placenta. PLoS ONE. 2017;12(9):e0184539.
Yu C, Qiu M, Zhang Z, Song X, Du H, Peng H, Li Q, Yang L, Xiong X, Xia B, et al. Transcriptome sequencing reveals genes involved in cadmium-triggered oxidative stress in the chicken heart. Poult Sci. 2021;100(3):100932.
Boo SY, Tan SW, Alitheen NB, Ho CL, Omar AR, Yeap SK. Transcriptome analysis of chicken intraepithelial lymphocyte natural killer cells infected with very virulent infectious bursal disease virus. Sci Rep. 2020;10(1):18348.
Ren L, Liu A, Wang Q, Wang H, Dong D, Liu L. Transcriptome analysis of embryonic muscle development in Chengkou Mountain Chicken. BMC Genomics. 2021;22(1):431.
Yang C, Han L, Li P, Ding Y, Zhu Y, Huang Z, Dan X, Shi Y, Kang X. Characterization and Duodenal Transcriptome Analysis of Chinese Beef Cattle With Divergent Feed Efficiency Using RNA-Seq. Front Genet. 2021;12:741878.
Wang J, Jia R, Gong H, Celi P, Zhuo Y, Ding X, Bai S, Zeng Q, Yin H, Xu S, et al. The Effect of Oxidative Stress on the Chicken Ovary: Involvement of Microbiota and Melatonin Interventions. Antioxidants (Basel). 2021;10(9):1422.
Attia YA, Al-Harthi MA, Abo El-Maaty HM. Calcium and Cholecalciferol Levels in Late-Phase Laying Hens: Effects on Productive Traits, Egg Quality, Blood Biochemistry, and Immune Responses. Front Vet Sci. 2020;7:389.
Moreno J, Osorno JL. Avian egg colour and sexual selection: does eggshell pigmentation reflect female condition and genetic quality? Ecol Lett. 2003;6(9):803–6.
Afonso S, Vanore G, Batlle A. Protoporphyrin IX and oxidative stress. Free Radic Res. 1999;31(3):161–70.
Schwartz-Albiez R, Monteiro RC, Rodriguez M, Binder CJ, Shoenfeld Y. Natural antibodies, intravenous immunoglobulin and their role in autoimmunity, cancer and inflammation. Clin Exp Immunol. 2009;158(Suppl 1):43–50.
Martínez-de la Puente J, Merino S, Moreno J, Tomás G, Morales J, Lobato E, García-Fraile S, Martínez J. Are eggshell spottiness and colour indicators of health and condition in blue tits Cyanistes caeruleus? J Avian Biol. 2007;38(3):377–84.
Fagarasan S, Honjo T. Intestinal IgA synthesis: regulation of front-line body defences. Nat Rev Immunol. 2003;3(1):63–72.
Zhen W, Shao Y, Wu Y, Li L, Pham VH, Abbas W, Wan Z, Guo Y, Wang Z. Dietary yeast beta-glucan supplementation improves eggshell color and fertile eggs hatchability as well as enhances immune functions in breeder laying hens. Int J Biol Macromol. 2020;159:607–21.
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–8.
Li J, Qin Q, Li YX, Leng XF, Wu YJ. Tri-ortho-cresyl phosphate exposure leads to low egg production and poor eggshell quality via disrupting follicular development and shell gland function in laying hens. Ecotoxicol Environ Saf. 2021;225:112771.
Fukaya M, Kamata A, Hara Y, Tamaki H, Katsumata O, Ito N, Takeda S, Hata Y, Suzuki T, Watanabe M, et al. SynArfGEF is a guanine nucleotide exchange factor for Arf6 and localizes preferentially at post-synaptic specializations of inhibitory synapses. J Neurochem. 2011;116(6):1122–37.
Um JW. Synaptic functions of the IQSEC family of ADP-ribosylation factor guanine nucleotide exchange factors. Neurosci Res. 2017;116:54–9.
Um JW, Choii G, Park D, Kim D, Jeon S, Kang H, Mori T, Papadopoulos T, Yoo T, Lee Y, et al. IQ Motif and SEC7 Domain-containing Protein 3 (IQSEC3) Interacts with Gephyrin to Promote Inhibitory Synapse Formation. J Biol Chem. 2016;291(19):10119–30.
Kim S, Kim H, Park D, Kim J, Hong J, Kim JS, Jung H, Kim D, Cheong E, Ko J, et al. Loss of IQSEC3 Disrupts GABAergic Synapse Maintenance and Decreases Somatostatin Expression in the Hippocampus. Cell Rep. 2020;30(6):1995-2005 e1995.
Erdö SL, Riesz M, Kárpáti E, Szporny L. GABAB receptor-mediated stimulation of the contractility of isolated rabbit oviduct. Eur J Pharmacol. 1984;99(4):333–6.
Roth FC, Draguhn A. GABA metabolism and transport: effects on synaptic efficacy. Neural Plast. 2012;2012:805830.
Erdő SL, László Á, Szporny L, Zsolnai B. High density of specific GABA binding sites in the human fallopian tube. Neurosci Lett. 1983;42(2):155–60.
Erdo SL, Rosdy B, Szporny L. Higher GABA concentrations in fallopian tube than in brain of the rat. J Neurochem. 1982;38(4):1174–6.
Sun X, Chen X, Zhao J, Ma C, Yan C, Liswaniso S, Xu R, Qin N. Transcriptome comparative analysis of ovarian follicles reveals the key genes and signaling pathways implicated in hen egg production. BMC Genomics. 2021;22(1):899.
Chen X, Sun X, Chimbaka IM, Qin N, Xu X, Liswaniso S, Xu R, Gonzalez JM. Transcriptome analysis of ovarian follicles reveals potential pivotal genes associated with increased and decreased rates of chicken egg production. Front Genet. 2021;12:622751.
Luan X, Liu D, Cao Z, Luo L, Liu M, Gao M, Zhang X. Transcriptome profiling identifies differentially expressed genes in Huoyan goose ovaries between the laying period and ceased period. PLoS ONE. 2014;9(11):e113211.
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.
Mishra SK, Chen B, Zhu Q, Xu Z, Ning C, Yin H, Wang Y, Zhao X, Fan X, Yang M, et al. Transcriptome analysis reveals differentially expressed genes associated with high rates of egg production in chicken hypothalamic-pituitary-ovarian axis. Sci Rep. 2020;10(1):5976.
Cox PR, Siddique T, Zoghbi HY. Genomic organization of Tropomodulins 2 and 4 and unusual intergenic and intraexonic splicing of YL-1 and Tropomodulin 4. BMC Genomics. 2001;2:7.
Cox PR, Zoghbi HY. Sequencing, expression analysis, and mapping of three unique human tropomodulin genes and their mouse orthologs. Genomics. 2000;63(1):97–107.
Yamashiro S, Gokhin DS, Kimura S, Nowak RB, Fowler VM. Tropomodulins: pointed-end capping proteins that regulate actin filament architecture in diverse cell types. Cytoskeleton (Hoboken). 2012;69(6):337–70.
Almenar-Queralt A, Lee A, Conley CA, Ribas de Pouplana L, Fowler VM. Identification of a novel tropomodulin isoform, skeletal tropomodulin, that caps actin filament pointed ends in fast skeletal muscle. J Biol Chem. 1999;274(40):28466–75.
Ren T, Li Z, Zhou Y, Liu X, Han R, Wang Y, Yan F, Sun G, Li H, Kang X. Sequencing and characterization of lncRNAs in the breast muscle of Gushi and Arbor Acres chickens. Genome. 2018;61(5):337–47.
Xu L, Zhao F, Ren H, Li L, Lu J, Liu J, Zhang S, Liu GE, Song J, Zhang L, et al. Co-expression analysis of fetal weight-related genes in ovine skeletal muscle during mid and late fetal development stages. Int J Biol Sci. 2014;10(9):1039–50.
Wu Y, Wang Y, Yin D, Mahmood T, Yuan J. Transcriptome analysis reveals a molecular understanding of nicotinamide and butyrate sodium on meat quality of broilers under high stocking density. BMC Genomics. 2020;21(1):412.
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.
Nys Y, Gautron J, Garcia-Ruiz JM, Hincke MT. Avian eggshell mineralization: biochemical and functional characterization of matrix proteins. CR Palevol. 2004;3(6–7):549–62.
Chen CS, Alonso JL, Ostuni E, Whitesides GM, Ingber DE. Cell shape provides global control of focal adhesion assembly. Biochem Biophys Res Commun. 2003;307(2):355–61.
Su H, Zhang H, Wei X, Pan D, Jing L, Zhao D, Zhao Y, Qi B. Comparative Proteomic Analysis of Rana chensinensis Oviduct. Molecules. 2018;23(6):1384.
Zhang T, Chen L, Han K, Zhang X, Zhang G, Dai G, Wang J, Xie K. Transcriptome analysis of ovary in relatively greater and lesser egg producing Jinghai Yellow Chicken. Anim Reprod Sci. 2019;208:106114.
Sun T, Xiao C, Deng J, Yang Z, Zou L, Du W, Li S, Huo X, Zeng L, Yang X. Transcriptome analysis reveals key genes and pathways associated with egg production in Nandan-Yao domestic chicken. Comp Biochem Physiol Part D Genomics Proteomics. 2021;40:100889.
Peter AT, Dhanasekaran N. Apoptosis of granulosa cells: a review on the role of MAPK-signalling modules. Reprod Domest Anim. 2003;38(3):209–13.
Jeong W, Kim J, Ahn SE, Lee SI, Bazer FW, Han JY, Song G. AHCYL1 is mediated by estrogen-induced ERK1/2 MAPK cell signaling and microRNA regulation to effect functional aspects of the avian oviduct. PLoS ONE. 2012;7(11):e49204.
Wang J, Huang X, Zhang K, Mao X, Ding X, Zeng Q, Bai S, Xuan Y, Peng H. Vanadate oxidative and apoptotic effects are mediated by the MAPK-Nrf2 pathway in layer oviduct magnum epithelial cells. Metallomics. 2017;9(11):1562–75.
Tao Z, Song W, Zhu C, Xu W, Liu H, Zhang S, Huifang L. Comparative transcriptomic analysis of high and low egg-producing duck ovaries. Poult Sci. 2017;96(12):4378–88.
Shioda K, Odajima J, Kobayashi M, Kobayashi M, Cordazzo B, Isselbacher KJ, Shioda T. Transcriptomic and Epigenetic Preservation of Genetic Sex Identity in Estrogen-feminized Male Chicken Embryonic Gonads. Endocrinology. 2021;162(1):bqaa208.
Huang YZ, Sun JJ, Zhang LZ, Li CJ, Womack JE, Li ZJ, Lan XY, Lei CZ, Zhang CL, Zhao X, et al. Genome-wide DNA methylation profiles and their relationships with mRNA and the microRNA transcriptome in bovine muscle tissue (Bos taurine). Sci Rep. 2014;4:6546.
Tan X, Liu R, Zhang Y, Wang X, Wang J, Wang H, Zhao G, Zheng M, Wen J. Integrated analysis of the methylome and transcriptome of chickens with fatty liver hemorrhagic syndrome. BMC Genomics. 2021;22(1):8.
Fu Y, Li J, Tang Q, Zou C, Shen L, Jin L, Li C, Fang C, Liu R, Li M, et al. Integrated analysis of methylome, transcriptome and miRNAome of three pig breeds. Epigenomics. 2018;10(5):597–612.
Lawrence M, Daujat S, Schneider R. Lateral thinking: how histone modifications regulate gene expression. Trends Genet. 2016;32(1):42–56.
Zhao BS, Roundtree IA, He C. Post-transcriptional gene regulation by mRNA modifications. Nat Rev Mol Cell Biol. 2017;18(1):31–42.
Chakravarti B, Abraham GN. Aging and T-cell-mediated immunity. Mech Ageing Dev. 1999;108(3):183–206.
Camougrand N, Rigoulet M. Aging and oxidative stress: studies of some genes involved both in aging and in response to oxidative stress. Respir Physiol. 2001;128(3):393–401.
Ribezzo F, Shiloh Y, Schumacher B. Systemic DNA damage responses in aging and diseases. Semin Cancer Biol. 2016;37–38:26–35.
Orr WC, Radyuk SN, Prabhudesai L, Toroser D, Benes JJ, Luchak JM, Mockett RJ, Rebrin I, Hubbard JG, Sohal RS. Overexpression of glutamate-cysteine ligase extends life span in Drosophila melanogaster. J Biol Chem. 2005;280(45):37331–8.
Park DS, Cohen AW, Frank PG, Razani B, Lee H, Williams TM, Chandra M, Shirani J, De Souza AP, Tang B, et al. Caveolin-1 null (-/-) mice show dramatic reductions in life span. Biochemistry. 2003;42(51):15124–31.
Park JW, Ji YI, Choi YH, Kang MY, Jung E, Cho SY, Cho HY, Kang BK, Joung YS, Kim DH, et al. Candidate gene polymorphisms for diabetes mellitus, cardiovascular disease and cancer are associated with longevity in Koreans. Exp Mol Med. 2009;41(11):772–81.
Jakobs PM, Hess JF, FitzGerald PG, Kramer P, Weleber RG, Litt M. Autosomal-dominant congenital cataract associated with a deletion mutation in the human beaded filament protein gene BFSP2. Am J Hum Genet. 2000;66(4):1432–6.
Conley YP, Erturk D, Keverline A, Mah TS, Keravala A, Barnes LR, Bruchis A, Hess JF, FitzGerald PG, Weeks DE, et al. A juvenile-onset, progressive cataract locus on chromosome 3q21-q22 is associated with a missense mutation in the beaded filament structural protein-2. Am J Hum Genet. 2000;66(4):1426–31.
Kumar P, Kasiviswanathan D, Sundaresan L, Kathirvel P, Veeriah V, Dutta P, Sankaranarayanan K, Gupta R, Chatterjee S. Harvesting clues from genome wide transcriptome analysis for exploring thalidomide mediated anomalies in eye development of chick embryo: Nitric oxide rectifies the thalidomide mediated anomalies by swinging back the system to normal transcriptome pattern. Biochimie. 2016;121:253–67.
Chen S, Zhou Y, Chen Y, Gu J. fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics. 2018;34(17):i884–90.
Kim D, Paggi JM, Park C, Bennett C, Salzberg SL. Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat Biotechnol. 2019;37(8):907–15.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R, Genome Project Data Processing S. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009;25(16):2078–9.
Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30(7):923–30.
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.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102(43):15545–50.
Krueger F, Andrews SR. Bismark: a flexible aligner and methylation caller for Bisulfite-Seq applications. Bioinformatics. 2011;27(11):1571–2.
Jones PA. Functions of DNA methylation: islands, start sites, gene bodies and beyond. Nat Rev Genet. 2012;13(7):484–92.
Akalin A, Kormaksson M, Li S, Garrett-Bakelman FE, Figueroa ME, Melnick A, Mason CE. methylKit: a comprehensive R package for the analysis of genome-wide DNA methylation profiles. Genome Biol. 2012;13(10):R87.
Zhou Y, Zhou B, Pache L, Chang M, Khodabakhshi AH, Tanaseichuk O, Benner C, Chanda SK. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. 2019;10(1):1523.
We are deeply grateful to all the participants involved in the study and all the research staff and students working on the project.
This research was funded by the China Agriculture Research System (CARS-40) and Beijing Agriculture Innovation Consortium for Poultry (BAIC06-2022-G01).
Ethics approval and consent to participate
All chicken were obtained from Rongde Breeding Company Ltd. (Hebei, China). Ethical and animal welfare issues were approved by the Ethics Committee of China Agricultural University (permit number: AW02602202-1–1). All experimental protocols were performed according to the Ministry of Science and Technology (Beijing, China). We declare that this study is reported in accordance with ARRIVE guidelines.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Supplementary Table 1.
Information on sequencing samples.
Additional file 2: Supplementary Table 2.
Specific primers for qRT-PCR.
Additional file 3: Supplementary Table 3.
Summary statistics for transcriptome sequence quality and alignment information of eggshell gland.
Additional file 4: Supplementary Table 4.
The gene set between between normal and speckle group.
Additional file 5: Supplementary Table 5.
Summary statistics for WGBS alignment and methylation information of eggshell gland.
Additional file 6: Supplementary Table 6.
Genes related to senescence and longevity.
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.
About this article
Cite this article
Cheng, X., Li, X., Liu, Y. et al. DNA methylome and transcriptome identified Key genes and pathways involved in Speckled Eggshell formation in aged laying hens. BMC Genomics 24, 31 (2023). https://0-doi-org.brum.beds.ac.uk/10.1186/s12864-022-09100-8
- Laying hens
- Speckled eggs
- DNA methylation