Search

The Comprehensive Antibiotic Resistance Database (CARD; card.mcmaster.ca) combines the Antibiotic Resistance Ontology (ARO) with curated AMR gene (ARG) sequences and resistance-conferring mutations to provide an informatics framework for annotation and interpretation of resistomes. As of version 3.2.4, CARD encompasses 6627 ontology terms, 5010 reference sequences, 1933 mutations, 3004 publications, and 5057 AMR detection models that can be used by the accompanying Resistance Gene Identifier (RGI) software to annotate genomic or metagenomic sequences. Focused curation enhancements since 2020 include expanded β-lactamase curation, incorporation of likelihood-based AMR mutations for Mycobacterium tuberculosis, addition of disinfectants and antiseptics plus their associated ARGs, and systematic curation of resistance-modifying agents. This expanded curation includes 180 new AMR gene families, 15 new drug classes, 1 new resistance mechanism, and two new ontological relationships: evolutionary_variant_of and is_small_molecule_inhibitor. In silico prediction of resistomes and prevalence statistics of ARGs has been expanded to 377 pathogens, 21,079 chromosomes, 2,662 genomic islands, 41,828 plasmids and 155,606 whole-genome shotgun assemblies, resulting in collation of 322,710 unique ARG allele sequences. New features include the CARD:Live collection of community submitted isolate resistome data and the introduction of standardized 15 character CARD Short Names for ARGs to support machine learning efforts.

Outbreaks of virulent and/or drug-resistant bacteria have a significant impact on human health and major economic consequences. Genomic islands (GIs; defined as clusters of genes of probable horizontal origin) are of high interest because they disproportionately encode virulence factors, some antimicrobial-resistance (AMR) genes, and other adaptations of medical or environmental interest. While microbial genome sequencing has become rapid and inexpensive, current computational methods for GI analysis are not amenable for rapid, accurate, user-friendly and scalable comparative analysis of sets of related genomes. To help fill this gap, we have developed IslandCompare, an open-source computational pipeline for GI prediction and comparison across several to hundreds of bacterial genomes. A dynamic and interactive visualization strategy displays a bacterial core-genome phylogeny, with bacterial genomes linearly displayed at the phylogenetic tree leaves. Genomes are overlaid with GI predictions and AMR determinants from the Comprehensive Antibiotic Resistance Database (CARD), and regions of similarity between the genomes are also displayed. GI predictions are performed using Sigi-HMM and IslandPath-DIMOB, the two most precise GI prediction tools based on nucleotide composition biases, as well as a novel blast-based consistency step to improve cross-genome prediction consistency. GIs across genomes sharing sequence similarity are grouped into clusters, further aiding comparative analysis and visualization of acquisition and loss of mobile GIs in specific sub-clades. IslandCompare is an open-source software that is containerized for local use, plus available via a user-friendly, web-based interface to allow direct use by bioinformaticians, biologists and clinicians (at https://islandcompare.ca).

Since its emergence in Wuhan, China, in December 2019, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has infected ≈6 million persons worldwide. As SARS-CoV-2 spreads across the planet, we explored the range of human cells that can be infected by this virus. We isolated SARS-CoV-2 from 2 infected patients in Toronto, Canada; determined the genomic sequences; and identified single-nucleotide changes in representative populations of our virus stocks. We also tested a wide range of human immune cells for productive infection with SARS-CoV-2. We confirm that human primary peripheral blood mononuclear cells are not permissive for SARS-CoV-2. As SARS-CoV-2 continues to spread globally, it is essential to monitor single-nucleotide polymorphisms in the virus and to continue to isolate circulating viruses to determine viral genotype and phenotype by using in vitro and in vivo infection models.

Enterococcus faecium is a ubiquitous opportunistic pathogen that is exhibiting increasing levels of antimicrobial resistance (AMR). Many of the genes that confer resistance and pathogenic functions are localized on mobile genetic elements (MGEs), which facilitate their transfer between lineages. Here, features including resistance determinants, virulence factors and MGEs were profiled in a set of 1273 E. faecium genomes from two disparate geographic locations (in the UK and Canada) from a range of agricultural, clinical and associated habitats. Neither lineages of E. faecium , type A and B, nor MGEs are constrained by geographic proximity, but our results show evidence of a strong association of many profiled genes and MGEs with habitat. Many features were associated with a group of clinical and municipal wastewater genomes that are likely forming a new human-associated ecotype within type A. The evolutionary dynamics of E. faecium make it a highly versatile emerging pathogen, and its ability to acquire, transmit and lose features presents a high risk for the emergence of new pathogenic variants and novel resistance combinations. This study provides a workflow for MGE-centric surveillance of AMR in Enterococcus that can be adapted to other pathogens.

Background The Public Health Alliance for Genomic Epidemiology (PHA4GE) (https://pha4ge.org) is a global coalition that is actively working to establish consensus standards, document and share best practices, improve the availability of critical bioinformatics tools and resources, and advocate for greater openness, interoperability, accessibility, and reproducibility in public health microbial bioinformatics. In the face of the current pandemic, PHA4GE has identified a need for a fit-for-purpose, open-source SARS-CoV-2 contextual data standard. Results As such, we have developed a SARS-CoV-2 contextual data specification package based on harmonizable, publicly available community standards. The specification can be implemented via a collection template, as well as an array of protocols and tools to support both the harmonization and submission of sequence data and contextual information to public biorepositories. Conclusions Well-structured, rich contextual data add value, promote reuse, and enable aggregation and integration of disparate datasets. Adoption of the proposed standard and practices will better enable interoperability between datasets and systems, improve the consistency and utility of generated data, and ultimately facilitate novel insights and discoveries in SARS-CoV-2 and COVID-19. The package is now supported by the NCBI’s BioSample database.

Obesity results from a caloric imbalance between energy intake, absorption and expenditure. In both rodents and humans, diet-induced thermogenesis contributes to energy expenditure and involves the activation of brown adipose tissue (BAT). We hypothesize that environmental toxicants commonly used as food additives or pesticides might reduce BAT thermogenesis through suppression of uncoupling protein 1 (UCP1) and this may contribute to the development of obesity. Using a step-wise screening approach, we discover that the organophosphate insecticide chlorpyrifos suppresses UCP1 and mitochondrial respiration in BAT at concentrations as low as 1 pM. In mice housed at thermoneutrality and fed a high-fat diet, chlorpyrifos impairs BAT mitochondrial function and diet-induced thermogenesis, promoting greater obesity, non-alcoholic fatty liver disease (NAFLD) and insulin resistance. This is associated with reductions in cAMP; activation of p38MAPK and AMPK; protein kinases critical for maintaining UCP1 and mitophagy, respectively in BAT. These data indicate that the commonly used pesticide chlorpyrifos, suppresses diet-induced thermogenesis and the activation of BAT, suggesting its use may contribute to the obesity epidemic., Chlorpyrifos is a widely-used pesticide and a common residue on vegetables and fruits. Here the authors show that at non-neurotoxic doses, chlorpyrifos reduces energy expenditure, by inhibiting diet induced thermogenesis, and promotes obesity and insulin resistance.

Metformin is the most commonly prescribed medication for type 2 diabetes, owing to its glucose-lowering effects, which are mediated through the suppression of hepatic glucose production (reviewed in refs. 1-3). However, in addition to its effects on the liver, metformin reduces appetite and in preclinical models exerts beneficial effects on ageing and a number of diverse diseases (for example, cognitive disorders, cancer, cardiovascular disease) through mechanisms that are not fully understood1-3. Given the high concentration of metformin in the liver and its many beneficial effects beyond glycemic control, we reasoned that metformin may increase the secretion of a hepatocyte-derived endocrine factor that communicates with the central nervous system4. Here we show, using unbiased transcriptomics of mouse hepatocytes and analysis of proteins in human serum, that metformin induces expression and secretion of growth differentiating factor 15 (GDF15). In primary mouse hepatocytes, metformin stimulates the secretion of GDF15 by increasing the expression of activating transcription factor 4 (ATF4) and C/EBP homologous protein (CHOP; also known as DDIT3). In wild-type mice fed a high-fat diet, oral administration of metformin increases serum GDF15 and reduces food intake, body mass, fasting insulin and glucose intolerance; these effects are eliminated in GDF15 null mice. An increase in serum GDF15 is also associated with weight loss in patients with type 2 diabetes who take metformin. Although further studies will be required to determine the tissue source(s) of GDF15 produced in response to metformin in vivo, our data indicate that the therapeutic benefits of metformin on appetite, body mass and serum insulin depend on GDF15.

Metformin is the most commonly prescribed medication for type 2 diabetes, owing to its glucose-lowering effects, which are mediated through the suppression of hepatic glucose production (reviewed in refs. 1-3). However, in addition to its effects on the liver, metformin reduces appetite and in preclinical models exerts beneficial effects on ageing and a number of diverse diseases (for example, cognitive disorders, cancer, cardiovascular disease) through mechanisms that are not fully understood1-3. Given the high concentration of metformin in the liver and its many beneficial effects beyond glycemic control, we reasoned that metformin may increase the secretion of a hepatocyte-derived endocrine factor that communicates with the central nervous system4. Here we show, using unbiased transcriptomics of mouse hepatocytes and analysis of proteins in human serum, that metformin induces expression and secretion of growth differentiating factor 15 (GDF15). In primary mouse hepatocytes, metformin stimulates the secretion of GDF15 by increasing the expression of activating transcription factor 4 (ATF4) and C/EBP homologous protein (CHOP; also known as DDIT3). In wild-type mice fed a high-fat diet, oral administration of metformin increases serum GDF15 and reduces food intake, body mass, fasting insulin and glucose intolerance; these effects are eliminated in GDF15 null mice. An increase in serum GDF15 is also associated with weight loss in patients with type 2 diabetes who take metformin. Although further studies will be required to determine the tissue source(s) of GDF15 produced in response to metformin in vivo, our data indicate that the therapeutic benefits of metformin on appetite, body mass and serum insulin depend on GDF15.

Genome sequencing of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is increasingly important to monitor the transmission and adaptive evolution of the virus. The accessibility of high-throughput methods and polymerase chain reaction (PCR) has facilitated a growing ecosystem of protocols. Two differing protocols are tiling multiplex PCR and bait capture enrichment. Each method has advantages and disadvantages but a direct comparison with different viral RNA concentrations has not been performed to assess the performance of these approaches. Here we compare Liverpool amplification, ARTIC amplification, and bait capture using clinical diagnostics samples. All libraries were sequenced using an Illumina MiniSeq with data analyzed using a standardized bioinformatics workflow (SARS-CoV-2 Illumina GeNome Assembly Line; SIGNAL). One sample showed poor SARS-CoV-2 genome coverage and consensus, reflective of low viral RNA concentration. In contrast, the second sample had a higher viral RNA concentration, which yielded good genome coverage and consensus. ARTIC amplification showed the highest depth of coverage results for both samples, suggesting this protocol is effective for low concentrations. Liverpool amplification provided a more even read coverage of the SARS-CoV-2 genome, but at a lower depth of coverage. Bait capture enrichment of SARS-CoV-2 cDNA provided results on par with amplification. While only two clinical samples were examined in this comparative analysis, both the Liverpool and ARTIC amplification methods showed differing efficacy for high and low concentration samples. In addition, amplification-free bait capture enriched sequencing of cDNA is a viable method for generating a SARS-CoV-2 genome sequence and for identification of amplification artifacts.

The Public Health Alliance for Genomic Epidemiology (PHA4GE) (https://pha4ge.org) is a global coalition that is actively working to establish consensus standards, document and share best practices, improve the availability of critical bioinformatic tools and resources, and advocate for greater openness, interoperability, accessibility and reproducibility in public health microbial bioinformatics. In the face of the current pandemic, PHA4GE has identified a clear and present need for a fit-for-purpose, open source SARS-CoV-2 contextual data standard. As such, we have developed an extension to the INSDC pathogen package, providing a SARS-CoV-2 contextual data specification based on harmonisable, publicly available, community standards. The specification is implementable via a collection template, as well as an array of protocols and tools to support the harmonisation and submission of sequence data and contextual information to public repositories. Well-structured, rich contextual data adds value, promotes reuse, and enables aggregation and integration of disparate data sets. Adoption of the proposed standard and practices will better enable interoperability between datasets and systems, improve the consistency and utility of generated data, and ultimately facilitate novel insights and discoveries in SARS-CoV-2 and COVID-19.

Vaccination has transformed public health, most notably including the eradication of smallpox. Despite its profound historical importance, little is known of the origins and diversity of the viruses used in smallpox vaccination. Prior to the twentieth century, the method, source and origin of smallpox vaccinations remained unstandardised and opaque. We reconstruct and analyse viral vaccine genomes associated with smallpox vaccination from historical artefacts. Significantly, we recover viral molecules through non-destructive sampling of historical materials lacking signs of biological residues. We use the authenticated ancient genomes to reveal the evolutionary relationships of smallpox vaccination viruses within the poxviruses as a whole.

SARS-CoV-2 emerged in December 2019 in Wuhan, China and has since infected over 1.5 million people, of which over 107,000 have died. As SARS-CoV-2 spreads across the planet, speculations remain about the range of human cells that can be infected by SARS-CoV-2. In this study, we report the isolation of SARS-CoV-2 from two COVID-19 patients in Toronto, Canada. We determined the genomic sequences of the two isolates and identified single nucleotide changes in representative populations of our virus stocks. More importantly, we tested a wide range of human immune cells for productive infection with SARS-CoV-2. Here we confirm that human primary peripheral blood mononuclear cells (PBMCs) are not permissive to SARS-CoV-2. As SARS-CoV-2 continues to spread globally, it is essential to monitor small nucleotide polymorphisms in the virus and to continue to isolate circulating viruses to determine cell susceptibility and pathogenicity using in vitro and in vivo infection models.

Obesity is a major risk factor for type 2 diabetes (T2D), non-alcoholic fatty liver disease (NAFLD) and cardiovascular disease that arises from a caloric surplus of as little as 10–30 kcal per day. And while increased consumption of energy dense foods and reduced physical activity are commonly thought to be the major contributors to this caloric imbalance, diet-induced thermogenesis is a quantitatively important component of the energy balance equation. In adult humans, recent studies have indicated that diet-induced thermogenesis requires the activation of brown adipose tissue (BAT), however, the determinants regulating this process and why they may differ between individuals are not fully understood. We hypothesized that environmental toxicants commonly used as food additives or pesticides might reduce diet-induced thermogenesis through suppression of uncoupling protein 1, the defining protein of human BAT thermogenesis. Through a screening approach of pesticides/toxicants chosen from the Toxcast chem Library, we discovered that the organophosphate insecticide chlorpyrifos potently suppressed the expression of uncoupling protein 1 (UCP1) and mitochondrial respiration in brown adipocytes at concentrations as low as 1 pM. Chloropyrifos-induced suppression of brown adipocyte thermogenesis was also observed in mice fed a diet high in fat and housed at thermoneutrality where it promoted greater obesity, non-alcoholic fatty liver disease and insulin resistance. Reductions in thermogenesis by chlorpyrifos were associated with impaired activation of the β3-adrenergic receptor and protein kinases critical for regulating UCP1 and mitophagy. These data indicate that the commonly used pesticide chlorpyrifos, at doses found within the food supply, suppresses the activation of brown adipose tissue, suggesting that its use may contribute to the obesity epidemic.

BACKGROUND: Physical inactivity impairs insulin sensitivity, which is exacerbated with aging. We examined the impact of 2 wk of acute inactivity and recovery on glycemic control, and integrated rates of muscle protein synthesis in older men and women. METHODS: Twenty-two overweight, prediabetic older adults (12 men, 10 women, 69 ± 4 y) undertook 7 d of habitual activity (baseline; BL), step reduction (SR

Additional file 8: Figure S4. Regression analyses of root-to-tip genetic distance on the ML phylogeny against either year of collection or year of sequencing.

Additional file 6: Figure S2. Maximum likelihood trees for Mutter vaccine strains called in reference to alternative VACV genomes.

The Comprehensive Antibiotic Resistance Database (CARD; https://card.mcmaster.ca) is a curated resource providing reference DNA and protein sequences, detection models and bioinformatics tools on the molecular basis of bacterial antimicrobial resistance (AMR). CARD focuses on providing high-quality reference data and molecular sequences within a controlled vocabulary, the Antibiotic Resistance Ontology (ARO), designed by the CARD biocuration team to integrate with software development efforts for resistome analysis and prediction, such as CARD’s Resistance Gene Identifier (RGI) software. Since 2017, CARD has expanded through extensive curation of reference sequences, revision of the ontological structure, curation of over 500 new AMR detection models, development of a new classification paradigm and expansion of analytical tools. Most notably, a new Resistomes & Variants module provides analysis and statistical summary of in silico predicted resistance variants from 82 pathogens and over 100 000 genomes. By adding these resistance variants to CARD, we are able to summarize predicted resistance using the information included in CARD, identify trends in AMR mobility and determine previously undescribed and novel resistance variants. Here, we describe updates and recent expansions to CARD and its biocuration process, including new resources for community biocuration of AMR molecular reference data.

Identification of the nucleotide sequences encoding antibiotic resistance elements and determination of their association with antibiotic resistance are critical to improve surveillance and monitor trends in antibiotic resistance. Current methods to study antibiotic resistance in various environments rely on extensive deep sequencing or laborious culturing of fastidious organisms, both of which are heavily time-consuming operations., Identification of the nucleotide sequences encoding antibiotic resistance elements and determination of their association with antibiotic resistance are critical to improve surveillance and monitor trends in antibiotic resistance. Current methods to study antibiotic resistance in various environments rely on extensive deep sequencing or laborious culturing of fastidious organisms, both of which are heavily time-consuming operations. An accurate and sensitive method to identify both rare and common resistance elements in complex metagenomic samples is needed. Referencing the sequences in the Comprehensive Antibiotic Resistance Database, we designed a set of 37,826 probes to specifically target over 2,000 nucleotide sequences associated with antibiotic resistance in clinically relevant bacteria. Testing of this probe set on DNA libraries generated from multidrug-resistant bacteria to selectively capture resistance genes reproducibly produced higher numbers of reads on target at a greater length of coverage than shotgun sequencing. We also identified additional resistance gene sequences from human gut microbiome samples that sequencing alone was not able to detect. Our method to capture the resistome enables a sensitive means of gene detection in diverse environments where genes encoding antibiotic resistance represent less than 0.1% of the metagenome.

Metformin is the most commonly prescribed medication for type 2 diabetes, owing to its glucose-lowering effects, which are mediated through the suppression of hepatic glucose production (reviewed in refs.

Type III Secretion Systems (T3SSs) are structurally conserved nanomachines that span the inner and outer bacterial membranes, and via a protruding needle complex contact host cell membranes and deliver type III effector proteins. T3SS are phylogenetically divided into several families based on structural basal body components. Here we have studied the evolutionary and functional conservation of four T3SS proteins from the Inv/Mxi-Spa family: a cytosolic chaperone, two hydrophobic translocators that form a plasma membrane-integral pore, and the hydrophilic “tip complex” translocator that connects the T3SS needle to the translocon pore. Salmonella enterica serovar Typhimurium (S. Typhimurium), a common cause of food-borne gastroenteritis, possesses two T3SSs, one belonging to the Inv/Mxi-Spa family. We used invasion-deficient S. Typhimurium mutants as surrogates for expression of translocator orthologs identified from an extensive phylogenetic analysis, and type III effector translocation and host cell invasion as a readout for complementation efficiency, and identified several Inv/Mxi-Spa orthologs that can functionally substitute for the S. Typhimurium chaperone and translocator proteins. Functional complementation correlates with amino acid sequence identity between orthologs, but varies considerably between the four proteins. This is the first in-depth survey of the functional interchangeability of Inv/Mxi-Spa T3SS proteins acting directly at the host-pathogen interface. This article is protected by copyright. All rights reserved.

Development is a complex and well-defined process characterized by rapid cell proliferation and apoptosis. At this stage in life, a developmentally young organism is more sensitive to toxicants as compared to an adult. In response to pro-oxidant exposure, members of the Cap’n’Collar (CNC) basic leucine zipper (b-ZIP) transcription factor family (including Nfe2 and Nfe2-related factors, Nrfs) activate the expression of genes whose protein products contribute to reduced toxicity. Here, we studied the role of the CNC protein, Nfe2, in the developmental response to pro-oxidant exposure in the zebrafish (Danio rerio). Following acute waterborne exposures to diquat or tert-buytlhydroperoxide (tBOOH) at one of three developmental stages, wildtype (WT) and nfe2 knockout (KO) embryos and larvae were morphologically scored and their transcriptomes sequenced. Early in development, KO animals suffered from hypochromia that was made more severe through exposure to pro-oxidants; this phenotype in the KO may be linked to decreased expression of alas2, a gene involved in heme synthesis. WT and KO eleutheroembryos and larvae were phenotypically equally affected by exposure to pro-oxidants, where tBOOH caused more pronounced phenotypes as compared to diquat. Comparing diquat and tBOOH exposed embryos relative to the WT untreated control, a greater number of genes were up-regulated in the tBOOH condition as compared to diquat (tBOOH: 304 vs diquat: 148), including those commonly found to be differentially regulated in the vertebrate oxidative stress response (OSR) (e.g. hsp70.2, txn1, and gsr). When comparing WT and KO across all treatments and times, there were 1170 genes that were differentially expressed, of which 33 are known targets of the Nrf proteins Nrf1 and Nrf2. More specifically, in animals exposed to pro-oxidants a total of 968 genes were differentially expressed between WT and KO across developmental time, representing pathways involved in coagulation, embryonic organ development, body fluid level regulation, erythrocyte differentiation, and oxidation-reduction, amongst others. The greatest number of genes that changed in expression between WT and KO occurred in animals exposed to diquat at 2 h post fertilization (hpf). Across time and treatment, there were six genes (dhx40, cfap70, dnajb9b, slc35f4, spi-c, and gpr19) that were significantly up-regulated in KO compared to WT and four genes (fhad1, cyp4v7, nlrp12, and slc16a6a) that were significantly down-regulated. None of these genes have been previously identified as targets of Nfe2 or the Nrf family. These results demonstrate that the zebrafish Nfe2 may be a regulator of both primitive erythropoiesis and the OSR during development.

The Comprehensive Antibiotic Resistance Database (CARD; http://arpcard.mcmaster.ca) is a manually curated resource containing high quality reference data on the molecular basis of antimicrobial resistance (AMR), with an emphasis on the genes, proteins and mutations involved in AMR. CARD is ontologically structured, model centric, and spans the breadth of AMR drug classes and resistance mechanisms, including intrinsic, mutation-driven and acquired resistance. It is built upon the Antibiotic Resistance Ontology (ARO), a custom built, interconnected and hierarchical controlled vocabulary allowing advanced data sharing and organization. Its design allows the development of novel genome analysis tools, such as the Resistance Gene Identifier (RGI) for resistome prediction from raw genome sequence. Recent improvements include extensive curation of additional reference sequences and mutations, development of a unique Model Ontology and accompanying AMR detection models to power sequence analysis, new visualization tools, and expansion of the RGI for detection of emergent AMR threats. CARD curation is updated monthly based on an interplay of manual literature curation, computational text mining, and genome analysis.

Type III Secretion Systems (T3SSs) are structurally conserved nanomachines that span the inner and outer bacterial membranes, and via a protruding needle complex contact host cell membranes and deliver type III effector proteins. T3SS are phylogenetically divided into several families based on structural basal body components. Here we have studied the evolutionary and functional conservation of four T3SS proteins from the Inv/Mxi-Spa family: a cytosolic chaperone, two hydrophobic translocators that form a plasma membrane-integral pore, and the hydrophilic 'tip complex' translocator that connects the T3SS needle to the translocon pore. Salmonella enterica serovar Typhimurium (S. Typhimurium), a common cause of food-borne gastroenteritis, possesses two T3SSs, one belonging to the Inv/Mxi-Spa family. We used invasion-deficient S. Typhimurium mutants as surrogates for expression of translocator orthologs identified from an extensive phylogenetic analysis, and type III effector translocation and host cell invasion as a readout for complementation efficiency, and identified several Inv/Mxi-Spa orthologs that can functionally substitute for the S. Typhimurium chaperone and translocator proteins. Functional complementation correlates with amino acid sequence identity between orthologs, but varies considerably between the four proteins. This is the first in-depth survey of the functional interchangeability of Inv/Mxi-Spa T3SS proteins acting directly at the host-pathogen interface.

Additional file 6: Figure S2. Maximum likelihood trees for Mutter vaccine strains called in reference to alternative VACV genomes.

Additional file 13: Figure S7. Comparison of terminal damage patterns of VK01.

Additional file 8: Figure S4. Regression analyses of root-to-tip genetic distance on the ML phylogeny against either year of collection or year of sequencing.

Additional file 7: Figure S3. GrapeTree analysis of 79 OPXV.

Additional file 13: Figure S7. Comparison of terminal damage patterns of VK01.

Additional file 3: Figure S1. Terminal damage rates for VK01, VK02, VK05, VK08 and VK12 libraries mapped to VACV strain Copenhagen reference (M35027) and human mitochondrial reference (rCRS).

Additional file 11: Figure S6. Mütter catalogue # 17090.29 before and after non-destructive sampling.

Additional file 15: Figures S8–10. VK01 coverage relative to Copenhagen VACV strain M35027, GRI90 CPXV strain X93455, and Horsepox strain DQ792504.

Additional file 1. Supplementary Material and Methods.

Additional file 16: Figure S11. Maximum likelihood phylogenetic analysis of historical vaccine strains in relation to other OPXV using largest de novo assembled contig for each sample.

Additional file 11: Figure S6. Mütter catalogue # 17090.29 before and after non-destructive sampling.

Additional file 7: Figure S3. GrapeTree analysis of 79 OPXV.