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.     

Molecular evolution of poxviruses

Previous restriction fragment length polymorphism analysis divided variola virus (VARV) strains into two subtypes, one of which included West African and South American isolates. This allowed a dating to be introduced for the first time in estimation of the VARV evolution rate. The results were used to analyze the molecular evolution of the total family Poxviridae. Comparisons of the known nucleotide sequences were performed for the extended conserved central genome region in 42 orthopoxvirus strains and for the eight genes of multisubunit RNA polymerase in 65 viruses belonging to various genera of the family Poxviridae. Using the Bayesian dating method, the mutation accumulation rate of poxviruses was estimated at (1.7-8.8) x 10(-6) nucleotide substitutions per site per year. Computations showed that the modem poxvirus genera started diverging from an ancestral virus more than 200 thousand years ago and that an ancestor of the genus Orthopoxvirus emerged 131 +/- 45 thousand years ago. The other genera of mammalian poxviruses with a low GC content diverged approximately 110-90 thousand years ago. The independent evolution of VARV started 3.4 +/- 0.8 thousand years ago. It was shown with the example of VARV and the monkeypox virus (MPXV) that divergent evolution of these orthopoxviruses started and the West African subtypes of VARV and MPXV were formed as geographical conditions changed to allow isolation of West African animals from other African regions.     

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.     

Molecular evolution of poxviruses

Previous restriction fragment length polymorphism analysis divided variola virus (VARV) strains into two subtypes, one of which included West African and South American isolates. This allowed a dating to be introduced for the first time in estimation of the VARV evolution rate. The results were used to analyze the molecular evolution of the total family Poxviridae. Comparisons of the known nucleotide sequences were performed for the extended conserved central genome region in 42 orthopoxvirus strains and for the eight genes of multisubunit RNA polymerase in 65 viruses belonging to various genera of the family Poxviridae. Using the Bayesian dating method, the mutation accumulation rate of poxviruses was estimated at (1.7–8.8) × 10^?6 nucleotide substitutions per site per year. Computations showed that the modern poxvirus genera started diverging from an ancestral virus more than 200 thousand years ago and that an ancestor of the genus Orthopoxvirus emerged 131 ± 45 thousand years ago. The other genera of mammalian poxviruses with a low GC content diverged approximately 110–90 thousand years ago. The independent evolution of VARV started 3.4 ± 0.8 thousand years ago. It was shown with the example of VARV and the monkeypox virus (MPXV) that divergent evolution of these orthopoxviruses started and the West African subtypes of VARV and MPXV were formed as geographical conditions changed to allow isolation of West African animals from other African regions.     

Time scale of Poxvirus evolution

Unlike in vertebrates and RNA viruses, the molecular clock has not been estimated so far for DNA viruses. The extended conserved central region (102 kb) of the orthopoxvirus genome and the DNA polymerase gene (3 kb) were analyzed in viruses representing several genera of the family Poxviridae. Analysis was based on the known dating of the variola virus (VARV) transfer from Western Africa to South America and previous data on the phylogenetic relatedness of modern West African and South American isolates of VARV. The mutation accumulation rate was for the first time estimated for these DNA viruses at (0.9–1.2) × 10^?6 substitutions per site per year. It was assumed that poxviruses diverged from an ancestor approximately 500,000 years ago to form the recent species and that the ancestor of the genus Orthopoxvirus emerged approximately 300,000 years ago and gave origin to the modern species approximately 14,000 years ago.     

Comparative analysis of variable regions in the genomes of variola virus

Nucleotide sequences of two extended segments of the terminal variable regions in variola virus genome were determined. The size of the left segment was 13.5 kbp and of the right, 10.5 kbp. Totally, over 540 kbp were sequenced for 22 variola virus strains. The conducted phylogenetic analysis and the data published earlier allowed us to find the interrelations between 70 variola virus isolates, the character of their clustering, and the degree of intergroup and intragroup variations of the clusters of variola virus strains. The most polymorphic loci of the genome segments studied were determined. It was demonstrated that that these loci are localized to either noncoding genome regions or to the regions of destroyed open reading frames, characteristic of the ancestor virus. These loci are promising for development of the strategy for genotyping variola virus strains. Analysis of recombination using various methods demonstrated that, with the only exception, no statistically significant recombinational events in the genomes of variola virus strains studied were detectable.     

Variola virus genome sequenced from an eighteenth-century museum specimen supports the recent origin of smallpox

Smallpox, caused by the variola virus (VARV), was a highly virulent disease with high mortality rates causing a major threat for global human health until its successful eradication in 1980. Despite previously published historic and modern VARV genomes, its past dissemination and diversity remain debated. To understand the evolutionary history of VARV with respect to historic and modern VARV genetic variation in Europe, we sequenced a VARV genome from a well-described eighteenth-century case from England (specimen P328). In our phylogenetic analysis, the new genome falls between the modern strains and another historic strain from Lithuania, supporting previous claims of larger diversity in early modern Europe compared to the twentieth century. Our analyses also resolve a previous controversy regarding the common ancestor between modern and historic strains by confirming a later date around the seventeenth century. Overall, our results point to the benefit of historic genomes for better resolution of past VARV diversity and highlight the value of such historic genomes from around the world to further understand the evolutionary history of smallpox as well as related diseases.This article is part of the theme issue ‘Insights into health and disease from ancient biomolecules’.     

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.     

Pandemic potential of poxviruses: From an ancient killer causing smallpox to the surge of monkeypox

Smallpox caused by the variola virus (VARV) was one of the greatest infectious killers of mankind. Historical records trace back smallpox for at least a millennium while phylogenetic analysis dated the ancestor of VARV circulating in the 20th century into the 19th century. The discrepancy was solved by the detection of distinct VARV sequences first in 17th-century mummies and then in human skeletons dated to the 7th century. The historical records noted marked variability in VARV virulence which scientists tentatively associated with gene losses occurring when broad-host poxviruses narrow their host range to a single host. VARV split from camel and gerbil poxviruses and had no animal reservoir, a prerequisite for its eradication led by WHO. The search for residual pockets of VARV led to the discovery of the monkeypox virus (MPXV); followed by the detection of endemic smallpox-like monkeypox (mpox) disease in Africa. Mpox is caused by less virulent clade 2 MPXV in West Africa and more virulent clade 1 MPXV in Central Africa. Exported clade 2 mpox cases associated with the pet animal trade were observed in 2003 in the USA. In 2022 a world-wide mpox epidemic infecting more than 80,000 people was noted, peaking in August 2022 although waning rapidly. The cases displayed particular epidemiological characteristics affecting nearly exclusively young men having sex with men (MSM). In contrast, mpox in Africa mostly affects children by non-sexual transmission routes possibly from uncharacterized animal reservoirs. While African children show a classical smallpox picture, MSM mpox cases show few mostly anogenital lesions, low-hospitalization rates and 140 fatal cases worldwide. MPXV strains from North America and Europe are closely related, derived from clade 2 African MPXV. Distinct transmission mechanisms are more likely causes for the epidemiological and clinical differences between endemic African cases and the 2022 epidemic cases than viral traits.     

Properties of the recombinant TNF-binding proteins from variola, monkeypox, and cowpox viruses are different

Tumor necrosis factor (TNF), a potent proinflammatory and antiviral cytokine, is a critical extracellular immune regulator targeted by poxviruses through the activity of virus-encoded family of TNF-binding proteins (CrmB, CrmC, CrmD, and CrmE). The only TNF-binding protein from variola virus (VARV), the causative agent of smallpox, infecting exclusively humans, is CrmB. Here we have aligned the amino acid sequences of CrmB proteins from 10 VARV, 14 cowpox virus (CPXV), and 22 monkeypox virus (MPXV) strains. Sequence analyses demonstrated a high homology of these proteins. The regions homologous to cd00185 domain of the TNF receptor family, determining the specificity of ligand-receptor binding, were found in the sequences of CrmB proteins. In addition, a comparative analysis of the C-terminal SECRET domain sequences of CrmB proteins was performed. The differences in the amino acid sequences of these domains characteristic of each particular orthopoxvirus species were detected. It was assumed that the species-specific distinctions between the CrmB proteins might underlie the differences in these physicochemical and biological properties. The individual recombinant proteins VARV-CrmB, MPXV-CrmB, and CPXV-CrmB were synthesized in a baculovirus expression system in insect cells and isolated. Purified VARV-CrmB was detectable as a dimer with a molecular weight of 90 kDa, while MPXV- and CPXV-CrmBs, as monomers when fractioned by non-reducing SDS-PAGE. The CrmB proteins of VARV, MPXV, and CPXV differed in the efficiencies of inhibition of the cytotoxic effects of human, mouse, or rabbit TNFs in L929 mouse fibroblast cell line. Testing of CrmBs in the experimental model of LPS-induced shock using SPF BALB/c mice detected a pronounced protective effect of VARV-CrmB. Thus, our data demonstrated the difference in anti-TNF activities of VARV-, MPXV-, and CPXV-CrmBs and efficiency of VARV-CrmB rather than CPXV- or MPXV-CrmBs against LPS-induced mortality in mice.     

Neutralizing Monoclonal Antibodies Cross-React with Fusion Proteins Encoded by <Emphasis Type="Italic">129L</Emphasis> of the Ectromelia Virus and <Emphasis Type="Italic">A30L</Emphasis> of the Variola Virus

PCR fragments containing the fusion protein genes 129L of the ectromelia virus (EV) and A30L of the variola virus (VARV) were cloned in pQE32. The expression products, recombinant prA30L and pr129L, were isolated from Escherichia coli cell lysates by metal-chelate affinity chromatography. The recombinant proteins retained the capability of oligomerization, characteristic of their natural analogs. ELISA and immunoblotting were used to test 22 monoclonal antibodies (mAbs) to orthopoxviruses (19 mAbs to EV, 2 mAbs to the vaccinia virus (VACV), and 1 mAb to the cowpox virus (CPXV)) for interaction with prA30L, pr129L, and orthopoxviruses. Twelve species-specific epitopes were found in the EV fusion protein 129L and its recombinant analog. Ten cross-reacting epitopes were found in the EV, CPXV, and VACV fusion proteins. Of these, nine epitopes were present both in prA30L and in the VARV fusion protein. Five mAbs interacting with cross-reacting epitopes were capable of efficient neutralization of VACV; two of these mAbs neutralized VARV. It was demonstrated that there are species-specific epitopes in EV 129L and cross-reacting epitopes in the EV, VARV, CPXV, and VACV fusion proteins, including epitopes that induced synthesis of virus-neutralizing antibodies against VACV and VARV.     

Evidence for Persistence of Ectromelia Virus in Inbred Mice,Recrudescence Following Immunosuppression and Transmission to Na?ve Mice

Orthopoxviruses (OPV), including variola, vaccinia, monkeypox, cowpox and ectromelia viruses cause acute infections in their hosts. With the exception of variola virus (VARV), the etiological agent of smallpox, other OPV have been reported to persist in a variety of animal species following natural or experimental infection. Despite the implications and significance for the ecology and epidemiology of diseases these viruses cause, those reports have never been thoroughly investigated. We used the mouse pathogen ectromelia virus (ECTV), the agent of mousepox and a close relative of VARV to investigate virus persistence in inbred mice. We provide evidence that ECTV causes a persistent infection in some susceptible strains of mice in which low levels of virus genomes were detected in various tissues late in infection. The bone marrow (BM) and blood appeared to be key sites of persistence. Contemporaneous with virus persistence, antiviral CD8 T cell responses were demonstrable over the entire 25-week study period, with a change in the immunodominance hierarchy evident during the first 3 weeks. Some virus-encoded host response modifiers were found to modulate virus persistence whereas host genes encoded by the NKC and MHC class I reduced the potential for persistence. When susceptible strains of mice that had apparently recovered from infection were subjected to sustained immunosuppression with cyclophosphamide (CTX), animals succumbed to mousepox with high titers of infectious virus in various organs. CTX treated index mice transmitted virus to, and caused disease in, co-housed naïve mice. The most surprising but significant finding was that immunosuppression of disease-resistant C57BL/6 mice several weeks after recovery from primary infection generated high titers of virus in multiple tissues. Resistant mice showed no evidence of a persistent infection. This is the strongest evidence that ECTV can persist in inbred mice, regardless of their resistance status.     

Comparison of the Interferon γ-Binding Proteins of the Variola and Monkeypox Viruses

PCR was used to amplify DNA fragments containing the genes for interferon γ (IFNγ)-binding proteins (IFNγBPs) of the variola virus (VARV) and monkeypox virus (MPXV). The genes were expressed from baculovirus DNAs in Sf21 insect cells. The recombinant proteins were isolated from the culture medium by affinity chromatography. PAGE and Western blot analysis of the culture media and affinity-purified recombinant proteins showed that, in contrast to their cell analogs, the viral IFNγBPs form dimers in the absence of the ligand, human IFNγ. The biological activity of recombinant IFNγBPs was inferred from suppression of the protective effect of human IFNγ on L68 cells infected with the mouse encephalomyocarditis virus. The viral proteins showed a dose-dependent IFNγ-neutralizing effect. The prospects of using IFNγBPs as IFNγ antagonists are discussed.     

Teaching a new mouse old tricks: Humanized mice as an infection model for Variola virus

Smallpox, caused by the solely human pathogen Variola virus (VARV), was declared eradicated in 1980. While known VARV stocks are secure, smallpox remains a bioterrorist threat agent. Recent U.S. Food and Drug Administration approval of the first smallpox anti-viral (tecovirimat) therapeutic was a successful step forward in smallpox preparedness; however, orthopoxviruses can become resistant to treatment, suggesting a multi-therapeutic approach is necessary. Animal models are required for testing medical countermeasures (MCMs) and ideally MCMs are tested directly against the pathogen of interest. Since VARV only infects humans, a representative animal model for testing therapeutics directly against VARV remains a challenge. Here we show that three different humanized mice strains are highly susceptible to VARV infection, establishing the first small animal model using VARV. In comparison, the non-humanized, immunosuppressed background mouse was not susceptible to systemic VARV infection. Following an intranasal VARV challenge that mimics the natural route for human smallpox transmission, the virus spread systemically within the humanized mouse before mortality (~ 13 days post infection), similar to the time from exposure to symptom onset for ordinary human smallpox. Our identification of a permissive/representative VARV animal model can facilitate testing of MCMs in a manner consistent with their intended use.     

The gene therapy of collagen-induced arthritis in rats by intramuscular administration of the plasmid encoding TNF-binding domain of variola virus CrmB protein

Wistar rats with collagen-induced arthritis were intramuscularly injected with the recombinant plasmid pcDNA/sTNF-BD encoding the sequence of the TNF-binding protein domain of variola virus CrmB protein (VARV sTNF-BD) or the pcDNA3.1 vector. Quantitative analysis showed that the histopathological changes in the hind-limb joints of rats were most severe in the animals injected with pcDNA3.1 and much less severe in the group of rats injected with pcDNA/sTNF-BD, which indicates that gene therapy of rheumatoid arthritis is promising in the case of local administration of plasmids governing the synthesis of VARV immunomodulatory proteins.     

Poxvirus-encoded gamma interferon binding protein dampens the host immune response to infection

Ectromelia virus (ECTV), a natural mouse pathogen and the causative agent of mousepox, is closely related to variola virus (VARV), which causes smallpox in humans. Mousepox is an excellent surrogate small-animal model for smallpox. Both ECTV and VARV encode a multitude of host response modifiers that target components of the immune system and that are thought to contribute to the high mortality rates associated with infection. Like VARV, ECTV encodes a protein homologous to the ectodomain of the host gamma interferon (IFN-gamma) receptor 1. We generated an IFN-gamma binding protein (IFN-gammabp) deletion mutant of ECTV to study the role of viral IFN-gammabp (vIFN-gammabp) in host-virus interaction and also to elucidate the contribution of this molecule to the outcome of infection. Our data show that the absence of vIFN-gammabp does not affect virus replication per se but does have a profound effect on virus replication and pathogenesis in mice. BALB/c mice, which are normally susceptible to infection with ECTV, were able to control replication of the mutant virus and survive infection. Absence of vIFN-gammabp from ECTV allowed the generation of an effective host immune response that was otherwise diminished by this viral protein. Mice infected with a vIFN-gammabp deletion mutant virus, designated ECTV-IFN-gammabp(Delta), produced increased levels of IFN-gamma and generated robust cell-mediated and antibody responses. Using several strains of mice that exhibit differential degrees of resistance to mousepox, we show that recovery or death from ECTV infection is determined by a balance between the host's ability to produce IFN-gamma and the virus' ability to dampen its effects.     

Molecular genetic methods of typing of the Bartonellae

The primer systems for the PCR detection of four house-keeping genes of bartonellae in clinical material were developed and tested. The tactics of the species RFLP typing was also developed and tested. The scheme of the species RFLP typing of bartonellae was tested using as an example two strains for the first time isolated in Russia from patients with endocarditis and fever of uncertain origin. The results of the typing were supported by sequencing of the amplicons obtained. According to the sequencing the isolates were attributed to the sub species Bartonella vinsonii, subsp. arupensis. The necessity of molecular epidemiological analysis of bartonelloses in Russia was substantiated.     

Determination of clonality and relatedness of Vibrio cholerae isolates by genomic fingerprinting, using long-range repetitive element sequence-based PCR

A high-throughput method which is applicable for rapid screening, identification, and delineation of isolates of Vibrio cholerae, sensitive to genome variation, and capable of providing phylogenetic inferences enhances environmental monitoring of this bacterium. We have developed and optimized a method for genomic fingerprinting of V. cholerae based on long-range PCR. The method uses a primer set directed to enterobacterial repetitive intergenic consensus sequences, a high-fidelity DNA polymerase, and analysis via conventional agarose gel electrophoresis. Long ( approximately 10 kb), highly reproducible amplicons were generated from V. cholerae isolates, including those from different geographical locations and historical strains isolated during the period 1931-2000. The amplicons yielded reduced variability in their densitometric band patterns to /=90% similarity, discriminating O serotypes and biotypes (classical versus El Tor) as well as pathogenic and nonpathogenic strains. Compared to genome similarity measured by DNA-DNA hybridization, the results showed good correlation (r = 0.7; P < 0.001), with five times less measurement error and without bias. The method permits both phylogenetic inference and clonal differentiation of individual V. cholerae strains, enables robust, high-throughput analysis, and does not require specialized equipment to perform. With access to a curated public database furnished with appropriate analytical software applications, the method should prove useful in large-scale multilaboratory surveys, especially those designed to detect specific pathogens in the natural environment.     

A Multiplex PCR/LDR Assay for the Simultaneous Identification of Category A Infectious Pathogens: Agents of Viral Hemorrhagic Fever and Variola Virus

CDC designated category A infectious agents pose a major risk to national security and require special action for public health preparedness. They include viruses that cause viral hemorrhagic fever (VHF) syndrome as well as variola virus, the agent of smallpox. VHF is characterized by hemorrhage and fever with multi-organ failure leading to high morbidity and mortality. Smallpox, a prior scourge, has been eradicated for decades, making it a particularly serious threat if released nefariously in the essentially non-immune world population. Early detection of the causative agents, and the ability to distinguish them from other pathogens, is essential to contain outbreaks, implement proper control measures, and prevent morbidity and mortality. We have developed a multiplex detection assay that uses several species-specific PCR primers to generate amplicons from multiple pathogens; these are then targeted in a ligase detection reaction (LDR). The resultant fluorescently-labeled ligation products are detected on a universal array enabling simultaneous identification of the pathogens. The assay was evaluated on 32 different isolates associated with VHF (ebolavirus, marburgvirus, Crimean Congo hemorrhagic fever virus, Lassa fever virus, Rift Valley fever virus, Dengue virus, and Yellow fever virus) as well as variola virus and vaccinia virus (the agent of smallpox and its vaccine strain, respectively). The assay was able to detect all viruses tested, including 8 sequences representative of different variola virus strains from the CDC repository. It does not cross react with other emerging zoonoses such as monkeypox virus or cowpox virus, or six flaviviruses tested (St. Louis encephalitis virus, Murray Valley encephalitis virus, Powassan virus, Tick-borne encephalitis virus, West Nile virus and Japanese encephalitis virus).     

Application of Multi Locus Sequence Typing to the analysis of the biodiversity of indigenous Saccharomyces cerevisiae wine yeasts from Lebanon

AIMS: To assess suitability of Multi Locus Sequence Typing (MLST) for investigating the biodiversity of wine yeast strains. This method was compared with established ones like microsatellite analysis or amplification of genomic regions flanked by repeated (delta) elements. METHODS AND RESULTS: DNA fragments were amplified and sequenced for 26 loci representing housekeeping genes, open reading frames (ORFs) of unknown functions or intergenic regions. A set of seven loci was tested on 84 Saccharomyces cerevisiae strains, including 65 strains isolated from traditional wineries in Lebanon, commercial wine strains and Asian isolates. An overall sequence diversity of 2.05% was observed, consisting of single nucleotide polymorphisms, 60% of them occurring in a heterozygous state. The number of polymorphic sites per locus varied between 4 and 14. The same set of strains was analysed by microsatellite typing on six polymorphic loci and by interdelta amplification. CONCLUSIONS: Clustering of MLST profiles clearly differentiated the Asian group of strains from Lebanese and European commercial strains that appear closely related. The current MLST scheme appears less discriminatory (92.27%) on closely related wine yeasts than microsatellite or interdelta typing (>99%). SIGNIFICANCE AND IMPACT OF THE STUDY: MLST is a highly reliable method for relatedness inference and promising for wine yeast typing.