Variation in piRNA and transposable element content in strains of Drosophila melanogaster

Transposable Element (TE) insertion predictions from RNAseq inputs and TE impact on RNA splicing and gene expression in Drosophila brain transcriptomes

Recent studies have suggested that Transposable Elements (TEs) residing in introns frequently splice into and alter primary gene-coding transcripts. To re-examine the exonization frequency of TEs into protein-coding gene transcripts, we re-analyzed a Drosophila neuron circadian rhythm RNAseq dataset and a deep long RNA fly midbrain RNAseq dataset using our Transposon Insertion and Depletion Analyzer (TIDAL) program. Our TIDAL results were able to predict several TE insertions from RNAseq data that were consistent with previous published studies. However, we also uncovered many discrepancies in TE-exonization calls, such as reads that mainly support intron retention of the TE and little support for chimeric mRNA spliced to the TE. We then deployed rigorous genomic DNA-PCR (gDNA-PCR) and RT-PCR procedures on TE-mRNA fusion candidates to see how many of bioinformatics predictions could be validated. By testing a w1118 strain from which the deeper long RNAseq data was derived and comparing to an OreR strain, only 9 of 23 TIDAL candidates (< 40%) could be validated as a novel TE insertion by gDNA-PCR, indicating that deeper study is needed when using RNAseq data as inputs into current TE-insertion prediction programs. Of these validated calls, our RT-PCR results only supported TE-intron retention. Lastly, in the Dscam2 and Bx genes of the w1118 strain that contained intronic TEs, gene expression was 23 times higher than the OreR genes lacking the TEs. This study's validation approach indicates that chimeric TE-mRNAs are infrequent and cautions that more optimization is required in bioinformatics programs to call TE insertions using RNAseq datasets.

The genomes of all lungfish inform on genome expansion and tetrapod evolution

The genomes of living lungfishes can inform on the molecular-developmental basis of the Devonian sarcopterygian fish-tetrapod transition. We de novo sequenced the genomes of the African (Protopterus annectens) and South American lungfishes (Lepidosiren paradoxa). The Lepidosiren genome (about 91 Gb, roughly 30 times the human genome) is the largest animal genome sequenced so far and more than twice the size of the Australian (Neoceratodus forsteri)1 and African2 lungfishes owing to enlarged intergenic regions and introns with high repeat content (about 90%). All lungfish genomes continue to expand as some transposable elements (TEs) are still active today. In particular, Lepidosiren's genome grew extremely fast during the past 100 million years (Myr), adding the equivalent of one human genome every 10 Myr. This massive genome expansion seems to be related to a reduction of PIWI-interacting RNAs and C2H2 zinc-finger and Krüppel-associated box (KRAB)-domain protein genes that suppress TE expansions. Although TE abundance facilitates chromosomal rearrangements, lungfish chromosomes still conservatively reflect the ur-tetrapod karyotype. Neoceratodus' limb-like fins still resemble those of their extinct relatives and remained phenotypically static for about 100 Myr. We show that the secondary loss of limb-like appendages in the Lepidosiren-Protopterus ancestor was probably due to loss of sonic hedgehog limb-specific enhancers.

piRNA-Guided Transposon Silencing and Response to Stress in Drosophila Germline

Transposons are integral genome constituents that can be domesticated for host functions, but they also represent a significant threat to genome stability. Transposon silencing is especially critical in the germline, which is dedicated to transmitting inherited genetic material. The small Piwi-interacting RNAs (piRNAs) have a deeply conserved function in transposon silencing in the germline. piRNA biogenesis and function are particularly well understood in Drosophila melanogaster, but some fundamental mechanisms remain elusive and there is growing evidence that the pathway is regulated in response to genotoxic and environmental stress. Here, we review transposon regulation by piRNAs and the piRNA pathway regulation in response to stress, focusing on the Drosophila female germline.

piOxi database: a web resource of germline and somatic tissue piRNAs identified by chemical oxidation

PIWI-interacting RNAs (piRNAs) are a class of small non-coding RNAs that are highly expressed and extensively studied from the germline. piRNAs associate with PIWI proteins to maintain DNA methylation for transposon silencing and transcriptional gene regulation for genomic stability. Mature germline piRNAs have distinct characteristics including a 24- to 32-nucleotide length and a 2'-O-methylation signature at the 3' end. Although recent studies have identified piRNAs in somatic tissues, they remain poorly characterized. For example, we recently demonstrated notable expression of piRNA in the murine soma, and while overall expression was lower than that of the germline, unique characteristics suggested tissue-specific functions of this class. While currently available databases commonly use length and association with PIWI proteins to identify piRNA, few have included a chemical oxidation method that detects piRNA based on its 3' modification. This method leads to reproducible and rigorous data processing when coupled with next-generation sequencing and bioinformatics analysis. Here, we introduce piOxi DB, a user-friendly web resource that provides a comprehensive analysis of piRNA, generated exclusively through sodium periodate treatment of small RNA. The current version of piOxi DB includes 435 749 germline and 9828 somatic piRNA sequences robustly identified from M. musculus, M. fascicularis and H. sapiens. The database provides species- and tissue-specific data that are further analyzed according to chromosome location and correspondence to gene and repetitive elements. piOxi DB is an informative tool to assist broad research applications in the fields of RNA biology, cancer biology, environmental toxicology and beyond. Database URL:  https://pioxidb.dcmb.med.umich.edu/.

The composition of piRNA clusters in Drosophila melanogaster deviates from expectations under the trap model

BACKGROUND

It is widely assumed that the invasion of a transposable element (TE) in mammals and invertebrates is stopped when a copy of the TE jumps into a piRNA cluster (i.e., the trap model). However, recent works, which for example showed that deletion of three major piRNA clusters has no effect on TE activity, cast doubt on the trap model.

RESULTS

Here, we test the trap model from a population genetics perspective. Our simulations show that the composition of regions that act as transposon traps (i.e., potentially piRNA clusters) ought to deviate from regions that have no effect on TE activity. We investigated TEs in five Drosophila melanogaster strains using three complementary approaches to test whether the composition of piRNA clusters matches these expectations. We found that the abundance of TE families inside and outside of piRNA clusters is highly correlated, although this is not expected under the trap model. Furthermore, the distribution of the number of TE insertions in piRNA clusters is also much broader than expected.

CONCLUSIONS

We found that the observed composition of piRNA clusters is not in agreement with expectations under the simple trap model. Dispersed piRNA producing TE insertions and temporal as well as spatial heterogeneity of piRNA clusters may account for these deviations.

Structural insights into the chromodomain of Oxpecker in complex with histone H3 lysine 9 trimethylation reveal a transposon silencing mechanism by heterodimerization

Oxpecker, the homolog of Rhino/HP1D, exclusively expressed in Drosophila ovaries, belongs to the Heterochromatin Protein 1 family, as does Rhino. Rhi recognizes piRNA clusters enriched with the heterochromatin marker H3K9me3 via its N-terminal chromodomain and recruits Deadlock via its C-terminal chromoshadow domain, further recruits Moonshiner, a paralog of the TATA box-binding protein-related factor 2 large subunits, to promote transcription of piRNA precursors, thereby protecting the genome. Despite Oxp possessing only the chromodomain, its loss leads to the upregulation of transposons in the female germline. In this study, we solved the crystal structure of the Oxp chromodomain in complex with the histone H3K9me3 peptide. As the Oxp chromodomain dimerizes, two H3K9me3 peptides bind to the Oxp chromodomain in an antiparallel manner. ITC experiments and site-directed mutagenesis experiments showed that E44 determines Oxp's five-fold stronger binding ability to H3K9me3 than that of Rhi. In addition, we found that Oxp and Rhi can form a heterodimer, which may shed light on the molecular mechanism by which Oxp regulates transposon silencing in the absence of CSD.

Epigenetic targeting of transposon relics: beating the dead horses of the genome?

Transposable elements (TEs) have been seen as selfish genetic elements that can propagate in a host genome. Their propagation success is however hindered by a combination of mechanisms such as mutations, selection, and their epigenetic silencing by the host genome. As a result, most copies of TEs in a given genome are dead relics: their sequence is too degenerated to allow any transposition. Nevertheless, these TE relics often, but not always, remain epigenetically silenced, and if not to prevent transposition anymore, one can wonder the reason for this phenomenon. The mere self-perpetuating loop inherent to epigenetic silencing could alone explain that even when inactive, TE copies remain silenced. Beyond this process, nevertheless, antagonistic selective forces are likely to act on TE relic silencing. Especially, without the benefit of preventing transposition, TE relic silencing may prove deleterious to the host fitness, suggesting that the maintenance of TE relic silencing is the result of a fine, and perhaps case-by-case, evolutionary trade-off between beneficial and deleterious effects. Ultimately, the release of TE relics silencing may provide a 'safe' ground for adaptive epimutations to arise. In this review, we provide an overview of these questions in both plants and animals.

High-quality, chromosome-scale genome assemblies: comparisons of three Diaphorina citri (Asian citrus psyllid) geographic populations

The Asian citrus psyllid, Diaphorina citri, is the insect vector of the causal agent of huanglongbing (HLB), a devastating bacterial disease of commercial citrus. Presently, few genomic resources exist for D. citri. In this study, we utilized PacBio HiFi and chromatin confirmation contact (Hi-C) sequencing to sequence, assemble, and compare three high-quality, chromosome-scale genome assemblies of D. citri collected from California, Taiwan, and Uruguay. Our assemblies had final sizes of 282.67 Mb (California), 282.89 Mb (Taiwan), and 266.67 Mb (Uruguay) assembled into 13 pseudomolecules—a reduction in assembly size of 41–45% compared with previous assemblies which we validated using flow cytometry. We identified the X chromosome in D. citri and annotated each assembly for repetitive elements, protein-coding genes, transfer RNAs, ribosomal RNAs, piwi-interacting RNA clusters, and endogenous viral elements. Between 19,083 and 20,357 protein-coding genes were predicted. Repetitive DNA accounts for 36.87–38.26% of each assembly. Comparative analyses and mitochondrial haplotype networks suggest that Taiwan and Uruguay D. citri are more closely related, while California D. citri are closely related to Florida D. citri. These high-quality, chromosome-scale assemblies provide new genomic resources to researchers to further D. citri and HLB research.

Transposable Element Expression and Regulation Profile in Gonads of Interspecific Hybrids of Drosophila&nbsp;arizonae and Drosophila mojavensis wrigleyi

Interspecific hybridization may lead to sterility and/or inviability through differential expression of genes and transposable elements (TEs). In Drosophila, studies have reported massive TE mobilization in hybrids from interspecific crosses of species presenting high divergence times. However, few studies have examined the consequences of TE mobilization upon hybridization in recently diverged species, such as Drosophila arizonae and D. mojavensis. We have sequenced transcriptomes of D. arizonae and the subspecies D. m. wrigleyi and their reciprocal hybrids, as well as piRNAs, to analyze the impact of genomic stress on TE regulation. Our results revealed that the differential expression in both gonadal tissues of parental species was similar. Globally, ovaries and testes showed few deregulated TEs compared with both parental lines. Analyses of small RNA data showed that in ovaries, the TE upregulation is likely due to divergence of copies inherited from parental genomes and lack of piRNAs mapping to them. Nevertheless, in testes, the divergent expression of genes associated with chromatin state and piRNA pathway potentially indicates that TE differential expression is related to the divergence of regulatory genes that play a role in modulating transcriptional and post-transcriptional mechanisms.

piRNA and Transposon Dynamics in Drosophila: A Female Story

The germlines of metazoans contain transposable elements (TEs) causing genetic instability and affecting fitness. To protect the germline from TE activity, gonads of metazoans produce TE-derived PIWI-interacting RNAs (piRNAs) that silence TE expression. In Drosophila, our understanding of piRNA biogenesis is mainly based on studies of the Drosophila melanogaster female germline. However, it is not known whether piRNA functions are also important in the male germline or whether and how piRNAs are affected by the global genomic context. To address these questions, we compared genome sequences, transcriptomes, and small RNA libraries extracted from entire testes and ovaries of two sister species: D. melanogaster and Drosophila simulans. We found that most TE-derived piRNAs were produced in ovaries and that piRNA pathway genes were strongly overexpressed in ovaries compared with testes, indicating that the silencing of TEs by the piRNA pathway mainly took place in the female germline. To study the relationship between host piRNAs and TE landscape, we analyzed TE genomic features and how they correlate with piRNA production in the two species. In D. melanogaster, we found that TE-derived piRNAs target recently active TEs. In contrast, although Drosophila simulans TEs do not display any features of recent activity, the host still intensively produced silencing piRNAs targeting old TE relics. Together, our results show that the piRNA silencing response mainly takes place in Drosophila ovaries and indicate that the host piRNA response is implemented following a burst of TE activity and could persist long after the extinction of active TE families.

A bird-like genome from a frog: Mechanisms of genome size reduction in the ornate burrowing frog, Platyplectrum ornatum

The diversity of genome sizes across the tree of life is of key interest in evolutionary biology. Various correlates of variation in genome size, such as accumulation of transposable elements (TEs) or rate of DNA gain and loss, are well known, but the underlying molecular mechanisms driving or constraining genome size are poorly understood. Here, we study one of the smallest genomes among frogs characterized thus far, that of the ornate burrowing frog (Platyplectrum ornatum) from Australia, and compare it to other published frog and vertebrate genomes to examine the forces driving reduction in genome size. At ∼1.06 gigabases (Gb), the P. ornatum genome is like that of birds, revealing four major mechanisms underlying TE dynamics: reduced abundance of all major classes of TEs; increased net deletion bias in TEs; drastic reduction in intron lengths; and expansion via gene duplication of the repertoire of TE-suppressing Piwi genes, accompanied by increased expression of Piwi-interacting RNA (piRNA)-based TE-silencing pathway genes in germline cells. Transcriptomes from multiple tissues in both sexes corroborate these results and provide insight into sex-differentiation pathways in Platyplectrum Genome skimming of two closely related frog species (Lechriodus fletcheri and Limnodynastes fletcheri) confirms a reduction in TEs as a major driver of genome reduction in Platyplectrum and supports a macroevolutionary scenario of small genome size in frogs driven by convergence in life history, especially rapid tadpole development and tadpole diet. The P. ornatum genome offers a model for future comparative studies on mechanisms of genome size reduction in amphibians and vertebrates generally.

A mosquito small RNA genomics resource reveals dynamic evolution and host responses to viruses and transposons

Although mosquitoes are major transmission vectors for pathogenic arboviruses, viral infection has little impact on mosquito health. This immunity is caused in part by mosquito RNA interference (RNAi) pathways that generate antiviral small interfering RNAs (siRNAs) and Piwi-interacting RNAs (piRNAs). RNAi also maintains genome integrity by potently repressing mosquito transposon activity in the germline and soma. However, viral and transposon small RNA regulatory pathways have not been systematically examined together in mosquitoes. Therefore, we developed an integrated mosquito small RNA genomics (MSRG) resource that analyzes the transposon and virus small RNA profiles in mosquito cell cultures and somatic and gonadal tissues across four medically important mosquito species. Our resource captures both somatic and gonadal small RNA expression profiles within mosquito cell cultures, and we report the evolutionary dynamics of a novel Mosquito-Conserved piRNA Cluster Locus (MCpiRCL) made up of satellite DNA repeats. In the larger culicine mosquito genomes we detected highly regular periodicity in piRNA biogenesis patterns coinciding with the expansion of Piwi pathway genes. Finally, our resource enables detection of cross talk between piRNA and siRNA populations in mosquito cells during a response to virus infection. The MSRG resource will aid efforts to dissect and combat the capacity of mosquitoes to tolerate and spread arboviruses.

Telomeric TART elements target the piRNA machinery in Drosophila

Coevolution between transposable elements (TEs) and their hosts can be antagonistic, where TEs evolve to avoid silencing and the host responds by reestablishing TE suppression, or mutualistic, where TEs are co-opted to benefit their host. The TART-A TE functions as an important component of Drosophila telomeres but has also reportedly inserted into the Drosophila melanogaster nuclear export factor gene nxf2. We find that, rather than inserting into nxf2, TART-A has actually captured a portion of nxf2 sequence. We show that TART-A produces abundant Piwi-interacting small RNAs (piRNAs), some of which are antisense to the nxf2 transcript, and that the TART-like region of nxf2 is evolving rapidly. Furthermore, in D. melanogaster, TART-A is present at higher copy numbers, and nxf2 shows reduced expression, compared to the closely related species Drosophila simulans. We propose that capturing nxf2 sequence allowed TART-A to target the nxf2 gene for piRNA-mediated repression and that these 2 elements are engaged in antagonistic coevolution despite the fact that TART-A is serving a critical role for its host genome.

Co-evolution between transposable elements (TEs) and their hosts can be antagonistic, where TEs evolve to avoid silencing and the host responds by re-establishing TE suppression, or mutualistic, where TEs are co-opted to benefit their host. This study shows that a specialized Drosophila retrotransposon that functions as a telomere has captured a portion of a host piRNA gene which may allow it to evade silencing.

A Transposon Story: From TE Content to TE Dynamic Invasion of Drosophila Genomes Using the Single-Molecule Sequencing Technology from Oxford Nanopore

Transposable elements (TEs) are the main components of genomes. However, due to their repetitive nature, they are very difficult to study using data obtained with short-read sequencing technologies. Here, we describe an efficient pipeline to accurately recover TE insertion (TEI) sites and sequences from long reads obtained by Oxford Nanopore Technology (ONT) sequencing. With this pipeline, we could precisely describe the landscapes of the most recent TEIs in wild-type strains of Drosophila melanogaster and Drosophila simulans. Their comparison suggests that this subset of TE sequences is more similar than previously thought in these two species. The chromosome assemblies obtained using this pipeline also allowed recovering piRNA cluster sequences, which was impossible using short-read sequencing. Finally, we used our pipeline to analyze ONT sequencing data from a D. melanogaster unstable line in which LTR transposition was derepressed for 73 successive generations. We could rely on single reads to identify new insertions with intact target site duplications. Moreover, the detailed analysis of TEIs in the wild-type strains and the unstable line did not support the trap model claiming that piRNA clusters are hotspots of TE insertions.

Taming the Turmoil Within: New Insights on the Containment of Transposable Elements

Transposable elements (TEs) are mobile genetic parasites that can exponentially increase their genomic abundance through self-propagation. Classic theoretical papers highlighted the importance of two potentially escalating forces that oppose TE spread: regulated transposition and purifying selection. Here, we review new insights into mechanisms of TE regulation and purifying selection, which reveal the remarkable foresight of these theoretical models. We further highlight emergent connections between transcriptional control enacted by small RNAs and the contribution of TE insertions to structural mutation and host-gene regulation. Finally, we call for increased comparative analysis of TE dynamics and fitness effects, as well as host control mechanisms, to reveal how interconnected forces shape the differential prevalence and distribution of TEs across the tree of life.

Rapid evolution of piRNA-mediated silencing of an invading transposable element was driven by abundant de novo mutations

The regulation of transposable element (TE) activity by small RNAs is a ubiquitous feature of germlines. However, despite the obvious benefits to the host in terms of ensuring the production of viable gametes and maintaining the integrity of the genomes they carry, it remains controversial whether TE regulation evolves adaptively. We examined the emergence and evolutionary dynamics of repressor alleles after P-elements invaded the Drosophila melanogaster genome in the mid-twentieth century. In many animals including Drosophila, repressor alleles are produced by transpositional insertions into piRNA clusters, genomic regions encoding the Piwi-interacting RNAs (piRNAs) that regulate TEs. We discovered that ∼94% of recently collected isofemale lines in the Drosophila melanogaster Genetic Reference Panel (DGRP) contain at least one P-element insertion in a piRNA cluster, indicating that repressor alleles are produced by de novo insertion at an exceptional rate. Furthermore, in our sample of approximately 200 genomes, we uncovered no fewer than 80 unique P-element insertion alleles in at least 15 different piRNA clusters. Finally, we observe no footprint of positive selection on P-element insertions in piRNA clusters, suggesting that the rapid evolution of piRNA-mediated repression in D. melanogaster was driven primarily by mutation. Our results reveal for the first time how the unique genetic architecture of piRNA production, in which numerous piRNA clusters can encode regulatory small RNAs upon transpositional insertion, facilitates the nonadaptive rapid evolution of repression.

Nuclear Ccr4-Not mediates the degradation of telomeric and transposon transcripts at chromatin in the Drosophila germline

Ccr4-Not is a highly conserved complex involved in cotranscriptional RNA surveillance pathways in yeast. In Drosophila, Ccr4-Not is linked to the translational repression of miRNA targets and the posttranscriptional control of maternal mRNAs during oogenesis and embryonic development. Here, we describe a new role for the Ccr4-Not complex in nuclear RNA metabolism in the Drosophila germline. Ccr4 depletion results in the accumulation of transposable and telomeric repeat transcripts in the fraction of chromatin-associated RNA; however, it does not affect small RNA levels or the heterochromatin state of the target loci. Nuclear targets of Ccr4 mainly comprise active full-length transposable elements (TEs) and telomeric and subtelomeric repeats. Moreover, Ccr4-Not foci localize at telomeres in a Piwi-dependent manner, suggesting a functional relationship between these pathways. Indeed, we detected interactions between the components of the Ccr4-Not complex and piRNA machinery, which indicates that these pathways cooperate in the nucleus to recognize and degrade TE transcripts at transcription sites. These data reveal a new layer of transposon control in the germline, which is critical for the maintenance of genome integrity.

Nanopore sequencing and Hi-C scaffolding provide insight into the evolutionary dynamics of transposable elements and piRNA production in wild strains of Drosophila melanogaster

Illumina sequencing has allowed for population-level surveys of transposable element (TE) polymorphism via split alignment approaches, which has provided important insight into the population dynamics of TEs. However, such approaches are not able to identify insertions of uncharacterized TEs, nor can they assemble the full sequence of inserted elements. Here, we use nanopore sequencing and Hi-C scaffolding to produce de novo genome assemblies for two wild strains of Drosophila melanogaster from the Drosophila Genetic Reference Panel (DGRP). Ovarian piRNA populations and Illumina split-read TE insertion profiles have been previously produced for both strains. We find that nanopore sequencing with Hi-C scaffolding produces highly contiguous, chromosome-length scaffolds, and we identify hundreds of TE insertions that were missed by Illumina-based methods, including a novel micropia-like element that has recently invaded the DGRP population. We also find hundreds of piRNA-producing loci that are specific to each strain. Some of these loci are created by strain-specific TE insertions, while others appear to be epigenetically controlled. Our results suggest that Illumina approaches reveal only a portion of the repetitive sequence landscape of eukaryotic genomes and that population-level resequencing using long reads is likely to provide novel insight into the evolutionary dynamics of repetitive elements.

The evolutionary arms race between transposable elements and piRNAs in Drosophila melanogaster

BACKGROUND

The piwi-interacting RNAs (piRNAs) are small non-coding RNAs that specifically repress transposable elements (TEs) in the germline of Drosophila. Despite our expanding understanding of TE:piRNA interaction, whether there is an evolutionary arms race between TEs and piRNAs was unclear.

RESULTS

Here, we studied the population genomics of TEs and piRNAs in the worldwide strains of D. melanogaster. By conducting a correlation analysis between TE contents and the abundance of piRNAs from ovaries of representative strains of D. melanogaster, we find positive correlations between TEs and piRNAs in six TE families. Our simulations further highlight that TE activities and the strength of purifying selection against TEs are important factors shaping the interactions between TEs and piRNAs. Our studies also suggest that the de novo generation of piRNAs is an important mechanism to repress the newly invaded TEs.

CONCLUSIONS

Our results revealed the existence of an evolutionary arms race between the copy numbers of TEs and the abundance of antisense piRNAs at the population level. Although the interactions between TEs and piRNAs are complex and many factors should be considered to impact their interaction dynamics, our results suggest the emergence, repression specificity and strength of piRNAs on TEs should be considered in studying the landscapes of TE insertions in Drosophila. These results deepen our understanding of the interactions between piRNAs and TEs, and also provide novel insights into the nature of genomic conflicts of other forms.

Host-transposon interactions: conflict, cooperation, and cooption

Transposable elements (TEs) are mobile DNA sequences that colonize genomes and threaten genome integrity. As a result, several mechanisms appear to have emerged during eukaryotic evolution to suppress TE activity. However, TEs are ubiquitous and account for a prominent fraction of most eukaryotic genomes. We argue that the evolutionary success of TEs cannot be explained solely by evasion from host control mechanisms. Rather, some TEs have evolved commensal and even mutualistic strategies that mitigate the cost of their propagation. These coevolutionary processes promote the emergence of complex cellular activities, which in turn pave the way for cooption of TE sequences for organismal function.

Rapid evolution of piRNA pathway and its transposon targets in Japanese flounder (Paralichthys olivaceus)

Piwi-interacting RNA (piRNA) pathway is essential for germline specification, gametogenesis, and genome integrity as defense against transposable elements (TEs). This pathway has been suggested to have undergone rapid adaptive evolution in spite of its conserved role in TE silencing. However, with diverse sexual development patterns, piRNA pathway evolution and its adaptation to transposon activity in teleost lineages remain less known. This study illustrated the evolutionary significance of piRNA pathway via a systematic comparative analysis on diverse teleosts, including flatfish lineages. Molecular evolution of piRNA pathway and microRNA/small interfering RNA pathway genes indicated a faster evolution of piRNA pathway in teleosts than in mammals. Positive selection was detected at the PAZ (Piwi-Argonaute-Zwille) domain involved in Piwi-piRNA interaction, thereby suggesting that the amino acid substitutions were adaptive to their functions in teleost piRNA pathway. Notably, Piwil1 evolved faster in Japanese flounder than in other teleosts, and the piRNA pathway genes expressed higher in testes than in ovaries. In addition, gonadal transcriptomic analysis revealed male under-represented TE families mainly from DNA transposons, which were the potential targets of the complex formed by male-biased Piwi genes and male over-represented piRNAs in Japanese flounder testes. The potential piRNA-TE regulatory relationships suggested that the rapidly evolved piRNA pathway in Japanese flounder was likely involved in the regulation of transposon activity in germlines and could play important roles in Japanese flounder gonadal development and spermatogenesis.

A unique cluster of roo insertions in the promoter region of a stress response gene in Drosophila melanogaster

Transposable elements (TEs) are not randomly distributed in the genome. A genome-wide analysis of the D. melanogaster genome found that differences in TE density across 50 kb genomic regions was due both to transposition and duplication. At smaller genomic scales, promoter regions of hsp genes and the promoter region of CG18446 have been shown to accumulate TE insertions. In this work, we have further analyzed the promoter region of CG18446. We screened 218 strains collected in 15 natural populations, and we found that the CG18446 promoter region contains 20 independent roo insertions. Based on phylogenetic analysis, we suggest that the presence of multiple roo insertions in this region is likely to be the result of several bursts of transposition. Moreover, we found that the roo insertional cluster in the CG18446 promoter region is unique: no other promoter region in the genome contains a similar number of roo insertions. We found that, similar to hsp gene promoters, chromatin accessibility could be one of the factors explaining the recurrent insertions of roo elements in CG18446 promoter region.

TEtools facilitates big data expression analysis of transposable elements and reveals an antagonism between their activity and that of piRNA genes

Over recent decades, substantial efforts have been made to understand the interactions between host genomes and transposable elements (TEs). The impact of TEs on the regulation of host genes is well known, with TEs acting as platforms of regulatory sequences. Nevertheless, due to their repetitive nature it is considerably hard to integrate TE analysis into genome-wide studies. Here, we developed a specific tool for the analysis of TE expression: TEtools. This tool takes into account the TE sequence diversity of the genome, it can be applied to unannotated or unassembled genomes and is freely available under the GPL3 (https://github.com/l-modolo/TEtools). TEtools performs the mapping of RNA-seq data obtained from classical mRNAs or small RNAs onto a list of TE sequences and performs differential expression analyses with statistical relevance. Using this tool, we analyzed TE expression from five Drosophila wild-type strains. Our data show for the first time that the activity of TEs is strictly linked to the activity of the genes implicated in the piwi-interacting RNA biogenesis and therefore fits an arms race scenario between TE sequences and host control genes.

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.

Characterization of the piRNA pathway during development of the sea anemone Nematostella vectensis

PIWI-interacting RNAs (piRNAs) and associated proteins comprise a conserved pathway for silencing transposons in metazoan germlines. piRNA pathway components are also expressed in multipotent somatic stem cells in various organisms. piRNA functions have been extensively explored in bilaterian model systems, however, comprehensive studies in non-bilaterian phyla remain limited. Here we investigate the piRNA pathway during the development of Nematostella vectensis, a well-established model system belonging to Cnidaria, the sister group to Bilateria. To date, no population of somatic stem cells has been identified in this organism, despite its long life-span and regenerative capacities that require a constant cell-renewal. We show that Nematostella piRNA pathway components are broadly expressed in early developmental stages, while piRNAs themselves show differential expression, suggesting specific developmental roles of distinct piRNA families. In adults, piRNA associated proteins are enriched in the germline but also expressed in somatic cells, indicating putative stem cell properties. Furthermore, we provide experimental evidence that Nematostella piRNAs cleave transposable elements as well as protein-coding genes. Our results demonstrate that somatic expression of piRNA associated proteins as well as the roles of piRNAs in transposon repression and gene regulation are likely ancestral features that evolved before the split between Cnidaria and Bilateria.

Genomic analysis of P elements in natural populations of Drosophila melanogaster

The Drosophila melanogaster P transposable element provides one of the best cases of horizontal transfer of a mobile DNA sequence in eukaryotes. Invasion of natural populations by the P element has led to a syndrome of phenotypes known as P-M hybrid dysgenesis that emerges when strains differing in their P element composition mate and produce offspring. Despite extensive research on many aspects of P element biology, many questions remain about the genomic basis of variation in P-M dysgenesis phenotypes across populations. Here we compare estimates of genomic P element content with gonadal dysgenesis phenotypes for isofemale strains obtained from three worldwide populations of D. melanogaster to illuminate the molecular basis of natural variation in cytotype status. We show that P element abundance estimated from genome sequences of isofemale strains is highly correlated across different bioinformatics approaches, but that abundance estimates are sensitive to method and filtering strategies as well as incomplete inbreeding of isofemale strains. We find that P element content varies significantly across populations, with strains from a North American population having fewer P elements but a higher proportion of full-length elements than strains from populations sampled in Europe or Africa. Despite these geographic differences in P element abundance and structure, neither the number of P elements nor the ratio of full-length to internally-truncated copies is strongly correlated with the degree of gonadal dysgenesis exhibited by an isofemale strain. Thus, variation in P element abundance and structure across different populations does not necessarily lead to corresponding geographic differences in gonadal dysgenesis phenotypes. Finally, we confirm that population differences in the abundance and structure of P elements that are observed from isofemale lines can also be observed in pool-seq samples from the same populations. Our work supports the view that genomic P element content alone is not sufficient to explain variation in gonadal dysgenesis across strains of D. melanogaster, and informs future efforts to decode the genomic basis of geographic and temporal differences in P element induced phenotypes.

Silencing of Transposable Elements by piRNAs in Drosophila: An Evolutionary Perspective

Transposable elements (TEs) are DNA sequences that can move within the genome. TEs have greatly shaped the genomes, transcriptomes, and proteomes of the host organisms through a variety of mechanisms. However, TEs generally disrupt genes and destabilize the host genomes, which substantially reduce fitness of the host organisms. Understanding the genomic distribution and evolutionary dynamics of TEs will greatly deepen our understanding of the TE-mediated biological processes. Most TE insertions are highly polymorphic in Drosophila melanogaster, providing us a good system to investigate the evolution of TEs at the population level. Decades of theoretical and experimental studies have well established “transposition-selection” population genetics model, which assumes that the equilibrium between TE replication and purifying selection determines the copy number of TEs in the genome. In the last decade, P-element-induced wimpy testis (PIWI)-interacting RNAs (piRNAs) were demonstrated to be master repressors of TE activities in Drosophila. The discovery of piRNAs revolutionized our understanding of TE repression, because it reveals that the host organisms have evolved an adaptive mechanism to defend against TE invasion. Tremendous progress has been made to understand the molecular mechanisms by which piRNAs repress active TEs, although many details in this process remain to be further explored. The interaction between piRNAs and TEs well explains the molecular mechanisms underlying hybrid dysgenesis for the I-R and P-M systems in Drosophila, which have puzzled evolutionary biologists for decades. The piRNA repression pathway provides us an unparalleled system to study the co-evolutionary process between parasites and host organisms.

Natural variation of piRNA expression affects immunity to transposable elements

In the Drosophila germline, transposable elements (TEs) are silenced by PIWI-interacting RNA (piRNA) that originate from distinct genomic regions termed piRNA clusters and are processed by PIWI-subfamily Argonaute proteins. Here, we explore the variation in the ability to restrain an alien TE in different Drosophila strains. The I-element is a retrotransposon involved in the phenomenon of I-R hybrid dysgenesis in Drosophila melanogaster. Genomes of R strains do not contain active I-elements, but harbour remnants of ancestral I-related elements. The permissivity to I-element activity of R females, called reactivity, varies considerably in natural R populations, indicating the existence of a strong natural polymorphism in defense systems targeting transposons. To reveal the nature of such polymorphisms, we compared ovarian small RNAs between R strains with low and high reactivity and show that reactivity negatively correlates with the ancestral I-element-specific piRNA content. Analysis of piRNA clusters containing remnants of I-elements shows increased expression of the piRNA precursors and enrichment by the Heterochromatin Protein 1 homolog, Rhino, in weak R strains, which is in accordance with stronger piRNA expression by these regions. To explore the nature of the differences in piRNA production, we focused on two R strains, weak and strong, and showed that the efficiency of maternal inheritance of piRNAs as well as the I-element copy number are very similar in both strains. At the same time, germline and somatic uni-strand piRNA clusters generate more piRNAs in strains with low reactivity, suggesting the relationship between the efficiency of primary piRNA production and variable response to TE invasions. The strength of adaptive genome defense is likely driven by naturally occurring polymorphisms in the rapidly evolving piRNA pathway proteins. We hypothesize that hyper-efficient piRNA production is contributing to elimination of a telomeric retrotransposon HeT-A, which we have observed in one particular transposon-resistant R strain.

Transposon activity in the germline is suppressed by the PIWI-interacting RNA (piRNA) pathway. The resistance of natural Drosophila strains to transposon invasion varies considerably, but the nature of this variability is unknown. We discovered that natural variation in the efficiency of primary piRNA production in the germline causes dramatic differences in the susceptibility to expansion of a newly invaded transposon. A high level of piRNA production in the germline is achieved by increased expression of piRNA precursors. In one of the most transposon-resistant strains, increased content of primary piRNA is observed in both the germline and ovarian somatic cells. We suggest that polymorphisms in piRNA pathway factors are responsible for increased piRNA production. piRNA pathway proteins have been shown to be evolving rapidly under selective pressure. Our data are the first to describe a phenotype that might be caused by this kind of polymorphism. We also demonstrate a likely explanation as to why an overly active piRNA pathway can cause more harm than good in Drosophila: Highly efficient piRNA processing leads to elimination of domesticated telomeric retrotransposons essential for telomere elongation, an effect which has been observed in a natural strain that is extremely resistant to transposon invasion.

Transposable Element Bm1645 is a Source of BmAGO2-associated Small RNAs that affect its expression in Bombyx mori

Background

A transposable element (TE) is a DNA fragment that can change its position within a genome. Transposable elements play important roles in maintaining the stability and diversity of organisms by transposition. Recent studies have shown that approximately half of the genes in Bombyx mori are TEs.

Results

We systematically identified and analyzed the BmAGO2-associated TEs, which exceed 100 in the B. mori genome. Additionally, we also mapped the small RNAs associated with BmAGO2 in B.mori. The transposon Bm1645 is the most abundant TE associated with BmAGO2, and Bm1645-derived small RNAs represent a small RNA pool. We determined the expression patterns of several Bm1645-derived small RNAs by northern blotting, and the results showed there was differential expression of multiple small RNAs in normal and BmNPV-infected BmN cells and silkworms from various developmental stages. We confirmed that four TE-siRNAs could bind to BmAGO2 using EMSA and also validated the recognition sites of these four TE-siRNAs in Bm1645 by dual-luciferase reporter assays. Furthermore, qRT-PCR analysis revealed the overexpression of the four TE-siRNAs could downregulate the expression of Bm1645 in BmN cells, and the transcription of Bm1645 was upregulated by the downregulation of BmAGO2.

Conclusions

Our results suggest Bm1645 functions as a source of small RNAs pool and this pool can produce many BmAGO2-associated small RNAs that regulate TE’s expression.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3598-5) contains supplementary material, which is available to authorized users.

Telomeric Retrotransposon HeT-A Contains a Bidirectional Promoter that Initiates Divergent Transcription of piRNA Precursors in Drosophila Germline

PIWI-interacting RNAs (piRNAs) provide the silencing of transposable elements in the germline. Drosophila telomeres are maintained by transpositions of specialized telomeric retroelements. piRNAs generated from sense and antisense transcripts of telomeric elements provide telomere length control in the germline. Previously, we have found that antisense transcription of the major telomeric retroelement HeT-A is initiated upstream of the HeT-A sense transcription start site. Here, we performed a deletion analysis of the HeT-A promoter and show that common regulatory elements are shared by sense and antisense promoters of HeT-A. Therefore, the HeT-A promoter is a bidirectional promoter capable of processive sense and antisense transcription. Ovarian small RNA data show that a solo HeT-A promoter within an euchromatic transgene initiates the divergent transcription of transgenic reporter genes and subsequent processing of these transcripts into piRNAs. These events lead to the formation of a divergent unistrand piRNA cluster at solo HeT-A promoters, in contrast to endogenous telomeres that represent strong dual-strand piRNA clusters. Solo HeT-A promoters are not immunoprecipitated with heterochromatin protein 1 (HP1) homolog Rhino, a marker of the dual-strand piRNA clusters, but are associated with HP1 itself, which provides piRNA-mediated transcriptional repression of the reporter genes. Unlike endogenous dual-strand piRNA clusters, the solo HeT-A promoter does not produce overlapping transcripts. In a telomeric context, however, bidirectional promoters of tandem HeT-A repeats provide a read-through transcription of both genomic strands, followed by Rhi binding. These data indicate that Drosophila telomeres share properties of unistrand and dual-strand piRNA clusters.

What Drives Positive Selection in the Drosophila piRNA Machinery? The Genomic Autoimmunity Hypothesis

In animals, PIWI-interacting RNAs (piRNAs) play a crucial role in genome defense. Moreover, because piRNAs can be maternally transmitted, they contribute to the epigenetic profile of inheritance. Multiple studies, especially in Drosophila, have demonstrated that the machinery of piRNA biogenesis is often the target of positive selection. Because transposable elements (TEs) are a form of genetic parasite, positive selection in the piRNA machinery is often explained by analogy to the signatures of positive selection commonly observed in genes that play a role in host-parasite dynamics. However, the precise mechanisms that drive positive selection in the piRNA machinery are not known. In this review, we outline several mechanistic models that might explain pervasive positive selection in the piRNA machinery of Drosophila species. We propose that recurrent positive selection in the piRNA machinery can be partly explained by an ongoing tension between selection for sensitivity required by genome defense and selection for specificity to avoid the off-target effects of maladaptive genic silencing by piRNA.

Unique transposon landscapes are pervasive across Drosophila melanogaster genomes

To understand how transposon landscapes (TLs) vary across animal genomes, we describe a new method called the Transposon Insertion and Depletion AnaLyzer (TIDAL) and a database of >300 TLs in Drosophila melanogaster (TIDAL-Fly). Our analysis reveals pervasive TL diversity across cell lines and fly strains, even for identically named sub-strains from different laboratories such as the ISO1 strain used for the reference genome sequence. On average, >500 novel insertions exist in every lab strain, inbred strains of the Drosophila Genetic Reference Panel (DGRP), and fly isolates in the Drosophila Genome Nexus (DGN). A minority (<25%) of transposon families comprise the majority (>70%) of TL diversity across fly strains. A sharp contrast between insertion and depletion patterns indicates that many transposons are unique to the ISO1 reference genome sequence. Although TL diversity from fly strains reaches asymptotic limits with increasing sequencing depth, rampant TL diversity causes unsaturated detection of TLs in pools of flies. Finally, we show novel transposon insertions negatively correlate with Piwi-interacting RNA (piRNA) levels for most transposon families, except for the highly-abundant roo retrotransposon. Our study provides a useful resource for Drosophila geneticists to understand how transposons create extensive genomic diversity in fly cell lines and strains.

P elements and the determinants of hybrid dysgenesis have different dynamics of propagation in Drosophila melanogaster populations

Intraspecific hybrid dysgenesis (HD) appears after some strains of D. melanogaster are crossed. The predominant idea is that the movement of transposable P elements causes HD. It is believed that P elements appeared in the D. melanogaster genome in the middle of the last century by horizontal transfer, simultaneously with the appearance of HD determinants. A subsequent simultaneous expansion of HD determinants and P elements occurred. We analyzed the current distribution of HD determinants in natural populations of D. melanogaster and found no evidence of their further spread. However, full-sized P elements were identified in the genomes of all analyzed natural D. melanogaster strains independent of their cytotypes. Thus, the expansion of P elements does not correlate with the expansion of HD determinants. We found that the ovaries of dysgenic females did not contain germ cells despite the equal number of primordial germ cells in early stages in dysgenic and non-dysgenic embryos. We propose that HD does not result from DNA damage caused by P element transposition, but it would be the disruption in the regulation of dysgenic ovarian formation that causes the dysgenic phenotypes.

When natural selection gives gene function the cold shoulder

It is tempting to invoke organismal selection as perpetually optimizing the function of any given gene. However, natural selection can drive genic functional change without improvement of biochemical activity, even to the extinction of gene activity. Detrimental mutations can creep in owing to linkage with other selectively favored loci. Selection can promote functional degradation, irrespective of genetic drift, when adaptation occurs by loss of gene function. Even stabilizing selection on a trait can lead to divergence of the underlying molecular constituents. Selfish genetic elements can also proliferate independent of any functional benefits to the host genome. Here we review the logic and evidence for these diverse processes acting in genome evolution. This collection of distinct evolutionary phenomena - while operating through easily understandable mechanisms - all contribute to the seemingly counterintuitive notion that maintenance or improvement of a gene's biochemical function sometimes do not determine its evolutionary fate.