Background: Previously, we developed a simple method for carrying out a restriction enzyme analysis of eukaryotic DNA in silico, based on the known DNA sequences of the genomes. This method allows the user to calculate lengths of all DNA fragments that are formed after a whole genome is digested at the theoretical recognition sites of a given restriction enzyme. A comparison of the observed peaks in distribution diagrams with the results from DNA cleavage using several restriction enzymes performed in vitro have shown good correspondence between the theoretical and experimental data in several cases. Here, we applied this approach to the annotated genome of Drosophila virilis which is extremely rich in various repeats. Results: Here we explored the combined approach to perform the restriction analysis of D. virilis DNA. This approach enabled to reveal three abundant medium-sized tandem repeats within the D. virilis genome. While the 225 bp repeats were revealed previously in intergenic non-transcribed spacers between ribosomal genes of D. virilis, two other families comprised of 154 bp and 172 bp repeats were not described. Tandem Repeats Finder search demonstrated that 154 bp and 172 bp units are organized in multiple clusters in the genome of D. virilis. Characteristically, only 154 bp repeats derived from Helitron transposon are transcribed. Conclusion: Using in silico digestion in combination with conventional restriction analysis and sequencing of repeated DNA fragments enabled us to isolate and characterize three highly abundant families of medium-sized repeats present in the D. virilis genome. These repeats comprise a significant portion of the genome and may have important roles in genome function and structural integrity. Therefore, we demonstrated an approach which makes possible to investigate in detail the gross arrangement and expression of medium-sized repeats basing on sequencing data even in the case of incompletely assembled and/or annotated genomes.

Pericentromeric heterochromatin in Drosophila generally consists of repetitive DNA, forming the environment associated with gene silencing. Despite the expanding knowledge of the impact of transposable elements (TEs) on the host genome, little is known about the evolution of pericentromeric heterochromatin, its structural composition, and age. During the evolution of the Drosophilidae, hundreds of genes have become embedded within pericentromeric regions yet retained activity. We investigated a pericentromeric heterochromatin fragment found in D. virilis and related species, describing the evolution of genes in this region and the age of TE invasion. Regardless of the heterochromatic environment, the amino acid composition of the genes is under purifying selection. However, the selective pressure affects parts of genes in varying degrees, resulting in expansion of gene introns due to TEs invasion. According to the divergence of TEs, the pericentromeric heterochromatin of the species of virilis group began to form more than 20 million years ago by invasions of retroelements, miniature inverted repeat transposable elements (MITEs), and Helitrons. Importantly, invasions into the heterochromatin continue to occur by TEs that fall under the scope of piRNA silencing. Thus, the pericentromeric heterochromatin, in spite of its ability to induce silencing, has the means for being dynamic, incorporating the regions of active transcription.

Piwi-interacting RNAs (piRNAs) play a crucial role in silencing transposable elements (TEs) in the germ cells of Metazoa by acting as sequence-specific guides. Originating from distinct genomic loci, called piRNA clusters, piRNA can trigger an epigenetic conversion of TE insertions into piRNA clusters by means of a paramutation-like process. However, the variability in piRNA clusters’ capacity to induce such conversions remains poorly understood. Here, we investigated two Drosophila virilis strains with differing capacities to produce piRNAs from the subtelomeric RhoGEF3 and Adar gene loci. We found that active piRNA generation correlates with high levels of the heterochromatic mark histone 3 lysine 9 trimethylation (H3K9me3) over genomic regions that give rise to piRNAs. Importantly, the maternal transmission of piRNAs drives their production in the progeny, even from homologous loci previously inactive in piRNA biogenesis. The RhoGEF3 locus, once epigenetically converted, maintained enhanced piRNA production in subsequent generations lacking the original allele carrying the active piRNA cluster. In contrast, piRNA expression from the converted Adar locus was lost in offspring lacking the inducer allele. The present findings suggest that the paramutation-like behavior of piRNA clusters may be influenced not only by piRNAs but also by structural features and the chromatin environment in the proximity to telomeres, providing new insights into the epigenetic regulation of the Drosophila genome.

Pericentric heterochromatin in Drosophila is generally composed of repetitive DNA forming a transcriptionally repressive environment. Nevertheless, dozens of genes were embedded into pericentric genome regions during evolution of Drosophilidae lineage and retained functional activity. However, factors that contribute to “immunity” of these gene loci to transcriptional silencing remain unknown. Here, we investigated molecular evolution of the essential Myb and Ranbp16 genes. These protein-coding genes reside in euchromatic loci of chromosome X in D. melanogaster and related species, while in other studied Drosophila species, including evolutionary distant ones, they are located in genomic regions highly enriched with the remnants of transposable elements (TEs), suggesting their heterochromatic nature and location. The promoter region of Myb exhibits a conserved structure throughout the Drosophila phylogeny and carries motifs for binding of chromatin remodeling factors, including insulator BEAF-32, regardless of eu- or heterochromatic surroundings. Importantly, BEAF-32 occupies not only the promoter region of Myb but is also found in the vicinity of transcriptional start sites (TSS) of Ranbp16 gene as well as in a wide range of genes located in the contrasting chromatin types in D. melanogaster and D. virilis, denoting the boundary of the nucleosome-free region available for RNA polymerase II recruitment and the surrounding heterochromatin. We also find that along with BEAF-32, insulators dCTCF and GAF are enriched at the TSS of heterochromatic genes in D. melanogaster . Thus, we propose that insulator proteins contribute to gene expression in the heterochromatic environment and, hence, facilitate the evolutionary repositioning of gene loci into heterochromatin. Author summary: Heterochromatin in Drosophila is generally associated with transcriptional silencing. Nevertheless, hundreds of essential genes have been identified in the pericentric heterochromatin of Drosophila melanogaster . Interestingly, genes embedded in pericentric heterochromatin of D. melanogaster may occupy different genomic loci, euchromatic or heterochromatic, due to repositioning in the course of evolution of Drosophila species. By surveying factors that contribute to the normal functioning of the relocated genes in distant Drosophila species, i.e. D. melanogaster and D. virilis , we identify certain insulator proteins (e.g.BEAF-32) that facilitate the expression of heterochromatic genes in spite of the repressive environment.