A pangenome graph reference of 30 chicken genomes allows genotyping of large and complex structural variants.

BACKGROUND: The red junglefowl, the wild outgroup of domestic chickens, has historically served as a reference for genomic studies of domestic chickens. These studies have provided insight into the etiology of traits of commercial importance. However, the use of a single reference genome does not capture diversity present among modern breeds, many of which have accumulated molecular changes due to drift and selection. While reference-based resequencing is well-suited to cataloging simple variants such as single-nucleotide changes and short insertions and deletions, it is mostly inadequate to discover more complex structural variation in the genome. METHODS: We present a pangenome for the domestic chicken consisting of thirty assemblies of chickens from different breeds and research lines. RESULTS: We demonstrate how this pangenome can be used to catalog structural variants present in modern breeds and untangle complex nested variation. We show that alignment of short reads from 100 diverse wild and domestic chickens to this pangenome reduces reference bias by 38%, which affects downstream genotyping results. This approach also allows for the accurate genotyping of a large and complex pair of structural variants at the K feathering locus using short reads, which would not be possible using a linear reference. CONCLUSIONS: We expect that this new paradigm of genomic reference will allow better pinpointing of exact mutations responsible for specific phenotypes, which will in turn be necessary for breeding chickens that meet new sustainability criteria and are resilient to quickly evolving pathogen threats.

Phenotypic Characterization of Recombinant Marek's Disease Virus in Live Birds Validates Polymorphisms Associated with Virulence.

Marek's disease (MD) is a highly infectious lymphoproliferative disease in chickens with a significant economic impact. Mardivirus gallidalpha 2, also known as Marek's disease virus (MDV), is the causative pathogen and has been categorized based on its virulence rank into four pathotypes: mild (m), virulent (v), very virulent (vv), and very virulent plus (vv+). A prior comparative genomics study suggested that several single-nucleotide polymorphisms (SNPs) and genes in the MDV genome are associated with virulence, including nonsynonymous (ns) SNPs in eight open reading frames (ORF): UL22, UL36, UL37, UL41, UL43, R-LORF8, R-LORF7, and ICP4. To validate the contribution of these nsSNPs to virulence, the vv+MDV strain 686 genome was modified by replacing nucleotides with those observed in the vMDV strains. Pathogenicity studies indicated that these substitutions reduced the MD incidence and increased the survival of challenged birds. Furthermore, using the best-fit pathotyping method to rank the virulence, the modified vv+MDV 686 viruses resulted in a pathotype similar to the vvMDV Md5 strain. Thus, these results support our hypothesis that SNPs in one or more of these ORFs are associated with virulence but, as a group, are not sufficient to result in a vMDV pathotype, suggesting that there are additional variants in the MDV genome associated with virulence, which is not surprising given this complex phenotype and our previous finding of additional variants and SNPs associated with virulence.

Targeted Ablation of Exon 2 of the Avian Leukosis Virus-A (ALV-A) Receptor Gene in a Chicken Fibroblast Cell Line by CRISPR Abrogates ALV-A Infection.

The U.S. Department of Agriculture Avian Disease and Oncology Laboratory currently relies on live birds of specific genetic backgrounds for producing chicken-embryo fibroblasts that are used for the diagnosis and subtyping of field isolates associated with avian leukosis virus (ALV) outbreaks. As an alternative to maintaining live animals for this purpose, we are currently developing cell lines capable of achieving the same result by ablation of the entry receptors utilized by ALV strains. We used CRISPR-Cas9 on the cell fibroblast-derived cell line DF-1 to disrupt the tva gene, which encodes the receptor required for binding and entry of ALV-A into cells. We ultimately identified seven DF-1 clones that had biallelic and homozygous indels at the Cas9 target site, exon 2 of tva. When tested in vitro for their ability to host ALV-A, the five clones that had frameshift mutations that disrupted the Tva protein were unable to support ALV-A replication. This result clearly demonstrates that modified cell lines can be used as part of a battery of tests to determine ALV subtype for isolate characterization, thus eliminating the need for live birds.

Nota de investigación- La ablación dirigida del exón 2 del gene del receptor del virus de la leucosis aviar A (ALV-A) en una línea celular de fibroblastos de pollo mediante CRISPR anula la infección por ALV-A. El Laboratorio de Oncología y Enfermedades Aviares del Departamento de Agricultura de los Estados Unidos. actualmente depende de aves vivas con antecedentes genéticos específicos para producir fibroblastos de embrión de pollo que se utilizan para el diagnóstico y la subtipificación de aislamientos de campo asociados con brotes del virus de la leucosis aviar (ALV). Como alternativa al mantenimiento de animales vivos para este propósito, actualmente se están desarrollando líneas celulares capaces de lograr el mismo resultado mediante la ablación de los receptores de entrada utilizados por las cepas ALV. Se utilizó el método repeticiones palindrómicas cortas agrupadas y regularmente interespaciadas o CRISPR-Cas9 en la línea celular DF-1 derivada de fibroblastos para interrumpir el gene Tva, que codifica el receptor requerido para la unión y entrada de ALV-A en las células. Finalmente, se identificaron siete clones de DF-1 que tenían inserciones y deleciones (indeles) bialélicos y homocigóticos en el sitio blanco Cas9, exón 2 del gene tva. Cuando se probó in vitro su capacidad para albergar ALV-A, los cinco clones que tenían mutaciones que involucraban al marco de lectura y que interrumpieron la proteína Tva no pudieron admitir la replicación de ALVA. Este resultado demuestra claramente que las líneas celulares modificadas se pueden utilizar como parte de una batería de pruebas para determinar el subtipo de ALV para la caracterización de los aislamientos, eliminando así la necesidad de aves vivas.

The immune cell landscape and response of Marek's disease resistant and susceptible chickens infected with Marek's disease virus.

Genetically resistant or susceptible chickens to Marek's disease (MD) have been widely used models to identify the molecular determinants of these phenotypes. However, these prior studies lacked the basic identification and understanding of immune cell types that could be translated toward improved MD control. To gain insights into specific immune cell types and their responses to Marek's disease virus (MDV) infection, we used single-cell RNA sequencing (scRNAseq) on splenic cells from MD resistant and susceptible birds. In total, 14,378 cells formed clusters that identified various immune cell types. Lymphocytes, specifically T cell subtypes, were the most abundant with significant proportional changes in some subtypes upon infection. The largest number of differentially expressed genes (DEG) response was seen in granulocytes, while macrophage DEGs differed in directionality by subtype and line. Among the most DEG in almost all immune cell types were granzyme and granulysin, both associated with cell-perforating processes. Protein interactive network analyses revealed multiple overlapping canonical pathways within both lymphoid and myeloid cell lineages. This initial estimation of the chicken immune cell type landscape and its accompanying response will greatly aid efforts in identifying specific cell types and improving our knowledge of host response to viral infection.

Contribution of the TCRß Repertoire to Marek's Disease Genetic Resistance in the Chicken.

Marek's disease (MD) is a lymphoproliferative disease of chickens induced by Marek's disease virus (MDV), an oncogenic α-herpesvirus. MDV has increased in virulence, prompting continued efforts in both improved vaccines and enhanced genetic resistance. Model pairs of genetically MD-resistant and MD-susceptible chickens that were either MHC-matched or MHC-congenic allowed characterization of T cell receptor (TCR) repertoires associated with MDV infection. MD-resistant chickens showed higher usage of Vß-1 TCRs than susceptible chickens in both the CD8 and CD4 subsets in the MHC-matched model, and in the CD8 subset only in the MHC-congenic model, with a shift towards Vß-1+ CD8 cells during MDV infection. Long and short read sequencing identified divergent TCRß loci between MHC-matched MD-resistant and MD-susceptible chickens, with MD-resistant chickens having more TCR Vß1 genes. TCR Vß1 CDR1 haplotype usage in MD-resistant x MD-susceptible F1 birds by RNAseq indicated that the most commonly used CDR1 variant was unique to the MD-susceptible line, suggesting that selection for MD resistance in the MHC-matched model optimized the TCR repertoire away from dominant recognition of one or more B2 haplotype MHC molecules. Finally, TCR downregulation during MDV infection in the MHC-matched model was strongest in the MD-susceptible line, and MDV reactivation downregulated TCR expression in a tumor cell line.

Fast-forwarding evolution-Accelerated adaptation in a proofreading-deficient hypermutator herpesvirus.

Evolution relies on the availability of genetic diversity for fitness-based selection. However, most deoxyribonucleic acid (DNA) viruses employ DNA polymerases (Pol) capable of exonucleolytic proofreading to limit mutation rates during DNA replication. The relative genetic stability produced by high-fidelity genome replication can make studying DNA virus adaptation and evolution an intensive endeavor, especially in slowly replicating viruses. Here, we present a proofreading-impaired Pol mutant (Y547S) of Marek's disease virus that exhibits a hypermutator phenotype while maintaining unimpaired growth in vitro and wild-type (WT)-like pathogenicity in vivo. At the same time, mutation frequencies observed in Y547S virus populations are 2-5-fold higher compared to the parental WT virus. We find that Y547S adapts faster to growth in originally non-permissive cells, evades pressure conferred by antiviral inhibitors more efficiently, and is more easily attenuated by serial passage in cultured cells compared to WT. Our results suggest that hypermutator viruses can serve as a tool to accelerate evolutionary processes and help identify key genetic changes required for adaptation to novel host cells and resistance to antiviral therapy. Similarly, the rapid attenuation achieved through adaptation of hypermutators to growth in cell culture enables identification of genetic changes underlying attenuation and virulence, knowledge that could practically exploited, e.g. in the rational design of vaccines.

Codon Deoptimization of UL54 in meq-Deleted Marek's Disease Vaccine Candidate Eliminates Lymphoid Atrophy but Reduces Vaccinal Protection.

Marek's disease (MD) is an oncogenic, lymphoproliferative, and highly contagious disease of chickens. Its etiologic agent is the alphaherpesvirus Marek's disease virus (MDV, Gallid alphaherpesvirus 2), and it is a chronic and ubiquitous problem for the poultry industry with significant economic impact in the United States and worldwide. We have previously demonstrated that MDV attenuated by dicodon deoptimization of the UL54 gene results in reduced gene product accumulation in vitro, with reduced viral genome copy number upon infection and reduced atrophy of bursa and thymus in vivo as well. In this report we detail our attempts to use the same attenuation strategy on a meq-deleted MDV mutant, rMd5B40ΔMeq. Unlike the wild-type rMd5B40 virus the rMd5B40ΔMeq is no longer oncogenic, but infected birds experience an unacceptable amount of bursa and thymus atrophy (BTA). We produced two meq-deleted MDV recombinants with a dicodon-deoptimized UL54 (rMd5B40ΔMeq/UL54deop1 and -deop2) and tested their tendency to cause BTA and to serve as a protective vaccine. We found that, although dicodon deoptimization of the UL54 gene results in a virus that spares the infected animal from atrophy of the bursa and thymus, the meq-deleted UL54-deoptimized recombinant is also less protective than the meq-deleted virus without UL54 deoptimization, the HVT + SB1 combination vaccine, or the Rispens (CVI988) vaccine.

Attenuation of Marek's disease virus by codon pair deoptimization of a core gene.

Marek's disease virus (MDV) is an oncogenic alphaherpesvirus of Gallus gallus, the domesticated chicken. Control strategies rely upon vaccination with live attenuated viruses of antigenically similar avian herpesviruses or attenuated strains of MDV. Recent studies in other viruses have shown that recoding certain viral genes to employ synonymous but rarely-used codon pairs resulted in viral attenuation. We deoptimized two MDV proteins, UL54/ICP27 and UL49/VP22, and demonstrate that the more severely deoptimized variant of UL54 accumulates significantly less gene product in vitro. Using these UL54 deoptimized mutants, we further demonstrate that animals infected with the UL54-recoded recombinant virus exhibited decreased viral genome copy number in lymphocytes, reduced lymphoid atrophy and reduced tumor incidence. This study demonstrates that codon pair deoptimization of a single viral gene can produce attenuated strains of MDV. This approach may be useful as a rational way of making novel live attenuated virus vaccines for MDV.