Exonic Short Interspersed Nuclear Element Insertion in FAM161A Is Associated with Autosomal Recessive Progressive Retinal Atrophy in the English Shepherd.

Progressive retinal atrophies (PRAs) are a genetically heterogeneous group of inherited eye diseases that affect over 100 breeds of dog. The initial clinical sign is visual impairment in scotopic conditions, as a consequence of rod photoreceptor cell degeneration. Photopic vision degeneration then follows, due to progression of the disease to the cone photoreceptors, and ultimately results in complete blindness. Two full-sibling English Shepherds were diagnosed with PRA at approximately 5 years old and tested clear of all published PRA genetic variants. This study sought to identify the novel PRA-associated variant segregating in the breed. We utilised a combined approach of whole genome sequencing of the probands and homozygosity mapping of four cases and 22 controls and identified a short interspersed nuclear element within an alternatively spliced exon in FAM161A. The XP_005626197.1 c.17929_ins210 variant was homozygous in six PRA cases and heterozygous or absent in control dogs, consistent with a recessive mode of inheritance. The insertion is predicted to extend exon 4 by 39 aberrant amino acids followed by an early termination stop codon. PRA is intractable to treatment, so the development of a genetic screening test, based on the associated variant, is significant, because it provides dog breeders/owners with a means of reducing the frequency of the disease variant within this breed as well as minimising the risk of breeding puppies that will develop this blinding disease.

Mapping of SINEs in the genome of Proechimys (Mammalia: Rodentia).

This study aimed to analyze the distribution of short interspersed elements (SINEs) in the chromosomes of five species of rodents of the genus Proechimys and in a variant karyotype of P. guyannensis. Molecular cytogenetic techniques were used to characterize the sequences of the B1, B4, MAR and THER SINEs, which were used as probes for hybridization in metaphase chromosomes. A wide distribution of SINEs was observed in the chromosomes of the Proechimys species examined, thus indicating differentiation of these retroelements. The signal of the B4 SINE was more evident than that of the B1 SINE, especially in P. echinothrix, P. longicaudatus, and P. cuvieri. Although the signal of the MAR SINE was more explosive than that of the THER SINE, in the species P. echinothrix, P. guyannensis (2n = 46) and P. longicaudatus, its distribution in the karyotypes was similar. The signals of these retroelements occurred at specific heterochromatic sites and were centromeric/pericentromeric and at the terminal regions in most chromosomes. This appears to be a typical distribution pattern of the SINEs and may indicate involvement with rearrangements during karyotypic diversification in Proechimys. The variation of the SINEs in the genome of Proechimys species demonstrates that these elements are distributed in a specific way in this genus and the preference for some sites, considered hotspots for chromosomal breakage, allows us to propose that these elements are related to the karyotypic evolution of Proechimys.

Variable patterns of retrotransposition in different HeLa strains provide mechanistic insights into SINE RNA mobilization processes.

Alu elements are non-autonomous Short INterspersed Elements (SINEs) derived from the 7SL RNA gene that are present at over one million copies in human genomic DNA. Alu mobilizes by a mechanism known as retrotransposition, which requires the Long INterspersed Element-1 (LINE-1) ORF2-encoded protein (ORF2p). Here, we demonstrate that HeLa strains differ in their capacity to support Alu retrotransposition. Human Alu elements retrotranspose efficiently in HeLa-HA and HeLa-CCL2 (Alu-permissive) strains, but not in HeLa-JVM or HeLa-H1 (Alu-nonpermissive) strains. A similar pattern of retrotransposition was observed for other 7SL RNA-derived SINEs and tRNA-derived SINEs. In contrast, mammalian LINE-1s, a zebrafish LINE, a human SINE-VNTR-Alu (SVA) element, and an L1 ORF1-containing mRNA can retrotranspose in all four HeLa strains. Using an in vitro reverse transcriptase-based assay, we show that Alu RNAs associate with ORF2p and are converted into cDNAs in both Alu-permissive and Alu-nonpermissive HeLa strains, suggesting that 7SL- and tRNA-derived SINEs use strategies to 'hijack' L1 ORF2p that are distinct from those used by SVA elements and ORF1-containing mRNAs. These data further suggest ORF2p associates with the Alu RNA poly(A) tract in both Alu-permissive and Alu-nonpermissive HeLa strains, but that Alu retrotransposition is blocked after this critical step in Alu-nonpermissive HeLa strains.

Deciphering the role of a SINE-VNTR-Alu retrotransposon polymorphism as a biomarker of Parkinson's disease progression.

SINE-VNTR-Alu (SVA) retrotransposons are transposable elements which represent a source of genetic variation. We previously demonstrated that the presence/absence of a human-specific SVA, termed SVA_67, correlated with the progression of Parkinson's disease (PD). In the present study, we demonstrate that SVA_67 acts as expression quantitative trait loci, thereby exhibiting a strong regulatory effect across the genome using whole genome and transcriptomic data from the Parkinson's progression markers initiative cohort. We further show that SVA_67 is polymorphic for its variable number tandem repeat domain which correlates with both regulatory properties in a luciferase reporter gene assay in vitro and differential expression of multiple genes in vivo. Additionally, this variation's utility as a biomarker is reflected in a correlation with a number of PD progression markers. These experiments highlight the plethora of transcriptomic and phenotypic changes associated with SVA_67 polymorphism which should be considered when investigating the missing heritability of neurodegenerative diseases.

High-Copy Transposons from a Pathogen Give Rise to a Conserved sRNA Family with a Novel Host Immunity Target.

Small RNAs (sRNAs) are involved in gene silencing in multiple ways, including through cross-kingdom transfers from parasites to their hosts. Little is known about the evolutionary mechanisms enabling eukaryotic microbes to evolve functional mimics of host small regulatory RNAs. Here, we describe the identification and functional characterization of SINE_sRNA1, an sRNA family derived from highly abundant short interspersed nuclear element (SINE) retrotransposons in the genome of the wheat powdery mildew pathogen. SINE_sRNA1 is encoded by a sequence motif that is conserved in multiple SINE families and corresponds to a functional plant microRNA (miRNA) mimic targeting Tae_AP1, a wheat gene encoding an aspartic protease only found in monocots. Tae_AP1 has a novel function enhancing both pattern-triggered immunity (PTI) and effector-triggered immunity (ETI), thereby contributing to the cross activation of plant defenses. We conclude that SINE_sRNA1 and Tae_AP1 are functional innovations, suggesting the contribution of transposons to the evolutionary arms race between a parasite and its host. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Decryption of sequence, structure, and functional features of SINE repeat elements in SINEUP non-coding RNA-mediated post-transcriptional gene regulation.

RNA structure folding largely influences RNA regulation by providing flexibility and functional diversity. In silico and in vitro analyses are limited in their ability to capture the intricate relationships between dynamic RNA structure and RNA functional diversity present in the cell. Here, we investigate sequence, structure and functional features of mouse and human SINE-transcribed retrotransposons embedded in SINEUPs long non-coding RNAs, which positively regulate target gene expression post-transcriptionally. In-cell secondary structure probing reveals that functional SINEs-derived RNAs contain conserved short structure motifs essential for SINEUP-induced translation enhancement. We show that SINE RNA structure dynamically changes between the nucleus and cytoplasm and is associated with compartment-specific binding to RBP and related functions. Moreover, RNA-RNA interaction analysis shows that the SINE-derived RNAs interact directly with ribosomal RNAs, suggesting a mechanism of translation regulation. We further predict the architecture of 18 SINE RNAs in three dimensions guided by experimental secondary structure data. Overall, we demonstrate that the conservation of short key features involved in interactions with RBPs and ribosomal RNA drives the convergent function of evolutionarily distant SINE-transcribed RNAs.

SINEs as Potential Expression Cassettes: Impact of Deletions and Insertions on Polyadenylation and Lifetime of B2 and Ves SINE Transcripts Generated by RNA Polymerase III.

Short Interspersed Elements (SINEs) are common in the genomes of most multicellular organisms. They are transcribed by RNA polymerase III from an internal promoter comprising boxes A and B. As transcripts of certain SINEs from mammalian genomes can be polyadenylated, such transcripts should contain the AATAAA sequence as well as those called ß- and τ-signals. One of the goals of this work was to evaluate how autonomous and independent other SINE parts are ß- and τ-signals. Extended regions outside of ß- and τ-signals were deleted from SINEs B2 and Ves and the derived constructs were used to transfect HeLa cells in order to evaluate the relative levels of their transcripts as well as their polyadenylation efficiency. If the deleted regions affected boxes A and B, the 5'-flanking region of the U6 RNA gene with the external promoter was inserted upstream. Such substitution of the internal promoter in B2 completely restored its transcription. Almost all tested deletions/substitutions did not reduce the polyadenylation capacity of the transcripts, indicating a weak dependence of the function of ß- and τ-signals on the neighboring sequences. A similar analysis of B2 and Ves constructs containing a 55-bp foreign sequence inserted between ß- and τ-signals showed an equal polyadenylation efficiency of their transcripts compared to those of constructs without the insertion. The acquired poly(A)-tails significantly increased the lifetime and thus the cellular level of such transcripts. The data obtained highlight the potential of B2 and Ves SINEs as cassettes for the expression of relatively short sequences for various applications.

Continuous synthesis of E. coli genome sections and Mb-scale human DNA assembly.

Whole-genome synthesis provides a powerful approach for understanding and expanding organism function1-3. To build large genomes rapidly, scalably and in parallel, we need (1) methods for assembling megabases of DNA from shorter precursors and (2) strategies for rapidly and scalably replacing the genomic DNA of organisms with synthetic DNA. Here we develop bacterial artificial chromosome (BAC) stepwise insertion synthesis (BASIS)-a method for megabase-scale assembly of DNA in Escherichia coli episomes. We used BASIS to assemble 1.1 Mb of human DNA containing numerous exons, introns, repetitive sequences, G-quadruplexes, and long and short interspersed nuclear elements (LINEs and SINEs). BASIS provides a powerful platform for building synthetic genomes for diverse organisms. We also developed continuous genome synthesis (CGS)-a method for continuously replacing sequential 100 kb stretches of the E. coli genome with synthetic DNA; CGS minimizes crossovers1,4 between the synthetic DNA and the genome such that the output for each 100 kb replacement provides, without sequencing, the input for the next 100 kb replacement. Using CGS, we synthesized a 0.5 Mb section of the E. coli genome-a key intermediate in its total synthesis1-from five episomes in 10 days. By parallelizing CGS and combining it with rapid oligonucleotide synthesis and episome assembly5,6, along with rapid methods for compiling a single genome from strains bearing distinct synthetic genome sections1,7,8, we anticipate that it will be possible to synthesize entire E. coli genomes from functional designs in less than 2 months.

Identification and functional characterization of the German cockroach, Blattella germanica, short interspersed nuclear elements.

In recent decades, experimental data has accumulated indicating that short interspersed nuclear elements (SINEs) can play a significant functional role in the regulation of gene expression in the host genome. In addition, molecular markers based on SINE insertion polymorphisms have been developed and are widely used for genetic differentiation of populations of eukaryotic organisms. Using routine bioinformatics analysis and publicly available genomic DNA and small RNA-seq data, we first described nine SINEs in the genome of the German cockroach, Blattella germanica. All described SINEs have tRNA promoters, and the start of their transcription begins 11 bp upstream of an "A" box of these promoters. The number of copies of the described SINEs in the B. germanica genome ranges from several copies to more than a thousand copies in a SINE-specific manner. Some of the described SINEs and their degenerate copies can be localized both in the introns of genes and loci known as piRNA clusters. piRNAs originating from piRNA clusters are shown to be mapped to seven of the nine types of SINEs described, including copies of SINEs localized in gene introns. We speculate that SINEs, localized in the introns of certain genes, may regulate the level of expression of these genes by a PIWI-related molecular mechanism.

Transcriptional Alterations in X-Linked Dystonia-Parkinsonism Caused by the SVA Retrotransposon.

X-linked dystonia-parkinsonism (XDP) is a severe neurodegenerative disorder that manifests as adult-onset dystonia combined with parkinsonism. A SINE-VNTR-Alu (SVA) retrotransposon inserted in an intron of the TAF1 gene reduces its expression and alters splicing in XDP patient-derived cells. As a consequence, increased levels of the TAF1 intron retention transcript TAF1-32i can be found in XDP cells as compared to healthy controls. Here, we investigate the sequence of the deep intronic region included in this transcript and show that it is also present in cells from healthy individuals, albeit in lower amounts than in XDP cells, and that it undergoes degradation by nonsense-mediated mRNA decay. Furthermore, we investigate epigenetic marks (e.g., DNA methylation and histone modifications) present in this intronic region and the spanning sequence. Finally, we show that the SVA evinces regulatory potential, as demonstrated by its ability to repress the TAF1 promoter in vitro. Our results enable a better understanding of the disease mechanisms underlying XDP and transcriptional alterations caused by SVA retrotransposons.

The non-LTR retrotransposons of Entamoeba histolytica: genomic organization and biology.

Genome sequence analysis of Entamoeba species revealed various classes of transposable elements. While E. histolytica and E. dispar are rich in non-long terminal repeat (LTR) retrotransposons, E. invadens contains predominantly DNA transposons. Non-LTR retrotransposons of E. histolytica constitute three families of long interspersed nuclear elements (LINEs), and their short, nonautonomous partners, SINEs. They occupy ~ 11% of the genome. The EhLINE1/EhSINE1 family is the most abundant and best studied. EhLINE1 is 4.8 kb, with two ORFs that encode functions needed for retrotransposition. ORF1 codes for the nucleic acid-binding protein, and ORF2 has domains for reverse transcriptase (RT) and endonuclease (EN). Most copies of EhLINEs lack complete ORFs. ORF1p is expressed constitutively, but ORF2p is not detected. Retrotransposition could be demonstrated upon ectopic over expression of ORF2p, showing that retrotransposition machinery is functional. The newly retrotransposed sequences showed a high degree of recombination. In transcriptomic analysis, RNA-Seq reads were mapped to individual EhLINE1 copies. Although full-length copies were transcribed, no full-length 4.8 kb transcripts were seen. Rather, sense transcripts mapped to ORF1, RT and EN domains. Intriguingly, there was strong antisense transcription almost exclusively from the RT domain. These unique features of EhLINE1 could serve to attenuate retrotransposition in E. histolytica.

Human Recombinant DNase I (Pulmozyme®) Inhibits Lung Metastases in Murine Metastatic B16 Melanoma Model That Correlates with Restoration of the DNase Activity and the Decrease SINE/LINE and c-Myc Fragments in Blood Cell-Free DNA.

Tumor-associated cell-free DNAs (cfDNA) play an important role in the promotion of metastases. Previous studies proved the high antimetastatic potential of bovine pancreatic DNase I and identified short interspersed nuclear elements (SINEs) and long interspersed nuclear elements (LINEs)and fragments of oncogenes in cfDNA as the main molecular targets of enzyme in the bloodstream. Here, recombinant human DNase I (commercial name Pulmozyme®), which is used for the treatment of cystic fibrosis in humans, was repurposed for the inhibition of lung metastases in the B16 melanoma model in mice. We found that Pulmozyme® strongly reduced migration and induced apoptosis of B16 cells in vitro and effectively inhibited metastases in lungs and liver in vivo. Pulmozyme® was shown to be two times more effective when administered intranasally (i.n.) than bovine DNase I, but intramuscular (i.m.) administration forced it to exhibit as high an antimetastatic activity as bovine DNase I. Both DNases administered to mice either i.m. or i.n. enhanced the DNase activity of blood serum to the level of healthy animals, significantly decreased cfDNA concentrations, efficiently degraded SINE and LINE repeats and c-Myc fragments in the bloodstream and induced apoptosis and disintegration of neutrophil extracellular traps in metastatic foci; as a result, this manifested as the inhibition of metastases spread. Thus, Pulmozyme®, which is already an approved drug, can be recommended for use in the treatment of lung metastases.

Anti-aging and anti-oxidant activities of murine short interspersed nuclear element antisense RNA.

Short interspersed nuclear elements (SINEs) play a key role in regulating gene expression, and SINE RNAs are involved in age-related diseases. We investigated the anti-aging effects of a genetically engineered murine SINE B1 antisense RNA (B1as RNA) and explored its mechanism of action in naturally senescent BALB/c (≥14 months) and moderately senscent C57BL/6N (≥9 months) mice. After tail vein injection, B1as RNA was available in the blood of mice for approximately 30 min, persisted for approximately 2-4 h in most detected tissues and persisted approximately 48 h in lungs. We found that treatment with B1as RNA improved stamina and promoted hair re-growth in aged mice. Treatment with B1as RNA also partially rescued the increase in mitochondrial DNA copy number in liver and spleen tissues observed in aged and moderately senescent mice. Finally, treatment with B1as RNA increased the activities of superoxide dismutase and glutathione peroxidase in aged and moderately senescent mice, reduced these animals' malondialdehyde and reactive oxygen species levels, and modulated the expression of several aging-associated genes, including Sirtuin 1, p21, p16Ink4a, p15Ink4b and p19Arf, and anti-oxidant genes (Sesn1 and Sesn 2). These data suggest that B1as RNA inhibits the aging process by enhancing antioxidant activity, promoting the scavenging of free radicals, and modulating the expression of aging-associated genes. This is the first report describing the anti-aging activity of SINE antisense RNA, which may serve as an effective nucleic acid drug for the treatment of age-related diseases.

A brain-specific pgc1α fusion transcript affects gene expression and behavioural outcomes in mice.

PGC1α is a transcriptional coactivator in peripheral tissues, but its function in the brain remains poorly understood. Various brain-specific Pgc1α isoforms have been reported in mice and humans, including two fusion transcripts (FTs) with non-coding repetitive sequences, but their function is unknown. The FTs initiate at a simple sequence repeat locus ∼570 Kb upstream from the reference promoter; one also includes a portion of a short interspersed nuclear element (SINE). Using publicly available genomics data, here we show that the SINE FT is the predominant form of Pgc1α in neurons. Furthermore, mutation of the SINE in mice leads to altered behavioural phenotypes and significant up-regulation of genes in the female, but not male, cerebellum. Surprisingly, these genes are largely involved in neurotransmission, having poor association with the classical mitochondrial or antioxidant programs. These data expand our knowledge on the role of Pgc1α in neuronal physiology and suggest that different isoforms may have distinct functions. They also highlight the need for further studies before modulating levels of Pgc1α in the brain for therapeutic purposes.

Functional characterization of T2D-associated SNP effects on baseline and ER stress-responsive ß cell transcriptional activation.

Genome-wide association studies (GWAS) have linked single nucleotide polymorphisms (SNPs) at >250 loci in the human genome to type 2 diabetes (T2D) risk. For each locus, identifying the functional variant(s) among multiple SNPs in high linkage disequilibrium is critical to understand molecular mechanisms underlying T2D genetic risk. Using massively parallel reporter assays (MPRA), we test the cis-regulatory effects of SNPs associated with T2D and altered in vivo islet chromatin accessibility in MIN6 ß cells under steady state and pathophysiologic endoplasmic reticulum (ER) stress conditions. We identify 1,982/6,621 (29.9%) SNP-containing elements that activate transcription in MIN6 and 879 SNP alleles that modulate MPRA activity. Multiple T2D-associated SNPs alter the activity of short interspersed nuclear element (SINE)-containing elements that are strongly induced by ER stress. We identify 220 functional variants at 104 T2D association signals, narrowing 54 signals to a single candidate SNP. Together, this study identifies elements driving ß cell steady state and ER stress-responsive transcriptional activation, nominates causal T2D SNPs, and uncovers potential roles for repetitive elements in ß cell transcriptional stress response and T2D genetics.

Analysis of SINE Families B2, Dip, and Ves with Special Reference to Polyadenylation Signals and Transcription Terminators.

Short Interspersed Elements (SINEs) are eukaryotic non-autonomous retrotransposons transcribed by RNA polymerase III (pol III). The 3'-terminus of many mammalian SINEs has a polyadenylation signal (AATAAA), pol III transcription terminator, and A-rich tail. The RNAs of such SINEs can be polyadenylated, which is unique for pol III transcripts. Here, B2 (mice and related rodents), Dip (jerboas), and Ves (vespertilionid bats) SINE families were thoroughly studied. They were divided into subfamilies reliably distinguished by relatively long indels. The age of SINE subfamilies can be estimated, which allows us to reconstruct their evolution. The youngest and most active variants of SINE subfamilies were given special attention. The shortest pol III transcription terminators are TCTTT (B2), TATTT (Ves and Dip), and the rarer TTTT. The last nucleotide of the terminator is often not transcribed; accordingly, the truncated terminator of its descendant becomes nonfunctional. The incidence of complete transcription of the TCTTT terminator is twice higher compared to TTTT and thus functional terminators are more likely preserved in daughter SINE copies. Young copies have long poly(A) tails; however, they gradually shorten in host generations. Unexpectedly, the tail shortening below A10 increases the incidence of terminator elongation by Ts thus restoring its efficiency. This process can be critical for the maintenance of SINE activity in the genome.

ADAR1 interaction with Z-RNA promotes editing of endogenous double-stranded RNA and prevents MDA5-dependent immune activation.

Loss of function of adenosine deaminase acting on double-stranded RNA (dsRNA)-1 (ADAR1) causes the severe autoinflammatory disease Aicardi-Goutières syndrome (AGS). ADAR1 converts adenosines into inosines within dsRNA. This process called A-to-I editing masks self-dsRNA from detection by the antiviral dsRNA sensor MDA5. ADAR1 binds to dsRNA in both the canonical A-form and the poorly defined Z conformation (Z-RNA). Mutations in the Z-RNA-binding Zα domain of ADAR1 are common in patients with AGS. How loss of ADAR1/Z-RNA interaction contributes to disease development is unknown. We demonstrate that abrogated binding of ADAR1 to Z-RNA leads to reduced A-to-I editing of dsRNA structures formed by base pairing of inversely oriented short interspersed nuclear elements. Preventing ADAR1 binding to Z-RNA triggers an MDA5/MAVS-mediated type I interferon response and leads to the development of lethal autoinflammation in mice. This shows that the interaction between ADAR1 and Z-RNA restricts sensing of self-dsRNA and prevents AGS development.

Whole-genome sequencing of H3K4me3 and DNA methylation in human sperm reveals regions of overlap linked to fertility and development.

The paternal environment has been linked to infertility and negative outcomes. Such effects may be transmitted via sperm through histone modifications. To date, in-depth profiling of the sperm chromatin in men has been limited. Here, we use deep sequencing to characterize the sperm profiles of histone H3 lysine 4 tri-methylation (H3K4me3) and DNA methylation in a representative reference population of 37 men. Our analysis reveals that H3K4me3 is localized throughout the genome and at genes for fertility and development. Remarkably, enrichment is also found at regions that escape epigenetic reprogramming in primordial germ cells, embryonic enhancers, and short-interspersed nuclear elements (SINEs). There is significant overlap in H3K4me3 and DNA methylation throughout the genome, suggesting a potential interplay between these marks previously reported to be mutually exclusive in sperm. Comparisons made between H3K4me3 marked regions in sperm and the embryonic transcriptome suggest an influence of paternal chromatin on embryonic gene expression.

Expression Quantitative Trait Loci (eQTLs) Associated with Retrotransposons Demonstrate their Modulatory Effect on the Transcriptome.

Transposable elements (TEs) are repetitive elements that belong to a variety of functional classes and have an important role in shaping genome evolution. Around 50% of the human genome contains TEs, and they have been termed the "dark matter" of the genome because relatively little is known about their function. While TEs have been shown to participate in aberrant gene regulation and the pathogenesis of diseases, only a few studies have explored the systemic effect of TEs on gene expression. In the present study, we analysed whole genome sequences and blood whole transcriptome data from 570 individuals within the Parkinson's Progressive Markers Initiative (PPMI) cohort to identify expression quantitative trait loci (eQTL) regulating genome-wide gene expression associated with TEs. We identified 2132 reference TEs that were polymorphic for their presence or absence in our study cohort. The presence or absence of the TE element could change the expression of the gene or gene clusters from zero to tens of thousands of copies of RNA. The main finding is that many TEs possess very strong regulatory effects, and they have the potential to modulate large genetic networks with hundreds of target genes over the genome. We illustrate the plethora of regulatory mechanisms using examples of their action at the HLA gene cluster and data showing different TEs' convergence to modulate WFS1 gene expression. In conclusion, the presence or absence of polymorphisms of TEs has an eminent genome-wide regulatory function with large effect size at the level of the whole transcriptome. The role of TEs in explaining, in part, the missing heritability for complex traits is convincing and should be considered.

Distinct Retrotransposon Evolution Profile in the Genome of Rabbit (Oryctolagus cuniculus).

Although the rabbit genome has already been annotated, it is mobilome remains largely unknown. Here, multiple pipelines were used to de novo mine and annotate the mobilome in rabbit. Four families and 19 subfamilies of LINE1s, two families and nine subfamilies of SINEs, and 12 ERV families were defined in rabbit based on sequence identity, structural organization, and phylogenetic tree. The analysis of insertion age and polymerase chain reaction suggests that a number of families are very young and may remain active, such as L1B, L1D, OcuSINEA, and OcuERV1. RepeatMasker annotation revealed a distinct transposable element landscape within the genome, with approximately two million copies of SINEs, representing the greatest proportion of the genome (19.61%), followed by LINEs (15.44%), and LTRs (4.11%), respectively, considerably different from most other mammal mobilomes except hedgehog and tree shrew, in which LINEs have the highest proportion. Furthermore, a very high rate of insertion polymorphisms (>85%) for the youngest subfamily (OcuSINEA1) was identified by polymerase chain reaction. The majority of retrotransposon insertions overlapped with protein-coding regions (>80%) and lncRNA (90%) genes. Genomic distribution bias was observed for retrotransposons, with those immediately upstream (-1 kb) and downstream (1 kb) of genes significantly depleted. Local GC content in 50-kb widows had significantly negative correlations with LINE (rs=-0.996) and LTR (rs=-0.829) insertions. The current study revealed a distinct mobilome landscape in rabbit, which will assist in the elucidation of the evolution of the genome of lagomorphs, and even other mammals.