Manipulating nucleosome disfavoring sequences allows fine-tune regulation of gene expression in yeast

Living in a noisy world-origins of gene expression noise and its impact on cellular decision-making

The expression level of a gene can vary between genetically identical cells under the same environmental condition-a phenomenon referred to as gene expression noise. Several studies have now elucidated a central role of transcription factors in the generation of expression noise. Transcription factors, as the key components of gene regulatory networks, drive many important cellular decisions in response to cellular and environmental signals. Therefore, a very relevant question is how expression noise impacts gene regulation and influences cellular decision-making. In this Review, we summarize the current understanding of the molecular origins of expression noise, highlighting the role of transcription factors in this process, and discuss the ways in which noise can influence cellular decision-making. As advances in single-cell technologies open new avenues for studying expression noise as well as gene regulatory circuits, a better understanding of the influence of noise on cellular decisions will have important implications for many biological processes.

Integrated strategies for efficient production of Streptomyces mobaraensis transglutaminase in Komagataella phaffii

Transglutaminase (TGase) from Streptomyces mobaraensis commonly used to improve protein-based foods due to its unique enzymatic reactions, which imply considerable attention in its production. Recently, TGase exhibit broad market potential in non-food industries. However, achieving efficient synthesis of TGase remains a significant challenge. Herein, we achieved a substantial amount of a fully functional and kinetically stable TGase produced by Komagataella phaffii (Pichia pastoris) using multiple strategies including Geneticin (G418) screening, combinatorial mutations, promoter optimization, and co-expression. The active TGase expression reached a maximum of 10.1 U mL-1 in shake flask upon 96 h of induction, which was 3.8-fold of the wild type. Also, the engineered strain exhibited a 6.4-fold increase in half-life and a 2-fold increase in specific activity, reaching 172.67 min at 60 °C (t1/2(60 °C)) and 65.3 U mg-1, respectively. Moreover, the high-cell density cultivation in 5-L fermenter was also applied to test the productivity at large scale. Following optimization at a fermenter, the secretory yield of TGase reached 47.96 U mL-1 in the culture supernatant. Given the complexity inherent in protein expression and secretion, our research is of great significance and offers a comprehensive guide for improving the production of a wide range of heterologous proteins.

MutLα suppresses error-prone DNA mismatch repair and preferentially protects noncoding DNA from mutations

The DNA mismatch repair (MMR) system promotes genome stability and protects humans from certain types of cancer. Its primary function is the correction of DNA polymerase errors. MutLα is an important eukaryotic MMR factor. We have examined the contributions of MutLα to maintaining genome stability. We show here that loss of MutLα in yeast increases the genome-wide mutation rate by ∼130-fold and generates a genome-wide mutation spectrum that consists of small indels and base substitutions. We also show that loss of yeast MutLα leads to error-prone MMR that produces T > C base substitutions in 5'-ATA-3' sequences. In agreement with this finding, our examination of human whole-genome DNA sequencing data has revealed that loss of MutLα in induced pluripotent stem cells triggers error-prone MMR that leads to the formation of T > C mutations in 5'-NTN-3' sequences. Our further analysis has shown that MutLα-independent MMR plays a role in suppressing base substitutions in N3 homopolymeric runs. In addition, we describe that MutLα preferentially protects noncoding DNA from mutations. Our study defines the contributions of MutLα-dependent and independent mechanisms to genome-wide MMR.

The scales, mechanisms, and dynamics of the genome architecture

Even when split into several chromosomes, DNA molecules that make up our genome are too long to fit into the cell nuclei unless massively folded. Such folding must accommodate the need for timely access to selected parts of the genome by transcription factors, RNA polymerases, and DNA replication machinery. Here, we review our current understanding of the genome folding inside the interphase nuclei. We consider the resulting genome architecture at three scales with a particular focus on the intermediate (meso) scale and summarize the insights gained from recent experimental observations and diverse computational models.

MutLα suppresses error-prone DNA mismatch repair and preferentially protects noncoding DNA from mutations

The DNA mismatch repair (MMR) system promotes genome stability and protects humans from certain types of cancer. Its primary function is the correction of DNA polymerase errors. MutLα is an important eukaryotic MMR factor. We have examined the contributions of MutLα to maintaining genome stability. We show here that loss of MutLα in yeast increases the genome-wide mutation rate by ~130-fold and generates a genome-wide mutation spectrum that consists of small indels and base substitutions. We also show that loss of yeast MutLα leads to error-prone MMR that produces T>C base substitutions in 5'-ATA-3' sequences. In agreement with this finding, our examination of human whole genome DNA sequencing data has revealed that loss of MutLα in induced pluripotent stem cells triggers error-prone MMR that leads to the formation of T>C mutations in 5'-NTN-3' sequences. Our further analysis has shown that MutLα-independent MMR plays a role in suppressing base substitutions in N3 homopolymeric runs. In addition, we describe that MutLα preferentially defends noncoding DNA from mutations. Our study defines the contributions of MutLα-dependent and independent mechanisms to genome-wide MMR.

Specific Patterns in Correlations of Super-Short Tandem Repeats (SSTRs) with G+C Content, Genic and Intergenic Regions, and Retrotransposons on All Human Chromosomes

The specific characteristics of k-mer words (2 ≤ k ≤ 11) regarding genomic distribution and evolutionary conservation were recently found. Among them are, in high abundance, words with a tandem repeat structure (repeat unit length of 1 bp to 3 bp). Furthermore, there seems to be a class of extremely short tandem repeats (≤12 bp), so far overlooked, that are non-random-distributed and, therefore, may play a crucial role in the functioning of the genome. In the following article, the positional distributions of these motifs we call super-short tandem repeats (SSTRs) were compared to other functional elements, like genes and retrotransposons. We found length- and sequence-dependent correlations between the local SSTR density and G+C content, and also between the density of SSTRs and genes, as well as correlations with retrotransposon density. In addition to many general interesting relations, we found that SINE Alu has a strong influence on the local SSTR density. Moreover, the observed connection of SSTR patterns to pseudogenes and -exons might imply a special role of SSTRs in gene expression. In summary, our findings support the idea of a special role and the functional relevance of SSTRs in the genome.

Systematic sequence engineering enhances the induction strength of the glucose-regulated GTH1 promoter of Komagataella phaffii

The promoter of the high-affinity glucose transporter Gth1 (PGTH1) is tightly repressed on glucose and glycerol surplus, and strongly induced in glucose-limitation, thus enabling regulated methanol-free production processes in the yeast production host Komagataella phaffii. To further improve this promoter, an intertwined approach of nucleotide diversification through random and rational engineering was pursued. Random mutagenesis and fluorescence activated cell sorting of PGTH1 yielded five variants with enhanced induction strength. Reverse engineering of individual point mutations found in the improved variants identified two single point mutations with synergistic action. Sequential deletions revealed the key promoter segments for induction and repression properties, respectively. Combination of the single point mutations and the amplification of key promoter segments led to a library of novel promoter variants with up to 3-fold higher activity. Unexpectedly, the effect of gaining or losing a certain transcription factor binding site (TFBS) was highly dependent on its context within the promoter. Finally, the applicability of the novel promoter variants for biotechnological production was proven for the secretion of different recombinant model proteins in fed batch cultivation, where they clearly outperformed their ancestors. In addition to advancing the toolbox for recombinant protein production and metabolic engineering of K. phaffii, we discovered single nucleotide positions and correspondingly affected TFBS that distinguish between glycerol- and glucose-mediated repression of the native promoter.

Short tandem repeats bind transcription factors to tune eukaryotic gene expression

Short tandem repeats (STRs) are enriched in eukaryotic cis-regulatory elements and alter gene expression, yet how they regulate transcription remains unknown. We found that STRs modulate transcription factor (TF)-DNA affinities and apparent on-rates by about 70-fold by directly binding TF DNA-binding domains, with energetic impacts exceeding many consensus motif mutations. STRs maximize the number of weakly preferred microstates near target sites, thereby increasing TF density, with impacts well predicted by statistical mechanics. Confirming that STRs also affect TF binding in cells, neural networks trained only on in vivo occupancies predicted effects identical to those observed in vitro. Approximately 90% of TFs preferentially bound STRs that need not resemble known motifs, providing a cis-regulatory mechanism to target TFs to genomic sites.

Using Synthetic DNA Libraries to Investigate Chromatin and Gene Regulation

Despite the recent explosion in genome-wide studies in chromatin and gene regulation, we are still far from extracting a set of genetic rules that can predict the function of the regulatory genome. One major reason for this deficiency is that gene regulation is a multi-layered process that involves an enormous variable space, which cannot be fully explored using native genomes. This problem can be partially solved by introducing synthetic DNA libraries into cells, a method that can test the regulatory roles of thousands to millions of sequences with limited variables. Here, we review recent applications of this method to study transcription factor (TF) binding, nucleosome positioning, and transcriptional activity. We discuss the design principles, experimental procedures, and major findings from these studies and compare the pros and cons of different approaches.

oFlowSeq: a quantitative approach to identify protein coding mutations affecting cell type enrichment using mosaic CRISPR-Cas9 edited cerebral organoids

Cerebral organoids are comprised of diverse cell types found in the developing human brain, and can be leveraged in the identification of critical cell types perturbed by genetic risk variants in common, neuropsychiatric disorders. There is great interest in developing high-throughput technologies to associate genetic variants with cell types. Here, we describe a high-throughput, quantitative approach (oFlowSeq) by utilizing CRISPR-Cas9, FACS sorting, and next-generation sequencing. Using oFlowSeq, we found that deleterious mutations in autism-associated gene KCTD13 resulted in increased proportions of Nestin+ cells and decreased proportions of TRA-1-60+ cells within mosaic cerebral organoids. We further identified that a locus-wide CRISPR-Cas9 survey of another 18 genes in the 16p11.2 locus resulted in most genes with > 2% maximum editing efficiencies for short and long indels, suggesting a high feasibility for an unbiased, locus-wide experiment using oFlowSeq. Our approach presents a novel method to identify genotype-to-cell type imbalances in an unbiased, high-throughput, quantitative manner.

Moderation of Structural DNA Properties by Coupled Dinucleotide Contents in Eukaryotes

Dinucleotides are known as determinants for various structural and physiochemical properties of DNA and for binding affinities of proteins to DNA. These properties (e.g., stiffness) and bound proteins (e.g., transcription factors) are known to influence important biological functions, such as transcription regulation and 3D chromatin organization. Accordingly, the question arises of how the considerable variations in dinucleotide contents of eukaryotic chromosomes could still provide consistent DNA properties resulting in similar functions and 3D conformations. In this work, we investigate the hypothesis that coupled dinucleotide contents influence DNA properties in opposite directions to moderate each other's influences. Analyzing all 2478 chromosomes of 155 eukaryotic species, considering bias from coding sequences and enhancers, we found sets of correlated and anti-correlated dinucleotide contents. Using computational models, we estimated changes of DNA properties resulting from this coupling. We found that especially pure A/T dinucleotides (AA, TT, AT, TA), known to influence histone positioning and AC/GT contents, are relevant moderators and that, e.g., the Roll property, which is known to influence histone affinity of DNA, is preferably moderated. We conclude that dinucleotide contents might indirectly influence transcription and chromatin 3D conformation, via regulation of histone occupancy and/or other mechanisms.

Mutation and selection processes regulating short tandem repeats give rise to genetic and phenotypic diversity across species

Short tandem repeats (STRs) are units of 1-6 bp that repeat in a tandem fashion in DNA. Along with single nucleotide polymorphisms and large structural variations, they are among the major genomic variants underlying genetic, and likely phenotypic, divergence. STRs experience mutation rates that are orders of magnitude higher than other well-studied genotypic variants. Frequent copy number changes result in a wide range of alleles, and provide unique opportunities for modulating complex phenotypes through variation in repeat length. While classical studies have identified key roles of individual STR loci, the advent of improved sequencing technology, high-quality genome assemblies for diverse species, and bioinformatics methods for genome-wide STR analysis now enable more systematic study of STR variation across wide evolutionary ranges. In this review, we explore mutation and selection processes that affect STR copy number evolution, and how these processes give rise to varying STR patterns both within and across species. Finally, we review recent examples of functional and adaptive changes linked to STRs.

Micro and macroevolution of sea anemone venom phenotype

Venom is a complex trait with substantial inter- and intraspecific variability resulting from strong selective pressures acting on the expression of many toxic proteins. However, understanding the processes underlying toxin expression dynamics that determine the venom phenotype remains unresolved. By interspecific comparisons we reveal that toxin expression in sea anemones evolves rapidly and that in each species different toxin family dictates the venom phenotype by massive gene duplication events. In-depth analysis of the sea anemone, Nematostella vectensis, revealed striking variation of the dominant toxin (Nv1) diploid copy number across populations (1-24 copies) resulting from independent expansion/contraction events, which generate distinct haplotypes. Nv1 copy number correlates with expression at both the transcript and protein levels with one population having a near-complete loss of Nv1 production. Finally, we establish the dominant toxin hypothesis which incorporates observations in other venomous lineages that animals have convergently evolved a similar strategy in shaping their venom.

Natural promoters and promoter engineering strategies for metabolic regulation in Saccharomyces cerevisiae

Sharomyces cerevisiae is currently one of the most important foreign gene expression systems. S. cerevisiae is an excellent host for high-value metabolite cell factories due to its advantages of simplicity, safety, and nontoxicity. A promoter, as one of the basic elements of gene transcription, plays an important role in regulating gene expression and optimizing metabolic pathways. Promoters control the direction and intensity of transcription, and the application of promoters with different intensities and performances will largely determine the effect of gene expression and ultimately affect the experimental results. Due to its significant role, there have been many studies on promoters for decades. While some studies have explored and analyzed new promoters with different functions, more studies have focused on artificially modifying promoters to meet their own scientific needs. Thus, this article reviews current research on promoter engineering techniques and related natural promoters in S. cerevisiae. First, we introduce the basic structure of promoters and the classification of natural promoters. Then, the classification of various promoter strategies is reviewed. Finally, by grouping related articles together using various strategies, this review anticipates the future development direction of promoter engineering.

FACS-Based Sequencing Approach to Evaluate Cell Type to Genotype Associations Using Cerebral Organoids

Recent technological developments have led to widespread applications of large-scale transcriptomics-based sequencing methods to identify genotype-to-cell type associations. Here we describe a fluorescence-activated cell sorting (FACS)-based sequencing method to utilize CRISPR/Cas9 edited mosaic cerebral organoids to identify or validate genotype-to-cell type associations. Our approach is high-throughput and quantitative and uses internal controls to enable comparisons of the results across different antibody markers and experiments.

Biologically relevant small variations of intra-cellular pH can have significant effect on stability of protein-DNA complexes, including the nucleosome

Stability of a protein-ligand complex may be sensitive to pH of its environment. Here we explore, computationally, stability of a set of protein-nucleic acid complexes using fundamental thermodynamic linkage relationship. The nucleosome, as well as an essentially random selection of 20 protein complexes with DNA or RNA, are included in the analysis. An increase in intra-cellular/intra-nuclear pH destabilizes most complexes, including the nucleosome. We propose to quantify the effect by ΔΔG^0.3-the change in the binding free energy due to pH increase of 0.3 units, corresponding to doubling of the H ⁺ activity; variations of pH of this amplitude can occur in living cells, including in the course of the cell cycle, and in cancer cells relative to normal ones. We suggest, based on relevant experimental findings, a threshold of biological significance of 1 2 k B T ( ∼ 0.3 k c a l / m o l ) for changes of stability of chromatin-related protein-DNA complexes: a change in the binding affinity above the threshold may have biological consequences. We find that for 70% of the examined complexes, Δ Δ G 0.3 > 1 2 k B T (for 10%, ΔΔG^0.3 is between 3 and 4 k _B T). Thus, small but relevant variations of intra-nuclear pH of 0.3 may have biological consequences for many protein-nucleic acid complexes. The binding affinity between the histone octamer and its DNA, which directly affects the DNA accessibility in the nucleosome, is predicted to be highly sensitive to intra-nuclear pH. A variation of 0.3 units results in ΔΔG^0.3 ∼ 10k _B T ( ∼ 6 k c a l / m o l ) ; for spontaneous unwrapping of 20 bp long entry/exit fragments of the nucleosomal DNA, ΔΔG^0.3 = 2.2k _B T; partial disassembly of the nucleosome into the tetrasome is characterized by ΔΔG^0.3 = 5.2k _B T. The predicted pH -induced modulations of the nucleosome stability are significant enough to suggest that they may have consequences relevant to the biological function of the nucleosome. Accessibility of the nucleosomal DNA is predicted to positively correlate with pH variations during the cell cycle; an increase in intra-cellular pH seen in cancer cells is predicted to lead to a more accessible nucleosomal DNA; a drop in pH associated with apoptosis is predicted to make nucleosomal DNA less accessible. We speculate that processes that depend on accessibility to the DNA in the nucleosomes, such as transcription or DNA replication, might become upregulated due to relatively small, but nevertheless realistic increases of intra-nuclear pH.

Construction and Optimization of Malonyl-CoA Sensors in Saccharomyces cerevisiae by Combining Promoter Engineering Strategies

Biosensors can be used for high-throughput screening, real-time monitoring of metabolites, and dynamic regulation of metabolic processes, which have been a popular research direction in recent years. Here, five promoters from Saccharomyces cerevisiae were selected to construct Malonyl-CoA sensors with the fapO/fapR system derived from Bacillus subtilis, and pCCW12 was finally selected for further optimization. Based on pCCW12, a series of sensors with different response sensitivities were obtained by selecting different fapO insertion sites and combining the best two or three of them. Then, through a combination of promoter hybrid, intron insertion, and transcription factor modification strategies, we obtained sensors with different effects, one of which, the H-pCCW12(TFBS)-Cti6~fapR sensor, had the lowest background noise, doubled response range and higher response sensitivity compared to the original sensor. Sensors with different characteristics constructed in this study, can be applied to Malonyl-CoA related high-throughput screening and finer regulation of metabolism. It also proves that the combined application of different promoter engineering strategies is a feasible idea for the precise construction and regulation of biosensors.

Manipulating chromatin architecture in C. elegans

BACKGROUND

Nucleosome-mediated chromatin compaction has a direct effect on the accessibility of trans-acting activators and repressors to DNA targets and serves as a primary regulatory agent of genetic expression. Understanding the nature and dynamics of chromatin is fundamental to elucidating the mechanisms and factors that epigenetically regulate gene expression. Previous work has shown that there are three types of canonical sequences that strongly regulate nucleosome positioning and thus chromatin accessibility: putative nucleosome-positioning elements, putative nucleosome-repelling sequences, and homopolymeric runs of A/T. It is postulated that these elements can be used to remodel chromatin in C. elegans. Here we show the utility of such elements in vivo, and the extreme efficacy of a newly discovered repelling sequence, PRS-322.

RESULTS

In this work, we show that it is possible to manipulate nucleosome positioning in C. elegans solely using canonical and putative positioning sequences. We have not only tested previously described sequences such as the Widom 601, but also have tested additional nucleosome-positioning sequences: the Trifonov sequence, putative repelling sequence-322 (PRS-322), and various homopolymeric runs of A and T nucleotides.

CONCLUSIONS

Using each of these types of putative nucleosome-positioning sequences, we demonstrate their ability to alter the nucleosome profile in C. elegans as evidenced by altered nucleosome occupancy and positioning in vivo. Additionally, we show the effect that PRS-322 has on nucleosome-repelling and chromatin remodeling.

A Poly(dA:dT) Tract in the IGF1 Gene Is a Genetic Marker for Growth Traits in Pigs

Insulin-like growth factor 1 (IGF1) is an important regulator of body growth, development, and metabolism. The poly(dA:dT) tract affects the accessibility of transcription factor binding sites to regulate transcription. Therefore, this study assessed the effects of two poly(dA:dT) tracts on the transcriptional activity of porcine IGF1. The luciferase assay results demonstrated that the poly(dA:dT) tract 2 (−264/−255) was a positive regulatory element for IGF1 gene expression, and the activities between the different lengths of the poly(dA:dT) tract 2 were significant (p<0.01). The transcription factor C/EBPα inhibited the transcription of IGF1 by binding to tract 2, and the expression levels between the lengths of tract 2 after C/EBPα binding were also statistically different (p<0.01). Only the alleles 10T and 11T were found in the tract 2 in commercial pig breeds, while the 9T, 10T, and 11T alleles were found in Chinese native pig breeds. The allele frequencies were in Hardy−Weinberg equilibrium in all pig breeds. The genotypes of tract 2 were significantly associated with the growth traits (days to 115 kg and average daily gain) (p<0.05) in commercial pig breeds. Based on these findings, it can be concluded that the tract 2 mutation could be applied as a candidate genetic marker for growth trait selection in pig breeding programs.

Nucleus size and its effect on nucleosome stability in living cells

DNA architectural proteins play a major role in organization of chromosomal DNA in living cells by packaging it into chromatin, whose spatial conformation is determined by an intricate interplay between the DNA-binding properties of architectural proteins and physical constraints applied to the DNA by a tight nuclear space. Yet, the exact effects of the nucleus size on DNA-protein interactions and chromatin structure currently remain obscure. Furthermore, there is even no clear understanding of molecular mechanisms responsible for the nucleus size regulation in living cells. To find answers to these questions, we developed a general theoretical framework based on a combination of polymer field theory and transfer-matrix calculations, which showed that the nucleus size is mainly determined by the difference between the surface tensions of the nuclear envelope and the endoplasmic reticulum membrane as well as the osmotic pressure exerted by cytosolic macromolecules on the nucleus. In addition, the model demonstrated that the cell nucleus functions as a piezoelectric element, changing its electrostatic potential in a size-dependent manner. This effect has been found to have a profound impact on stability of nucleosomes, revealing a previously unknown link between the nucleus size and chromatin structure. Overall, our study provides new insights into the molecular mechanisms responsible for regulation of the nucleus size, as well as the potential role of nuclear organization in shaping the cell response to environmental cues.

Nucleosome sliding can influence the spreading of histone modifications

Nucleosomes are the fundamental building blocks of chromatin that not only help in the folding of chromatin, but also in carrying epigenetic information. It is known that nucleosome sliding is responsible for dynamically organizing chromatin structure and the resulting gene regulation. Since sliding can move two neighboring nucleosomes physically close or away, can it play a role in the spreading of histone modifications? We investigate this by simulating a stochastic model that couples nucleosome dynamics with the kinetics of histone modifications. We show that the sliding of nucleosomes can affect the modification pattern as well as the time it takes to modify a given region of chromatin. Exploring different nucleosome densities and modification kinetic parameters, we show that nucleosome sliding can be important for creating histone modification domains. Our model predicts that nucleosome density coupled with sliding dynamics can create an asymmetric histone modification profile around regulatory regions. We also compute the probability distribution of modified nucleosomes and relaxation kinetics of modifications. Our predictions are comparable with known experimental results.

High-throughput techniques enable advances in the roles of DNA and RNA secondary structures in transcriptional and post-transcriptional gene regulation

The most stable structure of DNA is the canonical right-handed double helix termed B DNA. However, certain environments and sequence motifs favor alternative conformations, termed non-canonical secondary structures. The roles of DNA and RNA secondary structures in transcriptional regulation remain incompletely understood. However, advances in high-throughput assays have enabled genome wide characterization of some secondary structures. Here, we describe their regulatory functions in promoters and 3’UTRs, providing insights into key mechanisms through which they regulate gene expression. We discuss their implication in human disease, and how advances in molecular technologies and emerging high-throughput experimental methods could provide additional insights.

Human Virus Genomes Are Enriched in Conserved Adenine/Thymine/Uracil Multiple Tracts That Pause Polymerase Progression

The DNA secondary structures that deviate from the classic Watson and Crick base pairing are increasingly being reported to form transiently in the cell and regulate specific cellular mechanisms. Human viruses are cell parasites that have evolved mechanisms shared with the host cell to support their own replication and spreading. Contrary to human host cells, viruses display a diverse array of nucleic acid types, which include DNA or RNA in single-stranded or double-stranded conformations. This heterogeneity improves the possible occurrence of non-canonical nucleic acid structures. We have previously shown that human virus genomes are enriched in G-rich sequences that fold in four-stranded nucleic acid secondary structures, the G-quadruplexes.Here, by extensive bioinformatics analysis on all available genomes, we showed that human viruses are enriched in highly conserved multiple A (and T or U) tracts, with such an array that they could in principle form quadruplex structures. By circular dichroism, NMR, and Taq polymerase stop assays, we proved that, while A/T/U-quadruplexes do not form, these tracts still display biological significance, as they invariably trigger polymerase pausing within two bases from the A/T/U tract. “A” bases display the strongest effect. Most of the identified A-tracts are in the coding strand, both at the DNA and RNA levels, suggesting their possible relevance during viral translation. This study expands on the presence and mechanism of nucleic acid secondary structures in human viruses and provides a new direction for antiviral research.

Transcriptional neighborhoods regulate transcript isoform lengths and expression levels

Sequence features of genes and their flanking regulatory regions are determinants of RNA transcript isoform expression and have been used as context-independent plug-and-play modules in synthetic biology. However, genetic context-including the adjacent transcriptional environment-also influences transcript isoform expression levels and boundaries. We used synthetic yeast strains with stochastically repositioned genes to systematically disentangle the effects of sequence and context. Profiling 120 million full-length transcript molecules across 612 genomic perturbations, we observed sequence-independent alterations to gene expression levels and transcript isoform boundaries that were influenced by neighboring transcription. We identified features of transcriptional context that could predict these alterations and used these features to engineer a synthetic circuit where transcript length was controlled by neighboring transcription. This demonstrates how positional context can be leveraged in synthetic genome engineering.

The Active Mechanism of Nucleosome Depletion by Poly(dA:dT) Tracts In Vivo

Poly(dA:dT) tracts cause nucleosome depletion in many species, e.g., at promoters and replication origins. Their intrinsic biophysical sequence properties make them stiff and unfavorable for nucleosome assembly, as probed by in vitro nucleosome reconstitution. The mere correlation between nucleosome depletion over poly(dA:dT) tracts in in vitro reconstituted and in in vivo chromatin inspired an intrinsic nucleosome exclusion mechanism in vivo that is based only on DNA and histone properties. However, we compile here published and new evidence that this correlation does not reflect mechanistic causation. (1) Nucleosome depletion over poly(dA:dT) in vivo is not universal, e.g., very weak in S. pombe. (2) The energy penalty for incorporating poly(dA:dT) tracts into nucleosomes is modest (<10%) relative to ATP hydrolysis energy abundantly invested by chromatin remodelers. (3) Nucleosome depletion over poly(dA:dT) is much stronger in vivo than in vitro if monitored without MNase and (4) actively maintained in vivo. (5) S. cerevisiae promoters evolved a strand-biased poly(dA) versus poly(dT) distribution. (6) Nucleosome depletion over poly(dA) is directional in vivo. (7) The ATP dependent chromatin remodeler RSC preferentially and directionally displaces nucleosomes towards 5′ of poly(dA). Especially distribution strand bias and displacement directionality would not be expected for an intrinsic mechanism. Together, this argues for an in vivo mechanism where active and species-specific read out of intrinsic sequence properties, e.g., by remodelers, shapes nucleosome organization.

Saccharomyces cerevisiae Promoter Engineering before and during the Synthetic Biology Era

Simple Summary

Promoters are DNA sequences where the process of transcription starts. They can work constitutively or be controlled by environmental signals of different types. The quantity of proteins and RNA present in yeast genetic circuits highly depends on promoter strength. Hence, they have been deeply studied and modified over, at least, the last forty years, especially since the year 2000 when Synthetic Biology was born. Here, we present how promoter engineering changed over these four decades and discuss its possible future directions due to novel computational methods and technology.

Abstract

Synthetic gene circuits are made of DNA sequences, referred to as transcription units, that communicate by exchanging proteins or RNA molecules. Proteins are, mostly, transcription factors that bind promoter sequences to modulate the expression of other molecules. Promoters are, therefore, key components in genetic circuits. In this review, we focus our attention on the construction of artificial promoters for the yeast S. cerevisiae, a popular chassis for gene circuits. We describe the initial techniques and achievements in promoter engineering that predated the start of the Synthetic Biology epoch of about 20 years. We present the main applications of synthetic promoters built via different methods and discuss the latest innovations in the wet-lab engineering of novel promoter sequences.

Coupling between Sequence-Mediated Nucleosome Organization and Genome Evolution

The nucleosome is a major modulator of DNA accessibility to other cellular factors. Nucleosome positioning has a critical importance in regulating cell processes such as transcription, replication, recombination or DNA repair. The DNA sequence has an influence on the position of nucleosomes on genomes, although other factors are also implicated, such as ATP-dependent remodelers or competition of the nucleosome with DNA binding proteins. Different sequence motifs can promote or inhibit the nucleosome formation, thus influencing the accessibility to the DNA. Sequence-encoded nucleosome positioning having functional consequences on cell processes can then be selected or counter-selected during evolution. We review the interplay between sequence evolution and nucleosome positioning evolution. We first focus on the different ways to encode nucleosome positions in the DNA sequence, and to which extent these mechanisms are responsible of genome-wide nucleosome positioning in vivo. Then, we discuss the findings about selection of sequences for their nucleosomal properties. Finally, we illustrate how the nucleosome can directly influence sequence evolution through its interactions with DNA damage and repair mechanisms. This review aims to provide an overview of the mutual influence of sequence evolution and nucleosome positioning evolution, possibly leading to complex evolutionary dynamics.

Novel S. cerevisiae Hybrid Synthetic Promoters Based on Foreign Core Promoter Sequences

Promoters are fundamental components of synthetic gene circuits. They are DNA segments where transcription initiation takes place. New constitutive and regulated promoters are constantly engineered in order to meet the requirements for protein and RNA expression into different genetic networks. In this work, we constructed and optimized new synthetic constitutive promoters for the yeast Saccharomyces cerevisiae. We started from foreign (e.g., viral) core promoters as templates. They are, usually, unfunctional in yeast but can be activated by extending them with a short sequence, from the CYC1 promoter, containing various transcription start sites (TSSs). Transcription was modulated by mutating the TATA box composition and varying its distance from the TSS. We found that gene expression is maximized when the TATA box has the form TATAAAA or TATATAA and lies between 30 and 70 nucleotides upstream of the TSS. Core promoters were turned into stronger promoters via the addition of a short UAS. In particular, the 40 nt bipartite UAS from the GPD promoter can enhance protein synthesis considerably when placed 150 nt upstream of the TATA box. Overall, we extended the pool of S. cerevisiae promoters with 59 new samples, the strongest overcoming the native TEF2 promoter.

Molecular and evolutionary processes generating variation in gene expression

Heritable variation in gene expression is common within and between species. This variation arises from mutations that alter the form or function of molecular gene regulatory networks that are then filtered by natural selection. High-throughput methods for introducing mutations and characterizing their cis- and trans-regulatory effects on gene expression (particularly, transcription) are revealing how different molecular mechanisms generate regulatory variation, and studies comparing these mutational effects with variation seen in the wild are teasing apart the role of neutral and non-neutral evolutionary processes. This integration of molecular and evolutionary biology allows us to understand how the variation in gene expression we see today came to be and to predict how it is most likely to evolve in the future.

Capturing and Understanding the Dynamics and Heterogeneity of Gene Expression in the Living Cell

The regulation of gene expression is a fundamental process enabling cells to respond to internal and external stimuli or to execute developmental programs. Changes in gene expression are highly dynamic and depend on many intrinsic and extrinsic factors. In this review, we highlight the dynamic nature of transient gene expression changes to better understand cell physiology and development in general. We will start by comparing recent in vivo procedures to capture gene expression in real time. Intrinsic factors modulating gene expression dynamics will then be discussed, focusing on chromatin modifications. Furthermore, we will dissect how cell physiology or age impacts on dynamic gene regulation and especially discuss molecular insights into acquired transcriptional memory. Finally, this review will give an update on the mechanisms of heterogeneous gene expression among genetically identical individual cells. We will mainly focus on state-of-the-art developments in the yeast model but also cover higher eukaryotic systems.

Understanding the paradoxical mechanical response of in-phase A-tracts at different force regimes

A-tracts are A:T rich DNA sequences that exhibit unique structural and mechanical properties associated with several functions in vivo. The crystallographic structure of A-tracts has been well characterized. However, the mechanical properties of these sequences is controversial and their response to force remains unexplored. Here, we rationalize the mechanical properties of in-phase A-tracts present in the Caenorhabditis elegans genome over a wide range of external forces, using single-molecule experiments and theoretical polymer models. Atomic Force Microscopy imaging shows that A-tracts induce long-range (∼200 nm) bending, which originates from an intrinsically bent structure rather than from larger bending flexibility. These data are well described with a theoretical model based on the worm-like chain model that includes intrinsic bending. Magnetic tweezers experiments show that the mechanical response of A-tracts and arbitrary DNA sequences have a similar dependence with monovalent salt supporting that the observed A-tract bend is intrinsic to the sequence. Optical tweezers experiments reveal a high stretch modulus of the A-tract sequences in the enthalpic regime. Our work rationalizes the complex multiscale flexibility of A-tracts, providing a physical basis for the versatile character of these sequences inside the cell.

Quantitative analysis of transcription factor binding and expression using calling cards reporter arrays

We report a tool, Calling Cards Reporter Arrays (CCRA), that measures transcription factor (TF) binding and the consequences on gene expression for hundreds of synthetic promoters in yeast. Using Cbf1p and MAX, we demonstrate that the CCRA method is able to detect small changes in binding free energy with a sensitivity comparable to in vitro methods, enabling the measurement of energy landscapes in vivo. We then demonstrate the quantitative analysis of cooperative interactions by measuring Cbf1p binding at synthetic promoters with multiple sites. We find that the cooperativity between Cbf1p dimers varies sinusoidally with a period of 10.65 bp and energetic cost of 1.37 K_BT for sites that are positioned ‘out of phase’. Finally, we characterize the binding and expression of a group of TFs, Tye7p, Gcr1p and Gcr2p, that act together as a ‘TF collective’, an important but poorly characterized model of TF cooperativity. We demonstrate that Tye7p often binds promoters without its recognition site because it is recruited by other collective members, whereas these other members require their recognition sites, suggesting a hierarchy where these factors recruit Tye7p but not vice versa. Our experiments establish CCRA as a useful tool for quantitative investigations into TF binding and function.

Promoter Architecture and Promoter Engineering in Saccharomyces cerevisiae

Promoters play an essential role in the regulation of gene expression for fine-tuning genetic circuits and metabolic pathways in Saccharomyces cerevisiae (S. cerevisiae). However, native promoters in S. cerevisiae have several limitations which hinder their applications in metabolic engineering. These limitations include an inadequate number of well-characterized promoters, poor dynamic range, and insufficient orthogonality to endogenous regulations. Therefore, it is necessary to perform promoter engineering to create synthetic promoters with better properties. Here, we review recent advances related to promoter architecture, promoter engineering and synthetic promoter applications in S. cerevisiae. We also provide a perspective of future directions in this field with an emphasis on the recent advances of machine learning based promoter designs.

Transcription in Living Cells: Molecular Mechanisms of Bursting

Transcription in several organisms from certain bacteria to humans has been observed to be stochastic in nature: toggling between active and inactive states. Periods of active nascent RNA synthesis known as bursts represent individual gene activation events in which multiple polymerases are initiated. Therefore, bursting is the single locus illustration of both gene activation and repression. Although transcriptional bursting was originally observed decades ago, only recently have technological advances enabled the field to begin elucidating gene regulation at the single-locus level. In this review, we focus on how biochemical, genomic, and single-cell data describe the regulatory steps of transcriptional bursts.

Modulating Pathway Performance by Perturbing Local Genetic Context

Combinatorial engineering is a preferred strategy for attaining optimal pathway performance. Previous endeavors have been concentrated on regulatory elements (e.g., promoters, terminators, and ribosomal binding sites) and/or open reading frames. Accumulating evidence indicates that noncoding DNA sequences flanking a transcriptional unit on the genome strongly impact gene expression. Here, we sought to mimic the effect imposed on expression cassettes by the genome. We created variants of the model yeast Saccharomyces cerevisiae with significantly improved fluorescence or cellobiose consumption rate by randomizing the sequences adjacent to the GFP expression cassette or the cellobiose-utilization pathway, respectively. Interestingly, nucleotide specificity was observed at certain positions and showed to be essential for achieving optimal cellobiose assimilation. Further characterization suggested that the modulation effects of the short sequences flanking the expression cassettes could be potentially mediated by remodeling DNA packaging and/or recruiting transcription factors. Collectively, these results indicate that the often-overlooked contiguous DNA sequences can be exploited to rapidly achieve balanced pathway expression, and the corresponding approach could be easily stacked with other combinatorial engineering strategies.

Polymorphism of simple sequence repeats may quantitatively regulate gene transcription

The degree of polymorphism, i.e., DNA sequence divergence, of short AT-rich tandemly arranged simple sequence repeats at or near promoters and 5'- untranslated regions of mRNA may quantitatively regulate transcription of tissue-specific genes. Less polymorphic repeats allow greater gene expression. Preferential binding of hypophosphorylated H1 histone to these repeats may diminish binding of transcription factors. Preferential binding of hypophosphorylated high mobility group chromatin proteins would increase this binding. Shorter simple sequence repeats have undergone fewer point mutations than longer repeats, hence they are less polymorphic and more conserved. The role of transcribed simple sequence repeats in frog embryo germ layer determination is considered.

Sensitive Automated Measurement of Histone-DNA Affinities in Nucleosomes

The DNA of eukaryotes is wrapped around histone octamers to form nucleosomes. Although it is well established that the DNA sequence significantly influences nucleosome formation, its precise contribution has remained controversial, partially owing to the lack of quantitative affinity data. Here, we present a method to measure DNA-histone binding free energies at medium throughput and with high sensitivity. Competitive nucleosome formation is achieved through automation, and a modified epifluorescence microscope is used to rapidly and accurately measure the fractions of bound/unbound DNA based on fluorescence anisotropy. The procedure allows us to obtain full titration curves with high reproducibility. We applied this technique to measure the histone-DNA affinities for 47 DNA sequences and analyzed how the affinities correlate with relevant DNA sequence features. We found that the GC content has a significant impact on nucleosome-forming preferences, but 10 bp dinucleotide periodicities and the presence of poly(dA:dT) stretches do not.

•

Robotics permits full titration series to measure histone-DNA binding affinities

•

Fluorescence anisotropy used as a fast, sensitive readout of bound/unbound DNA

•

Free energies span three orders of magnitude, less for naturally occurring sequences

•

GC content is a major determinant of measured binding free energies

Engineering yeast artificial core promoter with designated base motifs

BACKGROUND

Synthetic biology requires toolbox of promoters to finely tune gene expression levels for building up efficient cell factories. Yeast promoters owned variable core promoter regions between the TATA-box and transcriptional starting site (TSS) at the length mostly around 20-80 bases. This region allowed flexible design of artificial promoter but potentially demand special base motifs to maintain or enhance the promoter's strength.

RESULTS

Here, we designed and screened the base motifs and tested the activities of yeast artificial core promoters. Different 30 bases of artificial sequences led to variable expression levels of CrtY enzyme which determined the lycopene-carotene compositions, represented in the colony-color spectrum of red-orange-yellow. The upstream sequences of two strong promoter P_EXP1 and P_GPD and two starting strains with distinguishable lycopene production levels were utilized to characterize the promoter sequences. Different partition designs of T-rich or G/C-rich base motifs led to distinguishable colony-color distributions. Finally, we screened a champion promoter with a highest 5.5-fold enhancement of lycopene-carotene transformation. Another selected promoter generated a highest beta-carotene production as 7.4 mg/g DCW.

CONCLUSIONS

This work offered an approach to redesign promoter with artificial sequences. We concluded that the core promoter region could be designated as 30 bases and different base motifs would enhance or weaken the promoter's strength. Generally, more T-rich elements, higher %T and lower G/C percentage were beneficial to enhance the strength of artificial core promoter.

SerraNA: a program to determine nucleic acids elasticity from simulation data

AT-rich motifs can generate extreme mechanical properties, which are critical for creating strong global bends when phased properly.

Nucleosome Dynamics: a new tool for the dynamic analysis of nucleosome positioning

We present Nucleosome Dynamics, a suite of programs integrated into a virtual research environment and created to define nucleosome architecture and dynamics from noisy experimental data. The package allows both the definition of nucleosome architectures and the detection of changes in nucleosomal organization due to changes in cellular conditions. Results are displayed in the context of genomic information thanks to different visualizers and browsers, allowing the user a holistic, multidimensional view of the genome/transcriptome. The package shows good performance for both locating equilibrium nucleosome architecture and nucleosome dynamics and provides abundant useful information in several test cases, where experimental data on nucleosome position (and for some cases expression level) have been collected for cells under different external conditions (cell cycle phase, yeast metabolic cycle progression, changes in nutrients or difference in MNase digestion level). Nucleosome Dynamics is a free software and is provided under several distribution models.

Modulating transcription through development of semi-synthetic yeast core promoters

Altering gene expression regulation by promoter engineering is a very effective way to fine-tune heterologous pathways in eukaryotic hosts. Typically, pathway building approaches in yeast still use a limited set of long, native promoters. With the today’s introduction of longer and more complex pathways, an expansion of this synthetic biology toolbox is necessary. In this study we elucidated the core promoter structure of the well-characterized yeast TEF1 promoter and determined the minimal length needed for sufficient protein expression. Furthermore, this minimal core promoter sequence was used for the creation of a promoter library covering different expression strengths. This resulted in a group of short, 69 bp promoters with an 8.0-fold expression range. One exemplar had a two and four times higher expression compared to the native CYC1 and ADH1 promoter, respectively. Additionally, as it was described that the protein expression range could be broadened by upstream activating sequences (UASs), we integrated earlier described single and multiple short, synthetic UASs in front of the strongest yeast core promoter. This approach resulted to further variation in protein expression and an overall promoter library spanning a 20-fold activity range and covering a length from 69 bp to maximally 129 bp. Furthermore, the robustness of this library was assessed on three alternative carbon sources besides glucose. As such, the suitability of short yeast core promoters for metabolic engineering applications on different media, either in an individual context or combined with UAS elements, was demonstrated.

The impact of short tandem repeat variation on gene expression

Short tandem repeats (STRs) have been implicated in a variety of complex traits in humans. However, genome-wide studies of the effects of STRs on gene expression thus far have had limited power to detect associations and provide insights into putative mechanisms. Here, we leverage whole-genome sequencing and expression data for 17 tissues from the Genotype-Tissue Expression Project to identify more than 28,000 STRs for which repeat number is associated with expression of nearby genes (eSTRs). We use fine-mapping to quantify the probability that each eSTR is causal and characterize the top 1,400 fine-mapped eSTRs. We identify hundreds of eSTRs linked with published genome-wide association study signals and implicate specific eSTRs in complex traits, including height, schizophrenia, inflammatory bowel disease and intelligence. Overall, our results support the hypothesis that eSTRs contribute to a range of human phenotypes, and our data should serve as a valuable resource for future studies of complex traits.

Clock-controlled rhythmic transcription: is the clock enough and how does it work?

Circadian clocks regulate the rhythmic expression of thousands of genes underlying the daily oscillations of biological functions. Here, we discuss recent findings showing that circadian clock rhythmic transcriptional outputs rely on additional mechanisms than just clock gene DNA binding, which may ultimately contribute to the plasticity of circadian transcriptional programs.

Chromatin accessibility and the regulatory epigenome

Physical access to DNA is a highly dynamic property of chromatin that plays an essential role in establishing and maintaining cellular identity. The organization of accessible chromatin across the genome reflects a network of permissible physical interactions through which enhancers, promoters, insulators and chromatin-binding factors cooperatively regulate gene expression. This landscape of accessibility changes dynamically in response to both external stimuli and developmental cues, and emerging evidence suggests that homeostatic maintenance of accessibility is itself dynamically regulated through a competitive interplay between chromatin-binding factors and nucleosomes. In this Review, we examine how the accessible genome is measured and explore the role of transcription factors in initiating accessibility remodelling; our goal is to illustrate how chromatin accessibility defines regulatory elements within the genome and how these epigenetic features are dynamically established to control gene expression.

Nucleosome positioning and spacing: from genome-wide maps to single arrays

The positioning of nucleosomes relative to DNA and their neighboring nucleosomes represents a fundamental layer of chromatin organization. Changes in nucleosome positioning and spacing affect the accessibility of DNA to regulatory factors and the formation of higher order chromatin structures. Sequencing of mononucleosomal fragments allowed mapping nucleosome positions on a genome-wide level in many organisms. This revealed that successions of evenly spaced and well-positioned nucleosomes-so called phased nucleosome arrays-occur at the 5' end of many active genes and in the vicinity of transcription factor and other protein binding sites. Phased arrays arise from the interplay of barrier elements on the DNA, which position adjacent nucleosomes, and the nucleosome spacing activity of ATP-dependent chromatin remodelers. A shortcoming of classic mononucleosomal mapping experiments is that they only reveal nucleosome spacing and array regularity at select sites in the genome with well-positioned nucleosomes. However, new technological approaches elucidate nucleosome array structure throughout the genome and with single-cell resolution. In the future, it will be interesting to see whether changes in nucleosome array regularity and spacing contribute to the formation of higher order chromatin structures and the spatial organization of the genome in vivo.

A k-mer grammar analysis to uncover maize regulatory architecture

BACKGROUND

Only a small percentage of the genome sequence is involved in regulation of gene expression, but to biochemically identify this portion is expensive and laborious. In species like maize, with diverse intergenic regions and lots of repetitive elements, this is an especially challenging problem that limits the use of the data from one line to the other. While regulatory regions are rare, they do have characteristic chromatin contexts and sequence organization (the grammar) with which they can be identified.

RESULTS

We developed a computational framework to exploit this sequence arrangement. The models learn to classify regulatory regions based on sequence features - k-mers. To do this, we borrowed two approaches from the field of natural language processing: (1) "bag-of-words" which is commonly used for differentially weighting key words in tasks like sentiment analyses, and (2) a vector-space model using word2vec (vector-k-mers), that captures semantic and linguistic relationships between words. We built "bag-of-k-mers" and "vector-k-mers" models that distinguish between regulatory and non-regulatory regions with an average accuracy above 90%. Our "bag-of-k-mers" achieved higher overall accuracy, while the "vector-k-mers" models were more useful in highlighting key groups of sequences within the regulatory regions.

CONCLUSIONS

These models now provide powerful tools to annotate regulatory regions in other maize lines beyond the reference, at low cost and with high accuracy.

Systematic interrogation of human promoters

Despite much research, our understanding of the architecture and cis-regulatory elements of human promoters is still lacking. Here, we devised a high-throughput assay to quantify the activity of approximately 15,000 fully designed sequences that we integrated and expressed from a fixed location within the human genome. We used this method to investigate thousands of native promoters and preinitiation complex (PIC) binding regions followed by in-depth characterization of the sequence motifs underlying promoter activity, including core promoter elements and TF binding sites. We find that core promoters drive transcription mostly unidirectionally and that sequences originating from promoters exhibit stronger activity than those originating from enhancers. By testing multiple synthetic configurations of core promoter elements, we dissect the motifs that positively and negatively regulate transcription as well as the effect of their combinations and distances, including a 10-bp periodicity in the optimal distance between the TATA and the initiator. By comprehensively screening 133 TF binding sites, we find that in contrast to core promoters, TF binding sites maintain similar activity levels in both orientations, supporting a model by which divergent transcription is driven by two distinct unidirectional core promoters sharing bidirectional TF binding sites. Finally, we find a striking agreement between the effect of binding site multiplicity of individual TFs in our assay and their tendency to appear in homotypic clusters throughout the genome. Overall, our study systematically assays the elements that drive expression in core and proximal promoter regions and sheds light on organization principles of regulatory regions in the human genome.