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

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.

Genetic variation in chicken interferon signalling pathway genes in research lines showing differential viral resistance

Avian viruses of economic interest are a significant burden on the poultry industry, affecting production traits and resulting in mortality. Furthermore, the zoonosis of avian viruses risks pandemics developing in humans. Vaccination is the most common method of controlling viruses; however current vaccines often lack cross-protection against multiple strains of each virus. The mutagenicity of these viruses has also led to virulent strains emerging that can overcome the protection offered by vaccines. Breeding chickens with a more robust innate immune response may help in tackling current and emerging viruses. Understanding the genetic evolution of different lines will thus provide a useful tool in helping the host in the fight against pathogens. This study focuses on the interferon genes and their receptors in different chicken lines that are known to be more resistant or susceptible to particular avian viruses. Comparing genotypic differences in these core immune genes between the chicken lines may explain the phenotypic differences observed and aid the identification of causative variations. The relative resistance/susceptibility of each line to viruses of interest (Marek's disease virus, infectious bursal disease, infectious bronchitis virus and avian influenza virus) has previously been determined. Here we identify single nucleotide polymorphisms in interferons and downstream genes. Functional prediction tools were used to identify variants that may be affecting protein structure, mRNA secondary structure or transcription factor and micro-RNA binding sites. These variants were then considered in the context of the research lines and their distribution between phenotypes. We highlight 60 variants of interest in the interferon pathway genes that may account for susceptibility/resistance to viral pathogens.

Gga-miR-181a modulates ANP32A expression and inhibits MDCC-MSB-1 cell

Marek's disease (MD), a highly contagious T cell lymphoid neoplasia disease of chickens, causes huge economic losses to the poultry industry. It is the only one tumor disease which can be prevented by vaccine in chickens; therefore, MD is considered to be an excellent model to study the pathogenesis of virus-induced cancer. Recently, abundant evidences have verified that miRNAs are regulators in the process of neoplastic transformation. In our previous study on miRNome analysis of MDV-induced lymphoma in chicken, we found that gga-miR-181a was downregulated drastically in MDV-infected spleens. To further investigate the role of gga-miR-181a in MDV-induced lymphomagenesis, we performed cell migration assay, and the results suggested that gga-miR-181a suppressed the migration of MDV-transformed lymphoid cell (MSB-1). Subsequently, luciferase reporter gene assay revealed that acidic nuclear phosphoprotein 32A (ANP32A) was a functional target gene of gga-miR181a. Real-time PCR and western blot assay showed that the mRNA and protein levels of ANP32A were downregulated in gga-miR-181a mimic group at 48-h and 96-h post-transfection, respectively, indicating that ANP32A was modulated by gga-miR-181a. All the results suggested that gga-miR-181a was an inhibitor in MSB-1 cell migration. ANP32A was a direct target gene of gga-miR-181a and they were implicated in MD lymphoma tumorigenesis.

Investigation of allele specific expression in various tissues of broiler chickens using the detection tool VADT

Differential abundance of allelic transcripts in a diploid organism, commonly referred to as allele specific expression (ASE), is a biologically significant phenomenon and can be examined using single nucleotide polymorphisms (SNPs) from RNA-seq. Quantifying ASE aids in our ability to identify and understand cis-regulatory mechanisms that influence gene expression, and thereby assist in identifying causal mutations. This study examines ASE in breast muscle, abdominal fat, and liver of commercial broiler chickens using variants called from a large sub-set of the samples (n = 68). ASE analysis was performed using a custom software called VCF ASE Detection Tool (VADT), which detects ASE of biallelic SNPs using a binomial test. On average ~ 174,000 SNPs in each tissue passed our filtering criteria and were considered informative, of which ~ 24,000 (~ 14%) showed ASE. Of all ASE SNPs, only 3.7% exhibited ASE in all three tissues, with ~ 83% showing ASE specific to a single tissue. When ASE genes (genes containing ASE SNPs) were compared between tissues, the overlap among all three tissues increased to 20.1%. Our results indicate that ASE genes show tissue-specific enrichment patterns, but all three tissues showed enrichment for pathways involved in translation.

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.

Genome-wide characterization of copy number variations in the host genome in genetic resistance to Marek's disease using next generation sequencing

Background

Marek’s disease (MD) is a highly neoplastic disease primarily affecting chickens, and remains as a chronic infectious disease that threatens the poultry industry. Copy number variation (CNV) has been examined in many species and is recognized as a major source of genetic variation that directly contributes to phenotypic variation such as resistance to infectious diseases. Two highly inbred chicken lines, 6₃ (MD-resistant) and 7₂ (MD-susceptible), as well as their F₁ generation and six recombinant congenic strains (RCSs) with varied susceptibility to MD, are considered as ideal models to identify the complex mechanisms of genetic and molecular resistance to MD.

Results

In the present study, to unravel the potential genetic mechanisms underlying resistance to MD, we performed a genome-wide CNV detection using next generation sequencing on the inbred chicken lines with the assistance of CNVnator. As a result, a total of 1649 CNV regions (CNVRs) were successfully identified after merging all the nine datasets, of which 90 CNVRs were overlapped across all the chicken lines. Within these shared regions, 1360 harbored genes were identified. In addition, 55 and 44 CNVRs with 62 and 57 harbored genes were specifically identified in line 6₃ and 7₂, respectively. Bioinformatics analysis showed that the nearby genes were significantly enriched in 36 GO terms and 6 KEGG pathways including JAK/STAT signaling pathway. Ten CNVRs (nine deletions and one duplication) involved in 10 disease-related genes were selected for validation by using quantitative real-time PCR (qPCR), all of which were successfully confirmed. Finally, qPCR was also used to validate two deletion events in line 7₂ that were definitely normal in line 6₃. One high-confidence gene, IRF2 was identified as the most promising candidate gene underlying resistance and susceptibility to MD in view of its function and overlaps with data from previous study.

Conclusions

Our findings provide valuable insights for understanding the genetic mechanism of resistance to MD and the identified gene and pathway could be considered as the subject of further functional characterization.

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.

Allele-specific expression is widespread in Bos indicus muscle and affects meat quality candidate genes

Differences between the expression of the two alleles of a gene are known as allele-specific expression (ASE), a common event in the transcriptome of mammals. Despite ASE being a source of phenotypic variation, its occurrence and effects on genetic prediction of economically relevant traits are still unexplored in bovines. Furthermore, as ASE events are likely driven by cis-regulatory mutations, scanning them throughout the bovine genome represents a significant step to elucidate the mechanisms underlying gene expression regulation. To address this question in a Bos indicus population, we built the ASE profile of the skeletal muscle tissue of 190 Nelore steers, using RNA sequencing data and SNPs genotypes from the Illumina BovineHD BeadChip (770 K bp). After quality control, 820 SNPs showed at least one sample with ASE. These SNPs were widespread among all autosomal chromosomes, being 32.01% found in 3'UTR and 31.41% in coding regions. We observed a considerable variation of ASE profile among individuals, which highlighted the need for biological replicates in ASE studies. Functional analysis revealed that ASE genes play critical biological functions in the development and maintenance of muscle tissue. Additionally, some of these genes were previously reported as associated with beef production and quality traits in livestock, thus indicating a possible source of bias on genomic predictions for these traits.

RNA Sequencing revealed differentially expressed genes functionally associated with immunity and tumor suppression during latent phase infection of a vv + MDV in chickens

Very virulent plus Marek's disease (MD) virus (vv + MDV) induces tumors in relatively resistant lines of chickens and early mortality in highly susceptible lines of chickens. The vv + MDV also triggers a series of cellular responses in both types of chickens. We challenged birds sampled from a highly inbred chicken line (line 6₃) that is relatively resistant to MD and from another inbred line (line 7₂) that is highly susceptible to MD with a vv + MDV. RNA-sequencing analysis was performed with samples extracted from spleen tissues taken at 10-day and 21-day post infection (dpi). A total of 64 and 106 differentially expressed genes was identified in response to the vv + MDV challenge at latent phase in the resistant and susceptible lines of chickens, respectively. Direct comparisons between samples of the two lines identified 90 and 126 differentially expressed genes for control and MDV challenged groups, respectively. The differentially expressed gene profiles illustrated that intensive defense responses were significantly induced by vv + MDV at 10 dpi and 21 dpi but with slight changes in the resistant line. In contrast, vv + MDV induced a measurable suppression of gene expression associated with host defense at 10 dpi but followed by an apparent activation of the defense response at 21 dpi in the susceptible line of chickens. The observed difference in gene expression between the two genetic lines of chickens in response to MDV challenge during the latent phase provided a piece of indirect evidence that time points for MDV reactivation differ between the genetic lines of chickens with different levels of genetic resistance to MD. Early MDV reactivation might be necessary and potent to host defense system readiness for damage control of tumorigenesis and disease progression, which consequently results in measurable differences in phenotypic characteristics including early mortality (8 to 20 dpi) and tumor incidence between the resistant and susceptible lines of chickens. Combining differential gene expression patterns with reported GO function terms and quantitative trait loci, a total of 27 top genes was selected as highly promising candidate genes for genetic resistance to MD. These genes are functionally involved with virus process (F13A1 and HSP90AB1), immunity (ABCB1LB, RGS5, C10ORF58, OSF-2, MMP7, CXCL12, GAL1, GAL2, GAL7, HVCN1, PDE4D, IL4I1, PARP9, EOMES, MPEG1, PDK4, CCLI10, K60 and FST), and tumor suppression (ADAMTS2, LXN, ARRDC3, WNT7A, CLDN1 and HPGD). It is anticipated that these findings will facilitate advancement in the fundamental understanding on mechanisms of genetic resistance to MD. In addition, such advancement may also provide insights on tumor virus-induced tumorigenesis in general and help the research community recognize MD study may serve as a good model for oncology study involving tumor viruses.

Allele-Specific Expression of CD4+ T Cells in Response to Marek's Disease Virus Infection

Marek's disease (MD) is a T cell lymphoma disease induced by Marek's disease virus (MDV), a highly oncogenic α herpesvirus primarily affecting chickens. MD is a chronic infectious disease that threatens the poultry industry. However, the mechanisms of genetic resistance for MD are complex and not completely understood. In this study, to identify high-confidence candidate genes of MD genetic resistance, high throughput sequencing (RNA-seq) was used to obtain transcriptomic data of CD4+ T cells isolated from MDV-infected and non-infected groups of two reciprocal crosses of individuals mating by two highly inbred chicken lines (63 MD-resistant and 72 MD-susceptible). After RNA-seq analysis with two biological replicates in each group, we identified 61 and 123 single nucleotide polymorphisms (SNPs) (false discovery rate (FDR) < 0.05) annotated in 39 and 132 genes in intercrosses 63 × 72 and 72 × 63, respectively, which exhibited allele-specific expression (ASE) in response to MDV infection. Similarly, we identified 62 and 79 SNPs annotated in 66 and 96 genes in infected and non-infected groups, respectively. We identified 534 and 1543 differentially expressed genes (DEGs) (FDR < 0.05) related to MDV infection in intercrosses 63 × 72 and 72 × 63, respectively. We also identified 328 and 20 DEGs in infected and non-infected groups, respectively. The qRT-PCR using seven DEGs further verified our results of RNA-seq analysis. The qRT-PCR of 11 important ASE genes was performed for gene functional validation in CD4+ T cells and tumors. Combining the analyses, six genes (MCL1, SLC43A2, PDE3B, ADAM33, BLB1, and DMB2), especially MCL1, were highlighted as the candidate genes with the potential to be involved in MDV infection. Gene-set enrichment analysis revealed that many ASE genes are linked to T cell activation, T cell receptor (TCR), B cell receptor (BCR), ERK/MAPK, and PI3K/AKT-mTOR signaling pathways, which play potentially important roles in MDV infection. Our approach underlines the importance of comprehensive functional studies for gaining valuable biological insight into the genetic factors behind MD and other complex traits, and our findings provide additional insights into the mechanisms of MD and disease resistance breeding in poultry.

Genetic Basis of Phenotypic Differences Between Chinese Yunling Black Goats and Nubian Goats Revealed by Allele-Specific Expression in Their F1 Hybrids

Chinese Yunling black goats and African Nubian goats are divergent breeds showing significant differences in body size, milk production, and environmental adaptation. However, the genetic mechanisms underlying these phenotypic differences remain to be elucidated. In this report, we provide a detailed portrait of allele-specific expression (ASE) from 54 RNA-Seq analyses across six tissues from nine F1 hybrid offspring generated by crossing the two breeds combined with 13 genomes of the two breeds. We identified a total of 524 genes with ASE, which are involved in bone development, muscle cell differentiation, and the regulation of lipid metabolic processes. We further found that 38 genes with ASE were also under directional selection by comparing 13 genomes of the two breeds; these 38 genes play important roles in metabolism, immune responses, and the adaptation to hot and humid environments. In conclusion, our study shows that the exploration of genes with ASE in F1 hybrids provides an efficient way to understand the genetic basis underlying the phenotypic differences of two diverse goat breeds.

Macrophages from Susceptible and Resistant Chicken Lines have Different Transcriptomes following Marek's Disease Virus Infection

Despite successful control by vaccination, Marek’s disease (MD) has continued evolving to greater virulence over recent years. To control MD, selection and breeding of MD-resistant chickens might be a suitable option. MHC-congenic inbred chicken lines, 6₁ and 7₂, are highly resistant and susceptible to MD, respectively, but the cellular and genetic basis for these phenotypes is unknown. Marek’s disease virus (MDV) infects macrophages, B-cells, and activated T-cells in vivo. This study investigates the cellular basis of resistance to MD in vitro with the hypothesis that resistance is determined by cells active during the innate immune response. Chicken bone marrow-derived macrophages from lines 6₁ and 7₂ were infected with MDV in vitro. Flow cytometry showed that a higher percentage of macrophages were infected in line 7₂ than in line 6₁. A transcriptomic study followed by in silico functional analysis of differentially expressed genes was then carried out between the two lines pre- and post-infection. Analysis supports the hypothesis that macrophages from susceptible and resistant chicken lines display a marked difference in their transcriptome following MDV infection. Resistance to infection, differential activation of biological pathways, and suppression of oncogenic potential are among host defense strategies identified in macrophages from resistant chickens.

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.

Application of omics technologies for a deeper insight into quali-quantitative production traits in broiler chickens: A review

The poultry industry is continuously facing substantial and different challenges such as the increasing cost of feed ingredients, the European Union's ban of antibiotic as growth promoters, the antimicrobial resistance and the high incidence of muscle myopathies and breast meat abnormalities. In the last decade, there has been an extraordinary development of many genomic techniques able to describe global variation of genes, proteins and metabolites expression level. Proper application of these cutting-edge omics technologies (mainly transcriptomics, proteomics and metabolomics) paves the possibility to understand much useful information about the biological processes and pathways behind different complex traits of chickens. The current review aimed to highlight some important knowledge achieved through the application of omics technologies and proteo-genomics data in the field of feed efficiency, nutrition, meat quality and disease resistance in broiler chickens.

Pattern recognition receptor genes expression profiling in indigenous chickens of India and White Leghorn

Pattern recognition receptors (PRR) such as Toll-like receptors, NOD-like receptors, RIG-I helicase receptors, and C-type lectin receptors play a critical role in innate immunity as a first line of defense against invading pathogens through recognition of pathogen and/or damage-associated molecular patterns. Genetic makeup of birds is known to play a role in resistance or susceptibility to various infectious diseases. Therefore, the present study was carried out to elucidate the differential expression of PRR and some of the cytokine genes in peripheral blood mononuclear cells of indigenous chicken breeds such as Ghagus and Nicobari and an exotic chicken breed, White Leghorn (WLH). The stability of expression of reference genes in peripheral blood mononuclear cells of 3 breeds was first determined using NormFinder and BestKeeper programs. NormFinder determined B2M and G6PDH reference genes as the best combination with stability value of 0.38. Out of total 14 genes studied, expression of ten genes was found to be significantly different among 3 breeds after normalization with these reference genes. Ghagus breed showed higher level of expression of TLR1LB, TLR7, NOD1, NOD5, B-Lec, IFNβ, IL1β, and IL8 genes when compared to Nicobari breed. Further, Ghagus showed higher expression of TLR1LB, MDA5, LGP2, B-Lec, IL1β, and IL8 genes as compared to WLH breed. Higher expression of LGP2 and MDA5 genes was observed in Nicobari compared to the WLH breed while higher expression of TLR7, NOD1, NOD5, and IFNβ genes was observed in WLH as compared to Nicobari breed. No difference was observed in the expression of TLR1LA, TLR3, B-NK, and IFNα genes among 3 breeds. Study revealed significant breed effect in expression profile of PRR and some of the cytokine genes and Ghagus breed seems to have better expression profile of these genes linked to the innate immunity when compared to the WLH and Nicobar breeds.

RNA-Seq Analyses Identify Frequent Allele Specific Expression and No Evidence of Genomic Imprinting in Specific Embryonic Tissues of Chicken

Epigenetic and genetic cis-regulatory elements in diploid organisms may cause allele specific expression (ASE) – unequal expression of the two chromosomal gene copies. Genomic imprinting is an intriguing type of ASE in which some genes are expressed monoallelically from either the paternal allele or maternal allele as a result of epigenetic modifications. Imprinted genes have been identified in several animal species and are frequently associated with embryonic development and growth. Whether genomic imprinting exists in chickens remains debatable, as previous studies have reported conflicting evidence. Albeit no genomic imprinting has been reported in the chicken embryo as a whole, we interrogated the existence or absence of genomic imprinting in the 12-day-old chicken embryonic brain and liver by examining ASE in F1 reciprocal crosses of two highly inbred chicken lines (Fayoumi and Leghorn). We identified 5197 and 4638 ASE SNPs, corresponding to 18.3% and 17.3% of the genes with a detectable expression in the embryonic brain and liver, respectively. There was no evidence detected of genomic imprinting in 12-day-old embryonic brain and liver. While ruling out the possibility of imprinted Z-chromosome inactivation, our results indicated that Z-linked gene expression is partially compensated between sexes in chickens.

Comparative transcriptomics of genetically divergent lines of chickens in response to Marek's disease virus challenge at cytolytic phase

Marek’s disease (MD), caused by Marek’s disease virus (MDV), remains an economically significant threat to the poultry industry worldwide. Genetic resistance to MD is a promising alternative strategy to augment current control measures (vaccination and management). However, only a few functional genes reportedly conferring MD resistance have been identified. Here, we performed a comparative transcriptomics analysis of two highly inbred yet genetically divergent lines of chickens (line 6₃ and 7₂) that are resistant and susceptible to MD, respectively, in response to a very virulent plus strain of MDV (vv+MDV) challenge at cytolytic phase. A total of 203 DEGs in response to MDV challenge were identified in the two lines. Of these, 96 DEGs were in common for both lines, in addition to 36 and 71 DEGs that were specific for line 6₃ and 7₂, respectively. Functional enrichment analysis results showed the DEGs were significantly enriched in GO terms and pathways associated with immune response. Especially, the four DEGs, FGA, ALB, FN1, and F13A1 that reportedly facilitate virus invasion or immunosuppression, were found to be significantly up-regulated in the susceptible line 7₂ but down-regulated in the resistant line 6₃ birds. These results provide new resources for future studies to further elucidate the genetic mechanism conferring MD resistance.

Differences in CD8αα and cecal microbiome community during proliferation and late cytolytic phases of Marek's disease virus infection are associated with genetic resistance to Marek's disease

Marek's disease (MD) is an important neoplastic disease of chickens caused by Marek ': s disease virus (MDV), a highly oncogenic alphaherpesvirus. In this study using two chicken lines, one resistant and another susceptible to MD, splenic T cells and cecal microbiome were profiled to gain a better understanding of primary differences in these lines. The percent of splenic CD4⁺ T cells were similar regardless of MDV challenge status in both bird lines. In contrast, CD8αα profiles were different (P < 0.005) between chicken lines under naïve status and under MDV challenge, suggesting that CD8αα T cells play a key role in mediating MDV infection. Microbiome composition was different between naïve resistant (Blautia spp.) and susceptible birds (Streptococcus spp.) (P < 0.05) during initial colonization. With MDV challenge, both chicken lines showed lower numbers of beneficial Faecalibacterium spp. and increased number of Lactobacillus spp. Metabolic profiles between naïve chicken types were similar but with MDV challenge, there were differences in metabolism in both chicken lines, with amino acid metabolism impacted in resistant birds and lipid metabolism in susceptible birds. These results provide insights into immune response and potential interplay with the microbiome during infection with an oncogenic virus.

High concentrations of atmospheric ammonia induce alterations of gene expression in the breast muscle of broilers (Gallus gallus) based on RNA-Seq

BACKGROUND

High concentrations of atmospheric ammonia are one of the key environmental stressors affecting broiler production performance, which causes remarkable economic losses as well as potential welfare problems of the broiler industry. Previous reports demonstrated that high levels of ammonia altered body fat distribution and meat quality of broilers. However, the molecular mechanisms and metabolic pathways in breast muscle altered by high concentrations of ambient ammonia exposure on broilers are still unknown.

RESULTS

This study utilized RNA-Seq to compare the transcriptomes of breast muscles to identify differentially enriched genes in broilers exposed to high and low concentrations of atmospheric ammonia. A total of 267 promising candidate genes were identified by differential expression analysis, among which 67 genes were up-regulated and 189 genes were down-regulated. Bioinformatics analysis suggested that the up and down-regulation of these genes were involved in the following two categories of cellular pathways and metabolisms: Steroid biosynthesis (gga00100) and peroxisome proliferator-activated receptor (PPAR) signaling pathway (gga03320), which both participated in the lipid metabolism processes.

CONCLUSIONS

This study suggests that longtime exposure to high concentrations of aerial ammonia can change fat content in breast muscle, meat quality and palatability via altering expression level of genes participating in important lipid metabolism pathways. These findings have provided novel insights into our understanding of molecular mechanisms of breast muscles exposed to ammonia in broilers. This study provides new information that could be used for genetic breeding and nutritional intervention in production practice of broilers industry in the future.

Genome-wide identification of copy number variations between two chicken lines that differ in genetic resistance to Marek's disease

Background

Copy number variation (CNV) is a major source of genome polymorphism that directly contributes to phenotypic variation such as resistance to infectious diseases. Lines 6₃ and 7₂ are two highly inbred experimental chicken lines that differ greatly in susceptibility to Marek’s disease (MD), and have been used extensively in efforts to identify the genetic and molecular basis for genetic resistance to MD. Using next generation sequencing, we present a genome-wide assessment of CNVs that are potentially associated with genetic resistance to MD.

Methods

Three chickens randomly selected from each line were sequenced to an average depth of 20×. Two popular software, CNVnator and Pindel, were used to call genomic CNVs separately. The results were combined to obtain a union set of genomic CNVs in the two chicken lines.

Results

A total of 5,680 CNV regions (CNVRs) were identified after merging the two datasets, of which 1,546 and 1,866 were specific to the MD resistant or susceptible line, respectively. Over half of the line-specific CNVRs were shared by 2 or more chickens, reflecting the reduced diversity in both inbred lines. The CNVRs fixed in the susceptible lines were significantly enriched in genes involved in MAPK signaling pathway. We also found 67 CNVRs overlapping with 62 genes previously shown to be strong candidates of the underlying genes responsible for the susceptibility to MD.

Conclusions

Our findings provide new insights into the genetic architecture of the two chicken lines and additional evidence that MAPK signaling pathway may play an important role in host response to MD virus infection. The rich source of line-specific CNVs is valuable for future disease-related association studies in the two chicken lines.

Electronic supplementary material

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

Fine mapping of QTL and genomic prediction using allele-specific expression SNPs demonstrates that the complex trait of genetic resistance to Marek's disease is predominantly determined by transcriptional regulation

Background

Marek’s disease (MD) is a lymphoproliferative disease of poultry induced by Marek’s disease virus (MDV), a highly oncogenic alphaherpesvirus. Identifying the underlying genes conferring MD genetic resistance is desired for more efficacious control measures including genomic selection, which requires accurately identified genetic markers throughout the chicken genome.

Methods

Hypothesizing that variants located in transcriptional regulatory regions are the main mechanism underlying this complex trait, a genome-wide association study was conducted by genotyping a ~1,000 bird MD resource population derived from experimental inbred layers with SNPs containing 1,824 previously identified allele-specific expression (ASE) SNPs in response to MDV infection as well as 3,097 random SNPs equally spaced throughout the chicken genome. Based on the calculated associations, genomic predictions were determined for 200 roosters and selected sires had their progeny tested for Marek’s disease incidence.

Results

Our analyses indicate that these ASE SNPs account for more than 83 % of the genetic variance and exhibit nearly all the highest associations. To validate these findings, 200 roosters had their genetic merit predicted from the ASE SNPs only, and the top 30 and bottom 30 ranked roosters were reciprocally mated to random hens. The resulting progeny showed that after only one generation of bidirectional selection, there was a 22 % difference in MD incidence and this approach gave a 125 % increase in accuracy compared to current pedigree-based estimates.

Conclusions

We conclude that variation in transcriptional regulation is the major driving cause for genetic resistance to MD, and ASE SNPs identify the underlying genes and are sufficiently linked to the causative polymorphisms that they can be used for accurate genomic prediction as well as help define the underlying molecular basis. Furthermore, this approach should be applicable to other complex traits.

Electronic supplementary material

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

Genome-wide identification of allele-specific expression in response to Streptococcus suis 2 infection in two differentially susceptible pig breeds

Although allele expression imbalance has been recognized in many species, and strongly linked to diseases, no whole transcriptome allele imbalance has been detected in pigs during pathogen infections. The pathogen Streptococcus suis 2 (SS2) causes serious zoonotic disease. Different pig breeds show differential susceptibility/resistance to pathogen infection, but the biological insight is little known. Here we analyzed allele-specific expression (ASE) using the spleen transcriptome of four pigs belonging to two phenotypically different breeds after SS2 infection. The comparative analysis of allele specific SNPs between control and infected animals revealed 882 and 1096 statistically significant differentially expressed allele SNPs (criteria: ratio ≧ 2 or ≦ 0.5) in Landrace and Enshi black pig, respectively. Twenty nine allelically imbalanced SNPs were further verified by Sanger sequencing, and later six SNPs were quantified by pyrosequencing assay. The pyrosequencing results are in agreement with the RNA-seq results, except two SNPs. Looking at the role of ASE in predisposition to diseases, the discovery of causative variants by ASE analysis might help the pig industry in long term to design breeding programs for improving SS2 resistance.

A hidden Markov approach for ascertaining cSNP genotypes from RNA sequence data in the presence of allelic imbalance by exploiting linkage disequilibrium

BACKGROUND

Allelic specific expression (ASE) increases our understanding of the genetic control of gene expression and its links to phenotypic variation. ASE testing is implemented through binomial or beta-binomial tests of sequence read counts of alternative alleles at a cSNP of interest in heterozygous individuals. This requires prior ascertainment of the cSNP genotypes for all individuals. To meet the needs, we propose hidden Markov methods to call SNPs from next generation RNA sequence data when ASE possibly exists.

RESULTS

We propose two hidden Markov models (HMMs), HMM-ASE and HMM-NASE that consider or do not consider ASE, respectively, in order to improve genotyping accuracy. Both HMMs have the advantages of calling the genotypes of several SNPs simultaneously and allow mapping error which, respectively, utilize the dependence among SNPs and correct the bias due to mapping error. In addition, HMM-ASE exploits ASE information to further improve genotype accuracy when the ASE is likely to be present. Simulation results indicate that the HMMs proposed demonstrate a very good prediction accuracy in terms of controlling both the false discovery rate (FDR) and the false negative rate (FNR). When ASE is present, the HMM-ASE had a lower FNR than HMM-NASE, while both can control the false discovery rate (FDR) at a similar level. By exploiting linkage disequilibrium (LD), a real data application demonstrate that the proposed methods have better sensitivity and similar FDR in calling heterozygous SNPs than the VarScan method. Sensitivity and FDR are similar to that of the BCFtools and Beagle methods. The resulting genotypes show good properties for the estimation of the genetic parameters and ASE ratios.

CONCLUSIONS

We introduce HMMs, which are able to exploit LD and account for the ASE and mapping errors, to simultaneously call SNPs from the next generation RNA sequence data. The method introduced can reliably call for cSNP genotypes even in the presence of ASE and under low sequencing coverage. As a byproduct, the proposed method is able to provide predictions of ASE ratios for the heterozygous genotypes, which can then be used for ASE testing.

Allele Workbench: transcriptome pipeline and interactive graphics for allele-specific expression

Sequencing the transcriptome can answer various questions such as determining the transcripts expressed in a given species for a specific tissue or condition, evaluating differential expression, discovering variants, and evaluating allele-specific expression. Differential expression evaluates the expression differences between different strains, tissues, and conditions. Allele-specific expression evaluates expression differences between parental alleles. Both differential expression and allele-specific expression have been studied for heterosis (hybrid vigor), where the hybrid has improved performance over the parents for one or more traits. The Allele Workbench software was developed for a heterosis study that evaluated allele-specific expression for a mouse F1 hybrid using libraries from multiple tissues with biological replicates. This software has been made into a distributable package, which includes a pipeline, a Java interface to build the database, and a Java interface for query and display of the results. The required input is a reference genome, annotation file, and one or more RNA-Seq libraries with optional replicates. It evaluates allelic imbalance at the SNP and transcript level and flags transcripts with significant opposite directional allele-specific expression. The Java interface allows the user to view data from libraries, replicates, genes, transcripts, exons, and variants, including queries on allele imbalance for selected libraries. To determine the impact of allele-specific SNPs on protein folding, variants are annotated with their effect (e.g., missense), and the parental protein sequences may be exported for protein folding analysis. The Allele Workbench processing results in transcript files and read counts that can be used as input to the previously published Transcriptome Computational Workbench, which has a new algorithm for determining a trimmed set of gene ontology terms. The software with demo files is available from https://code.google.com/p/allele-workbench. Additionally, all software is ready for immediate use from an Atmosphere Virtual Machine Image available from the iPlant Collaborative (www.iplantcollaborative.org).

Breed-specific transcriptome response of spleen from six to eight week old piglet after infection with Streptococcus suis type 2

Different pig breeds have shown differential susceptibility to the pathogen infection; however, molecular mechanisms of the infection susceptibility are not fully understood. Streptococcus suis type 2 (SS2) is an important zoonotic pathogen. To identify the genes responsible for infection susceptibility, pigs from two different breeds (Enshi black and Landrace) were inoculated with SS2 and their spleen transcriptome profiles were investigated in the present study. The differentially expressed genes (DEGs) were analyzed from infected versus control pigs in each breed, and then compared between both pig breeds. Enshi black pig showed more DEGs than Landrace (830 vs. 611) and most of these were due to down-regulated genes (543 vs. 387). However some DEGs were uniquely expressed in one breed, some were expressed in opposite direction in both breeds. A number of candidate genes and pathways are identified which might be involved in susceptibility to SS2, for example, MMP9 and Resistin were only significantly expressed in Landrace. NPG3 and PMAP23 were up-regulated in Landrace whereas down-regulated in Enshi black. LENG8 in control Landrace have inherently higher expression than control Enshi black. IGKV6 is down-regulated in Landrace but up-regulated in Enshi black. Overall, the transcriptome profiles are consistent with the clinical signs, i.e. the Enshi black is more susceptible to SS2 infection than Landrace. This is the first study to identify differential gene expression between indigenous and modern commercial pigs after in vivo SS2 infection using RNA-seq. The significant DEGs in splenic profiles between two pig breeds suggested considerable involvement of genetic background in susceptibility to the SS2 infection in pigs.

Implications of genomic signatures in the differential vulnerability to fetal alcohol exposure in C57BL/6 and DBA/2 mice

Maternal alcohol consumption inflicts a multitude of phenotypic consequences that range from undetectable changes to severe dysmorphology. Using tightly controlled murine studies that deliver precise amounts of alcohol at discrete developmental stages, our group and other labs demonstrated in prior studies that the C57BL/6 and DBA/2 inbred mouse strains display differential susceptibility to the teratogenic effects of alcohol. Since the phenotypic diversity extends beyond the amount, dosage and timing of alcohol exposure, it is likely that an individual's genetic background contributes to the phenotypic spectrum. To identify the genomic signatures associated with these observed differences in alcohol-induced dysmorphology, we conducted a microarray-based transcriptome study that also interrogated the genomic signatures between these two lines based on genetic background and alcohol exposure. This approach is called a gene x environment (GxE) analysis; one example of a GxE interaction would be a gene whose expression level increases in C57BL/6, but decreases in DBA/2 embryos, following alcohol exposure. We identified 35 candidate genes exhibiting GxE interactions. To identify cis-acting factors that mediated these interactions, we interrogated the proximal promoters of these 35 candidates and found 241 single nucleotide variants (SNVs) in 16 promoters. Further investigation indicated that 186 SNVs (15 promoters) are predicted to alter transcription factor binding. In addition, 62 SNVs created, removed or altered the placement of a CpG dinucleotide in 13 of the proximal promoters, 53 of which overlapped putative transcription factor binding sites. These 53 SNVs are also our top candidates for future studies aimed at examining the effects of alcohol on epigenetic gene regulation.

Differential expression of Toll-like receptor pathway genes in chicken embryo fibroblasts from chickens resistant and susceptible to Marek's disease

The Toll-like receptor (TLR) signaling pathway is one of the innate immune defense mechanisms against pathogens in vertebrates and invertebrates. However, the role of TLR in non-MHC genetic resistance or susceptibility to Marek's disease (MD) in the chicken is yet to be elucidated. Chicken embryo fibroblast (CEF) cells from MD susceptible and resistant lines were infected either with Marek's disease virus (MDV) or treated with polyionosinic-polycytidylic acid, a synthetic analog of dsRNA, and the expression of TLR and pro-inflammatory cytokines was studied at 8 and 36 h posttreatment by quantitative reverse transcriptase PCR. Findings of the present study reveal that MDV infection and polyionosinic-polycytidylic acid treatment significantly elevated the mRNA expression of TLR3, IL6, and IL8 in both susceptible and resistant lines. Furthermore, basal expression levels in uninfected CEF for TLR3, TLR7, and IL8 genes were significantly higher in resistant chickens compared with those of susceptible chickens. Our results suggest that TLR3 together with pro-inflammatory cytokines may play a significant role in genetic resistance to MD.

Transcriptional profiling of mEq-dependent genes in Marek's disease resistant and susceptible inbred chicken lines

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 ChIP-seq, 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.