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

Transposable elements as scaffold/matrix attachment regions: shaping organization and functions in genomes

The hierarchical structure of eukaryotic genomes has regulatory layers, one of them being epigenetic "indexing" of the genome that leads to cell-type-specific patterns of gene expression. By establishing loops and defining chromatin domains, cells can achieve coordinated control over multi-locus segments of the genome. This is thought to be achieved using scaffold/matrix attachment regions (S/MARs) that establish structural and functional loops and topologically associating domains (TADs) that define a self-interacting region of the genome. Large-scale genome-wide mapping of S/MARs has begun to uncover these aspects of genome organization. A recent genome-wide study showed the association of transposable elements (TEs) with a significant fraction of S/MARs, suggesting that the multitude of TE-derived repeats constitute a class of anchorage sites of chromatin loops to nuclear architecture. In this study, we provide an insight that TE-driven dispersal of S/MARs has the potential to restructure the chromosomes by creating novel loops and domains. The combination of TEs and S/MARs, as elements that can hop through the genome along with regulatory capabilities, may provide an active mechanism of genome evolution leading to the emergence of novel features in biological systems. The significance is that a genome-wide study mapping developmental S/MARs reveals an intriguing link between these elements and TEs. This article highlights the potential of the TE-S/MAR combination to drive evolution by restructuring and shaping the genome.

Selective DNA-binding of SP120 (rat ortholog of human hnRNP U) is mediated by arginine-glycine rich domain and modulated by RNA

A human protein heterogeneous ribonucleoprotein U (hnRNP U) also known as Scaffold attachment factor A (SAF-A) and its orthologous rat protein SP120 are abundant and multifunctional nuclear protein that directly binds to both DNA and RNA. The C-terminal region of hnRNP U enriched with arginine and glycine is essential for the interaction with RNA and the N-terminal region of SAF-A termed SAP domain has been ascribed to the DNA binding. We have reported that rat hnRNP U specifically and cooperatively binds to AT-rich DNA called nuclear scaffold/matrix-associated region (S/MAR) although its detailed mechanism remained unclear. In the present study analysis of hnRNP U deletion mutants revealed for the first time that a C-terminal domain enriched with Arg-Gly (defined here as 'RG domain') is predominantly important for the S/MAR-selective DNA binding activities. RG domain alone directly bound to S/MAR and coexistence with the SAP domain exerted a synergistic effect. The binding was inhibited by netropsin, a minor groove binder with preference to AT pairs that are enriched in S/MAR, suggesting that RG domain interacts with minor groove of S/MAR DNA. Interestingly, excess amounts of RNA attenuated the RG domain-dependent S/MAR-binding of hnRNP U. Taken together, hnRNP U may be the key element for the RNA-regulated recognition of S/MAR DNA and thus contributing to the dynamic structural changes of chromatin compartments.

Domain Model of Eukaryotic Genome Organization: From DNA Loops Fixed on the Nuclear Matrix to TADs

The article reviews the development of ideas on the domain organization of eukaryotic genome, with special attention on the studies of DNA loops anchored to the nuclear matrix and their role in the emergence of the modern model of eukaryotic genome spatial organization. Critical analysis of results demonstrating that topologically associated chromatin domains are structural-functional blocks of the genome supports the notion that these blocks are fundamentally different from domains whose existence was proposed by the domain hypothesis of eukaryotic genome organization formulated in the 1980s. Based on the discussed evidence, it is concluded that the model postulating that eukaryotic genome is built from uniformly organized structural-functional blocks has proven to be untenable.

Dynamics of nuclear matrix attachment regions during 5th instar posterior silk gland development in Bombyx mori

Background

Chromatin architecture is critical for gene expression during development. Matrix attachment regions (MARs) control and regulate chromatin dynamics. The position of MARs in the genome determines the expression of genes in the organism. In this study, we set out to elucidate how MARs temporally regulate the expression of the fibroin heavy chain (FIBH) gene during development. We addressed this by identifying MARs and studying their distribution and differentiation, in the posterior silk glands of Bombyx mori during 5th instar development.

Results

Of the MARs identified on three different days, 7.15% MARs were common to all 3 days, whereas, 1.41, 19.27 and 52.47% MARs were unique to day 1, day 5, and day 7, respectively highlighting the dynamic nature of the matrix associated DNA. The average chromatin loop length based on the chromosome wise distribution of MARs and the distances between these MAR regions decreased from day 1 (253.91 kb) to day 5 (73.54 kb) to day 7 (39.19 kb). Further significant changes in the MARs in the vicinity of the FIBH gene were found during different days of 5^th instar development which implied their role in the regulation and expression of the FIBH gene.

Conclusions

The presence of MARs in the flanking regions of genes found to exhibit differential expression during 5^th instar development indicates their possible role in the regulation of their expression. This reiterates the importance of MARs in the genomic functioning as regulators of the molecular mechanisms in the nucleus. This is the first study that takes into account the tissue specific genome-wide MAR association and the potential role of these MARs in developmentally regulated gene expression. The current study lays a foundation to understand the genome wide regulation of chromatin during development.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08446-3.

Specific subfamilies of transposable elements contribute to different domains of T lymphocyte enhancers

Transposable elements (TEs) compose nearly half of mammalian genomes and provide building blocks for cis-regulatory elements. Using high-throughput sequencing, we show that 84 TE subfamilies are overrepresented, and distributed in a lineage-specific fashion in core and boundary domains of CD8⁺ T cell enhancers. Endogenous retroviruses are most significantly enriched in core domains with accessible chromatin, and bear recognition motifs for immune-related transcription factors. In contrast, short interspersed elements (SINEs) are preferentially overrepresented in nucleosome-containing boundaries. A substantial proportion of these SINEs harbor a high density of the enhancer-specific histone mark H3K4me1 and carry sequences that match enhancer boundary nucleotide composition. Motifs with regulatory features are better preserved within enhancer-enriched TE copies compared to their subfamily equivalents located in gene deserts. TE-rich and TE-poor enhancers associate with both shared and unique gene groups and are enriched in overlapping functions related to lymphocyte and leukocyte biology. The majority of T cell enhancers are shared with other immune lineages and are accessible in common hematopoietic progenitors. A higher proportion of immune tissue-specific enhancers are TE-rich compared to enhancers specific to other tissues, correlating with higher TE occurrence in immune gene-associated genomic regions. Our results suggest that during evolution, TEs abundant in these regions and carrying motifs potentially beneficial for enhancer architecture and immune functions were particularly frequently incorporated by evolving enhancers. Their putative selection and regulatory cooption may have accelerated the evolution of immune regulatory networks.

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.

Comparative Analysis of Intra- and Inter-Specific Genomic Variability in the Peach Potato Aphid, Myzus persicae

The availability of genomic data in the last decade relating to different aphid species has allowed the analysis of the genomic variability occurring among such species, whereas intra-specific variability has hitherto very largely been neglected. In order to analyse the intra-genomic variability in the peach potato aphid, Myzus persicae, comparative analyses were performed revealing several clone-specific gene duplications, together with numerous deletions/rearrangements. Our comparative approach also allowed us to evaluate the synteny existing between the two M. persicae clones tested and between the peach potato aphid and the pea aphid, Acyrthosiphon pisum. Even if part of the observed rearrangements are related to a low quality of some assembled contigs and/or to the high number of contigs present in these aphid genomes, our evidence reveals that aphid clones are genetically more different than expected. These results suggest that the choice of performing genomes sequencing combining different biotypes/populations, as revealed in the case of the soybean aphid, Aphis glycines, is unlikely to be very informative in aphids. Interestingly, it is possible that the holocentric nature of aphid chromosomes favours genome rearrangements that can be successively inherited transgenerationally via the aphid's apomictic (parthenogenetic) mode of reproduction. Lastly, we evaluated the structure of the cluster of genes coding for the five histones (H1, H2A, H2B, H3 and H4) in order to better understand the quality of the two M. persicae genomes and thereby to improve our knowledge of this functionally important gene family.

Development of an improved lentiviral based vector system for the stable expression of monoclonal antibody in CHO cells

Therapeutic monoclonal antibodies (mAbs) have become the dominant products in biopharmaceutical industry. Mammalian cell expression systems including Chinese hamster ovary (CHO) cells are the most commonly used hosts for the production of complex recombinant proteins. However, development of stable, high producing CHO cell lines suffers from the low expression level and instability of the transgene. The increasing efforts in the development of novel therapeutic antibodies and the advent of biosimilars have revealed the necessity for the development of improved platforms for rapid production of products for initial characterization and testing. In line with this premise, vector design and engineering has been applied to improve the expression level and stability of the transgene. This study reports the application of an improved lentiviral vector system containing the human interferon-β scaffold attachment region (IFN-SAR) for the development of antibody producing stable CHO cells. mAb expressing clones producing 1100 µg/L of IgG1 monoclonal antibody were isolated without extensive screening of a large number of clones. Our results here indicate the positive effects of IFN-SAR on stable mAb expression using lentiviral based expression vectors. We also observed that although IFN-SAR can improve light chain (LC) and heavy chain (HC) gene copy numbers in stable cell pools, mAb expression in single cell clones was not affected by the transgene copy number.

Chimeric 3' flanking regions strongly enhance gene expression in plants

Plants represent a promising platform for the highly scalable production of recombinant proteins. Previously, we identified the tobacco extensin terminator lacking its intron as an element that reduced transcript read-through and improved recombinant protein production in a plant-based system. In this study, we systematically compared nonreplicating plant expression vectors containing over 20 commonly used or newly identified terminators from diverse sources. We found that eight gene terminators enhance reporter gene expression significantly more than the commonly used 35S and NOS terminators. The intronless extensin terminator provided a 13.6-fold increase compared with the NOS terminator. Combining terminators in tandem produced large synergistic effects, with many combinations providing a >25-fold increase in expression. Addition of the tobacco Rb7 or TM6 matrix attachment region (MAR) strongly enhanced protein production when added to most terminators, with the Rb7 MAR providing the greatest enhancement. Using deletion analysis, the full activity of the 1193 bp Rb7 MAR was found to require only a 463-bp region at its 3' end. Combined terminators and MAR together provided a >60-fold increase compared with the NOS terminator alone. These combinations were then placed in a replicating geminiviral vector, providing a total of >150-fold enhancement over the original NOS vector, corresponding to an estimated yield of 3-5 g recombinant protein per kg leaf fresh weight or around 50% of the leaf total soluble protein. These results demonstrate the importance of 3' flanking regions in optimizing gene expression and show great potential for 3' flanking regions to improve DNA-based recombinant protein production systems.

How mutations shape p53 interactions with the genome to promote tumorigenesis and drug resistance

The tumor suppressive transcription factor p53 regulates a wide array of cellular processes that confer upon cells an essential protection against cancer development. Wild-type p53 regulates gene expression by directly binding to DNA in a sequence-specific manner. p53 missense mutations are the most common mutations in malignant cells and can be regarded as synonymous with anticancer drug resistance and poor prognosis. The current review provides an overview of how the extraordinary variety of more than 2000 different mutant p53 proteins, known as the p53 mutome, affect the interaction of p53 with DNA. We discuss how the classification of p53 mutations to loss of function (LOF), gain of function (GOF), and dominant-negative (DN) inhibition of a remaining wild-type allele, hides a complex p53 mutation spectrum that depends on the distinctive nature of each mutant protein, requiring different therapeutic strategies for each mutant p53 protein. We propose to regard the different mutant p53 categories as continuous variables, that may not be independent of each other. In particular, we suggest here to consider GOF mutations as a special subset of LOF mutations, especially when mutant p53 binds to DNA through cooperation with other transcription factors, and we present a model for GOF mechanism that consolidates many observations on the GOF phenomenon. We review how novel mutant p53 targeting approaches aim to restore a wild-type-like DNA interaction and to overcome resistance to cancer therapy.

Sus scrofa matrix attachment region enhances expression of the PiggyBac system transfected into HEK293T cells

OBJECTIVES

To determine the effects of the Sus scrofa matrix attachment region (SusMAR) on transgene expression in HEK293T cells.

RESULTS

Three expression vectors with the MAR at different sites in the PiggyBac (PB) transposon vector backbone were compared: two MARs flanking the β-galactosidase (β-gal) expression cassette, and one at the upstream or downstream site. Bos taurus MAR (BosMAR) and a β-gal expression cassette without MARs were the positive and negative controls, respectively. Compared to the control, β-gal activity of all SusMAR and BosMAR vectors was significantly improved in the presence of PB transposase (PBase). However, only the downstream SusMAR and upstream BosMAR vectors showed increased expression in the absence of PBase. Expression was significantly increased in all vectors with the PBase group compared to those without the PBase group. Gene copy numbers were not increased compared to the negative control.

CONCLUSIONS

SusMAR enhanced recombinant gene expression levels and stability in HEK293T cells, was not increase transgene copy number. These results could facilitate the development of vectors for stable production of therapeutic proteins.

5' and 3' Untranslated Regions Strongly Enhance Performance of Geminiviral Replicons in Nicotiana benthamiana Leaves

We previously reported a recombinant protein production system based on a geminivirus replicon that yields high levels of vaccine antigens and monoclonal antibodies in plants. The bean yellow dwarf virus (BeYDV) replicon generates massive amounts of DNA copies, which engage the plant transcription machinery. However, we noticed a disparity between transcript level and protein production, suggesting that mRNAs could be more efficiently utilized. In this study, we systematically evaluated genetic elements from human, viral, and plant sources for their potential to improve the BeYDV system. The tobacco extensin terminator enhanced transcript accumulation and protein production compared to other commonly used terminators, indicating that efficient transcript processing plays an important role in recombinant protein production. Evaluation of human-derived 5' untranslated regions (UTRs) indicated that many provided high levels of protein production, supporting their cross-kingdom function. Among the viral 5' UTRs tested, we found the greatest enhancement with the tobacco mosaic virus omega leader. An analysis of the 5' UTRs from the Arabidopsis thaliana and Nicotinana benthamiana photosystem I K genes found that they were highly active when truncated to include only the near upstream region, providing a dramatic enhancement of transgene production that exceeded that of the tobacco mosaic virus omega leader. The tobacco Rb7 matrix attachment region inserted downstream from the gene of interest provided significant enhancement, which was correlated with a reduction in plant cell death. Evaluation of Agrobacterium strains found that EHA105 enhanced protein production and reduced cell death compared to LBA4301 and GV3101. We used these improvements to produce Norwalk virus capsid protein at >20% total soluble protein, corresponding to 1.8 mg/g leaf fresh weight, more than twice the highest level ever reported in a plant system. We also produced the monoclonal antibody rituximab at 1 mg/g leaf fresh weight.

A role for chromatin topology in imprinted domain regulation

Recently, many advancements in genome-wide chromatin topology and nuclear architecture have unveiled the complex and hidden world of the nucleus, where chromatin is organized into discrete neighbourhoods with coordinated gene expression. This includes the active and inactive X chromosomes. Using X chromosome inactivation as a working model, we utilized publicly available datasets together with a literature review to gain insight into topologically associated domains, lamin-associated domains, nucleolar-associating domains, scaffold/matrix attachment regions, and nucleoporin-associated chromatin and their role in regulating monoallelic expression. Furthermore, we comprehensively review for the first time the role of chromatin topology and nuclear architecture in the regulation of genomic imprinting. We propose that chromatin topology and nuclear architecture are important regulatory mechanisms for directing gene expression within imprinted domains. Furthermore, we predict that dynamic changes in chromatin topology and nuclear architecture play roles in tissue-specific imprint domain regulation during early development and differentiation.

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.

Chromatin structure and organization: the relation with gene expression during development and disease

The elementary level of chromatin fiber, namely the nucleofilament, is known to undergo a hierarchical compaction leading to local chromatin loops, then chromatin domains and ultimately chromosome territories. These successive folding levels rely on the formation of chromatin loops ranging from few kb to some Mb. In addition to a packaging and structural role, the high-order organization of genomes functionally impacts on gene expression program. This review summarises to which extent each level of chromatin compaction does affect gene regulation. In addition, we point out the structural and functional changes observed in diseases. Emphasis will be mainly placed on the large-scale organization of the chromatin.

Evolutionary force of AT-rich repeats to trap genomic and episomal DNAs into the rice genome: lessons from endogenous pararetrovirus

In plant genomes, the incorporation of DNA segments is not a common method of artificial gene transfer. Nevertheless, various segments of pararetroviruses have been found in plant genomes in recent decades. The rice genome contains a number of segments of endogenous rice tungro bacilliform virus-like sequences (ERTBVs), many of which are present between AT dinucleotide repeats (ATrs). Comparison of genomic sequences between two closely related rice subspecies, japonica and indica, allowed us to verify the preferential insertion of ERTBVs into ATrs. In addition to ERTBVs, the comparative analyses showed that ATrs occasionally incorporate repeat sequences including transposable elements, and a wide range of other sequences. Besides the known genomic sequences, the insertion sequences also represented DNAs of unclear origins together with ERTBVs, suggesting that ATrs have integrated episomal DNAs that would have been suspended in the nucleus. Such insertion DNAs might be trapped by ATrs in the genome in a host-dependent manner. Conversely, other simple mono- and dinucleotide sequence repeats (SSR) were less frequently involved in insertion events relative to ATrs. Therefore, ATrs could be regarded as hot spots of double-strand breaks that induce non-homologous end joining. The insertions within ATrs occasionally generated new gene-related sequences or involved structural modifications of existing genes. Likewise, in a comparison between Arabidopsis thaliana and Arabidopsis lyrata, the insertions preferred ATrs to other SSRs. Therefore ATrs in plant genomes could be considered as genomic dumping sites that have trapped various DNA molecules and may have exerted a powerful evolutionary force.

Relationship between DNA replication and the nuclear matrix

There is an extensive list of primary published work related to the nuclear matrix (NM). Here we review the aspects that are required to understand its relationship with DNA replication, while highlighting some of the difficulties in studying such a structure, and possible differences that arise from the choice of model system. We consider NM attachment regions of DNA and discuss their characteristics and potential function before reviewing data that deal specifically with functional interaction with DNA replication factors. Data have long existed indicating that newly synthesized DNA is associated with a nuclease-resistant NM, allowing the conclusion that the elongation step of DNA synthesis is immobilized within the nucleus. We review in more detail the emerging data that suggest that prereplication complex proteins and origins of replication are transiently recruited to the NM during late G1 and early S-phase. Collectively, these data suggest that the initiation step of the DNA replication process is also immobilized by attachment to the NM. We outline models that discuss the possible spatial relationships and highlight the emerging evidence that suggests there may be important differences between cell types.

Mechanisms of unphosphorylated STAT3 transcription factor binding to DNA

Phosphorylation of signal transducer and activator of transcription 3 (STAT3) on a single tyrosine residue in response to growth factors, cytokines, interferons, and oncogenes activates its dimerization, translocation to the nucleus, binding to the interferon γ (gamma)-activated sequence (GAS) DNA-binding site and activation of transcription of target genes. STAT3 is constitutively phosphorylated in various cancers and drives gene expression from GAS-containing promoters to promote tumorigenesis. Recently, roles for unphosphorylated STAT3 (U-STAT3) have been described in response to cytokine stimulation, in cancers, and in maintenance of heterochromatin stability. However, the mechanisms underlying U-STAT3 binding to DNA has not been fully investigated. Here, we explore STAT3-DNA interactions by atomic force microscopy (AFM) imaging. We observed that U-STAT3 molecules bind to the GAS DNA-binding site as dimers and monomers. In addition, we observed that U-STAT3 binds to AT-rich DNA sequence sites and recognizes specific DNA structures, such as 4-way junctions and DNA nodes, within negatively supercoiled plasmid DNA. These structures are important for chromatin organization and our data suggest a role for U-STAT3 as a chromatin/genome organizer. Unexpectedly, we found that a C-terminal truncated 67.5-kDa STAT3 isoform recognizes single-stranded spacers within cruciform structures that also have a role in chromatin organization and gene expression. This isoform appears to be abundant in the nuclei of cancer cells and, therefore, may have a role in regulation of gene expression. Taken together, our data highlight novel mechanisms by which U-STAT3 binds to DNA and supports U-STAT3 function as a transcriptional activator and a chromatin/genomic organizer.

CTCF binds to sites in the major histocompatibility complex that are rapidly reconfigured in response to interferon-gamma

Activation of the major histocompatibility complex (MHC) by interferon-gamma (IFN−γ) is a fundamental step in the adaptive immune response to pathogens. Here, we show that reorganization of chromatin loop domains in the MHC is evident within the first 30 min of IFN−γ treatment of fibroblasts, and that further dynamic alterations occur up to 6 h. These very rapid changes occur at genomic sites which are occupied by CTCF and are close to IFN−γ-inducible MHC genes. Early responses to IFN−γ are thus initiated independently of CIITA, the master regulator of MHC class II genes and prepare the MHC for subsequent induction of transcription.

High-level transgene expression by homologous recombination-mediated gene transfer

Gene transfer and expression in eukaryotes is often limited by a number of stably maintained gene copies and by epigenetic silencing effects. Silencing may be limited by the use of epigenetic regulatory sequences such as matrix attachment regions (MAR). Here, we show that successive transfections of MAR-containing vectors allow a synergistic increase of transgene expression. This finding is partly explained by an increased entry into the cell nuclei and genomic integration of the DNA, an effect that requires both the MAR element and iterative transfections. Fluorescence in situ hybridization analysis often showed single integration events, indicating that DNAs introduced in successive transfections could recombine. High expression was also linked to the cell division cycle, so that nuclear transport of the DNA occurs when homologous recombination is most active. Use of cells deficient in either non-homologous end-joining or homologous recombination suggested that efficient integration and expression may require homologous recombination-based genomic integration of MAR-containing plasmids and the lack of epigenetic silencing events associated with tandem gene copies. We conclude that MAR elements may promote homologous recombination, and that cells and vectors can be engineered to take advantage of this property to mediate highly efficient gene transfer and expression.

A mini review of MAR-binding proteins

Genomic DNA encompasses several levels of organization, the nuclear matrix mediates the formation of DNA loop domains that are anchored to matrix attachment regions (MARs). By means of specific interaction with MAR binding proteins (MARBPs), MAR plays an important regulation role in enhancing transgene expression, decreasing expression variation among individuals of different transformants and serving as the replication origin. Through these years, some MARBPs have been identified and characterized from humans, plants, animals and algae so far and the list is growing. Most of MARBPs exist in a co-repressor/co-activator complex and involve in chromosome folding, regulation of gene expression, influencing cell development and inducing cell apoptosis. This review covers recent advances that have contributed to our understanding of MARBPs.

Common chromosomal fragile sites (CFS) may be involved in normal and traumatic cognitive stress memory consolidation and altered nervous system immunity

Previous reports of specific patterns of increased fragility at common chromosomal fragile sites (CFS) found in association with certain neurobehavioural disorders did not attract attention at the time due to a shift towards molecular approaches to delineate neuropsychiatric disorder candidate genes. Links with miRNA, altered methylation and the origin of copy number variation indicate that CFS region characteristics may be part of chromatinomic mechanisms that are increasingly linked with neuroplasticity and memory. Current reports of large-scale double-stranded DNA breaks in differentiating neurons and evidence of ongoing DNA demethylation of specific gene promoters in adult hippocampus may shed new light on the dynamic epigenetic changes that are increasingly appreciated as contributing to long-term memory consolidation. The expression of immune recombination activating genes in key stress-induced memory regions suggests the adoption by the brain of this ancient pattern recognition and memory system to establish a structural basis for long-term memory through controlled chromosomal breakage at highly specific genomic regions. It is furthermore considered that these mechanisms for management of epigenetic information related to stress memory could be linked, in some instances, with the transfer of the somatically acquired information to the germline. Here, rearranged sequences can be subjected to further selection and possible eventual retrotranscription to become part of the more stable coding machinery if proven to be crucial for survival and reproduction. While linkage of cognitive memory with stress and fear circuitry and memory establishment through structural DNA modification is proposed as a normal process, inappropriate activation of immune-like genomic rearrangement processes through traumatic stress memory may have the potential to lead to undesirable activation of neuro-inflammatory processes. These theories could have a significant impact on the interpretation of risks posed by heredity and the environment and the search for neuropsychiatric candidate genes.

Conservation of the PRM1→PRM2→TNP2 Domain

In mouse and human, the genes encoding protamines PRM1, PRM2 and transition protein TNP2 are found clustered together on chromosome 16. In addition, these three genes lie in the same orientation to one another and are coordinately expressed in a haploid-specific manner during spermatogenesis. Previously, we have shown that the human PRM1 --> PRM2 --> TNP2 locus exists as a single chromatin domain bounded by two male germ cell-specific MARs, i.e. Matrix Attachment Regions. A third, somatic-specific MAR element lies immediately 3' of the PRM1 --> PRM2 --> TNP2 domain. This MAR maps to a conserved CpG island 5' of the human SOCS-1 gene. Similarly, two candidate MARs flank the mouse Prm1 --> Prm2 --> Tnp2 domain. Comparative analysis of the mouse and human promoter regions identified several conserved regulatory motifs for each of the genes of this cluster. This further establishes the synteny of this region. Global structural similarities and the functional relevance of the associated candidate regulatory elements are discussed.

Reconfiguration of genomic anchors upon transcriptional activation of the human major histocompatibility complex

The folding of chromatin into topologically constrained loop domains is essential for genomic function. We have identified genomic anchors that define the organization of chromatin loop domains across the human major histocompatibility complex (MHC). This locus contains critical genes for immunity and is associated with more diseases than any other region of the genome. Classical MHC genes are expressed in a cell type-specific pattern and can be induced by cytokines such as interferon-gamma (IFNG). Transcriptional activation of the MHC was associated with a reconfiguration of chromatin architecture resulting from the formation of additional genomic anchors. These findings suggest that the dynamic arrangement of genomic anchors and loops plays a role in transcriptional regulation.

Anchoring the genome

Although the principles governing chromosomal architecture are largely unresolved, there is evidence that higher-order chromatin folding is mediated by the anchoring of specific DNA sequences to the nuclear matrix. These genome anchors are also crucial regulators of gene expression and DNA replication, and play a role in pathogenesis.

Structural families of genomic microsatellites

We present an analysis of tandem repeats of short sequence motifs (microsatellites) in twelve eukaryotes for which a large part of the genome has been sequenced and assembled. The pattern of motif abundance varies significantly in different species, but it is very similar in different chromosomes of the same species. The most abundant repeats can be classified in two main families. The first family has a rigid conformation, with purines in one strand and pyrimidines in the complementary strand, mainly A(n)/T(n) and (AG)(n)/(CT)(n). The second family has alternating, flexible sequences, such as (AT)(n), (AC)(n) and related sequences. In the pluricellular organisms the relative frequency of both families is rather constant. These observations indicate that microsatellites have structural information and may be involved in the organization of chromatin fibers and in chromosome architecture in general. An additional intriguing finding is the absence of microsatellites with sequences which appear to be forbidden, such as (AATT)(n).

Genetics and smoking behavior

Accumulating data support the role of genetic factors in smoking initiation, progression to tobacco dependence, and smoking persistence. This review summarizes current research on the heritability of tobacco use phenotypes and genetic association studies of smoking-related behaviors. Although progress has been made in genetics research on smoking behavior, many studies have methodological limitations, including insufficient samples for detecting gene-gene and gene-environment interactions and use of less refined phenotypes. Pharmacogenetic investigations also are identifying variants in drug-metabolizing enzymes, receptors, and transporters that modify therapeutic response to smoking cessation medications; however, the field is relatively new, and most findings in this area have yet to be replicated. As this research advances, it will be important to study and address practical, economic, ethical, and social barriers to the translation of genetics research on tobacco use to clinical practice.

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.

MARs and MARBPs

The DNA in eukaryotic genome is compartmentalized into various domains by a series of loops tethered onto the base of nuclear matrix. Scaffold/ Matrix attachment regions (S/MAR) punctuate these attachment sites and govern the nuclear architecture by establishing chromatin boundaries. In this context, specific proteins that interact with and bind to MAR sequences called MAR binding proteins (MARBPs), are of paramount importance, as these sequences spool the proteins that regulate transcription, replication, repair and recombination. Recent evidences also suggest a role for these cis-acting elements in viral integration, replication and transcription, thereby affecting host immune system. Owing to the complex nature of these nucleotide sequences, less is known about the MARBPs that bind to and bring about diverse effects on chromatin architecture and gene function. Several MARBPs have been identified and characterized so far and the list is growing. The fact that most the MARBPs exist in a co-repressor/ co-activator complex and bring about gene regulation makes them quintessential for cellular processes. This participation in gene regulation means that any perturbation in the regulation and levels of MARBPs could lead to disease conditions, particularly those caused by abnormal cell proliferation, like cancer. In the present chapter, we discuss the role of MARs and MARBPs in eukaryotic gene regulation, recombination, transcription and viral integration by altering the local chromatin structure and their dysregulation in disease manifestation

Structure of two genes encoding parallel prothrombin activators in Tropidechis carinatus snake: gene duplication and recruitment of factor X gene to the venom gland

BACKGROUND

Proteins with new function originate through gene duplication followed by divergence. In nature, occurrence of structurally and functionally similar proteins performing highly diverse physiological roles within an organism is rare. Several Australian elapid snakes have two parallel prothrombin activating systems with distinct physiological roles. For example, in Tropidechis carinatus, trocarin D exists in the venom and acts as toxin, whereas coagulation factor (F) X exists in plasma and plays a role in hemostasis.

RESULTS

Here, we show that FX and the trocarin D genes are expressed in a highly tissue-specific manner in T. carinatus. To understand their origin, recruitment and evolution, we determined the complete structure of their genes. Both genes have eight exons with identical exon-intron boundaries. All the introns are 92-99% identical with the exception of intron 1, indicating a recent gene duplication event. The first intron of the trocarin D gene is also nearly identical to that of the FX gene, except for two deletions (255 and 1406 bp) and three insertions (214, 1975, and 2174 bp). The third insertion has a potential scaffold/matrix attached region. The putative promoter of the snake FX gene shares similar cis-elements compared with those of human and murine FX genes. Interestingly, the trocarin D promoter has a 264-bp insertion with core promoter sequences and cis-elements that are known to induce high-level expression. This insertion might be responsible for switching from constitutive expression of the FX gene to inducible expression of trocarin D. Thus, we named this segment as VERSE (Venom Recruitment/Switch Element).

CONCLUSION

To our knowledge, this is the first molecular evidence for the recruitment of a duplicated gene for expression in venom glands by a simple insertion.

Molecular evolution of toxin genes in Elapidae snakes

The venom of the sea krait, Laticauda semifasciata, consists primarily of two toxic proteins, phospholipase A(2) (PLA(2)) and a three-finger-structure toxin. We have cloned both toxic protein genes, including the upstream region. PLA(2) genes contain three types of inserted sequences: an AG-rich region, a chicken repeat 1-like long interspersed nucleotide element sequence and an intron II 3' side repeat sequence. The molecular divergence of L. semifasciata PLA(2) genes was defined on the basis of the inserted sequences and their sequence homology. The length of intron I in the three-finger-structure toxin genes differs from species to species. The alignment analysis of intron I of the three-finger-structure toxin genes revealed that the intron I sequence of the ancestral gene comprised ten genetic regions. A hypothetical evolutionary process for the three-finger-structure toxin genes has also been developed.

Large-scale oscillation of structure-related DNA sequence features in human chromosome 21

Human chromosome 21 is the only chromosome in the human genome that exhibits oscillation of the (G+C) content of a cycle length of hundreds kilobases (kb) ( 500 kb near the right telomere). We aim at establishing the existence of a similar periodicity in structure-related sequence features in order to relate this (G+C)% oscillation to other biological phenomena. The following quantities are shown to oscillate with the same 500 kb periodicity in human chromosome 21: binding energy calculated by two sets of dinucleotide-based thermodynamic parameters, AA/TT and AAA/TTT bi- and tri-nucleotide density, 5'-TA-3' dinucleotide density, and signal for 10- or 11-base periodicity of AA/TT or AAA/TTT. These intrinsic quantities are related to structural features of the double helix of DNA molecules, such as base-pair binding, untwisting or unwinding, stiffness, and a putative tendency for nucleosome formation.

Overview of the structure of all-AT oligonucleotides: organization in helices and packing interactions

We present the crystalline organization of 33 all-AT deoxyoligonucleotide duplexes, studied by x-ray diffraction. Most of them have very similar structures, with Watson-Crick basepairs and a standard average twist close to 36 degrees. The molecules are organized as parallel columns of stacked duplexes in a helical arrangement. Such organization of duplexes is very regular and repetitive: all sequences show the same pattern. It is mainly determined by the stacking of the terminal basepairs, so that the twist in the virtual TA base step between neighbor duplexes is always negative, approximately -22 degrees. The distance between the axes of parallel columns is practically identical in all cases, approximately 26 A. Interestingly, it coincides with that found in DNA viruses and fibers in their hexagonal phase. It appears to be a characteristic distance for ordered parallel DNA molecules. This feature is due to the absence of short range intermolecular forces, which are usually due to the presence of CG basepairs at the end of the oligonucleotide sequence. The duplexes apparently interact only through their diffuse ionic atmospheres. The results obtained can thus be considered as intermediate between liquid crystals, fibers, and standard crystal structures. They provide new information on medium range DNA-DNA interactions.

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.

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.

DNA coiled coils

We report the coiled-coil structure of DNA, which is generated by the dodecanucleotide d(ATATATATATAT). The structure has been determined by single-crystal x-ray crystallography. The molecules form duplexes with single-stranded overhangs, which associate with neighbor molecules and give rise to infinite double helices in a coiled-coil conformation, with staggered nicks in both strands. The coiled coils have six dodecamer duplexes per turn. Despite the presence of nicks, the structure is very rigid. Statistical disorder is present, which gives rise to continuous scattering along the layer lines. This observation can be interpreted by the classical theory of coiled coils developed for proteins, which we apply here to a DNA structure. A clear splitting of the original layer lines of the DNA double helix is detected. The structure we have found adds a previously unrecognized element to the architectures that can be built from DNA oligonucleotides. Any duplex with complementary single-strand overhangs should be expected to give rise to regular coiled-coil structures.

Isolation and characterisation of three moss-derived beta-tubulin promoters suitable for recombinant expression

The moss Physcomitrella patens is an excellent tool to study plant gene-function relationships due to its high rate of homologous recombination (HR). It has also been shown to be very useful in the production of recombinant proteins which are secreted into a simple medium. Thus, there is a need for suitable promoters functional in this well established model organism. We isolated genomic flanking regions of the beta-tubulin gene family from Physcomitrella, concentrating on those family members showing high transcript abundance integrated over gametophytic tissues. Using a novel, fast and reliable quantification assay based on the transient expression and secretion of a recombinant human protein, three genomic upstream regions were characterised in serial deletion constructs. Expression rates were up to three times higher than those obtained with the 35S cauliflower mosaic virus (35S) promoter, which served as a reference.

Reconstruction of putative DNA virus from endogenous rice tungro bacilliform virus-like sequences in the rice genome: implications for integration and evolution

Background

Plant genomes contain various kinds of repetitive sequences such as transposable elements, microsatellites, tandem repeats and virus-like sequences. Most of them, with the exception of virus-like sequences, do not allow us to trace their origins nor to follow the process of their integration into the host genome. Recent discoveries of virus-like sequences in plant genomes led us to set the objective of elucidating the origin of the repetitive sequences. Endogenous rice tungro bacilliform virus (RTBV)-like sequences (ERTBVs) have been found throughout the rice genome. Here, we reconstructed putative virus structures from RTBV-like sequences in the rice genome and characterized to understand evolutionary implication, integration manner and involvements of endogenous virus segments in the corresponding disease response.

Results

We have collected ERTBVs from the rice genomes. They contain rearranged structures and no intact ORFs. The identified ERTBV segments were shown to be phylogenetically divided into three clusters. For each phylogenetic cluster, we were able to make a consensus alignment for a circular virus-like structure carrying two complete ORFs. Comparisons of DNA and amino acid sequences suggested the closely relationship between ERTBV and RTBV. The Oryza AA-genome species vary in the ERTBV copy number. The species carrying low-copy-number of ERTBV segments have been reported to be extremely susceptible to RTBV. The DNA methylation state of the ERTBV sequences was correlated with their copy number in the genome.

Conclusions

These ERTBV segments are unlikely to have functional potential as a virus. However, these sequences facilitate to establish putative virus that provided information underlying virus integration and evolutionary relationship with existing virus. Comparison of ERTBV among the Oryza AA-genome species allowed us to speculate a possible role of endogenous virus segments against its related disease.

Nicotine dependence in a prospective population-based study of adolescents

Dopamine is a key neurotransmitter of the mesolimbic reward pathway in the human brain, and tyrosine hydroxylase (TH) is the rate-limiting enzyme in dopamine biosynthesis. Consequently, the gene encoding TH is a strong candidate for involvement in the genetic component of addiction. The importance of this gene in nicotine dependence is supported by many studies showing a link between nicotine administration and TH expression. A functional tetranucleotide repeat polymorphism within intron 1 of the TH gene (HUMTH01-VNTR) has been shown to modify tobacco use in two independent Caucasian samples from the USA and Australia. Using information drawn from an eight-wave Australian population-based longitudinal study of adolescent health, we tested the effect of the HUMTH01-VNTR on nicotine dependence. Comparisons were made between dependent smokers and non-dependent smokers. These data provide further support for a protective association between the K4 allele and dependent smoking (odds ratio 0.54, 95% confidence interval 0.28-1.0). No associations were observed at any of three other common TH polymorphisms (rs6356, rs6357 and HUMTH01-PstI). Including these data, three independent studies, two of which use identical phenotypes, have now identified a protective relationship between the K4 allele of the functional HUMTH01-VNTR polymorphism and high-level smoking.

Non-random features of loop-size chromatin fragmentation

Upon isolation of DNA from normal eukaryotic cells by standard methods involving extensive proteolytic treatment, a rather homogeneous population of loop-size, double-stranded DNA fragments is regularly obtained. These DNA molecules can be efficiently end-labeled by the DNA polymerase I Klenow fragment, as well as by a 3'- to -5'-exonuclease-free Klenow enzyme, but not by terminal transferase (TdT) unless the ends have been filled up by Klenow, suggesting that dominantly 5' protruding termini are generated upon fragmentation. The filled-up termini were used for cloning the distal parts of the approximately 50 kb fragments. BLAST analysis of the sequence of several clones allowed us to determine the sequence of the non-cloned side of the breakpoints. Comparison of 25, 600 bp-long breakpoint sequences demonstrated prevalence of repetitive elements. Consensus motives characteristic of the breakpoint sequences have been identified. Several sequences exhibit peculiar computed conformational characteristics, with sharp transition or center of symmetry located exactly at the breakpoint. Our data collectively suggest that chromatin fragmentation involves nucleolytic cleavages at fragile/hypersensitive sites delimiting loop-size fragments in a non-random manner. Interestingly, the sequence characteristics of the breakpoints are reminiscent of certain breakpoint cluster regions frequently subject to gene rearrangements.

Molecular evolution and diversification of snake toxin genes, revealed by analysis of intron sequences

The genes encoding erabutoxin (short chain neurotoxin) isoforms (Ea, Eb, and Ec), LsIII (long chain neurotoxin) and a novel long chain neurotoxin pseudogene were cloned from a Laticauda semifasciata genomic library. Short and long chain neurotoxin genes were also cloned from the genome of Laticauda laticaudata, a closely related species of L. semifasciata, by PCR. A putative matrix attached region (MAR) sequence was found in the intron I of the LsIII gene. Comparative analysis of 11 structurally relevant snake toxin genes (three-finger-structure toxins) revealed the molecular evolution of these toxins. Three-finger-structure toxin genes diverged from a common ancestor through two types of evolutionary pathways (long and short types), early in the course of evolution. At a later stage of evolution in each gene, the accumulation of mutations in the exons, especially exon II, by accelerated evolution may have caused the increased diversification in their functions. It was also revealed that the putative MAR sequence found in the LsIII gene was integrated into the gene after the species-level divergence.

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.