The Juan non-LTR retrotransposon in mosquitoes: genomic impact, vertical transmission and indications of recent and widespread activity

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.

Mosquitoes LTR retrotransposons: a deeper view into the genomic sequence of Culex quinquefasciatus

A set of 67 novel LTR-retrotransposon has been identified by in silico analyses of the Culex quinquefasciatus genome using the LTR_STRUC program. The phylogenetic analysis shows that 29 novel and putatively functional LTR-retrotransposons detected belong to the Ty3/gypsy group. Our results demonstrate that, by considering only families containing potentially autonomous LTR-retrotransposons, they account for about 1% of the genome of C. quinquefasciatus. In previous studies it has been estimated that 29% of the genome of C. quinquefasciatus is occupied by mobile genetic elements.The potential role of retrotransposon insertions strictly associated with host genes is described and discussed along with the possible origin of a retrotransposon with peculiar Primer Binding Site region. Finally, we report the presence of a group of 38 retrotransposons, carrying tandem repeated sequences but lacking coding potential, and apparently lacking "master copy" elements from which they could have originated. The features of the repetitive sequences found in these non-autonomous LTR retrotransposons are described, and their possible role discussed.These results integrate the existing data on the genomics of an important virus-borne disease vector.

Identifying insecticide resistance genes in mosquito by combining AFLP genome scans and 454 pyrosequencing

AFLP-based genome scans are widely used to study the genetics of adaptation and to identify genomic regions potentially under selection. However, this approach usually fails to detect the actual genes or mutations targeted by selection owing to the difficulties of obtaining DNA sequences from AFLP fragments. Here, we combine classical AFLP outlier detection with 454 sequencing of AFLP fragments to obtain sequences from outlier loci. We applied this approach to the study of resistance to Bacillus thuringiensis israelensis (Bti) toxins in the dengue vector Aedes aegypti. A genome scan of Bti-resistant and Bti-susceptible A. aegypti laboratory strains was performed based on 432 AFLP markers. Fourteen outliers were detected using two different population genetic algorithms. Out of these, 11 were successfully sequenced. Three contained transposable elements (TEs) sequences, and the 10 outliers that could be mapped at a unique location in the reference genome were located on different supercontigs. One outlier was in the vicinity of a gene coding for an aminopeptidase potentially involved in Bti toxin-binding. Patterns of sequence variability of this gene showed significant deviation from neutrality in the resistant strain but not in the susceptible strain, even after taking into account the known demographic history of the selected strain. This gene is a promising candidate for future functional analysis.

The mosquito Aedes aegypti has a large genome size and high transposable element load but contains a low proportion of transposon-specific piRNAs

BACKGROUND

The piRNA pathway has been shown in model organisms to be involved in silencing of transposons thereby providing genome stability. In D. melanogaster the majority of piRNAs map to these sequences. The medically important mosquito species Aedes aegypti has a large genome size, a high transposon load which includes Miniature Inverted repeat Transposable Elements (MITES) and an expansion of the piRNA biogenesis genes. Studies of transgenic lines of Ae. aegypti have indicated that introduced transposons are poorly remobilized and we sought to explore the basis of this. We wished to analyze the piRNA profile of Ae. aegypti and thereby determine if it is responsible for transposon silencing in this mosquito.

RESULTS

Estimated piRNA sequence diversity was comparable between Ae. aegypti and D. melanogaster, but surprisingly only 19% of mosquito piRNAs mapped to transposons compared to 51% for D. melanogaster. Ae. aegypti piRNA clusters made up a larger percentage of the total genome than those of D. melanogaster but did not contain significantly higher percentages of transposon derived sequences than other regions of the genome. Ae. aegypti contains a number of protein coding genes that may be sources of piRNA biogenesis with two, traffic jam and maelstrom, implicated in this process in model organisms. Several genes of viral origin were also targeted by piRNAs. Examination of six mosquito libraries that had previously been transformed with transposon derived sequence revealed that new piRNA sequences had been generated to the transformed sequences, suggesting that they may have stimulated a transposon inactivation mechanism.

CONCLUSIONS

Ae. aegypti has a large piRNA complement that maps to transposons but primarily gene sequences, including many viral-derived sequences. This, together the more uniform distribution of piRNA clusters throughout its genome, suggest that some aspects of the piRNA system differ between Ae. aegypti and D. melanogaster.

In silico evidence for the species-specific conservation of mosquito retroposons: implications as a molecular biomarker

BACKGROUND

Mosquitoes are the transmissive vectors for several infectious pathogens that affect man. However, the control of mosquitoes through insecticide and pesticide spraying has proved difficult in the past. We hypothesized that, by virtue of their reported vertical inheritance among mosquitoes, group II introns - a class of small coding ribonucleic acids (scRNAs) - may form a potential species-specific biomarker. Structurally, introns are a six-moiety complex. Depending on the function of the protein encoded within the IV moiety, the highly mobile class of group II introns or retroposons is sub-divided into two: Restriction Endonuclease (REase)-like and Apurinic aPyramydinic Endonuclease (APE)-like. REase-like retroposons are thought to be the ancestors of APE retroposons. Our aim in this study was to find evidence for the highly species-specific conservation of the APE subclass of mosquito retroposons.

METHODS AND RESULTS

In silico targeted sequence alignments were conducted across a 1,779-organism genome database (1,518 bacterial, 59 archeal, 201 eukaryotic, and the human), using three mosquito retroposon sequence tags (RST) as BLASTN queries [AJ970181 and AJ90201 of Culex pipien origin and AJ970301 of Anoplese sinensis origin]. At a calibration of E = 10, A & D = 100, default filtration and a homology cut-off of >95% identity, no hits were found on any of the 1,518 bacterial genomes. Eleven (100%) and 15 (100%) hits obtained on the 201-eukaryote genome database were homologs (>95% score) of C.pipien quinquefasciatus JHB retroposons, but none of An. sinensis. Twenty and 221 low score (30-43% identity) spurious hits were found at flanking ends of genes and contigs in the human genome with the C.pipien and An. sinensis RSTs respectively. Functional and positional inference revealed these to be possible relatives of human genomic spliceosomes. We advance two models for the application of mosquito RST: as precursors for developing molecular biomarkers for mosquitoes, and as RST-specific monoclonal antibody (MAb)-DDT immunoconjugates to enhance targeted toxicity.

CONCLUSION

We offer evidence to support the species-specific conservation of mosquito retroposons among lower taxa. Our findings suggest that retroposons may therefore constitute a unique biomarker for mosquito species that may be exploited in molecular entomology. Mosquito RST-specific MAbs may possibly permit synthesis of DDT immunoconjugates that could be used to achieve species-tailored toxicity.

The tempo and mode of evolution of transposable elements as revealed by molecular phylogenies reconstructed from mosquito genomes

Although many mathematical models exist predicting the dynamics of transposable elements (TEs), there is a lack of available empirical data to validate these models and inherent assumptions. Genomes can provide a snapshot of several TE families in a single organism, and these could have their demographics inferred by coalescent analysis, allowing for the testing of theories on TE amplification dynamics. Using the available genomes of the mosquitoes Aedes aegypti and Anopheles gambiae, we indicate that such an approach is feasible. Our analysis follows four steps: (1) mining the two mosquito genomes currently available in search of TE families; (2) fitting, to selected families found in (1), a phylogeny tree under the general time-reversible (GTR) nucleotide substitution model with an uncorrelated lognormal (UCLN) relaxed clock and a nonparametric demographic model; (3) fitting a nonparametric coalescent model to the tree generated in (2); and (4) fitting parametric models motivated by ecological theories to the curve generated in (3).

Comparative genomics of small RNA regulatory pathway components in vector mosquitoes

Background

Small RNA regulatory pathways (SRRPs) control key aspects of development and anti-viral defense in metazoans. Members of the Argonaute family of catalytic enzymes degrade target RNAs in each of these pathways. SRRPs include the microRNA, small interfering RNA (siRNA) and PIWI-type gene silencing pathways. Mosquitoes generate viral siRNAs when infected with RNA arboviruses. However, in some mosquitoes, arboviruses survive antiviral RNA interference (RNAi) and are transmitted via mosquito bite to a subsequent host. Increased knowledge of these pathways and functional components should increase understanding of the limitations of anti-viral defense in vector mosquitoes. To do this, we compared the genomic structure of SRRP components across three mosquito species and three major small RNA pathways.

Results

The Ae. aegypti, An. gambiae and Cx. pipiens genomes encode putative orthologs for all major components of the miRNA, siRNA, and piRNA pathways. Ae. aegypti and Cx. pipiens have undergone expansion of Argonaute and PIWI subfamily genes. Phylogenetic analyses were performed for these protein families. In addition, sequence pattern recognition algorithms MEME, MDScan and Weeder were used to identify upstream regulatory motifs for all SRRP components. Statistical analyses confirmed enrichment of species-specific and pathway-specific cis-elements over the rest of the genome.

Conclusion

Analysis of Argonaute and PIWI subfamily genes suggests that the small regulatory RNA pathways of the major arbovirus vectors, Ae. aegypti and Cx. pipiens, are evolving faster than those of the malaria vector An. gambiae and D. melanogaster. Further, protein and genomic features suggest functional differences between subclasses of PIWI proteins and provide a basis for future analyses. Common UCR elements among SRRP components indicate that 1) key components from the miRNA, siRNA, and piRNA pathways contain NF-kappaB-related and Broad complex transcription factor binding sites, 2) purifying selection has occurred to maintain common pathway-specific elements across mosquito species and 3) species-specific differences in upstream elements suggest that there may be differences in regulatory control among mosquito species. Implications for arbovirus vector competence in mosquitoes are discussed.