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

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.

Chromatin attachment to the nuclear matrix represses hypocotyl elongation in Arabidopsis thaliana

The nuclear matrix is a nuclear compartment that has diverse functions in chromatin regulation and transcription. However, how this structure influences epigenetic modifications and gene expression in plants is largely unknown. In this study, we show that a nuclear matrix binding protein, AHL22, together with the two transcriptional repressors FRS7 and FRS12, regulates hypocotyl elongation by suppressing the expression of a group of genes known as SMALL AUXIN UP RNAs (SAURs) in Arabidopsis thaliana. The transcriptional repression of SAURs depends on their attachment to the nuclear matrix. The AHL22 complex not only brings these SAURs, which contain matrix attachment regions (MARs), to the nuclear matrix, but it also recruits the histone deacetylase HDA15 to the SAUR loci. This leads to the removal of H3 acetylation at the SAUR loci and the suppression of hypocotyl elongation. Taken together, our results indicate that MAR-binding proteins act as a hub for chromatin and epigenetic regulators. Moreover, we present a mechanism by which nuclear matrix attachment to chromatin regulates histone modifications, transcription, and hypocotyl elongation.

Arabidopsis DNA Replication Initiates in Intergenic, AT-Rich Open Chromatin

The selection and firing of DNA replication origins play key roles in ensuring that eukaryotes accurately replicate their genomes. This process is not well documented in plants due in large measure to difficulties in working with plant systems. We developed a new functional assay to label and map very early replicating loci that must, by definition, include at least a subset of replication origins. Arabidopsis (Arabidopsis thaliana) cells were briefly labeled with 5-ethynyl-2'-deoxy-uridine, and nuclei were subjected to two-parameter flow sorting. We identified more than 5500 loci as initiation regions (IRs), the first regions to replicate in very early S phase. These were classified as strong or weak IRs based on the strength of their replication signals. Strong initiation regions were evenly spaced along chromosomal arms and depleted in centromeres, while weak initiation regions were enriched in centromeric regions. IRs are AT-rich sequences flanked by more GC-rich regions and located predominantly in intergenic regions. Nuclease sensitivity assays indicated that IRs are associated with accessible chromatin. Based on these observations, initiation of plant DNA replication shows some similarity to, but is also distinct from, initiation in other well-studied eukaryotic systems.

Genomic structure and transcript analysis of the Rapid Alkalinization Factor (RALF) gene family during host-pathogen crosstalk in Fragaria vesca and Fragaria x ananassa strawberry

Rapid Alkalinization Factors (RALFs) are cysteine-rich peptides ubiquitous within plant kingdom. They play multiple roles as hormonal signals in diverse processes, including root elongation, cell growth, pollen tube development, and fertilization. Their involvement in host-pathogen crosstalk as negative regulators of immunity in Arabidopsis has also been recognized. In addition, peptides homologous to RALF are secreted by different fungal pathogens as effectors during early stages of infection. Previous studies have identified nine RALF genes in the diploid strawberry (Fragaria vesca) genome. This work describes the genomic organization of the RALF gene families in commercial octoploid strawberry (Fragaria × ananassa) and the re-annotated genome of F. vesca, and then compares findings with orthologs in Arabidopsis thaliana. We reveal the presence of 15 RALF genes in F. vesca genotype Hawaii 4 and 50 in Fragaria x ananassa cv. Camarosa, showing a non-homogenous localization of genes among the different Fragaria x ananassa subgenomes. Expression analysis of Fragaria x ananassa RALF genes upon infection with Colletotrichum acutatum or Botrytis cinerea showed that FanRALF3-1 was the only fruit RALF gene upregulated after fungal infection. In silico analysis was used to identify distinct pathogen inducible elements upstream of the FanRALF3-1 gene. Agroinfiltration of strawberry fruit with deletion constructs of the FanRALF3-1 promoter identified a 5' region required for FanRALF3-1 expression in fruit, but failed to identify a region responsible for fungal induced expression.

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.

Benefits of using genomic insulators flanking transgenes to increase expression and avoid positional effects

For more than 20 years, plant biologists have tried to achieve complete control of transgene expression. Until the techniques to target transgenes to safe harbor sites in the genome become routine, flanking transgenes with genetic insulators, DNA sequences that create independent domains of gene expression, can help avoid positional effects and stabilize their expression. We have, for the first time, compared the effect of three insulator sequences previously described in the literature and one never tested before. Our results indicate that their use increases transgene expression, but only the last one reduces variability between lines and between individuals. We have analyzed the integration of insulator-flanked T-DNAs using whole genome re-sequencing (to our knowledge, also for the first time) and found data suggesting that chiMARs can shelter transgene insertions from neighboring repressive epigenetic states. Finally, we could also observe a loss of accuracy of the RB insertion in the lines harboring insulators, evidenced by a high frequency of truncation of T-DNAs and of insertion of vector backbone that, however, did not affect transgene expression. Our data supports that the effect of each genetic insulator is different and their use in transgenic constructs should depend on the needs of each specific experiment.

Entering the Next Dimension: Plant Genomes in 3D

After linear sequences of genomes and epigenomic landscape data, the 3D organization of chromatin in the nucleus is the next level to be explored. Different organisms present a general hierarchical organization, with chromosome territories at the top. Chromatin interaction maps, obtained by chromosome conformation capture (3C)-based methodologies, for eight plant species reveal commonalities, but also differences, among them and with animals. The smallest structures, found in high-resolution maps of the Arabidopsis genome, are single genes. Epigenetic marks (histone modification and DNA methylation), transcriptional activity, and chromatin interaction appear to be correlated, and whether structure is the cause or consequence of the function of interacting regions is being actively investigated.

A chimeric HS4 insulator-scaffold attachment region enhances transgene expression in transfected Chinese hamster ovary cells

Chinese hamster ovary (CHO) cells are one of the most commonly used expression systems for the production of recombinant proteins but low levels of transgene expression and transgene silencing are frequently encountered. Epigenetic regulatory elements such as the chicken β-globin locus control region hypersensitive site 4 (HS4) and scaffold/matrix attachment regions (S/MARs) have positive effects on transgene expression. In this study, a chimeric HS4-SAR was cloned upstream or downstream of an enhanced green fluorescent protein (eGFP) expression cassette in a eukaryotic vector, and the resulting vectors were transfected into CHO cells. eGFP was detected by flow cytometry. Real-time quantitative PCR (qPCR) was used to determine copy numbers of the stably transfected cells. And fluorescence in situ hybridization (FISH) was used to detect the status of vector in the host cell chromosome. The results showed that HS4-SAR positioned downstream of the expression cassette could enhance eGFP expression by 4.83-fold compared with the control vector. There may not be a relationship between transgene copy number and gene expression level. HS4-SAR did not appear to alter the integration of the transgene into the host cell chromosome or its position in the chromosome. We found a synthetic chimeric HS4-SAR positively increased transgene expression in CHO cells.

Nonrandom domain organization of the Arabidopsis genome at the nuclear periphery

The nuclear space is not a homogeneous biochemical environment. Many studies have demonstrated that the transcriptional activity of a gene is linked to its positioning within the nuclear space. Following the discovery of lamin-associated domains (LADs), which are transcriptionally repressed chromatin regions, the nonrandom positioning of chromatin at the nuclear periphery and its biological relevance have been studied extensively in animals. However, it remains unknown whether comparable chromatin organizations exist in plants. Here, using a strategy using restriction enzyme-mediated chromatin immunoprecipitation, we present genome-wide identification of nonrandom domain organization of chromatin at the peripheral zone of Arabidopsis thaliana nuclei. We show that in various tissues, 10%-20% of the regions on the chromosome arms are anchored at the nuclear periphery, and these regions largely overlap between different tissues. Unlike LADs in animals, the identified domains in plants are not gene-poor or A/T-rich. These domains are enriched with silenced protein-coding genes, transposable element genes, and heterochromatic marks, which collectively define a repressed environment. In addition, these domains strongly correlate with our genome-wide chromatin interaction data set (Hi-C) by largely explaining the patterns of chromatin compartments, revealed on Hi-C maps. Moreover, our results reveal a spatial compartment of different DNA methylation pathways that regulate silencing of transposable elements, where the CHH methylation of transposable elements located at the nuclear periphery and in the interior are preferentially mediated by CMT2 and DRM methyltransferases, respectively. Taken together, the results demonstrate functional partitioning of the Arabidopsis genome in the nuclear space.

Recurrence of Chromosome Rearrangements and Reuse of DNA Breakpoints in the Evolution of the Triticeae Genomes

Chromosomal rearrangements (CRs) play important roles in karyotype diversity and speciation. While many CR breakpoints have been characterized at the sequence level in yeast, insects, and primates, little is known about the structure of evolutionary CR breakpoints in plant genomes, which are much more dynamic in genome size and sequence organization. Here, we report identification of breakpoints of a translocation between chromosome arms 4L and 5L of Triticeae, which is fixed in several species, including diploid wheat and rye, by comparative mapping and analysis of the draft genome and chromosome survey sequences of the Triticeae species. The wheat translocation joined the ends of breakpoints downstream of a WD40 gene on 4AL and a gene of the PMEI family on 5AL. A basic helix-loop-helix transcription factor gene in 5AL junction was significantly restructured. Rye and wheat share the same position for the 4L breakpoint, but the 5L breakpoint positions are not identical, although very close in these two species, indicating the recurrence of 4L/5L translocations in the Triticeae. Although barley does not carry the translocation, collinearity across the breakpoints was violated by putative inversions and/or transpositions. Alignment with model grass genomes indicated that the translocation breakpoints coincided with ancient inversion junctions in the Triticeae ancestor. Our results show that the 4L/5L translocation breakpoints represent two CR hotspots reused during Triticeae evolution, and support breakpoint reuse as a widespread mechanism in all eukaryotes. The mechanisms of the recurrent translocation and its role in Triticeae evolution are also discussed.

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.

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.

Noise-plasticity correlations of gene expression in the multicellular organism Arabidopsis thaliana

Gene expression levels exhibit stochastic variations among genetically identical organisms under the same environmental conditions (called gene expression "noise" or phenotype "fluctuation"). In yeast and Escherichia coli, positive correlations have been found between such gene expression noise and "plasticity" with environmental variations. To determine the universality of such correlations in both unicellular and multicellular organisms, we focused on the relationships between gene expression "noise" and "plasticity" in Arabidopsis thaliana, a multicellular model organism. In recent studies on yeast and E. coli, only some gene groups with specific properties of promoter architecture, average expression levels, and functions exhibited strong noise-plasticity correlations. However, we found strong noise-plasticity correlations for most gene groups in Arabidopsis; additionally, promoter architecture, functional essentiality of genes, and circadian rhythm appeared to have only a weak influence on the correlation strength. The differences in the characteristics of noise-plasticity correlations may result from three-dimensional chromosomal structures and/or circadian rhythm.

Sequence structure of Lowary/Widom clones forming strong nucleosomes

Lowary and Widom selected from random sequences those which form exceptionally stable nucleosomes, including clone 601, the current champion of strong nucleosome (SN) sequences. This unique sequence database (LW sequences) carries sequence elements which confer stability on the nucleosomes formed on the sequences, and, thus, may serve as source of information on the structure of "ideal" or close to ideal nucleosome DNA sequence. An important clue is also provided by crystallographic study of Vasudevan and coauthors on clone 601 nucleosomes. It demonstrated that YR·YR dinucleotide stacks (primarily TA·TA) follow one another at distances 10 or 11 bases or multiples thereof, such that they all are located on the interface between DNA and histone octamer. Combining this important information with alignment of the YR-containing 10-mers and 11-mers from LW sequences, the bendability matrices of the stable nucleosome DNA are derived. The matrices suggest that the periodically repeated TA (YR), RR, and YY dinucleotides are the main sequence features of the SNs. This consensus coincides with the one for recently discovered SNs with visibly periodic DNA sequences. Thus, the experimentally observed stable LW nucleosomes and SNs derived computationally appear to represent the same entity - exceptionally stable SNs.

Insights into the evolution and diversification of the AT-hook Motif Nuclear Localized gene family in land plants

BACKGROUND

Members of the ancient land-plant-specific transcription factor AT-Hook Motif Nuclear Localized (AHL) gene family regulate various biological processes. However, the relationships among the AHL genes, as well as their evolutionary history, still remain unexplored.

RESULTS

We analyzed over 500 AHL genes from 19 land plant species, ranging from the early diverging Physcomitrella patens and Selaginella to a variety of monocot and dicot flowering plants. We classified the AHL proteins into three types (Type-I/-II/-III) based on the number and composition of their functional domains, the AT-hook motif(s) and PPC domain. We further inferred their phylogenies via Bayesian inference analysis and predicted gene gain/loss events throughout their diversification. Our analyses suggested that the AHL gene family emerged in embryophytes and further evolved into two distinct clades, with Type-I AHLs forming one clade (Clade-A), and the other two types together diversifying in another (Clade-B). The two AHL clades likely diverged before the separation of Physcomitrella patens from the vascular plant lineage. In angiosperms, Clade-A AHLs expanded into 5 subfamilies; while, the ones in Clade-B expanded into 4 subfamilies. Examination of their expression patterns suggests that the AHLs within each clade share similar expression patterns with each other; however, AHLs in one monophyletic clade exhibit distinct expression patterns from the ones in the other clade. Over-expression of a Glycine max AHL PPC domain in Arabidopsis thaliana recapitulates the phenotype observed when over-expressing its Arabidopsis thaliana counterpart. This result suggests that the AHL genes from different land plant species may share conserved functions in regulating plant growth and development. Our study further suggests that such functional conservation may be due to conserved physical interactions among the PPC domains of AHL proteins.

CONCLUSIONS

Our analyses reveal a possible evolutionary scenario for the AHL gene family in land plants, which will facilitate the design of new studies probing their biological functions. Manipulating the AHL genes has been suggested to have tremendous effects in agriculture through increased seedling establishment, enhanced plant biomass and improved plant immunity. The information gleaned from this study, in turn, has the potential to be utilized to further improve crop production.