Transcriptional coupling of neighboring genes and gene expression noise: evidence that gene orientation and noncoding transcripts are modulators of noise

Molecular models of bidirectional promoter regulation

Approximately 11% of human genes are transcribed by a bidirectional promoter (BDP), defined as two genes with <1 kb between their transcription start sites. Despite their evolutionary conservation and enrichment for housekeeping genes and oncogenes, the regulatory role of BDPs remains unclear. BDPs have been suggested to facilitate gene coregulation and/or decrease expression noise. This review discusses these potential regulatory functions through the context of six prospective underlying mechanistic models: a single nucleosome free region, shared transcription factor/regulator binding, cooperative negative supercoiling, bimodal histone marks, joint activation by enhancer(s), and RNA-mediated recruitment of regulators. These molecular mechanisms may act independently and/or cooperatively to facilitate the coregulation and/or decreased expression noise predicted of BDPs.

Cellular control of protein levels: A systems biology perspective

How cells regulate protein levels is a central question of biology. Over the past decades, molecular biology research has provided profound insights into the mechanisms and the molecular machinery governing each step of the gene expression process, from transcription to protein degradation. Recent advances in transcriptomics and proteomics have complemented our understanding of these fundamental cellular processes with a quantitative, systems-level perspective. Multi-omic studies revealed significant quantitative, kinetic and functional differences between the genome, transcriptome and proteome. While protein levels often correlate with mRNA levels, quantitative investigations have demonstrated a substantial impact of translation and protein degradation on protein expression control. In addition, protein-level regulation appears to play a crucial role in buffering protein abundances against undesirable mRNA expression variation. These findings have practical implications for many fields, including gene function prediction and precision medicine.

DNA supercoiling restricts the transcriptional bursting of neighboring eukaryotic genes

DNA supercoiling has emerged as a major contributor to gene regulation in bacteria, but how DNA supercoiling impacts transcription dynamics in eukaryotes is unclear. Here, using single-molecule dual-color nascent transcription imaging in budding yeast, we show that transcriptional bursting of divergent and tandem GAL genes is coupled. Temporal coupling of neighboring genes requires rapid release of DNA supercoils by topoisomerases. When DNA supercoils accumulate, transcription of one gene inhibits transcription at its adjacent genes. Transcription inhibition of the GAL genes results from destabilized binding of the transcription factor Gal4. Moreover, wild-type yeast minimizes supercoiling-mediated inhibition by maintaining sufficient levels of topoisomerases. Overall, we discover fundamental differences in transcriptional control by DNA supercoiling between bacteria and yeast and show that rapid supercoiling release in eukaryotes ensures proper gene expression of neighboring genes.

Higher-order modular regulation of the human proteome

Operons are transcriptional modules that allow bacteria to adapt to environmental changes by coordinately expressing the relevant set of genes. In humans, biological pathways and their regulation are more complex. If and how human cells coordinate the expression of entire biological processes is unclear. Here, we capture 31 higher-order co-regulation modules, which we term progulons, by help of supervised machine-learning on proteomics data. Progulons consist of dozens to hundreds of proteins that together mediate core cellular functions. They are not restricted to physical interactions or co-localisation. Progulon abundance changes are primarily controlled at the level of protein synthesis and degradation. Implemented as a web app at www.proteomehd.net/progulonFinder, our approach enables the targeted search for progulons of specific cellular processes. We use it to identify a DNA replication progulon and reveal multiple new replication factors, validated by extensive phenotyping of siRNA-induced knockdowns. Progulons provide a new entry point into the molecular understanding of biological processes.

Rational search of genetic design space for a heterologous terpene metabolic pathway in Streptomyces

Modern tools in DNA synthesis and assembly give genetic engineers control over the nucleotide-level design of complex, multi-gene systems. Systematic approaches to explore genetic design space and optimize the performance of genetic constructs are lacking. Here we explore the application of a five-level Plackett-Burman fractional factorial design to improve the titer of a heterologous terpene biosynthetic pathway in Streptomyces. A library of 125 engineered gene clusters encoding the production of diterpenoid ent-atiserenoic acid (eAA) via the methylerythritol phosphate pathway was constructed and introduced into Streptomyces albidoflavus J1047 for heterologous expression. The eAA production titer varied within the library by over two orders of magnitude and host strains showed unexpected and reproducible colony morphology phenotypes. Analysis of Plackett-Burman design identified expression of dxs, the gene encoding the first and the flux-controlling enzyme, having the strongest impact on eAA titer, but with a counter-intuitive negative correlation between dxs expression and eAA production. Finally, simulation modeling was performed to determine how several plausible sources of experimental error/noise and non-linearity impact the utility of Plackett-Burman analyses.

Additional Layer of Regulation via Convergent Gene Orientation in Yeasts

Convergent gene pairs can produce transcripts with complementary sequences. We had shown that mRNA duplexes form in vivo in Saccharomyces cerevisiae via interactions of mRNA overlapping 3′-ends and can lead to posttranscriptional regulatory events. Here we show that mRNA duplex formation is restricted to convergent genes separated by short intergenic distance, independently of their 3′-untranslated region (UTR) length. We disclose an enrichment in genes involved in biological processes related to stress among these convergent genes. They are markedly conserved in convergent orientation in budding yeasts, meaning that this mode of posttranscriptional regulation could be shared in these organisms, conferring an additional level for modulating stress response. We thus investigated the mechanistic advantages potentially conferred by 3′-UTR mRNA interactions. Analysis of genome-wide transcriptome data revealed that Pat1 and Lsm1 factors, having 3′-UTR binding preference and participating to the remodeling of messenger ribonucleoprotein particles, bind differently these messenger-interacting mRNAs forming duplexes in comparison to mRNAs that do not interact (solo mRNAs). Functionally, messenger-interacting mRNAs show limited translational repression upon stress. We thus propose that mRNA duplex formation modulates the regulation of mRNA expression by limiting their access to translational repressors. Our results thus show that posttranscriptional regulation is an additional factor that determines the order of coding genes.

PCNA-associated factor KIAA0101 transcriptionally induced by ELK1 controls cell proliferation and apoptosis in nasopharyngeal carcinoma: an integrated bioinformatics and experimental study

KIAA0101, previously identified as PCNA-associated factor, is overexpressed among almost majority of human cancers and has emerged as an important regulator of cancer progression; however, its function in human nasopharyngeal carcinoma (NPC) remain unknown. Integrated bioinformatics approaches were employed to determine the KIAA0101 expressions in the NPC samples. Lentiviral vectors carrying KIAA0101 shRNA were constructed and stable transfected cells were validated by qRT-PCR and western blot. Cellular functions were then evaluated by MTT, colony formation, Brdu staining, and flow cytometry. Mechanistic studies were systematically investigated by UCSC Genome Browser, GEO, UALCAN, QIAGEN, PROMO and JASPAR, ChIP, and the cBioPortal, et al. The results showed that KIAA0101 ranked top overexpressed gene lists in GSE6631 dataset. KIAA0101 was highly expressed in NPC tissues and cell lines. Furthermore, knockdown of KIAA0101 significantly inhibited cell proliferation and DNA replication, promoted apoptosis and cell cycle arrest in vitro. Meanwhile, the mechanistic study revealed that MAP kinase phosphorylation-dependent activation of ELK1 may enhance neighbor gene expressions of KIAA0101 and TRIP4 by binding both promotor regions in the NPC cells. Taken together, our findings indicate that overexpression of KIAA0101 activated by MAP kinase phosphorylation-dependent activation of ELK1 may play an important role in NPC progression.

Early fate of exogenous promoters in E. coli

Gene gain by horizontal gene transfer is a major pathway of genome innovation in bacteria. The current view posits that acquired genes initially need to be silenced and that a bacterial chromatin protein, H-NS, plays a role in this silencing. However, we lack direct observation of the early fate of a horizontally transferred gene to prove this theory. We combine sequencing, flow cytometry and sorting, followed by microscopy to monitor gene expression and its variability after large-scale random insertions of a reporter gene in a population of Escherichia coli bacteria. We find that inserted promoters have a wide range of gene-expression variability related to their location. We find that high-expression clones carry insertions that are not correlated with H-NS binding. Conversely, binding of H-NS correlates with silencing. Finally, while most promoters show a common level of extrinsic noise, some insertions show higher noise levels. Analysis of these high-noise clones supports a scenario of switching due to transcriptional interference from divergent ribosomal promoters. Altogether, our findings point to evolutionary pathways where newly-acquired genes are not necessarily silenced, but may immediately explore a wide range of expression levels to probe the optimal ones.

Long noncoding RNA repertoire and targeting by nuclear exosome, cytoplasmic exonuclease, and RNAi in fission yeast

Long noncoding RNAs (lncRNAs), which are longer than 200 nucleotides but often unstable, contribute a substantial and diverse portion to pervasive noncoding transcriptomes. Most lncRNAs are poorly annotated and understood, although several play important roles in gene regulation and diseases. Here we systematically uncover and analyze lncRNAs in Schizosaccharomyces pombe. Based on RNA-seq data from twelve RNA-processing mutants and nine physiological conditions, we identify 5775 novel lncRNAs, nearly 4× the previously annotated lncRNAs. The expression of most lncRNAs becomes strongly induced under the genetic and physiological perturbations, most notably during late meiosis. Most lncRNAs are cryptic and suppressed by three RNA-processing pathways: the nuclear exosome, cytoplasmic exonuclease, and RNAi. Double-mutant analyses reveal substantial coordination and redundancy among these pathways. We classify lncRNAs by their dominant pathway into cryptic unstable transcripts (CUTs), Xrn1-sensitive unstable transcripts (XUTs), and Dicer-sensitive unstable transcripts (DUTs). XUTs and DUTs are enriched for antisense lncRNAs, while CUTs are often bidirectional and actively translated. The cytoplasmic exonuclease, along with RNAi, dampens the expression of thousands of lncRNAs and mRNAs that become induced during meiosis. Antisense lncRNA expression mostly negatively correlates with sense mRNA expression in the physiological, but not the genetic conditions. Intergenic and bidirectional lncRNAs emerge from nucleosome-depleted regions, upstream of positioned nucleosomes. Our results highlight both similarities and differences to lncRNA regulation in budding yeast. This broad survey of the lncRNA repertoire and characteristics in S. pombe, and the interwoven regulatory pathways that target lncRNAs, provides a rich framework for their further functional analyses.

Transcriptional noise in intact and TGF-beta treated human kidney cells; the importance of time-series designs

The transforming growth factor (TGF)-β signaling pathway plays a key role in various cellular processes. However, insufficient knowledge about the complex and sometimes paradoxical functions of this pathway hinders its therapeutic targeting. In this study, the transcriptional profile of seven mediators and downstream elements of the TGF-β pathway were assessed in TGF-β treated and untreated human kidney derived cells for 2 weeks in a time course manner. As expected the up-regulation of ACTA2 and COL1A2 was evident in the treated cells. However, we observed remarkable fluctuations in gene expression, even in the supposedly steady states. The magnitude of noise was diverse in the examined genes. Our findings underscore the significance of time-course designs for gene expression analyses and clearly show that misleading data can be obtained in single point measurements. Furthermore, we propose specific considerations in the interpretation of time-course data in the context of noisy gene expression.

Non-equilibrium phase transition in a model for supercoiling-dependent DNA transcription

We study a variant of a recently proposed non-equilibrium stochastic model for supercoiling-dependent transcription in DNA. In the case of a circular DNA molecule with overall positive supercoiling, we find a non-equilibrium phase transition between an absorbing phase, where all genes are switched off due to the supercoiling, and an active phase with a non-zero transcription rate. Mean field theory predicts that the transition should be continuous at a critical value of the background supercoiling, and we focus our analysis on this case. Our simulations suggest that the switch between the transcribed and silent phases may actually be a fluctuation-induced discontinuous transition, where the jump in the transcription rate decreases with the number of genes in the system.

The Evolution of Gene-Specific Transcriptional Noise Is Driven by Selection at the Pathway Level

Biochemical reactions within individual cells result from the interactions of molecules, typically in small numbers. Consequently, the inherent stochasticity of binding and diffusion processes generates noise along the cascade that leads to the synthesis of a protein from its encoding gene. As a result, isogenic cell populations display phenotypic variability even in homogeneous environments. The extent and consequences of this stochastic gene expression have only recently been assessed on a genome-wide scale, owing, in particular, to the advent of single-cell transcriptomics. However, the evolutionary forces shaping this stochasticity have yet to be unraveled. Here, we take advantage of two recently published data sets for the single-cell transcriptome of the domestic mouse Mus musculus to characterize the effect of natural selection on gene-specific transcriptional stochasticity. We show that noise levels in the mRNA distributions (also known as transcriptional noise) significantly correlate with three-dimensional nuclear domain organization, evolutionary constraints on the encoded protein, and gene age. However, the position of the encoded protein in a biological pathway is the main factor that explains observed levels of transcriptional noise, in agreement with models of noise propagation within gene networks. Because transcriptional noise is under widespread selection, we argue that it constitutes an important component of the phenotype and that variance of expression is a potential target of adaptation. Stochastic gene expression should therefore be considered together with the mean expression level in functional and evolutionary studies of gene expression.

It's easier to get along with the quiet neighbours

Why are neighbouring genes co‐expressed at the RNA level? While it is tempting to think that this is to enable coordination of functionally related proteins, analysis of human proteomics data by Rappsilber and colleagues (Kustatscher et al , 2017 ) suggests this is the exception not the rule. Rather it might either be just something that happens or, in some instances, an epiphenomenon of coordination of expression to enable reduced gene expression noise.

Pervasive coexpression of spatially proximal genes is buffered at the protein level

Genes are not randomly distributed in the genome. In humans, 10% of protein-coding genes are transcribed from bidirectional promoters and many more are organised in larger clusters. Intriguingly, neighbouring genes are frequently coexpressed but rarely functionally related. Here we show that coexpression of bidirectional gene pairs, and closeby genes in general, is buffered at the protein level. Taking into account the 3D architecture of the genome, we find that co-regulation of spatially close, functionally unrelated genes is pervasive at the transcriptome level, but does not extend to the proteome. We present evidence that non-functional mRNA coexpression in human cells arises from stochastic chromatin fluctuations and direct regulatory interference between spatially close genes. Protein-level buffering likely reflects a lack of coordination of post-transcriptional regulation of functionally unrelated genes. Grouping human genes together along the genome sequence, or through long-range chromosome folding, is associated with reduced expression noise. Our results support the hypothesis that the selection for noise reduction is a major driver of the evolution of genome organisation.

Flipping between Polycomb repressed and active transcriptional states introduces noise in gene expression

Polycomb repressive complexes (PRCs) are important histone modifiers, which silence gene expression; yet, there exists a subset of PRC-bound genes actively transcribed by RNA polymerase II (RNAPII). It is likely that the role of Polycomb repressive complex is to dampen expression of these PRC-active genes. However, it is unclear how this flipping between chromatin states alters the kinetics of transcription. Here, we integrate histone modifications and RNAPII states derived from bulk ChIP-seq data with single-cell RNA-sequencing data. We find that Polycomb repressive complex-active genes have greater cell-to-cell variation in expression than active genes, and these results are validated by knockout experiments. We also show that PRC-active genes are clustered on chromosomes in both two and three dimensions, and interactions with active enhancers promote a stabilization of gene expression noise. These findings provide new insights into how chromatin regulation modulates stochastic gene expression and transcriptional bursting, with implications for regulation of pluripotency and development.Polycomb repressive complexes modify histones but it is unclear how changes in chromatin states alter kinetics of transcription. Here, the authors use single-cell RNAseq and ChIPseq to find that actively transcribed genes with Polycomb marks have greater cell-to-cell variation in expression.

A Slowed Cell Cycle Stabilizes the Budding Yeast Genome

During cell division, aberrant DNA structures are detected by regulators called checkpoints that slow division to allow error correction. In addition to checkpoint-induced delay, it is widely assumed, though rarely shown, that merely slowing the cell cycle might allow more time for error detection and correction, thus resulting in a more stable genome. Fidelity by a slowed cell cycle might be independent of checkpoints. Here we tested the hypothesis that a slowed cell cycle stabilizes the genome, independent of checkpoints, in the budding yeast Saccharomyces cerevisiae We were led to this hypothesis when we identified a gene (ERV14, an ER cargo membrane protein) that when mutated, unexpectedly stabilized the genome, as measured by three different chromosome assays. After extensive studies of pathways rendered dysfunctional in erv14 mutant cells, we are led to the inference that no particular pathway is involved in stabilization, but rather the slowed cell cycle induced by erv14 stabilized the genome. We then demonstrated that, in genetic mutations and chemical treatments unrelated to ERV14, a slowed cell cycle indeed correlates with a more stable genome, even in checkpoint-proficient cells. Data suggest a delay in G2/M may commonly stabilize the genome. We conclude that chromosome errors are more rarely made or are more readily corrected when the cell cycle is slowed (even ∼15 min longer in an ∼100-min cell cycle). And, some chromosome errors may not signal checkpoint-mediated responses, or do not sufficiently signal to allow correction, and their correction benefits from this "time checkpoint."

Sequence-Modified Antibiotic Resistance Genes Provide Sustained Plasmid-Mediated Transgene Expression in Mammals

Conventional plasmid vectors are incapable of achieving sustained levels of transgene expression in vivo even in quiescent mammalian tissues because the transgene expression cassette is silenced. Transcriptional silencing results from the presence of the bacterial plasmid backbone or virtually any DNA sequence of >1 kb in length placed outside of the expression cassette. Here, we show that transcriptional silencing can be substantially forestalled by increasing the An/Tn sequence composition in the plasmid bacterial backbone. Increasing numbers of An/Tn sequences increased sustained transcription of both backbone sequences and adjacent expression cassettes. In order to recapitulate these expression profiles in compact and portable plasmid DNA backbones, we engineered the standard kanamycin or ampicillin antibiotic resistance genes, optimizing the number of An/Tn sequence without altering the encoded amino acids. The resulting vector backbones yield sustained transgene expression from mouse liver, providing generic DNA vectors capable of sustained transgene expression without additional genes or mammalian regulatory elements.

Setd5 is essential for mammalian development and the co-transcriptional regulation of histone acetylation

SET domain-containing proteins play a vital role in regulating gene expression during development through modifications in chromatin structure. Here we show that SET domain-containing 5 (Setd5) is divergently transcribed with Gt(ROSA26)Sor, is necessary for mammalian development, and interacts with the PAF1 co-transcriptional complex and other proteins. Setd5-deficient mouse embryos exhibit severe defects in neural tube formation, somitogenesis and cardiac development, have aberrant vasculogenesis in embryos, yolk sacs and placentas, and die between embryonic day 10.5 and 11.5. Setd5-deficient embryonic stem cells have impaired cellular proliferation, increased apoptosis, defective cell cycle progression, a diminished ability to differentiate into cardiomyocytes and greatly perturbed gene expression. SETD5 co-immunoprecipitates with multiple components of the PAF1 and histone deacetylase-containing NCoR complexes and is not solely required for major histone lysine methylation marks. In the absence of Setd5, histone acetylation is increased at transcription start sites and near downstream regions. These findings suggest that SETD5 functions in a manner similar to yeast Set3p and Drosophila UpSET, and that it is essential for regulating histone acetylation during gene transcription.

The regulatory role of DNA supercoiling in nucleoprotein complex assembly and genetic activity

We argue that dynamic changes in DNA supercoiling in vivo determine both how DNA is packaged and how it is accessed for transcription and for other manipulations such as recombination. In both bacteria and eukaryotes, the principal generators of DNA superhelicity are DNA translocases, supplemented in bacteria by DNA gyrase. By generating gradients of superhelicity upstream and downstream of their site of activity, translocases enable the differential binding of proteins which preferentially interact with respectively more untwisted or more writhed DNA. Such preferences enable, in principle, the sequential binding of different classes of protein and so constitute an essential driver of chromatin organization.

Stochastic Model of Supercoiling-Dependent Transcription

We propose a stochastic model for gene transcription coupled to DNA supercoiling, where we incorporate the experimental observation that polymerases create supercoiling as they unwind the DNA helix and that these enzymes bind more favorably to regions where the genome is unwound. Within this model, we show that when the transcriptionally induced flux of supercoiling increases, there is a sharp crossover from a regime where torsional stresses relax quickly and gene transcription is random, to one where gene expression is highly correlated and tightly regulated by supercoiling. In the latter regime, the model displays transcriptional bursts, waves of supercoiling, and up regulation of divergent or bidirectional genes. It also predicts that topological enzymes which relax twist and writhe should provide a pathway to down regulate transcription.

GFOD1 and peejar are promising markers for clear-cell renal cell carcinoma disease progression

Clear cell renal cell carcinoma (ccRCC) is a common genitourinary malignancy. The molecular pathogenesis of ccRCC remains unclear and biomarkers for daily practice were still limited. We performed an integrative analysis of two public ccRCC microarray datasets, E-GEOD-22541 and E-MTAB-1050, The candidate differential expression genes (DEG) were then confirmed in the E-GEOD-53757 dataset. In addition, an independent cohort of 50 ccRCC and 36 non-tumor kidney tissues were analyzed to examine the selected DGEs by qRT-PCR. We identified and validated two DEGs, namely GFOD1 and peejar, which were significantly up-regulated in ccRCC compared with normal renal tissues (p < 0.001). Moreover, the expression of these two genes are related to histological grade and stage and decrease of their expression correlated with disease progression (p < 0.05). Furthermore, we found the expression of peejar was positively correlated with the expression of GFOD1 in ccRCC tissue, with Pearson correlation coefficiency reaching 0.939 (p < 0.001). GFOD1 and peejar were novel genes correlated with ccRCC disease progression and patients' poor prognosis.

Integration of multi-omics data of a genome-reduced bacterium: Prevalence of post-transcriptional regulation and its correlation with protein abundances

We developed a comprehensive resource for the genome-reduced bacterium Mycoplasma pneumoniae comprising 1748 consistently generated ‘-omics’ data sets, and used it to quantify the power of antisense non-coding RNAs (ncRNAs), lysine acetylation, and protein phosphorylation in predicting protein abundance (11%, 24% and 8%, respectively). These factors taken together are four times more predictive of the proteome abundance than of mRNA abundance. In bacteria, post-translational modifications (PTMs) and ncRNA transcription were both found to increase with decreasing genomic GC-content and genome size. Thus, the evolutionary forces constraining genome size and GC-content modify the relative contributions of the different regulatory layers to proteome homeostasis, and impact more genomic and genetic features than previously appreciated. Indeed, these scaling principles will enable us to develop more informed approaches when engineering minimal synthetic genomes.

Stochastic developmental variation, an epigenetic source of phenotypic diversity with far-reaching biological consequences

This article reviews the production of different phenotypes from the same genotype in the same environment by stochastic cellular events, nonlinear mechanisms during patterning and morphogenesis, and probabilistic self-reinforcing circuitries in the adult life. These aspects of phenotypic variation are summarized under the term 'stochastic developmental variation' (SDV) in the following. In the past, SDV has been viewed primarily as a nuisance, impairing laboratory experiments, pharmaceutical testing, and true-to-type breeding. This article also emphasizes the positive biological effects of SDV and discusses implications for genotype-to-phenotype mapping, biological individuation, ecology, evolution, and applied biology. There is strong evidence from experiments with genetically identical organisms performed in narrowly standardized laboratory set-ups that SDV is a source of phenotypic variation in its own right aside from genetic variation and environmental variation. It is obviously mediated by molecular and higher-order epigenetic mechanisms. Comparison of SDV in animals, plants, fungi, protists, bacteria, archaeans, and viruses suggests that it is a ubiquitous and phylogenetically old phenomenon. In animals, it is usually smallest for morphometric traits and highest for life history traits and behaviour. SDV is thought to contribute to phenotypic diversity in all populations but is particularly relevant for asexually reproducing and genetically impoverished populations, where it generates individuality despite genetic uniformity. In each generation, SDV produces a range of phenotypes around a well-adapted target phenotype, which is interpreted as a bet-hedging strategy to cope with the unpredictability of dynamic environments. At least some manifestations of SDV are heritable, adaptable, selectable, and evolvable, and therefore, SDV may be seen as a hitherto overlooked evolution factor. SDV is also relevant for husbandry, agriculture, and medicine because most pathogens are asexuals that exploit this third source of phenotypic variation to modify infectivity and resistance to antibiotics. Since SDV affects all types of organisms and almost all aspects of life, it urgently requires more intense research and a better integration into biological thinking.

Visualizing the activity of Escherichia coli divergent promoters and probing their dependence on superhelical density using dual-colour fluorescent reporter vector

Mosaic pattern of transcription in alternating directions is a common feature of prokaryotic and eukaryotic genomes which rationality and origin remain enigmatic. In Escherichia coli approximately 25% of genes comprise pairs of topologically linked divergently transcribed units. Given that transcriptional complex formation at each promoter in the pair induces topological changes and is itself sensitive to DNA structural perturbations, study of the functional anatomy in such areas requires special approaches. Here we suggested the dual-colour promoter probe vector which may become an ideal tool for divergent transcription profiling. The vector was used to characterize the specific genomic region nearby appY with multiple bidirectional promoters predicted in silico. Only three promoters of this region were shown to be engaged in the transcription initiation resulting in the expression of reporter genes. RNA product transcribed in antisense direction is suggested as a novel RNA. Nalidixin-induced topological modulation differentially affected transcription in sense and antisense directions thus exemplifying anticooperative mode in the response to topological alterations.

Oleic acid metabolism via a conserved cytochrome P450 system-mediated ω-hydroxylation in the bark beetle-associated fungus Grosmannia clavigera

The bark beetle-associated fungus Grosmannia clavigera participates in the large-scale destruction of pine forests. In the tree, it must tolerate saturating levels of toxic conifer defense chemicals (e.g. monoterpenes). The fungus can metabolize some of these compounds through the ß-oxidation pathway and use them as a source of carbon. It also uses carbon from pine triglycerides, where oleic acid is the most common fatty acid. High levels of free fatty acids, however, are toxic and can cause additional stress during host colonization. Fatty acids induce expression of neighboring genes encoding a cytochrome P450 (CYP630B18) and its redox partner, cytochrome P450 reductase (CPR2). The aim of this work was to study the function of this novel P450 system. Using LC/MS, we biochemically characterized CYP630 as a highly specific oleic acid ω-hydroxylase. We explain oleic acid specificity using protein interaction modeling. Our results underscore the importance of ω-oxidation when the main ß-oxidation pathway may be overwhelmed by other substrates such as host terpenoid compounds. Because this CYP-CPR gene cluster is evolutionarily conserved, our work has implications for metabolism studies in other fungi.

The parasite Trichomonas vaginalis expresses thousands of pseudogenes and long non-coding RNAs independently from functional neighbouring genes

Background

The human pathogen Trichomonas vaginalis is a parabasalian flagellate that is estimated to infect 3% of the world’s population annually. With a 160 megabase genome and up to 60,000 genes residing in six chromosomes, the parasite has the largest genome among sequenced protists. Although it is thought that the genome size and unusual large coding capacity is owed to genome duplication events, the exact reason and its consequences are less well studied.

Results

Among transcriptome data we found thousands of instances, in which reads mapped onto genomic loci not annotated as genes, some reaching up to several kilobases in length. At first sight these appear to represent long non-coding RNAs (lncRNAs), however, about half of these lncRNAs have significant sequence similarities to genomic loci annotated as protein-coding genes. This provides evidence for the transcription of hundreds of pseudogenes in the parasite. Conventional lncRNAs and pseudogenes are expressed in Trichomonas through their own transcription start sites and independently from flanking genes in Trichomonas. Expression of several representative lncRNAs was verified through reverse-transcriptase PCR in different T. vaginalis strains and case studies exclude the use of alternative start codons or stop codon suppression for the genes analysed.

Conclusion

Our results demonstrate that T. vaginalis expresses thousands of intergenic loci, including numerous transcribed pseudogenes. In contrast to yeast these are expressed independently from neighbouring genes. Our results furthermore illustrate the effect genome duplication events can have on the transcriptome of a protist. The parasite’s genome is in a steady state of changing and we hypothesize that the numerous lncRNAs could offer a large pool for potential innovation from which novel proteins or regulatory RNA units could evolve.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-906) contains supplementary material, which is available to authorized users.

Uncovering the functional constraints underlying the genomic organization of the odorant-binding protein genes

Animal olfactory systems have a critical role for the survival and reproduction of individuals. In insects, the odorant-binding proteins (OBPs) are encoded by a moderately sized gene family, and mediate the first steps of the olfactory processing. Most OBPs are organized in clusters of a few paralogs, which are conserved over time. Currently, the biological mechanism explaining the close physical proximity among OBPs is not yet established. Here, we conducted a comprehensive study aiming to gain insights into the mechanisms underlying the OBP genomic organization. We found that the OBP clusters are embedded within large conserved arrangements. These organizations also include other non-OBP genes, which often encode proteins integral to plasma membrane. Moreover, the conservation degree of such large clusters is related to the following: 1) the promoter architecture of the confined genes, 2) a characteristic transcriptional environment, and 3) the chromatin conformation of the chromosomal region. Our results suggest that chromatin domains may restrict the location of OBP genes to regions having the appropriate transcriptional environment, leading to the OBP cluster structure. However, the appropriate transcriptional environment for OBP and the other neighbor genes is not dominated by reduced levels of expression noise. Indeed, the stochastic fluctuations in the OBP transcript abundance may have a critical role in the combinatorial nature of the olfactory coding process.

Functional gene groups are concentrated within chromosomes, among chromosomes and in the nuclear space of the human genome

Genomes undergo changes in organization as a result of gene duplications, chromosomal rearrangements and local mutations, among other mechanisms. In contrast to prokaryotes, in which genes of a common function are often organized in operons and reside contiguously along the genome, most eukaryotes show much weaker clustering of genes by function, except for few concrete functional groups. We set out to check systematically if there is a relation between gene function and gene organization in the human genome. We test this question for three types of functional groups: pairs of interacting proteins, complexes and pathways. We find a significant concentration of functional groups both in terms of their distance within the same chromosome and in terms of their dispersal over several chromosomes. Moreover, using Hi-C contact map of the tendency of chromosomal segments to appear close in the 3D space of the nucleus, we show that members of the same functional group that reside on distinct chromosomes tend to co-localize in space. The result holds for all three types of functional groups that we tested. Hence, the human genome shows substantial concentration of functional groups within chromosomes and across chromosomes in space.

Torsion-mediated interaction between adjacent genes

DNA torsional stress is generated by virtually all biomolecular processes involving the double helix, in particular transcription where a significant level of stress propagates over several kilobases. If another promoter is located in this range, this stress may strongly modify its opening properties, and hence facilitate or hinder its transcription. This mechanism implies that transcribed genes distant of a few kilobases are not independent, but coupled by torsional stress, an effect for which we propose the first quantitative and systematic model. In contrast to previously proposed mechanisms of transcriptional interference, the suggested coupling is not mediated by the transcription machineries, but results from the universal mechanical features of the double-helix. The model shows that the effect likely affects prokaryotes as well as eukaryotes, but with different consequences owing to their different basal levels of torsion. It also depends crucially on the relative orientation of the genes, enhancing the expression of eukaryotic divergent pairs while reducing that of prokaryotic convergent ones. To test the in vivo influence of the torsional coupling, we analyze the expression of isolated gene pairs in the Drosophila melanogaster genome. Their orientation and distance dependence is fully consistent with the model, suggesting that torsional gene coupling may constitute a widespread mechanism of (co)regulation in eukaryotes.

Suppression of expression between adjacent genes within heterologous modules in yeast

Recent studies have shown that proximal arrangement of multiple genes can have complex effects on gene expression. For example, in the case of heterologous gene expression modules, certain arrangements of the selection marker and the gene expression cassette may have unintended consequences that limit the predictability and interpretability of module behaviors. The relationship between arrangement and expression has not been systematically characterized within heterologous modules to date. In this study, we quantitatively measured gene expression patterns of the selection marker (KlURA3 driven by the promoter, pKlURA) and the gene expression cassette (GFP driven by the galactose-inducible GAL1 promoter, pGAL1) in all their possible relative arrangements in Saccharomyces cerevisiae. First, we observed that pKlURA activity depends strongly on the relative arrangement and the activity of pGAL1. Most notably, we observed transcriptional suppression in the case of divergent arrangements: pKlURA activity was reduced when pGAL1 was inactive. Based on our nucleosome occupancy data, we attribute the observed transcriptional reduction to nucleosome repositioning. Second, we observed that pGAL1 activity also depends on the relative arrangement of pKlURA. In particular, strains with divergent promoters showed significantly different pGAL1 activation patterns from other strains, but only when their growth was compromised by lack of uracil. We reasoned that this difference in pGAL1 activation patterns arises from arrangement-dependent pKlURA activity that can affect the overall cell physiology (i.e., cell growth and survival in the uracil-depleted condition). Our results underscore the necessity to consider ramifications of promoter arrangement when using synthetic gene expression modules.

Divergent transcription: a driving force for new gene origination?

The mammalian genome is extensively transcribed, a large fraction of which is divergent transcription from promoters and enhancers that is tightly coupled with active gene transcription. Here, we propose that divergent transcription may shape the evolution of the genome by new gene origination.

Functional networking of human divergently paired genes (DPGs)

Divergently paired genes (DPGs), also known as bidirectional (head-to-head positioned) genes, are conserved across species and lineages, and thus deemed to be exceptional in genomic organization and functional regulation. Despite previous investigations on the features of their conservation and gene organization, the functional relationship among DPGs in a given species and lineage has not been thoroughly clarified. Here we report a network-based comprehensive analysis on human DPGs and our results indicate that the two members of the DPGs tend to participate in different biological processes while enforcing related functions as modules. Comparing to randomly paired genes as a control, the DPG pairs have a tendency to be clustered in similar “cellular components” and involved in similar “molecular functions”. The functional network bridged by DPGs consists of three major modules. The largest module includes many house-keeping genes involved in core cellular activities. This module also shows low variation in expression in both CNS (central nervous system) and non-CNS tissues. Based on analyses of disease transcriptome data, we further suggest that this particular module may play crucial roles in HIV infection and its disease mechanism.

Epigenetic regulatory elements associate with specific histone modifications to prevent silencing of telomeric genes

In eukaryotic cells, transgene expression levels may be limited by an unfavourable chromatin structure at the integration site. Epigenetic regulators are DNA sequences which may protect transgenes from such position effect. We evaluated different epigenetic regulators for their ability to protect transgene expression at telomeres, which are commonly associated to low or inconsistent expression because of their repressive chromatin environment. Although to variable extents, matrix attachment regions (MARs), ubiquitous chromatin opening element (UCOE) and the chicken cHS4 insulator acted as barrier elements, protecting a telomeric-distal transgene from silencing. MARs also increased the probability of silent gene reactivation in time-course experiments. Additionally, all MARs improved the level of expression in non-silenced cells, unlike other elements. MARs were associated to histone marks usually linked to actively expressed genes, especially acetylation of histone H3 and H4, suggesting that they may prevent the spread of silencing chromatin by imposing acetylation marks on nearby nucleosomes. Alternatively, an UCOE was found to act by preventing deposition of repressive chromatin marks. We conclude that epigenetic DNA elements used to enhance and stabilize transgene expression all have specific epigenetic signature that might be at the basis of their mode of action.

Physical linkage of metabolic genes in fungi is an adaptation against the accumulation of toxic intermediate compounds

Genomic analyses have proliferated without being tied to tangible phenotypes. For example, although coordination of both gene expression and genetic linkage have been offered as genetic mechanisms for the frequently observed clustering of genes participating in fungal metabolic pathways, elucidation of the phenotype(s) favored by selection, resulting in cluster formation and maintenance, has not been forthcoming. We noted that the cause of certain well-studied human metabolic disorders is the accumulation of toxic intermediate compounds (ICs), which occurs when the product of an enzyme is not used as a substrate by a downstream neighbor in the metabolic network. This raises the hypothesis that the phenotype favored by selection to drive gene clustering is the mitigation of IC toxicity. To test this, we examined 100 diverse fungal genomes for the simplest type of cluster, gene pairs that are both metabolic neighbors and chromosomal neighbors immediately adjacent to each other, which we refer to as "double neighbor gene pairs" (DNGPs). Examination of the toxicity of their corresponding ICs shows that, compared with chromosomally nonadjacent metabolic neighbors, DNGPs are enriched for ICs that have acutely toxic LD50 doses or reactive functional groups. Furthermore, DNGPs are significantly more likely to be divergently oriented on the chromosome; remarkably, ∼40% of these DNGPs have ICs known to be toxic. We submit that the structure of synteny in metabolic pathways of fungi is a signature of selection for protection against the accumulation of toxic metabolic intermediates.

Gene regulation by the act of long non-coding RNA transcription

Long non-protein-coding RNAs (lncRNAs) are proposed to be the largest transcript class in the mouse and human transcriptomes. Two important questions are whether all lncRNAs are functional and how they could exert a function. Several lncRNAs have been shown to function through their product, but this is not the only possible mode of action. In this review we focus on a role for the process of lncRNA transcription, independent of the lncRNA product, in regulating protein-coding-gene activity in cis. We discuss examples where lncRNA transcription leads to gene silencing or activation, and describe strategies to determine if the lncRNA product or its transcription causes the regulatory effect.

G-NEST: a gene neighborhood scoring tool to identify co-conserved, co-expressed genes

BACKGROUND

In previous studies, gene neighborhoods-spatial clusters of co-expressed genes in the genome-have been defined using arbitrary rules such as requiring adjacency, a minimum number of genes, a fixed window size, or a minimum expression level. In the current study, we developed a Gene Neighborhood Scoring Tool (G-NEST) which combines genomic location, gene expression, and evolutionary sequence conservation data to score putative gene neighborhoods across all possible window sizes simultaneously.

RESULTS

Using G-NEST on atlases of mouse and human tissue expression data, we found that large neighborhoods of ten or more genes are extremely rare in mammalian genomes. When they do occur, neighborhoods are typically composed of families of related genes. Both the highest scoring and the largest neighborhoods in mammalian genomes are formed by tandem gene duplication. Mammalian gene neighborhoods contain highly and variably expressed genes. Co-localized noisy gene pairs exhibit lower evolutionary conservation of their adjacent genome locations, suggesting that their shared transcriptional background may be disadvantageous. Genes that are essential to mammalian survival and reproduction are less likely to occur in neighborhoods, although neighborhoods are enriched with genes that function in mitosis. We also found that gene orientation and protein-protein interactions are partially responsible for maintenance of gene neighborhoods.

CONCLUSIONS

Our experiments using G-NEST confirm that tandem gene duplication is the primary driver of non-random gene order in mammalian genomes. Non-essentiality, co-functionality, gene orientation, and protein-protein interactions are additional forces that maintain gene neighborhoods, especially those formed by tandem duplicates. We expect G-NEST to be useful for other applications such as the identification of core regulatory modules, common transcriptional backgrounds, and chromatin domains. The software is available at http://docpollard.org/software.html.

Interplay of gene expression noise and ultrasensitive dynamics affects bacterial operon organization

Bacterial chromosomes are organized into polycistronic cotranscribed operons, but the evolutionary pressures maintaining them are unclear. We hypothesized that operons alter gene expression noise characteristics, resulting in selection for or against maintaining operons depending on network architecture. Mathematical models for 6 functional classes of network modules showed that three classes exhibited decreased noise and 3 exhibited increased noise with same-operon cotranscription of interacting proteins. Noise reduction was often associated with a decreased chance of reaching an ultrasensitive threshold. Stochastic simulations of the lac operon demonstrated that the predicted effects of transcriptional coupling hold for a complex network module. We employed bioinformatic analysis to find overrepresentation of noise-minimizing operon organization compared with randomized controls. Among constitutively expressed physically interacting protein pairs, higher coupling frequencies appeared at lower expression levels, where noise effects are expected to be dominant. Our results thereby suggest an important role for gene expression noise, in many cases interacting with an ultrasensitive switch, in maintaining or selecting for operons in bacterial chromosomes.

In some species, most notably bacteria, chromosomal genes are arranged into clusters called operons. In operons, the process of transcription is physically coupled: a single pass of the RNA polymerase enzyme reading that region of the chromosome simultaneously produces messenger RNA encoding multiple proteins. So far, we do not have a satisfying explanation for what evolutionary forces have maintained operons on bacterial chromosomes. We hypothesized that different types of interactions between operon-coded proteins affect how strongly operons are selected for between two genes. The proposed mechanism for this effect is that operons correlate gene expression noise, changing how it manifests in the post-translational network depending on the type of protein interaction. Mathematical models demonstrate that operons reduce noise for some types of interactions but not others. We found that operon-dependent noise reduction has an underlying dependence on surprisingly high sensitivity of the network to the ratio of proteins from each gene. Databases of genetic information show that E. coli has operons more frequently than random if operons reduce noise for the type of interaction various gene pairs have, but not otherwise. Our study thus provides an example of how the architecture of post-translational networks affects bacterial evolution.

Balancing noise and plasticity in eukaryotic gene expression

BACKGROUND

Coupling the control of expression stochasticity (noise) to the ability of expression change (plasticity) can alter gene function and influence adaptation. A number of factors, such as transcription re-initiation, strong chromatin regulation or genome neighboring organization, underlie this coupling. However, these factors do not necessarily combine in equivalent ways and strengths in all genes. Can we identify then alternative architectures that modulate in distinct ways the linkage of noise and plasticity?

RESULTS

Here we first show that strong chromatin regulation, commonly viewed as a source of coupling, can lead to plasticity without noise. The nature of this regulation is relevant too, with plastic but noiseless genes being subjected to general activators whereas plastic and noisy genes experience more specific repression. Contrarily, in genes exhibiting poor transcriptional control, it is translational efficiency what separates noise from plasticity, a pattern related to transcript length. This additionally implies that genome neighboring organization -as modifier- appears only effective in highly plastic genes. In this class, we confirm bidirectional promoters (bipromoters) as a configuration capable to reduce coupling by abating noise but also reveal an important trade-off, since bipromoters also decrease plasticity. This presents ultimately a paradox between intergenic distances and modulation, with short intergenic distances both associated and disassociated to noise at different plasticity levels.

CONCLUSIONS

Balancing the coupling among different types of expression variability appears as a potential shaping force of genome regulation and organization. This is reflected in the use of different control strategies at genes with different sets of functional constraints.

Long noncoding RNAs in C. elegans

Thousands of long noncoding RNAs (lncRNAs) have been found in vertebrate animals, a few of which have known biological roles. To better understand the genomics and features of lncRNAs in invertebrates, we used available RNA-seq, poly(A)-site, and ribosome-mapping data to identify lncRNAs of Caenorhabditis elegans. We found 170 long intervening ncRNAs (lincRNAs), which had single- or multiexonic structures that did not overlap protein-coding transcripts, and about sixty antisense lncRNAs (ancRNAs), which were complementary to protein-coding transcripts. Compared to protein-coding genes, the lncRNA genes tended to be expressed in a stage-dependent manner. Approximately 25% of the newly identified lincRNAs showed little signal for sequence conservation and mapped antisense to clusters of endogenous siRNAs, as would be expected if they serve as templates and targets for these siRNAs. The other 75% tended to be more conserved and included lincRNAs with intriguing expression and sequence features associating them with processes such as dauer formation, male identity, sperm formation, and interaction with sperm-specific mRNAs. Our study provides a glimpse into the lncRNA content of a nonvertebrate animal and a resource for future studies of lncRNA function.

Dynamics of transcription of closely spaced promoters in Escherichia coli, one event at a time

Many pairs of genes in Escherichia coli are driven by closely spaced promoters. We study the dynamics of expression of such pairs of genes driven by a model at the molecule and nucleotide level with delayed stochastic dynamics as a function of the binding affinity of the RNA polymerase to the promoter region, of the geometry of the promoter, of the distance between transcription start sites (TSSs) and of the repression mechanism. We find that the rate limiting steps of transcription at the TSS, the closed and open complex formations, strongly affect the kinetics of RNA production for all promoter configurations. Beyond a certain rate of transcription initiation events, we find that the interference between polymerases correlates the dynamics of production of the two RNA molecules from the two TSS and affects the distribution of intervals between consecutive productions of RNA molecules. The degree of correlation depends on the geometry, the distance between TSSs and repressors. Small changes in the distance between TSSs can cause abrupt changes in behavior patterns, suggesting that the sequence between adjacent promoters may be subject to strong selective pressure. The results provide better understanding on the sequence level mechanisms of transcription regulation in bacteria and may aid in the genetic engineering of artificial circuits based on closely spaced promoters.

Error prevention and mitigation as forces in the evolution of genes and genomes

Why are short introns rarely a multiple of three nucleotides long? Why do essential genes cluster? Why are genes in operons often lined up in the order in which they are needed in the encoded pathway? In this Opinion article, we argue that these and many other - ostensibly disparate - observations are all pieces of an emerging picture in which multiple aspects of gene anatomy and genome architecture have evolved in response to error-prone gene expression.

Coexpression of linked gene pairs persists long after their separation

In many eukaryotes, physically linked gene pairs tend to be coexpressed. However, it is still controversial to what extent this neighbor coexpression is maintained by selection and to what extent it is nonselective, purely mechanistic "leaky expression." Here, we analyze expression patterns of gene pairs that have lost their linkage in the evolution of Saccharomyces cerevisiae since its last common ancestor with Kluyveromyces waltii or that were never linked in the S. cerevisiae lineage but became neighbors in a related yeast. We demonstrate that coexpression of many linked genes is retained long after their separation and is thus likely to be functionally important. In addition, unlinked gene pairs that recently became neighbors in other yeast species tend to be coexpressed in S. cerevisiae. This suggests that natural selection often favors chromosomal rearrangements in which coexpressed genes become neighbors. Contrary to previous suggestions, selectively favorable coexpression appears not to be restricted to bidirectional promoters.