Generalized nuclear localization of retroelement transcripts

BACKGROUND

LINE-1s, Alus and SVAs are the only retrotransposition competent elements in humans. Their mobilization followed by insertional mutagenesis is often linked to disease. Apart from these rare integration events, accumulation of retrotransposition intermediates in the cytoplasm is potentially pathogenic due to induction of inflammatory response pathways. Although the retrotransposition of LINE-1 and Alu retroelements has been studied in considerable detail, there are mixed observations about the localization of their RNAs.

RESULTS

We undertook a comprehensive and unbiased approach to analyze retroelement RNA localization using common cell lines and publicly available datasets containing RNA-sequencing data from subcellular fractions. Using our customized analytic pipeline, we compared localization patterns of RNAs transcribed from retroelements and single-copy protein coding genes. Our results demonstrate a generalized characteristic pattern of retroelement RNA nuclear localization that is conserved across retroelement classes as well as evolutionarily young and ancient elements. Preferential nuclear enrichment of retroelement transcripts was consistently observed in cell lines, in vivo and across species. Moreover, retroelement RNA localization patterns were dynamic and subject to change during development, as seen in zebrafish embryos.

CONCLUSION

The pronounced nuclear localization of transcripts arising from ancient as well as de novo transcribed retroelements suggests that these transcripts are retained in the nucleus as opposed to being re-imported to the nucleus or degraded in the cytoplasm. This raises the possibility that there is adaptive value associated with this localization pattern to the host, the retroelements or possibly both.

Expression of Retroelements in Mammalian Gametes and Embryos

Retroelements are genetic mobile elements, expressed during male and female gamete differentiation. Retrotransposons are normally regulated by the methylation machinery, chromatin modifications, non-coding RNAs, and transcription factors, while retrotransposition control is of vital importance in cellular proliferation and differentiation process. Retrotransposition requires a transcription step, by a cellular RNA polymerase, followed by reverse transcription of an RNA intermediate to cDNA and its integration into a new genomic locus. Long interspersed elements (LINEs), human endogenous retroviruses (HERVs), short interspersed elements (SINEs) and SINE-VNTR-Alu elements (SVAs) constitute about half of the human genome, play a crucial role in genome organization, structure and function and interfere with several biological procedures. In this mini review, we discuss recent data regarding retroelement expression (LINE-1, HERVK-10, SVA and VL30) and retrotransposition events in mammalian oocytes and spermatozoa, as well as the importance of their impact on human and mouse preimplantation embryo development.

Haplotype Shuffling and Dimorphic Transposable Elements in the Human Extended Major Histocompatibility Complex Class II Region

The major histocompatibility complex (MHC) on chromosome 6p21 is one of the most single-nucleotide polymorphism (SNP)-dense regions of the human genome and a prime model for the study and understanding of conserved sequence polymorphisms and structural diversity of ancestral haplotypes/conserved extended haplotypes. This study aimed to follow up on a previous analysis of the MHC class I region by using the same set of 95 MHC haplotype sequences downloaded from a publicly available BioProject database at the National Center for Biotechnology Information to identify and characterize the polymorphic human leukocyte antigen (HLA)-class II genes, the MTCO3P1 pseudogene alleles, the indels of transposable elements as haplotypic lineage markers, and SNP-density crossover (XO) loci at haplotype junctions in DNA sequence alignments of different haplotypes across the extended class II region (∼1 Mb) from the telomeric PRRT1 gene in class III to the COL11A2 gene at the centromeric end of class II. We identified 42 haplotypic indels (20 Alu, 7 SVA, 13 LTR or MERs, and 2 indels composed of a mosaic of different transposable elements) linked to particular HLA-class II alleles. Comparative sequence analyses of 136 haplotype pairs revealed 98 unique XO sites between SNP-poor and SNP-rich genomic segments with considerable haplotype shuffling located in the proximity of putative recombination hotspots. The majority of XO sites occurred across various regions including in the vicinity of MTCO3P1 between HLA-DQB1 and HLA-DQB3, between HLA-DQB2 and HLA-DOB, between DOB and TAP2, and between HLA-DOA and HLA-DPA1, where most XOs were within a HERVK22 sequence. We also determined the genomic positions of the PRDM9-recombination suppression sequence motif ATCCATG/CATGGAT and the PRDM9 recombination activation partial binding motif CCTCCCCT/AGGGGAG in the class II region of the human reference genome (NC_ 000006) relative to published meiotic recombination positions. Both the recombination and anti-recombination PRDM9 binding motifs were widely distributed throughout the class II genomic regions with 50% or more found within repeat elements; the anti-recombination motifs were found mostly in L1 fragmented repeats. This study shows substantial haplotype shuffling between different polymorphic blocks and confirms the presence of numerous putative ancestral recombination sites across the class II region between various HLA class II genes.

Detection of Transposition Events from Next-Generation Sequencing Data

Detection of transposition events of a transposon from short reads of next-generation sequencing (NGS) is challenging because transposons are repetitive and difficult to be distinguished from already existing transposons in the genome. Many transposons generate target site duplication (TSD) as the result of chromosomal integration. Since TSDs flanking the 5'-end (head) and 3'-end (tail) of a transposon has the identical sequences which are absent from the reference copy, the short reads containing the head or tail sequences of the transposon following the same TSD sequence may reveal the evidence of transposition. Transposon Insertion Finder (TIF) focuses on the TSD with flanking sequence of transposon and detects transposition events from NGS data. TIF software is available at https://github.com/akiomiyao/tif .

SeqURE - a new copy-capture based method for sequencing of unknown Retroposition events

Background

Retroelements (REs) occupy a significant part of all eukaryotic genomes including humans. The majority of retroelements in the human genome are inactive and unable to retrotranspose. Dozens of active copies are repressed in most normal tissues by various cellular mechanisms. These copies can become active in normal germline and brain tissues or in cancer, leading to new retroposition events. The consequences of such events and their role in normal cell functioning and carcinogenesis are not yet fully understood. If new insertions occur in a small portion of cells they can be found only with the use of specific methods based on RE enrichment and high-throughput sequencing. The downside of the high sensitivity of such methods is the presence of various artifacts imitating real insertions, which in many cases cannot be validated due to lack of the initial template DNA. For this reason, adequate assessment of rare (< 1%) subclonal cancer specific RE insertions is complicated.

Results

Here we describe a new copy-capture technique which we implemented in a method called SeqURE for Sequencing Unknown of Retroposition Events that allows for efficient and reliable identification of new genomic RE insertions. The method is based on the capture of copies of target molecules (copy-capture), selective amplification and sequencing of genomic regions adjacent to active RE insertions from both sides. Importantly, the template genomic DNA remains intact and can be used for validation experiments. In addition, we applied a novel system for testing method sensitivity and precisely showed the ability of the developed method to reliably detect insertions present in 1 out of 100 cells and a substantial portion of insertions present in 1 out of 1000 cells. Using advantages of the method we showed the absence of somatic Alu insertions in colorectal cancer samples bearing tumor-specific L1HS insertions.

Conclusions

This study presents the first description and implementation of the copy-capture technique and provides the first methodological basis for the quantitative assessment of RE insertions present in a small portion of cells.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13100-020-00228-6.

A comparative study of long interspersed element-1 protein immunoreactivity in cutaneous malignancies

Background

Skin cancer is the most common cancer worldwide and commonly classified into malignant melanoma (MM) and Nonmelanoma skin cancers (NMSCs), which mainly include basal cell carcinoma (BCC) and squamous cell carcinoma (SCC). The extent to which Long Interspersed Element-1 (LINE-1, L1) ORF1p is expressed in cutaneous malignancies remains to be evaluated. This study aimed to assess LINE-1 ORF1p immunoreactivity in various skin cancer subtypes.

Method

The expression level of LINE-1 ORF1p was evaluated in 95 skin cancer specimens comprising 36 (37.9%) BCC, 28 (29.5%) SCC, and 31 (32.6%) melanoma using the tissue microarray (TMA) technique. Then the association between expression of LINE-1 encoded protein and clinicopathological parameters was analyzed.

Results

We showed that LINE-1 ORF1p expression level was substantially higher in BCC and SCC patients compared with melanoma samples (p < 0.001). BCC cases had a higher LINE-1 histochemical score (H-score) compared with SCC cases (p = 0.004). In SCC samples, a lower level of LINE-1 ORF1p expression was associated with age younger than the mean (p = 0.041). At the same time, no significant correlation was found between LINE-1 ORF1p expression and other clinicopathological parameters (all p > 0.05).

Conclusions

According to our observation, LINE-1 ORF1p immunoreactivity in various skin tumor subtypes extends previous studies of LINE-1 expression in different cancers. LINE-1ORF1p overexpression in NMSCs compared with MM can be considered with caution as a tumor-specific antigen for NMSCs.

Study of VIPER and TATE in kinetoplastids and the evolution of tyrosine recombinase retrotransposons

Background

Kinetoplastids are a flagellated group of protists, including some parasites, such as Trypanosoma and Leishmania species, that can cause diseases in humans and other animals. The genomes of these species enclose a fraction of retrotransposons including VIPER and TATE, two poorly studied transposable elements that encode a tyrosine recombinase (YR) and were previously classified as DIRS elements. This study investigated the distribution and evolution of VIPER and TATE in kinetoplastids to understand the relationships of these elements with other retrotransposons.

Results

We observed that VIPER and TATE have a discontinuous distribution among Trypanosomatidae, with several events of loss and degeneration occurring during a vertical transfer evolution. We were able to identify the terminal repeats of these elements for the first time, and we showed that these elements are potentially active in some species, including T. cruzi copies of VIPER. We found that VIPER and TATE are strictly related elements, which were named in this study as VIPER-like. The reverse transcriptase (RT) tree presented a low resolution, and the origin and relationships among YR groups remain uncertain. Conversely, for RH, VIPER-like grouped with Hepadnavirus, whereas for YR, VIPER-like sequences constituted two different clades that are closely allied to Crypton. Distinct topologies among RT, RH and YR trees suggest ancient rearrangements/exchanges in domains and a modular pattern of evolution with putative independent origins for each ORF.

Conclusions

Due to the presence of both elements in Bodo saltans, a nontrypanosomatid species, we suggested that VIPER and TATE have survived and remained active for more than 400 million years or were reactivated during the evolution of the host species. We did not find clear evidence of independent origins of VIPER-like from the other YR retroelements, supporting the maintenance of the DIRS group of retrotransposons. Nevertheless, according to phylogenetic findings and sequence structure obtained by this study and other works, we proposed separating DIRS elements into four subgroups: DIRS-like, PAT-like, Ngaro-like, and VIPER-like.

Electronic supplementary material

The online version of this article (10.1186/s13100-019-0175-2) contains supplementary material, which is available to authorized users.

Characterization of molecular attributes that influence LINE-1 restriction by all seven human APOBEC3 proteins

LINE-1 (L1) is a non-long terminal repeat (LTR) retrotransposon inserted throughout the human genome. APOBEC3 (A3) proteins are part of a network of host intrinsic defenses capable of restricting retroviruses and the replication of L1 retroelements. These enzymes inactivate retroviruses primarily through deamination of single-stranded viral DNA. In contrast, only A3A deaminates L1 DNA, while the other six A3 proteins restrict L1 to varying degrees through yet poorly defined mechanisms. Here we provide further insight into the molecular attributes of L1 restriction by A3 proteins. We specifically investigated the roles of A3 protein oligomerization, interactions with RNA and their binding to the various L1 proteins. Our results show that compromising the ability of A3 proteins to oligomerize or interact with a nucleic acid substrate diminished L1 restriction to varying degrees. However the efficiency of their binding to L1 proteins did not predict restriction or the potency of the restriction.

Genome-wide hypomethylation of LINE-1 and Alu retroelements in cell-free DNA of blood is an epigenetic biomarker of human aging

Aging associated DNA hypomethylation of LINE-1 and Alu retroelements may be a crucial determinant of loss of genomic integrity, deterioration and cancer. In peripheral blood LINE-1 hypomethylation has been reported to increase during aging, but other studies did not observe significant changes. We hypothesized that these apparently inconsistent reports might relate to differences between cellular and cell-free DNA. Using the technique of idiolocal normalization of real-time methylation-specific PCR (IDLN-MSP) for genetic imbalanced DNA specimens we obtained evidence that LINE-1 hypomethylation in cell-free DNA, but not cellular DNA from peripheral blood is an epigenetic biomarker for human aging. Furthermore, hypomethylation of cell-free DNA is more extensive in smokers, suggesting that it might be used as a surrogate marker for monitoring the improvement of smoking-induced adverse effects after cancelling smoking.

Role of the Unconventional Prefoldin Proteins URI and UXT in Transcription Regulation

The Unconventional prefoldin RPB5 interacting protein (URI), also known as RPB5-Mediating Protein (RMP) has been shown to play several regulatory roles in different cellular compartments including the mitochondria, as a phosphatase binding protein; in the cytoplasm, as a chaperone-like protein; and in the nucleus, as a transcriptional regulator through binding to RPB5 and RNA polymerase II (polII). This chapter focuses on the role URI plays in transcriptional regulation in the prostate cell. In prostate cells, URI is tightly bound to another prefoldin-like protein called UXT, a known androgen receptor (AR) cofactor. Part of a multiprotein complex, URI and UXT act as transcriptional repressors, and URI regulates KAP1 through PP2A phosphatase activity. The discovery of the interaction of URI and UXT with KAP1, AR, and PP2A, as well as the numerous interactions between URI and components of the R2TP/prefoldin-like complex, RPB5, and nuclear proteins involved in DNA damage response, chromatin remodeling and gene transcription, reveal a pleiotropic effect of the URI/UXT complex on nuclear processes. The mechanisms by which URI/UXT affect transcription, chromatin structure and regulation, and genome stability, remain to be elucidated but will be of fundamental importance considering the many processes affected by alterations of URI/UXT and other prefoldins and prefoldin-like proteins.

Chronic low-dose pro-oxidant treatment stimulates transcriptional activity of telomeric retroelements and increases telomere length in Drosophila

It has been proposed that oxidative stress, elicited by high levels of reactive oxygen species, accelerates telomere shortening by erosion of telomeric DNA repeats. While most eukaryotes counteract telomere shortening by telomerase-driven addition of these repeats, telomeric loss in Drosophila is compensated by retrotransposition of the telomeric retroelements HeT-A, TART and TAHRE to chromosome ends. In this study we tested the effect of chronic exposure of flies to non-/sub-lethal doses of paraquat, which is a redox cycling compound widely used to induce oxidative stress in various experimental paradigms including telomere length analyses. Indeed, chronic paraquat exposure for five generations resulted in elevated transcriptional activity of both telomeric and non-telomeric transposable elements, and extended telomeric length in the tested fly lines. We propose that low oxidative stress leads to increased telomere length within Drosophila populations. For a mechanistic understanding of the observed phenomenon we discuss two scenarios: adaption, acting through a direct stimulation of telomere extension, or positive selection favoring individuals with longer telomeres within the population.

HERV-K activation is strictly required to sustain CD133+ melanoma cells with stemness features

BACKGROUND

Melanoma is a heterogeneous tumor in which phenotype-switching and CD133 marker have been associated with metastasis promotion and chemotherapy resistance. CD133 positive (CD133+) subpopulation has also been suggested as putative cancer stem cell (CSC) of melanoma tumor. Human endogenous retrovirus type K (HERV-K) has been described to be aberrantly activated during melanoma progression and implicated in the etiopathogenesis of disease. Earlier, we reported that stress-induced HERV-K activation promotes cell malignant transformation and reduces the immunogenicity of melanoma cells. Herein, we investigated the correlation between HERV-K and the CD133+ melanoma cells during microenvironmental modifications.

METHODS

TVM-A12 cell line, isolated in our laboratory from a primary human melanoma lesion, and other commercial melanoma cell lines (G-361, WM-115, WM-266-4 and A375) were grown and maintained in the standard and stem cell media. RNA interference, Real-time PCR, flow cytometry analysis, self-renewal and migration/invasion assays were performed to characterize cell behavior and HERV-K expression.

RESULTS

Melanoma cells, exposed to stem cell media, undergo phenotype-switching and expansion of CD133+ melanoma cells, concomitantly promoted by HERV-K activation. Notably, the sorted CD133+ subpopulation showed stemness features, characterized by higher self-renewal ability, embryonic genes expression, migration and invasion capacities compared to the parental cell line. RNA interference-mediated downregulation experiments showed that HERV-K has a decisive role to expand and maintain the CD133+ melanoma subpopulation during microenvironmental modifications. Similarly, non nucleoside reverse transcriptase inhibitors (NNRTIs) efavirenz and nevirapine were effective to restrain the activation of HERV-K in melanoma cells, to antagonize CD133+ subpopulation expansion and to induce selective high level apoptosis in CD133+ cells.

CONCLUSIONS

HERV-K activation promotes melanoma cells phenotype-switching and is strictly required to expand and maintain the CD133+ melanoma cells with stemness features in response to microenvironmental modifications.

A survey of endogenous retrovirus (ERV) sequences in the vicinity of multiple sclerosis (MS)-associated single nucleotide polymorphisms (SNPs)

Although multiple sclerosis (MS) is one of the most common central nervous system diseases in young adults, little is known about its etiology. Several human endogenous retroviruses (ERVs) are considered to play a role in MS. We are interested in which ERVs can be identified in the vicinity of MS associated genetic marker to find potential initiators of MS. We analysed the chromosomal regions surrounding 58 single nucleotide polymorphisms (SNPs) that are associated with MS identified in one of the last major genome wide association studies. We scanned these regions for putative endogenous retrovirus sequences with large open reading frames (ORFs). We observed that more retrovirus-related putative ORFs exist in the relatively close vicinity of SNP marker indices in multiple sclerosis compared to control SNPs. We found very high homologies to HERV-K, HCML-ARV, XMRV, Galidia ERV, HERV-H/env62 and XMRV-like mouse endogenous retrovirus mERV-XL. The associated genes (CYP27B1, CD6, CD58, MPV17L2, IL12RB1, CXCR5, PTGER4, TAGAP, TYK2, ICAM3, CD86, GALC, GPR65 as well as the HLA DRB1*1501) are mainly involved in the immune system, but also in vitamin D regulation. The most frequently detected ERV sequences are related to the multiple sclerosis-associated retrovirus, the human immunodeficiency virus 1, HERV-K, and the Simian foamy virus. Our data shows that there is a relation between MS associated SNPs and the number of retroviral elements compared to control. Our data identifies new ERV sequences that have not been associated with MS, so far.

Bacterial Group II Introns: Identification and Mobility Assay

Group II introns are large catalytic RNAs and mobile retroelements that encode a reverse transcriptase. Here, we provide methods for their identification in bacterial genomes and further analysis of their splicing and mobility capacities.