- Research article
- Open Access
A biophysical model for transcription factories
BMC Biophysics volume 6, Article number: 2 (2013)
Transcription factories are nuclear domains where gene transcription takes placealthough the molecular basis for their formation and maintenance are unknown. In thisstudy, we explored how the properties of chromatin as a polymer may contribute to thestructure of transcription factories. We found that transcriptional active chromatincontains modifications like histone H4 acetylated at Lysine 16 (H4K16ac). Singlefibre analysis showed that this modification spans the entire body of the gene.Furthermore, H4K16ac genes cluster in regions up to 500 Kb alternating active andinactive chromatin. The introduction of H4K16ac in chromatin induces stiffness in thechromatin fibre. The result of this change in flexibility is that chromatin couldbehave like a multi-block copolymer with repetitions of stiff-flexible(active-inactive chromatin) components. Copolymers with such structure self-organizethrough spontaneous phase separation into microdomains. Consistent with such modelH4K16ac chromatin form foci that associates with nascent transcripts. We propose thattranscription factories are the result of the spontaneous concentration of H4K16acchromatin that are in proximity, mainly in cis.
Transcription in eukaryotes is organized in transcription factories (TFs), which arenuclear domains where several genes are grouped to be transcribed together [1, 2][3–5]. The current opinion is that the genes in a TF interact by a loopingmechanism [5–8]. It has been suggested that chromatin looping plays an important role incontrolling gene activity by bringing together promoters and enhancers or TFs . Some studies suggest that promoter-enhancer loops are maintained by theinteraction of proteins associated with these cis-regulatory elements . This interaction precedes chromatin activation, which is required for generelocation to the TF by an unknown mechanism . It has been proposed that TFs are maintained by depletion attraction forces(excluding volume effect) between RNA pol II molecules . However, experimental evidence has shown that genes remain at the factoryeven when active RNA pol II is not present . This makes the excluding volume model very improbable and suggests that thisstructure is not the result of transcription. Instead experimental evidence points tohistone acetylation as being responsible for loop formation . For these reasons we explored the possible contribution of chromatinacetylation in the formation of TFs.
Chromatin at the TF is decondensed  and contains active transcription marks like histone acetylation or H3K36me3 . Among all the possible Lysine residues that can be acetylated, H4K16Ac isvery special because it prevents the formation of compacted chromatin by inhibiting theinter-fibre interaction [15–18]. Moreover, H4K16 acetylation is associated with both active chromatin  and with the active transcription marker H3K4me3 [19–21].
Results and discussion
To confirm whether H4K16Ac is associated with active chromatin, we analysed thedistribution of H4K16Ac in the nucleus of the TFs of circulating lymphocytes. The TFswere visualised as sites of incorporation of Br-UTP into nascent RNA. TFs appeared asdiscrete foci distributed along the edge of condensed chromatin (Figure 1a) as previously described in other cell types [2, 22, 23]. H4K16Ac was scattered in foci overlapping or very close to these Br-RNAsites (Figure 1a). To study the extent and degree of thehyper-acetylated chromatin in individual transcription units (TUs), we deconstructed thenuclei of these cells by making chromatin spreads. Under these conditions active RNApolymerases and epigenetic modifications of the chromatin are preserved. This treatmentdisassembled nuclei and spread templates over a wide area. The DNA adopts a linearstructure with no visible nucleosomes  and 95% of active polymerases remain associated with the DNA . When we stained the chromatin spreads with antibodies against H4K16Ac, theyshowed almost continuous fluorescent tracks along the DNA fibres (Figure 1b). These H4K16Ac tracks corresponded with active chromatin, as demonstratedby co-localisation of H4K16Ac with the nascent transcripts that were labelled eitherin vivo or in vitro by using Bromo-Uridine (BrU) or BrUTPrespectively. The area covered by acetylated histones was larger than that stained bythe nascent transcripts (Figure 1c). This was to be expectedbecause histone acetylation extends over long stretches of genes, whilst only a few RNApol II molecules are ever found on a given gene . Nevertheless, to demonstrate that the distribution of active RNA pol IImolecules is not an artifact of non-natural nucleotide incorporation, we carried outchromatin spreads with cells where transcription was not labelled with BrU. Ourexperiments demonstrated a similar co-localisation of H4K16Ac with P-RNA pol II(hyper-phosphorylated Ser2) (Figure 1c).
The chromatin spreading technique allowed us to measure the length of H4K16Ac tracks.The distribution of H4K16Ac stretches showed a lognormal distribution with average sizeof ~15 Kb (Figure 1e). H4K16Ac tracks rarely appeared isolated,instead they tended to cluster, spanning several hundreds of Kb (348 ± 90; range235–530 Kb) (Figure 1d). The extension of the gaps betweentwo consecutive H4K16Ac tracks in the cluster showed a lognormal distribution with anaverage distance of ~30 Kb (Figure 1f). The analysis of thepolymerases loaded onto H4K16Ac tracks showed that not all the tracks were stained withBr-RNA or P-RNA pol II. The number of nascent transcripts or P-RNA pol II per track waslow (0.7 ± 1 transcripts/track and 0.8 ± 0.9 P-RNA pol II/track). This was inaccordance with our previous findings, suggesting that most of the TUs contain onemolecule of RNA pol II . The fact that some H4K16Ac tracks of chromatin were not associated to RNApol II or Br-RNA could reflect a temporal discrepancy between the transcription andacetylation processes of chromatin. Indeed, transcription by RNA pol II takes only a fewminutes [25–27] while deacetylation of active chromatin can take several hours , providing a molecular memory of recently-transcribed chromatin. On the otherhand, H4K16Ac tracks are not a special feature of lymphocytes as we were able to findthe same chromatin organisation in all the mammalian cell types tested including: Hela,Epstein Barr transformed lymphocytes, human lymphocytes, primary human fibroblasts,primary mouse fibroblasts and murine erythroleukemia cells (both differentiated andundifferentiated).
The clusters in all the different cell types analysed were identical with respect to thenumber of TUs (8 ± 2 TUs/Cluster), suggesting that co-linear active genes expressedat the same time, in agreement with the analysis of expression data using FDCP mix cells . The sliding window analysis (applying a window of 500 Kb) over the entiregenome showed that genes are active in clusters (Figure 1g), inaccordance with our chromatin spreads data. Moreover, our results are consistent withthe co-expression data after a Serial Analysis of Gene Expression where the cluster sizewas <500 Kb . From these data we can conclude that co-linear TUs are active at the sametime in the same cell.
How are these TUs organised in the cell nucleus?
Collinear active TUs are enriched in H4K16Ac which confers stiffness and inhibitsinter-fiber interaction [15–17]. In this way, chromatin appears as a multi-block copolymer with stiff andflexible monomers (rod-coil)n system, where the rod is the stiff activeTU. The multi-block copolymers function as amphiphiles whose components segregateinto domains to avoid unfavourable contact with each other. In these systems,complete phase separation is prevented by the covalent linkage between the components . The rod block does not have the same conformational entropy as the coilblock and this restricts homogeneous packaging. In consequence anisotropicinteractions occur between the stiff blocks ending in a liquid crystalline domainwhere the different TUs are aligned in a high order smectic phase  (Figure 2a). Multi-block copolymers can adopt manydifferent structures depending on the relative proportions of the rod and coilphases. For example, when the rod phase is lower than 20% the structures obtained aremicrospheres  (Figure 2b). A calculation of the amount of activechromatin in a Hela cell line gives a proportion of rods to coils of ~12%, which isconsistent with active chromatin separated in many microspheres. These spheres forH4K16Ac chromatin were observed in the cell nucleus of human lymphocytes (the shapefactor was 0.93 + 0.05). Microspheres are regularly distributed in artificialpolymers with regular coils and rods resulting from the repulsion forces of coilspushing in all directions. In the cell nucleus TFs are not regularly distributedbecause the sizes of genes and intergenic distances are not as regular as inartificial polymers. Moreover, H4K16Ac foci concentrated at the edge of condensedchromatin (Figure 1a). A possible reason for this discrepancyis that chromatin cohabits with the inter chromatin compartment (ICC), which iscomposed by RNPs and proteins. This results in a biphasic system where the ICC(inelastic phase) segregates from chromatin , creating an interphase between both components. Under these conditionsmicrospheres containing H4K16Ac may be pushed by the coil polymers to the interphasebetween chromatin and ICC (Figure 2c), resulting in thelocalisation along the edge of the chromatin as observed.
A prediction of the multi-block copolymer model is that microphase separation mustpersist as long as H4K16Ac is present in chromatin. In fact, H4K16Ac foci wereunperturbed by treatments like 2 M NaCl extraction, which disrupts chromatin;transcription inhibition by DRB (5,6-Dichlorobenzimidazole1-β-D-ribofuranoside), which reduces RNA pol II transcription by 98%; and heatshock (1h 45°C), which releases RNA pol II from the DNA  (Additional file 1: Figure S1). The only way todisrupt these foci was by formamide treatment, which works as a solvent for theelectrostatic self-assembled polymers (Additional file 2:Figure S2).
These experiments contradict the excluding volume model  and are in agreement with the multi-block microphase separation hypothesisproposed in this study.
The next question about the genes in a TF is where they come from. Several studieshave shown that genes in cis and in trans are able to interact in the same TF [5–7]. However, the analysis of chromatin spreads showed that collinear genesare active in the same cell at the same time. This guarantees that several H4K16Actracks are in close proximity. Therefore, most of these collinear TUs would probablyaggregate in the same microsphere, as occurs in similar situations with themulti-block copolymers . The experimental evidence from chromosome configuration capture analysissuggests that local chromatin is the primary source of interaction for any genomicloci . Nevertheless, we cannot exclude the possibility that some genes locatedfurther away in the same chromosome or in another chromosome can interact due toproximity or chromatin folding.
Finally, a remarkable feature of TFs is their constant size across species anddifferentiation stages . According to the multi-block copolymer model for chromatin organisation,the way to change the size of H4K16Ac foci (and consequently TFs) is by increasingthe number of active genes in a given region or by unrestricting the mobility of theactive chromatin. The latter has been reported in experiments using plasmids thatrendered larger TFs than the endogenous ones [35, 36].
In summary, we present evidence of the relationship between epigenetic marks and theTF structure. Our model proposes that active chromatin self-organises in the nucleusdue to the special physical properties of H4K16Ac modified chromatin. Therefore, ourmodel implies that chromatin becomes activated (H4K16Ac modified) before joining aTF. This is conceptually very different from current transcription factory model,which proposes that genes are targeted to TFs to “enhance production byconcentrating the relevant machines, resources, and expertise in one place” .
Materials and methods
Transcription in vivo and in vitro
For in vivo transcription, cells grown on coverslips were incubated in presence of2.5 mM BrU for several min.
For in vitro transcription, cells grown on coverslips were treated as described .
Cells (103 cells in 5 μl) were spotted onto a 22 × 50mm glass slide and 5 μl of lyses buffer were added (Lyses buffer: 1%sarkosyl, 25 U/ml ribonuclease inhibitor, 10 mM EDTA, and 100 mM Tris–HCl (pH7.4)). After 10 min at 20°C, the slide was tilted to allow the drop to run down.Samples were air-dried and fixed in 4% Paraformaldehyde for 10 min. Clusters weredefined as two o more hyper acetylated tracks in less than 100 Kb. For quantificationof clusters of hyper acetylated chromatin between 150 and 200 tracks of hyperacetylated chromatin were analysed.
After blocking for non-specific antibody binding, immunolabelling was carried out asdescribed . For detection of primary transcripts, we used mouse anti-IdU/BrdU (5mg/ml; Caltag Laboratories, Burlingame, CA). For detection of H4K16ac we usedantibodies raised in rabbit and mouse (Serotec, Kidlington, UK, Abcam). RNA pol IIhyperphosphorylated in Ser 2 was detected with H5 antibody (Covance). Secondaryantibodies donkey anti-mouse IgG or IgM tagged with Cy3 (1/200 dilution; JacksonImmunoResearch, Bar Harbor, ME) and donkey anti-rabbit IgG tagged with Alexa 488(1/200; prepared using a Molecular Probes kit, Inc., Eugene, OR). DNA staining wasperformed with 200 nM TOPRO-3 (Molecular Probes) for 5 min. Then coverslips weremounted on slides using Vectashield (Vector laboratories), and images were collectedusing a Radiance 2000 confocal microscope (Bio-Rad Laboratories, Hemel Hempstead,Herts, UK), Distances were measured using EasiVision software (Soft Imaging SystemsGmbH, Münster, Germany) and data exported to Excel (Microsoft) for analysis.
The degree of spreading of the chromatin was measured by hybridising the spreads witha fragment of DNA of 47.26 Kb; the spreading was 3.9 + 0.2 Kb/μm.
Microarrays and sliding window analysis
Mouse FCDP mix cells were used. cRNA synthesis and hybridisation to oligonucleotidearray were performed as described .
The sliding window analysis was performed by applying a window of 500 Kb over thechromosomes and moved at 5 Kb steps along a chromosome to know whether the genescontained in that window were more likely to be transcribed together than just bychance.
Iborra FJ, Jackson DA, Cook PR: The path of transcripts from extra-nucleolar synthetic sites to nuclear pores:transcripts in transit are concentrated in discrete structures containing SRproteins. Journal of cell science. 1998, 111: 2269-2282.
Iborra FJ, Pombo A, Jackson DA, Cook PR: Active RNA polymerases are localized within discrete transcription "factories' inhuman nuclei. Journal of cell science. 1996, 109: 1427-1436.
Iborra FJ, Escargueil AE, Kwek KY, Akoulitchev A, Cook PR: Molecular cross-talk between the transcription, translation, and nonsense-mediateddecay machineries. Journal of cell science. 2004, 117: 899-906. 10.1242/jcs.00933.
Iborra FJ: The path that RNA takes from the nucleus to the cytoplasm: a trip with somesurprises. Histochem Cell Biol. 2002, 118: 95-103.
Osborne CS, Chakalova L, Brown KE, Carter D, Horton A, Debrand E, Goyenechea B, Mitchell JA, Lopes S, Reik W, Fraser P: Active genes dynamically colocalize to shared sites of ongoing transcription. Nat Genet. 2004, 36: 1065-1071. 10.1038/ng1423.
Osborne CS, Chakalova L, Mitchell JA, Horton A, Wood AL, Bolland DJ, Corcoran AE, Fraser P: Myc dynamically and preferentially relocates to a transcription factory occupiedby Igh. PLoS Biol. 2007, 5: e192-10.1371/journal.pbio.0050192.
Chakalova L, Carter D, Debrand E, Goyenechea B, Horton A, Miles J, Osborne C, Fraser P: Developmental regulation of the beta-globin gene locus. Prog Mol Subcell Biol. 2005, 38: 183-206. 10.1007/3-540-27310-7_8.
Simonis M, Klous P, Splinter E, Moshkin Y, Willemsen R, de Wit E, van Steensel B, de Laat W: Nuclear organization of active and inactive chromatin domains uncovered bychromosome conformation capture-on-chip (4C). Nat Genet. 2006, 38: 1348-1354. 10.1038/ng1896.
Chambeyron S, Bickmore WA: Does looping and clustering in the nucleus regulate gene expression?. Curr Opin Cell Biol. 2004, 16: 256-262. 10.1016/j.ceb.2004.03.004.
Kim SI, Bultman SJ, Kiefer CM, Dean A, Bresnick EH: BRG1 Requirement for long-range interaction of a locus control region with adownstream promoter. Proc Natl Acad Sci U S A. 2009, 106: 2259-2264. 10.1073/pnas.0806420106.
Kim A, Dean A: Chromatin loop formation in the beta-globin locus and its role in globin genetranscription. Mol Cells. 2012,
Marenduzzo D, Micheletti C, Cook PR: Entropy-driven genome organization. Biophys J. 2006, 90: 3712-3721. 10.1529/biophysj.105.077685.
Mitchell JA, Fraser P: Transcription factories are nuclear subcompartments that remain in the absence oftranscription. Genes Dev. 2008, 22: 20-25. 10.1101/gad.454008.
Muller WG, Rieder D, Karpova TS, John S, Trajanoski Z, McNally JG: Organization of chromatin and histone modifications at a transcription site. J Cell Biol. 2007, 177: 957-967. 10.1083/jcb.200703157.
Shogren-Knaak M, Ishii H, Sun JM, Pazin MJ, Davie JR, Peterson CL: Histone H4-K16 acetylation controls chromatin structure and proteininteractions. Science. 2006, 311: 844-847. 10.1126/science.1124000.
Braunstein M, Sobel RE, Allis CD, Turner BM, Broach JR: Efficient transcriptional silencing in saccharomyces cerevisiae requires aheterochromatin histone acetylation pattern. Mol Cell Biol. 1996, 16: 4349-4356.
Grunstein M: Histone acetylation in chromatin structure and transcription. Nature. 1997, 389: 349-352. 10.1038/38664.
Robinson PJ, An W, Routh A, Martino F, Chapman L, Roeder RG, Rhodes D: 30 nm chromatin fibre decompaction requires both H4-K16 acetylation and linkerhistone eviction. J Mol Biol. 2008, 381: 816-825. 10.1016/j.jmb.2008.04.050.
Wang Z, Zang C, Rosenfeld JA, Schones DE, Barski A, Cuddapah S, Cui K, Roh TY, Peng W, Zhang MQ, Zhao K: Combinatorial patterns of histone acetylations and methylations in the humangenome. Nat Genet. 2008, 40: 897-903. 10.1038/ng.154.
Millar CB, Grunstein M: Genome-wide patterns of histone modifications in yeast. Nat Rev Mol Cell Biol. 2006, 7: 657-666. 10.1038/nrm1986.
O'Neill LP, VerMilyea MD, Turner BM: Epigenetic characterization of the early embryo with a chromatinimmunoprecipitation protocol applicable to small cell populations. Nat Genet. 2006, 38: 835-841. 10.1038/ng1820.
Jackson DA, Iborra FJ, Manders EM, Cook PR: Numbers and organization of RNA polymerases, nascent transcripts, andtranscription units in HeLa nuclei. Molecular biology of the cell. 1998, 9: 1523-1536.
Faro-Trindade I, Cook PR: A conserved organization of transcription during embryonic stem celldifferentiation and in cells with high C value. Mol Biol Cell. 2006, 17: 2910-2920. 10.1091/mbc.E05-11-1024.
Hawley DK, Roeder RG: Functional steps in transcription initiation and reinitiation from the major latepromoter in a HeLa nuclear extract. J Biol Chem. 1987, 262: 3452-3461.
Kimura H, Sugaya K, Cook PR: The transcription cycle of RNA polymerase II in living cells. The Journal of cell biology. 2002, 159: 777-782. 10.1083/jcb.200206019.
Hieda M, Winstanley H, Maini P, Iborra FJ, Cook PR: Different populations of RNA polymerase II in living mammalian cells. Chromosome research: an international journal on the molecular, supramolecularand evolutionary aspects of chromosome biology. 2005, 13: 135-144. 10.1007/s10577-005-7720-1.
das Neves RP, Jones NS, Andreu L, Gupta R, Enver T, Iborra FJ: Connecting variability in global transcription rate to mitochondrialvariability. PLoS Biol. 2010, 8: e1000560-10.1371/journal.pbio.1000560.
Kouskouti A, Talianidis I: Histone modifications defining active genes persist after transcriptional andmitotic inactivation. EMBO J. 2005, 24: 347-357. 10.1038/sj.emboj.7600516.
Bruno L, Hoffmann R, McBlane F, Brown J, Gupta R, Joshi C, Pearson S, Seidl T, Heyworth C, Enver T: Molecular signatures of self-renewal, differentiation, and lineage choice inmultipotential hemopoietic progenitor cells in vitro. Mol Cell Biol. 2004, 24: 741-756. 10.1128/MCB.24.2.741-756.2004.
Lercher MJ, Urrutia AO, Hurst LD: Clustering of housekeeping genes provides a unified model of gene order in thehuman genome. Nat Genet. 2002, 31: 180-183. 10.1038/ng887.
Olsen BDSR: Self-assembly of rod-coil block copolymers. Materials Science and Engineering R. 2008, 62: 37-66. 10.1016/j.mser.2008.04.001.
Matsen MW, Bates FS: Unifying weak- and strong-segregation block copolymer theories. Macromolecules. 1996, 29: 1091-1098. 10.1021/ma951138i.
Iborra FJ: Can visco-elastic phase separation, macromolecular crowding and colloidal physicsexplain nuclear organisation?. Theor Biol Med Model. 2007, 4: 15-10.1186/1742-4682-4-15.
Muthukumar M, Ober CK, Thomas EL: Competing interactions and levels of ordering in self-organizing polymericmaterials. Science. 1997, 277: 1225-1232. 10.1126/science.277.5330.1225.
Binnie A, Castelo-Branco P, Monks J, Proudfoot NJ: Homologous gene sequences mediate transcription-domain formation. J Cell Sci. 2006, 119: 3876-3887. 10.1242/jcs.03050.
Xu M, Cook PR: Similar active genes cluster in specialized transcription factories. J Cell Biol. 2008, 181: 615-623. 10.1083/jcb.200710053.
Papantonis A, Cook PR: Fixing the model for transcription: the DNA moves, not the polymerase. Transcription. 2011, 2: 41-44. 10.4161/trns.2.1.14275.
The authors would like to thank Elspeth McFarlane for helping with the English. To:Ministerio de Economia y Competitividad, (Spain) (Grant number: BFU2009-10792) andThe Medical Research Council (UK), for supporting this work. We also thankFundação Ciência e Tecnologia, Portugal for funding RPN.
AZC-H, RPN and JER performed some of the immunocytochemical experiments and acquireddata. SS and TE, performed the transcriptomic and the statistical analyses. CIparticipated in the draft the manuscript. VJB was in charge of the DNA fish experiments& FJI conceived of the study, and participated in its design and coordination andhelped to draft the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Figure S1: Stability of H4 K16Ac foci. Resistance of H4 K16Ac foci to various treatmentsthat disrupted transcription or chromatin structure. The aspect of H4 K16Acfoci did not change after DRB treatment (2h 150 μM) or heat shock(Hs) for 1h at 45°C. Both treatments led to the release of RNA pol II fromthe genes. These foci were also resistant to NaCl extraction (cellspermeabilised with 0.05% Triton X100 for 5 min in PBS at 4°C followed by10 min extraction with 2M NaCl for 10 min). The images were pseudo-coloured fordisplay. The bottom bar shows the scale of pseudo-colours used. (EPS 2 MB)
Additional file 2: Figure S2: H4 K16Ac foci are disassembled by formamide. The resistance of H4 K16Ac foci toformamide treatment. Cells were incubated for 5 min in PBS with differentconcentrations of formamide (0, 25, 50 and 100%) then fixed with 4%paraformaldehyde and immunolabelled with H4K16Ac antibodies. (a) TheH4K16Ac foci were disassembled by formamide treatment, as can be seen from thechange in the staining pattern, which is more diffuse and less intense than thecontrol. The images were pseudo-coloured for display. (b) Thedeconstruction of the foci was quantified by the change in the pixel intensityvariation coefficient (SD/mean). This analysis was performed by measuring themean intensity and the standard deviation (SD) of the H4K16Ac signal of thenuclear areas in at least 200 cells for each treatment. The images werepseudo-coloured for display. The bottom bar shows the scale of pseudo-coloursused. (EPS 2 MB)
About this article
Cite this article
Canals-Hamann, A.Z., das Neves, R.P., Reittie, J.E. et al. A biophysical model for transcription factories. BMC Biophys 6, 2 (2013) doi:10.1186/2046-1682-6-2
- Transcription Factories
- RNA pol II
- Nuclear organization