Comprehensive analysis of the Xya riparia genome uncovers the dominance of DNA transposons, LTR/Gypsy elements, and their evolutionary dynamics.

Transposable elements (TEs) are DNA sequences that can move or replicate within a genome, and their study has become increasingly important in understanding genome evolution and function. The Tridactylidae family, including Xya riparia (pygmy mole cricket), harbors a variety of transposable elements (TEs) that have been insufficiently investigated. Further research is required to fully understand their diversity and evolutionary characteristics. Hence, we conducted a comprehensive repeatome analysis of X. riparia species using the chromosome-level assembled genome. The study aimed to comprehensively analyze the abundance, distribution, and age of transposable elements (TEs) in the genome. The results indicated that the genome was 1.67 Gb, with 731.63 Mb of repetitive sequences, comprising 27% of Class II (443.25 Mb), 16% of Class I (268.45 Mb), and 1% of unknown TEs (19.92 Mb). The study found that DNA transposons dominate the genome, accounting for approximately 60% of the total repeat size, with retrotransposons and unknown elements accounting for 37% and 3% of the genome, respectively. The members of the Gypsy superfamily were the most abundant amongst retrotransposons, accounting for 63% of them. The transposable superfamilies (LTR/Gypsy, DNA/nMITE, DNA/hAT, and DNA/Helitron) collectively constituted almost 70% of the total repeat size of all six chromosomes. The study further unveiled a significant linear correlation (Pearson correlation: r = 0.99, p-value = 0.00003) between the size of the chromosomes and the repetitive sequences. The average age of DNA transposon and retrotransposon insertions ranges from 25 My (million years) to 5 My. The satellitome analysis discovered 13 satellite DNA families that comprise about 0.15% of the entire genome. In addition, the transcriptional analysis of TEs found that DNA transposons were more transcriptionally active than retrotransposons. Overall, the study suggests that the genome of X. riparia is complex, characterized by a substantial portion of repetitive elements. These findings not only enhance our understanding of TE evolution within the Tridactylidae family but also provide a foundation for future investigations into the genomic intricacies of related species.

The Ty1 retrotransposon harbors a DNA region that performs dual functions as both a gene silencing and chromatin insulator.

In various eukaryotic kingdoms, long terminal repeat (LTR) retrotransposons repress transcription by infiltrating heterochromatin generated within their elements. In contrast, the budding yeast LTR retrotransposon Ty1 does not itself undergo transcriptional repression, although it is capable of repressing the transcription of the inserted genes within it. In this study, we identified a DNA region within Ty1 that exerts its silencing effect via sequence orientation. We identified a DNA region within the Ty1 group-specific antigen (GAG) gene that causes gene silencing, termed GAG silencing (GAGsi), in which the silent chromatin in the GAGsi region is created by euchromatin-specific histone modifications. A characteristic inverted repeat (IR) sequence is present at the 5' end of this region, forming a chromatin boundary between promoter-specific chromatin upstream of the IR sequence and silent chromatin downstream of the IR sequence. In addition, Esc2 and Rad57, which are involved in DNA repair, were required for GAGsi silencing. Finally, the chromatin boundary was required for the transcription of Ty1 itself. Thus, the GAGsi sequence contributes to the creation of a chromatin environment that promotes Ty1 transcription.

The Families of Non-LTR Transposable Elements within Neritimorpha and Other Gastropoda.

Repeated sequences, especially transposable elements (TEs), are known to be abundant in some members of the important invertebrate class Gastropoda. TEs that do not have long terminal repeated sequences (non-LTR TEs) are frequently the most abundant type but have not been well characterised in any gastropod. Despite this, sequences in draft gastropod genomes are often described as non-LTR TEs, but without identification to family type. This study was conducted to characterise non-LTR TEs in neritimorph snails, using genomic skimming surveys of three species and the recently published draft genome of Theodoxus fluviatilis. Multiple families of non-LTR TEs from the I, Jockey, L1, R2 and RTE superfamilies were found, although there were notably few representatives of the first of these, which is nevertheless abundant in other Gastropoda. Phylogenetic analyses of amino acid sequences of the reverse transcriptase domain from the elements ORF2 regions found considerable interspersion of representatives of the four neritimorph taxa within non-LTR families and sub-families. In contrast, phylogenetic analyses of sequences from the elements' ORF1 region resolved the representatives from individual species as monophyletic. However, using either region, members of the two species of the Neritidae were closely related, suggesting their potential for investigation of phyletic evolution at the family level.

NF-κB1 deficiency promotes macrophage-derived adrenal tumors but decreases neurofibromas in HTLV-I LTR-Tax transgenic mice.

Human T-cell leukemia virus type I (HTLV-I) is an oncogenic virus whose infection can cause diverse diseases, most notably adult T-cell leukemia/lymphoma (ATL or ATLL), an aggressive and fatal malignancy of CD4 T cells. The oncogenic ability of HTLV-I is mostly attributed to the viral transcriptional transactivator Tax. Tax alone is sufficient to induce specific tumors in mice depending on the promotor used to drive Tax expression, thereby being used to understand HTLV-I tumorigenesis and model the tumor types developed in Tax transgenic mice. Tax exerts its oncogenic role predominantly by activating the cellular transcription factor NF-κB. Here, we report that genetic deletion of NF-κB1, the prototypic member of the NF-κB family, promotes adrenal medullary tumors but suppresses neurofibromas in mice with transgenic Tax driven by the HTLV-I Long Terminal Repeat (LTR) promoter. The adrenal tumors are derived from macrophages. Neoplastic macrophages also infiltrate the spleen and lymph nodes, causing splenomegaly and lymphadenopathy in mice. Nevertheless, the findings could be human relevant, because macrophages are important target cells of HTLV-I infection and serve as a virus reservoir in vivo. Moreover, the spleen, lymph nodes and adrenal glands are the most common sites of tumor cell infiltration in HTLV-I-infected patients. These data provide new mechanistic insights into the complex interaction between Tax and NF-κB, therefore improving our understanding of HTLV-I oncogenic pathogenesis. They also expand our knowledge and establish a new animal model of macrophage neoplasms and adrenal tumors.

KSHV genome harbors both constitutive and lytically induced enhancers.

Kaposi's sarcoma-associated herpesvirus (KSHV) belongs to the gamma-herpesvirus family and is a well-known human oncogenic virus. In infected cells, the viral genome of 165 kbp is circular DNA wrapped in chromatin. The tight control of gene expression is critical for latency, the transition into the lytic phase, and the development of viral-associated malignancies. Distal cis-regulatory elements, such as enhancers and silencers, can regulate gene expression in a position- and orientation-independent manner. Open chromatin is another characteristic feature of enhancers. To systematically search for enhancers, we cloned all the open chromatin regions in the KSHV genome downstream of the luciferase gene and tested their enhancer activity in infected and uninfected cells. A silencer was detected upstream of the latency-associated nuclear antigen promoter. Two constitutive enhancers were identified in the K12p-OriLyt-R and ORF29 Intron regions, where ORF29 Intron is a tissue-specific enhancer. The following promoters: OriLyt-L, PANp, ALTp, and the terminal repeats (TRs) acted as lytically induced enhancers. The expression of the replication and transcription activator (RTA), the master regulator of the lytic cycle, was sufficient to induce the activity of lytic enhancers in uninfected cells. We propose that the TRs that span about 24 kbp region serve as a "viral super-enhancer" that integrates the repressive effect of the latency-associated nuclear antigen (LANA) with the activating effect of RTA. Utilizing CRISPR activation and interference techniques, we determined the connections between these enhancers and their regulated genes. The silencer and enhancers described here provide an additional layer to the complex gene regulation of herpesviruses.IMPORTANCEIn this study, we performed a systematic functional assay to identify cis-regulatory elements within the genome of the oncogenic herpesvirus, Kaposi's sarcoma-associated herpesvirus (KSHV). Similar to other herpesviruses, KSHV presents both latent and lytic phases. Therefore, our assays were performed in uninfected cells, during latent infection, and under lytic conditions. We identified two constitutive enhancers, one of which seems to be a tissue-specific enhancer. In addition, four lytically induced enhancers, which are all responsive to the replication and transcription activator (RTA), were identified. Furthermore, a silencer was identified between the major latency promoter and the lytic gene locus. Utilizing CRISPR activation and interference techniques, we determined the connections between these enhancers and their regulated genes. The terminal repeats, spanning a region of about 24 kbp, seem like a "viral super-enhancer" that integrates the repressive effect of the latency-associated nuclear antigen (LANA) with the activating effect of RTA to regulate latency to lytic transition.

Anaerobic fungi contain abundant, diverse, and transcriptionally active Long Terminal Repeat retrotransposons.

Long Terminal Repeat (LTR) retrotransposons are a class of repetitive elements that are widespread in the genomes of plants and many fungi. LTR retrotransposons have been associated with rapidly evolving gene clusters in plants and virulence factor transfer in fungal-plant parasite-host interactions. We report here the abundance and transcriptional activity of LTR retrotransposons across several species of the early-branching Neocallimastigomycota, otherwise known as the anaerobic gut fungi (AGF). The ubiquity of LTR retrotransposons in these genomes suggests key evolutionary roles in these rumen-dwelling biomass degraders, whose genomes also contain many enzymes that are horizontally transferred from other rumen-dwelling prokaryotes. Up to 10% of anaerobic fungal genomes consist of LTR retrotransposons, and the mapping of sequences from LTR retrotransposons to transcriptomes shows that the majority of clusters are transcribed, with some exhibiting expression greater than 104 reads per kilobase million mapped reads (rpkm). Many LTR retrotransposons are strongly differentially expressed upon heat stress during fungal cultivation, with several exhibiting a nearly three-log10 fold increase in expression, whereas growth substrate variation modulated transcription to a lesser extent. We show that some LTR retrotransposons contain carbohydrate-active enzymes (CAZymes), and the expansion of CAZymes within genomes and among anaerobic fungal species may be linked to retrotransposon activity. We further discuss how these widespread sequences may be a source of promoters and other parts towards the bioengineering of anaerobic fungi.

rAAV capsid mutants eliminate leaky expression from DNA donor template for homologous recombination.

Precise genomic editing through the combination of CRISPR/Cas systems and recombinant adeno-associated virus (rAAV)-delivered homology directed repair (HDR) donor templates represents a powerful approach. However, the challenge of effectively suppressing leaky transcription from the rAAV vector, a phenomenon associated to cytotoxicity, persists. In this study, we demonstrated substantial promoter activities of various homology arms and inverted terminal repeats (ITR). To address this issue, we identified a novel rAAV variant, Y704T, which not only yields high-vector quantities but also effectively suppresses in cis mRNA transcription driven by a robust promoter. The Y704T variant maintains normal functionality in receptor interaction, intracellular trafficking, nuclear entry, uncoating, and second-strand synthesis, while specifically exhibiting defects in transcription. Importantly, this inhibitory effect is found to be independent of ITR, promoter types, and RNA polymerases. Mechanistic studies unveiled the involvement of Valosin Containing Protein (VCP/p97) in capsid-mediated transcription repression. Remarkably, the Y704T variant delivers HDR donor templates without compromising DNA replication ability and homologous recombination efficiency. In summary, our findings enhance the understanding of capsid-regulated transcription and introduce novel avenues for the application of the rAAV-CRISPR/Cas9 system in human gene therapy.

Karyotype and LTR-RTs analysis provide insights into oak genomic evolution.

BACKGROUND: Whole-genome duplication and long terminal repeat retrotransposons (LTR-RTs) amplification in organisms are essential factors that affect speciation, local adaptation, and diversification of organisms. Understanding the karyotype projection and LTR-RTs amplification could contribute to untangling evolutionary history. This study compared the karyotype and LTR-RTs evolution in the genomes of eight oaks, a dominant lineage in Northern Hemisphere forests. RESULTS: Karyotype projections showed that chromosomal evolution was relatively conservative in oaks, especially on chromosomes 1 and 7. Modern oak chromosomes formed through multiple fusions, fissions, and rearrangements after an ancestral triplication event. Species-specific chromosomal rearrangements revealed fragments preserved through natural selection and adaptive evolution. A total of 441,449 full-length LTR-RTs were identified from eight oak genomes, and the number of LTR-RTs for oaks from section Cyclobalanopsis was larger than in other sections. Recent amplification of the species-specific LTR-RTs lineages resulted in significant variation in the abundance and composition of LTR-RTs among oaks. The LTR-RTs insertion suppresses gene expression, and the suppressed intensity in gene regions was larger than in promoter regions. Some centromere and rearrangement regions indicated high-density peaks of LTR/Copia and LTR/Gypsy. Different centromeric regional repeat units (32, 78, 79 bp) were detected on different Q. glauca chromosomes. CONCLUSION: Chromosome fusions and arm exchanges contribute to the formation of oak karyotypes. The composition and abundance of LTR-RTs are affected by its recent amplification. LTR-RTs random retrotransposition suppresses gene expression and is enriched in centromere and chromosomal rearrangement regions. This study provides novel insights into the evolutionary history of oak karyotypes and the organization, amplification, and function of LTR-RTs.

Adaptive expansion of ERVK solo-LTRs is associated with Passeriformes speciation events.

Endogenous retroviruses (ERVs) are ancient retroviral remnants integrated in host genomes, and commonly deleted through unequal homologous recombination, leaving solitary long terminal repeats (solo-LTRs). This study, analysing the genomes of 362 bird species and their reptilian and mammalian outgroups, reveals an unusually higher level of solo-LTRs formation in birds, indicating evolutionary forces might have purged ERVs during evolution. Strikingly in the order Passeriformes, and especially the parvorder Passerida, endogenous retrovirus K (ERVK) solo-LTRs showed bursts of formation and recurrent accumulations coinciding with speciation events over past 22 million years. Moreover, our results indicate that the ongoing expansion of ERVK solo-LTRs in these bird species, marked by high transcriptional activity of ERVK retroviral genes in reproductive organs, caused variation of solo-LTRs between individual zebra finches. We experimentally demonstrated that cis-regulatory activity of recently evolved ERVK solo-LTRs may significantly increase the expression level of ITGA2 in the brain of zebra finches compared to chickens. These findings suggest that ERVK solo-LTRs expansion may introduce novel genomic sequences acting as cis-regulatory elements and contribute to adaptive evolution. Overall, our results underscore that the residual sequences of ancient retroviruses could influence the adaptive diversification of species by regulating host gene expression.

Transposition of HOPPLA in siRNA-deficient plants suggests a limited effect of the environment on retrotransposon mobility in Brachypodium distachyon.

Long terminal repeat retrotransposons (LTR-RTs) are powerful mutagens regarded as a major source of genetic novelty and important drivers of evolution. Yet, the uncontrolled and potentially selfish proliferation of LTR-RTs can lead to deleterious mutations and genome instability, with large fitness costs for their host. While population genomics data suggest that an ongoing LTR-RT mobility is common in many species, the understanding of their dual role in evolution is limited. Here, we harness the genetic diversity of 320 sequenced natural accessions of the Mediterranean grass Brachypodium distachyon to characterize how genetic and environmental factors influence plant LTR-RT dynamics in the wild. When combining a coverage-based approach to estimate global LTR-RT copy number variations with mobilome-sequencing of nine accessions exposed to eight different stresses, we find little evidence for a major role of environmental factors in LTR-RT accumulations in B. distachyon natural accessions. Instead, we show that loss of RNA polymerase IV (Pol IV), which mediates RNA-directed DNA methylation in plants, results in high transcriptional and transpositional activities of RLC_BdisC024 (HOPPLA) LTR-RT family elements, and that these effects are not stress-specific. This work supports findings indicating an ongoing mobility in B. distachyon and reveals that host RNA-directed DNA methylation rather than environmental factors controls their mobility in this wild grass model.

Systematic evaluation of retroviral LTRs as cis-regulatory elements in mouse embryos.

In mammals, many retrotransposons are de-repressed during zygotic genome activation (ZGA). However, their functions in early development remain elusive largely due to the challenge to simultaneously manipulate thousands of retrotransposon insertions in embryos. Here, we applied CRISPR interference (CRISPRi) to perturb the long terminal repeat (LTR) MT2_Mm, a well-known ZGA and totipotency marker that exists in ∼2,667 insertions throughout the mouse genome. CRISPRi robustly perturbed 2,485 (∼93%) MT2_Mm insertions and 1,090 (∼55%) insertions of the closely related MT2C_Mm in 2-cell embryos. Remarkably, such perturbation caused downregulation of hundreds of ZGA genes and embryonic arrest mostly at the morula stage. Mechanistically, MT2 LTRs are globally enriched for open chromatin and H3K27ac and function as promoters/enhancers downstream of OBOX/DUX proteins. Thus, we not only provide direct evidence to support the functional importance of MT2 activation in development but also systematically define cis-regulatory function of MT2 in embryos by integrating functional perturbation and multi-omic analyses.

Kaposi's sarcoma-associated herpesvirus terminal repeat regulates inducible lytic gene promoters.

The Kaposi's sarcoma-associated herpesvirus (KSHV) genome consists of an approximately 140-kb unique coding region flanked by 30-40 copies of a 0.8-kb terminal repeat (TR) sequence. A gene enhancer recruits transcription-related enzymes by having arrays of transcription factor binding sites. Here, we show that KSHV TR possesses transcription regulatory function with latency-associated nuclear antigen (LANA). Cleavage under targets and release using nuclease demonstrated that TR fragments were occupied by LANA-interacting histone-modifying enzymes in naturally infected cells. The TR was enriched with histone H3K27 acetylation (H3K27Ac) and H3K4 tri-methylation (H3K4me3) modifications and also expressed nascent RNAs. The sites of H3K27Ac and H3K4me3 modifications were also conserved in the KSHV unique region among naturally infected primary effusion lymphoma cells. KSHV origin of lytic replication (Ori-Lyt) showed similar protein and histone modification occupancies with that of TR. In the Ori-Lyt region, the LANA and LANA-interacting proteins colocalized with an H3K27Ac-modified nucleosome along with paused RNA polymerase II. The KSHV transactivator KSHV replication and transcription activator (K-Rta) recruitment sites franked the LANA-bound nucleosome, and reactivation evicted the LANA-bound nucleosome. Including TR fragments in reporter plasmid enhanced inducible viral gene promoter activities independent of the orientations. In the presence of TR in reporter plasmids, K-Rta transactivation was drastically increased, while LANA acquired the promoter repression function. KSHV TR, therefore, functions as an enhancer for KSHV inducible genes. However, in contrast to cellular enhancers bound by multiple transcription factors, perhaps the KSHV enhancer is predominantly regulated by the LANA nuclear body.IMPORTANCEEnhancers are a crucial regulator of differential gene expression programs. Enhancers are the cis-regulatory sequences determining target genes' spatiotemporal and quantitative expression. Here, we show that Kaposi's sarcoma-associated herpesvirus (KSHV) terminal repeats fulfill the enhancer definition for KSHV inducible gene promoters. The KSHV enhancer is occupied by latency-associated nuclear antigen (LANA) and its interacting proteins, such as CHD4. Neighboring terminal repeat (TR) fragments to lytic gene promoters drastically enhanced KSHV replication and transcription activator and LANA transcription regulatory functions. This study, thus, proposes a new latency-lytic switch model in which TR accessibility to the KSHV gene promoters regulates viral inducible gene expression.

TNRC18 engages H3K9me3 to mediate silencing of endogenous retrotransposons.

Trimethylation of histone H3 lysine 9 (H3K9me3) is crucial for the regulation of gene repression and heterochromatin formation, cell-fate determination and organismal development1. H3K9me3 also provides an essential mechanism for silencing transposable elements1-4. However, previous studies have shown that canonical H3K9me3 readers (for example, HP1 (refs. 5-9) and MPP8 (refs. 10-12)) have limited roles in silencing endogenous retroviruses (ERVs), one of the main transposable element classes in the mammalian genome13. Here we report that trinucleotide-repeat-containing 18 (TNRC18), a poorly understood chromatin regulator, recognizes H3K9me3 to mediate the silencing of ERV class I (ERV1) elements such as LTR12 (ref. 14). Biochemical, biophysical and structural studies identified the carboxy-terminal bromo-adjacent homology (BAH) domain of TNRC18 (TNRC18(BAH)) as an H3K9me3-specific reader. Moreover, the amino-terminal segment of TNRC18 is a platform for the direct recruitment of co-repressors such as HDAC-Sin3-NCoR complexes, thus enforcing optimal repression of the H3K9me3-demarcated ERVs. Point mutagenesis that disrupts the TNRC18(BAH)-mediated H3K9me3 engagement caused neonatal death in mice and, in multiple mammalian cell models, led to derepressed expression of ERVs, which affected the landscape of cis-regulatory elements and, therefore, gene-expression programmes. Collectively, we describe a new H3K9me3-sensing and regulatory pathway that operates to epigenetically silence evolutionarily young ERVs and exert substantial effects on host genome integrity, transcriptomic regulation, immunity and development.

m6A RNA demethylase AtALKBH9B promotes mobilization of a heat-activated long terminal repeat retrotransposon in Arabidopsis.

Transposons are mobile and ubiquitous DNA molecules that can cause vast genomic alterations. In plants, it is well documented that transposon mobilization is strongly repressed by DNA methylation; however, its regulation at the posttranscriptional level remains relatively uninvestigated. Here, we suggest that transposon RNA is marked by m6A RNA methylation and can be localized in stress granules (SGs). Intriguingly, SG-localized AtALKBH9B selectively demethylates a heat-activated retroelement, Onsen, and thereby releases it from spatial confinement, allowing for its mobilization. In addition, we show evidence that m6A RNA methylation contributes to transpositional suppression by inhibiting virus-like particle assembly and extrachromosomal DNA production. In summary, this study unveils a previously unknown role for m6A in the suppression of transposon mobility and provides insight into how transposons counteract the m6A-mediated repression mechanism by hitchhiking the RNA demethylase of the host.

Origin and Evolution of Plant Long Terminal Repeat Retrotransposons with Additional Ribonuclease H.

Retroviruses originated from long terminal repeat retrotransposons (LTR-RTs) through several structural adaptations. One such modification was the arrangement of an additional ribonuclease H (aRH) domain next to native RH, followed by degradation and subfunctionalization of the latter. We previously showed that this retrovirus-like structure independently evolved in Tat LTR-RTs in flowering plants, proposing its origin from sequential rearrangements of ancestral Tat structures identified in lycophytes and conifers. However, most nonflowering plant genome assemblies were not available at that time, therefore masking the history of aRH acquisition by Tat and challenging our hypothesis. Here, we revisited Tat's evolution scenario upon the aRH acquisition by covering most of the extant plant phyla. We show that Tat evolved and obtained aRH in an ancestor of land plants. Importantly, we found the retrovirus-like structure in clubmosses, hornworts, ferns, and gymnosperms, suggesting its ancient origin, broad propagation, and yet-to-be-understood benefit for the LTR-RTs' adaptation.

The Pleiotropic Effects of YBX1 on HTLV-1 Transcription.

HTLV-1 is an oncogenic human retrovirus and the etiologic agent of the highly aggressive ATL malignancy. Two viral genes, Tax and Hbz, are individually linked to oncogenic transformation and play an important role in the pathogenic process. Consequently, regulation of HTLV-1 gene expression is a central feature in the viral lifecycle and directly contributes to its pathogenic potential. Herein, we identified the cellular transcription factor YBX1 as a binding partner for HBZ. We found YBX1 activated transcription and enhanced Tax-mediated transcription from the viral 5' LTR promoter. Interestingly, YBX1 also interacted with Tax. shRNA-mediated loss of YBX1 decreased transcript and protein abundance of both Tax and HBZ in HTLV-1-transformed T-cell lines, as well as Tax association with the 5' LTR. Conversely, YBX1 transcriptional activation of the 5' LTR promoter was increased in the absence of HBZ. YBX1 was found to be associated with both the 5' and 3' LTRs in HTLV-1-transformed and ATL-derived T-cell lines. Together, these data suggest that YBX1 positively influences transcription from both the 5' and 3' promoter elements. YBX1 is able to interact with Tax and help recruit Tax to the 5' LTR. However, through interactions with HBZ, YBX1 transcriptional activation of the 5' LTR is repressed.

The blackcap (Sylvia atricapilla) genome reveals a recent accumulation of LTR retrotransposons.

Transposable elements (TEs) are mobile genetic elements that can move around the genome, and as such are a source of genomic variability. Based on their characteristics we can annotate TEs within the host genome and classify them into specific TE types and families. The increasing number of available high-quality genome references in recent years provides an excellent resource that will enhance the understanding of the role of recently active TEs on genetic variation and phenotypic evolution. Here we showcase the use of a high-quality TE annotation to understand the distinct effect of recent and ancient TE insertions on the evolution of genomic variation, within our study species the Eurasian blackcap (Sylvia atricapilla). We investigate how these distinct TE categories are distributed along the genome and evaluate how their coverage across the genome is correlated with four genomic features: recombination rate, gene coverage, CpG island coverage and GC content. We found within the recent TE insertions an accumulation of LTRs previously not seen in birds. While the coverage of recent TE insertions was negatively correlated with both GC content and recombination rate, the correlation with recombination rate disappeared and turned positive for GC content when considering ancient TE insertions.

Lineage-specific amplification and epigenetic regulation of LTR-retrotransposons contribute to the structure, evolution, and function of Fabaceae species.

BACKGROUND: Long terminal repeat (LTR)-retrotransposons (LTR-RTs) are ubiquitous and make up the majority of nearly all sequenced plant genomes, whereas their pivotal roles in genome evolution, gene expression regulation as well as their epigenetic regulation are still not well understood, especially in a large number of closely related species. RESULTS: Here, we analyzed the abundance and dynamic evolution of LTR-RTs in 54 species from an economically and agronomically important family, Fabaceae, and also selected two representative species for further analysis in expression of associated genes, transcriptional activity and DNA methylation patterns of LTR-RTs. Annotation results revealed highly varied proportions of LTR-RTs in these genomes (5.1%~68.4%) and their correlation with genome size was highly positive, and they were significantly contributed to the variance in genome size through species-specific unique amplifications. Almost all of the intact LTR-RTs were inserted into the genomes 4 Mya (million years ago), and more than 50% of them were inserted in the last 0.5 million years, suggesting that recent amplifications of LTR-RTs were an important force driving genome evolution. In addition, expression levels of genes with intronic, promoter, and downstream LTR-RT insertions of Glycine max and Vigna radiata, two agronomically important crops in Fabaceae, showed that the LTR-RTs located in promoter or downstream regions suppressed associated gene expression. However, the LTR-RTs within introns promoted gene expression or had no contribution to gene expression. Additionally, shorter and younger LTR-RTs maintained higher mobility and transpositional potential. Compared with the transcriptionally silent LTR-RTs, the active elements showed significantly lower DNA methylation levels in all three contexts. The distributions of transcriptionally active and silent LTR-RT methylation varied across different lineages due to the position of LTR-RTs located or potentially epigenetic regulation. CONCLUSION: Lineage-specific amplification patterns were observed and higher methylation level may repress the activity of LTR-RTs, further influence evolution in Fabaceae species. This study offers valuable clues into the evolution, function, transcriptional activity and epigenetic regulation of LTR-RTs in Fabaceae genomes.

Independent evolution of transposase and TIRs facilitated by recombination between Mutator transposons from divergent clades in maize.

Nearly all eukaryotes carry DNA transposons of the Robertson's Mutator (Mu) superfamily, a widespread source of genome instability and genetic variation. Despite their pervasive impact on host genomes, much remains unknown about the evolution of these transposons. Transposase recognition of terminal inverted repeats (TIRs) is thought to drive and constrain coevolution of MuDR transposase genes and TIRs. To address the extent of this relationship and its impact, we compared separate phylogenies of TIRs and MuDR gene sequences from Mu elements in the maize genome. Five major clades were identified. As expected, most Mu elements were bound by highly similar TIRs from the same clade (homomorphic type). However, a subset of elements contained dissimilar TIRs derived from divergent clades. These "heteromorphs" typically occurred in multiple copies indicating active transposition in the genome. In addition, analysis of internal sequences showed that exchanges between elements having divergent TIRs produced new mudra and mudrb gene combinations. In several instances, TIR homomorphs had been regenerated within a heteromorph clade with retention of distinctive internal MuDR sequence combinations. Results reveal that recombination between divergent clades facilitates independent evolution of transposase (mudra), transposase-binding targets (TIRs), and capacity for insertion (mudrb) of active Mu elements. This mechanism would be enhanced by the preference of Mu insertions for recombination-rich regions near the 5' ends of genes. We suggest that cycles of recombination give rise to alternating homo- and heteromorph forms that enhance the diversity on which selection for Mu fitness can operate.

Co-option of endogenous retroviruses through genetic escape from TRIM28 repression.

Endogenous retroviruses (ERVs) have rewired host gene networks. To explore the origins of co-option, we employed an active murine ERV, IAPEz, and an embryonic stem cell (ESC) to neural progenitor cell (NPC) differentiation model. Transcriptional silencing via TRIM28 maps to a 190 bp sequence encoding the intracisternal A-type particle (IAP) signal peptide, which confers retrotransposition activity. A subset of "escapee" IAPs (∼15%) exhibits significant genetic divergence from this sequence. Canonical repressed IAPs succumb to a previously undocumented demarcation by H3K9me3 and H3K27me3 in NPCs. Escapee IAPs, in contrast, evade repression in both cell types, resulting in their transcriptional derepression, particularly in NPCs. We validate the enhancer function of a 47 bp sequence within the U3 region of the long terminal repeat (LTR) and show that escapee IAPs convey an activating effect on nearby neural genes. In sum, co-opted ERVs stem from genetic escapees that have lost vital sequences required for both TRIM28 restriction and autonomous retrotransposition.