On the Problem of Genome Redundancy in Viruses and Prokaryotes

A specific index of nucleotide sequence redundancy, the specific restriction length of a finite frequency dictionary, was determined for a complete set of genes in some viral genomes and a genome of a bacterium, Bacillus subtilis. The distribution of the gene number over the specific restriction length was shown to be bimodal for viral genomes and unimodal for the Bac. subtilis genome. These results agree with earlier data.     

Analysis of a restriction endonuclease map for a rabbit papillomavirus DNA

A restriction endonuclease map was constructed for rabbit papillomavirus DNA. Analysis of the viral genome by cleavage with specific endonucleases, when compared with bovine papillomavirus types 1 and 2 genome maps, revealed similarities among the three genomes. A number of cleavage sites and their relative order on each genome were found to be conserved.     

Errant processing and structural alterations of genomes present in a varicella-zoster virus vaccine.

Five minority populations of aberrant, varicella-zoster virus (VZV)-derived genomes were identified among the encapsidated DNAs obtained from the nuclear and cytoplasmic fractions of an in vitro infection initiated with a lyophilized sample of the BIKEN VZV vaccine (strain Oka). These were (i) VZV genomes, present within nuclear but not cytoplasmic viral capsids, which had been cleaved at a specific site within the short segment and which were, therefore, 3.15 megadaltons (approximately 4% of the VZV genome length) short of full length; (ii) highly deleted, repetitive VZV genomes which contained the errant cleavage site but not the usual VZV genome terminal sequences; (iii) VZV genomes into which multiples of 1 through 5 defective genome repeat units had been inserted into a homologous site; (iv) VZV genomes with additions of 0.1 or 0.18 megadaltons of DNA at both the terminal and internal ends of the short segment; and (v) VZV DNA which had lost the HindIII restriction site at map position 0.11.     

Characterization of flat revertant cells isolated from simian virus 40-transformed mouse and rat cells which contain multiple copies of viral genomes.

Eight clones of flat revertants were isolated by negative selection from simian virus 40 (SV40)-transformed mouse and rat cell lines in which two and six viral genome equivalents per cell were integrated, respectively. These revertants showed either a normal cell phenotype or a phenotype intermediate between normal and transformed cells as to cellular morphology and saturation density and were unable to grow in soft agar medium. One revertant derived from SV40-transformed mouse cells was T antigen positive, whereas the other seven revertants were T antigen negative. SV40 could be rescued only from the T-antigen-positive revertant by fusion with permissive monkey cells. The susceptibility of the revertants to retransformation by wild-type SV40 was variable among these revertants. T-antigen-negative revertants from SV40-transformed mouse cells were retransformed at a frequency of 3 to 10 times higher than their grandparental untransformed cells. In contrast, T-antigen-negative revertants from SV40-transformed rat cells could not be retransformed. The arrangement of viral genomes was analyzed by digestion of cellular DNA with restriction enzymes of different specificity, followed by detection of DNA fragments containing a viral sequence and rat cells were serially arranged within the length of about 30 kilobases, with at least two intervening cellular sequences. A head-to-tail tandem array of unit length viral genomes was present in at least one insertion site in the transformed rat cells. All of the revertants had undergone a deletion(s), and only a part of the viral genome was retained in T-antigen-negative revertants. A relatively high frequency of reversion in the transformed rat cells suggests that reversion occurs by homologous recombination between the integrated viral genomes.     

Pathogenicity of diatraea saccharalis Densovirus to Host Insets and Characterization of its Viral Genome

Pathogenicity of the Diatraea saccharalis densovirus (DsDNV) was tested on its host larvae. The results showed that up to 4 days after inoculation, no larvae mortality was observed and the infected larvae started to exhibit the infection symptoms from the fourth day. After 5 days of infection, the cumulative mortality of infected larvae increased significantly and reached 60% after 12 days and 100% after 21 days of infection, whereas that of the control group was only 10% and 20%, respectively, after same periods of infection, suggesting that the high mortality of infected larvae groups was due to the high pathogenicity of DsDNV. The size of the DsDNA was determined by Electron microscopy visualization of viral DNA molecules and gel electrophoresis of both native and endonuclease digested DNA fragments. The total length of the native DsDNA was about 5.95 kb. The DsDNV DNA was digested with 16 restriction enzymes and a restriction map of those enzymes was constructed with 41 restriction sites. Comparison of the restriction map of the DsDNV genome with those of the genomes ofJunonia coenia densovirus (JcDNV) and Galleria mellonella densovirus (GmDNV) indicated that the three densovirus genomes were found to share many identical restriction sites. Thus, most of the restriction sites of the following endonucleases Bam H I, Hha I, Xba I, Cla I, Asp 700, Spe I, Nco I and Bcl I, were found to be conserved among the three densovirus genomes. Symmetrical cleavage sites mapped at the both ends of the genome suggested the presence of inverted terminal repeats (ITRs) whose size was estimated to be about 500 bp. The similar genome size, almost identical restriction sites and presence of an ITR of about 500 bp for these three densoviruses suggested that they belong to the same group of ambisense densoviruses.     

Pathogenicity of diatraea saccharalis densovirus to host insets and characterization of its viral genome

Pathogenicity of the Diatraea saccharalis densovirus (DsDNV) was tested on its host larvae.The results showed that up to 4 days after inoculation,no larvae mortality was observed and the infected larvae started to exhibit the infection symptoms from the fourth day.After 5 days of infection,the cumulative mortality of infected larvae increased significantly and reached 60% after 12 days and 100% after 21 days of infection,whereas that of the control group was only 10% and 20%,respectively,after same periods of infection,suggesting that the high mortality of infected larvae groups was due to the high pathogenicity of DsDNV.The size of the DsDNA was determined by Electron microscopy visualization of viral DNA molecules and gel electrophoresis of both native and endonuclease digested DNA fragments.The total length of the native DsDNA was about 5.95 kb.The DsDNV DNA was digested with 16 restriction enzymes and a restriction map of those enzymes was constructed with 41 restriction sites.Comparison of the restriction map of the DsDNV genome with those of the genomes ofJunonia coenia densovirus (JcDNV) and Galleria mellonella densovirus (GmDNV) indicated that the three densovirus genomes were found to share many identical restriction sites.Thus,most of the restriction sites of the following endonucleases Bam H Ⅰ,Hha Ⅰ,Xba Ⅰ,Cla Ⅰ,Asp 700,Spe Ⅰ,Nco Ⅰ and Bcl Ⅰ,were found to be conserved among the three densovirus genomes.Symmetrical cleavage sites mapped at the both ends of the genome suggested the presence of inverted terminal repeats (ITRs) whose size was estimated to be about 500 bp.The similar genome size,almost identical restriction sites and presence of an ITR of about 500 bp for these three densoviruses suggested that they belong to the same group of ambisense densoviruses.     

Viral genome segmentation can result from a trade-off between genetic content and particle stability

The evolutionary benefit of viral genome segmentation is a classical, yet unsolved question in evolutionary biology and RNA genetics. Theoretical studies anticipated that replication of shorter RNA segments could provide a replicative advantage over standard size genomes. However, this question has remained elusive to experimentalists because of the lack of a proper viral model system. Here we present a study with a stable segmented bipartite RNA virus and its ancestor non-segmented counterpart, in an identical genomic nucleotide sequence context. Results of RNA replication, protein expression, competition experiments, and inactivation of infectious particles point to a non-replicative trait, the particle stability, as the main driver of fitness gain of segmented genomes. Accordingly, measurements of the volume occupation of the genome inside viral capsids indicate that packaging shorter genomes involves a relaxation of the packaging density that is energetically favourable. The empirical observations are used to design a computational model that predicts the existence of a critical multiplicity of infection for domination of segmented over standard types. Our experiments suggest that viral segmented genomes may have arisen as a molecular solution for the trade-off between genome length and particle stability. Genome segmentation allows maximizing the genetic content without the detrimental effect in stability derived from incresing genome length.     

Integration site of polyoma virus DNA in the inducible LPT line of polyoma-transformed rat cells

The structure of the polyoma virus (Py) integration site in the inducible LPT line of Py-transformed rat cells was determined by biochemical methods of gene mapping. LPT cell DNA was digested with various restriction enzymes. The digestion products were electrophoresed in agarose gels and transferred onto nitrocellulose sheets by Southern blotting. Fragments containing viral or cell DNA sequences, or both, were identified by hybridization with Py DNA or with a cloned flanking cell DNA probe. Cleavage of LPT DNA with enzymes that restrict the Py genome once generated linear Py DNA molecules and two fragments containing both cell and viral DNA sequences. Cleavage of LPT DNA with enzymes which do not restrict Py DNA generated series of fragments whose lengths were found to differ by increments of a whole Py genome; the smallest fragment in each series was found to be longer than the viral genome. These data indicate that LPT cultures contain Py insertions of various lengths integrated into the same chromosomal site in all the cells. The length heterogeneity of the viral insertions is due to the presence of 0, 1, 2, 3. . . Py genomes arranged in a direct tandem repeat within invariable sequences of viral DNA. Double-digestion experiments were also carried out with the above enzymes and with enzymes that cleave the Py genome at multiple sites. The data obtained in these experiments were used to construct a physical map of the integration site. This map showed that the early region of the virus remained intact even in the smallest insertion (which contains no whole duplicated genomes), whereas the late region was partially duplicated and split during integration. The smallest insertion is colinear with the Py physical map over a region including the entire Py genome and at least a part of the duplicated segment. This structure could give rise to nondefective circular viral DNA molecules by single homologous recombination events. Similar recombination events may occur at a higher frequency in the longer insertions, which include longer regions of homology, and may yield many more free viral genomes. The presence of these insertions in LPT cells could thus be one of the factors which account for the high inducibility of the LPT line.     

Comparison of the physical maps and redundant ends of the chromosomes of phages 2C, SP01, SP82 and phi e

The physical map of 2C DNA (cf. following paper in this journal) was compared to the maps of SP01, SP82 and phi e (three other Bacillus subtilis phages containing hydroxymethyluracil in place of thymine in their DNA). The overall organization of the four genomes was remarkably similar, as indicated by the topology of HaeIII and SalI cleavage segments. The proof was gathered for the presence in the four phage DNAs of large redundant ends carrying a single HaeIII recognition site. The location of the latter proved identical for 2C and SP01, but was shifted in the DNAs of SP82 and phi e. Since the redundant end components of these hydroxymethyluracil genomes are colinear, as shown by cross-hybridization studies, the shifting of the HaeIII cleavage site is presumably due to two base substitutions, suppressing an endonuclease recognition site and establishing a new site elsewhere. Relatedness between the genomes of this family of viruses was evaluated from the fraction of conserved restriction fragments. According to these calculations, 6% base substitutions have occurred within the four viral DNAs, in the course of evolution. However, specific segments of 2C DNA were not present in SP01 and phi e DNA, as shown by cross-hybridization with restriction fragments. These data indicate the occurrence of deletions, in addition to base substitutions, as evolutionary mechanisms prevailing in the genomes of this family of phages.     

鳞翅目基因组IIB型内切酶位点的统计分析

IIB型限制内切酶能够识别并切割特异酶切位点两端特定距离的DNA,形成粘性末端的30 bp左右的等长DNA片段。利用其特性与限制性酶切位点关联测序技术(RAD)相结合发展出2b-RAD简化基因组测序技术,应用于遗传图谱构建、种群遗传结构分析、性状定位以及细菌分型等多种研究领域。构建2b-RAD测序文库之前,需要对基因组中的IIB型限制内切酶位点进行预测与统计分析,制定有效的测序文库构建方案。本文利用Python语言构建分析基因组中IIB型限制内切酶位点的流程,预测并统计6个鳞翅目代表物种基因组含有的8个商业化IIB型限制内切酶的酶切位点,比较了各个基因组与IIB型限制内切酶之间含有的酶切位点总量、重复序列数量以及酶切间隔长度的关系,为在昆虫基因组中进一步试行2b-RAD研究提供了参考。     

Identification of the avian myeloblastosis virus genome. II. Restriction endonuclease analysis of DNA from lambda proviral recombinants and leukemic myeoblast clones.

Two lambda proviral DNA recombinants were characterized with a number of restriction endonucleases. One recombinant contained a complete presumptive avian myeloblastosis virus (AMV) provirus flanked by cellular sequences on either side, and the second recombinant contained 85% of a myeloblastosis-associated virus type 1 (MAV-1)-like provirus with cellular sequences adjacent to the 5' end of the provirus. Comparing the restriction maps for the proviral DNAs contained in each lambda hybrid showed that the putative AMV and MAV-1-like genomes shared identical enzyme sites for 3.6 megadaltons beginning at the 5' termini of the proviruses with respect to viral RNA. Two enzyme sites near the 3'-end of the MAV-1-like provirus were not present in the putative AMV genome. We also examined a number of leukemic myeloblast clones for proviral content and cell-provirus integration sites. The presumptive AMV provirus was present in all the leukemic myeloblast clones regardless of the endogenous proviral content of the target cells or the AMV pseudotype used for conversion. Multiple cellular sites were suitable for integration of the putative AMV genome and the helper genomes. The proviral genomes were all integrated colinearly with respect to linear viral DNA.     

Molecular evolution of bacteriophages: evidence of selection against the recognition sites of host restriction enzymes

Restriction enzymes produced by bacteria serve as a defense against invading bacteriophages, and so phages without other protection would be expected to undergo selection to eliminate recognition sites for these enzymes from their genomes. The observed frequencies of all restriction sites in the genomes of all completely sequenced DNA phages (T7, lambda, phi X174, G4, M13, f1, fd, and IKe) have been compared to expected frequencies derived from trinucleotide frequencies. Attention was focused on 6-base palindromes since they comprise the typical recognition sites for type II restriction enzymes. All of these coliphages, with the exception of lambda and G4, exhibit significant avoidance of the particular sequences that are enterobacterial restriction sites. As expected, the sequenced fraction of the genome of phi 29, a Bacillus subtilis phage, lacks Bacillus restriction sites. By contrast, the RNA phage MS2, several viruses that infect eukaryotes (EBV, adenovirus, papilloma, and SV40), and three mitochondrial genomes (human, mouse, and cow) were found not to lack restriction sites. Because the particular palindromes avoided correspond closely with the recognition sites for host enzymes and because other viruses and small genomes do not show this avoidance, it is concluded that the effect indeed results from natural selection.      

Regulation of expression and chromosomal subunit conformation of avian retrovirus genomes.

We have investigated the copy number, chromosomal subunit conformation and regulation of expression of integrated avian retrovirus genomes. Our results indicate that there are approximately two copies of the endogenous viral genomes (RAV-O) per haploid cell genome in uninfected chick embryo fibroblasts (CEF) and red blood cells (RBC). The copy number and subunit conformation (as measured by DNAasel sensitivity) of the RAV-O genomes are independent of the level of expression of these viral DNA sequences. In cells isolated from embryos of the V+, gs-chf- and gs+chf+ phenotypes, approximately one of the two viral genomes is in a DNAase l-sensitive conformation. Upon infection with an exogenous Rous sarcoma virus (PR-RSV-C), one new viral genome is integrated per haploid CEF genome. The newly integrated RSV genome is completely sensitive to DNAase l, and the subunit conformation of the endogenous viral genomes is not altered by the integration of additional exogenous proviruses. Both the endogenous and newly integrated exogenous viral genomes are present in "nu-body" structures, and the selective sensitivity of these proviral DNA sequences to DNAase l is maintained in isolated nucleosomes. Our experiments revealing the DNAase l sensitivity of one of the two RAV-O genomes in gs-chf-CEF led us to reexamine the level of viral specific RNA in CEF of various GS genotypes. We find that GS/GS CEF contain approximately 100 copies of viral RNA per cell, gs/gs CEF contain no detectable viral RNA, and the heterozygote GS/gs CEF contain approximately 50 copies of viral specific RNA per cell. These results suggest that the GS gene controls production of RAV-O RNA sequences in CEF in a "cis" fashion. In RBCs, however, the expression of the RAV-O genome is independent of the GS gene, with both GS/GS and gs/gs RBCs containing roughly equivalent amounts of viral specific RNA. Our results suggest that the chromosomal structure of the endogenous viral genes is independent of the GS gene, and that the GS gene is cis-acting and tissue-specific.     

Engineering large viral DNA genomes using the CRISPR‐Cas9 system

Manipulation of viral genomes is essential for studying viral gene function and utilizing viruses for therapy. Several techniques for viral genome engineering have been developed. Homologous recombination in virus‐infected cells has traditionally been used to edit viral genomes; however, the frequency of the expected recombination is quite low. Alternatively, large viral genomes have been edited using a bacterial artificial chromosome (BAC) plasmid system. However, cloning of large viral genomes into BAC plasmids is both laborious and time‐consuming. In addition, because it is possible for insertion into the viral genome of drug selection markers or parts of BAC plasmids to affect viral function, artificial genes sometimes need to be removed from edited viruses. Herpes simplex virus (HSV), a common DNA virus with a genome length of 152 kbp, causes labialis, genital herpes and encephalitis. Mutant HSV is a candidate for oncotherapy, in which HSV is used to kill tumor cells. In this study, the clustered regularly interspaced short palindromic repeat‐Cas9 system was used to very efficiently engineer HSV without inserting artificial genes into viral genomes. Not only gene‐ablated HSV but also gene knock‐in HSV were generated using this method. Furthermore, selection with phenotypes of edited genes promotes the isolation efficiencies of expectedly mutated viral clones. Because our method can be applied to other DNA viruses such as Epstein–Barr virus, cytomegaloviruses, vaccinia virus and baculovirus, our system will be useful for studying various types of viruses, including clinical isolates.     

Comparative genome analysis of Bacillus cereus group genomes with Bacillus subtilis

Genome features of the Bacillus cereus group genomes (representative strains of Bacillus cereus, Bacillus anthracis and Bacillus thuringiensis sub spp. israelensis) were analyzed and compared with the Bacillus subtilis genome. A core set of 1381 protein families among the four Bacillus genomes, with an additional set of 933 families common to the B. cereus group, was identified. Differences in signal transduction pathways, membrane transporters, cell surface structures, cell wall, and S-layer proteins suggesting differences in their phenotype were identified. The B. cereus group has signal transduction systems including a tyrosine kinase related to two-component system histidine kinases from B. subtilis. A model for regulation of the stress responsive sigma factor sigmaB in the B. cereus group different from the well studied regulation in B. subtilis has been proposed. Despite a high degree of chromosomal synteny among these genomes, significant differences in cell wall and spore coat proteins that contribute to the survival and adaptation in specific hosts has been identified.     

Fingerprinting of single viral genomes

We demonstrate the use of technology developed for optical mapping to acquire DNA fingerprints from single genomes for the purpose of discrimination and identification of bacteria and viruses. Single genome fingerprinting (SGF) provides not only the size but also the order of the restriction fragments, which adds another dimension to the information that can be used for discrimination. Analysis of single organisms may eliminate the need to culture cells and thereby significantly reduce analysis time. In addition, samples containing mixtures of several organisms can be analyzed. For analysis, cells are embedded in an agarose matrix, lysed, and processed to yield intact DNA. The DNA is then deposited on a derivatized glass substrate. The elongated genome is digested with a restriction enzyme and stained with the intercalating dye YOYO-1. DNA is then quantitatively imaged with a fluorescence microscope and the fragments are sized to an accuracy >or=90% by their fluorescence intensity and contour length. Single genome fingerprints were obtained from pure samples of adenovirus, from bacteriophages lambda and T4 GT7, and from a mixture of the three viral genomes. SGF will enable the fingerprinting of uncultured and unamplified samples and allow rapid identification of microorganisms with applications in forensics, medicine, public health, and environmental microbiology.     

Genome homology and divergence in the SP beta-related bacteriophages of Bacillus subtilis

Chromosomal organization in related temperate Bacillus subtilis bacteriophages SP beta, phi 3T, rho 11, Z, and E was compared. DNA-DNA hybridization studies done in conjunction with available restriction fragment maps of SP beta, phi 3T, and rho 11 demonstrated that DNA homology between these three phages extended over most of their respective genomes, although each contained unique chromosomal segments, phi 3T, rho 11, Z, and E, but not SP beta, possessed apparently homologous structural genes (thyP) for thymidylate synthetase. DNA from all thyP-containing phages transformed thymine auxotrophs of B. subtilis SP beta lysogens to prototrophy. This transformation commonly involved incorporation of the thyP gene into SP beta prophage within a region corresponding to the middle of the viral chromosome. Chimeric plasmids containing the thyP gene from phi 3T or cloned fragments of SP beta DNA were used in DNA-DNA hybridization studies to locate the thymidylate synthetase gene near the center of the phi 3T chromosome, and to demonstrate that the organization of this region resembled the analogous portion of the SP beta genome. Profiles of virion structural proteins from the five phages were also very similar, further suggesting functional homology between these viruses. However, despite these evidences of relatedness, populations of fragments generated by digesting SP beta, phi 3T, rho 11, Z, and E DNA with restriction enzymes were quite dissimilar.     

Cloning of human papilloma virus genomic DNAs and analysis of homologous polynucleotide sequences.

The complete DNA genomes of four distinct human papilloma viruses (human papilloma virus subtype 1a [HPV-1a], HPV-1b, HPV-2a, and HPV-4) were molecularly cloned in Escherichia coli, using the certified plasmid vector pBR322. The restriction endonuclease patterns of the cloned HPV-1a and HPV-1b DNAs were similar to those already published for uncloned DNAs. Physical maps were constructed for HPV-2a DNA and HPV-4 DNA, since these viral DNAs had not been previously mapped. By using the cloned DNAs, the genomes of HPV-1a, HPV-2a, and HPV-4 were analyzed for nucleotide sequence homology. Under standard hybridization conditions (Tm = --28 degrees C), no homology was detectable among the genomes of these papilloma viruses, in agreement with previous reports. However, under less stringent conditions (i.e., Tm = --50 degrees C), stable DNA hybrids could be detected between these viral DNAs, indicating homologous segments in the genomes with approximately 30% base mismatch. By using specific DNA fragments immobilized on nitrocellulose filters, these regions of homology were mapped. Hybridization experiments between radiolabeled bovine papilloma virus type 1 (BPV-1) DNA and the unlabeled HPV-1a, HPV-2a, or HPV-4 DNA restriction fragments under low-stringency conditions indicated that the regions of homology among the HPV DNAs are also conserved in the BPV-1 genome with approximately the same degree of base mismatch.     

Targeted isolation of a designated region of the Bacillus subtilis genome by recombinational transfer

A method for positional cloning of the Bacillus subtilis genome was developed. The method requires a set of two small DNA fragments that flank the region to be copied. A 38-kb segment that carries genes ppsABCDE encoding five enzymes for antibiotic plipastatin synthesis and another genome locus as large as 100 kb including one essential gene were examined for positional cloning. The positional cloning vector for ppsABCDE was constructed using a B. subtilis low-copy-number plasmid that faithfully copied the precise length of the 38-kb DNA in vivo via the recombinational transfer system of this bacterium. Structure of the copied DNA was confirmed by restriction enzyme analyses. Furthermore, the unaltered structure of the 38-kb DNA was demonstrated by complementation of a ppsABCDE deletion mutant.     

Repair and polyadenylation of a naturally occurring hepatitis C virus 3' nontranslated region-shorter variant in selectable replicon cell lines

The 3' nontranslated region (NTR) of the hepatitis C virus (HCV) genome is highly conserved and contains specific cis-acting RNA motifs that are essential in directing the viral replication machinery to initiate at the correct 3' end of the viral genome. Since the ends of viral genomes may be damaged by cellular RNases, preventing the initiation of viral RNA replication, stable RNA hairpin structures in the 3' NTR may also be essential in host defense against exoribonucleases. During 3'-terminal sequence analysis of serum samples of a patient with chronic hepatitis related to an HCV1b infection, a number of clones were obtained that were several nucleotides shorter at the extreme 3' end of the genome. These shorter 3' ends were engineered in selectable HCV replicons in order to enable the study of RNA replication in cell culture. When in vitro-transcribed subgenomic RNAs, containing shorter 3' ends, were introduced into Huh-7 cells, a few selectable colonies were obtained, and the 3' terminus of these subgenomic RNAs was sequenced. Interestingly, most genomes recovered from these colonies had regained the wild-type 3' ends, showing that HCV, like several other positive-stranded RNA viruses, has developed a strategy to repair deleted 3' end nucleotides. Furthermore, we found several genomes in these replicon colonies that contained a poly(A) tail and a short linker sequence preceding the poly(A) tail. After recloning and subsequent passage in Huh-7 cells, these poly(A) tails persisted and varied in length. In addition, the connecting linker became highly diverse in sequence and length, suggesting that these tails are actively replicated. The possible terminal repair mechanisms, including roles for the poly(A) tail addition, are discussed.