Molecular mechanisms of multiple toxin-antitoxin systems are coordinated to govern the persister phenotype

HigA2 (Rv2021c) Is a Transcriptional Regulator with Multiple Regulatory Targets in Mycobacterium tuberculosis

Toxin-antitoxin (TA) systems are the major mechanism for persister formation in Mycobacterium tuberculosis (Mtb). Previous studies found that HigBA2 (Rv2022c-Rv2021c), a predicted type II TA system of Mtb, could be activated for transcription in response to multiple stresses such as anti-tuberculosis drugs, nutrient starvation, endure hypoxia, acidic pH, etc. In this study, we determined the binding site of HigA2 (Rv2021c), which is located in the coding region of the upstream gene higB2 (Rv2022c), and the conserved recognition motif of HigA2 was characterized via oligonucleotide mutation. Eight binding sites of HigA2 were further found in the Mtb genome according to the conserved motif. RT-PCR showed that HigA2 can regulate the transcription level of all eight of these genes and three adjacent downstream genes. DNA pull-down experiments showed that twelve functional regulators sense external regulatory signals and may regulate the transcription of the HigBA2 system. Of these, Rv0903c, Rv0744c, Rv0474, Rv3124, Rv2603c, and Rv3583c may be involved in the regulation of external stress signals. In general, we identified the downstream target genes and possible upstream regulatory genes of HigA2, which paved the way for the illustration of the persistence establishment mechanism in Mtb.

Leveraging Cancer Phenotypic Plasticity for Novel Treatment Strategies

Cancer cells, like all other organisms, are adept at switching their phenotype to adjust to the changes in their environment. Thus, phenotypic plasticity is a quantitative trait that confers a fitness advantage to the cancer cell by altering its phenotype to suit environmental circumstances. Until recently, new traits, especially in cancer, were thought to arise due to genetic factors; however, it is now amply evident that such traits could also emerge non-genetically due to phenotypic plasticity. Furthermore, phenotypic plasticity of cancer cells contributes to phenotypic heterogeneity in the population, which is a major impediment in treating the disease. Finally, plasticity also impacts the group behavior of cancer cells, since competition and cooperation among multiple clonal groups within the population and the interactions they have with the tumor microenvironment also contribute to the evolution of drug resistance. Thus, understanding the mechanisms that cancer cells exploit to tailor their phenotypes at a systems level can aid the development of novel cancer therapeutics and treatment strategies. Here, we present our perspective on a team medicine-based approach to gain a deeper understanding of the phenomenon to develop new therapeutic strategies.

Preventing toxicity in toxin-antitoxin systems: An overview of regulatory mechanisms

Toxin-antitoxin systems (TAs) are generally two-component genetic modules present in almost every prokaryotic genome. The production of the free and active toxin is able to disrupt key cellular processes leading to the growth inhibition or death of its host organism in absence of its cognate antitoxin. The functions attributed to TAs rely on this lethal phenotype ranging from mobile genetic elements stabilization to phage defense. Their abundance in prokaryotic genomes as well as their lethal potential make them attractive targets for new antibacterial strategies. The hijacking of TAs requires a deep understanding of their regulation to be able to design such approach. In this review, we summarize the accumulated knowledge on how bacteria cope with these toxic genes in their genome. The characterized TAs can be grouped based on the way they prevent toxicity. Some systems rely on a tight control of the expression to prevent the production of the toxin while others control the activity of the toxin at the post-translational level.

Genes and genome-resolved metagenomics reveal the microbial functional make up of Amazon peatlands under geochemical gradients

The Pastaza-Marañón Foreland Basin (PMFB) holds the most extensive tropical peatland area in South America. PMFB peatlands store ~7.07 Gt of organic carbon interacting with multiple microbial heterotrophic, methanogenic, and other aerobic/anaerobic respirations. Little is understood about the contribution of distinct microbial community members inhabiting tropical peatlands. Here, we studied the metagenomes of three geochemically distinct peatlands spanning minerotrophic, mixed, and ombrotrophic conditions. Using gene- and genome-centric approaches, we evaluate the functional potential of the underlying microbial communities. Abundance analyses show significant differences in C, N, P, and S acquisition genes. Furthermore, community interactions mediated by toxin-antitoxin and CRISPR-Cas systems were enriched in oligotrophic soils, suggesting that non-metabolic interactions may exert additional controls in low-nutrient environments. Additionally, we reconstructed 519 metagenome-assembled genomes spanning 28 phyla. Our analyses detail key differences across the geochemical gradient in the predicted microbial populations involved in degradation of organic matter, and the cycling of N and S. Notably, we observed differences in the nitric oxide (NO) reduction strategies between sites with high and low N2 O fluxes and found phyla putatively capable of both NO and sulfate reduction. Our findings detail how gene abundances and microbial populations are influenced by geochemical differences in tropical peatlands.

Bacterial toxin-antitoxin systems: Novel insights on toxin activation across populations and experimental shortcomings

The alarming rise in hard-to-treat bacterial infections is of great concern to human health. Thus, the identification of molecular mechanisms that enable the survival and growth of pathogens is of utmost urgency for the development of more efficient antimicrobial therapies. In challenging environments, such as presence of antibiotics, or during host infection, metabolic adjustments are essential for microorganism survival and competitiveness. Toxin-antitoxin systems (TASs) consisting of a toxin with metabolic modulating activity and a cognate antitoxin that antagonizes that toxin are important elements in the arsenal of bacterial stress defense. However, the exact physiological function of TA systems is highly debatable and with the exception of stabilization of mobile genetic elements and phage inhibition, other proposed biological functions lack a broad consensus. This review aims at gaining new insights into the physiological effects of TASs in bacteria and exploring the experimental shortcomings that lead to discrepant results in TAS research. Distinct control mechanisms ensure that only subsets of cells within isogenic cultures transiently develop moderate levels of toxin activity. As a result, TASs cause phenotypic growth heterogeneity rather than cell stasis in the entire population. It is this feature that allows bacteria to thrive in diverse environments through the creation of subpopulations with different metabolic rates and stress tolerance programs.

Metabolome shifts triggered by chlorine sanitisation induce Escherichia coli on fresh produce into the viable but nonculturable state

Facing the increasing occurrence of "big six" Escherichia coli outbreaks linked to fresh produce, chlorine-based sanitisers are widely used for fresh produce decontamination in recent years. However, latest finding that chlorine may induce E. coli cells into a viable not nonculturable (VBNC) state is bringing a new challenge to the fresh produce industry. VBNC cells are undetectable by the plate count test, and yet they retain pathogenicity and are more antibiotic-resistant than culturable cells. As a result, their eradication is critical to ensure the safety of fresh produce. Understanding VBNC cells at the metabolic level may provide a breakthrough for their eradication. Therefore, this study was carried out to collect the VBNC pathogenic E. coli (O26:H11, O121:H19, and O157:H7) cells from chlorine-treated pea sprouts and characterise them using NMR-based metabolomics. From the globally increased metabolite contents detected in the VBNC E. coli cells as compared to the culturable cells, mechanisms underlying E. coli's VBNC induction were elucidated. These include rendering the energy generation scheme to become more compatible with the lowered energy needs, disaggregating protein aggregates to release amino acids for osmoprotection and later resuscitation, as well as increasing cAMP content to downregulate RpoS. These identified metabolic characteristics can inspire future development of targeted measures for VBNC E. coli cell inhibition. Our methods can also be applied to other pathogens to help lower the risk of overall foodborne diseases.

Impact of Genomic Deletion RD16 on the Expression of the Mycobacterium bovis BCG Moreau VapBC47 Toxin-Antitoxin System

Mycobacterium bovis BCG is the only vaccine against tuberculosis. The variable forms of cultivation throughout the years, before seed-lots were developed, allowed in vitro evolution of the original strain, generating a family of vaccines with different phenotypic and genotypic characteristics. Molecular studies revealed regions of difference (RDs) in the genomes of the various BCG strains. This work aims to characterize the gene pair rv3407-rv3408 (vapB47-vapC47), coding for a toxin-antitoxin system of the VapBC family, and to evaluate possible transcriptional effects due to the adjacent BCG Moreau-specific genomic deletion RD16. We show that these genes are co-transcribed in BCG strains Moreau and Pasteur, and that the inactivation of an upstream transcriptional repressor (Rv3405c) due to RD16 has a polar effect, leading to increased vapBC47 expression. Furthermore, we detect VapB47 DNA binding in vitro, dependent on a 5' vapB47 sequence that contributes to a palindrome, spanning the promoter and coding region. Our data shed light on the regulation of VapBC systems and on the impact of the BCG Moreau RD16 deletion in the expression of adjacent genes, contributing to a better understanding of BCG Moreau physiology.

Comparative analysis of five type II TA systems identified in Pseudomonas aeruginosa reveals their contributions to persistence and intracellular survival

Background

Pseudomonas aeruginosa is a grave nosocomial pathogen that persistently inhabits the lungs of patients with cystic fibrosis (CF) and causes various chronic infections. The bacterial toxin-antitoxin (TA) system is associated with latent and long-term infections, but the underlying mechanisms remain to be fully characterized.

Methods

We here investigated the diversity and function of five genomic type II TA systems widely distributed among P. aeruginosa clinical isolates. We also examined the distinct structural features of the toxin protein from different TA systems and characterized their contributions to persistence, invasion ability, and intracellular infection caused by P. aeruginosa.

Results

ParDE, PA1030/PA1029, and HigBA could modulate persister cell formation under treatment with specific antibiotics. Furthermore, cell-based transcriptional and invasion assays revealed that PA1030/PA1029 and HigBA TA systems were critical for intracellular survival.

Discussion

Our results highlight the prevalence and diverse roles of type II TA systems in P. aeruginosa and evaluate the possibility of using PA1030/PA1029 and HigBA TA pairs as targets for novel antibiotic treatments.

The High Mutational Sensitivity of ccdA Antitoxin Is Linked to Codon Optimality

Deep mutational scanning studies suggest that synonymous mutations are typically silent and that most exposed, nonactive-site residues are tolerant to mutations. Here, we show that the ccdA antitoxin component of the Escherichia coli ccdAB toxin-antitoxin system is unusually sensitive to mutations when studied in the operonic context. A large fraction (∼80%) of single-codon mutations, including many synonymous mutations in the ccdA gene shows inactive phenotype, but they retain native-like binding affinity towards cognate toxin, CcdB. Therefore, the observed phenotypic effects are largely not due to alterations in protein structure/stability, consistent with a large region of CcdA being intrinsically disordered. E. coli codon preference and strength of ribosome-binding associated with translation of downstream ccdB gene are found to be major contributors of the observed ccdA mutant phenotypes. In select cases, proteomics studies reveal altered ratios of CcdA:CcdB protein levels in vivo, suggesting that the ccdA mutations likely alter relative translation efficiencies of the two genes in the operon. We extend these results by studying single-site synonymous mutations that lead to loss of function phenotypes in the relBE operon upon introduction of rarer codons. Thus, in their operonic context, genes are likely to be more sensitive to both synonymous and nonsynonymous point mutations than inferred previously.

Phenomenological Approach to Cancer Cell Persistence

Drug persistence is a phenomenon by which a small percentage of cancer cells survive the presentation of targeted therapy by transitioning to a quiescent state. Eventually some of these persister cells can transition back to an active growing state and give rise to resistant tumors. Here we introduce a quantitative genetics approach to drug-exposed populations of cancer cells in order to interpret recent experimental data regarding inheritance of persister probability. Our results indicate that alternating periods of drug treatment and drug removal may not be an effective strategy for eliminating persisters.

Integrative biology of persister cell formation: molecular circuitry, phenotypic diversification and fitness effects

Microbial populations often contain persister cells, which reduce the extinction risk upon sudden stresses. Persister cell formation is deeply intertwined with physiology. Due to this complexity, it cannot be satisfactorily understood by focusing only on mechanistic, physiological or evolutionary aspects. In this review, we take an integrative biology perspective to identify common principles of persister cell formation, which might be applicable across evolutionary-distinct microbes. Persister cells probably evolved to cope with a fundamental trade-off between cellular stress and growth tasks, as any biosynthetic resource investment in growth-supporting proteins is at the expense of stress tasks and vice versa. Natural selection probably favours persister cell subpopulation formation over a single-phenotype strategy, where each cell is prepared for growth and stress to a suboptimal extent, since persister cells can withstand harsher environments and their coexistence with growing cells leads to a higher fitness. The formation of coexisting phenotypes requires bistable molecular circuitry. Bistability probably emerges from growth-modulated, positive feedback loops in the cell's growth versus stress control network, involving interactions between sigma factors, guanosine pentaphosphate and toxin-antitoxin (TA) systems. We conclude that persister cell formation is most likely a response to a sudden reduction in growth rate, which can be achieved by antibiotic addition, nutrient starvation, sudden stresses, nutrient transitions or activation of a TA system.

Comparative Analysis of Diverse Acetyltransferase-Type Toxin-Antitoxin Loci in Klebsiella pneumoniae

Toxin-antitoxin (TA) modules containing a Gcn5-related N-acetyltransferase (GNAT) toxin domain regulate bacterial physiology under adverse environmental stresses. Multiple GNAT-ribbon-helix-helix domain (RHH) TA loci have been identified in single bacterial genomes. However, their diversity and interactions are still obscure. Our previous analysis showed that the GNAT toxin of Klebsiella pneumoniae, KacT, introduces antibiotic tolerance and the toxicity of GNAT is neutralized by KacA, an RHH antitoxin. We here present a phylogenetic analysis of GNAT toxins of more than 1,000 GNAT-RHH pairs detected in completely sequenced K. pneumoniae genomes, revealing that the GNAT toxins are diverse and grouped into four distinct clades. Overexpression of GNAT toxins representative of each of the four clades halts the cell growth of K. pneumoniae, while the coexpression of the cognate RHH antitoxin neutralizes GNAT toxicity. We also identify point mutations that inactivate the GNAT toxins. Moreover, we observe a cross-interaction between GNAT-RHH pairs encoded by different replicons, where a chromosomal toxin (KacT2) can be neutralized by its cognate RHH antitoxin (KacA2 on a chromosome) and another antitoxin (KacA3 on a plasmid). Finally, statistical analysis of the distribution of GNAT-RHH loci in K. pneumoniae strains shows pronounced deviation from random distribution within the same clades. Moreover, we also obtain statistically significant correlations between different clades, which we discuss in terms of the experimental results.

IMPORTANCE Elucidating the roles of multifaceted GNAT-RHH TA loci is essential for understanding how these TAs interact among themselves. Recently, the reaction mechanisms and structures of several GNAT-RHH pairs have been reported. While bacterial strains can carry multiple GNAT-RHH loci with diverse origins, studies on the possible cross-interactions of these TA pairs are still limited. Here, we find that 1,000 predicted GNAT toxins of K. pneumoniae can be grouped into four distinct clades. The distributions of TA loci within these clades in K. pneumoniae strains are highly nonrandom, with the presence of a single locus of each clade per strain being highly overrepresented. Moreover, the toxicity of a GNAT toxin encoded by a chromosome was alleviated by a noncognate RHH antitoxin on a plasmid. These results might yield a profound understanding of the widespread GNAT-RHH TA pairs and the cross-interactions between noncognate TA pairs located on different replicons.

Disrupting Mechanisms that Regulate Genomic Repeat Elements to Combat Cancer and Drug Resistance

Despite advancements in understanding cancer pathogenesis and the development of many effective therapeutic agents, resistance to drug treatment remains a widespread challenge that substantially limits curative outcomes. The historical focus on genetic evolution under drug “pressure” as a key driver of resistance has uncovered numerous mechanisms of therapeutic value, especially with respect to acquired resistance. However, recent discoveries have also revealed a potential role for an ancient evolutionary balance between endogenous “viral” elements in the human genome and diverse factors involved in their restriction in tumor evolution and drug resistance. It has long been appreciated that the stability of genomic repeats such as telomeres and centromeres affect tumor fitness, but recent findings suggest that de-regulation of other repetitive genome elements, including retrotransposons, might also be exploited as cancer therapy. This review aims to present an overview of these recent findings.

Overexpression of mqsR in Xylella fastidiosa Leads to a Priming Effect of Cells to Copper Stress Tolerance

Copper-based compounds are widely used in agriculture as a chemical strategy to limit the spread of multiple plant diseases; however, the continuous use of this heavy metal has caused environmental damage as well as the development of copper-resistant strains. Thus, it is important to understand how the bacterial phytopathogens evolve to manage with this metal in the field. The MqsRA Toxin-Antitoxin system has been recently described for its function in biofilm formation and copper tolerance in Xylella fastidiosa, a plant-pathogen bacterium responsible for economic damage in several crops worldwide. Here we identified differentially regulated genes by X. fastidiosa MqsRA by assessing changes in global gene expression with and without copper. Results show that mqsR overexpression led to changes in the pattern of cell aggregation, culminating in a global phenotypic heterogeneity, indicative of persister cell formation. This phenotype was also observed in wild-type cells but only in the presence of copper. This suggests that MqsR regulates genes that alter cell behavior in order to prime them to respond to copper stress, which is supported by RNA-Seq analysis. To increase cellular tolerance, proteolysis and efflux pumps and regulator related to multidrug resistance are induced in the presence of copper, in an MqsR-independent response. In this study we show a network of genes modulated by MqsR that is associated with induction of persistence in X. fastidiosa. Persistence in plant-pathogenic bacteria is an important genetic tolerance mechanism still neglected for management of phytopathogens in agriculture, for which this work expands the current knowledge and opens new perspectives for studies aiming for a more efficient control in the field.

A minimal model for gene expression dynamics of bacterial type II toxin-antitoxin systems

Toxin-antitoxin (TA) modules are part of most bacteria's regulatory machinery for stress responses and general aspects of their physiology. Due to the interplay of a long-lived toxin with a short-lived antitoxin, TA modules have also become systems of interest for mathematical modelling. Here we resort to previous modelling efforts and extract from these a minimal model of type II TA system dynamics on a timescale of hours, which can be used to describe time courses derived from gene expression data of TA pairs. We show that this model provides a good quantitative description of TA dynamics for the 11 TA pairs under investigation here, while simpler models do not. Our study brings together aspects of Biophysics with its focus on mathematical modelling and Computational Systems Biology with its focus on the quantitative interpretation of 'omics' data. This mechanistic model serves as a generic transformation of time course information into kinetic parameters. The resulting parameter vector can, in turn, be mechanistically interpreted. We expect that TA pairs with similar mechanisms are characterized by similar vectors of kinetic parameters, allowing us to hypothesize on the mode of action for TA pairs still under discussion.

The Staphylococcus aureus toxin-antitoxin system YefM-YoeB is associated with antibiotic tolerance and extracellular dependent biofilm formation

The high antibiotic tolerance of Staphylococcus aureus biofilms is associated with challenges for treating periprosthetic joint infection. The toxin–antitoxin system, YefM–YoeB, is thought to be a regulator for antibiotic tolerance, but its physiological role is unknown. The objective of this study was to determine the biofilm and antibiotic susceptibility phenotypes associated with S. aureus yoeB homologs. We hypothesized the toxin–antitoxin yoeB homologs contribute to biofilm formation and antibiotic susceptibility. Disruption of yoeB1 and yoeB2 resulted in decreased biofilm formation in comparison to Newman and JE2 wild-type (WT) S. aureus strains. In comparison to yoeB mutants, both Newman and JE2 WT strains had higher polysaccharide intercellular adhesin (PIA) production. Treatment with sodium metaperiodate increased biofilm formation in Newman WT, indicating biofilm formation may be increased under conditions of oxidative stress. DNase I treatment decreased biofilm formation in Newman WT but not in the absence of yoeB1 or yoeB2. Additionally, WT strains had a higher extracellular DNA (eDNA) content in comparison to yoeB mutants but no differences in biofilm protein content. Moreover, loss of yoeB1 and yoeB2 decreased biofilm survival in both Newman and JE2 strains. Finally, in a neutropenic mouse abscess model, deletion of yoeB1 and yoeB2 resulted in reduced bacterial burden. In conclusion, our data suggest that yoeB1 and yoeB2 are associated with S. aureus planktonic growth, extracellular dependent biofilm formation, antibiotic tolerance, and virulence.

Impacts of type II toxin-antitoxin systems on cell physiology and environmental behavior in acetic acid bacteria

Acetic acid bacteria (AAB) are a group of Gram-negative and strictly aerobic microorganisms widely used in vinegar industry, especially the species belonging to the genera Acetobacter and Komagataeibacter. The environments inhabited by AAB during the vinegar fermentation, in particular those natural traditional bioprocesses, are complex and dynamically changed, usually accompanied by diverse microorganisms, bacteriophages, and the increasing acetic acid concentration. For this reason, how AAB survive to such harsh niches has always been an interesting research field. Previous omic analyses (e.g., genomics, proteomics, and transcriptomics) have provided abundant clues for the metabolic pathways and bioprocesses indispensable for the acid stress adaptation of AAB. Nevertheless, it is far from fully understanding what factors regulate these modular mechanisms overtly and covertly upon shifting environments. Bacterial toxin-antitoxin systems (TAS), usually consisting of a pair of genes encoding a stable toxin and an unstable antitoxin that is capable of counteracting the toxin, have been uncovered to have a variety of biological functions. Recent studies focusing on the role of TAS in Acetobacter pasteurianus suggest that TAS contribute substantially to the acid stress resistance. In this mini review, we discuss the biological functions of type II TAS in the context of AAB with regard to the acid stress resistance, persister formation and resuscitation, genome stability, and phage immunity. KEY POINTS: • Type II TAS act as regulators in the acid stress resistance of AAB. • Type II TAS are implicated in the formation of acid-tolerant persister cells in AAB. • Type II TAS are potential factors responsible for phage immunity and genome stability.

Evaluation of Putative Type II Toxin-Antitoxin Systems and Lon Protease Expression in Shigella flexneri Following Infection of Caco-2 Cells

: Shigella flexneri causes bacillary dysentery in developing countries. Due to recent reports regarding antimicrobial resistance in human S. flexneri, finding alternative therapeutics is of vital importance. Toxin-antitoxin (TA) systems have recently been introduced as antimicrobial targets owing to their involvement in bacterial survival in stress conditions and “persister” cell formation. In this study, the presence of four TA loci were studied in S. flexneri ATCC 12022. The presence of genes coding for the identified TA loci and Lon protease were confirmed by the PCR method using specific primers. Caco-2 cell lines were then infected with this standard strain, and 8 and 24 h post-infection, expression levels of genes coding for the studied TA loci, and Lon protease were evaluated using a real-time PCR method. Expression of mazF, GNAT (Gcn5-related N-acetyltransferase), yeeU, pfam13975, and Lon genes showed 5.4, 9.8, 2.3, 2.7, and 13.8-fold increase, respectively, 8 h after bacterial invasion of the Caco-2 cell line. In addition, the expression of the aforementioned genes showed 4.8, 10.8, 2.3, 3.7, and 16.8-fold increase after 24 h. The GNAT and lon genes showed significantly higher expression levels compared to the control (P value < 0.05). However, the increase in the expression level of yeeU was the same at 8 h and 24 h post-infection. In addition, mazF expression level showed a slight decrease at 24 h compared to 8h post-infection. Genes coding for GNAT and Lon protease showed a significantly higher expression after invading the Caco-2 cell line. Therefore, targeting GNAT or Lon protease can be taken into consideration for finding novel antimicrobial drug strategies. The exact functions and mechanisms of TA systems in S. flexneri isolates are suggested to be experimentally determined.

Mechanistic Modeling of Biochemical Systems without A Priori Parameter Values Using the Design Space Toolbox v.3.0

Mechanistic models of biochemical systems provide a rigorous description of biological phenomena. They are indispensable for making predictions and elucidating biological design principles. To date, mathematical analysis and characterization of these models encounter a bottleneck consisting of large numbers of unknown parameter values. Here, we introduce the Design Space Toolbox v.3.0 (DST3), a software implementation of the Design Space formalism enabling mechanistic modeling without requiring previous knowledge of parameter values. This is achieved by using a phenotype-centric modeling approach, in which the system is first decomposed into a series of biochemical phenotypes. Parameter values realizing phenotypes of interest are subsequently predicted. DST3 represents the most generally applicable implementation of the Design Space formalism and offers unique advantages over earlier versions. By expanding the Design Space formalism and streamlining its distribution, DST3 represents a valuable tool for elucidating biological design principles and designing novel synthetic circuits.

•

DST3 extends the Design Space formalism by identifying additional phenotypes

•

Additional phenotypes arise from cycles, conservations, and metabolic imbalances

•

DST3 enables mechanistic modeling without previous knowledge of parameter values

•

It fully unlocks the potential of the novel phenotype-centric modeling strategy

Role of antibiotic stress in phenotypic switching to persister cells of antibiotic-resistant Staphylococcus aureus

Purpose

This study was designed to evaluate phenotypic and genotypic properties of persister cells formed by Staphylococcus aureus ATCC 15564 (SAWT), oxacillin-induced S. aureus (SAOXA), ciprofloxacin-induced S. aureus (SACIP), and clinically isolated multidrug-resistant S. aureus CCARM 3080 (SAMDR).

Methods

The dose-dependent biphasic killing patterns were observed for SAWT, SAOXA, SACIP, and SAMDR in response to twofold minimum inhibitory concentrate (MIC) of ciprofloxacin. The surviving cells of SAWT, SAOXA, SACIP, and SAMDR after twofold MIC of ciprofloxacin treatment were analyzed using a metabolic-based assay to estimate the fractions of persister cells.

Results

The least persister formation was induced in SACIP after twofold MIC of ciprofloxacin treatment, showing 58% of persistence. The lowest fitness cost of resistance was observed for the recovered persister cells of SACIP (relative fitness = 0.95), followed by SAMDR (relative fitness = 0.70), while the highest fitness cost was observed for SAWT (relative fitness = 0.26). The mRNA transcripts were analyzed by RT-PCR assay in recovered persister cells pre-incubated with ciprofloxacin. The highest expression levels of stress-related genes (dnaK and groEL) and efflux pump-related genes (mepR, norA, and norB) were observed in the recovered persister cells of SAOXA and SAMDR.

Conclusion

This study provides valuable information for understanding crosstalk between antibiotic resistance, tolerance, and persistence in different antibiotic-resistant S. aureus strains.

Deciphering the Antitoxin-Regulated Bacterial Stress Response via Single-Cell Analysis

Bacterial toxin-antitoxin (TA) systems, which are diverse and widespread among prokaryotes, are responsible for tolerance to drugs and environmental stresses. However, the low abundance of toxin and antitoxin proteins renders their quantitative measurement in single bacteria challenging. Employing a laboratory-built nano-flow cytometer (nFCM) to monitor a tetracysteine (TC)-tagged TA system labeled with the biarsenical dye FlAsH, we here report the development of a sensitive method that enables the detection of basal-level expression of antitoxin. Using the Escherichia coli MqsR/MqsA as a model TA system, we reveal for the first time that under its native promoter and in the absence of environmental stress, there exist two populations of bacteria with high or low levels of antitoxin MqsA. Under environmental stress, such as bile acid stress, heat shock, and amino acid starvation, the two populations of bacteria responded differently in terms of MqsA degradation and production. Subsequently, resumed production of MqsA after amino acid stress was observed for the first time. Taking advantage of the multiparameter capability of nFCM, bacterial growth rate and MqsA production were analyzed simultaneously. We found that under environmental stress, the response of bacterial growth was consistent with MqsA production but with an approximate 60 min lag. Overall, the results of the present study indicate that stochastic elevation of MqsA level facilitates bacterial survival, and the two populations with distinct phenotypes empower bacteria to deal with fluctuating environments. This analytical method will help researchers gain deeper insight into the heterogeneity and fundamental role of TA systems.

Antibiotic resilience: a necessary concept to complement antibiotic resistance?

Resilience is the capacity of systems to recover their initial state or functions after a disturbance. The concepts of resilience and resistance are complementary in ecology and both represent different aspects of the stability of ecosystems. However, antibiotic resilience is not used in clinical bacteriology whereas antibiotic resistance is a recognized major problem. To join the fields of ecology and clinical bacteriology, we first review the resilience concept from ecology, socio-ecological systems and microbiology where it is widely developed. We then review resilience-related concepts in microbiology, including bacterial tolerance and persistence, phenotypic heterogeneity and collective tolerance and resistance. We discuss how antibiotic resilience could be defined and argue that the use of this concept largely relies on its experimental measure and its clinical relevance. We review indicators in microbiology which could be used to reflect antibiotic resilience and used as valuable indicators to anticipate the capacity of bacteria to recover from antibiotic treatments.

Stringent response protein as a potential target to intervene persistent bacterial infection

More than half of the world's population is infected with persistent bacterial infections, consequently, persisters are gradually becoming a major public health concern. During the persistent phase, bacterial pathogens deploy many regulatory strategies to compensate unfavorable host environmental conditions. The stringent response is one of such gene regulatory mechanisms which is stimulated by nutrient starvation. It is regulated by the synthesis of highly phosphorylated signaling nucleotides, (p)ppGpp or alarmone. (p)ppGpp is synthesized by ppGpp synthetases, and these proteins are classified as RelA/SpoT homolog (RSH) proteins. Subsequently, (p)ppGpp modulate several molecular and biochemical processes ranging from transcription to metabolism. Imperativeness of (p)ppGpp synthetases has been investigated by numerous approaches including microbiology and animal studies, thereby establishing that Rel enzyme deleted strains of pathogenic bacteria were unable to transform in persister form. In this review, we summarize recent findings to corroborate the rationality to consider (p)ppGpp synthetase as a potential target in discovering a novel class of antimicrobial agents to combat persistent infections. Moreover, inhibition studies on Mycobacterium tuberculosis (p)ppGpp synthetase shows that these inhibitors prevent dormant state transition and biofilm formation. Also, we have highlighted the structural biology of (p)ppGpp synthetases, which may provide significant information that could be used in structure-based inhibitor design.

Identification and first characterization of DinJ-YafQ toxin-antitoxin systems in Lactobacillus species of biotechnological interest

DinJ-YafQ is a type II TA system comprising the ribosome-dependent RNase YafQ toxin and the DinJ antitoxin protein. Although the module has been extensively characterized in Escherichia coli, little information is available for homologous systems in lactic acid bacteria. In this study, we employed bioinformatics tools to identify DinJ-YafQ systems in Lactobacillus casei, Lactobacillus paracasei and Lactobacillus rhamnosus species, commonly used in biotechnological processes. Among a total of nineteen systems found, two TA modules from Lactobacillus paracasei and two modules from Lactobacillus rhamnosus wild strains were isolated and their activity was verified by growth assays in Escherichia coli either in liquid and solid media. The RNase activity of the YafQ toxins was verified in vivo by probing mRNA dynamics and metabolism with single-cell Thioflavin T fluorescence. Our findings demonstrate that, albeit DinJ-YafQ TA systems are widely distributed in lactic acid bacteria, only few are fully functional, while others have lost toxicity even though they maintain high sequence identity with wild type YafQ and a likely functional antitoxin protein.

The relevance of persisters in tuberculosis drug discovery

Bacterial persisters are a subpopulation of cells that exhibit phenotypic resistance during exposure to a lethal dose of antibiotics. They are difficult to target and thought to contribute to the long treatment duration required for tuberculosis. Understanding the molecular and cellular biology of persisters is critical to finding new tuberculosis drugs that shorten treatment. This review focuses on mycobacterial persisters and describes the challenges they pose in tuberculosis therapy, their characteristics and formation, how persistence leads to resistance, and the current approaches being used to target persisters within mycobacterial drug discovery.

A Systematic Overview of Type II and III Toxin-Antitoxin Systems with a Focus on Druggability

Toxin-antitoxin (TA) systems are known to play various roles in physiological processes, such as gene regulation, growth arrest and survival, in bacteria exposed to environmental stress. Type II TA systems comprise natural complexes consisting of protein toxins and antitoxins. Each toxin and antitoxin participates in distinct regulatory mechanisms depending on the type of TA system. Recently, peptides designed by mimicking the interfaces between TA complexes showed its potential to activate the activity of toxin by competing its binding counterparts. Type II TA systems occur more often in pathogenic bacteria than in their nonpathogenic kin. Therefore, they can be possible drug targets, because of their high abundance in some pathogenic bacteria, such as Mycobacterium tuberculosis. In addition, recent bioinformatic analyses have shown that type III TA systems are highly abundant in the intestinal microbiota, and recent clinical studies have shown that the intestinal microbiota is linked to inflammatory diseases, obesity and even several types of cancer. We therefore focused on exploring the putative relationship between intestinal microbiota-related human diseases and type III TA systems. In this paper, we review and discuss the development of possible druggable materials based on the mechanism of type II and type III TA system.

Persister Cells - a Plausible Outcome of Neutral Coevolutionary Drift

The phenomenon of bacterial persistence - a non-inherited antibiotic tolerance in a minute fraction of the bacterial population, was observed more than 70 years ago. Nowadays, it is suggested that "persister cells" undergo an alternative scenario of the cell cycle; however, pathways involved in its emergence are still not identified. We present a mathematically grounded scenario of such possibility. We have determined that population drift in the space of multiple neutrally coupled mutations, which we called "neutrally coupled co-evolution" (NCCE), leads to increased dynamic complexity of bacterial populations via appearance of cells capable of carrying out a single cell cycle in two or more alternative ways and that universal properties of the coupled transcription-translation system underlie this phenotypic multiplicity. According to our hypothesis, modern persister cells have derived from such cells and regulatory mechanisms that govern the consolidation of this phenomenon represented the trigger. We assume that the described type of neutrally coupled co-evolution could play an important role in the origin of extremophiles, both in bacteria and archaea.

Dynamics of Bacterial Gene Regulatory Networks

The ability of bacterial cells to adjust their gene expression program in response to environmental perturbation is often critical for their survival. Recent experimental advances allowing us to quantitatively record gene expression dynamics in single cells and in populations coupled with mathematical modeling enable mechanistic understanding on how these responses are shaped by the underlying regulatory networks. Here, we review how the combination of local and global factors affect dynamical responses of gene regulatory networks. Our goal is to discuss the general principles that allow extrapolation from a few model bacteria to less understood microbes. We emphasize that, in addition to well-studied effects of network architecture, network dynamics are shaped by global pleiotropic effects and cell physiology.

Phenotypic Plasticity, Bet-Hedging, and Androgen Independence in Prostate Cancer: Role of Non-Genetic Heterogeneity

It is well known that genetic mutations can drive drug resistance and lead to tumor relapse. Here, we focus on alternate mechanisms-those without mutations, such as phenotypic plasticity and stochastic cell-to-cell variability that can also evade drug attacks by giving rise to drug-tolerant persisters. The phenomenon of persistence has been well-studied in bacteria and has also recently garnered attention in cancer. We draw a parallel between bacterial persistence and resistance against androgen deprivation therapy in prostate cancer (PCa), the primary standard care for metastatic disease. We illustrate how phenotypic plasticity and consequent mutation-independent or non-genetic heterogeneity possibly driven by protein conformational dynamics can stochastically give rise to androgen independence in PCa, and suggest that dynamic phenotypic plasticity should be considered in devising therapeutic dosing strategies designed to treat and manage PCa.

Division rate, cell size and proteome allocation: impact on gene expression noise and implications for the dynamics of genetic circuits

The cell division rate, size and gene expression programmes change in response to external conditions. These global changes impact on average concentrations of biomolecule and their variability or noise. Gene expression is inherently stochastic, and noise levels of individual proteins depend on synthesis and degradation rates as well as on cell-cycle dynamics. We have modelled stochastic gene expression inside growing and dividing cells to study the effect of division rates on noise in mRNA and protein expression. We use assumptions and parameters relevant to Escherichia coli, for which abundant quantitative data are available. We find that coupling of transcription, but not translation rates to the rate of cell division can result in protein concentration and noise homeostasis across conditions. Interestingly, we find that the increased cell size at fast division rates, observed in E. coli and other unicellular organisms, buffers noise levels even for proteins with decreased expression at faster growth. We then investigate the functional importance of these regulations using gene regulatory networks that exhibit bi-stability and oscillations. We find that network topology affects robustness to changes in division rate in complex and unexpected ways. In particular, a simple model of persistence, based on global physiological feedback, predicts increased proportion of persister cells at slow division rates. Altogether, our study reveals how cell size regulation in response to cell division rate could help controlling gene expression noise. It also highlights that understanding circuits' robustness across growth conditions is key for the effective design of synthetic biological systems.

Division rate, cell size and proteome allocation: impact on gene expression noise and implications for the dynamics of genetic circuits

The cell division rate, size and gene expression programmes change in response to external conditions. These global changes impact on average concentrations of biomolecule and their variability or noise. Gene expression is inherently stochastic, and noise levels of individual proteins depend on synthesis and degradation rates as well as on cell-cycle dynamics. We have modelled stochastic gene expression inside growing and dividing cells to study the effect of division rates on noise in mRNA and protein expression. We use assumptions and parameters relevant to Escherichia coli, for which abundant quantitative data are available. We find that coupling of transcription, but not translation rates to the rate of cell division can result in protein concentration and noise homeostasis across conditions. Interestingly, we find that the increased cell size at fast division rates, observed in E. coli and other unicellular organisms, buffers noise levels even for proteins with decreased expression at faster growth. We then investigate the functional importance of these regulations using gene regulatory networks that exhibit bi-stability and oscillations. We find that network topology affects robustness to changes in division rate in complex and unexpected ways. In particular, a simple model of persistence, based on global physiological feedback, predicts increased proportion of persister cells at slow division rates. Altogether, our study reveals how cell size regulation in response to cell division rate could help controlling gene expression noise. It also highlights that understanding circuits' robustness across growth conditions is key for the effective design of synthetic biological systems.

Elimination of Antibiotic Resistant Surgical Implant Biofilms Using an Engineered Cationic Amphipathic Peptide WLBU2

Antibiotics are unable to remove biofilms from surgical implants. This high antibiotic tolerance is related to bacterial persisters, a sub-population of bacteria phenotypically tolerant to antibiotics secondary to a reduced metabolic state. WLBU2 is an engineered cationic amphipathic peptide designed to maximize antimicrobial activity with minimal mammalian cell toxicity. The objective of this study was to test the ability of WLBU2 to remove Staphylococcus aureus surgical implant biofilms. WLBU2 effectively treated S. aureus biofilms formed by a variety of clinical MSSA and MRSA strains and created culture-negative implants in the in vitro biofilm model. Blocking bacterial metabolism by inhibiting oxidative phosphorylation did not affect WLBU2 killing compared to decreased killing by cefazolin. In the surgical implant infection animal model, WLBU2 decreased biofilm mass as compared to control, untreated samples. WLBU2 could rapidly eliminate implants in vitro and had sufficient efficacy in vivo with minimal systemic toxicity.

msaABCR operon is involved in persister cell formation in Staphylococcus aureus

Background

Persister cells comprise a phenotypic variant that shows extreme antibiotic tolerance resulting in treatment failures of bacterial infections. While this phenomenon has posed a great threat in public health, mechanisms underlying their formation in Staphylococcus aureus remain largely unknown. Increasing evidences of the presence of persister cells in recalcitrant infections underscores the great urgency to unravel the mechanism by which these cells develop. Previously, we characterized msaABCR operon that plays roles in regulation of virulence, biofilm development and antibiotic resistance. We also characterized the function of MsaB protein and showed that MsaB is a putative transcription factor that binds target DNA in response to nutrients availability.

Results

In this study, we compared the number of persister cell in wild type, msaABCR deletion mutant and the complemented strain in two backgrounds USA300 LAC and Mu50. Herein, we report that msaABCR deletion mutant forms significantly less number of persister cells relative to wild type after challenge with various antibiotics in planktonic and biofilm growth conditions. Complementation of the msaABCR operon restored wild type phenotype. Combined antibiotic therapy along with msaABCR deletion significantly improves the killing kinetics of stationary phase and biofilm S. aureus cells. Transcriptomics analysis showed that msaABCR regulates several metabolic genes, transcription factors, transporters and enzymes that may play role in persister cells formation, which we seek to define in the future.

Conclusions

This study presented a new regulator, msaABCR operon, that is involved in the persister cells formation, which is a poorly understood in S. aureus. Indeed, we showed that msaABCR deletion significantly reduces the persister cells formation in all growth phases tested. Although, we have not yet defined the mechanism, we have shown that msaABCR regulates several metabolic, transporters, and extracellular proteases genes that have been previously linked with persister cells formation in other bacterial systems. Taken together, this study showed that inactivation of the msaABCR operon enhances the effectiveness of antibiotics for the treatment of S. aureus infections, especially in context of persister cells.

Electronic supplementary material

The online version of this article (10.1186/s12866-017-1129-9) contains supplementary material, which is available to authorized users.

Synchronized switching of multiple toxin-antitoxin modules by (p)ppGpp fluctuation

Toxin–antitoxin (TA) loci are widespread in bacteria including important pathogenic species. Recent studies suggest that TA systems play a key role in persister formation. However, the persistence phenotype shows only weak dependence on the number of TA systems, i.e. they are functionally redundant. We use a mathematical model to investigate the interaction of multiple TA systems in the switching between growth and persistence. We explore two scenarios: (i) TA systems are bistable and each TA system experiences its own noise and (ii) the noise in the level of common stress signal (e.g. (p)ppGpp) coordinates all TA systems simultaneously. We find that in the first scenario the exit from the persister state strongly depends on the number of TA systems. However in the second case, we could reproduce the weak dependence. The duration of the high (p)ppGpp state was found to be the key parameter for persistence. The (p)ppGpp-driven synchronized transition of all TA systems results in the redundancy.

Fic Proteins of Campylobacter fetus subsp. venerealis Form a Network of Functional Toxin-Antitoxin Systems

Enzymes containing the FIC (filamentation induced by cyclic AMP) domain catalyze post-translational modifications of target proteins. In bacteria the activity of some Fic proteins resembles classical toxin–antitoxin (TA) systems. An excess of toxin over neutralizing antitoxin can enable bacteria to survive some stress conditions by slowing metabolic processes and promoting dormancy. The cell can return to normal growth when sufficient antitoxin is present to block toxin activity. Fic genes of the human and animal pathogen Campylobacter fetus are significantly associated with just one subspecies, which is specifically adapted to the urogenital tract. Here, we demonstrate that the fic genes of virulent isolate C. fetus subsp. venerealis 84-112 form multiple TA systems. Expression of the toxins in Escherichia coli caused filamentation and growth inhibition phenotypes reversible by concomitant antitoxin expression. Key active site residues involved in adenylylation by Fic proteins are conserved in Fic1, Fic3 and Fic4, but degenerated in Fic2. We show that both Fic3 and the non-canonical Fic2 disrupt assembly and function of E. coli ribosomes when expressed independently of a trans-acting antitoxin. Toxicity of the Fic proteins is controlled by different mechanisms. The first involves intramolecular regulation by an inhibitory helix typical for Fic proteins. The second is an unusual neutralization by heterologous Fic–Fic protein interactions. Moreover, a small interacting antitoxin called Fic inhibitory protein 3, which appears unrelated to known Fic antitoxins, has the novel capacity to bind and neutralize Fic toxins encoded in cis and at distant sites. These findings reveal a remarkable system of functional crosstalk occurring between Fic proteins expressed from chromosomal and extrachromosomal modules. Conservation of fic genes in other bacteria that either inhabit or establish pathology in the urogenital tract of humans and animals underscores the significance of these factors for niche-specific adaptation and virulence.

Toxin-Antitoxin Systems: Implications for Plant Disease

Toxin-antitoxin (TA) systems are gene modules that are ubiquitous in free-living prokaryotes. Diverse in structure, cellular function, and fitness roles, TA systems are defined by the presence of a toxin gene that suppresses bacterial growth and a toxin-neutralizing antitoxin gene, usually encoded in a single operon. Originally viewed as DNA maintenance modules, TA systems are now thought to function in many roles, including bacterial stress tolerance, virulence, phage defense, and biofilm formation. However, very few studies have investigated the significance of TA systems in the context of plant-microbe interactions. This review discusses the potential impact and application of TA systems in plant-associated bacteria, guided by insights gained from animal-pathogenic model systems.

Phenotypic Heterogeneity in Mycobacterium tuberculosis

The interaction between the host and the pathogen is extremely complex and is affected by anatomical, physiological, and immunological diversity in the microenvironments, leading to phenotypic diversity of the pathogen. Phenotypic heterogeneity, defined as nongenetic variation observed in individual members of a clonal population, can have beneficial consequences especially in fluctuating stressful environmental conditions. This is all the more relevant in infections caused by Mycobacterium tuberculosis wherein the pathogen is able to survive and often establish a lifelong persistent infection in the host. Recent studies in tuberculosis patients and in animal models have documented the heterogeneous and diverging trajectories of individual lesions within a single host. Since the fate of the individual lesions appears to be determined by the local tissue environment rather than systemic response of the host, studying this heterogeneity is very relevant to ensure better control and complete eradication of the pathogen from individual lesions. The heterogeneous microenvironments greatly enhance M. tuberculosis heterogeneity influencing the growth rates, metabolic potential, stress responses, drug susceptibility, and eventual lesion resolution. Single-cell approaches such as time-lapse microscopy using microfluidic devices allow us to address cell-to-cell variations that are often lost in population-average measurements. In this review, we focus on some of the factors that could be considered as drivers of phenotypic heterogeneity in M. tuberculosis as well as highlight some of the techniques that are useful in addressing this issue.

Relating switching rates between normal and persister cells to substrate and antibiotic concentrations: a mathematical modelling approach supported by experiments

We developed and compared two mathematical models of variable phenotypic switching rates between normal and persister cells that depend on substrate concentration and antibiotic presence. They could be used to simulate the formation of persisters in environments with concentration gradients such as biofilms. Our models are extensions of a previous model of the dynamics of normal and persistent cell populations developed by Balaban et al. (2004, Science 305: 1622). We calibrated the models’ parameters with experimental killing curves obtained after ciprofloxacin treatment of samples regularly harvested from planktonic batch cultures of Klebsiella pneumoniae. Our switching models accurately reproduced the dynamics of normal and persistent populations in planktonic batch cultures and under antibiotic treatment. Results showed that the models are valid for a large range of substrate concentrations and for zero or high doses of antibiotics.

Diversity, Prevalence, and Longitudinal Occurrence of Type II Toxin-Antitoxin Systems of Pseudomonas aeruginosa Infecting Cystic Fibrosis Lungs

Type II toxin-antitoxin (TA) systems are most commonly composed of two genes encoding a stable toxin, which harms the cell, and an unstable antitoxin that can inactivate it. TA systems were initially characterized as selfish elements, but have recently gained attention for regulating general stress responses responsible for pathogen virulence, formation of drug-tolerant persister cells and biofilms—all implicated in causing recalcitrant chronic infections. We use a bioinformatics approach to explore the distribution and evolution of type II TA loci of the opportunistic pathogen, Pseudomonas aeruginosa, across longitudinally sampled isolates from cystic fibrosis lungs. We identify their location in the genome, mutations, and gain/loss during infection to elucidate their function(s) in stabilizing selfish elements and pathogenesis. We found (1) 26 distinct TA systems, where all isolates harbor four in their core genome and a variable number of the remaining 22 on genomic islands; (2) limited mutations in core genome TA loci, suggesting they are not under negative selection; (3) no evidence for horizontal transmission of elements with TA systems between clone types within patients, despite their ability to mobilize; (4) no gain and limited loss of TA-bearing genomic islands, and of those elements partially lost, the remnant regions carry the TA systems supporting their role in genomic stabilization; (5) no significant correlation between frequency of TA systems and strain ability to establish as chronic infection, but those with a particular TA, are more successful in establishing a chronic infection.

Structural conservation of the PIN domain active site across all domains of life

The PIN (PilT N-terminus) domain is a compact RNA-binding protein domain present in all domains of life. This 120-residue domain consists of a central and parallel β sheet surrounded by α helices, which together organize 4-5 acidic residues in an active site that binds one or more divalent metal ions and in many cases has endoribonuclease activity. In bacteria and archaea, the PIN domain is primarily associated with toxin-antitoxin loci, consisting of a toxin (the PIN domain nuclease) and an antitoxin that inhibits the function of the toxin under normal growth conditions. During nutritional or antibiotic stress, the antitoxin is proteolytically degraded causing activation of the PIN domain toxin leading to a dramatic reprogramming of cellular metabolism to cope with the new situation. In eukaryotes, PIN domains are commonly found as parts of larger proteins and are involved in a range of processes involving RNA cleavage, including ribosomal RNA biogenesis and nonsense-mediated mRNA decay. In this review, we provide a comprehensive overview of the structural characteristics of the PIN domain and compare PIN domains from all domains of life in terms of structure, active site architecture, and activity.

A Nonsynonymous SNP Catalog of Mycobacterium tuberculosis Virulence Genes and Its Use for Detecting New Potentially Virulent Sublineages

Mycobacterium tuberculosis is divided into several distinct lineages, and various genetic markers such as IS-elements, VNTR, and SNPs are used for lineage identification. We propose an M. tuberculosis classification approach based on functional polymorphisms in virulence genes. An M. tuberculosis virulence genes catalog has been established, including 319 genes from various protein groups, such as proteases, cell wall proteins, fatty acid and lipid metabolism proteins, sigma factors, toxin–antitoxin systems. Another catalog of 1,573 M. tuberculosis isolates of different lineages has been developed. The developed SNP-calling program has identified 3,563 nonsynonymous SNPs. The constructed SNP-based phylogeny reflected the evolutionary relationship between lineages and detected new sublineages. SNP analysis of sublineage F15/LAM4/KZN revealed four lineage-specific mutations in cyp125, mce3B, vapC25, and vapB34. The Ural lineage has been divided into two geographical clusters based on different SNPs in virulence genes. A new sublineage, B0/N-90, was detected inside the Beijing-B0/W-148 by SNPs in irtB, mce3F and vapC46. We have found 27 members of B0/N-90 among the 227 available genomes of the Beijing-B0/W-148 sublineage. Whole-genome sequencing of strain B9741, isolated from an HIV-positive patient, was demonstrated to belong to the new B0/N-90 group. A primer set for PCR detection of B0/N-90 lineage-specific mutations has been developed. The prospective use of mce3 mutant genes as genetically engineered vaccine is discussed.

Physical and Functional Interplay between MazF1Bif and Its Noncognate Antitoxins from Bifidobacterium longum

Bifidobacterium longum strain JDM301, a widely used commercial strain in China, encodes at least two MazEF-like modules and one RelBE-like toxin-antitoxin (TA) system in its chromosome, designated MazE₁F₁^Bif, MazE₂F₂^Bif, and RelBE^Bif, respectively. Bacterial TA systems play an important role in several stress responses, but the relationship between these TA systems is largely unknown. In this study, the interactions between MazF₁^Bif and MazE₂^Bif or RelB^Bif were assessed in B. longum strain JDM301. MazF₁^Bif caused the degradation of tufA^Bif mRNA, and its toxicity was inhibited by forming a protein complex with its cognate antitoxin, MazE₁^Bif Notably, MazF₁^Bif toxicity was also partially neutralized when jointly expressed with noncognate antitoxin MazE₂^Bif or RelB^Bif Our results show that the two noncognate antitoxins also inhibited mRNA degradation caused by MazF₁^Bif toxin. Furthermore, the physical interplay between MazF₁^Bif and its noncognate antitoxins was confirmed by immunoprecipitation. These results suggest that MazF₁^Bif can arrest cell growth and that MazF₁^Bif toxicity can be neutralized by its cognate and noncognate antitoxins. These results imply that JDM301 uses a sophisticated toxin-antitoxin interaction network to alter its physiology when coping with environmental stress.IMPORTANCE Although toxin-antitoxin (TA) systems play an important role in several stress responses, the regulatory mechanisms of multiple TA system homologs in the bacterial genome remain largely unclear. In this study, the relationships between MazE₁F₁^Bif and the other two TA systems of Bifidobacterium longum strain JDM301 were explored, and the interactions between MazF₁^Bif and MazE₂^Bif or RelB^Bif were characterized. In addition, the mRNA degradation activity of MazF₁^Bif was demonstrated. In particular, the interaction of the toxin with noncognate antitoxins was shown, even between different TA families (MazF₁^Bif toxin and RelB^Bif antitoxin) in JDM301. This work provides insight into the regulatory mechanisms of TA systems implicated in the stress responses of bifidobacteria.

Bifurcation of transition paths induced by coupled bistable systems

We discuss the transition paths in a coupled bistable system consisting of interacting multiple identical bistable motifs. We propose a simple model of coupled bistable gene circuits as an example and show that its transition paths are bifurcating. We then derive a criterion to predict the bifurcation of transition paths in a generalized coupled bistable system. We confirm the validity of the theory for the example system by numerical simulation. We also demonstrate in the example system that, if the steady states of individual gene circuits are not changed by the coupling, the bifurcation pattern is not dependent on the number of gene circuits. We further show that the transition rate exponentially decreases with the number of gene circuits when the transition path does not bifurcate, while a bifurcation facilitates the transition by lowering the quasi-potential energy barrier.