Nuclear Pore Proteins Nup153 and Megator Define Transcriptionally Active Regions in the Genome
Transcriptional regulation is one of the most important processes for modulating gene expression. Though much of this control is attributed to transcription factors, histones, and associated enzymes, it is increasingly apparent that the spatial organization of chromosomes within the nucleus has a profound effect on transcriptional activity. Studies in yeast indicate that the nuclear pore complex might promote transcription by recruiting chromatin to the nuclear periphery. In higher eukaryotes, however, it is not known whether such regulation has global significance. Here we establish nucleoporins as a major class of global regulators for gene expression in Drosophila melanogaster. Using chromatin-immunoprecipitation combined with microarray hybridisation, we show that Nup153 and Megator (Mtor) bind to 25% of the genome in continuous domains extending 10 kb to 500 kb. These Nucleoporin-Associated Regions (NARs) are dominated by markers for active transcription, including high RNA polymerase II occupancy and histone H4K16 acetylation. RNAi–mediated knock-down of Nup153 alters the expression of ∼5,700 genes, with a pronounced down-regulatory effect within NARs. We find that nucleoporins play a central role in coordinating dosage compensation—an organism-wide process involving the doubling of expression of the male X chromosome. NARs are enriched on the male X chromosome and occupy 75% of this chromosome. Furthermore, Nup153-depletion abolishes the normal function of the male-specific dosage compensation complex. Finally, by extensive 3D imaging, we demonstrate that NARs contribute to gene expression control irrespective of their sub-nuclear localization. Therefore, we suggest that NAR–binding is used for chromosomal organization that enables gene expression control.
Published in the journal:
. PLoS Genet 6(2): e32767. doi:10.1371/journal.pgen.1000846
Category:
Research Article
doi:
https://doi.org/10.1371/journal.pgen.1000846
Summary
Transcriptional regulation is one of the most important processes for modulating gene expression. Though much of this control is attributed to transcription factors, histones, and associated enzymes, it is increasingly apparent that the spatial organization of chromosomes within the nucleus has a profound effect on transcriptional activity. Studies in yeast indicate that the nuclear pore complex might promote transcription by recruiting chromatin to the nuclear periphery. In higher eukaryotes, however, it is not known whether such regulation has global significance. Here we establish nucleoporins as a major class of global regulators for gene expression in Drosophila melanogaster. Using chromatin-immunoprecipitation combined with microarray hybridisation, we show that Nup153 and Megator (Mtor) bind to 25% of the genome in continuous domains extending 10 kb to 500 kb. These Nucleoporin-Associated Regions (NARs) are dominated by markers for active transcription, including high RNA polymerase II occupancy and histone H4K16 acetylation. RNAi–mediated knock-down of Nup153 alters the expression of ∼5,700 genes, with a pronounced down-regulatory effect within NARs. We find that nucleoporins play a central role in coordinating dosage compensation—an organism-wide process involving the doubling of expression of the male X chromosome. NARs are enriched on the male X chromosome and occupy 75% of this chromosome. Furthermore, Nup153-depletion abolishes the normal function of the male-specific dosage compensation complex. Finally, by extensive 3D imaging, we demonstrate that NARs contribute to gene expression control irrespective of their sub-nuclear localization. Therefore, we suggest that NAR–binding is used for chromosomal organization that enables gene expression control.
Introduction
The spatial organisation of DNA, both at the nucleotide and chromosomal levels, allows efficient storage of genetic information inside the nucleus. However, DNA-dependent processes such as transcription, require the chromosomal structure to be modified in order to allow access to this information.
The regulation of chromatin accessibility is an intensely studied subject [1],[2]. Molecular and genomic investigations have examined how nucleotide sequences and ATP-dependent chromatin-remodelling enzymes specify the locations for nucleosomal-binding, and how histone-modifying enzymes modulate the stability of histone-nucleic acid interactions. These enzymes are recruited to precise genomic loci with the aid of sequence-specific DNA-binding transcription factors. In turn, particular histone modifications influence transcription factor-binding to target sites on the genome, so controlling transcriptional initiation. Despite the importance of these cis- and trans-acting factors on the local chromosomal environment and the transcription of nearby genes, it has become increasingly clear that they explain just one level at which chromatin is regulated [3],[4].
The eukaryotic genome is spatially distributed in a highly organised manner, with entire chromosomal regions localising to well-defined sub-nuclear positions [5]. This organisation has a profound effect on chromatin accessibility and transcriptional activity on a genome-wide level [6]–[8]. For instance, chromosomal regions at the nuclear envelope tend to form closed heterochromatin, a structure that is generally indicative of transcriptional repression [9]. Genomic studies in Drosophila melanogaster and humans established that lamins—proteins lining the nuclear membrane [10]—are major contributors to sub-nuclear localisation and gene regulation [11],[12]. Comparisons of binding profiles with gene expression data and histone marker information showed that chromosomal regions containing dense lamin-binding were transcriptionally repressed.
Although the nuclear periphery has been primarily associated with repression, recent evidence has also suggested a role for membrane components in transcriptional activation [9], [13]–[16]. The nuclear pore complex is a large structure comprising about 30 protein subunits, and it is the primary channel through which macromolecules traverse the nuclear envelope [17]. Interestingly, investigations in Saccharomyces cerevisiae identified subunits of the nuclear pore complex that preferentially bound transcriptionally active genes [18]. Moreover, several target loci such as GAL2 and INO1 were found to relocate from the interior to the periphery upon activation [13], although there were exceptions to this behaviour [19]–[22]. Thus, it is becoming increasingly clear that nuclear periphery components can have both positive and negative influence on gene regulation.
Since there are differences in the composition of the nuclear envelope—such as the lack of lamins—it is important to also study the contribution of nuclear envelope components in gene regulation in higher organisms [9], [17], [23]–[25]. So far just one study has explored the global interactions of nucleoporin subunit Nup93 with human chromosomes 5, 7 and 16 [26]; the publication reported only a low density of binding sites, and their influence on gene regulation was inconclusive.
Recently, we revealed a biochemical association between nucleoporins and the dosage compensation apparatus in higher eukaryotes including humans [27]. In Drosophila, the Male Specific Lethal (MSL) complex offsets the imbalance in the number of sex chromosomes in males and females by doubling the expression of genes on the male X chromosome [28],[29]. By purifying enzymatically active MOF complexes, we identified interactions with the nucleoporins Nup153 and Megator (Mtor). Strikingly, depletion of either subunit resulted in the loss of dosage compensation in male cells. Therefore, our work suggested a vital role for nucleoporins in promoting transcriptional activation on a large-scale.
Here, we present the first genome-wide study of nucleoporin-binding in a higher eukaryote. Using chromatin immunoprecipitation followed by hybridisation to high resolution tiling microarrays, we show that Nup153 and Mtor interact with 25% of the Drosophila genome in large domains spanning 10–500 kb in size. These regions—which we term nucleoporin associated regions (NARs)—contain large numbers of highly expressed genes, and are enriched for markers of active transcription including RNA polymerase-binding and histone H4 lysine 16 acetylation. Additionally, we reveal a remarkably high density of NARs on the male X chromosome, which correlate extremely well with the binding pattern of the dosage compensation complex. Finally, we demonstrate that chromosomal regions bound by these nucleoporins are composed of peripheral as well as non-peripheral pools of these proteins but interestingly the X chromosomal target regions are preferentially localised closer to the nuclear periphery. In summary, we firmly establish nucleoporins as a major class of chromatin-binding proteins in higher eukaryotes, with a general role in transcriptional regulation and three-dimensional chromosomal organisation. Finally we show for the first time, the importance of nucleoporin-binding not only as a mechanism for transcriptional control, but also in maintaining a complex organism-level biological system namely dosage compensation.
Results
Nup153 and Mtor bind chromatin in a genome-wide fashion
We produced DNA-binding profiles for nuclear pore components Mtor and Nup153 in Drosophila male SL-2 and female KC cell lines using chromatin immunoprecipitation followed by hybridisation to Affymetrix tiling arrays [30],[31] (Figure 1). Raw data were processed as in Kind et al (2008) to minimise false-positive signals from aberrant array probes (Figure S1) [32].
The ChIP-chip profiles for the two proteins strongly correlate, indicating they bind to similar locations throughout the genome (r = 0.77 and 0.88 for SL-2 and KC cells respectively; Figure 1D, Figure S4). We confirmed the reproducibility of results by performing three biological replicates for each condition (r = 0.73), and we validated binding at 18 control genes by real-time PCR in triplicate (Figure S2).
Both Mtor and Nup153 exhibit extensive binding across the whole genome, and together they bind to 42% of the Drosophila genome (calculated as a fraction of base-pairs covered with two-fold cut-off). Thus nucleoporins represent a new class of global chromatin-binding proteins for higher eukaryotes.
Nucleoporin-binding occurs in large chromosomal domains
Visual inspection of the ChIP-chip profiles reveals that Nup153 and Mtor interact with the genome in a manner not observed for traditional transcription factors (Figure 1B and 1C) [33]. Instead of associating with discrete loci, nucleoporins bind extended chromosomal regions that alternate between domains of high-density binding with those of low occupancy.
In order to analyse the visual observations in a statistically rigorous fashion, we quantified binding that takes place within a 10 kb sliding window that was scanned along the genome (see Materials and Methods). Windows containing more than 70% binding (as a proportion of array probes with positive binding signal) were classified as Nucleoporin Associated Regions (NARs), and neighbouring windows reaching this threshold were grouped together as continuous NARs. The detection method is robust: the 70% threshold ensures that no NARs are found when binding sites are randomly distributed across the genome and we identify very similar sets of NARs for windows ranging 5 kb to 500 kb in size. Moreover, application of the domain-finding approach described by Guelen et al [11] returns over 80% agreement with our method (in terms of base-pairs classified as NARs).
There is considerable NAR-occurrence (Figure 1A–1C); in male SL-2 cells, a total of 1,384 NARs cover a quarter of the entire Drosophila genome (25Mb and 29Mb for Nup153 and Mtor respectively) and in female Kc cells 1,865 NARs occupy a similar proportion of the genome (33Mb and 35Mb for Nup153 and Mtor respectively; Figure S3). Most domains range in size from 10 kb to 100 kb, although some even extend to over 500 kb (Figure 1F, Figure S4). Most nucleotide positions within NARs are occupied by both Nup153 and Mtor. Moreover, even where the overlap is not perfect, NARs tend to occur in similar genomic loci (Figure 1E; Figure 1B chromosomal positions 560,000–600,000). Most importantly, NARs occur in gene-rich areas that encompass over 4,700 protein-coding genes whose activities might be affected by nucleoporin-binding.
Nucleoporin-binding demarcates actively transcribed chromosomal regions
A direct relationship between nucleoporin-binding and gene expression has not been established so far in higher eukaryotes. Therefore, we explored the impact of NARs on transcriptional regulation by examining the activity of genes encoded within these regions (Figure 2; Tables S1, S2).
We measured gene expression levels using Affymetrix GeneChips (see Materials and Methods). Using the present-absence calls defined by the MAS5.0 algorithm [34], we detected the expression of 6,478 and 6,219 genes in SL-2 and Kc cells respectively. These genes are preferentially located within NARs: 63% of genes inside NARs are expressed compared with just 40% outside, indicating a significantly elevated transcriptional activity in the former (p-value <2.2e−16). This observation is supported by data quantifying RNA polymerase II-occupancy (Figure 2; Tables S1, S2); by mapping publicly available ChIP-chip data [35], we find the Pol II-binding is highly enriched at the promoters of genes inside NARs compared with those outside (p-value <2.2e−16).
Recent publications demonstrated that histone modifications, MOF acetyltransferase- and lamin-binding are robust genome-wide indicators of transcriptional activity. In both SL-2 and KC cells, acetylated histone H4 lysine (H4K16Ac) and MOF-binding [32]—strong markers for active transcription—are extremely prominent within NARs (Figure 2; Tables S1, S2; p-value <2.2×10−16). In contrast, histone H3 lysine 27 tri-methylation [36] and lamin-binding [12]—markers of transcriptional repression—are enriched outside NARs (Figure 2, Figure S5; Tables S1, S2; p-value <2.2e−16).
Finally, we confirmed a causal link between nucleoporin-binding and transcriptional regulation by measuring gene expression levels following RNAi-mediated knock-down of Nup153 (Figure 2, Figure S7; Tables S1, S2). The depletion results in large and wide-spread transcriptional changes in cells collected after seven days: 5,684 genes −40% of Drosophila genes represented on the array—are differentially expressed in SL-2 cells (p-value <0.05). Moreover, there is a large enrichment of down-regulated genes within NARs (29% of all genes; 40% of ‘present’ genes) compared with non-NARs (19% of all genes; p-value <2.2e−16). We obtain similar enrichments for cells collected five days after RNAi-treatment, and also upon Mtor-depletion (data not shown). These observations strongly indicate that nucleoporin-binding promotes a high-level of transcriptional activity, which may be due to the formation of an open chromatin environment.
NARs are enriched on the male X chromosome
One of the most important manifestations of gene expression control in higher eukaryotes is dosage compensation for different number of sex chromosomes between the two sexes. In Drosophila—in which females have two X chromosomes but males possess only a single X—the dosage compensation complex offsets the imbalance in gene content by doubling the expression of the male X chromosome. Thus, the chromosome represents an outstanding example of an exceptionally highly transcribed genomic region.
In order to explore the association of Nup153 and Mtor with the dosage compensation complex further, we compared the patterns of nucleoporin-binding in male SL-2 and female Kc cells (Figure 1A, Figure 3A–3D, Figure S3). There is a dramatic difference between the two sexes: in females, NARs are evenly distributed throughout the entire genome with only a 1.2-fold difference in % NAR occupancy between chromosome X (7.4Mb and 33% for Nup153; 8.0Mb and 36% for Mtor) and autosomes (26.0Mb and 27% for Nup153; 27.1Mb and 28% for Mtor); but in males, NARs are overwhelmingly biased towards the X chromosome (14.9Mb and 67% for Nup153; 16.6Mb and 75% for Mtor) compared with the autosomes (9.7Mb and 10% for Nup153; 12.0Mb and 12% for Mtor) with a 6-fold difference in occupancy. Further, domains on the male X chromosome (median length = 62Kb, 94Kb for Nup153 and Mtor respectively) are much longer than those found on any other chromosomes (median length = 22Kb for Nup153 and Mtor in male autosomes, ∼35Kb for female autosomes and X chromosome).
Having established that the nucleoporins are enriched on the male X chromosome, we explored the association with the dosage compensation system further. Recently, we demonstrated that the members of the dosage compensation complex—MSL1, MSL3 and MOF—preferentially bind to the male X chromosome [32]. A comparison of this previously published dataset with our current analysis shows that NARs on the male X chromosome coincide very well with the binding sites of the dosage compensation complex (Figure 3E).
We also tested the effects of Nup153-depletion on MSL1 and MOF-binding to 10 known target loci using chromatin-immunoprecipitation followed by qPCR. X-chromosomal binding is severely reduced for both proteins (Figure 3F), and the additional binding to autosomal targets is lost for MOF (Figure S8). The effects are clearly specific to Nup153, as depleting another nucleoporin, Nup50 does not influence MSL1 and MOF-localisation and binding (Figure S9; data not shown). Moreover, the observations are not due to an effect on MSL protein concentrations or defects in the RNA export pathway [27]: we previously showed that MSL levels remain unaffected in Nup153 and Mtor-depleted cells; and impairment of the major export pathways through NFX1-depletion does not disrupt the localisation of the MSL complex to the X chromosomes.
Spatial localisation of NARs versus non–NARs in the nucleus
Although nucleoporins are primarily located at the nuclear periphery, some display dynamic association with the nuclear pore complex [37], and it remains unclear whether nucleoporin-chromatin interactions would affect transcription at the periphery or within the nucleoplasm. Therefore, we assessed the spatial localisation of different chromosomal regions within the nucleus using three-dimensional imaging of Fluorescence In Situ Hybridisation (3D-FISH) in male and female cells (Figure 4). We selected 26 chromosomal regions of average length 15–20 kb for analysis (Table S3), comprising 18 NAR (targets T1-18) and 8 non-NAR loci (targets N1-8). An independent lamin-bound locus (target L105) was used as a positive control representing a region previously shown to localise at the nuclear periphery [12].
First we checked the localisation of Nup153 and Mtor themselves (Figure S6). Immunostaining of SL-2 cells and salivary glands from male larvae confirm that both proteins predominantly reside in the nuclear periphery, although we also detected some staining within the nucleus. This is consistent with earlier reports that these proteins are dynamic components of the nuclear pore complex, with the capacity to shuttle between different sub-nuclear locations [25],[37].
Next, we used DAPI and lamin protein-immunostaining to assess the nuclear localisation of our target loci. We display a selection of images in Figure 4A: the lamin protein in green defines the nuclear boundary, the DAPI in blue the distribution of genomic DNA, and the FISH signal in red specifies the position of the target locus. In order to account for cell-to-cell variation in localisation that results from the dynamic behaviour of chromatin, we measured the distance between the FISH signal and nuclear boundary for a large number of samples (44<n<91). Size differences between nuclei were normalised by representing distances as a percentage of the nuclear radius. In Figure 4B, we show the expected distribution of distances for a simulated locus situated at the periphery; for a FISH signal with 30% radius, we find that most measurements lie between 0% and 30% of the distance to the centre of the nucleus. In contrast simulations for a signal positioned halfway between the periphery and the centre results in a distinct, more symmetrically shaped distribution, with most measurements falling between 20% and 60% of the distance to the centre (Figure 4C; Figures S10, S11; Videos S1, S2, S3, S4).
The lamin-bound L105 locus displays a distribution that is heavily skewed towards the periphery (Figure 4D); however the profile is broader than the simulation, signifying that the locus is present at the interior of the nucleus at least part of the time. On the other hand, target N2 resembles that of the non-peripheral simulation (Figure 4E), albeit with a broader distribution, which indicates that the locus predominantly resides in the interior. Since both loci are NAR-independent, they were assigned as in vivo controls representing peripheral and non-peripheral localisation.
Many NAR-target distributions show almost perfect overlap with L105, demonstrating that they are preferentially situated at the periphery (Figure 4F–4G; see Materials and Methods); interestingly however a subset of NAR loci displays distributions that are indicative of non-peripheral localisation (Figure 4I). For non-NARs, targets such as N1 display good overlap with the negative control N2 (Figure 4H), but some are found at the periphery. It is clear, therefore, that many targets regions tested here do not conform to the behaviour expected from NPC-binding.
In fact, we find that NARs from chromosome X tend to reside at the periphery (6 out of 10 targets; Table S4), whereas only a small number of autosomal NARs do so (1 out of 8; Table S4). This is reflected in the aggregate distributions, in which X-chromosomal loci display the characteristic skewed profiles compared with autosomal regions (Figure 4J–4L). Among non-NARs (Figure 4K–4M), autosomal loci are invariably non-peripheral, whereas the X chromosomal targets display a tendency for peripheral localisation; the positioning of the latter is probably influenced by neighbouring NARs as there is such a large amount of binding on the X chromosome. For comparison, peripheral localisation of the X chromosome is absent in female Kc cells (data not shown). Thus in striking contrast to prior expectations, we reveal that interior as well as peripheral populations of nucleoporins bind chromatin and mediate transcriptional activity at NARs. Furthermore, interactions with the X chromosome promotes peripheral localisation of the chromosome—most likely as a result of the overwhelming amount of binding in males—but this is generally not the case for autosomes.
Finally to confirm the influence of nucleoporins on localisation, we tested the effects of RNAi-mediated Nup153-knockdown for six loci: three peripheral X chromosomal NARs (T4, T5, T7), a non-peripheral X chromosomal NAR (T11), a non-peripheral autosomal NAR (T9) and the non-peripheral control (N2). For each we compared the distribution of Nup153-depleted samples against a mock EGFP RNAi-treatment (Figure 5, Figure S7). All three peripheral targets on the X chromosome displace to a more intra-nuclear position upon loss of Nup153 (Figure 5A–5C; p-value <0.05), but in contrast there was no significant change for any of the non-peripheral loci (Figure 5D–5F; p-value >0.05). These data suggest that the sub-nuclear positioning of peripheral NARs—specifically those on the male X—depends on the presence of Nup153, whereas the localisation of intra-nuclear loci is independent regardless of whether they are bound by nucleoporins.
Discussion
The classical view of transcriptional regulation describes the interplay of transcription factors, histones and associated enzymes with DNA in order to recruit the transcriptional machinery to the appropriate genomic loci. However, it has become increasingly clear that these interactions explain only one level at which gene expression is controlled. At a genome-wide level, the spatial organisation of chromosomes within the nucleus is increasingly considered to have a profound effect on chromatin structure and transcriptional activity [5]. In particular, studies in yeast indicate that members of the nuclear pore complex might promote transcription by recruiting chromatin to the nuclear periphery [14],[18]. However, the importance of such regulation in higher eukaryotes has remained unresolved [26].
In this study, we established conclusively that nucleoporins play a central role in mediating transcriptional regulation in a complex, multicellular organism. For the first time in any higher eukaryote, we generated a genome-wide profile of nucleoporin-binding; contrary to preliminary observations, binding is widespread, occurring across 40% of the genome. Thus, we reveal that nucleoporins—Nup153 and Mtor in particular—represent a major new class of global chromatin-binding proteins.
Intriguingly, these proteins interact with the genome differently to traditional transcription factors. Rather than associate with individual loci, nucleoporins bind continuous sections of chromosomes at very high density. Termed NARs, these regions extend up to 500kb in length and occupy 25% of the entire Drosophila genome. Moreover, NARs are functionally important as they demarcate regions of open chromatin and transcriptional activity, which is lost on depletion of Nup153. It is significant that the male X chromosome—a prime example for hyper-transcription—is almost entirely occupied by NARs. Therefore, we suggest that Nup153 and Mtor may stimulate transcription by promoting the formation of an open chromatin environment.
In dramatic contrast to expectations, nucleoporin-binding does not automatically lead to localisation at the nuclear periphery, though the male X chromosome is an exception in this regard. Since Nup153 and Mtor are known to be dynamic components of the nuclear pore complex, it appears likely that both peripheral and intra-nuclear pools of nucleoporins contribute to chromatin-binding. Given the dynamic nature of chromatin-localisation, it is also possible that NARs are located at the periphery in a very transient manner, and further developments in imaging techniques will help clarify this. Where NAR-formation and peripheral localisation do coincide however, Nup153 is necessary for sustained positioning.
Chromosomal domains have been implicated in the formation of three-dimensional sub-nuclear structures to coordinate the expression of otherwise distant loci [38] such as the human beta-globin genes [39],[40]. We speculate that NARs may indicate the genomic regions required for the assembly of these transcription factories on a very large scale. Within this context, the dynamic nature of Nup153 and Mtor is significant, as re-localisation of these proteins might allow a basis for global transcriptional control in response to cellular cues. Additionally, given the primary function of the nuclear pore complex in transporting macromolecules to and from the nucleus, Nup153 and Mtor may provide a means to couple transcriptional control with post-transcriptional events. We stress however that the mechanisms behind such processes are the subject of intense research activity and many controversies remain.
Finally, the special link with dosage compensation confirms the importance of nucleoporin-binding not only as a molecular mechanism for transcriptional control, but also in maintaining a complex, organism-level biological system.
Materials and Methods
ChIP–chip and qPCR analysis
Chromatin immunoprecipitation combined with microarray hybridisation (ChIP-chip), and qPCR experiments were performed as described previously in Kind et al [32]. Primer sequences are provided in Text S1.
Numerical data from Affymetrix Drosophila Tiling 2.0R Arrays (Dm35b_MR_v02) were processed as in Kind et al [32]. Briefly, array data were background corrected using GCRMA and quantile normalised [41]. Log2 (ChIP/input) ratios were calculated using the average from three replicate experiments. Log2 ratios were then smoothed by averaging the signal within a 500 bp window centred on each probe (Figure S1).
Identification of Nucleoporin Associated Regions (NARs)
Chromosomal regions with high densities of Nup153- and Mtor-binding were identified by sliding a 10 kb window along each chromosome, centred on the start position of each probe. NARs were defined as continuous chromosomal regions containing positive binding signal (ie, log2 ratio >0) for more than 70% of probes. We also implemented the two-stage domain-finding method described by Guelen et al [11]. Our method recovered at least 80% of all probes defined as domains by the Guelen approach.
RNAi on cultured cells
Nup153 and Nup50 were depleted as previously described in Mendjan et al [27]. Briefly, cells were incubated with dsRNA for five or seven days with a boost on day two. Cells were subsequently harvested for Western blot analysis, ChIP, gene expression profiling, or immunofluorescence experiments. Control experiments were performed using mock treatment (EGFP RNAi).
Gene expression profiling
Gene expression was measured using Affymetrix Drosophila2 GeneChips in triplicate for each condition. Data analysis was performed using publicly available packages in the BioConductor Software Suite [43]. Raw .CEL files were processed using RMA [44] and probe-sets were mapped to genes using the annotation from the Ensembl database (v41) [45].
In control (EGFP-treated) cells, expressed genes were identified as those outputting MAS5.0 ‘present’ cells in all three replicates [34]. For comparisons of Nup153-depleted and mock-treated cells, differentially expressed genes were determined using the Limma package [46]; p-values were corrected for multiple-testing using FDR [47] and a significance threshold of p-value<0.05 was selected.
Overlap of NARs with markers for transcriptional activity
We compared the overlap between NARs and genomic features. For ease of comparison, all data were mapped onto the Drosophila genome provided by the Ensembl database (v. 41) [45]. Accompanying each entry is the statistical significance of the difference in the amount of genomic feature found within NARs and non-NARs.
(i) Histone H4 lysine K16 acetylation (H4K16Ac; p<2.2e−16; t-test): processed ChIP-chip profiles obtained from Kind et al [32]. (ii) MOF-binding (p<2.2e−16; Fisher test): processed ChIP-chip profiles obtained from Kind et al [32]. (iii) RNA PolII-occupancy (p<2.2e−16; Fisher test): PolII-bound genes obtained from Muse et al [35]. For visualisation purposes in Figure 2, bound genes were represented as 1kb windows centred on the transcription start site. (iv) Gene density (p-value <2.2e−16; Wilcoxon test): number of genes as annotated by the Ensembl database within a 20kb sliding window with a 1 kb offset. (v) Expressed genes (p-value <2.2e−16; Fisher test): gene expression measured using Affymetrix Drosophila2 GeneChips as described above. (vi) Down-regulated genes upon Nup153-depletion (p-value <2.2e−16; Fisher test): differentially expressed genes in RNAi-treated cells compared with untreated cells as described above. (vii) Lamin-binding (p-value <2.2e−16; Fisher test): processed ChIP-chip data were obtained from Pickersgill et al [12]. Note that the study used low-resolution cDNA arrays, and therefore unlike the human study, the authors were unable to detect high-density lamin-associated domains. (viii) Histone H3 lysine 27 tri-methylation (H3K27me3; p-value <2.2e−16; Fisher test): processed ChIP-chip profiles obtained from Schwartz et al [36].
Fluorescent In Situ Hybridisation on cultured cells
DNA FISH on SL-2 cells was performed as previously described by Lanzuolo et al [48]. Briefly for DNA FISH 1×106 cells were centrifuged, re-suspended in 0.4 ml of medium and placed for 30 min at room temperature on a poly-lysine-coated slide (10 mm diameter). After rinsing with PBS, the cells were fixed with 4% paraformaldehyde in PBT (PBS, 0.1% Tween 20) for 10 min at room temperature. Cells were then washed three times with PBT, incubated for 1 h at room temperature with RNAse A (100 µg/ml in PBT). After rinsing with PBS, cells were incubated with 0.5% Triton in PBS for 10 min at room temperature. Cells were rinsed again with PBS and incubated with 20% glycerol in PBS for 30 min at room temperature. Cells were then frozen in liquid nitrogen, thawed at room temperature and soaked in 20% glycerol in PBS, repeatedly four times. After washing the cells again with PBS three times, they were incubated for 5 min in 0.1N HCl, briefly rinsed in 2XSSC twice, and stored in 50% formamide, 2XSSC, 10% dextransulphate, pH 7.0. Fluorescent probes were prepared with the FISH Tag DNA Kit (Invitrogen, Carlsbad, CA), dissolved in the hybridization mixture (50% deionized formamide, 2XSSC, 10% dextransulphate, salmon sperm DNA at 0.5 mg/ml), applied to cells and sealed under coverslips with rubber cement. Probe and cellular DNA were denatured simultaneously on a hot block at 78°C for 3 min. Hybridization was carried out in a humid atmosphere at 37°C for 1 d. After hybridization, slides were washed in 2XSSC three times for 5 min at 37°C, and in 0.1XSSC three times for 5 min at 45°C, rinsed in PBS twice and counter-stained with DAPI.
For immuno-FISH, the following procedure is added after washing with 0.1XSSC at 45°C. Wash twice with 2XSSC 5 min each at RT. Blocking with (TNT buffer; 0.1M Tris-HCl pH 7.5, 0.15M NaCl, 5% BSA) for 1 h at RT. Anti-lamin antibody is incubated for overnight at 4°C in TNT buffer, wash with wash buffer three times for 5 min. Second antibody is applied in TNT buffer for 2–3 h at RT, wash with wash buffer (0.1M Tris-HCl pH 7.5, 0.15M NaCl), including DAPI staining as described above. Cells were mounted on the glass slide with FluoromountG (Southern Biotech. Birmingham, AL). Three-dimensional image stacks were taken with Leica SP5 confocal microscope (Leica Microsystems, Exton, PA) using an x63 oil immersion objective with a numerical aperture of 1.4, and zoom 3.2±0.2.
To perform DNA FISH on target and non-target probes, approximately 15 kb region were chosen, except for the repeated sequence, in the genome and amplified by PCR from genomic DNA with 5–10 primers pairs, each covering around 0.5–3 kb. Primer sequences are available on request.
Image analysis of FISH localisation
To determine quantitatively the three-dimensional position of the FISH signal within the nucleus, we used the ImageJ software [49]. The nuclear envelope was initially defined by segmentation of the DAPI image using the automated Otsu thresholding algorithm. The boundary definition was then refined against the lamin-staining, flagging significant deviations between the two signals if necessary. Figure S10 shows a schematic diagram of the procedure. We also display a distribution of radii calculated for 62 nuclei, demonstrating that the DAPI and lamin signals provide very consistent definitions of the nuclear boundary. Segmented images were then stacked in order to recreate the three-dimensional nucleus.
Next we calculated the distances between the FISH signal and the nuclear boundary (Figure S10, S11, Videos S1, S2, S3, S4). The segmented three-dimensional images of the nucleus were converted into a three-dimensional distance map using the Local Thickness plug-in (http://www.optinav.com/Local_Thickness.htm). We thresholded the FISH images to identify voxels within the nucleus that corresponded to the FISH signal and we measured the distances between all such voxels and the closest point on the nuclear boundary. For each nucleus we calculated the mean distance, and then for each test locus, we use the set of mean distances for all nuclei to plot the distance distribution. Similar results were obtained when we used the centre of mass of the FISH signal as the reference point instead of the mean distances for individual voxels (data not shown).
In total, we examined 1,712 nuclei (35–91 samples for each target locus; total 1,172 nuclei for NAR; total 540 nuclei for non-NARs). For a given target, we compiled all distance measurements from all relevant nuclei to produce a distribution of distances as shown in Figure 4 and Table S4.
The lamin L105 and N2 non-NAR targets were selected as in vivo controls with representative distributions for peripheral and non-peripheral sub-nuclear localisation. We compared the localisation of each target locus by comparing its distance measurements against the L105 and N2 controls separately. Statistical significance was calculated using the Wilcoxon test, with a FDR-corrected threshold of p <0.05. Briefly, a non-significant p-value (ie p-value >0.05) compared with the L105 distribution is indicative of peripheral localisation, whereas a non-significant p-value (i.e. p-value >0.05) compared with the N2 distribution is indicative of non-peripheral localisation.
Accession numbers
Microarray data are available in the ArrayExpress databaset [42] under accession numbers E-MEXP-2523 (gene expression data) and E-MEXP-2525 (ChIP-chip data).
Supporting Information
Zdroje
1. KouzaridesT
2007 Chromatin modifications and their function. Cell 128 693 705
2. LiB
CareyM
WorkmanJL
2007 The role of chromatin during transcription. Cell 128 707 719
3. ClapierCR
CairnsBR
2009 The biology of chromatin remodelling complexes. Annu Rev Biochem 78 273 304
4. SuganumaT
WorkmanJL
2008 Crosstalk among histone modifications. Cell 135 604 607
5. LanctôtC
CheutinT
CremerM
CavalliG
CremerT
2007 Dynamic genome architecture in the nuclear space: regulation of gene expression in three dimensions. Nat Rev Genet 8 104 115
6. BrancoMR
PomboA
2006 Intermingling of chromosome territories in interphase suggests role in translocations and transcription-dependent associations. PLoS Biol 4 e138 doi:10.1371/journal.pbio.0040138
7. BrancoMR
PomboA
2007 Chromosome organization: new facts, new models. Trends Cell Biol 17 127 134
8. CremerT
CremerM
DietzelS
MüllerS
SoloveiI
2006 Chromosome territories - a functional nuclear landscape. Curr Opin Cell Biol 18 307 316
9. ShaklaiS
AmariglioN
RechaviG
SimonAJ
2007 Gene silencing at the nuclear periphery. FEBS J 274 1383 1392
10. DechatT
PfleghaarK
SenguptaK
ShimiT
ShumakerDK
2008 Nuclear lamins: major factors in the structural organization and function of the nucleus and chromatin. Genes Dev 22 832 853
11. GuelenL
PagieL
BrassetE
MeulemanW
FazaMB
2008 Domain organization of human chromosomes revealed by mapping of nuclear lamina interactions. Nature 453 948 951
12. PickersgillH
KalverdaB
de WitE
TalhoutW
FornerodM
2006 Characterization of the Drosophila melanogaster genome at the nuclear lamina. Nat Genet 38 1005 1014
13. BricknerJH
WalterP
2004 Gene recruitment of the activated INO1 locus to the nuclear membrane. PLoS Biol 2 e342 doi:10.1371/journal.pbio.0020342
14. BrownCR
SilverPA
2007 Transcriptional regulation at the nuclear pore complex. Curr Opin Genet Dev 17 100 106
15. CabalGG
GenovesioA
Rodriguez-NavarroS
ZimmerC
GadalO
2006 SAGA interacting factors confine sub-diffusion of transcribed genes to the nuclear envelope. Nature 441 770 773
16. TaddeiA
Van HouweG
HedigerF
KalckV
CubizollesF
2006 Nuclear pore association confers optimal expression levels for an inducible yeast gene. Nature 441 774 778
17. TranEJ
WenteSR
2006 Dynamic nuclear pore complexes: life on the edge. Cell 125 1041 1053
18. CasolariJM
BrownCR
KomiliS
WestJ
HieronymusH
2004 Genome-wide localization of the nuclear transport machinery couples transcriptional status and nuclear organization. Cell 117 427 439
19. ChambeyronS
BickmoreWA
2004 Chromatin decondensation and nuclear reorganization of the HoxB locus upon induction of transcription. Genes Dev 18 1119 1130
20. FerreiraJ
PaolellaG
RamosC
LamondAI
1997 Spatial organization of large-scale chromatin domains in the nucleus: a magnified view of single chromosome territories. J Cell Biol 139 1597 1610
21. GilbertDM
2001 Nuclear position leaves its mark on replication timing. J Cell Biol 152 F11 F15
22. ZinkD
AmaralMD
EnglmannA
LangS
ClarkeLA
2004 Transcription-dependent spatial arrangements of CFTR and adjacent genes in human cell nuclei. J Cell Biol 166 815 825
23. AkhtarA
GasserSM
2007 The nuclear envelope and transcriptional control. Nat Rev Genet 8 507 517
24. GoldmanRD
ShumakerDK
ErdosMR
ErikssonM
GoldmanAE
2004 Accumulation of mutant lamin A causes progressive changes in nuclear architecture in Hutchinson-Gilford progeria syndrome. Proc Natl Acad Sci U S A 101 8963 8968
25. PaddyMR
BelmontAS
SaumweberH
AgardDA
SedatJW
1990 Interphase nuclear envelope lamins form a discontinuous network that interacts with only a fraction of the chromatin in the nuclear periphery. Cell 62 89 106
26. BrownCR
KennedyCJ
DelmarVA
ForbesDJ
SilverPA
2008 Global histone acetylation induces functional genomic reorganization at mammalian nuclear pore complexes. Genes Dev 22 627 639
27. MendjanS
TaipaleM
KindJ
HolzH
GebhardtP
2006 Nuclear pore components are involved in the transcriptional regulation of dosage compensation in Drosophila. Mol Cell 21 811 823
28. LucchesiJC
KellyWG
PanningB
2005 Chromatin remodeling in dosage compensation. Annu Rev Genet 39 615 651
29. StraubT
BeckerPB
2007 Dosage compensation: the beginning and end of generalization. Nat Rev Genet 8 47 57
30. IyerVR
HorakCE
ScafeCS
BotsteinD
SnyderM
2001 Genomic binding sites of the yeast cell-cycle transcription factors SBF and MBF. Nature 409 533 538
31. RenB
RobertF
WyrickJJ
AparicioO
JenningsEG
2000 Genome-wide location and function of DNA binding proteins. Science 290 2306 2309
32. KindJ
VaquerizasJM
GebhardtP
GentzelM
LuscombeNM
2008 Genome-wide analysis reveals MOF as a key regulator of dosage compensation and gene expression in Drosophila. Cell 133 813 828
33. Farnham
PJ
2009 Insights from genomic profiling of transcription factors. Nat Rev Genet 10 605 616
34. HubbellE
LiuWM
MeiR
2002 Robust estimators for expression analysis. Bioinformatics 18 1585 1592
35. MuseGW
GilchristDA
NechaevS
ShahR
ParkerJS
2007 RNA polymerase is poised for activation across the genome. Nat Genet 39 1507 1511
36. SchwartzYB
KahnTG
NixDA
LiXY
BourgonR
2006 Genome-wide analysis of Polycomb targets in Drosophila melanogaster. Nat Genet 38 700 705
37. RabutG
DoyeV
EllenbergJ
2004 Mapping the dynamic organization of the nuclear pore complex inside single living cells. Nat Cell Biol 6 1114 1121
38. SutherlandH
BickmoreWA
2009 Transcription factories: gene expression in unions? Nat Rev Genet 10 457 466
39. OsborneCS
ChakalovaL
BrownKE
CarterD
HortonA
2004 Active genes dynamically colocalize to shared sites of ongoing transcription. Nat Genet 36 1065 1071
40. SimonisM
KlousP
SplinterE
MoshkinY
WillemsenR
2006 Nuclear organization of active and inactive chromatin domains uncovered by chromosome conformation capture-on-chip (4C). Nat Genet 38 1348 1354
41. WuZ
IrizarryRA
GentlemanR
Martinez-MurilloF
SpencerF
2004 A model-based background adjustment for oligonucleotide expression arrays. Journal of the American Statistical Association 99 909 917
42. ParkinsonH
KapusheskyM
KolesnikovN
RusticiG
ShojatalabM
2009 ArrayExpress update - from an archive of functional genomics experiments to the atlas of gene expression. Nucleic Acids Res 37 D868 D872
43. GentlemanRC
CareyVJ
BatesDM
BolstadB
DettlingM
2004 Bioconductor: open software development for computational biology and bioinformatics. Genome Biol 5 R80
44. IrizarryRA
BolstadBM
CollinF
CopeLM
HobbsB
2003 Summaries of Affymetrix GeneChip probe level data. Nucleic Acids Res 31 e15
45. FlicekP
AkenBL
BallesterB
BealK
BraginE
2010 Ensembl's 10th year. Nucleic Acids Res 38 D557 D562
46. SmythGK
2004 Linear models and empirical Bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol 3 Article3
47. BenjaminiY
HochbergY
1995 Controlling the false discovery rate: a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society Series B 57 289 300
48. LanzuoloC
RoureV
DekkerJ
BantigniesF
OrlandoV
2007 Polycomb response elements mediate the formation of chromosome higher-order structures in the bithorax complex. Nat Cell Biol 9 1167 1174
49. DoughertyR
KunzelmannKH
2007 Computing local thickness of 3D structures with ImageJ. Microscopy and Microanalysis 13 1678 1679
Štítky
Genetika Reprodukční medicínaČlánek vyšel v časopise
PLOS Genetics
2010 Číslo 2
- Primární hyperoxalurie – aktuální možnosti diagnostiky a léčby
- Srdeční frekvence embrya může být faktorem užitečným v předpovídání výsledku IVF
- Akutní intermitentní porfyrie
- Vztah užívání alkoholu a mužské fertility
- Šanci na úspěšný průběh těhotenství snižují nevhodné hladiny progesteronu vznikající při umělém oplodnění
Nejčtenější v tomto čísle
- Genome-Wide Association Study in Asian Populations Identifies Variants in and Associated with Systemic Lupus Erythematosus
- Nuclear Pore Proteins Nup153 and Megator Define Transcriptionally Active Regions in the Genome
- The Genetic Interpretation of Area under the ROC Curve in Genomic Profiling
- Nucleoporins and Transcription: New Connections, New Questions