Target Site Duplications Research Articles

Creating Meiotic Recombination-Regulating DNA Sites by SpEDIT in Fission Yeast Reveals Inefficiencies, Target-Site Duplications, and Ectopic Insertions.

Recombination hotspot-activating DNA sites (e.g., M26, CCAAT, Oligo-C) and their binding proteins (e.g., Atf1-Pcr1 heterodimer; Php2-Php3-Php5 complex, Rst2, Prdm9) regulate the distribution of Spo11 (Rec12)-initiated meiotic recombination. We sought to create 14 different candidate regulatory DNA sites via bp substitutions in the ade6 gene of Schizosaccharomyces pombe. We used a fission yeast-optimized CRISPR-Cas9 system (SpEDIT) and 196 bp-long dsDNA templates with centrally located bp substitutions designed to ablate the genomic PAM site, create specific 15 bp-long DNA sequences, and introduce a stop codon. After co-transformation with a plasmid that encoded both the guide RNA and Cas9 enzyme, about one-third of colonies had a phenotype diagnostic for DNA sequence changes at ade6. PCR diagnostics and DNA sequencing revealed a diverse collection of alterations at the target locus, including: (A) complete or (B) partial template-directed substitutions; (C) non-homologous end joinings; (D) duplications; (E) bp mutations, and (F) insertions of ectopic DNA. We concluded that SpEDIT can be used successfully to generate a diverse collection of DNA sequence elements within a reporter gene of interest. However, its utility is complicated by low efficiency, incomplete template-directed repair events, and undesired alterations to the target locus.

The identification of retro-DNAs in primate genomes as DNA transposons mobilizing via retrotransposition.

Mobile elements (MEs) constitute a major portion of the genome in primates and other higher eukaryotes, and they play important role in genome evolution and gene function. MEs can be divided into two fundamentally different classes: DNA transposons which transpose in the genome in a "cut-and-paste" style, and retrotransposons which propagate in a "copy-and-paste" fashion via a process involving transcription and reverse-transcription. In primate genomes, DNA transposons are mostly dead, while many retrotransposons are still highly active. We report here the identification of a unique group of MEs, which we call "retro-DNAs", for their combined characteristics of these two fundamentally different ME classes. A comparative computational genomic approach was used to analyze the reference genome sequences of 10 primate species consisting of five apes, four monkeys, and marmoset. From our analysis, we identified a total of 1,750 retro-DNAs, representing 748 unique insertion events in the genomes of ten primate species including human. These retro-DNAs contain sequences of DNA transposons but lack the terminal inverted repeats (TIRs), the hallmark of DNA transposons. Instead, they show characteristics of retrotransposons, such as polyA tails, longer target-site duplications (TSDs), and the "TT/AAAA" insertion site motif, suggesting the use of the L1-based target- primed reverse transcription (TPRT) mechanism. At least 40% of these retro-DNAs locate into genic regions, presenting potentials for impacting gene function. More interestingly, some retro-DNAs, as well as their parent sites, show certain levels of expression, suggesting that they have the potential to create more retro-DNA copies in the present primate genomes. Although small in number, the identification of these retro-DNAs reveals a new mean for propagating DNA transposons in primate genomes without active canonical DNA transposon activity. Our data also suggest that the TPRT machinery may transpose a wider variety of DNA sequences in the genomes.

The identification of retro-DNAs in primate genomes as DNA transposons mobilizing via retrotransposition

Background Mobile elements (MEs) constitute a major portion of the genome in primates and other higher eukaryotes, and they play important role in genome evolution and gene function. MEs can be divided into two fundamentally different classes: DNA transposons which transpose in the genome in a “cut-and-paste” style, and retrotransposons which propagate in a “copy-and-paste” fashion via a process involving transcription and reverse-transcription. In primate genomes, DNA transposons are mostly dead, while many retrotransposons are still highly active. We report here the identification of a new type of MEs, which we call “retro-DNAs”, for their combined characteristics of these two fundamentally different ME classes. Methods A comparative computational genomic approach was used to analyze the reference genome sequences of 10 primate species consisting of five apes, four monkeys, and marmoset. Results From our analysis, we identified a total of 1,750 retro-DNAs, representing 748 unique insertion events in the genomes of ten primate species including human. These retro-DNAs contain sequences of DNA transposons but lack the terminal inverted repeats (TIRs), the hallmark of DNA transposons. Instead, they show characteristics of retrotransposons, such as polyA tails, longer target-site duplications (TSDs), and the “TT/AAAA” insertion site motif, suggesting the use of the L1-based target-primed reverse transcription (TPRT) mechanism. At least 40% of these retro-DNAs locate into genic regions, presenting potentials for impacting gene function. More interestingly, some retro-DNAs, as well as their parent sites, show certain levels of expression, suggesting that they have the potential to create more retro-DNA copies in the present primate genomes. Conclusions Although small in number, the identification of these retro-DNAs reveals a new mechanism for propagating DNA transposons in primate genomes without active canonical DNA transposon activity. Our data also suggest that the TPRT machinery may transpose a wider variety of DNA sequences in the genomes.

Insights into the Evolution of IncR Plasmids Found in the Southern European Clone of the Monophasic Variant of Salmonella enterica Serovar Typhimurium.

Salmonella enterica subspecies enterica serovar 4,[5],12:i:- is a monophasic variant of S. Typhimurium which has emerged as a world-wide distributed pathogen in the last decades. Several clones have been identified within this variant, the European clone, the Spanish clone, the Southern European clone and the U.S./American clone. The present study focused on isolates of the Southern European clone that were obtained from clinical samples at Spanish hospitals. The selected isolates were multidrug resistant, with most resistance genes residing on IncR plasmids that also carried virulence genes. These plasmids had a mosaic structure, comprising a highly reduced IncR backbone, which has acquired a large amount of exogenous DNA mostly derived from pSLT and IncI1-I(alfa) plasmids. Although composed of approximately the same elements, the investigated plasmids displayed a high diversity, consistent with active evolution driven by a wealth of mobile genetic elements. They comprise multiple intact or truncated insertion sequences, transposons, pseudo-compound transposons and integrons. Particularly relevant was the role of IS26 (with six to nine copies per plasmid) in generating insertions, deletions and inversions, with many of the rearrangements uncovered by tracking the patterns of eight bp target site duplications. Most of the resistance genes detected in the analyzed isolates have been previously associated with the Southern European clone. However, erm(B), lnu(G) and blaTEM-1B are novel, with the last two carried by a second resistance plasmid found in one of the IncR-positive isolates. Thus, evolution of resistance in the Southern European clone is not only mediated by diversification of the IncR plasmids, but also through acquisition of additional plasmids. All isolates investigated in the present study have the large deletion affecting the fljBA region previously found to justify the monophasic phenotype in the Southern European and U.S./American clones. An SNP-based phylogenetic analysis revealed the close relationship amongst our isolates, and support that those sharing the large fljBA deletion could be more heterogeneous than previously anticipated.

Aminoglycoside resistance genes in early members of the Acinetobacter baumannii ST78A (SMAL, Italian clone) reside in an IS26-bounded island in the chromosome.

The Acinetobacter baumannii isolate called SMAL, previously used to determine the structures of capsular polysaccharide and lipooligosaccharide, was recovered in Pavia, Italy in 2002 among the collection of aminoglycoside-resistant isolates designated as SMAL type. This type was later called the Italian clone, then ST78. ST78 isolates are now widely distributed. To establish the resistance gene complement and the location and structure of acquired resistance regions in early members of the Italian/ST78 clone. The draft genome of SMAL2002 was assembled from Illumina MiSeq reads. Contigs containing resistance genes were joined and located in the chromosome using PCR with custom primers. The resistance profile was determined using disc diffusion. SMAL2002 is an ST78A isolate and includes three aminoglycoside resistance genes, aadB (gentamicin, kanamycin, tobramycin) aphA1 (kanamycin, neomycin) and aac(6')-Ian (amikacin, kanamycin, tobramycin). The aadB gene cassette is incorporated at a secondary site in a relative of the aphA1-containing, IS26-bounded pseudo-compound transposon, PTn6020. The aac(6')-Ian gene is in an adjacent IS26-bounded structure that includes sul2 (sulphonamide) and floR (florfenicol) resistance genes. The two pseudo-compound transposons overlap and are in the chromosomal hutU gene flanked by an 8 bp target site duplication. Although aac(6')-Ian was not noticed previously, the same genes and structures were found in several available draft genomes of early ST78A isolates. This study highlights the importance of correlating resistance profiles with resistance gene content. The location of acquired resistance genes in the SMAL2002 chromosome represents the original location in the ST78A lineage of ST78.

Role of amino acid 159 in carbapenem and temocillin hydrolysis of OXA-933, a novel OXA-48 variant.

OXA-48 has rapidly disseminated worldwide and become one of the most common carbapenemases in many countries with more than 45 variants reported with, in some cases, significant differences in their hydrolysis profiles. The R214 residue, located in the ß5-ß6 loop, is crucial for the carbapenemase activity, as it stabilizes carbapenems in the active site and maintains the shape of the active site through interactions with D159. In this study, we have characterized a novel variant of OXA-48, OXA-933 with a single D159N change. To evaluate the importance of this residue, point mutations were generated (D159A, D159G, D159K, and D159W), kinetic parameters of OXA-933, OXA-48 D159G, and OXA-48 D159K were determined and compared to those of OXA-48 and OXA-244. The blaOXA-933 gene was borne on Tn2208, a 2,696-bp composite transposon made of two IS1 elements surrounded by 9 bp target site duplications and inserted into a non-self-transmissible plasmid pOXA-933 of 7,872 bp in size. Minimal inhibitory concentration values of E. coli expressing the blaOXA-933 gene or of its point mutant derivatives were lower for carbapenems (except for D159G) as compared to those expressing the blaOXA-48 gene. Steady-state kinetic parameters revealed lower catalytic efficiencies for expanded spectrum cephalosporins and carbapenems. A detailed structural analysis confirmed the crucial role of D159 in shaping the active site of OXA-48 enzymes by interacting with R214. Our work further illustrates the remarkable propensity of OXA-48-like carbapenemases to evolve through mutations at positions outside the β5-β6 loop, but interacting with key residues of it.

Helenus and Ajax, Two Groups of Non-Autonomous LTR Retrotransposons, Represent a New Type of Small RNA Gene-Derived Mobile Elements.

Terminal repeat retrotransposons in miniature (TRIMs) are short non-autonomous long terminal repeat (LTR) retrotransposons found from various eukaryotes. Cassandra is a unique TRIM lineage which contains a 5S rRNA-derived sequence in its LTRs. Here, two new groups of TRIMs, designated Helenus and Ajax, are reported based on bioinformatics analysis and the usage of Repbase. Helenus is found from fungi, animals, and plants, and its LTRs contain a tRNA-like sequence. It includes two LTRs and between them, a primer-binding site (PBS) and polypurine tract (PPT) exist. Fungal and plant Helenus generate 5 bp target site duplications (TSDs) upon integration, while animal Helenus generates 4 bp TSDs. Ajax includes a 5S rRNA-derived sequence in its LTR and is found from two nemertean genomes. Ajax generates 5 bp TSDs upon integration. These results suggest that despite their unique promoters, Helenus and Ajax are TRIMs whose transposition is dependent on autonomous LTR retrotransposon. These TRIMs can originate through an insertion of SINE in an LTR of TRIM. The discovery of Helenus and Ajax suggests the presence of TRIMs with a promoter for RNA polymerase III derived from a small RNA gene, which is here collectively termed TRIMp3.

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.

The dynamic evolution and IS26-mediated interspecies transfer of a blaNDM-1-bearing fusion plasmid leading to a hypervirulent carbapenem-resistant Klebsiella pneumoniae strain harbouring blaKPC-2 in a single patient

To characterize the evolution and interspecies transfer of plasmids between Klebsiella pneumoniae and Escherichia coli within a single patient. Minimum inhibitory concentrations were measured using broth microdilution assays. Conjugation assays, string tests, and Galleria mellonella infection model experiments were also conducted. Whole-genome sequencing was performed on the Illumina and Nanopore platforms. Antimicrobial resistance determinants, insertion sequences, and virulence factors were identified using ABRicate/ResFinder database, ISFinder, and virulence factor database. Wzi and capsular polysaccharide (KL) were typed using Kleborate and Kaptive. Multi-locus sequence typing (MLST), replicon typing, and single nucleotide polymorphism analyses were conducted using the BacWGSTdb server. The carbapenem-resistant K. pneumoniae 2111KP was characterized as ST11, wzi64, and KL64, with a positive string test result and a relatively high virulence phenotype. Analysis of the 2111KP genome revealed that blaNDM-1 was located in a 268,400-bp IncFIB/IncHI1B/IncX3 conjugative plasmid (p2111KP-1), regulated by IS26, IS5, and ISKox3. p2111KP-1 was also a rmpA2-associated virulence plasmid with an iutA-iucABCD gene cluster and a IS26-mediated multidrug-resistant fusion plasmid, which contained 8-bp (AGCTGCAC or GGCCTTTG) target site duplications. Segments flanked by IS26 of p2111KP-1 were 99.99% identical to a 49,016-bp E. coli plasmid. This study provided direct evidence of plasmid fusion via IS26 between two different bacterial species within one patient and revealed the process by which genetic elements conferring carbapenem resistance and virulence were simultaneously transferred between these species. It highlights the need for strategic antibiotic use and rigorous monitoring to prevent the plasmid-mediated fusion and transmission of drug-resistance/virulence factors.

Emergence of IS26-mediated pLVPK-like virulence and NDM-1 conjugative fusion plasmid in hypervirulent carbapenem-resistant Klebsiella pneumoniae.

Hypervirulent carbapenem-resistant Klebsiella pneumoniae (hv-CRKP) has been widely reported and poses a global threat. However, the comprehensive genetic structure of ST11-KL64 hv-CRKP and the possible evolutionary mechanisms from a genetic structure perspective of this high-risk clone remain unclear. Here, a blaKPC-2-blaNDM-1-positive ST11-KL64 hv-CRKP isolate was obtained from a human bloodstream infection (BSI). Whole-genome sequencing and bioinformatics analyses revealed that it contained a fusion plasmid, pKPTCM2-1. pKPTCM2-1 is a conjugative plasmid composed of an oriT-positive pLVPK-like virulence plasmid and a type IV secretion system-produced blaNDM-1-bearing IncX3 plasmid mediated by IS26-based co-integration. This progress generated 8-bp target site duplications (TGAAAACC) on both sides. The fusion plasmid possessed self-transferability and could be transferred to blaKPC-2-harboring ST11-KL64 CRKP to form the ST11-KL64 hv-CRKP clone. The pLVPK-like-positive ST11-KL64 strain exhibited virulence levels similar to those of the typical hypervirulent K. pneumoniae NTUH-2044. The mutation, Tet(A) (A276S), which was believed to lead to tigecycline resistance was observed. Overall, this high-risk clone has emerged as a tremendous threat in fatal BSIs and thus, targeted surveillance is an urgent need to contain the hv-CRKP clones.

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.

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.

Active in vivo translocation of the Methanosarcina mazei Gö1 Casposon.

Casposons are transposable elements containing the CRISPR associated gene Cas1solo. Identified in many archaeal genomes, casposons are discussed as the origin of CRISPR-Cas systems due to their proposed Cas1solo-dependent translocation. However, apart from bioinformatic approaches and the demonstration of Cas1solo integrase and endonuclease activity in vitro, casposon transposition has not yet been shown in vivo. Here, we report on active casposon translocations in Methanosarcina mazei Gö1 using two independent experimental approaches. First, mini-casposons, consisting of a R6Kγ origin and two antibiotic resistance cassettes, flanked by target site duplications (TSDs) and terminal inverted repeats (TIRs), were generated, and shown to actively translocate from a suicide plasmid and integrate into the chromosomal MetMaz-C1 TSD IS1a. Second, casposon excision activity was confirmed in a long-term evolution experiment using a Cas1solo overexpression strain in comparison to an empty vector control under four different treatments (native, high temperature, high salt, mitomycin C) to study stress-induced translocation. Analysis of genomic DNA using a nested qPCR approach provided clear evidence of casposon activity in single cells and revealed significantly different casposon excision frequencies between treatments and strains. Our results, providing the first experimental evidence for in vivo casposon activity are summarized in a modified hypothetical translocation model.

Evolution of piggyBac Transposons in Apoidea

In this study, we investigated the presence of piggyBac (PB) transposons in 44 bee genomes from the Apoidea order, which is a superfamily within the Hymenoptera, which includes a large number of bee species crucial for pollination. We annotated the PB transposons in these 44 bee genomes and examined their evolution profiles, including structural characteristics, distribution, diversity, activity, and abundance. The mined PB transposons were divided into three clades, with uneven distribution in each genus of PB transposons in Apoidea. The complete PB transposons we discovered are around 2.23-3.52 kb in length and encode transposases of approximately 580 aa, with terminal inverted repeats (TIRs) of about 14 bp and 4 bp (TTAA) target-site duplications. Long TIRs (200 bp, 201 bp, and 493 bp) were also detected in some species of bees. The DDD domains of the three transposon types were more conserved, while the other protein domains were less conserved. Generally, most PB transposons showed low abundance in the genomes of Apoidea. Divergent evolution dynamics of PB were observed in the genomes of Apoidea. PB transposons in some identified species were relatively young, whiles others were older and with some either active or inactive. In addition, multiple invasions of PB were also detected in some genomes of Apoidea. Our findings highlight the contribution of PB transposons to genomic variation in these species and suggest their potential as candidates for future gene transfer tools.

Characterization of miniature inverted repeat transposable elements inserted in the CitRWP gene controlling nucellar embryony and development of molecular markers for reliable genotyping of CitRWP in Citrus species

Nucellar embryony, a form of apomixis, is one of the most important traits in citrus breeding. Transposition of a miniature inverted repeat transposable element (MITE) in the promoter of CitRWP gene has been previously reported as the causal mutation. However, features of the MITE such as sequences of terminal inverted repeats (TIRs) and target site duplications (TSDs) remain unclear. An MITE family, DTM10, showing the highest homology was identified from 110 citrus MITE families. At least 80 DTM10 elements were identified from each of 10 genome sequences of Citrinae species. Analysis of 114 sweet orange elements showed that DTM10 contained 56-bp TIRs and mostly 9-bp TSDs. Phylogenetic analyses of DTM10 elements identified from five Citrinae species showed that MITEs in CitRWP might have transposed in a common ancestor of Citrus species. Two inverted MITEs were identified from three mutant alleles analyzed. In addition, three tandem MITEs were identified in the mutant allele of Fortunella CitRWP. A strong DNA secondary structure formed between MITEs was predicted. Due to a strong stem-loop structure, a phenomenon of polymerase jumping was observed during PCR amplification. To develop reliable molecular markers, sequences of mutant and wild-type CitRWP alleles were obtained from diverse citrus accessions. Phylogenetic analysis of CitRWP alleles showed that mutant alleles might have originated from a common ancestor of Citrus species. Based on MITE insertion and a single nucleotide polymorphism in the genic region, three complementary molecular markers were developed. Genotyping results of 241 citrus accessions were found to be identical among three molecular markers.

The identification of retro-DNAs in primate genomes as DNA transposons mobilizing via retrotransposition

Background: Mobile elements (MEs) constitute a major portion of the genome in primates and other higher eukaryotes, and they play important role in genome evolution and gene function. MEs can be divided into two fundamentally different classes: DNA transposons which transpose in the genome in a “cut-and-paste” style, and retrotransposons which propagate in a “copy-and-paste” fashion via a process involving transcription and reverse-transcription. In primate genomes, DNA transposons are mostly dead, while many retrotransposons are still highly active. We report here the identification of a new type of MEs, which we call “retro-DNAs”, for their combined characteristics of these two fundamentally different ME classes. Methods: A comparative computational genomic approach was used to analyze the reference genome sequences of 10 primate species consisting of five apes, four monkeys, and marmoset. Results: From our analysis, we identified a total of 1,750 retro-DNAs, representing 748 unique insertion events in the genomes of ten primate species including human. These retro-DNAs contain sequences of DNA transposons but lack the terminal inverted repeats (TIRs), the hallmark of DNA transposons. Instead, they show characteristics of retrotransposons, such as polyA tails, longer target-site duplications (TSDs), and the “TT/AAAA” insertion site motif, suggesting the use of the L1-based target-primed reverse transcription (TPRT) mechanism. At least 40% of these retro-DNAs locate into genic regions, presenting potentials for impacting gene function. More interestingly, some retro-DNAs, as well as their parent sites, show certain levels of expression, suggesting that they have the potential to create more retro-DNA copies in the present primate genomes. Conclusions: Although small in number, the identification of these retro-DNAs reveals a new mechanism for propagating DNA transposons in primate genomes without active canonical DNA transposon activity. Our data also suggest that the TPRT machinery may transpose a wider variety of DNA sequences in the genomes.

Insertion sequences related to ISAjo2 target pdif and dif sites and belong to a new IS family, the IS1202 family.

Several insertion sequences (IS) found in various Acinetobacter species exhibit target specificity. They are found, in the same orientation, 5 bp from the XerC binding site of the pdif sites associated with dif modules in Acinetobacter plasmids, and searches revealed they are also found near chromosomal dif sites of Acinetobacter species. These IS are 1.5 kb long, bounded by 24-26 bp imperfect terminal inverted repeats (TIRs) and encode a large transposase of 441-457 aa. They generate 5 bp target site duplications (TSDs). Structural predictions of the ISAjo2 transposase, TnpAjo2, modelled on TnsB of Tn7 revealed two N-terminal HTH domains followed by an RNaseH fold (DDE domain), a β barrel and a C-terminal domain. Similar to Tn7, the outer IS ends are 5'-TGT and ACA-3', and an additional Tnp binding site, corresponding to the internal portion of the IR, is found near each end. However, the Acinetobacter IS do not encode further proteins related to those required by Tn7 for targeted transposition, and the transposase may interact directly with XerC bound to a dif-like site. We propose that these IS, currently in the IS1202 group in the not characterized yet (NCY) category in ISFinder, are part of a distinct IS1202 family. Other IS listed as in the IS1202 group encode transposases related to TnpAjo2 (25-56 % amino acid identity) and have similar TIRs but fall into three groups based on the TSD length (3-5, >15, 0 bp). Those with 3-5 bp TSDs may also target dif-like sites but targets were not found for the other groups.

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.

Structural Analysis of the A Mating Type Locus and Development of the Mating Type Marker of Agaricus bisporus var. bisporus

Karyotyping in Agaricus bisporus is crucial for both the isolation of homokaryotic strains and the confirmation of dikaryon establishment. For the verification of the karyotype, the A mating type loci of two homokaryotic strains, H39 and H97, were analyzed through comparative sequence analysis. The two loci showed major differences in two sequence regions designated as Region 1 and Region 2. H97 had a putative DNA transposon in Region 1 that had target site duplications (TSDs), terminal inverted repeats (TIRs), and a loop sequence, in contrast to H39, which only had the insertional target sequence. Homologous sequences of the transposon were discovered in the two different chromosomes of H97 and in one of H39, all of which have different TSDs but share high sequence homology in TIR. Region 2 shared three consensus sequences between H97 and H39. However, it was only from H97 that a large insertional sequence of unknown origin was discovered between the first and second consensus sequences. The difference in length in Region 1, employed for the verification of the A mating type, resulted in the successful verification of mating types in the heterokaryotic and homokaryotic strains. This length difference enables the discrimination between homo- and heterokaryotic spores by PCR. The present study suggests that the A mating type locus in A. bisporus H97 has evolved through transposon insertion, allowing the discrimination of the mating type, and thus the nuclear type, between A. bisporus H97 and H39.

The Transcriptomic Response of the Boll Weevil, Anthonomus grandis grandis Boheman (Coleoptera: Curculionidae), following Exposure to the Organophosphate Insecticide Malathion.

Insecticide tolerance and resistance have evolved countless times in insect systems. Molecular drivers of resistance include mutations in the insecticide target site and/or gene duplication, and increased gene expression of detoxification enzymes. The boll weevil, Anthonomus grandis grandis Boheman (Coleoptera: Curculionidae), is a pest of commercial cotton and has developed resistance in the field to several insecticides; however, the current organophosphate insecticide, malathion, used by USA eradication programs remains effective despite its long-term use. Here, we present findings from an RNA-seq experiment documenting gene expression post-exposure to field-relevant concentrations of malathion, which was used to provide insight on the boll weevil's continued susceptibility to this insecticide. Additionally, we incorporated a large collection of boll weevil whole-genome resequencing data from nearly 200 individuals collected from three geographically distinct areas to determine SNP allele frequency of the malathion target site, as a proxy for directional selection in response to malathion exposure. No evidence was found in the gene expression data or SNP data consistent with a mechanism of enhanced tolerance or resistance adaptation to malathion in the boll weevil. Although this suggests continued effectiveness of malathion in the field, we identified important temporal and qualitative differences in gene expression between weevils exposed to two different concentrations of malathion. We also identified several tandem isoforms of the detoxifying esterase B1 and glutathione S-transferases, which are putatively associated with organophosphate resistance.