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

Epigenetic Factor MicroRNAs Likely Mediate Vaccine Protection Efficacy against Lymphomas in Response to Tumor Virus Infection in Chickens through Target Gene Involved Signaling Pathways

Epigenetic factors, including microRNAs (miRNAs), play an important role in affecting gene expression and, therefore, are involved in various biological processes including immunity protection against tumors. Marek's disease (MD) is a highly contagious disease of chickens caused by the MD virus (MDV). MD has been primarily controlled by vaccinations. MD vaccine efficacy might, in part, be dependent on modulations of a complex set of factors including host epigenetic factors. This study was designed to identify differentially expressed miRNAs in the primary lymphoid organ, bursae of Fabricius, in response to MD vaccination followed by MDV challenge in two genetically divergent inbred lines of White Leghorns. Small RNA sequencing and bioinformatic analyses of the small RNA sequence reads identified hundreds of miRNAs among all the treatment groups. A small portion of the identified miRNAs was differentially expressed within each of the four treatment groups, which were HVT or CVI988/Rispens vaccinated line 63-resistant birds and line 72-susceptible birds. A direct comparison between the resistant line 63 and susceptible line 72 groups vaccinated with HVT followed by MDV challenge identified five differentially expressed miRNAs. Gene Ontology analysis of the target genes of those five miRNAs revealed that those target genes, in addition to various GO terms, are involved in multiple signaling pathways including MAPK, TGF-β, ErbB, and EGFR1 signaling pathways. The general functions of those pathways reportedly play important roles in oncogenesis, anti-cancer immunity, cancer cell migration, and metastatic progression. Therefore, it is highly likely that those miRNAs may, in part, influence vaccine protection through the pathways.

Integrated analysis of methylation profiles and transcriptome of MDV-infected chicken spleens reveal hypomethylation of CD4 and HMGB1 genes might promote MD tumorigenesis

Marek's disease (MD) is a lymphoproliferative neoplastic disease caused by Marek's disease virus (MDV). Previous studies have showed that DNA methylation was involved in MD development, but systematic studies are still lacking. Herein, we performed whole genome bisulfite sequencing (WGBS) and RNA-seq in MDV-infected tumorous spleens (IN), noninfected spleens (NoIN), and survivor (SUR) spleens of chickens to identify the genes playing important roles in MD tumor transformation. We generated the first genome-wide DNA methylation profile of MDV-infected, noninfected, and survivor chickens. Combined the WGBS and RNA-Seq, we found that the expression of 25% differential expression genes (DEGs) were significantly correlated with methylation of CpG sites in their gene bodies or promoters. Further, we focused on the DEGs with differentially methylated regions (DMRs) on genes' body and promoter, and it showed the expression of 60% DEGs were significantly correlated with methylation of CpG sites in DMRs. Finally, we identified 8 genes, including CD4, CTLA4, DTL, HMGB1, LGMN, NUP210, RAD52, and ZAP70, and their expression was negatively correlated with methylation of DMRs in their promoters in both IN vs. NoIN and IN vs. SUR. These 8 genes showed specifically high expression in IN groups and clustered in module turquoise analyzed by WGCNA. Out of 8 genes, CD4 and HMGB1 were drop in QTLs associated with MD resistance. Thus, we overexpressed the 2 genes to simulate their high expression in the IN group and found they significantly promoted MDCC-MSB-1 cell proliferation, which revealed they might play promoting roles in MD tumorigenesis in IN due to their high expression induced by hypomethylation.

Pre-hatching and post-hatching environmental factors related to epigenetic mechanisms in poultry

Epigenetic modifications are phenotypic changes unrelated to the modification of the DNA sequence. These modifications are essential for regulating cellular differentiation and organism development. In this case, epigenetics controls how the animal's genetic potential is used. The main epigenetic mechanisms are microRNA activity, DNA methylation and histone modification. The literature has repeatedly shown that environmental modulation has a significant influence on the regulation of epigenetic mechanisms in poultry. The aim of this review is to give an overview of the current state of the knowledge in poultry epigenetics in terms of issues relevant to overall poultry production and the improvement of the health status in chickens and other poultry species. One of the main differences between birds and mammals is the stage of embryonic development. The bird's embryo develops outside its mother, so an optimal environment of egg incubation before hatching is crucial for development. It is also the moment when many factors influence the activation of epigenetic mechanisms, i.e., incubation temperature, humidity, light, as well as in ovo treatments. Epigenome of the adult birds, might be modulated by: nutrition, supplementation and treatment, as well as modification of the intestinal microbiota. In addition, the activation of epigenetic mechanisms is influenced by pathogens (i.e., pathogenic bacteria, toxins, viruses and fungi) as well as, the maintenance conditions. Farm animal epigenetics is still a big challenge for scientists. This is a research area with many open questions. Modern methods of epigenetic analysis can serve both in the analysis of biological mechanisms and in the research and applied to production system, poultry health and welfare.

Exploration of Alternative Splicing (AS) Events in MDV-Infected Chicken Spleens

Marek’s disease (MD) was an immunosuppression disease induced by Marek’s disease virus (MDV). MD caused huge economic loss to the global poultry industry, but it also provided an ideal model for studying diseases induced by the oncogenic virus. Alternative splicing (AS) simultaneously produced different isoform transcripts, which are involved in various diseases and individual development. To investigate AS events in MD, RNA-Seq was performed in tumorous spleens (TS), spleens from the survivors (SS) without any lesion after MDV infection, and non-infected chicken spleens (NS). In this study, 32,703 and 25,217 AS events were identified in TS and SS groups with NS group as the control group, and 1198, 1204, and 348 differently expressed (DE) AS events (p-value < 0.05 and FDR < 0.05) were identified in TS vs. NS, TS vs. SS, SS vs. NS, respectively. Additionally, Function enrichment analysis showed that ubiquitin-mediated proteolysis, p53 signaling pathway, and phosphatidylinositol signaling system were significantly enriched (p-value < 0.05). Small structural variations including SNP and indel were analyzed based on RNA-Seq data, and it showed that the TS group possessed more variants on the splice site region than those in SS and NS groups, which might cause more AS events in the TS group. Combined with previous circRNA data, we found that 287 genes could produce both circular and linear RNAs, which suggested these genes were more active in MD lymphoma transformation. This study has expanded the understanding of the MDV infection process and provided new insights for further analysis of resistance/susceptibility mechanisms.

Chicken embryo as a model in epigenetic research

Epigenetics is defined as the study of changes in gene function that are mitotically or meiotically heritable and do not lead to a change in DNA sequence. Epigenetic modifications are important mechanisms that fine tune the expression of genes in response to extracellular signals and environmental changes. In vertebrates, crucial epigenetic reprogramming events occur during early embryogenesis and germ cell development. Chicken embryo, which develops external to the mother's body, can be easily manipulated in vivo and in vitro, and hence, it is an excellent model for performing epigenetic studies. Environmental factors such as temperature can affect the development of an embryo into the phenotype of an adult. A better understanding of the environmental impact on embryo development can be achieved by analyzing the direct effects of epigenetic modifications as well as their molecular background and their intergenerational and transgenerational inheritance. In this overview, the current possibility of epigenetic changes during chicken embryonic development and their effects on long-term postembryonic development are discussed.

Differential DNA Methylation and Gene Expression Between ALV-J-Positive and ALV-J-Negative Chickens

Background: Avian leukosis virus subgroup J (ALV-J) is an oncogenic virus that causes serious economic losses in the poultry industry; unfortunately, there is no effective vaccine against ALV-J. DNA methylation plays a crucial role in several biological processes, and an increasing number of diseases have been proven to be related to alterations in DNA methylation. In this study, we screened ALV-J-positive and -negative chickens. Subsequently, we generated and provided the genome-wide gene expression and DNA methylation profiles by MeDIP-seq and RNA-seq of ALV-J-positive and -negative chicken samples; 8,304 differentially methylated regions (DMRs) were identified by MeDIP-seq analysis (p ≤ 0.005) and 515 differentially expressed genes were identified by RNA-seq analysis (p ≤ 0.05). As a result of an integration analysis, we screened six candidate genes to identify ALV-J-negative chickens that possessed differential methylation in the promoter region. Furthermore, TGFB2 played an important role in tumorigenesis and cancer progression, which suggested TGFB2 may be an indicator for identifying ALV-J infections.

Chicken cecal DNA methylome alteration in the response to Salmonella enterica serovar Enteritidis inoculation

BACKGROUND

Salmonella enterica serovar Enteritidis (SE) is one of the pathogenic bacteria, which affects poultry production and poses a severe threat to public health. Chicken meat and eggs are the main sources of human salmonellosis. DNA methylation is involved in regulatory processes including gene expression, chromatin structure and genomic imprinting. To understand the methylation regulation in the response to SE inoculation in chicken, the genome-wide DNA methylation profile following SE inoculation was analyzed through whole-genome bisulfite sequencing in the current study.

RESULTS

There were 185,362,463 clean reads and 126,098,724 unique reads in the control group, and 180,530,750 clean reads and 126,782,896 unique reads in the inoculated group. The methylation density in the gene body was higher than that in the upstream and downstream regions of the gene. There were 8946 differentially methylated genes (3639 hypo-methylated genes, 5307 hyper-methylated genes) obtained between inoculated and control groups. Methylated genes were mainly enriched in immune-related Gene Ontology (GO) terms and metabolic process terms. Cytokine-cytokine receptor interaction, TGF-beta signaling pathway, FoxO signaling pathway, Wnt signaling pathway and several metabolism-related pathways were significantly enriched. The density of differentially methylated cytosines in miRNAs was the highest. HOX genes were widely methylated.

CONCLUSIONS

The genome-wide DNA methylation profile in the response to SE inoculation in chicken was analyzed. SE inoculation promoted the DNA methylation in the chicken cecum and caused methylation alteration in immune- and metabolic- related genes. Wnt signal pathway, miRNAs and HOX gene family may play crucial roles in the methylation regulation of SE inoculation in chicken. The findings herein will deepen the understanding of epigenetic regulation in the response to SE inoculation in chicken.

Marek's Disease Virus Infection Induced Mitochondria Changes in Chickens

Mitochondria are crucial cellular organelles in eukaryotes and participate in many cell processes including immune response, growth development, and tumorigenesis. Marek’s disease (MD), caused by an avian alpha-herpesvirus Marek’s disease virus (MDV), is characterized with lymphomas and immunosuppression. In this research, we hypothesize that mitochondria may play roles in response to MDV infection. To test it, mitochondrial DNA (mtDNA) abundance and gene expression in immune organs were examined in two well-defined and highly inbred lines of chickens, the MD-susceptible line 7₂ and the MD-resistant line 6₃. We found that mitochondrial DNA contents decreased significantly at the transformation phase in spleen of the MD-susceptible line 7₂ birds in contrast to the MD-resistant line 6₃. The mtDNA-genes and the nucleus-genes relevant to mtDNA maintenance and transcription, however, were significantly up-regulated. Interestingly, we found that POLG2 might play a potential role that led to the imbalance of mtDNA copy number and gene expression alteration. MDV infection induced imbalance of mitochondrial contents and gene expression, demonstrating the indispensability of mitochondria in virus-induced cell transformation and subsequent lymphoma formation, such as MD development in chicken. This is the first report on relationship between virus infection and mitochondria in chicken, which provides important insights into the understanding on pathogenesis and tumorigenesis due to viral infection.

Identification of DNA methylation regulated novel host genes relevant to inhibition of virus replication in porcine PK15 cell using double stranded RNA mimics and DNA methyltransferase inhibitor

During RNA viruses's replication, double-stranded RNA (dsRNA) is normally produced and induce host innate immune response. Most of gene activation due cytokine mediated but which are due to methylation mediated is still unknown. In the study, DNA methylome was integrated with our previous transcriptome data to investigate the differentially methylated regions and genes using MeDIP-chip technology. We found that the transcriptional expressions of 15, 37 and 18 genes were negatively related with their promoter DNA methylation levels in the cells treated by PolyI:C, Aza-CdR, as well as PolyI:C plus Aza-CdR, respectively, compared with the untreated cells. GO analysis revealed hypo-methylated genes (BNIP3L and CDK9) and a hyper-methylated gene (ZC3HAV1) involved in the host response to viral replication. Our results suggest that these novel genes targeted by DNA methylation can be potential markers relevant to virus replication and host innate immune response to set up a medical model of infectious diseases.

Gga-miR-130b-3p inhibits MSB1 cell proliferation, migration, invasion, and its downregulation in MD tumor is attributed to hypermethylation

Marek's disease is an oncogenic and lymphoproliferative disease of chickens caused by Marek's disease virus. Hypermethylation or hypomethylation of CpG islands in gene promoter region are involved in the initiation and progression of carcinogenesis. In this study, we analyzed differential methylation levels of upstream region of gga-miR-130b-3p gene between Marek's disease virus-infected tumorous and non-infected spleens. Around the upstream 1 kb of gga-miR-130b-3p gene, two amplicons were designed that covered 616 bp. There were forty-eight CpG sites in this region. CpG sites in this region presented higher methylation level in tumorous spleens compared with that in non-infected ones. There were eight CpG sites significantly hypermethylated in tumorous spleens. The expression level of three DNA methyltransferases including DNMT1, DNMT3a and DNMT3b increased and the expression level of Tet ten-eleven translocation protein 2 remarkably decreased in tumorous spleens. Hypermethylation in the upstream region of gga-miR-130b-3p gene might be a direct reason for its downregulation in MD tumorous tissues. Moreover, cell proliferation of Marek's disease lymphoblastoid cell line MDCC-MSB1 was remarkably inhibited at 24, 36, 48, 60 and 72 h post-gga-miR-130b-3p-agomir transfection. The transwell migration assay indicated cell number of migration was significantly lower in miRNA agomir transfection group. Matrix metalloproteinases MMP2 and MMP9 are involved in tumor invasion, and their protein levels were significantly downregulated at 72 h post-miRNA-agomir transfection. Collectively, these results indicated that hypermethylation in upstream region of gga-miR-130b-3p gene contributed to its downregulation in tumorous tissues. Gga-miR-130b-3p plays an inhibitory role in lymphomatous cell transformation.

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.

Genome-wide mapping of DNase I hypersensitive sites and association analysis with gene expression in MSB1 cells

DNase I hypersensitive sites (DHSs) mark diverse classes of cis-regulatory regions, such as promoters and enhancers. MSB-1 derived from chicken Marek's disease (MD) lymphomas is an MDV-transformed CD4+ T-cell line for MD study. Previously, DNase I HS sites were studied mainly in human cell types for mammalian. To capture the regulatory elements specific to MSB1 cells and explore the molecular mechanisms of T-cell transformation caused by MDV in MD, we generated high-quality of DHSs map and gene expression profile for functional analysis in MSB1 cell line. The total of 21,724 significant peaks of DHSs was identified from around 40 million short reads. DHSs distribution varied between chromosomes and they preferred to enrich in the gene-rich chromosomes. More interesting, DHSs enrichments appeared to be scarce on regions abundant in CpG islands. Besides, we integrated DHSs into the gene expression data and found that DHSs tended to enrich on high expressed genes throughout whole gene regions while DHSs did not show significant changes for low and silent expressed genes. Furthermore, the correlation of DHSs with lincRNAs expression was also calculated and it implied that enhancer-associated lincRNAs probably originated from enhancer-like regions of DHSs. Together, our results indicated that DNase I HS sites highly correlate with active genes expression in MSB1 cells, suggesting DHSs can be considered as markers to identify the cis-regulatory elements associated with chicken Marek's disease.

DNA methylome in spleen of avian pathogenic Escherichia coli-challenged broilers and integration with mRNA expression

Avian pathogenic Escherichia coli (APEC) are responsible for heavy economic losses in poultry industry. Here we investigate DNA methylome of spleen and identify functional DNA methylation changes related to host response to APEC among groups of non-challenged chickens (NC), challenged with mild (MD) and severe pathology (SV). DNA methylation was enriched in the gene bodies and repeats. Promoter and CGIs are hypomethylated. Integration analysis revealed 22, 87, and 9 genes exhibiting inversely changed DNA methylation and gene expression in NC vs. MD, NC vs. SV, and MD vs. SV, respectively. IL8, IL2RB, and IL1RAPL1 were included. Gene network analysis suggested that besides inflammatory response, other networks and pathways such as organismal injury and abnormalities, cell signaling and molecular transport, are probably related to host response to APEC infection. Moreover, methylation changes in cell cycle processes might contribute to the lesion phenotype differences between MD and SV.

Genome-wide effects of DNA methyltransferase inhibitor on gene expression in double-stranded RNA transfected porcine PK15 cells

Double-stranded RNA (dsRNA) is produced in host cells during viral replication. The effects of DNA demethylation on gene expression in dsRNA transfected swine cells are unclear. The study aims to profile the transcriptome changes which are induced by DNA methyltransferase inhibitor (Aza-CdR) in porcine PK15 cells transfected with viral-like dsRNA (Poly(I:C)). A total of 44, 76 and 952 differentially expressed genes (DEGs) were detected in the cells treated by Poly(I:C) plus Aza-CdR (P+A), Poly(I:C) (P) or Aza-CdR (A) alone compared to the controls (C). Immune response-related pathways are observed in the comparison of A vs. C and P vs. C, and the genes in the pathways were recovered in the comparison of (P+A) vs. C. GO analysis indicated that Aza-CdR has negative regulatory effects on viral reproduction. The results suggest that the stimulant of Poly(I:C) could be regressed by Aza-CdR. These observations provide new insights into the epigenetic regulatory effects on viral replication.

Mosaicism of HTLV-1 5' LTR CpG methylation in the absence of malignancy

Viral expression varies widely between untransformed HTLV-1 positive clones derived from infected individuals without malignancy. Here we show that, in the absence of malignancy, 68% of HTLV-1 positive clones carry deleted (10%) or methylated (58%) 5' LTR. These changes were found to contribute to the fluctuation of viral expression between clones. 5' LTR deletions strongly impaired tax expression and thereby clonal expansion, telomere homeostasis, and genetic instability. The effects of CpG methylation on viral transcription were qualitatively and quantitatively different from those of LTR deletions and preserved the preleukemic features of HTLV-1(+)CD4(+) clones. 5' LTR methylation varied not only between clones but also between cells belonging to the same clones, between distinct integrated sequences within the same cells, and between CpG dyads along the 5' LTR. Such distributions suggest that a dynamic methylation spectrum might help sustain persistent infection.

Genomic instability in liver cells caused by an LPS-induced bystander-like effect

Bacterial infection has been linked to carcinogenesis, however, there is lack of knowledge of molecular mechanisms that associate infection with the development of cancer. We analyzed possible effects of the consumption of heat-killed E. coli O157:H7 cells or its cellular components, DNA, RNA, protein or lipopolysaccharides (LPS) on gene expression in naïve liver cells. Four week old mice were provided water supplemented with whole heat-killed bacteria or bacterial components for a two week period. One group of animals was sacrificed immediately, whereas another group was allowed to consume uncontaminated tap water for an additional two weeks, and liver samples were collected, post mortem. Liver cells responded to exposure of whole heat-killed bacteria and LPS with alteration in γH2AX levels and levels of proteins involved in proliferation, DNA methylation (MeCP2, DNMT1, DNMT3A and 3B) or DNA repair (APE1 and KU70) as well as with changes in the expression of genes involved in stress response, cell cycle control and bile acid biosynthesis. Other bacterial components analysed in this study did not lead to any significant changes in the tested molecular parameters. This study suggests that lipopolysaccharides are a major component of Gram-negative bacteria that induce molecular changes within naïve cells of the host.

Epigenetics and phenotypic variability: some interesting insights from birds

Little is known about epigenetic mechanisms in birds with the exception of the phenomenon of dosage compensation of sex chromosomes, although such mechanisms could be involved in the phenotypic variability of birds, as in several livestock species. This paper reviews the literature on epigenetic mechanisms that could contribute significantly to trait variability in birds, and compares the results to the existing knowledge of epigenetic mechanisms in mammals. The main issues addressed in this paper are: (1) Does genomic imprinting exist in birds? (2) How does the embryonic environment influence the adult phenotype in avian species? (3) Does the embryonic environment have an impact on phenotypic variability across several successive generations? The potential for epigenetic studies to improve the performance of individual animals through the implementation of limited changes in breeding conditions or the addition of new parameters in selection models is still an open question.

DNMT gene expression and methylome in Marek's disease resistant and susceptible chickens prior to and following infection by MDV

Marek's disease (MD) is characterized as a T cell lymphoma induced by a cell-associated α-herpesvirus, Marek's disease virus type 1 (MDV1). As with many viral infectious diseases, DNA methylation variations were observed in the progression of MD; these variations are thought to play an important role in host-virus interactions. We observed that DNA methyltransferase 3a (DNMT3a) and 3b (DNMT3b) were differentially expressed in chicken MD-resistant line 6 3 and MD-susceptible line 7 2 at 21 d after MDV infection. To better understand the role of methylation variation induced by MDV infection in both chicken lines, we mapped the genome-wide DNA methylation profiles in each line using Methyl-MAPS (methylation mapping analysis by paired-end sequencing). Collectively, the data sets collected in this study provide a more comprehensive picture of the chicken methylome. Overall, methylation levels were reduced in chickens from the resistant line 6 3 after MDV infection. We identified 11,512 infection-induced differential methylation regions (iDMRs). The number of iDMRs was larger in line 7 2 than in line 6 3, and most of iDMRs found in line 6 3 were overlapped with the iDMRs found in line 7 2. We further showed that in vitro methylation levels were associated with MDV replication, and found that MDV propagation in the infected cells was restricted by pharmacological inhibition of DNA methylation. Our results suggest that DNA methylation in the host may be associated with disease resistance or susceptibility. The methylation variations induced by viral infection may consequentially change the host transcriptome and result in diverse disease outcomes.