A survey of transposable element classification systems--a call for a fundamental update to meet the challenge of their diversity and complexity

Systematic identification of cargo-mobilizing genetic elements reveals new dimensions of eukaryotic diversity

Cargo-mobilizing mobile elements (CMEs) are genetic entities that faithfully transpose diverse protein coding sequences. Although common in bacteria, we know little about eukaryotic CMEs because no appropriate tools exist for their annotation. For example, Starships are giant fungal CMEs whose functions are largely unknown because they require time-intensive manual curation. To address this knowledge gap, we developed starfish, a computational workflow for high-throughput eukaryotic CME annotation. We applied starfish to 2 899 genomes of 1 649 fungal species and found that starfish recovers known Starships with 95% combined precision and recall while expanding the number of annotated elements ten-fold. Extant Starship diversity is partitioned into 11 families that differ in their enrichment patterns across fungal classes. Starship cargo changes rapidly such that elements from the same family differ substantially in their functional repertoires, which are predicted to contribute to diverse biological processes such as metabolism. Many elements have convergently evolved to insert into 5S rDNA and AT-rich sequence while others integrate into random locations, revealing both specialist and generalist strategies for persistence. Our work establishes a framework for advancing mobile element biology and provides the means to investigate an emerging dimension of eukaryotic genetic diversity, that of genomes within genomes.

FishTEDB 2.0: an update fish transposable element (TE) database with new functions to facilitate TE research

We launched the initial version of FishTEDB in 2018, which aimed to establish an open-source, user-friendly, data-rich transposable element (TE) database. Over the past 5 years, FishTEDB 1.0 has gained approximately 10 000 users, accumulating more than 450 000 interactions. With the unveiling of extensive fish genome data and the increasing emphasis on TE research, FishTEDB needs to extend the richness of data and functions. To achieve the above goals, we introduced 33 new fish species to FishTEDB 2.0, encompassing a wide array of fish belonging to 48 orders. To make the updated database more functional, we added a genome browser to visualize the positional relationship between TEs and genes and the estimated TE insertion time in different species. In conclusion, we released a new version of the fish TE database, FishTEDB 2.0, designed to assist researchers in the future study of TE functions and promote the progress of biological theories related to TEs. Database URL: https://www.fishtedb.com/.

A metagenomics pipeline reveals insertion sequence-driven evolution of the microbiota

Insertion sequence (IS) elements are mobile genetic elements in bacterial genomes that support adaptation. We developed a database of IS elements coupled to a computational pipeline that identifies IS element insertions in the microbiota. We discovered that diverse IS elements insert into the genomes of intestinal bacteria regardless of human host lifestyle. These insertions target bacterial accessory genes that aid in their adaptation to unique environmental conditions. Using IS expansion in Bacteroides, we show that IS activity leads to the insertion of "hot spots" in accessory genes. We show that IS insertions are stable and can be transferred between humans. Extreme environmental perturbations force IS elements to fall out of the microbiota, and many fail to rebound following homeostasis. Our work shows that IS elements drive bacterial genome diversification within the microbiota and establishes a framework for understanding how strain-level variation within the microbiota impacts human health.

Mobile genetic element-based gene editing and genome engineering: Recent advances and applications

Genome engineering has revolutionized several scientific fields, ranging from biochemistry and fundamental research to therapeutic uses and crop development. Diverse engineering toolkits have been developed and used to effectively modify the genome sequences of organisms. However, there is a lack of extensive reviews on genome engineering technologies based on mobile genetic elements (MGEs), which induce genetic diversity within host cells by changing their locations in the genome. This review provides a comprehensive update on the versatility of MGEs as powerful genome engineering tools that offers efficient solutions to challenges associated with genome engineering. MGEs, including DNA transposons, retrotransposons, retrons, and CRISPR-associated transposons, offer various advantages, such as a broad host range, genome-wide mutagenesis, efficient large-size DNA integration, multiplexing capabilities, and in situ single-stranded DNA generation. We focused on the components, mechanisms, and features of each MGE-based tool to highlight their cellular applications. Finally, we discussed the current challenges of MGE-based genome engineering and provided insights into the evolving landscape of this transformative technology. In conclusion, the combination of genome engineering with MGE demonstrates remarkable potential for addressing various challenges and advancing the field of genetic manipulation, and promises to revolutionize our ability to engineer and understand the genomes of diverse organisms.

Repeat turnover meets stable chromosomes: repetitive DNA sequences mark speciation and gene pool boundaries in sugar beet and wild beets

Sugar beet and its wild relatives share a base chromosome number of nine and similar chromosome morphologies. Yet, interspecific breeding is impeded by chromosome and sequence divergence that is still not fully understood. Since repetitive DNAs are among the fastest evolving parts of the genome, we investigated, if repeatome innovations and losses are linked to chromosomal differentiation and speciation. We traced genome and chromosome-wide evolution across 13 beet species comprising all sections of the genera Beta and Patellifolia. For this, we combined short and long read sequencing, flow cytometry, and cytogenetics to build a comprehensive framework that spans the complete scale from DNA to chromosome to genome. Genome sizes and repeat profiles reflect the separation into three gene pools with contrasting evolutionary patterns. Among all repeats, satellite DNAs harbor most genomic variability, leading to fundamentally different centromere architectures, ranging from chromosomal uniformity in Beta and Patellifolia to the formation of patchwork chromosomes in Corollinae/Nanae. We show that repetitive DNAs are causal for the genome expansions and contractions across the beet genera, providing insights into the genomic underpinnings of beet speciation. Satellite DNAs in particular vary considerably between beet genomes, leading to the evolution of distinct chromosomal setups in the three gene pools, likely contributing to the barriers in beet breeding. Thus, with their isokaryotypic chromosome sets, beet genomes present an ideal system for studying the link between repeats, genomic variability, and chromosomal differentiation and provide a theoretical fundament for understanding barriers in any crop breeding effort.

Orthoptera-TElib: a library of Orthoptera transposable elements for TE annotation

Transposable elements (TEs) are a major component of eukaryotic genomes and are present in almost all eukaryotic organisms. TEs are highly dynamic between and within species, which significantly affects the general applicability of the TE databases. Orthoptera is the only known group in the class Insecta with a significantly enlarged genome (0.93-21.48 Gb). When analyzing the large genome using the existing TE public database, the efficiency of TE annotation is not satisfactory. To address this limitation, it becomes imperative to continually update the available TE resource library and the need for an Orthoptera-specific library as more insect genomes are publicly available. Here, we used the complete genome data of 12 Orthoptera species to de novo annotate TEs, then manually re-annotate the unclassified TEs to construct a non-redundant Orthoptera-specific TE library: Orthoptera-TElib. Orthoptera-TElib contains 24,021 TE entries including the re-annotated results of 13,964 unknown TEs. The naming of TE entries in Orthoptera-TElib adopts the same naming as RepeatMasker and Dfam and is encoded as the three-level form of "level1/level2-level3". Orthoptera-TElib can be directly used as an input reference database and is compatible with mainstream repetitive sequence analysis software such as RepeatMasker and dnaPipeTE. When analyzing TEs of Orthoptera species, Orthoptera-TElib performs better TE annotation as compared to Dfam and Repbase regardless of using low-coverage sequencing or genome assembly data. The most improved TE annotation result is Angaracris rhodopa, which has increased from 7.89% of the genome to 53.28%. Finally, Orthoptera-TElib is stored in Sqlite3 for the convenience of data updates and user access.

A metagenomics pipeline reveals insertion sequence-driven evolution of the microbiota

Insertion sequence (IS) elements are mobile genetic elements in bacterial genomes that support adaptation. We developed a database of IS elements coupled to a computational pipeline that identifies IS element insertions in the microbiota. We discovered that diverse IS elements insert into the genomes of intestinal bacteria regardless of human host lifestyle. These insertions target bacterial accessory genes that aid in their adaptation to unique environmental conditions. Using IS expansion in Bacteroides, we show that IS activity leads to insertion "hot spots" in accessory genes. We show that IS insertions are stable and can be transferred between humans. Extreme environmental perturbations force IS elements to fall out of the microbiota and many fail to rebound following homeostasis. Our work shows that IS elements drive bacterial genome diversification within the microbiota and establishes a framework for understanding how strain level variation within the microbiota impacts human health.

Evolutionary dynamics of the LTR-retrotransposon crapaud in the Podospora anserina species complex and the interaction with repeat-induced point mutations

BACKGROUND

The genome of the filamentous ascomycete Podospora anserina shows a relatively high abundance of retrotransposons compared to other interspersed repeats. The LTR-retrotransposon family crapaud is particularly abundant in the genome, and consists of multiple diverged sequence variations specifically localized in the 5' half of both long terminal repeats (LTRs). P. anserina is part of a recently diverged species-complex, which makes the system ideal to classify the crapaud family based on the observed LTR variation and to study the evolutionary dynamics, such as the diversification and bursts of the elements over recent evolutionary time.

RESULTS

We developed a sequence similarity network approach to classify the crapaud repeats of seven genomes representing the P. anserina species complex into 14 subfamilies. This method does not utilize a consensus sequence, but instead it connects any copies that share enough sequence similarity over a set sequence coverage. Based on phylogenetic analyses, we found that the crapaud repeats likely diversified in the ancestor of the complex and have had activity at different time points for different subfamilies. Furthermore, while we hypothesized that the evolution into multiple subfamilies could have been a direct effect of escaping the genome defense system of repeat induced point mutations, we found this not to be the case.

CONCLUSIONS

Our study contributes to the development of methods to classify transposable elements in fungi, and also highlights the intricate patterns of retrotransposon evolution over short timescales and under high mutational load caused by nucleotide-altering genome defense.

Selective control of transposable element expression during T cell exhaustion and anti-PD-1 treatment

In chronic infections and cancer, T cells are exposed to prolonged antigen stimulation, resulting in loss of function (or exhaustion) and impairment of effective immunological protection. Exhausted T cells are heterogeneous and include early progenitors (Tpex) and terminally exhausted cells (Tex). Here, we used bulk and single-cell transcriptomics to analyze expression of transposable elements (TEs) in subpopulations of mouse and human CD8+ tumor-infiltrating T lymphocytes (TILs). We show that in mice, members of the virus-like murine VL30 TE family (mostly intact, evolutionary young ERV1s) are strongly repressed in terminally exhausted CD8+ T cells in both tumor and viral models of exhaustion. Tpex expression of these VL30s, which are mainly intergenic and transcribed independently of their closest gene neighbors, was driven by Fli1, a transcription factor involved in progression from Tpex to Tex. Immune checkpoint blockade (ICB) in both mice and patients with cancer increased TE expression (including VL30 in mice), demonstrating that TEs may be applicable as ICB response biomarkers. We conclude that expression of TEs is tightly regulated in TILs during establishment of exhaustion and reprogramming by ICB. Analyses of TE expression on single cells and bulk populations open opportunities for understanding immune cell identity and heterogeneity, as well as for defining cellular gene expression signatures and disease biomarkers.

Analyses of 600+ insect genomes reveal repetitive element dynamics and highlight biodiversity-scale repeat annotation challenges

Repetitive elements (REs) are integral to the composition, structure, and function of eukaryotic genomes, yet remain understudied in most taxonomic groups. We investigated REs across 601 insect species and report wide variation in RE dynamics across groups. Analysis of associations between REs and protein-coding genes revealed dynamic evolution at the interface between REs and coding regions across insects, including notably elevated RE-gene associations in lineages with abundant long interspersed nuclear elements (LINEs). We leveraged this large, empirical data set to quantify impacts of long-read technology on RE detection and investigate fundamental challenges to RE annotation in diverse groups. In long-read assemblies, we detected ∼36% more REs than short-read assemblies, with long terminal repeats (LTRs) showing 162% increased detection, whereas DNA transposons and LINEs showed less respective technology-related bias. In most insect lineages, 25%-85% of repetitive sequences were "unclassified" following automated annotation, compared with only ∼13% in Drosophila species. Although the diversity of available insect genomes has rapidly expanded, we show the rate of community contributions to RE databases has not kept pace, preventing efficient annotation and high-resolution study of REs in most groups. We highlight the tremendous opportunity and need for the biodiversity genomics field to embrace REs and suggest collective steps for making progress toward this goal.

From parasites to partners: exploring the intricacies of host-transposon dynamics and coevolution

Transposable elements, often referred to as "jumping genes," have long been recognized as genomic parasites due to their ability to integrate and disrupt normal gene function and induce extensive genomic alterations, thereby compromising the host's fitness. To counteract this, the host has evolved a plethora of mechanisms to suppress the activity of the transposons. Recent research has unveiled the host-transposon relationships to be nuanced and complex phenomena, resulting in the coevolution of both entities. Transposition increases the mutational rate in the host genome, often triggering physiological pathways such as immune and stress responses. Current gene transfer technologies utilizing transposable elements have potential drawbacks, including off-target integration, induction of mutations, and modifications of cellular machinery, which makes an in-depth understanding of the host-transposon relationship imperative. This review highlights the dynamic interplay between the host and transposable elements, encompassing various factors and components of the cellular machinery. We provide a comprehensive discussion of the strategies employed by transposable elements for their propagation, as well as the mechanisms utilized by the host to mitigate their parasitic effects. Additionally, we present an overview of recent research identifying host proteins that act as facilitators or inhibitors of transposition. We further discuss the evolutionary outcomes resulting from the genetic interactions between the host and the transposable elements. Finally, we pose open questions in this field and suggest potential avenues for future research.

Comparative genomics reveals insights into anuran genome size evolution

BACKGROUND

Amphibians, particularly anurans, display an enormous variation in genome size. Due to the unavailability of whole genome datasets in the past, the genomic elements and evolutionary causes of anuran genome size variation are poorly understood. To address this, we analyzed whole-genome sequences of 14 anuran species ranging in size from 1.1 to 6.8 Gb. By annotating multiple genomic elements, we investigated the genomic correlates of anuran genome size variation and further examined whether the genome size relates to habitat types.

RESULTS

Our results showed that intron expansions or contraction and Transposable Elements (TEs) diversity do not contribute significantly to genome size variation. However, the recent accumulation of transposable elements (TEs) and the lack of deletion of ancient TEs primarily accounted for the evolution of anuran genome sizes. Our study showed that the abundance and density of simple repeat sequences positively correlate with genome size. Ancestral state reconstruction revealed that genome size exhibits a taxon-specific pattern of evolution, with families Bufonidae and Pipidae experiencing extreme genome expansion and contraction events, respectively. Our result showed no relationship between genome size and habitat types, although large genome-sized species are predominantly found in humid habitats.

CONCLUSIONS

Overall, our study identified the genomic element and their evolutionary dynamics accounting for anuran genome size variation, thus paving a path to a greater understanding of the size evolution of the genome in amphibians.

Transposable elements and their role in aging

Transposable elements (TEs) are an important part of eukaryotic genomes. The role of somatic transposition in aging, carcinogenesis, and other age-related diseases has been determined. This review discusses the fundamental properties of TEs and their complex interactions with cellular processes, which are crucial for understanding the diverse effects of their activity on the genetics and epigenetics of the organism. The interactions of TEs with recombination, replication, repair, and chromosomal regulation; the ability of TEs to maintain a balance between their own activity and repression, the involvement of TEs in the creation of new or alternative genes, the expression of coding/non-coding RNA, and the role in DNA damage and modification of regulatory networks are reviewed. The contribution of the derepressed TEs to age-dependent effects in individual cells/tissues in different organisms was assessed. Conflicting information about TE activity under stress as well as theories of aging mechanisms related to TEs is discussed. On the one hand, transposition activity in response to stressors can lead to organisms acquiring adaptive innovations of great importance for evolution at the population level. On the other hand, the TE expression can cause decreased longevity and stress tolerance at the individual level. The specific features of TE effects on aging processes in germline and soma and the ways of their regulation in cells are highlighted. Recent results considering somatic mutations in normal human and animal tissues are indicated, with the emphasis on their possible functional consequences. In the context of aging, the correlation between somatic TE activation and age-related changes in the number of proteins required for heterochromatin maintenance and longevity regulation was analyzed. One of the original features of this review is a discussion of not only effects based on the TEs insertions and the associated consequences for the germline cell dynamics and somatic genome, but also the differences between transposon- and retrotransposon-mediated structural genome changes and possible phenotypic characteristics associated with aging and various age-related pathologies. Based on the analysis of published data, a hypothesis about the influence of the species-specific features of number, composition, and distribution of TEs on aging dynamics of different animal genomes was formulated.

Transposable elements in the transcriptome of the velvetbean caterpillar Anticarsia gemmatalis Hübner, 1818 (Lepidoptera: Erebidae)

Transposable elements (TEs) are DNA sequences that possess the ability to move from one genomic location to another. These sequences contribute to a significant fraction of the genomes of most eukaryotes and can impact their architecture and regulation. In this paper, we present the first data related to the identification and characterization of TEs present in the transcriptome of Anticarsia gemmatalis. Approximately, 835 transcripts showed significant similarity to TEs and (or) characteristic domains. Retrotransposons accounted for 71.2% (595 sequences) of the identified elements, while DNA transposons were less abundant, with 240 annotations (28.8%). TEs were classified into 30 superfamilies, with SINE3/5S and Gypsy being the most abundant. Based on the sequences of TEs found in the transcriptome, we were able to locate conserved regions in the chromosomes of this species. The analysis of differential expression of TEs in susceptible and resistant strains, challenged and not challenged with Bacillus thuringiensis ( Bt) from in silico analysis, indicated that exposure to Bt can regulate the transcription of mobile genetic elements in the velvetbean caterpillar. Thus, these data contribute significantly to the knowledge of the structure and composition of these elements in the genome of this species, and suggest the role of stress on their expression.

Transposable Elements Co-Option in Genome Evolution and Gene Regulation

The genome is no longer deemed as a fixed and inert item but rather as a moldable matter that is continuously evolving and adapting. Within this frame, Transposable Elements (TEs), ubiquitous, mobile, repetitive elements, are considered an alive portion of the genomes to date, whose functions, although long considered "dark", are now coming to light. Here we will review that, besides the detrimental effects that TE mobilization can induce, TEs have shaped genomes in their current form, promoting genome sizing, genomic rearrangements and shuffling of DNA sequences. Although TEs are mostly represented in the genomes by evolutionarily old, short, degenerated, and sedentary fossils, they have been thoroughly co-opted by the hosts as a prolific and original source of regulatory instruments for the control of gene transcription and genome organization in the nuclear space. For these reasons, the deregulation of TE expression and/or activity is implicated in the onset and progression of several diseases. It is likely that we have just revealed the outermost layers of TE functions. Further studies on this portion of the genome are required to unlock novel regulatory functions that could also be exploited for diagnostic and therapeutic approaches.

Chemotherapy Drugs Act Differently in the Expression and Somatic Mobilization of the mariner Transposable Element in Drosophila simulans

Transposable elements (TEs) are abundant in genomes. Their mobilization can lead to genetic variability that is useful for evolution, but can also have deleterious biological effects. Somatic mobilization (SM) has been linked to degenerative diseases, such as Alzheimer's disease and cancer. We used a Drosophila simulans strain, in which SM can be measured by counting red spots in the eyes, to investigate how chemotherapeutic agents affect expression and SM of the mariner TE. Flies were treated with Cisplatin, Dacarbazine, and Daunorubicin. After acute exposure, relative expression of mariner was quantified by RT-qPCR and oxidative stress was measured by biochemical assays. Exposure to 50 and 100 µg/mL Cisplatin increased mariner expression and ROS levels; catalase activity increased at 100 µg/mL. With chronic exposure, the number of spots also increased, indicating higher mariner SM. Dacarbazine (50 and 100 µg/mL) did not significantly alter mariner expression or mobilization or ROS levels, but decreased catalase activity (100 µg/mL). Daunorubicin (25 and 50 µM) increased mariner expression, but decreased mariner SM. ROS and catalase activity were also reduced. Our data suggest that stress factors may differentially affect the expression and SM of TEs. The increase in mariner transposase gene expression is necessary, but not sufficient for mariner SM.

Variation in the Number and Position of rDNA Loci Contributes to the Diversification and Speciation in Nigella (Ranunculaceae)

Nigella is a small genus belonging to the Ranunculaceae family which is presumably originated and distributed in Aegean and the adjacent Western-Irano-Turanian region. Comparative repeat analysis of N. sativa, N. damascena and N. bucharica was performed using low-pass Illumina genomic reads followed by karyotyping and FISH mapping of seven Nigella species using the in silico identified repeats and ribosomal DNA (rDNA) probes. High- and moderate-copy repeat sequences occupy 57.52, 59.01, and 64.73% of N. sativa, N. damascena and N. bucharica genomes, respectively. Roughly, half of the genomes are retrotransposons (class I transposons), while DNA transposons (class II transposons) contributed to only about 2% of the genomes. The analyzed Nigella species possess large genomes of about 7.4 to 12.4 Gbp/1C. Only two satellite repeats in N. sativa, one in N. damascena and four in N. bucharica were identified, which were mostly (peri)centromeric and represented about 1% of each genome. A high variation in number and position of 45S rDNA loci were found among Nigella species. Interestingly, in N. hispanica, each chromosome revealed at least one 45S rDNA site and one of them occurs in hemizygous condition. Based on the chromosome numbers, genome size and (peri)centromeric satellites, three karyotype groups were observed: Two with 2n = 2x = 12 and a karyotype formula of 10m + 2t (including N. sativa, N. arvensis, N. hispanica as the first group and N. damascena and N. orientalis as the second group) and a more distant group with 2n = 2x = 14 and a karyotype formula of 8m + 2st + 4t (including N. integrifolia and N. bucharica). These karyotype groups agreed with the phylogenetic analysis using ITS and rbcL sequences. We conclude that variation in (peri)centromeric sequences, number and localization of rDNA sites as well as chromosome number (dysploidy) are involved in the diversification of the genus Nigella.

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.

From telomere to telomere: The transcriptional and epigenetic state of human repeat elements

Mobile elements and repetitive genomic regions are sources of lineage-specific genomic innovation and uniquely fingerprint individual genomes. Comprehensive analyses of such repeat elements, including those found in more complex regions of the genome, require a complete, linear genome assembly. We present a de novo repeat discovery and annotation of the T2T-CHM13 human reference genome. We identified previously unknown satellite arrays, expanded the catalog of variants and families for repeats and mobile elements, characterized classes of complex composite repeats, and located retroelement transduction events. We detected nascent transcription and delineated CpG methylation profiles to define the structure of transcriptionally active retroelements in humans, including those in centromeres. These data expand our insight into the diversity, distribution, and evolution of repetitive regions that have shaped the human genome.

A beginner's guide to manual curation of transposable elements

Background

In the study of transposable elements (TEs), the generation of a high confidence set of consensus sequences that represent the diversity of TEs found in a given genome is a key step in the path to investigate these fascinating genomic elements. Many algorithms and pipelines are available to automatically identify putative TE families present in a genome. Despite the availability of these valuable resources, producing a library of high-quality full-length TE consensus sequences largely remains a process of manual curation. This know-how is often passed on from mentor-to-mentee within research groups, making it difficult for those outside the field to access this highly specialised skill.

Results

Our manuscript attempts to fill this gap by providing a set of detailed computer protocols, software recommendations and video tutorials for those aiming to manually curate TEs. Detailed step-by-step protocols, aimed at the complete beginner, are presented in the Supplementary Methods.

Conclusions

The proposed set of programs and tools presented here will make the process of manual curation achievable and amenable to all researchers and in special to those new to the field of TEs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13100-021-00259-7.

Epigenetic influences of mobile genetic elements on ciliate genome architecture and evolution

Mobile genetic elements (MGEs) are transient genetic material that can move either within a single organism's genome or between individuals or species. While historically considered 'junk' DNA (i.e. deleterious or at best neutral), more recent studies reveal the adaptive advantages MGEs provide in lineages across the tree of life. Ciliates, a group of single-celled microbial eukaryotes characterized by nuclear dimorphism, exemplify how epigenetic influences from MGEs shape genome architecture and patterns of molecular evolution. Ciliate nuclear dimorphism may have evolved as a response to transposon invasion and ciliates have since co-opted transposons to carry out programmed DNA deletion. Another example of the effect of MGEs is in providing mechanisms for lateral gene transfer from bacteria, which introduces genetic diversity and, in several cases, drives ecological specialization in ciliates. As a third example, the integration of viral DNA, likely through transduction, provides new genetic material and can change the way host cells defend themselves against other viral pathogens. We argue that the acquisition of MGEs through non-Mendelian patterns of inheritance, coupled with their effects on ciliate genome architecture and expression and persistence throughout evolutionary history, exemplify how the transmission of mobile elements should be considered a mechanism of transgenerational epigenetic inheritance.

inPOSE: A Flexible Toolbox for Chromosomal Cloning and Amplification of Bacterial Transgenes

Cloning the genes and operons encoding heterologous functions in bacterial hosts is now almost exclusively carried out using plasmid vectors. This has multiple drawbacks, including the need for constant selection and variation in copy numbers. The chromosomal integration of transgenes has always offered a viable alternative; however, to date, it has been of limited use due to its tedious nature and often being limited to a single copy. We introduce here a strategy that uses bacterial insertion sequences, which are the simplest autonomous transposable elements to insert and amplify genetic cargo into a bacterial chromosome. Transgene insertion can take place either as transposition or homologous recombination, and copy number amplification is achieved using controlled copy-paste transposition. We display the successful use of IS1 and IS3 for this purpose in Escherichia coli cells using various selection markers. We demonstrate the insertion of selectable genes, an unselectable gene and a five-gene operon in up to two copies in a single step. We continue with the amplification of the inserted cassette to double-digit copy numbers within two rounds of transposase induction and selection. Finally, we analyze the stability of the cloned genetic constructs in the lack of selection and find it to be superior to all investigated plasmid-based systems. Due to the ubiquitous nature of transposable elements, we believe that with proper design, this strategy can be adapted to numerous other bacterial species.

Genome-wide characterization of Mariner-like transposons and their derived MITEs in the Whitefly Bemisia tabaci (Hemiptera: Aleyrodidae)

The whitefly, Bemisia tabaci is a hemipteran pest of vegetable crops vectoring a broad category of viruses. Currently, this insect pest showed a high adaptability and resistance to almost all the chemical compounds commonly used for its control. In many cases, transposable elements (TEs) contributed to the evolution of host genomic plasticity. This study focuses on the annotation of Mariner-like elements (MLEs) and their derived Miniature Inverted repeat Transposable Elements (MITEs) in the genome of B. tabaci. Two full-length MLEs belonging to mauritiana and irritans subfamilies were detected and named Btmar1.1 and Btmar2.1, respectively. Additionally, 548 defective MLE sequences clustering mainly into 19 different Mariner lineages of mauritiana and irritans subfamilies were identified. Each subfamily showed a significant variation in MLE copy number and size. Furthermore, 71 MITEs were identified as MLEs derivatives that could be mobilized via the potentially active transposases encoded by Btmar 1.1 and Btmar2.1. The vast majority of sequences detected in the whitefly genome present unusual terminal inverted repeats (TIRs) of up to 400 bp in length. However, some exceptions are sequences without TIRs. This feature of the MLEs and their derived MITEs in B. tabaci genome that distinguishes them from all the other MLEs so far described in insects, which have TIRs size ranging from 20 to 40 bp. Overall, our study provides an overview of MLEs, especially those with large TIRs, and their related MITEs, as well as diversity of their families, which will provide a better understanding of the evolution and adaptation of the whitefly genome.

Transposable elements and Piwi‑interacting RNAs in hemato‑oncology with a focus on myelodysplastic syndrome (Review)

Our current understanding of hematopoietic stem cell differentiation and the abnormalities that lead to leukemogenesis originates from the accumulation of knowledge regarding protein‑coding genes. However, the possible impact of transposable element (TE) mobilization and the expression of P‑element‑induced WImpy testis‑interacting RNAs (piRNAs) on leukemogenesis has been beyond the scope of scientific interest to date. The expression profiles of these molecules and their importance for human health have only been characterized recently due to the rapid progress of high‑throughput sequencing technology development. In the present review, current knowledge on the expression profile and function of TEs and piRNAs was summarized, with specific focus on their reported involvement in leukemogenesis and pathogenesis of myelodysplastic syndrome.

Engines of change: Transposable element mutation rates are high and variable within Daphnia magna

Transposable elements (TEs) represent a major portion of most eukaryotic genomes, yet little is known about their mutation rates or how their activity is shaped by other evolutionary forces. Here, we compare short- and long-term patterns of genome-wide mutation accumulation (MA) of TEs among 9 genotypes from three populations of Daphnia magna from across a latitudinal gradient. While the overall proportion of the genome comprised of TEs is highly similar among genotypes from Finland, Germany, and Israel, populations are distinguishable based on patterns of insertion site polymorphism. Our direct rate estimates indicate TE movement is highly variable (net rates ranging from -11.98 to 12.79 x 10-5 per copy per generation among genotypes), differing both among populations and TE families. Although gains outnumber losses when selection is minimized, both types of events appear to be highly deleterious based on their low frequency in control lines where propagation is not limited to random, single-progeny descent. With rate estimates 4 orders of magnitude higher than base substitutions, TEs clearly represent a highly mutagenic force in the genome. Quantifying patterns of intra- and interspecific variation in TE mobility with and without selection provides insight into a powerful mechanism generating genetic variation in the genome.

The Dynamism of Transposon Methylation for Plant Development and Stress Adaptation

Plant development processes are regulated by epigenetic alterations that shape nuclear structure, gene expression, and phenotypic plasticity; these alterations can provide the plant with protection from environmental stresses. During plant growth and development, these processes play a significant role in regulating gene expression to remodel chromatin structure. These epigenetic alterations are mainly regulated by transposable elements (TEs) whose abundance in plant genomes results in their interaction with genomes. Thus, TEs are the main source of epigenetic changes and form a substantial part of the plant genome. Furthermore, TEs can be activated under stress conditions, and activated elements cause mutagenic effects and substantial genetic variability. This introduces novel gene functions and structural variation in the insertion sites and primarily contributes to epigenetic modifications. Altogether, these modifications indirectly or directly provide the ability to withstand environmental stresses. In recent years, many studies have shown that TE methylation plays a major role in the evolution of the plant genome through epigenetic process that regulate gene imprinting, thereby upholding genome stability. The induced genetic rearrangements and insertions of mobile genetic elements in regions of active euchromatin contribute to genome alteration, leading to genomic stress. These TE-mediated epigenetic modifications lead to phenotypic diversity, genetic variation, and environmental stress tolerance. Thus, TE methylation is essential for plant evolution and stress adaptation, and TEs hold a relevant military position in the plant genome. High-throughput techniques have greatly advanced the understanding of TE-mediated gene expression and its associations with genome methylation and suggest that controlled mobilization of TEs could be used for crop breeding. However, development application in this area has been limited, and an integrated view of TE function and subsequent processes is lacking. In this review, we explore the enormous diversity and likely functions of the TE repertoire in adaptive evolution and discuss some recent examples of how TEs impact gene expression in plant development and stress adaptation.

The genomic ecosystem of transposable elements in maize

Transposable elements (TEs) constitute the majority of flowering plant DNA, reflecting their tremendous success in subverting, avoiding, and surviving the defenses of their host genomes to ensure their selfish replication. More than 85% of the sequence of the maize genome can be ascribed to past transposition, providing a major contribution to the structure of the genome. Evidence from individual loci has informed our understanding of how transposition has shaped the genome, and a number of individual TE insertions have been causally linked to dramatic phenotypic changes. Genome-wide analyses in maize and other taxa have frequently represented TEs as a relatively homogeneous class of fragmentary relics of past transposition, obscuring their evolutionary history and interaction with their host genome. Using an updated annotation of structurally intact TEs in the maize reference genome, we investigate the family-level dynamics of TEs in maize. Integrating a variety of data, from descriptors of individual TEs like coding capacity, expression, and methylation, as well as similar features of the sequence they inserted into, we model the relationship between attributes of the genomic environment and the survival of TE copies and families. In contrast to the wholesale relegation of all TEs to a single category of junk DNA, these differences reveal a diversity of survival strategies of TE families. Together these generate a rich ecology of the genome, with each TE family representing the evolution of a distinct ecological niche. We conclude that while the impact of transposition is highly family- and context-dependent, a family-level understanding of the ecology of TEs in the genome can refine our ability to predict the role of TEs in generating genetic and phenotypic diversity.

Transposable Elements Contribute to Genome Dynamics and Gene Expression Variation in the Fungal Plant Pathogen Verticillium dahliae

Transposable elements (TEs) are a major source of genetic and regulatory variation in their host genome and are consequently thought to play important roles in evolution. Many fungal and oomycete plant pathogens have evolved dynamic and TE-rich genomic regions containing genes that are implicated in host colonization and adaptation. TEs embedded in these regions have typically been thought to accelerate the evolution of these genomic compartments, but little is known about their dynamics in strains that harbor them. Here, we used whole-genome sequencing data of 42 strains of the fungal plant pathogen Verticillium dahliae to systematically identify polymorphic TEs that may be implicated in genomic as well as in gene expression variation. We identified 2,523 TE polymorphisms and characterize a subset of 8% of the TEs as polymorphic elements that are evolutionary younger, less methylated, and more highly expressed when compared with the remaining 92% of the total TE complement. As expected, the polyrmorphic TEs are enriched in the adaptive genomic regions. Besides, we observed an association of polymorphic TEs with pathogenicity-related genes that localize nearby and that display high expression levels. Collectively, our analyses demonstrate that TE dynamics in V. dahliae contributes to genomic variation, correlates with expression of pathogenicity-related genes, and potentially impacts the evolution of adaptive genomic regions.

TE Hub: A community-oriented space for sharing and connecting tools, data, resources, and methods for transposable element annotation

Transposable elements (TEs) play powerful and varied evolutionary and functional roles, and are widespread in most eukaryotic genomes. Research into their unique biology has driven the creation of a large collection of databases, software, classification systems, and annotation guidelines. The diversity of available TE-related methods and resources raises compatibility concerns and can be overwhelming to researchers and communicators seeking straightforward guidance or materials. To address these challenges, we have initiated a new resource, TE Hub, that provides a space where members of the TE community can collaborate to document and create resources and methods. The space consists of (1) a website organized with an open wiki framework,  https://tehub.org , (2) a conversation framework via a Twitter account and a Slack channel, and (3) bi-monthly Hub Update video chats on the platform's development. In addition to serving as a centralized repository and communication platform, TE Hub lays the foundation for improved integration, standardization, and effectiveness of diverse tools and protocols. We invite the TE community, both novices and experts in TE identification and analysis, to join us in expanding our community-oriented resource.

In and out: Evolution of viral sequences in the mitochondrial genomes of legumes (Fabaceae)

Plant specific mitoviruses in the 'genus' Mitovirus (Narnaviridae) and their integrated sequences (non-retroviral endogenous RNA viral elements or NERVEs) have been recently identified in various plant lineages. However, the sparse phylogenetic coverage of complete plant mitochondrial genome (mitogenome) sequences and the non-conserved nature of mitochondrial intergenic regions have hindered comparative studies on mitovirus NERVEs in plants. In this study, 10 new mitogenomes were sequenced from legumes (Fabaceae). Based on comparative genomic analysis of 27 total mitogenomes, we identified mitovirus NERVEs and transposable elements across the family. All legume mitogenomes included NERVEs and total NERVE length varied from ca. 2 kb in the papilionoid Trifolium to 35 kb in the mimosoid Acacia. Most of the NERVE integration sites were in highly variable intergenic regions, however, some were positioned in six cis-spliced mitochondrial introns. In the Acacia mitogenome, there were L1-like transposon sequences including an almost full-length copy with target site duplications (TSDs). The integration sites of NERVEs in four introns showed evidence of L1-like retrotransposition events. Phylogenetic analysis revealed that there were multiple instances of precise deletion of NERVEs between TSDs. This study provides clear evidence that a L1-like retrotransposition mechanism has a long history of contributing to the integration of viral RNA into plant mitogenomes while microhomology-mediated deletion can restore the integration site.

Curation Guidelines for de novo Generated Transposable Element Families

Transposable elements (TEs) have the ability to alter individual genomic landscapes and shape the course of evolution for species in which they reside. Such profound changes can be understood by studying the biology of the organism and the interplay of the TEs it hosts. Characterizing and curating TEs across a wide range of species is a fundamental first step in this endeavor. This protocol employs techniques honed while developing TE libraries for a wide range of organisms and specifically addresses: (1) the extension of truncated de novo results into full-length TE families; (2) the iterative refinement of TE multiple sequence alignments; and (3) the use of alignment visualization to assess model completeness and subfamily structure. © 2021 Wiley Periodicals LLC. Basic Protocol: Extension and edge polishing of consensi and seed alignments derived from de novo repeat finders Support Protocol: Generating seed alignments using a library of consensi and a genome assembly.

Expression Quantitative Trait Loci (eQTLs) Associated with Retrotransposons Demonstrate their Modulatory Effect on the Transcriptome

Transposable elements (TEs) are repetitive elements that belong to a variety of functional classes and have an important role in shaping genome evolution. Around 50% of the human genome contains TEs, and they have been termed the "dark matter" of the genome because relatively little is known about their function. While TEs have been shown to participate in aberrant gene regulation and the pathogenesis of diseases, only a few studies have explored the systemic effect of TEs on gene expression. In the present study, we analysed whole genome sequences and blood whole transcriptome data from 570 individuals within the Parkinson's Progressive Markers Initiative (PPMI) cohort to identify expression quantitative trait loci (eQTL) regulating genome-wide gene expression associated with TEs. We identified 2132 reference TEs that were polymorphic for their presence or absence in our study cohort. The presence or absence of the TE element could change the expression of the gene or gene clusters from zero to tens of thousands of copies of RNA. The main finding is that many TEs possess very strong regulatory effects, and they have the potential to modulate large genetic networks with hundreds of target genes over the genome. We illustrate the plethora of regulatory mechanisms using examples of their action at the HLA gene cluster and data showing different TEs' convergence to modulate WFS1 gene expression. In conclusion, the presence or absence of polymorphisms of TEs has an eminent genome-wide regulatory function with large effect size at the level of the whole transcriptome. The role of TEs in explaining, in part, the missing heritability for complex traits is convincing and should be considered.

Two repeated motifs enriched within some enhancers and origins of replication are bound by SETMAR isoforms in human colon cells

Setmar is a gene specific to simian genomes. The function(s) of its isoforms are poorly understood and their existence in healthy tissues remains to be validated. Here we profiled SETMAR expression and its genome-wide binding landscape in colon tissue. We found isoforms V3 and V6 in healthy and tumour colon tissues as well as incell lines. In two colorectal cell lines SETMAR binds to several thousand Hsmar1 and MADE1 terminal ends, transposons mostly located in non-genic regions of active chromatin including in enhancers. It also binds to a 12-bp motifs similar to an inner motif in Hsmar1 and MADE1 terminal ends. This motif is interspersed throughout the genome and is enriched in GC-rich regions as well as in CpG islands that contain constitutive replication origins. It is also found in enhancers other than those associated with Hsmar1 and MADE1. The role of SETMAR in the expression of genes, DNA replication and in DNA repair are discussed.

Retrotransposable Elements: DNA Fingerprinting and the Assessment of Genetic Diversity

Retrotransposable elements (RTEs) are highly common mobile genetic elements that are composed of several classes and make up the majority of eukaryotic genomes. The "copy-out and paste-in" life cycle of replicative transposition in these dispersive and ubiquitous RTEs leads to new genome insertions without excision of the original element. RTEs are important drivers of species diversity; they exhibit great variety in structure, size, and mechanisms of transposition, making them important putative components in genome evolution. Accordingly, various applications have been developed to explore the polymorphisms in RTE insertion patterns. These applications include conventional or anchored polymerase chain reaction (PCR) and quantitative or digital PCR with primers designed for the 5' or 3' junction. Marker systems exploiting these PCR methods can be easily developed and are inexpensively used in the absence of extensive genome sequence data. The main inter-repeat amplification polymorphism techniques include inter-retrotransposon amplified polymorphism (IRAP), retrotransposon microsatellite amplified polymorphism (REMAP), and Inter-Primer Binding Site (iPBS) for PCR amplification with a single or two primers.

Opinion: Genetic Conflict With Mobile Elements Drives Eukaryotic Genome Evolution, and Perhaps Also Eukaryogenesis

Through analyses of diverse microeukaryotes, we have previously argued that eukaryotic genomes are dynamic systems that rely on epigenetic mechanisms to distinguish germline (i.e., DNA to be inherited) from soma (i.e., DNA that undergoes polyploidization, genome rearrangement, etc.), even in the context of a single nucleus. Here, we extend these arguments by including two well-documented observations: (1) eukaryotic genomes interact frequently with mobile genetic elements (MGEs) like viruses and transposable elements (TEs), creating genetic conflict, and (2) epigenetic mechanisms regulate MGEs. Synthesis of these ideas leads to the hypothesis that genetic conflict with MGEs contributed to the evolution of a dynamic eukaryotic genome in the last eukaryotic common ancestor (LECA), and may have contributed to eukaryogenesis (i.e., may have been a driver in the evolution of FECA, the first eukaryotic common ancestor). Sex (i.e., meiosis) may have evolved within the context of the development of germline-soma distinctions in LECA, as this process resets the germline genome by regulating/eliminating somatic (i.e., polyploid, rearranged) genetic material. Our synthesis of these ideas expands on hypotheses of the origin of eukaryotes by integrating the roles of MGEs and epigenetics.

A genomic survey of Tc1-mariner transposons in nematodes suggests extensive horizontal transposon transfer events

The number of reports concerning horizontal transposon transfers (HTT) in metazoan species is considerably increased, alongside with the exponential growth of genomic sequence data However, our understanding of the mechanisms of such phenomenon is still at an early stage. Nematodes constitute an animal phylum successfully adapted to almost every ecosystem and for this reason could potentially contribute to spreading the genetic information through horizontal transfer. To date, few studies describe HTT of nematode retrotransposons. This is due to the lack of annotation of transposable elements in the sequenced nematode genomes, especially DNA transposons, which are acknowledged as the best horizontal travelers among mobile sequences. We have therefore started a survey of DNA transposons and their possible involvement in HTT in sequenced nematode genomes. Here, we describe 83 new Tc1/mariner elements distributed in 17 nematode species. Among them, nine families were possibly horizontally transferred between nematodes and the most diverse animal species, including ants as preferred partner of HTT. The results obtained suggest that HTT events involving nematodes Tc1/mariner elements are not uncommon, and that nematodes could have a possible role as transposon reservoir that, in turn, can be redistributed among animal genomes. Overall, this could be relevant to understand how the inter-species genetic flows shape the landscape of genetic variation of organisms inhabiting specific environmental communities.

The piggyBac-derived protein 5 (PGBD5) transposes both the closely and the distantly related piggyBac-like elements Tcr-pble and Ifp2

The vertebrate piggyBac derived transposase 5 (PGBD5) encodes a domesticated transposase, which is active and able to transpose its distantly related piggyBac-like element (pble), Ifp2. This raised the question whether PGBD5 would be more effective at mobilizing a phylogenetically closely related pble element. We aimed to identify the pble most closely related to the pgbd5 gene. We updated the landscape of vertebrate pgbd genes to develop efficient filters and identify the most closely related pble to each of these genes. We found that Tcr-pble is phylogenetically the closest pble to the pgbd5 gene. Furthermore, we evaluated the capacity of two murine and human PGBD5 isoforms, Mm523 and Hs524, to transpose both Tcr-pble and Ifp2 elements. We found that both pbles could be transposed by Mm523 with similar efficiency. However, integrations of both pbles occurred through both proper transposition and improper PGBD5-dependent recombination. This suggested that the ability of PGBD5 to bind both pbles may not be based on the primary sequence of element ends, but may involve recognition of inner DNA motifs, possibly related to palindromic repeats. In agreement with this hypothesis, we identified internal palindromic repeats near the end of 24 pble sequences, which display distinct sequences.

The C-terminal Domain of piggyBac Transposase Is Not Required for DNA Transposition

PiggyBac(PB)-like elements (pble) are members of a eukaryotic DNA transposon family. This family is of interest to evolutionary genomics because pble transposases have been domesticated at least 9 times in vertebrates. The amino acid sequence of pble transposases can be split into three regions: an acidic N-terminal domain (~100 aa), a central domain (~400 aa) containing a DD[D/E] catalytic triad, and a cysteine-rich domain (CRD; ~90 aa). Two recent reports suggested that a functional CRD is required for pble transposase activity. Here we found that two CRD-deficient pble transposases, a PB variant and an isoform encoded by the domesticated PB-derived vertebrate transposase gene 5 (pgbd5) trigger transposition of the Ifp2 pble. When overexpressed in HeLa cells, these CRD-deficient transposases can insert Ifp2 elements with proper and improper transposon ends, associated with deleterious effects on cells. Finally, we found that mouse CRD-deficient transposase Pgbd5, as well as PB, do not insert pbles at random into chromosomes. Transposition events occurred more often in genic regions, in the neighbourhood of the transcription start sites and were often found in genes predominantly expressed in the human central nervous system.

The Dfam community resource of transposable element families, sequence models, and genome annotations

Dfam is an open access database of repetitive DNA families, sequence models, and genome annotations. The 3.0-3.3 releases of Dfam ( https://dfam.org ) represent an evolution from a proof-of-principle collection of transposable element families in model organisms into a community resource for a broad range of species, and for both curated and uncurated datasets. In addition, releases since Dfam 3.0 provide auxiliary consensus sequence models, transposable element protein alignments, and a formalized classification system to support the growing diversity of organisms represented in the resource. The latest release includes 266,740 new de novo generated transposable element families from 336 species contributed by the EBI. This expansion demonstrates the utility of many of Dfam's new features and provides insight into the long term challenges ahead for improving de novo generated transposable element datasets.

Genome wide annotation and characterization of young, intact long terminal repeat retrotransposons (In-LTR-RTs) of seven legume species

Availability of genome sequence of different legume species has provided an opportunity to characterize the abundance, distribution, and divergence of canonical intact long terminal retrotransposons (In-LTR-RT) superfamilies. Among seven legume species, Arachis ipaensis (Aip) showed the highest number of full-length canonical In-LTR-RTs (3325), followed by Glycine max (Gma, 2328), Vigna angularis (Van, 1625), Arachis durensis (Adu, 1348), Lotus japonicus (Lja, 1294), Medicago truncatula (Mtr, 788), and Circer arietinum (Car, 124). Divergence time analysis demonstrated that the amplification timeframe of LTR-RTs dramatically varied in different families. The average insertion time of Copia element varied from 0.51 (Van) to 1.37 million years ago (Mya) (Adu, and Aip), whereas that of Gypsy was between 0.22 (Mtr) and 1.82 Mya (Adu). Bayesian phylogenetic tree analysis suggested that the 1397 and 1917 reverse transcriptase (RT) domains of Copia and Gypsy families of the seven legume species were clustered into 7 and 14 major groups, respectively. The highest proportion (approximately 94.79-100%) of transposable element (TE)-associated genes assigned to pathways was mapped to metabolism-related pathways in all species. The results enabled the structural understanding of full-length In-LTR-RTs and will be valuable resource for the further study of the impact of TEs on gene structure and expression in legume species.

Measuring and interpreting transposable element expression

Transposable elements (TEs) are insertional mutagens that contribute greatly to the plasticity of eukaryotic genomes, influencing the evolution and adaptation of species as well as physiology or disease in individuals. Measuring TE expression helps to understand not only when and where TE mobilization can occur but also how this process alters gene expression, chromatin accessibility or cellular signalling pathways. Although genome-wide gene expression assays such as RNA sequencing include transposon-derived transcripts, most computational analytical tools discard or misinterpret TE-derived reads. Emerging approaches are improving the identification of expressed TE loci and helping to discriminate TE transcripts that permit TE mobilization from chimeric gene-TE transcripts or pervasive transcription. Here we review the main challenges associated with the detection of TE expression, including mappability, insertional and internal sequence polymorphisms, and the diversity of the TE transcriptional landscape, as well as the different experimental and computational strategies to solve them.

RepeatModeler2 for automated genomic discovery of transposable element families

The accelerating pace of genome sequencing throughout the tree of life is driving the need for improved unsupervised annotation of genome components such as transposable elements (TEs). Because the types and sequences of TEs are highly variable across species, automated TE discovery and annotation are challenging and time-consuming tasks. A critical first step is the de novo identification and accurate compilation of sequence models representing all of the unique TE families dispersed in the genome. Here we introduce RepeatModeler2, a pipeline that greatly facilitates this process. This program brings substantial improvements over the original version of RepeatModeler, one of the most widely used tools for TE discovery. In particular, this version incorporates a module for structural discovery of complete long terminal repeat (LTR) retroelements, which are widespread in eukaryotic genomes but recalcitrant to automated identification because of their size and sequence complexity. We benchmarked RepeatModeler2 on three model species with diverse TE landscapes and high-quality, manually curated TE libraries: Drosophila melanogaster (fruit fly), Danio rerio (zebrafish), and Oryza sativa (rice). In these three species, RepeatModeler2 identified approximately 3 times more consensus sequences matching with >95% sequence identity and sequence coverage to the manually curated sequences than the original RepeatModeler. As expected, the greatest improvement is for LTR retroelements. Thus, RepeatModeler2 represents a valuable addition to the genome annotation toolkit that will enhance the identification and study of TEs in eukaryotic genome sequences. RepeatModeler2 is available as source code or a containerized package under an open license (https://github.com/Dfam-consortium/RepeatModeler, http://www.repeatmasker.org/RepeatModeler/).

Repetitive DNA landscape in essential A and supernumerary B chromosomes of Festuca pratensis Huds

Here, we characterized the basic properties of repetitive sequences in essential A and supernumerary B chromosomes of Festuca pratensis Huds. This was performed by comparative analysis of low-pass Illumina sequence reads of B chromosome lacking (-B) and B chromosome containing (+B) individuals of F. pratensis. 61% of the nuclear genome is composed of repetitive sequences. 43.1% of the genome are transposons of which DNA transposons and retrotransposons made up 2.3% and 40.8%, respectively. LTR retrotransposons are the most abundant mobile elements and contribute to 40.7% of the genome and divided into Ty3-gypsy and Ty1-copia super families with 32.97% and 7.78% of the genome, respectively. Eighteen different satellite repeats were identified making up 3.9% of the genome. Five satellite repeats were used as cytological markers for chromosome identification and genome analysis in the genus Festuca. Four satellite repeats were identified on B chromosomes among which Fp-Sat48 and Fp-Sat253 were specific to the B chromosome of F. pratensis.

A systematic review of the application of machine learning in the detection and classification of transposable elements

Background

Transposable elements (TEs) constitute the most common repeated sequences in eukaryotic genomes. Recent studies demonstrated their deep impact on species diversity, adaptation to the environment and diseases. Although there are many conventional bioinformatics algorithms for detecting and classifying TEs, none have achieved reliable results on different types of TEs. Machine learning (ML) techniques can automatically extract hidden patterns and novel information from labeled or non-labeled data and have been applied to solving several scientific problems.

Methodology

We followed the Systematic Literature Review (SLR) process, applying the six stages of the review protocol from it, but added a previous stage, which aims to detect the need for a review. Then search equations were formulated and executed in several literature databases. Relevant publications were scanned and used to extract evidence to answer research questions.

Results

Several ML approaches have already been tested on other bioinformatics problems with promising results, yet there are few algorithms and architectures available in literature focused specifically on TEs, despite representing the majority of the nuclear DNA of many organisms. Only 35 articles were found and categorized as relevant in TE or related fields.

Conclusions

ML is a powerful tool that can be used to address many problems. Although ML techniques have been used widely in other biological tasks, their utilization in TE analyses is still limited. Following the SLR, it was possible to notice that the use of ML for TE analyses (detection and classification) is an open problem, and this new field of research is growing in interest.

Stowaway miniature inverted repeat transposable elements are important agents driving recent genomic diversity in wild and cultivated carrot

BACKGROUND

Miniature inverted repeat transposable elements (MITEs) are small non-autonomous DNA transposons that are ubiquitous in plant genomes, and are mobilised by their autonomous relatives. Stowaway MITEs are derived from and mobilised by elements from the mariner superfamily. Those elements constitute a significant portion of the carrot genome; however the variation caused by Daucus carota Stowaway MITEs (DcStos), their association with genes and their putative impact on genome evolution has not been comprehensively analysed.

RESULTS

Fourteen families of Stowaway elements DcStos occupy about 0.5% of the carrot genome. We systematically analysed 31 genomes of wild and cultivated Daucus carota, yielding 18.5 thousand copies of these elements, showing remarkable insertion site polymorphism. DcSto element demography differed based on the origin of the host populations, and corresponded with the four major groups of D. carota, wild European, wild Asian, eastern cultivated and western cultivated. The DcStos elements were associated with genes, and most frequently occurred in 5' and 3' untranslated regions (UTRs). Individual families differed in their propensity to reside in particular segments of genes. Most importantly, DcSto copies in the 2 kb regions up- and downstream of genes were more frequently associated with open reading frames encoding transcription factors, suggesting their possible functional impact. More than 1.5% of all DcSto insertion sites in different host genomes contained different copies in exactly the same position, indicating the existence of insertional hotspots. The DcSto7b family was much more polymorphic than the other families in cultivated carrot. A line of evidence pointed at its activity in the course of carrot domestication, and identified Dcmar1 as an active carrot mariner element and a possible source of the transposition machinery for DcSto7b.

CONCLUSION

Stowaway MITEs have made a substantial contribution to the structural and functional variability of the carrot genome.

Evolutionary history and classification of Micropia retroelements in Drosophilidae species

Transposable elements (TEs) have the main role in shaping the evolution of genomes and host species, contributing to the creation of new genes and promoting rearrangements frequently associated with new regulatory networks. Support for these hypotheses frequently results from studies with model species, and Drosophila provides a great model organism to the study of TEs. Micropia belongs to the Ty3/Gypsy group of long terminal repeats (LTR) retroelements and comprises one of the least studied Drosophila transposable elements. In this study, we assessed the evolutionary history of Micropia within Drosophilidae, while trying to assist in the classification of this TE. At first, we performed searches of Micropia presence in the genome of natural populations from several species. Then, based on searches within online genomic databases, we retrieved Micropia-like sequences from the genomes of distinct Drosophilidae species. We expanded the knowledge of Micropia distribution within Drosophila species. The Micropia retroelements we detected consist of an array of divergent sequences, which we subdivided into 20 subfamilies. Even so, a patchy distribution of Micropia sequences within the Drosophilidae phylogeny could be identified, with incongruences between the species phylogeny and the Micropia phylogeny. Comparing the pairwise synonymous distance (dS) values between Micropia and three host nuclear sequences, we found several cases of unexpectedly high levels of similarity between Micropia sequences in divergent species. All these findings provide a hypothesis to the evolution of Micropia within Drosophilidae, which include several events of vertical and horizontal transposon transmission, associated with ancestral polymorphisms and recurrent Micropia sequences diversification.

Characterization of Herves-like transposable elements (hATs) in Drosophila species and their evolutionary scenario

A monophyletic group of Drosophila hAT transposable elements, referred to as Herves-like, was characterized and found to be present in 46% of 57 screened Drosophila species. A remarkable characteristic of these elements is the presence of a long array of minisatellite repeats (MnRs) in both subterminal extremities of the elements. The copy number of these minisatellites was highly variable between and within populations. Twenty-three strains of Drosophila willistoni, covering its geographic distribution, were screened for polymorphism in the copy number of 5' MnRs, showing a variation from 7 to 20 repeat copies. These MnRs are well conserved among Drosophila species and probably function as transposase binding sequences, as provided by short subterminal repeats in other hAT elements. Miniature inverted repeat transposable elements were found in 27% of species carrying Herves-like elements. Phylogenetic analysis showed incongruences between transposable elements and species phylogenies, suggesting that at least four horizontal transfer events have occurred.

Retrotransposons in Plant Genomes: Structure, Identification, and Classification through Bioinformatics and Machine Learning

Transposable elements (TEs) are genomic units able to move within the genome of virtually all organisms. Due to their natural repetitive numbers and their high structural diversity, the identification and classification of TEs remain a challenge in sequenced genomes. Although TEs were initially regarded as "junk DNA", it has been demonstrated that they play key roles in chromosome structures, gene expression, and regulation, as well as adaptation and evolution. A highly reliable annotation of these elements is, therefore, crucial to better understand genome functions and their evolution. To date, much bioinformatics software has been developed to address TE detection and classification processes, but many problematic aspects remain, such as the reliability, precision, and speed of the analyses. Machine learning and deep learning are algorithms that can make automatic predictions and decisions in a wide variety of scientific applications. They have been tested in bioinformatics and, more specifically for TEs, classification with encouraging results. In this review, we will discuss important aspects of TEs, such as their structure, importance in the evolution and architecture of the host, and their current classifications and nomenclatures. We will also address current methods and their limitations in identifying and classifying TEs.

Transposable Elements Adaptive Role in Genome Plasticity, Pathogenicity and Evolution in Fungal Phytopathogens

Transposable elements (TEs) are agents of genetic variability in phytopathogens as they are a source of adaptive evolution through genome diversification. Although many studies have uncovered information on TEs, the exact mechanism behind TE-induced changes within the genome remains poorly understood. Furthermore, convergent trends towards bigger genomes, emergence of novel genes and gain or loss of genes implicate a TE-regulated genome plasticity of fungal phytopathogens. TEs are able to alter gene expression by revamping the cis-regulatory elements or recruiting epigenetic control. Recent findings show that TEs recruit epigenetic control on the expression of effector genes as part of the coordinated infection strategy. In addition to genome plasticity and diversity, fungal pathogenicity is an area of economic concern. A survey of TE distribution suggests that their proximity to pathogenicity genes TEs may act as sites for emergence of novel pathogenicity factors via nucleotide changes and expansion or reduction of the gene family. Through a systematic survey of literature, we were able to conclude that the role of TEs in fungi is wide: ranging from genome plasticity, pathogenicity to adaptive behavior in evolution. This review also identifies the gaps in knowledge that requires further elucidation for a better understanding of TEs' contribution to genome architecture and versatility.