Gene amplification as a form of population-level gene expression regulation

Population-level amplification of gene regulation by programmable gene transfer

Engineering cells to sense and respond to environmental cues often focuses on maximizing gene regulation at the single-cell level. Inspired by population-level control mechanisms like the immune response, we demonstrate dynamic control and amplification of gene regulation in bacterial populations using programmable plasmid-mediated gene transfer. By regulating plasmid loss rate, transfer rate and fitness effects via Cas9 endonuclease, F conjugation machinery and antibiotic selection, we modulate the fraction of plasmid-carrying cells, serving as an amplification factor for single-cell-level regulation. This approach expands the dynamic range of gene expression and allows orthogonal control across populations. Our platform offers a versatile strategy for dynamically regulating gene expression in engineered microbial communities.

Engineering plasmid copy number heterogeneity for dynamic microbial adaptation

Natural microbial populations exploit phenotypic heterogeneity for survival and adaptation. However, in engineering biology, limiting the sources of variability is a major focus. Here we show that intentionally coupling distinct plasmids via shared replication mechanisms enables bacterial populations to adapt to their environment. We demonstrate that plasmid coupling of carbon-metabolizing operons facilitates copy number tuning of an essential but burdensome construct through the action of a stably maintained, non-essential plasmid. For specific cost-benefit situations, incompatible two-plasmid systems can stably persist longer than compatible ones. We also show using microfluidics that plasmid coupling of synthetic constructs generates population-state memory of previous environmental adaptation without additional regulatory control. This work should help to improve the design of synthetic populations by enabling adaptive engineered strains to function under changing growth conditions without strict fine-tuning of the genetic circuitry.

Shining a light on Ohno's dilemma

Laboratory experiments on a fluorescent protein in E. coli reveal how duplicate genes are rapidly inactivated by mutations during evolution.

Large-scale duplication events underpin population-level flexibility in tRNA gene copy number in Pseudomonas fluorescens SBW25

The complement of tRNA genes within a genome is typically considered to be a (relatively) stable characteristic of an organism. Here, we demonstrate that bacterial tRNA gene set composition can be more flexible than previously appreciated, particularly regarding tRNA gene copy number. We report the high-rate occurrence of spontaneous, large-scale, tandem duplication events in laboratory populations of the bacterium Pseudomonas fluorescens SBW25. The identified duplications are up to ∼1 Mb in size (∼15% of the wildtype genome) and are predicted to change the copy number of up to 917 genes, including several tRNA genes. The observed duplications are inherently unstable: they occur, and are subsequently lost, at extremely high rates. We propose that this unusually plastic type of mutation provides a mechanism by which tRNA gene set diversity can be rapidly generated, while simultaneously preserving the underlying tRNA gene set in the absence of continued selection. That is, if a tRNA set variant provides no fitness advantage, then high-rate segregation of the duplication ensures the maintenance of the original tRNA gene set. However, if a tRNA gene set variant is beneficial, the underlying duplication fragment(s) may persist for longer and provide raw material for further, more stable, evolutionary change.

Duplicated antibiotic resistance genes reveal ongoing selection and horizontal gene transfer in bacteria

Horizontal gene transfer (HGT) and gene duplication are often considered as separate mechanisms driving the evolution of new functions. However, the mobile genetic elements (MGEs) implicated in HGT can copy themselves, so positive selection on MGEs could drive gene duplications. Here, we use a combination of modeling and experimental evolution to examine this hypothesis and use long-read genome sequences of tens of thousands of bacterial isolates to examine its generality in nature. Modeling and experiments show that antibiotic selection can drive the evolution of duplicated antibiotic resistance genes (ARGs) through MGE transposition. A key implication is that duplicated ARGs should be enriched in environments associated with antibiotic use. To test this, we examined the distribution of duplicated ARGs in 18,938 complete bacterial genomes with ecological metadata. Duplicated ARGs are highly enriched in bacteria isolated from humans and livestock. Duplicated ARGs are further enriched in an independent set of 321 antibiotic-resistant clinical isolates. Our findings indicate that duplicated genes often encode functions undergoing positive selection and horizontal gene transfer in microbial communities.

An overview of the co-transcription factor NACC1: Beyond its pro-tumor effects

Nucleus accumbens-associated protein 1 (NACC1) is a member of the broad complex, tramtrack, bric-a-brac/poxvirus and zinc finger (BTB/POZ) protein families, mainly exerting its biological functions as a transcription co-regulator. NACC1 forms homo- or hetero-dimers through the BTB/POZ or BANP, E5R, and NACC1 (BEN) domain with other transcriptional regulators to regulate downstream signals. Recently, the overexpression of NACC1 has been observed in various tumors and is positively associated with tumor progression, high recurrence rate, indicating poor prognosis. NACC1 also regulates biological processes such as embryonic development, stem cell pluripotency, innate immunity, and related diseases. Our review combines recent research to summarize advancements in the structure, biological functions, and relative molecular mechanisms of NACC1. The future development of NACC1 clinical appliances is also discussed.

The genome formula of a multipartite virus is regulated both at the individual segment and the segment group levels

Differential accumulation of the distinct genome segments is a common feature of viruses with segmented genomes. The reproducible and specific pattern of genome segment accumulation within the host is referred to as the "genome formula". There is speculation and some experimental support for a functional role of the genome formula by modulating gene expression through copy number variations. However, the mechanisms of genome formula regulation have not yet been identified. In this study, we investigated whether the genome formula of the octopartite nanovirus faba bean necrotic stunt virus (FBNSV) is regulated by processes acting at the individual segment vs. viral population levels. We used a leaf infiltration system to show that the two most accumulated genome segments of the FBNSV possess a greater intrinsic accumulation capacity in Vicia faba tissues than the other segments. Nevertheless, processes acting at the individual segment level are insufficient to generate the genome formula, suggesting the involvement of additional mechanisms acting at the supra-segment level. Indeed, the absence of segments with important functions during systemic infection strongly modifies the relative frequency of the others, indicating that the genome formula is a property of the segment group. Together, these results demonstrate that the FBNSV genome formula is shaped by a complex process acting at both the individual segment and the segment group levels.

Replication mechanisms of circular ssDNA plant viruses and their potential implication in viral gene expression regulation

The replication of members of the two circular single-stranded DNA (ssDNA) virus families Geminiviridae and Nanoviridae, the only ssDNA viruses infecting plants, is believed to be processed by rolling-circle replication (RCR) and recombination-dependent replication (RDR) mechanisms. RCR is a ubiquitous replication mode for circular ssDNA viruses and involves a virus-encoded Replication-associated protein (Rep) which fulfills multiple functions in the replication mechanism. Two key genomic elements have been identified for RCR in Geminiviridae and Nanoviridae: (i) short iterative sequences called iterons which determine the specific recognition of the viral DNA by the Rep and (ii) a sequence enabling the formation of a stem-loop structure which contains a conserved motif and constitutes the origin of replication. In addition, studies in Geminiviridae provided evidence for a second replication mode, RDR, which has also been documented in some double-stranded DNA viruses. Here, we provide a synthesis of the current understanding of the two presumed replication modes of Geminiviridae and Nanoviridae, and we identify knowledge gaps and discuss the possibility that these replication mechanisms could regulate viral gene expression through modulation of gene copy number.

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.

Gene amplifications cause high-level resistance against albicidin in gram-negative bacteria

Antibiotic resistance is a continuously increasing concern for public healthcare. Understanding resistance mechanisms and their emergence is crucial for the development of new antibiotics and their effective use. The peptide antibiotic albicidin is such a promising candidate that, as a gyrase poison, shows bactericidal activity against a wide range of gram-positive and gram-negative bacteria. Here, we report the discovery of a gene amplification-based mechanism that imparts an up to 1000-fold increase in resistance levels against albicidin. RNA sequencing and proteomics data show that this novel mechanism protects Salmonella Typhimurium and Escherichia coli by increasing the copy number of STM3175 (YgiV), a transcription regulator with a GyrI-like small molecule binding domain that traps albicidin with high affinity. X-ray crystallography and molecular docking reveal a new conserved motif in the binding groove of the GyrI-like domain that can interact with aromatic building blocks of albicidin. Phylogenetic studies suggest that this resistance mechanism is ubiquitous in gram-negative bacteria, and our experiments confirm that STM3175 homologs can confer resistance in pathogens such as Vibrio vulnificus and Pseudomonas aeruginosa.

Microfluidics for long-term single-cell time-lapse microscopy: Advances and applications

Cells are inherently dynamic, whether they are responding to environmental conditions or simply at equilibrium, with biomolecules constantly being made and destroyed. Due to their small volumes, the chemical reactions inside cells are stochastic, such that genetically identical cells display heterogeneous behaviors and gene expression profiles. Studying these dynamic processes is challenging, but the development of microfluidic methods enabling the tracking of individual prokaryotic cells with microscopy over long time periods under controlled growth conditions has led to many discoveries. This review focuses on the recent developments of one such microfluidic device nicknamed the mother machine. We overview the original device design, experimental setup, and challenges associated with this platform. We then describe recent methods for analyzing experiments using automated image segmentation and tracking. We further discuss modifications to the experimental setup that allow for time-varying environmental control, replicating batch culture conditions, cell screening based on their dynamic behaviors, and to accommodate a variety of microbial species. Finally, this review highlights the discoveries enabled by this technology in diverse fields, such as cell-size control, genetic mutations, cellular aging, and synthetic biology.

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.

Positively selected genes in the hoary bat (Lasiurus cinereus) lineage: prominence of thymus expression, immune and metabolic function, and regions of ancient synteny

Background

Bats of the genus Lasiurus occur throughout the Americas and have diversified into at least 20 species among three subgenera. The hoary bat (Lasiurus cinereus) is highly migratory and ranges farther across North America than any other wild mammal. Despite the ecological importance of this species as a major insect predator, and the particular susceptibility of lasiurine bats to wind turbine strikes, our understanding of hoary bat ecology, physiology, and behavior remains poor.

Methods

To better understand adaptive evolution in this lineage, we used whole-genome sequencing to identify protein-coding sequence and explore signatures of positive selection. Gene models were predicted with Maker and compared to seven well-annotated and phylogenetically representative species. Evolutionary rate analysis was performed with PAML.

Results

Of 9,447 single-copy orthologous groups that met evaluation criteria, 150 genes had a significant excess of nonsynonymous substitutions along the L. cinereus branch (P < 0.001 after manual review of alignments). Selected genes as a group had biased expression, most strongly in thymus tissue. We identified 23 selected genes with reported immune functions as well as a divergent paralog of Steep1 within suborder Yangochiroptera. Seventeen genes had roles in lipid and glucose metabolic pathways, partially overlapping with 15 mitochondrion-associated genes; these adaptations may reflect the metabolic challenges of hibernation, long-distance migration, and seasonal variation in prey abundance. The genomic distribution of positively selected genes differed significantly from background expectation by discrete Kolmogorov-Smirnov test (P < 0.001). Remarkably, the top three physical clusters all coincided with islands of conserved synteny predating Mammalia, the largest of which shares synteny with the human cat-eye critical region (CECR) on 22q11. This observation coupled with the expansion of a novel Tbx1-like gene family may indicate evolutionary innovation during pharyngeal arch development: both the CECR and Tbx1 cause dosage-dependent congenital abnormalities in thymus, heart, and head, and craniodysmorphy is associated with human orthologs of other positively selected genes as well.

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.

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.

A Bacillus subtilis ΔpdxT mutant suppresses vitamin B6 limitation by acquiring mutations enhancing pdxS gene dosage and ammonium assimilation

Pyridoxal-5'-phosphate (PLP), the biologically active form of vitamin B6, serves as a cofactor for many enzymes. The Gram-positive model bacterium Bacillus subtilis synthesizes PLP via the PdxST enzyme complex, consisting of the PdxT glutaminase and the PdxS PLP synthase subunits, respectively. PdxT converts glutamine to glutamate and ammonia of which the latter is channelled to PdxS. At high extracellular ammonium concentrations, the PdxS PLP synthase subunit does not depend on PdxT. Here, we assessed the potential of a B. subtilis ΔpdxT mutant to adapt to PLP limitation at the genome level. The majority of ΔpdxT suppressors had amplified a genomic region containing the pdxS gene. We also identified mutants having acquired as yet undescribed mutations in ammonium assimilation genes, indicating that the overproduction of PdxS and the NrgA ammonium transporter partially relieve vitamin B6 limitation in a ΔpdxT mutant when extracellular ammonium is scarce. Furthermore, we found that PdxS positively affects complex colony formation in B. subtilis. The catalytic mechanism of the PdxS PLP synthase subunit could be the reason for the limited evolution of the enzyme and why we could not identify a PdxS variant producing PLP independently of PdxT at low ammonium concentrations.

Local genetic context shapes the function of a gene regulatory network

Gene expression levels are influenced by multiple coexisting molecular mechanisms. Some of these interactions such as those of transcription factors and promoters have been studied extensively. However, predicting phenotypes of gene regulatory networks (GRNs) remains a major challenge. Here, we use a well-defined synthetic GRN to study in Escherichia coli how network phenotypes depend on local genetic context, i.e. the genetic neighborhood of a transcription factor and its relative position. We show that one GRN with fixed topology can display not only quantitatively but also qualitatively different phenotypes, depending solely on the local genetic context of its components. Transcriptional read-through is the main molecular mechanism that places one transcriptional unit (TU) within two separate regulons without the need for complex regulatory sequences. We propose that relative order of individual TUs, with its potential for combinatorial complexity, plays an important role in shaping phenotypes of GRNs.

The Antibiotic Dosage of Fastest Resistance Evolution: gene amplifications underpinning the inverted-U

To determine the dosage at which antibiotic resistance evolution is most rapid, we treated Escherichia coli in vitro, deploying the antibiotic erythromycin at dosages ranging from zero to high. Adaptation was fastest just below erythromycin’s minimal inhibitory concentration (MIC) and genotype-phenotype correlations determined from whole genome sequencing revealed the molecular basis: simultaneous selection for copy number variation in three resistance mechanisms which exhibited an “inverted-U” pattern of dose-dependence, as did several insertion sequences and an integron. Many genes did not conform to this pattern, however, reflecting changes in selection as dose increased: putative media adaptation polymorphisms at zero antibiotic dosage gave way to drug target (ribosomal RNA operon) amplification at mid dosages whereas prophage-mediated drug efflux amplifications dominated at the highest dosages. All treatments exhibited E. coli increases in the copy number of efflux operons acrAB and emrE at rates that correlated with increases in population density. For strains where the inverted-U was no longer observed following the genetic manipulation of acrAB, it could be recovered by prolonging the antibiotic treatment at subMIC dosages.

Beyond horizontal gene transfer: the role of plasmids in bacterial evolution

Plasmids have a key role in bacterial ecology and evolution because they mobilize accessory genes by horizontal gene transfer. However, recent studies have revealed that the evolutionary impact of plasmids goes above and beyond their being mere gene delivery platforms. Plasmids are usually kept at multiple copies per cell, producing islands of polyploidy in the bacterial genome. As a consequence, the evolution of plasmid-encoded genes is governed by a set of rules different from those affecting chromosomal genes, and these rules are shaped by unusual concepts in bacterial genetics, such as genetic dominance, heteroplasmy or segregational drift. In this Review, we discuss recent advances that underscore the importance of plasmids in bacterial ecology and evolution beyond horizontal gene transfer. We focus on new evidence that suggests that plasmids might accelerate bacterial evolution, mainly by promoting the evolution of plasmid-encoded genes, but also by enhancing the adaptation of their host chromosome. Finally, we integrate the most relevant theoretical and empirical studies providing a global understanding of the forces that govern plasmid-mediated evolution in bacteria.

Diversity of Omega Glutathione Transferases in mushroom-forming fungi revealed by phylogenetic, transcriptomic, biochemical and structural approaches

The Omega class of glutathione transferases (GSTs) forms a distinct class within the cytosolic GST superfamily because most of them possess a catalytic cysteine residue. The human GST Omega 1 isoform was first characterized twenty years ago, but it took years of work to clarify the roles of the human isoforms. Concerning the kingdom of fungi, little is known about the cellular functions of Omega glutathione transferases (GSTOs), although they are widely represented in some of these organisms. In this study, we re-assess the phylogeny and the classification of GSTOs based on 240 genomes of mushroom-forming fungi (Agaricomycetes). We observe that the number of GSTOs is not only extended in the order of Polyporales but also in other orders such as Boletales. Our analysis leads to a new classification in which the fungal GSTOs are divided into two Types A and B. The catalytic residue of Type-A is either cysteine or serine, while that of Type-B is cysteine. The present study focuses on Trametes versicolor GSTO isoforms that possess a catalytic cysteine residue. Transcriptomic data show that Type-A GSTOs are constitutive enzymes while Type-B are inducible ones. The crystallographic analysis reveals substantial structural differences between the two types while they have similar biochemical profiles in the tested conditions. Additionally, these enzymes have the ability to bind antioxidant molecules such as wood polyphenols in two possible binding sites as observed from X-ray structures. The multiplication of GSTOs could allow fungal organisms to adapt more easily to new environments.