scRNA seq of an F1 cross of Marek’s disease resistant and susceptible chickens identifies allele specific expression signatures enriched in transcription modulators

Abstract Marek’s disease (MD), a T cell lymphoma disease in chickens, is caused by the Marek’s disease virus (MDV) found ubiquitously in the poultry industry. Genetically resistant Line 63 (L6) and susceptible Line 72 (L7) chickens have been instrumental to research on avian immune system response t...

Full description

Saved in:
Bibliographic Details
Main Authors: Deborah Velez-Irizarry, Hans Cheng, Cari Hearn
Format: Article
Language:English
Published: Nature Portfolio 2025-01-01
Series:Scientific Reports
Subjects:
Online Access:https://doi.org/10.1038/s41598-025-86174-w
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1832571760059351040
author Deborah Velez-Irizarry
Hans Cheng
Cari Hearn
author_facet Deborah Velez-Irizarry
Hans Cheng
Cari Hearn
author_sort Deborah Velez-Irizarry
collection DOAJ
description Abstract Marek’s disease (MD), a T cell lymphoma disease in chickens, is caused by the Marek’s disease virus (MDV) found ubiquitously in the poultry industry. Genetically resistant Line 63 (L6) and susceptible Line 72 (L7) chickens have been instrumental to research on avian immune system response to MDV infection. In this study we characterized molecular signatures unique to splenic immune cell types across different genetic backgrounds 6 days after infection. Using three populations, L6, L7, and an F1 cross between L6xL7, we evaluated the immune cell transcriptome of responding cell types using single cell RNA sequencing. Several MDV genes were found expressed mainly in cytotoxic T cells while ICP4 and MEQ MDV genes were expressed across infected cell types. Using the F1 we quantified allele specific expression (ASE) of biallelic SNPs and found biased expression of parental alleles specific to immune cell subtypes. We identified 22 SNPs with ASE in response to MDV infection mapped to gene rich regions surrounding 59 genes of critical importance for chromatin remodeling and transcriptional regulation. Histone deacetylase genes (HDAC1 and HDAC8) had increased expression of L6 alleles, while small nuclear RNA genes (SNORA68 and SNORA72) expressed higher levels of L7 alleles with infection in T cell subsets. SNPs with ASE also mapped genes important for an adequate immune response including GNLY (cytotoxic activity) and PDIA3 (component of MHC class I peptide loading complex), and genes known to promote viral replication (MCM5 and EIF3M). These results show that functional variants associated with susceptibility to MD may have a bigger impact in subsets of immune cell types, and by characterizing the transcriptomes of these subtypes we can unravel molecular signatures specific to MD genomic resistance.
format Article
id doaj-art-e3e704888d7c44c7b6b386732106b572
institution Kabale University
issn 2045-2322
language English
publishDate 2025-01-01
publisher Nature Portfolio
record_format Article
series Scientific Reports
spelling doaj-art-e3e704888d7c44c7b6b386732106b5722025-02-02T12:21:11ZengNature PortfolioScientific Reports2045-23222025-01-0115111610.1038/s41598-025-86174-wscRNA seq of an F1 cross of Marek’s disease resistant and susceptible chickens identifies allele specific expression signatures enriched in transcription modulatorsDeborah Velez-Irizarry0Hans Cheng1Cari Hearn2USDA, Agricultural Research Service, US National Poultry Research CenterUSDA, Agricultural Research Service, US National Poultry Research CenterUSDA, Agricultural Research Service, US National Poultry Research CenterAbstract Marek’s disease (MD), a T cell lymphoma disease in chickens, is caused by the Marek’s disease virus (MDV) found ubiquitously in the poultry industry. Genetically resistant Line 63 (L6) and susceptible Line 72 (L7) chickens have been instrumental to research on avian immune system response to MDV infection. In this study we characterized molecular signatures unique to splenic immune cell types across different genetic backgrounds 6 days after infection. Using three populations, L6, L7, and an F1 cross between L6xL7, we evaluated the immune cell transcriptome of responding cell types using single cell RNA sequencing. Several MDV genes were found expressed mainly in cytotoxic T cells while ICP4 and MEQ MDV genes were expressed across infected cell types. Using the F1 we quantified allele specific expression (ASE) of biallelic SNPs and found biased expression of parental alleles specific to immune cell subtypes. We identified 22 SNPs with ASE in response to MDV infection mapped to gene rich regions surrounding 59 genes of critical importance for chromatin remodeling and transcriptional regulation. Histone deacetylase genes (HDAC1 and HDAC8) had increased expression of L6 alleles, while small nuclear RNA genes (SNORA68 and SNORA72) expressed higher levels of L7 alleles with infection in T cell subsets. SNPs with ASE also mapped genes important for an adequate immune response including GNLY (cytotoxic activity) and PDIA3 (component of MHC class I peptide loading complex), and genes known to promote viral replication (MCM5 and EIF3M). These results show that functional variants associated with susceptibility to MD may have a bigger impact in subsets of immune cell types, and by characterizing the transcriptomes of these subtypes we can unravel molecular signatures specific to MD genomic resistance.https://doi.org/10.1038/s41598-025-86174-wMarek’s Disease VirusSingle cell RNA sequencingAllele specific expressionGenetic resistanceAvian immune response
spellingShingle Deborah Velez-Irizarry
Hans Cheng
Cari Hearn
scRNA seq of an F1 cross of Marek’s disease resistant and susceptible chickens identifies allele specific expression signatures enriched in transcription modulators
Scientific Reports
Marek’s Disease Virus
Single cell RNA sequencing
Allele specific expression
Genetic resistance
Avian immune response
title scRNA seq of an F1 cross of Marek’s disease resistant and susceptible chickens identifies allele specific expression signatures enriched in transcription modulators
title_full scRNA seq of an F1 cross of Marek’s disease resistant and susceptible chickens identifies allele specific expression signatures enriched in transcription modulators
title_fullStr scRNA seq of an F1 cross of Marek’s disease resistant and susceptible chickens identifies allele specific expression signatures enriched in transcription modulators
title_full_unstemmed scRNA seq of an F1 cross of Marek’s disease resistant and susceptible chickens identifies allele specific expression signatures enriched in transcription modulators
title_short scRNA seq of an F1 cross of Marek’s disease resistant and susceptible chickens identifies allele specific expression signatures enriched in transcription modulators
title_sort scrna seq of an f1 cross of marek s disease resistant and susceptible chickens identifies allele specific expression signatures enriched in transcription modulators
topic Marek’s Disease Virus
Single cell RNA sequencing
Allele specific expression
Genetic resistance
Avian immune response
url https://doi.org/10.1038/s41598-025-86174-w
work_keys_str_mv AT deborahvelezirizarry scrnaseqofanf1crossofmareksdiseaseresistantandsusceptiblechickensidentifiesallelespecificexpressionsignaturesenrichedintranscriptionmodulators
AT hanscheng scrnaseqofanf1crossofmareksdiseaseresistantandsusceptiblechickensidentifiesallelespecificexpressionsignaturesenrichedintranscriptionmodulators
AT carihearn scrnaseqofanf1crossofmareksdiseaseresistantandsusceptiblechickensidentifiesallelespecificexpressionsignaturesenrichedintranscriptionmodulators