Genome Sequences of 17 Diverse Pseudomonas aeruginosa Phages.

Here, we describe genome sequences of 17 Pseudomonas aeruginosa phages, including therapeutic candidates. They belong to the families Myoviridae, Podoviridae, and Siphoviridae and six different genera. The genomes ranged in size from 42,788 to 88,805 bp, with G+C contents of 52.5% to 64.3% and numbers of coding sequences from 58 to 179.

Complete Genome Sequences of 10 Phages Lytic against Multidrug-Resistant Pseudomonas aeruginosa.

We report the genome sequences of 10 Pseudomonas aeruginosa phages studied for their potential for formulation of a therapeutic cocktail; they represent the families Myoviridae, Podoviridae, and Siphoviridae Genome sizes ranged from 43,299 to 88,728 nucleotides, with G+C contents of 52.1% to 62.2%. The genomes contained 68 to 168 coding sequences.

Comparative genomic analysis of four multidrug-resistant isolates of Acinetobacter baumannii from Georgia.

OBJECTIVES: This study reports the draft genomes of four newly isolated multidrug-resistant (MDR) Acinetobacter baumannii (A. baumannii) isolates (0830, 0365, 4022, and 2846) from western Georgia to identify putative antimicrobial resistance genes (ARGs) and to determine the clonal subtypes of local clinical isolates. METHODS: An Illumina MiSeq sequencer was used to perform whole-genome sequencing (WGS). The Vitek 2 automated system was used for microbial identification and antimicrobial resistance profiling. RESULTS: Taxonomical identification as A. baumannii was confirmed by WGS. In silico analyses resolved their ARG content and clonal relatedness using the Oxford (Oxf) and Pasteur (Pas) multi-locus sequence typing schemes. Isolates 0365 and 4022 displayed similar allelic profiles corresponding to ST944Oxf/ST78Pas. Isolate 2846 displayed a different allelic profile consistent with ST19Pas/IC 1 (International or European Clone I) and exhibited a novel Oxford ST that was designated as 1868. Isolate 0830 displayed the ST78Pas allelic profile, similar to isolates 0365 and 4022, and also possessed a single allelic mismatch in the gpi gene, resulting in an ST1104Oxf allele profile in the Oxford typing scheme. CONCLUSION: Circulating MDR A. baumannii exhibited genetic heterogeneity with variations in the structure and content of genomic A. baumannii resistance islands and encoded multiple putative ARGs. This report represents the first clonal subtype information and genomic characterization of MDR A. baumannii in Georgia and may inform future epidemiological investigations.

Correlation of Host Range Expansion of Therapeutic Bacteriophage Sb-1 with Allele State at a Hypervariable Repeat Locus.

Staphylococci are frequent agents of health care-associated infections and include methicillin-resistant Staphylococcus aureus (MRSA), which is resistant to first-line antibiotic treatments. Bacteriophage (phage) therapy is a promising alternative antibacterial option to treat MRSA infections. S. aureus-specific phage Sb-1 has been widely used in Georgia to treat a variety of human S. aureus infections. Sb-1 has a broad host range within S. aureus, including MRSA strains, and its host range can be further expanded by adaptation to previously resistant clinical isolates. The susceptibilities of a panel of 25 genetically diverse clinical MRSA isolates to Sb-1 phage were tested, and the phage had lytic activity against 23 strains (92%). The adapted phage stock (designated Sb-1A) was tested in comparison with the parental phage (designated Sb-1P). Sb-1P had lytic activity against 78/90 strains (87%) in an expanded panel of diverse global S. aureus isolates, while eight additional strains in this panel were susceptible to Sb-1A (lytic against 86/90 strains [96%]). The Sb-1A stock was shown to be a mixed population of phage clones, including approximately 4% expanded host range mutants, designated Sb-1M. In an effort to better understand the genetic basis for this host range expansion, we sequenced the complete genomes of the parental Sb-1P and two Sb-1M mutants. Comparative genomic analysis revealed a hypervariable complex repeat structure in the Sb-1 genome that had a distinct allele that correlated with the host range expansion. This hypervariable region was previously uncharacterized in Twort-like phages and represents a novel putative host range determinant.IMPORTANCE Because of limited therapeutic options, infections caused by methicillin-resistant Staphylococcus aureus represent a serious problem in both civilian and military health care settings. Phages have potential as alternative antibacterial agents that can be used in combination with antibiotic drugs. For decades, phage Sb-1 has been used in former Soviet Union countries for antistaphylococcal treatment in humans. The therapeutic spectrum of activity of Sb-1 can be increased by selecting mutants of the phage with expanded host ranges. In this work, the host range of phage Sb-1 was expanded in the laboratory, and a hypervariable region in its genome was identified with a distinct allele state that correlated with this host range expansion. These results provide a genetic basis for better understanding the mechanisms of phage host range expansion.

Intra-epidemic genome variation in highly pathogenic African swine fever virus (ASFV) from the country of Georgia.

BACKGROUND: African swine fever virus (ASFV) causes an acute hemorrhagic infection in suids with a mortality rate of up to 100%. No vaccine is available and the potential for catastrophic disease in Europe remains elevated due to the ongoing ASF epidemic in Russia and Baltic countries. To date, intra-epidemic whole-genome variation for ASFV has not been reported. To provide a more comprehensive baseline for genetic variation early in the ASF outbreak, we sequenced two Georgian ASFV samples, G-2008/1 and G-2008/2, derived from domestic porcine blood collected in 2008. METHODS: Genomic DNA was extracted directly from low-volume ASFV PCR-positive porcine blood samples and subjected to next generation sequencing on the Illumina Miseq platform. De novo and mapped sequence assemblies were performed using CLCBio software. Genomic illustrations, sequence alignments and assembly figures were generated using Geneious v10.2.4. Sequence repeat architecture was analyzed using DNASTAR GeneQuest 14.1.0. RESULTS: The G-2008/1 and G-2008/2 genomes were distinguished from each other by coding changes in seven genes, including MGF 110-1 L, X69R, MGF 505-10R, EP364R, H233R, E199L, and MGF 360-21R in addition to eight homopolymer tract variations. The 2008/2 genome possessed a novel allele state at a previously undescribed intergenic repeat locus between genes C315R and C147L. The C315R/C147L locus represents the earliest observed variable repeat sequence polymorphism reported among isolates from this epidemic. No sequence variation was observed in conventional ASFV subtyping markers. The two genomes exhibited complete collinearity and identical gene content with the Georgia 2007/1 reference genome. Approximately 56 unique homopolymer A/T-tract variations were identified that were unique to the Georgia 2007/1 genome. In both 2008 genomes, within-sample sequence read heterogeneity was evident at six homopolymeric G/C-tracts confined to the known hypervariable ~ 7 kb region in the left terminal region of the genome. CONCLUSIONS: This is the first intra-epidemic comparative genomic analysis reported for ASFV and provides insight into the intra-epidemic microevolution of ASFV. The genomes reported here, in addition to the G-2007/1 genome, provide an early baseline for future genome-level comparisons and epidemiological tracing efforts.

Genome Sequence of the Soviet/Russian Bacillus anthracis Vaccine Strain 55-VNIIVViM.

Bacillus anthracis strain 55-VNIIVViM is a live-attenuated nonencapsulated Soviet/Russian veterinary anthrax vaccine strain. We report here the genome of 55-VNIIVViM and confirm its phylogenetic placement in the global population structure of B. anthracis.

Seroepidemiology and molecular diversity of Leishmania donovani complex in Georgia.

BACKGROUND: Leishmaniasis includes multiple clinical syndromes, most notably visceral, cutaneous, and mucosal forms. Visceral leishmaniasis (VL), also known as kala-azar, is a potentially fatal disease endemic to large parts of Africa and Asia, and in South-Eastern Europe (Greece, Turkey, Georgia). Visceral leishmaniasis is a parasitic zoonosis caused by species of the L. donovani complex. In the classical epidemiological model the main reservoir for VL are canines. METHODS: The study included a cohort of 513 individuals of both genders (190 males and 323 females) from the ages of 1 to 70 years that were screened in ten villages across two districts in Kakheti using the Kalazar Detect™ rK39 rapid diagnostic test. The phylogenetic diversity patterns of local strains, based on the rDNA internal transcribed spacer (ITS) sequences, were assessed for samples obtained from patients with suspected L. donovani infection, from canine reservoirs and from Phlebotomus sand flies obtained from different geographical areas of Georgia and from Azerbaijan. RESULTS: Out of a total of 600 domestic dog blood samples 95 (15.8 %) were positive by rK39 rapid diagnostic tests. For symptomatic domestic dogs, the testing of conjunctival swabs or bone marrow aspirates revealed a higher VL incidence in Kvareli District (Kvareli; 19.4 %, n = 329) compared with that observed for Sagarejo District (Sagarejo; 11.4 %, n = 271). A total of 231 sand flies of both genders were collected during the 2-month period; of the 114 females, 1.75 % were PCR positive for the presence of Leishmania spp. CONCLUSIONS: VL infection rates remain high in both canines and humans in Georgia, with disease in several known natural foci. The genetic relationships derived from rDNA internal transcribed spacer (ITS) sequence comparisons identified genetic subgroups, revealing preliminary insights into the genetic structure of L. donovani complex members currently circulating in the South Caucasus and demonstrates the utility of ITS-based genotyping in the resource-limited country of Georgia.

Whole genome sequence comparison of ten diagnostic brucellaphages propagated on two Brucella abortus hosts.

BACKGROUND: Recently the genome sequences of two brucellaphages, isolated in Georgia (Tb) and Mexico (Pr) were reported revealing pronounced sequence homogeneity and the presence of two major indels discriminating the two phages. Subsequent genome sequencing of six diagnostic brucellaphages: Tbilisi (Tb), Firenze (Fz), Weybridge (Wb), S708, Berkeley (Bk) and R/C phages identified three major genetic groups. However, the propensity for fine-scale genetic variability of diverse brucellaphages grown on multiple hosts within a single Brucella species remains unknown. METHODS: We sequenced the complete genomes of ten brucellaphages following initial propagation on B. abortus strain 141 and after subsequent propagation on B. abortus strain S19. All brucellaphages were isolated and propagated at the Eliava Institute including the original Tb phage. Genomic libraries were quantified using the Qbit and sheared on the Covaris M220. QC for fragmentation was performed on the BioAnalyzer 2100. DNA libraries were prepared using an Illumina Paired-End protocol and sequenced on the Illumina MiSeq. Sequence analysis was performed using Geneious and MEGA software. RESULTS: Comparative whole genome sequence analysis revealed genetic homogeneity consistent with previously published data as well as multiple nucleotide variations. Genomic changes as a result of passages were observed in similar genes and predominantly occurred at identical sites in separate phages. Multiple instances of within-sample genetic heterogeneity were observed often at shared genomics positions across phages. Positive selection was detected in the tail collar protein gene. We also identified a Staphylothermus marinus F1-like CRISPR spacer and sequences orthologous to both prophage antirepressors of Brucella spp. and intergenic sequences encoded by Ochrobactrum anthropi. CONCLUSION: We surveyed whole genome level diversity in phage lytic for B. abortus as they are propagated on alternate vaccine strains within the species. Our data extend previous results indicating select variable hotspots and broad genomic homogeneity as well as multiple common polymorphisms and within-sample variation. These data also provide additional genomes for future reference in comparative studies involving the molecular evolution and host specificity of brucellaphages.

Comparative whole genome analysis of six diagnostic brucellaphages.

Whole genome sequencing of six diagnostic brucellaphages, Tbilisi (Tb), Firenze (Fz), Weybridge (Wb), S708, Berkeley (Bk) and R/C, was followed with genomic comparisons including recently described genomes of the Tb phage from Mexico (TbM) and Pr phage to elucidate genomic diversity and candidate host range determinants. Comparative whole genome analysis revealed high sequence homogeneity among these brucellaphage genomes and resolved three genetic groups consistent with defined host range phenotypes. Group I was composed of Tb and Fz phages that are predominantly lytic for Brucella abortus and Brucella neotomae; Group II included Bk, R/C, and Pr phages that are lytic mainly for B. abortus, Brucella melitensis and Brucella suis; Group III was composed of Wb and S708 phages that are lytic for B. suis, B. abortus and B. neotomae. We found that the putative phage collar protein is a variable locus with features that may be contributing to the host specificities exhibited by different brucellaphage groups. The presence of several candidate host range determinants is illustrated herein for future dissection of the differential host specificity observed among these phages.

Phylogeography of Francisella tularensis subspecies holarctica from the country of Georgia.

BACKGROUND: Francisella tularensis, the causative agent of tularemia, displays subspecies-specific differences in virulence, geographic distribution, and genetic diversity. F. tularensis subsp. holarctica is widely distributed throughout the Northern Hemisphere. In Europe, F. tularensis subsp. holarctica isolates have largely been assigned to two phylogenetic groups that have specific geographic distributions. Most isolates from Western Europe are assigned to the B.Br.FTNF002-00 group, whereas most isolates from Eastern Europe are assigned to numerous lineages within the B.Br.013 group. The eastern geographic extent of the B.Br.013 group is currently unknown due to a lack of phylogenetic knowledge about populations at the European/Asian juncture and in Asia. In this study, we address this knowledge gap by describing the phylogenetic structure of F. tularensis subsp. holarctica isolates from the country of Georgia, and by placing these isolates into a global phylogeographic context. RESULTS: We identified a new genetic lineage of F. tularensis subsp. holarctica from Georgia that belongs to the B.Br.013 group. This new lineage is genetically and geographically distinct from lineages previously described from the B.Br.013 group from Central-Eastern Europe. Importantly, this new lineage is basal within the B.Br.013 group, indicating the Georgian lineage diverged before the diversification of the other known B.Br.013 lineages. Although two isolates from the Georgian lineage were collected nearby in the Ukrainian region of Crimea, all other global isolates assigned to this lineage were collected in Georgia. This restricted geographic distribution, as well as the high levels of genetic diversity within the lineage, is consistent with a relatively older origin and localized differentiation. CONCLUSIONS: We identified a new lineage of F. tularensis subsp. holarctica from Georgia that appears to have an older origin than any other diversified lineages previously described from the B.Br.013 group. This finding suggests that additional phylogenetic studies of F. tularensis subsp. holarctica populations in Eastern Europe and Asia have the potential to yield important new insights into the evolutionary history and phylogeography of this broadly dispersed F. tularensis subspecies.

Comparative whole genome sequence analysis of wild-type and cidofovir-resistant monkeypoxvirus.

We performed whole genome sequencing of a cidofovir {[(S)-1-(3-hydroxy-2-phosphonylmethoxy-propyl) cytosine] [HPMPC]}-resistant (CDV-R) strain of Monkeypoxvirus (MPV). Whole-genome comparison with the wild-type (WT) strain revealed 55 single-nucleotide polymorphisms (SNPs) and one tandem-repeat contraction. Over one-third of all identified SNPs were located within genes comprising the poxvirus replication complex, including the DNA polymerase, RNA polymerase, mRNA capping methyltransferase, DNA processivity factor, and poly-A polymerase. Four polymorphic sites were found within the DNA polymerase gene. DNA polymerase mutations observed at positions 314 and 684 in MPV were consistent with CDV-R loci previously identified in Vaccinia virus (VACV). These data suggest the mechanism of CDV resistance may be highly conserved across Orthopoxvirus (OPV) species. SNPs were also identified within virulence genes such as the A-type inclusion protein, serine protease inhibitor-like protein SPI-3, Schlafen ATPase and thymidylate kinase, among others. Aberrant chain extension induced by CDV may lead to diverse alterations in gene expression and viral replication that may result in both adaptive and attenuating mutations. Defining the potential contribution of substitutions in the replication complex and RNA processing machinery reported here may yield further insight into CDV resistance and may augment current therapeutic development strategies.

Global genetic population structure of Bacillus anthracis.

Anthrax, caused by the bacterium Bacillus anthracis, is a disease of historical and current importance that is found throughout the world. The basis of its historical transmission is anecdotal and its true global population structure has remained largely cryptic. Seven diverse B. anthracis strains were whole-genome sequenced to identify rare single nucleotide polymorphisms (SNPs), followed by phylogenetic reconstruction of these characters onto an evolutionary model. This analysis identified SNPs that define the major clonal lineages within the species. These SNPs, in concert with 15 variable number tandem repeat (VNTR) markers, were used to subtype a collection of 1,033 B. anthracis isolates from 42 countries to create an extensive genotype data set. These analyses subdivided the isolates into three previously recognized major lineages (A, B, and C), with further subdivision into 12 clonal sub-lineages or sub-groups and, finally, 221 unique MLVA15 genotypes. This rare genomic variation was used to document the evolutionary progression of B. anthracis and to establish global patterns of diversity. Isolates in the A lineage are widely dispersed globally, whereas the B and C lineages occur on more restricted spatial scales. Molecular clock models based upon genome-wide synonymous substitutions indicate there was a massive radiation of the A lineage that occurred in the mid-Holocene (3,064-6,127 ybp). On more recent temporal scales, the global population structure of B. anthracis reflects colonial-era importation of specific genotypes from the Old World into the New World, as well as the repeated industrial importation of diverse genotypes into developed countries via spore-contaminated animal products. These findings indicate humans have played an important role in the evolution of anthrax by increasing the proliferation and dispersal of this now global disease. Finally, the value of global genotypic analysis for investigating bioterrorist-mediated outbreaks of anthrax is demonstrated.

Francisella tularensis in the United States.

The causative agent of tularemia, Francisella tularensis, is a formidable biologic agent that occurs naturally throughout North America. We examined genetic and spatial diversity patterns among 161 US F. tularensis isolates by using a 24-marker multiple-locus variable-number tandem repeat analysis (MLVA) system. MLVA identified 126 unique genotypes. Phylogenetic analyses showed patterns similar to recently reported global-scale analyses. We observed clustering by subspecies, low genetic diversity within F. tularensis subsp. holarctica, and division of F. tularensis subsp. tularensis into 2 distinct subpopulations: A.I. and A.II. The 2 F. tularensis subsp. tularensis subpopulations also represent geographically distinct groups; A.I. occurs primarily in the central United States, and A.II. occurs primarily in the western United States. These spatial distributions are correlated with geographic ranges of particular vectors, hosts of tularemia, and abiotic factors. These correlates provide testable hypotheses regarding ecologic factors associated with maintaining tularemia foci.

Worldwide genetic relationships among Francisella tularensis isolates determined by multiple-locus variable-number tandem repeat analysis.

The intracellular bacterium Francisella tularensis is the causative agent of tularemia and poses a serious threat as an agent of bioterrorism. We have developed a highly effective molecular subtyping system from 25 variable-number tandem repeat (VNTR) loci. In our study, multiple-locus VNTR analysis (MLVA) was used to analyze genetic relationships and potential population structure within a global collection of 192 F. tularensis isolates, including representatives from each of the four subspecies. The VNTR loci displayed between 2 and 31 alleles with Nei's diversity values between 0.05 and 0.95. Neighbor-joining cluster analysis of VNTR data revealed 120 genotypes among the 192 F. tularensis isolates, including accurate subspecies identification. F. tularensis subsp. tularensis (type A) isolates showed great diversity at VNTR loci, while F. tularensis subsp. holarctica (type B) isolates showed much lower levels despite a much broader geographical prevalence. The resolution of two distinct clades within F. tularensis subsp. tularensis (designated A.I and A.II) revealed a previously unrecognized genetic division within this highly virulent subspecies. F. tularensis subsp. holarctica appears to have recently spread globally across continents from a single origin, while F. tularensis subsp. tularensis has a long and complex evolutionary history almost exclusively in North America. The sole non-North American type A isolates (Slovakian) were closely related to the SCHU S4 strain. Significant linkage disequilibrium was detected among VNTR loci of F. tularensis consistent with a clonal population structure. Overall, this work greatly augments the study of tularemia ecology and epidemiology, while providing a framework for future forensic analysis of F. tularensis isolates.

Strain typing of Borrelia burgdorferi, Borrelia afzelii, and Borrelia garinii by using multiple-locus variable-number tandem repeat analysis.

Human Lyme borreliosis (LB) is the most prevalent arthropod-borne infection in temperate climate zones around the world and is caused by Borrelia spirochetes. We have identified 10 variable-number tandem repeat (VNTR) loci present within the genome of Borrelia burgdorferi and subsequently developed a multiple-locus VNTR analysis (MLVA) typing system for this disease agent. We report here the successful application of MLVA for strain discrimination among a group of 41 globally diverse Borrelia isolates including B. burgdorferi, B. afzelii, and B. garinii. PCR assays displayed diversity at these loci, with total allele numbers ranging from two to nine and Nei's diversity (D) values ranging from 0.10 to 0.87. The average D value was 0.53 across all VNTR loci. A clear correlation exists between the repeat copy number and the D value (r = 0.62) or the number of alleles (r = 0.93) observed across diverse strains. Cluster analysis by the unweighted pair-group method with arithmetic means resolved the 30 observed unique Borrelia genotypes into five distinct groups. B. burgdorferi, B. afzelii, and B. garinii clustered into distinct affiliations, consistent with current 16S rRNA phylogeny studies. Genetic similarity and diversity suggest that B. afzelii and B. garinii are close relatives and were perhaps recently derived from B. burgdorferi. MLVA provides both phylogenetic relationships and additional resolution to discriminate among strains of Borrelia species. This new level of strain identification and discrimination will allow more detailed epidemiological and phylogenetic analysis in future studies.