Short Transposable Elements Research Articles

A 21L/BA.2-21K/BA.1 "MixOmicron" SARS-CoV-2 hybrid undetected by qPCR that screen for variant in routine diagnosis.

Among the multiple SARS-CoV-2 variants identified since summer 2020, several have co-circulated, creating opportunities for coinfections and potentially genetic recombinations that are common in coronaviruses. Viral recombinants are indeed beginning to be reported more frequently. Here, we describe a new SARS-CoV-2 recombinant genome that is mostly that of a Omicron 21L/BA.2 variant but with a 3′ tip originating from a Omicron 21K/BA.1 variant. Two such genomes were obtained in our institute from adults sampled in February 2022 in university hospitals of Marseille, southern France, by next-generation sequencing carried out with the Illumina or Nanopore technologies. The recombination site was located between nucleotides 26,858-27,382. In the two genomic assemblies, mean sequencing depth at mutation-harboring positions was 271 and 1362 reads and mean prevalence of the majoritary nucleotide was 99.3 ± 2.2% and 98.8 ± 1.6%, respectively. Phylogeny generated trees with slightly different topologies according to whether genomes were depleted or not of the 3′ tip. This 3′ terminal end brought in the Omicron 21L/BA.2 genome a short transposable element of 41 nucleotides named S2m that is present in most SARS-CoV-2 except a few variants among which the Omicron 21L/BA.2 variant and may be involved in virulence. Importantly, this recombinant is not detected by currently used qPCR that screen for variants in routine diagnosis. The present observation emphasizes the need to survey closely the genetic pathways of SARS-CoV-2 variability by whole genome sequencing, and it could contribute to gain a better understanding of factors that lead to observed differences between epidemic potentials of the different variants.

Disentangling the determinants of transposable elements dynamics in vertebrate genomes using empirical evidences and simulations.

The interactions between transposable elements (TEs) and their hosts constitute one of the most profound co-evolutionary processes found in nature. The population dynamics of TEs depends on factors specific to each TE families, such as the rate of transposition and insertional preference, the demographic history of the host and the genomic landscape. How these factors interact has yet to be investigated holistically. Here we are addressing this question in the green anole (Anolis carolinensis) whose genome contains an extraordinary diversity of TEs (including non-LTR retrotransposons, SINEs, LTR-retrotransposons and DNA transposons). We observed a positive correlation between recombination rate and frequency of TEs and densities for LINEs, SINEs and DNA transposons. For these elements, there was a clear impact of demography on TE frequency and abundance, with a loss of polymorphic elements and skewed frequency spectra in recently expanded populations. On the other hand, some LTR-retrotransposons displayed patterns consistent with a very recent phase of intense amplification. To determine how demography, genomic features and intrinsic properties of TEs interact we ran simulations using SLiM3. We determined that i) short TE insertions are not strongly counter-selected, but long ones are, ii) neutral demographic processes, linked selection and preferential insertion may explain positive correlations between average TE frequency and recombination, iii) TE insertions are unlikely to have been massively recruited in recent adaptation. We demonstrate that deterministic and stochastic processes have different effects on categories of TEs and that a combination of empirical analyses and simulations can disentangle these mechanisms.

Transcriptional regulation of genes bearing intronic heterochromatin in the rice genome.

Intronic regions of eukaryotic genomes accumulate many Transposable Elements (TEs). Intronic TEs often trigger the formation of transcriptionally repressive heterochromatin, even within transcription-permissive chromatin environments. Although TE-bearing introns are widely observed in eukaryotic genomes, their epigenetic states, impacts on gene regulation and function, and their contributions to genetic diversity and evolution, remain poorly understood. In this study, we investigated the genome-wide distribution of intronic TEs and their epigenetic states in the Oryza sativa genome, where TEs comprise 35% of the genome. We found that over 10% of rice genes contain intronic heterochromatin, most of which are associated with TEs and repetitive sequences. These heterochromatic introns are longer and highly enriched in promoter-proximal positions. On the other hand, introns also accumulate hypomethylated short TEs. Genes with heterochromatic introns are implicated in various biological functions. Transcription of genes bearing intronic heterochromatin is regulated by an epigenetic mechanism involving the conserved factor OsIBM2, mutation of which results in severe developmental and reproductive defects. Furthermore, we found that heterochromatic introns evolve rapidly compared to non-heterochromatic introns. Our study demonstrates that heterochromatin is a common epigenetic feature associated with actively transcribed genes in the rice genome.

MUSTv2: An Improved De Novo Detection Program for Recently Active Miniature Inverted Repeat Transposable Elements (MITEs).

Background Miniature inverted repeat transposable element (MITE) is a short transposable element, carrying no protein-coding regions. However, its high proliferation rate and sequence-specific insertion preference renders it as a good genetic tool for both natural evolution and experimental insertion mutagenesis. Recently active MITE copies are those with clear signals of Terminal Inverted Repeats (TIRs) and Direct Repeats (DRs), and are recently translocated into their current sites. Their proliferation ability renders them good candidates for the investigation of genomic evolution.Results This study optimizes the C++ code and running pipeline of the MITE Uncovering SysTem (MUST) by assuming no prior knowledge of MITEs required from the users, and the current version, MUSTv2, shows significantly increased detection accuracy for recently active MITEs, compared with similar programs. The running speed is also significantly increased compared with MUSTv1. We prepared a benchmark dataset, the simulated genome with 150 MITE copies for researchers who may be of interest.Conclusions MUSTv2 represents an accurate detection program of recently active MITE copies, which is complementary to the existing template-based MITE mapping programs. We believe that the release of MUSTv2 will greatly facilitate the genome annotation and structural analysis of the bioOMIC big data researchers.

DNA demethylases target promoter transposable elements to positively regulate stress responsive genes in Arabidopsis.

BackgroundDNA demethylases regulate DNA methylation levels in eukaryotes. Arabidopsis encodes four DNA demethylases, DEMETER (DME), REPRESSOR OF SILENCING 1 (ROS1), DEMETER-LIKE 2 (DML2), and DML3. While DME is involved in maternal specific gene expression during seed development, the biological function of the remaining DNA demethylases remains unclear.ResultsWe show that ROS1, DML2, and DML3 play a role in fungal disease resistance in Arabidopsis. A triple DNA demethylase mutant, rdd (ros1 dml2 dml3), shows increased susceptibility to the fungal pathogen Fusarium oxysporum. We identify 348 genes differentially expressed in rdd relative to wild type, and a significant proportion of these genes are downregulated in rdd and have functions in stress response, suggesting that DNA demethylases maintain or positively regulate the expression of stress response genes required for F. oxysporum resistance. The rdd-downregulated stress response genes are enriched for short transposable element sequences in their promoters. Many of these transposable elements and their surrounding sequences show localized DNA methylation changes in rdd, and a general reduction in CHH methylation, suggesting that RNA-directed DNA methylation (RdDM), responsible for CHH methylation, may participate in DNA demethylase-mediated regulation of stress response genes. Many of the rdd-downregulated stress response genes are downregulated in the RdDM mutants nrpd1 and nrpe1, and the RdDM mutants nrpe1 and ago4 show enhanced susceptibility to F. oxysporum infection.ConclusionsOur results suggest that a primary function of DNA demethylases in plants is to regulate the expression of stress response genes by targeting promoter transposable element sequences.Electronic supplementary materialThe online version of this article (doi:10.1186/s13059-014-0458-3) contains supplementary material, which is available to authorized users.

DNA demethylases target promoter transposable elements to positively regulate stress responsive genes in Arabidopsis

BackgroundDNA demethylases regulate DNA methylation levels in eukaryotes. Arabidopsis encodes four DNA demethylases, DEMETER (DME), REPRESSOR OF SILENCING 1 (ROS1), DEMETER-LIKE 2 (DML2), and DML3. While DME is involved in maternal specific gene expression during seed development, the biological function of the remaining DNA demethylases remains unclear.ResultsWe show that ROS1, DML2, and DML3 play a role in fungal disease resistance in Arabidopsis. A triple DNA demethylase mutant, rdd (ros1 dml2 dml3), shows increased susceptibility to the fungal pathogen Fusarium oxysporum. We identify 348 genes differentially expressed in rdd relative to wild type, and a significant proportion of these genes are downregulated in rdd and have functions in stress response, suggesting that DNA demethylases maintain or positively regulate the expression of stress response genes required for F. oxysporum resistance. The rdd-downregulated stress response genes are enriched for short transposable element sequences in their promoters. Many of these transposable elements and their surrounding sequences show localized DNA methylation changes in rdd, and a general reduction in CHH methylation, suggesting that RNA-directed DNA methylation (RdDM), responsible for CHH methylation, may participate in DNA demethylase-mediated regulation of stress response genes. Many of the rdd-downregulated stress response genes are downregulated in the RdDM mutants nrpd1 and nrpe1, and the RdDM mutants nrpe1 and ago4 show enhanced susceptibility to F. oxysporum infection.ConclusionsOur results suggest that a primary function of DNA demethylases in plants is to regulate the expression of stress response genes by targeting promoter transposable element sequences.

TSIDER1, a short and non-autonomous Salivarian trypanosome-specific retroposon related to the ingi6 subclade

Retroposons of the ingi clade are the most abundant transposable elements identified in the trypanosomatid genomes. Some are long autonomous elements (ingi, L1Tc) while others, such as RIME and NARTc, are short non-coding elements that parasitize the retrotransposition machinery of the active autonomous ones for their own mobilization. Here, we identified a new family of short non-autonomous retroposons of the ingi clade, called TSIDER1, which are present in the genome of Salivarian (African) trypanosomes, Trypanosoma brucei, T. congolense and T. vivax, but absent in the T. cruzi and Leishmania spp. genomes and, as such, TSIDER1 is the only retroposon subfamily conserved at the nucleotide level between African trypanosome species. We identified three TvSIDER1 families within the genome of T. vivax and the high level of sequence conservation within the TvSIDER1a and TvSIDER1b groups suggests that they are still active. We propose that TvSIDER1a/b elements are using the Tvingi retrotransposition machinery, as they are preceded by the same conserved pattern characteristic of the ingi6 subclade, which corresponds to the retroposon-encoded endonuclease binding site. In contrast, TcoSIDER1, TbSIDER1 and TvSIDER1c are too divergent to be considered as active retroposons. The relatively low number of SIDER elements identified in the T. congolense (70 copies), T. vivax (32 copies) and T. brucei (22 copies) genomes confirms that trypanosomes have not expanded short transposable elements, which is in contrast to Leishmania spp. (∼2000 copies), where SIDER play a role in the regulation of gene expression.

Short- and Long-term Evolutionary Dynamics of Bacterial Insertion Sequences: Insights from Wolbachia Endosymbionts

Transposable elements (TE) are one of the major driving forces of genome evolution, raising the question of the long-term dynamics underlying their evolutionary success. Long-term TE evolution can readily be reconstructed in eukaryotes, thanks to many degraded copies constituting genomic fossil records of past TE proliferations. By contrast, bacterial genomes usually experience high sequence turnover and short TE retention times, thereby obscuring ancient TE evolutionary patterns. We found that Wolbachia bacterial genomes contain 52–171 insertion sequence (IS) TEs. IS account for 11% of Wolbachia wRi, which is one of the highest IS genomic coverage reported in prokaryotes to date. We show that many IS groups are currently expanding in various Wolbachia genomes and that IS horizontal transfers are frequent among strains, which can explain the apparent synchronicity of these IS proliferations. Remarkably, >70% of Wolbachia IS are nonfunctional. They constitute an unusual bacterial IS genomic fossil record providing direct empirical evidence for a long-term IS evolutionary dynamics following successive periods of intense transpositional activity. Our results show that comprehensive IS annotations have the potential to provide new insights into prokaryote TE evolution and, more generally, prokaryote genome evolution. Indeed, the identification of an important IS genomic fossil record in Wolbachia demonstrates that IS elements are not always of recent origin, contrary to the conventional view of TE evolution in prokaryote genomes. Our results also raise the question whether the abundance of IS fossils is specific to Wolbachia or it may be a general, albeit overlooked, feature of prokaryote genomes.

Identification of miniature inverted-repeat transposable elements (MITEs) and biogenesis of their siRNAs in the Solanaceae: new functional implications for MITEs.

Small RNAs regulate the genome by guiding transcriptional and post-transcriptional silencing machinery to specific target sequences, including genes and transposable elements (TEs). Although miniature inverted-repeat transposable elements (MITEs) are closely associated with euchromatic genes, the broader functional impact of these short TE insertions in genes is largely unknown. We identified 22 families of MITEs in the Solanaceae (MiS1-MiS22) and found abundant MiS insertions in Solanaceae genomic DNA and expressed sequence tags (EST). Several Solanaceae MITEs generate genome changes that potentially affect gene function and regulation, most notably, a MiS insertion that provides a functionally indispensable alternative exon in the tobacco mosaic virus N resistance gene. We show that MITEs generate small RNAs that are primarily 24 nt in length, as detected by Northern blot hybridization and by sequencing small RNAs of Solanum demissum, Nicotiana glutinosa, and Nicotiana benthamiana. Additionally, we show that stable RNAi lines silencing DICER-LIKE3 (DCL3) in tobacco and RNA-dependent RNA polymerase 2 (RDR2) in potato cause a reduction in 24-nt MITE siRNAs, suggesting that, as in Arabidopsis, TE-derived siRNA biogenesis is DCL3 and RDR2 dependent. We provide evidence that DICER-LIKE4 (DCL4) may also play a role in MITE siRNA generation in the Solanaceae.

Evolution and phylogenetic utility of the PHOT gene duplicates in the Verbena complex (Verbenaceae): dramatic intron size variation and footprint of ancestral recombination

A well-resolved species level phylogeny is critically important in studying organismal evolution (e.g., hybridization, polyploidization, adaptive speciation). Lack of appropriate molecular markers that give sufficient resolution to gene trees is one of the major impediments to inferring species level phylogenies. In addition, sampling multiple independent loci is essential to overcome the lineage sorting problem. The availability of nuclear loci has often been a limiting factor in plant species-level phylogenetic studies. Here the two PHOT loci were developed as new sources of nuclear gene trees. The PHOT1 and PHOT2 gene trees of the Verbena complex (Verbenaceae) are well resolved and have good clade support. These gene trees are consistent with each other and previously generated chloroplast and nuclear waxy gene trees in most of the phylogenetic backbone as well as some terminal relationships, but are incongruent in some other relationships. Locus-specific primers were optimized for amplifying and sequencing these two loci in all Lamiales. Comparing intron size in the context of the gene trees shows dramatic variation within the Verbena complex, particularly at the PHOT1 locus. These variations are largely caused by invasions of short transposable elements and frequent long deletions and insertions of unknown causes. In addition, inspection of DNA sequences and phylogenetic analyses unmask a clear footprint of ancestral recombination in one species.

Peas-Mea1-Ppp2r5d overlapping gene complex: a transposon mediated-gene formation in mammals.

Human and mouse MEA1/Mea1 is flanked by two overlapping genes, a novel PEAS/Peas in a head-to-head orientation and PPP2R5D/Ppp2r5d in a tail-to-tail orientation making a Peas-Mea1-Ppp2r5d overlapping gene complex (PMP-complex). Genomic zoo blot analyses and database searching revealed that Mea1 exists only in mammals, while Peas and Ppp2r5d are conserved in eukaryotes. Mea1 and Peas are transcribed from a testis-expressed bidirectional promoter. Mea1-Ppp2r5d overlapping segment (MPOS) contains polyadenylation signals for both genes and shows marked conservation throughout mammals. Furthermore, the MPOS occupies 3'-region of transcripts of both genes is expected to form a clover-like intramolecular structure. Mouse genomic library Screening and database searches identified two MPOS-derived sequences in Odf2 gene and RP23-86H7 cosmid clone, respectively, in which MPOS might be a core segment for the retropositions. Thus, a key role of MPOS, a short transposable element containing polyadenylation signals on both strands, in the formation of the Mea1 during mammalian evolution is suggested.

One member of the tRNA(Glu) gene family in yeast codes for a minor GAGtRNA(Glu) species and is associated with several short transposable elements.

During characterization of the whole tRNA-(Glu) family from the yeast, Saccharomyces cerevisiae, we isolated one cosmid clone bearing a tRNA(Glu) gene copy that is deviant from the major tRNA(Glu3) gene members in only five positions. This divergent tRNA-(Glu) is a minor species and is represented by a single gene copy. One of the nucleotide exchanges concerns the anticodon which is modified from T-T-C in the tRNA(Glu3) gene to C-T-C which implies that this tRNA serves the codon triplet G-A-G. Two other minor yeast tRNA species have been reported which appear to be particularly designed for the translation of those codons that have a G in its third (Wobble) position. The low abundance of such minor tRNA species correlates positively to the low occurrence of most of the N-N-G codons in yeast. Furthermore, the GAGtRNA-(Glu) locus represents another case of the general phenomenon in which the majority of the tRNA genes in yeast are associated with one or several transposable elements forming complex patterns. In this particular case, divergent segments of delta and tau are present in the 5' flanking region of the tRNA gene and arranged in a novel configuration. The sequence data lend support to the view that tau is not an evolutionary young element as was earlier anticipated.