- Volume 6, Issue 1, 2020

Tuberculosis (TB) surveillance is scarce in most African countries, even though it is the continent with the greatest disease incidence according to the World Health Organization. Liberia is within the 30 countries with the highest TB burden, probably as a consequence of the long civil war and the recent Ebola outbreak, both crippling the health system and depreciating the TB prevention and control programmes. Due to difficulties working in the country, there is a lack of resistance surveys and bacillus characterization. Here, we use genome sequencing of Mycobacterium tuberculosis clinical isolates to fill this gap. Our results highlight that the bacillus population structure is dominated by lineage 4 strains that harbour an outstanding genetic diversity, higher than in the rest of Africa as a whole. Coalescent analyses demonstrate that strains currently circulating in Liberia were introduced several times beginning in the early year 600 CE until very recently coinciding with migratory movements associated with the civil war and Ebola epidemics. A higher multidrug-resistant (MDR)-TB frequency (23.5 %) than current estimates was obtained together with non-catalogued drug-resistance mutations. Additionally, 39 % of strains were in genomic clusters revealing that ongoing transmission is a major contribution to the TB burden in the country. Our report emphasizes the importance of TB surveillance and control in African countries where bacillus diversity, MDR-TB prevalence and transmission are coalescing to jeopardize TB control programmes.

The majority of bacterial genomes have high coding efficiencies, but there are some genomes of intracellular bacteria that have low gene density. The genome of the endosymbiont Sodalis glossinidius contains almost 50 % pseudogenes containing mutations that putatively silence them at the genomic level. We have applied multiple ‘omic’ strategies, combining Illumina and Pacific Biosciences Single-Molecule Real-Time DNA sequencing and annotation, stranded RNA sequencing and proteome analysis to better understand the transcriptional and translational landscape of Sodalis pseudogenes, and potential mechanisms for their control. Between 53 and 74 % of the Sodalis transcriptome remains active in cell-free culture. The mean sense transcription from coding domain sequences (CDSs) is four times greater than that from pseudogenes. Comparative genomic analysis of six Illumina-sequenced Sodalis isolates from different host Glossina species shows pseudogenes make up ~40 % of the 2729 genes in the core genome, suggesting that they are stable and/or that Sodalis is a recent introduction across the genus Glossina as a facultative symbiont. These data shed further light on the importance of transcriptional and translational control in deciphering host–microbe interactions. The combination of genomics, transcriptomics and proteomics gives a multidimensional perspective for studying prokaryotic genomes with a view to elucidating evolutionary adaptation to novel environmental niches.

Sequence type (ST)73 has emerged as one of the most frequently isolated extraintestinal pathogenic Escherichia coli. To examine the localized diversity of ST73 clonal groups, including their mobile genetic element profile, we sequenced the genomes of 16 multiple-drug resistant ST73 isolates from patients with urinary tract infection from a single hospital in Sydney, Australia, between 2009 and 2011. Genome sequences were used to generate a SNP-based phylogenetic tree to determine the relationship of these isolates in a global context with ST73 sequences (n=210) from public databases. There was no evidence of a dominant outbreak strain of ST73 in patients from this hospital, rather we identified at least eight separate groups, several of which reoccurred, over a 2 year period. The inferred phylogeny of all ST73 strains (n=226) including the ST73 clone D i2 reference genome shows high bootstrap support and clusters into four major groups that correlate with serotype. The Sydney ST73 strains carry a wide variety of virulence-associated genes, but the presence of iss, pic and several iron-acquisition operons was notable.

Carbapenemase-producing Enterobacteriaceae (CPE) are an increasingly common cause of healthcare-associated infections and may occasionally be identified in patients without extensive healthcare exposure. bla IMP-4 is the most frequently detected carbapenemase gene in Enterobacteriaceae within Australia, but little is known about the mechanisms behind its persistence. Here we used whole genome sequencing (WGS) to investigate the molecular epidemiology of bla IMP-4 in Queensland, Australia. In total, 107 CPE were collected between 2014 and 2017 and sent for WGS on an Illumina NextSeq500. Resistance genes and plasmid types were detected using a combination of read mapping and nucleotide comparison of de novo assemblies. Six isolates were additionally sequenced using Oxford Nanopore MinION to generate long-reads and fully characterize the context of the bla IMP-4 gene. Of 107 CPE, 93 carried the bla IMP-4 gene; 74/107 also carried an IncHI2 plasmid, suggesting carriage of the bla IMP-4 gene on an IncHI2 plasmid. Comparison of these isolates to a previously characterized IncHI2 plasmid pMS7884A (isolated from an Enterobacter hormaechei strain in Brisbane) suggested that all isolates carried a similar plasmid. Five of six representative isolates sequenced using Nanopore long-read technology carried IncHI2 plasmids harbouring the bla IMP-4 gene. While the vast majority of isolates represented  E. hormaechei , several other species were also found to carry the IncHI2 plasmid, including Klebsiella species, Escherichia coli and Citrobacter species. Several clonal groups of E. hormaechei were also identified, suggesting that persistence of bla IMP-4 is driven by both presence on a common plasmid and clonal spread of certain E. hormaechei lineages.

Phages and plasmids play important roles in bacterial evolution and diversification. Although many draft genomes have been generated, phage and plasmid genomes are usually fragmented, limiting our understanding of their dynamics. Here, we performed a systematic analysis of 239 draft genomes and 7 complete genomes of Shiga toxin (Stx)-producing Escherichia coli O145:H28, the major virulence factors of which are encoded by prophages (PPs) or plasmids. The results indicated that PPs are more stably maintained than plasmids. A set of ancestrally acquired PPs was well conserved, while various PPs, including Stx phages, were acquired by multiple sublineages. In contrast, gains and losses of a wide range of plasmids have frequently occurred across the O145:H28 lineage, and only the virulence plasmid was well conserved. The different dynamics of PPs and plasmids have differentially impacted the pangenome of O145:H28, with high proportions of PP- and plasmid-associated genes in the variably present and rare gene fractions, respectively. The dynamics of PPs and plasmids have also strongly impacted virulence gene repertoires, such as the highly variable distribution of stx genes and the high conservation of a set of type III secretion effectors, which probably represents the core effectors of O145:H28 and the genes on the virulence plasmid in the entire O145:H28 population. These results provide detailed insights into the dynamics of PPs and plasmids, and show the application of genomic analyses using a large set of draft genomes and appropriately selected complete genomes.

Cryptococcus neoformans is an opportunistic fungal pathogen that at its peak epidemic levels caused an estimated million cases of cryptococcal meningitis per year worldwide. This species can grow in diverse environmental (trees, soil and bird excreta) and host niches (intracellular microenvironments of phagocytes and free-living in host tissues). The genetic basic for adaptation to these different conditions is not well characterized, as most experimental work has relied on a single reference strain of C. neoformans. To identify genes important for yeast infection and disease progression, we profiled the gene expression of seven C. neoformans isolates grown in five representative in vitro environmental and in vivo conditions. We characterized gene expression differences using RNA-Seq (RNA sequencing), comparing clinical and environmental isolates from two of the major lineages of this species, VNI and VNBI. These comparisons highlighted genes showing lineage-specific expression that are enriched in subtelomeric regions and in lineage-specific gene clusters. By contrast, we find few expression differences between clinical and environmental isolates from the same lineage. Gene expression specific to in vivo stages reflects available nutrients and stresses, with an increase in fungal metabolism within macrophages, and an induction of ribosomal and heat-shock gene expression within the subarachnoid space. This study provides the widest view to date of the transcriptome variation of C. neoformans across natural isolates, and provides insights into genes important for in vitro and in vivo growth stages.

As sequencing read length has increased, researchers have quickly adopted longer reads for their experiments. Here, we examine 14 pathogen or host–pathogen differential gene expression data sets to assess whether using longer reads is warranted. A variety of data sets was used to assess what genomic attributes might affect the outcome of differential gene expression analysis including: gene density, operons, gene length, number of introns/exons and intron length. No genome attribute was found to influence the data in principal components analysis, hierarchical clustering with bootstrap support, or regression analyses of pairwise comparisons that were undertaken on the same reads, looking at all combinations of paired and unpaired reads trimmed to 36, 54, 72 and 101 bp. Read pairing had the greatest effect when there was little variation in the samples from different conditions or in their replicates (e.g. little differential gene expression). But overall, 54 and 72 bp reads were typically most similar. Given differences in costs and mapping percentages, we recommend 54 bp reads for organisms with no or few introns and 72 bp reads for all others. In a third of the data sets, read pairing had absolutely no effect, despite paired reads having twice as much data. Therefore, single-end reads seem robust for differential-expression analyses, but in eukaryotes paired-end reads are likely desired to analyse splice variants and should be preferred for data sets that are acquired with the intent to be community resources that might be used in secondary data analyses.