Histone methylation analysis and pathway predictions in chickens after MDV infection

Transcriptomic analyses of host-virus interactions during in vitro infection with wild-type and glycoprotein g-deficient (ΔgG) strains of ILTV in primary and continuous cell cultures

Infectious laryngotracheitis (ILT) remains a significant concern for the poultry industry worldwide due to its impact on animal welfare and its substantial economic consequences. The disease is caused by the alphaherpesvirus, infectious laryngotracheitis virus (ILTV). This study investigated in vitro host-virus interactions of a glycoprotein G (gG) deletion mutant vaccine strain of ILTV (ΔgG ILTV), and its parent wild-type strain (CSW-1 ILTV). Inoculations were performed separately for the two strains of ILTV using both a primary (chicken embryonic kidney, CEK) and a continuous culture (leghorn male hepatoma, LMH) of chicken cells. Transcriptome analysis was performed at 12 hours post infection. Each cell-type displayed distinct effects on host and viral gene transcription, with a greater number of viral and host genes differentially transcribed in CEK cells and LMH cells, respectively. Both cell-types infected with either strain demonstrated enrichment of pathways related to signalling, and gene ontologies (GO) associated with chemotaxis. Infection with either strain upregulated both SOCS proteins and certain proto-oncogenes, which may contribute to prolonged viral persistence by promoting immunosuppression and preventing apoptosis, respectively. Patterns of gene transcription related to cytokines, chemokines, endosomal TLRs, and interferon responses, as well as pathways associated with histone acetylation, transport, and extracellular matrix organization were similar within each cell type, regardless of the viral strain. In CEK cells, GO terms and pathways were downregulated uniquely after CSW-1 ILTV infection, indicating a viral-strain specific effect in this cell-type. Overall, this study highlights that the observed differences in host and ILTV gene transcription in vitro were more strongly influenced by the cell-types used rather than the presence or absence of gG. This underscores the importance of cell-line selection in studying host-virus interactions and interpreting experimental results.

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.

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.

Non-additive effect of the DNA methylation inhibitor, 5-Aza-dC, and glass as a culture surface on osteogenic differentiation

The clinical need for bone regenerative solutions is expanding with increasing life expectancy and escalating incidence of accidents. Several strategies are being investigated to enhance the osteogenic differentiation of stem cells. We previously reported two different approaches for this purpose, in monolayer and three-dimensional cell culture. The first approach was based on pretreating cells with 5-Aza-dC, a DNA methylation inhibitor, before the applying the differentiation media. The second approach was based on culturing cells on a glass surface during differentiation. In this study, we investigated the potential effect of combining both methods. Our results suggested that both approaches were associated with decreasing global DNA methylation levels. Cells cultured as a monolayer on glass surface showed enhancement in alkaline phosphatase activity at day 10, while 5-Aza-dC pretreatment enhanced the activity at day 5, irrespective of the culture surface. In three-dimensional pellet culture, 5-Aza-dC pretreatment enhanced osteogenesis through Runx-2 and TGF-β1 upregulation while the glass surface induced Osterix. Furthermore, pellets cultured on glass showed upregulation of a group of miRNAs, including pro-osteogenesis miR- 20a and miR -148b and anti-osteogenesis miR -125b, miR -31, miR -138, and miR -133a. Interestingly, 5-Aza-dC was not associated with a change of miRNAs in cells cultured on tissue culture plastic but reverted the upregulated miRNAs on the glass to the basal level. This study confirms the two approaches for enhancing osteogenic differentiation and contradicts their combination.

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.

Effects of acute heat stress on protein expression and histone modification in the adrenal gland of male layer-type country chickens

The adrenal gland responds to heat stress by epinephrine and glucocorticoid release to alleviate the adverse effects. This study investigated the effect of acute heat stress on the protein profile and histone modification in the adrenal gland of layer-type country chickens. A total of 192 roosters were subject to acute heat stress and thereafter classified into a resistant or susceptible group according to body temperature change. The iTRAQ analysis identified 80 differentially expressed proteins, in which the resistant group had a higher level of somatostatin and hydroxy-δ-5-steroid dehydrogenase but a lower parathymosin expression in accordance with the change of serum glucocorticoid levels. Histone modification analysis identified 115 histone markers. The susceptible group had a higher level of tri-methylation of histone H3 lysine 27 (H3K27me3) and showed a positive crosstalk with K36me and K37me in the H3 tails. The differential changes of body temperature projected in physiological regulation at the hypothalamus-pituitary-adrenal axis suggest the genetic heterogeneity in basic metabolic rate and efficiency for heat dissipation to acclimate to thermal stress and maintain body temperature homeostasis. The alteration of adrenal H3K27me3 level was associated with the endocrine function of adrenal gland and may contribute to the thermotolerance of chickens.

Mapping QTL Associated with Resistance to Avian Oncogenic Marek's Disease Virus (MDV) Reveals Major Candidate Genes and Variants

Marek’s disease (MD) represents a significant global economic and animal welfare issue. Marek’s disease virus (MDV) is a highly contagious oncogenic and highly immune-suppressive α-herpes virus, which infects chickens, causing neurological effects and tumour formation. Though partially controlled by vaccination, MD continues to have a profound impact on animal health and on the poultry industry. Genetic selection provides an alternative and complementary method to vaccination. However, even after years of study, the genetic mechanisms underlying resistance to MDV remain poorly understood. The Major Histocompatability Complex (MHC) is known to play a role in disease resistance, along with a handful of other non-MHC genes. In this study, one of the largest to date, we used a multi-facetted approach to identify quantitative trait locus regions (QTLR) influencing resistance to MDV, including an F₆ population from a full-sib advanced intercross line (FSIL) between two elite commercial layer lines differing in resistance to MDV, RNA-seq information from virus challenged chicks, and genome wide association study (GWAS) from multiple commercial lines. Candidate genomic elements residing in the QTLR were further tested for association with offspring mortality in the face of MDV challenge in eight pure lines of elite egg-layer birds. Thirty-eight QTLR were found on 19 chicken chromosomes. Candidate genes, microRNAs, long non-coding RNAs and potentially functional mutations were identified in these regions. Association tests were carried out in 26 of the QTLR, using eight pure lines of elite egg-layer birds. Numerous candidate genomic elements were strongly associated with MD resistance. Genomic regions significantly associated with resistance to MDV were mapped and candidate genes identified. Various QTLR elements were shown to have a strong genetic association with resistance. These results provide a large number of significant targets for mitigating the effects of MDV infection on both poultry health and the economy, whether by means of selective breeding, improved vaccine design, or gene-editing technologies.

Adiponectin and its receptor genes' expression in response to Marek's disease virus infection of White Leghorns

Marek’s disease virus (MDV) causes T-cell lymphoma in susceptible chicken and is also related to an imbalance of the lipid metabolism. Adiponectin is a circulatory cytokine secreted from adipose tissue and exerts critical metabolic functions. Although the associations between adiponectin and diseases, including lipid disorder and noncardiac vascular diseases, have been reported, little is known about the relationship between MDV infection and adiponectin. Here, we challenged white Leghorns from Marek’s disease (MD)-susceptible and MD-resistant lines with MDV at 7 D of age and then explored the body weight and plasma lipoprotein levels at 21 D after MDV infection. Meanwhile, adiponectin and the expression of its receptors were detected using quantitative real-time PCR and Western blot. The results showed that MDV infection induced body weight loss in all the experimental birds. Meanwhile, the concentrations of total cholesterol and high-density lipoprotein were lower after the infection, although there was no significant difference (P > 0.05). However, the infection did not affect adiponectin circulating levels in plasma. MD-susceptible birds had much lower plasma adiponectin than MD-resistant birds (P < 0.01). In abdominal fat, there was no significant difference in adiponectin mRNA level. Still, we observed a significant decrease in adiponectin protein concentration, as well as adipoR1 and adipoR2, at both mRNA and protein levels in the infected compared with the noninfected MD-susceptible chickens. In the spleen, MDV infection significantly reduced the adiponectin mRNA expression but increased the protein in MD-susceptible chickens, which decreased both adipoR1 mRNA expression and protein levels. Also interestingly, the adipoR1 mRNA expression level was significantly increased in MD-susceptible chickens in the liver after MDV infection. All findings in the present study provided interesting insights into adiponectin metabolism in chickens after MDV infection, which helps to advance the understanding of lipid metabolism in response to herpesvirus infection.

Genetic assessment of inbred chicken lines indicates genomic signatures of resistance to Marek's disease

Background

Marek’s disease (MD) is a highly contagious pathogenic and oncogenic disease primarily affecting chickens. However, the mechanisms of genetic resistance for MD are complex and not fully understood. MD-resistant line 6₃ and MD-susceptible line 7₂ are two highly inbred progenitor lines of White Leghorn. Recombinant Congenic Strains (RCS) were developed from these two lines, which show varied susceptibility to MD.

Results

We investigated genetic structure and genomic signatures across the genome, including the line 6₃ and line 7₂, six RCSs, and two reciprocally crossed flocks between the lines 6₃ and 7₂ (F1 6₃ × 7₂ and F1 7₂ × 6₃) using Affymetrix® Axiom® HD 600 K genotyping array. We observed 18 chickens from RCS lines were specifically clustered into resistance sub-groups distributed around line 6₃. Additionally, homozygosity analysis was employed to explore potential genetic components related to MD resistance, while runs of homozygosity (ROH) are regions of the genome where the identical haplotypes are inherited from each parent. We found several genes including SIK, SOX1, LIG4, SIK1 and TNFSF13B were contained in ROH region identified in resistant group (line 6₃ and RCS), and these genes have been reported that are contribute to immunology and survival. Based on F_ST based population differential analysis, we also identified important genes related to cell death and anti-apoptosis, including AKT1, API5, CDH13, CFDP and USP15, which could be involved in divergent selection during inbreeding process.

Conclusions

Our findings offer valuable insights for understanding the genetic mechanism of resistance to MD and the identified genes could be considered as candidate biomarkers in further evaluation.

Electronic supplementary material

The online version of this article (10.1186/s40104-018-0281-x) contains supplementary material, which is available to authorized users.

A demonstration of the H3 trimethylation ChIP-seq analysis of galline follicular mesenchymal cells and male germ cells

OBJECTIVE

Trimethylation of histone 3 (H3) at 4th lysine N-termini (H3K4me3) in gene promoter region was the universal marker of active genes specific to cell lineage. On the contrary, coexistence of trimethylation at 27th lysine (H3K27me3) in the same loci-the bivalent H3K4m3/H3K27me3 was known to suspend the gene transcription in germ cells, and could also be inherited to the developed stem cell. In galline species, throughout example of H3K4m3 and H3K27me3 ChIP-seq analysis was still not provided. We therefore designed and demonstrated such procedures using ChIP-seq and mRNA-seq data of chicken follicular mesenchymal cells and male germ cells.

METHODS

Analytical workflow was designed and provided in this study. ChIP-seq and RNA-seq datasets of follicular mesenchymal cells and male germ cells were acquired and properly preprocessed. Peak calling by Model-based analysis of ChIP-seq 2 was performed to identify H3K4m3 or H3K27me3 enriched regions (Fold-change≥2, FDR≤0.01) in gene promoter regions. Integrative genomics viewer was utilized for cellular retinoic acid binding protein 1 (CRABP1), growth differentiation factor 10 (GDF10), and gremlin 1 (GREM1) gene explorations.

RESULTS

The acquired results indicated that follicular mesenchymal cells and germ cells shared several unique gene promoter regions enriched with H3K4me3 (5,704 peaks) and also unique regions of bivalent H3K4m3/H3K27me3 shared between all cell types and germ cells (1,909 peaks). Subsequent observation of follicular mesenchyme-specific genes-CRABP1, GDF10, and GREM1 correctly revealed vigorous transcriptions of these genes in follicular mesenchymal cells. As expected, bivalent H3K4m3/H3K27me3 pattern was manifested in gene promoter regions of germ cells, and thus suspended their transcriptions.

CONCLUSION

According the results, an example of chicken H3K4m3/H3K27me3 ChIP-seq data analysis was successfully demonstrated in this study. Hopefully, the provided methodology should hereby be useful for galline ChIP-seq data analysis in the future.

An Assessment of Fixed and Native Chromatin Preparation Methods to Study Histone Post-Translational Modifications at a Whole Genome Scale in Skeletal Muscle Tissue

Background

Genomic loci associated with histone marks are typically analyzed by immunoprecipitation of the chromatin followed by quantitative-PCR (ChIP-qPCR) or high throughput sequencing (ChIP-seq). Chromatin can be either cross-linked (X-ChIP) or used in the native state (N-ChIP). Cross-linking of DNA and proteins helps stabilizing their interactions before analysis. Despite X-ChIP is the most commonly used method, muscle tissue fixation is known to be relatively inefficient. Moreover, no protocol described a simple and reliable preparation of skeletal muscle chromatin of sufficient quality for subsequent high-throughput sequencing. Here we aimed to set-up and compare both chromatin preparation methods for a genome-wide analysis of H3K27me3, a broad-peak histone mark, using chicken P. major muscle tissue.

Results

Fixed and unfixed chromatin were prepared from chicken muscle tissues (Pectoralis major). Chromatin fixation, shearing by sonication or digestion and immunoprecipitation performed equivalently. High-quality Illumina reads were obtained (q30 > 93%). The bioinformatic analysis of the data was performed using epic, a tool based on SICER, and MACS2. Forty millions of reads were analyzed for both X-ChIP-seq and N-ChIP-seq experiments. Surprisingly, H3K27me3 X-ChIP-seq analysis led to the identification of only 2000 enriched regions compared to about 15,000 regions identified in the case of N-ChIP-seq. N-ChIP-seq peaks were more consistent between replicates compared to X-ChIP-seq. Higher N-ChIP-seq enrichments were confirmed by ChIP-qPCR at the PAX5 and SOX2 loci known to be enriched for H3K27me3 in myotubes and at the loci of common regions of enrichment identified in this study.

Conclusions

Our findings suggest that the preparation of muscle chromatin for ChIP-seq in cross-linked conditions can compromise the systematic analysis of broad histone marks. Therefore, native chromatin preparation should be preferred to cross-linking when a ChIP experiment has to be performed on skeletal muscle tissue, particularly when a broad source signal is considered.

Electronic supplementary material

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

Characterization of Copy Number Variation's Potential Role in Marek's Disease

Marek’s Disease (MD) is a highly contagious pathogenic and oncogenic disease primarily affecting chickens. Chicken Lines 63 and 72, as well as their recombinant congenic strains (RCS) with varied susceptibility to MD, are ideal models to study the complex mechanisms of genetic resistance to MD. In this study, we investigated copy number variation (CNV) in these inbred chicken lines using the Affymetrix Axiom HD 600 K SNP genotyping array. We detected 393 CNV segments across all ten chicken lines, of which 12 CNVs were specifically identified in Line 72. We then assessed genetic structure based on CNV and observed markedly different patterns. Finally, we validated two deletion events in Line 72 and correlated them with genes expression using qPCR and RNA-seq, respectively. Our combined results indicated that these two CNV deletions were likely to contribute to MD susceptibility.

Genome-Wide H3K4me3 Analysis in Angus Cattle with Divergent Tenderness

Tenderness is one of the most important properties of meat quality, which is influenced by genetic and environmental factors. As an intensively studied epigenetic marker, histone methylation, occurring on arginine and lysine residues, has pivotal regulatory functions on gene expression. To examine whether histone methylation involves in beef tenderness variation, we analyzed the transcriptome and H3K4me3 enrichment profiles of muscle strips obtained from the longissimus dorsi (LD) of Angus steers previously classify to the tender or tough group. We first plotted a global bovine H3K4me3 map on chromosomes and called peak-enriched regions and genes. We found that majorities of H3K4me3 on genes were occupying the first intron and intergenic regions and its maps displayed similar patterns in tender and tough groups, with high H3K4me3 enrichment surrounding the transcription start site (TSS). We also explored the relationship of H3K4me3 and gene expression. The results showed that H3K4me3 enrichment is highly positively correlated with gene expression across the whole genome. Cluster analysis results confirmed the relationship of H3K4me3 enrichment and gene expression. By using a pathway-based approach in genes with H3K4me3 enrichment in promoter regions from the tender cluster, we revealed that those genes involved in the development of different tissues-connective tissue, skeletal and muscular system and functional tissues-; while in tough group those genes engaged in cell death, lipid metabolism and small molecule biochemistry. The results from this study provide a deep insight into understanding of the mechanisms of epigenetic regulations in meat quality and beef tenderness.

Histone modifications induced by MDV infection at early cytolytic and latency phases

Background

Marek’s disease (MD) is a highly contagious, lymphomatous disease of chickens induced by a herpesvirus, Marek’s disease virus (MDV) that is the cause of major annual losses to the poultry industry. MD pathogenesis involves multiple stages including an early cytolytic phase and latency, and transitions between these stages are governed by several host and environmental factors. The success of vaccination strategies has led to the increased virulence of MDV and selective breeding of naturally resistant chickens is seen as a viable alternative. While multiple gene expression studies have been performed in resistant and susceptible populations, little is known about the epigenetic effects of infection.

Results

In this study, we investigated temporal chromatin signatures induced by MDV by analyzing early cytolytic and latent phases of infection in the bursa of Fabricius of MD-resistant and –susceptible birds. Major global variations in chromatin marks were observed at different stages of MD in the two lines. Differential H3K27me3 marks were associated with immune-related pathways, such as MAP kinase signaling, focal adhesion and neuroactive ligand receptor interaction, and suggested varying degrees of silencing in response to infection. Immune-related microRNAs, e.g. gga-miR-155 and gga-miR-10b, bore chromatin signatures, which suggested their contribution to MD-susceptibility. Finally, several members of the focal adhesion pathway, e.g. THBS4 and ITGA1, showed marked concordance between gene expression and chromatin marks indicating putative epigenetic regulation in response to MDV infection.

Conclusion

Our comprehensive analysis of chromatin signatures, therefore, revealed further clues about the epigenetic effects of MDV infection although further studies are necessary to elucidate the functional implications of the observed variations in histone modifications.

Electronic supplementary material

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

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.

Marek's disease virus infection induces widespread differential chromatin marks in inbred chicken lines

Background

Marek’s disease (MD) is a neoplastic disease in chickens caused by the MD virus (MDV). Successful vaccine development against MD has resulted in increased virulence of MDV and the understanding of genetic resistance to the disease is, therefore, crucial to long-term control strategies. Also, epigenetic factors are believed to be one of the major determinants of disease response.

Results

Here, we carried out comprehensive analyses of the epigenetic landscape induced by MDV, utilizing genome-wide histone H3 lysine 4 and lysine 27 trimethylation maps from chicken lines with varying resistance to MD. Differential chromatin marks were observed on genes previously implicated in the disease such as MX1 and CTLA-4 and also on genes reported in other cancers including IGF2BP1 and GAL. We detected bivalent domains on immune-related transcriptional regulators BCL6, CITED2 and EGR1, which underwent dynamic changes in both lines as a result of MDV infection. In addition, putative roles for GAL in the mechanism of MD progression were revealed.

Conclusion

Our results confirm the presence of widespread epigenetic differences induced by MD in chicken lines with different levels of genetic resistance. A majority of observed epigenetic changes were indicative of increased levels of viral infection in the susceptible line symptomatic of lowered immunocompetence in these birds caused by early cytolytic infection. The GAL system that has known anti-proliferative effects in other cancers is also revealed to be potentially involved in MD progression. Our study provides further insight into the mechanisms of MD progression while revealing a complex landscape of epigenetic regulatory mechanisms that varies depending on host factors.