We investigated genome folding across the eukaryotic tree of life. We find two types of three-dimensional (3D) genome architectures at the chromosome scale. Each type appears and disappears repeatedly during eukaryotic evolution. The type of genome architecture that an organism exhibits correlates with the absence of condensin II subunits. Moreover, condensin II depletion converts the architecture of the human genome to a state resembling that seen in organisms such as fungi or mosquitoes. In this state, centromeres cluster together at nucleoli, and heterochromatin domains merge. We propose a physical model in which lengthwise compaction of chromosomes by condensin II during mitosis determines chromosome-scale genome architecture, with effects that are retained during the subsequent interphase. This mechanism likely has been conserved since the last common ancestor of all eukaryotes.

Endogenous viral elements (EVEs) are genetic remnants of viruses that have integrated into host genomes millions of years ago and retained as heritable elements passed on to offspring until present-day. As a result, EVEs provide an opportunity to analyze the genomes of extinct viruses utilizing these genomic viral fossils to study evolution of viruses over large timescales. Analysis of sequences from near full-length EVEs of dependoparvoviral origin identified within three mammalian taxa, Whippomorpha (whales and hippos), Vespertillionade (smooth-nosed bats), and Lagomorpha (rabbits, hares, and pikas), indicates that distinct ancestral dependoparvovirus species integrated into these host genomes approximately 77-23 million years ago (MYA). These ancestral viruses are unique relative to modern adeno-associated viruses (AAVs), and distinct from extant species of genus Dependoparvovirus. These EVE sequences show characteristics previously unseen in modern, mammalian AAVs, but instead appear more similar to the more primitive, autonomously-replicating and pathogenic waterfowl dependoparvoviruses. Phylogeny reconstruction suggests that the whippomorph EVE orthologue derives from exogenous ancestors of autonomous and highly pathogenic dependoparvovirus lineages, believed to have uniquely co-evolved with waterfowl birds to present date. In contrast, ancestors of the two other mammalian orthologues (Lagomorpha and Vespertilionidae) likely shared the same lineage as all other known mammalian exogenous AAVs. Comparative in silico analysis of the EVE genomes revealed remarkable overall conservation of AAV rep and cap genes, despite millions of years of integration within the host germline. Modeling these proteins identified unexpected variety, even between orthologues, in previously defined capsid viral protein (VP) variable regions (VRs), especially in those related to the threefold and fivefold symmetry axes of the capsid. Moreover, the normally well-conserved phospholipase A2 (PLA2) domain of the predicted minor VP1 also exhibited a high degree of sequence variance. These findings may indicate unique biological properties for these virus “fossils” relative to extant dependoparvoviruses and suggest key regions to explore within capsid sequences that may confer novel properties for engineered gene therapy vectors based upon paleovirology data.

Marek's disease (MD) is a lymphoproliferative disease of chickens caused by the oncogenic Gallid herpesvirus 2, commonly known as Marek's disease virus (MDV). MD vaccines, the primary control method, are often generated by repeated in vitro serial passage of this highly cell-associated virus to attenuate virulent MDV strains. To understand the genetic basis of attenuation, we used experimental evolution by serially passing three virulent MDV replicates generated from an infectious bacterial artificial chromosome (BAC) clone. All replicates became completely or highly attenuated, indicating that de novo mutation, and not selection among quasispecies existing in a strain, is the primary driving force for the reduction in virulence. Sequence analysis of the attenuated replicates revealed 41 to 95 single-nucleotide variants (SNVs) at 2% or higher frequency in each population and several candidate genes containing high-frequency, nonsynonymous mutations. Five candidate mutations were incorporated into recombinant viruses to determine their in vivo effect. SNVs within UL42 (DNA polymerase auxiliary subunit) and UL46 (tegument) had no measurable influence, while two independent mutations in LORF2 (a gene of unknown function) improved survival time of birds but did not alter disease incidence. A fifth SNV located within UL5 (helicase-primase subunit) greatly reduced in vivo viral replication, increased survival time of birds, and resulted in only 0 to 11% disease incidence. This study shows that multiple genes, often within pathways involving DNA replication and transcriptional regulation, are involved in de novo attenuation of MDV and provides targets for the rational design of future MD vaccines. IMPORTANCEMarek's disease virus (MDV) is a very important pathogen in chickens that costs the worldwide poultry industry $1 billion to $2 billion annually. Marek's disease (MD) vaccines, the primary control method, are often produced by passing virulent strains in cell culture until attenuated. To understand this process, we identified all the changes in the viral genome that occurred during repeated cell passage. We find that a single mutation in the UL5 gene, which encodes a viral protein necessary for DNA replication, reduces disease incidence by 90% or more. In addition, other candidate genes were identified. This information should lead to the development of more effective and rationally designed MD vaccines leading to improved animal health and welfare and lower costs to consumers.

Marek’s disease virus (MDV) is the causative agent for Marek’s disease (MD), which is characterized by T-cell lymphomas in chickens. While the viral Meq oncogene is necessary for transformation, it is insufficient, as not every bird infected with virulent MDV goes on to develop a gross tumor. Thus, we postulated that the chicken genome contains cancer driver genes; i.e., ones with somatic mutations that promote tumors, as is the case for most human cancers. To test this hypothesis, MD tumors and matching control tissues were sequenced. Using a custom bioinformatics pipeline, 9 of the 22 tumors analyzed contained one or more somatic mutation in Ikaros (IKFZ1), a transcription factor that acts as the master regulator of lymphocyte development. The mutations found were in key Zn-finger DNA-binding domains that also commonly occur in human cancers such as B-cell acute lymphoblastic leukemia (B-ALL). To validate that IKFZ1 was a cancer driver gene, recombinant MDVs that expressed either wild-type or a mutated Ikaros allele were used to infect chickens. As predicted, birds infected with MDV expressing the mutant Ikaros allele had high tumor incidences (~90%), while there were only a few minute tumors (~12%) produced in birds infected with the virus expressing wild-type Ikaros. Thus, in addition to Meq, key somatic mutations in Ikaros or other potential cancer driver genes in the chicken genome are necessary for MDV to induce lymphomas.