Structural properties of promoters: similarities and differences between prokaryotes and eukaryotes

Augusta: From RNA-Seq to gene regulatory networks and Boolean models

Computational models of gene regulations help to understand regulatory mechanisms and are extensively used in a wide range of areas, e.g., biotechnology or medicine, with significant benefits. Unfortunately, there are only a few computational gene regulatory models of whole genomes allowing static and dynamic analysis due to the lack of sophisticated tools for their reconstruction. Here, we describe Augusta, an open-source Python package for Gene Regulatory Network (GRN) and Boolean Network (BN) inference from the high-throughput gene expression data. Augusta can reconstruct genome-wide models suitable for static and dynamic analyses. Augusta uses a unique approach where the first estimation of a GRN inferred from expression data is further refined by predicting transcription factor binding motifs in promoters of regulated genes and by incorporating verified interactions obtained from databases. Moreover, a refined GRN is transformed into a draft BN by searching in the curated model database and setting logical rules to incoming edges of target genes, which can be further manually edited as the model is provided in the SBML file format. The approach is applicable even if information about the organism under study is not available in the databases, which is typically the case for non-model organisms including most microbes. Augusta can be operated from the command line and, thus, is easy to use for automated prediction of models for various genomes. The Augusta package is freely available at github.com/JanaMus/Augusta. Documentation and tutorials are available at augusta.readthedocs.io.

From straight to curved: A historical perspective of DNA shape

DNA is the macromolecule responsible for storing the genetic information of a cell and it has intrinsic properties such as deformability, stability and curvature. DNA Curvature plays an important role in gene transcription and, consequently, in the subsequent production of proteins, a fundamental process of cells. With recent advances in bioinformatics and theoretical biology, it became possible to analyze and understand the involvement of DNA Curvature as a discriminatory characteristic of gene-promoting regions. These regions act as sites where RNAp (ribonucleic acid-polymerase) binds to initiate transcription. This review aims to describe the formation of Curvature, as well as highlight its importance in predicting promoters. Furthermore, this article provides the potential of DNA Curvature as a distinguishing feature for promoter prediction tools, as well as outlining the calculation procedures that have been described by other researchers. This work may support further studies directed towards the enhancement of promoter prediction software.

Identification of an Endogenous Strong Promoter in Burkholderia sp. JP2-270

Burkholderia is the second largest source of natural product bacteria after Actinomyces and can produce many secondary metabolites including pyrrolnitrin (PRN). Natural products of microbial origin are usually found in trace amounts, so in metabolic engineering, promoter engineering is often used to regulate gene expression to increase yield. In this study, an endogenous strong promoter was identified based on RNA-seq to overexpress biosynthetic genes to increase the production of PRN. By analyzing the transcriptomic data of the antagonistic bacterium Burkholderia sp. JP2-270 in three different development periods, we screened 50 endogenous promoters with high transcriptional activity, nine of which were verified by an obvious fluorescent signal via fluorescence observation. Then, combined with RT-qPCR analysis, Php, the promoter of a hypothetical protein, was found to be significantly expressed in all three periods. In order to increase the suitability of endogenous promoters, the promoter Php was shortened at different lengths, and the results show that a sequence length of 173 bp was necessary for its activity. Moreover, this promoter was used to overexpress the PRN biosynthesis genes (prnA, prnB, prnC and prnD) in JP2-270, resulting in a successful increase in gene expression levels by 40-80 times. Only the overexpression of the prnB gene successfully increased PRN production to 1.46 times that of the wild type. Overall, the endogenous strong promoters screened in this study can improve gene expression and increase the production of secondary metabolites in JP2-270 and other strains.

DNASCANNER v2: A Web-Based Tool to Analyze the Characteristic Properties of Nucleotide Sequences

Throughout the process of evolution, DNA undergoes the accumulation of distinct mutations, which can often result in highly organized patterns that serve various essential biological functions. These patterns encompass various genomic elements and provide valuable insights into the regulatory and functional aspects of DNA. The physicochemical, mechanical, thermodynamic, and structural properties of DNA sequences play a crucial role in the formation of specific patterns. These properties contribute to the three-dimensional structure of DNA and influence their interactions with proteins, regulatory elements, and other molecules. In this study, we introduce DNASCANNER v2, an advanced version of our previously published algorithm DNASCANNER for analyzing DNA properties. The current tool is built using the FLASK framework in Python language. Featuring a user-friendly interface tailored for nonspecialized researchers, it offers an extensive analysis of 158 DNA properties, including mono/di/trinucleotide frequencies, structural, physicochemical, thermodynamics, and mechanical properties of DNA sequences. The tool provides downloadable results and offers interactive plots for easy interpretation and comparison between different features. We also demonstrate the utility of DNASCANNER v2 in analyzing splice-site junctions, casposon insertion sequences, and transposon insertion sites (TIS) within the bacterial and human genomes, respectively. We also developed a deep learning module for the prediction of potential TIS in a given nucleotide sequence. In the future, we aim to optimize the performance of this prediction model through extensive training on larger data sets.

Tailored UPRE2 variants for dynamic gene regulation in yeast

Genetic elements are foundational in synthetic biology serving as vital building blocks. They enable programming host cells for efficient production of valuable chemicals and recombinant proteins. The unfolded protein response (UPR) is a stress pathway in which the transcription factor Hac1 interacts with the upstream unfolded protein response element (UPRE) of the promoter to restore endoplasmic reticulum (ER) homeostasis. Here, we created a UPRE2 mutant (UPRE2m) library. Several rounds of screening identified many elements with enhanced responsiveness and a wider dynamic range. The most active element m84 displayed a response activity 3.72 times higher than the native UPRE2. These potent elements are versatile and compatible with various promoters. Overexpression of HAC1 enhanced stress signal transduction, expanding the signal output range of UPRE2m. Through molecular modeling and site-directed mutagenesis, we pinpointed the DNA-binding residue Lys60 in Hac1(Hac1-K60). We also confirmed that UPRE2m exhibited a higher binding affinity to Hac1. This shed light on the mechanism underlying the Hac1-UPRE2m interaction. Importantly, applying UPRE2m for target gene regulation effectively increased both recombinant protein production and natural product synthesis. These genetic elements provide valuable tools for dynamically regulating gene expression in yeast cell factories.

Towards a rational approach to promoter engineering: understanding the complexity of transcription initiation in prokaryotes

Promoter sequences are important genetic control elements. Through their interaction with RNA polymerase they determine transcription strength and specificity, thereby regulating the first step in gene expression. Consequently, they can be targeted as elements to control predictability and tuneability of a genetic circuit, which is essential in applications such as the development of robust microbial cell factories. This review considers the promoter elements implicated in the three stages of transcription initiation, detailing the complex interplay of sequence-specific interactions that are involved, and highlighting that DNA sequence features beyond the core promoter elements work in a combinatorial manner to determine transcriptional strength. In particular, we emphasize that, aside from promoter recognition, transcription initiation is also defined by the kinetics of open complex formation and promoter escape, which are also known to be highly sequence specific. Significantly, we focus on how insights into these interactions can be manipulated to lay the foundation for a more rational approach to promoter engineering.

Single-molecule visualization of twin-supercoiled domains generated during transcription

Transcription-coupled supercoiling of DNA is a key factor in chromosome compaction and the regulation of genetic processes in all domains of life. It has become common knowledge that, during transcription, the DNA-dependent RNA polymerase (RNAP) induces positive supercoiling ahead of it (downstream) and negative supercoils in its wake (upstream), as rotation of RNAP around the DNA axis upon tracking its helical groove gets constrained due to drag on its RNA transcript. Here, we experimentally validate this so-called twin-supercoiled-domain model with in vitro real-time visualization at the single-molecule scale. Upon binding to the promoter site on a supercoiled DNA molecule, RNAP merges all DNA supercoils into one large pinned plectoneme with RNAP residing at its apex. Transcription by RNAP in real time demonstrates that up- and downstream supercoils are generated simultaneously and in equal portions, in agreement with the twin-supercoiled-domain model. Experiments carried out in the presence of RNases A and H, revealed that an additional viscous drag of the RNA transcript is not necessary for the RNAP to induce supercoils. The latter results contrast the current consensus and simulations on the origin of the twin-supercoiled domains, pointing at an additional mechanistic cause underlying supercoil generation by RNAP in transcription.

Identification and characterization of a strong constitutive promoter stnYp for activating biosynthetic genes and producing natural products in streptomyces

BACKGROUND

Streptomyces are well known for their potential to produce various pharmaceutically active compounds, the commercial development of which is often limited by the low productivity and purity of the desired compounds expressed by natural producers. Well-characterized promoters are crucial for driving the expression of target genes and improving the production of metabolites of interest.

RESULTS

A strong constitutive promoter, stnYp, was identified in Streptomyces flocculus CGMCC4.1223 and was characterized by its effective activation of silent biosynthetic genes and high efficiency of heterologous gene expression. The promoter stnYp showed the highest activity in model strains of four Streptomyces species compared with the three frequently used constitutive promoters ermEp*, kasOp*, and SP44. The promoter stnYp could efficiently activate the indigoidine biosynthetic gene cluster in S. albus J1074, which is thought to be silent under routine laboratory conditions. Moreover, stnYp was found suitable for heterologous gene expression in different Streptomyces hosts. Compared with the promoters ermEp*, kasOp*, and SP44, stnYp conferred the highest production level of diverse metabolites in various heterologous hosts, including the agricultural-bactericide aureonuclemycin and the antitumor compound YM-216391, with an approximately 1.4 - 11.6-fold enhancement of the yields. Furthermore, the purity of tylosin A was greatly improved by overexpressing rate-limiting genes through stnYp in the industrial strain. Further, the yield of tylosin A was significantly elevated to 10.30 ± 0.12 g/L, approximately 1.7-fold higher than that of the original strain.

CONCLUSIONS

The promoter stnYp is a reliable, well-defined promoter with strong activity and broad suitability. The findings of this study can expand promoter diversity, facilitate genetic manipulation, and promote metabolic engineering in multiple Streptomyces species.

Evolutionary Invariant of the Structure of DNA Double Helix in RNAP II Core Promoters

Eukaryotic and archaeal RNA polymerase II (POL II) machinery is highly conserved, regardless of the extreme changes in promoter sequences in different organisms. The goal of our work is to find the cause of this conservatism. The representative sets of aligned promoter sequences of fifteen organisms belonging to different evolutional stages were studied. Their textual profiles, as well as profiles of the indexes that characterize the secondary structure and the mechanical and physicochemical properties, were analyzed. The evolutionarily stable, extremely heterogeneous special secondary structure of POL II core promoters was revealed, which includes two singular regions—hexanucleotide “INR” around TSS and octanucleotide “TATA element” of about −28 bp upstream. Such structures may have developed at some stage of evolution. It turned out to be so well matched for the pre-initiation complex formation and the subsequent initiation of transcription for POL II machinery that in the course of evolution there were selected only those nucleotide sequences that were able to reproduce these structural properties. The individual features of specific sequences representing the singular region of the promoter of each gene can affect the kinetics of DNA-protein complex formation and facilitate strand separation in double-stranded DNA at the TSS position.

Genome-Wide Evolutionary Analysis of Putative Non-Specific Herbicide Resistance Genes and Compilation of Core Promoters between Monocots and Dicots

Herbicides are key weed-control tools, but their repeated use across large areas has favored the evolution of herbicide resistance. Although target-site has been the most prevalent and studied type of resistance, non-target-site resistance (NTSR) is increasing. However, the genetic factors involved in NTSR are widely unknown. In this study, four gene groups encoding putative NTSR enzymes, namely, cytochrome-P450, glutathione-S-transferase (GST), uridine 5'-diphospho-glucuronosyltransferase (UDPGT), and nitronate monooxygenase (NMO) were analyzed. The monocot and dicot gene sequences were downloaded from publicly available databases. Phylogenetic trees revealed that most of the CYP450 resistance-related sequences belong to CYP81 (5), and in GST, most of the resistance sequences belonged to GSTU18 (9) and GSTF6 (8) groups. In addition, the study of upstream promoter sequences of these NTSR genes revealed stress-related cis-regulatory motifs, as well as eight transcription factor binding sites (TFBS) were identified. The discovered TFBS were commonly present in both monocots and dicots, and the identified motifs are known to play key roles in countering abiotic stress. Further, we predicted the 3D structure for the resistant CYP450 and GST protein and identified the substrate recognition site through the homology approach. Our description of putative NTSR enzymes may be used to develop innovative weed control techniques to delay the evolution of NTSR.

Bio-informatic analysis of CRISPR protospacer adjacent motifs (PAMs) in T4 genome

BACKGROUND

The existence of protospacer adjacent motifs (PAMs) sequences in bacteriophage genome is critical for the recognition and function of the clustered regularly interspaced short palindromic repeats-Cas (CRISPR-Cas) machinery system. We further elucidate the significance of PAMs and their function, particularly as a part of transcriptional regulatory regions in T4 bacteriophages.

METHODS

A scripting language was used to analyze a sequence of T4 phage genome, and a list of few selected PAMs. Mann-Whitney Wilcoxon (MWW) test was used to compare the sequence hits for the PAMs versus the hits of all the possible sequences of equal lengths.

RESULTS

The results of MWW test show that certain PAMs such as: 'NGG' and 'TATA' are preferably located at the core of phage promoters: around -10 position, whereas the position around -35 appears to have no detectable count variation of any of the tested PAMs. Among all tested PAMs, the following three sequences: 5'-GCTV-3', 5'-TTGAAT-3' and 5'-TTGGGT-3' have higher prevalence in essential genes. By analyzing all the possible ways of reading PAM sequences as codons for the corresponding amino acids, it was found that deduced amino acids of some PAMs have a significant tendency to prefer the surface of proteins.

CONCLUSION

These results provide novel insights into the location and the subsequent identification of the role of PAMs as transcriptional regulatory elements. Also, CRISPR targeting certain PAM sequences is somehow likely to be connected to the hydrophilicity (water solubility) of amino acids translated from PAM's triplets. Therefore, these amino acids are found at the interacting unit at protein-protein interfaces.

Machine learning and statistics shape a novel path in archaeal promoter annotation

Background

Archaea are a vast and unexplored domain. Bioinformatic techniques might enlighten the path to a higher quality genome annotation in varied organisms. Promoter sequences of archaea have the action of a plethora of proteins upon it. The conservation found in a structural level of the binding site of proteins such as TBP, TFB, and TFE aids RNAP-DNA stabilization and makes the archaeal promoter prone to be explored by statistical and machine learning techniques.

Results and discussions

In this study, experimentally verified promoter sequences of the organisms Haloferax volcanii, Sulfolobus solfataricus, and Thermococcus kodakarensis were converted into DNA duplex stability attributes (i.e. numerical variables) and were classified through Artificial Neural Networks and an in-house statistical method of classification, being tested with three forms of controls. The recognition of these promoters enabled its use to validate unannotated promoter sequences in other organisms. As a result, the binding site of basal transcription factors was located through a DNA duplex stability codification. Additionally, the classification presented satisfactory results (above 90%) among varied levels of control.

Concluding remarks

The classification models were employed to perform genomic annotation into the archaea Aciduliprofundum boonei and Thermofilum pendens, from which potential promoters have been identified and uploaded into public repositories.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04714-x.

Does the Pachytene Checkpoint, a Feature of Meiosis, Filter Out Mistakes in Double-Strand DNA Break Repair and as a side-Effect Strongly Promote Adaptive Speciation?

This essay aims to explain two biological puzzles: why eukaryotic transcription units are composed of short segments of coding DNA interspersed with long stretches of non-coding (intron) DNA, and the near ubiquity of sexual reproduction. As is well known, alternative splicing of its coding sequences enables one transcription unit to produce multiple variants of each encoded protein. Additionally, padding transcription units with non-coding DNA (often many thousands of base pairs long) provides a readily evolvable way to set how soon in a cell cycle the various mRNAs will begin being expressed and the total amount of mRNA that each transcription unit can make during a cell cycle. This regulation complements control via the transcriptional promoter and facilitates the creation of complex eukaryotic cell types, tissues, and organisms. However, it also makes eukaryotes exceedingly vulnerable to double-strand DNA breaks, which end-joining break repair pathways can repair incorrectly. Transcription units cover such a large fraction of the genome that any mis-repair producing a reorganized chromosome has a high probability of destroying a gene. During meiosis, the synaptonemal complex aligns homologous chromosome pairs and the pachytene checkpoint detects, selectively arrests, and in many organisms actively destroys gamete-producing cells with chromosomes that cannot adequately synapse; this creates a filter favoring transmission to the next generation of chromosomes that retain the parental organization, while selectively culling those with interrupted transcription units. This same meiotic checkpoint, reacting to accidental chromosomal reorganizations inflicted by error-prone break repair, can, as a side effect, provide a mechanism for the formation of new species in sympatry. It has been a long-standing puzzle how something as seemingly maladaptive as hybrid sterility between such new species can arise. I suggest that this paradox is resolved by understanding the adaptive importance of the pachytene checkpoint, as outlined above.

Programmable Synthetic Upstream Activating Sequence Library for Fine-Tuning Gene Expression Levels in Saccharomyces cerevisiae

A wide dynamic range of promoters is necessary for fine-tuning transcription levels. However, weak intensity and narrow dynamic range limit transcriptional regulation via constitutive promoters. The upstream activation sequence (UAS) located upstream of the core promoter is a crucial region that could obviously enhance promoter strength. Herein, we created a random mutagenesis library consisting of 330 different variants based on the UAS of the TDH3 promoter with an ∼37-fold dynamic range by error-prone polymerase chain reaction (PCR) and obtained strong intensity mutant UAS, which was ∼12-fold greater than the wild-type UAS_TDH3. Analysis of the mutant library revealed 15 strength-enhancing sites and their corresponding bases of the UAS_TDH3 regions, which provided the impetus for a synthetic library. The resulting 32 768 mutant UAS library was constructed by permutation and combination of the bases of the 15 enhancing sites. To characterize the library, a strength prediction model was built by correlating DNA structural features and UAS strength, which provided a model between UAS sequence and intensity. Following characterization, the UAS library was applied to precisely regulate gene expression in the production of β-carotene, proving that the UAS library would be a useful tool for gene tuning in metabolic engineering. In summary, we designed, constructed, and characterized a UAS library that facilitated precise tuning of transcription levels of target proteins.

Advancing Yarrowia lipolytica as a superior biomanufacturing platform by tuning gene expression using promoter engineering

Yarrowia lipolytica is recognized as an excellent non-conventional yeast in the field of biomanufacturing, where it is used as a host to produce oleochemicals, terpenes, organic acids, polyols and recombinant proteins. Consequently, metabolic engineering of this yeast is becoming increasingly popular to advance it as a superior biomanufacturing platform, of which promoters are the most basic elements for tuning gene expression. Endogenous promoters of Yarrowia lipolytica were reviewed, which are the basis for promoter engineering. The engineering strategies, such as hybrid promoter engineering, intron enhancement promoter engineering, and transcription factor-based inducible promoter engineering are described. Additionally, the applications of Yarrowia lipolytica promoter engineering to rationally reconstruct biosynthetic gene clusters and improve the genome-editing efficiency of the CRISPR-Cas systems were reviewed. Finally, research needs and future directions for promoter engineering are also discussed in this review.

What Have We Learned in 30 Years of Investigations on Bari Transposons?

Transposable elements (TEs) have been historically depicted as detrimental genetic entities that selfishly aim at perpetuating themselves, invading genomes, and destroying genes. Scientists often co-opt "special" TEs to develop new and powerful genetic tools, that will hopefully aid in changing the future of the human being. However, many TEs are gentle, rarely unleash themselves to harm the genome, and bashfully contribute to generating diversity and novelty in the genomes they have colonized, yet they offer the opportunity to develop new molecular tools. In this review we summarize 30 years of research focused on the Bari transposons. Bari is a "normal" transposon family that has colonized the genomes of several Drosophila species and introduced genomic novelties in the melanogaster species. We discuss how these results have contributed to advance the field of TE research and what future studies can still add to the current knowledge.

Genome-Wide Prediction of Transcription Start Sites in Conifers

The identification of promoters is an essential step in the genome annotation process, providing a framework for gene regulatory networks and their role in transcription regulation. Despite considerable advances in the high-throughput determination of transcription start sites (TSSs) and transcription factor binding sites (TFBSs), experimental methods are still time-consuming and expensive. Instead, several computational approaches have been developed to provide fast and reliable means for predicting the location of TSSs and regulatory motifs on a genome-wide scale. Numerous studies have been carried out on the regulatory elements of mammalian genomes, but plant promoters, especially in gymnosperms, have been left out of the limelight and, therefore, have been poorly investigated. The aim of this study was to enhance and expand the existing genome annotations using computational approaches for genome-wide prediction of TSSs in the four conifer species: loblolly pine, white spruce, Norway spruce, and Siberian larch. Our pipeline will be useful for TSS predictions in other genomes, especially for draft assemblies, where reliable TSS predictions are not usually available. We also explored some of the features of the nucleotide composition of the predicted promoters and compared the GC properties of conifer genes with model monocot and dicot plants. Here, we demonstrate that even incomplete genome assemblies and partial annotations can be a reliable starting point for TSS annotation. The results of the TSS prediction in four conifer species have been deposited in the Persephone genome browser, which allows smooth visualization and is optimized for large data sets. This work provides the initial basis for future experimental validation and the study of the regulatory regions to understand gene regulation in gymnosperms.

Analysis of nucleoid-associated protein-binding regions reveals DNA structural features influencing genome organization in Mycobacterium tuberculosis

Nucleoid-associated proteins (NAPs) maintain bacterial nucleoid configuration through their architectural properties of DNA bending, wrapping, and bridging. However, the contribution of DNA structural alterations to DNA-NAP recognition at the genomic scale remains unresolved. Present work dissects the DNA sequence, shape and altered structural preferences at a genomic scale for six NAPs in Mycobacterium tuberculosis. Results suggest narrower minor groove width (MGW) and higher DNA rigidity are marked for the binding sites of EspR and Lsr2, while mIHF, MtHU and NapM have heterogeneous DNA structural predilections. In contrast, WhiB4-DNA-binding sites were characterized by wider MGW, highly deformable and less curved DNA. This work provides systematic insight into NAP-mediated genome organization as a function of DNA structural features.

Low-level expression of the Type II restriction-modification system confers potent bacteriophage resistance in Escherichia coli

Restriction–modification systems (R–M) are one of the antiviral defense tools used by bacteria, and those of the Type II family are composed of a restriction endonuclease (REase) and a DNA methyltransferase (MTase). Most entering DNA molecules are usually cleaved by the REase before they can be methylated by MTase, although the observed level of fragmented DNA may vary significantly. Using a model EcoRI R–M system, we report that the balance between DNA methylation and cleavage may be severely affected by transcriptional signals coming from outside the R–M operon. By modulating the activity of the promoter, we obtained a broad range of restriction phenotypes for the EcoRI R–M system that differed by up to 4 orders of magnitude in our biological assays. Surprisingly, we found that high expression levels of the R–M proteins were associated with reduced restriction of invading bacteriophage DNA. Our results suggested that the regulatory balance of cleavage and methylation was highly sensitive to fluctuations in transcriptional signals both up- and downstream of the R–M operon. Our data provided further insights into Type II R–M system maintenance and the potential conflict within the host bacterium.

A supernumerary designer chromosome for modular in vivo pathway assembly in Saccharomyces cerevisiae

The construction of microbial cell factories for sustainable production of chemicals and pharmaceuticals requires extensive genome engineering. Using Saccharomyces cerevisiae, this study proposes synthetic neochromosomes as orthogonal expression platforms for rewiring native cellular processes and implementing new functionalities. Capitalizing the powerful homologous recombination capability of S. cerevisiae, modular neochromosomes of 50 and 100 kb were fully assembled de novo from up to 44 transcriptional-unit-sized fragments in a single transformation. These assemblies were remarkably efficient and faithful to their in silico design. Neochromosomes made of non-coding DNA were stably replicated and segregated irrespective of their size without affecting the physiology of their host. These non-coding neochromosomes were successfully used as landing pad and as exclusive expression platform for the essential glycolytic pathway. This work pushes the limit of DNA assembly in S. cerevisiae and paves the way for de novo designer chromosomes as modular genome engineering platforms in S. cerevisiae.

DIANA-miRGen v4: indexing promoters and regulators for more than 1500 microRNAs

Deregulation of microRNA (miRNA) expression plays a critical role in the transition from a physiological to a pathological state. The accurate miRNA promoter identification in multiple cell types is a fundamental endeavor towards understanding and characterizing the underlying mechanisms of both physiological as well as pathological conditions. DIANA-miRGen v4 (www.microrna.gr/mirgenv4) provides cell type specific miRNA transcription start sites (TSSs) for over 1500 miRNAs retrieved from the analysis of >1000 cap analysis of gene expression (CAGE) samples corresponding to 133 tissues, cell lines and primary cells available in FANTOM repository. MiRNA TSS locations were associated with transcription factor binding site (TFBSs) annotation, for >280 TFs, derived from analyzing the majority of ENCODE ChIP-Seq datasets. For the first time, clusters of cell types having common miRNA TSSs are characterized and provided through a user friendly interface with multiple layers of customization. DIANA-miRGen v4 significantly improves our understanding of miRNA biogenesis regulation at the transcriptional level by providing a unique integration of high-quality annotations for hundreds of cell specific miRNA promoters with experimentally derived TFBSs.

Modeling SELEX for regulatory regions using Royal Road and Royal Staircase fitness functions

The field of evolutionary algorithms (EAs) emerged in the area of computer science due to transfer of ideas from biology and developed independently for several decades, enriched with techniques from probability theory, complexity theory and optimization methods. In this paper, we consider some recent results form the EAs theory transferred back into biology. The well-known biotechnological procedure SELEX (Systematic Evolution of Ligands by EXponential enrichment) is viewed as an experimental implementation of an evolutionary algorithm. Theoretical bounds on EAs runtime are applied to model SELEX search for a regulatory region consisting of promoter and enhancer sequences. A comparison of theoretical bounds to the results of computational simulation indicates some cases where the theoretical bounds give favorable prediction, while simulation requires prohibitive computational resource.

Delving into Eukaryotic Origins of Replication Using DNA Structural Features

DNA replication in eukaryotes is an intricate process, which is precisely synchronized by a set of regulatory proteins, and the replication fork emanates from discrete sites on chromatin called origins of replication (Oris). These spots are considered as the gateway to chromosomal replication and are stereotyped by sequence motifs. The cognate sequences are noticeable in a small group of entire origin regions or totally absent across different metazoans. Alternatively, the use of DNA secondary structural features can provide additional information compared to the primary sequence. In this article, we report the trends in DNA sequence-based structural properties of origin sequences in nine eukaryotic systems representing different families of life. Biologically relevant DNA secondary structural properties, namely, stability, propeller twist, flexibility, and minor groove shape were studied in the sequences flanking replication start sites. Results indicate that Oris in yeasts show lower stability, more rigidity, and narrow minor groove preferences compared to genomic sequences surrounding them. Yeast Oris also show preference for A-tracts and the promoter element TATA box in the vicinity of replication start sites. On the contrary, Drosophila melanogaster, humans, and Arabidopsis thaliana do not have such features in their Oris, and instead, they show high preponderance of G-rich sequence motifs such as putative G-quadruplexes or i-motifs and CpG islands. Our extensive study applies the DNA structural feature computation to delve into origins of replication across organisms ranging from yeasts to mammals and including a plant. Insights from this study would be significant in understanding origin architecture and help in designing new algorithms for predicting DNA trans-acting factor recognition events.

Antibiotic drug discovery: Challenges and perspectives in the light of emerging antibiotic resistance

Amid a rising threat of antimicrobial resistance in a global scenario, our huge investments and high-throughput technologies injected for rejuvenating the key therapeutic scaffolds to suppress these rising superbugs has been diminishing severely. This has grasped world-wide attention, with increased consideration being given to the discovery of new chemical entities. Research has now proven that the relatively tiny and simpler microbes possess enhanced capability of generating novel and diverse chemical constituents with huge therapeutic leads. The usage of these beneficial organisms could help in producing new chemical scaffolds that govern the power to suppress the spread of obnoxious superbugs. Here in this review, we have explicitly focused on several appealing strategies employed for the generation of new chemical scaffolds. Also, efforts on providing novel insights on some of the unresolved questions in the production of metabolites, metabolic profiling and also the serendipity of getting "hit molecules" have been rigorously discussed. However, we are highly aware that biosynthetic pathway of different classes of secondary metabolites and their biosynthetic route is a vast topic, thus we have avoided discussion on this topic.

Solving the transcription start site identification problem with ADAPT-CAGE: a Machine Learning algorithm for the analysis of CAGE data

Cap Analysis of Gene Expression (CAGE) has emerged as a powerful experimental technique for assisting in the identification of transcription start sites (TSSs). There is strong evidence that CAGE also identifies capping sites along various other locations of transcribed loci such as splicing byproducts, alternative isoforms and capped molecules overlapping introns and exons. We present ADAPT-CAGE, a Machine Learning framework which is trained to distinguish between CAGE signal derived from TSSs and transcriptional noise. ADAPT-CAGE provides highly accurate experimentally derived TSSs on a genome-wide scale. It has been specifically designed for flexibility and ease-of-use by only requiring aligned CAGE data and the underlying genomic sequence. When compared to existing algorithms, ADAPT-CAGE exhibits improved performance on every benchmark that we designed based on both annotation- and experimentally-driven strategies. This performance boost brings ADAPT-CAGE in the spotlight as a computational framework that is able to assist in the refinement of gene regulatory networks, the incorporation of accurate information of gene expression regulators and alternative promoter usage in both physiological and pathological conditions.

Variation of gene expression in plants is influenced by gene architecture and structural properties of promoters

In higher eukaryotes, gene architecture and structural properties of promoters have emerged as significant factors influencing variation in number of transcripts (expression level) and specificity of gene expression in a tissue (expression breadth), which eventually shape the phenotype. In this study, transcriptome data of different tissue types at various developmental stages of A. thaliana, O. sativa, S. bicolor and Z. mays have been used to understand the relationship between properties of gene components and its expression. Our findings indicate that in plants, among all gene architecture and structural properties of promoters, compactness of genes in terms of intron content is significantly linked to gene expression level and breadth, whereas in human an exactly opposite scenario is seen. In plants, for the first time we have carried out a quantitative estimation of effect of a particular trait on expression level and breadth, by using multiple regression analysis and it confirms that intron content of primary transcript (as %) is a powerful determinant of expression breadth. Similarly, further regression analysis revealed that among structural properties of the promoters, stability is negatively linked to expression breadth, while DNase1 sensitivity strongly governs gene expression breadth in monocots and gene expression level in dicots. In addition, promoter regions of tissue specific genes are found to be enriched with TATA box and Y-patch motifs. Finally, multi copy orthologous genes in plants are found to be longer, highly regulated and tissue specific.

Truncated, strong inducible promoter Pmcl1 from Metarhizium anisopliae

In this study, Metarhizium collagen -like protein (MCL1) promoter from Metarhizium anisopliae was analysed and truncated into different sizes through series of targeted and random deletions based on the presence of various transcription factor-binding sites. Synthetic Green Fluorescent Protein (sGFP) was being utilized as a reporter gene to study the relative expression driving capability of unmodified and truncated promoters. Conserved promoter sequence analysis revealed similarity between the paralogous promoters from M. brunneum and M. acridum. sGFP expression in the haemolymph was directed with the help of mcl1 signal peptide sequence. Deleting the promoter region from - 2764 to - 1583 bp increases the promoter mcl1 (Pmcl1) activity by twofolds, while deletions of the regions upstream of - 1150 bp and - 840 bp caused a decrease of sGFP expression level (80% and 70%, respectively). Transcriptional binding sites predicted for the deleted region revealed the loss of upstream repressing sequences such as Matalpha2 along with ROX1 and Rap1 repressor-binding sites located - 2234 bp, - 1754 bp and - 1724 bp from the TSS. Compared with Pmcl1-wild type (2.7 kbp), Pmcl1-1583 bp had a shorter sequence and showed statistically significant expression in M. anisopliae. This study introduces a highly efficient strong inducible promoter for over-expression of target genes in M. anisopliae.

DNA sequence encodes the position of DNA supercoils

The three-dimensional organization of DNA is increasingly understood to play a decisive role in vital cellular processes. Many studies focus on the role of DNA-packaging proteins, crowding, and confinement in arranging chromatin, but structural information might also be directly encoded in bare DNA itself. Here, we visualize plectonemes (extended intertwined DNA structures formed upon supercoiling) on individual DNA molecules. Remarkably, our experiments show that the DNA sequence directly encodes the structure of supercoiled DNA by pinning plectonemes at specific sequences. We develop a physical model that predicts that sequence-dependent intrinsic curvature is the key determinant of pinning strength and demonstrate this simple model provides very good agreement with the data. Analysis of several prokaryotic genomes indicates that plectonemes localize directly upstream of promoters, which we experimentally confirm for selected promotor sequences. Our findings reveal a hidden code in the genome that helps to spatially organize the chromosomal DNA.

ShapeGTB: the role of local DNA shape in prioritization of functional variants in human promoters with machine learning

Motivation

The identification of functional sequence variations in regulatory DNA regions is one of the major challenges of modern genetics. Here, we report results of a combined multifactor analysis of properties characterizing functional sequence variants located in promoter regions of genes.

Results

We demonstrate that GC-content of the local sequence fragments and local DNA shape features play significant role in prioritization of functional variants and outscore features related to histone modifications, transcription factors binding sites, or evolutionary conservation descriptors. Those observations allowed us to build specialized machine learning classifier identifying functional single nucleotide polymorphisms within promoter regions—ShapeGTB. We compared our method with more general tools predicting pathogenicity of all non-coding variants. ShapeGTB outperformed them by a wide margin (average precision 0.93 vs. 0.47–0.55). On the external validation set based on ClinVar database it displayed worse performance but was still competitive with other methods (average precision 0.47 vs. 0.23–0.42). Such results suggest unique characteristics of mutations located within promoter regions and are a promising signal for the development of more accurate variant prioritization tools in the future.

Ambiguity in logic-based models of gene regulatory networks: An integrative multi-perturbation analysis

Most studies of gene regulatory network (GRN) inference have focused extensively on identifying the interaction map of the GRNs. However, in order to predict the cellular behavior, modeling the GRN in terms of logic circuits, i.e., Boolean networks, is necessary. The perturbation techniques, e.g., knock-down and over-expression, should be utilized for identifying the underlying logic behind the interactions. However, we will show that by using only transcriptomic data obtained by single-perturbation experiments, we cannot observe all regulatory interactions, and this invisibility causes ambiguity in our model. Consequently, we need to employ the data of multiple omics layers (genome, transcriptome, and proteome) as well as multiple perturbation experiments to reduce or eliminate ambiguity in our modeling. In this paper, we introduce a multi-step perturbation experiment to deal with ambiguity. Moreover, we perform a thorough analysis to investigate which types of perturbations and omics layers play the most important role in the unambiguous modeling of the GRNs and how much ambiguity will be eliminated by considering more perturbations and more omics layers. Our analysis shows that performing both knock-down and over-expression is necessary in order to achieve the least ambiguous model. Moreover, the more steps of the perturbation are taken, the more ambiguity is eliminated. In addition, we can even achieve an unambiguous model of the GRN by using multi-step perturbation and integrating transcriptomic, protein-protein interaction, and cis-element data. Finally, we demonstrate the effect of utilizing different types of perturbation experiment and integrating multi-omics data on identifying the logic behind the regulatory interactions in a synthetic GRN. In conclusion, relying on the results of only knock-down experiments and not including as many omics layers as possible in the GRN inference, makes the results ambiguous, unreliable, and less accurate.

RAIL: A new tool for defining bacterial promoter regions

In this issue of the Journal of Bacteriology, Hustmyer and colleagues describe a new method for rapidly generating reporter libraries (Hustmyer citation). This RAIL technique (Rapid Arbitrary PCR Insertion Libraries) uses arbitrary PCR and isothermal DNA assembly to insert random fragments of promoter regions into reporter plasmids, resulting in libraries that can be screened to identify regions required for gene expression. This technique will likely be useful for a number of different genetic applications.

Synthetic promoter design for new microbial chassis

The judicious choice of promoter to drive gene expression remains one of the most important considerations for synthetic biology applications. Constitutive promoter sequences isolated from nature are often used in laboratory settings or small-scale commercial production streams, but unconventional microbial chassis for new synthetic biology applications require well-characterized, robust and orthogonal promoters. This review provides an overview of the opportunities and challenges for synthetic promoter discovery and design, including molecular methodologies, such as saturation mutagenesis of flanking regions and mutagenesis by error-prone PCR, as well as the less familiar use of computational and statistical analyses for de novo promoter design.

Pgas, a Low-pH-Induced Promoter, as a Tool for Dynamic Control of Gene Expression for Metabolic Engineering of Aspergillus niger

The dynamic control of gene expression is important for adjusting fluxes in order to obtain desired products and achieve appropriate cell growth, particularly when the synthesis of a desired product drains metabolites required for cell growth. For dynamic gene expression, a promoter responsive to a particular environmental stressor is vital. Here, we report a low-pH-inducible promoter, Pgas, which promotes minimal gene expression at pH values above 5.0 but functions efficiently at low pHs, such as pH 2.0. First, we performed a transcriptional analysis of Aspergillus niger, an excellent platform for the production of organic acids, and we found that the promoter Pgas may act efficiently at low pH. Then, a gene for synthetic green fluorescent protein (sGFP) was successfully expressed by Pgas at pH 2.0, verifying the results of the transcriptional analysis. Next, Pgas was used to express the cis-aconitate decarboxylase (cad) gene of Aspergillus terreus in A. niger, allowing the production of itaconic acid at a titer of 4.92 g/liter. Finally, we found that Pgas strength was independent of acid type and acid ion concentration, showing dependence on pH only.IMPORTANCE The promoter Pgas can be used for the dynamic control of gene expression in A. niger for metabolic engineering to produce organic acids. This promoter may also be a candidate tool for genetic engineering.

DNA structural features of eukaryotic TATA-containing and TATA-less promoters

Eukaryotic genes can be broadly classified as TATA‐containing and TATA‐less based on the presence of TATA box in their promoters. Experiments on both classes of genes have revealed a disparity in the regulation of gene expression and cellular functions between the two classes. In this study, we report characteristic differences in promoter sequences and associated structural properties of the two categories of genes in six different eukaryotes. We have analyzed three structural features, DNA duplex stability, bendability, and curvature along with the distribution of A‐tracts, G‐quadruplex motifs, and CpG islands. The structural feature analyses reveal that while the two classes of gene promoters are distinctly different from each other, the properties are also distinguishable across the six organisms.

Unveiling DNA structural features of promoters associated with various types of TSSs in prokaryotic transcriptomes and their role in gene expression

Next-generation sequencing studies have revealed that a variety of transcripts are present in the prokaryotic transcriptome and a significant fraction of them are functional, being involved in various regulatory activities apart from coding for proteins. Identification of promoters associated with different transcripts is necessary for characterization of the transcriptome. Promoter regions have been shown to have unique structural features as compared with their flanking region, in organisms covering all domains of life. Here we report an in silico analysis of DNA sequence dependent structural properties like stability, bendability and curvature in the promoter region of six different prokaryotic transcriptomes. Using these structural features, we predicted promoters associated with different categories of transcripts (mRNA, internal, antisense and non-coding), which constitute the transcriptome. Promoter annotation using structural features is fairly accurate and reliable with about 50% of the primary promoters being characterized by all three structural properties while at least one property identifies 95%. We also studied the relative differences of these structural features in terms of gene expression and found that the features, viz. lower stability, lesser bendability and higher curvature are more prominent in the promoter regions which are associated with high gene expression as compared with low expression genes. Hence, promoters, which are associated with higher gene expression, get annotated well using DNA structural features as compared with those, which are linked to lower gene expression.

Structural features of DNA that determine RNA polymerase II core promoter

Background

The general structure and action of all eukaryotic and archaeal RNA polymerases machinery have an astonishing similarity despite the diversity of core promoter sequences in different species. The goal of our work is to find common characteristics of DNA region that define it as a promoter for the RNA polymerase II (Pol II).

Results

The profiles of a large number of physical and structural characteristics, averaged over representative sets of the Pol II minimal core promoters of the evolutionary divergent species from animals, plants and unicellular fungi were analysed. In addition to the characteristics defined at the base-pair steps, we, for the first time, use profiles of the ultrasonic cleavage and DNase I cleavage indexes, informative for internal properties of each complementary strand.

Conclusions

DNA of the core promoters of metazoans and Schizosaccharomyces pombe has similar structural organization. Its mechanical and 3D structural characteristics have singular properties at the positions of TATA-box. The minor groove is broadened and conformational motion is decreased in that region. Special characteristics of conformational behavior are revealed in metazoans at the region, which connects the end of TATA-box and the transcription start site (TSS). The intensities of conformational motions in the complementary strands are periodically changed in opposite phases. They are noticeable, best of all, in mammals. Such conformational features are lacking in the core promoters of S. pombe. The profiles of Saccharomyces cerevisiae core promoters significantly differ: their singular region is shifted down thus pointing to the uniqueness of their structural organization. Obtained results may be useful in genetic engineering for artificial modulation of the promoter strength.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3292-z) contains supplementary material, which is available to authorized users.

MethSMRT: an integrative database for DNA N6-methyladenine and N4-methylcytosine generated by single-molecular real-time sequencing

DNA methylation is an important type of epigenetic modifications, where 5- methylcytosine (5mC), 6-methyadenine (6mA) and 4-methylcytosine (4mC) are the most common types. Previous efforts have been largely focused on 5mC, providing invaluable insights into epigenetic regulation through DNA methylation. Recently developed single-molecule real-time (SMRT) sequencing technology provides a unique opportunity to detect the less studied DNA 6mA and 4mC modifications at single-nucleotide resolution. With a rapidly increased amount of SMRT sequencing data generated, there is an emerging demand to systematically explore DNA 6mA and 4mC modifications from these data sets. MethSMRT is the first resource hosting DNA 6mA and 4mC methylomes. All the data sets were processed using the same analysis pipeline with the same quality control. The current version of the database provides a platform to store, browse, search and download epigenome-wide methylation profiles of 156 species, including seven eukaryotes such as Arabidopsis, C. elegans, Drosophila, mouse and yeast, as well as 149 prokaryotes. It also offers a genome browser to visualize the methylation sites and related information such as single nucleotide polymorphisms (SNP) and genomic annotation. Furthermore, the database provides a quick summary of statistics of methylome of 6mA and 4mC and predicted methylation motifs for each species. MethSMRT is publicly available at http://sysbio.sysu.edu.cn/methsmrt/ without use restriction.

Genome-Wide Targets Regulated by the OsMADS1 Transcription Factor Reveals Its DNA Recognition Properties

OsMADS1 controls rice (Oryza sativa) floral fate and organ development. Yet, its genome-wide targets and the mechanisms underlying its role as a transcription regulator controlling developmental gene expression are unknown. We identify 3112 gene-associated OsMADS1-bound sites in the floret genome. These occur in the vicinity of transcription start sites, within gene bodies, and in intergenic regions. Majority of the bound DNA contained CArG motif variants or, in several cases, only A-tracts. Sequences flanking the binding peak had a higher AT nucleotide content, implying that broader DNA structural features may define in planta binding. Sequences for binding by other transcription factor families like MYC, AP2/ERF, bZIP, etc. are enriched in OsMADS1-bound DNAs. Target genes implicated in transcription, chromatin remodeling, cellular processes, and hormone metabolism were enriched. Combining expression data from OsMADS1 knockdown florets with these DNA binding data, a snapshot of a gene regulatory network was deduced where targets, such as AP2/ERF and bHLH transcription factors and chromatin remodelers form nodes. We show that the expression status of these nodal factors can be altered by inducing the OsMADS1-GR fusion protein and present a model for a regulatory cascade where the direct targets of OsMADS1, OsbHLH108/SPT, OsERF034, and OsHSF24, in turn control genes such as OsMADS32 and OsYABBY5 This cascade, with other similar relationships, cumulatively contributes to floral organ development. Overall, OsMADS1 binds to several regulatory genes and, probably in combination with other factors, controls a gene regulatory network that ensures rice floret development.

Bioinformatic Identification of Conserved Cis-Sequences in Coregulated Genes

Bioinformatics tools can be employed to identify conserved cis-sequences in sets of coregulated plant genes because more and more gene expression and genomic sequence data become available. Knowledge on the specific cis-sequences, their enrichment and arrangement within promoters, facilitates the design of functional synthetic plant promoters that are responsive to specific stresses. The present chapter illustrates an example for the bioinformatic identification of conserved Arabidopsis thaliana cis-sequences enriched in drought stress-responsive genes. This workflow can be applied for the identification of cis-sequences in any sets of coregulated genes. The workflow includes detailed protocols to determine sets of coregulated genes, to extract the corresponding promoter sequences, and how to install and run a software package to identify overrepresented motifs. Further bioinformatic analyses that can be performed with the results are discussed.

Diversity, expansion, and evolutionary novelty of plant DNA-binding transcription factor families

Plant transcription factors (TFs) that interact with specific sequences via DNA-binding domains are crucial for regulating transcriptional initiation and are fundamental to plant development and environmental response. In addition, expansion of TF families has allowed functional divergence of duplicate copies, which has contributed to novel, and in some cases adaptive, traits in plants. Thus, TFs are central to the generation of the diverse plant species that we see today. Major plant agronomic traits, including those relevant to domestication, have also frequently arisen through changes in TF coding sequence or expression patterns. Here our goal is to provide an overview of plant TF evolution by first comparing the diversity of DNA-binding domains and the sizes of these domain families in plants and other eukaryotes. Because TFs are among the most highly expanded gene families in plants, the birth and death process of TFs as well as the mechanisms contributing to their retention are discussed. We also provide recent examples of how TFs have contributed to novel traits that are important in plant evolution and in agriculture.This article is part of a Special Issue entitled: Plant Gene Regulatory Mechanisms and Networks, edited by Dr. Erich Grotewold and Dr. Nathan Springer.

Blue light-mediated transcriptional activation and repression of gene expression in bacteria

Light-regulated modules offer unprecedented new ways to control cellular behavior in precise spatial and temporal resolution. The availability of such tools may dramatically accelerate the progression of synthetic biology applications. Nonetheless, current optogenetic toolbox of prokaryotes has potential issues such as lack of rapid and switchable control, less portable, low dynamic expression and limited parts. To address these shortcomings, we have engineered a novel bidirectional promoter system for Escherichia coli that can be induced or repressed rapidly and reversibly using the blue light dependent DNA-binding protein EL222. We demonstrated that by modulating the dosage of light pulses or intensity we could control the level of gene expression precisely. We show that both light-inducible and repressible system can function in parallel with high spatial precision in a single cell and can be switched stably between ON- and OFF-states by repetitive pulses of blue light. In addition, the light-inducible and repressible expression kinetics were quantitatively analysed using a mathematical model. We further apply the system, for the first time, to optogenetically synchronize two receiver cells performing different logic behaviors over time using blue light as a molecular clock signal. Overall, our modular approach layers a transformative platform for next-generation light-controllable synthetic biology systems in prokaryotes.

Fast Prediction of DNA Melting Bubbles Using DNA Thermodynamic Stability

DNA melting bubbles are the basis of many DNA-protein interactions, such as those in regulatory DNA regions driving gene expression, DNA replication and bacterial horizontal gene transfer. Bubble formation is affected by DNA duplex stability and thermally induced duplex destabilization (TIDD). Although prediction of duplex stability with the nearest neighbor (NN) method is much faster than prediction of TIDD with the Peyrard-Bishop-Dauxois (PBD) model, PBD predicted TIDD defines regulatory DNA regions with higher accuracy and detail. Here, we considered that PBD predicted TIDD is inherently related to the intrinsic duplex stabilities of destabilization regions. We show by regression modeling that NN duplex stabilities can be used to predict TIDD almost as accurately as is predicted with PBD. Predicted TIDD is in fact ascribed to non-linear transformation of NN duplex stabilities in destabilization regions as well as effects of neighboring regions relative to destabilization size. Since the prediction time of our models is over six orders of magnitude shorter than that of PBD, the models present an accessible tool for researchers. TIDD can be predicted on our webserver at http://tidd.immt.eu.

DNA topoisomerases beyond the standard role

Chromatin is dynamically changing its structure to accommodate and control DNA-dependent processes inside of eukaryotic cells. These changes are necessarily linked to changes of DNA topology, which might itself serve as a regulatory signal to be detected by proteins. Thus, DNA Topoisomerases may contribute to the regulation of many events occurring during the transcription cycle. In this review we will focus on DNA Topoisomerase functions in transcription, with particular emphasis on the multiplicity of tasks beyond their widely appreciated role in solving topological problems associated with transcription elongation.