The transposable elements of the Drosophila melanogaster euchromatin: a genomics perspective

Expansion of three types of transposon superfamilies within 25 Mya lead to large genome size of a rice insect pest

The brown planthoppers (BPH, Nilaparvata lugens), white backed planthopper (WBPH, Sogatella furcifera) and small brown planthopper (SBPH, Laodelphax striatellus) are widely distributed rice insect pests, causing huge annual yield loss of rice production. Though these three planthoppers belong to the same family, Delphacidae of Hemiptera, their genome sizes (GS) are very different, ranging from 541 to 1088 Mb. To uncover the main factors driving GS changes of three planthoppers, we first estimated the GS of their ancestor Fulgoroidea, to be 794.33 Mb, indicating GS expansion in BPH but contraction in SBPH and WBPH. Next, we identified repetitive sequences and compared the TE landscapes, showed that three types of transposon superfamilies, hAT, Tc1-Mariner and Gypsy, expanded within 25 Mya in BPH. In addition, BPH kept ancient TEs of Fulgoroidea dated back to 175 Mya, while SBPH and WBPH have lost most of these ancient TEs. Here, we present evidence that the gain of recently expanded TEs driving the GS expansion and loss of ancient TEs leading to the GS contraction, providing new insights into the mechanism of GS variation.

MCHelper automatically curates transposable element libraries across eukaryotic species

The number of species with high-quality genome sequences continues to increase, in part due to the scaling up of multiple large-scale biodiversity sequencing projects. While the need to annotate genic sequences in these genomes is widely acknowledged, the parallel need to annotate transposable element (TE) sequences that have been shown to alter genome architecture, rewire gene regulatory networks, and contribute to the evolution of host traits is becoming ever more evident. However, accurate genome-wide annotation of TE sequences is still technically challenging. Several de novo TE identification tools are now available, but manual curation of the libraries produced by these tools is needed to generate high-quality genome annotations. Manual curation is time-consuming, and thus impractical for large-scale genomic studies, and lacks reproducibility. In this work, we present the Manual Curator Helper tool MCHelper, which automates the TE library curation process. By leveraging MCHelper's fully automated mode with the outputs from three de novo TE identification tools, RepeatModeler2, EDTA, and REPET, in the fruit fly, rice, hooded crow, zebrafish, maize, and human, we show a substantial improvement in the quality of the TE libraries and genome annotations. MCHelper libraries are less redundant, with up to 65% reduction in the number of consensus sequences, have up to 11.4% fewer false positive sequences, and up to ∼48% fewer "unclassified/unknown" TE consensus sequences. Genome-wide TE annotations are also improved, including larger unfragmented insertions. Moreover, MCHelper is an easy-to-install and easy-to-use tool.

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.

Chromosome-scale genome of the polyphagous pest Anastrepha ludens (Diptera: Tephritidae) provides insights on sex chromosome evolution in Anastrepha

The Mexican fruit fly, Anastrepha ludens, is a polyphagous true fruit fly (Diptera: Tephritidae) considered one of the most serious insect pests in Central and North America to various economically relevant fruits. Despite its agricultural relevance, a high-quality genome assembly has not been reported. Here, we described the generation of a chromosome-level genome for the A. ludens using a combination of PacBio high fidelity long-reads and chromatin conformation capture sequencing data. The final assembly consisted of 140 scaffolds (821 Mb, N50 = 131 Mb), containing 99.27% complete conserved orthologs (BUSCO) for Diptera. We identified the sex chromosomes using three strategies: 1) visual inspection of Hi-C contact map and coverage analysis using the HiFi reads, 2) synteny with Drosophila melanogaster, and 3) the difference in the average read depth of autosomal versus sex chromosomal scaffolds. The X chromosome was found in one major scaffold (100 Mb) and eight smaller contigs (1.8 Mb), and the Y chromosome was recovered in one large scaffold (6.1 Mb) and 35 smaller contigs (4.3 Mb). Sex chromosomes and autosomes showed considerable differences of transposable elements and gene content. Moreover, evolutionary rates of orthologs of A. ludens and Anastrepha obliqua revealed a faster evolution of X-linked, compared to autosome-linked, genes, consistent with the faster-X effect, leading us to new insights on the evolution of sex chromosomes in this diverse group of flies. This genome assembly provides a valuable resource for future evolutionary, genetic, and genomic translational research supporting the management of this important agricultural pest.

Whole-Genome and Poly(A)+Transcriptome Analysis of the Drosophila Mutant agnts3 with Cognitive Dysfunctions

The temperature-sensitive Drosophila mutant agnts3 exhibits the restoration of learning defects both after heat shock (HS) and under hypomagnetic conditions (HMC). Previously, agnts3 was shown to have an increased level of LIM kinase 1 (LIMK1). However, its limk1 sequence did not significantly differ from that of the wild-type strain Canton-S (CS). Here, we performed whole-genome and poly(A)-enriched transcriptome sequencing of CS and agnts3 males normally, after HMC, and after HS. Several high-effect agnts3-specific mutations were identified, including MED23 (regulation of HS-dependent transcription) and Spn42De, the human orthologs of which are associated with intellectual disorders. Pronounced interstrain differences between the transcription profiles were revealed. Mainly, they included the genes of defense and stress response, long non-coding RNAs, and transposons. After HS, the differences between the transcriptomes became less pronounced. In agnts3, prosalpha1 was the only gene whose expression changed after both HS and HMC. The normal downregulation of prosalpha1 and Spn42De in agnts3 was confirmed by RT-PCR. Analysis of limk1 expression did not reveal any interstrain differences or changes after stress. Thus, behavioral differences between CS and agnts3 both under normal and stressed conditions are not due to differences in limk1 transcription. Instead, MED23, Spn42De, and prosalpha1 are more likely to contribute to the agnts3 phenotype.

Genomic context-dependent histone H3K36 methylation by three Drosophila methyltransferases and implications for dedicated chromatin readers

Methylation of histone H3 at lysine 36 (H3K36me3) marks active chromatin. The mark is interpreted by epigenetic readers that assist transcription and safeguard the integrity of the chromatin fiber. The chromodomain protein MSL3 binds H3K36me3 to target X-chromosomal genes in male Drosophila for dosage compensation. The PWWP-domain protein JASPer recruits the JIL1 kinase to active chromatin on all chromosomes. Unexpectedly, depletion of K36me3 had variable, locus-specific effects on the interactions of those readers. This observation motivated a systematic and comprehensive study of K36 methylation in a defined cellular model. Contrasting prevailing models, we found that K36me1, K36me2 and K36me3 each contribute to distinct chromatin states. A gene-centric view of the changing K36 methylation landscape upon depletion of the three methyltransferases Set2, NSD and Ash1 revealed local, context-specific methylation signatures. Set2 catalyzes K36me3 predominantly at transcriptionally active euchromatin. NSD places K36me2/3 at defined loci within pericentric heterochromatin and on weakly transcribed euchromatic genes. Ash1 deposits K36me1 at regions with enhancer signatures. The genome-wide mapping of MSL3 and JASPer suggested that they bind K36me2 in addition to K36me3, which was confirmed by direct affinity measurement. This dual specificity attracts the readers to a broader range of chromosomal locations and increases the robustness of their actions.

Rapid evolution of piRNA clusters in the Drosophila melanogaster ovary

The piRNA pathway is a highly conserved mechanism to repress transposable element (TE) activity in the animal germline via a specialized class of small RNAs called piwi-interacting RNAs (piRNAs). piRNAs are produced from discrete genomic regions called piRNA clusters (piCs). Although the molecular processes by which piCs function are relatively well understood in Drosophila melanogaster, much less is known about the origin and evolution of piCs in this or any other species. To investigate piC origin and evolution, we use a population genomic approach to compare piC activity and sequence composition across eight geographically distant strains of D. melanogaster with high-quality long-read genome assemblies. We perform annotations of ovary piCs and genome-wide TE content in each strain. Our analysis uncovers extensive variation in piC activity across strains and signatures of rapid birth and death of piCs. Most TEs inferred to be recently active show an enrichment of insertions into old and large piCs, consistent with the previously proposed "trap" model of piC evolution. In contrast, a small subset of active LTR families is enriched for the formation of new piCs, suggesting that these TEs have higher proclivity to form piCs. Thus, our findings uncover processes leading to the origin of piCs. We propose that piC evolution begins with the emergence of piRNAs from individual insertions of a few select TE families prone to seed new piCs that subsequently expand by accretion of insertions from most other TE families during evolution to form larger "trap" clusters. Our study shows that TEs themselves are the major force driving the rapid evolution of piCs.

The impact of differential transposition activities of autonomous and nonautonomous hAT transposable elements on genome architecture and gene expression in Caenorhabditis inopinata

Transposable elements are DNA sequences capable of moving within genomes and significantly influence genomic evolution. The nematode Caenorhabditis inopinata exhibits a much higher transposable element copy number than its sister species, Caenorhabditis elegans. In this study, we identified a novel autonomous transposable element belonging to the hAT superfamily from a spontaneous transposable element-insertion mutant in C. inopinata and named this transposon Ci-hAT1. Further bioinformatic analyses uncovered 3 additional autonomous hAT elements-Ci-hAT2, Ci-hAT3, and Ci-hAT4-along with over 1,000 copies of 2 nonautonomous miniature inverted-repeat transposable elements, mCi-hAT1 and mCi-hAT4, likely derived from Ci-hAT1 and Ci-hAT4 through internal deletion. We tracked at least 3 sequential transpositions of Ci-hAT1 over several years. However, the transposition rates of the other 3 autonomous hAT elements were lower, suggesting varying activity levels. Notably, the distribution patterns of the 2 miniature inverted-repeat transposable element families differed significantly: mCi-hAT1 was primarily located in the chromosome arms, a pattern observed in the transposable elements of other Caenorhabditis species, whereas mCi-hAT4 was more evenly distributed across chromosomes. Additionally, interspecific transcriptome analysis indicated that C. inopinata genes with upstream or intronic these miniature inverted-repeat transposable element insertions tend to be more highly expressed than their orthologous genes in C. elegans. These findings highlight the significant role of de-silenced transposable elements in driving the evolution of genomes and transcriptomes, leading to species-specific genetic diversity.

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.

Piwi expressed in Drosophila adipose tissues regulates systemic IGF signaling and growth via IGF-binding protein

Piwi and its partner, Piwi-interacting RNA (piRNA), are pivotal in suppressing the harmful effects of transposable elements (TEs) linked to genomic insertional mutagenesis. While primarily active in Drosophila's adult gonadal tissues, causing sterility in its absence, Piwi's role in post-embryonic development remains unclear. Our study reveals Piwi's functional presence in the larval fat body, where it governs developmental growth through systemic insulin/insulin-like growth factor (IGF) signaling (IIS). Piwi knockdown in the fat body resulted in dysregulated TE expression, reduced developmental rate and body growth, and diminished systemic IIS activity. Notably, Piwi knockdown increased Imaginal Morphogenic Protein Late 2 (Imp-L2) expression, akin to insulin-like growth factor-binding protein 7 (IGFBP7), reducing systemic IIS and inhibiting body growth. This unveils a novel role for Piwi in larval adipose tissues, emphasizing its importance in regulating systemic IIS and overall organismal growth.

ChimeraTE: a pipeline to detect chimeric transcripts derived from genes and transposable elements

Transposable elements (TEs) produce structural variants and are considered an important source of genetic diversity. Notably, TE-gene fusion transcripts, i.e. chimeric transcripts, have been associated with adaptation in several species. However, the identification of these chimeras remains hindered due to the lack of detection tools at a transcriptome-wide scale, and to the reliance on a reference genome, even though different individuals/cells/strains have different TE insertions. Therefore, we developed ChimeraTE, a pipeline that uses paired-end RNA-seq reads to identify chimeric transcripts through two different modes. Mode 1 is the reference-guided approach that employs canonical genome alignment, and Mode 2 identifies chimeras derived from fixed or insertionally polymorphic TEs without any reference genome. We have validated both modes using RNA-seq data from four Drosophila melanogaster wild-type strains. We found ∼1.12% of all genes generating chimeric transcripts, most of them from TE-exonized sequences. Approximately ∼23% of all detected chimeras were absent from the reference genome, indicating that TEs belonging to chimeric transcripts may be recent, polymorphic insertions. ChimeraTE is the first pipeline able to automatically uncover chimeric transcripts without a reference genome, consisting of two running Modes that can be used as a tool to investigate the contribution of TEs to transcriptome plasticity.

A genetic basis for facultative parthenogenesis in Drosophila

Facultative parthenogenesis enables sexually reproducing organisms to switch between sexual and asexual parthenogenetic reproduction. To gain insights into this phenomenon, we sequenced the genomes of sexually reproducing and parthenogenetic strains of Drosophila mercatorum and identified differences in the gene expression in their eggs. We then tested whether manipulating the expression of candidate gene homologs identified in Drosophila mercatorum could lead to facultative parthenogenesis in the non-parthenogenetic species Drosophila melanogaster. This identified a polygenic system whereby increased expression of the mitotic protein kinase polo and decreased expression of a desaturase, Desat2, caused facultative parthenogenesis in the non-parthenogenetic species that was enhanced by increased expression of Myc. The genetically induced parthenogenetic Drosophila melanogaster eggs exhibit de novo centrosome formation, fusion of the meiotic products, and the onset of development to generate predominantly triploid offspring. Thus, we demonstrate a genetic basis for sporadic facultative parthenogenesis in an animal.

Generation of LexA enhancer-trap lines in Drosophila by an international scholastic network

Conditional gene regulation in Drosophila through binary expression systems like the LexA-LexAop system provides a superb tool for investigating gene and tissue function. To increase the availability of defined LexA enhancer trap insertions, we present molecular, genetic, and tissue expression studies of 301 novel Stan-X LexA enhancer traps derived from mobilization of the index SX4 line. This includes insertions into distinct loci on the X, II, and III chromosomes that were not previously associated with enhancer traps or targeted LexA constructs, an insertion into ptc, and seventeen insertions into natural transposons. A subset of enhancer traps was expressed in CNS neurons known to produce and secrete insulin, an essential regulator of growth, development, and metabolism. Fly lines described here were generated and characterized through studies by students and teachers in an international network of genetics classes at public, independent high schools, and universities serving a diversity of students, including those underrepresented in science. Thus, a unique partnership between secondary schools and university-based programs has produced and characterized novel resources in Drosophila, establishing instructional paradigms devoted to unscripted experimental science.

Third-Generation Sequencing Reveals the Adaptive Role of the Epigenome in Three Deep-Sea Polychaetes

The roles of DNA methylation in invertebrates are poorly characterized, and critical data are missing for the phylum Annelida. We fill this knowledge gap by conducting the first genome-wide survey of DNA methylation in the deep-sea polychaetes dominant in deep-sea vents and seeps: Paraescarpia echinospica, Ridgeia piscesae, and Paralvinella palmiformis. DNA methylation calls were inferred from Oxford Nanopore sequencing after assembling high-quality genomes of these animals. The genomes of these worms encode all the key enzymes of the DNA methylation metabolism and possess a mosaic methylome similar to that of other invertebrates. Transcriptomic data of these polychaetes support the hypotheses that gene body methylation strengthens the expression of housekeeping genes and that promoter methylation acts as a silencing mechanism but not the hypothesis that DNA methylation suppresses the activity of transposable elements. The conserved epigenetic profiles of genes responsible for maintaining homeostasis under extreme hydrostatic pressure suggest DNA methylation plays an important adaptive role in these worms.

Retrotransposons hijack alt-EJ for DNA replication and eccDNA biogenesis

Retrotransposons are highly enriched in the animal genome1-3. The activation of retrotransposons can rewrite host DNA information and fundamentally impact host biology1-3. Although developmental activation of retrotransposons can offer benefits for the host, such as against virus infection, uncontrolled activation promotes disease or potentially drives ageing1-5. After activation, retrotransposons use their mRNA as templates to synthesize double-stranded DNA for making new insertions in the host genome1-3,6. Although the reverse transcriptase that they encode can synthesize the first-strand DNA1-3,6, how the second-strand DNA is generated remains largely unclear. Here we report that retrotransposons hijack the alternative end-joining (alt-EJ) DNA repair process of the host for a circularization step to synthesize their second-strand DNA. We used Nanopore sequencing to examine the fates of replicated retrotransposon DNA, and found that 10% of them achieve new insertions, whereas 90% exist as extrachromosomal circular DNA (eccDNA). Using eccDNA production as a readout, further genetic screens identified factors from alt-EJ as essential for retrotransposon replication. alt-EJ drives the second-strand synthesis of the long terminal repeat retrotransposon DNA through a circularization process and is therefore necessary for eccDNA production and new insertions. Together, our study reveals that alt-EJ is essential in driving the propagation of parasitic genomic retroelements. Our study uncovers a conserved function of this understudied DNA repair process, and provides a new perspective to understand-and potentially control-the retrotransposon life cycle.

Dicentric chromosome breakage in Drosophila melanogaster is influenced by pericentric heterochromatin and occurs in nonconserved hotspots

Chromosome breakage plays an important role in the evolution of karyotypes and can produce deleterious effects within a single individual, such as aneuploidy or cancer. Forces that influence how and where chromosomes break are not fully understood. In humans, breakage tends to occur in conserved hotspots called common fragile sites (CFS), especially during replication stress. By following the fate of dicentric chromosomes in Drosophila melanogaster, we find that breakage under tension also tends to occur in specific hotspots. Our experimental approach was to induce sister chromatid exchange in a ring chromosome to generate a dicentric chromosome with a double chromatid bridge. In the following cell division, the dicentric bridges may break. We analyzed the breakage patterns of 3 different ring-X chromosomes. These chromosomes differ by the amount and quality of heterochromatin they carry as well as their genealogical history. For all 3 chromosomes, breakage occurs preferentially in several hotspots. Surprisingly, we found that the hotspot locations are not conserved between the 3 chromosomes: each displays a unique array of breakage hotspots. The lack of hotspot conservation, along with a lack of response to aphidicolin, suggests that these breakage sites are not entirely analogous to CFS and may reveal new mechanisms of chromosome fragility. Additionally, the frequency of dicentric breakage and the durability of each chromosome's spindle attachment vary significantly between the 3 chromosomes and are correlated with the origin of the centromere and the amount of pericentric heterochromatin. We suggest that different centromere strengths could account for this.

The Landscape of the DNA Transposons in the Genome of the Horezu_LaPeri Strain of Drosophila melanogaster

Natural transposons (NTs) represent mobile DNA sequences found in both prokaryotic and eukaryotic genomes. Drosophila melanogaster (the fruit fly) is a eukaryotic model organism with NTs standing for about 20% of its genome and has contributed significantly to the understanding of various aspects of transposon biology. Our study describes an accurate approach designed to map class II transposons (DNA transposons) in the genome of the Horezu_LaPeri fruit fly strain, consecutive to Oxford Nanopore Technology sequencing. A whole genome bioinformatics analysis was conducted using Genome ARTIST_v2, LoRTE and RepeatMasker tools to identify DNA transposons insertions. Then, a gene ontology enrichment analysis was performed in order to evaluate the potential adaptive role of some DNA transposons insertions. Herein, we describe DNA transposon insertions specific for the Horezu_LaPeri genome and a predictive functional analysis of some insertional alleles. The PCR validation of P-element insertions specific for this fruit fly strain, along with a putative consensus sequence for the KP element, is also reported. Overall, the genome of the Horezu_LaPeri strain contains several insertions of DNA transposons associated with genes known to be involved in adaptive processes. For some of these genes, insertional alleles obtained via mobilization of the artificial transposons were previously reported. This is a very alluring aspect, as it suggests that insertional mutagenesis experiments conducting adaptive predictions for laboratory strains may be confirmed by mirroring insertions which are expected to be found at least in some natural fruit fly strains.

Rapid evolution of piRNA clusters in the Drosophila melanogaster ovary

Animal genomes are parasitized by a horde of transposable elements (TEs) whose mutagenic activity can have catastrophic consequences. The piRNA pathway is a conserved mechanism to repress TE activity in the germline via a specialized class of small RNAs associated with effector Piwi proteins called piwi-associated RNAs (piRNAs). piRNAs are produced from discrete genomic regions called piRNA clusters (piCs). While piCs are generally enriched for TE sequences and the molecular processes by which they are transcribed and regulated are relatively well understood in Drosophila melanogaster, much less is known about the origin and evolution of piCs in this or any other species. To investigate piC evolution, we use a population genomics approach to compare piC activity and sequence composition across 8 geographically distant strains of D. melanogaster with high quality long-read genome assemblies. We perform extensive annotations of ovary piCs and TE content in each strain and test predictions of two proposed models of piC evolution. The 'de novo' model posits that individual TE insertions can spontaneously attain the status of a small piC to generate piRNAs silencing the entire TE family. The 'trap' model envisions large and evolutionary stable genomic clusters where TEs tend to accumulate and serves as a long-term "memory" of ancient TE invasions and produce a great variety of piRNAs protecting against related TEs entering the genome. It remains unclear which model best describes the evolution of piCs. Our analysis uncovers extensive variation in piC activity across strains and signatures of rapid birth and death of piCs in natural populations. Most TE families inferred to be recently or currently active show an enrichment of strain-specific insertions into large piCs, consistent with the trap model. By contrast, only a small subset of active LTR retrotransposon families is enriched for the formation of strain-specific piCs, suggesting that these families have an inherent proclivity to form de novo piCs. Thus, our findings support aspects of both 'de novo' and 'trap' models of piC evolution. We propose that these two models represent two extreme stages along an evolutionary continuum, which begins with the emergence of piCs de novo from a few specific LTR retrotransposon insertions that subsequently expand by accretion of other TE insertions during evolution to form larger 'trap' clusters. Our study shows that piCs are evolutionarily labile and that TEs themselves are the major force driving the formation and evolution of piCs.

Epigenetic Silencing of P-Element Reporter Genes Induced by Transcriptionally Active Domains of Constitutive Heterochromatin in Drosophila melanogaster

Reporter genes inserted via P-element integration into different locations of the Drosophila melanogaster genome have been routinely used to monitor the functional state of chromatin domains. It is commonly thought that P-element-derived reporter genes are subjected to position effect variegation (PEV) when transposed into constitutive heterochromatin because they acquire heterochromatin-like epigenetic modifications that promote silencing. However, sequencing and annotation of the D. melanogaster genome have shown that constitutive heterochromatin is a genetically and molecularly heterogeneous compartment. In fact, in addition to repetitive DNAs, it harbors hundreds of functional genes, together accounting for a significant fraction of its entire genomic territory. Notably, most of these genes are actively transcribed in different developmental stages and tissues, irrespective of their location in heterochromatin. An open question in the genetic and molecular studies on PEV in D. melanogaster is whether functional heterochromatin domains, i.e., heterochromatin harboring active genes, are able to silence reporter genes therein transposed or, on the contrary, can drive their expression. In this work, we provide experimental evidence showing that strong silencing of the Pw⁺ reporters is induced even when they are integrated within or near actively transcribed loci in the pericentric regions of chromosome 2. Interestingly, some Pw⁺ reporters were found insensitive to the action of a known PEV suppressor. Two of them are inserted within Yeti, a gene expressed in the deep heterochromatin of chromosome 2 which carries active chromatin marks. The difference sensitivity to suppressors-exhibited Pw⁺ reporters supports the view that different epigenetic regulators or mechanisms control different regions of heterochromatin. Together, our results suggest that there may be more complexity regarding the molecular mechanisms underlying PEV.

Retrotransposon activation during Drosophila metamorphosis conditions adult antiviral responses

Retrotransposons are one type of mobile genetic element that abundantly reside in the genomes of nearly all animals. Their uncontrolled activation is linked to sterility, cancer and other pathologies, thereby being largely considered detrimental. Here we report that, within a specific time window of development, retrotransposon activation can license the host's immune system for future antiviral responses. We found that the mdg4 (also known as Gypsy) retrotransposon selectively becomes active during metamorphosis at the Drosophila pupal stage. At this stage, mdg4 activation educates the host's innate immune system by inducing the systemic antiviral function of the nuclear factor-κB protein Relish in a dSTING-dependent manner. Consequently, adult flies with mdg4, Relish or dSTING silenced at the pupal stage are unable to clear exogenous viruses and succumb to viral infection. Altogether, our data reveal that hosts can establish a protective antiviral response that endows a long-term benefit in pathogen warfare due to the developmental activation of mobile genetic elements.

Genetic variation in P-element dysgenic sterility is associated with double-strand break repair and alternative splicing of TE transcripts

The germline mobilization of transposable elements (TEs) by small RNA mediated silencing pathways is conserved across eukaryotes and critical for ensuring the integrity of gamete genomes. However, genomes are recurrently invaded by novel TEs through horizontal transfer. These invading TEs are not targeted by host small RNAs, and their unregulated activity can cause DNA damage in germline cells and ultimately lead to sterility. Here we use hybrid dysgenesis-a sterility syndrome of Drosophila caused by transposition of invading P-element DNA transposons-to uncover host genetic variants that modulate dysgenic sterility. Using a panel of highly recombinant inbred lines of Drosophila melanogaster, we identified two linked quantitative trait loci (QTL) that determine the severity of dysgenic sterility in young and old females, respectively. We show that ovaries of fertile genotypes exhibit increased expression of splicing factors that suppress the production of transposase encoding transcripts, which likely reduces the transposition rate and associated DNA damage. We also show that fertile alleles are associated with decreased sensitivity to double-stranded breaks and enhanced DNA repair, explaining their ability to withstand high germline transposition rates. Together, our work reveals a diversity of mechanisms whereby host genotype modulates the cost of an invading TE, and points to genetic variants that were likely beneficial during the P-element invasion.

ONT-Based Alternative Assemblies Impact on the Annotations of Unique versus Repetitive Features in the Genome of a Romanian Strain of Drosophila melanogaster

To date, different strategies of whole-genome sequencing (WGS) have been developed in order to understand the genome structure and functions. However, the analysis of genomic sequences obtained from natural populations is challenging and the biological interpretation of sequencing data remains the main issue. The MinION device developed by Oxford Nanopore Technologies (ONT) is able to generate long reads with minimal costs and time requirements. These valuable assets qualify it as a suitable method for performing WGS, especially in small laboratories. The long reads resulted using this sequencing approach can cover large structural variants and repetitive sequences commonly present in the genomes of eukaryotes. Using MinION, we performed two WGS assessments of a Romanian local strain of Drosophila melanogaster, referred to as Horezu_LaPeri (Horezu). In total, 1,317,857 reads with a size of 8.9 gigabytes (Gb) were generated. Canu and Flye de novo assembly tools were employed to obtain four distinct assemblies with both unfiltered and filtered reads, achieving maximum reference genome coverages of 94.8% (Canu) and 91.4% (Flye). In order to test the quality of these assemblies, we performed a two-step evaluation. Firstly, we considered the BUSCO scores and inquired for a supplemental set of genes using BLAST. Subsequently, we appraised the total content of natural transposons (NTs) relative to the reference genome (ISO1 strain) and mapped the mdg1 retroelement as a resolution assayer. Our results reveal that filtered data provide only slightly enhanced results when considering genes identification, but the use of unfiltered data had a consistent positive impact on the global evaluation of the NTs content. Our comparative studies also revealed differences between Flye and Canu assemblies regarding the annotation of unique versus repetitive genomic features. In our hands, Flye proved to be moderately better for gene identification, while Canu clearly outperformed Flye for NTs analysis. Data concerning the NTs content were compared to those obtained with ONT for the D. melanogaster ISO1 strain, revealing that our strategy conducted to better results. Additionally, the parameters of our ONT reads and assemblies are similar to those reported for ONT experiments performed on various model organisms, revealing that our assembly data are appropriate for a proficient annotation of the Horezu genome.

Transposable element finder (TEF): finding active transposable elements from next generation sequencing data

BACKGROUND

Detection of newly transposed events by transposable elements (TEs) from next generation sequence (NGS) data is difficult, due to their multiple distribution sites over the genome containing older TEs. The previously reported Transposon Insertion Finder (TIF) detects TE transpositions on the reference genome from NGS short reads using end sequences of target TE. TIF requires the sequence of target TE and is not able to detect transpositions for TEs with an unknown sequence.

RESULT

The new algorithm Transposable Element Finder (TEF) enables the detection of TE transpositions, even for TEs with an unknown sequence. TEF is a finding tool of transposed TEs, in contrast to TIF as a detection tool of transposed sites for TEs with a known sequence. The transposition event is often accompanied with a target site duplication (TSD). Focusing on TSD, two algorithms to detect both ends of TE, TSDs and target sites are reported here. One is based on the grouping with TSDs and direct comparison of k-mers from NGS without similarity search. The other is based on the junction mapping of TE end sequence candidates. Both methods succeed to detect both ends and TSDs of known active TEs in several tests with rice, Arabidopsis and Drosophila data and discover several new TEs in new locations. PCR confirmed the detected transpositions of TEs in several test cases in rice.

CONCLUSIONS

TEF detects transposed TEs with TSDs as a result of TE transposition, sequences of both ends and their inserted positions of transposed TEs by direct comparison of NGS data between two samples. Genotypes of transpositions are verified by counting of junctions of head and tail, and non-insertion sequences in NGS reads. TEF is easy to run and independent of any TE library, which makes it useful to detect insertions from unknown TEs bypassed by common TE annotation pipelines.

Lysine 27 dimethylation of Drosophila linker histone dH1 contributes to heterochromatin organization independently of H3K9 methylation

Post-translational modifications (PTMs) of core histones are important epigenetic determinants that correlate with functional chromatin states. However, despite multiple linker histone H1s PTMs have been identified, little is known about their genomic distribution and contribution to the epigenetic regulation of chromatin. Here, we address this question in Drosophila that encodes a single somatic linker histone, dH1. We previously reported that dH1 is dimethylated at K27 (dH1K27me2). Here, we show that dH1K27me2 is a major PTM of Drosophila heterochromatin. At mitosis, dH1K27me2 accumulates at pericentromeric heterochromatin, while, in interphase, it is also detected at intercalary heterochromatin. ChIPseq experiments show that >98% of dH1K27me2 enriched regions map to heterochromatic repetitive DNA elements, including transposable elements, simple DNA repeats and satellite DNAs. Moreover, expression of a mutated dH1K27A form, which impairs dH1K27me2, alters heterochromatin organization, upregulates expression of heterochromatic transposable elements and results in the accumulation of RNA:DNA hybrids (R-loops) in heterochromatin, without affecting H3K9 methylation and HP1a binding. The pattern of dH1K27me2 is H3K9 methylation independent, as it is equally detected in flies carrying a H3K9R mutation, and is not affected by depletion of Su(var)3-9, HP1a or Su(var)4-20. Altogether these results suggest that dH1K27me2 contributes to heterochromatin organization independently of H3K9 methylation.

Species-specific chromatin landscape determines how transposable elements shape genome evolution

Transposable elements (TEs) are selfish genetic parasites that increase their copy number at the expense of host fitness. The ‘success’, or genome-wide abundance, of TEs differs widely between species. Deciphering the causes for this large variety in TE abundance has remained a central question in evolutionary genomics. We previously proposed that species-specific TE abundance could be driven by the inadvertent consequences of host-direct epigenetic silencing of TEs—the spreading of repressive epigenetic marks from silenced TEs into adjacent sequences. Here, we compared this TE-mediated local enrichment of repressive marks, or ‘the epigenetic effect of TEs’, in six species in the Drosophila melanogaster subgroup to dissect step-by-step the role of such effect in determining genomic TE abundance. We found that TE-mediated local enrichment of repressive marks is prevalent and substantially varies across and even within species. While this TE-mediated effect alters the epigenetic states of adjacent genes, we surprisingly discovered that the transcription of neighboring genes could reciprocally impact this spreading. Importantly, our multi-species analysis provides the power and appropriate phylogenetic resolution to connect species-specific host chromatin regulation, TE-mediated epigenetic effects, the strength of natural selection against TEs, and genomic TE abundance unique to individual species. Our findings point toward the importance of host chromatin landscapes in shaping genome evolution through the epigenetic effects of a selfish genetic parasite.

All the instructions required for life are encoded in the set of DNA present in a cell. It therefore seems natural to think that every bit of this genetic information should serve the organism. And yet most species carry parasitic ‘transposable’ sequences, or transposons, whose only purpose is to multiply and insert themselves at other positions in the genome.

It is possible for cells to suppress these selfish elements. Chemical marks can be deposited onto the DNA to temporarily ‘silence’ transposons and prevent them from being able to move and replicate. However, this sometimes comes at a cost: the repressive chemical modifications can spread to nearby genes that are essential for the organism and perturb their function.

Strangely, the prevalence of transposons varies widely across the tree of life. These sequences form the majority of the genome of certain species – in fact, they represent about half of the human genetic information. But their abundance is much lower in other organisms, forming a measly 6% of the genome of puffer fish for instance. Even amongst fruit fly species, the prevalence of transposable elements can range between 2% and 25%. What explains such differences?

Huang et al. set out to examine this question through the lens of transposon silencing, systematically comparing how this process impacts nearby regions in six species of fruit flies. This revealed variations in the strength of the side effects associated with transposon silencing, resulting in different levels of perturbation on neighbouring genes. A stronger impact was associated with the species having fewer transposons in its genome, suggesting that an evolutionary pressure is at work to keep the abundance of transposons at a low level in these species. Further analyses showed that the genes which determine how silencing marks are distributed may also be responsible for the variations in the impact of transposon silencing. They could therefore be the ones driving differences in the abundance of transposons between species.

Overall, this work sheds light on the complex mechanisms shaping the evolution of genomes, and it may help to better understand how transposons are linked to processes such as aging and cancer.

Artificially stimulating retrotransposon activity increases mortality and accelerates a subset of aging phenotypes in Drosophila

Transposable elements (TEs) are mobile sequences of DNA that can become transcriptionally active as an animal ages. Whether TE activity is simply a by-product of heterochromatin breakdown or can contribute toward the aging process is not known. Here, we place the TE gypsy under the control of the UAS GAL4 system to model TE activation during aging. We find that increased TE activity shortens the life span of male Drosophila melanogaster. The effect is only apparent in middle-aged animals. The increase in mortality is not seen in young animals. An intact reverse transcriptase is necessary for the decrease in life span, implicating a DNA-mediated process in the effect. The decline in life span in the active gypsy flies is accompanied by the acceleration of a subset of aging phenotypes. TE activity increases sensitivity to oxidative stress and promotes a decline in circadian rhythmicity. The overexpression of the Forkhead-box O family (FOXO) stress response transcription factor can partially rescue the detrimental effects of increased TE activity on life span. Our results provide evidence that active TEs can behave as effectors in the aging process and suggest a potential novel role for dFOXO in its promotion of longevity in D. melanogaster.

Stress does not induce a general transcription of transposable elements in Drosophila

Transposable elements, also known as "jumping genes," have the ability to hop within the host genome. Nonetheless, this capacity is kept in check by the host cell defense systems to avoid unbridled TE mobilization. Different types of stressors can activate TEs in Drosophila, suggesting that TEs may play an adaptive role in the stress response, especially in generating genetic variability for adaptive evolution. TE activation by stressors may also lead to the notion, usually found in the literature, that any form of stress could activate all or the majority of TEs. In this review, we define what stress is. We then present and discuss RNA sequencing results from several studies demonstrating that stress does not trigger TE transcription broadly in Drosophila. An explanation for the LTR order of TEs being the most overexpressed is also proposed.

Epistatic selection on a selfish Segregation Distorter supergene - drive, recombination, and genetic load

Meiotic drive supergenes are complexes of alleles at linked loci that together subvert Mendelian segregation resulting in preferential transmission. In males, the most common mechanism of drive involves the disruption of sperm bearing one of a pair of alternative alleles. While at least two loci are important for male drive-the driver and the target-linked modifiers can enhance drive, creating selection pressure to suppress recombination. In this work, we investigate the evolution and genomic consequences of an autosomal, multilocus, male meiotic drive system, Segregation Distorter (SD) in the fruit fly, Drosophila melanogaster. In African populations, the predominant SD chromosome variant, SD-Mal, is characterized by two overlapping, paracentric inversions on chromosome arm 2R and nearly perfect (~100%) transmission. We study the SD-Mal system in detail, exploring its components, chromosomal structure, and evolutionary history. Our findings reveal a recent chromosome-scale selective sweep mediated by strong epistatic selection for haplotypes carrying Sd, the main driving allele, and one or more factors within the double inversion. While most SD-Mal chromosomes are homozygous lethal, SD-Mal haplotypes can recombine with other, complementing haplotypes via crossing over, and with wildtype chromosomes via gene conversion. SD-Mal chromosomes have nevertheless accumulated lethal mutations, excess non-synonymous mutations, and excess transposable element insertions. Therefore, SD-Mal haplotypes evolve as a small, semi-isolated subpopulation with a history of strong selection. These results may explain the evolutionary turnover of SD haplotypes in different populations around the world and have implications for supergene evolution broadly.

Characterization of transposable elements within the Bemisia tabaci species complex

Background

Whiteflies are agricultural pests that cause negative impacts globally to crop yields resulting at times in severe economic losses and food insecurity. The Bemisia tabaci whitefly species complex is the most damaging in terms of its broad crop host range and its ability to serve as vector for over 400 plant viruses. Genomes of whiteflies belonging to this species complex have provided valuable genomic data; however, transposable elements (TEs) within these genomes remain unexplored. This study provides the first accurate characterization of TE content within the B. tabaci species complex.

Results

This study identified that an average of 40.61% of the genomes of three whitefly species (MEAM1, MEDQ, and SSA-ECA) consists of TEs. The majority of the TEs identified were DNA transposons (22.85% average) while SINEs (0.14% average) were the least represented. This study also compared the TE content of the three whitefly genomes with three other hemipteran genomes and found significantly more DNA transposons and less LINEs in the whitefly genomes. A total of 63 TE superfamilies were identified to be present across the three whitefly species (39 DNA transposons, six LTR, 16 LINE, and two SINE). The sequences of the identified TEs were clustered which generated 5766 TE clusters. A total of 2707 clusters were identified as uniquely found within the whitefly genomes while none of the generated clusters were from both whitefly and non-whitefly TE sequences.

This study is the first to characterize TEs found within different B. tabaci species and has created a standardized annotation workflow that could be used to analyze future whitefly genomes.

Conclusion

This study is the first to characterize the landscape of TEs within the B. tabaci whitefly species complex. The characterization of these elements within the three whitefly genomes shows that TEs occupy significant portions of B. tabaci genomes, with DNA transposons representing the vast majority. This study also identified TE superfamilies and clusters of TE sequences of potential interest, providing essential information, and a framework for future TE studies within this species complex.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13100-022-00270-6.

Population-scale long-read sequencing uncovers transposable elements associated with gene expression variation and adaptive signatures in Drosophila

High quality reference genomes are crucial to understanding genome function, structure and evolution. The availability of reference genomes has allowed us to start inferring the role of genetic variation in biology, disease, and biodiversity conservation. However, analyses across organisms demonstrate that a single reference genome is not enough to capture the global genetic diversity present in populations. In this work, we generate 32 high-quality reference genomes for the well-known model species D. melanogaster and focus on the identification and analysis of transposable element variation as they are the most common type of structural variant. We show that integrating the genetic variation across natural populations from five climatic regions increases the number of detected insertions by 58%. Moreover, 26% to 57% of the insertions identified using long-reads were missed by short-reads methods. We also identify hundreds of transposable elements associated with gene expression variation and new TE variants likely to contribute to adaptive evolution in this species. Our results highlight the importance of incorporating the genetic variation present in natural populations to genomic studies, which is essential if we are to understand how genomes function and evolve.

Satellitome Analysis and Transposable Elements Comparison in Geographically Distant Populations of Spodoptera frugiperda

Spodoptera frugiperda (fall armyworm) is a member of the superfamily Noctuoidea that accounts for more than a third of all Lepidoptera and includes a considerable number of agricultural and forest pest species. Spodoptera frugiperda is a polyphagous species that is a significant agricultural pest worldwide, emphasizing its economic importance. Spodoptera frugiperda's genome size, assembly, phylogenetic classification, and transcriptome analysis have all been previously described. However, the different studies reported different compositions of repeated DNA sequences that occupied the whole assembled genome, and the Spodoptera frugiperda genome also lacks the comprehensive study of dynamic satellite DNA. We conducted a comparative analysis of repetitive DNA across geographically distant populations of Spodoptera frugiperda, particularly satellite DNA, using publicly accessible raw genome data from eight different geographical regions. Our results showed that most transposable elements (TEs) were commonly shared across all geographically distant samples, except for the Maverick and PIF/Harbinger elements, which have divergent repeat copies. The TEs age analysis revealed that most TEs families consist of young copies 1-15 million years old; however, PIF/Harbinger has some older/degenerated copies of 30-35 million years old. A total of seven satellite DNA families were discovered, accounting for approximately 0.65% of the entire genome of the Spodoptera frugiperda fall armyworm. The repeat profiling analysis of satellite DNA families revealed differential read depth coverage or copy numbers. The satellite DNA families range in size from the lowest 108 bp SfrSat06-108 families to the largest (1824 bp) SfrSat07-1824 family. We did not observe a statistically significant correlation between monomer length and K2P divergence, copy number, or abundance of each satellite family. Our findings suggest that the satellite DNA families identified in Spodoptera frugiperda account for a considerable proportion of the genome's repetitive fraction. The satellite DNA families' repeat profiling revealed a point mutation along the reference sequences. Limited TEs differentiation exists among geographically distant populations of Spodoptera frugiperda.

Meiotic, genomic and evolutionary properties of crossover distribution in Drosophila yakuba

The number and location of crossovers across genomes are highly regulated during meiosis, yet the key components controlling them are fast evolving, hindering our understanding of the mechanistic causes and evolutionary consequences of changes in crossover rates. Drosophila melanogaster has been a model species to study meiosis for more than a century, with an available high-resolution crossover map that is, nonetheless, missing for closely related species, thus preventing evolutionary context. Here, we applied a novel and highly efficient approach to generate whole-genome high-resolution crossover maps in D. yakuba to tackle multiple questions that benefit from being addressed collectively within an appropriate phylogenetic framework, in our case the D. melanogaster species subgroup. The genotyping of more than 1,600 individual meiotic events allowed us to identify several key distinct properties relative to D. melanogaster. We show that D. yakuba, in addition to higher crossover rates than D. melanogaster, has a stronger centromere effect and crossover assurance than any Drosophila species analyzed to date. We also report the presence of an active crossover-associated meiotic drive mechanism for the X chromosome that results in the preferential inclusion in oocytes of chromatids with crossovers. Our evolutionary and genomic analyses suggest that the genome-wide landscape of crossover rates in D. yakuba has been fairly stable and captures a significant signal of the ancestral crossover landscape for the whole D. melanogaster subgroup, even informative for the D. melanogaster lineage. Contemporary crossover rates in D. melanogaster, on the other hand, do not recapitulate ancestral crossovers landscapes. As a result, the temporal stability of crossover landscapes observed in D. yakuba makes this species an ideal system for applying population genetic models of selection and linkage, given that these models assume temporal constancy in linkage effects. Our studies emphasize the importance of generating multiple high-resolution crossover rate maps within a coherent phylogenetic context to broaden our understanding of crossover control during meiosis and to improve studies on the evolutionary consequences of variable crossover rates across genomes and time.

Transposable element landscapes in aging Drosophila

Genetic mechanisms that repress transposable elements (TEs) in young animals decline during aging, as reflected by increased TE expression in aged animals. Does increased TE expression during aging lead to more genomic TE copies in older animals? To address this question, we quantified TE Landscapes (TLs) via whole genome sequencing of young and aged Drosophila strains of wild-type and mutant backgrounds. We quantified TLs in whole flies and dissected brains and validated the feasibility of our approach in detecting new TE insertions in aging Drosophila genomes when small RNA and RNA interference (RNAi) pathways are compromised. We also describe improved sequencing methods to quantify extra-chromosomal DNA circles (eccDNAs) in Drosophila as an additional source of TE copies that accumulate during aging. Lastly, to combat the natural progression of aging-associated TE expression, we show that knocking down PAF1, a conserved transcription elongation factor that antagonizes RNAi pathways, may bolster suppression of TEs during aging and extend lifespan. Our study suggests that in addition to a possible influence by different genetic backgrounds, small RNA and RNAi mechanisms may mitigate genomic TL expansion despite the increase in TE transcripts during aging.

Transposable elements, also called transposons, are genetic parasites found in all animal genomes. Normally, transposons are compacted away in silent chromatin in young animals. But, as animals age and transposon-silencing defense mechanisms break down, transposon RNAs accumulate to significant levels in old animals like fruit flies. An open question is whether the increased levels of transposon RNAs in older animals also correspond to increased genomic copies of transposons. This study approached this question by sequencing the whole genomes of young and old wild-type and mutant flies lacking a functional RNA interference (RNAi) pathway, which naturally silences transposon RNAs. Although the wild-type flies with intact RNAi activity had little new accumulation of transposon copies, the sequencing approach was able to detect several transposon accumulation occurrences in some RNAi mutants. In addition, we found that some fly transposon families can also accumulate as extra-chromosomal circular DNA copies. Lastly, we showed that genetically augmenting the expression of RNAi factors can counteract the rising transposon RNA levels in aging and promote longevity. This study improves our understanding of the animal host genome relationship with transposons during natural aging processes.

Constitutive Heterochromatin in Eukaryotic Genomes: A Mine of Transposable Elements

Transposable elements (TEs) are abundant components of constitutive heterochromatin of the most diverse evolutionarily distant organisms. TEs enrichment in constitutive heterochromatin was originally described in the model organism Drosophila melanogaster, but it is now considered as a general feature of this peculiar portion of the genomes. The phenomenon of TE enrichment in constitutive heterochromatin has been proposed to be the consequence of a progressive accumulation of transposable elements caused by both reduced recombination and lack of functional genes in constitutive heterochromatin. However, this view does not take into account classical genetics studies and most recent evidence derived by genomic analyses of heterochromatin in Drosophila and other species. In particular, the lack of functional genes does not seem to be any more a general feature of heterochromatin. Sequencing and annotation of Drosophila melanogaster constitutive heterochromatin have shown that this peculiar genomic compartment contains hundreds of transcriptionally active genes, generally larger in size than that of euchromatic ones. Together, these genes occupy a significant fraction of the genomic territory of heterochromatin. Moreover, transposable elements have been suggested to drive the formation of heterochromatin by recruiting HP1 and repressive chromatin marks. In addition, there are several pieces of evidence that transposable elements accumulation in the heterochromatin might be important for centromere and telomere structure. Thus, there may be more complexity to the relationship between transposable elements and constitutive heterochromatin, in that different forces could drive the dynamic of this phenomenon. Among those forces, preferential transposition may be an important factor. In this article, we present an overview of experimental findings showing cases of transposon enrichment into the heterochromatin and their positive evolutionary interactions with an impact to host genomes.

piRNAs of Caenorhabditis elegans broadly silence nonself sequences through functionally random targeting

Small noncoding RNAs such as piRNAs are guides for Argonaute proteins, enabling sequence-specific, post-transcriptional regulation of gene expression. The piRNAs of Caenorhabditis elegans have been observed to bind targets with high mismatch tolerance and appear to lack specific transposon targets, unlike piRNAs in Drosophila melanogaster and other organisms. These observations support a model in which C. elegans piRNAs provide a broad, indiscriminate net of silencing, competing with siRNAs associated with the CSR-1 Argonaute that specifically protect self-genes from silencing. However, the breadth of piRNA targeting has not been subject to in-depth quantitative analysis, nor has it been explained how piRNAs are distributed across sequence space to achieve complete coverage. Through a bioinformatic analysis of piRNA sequences, incorporating an original data-based metric of piRNA-target distance, we demonstrate that C. elegans piRNAs are functionally random, in that their coverage of sequence space is comparable to that of random sequences. By possessing a sufficient number of distinct, essentially random piRNAs, C. elegans is able to target arbitrary nonself sequences with high probability. We extend this approach to a selection of other nematodes, finding results which elucidate the mechanism by which nonself mRNAs are silenced, and have implications for piRNA evolution and biogenesis.

The Ribosomal Protein RpL22 Interacts In Vitro with 5′-UTR Sequences Found in Some Drosophila melanogaster Transposons

Mobility of eukaryotic transposable elements (TEs) are finely regulated to avoid an excessive mutational load caused by their movement. The transposition of retrotransposons is usually regulated through the interaction of host- and TE-encoded proteins, with non-coding regions (LTR and 5′-UTR) of the transposon. Examples of new potent cis-acting sequences, identified and characterized in the non-coding regions of retrotransposons, include the insulator of gypsy and Idefix, and the enhancer of ZAM of Drosophila melanogaster. Recently we have shown that in the 5′-UTR of the LTR-retrotransposon ZAM there is a sequence structured in tandem-repeat capable of operating as an insulator both in Drosophila (S2R⁺) and human cells (HEK293). Here, we test the hypothesis that tandem repeated 5′-UTR of a different LTR-retrotransposon could accommodate similar regulatory elements. The comparison of the 5′-UTR of some LTR-transposons allowed us to identify a shared motif of 13 bp, called Transposable Element Redundant Motif (TERM). Surprisingly, we demonstrated, by Yeast One-Hybrid assay, that TERM interacts with the D. melanogaster ribosomal protein RpL22. The Drosophila RpL22 has additional Ala-, Lys- and Pro-rich sequences at the amino terminus, which resembles the carboxy-terminal portion of histone H1 and histone H5. For this reason, it has been hypothesized that RpL22 might have two functions, namely the role in organizing the ribosome, and a potential regulatory role involving DNA-binding similar to histone H1, which represses transcription in Drosophila. In this paper, we show, by two independent sets of experiments, that DmRpL22 is able to directly and specifically bind DNA of Drosophila melanogaster.

Synergistic epistasis of the deleterious effects of transposable elements

The replicative nature and generally deleterious effects of transposable elements (TEs) raise an outstanding question about how TE copy number is stably contained in host populations. Classic theoretical analyses predict that, when the decline in fitness due to each additional TE insertion is greater than linear, or when there is synergistic epistasis, selection against TEs can result in a stable equilibrium of TE copy number. While several mechanisms are predicted to yield synergistic deleterious effects of TEs, we lack empirical investigations of the presence of such epistatic interactions. Purifying selection with synergistic epistasis generates repulsion linkage between deleterious alleles. We investigated this population genetic signal in the likely ancestral Drosophila melanogaster population and found evidence supporting the presence of synergistic epistasis among TE insertions, especially TEs expected to exert large fitness impacts. Even though synergistic epistasis of TEs has been predicted to arise through ectopic recombination and TE-mediated epigenetic silencing mechanisms, we only found mixed support for the associated predictions. We observed signals of synergistic epistasis for a large number of TE families, which is consistent with the expectation that such epistatic interaction mainly happens among copies of the same family. Curiously, significant repulsion linkage was also found among TE insertions from different families, suggesting the possibility that synergism of TEs' deleterious fitness effects could arise above the family level and through mechanisms similar to those of simple mutations. Our findings set the stage for investigating the prevalence and importance of epistatic interactions in the evolutionary dynamics of TEs.

A genomic survey of LINE elements in Pipidae aquatic frogs shed light on Rex-elements evolution in these genomes

The transposable elements (TE) represent a large portion of anuran genomes that act as components of genetic diversification. The LINE order of retrotransposons is among the most representative and diverse TEs and is poorly investigated in anurans. Here we explored the LINE diversity with an emphasis on the elements generically called Rex in Pipidae species, more specifically, in the genomes ofXenopus tropicalis, used as a model genome in the study of anurans,the allotetraploid sister species Xenopus laevis and theAmerican species Pipa carvalhoi. We were able to identify a great diversity of LINEs from five clades, Rex1, L2, CR1, L1 and Tx1, in these three species, and the RTE clade was lost in X. tropicalis. It is clear that elements classified as Rex are distributed in distinct clades. The evolutionary pattern of Rex1 elements denote a complex evolution with independent losses of families and some horizontal transfer events between fishes and amphibians which were supported not only by the phylogenetic inconsistencies but also by the very low Ks values found for the TE sequences. The data obtained here update the knowledge of the LINEs diversity in X. laevis and represent the first study of TEs in P. carvalhoi.

Patterns of piRNA Regulation in Drosophila Revealed through Transposable Element Clade Inference

Transposable elements (TEs) are self-replicating "genetic parasites" ubiquitous to eukaryotic genomes. In addition to conflict between TEs and their host genomes, TEs of the same family are in competition with each other. They compete for the same genomic niches while experiencing the same regime of copy-number selection. This suggests that competition among TEs may favor the emergence of new variants that can outcompete their ancestral forms. To investigate the sequence evolution of TEs, we developed a method to infer clades: collections of TEs that share SNP variants and represent distinct TE family lineages. We applied this method to a panel of 85 Drosophila melanogaster genomes and found that the genetic variation of several TE families shows significant population structure that arises from the population-specific expansions of single clades. We used population genetic theory to classify these clades into younger versus older clades and found that younger clades are associated with a greater abundance of sense and antisense piRNAs per copy than older ones. Further, we find that the abundance of younger, but not older clades, is positively correlated with antisense piRNA production, suggesting a general pattern where hosts preferentially produce antisense piRNAs from recently active TE variants. Together these findings suggest a pattern whereby new TE variants arise by mutation and then increase in copy number, followed by the host producing antisense piRNAs that may be used to silence these emerging variants.

Chromatin remodelling and retrotransposons activities during regeneration in Drosophila

Regeneration is a response mechanism aiming to reconstruct lost or damaged structures. To achieve this, the cells repopulating the lost tissue often have to change their original identity, a process that involves chromatin remodelling.We have analysed the issue of chromatin remodelling during regeneration in the wing disc of Drosophila . In this disc the ablation of the central region (the pouch) induces the regenerative response of the cells from the lateral region (the hinge), which reconstitute the wing pouch. We have examined euchromatin and heterochromatin histone marks during the process and find that heterochromatin marks disappear but are recovered when regeneration is complete. Euchromatin marks are not modified. We also describe the transcription of two retrotransposons, Roo and F-element in the regenerating cells. We have established a temporal correlation between the alterations of heterochromatin marks and the levels of transcription of two retrotransposons, Roo and F-element, both during embryonic development and in the regeneration process.

Large Drosophila germline piRNA clusters are evolutionarily labile and dispensable for transposon regulation

PIWI proteins and their guiding Piwi-interacting small RNAs (piRNAs) are crucial for fertility and transposon defense in the animal germline. In most species, the majority of piRNAs are produced from distinct large genomic loci, called piRNA clusters. It is assumed that germline-expressed piRNA clusters, particularly in Drosophila, act as principal regulators to control transposons dispersed across the genome. Here, using synteny analysis, we show that large clusters are evolutionarily labile, arise at loci characterized by recurrent chromosomal rearrangements, and are mostly species-specific across the Drosophila genus. By engineering chromosomal deletions in D. melanogaster, we demonstrate that the three largest germline clusters, which account for the accumulation of >40% of all transposon-targeting piRNAs in ovaries, are neither required for fertility nor for transposon regulation in trans. We provide further evidence that dispersed elements, rather than the regulatory action of large Drosophila germline clusters in trans, may be central for transposon defense.

The Drosophila ribosome protein S5 paralog RpS5b promotes germ cell and follicle cell differentiation during oogenesis

Emerging evidence suggests that ribosome heterogeneity may have important functional consequences in the translation of specific mRNAs within different cell types and under various conditions. Ribosome heterogeneity comes in many forms, including post-translational modification of ribosome proteins (RPs), absence of specific RPs and inclusion of different RP paralogs. The Drosophila genome encodes two RpS5 paralogs: RpS5a and RpS5b. While RpS5a is ubiquitously expressed, RpS5b exhibits enriched expression in the reproductive system. Deletion of RpS5b results in female sterility marked by developmental arrest of egg chambers at stages 7-8, disruption of vitellogenesis and posterior follicle cell (PFC) hyperplasia. While transgenic rescue experiments suggest functional redundancy between RpS5a and RpS5b, molecular, biochemical and ribo-seq experiments indicate that RpS5b mutants display increased rRNA transcription and RP production, accompanied by increased protein synthesis. Loss of RpS5b results in microtubule-based defects and in mislocalization of Delta and Mindbomb1, leading to failure of Notch pathway activation in PFCs. Together, our results indicate that germ cell-specific expression of RpS5b promotes proper egg chamber development by ensuring the homeostasis of functional ribosomes.

Diversity and Evolution of pogo and Tc1/mariner Transposons in the Apoidea Genomes

Bees (Apoidea), the largest and most crucial radiation of pollinators, play a vital role in the ecosystem balance. Transposons are widely distributed in nature and are important drivers of species diversity. However, transposons are rarely reported in important pollinators such as bees. Here, we surveyed 37 bee genomesin Apoidea, annotated the pogo and Tc1/mariner transposons in the genome of each species, and performed a phylogenetic analysis and determined their overall distribution. The pogo and Tc1/mariner families showed high diversity and low abundance in the 37 species, and their proportion was significantly higher in solitary bees than in social bees. DD34D/mariner was found to be distributed in almost all species and was found in Apis mellifera, Apis mellifera carnica, Apis mellifera caucasia, and Apis mellifera mellifera, and Euglossa dilemma may still be active. Using horizontal transfer analysis, we found that DD29-30D/Tigger may have experienced horizontal transfer (HT) events. The current study displayed the evolution profiles (including diversity, activity, and abundance) of the pogo and Tc1/mariner transposons across 37 species of Apoidea. Our data revealed their contributions to the genomic variations across these species and facilitated in understanding of the genome evolution of this lineage.

Phenotypic and Transcriptomic Responses to Stress Differ According to Population Geography in an Invasive Species

Adaptation to rapid environmental changes must occur within a short-time scale. In this context, studies of invasive species may provide insights into the underlying mechanisms of rapid adaptation as these species have repeatedly encountered and adapted to novel environmental conditions. We investigated how invasive and noninvasive genotypes of Drosophila suzukii deal with oxidative stress at the phenotypic and molecular levels. We also studied the impact of transposable element (TE) insertions on the gene expression in response to stress. Our results show that flies from invasive areas (France and the United States) live longer in natural conditions than the ones from native Japanese areas. As expected, lifespan for all genotypes was significantly reduced following exposure to paraquat, but this reduction varied among genotypes (genotype-by-environment interaction) with invasive genotypes appearing more affected by exposure than noninvasive ones. A transcriptomic analysis of genotypes upon paraquat treatment detected many genes differentially expressed (DE). Although a small core set of genes were DE in all genotypes following paraquat exposure, much of the response of each genotype was unique. Moreover, we showed that TEs were not activated after oxidative stress and DE genes were significantly depleted of TEs. In conclusion, it is likely that transcriptomic changes are involved in the rapid adaptation to local environments. We provide new evidence that in the decade since the invasion from Asia, the sampled genotypes in Europe and the United States of D. suzukii diverged from the ones from the native area regarding their phenotypic and genomic response to oxidative stress.

piRNAs and endo-siRNAs: Small molecules with large roles in the nervous system

Since their discovery, small non-coding RNAs have emerged as powerhouses in the regulation of numerous cellular processes. In addition to guarding the integrity of the reproductive system, small non-coding RNAs play critical roles in the maintenance of the soma. Accumulating evidence indicates that small non-coding RNAs perform vital functions in the animal nervous system such as restricting the activity of deleterious transposable elements, regulating nerve regeneration, and mediating learning and memory. In this review, we provide an overview of the current understanding of the contribution of two major classes of small non-coding RNAs, piRNAs and endo-siRNAs, to the nervous system development and function, and present highlights on how the dysregulation of small non-coding RNA pathways can assist in understanding the neuropathology of human neurological disorders.

Transposons-Based Clonal Diversity in Trematode Involves Parts of CR1 (LINE) in Eu- and Heterochromatin

Trematode parthenitae have long been believed to form clonal populations, but clonal diversity has been discovered in this asexual stage of the lifecycle. Clonal polymorphism in the model species Himasthla elongata has been previously described, but the source of this phenomenon remains unknown. In this work, we traced cercarial clonal diversity using a simplified amplified fragment length polymorphism (SAFLP) method and characterised the nature of fragments in diverse electrophoretic bands. The repetitive elements were identified in both the primary sequence of the H. elongata genome and in the transcriptome data. Long-interspersed nuclear elements (LINEs) and long terminal repeat retrotransposons (LTRs) were found to represent an overwhelming majority of the genome and the transposon transcripts. Most sequenced fragments from SAFLP pattern contained the reverse transcriptase (RT, ORF2) domains of LINEs, and only a few sequences belonged to ORFs of LTRs and ORF1 of LINEs. A fragment corresponding to a CR1-like (LINE) spacer region was discovered and named CR1-renegade (CR1-rng). In addition to RT-containing CR1 transcripts, we found short CR1-rng transcripts in the redia transcriptome and short contigs in the mobilome. Probes against CR1-RT and CR1-rng presented strikingly different pictures in FISH mapping, despite both being fragments of CR1. In silico data and Southern blotting indicated that CR1-rng is not tandemly organised. CR1 involvement in clonal diversity is discussed.

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.

Genome-wide characterization of LTR retrotransposons in the non-model deep-sea annelid Lamellibrachia luymesi

Background

Long Terminal Repeat retrotransposons (LTR retrotransposons) are mobile genetic elements composed of a few genes between terminal repeats and, in some cases, can comprise over half of a genome’s content. Available data on LTR retrotransposons have facilitated comparative studies and provided insight on genome evolution. However, data are biased to model systems and marine organisms, including annelids, have been underrepresented in transposable elements studies. Here, we focus on genome of Lamellibrachia luymesi, a vestimentiferan tubeworm from deep-sea hydrocarbon seeps, to gain knowledge of LTR retrotransposons in a deep-sea annelid.

Results

We characterized LTR retrotransposons present in the genome of L. luymesi using bioinformatic approaches and found that intact LTR retrotransposons makes up about 0.1% of L. luymesi genome. Previous characterization of the genome has shown that this tubeworm hosts several known LTR-retrotransposons. Here we describe and classify LTR retrotransposons in L. luymesi as within the Gypsy, Copia and Bel-pao superfamilies. Although, many elements fell within already recognized families (e.g., Mag, CSRN1), others formed clades distinct from previously recognized families within these superfamilies. However, approximately 19% (41) of recovered elements could not be classified. Gypsy elements were the most abundant while only 2 Copia and 2 Bel-pao elements were present. In addition, analysis of insertion times indicated that several LTR-retrotransposons were recently transposed into the genome of L. luymesi, these elements had identical LTR’s raising possibility of recent or ongoing retrotransposon activity.

Conclusions

Our analysis contributes to knowledge on diversity of LTR-retrotransposons in marine settings and also serves as an important step to assist our understanding of the potential role of retroelements in marine organisms. We find that many LTR retrotransposons, which have been inserted in the last few million years, are similar to those found in terrestrial model species. However, several new groups of LTR retrotransposons were discovered suggesting that the representation of LTR retrotransposons may be different in marine settings. Further study would improve understanding of the diversity of retrotransposons across animal groups and environments.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07749-1.

Behavioral and Genomic Sensory Adaptations Underlying the Pest Activity of Drosophila suzukii

Studying how novel phenotypes originate and evolve is fundamental to the field of evolutionary biology as it allows us to understand how organismal diversity is generated and maintained. However, determining the basis of novel phenotypes is challenging as it involves orchestrated changes at multiple biological levels. Here, we aim to overcome this challenge by using a comparative species framework combining behavioral, gene expression, and genomic analyses to understand the evolutionary novel egg-laying substrate-choice behavior of the invasive pest species Drosophila suzukii. First, we used egg-laying behavioral assays to understand the evolution of ripe fruit oviposition preference in D. suzukii compared with closely related species D. subpulchrella and D. biarmipes as well as D. melanogaster. We show that D. subpulchrella and D. biarmipes lay eggs on both ripe and rotten fruits, suggesting that the transition to ripe fruit preference was gradual. Second, using two-choice oviposition assays, we studied how D. suzukii, D. subpulchrella, D. biarmipes, and D. melanogaster differentially process key sensory cues distinguishing ripe from rotten fruit during egg-laying. We found that D. suzukii's preference for ripe fruit is in part mediated through a species-specific preference for stiff substrates. Last, we sequenced and annotated a high-quality genome for D. subpulchrella. Using comparative genomic approaches, we identified candidate genes involved in D. suzukii's ability to seek out and target ripe fruits. Our results provide detail to the stepwise evolution of pest activity in D. suzukii, indicating important cues used by this species when finding a host, and the molecular mechanisms potentially underlying their adaptation to a new ecological niche.