Phenotypic heterogeneity promotes adaptive evolution

Climate warming promotes collateral antibiotic resistance development in cyanobacteria

Both cyanobacterial blooms and antibiotic resistance have aggravated worldwide and posed a great threat to public health in recent years. As a significant source and reservoir of water environmental resistome, cyanobacteria exhibit confusing discrepancy between their reduced susceptibility and their chronic exposure to antibiotic mixtures at sub-inhibitory concentrations. How the increasing temperature affects the adaptive evolution of cyanobacteria-associated antibiotic resistance in response to low-level antibiotic combinations under climate change remains unclear. Here we profiled the antibiotic interaction and collateral susceptibility networks among 33 commonly detected antibiotics in 600 cyanobacterial strains isolated from 50 sites across four eutrophicated lakes in China. Cyanobacteria-associated antibiotic resistance level was found positively correlated to antibiotic heterogeneity across all sites. Among 528 antibiotic combinations, antagonism was observed for 62 % interactions and highly conserved within cyanobacterial species. Collateral resistance was detected in 78.5 % of pairwise antibiotic interaction, leading to a widened or shifted upwards mutant selection window for increased opportunity of acquiring second-step mutations. We quantified the interactive promoting effect of collateral resistance and increasing temperature on the evolution of both phenotypic and genotypic cyanobacteria-associated resistance under chronic exposure to environmental level of antibiotic combinations. With temperature increasing from 16 °C to 36 °C, the evolvability index and genotypic resistance level increased by 1.25 - 2.5 folds and 3 - 295 folds in the collateral-resistance-informed lineages, respectively. Emergence of resistance mutation pioneered by tolerance, which was jointly driven by mutation rate and persister fraction, was found to be accelerated by increased temperature and antibiotic switching rate. Our findings provided mechanic insights into the boosting effect of climate warming on the emergence and development of cyanobacteria-associated resistance against collateral antibiotic phenotypes.

Relationships between genomic dissipation and de novo SNP evolution

Patterns of single nucleotide polymorphisms (SNPs) in eukaryotic DNA are traditionally attributed to selective pressure, drift, identity descent, or related factors-without accounting for ways in which bias during de novo SNP formation, itself, might contribute. A functional and phenotypic analysis based on evolutionary resilience of DNA points to decreased numbers of non-synonymous SNPs in human and other genomes, with a predominant component of SNP depletion in the human gene pool caused by robust preferences during de novo SNP formation (rather than selective constraint). Ramifications of these findings are broad, belie a number of concepts regarding human evolution, and point to a novel interpretation of evolving DNA across diverse species.

Nonlethal deleterious mutation-induced stress accelerates bacterial aging

Random mutagenesis, including when it leads to loss of gene function, is a key mechanism enabling microorganisms' long-term adaptation to new environments. However, loss-of-function mutations are often deleterious, triggering, in turn, cellular stress and complex homeostatic stress responses, called "allostasis," to promote cell survival. Here, we characterize the differential impacts of 65 nonlethal, deleterious single-gene deletions on Escherichia coli growth in three different growth environments. Further assessments of select mutants, namely, those bearing single adenosine triphosphate (ATP) synthase subunit deletions, reveal that mutants display reorganized transcriptome profiles that reflect both the environment and the specific gene deletion. We also find that ATP synthase α-subunit deleted (ΔatpA) cells exhibit elevated metabolic rates while having slower growth compared to wild-type (wt) E. coli cells. At the single-cell level, compared to wt cells, individual ΔatpA cells display near normal proliferation profiles but enter a postreplicative state earlier and exhibit a distinct senescence phenotype. These results highlight the complex interplay between genomic diversity, adaptation, and stress response and uncover an "aging cost" to individual bacterial cells for maintaining population-level resilience to environmental and genetic stress; they also suggest potential bacteriostatic antibiotic targets and -as select human genetic diseases display highly similar phenotypes, - a bacterial origin of some human diseases.

Deep-learning-assisted Sort-Seq enables high-throughput profiling of gene expression characteristics with high precision

Owing to the nondeterministic and nonlinear nature of gene expression, the steady-state intracellular protein abundance of a clonal population forms a distribution. The characteristics of this distribution, including expression strength and noise, are closely related to cellular behavior. However, quantitative description of these characteristics has so far relied on arrayed methods, which are time-consuming and labor-intensive. To address this issue, we propose a deep-learning-assisted Sort-Seq approach (dSort-Seq) in this work, enabling high-throughput profiling of expression properties with high precision. We demonstrated the validity of dSort-Seq for large-scale assaying of the dose-response relationships of biosensors. In addition, we comprehensively investigated the contribution of transcription and translation to noise production in Escherichia coli, from which we found that the expression noise is strongly coupled with the mean expression level. We also found that the transcriptional interference caused by overlapping RpoD-binding sites contributes to noise production, which suggested the existence of a simple and feasible noise control strategy in E. coli.

Polygenic pathogen networks influence transcriptional plasticity in the Arabidopsis-Botrytis pathosystem

Bidirectional flow of information shapes the outcome of the host-pathogen interactions and depends on the genetics of each organism. Recent work has begun to use co-transcriptomic studies to shed light on this bidirectional flow, but it is unclear how plastic the co-transcriptome is in response to genetic variation in both the host and pathogen. To study co-transcriptome plasticity, we conducted transcriptomics using natural genetic variation in the pathogen, Botrytis cinerea, and large-effect genetic variation abolishing defense signaling pathways within the host, Arabidopsis thaliana. We show that genetic variation in the pathogen has a greater influence on the co-transcriptome than mutations that abolish defense signaling pathways in the host. Genome-wide association mapping using the pathogens' genetic variation and both organisms' transcriptomes allowed an assessment of how the pathogen modulates plasticity in response to the host. This showed that the differences in both organism's responses were linked to trans-expression quantitative trait loci (eQTL) hotspots within the pathogen's genome. These hotspots control gene sets in either the host or pathogen and show differential allele sensitivity to the host's genetic variation rather than qualitative host specificity. Interestingly, nearly all the trans-eQTL hotspots were unique to the host or pathogen transcriptomes. In this system of differential plasticity, the pathogen mediates the shift in the co-transcriptome more than the host.

Nonmonotone invasion landscape by noise-aware control of metastasis activator levels

A major pharmacological assumption is that lowering disease-promoting protein levels is generally beneficial. For example, inhibiting metastasis activator BACH1 is proposed to decrease cancer metastases. Testing such assumptions requires approaches to measure disease phenotypes while precisely adjusting disease-promoting protein levels. Here we developed a two-step strategy to integrate protein-level tuning, noise-aware synthetic gene circuits into a well-defined human genomic safe harbor locus. Unexpectedly, engineered MDA-MB-231 metastatic human breast cancer cells become more, then less and then more invasive as we tune BACH1 levels up, irrespective of the native BACH1. BACH1 expression shifts in invading cells, and expression of BACH1's transcriptional targets confirm BACH1's nonmonotone phenotypic and regulatory effects. Thus, chemical inhibition of BACH1 could have unwanted effects on invasion. Additionally, BACH1's expression variability aids invasion at high BACH1 expression. Overall, precisely engineered, noise-aware protein-level control is necessary and important to unravel disease effects of genes to improve clinical drug efficacy.

Dynamic proteome trade-offs regulate bacterial cell size and growth in fluctuating nutrient environments

Bacteria dynamically regulate cell size and growth to thrive in changing environments. While previous studies have characterized bacterial growth physiology at steady-state, a quantitative understanding of bacterial physiology in time-varying environments is lacking. Here we develop a quantitative theory connecting bacterial growth and division rates to proteome allocation in time-varying nutrient environments. In such environments, cell size and growth are regulated by trade-offs between prioritization of biomass accumulation or division, resulting in decoupling of single-cell growth rate from population growth rate. Specifically, bacteria transiently prioritize biomass accumulation over production of division machinery during nutrient upshifts, while prioritizing division over growth during downshifts. When subjected to pulsatile nutrient concentration, we find that bacteria exhibit a transient memory of previous metabolic states due to the slow dynamics of proteome reallocation. This allows for faster adaptation to previously seen environments and results in division control which is dependent on the time-profile of fluctuations.

Frank Beach Award Winner: The centrality of the hypothalamic-pituitary-adrenal axis in dealing with environmental change across temporal scales

Understanding if and how individuals and populations cope with environmental change is an enduring question in evolutionary ecology that has renewed importance given the pace of change in the Anthropocene. Two evolutionary strategies of coping with environmental change may be particularly important in rapidly changing environments: adaptive phenotypic plasticity and/or bet hedging. Adaptive plasticity could enable individuals to match their phenotypes to the expected environment if there is an accurate cue predicting the selective environment. Diversifying bet hedging involves the production of seemingly random phenotypes in an unpredictable environment, some of which may be adaptive. Here, I review the central role of the hypothalamic-pituitary-adrenal (HPA) axis and glucocorticoids (GCs) in enabling vertebrates to cope with environmental change through adaptive plasticity and bet hedging. I first describe how the HPA axis mediates three types of adaptive plasticity to cope with environmental change (evasion, tolerance, recovery) over short timescales (e.g., 1-3 generations) before discussing how the implications of GCs on phenotype integration may depend upon the timescale under consideration. GCs can promote adaptive phenotypic integration, but their effects on phenotypic co-variation could also limit the dimensions of phenotypic space explored by animals over longer timescales. Finally, I discuss how organismal responses to environmental stressors can act as a bet hedging mechanism and therefore enhance evolvability by increasing genetic or phenotypic variability or reducing patterns of genetic and phenotypic co-variance. Together, this emphasizes the crucial role of the HPA axis in understanding fundamental questions in evolutionary ecology.

Metabolic heterogeneity and cross-feeding within isogenic yeast populations captured by DILAC

Genetically identical cells are known to differ in many physiological parameters such as growth rate and drug tolerance. Metabolic specialization is believed to be a cause of such phenotypic heterogeneity, but detection of metabolically divergent subpopulations remains technically challenging. We developed a proteomics-based technology, termed differential isotope labelling by amino acids (DILAC), that can detect producer and consumer subpopulations of a particular amino acid within an isogenic cell population by monitoring peptides with multiple occurrences of the amino acid. We reveal that young, morphologically undifferentiated yeast colonies contain subpopulations of lysine producers and consumers that emerge due to nutrient gradients. Deconvoluting their proteomes using DILAC, we find evidence for in situ cross-feeding where rapidly growing cells ferment and provide the more slowly growing, respiring cells with ethanol. Finally, by combining DILAC with fluorescence-activated cell sorting, we show that the metabolic subpopulations diverge phenotypically, as exemplified by a different tolerance to the antifungal drug amphotericin B. Overall, DILAC captures previously unnoticed metabolic heterogeneity and provides experimental evidence for the role of metabolic specialization and cross-feeding interactions as a source of phenotypic heterogeneity in isogenic cell populations.

Exploring the expanse between theoretical questions and experimental approaches in the modern study of evolvability

Despite several decades of computational and experimental work across many systems, evolvability remains on the periphery with regards to its status as a widely accepted and regularly applied theoretical concept. Here we propose that its marginal status is partly a result of large gaps between the diverse but disconnected theoretical treatments of evolvability and the relatively narrower range of studies that have tested it empirically. To make this case, we draw on a range of examples-from experimental evolution in microbes, to molecular evolution in proteins-where attempts have been made to mend this disconnect. We highlight some examples of progress that has been made and point to areas where synthesis and translation of existing theory can lead to further progress in the still-new field of empirical measurements of evolvability.

Phenotypic plasticity as a facilitator of microbial evolution

Tossed about by the tides of history, the inheritance of acquired characteristics has found a safe harbor at last in the rapidly expanding field of epigenetics. The slow pace of genetic variation and high opportunity cost associated with maintaining a diverse genetic pool are well-matched by the flexibility of epigenetic traits, which can enable low-cost exploration of phenotypic space and reactive tuning to environmental pressures. Aiding in the generation of a phenotypically plastic population, epigenetic mechanisms often provide a hotbed of innovation for countering environmental pressures, while the potential for genetic fixation can lead to strong epigenetic-genetic evolutionary synergy. At the level of cells and cellular populations, we begin this review by exploring the breadth of mechanisms for the storage and intergenerational transmission of epigenetic information, followed by a brief review of common and exotic epigenetically regulated phenotypes. We conclude by offering an in-depth coverage of recent papers centered around two critical issues: the evolvability of epigenetic traits through Baldwinian adaptive phenotypic plasticity and the potential for synergy between epigenetic and genetic evolution.

Horizontal gene transfer enables programmable gene stability in synthetic microbiota

The functions of many microbial communities exhibit remarkable stability despite fluctuations in the compositions of these communities. To date, a mechanistic understanding of this function-composition decoupling is lacking. Statistical mechanisms have been commonly hypothesized to explain such decoupling. Here, we proposed that dynamic mechanisms, mediated by horizontal gene transfer (HGT), also enable the independence of functions from the compositions of microbial communities. We combined theoretical analysis with numerical simulations to illustrate that HGT rates can determine the stability of gene abundance in microbial communities. We further validated these predictions using engineered microbial consortia of different complexities transferring one or more than a dozen clinically isolated plasmids, as well as through the reanalysis of data from the literature. Our results demonstrate a generalizable strategy to program the gene stability of microbial communities.

The role of the Aspergillus nidulans high mobility group B protein HmbA, the orthologue of Saccharomyces cerevisiae Nhp6p

The mammalian HMGB1 is a high-mobility-group B protein, which is both an architectural and functional element of chromatin. Nhp6p, the extensively studied fungal homologue of HMGB1 in Saccharomyces cerevisiae has pleiotropic physiological functions. Despite the existence of Nhp6p orthologues in filamentous ascomycetes, little is known about their physiological roles besides their contribution to sexual development. Here we study the function of HmbA, the Aspergillus nidulans orthologue of Nhp6p. We show that HmbA influences the utilization of various carbon- and nitrogen sources, stress tolerance, secondary metabolism, hyphae elongation and maintenance of polarized growth. Additionally, by conducting heterologous expression studies, we demonstrate that HmbA and Nhp6p are partially interchangeable. HmbA restores SNR6 transcription and fitness of nhp6AΔBΔ mutant and reverses its heat sensitivity. Nhp6Ap complements several phenotypes of hmbAΔ, including ascospore formation, utilization of various carbon- and nitrogen-sources, radial growth rate, hypha elongation by polarized growth. However, Nhp6Ap does not complement sterigmatocystin production in a hmbAΔ strain. Finally, we also show that HmbA is necessary for the normal expression of the endochitinase chiA, a cell wall re-modeller that is pivotal for the normal mode of maintenance of polar growth.

Phenotypic selection during laboratory evolution of yeast populations leads to a genome-wide sustainable chromatin compaction shift

In a previous study, we have shown how microbial evolution has resulted in a persistent reduction in expression after repeatedly selecting for the lowest PGAL1-YFP-expressing cells. Applying the ATAC-seq assay on samples collected from this 28-day evolution experiment, here we show how genome-wide chromatin compaction changes during evolution under selection pressure. We found that the chromatin compaction was altered not only on GAL network genes directly impacted by the selection pressure, showing an example of selection-induced non-genetic memory, but also at the whole-genome level. The GAL network genes experienced chromatin compaction accompanying the reduction in PGAL1-YFP reporter expression. Strikingly, the fraction of global genes with differentially compacted chromatin states accounted for about a quarter of the total genome. Moreover, some of the ATAC-seq peaks followed well-defined temporal dynamics. Comparing peak intensity changes on consecutive days, we found most of the differential compaction to occur between days 0 and 3 when the selection pressure was first applied, and between days 7 and 10 when the pressure was lifted. Among the gene sets enriched for the differential compaction events, some had increased chromatin availability once selection pressure was applied and decreased availability after the pressure was lifted (or vice versa). These results intriguingly show that, despite the lack of targeted selection, transcriptional availability of a large fraction of the genome changes in a very diverse manner during evolution, and these changes can occur in a relatively short number of generations.

A Trait-Based Approach to Predicting Viral Host-Range Evolvability

Predicting the evolution of virus host range has proven to be extremely difficult, in part because of the sheer diversity of viruses, each with unique biology and ecological interactions. We have not solved this problem, but to make the problem more tractable, we narrowed our focus to three traits intrinsic to all viruses that may play a role in host-range evolvability: mutation rate, recombination rate, and phenotypic heterogeneity. Although each trait should increase evolvability, they cannot do so unbounded because fitness trade-offs limit the ability of all three traits to maximize evolvability. By examining these constraints, we can begin to identify groups of viruses with suites of traits that make them especially concerning, as well as ecological and environmental conditions that might push evolution toward accelerating host-range expansion.

Influence of Different Aggregation States on Volatile Organic Compounds Released by Dairy Kluyveromyces marxianus Strains

Kluyveromyces marxianus has the ability to contribute to the aroma profile of foods and beverages since it is able to produce several volatile organic compounds (VOCs). In this study, 8 dairy K. marxianus strains, previously selected for their adhesion properties, were tested for VOCs production when grown in different conditions: planktonic, biofilm-detached, and MATS forming-cells. It was shown that biofilm-detached cells were mainly able to produce higher alcohols (64.57 mg/L), while esters were mainly produced by planktonic and MATS forming-cells (117.86 and 94.90 mg/L, respectively). Moreover, K. marxianus biofilm-detached cells were able to produce VOCs with flavor and odor impacts, such as ketons, phenols, and terpenes, which were not produced by planktonic cells. In addition, specific unique compounds were associated to the different conditions tested. Biofilm-detached cells were characterized by the production of 9 unique compounds, while planktonic and MATS forming-cells by 7 and 12, respectively. The obtained results should be exploited to modulate the volatilome of foods and beverages and improve the production of certain compounds at the industrial level. Further studies will be carried out to better understand the genetic mechanisms underlying the metabolic pathways activated under different conditions.

Inheritance of somatic mutations by animal offspring

Since 1892, it has been widely assumed that somatic mutations are evolutionarily irrelevant in animals because they cannot be inherited by offspring. However, some nonbilaterians segregate the soma and germline late in development or never, leaving the evolutionary fate of their somatic mutations unknown. By investigating uni- and biparental reproduction in the coral Acropora palmata (Cnidaria, Anthozoa), we found that uniparental, meiotic offspring harbored 50% of the 268 somatic mutations present in their parent. Thus, somatic mutations accumulated in adult coral animals, entered the germline, and were passed on to swimming larvae that grew into healthy juvenile corals. In this way, somatic mutations can increase allelic diversity and facilitate adaptation across habitats and generations in animals.

Non-genetic resistance facilitates survival while hindering the evolution of drug resistance due to intraspecific competition

Rising rates of resistance to antimicrobial drugs threaten the effective treatment of infections across the globe. Drug resistance has been established to emerge from non-genetic mechanisms as well as from genetic mechanisms. However, it is still unclear how non-genetic resistance affects the evolution of genetic drug resistance. We develop deterministic and stochastic population models that incorporate resource competition to quantitatively investigate the transition from non-genetic to genetic resistance during the exposure to static and cidal drugs. We find that non-genetic resistance facilitates the survival of cell populations during drug treatment while hindering the development of genetic resistance due to competition between the non-genetically and genetically resistant subpopulations. Non-genetic resistance in the presence of subpopulation competition increases the fixation times of drug resistance mutations, while increasing the probability of mutation before population extinction during cidal drug treatment. Intense intraspecific competition during drug treatment leads to extinction of susceptible and non-genetically resistant subpopulations. Alternating between drug and no drug conditions results in oscillatory population dynamics, increased resistance mutation fixation timescales, and reduced population survival. These findings advance our fundamental understanding of the evolution of resistance and may guide novel treatment strategies for patients with drug-resistant infections.

Low protein expression enhances phenotypic evolvability by intensifying selection on folding stability

Protein abundance affects the evolution of protein genotypes, but we do not know how it affects the evolution of protein phenotypes. Here we investigate the role of protein abundance in the evolvability of green fluorescent protein (GFP) towards the novel phenotype of cyan fluorescence. We evolve GFP in E. coli through multiple cycles of mutation and selection and show that low GFP expression facilitates the evolution of cyan fluorescence. A computational model whose predictions we test experimentally helps explain why: lowly expressed proteins are under stronger selection for proper folding, which facilitates their evolvability on short evolutionary time scales. The reason is that high fluorescence can be achieved by either few proteins that fold well or by many proteins that fold less well. In other words, we observe a synergy between a protein's scarcity and its stability. Because many proteins meet the essential requirements for this scarcity-stability synergy, it may be a widespread mechanism by which low expression helps proteins evolve new phenotypes and functions.

The ecological roles of bacterial chemotaxis

How bacterial chemotaxis is performed is much better understood than why. Traditionally, chemotaxis has been understood as a foraging strategy by which bacteria enhance their uptake of nutrients and energy, yet it has remained puzzling why certain less nutritious compounds are strong chemoattractants and vice versa. Recently, we have gained increased understanding of alternative ecological roles of chemotaxis, such as navigational guidance in colony expansion, localization of hosts or symbiotic partners and contribution to microbial diversity by the generation of spatial segregation in bacterial communities. Although bacterial chemotaxis has been observed in a wide range of environmental settings, insights into the phenomenon are mostly based on laboratory studies of model organisms. In this Review, we highlight how observing individual and collective migratory behaviour of bacteria in different settings informs the quantification of trade-offs, including between chemotaxis and growth. We argue that systematically mapping when and where bacteria are motile, in particular by transgenerational bacterial tracking in dynamic environments and in situ approaches from guts to oceans, will open the door to understanding the rich interplay between metabolism and growth and the contribution of chemotaxis to microbial life.

Narrow mutational signatures drive acquisition of multidrug resistance in the fungal pathogen Candida glabrata

Fungal infections are a growing medical concern, in part due to increased resistance to one or multiple antifungal drugs. However, the evolutionary processes underpinning the acquisition of antifungal drug resistance are poorly understood. Here, we used experimental microevolution to study the adaptation of the yeast pathogen Candida glabrata to fluconazole and anidulafungin, two widely used antifungal drugs with different modes of action. Our results show widespread ability of rapid adaptation to one or both drugs. Resistance, including multidrug resistance, is often acquired at moderate fitness costs and mediated by mutations in a limited set of genes that are recurrently and specifically mutated in strains adapted to each of the drugs. Importantly, we uncover a dual role of ERG3 mutations in resistance to anidulafungin and cross-resistance to fluconazole in a subset of anidulafungin-adapted strains. Our results shed light on the mutational paths leading to resistance and cross-resistance to antifungal drugs.

The Origins of Phenotypic Heterogeneity in Cancer

Heterogeneity is a pervasive feature of cancer, and understanding the sources and regulatory mechanisms underlying heterogeneity could provide key insights to help improve the diagnosis and treatment of cancer. In this review, we discuss the origin of heterogeneity in the phenotype of individual cancer cells. Genotype-phenotype (G-P) maps are widely used in evolutionary biology to represent the complex interactions of genes and the environment that lead to phenotypes that impact fitness. Here, we present the rationale of an extended G-P (eG-P) map with a cone structure in cancer. The eG-P cone is formed by cells that are similar at the genome layer but gradually increase variability in the epigenome, transcriptome, proteome, metabolome and signalome layers to produce large variability at the phenome layer. Experimental evidence from single-cell -omics analyses supporting the cancer eG-P cone concept is presented, and the impact of epimutations and the interaction of cancer and tumor microenvironmental eG-P cones are integrated with the current understanding of cancer biology. The eG-P cone concept uncovers potential therapeutic strategies to reduce cancer evolution and improve cancer treatment. More methods to study phenotypes in single cells will be key to better understand cancer cell fitness in tumor biology and therapeutics.

Bridging Tumorigenesis and Therapy Resistance With a Non-Darwinian and Non-Lamarckian Mechanism of Adaptive Evolution

Although neo-Darwinian (and less often Lamarckian) dynamics are regularly invoked to interpret cancer's multifarious molecular profiles, they shine little light on how tumorigenesis unfolds and often fail to fully capture the frequency and breadth of resistance mechanisms. This uncertainty frames one of the most problematic gaps between science and practice in modern times. Here, we offer a theory of adaptive cancer evolution, which builds on a molecular mechanism that lies outside neo-Darwinian and Lamarckian schemes. This mechanism coherently integrates non-genetic and genetic changes, ecological and evolutionary time scales, and shifts the spotlight away from positive selection towards purifying selection, genetic drift, and the creative-disruptive power of environmental change. The surprisingly simple use-it or lose-it rationale of the proposed theory can help predict molecular dynamics during tumorigenesis. It also provides simple rules of thumb that should help improve therapeutic approaches in cancer.

Imaging and Quantification of mRNA Molecules at Single-Cell Resolution in the Human Fungal Pathogen Candida albicans

The study of gene expression in fungi has typically relied on measuring transcripts in populations of cells. A major disadvantage of this approach is that the transcripts' spatial distribution and stochastic variation among individual cells within a clonal population is lost. Traditional fluorescence in situ hybridization techniques have been of limited use in fungi due to poor specificity and high background signal. Here, we report that in situ hybridization chain reaction (HCR), a method that employs split-initiator probes to trigger signal amplification upon mRNA-probe hybridization, is ideally suited for the imaging and quantification of low-abundance transcripts at single-cell resolution in the fungus Candida albicans. We show that HCR allows the absolute quantification of transcripts within a cell by microscopy as well as their relative quantification by flow cytometry. mRNA imaging also revealed the subcellular localization of specific transcripts. Furthermore, we establish that HCR is amenable to multiplexing by visualizing different transcripts in the same cell. Finally, we combine HCR with immunostaining to image specific mRNAs and proteins simultaneously within a single C. albicans cell. The fungus is a major pathogen in humans where it can colonize and invade mucosal surfaces and most internal organs. The technical development that we introduce, therefore, paves the way to study the patterns of expression of pathogenesis-associated C. albicans genes in infected organs at single-cell resolution. IMPORTANCE Tools to visualize and quantify transcripts at single-cell resolution have enabled the dissection of spatiotemporal patterns of gene expression in animal cells and tissues. Yet the accurate quantification of transcripts at single-cell resolution remains challenging for the much smaller microbial cells. Widespread phenomena such as stochastic variation in transcript levels among cells-even within a clonal population-seem to play important roles in the biology of many microorganisms. Investigating this process requires microbial cell-optimized procedures to image and measure mRNAs at single-molecule resolution. In this report, we adapt and expand in situ hybridization chain reaction (HCR) combined with split-initiator probes to visualize transcripts in the human-pathogenic fungus Candida albicans at high resolution.

Non-Genetic Diversity in Chemosensing and Chemotactic Behavior

Non-genetic phenotypic diversity plays a significant role in the chemotactic behavior of bacteria, influencing how populations sense and respond to chemical stimuli. First, we review the molecular mechanisms that generate phenotypic diversity in bacterial chemotaxis. Next, we discuss the functional consequences of phenotypic diversity for the chemosensing and chemotactic performance of single cells and populations. Finally, we discuss mechanisms that modulate the amount of phenotypic diversity in chemosensory parameters in response to changes in the environment.

Do Multiple Drug Resistance Transporters Interfere with Cell Functioning under Normal Conditions?

Eukaryotic cells rely on multiple mechanisms to protect themselves from exogenous toxic compounds. For instance, cells can limit penetration of toxic molecules through the plasma membrane or sequester them within the specialized compartments. Plasma membrane transporters with broad substrate specificity confer multiple drug resistance (MDR) to cells. These transporters efflux toxic compounds at the cost of ATP hydrolysis (ABC-transporters) or proton influx (MFS-transporters). In our review, we discuss the possible costs of having an active drug-efflux system using yeast cells as an example. The pleiotropic drug resistance (PDR) subfamily ABC-transporters are known to constitutively hydrolyze ATP even without any substrate stimulation or transport across the membrane. Besides, some MDR-transporters have flippase activity allowing transport of lipids from inner to outer lipid layer of the plasma membrane. Thus, excessive activity of MDR-transporters can adversely affect plasma membrane properties. Moreover, broad substrate specificity of ABC-transporters also suggests the possibility of unintentional efflux of some natural metabolic intermediates from the cells. Furthermore, in some microorganisms, transport of quorum-sensing factors is mediated by MDR transporters; thus, overexpression of the transporters can also disturb cell-to-cell communications. As a result, under normal conditions, cells keep MDR-transporter genes repressed and activate them only upon exposure to stresses. We speculate that exploiting limitations of the drug-efflux system is a promising strategy to counteract MDR in pathogenic fungi.

Does transcriptional heterogeneity facilitate the development of genetic drug resistance?

Non-genetic forms of antimicrobial (drug) resistance can result from cell-to-cell variability that is not encoded in the genetic material. Data from recent studies also suggest that non-genetic mechanisms can facilitate the development of genetic drug resistance. We speculate on how the interplay between non-genetic and genetic mechanisms may affect microbial adaptation and evolution during drug treatment. We argue that cellular heterogeneity arising from fluctuations in gene expression, epigenetic modifications, as well as genetic changes contribute to drug resistance at different timescales, and that the interplay between these mechanisms enhance pathogen resistance. Accordingly, developing a better understanding of the role of non-genetic mechanisms in drug resistance and how they interact with genetic mechanisms will enhance our ability to combat antimicrobial resistance. Also see the video abstract here: https://youtu.be/aefGpdh-bgU.

Optimizing mitochondrial maintenance in extended neuronal projections

Neurons rely on localized mitochondria to fulfill spatially heterogeneous metabolic demands. Mitochondrial aging occurs on timescales shorter than the neuronal lifespan, necessitating transport of fresh material from the soma. Maintaining an optimal distribution of healthy mitochondria requires an interplay between a stationary pool localized to sites of high metabolic demand and a motile pool capable of delivering new material. Interchange between these pools can occur via transient fusion / fission events or by halting and restarting entire mitochondria. Our quantitative model of neuronal mitostasis identifies key parameters that govern steady-state mitochondrial health at discrete locations. Very infrequent exchange between stationary and motile pools optimizes this system. Exchange via transient fusion allows for robust maintenance, which can be further improved by selective recycling through mitophagy. These results provide a framework for quantifying how perturbations in organelle transport and interactions affect mitochondrial homeostasis in neurons, a key aspect underlying many neurodegenerative disorders.

Tracking bacterial lineages in complex and dynamic environments with applications for growth control and persistence

As bacteria transition from exponential to stationary phase, they change substantially in size, morphology, growth and expression profiles. These responses also vary between individual cells, but it has proved difficult to track cell lineages along the growth curve to determine the progression of events or correlations between how individual cells enter and exit dormancy. Here, we developed a platform for tracking more than 10⁵ parallel cell lineages in dense and changing cultures, independently validating that the imaged cells closely track batch populations. Initial applications show that for both Escherichia coli and Bacillus subtilis, growth changes from an 'adder' mode in exponential phase to mixed 'adder-timers' entering stationary phase, and then a near-perfect 'sizer' upon exit-creating broadly distributed cell sizes in stationary phase but rapidly returning to narrowly distributed sizes upon exit. Furthermore, cells that undergo more divisions when entering stationary phase suffer reduced survival after long periods of dormancy but are the only cells observed that persist following antibiotic treatment.

How to characterize a strain? Clonal heterogeneity in industrial Saccharomyces influences both phenotypes and heterogeneity in phenotypes

Populations of microbes are constantly evolving heterogeneity that selection acts upon, yet heterogeneity is nontrivial to assess methodologically. The necessary practice of isolating single-cell colonies and thus subclone lineages for establishing, transferring, and using a strain results in single-cell bottlenecks with a generally neglected effect on the characteristics of the strain itself. Here, we present evidence that various subclone lineages for industrial yeasts sequenced for recent genomic studies show considerable differences, ranging from loss of heterozygosity to aneuploidies. Subsequently, we assessed whether phenotypic heterogeneity is also observable in industrial yeast, by individually testing subclone lineages obtained from products. Phenotyping of industrial yeast samples and their newly isolated subclones showed that single-cell bottlenecks during isolation can indeed considerably influence the observable phenotype. Next, we decoupled fitness distributions on the level of individual cells from clonal interference by plating single-cell colonies and quantifying colony area distributions. We describe and apply an approach using statistical modeling to compare the heterogeneity in phenotypes across samples and subclone lineages. One strain was further used to show how individual subclonal lineages are remarkably different not just in phenotype but also in the level of heterogeneity in phenotype. With these observations, we call attention to the fact that choosing an initial clonal lineage from an industrial yeast strain may vastly influence downstream performances and observations on karyotype, on phenotype, and also on heterogeneity.

Empirical evidence for epigenetic inheritance driving evolutionary adaptation

The cellular machinery that regulates gene expression can be self-propagated across cell division cycles and even generations. This renders gene expression states and their associated phenotypes heritable, independently of genetic changes. These phenotypic states, in turn, can be subject to selection and may influence evolutionary adaptation. In this review, we will discuss the molecular basis of epigenetic inheritance, the extent of its transmission and mechanisms of evolutionary adaptation. The current work shows that heritable gene expression can facilitate the process of adaptation through the increase of survival in a novel environment and by enlarging the size of beneficial mutational targets. Moreover, epigenetic control of gene expression enables stochastic switching between different phenotypes in populations that can potentially facilitate adaptation in rapidly fluctuating environments. Ecological studies of the variation of epigenetic markers (e.g. DNA methylation patterns) in wild populations show a potential contribution of this mode of inheritance to local adaptation in nature. However, the extent of the adaptive contribution of the naturally occurring variation in epi-alleles compared to genetic variation remains unclear. This article is part of the theme issue 'How does epigenetics influence the course of evolution?'

Topological signatures in regulatory network enable phenotypic heterogeneity in small cell lung cancer

Phenotypic (non-genetic) heterogeneity has significant implications for the development and evolution of organs, organisms, and populations. Recent observations in multiple cancers have unraveled the role of phenotypic heterogeneity in driving metastasis and therapy recalcitrance. However, the origins of such phenotypic heterogeneity are poorly understood in most cancers. Here, we investigate a regulatory network underlying phenotypic heterogeneity in small cell lung cancer, a devastating disease with no molecular targeted therapy. Discrete and continuous dynamical simulations of this network reveal its multistable behavior that can explain co-existence of four experimentally observed phenotypes. Analysis of the network topology uncovers that multistability emerges from two teams of players that mutually inhibit each other, but members of a team activate one another, forming a 'toggle switch' between the two teams. Deciphering these topological signatures in cancer-related regulatory networks can unravel their 'latent' design principles and offer a rational approach to characterize phenotypic heterogeneity in a tumor.

Microfluidic Evolution-On-A-Chip Reveals New Mutations that Cause Antibiotic Resistance

Microfluidic devices can mimic naturally occurring microenvironments and create microbial population heterogeneities ranging from planktonic cells to biofilm states. The exposure of such populations to spatially organized stress gradients can promote their adaptation into complex phenotypes, which are otherwise difficult to achieve with conventional experimental setups. Here a microfluidic chip that employs precise chemical gradients in consecutive microcompartments to perform microbial adaptive laboratory evolution (ALE), a key tool to study evolution in fundamental and applied contexts is described. In the chip developed here, microbial cells can be exposed to a defined profile of stressors such as antibiotics. By modulating this profile, stress adaptation in the chip through resistance or persistence can be specifically controlled. Importantly, chip-based ALE leads to the discovery of previously unknown mutations in Escherichia coli that confer resistance to nalidixic acid. The microfluidic device presented here can enhance the occurrence of mutations employing defined micro-environmental conditions to generate data to better understand the parameters that influence the mechanisms of antibiotic resistance.

Assessing the virulence of Cryptococcus neoformans causing meningitis in HIV infected and uninfected patients in Vietnam

We previously observed a substantial burden of cryptococcal meningitis in Vietnam atypically arising in individuals who are uninfected with human immunodeficiency virus (HIV). This disease was associated with a single genotype of Cryptococcus neoformans (sequence type [ST]5), which was significantly less common in HIV-infected individuals. Aiming to compare the phenotypic characteristics of ST5 and non-ST5 C. neoformans, we selected 30 representative Vietnamese isolates and compared their in vitro pathogenic potential and in vivo virulence. ST5 and non-ST5 organisms exhibited comparable characteristics with respect to in vitro virulence markers including melanin production, replication at 37°C, and growth in cerebrospinal fluid. However, the ST5 isolates had significantly increased variability in cellular and capsular sizing compared with non-ST5 organisms (P < .001). Counterintuitively, mice infected with ST5 isolates had significantly longer survival with lower fungal burdens at day 7 than non-ST5 isolates. Notably, ST5 isolates induced significantly greater initial inflammatory responses than non-ST5 strains, measured by TNF-α concentrations (P < .001). Despite being generally less virulent in the mouse model, we hypothesize that the significant within strain variation seen in ST5 isolates in the tested phenotypes may represent an evolutionary advantage enabling adaptation to novel niches including apparently immunocompetent human hosts.

Gene expression variability in human and chimpanzee populations share common determinants

Inter-individual variation in gene expression has been shown to be heritable and is often associated with differences in disease susceptibility between individuals. Many studies focused on mapping associations between genetic and gene regulatory variation, yet much less attention has been paid to the evolutionary processes that shape the observed differences in gene regulation between individuals in humans or any other primate. To begin addressing this gap, we performed a comparative analysis of gene expression variability and expression quantitative trait loci (eQTLs) in humans and chimpanzees, using gene expression data from primary heart samples. We found that expression variability in both species is often determined by non-genetic sources, such as cell-type heterogeneity. However, we also provide evidence that inter-individual variation in gene regulation can be genetically controlled, and that the degree of such variability is generally conserved in humans and chimpanzees. In particular, we found a significant overlap of orthologous genes associated with eQTLs in both species. We conclude that gene expression variability in humans and chimpanzees often evolves under similar evolutionary pressures.

The Boggarts of biology: how non-genetic changes influence the genotype

The notion that there is a one-one mapping from genotype to phenotype was overturned a long time ago. Along with genotype and environment, 'non-genetic changes' orchestrated by altered RNA and protein molecules also guide the development of phenotype. The idea that there is a route through which changes in phenotype can lead to changes in genotype impinges on several phenomena of molecular, developmental, evolutionary and applied interest. Phenotypic changes that do not alter the underlying DNA sequence have been studied across model systems (eg: DNA and histone modifications, RNA editing, prion formation) and are known to play an important role in short-term adaptation. However, because of their transient nature and unstable inheritance, the role of such changes in long-term evolution has remained controversial. I classify and review three ways in which non-genetic changes can influence genotype and impact cellular fitness across generations, with an emphasis on the enticing idea that they may act as stepping stones for genetic adaptation. I focus on work from microbial systems and attempt to highlight recent experiments and models that bear on this idea. Overall, I review evidence which suggests that non-genetic changes can impact phenotype via their influence on the genotype, and thus play a role in evolutionary change.

Gene expression noise can promote the fixation of beneficial mutations in fluctuating environments

Nongenetic phenotypic variation can either speed up or slow down adaptive evolution. We show that it can speed up evolution in environments where available carbon and energy sources change over time. To this end, we use an experimentally validated model of Escherichia coli growth on two alternative carbon sources, glucose and acetate. On the superior carbon source (glucose), all cells achieve high growth rates, while on the inferior carbon source (acetate) only a small fraction of the population manages to initiate growth. Consequently, populations experience a bottleneck when the environment changes from the superior to the inferior carbon source. Growth on the inferior carbon source depends on a circuit under the control of a transcription factor that is repressed in the presence of the superior carbon source. We show that noise in the expression of this transcription factor can increase the probability that cells start growing on the inferior carbon source. In doing so, it can decrease the severity of the bottleneck and increase mean population fitness whenever this fitness is low. A modest amount of noise can also enhance the fitness effects of a beneficial allele that increases the fraction of a population initiating growth on acetate. Additionally, noise can protect this allele from extinction, accelerate its spread, and increase its likelihood of going to fixation. Central to the adaptation-enhancing principle we identify is the ability of noise to mitigate population bottlenecks, particularly in environments that fluctuate periodically. Because such bottlenecks are frequent in fluctuating environments, and because periodically fluctuating environments themselves are common, this principle may apply to a broad range of environments and organisms.

Phenotypic noise and the cost of complexity

Experimental studies demonstrate the existence of phenotypic diversity despite constant genotype and environment. Theoretical models based on a single phenotypic character predict that during an adaptation event, phenotypic noise should be positively selected far from the fitness optimum because it increases the fitness of the genotype, and then be selected against when the population reaches the optimum. It is suggested that because of this fitness gain, phenotypic noise should promote adaptive evolution. However, it is unclear how the selective advantage of phenotypic noise is linked to the rate of evolution, and whether any advantage would hold for more realistic, multidimensional phenotypes. Indeed, complex organisms suffer a cost of complexity, where beneficial mutations become rarer as the number of phenotypic characters increases. Using a quantitative genetics approach, we first show that for a one-dimensional phenotype, phenotypic noise promotes adaptive evolution on plateaus of positive fitness, independently from the direct selective advantage on fitness. Second, we show that for multidimensional phenotypes, phenotypic noise evolves to a low-dimensional configuration, with elevated noise in the direction of the fitness optimum. Such a dimensionality reduction of the phenotypic noise promotes adaptive evolution and numerical simulations show that it reduces the cost of complexity.

Current Status and Applications of Adaptive Laboratory Evolution in Industrial Microorganisms

Adaptive laboratory evolution (ALE) is an evolutionary engineering approach in artificial conditions that improves organisms through the imitation of natural evolution. Due to the development of multi-level omics technologies in recent decades, ALE can be performed for various purposes at the laboratory level. This review delineates the basics of the experimental design of ALE based on several ALE studies of industrial microbial strains and updates current strategies combined with progressed metabolic engineering, in silico modeling and automation to maximize the evolution efficiency. Moreover, the review sheds light on the applicability of ALE as a strain development approach that complies with non-recombinant preferences in various food industries. Overall, recent progress in the utilization of ALE for strain development leading to successful industrialization is discussed.

For the Greater (Bacterial) Good: Heterogeneous Expression of Energetically Costly Virulence Factors

Bacterial populations are phenotypically heterogeneous, which allows subsets of cells to survive and thrive following changes in environmental conditions. For bacterial pathogens, changes within the host environment occur over the course of the immune response to infection and can result in exposure to host-derived, secreted antimicrobials or force direct interactions with immune cells. Many recent studies have shown host cell interactions promote virulence factor expression, forcing subsets of bacterial cells to battle the host response, while other bacteria reap the benefits of this pacification. It still remains unclear whether virulence factor expression is truly energetically costly within host tissues and whether expression is sufficient to impact the growth kinetics of virulence factor-expressing cells. However, it is clear that slow-growing subsets of bacteria emerge during infection and that these subsets are particularly difficult to eliminate with antibiotics. This minireview will focus on our current understanding of heterogenous virulence factor expression and discuss the evidence that supports or refutes the hypothesis that virulence factor expression is linked to slowed growth and antibiotic tolerance.

Compounding Effects of Climate Warming and Antibiotic Resistance

Bacteria have evolved diverse mechanisms to survive environments with antibiotics. Temperature is both a key factor that affects the survival of bacteria in the presence of antibiotics and an environmental trait that is drastically increasing due to climate change. Therefore, it is timely and important to understand links between temperature changes and selection of antibiotic resistance. This review examines these links by synthesizing results from laboratories, hospitals, and environmental studies. First, we describe the transient physiological responses to temperature that alter cellular behavior and lead to antibiotic tolerance and persistence. Second, we focus on the link between thermal stress and the evolution and maintenance of antibiotic resistance mutations. Finally, we explore how local and global changes in temperature are associated with increases in antibiotic resistance and its spread. We suggest that a multidisciplinary, multiscale approach is critical to fully understand how temperature changes are contributing to the antibiotic crisis.

Rapid evolution in plant-microbe interactions - a molecular genomics perspective

Rapid (co-)evolution at multiple timescales is a hallmark of plant-microbe interactions. The mechanistic basis for the rapid evolution largely rests on the features of the genomes of the interacting partners involved. Here, we review recent insights into genomic characteristics and mechanisms that enable rapid evolution of both plants and phytopathogens. These comprise fresh insights in allelic series of matching pairs of resistance and avirulence genes, the generation of novel pathogen effectors, the recently recognised small RNA warfare, and genomic aspects of secondary metabolite biosynthesis. In addition, we discuss the putative contributions of permissive host environments, transcriptional plasticity and the role of ploidy on the interactions. We conclude that the means underlying the rapid evolution of plant-microbe interactions are multifaceted and depend on the particular nature of each interaction.

Quantitation and Comparison of Phenotypic Heterogeneity Among Single Cells of Monoclonal Microbial Populations

Phenotypic heterogeneity within microbial populations arises even when the cells are exposed to putatively constant and homogeneous conditions. The outcome of this phenomenon can affect the whole function of the population, resulting in, for example, new "adapted" metabolic strategies and impacting its fitness at given environmental conditions. Accounting for phenotypic heterogeneity becomes thus necessary, due to its relevance in medical and applied microbiology as well as in environmental processes. Still, a comprehensive evaluation of this phenomenon requires a common and unique method of quantitation, which allows for the comparison between different studies carried out with different approaches. Consequently, in this study, two widely applicable indices for quantitation of heterogeneity were developed. The heterogeneity coefficient (HC) is valid when the population follows unimodal activity, while the differentiation tendency index (DTI) accounts for heterogeneity implying outbreak of subpopulations and multimodal activity. We demonstrated the applicability of HC and DTI for heterogeneity quantitation on stable isotope probing with nanoscale secondary ion mass spectrometry (SIP-nanoSIMS), flow cytometry, and optical microscopy datasets. The HC was found to provide a more accurate and precise measure of heterogeneity, being at the same time consistent with the coefficient of variation (CV) applied so far. The DTI is able to describe the differentiation in single-cell activity within monoclonal populations resolving subpopulations with low cell abundance, individual cells with similar phenotypic features (e.g., isotopic content close to natural abundance, as detected with nanoSIMS). The developed quantitation approach allows for a better understanding on the impact and the implications of phenotypic heterogeneity in environmental, medical and applied microbiology, microbial ecology, cell biology, and biotechnology.

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.