Computer-assisted search for sites of nuclear matrix attachment

Large-scale organization of ribosomal DNA chromatin is regulated by Tip5

The DNase I accessibility and chromatin organization of genes within the nucleus do correlate to their transcriptional activity. Here, we show that both serum starvation and overexpression of Tip5, a key regulator of ribosomal RNA gene (rDNA) repression, dictate DNase I accessibility, facilitate the association of rDNA with the nuclear matrix and thus regulate large-scale rDNA chromatin organization. Tip5 contains four AT-hooks and a TAM (Tip5/ARBP/MBD) domain, which were proposed to bind matrix-attachment regions (MARs) of the genome. Remarkably, the TAM domain of Tip5 functions as nucleolar localization and nuclear matrix targeting module, whereas AT-hooks do not mediate association with the nuclear matrix, but they are required for nucleolar targeting. These findings suggest a dual role for Tip5's AT-hooks and TAM domain, targeting the nucleolus and anchoring to the nuclear matrix, and suggest a function for Tip5 in the regulation of higher-order rDNA chromatin structure.

Cis-regulation of microRNA expression by scaffold/matrix-attachment regions

microRNAs (miRNAs) spatio-temporally modulate gene expression; however, very little is known about the regulation of their expression. Here, we hypothesized that the well-known cis-regulatory elements of gene expression, scaffold/matrix-attachment regions (MARs) could modulate miRNA expression. Accordingly, we found MARs to be enriched in the upstream regions of miRNA genes. To determine their role in cell type-specific expression of miRNAs, we examined four individual miRNAs (let-7b, miR-17, miR-93 and miR-221) and the miR-17-92 cluster, known to be overexpressed in neuroblastoma. Our results show that MARs indeed define the cell-specific expression of these miRNAs by tethering the chromatin to nuclear matrix. This is brought about by cell type-specific binding of HMG I/Y protein to MARs that then promotes the local acetylation of histones, serving as boundary elements for gene activation. The binding, chromatin tethering and gene activation by HMG I/Y was not observed in fibroblast control cells but were restricted to neuroblastoma cells. This study implies that the association of MAR binding proteins to MARs could dictate the tissue/context specific regulation of miRNA genes by serving as a boundary element signaling the transcriptional activation.

Decondensing the protamine domain for transcription

Potentiation is the transition from higher-order, transcriptionally silent chromatin to a less condensed state requisite to accommodating the molecular elements required for transcription. To examine the underlying mechanism of potentiation an approximately 13.7-kb mouse protamine domain of increased nuclease sensitivity flanked by 5' and 3' nuclear matrix attachment regions was defined. The potentiated DNase I-sensitive region is formed at the pachytene spermatocyte stage with the recruitment to the nuclear matrix of a large approximately 9.6-kb region just upstream of the domain. Attachment is then specified in the transcribing round spermatid, recapitulating the organization of the human cluster. In comparison to other modifiers that have no effect, i.e., histone methylation, HP1, and SATB1, topoisomerase engages nuclear matrix binding as minor marks of histone acetylation appear. Reorganization is marked by specific sites of topoisomerase II activity that are initially detected in leptotene-zygotene spermatocytes just preceding the formation of the DNase I-sensitive domain. This has provided a likely model of the events initiating potentiation, i.e., the opening of a chromatin domain.

DNA bending in the replication zone of the C3 DNA puff amplicon of Rhynchosciara americana (Diptera: Sciaridae)

Intrinsic bent DNA sites were identified in the 4289 bp segment encompassing the replication zone which directs DNA amplification and transcription of the C3-22 gene of Rhynchosciara americana. Restriction fragments showed reduced electrophoretic mobility in polyacrylamide gels. The 2D modeling of the 3D DNA path and the ENDS ratio values obtained from the dinucleotide wedge model of Trifonov revealed the presence of four major bent sites, positioned at nucleotides -6753, -5433, -5133 and -4757. Sequence analysis showed that these bends are composed of 2-6 bp dA.dT tracts in phase with the DNA helical repeat. The circular permutation analysis permitted the verification that the fragments containing the bending sites promote curvature in other sequence contexts. Computer analyses of the 4289 bp sequence revealed low helical stability (DeltaG values), negative roll angles indicating a narrow minor groove and a putative matrix attachment region. The data presented in this paper add to information about the structural features involved in this amplified segment.

Correlations between scaffold/matrix attachment region (S/MAR) binding activity and DNA duplex destabilization energy

Scaffold or matrix-attachment regions (S/MARs) are thought to be involved in the organization of eukaryotic chromosomes and in the regulation of several DNA functions. Their characteristics are conserved between plants and humans, and a variety of biological activities have been associated with them. The identification of S/MARs within genomic sequences has proved to be unexpectedly difficult, as they do not appear to have consensus sequences or sequence motifs associated with them. We have shown that S/MARs do share a characteristic structural property, they have a markedly high predicted propensity to undergo strand separation when placed under negative superhelical tension. This result agrees with experimental observations, that S/MARs contain base-unpairing regions (BURs). Here, we perform a quantitative evaluation of the association between the ease of stress-induced DNA duplex destabilization (SIDD) and S/MAR binding activity. We first use synthetic oligomers to investigate how the arrangement of localized unpairing elements within a base-unpairing region affects S/MAR binding. The organizational properties found in this way are applied to the investigation of correlations between specific measures of stress-induced duplex destabilization and the binding properties of naturally occurring S/MARs. For this purpose, we analyze S/MAR and non-S/MAR elements that have been derived from the human genome or from the tobacco genome. We find that S/MARs exhibit long regions of extensive destabilization. Moreover, quantitative measures of the SIDD attributes of these fragments calculated under uniform conditions are found to correlate very highly (r2>0.8) with their experimentally measured S/MAR-binding strengths. These results suggest that duplex destabilization may be involved in the mechanisms by which S/MARs function. They suggest also that SIDD properties may be incorporated into an improved computational strategy to search genomic DNA sequences for sites having the necessary attributes to function as S/MARs, and even to estimate their relative binding strengths.

Matrix attachment region (MAR) properties and abnormal expansion of AT island minisatellites in FRA16B fragile sites in leukemic CEM cells

AT-rich minisatellites (AT islands) are sites of genomic instability in cancer cells and targets for extremely lethal AT-specific drugs, such as bizelesin. Here we investigated the AT islands in the FRA16B fragile site region for their possible roles in the organization of DNA on the nuclear matrix. The FRA16B AT island nominally spans approximately 3 kb of mostly >90% A/T DNA. In silico analysis indicates that this domain exhibits characteristics of nuclear matrix attachment regions (MARs): an exceptionally intense computed 'MAR potential' and profound duplex destabilization and flexibility. FRA16B repeats specifically bind to isolated nuclear matrices, which indicates their in vitro MAR function. This binding is several-fold greater than that of a known MAR in the c-myc gene. AT islands in fragile sites FRA16B and FRA16D are significantly more abundant in CEM cells that are hypersensitive to bizelesin compared to normal WI-38 cells. FRA16B overabundance in CEM is due to an approximately 10-fold expansion of FRA16B repeats. The expanded FRA16B minisatellites in CEM cells preferentially localize to the nuclear matrix-associated DNA indicating their in vivo MAR function. The unexpanded repeats in WI-38 cells localize to the loop DNA. The c-myc MAR is also matrix-associated in CEM cells while localizing to loop DNA in WI-38 cells. These results are the first to demonstrate that AT islands in fragile sites can function as MARs both in vitro and in vivo. The ability of FRA16B-mediated MAR sites to rearrange depending on the repeat expansion status could be relevant to both genomic instability of cancer cells and their sensitivity to AT-island targeting drugs.

Structure-specific DNA-binding proteins as the foundation for three-dimensional chromatin organization

Any functions of tandem repetitive sequences need proteins that specifically bind to them. Telomere-binding TRF2/MTBP attaches telomeres to the nuclear envelope in interphase due to its rod-domain-like motif. Interphase nuclei organized as a number of sponge-like ruffly round chromosome territories that could be rotated from outside. SAF-A/hnRNP-U and p68-helicase are proteins suitable to do that. Their location in the interchromosome territory space, ATPase domains, and the ability to be bound by satellite DNAs (satDNA) make them part of the wires used to help chromosome territory rotates. In case of active transcription p68-helicase can be involved in the formation of local "gene expression matrices" and due to its satDNA-binding specificity cause the rearrangement of the local chromosome territory. The marks of chromatin rearrangement, which have to be heritable, could be provided by SAF-A/hnRNP-U. During telophase unfolding the proper chromatin arrangement is restored according to these marks. The structural specificity of both proteins to the satDNAs provides a regulative but relatively stable mode of binding. The structural specificity of protein binding could help to find the "magic" centromeric sequence. With future investigations of proteins with the structural specificity of binding during early embryogenesis, when heterochromatin formation goes on, the molecular mechanisms of the "gene gating" hypothesis (Blobel, 1985) will be confirmed.

Functional analysis of two matrix attachment region (MAR) elements in transgenic maize plants

Matrix attachment regions (MARs) are binding sites for nuclear scaffold proteins in vitro, and are proposed to mediate the attachment of chromatin to the nuclear scaffold in vivo. Previous reports suggest that MAR elements may stabilize transgene expression. Here, we tested the effects of two maize MAR elements (P-MAR from the P1-rr gene, and Adh1-MAR from the adh1 gene) on the expression of a gusA reporter gene driven by three different promoters: the maize p1 gene promoter, a wheat peroxidase (WP) gene promoter, or a synthetic promoter (Rsyn7). The inclusion of P-MAR or Adh1-MAR on P::GUS transgene constructs did not reduce variation in the levels of GUS activity among independent transformation events, nor among the progeny derived from each event. The Adh1-MAR element did not affect GUS expression driven by the WP promoter, but did modify the spatial pattern of expression of the Rsyn7::GUS transgene. These results indicate that, in transgenic maize plants, the effects of MAR elements can vary significantly depending upon the promoter used to drive the transgene.

Chromatin and CD4, CD8A and CD8B gene expression during thymic differentiation

The regulation of gene expression during thymocyte development provides an ideal experimental system to study lineage-commitment processes. In particular, expression of the CD4, CD8A and CD8B genes seems to correlate well with the cell-fate decisions that are taken by thymocytes, and elucidating the molecular mechanisms that underlie the differential expression of these genes could reveal key events in differentiation processes. Here, we review examples of how gene cis elements (such as promoters, enhancers and locus control regions) and trans elements (such as transcription factors, chromatin-remodelling complexes and histone-modification enzymes) come together to orchestrate a finely tuned sequence of events that results in the complex pattern of CD4, CD8A and CD8B gene expression that is observed during thymocyte development.

Comparison of orthologous and paralogous DNA flanking the wheat high molecular weight glutenin genes: sequence conservation and divergence, transposon distribution, and matrix-attachment regions

Extended flanking DNA sequences were characterized for five members of the wheat high molecular weight (HMW) glutenin gene family to understand more of the structure, control, and evolution of these genes. Analysis revealed more sequence conservation among orthologous regions than between paralogous regions, with differences mainly owing to transposition events involving putative retrotransposons and several miniature inverted transposable elements (MITEs). Both gyspy-like long terminal repeat (LTR) and non-LTR retrotransposon sequences are represented in the flanking DNAs. One of the MITEs is a novel class, but another MITE is related to the maize Stowaway family and is widely represented in Triticeae express sequence tags (ESTs). Flanking DNA of the longest sequence, a 20 425-bp fragment including and surrounding the HMW-glutenin Bx7 gene, showed additional cereal gene-like sequences both immediately 5' and 3' to the HMW-glutenin coding region. The transcriptional activities of sequences related to these flanking putative genes and the retrotransposon-related regions were indicated by matches to wheat and other Triticeae ESTs. Predictive analysis of matrix-attachment regions (MARs) of the HMW glutenin and several alpha-, gamma-, and omega-gliadin flanking DNAs indicate potential MARs immediately flanking each of the genes. Matrix binding activity in the predicted regions was confirmed for two of the HMW-glutenin genes.

Breakpoints of t(4;11) translocations in the human MLL and AF4 genes in ALL patients are preferentially clustered outside of high-affinity matrix attachment regions

Chromosomal translocations t(4;11) are based on illegitimate recombinations between the human MLL and AF4 genes, and are associated with high-risk acute leukemias of infants and young children. Here, the question was asked, whether a correlation exists between the location of translocation breakpoints within both genes and the location of S/MARs. In "halo mapping experiments" (to define SARs), about 20 kb of MLL DNA was found to be attached to the nuclear matrix. Similar experiments performed for the translocation partner gene AF4 revealed that SARs are spanning nearly the complete breakpoint cluster region of the AF4 gene. By using short DNA fragments in "scaffold reassociation experiments" (to define MARs), similar results were obtained for both genes. However, Distamycin A competition experiments in combination with "scaffold reassociation experiments" revealed specific differences in the affinity of each tested DNA fragment to bind the isolated nuclear matrix proteins. When the latter data were compared with the known location of chromosomal breakpoints for both genes, an unexpected correlation was observed. DNA areas with strong MAR affinity contained fewer translocation breakpoints, while areas with weak or absent MAR affinity showed a higher density of chromosomal breakpoints.

AT-rich islands in genomic DNA as a novel target for AT-specific DNA-reactive antitumor drugs

Interstrand cross-links at T(A/T)4A sites in cellular DNA are associated with hypercytotoxicity of an anticancer drug, bizelesin. Here we evaluated whether these lethal effects reflect targeting critical genomic regions. An in silico analysis of human sequences showed that T(A/T)4A motifs are on average scarce and scattered. However, significantly higher local motif densities were identified in distinct minisatellite regions (200-1000 base pairs of approximately 85-100% AT), herein referred to as "AT islands." Experimentally detected bizelesin lesions agree with these in silico predictions. Actual bizelesin adducts clustered within the model AT island naked DNA, whereas motif-poor sequences were only sparsely adducted. In cancer cells, bizelesin produced high levels of lesions (approximately 4.7-7.1 lesions/kilobase pair/microM drug) in several prominent AT islands, compared with markedly lower lesion levels in several motif-poor loci and in bulk cellular DNA (approximately 0.8-1.3 and approximately 0.9 lesions/kilobase pair/microM drug, respectively). The identified AT islands exhibit sequence attributes of matrix attachment regions (MARs), domains that organize DNA loops on the nuclear matrix. The computed "MAR potential" and propensity for supercoiling-induced duplex destabilization (both predictive of strong MARs) correlate with the total number of bizelesin binding sites. Hence, MAR-like AT-rich non-coding domains can be regarded as a novel class of critical targets for anticancer drugs.

The transgeneticist's toolbox: novel methods for the targeted modification of eukaryotic genomes

Classical techniques for gene transfer into mammalian cells involve tedious screening procedures to identify transgenic clones or animals with the appropriate level and stability of expression or with the correct developmental patterns. These first generation technologies are clearly inadequate for complex genetic strategies by which gene regulation can be studied in its entire complexity. While site-specific insertions can principally be achieved by homologous recombination or by adapting the recombination apparatus from phages or yeast, these methods usually lack the required efficiency or they perturb expression patterns by the co-insertion of prokaryotic vector parts. Virtually all of these problems can be overcome by recombinase-mediated cassette exchange (RMCE) techniques which cleanly replace a resident cassette that is flanked by two hetero-specific recombination target sites for a second cassette with the analogous design, presented on a targeting vector. After illustrating the fundamentals of site-specific recombination by selected experiments, the authors (arranged in the chronological order of their contribution) will describe their efforts to develop RMCE into a method of wide applicability. Further developments that have been initiated utilizing the particular potential of the RMCE principle will be outlined.

Conserved characteristics of heterochromatin-forming DNA at the 15q11-q13 imprinting center

Nuclear matrix binding assays (NMBAs) define certain DNA sequences as matrix attachment regions (MARs), which often have cis-acting epigenetic regulatory functions. We used NMBAs to analyze the functionally important 15q11-q13 imprinting center (IC). We find that the IC is composed of an unusually high density of MARs, located in close proximity to the germ line elements that are proposed to direct imprint switching in this region. Moreover, we find that the organization of MARs is the same at the homologous mouse locus, despite extensive divergence of DNA sequence. MARs of this size are not usually associated with genes but rather with heterochromatin-forming areas of the genome. In contrast, the 15q11-q13 region contains multiple transcribed genes and is unusual for being subject to genomic imprinting, causing the maternal chromosome to be more transcriptionally silent, methylated, and late replicating than the paternal chromosome. We suggest that the extensive MAR sequences at the IC are organized as heterochromatin during oogenesis, an organization disrupted during spermatogenesis. Consistent with this model, multicolor fluorescence in situ hybridization to halo nuclei demonstrates a strong matrix association of the maternal IC, whereas the paternal IC is more decondensed, extending into the nuclear halo. This model also provides a mechanism for spreading of the imprinting signal, because heterochromatin at the IC on the maternal chromosome may exert a suppressive position effect in cis. We propose that the germ line elements at the 15q11-q13 IC mediate their effects through the candidate heterochromatin-forming DNA identified in this study.

Nuclear matrix attachment regions of human papillomavirus type 16 point toward conservation of these genomic elements in all genital papillomaviruses

The gene functions, transcriptional regulation, and genome replication of human papillomaviruses (HPVs) have been extensively studied. Thus far, however, there has been little research on the organization of HPV genomes in the nuclei of infected cells. As a first step to understand how chromatin and suprachromatin structures may modulate the life cycles of these viruses, we have identified and mapped interactions of HPV DNAs with the nuclear matrix. The endogenous genomes of HPV type 16 (HPV-16) which are present in SiHa, HPKI, and HPKII cells, adhere in vivo to the nuclear matrixes of these cell lines. A tight association with the nuclear matrix in vivo may be common to all genital HPV types, as the genomes of HPV-11, HPV-16, HPV-18, and HPV-33 showed high affinity in vitro to preparations of the nuclear matrix of C33A cells, as did the well-known nuclear matrix attachment region (MAR) of the cellular beta interferon gene. Affinity to the nuclear matrix is not evenly spread over the HPV-16 genome. Five genomic segments have strong MAR properties, while the other parts of the genome have low or no affinity. Some of the five MARs correlate with known cis-responsive elements: a strong MAR lies in the 5' segment of the long control region (LCR), and another one lies in the E6 gene, flanking the HPV enhancer, the replication origin, and the E6 promoter. The strongest MAR coincides with the E5 gene and the early-late intergenic region. Weak MAR activity is present in the E1 and E2 genes and in the 3' part of L2. The in vitro map of MAR activity appears to reflect MAR properties in vivo, as we found for two selected fragments with and without MAR activity. As is typical for many MARs, the two segments with highest affinity, namely, the 5' LCR and the early-late intergenic region, have an extraordinarily high A-T content (up to 85%). It is likely that these MARs have specific functions in the viral life cycle, as MARs predicted by nucleotide sequence analysis, patterns of A-T content, transcription factor YY1 binding sites, and likely topoisomerase II cleavage sites are conserved in similar positions throughout all genital HPVs.

Stress-induced duplex DNA destabilization in scaffold/matrix attachment regions

S/MARs are DNA elements 300 to several thousand base-pairs long, which are operationally defined by their affinity for the nuclear scaffold or matrix. S/MARs occur exclusively in eukaryotic genomes, where they mediate several functions. Because S/MARs do not have a clearcut consensus sequence, the characteristics that define their activity are thought to be structural. Ubiquitous S/MAR binding proteins have been identified, but to date no unique binding sequence or structural motif has been found. Here we show by computational analysis that S/MARs conform to a specific design whose essential attribute is the presence of stress-induced base-unpairing regions (BURs). Stress-induced destabilization (SIDD) profiles are calculated using a previously developed statistical mechanical procedure in which the superhelical deformation is partitioned between strand separation, twisting within denatured regions, and residual superhelicity. The results of these calculations show that BURs exhibit a succession of evenly spaced destabilized sites that would render part or all of the S/MAR sequence single stranded at sufficient superhelicity. These analyses are performed for a range of sequenced S/MAR elements from the borders of eukaryotic gene domains, from centromeres, and from positions where S/MARs are known to support the action of an enhancer. The results reported here are in excellent agreement with earlier in vitro chemical reactivity studies. This approach demonstrates the potential for computational analysis to predict the points of division of the eukaryotic genome into functional units (domains), and also to locate certain cis-regulatory sequences.

A novel method to develop highly specific models for regulatory units detects a new LTR in GenBank which contains a functional promoter

Functional promoters are composed of individual modules (e.g. transcription factor binding sites, secondary structure elements, repeats) arranged in distinct patterns. Recognition of such patterns is essential for identification of promoters in non-coding sequences. However, this is difficult due to the absence of overall sequence similarity in promoters even if they are regulated in a similar way. We implemented simple formal representations of general features of regulatory regions into an algorithm capable of developing complex models reflecting both the element composition and the functional organization of individual elements (ModelGenerator). Though ModelGenerator requires a very simple initial model (e.g. two modules and their relative order) it will generate a much more sophisticated model by analysis of the training set of at least ten sequences. We show ModelGenerator to successfully model different retroviral long terminal repeat (LTR) classes (Lentivirus as well as avian and mammalian C-type) which contain functional promoters. Database searches with the program ModelInspector demonstrated the high specificity of these models and no apparent false negatives were detected. We also verified one match from GenBank to the mammalian C-type LTR model experimentally and showed this sequence to contain an active promoter. Thus, the concept of modular organization of functional regulatory DNA regions (e.g. promoters) could be successfully implemented into a set of computer tools which might be flexible and specific enough to be suitable for prospective analysis of new genomic DNA sequences.

Mathematical model to predict regions of chromatin attachment to the nuclear matrix

The potentiation and subsequent initiation of transcription are complex biological phenomena. The region of attachment of the chromatin fiber to the nuclear matrix, known as the matrix attachment region or scaffold attachment region (MAR or SAR), are thought to be requisite for the transcriptional regulation of the eukaryotic genome. As expressed sequences should be contained in these regions, it becomes significant to answer the following question: can these regions be identified from the primary sequence data alone and subsequently used as markers for expressed sequences? This paper represents an effort toward achieving this goal and describes a mathematical model for the detection of MARs. The location of matrix associated regions has been linked to a variety of sequence patterns. Consequently, a list of these patterns is compiled and represented as a set of decision rules using an AND-OR formulation. The DNA sequence was then searched for the presence of these patterns and a statistical significance was associated with the frequency of occurrence of the various patterns. Subsequently, a mathematical potential value,MAR-Potential, was assigned to a sequence region as the inverse proportion to the probability that the observed pattern population occurred at random. Such a MAR detection process was applied to the analysis of a variety of known MAR containing sequences. Regions of matrix association predicted by the software essentially correspond to those determined experimentally. The human T-cell receptor and the DNA sequence from the Drosophila bithorax region were also analyzed. This demonstrates the usefulness of the approach described as a means to direct experimental resources.

Matrix-attachment regions in the mouse chromosome 7F imprinted domain

We have mapped the matrix-attachment regions (MARs) in 200 kilobases of the mouse Chromosome (Chr) 7F imprinted domain. MARs are genetic elements known to have effects in cis on methylation at nonimprinted loci. The imprinting of the Igf2 and Ins2 genes is dependent on the transcription of the downstream H19 gene. The transcription of H19 is dependent in turn on its methylation status. The cis-acting regulators of methylation at this site are not known. As MARs are potential regulators not only of methylation but also other elements of genomic imprinting, we mapped the MARs within the 200 kilobases around H19. This report describes the mapping of four MARs from this region.