Molecular and Chromosomal Evolution of DsTc1_5 and DsPogo_8 Elements (Tc1-Mariner) in Species of the Dichotomius Genus (Coleoptera: Scarabaeidae): Contributions to the Speciation Process.

In Dichotomius genus, transposable elements (TE) have been related to chromosome remodeling, genomic evolution, and, possibly, to the speciation process. The objective of this study was to verify the interpopulational and interspecific conservation/variation of Tc1-Mariner elements (possibly autonomous) in Dichotomius species, aiming to identify possible contributions in the speciation process of this group. The analysis was performed on four species of Dichotomius, belonging to the Selenocopris subgenus. We verified the presence of the DsPogo_8 and DsTc1_5 elements by PCR and sequencing. We also isolated and sequenced the 28S and 16S rRNA genes aiming at the phylogenetic reconstruction of the analyzed species. Chromosomal mapping of TEs DsTc1_5 and DsPogo_8 was performed by fluorescent in situ hybridization. The results revealed the presence of the elements in the different species analyzed, except for DsTc1_5 in D. (S.) geminatus. These results suggest a vertical inheritance, with the presence of these elements in the common ancestor of these species. In the analyzed species, the nucleotide similarity of DsTc1_5 was higher than that of the 28S and 16S rRNA genes, suggesting the occurrence of horizontal transfer. The phylogenetic tree indicated that the absence of DsTc1_5 in D. (S.) geminatus is related to stochastic loss of this TE. Chromosomal mapping revealed dispersed signals, with predominance in euchromatic regions and wide variation in the chromosomal localization pattern of DsTc1_5 and DsPogo_8, both interpopulational and interspecific. This variation indicates that DsTc1_5 and DsPogo_8 may have contributed to prezygotic and postzygotic isolation, thus contributing to the speciation of these species.

Concurrent transposon engineering and CRISPR/Cas9 genome editing of primary CLL-1 chimeric antigen receptor-natural killer cells.

BACKGROUND: Natural killer (NK) cell genome editing promises to enhance the innate and alloreactive anti-tumor potential of NK cell adoptive transfer. DNA transposons are versatile non-viral gene vectors now being adapted to primary NK cells, representing important tools for research and clinical product development. AIMS AND METHODS: We set out to generate donor-derived, primary chimeric antigen receptor (CAR)-NK cells by combining the TcBuster transposon system with Epstein-Barr virus-transformed lymphoblastoid feeder cell-mediated activation and expansion. RESULTS: This approach allowed for clinically relevant NK-cell expansion capability and CAR expression, which was further enhanced by immunomagnetic selection based on binding to the CAR target protein.The resulting CAR-NK cells targeting the myeloid associated antigen CLL-1 efficiently targeted CLL-1-positive AML cell lines and primary AML populations, including a population enriched for leukemia stem cells. Subsequently, concurrent delivery of CRISPR/Cas9 cargo was applied to knockout the NK cell cytokine checkpoint cytokine-inducible SH2-containing protein (CIS, product of the CISH gene), resulting in enhanced cytotoxicity and an altered NK cell phenotype. CONCLUSIONS: This report contributes a promising application of transposon engineering to donor-derived NK cells and emphasizes the importance of feeder mediated NK cell activation and expansion to current protocols.

Functional Characterization of the N-Terminal Disordered Region of the piggyBac Transposase.

The piggyBac DNA transposon is an active element initially isolated from the cabbage looper moth, but members of this superfamily are also present in most eukaryotic evolutionary lineages. The functionally important regions of the transposase are well described. There is an RNase H-like fold containing the DDD motif responsible for the catalytic DNA cleavage and joining reactions and a C-terminal cysteine-rich domain important for interaction with the transposon DNA. However, the protein also contains a ~100 amino acid long N-terminal disordered region (NTDR) whose function is currently unknown. Here we show that deletion of the NTDR significantly impairs piggyBac transposition, although the extent of decrease is strongly cell-type specific. Moreover, replacing the NTDR with scrambled but similarly disordered sequences did not rescue transposase activity, indicating the importance of sequence conservation. Cell-based transposon excision and integration assays reveal that the excision step is more severely affected by NTDR deletion. Finally, bioinformatic analyses indicated that the NTDR is specific for the piggyBac superfamily and is also present in domesticated, transposase-derived proteins incapable of catalyzing transposition. Our results indicate an essential role of the NTDR in the "fine-tuning" of transposition and its significance in the functions of piggyBac-originated co-opted genes.

Identification of an active miniature inverted-repeat transposable element mJing in rice.

Miniature inverted-repeat transposable elements (MITEs) are structurally homogeneous non-autonomous DNA transposons with high copy numbers that play important roles in genome evolution and diversification. Here, we analyzed the rice high-tillering dwarf (htd) mutant in an advanced backcross population between cultivated and wild rice, and identified an active MITE named miniature Jing (mJing). The mJing element belongs to the PIF/Harbinger superfamily. japonica rice var. Nipponbare and indica var. 93-11 harbor 72 and 79 mJing family members, respectively, have undergone multiple rounds of amplification bursts during the evolution of Asian cultivated rice (Oryza sativa L.). A heterologous transposition experiment in Arabidopsis thaliana indicated that the autonomous element Jing is likely to have provides the transposase needed for mJing mobilization. We identified 297 mJing insertion sites and their presence/absence polymorphism among 71 rice samples through targeted high-throughput sequencing. The results showed that the copy number of mJing varies dramatically among Asian cultivated rice (O. sativa), its wild ancestor (O. rufipogon), and African cultivated rice (O. glaberrima) and that some mJing insertions are subject to directional selection. These findings suggest that the amplification and removal of mJing elements have played an important role in rice genome evolution and species diversification.

Rigidity and flexibility characteristics of DD[E/D]-transposases Mos1 and Sleeping Beauty.

DD[E/D]-transposases catalyze the multistep reaction of cut-and-paste DNA transposition. Structurally, several DD[E/D]-transposases have been characterized, revealing a multi-domain structure with the catalytic domain possessing the RNase H-like structural motif that brings three catalytic residues (D, D, and E or D) into close proximity for the catalysis. However, the dynamic behavior of DD[E/D]-transposases during transposition remains poorly understood. Here, we analyze the rigidity and flexibility characteristics of two representative DD[E/D]-transposases Mos1 and Sleeping Beauty (SB) using the minimal distance constraint model (mDCM). We find that the catalytic domain of both transposases is globally rigid, with the notable exception of the clamp loop being flexible in the DNA-unbound form. Within this globally rigid structure, the central ß-sheet of the RNase H-like motif is much less rigid in comparison to its surrounding α-helices, forming a cage-like structure. The comparison of the original SB transposase to its hyperactive version SB100X reveals the region where the change in flexibility/rigidity correlates with increased activity. This region is found to be within the RNase H-like structural motif and comprise the loop leading from beta-strand B3 to helix H1, helices H1 and H2, which are located on the same side of the central beta-sheet, and the loop between helix H3 and beta-strand B5. We further identify the RKEN214-217DAVQ mutations of the set of hyperactive mutations within the catalytic domain of SB transposase to be the driving factor that induces change in residue-pair rigidity correlations within SB transposase. Given that a signature RNase H-like structural motif is found in DD[E/D]-transposases and, more broadly, in a large superfamily of polynucleotidyl transferases, our results are relevant to these proteins as well.

LARGE GRAIN Encodes a Putative RNA-Binding Protein that Regulates Spikelet Hull Length in Rice.

Grain size is a key determiner of grain weight, one of the yield components in rice (Oryza sativa). Therefore, to increase grain yield, it is important to elucidate the detailed mechanisms regulating grain size. The Large grain (Lgg) mutant, found in the nonautonomous DNA-based active rice transposon1 (nDart1)-tagged lines of Koshihikari, is caused by a truncated nDart1-3 and 355 bp deletion in the 5' untranslated region of LGG, which encodes a putative RNA-binding protein, through transposon display and cosegregation analysis between grain length and LGG genotype in F2 and F3. Clustered regularly interspaced short palindromic repeats/CRISPR-associated 9-mediated knockout and overexpression of LGG led to longer and shorter grains than wild type, respectively, showing that LGG regulates spikelet hull length. Expression of LGG was highest in the 0.6-mm-long young panicle and gradually decreased as the panicle elongated. LGG was also expressed in roots and leaves. These results show that LGG functions at the very early stage of panicle development. Longitudinal cell numbers of spikelet hulls of Lgg, knockout and overexpressed plants were significantly different from those of the wild type, suggesting that LGG might regulate longitudinal cell proliferation in the spikelet hull. RNA-Seq analysis of 1-mm-long young panicles from LGG knockout and overexpressing plants revealed that the expressions of many cell cycle-related genes were reduced in knockout plants relative to LGG-overexpressing plants and wild type, whereas some genes for cell proliferation were highly expressed in knockout plants. Taken together, these results suggest that LGG might be a regulator of cell cycle and cell division in the rice spikelet hull.

Transgenerational activation of an autonomous DNA transposon, Dart1-24, by 5-azaC treatment in rice.

KEY MESSAGE: Dart1-24, one of the 37 autonomous DNA transposon Dart1s, was heritably activated by the demethylation of the 5' region following 5-azaC treatment of rice seeds. Transposons are controlled by epigenetic regulations. To obtain newly activated autonomous elements of Dart1, a DNA transposon, in rice, seeds of a stable pale yellow leaf (pyl-stb) mutant caused by the insertion of nDart1-0, a nonautonomous element in OsClpP5, were treated with 5-azaC, a demethylating agent. In the 5-azaC-treated M1 plants, 60-70% of the plants displayed variegated pale yellow leaf (pyl-v) phenotype, depending on the concentration of 5-azaC used, suggesting that inactivated Dart1 might become highly activated by 5-azaC treatment and nDart1-0 was excised from OsClpP5 by the activated Dart1s. Although the M2 plants derived from most of these pyl-v plants showed stable pyl phenotypes, some variegated M1 plants generated pyl-v M2 progeny. These results indicated that most M1 pyl-v phenotypes at M1 were not heritable. Dart1-24, 1-27 and 1-28 were expressed in the M2 pyl-v plants, and mapping analysis confirmed that Dart1-24 was newly activated. Further, the transgenerational activation of Dart1-24 was demonstrated to be caused by the demethylation of nucleotides in its 5' region.

Establishment of nDart1-tagged lines of Koshihikari, an elite variety of rice in Japan.

To utilize a transposon-tagged mutant as a breeding material in rice, an endogenous DNA transposon, nDart1-0, was introduced into Koshihikari by successive backcrossing together with aDart1-27, an active autonomous element. The founder line for nDart1-tagged lines of Koshihikari carried nDart1-0 on chromosome 9 and transposed nDart1-12s on chromosomes 1 and 8 and nDart1-3 on chromosome 11. In nDart1-tagged lines, there were the most abnormal phenotypic mutants and many aberrant chlorophyll mutants at seedling stage. At mature stage, many semi-sterile mutants were observed. Dwarf, reduced culm number and lesion mimic mutants were also found. In total, 43.2% of the lines segregated some phenotypic mutants. Thus, the nDart1-tagged lines of Koshihikari are expected to be potentially useful for screening stress-tolerant mutants under abiotic or biotic stress conditions.

Transcriptional activity of PIF and Pong-like Class II transposable elements in Triticeae.

BACKGROUND: Transposable elements are major contributors to genome size and variability, accounting for approximately 70-80% of the maize, barley, and wheat genomes. PIF and Pong-like elements belong to two closely-related element families within the PIF/Harbinger superfamily of Class II (DNA) transposons. Both elements contain two open reading frames; one encodes a transposase (ORF2) that catalyzes transposition of the functional elements and their related non-autonomous elements, while the function of the second is still debated. In this work, we surveyed for PIF- and Pong-related transcriptional activity in 13 diploid Triticeae species, all of which have been previously shown to harbor extensive within-genome diversity of both groups of elements. RESULTS: The results revealed that PIF elements have considerable transcriptional activity in Triticeae, suggesting that they can escape the initial levels of plant cell control and are regulated at the post-transcriptional level. Phylogenetic analysis of 156 PIF cDNA transposase fragments along with 240 genomic partial transposase sequences showed that most, if not all, PIF clades are transcriptionally competent, and that multiple transposases coexisting within a single genome have the potential to act simultaneously. In contrast, we did not detect any transcriptional activity of Pong elements in any sample. CONCLUSIONS: The lack of Pong element transcription shows that even closely related transposon families can exhibit wide variation in their transposase transcriptional activity within the same genome.

Evolutionary histories of transposable elements in the genome of the largest living marsupial carnivore, the Tasmanian devil.

The largest living carnivorous marsupial, the Tasmanian devil (Sarcophilus harrisii), is the sole survivor of a lineage originating about 12 Ma. We set out to investigate the spectrum of transposable elements found in the Tasmanian devil genome, the first high-coverage genome of an Australian marsupial. Marsupial genomes have been shown to have the highest amount of transposable elements among vertebrates. We analyzed the horizontally transmitted DNA transposons OC1 and hAT-1_MEu in the Tasmanian devil genome. OC1 is present in all carnivorous marsupials, while having a very limited distribution among the remaining Australian marsupial orders. In contrast, hAT-1_MEu is present in all Australian marsupial orders, and has so far only been identified in a few placental mammals. We screened 158 introns for phylogenetically informative retrotransposons in the order Dasyuromorphia, and found that the youngest SINE (Short INterspersed Element), WSINE1, is no longer active in the subfamily Dasyuridae. The lack of detectable WSINE1 activity in this group may be due to a retrotransposon inactivation event approximately 30 Ma. We found that the Tasmanian devil genome contains a relatively low number of continuous full-length LINE-1 (Long INterspersed Element 1, L1) retrotransposons compared with the opossum genome. Furthermore, all L1 elements in the Tasmanian devil appeared to be nonfunctional. Hidden Markov Model approaches suggested that other potential sources of functional reverse transcriptase are absent from the genome. We discuss the issues associated with assembling long, highly similar L1 copies from short read Illumina data and describe how assembly artifacts can potentially lead to erroneous conclusions.

Diversity, abundance, and evolutionary dynamics of Pong-like transposable elements in Triticeae.

Pong-like elements are members of the PIF/Harbinger superfamily of DNA transposons that has been described in many plants, animals, and fungi. Most Pong elements contain two open reading frames (ORFs). One encodes a transposase (ORF2) that catalyzes transposition of Pong and related non-autonomous elements, while the function of the second is unknown. Little is known about the evolutionary history of Pong elements in flowering plants. In this work, we present the first comprehensive analysis of the diversity, abundance, and evolution of the Pong-like transposase gene in the genomes of 21 diploid species from the wheat tribe, Triticeae, and we present the first convincing evidence of horizontal transfer of nuclear-encoded Pong elements in any organism. A phylogenetic analysis of nearly 300 Pong sequences based on a conserved region of the transposase domain revealed a complex evolutionary history of Pong elements that can be best explained by ancestral polymorphism, followed by differential evolutionary success of some transposase lineages, and by occasional horizontal transfer between phylogenetically distant genera. In addition, we used transposon display to estimate the abundance of the transposase gene within Triticeae genomes, and our results revealed varying levels of Pong proliferation, with numbers of transposase copies ranging from 22 to 92. Comparisons of Pong transposase abundance to flow cytometry estimates of genome size revealed that larger Triticeae genome size was not correlated with transposase abundance.

A mutable albino allele in rice reveals that formation of thylakoid membranes requires the SNOW-WHITE LEAF1 gene.

Active DNA transposons are important tools for gene functional analysis. The endogenous non-autonomous transposon, nDart1-0, in rice (Oryza sativa L.) is expected to generate various transposon-insertion mutants because nDart1-0 elements tend to insert into genic regions under natural growth conditions. We have developed a specific method (nDart1-0-iPCR) for efficient detection of nDart1-0 insertions and successfully identified the SNOW-WHITE LEAF1 (SWL1) gene in a variegated albino (swl1-v) mutant obtained from the nDart1-promoted rice tagging line. The variegated albino phenotype was caused by insertion and excision of nDart1-0 in the 5'-untranslated region of the SWL1 gene predicted to encode an unknown protein with the N-terminal chloroplast transit peptide. SWL1 expression was detected in various rice tissues at different developmental stages. However, immunoblot analysis indicated that SWL1 protein accumulation was strictly regulated in a tissue-specific manner. In the swl1 mutant, formations of grana and stroma thylakoids and prolamellar bodies were inhibited. This study revealed that SWL1 is essential for the beginning of thylakoid membrane organization during chloroplast development. Furthermore, we provide a developmental perspective on the nDart1-promoted tagging line to characterize unidentified gene functions in rice.

Scatter: a novel family of miniature inverted-repeat transposable elements in the fungus Botrytis cinerea.

Miniature inverted-repeat transposable elements (MITEs) are short non-autonomous DNA transposons that play an important role in genome structure and function. Here, we described a novel family of MITEs, named Scatter, identified from the genomes of three strains of the fungus Botrytis cinerea (T4, B05.10, and TBC-A). Intact Scatter elements are typically an average of 247 bp, and contain 41 bp terminal inverted repeats (TIRs) and 2-bp "TA" target site duplications (TSDs). Based a search against the transposable elements database and GenBank, Scatter is a novel and potentially species-specific family of MITEs. Moderate heterogeneity in sequence and size of individual Scatter copies suggests that Scatter elements were not recently proliferated. Most integrated sites were conserved across all three strains tested and elements inserted at equivalent sites shared high identity at the nucleotide level. This conservation, in combination with the presence of a similar copy number (22-24), in B. cinerea strains tested suggests that Scatter may be a relic of an ancient transposition developed prior to the strain divergence of B. cinerea. Two unique insertion instances were observed, indicating that some copies of Scatter may have remained active following strain divergence of B. cinerea. Because only a few subtle insertion differences among B. cinerea strains were observed, Scatter may play only a minor role in the genetic diversity in B. cinerea species. Most Scatter elements appear to be inserted in potential regulatory regions of adjacent coding regions, highlighting their role in transcriptional regulation. The origin of Scatter remains to be addressed. Scatter is the first well-characterized family of MITEs in B. cinerea.

Fake IDs? Widespread misannotation of DNA transposons as a general transcription factor.

Accurate annotation of genes and transposable elements (TEs) is vital for understanding genomes, but current annotation pipelines often misannotate TEs as genes. This study reveals how the general transcription factor II-I repeat domain-containing protein 2 (GTF2IRD2) erroneously annotated DNA transposons in non-mammalian species, as it contains a 3' fused hAT transposase domain. We also demonstrate the generality of this problem by identifying misannotated TEs as genes in other vertebrate genomes. Such misannotations can lead to errors in phylogenetic analyses and wasted time for investigators. The study proposes adding a final TE-check to gene annotation pipelines to mitigate this problem.

TEMM: A Curated Data Resource for Transposon Element-Based Molecular Markers in Plants.

Transposon elements (TEs) are mobile genetic elements that can insert themselves into new locations and modify the plant genome. In recent years, they have been used as molecular markers in plant breeding programs. TE-based molecular markers (TE-markers) are divided into two categories depending on the transcription mechanism of the TEs. The first category is retrotransposon-based molecular markers, which include RBIP, IRAP, REMAP, and iPBS. The second group is DNA-based-TE-markers, which include MITE, TE-junction, and CACTA TE-markers. These markers are a good tool for studying genetic diversity and can provide information on plants' phylogenetic and evolutionary history. They can help improve breeding programs to increase agronomic traits and develop new varieties. Overall, TE-markers play an important role in plant genetics and plant breeding and contribute to a better understanding of plant biology. Here, we present TEMM, a curated data resource for TE-markers in plants. Relevant research articles were screened to collect primer sequences and related information. Only articles containing primer sequences are added to the present data resource. TEMM contains 784 primers with their associated PCR reaction programs and their applications in various crops. These include 203 IPBS, 191 RBIP, 140 IRAP, 78 TE-junction, 76 IRAPS, 47 RBIP-IRAP, 16 IRAP-REMAP, 12 REMAP, 12 REMA-IRAP, 6 REMA, and 3 ISBP primers. The data resource is freely available at https://bioinformatics.um6p.ma/TEMM .

Prolonged activity of the transposase helper may raise safety concerns during DNA transposon-based gene therapy.

DNA transposon-based gene delivery vectors represent a promising new branch of randomly integrating vector development for gene therapy. For the side-by-side evaluation of the piggyBac and Sleeping Beauty systems-the only DNA transposons currently employed in clinical trials-during therapeutic intervention, we treated the mouse model of tyrosinemia type I with liver-targeted gene delivery using both transposon vectors. For genome-wide mapping of transposon insertion sites we developed a new next-generation sequencing procedure called streptavidin-based enrichment sequencing, which allowed us to identify approximately one million integration sites for both systems. We revealed that a high proportion of piggyBac integrations are clustered in hot regions and found that they are frequently recurring at the same genomic positions among treated animals, indicating that the genome-wide distribution of Sleeping Beauty-generated integrations is closer to random. We also revealed that the piggyBac transposase protein exhibits prolonged activity, which predicts the risk of oncogenesis by generating chromosomal double-strand breaks. Safety concerns associated with prolonged transpositional activity draw attention to the importance of squeezing the active state of the transposase enzymes into a narrower time window.

IS481EU Shows a New Connection between Eukaryotic and Prokaryotic DNA Transposons.

DDD/E transposase gene is the most abundant gene in nature and many DNA transposons in all three domains of life use it for their transposition. A substantial number of eukaryotic DNA transposons show similarity to prokaryotic insertion sequences (ISs). The presence of IS481-like DNA transposons was indicated in the genome of Trichomonas vaginalis. Here, we surveyed IS481-like eukaryotic sequences using a bioinformatics approach and report a group of eukaryotic IS481-like DNA transposons, designated IS481EU, from parabasalids including T. vaginalis. The lengths of target site duplications (TSDs) of IS481EU are around 4 bps, around 15 bps, or around 25 bps, and strikingly, these discrete lengths of TSDs can be observed even in a single IS481EU family. Phylogenetic analysis indicated the close relationships of IS481EU with some of the prokaryotic IS481 family members. IS481EU was not well separated from IS3EU/GingerRoot in the phylogenetic analysis, but was distinct from other eukaryotic DNA transposons including Ginger1 and Ginger2. The unique characteristics of IS481EU in protein sequences and the distribution of TSD lengths support its placement as a new superfamily of eukaryotic DNA transposons.

Diversity and Evolution of DNA Transposons Targeting Multicopy Small RNA Genes from Actinopterygian Fish.

Dada is a unique superfamily of DNA transposons, inserted specifically in multicopy RNA genes. The zebrafish genome harbors five families of Dada transposons, whose targets are U6 and U1 snRNA genes, and tRNA-Ala and tRNA-Leu genes. Dada-U6, which is inserted specifically in U6 snRNA genes, is found in four animal phyla, but other target-specific lineages have been reported only from one or two species. Here, vertebrate genomes and transcriptomes were surveyed to characterize Dada families with new target specificities, and over 120 Dada families were characterized from the genomes of actinopterygian fish. They were classified into 12 groups with confirmed target specificities. Newly characterized Dada families target tRNA genes for Asp, Asn, Arg, Gly, Lys, Ser, Tyr, and Val, and 5S rRNA genes. Targeted positions inside of tRNA genes are concentrated in two regions: around the anticodon and the A box of RNA polymerase III promoter. Phylogenetic analysis revealed the relationships among actinopterygian Dada families, and one domestication event in the common ancestor of carps and minnows belonging to Cyprinoidei, Cypriniformes. Sequences targeted by phylogenetically related Dada families show sequence similarities, indicating that the target specificity of Dada is accomplished through the recognition of primary nucleotide sequences.

CRISPR-Cas9-directed gene tagging using a single integrase-defective lentiviral vector carrying a transposase-based Cas9 off switch.

Locus-directed DNA cleavage induced by the CRISPR-Cas9 system triggers DNA repair mechanisms allowing gene repair or targeted insertion of foreign DNA. For gene insertion to be successful, availability of a homologous donor template needs to be timed with cleavage of the DNA by the Cas9 endonuclease guided by a target-specific single guide RNA (sgRNA). We present a novel approach for targeted gene insertion based on a single integrase-defective lentiviral vector (IDLV) carrying a Cas9 off switch. Gene insertion using this approach benefits from transposon-based stable Cas9 expression, which is switched off by excision-only transposase protein co-delivered in IDLV particles carrying a combined sgRNA/donor vector. This one-vector approach supports potent (up to >80%) knockin of a full-length EGFP gene sequence. This traceless cell engineering method benefits from high stable levels of Cas9, timed intracellular availability of the molecular tools, and a built-in feature to turn off Cas9 expression after DNA cleavage. The simple technique is based on transduction with a single IDLV, which holds the capacity to transfer larger donor templates, allowing robust gene knockin or tagging of genes in a single step.

A unique eukaryotic lineage of composite-like DNA transposons encoding a DDD/E transposase and a His-Me finger homing endonuclease.

BACKGROUND: DNA transposons are ubiquitous components of eukaryotic genomes. A major group of them encode a DDD/E transposase and contain terminal inverted repeats (TIRs) of varying lengths. The Kolobok superfamily of DNA transposons has been found in a wide spectrum of organisms. RESULTS: Here we report a new Kolobok lineage, designated KolobokP. They were identified in 7 animal phyla (Mollusca, Phoronida, Annelida, Nemertea, Bryozoa, Chordata, and Echinodermata), and are especially rich in bivalves. Unlike other Kolobok families, KolobokP adopts a composite-like architecture: an internal region (INT) flanked by two long terminal direct repeats (LTDRs), which exhibit their own short terminal inverted repeats ranging up to 18 bps. The excision of LTDRs was strongly suggested. The LTDR lengths seem to be constrained to be either around 450-bp or around 660-bp. The internal region encodes a DDD/E transposase and a small His-Me finger nuclease, which likely originated from the homing endonuclease encoded by a group I intron from a eukaryotic species. The architecture of KolobokP resembles composite DNA transposons, usually observed in bacterial genomes, and long terminal repeat (LTR) retrotransposons. In addition to this monomeric LTDR-INT-LTDR structure, plenty of solo LTDRs and multimers represented as (LTDR-INT)n-LTDR are also observed. Our structural and phylogenetic analysis supported the birth of KolobokP in the late stage of the Kolobok evolution. We propose KolobokP families propagate themselves in two ways: the canonical transposition catalyzed by their transposase and the sequence-specific cleavage by their endonuclease followed by the multimerization through the unequal crossover. CONCLUSIONS: The presence of homing endonuclease and long terminal direct repeats of KolobokP families suggest their unique dual replication mechanisms: transposition and induced unequal crossover.