A new subpopulation of polar bears (Ursus maritimus) was recently discovered in the South-East of Greenland. This isolated colony inhabits a warmer climate zone, akin to the predicted future environments of polar bears with vastly reduced sea ice habitats, rendering this population of bears particularly important. Over two-thirds of polar bears will be extinct by 2050 with total extinction predicted by the end of this century, therefore understanding possible mechanisms of adaptation via genomic analyses and preservation are critical. Transposable elements (TEs) are mobile elements that may play a role in an adaptive response to environmental challenges. We analysed transcriptome data from polar bear sub-populations in cooler North-East (NEG) and warmer South-East Greenland (SEG) to compare TE activity between the two populations and its correlation with temperature and associated changes in gene expression. We identified activity hotspots in the genome of regions with significantly differentially expressed TEs. LINE family TEs were the most abundant, and most differentially expressed and divergent in the SEG population compared to reference TEs. We report a significant shift in TE activity and age, with younger more abundant TEs in the SEG populations. Differentially expressed genes in SEG populations were linked to Foxo signalling, ageing and metabolic pathways. Our results provide insights into how a genomic response at the TE level may allow the SEG subpopulations to adapt and survive to climate change and provides a useful resource for conservation in polar bears.

BackgroundThe recognition that transposable elements (TEs) play important roles in many biological processes has elicited growing interest in analyzing sequencing data derived from this ‘dark genome’. This goal is complicated by the highly repetitive nature of these sequences in genomes, requiring the deployment of several problem-specific tools as well as the curation of appropriate genome annotations. This pipeline aims to make the analysis of TE sequences and their expression more generally accessible.ResultsThe TE-Seq pipeline conducts an end-to-end analysis of RNA sequencing data, examining both genes and TEs, and is compatible with most eukaryotic species. It implements computational methods tailor-made for TEs, and produces a comprehensive analysis of TE expression at both the level of the individual element and at the TE clade level. Furthermore, if supplied with long-read DNA sequencing data, it is able to assess TE expression from non-reference (polymorphic) loci. As a demonstration, we analyzed proliferating, early senescent, and late senescent human lung fibroblast RNA-Seq data, and created a custom reference genome and annotations for this cell strain using Nanopore sequencing data. We found that several retrotransposable element clades were upregulated in senescence, which included non-reference, intact, and potentially active elements.ConclusionsTE-Seq is made available as a Snakemake pipeline which can be obtained at https://github.com/maxfieldk/TE-Seq.Supplementary InformationThe online version contains supplementary material available at 10.1186/s13100-025-00381-w.

Tn7 mobile genetic elements are known for their sophisticated target-site selection mechanisms and, in some cases, programmability. Recognition of target sites is mediated by designated transposon-encoded proteins and modulated by host factor proteins. In the case of the CRISPR-associated Tn7 elements from the type V-K, the ribosomal protein uS15 is an integral component of recruitment complex that promotes R-loop completion. Previous biochemical work also revealed that the ribosomal protein uL29 and the acyl carrier protein (ACP) influence Tn7 transposition frequency in vitro. However, how uL29 and ACP regulate the formation of the Tn7 targeting complex remains unclear. The prototypical Tn7 element encodes a heteromeric transposase (TnsA, TnsB), a AAA + adaptor (TnsC), and two target-site selection proteins (TnsD and TnsE). TnsD targets a highly conserved site at the end of the glmS gene (attTn7). However, poor protein stability has precluded the molecular characterization of how TnsD recognizes its target site. Here, we show that ACP and uL29 interact with the C-terminal region of TnsD through reciprocal electrostatic interactions, in turn, mitigating its tendency to aggregate. Additionally, we identify the uL29 and ACP residues that mediate the interaction with TnsD and stimulate DNA binding. These results unveil unique features of the TnsD-mediated target-site selection complex.Supplementary InformationThe online version contains supplementary material available at 10.1186/s13100-025-00369-6.

BackgroundTransposable Elements (TEs) represent a sizeable amount of mammalian genomes, providing regulatory sequences involved in shaping gene expression patterns. NF-Y is a Transcription factor -TF- trimer that binds to the CCAAT box, belonging to a selected group implicated in determining initiation of coding and noncoding RNAs.ResultsWe focus on NF-Y TE locations in 8 human and 8 mouse cells. Binding is exclusive for retroviral LTR12, MLT1 and MER in human and RLTR10 and IAPLTR in mouse cells. Cobinding and analysis of the DNA matrices signal enrichment of distinct TFs neighboring CCAAT in the three TE classes: MAFK/F/G in LTR12 and USF1/2 in MLT1 with precise alignment of sites, PKNOX1, MEIS2, PBX2/3 TALE TFs in MER57. The presence of “epigenetic” marks in human cells indicate prevalent co-association with open chromatin in MER, closed in LTR12 and mixed in MLT1. Based on chromatin features, these locations are mostly marked as enhancers, as confirmed by analysis of loci predicted to generate eRNAs.ConclusionsThese results are discussed in the context of functional data, suggesting a complex -positive and potentially-negative role of NF-Y on distinct classes of repetitive sequences.

Transposable elements (TEs) drive genome evolution and can affect gene expression through diverse mechanisms. In breast cancer, disrupted regulation of TE sequences may facilitate tumor-specific transcriptomic alterations. We examine 142,514 full-length isoforms derived from long-read RNA sequencing (LR-seq) of 30 breast samples to investigate the effects of TEs on the breast cancer transcriptome. Approximately half of these isoforms contain TE sequences, and these contribute to half of the novel annotated splice junctions. We quantify splicing of these LR-seq derived isoforms in 1,135 breast tumors from The Cancer Genome Atlas (TCGA) and 1,329 healthy tissue samples from the Genotype-Tissue Expression (GTEx), and find 300 TE-overlapping tumor-specific splicing events. Some splicing events are enriched in specific breast cancer subtypes – for example, a TE-driven transcription start site upstream of ERBB2 in HER2 + tumors, and several TE-mediated splicing events are associated with patient survival and poor prognosis. The full-length sequences we capture with LR-seq reveal thousands of isoforms with signatures of RNA editing, including a novel isoform belonging to RHOA; a gene previously implicated in tumor progression. We utilize our full-length isoforms to discover polymorphic TE insertions that alter splicing and validate one of these events in breast cancer cell lines. Together, our results demonstrate the widespread effects of dysregulated TEs on breast cancer transcriptomes and highlight the advantages of long-read isoform sequencing for understanding TE biology. TE-derived isoforms may alter the expression of genes important in cancer and can potentially be used as novel, disease-specific therapeutic targets or biomarkers.One sentence summary: Transposable elements generate alternative isoforms and alter post-transcriptional regulation in human breast cancer.

BackgroundSystemic lupus erythematosus (SLE) is a chronic autoimmune disease with an unpredictable course of recurrent exacerbations alternating with more stable disease. SLE is characterized by broad immune activation and autoantibodies against double-stranded DNA and numerous proteins that exist in cells as aggregates with nucleic acids, such as Ro60, MOV10, and the L1 retrotransposon-encoded ORF1p.ResultsHere we report that these 3 proteins are co-expressed and co-localized in a subset of SLE granulocytes and are concentrated in cytosolic dots that also contain DNA: RNA heteroduplexes and the DNA sensor ZBP1, but not cGAS. The DNA: RNA heteroduplexes vanished from the neutrophils when they were treated with a selective inhibitor of the L1 reverse transcriptase. We also report that ORF1p granules escape neutrophils during the extrusion of neutrophil extracellular traps (NETs) and, to a lesser degree, from neutrophils dying by pyroptosis, but not apoptosis.ConclusionsThese results bring new insights into the composition of ORF1p granules in SLE neutrophils and may explain, in part, why proteins in these granules become targeted by autoantibodies in this disease.

Eukaryotic retroelements are generally divided into two classes: long terminal repeat (LTR) retrotransposons and non-LTR retrotransposons. A third class of eukaryotic retroelement, the Penelope-like elements (PLEs), has been well-characterized bioinformatically, but relatively little is known about the transposition mechanism of these elements. PLEs share some features with the R2 retrotransposon from Bombyx mori, which uses a target-primed reverse transcription (TPRT) mechanism, but their distinct phylogeny suggests PLEs may utilize a novel mechanism of mobilization. Using protein purified from E. coli, we report unique in vitro properties of a PLE from the green anole (Anolis carolinensis), revealing mechanistic aspects not shared by other retrotransposons. We found that reverse transcription is initiated at two adjacent sites within the transposon RNA that is not homologous to the cleaved DNA, a feature that is reflected in the genomic “tail” signature shared between and unique to PLEs. Our results for the first active PLE in vitro provide a starting point for understanding PLE mobilization and biology.

How novel protein functions are acquired is a central question in molecular biology. Key paths to novelty include gene duplications, recombination or horizontal acquisition. Transposable elements (TEs) are increasingly recognized as a major source of novel domain-encoding sequences. However, the impact of TE coding sequences on the evolution of the proteome remains understudied. Here, we analyzed 1237 genomes spanning the phylogenetic breadth of the fungal kingdom. We scanned proteomes for evidence of co-occurrence of TE-derived domains along with other conventional protein functional domains. We detected more than 13,000 predicted proteins containing potentially TE-derived domain, of which 825 were identified in more than five genomes, indicating that many host-TE fusions may have persisted over long evolutionary time scales. We used the phylogenetic context to identify the origin and retention of individual TE-derived domains. The most common TE-derived domains are helicases derived from Academ, Kolobok or Helitron. We found putative TE co-options at a higher rate in genomes of the Saccharomycotina, providing an unexpected source of protein novelty in these generally TE depleted genomes. We investigated in detail a candidate host-TE fusion with a heterochromatic transcriptional silencing function that may play a role in TE and gene regulation in ascomycetes. The affected gene underwent multiple full or partial losses within the phylum. Overall, our work establishes a kingdom-wide view of putative host-TE fusions and facilitates systematic investigations of candidate fusion proteins.

In November 2022 the first Dark Genome Symposium was held in Boston, USA. The meeting was hosted by Rome Therapeutics and Enara Bio, two biotechnology companies working on translating our growing understanding of this vast genetic landscape into therapies for human disease. The spirit and ambition of the meeting was one of shared knowledge, looking to strengthen the network of researchers engaged in the field. The meeting opened with a welcome from Rosana Kapeller and Kevin Pojasek followed by a first session of field defining talks from key academics in the space. A series of panels, bringing together academia and industry views, were then convened covering a wide range of pertinent topics. Finally, Richard Young and David Ting gave their views on the future direction and promise for patient impact inherent in the growing understanding of the Dark Genome.

Accumulating evidence suggests that endogenous retroviruses (ERVs) play an important role in the host response to infection and the development of disease. By analyzing ChIP-sequencing data sets, we show that SARS-CoV-2 infection induces H3K27 acetylation of several loci within the LTR69 subfamily of ERVs. Using functional assays, we identified one SARS-CoV-2-activated LTR69 locus, termed Dup69, which exhibits regulatory activity and is responsive to the transcription factors IRF3 and p65/RELA. LTR69_Dup69 is located about 500 bp upstream of a long non-coding RNA gene (ENSG00000289418) and within the PTPRN2 gene encoding a diabetes-associated autoantigen. Both ENSG00000289418 and PTPRN2 showed a significant increase in expression upon SARS-CoV-2 infection. Thus, our study sheds light on the interplay of exogenous with endogenous viruses and helps to understand how ERVs regulate gene expression during infection.