Laszlo Zsak

The Genome of Fowlpox Virus

ABSTRACT Here we present the genomic sequence, with analysis, of a pathogenic fowlpox virus (FPV). The 288-kbp FPV genome consists of a central coding region bounded by identical 9.5-kbp inverted terminal repeats and contains 260 open reading frames, of which 101 exhibit similarity to genes of known function. Comparison of the FPV genome with those of other chordopoxviruses (ChPVs) revealed 65 conserved gene homologues, encoding proteins involved in transcription and mRNA biogenesis, nucleotide metabolism, DNA replication and repair, protein processing, and virion structure. Comparison of the FPV genome with those of other ChPVs revealed extensive genome colinearity which is interrupted in FPV by a translocation and a major inversion, the presence of multiple and in some cases large gene families, and novel cellular homologues. Large numbers of cellular homologues together with 10 multigene families largely account for the marked size difference between the FPV genome (260 to 309 kbp) and other known ChPV genomes (178 to 191 kbp). Predicted proteins with putative functions involving immune evasion included eight natural killer cell receptors, four CC chemokines, three G-protein-coupled receptors, two β nerve growth factors, transforming growth factor β, interleukin-18-binding protein, semaphorin, and five serine proteinase inhibitors (serpins). Other potential FPV host range proteins included homologues of those involved in apoptosis (e.g., Bcl-2 protein), cell growth (e.g., epidermal growth factor domain protein), tissue tropism (e.g., ankyrin repeat-containing gene family, N1R/p28 gene family, and a T10 homologue), and avian host range (e.g., a protein present in both fowl adenovirus and Mareks disease virus). The presence of homologues of genes encoding proteins involved in steroid biogenesis (e.g., hydroxysteroid dehydrogenase), antioxidant functions (e.g., glutathione peroxidase), vesicle trafficking (e.g., two α-type soluble NSF attachment proteins), and other, unknown conserved cellular processes (e.g., Hal3 domain protein and GSN1/SUR4) suggests that significant modification of host cell function occurs upon viral infection. The presence of a cyclobutane pyrimidine dimer photolyase homologue in FPV suggests the presence of a photoreactivation DNA repair pathway. This diverse complement of genes with likely host range functions in FPV suggests significant viral adaptation to the avian host.

The Genome of a Very Virulent Marek's Disease Virus

ABSTRACT Here we present the first complete genomic sequence, with analysis, of a very virulent strain of Mareks disease virus serotype 1 (MDV1), Md5. The genome is 177,874 bp and is predicted to encode 103 proteins. MDV1 is colinear with the prototypic alphaherpesvirus herpes simplex virus type 1 (HSV-1) within the unique long (UL) region, and it is most similar at the amino acid level to MDV2, herpesvirus of turkeys (HVT), and nonavian herpesviruses equine herpesviruses 1 and 4. MDV1 encodes 55 HSV-1 UL homologues together with 6 additional UL proteins that are absent in nonavian herpesviruses. The unique short (US) region is colinear with and has greater than 99% nucleotide identity to that of MDV1 strain GA; however, an extra nucleotide sequence at the Md5 US/short terminal repeat boundary results in a shorter US region and the presence of a second gene (encoding MDV097) similar to the SORF2 gene. MD5, like HVT, encodes an ICP4 homologue that contains a 900-amino-acid amino-terminal extension not found in other herpesviruses. Putative virulence and host range gene products include the oncoprotein MEQ, oncogenicity-associated phosphoproteins pp38 and pp24, a lipase homologue, a CxC chemokine, and unique proteins of unknown function MDV087 and MDV097 (SORF2 homologues) and MDV093 (SORF4). Consistent with its virulent phenotype, Md5 contains only two copies of the 132-bp repeat which has previously been associated with viral attenuation and loss of oncogenicity.

The Genomes of Sheeppox and Goatpox Viruses

ABSTRACT Sheeppox virus (SPPV) and goatpox virus (GTPV), members of the Capripoxvirus genus of the Poxviridae, are etiologic agents of important diseases of sheep and goats in northern and central Africa, southwest and central Asia, and the Indian subcontinent. Here we report the genomic sequence and comparative analysis of five SPPV and GTPV isolates, including three pathogenic field isolates and two attenuated vaccine viruses. SPPV and GTPV genomes are approximately 150 kbp and are strikingly similar to each other, exhibiting 96% nucleotide identity over their entire length. Wild-type genomes share at least 147 putative genes, including conserved poxvirus replicative and structural genes and genes likely involved in virulence and host range. SPPV and GTPV genomes are very similar to that of lumpy skin disease virus (LSDV), sharing 97% nucleotide identity. All SPPV and GTPV genes are present in LSDV. Notably in both SPPV and GTPV genomes, nine LSDV genes with likely virulence and host range functions are disrupted, including a gene unique to LSDV (LSDV132) and genes similar to those coding for interleukin-1 receptor, myxoma virus M003.2 and M004.1 genes (two copies each), and vaccinia virus F11L, N2L, and K7L genes. The absence of these genes in SPPV and GTPV suggests a significant role for them in the bovine host range. SPPV and GTPV genomes contain specific nucleotide differences, suggesting they are phylogenetically distinct. Relatively few genomic changes in SPPV and GTPV vaccine viruses account for viral attenuation, because they contain 71 and 7 genomic changes compared to their respective field strains. Notable genetic changes include mutation or disruption of genes with predicted functions involving virulence and host range, including two ankyrin repeat proteins in SPPV and three kelch-like proteins in GTPV. These comparative genomic data indicate the close genetic relationship among capripoxviruses, and they suggest that SPPV and GTPV are distinct and likely derived from an LSDV-like ancestor.

The Genome of Canarypox Virus

ABSTRACT Here we present the genomic sequence, with analysis, of a canarypox virus (CNPV). The 365-kbp CNPV genome contains 328 potential genes in a central region and in 6.5-kbp inverted terminal repeats. Comparison with the previously characterized fowlpox virus (FWPV) genome revealed avipoxvirus-specific genomic features, including large genomic rearrangements relative to other chordopoxviruses and novel cellular homologues and gene families. CNPV also contains many genomic differences with FWPV, including over 75 kbp of additional sequence, 39 genes lacking FWPV homologues, and an average of 47% amino acid divergence between homologues. Differences occur primarily in terminal and, notably, localized internal genomic regions and suggest significant genomic diversity among avipoxviruses. Divergent regions contain gene families, which overall comprise over 49% of the CNPV genome and include genes encoding 51 proteins containing ankyrin repeats, 26 N1R/p28-like proteins, and potential immunomodulatory proteins, including those similar to transforming growth factor β and β-nerve growth factor. CNPV genes lacking homologues in FWPV encode proteins similar to ubiquitin, interleukin-10-like proteins, tumor necrosis factor receptor, PIR1 RNA phosphatase, thioredoxin binding protein, MyD116 domain proteins, circovirus Rep proteins, and the nucleotide metabolism proteins thymidylate kinase and ribonucleotide reductase small subunit. These data reveal genomic differences likely affecting differences in avipoxvirus virulence and host range, and they will likely be useful for the design of improved vaccine vectors.

Genome of Lumpy Skin Disease Virus

ABSTRACT Lumpy skin disease virus (LSDV), a member of the capripoxvirus genus of the Poxviridae, is the etiologic agent of an important disease of cattle in Africa. Here we report the genomic sequence of LSDV. The 151-kbp LSDV genome consists of a central coding region bounded by identical 2.4 kbp-inverted terminal repeats and contains 156 putative genes. Comparison of LSDV with chordopoxviruses of other genera reveals 146 conserved genes which encode proteins involved in transcription and mRNA biogenesis, nucleotide metabolism, DNA replication, protein processing, virion structure and assembly, and viral virulence and host range. In the central genomic region, LSDV genes share a high degree of colinearity and amino acid identity (average of 65%) with genes of other known mammalian poxviruses, particularly suipoxvirus, yatapoxvirus, and leporipoxviruses. In the terminal regions, colinearity is disrupted and poxvirus homologues are either absent or share a lower percentage of amino acid identity (average of 43%). Most of these differences involve genes and gene families with likely functions involving viral virulence and host range. Although LSDV resembles leporipoxviruses in gene content and organization, it also contains homologues of interleukin-10 (IL-10), IL-1 binding proteins, G protein-coupled CC chemokine receptor, and epidermal growth factor-like protein which are found in other poxvirus genera. These data show that although LSDV is closely related to other members of the Chordopoxvirinae, it contains a unique complement of genes responsible for viral host range and virulence.

The Genome of Turkey Herpesvirus

ABSTRACT Here we present the first complete genomic sequence of Mareks disease virus serotype 3 (MDV3), also known as turkey herpesvirus (HVT). The 159,160-bp genome encodes an estimated 99 putative proteins and resembles alphaherpesviruses in genomic organization and gene content. HVT is very similar to MDV1 and MDV2 within the unique long (UL) and unique short (US) genomic regions, where homologous genes share a high degree of colinearity and their proteins share a high level of amino acid identity. Within the UL region, HVT contains 57 genes with homologues found in herpes simplex virus type 1 (HSV-1), six genes with homologues found only in MDV, and two genes (HVT068 and HVT070 genes) which are unique to HVT. The HVT US region is 2.2 kb shorter than that of MDV1 (Md5 strain) due to the absence of an MDV093 (SORF4) homologue and to differences at the UL/short repeat (RS) boundary. HVT lacks a homologue of MDV087, a protein encoded at the UL/RS boundary of MDV1 (Md5), and it contains two homologues of MDV096 (glycoprotein E) in the RS. HVT RS are 1,039 bp longer than those in MDV1, and with the exception of an ICP4 gene homologue, the gene content is different from that of MDV1. Six unique genes, including a homologue of the antiapoptotic gene Bcl-2, are found in the RS. This is the first reported Bcl-2 homologue in an alphaherpesvirus. HVT long repeats (RL) are 7,407 bp shorter than those in MDV1 and do not contain homologues of MDV1 genes with functions involving virulence, oncogenicity, and immune evasion. HVT lacks homologues of MDV1 oncoprotein MEQ, CxC chemokine, oncogenicity-associated phosphoprotein pp24, and conserved domains of phosphoprotein pp38. These significant genomic differences in and adjacent to RS and RL regions likely account for the differences in host range, virulence, and oncogenicity between nonpathogenic HVT and highly pathogenic MDV1.

The Genome of Swinepox Virus

ABSTRACT Swinepox virus (SWPV), the sole member of the Suipoxvirus genus of the Poxviridae, is the etiologic agent of a worldwide disease specific for swine. Here we report the genomic sequence of SWPV. The 146-kbp SWPV genome consists of a central coding region bounded by identical 3.7-kbp inverted terminal repeats and contains 150 putative genes. Comparison of SWPV with chordopoxviruses reveals 146 conserved genes encoding proteins involved in basic replicative functions, viral virulence, host range, and immune evasion. Notably, these include genes with similarity to genes for gamma interferon (IFN-γ) receptor, IFN resistance protein, interleukin-18 binding protein, IFN-α/β binding protein, extracellular enveloped virus host range protein, dUTPase, hydroxysteroid dehydrogenase, superoxide dismutase, serpin, herpesvirus major histocompatibility complex inhibitor, ectromelia virus macrophage host range protein, myxoma virus M011L, variola virus B22R, four ankyrin repeat proteins, three kelch-like proteins, five vaccinia virus (VV) A52R-like family proteins, and two G protein-coupled receptors. The most conserved genomic region is centrally located and corresponds to the VV region located between genes F9L and A38L. Within the terminal 13 kbp, colinearity is disrupted and multiple poxvirus gene homologues are absent or share a lower percentage of amino acid identity. Most of these differences involve genes and gene families with likely functions involving viral virulence and host range. Three open reading frames (SPV018, SPV019. and SPV020) are unique for SWPV. Phylogenetic analysis, genome organization, and amino acid identity indicate that SWPV is most closely related to the capripoxvirus lumpy skin disease virus, followed by the yatapoxvirus yaba-like disease virus and the leporipoxviruses. The gene complement of SWPV better defines Suipoxvirus within the Chordopoxvirinae subfamily and provides a basis for future genetic comparisons.

Genome of Horsepox Virus

ABSTRACT Here we present the genomic sequence of horsepox virus (HSPV) isolate MNR-76, an orthopoxvirus (OPV) isolated in 1976 from diseased Mongolian horses. The 212-kbp genome contained 7.5-kbp inverted terminal repeats and lacked extensive terminal tandem repetition. HSPV contained 236 open reading frames (ORFs) with similarity to those in other OPVs, with those in the central 100-kbp region most conserved relative to other OPVs. Phylogenetic analysis of the conserved region indicated that HSPV is closely related to sequenced isolates of vaccinia virus (VACV) and rabbitpox virus, clearly grouping together these VACV-like viruses. Fifty-four HSPV ORFs likely represented fragments of 25 orthologous OPV genes, including in the central region the only known fragmented form of an OPV ribonucleotide reductase large subunit gene. In terminal genomic regions, HSPV lacked full-length homologues of genes variably fragmented in other VACV-like viruses but was unique in fragmentation of the homologue of VACV strain Copenhagen B6R, a gene intact in other known VACV-like viruses. Notably, HSPV contained in terminal genomic regions 17 kbp of OPV-like sequence absent in known VACV-like viruses, including fragments of genes intact in other OPVs and approximately 1.4 kb of sequence present only in cowpox virus (CPXV). HSPV also contained seven full-length genes fragmented or missing in other VACV-like viruses, including intact homologues of the CPXV strain GRI-90 D2L/I4R CrmB and D13L CD30-like tumor necrosis factor receptors, D3L/I3R and C1L ankyrin repeat proteins, B19R kelch-like protein, D7L BTB/POZ domain protein, and B22R variola virus B22R-like protein. These results indicated that HSPV contains unique genomic features likely contributing to a unique virulence/host range phenotype. They also indicated that while closely related to known VACV-like viruses, HSPV contains additional, potentially ancestral sequences absent in other VACV-like viruses.

African Swine Fever Virus Multigene Family 360 and 530 Genes Affect Host Interferon Response

ABSTRACT African swine fever virus (ASFV) multigene family 360 and 530 (MGF360/530) genes affect viral growth in macrophage cell cultures and virulence in pigs (L. Zsak, Z. Lu, T. G. Burrage, J. G. Neilan, G. F. Kutish, D. M. Moore, and D. L. Rock, J. Virol. 75:3066-3076, 2001). The mechanism by which these novel genes affect virus-host interactions is unknown. To define MGF360/530 gene function, we compared macrophage transcriptional responses following infection with parental ASFV (Pr4) and an MGF360/530 deletion mutant (Pr4Δ35). A swine cDNA microarray containing 7,712 macrophage cDNA clones was used to compare the transcriptional profiles of swine macrophages infected with Pr4 and Pr4Δ35 at 3 and 6 h postinfection (hpi). While at 3 hpi most (7,564) of the genes had similar expression levels in cells infected with either virus, 38 genes had significantly increased (>2.0-fold, P < 0.05) mRNA levels in Pr4Δ35-infected macrophages. Similar up-regulation of these genes was observed at 6 hpi. Viral infection was required for this induced transcriptional response. Most Pr4Δ35 up-regulated genes were part of a type I interferon (IFN) response or were genes that are normally induced by double-stranded RNA and/or viral infection. These included monocyte chemoattractant protein, transmembrane protein 3, tetratricopeptide repeat protein 1, a ubiquitin-like 17-kDa protein, ubiquitin-specific protease ISG43, an RNA helicase DEAD box protein, GTP-binding MX protein, the cytokine IP-10, and the PKR activator PACT. Differential expression of IFN early-response genes in Pr4Δ35 relative to Pr4 was confirmed by Northern blot analysis and real-time PCR. Analysis of IFN-α mRNA and secreted IFN-α levels at 3, 8, and 24 hpi revealed undetectable IFN-α in mock- and Pr4-infected macrophages but significant IFN-α levels at 24 hpi in Pr4Δ35-infected macrophages. The absence of IFN-α in Pr4-infected macrophages suggests that MGF360/530 genes either directly or indirectly suppress a type I IFN response. An inability to suppress host type I IFN responses may account for the growth defect of Pr4Δ35 in macrophages and its attenuation in swine.

Metagenomic analysis of the turkey gut RNA virus community

Viral enteric disease is an ongoing economic burden to poultry producers worldwide, and despite considerable research, no single virus has emerged as a likely causative agent and target for prevention and control efforts. Historically, electron microscopy has been used to identify suspect viruses, with many small, round viruses eluding classification based solely on morphology. National and regional surveys using molecular diagnostics have revealed that suspect viruses continuously circulate in United States poultry, with many viruses appearing concomitantly and in healthy birds. High-throughput nucleic acid pyrosequencing is a powerful diagnostic technology capable of determining the full genomic repertoire present in a complex environmental sample. We utilized the Roche/454 Life Sciences GS-FLX platform to compile an RNA virus metagenome from turkey flocks experiencing enteric disease. This approach yielded numerous sequences homologous to viruses in the BLAST nr protein database, many of which have not been described in turkeys. Our analysis of this turkey gut RNA metagenome focuses in particular on the turkey-origin members of the Picornavirales, the Caliciviridae, and the turkey Picobirnaviruses.