In silico prediction of scaffold/matrix attachment regions in large genomic sequences

Mapping of scaffold/matrix attachment regions in human genome: a data mining exercise

Scaffold/matrix attachment regions (S/MARs) are DNA elements that serve to compartmentalize the chromatin into structural and functional domains. These elements are involved in control of gene expression which governs the phenotype and also plays role in disease biology. Therefore, genome-wide understanding of these elements holds great therapeutic promise. Several attempts have been made toward identification of S/MARs in genomes of various organisms including human. However, a comprehensive genome-wide map of human S/MARs is yet not available. Toward this objective, ChIP-Seq data of 14 S/MAR binding proteins were analyzed and the binding site coordinates of these proteins were used to prepare a non-redundant S/MAR dataset of human genome. Along with co-ordinate (location) details of S/MARs, the dataset also revealed details of S/MAR features, namely, length, inter-SMAR length (the chromatin loop size), nucleotide repeats, motif abundance, chromosomal distribution and genomic context. S/MARs identified in present study and their subsequent analysis also suggests that these elements act as hotspots for integration of retroviruses. Therefore, these data will help toward better understanding of genome functioning and designing effective anti-viral therapeutics. In order to facilitate user friendly browsing and retrieval of the data obtained in present study, a web interface, MARome (http://bioinfo.net.in/MARome), has been developed.

Chromosome territories and the global regulation of the genome

Spatial positioning is a fundamental principle governing nuclear processes. Chromatin is organized as a hierarchy from nucleosomes to Mbp chromatin domains (CD) or topologically associating domains (TADs) to higher level compartments culminating in chromosome territories (CT). Microscopic and sequencing techniques have substantiated chromatin organization as a critical factor regulating gene expression. For example, enhancers loop back to interact with their target genes almost exclusively within TADs, distally located coregulated genes reposition into common transcription factories upon activation, and Mbp CDs exhibit dynamic motion and configurational changes in vivo. A longstanding question in the nucleus field is whether an interactive nuclear matrix provides a direct link between structure and function. The findings of nonrandom radial positioning of CT within the nucleus suggest the possibility of preferential interaction patterns among populations of CT. Sequential labeling up to 10 CT followed by application of computer imaging and geometric graph mining algorithms revealed cell-type specific interchromosomal networks (ICN) of CT that are altered during the cell cycle, differentiation, and cancer progression. It is proposed that the ICN correlate with the global level of genome regulation. These approaches also demonstrated that the large scale 3-D topology of CT is specific for each CT. The cell-type specific proximity of certain chromosomal regions in normal cells may explain the propensity of distinct translocations in cancer subtypes. Understanding how genes are dysregulated upon disruption of the normal "wiring" of the nucleus by translocations, deletions, and amplifications that are hallmarks of cancer, should enable more targeted therapeutic strategies.

Matrix association region/scaffold attachment region: the crucial player in defining the positions of chromosome breaks mediated by bile acid-induced apoptosis in nasopharyngeal epithelial cells

Background

It has been found that chronic rhinosinusitis (CRS) increases the risk of developing nasopharyngeal carcinoma (NPC). CRS can be caused by gastro-oesophageal reflux (GOR) that may reach nasopharynx. The major component of refluxate, bile acid (BA) has been found to be carcinogenic and genotoxic. BA-induced apoptosis has been associated with various cancers. We have previously demonstrated that BA induced apoptosis and gene cleavages in nasopharyngeal epithelial cells. Chromosomal cleavage occurs at the early stage of both apoptosis and chromosome rearrangement. It was suggested that chromosome breaks tend to cluster in the region containing matrix association region/scaffold attachment region (MAR/SAR). This study hypothesised that BA may cause chromosome breaks at MAR/SAR leading to chromosome aberrations in NPC. This study targeted the AF9 gene located at 9p22 because 9p22 is a deletion hotspot in NPC.

Methods

Potential MAR/SAR sites were predicted in the AF9 gene by using MAR/SAR prediction tools. Normal nasopharyngeal epithelial cells (NP69) and NPC cells (TWO4) were treated with BA at neutral and acidic pH. Inverse-PCR (IPCR) was used to identify chromosome breaks in SAR region (contains MAR/SAR) and non-SAR region (does not contain MAR/SAR). To map the chromosomal breakpoints within the AF9 SAR and non-SAR regions, DNA sequencing was performed.

Results

In the AF9 SAR region, the gene cleavage frequencies of BA-treated NP69 and TWO4 cells were significantly higher than those of untreated control. As for the AF9 non-SAR region, no significant difference in cleavage frequency was detected between untreated and BA-treated cells. A few breakpoints detected in the SAR region were mapped within the AF9 region that was previously reported to translocate with the mixed lineage leukaemia (MLL) gene in an acute lymphoblastic leukaemia (ALL) patient.

Conclusions

Our findings suggest that MAR/SAR may be involved in defining the positions of chromosomal breakages induced by BA. Our report here, for the first time, unravelled the relation of these BA-induced chromosomal breakages to the AF9 chromatin structure.

Electronic supplementary material

The online version of this article (10.1186/s12920-018-0465-4) contains supplementary material, which is available to authorized users.

Genome-Wide Anaplasma phagocytophilum AnkA-DNA Interactions Are Enriched in Intergenic Regions and Gene Promoters and Correlate with Infection-Induced Differential Gene Expression

Anaplasma phagocytophilum, an obligate intracellular prokaryote, infects neutrophils, and alters cardinal functions via reprogrammed transcription. Large contiguous regions of neutrophil chromosomes are differentially expressed during infection. Secreted A. phagocytophilum effector AnkA transits into the neutrophil or granulocyte nucleus to complex with DNA in heterochromatin across all chromosomes. AnkA binds to gene promoters to dampen cis-transcription and also has features of matrix attachment region (MAR)-binding proteins that regulate three-dimensional chromatin architecture and coordinate transcriptional programs encoded in topologically-associated chromatin domains. We hypothesize that identification of additional AnkA binding sites will better delineate how A. phagocytophilum infection results in reprogramming of the neutrophil genome. Using AnkA-binding ChIP-seq, we showed that AnkA binds broadly throughout all chromosomes in a reproducible pattern, especially at: (i) intergenic regions predicted to be MARs; (ii) within predicted lamina-associated domains; and (iii) at promoters ≤ 3000 bp upstream of transcriptional start sites. These findings provide genome-wide support for AnkA as a regulator of cis-gene transcription. Moreover, the dominant mark of AnkA in distal intergenic regions known to be AT-enriched, coupled with frequent enrichment in the nuclear lamina, provides strong support for its role as a MAR-binding protein and genome “re-organizer.” AnkA must be considered a prime candidate to promote neutrophil reprogramming and subsequent functional changes that belie improved microbial fitness and pathogenicity.

Genome-wide mapping of matrix attachment regions in Drosophila melanogaster

Background

Eukaryotic genome acquires functionality upon proper packaging within the nucleus. This process is facilitated by the structural framework of Nuclear Matrix, a nucleo-proteinaceous meshwork. Matrix Attachment Regions (MARs) in the genome serve as anchoring sites to this framework.

Results

Here we report direct sequencing of the MAR preparation from Drosophila melanogaster embryos and identify >7350 MARs. This amounts to ~2.5% of the fly genome and often coincide with AT rich non-coding regions. We find significant association of MARs with the origins of replication, transcription start sites, paused RNA Polymerase II sites and exons, but not introns, of highly expressed genes. We also identified sequence motifs and repeats that constitute MARs.

Conclusion

Our data reveal the contact points of genome to the nuclear architecture and provide a link between nuclear functions and genomic packaging.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1022) contains supplementary material, which is available to authorized users.

Dynamic epigenetic control of highly conserved noncoding elements

Background

Many noncoding genomic loci have remained constant over long evolutionary periods, suggesting that they are exposed to strong selective pressures. The molecular functions of these elements have been partially elucidated, but the fundamental reason for their extreme conservation is still unknown.

Results

To gain new insights into the extreme selection of highly conserved noncoding elements (HCNEs), we used a systematic analysis of multi-omic data to study the epigenetic regulation of such elements during the development of Drosophila melanogaster. At the sequence level, HCNEs are GC-rich and have a characteristic oligomeric composition. They have higher levels of stable nucleosome occupancy than their flanking regions, and lower levels of mononucleosomes and H3.3, suggesting that these regions reside in compact chromatin. Furthermore, these regions showed remarkable modulations in histone modification and the expression levels of adjacent genes during development. Although HCNEs are primarily initiated late in replication, about 10% were related to early replication origins. Finally, HCNEs showed strong enrichment within lamina-associated domains.

Conclusion

HCNEs have distinct and protective sequence properties, undergo dynamic epigenetic regulation, and appear to be associated with the structural components of the chromatin, replication origins, and nuclear matrix. These observations indicate that such elements are likely to have essential cellular functions, and offer insights into their epigenetic properties.

Genome analysis of a major urban malaria vector mosquito, Anopheles stephensi

Background

Anopheles stephensi is the key vector of malaria throughout the Indian subcontinent and Middle East and an emerging model for molecular and genetic studies of mosquito-parasite interactions. The type form of the species is responsible for the majority of urban malaria transmission across its range.

Results

Here, we report the genome sequence and annotation of the Indian strain of the type form of An. stephensi. The 221 Mb genome assembly represents more than 92% of the entire genome and was produced using a combination of 454, Illumina, and PacBio sequencing. Physical mapping assigned 62% of the genome onto chromosomes, enabling chromosome-based analysis. Comparisons between An. stephensi and An. gambiae reveal that the rate of gene order reshuffling on the X chromosome was three times higher than that on the autosomes. An. stephensi has more heterochromatin in pericentric regions but less repetitive DNA in chromosome arms than An. gambiae. We also identify a number of Y-chromosome contigs and BACs. Interspersed repeats constitute 7.1% of the assembled genome while LTR retrotransposons alone comprise more than 49% of the Y contigs. RNA-seq analyses provide new insights into mosquito innate immunity, development, and sexual dimorphism.

Conclusions

The genome analysis described in this manuscript provides a resource and platform for fundamental and translational research into a major urban malaria vector. Chromosome-based investigations provide unique perspectives on Anopheles chromosome evolution. RNA-seq analysis and studies of immunity genes offer new insights into mosquito biology and mosquito-parasite interactions.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0459-2) contains supplementary material, which is available to authorized users.

BAC-based cellular model for screening regulators of BDNF gene transcription

Background

Brain derived neurotrophic factor (BDNF) belongs to a family of structurally related proteins called neurotrophins that have been shown to regulate survival and growth of neurons in the developing central and peripheral nervous system and also to take part in synaptic plasticity related processes in adulthood. Since BDNF is associated with several nervous system disorders it would be beneficial to have cellular reporter system for studying its expression regulation.

Methods

Using modified bacterial artificial chromosome (BAC), we generated several transgenic cell lines expressing humanised Renilla luciferase (hRluc)-EGFP fusion reporter gene under the control of rat BDNF gene regulatory sequences (rBDNF-hRluc-EGFP) in HeLa background. To see if the hRluc-EGFP reporter was regulated in response to known regulators of BDNF expression we treated cell lines with substances known to regulate BDNF and also overexpressed transcription factors known to regulate BDNF gene in established cell lines.

Results

rBDNF-hRluc-EGFP cell lines had high transgene copy numbers when assayed with qPCR and FISH analysis showed that transgene was maintained episomally in all cell lines. Luciferase activity in transgenic cell lines was induced in response to ionomycin-mediated rise of intracellular calcium levels, treatment with HDAC inhibitors and by over-expression of transcription factors known to increase BDNF expression, indicating that transcription of the transgenic reporter is regulated similarly to the endogenous BDNF gene.

Conclusions

Generated rBDNF-hRluc-EGFP BAC cell lines respond to known modulators of BDNF expression and could be used for screening of compounds/small molecules or transcription factors altering BDNF expression.

In vivo mapping of arabidopsis scaffold/matrix attachment regions reveals link to nucleosome-disfavoring poly(dA:dT) tracts

Scaffold or matrix attachment regions (S/MARs) are found in all eukaryotes. The pattern of distribution and genomic context of S/MARs is thought to be important for processes such as chromatin organization and modulation of gene expression. Despite the importance of such processes, much is unknown about the large-scale distribution and sequence content of S/MARs in vivo. Here, we report the use of tiling microarrays to map 1358 S/MARs on Arabidopsis thaliana chromosome 4 (chr4). S/MARs occur throughout chr4, spaced much more closely than in the large plant and animal genomes that have been studied to date. Arabidopsis S/MARs can be divided into five clusters based on their association with other genomic features, suggesting a diversity of functions. While some Arabidopsis S/MARs may define structural domains, most occur near the transcription start sites of genes. Genes associated with these S/MARs have an increased probability of expression, which is particularly pronounced in the case of transcription factor genes. Analysis of sequence motifs and 6-mer enrichment patterns show that S/MARs are preferentially enriched in poly(dA:dT) tracts, sequences that resist nucleosome formation, and the majority of S/MARs contain at least one nucleosome-depleted region. This global view of S/MARs provides a framework to begin evaluating genome-scale models for S/MAR function.

HMGA regulates the global chromatin state and neurogenic potential in neocortical precursor cells

Neural precursor cells (NPCs) in the mouse neocortex generate various neuronal and glial cell types in a developmental stage–dependent manner. Most NPCs lose their neurogenic potential during development, although the underlying mechanisms of this process are not fully understood. We found that the chromatin of mouse NPCs gradually becomes more condensed and less dynamic on a global scale during neocortical development. Furthermore, we found high mobility group A (HMGA) proteins to be essential for the open chromatin state of NPCs at early developmental stages. Knockdown of HMGA proteins in early-stage NPCs reduced their neurogenic potential. Conversely, overexpression of HMGA proteins conferred neurogenic potential on late-stage NPCs, an effect that was antagonized by coexpression of a histone H1 mutant that inhibits chromatin opening. Thus, HMGA proteins contribute to the neurogenic potential of NPCs in the early stages of neocortical development, possibly through induction of an open chromatin state.

The zebrafish IgH locus contains multiple transcriptional regulatory regions

Many fish have, in addition to IgM and IgD, a third isotype called IgZ or IgT. The ζ-chain locus is embedded among the Ig heavy chain V-, D- and J-elements in a manner reminiscent of the TcR δ/α locus. Isotype selection thus occurs during VDJ recombination, a process that is facilitated by intralocus transcription. Using in silico analyses and enhancer reporter vectors we identified 3 new regions within the zebrafish IgH locus through which transcription can be activated in catfish B-cell lines. Two of these, termed Eζi (Jζ to Cζ1 intronic) and Eζ3' regions flank the ζ-chain constant domain exons. A third region, Eδ3', resides downstream of the δ-chain exons. All regions contain predicted binding sites for transcription factors that contribute to B-cell specific transcription in fish and mammals. Each region also has proximal matrix attachment regions, which may further contribute to transcriptional activation and chromatin remodeling. We discuss possible roles for these regions during VDJ recombination.

Genome mapping and characterization of the Anopheles gambiae heterochromatin

BACKGROUND

Heterochromatin plays an important role in chromosome function and gene regulation. Despite the availability of polytene chromosomes and genome sequence, the heterochromatin of the major malaria vector Anopheles gambiae has not been mapped and characterized.

RESULTS

To determine the extent of heterochromatin within the An. gambiae genome, genes were physically mapped to the euchromatin-heterochromatin transition zone of polytene chromosomes. The study found that a minimum of 232 genes reside in 16.6 Mb of mapped heterochromatin. Gene ontology analysis revealed that heterochromatin is enriched in genes with DNA-binding and regulatory activities. Immunostaining of the An. gambiae chromosomes with antibodies against Drosophila melanogaster heterochromatin protein 1 (HP1) and the nuclear envelope protein lamin Dm0 identified the major invariable sites of the proteins' localization in all regions of pericentric heterochromatin, diffuse intercalary heterochromatin, and euchromatic region 9C of the 2R arm, but not in the compact intercalary heterochromatin. To better understand the molecular differences among chromatin types, novel Bayesian statistical models were developed to analyze genome features. The study found that heterochromatin and euchromatin differ in gene density and the coverage of retroelements and segmental duplications. The pericentric heterochromatin had the highest coverage of retroelements and tandem repeats, while intercalary heterochromatin was enriched with segmental duplications. We also provide evidence that the diffuse intercalary heterochromatin has a higher coverage of DNA transposable elements, minisatellites, and satellites than does the compact intercalary heterochromatin. The investigation of 42-Mb assembly of unmapped genomic scaffolds showed that it has molecular characteristics similar to cytologically mapped heterochromatin.

CONCLUSIONS

Our results demonstrate that Anopheles polytene chromosomes and whole-genome shotgun assembly render the mapping and characterization of a significant part of heterochromatic scaffolds a possibility. These results reveal the strong association between characteristics of the genome features and morphological types of chromatin. Initial analysis of the An. gambiae heterochromatin provides a framework for its functional characterization and comparative genomic analyses with other organisms.

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.

Genome landscape and evolutionary plasticity of chromosomes in malaria mosquitoes

BACKGROUND

Nonrandom distribution of rearrangements is a common feature of eukaryotic chromosomes that is not well understood in terms of genome organization and evolution. In the major African malaria vector Anopheles gambiae, polymorphic inversions are highly nonuniformly distributed among five chromosomal arms and are associated with epidemiologically important adaptations. However, it is not clear whether the genomic content of the chromosomal arms is associated with inversion polymorphism and fixation rates.

METHODOLOGY/PRINCIPAL FINDINGS

To better understand the evolutionary dynamics of chromosomal inversions, we created a physical map for an Asian malaria mosquito, Anopheles stephensi, and compared it with the genome of An. gambiae. We also developed and deployed novel Bayesian statistical models to analyze genome landscapes in individual chromosomal arms An. gambiae. Here, we demonstrate that, despite the paucity of inversion polymorphisms on the X chromosome, this chromosome has the fastest rate of inversion fixation and the highest density of transposable elements, simple DNA repeats, and GC content. The highly polymorphic and rapidly evolving autosomal 2R arm had overrepresentation of genes involved in cellular response to stress supporting the role of natural selection in maintaining adaptive polymorphic inversions. In addition, the 2R arm had the highest density of regions involved in segmental duplications that clustered in the breakpoint-rich zone of the arm. In contrast, the slower evolving 2L, 3R, and 3L, arms were enriched with matrix-attachment regions that potentially contribute to chromosome stability in the cell nucleus.

CONCLUSIONS/SIGNIFICANCE

These results highlight fundamental differences in evolutionary dynamics of the sex chromosome and autosomes and revealed the strong association between characteristics of the genome landscape and rates of chromosomal evolution. We conclude that a unique combination of various classes of genes and repetitive DNA in each arm, rather than a single type of repetitive element, is likely responsible for arm-specific rates of rearrangements.

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.

AGAMOUS controls GIANT KILLER, a multifunctional chromatin modifier in reproductive organ patterning and differentiation

The floral, homeotic protein AGAMOUS coordinates multiple downstream genes through direct transcriptional regulation of the nuclear matrix attachment region binding protein GIANT KILLER.

The Arabidopsis homeotic protein AGAMOUS (AG), a MADS domain transcription factor, specifies reproductive organ identity during flower development. Using a binding assay and expression analysis, we identified a direct target of AG, GIANT KILLER (GIK), which fine-tunes the expression of multiple genes downstream of AG. The GIK protein contains an AT-hook DNA binding motif that is widely found in chromosomal proteins and that binds to nuclear matrix attachment regions of DNA elements. Overexpression and loss of function of GIK cause wide-ranging defects in patterning and differentiation of reproductive organs. GIK directly regulates the expression of several key transcriptional regulators, including ETTIN/AUXIN RESPONSE FACTOR 3 (ETT/ARF3) that patterns the gynoecium, by binding to the matrix attachment regions of target promoters. Overexpression of GIK causes a swift and dynamic change in repressive histone modification in the ETT promoter. We propose that GIK acts as a molecular node downstream of the homeotic protein AG, regulating patterning and differentiation of reproductive organs through chromatin organization.

Multicellular development depends on proper expression of thousands of genes. Master regulators, such as homeotic proteins, code for transcription factors in both plants and animals and are thought to act by regulating other genes. Recent genomic studies in the plant Arabidopsis have shown that over 1,000 genes are regulated by homeotic proteins that directly control various target genes, including different classes of transcriptional regulators. It is not known, however, how expression of so many genes is coordinated by a single homeotic gene to form functional organs and tissues. Here we identified a transcriptional target of the plant homeotic protein AGAMOUS using bioinformatics analysis and showed that AGAMOUS directly controls GIANT KILLER, a multifunctional chromatin modifier. GIANT KILLER then binds to the upstream regions of multiple genes involved in patterning and differentiation in the AGAMOUS pathway and fine-tunes the expression of these genes. These data therefore provide a possible mechanism by which a homeotic gene coordinates multiple downstream targets in plants.

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.

In silico comparative genomic analysis of GABAA receptor transcriptional regulation

Background

Subtypes of the GABA_A receptor subunit exhibit diverse temporal and spatial expression patterns. In silico comparative analysis was used to predict transcriptional regulatory features in individual mammalian GABA_A receptor subunit genes, and to identify potential transcriptional regulatory components involved in the coordinate regulation of the GABA_A receptor gene clusters.

Results

Previously unreported putative promoters were identified for the β2, γ1, γ3, ε, θ and π subunit genes. Putative core elements and proximal transcriptional factors were identified within these predicted promoters, and within the experimentally determined promoters of other subunit genes. Conserved intergenic regions of sequence in the mammalian GABA_A receptor gene cluster comprising the α1, β2, γ2 and α6 subunits were identified as potential long range transcriptional regulatory components involved in the coordinate regulation of these genes. A region of predicted DNase I hypersensitive sites within the cluster may contain transcriptional regulatory features coordinating gene expression. A novel model is proposed for the coordinate control of the gene cluster and parallel expression of the α1 and β2 subunits, based upon the selective action of putative Scaffold/Matrix Attachment Regions (S/MARs).

Conclusion

The putative regulatory features identified by genomic analysis of GABA_A receptor genes were substantiated by cross-species comparative analysis and now require experimental verification. The proposed model for the coordinate regulation of genes in the cluster accounts for the head-to-head orientation and parallel expression of the α1 and β2 subunit genes, and for the disruption of transcription caused by insertion of a neomycin gene in the close vicinity of the α6 gene, which is proximal to a putative critical S/MAR.

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.

The influence of matrix attachment regions on transgene expression in Arabidopsis thaliana wild type and gene silencing mutants

Many studies in both animal and plant systems have shown that matrix attachment regions (MARs) can increase the expression of flanking transgenes. However, our previous studies revealed no effect of the chicken lysozyme MARs (chiMARs) on transgene expression in the first generation transgenic Arabidopsis thaliana plants transformed with a beta-glucuronidase gene (uidA) unless gene silencing mutants were used as genetic background for transformation. In the present study, we investigated why chiMARs do not influence transgene expression in transgenic wild-type Arabidopsis plants. We first studied the effect of chiMARs on transgene expression in the progeny of primary transformants harboring chiMAR-flanked T-DNAs. Our data indicate that chiMARs do not affect transgene expression in consecutive generations of wild-type A. thaliana plants. Next, we examined whether these observed results in A. thaliana transformants are influenced by the applied transformation method. The results from in vitro transformed A. thaliana plants are in accordance with those from in planta transformed A. thaliana plants and again reveal no influence of chiMARs on transgene expression in A. thaliana wild-type transformants. The effect of chi-MARs on transgene expression is also examined in in vitro transformed Nicotiana tabacum plants, but as for A. thaliana, the transgene expression in tobacco transformants is not altered by the presence of chi-MARs. Taken together, our results show that the applied method or the plant species used for transformation does not influence whether and how chiMARs have an effect on transgene expression. Finally, we studied the effect of MARs (tabMARs) of plant origin (tobacco) on the transgene expression in A. thaliana wild-type plants and suppressed gene silencing (sgs2) mutants. Our results clearly show that similar to chiMARs, the tobacco-derived MARs do not enhance transgene expression in a wild-type background but can be used to enhance transgene expression in a mutant impaired in gene silencing.

Chimeric Mos1 and piggyBac transposases result in site-directed integration

Genetic transformation systems based on Mos1 and piggyBac transposable elements are used to achieve stable chromosomal integration. However, integration sites are randomly distributed in the genome and transgene expression can be influenced by position effects. We developed a novel technology that utilizes chimeric transposases to direct integration into specific sites on a target DNA molecule. The Gal4 DNA binding domain was fused to the NH(2) terminus of the Mos1 and piggyBac transposases and a target plasmid was created that contained upstream activating sequences (UAS), to which the Gal4 DBD binds with high affinity. The transpositional activity of the Gal4-Mos1 transposase was 12.7-fold higher compared to controls where the Gal4-UAS interaction was absent and 96% of the recovered transposition products were identical, with integration occurring at the same TA site. In a parallel experiment, a Gal4-piggyBac transposase resulted in an 11.6-fold increase in transpositional activity compared to controls, with 67% of the integrations occurring at a single TTAA site. This technology has the potential to minimize nonspecific integration events that may result in insertional mutagenesis and reduced fitness. Site-directed integration will be advantageous to the manipulation of genomes, study of gene function, and for the development of gene therapy techniques.

Spatiotemporal expression control correlates with intragenic scaffold matrix attachment regions (S/MARs) in Arabidopsis thaliana

Scaffold/matrix attachment regions (S/MARs) are essential for structural organization of the chromatin within the nucleus and serve as anchors of chromatin loop domains. A significant fraction of genes in Arabidopsis thaliana contains intragenic S/MAR elements and a significant correlation of S/MAR presence and overall expression strength has been demonstrated. In this study, we undertook a genome scale analysis of expression level and spatiotemporal expression differences in correlation with the presence or absence of genic S/MAR elements. We demonstrate that genes containing intragenic S/MARs are prone to pronounced spatiotemporal expression regulation. This characteristic is found to be even more pronounced for transcription factor genes. Our observations illustrate the importance of S/MARs in transcriptional regulation and the role of chromatin structural characteristics for gene regulation. Our findings open new perspectives for the understanding of tissue- and organ-specific regulation of gene expression.

Scaffold/matrix attachment regions (S/MARs) are AT-rich DNA sequences that mediate structural organization of the chromatin within the nucleus. These elements constitute anchor points of the DNA for the chromatin scaffold and serve to organize the chromatin into structural domains. Studies on individual genes led to the conclusion that the dynamic and complex organization of the chromatin mediated by S/MAR elements plays an important role in the regulation of gene expression. In addition to intergenic S/MARs, which likely exert import insulator effects, more than 2,000 intragenic S/MARs have been shown to be present within the Arabidopsis genome. In this study, the authors set out to analyze the effects of these intragenic S/MAR elements on the regulation of the genes affected. Making use of exhaustive and multidimensional expression datasets available for Arabidopsis, the authors analyzed overall expression differences and correlation of intragenic S/MARs with spatiotemporal expression of genes. On a genome scale, pronounced tissue- and organ-specific and developmental expression patterns of S/MAR-containing genes have been detected. Notably, transcription factor genes contain a significant higher portion of S/MARs. The pronounced difference in expression characteristics of S/MAR-containing genes emphasizes their functional importance and the importance of structural chromosomal characteristics for gene regulation in plants as well as within other eukaryotes.

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.

The pattern of chromosome folding in interphase is outlined by the linear gene density profile

Spatial organization of chromatin in the interphase nucleus plays a role in gene expression and inheritance. Although it appears not to be random, the principles of this organization are largely unknown. In this work, we show an explicit relationship between the intranuclear localization of various chromosome segments and the pattern of gene distribution along the genome sequence. Using a 7-megabase-long region of the Drosophila melanogaster chromosome 2 as a model, we observed that the six gene-poor chromosome segments identified in the region interact with components of the nuclear matrix to form a compact stable cluster. The six gene-rich segments form a spatially segregated unstable cluster dependent on nonmatrix nuclear proteins. The resulting composite structure formed by clusters of gene-rich and gene-poor regions is reproducible between the nuclei. We suggest that certain aspects of chromosome folding in interphase are predetermined and can be inferred through in silico analysis of chromosome sequence, using gene density profile as a manifestation of "folding code."

Osteo-Promoter Database (OPD) -- promoter analysis in skeletal cells

Background

Increasing our knowledge about the complex expression of genes in skeletal tissue will provide a better understanding of the physiology of skeletal cells. The study summarizes transcriptional regulation factors interacting and cooperating at promoter regions that regulate gene expression. Specifically, we analyzed A/T rich elements along the promoter sequences.

Description

The Osteo-Promoter Database (OPD) is a collection of genes and promoters expressed in skeletal cells. We have compiled a new viewer, OPD, as unique database developed and created as an accessible tool for skeletal promoter sequences. OPD can navigate to identify genes specific to skeletal cDNA databases and promoter analysis sites. OPD offers exclusive access to facilitate a dynamic extraction of promoters' gene-specific analyses in skeletal tissue. The data on promoters included in OPD contains cloned promoters or predicted promoters that were analyzed by bioinformatics tools. OPD offers MAR-analysis, which allocates A/T rich elements along these promoter sequences.

Conclusion

The analysis leads to a better insight of proteins that bind to DNA, regulate DNA, and function in chromatin remodeling. The OPD is a distinctive tool for understanding the complex function of chromatin remodeling and transcriptional regulation of specific gene expression in skeletal tissue.

Evolution of transcriptional enhancers in the immunoglobulin heavy-chain gene: functional characteristics of the zebrafish Emu3' enhancer

Past studies of the channel catfish immunoglobulin heavy-chain (IgH) locus indicates that it lacks an Emu enhancer in the J(H)-Cmu1 intron but does have an enhancer, termed Emu3', in the mu-delta intergenic region. The positioning of the catfish enhancer downstream of the mu-chain exons is predicted to be unfavorable for antibody-affinity maturation in catfish, and would also have been an impediment to the evolution of class switch recombination, had it existed in early tetrapods. To determine if this downstream enhancer is a general feature of teleost fish, we have identified the location of the transcriptional enhancer in the zebrafish IgH locus. We find that zebrafish, like catfish, only have an Emu3'-like enhancer that has cross-species activity, but which is B-cell-specific in its activity. A 300-bp region of the zebrafish enhancer shares sequence homology with the core of the catfish Emu3', although there has been loss and gain of specific octamer enhancer motifs. Mutagenesis studies demonstrate that the zebrafish IgH enhancer depends on a pair of E-box motifs that are found in the enhancer core. Similarly spaced E-box motifs appear to exist in the Emu3' enhancer regions of other teleost fish, suggesting this is a common feature among fish IgH enhancers. We discuss how this distal positioning of the enhancer may influence affinity maturation in extant teleosts as well as the evolution of this process in the early tetrapods.

Towards understanding the epigenetics of transcription by chromatin structure and the nuclear matrix

The eukaryotic nucleus houses a significant amount of information that is carefully ordered to ensure that genes can be transcribed as needed throughout development and differentiation. The genome is partitioned into regions containing functional transcription units, providing the means for the cell to selectively activate some, while keeping other regions of the genome silent. Over the last quarter of a century the structure of chromatin and how it is influenced by epigenetics has come into the forefront of modern biology. However, it has thus far failed to identify the mechanism by which individual genes or domains are selected for expression. Through covalent and structural modification of the DNA and chromatin proteins, epigenetics maintains both active and silent chromatin states. This is the "other" genetic code, often superseding that dictated by the nucleotide sequence. The nuclear matrix is rich in many of the factors that govern nuclear processes. It includes a host of unknown factors that may provide our first insight into the structural mechanism responsible for the genetic selectivity of a differentiating cell. This review will consider the nuclear matrix as an integral component of the epigenetic mechanism.

Genome-wide in silico mapping of scaffold/matrix attachment regions in Arabidopsis suggests correlation of intragenic scaffold/matrix attachment regions with gene expression

We carried out a genome-wide prediction of scaffold/matrix attachment regions (S/MARs) in Arabidopsis. Results indicate no uneven distribution on the chromosomal level but a clear underrepresentation of S/MARs inside genes. In cases where S/MARs were predicted within genes, these intragenic S/MARs were preferentially located within the 5'-half, most prominently within introns 1 and 2. Using Arabidopsis whole-genome expression data generated by the massively parallel signature sequencing methodology, we found a negative correlation between S/MAR-containing genes and transcriptional abundance. Expressed sequence tag data correlated the same way with S/MAR-containing genes. Thus, intragenic S/MARs show a negative correlation with transcription level. For various genes it has been shown experimentally that S/MARs can function as transcriptional regulators and that they have an implication in stabilizing expression levels within transgenic plants. On the basis of a genome-wide in silico S/MAR analysis, we found a significant correlation between the presence of intragenic S/MARs and transcriptional down-regulation.

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.

Structure-function analysis of the estrogen receptor alpha corepressor scaffold attachment factor-B1: identification of a potent transcriptional repression domain

Scaffold attachment factor-B1 (SAFB1) is a nuclear matrix protein that has been proposed to couple chromatin structure, transcription, and RNA processing. We have previously shown that SAFB1 can repress estrogen receptor (ERalpha)-mediated transactivation. Here we present a structure-function study showing that transactivation is mediated via an intrinsic and transferable C-terminal repression domain (RD). A similar C-terminal RD was found in the family member SAFB2. Removal of the RD from SAFB1 resulted in a dominant-negative SAFB1 protein that increased ligand-dependent and -independent ERalpha activity. SAFB1RD-mediated repression was partly blocked by histone deacetylase inhibitors; however, no histone deacetylase inhibitors were identified in a yeast two-hybrid screen using the RD as bait. Instead, SAFB1RD was found to interact with TAFII68, a member of the basal transcription machinery. We propose a model in which SAFB1 represses ERalpha activity via indirect association with histone deacetylation and interaction with the basal transcription machinery.

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.

Nuclear matrix association of the human beta-globin locus utilizing a novel approach to quantitative real-time PCR

The human beta-globin locus is home to five genes that are regulated in a tissue-specific and developmental stage-specific manner. While the exact mode of expression remains somewhat enigmatic, a significant effort has been focused at the locus control region (LCR). The LCR is marked by five DNase I-hypersensitive sites (HS) approximately 15 kb upstream of the epsilon-globin gene. Nuclear matrix-associated regions (MARs) organize chromatin into functional domains and at least one of the HS appears bound to the nuclear matrix. We have employed an in vivo based PCR MAR assay to investigate the role of MAR-mediated regulation of the beta-globin locus. This was facilitated with a novel reaction efficiency based quantitative real-time PCR analysis software tool, Target Analysis Quantification. Using a log-linear regression strategy, discordances were eliminated. This allowed us to reliably estimate the relative amount of initial template associated with the nuclear matrix at 15 unique regions spanning the beta-globin locus in both non-expressing and expressing cell lines. A dynamic association dependent on expression status was revealed both at the LCR/5'HS region and within the second intron of the beta-globin gene. These results provide the first evidence that nuclear matrix association dynamically mediates the looping of the beta-globin locus to achieve transcriptional control.

In silico identification of metazoan transcriptional regulatory regions

Transcriptional regulation remains one of the most intriguing and challenging subjects in biomedical research. The catalysis of transcription is a clear example of multiple proteins interacting to orchestrate a biological process, offering a starting point for the study of biological systems. Transcriptional regulation is viewed as one of the principal mechanisms governing the spatial and temporal distribution of gene expression, thus the field of transcriptional regulation provides a natural stage for quantitative studies of multiple gene systems. Building on the body of focused experimental studies and new genomics-driven data, computational biologists are making significant strides in accelerating our understanding of the transcriptional regulatory process in metazoan cells. Recent advances in the computational analysis of the interplay between factors have been fueled by well-defined computational methods for the modeling of the binding of individual transcription factors. We present here an overview of advances in the analysis of regulatory systems and the fundamental methods that underlie the recent developments.

Molecular structure and regulatory potential of a T-DNA integration site in petunia

The genomic structure surrounding a T-DNA integration site in a transgenic petunia plant, which shows deregulation of a root-specific promoter, was investigated. We have already demonstrated that T-DNA integration in this transformant (P13) had occurred close to a scaffold/matrix attachment region (S/MAR). A major question regarding the observed promoter leakiness was whether the T-DNA had integrated into the centre or at the border of the Petun-SAR and whether other regulatory elements are located within this genomic region. While small rearrangements were shown to occur during T-DNA integration in agreement with other reports, we find indications of the presence of a SINE retroposon--an apparent landmark for recombinogenic targets--at the integration site. Binding assays to both plant and animal nuclear scaffolds, supported by biomathematical analyses, reveal that the T-DNA is definitely located at the border of a strong S/MAR, which is in agreement with current models on the structure of integration sites. These results, together with a developmentally regulated leaf-specific enhancer effect of the Petun-SAR on gene expression in transgenic tobacco plants, indicate that the Petun-SAR demarcates the right border of a chromatin domain with genes predominantly active in leaves.

Evaluation of sequence motifs found in scaffold/matrix-attached regions (S/MARs)

Based on the contents of the database S/MARt DB, the most comprehensive data collection of scaffold/matrix-attached regions (S/MARs) publicly available thus far, we initiated a systematic evaluation of the stored data. By analyzing the 245 S/MAR sequences presently described in this database, we found that the S/MARs contained in this collection are generally AT-rich, with certain significant exceptions. Comparative analyses showed that most of the AT-rich motifs which were found to be enriched in S/MARs are also enriched in randomized S/MAR sequences of the same AT content. Some sequence patterns previously suggested to be characteristic for S/MARs were also investigated, among them potential binding sites for homeodomain transcription factors. Even though hexanucleotides containing the core motif of homeodomain factors were frequently observed in S/MARs, only a few potential binding sites for these factors were found enriched when compared with regulatory regions or exon sequences. All our analyses indicated that, on average, the observed frequency of motifs in S/MAR elements is largely influenced by the AT content. Our results can serve as a guideline for further improvements in the definition of S/MARs, which are now believed to constitute the functional coordinate system for genomic regulatory regions.