Evolution of pogo, a separate superfamily of IS630-Tc1-mariner transposons, revealing recurrent domestication events in vertebrates

Unveiling the Genetic Blueprint of a Desert Scorpion: A Chromosome-level Genome of Hadrurus arizonensis Provides the First Reference for Parvorder Iurida

Over 400 million years old, scorpions represent an ancient group of arachnids and one of the first animals to adapt to life on land. Presently, the lack of available genomes within scorpions hinders research on their evolution. This study leverages ultralong nanopore sequencing and Pore-C to generate the first chromosome-level assembly and annotation for the desert hairy scorpion, Hadrurus arizonensis. The assembled genome is 2.23 Gb in size with an N50 of 280 Mb. Pore-C scaffolding reoriented 99.6% of bases into nine chromosomes and BUSCO identified 998 (98.6%) complete arthropod single copy orthologs. Repetitive elements represent 54.69% of the assembled bases, including 872,874 (29.39%) LINE elements. A total of 18,996 protein-coding genes and 75,256 transcripts were predicted, and extracted protein sequences yielded a BUSCO score of 97.2%. This is the first genome assembled and annotated within the family Hadruridae, representing a crucial resource for closing gaps in genomic knowledge of scorpions, resolving arachnid phylogeny, and advancing studies in comparative and functional genomics.

PiggyBac Transposon Mining in the Small Genomes of Animals

TEs, including DNA transposons, are major contributors of genome expansions, and have played a very significant role in shaping the evolution of animal genomes, due to their capacity to jump from one genomic position to the other. In this study, we investigated the evolution landscapes of PB transposons, including their distribution, diversity, activity and structure organization in 79 species of small (compact) genomes of animals comprising both vertebrate and invertebrates. Overall, 212 PB transposon types were detected from almost half (37) of the total number of the small genome species (79) investigated. The detected PB transposon types, which were unevenly distributed in various genera and phyla, have been classified into seven distinct clades or families with good bootstrap support (>80%). The PB transposon types that were identified have a length ranging from 1.23 kb to 9.51 kb. They encode transposases of approximately ≥500 amino acids in length, and possess terminal inverted repeats (TIRs) ranging from 4 bp to 24 bp. Though some of the transposon types have long TIRs (528 bp), they still maintain the consistent and reliable 4 bp target site duplication (TSD) of TTAA. However, PiggyBac-2_Cvir transposon originating from the Crassostrea virginica species exhibits a unique TSD of TATG. The TIRs of the transposons in all the seven families display high divergence, with a highly conserved 5' end motif. The core transposase domains (DDD) were better conserved among the seven different families compared to the other protein domains, which were less prevalent in the vertebrate genome. The divergent evolution dynamics analysis also indicated that the majority of the PB transposon types identified in this study are either relatively young or old, with some being active. Additionally, numerous invasions of PB transposons were found in the genomes of both vertebrate and invertebrate animals. The data reveals that the PB superfamily is widely distributed in these species. PB transposons exhibit high diversity and activity in the small genomes of animals, and might play a crucial role in shaping the evolution of these small genomes of animals.

Hiker, a new family of DNA transposons encoding transposases with DD35E motifs, displays a distinct phylogenetic relationship with most known DNA transposon families of IS630-Tc1-mariner (ITm)

DNA transposons play a crucial role in determining the size and structure of eukaryotic genomes. In this study, a new family of IS630-Tc1-mariner (ITm) DNA transposons, named Hiker (HK), was identified. HK is characterized by a DD35E catalytic domain and is distinct from all previously known families of the ITm group. Phylogenetic analyses showed that DD35E/Hiker forms a monophyletic clade with DD34E/Gambol, indicating that they may represent a separate superfamily of ITm. A total of 178 Hiker species were identified, with 170 found mainly in Actinopterygii, one in Chondrichthyes, six in Anura and one in Mollusca. Gambol (GM), on the other hand, are found in invertebrates, with 18 in Arthropoda and one in Platyhelminthes. Hiker transposons have a total length ranging from 2.14 to 3.67 kb and contain a single open reading frame that encodes a protein of approximately 370 amino acids (range 311-413 aa). They are flanked by short terminal inverted repeats (TIRs) of 16-30 base pairs and two base pair (TA) target-site duplications. In contrast, most transposons of the Gambol family have a total length of 1.35-5.96 kb, encode a transposase protein of approximately 350 amino acids (range 306-374 aa), and are flanked by TIRs that range from 32 to 1097 bp in length. Both Hiker and Gambol transposases have several conserved motifs, including helix-turn-helix (HTH) motifs and a DDE domain. Our study observed multiple amplification waves and repeated horizontal transfer (HT) events of HK transposons in vertebrate genomes, indicating their role in diversifying and shaping the genomes of Actinopterygii, Chondrichthyes, and Anura. Conversely, GM transposons showed few Horizontal transfer events. According to cell-based transposition assays, most HK transposons are likely inactive due to the truncated DNA binding domains of their transposases. We present an updated classification of the ITm group based on these findings, which will enhance the understanding of both the evolution of ITm transposons and that of their hosts.

Whole Genome Assembly and Annotation of Blackstripe Livebearer Poeciliopsis prolifica

The blackstripe livebearer Poeciliopsis prolifica is a live-bearing fish belonging to the family Poeciliidae with high level of postfertilization maternal investment (matrotrophy). This viviparous matrotrophic species has evolved a structure similarly to the mammalian placenta. Placentas have independently evolved multiple times in Poeciliidae from nonplacental ancestors, which provide an opportunity to study the placental evolution. However, there is a lack of high-quality reference genomes for the placental species in Poeciliidae. In this study, we present a 674 Mb assembly of P. prolifica in 504 contigs with excellent continuity (contig N50 7.7 Mb) and completeness (97.2% Benchmarking Universal Single-Copy Orthologs [BUSCO] completeness score, including 92.6% single-copy and 4.6% duplicated BUSCO score). A total of 27,227 protein-coding genes were annotated from the merged datasets based on bioinformatic prediction, RNA sequencing and homology evidence. Phylogenomic analyses revealed that P. prolifica diverged from the guppy (Poecilia reticulata) ∼19 Ma. Our research provides the necessary resources and the genomic toolkit for investigating the genetic underpinning of placentation.

Mosquito (MS), a DD37E Family of Tc1/Mariner, Displaying a Distinct Evolution Profile from DD37E/TRT and DD37E/L18

Diverse Tc1/mariner elements with the DD37E signature have been detected. However, their evolutionary relationship and profiles are largely unknown. Using bioinformatics methods, we defined the evolution profile of a Tc1/Mariner family, which harbors the catalytic domain with the DD37E signature, and renamed it DD37E/Mosquito (MS). MS transposons form a separate monophyletic clade in the phylogenetic tree, distinct from the other two groups of elements with the DD37E signature, DD37E/L18 and DD37E/TRT (transposon related to Tc1), and represent a very different taxonomic distribution from that of DD37E/TRT. MS is only detected in invertebrate and is mostly present in Arthropoda, as well as in Cnidaria, Ctenophora, Mollusca, Nematoda, and Platyhelminthes, with a total length of about 1.3 kb, containing an open reading frame (ORF) encoding about 340 amino acids transposases, with a conserved DD37E catalytic domain. The terminal inverted repeat (TIR) lengths range from 19 bp to 203 bp, and the target site duplication (TSD) is TA. We also identified few occurrences of MS horizontal transfers (HT) across lineages of diptera. In this paper, the distribution characteristics, structural characteristics, phylogenetic evolution, and horizontal transfer of the MS family are fully analyzed, which is conducive to supplementing and improving the Tc1/Mariner superfamily and excavating active transposons.

A high-quality reference genome for the fission yeast Schizosaccharomyces osmophilus

Fission yeasts are an ancient group of fungal species that diverged from each other from tens to hundreds of million years ago. Among them is the preeminent model organism Schizosaccharomyces pombe, which has significantly contributed to our understandings of molecular mechanisms underlying fundamental cellular processes. The availability of the genomes of S. pombe and 3 other fission yeast species S. japonicus, S. octosporus, and S. cryophilus has enabled cross-species comparisons that provide insights into the evolution of genes, pathways, and genomes. Here, we performed genome sequencing on the type strain of the recently identified fission yeast species S. osmophilus and obtained a complete mitochondrial genome and a nuclear genome assembly with gaps only at rRNA gene arrays. A total of 5,098 protein-coding nuclear genes were annotated and orthologs for more than 95% of them were identified. Genome-based phylogenetic analysis showed that S. osmophilus is most closely related to S. octosporus and these 2 species diverged around 16 million years ago. To demonstrate the utility of this S. osmophilus reference genome, we conducted cross-species comparative analyses of centromeres, telomeres, transposons, the mating-type region, Cbp1 family proteins, and mitochondrial genomes. These analyses revealed conservation of repeat arrangements and sequence motifs in centromere cores, identified telomeric sequences composed of 2 types of repeats, delineated relationships among Tf1/sushi group retrotransposons, characterized the evolutionary origins and trajectories of Cbp1 family domesticated transposases, and discovered signs of interspecific transfer of 2 types of mitochondrial selfish elements.

A draft of the genome of the Gulf Coast tick, Amblyomma maculatum

The Gulf Coast tick, Amblyomma maculatum, inhabits the Southeastern states of the USA bordering the Gulf of Mexico, Mexico, and other Central and South American countries. More recently, its U.S. range has extended West to Arizona and Northeast to New York state and Connecticut. It is a vector of Rickettsia parkeri and Hepatozoon americanum. This tick species has become a model to study tick/Rickettsia interactions. To increase our knowledge of the basic biology of A. maculatum we report here a draft genome of this tick and an extensive functional classification of its proteome. The DNA from a single male tick was used as a genomic source, and a 10X genomics protocol determined 28,460 scaffolds having equal or more than 10 Kb, totaling 1.98 Gb. The N50 scaffold size was 19,849 Kb. The BRAKER pipeline was used to find the protein-coding gene boundaries on the assembled A. maculatum genome, discovering 237,921 CDS. After trimming and classifying the transposable elements, bacterial contaminants, and truncated genes, a set of 25,702 were annotated and classified as the core gene products. A BUSCO analysis revealed 83.4% complete BUSCOs. A hyperlinked spreadsheet is provided, allowing browsing of the individual gene products and their matches to several databases.

Passer, a highly active transposon from a fish genome, as a potential new robust genetic manipulation tool

The discovery of new, active DNA transposons can expand the range of genetic tools and provide more options for genomic manipulation. In this study, a bioinformatics analysis suggested that Passer (PS) transposons, which are members of the pogo superfamily, show signs of recent and current activity in animals and may be active in some species. Cell-based transposition assays revealed that the native PS transposases from Gasterosteus aculeatus and Danio rerio displayed very high activity in human cells relative to the Sleeping Beauty transposon. A typical overproduction inhibition phenomenon was observed for PS, and transposition capacity was decreased by ∼12% with each kilobase increase in the insertion size. Furthermore, PS exhibited a pronounced integration preference for genes and their transcriptional regulatory regions. We further show that two domesticated human proteins derived from PS transposases have lost their transposition activity. Overall, PS may represent an alternative with a potentially efficient genetic manipulation tool for transgenesis and mutagenesis applications.

MITE infestation accommodated by genome editing in the germline genome of the ciliate Blepharisma

During their development following sexual conjugation, ciliates excise numerous internal eliminated sequences (IESs) from a copy of the germline genome to produce the functional somatic genome. Most IESs are thought to have originated from transposons, but the presumed homology is often obscured by sequence decay. To obtain more representative perspectives on the nature of IESs and ciliate genome editing, we assembled 40,000 IESs of Blepharisma stoltei, a species belonging to a lineage (Heterotrichea) that diverged early from those of the intensively studied model ciliate species. About a quarter of IESs were short (<115 bp), largely nonrepetitive, and with a pronounced ~10 bp periodicity in length; the remainder were longer (up to 7 kbp) and nonperiodic and contained abundant interspersed repeats. Contrary to the expectation from current models, the assembled Blepharisma germline genome encodes few transposases. Instead, its most abundant repeat (8,000 copies) is a Miniature Inverted-repeat Transposable Element (MITE), apparently a deletion derivative of a germline-limited Pogo-family transposon. We hypothesize that MITEs are an important source of IESs whose proliferation is eventually self-limiting and that rather than defending the germline genomes against mobile elements, transposase domestication actually facilitates the accumulation of junk DNA.

Transposon-derived transcription factors across metazoans

Transposable elements (TE) could serve as sources of new transcription factors (TFs) in plants and some other model species, but such evidence is lacking for most animal lineages. Here, we discovered multiple independent co-options of TEs to generate 788 TFs across Metazoa, including all early-branching animal lineages. Six of ten superfamilies of DNA transposon-derived conserved TF families (ZBED, CENPB, FHY3, HTH-Psq, THAP, and FLYWCH) were identified across nine phyla encompassing the entire metazoan phylogeny. The most extensive convergent domestication of potentially TE-derived TFs occurred in the hydroid polyps, polychaete worms, cephalopods, oysters, and sea slugs. Phylogenetic reconstructions showed species-specific clustering and lineage-specific expansion; none of the identified TE-derived TFs revealed homologs in their closest neighbors. Together, our study established a framework for categorizing TE-derived TFs and informing the origins of novel genes across phyla.

Horizontal Transfer and Evolutionary Profiles of Two Tc1/DD34E Transposons (ZB and SB) in Vertebrates

Both ZeBrafish (ZB), a recently identified DNA transposon in the zebrafish genome, and SB, a reconstructed transposon originally discovered in several fish species, are known to exhibit high transposition activity in vertebrate cells. Although a similar structural organization was observed for ZB and SB transposons, the evolutionary profiles of their homologs in various species remain unknown. In the present study, we compared their taxonomic ranges, structural arrangements, sequence identities, evolution dynamics, and horizontal transfer occurrences in vertebrates. In total, 629 ZB and 366 SB homologs were obtained and classified into four distinct clades, named ZB, ZB-like, SB, and SB-like. They displayed narrow taxonomic distributions in eukaryotes, and were mostly found in vertebrates, Actinopterygii in particular tended to be the major reservoir hosts of these transposons. Similar structural features and high sequence identities were observed for transposons and transposase, notably homologous to the SB and ZB elements. The genomic sequences that flank the ZB and SB transposons in the genomes revealed highly conserved integration profiles with strong preferential integration into AT repeats. Both SB and ZB transposons experienced horizontal transfer (HT) events, which were most common in Actinopterygii. Our current study helps to increase our understanding of the evolutionary properties and histories of SB and ZB transposon families in animals.

High Expression of POGK Predicts Poor Prognosis in Patients with Hepatocellular Carcinoma

OBJECTIVE

Kruppel-associated box (KRAB) proteins reportedly play a dual role in neoplastic transformation. At present, little is known about the function of the proteins encoded by the human pogo transposable element derived with KRAB domain (POGK) gene. Herein, we evaluated the prognostic significance of POGK expression in patients with hepatocellular carcinoma (HCC).

METHODS

The data of HCC patients was downloaded from The Cancer Genome Atlas (TCGA) database. To determine the relationship between POGK and clinical features, logistic regression was applied. Cox regression and Kaplan-Meier analyses were used to evaluate the correlation between POGK and survival rates. Gene ontology (GO) analysis and Gene set enrichment analysis (GSEA) were conducted to identify the enriched pathways and functions associated with POGK.

RESULTS

A total of 374 HCC patients were identified in TCGA. POGK was significantly upregulated in HCC and correlated with tumor status (p = 0.036), race (p = 0.025), weight (p = 0.002), body mass index (p = 0.033), histologic grade (p < 0.001), and alpha-fetoprotein (p < 0.001). High POGK expression in HCC patients correlated with a poor outcome in terms of overall survival (p = 0.0018), progression-free survival (p = 0.0087), relapse-free survival (p = 0.045), and disease-specific survival (p = 0.014), according to Kaplan-Meier analysis. Receiver operating characteristic curve analysis showed that the area under the curve of POGK expression for HCC diagnosis was 0.891. GSEA showed that high POGK expression might activate mitotic prometaphase, kinesins, homologous DNA pairing and strand exchange, MET activates PTK2 signaling pathway, G1 to S cell cycle control, Aurora B pathway, ncRNAs involved in WNT signaling pathway, hepatitis C, and ncRNAs involved in the STAT3 signaling pathway. POGK expression correlated with the abundance of adaptive and innate immunocytes in HCC.

CONCLUSION

High expression of POGK has high diagnostic and prognostic values in patients with HCC. Moreover, POGK expression is correlated with immune infiltration in HCC.

The TIGD5 gene located in 8q24 and frequently amplified in ovarian cancers is a tumor suppressor

Ovarian cancer (OC) is the fifth most common cancer of female cancer death and leading cause of lethal gynecological cancers. High-grade serous ovarian carcinoma (HGSOC) is an aggressive malignancy that is rapidly fatal. Many cases of OC show amplification of the 8q24 chromosomal region, which contains the well-known oncogene MYC. Although MYC amplification is more frequently observed in OCs than in other tumor types, due to the large size of the 8q24 amplicon, the functions of the vast majority of the genes it contains are still unknown. The TIGD5 gene is located at 8q24.3 and encodes a nuclear protein with a DNA-binding motif, but its precise role is obscure. We show here that TIGD5 often co-amplifies with MYC in OCs, and that OC patients with high TIGD5 mRNA expression have a poor prognosis. However, we also found that TIGD5 overexpression in ovarian cancer cell lines unexpectedly suppressed their growth, adhesion, and invasion in vitro, and also reduced tumor growth in xenografted nude mice in vivo. Thus, our work suggests that TIGD5 may in fact operate as a tumor suppressor in OCs rather than as an oncogene.

Identification of putative enhancer-like elements predicts regulatory networks active in planarian adult stem cells

Planarians have become an established model system to study regeneration and stem cells, but the regulatory elements in the genome remain almost entirely undescribed. Here, by integrating epigenetic and expression data we use multiple sources of evidence to predict enhancer elements active in the adult stem cell populations that drive regeneration. We have used ChIP-seq data to identify genomic regions with histone modifications consistent with enhancer activity, and ATAC-seq data to identify accessible chromatin. Overlapping these signals allowed for the identification of a set of high-confidence candidate enhancers predicted to be active in planarian adult stem cells. These enhancers are enriched for predicted transcription factor (TF) binding sites for TFs and TF families expressed in planarian adult stem cells. Footprinting analyses provided further evidence that these potential TF binding sites are likely to be occupied in adult stem cells. We integrated these analyses to build testable hypotheses for the regulatory function of TFs in stem cells, both with respect to how pluripotency might be regulated, and to how lineage differentiation programs are controlled. We found that our predicted GRNs were independently supported by existing TF RNAi/RNA-seq datasets, providing further evidence that our work predicts active enhancers that regulate adult stem cells and regenerative mechanisms.

Evolution of eukaryotic centromeres by drive and suppression of selfish genetic elements

Despite the universal requirement for faithful chromosome segregation, eukaryotic centromeres are rapidly evolving. It is hypothesized that rapid centromere evolution represents an evolutionary arms race between selfish genetic elements that drive, or propagate at the expense of organismal fitness, and mechanisms that suppress fitness costs. Selfish centromere DNA achieves preferential inheritance in female meiosis by recruiting more effector proteins that alter spindle microtubule interaction dynamics. Parallel pathways for effector recruitment are adaptively evolved to suppress functional differences between centromeres. Opportunities to drive are not limited to female meiosis, and selfish transposons, plasmids and B chromosomes also benefit by maximizing their inheritance. Rapid evolution of selfish genetic elements can diversify suppressor mechanisms in different species that may cause hybrid incompatibility.

Chromosome-level assembly and annotation of the blue catfish Ictalurus furcatus, an aquaculture species for hybrid catfish reproduction, epigenetics, and heterosis studies

Background

The blue catfish is of great value in aquaculture and recreational fisheries. The F₁ hybrids of female channel catfish (Ictalurus punctatus) × male blue catfish (Ictalurusfurcatus) have been the primary driver of US catfish production in recent years because of superior growth, survival, and carcass yield. The channel–blue hybrid also provides an excellent model to investigate molecular mechanisms of environment-dependent heterosis. However, transcriptome and methylome studies suffered from low alignment rates to the channel catfish genome due to divergence, and the genome resources for blue catfish are not publicly available.

Results

The blue catfish genome assembly is 841.86 Mbp in length with excellent continuity (8.6 Mbp contig N50, 28.2 Mbp scaffold N50) and completeness (98.6% Eukaryota and 97.0% Actinopterygii BUSCO). A total of 30,971 protein-coding genes were predicted, of which 21,781 were supported by RNA sequencing evidence. Phylogenomic analyses revealed that it diverged from channel catfish approximately 9 million years ago with 15.7 million fixed nucleotide differences. The within-species single-nucleotide polymorphism (SNP) density is 0.32% between the most aquaculturally important blue catfish strains (D&B and Rio Grande). Gene family analysis discovered significant expansion of immune-related families in the blue catfish lineage, which may contribute to disease resistance in blue catfish.

Conclusions

We reported the first high-quality, chromosome-level assembly of the blue catfish genome, which provides the necessary genomic tool kit for transcriptome and methylome analysis, SNP discovery and marker-assisted selection, gene editing and genome engineering, and reproductive enhancement of the blue catfish and hybrid catfish.

Leafy and weedy seadragon genomes connect genic and repetitive DNA features to the extravagant biology of syngnathid fishes

Seadragons are widely recognized for their derived traits, which include leaf-like appendages and extreme spinal curvature. Efforts to understand the genetic basis of these unique traits and conserve these species and their relatives have been limited by genomic resource gaps. In this paper we present full, annotated genomes of leafy and weedy seadragons, which we use to uncover surprising features of gene family and genome architecture evolution that likely relate to the extravagant phenotypic traits of seadragons and their pipefish and seahorse relatives. These genomes and their analyses are important advances for the study of elaborate vertebrate traits, leveraging this diverse, morphologically exceptional group of fishes.

Seadragons are a remarkable lineage of teleost fishes in the family Syngnathidae, renowned for having evolved male pregnancy. Comprising three known species, seadragons are widely recognized and admired for their fantastical body forms and coloration, and their specific habitat requirements have made them flagship representatives for marine conservation and natural history interests. Until recently, a gap has been the lack of significant genomic resources for seadragons. We have produced gene-annotated, chromosome-scale genome models for the leafy and weedy seadragon to advance investigations of evolutionary innovation and elaboration of morphological traits in seadragons as well as their pipefish and seahorse relatives. We identified several interesting features specific to seadragon genomes, including divergent noncoding regions near a developmental gene important for integumentary outgrowth, a high genome-wide density of repetitive DNA, and recent expansions of transposable elements and a vesicular trafficking gene family. Surprisingly, comparative analyses leveraging the seadragon genomes and additional syngnathid and outgroup genomes revealed striking, syngnathid-specific losses in the family of fibroblast growth factors (FGFs), which likely involve reorganization of highly conserved gene regulatory networks in ways that have not previously been documented in natural populations. The resources presented here serve as important tools for future evolutionary studies of developmental processes in syngnathids and hold value for conservation of the extravagant seadragons and their relatives.

Revisiting the Tigger Transposon Evolution Revealing Extensive Involvement in the Shaping of Mammal Genomes

Simple Summary

Despite the discovery of the Tigger family of pogo transposons in the mammalian genome, the evolution profile of this family is still incomplete. Here, we conducted a systematic evolution analysis for Tigger in nature. The data revealed that Tigger was found in a broad variety of animals, and extensive invasion of Tigger was observed in mammal genomes. Common horizontal transfer events of Tigger elements were observed across different lineages of animals, including mammals, that may have led to their widespread distribution, while parasites and invasive species may have promoted Tigger HT events. Our results also indicate that the activity of Tigger transposons tends to be low in vertebrates; only one mammalian genome and fish genome may harbor active Tigger.

Abstract

The data of this study revealed that Tigger was found in a wide variety of animal genomes, including 180 species from 36 orders of invertebrates and 145 species from 29 orders of vertebrates. An extensive invasion of Tigger was observed in mammals, with a high copy number. Almost 61% of those species contain more than 50 copies of Tigger; however, 46% harbor intact Tigger elements, although the number of these intact elements is very low. Common HT events of Tigger elements were discovered across different lineages of animals, including mammals, that may have led to their widespread distribution, whereas Helogale parvula and arthropods may have aided Tigger HT incidences. The activity of Tigger seems to be low in the kingdom of animals, most copies were truncated in the mammal genomes and lost their transposition activity, and Tigger transposons only display signs of recent and current activities in a few species of animals. The findings suggest that the Tigger family is important in structuring mammal genomes.

The Role of Transposable Elements of the Human Genome in Neuronal Function and Pathology

Transposable elements (TEs) have been extensively studied for decades. In recent years, the introduction of whole-genome and whole-transcriptome approaches, as well as single-cell resolution techniques, provided a breakthrough that uncovered TE involvement in host gene expression regulation underlying multiple normal and pathological processes. Of particular interest is increased TE activity in neuronal tissue, and specifically in the hippocampus, that was repeatedly demonstrated in multiple experiments. On the other hand, numerous neuropathologies are associated with TE dysregulation. Here, we provide a comprehensive review of literature about the role of TEs in neurons published over the last three decades. The first chapter of the present review describes known mechanisms of TE interaction with host genomes in general, with the focus on mammalian and human TEs; the second chapter provides examples of TE exaptation in normal neuronal tissue, including TE involvement in neuronal differentiation and plasticity; and the last chapter lists TE-related neuropathologies. We sought to provide specific molecular mechanisms of TE involvement in neuron-specific processes whenever possible; however, in many cases, only phenomenological reports were available. This underscores the importance of further studies in this area.

Taming, Domestication and Exaptation: Trajectories of Transposable Elements in Genomes

During evolution, several types of sequences pass through genomes. Along with mutations and internal genetic tinkering, they are a useful source of genetic variability for adaptation and evolution. Most of these sequences are acquired by horizontal transfers (HT), but some of them may come from the genomes themselves. If they are not lost or eliminated quickly, they can be tamed, domesticated, or even exapted. Each of these processes results from a series of events, depending on the interactions between these sequences and the host genomes, but also on environmental constraints, through their impact on individuals or population fitness. After a brief reminder of the characteristics of each of these states (taming, domestication, exaptation), the evolutionary trajectories of these new or acquired sequences will be presented and discussed, emphasizing that they are not totally independent insofar as the first can constitute a step towards the second, and the second is another step towards the third.

Prokaryotic and Eukaryotic Horizontal Transfer of Sailor (DD82E), a New Superfamily of IS630-Tc1-Mariner DNA Transposons

Simple Summary

Transposable elements, including DNA transposons, play a significant role in genetic material exchanges between prokaryotes and eukaryotes. Comparative profiling of the evolution pattern of DNA transposons between prokaryotes and eukaryotes may identify potential genetic material exchanges between them and provide insights into the evolutionary history of prokaryotic and eukaryotic genomes. The members of the IS630-Tc1-mariner (ITm) group may represent the most diverse and widely distributed DNA transposons in nature, and the discovery of new members of this group is highly expected based on the increasing availability of genome sequencing data. We discovered a new superfamily (termed Sailor) belonging to the ITm hyperfamily, which differed from the known superfamilies of Tc1/mariner, DDxD/pogo and DD34E/Gambol, regarding phylogenetic position and catalytic domain. Our data revealed that Sailor was distributed in both prokaryotes and eukaryotes and suggested that horizontal transfer (HT) events of Sailor may occur from prokaryotic to eukaryotic genomes. Finally, internal transmissions of Sailor in prokaryotes and eukaryotes were also detected.

Abstract

Here, a new superfamily of IS630-Tc1-mariner (ITm) DNA transposons, termed Sailor, is identified, that is characterized by a DD82E catalytic domain and is distinct from all previously known superfamilies of the ITm group. Phylogenetic analyses revealed that Sailor forms a monophyletic clade with a more intimate link to the clades of Tc1/mariner and DD34E/Gambol. Sailor was detected in both prokaryotes and eukaryotes and invaded a total of 256 species across six kingdoms. Sailor is present in nine species of bacteria, two species of plantae, four species of protozoa, 23 species of Chromista, 12 species of Fungi and 206 species of animals. Moreover, Sailor is extensively distributed in invertebrates (a total of 206 species from six phyla) but is absent in vertebrates. Sailor transposons are 1.38–6.98 kb in total length and encoded transposases of ~676 aa flanked by TIRs with lengths between 18, 1362 and 4 bp (TATA) target-site duplications. Furthermore, our analysis provided strong evidence of Sailor transmissions from prokaryotes to eukaryotes and internal transmissions in both. These data update the classification of the ITm group and will contribute to the understanding of the evolution of ITm transposons and that of their hosts.

Clinical and Molecular Features of Anti-CENP-B Autoantibodies

Centromeric proteins are the foundation for assembling the kinetochore, a macromolecular complex that is essential for accurate chromosome segregation during mitosis. Anti-centromere antibodies (ACAs) are polyclonal autoantibodies targeting centromeric proteins (CENP-A, CENP-B, CENP-C), predominantly CENP-B, and are highly associated with rheumatologic disease (lcSSc/CREST syndrome). CENP-B autoantibodies have also been reported in cancer patients without symptoms of rheumatologic disease. The rise of oncoimmunotherapy stimulates inquiry into how and why anti-CENP-B autoantibodies are formed. In this review, we describe the clinical correlations between anti-CENP-B autoantibodies, rheumatologic disease, and cancer; the molecular features of CENP-B; possible explanations for autoantigenicity; and, finally, a possible mechanism for induction of autoantibody formation.

Diversity and Evolution of pogo and Tc1/mariner Transposons in the Apoidea Genomes

Bees (Apoidea), the largest and most crucial radiation of pollinators, play a vital role in the ecosystem balance. Transposons are widely distributed in nature and are important drivers of species diversity. However, transposons are rarely reported in important pollinators such as bees. Here, we surveyed 37 bee genomesin Apoidea, annotated the pogo and Tc1/mariner transposons in the genome of each species, and performed a phylogenetic analysis and determined their overall distribution. The pogo and Tc1/mariner families showed high diversity and low abundance in the 37 species, and their proportion was significantly higher in solitary bees than in social bees. DD34D/mariner was found to be distributed in almost all species and was found in Apis mellifera, Apis mellifera carnica, Apis mellifera caucasia, and Apis mellifera mellifera, and Euglossa dilemma may still be active. Using horizontal transfer analysis, we found that DD29-30D/Tigger may have experienced horizontal transfer (HT) events. The current study displayed the evolution profiles (including diversity, activity, and abundance) of the pogo and Tc1/mariner transposons across 37 species of Apoidea. Our data revealed their contributions to the genomic variations across these species and facilitated in understanding of the genome evolution of this lineage.

Parallel pathways for recruiting effector proteins determine centromere drive and suppression

Selfish centromere DNA sequences bias their transmission to the egg in female meiosis. Evolutionary theory suggests that centromere proteins evolve to suppress costs of this "centromere drive." In hybrid mouse models with genetically different maternal and paternal centromeres, selfish centromere DNA exploits a kinetochore pathway to recruit microtubule-destabilizing proteins that act as drive effectors. We show that such functional differences are suppressed by a parallel pathway for effector recruitment by heterochromatin, which is similar between centromeres in this system. Disrupting the kinetochore pathway with a divergent allele of CENP-C reduces functional differences between centromeres, whereas disrupting heterochromatin by CENP-B deletion amplifies the differences. Molecular evolution analyses using Murinae genomes identify adaptive evolution in proteins in both pathways. We propose that centromere proteins have recurrently evolved to minimize the kinetochore pathway, which is exploited by selfish DNA, relative to the heterochromatin pathway that equalizes centromeres, while maintaining essential functions.

Characterization and expression pattern of ZB and PS transposons in zebrafish

Despite comprising much of the genome, transposons were once thought of as junk. However, transposons play many roles in the eukaryotic genome, such as providing new proteins as domesticated genes, expressing during germline-soma differentiation, function in DNA rearrangement in the offspring, and so on. We sought to describe the distribution and structural organization of the two autonomous transposons (ZB and PS) in the zebrafish genome and examine their expression patterns in embryos and adult tissues. The intact copy of ZB and PS was queried by BLAST on NCBI and ENSEMBL using default parameters. Of the copies with coverage and identity, more than 90 % were downloaded to do structural analysis. Spatial and temporal expression patterns were detected by qRT-PCR and Whole-mount in situ hybridization (WISH). There are 19 intact copies of ZB, encoding 341 amino acid residues with DD34E catalytic domain and flanked by 201bp TIRs, and seven intact PS copies, containing 425 amino acid residues with DD35D catalytic domain flanked by 28bp TIRs, were detected in the genome of zebrafish respectively. Analysis of genomic insertions indicated that both ZB and PS transposons are prone to be retained in the intron and intergenic regions of the zebrafish genome. The sense and antisense transcripts of ZB and PS were detected during embryonic development stages and exhibited similar expression patterns. The difference is that the sense strand transcript of ZB was explicitly expressed in midbrain-hindbrain boundary (MHB) and otic vesicle (OV), and pharyngeal arches and pharyngeal pouches (PA&PP) at 48 hpf. In adult zebrafish, the expressions of ZB and PS in muscle and brain are much higher than in other tissues. Our study results indicate that ZB and PS transposons may be involved in the embryonic development and regulation of somatic cells of certain adult tissues, such as the brain and muscle.

Characterization of Mariner transposons in seven species of Rhus gall aphids

Transposable elements (TEs), also known as jumping genes, are widely spread in the genomes of insects and play a considerable role in genomic evolution. Mariner/DD34D family belongs to class II transposable elements which is widely spread in the genomes of insects and have considerable role in genomic evolution. Mariner like elements (MLEs) were searched in the genomes of seven species of Rhus gall aphids belonging to six genera. In total, 121 MLEs were detected in the genomes of the seven investigated species of Rhus gall aphids, which showed a wide distribution in both close and distant related species. The sequences of MLEs ranged from 1 to 1.4 kb in length and the structural analysis of the MLEs showed that only five copies were potentially active with intact open reading frame (ORF) and terminal inverted repeats (TIRs). Phylogenetic analysis showed that all the 121 MLE sequences belonged to four subfamilies, i.e., Mauritiana, Drosophila, Vertumana and Irritans, among which Drosophila and Vertumana subfamilies were reported in aphids for the first time. Our present report revealed the diversity and distribution of MLEs in Rhus gall aphid genomes and expanded our understandings on the characterization of transposable elements in aphid genomes, which might be useful as genetic markers and tools and would play an important role in genomic evolution and adaptation of aphids.

The IS630/Tc1/mariner transposons in three ctenophore genomes

Transposable elements (TEs) exert a significant effect on the structure and functioning of the genomes and also serve as a source of the new genes. The study of the TE diversity and evolution in different taxa is indispensable for the fundamental understanding of their roles in the genomes. IS630/Tc1/mariner (ITm) transposable elements represent the most prevalent and diverse group of DNA transposons. In this work, we studied the diversity, evolutionary dynamics and the phylogenetic relationships of the ITm transposons found in three ctenophore species: Mnemiopsis leidyi, Pleurobrachia bachei, Beroe ovata. We identified 29 ITm transposons, seven of which possess the terminal inverted repeats (TIRs) and an intact transposase, and, thus, are, presumably, active. Four other ITm transposons have the features of domesticated TEs. According to the results of the phylogenetic analysis, the ITm transposons of the ctenophores represent five groups - MLE/DD34D, TLE/DD34-38E, mosquito/DD37E, Visiror/DD41D and pogo/DDxD. Pogo/DDxD superfamily turnes out to be the most diverse and prevalent, since it accounts for more than 40% of the TEs identified. The data obtained in this research will fill the gap of knowledge of the diversity and evolution of the ITm transposons in the multicellular genomes and will lay the ground for the study of the TE effects on the evolution of the ctenophores.

Distribution of Merlin in eukaryotes and first report of DNA transposons in kinetoplastid protists

DNA transposons are defined as repeated DNA sequences that can move within the host genome through the action of transposases. The transposon superfamily Merlin was originally found mainly in animal genomes. Here, we describe a global distribution of the Merlin in animals, fungi, plants and protists, reporting for the first time their presence in Rhodophyceae, Metamonada, Discoba and Alveolata. We identified a great variety of potentially active Merlin families, some containing highly imperfect terminal inverted repeats and internal tandem repeats. Merlin-related sequences with no evidence of mobilization capacity were also observed and may be products of domestication. The evolutionary trees support that Merlin is likely an ancient superfamily, with early events of diversification and secondary losses, although repeated re-invasions probably occurred in some groups, which would explain its diversity and discontinuous distribution. We cannot rule out the possibility that the Merlin superfamily is the product of multiple horizontal transfers of related prokaryotic insertion sequences. Moreover, this is the first account of a DNA transposon in kinetoplastid flagellates, with conserved Merlin transposase identified in Bodo saltans and Perkinsela sp., whereas it is absent in trypanosomatids. Based on the level of conservation of the transposase and overlaps of putative open reading frames with Merlin, we propose that in protists it may serve as a raw material for gene emergence.

Divergent evolution profiles of DD37D and DD39D families of Tc1/mariner transposons in eukaryotes

DNA transposons play a significant role in shaping the size and structure of eukaryotic genomes. The Tc1/mariner transposons are the most diverse and widely distributed superfamily of DNA transposons and the structure and distribution of several Tc1/mariner families, such as DD35E/TR, DD36E/IC, DD37E/TRT, and DD41D/VS, have been well studied. Nonetheless, a greater understanding of the structure and diversity of Tc1/mariner transposons will provide insight into the evolutionary history of eukaryotic genomes. Here, we conducted further analysis of DD37D/maT and DD39D (named Guest, GT), which were identified by the specific catalytic domains DD37D and DD39D. Most transposons of the maT family have a total length of approximately 1.3 kb and harbor a single open reading frame encoding a ~ 346 amino acid (range 302-398 aa) transposase protein, flanked by short terminal inverted repeats (TIRs) (13-48 base pairs, bp). In contrast, GTs transposons were longer (2.0-5.8 kb), encoded a transposase protein of ~400 aa (range 140-592 aa), and were flanked by short TIRs (19-41 bp). Several conserved motifs, including two helix-turn-helix (HTH) motifs, a GRPR (GRKR) motif, a nuclear localization sequence, and a DDD domain, were also identified in maT and GT transposases. Phylogenetic analyses of the DDD domain showed that the maT and GT families each belong to a monophyletic clade and appear to be closely related to DD41D/VS and DD34D/mariner. In addition, maTs are mainly distributed in invertebrates (144 species), whereas GTs are mainly distributed in land plants through a small number of GTs are present in Chromista and animals. Sequence identity and phylogenetic analysis revealed that horizontal transfer (HT) events of maT and GT might occur between kingdoms and phyla of eukaryotes; however, pairwise distance comparisons between host genes and transposons indicated that HT events involving maTs might be less frequent between invertebrate species and HT events involving GTs may be less frequent between land plant species. Overall, the DD37D/maT and DD39D/GT families display significantly different distribution and tend to be identified in more ancient evolutionary families. The discovery of intact transposases, perfect TIRs, and target site duplications (TSD) of maTs and GTs illustrates that the DD37D/maT and DD39D/GT families may be active. Together, these findings improve our understanding of the diversity of Tc1/mariner transposons and their impact on eukaryotic genome evolution.

Intruder (DD38E), a recently evolved sibling family of DD34E/Tc1 transposons in animals

BACKGROUND

A family of Tc1/mariner transposons with a characteristic DD38E triad of catalytic amino acid residues, named Intruder (IT), was previously discovered in sturgeon genomes, but their evolutionary landscapes remain largely unknown.

RESULTS

Here, we comprehensively investigated the evolutionary profiles of ITs, and evaluated their cut-and-paste activities in cells. ITs exhibited a narrow taxonomic distribution pattern in the animal kingdom, with invasions into two invertebrate phyla (Arthropoda and Cnidaria) and three vertebrate lineages (Actinopterygii, Agnatha, and Anura): very similar to that of the DD36E/IC family. Some animal orders and species seem to be more hospitable to Tc1/mariner transposons, one order of Amphibia and seven Actinopterygian orders are the most common orders with horizontal transfer events and have been invaded by all four families (DD38E/IT, DD35E/TR, DD36E/IC and DD37E/TRT) of Tc1/mariner transposons, and eight Actinopterygii species were identified as the major hosts of these families. Intact ITs have a total length of 1.5-1.7 kb containing a transposase gene flanked by terminal inverted repeats (TIRs). The phylogenetic tree and sequence identity showed that IT transposases were most closely related to DD34E/Tc1. ITs have been involved in multiple events of horizontal transfer in vertebrates and have invaded most lineages recently (< 5 million years ago) based on insertion age analysis. Accordingly, ITs presented high average sequence identity (86-95%) across most vertebrate species, suggesting that some are putatively active. ITs can transpose in human HeLa cells, and the transposition efficiency of consensus TIRs was higher than that of the TIRs of natural isolates.

CONCLUSIONS

We conclude that DD38E/IT originated from DD34E/Tc1 and can be detected in two invertebrate phyla (Arthropoda and Cnidaria), and in three vertebrate lineages (Actinopterygii, Agnatha and Anura). IT has experienced multiple HT events in animals, dominated by recent amplifications in most species and has high identity among vertebrate taxa. Our reconstructed IT transposon vector designed according to the sequence from the "cat" genome showed high cut-and-paste activity. The data suggest that IT has been acquired recently and is active in many species. This study is meaningful for understanding the evolution of the Tc1/mariner superfamily members and their hosts.