Genome-Wide Identification of Allele-Specific Expression (ASE) in Response to Marek’s Disease Virus Infection Using Next Generation Sequencing.

Background Marek’s disease (MD), a T cell lymphoma induced by the highly oncogenic α-herpesvirus Marek’s disease virus (MDV), is the main chronic infectious disease concern threatening the poultry industry. Enhancing genetic resistance to MD in commercial poultry is an attractive method to augment MD vaccines, which is currently the control method of choice. In order to optimally implement this control strategy through marker-assisted selection (MAS) and to gain biological information, it is necessary to identify specific genes that influence MD incidence. Methods A genome-wide screen for allele-specific expression (ASE) in response to MDV infection was conducted. The highly inbred ADOL chicken lines 6 (MD resistant) and 7 (MD susceptible) were inter-mated in reciprocal crosses and half of the progeny challenged with MDV. Splenic RNA pools at a single time after infection for each treatment group point were generated, sequenced using a next generation sequencer, then analyzed for allele-specific expression (ASE). To validate and extend the results, Illumina GoldenGate assays for selected cSNPs were developed and used on all RNA samples from all 6 time points following MDV challenge. Results RNA sequencing resulted in 11-13+ million mappable reads per treatment group, 1.7+ Gb total sequence, and 22,655 high-confidence cSNPs. Analysis of these cSNPs revealed that 5360 cSNPs in 3773 genes exhibited statistically significant allelic imbalance. Of the 1536 GoldenGate assays, 1465 were successfully scored with all but 19 exhibiting evidence for allelic imbalance. Conclusions ASE is an efficient method to identify potentially all or most of the genes influencing this complex trait. The identified cSNPs can be further evaluated in resource populations to determine their allelic direction and size of effect on genetic resistance to MD as well as being directly implemented in genomic selection programs. The described method, although demonstrated in inbred chicken lines, is applicable to all traits in any diploid species, and should prove to be a simple method to identify the majority of genes controlling any complex trait

Genome-Wide Identification and Quantification of cis- and trans-Regulated Genes Responding to Marek’s Disease Virus Infection via Analysis of Allele-Specific Expression

Marek’s disease (MD) is a commercially important neoplastic disease of chickens caused by Marek’s disease virus (MDV), a naturally occurring oncogenic alphaherpesvirus. Selecting for increased genetic resistance to MD is a control strategy that can augment vaccinal control measures. To identify high-confidence candidate MD resistance genes, we conducted a genome-wide screen for allele-specific expression (ASE) amongst F1 progeny of two inbred chicken lines that differ substantially in MD resistance. High throughput sequencing was initially used to profile transcriptomes from pools of uninfected and infected individuals at 4 days post-infection to identify any genes showing ASE in response to MDV infection. RNA sequencing identified 22,655 single nucleotide polymorphisms (SNPs) of which 5,360 in 3,773 genes exhibited significant allelic imbalance. Illumina GoldenGate assays were subsequently used to quantify regulatory variation controlled at the gene (cis) and elsewhere in the genome (trans) by examining differences in expression between F1 individuals and artificial F1 RNA pools over six time periods in 1,536 of the most significant SNPs identified by RNA sequencing. Allelic imbalance as a result of cis-regulatory changes was confirmed in 861 of the 1,233 GoldenGate assays successfully examined. Furthermore we have identified seven genes that display trans-regulation only in infected animals and ∼500 SNP that show a complex interaction between cis- and trans-regulatory changes. Our results indicate ASE analyses are a powerful approach to identify regulatory variation responsible for differences in transcript abundance in genes underlying complex traits. And the genes with SNPs exhibiting ASE provide a strong foundation to further investigate the causative polymorphisms and genetic mechanisms for MD resistance. Finally, the methods used here for identifying specific genes and SNPs have practical implications for applying marker-assisted selection to complex traits that are difficult to measure in agricultural species, when expression differences are expected to control a portion of the phenotypic variance

Chromosomal integration of an avian oncogenic herpesvirus reveals telomeric preferences and evidence for lymphoma clonality

Abstract Background Herpesviruses are a major health concern for numerous organisms, including humans, causing both acute and chronic infections recurrent over an individual's lifespan. Marek's disease virus (MDV) is a highly contagious herpesvirus which causes a neoplastic condition in chicken populations. Several vertebrate-infecting herpesviruses have been shown to exist in an integrated state during latent periods of infection. However the status of MDV during latency has been a topic of debate. Results Here we employed high-resolution multi-color fluorescence in situ hybridization (FISH) to show integration of MDV at the telomeres of chicken chromosomes. Cytogenomic mapping of the chromosomal integrations allowed us to examine the clonal relationships among lymphomas within individuals, whereas analysis of tumors from multiple individuals indicated the potential for chromosomal preferences. Conclusions Our data highlight that substantive genome-level interactions between the virus and host exist, and merit consideration for their potential impact and role in key aspects of herpesvirus pathobiology including infection, latency, cellular transformation, latency-breaks and viral evolution

Improving Run Time in Three-Dimensional Reservoir Hydrodynamics and Water Quality Modeling

Source: ICHE Conference Archive - https://mdi-de.baw.de/icheArchiv

Comparison and Contrast of Genes and Biological Pathways Responding to Marek’s Disease Virus Infection Using Allele-Specific Expression and Differential Expression in Broiler and Layer Chickens.

Background Marek’s disease (MD) is a commercially important neoplastic disease of chickens caused by the Marek’s disease virus (MDV), a naturally occurring oncogenic alphaherpesvirus. Enhancing MD genetic resistance is desirable to augment current vaccines and other MD control measures. High throughput sequencing was used to profile splenic transcriptomes from individual F1 progeny infected with MDV at 4 days of age from both outbred broilers (meat-type) and inbred layer (egg-type) chicken lines that differed in MD genetic resistance. The resulting information was used to identify SNPs, genes, and biological pathways exhibiting allele-specific expression (ASE) in response to MDV infection in each type of chicken. In addition, we compared and contrasted the results of pathway analyses (ASE and differential expression (DE)) between chicken types to help inform on the biological response to MDV infection. Results With 7 individuals per line and treatment group providing high power, we identified 6,132 single nucleotide polymorphisms (SNPs) in 4,768 genes and 4,528 SNPs in 3,718 genes in broilers and layers, respectively, that exhibited ASE in response to MDV infection. Furthermore, 548 and 434 genes in broilers and layers, respectively, were found to show DE following MDV infection. Comparing the datasets, only 72 SNPs and 850 genes for ASE and 20 genes for DE were common between the two bird types. Although the chicken types used in this study were genetically different, at the pathway level, both TLR receptor and JAK/STAT signaling pathways were enriched as well as exhibiting a high proportion of ASE genes, especially at the beginning of both above mentioned regulatory pathways. Conclusions RNA sequencing with adequate biological replicates is a powerful approach to identify high confidence SNPs, genes, and pathways that are associated with transcriptional response to MDV infection. In addition, the SNPs exhibiting ASE in response to MDV infection provide a strong foundation for determining the extent to which variation in expression influences MD incidence plus yield genetic markers for genomic selection. However, given the paucity of overlap among ASE SNP sets (broilers vs. layers), it is likely that separate screens need to be incorporated for each population. Finally, comparison of gene lists obtained between these two diverse chicken types indicate the TLR and JAK/STAT signaling are conserved when responding to MDV infection and may be altered by selection of genes exhibiting ASE found at the start of each pathway

Transcriptional Profiling of MEq-Dependent Genes in Marek's Disease Resistant and Susceptible Inbred Chicken Lines

Abstract Marek's disease (MD) is an economically significant disease in chickens caused by the highly oncogenic Marek's disease virus (MDV). Understanding the genes and biological pathways that confer MD genetic resistance should lead towards the development of more disease resistant commercial poultry flocks or improved MD vaccines. MDV mEq, a bZIP transcription factor, is largely attributed to viral oncogenicity though only a few host target genes have been described, which has impeded our understanding of MDV-induced tumorigenesis. Given the importance of mEq in MDV-induced pathogenesis, we explored the role of mEq in genetic resistance to MDV. Using global transcriptome analysis and cells from MD resistant or susceptible birds, we compared the response to infection with either wild type MDV or a nononcogenic recombinant lacking mEq. As a result, we identified a number of specific genes and pathways associated with either MD resistance or susceptibility. Additionally, integrating prior information from ChIPseq, microarray analysis, and SNPs exhibiting allele-specific expression (ASE) in response to MDV infection, we were able to provide evidence for 24 genes that are polymorphic within mEq binding sites are likely to account for gene expression in an allele-specific manner and potentially for the underlying genetic differences in MD incidence

Pathogen transmission from vaccinated hosts can cause dose-dependent reduction in virulence

Many livestock and human vaccines are leaky because they block symptoms but do not prevent infection or onward transmission. This leakiness is concerning because it increases vaccination coverage required to prevent disease spread and can promote evolution of increased pathogen virulence. Despite leakiness, vaccination may reduce pathogen load, affecting disease transmission dynamics. However, the impacts on post-transmission disease development and infectiousness in contact individuals are unknown. Here, we use transmission experiments involving Marek disease virus (MDV) in chickens to show that vaccination with a leaky vaccine substantially reduces viral load in both vaccinated individuals and unvaccinated contact individuals they infect. Consequently, contact birds are less likely to develop disease symptoms or die, show less severe symptoms, and shed less infectious virus themselves, when infected by vaccinated birds. These results highlight that even partial vaccination with a leaky vaccine can have unforeseen positive consequences in controlling the spread and symptoms of disease

Genetic assessment of inbred chicken lines indicates genomic signatures of resistance to Marek\u27s 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 63 and MD-susceptible line 72 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 63 and line 72, six RCSs, and two reciprocally crossed flocks between the lines 63 and 72 (F1 63 × 72 and F1 72 × 63) 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 63. 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 63 and RCS), and these genes have been reported that are contribute to immunology and survival. Based on FST 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