Late replication domains are evolutionary conserved in the Drosophila genome

Dynamic changes in ORC localization and replication fork progression during tissue differentiation

Background

Genomic regions repressed for DNA replication, resulting in either delayed replication in S phase or underreplication in polyploid cells, are thought to be controlled by inhibition of replication origin activation. Studies in Drosophila polytene cells, however, raised the possibility that impeding replication fork progression also plays a major role.

Results

We exploited genomic regions underreplicated (URs) with tissue specificity in Drosophila polytene cells to analyze mechanisms of replication repression. By localizing the Origin Recognition Complex (ORC) in the genome of the larval fat body and comparing this to ORC binding in the salivary gland, we found that sites of ORC binding show extensive tissue specificity. In contrast, there are common domains nearly devoid of ORC in the salivary gland and fat body that also have reduced density of ORC binding sites in diploid cells. Strikingly, domains lacking ORC can still be replicated in some polytene tissues, showing absence of ORC and origins is insufficient to repress replication. Analysis of the width and location of the URs with respect to ORC position indicates that whether or not a genomic region lacking ORC is replicated is controlled by whether replication forks formed outside the region are inhibited.

Conclusions

These studies demonstrate that inhibition of replication fork progression can block replication across genomic regions that constitutively lack ORC. Replication fork progression can be inhibited in both tissue-specific and genome region-specific ways. Consequently, when evaluating sources of genome instability it is important to consider altered control of replication forks in response to differentiation.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4992-3) contains supplementary material, which is available to authorized users.

Banding Pattern of Polytene Chromosomes as a Representation of Universal Principles of Chromatin Organization into Topological Domains

Drosophila polytene chromosomes are widely used as a model of eukaryotic interphase chromosomes. The most noticeable feature of polytene chromosome is transverse banding associated with alternation of dense stripes (dark or black bands) and light diffuse areas that encompass alternating less compact gray bands and interbands visible with an electron microscope. In recent years, several approaches have been developed to predict location of morphological structures of polytene chromosomes based on the distribution of proteins on the molecular map of Drosophila genome. Comparison of these structures with the results of analysis of the three-dimensional chromatin organization by the Hi-C method indicates that the morphology of polytene chromosomes represents direct visualization of the interphase nucleus spatial organization into topological domains. Compact black bands correspond to the extended topological domains of inactive chromatin, while interbands are the barriers between the adjacent domains. Here, we discuss the prospects of using polytene chromosomes to study mechanisms of spatial organization of interphase chromosomes, as well as their dynamics and evolution.

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

Morphologically, polytene chromosomes of Drosophila melanogaster consist of compact “black” bands alternating with less compact “grey” bands and interbands. We developed a comprehensive approach that combines cytological mapping data of FlyBase-annotated genes and novel tools for predicting cytogenetic features of chromosomes on the basis of their protein composition and determined the genomic coordinates for all black bands of polytene chromosome 2R. By a PCNA immunostaining assay, we obtained the replication timetable for all the bands mapped. The results allowed us to compare replication timing between polytene chromosomes in salivary glands and chromosomes from cultured diploid cell lines and to observe a substantial similarity in the global replication patterns at the band resolution level. In both kinds of chromosomes, the intervals between black bands correspond to early replication initiation zones. Black bands are depleted of replication initiation events and are characterized by a gradient of replication timing; therefore, the time of replication completion correlates with the band length. The bands are characterized by low gene density, contain predominantly tissue-specific genes, and are represented by silent chromatin types in various tissues. The borders of black bands correspond well to the borders of topological domains as well as to the borders of the zones showing H3K27me3, SUUR, and LAMIN enrichment. In conclusion, the characteristic pattern of polytene chromosomes reflects partitioning of the Drosophila genome into two global types of domains with contrasting properties. This partitioning is conserved in different tissues and determines replication timing in Drosophila.

Chromatin Heterogeneity and Distribution of Regulatory Elements in the Late-Replicating Intercalary Heterochromatin Domains of Drosophila melanogaster Chromosomes

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 late-replicating 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.

Genetic organization of interphase chromosome bands and interbands in Drosophila melanogaster

Drosophila melanogaster polytene chromosomes display specific banding pattern; the underlying genetic organization of this pattern has remained elusive for many years. In the present paper, we analyze 32 cytology-mapped polytene chromosome interbands. We estimated molecular locations of these interbands, described their molecular and genetic organization and demonstrate that polytene chromosome interbands contain the 5' ends of housekeeping genes. As a rule, interbands display preferential "head-to-head" orientation of genes. They are enriched for "broad" class promoters characteristic of housekeeping genes and associate with open chromatin proteins and Origin Recognition Complex (ORC) components. In two regions, 10A and 100B, coding sequences of genes whose 5'-ends reside in interbands map to constantly loosely compacted, early-replicating, so-called "grey" bands. Comparison of expression patterns of genes mapping to late-replicating dense bands vs genes whose promoter regions map to interbands shows that the former are generally tissue-specific, whereas the latter are represented by ubiquitously active genes. Analysis of RNA-seq data (modENCODE-FlyBase) indicates that transcripts from interband-mapping genes are present in most tissues and cell lines studied, across most developmental stages and upon various treatment conditions. We developed a special algorithm to computationally process protein localization data generated by the modENCODE project and show that Drosophila genome has about 5700 sites that demonstrate all the features shared by the interbands cytologically mapped to date.

Underreplicated regions in Drosophila melanogaster are enriched with fast-evolving genes and highly conserved noncoding sequences

Many late replicating regions are underreplicated in polytene chromosomes of Drosophila melanogaster. These regions contain silenced chromatin and overlap long syntenic blocks of conserved gene order in drosophilids. In this report we show that in D. melanogaster the underreplicated regions are enriched with fast-evolving genes lacking homologs in distant species such as mosquito or human, indicating that the phylogenetic conservation of genes correlates with replication timing and chromatin status. Drosophila genes without human homologs located in the underreplicated regions have higher nonsynonymous substitution rate and tend to encode shorter proteins when compared with those in the adjacent regions. At the same time, the underreplicated regions are enriched with ultraconserved elements and highly conserved noncoding sequences, especially in introns of very long genes indicating the presence of an extensive regulatory network that may be responsible for the conservation of gene order in these regions. The regions have a modest preference for long noncoding RNAs but are depleted for small nucleolar RNAs, microRNAs, and transfer RNAs. Our results demonstrate that the underreplicated regions have a specific genic composition and distinct pattern of evolution.