MyoD1 promotes the transcription of BIK and plays an apoptosis-promoting role in the development of gastric cancer.

Myogenic differentiation (MyoD) 1, which is known as a pivotal transcription factor during myogenesis, has been proven dysregulated in several cancers. However, litter is known about the precise role and downstream genes of MyoD1 in gastric cancer (GC) cells. Here, we report that MyoD1 is lowly expressed in primary GC tissues and cells. In our experiments, overexpression of MyoD1 inhibited cell proliferation. Downstream genes of MyoD1 regulation were investigated using RNA-Seq. As a result, 138 up-regulated genes and 20 down-regulated genes and 27 up-regulated lncRNAs and 20 down-regulated lncRNAs were identified in MyoD1 overexpressed MKN-45 cells, which participated in epithelial cell signaling in Helicobacter pylori infection, glycosaminoglycan biosynthesis (keratan sulfate), notch signaling pathway, and others. Among these genes, BIK was directly regulated by MyoD1 in GC cells and inhibited cancer cell proliferation. The BIK knockdown rescued the effects of MyoD1 overexpression on GC cells. In conclusion, MyoD1 inhibited cell proliferation via 158 genes and 47 lncRNAs downstream directly or indirectly that participated in multiple signaling pathways in GC, and among these, MyoD1 promotes BIK transcription by binding to its promoter, then promotes BIK-Bcl2-caspase 3 axis and regulates GC cell apoptosis.

Expression of human RECQL5 in Saccharomyces cerevisiae causes transcription defects and transcription-associated genome instability.

RECQL5 is a member of the conserved RecQ family of DNA helicases involved in the maintenance of genome stability that is specifically found in higher eukaryotes and associates with the elongating RNA polymerase II. To expand our understanding of its function we expressed human RECQL5 in the yeast Saccharomyces cerevisiae, which does not have a RECQL5 ortholog. We found that RECQL5 expression leads to cell growth inhibition, increased genotoxic sensitivity and transcription-associated hyperrecombination. Chromatin immunoprecipitation and transcriptomic analysis of yeast cells expressing human RECQL5 shows that this is recruited to transcribed genes and although it causes only a weak impact on gene expression, in particular at G + C-rich genes, it leads to a transcription termination defect detected as readthrough transcription. The data indicate that the interaction between RNAPII and RECQL5 is conserved from yeast to humans. Unexpectedly, however, the RECQL5-ID mutant, previously shown to have reduced the association with RNAPII in vitro, associates with the transcribing polymerase in cells. As a result, expression of RECQL5-ID leads to similar although weaker phenotypes than wild-type RECQL5 that could be transcription-mediated. Altogether, the data suggests that RECQL5 has the intrinsic ability to function in transcription-dependent and independent genome dynamics in S. cerevisiae.

SUMOylation of zebrafish transcription factor Zbtb21 affects its transcription activity.

Background: Post-translational modification by Small Ubiquitin-like MOdifier (SUMO) is an important mechanism to regulate protein activity, protein stability, and localization of substrates. Zbtb21 is a zinc finger and BTB (Broad-complex, Tram-track and Bric à brac) domain-containing transcription factor. Bioinformatic prediction suggests several putative SUMOylated sites in Zbtb21 protein. Methods: Two evolutionarily conserved lysine residues in Zbtb21 protein were mutated alone or in combination to disrupt the binding with SUMO molecules. Western blot and co-immunoprecipitation analyses were performed to detect the SUMOylation state of wild type and mutant Zbtb21 proteins, respectively. Luciferase reporter assays were conducted to evaluate their transcription activities. Meanwhile, immunofluorescence staining was carried out to show their sub-nuclear localizations. Finally, co-immunoprecipitation was performed to detect the interaction between Zbtb21 and its partners. Results: Phylogenetically conserved lysines 419 and 845 of zebrafish Zbtb21 protein can be conjugated with SUMO molecules. SUMOylation does not affect the subcellular localization and protein stability of Zbtb21, as well as the interaction with Zbtb14 or Zbtb21. Nevertheless, luciferase reporter assays revealed that Zbtb21 is a dual-function transcription factor which exerts activation or repression effect on different promoters, and SUMOylation can modulate the transcriptional activity of Zbtb21 in regulating downstream target genes. Hence, Zbtb21 is identified as a novel substrate of SUMOylation, which would be important for its function. Conclusions: Zebrafish Zbtb21 protein can be SUMOylated on lysines 419 and 845, which is evolutionary conserved. SUMOylation affects the dual role of Zbtb21 on transcription.

Temperature-driven coordination of circadian transcriptional regulation.

The circadian clock is an evolutionarily-conserved molecular oscillator that enables species to anticipate rhythmic changes in their environment. At a molecular level, the core clock genes induce circadian oscillations in thousands of genes in a tissue-specific manner, orchestrating myriad biological processes. While previous studies have investigated how the core clock circuit responds to environmental perturbations such as temperature, the downstream effects of such perturbations on circadian regulation remain poorly understood. By analyzing bulk-RNA sequencing of Drosophila fat bodies harvested from flies subjected to different environmental conditions, we demonstrate a highly condition-specific circadian transcriptome: genes are cycling in a temperature-specific manner, and the distributions of their phases also differ between the two conditions. Further employing a reference-based gene regulatory network (Reactome), we find evidence of increased gene-gene coordination at low temperatures and synchronization of rhythmic genes that are network neighbors. We report that the phase differences between cycling genes increase as a function of geodesic distance in the low temperature condition, suggesting increased coordination of cycling on the gene regulatory network. Our results suggest a potential mechanism whereby the circadian clock mediates the fly's response to seasonal changes in temperature.

Non-canonical functions of enhancers: regulation of RNA polymerase III transcription, DNA replication, and V(D)J recombination.

Enhancers are the key regulators of other DNA-based processes by virtue of their unique ability to generate nucleosome-depleted regions in a highly regulated manner. Enhancers regulate cell-type-specific transcription of tRNA genes by RNA polymerase III (Pol III). They are also responsible for the binding of the origin replication complex (ORC) to DNA replication origins, thereby regulating origin utilization, replication timing, and replication-dependent chromosome breaks. Additionally, enhancers regulate V(D)J recombination by increasing access of the recombination-activating gene (RAG) recombinase to target sites and by generating non-coding enhancer RNAs and localized regions of trimethylated histone H3-K4 recognized by the RAG2 PHD domain. Thus, enhancers represent the first step in decoding the genome, and hence they regulate biological processes that, unlike RNA polymerase II (Pol II) transcription, do not have dedicated regulatory proteins.

Cancer cells forgo translating mRNA transcribed from genes of nonspecialized tasks.

The coupling of transcription and translation enables prokaryotes to regulate mRNA stability and reduce nonfunctional transcripts. Eukaryotes evolved other means to perform these functions. Here, we quantify the disparity between gene expression and protein levels and attempt to explain its origins. We collected publicly available simultaneous measurements of gene expression, protein level, division rate, and growth inhibition of breast cancer cells under drug perturbation. We used the cell lines as entities with shared origin, different evolutionary trajectories, and cancer hallmarks to define tasks subject to specializing and trading-off. We observed varying average mRNA and protein correlation across cell lines, and it was consistently higher for the gene products in the cancer hallmarks. The enrichment of hallmark gene products signifies the resources invested in it as a task. Enrichment based on mRNA or protein abundance corresponds to the relative resources dedicated to transcription and translation. The differences in gene- and protein-based enrichment correlated with nominal division rates but not growth inhibition under drug perturbations. Comparing the range of enrichment scores of the hallmarks within each cell signifies the resources dedicated to each. Cells appear to have a wider range of enrichment in protein synthesis relative to gene transcription. The difference and range of enrichment of the hallmark genes and proteins correlated with cell division and inhibition in response to drug treatments. We posit that cancer cells may express the genes coding for seemingly nonspecialized tasks but do not translate them to the corresponding proteins. This trade-off may cost the cells under normal conditions but confer benefits during stress.

ZEB1 controls a lineage-specific transcriptional program essential for melanoma cell state transitions.

Cell plasticity sustains intra-tumor heterogeneity and treatment resistance in melanoma. Deciphering the transcriptional mechanisms governing reversible phenotypic transitions between proliferative/differentiated and invasive/stem-like states is required. Expression of the ZEB1 transcription factor is frequently activated in melanoma, where it fosters adaptive resistance to targeted therapies. Here, we performed a genome-wide characterization of ZEB1 transcriptional targets, by combining ChIP-sequencing and RNA-sequencing, upon phenotype switching in melanoma models. We identified and validated ZEB1 binding peaks in the promoter of key lineage-specific genes crucial for melanoma cell identity. Mechanistically, ZEB1 negatively regulates SOX10-MITF dependent proliferative/melanocytic programs and positively regulates AP-1 driven invasive and stem-like programs. Comparative analyses with breast carcinoma cells revealed lineage-specific ZEB1 binding, leading to the design of a more reliable melanoma-specific ZEB1 regulon. We then developed single-cell spatial multiplexed analyses to characterize melanoma cell states intra-tumoral heterogeneity in human melanoma samples. Combined with scRNA-Seq analyses, our findings confirmed increased ZEB1 expression in Neural-Crest-like cells and mesenchymal cells, underscoring its significance in vivo in both populations. Overall, our results define ZEB1 as a major transcriptional regulator of cell states transitions and provide a better understanding of lineage-specific transcriptional programs sustaining intra-tumor heterogeneity in melanoma.

RNA polymerase collisions and their role in transcription.

RNA polymerases are the central enzymes of gene expression and function frequently in either a head-on or co-directional manner on the busy DNA track. Whether and how these collisions between RNA polymerases contribute to transcriptional regulation is mysterious. Increasing evidence from biochemical and single-molecule studies suggests that RNA polymerase collisions function as an important regulator to fine-tune transcription, rather than creating deleterious "traffic jams". This review summarizes the recent progress on elucidating the consequences of RNA polymerase collisions during transcription and highlights the significance of cooperation and coordination between RNA polymerases.

An emerging paradigm in epigenetic marking: coordination of transcription and replication.

DNA replication and RNA transcription both utilize DNA as a template and therefore need to coordinate their activities. The predominant theory in the field is that in order for the replication fork to proceed, transcription machinery has to be evicted from DNA until replication is complete. If that does not occur, these machineries collide, and these collisions elicit various repair mechanisms which require displacement of one of the enzymes, often RNA polymerase, in order for replication to proceed. This model is also at the heart of the epigenetic bookmarking theory, which implies that displacement of RNA polymerase during replication requires gradual re-building of chromatin structure, which guides recruitment of transcriptional proteins and resumption of transcription. We discuss these theories but also bring to light newer data that suggest that these two processes may not be as detrimental to one another as previously thought. This includes findings suggesting that these processes can occur without fork collapse and that RNA polymerase may only be transiently displaced during DNA replication. We discuss potential mechanisms by which RNA polymerase may be retained at the replication fork and quickly rebind to DNA post-replication. These discoveries are important, not only as new evidence as to how these two processes are able to occur harmoniously but also because they have implications on how transcriptional programs are maintained through DNA replication. To this end, we also discuss the coordination of replication and transcription in light of revising the current epigenetic bookmarking theory of how the active gene status can be transmitted through S phase.

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.

Transcription-replication interactions reveal bacterial genome regulation.

Organisms determine the transcription rates of thousands of genes through a few modes of regulation that recur across the genome1. In bacteria, the relationship between the regulatory architecture of a gene and its expression is well understood for individual model gene circuits2,3. However, a broader perspective of these dynamics at the genome scale is lacking, in part because bacterial transcriptomics has hitherto captured only a static snapshot of expression averaged across millions of cells4. As a result, the full diversity of gene expression dynamics and their relation to regulatory architecture remains unknown. Here we present a novel genome-wide classification of regulatory modes based on the transcriptional response of each gene to its own replication, which we term the transcription-replication interaction profile (TRIP). Analysing single-bacterium RNA-sequencing data, we found that the response to the universal perturbation of chromosomal replication integrates biological regulatory factors with biophysical molecular events on the chromosome to reveal the local regulatory context of a gene. Whereas the TRIPs of many genes conform to a gene dosage-dependent pattern, others diverge in distinct ways, and this is shaped by factors such as intra-operon position and repression state. By revealing the underlying mechanistic drivers of gene expression heterogeneity, this work provides a quantitative, biophysical framework for modelling replication-dependent expression dynamics.

Transcription regulation by long non-coding RNAs: mechanisms and disease relevance.

Long non-coding RNAs (lncRNAs) outnumber protein-coding transcripts, but their functions remain largely unknown. In this Review, we discuss the emerging roles of lncRNAs in the control of gene transcription. Some of the best characterized lncRNAs have essential transcription cis-regulatory functions that cannot be easily accomplished by DNA-interacting transcription factors, such as XIST, which controls X-chromosome inactivation, or imprinted lncRNAs that direct allele-specific repression. A growing number of lncRNA transcription units, including CHASERR, PVT1 and HASTER (also known as HNF1A-AS1) act as transcription-stabilizing elements that fine-tune the activity of dosage-sensitive genes that encode transcription factors. Genetic experiments have shown that defects in such transcription stabilizers often cause severe phenotypes. Other lncRNAs, such as lincRNA-p21 (also known as Trp53cor1) and Maenli (Gm29348) contribute to local activation of gene transcription, whereas distinct lncRNAs influence gene transcription in trans. We discuss findings of lncRNAs that elicit a function through either activation of their transcription, transcript elongation and processing or the lncRNA molecule itself. We also discuss emerging evidence of lncRNA involvement in human diseases, and their potential as therapeutic targets.

A mutation rate model at the basepair resolution identifies the mutagenic effect of polymerase III transcription.

De novo mutations occur at substantially different rates depending on genomic location, sequence context and DNA strand. The success of methods to estimate selection intensity, infer demographic history and map rare disease genes, depends strongly on assumptions about the local mutation rate. Here we present Roulette, a genome-wide mutation rate model at basepair resolution that incorporates known determinants of local mutation rate. Roulette is shown to be more accurate than existing models. We use Roulette to refine the estimates of population growth within Europe by incorporating the full range of human mutation rates. The analysis of significant deviations from the model predictions revealed a tenfold increase in mutation rate in nearly all genes transcribed by polymerase III (Pol III), suggesting a new mutagenic mechanism. We also detected an elevated mutation rate within transcription factor binding sites restricted to sites actively used in testis and residing in promoters.

Elongation rate of RNA polymerase II affects pausing patterns across 3' UTRs.

Yeast mRNAs are polyadenylated at multiple sites in their 3' untranslated regions (3' UTRs), and poly(A) site usage is regulated by the rate of transcriptional elongation by RNA polymerase II (Pol II). Slow Pol II derivatives favor upstream poly(A) sites, and fast Pol II derivatives favor downstream poly(A) sites. Transcriptional elongation and polyadenylation are linked at the nucleotide level, presumably reflecting Pol II dwell time at each residue that influences the level of polyadenylation. Here, we investigate the effect of Pol II elongation rate on pausing patterns and the relationship between Pol II pause sites and poly(A) sites within 3' UTRs. Mutations that affect Pol II elongation rate alter sequence preferences at pause sites within 3' UTRs, and pausing preferences differ between 3' UTRs and coding regions. In addition, sequences immediately flanking the pause sites show preferences that are largely independent of Pol II speed. In wild-type cells, poly(A) sites are preferentially located < 50 nucleotides upstream from Pol II pause sites, but this spatial relationship is diminished in cells harboring Pol II speed mutants. Based on a random forest classifier, Pol II pause sites are modestly predicted by the distance to poly(A) sites but are better predicted by the chromatin landscape in Pol II speed derivatives. Transcriptional regulatory proteins can influence the relationship between Pol II pausing and polyadenylation but in a manner distinct from Pol II elongation rate derivatives. These results indicate a complex relationship between Pol II pausing and polyadenylation.

Transcriptional linkage analysis with in vivo AAV-Perturb-seq.

The ever-growing compendium of genetic variants associated with human pathologies demands new methods to study genotype-phenotype relationships in complex tissues in a high-throughput manner1,2. Here we introduce adeno-associated virus (AAV)-mediated direct in vivo single-cell CRISPR screening, termed AAV-Perturb-seq, a tuneable and broadly applicable method for transcriptional linkage analysis as well as high-throughput and high-resolution phenotyping of genetic perturbations in vivo. We applied AAV-Perturb-seq using gene editing and transcriptional inhibition to systematically dissect the phenotypic landscape underlying 22q11.2 deletion syndrome3,4 genes in the adult mouse brain prefrontal cortex. We identified three 22q11.2-linked genes involved in known and previously undescribed pathways orchestrating neuronal functions in vivo that explain approximately 40% of the transcriptional changes observed in a 22q11.2-deletion mouse model. Our findings suggest that the 22q11.2-deletion syndrome transcriptional phenotype found in mature neurons may in part be due to the broad dysregulation of a class of genes associated with disease susceptibility that are important for dysfunctional RNA processing and synaptic function. Our study establishes a flexible and scalable direct in vivo method to facilitate causal understanding of biological and disease mechanisms with potential applications to identify genetic interventions and therapeutic targets for treating disease.

High-throughput functional characterization of combinations of transcriptional activators and repressors.

Despite growing knowledge of the functions of individual human transcriptional effector domains, much less is understood about how multiple effector domains within the same protein combine to regulate gene expression. Here, we measure transcriptional activity for 8,400 effector domain combinations by recruiting them to reporter genes in human cells. In our assay, weak and moderate activation domains synergize to drive strong gene expression, whereas combining strong activators often results in weaker activation. In contrast, repressors combine linearly and produce full gene silencing, and repressor domains often overpower activation domains. We use this information to build a synthetic transcription factor whose function can be tuned between repression and activation independent of recruitment to target genes by using a small-molecule drug. Altogether, we outline the basic principles of how effector domains combine to regulate gene expression and demonstrate their value in building precise and flexible synthetic biology tools. A record of this paper's transparent peer review process is included in the supplemental information.

Butt-seq: a new method for facile profiling of transcription.

A wide range of sequencing methods has been developed to assess nascent RNA transcription and resolve the single-nucleotide position of RNA polymerase genome-wide. These techniques are often burdened with high input material requirements and lengthy protocols. We leveraged the template-switching properties of thermostable group II intron reverse transcriptase (TGIRT) and developed Butt-seq (bulk analysis of nascent transcript termini sequencing), which can produce libraries from purified nascent RNA in 6 h and from as few as 10,000 cells-an improvement of at least 10-fold over existing techniques. Butt-seq shows that inhibition of the superelongation complex (SEC) causes promoter-proximal pausing to move upstream in a fashion correlated with subnucleosomal fragments. To address transcriptional regulation in a tissue, Butt-seq was used to measure the circadian regulation of transcription from fly heads. All the results indicate that Butt-seq is a simple and powerful technique to analyze transcription at a high level of resolution.

What is an enhancer?

Tight control of the transcription process is essential for the correct spatial and temporal gene expression pattern during development and in homeostasis. Enhancers are at the core of correct transcriptional activation. The original definition of an enhancer is straightforward: a DNA sequence that activates transcription independent of orientation and direction. Dissection of numerous enhancer loci has shown that many enhancer-like elements might not conform to the original definition, suggesting that enhancers and enhancer-like elements might use multiple different mechanisms to contribute to transcriptional activation. Here, we review methodologies to identify enhancers and enhancer-like elements and discuss pitfalls and consequences for our understanding of transcriptional regulation.