Yanghua He

The conservation and signatures of lincRNAs in Marek's disease of chicken.

Long intergenic non-coding RNAs (lincRNAs) associated with a number of cancers and other diseases have been identified in mammals, but they are still formidable to be comprehensively identified and characterized. Marek’s disease (MD) is a T cell lymphoma of chickens induced by Marek’s disease virus (MDV). Here, we used a MD chicken model to develop a precise pipeline for identifying lincRNAs and to determine the roles of lincRNAs in T cell tumorigenesis. More than 1,000 lincRNA loci were identified in chicken bursa. Computational analyses demonstrated that lincRNAs are conserved among different species such as human, mouse and chicken. The putative lincRNAs were found to be associated with a wide range of biological functions including immune responses. Interestingly, we observed distinct lincRNA expression signatures in bursa between MD resistant and susceptible lines of chickens. One of the candidate lincRNAs, termed linc-satb1, was found to play a crucial role in MD immune response by regulating a nearby protein-coding gene SATB1. Thus, our results manifested that lincRNAs may exert considerable influence on MDV-induced T cell tumorigenesis and provide a rich resource for hypothesis-driven functional studies to reveal genetic mechanisms underlying susceptibility to tumorigenesis.

Integrated metabolomic and transcriptome analyses reveal finishing forage affects metabolic pathways related to beef quality and animal welfare.

Beef represents a major dietary component and source of protein in many countries. With an increasing demand for beef, the industry is currently undergoing changes towards naturally produced beef. However, the true differences between the feeding systems, especially the biochemical and nutritional aspects, are still unclear. Using transcriptome and metabolome profiles, we identified biological pathways related to the differences between grass- and grain-fed Angus steers. In the latissimus dorsi muscle, we have recognized 241 differentially expressed genes (FDR < 0.1). The metabolome examinations of muscle and blood revealed 163 and 179 altered compounds in each tissue (P < 0.05), respectively. Accordingly, alterations in glucose metabolism, divergences in free fatty acids and carnitine conjugated lipid levels, and altered β-oxidation have been observed. The anti-inflammatory n3 polyunsaturated fatty acids are enriched in grass finished beef, while higher levels of n6 PUFAs in grain finished animals may promote inflammation and oxidative stress. Furthermore, grass-fed animals produce tender beef with lower total fat and a higher omega3/omega6 ratio than grain-fed ones, which could potentially benefit consumer health. Most importantly, blood cortisol levels strongly indicate that grass-fed animals may experience less stress than the grain-fed individuals. These results will provide deeper insights into the merits and mechanisms of muscle development.

Ruminal Transcriptomic Analysis of Grass-Fed and Grain-Fed Angus Beef Cattle

Beef represents a major diet component and one of the major sources of protein in human. The beef industry in the United States is currently undergoing changes and is facing increased demands especially for natural grass-fed beef. The grass-fed beef obtained their nutrients directly from pastures, which contained limited assimilable energy but abundant amount of fiber. On the contrary, the grain-fed steers received a grain-based regime that served as an efficient source of high-digestible energy. Lately, ruminant animals have been accused to be a substantial contributor for the green house effect. Therefore, the concerns from environmentalism, animal welfare and public health have driven consumers to choose grass-fed beef. Rumen is one of the key workshops to digest forage constituting a critical step to supply enough nutrients for animals’ growth and production. We hypothesize that rumen may function differently in grass- and grain-fed regimes. The objective of this study was to find the differentially expressed genes in the ruminal wall of grass-fed and grain-fed steers, and then explore the potential biopathways. In this study, the RNA Sequencing (RNA-Seq) method was used to measure the gene expression level in the ruminal wall. The total number of reads per sample ranged from 24,697,373 to 36,714,704. The analysis detected 342 differentially expressed genes between ruminal wall samples of animals raised under different regimens. The Fisher’s exact test performed in the Ingenuity Pathway Analysis (IPA) software found 16 significant molecular networks. Additionally, 13 significantly enriched pathways were identified, most of which were related to cell development and biosynthesis. Our analysis demonstrated that most of the pathways enriched with the differentially expressed genes were related to cell development and biosynthesis. Our results provided valuable insights into the molecular mechanisms resulting in the phenotype difference between grass-fed and grain-fed cattle.

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 Mareks 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 Mareks disease.

Methylome Analysis in Chickens Immunized with Infectious Laryngotracheitis Vaccine

In this study we investigated the methylome of chickens immunized with Infectious laryngotracheitis (ILT) vaccine derived from chicken embryos. Methyl-CpG binding domain protein-enriched genome sequencing (MBD-Seq) method was employed in the detection of the 1,155 differentially methylated regions (DMRs) across the entire genome. After validation, we ascertained the genomic DMRs distribution and annotated them regarding genes, transcription start sites (TSS) and CpG islands. We found that global DNA methylation decreased in vaccinated birds, presenting 704 hypomethylated and 451 hypermethylated DMRs, respectively. Additionally, we performed an enrichment analysis detecting gene networks, in which cancer and RNA post-transcriptional modification appeared in the first place, followed by humoral immune response, immunological disease and inflammatory disease. The top four identified canonical pathways were EIF2 signaling, regulation of EIF4 and p70S6K signaling, axonal guidance signaling and mTOR signaling, providing new insight regarding the mechanisms of ILT etiology. Lastly, the association between DNA methylation and differentially expressed genes was examined, and detected negative correlation in seventeen of the eighteen genes.

Combined analysis of DNA methylome and transcriptome reveal novel candidate genes with susceptibility to bovine Staphylococcus aureus subclinical mastitis

Subclinical mastitis is a widely spread disease of lactating cows. Its major pathogen is Staphylococcus aureus (S. aureus). In this study, we performed genome-wide integrative analysis of DNA methylation and transcriptional expression to identify candidate genes and pathways relevant to bovine S. aureus subclinical mastitis. The genome-scale DNA methylation profiles of peripheral blood lymphocytes in cows with S. aureus subclinical mastitis (SA group) and healthy controls (CK) were generated by methylated DNA immunoprecipitation combined with microarrays. We identified 1078 differentially methylated genes in SA cows compared with the controls. By integrating DNA methylation and transcriptome data, 58 differentially methylated genes were shared with differently expressed genes, in which 20.7% distinctly hypermethylated genes showed down-regulated expression in SA versus CK, whereas 14.3% dramatically hypomethylated genes showed up-regulated expression. Integrated pathway analysis suggested that these genes were related to inflammation, ErbB signalling pathway and mismatch repair. Further functional analysis revealed that three genes, NRG1, MST1 and NAT9, were strongly correlated with the progression of S. aureus subclinical mastitis and could be used as powerful biomarkers for the improvement of bovine mastitis resistance. Our studies lay the groundwork for epigenetic modification and mechanistic studies on susceptibility of bovine mastitis.

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.

Corrigendum: The conservation and signatures of lincRNAs in Marek's disease of chicken.

Long intergenic non-coding RNAs (lincRNAs) associated with a number of cancers and other diseases have been identified in mammals, but they are still formidable to be comprehensively identified and characterized. Marek’s disease (MD) is a T cell lymphoma of chickens induced by Marek’s disease virus (MDV). Here, we used a MD chicken model to develop a precise pipeline for identifying lincRNAs and to determine the roles of lincRNAs in T cell tumorigenesis. More than 1,000 lincRNA loci were identified in chicken bursa. Computational analyses demonstrated that lincRNAs are conserved among different species such as human, mouse and chicken. The putative lincRNAs were found to be associated with a wide range of biological functions including immune responses. Interestingly, we observed distinct lincRNA expression signatures in bursa between MD resistant and susceptible lines of chickens. One of the candidate lincRNAs, termed linc-satb1, was found to play a crucial role in MD immune response by regulating a nearby protein-coding gene SATB1. Thus, our results manifested that lincRNAs may exert considerable influence on MDV-induced T cell tumorigenesis and provide a rich resource for hypothesis-driven functional studies to reveal genetic mechanisms underlying susceptibility to tumorigenesis.

Transcriptomic Profiling of Spleen in Grass-Fed and Grain-Fed Angus Cattle

The grass-fed cattle obtain nutrients directly from pastures containing limited assimilable energy but abundant amount of fiber; by contrast, grain-fed steers receive a diet that is comprised mainly of grains and serves as an efficient source of high-digestible energy. Besides energy, these two types of diet differ in a large number of nutritional components. Additionally, animals maintained on rich-energy regimen are more likely to develop metabolic disorders and infectious diseases than pasture raised individuals. Thus, we hypothesize that spleen–a relevant immune organ–may function differently under disparate regimes. The objective of this study was to find the differentially expressed genes in the spleen of grass-fed and grain-fed steers, and furtherly explore the potential involved biopathways. Through RNA sequencing (RNA-Seq), we detected 123 differentially expressed genes. Based on these genes, we performed an Ingenuity Pathway Analysis (IPA) and identified 9 significant molecular networks and 13 enriched biological pathways. Two of the pathways, Nur77 signaling in T lymphocytes and calcium-induced T lymphocyte apoptosis which are immune related, contain a pair of genes HLA-DRA and NR4A1 with dramatically altered expression level. Collectively, our results provided valuable insights into understanding the molecular mechanism of spleen under varied feeding regimens.

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

BackgroundMarek’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.ResultsIn 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.ConclusionOur 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.