Evolutionary Dynamics of the Pericentromeric Heterochromatin in Drosophila virilis and Related Species.

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.

Adaptation of gene loci to heterochromatin in the course of Drosophila evolution is associated with insulator proteins.

Pericentromeric heterochromatin is generally composed of repetitive DNA forming a transcriptionally repressive environment. Dozens of genes were embedded into pericentromeric heterochromatin during evolution of Drosophilidae lineage while retaining activity. However, factors that contribute to insusceptibility of gene loci to transcriptional silencing remain unknown. Here, we find that the promoter region of genes that can be embedded in both euchromatin and heterochromatin exhibits a conserved structure throughout the Drosophila phylogeny and carries motifs for binding of certain chromatin remodeling factors, including insulator proteins. Using ChIP-seq data, we demonstrate that evolutionary gene relocation between euchromatin and pericentric heterochromatin occurred with preservation of sites of insulation of BEAF-32 in evolutionarily distant species, i.e. D. melanogaster and D. virilis. Moreover, promoters of virtually all protein-coding genes located in heterochromatin in D. melanogaster are enriched with insulator proteins BEAF-32, GAF and dCTCF. Applying RNA-seq of a BEAF-32 mutant, we show that the impairment of BEAF-32 function has a complex effect on gene expression in D. melanogaster, affecting even those genes that lack BEAF-32 association in their promoters. We propose that conserved intrinsic properties of genes, such as sites of insulation near the promoter regions, may contribute to adaptation of genes to the heterochromatic environment and, hence, facilitate the evolutionary relocation of genes loci between euchromatin and heterochromatin.

Spontaneous gain of susceptibility suggests a novel mechanism of resistance to hybrid dysgenesis in Drosophila virilis.

Syndromes of hybrid dysgenesis (HD) have been critical for our understanding of the transgenerational maintenance of genome stability by piRNA. HD in D. virilis represents a special case of HD since it includes simultaneous mobilization of a set of TEs that belong to different classes. The standard explanation for HD is that eggs of the responder strains lack an abundant pool of piRNAs corresponding to the asymmetric TE families transmitted solely by sperm. However, there are several strains of D. virilis that lack asymmetric TEs, but exhibit a "neutral" cytotype that confers resistance to HD. To characterize the mechanism of resistance to HD, we performed a comparative analysis of the landscape of ovarian small RNAs in strains that vary in their resistance to HD mediated sterility. We demonstrate that resistance to HD cannot be solely explained by a maternal piRNA pool that matches the assemblage of TEs that likely cause HD. In support of this, we have witnessed a cytotype shift from neutral (N) to susceptible (M) in a strain devoid of all major TEs implicated in HD. This shift occurred in the absence of significant change in TE copy number and expression of piRNAs homologous to asymmetric TEs. Instead, this shift is associated with a change in the chromatin profile of repeat sequences unlikely to be causative of paternal induction. Overall, our data suggest that resistance to TE-mediated sterility during HD may be achieved by mechanisms that are distinct from the canonical syndromes of HD.

Activity of heat shock genes' promoters in thermally contrasting animal species.

Heat shock gene promoters represent a highly conserved and universal system for the rapid induction of transcription after various stressful stimuli. We chose pairs of mammalian and insect species that significantly differ in their thermoresistance and constitutive levels of Hsp70 to compare hsp promoter strength under normal conditions and after heat shock (HS). The first pair includes the HSPA1 gene promoter of camel (Camelus dromedarius) and humans. It was demonstrated that the camel HSPA1A and HSPA1L promoters function normally in vitro in human cell cultures and exceed the strength of orthologous human promoters under basal conditions. We used the same in vitro assay for Drosophila melanogaster Schneider-2 (S2) cells to compare the activity of the hsp70 and hsp83 promoters of the second species pair represented by Diptera, i.e., Stratiomys singularior and D. melanogaster, which dramatically differ in thermoresistance and the pattern of Hsp70 accumulation. Promoter strength was also monitored in vivo in D. melanogaster strains transformed with constructs containing the S. singularior hsp70 ORF driven either by its own promoter or an orthologous promoter from the D. melanogaster hsp70Aa gene. Analysis revealed low S. singularior hsp70 promoter activity in vitro and in vivo under basal conditions and after HS in comparison with the endogenous promoter in D. melanogaster cells, which correlates with the absence of canonical GAGA elements in the promoters of the former species. Indeed, the insertion of GAGA elements into the S. singularior hsp70 regulatory region resulted in a dramatic increase in promoter activity in vitro but only modestly enhanced the promoter strength in the larvae of the transformed strains. In contrast with hsp70 promoters, hsp83 promoters from both of the studied Diptera species demonstrated high conservation and universality.

Medium-sized tandem repeats represent an abundant component of the Drosophila virilis genome.

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.

Evolution and dynamics of small RNA response to a retroelement invasion in Drosophila.

Although small RNAs efficiently control transposition activity of most transposons in the host genome, such an immune system is not always applicable against a new transposon's invasions. Here, we explored a possibility to introduce potentially mobile copy of the Penelope retroelement previously implicated in hybrid dysgenesis syndrome in Drosophila virilis into the genomes of two distant Drosophila species. The consequences of such introduction were monitored at different phases after experimental colonization as well as in D. virilis species, which is apparently in the process of ongoing Penelope invasion. We investigated the expression of Penelope and biogenesis of Penelope-derived small RNAs in D. virilis and D. melanogaster strains originally lacking active copies of this element after experimental Penelope invasion. These strains were transformed by constructs containing intact Penelope copies. We show that immediately after transformation, which imitates the first stage of retroelement invasion, Penelope undergoes transposition predominantly in somatic tissues, and may produce siRNAs that are apparently unable to completely silence its activity. However, at the later stages of colonization Penelope copies may jump into one of the piRNA-clusters, which results in production of homologous piRNAs that are maternally deposited and can silence euchromatic transcriptionally active copies of Penelope in trans and, hence, prevent further amplification of the invader in the host genome. Intact Penelope copies and different classes of Penelope-derived small RNAs were found in most geographical strains of D. virilis collected throughout the world. Importantly, all strains of this species containing full-length Penelope tested do not produce gonadal sterility in dysgenic crosses and, hence, exhibit neutral cytotype. To understand whether RNA interference mechanism able to target Penelope operates in related species of the virilis group, we correlated the presence of full-length and potentially active Penelope with the occurrence of piRNAs homologous to this transposable element in the ovaries of species comprising the group. It was demonstrated that Penelope-derived piRNAs are present in all virilis group species containing full-length but transcriptionally silent copies of this element that probably represent the remnants of its previous invasions taking place in the course of the virilis species divergent evolution.

Ontogenetic consequences of dysgenic crosses in Drosophila virilis.

Hybrid dysgenesis (HD) syndrome in Drosophila virilis presumably results from the mobilization of several unrelated mobile genetic elements in dysgenic hybrids. Morphogenetic events during oogenesis and spermatogenesis were investigated in detail in the progeny of D. virilis dysgenic crosses. Using germ-cell specific anti-Vasa staining, we monitored the fate of germline cells at different ontogenetic stages in strains of D. virilis and their hybrids. Anti-Vasa staining indicated that the major loss of pole cells occurs in dysgenic embryos at stage 11-14 after primordial germ cells (PGC) pass the midgut wall. At later ontogenetic stages, including larvae, pupae and imagoes, we often observed an abnormal development of gonads in dysgenic individuals with a frequent occurrence of unilateral and bilateral gonadal atrophy. Dysgenic females were characterized by the presence of various sterile ovarian phenotypes that predominantly include agametic ovarioles, while other atypical forms such as tumor-like ovarioles and dorsalized ovariolar follicles may also be present. Testis abnormalities were also frequently observed in dysgenic males. The sterility manifestations depended on the strain, the growing temperature and the age of the flies used in crosses. The observed gonadal sterility and other HD manifestations correlated with the absence of maternal piRNAs homologous to Penelope and other transposons in the early dysgenic embryos. We speculate that gonadal abnormalities mimicking several known sterility mutations probably result from the disturbance of developmental gene expression machinery due to the activation of unrelated families of transposons in early dysgenic embryos.

Expression of Drosophila virilis retroelements and role of small RNAs in their intrastrain transposition.

Transposition of two retroelements (Ulysses and Penelope) mobilized in the course of hybrid dysgenesis in Drosophila virilis has been investigated by in situ hybridization on polytene chromosomes in two D. virilis strains of different cytotypes routinely used to get dysgenic progeny. The analysis has been repeatedly performed over the last two decades, and has revealed transpositions of Penelope in one of the strains, while, in the other strain, the LTR-containing element Ulysses was found to be transpositionally active. The gypsy retroelement, which has been previously shown to be transpositionally inactive in D. virilis strains, was also included in the analysis. Whole mount is situ hybridization with the ovaries revealed different subcellular distribution of the transposable elements transcripts in the strains studied. Ulysses transpositions occur only in the strain where antisense piRNAs homologous to this TE are virtually absent and the ping-pong amplification loop apparently does not take place. On the other hand small RNAs homologous to Penelope found in the other strain, belong predominantly to the siRNA category (21nt), and consist of sense and antisense species observed in approximately equal proportion. The number of Penelope copies in the latter strain has significantly increased during the last decades, probably because Penelope-derived siRNAs are not maternally inherited, while the low level of Penelope-piRNAs, which are faithfully transmitted from mother to the embryo, is not sufficient to silence this element completely. Therefore, we speculate that intrastrain transposition of the three retroelements studied is controlled predominantly at the post-transcriptional level.

Small RNA-based silencing strategies for transposons in the process of invading Drosophila species.

Colonization of a host by an active transposon can increase mutation rates or cause sterility, a phenotype termed hybrid dysgenesis. As an example, intercrosses of certain Drosophila virilis strains can produce dysgenic progeny. The Penelope element is present only in a subset of laboratory strains and has been implicated as a causative agent of the dysgenic phenotype. We have also introduced Penelope into Drosophila melanogaster, which are otherwise naive to the element. We have taken advantage of these natural and experimentally induced colonization processes to probe the evolution of small RNA pathways in response to transposon challenge. In both species, Penelope was predominantly targeted by endo-small-interfering RNAs (siRNAs) rather than by piwi-interacting RNAs (piRNAs). Although we do observe correlations between Penelope transcription and dysgenesis, we could not correlate differences in maternally deposited Penelope piRNAs with the sterility of progeny. Instead, we found that strains that produced dysgenic progeny differed in their production of piRNAs from clusters in subtelomeric regions, possibly indicating that changes in the overall piRNA repertoire underlie dysgenesis. Considered together, our data reveal unexpected plasticity in small RNA pathways in germ cells, both in the character of their responses to invading transposons and in the piRNA clusters that define their ability to respond to mobile elements.

Penelope retroelements from Drosophila virilis are active after transformation of Drosophila melanogaster.

The Penelope family of retroelements was first described in species of the Drosophila virilis group. Intact elements encode a reverse transcriptase and an endonuclease of the UvrC type, which may play a role in Penelope integration. Penelope is a key element in the induction of D. virilis hybrid dysgenesis, which involves the mobilization of several unrelated families of transposable elements. We here report the successful introduction of Penelope into the germ line of Drosophila melanogaster by P element-mediated transformation with three different constructs. Penelope is actively transcribed in the D. melanogaster genome only in lines transformed with a construct containing a full-length Penelope clone. The transcript is identical to that detected in D. virilis dysgenic hybrids. Most newly transposed Penelope elements have a very complex organization. Significant proliferation of Penelope copy number occurred in some lines during the 24-month period after transformation. The absence of copy number increase with two other constructs suggests that the 5' andor 3' UTRs of Penelope are required for successful transposition in D. melanogaster. No insect retroelement has previously been reported to be actively transcribed and to increase in copy number after interspecific transformation.