IS1096-mediated DNA rearrangements play a key role in genome evolution of Mycobacterium smegmatis

Insertion sequence transposition activates antimycobacteriophage immunity through an lsr2-silenced lipid metabolism gene island

Insertion sequences (ISs) exist widely in bacterial genomes, but their roles in the evolution of bacterial antiphage defense remain to be clarified. Here, we report that, under the pressure of phage infection, the IS1096 transposition of Mycobacterium smegmatis into the lsr2 gene can occur at high frequencies, which endows the mutant mycobacterium with a broad-spectrum antiphage ability. Lsr2 functions as a negative regulator and directly silences expression of a gene island composed of 11 lipid metabolism-related genes. The complete or partial loss of the gene island leads to a significant decrease of bacteriophage adsorption to the mycobacterium, thus defending against phage infection. Strikingly, a phage that has evolved mutations in two tail-filament genes can re-escape from the lsr2 inactivation-triggered host defense. This study uncovered a new signaling pathway for activating antimycobacteriophage immunity by IS transposition and provided insight into the natural evolution of bacterial antiphage defense.

Assembling the perfect bacterial genome using Oxford Nanopore and Illumina sequencing

A perfect bacterial genome assembly is one where the assembled sequence is an exact match for the organism's genome-each replicon sequence is complete and contains no errors. While this has been difficult to achieve in the past, improvements in long-read sequencing, assemblers, and polishers have brought perfect assemblies within reach. Here, we describe our recommended approach for assembling a bacterial genome to perfection using a combination of Oxford Nanopore Technologies long reads and Illumina short reads: Trycycler long-read assembly, Medaka long-read polishing, Polypolish short-read polishing, followed by other short-read polishing tools and manual curation. We also discuss potential pitfalls one might encounter when assembling challenging genomes, and we provide an online tutorial with sample data (github.com/rrwick/perfect-bacterial-genome-tutorial).

Multiple genetic paths including massive gene amplification allow Mycobacterium tuberculosis to overcome loss of ESX-3 secretion system substrates

The Mycobacterium tuberculosis (Mtb) ESX-3 type VII secretion system plays a critical role in iron acquisition. Infection of mice with highly attenuated Mtb deletion mutants lacking esxG or esxH, genes encoding key ESX-3 substrates, unexpectedly yielded suppressor mutants with restored capacity to grow in vivo and in vitro in the absence of iron supplementation. Whole-genome sequencing identified two mechanisms of suppression, the disruption of a transcriptional repressor that regulates expression of an ESX-3 paralogous region encoding EsxR and EsxS, and a massive 38- to 60-fold gene amplification of this same region. These data are significant because they reveal a previously unrecognized iron acquisition regulon and inform mechanisms of Mtb chromosome evolution.

Mycobacterium tuberculosis (Mtb) possesses five type VII secretion systems (T7SS), virulence determinants that include the secretion apparatus and associated secretion substrates. Mtb strains deleted for the genes encoding substrates of the ESX-3 T7SS, esxG or esxH, require iron supplementation for in vitro growth and are highly attenuated in vivo. In a subset of infected mice, suppressor mutants of esxG or esxH deletions were isolated, which enabled growth to high titers or restored virulence. Suppression was conferred by mechanisms that cause overexpression of an ESX-3 paralogous region that lacks genes for the secretion apparatus but encodes EsxR and EsxS, apparent ESX-3 orphan substrates that functionally compensate for the lack of EsxG or EsxH. The mechanisms include the disruption of a transcriptional repressor and a massive 38- to 60-fold gene amplification. These data identify an iron acquisition regulon, provide insight into T7SS, and reveal a mechanism of Mtb chromosome evolution involving “accordion-type” amplification.

Elimination of PknL and MSMEG_4242 in Mycobacterium smegmatis alters the character of the outer cell envelope and selects for mutations in Lsr2

•

Elimination of pknL and adjacent gene MSMEG_4242 in M. smegmatis produces rough & smooth colonies.

•

All smooth colony mutants have inactivated lsr2 genes.

•

Inactivated lsr2 leads to high expression of MSMEG_4727 and increased synthesis of LOS.

•

Smooth mutants don’t form biofilms, have longer bacilli and increased sliding motility.

•

Rough strains deleted for pknL and MSMEG_4242 form biofilms but have aberrant sliding motility.

•

We propose a phosphorylation cascade of PknB phosphorylating PknL that then phosphorylates Lsr2.

Four serine/threonine kinases are present in all mycobacteria: PknA, PknB, PknG and PknL. PknA and PknB are essential for growth and replication, PknG regulates metabolism, but little is known about PknL. Inactivation of pknL and adjacent regulator MSMEG_4242 in rough colony M. smegmatis mc²155 produced both smooth and rough colonies. Upon restreaking rough colonies, smooth colonies appeared at a frequency of ~ 1/250. Smooth mutants did not form biofilms, showed increased sliding motility and anomalous lipids on thin-layer chromatography, identified by mass spectrometry as lipooligosaccharides and perhaps also glycopeptidolipids. RNA-seq and Sanger sequencing revealed that all smooth mutants had inactivated lsr2 genes due to mutations and different IS1096 insertions. When complemented with lsr2, the colonies became rough, anomalous lipids disappeared and sliding motility decreased. Smooth mutants showed increased expression of IS1096 transposase TnpA and MSMEG_4727, which encodes a protein similar to PKS5. When MSMEG_4727 was deleted, smooth pknL/MSMEG_4242/lsr2 mutants reverted to rough, formed good biofilms, their motility decreased slightly and their anomalous lipids disappeared. Rough delpknL/del4242 mutants formed poor biofilms and showed decreased, aberrant sliding motility and both phenotypes were complemented with the two deleted genes. Inactivation of lsr2 changes colony morphology from rough to smooth, augments sliding motility and increases expression of MSMEG_4727 and other enzymes synthesizing lipooligosaccharides, apparently preventing biofilm formation. Similar morphological phase changes occur in other mycobacteria, likely reflecting environmental adaptations. PknL and MSMEG_4242 regulate lipid components of the outer cell envelope and their absence selects for lsr2 inactivation. A regulatory, phosphorylation cascade model is proposed.

A cytoplasmic peptidoglycan amidase homologue controls mycobacterial cell wall synthesis

Regulation of cell wall assembly is essential for bacterial survival and contributes to pathogenesis and antibiotic tolerance in Mycobacterium tuberculosis (Mtb). However, little is known about how the cell wall is regulated in stress. We found that CwlM, a protein homologous to peptidoglycan amidases, coordinates peptidoglycan synthesis with nutrient availability. Surprisingly, CwlM is sequestered from peptidoglycan (PG) by localization in the cytoplasm, and its enzymatic function is not essential. Rather, CwlM is phosphorylated and associates with MurA, the first enzyme in PG precursor synthesis. Phosphorylated CwlM activates MurA ~30 fold. CwlM is dephosphorylated in starvation, resulting in lower MurA activity, decreased cell wall metabolism, and increased tolerance to multiple antibiotics. A phylogenetic analysis of cwlM implies that localization in the cytoplasm drove the evolution of this factor. We describe a system that controls cell wall metabolism in response to starvation, and show that this regulation contributes to antibiotic tolerance.

Origins of a 350-kilobase genomic duplication in Mycobacterium tuberculosis and its impact on virulence

In the present study, we have investigated the evolution and impact on virulence of a 350-kb genomic duplication present in the most recently evolved members of the Mycobacterium tuberculosis East Asian lineage. In a mouse model of infection, comparing HN878 subclones HN878-27 (no duplication) and HN878-45 (with the 350-kb duplication) revealed that the latter is impaired for in vivo growth during the initial 3 weeks of infection. Furthermore, the median survival time of mice infected with isolate HN878-45 is significantly longer (77 days) than that of mice infected with HN878-27. Whole-genome sequencing of both isolates failed to reveal any mutational events other than the duplication that could account for such a substantial difference in virulence. Although we and others had previously speculated that the 350-kb duplication arose in response to some form of host-applied selective pressure (P. Domenech, G. S. Kolly, L. Leon-Solis, A. Fallow, M. B. Reed, J. Bacteriol. 192: 4562-4570, 2010, and B. Weiner, J. Gomez, T. C. Victor, R. M. Warren, A. Sloutsky, B. B. Plikaytis, J. E. Posey, P. D. van Helden, N. C. Gey van Pittius, M. Koehrsen, P. Sisk, C. Stolte, J. White, S. Gagneux, B. Birren, D. Hung, M. Murray, J. Galagan, PLoS One 7: e26038, 2012), here we show that this large chromosomal amplification event is very rapidly selected within standard in vitro broth cultures in a range of isolates. Indeed, subclones harboring the duplication were detectable after just five rounds of in vitro passage. In contrast, the duplication appears to be highly unstable in vivo and is negatively selected during the later stages of infection in mice. We believe that the rapid in vitro evolution of M. tuberculosis is an underappreciated aspect of its biology that is often ignored, despite the fact that it has the potential to confound the data and conclusions arising from comparative studies of isolates at both the genotypic and phenotypic levels.

Large-scale genomic analysis suggests a neutral punctuated dynamics of transposable elements in bacterial genomes

Insertion sequences (IS) are the simplest and most abundant form of transposable DNA found in bacterial genomes. When present in multiple copies, it is thought that they can promote genomic plasticity and genetic exchange, thus being a major force of evolutionary change. The main processes that determine IS content in genomes are, though, a matter of debate. In this work, we take advantage of the large amount of genomic data currently available and study the abundance distributions of 33 IS families in 1811 bacterial chromosomes. This allows us to test simple models of IS dynamics and estimate their key parameters by means of a maximum likelihood approach. We evaluate the roles played by duplication, lateral gene transfer, deletion and purifying selection. We find that the observed IS abundances are compatible with a neutral scenario where IS proliferation is controlled by deletions instead of purifying selection. Even if there may be some cases driven by selection, neutral behavior dominates over large evolutionary scales. According to this view, IS and hosts tend to coexist in a dynamic equilibrium state for most of the time. Our approach also allows for a detection of recent IS expansions, and supports the hypothesis that rapid expansions constitute transient events—punctuations—during which the state of coexistence of IS and host becomes perturbated.

Insertion sequences (IS) are mobile genetic elements found in most prokaryotic genomes. They are able to autonomously change position and proliferate in chromosomes. The nature of the coevolutionary dynamics of IS with the genome that hosts them is a matter of debate: Do IS proliferate to the point of causing the extinction of the host? Is it possible that IS and hosts stably coexist? Can environmental perturbations cause IS expansions? What is the role of selection in controlling IS copy number? In this study, we have analysed abundance patterns of IS families to test two different evolutionary hypotheses: in the first one IS evolve neutrally, while in the second case they are affected by selection. Our results indicate that, most of the time, IS and their hosts coexist stably in a neutral scenario where the proliferation of IS through duplications and lateral gene transfer is balanced by regular deletions. Occasionally, though, this balance may be disrupted, causing temporary explosions of IS abundance.

Population level analysis of evolved mutations underlying improvements in plant hemicellulose and cellulose fermentation by Clostridium phytofermentans

Background

The complexity of plant cell walls creates many challenges for microbial decomposition. Clostridium phytofermentans, an anaerobic bacterium isolated from forest soil, directly breaks down and utilizes many plant cell wall carbohydrates. The objective of this research is to understand constraints on rates of plant decomposition by Clostridium phytofermentans and identify molecular mechanisms that may overcome these limitations.

Results

Experimental evolution via repeated serial transfers during exponential growth was used to select for C. phytofermentans genotypes that grow more rapidly on cellobiose, cellulose and xylan. To identify the underlying mutations an average of 13,600,000 paired-end reads were generated per population resulting in ∼300 fold coverage of each site in the genome. Mutations with allele frequencies of 5% or greater could be identified with statistical confidence. Many mutations are in carbohydrate-related genes including the promoter regions of glycoside hydrolases and amino acid substitutions in ABC transport proteins involved in carbohydrate uptake, signal transduction sensors that detect specific carbohydrates, proteins that affect the export of extracellular enzymes, and regulators of unknown specificity. Structural modeling of the ABC transporter complex proteins suggests that mutations in these genes may alter the recognition of carbohydrates by substrate-binding proteins and communication between the intercellular face of the transmembrane and the ATPase binding proteins.

Conclusions

Experimental evolution was effective in identifying molecular constraints on the rate of hemicellulose and cellulose fermentation and selected for putative gain of function mutations that do not typically appear in traditional molecular genetic screens. The results reveal new strategies for evolving and engineering microorganisms for faster growth on plant carbohydrates.

Unusual large-scale chromosomal rearrangements in Mycobacterium tuberculosis Beijing B0/W148 cluster isolates

The Mycobacterium tuberculosis (MTB) Beijing family isolates are geographically widespread, and there are examples of Beijing isolates that are hypervirulent and associated with drug resistance. One-fourth of Beijing genotype isolates found in Russia belong to the B0/W148 group. The aim of the present study was to investigate features of these endemic strains on a genomic level. Four Russian clinical isolates of this group were sequenced, and the data obtained was compared with published sequences of various MTB strain genomes, including genome of strain W-148 of the same B0/W148 group. The comparison of the W-148 and H37Rv genomes revealed two independent inversions of large segments of the chromosome. The same inversions were found in one of the studied strains after deep sequencing using both the fragment and mate-paired libraries. Additionally, inversions were confirmed by RFLP hybridization analysis. The discovered rearrangements were verified by PCR in all four newly sequenced strains in the study and in four additional strains of the same Beijing B0/W148 group. The other 32 MTB strains from different phylogenetic lineages were tested and revealed no inversions. We suggest that the initial largest inversion changed the orientation of the three megabase (Mb) segment of the chromosome, and the second one occurred in the previously inverted region and partly restored the orientation of the 2.1 Mb inner segment of the region. This is another remarkable example of genomic rearrangements in the MTB in addition to the recently published of large-scale duplications. The described cases suggest that large-scale genomic rearrangements in the currently circulating MTB isolates may occur more frequently than previously considered, and we hope that further studies will help to determine the exact mechanism of such events.

Direct and inverted repeats elicit genetic instability by both exploiting and eluding DNA double-strand break repair systems in mycobacteria

Repetitive DNA sequences with the potential to form alternative DNA conformations, such as slipped structures and cruciforms, can induce genetic instability by promoting replication errors and by serving as a substrate for DNA repair proteins, which may lead to DNA double-strand breaks (DSBs). However, the contribution of each of the DSB repair pathways, homologous recombination (HR), non-homologous end-joining (NHEJ) and single-strand annealing (SSA), to this sort of genetic instability is not fully understood. Herein, we assessed the genome-wide distribution of repetitive DNA sequences in the Mycobacterium smegmatis, Mycobacterium tuberculosis and Escherichia coli genomes, and determined the types and frequencies of genetic instability induced by direct and inverted repeats, both in the presence and in the absence of HR, NHEJ, and SSA. All three genomes are strongly enriched in direct repeats and modestly enriched in inverted repeats. When using chromosomally integrated constructs in M. smegmatis, direct repeats induced the perfect deletion of their intervening sequences ~1,000-fold above background. Absence of HR further enhanced these perfect deletions, whereas absence of NHEJ or SSA had no influence, suggesting compromised replication fidelity. In contrast, inverted repeats induced perfect deletions only in the absence of SSA. Both direct and inverted repeats stimulated excision of the constructs from the attB integration sites independently of HR, NHEJ, or SSA. With episomal constructs, direct and inverted repeats triggered DNA instability by activating nucleolytic activity, and absence of the DSB repair pathways (in the order NHEJ>HR>SSA) exacerbated this instability. Thus, direct and inverted repeats may elicit genetic instability in mycobacteria by 1) directly interfering with replication fidelity, 2) stimulating the three main DSB repair pathways, and 3) enticing L5 site-specific recombination.

Polycyclic aromatic hydrocarbon degrading gene islands in five pyrene-degrading Mycobacterium isolates from different geographic locations

Mycobacterium sp. strain KMS utilizes pyrene, a high-molecular weight polycyclic aromatic hydrocarbon (PAH), as a sole carbon source. Bioinformatic analysis of the genome of isolate KMS predicted 25 genes with the potential to encode 17 pyrene-induced proteins identified by proteomics; these genes were clustered on both the chromosome and a circular plasmid. RT-PCR analysis of total RNA isolated from KMS cells grown with or without pyrene showed that the presence of pyrene increased the transcript accumulation of 20 of the predicted chromosome- and plasmid-located genes encoding pyrene-induced proteins. The transcribed genes from both the chromosome and a circular plasmid were within larger regions containing genes required for PAH degradation constituting PAH-degrading gene islands. Genes encoding integrases and transposases were found within and outside the PAH-degrading gene islands. The lower GC content of the genes within the gene island (61%-64%) compared with the average genome content (68%) suggested that these mycobacteria initially acquired these genes by horizontal gene transfer. Synteny was detected for the PAH-degrading islands in the genomes of two additional Mycobacterium isolates from the same PAH-polluted site and of two other pyrene-degrading Mycobacterium from different sites in the United States of America. Consequently, the gene islands have been conserved from a common ancestral strain.

Analysis of insertion sequences in thermophilic cyanobacteria: exploring the mechanisms of establishing, maintaining, and withstanding high insertion sequence abundance

Insertion sequences (ISs) are simple mobile genetic elements capable of relocating within a genome. Through this transposition activity, they are known to create mutations which are mostly deleterious to the cell, although occasionally they are beneficial. Two closely related isolates of thermophilic Synechococcus species from hot spring microbial mats are known to harbor a large number of diverse ISs. To explore the mechanism of IS acquisition within natural populations and survival in the face of high IS abundance, we examined IS content and location in natural populations of Synechococcus by comparing metagenomic data to the genomes of fully sequenced cultured isolates. The observed IS distribution in the metagenome was equivalent to the distribution in the isolates, indicating that the cultured isolates are appropriate models for the environmental population. High sequence conservation between IS families shared between the two isolates suggests that ISs are able to move between individuals within populations and between species via lateral gene transfer, consistent with models for IS family accumulation. Most IS families show evidence of recent activity, and interruption of critical genes in some individuals was observed, demonstrating that transposition is an ongoing mutational force in the populations.

Massive gene duplication event among clinical isolates of the Mycobacterium tuberculosis W/Beijing family

As part of our effort to uncover the molecular basis for the phenotypic variation among clinical Mycobacterium tuberculosis isolates, we have previously reported that isolates belonging to the W/Beijing lineage constitutively overexpress the DosR-regulated transcriptional program. While generating dosR knockouts in two independent W/Beijing sublineages, we were surprised to discover that they possess two copies of dosR. This dosR amplification is part of a massive genomic duplication spanning 350 kb and encompassing >300 genes. In total, this equates to 8% of the genome being present as two copies. The presence of IS6110 elements at both ends of the region of duplication, and in the novel junction region, suggests that it arose through unequal homologous recombination of sister chromatids at the IS6110 sequences. Analysis of isolates representing the major M. tuberculosis lineages has revealed that the 350-kb duplication is restricted to the most recently evolved sublineages of the W/Beijing family. Within these isolates, the duplication is partly responsible for the constitutive dosR overexpression phenotype. Although the nature of the selection event giving rise to the duplication remains unresolved, its evolution is almost certainly the result of specific selective pressure(s) encountered inside the host. A preliminary in vitro screen has failed to reveal a role of the duplication in conferring resistance to common antitubercular drugs, a trait frequently associated with W/Beijing isolates. Nevertheless, this first description of a genetic remodeling event of this nature for M. tuberculosis further highlights the potential for the evolution of diversity in this important global pathogen.

Molecular characterization of heterologous HIV-1gp120 gene expression disruption in mycobacterium bovis BCG host strain: a critical issue for engineering mycobacterial based-vaccine vectors

Mycobacterium bovis Bacillus Calmette-Guérin (BCG) as a live vector of recombinant bacterial vaccine is a promising system to be used. In this study, we evaluate the disrupted expression of heterologous HIV-1gp120 gene in BCG Pasteur host strain using replicative vectors pMV261 and pJH222. pJH222 carries a lysine complementing gene in BCG lysine auxotrophs. The HIV-1 gp120 gene expression was regulated by BCG hsp60 promoter (in plasmid pMV261) and Mycobacteria spp. α-antigen promoter (in plasmid pJH222). Among 14 rBCG:HIV-1gp120 (pMV261) colonies screened, 12 showed a partial deletion and two showed a complete deletion. However, deletion was not observed in all 10 rBCG:HIV-1gp120 (pJH222) colonies screened. In this study, we demonstrated that E. coli/Mycobacterial expression vectors bearing a weak promoter and lysine complementing gene in a recombinant lysine auxotroph of BCG could prevent genetic rearrangements and disruption of HIV 1gp120 gene expression, a key issue for engineering Mycobacterial based vaccine vectors.