Characterization and functional annotation of nested transposable elements in eukaryotic genomes

Modeling early germline immunization after horizontal transfer of transposable elements reveals internal piRNA cluster heterogeneity

BACKGROUND

A fraction of all genomes is composed of transposable elements (TEs) whose mobility needs to be carefully controlled. In gonads, TE activity is repressed by PIWI-interacting RNAs (piRNAs), a class of small RNAs synthesized by heterochromatic loci enriched in TE fragments, called piRNA clusters. Maintenance of active piRNA clusters across generations is secured by maternal piRNA inheritance providing the memory for TE repression. On rare occasions, genomes encounter horizontal transfer (HT) of new TEs with no piRNA targeting them, threatening the host genome integrity. Naïve genomes can eventually start to produce new piRNAs against these genomic invaders, but the timing of their emergence remains elusive.

RESULTS

Using a set of TE-derived transgenes inserted in different germline piRNA clusters and functional assays, we have modeled a TE HT in Drosophila melanogaster. We have found that the complete co-option of these transgenes by a germline piRNA cluster can occur within four generations associated with the production of new piRNAs all along the transgenes and the germline silencing of piRNA sensors. Synthesis of new transgenic TE piRNAs is linked to piRNA cluster transcription dependent on Moonshiner and heterochromatin mark deposition that propagates more efficiently on short sequences. Moreover, we found that sequences located within piRNA clusters can have different piRNA profiles and can influence transcript accumulation of nearby sequences.

CONCLUSIONS

Our study reveals that genetic and epigenetic properties, such as transcription, piRNA profiles, heterochromatin, and conversion efficiency along piRNA clusters, could be heterogeneous depending on the sequences that compose them. These findings suggest that the capacity of transcriptional signal erasure induced by the chromatin complex specific of the piRNA cluster can be incomplete through the piRNA cluster loci. Finally, these results have revealed an unexpected level of complexity that highlights a new magnitude of piRNA cluster plasticity fundamental for the maintenance of genome integrity.

Transposable elements and their role in aging

Transposable elements (TEs) are an important part of eukaryotic genomes. The role of somatic transposition in aging, carcinogenesis, and other age-related diseases has been determined. This review discusses the fundamental properties of TEs and their complex interactions with cellular processes, which are crucial for understanding the diverse effects of their activity on the genetics and epigenetics of the organism. The interactions of TEs with recombination, replication, repair, and chromosomal regulation; the ability of TEs to maintain a balance between their own activity and repression, the involvement of TEs in the creation of new or alternative genes, the expression of coding/non-coding RNA, and the role in DNA damage and modification of regulatory networks are reviewed. The contribution of the derepressed TEs to age-dependent effects in individual cells/tissues in different organisms was assessed. Conflicting information about TE activity under stress as well as theories of aging mechanisms related to TEs is discussed. On the one hand, transposition activity in response to stressors can lead to organisms acquiring adaptive innovations of great importance for evolution at the population level. On the other hand, the TE expression can cause decreased longevity and stress tolerance at the individual level. The specific features of TE effects on aging processes in germline and soma and the ways of their regulation in cells are highlighted. Recent results considering somatic mutations in normal human and animal tissues are indicated, with the emphasis on their possible functional consequences. In the context of aging, the correlation between somatic TE activation and age-related changes in the number of proteins required for heterochromatin maintenance and longevity regulation was analyzed. One of the original features of this review is a discussion of not only effects based on the TEs insertions and the associated consequences for the germline cell dynamics and somatic genome, but also the differences between transposon- and retrotransposon-mediated structural genome changes and possible phenotypic characteristics associated with aging and various age-related pathologies. Based on the analysis of published data, a hypothesis about the influence of the species-specific features of number, composition, and distribution of TEs on aging dynamics of different animal genomes was formulated.

Unravelling complex transposable elements surrounding blaGES-16 in a Pseudomonas aeruginosa ExoU strain

OBJECTIVES

We characterised the complex surrounding regions of bla_GES-16 in a Pseudomonas aeruginosa exoU+ strain (P-10.226) in Brazil.

METHODS

Species identification was performed by MALDI-TOF MS, and the antimicrobial susceptibility profile was determined by broth microdilution based on European Committee on Antimicrobial Susceptibility Testing (EUCAST) breakpoints. The whole genome sequencing (WGS) of P-10.226 strain was performed using both short-read paired-end sequencing on the Illumina MiSeq platform as well as the long-read Oxford Nanopore MinION.

RESULTS

WGS analysis showed that P-10.226 carried bla_GES-16, which was found as a gene cassette inserted into a novel class I integron, In1992 (aadB-bla_OXA-56-bla_GES-16-aadB-aadA6c), whose 3'-CS was truncated by a nested transposable element, IS5564::ISPa157. The structure was even more complex since IS6100-ΔIS6100 structure and a TnAs2-like harbouring the operon merRTPADE was found downstream In1992. Fragments of TnAs3 harbouring 25-bp imperfect inverted repeats were identified bordering the intl1 of In1992 and also flanking IS6100-ΔIS6100, which might be genetic marks of its previous presence in the genome. Interestingly, In1992 also shows a distinct cassette array from In581 (bla_GES-16-dfrA22-aacA27-aadA1), which was previously reported in Serratia marcescens strains recovered in Brazil. Finally, exoU gene, which encodes a potent cytotoxin of type III secretion systems (T3SS) effector proteins from P. aeruginosa and is associated to severe infections, was also detected.

CONCLUSION

We described the novel In1992 carrying bla_GES-16 surrounded by complex transposition events in a XDR P. aeruginosa strain. The identification of many sets of direct repeats adjacent to TnAs3 fragments indicates a major past of transposition events that shaped the current genetic environment of In1992.

Automatic curation of LTR retrotransposon libraries from plant genomes through machine learning

Transposable elements are mobile sequences that can move and insert themselves into chromosomes, activating under internal or external stimuli, giving the organism the ability to adapt to the environment. Annotating transposable elements in genomic data is currently considered a crucial task to understand key aspects of organisms such as phenotype variability, species evolution, and genome size, among others. Because of the way they replicate, LTR retrotransposons are the most common transposable elements in plants, accounting in some cases for up to 80% of all DNA information. To annotate these elements, a reference library is usually created, a curation process is performed, eliminating TE fragments and false positives and then annotated in the genome using the homology method. However, the curation process can take weeks, requires extensive manual work and the execution of multiple time-consuming bioinformatics software. Here, we propose a machine learning-based approach to perform this process automatically on plant genomes, obtaining up to 91.18% F1-score. This approach was tested with four plant species, obtaining up to 93.6% F1-score (Oryza granulata) in only 22.61 s, where bioinformatics methods took approximately 6 h. This acceleration demonstrates that the ML-based approach is efficient and could be used in massive sequencing projects.

Software evaluation for de novo detection of transposons

Transposable elements (TEs) are major genomic components in most eukaryotic genomes and play an important role in genome evolution. However, despite their relevance the identification of TEs is not an easy task and a number of tools were developed to tackle this problem. To better understand how they perform, we tested several widely used tools for de novo TE detection and compared their performance on both simulated data and well curated genomic sequences. As expected, tools that build TE-models performed better than k-mer counting ones, with RepeatModeler beating competitors in most datasets. However, there is a tendency for most tools to identify TE-regions in a fragmented manner and it is also frequent that small TEs or fragmented TEs are not detected. Consequently, the identification of TEs is still a challenging endeavor and it requires a significant manual curation by an experienced expert. The results will be helpful for identifying common issues associated with TE-annotation and for evaluating how comparable are the results obtained with different tools.

A large transposable element mediates metal resistance in the fungus Paecilomyces variotii

The horizontal transfer of large gene clusters by mobile elements is a key driver of prokaryotic adaptation in response to environmental stresses. Eukaryotic microbes face similar stresses; however, a parallel role for mobile elements has not been established. A stress faced by many microorganisms is toxic metal ions in their environment. In fungi, identified mechanisms for protection against metals generally rely on genes that are dispersed within an organism's genome. Here, we discover a large (∼85 kb) region that confers tolerance to five metal/metalloid ions (arsenate, cadmium, copper, lead, and zinc) in the genomes of some, but not all, strains of a fungus, Paecilomyces variotii. We name this region HEPHAESTUS (Hφ) and present evidence that it is mobile within the P. variotii genome with features characteristic of a transposable element. HEPHAESTUS contains the greatest complement of host-beneficial genes carried by a transposable element in eukaryotes, suggesting that eukaryotic transposable elements might play a role analogous to bacteria in the horizontal transfer of large regions of host-beneficial DNA. Genes within HEPHAESTUS responsible for individual metal tolerances include those encoding a P-type ATPase transporter-PcaA-required for cadmium and lead tolerance, a transporter-ZrcA-providing tolerance to zinc, and a multicopper oxidase-McoA-conferring tolerance to copper. In addition, a subregion of Hφ confers tolerance to arsenate. The genome sequences of other fungi in the Eurotiales contain further examples of HEPHAESTUS, suggesting that it is responsible for independently assembling tolerance to a diverse array of ions, including chromium, mercury, and sodium.

TE-greedy-nester: structure-based detection of LTR retrotransposons and their nesting

Motivation

Transposable elements (TEs) in eukaryotes often get inserted into one another, forming sequences that become a complex mixture of full-length elements and their fragments. The reconstruction of full-length elements and the order in which they have been inserted is important for genome and transposon evolution studies. However, the accumulation of mutations and genome rearrangements over evolutionary time makes this process error-prone and decreases the efficiency of software aiming to recover all nested full-length TEs.

Results

We created software that uses a greedy recursive algorithm to mine increasingly fragmented copies of full-length LTR retrotransposons in assembled genomes and other sequence data. The software called TE-greedy-nester considers not only sequence similarity but also the structure of elements. This new tool was tested on a set of natural and synthetic sequences and its accuracy was compared to similar software. We found TE-greedy-nester to be superior in a number of parameters, namely computation time and full-length TE recovery in highly nested regions.

Availability and implementation

http://gitlab.fi.muni.cz/lexa/nested.

Supplementary information

Supplementary data are available at Bioinformatics online.

Co-option of the lineage-specific LAVA retrotransposon in the gibbon genome

Co-option of transposable elements (TEs) to become part of existing or new enhancers is an important mechanism for evolution of gene regulation. However, contributions of lineage-specific TE insertions to recent regulatory adaptations remain poorly understood. Gibbons present a suitable model to study these contributions as they have evolved a lineage-specific TE called LAVA (LINE-AluSz-VNTR-Alu _LIKE), which is still active in the gibbon genome. The LAVA retrotransposon is thought to have played a role in the emergence of the highly rearranged structure of the gibbon genome by disrupting transcription of cell cycle genes. In this study, we investigated whether LAVA may have also contributed to the evolution of gene regulation by adopting enhancer function. We characterized fixed and polymorphic LAVA insertions across multiple gibbons and found 96 LAVA elements overlapping enhancer chromatin states. Moreover, LAVA was enriched in multiple transcription factor binding motifs, was bound by an important transcription factor (PU.1), and was associated with higher levels of gene expression in cis We found gibbon-specific signatures of purifying/positive selection at 27 LAVA insertions. Two of these insertions were fixed in the gibbon lineage and overlapped with enhancer chromatin states, representing putative co-opted LAVA enhancers. These putative enhancers were located within genes encoding SETD2 and RAD9A, two proteins that facilitate accurate repair of DNA double-strand breaks and prevent chromosomal rearrangement mutations. Co-option of LAVA in these genes may have influenced regulation of processes that preserve genome integrity. Our findings highlight the importance of considering lineage-specific TEs in studying evolution of gene regulatory elements.

Divergence of 3' ends as a driver of short interspersed nuclear element (SINE) evolution in the Salicaceae

Short interspersed nuclear elements (SINEs) are small, non-autonomous and heterogeneous retrotransposons that are widespread in plants. To explore the amplification dynamics and evolutionary history of SINE populations in representative deciduous tree species, we analyzed the genomes of the six following Salicaceae species: Populus deltoides, Populus euphratica, Populus tremula, Populus tremuloides, Populus trichocarpa, and Salix purpurea. We identified 11 Salicaceae SINE families (SaliS-I to SaliS-XI), comprising 27 077 full-length copies. Most of these families harbor segmental similarities, providing evidence for SINE emergence by reshuffling or heterodimerization. We observed two SINE groups, differing in phylogenetic distribution pattern, similarity and 3' end structure. These groups probably emerged during the 'salicoid duplication' (~65 million years ago) in the Salix-Populus progenitor and during the separation of the genus Salix (45-65 million years ago), respectively. In contrast to conserved 5' start motifs across species and SINE families, the 3' ends are highly variable in sequence and length. This extraordinary 3'-end variability results from mutations in the poly(A) tail, which were fixed by subsequent amplificational bursts. We show that the dissemination of newly evolved 3' ends is accomplished by a displacement of older motifs, leading to various 3'-end subpopulations within the SaliS families.

DeepTE: a computational method for de novo classification of transposons with convolutional neural network

Motivation

Transposable elements (TEs) classification is an essential step to decode their roles in genome evolution. With a large number of genomes from non-model species becoming available, accurate and efficient TE classification has emerged as a new challenge in genomic sequence analysis.

Results

We developed a novel tool, DeepTE, which classifies unknown TEs using convolutional neural networks (CNNs). DeepTE transferred sequences into input vectors based on k-mer counts. A tree structured classification process was used where eight models were trained to classify TEs into super families and orders. DeepTE also detected domains inside TEs to correct false classification. An additional model was trained to distinguish between non-TEs and TEs in plants. Given unclassified TEs of different species, DeepTE can classify TEs into seven orders, which include 15, 24 and 16 super families in plants, metazoans and fungi, respectively. In several benchmarking tests, DeepTE outperformed other existing tools for TE classification. In conclusion, DeepTE successfully leverages CNN for TE classification, and can be used to precisely classify TEs in newly sequenced eukaryotic genomes.

Availability and implementation

DeepTE is accessible at https://github.com/LiLabAtVT/DeepTE.

Supplementary information

Supplementary data are available at Bioinformatics online.

Nested plant LTR retrotransposons target specific regions of other elements, while all LTR retrotransposons often target palindromes and nucleosome-occupied regions: in silico study

Background

Nesting is common in LTR retrotransposons, especially in large genomes containing a high number of elements.

Results

We analyzed 12 plant genomes and obtained 1491 pairs of nested and original (pre-existing) LTR retrotransposons. We systematically analyzed mutual nesting of individual LTR retrotransposons and found that certain families, more often belonging to the Ty3/gypsy than Ty1/copia superfamilies, showed a higher nesting frequency as well as a higher preference for older copies of the same family ("autoinsertions"). Nested LTR retrotransposons were preferentially located in the 3'UTR of other LTR retrotransposons, while coding and regulatory regions (LTRs) are not commonly targeted. Insertions displayed a weak preference for palindromes and were associated with a strong positional pattern of higher predicted nucleosome occupancy. Deviation from randomness in target site choice was also found in 13,983 non-nested plant LTR retrotransposons.

Conclusions

We reveal that nesting of LTR retrotransposons is not random. Integration is correlated with sequence composition, secondary structure and the chromatin environment. Insertion into retrotransposon positions with a low negative impact on family fitness supports the concept of the genome being viewed as an ecosystem of various elements.

Expression of A New Endogenous Retrovirus-Associated Transcript in Hodgkin Lymphoma Cells

During characterization of a cDNA library from the Hodgkin lymphoma (HL) cell line L-1236, we discovered a new transcript derived from chromosome 1 at the long intergenic non-protein coding RNA 1768 (LINC01768)/colony stimulating factor 1 (CSF1) region. The first exon of this transcript from Hodgkin lymphoma cells (THOLE) starts in the predicted exon 4 of LINC01768 and is part of an endogenous retrovirus (ERV) from the HUERS-P1/LTR8 family. High expression of THOLE was only detectable in HL cell line L-1236. The expression of THOLE in L-1236 cell is another example for ERV/LTR-associated gene expression in HL cells. At the genome level, the HUERS-P1/LTR8 region including THOLE is only present in Hominoidea. The influence of ERV/LTRs on gene expression might explain the characteristic phenotype of human HL.

Environmentally-induced epigenetic conversion of a piRNA cluster

Transposable element (TE) activity is repressed in animal gonads by PIWI-interacting RNAs (piRNAs) produced by piRNA clusters. Current models in flies propose that germinal piRNA clusters are functionally defined by the maternal inheritance of piRNAs produced during the previous generation. Taking advantage of an inactive, but ready to go, cluster of P-element derived transgene insertions in Drosophila melanogaster, we show here that raising flies at high temperature (29°C) instead of 25°C triggers the stable conversion of this locus from inactive into actively producing functional piRNAs. The increase of antisense transcripts from the cluster at 29°C combined with the requirement of transcription of euchromatic homologous sequences, suggests a role of double stranded RNA in the production of de novo piRNAs. This report describes the first case of the establishment of an active piRNA cluster by environmental changes in the absence of maternal inheritance of homologous piRNAs.

Editorial note

This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).

FishTEDB: a collective database of transposable elements identified in the complete genomes of fish

Transposable elements (TEs) are important for host gene regulation and genome evolution. Consensus sequences of TEs can assist investigators in accelerating studies on TE origins, amplification, functions and evolution, as well as comparative analyses and prediction of TEs in different species. In evolution, physiology, ecology and heredity research, fish are important models. However, to date, no comprehensive resource for TE consensus sequences exists for fish. Here, we collected genome-wide data and developed a novel database, FishTEDB, including 27 bony fishes, 1 cartilaginous fish, 1 lamprey and 1 lancelet. De novo, structure-based and homology-based approaches were combined to detect TEs. The database is open-source and user-friendly, and users can browse, search and download all data. FishTEDB also provides GetORF, BLAST and HMMER tools to analyze sequences.

Database URL: http://www.fishtedb.org/

Ty3/Gypsy retrotransposons in the Pacific abalone Haliotis discus hannai: characterization and use for species identification in the genus Haliotis

Transposable elements are highly abundant elements that are present in all eukaryotic species. Here, we present a molecular description of abalone retrotransposon (Abret) elements. The genome of Haliotis discus hannai contains 130 Abret elements which were all Ty3/Gypsy retrotransposons. The Ty1/Copia elements were absent in the H. discus hannai genome. Most of the elements were not complete due to sequence truncation or coding region decay. However, three elements Abret-296, Abret-935, and Abret-3259 had most of the canonical features of LTR (long terminal repeat)-retrotransposons. There were several reading frame shifts in Abret-935 and Abret-3259 elements. Surprisingly, phylogenetic analysis indicated that all of the elements belonged to the Osvaldo lineage. The sequence divergence between LTRs revealed that the Abret elements were mostly active within 2 million years ago. Abret elements were used as molecular markers in SSAP analyses, which allowed clear distinction of different species in the genus Haliotis. The polymorphic markers were converted into SCAR markers for use in species identification by simple PCR in the Haliotis genus.

Russian Doll Genes and Complex Chromosome Rearrangements in Oxytricha trifallax

Ciliates have two different types of nuclei per cell, with one acting as a somatic, transcriptionally active nucleus (macronucleus; abbr. MAC) and another serving as a germline nucleus (micronucleus; abbr. MIC). Furthermore, Oxytricha trifallax undergoes extensive genome rearrangements during sexual conjugation and post-zygotic development of daughter cells. These rearrangements are necessary because the precursor MIC loci are often both fragmented and scrambled, with respect to the corresponding MAC loci. Such genome architectures are remarkably tolerant of encrypted MIC loci, because RNA-guided processes during MAC development reorganize the gene fragments in the correct order to resemble the parental MAC sequence. Here, we describe the germline organization of several nested and highly scrambled genes in Oxytricha trifallax These include cases with multiple layers of nesting, plus highly interleaved or tangled precursor loci that appear to deviate from previously described patterns. We present mathematical methods to measure the degree of nesting between precursor MIC loci, and revisit a method for a mathematical description of scrambling. After applying these methods to the chromosome rearrangement maps of O. trifallax we describe cases of nested arrangements with up to five layers of embedded genes, as well as the most scrambled loci in O. trifallax.

Genome-wide analysis of transposable elements in the coffee berry borer Hypothenemus hampei (Coleoptera: Curculionidae): description of novel families

The coffee berry borer (CBB) Hypothenemus hampei is the most limiting pest of coffee production worldwide. The CBB genome has been recently sequenced; however, information regarding the presence and characteristics of transposable elements (TEs) was not provided. Using systematic searching strategies based on both de novo and homology-based approaches, we present a library of TEs from the draft genome of CBB sequenced by the Colombian Coffee Growers Federation. The library consists of 880 sequences classified as 66% Class I (LTRs: 46%, non-LTRs: 20%) and 34% Class II (DNA transposons: 8%, Helitrons: 16% and MITEs: 10%) elements, including families of the three main LTR (Gypsy, Bel-Pao and Copia) and non-LTR (CR1, Daphne, I/Nimb, Jockey, Kiri, R1, R2 and R4) clades and DNA superfamilies (Tc1-mariner, hAT, Merlin, P, PIF-Harbinger, PiggyBac and Helitron). We propose the existence of novel families: Hypo, belonging to the LTR Gypsy superfamily; Hamp, belonging to non-LTRs; and rosa, belonging to Class II or DNA transposons. Although the rosa clade has been previously described, it was considered to be a basal subfamily of the mariner family. Based on our phylogenetic analysis, including Tc1, mariner, pogo, rosa and Lsra elements from other insects, we propose that rosa and Lsra elements are subfamilies of an independent family of Class II elements termed rosa. The annotations obtained indicate that a low percentage of the assembled CBB genome (approximately 8.2%) consists of TEs. Although these TEs display high diversity, most sequences are degenerate, with few full-length copies of LTR and DNA transposons and several complete and putatively active copies of non-LTR elements. MITEs constitute approximately 50% of the total TEs content, with a high proportion associated with DNA transposons in the Tc1-mariner superfamily.

Characterization of new transposable element sub-families from white clover (Trifolium repens) using PCR amplification

Transposable elements (TEs) dominate the landscapes of most plant and animal genomes. Once considered junk DNA and genetic parasites, these interspersed, repetitive DNA elements are now known to play major roles in both genetic and epigenetic processes that sponsor genome variation and regulate gene expression. Knowledge of TE consensus sequences from elements in species whose genomes have not been sequenced is limited, and the individual TEs that are encountered in clones or short-reads rarely represent potentially canonical, let alone, functional representatives. In this study, we queried the Repbase database with eight BAC clones from white clover (Trifolium repens), identified a large number of candidate TEs, and used polymerase chain reaction and Sanger sequencing to create consensus sequences for three new TE families. The results show that TE family consensus sequences can be obtained experimentally in species for which just a single, full-length member of a TE family has been sequenced.

Development of a D genome specific marker resource for diploid and hexaploid wheat

BACKGROUND

Mapping and map-based cloning of genes that control agriculturally and economically important traits remain great challenges for plants with complex highly repetitive genomes such as those within the grass tribe, Triticeae. Mapping limitations in the Triticeae are primarily due to low frequencies of polymorphic gene markers and poor genetic recombination in certain genetic regions. Although the abundance of repetitive sequence may pose common problems in genome analysis and sequence assembly of large and complex genomes, they provide repeat junction markers with random and unbiased distribution throughout chromosomes. Hence, development of a high-throughput mapping technology that combine both gene-based and repeat junction-based markers is needed to generate maps that have better coverage of the entire genome.

RESULTS

In this study, the available genomics resource of the diploid Aegilop tauschii, the D genome donor of bread wheat, were used to develop genome specific markers that can be applied for mapping in modern hexaploid wheat. A NimbleGen array containing both gene-based and repeat junction probe sequences derived from Ae. tauschii was developed and used to map the Chinese Spring nullisomic-tetrasomic lines and deletion bin lines of the D genome chromosomes. Based on these mapping data, we have now anchored 5,171 repeat junction probes and 10,892 gene probes, corresponding to 5,070 gene markers, to the delineated deletion bins of the D genome. The order of the gene-based markers within the deletion bins of the Chinese Spring can be inferred based on their positions on the Ae. tauschii genetic map. Analysis of the probe sequences against the Chinese Spring chromosome sequence assembly database facilitated mapping of the NimbleGen probes to the sequence contigs and allowed assignment or ordering of these sequence contigs within the deletion bins. The accumulated length of anchored sequence contigs is about 155 Mb, representing ~ 3.2 % of the D genome. A specific database was developed to allow user to search or BLAST against the probe sequence information and to directly download PCR primers for mapping specific genetic loci.

CONCLUSIONS

In bread wheat, aneuploid stocks have been extensively used to assign markers linked with genes/traits to chromosomes, chromosome arms, and their specific bins. Through this study, we added thousands of markers to the existing wheat chromosome bin map, representing a significant step forward in providing a resource to navigate the wheat genome. The database website ( http://probes.pw.usda.gov/ATRJM/ ) provides easy access and efficient utilization of the data. The resources developed herein can aid map-based cloning of traits of interest and the sequencing of the D genome of hexaploid wheat.

Highly expressed captured genes and cross-kingdom domains present in Helitrons create novel diversity in Pleurotus ostreatus and other fungi

BACKGROUND

Helitrons are class-II eukaryotic transposons that transpose via a rolling circle mechanism. Due to their ability to capture and mobilize gene fragments, they play an important role in the evolution of their host genomes. We have used a bioinformatics approach for the identification of helitrons in two Pleurotus ostreatus genomes using de novo detection and homology-based searching. We have analyzed the presence of helitron-captured genes as well as the expansion of helitron-specific helicases in fungi and performed a phylogenetic analysis of their conserved domains with other representative eukaryotic species.

RESULTS

Our results show the presence of two helitron families in P. ostreatus that disrupt gene colinearity and cause a lack of synteny between their genomes. Both putative autonomous and non-autonomous helitrons were transcriptionally active, and some of them carried highly expressed captured genes of unknown origin and function. In addition, both families contained eukaryotic, bacterial and viral domains within the helitron's boundaries. A phylogenetic reconstruction of RepHel helicases using the Helitron-like and PIF1-like helicase conserved domains revealed a polyphyletic origin for eukaryotic helitrons.

CONCLUSION

P. ostreatus helitrons display features similar to other eukaryotic helitrons and do not tend to capture host genes or gene fragments. The occurrence of genes probably captured from other hosts inside the helitrons boundaries pose the hypothesis that an ancient horizontal transfer mechanism could have taken place. The viral domains found in some of these genes and the polyphyletic origin of RepHel helicases in the eukaryotic kingdom suggests that virus could have played a role in a putative lateral transfer of helitrons within the eukaryotic kingdom. The high similarity of some helitrons, along with the transcriptional activity of its RepHel helicases indicates that these elements are still active in the genome of P. ostreatus.

Genome-wide development of transposable elements-based markers in foxtail millet and construction of an integrated database

Transposable elements (TEs) are major components of plant genome and are reported to play significant roles in functional genome diversity and phenotypic variations. Several TEs are highly polymorphic for insert location in the genome and this facilitates development of TE-based markers for various genotyping purposes. Considering this, a genome-wide analysis was performed in the model plant foxtail millet. A total of 30,706 TEs were identified and classified as DNA transposons (24,386), full-length Copia type (1,038), partial or solo Copia type (10,118), full-length Gypsy type (1,570), partial or solo Gypsy type (23,293) and Long- and Short-Interspersed Nuclear Elements (3,659 and 53, respectively). Further, 20,278 TE-based markers were developed, namely Retrotransposon-Based Insertion Polymorphisms (4,801, ∼24%), Inter-Retrotransposon Amplified Polymorphisms (3,239, ∼16%), Repeat Junction Markers (4,451, ∼22%), Repeat Junction-Junction Markers (329, ∼2%), Insertion-Site-Based Polymorphisms (7,401, ∼36%) and Retrotransposon-Microsatellite Amplified Polymorphisms (57, 0.2%). A total of 134 Repeat Junction Markers were screened in 96 accessions of Setaria italica and 3 wild Setaria accessions of which 30 showed polymorphism. Moreover, an open access database for these developed resources was constructed (Foxtail millet Transposable Elements-based Marker Database; http://59.163.192.83/ltrdb/index.html). Taken together, this study would serve as a valuable resource for large-scale genotyping applications in foxtail millet and related grass species.

Identification of a recently active Prunus-specific non-autonomous Mutator element with considerable genome shaping force

Miniature inverted-repeat transposable elements (MITEs) are known to contribute to the evolution of plants, but only limited information is available for MITEs in the Prunus genome. We identified a MITE that has been named Falling Stones, FaSt. All structural features (349-bp size, 82-bp terminal inverted repeats and 9-bp target site duplications) are consistent with this MITE being a putative member of the Mutator transposase superfamily. FaSt showed a preferential accumulation in the short AT-rich segments of the euchromatin region of the peach genome. DNA sequencing and pollination experiments have been performed to confirm that the nested insertion of FaSt into the S-haplotype-specific F-box gene of apricot resulted in the breakdown of self-incompatibility (SI). A bioinformatics-based survey of the known Rosaceae and other genomes and a newly designed polymerase chain reaction (PCR) assay verified the Prunoideae-specific occurrence of FaSt elements. Phylogenetic analysis suggested a recent activity of FaSt in the Prunus genome. The occurrence of a nested insertion in the apricot genome further supports the recent activity of FaSt in response to abiotic stress conditions. This study reports on a presumably active non-autonomous Mutator element in Prunus that exhibits a major indirect genome shaping force through inducing loss-of-function mutation in the SI locus.

Genomic landscape of human, bat, and ex vivo DNA transposon integrations

The integration and fixation preferences of DNA transposons, one of the major classes of eukaryotic transposable elements, have never been evaluated comprehensively on a genome-wide scale. Here, we present a detailed study of the distribution of DNA transposons in the human and bat genomes. We studied three groups of DNA transposons that integrated at different evolutionary times: 1) ancient (>40 My) and currently inactive human elements, 2) younger (<40 My) bat elements, and 3) ex vivo integrations of piggyBat and Sleeping Beauty elements in HeLa cells. Although the distribution of ex vivo elements reflected integration preferences, the distribution of human and (to a lesser extent) bat elements was also affected by selection. We used regression techniques (linear, negative binomial, and logistic regression models with multiple predictors) applied to 20-kb and 1-Mb windows to investigate how the genomic landscape in the vicinity of DNA transposons contributes to their integration and fixation. Our models indicate that genomic landscape explains 16-79% of variability in DNA transposon genome-wide distribution. Importantly, we not only confirmed previously identified predictors (e.g., DNA conformation and recombination hotspots) but also identified several novel predictors (e.g., signatures of double-strand breaks and telomere hexamer). Ex vivo integrations showed a bias toward actively transcribed regions. Older DNA transposons were located in genomic regions scarce in most conserved elements-likely reflecting purifying selection. Our study highlights how DNA transposons are integral to the evolution of bat and human genomes, and has implications for the development of DNA transposon assays for gene therapy and mutagenesis applications.

Sequence-Based Analysis of Structural Organization and Composition of the Cultivated Sunflower (Helianthus annuus L.) Genome

Sunflower is an important oilseed crop, as well as a model system for evolutionary studies, but its 3.6 gigabase genome has proven difficult to assemble, in part because of the high repeat content of its genome. Here we report on the sequencing, assembly, and analyses of 96 randomly chosen BACs from sunflower to provide additional information on the repeat content of the sunflower genome, assess how repetitive elements in the sunflower genome are organized relative to genes, and compare the genomic distribution of these repeats to that found in other food crops and model species. We also examine the expression of transposable element-related transcripts in EST databases for sunflower to determine the representation of repeats in the transcriptome and to measure their transcriptional activity. Our data confirm previous reports in suggesting that the sunflower genome is >78% repetitive. Sunflower repeats share very little similarity to other plant repeats such as those of Arabidopsis, rice, maize and wheat; overall 28% of repeats are “novel” to sunflower. The repetitive sequences appear to be randomly distributed within the sequenced BACs. Assuming the 96 BACs are representative of the genome as a whole, then approximately 5.2% of the sunflower genome comprises non TE-related genic sequence, with an average gene density of 18kbp/gene. Expression levels of these transposable elements indicate tissue specificity and differential expression in vegetative and reproductive tissues, suggesting that expressed TEs might contribute to sunflower development. The assembled BACs will also be useful for assessing the quality of several different draft assemblies of the sunflower genome and for annotating the reference sequence.

A MATLAB-based tool for accurate detection of perfect overlapping and nested inverted repeats in DNA sequences

SUMMARY

Palindromic sequences, or inverted repeats (IRs), in DNA sequences involve important biological processes such as DNA-protein binding, DNA replication and DNA transposition. Development of bioinformatics tools that are capable of accurately detecting perfect IRs can enable genome-wide studies of IR patterns in both prokaryotes and eukaryotes. Different from conventional string-comparison approaches, we propose a novel algorithm that uses a cumulative score system based on a prime number representation of nucleotide bases. We then implemented this algorithm as a MATLAB-based program for perfect IR detection. In comparison with other existing tools, our program demonstrates a high accuracy in detecting nested and overlapping IRs.

AVAILABILITY AND IMPLEMENTATION

The source code is freely available on (http://bioinfolab.miamioh.edu/bioinfolab/palindrome.php)

CONTACT

liangc@miamioh.edu or karroje@miamioh.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Internal deletions of transposable elements: the case of Lemi elements

Mobile elements using a "cut and paste" mechanism of transposition (Class II) are frequently prone to internal deletions and the question of the origin of these copies remains elusive. In this study, we looked for copies belonging to the Lemi Family (Tc1-mariner-IS630 SuperFamily) in the plant genomes, and copies within internal deletions were analyzed in detail. Lemi elements are found exclusively in Eudicots, and more than half of the copies have been deleted. All deletions occur between microhomologies (direct repeats from 2 to 13 bp). Copies less than 500 bp long, similar to MITEs, are frequent. These copies seem to result from large deletions occurring between microhomologies present within a region of 300 bp at both extremities of the element. These regions are particularly A/T rich, compared to the internal part of the element, which increases the probability of observing short direct repeats. Most of the molecular mechanisms responsible for double strand break repair are able to induce deletions between microhomologies during the repair process. This could be a quick way to reduce the population of active copies within a genome and, more generally, to reduce the overall activity of the element after it has entered a naive genome.