CRISPR-Cas systems are widespread accessory elements across bacterial and archaeal plasmids.

Many prokaryotes encode CRISPR-Cas systems as immune protection against mobile genetic elements (MGEs), yet a number of MGEs also harbor CRISPR-Cas components. With a few exceptions, CRISPR-Cas loci encoded on MGEs are uncharted and a comprehensive analysis of their distribution, prevalence, diversity, and function is lacking. Here, we systematically investigated CRISPR-Cas loci across the largest curated collection of natural bacterial and archaeal plasmids. CRISPR-Cas loci are widely but heterogeneously distributed across plasmids and, in comparison to host chromosomes, their mean prevalence per Mbp is higher and their distribution is distinct. Furthermore, the spacer content of plasmid CRISPRs exhibits a strong targeting bias towards other plasmids, while chromosomal arrays are enriched with virus-targeting spacers. These contrasting targeting preferences highlight the genetic independence of plasmids and suggest a major role for mediating plasmid-plasmid conflicts. Altogether, CRISPR-Cas are frequent accessory components of many plasmids, which is an overlooked phenomenon that possibly facilitates their dissemination across microbiomes.

Type IV CRISPR-Cas systems are highly diverse and involved in competition between plasmids.

CRISPR-Cas systems provide prokaryotes with adaptive immune functions against viruses and other genetic parasites. In contrast to all other types of CRISPR-Cas systems, type IV has remained largely overlooked. Here, we describe a previously uncharted diversity of type IV gene cassettes, primarily encoded by plasmid-like elements from diverse prokaryotic taxa. Remarkably, via a comprehensive analysis of their CRISPR spacer content, these systems were found to exhibit a strong bias towards the targeting of other plasmids. Our data indicate that the functions of type IV systems have diverged from those of other host-related CRISPR-Cas immune systems to adopt a role in mediating conflicts between plasmids. Furthermore, we find evidence for cross-talk between certain type IV and type I CRISPR-Cas systems that co-exist intracellularly, thus providing a simple answer to the enigmatic absence of type IV adaptation modules. Collectively, our results lead to the expansion and reclassification of type IV systems and provide novel insights into the biological function and evolution of these elusive systems.

Comprehensive search for accessory proteins encoded with archaeal and bacterial type III CRISPR-cas gene cassettes reveals 39 new cas gene families.

A study was undertaken to identify conserved proteins that are encoded adjacent to cas gene cassettes of Type III CRISPR-Cas (Clustered Regularly Interspaced Short Palindromic Repeats - CRISPR associated) interference modules. Type III modules have been shown to target and degrade dsDNA, ssDNA and ssRNA and are frequently intertwined with cofunctional accessory genes, including genes encoding CRISPR-associated Rossman Fold (CARF) domains. Using a comparative genomics approach, and defining a Type III association score accounting for coevolution and specificity of flanking genes, we identified and classified 39 new Type III associated gene families. Most archaeal and bacterial Type III modules were seen to be flanked by several accessory genes, around half of which did not encode CARF domains and remain of unknown function. Northern blotting and interference assays in Synechocystis confirmed that one particular non-CARF accessory protein family was involved in crRNA maturation. Non-CARF accessory genes were generally diverse, encoding nuclease, helicase, protease, ATPase, transporter and transmembrane domains with some encoding no known domains. We infer that additional families of non-CARF accessory proteins remain to be found. The method employed is scalable for potential application to metagenomic data once automated pipelines for annotation of CRISPR-Cas systems have been developed. All accessory genes found in this study are presented online in a readily accessible and searchable format for researchers to audit their model organism of choice: http://accessory.crispr.dk .

Stable maintenance of the rudivirus SIRV3 in a carrier state in Sulfolobus islandicus despite activation of the CRISPR-Cas immune response by a second virus SMV1.

Carrier state viral infection constitutes an equilibrium state in which a limited fraction of a cellular population is infected while the remaining cells are transiently resistant to infection. This type of infection has been characterized for several bacteriophages but not, to date, for archaeal viruses. Here we demonstrate that the rudivirus SIRV3 can produce a host-dependent carrier state infection in the model crenarchaeon Sulfolobus. SIRV3 only infected a fraction of a Sulfolobus islandicus REY15A culture over several days during which host growth was unimpaired and no chromosomal DNA degradation was observed. CRISPR spacer acquisition from SIRV3 DNA was induced by coinfecting with the monocaudavirus SMV1 and it was coincident with increased transcript levels from subtype I-A adaptation and interference cas genes. However, this response did not significantly affect the carrier state infection of SIRV3 and both viruses were maintained in the culture over 12 days during which SIRV3 anti-CRISPR genes were shown to be expressed. Transcriptome and proteome analyses demonstrated that most SIRV3 genes were expressed at varying levels over time whereas SMV1 gene expression was generally low. The study yields insights into the basis for the stable infection of SIRV3 and the resistance to the different host CRISPR-Cas interference mechanisms. It also provides a rationale for the commonly observed coinfection of archaeal cells by different viruses in natural environments.

Diverse CRISPR-Cas responses and dramatic cellular DNA changes and cell death in pKEF9-conjugated Sulfolobus species.

The Sulfolobales host a unique family of crenarchaeal conjugative plasmids some of which undergo complex rearrangements intracellularly. Here we examined the conjugation cycle of pKEF9 in the recipient strain Sulfolobus islandicus REY15A. The plasmid conjugated and replicated rapidly generating high average copy numbers which led to strong growth retardation that was coincident with activation of CRISPR-Cas adaptation. Simultaneously, intracellular DNA was extensively degraded and this also occurred in a conjugated Δcas6 mutant lacking a CRISPR-Cas immune response. Furthermore, the integrated forms of pKEF9 in the donor Sulfolobus solfataricus P1 and recipient host were specifically corrupted by transposable orfB elements, indicative of a dual mechanism for inactivating free and integrated forms of the plasmid. In addition, the CRISPR locus of pKEF9 was progressively deleted when conjugated into the recipient strain. Factors influencing activation of CRISPR-Cas adaptation in the recipient strain are considered, including the first evidence for a possible priming effect in Sulfolobus The 3-Mbp genome sequence of the donor P1 strain is presented.

Transcriptome changes in STSV2-infected Sulfolobus islandicus REY15A undergoing continuous CRISPR spacer acquisition.

A transcriptome study was performed on Sulfolobus islandicus REY15A actively undergoing CRISPR spacer acquisition from the crenarchaeal monocaudavirus STSV2 in rich and basal media over a 6 day period. Spacer acquisition preceded strong host growth retardation, altered transcriptional activity of four different CRISPR-Cas modules and changes in viral copy numbers, and with significant differences in the two media. Transcript levels of proteins involved in the cell cycle were reduced, whereas those of DNA replication, DNA repair, transcriptional regulation and some antitoxin-toxin pairs and transposases were unchanged or enhanced. Antisense RNAs were implicated in the transcriptional regulation of adaptation and interference modules of the type I-A CRISPR-Cas system, and evidence was found for the occurrence of functional co-ordination between the single CRISPR-Cas adaptation module and the functionally diverse interference modules.

Characterizing leader sequences of CRISPR loci.

MOTIVATION: The CRISPR-Cas system is an adaptive immune system in many archaea and bacteria, which provides resistance against invading genetic elements. The first phase of CRISPR-Cas immunity is called adaptation, in which small DNA fragments are excised from genetic elements and are inserted into a CRISPR array generally adjacent to its so called leader sequence at one end of the array. It has been shown that transcription initiation and adaptation signals of the CRISPR array are located within the leader. However, apart from promoters, there is very little knowledge of sequence or structural motifs or their possible functions. Leader properties have mainly been characterized through transcriptional initiation data from single organisms but large-scale characterization of leaders has remained challenging due to their low level of sequence conservation. RESULTS: We developed a method to successfully detect leader sequences by focusing on the consensus repeat of the adjacent CRISPR array and weak upstream conservation signals. We applied our tool to the analysis of a comprehensive genomic database and identified several characteristic properties of leader sequences specific to archaea and bacteria, ranging from distinctive sizes to preferential indel localization. CRISPRleader provides a full annotation of the CRISPR array, its strand orientation as well as conserved core leader boundaries that can be uploaded to any genome browser. In addition, it outputs reader-friendly HTML pages for conserved leader clusters from our database. AVAILABILITY AND IMPLEMENTATION: CRISPRleader and multiple sequence alignments for all 195 leader clusters are available at http://www.bioinf.uni-freiburg.de/Software/CRISPRleader/ CONTACT: costa@informatik.uni-freiburg.de or backofen@informatik.uni-freiburg.de SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Major and minor crRNA annealing sites facilitate low stringency DNA protospacer binding prior to Type I-A CRISPR-Cas interference in Sulfolobus.

The stringency of crRNA-protospacer DNA base pair matching required for effective CRISPR-Cas interference is relatively low in crenarchaeal Sulfolobus species in contrast to that required in some bacteria. To understand its biological significance we studied crRNA-protospacer interactions in Sulfolobus islandicus REY15A which carries multiple, and functionally diverse, interference complexes. A range of mismatches were introduced into a vector-borne protospacer that was identical to spacer 1 of CRISPR locus 2, with a cognate CCN PAM sequence. Two important crRNA annealing regions were identified on the 39 bp protospacer, a strong primary site centered on nucleotides 3 - 7 and a weaker secondary site at nucleotides 21 - 25. Multiple mismatches introduced into remaining protospacer regions did not seriously impair interference. Extending the study to different protospacers demonstrated that the efficacy of the secondary site was greatest for protospacers with higher G+C contents. In addition, the interference effects were assigned specifically to the type I-A dsDNA-targeting module by repeating the experiments with mutated protospacer constructs that were transformed into an S. islandicus mutant lacking type III-Bα and III-Bß interference gene cassettes, which showed similar interference levels to those of the wild-type strain. Parallels are drawn to the involvement of 2 annealing sites for microRNAs on some eukaryal mRNAs which provide enhanced binding capacity and specificity. A biological rationale for the relatively low crRNA-protospacer base pairing stringency among the Sulfolobales is considered.

Inter-viral conflicts that exploit host CRISPR immune systems of Sulfolobus.

Infection of Sulfolobus islandicus REY15A with mixtures of different Sulfolobus viruses, including STSV2, did not induce spacer acquisition by the host CRISPR immune system. However, coinfection with the tailed fusiform viruses SMV1 and STSV2 generated hyperactive spacer acquisition in both CRISPR loci, exclusively from STSV2, with the resultant loss of STSV2 but not SMV1. SMV1 was shown to activate adaptation while itself being resistant to CRISPR-mediated adaptation and DNA interference. Exceptionally, a single clone S-1 isolated from an SMV1 + STSV2-infected culture, that carried STSV2-specific spacers and had lost STSV2 but not SMV1, acquired spacers from SMV1. This effect was also reproducible on reinfecting wild-type host cells with a variant SMV1 isolated from the S-1 culture. The SMV1 variant lacked a virion protein ORF114 that was shown to bind DNA. This study also provided evidence for: (i) limits on the maximum sizes of CRISPR loci; (ii) spacer uptake strongly retarding growth of infected cultures; (iii) protospacer selection being essentially random and non-directional, and (iv) the reversible uptake of spacers from STSV2 and SMV1. A hypothesis is presented to explain the interactive conflicts between SMV1 and the host CRISPR immune system.

CRISPRstrand: predicting repeat orientations to determine the crRNA-encoding strand at CRISPR loci.

MOTIVATION: The discovery of CRISPR-Cas systems almost 20 years ago rapidly changed our perception of the bacterial and archaeal immune systems. CRISPR loci consist of several repetitive DNA sequences called repeats, inter-spaced by stretches of variable length sequences called spacers. This CRISPR array is transcribed and processed into multiple mature RNA species (crRNAs). A single crRNA is integrated into an interference complex, together with CRISPR-associated (Cas) proteins, to bind and degrade invading nucleic acids. Although existing bioinformatics tools can recognize CRISPR loci by their characteristic repeat-spacer architecture, they generally output CRISPR arrays of ambiguous orientation and thus do not determine the strand from which crRNAs are processed. Knowledge of the correct orientation is crucial for many tasks, including the classification of CRISPR conservation, the detection of leader regions, the identification of target sites (protospacers) on invading genetic elements and the characterization of protospacer-adjacent motifs. RESULTS: We present a fast and accurate tool to determine the crRNA-encoding strand at CRISPR loci by predicting the correct orientation of repeats based on an advanced machine learning approach. Both the repeat sequence and mutation information were encoded and processed by an efficient graph kernel to learn higher-order correlations. The model was trained and tested on curated data comprising >4500 CRISPRs and yielded a remarkable performance of 0.95 AUC ROC (area under the curve of the receiver operator characteristic). In addition, we show that accurate orientation information greatly improved detection of conserved repeat sequence families and structure motifs. We integrated CRISPRstrand predictions into our CRISPRmap web server of CRISPR conservation and updated the latter to version 2.0. AVAILABILITY: CRISPRmap and CRISPRstrand are available at http://rna.informatik.uni-freiburg.de/CRISPRmap. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Solution properties of the archaeal CRISPR DNA repeat-binding homeodomain protein Cbp2.

Clustered regularly interspaced short palindromic repeats (CRISPR) form the basis of diverse adaptive immune systems directed primarily against invading genetic elements of archaea and bacteria. Cbp1 of the crenarchaeal thermoacidophilic order Sulfolobales, carrying three imperfect repeats, binds specifically to CRISPR DNA repeats and has been implicated in facilitating production of long transcripts from CRISPR loci. Here, a second related class of CRISPR DNA repeat-binding protein, denoted Cbp2, is characterized that contains two imperfect repeats and is found amongst members of the crenarchaeal thermoneutrophilic order Desulfurococcales. DNA repeat-binding properties of the Hyperthermus butylicus protein Cbp2Hb were characterized and its three-dimensional structure was determined by NMR spectroscopy. The two repeats generate helix-turn-helix structures separated by a basic linker that is implicated in facilitating high affinity DNA binding of Cbp2 by tethering the two domains. Structural studies on mutant proteins provide support for Cys(7) and Cys(28) enhancing high thermal stability of Cbp2Hb through disulphide bridge formation. Consistent with their proposed CRISPR transcriptional regulatory role, Cbp2Hb and, by inference, other Cbp1 and Cbp2 proteins are closely related in structure to homeodomain proteins with linked helix-turn-helix (HTH) domains, in particular the paired domain Pax and Myb family proteins that are involved in eukaryal transcriptional regulation.

A novel interference mechanism by a type IIIB CRISPR-Cmr module in Sulfolobus.

Recent studies on CRISPR-based adaptive immune systems have revealed extensive structural and functional diversity of the interference complexes which often coexist intracellularly. The archaeon Sulfolobus islandicus REY15A encodes three interference modules, one of type IA and two of type IIIB. Earlier we showed that type IA activity eliminated plasmid vectors carrying matching protospacers with specific CCN PAM sequences. Here we demonstrate that interference-mediated by one type IIIB module Cmr-α, and a Csx1 protein, efficiently eliminated plasmid vectors carrying matching protospacers but lacking PAM motifs. Moreover, Cmr-α-mediated interference was dependent on directional transcription of the protospacer, in contrast to the transcription-independent activities of the type IA and type IIIA DNA interference. We infer that the interference mechanism involves transcription-dependent DNA targeting. A rationale is provided for the intracellular coexistence of the different interference systems in S. islandicus REY15A which cooperate functionally by sharing a single Cas6 protein for crRNA processing and utilize crRNA products from identical CRISPR spacers.

A novel single-tailed fusiform Sulfolobus virus STSV2 infecting model Sulfolobus species.

A newly isolated single-tailed fusiform virus, Sulfolobus tengchongensis spindle-shaped virus STSV2, from Hamazui, China, is characterised. It contains a double-stranded modified DNA genome of 76,107 bp and is enveloped by a lipid membrane structure. Virions exhibit a single coat protein that forms oligomers when isolated. STSV2 is related to the single-tailed fusiform virus STSV1 and, more distantly, to the two-tailed bicaudavirus ATV. The virus can be stably cultured over long periods in laboratory strains of Sulfolobus and no evidence was found for cell lysis under different stress conditions. Therefore, it constitutes an excellent model virus for archaeal virus-host studies.

A backward view from 16S rRNA to archaea to the universal tree of life to progenotes: reminiscences of Carl Woese.

I first became aware of Carl Woese in the mid-1970s when he and George Fox criticized a few of the 16S rRNA oligonucleotide sequences emerging from Strasbourg in the 10-12 y RNA sequencing project of the first 16S rRNA from Escherichia coli, some of which we were using for assembling RNA binding sites of ribosomal proteins. When I realized that they were attempting to sequence 16S rRNAs from a range of bacteria to classify them phylogenetically, I seriously questioned their sanity. Not because of the goal, which was admirable, but because of the sheer technical difficulty, and slowness, of sequencing large RNA molecules using the original Sanger RNA sequencing method, not to mention the health hazards of regularly preparing rRNA using 20-30 mCi [ (32)P]. My view changed radically, however, with their subsequent prediction of 5S rRNA secondary structures using a phylogenetic approach. Previously, the molecular biology community had been competing to generate the maximum numbers of base pairs in the model RNA molecule E. coli 5S RNA when Fox and Woese introduced the concept of compensatory base changes based on phylogeny for defining secondary structure and applied it to 5S RNA, they found evidence for only about 50% base pairing. This approach had previously been used for tRNA secondary structure predictions but its more general significance had never been acknowledged. Carl subsequently persuaded Harry Noller to apply the same method to predicting secondary structures of the large rRNAs.

CRISPR adaptive immune systems of Archaea.

CRISPR adaptive immune systems were analyzed for all available completed genomes of archaea, which included representatives of each of the main archaeal phyla. Initially, all proteins encoded within, and proximal to, CRISPR-cas loci were clustered and analyzed using a profile-profile approach. Then cas genes were assigned to gene cassettes and to functional modules for adaptation and interference. CRISPR systems were then classified primarily on the basis of their concatenated Cas protein sequences and gene synteny of the interference modules. With few exceptions, they could be assigned to the universal Type I or Type III systems. For Type I, subtypes I-A, I-B, and I-D dominate but the data support the division of subtype I-B into two subtypes, designated I-B and I-G. About 70% of the Type III systems fall into the universal subtypes III-A and III-B but the remainder, some of which are phyla-specific, diverge significantly in Cas protein sequences, and/or gene synteny, and they are classified separately. Furthermore, a few CRISPR systems that could not be assigned to Type I or Type III are categorized as variant systems. Criteria are presented for assigning newly sequenced archaeal CRISPR systems to the different subtypes. Several accessory proteins were identified that show a specific gene linkage, especially to Type III interference modules, and these may be cofunctional with the CRISPR systems. Evidence is presented for extensive exchange having occurred between adaptation and interference modules of different archaeal CRISPR systems, indicating the wide compatibility of the functionally diverse interference complexes with the relatively conserved adaptation modules.

Modulation of CRISPR locus transcription by the repeat-binding protein Cbp1 in Sulfolobus.

CRISPR loci are essential components of the adaptive immune system of archaea and bacteria. They consist of long arrays of repeats separated by DNA spacers encoding guide RNAs (crRNA), which target foreign genetic elements. Cbp1 (CRISPR DNA repeat binding protein) binds specifically to the multiple direct repeats of CRISPR loci of members of the acidothermophilic, crenarchaeal order Sulfolobales. cbp1 gene deletion from Sulfolobus islandicus REY15A produced a strong reduction in pre-crRNA yields from CRISPR loci but did not inhibit the foreign DNA targeting capacity of the CRISPR/Cas system. Conversely, overexpression of Cbp1 in S. islandicus generated an increase in pre-crRNA yields while the level of reverse strand transcripts from CRISPR loci remained unchanged. It is proposed that Cbp1 modulates production of longer pre-crRNA transcripts from CRISPR loci. A possible mechanism is that it minimizes interference from potential transcriptional signals carried on spacers deriving from A-T-rich genetic elements and, occasionally, on DNA repeats. Supporting evidence is provided by microarray and northern blotting analyses, and publicly available whole-transcriptome data for S. solfataricus P2.

Selective and hyperactive uptake of foreign DNA by adaptive immune systems of an archaeon via two distinct mechanisms.

Central to the disparate adaptive immune systems of archaea and bacteria are clustered regularly interspaced short palindromic repeats (CRISPR). The spacer regions derive from invading genetic elements and, via RNA intermediates and associated proteins, target and cleave nucleic acids of the invader. Here we demonstrate the hyperactive uptake of hundreds of unique spacers within CRISPR loci associated with type I and IIIB immune systems of a hyperthermophilic archaeon. Infection with an environmental virus mixture resulted in the exclusive uptake of protospacers from a co-infecting putative conjugative plasmid. Spacer uptake occurred by two distinct mechanisms in only one of two CRISPR loci subfamilies present. In two loci, insertions, often multiple, occurred adjacent to the leader while in a third locus single spacers were incorporated throughout the array. Protospacer DNAs were excised from the invading genetic element immediately after CCN motifs, on either strand, with the secondary cut apparently produced by a ruler mechanism. Over a 10-week period, there was a gradual decrease in the number of wild-type cells present in the culture but the virus and putative conjugative plasmid were still propagating. The results underline the complex dynamics of CRISPR-based immune systems within a population infected with genetic elements.

SMV1 virus-induced CRISPR spacer acquisition from the conjugative plasmid pMGB1 in Sulfolobus solfataricus P2.

Organisms of the crenarchaeal order Sulfolobales carry complex CRISPR (clustered regularly interspaced short palindromic repeats) adaptive immune systems. These systems are modular and show extensive structural and functional diversity, especially in their interference complexes. The primary targets are an exceptional range of diverse viruses, many of which propagate stably within cells and follow lytic life cycles without producing cell lysis. These properties are consistent with the difficulty of activating CRISPR spacer uptake in the laboratory, but appear to conflict with the high complexity and diversity of the CRISPR immune systems that are found among the Sulfolobales. In the present article, we re-examine the first successful induction of archaeal spacer acquisition in our laboratory that occurred exclusively for the conjugative plasmid pMGB1 in Sulfolobus solfataricus P2 that was co-infected with the virus SMV1 (Sulfolobus monocaudavirus 1). Although we reaffirm that protospacer selection is essentially a random process with respect to the pMGB1 genome, we identified single spacer sequences specific for each of CRISPR loci C, D and E that, exceptionally, occurred in many sequenced clones. Moreover, the same sequence was reproducibly acquired for a given locus in independent experiments, consistent with it being the first protospacer to be selected. There was also a small protospacer bias (1.6:1) to the antisense strand of protein genes. In addition, new experiments demonstrated that spacer acquisition in the previously inactive CRISPR locus A could be induced on freeze-thawing of the infected cells, suggesting that environmental stress can facilitate activation. Coincidentally with spacer acquisition, a mobile OrfB element was deleted from pMGB1, suggesting that interplay can occur between spacer acquisition and transposition.

Protospacer recognition motifs: mixed identities and functional diversity.

Protospacer adjacent motifs (PAMs) were originally characterized for CRISPR-Cas systems that were classified on the basis of their CRISPR repeat sequences. A few short 2-5 bp sequences were identified adjacent to one end of the protospacers. Experimental and bioinformatical results linked the motif to the excision of protospacers and their insertion into CRISPR loci. Subsequently, evidence accumulated from different virus- and plasmid-targeting assays, suggesting that these motifs were also recognized during DNA interference, at least for the recently classified type I and type II CRISPR-based systems. The two processes, spacer acquisition and protospacer interference, employ different molecular mechanisms, and there is increasing evidence to suggest that the sequence motifs that are recognized, while overlapping, are unlikely to be identical. In this article, we consider the properties of PAM sequences and summarize the evidence for their dual functional roles. It is proposed to use the terms protospacer associated motif (PAM) for the conserved DNA sequence and to employ spacer acqusition motif (SAM) and target interference motif (TIM), respectively, for acquisition and interference recognition sites.

Novel insights into gene regulation of the rudivirus SIRV2 infecting Sulfolobus cells.

Microarray analysis of infection by a lytic Sulfolobus rudivirus, SIRV2, revealed both the temporal expression of viral genes and the differential regulation of host genes. A highly susceptible strain derived from Sulfolobus solfataricus P2 with a large genomic deletion spanning CRISPR clusters A to D was infected with SIRV2, and subjected to a microarray analysis. Transcripts from a few viral genes were detected at 15 min post-infection and all except one were expressed within 2 h. The earliest expressed genes were located mainly at the termini of the linear viral genome while later expressed genes were concentrated in the central region. Timing of the expression correlated with the known or predicted functions of the viral gene products and, thus, should facilitate functional characterization of many hypothetical viral genes. Evaluation of the microarray data with quantitative reverse-transcription PCR analyses of a few selected viral genes revealed a good correlation between the two methods. Expression of about 3,000 host genes was examined. Seventy-two were downregulated>2-fold that were mainly associated with stress response and vesicle formation, as well as chromosome structure maintenance, which appears to contribute to host chromosome degradation and cellular collapse. A further 76 host genes were upregulated>2-fold and they were dominated by genes associated with metabolism and membrane transport, including phosphate transport and DNA precursor synthesis. The altered transcriptional patterns suggest that the virus reprograms the host cellular machinery to facilitate its own DNA replication and to inhibit cellular processes required for defense against viruses.