Screening and Genomic Characterization of Filamentous Hemagglutinin-Deficient Bordetella pertussis.

Despite high vaccine coverage, pertussis cases in the United States have increased over the last decade. Growing evidence suggests that disease resurgence results, in part, from genetic divergence of circulating strain populations away from vaccine references. The United States employs acellular vaccines exclusively, and current Bordetella pertussis isolates are predominantly deficient in at least one immunogen, pertactin (Prn). First detected in the United States retrospectively in a 1994 isolate, the rapid spread of Prn deficiency is likely vaccine driven, raising concerns about whether other acellular vaccine immunogens experience similar pressures, as further antigenic changes could potentially threaten vaccine efficacy. We developed an electrochemiluminescent antibody capture assay to monitor the production of the acellular vaccine immunogen filamentous hemagglutinin (Fha). Screening 722 U.S. surveillance isolates collected from 2010 to 2016 identified two that were both Prn and Fha deficient. Three additional Fha-deficient laboratory strains were also identified from a historic collection of 65 isolates dating back to 1935. Whole-genome sequencing of deficient isolates revealed putative, underlying genetic changes. Only four isolates harbored mutations to known genes involved in Fha production, highlighting the complexity of its regulation. The chromosomes of two Fha-deficient isolates included unexpected structural variation that did not appear to influence Fha production. Furthermore, insertion sequence disruption of fhaB was also detected in a previously identified pertussis toxin-deficient isolate that still produced normal levels of Fha. These results demonstrate the genetic potential for additional vaccine immunogen deficiency and underscore the importance of continued surveillance of circulating B. pertussis evolution in response to vaccine pressure.

Revisiting the taxonomy of the genus Elizabethkingia using whole-genome sequencing, optical mapping, and MALDI-TOF, along with proposal of three novel Elizabethkingia species: Elizabethkingia bruuniana sp. nov., Elizabethkingia ursingii sp. nov., and Elizabethkingia occulta sp. nov.

The genus Elizabethkingia is genetically heterogeneous, and the phenotypic similarities between recognized species pose challenges in correct identification of clinically derived isolates. In addition to the type species Elizabethkingia meningoseptica, and more recently proposed Elizabethkingia miricola, Elizabethkingia anophelis and Elizabethkingia endophytica, four genomospecies have long been recognized. By comparing historic DNA-DNA hybridization results with whole genome sequences, optical maps, and MALDI-TOF mass spectra on a large and diverse set of strains, we propose a comprehensive taxonomic revision of this genus. Genomospecies 1 and 2 contain the type strains E. anophelis and E. miricola, respectively. Genomospecies 3 and 4 are herein proposed as novel species named as Elizabethkingia bruuniana sp. nov. (type strain, G0146T = DSM 2975T = CCUG 69503T = CIP 111191T) and Elizabethkingia ursingii sp. nov. (type strain, G4122T = DSM 2974T = CCUG 69496T = CIP 111192T), respectively. Finally, the new species Elizabethkingia occulta sp. nov. (type strain G4070T = DSM 2976T = CCUG 69505T = CIP 111193T), is proposed.

The History of Bordetella pertussis Genome Evolution Includes Structural Rearrangement.

Despite high pertussis vaccine coverage, reported cases of whooping cough (pertussis) have increased over the last decade in the United States and other developed countries. Although Bordetella pertussis is well known for its limited gene sequence variation, recent advances in long-read sequencing technology have begun to reveal genomic structural heterogeneity among otherwise indistinguishable isolates, even within geographically or temporally defined epidemics. We have compared rearrangements among complete genome assemblies from 257 B. pertussis isolates to examine the potential evolution of the chromosomal structure in a pathogen with minimal gene nucleotide sequence diversity. Discrete changes in gene order were identified that differentiated genomes from vaccine reference strains and clinical isolates of various genotypes, frequently along phylogenetic boundaries defined by single nucleotide polymorphisms. The observed rearrangements were primarily large inversions centered on the replication origin or terminus and flanked by IS481, a mobile genetic element with >240 copies per genome and previously suspected to mediate rearrangements and deletions by homologous recombination. These data illustrate that structural genome evolution in B. pertussis is not limited to reduction but also includes rearrangement. Therefore, although genomes of clinical isolates are structurally diverse, specific changes in gene order are conserved, perhaps due to positive selection, providing novel information for investigating disease resurgence and molecular epidemiology.IMPORTANCE Whooping cough, primarily caused by Bordetella pertussis, has resurged in the United States even though the coverage with pertussis-containing vaccines remains high. The rise in reported cases has included increased disease rates among all vaccinated age groups, provoking questions about the pathogen's evolution. The chromosome of B. pertussis includes a large number of repetitive mobile genetic elements that obstruct genome analysis. However, these mobile elements facilitate large rearrangements that alter the order and orientation of essential protein-encoding genes, which otherwise exhibit little nucleotide sequence diversity. By comparing the complete genome assemblies from 257 isolates, we show that specific rearrangements have been conserved throughout recent evolutionary history, perhaps by eliciting changes in gene expression, which may also provide useful information for molecular epidemiology.

Evaluation of a TaqMan Array Card for Detection of Central Nervous System Infections.

Infections of the central nervous system (CNS) are often acute, with significant morbidity and mortality. Routine diagnosis of such infections is limited in developing countries and requires modern equipment in advanced laboratories that may be unavailable to a number of patients in sub-Saharan Africa. We developed a TaqMan array card (TAC) that detects multiple pathogens simultaneously from cerebrospinal fluid. The 21-pathogen CNS multiple-pathogen TAC (CNS-TAC) assay includes two parasites (Balamuthia mandrillaris and Acanthamoeba), six bacterial pathogens (Streptococcus pneumoniae, Haemophilus influenzae, Neisseria meningitidis, Mycoplasma pneumoniae, Mycobacterium tuberculosis, and Bartonella), and 13 viruses (parechovirus, dengue virus, Nipah virus, varicella-zoster virus, mumps virus, measles virus, lyssavirus, herpes simplex viruses 1 and 2, Epstein-Barr virus, enterovirus, cytomegalovirus, and chikungunya virus). The card also includes human RNase P as a nucleic acid extraction control and an internal manufacturer control, GAPDH (glyceraldehyde-3-phosphate dehydrogenase). This CNS-TAC assay can test up to eight samples for all 21 agents within 2.5 h following nucleic acid extraction. The assay was validated for linearity, limit of detection, sensitivity, and specificity by using either live viruses (dengue, mumps, and measles viruses) or nucleic acid material (Nipah and chikungunya viruses). Of 120 samples tested by individual real-time PCR, 35 were positive for eight different targets, whereas the CNS-TAC assay detected 37 positive samples across nine different targets. The CNS-TAC assays showed 85.6% sensitivity and 96.7% specificity. Therefore, the CNS-TAC assay may be useful for outbreak investigation and surveillance of suspected neurological disease.

Enterococcus crotali sp. nov., isolated from faecal material of a timber rattlesnake.

A facultatively anaerobic, Gram-stain-positive bacterium, designated ETRF1T, was found in faecal material of a timber rattlesnake (Crotalus horridus). Based on a comparative 16S rRNA gene sequence analysis, the isolate was assigned to the genus Enterococcus. The 16S rRNA gene sequence of strain ETRF1T showed >97 % similarity to that of the type strains of Enterococcus rotai, E. caccae, E. silesiacus, E haemoperoxidus, E. ureasiticus, E. moraviensis, E. plantarum, E. quebecensis, E. ureilyticus, E. termitis, E. rivorum and E. faecalis. The organism could be distinguished from these 12 phylogenetically related enterococci using conventional biochemical testing, the Rapid ID32 Strep system, comparative pheS and rpoA gene sequence analysis, and comparative whole genome sequence analysis. The estimated in silico DNA-DNA hybridization values were <70 %, and average nucleotide identity values were <96 %, when compared to these 12 species, further validating that ETRF1T represents a unique species within the genus Enterococcus. On the basis of these analyses, strain ETRF1T (=CCUG 65857T=LMG 28312T) is proposed as the type strain of a novel species, Enterococcus crotali sp. nov.

Molecular detection of airborne Coccidioides in Tucson, Arizona.

Environmental surveillance of the soil-dwelling fungus Coccidioides is essential for the prevention of Valley fever, a disease primarily caused by inhalation of the arthroconidia. Methods for collecting and detecting Coccidioides in soil samples are currently in use by several laboratories; however, a method utilizing current air sampling technologies has not been formally demonstrated for the capture of airborne arthroconidia. In this study, we collected air/dust samples at two sites (Site A and Site B) in the endemic region of Tucson, Arizona, and tested a variety of air samplers and membrane matrices. We then employed a single-tube nested qPCR assay for molecular detection. At both sites, numerous soil samples (n = 10 at Site A and n = 24 at Site B) were collected and Coccidioides was detected in two samples (20%) at Site A and in eight samples (33%) at Site B. Of the 25 air/dust samples collected at both sites using five different air sampling methods, we detected Coccidioides in three samples from site B. All three samples were collected using a high-volume sampler with glass-fiber filters. In this report, we describe these methods and propose the use of these air sampling and molecular detection strategies for environmental surveillance of Coccidioides.

Campylobacter fetus subsp. testudinum subsp. nov., isolated from humans and reptiles.

A polyphasic study was undertaken to determine the taxonomic position of 13 Campylobacter fetus-like strains from humans (n = 8) and reptiles (n = 5). The results of matrix-assisted laser desorption ionization time-of-flight (MALDI-TOF) MS and genomic data from sap analysis, 16S rRNA gene and hsp60 sequence comparison, pulsed-field gel electrophoresis, amplified fragment length polymorphism analysis, DNA-DNA hybridization and whole genome sequencing demonstrated that these strains are closely related to C. fetus but clearly differentiated from recognized subspecies of C. fetus. Therefore, this unique cluster of 13 strains represents a novel subspecies within the species C. fetus, for which the name Campylobacter fetus subsp. testudinum subsp. nov. is proposed, with strain 03-427(T) ( = ATCC BAA-2539(T) = LMG 27499(T)) as the type strain. Although this novel taxon could not be differentiated from C. fetus subsp. fetus and C. fetus subsp. venerealis using conventional phenotypic tests, MALDI-TOF MS revealed the presence of multiple phenotypic biomarkers which distinguish Campylobacter fetus subsp. testudinum subsp. nov. from recognized subspecies of C. fetus.

The first complete genome of the simian malaria parasite Plasmodium brasilianum.

Naturally occurring human infections by zoonotic Plasmodium species have been documented for P. knowlesi, P. cynomolgi, P. simium, P. simiovale, P. inui, P. inui-like, P. coatneyi, and P. brasilianum. Accurate detection of each species is complicated by their morphological similarities with other Plasmodium species. PCR-based assays offer a solution but require prior knowledge of adequate genomic targets that can distinguish the species. While whole genomes have been published for P. knowlesi, P. cynomolgi, P. simium, and P. inui, no complete genome for P. brasilianum has been available. Previously, we reported a draft genome for P. brasilianum, and here we report the completed genome for P. brasilianum. The genome is 31.4 Mb in size and comprises 14 chromosomes, the mitochondrial genome, the apicoplast genome, and 29 unplaced contigs. The chromosomes consist of 98.4% nucleotide sites that are identical to the P. malariae genome, the closest evolutionarily related species hypothesized to be the same species as P. brasilianum, with 41,125 non-synonymous SNPs (0.0722% of genome) identified between the two genomes. Furthermore, P. brasilianum had 4864 (82.1%) genes that share 80% or higher sequence similarity with 4970 (75.5%) P. malariae genes. This was demonstrated by the nearly identical genomic organization and multiple sequence alignments for the merozoite surface proteins msp3 and msp7. We observed a distinction in the repeat lengths of the circumsporozoite protein (CSP) gene sequences between P. brasilianum and P. malariae. Our results demonstrate a 97.3% pairwise identity between the P. brasilianum and the P. malariae genomes. These findings highlight the phylogenetic proximity of these two species, suggesting that P. malariae and P. brasilianum are strains of the same species, but this could not be fully evaluated with only a single genomic sequence for each species.

Draft Chromosome Sequences of a Clinical Isolate of the Free-Living Ameba Naegleria fowleri.

We present the chromosome sequences of a Naegleria fowleri isolate from a human primary amebic meningoencephalitis (PAM) case. The genome sequences were assembled from Illumina HiSeq and PacBio sequencing data and verified with the optical mapping data. This led to the identification of 37 contigs representing 37 chromosomes in N. fowleri.

Evaluation of laboratory methods for diagnosis of varicella.

BACKGROUND: The incidence of varicella disease is declining as a result of vaccination, making clinical diagnosis more challenging, particularly for vaccine-modified cases. We conducted a comprehensive evaluation of laboratory tests and specimen types to assess diagnostic performance and determine what role testing can play after skin lesions have resolved. METHODS: We enrolled patients with suspected varicella disease in 2 communities. Enrollees were visited at the time of rash onset and 2 weeks later. Multiple skin lesion, oral, urine, and blood or serum specimens were requested at each visit and tested for varicella zoster virus (VZV) immunoglobulin (Ig) G, IgM, and IgA antibody by enzyme-linked immunoassay; for VZV antigen by direct fluorescent antibody; and/or for VZV DNA by polymerase chain reaction (PCR). Clinical certainty of the diagnosis of varicella disease was scored. PCR results from first-visit vesicles or scab specimens served as the gold standard in assessing test performance. RESULTS: Of 93 enrollees, 53 were confirmed to have varicella disease. Among 20 unmodified cases, PCR testing was 95%-100% sensitive for macular and/or papular lesions and for oral specimens collected at the first visit; most specimens from the second visit yielded negative results. Among 27 vaccine-modified cases, macular and/or papular lesions collected at the first visit were also 100% sensitive; yields from other specimens were poorer, and few specimens from the second visit tested positive. Clinical diagnosis was 100% and 85% sensitive for diagnosing unmodified and vaccine-modified varicella cases, respectively. CONCLUSIONS: PCR testing of skin lesion specimens remains convenient and accurate for diagnosing varicella disease in vaccinated and unvaccinated persons. PCR of oral specimens can sometimes aid in diagnosis of varicella disease, even after rash resolves.

Rapid Inactivation of Non-Endospore-Forming Bacterial Pathogens by Heat Stabilization is Compatible with Downstream Next-Generation Sequencing.

Introduction: Heat stabilization treatment preserves the in vivo state of biological samples by rapidly inactivating enzymes that cause degradation of proteins and nucleic acids. Historically, proteomics studies used this technique as an alternative to chemical fixation. More recently, microbiologists discovered that heat stabilization treatment rapidly inactivates pathogens present in tissue samples and preserves deoxyribonucleic acid (DNA) in the tissue. However, these recent studies did not investigate the inactivation of high-density bacterial suspensions and the quality of bacterial DNA. Methods and Results: High-density suspensions of Escherichia coli (>109 cfu/mL) were completely inactivated by heat stabilization treatment using the Denator Stabilizor T1 instrument at 72°C and 95°C for 45 seconds. Using the heat stabilization instrument, a panel of 30 species, 20 Gram-negative and 10 non-endospore-forming Gram-positive species, were fully inactivated by treatment (95°C for 45 seconds). DNA was isolated from bacterial suspensions of Gram-negative bacteria, including E. albertii, E. coli, Shigella dysenteriae, and S. flexneri, following inactivation via heat stabilization treatment and without treatment. DNA isolated following heat stabilization treatment was fully compatible with all downstream molecular applications tested, including next-generation sequencing, pulsed-field gel electrophoresis, multiplex polymerase chain reaction (PCR), and real-time PCR. Conclusions and Discussion: Heat stabilization treatment of Gram-negative and non-endospore-forming Gram-positive pathogens completely inactivates high-density bacterial suspensions. This treatment is compatible with downstream DNA molecular assays, including next-generation sequencing, pulsed-field gel electrophoresis, and PCR. Inactivation by heat stabilization is a rapid process that may increase safety by decreasing risks for laboratory-associated infections and risks associated with transportation of infectious materials.

Draft Genome Sequences of Leishmania (Leishmania) amazonensis, Leishmania (Leishmania) mexicana, and Leishmania (Leishmania) aethiopica, Potential Etiological Agents of Diffuse Cutaneous Leishmaniasis.

We present here the draft genome sequences of Leishmania (Leishmania) amazonensis, Leishmania (Leishmania) mexicana, and Leishmania (Leishmania) aethiopica, potential etiological agents of diffuse cutaneous leishmaniasis (DCL). Sequence data were obtained using PacBio and MiSeq platforms. The PacBio assemblies generated using Canu v1.6 are more contiguous than are those in the available data.

First Draft Genome Sequence of Leishmania (Viannia) lainsoni Strain 216-34, Isolated from a Peruvian Clinical Case.

We present here the first draft genome sequence of Leishmania (Viannia) lainsoni strain 216-34, sequenced using PacBio and MiSeq platforms. PacBio contigs were generated from de novo assemblies using CANU version 1.6 and polished using Illumina reads.

Conserved Patterns of Symmetric Inversion in the Genome Evolution of Bordetella Respiratory Pathogens.

Whooping cough (pertussis), primarily caused by Bordetella pertussis, has resurged in the United States, and circulating strains exhibit considerable chromosome structural fluidity in the form of rearrangement and deletion. The genus Bordetella includes additional pathogenic species infecting various animals, some even causing pertussis-like respiratory disease in humans; however, investigation of their genome evolution has been limited. We studied chromosome structure in complete genome sequences from 167 Bordetella species isolates, as well as 469 B. pertussis isolates, to gain a generalized understanding of rearrangement patterns among these related pathogens. Observed changes in gene order primarily resulted from large inversions and were only detected in species with genomes harboring multicopy insertion sequence (IS) elements, most notably B. holmesii and B. parapertussis While genomes of B. pertussis contain >240 copies of IS481, IS elements appear less numerous in other species and yield less chromosome structural diversity through rearrangement. These data were further used to predict all possible rearrangements between IS element copies present in Bordetella genomes, revealing that only a subset is observed among circulating strains. Therefore, while it appears that rearrangement occurs less frequently in other species than in B. pertussis, these clinically relevant respiratory pathogens likely experience similar mutation of gene order. The resulting chromosome structural fluidity presents both challenges and opportunity for the study of Bordetella respiratory pathogens.IMPORTANCE Bordetella pertussis is the primary agent of whooping cough (pertussis). The Bordetella genus includes additional pathogens of animals and humans, including some that cause pertussis-like respiratory illness. The chromosome of B. pertussis has previously been shown to exhibit considerable structural rearrangement, but insufficient data have prevented comparable investigation in related species. In this study, we analyze chromosome structure variation in several Bordetella species to gain a generalized understanding of rearrangement patterns in this genus. Just as in B. pertussis, we observed inversions in other species that likely result from common mutational processes. We used these data to further predict additional, unobserved inversions, suggesting that specific genome structures may be preferred in each species.

Identification of five major and two minor genotypes of varicella-zoster virus strains: a practical two-amplicon approach used to genotype clinical isolates in Australia and New Zealand.

Whole genome phylogenetic analysis in this study resolved a total of five major genotypes among the 22 varicella-zoster virus (VZV) strains or isolates for which complete genomic sequences are available. Consistent with earlier publications we have designated these genotypes European 1 (E1), European 2 (E2), Japanese (J), mosaic 1 (M1), and mosaic 2 (M2). Single nucleotide polymorphism (SNP) analysis performed in a whole-genome alignment revealed that VZV isolates of all five genotypes can be accurately genotyped using SNPs from two amplicons: open reading frame 22 (ORF22) and either ORF21 or ORF50. This modified approach identifies all of the genotypes observed using any of the published genotyping protocols. Of 165 clinical varicella and zoster isolates from Australia and New Zealand typed using this approach, 67 of 127 eastern Australian isolates were E1, 30 were E2, 16 were J, 10 were M1, and 4 were M2; 25 of 38 New Zealand isolates were E1, 8 were E2, and 5 were M1. VZV strain diversity in eastern Australia is thus broader than has been described for any other region, including Europe, Africa, and North America. J strains were far more prevalent than previously observed in countries other than Japan. Two-amplicon typing was in complete accord with genotypes derived using SNP in multiple ORFs (ORFs 1, 21, 22, 38, 50, 54, and 62). Two additional minor genotypes, M3 and M4, could also be resolved using two-amplicon typing.

Draft Genome Sequence of French Guiana Leishmania (Viannia) guyanensis Strain 204-365, Assembled Using Long Reads.

We present here the draft genome sequence for Leishmania (Viannia) guyanensis. The isolate was obtained from a clinical case of cutaneous leishmaniasis in French Guiana. Genomic DNA was sequenced using PacBio and MiSeq platforms.

Genome Sequence of the Amphotericin B-Resistant Candida duobushaemulonii Strain B09383.

Candida duobushaemulonii is a drug-resistant yeast that can cause invasive candidiasis. Here, we report the first genome sequence of C. duobushaemulonii, isolate B09383, generated using PacBio sequencing technology. The estimated genome size was 12.5 Mb with a GC content of 46.84%.

PacBio Genome Sequences of Escherichia coli Serotype O157:H7, Diffusely Adherent E. coli, and Salmonella enterica Strains, All Carrying Plasmids with an mcr-1 Resistance Gene.

We report here Illumina-corrected PacBio whole-genome sequences of an Escherichia coli serotype O157:H7 strain (2017C-4109), an E. coli serotype O[undetermined]:H2 strain (2017C-4173W12), and a Salmonella enterica subsp. enterica serovar Enteritidis strain (2017K-0021), all of which carried the mcr-1 resistance gene on an IncI2 or IncX4 plasmid. We also determined that pMCR-1-CTSe is identical to a previously published plasmid, pMCR-1-CT.

Genomic insights into multidrug-resistance, mating and virulence in Candida auris and related emerging species.

Candida auris is an emergent multidrug-resistant fungal pathogen causing increasing reports of outbreaks. While distantly related to C. albicans and C. glabrata, C. auris is closely related to rarely observed and often multidrug-resistant species from the C. haemulonii clade. Here, we analyze near complete genome assemblies for the four C. auris clades and three related species, and map intra- and inter-species rearrangements across the seven chromosomes. Using RNA-Seq-guided gene predictions, we find that most mating and meiosis genes are conserved and that clades contain either the MTLa or MTLα mating loci. Comparing the genomes of these emerging species to those of other Candida species identifies genes linked to drug resistance and virulence, including expanded families of transporters and lipases, as well as mutations and copy number variants in ERG11. Gene expression analysis identifies transporters and metabolic regulators specific to C. auris and those conserved with related species which may contribute to differences in drug response in this emerging fungal clade.

High-Quality Complete and Draft Genome Sequences for Three Escherichia spp. and Three Shigella spp. Generated with Pacific Biosciences and Illumina Sequencing and Optical Mapping.

Escherichia spp., including E. albertii and E. coli, Shigella dysenteriae, and S. flexneri are causative agents of foodborne disease. We report here reference-level whole-genome sequences of E. albertii (2014C-4356), E. coli (2011C-4315 and 2012C-4431), S. dysenteriae (BU53M1), and S. flexneri (94-3007 and 71-2783).