The genomes of three North American orthopoxviruses.

The complete genomes of a skunkpox, volepox, and raccoonpox virus were sequenced and annotated. Phylogenetic analysis of these genomes indicates that although these viruses are all orthopoxviruses, they form a distinct clade to the other known species. This supports the ancient divergence of the North American orthopoxviruses from other members of the orthopoxviruses. Only two open reading frames appear to be unique to this group of viruses, but a relatively small number of insertions/deletions contribute to the varied gene content of this clade. The availability of these genomes will help determine whether skunkpox and volepox viruses share the characteristics that make raccoonpox a useful vaccine vector.

Isolation and enrichment of Cryptosporidium DNA and verification of DNA purity for whole-genome sequencing.

Whole-genome sequencing of Cryptosporidium spp. is hampered by difficulties in obtaining sufficient, highly pure genomic DNA from clinical specimens. In this study, we developed procedures for the isolation and enrichment of Cryptosporidium genomic DNA from fecal specimens and verification of DNA purity for whole-genome sequencing. The isolation and enrichment of genomic DNA were achieved by a combination of three oocyst purification steps and whole-genome amplification (WGA) of DNA from purified oocysts. Quantitative PCR (qPCR) analysis of WGA products was used as an initial quality assessment of amplified genomic DNA. The purity of WGA products was assessed by Sanger sequencing of cloned products. Next-generation sequencing tools were used in final evaluations of genome coverage and of the extent of contamination. Altogether, 24 fecal specimens of Cryptosporidium parvum, C. hominis, C. andersoni, C. ubiquitum, C. tyzzeri, and Cryptosporidium chipmunk genotype I were processed with the procedures. As expected, WGA products with low (<16.0) threshold cycle (CT) values yielded mostly Cryptosporidium sequences in Sanger sequencing. The cloning-sequencing analysis, however, showed significant contamination in 5 WGA products (proportion of positive colonies derived from Cryptosporidium genomic DNA, ≤25%). Following this strategy, 20 WGA products from six Cryptosporidium species or genotypes with low (mostly <14.0) CT values were submitted to whole-genome sequencing, generating sequence data covering 94.5% to 99.7% of Cryptosporidium genomes, with mostly minor contamination from bacterial, fungal, and host DNA. These results suggest that the described strategy can be used effectively for the isolation and enrichment of Cryptosporidium DNA from fecal specimens for whole-genome sequencing.

A distinct lineage of influenza A virus from bats.

Influenza A virus reservoirs in animals have provided novel genetic elements leading to the emergence of global pandemics in humans. Most influenza A viruses circulate in waterfowl, but those that infect mammalian hosts are thought to pose the greatest risk for zoonotic spread to humans and the generation of pandemic or panzootic viruses. We have identified an influenza A virus from little yellow-shouldered bats captured at two locations in Guatemala. It is significantly divergent from known influenza A viruses. The HA of the bat virus was estimated to have diverged at roughly the same time as the known subtypes of HA and was designated as H17. The neuraminidase (NA) gene is highly divergent from all known influenza NAs, and the internal genes from the bat virus diverged from those of known influenza A viruses before the estimated divergence of the known influenza A internal gene lineages. Attempts to propagate this virus in cell cultures and chicken embryos were unsuccessful, suggesting distinct requirements compared with known influenza viruses. Despite its divergence from known influenza A viruses, the bat virus is compatible for genetic exchange with human influenza viruses in human cells, suggesting the potential capability for reassortment and contributions to new pandemic or panzootic influenza A viruses.

Phylogenetic and ecologic perspectives of a monkeypox outbreak, southern Sudan, 2005.

Identification of human monkeypox cases during 2005 in southern Sudan (now South Sudan) raised several questions about the natural history of monkeypox virus (MPXV) in Africa. The outbreak area, characterized by seasonally dry riverine grasslands, is not identified as environmentally suitable for MPXV transmission. We examined possible origins of this outbreak by performing phylogenetic analysis of genome sequences of MPXV isolates from the outbreak in Sudan and from differing localities. We also compared the environmental suitability of study localities for monkeypox transmission. Phylogenetically, the viruses isolated from Sudan outbreak specimens belong to a clade identified in the Congo Basin. This finding, added to the political instability of the area during the time of the outbreak, supports the hypothesis of importation by infected animals or humans entering Sudan from the Congo Basin, and person-to-person transmission of virus, rather than transmission of indigenous virus from infected animals to humans.

Comparative genomics of Vibrio cholerae from Haiti, Asia, and Africa.

Cholera was absent from the island of Hispaniola at least a century before an outbreak that began in Haiti in the fall of 2010. Pulsed-field gel electrophoresis (PFGE) analysis of clinical isolates from the Haiti outbreak and recent global travelers returning to the United States showed indistinguishable PFGE fingerprints. To better explore the genetic ancestry of the Haiti outbreak strain, we acquired 23 whole-genome Vibrio cholerae sequences: 9 isolates obtained in Haiti or the Dominican Republic, 12 PFGE pattern-matched isolates linked to Asia or Africa, and 2 nonmatched outliers from the Western Hemisphere. Phylogenies for whole-genome sequences and core genome single-nucleotide polymorphisms showed that the Haiti outbreak strain is genetically related to strains originating in India and Cameroon. However, because no identical genetic match was found among sequenced contemporary isolates, a definitive genetic origin for the outbreak in Haiti remains speculative.

A computational genomics pipeline for prokaryotic sequencing projects.

MOTIVATION: New sequencing technologies have accelerated research on prokaryotic genomes and have made genome sequencing operations outside major genome sequencing centers routine. However, no off-the-shelf solution exists for the combined assembly, gene prediction, genome annotation and data presentation necessary to interpret sequencing data. The resulting requirement to invest significant resources into custom informatics support for genome sequencing projects remains a major impediment to the accessibility of high-throughput sequence data. RESULTS: We present a self-contained, automated high-throughput open source genome sequencing and computational genomics pipeline suitable for prokaryotic sequencing projects. The pipeline has been used at the Georgia Institute of Technology and the Centers for Disease Control and Prevention for the analysis of Neisseria meningitidis and Bordetella bronchiseptica genomes. The pipeline is capable of enhanced or manually assisted reference-based assembly using multiple assemblers and modes; gene predictor combining; and functional annotation of genes and gene products. Because every component of the pipeline is executed on a local machine with no need to access resources over the Internet, the pipeline is suitable for projects of a sensitive nature. Annotation of virulence-related features makes the pipeline particularly useful for projects working with pathogenic prokaryotes. AVAILABILITY AND IMPLEMENTATION: The pipeline is licensed under the open-source GNU General Public License and available at the Georgia Tech Neisseria Base (http://nbase.biology.gatech.edu/). The pipeline is implemented with a combination of Perl, Bourne Shell and MySQL and is compatible with Linux and other Unix systems.

Smallpox virus resequencing GeneChips can also rapidly ascertain species status for some zoonotic non-variola orthopoxviruses.

We recently developed a set of seven resequencing GeneChips for the rapid sequencing of Variola virus strains in the WHO Repository of the Centers for Disease Control and Prevention. In this study, we attempted to hybridize these GeneChips with some known non-Variola orthopoxvirus isolates, including monkeypox, cowpox, and vaccinia viruses, for rapid detection.

Distinct evolutionary patterns of Neisseria meningitidis serogroup B disease outbreaks at two universities in the USA.

Neisseria meningitidis serogroup B (MnB) was responsible for two independent meningococcal disease outbreaks at universities in the USA during 2013. The first at University A in New Jersey included nine confirmed cases reported between March 2013 and March 2014. The second outbreak occurred at University B in California, with four confirmed cases during November 2013. The public health response to these outbreaks included the approval and deployment of a serogroup B meningococcal vaccine that was not yet licensed in the USA. This study investigated the use of whole-genome sequencing(WGS) to examine the genetic profile of the disease-causing outbreak isolates at each university. Comparative WGS revealed differences in evolutionary patterns between the two disease outbreaks. The University A outbreak isolates were very closely related, with differences primarily attributed to single nucleotide polymorphisms/insertion-deletion (SNP/indel) events. In contrast, the University B outbreak isolates segregated into two phylogenetic clades, differing in large part due to recombination events covering extensive regions (>30 kb) of the genome including virulence factors. This high-resolution comparison of two meningococcal disease outbreaks further demonstrates the genetic complexity of meningococcal bacteria as related to evolution and disease virulence.

Detection of minority drug resistant mutations in Malawian HIV-1 subtype C-positive patients initiating and on first-line antiretroviral therapy.

BACKGROUND: Minority drug resistance mutations (DRMs) that are often missed by Sanger sequencing are clinically significant, as they can cause virologic failure in individuals treated with antiretroviral therapy (ART) drugs. OBJECTIVE: This study aimed to estimate the prevalence of minor DRMs among patients enrolled in a Malawi HIV drug resistance monitoring survey at baseline and at one year after initiation of ART. METHODS: Forty-one plasma specimens collected from HIV-1 subtype C-positive patients and seven clonal control samples were analysed using ultra-deep sequencing technology. RESULTS: Deep sequencing identified all 72 DRMs detected by Sanger sequencing at the level of ≥20% and 79 additional minority DRMs at the level of < 20% from the 41 Malawian clinical specimens. Overall, DRMs were detected in 85% of pre-ART and 90.5% of virologic failure patients by deep sequencing. Among pre-ART patients, deep sequencing identified a statistically significant higher prevalence of DRMs to nucleoside reverse transcriptase inhibitors (NRTIs) compared with Sanger sequencing. The difference was mainly due to the high prevalence of minority K65R and M184I mutations. Most virologic failure patients harboured DRMs against both NRTIs and non-nucleoside reverse transcriptase inhibitors (NNRTIs). These minority DRMs contributed to the increased or enhanced virologic failures in these patients. CONCLUSION: The results revealed the presence of minority DRMs to NRTIs and NNRTIs in specimens collected at baseline and virologic failure time points. These minority DRMs not only increased resistance levels to NRTIs and NNRTIs for the prescribed ART, but also expanded resistance to additional major first-line ART drugs. This study suggested that drug resistance testing that uses more sensitive technologies, is needed in this setting.

Genome Sequence of Vaccinia virus Strain Lister-Butantan, a Lister Vaccine Variant Used during a Smallpox Eradication Campaign in Brazil.

Here, we report the 187.8-kb genome sequence of Vaccinia virus Lister-Butantan, which was used in Brazil during the WHO smallpox eradication campaign. Its genome showed an average similarity of 98.18% with the original Lister isolate, highlighting the low divergence among related Vaccinia virus vaccine strains, even after several passages in animals and cell culture.

Genome Structural Diversity among 31 Bordetella pertussis Isolates from Two Recent U.S. Whooping Cough Statewide Epidemics.

During 2010 and 2012, California and Vermont, respectively, experienced statewide epidemics of pertussis with differences seen in the demographic affected, case clinical presentation, and molecular epidemiology of the circulating strains. To overcome limitations of the current molecular typing methods for pertussis, we utilized whole-genome sequencing to gain a broader understanding of how current circulating strains are causing large epidemics. Through the use of combined next-generation sequencing technologies, this study compared de novo, single-contig genome assemblies from 31 out of 33 Bordetella pertussis isolates collected during two separate pertussis statewide epidemics and 2 resequenced vaccine strains. Final genome architecture assemblies were verified with whole-genome optical mapping. Sixteen distinct genome rearrangement profiles were observed in epidemic isolate genomes, all of which were distinct from the genome structures of the two resequenced vaccine strains. These rearrangements appear to be mediated by repetitive sequence elements, such as high-copy-number mobile genetic elements and rRNA operons. Additionally, novel and previously identified single nucleotide polymorphisms were detected in 10 virulence-related genes in the epidemic isolates. Whole-genome variation analysis identified state-specific variants, and coding regions bearing nonsynonymous mutations were classified into functional annotated orthologous groups. Comprehensive studies on whole genomes are needed to understand the resurgence of pertussis and develop novel tools to better characterize the molecular epidemiology of evolving B. pertussis populations. IMPORTANCE Pertussis, or whooping cough, is the most poorly controlled vaccine-preventable bacterial disease in the United States, which has experienced a resurgence for more than a decade. Once viewed as a monomorphic pathogen, B. pertussis strains circulating during epidemics exhibit diversity visible on a genome structural level, previously undetectable by traditional sequence analysis using short-read technologies. For the first time, we combine short- and long-read sequencing platforms with restriction optical mapping for single-contig, de novo assembly of 31 isolates to investigate two geographically and temporally independent U.S. pertussis epidemics. These complete genomes reshape our understanding of B. pertussis evolution and strengthen molecular epidemiology toward one day understanding the resurgence of pertussis.

Three Genome Sequences of Legionella pneumophila subsp. pascullei Associated with Colonization of a Health Care Facility.

Here, we report the complete genome sequences of three Legionella pneumophila subsp. pascullei strains (including both serogroup 1 and 5 strains) that were found in the same health care facility in 1982 and 2012.

Serro 2 Virus Highlights the Fundamental Genomic and Biological Features of a Natural Vaccinia Virus Infecting Humans.

Vaccinia virus (VACV) has been implicated in infections of dairy cattle and humans, and outbreaks have substantially impacted local economies and public health in Brazil. During a 2005 outbreak, a VACV strain designated Serro 2 virus (S2V) was collected from a 30-year old male milker. Our aim was to phenotypically and genetically characterize this VACV Brazilian isolate. S2V produced small round plaques without associated comets when grown in BSC40 cells. Furthermore, S2V was less virulent than the prototype strain VACV-Western Reserve (WR) in a murine model of intradermal infection, producing a tiny lesion with virtually no surrounding inflammation. The genome of S2V was sequenced by primer walking. The coding region spans 184,572 bp and contains 211 predicted genes. Mutations in envelope genes specifically associated with small plaque phenotypes were not found in S2V; however, other alterations in amino acid sequences within these genes were identified. In addition, some immunomodulatory genes were truncated in S2V. Phylogenetic analysis using immune regulatory-related genes, besides the hemagglutinin gene, segregated the Brazilian viruses into two clusters, grouping the S2V into Brazilian VACV group 1. S2V is the first naturally-circulating human-associated VACV, with a low passage history, to be extensively genetically and phenotypically characterized.

Draft Genome Sequence of Legionella pneumophila D-5864, a Serogroup 6 Strain.

Legionella pneumophila is the leading etiology of legionellosis infections in North America and Europe. Here we report the draft genome sequence of L. pneumophila D-5864, a serogroup 6 strain, which was isolated from a bronchial alveolar lavage specimen of a male patient from Arizona in 2009. Genes within the lipopolysaccharide (LPS)-biosynthesis region could potentially be determinants of serogroup specificity.

Identification and Characterization of Microsatellite Markers Derived from the Whole Genome Analysis of Taenia solium.

BACKGROUND: Infections with Taenia solium are the most common cause of adult acquired seizures worldwide, and are the leading cause of epilepsy in developing countries. A better understanding of the genetic diversity of T. solium will improve parasite diagnostics and transmission pathways in endemic areas thereby facilitating the design of future control measures and interventions. Microsatellite markers are useful genome features, which enable strain typing and identification in complex pathogen genomes. Here we describe microsatellite identification and characterization in T. solium, providing information that will assist in global efforts to control this important pathogen. METHODS: For genome sequencing, T. solium cysts and proglottids were collected from Huancayo and Puno in Peru, respectively. Using next generation sequencing (NGS) and de novo assembly, we assembled two draft genomes and one hybrid genome. Microsatellite sequences were identified and 36 of them were selected for further analysis. Twenty T. solium isolates were collected from Tumbes in the northern region, and twenty from Puno in the southern region of Peru. The size-polymorphism of the selected microsatellites was determined with multi-capillary electrophoresis. We analyzed the association between microsatellite polymorphism and the geographic origin of the samples. RESULTS: The predicted size of the hybrid (proglottid genome combined with cyst genome) T. solium genome was 111 MB with a GC content of 42.54%. A total of 7,979 contigs (>1,000 nt) were obtained. We identified 9,129 microsatellites in the Puno-proglottid genome and 9,936 in the Huancayo-cyst genome, with 5 or more repeats, ranging from mono- to hexa-nucleotide. Seven microsatellites were polymorphic and 29 were monomorphic within the analyzed isolates. T. solium tapeworms were classified into two genetic groups that correlated with the North/South geographic origin of the parasites. CONCLUSIONS/SIGNIFICANCE: The availability of draft genomes for T. solium represents a significant step towards the understanding the biology of the parasite. We report here a set of T. solium polymorphic microsatellite markers that appear promising for genetic epidemiology studies.

Finished Annotated Genome Sequence of Burkholderia pseudomallei Strain Bp1651, a Multidrug-Resistant Clinical Isolate.

Burkholderia pseudomallei strain Bp1651, a human isolate, is resistant to all clinically relevant antibiotics. We report here on the finished genome sequence assembly and annotation of the two chromosomes of this strain. This genome sequence may assist in understanding the mechanisms of antimicrobial resistance for this pathogenic species.

Draft Genome Sequence of Bacillus cereus Strain BcFL2013, a Clinical Isolate Similar to G9241.

Bacillus cereus strains, such as G9241, causing anthrax-like illnesses have recently been discovered. We report the genome sequence of a clinical strain, B. cereus BcFL2013, which is similar to G9241, recovered from a patient in Florida.

Draft Whole-Genome Sequences of Nine Non-O157 Shiga Toxin-Producing Escherichia coli Strains.

Shiga toxin-producing Escherichia coli (STEC) is an important food-borne pathogen. Here, we report the draft whole-genome sequences of nine STEC strains isolated from clinical cases in the United States. This is the first report of such information for STEC of serotypes O69, H11, O145:H25, O118:H16, O91:H21, O146:H21, O45:H2, O128:H2, and O121:H19.

vanG element insertions within a conserved chromosomal site conferring vancomycin resistance to Streptococcus agalactiae and Streptococcus anginosus.

Three vancomycin-resistant streptococcal strains carrying vanG elements (two invasive Streptococcus agalactiae isolates [GBS-NY and GBS-NM, both serotype II and multilocus sequence type 22] and one Streptococcus anginosus [Sa]) were examined. The 45,585-bp elements found within Sa and GBS-NY were nearly identical (together designated vanG-1) and shared near-identity over an ~15-kb overlap with a previously described vanG element from Enterococcus faecalis. Unexpectedly, vanG-1 shared much less homology with the 49,321-bp vanG-2 element from GBS-NM, with widely different levels (50% to 99%) of sequence identity shared among 44 related open reading frames. Immediately adjacent to both vanG-1 and vanG-2 were 44,670-bp and 44,680-bp integrative conjugative element (ICE)-like sequences, designated ICE-r, that were nearly identical in the two group B streptococcal (GBS) strains. The dual vanG and ICE-r elements from both GBS strains were inserted at the same position, between bases 1328 and 1329, within the identical RNA methyltransferase (rumA) genes. A GenBank search revealed that although most GBS strains contained insertions within this specific site, only sequence type 22 (ST22) GBS strains contained highly related ICE-r derivatives. The vanG-1 element in Sa was also inserted within this position corresponding to its rumA homolog adjacent to an ICE-r derivative. vanG-1 insertions were previously reported within the same relative position in the E. faecalis rumA homolog. An ICE-r sequence perfectly conserved with respect to its counterpart in GBS-NY was apparent within the same site of the rumA homolog of a Streptococcus dysgalactiae subsp. equisimilis strain. Additionally, homologous vanG-like elements within the conserved rumA target site were evident in Roseburia intestinalis. Importance: These three streptococcal strains represent the first known vancomycin-resistant strains of their species. The collective observations made from these strains reveal a specific hot spot for insertional elements that is conserved between streptococci and different Gram-positive species. The two GBS strains potentially represent a GBS lineage that is predisposed to insertion of vanG elements.