"Gene accordions" cause genotypic and phenotypic heterogeneity in clonal populations of Staphylococcus aureus

Evolution of a bistable genetic system in fluctuating and nonfluctuating environments

Epigenetic mechanisms can generate bacterial lineages capable of spontaneously switching between distinct phenotypes. Currently, mathematical models and simulations propose epigenetic switches as a mechanism of adaptation to deal with fluctuating environments. However, bacterial evolution experiments for testing these predictions are lacking. Here, we exploit an epigenetic switch in Salmonella enterica, the opvAB operon, to show clear evidence that OpvAB bistability persists in changing environments but not in stable conditions. Epigenetic control of transcription in the opvAB operon produces OpvABOFF (phage-sensitive) and OpvABON (phage-resistant) cells in a reversible manner and may be interpreted as an example of bet-hedging to preadapt Salmonella populations to the encounter with phages. Our experimental observations and computational simulations illustrate the adaptive value of epigenetic variation as an evolutionary strategy for mutation avoidance in fluctuating environments. In addition, our study provides experimental support to game theory models predicting that phenotypic heterogeneity is advantageous in changing and unpredictable environments.

SLUSH peptides of the PSMβ family enable Staphylococcus lugdunensis to use erythrocytes as a sole source of nutrient iron

During infection, the host employs nutritional immunity to restrict access to iron. Staphylococcus lugdunensis has been recognized for its ability to utilize host-derived heme to overcome iron restriction. However, the mechanism behind this process involves the release of hemoglobin from erythrocytes, and the hemolytic factors of S. lugdunensis remain poorly understood. S. lugdunensis encodes four phenol-soluble modulins (PSMs), short peptides with hemolytic activity. The peptides SLUSH A, SLUSH B, and SLUSH C are β-type PSMs, and OrfX is an α-type PSM. Our study shows the SLUSH locus to be essential for the hemolytic phenotype of S. lugdunensis. All four peptides individually exhibited hemolytic activity against human and sheep erythrocytes, but synergism with sphingomyelinase was observed exclusively against sheep erythrocytes. Furthermore, our findings demonstrate that SLUSH is crucial for allowing the utilization of erythrocytes as the sole source of nutritional iron and confirm the transcriptional regulation of SLUSH by Agr. Additionally, our study reveals that SLUSH peptides stimulate the human immune system. Our analysis identifies SLUSH as a pivotal hemolytic factor of S. lugdunensis and demonstrates its concerted action with heme acquisition systems to overcome iron limitation in the presence of host erythrocytes.

The Staphylococcus aureus regulatory program in a human skin-like environment

Staphylococcus aureus is a Gram-positive pathogen responsible for the majority of skin and soft tissue infections (SSTIs). S. aureus colonizes the anterior nares of approximately 20%-30% of the population and transiently colonizes the skin, thereby increasing the risk of developing SSTIs and more serious infections. Current laboratory models that mimic the skin surface environment are expensive, require substantial infrastructure, and limit the scope of bacterial physiology studies under human skin conditions. To overcome these limitations, we developed a cost-effective, open-source, chemically defined media recipe termed skin-like medium (SLM) that incorporates key aspects of the human skin surface environment and supports growth of several staphylococcal species. We utilized SLM to investigate the transcriptional response of methicillin-resistant Staphylococcus aureus (MRSA) following growth in SLM compared to a commonly used laboratory media. Through RNA-seq analysis, we observed the upregulation of several virulence factors, including genes encoding functions involved in adhesion, proteolysis, and cytotoxicity. To further explore these findings, we conducted quantitative reverse transcription-PCR (qRT-PCR) experiments to determine the influence of media composition, pH, and temperature on the transcriptional response of key factors involved in adhesion and virulence. We also demonstrated that MRSA primed in SLM adhered better to human corneocytes and demonstrated adhesin-specific phenotypes that previously required genetic manipulation. This improved adherence to corneocytes was dependent on both acidic pH and growth in SLM. These results support the potential utility of SLM as an in vitro model for assessing staphylococcal physiology and metabolism on human skin.

IMPORTANCE

Staphylococcus aureus is the major cause of skin diseases, and its increased prevalence in skin colonization and infections present a need to understand its physiology in this environment. The work presented here outlines S. aureus upregulation of colonization and virulence factors using a newly developed medium that strives to replicate the human skin surface environment and demonstrates roles for adhesins clumping factor A (ClfA), serine-rich repeat glycoprotein adhesin (SraP), and the fibronectin binding proteins (Fnbps) in human corneocyte adherence.

Direct long read visualization reveals metabolic interplay between two antimalarial drug targets

Increases in the copy number of large genomic regions, termed genome amplification, are an important adaptive strategy for malaria parasites. Numerous amplifications across the Plasmodium falciparum genome contribute directly to drug resistance or impact the fitness of this protozoan parasite. During the characterization of parasite lines with amplifications of the dihydroorotate dehydrogenase (DHODH) gene, we detected increased copies of an additional genomic region that encompassed 3 genes (~5 kb) including GTP cyclohydrolase I (GCH1 amplicon). While this gene is reported to increase the fitness of antifolate resistant parasites, GCH1 amplicons had not previously been implicated in any other antimalarial resistance context. Here, we further explored the association between GCH1 and DHODH copy number. Using long read sequencing and single read visualization, we directly observed a higher number of tandem GCH1 amplicons in parasites with increased DHODH copies (up to 9 amplicons) compared to parental parasites (3 amplicons). While all GCH1 amplicons shared a consistent structure, expansions arose in 2-unit steps (from 3 to 5 to 7, etc copies). Adaptive evolution of DHODH and GCH1 loci was further bolstered when we evaluated prior selection experiments; DHODH amplification was only successful in parasite lines with pre-existing GCH1 amplicons. These observations, combined with the direct connection between metabolic pathways that contain these enzymes, lead us to propose that the GCH1 locus is beneficial for the fitness of parasites exposed to DHODH inhibitors. This finding highlights the importance of studying variation within individual parasite genomes as well as biochemical connections of drug targets as novel antimalarials move towards clinical approval.

A type VII-secreted lipase toxin with reverse domain arrangement

The type VII protein secretion system (T7SS) is found in many Gram-positive bacteria and in pathogenic mycobacteria. All T7SS substrate proteins described to date share a common helical domain architecture at the N-terminus that typically interacts with other helical partner proteins, forming a composite signal sequence for targeting to the T7SS. The C-terminal domains are functionally diverse and in Gram-positive bacteria such as Staphylococcus aureus often specify toxic anti-bacterial activity. Here we describe the first example of a class of T7 substrate, TslA, that has a reverse domain organisation. TslA is widely found across Bacillota including Staphylococcus, Enterococcus and Listeria. We show that the S. aureus TslA N-terminal domain is a phospholipase A with anti-staphylococcal activity that is neutralised by the immunity lipoprotein TilA. Two small helical partner proteins, TlaA1 and TlaA2 are essential for T7-dependent secretion of TslA and at least one of these interacts with the TslA C-terminal domain to form a helical stack. Cryo-EM analysis of purified TslA complexes indicate that they share structural similarity with canonical T7 substrates. Our findings suggest that the T7SS has the capacity to recognise a secretion signal present at either end of a substrate.

Genomic insights into linezolid-resistant Enterococci revealed its evolutionary diversity and poxtA copy number heterogeneity

OBJECTIVES

This study aimed to determine the molecular mechanisms of linezolid-resistant enterococci (LRE) in swine slaughterhouses in China and apply the "One Health" perspective to analyse the evolutionary dynamics of poxtA-positive E. faecium in clinical and non-clinical settings worldwide.

METHODS

The phenotypic and genomic characteristics of multiple LRE isolates were systematically investigated using antimicrobial susceptibility testing, transfer assays, evolutionary experiments, quantitative RT-PCR assays, whole-genome sequencing, and bioinformatics analyses.

RESULTS

Swine faeces served as a significant reservoir for LRE isolates, and optrA and poxtA were the primary contributors to linezolid resistance. Co-occurrence network analysis revealed a significant interconnection between optrA and several other ARGs. The poxtA copy number heterogeneity and polymorphism were initially observed in E. faecium parental and evolved isolates. The poxtA-carrying tandem repeat region exhibits high mobility and has undergone extensive duplication owing to linezolid pressure. The poxtA copy number varies from four copies on the plasmid of E. faecium IC25 to 11 copies on the plasmid and six copies on the chromosome in the evolved isolate IC25-50_poxtA. Furthermore, phylogenetic analysis of 185 poxtA-positive E. faecium strains worldwide found that one isolate from a French patient in 2018 shared only two SNPs with CC17 E. faecium isolates IC25 and IC7-2 from this study, highlighting the potential global transmission of CC17 poxtA-positive E. faecium between humans and animals.

CONCLUSION

This study identified amplification of poxtA as a response of E. faecium to linezolid pressure. Phylogenetic analysis shed light on the potential global transmission of hospital-associated CC17 poxtA-positive E. faecium in clinical and non-clinical settings.

Diverse mechanisms of bioproduction heterogeneity in fermentation and their control strategies

Microbial bioproduction often faces challenges related to populational heterogeneity, where cells exhibit varying biosynthesis capabilities. Bioproduction heterogeneity can stem from genetic and non-genetic factors, resulting in decreased titer, yield, stability, and reproducibility. Consequently, understanding and controlling bioproduction heterogeneity are crucial for enhancing the economic competitiveness of large-scale biomanufacturing. In this review, we provide a comprehensive overview of current understandings of the various mechanisms underlying bioproduction heterogeneity. Additionally, we examine common strategies for controlling bioproduction heterogeneity based on these mechanisms. By implementing more robust measures to mitigate heterogeneity, we anticipate substantial enhancements in the scalability and stability of bioproduction processes.

ONE-SENTENCE SUMMARY

This review summarizes current understandings of different mechanisms of bioproduction heterogeneity and common control strategies based on these mechanisms.

Antibiotics Limit Adaptation of Drug-Resistant Staphylococcus aureus to Hypoxia

Bacterial pathogens are confronted with a range of challenges at the site of infection, including exposure to antibiotic treatment and harsh physiological conditions, that can alter the fitness benefits and costs of acquiring antibiotic resistance. Here, we develop an experimental system to recapitulate resistance gene acquisition by Staphylococcus aureus and test how the subsequent evolution of the resistant bacterium is modulated by antibiotic treatment and oxygen levels, both of which are known to vary extensively at sites of infection. We show that acquiring tetracycline resistance was costly, reducing competitive growth against the isogenic strain without the resistance gene in the absence of the antibiotic, for S. aureus under hypoxic but not normoxic conditions. Treatment with tetracycline or doxycycline drove the emergence of enhanced resistance through mutations in an RluD-like protein-encoding gene and duplications of tetL, encoding the acquired tetracycline-specific efflux pump. In contrast, evolutionary adaptation by S. aureus to hypoxic conditions, which evolved in the absence of antibiotics through mutations affecting gyrB, was impeded by antibiotic treatment. Together, these data suggest that the horizontal acquisition of a new resistance mechanism is merely a starting point for the emergence of high-level resistance under antibiotic selection but that antibiotic treatment constrains pathogen adaptation to other important environmental selective forces such as hypoxia, which in turn could limit the survival of these highly resistant but poorly adapted genotypes after antibiotic treatment is ended.

From genome structure to function: insights into structural variation in microbiology

Structural variation in bacterial genomes is an important evolutionary driver. Genomic rearrangements, such as inversions, duplications, and insertions, can regulate gene expression and promote niche adaptation. Importantly, many of these variations are reversible and preprogrammed to generate heterogeneity. While many tools have been developed to detect structural variation in eukaryotic genomes, variation in bacterial genomes and metagenomes remains understudied. However, recent advances in genome sequencing technology and the development of new bioinformatic pipelines hold promise in further understanding microbial genomics.

Gene copy number variations at the within-host population level modulate gene expression in a multipartite virus

Multipartite viruses have a segmented genome, with each segment encapsidated separately. In all multipartite virus species for which the question has been addressed, the distinct segments reproducibly accumulate at a specific and host-dependent relative frequency, defined as the 'genome formula'. Here, we test the hypothesis that the multipartite genome organization facilitates the regulation of gene expression via changes of the genome formula and thus via gene copy number variations. In a first experiment, the faba bean necrotic stunt virus (FBNSV), whose genome is composed of eight DNA segments each encoding a single gene, was inoculated into faba bean or alfalfa host plants, and the relative concentrations of the DNA segments and their corresponding messenger RNAs (mRNAs) were monitored. In each of the two host species, our analysis consistently showed that the genome formula variations modulate gene expression, the concentration of each genome segment linearly and positively correlating to that of its cognate mRNA but not of the others. In a second experiment, twenty parallel FBNSV lines were transferred from faba bean to alfalfa plants. Upon host switching, the transcription rate of some genome segments changes, but the genome formula is modified in a way that compensates for these changes and maintains a similar ratio between the various viral mRNAs. Interestingly, a deep-sequencing analysis of these twenty FBNSV lineages demonstrated that the host-related genome formula shift operates independently of DNA-segment sequence mutation. Together, our results indicate that nanoviruses are plastic genetic systems, able to transiently adjust gene expression at the population level in changing environments, by modulating the copy number but not the sequence of each of their genes.

A high-quality reference genome for the fish pathogen Streptococcus iniae

Fish mortality caused by Streptococcus iniae is a major economic problem in aquaculture in warm and temperate regions globally. There is also risk of zoonotic infection by S. iniae through handling of contaminated fish. In this study, we present the complete genome sequence of S. iniae strain QMA0248, isolated from farmed barramundi in South Australia. The 2.12 Mb genome of S. iniae QMA0248 carries a 32 kb prophage, a 12 kb genomic island and 92 discrete insertion sequence (IS) elements. These include nine novel IS types that belong mostly to the IS3 family. Comparative and phylogenetic analysis between S. iniae QMA0248 and publicly available complete S. iniae genomes revealed discrepancies that are probably due to misassembly in the genomes of isolates ISET0901 and ISNO. Long-range PCR confirmed five rRNA loci in the PacBio assembly of QMA0248, and, unlike S. iniae 89353, no tandemly repeated rRNA loci in the consensus genome. However, we found sequence read evidence that the tandem rRNA repeat existed within a subpopulation of the original QMA0248 culture. Subsequent nanopore sequencing revealed that the tandem rRNA repeat was the most prevalent genotype, suggesting that there is selective pressure to maintain fewer rRNA copies under uncertain laboratory conditions. Our study not only highlights assembly problems in existing genomes, but provides a high-quality reference genome for S. iniae QMA0248, including manually curated mobile genetic elements, that will assist future S. iniae comparative genomic and evolutionary studies.

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.

Genomic analysis of the progenitor strains of Staphylococcus aureus RN6390

RN6390 is a commonly used laboratory strain of Staphylococcus aureus derived from NCTC8325. In this study, we sequenced the RN6390 genome and compared it to available genome sequences for NCTC8325. We confirmed that three prophages, Φ11, Φ12 and Φ13, which are present in NCTC8325 are absent from the genome of RN6390, consistent with the successive curing events leading to the generation of this strain. However, we noted that a separate prophage is present in RN6390 that is not found in NCTC8325. Two separate genome sequences have been deposited for the parental strain, NCTC8325. Analysis revealed several differences between these sequences, in particular, between the copy number of esaG genes, which encode immunity proteins to the type VII secreted anti-bacterial toxin, EsaD. Single nucleotide polymorphisms were also detected in ribosomal RNA genes and genes encoding microbial surface components recognizing adhesive matrix molecules (MSCRAMM) between the two NCTC8325 sequences. Comparing each NCTC8325 sequence to other strains in the RN6390 lineage confirmed that sequence assembly errors in the earlier NCTC8325 sequence are the most likely explanation for most of the differences observed.

Adaptation dynamics between copy-number and point mutations

Together, copy-number and point mutations form the basis for most evolutionary novelty, through the process of gene duplication and divergence. While a plethora of genomic data reveals the long-term fate of diverging coding sequences and their cis-regulatory elements, little is known about the early dynamics around the duplication event itself. In microorganisms, selection for increased gene expression often drives the expansion of gene copy-number mutations, which serves as a crude adaptation, prior to divergence through refining point mutations. Using a simple synthetic genetic reporter system that can distinguish between copy-number and point mutations, we study their early and transient adaptive dynamics in real time in Escherichia coli. We find two qualitatively different routes of adaptation, depending on the level of functional improvement needed. In conditions of high gene expression demand, the two mutation types occur as a combination. However, under low gene expression demand, copy-number and point mutations are mutually exclusive; here, owing to their higher frequency, adaptation is dominated by copy-number mutations, in a process we term amplification hindrance. Ultimately, due to high reversal rates and pleiotropic cost, copy-number mutations may not only serve as a crude and transient adaptation, but also constrain sequence divergence over evolutionary time scales.

Effect of genetic background on the evolution of Vancomycin-Intermediate Staphylococcus aureus (VISA)

Vancomycin-intermediate Staphylococcus aureus (VISA) typically arises through accumulation of chromosomal mutations that alter cell-wall thickness and global regulatory pathways. Genome-based prediction of VISA requires understanding whether strain background influences patterns of mutation that lead to resistance. We used an iterative method to experimentally evolve three important methicillin-resistant S. aureus (MRSA) strain backgrounds-(CC1, CC5 and CC8 (USA300)) to generate a library of 120 laboratory selected VISA isolates. At the endpoint, isolates had vancomycin MICs ranging from 4 to 10 μg/mL. We detected mutations in more than 150 genes, but only six genes (already known to be associated with VISA from prior studies) were mutated in all three background strains (walK, prs, rpoB, rpoC, vraS, yvqF). We found evidence of interactions between loci (e.g., vraS and yvqF mutants were significantly negatively correlated) and rpoB, rpoC, vraS and yvqF were more frequently mutated in one of the backgrounds. Increasing vancomycin resistance was correlated with lower maximal growth rates (a proxy for fitness) regardless of background. However, CC5 VISA isolates had higher MICs with fewer rounds of selection and had lower fitness costs than the CC8 VISA isolates. Using multivariable regression, we found that genes differed in their contribution to overall MIC depending on the background. Overall, these results demonstrated that VISA evolved through mutations in a similar set of loci in all backgrounds, but the effect of mutation in common genes differed with regard to fitness and contribution to resistance in different strains.

Waddington's Landscapes in the Bacterial World

Conrad Waddington’s epigenetic landscape, a visual metaphor for the development of multicellular organisms, is appropriate to depict the formation of phenotypic variants of bacterial cells. Examples of bacterial differentiation that result in morphological change have been known for decades. In addition, bacterial populations contain phenotypic cell variants that lack morphological change, and the advent of fluorescent protein technology and single-cell analysis has unveiled scores of examples. Cell-specific gene expression patterns can have a random origin or arise as a programmed event. When phenotypic cell-to-cell differences are heritable, bacterial lineages are formed. The mechanisms that transmit epigenetic states to daughter cells can have strikingly different levels of complexity, from the propagation of simple feedback loops to the formation of complex DNA methylation patterns. Game theory predicts that phenotypic heterogeneity can facilitate bacterial adaptation to hostile or unpredictable environments, serving either as a division of labor or as a bet hedging that anticipates future challenges. Experimental observation confirms the existence of both types of strategies in the bacterial world.