Transposable Elements In Genomes Research Articles

Transposable element diversity and activity patterns in neotropical salamanders.

Transposable elements (TEs) compose a substantial proportion of the largest eukaryotic genomes. TE diversity has been hypothesized to be negatively correlated with genome size, yet empirical demonstrations of such a relationship in a phylogenetic context are largely lacking. Furthermore, the most abundant type of TEs in genomes varies across groups, and it is not clear if there are patterns of TE activity consistent with genome size among different taxa with large genome sizes. We use low-coverage sequencing of 16 species of Neotropical salamanders, which vary approximately seven-fold in genome size, to estimate TE relative abundance and diversity for each species. We also compare divergence of copies of each TE superfamily to estimate patterns of TE activity in each species. We find a negative relationship between TE diversity and genome size, which is consistent with the hypothesis that either competition among TEs or reduced selection against ectopic recombination may result in lower diversity in the largest genomes. We also find divergent activity patterns in the largest vs. the smallest genomes, suggesting that the history of TE activity may explain differences in genome size. Our results suggest that both TE diversity and relative abundance may be predictable, at least within taxonomic groups.

Prevalent fast evolution of genes involved in heterochromatin functions.

Heterochromatin is a gene-poor and repeat-rich genomic compartment universally found in eukaryotes. Despite its low transcriptional activity, heterochromatin plays important roles in maintaining genome stability, organizing chromosomes, and suppressing transposable elements (TEs). Given the importance of these functions, it is expected that the genes involved in heterochromatin regulation would be highly conserved. Yet, a handful of these genes were found to evolve rapidly. To investigate whether these previous findings are anecdotal or general to genes modulating heterochromatin, we compile an exhaustive list of 106 candidate genes involved in heterochromatin functions and investigate their evolution over short and long evolutionary time scales in Drosophila. Our analyses find that these genes exhibit significantly more frequent evolutionary changes, both in the forms of amino acid substitutions and gene copy number change, when compared to genes involved in Polycomb-based repressive chromatin. While positive selection drives amino acid changes within both structured domains with diverse functions and intrinsically disordered regions (IDRs), purifying selection may have maintained the proportions of IDRs of these proteins. Together with the observed negative associations between evolutionary rates of these genes and genomic TE abundance, we propose an evolutionary model where the fast evolution of genes involved in heterochromatin functions is an inevitable outcome of the unique functional roles of heterochromatin, while the rapid evolution of TEs may be an effect rather than cause. Our study provides an important global view of the evolution of genes involved in this critical cellular domain and provides insights into the factors driving the distinctive evolution of heterochromatin.

Technology to the rescue: how to uncover the role of transposable elements in preimplantation development.

Transposable elements (TEs) are highly expressed in preimplantation development. Preimplantation development is the phase when the cells of the early embryo undergo the first cell fate choice and change from being totipotent to pluripotent. A range of studies have advanced our understanding of TEs in preimplantation, as well as their epigenetic regulation and functional roles. However, many questions remain about the implications of TE expression during early development. Challenges originate first due to the abundance of TEs in the genome, and second because of the limited cell numbers in preimplantation. Here we review the most recent technological advancements promising to shed light onto the role of TEs in preimplantation development. We explore novel avenues to identify genomic TE insertions and improve our understanding of the regulatory mechanisms and roles of TEs and their RNA and protein products during early development.

Endogenous Retroviruses Unveiled: A Comprehensive Review of Inflammatory Signaling/Senescence-Related Pathways and Therapeutic Strategies.

Endogenous retroviruses (ERVs), a subset of genomic transposable elements (TEs) in a broader sense, have remained latent within mammalian genomes for tens of millions of years. These genetic elements are typically in a silenced state due to stringent regulatory mechanisms. However, under specific conditions, they can become activated, triggering inflammatory responses through diverse mechanisms. This activation has been shown to play a potential role in various neurological disorders, tumors, and cellular senescence. Consequently, the regulation of ERV expression through various methods holds promise for clinical applications in disease treatment. ERVs also engage in interactions with a variety of exogenous viruses, thereby influencing the outcomes of viral infectious diseases. This article comprehensively reviews the pathogenic cascade of ERVs, encompassing activation, inflammation, associated diseases, senescence, and interplay with viruses. Additionally, it outlines therapeutic strategies targeting ERVs with the aim of offering novel research directions for understanding the relationship between ERVs and diseases, along with corresponding treatment modalities.

Insights into the Evolution of Ohnologous Sequences and Their Epigenetic Marks Post-WGD in Malus Domestica.

A Whole Genome Duplication (WGD) event occurred several Ma in a Rosaceae ancestor, giving rise to the Maloideae subfamily which includes today many pome fruits such as pear (Pyrus communis) and apple (Malus domestica). This complete and well-conserved genome duplication makes the apple an organism of choice to study the early evolutionary events occurring to ohnologous chromosome fragments. In this study, we investigated gene sequence evolution and expression, transposable elements (TE) density, and DNA methylation level. Overall, we identified 16,779 ohnologous gene pairs in the apple genome, confirming the relatively recent WGD. We identified several imbalances in QTL localization among duplicated chromosomal fragments and characterized various biases in genome fractionation, gene transcription, TE densities, and DNA methylation. Our results suggest a particular chromosome dominance in this autopolyploid species, a phenomenon that displays similarities with subgenome dominance that has only been described so far in allopolyploids.

Perspective for Studying the Relationship of miRNAs with Transposable Elements.

Transposable elements are important sources of miRNA, long non-coding RNAs genes, and their targets in the composition of protein-coding genes in plants and animals. Therefore, the detection of expression levels of specific non-coding RNAs in various tissues and cells in normal and pathological conditions may indicate a programmed pattern of transposable elements' activation. This reflects the species-specific composition and distribution of transposable elements in genomes, which underlie gene regulation in every cell division, including during aging. TEs' expression is also regulated by epigenetic factors (DNA methylation, histone modifications), SIRT6, cytidine deaminases APOBEC3, APOBEC1, and other catalytic proteins, such as ERCC, TREX1, RB1, HELLS, and MEGP2. In evolution, protein-coding genes and their regulatory elements are derived from transposons. As part of non-coding regions and introns of genes, they are sensors for transcriptional and post-transcriptional control of expression, using miRNAs and long non-coding RNAs, that arose from transposable elements in evolution. Methods (Orbld, ncRNAclassifier) and databases have been created for determining the occurrence of miRNAs from transposable elements in plants (PlanTE-MIR DB, PlaNC-TE), which can be used to design epigenetic gene networks in ontogenesis. Based on the data accumulated in the scientific literature, the presence of 467 transposon-derived miRNA genes in the human genome has been reliably established. It was proposed to create an updated and controlled online bioinformatics database of miRNAs derived from transposable elements in healthy individuals, as well as expression changes of these miRNAs during aging and various diseases, such as cancer and difficult-to-treat diseases. The use of the information obtained can open new horizons in the management of tissue and organ differentiation to aging slow down. In addition, the created database could become the basis for clarifying the mechanisms of pathogenesis of various diseases (imbalance in the activity of transposable elements, reflected in changes in the expression of miRNAs) and designing their targeted therapy using specific miRNAs as targets. This article provides examples of the detection of transposable elements-derived miRNAs involved in the development of specific malignant neoplasms, aging, and idiopathic pulmonary fibrosis.

A population genetics theory for piRNA-regulated transposable elements

Transposable elements (TEs) are self-reproducing selfish DNA sequences that can invade the genome of virtually all living species. Population genetics models have shown that TE copy numbers generally reach a limit, either because the transposition rate decreases with the number of copies (transposition regulation) or because TE copies are deleterious, and thus purged by natural selection. Yet, recent empirical discoveries suggest that TE regulation may mostly rely on piRNAs, which require a specific mutational event (the insertion of a TE copy in a piRNA cluster) to be activated — the so-called TE regulation “trap model”. We derived new population genetics models accounting for this trap mechanism, and showed that the resulting equilibria differ substantially from previous expectations based on a transposition–selection equilibrium. We proposed three sub-models, depending on whether or not genomic TE copies and piRNA cluster TE copies are selectively neutral or deleterious, and we provide analytical expressions for maximum and equilibrium copy numbers, as well as cluster frequencies for all of them. In the full neutral model, the equilibrium is achieved when transposition is completely silenced, and this equilibrium does not depend on the transposition rate. When genomic TE copies are deleterious but not cluster TE copies, no long-term equilibrium is possible, and active TEs are eventually eliminated after an active incomplete invasion stage. When all TE copies are deleterious, a transposition–selection equilibrium exists, but the invasion dynamics is not monotonic, and the copy number peaks before decreasing. Mathematical predictions were in good agreement with numerical simulations, except when genetic drift and/or linkage disequilibrium dominates. Overall, the trap-model dynamics appeared to be substantially more stochastic and less repeatable than traditional regulation models.

Transposable element and host silencing activity in gigantic genomes.

Transposable elements (TEs) and the silencing machinery of their hosts are engaged in a germline arms-race dynamic that shapes TE accumulation and, therefore, genome size. In animal species with extremely large genomes (>10Gb), TE accumulation has been pushed to the extreme, prompting the question of whether TE silencing also deviates from typical conditions. To address this question, we characterize TE silencing via two pathways-the piRNA pathway and KRAB-ZFP transcriptional repression-in the male and female gonads of Ranodon sibiricus, a salamander species with a ∼21Gb genome. We quantify 1) genomic TE diversity, 2) TE expression, and 3) small RNA expression and find a significant relationship between the expression of piRNAs and TEs they target for silencing in both ovaries and testes. We also quantified TE silencing pathway gene expression in R. sibiricus and 14 other vertebrates with genome sizes ranging from 1 to 130Gb and find no association between pathway expression and genome size. Taken together, our results reveal that the gigantic R. sibiricus genome includes at least 19 putatively active TE superfamilies, all of which are targeted by the piRNA pathway in proportion to their expression levels, suggesting comprehensive piRNA-mediated silencing. Testes have higher TE expression than ovaries, suggesting that they may contribute more to the species' high genomic TE load. We posit that apparently conflicting interpretations of TE silencing and genomic gigantism in the literature, as well as the absence of a correlation between TE silencing pathway gene expression and genome size, can be reconciled by considering whether the TE community or the host is currently "on the attack" in the arms race dynamic.

PlanTEnrichment: A How-to Guide on Rapid Identification of Transposable Elements Associated with Regions of Interest in Select Plant Genomes.

Transposable elements (TEs) are repeat elements that can relocate or create novel copies of themselves in the genome and contribute to genomic complexity and expansion, via events such as chromosome recombination or regulation of gene expression. However, given the large number of such repeats across the genome, identifying repeats of interest can be a challenge in even well-annotated genomes, especially in more complex, TE-rich plant genomes. Here, we describe a protocol for PlanTEnrichment, a database we created comprising information on 11 plant genomes to analyze stress-associated TEs using publicly available data. By selecting a genome and providing a list of genes or genomic regions whose TE associations the user wants to identify, the user can rapidly obtain TE subfamilies found near the provided regions, as well as their superfamily and class, and the enrichment values of the repeats. The results also provide the locations of individual repeat instances found, alongside the input regions or genes they are associated with, and a bar graph of the top ten most significant repeat subfamilies identified. PlanTEnrichment is freely available at http://tools.ibg.deu.edu.tr/plantenrichment/ and can be used by researchers with rudimentary or no proficiency in computational analysis of TE elements, allowing for expedience in the identification of TEs of interest and helping further our understanding of the potential contributions of TEs in plant genomes.

Genomic Diversity and Chromosomal Rearrangements in Neisseria gonorrhoeae and Neisseria meningitidis

Chromosomal rearrangements in N. gonorrhoeae and N. meningitidis were studied with the determination of mobile elements and their role in rearrangements. The results of whole-genome sequencing and de novo genome assembly for 50 N. gonorrhoeae isolates collected in Russia were compared with 96 genomes of N. gonorrhoeae and 138 genomes of N. meningitidis from the databases. Rearrangement events with the determination of the coordinates of syntenic blocks were analyzed using the SibeliaZ software v.1.2.5, the minimum number of events that allow one genome to pass into another was calculated using the DCJ-indel model using the UniMoG program v.1.0. Population-level analysis revealed a stronger correlation between changes in the gene order and phylogenetic proximity for N. meningitidis in contrast to N. gonorrhoeae. Mobile elements were identified, including Correa elements; Spencer-Smith elements (in N. gonorrhoeae); Neisserial intergenic mosaic elements; IS elements of IS5, IS30, IS110, IS1595 groups; Nf1-Nf3 prophages; NgoФ1-NgoФ9 prophages; and Mu-like prophages Pnm1, Pnm2, MuMenB (in N. meningitidis). More than 44% of the observed rearrangements most likely occurred with the participation of mobile elements, including prophages. No differences were found between the Russian and global N. gonorrhoeae population both in terms of rearrangement events and in the number of transposable elements in genomes.

Non-reference genome transposable elements (TEs) have a significant impact on the progression of the Parkinson's disease.

The pathophysiology of Parkinson's disease (PD) is a complex process of the interaction between genetic and environmental factors. Studies on the genetic component of PD have predominantly focused on single nucleotide polymorphisms (SNPs) using a cross-sectional case-control design in large genome-wide association studies. This approach while giving insight into a significant portion of the genetics of PD does not fully account for all the genetic components resulting in missing heritability. In this study, we approached this problem by focusing on the non-reference genome transposable elements (TEs) and their impact on the progression of PD using a longitudinal study design within the Parkinson's progression markers initiative (PPMI) cohort. We analyzed 2886 Alu repeats, 360 LINE1 and 128 SINE-VNTR-Alus (SVAs) that were called from the whole-genome sequence data which are not within the reference genome. The presence or absence of these non-reference TE variants is known as a retrotransposon insertion polymorphism, and measuring this polymorphism describes the impact of TEs on the traits. The variations for the presence or absence of the non-reference TE elements were modeled to align with the changes in the 114 outcome measures during the five-year follow-up period of the PPMI cohort. Linear mixed-effects models were used, and many TEs were found to have a highly significant effect on the longitudinal changes in the clinically important PD outcomes such as UPDRS subscale II, UPDRS total scores, and modified Schwab and England ADL scale. In addition, the progression of several imaging and functional measures, including the Caudate/Putamen ratio and levodopa equivalent daily dose (LEDD) were also significantly affected by the TEs. In conclusion, this study identified the overwhelming effect of the non-reference TEs on the progression of PD and is a good example of the impact the variations in the "junk DNA" have on complex diseases.

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

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

Repeat Dynamics across Timescales: A Perspective from Sibling Allotetraploid Marsh Orchids (Dactylorhiza majalis s.l.).

To provide insights into the fate of transposable elements (TEs) across timescales in a post-polyploidization context, we comparatively investigate five sibling Dactylorhiza allotetraploids (Orchidaceae) formed independently and sequentially between 500 and 100K generations ago by unidirectional hybridization between diploids D. fuchsii and D. incarnata. Our results first reveal that the paternal D. incarnata genome shows a marked increased content of LTR retrotransposons compared to the maternal species, reflected in its larger genome size and consistent with a previously hypothesized bottleneck. With regard to the allopolyploids, in the youngest D. purpurella both genome size and TE composition appear to be largely additive with respect to parents, whereas for polyploids of intermediate ages we uncover rampant genome expansion on a magnitude of multiple entire genomes of some plants such as Arabidopsis. The oldest allopolyploids in the series are not larger than the intermediate ones. A putative tandem repeat, potentially derived from a non-autonomous miniature inverted-repeat TE (MITE) drives much of the genome dynamics in the allopolyploids. The highly dynamic MITE-like element is found in higher proportions in the maternal diploid, D. fuchsii, but is observed to increase in copy number in both subgenomes of the allopolyploids. Altogether, the fate of repeats appears strongly regulated and therefore predictable across multiple independent allopolyploidization events in this system. Apart from the MITE-like element, we consistently document a mild genomic shock following the allopolyploidizations investigated here, which may be linked to their relatively large genome sizes, possibly associated with strong selection against further genome expansions.

Transposable element landscapes in aging Drosophila.

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

Taming, Domestication and Exaptation: Trajectories of Transposable Elements in Genomes.

During evolution, several types of sequences pass through genomes. Along with mutations and internal genetic tinkering, they are a useful source of genetic variability for adaptation and evolution. Most of these sequences are acquired by horizontal transfers (HT), but some of them may come from the genomes themselves. If they are not lost or eliminated quickly, they can be tamed, domesticated, or even exapted. Each of these processes results from a series of events, depending on the interactions between these sequences and the host genomes, but also on environmental constraints, through their impact on individuals or population fitness. After a brief reminder of the characteristics of each of these states (taming, domestication, exaptation), the evolutionary trajectories of these new or acquired sequences will be presented and discussed, emphasizing that they are not totally independent insofar as the first can constitute a step towards the second, and the second is another step towards the third.

Roles of Noncoding RNA in Reproduction.

The World Health Organization predicts that infertility will be the third major health threat after cancer and cardiovascular disease, and will become a hot topic in medical research. Studies have shown that epigenetic changes are an important component of gametogenesis and related reproductive diseases. Epigenetic regulation of noncoding RNA (ncRNA) is appropriate and is a research hotspot in the biomedical field; these include long noncoding RNA (lncRNA), microRNA (miRNA), and PIWI-interacting RNA (piRNA). As vital members of the intracellular gene regulatory network, they affect various life activities of cells. LncRNA functions as a molecular bait, molecular signal and molecular scaffold in the body through molecular guidance. miRNAs are critical regulators of gene expression; they mainly control the stability or translation of their target mRNA after transcription. piRNA functions mainly through silencing genomic transposable elements and the post-transcriptional regulation of mRNAs in animal germ cells. Current studies have shown that these ncRNAs also play significant roles in the reproductive system and are involved in the regulation of essential cellular events in spermatogenesis and follicular development. The abnormal expression of ncRNA is closely linked to testicular germ cell tumors, poly cystic ovary syndrome and other diseases. This paper briefly presents the research on the reproductive process and reproductive diseases involving ncRNAs.

An Atlas of Plant Transposable Elements.

Advances in genomic sequencing have recently offered vast opportunities for biological exploration, unraveling the evolution and improving our understanding of Earth biodiversity. Due to distinct plant species characteristics in terms of genome size, ploidy and heterozygosity, transposable elements (TEs) are common characteristics of many genomes. TEs are ubiquitous and dispersed repetitive DNA sequences that frequently impact the evolution and composition of the genome, mainly due to their redundancy and rearrangements. For this study, we provided an atlas of TE data by employing an easy-to-use portal ( APTE website ). To our knowledge, this is the most extensive and standardized analysis of TEs in plant genomes. We evaluated 67 plant genomes assembled at chromosome scale, recovering a total of 49,802,023 TE records, representing a total of 47,992,091,043 (~47,62%) base pairs (bp) of the total genomic space. We observed that new types of TEs were identified and annotated compared to other data repositories. By establishing a standardized catalog of TE annotation on 67 genomes, new hypotheses, exploration of TE data and their influences on the genomes may allow a better understanding of their function and processes. All original code and an example of how we developed the TE annotation strategy is available on GitHub ( Extended data).

Identifying piRNA targets on mRNAs in C. elegans using a deep multi-head attention network

BackgroundPiwi-interacting RNAs (piRNAs) are the small non-coding RNAs (ncRNAs) that silence genomic transposable elements. And researchers found out that piRNA also regulates various endogenous transcripts. However, there is no systematic understanding of the piRNA binding patterns and how piRNA targets genes. While various prediction methods have been developed for other similar ncRNAs (e.g., miRNAs), piRNA holds distinctive characteristics and requires its own computational model for binding target prediction.ResultsRecently, transcriptome-wide piRNA binding events in C. elegans were probed by PRG-1 CLASH experiments. Based on the probed piRNA-messenger RNAs (mRNAs) binding pairs, in this research, we devised the first deep learning architecture based on multi-head attention to computationally identify piRNA targeting mRNA sites. In the devised deep network, the given piRNA and mRNA segment sequences are first one-hot encoded and undergo a combined operation of convolution and squeezing-extraction to unravel motif patterns. And we incorporate a novel multi-head attention sub-network to extract the hidden piRNA binding rules that can simulate the biological piRNA target recognition process. Finally, the true piRNA–mRNA binding pairs are identified by a deep fully connected sub-network. Our model obtains a supreme discriminatory power of AUC = 93.3% on an independent test set and successfully extracts the verified binding pattern of a synthetic piRNA. These results demonstrated that the devised model achieves high prediction performance and suggests testable potential biological piRNA binding rules.ConclusionsIn this research, we developed the first deep learning method to identify piRNA targeting sites on C. elegans mRNAs. And the developed deep learning method is demonstrated to be of high accuracy and can provide biological insights into piRNA–mRNA binding patterns. The piRNA binding target identification network can be downloaded from http://cosbi2.ee.ncku.edu.tw/data_download/piRNA_mRNA_binding.

Monitoring Insect Transposable Elements in Large Double-Stranded DNA Viruses Reveals Host-to-Virus and Virus-to-Virus Transposition.

The mechanisms by which transposable elements (TEs) can be horizontally transferred between animals are unknown, but viruses are possible candidate vectors. Here, we surveyed the presence of host-derived TEs in viral genomes in 35 deep sequencing data sets produced from 11 host–virus systems, encompassing nine arthropod host species (five lepidopterans, two dipterans, and two crustaceans) and six different double-stranded (ds) DNA viruses (four baculoviruses and two iridoviruses). We found evidence of viral-borne TEs in 14 data sets, with frequencies of viral genomes carrying a TE ranging from 0.01% to 26.33% for baculoviruses and from 0.45% to 7.36% for iridoviruses. The analysis of viral populations separated by a single replication cycle revealed that viral-borne TEs originating from an initial host species can be retrieved after viral replication in another host species, sometimes at higher frequencies. Furthermore, we detected a strong increase in the number of integrations in a viral population for a TE absent from the hosts’ genomes, indicating that this TE has undergone intense transposition within the viral population. Finally, we provide evidence that many TEs found integrated in viral genomes (15/41) have been horizontally transferred in insects. Altogether, our results indicate that multiple large dsDNA viruses have the capacity to shuttle TEs in insects and they underline the potential of viruses to act as vectors of horizontal transfer of TEs. Furthermore, the finding that TEs can transpose between viral genomes of a viral species sets viruses as possible new niches in which TEs can persist and evolve.

Comprehensive Characterization of Transposable Elements in the Acute Myeloid Leukemia Transcriptome

Over half of the human genome is comprised of genomic transposable elements (TEs). TE expression has been implicated in cancer. It has been shown to induce genomic instability and more recently, a potential beneficial effect by causing cancer cell death by activating the viral recognition pathways. Recently we showed that the expression of distinct classes of TE such as SINEs and LTRs, but not LINEs, are suppressed in acute myeloid leukemia (AML) stem cells. Hence, it is likely that different classes and subtypes of TE perform diverse functions in cancer.Several large-scale studies characterized the transcriptome of numerous cancer types, and improved our understanding of oncogenesis in cancer. However, the expression of TEs has not been analyzed in these impactful studies. We developed a tool, Arkas, to comprehensively characterize the expression of TEs, accompanying coding gene expression, and gene pathways in large transcriptome databases. Using this we investigated the dysregulation of TEs in 178 AML transcriptome from The Cancer Genome Atlas.We first investigated the effect of recurrent AML mutations on TE dysregulation. In order to disentangle the effect of recurrent mutations on the TE expression, we combined mutation with TE expression in a multivariate linear model. Using a mutation-centric view, we observed that the number TEs significantly differentially expressed with respect to mutational drivers varied significantly with different mutations (Figure A). TE dysregulation was most affected by TP53 mutation, showing downregulation. MLL gene rearrangement was the second most, with mostly upregulated TE expression.Regulation of TE expression and its downstream effects is not fully understood. In the past, assumptions have been made generalizing the function of TE as a whole without understanding the diversity in TE. In addition to directly activating the interferon pathway and DNA damage, TEs are key regulators of coding gene expression, hence likely exhibiting diverse function. In order to gain insight into this, we correlated the expression of TE biotypes with coding gene networks. The coding genes were first clustered together (to form modules) based on their co-expression profile. 50 modules were formed using this method (Figure D). These modules are likely co-regulated and/or functionally related. These modules were then correlated with the expression of specific TE types. Distinct correlation pattern emerged with coding gene networks; positively associated modules with LINE1 were negatively associated with LTR elements including ERVs. Within the LTR class, distinct modules were associated with each ERV types. Though these findings do not prove direct functional relationship between TE types and the coding gene networks, they support diverse functions of exhibited by each TE type and subtype.Finally, we investigated the role of TE in predicting prognosis in AML. Using a penalized cox regression analysis, which included nested 10-fold cross validation, we examined the TE expression association to survival prediction. We identified 17 TE transcripts able to predict prognosis in AML (Figure C), classifying them into good-risk (N=77) and bad risk AML (N=101) (P=0.0011) (Figure B). Interestingly, this signature was able to further divide good-risk AML (N=87, identified based on gene mutations) into good-risk (N=41) and bad-risk AML (N=46) (P=0.026), suggesting TE predicts prognosis independently from gene mutations. 12 (9 of which were ERVs) of the 17 TEs predicting prognosis were associated with good prognosis. 5 TEs including L1M3B_5, which is LINE1, predicted worse prognosis. Interestingly, 3 of the good risk TEs were suppressed in TP53 mutation. In summary, this is the first study to comprehensively characterize TE expression in cancer. TE expression in AML shows distinct expression pattern associated with specific mutations and coding gene networks. TE expression offers independent prognostic information. These finding suggest the need for future mechanistic characterization of the role of distinct TEs in AML. [Display omitted] DisclosuresNo relevant conflicts of interest to declare.