Microbial Genomics - Online First

Research into the lower urinary tract (LUT) microbiota has primarily focused on its relationship to LUT symptoms (LUTS), taking snapshots of these communities in individuals with and without LUTS. While certain bacterial taxa have been associated with LUTS, or the lack thereof, the temporal dynamics of this community were largely unknown. Recently, we conducted a longitudinal study and found that vaginal intercourse resulted in a shift in species richness and diversity within the LUT microbiota. This is particularly relevant as frequent vaginal intercourse is a major risk factor for urinary tract infection (UTI) in premenopausal women (Aydin et al. Int Urogynecol J 2015;26:795–804). To further investigate the relationship between vaginal intercourse and LUT microbiota, here we present the results of a 3 week study in which daily urogenital specimens were collected from a female participant and her male sexual partner. Consistent with our previous findings, the LUT microbiota changed after vaginal intercourse, most notably a high abundance of Streptococcus mitis was observed post-coitus. We isolated and sequenced S. mitis from both sexual partners finding that: (i) the S. mitis isolates from the female partner’s urogenital tract were genomically similar throughout the duration of the study, and (ii) they were related to one isolate from the male partner’s oral cavity collected at the end of the study, suggesting transmission between the two individuals. We hypothesize that blooms in S. mitis after vaginal intercourse may play a role in coitus-related UTI. We found that a S. mitis isolate, in contrast to a Lactobacillus jensenii isolate displaced after vaginal intercourse, cannot inhibit the growth of uropathogenic Escherichia coli . Thus, this bloom in S. mitis may provide a window of opportunity for a uropathogen to colonize the LUT.

In 2018, a cluster of two cases of cryptococcosis occurred at the Queen Elizabeth University Hospital (QEUH) in Glasgow, Scotland (UK). It was postulated that these cases may have been linked to pigeon droppings found on the hospital site, given there have been previous reports of Cryptococcus neoformans associated with pigeon guano. Although some samples of pigeon guano taken from the site yielded culturable yeast from genera related to Cryptococcus, they have since been classified as Naganishia or Papiliotrema spp., and no isolates of C. neoformans were recovered from either the guano or subsequent widespread air sampling. In an attempt to further elucidate any possible shared source of the clinical isolates, we used whole-genome sequencing and phylogenetic analysis to examine the relationship of the two Cryptococcus isolates from the QEUH cases, along with two isolates from sporadic cases treated at a different Glasgow hospital earlier in 2018. Our work demonstrated that these four clinical isolates were not clonally related; while all isolates were from the VNI global lineage and of the same mating type (MATα), the genotypes of the two QEUH isolates were separated by 1885 base changes and belonged to different sub-lineages, recently described as the intercontinental sub-clades VNIa-93 and VNIa-5. In contrast, one of the two sporadic 2018 clinical isolates was determined to belong to the VNIb sub-lineage and the other classified as a VNIV/VNI hybrid. Our work demonstrated that the two 2018 QEUH isolates and the two prior C. neoformans clinical isolates were all genetically distinct. It was not possible to determine whether the QEUH genotypes stemmed from independent sources or from the same source, i.e. pigeons carrying different genotypes, but it should be noted that whilst members of allied genera within the Tremellomycetes were isolated from the hospital environment, there were no environmental isolations of C. neoformans.

Background. Infections caused by carbapenem-resistant Acinetobacter baumannii (CR-Ab) have become increasingly prevalent in clinical settings and often result in significant morbidity and mortality due to their multidrug resistance (MDR). Here we present an integrated whole-genome sequencing (WGS) response to a persistent CR-Ab outbreak in a Brisbane hospital between 2016–2018.

Methods. A. baumannii, Klebsiella pneumoniae, Serratia marcescens and Pseudomonas aeruginosa isolates were sequenced using the Illumina platform primarily to establish isolate relationships based on core-genome SNPs, MLST and antimicrobial resistance gene profiles. Representative isolates were selected for PacBio sequencing. Environmental metagenomic sequencing with Illumina was used to detect persistence of the outbreak strain in the hospital.

Results. In response to a suspected polymicrobial outbreak between May to August of 2016, 28 CR-Ab (and 21 other MDR Gram-negative bacilli) were collected from Intensive Care Unit and Burns Unit patients and sent for WGS with a 7 day turn-around time in clinical reporting. All CR-Ab were sequence type (ST)1050 (Pasteur ST2) and within 10 SNPs apart, indicative of an ongoing outbreak, and distinct from historical CR-Ab isolates from the same hospital. Possible transmission routes between patients were identified on the basis of CR-Ab and K. pneumoniae SNP profiles. Continued WGS surveillance between 2016 to 2018 enabled suspected outbreak cases to be refuted, but a resurgence of the outbreak CR-Ab mid-2018 in the Burns Unit prompted additional screening. Environmental metagenomic sequencing identified the hospital plumbing as a potential source. Replacement of the plumbing and routine drain maintenance resulted in rapid resolution of the secondary outbreak and significant risk reduction with no discernable transmission in the Burns Unit since.

Conclusion. We implemented a comprehensive WGS and metagenomics investigation that resolved a persistent CR-Ab outbreak in a critical care setting.

Although RNA helicases are essentially ubiquitous and perform roles in all stages of RNA metabolism, phylogenetic analysis of the DEAD (Asp-Glu-Ala-Asp)-box RNA helicase family in a single phylum has not been performed. Here, we performed a phylogenetic analysis on DEAD-box helicases from all currently available cyanobacterial genomes, comprising a total of 362 helicase protein sequences from 280 strains. DEAD-box helicases belonging to three distinct clades were observed. Two clades, the CsdA (cold shock DEAD-box A)-like and RhlE (RNA helicase E)-like helicases, cluster with the homologous proteins from Escherichia coli . The third clade, the CrhR (cyanobacterial RNA helicase Redox)-like helicases, is unique to cyanobacteria and characterized by a conserved sequence motif in the C-terminal extension. Restricted distribution is observed across cyanobacterial diversity with respect to both helicase type and strain. CrhR-like and CsdA-like helicases essentially never occur together, while RhlE always occurs with either a CrhR-like or CsdA-like helicase. CrhR-like and RhlE-like proteins occurred in filamentous cyanobacteria of the orders Nostocales , Oscillatoriales and Synechococcales . Similarly, CsdA- and RhlE-like proteins are restricted to unicellular cyanobacteria of the genera Cyanobium and Synechococcus . In addition, the unexpected occurrence of RhlE in two Synechococcus strains suggests recent acquisition and evolutionary divergence. This study, therefore, raises physiological and evolutionary questions as to why DEAD-box RNA helicases encoded in cyanobacterial lineages display restricted distributions, suggesting niches that require either CrhR or CsdA RNA helicase activity but not both. Extensive conservation of gene synteny surrounding the previously described rimO–crhR operon is also observed, indicating a role in the maintenance of photosynthesis. The analysis provides insights into the evolution, origin and dissemination of sequences within a single gene family to yield divergent functional roles.

Whole-genome sequencing (WGS) is fundamental to Mycobacterium tuberculosis basic research and many clinical applications. Coverage across Illumina-sequenced M. tuberculosis genomes is known to vary with sequence context, but this bias is poorly characterized. Here, through a novel application of phylogenomics that distinguishes genuine coverage bias from deletions, we discern Illumina ‘blind spots’ in the M. tuberculosis reference genome for seven sequencing workflows. We find blind spots to be widespread, affecting 529 genes, and provide their exact coordinates, enabling salvage of unaffected regions. Fifty-seven pe/ppe genes (the primary families assumed to exhibit Illumina bias) lack blind spots entirely, while the remaining pe/ppe genes account for 55.1 % of blind spots. Surprisingly, we find coverage bias persists in homopolymers as short as 6 bp, shorter tracts than previously reported. While G+C-rich regions challenge all Illumina sequencing workflows, a modified Nextera library preparation that amplifies DNA with a high-fidelity polymerase markedly attenuates coverage bias in G+C-rich and homopolymeric sequences, expanding the ‘Illumina-sequenceable’ genome. Through these findings, and by defining workflow-specific exclusion criteria, we spotlight effective strategies for handling bias in M. tuberculosis Illumina WGS. This empirical analysis framework may be used to systematically evaluate coverage bias in other species using existing sequencing data.

Species belonging to the family Lactobacillaceae are found in highly diverse environments and play an important role in fermented foods and probiotic products. Many of these species have been individually reported to harbour plasmids that encode important genes. In this study, we performed comparative genomic analysis of publicly available data for 512 plasmids from 282 strains represented by 51 species of this family and correlated the genomic features of plasmids with the ecological niches in which these species are found. Two-thirds of the species had at least one plasmid-harbouring strain. Plasmid abundance and GC content were significantly lower in vertebrate-adapted species as compared to nomadic and free-living species. Hierarchical clustering highlighted the distinct nature of plasmids from the nomadic and free-living species than those from the vertebrate-adapted species. EggNOG-assisted functional annotation revealed that genes associated with transposition, conjugation, DNA repair and recombination, exopolysaccharide production, metal ion transport, toxin–antitoxin system, and stress tolerance were significantly enriched on the plasmids of the nomadic and in some cases nomadic and free-living species. On the other hand, genes related to anaerobic metabolism, ABC transporters and the major facilitator superfamily were overrepresented on the plasmids of the vertebrate-adapted species. These genomic signatures correlate with the comparatively nutrient-depleted, stressful and dynamic environments of nomadic and free-living species and nutrient-rich and anaerobic environments of vertebrate-adapted species. Thus, these results indicate the contribution of the plasmids in the adaptation of lactobacilli to their respective habitats. This study also underlines the potential application of these plasmids in improving the technological and probiotic properties of lactic acid bacteria.

Salmonella enterica serovar Typhimurium is the leading cause of salmonellosis in Australia, and the ability to identify outbreaks and their sources is vital to public health. Here, we examined the utility of whole-genome sequencing (WGS), including complete genome sequencing with Oxford Nanopore technologies, in examining 105 isolates from an endemic multi-locus variable number tandem repeat analysis (MLVA) type over 5 years. The MLVA type was very homogeneous, with 90 % of the isolates falling into groups with a five SNP cut-off. We developed a new two-step approach for outbreak detection using WGS. The first clustering at a zero single nucleotide polymorphism (SNP) cut-off was used to detect outbreak clusters that each occurred within a 4 week window and then a second clustering with dynamically increased SNP cut-offs were used to generate outbreak investigation clusters capable of identifying all outbreak cases. This approach offered optimal specificity and sensitivity for outbreak detection and investigation, in particular of those caused by endemic MLVA types or clones with low genetic diversity. We further showed that inclusion of complete genome sequences detected no additional mutational events for genomic outbreak surveillance. Phylogenetic analysis found that the MLVA type was likely to have been derived recently from a single source that persisted over 5 years, and seeded numerous sporadic infections and outbreaks. Our findings suggest that SNP cut-offs for outbreak cluster detection and public-health surveillance should be based on the local diversity of the relevant strains over time. These findings have general applicability to outbreak detection of bacterial pathogens.