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.

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.

The annotation of transposable elements (transposons) is a very dynamic field of genomics and various tools assigned to support this bioinformatics endeavor have been developed and described. Genome ARTIST v1.19 (GA_v1.19) software was conceived for mapping artificial transposons mobilized during insertional mutagenesis projects, but the new functions of GA_v2 qualify it as a tool for the mapping and annotation of natural transposons (NTs) in long reads, contigs and assembled genomes. The tabular export of mapping and annotation data for high-throughput data analysis, the generation of a list of flanking sequences around the coordinates of insertion or around the target site duplications and the computing of a consensus sequence for the flanking sequences are all key assets of GA_v2. Additionally, we developed a set of scripts that enable the user to annotate NTs, to harness annotations offered by FlyBase for Drosophila melanogaster genome, to convert sequence files from .fasta to .raw, and to extract junction query sequences essential for NTs mapping. Herein, we present the applicability of GA_v2 for a preliminary annotation of P-element and hobo class II NTs and copia retrotransposon in the genome of D. melanogaster strain Horezu_LaPeri (Horezu), Romania, which was sequenced with Nanopore technology in our laboratory. We used contigs assembled with Flye tool and a Q10 quality filter of the reads. Our results suggest that GA_v2 is a reliable autonomous tool able to perform mapping and annotation of NTs in genomes sequenced by long sequencing technology. GA_v2 is open-source software compatible with Linux, Mac OS and Windows and is available at GitHub repository and dedicated website.

Drosophila melanogaster (the fruit fly) is a valuable experimental platform for modeling host–pathogen interactions. It is also commonly used to define innate immunity pathways and to understand the mechanisms of both host tolerance to commensal microbiota and response to pathogenic agents. Herein, we investigate how the host response to bacterial infection is mirrored in the expression of genes of Imd and Toll pathways when D. melanogaster strains with different γCOP genetic backgrounds are infected with Pseudomonas aeruginosa ATCC 27853. Using microarray technology, we have interrogated the whole-body transcriptome of infected versus uninfected fruit fly males with three specific genotypes, namely wild-type Oregon, γCOP^S057302/TM6B and γCOP^14a/γCOP^14a. While the expression of genes pertaining to Imd and Toll is not significantly modulated by P. aeruginosa infection in Oregon males, many of the components of these cascades are up- or downregulated in both infected and uninfected γCOP^S057302/TM6B and γCOP^14a/γCOP^14a males. Thus, our results suggest that a γCOP genetic background modulates the gene expression profiles of Imd and Toll cascades involved in the innate immune response of D. melanogaster, inducing the occurrence of immunological dysfunctions in γCOP mutants.

Transposons are mobile DNA sequences, known for their ability to insert into other locations in the genome. Genome sequencing allowed the identification of the high content of transposons in various model organisms. Here, we present a bioinformatics pipeline developed to estimate the evolutionary relationship between Drosophila melanogaster and other Drosophilidae based on the comparative analysis of the presence and distribution of class II transposons. Our study reveals that the presence and distribution of transposons hobo, HB, Tc1, Tc1-2, hopper and Bari1 points to close evolutionary relationship among D. melanogaster, D. simulans, D. sechellia and D. yakuba , which is in accordance with other data available in literature.