Quantitative evaluation of DNA methylation patterns for ALVE and TVB genes in a neoplastic disease susceptible and resistant chicken model

Transgenerational epigenetic inheritance and immunity in chickens that vary in Marek's disease resistance

Marek's disease virus (MDV), a naturally oncogenic, highly contagious alpha herpesvirus, induces a T cell lymphoma in chickens that causes severe economic loss. Marek's disease (MD) outcome in an individual is attributed to genetic and environmental factors. Further investigation of the host-virus interaction mechanisms that impact MD resistance is needed to achieve greater MD control. This study analyzed genome-wide DNA methylation patterns in 2 highly inbred parental lines 63 and 72 and 5 recombinant congenic strains (RCS) C, L, M, N, and X strains from those parents. Lines 63 and 72, are MD resistant and susceptible, respectively, whereas the RCS have different combinations of 87.5% Line 63 and 12.5% Line 72. Our DNA methylation cluster showed a strong association with MD incidence. Differentially methylated regions (DMRs) between the parental lines and the 5 RCS were captured. MD-resistant and MD-susceptible markers of DNA methylation were identified as transgenerational epigenetic inheritable. In addition, the growth of v-src DNA tumors and antibody response against sheep red blood cells differed among the 2 parental lines and the RCS. Overall, our results provide very solid evidence that DNA methylation patterns are transgenerational epigenetic inheritance (TEI) in chickens and also play a vital role in MD tumorigenesis and other immune responses; the specific methylated regions may be important modulators of general immunity.

Association of Aberrant DNA Methylation Level in the CD4 and JAK-STAT-Pathway-Related Genes with Mastitis Indicator Traits in Chinese Holstein Dairy Cattle

The present study was designed to evaluate the gene expression and DNA methylation level in the promoter region of the CD4 and the JAK-STAT-pathway-related genes. A total of 24 samples were deployed in the gene expression and 118 samples were used in the DNA methylation study. Student's t-tests were used to analyze the gene expression and DNA methylation. The evaluation of DNA methylation in promoter regions of JAK2 and STAT5A revealed hypo-methylation levels of CpG sites and higher gene expression in cows diagnosed with mastitis as compared to the healthy control, and vice versa in those with CD4. DNA methylation was negatively correlated with gene expression in JAK2, STAT5A, and CD4 genes. Six, two, and four active transcription factors were identified on the CpG sites in the promoter regions of JAK2, STAT5A, and CD4 genes, respectively. Regarding correlation analysis, the DNA methylation levels of CD4 showed significantly higher positive correlations with somatic cell counts (p < 0.05). Findings of the current study inferred that aberrant DNA methylation in the CpG sites at the 1 kb promoter region in JAK2, STAT5A, and CD4 genes due to mastitis in cows can be used as potential epigenetic markers to estimate bovine mastitis susceptibility in dairy cattle.

The impact of endogenous Avian Leukosis Viruses (ALVE) on production traits in elite layer lines

Avian Leukosis Virus subgroup E (ALVE) integrations are endogenous retroviral elements found in the chicken genome. The presence of ALVE has been reported to have negative impacts on multiple traits, including egg production and body weight. The recent development of rapid, inexpensive and specific ALVE detection methods has facilitated their characterization in elite commercial egg production lines across multiple generations. The presence of 20 ALVE was examined in 8 elite lines, from 3 different breeds. Seventeen of these ALVE (85%) were informative and found to be segregating in at least one of the lines. To test for an association between specific ALVE inserts and traits, a large genotype by phenotype study was undertaken. Genotypes were obtained for 500 to 1500 males per line, and the phenotypes used were sire-daughter averages. Phenotype data were analyzed by line with a linear model that included the effects of generation, ALVE genotype and their interaction. If genotype effect was significant, the number of ALVE copies was fitted as a regression to estimate additive ALVE gene substitution effect. Significant associations between the presence of specific ALVE inserts and 18 commercially relevant performance and egg quality traits, including egg production, egg weight and albumen height, were observed. When an ALVE was segregating in more than one line, these associations did not always have the same impact (negative, positive or none) in each line. It is hypothesized that the presence of ALVE in the chicken genome may influence production traits by 3 mechanisms: viral protein production may modulate the immune system and impact overall production performance (virus effect); insertional mutagenesis caused by viral integration may cause direct gene alterations or affect gene regulation (gene effect); or the integration site may be within or adjacent to a quantitative trait region which impacts a performance trait (linkage disequilibrium, marker effect).

Identification of a novel differentially methylated region adjacent to ATG16L2 in lung cancer cells using methyl-CpG binding domain protein-enriched genome sequencing

Lung cancer is the most common cancer worldwide. Epigenetic modifications like DNA methylation play fundamental roles in the dynamic process of lung cancer. The objective of this study was to use methyl-CpG binding domain protein-enriched genome sequencing (MBD-Seq) to identify novel and high-confidence DNA methylation in lung tumor. We first compared the whole-genome DNA methylation of three lung cancer cell lines, including A549, H1299, and SK-MES-1, against BEAS-2B, a lung/bronchial normal epithelial cell line. We then used pyrosequencing and OneStep qMethyl kit methods to verify the results in the cell line specimens. MBD-Seq identified 279, 8046, and 22 887 differentially methylated regions (DMRs), respectively, with 120 common DMRs among three comparison groups. Three DMRs were consistent with the MBD-Seq results by both pyrosequencing and OneStep qMethyl validations. Furthermore, OneStep qMethyl kit was also performed for functional validation of these three potential DMRs in sputum DNA from clinical participants. We successfully identified one new DMR adjacent to ATG16L2. The novel DMR might have an important function in lung carcinogenesis. Further validation of the finding in clinical specimens of lung cancer patients and functional analysis of this novel DMR in the development of lung cancer through transcriptional silencing of ATG16L2 are warranted.

Genome-wide DNA methylation pattern in a mouse model reveals two novel genes associated with Staphylococcus aureus mastitis

Objective

Staphylococcus aureus (S. aureus) is one of the major microorganisms responsible for subclinical mastitis in dairy cattle. The present study was designed with the aim to explore the DNA methylation patterns using the Fluorescence-labeled methylation-sensitive amplified polymorphism (F-MSAP) techniques in a S. aureus-infected mouse model.

Methods

A total of 12 out-bred Institute of Cancer Research female mice ranging from 12 to 13 weeks-old were selected to construct a mastitis model. F-MSAP analysis was carried out to detect fluctuations of DNA methylation between control group and S. aureus mastitis group.

Results

Visible changes were observed in white cell counts in milk, percentage of granulocytes, percentage of lymphocytes, CD4⁺/CD8⁺ ratio (CD4⁺/CD8⁺), and histopathology of mice pre- and post-challenge with S. aureus. These findings showed the suitability of the S. aureus-infected mouse model. A total of 369 fragments was amplified from udder tissue samples from the two groups (S. aureus-infected mastitis group and control group) using eight pairs of selective primers. Results indicated that the methylation level of mastitis mouse group was higher than that in the control group. In addition, NCK-associated protein 5 (Nckap5) and transposon MTD were identified to be differentially methylated through secondary polymerase chain reaction and sequencing in the mastitis group. These observations might play an important role in the development of S. aureus mastitis.

Conclusion

Collectively, our study suggests that the methylation modification in Nckap5 and transposon MTD might be considered as epigenetic markers in resistance to S. aureus-infected mastitis and provided a new insight into S. aureus mastitis research in dairy industry and public health.

High-throughput sequencing of pituitary and hypothalamic microRNA transcriptome associated with high rate of egg production

Background

MicroRNAs exist widely in viruses, plants and animals. As endogenous small non-coding RNAs, miRNAs regulate a variety of biological processes. Tissue miRNA expression studies have discovered numerous functions for miRNAs in various tissues of chicken, but the regulation of miRNAs in chicken pituitary and hypothalamic development related to high and low egg-laying performance has remained unclear.

Results

In this study, using high-throughput sequencing technology, we sequenced two tissues (pituitary and hypothalamus) in 3 high- and 3 low-rate egg production Luhua chickens at the age of 300 days. By comparing low- and high-rate egg production chickens, 46 known miRNAs and 27 novel miRNAs were identified as differentially expressed (P < 0.05). Six differentially expressed known miRNAs, which are expressed in both tissues, were used in RT-qPCR validation and SNP detection. Among them, seven SNPs in two miRNA precursors (gga-miR-1684a and gga-miR-1434) were found that might enhance or reduce the production of the mature miRNAs. In addition, 124 and 30 reciprocally expressed miRNA-target pairs were identified by RNA-seq in pituitary and hypothalamic tissues, respectively and randomly selected candidate miRNA and miRNA-target pairs were validated by RT-qPCR in Jiuyuan black fowl. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway annotation illustrated that a large number of egg laying-related pathways were enriched in the high-rate egg production chickens, including ovarian steroidogenesis and steroid hormone biosynthesis.

Conclusions

These differentially expressed miRNAs and their predicted target genes, especially identified reciprocally expressed miRNA-target pairs, advance the study of miRNA function and egg production associated miRNA identification. The analysis of the miRNA-related SNPs and their effects provided insights into the effects of SNPs on miRNA biogenesis and function. The data generated in this study will further our understanding of miRNA regulation mechanisms in the chicken egg-laying process.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3644-3) contains supplementary material, which is available to authorized users.

Expressed microRNA associated with high rate of egg production in chicken ovarian follicles

MicroRNA (miRNA) is a highly conserved class of small noncoding RNA about 19-24 nucleotides in length that function in a specific manner to post-transcriptionally regulate gene expression in organisms. Tissue miRNA expression studies have discovered a myriad of functions for miRNAs in various aspects, but a role for miRNAs in chicken ovarian tissue at 300 days of age has not hitherto been reported. In this study, we performed the first miRNA analysis of ovarian tissues in chickens with low and high rates of egg production using high-throughput sequencing. By comparing low rate of egg production chickens with high rate of egg production chickens, 17 significantly differentially expressed miRNAs were found (P < 0.05), including 11 known and six novel miRNAs. We found that all 11 known miRNAs were involved mainly in pathways of reproduction regulation, such as steroid hormone biosynthesis and dopaminergic synapse. Additionally, expression profiling of six randomly selected differentially regulated miRNAs were validated by quantitative real-time polymerase chain reaction (RT-qPCR). Some miRNAs, such as gga-miR-34b, gga-miR-34c and gga-miR-216b, were reported to regulate processes such as proliferation, cell cycle, apoptosis and metastasis and were expressed differentially in ovaries of chickens with high rates of egg production, suggesting that these miRNAs have an important role in ovary development and reproductive management of chicken. Furthermore, we uncovered that a significantly up-regulated miRNA-gga-miR-200a-3p-is ubiquitous in reproduction-regulation-related pathways. This miRNA may play a special central role in the reproductive management of chicken, and needs to be further studied for confirmation.

Expression patterns of endogenous avian retrovirus ALVE1 and its response to infection with exogenous avian tumour viruses

Endogenous retroviruses (ERVs) are genomic elements that are present in a wide range of vertebrates and have been implicated in a variety of human diseases, including cancer. However, the characteristic expression patterns of ERVs, particularly in virus-induced tumours, is not fully clear. DNA methylation was analysed by bisulfite pyrosequencing, and gene expression was analysed by RT-qPCR. In this study, we first found that the endogenous avian retrovirus ALVE1 was highly expressed in some chicken tissues (including the heart, bursa, thymus, and spleen) at 2 days of age, but its expression was markedly decreased at 35 days of age. In contrast, the CpG methylation level of ALVE1 was significantly lower in heart and bursa at 2 days than at 35 days of age. Moreover, we found that the expression of ALVE1 was significantly inhibited in chicken embryo fibroblast cells (CEFs) and MSB1 cells infected with avian leukosis virus subgroup J (ALVJ) and reticuloendotheliosis virus (REV) at the early stages of infection. In contrast, the expression of the ALVE1 env gene was significantly induced in CEFs and MSB1 cells infected with Marek's disease virus (MDV). However, the methylation and expression levels of the ALVE1 long terminal repeat (LTR) did not show obvious alterations in response to viral infection. The present study revealed the expression patterns of ALVE1 in a variety of chicken organs and tissues and in chicken cells in response to avian tumour virus infection. These findings may be of significance for understanding the role and function of ERVs that are present in the host genome.

Epigenetics and inheritance of phenotype variation in livestock

Epigenetic inheritance plays a crucial role in many biological processes, such as gene expression in early embryo development, imprinting and the silencing of transposons. It has recently been established that epigenetic effects can be inherited from one generation to the next. Here, we review examples of epigenetic mechanisms governing animal phenotype and behaviour, and we discuss the importance of these findings in respect to animal studies, and livestock in general. Epigenetic parameters orchestrating transgenerational effects, as well as heritable disorders, and the often-overlooked areas of livestock immunity and stress, are also discussed. We highlight the importance of nutrition and how it is linked to epigenetic alteration. Finally, we describe how our understanding of epigenetics is underpinning the latest cancer research and how this can be translated into directed efforts to improve animal health and welfare.

Epigenetics and inheritance of phenotype variation in livestock

Epigenetic inheritance plays a crucial role in many biological processes, such as gene expression in early embryo development, imprinting and the silencing of transposons. It has recently been established that epigenetic effects can be inherited from one generation to the next. Here, we review examples of epigenetic mechanisms governing animal phenotype and behaviour, and we discuss the importance of these findings in respect to animal studies, and livestock in general. Epigenetic parameters orchestrating transgenerational effects, as well as heritable disorders, and the often-overlooked areas of livestock immunity and stress, are also discussed. We highlight the importance of nutrition and how it is linked to epigenetic alteration. Finally, we describe how our understanding of epigenetics is underpinning the latest cancer research and how this can be translated into directed efforts to improve animal health and welfare.

Temporal expression and DNA hypomethylation profile of CD30 in Marek's disease virus-infected chicken spleens

Marek's disease (MD) is a viral neoplastic disease of chickens caused by Marek's disease virus (MDV), which is serious threat to worldwide poultry industry. Our previous studies showed that the CD30 gene was hypomethylated in MD lymphoma. In this study, we further analyzed differential expression patterns and methylation levels of the CD30 gene between MDV-infected and noninfected spleens at 4, 7, 14, 21, and 28 d postinfection (dpi). The results showed that the expression of CD30 in MDV-infected spleens was significantly lower than that in noninfected spleens at 4 dpi. The expression of CD30 did not present significant difference between MDV-infected and noninfected spleens at 7 and 14 dpi. However, an increased expression of CD30 was presented in MDV-infected spleens at both 21 and 28 dpi. Simultaneously, CD30 showed a lower DNA methylation level in MDV-infected spleens at 14, 21, and 28 dpi. The results indicated that CD30 gene was involved in the whole process of MD tumorigenesis and upregulated expression of CD30 in MDV-infected spleens might be attributed to the hypomethylation of promoter of CD30 gene.

Transcriptional profiles of host-pathogen responses to necrotic enteritis and differential regulation of immune genes in two inbreed chicken lines showing disparate disease susceptibility

Necrotic enteritis (NE) is an important intestinal infectious disease of commercial poultry flocks caused by Clostridium perfringens. Using an experimental model of NE involving co-infection with C. perfringens and Eimeria maxima, transcriptome profiling and functional genomics approaches were applied to identify the genetic mechanisms that might regulate the host response to this disease. Microarray hybridization identified 1,049 transcripts whose levels were altered (601 increased, 448 decreased) in intestinal lymphocytes from C. perfringens/E. maxima co-infected Ross chickens compared with uninfected controls. Five biological functions, all related to host immunity and inflammation, and 11 pathways were identified from this dataset. To further elucidate the role of host genetics in NE susceptibility, two inbred chicken lines, ADOL line 6 and line 7 which share an identical B² major histocompatibility complex haplotype but differ in their susceptibility to virus infection, were compared for clinical symptoms and the expression levels of a panel of immune-related genes during experimental NE. Line 6 chickens were more susceptible to development of experimental NE compared with line 7, as revealed by decreased body weight gain and increased E. maxima oocyst shedding. Of 21 immune-related genes examined, 15 were increased in C. perfringens/E. maxima co-infected line 6 vs. line 7 chickens. These results suggest that immune pathways are activated in response to experimental NE infection and that genetic determinants outside of the chicken B complex influence resistance to this disease.

DNA methylation and regulation of the CD8A after duck hepatitis virus type 1 infection

Background

Cluster of differentiation 8 (CD8) is expressed in cytotoxic T cells, where it functions as a co-receptor for the T-cell receptor by binding to major histocompatibility complex class I (MHCI) proteins, which present peptides on the cell surface. CD8A is critical for cell-mediated immune defense and T-cell development. CD8A transcription is controlled by several cis-acting elements and trans-acting elements and is also regulated by DNA methylation. However, the epigenetic regulation of CD8A in the duck and its relationship with virus infection are still unclear.

Results

We investigated the epigenetic transcriptional regulatory mechanisms, such as DNA methylation, for the expression of the CD8A and further evaluated the contribution of such epigenetic regulatory mechanisms to DHV-I infection in the duck. Real-time quantitative polymerase chain reaction (RT-qPCR) revealed the highest level of CD8A expression to be in the thymus, followed by the lungs, spleen, and liver, and the levels of CD8A expression were very low in the kidney, cerebrum, cerebellum, and muscle in the duck. RT-qPCR also demonstrated that the CD8A mRNA was down-regulated significantly in morbid ducklings treated with DHV-1 and up-regulated significantly in non-morbid ducklings in all the tissues tested. In addition, hypermethylation of CD8A was detected in the morbid ducklings, whereas relatively low methylation of CD8A was evident in the non-morbid ducklings. The CD8A mRNA level was negatively associated with the CpG methylation level of CD8A and global methylation status.

Conclusions

We concluded that the mRNA level of the CD8A was negatively associated with the CpG methylation level of CD8A and global methylation status in the duck, suggesting that the hypermethylation of CD8A may be associated with DHV-1 infection. The first two CpG sites of the CD8A promoter region could be considered as epigenetic biomarkers for resistance breeding against duckling hepatitis disease in the duck.

Histone methylation analysis and pathway predictions in chickens after MDV infection

Marek's disease (MD) is a lymphoproliferative disease in chicken induced by Marek's disease virus (MDV). Although studies have focused on the genetic differences between the resistant and susceptible chicken, less is known about the role of epigenetic factors in MD. In this study, genome-wide histone modifications in the non-MHC-associated resistant and susceptible chicken lines were examined. We found that tri-methylation at histone H3 Lys4 (H3K4me3) enrichment is positively correlated with the expression of protein coding genes as well as microRNA (miRNA) genes, whereas tri-methylation at histone H3 Lys27 (H3K27me3) exhibits a negative correlation. By identifying line-specific histone modifications in MDV infection, we found unique H3K4me3 islands in the resistant chicken activated genes, which are related to immune response and cell adhesion. Interestingly, we also found some miRNAs from unique H3K27me3 patterns in the susceptible chickens that targeted genes involved in 5-hydroxytryptamine (5-HT)-receptor and adrenergic receptor pathways. In conclusion, dynamic line-specific histone modifications in response to MDV infection suggested that intrinsic epigenetic mechanisms may play a role in MD-resistance and -susceptibility.

Functional genomic analysis of variation on beef tenderness induced by acute stress in angus cattle

Beef is one of the leading sources of protein, B vitamins, iron, and zinc in human food. Beef palatability is based on three general criteria: tenderness, juiciness, and flavor, of which tenderness is thought to be the most important factor. In this study, we found that beef tenderness, measured by the Warner-Bratzler shear force (WBSF), was dramatically increased by acute stress. Microarray analysis and qPCR identified a variety of genes that were differentially expressed. Pathway analysis showed that these genes were involved in immune response and regulation of metabolism process as activators or repressors. Further analysis identified that these changes may be related with CpG methylation of several genes. Therefore, the results from this study provide an enhanced understanding of the mechanisms that genetic and epigenetic regulations control meat quality and beef tenderness.

DNA Methylation Fluctuation Induced by Virus Infection Differs between MD-resistant and -susceptible Chickens

Marek’s disease (MD) is a lymphoproliferative disease induced by Marek’s disease virus (MDV) infection. To augment vaccination measures in MD control, host genetic resistant to MD becomes obviously more and more important. To elucidate the mechanism of MD-resistance, most of researches were focused on the genetic differences between resistant and susceptible chickens. However, epigenetic features between MD resistant and susceptible chickens are poorly characterized. Using bisulfite pyrosequencing method, we found some candidate genes have higher promoter methylation in the MD-susceptible (L7₂) chickens than in the MD-resistant (L6₃) chickens. The hypermethylated genes, involved in cellular component organization, responding to stimulus, cell adhesion, and immune system process, may play important role in susceptibility to disease by deregulation of these genes. MDV infection induced the expression changes of all three methyltransferases genes (DNMT1, DNMT3a, and DNMT3b) in both lines of chickens. The DNMT1 was up-regulated in L7₂, whereas the DNMT3b was down-regulated in L6₃ at 21 dpi. Interestingly, a dynamic change of promoter methylation was observed during MDV life cycle. Some genes, including HDAC9, GH, STAT1, CIITA, FABP3, LATS2, and H2Ac, showed differential methylation behaviors between the two lines of chickens. In summary, the findings from this study suggested that DNA methylation heterogeneity and MDV infection induced methylation alterations differences existed between the two lines of chickens. Therefore, it is suggested that epigenetic mechanisms may be involved in modulating the resistance and/or susceptibility to MD in chickens.

Molecular identification of avian leukosis virus subgroup E loci and tumor virus B locus in Chinese indigenous chickens

Avian leukosis virus (ALV) subgroup E (ALVE) is an endogenous retrovirus in the chicken genome. The chickens carrying ALVE locus 3 (ALVE3), 6 (ALVE6), 9 (ALVE9), and 21 (ALVE21) have been proved to be susceptible to ALV. Tumor virus locus B (TVB) encodes the cellular receptor for ALV subgroups B, D, and E. The insertions of the 4 ALVE loci and the genotypes of TVB have not been demonstrated in Chinese indigenous chicken breeds. In the present study, the existence of ALVE3, ALVE6, ALVE9, and ALVE21 were detected in 10 native breeds of Chinese chickens and an introduced breed, the White Leghorn (2 populations in this study, WL1 and WL2), by locus-specific PCR. The PCR products of ALVE were further confirmed by sequencing assay. We also surveyed the status of genotypes of TVB in Silkie, Beijing You, and White Leghorn (WL1 and WL2) chickens with pyrosequencing assays. The results showed that the carrier frequency of ALVE3 was 1.3% in the Chinese chicken population, and was 10.3 in WL1 and 49.2% in WL2. The carrier frequency of ALVE6 was 5.4% in native breeds of Chinese birds, in contrast with 0% in WL1 and 6.8% in WL2. The carrier frequency of ALVE9 was 0.1% in the Chinese indigenous population, and was 16.0% in WL1 and 11.9% in WL2. The carrier frequency of ALVE21 was 10.4% in Chinese chickens, whereas ALVE21 was detected with a frequency of 0% in WL1 and 50% in WL2. The frequency of the TVB resistance allele (TVB*R and TVB*R') was 0.4% in Beijing You chickens, whereas it was 70.5% in WL1 and 54.5% in WL2. No carriers of ALVE3, ALVE9, and ALVE21 were detected in Silkie fowl, a famous Chinese native breed that has been used as a source for alternative medicine. These results present molecular evidence of ALVE3, ALVE6, ALVE9, and ALVE21 insertions and TVB genotypes in Chinese indigenous chickens and could provide potential molecular insights into anti-ALV breeding in chickens.

An integrated epigenetic and genetic analysis of DNA methyltransferase genes (DNMTs) in tumor resistant and susceptible chicken lines

Both epigenetic alterations and genetic variations play essential roles in tumorigenesis. The epigenetic modification of DNA methylation is catalyzed and maintained by the DNA methyltransferases (DNMT3a, DNMT3b and DNMT1). DNA mutations and DNA methylation profiles of DNMTs themselves and their relationships with chicken neoplastic disease resistance and susceptibility are not yet defined. In the present study, we analyzed the complexity of the DNA methylation variations and DNA mutations in the first exon of three DNMTs genes over generations, tissues, and ages among chickens of two highly inbred White Leghorn lines, Marek's disease-resistant line 6(3) and -susceptible line 7(2), and six recombinant congenic strains (RCSs). Among them, tissue-specific methylation patterns of DNMT3a were disclosed in spleen, liver, and hypothalamus in lines 6(3) and 7(2). The methylation level of DNMT3b on four CpG sites was not significantly different among four tissues of the two lines. However, two line-specific DNA transition mutations, CpG-->TpG (Chr20:10203733 and 10203778), were discovered in line 7(2) compared to the line 6(3) and RCSs. The methylation contents of DNMT1 in blood cell showed significant epimutations in the first CpG site among the two inbred lines and the six RCSs (P<0.05). Age-specific methylation of DNMT1 was detected in comparisons between 15 month-old and 2 month-old chickens in both lines except in spleen samples from line 7(2). No DNA mutations were discovered on the studied regions of DNMT1 and DNMT3a among the two lines and the six RCSs. Moreover, we developed a novel method that can effectively test the significance of DNA methylation patterns consisting of continuous CpG sites. Taken together, these results highlight the potential of epigenetic alterations in DNMT1 and DNMT3a, as well as the DNA mutations in DNMT3b, as epigenetic and genetic factors to neoplastic diseases of chickens.