Comparative study and prediction of DNA fragments associated with various elements of the nuclear matrix

Computational prediction of polycomb-associated long non-coding RNAs

Among thousands of long non-coding RNAs (lncRNAs) only a small subset is functionally characterized and the functional annotation of lncRNAs on the genomic scale remains inadequate. In this study we computationally characterized two functionally different parts of human lncRNAs transcriptome based on their ability to bind the polycomb repressive complex, PRC2. This classification is enabled by the fact that while all lncRNAs constitute a diverse set of sequences, the classes of PRC2-binding and PRC2 non-binding lncRNAs possess characteristic combinations of sequence-structure patterns and, therefore, can be separated within the feature space. Based on the specific combination of features, we built several machine-learning classifiers and identified the SVM-based classifier as the best performing. We further showed that the SVM-based classifier is able to generalize on the independent data sets. We observed that this classifier, trained on the human lncRNAs, can predict up to 59.4% of PRC2-binding lncRNAs in mice. This suggests that, despite the low degree of sequence conservation, many lncRNAs play functionally conserved biological roles.

Identification of Scaffold/Matrix Attachment (S/MAR) like DNA element from the gastrointestinal protozoan parasite Giardia lamblia

BACKGROUND

Chromatin in the nucleus of all eukaryotes is organized into a system of loops and domains. These loops remain fastened at their bases to the fundamental framework of the nucleus, the matrix or the scaffold. The DNA sequences which anchor the bases of the chromatin loops to the matrix are known as Scaffold/Matrix Attachment Regions or S/MARs. Though S/MARs have been studied in yeast and higher eukaryotes and they have been found to be associated with gene organization and regulation of gene expression, they have not been reported in protists like Giardia. Several tools have been discovered and formulated to predict S/MARs from a genome of a higher eukaryote which take into account a number of features. However, the lack of a definitive consensus sequence in S/MARs and the randomness of the protozoan genome in general, make it a challenge to predict and identify such sequences from protists.

RESULTS

Here, we have analysed the Giardia genome for the probable S/MARs predicted by the available computational tools; and then shown these sequences to be physically associated with the nuclear matrix. Our study also reflects that while no single computational tool is competent to predict such complex elements from protist genomes, a combination of tools followed by experimental verification is the only way to confirm the presence of these elements from these organisms.

CONCLUSION

This is the first report of S/MAR elements from the protozoan parasite Giardia lamblia. This initial work is expected to lay a framework for future studies relating to genome organization as well as gene regulatory elements in this parasite.

Determination of the in vivo structural DNA loop organization in the genomic region of the rat albumin locus by means of a topological approach

Nuclear DNA of metazoans is organized in supercoiled loops anchored to a proteinaceous substructure known as the nuclear matrix (NM). DNA is anchored to the NM by non-coding sequences known as matrix attachment regions (MARs). There are no consensus sequences for identification of MARs and not all potential MARs are actually bound to the NM constituting loop attachment regions (LARs). Fundamental processes of nuclear physiology occur at macromolecular complexes organized on the NM; thus, the topological organization of DNA loops must be important. Here, we describe a general method for determining the structural DNA loop organization in any large genomic region with a known sequence. The method exploits the topological properties of loop DNA attached to the NM and elementary topological principles such as that points in a deformable string (DNA) can be positionally mapped relative to a position-reference invariant (NM), and from such mapping, the configuration of the string in third dimension can be deduced. Therefore, it is possible to determine the specific DNA loop configuration without previous characterization of the LARs involved. We determined in hepatocytes and B-lymphocytes of the rat the DNA loop organization of a genomic region that contains four members of the albumin gene family.

Association of the matrix attachment region recognition signature with coding regions in Caenorhabditis elegans

Background

Matrix attachment regions (MAR) are the sites on genomic DNA that interact with the nuclear matrix. There is increasing evidence for the involvement of MAR in regulation of gene expression. The unsuitability of experimental detection of MAR for genome-wide analyses has led to the development of computational methods of detecting MAR. The MAR recognition signature (MRS) has been reported to be associated with a significant fraction of MAR in C. elegans and has also been found in MAR from a wide range of other eukaryotes. However the effectiveness of the MRS in specifically and sensitively identifying MAR remains unresolved.

Results

Using custom software, we have mapped the occurrence of MRS across the entire C. elegans genome. We find that MRS have a distinctive chromosomal distribution, in which they appear more frequently in the gene-rich chromosome centres than in arms. Comparison to distributions of MRS estimated from chromosomal sequences randomised using mono-, di- tri- and tetra-nucleotide frequency patterns showed that, while MRS are less common in real sequence than would be expected from nucleotide content alone, they are more frequent than would be predicted from short-range nucleotide structure. In comparison to the rest of the genome, MRS frequency was elevated in 5' and 3' UTRs, and striking peaks of average MRS frequency flanked C. elegans coding sequence (CDS). Genes associated with MRS were significantly enriched for receptor activity annotations, but not for expression level or other features.

Conclusion

Through a genome-wide analysis of the distribution of MRS in C. elegans we have shown that they have a distinctive distribution, particularly in relation to genes. Due to their association with untranslated regions, it is possible that MRS could have a post-transcriptional role in the control of gene expression. A role for MRS in nuclear scaffold attachment is not supported by these analyses.

A comparative study of S/MAR prediction tools

BACKGROUND

S/MARs are regions of the DNA that are attached to the nuclear matrix. These regions are known to affect substantially the expression of genes. The computer prediction of S/MARs is a highly significant task which could contribute to our understanding of chromatin organisation in eukaryotic cells, the number and distribution of boundary elements, and the understanding of gene regulation in eukaryotic cells. However, while a number of S/MAR predictors have been proposed, their accuracy has so far not come under scrutiny.

RESULTS

We have selected S/MARs with sufficient experimental evidence and used these to evaluate existing methods of S/MAR prediction. Our main results are: 1.) all existing methods have little predictive power, 2.) a simple rule based on AT-percentage is generally competitive with other methods, 3.) in practice, the different methods will usually identify different sub-sequences as S/MARs, 4.) more research on the H-Rule would be valuable.

CONCLUSION

A new insight is needed to design a method which will predict S/MARs well. Our data, including the control data, has been deposited as additional material and this may help later researchers test new predictors.

Scaffold/Matrix Attachment Regions and intrinsic DNA curvature

Recent approaches have failed to detect nucleotide sequence motifs in Scaffold/Matrix Attachment Regions (S/MARs). The lack of any known motifs, together with the confirmation that some S/MARs are not associated to any peculiar sequence, indicates that some structural elements, such as DNA curvature, have a role in chromatin organization and on their efficiency in protein binding. Similar to DNA curvature, S/MARs are located close to promoters, replication origins, and multiple nuclear processes like recombination and breakpoint sites. The chromatin structure in these regulatory regions is important to chromosome organization for accurate regulation of nuclear processes. In this article we review the biological importance of the co-localization between bent DNA sites and S/MARs.

In-silico prediction and observations of nuclear matrix attachment

The nuclear matrix is a functionally adaptive structural framework interior to the nuclear envelope. The nature and function of this nuclear organizer remains the subject of widespread discussion in the epigenetic literature. To draw this discussion together with a view to suggest a way forward we summarize the biochemical evidence for the modalities of DNA-matrix binding alongside the in-silico predictions. Concordance is exhibited at various, but not all levels. On the one hand, both the reiteration and sequence similarity of some elements of Matrix Attachment Regions suggest conservation. On the other hand, in-silico predictions suggest additional unique components. In bringing together biological and sequence evidence we conclude that binding may be hierarchical in nature, reflective of a biological role in replicating, transcribing and potentiating chromatin. Nuclear matrix binding may well be more complex than the widely accepted simple loop model.

Matrix attachment regions as targets for retroviral integration

BACKGROUND

The randomness of retroviral integration has been debated for many years. Recent evidence indicates that integration site selection is not random, and that it is influenced by both viral and cellular factors. To study the role of DNA structure in site selection, retroviral integration near matrix attachment regions (MARs) was analyzed for three different groups of retroviruses. The objective was to assess whether integration near MARs may be a factor for integration site selection.

RESULTS

Results indicated that MLV, SL3-3 MuLV, HIV-1 and HTLV-1 integrate preferentially near MARs, specifically within 2-kilobases (kb). In addition, a preferential position and orientation relative to the adjacent MAR was observed for each virus. Further analysis of SL3-3 MuLV insertions in common integration sites (CISs) demonstrated a higher frequency of integration near MARs and an orientation preference that was not observed for integrations outside CISs.

CONCLUSION

These findings contribute to a growing body of evidence indicating that retroviral integration is not random, that MARs influence integration site selection for some retroviruses, and that integration near MARs may have a role in the insertional activation of oncogenes by gammaretroviruses.

Comparative analysis of chromatin landscape in regulatory regions of human housekeeping and tissue specific genes

BACKGROUND

Global regulatory mechanisms involving chromatin assembly and remodelling in the promoter regions of genes is implicated in eukaryotic transcription control especially for genes subjected to spatial and temporal regulation. The potential to utilise global regulatory mechanisms for controlling gene expression might depend upon the architecture of the chromatin in and around the gene. In-silico analysis can yield important insights into this aspect, facilitating comparison of two or more classes of genes comprising of a large number of genes within each group.

RESULTS

In the present study, we carried out a comparative analysis of chromatin characteristics in terms of the scaffold/matrix attachment regions, nucleosome formation potential and the occurrence of repetitive sequences, in the upstream regulatory regions of housekeeping and tissue specific genes. Our data show that putative scaffold/matrix attachment regions are more abundant and nucleosome formation potential is higher in the 5' regions of tissue specific genes as compared to the housekeeping genes.

CONCLUSION

The differences in the chromatin features between the two groups of genes indicate the involvement of chromatin organisation in the control of gene expression. The presence of global regulatory mechanisms mediated through chromatin organisation can decrease the burden of invoking gene specific regulators for maintenance of the active/silenced state of gene expression. This could partially explain the lower number of genes estimated in the human genome.

Human matrix attachment regions are necessary for the establishment but not the maintenance of transgene insulation in Drosophila melanogaster

Human matrix attachment regions (MARs) can insulate transgene expression from chromosomal position effects in Drosophila melanogaster. To gain insight into the mechanism(s) by which chromosomal insulation occurs, we studied the expression phenotypes of Drosophila transformants expressing mini-white transgenes in which MAR sequences from the human apoB gene were arranged in a variety of ways. In agreement with previous reports, we found that a single copy of the insulating element was not sufficient for position-independent transgene expression; rather, two copies were required. However, the arrangement of the two elements within the transgene was unimportant, since chromosomal insulation was equally apparent when both copies of the insulator were upstream of the mini-white reporter as when the transcription unit was flanked by insulator elements. Moreover, experiments in which apoB 3' MAR sequences were removed from integrated transgenes in vivo by site-specific recombination demonstrated that MAR sequences were required for the establishment but not for the maintenance of chromosomal insulation. These observations are not compatible with the chromosomal loop model in its simplest form. Alternate mechanisms for MAR function in this system are proposed.

Comparative analysis of the PDCD2–TBP–PSMB1 region in vertebrates

Three orthologous genes encoding programmed cell death 2 (PDCD2), TATA-binding protein (TBP), and proteasomal subunit C5 (PSMB1) proteins have been shown previously to be nonrandomly distributed in both mammalian and invertebrate genomes. Here we analyze a conserved synteny of the PDCD2, TBP, and PSMB1 orthologs in four nonmammalian vertebrates. Homologous genes of the chicken, zebrafish, fugu, and Tetraodon nigroviridis were identified. A chicken cosmid harboring the orthologs of these three genes was completely sequenced. The fish genes were analyzed in silico. In all seven vertebrates thus far investigated, the PDCD2 and TBP genes are located tail-to-tail. In six tested species but the zebrafish, the PSMB1 gene mapped head-to-head or in the close vicinity to the TBP, but even in the zebrafish, all three genes were syntenic. In contrast, a three times reused synteny breakpoint in the 5'-region from PDCD2 was detected. A comparative analysis revealed the distribution of putative matrix-attached regions (MARs), which may affect the synteny conservation.

A short synthetic MAR positively affects transgene expression in rice and Arabidopsis

Matrix Attachment Regions (MARs) are DNA elements that are thought to influence gene expression by anchoring active chromatin domains to the nuclear matrix. When flanking a construct in transgenic plants, MARs could be useful for enhancing transgene expression. Naturally occurring MARs have a number of sequence features and DNA elements in common, and using different subsets of these sequence elements, three independent synthetic MARs were created. Although short, these MARs were able to bind nuclear scaffold preparations with an affinity equal to or greater than naturally occurring plant MARs. One synthetic MAR was extensively tested for its effect on transgene expression, using different MAR orientations, plant promoters, transformation methods and plant species. This MAR was able to increase average transgene expression and produced integration patterns of lower complexity. These data show the potential of making well defined synthetic MARs and using them to improve transgene expression.

A rapid method of genomic array analysis of scaffold/matrix attachment regions (S/MARs) identifies a 2.5-Mb region of enhanced scaffold/matrix attachment at a human neocentromere

Human neocentromeres are fully functional centromeres that arise at previously noncentromeric regions of the genome. We have tested a rapid procedure of genomic array analysis of chromosome scaffold/matrix attachment regions (S/MARs), involving the isolation of S/MAR DNA and hybridization of this DNA to a genomic BAC/PAC array. Using this procedure, we have defined a 2.5-Mb domain of S/MAR-enriched chromatin that fully encompasses a previously mapped centromere protein-A (CENP-A)-associated domain at a human neocentromere. We have independently verified this procedure using a previously established fluorescence in situ hybridization method on salt-treated metaphase chromosomes. In silico sequence analysis of the S/MAR-enriched and surrounding regions has revealed no outstanding sequence-related predisposition. This study defines the S/MAR-enriched domain of a higher eukaryotic centromere and provides a method that has broad application for the mapping of S/MAR attachment sites over large genomic regions or throughout a genome.

A significant fraction of conserved noncoding DNA in human and mouse consists of predicted matrix attachment regions

Noncoding DNA in the human-mouse orthologous intergenic regions contains "islands" of conserved sequences, the functions of which remain largely unknown. We hypothesized that some of these regions might be matrix-scaffold attachment regions, MARs (or S/MARs). MARs comprise one of the few classes of eukaryotic noncoding DNA with an experimentally characterized function, being involved in the attachment of chromatin to the nuclear matrix, chromatin remodeling and transcription regulation. To test our hypothesis, we analyzed the co-occurrence of predicted MARs with highly conserved noncoding DNA regions in human-mouse genomic alignments. We found that 11% of the conserved noncoding DNA consists of predicted MARs. Conversely, more than half of the predicted MARs co-occur with one or more independently identified conserved sequence blocks. An excess of conserved predicted MARs is seen in intergenic regions preceding 5' ends of genes, suggesting that these MARs are primarily involved in transcriptional control.

High frequency of matrix attachment regions and cut-like protein x/CCAAT-displacement protein and B cell regulator of IgH transcription binding sites flanking Ig V region genes

A major component in controlling V(D)J recombination is differential accessibility through localized changes in chromatin structure. Attachment of DNA to the nuclear matrix via matrix attachment region (MAR) sequences, and interaction with MAR-binding proteins have been shown to alter chromatin conformation, promote histone acetylation, and influence gene transcription. In this study, the flanking regions of several human and mouse Ig V(H) and Ig Vkappa genes were analyzed extensively for the presence of MARs by in vitro matrix-binding assay, and for interaction with the MAR-binding proteins cut-like protein x/CCAAT-displacement protein (Cux/CDP), B cell regulator of IgH transcription (Bright), and special AT-rich sequence-binding protein (SATB1) by EMSA. Cux/CDP and SATB1 are associated with repression, while Bright is an activator of Ig transcription. Binding sites were identified in the vicinity of all analyzed Ig V genes, and were also found flanking TCR Vbeta genes. We also show that the binding sites of the different factors do not always occur at MAR sequences. MAR sequences were also found within the Ig V loci at a much higher frequency than throughout the rest of the genome. Overall, the frequency and location of binding sites relative to the coding regions, and the strength of DNA-protein interaction showed much heterogeneity. Thus, variations in factor binding and MAR activity could potentially influence the extent of localized accessibility to V(D)J recombination and thus could play a role in unequal rearrangement of individual V genes. These sites could also contribute to effective transcription of Ig genes in mature and/or activated B cells, bringing both the promoter as well as the enhancer regions into close proximity at the nuclear matrix.

Cytoplasmic intermediate filaments are stably associated with nuclear matrices and potentially modulate their DNA-binding function

The tight association of cytoplasmic intermediate filaments (cIFs) with the nucleus and the isolation of crosslinkage products of vimentin with genomic DNA fragments, including nuclear matrix attachment regions (MARs) from proliferating fibroblasts, point to a participation of cIFs in nuclear activities. To test the possibility that cIFs are complementary nuclear matrix elements, the nuclei of a series of cultured cells were subjected to the Li-diiodosalicylate (LIS) extraction protocol developed for the preparation of nuclear matrices and analyzed by immunofluorescence microscopy and immunoblotting with antibodies directed against lamin B and cIF proteins. When nuclei released from hypotonically swollen L929 suspension cells in the presence of digitonin or Triton X-100 were exposed to such strong shearing forces that a considerable number were totally disrupted, a thin, discontinuous layer of vimentin IFs remained tenaciously adhering to still intact nuclei, in apparent coalignment with the nuclear lamina. Even in broken nuclei, the distribution of vimentin followed that of lamin B in areas where the lamina still appeared intact. The same retention of vimentin together with desmin and glial IFs was observed on the nuclei isolated from differentiating C2C12 myoblast and U333 glioma cells, respectively. Nuclei from epithelial cells shed their residual perinuclear IF layers as coherent cytoskeletal ghosts, except for small fractions of vimentin and cytokeratin IFs, which remained in a dot-to cap-like arrangement on the nuclear surface, in apparent codistribution with lamin B. LIS extraction did not bring about a reduction in the cIF protein contents of such nuclei upon their transformation into nuclear matrices. Moreover, in whole mount preparations of mouse embryo fibroblasts, DNA/chromatin emerging from nuclei during LIS extraction mechanically and chemically cleaned the nuclear surface and perinuclear area from loosely anchored cytoplasmic material with the production of broad, IF-free annular spaces, but left substantial fractions of the vimentin IFs in tight association with the nuclear surface. Accordingly, double-immunogold electron microscopy of fixed and permeabilized fibroblasts disclosed a close neighborhood of vimentin IFs and lamin B, with a minimal distance between the nanogold particles of ca. 30 nm. These data indicate an extremely solid interconnection of cIFs with structural elements of the nuclear matrix, and make them, together with their susceptibility to crosslinkage to MARs and other genomic DNA sequences under native conditions, complementary or even integral constituents of the karyoskeleton.

Small tumor virus genomes are integrated near nuclear matrix attachment regions in transformed cells

More than 15% of human cancers have a viral etiology. In benign lesions induced by the small DNA tumor viruses, viral genomes are typically maintained extrachromosomally. Malignant progression is often associated with viral integration into host cell chromatin. To study the role of viral integration in tumorigenesis, we analyzed the positions of integrated viral genomes in tumors and tumor cell lines induced by the small oncogenic viruses, including the high-risk human papillomaviruses, hepatitis B virus, simian virus 40, and human T-cell leukemia virus type 1. We show that viral integrations in tumor cells lie near cellular sequences identified as nuclear matrix attachment regions (MARs), while integrations in nonneoplastic cells show no significant correlation with these regions. In mammalian cells, the nuclear matrix functions in gene expression and DNA replication. MARs play varied but poorly understood roles in eukaryotic gene expression. Our results suggest that integrated tumor virus genomes are subject to MAR-mediated transcriptional regulation, providing insight into mechanisms of viral carcinogenesis. Furthermore, the viral oncoproteins serve as invaluable tools for the study of mechanisms controlling cellular growth. Similarly, our demonstration that integrated viral genomes may be subject to MAR-mediated transcriptional effects should facilitate elucidation of fundamental mechanisms regulating eukaryotic gene expression.