SEARCH · Search PubMed

Results for “LTR retrotransposons”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

15 recordsLinked to original sources

PubMed2026

Hide and seek: de novo identification in sugar beet reveals impact of non-autonomous LTR retrotransposons.

Plant genomes are filled with retrotransposons and their derivatives, constantly undergoing sequence diversification and structural rearrangement. Among them, short, non-autonomous retrotransposons lack full coding capacity and often form subfamilies. As a result, non-autonomous retrotransposons are incompletely identified in most to all genome assemblies.Here, we capitalize on our comprehensive understanding of the transposable element (TE) landscape in sugar beet (Beta vulgaris) to assess the extent of the blind spot for non-autonomous long terminal repeat (LTR) retrotransposons. This use case serves to answer if all of these sequences are derivatives of easier-to-identify full-length elements or if there is more variability that is currently overlooked.For this we applied a semi-automated structural discovery workflow followed by in-depth manual verification to characterize non-autonomous LTR retrotransposons in sugar beet. We retrieve more than 100 non-autonomous LTR retrotransposon families that lack complete autonomous coding capacity, including canonical terminal-repeat retrotransposons in miniature (TRIMs), elongated non-coding derivatives and families retaining fragmented coding remnants. The identified families span a broad range, including elements exceeding 15,000 bp in length and display evidence for reshuffling and modular evolution. Only a subset of families could be confidently linked to autonomous retrotransposons, showing sequence diversification within the non-autonomous LTR retrotransposon fraction beyond the autonomous genomic templates.We highlight that a large fraction of non-autonomous LTR retrotransposons is incompletely recovered with the current TE identification workflows, even if the output is well-curated and condensed into TE libraries and suggest procedures to remedy this gap. This study gives a genome-wide view into the non-autonomous LTR retrotransposon landscape of a single plant genome and highlights the importance of structure-based approaches for their identification and classification.

LTR retrotransposons↗

PubMed2026

Hijacking pre-tRNA enables LTR-retrotransposon-initiated constitutive heterochromatin formation.

Pericentric heterochromatin serves as a fundamental component of eukaryotic chromosomes, endowing specialized genomic architecture with broad functional consequences. Although it is universally marked by H3K9me3 modification, the underlying pericentric DNA sequences diverge substantially across species. Here, by leveraging a transposition reporter system combined with a genome-wide RNA interference (RNAi) screen, we identified a specialized mechanism for recruiting SUV39H methyltransferase to initiate pericentric heterochromatin formation. This pathway depends on a highly ordered complex comprising the Puf68, pre-transfer RNAs (tRNAs), and the primer binding site (PBS). Puf68 binds with high affinity to poly-U tracts in pre-tRNA 3' trailer, forming a Puf68/pre-tRNA complex that subsequently base-pairs with the PBS of nascent long terminal repeat (LTR)-retrotransposons. Through direct interaction, Puf68 recruits Su(var)3-9 to these regions, catalyzing H3K9 trimethylation. Notably, Puf68 is sufficient to initiate de novo heterochromatin assembly both at pericentric and ectopically integrated LTR-retrotransposon regions. Our findings not only uncover a previously unrecognized mechanism of heterochromatin initiation but also resolve a long-standing question of how hosts harness nascent LTR-retrotransposon transcripts.

Heterochromatin↗

PubMed2026

The piRNA pathway mediates transcriptional silencing of LTR retrotransposons in ovaries and somatic tissues of Aedes mosquitoes.

The PIWI-interacting RNA (piRNA) pathway preserves genomic integrity by suppressing transposable elements in animal germlines. Despite its well-established function in the animal germline, piRNAs and PIWI proteins are expressed in somatic tissues across arthropod species, and their functions outside the gonads remain poorly understood. Aedes albopictus mosquitoes express four PIWI genes, Piwi4, Piwi5, Piwi6, and Ago3, in both gonadal and somatic tissues. Here, we generated Piwi6 knockout (KO) Ae. albopictus cell lines and observed a substantial upregulation of long terminal repeat retrotransposons, including a full-length endogenous retrovirus that we named Aedes albopictus Endogenous Retrovirus-1 (AalERV1). Nascent RNA sequencing and Cleavage Under Targets and Tagmentation (CUT&Tag) analyses revealed that Piwi6 silences AalERV1 transcriptionally by guiding the deposition of the repressive H3K9me3 histone mark. Consistently, Piwi6 localized to both the cytoplasm and nucleus, with sequences in the intrinsically disordered region guiding nuclear translocation. Reintroduction of full-length GFP-Piwi6, but not a mutant GFP-Piwi6 defective in nuclear localization, rescued AalERV1 repression in Piwi6 KO cells. Importantly, Piwi6-mediated control of AalERV1 was recapitulated in vivo as Piwi6 knockdown increased AalERV1 expression in both ovaries and somatic tissues of Ae. albopictus mosquitoes. These results establish Aedes mosquitoes as a model to study nuclear PIWI functions and suggest that somatic piRNA-mediated transposon silencing is evolutionarily conserved across arthropod species.

Animals↗

PubMed2025

Recent Non-LTR Retrotransposon Activity Predicts Cancer Prevalence in Mammals.

Non-long terminal repeat retrotransposons (nLTRs), including long and short interspersed nuclear elements (L1 and SINEs), are the most abundant and active mobile elements in mammals. NLTRs play critical mutagenic and regulatory roles during oncogenesis in humans and model species. However, it is not known whether recent nLTR activity in the genome is related to the lifetime cancer risk of a species beyond humans and conventional model organisms. We examined whether recent nLTR activity predicts cancer prevalence across mammals using comparative analyses of de novo whole-genome repeat annotations from 55 species, each with over 20 published zoo pathology records. We quantified nLTR activity as the number of potentially active elements, their proximity to protein-coding genes and cancer gene orthologs (CGOs), and insertions within these genes. Across all three metrics, neoplasia prevalence was associated with both L1 and combined L1-SINE activity, while malignancy was linked exclusively to the L1-SINE predictors. This pattern suggests a complementary and escalating trajectory, where L1s contribute to early tumorigenic events, while SINE activity, driven by L1s, amplifies their impact and fuels the transition to malignancy. Moreover, genomes harboring more CGOs tended to exhibit higher neoplasia prevalence, and the number of fusion cancer genes was strongly correlated with the number of potentially active L1s across species. Our results further revealed a pattern wherein species with minimal cancer prevalence exhibit restricted activity of at least one major nLTR superfamily, suggesting that preserving genome stability through limited retrotransposition may serve as a protective mechanism against cancer.

Cancer Genes↗

PubMed2026

Retrotransposon activation during spermatogenesis achieves massive ecDNA biogenesis but rare integration.

Retrotransposon mobilization in germline cells enables the rewriting of genetic information to drive genome innovation, species evolution, and adaptation through the generation of de novo mutations. However, uncontrolled mobilization can cause DNA breaks and genome instability, often leading to sterility. How retrotransposon mobilization that can be retained for genome evolution persists despite negative outcomes of retrotransposon activity remains poorly understood. Here, we used Drosophila spermatogenesis as a model to investigate retrotransposon mobilization dynamics. Although many retrotransposon families are transcriptionally active, we found that the LTR retrotransposon nomad completes the full mobilization cascade (including mRNA export, protein translation, and reverse transcription) to produce double-stranded DNA (dsDNA) the most efficiently. Strikingly, despite successfully generating dsDNA, nomad rarely achieves genomic reintegration. Instead, its newly synthesized DNA predominantly forms extrachromosomal circular DNA (ecDNA). These findings show that retrotransposon-derived DNA largely remains as ecDNA. This could prevent widespread genomic integration during spermatogenesis, potentially preserving genome stability with the presence of limited retrotransposon activity.

Animals↗

PubMed2026

Revisiting the genome assembly of Lupinus species reveals differential diploidization after a shared whole-genome duplication.

Accurate genome assemblies are essential for comparative genomics, yet Hi-C-guided scaffolding can introduce structural errors that misrepresent chromosome architecture and bias evolutionary inferences. Here, we identified pervasive scaffolding errors-including artificial fusions, internal inversions, and incomplete contig mounting-in 2 previously published Lupinus genomes (L. cosentinii and L. digitatus) using a segmentation method based on long terminal repeat (LTR) retrotransposon density. We reassembled both genomes, producing chromosome-level references of 472.7 Mb (16 chromosomes) and 427.2 Mb (21 chromosomes), with BUSCO completeness >98.5%. Synteny validation and reapplication of LTR profiling confirmed that all prior errors were resolved. Using these corrected genomes together with 4 additional Lupinus species and 2 outgroup legumes, we investigated postpolyploid evolution. Synonymous substitution rate (Ks) analysis revealed a genus-specific whole-genome duplication (WGD) event (Ks = 0.17) shared by all 6 Lupinus species. The proportion of WGD-derived genes varied markedly, from 60% in L. digitatus to only 36% in L. mutabilis, indicating differential diploidization. While all species retained a core set of WGD duplicates enriched in cytoskeleton organization, ion transport, and defense responses, each exhibited lineage-specific functional trajectories: cell wall modification in L. cosentinii and L. digitatus, nitrogen metabolism in L. albus and L. angustifolius, flower development in L. luteus, and stress/lipid metabolism in L. mutabilis. Our corrected assemblies provide optimal references for Lupinus comparative genomics, and our findings demonstrate that a shared WGD event can lead to both conserved and highly divergent postpolyploid fates, likely underpinning adaptive diversification within the genus.

Lupinus↗

PubMed2026

Heat-responsive ONSEN long terminal repeats integrate heat shock factor motifs, DNA methylation and natural sequence variation in Arabidopsis.

ONSEN is a heat-activated Ty1/copia retrotransposon in Arabidopsis thaliana controlled by heat shock factors (HSFs) and epigenetic silencing. Heat shock element (HSE)-like sequences in ONSEN long terminal repeats (LTRs) contribute to heat responsiveness, but relationships among sequence architecture, basal DNA methylation and natural variation remain unclear. We combined transcription-factor motif prediction, transposable-element comparisons, methylome and RNA sequencing (RNA-seq) data, and Arabidopsis genome assemblies. In silico disruption of five HSE cores eliminated HSF-family motif compatibility in the selected design and all 5119 exact-guanine-cytosine (GC) alternatives. Across 16 curated Columbia-0 terminal windows, ONSEN contained 33-49 non-redundant HSF motif-coordinate placements per 800 bp window and was strongly enriched relative to 1930 non-ONSEN transposable elements across score thresholds and continuous metrics. Direct comparison with 779 non-ONSEN LTR retrotransposons showed selectively elevated basal CHH methylation (where H = A, C or T) at ONSEN termini. Genome-wide RNA-seq analysis revealed broad heat-responsive gene and transposable-element changes, including strong ONSEN induction, whereas candidate-window analysis distinguished ONSEN from most HSF-rich non-ONSEN outliers. ONSEN-like variants across eight accessions generally retained HSF-compatible motifs while altering predicted DNA binding with one finger-family motif composition. Together, these findings define ONSEN terminal regions as HSF-rich regulatory sequences that retain heat-responsive potential within a methylated chromatin context and identify candidates for functional analysis.

DNA Methylation↗

PubMed2026

Complete telomere-to-telomere genome assembly of Guazuma ulmifolia uncovers evolutionary mechanisms, drought adaptation, and flavonoid biosynthesis.

The first T2T reference genome of Guazuma ulmifolia is reported, which serves as a core genomic resource for stress adaptation research and stress-tolerant breeding in cacao wild relatives. Climate change, particularly increased incidence of drought, poses a major threat to food security. Understanding the genomic basis of environmental adaptation in crop wild relatives can provide valuable resources for improving stress resilience. Guazuma ulmifolia, a wild relative of Theobroma cacao with important ecological and medicinal value, lacks high-quality reference genomic resources. Here, we report the first telomere-to-telomere (T2T) chromosome-level genome assembly of G. ulmifolia, with a genome size of 311.31 Mb, contig N50 of 35.19 Mb, and 98.70% BUSCO completeness. Repetitive sequences constitute 27.43% of the G. ulmifolia genome, with LTR retrotransposons as the predominant class. Comparative genomic analyses revealed that genome-size variation among Malvaceae species is associated with differences in polyploidization history and TE dynamics. Ancestral karyotype reconstruction identified five lineage-specific chromosome fusion events distinguishing G. ulmifolia from T. cacao. Comparative analyses further identified tandem duplication-associated expansion of stress-related LEA and GST gene families, suggesting potential genomic features associated with stress responses. Flavonoid biosynthesis genes were largely conserved in copy number but showed tissue-specific expression patterns, providing candidate genes for investigating secondary metabolism. Together, this study establishes a high-quality T2T genome resource for exploring genome evolution, chromosome organization, and stress-related genomic features in Malvaceae.

Genome, Plant↗

PubMed2026

Genome assembly of Astatotilapia latifasciata uncovers B chromosome-linked chromatin reorganization.

B chromosomes (Bs) are supernumerary genomic elements found in many eukaryotes, yet their full sequence composition, functional potential, and regulatory impact on the host genome remain unclear. Here, we present a chromosome-level genome assembly of the cichlid fish Astatotilapia latifasciata, integrating PacBio long reads, Illumina short reads, and Hi-C chromatin contact maps to resolve both A and B chromosomes. The 0.93 Gb assembly (N50 = 36.2 Mb) includes a 34 Mb B chromosome containing 789 predicted protein-coding genes and a markedly higher density of transposable elements (TEs), especially long terminal repeats (LTR) retrotransposons. Transcriptome profiling revealed that B-linked genes are predominantly transcriptionally repressed relative to their A chromosome paralogs. Hi-C-based chromatin modeling uncovered distinct 3D structural configurations associated with the B chromosome, including fewer topologically associating domains (TADs), reduced loop formation, and altered compartmentalization. These changes are linked to long-range chromatin interactions and genomic rearrangements, suggesting that the B chromosome reshapes the nuclear architecture of the host genome. Our study proposes a potential regulatory role of Bs in genome and provides a genomic resource for investigating chromosome evolution in cichlids.

Animals↗

PubMed2026

Chromosome-level genome assembly of Manglietia pachyphylla.

Manglietia pachyphylla, an endangered evergreen tree within the Magnoliaceae family, is renowned for its exceptional ornamental value in landscape horticulture. Despite its classification as a Category II nationally protected plant species in China, the genetic basis of its adaptive traits and conservation priorities remains poorly understood. To address this, we present the first chromosome-scale genome assembly of M. pachyphylla utilizing an integrated approach combining PacBio HiFi long-read and Hi-C chromosome conformation capture sequencing technologies. The assembled genome spans 2.15 Gb (contig N50 = 43.57 Mb), exhibiting a heterozygosity rate of 0.78% and repeat content of 78.64%, predominantly comprising long terminal repeat (LTR) retrotransposons (52.86%). Hi-C scaffolding anchored 99.57% of the assembly to 19 pseudochromosomes, achieving a BUSCO completeness score of 96.4%. Annotation revealed 42,505 putative protein-coding genes, with 84.46% of predicted genes were functionally annotated. Phylogenomic analysis positioned M. pachyphylla and Oyama sieboldii clustered together in a well-supported group. This high-contiguity genome assembly enables future investigations into adaptive evolution, functional genomics, and evidence-based conservation strategies for this endangered species.

Chromosomes, Plant↗

PubMed2026

Chromosome-level genome assembly of the hemiparasitic Taxillus sutchuenensis (Loranthaceae).

Taxillus sutchuenensis, an ecologically and medicinally important hemiparasitic plant that parasitizes diverse woody hosts, was sequenced to generate a high-quality chromosome-level genome assembly. PacBio HiFi long reads, RNA-seq transcriptome data, and Hi-C data were used to assemble a 406.32 Mb genome anchored onto nine pseudo-chromosomes, with a scaffold N50 of 45.59 Mb. The assembly showed high completeness and accuracy, supported by BUSCO (93.6%) and Merqury QV (70.6) assessments. The LTR Assembly Index (LAI) of 13.98 indicated excellent continuity. A total of 21,795 protein-coding genes were predicted, with 94.46% functionally annotated. Repetitive sequences accounted for 50.05% of the genome, primarily LTR retrotransposons. This genome provides a valuable resource for investigating the evolution, functional genomics, and parasitic mechanisms of hemiparasitic plants.

Genome, Plant↗

PubMed2026

Genomes of Conopholis americana and Epifagus virginiana: two holoparasitic plants (Orobanchaceae).

Conopholis americana (American cancer-root) and Epifagus virginiana (beechdrops) are sister genera of holoparasitic plants (Orobanchaceae) native to eastern North America, parasitizing oaks and American beech, respectively. Both have served as models for plastid genome reduction, yet no nuclear genomes exist for either genus or any New World holoparasitic Orobanchaceae. Here we present the first nuclear genome assemblies for both species using PacBio HiFi sequencing. The C. americana assembly totals 1.82 Gb and E. virginiana totals 440 Mb, representing an approximately 4-fold difference in genome size between these sister genera. We observed a BUSCO completeness of 79% to 80% in both species, which is typical of holoparasites. While gene prediction identified 33,889 genes in C. americana and 21,031 in E. virginiana, repeat annotation revealed that LTR retrotransposons account for 78% of the genome size difference. These assemblies reveal contrasting mechanisms of genome evolution in sister holoparasitic genera and provide foundational resources for comparative genomics of parasitic plants.

Genome, Plant↗

PubMed2025

A transposable element insertion in AUX/IAA16 disrupts splicing and causes auxin resistance in Bassia scoparia.

A dicamba-resistant population of kochia (Bassia scoparia) identified in Colorado, USA in 2012 was used to generate a synthetic mapping population that segregated for dicamba resistance. Linkage mapping associating dicamba injury with genotype derived from restriction-site-associated DNA sequencing identified a single locus in the kochia genome associated with resistance on chromosome 4. A mutant version of Auxin/Indole-3-Acetic Acid 16 (AUX/IAA16; a gene previously implicated in dicamba resistance in kochia) was found near the middle of this locus in resistant plants. Long-read sequencing of dicamba-resistant plants identified a recently inserted long-terminal repeat (LTR) retrotransposon TRIM element near the beginning of the second exon of AUX/IAA16, leading to disruption of normal splicing and a mutated degron domain. Stable transgenic lines of Arabidopsis thaliana ectopically expressing the mutant and wild-type alleles of AUX/IAA16 were developed. Arabidopsis thaliana plants expressing the mutant AUX/IAA16 allele grew shorter roots on control media. However, transgenic root growth was less inhibited on media containing either dicamba (5 μM) or IAA (0.5 μM) when compared with non-transgenic plants or those expressing the wild-type allele of AUX/IAA16. In vitro assays indicate reduced binding affinity and more rapid dissociation of the mutant AUX/IAA16 with TIR1 in the presence of several auxins, and protein modeling suggests the substitution of the glycine residue in the degron domain of AUX/IAA16 is especially important for resistance. A fitness cost associated with the mutant allele of AUX/IAA16 has implications for resistance evolution and management of kochia populations with this resistance mechanism.

Indoleacetic Acids↗

PubMed2025

Transposable elements drive evolution and perturb gene expression in Brassica rapa and B. oleracea.

Transposable elements (TEs) significantly influence genomic diversity and gene regulation in plants. Brassica rapa and B. oleracea, with their distinct domestication histories, offer excellent models to explore TE dynamics. Here, we developed a refined TE classification method and systematically analyzed TEs across 12 B. rapa and B. oleracea genomes, identifying 1878 TE families. Approximately half (49.5%) of these TE families were shared between the two species, reflecting a common evolutionary origin, whereas species-specific expansions, particularly among long-terminal repeat (LTR) retrotransposons, underscore their roles in genomic differentiation. We notably characterized a heat-responsive Ty1-copia family (Copia0035) in B. oleracea roots, distinguished by low GC content and the absence of CG and CHG methylation motifs, sharing regulatory similarities with the Arabidopsis heat-induced ONSEN element. Syntenic analyses of gene-TE associations highlighted significant intraspecies TE insertion variability, with more accession-specific insertions in B. rapa and more conserved insertions, often associated with distinct morphotypes in B. oleracea. Gene ontology enrichment indicated TE involvement in developmental, reproductive, and stress response pathways. Transcriptome analysis across diverse accessions revealed that genes proximal to TEs, particularly those regulating floral development and flowering time, exhibit increased expression variability. These findings advance our understanding of TE-mediated genome evolution in Brassica species and underscore their potential utility in breeding and genome engineering strategies for crop improvement.

DNA Transposable Elements↗

PubMed2026

The complete sequence of the silkworm W chromosome uncovers its rapid evolution by large-scale duplications/deletions and translocation of W-linked genes.

The complete sequence of the W chromosome, which carries feminization activity in the silkworm, is crucial for understanding the sex-determination system in Lepidoptera. However, extensive accumulation of transposons due to lack of recombination, the very rare protein-coding genes and almost no information about molecular markers has hindered full W sequencing. We report the first complete silkworm W sequence (T2T_W, 11683305 bp) obtained by combining sequencing-assembly technologies and newly developed error detection methods, evaluated with genetically mapped W-RAPD markers, W-mutants, and W-derived BAC clones. The T2T_W sequence showed that the W is composed of a massive 92% accumulation of transposons and repeat sequences, among which the main constituents are intact LTR/LINE retrotransposons indicating recent expansions. In addition to Fem clusters producing Fem piRNA (Feminizer-derived PIWI-interacting RNA), we found 26 protein-coding genes in the W sequence. These include four gene pairs encoding zinc-finger motifs designated z1:z20 and a gene encoding serine/arginine repetitive matrix protein 1-like (SRRM1-like). To identify candidate genes for female sex-determination and differentiation we also sequenced the shortest W (3.8 Mb) from a translocation mutant with feminizing activity, which harbored four conventional genes: a Fem cluster, a pair of z1:z20 isoforms, z20-S, and a SRRM1-like gene. Phylogenetic analysis revealed that z1:z20 originated from a copy of an autosomal zinc-finger gene pair, z2:z21, translocated onto the W around 2.43 Mya and subsequently amplified to yield 4 W-linked zinc-finger gene pairs. The complete W sequence revealed that large-scale deletions and amplifications played a significant role in W chromosome evolution.

Animals↗