Gene density profile reveals the marking of late replicated domains in the Drosophila melanogaster genome

Replication timing in Drosophila and its peculiarities in polytene chromosomes

Drosophila melanogaster is one of the popular model organisms in DNA replication studies. Since the 1960s, DNA replication of polytene chromosomes has been extensively studied by cytological methods. In the recent two decades, the progress in our understanding of DNA replication was associated with new techniques. Use of fluorescent dyes increased the resolution of cytological methods significantly. High-throughput methods allowed analysis of DNA replication on a genome scale, as well as its correlation with chromatin structure and gene activi ty. Precise mapping of the cytological structures of polytene chromosomes to the genome assembly allowed comparison of replication between polytene chromosomes and chromosomes of diploid cells. New features of replication characteristic for D. melanogaster were described for both diploid and polytene chromosomes. Comparison of genomic replication profiles revealed a significant similarity between Drosophila and other well-studi ed eukaryotic species, such as human. Early replication is often confined to intensely transcribed gene-dense regions characterized by multiple replication initiation sites. Features of DNA replication in Drosophila might be explained by a compact genome. The organization of replication in polytene chromosomes has much in common with the organization of replication in chromosomes in diploid cells. The most important feature of replication in polytene chromosomes is its low rate and the dependence of S-phase duration on many factors: external and internal, local and global. The speed of replication forks in D. melanogaster polytene chromosomes is affected by SUUR and Rif1 proteins. It is not known yet how universal the mechanisms associated with these factors are, but their study is very promising.

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.

Functional dissection of Drosophila melanogaster SUUR protein influence on H3K27me3 profile

BACKGROUND

In eukaryotes, heterochromatin replicates late in S phase of the cell cycle and contains specific covalent modifications of histones. SuUR mutation found in Drosophila makes heterochromatin replicate earlier than in wild type and reduces the level of repressive histone modifications. SUUR protein was shown to be associated with moving replication forks, apparently through the interaction with PCNA. The biological process underlying the effects of SUUR on replication and composition of heterochromatin remains unknown.

RESULTS

Here we performed a functional dissection of SUUR protein effects on H3K27me3 level. Using hidden Markow model-based algorithm we revealed SuUR-sensitive chromosomal regions that demonstrated unusual characteristics: They do not contain Polycomb and require SUUR function to sustain H3K27me3 level. We tested the role of SUUR protein in the mechanisms that could affect H3K27me3 histone levels in these regions. We found that SUUR does not affect the initial H3K27me3 pattern formation in embryogenesis or Polycomb distribution in the chromosomes. We also ruled out the possible effect of SUUR on histone genes expression and its involvement in DSB repair.

CONCLUSIONS

Obtained results support the idea that SUUR protein contributes to the heterochromatin maintenance during the chromosome replication. A model that explains major SUUR-associated phenotypes is proposed.

Computer Analysis of Glioma Transcriptome Profiling: Alternative Splicing Events

Here we present the analysis of alternative splicing events on an example of glioblastoma cell culture samples using a set of computer tools in combination with database integration. The gene expression profiles of glioblastoma were obtained from cell culture samples of primary glioblastoma which were isolated and processed for RNA extraction. Transcriptome profiling of normal brain samples and glioblastoma were done by Illumina sequencing. The significant differentially expressed exon-level probes and their corresponding genes were identified using a combination of the splicing index method. Previous studies indicated that tumor-specific alternative splicing is important in the regulation of gene expression and corresponding protein functions during cancer development. Multiple alternative splicing transcripts have been identified as progression markers, including generalized splicing abnormalities and tumor- and stage-specific events. We used a set of computer tools which were recently applied to analysis of gene expression in laboratory animals to study differential splicing events. We found 69 transcripts that are differentially alternatively spliced. Three cancer-associated genes were considered in detail, in particular: APP (amyloid beta precursor protein), CASC4 (cancer susceptibility candidate 4) and TP53. Such alternative splicing opens new perspectives for cancer research.

Protein and Genetic Composition of Four Chromatin Types in Drosophila melanogaster Cell Lines

BACKGROUND

Recently, we analyzed genome-wide protein binding data for the Drosophila cell lines S2, Kc, BG3 and Cl.8 (modENCODE Consortium) and identified a set of 12 proteins enriched in the regions corresponding to interbands of salivary gland polytene chromosomes. Using these data, we developed a bioinformatic pipeline that partitioned the Drosophila genome into four chromatin types that we hereby refer to as aquamarine, lazurite, malachite and ruby.

RESULTS

Here, we describe the properties of these chromatin types across different cell lines. We show that aquamarine chromatin tends to harbor transcription start sites (TSSs) and 5' untranslated regions (5'UTRs) of the genes, is enriched in diverse "open" chromatin proteins, histone modifications, nucleosome remodeling complexes and transcription factors. It encompasses most of the tRNA genes and shows enrichment for non-coding RNAs and miRNA genes. Lazurite chromatin typically encompasses gene bodies. It is rich in proteins involved in transcription elongation. Frequency of both point mutations and natural deletion breakpoints is elevated within lazurite chromatin. Malachite chromatin shows higher frequency of insertions of natural transposons. Finally, ruby chromatin is enriched for proteins and histone modifications typical for the "closed" chromatin. Ruby chromatin has a relatively low frequency of point mutations and is essentially devoid of miRNA and tRNA genes. Aquamarine and ruby chromatin types are highly stable across cell lines and have contrasting properties. Lazurite and malachite chromatin types also display characteristic protein composition, as well as enrichment for specific genomic features. We found that two types of chromatin, aquamarine and ruby, retain their complementary protein patterns in four Drosophila cell lines.

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.

DNA replication and transcription programs respond to the same chromatin cues

DNA replication is a dynamic process that occurs in a temporal order along each of the chromosomes. A consequence of the temporally coordinated activation of replication origins is the establishment of broad domains (>100 kb) that replicate either early or late in S phase. This partitioning of the genome into early and late replication domains is important for maintaining genome stability, gene dosage, and epigenetic inheritance; however, the molecular mechanisms that define and establish these domains are poorly understood. The modENCODE Project provided an opportunity to investigate the chromatin features that define the Drosophila replication timing program in multiple cell lines. The majority of early and late replicating domains in the Drosophila genome were static across all cell lines; however, a small subset of domains was dynamic and exhibited differences in replication timing between the cell lines. Both origin selection and activation contribute to defining the DNA replication program. Our results suggest that static early and late replicating domains were defined at the level of origin selection (ORC binding) and likely mediated by chromatin accessibility. In contrast, dynamic domains exhibited low ORC densities in both cell types, suggesting that origin activation and not origin selection governs the plasticity of the DNA replication program. Finally, we show that the male-specific early replication of the X chromosome is dependent on the dosage compensation complex (DCC), suggesting that the transcription and replication programs respond to the same chromatin cues. Specifically, MOF-mediated hyperacetylation of H4K16 on the X chromosome promotes both the up-regulation of male-specific transcription and origin activation.

DNA replication in nurse cell polytene chromosomes of Drosophila melanogaster otu mutants

Drosophila cell lines are used extensively to study replication timing, yet data about DNA replication in larval and adult tissues are extremely limited. To address this gap, we traced DNA replication in polytene chromosomes from nurse cells of Drosophila melanogaster otu mutants using bromodeoxyuridine incorporation. Importantly, nurse cells are of female germline origin, unlike the classical larval salivary glands, that are somatic. In contrast to salivary gland polytene chromosomes, where replication begins simultaneously across all puffs and interbands, replication in nurse cells is first observed at several specific chromosomal regions. For instance, in the chromosome 2L, these include the regions 31B-E and 37E and proximal parts of 34B and 35B, with the rest of the decondensed chromosomal regions joining replication process a little later. We observed that replication timing of pericentric heterochromatin in nurse cells was shifted from late S phase to early and mid stages. Curiously, chromosome 4 may represent a special domain of the genome, as it replicates on its own schedule which is uncoupled from the rest of the chromosomes. Finally, we report that SUUR protein, an established marker of late replication in salivary gland polytene chromosomes, does not always colocalize with late-replicating regions in nurse cells.

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.

Late replication domains are evolutionary conserved in the Drosophila genome

Drosophila chromosomes are organized into distinct domains differing in their predominant chromatin composition, replication timing and evolutionary conservation. We show on a genome-wide level that genes whose order has remained unaltered across 9 Drosophila species display late replication timing and frequently map to the regions of repressive chromatin. This observation is consistent with the existence of extensive domains of repressive chromatin that replicate extremely late and have conserved gene order in the Drosophila genome. We suggest that such repressive chromatin domains correspond to a handful of regions that complete replication at the very end of S phase. We further demonstrate that the order of genes in these regions is rarely altered in evolution. Substantial proportion of such regions significantly coincide with large synteny blocks. This indicates that there are evolutionary mechanisms maintaining the integrity of these late-replicating chromatin domains. The synteny blocks corresponding to the extremely late-replicating regions in the D. melanogaster genome consistently display two-fold lower gene density across different Drosophila species.

Comparison of ultra-conserved elements in drosophilids and vertebrates

Metazoan genomes contain many ultra-conserved elements (UCEs), long sequences identical between distant species. In this study we identified UCEs in drosophilid and vertebrate species with a similar level of phylogenetic divergence measured at protein-coding regions, and demonstrated that both the length and number of UCEs are larger in vertebrates. The proportion of non-exonic UCEs declines in distant drosophilids whilst an opposite trend was observed in vertebrates. We generated a set of 2,126 Sophophora UCEs by merging elements identified in several drosophila species and compared these to the eutherian UCEs identified in placental mammals. In contrast to vertebrates, the Sophophora UCEs are depleted around transcription start sites. Analysis of 52,954 P-element, piggyBac and Minos insertions in the D. melanogaster genome revealed depletion of the P-element and piggyBac insertions in and around the Sophophora UCEs. We examined eleven fly strains with transposon insertions into the intergenic UCEs and identified associated phenotypes in five strains. Four insertions behave as recessive lethals, and in one case we observed a suppression of the marker gene within the transgene, presumably by silenced chromatin around the integration site. To confirm the lethality is caused by integration of transposons we performed a phenotype rescue experiment for two stocks and demonstrated that the excision of the transposons from the intergenic UCEs restores viability. Sequencing of DNA after the transposon excision in one fly strain with the restored viability revealed a 47 bp insertion at the original transposon integration site suggesting that the nature of the mutation is important for the appearance of the phenotype. Our results suggest that the UCEs in flies and vertebrates have both common and distinct features, and demonstrate that a significant proportion of intergenic drosophila UCEs are sensitive to disruption.

[Regulation of DNA replication timing]

Although distinct chromatin types have been long known to replicate at different timepoints of S phase, fine replication control has only recently become considered as an epigenetic phenomenon. It is now clear that in course of differentiation significant changes in genome replication timing occur, and these changes are intimately linked with the changes in transcriptional activity and nuclear architecture. Temporally coordinate replication is organized spatially into discrete units having specific chromosomal organization and function. Even though the functional aspects of such tight control of replication timing remain to be explored, one can confidently consider the replication program as yet another fundamental feature characteristic of the given differentiation state. The present review touches upon the molecular mechanisms of spatial and temporal control of replication timing, involving individual replication origins as well as large chromatin domains.

Induced transcription results in local changes in chromatin structure, replication timing, and DNA polytenization in a site of intercalary heterochromatin

In salivary gland polytene chromosomes of Drosophila melanogaster, the regions of intercalary heterochromatin are characterized by late replication, under-replication, and genetic silencing. Using Gal4/UAS system, we induced transcription of sequences adjacent to transgene insertions in the band 11A6-9. This activation resulted in a loss of "silent" and appearance of "active" epigenetic marks, recruitment of RNA polymerase II, and formation of a puff. The activated region is now early replicating and shows increased level of DNA polytenization. Notably, all these changes are restricted to the area around the inserts, whereas the rest of the band remains inactive and late replicating. Although only a short area near the insertion site is transcribed, it results in an "open" chromatin conformation in a much broader region. We conclude that regions of intercalary heterochromatin do not form stand-alone units of late replication and under-replication. Every part of such regions can be activated and polytenized independently of other parts.

Banding patterns in Drosophila melanogaster polytene chromosomes correlate with DNA-binding protein occupancy

The most enigmatic feature of polytene chromosomes is their banding pattern, the genetic organization of which has been a very attractive puzzle for many years. Recent genome-wide protein mapping efforts have produced a wealth of data for the chromosome proteins of Drosophila cells. Based on their specific protein composition, the chromosomes comprise two types of bands, as well as interbands. These differ in terms of time of replication and specific types of proteins. The interbands are characterized by their association with "active" chromatin proteins, nucleosome remodeling, and origin recognition complexes, and so they have three functions: acting as binding sites for RNA pol II, initiation of replication and nucleosome remodeling of short fragments of DNA. The borders and organization of the same band and interband regions are largely identical, irrespective of the cell type studied. This demonstrates that the banding pattern is a universal principle of the organization of interphase polytene and non-polytene chromosomes.

Late replication domains in polytene and non-polytene cells of Drosophila melanogaster

In D. melanogaster polytene chromosomes, intercalary heterochromatin (IH) appears as large dense bands scattered in euchromatin and comprises clusters of repressed genes. IH displays distinctly low gene density, indicative of their particular regulation. Genes embedded in IH replicate late in the S phase and become underreplicated. We asked whether localization and organization of these late-replicating domains is conserved in a distinct cell type. Using published comprehensive genome-wide chromatin annotation datasets (modENCODE and others), we compared IH organization in salivary gland cells and in a Kc cell line. We first established the borders of 60 IH regions on a molecular map, these regions containing underreplicated material and encompassing ∼12% of Drosophila genome. We showed that in Kc cells repressed chromatin constituted 97% of the sequences that corresponded to IH bands. This chromatin is depleted for ORC-2 binding and largely replicates late. Differences in replication timing between the cell types analyzed are local and affect only sub-regions but never whole IH bands. As a rule such differentially replicating sub-regions display open chromatin organization, which apparently results from cell-type specific gene expression of underlying genes. We conclude that repressed chromatin organization of IH is generally conserved in polytene and non-polytene cells. Yet, IH domains do not function as transcription- and replication-regulatory units, because differences in transcription and replication between cell types are not domain-wide, rather they are restricted to small “islands” embedded in these domains. IH regions can thus be defined as a special class of domains with low gene density, which have narrow temporal expression patterns, and so displaying relatively conserved organization.

Late-replicating domains have higher divergence and diversity in Drosophila melanogaster

Several reports from mammals indicate that an increase in the mutation rate in late-replicating regions may, in part, be responsible for the observed genomic heterogeneity in neutral substitution rates and levels of diversity, although the mechanisms for this remain poorly understood. Recent evidence also suggests that late replication is associated with high mutability in yeast. This then raises the question as to whether a similar effect is operating across all eukaryotes. Limited evidence from one chromosome arm in Drosophila melanogaster suggests the opposite pattern, with regions overlapping early-firing origins showing increased levels of diversity and divergence. Given the availability of genome-wide replication timing profiles for D. melanogaster, we now return to this issue. Consistent with what is seen in other taxa, we find that divergence at synonymous sites in exon cores, as well as divergence at putatively unconstrained intronic sites, is elevated in late-replicating regions. Analysis of genes with low codon usage bias suggests a ∼30% difference in mutation rate between the earliest and the latest replicating sequence. Intronic sequence suggests a more modest difference. We additionally show that an increase in diversity in late-replicating sequences is not owing to replication timing covarying with the local recombination rate. If anything, the effects of recombination mask the impact of replication timing. We conclude that, contrary to prior reports and consistent with what is seen in mammals and yeast, there is indeed a relationship between rates of nucleotide divergence and diversity and replication timing that is consistent with an increase in the mutation rate during late S-phase in D. melanogaster. It is therefore plausible that such an effect might be common among eukaryotes. The result may have implications for the inference of positive selection.

Identical functional organization of nonpolytene and polytene chromosomes in Drosophila melanogaster

Salivary gland polytene chromosomes demonstrate banding pattern, genetic meaning of which is an enigma for decades. Till now it is not known how to mark the band/interband borders on physical map of DNA and structures of polytene chromosomes are not characterized in molecular and genetic terms. It is not known either similar banding pattern exists in chromosomes of regular diploid mitotically dividing nonpolytene cells. Using the newly developed approach permitting to identify the interband material and localization data of interband-specific proteins from modENCODE and other genome-wide projects, we identify physical limits of bands and interbands in small cytological region 9F13-10B3 of the X chromosome in D. melanogaster, as well as characterize their general molecular features. Our results suggests that the polytene and interphase cell line chromosomes have practically the same patterns of bands and interbands reflecting, probably, the basic principle of interphase chromosome organization. Two types of bands have been described in chromosomes, early and late-replicating, which differ in many aspects of their protein and genetic content. As appeared, origin recognition complexes are located almost totally in the interbands of chromosomes.