eRNA co-expression network uncovers TF dependency and convergent cooperativity

Enhancer RNAs (eRNAs) are non-coding RNAs produced by transcriptional enhancers that are highly correlated with their activity. Using a capped nascent RNA sequencing (PRO-cap) dataset in human lymphoblastoid cell lines across 67 individuals, we identified inter-individual variation in the expression of over 80 thousand transcribed transcriptional regulatory elements (tTREs), in both enhancers and promoters. Co-expression analysis of eRNAs from tTREs across individuals revealed how enhancers are associated with each other and with promoters. Mid- to long-range co-expression showed a distance-dependent decay that was modified by TF occupancy. In particular, we found a class of "bivalent" TFs, including Cohesin, that both facilitate and isolate the interaction between enhancers and/or promoters, depending on their topology. At short distances, we observed strand-specific correlations between nearby eRNAs in both convergent and divergent orientations. Our results support a cooperative model of convergent eRNAs, consistent with eRNAs facilitating adjacent enhancers rather than interfering with each other. Therefore, our approach to infer functional interactions from co-expression analyses provided novel insights into the principles of enhancer interactions as a function of distance, orientation, and binding landscapes of TFs.

Revealing eRNA interactions: TF dependency and convergent cooperativity

Enhancer RNAs (eRNAs) are non-coding RNAs produced from transcriptional enhancers that are highly correlated with their activities. Using capped nascent RNA sequencing (PRO-cap) dataset in human lymphoblastoid cell lines across individuals, we identified inter-individual variation of expression in over 80 thousand transcribed transcriptional regulatory elements (tTREs), in both enhancers and promoters. Co-expression analysis of eRNAs from tTREs across individuals revealed how enhancers interact with each other and with promoters. Mid-to-long range interactions showed distance-dependent decay, which was modified by TF occupancy. In particular, we found a class of 'bivalent' TFs, including Cohesin, which both facilitates and insulates the interaction between enhancers and/or promoters depending on the topology. In short ranges, we observed strand specific interactions between nearby eRNAs in both convergent or divergent orientations. Our finding supports a cooperative convergent eRNA model, which is compatible with eRNA remodeling neighboring enhancers rather than interfering with each other. Therefore, our approach to infer functional interactions from co-expression analyses provided novel insights into the principles of enhancer interactions depending on the distance, orientation, and the binding landscapes of TFs.

RANKL-responsive epigenetic mechanism reprograms macrophages into bone-resorbing osteoclasts

Monocyte/macrophage lineage cells are highly plastic and can differentiate into various cells under different environmental stimuli. Bone-resorbing osteoclasts are derived from the monocyte/macrophage lineage in response to receptor activator of NF-κB ligand (RANKL). However, the epigenetic signature contributing to the fate commitment of monocyte/macrophage lineage differentiation into human osteoclasts is largely unknown. In this study, we identified RANKL-responsive human osteoclast-specific superenhancers (SEs) and SE-associated enhancer RNAs (SE-eRNAs) by integrating data obtained from ChIP-seq, ATAC-seq, nuclear RNA-seq and PRO-seq analyses. RANKL induced the formation of 200 SEs, which are large clusters of enhancers, while suppressing 148 SEs in macrophages. RANKL-responsive SEs were strongly correlated with genes in the osteoclastogenic program and were selectively increased in human osteoclasts but marginally presented in osteoblasts, CD4+ T cells, and CD34+ cells. In addition to the major transcription factors identified in osteoclasts, we found that BATF binding motifs were highly enriched in RANKL-responsive SEs. The depletion of BATF1/3 inhibited RANKL-induced osteoclast differentiation. Furthermore, we found increased chromatin accessibility in SE regions, where RNA polymerase II was significantly recruited to induce the extragenic transcription of SE-eRNAs, in human osteoclasts. Knocking down SE-eRNAs in the vicinity of the NFATc1 gene diminished the expression of NFATc1, a major regulator of osteoclasts, and osteoclast differentiation. Inhibiting BET proteins suppressed the formation of some RANKL-responsive SEs and NFATc1-associated SEs, and the expression of SE-eRNA:NFATc1. Moreover, SE-eRNA:NFATc1 was highly expressed in the synovial macrophages of rheumatoid arthritis patients exhibiting high-osteoclastogenic potential. Our genome-wide analysis revealed RANKL-inducible SEs and SE-eRNAs as osteoclast-specific signatures, which may contribute to the development of osteoclast-specific therapeutic interventions.

Dynamic and widespread control of poly(A) tail length during macrophage activation

The poly(A) tail enhances translation and transcript stability, and tail length is under dynamic control during cell state transitions. Tail regulation plays essential roles in translational timing and fertilization in early development, but poly(A) tail dynamics have not been fully explored in post-embryonic systems. Here, we examined the landscape and impact of tail length control during macrophage activation. Upon activation, more than 1500 mRNAs, including proinflammatory genes, underwent distinctive changes in tail lengths. Increases in tail length correlated with mRNA levels regardless of transcriptional activity, and many mRNAs that underwent tail extension encode proteins necessary for immune function and post-transcriptional regulation. Strikingly, we found that ZFP36, whose protein product destabilizes target transcripts, undergoes tail extension. Our analyses indicate that many mRNAs undergoing tail lengthening are, in turn, degraded by elevated levels of ZFP36, constituting a post-transcriptional feedback loop that ensures transient regulation of transcripts integral to macrophage activation. Taken together, this study establishes the complexity, relevance, and widespread nature of poly(A) tail dynamics, and the resulting post-transcriptional regulation during macrophage activation.

Genome-Wide Identification of Polyadenylation Dynamics with TED-Seq

Polyadenylation and deadenylation of mRNA are major RNA modifications associated with nucleus-to-cytoplasm translocation, mRNA stability, translation efficiency, and mRNA decay pathways. Our current knowledge of polyadenylation and deadenylation has been expanded due to recent advances in transcriptome-wide poly(A) tail length assays. Whereas these methods measure poly(A) length by quantifying the adenine (A) base stretch at the 3' end of mRNA, we developed a more cost-efficient technique that does not rely on A-base counting, called tail-end-displacement sequencing (TED-seq). Through sequencing highly size-selected 3' RNA fragments including the poly(A) tail pieces, TED-seq provides accurate measure of transcriptome-wide poly(A)-tail lengths in high resolution, economically suitable for larger scale analysis under various biologically transitional contexts.

Single-Cell Transcriptome Analysis of Colon Cancer Cell Response to 5-Fluorouracil-Induced DNA Damage

DNA damage often induces heterogeneous cell-fate responses, such as cell-cycle arrest and apoptosis. Through single-cell RNA sequencing (scRNA-seq), we characterize the transcriptome response of cultured colon cancer cell lines to 5-fluorouracil (5FU)-induced DNA damage. After 5FU treatment, a single population of colon cancer cells adopts three distinct transcriptome phenotypes, which correspond to diversified cell-fate responses: apoptosis, cell-cycle checkpoint, and stress resistance. Although some genes are regulated uniformly across all groups of cells, many genes showed group-specific expression patterns mediating DNA damage responses specific to the corresponding cell fate. Some of these observations are reproduced at the protein level by flow cytometry and are replicated in cells treated with other 5FU-unrelated genotoxic drugs, camptothecin and etoposide. This work provides a resource for understanding heterogeneous DNA damage responses involving fractional killing and chemoresistance, which are among the major challenges in current cancer chemotherapy.

TED-Seq Identifies the Dynamics of Poly(A) Length during ER Stress

Post-transcriptional RNA processing is a core mechanism of gene expression control in cell stress response. The poly(A) tail influences mRNA translation and stability, but it is unclear whether there are global roles of poly(A)-tail lengths in cell stress. To address this, we developed tail-end displacement sequencing (TED-seq) for an efficient transcriptome-wide profiling of poly(A) lengths and applied it to endoplasmic reticulum (ER) stress in human cells. ER stress induced increases in the poly(A) lengths of certain mRNAs, including known ER stress regulators, XBP1, DDIT3, and HSPA5. Importantly, the mRNAs with increased poly(A) lengths are both translationally de-repressed and stabilized. Furthermore, mRNAs in stress-induced RNA granules have shorter poly(A) tails than in the cytoplasm, supporting the view that RNA processing is compartmentalized. In conclusion, TED-seq reveals that poly(A) length is dynamically regulated upon ER stress, with potential consequences for both translation and mRNA turnover.

m6A enhances the phase separation potential of mRNA

N⁶-methyladenosine (m⁶A) is the most prevalent modified nucleotide in mRNA^1,2, with ~25% of mRNAs containing at least one m⁶A. Methylation of mRNA to form m⁶A is required for diverse cellular and physiological processes³. Although the presence of m⁶A in an mRNA can affect its fate in different ways, it is unclear how m⁶A directs this process and why the effects of m⁶A can vary in different cellular contexts. Here we show that the cytosolic m⁶A-binding proteins, YTHDF1–3, undergo liquid-liquid phase separation (LLPS) in vitro and in cells. This LLPS is markedly enhanced by mRNAs that contain multiple, but not single, m⁶A residues. Polymethylated mRNAs act as a multivalent scaffold for binding YTHDF proteins, juxtaposing their low-complexity domains, leading to phase separation. The resulting mRNA-YTHDF complexes then partition into different endogenous phase-separated compartments, such as P-bodies, stress granules, or neuronal RNA granules. m⁶A-mRNA is subject to compartment-specific regulation, including reduced mRNA stability and translation. These studies reveal that the number and distribution of m⁶A sites in cellular mRNAs can regulate and influence the composition of the phase-separated transcriptome. Additionally, these findings indicate that the cellular properties of m⁶A-modified mRNAs are governed by liquid-liquid phase separation principles.

Chromatin run-on and sequencing maps the transcriptional regulatory landscape of glioblastoma multiforme

The human genome encodes a variety of poorly understood RNA species that remain challenging to identify using existing genomic tools. We developed chromatin run-on and sequencing (ChRO-seq) to map the location of RNA polymerase for almost any input sample, including samples with degraded RNA that are intractable to RNA sequencing. We used ChRO-seq to map nascent transcription in primary human glioblastoma (GBM) brain tumors. Enhancers identified in primary GBMs resemble open chromatin in the normal human brain. Rare enhancers that are activated in malignant tissue drive regulatory programs similar to the developing nervous system. We identified enhancers that regulate groups of genes that are characteristic of each known GBM subtype and transcription factors that drive them. Finally we discovered a core group of transcription factors that control the expression of genes associated with clinical outcomes. This study characterizes the transcriptional landscape of GBM and introduces ChRO-seq as a method to map regulatory programs that contribute to complex diseases.

Aspm knockout ferret reveals an evolutionary mechanism governing cerebral cortical size

The human cerebral cortex is distinguished by its large size and abundant gyrification, or folding, yet the evolutionary mechanisms driving cortical size and structure are unknown. While genes essential for cortical developmental expansion have been identified from the genetics of human primary microcephaly (“small head”, associated with reduced brain size and intellectual disability)¹, studies of these genes in mice, whose smooth cortex is one thousand times smaller than that of humans, have provided limited insight. Mutations of abnormal spindle-like microcephaly-associated (ASPM), the most common recessive microcephaly gene, reduce cortical volume by ≥50% in humans^2–4, but have little effect in mice^5–9, likely reflecting evolutionarily divergent functions of ASPM^10,11. We used genome editing to create a germline knockout (KO) of Aspm in the ferret (Mustela putorius furo), a species with a larger, gyrified cortex and greater neural progenitor cell (NPC) diversity^12–14 than mice, and closer Aspm protein sequence homology to human. Aspm KO ferrets exhibit severe microcephaly (25–40% decreases in brain weight), reflecting reduced cortical surface area without significant change in cortical thickness, as in human patients^3,4, suggesting loss of “cortical units”. The mutant ferret fetal cortex displays a massive premature displacement of ventricular radial glial cells (VRG) to the outer subventricular zone (OSVZ), where many resemble outer radial glia (ORG), an NPC subtype essentially absent in mice and implicated in cerebral cortical expansion in primates^12–16. These data suggest an evolutionary mechanism whereby Aspm regulates cortical expansion by controlling the affinity of VRG for the ventricular surface, thus modulating the ratio of VRG, the most undifferentiated cell type, to ORG, a more differentiated progenitor.

Systematic Characterization of Stress-Induced RNA Granulation

Upon stress, cytoplasmic mRNA is sequestered to insoluble ribonucleoprotein (RNP) granules, such as the stress granule (SG). Partially due to the belief that translationally suppressed mRNAs are recruited to SGs in bulk, stress-induced dynamic redistribution of mRNA has not been thoroughly characterized. Here, we report that endoplasmic reticulum (ER) stress targets only a small subset of translationally suppressed mRNAs into the insoluble RNP granule fraction (RG). This subset, characterized by extended length and adenylate-uridylate (AU)-rich motifs, is highly enriched with genes critical for cell survival and proliferation. This pattern of RG targeting was conserved for two other stress types, heat shock and arsenite toxicity, which induce distinct responses in the total cytoplasmic transcriptome. Nevertheless, stress-specific RG-targeting motifs, such as guanylate-cytidylate (GC)-rich motifs in heat shock, were also identified. Previously underappreciated, transcriptome profiling in the RG may contribute to understanding human diseases associated with RNP dysfunction, such as cancer and neurodegeneration.

RNA-DNA sequence differences in Saccharomyces cerevisiae

Alterations of RNA sequences and structures, such as those from editing and alternative splicing, result in two or more RNA transcripts from a DNA template. It was thought that in yeast, RNA editing only occurs in tRNAs. Here, we found that Saccharomyces cerevisiae have all 12 types of RNA–DNA sequence differences (RDDs) in the mRNA. We showed these sequence differences are propagated to proteins, as we identified peptides encoded by the RNA sequences in addition to those by the DNA sequences at RDD sites. RDDs are significantly enriched at regions with R-loops. A screen of yeast mutants showed that RDD formation is affected by mutations in genes regulating R-loops. Loss-of-function mutations in ribonuclease H, senataxin, and topoisomerase I that resolve RNA–DNA hybrids lead to increases in RDD frequency. Our results demonstrate that RDD is a conserved process that diversifies transcriptomes and proteomes and provide a mechanistic link between R-loops and RDDs.

High-Resolution Mapping of RNA Polymerases Identifies Mechanisms of Sensitivity and Resistance to BET Inhibitors in t(8;21) AML

Bromodomain and extra-terminal domain (BET) family inhibitors offer an approach to treating hematological malignancies. We used precision nuclear run-on transcription sequencing (PRO-seq) to create high-resolution maps of active RNA polymerases across the genome in t(8;21) acute myeloid leukemia (AML), as these polymerases are exceptionally sensitive to BET inhibitors. PRO-seq identified over 1,400 genes showing impaired release of promoter-proximal paused RNA polymerases, including the stem cell factor receptor tyrosine kinase KIT that is mutated in t(8;21) AML. PRO-seq also identified an enhancer 3' to KIT. Chromosome conformation capture confirmed contacts between this enhancer and the KIT promoter, while CRISPRi-mediated repression of this enhancer impaired cell growth. PRO-seq also identified microRNAs, including MIR29C and MIR29B2, that target the anti-apoptotic factor MCL1 and were repressed by BET inhibitors. MCL1 protein was upregulated, and inhibition of BET proteins sensitized t(8:21)-containing cells to MCL1 inhibition, suggesting a potential mechanism of resistance to BET-inhibitor-induced cell death.

Base-pair-resolution genome-wide mapping of active RNA polymerases using precision nuclear run-on (PRO-seq)

We provide a protocol for precision nuclear run-on sequencing (PRO-seq) and its variant, PRO-cap, which map the location of active RNA polymerases (PRO-seq) or transcription start sites (TSSs) (PRO-cap) genome-wide at high resolution. The density of RNA polymerases at a particular genomic locus directly reflects the level of nascent transcription at that region. Nuclei are isolated from cells and, under nuclear run-on conditions, transcriptionally engaged RNA polymerases incorporate one or, at most, a few biotin-labeled nucleotide triphosphates (biotin-NTPs) into the 3' end of nascent RNA. The biotin-labeled nascent RNA is used to prepare sequencing libraries, which are sequenced from the 3' end to provide high-resolution positional information for the RNA polymerases. PRO-seq provides much higher sensitivity than ChIP-seq, and it generates a much larger fraction of usable sequence reads than ChIP-seq or NET-seq (native elongating transcript sequencing). Similarly to NET-seq, PRO-seq maps the RNA polymerase at up to base-pair resolution with strand specificity, but unlike NET-seq it does not require immunoprecipitation. With the protocol provided here, PRO-seq (or PRO-cap) libraries for high-throughput sequencing can be generated in 4-5 working days. The method has been applied to human, mouse, Drosophila melanogaster and Caenorhabditis elegans cells and, with slight modifications, to yeast.

Chromatin Architecture of the Pitx2 Locus Requires CTCF- and Pitx2-Dependent Asymmetry that Mirrors Embryonic Gut Laterality

Expression of Pitx2 on the left side of the embryo patterns left-right (LR) organs including the dorsal mesentery (DM), whose asymmetric cell behavior directs gut looping. Despite the importance of organ laterality, chromatin-level regulation of Pitx2 remains undefined. Here, we show that genes immediately neighboring Pitx2 in chicken and mouse, including a long noncoding RNA (Pitx2 locus-asymmetric regulated RNA or Playrr), are expressed on the right side and repressed by Pitx2. CRISPR/Cas9 genome editing of Playrr, 3D fluorescent in situ hybridization (FISH), and variations of chromatin conformation capture (3C) demonstrate that mutual antagonism between Pitx2 and Playrr is coordinated by asymmetric chromatin interactions dependent on Pitx2 and CTCF. We demonstrate that transcriptional and morphological asymmetries driving gut looping are mirrored by chromatin architectural asymmetries at the Pitx2 locus. We propose a model whereby Pitx2 auto-regulation directs chromatin topology to coordinate LR transcription of this locus essential for LR organogenesis.

Genome-wide dynamics of Pol II elongation and its interplay with promoter proximal pausing, chromatin, and exons

Production of mRNA depends critically on the rate of RNA polymerase II (Pol II) elongation. To dissect Pol II dynamics in mouse ES cells, we inhibited Pol II transcription at either initiation or promoter-proximal pause escape with Triptolide or Flavopiridol, and tracked Pol II kinetically using GRO-seq. Both inhibitors block transcription of more than 95% of genes, showing that pause escape, like initiation, is a ubiquitous and crucial step within the transcription cycle. Moreover, paused Pol II is relatively stable, as evidenced from half-life measurements at ∼3200 genes. Finally, tracking the progression of Pol II after drug treatment establishes Pol II elongation rates at over 1000 genes. Notably, Pol II accelerates dramatically while transcribing through genes, but slows at exons. Furthermore, intergenic variance in elongation rates is substantial, and is influenced by a positive effect of H3K79me2 and negative effects of exon density and CG content within genes.DOI: http://dx.doi.org/10.7554/eLife.02407.001.

The Hmr and Lhr hybrid incompatibility genes suppress a broad range of heterochromatic repeats

Hybrid incompatibilities (HIs) cause reproductive isolation between species and thus contribute to speciation. Several HI genes encode adaptively evolving proteins that localize to or interact with heterochromatin, suggesting that HIs may result from co-evolution with rapidly evolving heterochromatic DNA. Little is known, however, about the intraspecific function of these HI genes, the specific sequences they interact with, or the evolutionary forces that drive their divergence. The genes Hmr and Lhr genetically interact to cause hybrid lethality between Drosophila melanogaster and D. simulans, yet mutations in both genes are viable. Here, we report that Hmr and Lhr encode proteins that form a heterochromatic complex with Heterochromatin Protein 1 (HP1a). Using RNA-Seq analyses we discovered that Hmr and Lhr are required to repress transcripts from satellite DNAs and many families of transposable elements (TEs). By comparing Hmr and Lhr function between D. melanogaster and D. simulans we identify several satellite DNAs and TEs that are differentially regulated between the species. Hmr and Lhr mutations also cause massive overexpression of telomeric TEs and significant telomere lengthening. Hmr and Lhr therefore regulate three types of heterochromatic sequences that are responsible for the significant differences in genome size and structure between D. melanogaster and D. simulans and have high potential to cause genetic conflicts with host fitness. We further find that many TEs are overexpressed in hybrids but that those specifically mis-expressed in lethal hybrids do not closely correlate with Hmr function. Our results therefore argue that adaptive divergence of heterochromatin proteins in response to repetitive DNAs is an important underlying force driving the evolution of hybrid incompatibility genes, but that hybrid lethality likely results from novel epistatic genetic interactions that are distinct to the hybrid background.

Sister species capable of mating often produce hybrids that are sterile or die during development. This reproductive isolation is caused by incompatibilities between the two sister species' genomes. Some hybrid incompatibilities involve genes that encode rapidly evolving proteins that localize to heterochromatin. Heterochromatin is largely made up of highly repetitive transposable elements and satellite DNAs. It has been hypothesized that rapid changes in heterochromatic DNA drives the changes in these HI genes and thus the evolution of reproductive isolation. In support of this model, we show that two rapidly evolving HI proteins, Lhr and Hmr, which reproductively isolate the fruit fly sister species D. melanogaster and D. simulans, repress transposable elements and satellite DNAs. These proteins also help regulate the length of the atypical Drosophila telomeres, which are themselves made of domesticated transposable elements. Our data suggest that these proteins are part of the adaptive machinery that allows the host to respond to changes and increases in heterochromatin and to maintain the activity of genes located within or adjacent to heterochromatin.

RNA-DNA differences are generated in human cells within seconds after RNA exits polymerase II

RNA sequences are expected to be identical to their corresponding DNA sequences. Here, we found all 12 types of RNA-DNA sequence differences (RDDs) in nascent RNA. Our results show that RDDs begin to occur in RNA chains ~55 nt from the RNA polymerase II (Pol II) active site. These RDDs occur so soon after transcription that they are incompatible with known deaminase-mediated RNA-editing mechanisms. Moreover, the 55 nt delay in appearance indicates that they do not arise during RNA synthesis by Pol II or as a direct consequence of modified base incorporation. Preliminary data suggest that RDD and R-loop formations may be coupled. These findings identify sequence substitution as an early step in cotranscriptional RNA processing.

Kinetics of promoter Pol II on Hsp70 reveal stable pausing and key insights into its regulation

The kinetics with which promoter-proximal paused RNA polymerase II (Pol II) undergoes premature termination versus productive elongation is central to understanding underlying mechanisms of metazoan transcription regulation. To assess the fate of Pol II quantitatively, we tracked photoactivatable GFP-tagged Pol II at uninduced Hsp70 on polytene chromosomes and showed that Pol II is stably paused with a half-life of 5 min. Biochemical analysis of short nascent RNA from Hsp70 reveals that this half-life is determined by two comparable rates of productive elongation and premature termination of paused Pol II. Importantly, heat shock dramatically increases elongating Pol II without decreasing termination, indicating that regulation acts at the step of paused Pol II entry to productive elongation.

Defining NELF-E RNA binding in HIV-1 and promoter-proximal pause regions

The four-subunit Negative Elongation Factor (NELF) is a major regulator of RNA Polymerase II (Pol II) pausing. The subunit NELF-E contains a conserved RNA Recognition Motif (RRM) and is proposed to facilitate Poll II pausing through its association with nascent transcribed RNA. However, conflicting ideas have emerged for the function of its RNA binding activity. Here, we use in vitro selection strategies and quantitative biochemistry to identify and characterize the consensus NELF-E binding element (NBE) that is required for sequence specific RNA recognition (NBE: CUGAGGA(U) for Drosophila). An NBE-like element is present within the loop region of the transactivation-response element (TAR) of HIV-1 RNA, a known regulatory target of human NELF-E. The NBE is required for high affinity binding, as opposed to the lower stem of TAR, as previously claimed. We also identify a non-conserved region within the RRM that contributes to the RNA recognition of Drosophila NELF-E. To understand the broader functional relevance of NBEs, we analyzed promoter-proximal regions genome-wide in Drosophila and show that the NBE is enriched +20 to +30 nucleotides downstream of the transcription start site. Consistent with the role of NELF in pausing, we observe a significant increase in NBEs among paused genes compared to non-paused genes. In addition to these observations, SELEX with nuclear run-on RNA enrich for NBE-like sequences. Together, these results describe the RNA binding behavior of NELF-E and supports a biological role for NELF-E in promoter-proximal pausing of both HIV-1 and cellular genes.

RNA polymerase II (Pol II) is a molecular machine that is responsible for transcribing all protein coding genes in the eukaryotic genome. Transcription by Pol II is a highly regulated process consisting of several rate-limiting steps. During transcription elongation, a number of transcription factors are essential to modulate Pol II activity. One of these factors is the Negative Elongation Factor (NELF), and it plays a major role in promoter-proximal pausing, a widespread phenomenon during early transcription elongation. NELF-E, a protein subunit of the NELF complex contains a conserved RNA binding domain that is thought to regulate transcription through its interaction with newly transcribed RNA made by Pol II. However, the function of the RNA binding activity of NELF-E remains unresolved due to prior conflicting studies. Here, we clarify the RNA binding properties of NELF-E and provide insight into how this protein might facilitate promoter-proximal pausing of Pol II in transcription. Moreover, we identify the precise region of NELF-E binding in one of its known regulatory targets, HIV-1. Taken together, the results presented indicate a dynamic interplay between NELF and specific RNA sequences around the promoter pause region to modulate early transcription elongation.

Control of transcriptional elongation

Elongation is becoming increasingly recognized as a critical step in eukaryotic transcriptional regulation. Although traditional genetic and biochemical studies have identified major players of transcriptional elongation, our understanding of the importance and roles of these factors is evolving rapidly through the recent advances in genome-wide and single-molecule technologies. Here, we focus on how elongation can modulate the transcriptional outcome through the rate-liming step of RNA polymerase II (Pol II) pausing near promoters and how the participating factors were identified. Among the factors we describe are the pausing factors--NELF (negative elongation factor) and DSIF (DRB sensitivity-inducing factor)--and P-TEFb (positive elongation factor b), which is the key player in pause release. We also describe the high-resolution view of Pol II pausing and propose nonexclusive models for how pausing is achieved. We then discuss Pol II elongation through the bodies of genes and the roles of FACT and SPT6, factors that allow Pol II to move through nucleosomes.

Precise maps of RNA polymerase reveal how promoters direct initiation and pausing

Transcription regulation occurs frequently through promoter-associated pausing of RNA polymerase II (Pol II). We developed a precision nuclear run-on and sequencing (PRO-seq) assay to map the genome-wide distribution of transcriptionally engaged Pol II at base pair resolution. Pol II accumulates immediately downstream of promoters, at intron-exon junctions that are efficiently used for splicing, and over 3' polyadenylation sites. Focused analyses of promoters reveal that pausing is not fixed relative to initiation sites, nor is it specified directly by the position of a particular core promoter element or the first nucleosome. Core promoter elements function beyond initiation, and when optimally positioned they act collectively to dictate the position and strength of pausing. This "complex interaction" model was tested with insertional mutagenesis of the Drosophila Hsp70 core promoter.

Genome-wide control of RNA polymerase II activity by cohesin

Cohesin is a well-known mediator of sister chromatid cohesion, but it also influences gene expression and development. These non-canonical roles of cohesin are not well understood, but are vital: gene expression and development are altered by modest changes in cohesin function that do not disrupt chromatid cohesion. To clarify cohesin's roles in transcription, we measured how cohesin controls RNA polymerase II (Pol II) activity by genome-wide chromatin immunoprecipitation and precision global run-on sequencing. On average, cohesin-binding genes have more transcriptionally active Pol II and promoter-proximal Pol II pausing than non-binding genes, and are more efficient, producing higher steady state levels of mRNA per transcribing Pol II complex. Cohesin depletion frequently decreases gene body transcription but increases pausing at cohesin-binding genes, indicating that cohesin often facilitates transition of paused Pol II to elongation. In many cases, this likely reflects a role for cohesin in transcriptional enhancer function. Strikingly, more than 95% of predicted extragenic enhancers bind cohesin, and cohesin depletion can reduce their association with Pol II, indicating that cohesin facilitates enhancer-promoter contact. Cohesin depletion decreases the levels of transcriptionally engaged Pol II at the promoters of most genes that don't bind cohesin, suggesting that cohesin controls expression of one or more broadly acting general transcription factors. The multiple transcriptional roles of cohesin revealed by these studies likely underlie the growth and developmental deficits caused by minor changes in cohesin activity.

Defining the status of RNA polymerase at promoters

Recent genome-wide studies in metazoans have shown that RNA polymerase II (Pol II) accumulates to high densities on many promoters at a rate-limited step in transcription. However, the status of this Pol II remains an area of debate. Here, we compare quantitative outputs of a global run-on sequencing assay and chromatin immunoprecipitation sequencing assays and demonstrate that the majority of the Pol II on Drosophila promoters is transcriptionally engaged; very little exists in a preinitiation or arrested complex. These promoter-proximal polymerases are inhibited from further elongation by detergent-sensitive factors, and knockdown of negative elongation factor, NELF, reduces their levels. These results not only solidify the notion that pausing occurs at most promoters, but demonstrate that it is the major rate-limiting step in early transcription at these promoters. Finally, the divergent elongation complexes seen at mammalian promoters are far less prevalent in Drosophila, and this specificity in orientation correlates with directional core promoter elements, which are abundant in Drosophila.

First-principles study of confinement effects on the Raman spectra of Si nanocrystals

The Raman spectra of Si nanocrystals are studied as a function of nanocrystal diameter using pseudopotential density functional theory and the Placzek approximation. Our calculations reproduce the redshift and broadening of the optical Raman peak with decreasing nanocrystal size, and calculated peak frequencies show good agreement with experimental values. We also find that a surface induced softening of vibrational modes is largely responsible for the Raman redshift, with relaxation of momentum conservation playing only a minor role.

Chemical inhibitors destabilize HuR binding to the AU-rich element of TNF-alpha mRNA

Hu protein R (HuR) binds to the AU-rich element (ARE) in the 3UTR to stabilize TNF-alpha mRNA. Here, we identified chemical inhibitors of the interaction between HuR and the ARE of TNF-alpha mRNA using RNA electrophoretic mobility gel shift assay (EMSA) and filter binding assay. Of 179 chemicals screened, we identified three with a half-maximal inhibitory concentration (IC(50)) below 10 microM. The IC(50) of quercetin, b-40, and b-41 were 1.4, 0.38, and 6.21 microM, respectively, for binding of HuR protein to TNF-alpha mRNA. Quercetin and b-40 did not inhibit binding of tristetraprolin to the ARE of TNF-alpha mRNA. When LPS-treated RAW264.7 cells were treated with quercetin and b-40, we observed decreased stability of TNF-alpha mRNA and decreased levels of secreted TNF-alpha. From these results, we could find inhibitors for the TNF-alpha mRNA stability, which might be used advantageously for both the study for post-transcriptional regulation and the discovery of new anti-inflammation drugs.

Duodenal polypectomy is a rare cause of acute pancreatitis

Endoscopic duodenal polypectomy is a routine procedure particularly useful for obtaining histological diagnosis but it is not without serious complications. This is a case report of severe necrotising pancreatitis after duodenal polypectomy. We suggest that experienced endoscopists should carry out polypectomies and that clear guidelines for the management of duodenal polyps are required. Patients undergoing endoscopic duodenal polypectomies should be placed at the beginning of the endoscopy list and observed for at least 4 h.

Arginine-genotype interactions and immune status

The effects of arginine on selective immune responses were investigated in a high arginine-requiring (HA) and low arginine-requiring (LA) strain of chickens. Female chickens from these strains were fed diet containing a nutritionally inadequate level of arginine (0.53% arginine diet) or a surfeit level of arginine (1.53% arginine diet) for 2 weeks. Compared to LA chickens, HA chickens showed a higher feed efficiency, body weight gain, and relative thymus and spleen weights with L-arginine supplementation (p < 0.05). In both HA and LA chickens, a deficiency of arginine significantly decreased the delayed-type hypersensitivity response (p < 0.05) and nitric oxide (NO) production from macrophages. Chickens of the HA strain had higher NO production than those of LA strain with E. coli lipopolysaccharide (LPS) activation. This study indicates that dietary arginine concentration influences the immune status of chickens and that strains that differ in arginine requirements for growth may differ in their arginine needs for immune function.

Interrelationship of apoptosis, mutation, and cell proliferation in N-methyl-N'-nitro-N-nitrosoguanidine (MNNG)-induced medaka carcinogenesis model

The present study examined the interrelationship of GSH depletion, apoptosis, mutation, and cell proliferation following carcinogen exposure. Medaka (Oryzias latipes) were investigated following a 28 day, three times/week pulse exposure to N-methyl-N'-nitro-N-nitrosoguanidine (MNNG). Fish (5 weeks old) were exposed to MNNG at concentrations of 0, 0.5, and 1 mg l(-1) and reared for 3, 5 and 7 more months after the last day of exposure. GSH levels were decreased in the higher concentration groups and longer-reared groups. Flow cytometric analysis revealed that fish from the groups reared 3 and 5 months showed active apoptotic changes in the dose- and time-dependent manner, but the group reared 7 months had fewer apoptotic, rather showed more necrotic and carcinogenic alterations. Mutational responses were detected by an arbitrarily primed polymerase chain reaction (AP-PCR) fingerprinting method using whole body DNA samples as templates and pBR primer. A mutational change was expressed by a loss or gain of a band. There was a time-dependent mutational change, but no distinctive concentration-dependent one. A band from normal fish sample that disappeared after treatment of MNNG was excised and sequenced. The band had an 869 base pair-long sequence, however, there was no putative protein-coding region based on an analysis by DNAsis. Spindle cell sarcomas invading muscle were detected on the whole body sections from three of ten fish examined, and immunohistochemical analysis with PCNA showed that tumor cells were actively proliferating. However, terminal deoxynucleotidyl transferase (TdT) assay showed that tumored fish still had active apoptotic cell changes in the tissues without tumor. This study shows not only the interrelationship of GSH depletion, apoptosis, mutation and cell proliferation, but also indicates that medaka is appropriate as a fish model for research on the passage of carcinogenesis.

Mutagenicity of recombinant antihemophilic factor (GC-gamma AHF)

This study was carried out to evaluate the mutagenic potential of recombinant antihemophilic factor VIII (GC-gamma AHF). Salmonella typhimurium (S. typhimurium) reversion assay with/without histidine moiety, chromosomal aberration assay on Chinese hamster lung (CHL) fibroblast cells and in vivo micronucleus assay using mouse bone marrow cells and supravital micronucleus assay using peripheral blood were performed. GC-gamma AHF containing histidine did show inconsistent and irregular mutagenic effects on S. typhimurium TA98, TA100, TA1535 and TA1537 both in the absence and presence of the metabolic activation system, however, GC-gamma AHF without histidine showed no mutagenic effects regardless of the metabolic activation system, thus suggesting that the histidine moiety in GC-gamma AHF might cause inconsistent mutagenic effect. Also GC-gamma AHF did not increase the number of cells having structural or numerical chromosome aberration in the cytogenetic test. In classical and supravital micronucleus assay, no significant increases were observed in the occurrence of micronucleated polychromatic erythrocytes and micronucleated peripheral lymphocytes in male ICR mice. These results strongly indicate that GC-gamma AHF has no genetic toxicity under these experimental conditions.

Influence of dietary arginine concentration on lymphoid organ growth in chickens

In vivo effects of graded dietary levels of arginine on the body and lymphoid organs were investigated using Cornell K strain chickens of the B15/B15 haplotype. Two-week-old birds were fed an arginine-deficient basal diet (0.53% arginine) supplemented with additional arginine (up to 1.0% L-arginine to the diet). At four weeks of age, body weight, lymphoid organ weight, and concentrations of amino acids in plasma were measured. Arginine supplementation produced significant increases in plasma arginine (from 200 nM in chicks fed the basal diet to 2,000 nM in chicks receiving the 1.5% arginine diet) and ornithine concentrations (from 17 nM in chicks fed the basal diet to 500 nM in chicks receiving the 1.5% arginine diet). The arginine-deficient diet reduced body weight gain (P < 0.0001) and thymus, spleen, and bursa of Fabricius weights (P < 0.05). In contrast to the bursa weight, the thymus and spleen weights, as percentages of body weight, were also decreased (P < 0.05). This study suggests that arginine markedly influences lymphoid organ development, with a more pronounced effect on the thymus and spleen than on the bursa of Fabricius.