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.

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.

To investigate the genetic basis of differing thermotolerance in the closely related species Drosophila virilis and Drosophila lummei, which replace one another along a latitudinal cline, we characterized the hsp70 gene cluster in multiple strains of both species. In both species, all hsp70 copies cluster in a single chromosomal locus, 29C1, and each cluster includes two hsp70 genes arranged as an inverted pair, the ancestral condition. The total number of hsp70 copies is maximally seven in the more thermotolerant D. virilis and five in the less tolerant D. lummei, with some strains of each species exhibiting lower copy numbers. Thus, maximum hsp70 copy number corresponds to hsp70 mRNA and Hsp70 protein levels reported previously and the size of heat-induced puffs at 29C1. The nucleotide sequence and spacing of the hsp70 copies are consistent with tandem duplication of the hsp70 genes in a common ancestor of D. virilis and D. lummei followed by loss of hsp70 genes in D. lummei. These and other data for hsp70 in Drosophila suggest that evolutionary adaptation has repeatedly modified hsp70 copy number by several different genetic mechanisms.

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.

No abstract