Next Article in Journal
Autophosphorylation and Self-Activation of DNA-Dependent Protein Kinase
Previous Article in Journal
A Large-Scale Genome-Wide Association Study of Epistasis Effects of Production Traits and Daughter Pregnancy Rate in U.S. Holstein Cattle
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparison of the Transcriptomic and Epigenetic Profiles of Gonadal Primordial Germ Cells of White Leghorn and Green-Legged Partridgelike Chicken Embryos

Department of Animal Biotechnology and Genetics, UTP University of Science and Technology, 85-796 Bydgoszcz, Poland
*
Author to whom correspondence should be addressed.
Submission received: 15 June 2021 / Revised: 11 July 2021 / Accepted: 17 July 2021 / Published: 19 July 2021
(This article belongs to the Section Animal Genetics and Genomics)

Abstract

:
The Green-legged Partridgelike fowl is a native, dual-purpose Polish chicken. The White Leghorn has been intensively selected for several decades to mainly improve reproductive traits. Primordial germ cells (PGCs) represent the germline stem cells in chickens and are the only cells that can transfer the information stored in the genetic material from generation to generation. The aim of the study was to carry out a transcriptomic and an epigenetic comparison of the White Leghorn and Green-legged Partridgelike gonadal PGCs (gPGCs) at three developmental stages: days 4.5, 8, and 12 of the embryonic development. RNA and DNA were isolated from collected gPGCs. The RNA was further subjected to microarray analysis. An epigenetic analysis was performed based on the global methylation analysis and qMSP method for the particular silenced genes demonstrated in transcriptomic analysis. Statistically significant differences between the gPGCs from both breeds were detected on the day 8 of embryonic development. Global methylation analysis showed significant changes at the methylation level in the White Leghorn gPGCs on day 8 of embryonic development. The results suggest faster development of Green-legged Partridgelike embryos as compared to White Leghorn embryos. Changes in the levels of gene expression during embryonic development are determined by genetic and environmental factors, and this variability is influenced by breed and gender.

1. Introduction

Intensive genetic selection of chicken for the economically important production traits, which has been carried out for almost a century, has resulted in the considerable differences in the mechanisms of growth and development and, as a consequence, in the metabolism and reproduction of avians [1]. These differences might be identified at the early stages of avian embryogenesis [2,3,4]. Domestic chicken is also an important model for research in experimental embryology [5] focusing on the molecular basis of development and cell–cell interactions [6]. It is also used for studies in the fields of genomics [7], experimental medicine [8], and epigenetics [9].
Primordial germ cells (PGCs) represent the germline stem cells in chicken, which can be identified at an early stage X of HamburgerHamilton (HH). Further differentiation of these cells depends on the formation of the embryo [10]. Initially, they are attached to the posterior side of the hypoblast, but as growth progresses, they move and finally reach the anterior part of the blastodisc, the so-called germinal crescent. When blood vessels are formed, PGCs, called circulating PGCs (cPGCs), pass into the embryonic bloodstream, through which they enter the genital ridges. After settling into the future gonadal region, gonadal PGCs (gPGCs) proliferate and begin to differentiate into male and female reproductive cells [11]. After a period of mitotic proliferation, gPGCs undergo meiosis and differentiate into functional gametes—either ova or sperms. Post-mating, the fertilization of ova and sperm initiates embryogenesis. Subsequently, new PGCs begin to be produced in the embryo. Thus, PGCs are the only cells that can transfer the information stored in the genetic material from generation to generation. They are potentially of a great importance in understanding the epigenetic regulatory mechanisms in vertebrates. However, due to limited cell numbers, the epigenetic regulation of gene expression during the PGC determination and cell migration in vivo has been rarely studied [12].
In our previous study on gPGCs isolated from an intensively selected breed, White Leghorn (WL), we investigated the dynamics of the transcriptome during chicken embryo development [13]. We concluded that the changes in the gene expression depend on the time at which the PGCs were isolated during embryonic development and the highest number of genes was regulated between days 8 and 12 of an incubation. However, the study lacked information about the effect of sex and breed of the embryo on the molecular level in chicken PGCs. In the current study, we compared the leading breed in egg production—WL chicken—with a common native dual-purpose chicken (Green-legged Partridgelike—GP) in an identical experimental set-up of transcriptomic analysis and performed global and gene-specific methylation analysis in both the breeds of chicken considering the effect of gender. It is an innovative examination of the molecular changes occurring at the level of embryonic development related to gene expression.
This line of research was chosen based on an analysis of the WL transcriptome, which showed that intensive gene expression silencing also occurs. We assume that this process may be related to the DNA methylation. The impact of the DNA methylation on host genome depends on the gene region in which it occurs. DNA methylation of regulatory elements suppresses gene expression. Hence, DNA methylation of gene-deficient regions is considered crucial for the maintenance of chromosome structure and integrity. Methylation of CpG dinucleotides might be inherited through cell division and transmitted to the next generation by germ cells. DNA hypermethylation of these dinucleotides causes silencing of gene expression, whereas DNA hypomethylation of these is related to gene expression [14].
In this study, we hypothesize that changes at the level of gene expression during embryonic development are determined by genetic factors, and the time variation in their occurrence is determined by breed and gender. The study is a follow-up transcriptomic and epigenetic analysis of gPGCs in two chicken lines with different genetic backgrounds.

2. Materials and Methods

Experimental set-up is presented in Figure 1.

2.1. Gonadal PGC Isolation

Fertilized eggs of GP were incubated at 37.8 °C for 4.5, 8, and 12 days to obtain embryos of suitable developmental stage. Gonadal PGCs (n = 6/each group) were collected based on the method described by Nakajima et al. [15]. Gonads were obtained from the embryos on day 4.5, 8, and 12 of incubation (stage: 26, 34, and 38, respectively, according to Hamburger and Hamilton, 1951) [16]. The gonads were cultured for up to 90 min at 37.8 °C in phosphate-buffered saline without Ca2+ and Mg2+ (PBS [–]) to obtain live gonadal germ cells. The discharged, viable PGCs were identified by the morphological criteria and qPCR identification [13]. The PGCs have a typically spherical shape with a nucleus as the pale rounded structure. PGCs have numerous lipid droplets in their cytoplasm.

2.2. RNA and DNA Isolation

The PGCs were collected in a lysis buffer, and RNA was isolated by using a commercial kit for single-cell application according to the manufacturer’s protocol (GenElute Single Cell RNA Purification Kit, Sigma-Aldrich, St. Louis, MO, USA). The quality and quantity of the RNA were estimated with a 2100 Bioanalyzer instrument using the Agilent RNA 6000 Nano Kit (Agilent Technologies, Santa Clara, CA, USA). The RNA integrity number was assessed during the analysis. RNA of only high quality and purity was subjected to further molecular analysis.
For the methylation analysis, DNA was isolated using a commercial kit for the isolation of gDNA from materials that contained its trace amounts (Sherlock AX, A&A Biotechnology, Gdynia, Poland). The DNA was assessed using a 2% agarose gel and a Nanodrop2000 spectrophotometer (Thermo-Fisher Scientific, Waltham, MA, USA). PGC-related markers were examined to show that the isolated nucleic acid was derived from PGCs, as described by Dunislawska et al. [13].

2.3. Sex Determination of Embryos

Sexing of embryos for methylation analysis was performed by PCR. The primer sequences and reaction conditions for detecting W chromosome-specific repeating sequences and control PCR to amplify the chicken 18S ribosomal gene were chosen based on Clinton et al. [17]. W chromosome sequence was chosen to amplify a 415 bp product (primer sequences—F: CCCAAATATAACACGCTTCACT; R: GAAATGAATTATTTTCTGGCGAC) and ribosomal gene sequence was selected to amplify a 256 bp product (primer sequences—F: AGCTCTTTCTCGATTCCGTG; R: GGGTAGACACAAGCTGAGCC). The reaction was performed in the presence of Color Taq DNA polymerase (EURx, Gdansk, Poland). After amplification, the reaction product was separated on a 2% agarose gel and visualized in a gBOX instrument (fluorescence gel imaging, Syngene, Frederick, MD, USA) following staining with SimplySafe dye (EURx, Gdansk, Poland). Females were detected based on the presence of a band of 415 bp, whereas a band of 256 bp was detected in the control reaction.

2.4. Transcriptome Analysis

Microarray analysis was carried out without gender division (samples of randomly selected individuals) with 100 ng RNA using SurePrint G3 Custom GE 8 × 60 k microarrays for Gallus (n = 4/each group). Single-color microarrays were chosen because they allow a comparison of each sample with others and provide an immediate reference at any time. For this purpose, three slides were used with four matrices on each. The analysis within a time point was performed in four independent biological replicates. The microarray preparation procedure was performed following the manufacturer’s protocol for one-color microarray-based gene expression analysis (Agilent Technologies, Santa Clara, CA, USA), which includes the following stages: (1) sample preparation—template RNA with a spike-in, (2) cDNA synthesis, (3) cRNA synthesis and amplification, (4) cRNA purification, (5) preparation of the hybridization sample, (6) 17 h of hybridization at 65 °C, (7) microarray wash, (8) scan in Agilent Microarray Scanner System—SureScan (Agilent G3 GX1 color protocol), and (9) extraction of data using Agilent Feature Extraction Software.
After scanning and feature extraction, the obtained data were analyzed by GeneSpring GX software, version 14.9 (Agilent Technologies, Santa Clara, CA, USA). The data were normalized and subjected to quality control in the software. The following criteria were used to analyze the obtained list of genes: statistical significance (p value) higher than 0.05 and cut-off value greater than 2.0 indicating upregulated genes and lower than—2.0 indicating downregulated genes (cut-off of fold change value > 2 or <−2). Statistical analysis was performed using one-way analysis of variance. The selected time points were also compared with each other and between genotypes (GP vs. WL) as described [13]. The Agilent microarray was validated in Dunislawska et al. for WL [13]. For qualitative analysis, a heat map was used to present the general direction of changes in expression. GeneSpring software was also used for the qualitative assessment to assign genes to the main terms of gene ontology (GO).

2.5. Global Methylation Analysis

DNA for global methylation analysis was intended for analysis taking into account the division into sex of individuals after prior PCR identification and prepared using a commercial kit for methylated DNA quantification (MDQ1, Imprint Methylated DNA Quantification Kit, Sigma-Aldrich) according to the manufacturer’s protocol. DNA isolated from gPGCs (n = 5 for each time point and each sex) was diluted in a binding solution to a final concentration of 150 ng/µL. Then, it was used for estimating the level of methylated DNA based on enzyme-linked immunosorbent assay on 96-well plates. Positive (methylated) and blank controls were analyzed together with the DNA samples. The absorbance was measured at 450 nm by using the Multiscan (Thermo-Fisher Scientific, Waltham, MA, USA) microplate reader and read by SkanIt software 6.0.2. (Thermo-Fisher Scientific, Waltham, MA, USA). For each time point (4.5, 8, and 12 days), five samples per group (male/female), each derived from a different individual, were analyzed. The absorbance measurements were taken in two technical repetitions, each time using the same amount of DNA. The two measurements were averaged, and the mean value was used for the further analysis. Global DNA methylation levels were expressed as percentages relative to the methylated control and were calculated using the following formula: A 450 S A 450 B A 450 M C A 450 B × 100 % , where A 450 S is the average absorbance of the sample, A 450 B is the average absorbance of the blank, and A 450 M C is the average absorbance of the methylated control. Statistical analysis was carried out by SAS Enterprise Guide 8.2 (SAS Institute Inc., Cary, NC, USA). The quantitative values were first analyzed for normality using the Shapiro–Wilk test. Obtained values failed normal distribution assumption; therefore, the Kruskal–Wallis test was used (the influence of sex and the influence of developmental day; p < 0.05).

2.6. Gene-Specific Methylation Analysis

The isolated DNA was also subjected to the methylation analysis using the real-time quantitative methylation-specific polymerase chain reaction (qMSP-PCR) as described by Dunislawska et al. [9].
Based on the microarray data (presented in the current publication and our previous manuscript [13]), genes that have been silenced and whose functions were assumed to be important for embryonic development were selected for qMSP.
The conversion was carried out using the EpiJet Bisulfite Conversion Kit (Thermo-Fisher Scientific/Fermentas, Vilnius, Lithuania) according to the manufacturer’s instructions. qPCR was performed for the selected genes. For each gene, two primer pairs were designed—specific for methylated and nonmethylated DNA—using the MethPrimer tool according to Dunislawska et al. [9]. Primers for qMSP were complementary to the gene promoter region and were designed using the following criteria: island size > 100, GC% > 50.0, and obs./exp > 0.60. DNA oligonucleotides were synthesized by Sigma-Aldrich (St. Louis, MO, USA) (sequences presented in Table 1). The qPCR analysis was performed in the LightCycler 480 (Roche Diagnostics, Risch-Rothreuz, Switzerland) thermal cycler. The reaction mixture contained the Maxima SYBR Green qPCR Master Mix intercalating dye (Thermo-Fisher Scientific/Fermentas, Vilnius, Lithuania). The optimized melting temperature was 58 °C. The relative level of DNA methylation (%) was calculated from the melting curves (red fluorescence level) for each DNA sequence according to the formula described by Dunislawska et al. [9]. The quantitative values were first analyzed for normality using the Shapiro–Wilk test. Obtained values failed normal distribution assumption; therefore, the Kruskal–Wallis test was used (the influence of sex, influence of genotype, and the influence of embryonic developmental day; p < 0.05).

3. Results

3.1. Transcriptome Analysis

Comparative analysis showed changes in the expression (cut-off of fold change > 2 or <−2) in the GP only on day 8 of embryonic development. However, these differences between various time points of the GO embryo development were not statistically significant. The Volcano plot of the gene expression profiles at the studied time points are shown in Figure 2. The complete list of genes obtained in the microarray analysis is presented in Table S1 in supplementary file.
The numerical analysis of the gene expression levels data obtained from the microarray at the three time points is presented by heat maps in Figure 3.
Statistically significant differences (p < 0.05) between PGCs from both breeds, GP and WL (data available in [13]) were detected on day 8 of embryonic development. Strong regulation of the gene expression (cut-off >2 or <−2) was demonstrated in the PGCs of GP as compared to those of WL. Differences in the global expression on day 8 of embryo development are presented in Figure 4.
Genes significantly differentiating GP and WL were found to be mostly involved in the pathways related to, e.g., the regulation of biological and cellular processes, developmental processes, anatomical structure development, and multicellular organism development and organ development. Table S2 in supplementary file shows a list of pathways in which regulated genes are involved and the number of genes involved categorized by GeneSpring GX software. In addition, Table S3 provides a list of genes that have been included in the reported GO terms.

3.2. Global Methylation Analysis

There were no significant differences, but there were numerical differences between individual days and gender in the methylation analysis. In the case of the GP, the level of global methylation decreased with the following days of the embryonic development, which is especially noticeable in females. In WL, a decrease in methylation was observed on day 8, with a further increase on day 12. The results of global methylation analysis in gPGCs are presented in Figure 5.

3.3. Methylation of Silenced Genes

The selected genes for methylation analysis based on microarray data with statistically significant changes in the level of their expression are presented in Table 2.
There are statistically significant changes in the CYR61 gene methylation observed between genotypes (WL vs. GP) (p < 0.05) and between genders in WL (p < 0.01). CYR61 analyzed in the GP breed showed only a numerical (not statistical) increase in the methylation, especially in males. Statistically significant changes (p < 0.05) were also observed in the MIEN1 gene in WL between the different days of embryonic development. The overall changes in gene methylation levels are shown in Figure 6. These results are consistent with the gene expression levels. The greatest negative regulation was observed for the expression of the CYR61 gene both for the genotype and for the embryonic development day. In the case of methylation of the MIEN1 gene, a significant decrease in expression on day 12 in WL is consistent with an increase in methylation.

4. Discussion

In this study, we aimed to verify the hypothesis that changes at the level of gene expression during embryonic development may be nongenetic and are determined by the breed and gender. We made such a hypothesis on the basis of our previous analyses, where it was shown that the day of incubation is important for the transcriptomic profile between individual embryonic stages in WL. On the basis of the obtained results, we decided to verify not only whether such a relationship is also observed in other breeds but also whether silencing the expression of individual genes may be related to the level of their methylation.
Two very different chicken breeds were subjected to this study. One of them is a GP, an old native, dual-purpose Polish chicken, which became recognized as a breed at the end of the nineteenth century. These chickens can adapt well to extreme environmental conditions and hence are considered to be more robust and disease-resistant than other breeds. Currently, the GP chicken breed is protected from extinction by in situ preservation. GP chickens have been kept in a flock without subjecting them to any selection process for more than 60 generations [18]. Our previous studies suggested that a GP can be an excellent model for research, for example, as a good source of PGCs [19] and hence has been extensively used for biotechnological manipulations [20,21].
The WL breed, on the other hand, has been intensively selected for several decades mainly to improve reproductive traits and has been used widely to create highly productive egg-laying hybrids for commercial and industrial operations. Literature reports suggest different numbers of cPGCs and gPGCs in both genotypes.
Kuwana et al. indicated a large variation in the PGC concentration even within the same breed [22]. The results of Szczerba et al. [19] showed that the number of gPGCs was probably influenced by the stage of embryonic development. In the case of inbred lines or endangered breeds, the number of total PGCs seems to be lower compared to commercial chicken breeds [23,24]. This lower number of PGCs seems likely to be correlated with the lower productivity of the local chicken. This is in line with the report of Zhao et al. [25] that the number of chicken PGCs correlated with the number of eggs. Several studies of the past have indicated that the consequences of intensive genetic selection of chicken for economically important production traits, which has been carried out for almost a century, can be detected in the early stages of avian embryogenesis [2,4,26].

4.1. Transcriptome Analysis

Differences in gene expression were detected on day 8 of embryonic development. This is a very specific time point of the development of an embryo when the female PGCs in the germinal epithelium are more numerous than in the earlier stages. At this stage, the PGCs undergo a rapid mitotic cell division and gonadal sex differentiation in the ovarian cortex [27]. Moreover, day 8 is marked by a rapid growth of external tissues in the chicken embryo [28]. The results of transcriptomic analysis obtained from the present and previous studies [13] suggest a faster development of the GP embryos compared to the WL embryos.

4.2. Global Methylation Analysis

Research shows that during embryogenesis, right after zygote formation, DNA demethylation takes place (the reverse process of methylation, in which a methyl group is removed from the molecule), and a new methylation profile begins to be established de novo from the blastocyst stage [29]. Gryzinska et al. showed that methylation is age related. An increase in the level of global DNA methylation in Polbar embryos was demonstrated during days 6–18 of embryonic development [30]. These authors also proved that global methylation levels decrease in individuals in the postembryonic phase. In our study, this relationship was confirmed in WL breed. The literature shows that the level of methylation is species and tissue specific [9,31], whereas the pattern of this process is determined during embryogenesis.
In our study, global methylation analysis confirmed that day 8 is crucial and that DNA methylation may be breed dependent. The results of the analysis obtained for both WL and GP are in line with Jang et al.’s study [14]. These authors performed DNA methylation chip analysis of PGCs for both genders and reported that DNA methylation is evenly distributed among male and female PGCs.
In our previous studies, we also found differences in global methylation levels between two different in ovo-stimulated chicken genotypes in the early stage of embryonic development. The analyses were performed on two genotypes of different origins and selection history. The difference observed in the analyses suggested that selection history and genetic diversity play a key role in the tissue-specific methylation process [32].

4.3. Gene-Specific Methylation Analysis

Analysis of individual gene methylation between WL and GP clearly shows breed-related changes in the gene profile. In WL, numerically greater differences in the levels of gene methylation were observed between the different days of embryonic development. Interestingly, in the transcriptome analysis in GP, no significant differences in gene expression were observed between different days, which may also translate into the lack of differences between days in the context of methylation.
The CYR61 (cysteine-rich angiogenic inducer 61) gene is involved in the regulation of cell signaling and cell adhesion and in the developmental processes of the heart and skeleton. So far, the expression of this gene has been noted only in somatic cells [33], and its role in the development of gonads has not been reported. Our previous research has shown that this gene is silenced in the metabolic and immune tissues of adult chickens after stimulation of the intestinal microbiota during embryonic development [34]. However, it was observed that this silencing was not related to its methylation in the liver, spleen, or cecum tonsils in adult GP and broiler chicken [32]. Thus, silencing of gene expression in these cases may be related to RNA interference, and not DNA methylation. The present study showed numerically significant changes in the methylation levels of the CYR61 gene of PGCs in GP while also considering the influence of gender. This result may suggest that these changes are closely related to the previously discussed mechanism of methylation decrease during embryonic development. It also may have a biological basis, given the role of this gene in the development of the heart. In the development of chicken embryo, the heart is the first organ to function, and it starts to develop in the HH2 stage [35].
The MYH7 (myosin heavy-chain 7) gene encodes the β heavy-chain subunit of cardiac myosin. In our study, downregulation of the MYH7 expression in WL was observed in the analyzed days of embryonic development. Interestingly, we observed positive regulation of its expression in GP on the eighth day of development and a decrease in its methylation in GP, regardless of the gender. MYH7 has been described primarily in the context of hypertrophic cardiomyopathy in humans. However, chicken MYH7 has been shown to have sequence similarity and evolutionary association with mutant human MYH7. Karunanithi et al. have pointed out the anatomical similarities between the heart of chicken and humans. These authors highlighted that MYH7 is one of the proteins responsible for heart contractions, and, therefore, abnormalities in the gene encoding this protein can lead to excessive thickening of the myocardium [36]. It is necessary to conduct further analyses in this direction, because in the conducted experiment no phenotypic analyses in this area were performed.
The NEFL (neurofilament light) gene showed slightly upregulated expression in GP on day 8 in the present study, thus showing no significant differences in methylation. In WL, its expression was clearly negative whereas methylation differed between days—especially in terms of sex (a significant increase in expression in the following days in males). This result suggests that the expression-silencing mechanism may be closely related to gene methylation. The significant function of this gene has been described as a typical neuronal marker of brain development based on the analysis of chicken embryos by high-throughput sequencing [37].
The MIEN1 (Migration and invasion enhancer 1) gene is mainly involved in the regulation of apoptotic processes and cell migration. In WL, this gene shows a strong negative regulation during embryonic development between days 8 and 12, with a simultaneous increase in methylation from day 8. The level of its methylation in GP was found to be high for 4.5 days. Our previous studies showed that at the 14th stage of HH (i.e., about day 2) of embryonic development of GP, the peak total number of cPGCs and the maximum concentration were observed [19]. There was a decrease in concentration at a later stage, suggesting that cPGCs may be gradually trapped by capillary networks.
Our hypothesis that the GP embryo develops faster than the WL is supported by the statistically significant increase in methylation in WL observed between 4.5 and 8 days, while indicating a high level of methylation in GP from 4.5 days. In this case, methylation may be an important mechanism for silencing the expression of the genes that are responsible for cell migration.

5. Conclusions

PGCs allow a proper cell programming at the epigenetic level. Our research is the first step towards a comparative analysis of epigenetic mechanisms and their impact on gene expression silencing in two different genotypes of chickens. Based on our research, we can conclude that the transcriptomic profile differs according to the breed and stage of embryonic development. The analyses also confirmed that changes in the level of gene expression during embryonic development of WL and GP may be epigenetic and dependent on methylation. These changes are also determined by the sex and breed of the animal.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/genes12071090/s1, Table S1: List of genes whose differentiation was demonstrated on the expression microarray by comparing GP and WL on day 8 of embryonic development; Table S2: Biological pathways (GO terms) involving genes whose different expression was observed on day 8 of embryonic development between GP and WL; Table S3: List of genes that have been classified into GO terms by the GeneSpring software.

Author Contributions

Conceptualization, M.B.; methodology, A.D. and M.B., software, A.D.; validation, A.D., formal analysis, A.D.; investigation, A.D., resources, M.B.; data curation, A.D.; writing—original draft preparation, A.D., M.B. and M.S.; writing—review and editing, M.B., M.S. and K.S.; visualization, A.D.; supervision, M.B.; project administration, M.B.; funding acquisition, M.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by grant UMO-2017/27/B/NZ9/01510 funded by the National Science Center (Cracow, Poland).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

New data were created, analyzed, and presented in this study. Data sharing is not applicable to this article.

Acknowledgments

The authors wish to acknowledge Agata Szczerba, Marta Plewa, and Elzbieta Pietrzak for their excellent laboratory work and support with the isolation of PGCs and DNA.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Buzała, M.; Janicki, B.; Czarnecki, R. Consequences of different growth rates in broiler breeder and layer hens on embryogenesis, metabolism and metabolic rate: A review. Poult. Sci. 2014, 94, 728–733. [Google Scholar] [CrossRef]
  2. Bednarczyk, M.; Rosinski, A. Comparison of egg hatchability and in vitro survival of goose embryos of various origins. Poult. Sci. 1999, 78, 579–585. [Google Scholar] [CrossRef]
  3. Sawicka, D.; Samek, K.; Chojnacka-Puchta, L.; Witkowski, A.; Knaga, S.; Debowska, M.; Bednarczyk, M. Changes in quail blastodermal cell status as a result of selection. Folia Biol. 2015, 63, 63–67. [Google Scholar] [CrossRef] [PubMed]
  4. Sadeghi, R.; Pakdel, M.; Moradi Shabrbabak, M.; Moradi Shabrbabak, H. Divergent Selection Effect on Reproductive Trait in Japanese Quails. Res. J. Poult. Sci. 2013, 6, 18–22. [Google Scholar]
  5. Rashidi, H.; Sottile, V. The chick embryo: Hatching a model for contemporary biomedical research. BioEssays 2009, 31, 459–465. [Google Scholar] [CrossRef]
  6. Weeke-Klimp, A.; Bax, N.A.M.; Bellu, A.R.; Winter, E.M.; Vrolijk, J.; Plantinga, J.; Maas, S.; Brinker, M.; Mahtab, E.A.F.; Gittenberger-de Groot, A.C.; et al. Epicardium-derived cells enhance proliferation, cellular maturation and alignment of cardiomyocytes. J. Mol. Cell. Cardiol. 2010, 49, 606–616. [Google Scholar] [CrossRef] [PubMed]
  7. Cogburn, L.A.; Wang, X.; Carre, W.; Rejto, L.; Porter, T.E.; Aggrey, S.E.; Simon, J. Systems-wide chicken DNA microarrays, gene expression profiling, and discovery of functional genes. Poult. Sci. 2003, 82, 939–951. [Google Scholar] [CrossRef] [PubMed]
  8. Ribatti, D. Chicken chorioallantoic membrane angiogenesis model. Methods Mol. Biol. 2012, 843, 47–57. [Google Scholar]
  9. Dunislawska, A.; Slawinska, A.; Siwek, M. Hepatic DNA Methylation in Response to Early Stimulation of Microbiota with Lactobacillus Synbiotics in Broiler Chickens. Genes 2020, 11, 579. [Google Scholar] [CrossRef] [PubMed]
  10. Ginsburg, M.; Eyal-Giladi, H. Primordial Germ Cells of the Young Chick Blastoderm Originate from the Central Zone of the Area Pellucida Irrespective of the Embryo-Forming Process. Development 1987, 101, 209–219. [Google Scholar] [CrossRef]
  11. Chojnacka-Puchta, L.; Kasperczyk, K.; Płucienniczak, G.; Sawicka, D.; Bednarczyk, M. Primordial germ cells (PGCs) as a tool for creating transgenic chickens. Pol. J. Vet. Sci. 2012, 15, 181–188. [Google Scholar] [CrossRef] [Green Version]
  12. Mochizuki, K.; Matsui, Y. Epigenetic profiles in primordial germ cells: Global modulation and fine tuning of the epigenome for acquisition of totipotency. Dev. Growth Differ. 2010, 52, 517–525. [Google Scholar] [CrossRef]
  13. Dunislawska, A.; Szczerba, A.; Siwek, M.; Bednarczyk, M. Dynamics of the transcriptome during chicken embryo development based on primordial germ cells. BMC Res. Notes 2020, 13, 441. [Google Scholar] [CrossRef] [PubMed]
  14. Jang, H.J.; Seo, H.W.; Lee, B.R.; Yoo, M.; Womack, J.E.; Han, J.Y. Gene expression and DNA methylation status of chicken primordial germ cells. Mol. Biotechnol. 2013, 54, 177–186. [Google Scholar] [CrossRef]
  15. Nakajima, Y.; Minematsu, T.; Naito, M.; Tajima, A. A New Method for Isolating Viable Gonadal Germ Cells from 7-day-old Chick Embryos. J. Poult. Sci. 2011, 48, 106–111. [Google Scholar] [CrossRef] [Green Version]
  16. Hamburger, V.; Hamilton, H.L. A series of normal stages in the development of the chick embryo. J. Morphol. 1951, 88, 49–92. [Google Scholar] [CrossRef] [PubMed]
  17. Clinton, M.; Haines, L.; Belloir, B.; McBride, D. Sexing chick embryos: A rapid and simple protocol. Br. Poult. Sci. 2001, 42, 134–138. [Google Scholar] [CrossRef] [PubMed]
  18. Jawor, M.; Knaga, S.; Kozłowska, I.; Barna, J.; Váradi, É.; Kasperek, K.; Drobnyák, Á.; Bodzsár, N.; Várkonyi, E.P.; Jeżewska-Witkowska, G.; et al. Population Structure of Four Indigenous Chicken Breeds Undergoing In Situ Conservation. Anim. Sci. Pap. 2020, 38, 167–179. [Google Scholar]
  19. Szczerba, A.; Kuwana, T.; Bednarczyk, M. Concentration and total number of circulating primordial germ cells in Green-legged Partridgelike chicken embryos. Poult. Sci. 2021, 100, 319–324. [Google Scholar] [CrossRef]
  20. Siwek, M.; Wragg, D.; Sławińska, A.; Malek, M.; Hanotte, O.; Mwacharo, J.M. Insights into the genetic history of Green-legged Partridgelike fowl: MtDNA and genome-wide SNP analysis. Anim. Genet. 2013, 44, 522–532. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Bednarczyk, M.; Lakota, P.; Slomski, R.; Plawski, A.; Lipinski, D.; Siemieniako, B.; Lisowski, M.; Czekalski, P.; Grajewski, B.; Dluzniewska, P. Reconstitution of a chicken breed by inter se mating of germline chimeric birds. Poult. Sci. 2002, 81, 1347–1353. [Google Scholar] [CrossRef]
  22. Kuwana, T.; Kawashima, T.; Naito, M.; Yamashita, H.; Matsuzaki, M.; Takano, T. Conservation of a Threatened Indigenous Fowl (Kureko Dori) Using the Germline Chimeras Transplanted from Primordial Germ Cells. J. Poult. Sci. 2006, 43, 60–66. [Google Scholar] [CrossRef] [Green Version]
  23. Kostaman, T.; Yusuf, T.L.; Fahrudin, M.; Setiadi, M.A. Isolation and number of circulated primordial germ cells (circulated-PGCs) on stages of embryonic development of Gaok chicken. Jurnal Ilmu Ternak dan Veteriner 2013, 18, 54–62. [Google Scholar] [CrossRef]
  24. Sopiyana, S.; Setiadi, M.A.; Fahrudin, M.; Supriatna, I. Isolation and number of gonadal primordial germ cells (gonadal PGCs) on the stages of early embryonic development of KUB chicken. Media Peternak. 2017, 40, 1–6. [Google Scholar] [CrossRef] [Green Version]
  25. Zhao, D.-F.; Yamashita, H.; Matsuzaki, M.; Takano, T.; Abe, S.-I.; Naito, M.; Kuwana, T. Genetic Factors Affect the Number of Circulating Primordial Germ Cells in Early Chick Embryos. J. Poult. Sci. 2003, 40, 101–113. [Google Scholar] [CrossRef] [Green Version]
  26. Chojnacka-Puchta, L.; Sawicka, D.; Lakota, P.; Plucienniczak, G.; Bednarczyk, M.; Plucienniczak, A. Obtaining chicken primordial germ cells used for gene transfer: In vitro and in vivo results. J. Appl. Genet. 2015, 56, 493–504. [Google Scholar] [CrossRef] [Green Version]
  27. Kim, Y.M.; Han, J.Y. The early development of germ cells in chicken. Int. J. Dev. Biol. 2018, 62, 145–152. [Google Scholar] [CrossRef] [Green Version]
  28. Li, S.; Zhu, Y.; Zhi, L.; Han, X.; Shen, J.; Liu, Y.; Yao, J.; Yang, X. DNA Methylation Variation Trends during the Embryonic Development of Chicken. PLoS ONE 2016, 11, e0159230. [Google Scholar] [CrossRef] [Green Version]
  29. Messerschmidt, D.M.; Knowles, B.B.; Solter, D. DNA methylation dynamics during epigenetic reprogramming in the germline and preimplantation embryos. Genes Dev. 2014, 28, 812–828. [Google Scholar] [CrossRef] [Green Version]
  30. Gryzinska, M.; Blaszczak, E.; Strachecka, A.; Jezewska-Witkowska, G. Analysis of age-related global DNA methylation in chicken. Biochem. Genet. 2013, 51, 554–563. [Google Scholar] [CrossRef] [Green Version]
  31. Romanov, G.A.; Vanyushin, B.F. Methylation of reiterated sequences in mammalian DNAs Effects of the tissue type, age, malignancy and hormonal induction. BBA Sect. Nucleic Acids Protein Synth. 1981, 653, 204–218. [Google Scholar] [CrossRef]
  32. Dunislawska, A.; Slawinska, A.; Gryzinska, M.; Siwek, M. Interaction between early in ovo stimulation of the gut microbiota and chicken host—Splenic changes in gene expression and methylation. J. Anim. Sci. Biotechnol. 2021, 12, 73. [Google Scholar] [CrossRef]
  33. Piprek, R.P.; Kolasa, M.; Podkowa, D.; Kloc, M.; Kubiak, J.Z. Transcriptional profiling validates involvement of extracellular matrix and proteinases genes in mouse gonad development. Mech. Dev. 2018, 149, 9–19. [Google Scholar] [CrossRef] [PubMed]
  34. Dunislawska, A.; Slawinska, A.; Bednarczyk, M.; Siwek, M. Transcriptome modulation by in ovo delivered Lactobacillus synbiotics in a range of chicken tissues. Gene 2019, 698, 27–33. [Google Scholar] [CrossRef] [PubMed]
  35. Martinsen, B.J. Reference guide to the stages of chick heart embryology. Dev. Dyn. 2005, 233, 1217–1237. [Google Scholar] [CrossRef] [PubMed]
  36. Karunanithi, S.; Venkatachalapathy, L.; Sudhalakshmi, Y. Comparative In Silico Analysis of Hypertrophic Cardiomyopathy Heart in Human and Normal Chicken Heart. J. Adv. Lab. Res. Biol. 2010, 1, 60–63. [Google Scholar]
  37. Xu, Z.; Che, T.; Li, F.; Tian, K.; Zhu, Q.; Mishra, S.K.; Dai, Y.; Li, M.; Li, D. The temporal expression patterns of brain transcriptome during chicken development and ageing. BMC Genom. 2018, 19, 917. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Experimental set-up (created with BioRender.com).
Figure 1. Experimental set-up (created with BioRender.com).
Genes 12 01090 g001
Figure 2. Volcano plot of genes expressed at the three time points of embryonic development of Green-legged Partridge-like (GP) chicken: 4.5, 8, and 12 days (red dots: upregulated genes; blue dots: downregulated genes statistically significant; gray dots: not statistically significant).
Figure 2. Volcano plot of genes expressed at the three time points of embryonic development of Green-legged Partridge-like (GP) chicken: 4.5, 8, and 12 days (red dots: upregulated genes; blue dots: downregulated genes statistically significant; gray dots: not statistically significant).
Genes 12 01090 g002
Figure 3. Heat maps of genes expressed at three time points of embryonic development of GP chicken: 4.5, 8, and 12 days.
Figure 3. Heat maps of genes expressed at three time points of embryonic development of GP chicken: 4.5, 8, and 12 days.
Genes 12 01090 g003
Figure 4. Heat maps of statistically significant differences in gene expression between two genotypes (comparison of GP to White Leghorn (WL)) on day 8 of embryonic development.
Figure 4. Heat maps of statistically significant differences in gene expression between two genotypes (comparison of GP to White Leghorn (WL)) on day 8 of embryonic development.
Genes 12 01090 g004
Figure 5. Global methylation levels (%) in three stages of embryonic development (days 4.5, 8, and 12) of two chicken breeds: WL and GP and two sexes (male and female).
Figure 5. Global methylation levels (%) in three stages of embryonic development (days 4.5, 8, and 12) of two chicken breeds: WL and GP and two sexes (male and female).
Genes 12 01090 g005
Figure 6. DNA methylation of the CYR61, MYH7, NEFL, and MIEN1 genes in gonadal PGCs of WL and GP in three stages of embryonic development (days 4.5, 8, and 12) and categorized by sex. X-axis: days of embryonic development on which PGCs are collected; Y-axis: percentage of gene methylation.
Figure 6. DNA methylation of the CYR61, MYH7, NEFL, and MIEN1 genes in gonadal PGCs of WL and GP in three stages of embryonic development (days 4.5, 8, and 12) and categorized by sex. X-axis: days of embryonic development on which PGCs are collected; Y-axis: percentage of gene methylation.
Genes 12 01090 g006
Table 1. Sequences of the primers designed for qMSP reaction by using the MethPrimer tool.
Table 1. Sequences of the primers designed for qMSP reaction by using the MethPrimer tool.
Gene (NCBI No.)Primer SequencesGC%Amplicon
Size
Reference
CYR61
(429089)
MF: TTTGGTTTTAGTGTTTAAAGACGT
R: TTATATTTACCTTCAAAAAAACGTA
58.33
44.00
150[9]
UF: TTTTGGTTTTAGTGTTTAAAGATGT
R: TATTTATATTTACCTTCAAAAAAACATA
56.00
42.86
154
MYH7
(395350)
MF: AGGGTTTTGTTTCGTGTTTTATTC
R: CTCCCCCATCTCTATAATAACGAT
70.83
62.50
100This study
UF: GGGTTTTGTTTTGTGTTTTATTTGT
R: CCTCCCCCATCTCTATAATAACAAT
76.00
64.00
100
NEFL
(419528)
MF: TTTTTTGTATTCGGTGGATAGTTTC
R: TAAAATCCTACAACTAAACCCGCT
68.00
70.83
104This study
UF: TTTTTGTATTTGGTGGATAGTTTTG
R: TTAAAATCCTACAACTAAACCCACT
68.00
68.00
104
MIEN1
(100858225)
MF: GGGGTAGTTGAGAGTTATACGT
R:TACAAAATAATACAAAAAAAACGAC
68.18
68.00
125This study
UF: TGTGGGGTAGTTGAGAGTTATATGT
R: TACAAAATAATACAAAAAAAACAAC
64.00
68.00
128
M: specific for methylated DNA; U: specific for unmethylated DNA.
Table 2. Statistically significant gene expression changes detected on the microarray (GeneSpring software) in particular points of embryonic development of two chicken breeds—WL and GP.
Table 2. Statistically significant gene expression changes detected on the microarray (GeneSpring software) in particular points of embryonic development of two chicken breeds—WL and GP.
BreedDayCYR61MYH7NEFLMIEN1
WL8 vs. 4.5–3.26–2.37ndnd
12 vs. 4.5–2.34–3.91–4.19–2.52
GP8 vs. 4.5–1.781.300.37nd
WL: White Leghorn; GP: Green-legged Partridgelike; nd: no statistically significant data.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Dunislawska, A.; Siwek, M.; Stadnicka, K.; Bednarczyk, M. Comparison of the Transcriptomic and Epigenetic Profiles of Gonadal Primordial Germ Cells of White Leghorn and Green-Legged Partridgelike Chicken Embryos. Genes 2021, 12, 1090. https://0-doi-org.brum.beds.ac.uk/10.3390/genes12071090

AMA Style

Dunislawska A, Siwek M, Stadnicka K, Bednarczyk M. Comparison of the Transcriptomic and Epigenetic Profiles of Gonadal Primordial Germ Cells of White Leghorn and Green-Legged Partridgelike Chicken Embryos. Genes. 2021; 12(7):1090. https://0-doi-org.brum.beds.ac.uk/10.3390/genes12071090

Chicago/Turabian Style

Dunislawska, Aleksandra, Maria Siwek, Katarzyna Stadnicka, and Marek Bednarczyk. 2021. "Comparison of the Transcriptomic and Epigenetic Profiles of Gonadal Primordial Germ Cells of White Leghorn and Green-Legged Partridgelike Chicken Embryos" Genes 12, no. 7: 1090. https://0-doi-org.brum.beds.ac.uk/10.3390/genes12071090

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop