Chromatin heterogeneity and distribution of regulatory elements in the late-replicating intercalary heterochromatin domains of drosophila melanogaster chromosomes

19Citations
Citations of this article
16Readers
Mendeley users who have this article in their library.

Abstract

Late-replicating domains (intercalary heterochromatin) in the Drosophila genome display a number of features suggesting their organization is quite unique. Typically, they are quite large and encompass clusters of functionally unrelated tissue-specific genes. They correspond to the topologically associating domains and conserved microsynteny blocks. Our study aims at exploring further details of molecular organization of intercalary heterochromatin and has uncovered surprising heterogeneity of chromatin composition in these regions. Using the 4HMM model developed in our group earlier, intercalary heterochromatin regions were found to host chromatin fragments with a particular epigenetic profile. Aquamarine chromatin fragments (spanning 0.67% of late-replicating regions) are characterized as a class of sequences that appear heterogeneous in terms of their decompactization. These fragments are enriched with enhancer sequences and binding sites for insulator proteins. They likely mark the chromatin state that is related to the binding of cis-regulatory proteins. Malachite chromatin fragments (11% of late-replicating regions) appear to function as universal transitional regions between two contrasting chromatin states. Namely, they invariably delimit intercalary heterochromatin regions from the adjacent active chromatin of interbands. Malachite fragments also flank aquamarine fragments embedded in the repressed chromatin of latereplicating regions. Significant enrichment of insulator proteins CP190, SU(HW), and MOD2.2 was observed in malachite chromatin. Neither aquamarine nor malachite chromatin types appear to correlate with the positions of highly conserved non-coding elements (HCNE) that are typically replete in intercalary heterochromatin. Malachite chromatin found on the flanks of intercalary heterochromatin regions tends to replicate earlier than the malachite chromatin embedded in intercalary heterochromatin. In other words, there exists a gradient of replication progressing from the flanks of intercalary heterochromatin regions center-wise. The peculiar organization and features of replication in large late-replicating regions are discussed as possible factors shaping the evolutionary stability of intercalary heterochromatin.

Figures

  • Fig 1. Two representative chromatin passports for the IH regions 9A3 (A) and 89E1-4 (BXC) (B). 4HMMmap vs compactization map (3CM) by Milon et al., [10], 5 chromatin states by Filion et al., [4], 9 state map by Kharchenko et al., [8] and distribution of enhancer sites, sites of insulator protein binding,
  • Fig 2. Percentages of 4HMM-predicted chromatin types in IH regions (A) and genome-wide (B). The length of each chromatin type (kb) is indicated inside the bars (the length of gaps was not taken into account).
  • Fig 3. Contacts between the 4HMM states. The width of connecting lines is proportional to the frequency of contacts observed.
  • Fig 4. Overlap between 4HMM fragments and chromatin states and types. 5 principal chromatin colors reported in Filion et al., [4] (A); 9 chromatin states by Kharchenko et al., [8], S2 cells (B), BG3 cells (C); and 3 chromatin compactization classes [10] (D).
  • Fig 5. Gene expression in different 4HMM chromatin types. (A) The number of tissues where genes are active (RPKM>3). (B) Magnitude of gene expression summed for 29 tissues. Expression range in S2 (C) and Kc167 cells (D). Quartiles computed for the RPKM values are classified by the chromatin type. The distribution of the first, second and third quartiles of RPKM values for the datasets of transcripts are restricted by the chromatin types. For each chromatin color, the bottom part of the bar denotes the interval from the first to the second quartile; the top part denotes the interval from the second to the third quartile. In all panels, various chromatin types are shown on the X axis (from left to right: border malachite IH, internal malachite IH, internal aquamarine IH, aquamarine genome, ruby genome. Whiskers below/above the 1st/3rd quartile correspond to the 12.5th and 87.5th percentiles.
  • Fig 6. Enhancers and protein distributions across different 4HMM chromatin types. (A) Ratio of the observed fraction of overlapping fragments to the expected one. Observed fraction means the ratio of the total length of genomic regions associated with a protein of interest to the total length of the chromatin type in IH. Expected fraction is the fraction of overlap expected by chance (under random distribution model). Only the values above the “expected” threshold are shown. Asterisks denote probabilities of occurrence by chance *–p<0.05; **–p<1E-3; ***–p<1E-25. (B) Probability values that the observed overlap happened by chance. Bar height (-log10[P]) shows the significance levels for the enrichment of a chromatin type with regulatory elements or proteins indicated on the X axis. The probabilities were computed by the permutation Monte-Carlo test, as described in Materials and Methods section. Enhancers (1) and (2) are taken from RedFly [36] and Kvon et al. [29], respectively.
  • Fig 7. Enrichment levels of the paused RNA pol II in the four chromatin types. X axis shows the number of 5’-forward short non-polyadenylated transcripts (reads) produced by the paused RNA pol II (peaks). The density of peaks overlapping with a particular chromatin type is shown on the Y axis.
  • Fig 8. Replication timing of various 4HMM statesmeasured genome-wide in Kc cells.Replication timing score is shown on the X axis with positive and negative values corresponding to early and late replication timing, respectively. The proportion of probes in each bin is shown on the Y axis.

References Powered by Scopus

BEDTools: A flexible suite of utilities for comparing genomic features

16333Citations
N/AReaders
Get full text

The biology of chromatin remodeling complexes

1788Citations
N/AReaders
Get full text

Three-dimensional folding and functional organization principles of the Drosophila genome

1463Citations
N/AReaders
Get full text

Cited by Powered by Scopus

REDfly: The transcriptional regulatory element database for Drosophila

43Citations
N/AReaders
Get full text

Similarity in replication timing between polytene and diploid cells is associated with the organization of the drosophila genome

17Citations
N/AReaders
Get full text

Arp2/3 and Unc45 maintain heterochromatin stability in Drosophila polytene chromosomes

12Citations
N/AReaders
Get full text

Register to see more suggestions

Mendeley helps you to discover research relevant for your work.

Already have an account?

Cite

CITATION STYLE

APA

Khoroshko, V. A., Levitsky, V. G., Zykova, T. Y., Antonenko, O. V., Belyaeva, E. S., & Zhimulev, I. F. (2016). Chromatin heterogeneity and distribution of regulatory elements in the late-replicating intercalary heterochromatin domains of drosophila melanogaster chromosomes. PLoS ONE, 11(6). https://doi.org/10.1371/journal.pone.0157147

Readers over time

‘16‘17‘18‘19‘20‘21‘2302468

Readers' Seniority

Tooltip

PhD / Post grad / Masters / Doc 6

55%

Researcher 4

36%

Professor / Associate Prof. 1

9%

Readers' Discipline

Tooltip

Agricultural and Biological Sciences 7

50%

Biochemistry, Genetics and Molecular Bi... 5

36%

Environmental Science 1

7%

Neuroscience 1

7%

Article Metrics

Tooltip
Social Media
Shares, Likes & Comments: 1

Save time finding and organizing research with Mendeley

Sign up for free
0