The Deinococcus protease PprI senses DNA damage by directly interacting with single-stranded DNA

Bacteria have evolved various response systems to adapt to environmental stress. A protease-based derepression mechanism in response to DNA damage was characterized in Deinococcus, which is controlled by the specific cleavage of repressor DdrO by metallopeptidase PprI (also called IrrE). Despite the efforts to document the biochemical, physiological, and downstream regulation of PprI-DdrO, the upstream regulatory signal activating this system remains unclear. Here, we show that single-stranded DNA physically interacts with PprI protease, which enhances the PprI-DdrO interactions as well as the DdrO cleavage in a length-dependent manner both in vivo and in vitro. Structures of PprI, in its apo and complexed forms with single-stranded DNA, reveal two DNA-binding interfaces shaping the cleavage site. Moreover, we show that the dynamic monomer-dimer equilibrium of PprI is also important for its cleavage activity. Our data provide evidence that single-stranded DNA could serve as the signal for DNA damage sensing in the metalloprotease/repressor system in bacteria. These results also shed light on the survival and acquired drug resistance of certain bacteria under antimicrobial stress through a SOS-independent pathway.

TcrXY is an acid-sensing two-component transcriptional regulator of Mycobacterium tuberculosis required for persistent infection

The ability of Mycobacterium tuberculosis (Mtb) to persist in the host complicates and prolongs tuberculosis (TB) patient chemotherapy. Here we demonstrate that a neglected two-component system (TCS) of Mtb, TcrXY, is an autoregulated acid-sensing TCS that controls a functionally diverse 70-gene regulon required for bacterial persistence. Characterisation of two representatives of this regulon, Rv3706c and Rv3705A, implicate these genes as key determinants for the survival of Mtb in vivo by serving as important effectors to mitigate redox stress at acidic pH. We show that genetic silencing of the response regulator tcrX using CRISPR interference attenuates the persistence of Mtb during chronic mouse infection and improves treatment with the two front-line anti-TB drugs, rifampicin and isoniazid. We propose that targeting TcrXY signal transduction blocks the ability of Mtb to sense and respond to acid stress, resulting in a disordered program of persistence to render the organism vulnerable to existing TB chemotherapy.

Commensal lifestyle regulated by a negative feedback loop between Arabidopsis ROS and the bacterial T2SS

Despite the plant health-promoting effects of plant microbiota, these assemblages also comprise potentially detrimental microbes. How plant immunity controls its microbiota to promote plant health under these conditions remains largely unknown. We find that commensal bacteria isolated from healthy Arabidopsis plants trigger diverse patterns of reactive oxygen species (ROS) production dependent on the immune receptors and completely on the NADPH oxidase RBOHD that selectively inhibited specific commensals, notably Xanthomonas L148. Through random mutagenesis, we find that L148 gspE, encoding a type II secretion system (T2SS) component, is required for the damaging effects of Xanthomonas L148 on rbohD mutant plants. In planta bacterial transcriptomics reveals that RBOHD suppresses most T2SS gene expression including gspE. L148 colonization protected plants against a bacterial pathogen, when gspE was inhibited by ROS or mutation. Thus, a negative feedback loop between Arabidopsis ROS and the bacterial T2SS tames a potentially detrimental leaf commensal and turns it into a microbe beneficial to the host.

σE of Streptomyces coelicolor can function both as a direct activator or repressor of transcription

σ factors are considered as positive regulators of gene expression. Here we reveal the opposite, inhibitory role of these proteins. We used a combination of molecular biology methods and computational modeling to analyze the regulatory activity of the extracytoplasmic σE factor from Streptomyces coelicolor. The direct activator/repressor function of σE was then explored by experimental analysis of selected promoter regions in vivo. Additionally, the σE interactome was defined. Taken together, the results characterize σE, its regulation, regulon, and suggest its direct inhibitory function (as a repressor) in gene expression, a phenomenon that may be common also to other σ factors and organisms.

Carbapenem-resistant Escherichia coli exhibit diverse spatiotemporal epidemiological characteristics across the globe

Carbapenem-resistant Escherichia coli (CREC) poses a severe global public health risk. This study reveals the worldwide geographic spreading patterns and spatiotemporal distribution characteristics of resistance genes in 7918 CREC isolates belonging to 497 sequence types (ST) and originating from 75 countries. In the last decade, there has been a transition in the prevailing STs from highly virulent ST131 and ST38 to higher antibiotic-resistant ST410 and ST167. The rise of multi-drug resistant strains of CREC carrying plasmids with extended-spectrum beta-lactamase (ESBL) resistance genes could be attributed to three important instances of host-switching events. The spread of CREC was associated with the changing trends in blaNDM-5, blaKPC-2, and blaOXA-48, as well as the plasmids IncFI, IncFII, and IncI. There were intercontinental geographic transfers of major CREC strains. Various crucial transmission hubs and patterns have been identified for ST131 in the United Kingdom, Italy, the United States, and China, ST167 in India, France, Egypt, and the United States, and ST410 in Thailand, Israel, the United Kingdom, France, and the United States. This work is valuable in managing CREC infections and preventing CREC occurrence and transmission inside healthcare settings and among diverse hosts.

CRISPR-Cas-based identification of a sialylated human milk oligosaccharides utilization cluster in the infant gut commensal Bacteroides dorei

The infant gut microbiome is impacted by early-life feeding, as human milk oligosaccharides (HMOs) found in breastmilk cannot be digested by infants and serve as nutrients for their gut bacteria. While the vast majority of HMO-utilization research has focused on Bifidobacterium species, recent studies have suggested additional HMO-utilizers, mostly Bacteroides, yet their utilization mechanism is poorly characterized. Here, we investigate Bacteroides dorei isolates from breastfed-infants and identify that polysaccharide utilization locus (PUL) 33 enables B. dorei to utilize sialylated HMOs. We perform transcriptional profiling and identity upregulated genes when growing on sialylated HMOs. Using CRISPR-Cas12 to knock-out four PUL33 genes, combined with complementation assays, we identify GH33 as the critical gene in PUL33 for sialylated HMO-utilization. This demonstration of an HMO-utilization system by Bacteroides species isolated from infants opens the way to further characterization of additional such systems, to better understand HMO-utilization in the infant gut.

The CRISPR effector Cam1 mediates membrane depolarization for phage defence

Prokaryotic type III CRISPR-Cas systems provide immunity against viruses and plasmids using CRISPR-associated Rossman fold (CARF) protein effectors1-5. Recognition of transcripts of these invaders with sequences that are complementary to CRISPR RNA guides leads to the production of cyclic oligoadenylate second messengers, which bind CARF domains and trigger the activity of an effector domain6,7. Whereas most effectors degrade host and invader nucleic acids, some are predicted to contain transmembrane helices without an enzymatic function. Whether and how these CARF-transmembrane helix fusion proteins facilitate the type III CRISPR-Cas immune response remains unknown. Here we investigate the role of cyclic oligoadenylate-activated membrane protein 1 (Cam1) during type III CRISPR immunity. Structural and biochemical analyses reveal that the CARF domains of a Cam1 dimer bind cyclic tetra-adenylate second messengers. In vivo, Cam1 localizes to the membrane, is predicted to form a tetrameric transmembrane pore, and provides defence against viral infection through the induction of membrane depolarization and growth arrest. These results reveal that CRISPR immunity does not always operate through the degradation of nucleic acids, but is instead mediated via a wider range of cellular responses.

The Helicobacter pylori Genome Project: insights into H. pylori population structure from analysis of a worldwide collection of complete genomes

Helicobacter pylori, a dominant member of the gastric microbiota, shares co-evolutionary history with humans. This has led to the development of genetically distinct H. pylori subpopulations associated with the geographic origin of the host and with differential gastric disease risk. Here, we provide insights into H. pylori population structure as a part of the Helicobacter pylori Genome Project (HpGP), a multi-disciplinary initiative aimed at elucidating H. pylori pathogenesis and identifying new therapeutic targets. We collected 1011 well-characterized clinical strains from 50 countries and generated high-quality genome sequences. We analysed core genome diversity and population structure of the HpGP dataset and 255 worldwide reference genomes to outline the ancestral contribution to Eurasian, African, and American populations. We found evidence of substantial contribution of population hpNorthAsia and subpopulation hspUral in Northern European H. pylori. The genomes of H. pylori isolated from northern and southern Indigenous Americans differed in that bacteria isolated in northern Indigenous communities were more similar to North Asian H. pylori while the southern had higher relatedness to hpEastAsia. Notably, we also found a highly clonal yet geographically dispersed North American subpopulation, which is negative for the cag pathogenicity island, and present in 7% of sequenced US genomes. We expect the HpGP dataset and the corresponding strains to become a major asset for H. pylori genomics.

Post-crotonylation oxidation by a monooxygenase promotes acetyl-CoA synthetase degradation in Streptomyces roseosporus

Protein post-translational modifications (PTMs) with various acyl groups play central roles in Streptomyces. But whether these acyl groups can be further modified, and the influences of these potential modifications on bacterial physiology have not been addressed. Here in Streptomyces roseosporus with rich crotonylation, a luciferase monooxygenase LimB is identified to elaborately regulate the crotonylation level, morphological development and antibiotic production by oxidation on the crotonyl groups of an acetyl-CoA synthetase Acs. This chemical modification on crotonylation leads to Acs degradation via the protease ClpP1/2 pathway and lowered intracellular crotonyl-CoA pool. Thus, we show that acyl groups after PTMs can be further modified, herein named post-PTM modification (PPM), and LimB is a PTM modifier to control the substrate protein turnover for cell development of Streptomyces. These findings expand our understanding of the complexity of chemical modifications on proteins for physiological regulation, and also suggest that PPM would be widespread.

Ongoing diversification of the global fish pathogen Piscirickettsia salmonis through genetic isolation and transposition bursts

The management of bacterial pathogens remains a key challenge of aquaculture. The marine gammaproteobacterium Piscirickettsia salmonis is the etiological agent of piscirickettsiosis and causes multi-systemic infections in different salmon species, resulting in considerable mortality and substantial commercial losses. Here, we elucidate its global diversity, evolution, and selection during human interventions. Our comprehensive analysis of 73 closed, high quality genome sequences covered strains from major outbreaks and was supplemented by an analysis of all P. salmonis 16S rRNA gene sequences and metagenomic reads available in public databases. Genome comparison showed that Piscirickettsia comprises at least three distinct, genetically isolated species of which two showed evidence for continuing speciation. However, at least twice the number of species exist in marine fish or seawater. A hallmark of Piscirickettsia diversification is the unprecedented amount and diversity of transposases which are particularly active in subgroups undergoing rapid speciation and are key to the acquisition of novel genes and to pseudogenization. Several group-specific genes are involved in surface antigen synthesis and may explain the differences in virulence between strains. However, the frequent failure of antibiotic treatment of piscirickettsiosis outbreaks cannot be explained by horizontal acquisition of resistance genes which so far occurred only very rarely. Besides revealing a dynamic diversification of an important pathogen, our study also provides the data for improving its surveillance, predicting the emergence of novel lineages, and adapting aquaculture management, and thereby contributes towards the sustainability of salmon farming.

Soil bacterium manipulates antifungal weapons by sensing intracellular type IVA secretion system effectors of a competitor

Soil beneficial bacteria can effectively inhibit bacterial pathogens by assembling contact-dependent killing weapons, such as the type IVA secretion system (T4ASS). It's not clear whether these antibacterial weapons are involved in biotrophic microbial interactions in soil. Here we showed that an antifungal antibiotic 2,4-DAPG production of the soil bacterium, Pseudomonas protegens can be triggered by another soil bacterium, Lysobacter enzymogenes, via T4ASS by co-culturing on agar plates to mimic cell-to-cell contact. We demonstrated that the induced 2,4-DAPG production of P. protegens is achieved by intracellular detection of the T4ASS effector protein Le1519 translocated from L. enzymogenes. We defined Le1519 as LtaE (Lysobacter T4E triggering antifungal effects), which specifically stimulates the expression of 2,4-DAPG biosynthesis genes in P. protegens, thereby protecting soybean seedlings from infection by the fungus Rhizoctonia solani. We further found that LtaE directly bound to PhlF, a pathway-specific transcriptional repressor of the 2,4-DAPG biosynthesis, then activated the 2,4-DAPG production. Our results highlight a novel pattern of microbial interspecies and interkingdom interactions, providing a unique case for expanding the diversity of soil microbial interactions.

Phylogenetic and ecophysiological novelty of subsurface mercury methylators in mangrove sediments

Mangrove sediment is a crucial component in the global mercury (Hg) cycling and acts as a hotspot for methylmercury (MeHg) production. Early evidence has documented the ubiquity of well-studied Hg methylators in mangrove superficial sediments; however, their diversity and metabolic adaptation in the more anoxic and highly reduced subsurface sediments are lacking. Through MeHg biogeochemical assay and metagenomic sequencing, we found that mangrove subsurface sediments (20-100 cm) showed a less hgcA gene abundance but higher diversity of Hg methylators than superficial sediments (0-20 cm). Regional-scale investigation of mangrove subsurface sediments spanning over 1500 km demonstrated a prevalence and family-level novelty of Hg-methylating microbial lineages (i.e., those affiliated to Anaerolineae, Phycisphaerae, and Desulfobacterales). We proposed the candidate phylum Zixibacteria lineage with sulfate-reducing capacity as a currently understudied Hg methylator across anoxic environments. Unlike other Hg methylators, the Zixibacteria lineage does not use the Wood-Ljungdahl pathway but has unique capabilities of performing methionine synthesis to donate methyl groups. The absence of cobalamin biosynthesis pathway suggests that this Hg-methylating lineage may depend on its syntrophic partners (i.e., Syntrophobacterales members) for energy in subsurface sediments. Our results expand the diversity of subsurface Hg methylators and uncover their unique ecophysiological adaptations in mangrove sediments.

When less is more: shortening the Lpp protein leads to increased vancomycin resistance in Escherichia coli

Vancomycin is a naturally occurring cell-wall-targeting glycopeptide antibiotic. Due to the low potency of this antibiotic against Gram-negative pathogens, such as Escherichia coli, there is a limited knowledge about interactions between vancomycin and this group of bacteria. Here, we show that an in-frame 63 bp deletion of the lpp gene caused a fourfold increase in vancomycin resistance in E. coli. The resulting protein, LppΔ21, is 21 amino acids shorter than the wild-type Lpp, a helical structural lipoprotein that controls the width of the periplasmic space through its length. The mutant remains susceptible to synergistic growth inhibition by combination of furazolidone and vancomycin; with furazolidone decreasing the vancomycin MIC by eightfold. These findings have clinical relevance, given that the vancomycin concentration required to select the lpp mutation is reachable during typical vancomycin oral administration for treating Clostridioides difficile infections. Combination therapy with furazolidone, however, is likely to prevent emergence and outgrowth of the lpp-mutated Gram-negative coliforms, avoiding exacerbation of the patient's condition during the treatment.

Strain-level diversity of symbiont communities between individuals and populations of a bioluminescent fish

The bioluminescent symbiosis involving the urchin cardinalfish, Siphamia tubifer, and Photobacterium mandapamensis, a luminous member of the Vibrionaceae, is highly specific compared to other bioluminescent fish-bacteria associations. Despite this high degree of specificity, patterns of genetic diversity have been observed for the symbionts from hosts sampled over relatively small spatial scales. We characterized and compared sub-species, strain-level symbiont diversity within and between S. tubifer hosts sampled from the Philippines and Japan using PCR fingerprinting. We then carried out whole genome sequencing of the unique symbiont genotypes identified to characterize the genetic diversity of the symbiont community and the symbiont pangenome. We determined that an individual light organ contains six symbiont genotypes on average, but varied between 1-13. Additionally, we found that there were few genotypes shared between hosts from the same location. A phylogenetic analysis of the unique symbiont strains indicated location-specific clades, suggesting some genetic differentiation in the symbionts between host populations. We also identified symbiont genes that were variable between strains, including luxF, a member of the lux operon, which is responsible for light production. We quantified the light emission and growth rate of two strains missing luxF along with the other strains isolated from the same light organs and determined that strains lacking luxF were dimmer but grew faster than most of the other strains, suggesting a potential metabolic trade-off. This study highlights the importance of strain-level diversity in microbial associations and provides new insight into the underlying genetic architecture of intraspecific symbiont communities within a host.

TRS: a method for determining transcript termini from RNAtag-seq sequencing data

In bacteria, determination of the 3' termini of transcripts plays an essential role in regulation of gene expression, affecting the functionality and stability of the transcript. Several experimental approaches were developed to identify the 3' termini of transcripts, however, these were applied only to a limited number of bacteria and growth conditions. Here we present a straightforward approach to identify 3' termini from widely available RNA-seq data without the need for additional experiments. Our approach relies on the observation that the RNAtag-seq sequencing protocol results in overabundance of reads mapped to transcript 3' termini. We present TRS (Termini by Read Starts), a computational pipeline exploiting this property to identify 3' termini in RNAtag-seq data, and show that the identified 3' termini are highly reliable. Since RNAtag-seq data are widely available for many bacteria and growth conditions, our approach paves the way for studying bacterial transcription termination in an unprecedented scope.

Diurnal cycles drive rhythmic physiology and promote survival in facultative phototrophic bacteria

Bacteria have evolved many strategies to spare energy when nutrients become scarce. One widespread such strategy is facultative phototrophy, which helps heterotrophs supplement their energy supply using light. Our knowledge of the impact that such behaviors have on bacterial fitness and physiology is, however, still limited. Here, we study how a representative of the genus Porphyrobacter, in which aerobic anoxygenic phototrophy is ancestral, responds to different light regimes under nutrient limitation. We show that bacterial survival in stationary phase relies on functional reaction centers and varies depending on the light regime. Under dark-light alternance, our bacterial model presents a diphasic life history dependent on phototrophy: during dark phases, the cells inhibit DNA replication and part of the population lyses and releases nutrients, while subsequent light phases allow for the recovery and renewed growth of the surviving cells. We correlate these cyclic variations with a pervasive pattern of rhythmic transcription which reflects global changes in diurnal metabolic activity. Finally, we demonstrate that, compared to either a phototrophy mutant or a bacteriochlorophyll a overproducer, the wild type strain is better adapted to natural environments, where regular dark-light cycles are interspersed with additional accidental dark episodes. Overall, our results highlight the importance of light-induced biological rhythms in a new model of aerobic anoxygenic phototroph representative of an ecologically important group of environmental bacteria.

Global pathogenomic analysis identifies known and candidate genetic antimicrobial resistance determinants in twelve species

Surveillance programs for managing antimicrobial resistance (AMR) have yielded thousands of genomes suited for data-driven mechanism discovery. We present a workflow integrating pangenomics, gene annotation, and machine learning to identify AMR genes at scale. When applied to 12 species, 27,155 genomes, and 69 drugs, we 1) find AMR gene transfer mostly confined within related species, with 925 genes in multiple species but just eight in multiple phylogenetic classes, 2) demonstrate that discovery-oriented support vector machines outperform contemporary methods at recovering known AMR genes, recovering 263 genes compared to 145 by Pyseer, and 3) identify 142 AMR gene candidates. Validation of two candidates in E. coli BW25113 reveals cases of conditional resistance: ΔcycA confers ciprofloxacin resistance in minimal media with D-serine, and frdD V111D confers ampicillin resistance in the presence of ampC by modifying the overlapping promoter. We expect this approach to be adaptable to other species and phenotypes.

Distributed genotyping and clustering of Neisseria strains reveal continual emergence of epidemic meningococcus over a century

Core genome multilocus sequence typing (cgMLST) is commonly used to classify bacterial strains into different types, for taxonomical and epidemiological applications. However, cgMLST schemes require central databases for the nomenclature of new alleles and sequence types, which must be synchronized worldwide and involve increasingly intensive calculation and storage demands. Here, we describe a distributed cgMLST (dcgMLST) scheme that does not require a central database of allelic sequences and apply it to study evolutionary patterns of epidemic and endemic strains of the genus Neisseria. We classify 69,994 worldwide Neisseria strains into multi-level clusters that assign species, lineages, and local disease outbreaks. We divide Neisseria meningitidis into 168 endemic lineages and three epidemic lineages responsible for at least 9 epidemics in the past century. According to our analyses, the epidemic and endemic lineages experienced very different population dynamics in the past 100 years. Epidemic lineages repetitively emerged from endemic lineages, disseminated worldwide, and apparently disappeared rapidly afterward. We propose a stepwise model for the evolutionary trajectory of epidemic lineages in Neisseria, and expect that the development of similar dcgMLST schemes will facilitate epidemiological studies of other bacterial pathogens.

Depletion of tRNA CCA-adding enzyme in Mycobacterium tuberculosis leads to polyadenylation of transcripts and precursor tRNAs

In reference to gene annotation, more than half of the tRNA species synthesized by Mycobacterium tuberculosis require the enzymatic addition of the cytosine-cytosine-adenine (CCA) tail, which is indispensable for amino acid charging and tRNA functionality. It makes the mycobacterial CCA-adding enzyme essential for survival of the bacterium and a potential target for novel pipelines in drug discovery avenues. Here, we described the rv3907c gene product, originally annotated as poly(A)polymerase (rv3907c, PcnA) as a functional CCA-adding enzyme (CCAMtb) essential for viability of M. tuberculosis. The depletion of the enzyme affected tRNAs maturation, inhibited bacilli growth, and resulted in abundant accumulation of polyadenylated RNAs. We determined the enzymatic activities displayed by the mycobacterial CCAMtb in vitro and studied the effects of inhibiting of its transcription in bacterial cells. We are the first to properly confirm the existence of RNA polyadenylation in mycobacteria, a previously controversial phenomenon, which we found promoted upon CCA-adding enzyme downexpression.

The genomic epidemiology of shigellosis in South Africa

Shigellosis, a leading cause of diarrhoeal mortality and morbidity globally, predominantly affects children under five years of age living in low- and middle-income countries. While whole genome sequence analysis (WGSA) has been effectively used to further our understanding of shigellosis epidemiology, antimicrobial resistance, and transmission, it has been under-utilised in sub-Saharan Africa. In this study, we applied WGSA to large sub-sample of surveillance isolates from South Africa, collected from 2011 to 2015, focussing on Shigella flexneri 2a and Shigella sonnei. We find each serotype is epidemiologically distinct. The four identified S. flexneri 2a clusters having distinct geographical distributions, and antimicrobial resistance (AMR) and virulence profiles, while the four sub-Clades of S. sonnei varied in virulence plasmid retention. Our results support serotype specific lifestyles as a driver for epidemiological differences, show AMR is not required for epidemiological success in S. flexneri, and that the HIV epidemic may have promoted Shigella population expansion.

Regulation of the physiology and virulence of Ralstonia solanacearum by the second messenger 2',3'-cyclic guanosine monophosphate

Previous studies have demonstrated that bis-(3',5')-cyclic diguanosine monophosphate (bis-3',5'-c-di-GMP) is a ubiquitous second messenger employed by bacteria. Here, we report that 2',3'-cyclic guanosine monophosphate (2',3'-cGMP) controls the important biological functions, quorum sensing (QS) signaling systems and virulence in Ralstonia solanacearum through the transcriptional regulator RSp0980. This signal specifically binds to RSp0980 with high affinity and thus abolishes the interaction between RSp0980 and the promoters of target genes. In-frame deletion of RSp0334, which contains an evolved GGDEF domain with a LLARLGGDQF motif required to catalyze 2',3'-cGMP to (2',5')(3',5')-cyclic diguanosine monophosphate (2',3'-c-di-GMP), altered the abovementioned important phenotypes through increasing the intracellular 2',3'-cGMP levels. Furthermore, we found that 2',3'-cGMP, its receptor and the evolved GGDEF domain with a LLARLGGDEF motif also exist in the human pathogen Salmonella typhimurium. Together, our work provides insights into the unusual function of the GGDEF domain of RSp0334 and the special regulatory mechanism of 2',3'-cGMP signal in bacteria.

High-resolution temporal profiling of E. coli transcriptional response

Understanding how cells dynamically adapt to their environment is a primary focus of biology research. Temporal information about cellular behavior is often limited by both small numbers of data time-points and the methods used to analyze this data. Here, we apply unsupervised machine learning to a data set containing the activity of 1805 native promoters in E. coli measured every 10 minutes in a high-throughput microfluidic device via fluorescence time-lapse microscopy. Specifically, this data set reveals E. coli transcriptome dynamics when exposed to different heavy metal ions. We use a bioinformatics pipeline based on Independent Component Analysis (ICA) to generate insights and hypotheses from this data. We discovered three primary, time-dependent stages of promoter activation to heavy metal stress (fast, intermediate, and steady). Furthermore, we uncovered a global strategy E. coli uses to reallocate resources from stress-related promoters to growth-related promoters following exposure to heavy metal stress.

A smooth tubercle bacillus from Ethiopia phylogenetically close to the Mycobacterium tuberculosis complex

The Mycobacterium tuberculosis complex (MTBC) includes several human- and animal-adapted pathogens. It is thought to have originated in East Africa from a recombinogenic Mycobacterium canettii-like ancestral pool. Here, we describe the discovery of a clinical tuberculosis strain isolated in Ethiopia that shares archetypal phenotypic and genomic features of M. canettii strains, but represents a phylogenetic branch much closer to the MTBC clade than to the M. canettii strains. Analysis of genomic traces of horizontal gene transfer in this isolate and previously identified M. canettii strains indicates a persistent albeit decreased recombinogenic lifestyle near the emergence of the MTBC. Our findings support that the MTBC emergence from its putative free-living M. canettii-like progenitor is evolutionarily very recent, and suggest the existence of a continuum of further extant derivatives from ancestral stages, close to the root of the MTBC, along the Great Rift Valley.

The metabolic, virulence and antimicrobial resistance profiles of colonising Streptococcus pneumoniae shift after PCV13 introduction in urban Malawi

Streptococcus pneumoniae causes substantial mortality among children under 5-years-old worldwide. Polysaccharide conjugate vaccines (PCVs) are highly effective at reducing vaccine serotype disease, but emergence of non-vaccine serotypes and persistent nasopharyngeal carriage threaten this success. We investigated the hypothesis that following vaccine, adapted pneumococcal genotypes emerge with the potential for vaccine escape. We genome sequenced 2804 penumococcal isolates, collected 4-8 years after introduction of PCV13 in Blantyre, Malawi. We developed a pipeline to cluster the pneumococcal population based on metabolic core genes into "Metabolic genotypes" (MTs). We show that S. pneumoniae population genetics are characterised by emergence of MTs with distinct virulence and antimicrobial resistance (AMR) profiles. Preliminary in vitro and murine experiments revealed that representative isolates from emerging MTs differed in growth, haemolytic, epithelial infection, and murine colonisation characteristics. Our results suggest that in the context of PCV13 introduction, pneumococcal population dynamics had shifted, a phenomenon that could further undermine vaccine control and promote spread of AMR.