SEARCH · Search PubMed

Results for “Retrotransposon”

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.

At least 37 records · Page 2Linked to original sources

PubMed2025

Distinct evolutionary trajectories of subgenomic centromeres in polyploid wheat.

BACKGROUND: Centromeres are crucial for precise chromosome segregation and maintaining genome stability during cell division. However, their evolutionary dynamics, particularly in polyploid organisms with complex genomic architectures, remain largely enigmatic. Allopolyploid wheat, with its well-defined hierarchical ploidy series and recent polyploidization history, serves as an excellent model to explore centromere evolution. RESULTS: In this study, we perform a systematic comparative analysis of centromeres in common wheat and its corresponding ancestral species, utilizing the latest comprehensive reference genome assembly available. Our findings reveal that wheat centromeres predominantly consist of five types of centromeric-specific retrotransposon elements (CRWs), with CRW1 and CRW2 being the most prevalent. We identify distinct evolutionary trajectories in the functional centromeres of each subgenome, characterized by variations in copy number, insertion age, and CRW composition. By utilizing CENH3-ChIP data across various ploidy levels, we uncover a series of CRW invasion events that have shaped the evolution of AA subgenome centromeres. Conversely, the evolutionary process of the DD subgenome centromeres involves their expansion from diploid to hexaploid wheat, facilitating adaptation to a larger genomic context. Integration of complete einkorn centromere assemblies and Aegilops tauschii pan-genomes further revealed subgenome-specific centromere evolutionary trajectories. By inclusion of synthetic hexaploid from S2-S3 generations, alongside 2x/6 × natural accessions, we demonstrate that DD subgenome centromere expansion represents a gradual evolutionary process rather than an immediate response to polyploidization. CONCLUSIONS: Our study provides a comprehensive landscape of centromere adaptation, evolution, and maturation, along with insights into how retrotransposon invasions drive centromere evolution in polyploid wheat.

Centromere↗

PubMed2024

Structure and Methyl-lysine Binding Selectivity of the HUSH Complex Subunit MPP8.

The Human Silencing Hub (HUSH) guards the genome from the pathogenic effects of retroelement expression. Composed of MPP8, TASOR, and Periphilin-1, HUSH recognizes actively transcribed retrotransposed sequences by the presence of long (>1.5-kb) nascent transcripts without introns. HUSH recruits effectors that alter chromatin structure, degrade transcripts, and deposit transcriptionally repressive epigenetic marks. Here, we report the crystal structure of the C-terminal domain (CTD) of MPP8 necessary for HUSH activity. The MPP8 CTD consists of five ankyrin repeats followed by a domain with structural homology to the PINIT domains of Siz/PIAS-family SUMO E3 ligases. AlphaFold3 modeling of the MPP8-TASOR complex predicts that a SPOC domain and a domain with a novel fold in TASOR form extended interaction interfaces with the MPP8 CTD. Point mutations at these interfaces resulted in loss of HUSH-dependent transcriptional repression in a cell-based reporter assay, validating the AlphaFold3 model. The MPP8 chromodomain, known to bind the repressive mark H3K9me3, bound with similar or higher affinity to sequences in the H3K9 methyltransferase subunits SETDB1, ATF7IP, G9a, and GLP. Hence, MPP8 promotes heterochromatinization by recruiting H3K9 methyltransferases. Our work identifies novel structural elements in MPP8 required for HUSH complex assembly and silencing, thereby fulfilling vital functions in controlling retrotransposons.

Humans↗

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↗

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↗

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↗

PubMed2025

Genomic insights into karyotype evolution and adaptive mechanisms in Polygonaceae species.

Polygonaceae, with ecological versatility and global distribution, is an ideal system for investigating plant adaptation. However, the genomic mechanisms underlying its karyotype evolution and environmental resilience remain unclear. We herein present chromosome-level genomes of 11 species from 10 Polygonaceae genera. Our analyses reveal that Gypsy retrotransposons are key drivers of genome size variations in Polygonaceae. We reconstructed a Polygonaceae ancestral karyotype comprising 28 proto-chromosomes and elucidated evolutionary trajectories via extensive chromosomal rearrangements. Furthermore, we constructed a cross-genus super pan-genome for Polygonaceae, identifying 80,055 gene families, of which 9,845 (12.30%) are core gene families. Private genes are found to contribute significantly to interspecific differences in adaptability. Notably, gene copy number variations are identified as a critical factor influencing adaptations to diverse niches involving species-specific increases in metabolic pathways. This study provides a genomic framework for Polygonaceae karyotype plasticity and adaptive innovation, offering insights into plant evolution under environmental challenges.

Karyotype↗

PubMed2026

Epigenetic drift and LINE-1 activation in aging brain: Implications for neurodegenerative disease.

Brain aging and age-associated neurological diseases, such as Alzheimer's Disease (AD), Parkinson's Disease (PD), and Amyotrophic Lateral Sclerosis (ALS), are largely attributed to epigenetic drift which is characterized by the gradual accumulation of alterations in neural cell methylation patterns over time. These methylation changes are particularly evident in transposable element (TE)-derived sequences such as Long interspersed element-1 (LINE-1) which comprises approximately 17% of the human genome. During aging, LINE-1 elements gradually lose their methylation, as well as the regulatory safeguard mechanisms that usually keep them inactive. This repression loss can lead to LINE-1 reactivation, contributing to harmful effects including genomic instability, neuroinflammation, and more. Together these findings indicate that impaired epigenetic maintenance, especially in repetitive genome regions, plays a key role in biological aging of neurons and glial cells. In this narrative review, we discuss the methylation dynamics and regulatory mechanisms of LINE-1 retrotransposons, their activation processes during aging, and contribution to age-associated neurological diseases. We also highlight the potential of targeting LINE-1 methylation to restore methylation homeostasis, epigenetic stability and delay brain aging.

Humans↗

PubMed2026

THE1 repeats: Ancient endogenous retroviruses rampaging behind sarcoid myopathy.

Transposon-like human element (THE1) repeats are classified as part of the ERVL (endogenous retrovirus-like) -MaLR family and were previously thought to be remnants derived from retroviral infections in the germ cells of ancient primates. At present, at least 60,561 THE1 repeats have been mapped in the human genome, and they are divided into four subfamilies: THE1A, THE1B, THE1C, and THE1D. Although these repeats are largely silenced in somatic cells through multiple mechanisms, some of them are known to influence gene functions and genomic structures, particularly under disease conditions. Recently, we happened to identify the CIR1 driven from THE1B in the individuals with sarcoid myopathy (SM), which is characterized by the non-caseating granulomatous inflammation recognized in skeletal muscles biopsy. Retrotrans-genomics based on the COFFEE method showed that various THE1B fusion transcripts (THE1B/FTs) and retrotransposon-associated genomic dynamics in muscle of SM, including (1) splicing of THE1B to neighborhood exons, (2) expression of THE1B associated long noncoding RNAs in the intergenic regions, (3) ectopic expression of tissue specific transcripts driven by THE1B in SM, and (4) identification of a read through transcript, SIRPB1-SIRPD. Seven of eight THE1B/FTs downregulated by tofacitinib in cutaneous sarcoidosis were confined to the granuloma associated macrophage in sarcoidosis. The remaining transcript, SMFT4 exhibited stronger expression in mesenchymal stem cells. In addition, ectopic co-expression of canonical and LINE1 and THE1A-exonizing IL23R transcripts were identified in SM. This directory of THE1B/FTs in SM may serve as a fundamental molecular basis for understanding the pathogenesis of sarcoidosis.

Journal Article↗

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-scale genomes and population resequencing resolve subgenome diversity and halophyte adaptation in Salicornia.

Amid escalating water scarcity and groundwater depletion, halophytes such as Salicornia (Amaranthaceae) represent valuable models for extreme salt tolerance and hold promise for saltwater-based agriculture. Here, we show chromosome-scale genome assemblies for six Salicornia species, revealing four distinct subgenomes, reconciling our assemblies with two existing reference genomes (S. ramosissima UK and S. europaea China), correcting chromosome numbering and orientation. Comparative analyses across ploidy levels demonstrate genome expansion in North American lineages driven by Gypsy retrotransposons, and lineage-specific expansions of two gene families implicated in stress metabolism. Phylogenetic and population-structure analyses of a global resequencing panel of 318 accessions resolve interspecific relationships and establish curated germplasm collections for future crop breeding. Genetic analyses uncover a contrasting population-genetic signal on chromosome 6A between two species, highlighting an OSCA calcium-permeable channel gene as a candidate locus for osmotic adaptation. Together, these resources establish a genomic framework for Salicornia that supports evolutionary studies of halophyte adaptation and crop development.

Chenopodiaceae↗

PubMed2025

A chromosome-level genome assembly of Coffea arabica L. var. 'Kona Typica'.

Coffea arabica L. var. 'Kona Typica' is renowned for its premium cup quality, but its vulnerability to pests and diseases limits production. To accelerate cultivar improvement, we generated a chromosome-level genome assembly of 'Kona Typica' using PacBio HiFi sequencing and Hi-C scaffolding technology. The final assembly spans 1.13 Gb, with a scaffold N50 of 50.50 Mb, organized into 22 chromosomes. BUSCO assessment indicated a high completeness at 99.1%. We annotated 65,458 protein-coding genes and identified 1,073,545 interspersed repeats, accounting for 65.16% of the genome. Analysis of transposon insertion ages revealed that most long terminal repeat retrotransposons proliferated after the polyploidization event. This high-quality genome assembly of 'Kona Typica' provides a valuable resource for exploring coffee genomic evolution and genetic mechanisms of complex traits, facilitating genomics studies and the development of improved coffee cultivars with enhanced disease resistance and quality traits.

Coffea↗

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

Cleavage region organizes the structural architecture of the SINE-derived B2 repressive ribozyme.

The SINE-encoded B2 retrotransposon is an RNA Polymerase III (POL-III)-derived transcript whose expression is substantially upregulated during various cellular stress responses. Beyond retrotransposition, the B2 non-coding RNA can directly bind and repress the activity of RNA Polymerase II (POL-II), leading to a significant downregulation of transcripts during stress. Notably, our recent findings have shown that B2 is a self-cleaving ribozyme whose activity can be induced by interactions with chromatin-modifying factors through non-canonical epigenetic mechanisms that co-regulate its function across distinct chromatin-binding target loci. Here, by integrating RNA chemical probing, small-angle X-ray scattering, and 3D motif modeling, we determine structural ensemble-to-function relations for the B2 SINE ribozyme RNA. Genetic perturbations of the RNA suggest that the B2 SINE ribozyme has a well-defined secondary and dynamic tertiary structure that depends on the integrity of the critical region, which confers ribozymatic activity and repressive extent by POL-II. Using an RNA engineering approach, we examine the effects of point mutations, deletions of the main cleavage site, and deletions of the cleavage domain on the structural ensemble of the RNA. Combining this approach with in vitro and in vivo functional perturbation methods highlights the relationships between structural ensembles and various biologically relevant functional outcomes.

RNA, Catalytic↗

PubMed2025

STING inhibits LINE-1 retrotransposition through sorting ORF1p to lysosomes for degradation.

The cyclic dinucleotide sensor stimulator of interferon (IFN) genes (STING) is known for its critical role in interferon and inflammatory responses. In addition, STING also has functions independent of interferon induction. In this study, we report that STING restricts the mobilization of the cellular retrotransposon long interspersed nuclear element 1 (LINE-1) independent of cGAS and interferon induction. LINE-1 is the only active autonomous retrotransposable element in the human genome and its transposition can cause genetic and autoimmune diseases. STING inhibition of LINE-1 requires its dimerization. Mechanistically, STING interacts with LINE-1 ORF1p, then the complex translocates to the ER-Golgi intermediate compartment (ERGIC) and the Golgi followed by sorting to Rab7-positive lysosomes for degradation. Our data unveil a function of STING in maintaining host genome integrity by restricting LINE-1 retrotransposition via an IFN-independent mechanism.

Humans↗

PubMed2026

Transcriptional perturbation of LINE-1 elements reveals their cis-regulatory potential.

Long interspersed element-1 (LINE-1 or L1) retrotransposons constitute the largest transposable element family in mammalian genomes and contribute prominently to inter- and intra-individual genetic variation. Although most L1 elements are inactive, some evolutionary younger elements remain intact and genetically competent for transcription and occasionally retrotransposition. Despite being generally more abundant in gene-poor regions, intact or full-length L1s (FL-L1) are also enriched around specific classes of genes and on the eutherian X chromosome. How proximal FL-L1 may affect nearby gene expression remains unclear. Here, we examine this systematically using engineered mouse embryonic stem cells (ESCs) in which expression of one active L1 subfamily is perturbed. We find that FL-L1 activation leads to the misregulation of ~1024 genes, whereas FL-L1 repression affects ~81 genes. In most cases (68%), misexpressed genes contain an intronic FL-L1 or lie near a FL-L1 (< 260 kb). Gene ontology analysis shows that upon L1 activation, upregulated genes are enriched for neuronal function-related terms, suggesting that some L1 elements may have evolved to control neuronal gene networks. These results illustrate the cis-regulatory potential of FL-L1 elements and suggest a broader role for L1s than originally anticipated.

Animals↗

PubMed2025

Correlation of extracellular vesicle Alu RNA with brain aging and neuronal injury: a potential biomarker for brain aging.

BACKGROUND: Extracellular vesicles (EVs) are promising biomarkers for neurodegeneration. Alu elements are retrotransposons increasingly expressed with age and may be involved in aging-related diseases. OBJECTIVE: To determine the potential of Alu RNA in plasma-derived EVs as a biomarker for brain aging and neuronal injury. METHODS: EVs were isolated from plasma samples across different age groups. EV Alu RNA levels were measured and their associations with biomarkers of brain aging, including plasma neurofilament light chain (NfL), plasma amyloid-beta (Aβ42 and Aβ40), and plasma phosphorylated tau (p-Tau181), were analyzed. RESULTS: EV Alu RNA levels were increased significantly with age and were strongly correlated with plasma NfL, suggesting a strong association between EV Alu RNA and neuronal injury. Significant correlations were also found between EV Alu RNA and plasma amyloid-beta levels, while no significant association was observed with tau pathology. CONCLUSIONS: EV Alu RNA levels are elevated with age and associated with neuronal injury, highlighting their potential as a novel, non-invasive biomarker for brain aging and neurodegeneration.

Humans↗

PubMed2026

Molecular investigation of the progenitors, origin and domestication patterns of diploid Chinese old garden roses.

BACKGROUND AND AIMS: Chinese old garden roses are major contributors to the genetic development of modern roses. The RoKSN gene is associated with continuous flowering in roses and is proposed to have originated from Chinese wild roses. However, the wild roses that are implicated in the breeding of Chinese old garden roses and the origin of the RoKSN locus remain unidentified. We collected 25 of the most renowned and classic diploid Chinese old garden roses along with all related wild roses from East Asia. These roses were analysed with the aim of identifying the wild species that contributed to the genetic composition of Chinese old garden roses. In addition, we aimed to infer the geographical origin of the RoKSN gene and to develop a schematic overview of hybrid domestication of Chinese old garden roses. METHODS: We compared the haplotypes of internal transcribed spacers (nrITS), six nuclear single-copy genes and three chloroplast genes between Chinese old garden roses and wild roses. Additionally, we assessed genetic organization using 21 expressed sequence tag-simple sequence repeats to identify potential donor species that contributed to the emergence of these cultivars. Primers were designed for RoKSN to allow comparison of the gene across the entire distribution range of Rosa sect. Chinenses. KEY RESULTS: Our findings confirmed that the majority of rose cultivars are descendants of early hybridization events. Rosa chinensis var. spontanea, R. odorata var. gigantea and R. multiflora var. cathayensis were the primary donors for the 25 cultivar roses. Chinese old garden roses were categorized into four groups. Ten cultivars were hybrids between R. chinensis var. spontanea and R. multiflora var. cathayensis, thereby forming the 'Old Blush' group. Five cultivars were hybrids between 'Old Blush' and the R. kwangtungensis species complex, thereby forming the 'Slater's crimson' group. Six cultivars were hybrids between 'Old Blush' and R. odorata var. gigantea, thereby forming the 'Tea Rose' group, and three cultivars were hybrids that evolved from more than three donors. Moreover, we observed relatively close genetic proximity among Chinese old garden roses with an identical RoKSN-copia gene that is responsible for continuous flowering, which indicates a single origin for this retrotransposon-containing allele. Additionally, we determined that the haplotypes of the RoKSN-copia gene predominantly occurred in the Sichuan Basin region. In contrast, R. chinensis cultivated in the Ya'an region showed no markers of hybridization and displayed a genetic composition that was close to that of the wild species R. chinensis var. spontanea. This cultivar may represent the earliest mutated individual that bears the RoKSN-copia gene and may have served as a bridge from wild species to continuous-flowering old rose cultivars. CONCLUSIONS: The study provides crucial evidence that elucidates the origin of cultivated roses and lays the groundwork for further analysis of the breeding history of Chinese old garden roses using genomic data.

Domestication↗