Smallpox holds a unique position in the history of medicine. It was the first disease for which a vaccine was developed and remains the only human disease eradicated by vaccination. Although there have been claims of smallpox in Egypt, India, and China dating back millennia [1-4], the timescale of emergence of the causative agent, variola virus (VARV), and how it evolved in the context of increasingly widespread immunization, have proven controversial [4-9]. In particular, some molecular-clock-based studies have suggested that key events in VARV evolution only occurred during the last two centuries [4-6] and hence in apparent conflict with anecdotal historical reports, although it is difficult to distinguish smallpox from other pustular rashes by description alone. To address these issues, we captured, sequenced, and reconstructed a draft genome of an ancient strain of VARV, sampled from a Lithuanian child mummy dating between 1643 and 1665 and close to the time of several documented European epidemics [1, 2, 10]. When compared to vaccinia virus, this archival strain contained the same pattern of gene degradation as 20th century VARVs, indicating that such loss of gene function had occurred before ca. 1650. Strikingly, the mummy sequence fell basal to all currently sequenced strains of VARV on phylogenetic trees. Molecular-clock analyses revealed a strong clock-like structure and that the timescale of smallpox evolution is more recent than often supposed, with the diversification of major viral lineages only occurring within the 18th and 19th centuries, concomitant with the development of modern vaccination.

Variola virus, the agent of smallpox, has a severely restricted host range (humans) but a devastatingly high mortality rate. Although smallpox has been eradicated by a World Health Organization vaccination program, knowledge of the evolutionary processes by which human super-pathogens such as variola virus arise is important. By analyzing the evolution of variola and other closely related poxviruses at the level of single nucleotide polymorphisms we detected a hotspot of genome variation within the smallpox ortholog of the vaccinia virus O1L gene, which is known to be necessary for efficient replication of vaccinia virus in human cells. These mutations in the variola virus ortholog and the subsequent loss of the functional gene from camelpox virus and taterapox virus, the two closest relatives of variola virus, strongly suggest that changes within this region of the genome may have played a key role in the switch to humans as a host for the ancestral virus and the subsequent host-range restriction that must have occurred to create the phenotype exhibited by smallpox.

We report a major improvement to the assembly of published short read sequencing data from an ancient variola virus (VARV) genome by the removal of contig-capping sequencing tags and manual searches for gap-spanning reads. The new assembly, together with camelpox and taterapox genomes, permitted new dates to be calculated for the last common ancestor of all VARV genomes. The analysis of recently sequenced VARV-like cowpox virus genomes showed that single nucleotide polymorphisms (SNPs) and amino acid changes in the vaccinia virus (VACV)-Cop-O1L ortholog, predicted to be associated with VARV host specificity and virulence, were introduced into the lineage before the divergence of these viruses. A comparison of the ancient and modern VARV genome sequences also revealed a measurable drift towards adenine + thymine (A + T) richness.

Due to the high degree of sequence identity between Orthopoxvirus species, the specific B and T cell responses raised against these viruses are largely cross-reactive and poorly selective. We therefore searched for CD4 T cell epitopes present in the conserved parts of the Vaccinia genome (VACV) but absent from Variola viruses (VARV), with a view to identifying immunogenic sequences selective for VACV. We identified three long peptide fragments from the B7R, B10R and E7R proteins by in silico comparisons of the poxvirus genomes, and evaluated the recognition of these fragments by VACV-specific T cell lines derived from healthy donors. For the 12 CD4 T cell epitopes identified, we assessed their binding to common HLA-DR allotypes and their capacity to induce peptide-specific CD4 T-cell lines. Four peptides from B7R and B10R displayed a broad binding specificity for HLA-DR molecules and induced multiple T cell lines from healthy donors. Besides their absence from VARV, the two B10R peptide sequences were mutated in the Cowpox virus and completely absent from the Monkeypox genome. This work contributes to the development of differential diagnosis of poxvirus infections.

Continuing the work developed by our research group, in the present manuscript, we performed a theoretical study of 10 new structures derived from the antivirals cidofovir and ribavirin, as inhibitor prototypes for the enzyme thymidylate kinase from Variola virus (VarTMPK). The proposed structures were subjected to docking calculations, molecular dynamics simulations, and free energy calculations, using the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method, inside the active sites of VarTMPK and human TMPK (HssTMPK). The docking and molecular dynamic studies pointed to structures 2, 3, 4, 6, and 9 as more selective towards VarTMPK. In addition, the free energy data calculated through the MM-PBSA method, corroborated these results. This suggests that these compounds are potential selective inhibitors of VarTMPK and, thus, can be considered as template molecules to be synthesized and experimentally evaluated against smallpox.

Two distinct phenomena of airborne transmission of variola virus (smallpox) were described in the pre-eradication era-direct respiratory transmission, and a unique phenomenon of transmission over greater distances, referred to as "aerial convection". We conducted an analysis of data obtained from a systematic review following the PRISMA criteria, on the long-distance transmission of smallpox. Of 8179 studies screened, 22 studies of 17 outbreaks were identified-12 had conclusive evidence of aerial convection and five had partially conclusive evidence. Aerial convection was first documented in 1881 in England, when smallpox incidence had waned substantially following mass vaccination, making unusual transmissions noticeable. National policy at the time stipulated spatial separation of smallpox hospitals from other buildings and communities. The evidence supports the transmission of smallpox through aerial convection at distances ranging from 0.5 to 1 mile, and one instance of 15 km related to bioweapons testing. Other explanations are also possible, such as missed chains of transmission, fomites or secondary aerosolization from contaminated material such as bedding. The window of observation of aerial convection was within the 100 years prior to eradication. Aerial convection appears unique to the variola virus and is not considered in current hospital infection control protocols. Understanding potential aerial convection of variola should be an important consideration in planning for smallpox treatment facilities and protecting potential contacts and surrounding communities.

Archaeovirology efforts provided a rich portrait of the evolutionary history of variola virus (VARV, the cause of smallpox), which was characterized by lineage extinctions and a relatively recent origin of the virus as a human pathogen (~1700 years ago, ya). This contrasts with historical records suggesting the presence of smallpox as early as 3500 ya. By performing an analysis of ancestry components in modern, historic, and ancient genomes, we unveil the progressive drifting of VARV lineages from a common ancestral population and we show that a small proportion of Viking Age ancestry persisted until the 18th century. After the split of the P-I and P-II lineages, the former experienced a severe bottleneck. With respect to the emergence of VARV as a human pathogen, we revise time estimates by accounting for the time-dependent rate phenomenon. We thus estimate that VARV emerged earlier than 3800 ya, supporting its presence in ancient societies, as pockmarked Egyptian mummies suggest.

Zoonotic infections caused by several orthopoxviruses (OPV) like monkeypox virus or vaccinia virus have a significant impact on human health. In Europe, the number of diagnosed infections with cowpox viruses (CPXV) is increasing in animals as well as in humans. CPXV used to be enzootic in cattle; however, such infections were not being diagnosed over the last decades. Instead, individual cases of cowpox are being found in cats or exotic zoo animals that transmit the infection to humans. Both animals and humans reveal local exanthema on arms and legs or on the face. Although cowpox is generally regarded as a self-limiting disease, immunosuppressed patients can develop a lethal systemic disease resembling smallpox. To date, only limited information on the complex and, compared to other OPV, sparsely conserved CPXV genomes is available. Since CPXV displays the widest host range of all OPV known, it seems important to comprehend the genetic repertoire of CPXV which in turn may help elucidate specific mechanisms of CPXV pathogenesis and origin. Therefore, 22 genomes of independent CPXV strains from clinical cases, involving ten humans, four rats, two cats, two jaguarundis, one beaver, one elephant, one marah and one mongoose, were sequenced by using massive parallel pyrosequencing. The extensive phylogenetic analysis showed that the CPXV strains sequenced clearly cluster into several distinct clades, some of which are closely related to Vaccinia viruses while others represent different clades in a CPXV cluster. Particularly one CPXV clade is more closely related to Camelpox virus, Taterapox virus and Variola virus than to any other known OPV. These results support and extend recent data from other groups who postulate that CPXV does not form a monophyletic clade and should be divided into multiple lineages.

Smallpox was one of the most devastating diseases in the human history and still represents a serious menace today due to its potential use by bioterrorists. Considering this threat and the non-existence of effective chemotherapy, we propose the enzyme thymidylate kinase from Variola virus (VarTMPK) as a potential target to the drug design against smallpox. We first built a homology model for VarTMPK and performed molecular docking studies on it in order to investigate the interactions with inhibitors of Vaccinia virus TMPK (VacTMPK). Subsequently, molecular dynamics (MD) simulations of these compounds inside VarTMPK and human TMPK (HssTMPK) were carried out in order to select the most promising and selective compounds as leads for the design of potential VarTMPK inhibitors. Results of the docking and MD simulations corroborated to each other, suggesting selectivity towards VarTMPK and, also, a good correlation with the experimental data.

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.

Subversion of host cell apoptosis is an important survival strategy for viruses to ensure their own proliferation and survival. Certain viruses express proteins homologous in sequence, structure and function to mammalian pro-survival B-cell lymphoma 2 (Bcl-2) proteins, which prevent rapid clearance of infected host cells. In vaccinia virus (VV), the virulence factor F1L was shown to be a potent inhibitor of apoptosis that functions primarily be engaging pro-apoptotic Bim. Variola virus (VAR), the causative agent of smallpox, harbors a homolog of F1L of unknown function. We show that VAR F1L is a potent inhibitor of apoptosis, and unlike all other characterized anti-apoptotic Bcl-2 family members lacks affinity for the Bim Bcl-2 homology 3 (BH3) domain. Instead, VAR F1L engages Bid BH3 as well as Bak and Bax BH3 domains. Unlike its VV homolog, variola F1L only protects against Bax-mediated apoptosis in cellular assays. Crystal structures of variola F1L bound to Bid and Bak BH3 domains reveal that variola F1L forms a domain-swapped Bcl-2 fold, which accommodates Bid and Bak BH3 in the canonical Bcl-2-binding groove, in a manner similar to VV F1L. Despite the observed conservation of structure and sequence, variola F1L inhibits apoptosis using a startlingly different mechanism compared with its VV counterpart. Our results suggest that unlike during VV infection, Bim neutralization may not be required during VAR infection. As molecular determinants for the human-specific tropism of VAR remain essentially unknown, identification of a different mechanism of action and utilization of host factors used by a VAR virulence factor compared with its VV homolog suggest that studying VAR directly may be essential to understand its unique tropism.

Responding to viral infection, the interferon-induced, double-stranded RNA (dsRNA)-activated protein kinase PKR phosphorylates translation initiation factor eIF2α to inhibit cellular and viral protein synthesis. To overcome this host defense mechanism, many poxviruses express the protein E3L, containing an N-terminal Z-DNA binding (Zα) domain and a C-terminal dsRNA-binding domain (dsRBD). While E3L is thought to inhibit PKR activation by sequestering dsRNA activators and by directly binding the kinase, the role of the Zα domain in PKR inhibition remains unclear. Here, we show that the E3L Zα domain is required to suppress the growth-inhibitory properties associated with expression of human PKR in yeast, to inhibit PKR kinase activity in vitro, and to reverse the inhibitory effects of PKR on reporter gene expression in mammalian cells treated with dsRNA. Whereas previous studies revealed that the Z-DNA binding activity of E3L is critical for viral pathogenesis, we identified point mutations in E3L that functionally uncouple Z-DNA binding and PKR inhibition. Thus, our studies reveal a molecular distinction between the nucleic acid binding and PKR inhibitory functions of the E3L Zα domain, and they support the notion that E3L contributes to viral pathogenesis by targeting PKR and other components of the cellular anti-viral defense pathway.

Concerns that variola viruses might be used as bioweapons have renewed the interest in developing new and safer smallpox vaccines. Variola virus genomes are now widely available, allowing computational characterization of the entire T-cell epitome and the use of such information to develop safe and yet effective vaccines. To this end, we identified 124 proteins shared between various species of pathogenic orthopoxviruses including variola minor and major, monkeypox, cowpox, and vaccinia viruses, and we targeted them for T-cell epitope prediction. We recognized 8,106, and 8,483 unique class I and class II MHC-restricted T-cell epitopes that are shared by all mentioned orthopoxviruses. Subsequently, we developed an immunological resource, EPIPOX, upon the predicted T-cell epitome. EPIPOX is freely available online and it has been designed to facilitate reverse vaccinology. Thus, EPIPOX includes key epitope-focused protein annotations: time point expression, presence of leader and transmembrane signals, and known location on outer membrane structures of the infective viruses. These features can be used to select specific T-cell epitopes suitable for experimental validation restricted by single MHC alleles, as combinations thereof, or by MHC supertypes.

Recently we constructed a homology model of the enzyme thymidylate kinase from Variola virus (VarTMPK) and proposed it as a new target to the drug design against smallpox. In the present work, we used the antivirals cidofovir and acyclovir as reference compounds to choose eleven compounds as leads to the drug design of inhibitors for VarTMPK. Docking and molecular dynamics (MD) studies of the interactions of these compounds inside VarTMPK and human TMPK (HssTMPK) suggest that they compete for the binding region of the substrate and were used to propose the structures of ten new inhibitors for VarTMPK. Further docking and MD simulations of these compounds, inside VarTMPK and HssTMPK, suggest that nine among ten are potential selective inhibitors of VarTMPK.

Histone deacetylases (HDACs) regulate gene expression and innate immunity. Previously, we showed that HDAC5 is degraded during Vaccinia virus (VACV) infection and is a restriction factor for VACV and herpes simplex virus type 1. Here, we report that HDAC5 promotes interferon regulatory factor 3 (IRF3) activation downstream of Toll-IL-1 receptor (TIR) domain-containing adaptor molecule-1 or Sendai virus-mediated stimulation without requiring HDAC activity. Loss of HDAC5-mediated IRF3 activation is restored by re-introduction of HDAC5 but not HDAC1 or HDAC4. The antiviral activity of HDAC5 is antagonized by VACV protein C6 and orthologs from the orthopoxviruses cowpox, rabbitpox, camelpox, monkeypox, and variola. Infection by many of these viruses induces proteasomal degradation of HDAC5, and expression of C6 alone can induce HDAC5 degradation. Mechanistically, C6 binds to the dimerization domain of HDAC5 and prevents homodimerization and heterodimerization with HDAC4. Overall, this study describes HDAC5 as a positive regulator of IRF3 activation and provides mechanistic insight into how the poxviral protein C6 binds to HDAC5 to antagonize its function.

After the recent summary of World Health Organization-authorized research on smallpox, several clinical issues remain. This policy review addresses whether early hemorrhagic smallpox is disseminated intravascular coagulation and speculates about the cause of the high mortality rate among pregnant women and whether ocular smallpox is partly the result of trachoma or vitamin A deficiency. The joint destruction common in children with smallpox might be prevented by antiviral drugs, but intraarticular infusion of antiviral drugs is unprecedented. Development of highly effective antiviral drugs against smallpox raises the issue of whether postexposure vaccination can be performed without interference by an antiviral drug. Clinicians should consider whether patients with smallpox should be admitted to general hospitals. Although an adequate supply of second-generation smallpox vaccine exists in the United States, its use is unclear. Finally, political and ethical forces suggest that destruction of the remaining stocks of live smallpox virus is now appropriate.

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'.

Acknowledging the importance of studies toward the development of measures against terrorism and bioterrorism, this study aims to contribute to the design of new prototypes of potential drugs against smallpox. Based on a former study, nine synthetic feasible prototypes of selective inhibitors for thymidylate kinase from Variola virus (VarTMPK) were designed and submitted to molecular docking, molecular dynamics simulations and binding energy calculations. The compounds are simplifications of two more complex scaffolds, with a guanine connected to an amide or alcohol through a spacer containing ether and/or amide groups, formerly suggested as promising for the design of selective inhibitors of VarTMPK. Our study showed that, despite the structural simplifications, the compounds presented effective energy values in interactions with VarTMPK and HssTMPK and that the guanine could be replaced by a simpler imidazole ring linked to a -NH2 group, without compromising the affinity for VarTMPK. It was also observed that a positive charge in the imidazole ring is important for the selectivity toward VarTMPK and that an amide group in the spacer does not contribute to selectivity. Finally, prototype 3 was pointed as the most promising to be synthesized and experimentally evaluated. Communicated by Ramaswamy H. Sarma.

Variola virus (VARV) the causative agent of smallpox, eradicated in 1980, have wide spectrum of immunomodulatory proteins to evade host immunity. Recently additional biological activity was discovered for VARV CrmB protein, known to bind and inhibit tumour necrosis factor (TNF) through its N-terminal domain homologous to cellular TNF receptors. Besides binding TNF, this protein was also shown to bind with high affinity several chemokines which recruit B- and T-lymphocytes and dendritic cells to sites of viral entry and replication. Ability to bind chemokines was shown to be associated with unique C-terminal domain of CrmB protein. This domain named SECRET (Smallpox virus-Encoded Chemokine Receptor) is unrelated to the host proteins and lacks significant homology with other known viral chemokine-binding proteins or any other known protein.

Numerous animal models of systemic orthopoxvirus disease have been developed to evaluate therapeutics against variola virus (VARV), the causative agent of smallpox. These animal models do not resemble the disease presentation in human smallpox and most used surrogate Orthopoxviruses. A rodent model using VARV has a multitude of advantages, and previous investigations identified the CAST/EiJ mouse as highly susceptible to monkeypox virus infection, making it of interest to determine if these rodents are also susceptible to VARV infection. In this study, we inoculated CAST/EiJ mice with a range of VARV doses (102-106 plaque forming units). Some animals had detectable viable VARV from the oropharynx between days 3 and 12 post inoculation. Despite evidence of disease, the CAST/EiJ mouse does not provide a model for clinical smallpox due to mild signs of morbidity and limited skin lesions. However, in contrast to previous rodent models using VARV challenge (i.e. prairie dogs and SCID mice), a robust immune response was observed in the CAST/EiJ mice (measured by Immunoglobulin G enzyme-linked immunosorbent assay). This is an advantage of this model for the study of VARV and presents a unique potential for the study of the immunomodulatory pathways following VARV infection.