A Synthetic Oligo Library and Sequencing Approach Reveals an Insulation Mechanism Encoded within Bacterial σ 54 Promoters

Toolbox of Characterized Genetic Parts for Staphylococcus aureus

Staphylococcus aureus is an important clinical bacterium prevalent in human-associated microbiomes and the cause of many diseases. However, S. aureus has been intractable to synthetic biology approaches due to limited characterized genetic parts for this nonmodel Gram-positive bacterium. Moreover, genetic manipulation of S. aureus has relied on cumbersome and inefficient cloning strategies. Here, we report the first standardized genetic parts toolbox for S. aureus, which includes characterized promoters, ribosome binding sites, terminators, and plasmid replicons from a variety of bacteria for precise control of gene expression. We established a standard relative expression unit (REU) for S. aureus using a plasmid reference and characterized genetic parts in standardized REUs using S. aureus ATCC 12600. We constructed promoter and terminator part plasmids that are compatible with an efficient Type IIS DNA assembly strategy to effectively build multipart DNA constructs. A library of 24 constitutive promoters was built and characterized in S. aureus, which showed a 380-fold activity range. This promoter library was also assayed in Bacillus subtilis (122-fold activity range) to demonstrate the transferability of the constitutive promoters between these Gram-positive bacteria. By applying an iterative design-build-test-learn cycle, we demonstrated the use of our toolbox for the rational design and engineering of a tetracycline sensor in S. aureus using the PXyl-TetO aTc-inducible promoter that achieved 25.8-fold induction. This toolbox greatly expands the growing number of genetic parts for Gram-positive bacteria and will allow researchers to leverage synthetic biology approaches to study and engineer cellular processes in S. aureus.

AP profiling resolves co-translational folding pathway and chaperone interactions in vivo

Natural proteins have evolved to fold robustly along specific pathways. Folding begins during synthesis, guided by interactions of the nascent protein with the ribosome and molecular chaperones. However, the timing and progression of co-translational folding remain largely elusive, in part because the process is difficult to measure in the natural environment of the cytosol. We developed a high-throughput method to quantify co-translational folding in live cells that we term Arrest Peptide profiling (AP profiling). We employed AP profiling to delineate co-translational folding for a set of GTPase domains with very similar structures, defining how topology shapes folding pathways. Genetic ablation of major nascent chain-binding chaperones resulted in localized folding changes that suggest how functional redundancies among chaperones are achieved by distinct interactions with the nascent protein. Collectively, our studies provide a window into cellular folding pathways of complex proteins and pave the way for systematic studies on nascent protein folding at unprecedented resolution and throughput.

Challenges and advances towards the rational design of microalgal synthetic promoters in Chlamydomonas reinhardtii

Microalgae hold enormous potential to provide a safe and sustainable source of high-value compounds, acting as carbon-fixing biofactories that could help to mitigate rapidly progressing climate change. Bioengineering microalgal strains will be key to optimizing and modifying their metabolic outputs, and to render them competitive with established industrial biotechnology hosts, such as bacteria or yeast. To achieve this, precise and tuneable control over transgene expression will be essential, which would require the development and rational design of synthetic promoters as a key strategy. Among green microalgae, Chlamydomonas reinhardtii represents the reference species for bioengineering and synthetic biology; however, the repertoire of functional synthetic promoters for this species, and for microalgae generally, is limited in comparison to other commercial chassis, emphasizing the need to expand the current microalgal gene expression toolbox. Here, we discuss state-of-the-art promoter analyses, and highlight areas of research required to advance synthetic promoter development in C. reinhardtii. In particular, we exemplify high-throughput studies performed in other model systems that could be applicable to microalgae, and propose novel approaches to interrogating algal promoters. We lastly outline the major limitations hindering microalgal promoter development, while providing novel suggestions and perspectives for how to overcome them.

Molecular and experimental tools to design synthetic enhancers

Understanding the grammar of enhancers and how they regulate gene expression is key for both basic research and for the pharma and biotech industries. The design and characterization of synthetic enhancers can expand the known regulatory space. This is achieved by the utilization of DNA Oligo Libraries (OLs), which facilitates screening of as many as millions of synthetic enhancer variants simultaneously. This review includes the latest commercial DNA OL synthesis technology and its capabilities, and a general 'know-how' guide for the design, construction, and analysis of OL-based synthetic enhancer characterization experiments. Specifically, we focus on synthetic-enhancer-based massively parallel reporter assay, Sort-seq methodologies (e.g. flow cytometry, deep sequencing), and a brief description of machine learning-based attempts for OL-analysis and follow-up validation experiments.

Elucidation of Sequence-Function Relationships for an Improved Biobutanol In Vivo Biosensor in E. coli

Transcription factor (TF)–promoter pairs have been repurposed from native hosts to provide tools to measure intracellular biochemical production titer and dynamically control gene expression. Most often, native TF–promoter systems require rigorous screening to obtain desirable characteristics optimized for biotechnological applications. High-throughput techniques may provide a rational and less labor-intensive strategy to engineer user-defined TF–promoter pairs using fluorescence-activated cell sorting and deep sequencing methods (sort-seq). Based on the designed promoter library’s distribution characteristics, we elucidate sequence–function interactions between the TF and DNA. In this work, we use the sort-seq method to study the sequence–function relationship of a σ⁵⁴-dependent, butanol-responsive TF–promoter pair, BmoR-P_BMO derived from Thauera butanivorans, at the nucleotide level to improve biosensor characteristics, specifically an improved dynamic range. Activities of promoters from a mutagenized P_BMO library were sorted based on gfp expression and subsequently deep sequenced to correlate site-specific sequences with changes in dynamic range. We identified site-specific mutations that increase the sensor output. Double mutant and a single mutant, CA(129,130)TC and G(205)A, in P_BMO promoter increased dynamic ranges of 4-fold and 1.65-fold compared with the native system, respectively. In addition, sort-seq identified essential sites required for the proper function of the σ⁵⁴-dependent promoter biosensor in the context of the host. This work can enable high-throughput screening methods for strain development.

SpeedyGenesXL: an Automated, High-Throughput Platform for the Preparation of Bespoke Ultralarge Variant Libraries for Directed Evolution

Directed evolution of proteins is a highly effective strategy for tailoring biocatalysts to a particular application, and is capable of engineering improvements such as k_cat, thermostability and organic solvent tolerance. It is recognized that large and systematic libraries are required to navigate a protein's vast and rugged sequence landscape effectively, yet their preparation is nontrivial and commercial libraries are extremely costly. To address this, we have developed SpeedyGenesXL, an automated, high-throughput platform for the production of wild-type genes, Boolean OR, combinatorial, or combinatorial-OR-type libraries based on the SpeedyGenes methodology. Together this offers a flexible platform for library synthesis, capable of generating many different bespoke, diverse libraries simultaneously.

A broad analysis of splicing regulation in yeast using a large library of synthetic introns

RNA splicing is a key process in eukaryotic gene expression, in which an intron is spliced out of a pre-mRNA molecule to eventually produce a mature mRNA. Most intron-containing genes are constitutively spliced, hence efficient splicing of an intron is crucial for efficient regulation of gene expression. Here we use a large synthetic oligo library of ~20,000 variants to explore how different intronic sequence features affect splicing efficiency and mRNA expression levels in S. cerevisiae. Introns are defined by three functional sites, the 5’ donor site, the branch site, and the 3’ acceptor site. Using a combinatorial design of synthetic introns, we demonstrate how non-consensus splice site sequences in each of these sites affect splicing efficiency. We then show that S. cerevisiae splicing machinery tends to select alternative 3’ splice sites downstream of the original site, and we suggest that this tendency created a selective pressure, leading to the avoidance of cryptic splice site motifs near introns’ 3’ ends. We further use natural intronic sequences from other yeast species, whose splicing machineries have diverged to various extents, to show how intron architectures in the various species have been adapted to the organism’s splicing machinery. We suggest that the observed tendency for cryptic splicing is a result of a loss of a specific splicing factor, U2AF1. Lastly, we show that synthetic sequences containing two introns give rise to alternative RNA isoforms in S. cerevisiae, demonstrating that merely a synthetic fusion of two introns might be suffice to facilitate alternative splicing in yeast. Our study reveals novel mechanisms by which introns are shaped in evolution to allow cells to regulate their transcriptome. In addition, it provides a valuable resource to study the regulation of constitutive and alternative splicing in a model organism.

RNA splicing is a process in which parts of a new pre-mRNA are spliced out of the mRNA molecule to produce eventually a mature mRNA. Those RNA segments that are spliced out are termed introns, and they are found in most genes in eukaryotic organisms. Hence regulation of this process has a major role in the control of gene expression.

The budding yeast S. cerevisiae is a popular model organism for eukaryotic cell biology, but in terms of splicing it differs, as it has only few intron-containing genes. Nevertheless, this species has been used to study basic principles of splicing regulation based on its ~300 introns. Here we used the technology of a large synthetic genetic library to introduce many new intron-containing genes to the yeast genome, to explore splicing regulation at a wider scope than was possible so far.

Reassuringly, our results confirm known regulatory mechanisms, and further expand our understanding of splicing regulation, specifically how the yeast splicing machinery interacts with the end of introns, and how through evolution introns have evolved to avoid unwanted misidentifications of this end. We further demonstrate the potential of the yeast splicing machinery to alternatively splice a two-intron gene, which is common in other eukaryotes but rare in yeast. Our work presents a first-of-its-kind resource for the systematic study of splicing in live cells.

SOLQC: Synthetic Oligo Library Quality Control tool

MOTIVATION

Recent years have seen a growing number and an expanding scope of studies using synthetic oligo libraries for a range of applications in synthetic biology. As experiments are growing by numbers and complexity, analysis tools can facilitate quality control and support better assessment and inference.

RESULTS

We present a novel analysis tool, called SOLQC, which enables fast and comprehensive analysis of synthetic oligo libraries, based on NGS analysis performed by the user. SOLQC provides statistical information such as the distribution of variant representation, different error rates and their dependence on sequence or library properties. SOLQC produces graphical reports from the analysis, in a flexible format. We demonstrate SOLQC by analyzing literature libraries. We also discuss the potential benefits and relevance of the different components of the analysis.

AVAILABILITY AND IMPLEMENTATION

SOLQC is a free software for non-commercial use, available at https://app.gitbook.com/@yoav-orlev/s/solqc/. For commercial use please contact the authors.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Overcoming the design, build, test bottleneck for synthesis of nonrepetitive protein-RNA cassettes

We apply an oligo-library and machine learning-approach to characterize the sequence and structural determinants of binding of the phage coat proteins (CPs) of bacteriophages MS2 (MCP), PP7 (PCP), and Qβ (QCP) to RNA. Using the oligo library, we generate thousands of candidate binding sites for each CP, and screen for binding using a high-throughput dose-response Sort-seq assay (iSort-seq). We then apply a neural network to expand this space of binding sites, which allowed us to identify the critical structural and sequence features for binding of each CP. To verify our model and experimental findings, we design several non-repetitive binding site cassettes and validate their functionality in mammalian cells. We find that the binding of each CP to RNA is characterized by a unique space of sequence and structural determinants, thus providing a more complete description of CP-RNA interaction as compared with previous low-throughput findings. Finally, based on the binding spaces we demonstrate a computational tool for the successful design and rapid synthesis of functional non-repetitive binding-site cassettes.

Phage-coat proteins can be used to build synthetic biology parts. Here the authors use an oligo library and machine learning to predict and verify sequences based on binding scores.

Evolution of Regulated Transcription

The genomes of all organisms abound with various cis-regulatory elements, which control gene activity. Transcriptional enhancers are a key group of such elements in eukaryotes and are DNA regions that form physical contacts with gene promoters and precisely orchestrate gene expression programs. Here, we follow gradual evolution of this regulatory system and discuss its features in different organisms. In eubacteria, an enhancer-like element is often a single regulatory element, is usually proximal to the core promoter, and is occupied by one or a few activators. Activation of gene expression in archaea is accompanied by the recruitment of an activator to several enhancer-like sites in the upstream promoter region. In eukaryotes, activation of expression is accompanied by the recruitment of activators to multiple enhancers, which may be distant from the core promoter, and the activators act through coactivators. The role of the general DNA architecture in transcription control increases in evolution. As a whole, it can be seen that enhancers of multicellular eukaryotes evolved from the corresponding prototypic enhancer-like regulatory elements with the gradually increasing genome size of organisms.

Identification of a complete dibenzothiophene biodesulfurization operon and its regulator by functional metagenomics

Functional screening for aromatic ring oxygenases of an oil contaminated soil metagenome identified 25 different clones bearing monooxygenases coding genes. One fosmid bore an operon containing four tightly linked genes coding for a complete dibenzothiophene biodesulfurization pathway, which included the predicted monooxygenases DszC and DszA, the desulfinase DszB, and an FMN-oxidoreductase designated DszE. The dszEABC operon provided Escherichia coli with the ability to use dibenzothiophene as the only sulfur source. Transcription of the operon is driven from a σ^N -dependent promoter and regulated by an activator that was designated dszR. DszR has been purified and characterized in vitro and shown to be a constitutively active σ^N -dependent activator of the group IV, which binds to two contiguous sequences located upstream of the promoter. The dsz promoter and dszE and dszR genes have apparently been recruited from an aliphatic sulfonate biodegradation pathway. If transcribed from a heterologous upstream promoter, the σ^N -dependent promoter region functions as an 'insulator' that prevents translation of dszE, by binding with its ribosome binding site. Translational coupling, in turn, prevents translation of the downstream dszABC genes. The silencer combined with translational coupling thus represents an effective way of preventing expression of operons when spuriously transcribed from upstream promoters.

Data storage in DNA with fewer synthesis cycles using composite DNA letters

The density and long-term stability of DNA make it an appealing storage medium, particularly for long-term data archiving. Existing DNA storage technologies involve the synthesis and sequencing of multiple nominally identical molecules in parallel, resulting in information redundancy. We report the development of encoding and decoding methods that exploit this redundancy using composite DNA letters. A composite DNA letter is a representation of a position in a sequence that consists of a mixture of all four DNA nucleotides in a predetermined ratio. Our methods encode data using fewer synthesis cycles. We encode 6.4 MB into composite DNA, with distinguishable composition medians, using 20% fewer synthesis cycles per unit of data, as compared to previous reports. We also simulate encoding with larger composite alphabets, with distinguishable composition deciles, to show that 75% fewer synthesis cycles are potentially sufficient. We describe applicable error-correcting codes and inference methods, and investigate error patterns in the context of composite DNA letters.

Synthetic 5' UTRs Can Either Up- or Downregulate Expression upon RNA-Binding Protein Binding

The construction of complex gene-regulatory networks requires both inhibitory and upregulatory modules. However, the vast majority of RNA-based regulatory "parts" are inhibitory. Using a synthetic biology approach combined with SHAPE-seq, we explored the regulatory effect of RNA-binding protein (RBP)-RNA interactions in bacterial 5' UTRs. By positioning a library of RNA hairpins upstream of a reporter gene and co-expressing them with the matching RBP, we observed a set of regulatory responses, including translational stimulation, translational repression, and cooperative behavior. Our combined approach revealed three distinct states in vivo: in the absence of RBPs, the RNA molecules can be found in either a molten state that is amenable to translation or a structured phase that inhibits translation. In the presence of RBPs, the RNA molecules are in a semi-structured phase with partial translational capacity. Our work provides new insight into RBP-based regulation and a blueprint for designing complete gene-regulatory circuits at the post-transcriptional level.

Mapping DNA sequence to transcription factor binding energy in vivo

Despite the central importance of transcriptional regulation in biology, it has proven difficult to determine the regulatory mechanisms of individual genes, let alone entire gene networks. It is particularly difficult to decipher the biophysical mechanisms of transcriptional regulation in living cells and determine the energetic properties of binding sites for transcription factors and RNA polymerase. In this work, we present a strategy for dissecting transcriptional regulatory sequences using in vivo methods (massively parallel reporter assays) to formulate quantitative models that map a transcription factor binding site’s DNA sequence to transcription factor-DNA binding energy. We use these models to predict the binding energies of transcription factor binding sites to within 1 k_BT of their measured values. We further explore how such a sequence-energy mapping relates to the mechanisms of trancriptional regulation in various promoter contexts. Specifically, we show that our models can be used to design specific induction responses, analyze the effects of amino acid mutations on DNA sequence preference, and determine how regulatory context affects a transcription factor’s sequence specificity.

It has been said that we live in the “genomic era,” a time where we can readily sequence full genomes at will. However, it remains difficult to interpret much of the information within a genome. This is especially true of non-coding sequences such as promoters, which contain a number of features such as transcription factor binding sites that determine how genes are regulated. There is no straightforward regulatory “code” that tells us how transcription factor binding sites are organized within a promoter. In this work we examine how DNA sequence determines one of the most important features of a promoter, the strength with which a transcription factor binds to its DNA binding site. We discuss an approach to modeling DNA sequence-specific transcription factor binding energies in vivo using a massively parellel reporter assay. We develop models that allow us to predict the binding energy between a transcription factor and a mutated version of its binding site. We then show that this modeling technique can be used to address a number of scientific and design questions, such as engineering the behavior of genetic circuit elements or examining how transcription factors and their binding sites co-evolve.

An in Vivo Binding Assay for RNA-Binding Proteins Based on Repression of a Reporter Gene

We study translation repression in bacteria by engineering a regulatory circuit that functions as a binding assay for RNA binding proteins (RBP) in vivo. We do so by inducing expression of a fluorescent protein-RBP chimera, together with encoding its binding site at various positions within the ribosomal initiation region (+11-13 nt from the AUG) of a reporter module. We show that when bound by their cognate RBPs, the phage coat proteins for PP7 (PCP) and Qβ (QCP), strong repression is observed for all hairpin positions within the initiation region. Yet, a sharp transition to no-effect is observed when positioned in the elongation region, at a single-nucleotide resolution. Employing in vivo Selective 2'-hydroxyl acylation analyzed by primer extension followed by sequencing (SHAPE-seq) for a representative construct, established that in the translationally active state the mRNA molecule is nonstructured, while in the repressed state a structured signature was detected. We then utilize this regulatory phenomena to quantify the binding affinity of the coat proteins of phages MS2, PP7, GA, and Qβ to 14 cognate and noncognate binding sites in vivo. Using our circuit, we demonstrate qualitative differences between in vitro to in vivo binding characteristics for various variants when comparing to past studies. Furthermore, by introducing a simple mutation to the loop region for the Qβ-wt site, MCP binding is abolished, creating the first high-affinity QCP site that is completely orthogonal to MCP. Consequently, we demonstrate that our hybrid transcriptional-post-transcriptional circuit can be utilized as a binding assay to quantify RNA-RBP interactions in vivo.

Hidden Secrets of Sigma54 Promoters Revealed

Bacterial sigma54 (σ⁵⁴) promoters are the DNA-binding motif for σ⁵⁴-containing RNA polymerase (RNAP) holoenzymes. A recent study using a combination of synthetic oligonucleotide library screening, biochemical characterization, and bioinformatics has uncovered a new and unexpected role for σ⁵⁴ promoters, encoding a form of bacterial 'insulator sequence' to dampen unwanted translation.