Transposable elements contribute to the genomic response to insecticides in Drosophila melanogaster

Comparative transcriptional analysis between susceptible and resistant populations of Aedes (Stegomyia) aegypti (Linnaeus, 1762) after malathion exposure

Aedes aegypti is an important vector of arboviruses, including dengue, chikungunya and Zika. The application of synthetic insecticides is a frequently used strategy to control this insect. Malathion is an organophosphate insecticide that was widely used in Brazil in the 1980s and 1990s to control the adult form of A. aegypti. In situations where resistance to currently used insecticides is detected, the use of malathion may be resumed as a control measure. Many studies have confirmed resistance to malathion, however, comparative studies of differential gene expression of the entire transcriptome of resistant and susceptible insects are scarce. Therefore, understanding the molecular basis of resistance to this insecticide in this species is extremely important. In this paper, we present the first transcriptomic description of susceptible and resistant strains of A. aegypti challenged with malathion. Guided transcriptome assembly resulted in 39,904 transcripts, where 2133 differentially expressed transcripts were detected, and three were validated by RT-qPCR. Enrichment analysis for these identified transcripts resulted in 13 significant pathways (padj < 0.05), 8 associated with down-regulated and 5 with up-regulated transcripts in treated resistant insects. It was possible to divide the transcripts according to the mechanism of action into three main groups: (i) genes involved in detoxification metabolic pathways; (ii) genes of proteins located in the membrane/extracellular region; and (iii) genes related to DNA integration/function. These results are important in advancing knowledge of genes related to resistance mechanisms in this insect, enabling the development of effective technologies and strategies for managing insecticide resistance.

Transposable elements in Drosophila montana from harsh cold environments

BACKGROUND

Substantial discoveries during the past century have revealed that transposable elements (TEs) can play a crucial role in genome evolution by affecting gene expression and inducing genetic rearrangements, among other molecular and structural effects. Yet, our knowledge on the role of TEs in adaptation to extreme climates is still at its infancy. The availability of long-read sequencing has opened up the possibility to identify and study potential functional effects of TEs with higher precision. In this work, we used Drosophila montana as a model for cold-adapted organisms to study the association between TEs and adaptation to harsh climates.

RESULTS

Using the PacBio long-read sequencing technique, we de novo identified and manually curated TE sequences in five Drosophila montana genomes from eco-geographically distinct populations. We identified 489 new TE consensus sequences which represented 92% of the total TE consensus in D. montana. Overall, 11-13% of the D. montana genome is occupied by TEs, which as expected are non-randomly distributed across the genome. We identified five potentially active TE families, most of them from the retrotransposon class of TEs. Additionally, we found TEs present in the five analyzed genomes that were located nearby previously identified cold tolerant genes. Some of these TEs contain promoter elements and transcription binding sites. Finally, we detected TEs nearby fixed and polymorphic inversion breakpoints.

CONCLUSIONS

Our research revealed a significant number of newly identified TE consensus sequences in the genome of D. montana, suggesting that non-model species should be studied to get a comprehensive view of the TE repertoire in Drosophila species and beyond. Genome annotations with the new D. montana library allowed us to identify TEs located nearby cold tolerant genes, and present at high population frequencies, that contain regulatory regions and are thus good candidates to play a role in D. montana cold stress response. Finally, our annotations also allow us to identify for the first time TEs present in the breakpoints of three D. montana inversions.

Protocol to isolate nuclei from Chlamydomonas reinhardtii for ATAC sequencing

The isolation of sufficient amounts of intact nuclei is essential to obtain high-resolution maps of chromatin accessibility via assay for transposase-accessible chromatin using sequencing (ATAC-seq). Here, we present a protocol for tag-free isolation of nuclei from both cell walled and cell wall-deficient strains of the green model alga Chlamydomonas reinhardtii at a suitable quality for ATAC-seq. We describe steps for nuclei isolation, quantification, and downstream ATAC-seq. This protocol is optimized to shorten the time of isolation and quantification of nuclei.

Stress Induced Activation of LTR Retrotransposons in the Drosophila melanogaster Genome

Background: Retrotransposons with long terminal repeats (LTR retrotransposons) are widespread in all groups of eukaryotes and are often both the cause of new mutations and the source of new sequences. Apart from their high activity in generative and differentiation-stage tissues, LTR retrotransposons also become more active in response to different stressors. The precise causes of LTR retrotransposons' activation in response to stress, however, have not yet been thoroughly investigated. Methods: We used RT-PCR to investigate the transcriptional profile of LTR retrotransposons and piRNA clusters in response to oxidative and chronic heat stresses. We used Oxford Nanopore sequencing to investigate the genomic environment of new insertions of the retrotransposons. We used bioinformatics methods to find the stress-induced transcription factor binding sites in LTR retrotransposons. Results: We studied the transposition activity and transcription level of LTR retrotransposons in response to oxidative and chronic heat stress and assessed the contribution of various factors that can affect the increase in their expression under stress conditions: the state of the piRNA-interference system, the influence of the genomic environment on individual copies, and the presence of the stress-induced transcription factor binding sites in retrotransposon sequences. Conclusions: The main reason for the activation of LTR retrotransposons under stress conditions is the presence of transcription factor binding sites in their regulatory sequences, which are triggered in response to stress and are necessary for tissue regeneration processes. Stress-induced transposable element activation can function as a trigger mechanism, triggering multiple signal pathways and resulting in a polyvariant cell response.

Transposon accumulation at xenobiotic gene family loci in aphids

The evolution of resistance is a major challenge for the sustainable control of pests and pathogens. Thus, a deeper understanding of the evolutionary and genomic mechanisms underpinning resistance evolution is required to safeguard health and food production. Several studies have implicated transposable elements (TEs) in xenobiotic-resistance evolution in insects. However, analyses are generally restricted to one insect species and/or one or a few xenobiotic gene families (XGFs). We examine evidence for TE accumulation at XGFs by performing a comparative genomic analysis across 20 aphid genomes, considering major subsets of XGFs involved in metabolic resistance to insecticides: cytochrome P450s, glutathione S-transferases, esterases, UDP-glucuronosyltransferases, and ABC transporters. We find that TEs are significantly enriched at XGFs compared with other genes. XGFs show similar levels of TE enrichment to those of housekeeping genes. But unlike housekeeping genes, XGFs are not constitutively expressed in germline cells, supporting the selective enrichment of TEs at XGFs rather than enrichment owing to chromatin availability. Hotspots of extreme TE enrichment occur around certain XGFs. We find, in aphids of agricultural importance, particular enrichment of TEs around cytochrome P450 genes with known functions in the detoxification of synthetic insecticides. Our results provide evidence supporting a general role for TEs as a source of genomic variation at host XGFs and highlight the existence of considerable variability in TE content across XGFs and host species. These findings show the need for detailed functional verification analyses to clarify the significance of individual TE insertions and elucidate underlying mechanisms at TE-XGF hotspots.

Genetic variation in chromatin state across multiple tissues in Drosophila melanogaster

We use ATAC-seq to examine chromatin accessibility for four different tissues in Drosophila melanogaster: adult female brain, ovaries, and both wing and eye-antennal imaginal discs from males. Each tissue is assayed in eight different inbred strain genetic backgrounds, seven associated with a reference quality genome assembly. We develop a method for the quantile normalization of ATAC-seq fragments and test for differences in coverage among genotypes, tissues, and their interaction at 44099 peaks throughout the euchromatic genome. For the strains with reference quality genome assemblies, we correct ATAC-seq profiles for read mis-mapping due to nearby polymorphic structural variants (SVs). Comparing coverage among genotypes without accounting for SVs results in a highly elevated rate (55%) of identifying false positive differences in chromatin state between genotypes. After SV correction, we identify 1050, 30383, and 4508 regions whose peak heights are polymorphic among genotypes, among tissues, or exhibit genotype-by-tissue interactions, respectively. Finally, we identify 3988 candidate causative variants that explain at least 80% of the variance in chromatin state at nearby ATAC-seq peaks.

The role of structural variants in pest adaptation and genome evolution of the Colorado potato beetle, Leptinotarsa decemlineata (Say)

Structural variation has been associated with genetic diversity and adaptation. Despite these observations, it is not clear what their relative importance is for evolution, especially in rapidly adapting species. Here, we examine the significance of structural polymorphisms in pesticide resistance evolution of the agricultural super-pest, the Colorado potato beetle, Leptinotarsa decemlineata. By employing a parent offspring trio sequencing procedure, we develop highly contiguous reference genomes to characterize structural variation. These updated assemblies represent >100-fold improvement of contiguity and include derived pest and ancestral nonpest individuals. We identify >200,000 structural variations, which appear to be nonrandomly distributed across the genome as they co-occur with transposable elements and genes. Structural variations intersect with exons in a large proportion of gene annotations (~20%) that are associated with insecticide resistance (including cytochrome P450s), development, and transcription. To understand the role structural variations play in adaptation, we measure their allele frequencies among an additional 57 individuals using whole genome resequencing data, which represents pest and nonpest populations of North America. Incorporating multiple independent tests to detect the signature of natural selection using SNP data, we identify 14 genes that are probably under positive selection, include structural variations, and SNPs of elevated frequency within the pest lineages. Among these, three are associated with insecticide resistance based on previous research. One of these genes, CYP4g15, is coinduced during insecticide exposure with glycosyltransferase-13, which is a duplicated gene enclosed within a structural variant adjacent to the CYP4g15 genic region. These results demonstrate the significance of structural variations as a genomic feature to describe species history, genetic diversity, and adaptation.

Discovery of malathion resistance QTL in Drosophila melanogaster using a bulked phenotyping approach

Drosophila melanogaster has proved an effective system with which to understand the evolutionary genetics and molecular mechanisms of insecticide resistance. Insecticide use has left signatures of selection in the fly genome, and both functional and quantitative genetic studies in the system have identified genes and variants associated with resistance. Here, we use D. melanogaster and leverage a bulk phenotyping and pooled sequencing "extreme quantitative trait loci" approach to genetically dissect variation in resistance to malathion, an organophosphate insecticide. We resolve 2 quantitative trait loci, one of which implicates allelic variation at the cytochrome P450 gene Cyp6g1, a strong candidate based on previous work. The second shows no overlap with hits from a previous genome-wide association study for malathion resistance, recapitulating other studies showing that different strategies for complex trait dissection in flies can yield apparently different architectures. Notably, we see no genetic signal at the Ace gene. Ace encodes the target of organophosphate insecticide inhibition, and genome-wide association studies have identified strong Ace-linked associations with resistance in flies. The absence of quantitative trait locus implicating Ace here is most likely because our mapping population does not segregate for several of the known functional polymorphisms impacting resistance at Ace, perhaps because our population is derived from flies collected prior to the widespread use of organophosphate insecticides. Our fundamental approach can be an efficient, powerful strategy to dissect genetic variation in resistance traits. Nonetheless, studies seeking to interrogate contemporary insecticide resistance variation may benefit from deriving mapping populations from more recently collected strains.

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.

Genome-Wide Screening of Transposable Elements in the Whitefly, Bemisia tabaci (Hemiptera: Aleyrodidae), Revealed Insertions with Potential Insecticide Resistance Implications

Simple Summary

Transposable elements (TEs) are mobile DNA sequences hosted in the genomes of various organisms. These elements have the ability to mediate regulatory changes, which can result in changes in gene expression. Bemisia tabaci is an important agricultural pest that has been linked to several cases of insecticide resistance. In this study, we conducted a genome-wide screening of TEs in the B. tabaci genome using bioinformatics tools. Results revealed a total of 1,292,393 TE copies clustered into 4872 lineages. The TE insertion site analysis revealed 94 insertions within or near defensome genes.

Abstract

Transposable elements (TEs) are genetically mobile units that move from one site to another within a genome. These units can mediate regulatory changes that can result in massive changes in genes expression. In fact, a precise identification of TEs can allow the detection of the mechanisms involving these elements in gene regulation and genome evolution. In the present study, a genome-wide analysis of the Hemipteran pest Bemisia tabaci was conducted using bioinformatics tools to identify, annotate and estimate the age of TEs, in addition to their insertion sites, within or near of the defensome genes involved in insecticide resistance. Overall, 1,292,393 TE copies were identified in the B. tabaci genome grouped into 4872 lineages. A total of 699 lineages were found to belong to Class I of TEs, 1348 belong to Class II, and 2825 were uncategorized and form the largest part of TEs (28.81%). The TE age estimation revealed that the oldest TEs invasion happened 14 million years ago (MYA) and the most recent occurred 0.2 MYA with the insertion of Class II TE elements. The analysis of TE insertion sites in defensome genes revealed 94 insertions. Six of these TE insertions were found within or near previously identified differentially expressed insecticide resistance genes. These insertions may have a potential role in the observed insecticide resistance in these pests.

Transposable element variants and their potential adaptive impact in urban populations of the malaria vector Anopheles coluzzii

Anopheles coluzzii is one of the primary vectors of human malaria in sub-Saharan Africa. Recently, it has spread into the main cities of Central Africa threatening vector control programs. The adaptation of An. coluzzii to urban environments partly results from an increased tolerance to organic pollution and insecticides. Some of the molecular mechanisms for ecological adaptation are known, but the role of transposable elements (TEs) in the adaptive processes of this species has not been studied yet. As a first step toward assessing the role of TEs in rapid urban adaptation, we sequenced using long reads six An. coluzzii genomes from natural breeding sites in two major Central Africa cities. We de novo annotated TEs in these genomes and in an additional high-quality An. coluzzii genome, and we identified 64 new TE families. TEs were nonrandomly distributed throughout the genome with significant differences in the number of insertions of several superfamilies across the studied genomes. We identified seven putatively active families with insertions near genes with functions related to vectorial capacity, and several TEs that may provide promoter and transcription factor binding sites to insecticide resistance and immune-related genes. Overall, the analysis of multiple high-quality genomes allowed us to generate the most comprehensive TE annotation in this species to date and identify several TE insertions that could potentially impact both genome architecture and the regulation of functionally relevant genes. These results provide a basis for future studies of the impact of TEs on the biology of An. coluzzii.

Molecular innovations underlying resistance to nicotine and neonicotinoids in the aphid Myzus persicae

The green peach aphid, Myzus persicae, is a globally distributed highly damaging crop pest. This species has demonstrated an exceptional ability to evolve resistance to both synthetic insecticides used for control, and natural insecticides produced by certain plants as a chemical defense against insect attack. Here we review work characterizing the evolution of resistance in M. persicae to the natural insecticide nicotine and the structurally related class of synthetic neonicotinoid insecticides. We outline how research on this topic has provided insights into long-standing questions of both evolutionary and applied importance. These include questions pertaining to the origins of novel traits, the number and nature of mutational events or 'adaptive steps' underlying the evolution of new phenotypes, and whether host plant adaptations can be co-opted to confer resistance to synthetic insecticides. Finally, research on the molecular mechanisms underlying insecticide resistance in M. persicae has generated several outstanding questions on the genetic architecture of resistance to both natural and synthetic xenobiotics, and we conclude by identifying key knowledge gaps for future research. © 2021 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Regulatory regions in natural transposable element insertions drive interindividual differences in response to immune challenges in Drosophila

Background

Variation in gene expression underlies interindividual variability in relevant traits including immune response. However, the genetic variation responsible for these gene expression changes remains largely unknown. Among the non-coding variants that could be relevant, transposable element insertions are promising candidates as they have been shown to be a rich and diverse source of cis-regulatory elements.

Results

In this work, we use a population genetics approach to identify transposable element insertions likely to increase the tolerance of Drosophila melanogaster to bacterial infection by affecting the expression of immune-related genes. We identify 12 insertions associated with allele-specific expression changes in immune-related genes. We experimentally validate three of these insertions including one likely to be acting as a silencer, one as an enhancer, and one with a dual role as enhancer and promoter. The direction in the change of gene expression associated with the presence of several of these insertions is consistent with an increased survival to infection. Indeed, for one of the insertions, we show that this is the case by analyzing both natural populations and CRISPR/Cas9 mutants in which the insertion is deleted from its native genomic context.

Conclusions

We show that transposable elements contribute to gene expression variation in response to infection in D. melanogaster and that this variation is likely to affect their survival capacity. Because the role of transposable elements as regulatory elements is not restricted to Drosophila, transposable elements are likely to play a role in immune response in other organisms as well.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-021-02471-3.

Transposable elements expression in Rhinella marina (cane toad) specimens submitted to immune and stress challenge

Transposable elements (TEs) are important components of eukaryotic genomes and compose around 30% of the genome of Rhinella marina, an invasive toad species. Considering the possible role of TEs in the adaptation of populations, we have analyzed the expression of TEs in publicly available spleen tissue transcriptomic data generated for this species after immune and stress challenge. By analyzing the transcriptome assembly, we detected a high number of TE segments. Moreover, some distinct TE families were differentially expressed in some conditions. Our result shows that several TEs are capable of being transcribed in R. marina and they could help to generate a rapid response of specimens to the environment. Also, we can suggest that these TEs could be activated in the germinative cells as well producing variability to be selected and shaped by the evolutionary processes behind the success of this invasive species. Thus, the TEs are important targets for investigation in the context of R. marina adaptation.

Predicted Cellular Interactors of the Endogenous Retrovirus-K Integrase Enzyme

Integrase (IN) enzymes are found in all retroviruses and are crucial in the retroviral integration process. Many studies have revealed how exogenous IN enzymes, such as the human immunodeficiency virus (HIV) IN, contribute to altered cellular function. However, the same consideration has not been given to viral IN originating from symbionts within our own DNA. Endogenous retrovirus-K (ERVK) is pathologically associated with neurological and inflammatory diseases along with several cancers. The ERVK IN interactome is unknown, and the question of how conserved the ERVK IN protein-protein interaction motifs are as compared to other retroviral integrases is addressed in this paper. The ERVK IN protein sequence was analyzed using the Eukaryotic Linear Motif (ELM) database, and the results are compared to ELMs of other betaretroviral INs and similar eukaryotic INs. A list of putative ERVK IN cellular protein interactors was curated from the ELM list and submitted for STRING analysis to generate an ERVK IN interactome. KEGG analysis was used to identify key pathways potentially influenced by ERVK IN. It was determined that the ERVK IN potentially interacts with cellular proteins involved in the DNA damage response (DDR), cell cycle, immunity, inflammation, cell signaling, selective autophagy, and intracellular trafficking. The most prominent pathway identified was viral carcinogenesis, in addition to select cancers, neurological diseases, and diabetic complications. This potentiates the role of ERVK IN in these pathologies via protein-protein interactions facilitating alterations in key disease pathways.

Read between the lines: Diversity of non-translational selection pressures on local codon usage

Protein coding genes can contain specific motifs within their nucleotide sequence that function as a signal for various biological pathways. The presence of such sequence motifs within a gene can have beneficial or detrimental effects on the phenotype and fitness of an organism, and this can lead to the enrichment or avoidance of this sequence motif. The degeneracy of the genetic code allows for the existence of alternative synonymous sequences that exclude or include these motifs, while keeping the encoded amino acid sequence intact. This implies that locally, there can be a selective pressure for preferentially using a codon over its synonymous alternative in order to avoid or enrich a specific sequence motif. This selective pressure could -in addition to mutation, drift and selection for translation efficiency and accuracy- contribute to shape the codon usage bias. In this review, we discuss patterns of avoidance of (or enrichment for) the various biological signals contained in specific nucleotide sequence motifs: transcription and translation initiation and termination signals, mRNA maturation signals, and antiviral immune system targets. Experimental data on the phenotypic or fitness effects of synonymous mutations in these sequence motifs confirm that they can be targets of local selection pressures on codon usage. We also formulate the hypothesis that transposable elements could have a similar impact on codon usage through their preferred integration sequences. Overall, selection on codon usage appears to be a combination of a global selection pressure imposed by the translation machinery, and a patchwork of local selection pressures related to biological signals contained in specific sequence motifs.

Transposable Element Landscape in Drosophila Populations Selected for Longevity

Transposable elements (TEs) inflict numerous negative effects on health and fitness as they replicate by integrating into new regions of the host genome. Even though organisms employ powerful mechanisms to demobilize TEs, transposons gradually lose repression during aging. The rising TE activity causes genomic instability and was implicated in age-dependent neurodegenerative diseases, inflammation, and the determination of lifespan. It is therefore conceivable that long-lived individuals have improved TE silencing mechanisms resulting in reduced TE expression relative to their shorter-lived counterparts and fewer genomic insertions. Here, we test this hypothesis by performing the first genome-wide analysis of TE insertions and expression in populations of Drosophila melanogaster selected for longevity through late-life reproduction for 50–170 generations from four independent studies. Contrary to our expectation, TE families were generally more abundant in long-lived populations compared with nonselected controls. Although simulations showed that this was not expected under neutrality, we found little evidence for selection driving TE abundance differences. Additional RNA-seq analysis revealed a tendency for reducing TE expression in selected populations, which might be more important for lifespan than regulating genomic insertions. We further find limited evidence of parallel selection on genes related to TE regulation and transposition. However, telomeric TEs were genomically and transcriptionally more abundant in long-lived flies, suggesting improved telomere maintenance as a promising TE-mediated mechanism for prolonging lifespan. Our results provide a novel viewpoint indicating that reproduction at old age increases the opportunity of TEs to be passed on to the next generation with little impact on longevity.

New Environment, New Invaders-Repeated Horizontal Transfer of LINEs to Sea Snakes

Although numerous studies have found horizontal transposon transfer (HTT) to be widespread across metazoans, few have focused on HTT in marine ecosystems. To investigate potential recent HTTs into marine species, we searched for novel repetitive elements in sea snakes, a group of elapids which transitioned to a marine habitat at most 18 Ma. Our analysis uncovered repeated HTTs into sea snakes following their marine transition. The seven subfamilies of horizontally transferred LINE retrotransposons we identified in the olive sea snake (Aipysurus laevis) are transcribed, and hence are likely still active and expanding across the genome. A search of 600 metazoan genomes found all seven were absent from other amniotes, including terrestrial elapids, with the most similar LINEs present in fish and marine invertebrates. The one exception was a similar LINE found in sea kraits, a lineage of amphibious elapids which independently transitioned to a marine environment 25 Ma. Our finding of repeated horizontal transfer events into marine snakes greatly expands past findings that the marine environment promotes the transfer of transposons. Transposons are drivers of evolution as sources of genomic sequence and hence genomic novelty. We identified 13 candidate genes for HTT-induced adaptive change based on internal or neighboring HTT LINE insertions. One of these, ADCY4, is of particular interest as a part of the KEGG adaptation pathway “Circadian Entrainment.” This provides evidence of the ecological interactions between species influencing evolution of metazoans not only through specific selection pressures, but also by contributing novel genomic material.

Crossroads between transposons and gene regulation

Transposons are mobile genetic elements that have made a large contribution to genome evolution in a largely species-specific manner. A wide variety of different transposons have invaded genomes throughout evolution, acting in a first instance as 'selfish' elements, whose success was determined by their ability to self-replicate and expand within the host genome. However, their coevolution with the host has created many crossroads between transposons and the regulation of host gene expression. Transposons are an abundant source of transcriptional modulatory elements, such as gene promoters and enhancers, splicing and termination sites, and regulatory non-coding RNAs. Moreover, transposons have driven the evolution of host defence mechanisms that have been repurposed for gene regulation. However, dissecting the potential functional roles of transposons remains challenging owing to their evolutionary path, as well as their repetitive nature, which requires the development of specialized analytical tools. In this special issue, we present a collection of articles that lay out current paradigms in the field and discuss a vision for future research. This article is part of a discussion meeting issue 'Crossroads between transposons and gene regulation'.