Comparative analysis of variable regions in the variola virus genome

Two segments of the variable terminal regions of the variola virus (VARV) genome were sequenced in 22 strains from the Russian collection, including about 13.5 kb of the left segment and about 10.5 kb of the right segment. The total length of the sequences was over 540 kb. Phylogenetic analysis of the new and published data determined the relationships among 70 VARV strains, the character of their clustering, and the intergroup and intragroup variation of the strain clusters. Loci with the highest polymorphism rate were identified and proved to map to noncoding regions or to regions of damaged open reading frames, characteristic of the ancestral virus. These loci offer attractive possibilities for developing a strategy of VARV strain genotyping. Recombination analysis by different methods did not detect, except for a single case, significant recombination events in the VARV strains examined.     

Comparative restriction enzyme analysis of the genome in variola virus strains from the Russian collection

Comparative RFLP analysis was for the first time performed for 21 variola virus (VARV) strains of the Russian collection with 20 amplicons covering the total VARV genome. The amplicons were synthesized in the long polymerase chain reaction. A database useful as a reference for identifying VARV strains was generated. VARV strains isolated in different geographical regions were compared and proved to vary mostly in variable genome regions. Each of the dendrograms constructed included three clusters of African, Asian, and VARV-alastrim isolates. The VARV-alastrim isolates differed to the greatest extent from the other strains. VARV strains isolated during an ecdemic variola burst in Moscow (1960) grouped with Asian isolates. Polymorphism of VARV strains was for the first time observed for a single variola burst with a few affected patients.     

Comparative Restriction Enzyme Analysis of the Genome in Variola Virus Strains from the Russian Collection

Comparative RFLP analysis was for the first time performed for 21 variola virus (VARV) strains of the Russian collection with 20 amplicons covering the total VARV genome. The amplicons were synthesized in the long polymerase chain reaction. A database useful as a reference for identifying VARV strains was generated. VARV strains isolated in different geographical regions were compared and proved to vary mostly in variable genome regions. Each of the dendrograms constructed included three clusters of African, Asian, and VARV-alastrim isolates. The VARV-alastrim isolates differed to the greatest extent from the other strains. VARV strains isolated during an ecdemic variola burst in Moscow (1960) grouped with Asian isolates. Polymorphism of VARV strains was for the first time observed for a single variola burst with a few affected patients.     

The structure of the human c-fes/fps proto-oncogene.

We have determined the complete nucleotide sequence of a human DNA fragment of approximately 13 kbp, which was shown by Southern blot analysis to contain the entire v-fes/fps cellular homolog. The v-fes/fps homologous sequences were dispersed over 11 kbp in 18 interspersed segments which were flanked by splice junctions. Fusion of these segments created a DNA fragment in which coding regions similar to those observed in the viral oncogenes v-fes of the Gardner-Arnstein (GA) and Snyder-Theilen (ST) strains of feline sarcoma virus and v-fps found in Fujinami sarcoma virus could be identified. A potential initiation site in the first exon was found. About 200 nucleotides downstream of a translational stop codon in the v-fes/fps homologous region, a poly(A) addition signal was identified. The deduced amino acid sequence has a molecular weight of 93 390 dalton resembling NCP92, the recently described human c-fes/fps product. The topography of human c-fes/fps appeared to resemble that of chicken c-fps.     

Controlling the false-positive rate in multilocus genome scans for selection

Rapid typing of genetic variation at many regions of the genome is an efficient way to survey variability in natural populations in an effort to identify segments of the genome that have experienced recent natural selection. Following such a genome scan, individual regions may be chosen for further sequencing and a more detailed analysis of patterns of variability, often to perform a parametric test for selection and to estimate the strength of a recent selective sweep. We show here that not accounting for the ascertainment of loci in such analyses leads to false inference of natural selection when the true model is selective neutrality, because the procedure of choosing unusual loci (in comparison to the rest of the genome-scan data) selects regions of the genome with genealogies similar to those expected under models of recent directional selection. We describe a simple and efficient correction for this ascertainment bias, which restores the false-positive rate to near-nominal levels. For the parameters considered here, we find that obtaining a test with the expected distribution of P-values depends on accurately accounting both for ascertainment of regions and for demography. Finally, we use simulations to explore the utility of relying on outlier loci to detect recent selective sweeps. We find that measures of diversity and of population differentiation are more effective than summaries of the site-frequency spectrum and that sequencing larger regions (2.5 kbp) in genome-scan studies leads to more power to detect recent selective sweeps.     

Genomic maps of some strains within the Mycoplasma mycoides cluster

Genomic restriction maps for the small colony (SC) strains (PG1, KH3J, Gladysdale, and V5) of Mycoplasma mycoides subsp. mycoides (the agent of contagious bovine pleuropneumonia) and for Mycoplasma strain PG50 (classified as bovine serogroup 7), with respective sizes of 1,280, 1,280, 1,260, 1,230, and 1,040 kbp, were compared with the map (1,200 kbp) for a large colony strain (Y goat) of M. mycoides subsp. mycoides. The number and order of all mapped restriction sites were fully conserved in the SC genomes, as were the approximate positions of mapped loci. A number of these restriction sites in the Y genome and some, but fewer, in the PG50 genome appeared to be conserved. The SC and large colony strains shared conservation in the relative positions of the mapped loci, except for rpoC.     

Human cytomegalovirus clinical isolates carry at least 19 genes not found in laboratory strains.

Nucleotide sequence comparisons were performed on a highly heterogeneous region of three human cytomegalovirus strains, Toledo, Towne, and AD169. The low-passage, virulent Toledo genome contained a DNA segment of approximately 13 kbp that was not found in the Towne genome and a segment of approximately 15 kbp that was not found in the AD169 genome. The Towne strain contained approximately 4.7 kbp of DNA that was absent from the AD169 genome, and only about half of this segment was present, arranged in an inverted orientation, in the Toledo genome. These additional sequences were located at the unique long (UL)/b' (IRL) boundary within the L component of the viral genome. A region representing nucleotides 175082 to 178221 of the AD169 genome was conserved in all three strains; however, substantial reduction in the size of the adjacent b' sequence was found. The additional DNA segment within the Toledo genome contained 19 open reading frames not present in the AD169 genome. The additional DNA segment within the Towne genome contained four new open reading frames, only one of which shared homology with the Toledo genome. This comparison was extended to five additional clinical isolates, and the additional Toledo sequence was conserved in all. These findings reveal a dramatic level of genome sequence complexity that may explain the differences that these strains exhibit in virulence and tissue tropism. Although the additional sequences have not altered the predicted size of the viral genome (230 to 235 kbp), a total of 22 new open reading frames (denoted UL133 to UL154), many of which have sequence characteristics of glycoproteins, are now defined as cytomegalovirus specific. Our work suggests that wild-type virus carries more than 220 genes, some of which are lost by large-scale deletion and rearrangement of the UL/b' region during laboratory passage.     

Genetic analysis of endogenous xenotropic murine leukemia viruses: association with two common mouse mutations and the viral restriction locus Fv-1.

We have defined 40 endogenous xenotropic virus (Xmv) loci from several common inbred strains of mice by examining provirus-cell DNA junction fragments in recombinant inbred mice. Some inbred strains carried unique proviruses, but most Xmv loci were present in several strains, indicating that many Xmv integration events preexisted modern inbreeding. It was also clear that most Xmv junction fragment variation between inbred strains resulted from independent integration events and not modification or restriction site polymorphism following integration. Chromosomal assignments were determined for 32 Xmv loci by comparing their recombinant inbred strain distribution patterns to those of known genetic markers. The Xmv loci were generally dispersed throughout the genome, but several chromosomal regions contained more than one provirus. Furthermore, several close genetic associations with cellular genes were discovered. Four Xmv loci were closely linked to Fv-1b, a dominant viral resistance gene present in C57BL/6J, BALB/cJ, A/J, and several other strains. Xmv-28 was closely linked to rd (retinal degeneration), and Xmv-10 was closely linked to a (non-agouti), both of which are old mutations as inferred from their broad distribution in mice. We suggest that Xmv integration contributed to genetic diversity in the past and that much of this diversity exists today in common laboratory strains.     

Serial assembly of Thermus megaplasmid DNA in the genome of Bacillus subtilis 168: a BAC-based domino method applied to DNA with a high GC content

Bacillus subtilis is the only bacterium-based host able to clone giant DNA above 1000 kbp. DNA previously handled by this host was limited to that with GC content similar to or lower than that of the B. subtilis genome. To expand the target DNA range to higher GC content, we tried to clone a pTT27 megaplasmid (257 kbp, 69% of G+C) from Thermus thermophilus. To facilitate the reconstruction process, we subcloned pTT27 in a bacterial artificial chromosome (BAC) vector of Escherichia coli. Owing to the ability of BAC to carry around 100 kbp DNA, only 4 clones were needed to cover the pTT27 and conduct step-by-step assembly in the B. subtilis genome. The full length of 257 kbp was reconstructed through 3 intermediary lengths (108, 153, and 226 kbp), despite an unexpected difficulty in the maintenance of DNA >200 kbp. Retrieval of these four pTT27 segments from the B. subtilis genome by genetic transfer to a plasmid pLS20 was attempted. A stable plasmid clone was obtained only for the 108 and 153 kbp intermediates. The B. subtilis genome was demonstrated to accommodate large DNA with a high GC content, but may be restricted to less than 200 kbp by unidentified mechanisms.     

Pulsed-field gel electrophoresis-fingerprinting, genome size estimation and rrn loci number of Rhizobium galegae

I. HUBER AND S. SELENSKA-POBELL. 1994. The genomes of several Rhizobium galegae (of.) strains, which effectively nodulate Galega officinalis host, were analysed by pulsed-field gel electrophoresis (PFGE). Individual PFGE-fingerprints were obtained for every particular strain when the rarely cutting restriction endonucleases Spe I and Asn I were applied. In hybridization experiments, where a DNA fragment carrying the rrnB ribosomal RNA operon of Escherichia coli was used as a probe, the number of the resulting strain-specific Spe I and Asn I bands was reduced to three for all of the strains studied. This suggests that in Rh. galegae (of.) there are at least three rrn loci. On the basis of the lengths of the Spe I fragments, separated by PFGE, the genome size of five Rh. galegae (of.) strains was estimated to be 5852 ± 198 kbp.     

Unidirectional gene conversions in the chloroplast of Chlamydomonas interspecific hybrids

Summary With the goal of studying directly the inheritance and recombination of physically mapped markers on the chloroplast genome, we have recently identified and localized physical differences between the chloroplast DNAs (cpDNAs) of the interfertile algae Chlamydomonas eugametos and C. moewusii. Here we report the inheritance patterns of 24 polymorphic loci mapping throughout the chloroplast genome in hybrids recovered from reciprocal crosses between the two algae. Most polymorphic loci were found to be inherited mainly from the mt ⁺ parent, with no apparent preference for one or the other parental alternatives in reciprocal crosses. Virtually all hybrids, however, inherited exclusively the long alleles of three loci; i.e. an intron in the large subunit ribosomal RNA gene of C. eugametos, a 21 kbp sequence addition in the inverted repeat of the C. moewusii cpDNA and a 5.8 kbp sequence addition in one of the single-copy regions of C. moewusii cpDNA. As these alleles are derived from opposite parental strains, their unidirectional inheritance in hybrids results necessarily from interspecific recombination of cpDNA molecules. We propose that gene conversion events led to the spreading of the long alleles of the three loci.     

Polymorphism of Avian Leukosis Virus Subgroup E Loci Showing Selective Footprints in Chicken

Avian leukosis virus subgroup E (ALVE) is a family of endogenous retroviruses in the chicken genome. To investigate the genetic consequences of chicken domestication, we analyzed 18 ALVE loci in red jungle fowls, layers, broilers, and Chinese indigenous chickens. None of the ALVE loci tested were found in red jungle fowls, but 12 were present in domestic chickens. ALVE1 and ALVE16 are found in regions of the genome that harbor quantitative trait loci (QTL) affecting egg production traits. ALVE1 was fixed and ALVE16 was detected only in layers. By contrast, ALVE-b1, ALVE-b5, ALVE-b6, and ALVE-b8 integrated into regions of the genome that harbor QTL affecting meat production traits. Carrier frequencies of these four ALVE loci were high in broilers and low in Chinese local chickens; the loci were not found in the layers. This study demonstrated that insertionally polymorphic ALVE loci can illustrate the selective footprints in the chicken genome.     

Diagnosis of rotavirus infection with cloned cDNA copies of viral genome segments.

The diagnostic potential of cloned cDNA copies of human rotavirus (strain WA) genome segments for the detection of rotavirus in clinical specimens has been determined. A hybridization assay in which a mixture of 32P-labeled cDNAs representing the 11 rotavirus segments was used as a probe compared favorably with three frequently used diagnostic tests for rotavirus in terms of both specificity and sensitivity. Significantly, clinical isolates could be readily distinguished when cloned cDNA copies of individual genome segments were used independently as a probe. In assays in which genome RNA from rotaviruses of known subgroups and serotypes were tested, cloned probes that encode nonstructural viral proteins hybridized efficiently to genome RNAs of all strains, whereas cloned probes corresponding to genome segments 6 and 9 exhibited the potential for differentiating strains of different subgroups and serotypes. Cloned cDNA copies of rotavirus genome segments therefore offer considerable potential for improved general diagnosis of rotavirus in clinical specimens, as well as for epidemiological studies in which virus isolates can be distinguished on the basis of nucleotide sequence homology of individual genome segments.     

Isolation of a repeated DNA sequence from Bordetella pertussis

A repeated DNA sequence in the genome of Bordetella pertussis has been demonstrated. At least 20 copies of this sequence could be observed in either BamHI or EcoRI restriction enzyme digests of chromosomal DNA; fragments carrying the repeated DNA sequence ranged in size from about 1.5 to 20 kbp. The repeated DNA sequence was cloned from two separate regions of the B. pertussis genome, as shown by restriction enzyme site maps of the two clones and by hybridization studies. A small number of differences in the pattern of hybridization of the repeated DNA sequence to chromosomal DNA from several strains of B. pertussis was observed. No repeated DNA sequences were observed in one strain each of B. parapertussis and B. bronchiseptica, and there was no hybridization of B. pertussis DNA to Escherichia coli chromosomal DNA. The repeated DNA sequence was subcloned on a 2.54 kbp BamHI fragment from one of the two original clones. Restriction enzyme digests and hybridization studies showed that the repeated DNA sequence was about 1 kbp in size and had a single, internal ClaI site.     

Sequence organization of the rat genome by electron microscopy.

The size and arrangement of repetitive and inverted repeat (foldback) sequences in rat DNA were studied by visualization of hybrid and heteroduplex structures in the electron microscope. The self-reassociation of repetitive sequence-bearing DNA strands often results in the formation of four-ended "H" structures, whose duplex regions equal the repetitive sequence length and can be measured in the electron microscope. In this way, it was determined that the average size of the class of numerous short repetitive sequences is 0.40 +/- 0.15 kbp. Heteroduplex structures were prepared between long whole DNA single strands and short repeat-sequence-bearing strands. The analysis of these structures confirms that the size of the repetitive sequences in 0.4 kbp on average. Length measurements between adjacent duplexes show that the average spacing between two interspersed repeats is at least 1.5-1.8 kbp. By examining 29.4-kbp single strands after brief renaturation, the size and distribution of foldback sequences were determined. There are 1.9 X 10(5) foldback apirs per rat genome, spaced an average of 9.7 kbp apart according to our measurement. Repetitive, inverted repeat and unique sequences are interspersed with each other in at least half the genome.     

Genomic Sequence and Virulence of Clonal Isolates of Vaccinia Virus Tiantan,the Chinese Smallpox Vaccine Strain

Despite the worldwide eradication of smallpox in 1979, the potential bioterrorism threat from variola virus and the ongoing use of vaccinia virus (VACV) as a vector for vaccine development argue for continued research on VACV. In China, the VACV Tiantan strain (TT) was used in the smallpox eradication campaign. Its progeny strain is currently being used to develop a human immunodeficiency virus (HIV) vaccine. Here we sequenced the full genomes of five TT clones isolated by plaque purification from the TT (752-1) viral stock. Phylogenetic analysis with other commonly used VACV strains showed that TT (752-1) and its clones clustered and exhibited higher sequence diversity than that found in Dryvax clones. The ∼190 kbp genomes of TT appeared to encode 273 open reading frames (ORFs). ORFs located in the middle of the genome were more conserved than those located at the two termini, where many virulence and immunomodulation associated genes reside. Several patterns of nucleotide changes including point mutations, insertions and deletions were identified. The polymorphisms in seven virulence-associated proteins and six immunomodulation-related proteins were analyzed. We also investigated the neuro- and skin- virulence of TT clones in mice and rabbits, respectively. The TT clones exhibited significantly less virulence than the New York City Board of Health (NYCBH) strain, as evidenced by less extensive weight loss and morbidity in mice as well as produced smaller skin lesions and lower incidence of putrescence in rabbits. The complete genome sequences, ORF annotations, and phenotypic diversity yielded from this study aid our understanding of the Chinese historic TT strain and are useful for HIV vaccine projects employing TT as a vector.     

Structural organization and biological activity of molecular clones of the integrated genome of a BALB/c mouse sarcoma virus.

BALB/c mouse sarcoma virus (BALB-MSV) is a spontaneously occurring transforming retrovirus of mouse origin. The integrated form of the viral genome was cloned from the DNA of a BALB-MSV-transformed nonproducer NRK cell line in the Charon 9 strain of bacteriophage lambda. In transfection assays, the 19-kilobase-pair (kbp) recombinant DNA clone transformed NIH/3T3 mouse cells with an efficiency of 3 X 10(4) focus-forming units per pmol. Such transformants possessed typical BALB-MSV morphology and released BALB-MSV after helper virus superinfection. A 6.8-kbp DNA segment within the 19-kbp DNA possessed restriction enzyme sites identical to those of the linear BALB-MSV genome. Long terminal repeats of approximately 0.6 kbp were localized at either end of the viral genome by the presence of a repeated constellation of restriction sites and by hybridization of segments containing these sites with nick-translated Moloney murine leukemia virus long terminal repeat DNA. A continuous segment of at least 0.6 and no more than 0.9 kbp of helper virus-unrelated sequences was localized toward the 3' end of the viral genome in relation to viral RNA. A probe composed of these sequences detected six EcoRI-generated DNA bands in normal mouse cell DNA as well as a smaller number of bands in rat and human DNAs. These studies demonstrate that BALB-MSV, like previously characterized avian and mammalian transforming retroviruses, arose by recombination of a type C helper virus with a well-conserved cellular gene.     

The time scale in poxvirus evolution

Unlike vertebrates and RNA-containing viruses, the objective estimate of molecular clock for DNA-containing viruses was so far absent. An extended central conservative genomic region of orthopoxviruses (about 102 kbp) and the sequence of DNA polymerase gene (about 3 kbp) of the viruses belonging to various genera from the family Poxviridae were analyzed. During this analysis, the known dating of variola virus (VARV) transfer from West Africa to South America (XVI century) and our own data on close phylogenetic relations between the modem West African and South American VARV isolates were used. As a result of this work, it was calculated for the first time that the rate of mutation accumulation in these DNA-containing viruses amounted to 0.9-1.2 x 10(-6) substitutions per site per year. The poxviruses started separating from the ancestor virus to form the modem genera approximately 500 thousand years ago; the ancestor of the genus Orthopoxvirus separated about 300 thousand years ago; and its division into the modem studied species took place approximately 14 thousand years ago.     

Variola virus immune evasion proteins

Variola virus, the causative agent of smallpox, encodes approximately 200 proteins. Over 80 of these proteins are located in the terminal regions of the genome, where proteins associated with host immune evasion are encoded. To date, only two variola proteins have been characterized. Both are located in the terminal regions and demonstrate immunoregulatory functions. One protein, the smallpox inhibitor of complement enzymes (SPICE), is homologous to a vaccinia virus virulence factor, the vaccinia virus complement-control protein (VCP), which has been found experimentally to be expressed early in the course of vaccinia infection. Both SPICE and VCP are similar in structure and function to the family of mammalian complement regulatory proteins, which function to prevent inadvertent injury to adjacent cells and tissues during complement activation. The second variola protein is the variola virus high-affinity secreted chemokine-binding protein type II (CKBP-II, CBP-II, vCCI), which binds CC-chemokine receptors. The vaccinia homologue of CKBP-II is secreted both early and late in infection. CKBP-II proteins are highly conserved among orthopoxviruses, sharing approximately 85% homology, but are absent in eukaryotes. This characteristic sets it apart from other known virulence factors in orthopoxviruses, which share sequence homology with known mammalian immune regulatory gene products. Future studies of additional variola proteins may help illuminate factors associated with its virulence, pathogenesis and strict human tropism. In addition, these studies may also assist in the development of targeted therapies for the treatment of both smallpox and human immune-related diseases.     

Multiple double-stranded RNA segments are associated with virus particles infecting Trichomonas vaginalis.

Previous studies demonstrated that some isolates of the sexually transmitted protozoan Trichomonas vaginalis are infected with a nonsegmented, double-stranded RNA (dsRNA) virus. A reexamination of the total dsRNA extracted from several virus-harboring isolates indicated the presence of at least three dsRNAs with sizes ranging from 4.8 to 4.3 kbp. The double-stranded nature of each of the three segments was determined by hybridization experiments using riboprobes of opposite polarities obtained from cDNA generated to each of the segments. All three segments were present in agar clones originating from single organisms of T. vaginalis isolates, suggesting that the three segments were not the result of a mixed population of trichomonads harboring different sizes of dsRNA. The three segments were associated with CsCl-purified virus particles, as evidenced by electron microscopy, and RNAse treatment of the preparation containing virus particles did not destroy the dsRNAs. Finally, the individual dsRNA segments were purified for use as probes to determine whether the three dsRNAs shared any sequence homology. Each end-labeled dsRNA segment did not cross-hybridize to any of the other two segments, a finding consistent with the hybridization of labeled cDNAs to only the segments from which they were derived. These results show that the coding capacity of the dsRNA virus may be at least three times greater than that estimated earlier and illustrates further the complexity of this virus-parasite interrelationship.