Active human retrotransposons: variation and disease

Transposable Elements: Major Players in Shaping Genomic and Evolutionary Patterns

Transposable elements (TEs) are ubiquitous genetic elements, able to jump from one location of the genome to another, in all organisms. For this reason, on the one hand, TEs can induce deleterious mutations, causing dysfunction, disease and even lethality in individuals. On the other hand, TEs can increase genetic variability, making populations better equipped to respond adaptively to environmental change. To counteract the deleterious effects of TEs, organisms have evolved strategies to avoid their activation. However, their mobilization does occur. Usually, TEs are maintained silent through several mechanisms, but they can be reactivated during certain developmental windows. Moreover, TEs can become de-repressed because of drastic changes in the external environment. Here, we describe the ‘double life’ of TEs, being both ‘parasites’ and ‘symbionts’ of the genome. We also argue that the transposition of TEs contributes to two important evolutionary processes: the temporal dynamic of evolution and the induction of genetic variability. Finally, we discuss how the interplay between two TE-dependent phenomena, insertional mutagenesis and epigenetic plasticity, plays a role in the process of evolution.

Retrotransposons in the Mammalian Male Germline

Retrotransposons are a subset of DNA sequences that constitute a large part of the mammalian genome. They can translocate autonomously or non-autonomously, potentially jeopardizing the heritable germline genome. Retrotransposons coevolved with the host genome, and the germline is the prominent battlefield between retrotransposons and the host genome to maximize their mutual fitness. Host genomes have developed various mechanisms to suppress and control retrotransposons, including DNA methylation, histone modifications, and Piwi-interacting RNA (piRNA), for their own benefit. Thus, rapidly evolved retrotransposons often acquire positive functions, including gene regulation within the germline, conferring reproductive fitness in a species over the course of evolution. The male germline serves as an ideal model to examine the regulation and evolution of retrotransposons, resulting in genomic co-evolution with the host genome. In this review, we summarize and discuss the regulatory mechanisms of retrotransposons, stage-by-stage, during male germ cell development, with a particular focus on mice as an extensively studied mammalian model, highlighting suppression mechanisms and emerging functions of retrotransposons in the male germline.

Intracellular Reverse Transcription of Pfizer BioNTech COVID-19 mRNA Vaccine BNT162b2 In Vitro in Human Liver Cell Line

Preclinical studies of COVID-19 mRNA vaccine BNT162b2, developed by Pfizer and BioNTech, showed reversible hepatic effects in animals that received the BNT162b2 injection. Furthermore, a recent study showed that SARS-CoV-2 RNA can be reverse-transcribed and integrated into the genome of human cells. In this study, we investigated the effect of BNT162b2 on the human liver cell line Huh7 in vitro. Huh7 cells were exposed to BNT162b2, and quantitative PCR was performed on RNA extracted from the cells. We detected high levels of BNT162b2 in Huh7 cells and changes in gene expression of long interspersed nuclear element-1 (LINE-1), which is an endogenous reverse transcriptase. Immunohistochemistry using antibody binding to LINE-1 open reading frame-1 RNA-binding protein (ORFp1) on Huh7 cells treated with BNT162b2 indicated increased nucleus distribution of LINE-1. PCR on genomic DNA of Huh7 cells exposed to BNT162b2 amplified the DNA sequence unique to BNT162b2. Our results indicate a fast up-take of BNT162b2 into human liver cell line Huh7, leading to changes in LINE-1 expression and distribution. We also show that BNT162b2 mRNA is reverse transcribed intracellularly into DNA in as fast as 6 h upon BNT162b2 exposure.

A framework to score the effects of structural variants in health and disease

While technological advances improved the identification of structural variants (SVs) in the human genome, their interpretation remains challenging. Several methods utilize individual mechanistic principles like the deletion of coding sequence or 3D genome architecture disruptions. However, a comprehensive tool using the broad spectrum of available annotations is missing. Here, we describe CADD-SV, a method to retrieve and integrate a wide set of annotations to predict the effects of SVs. Previously, supervised learning approaches were limited due to a small number and biased set of annotated pathogenic or benign SVs. We overcome this problem by using a surrogate training-objective, the Combined Annotation Dependent Depletion (CADD) of functional variants. We use human and chimpanzee derived SVs as proxy-neutral and contrast them with matched simulated variants as proxy-deleterious, an approach that has proven powerful for short sequence variants. Our tool computes summary statistics over diverse variant annotations and uses random forest models to prioritize deleterious structural variants. The resulting CADD-SV scores correlate with known pathogenic and rare population variants. We further show that we can prioritize somatic cancer variants as well as noncoding variants known to affect gene expression. We provide a website and offline-scoring tool for easy application of CADD-SV.

FNC: An Advanced Anticancer Therapeutic or Just an Underdog?

Azvudine (FNC) is a novel cytidine analogue that has both antiviral and anticancer activities. This minireview focuses on its underlying molecular mechanisms of suppressing viral life cycle and cancer cell growth and discusses applications of this nucleoside drug for advanced therapy of tumors and malignant blood diseases. FNC inhibits positive-stand RNA viruses, like HCV, EV, SARS-COV-2, HBV, and retroviruses, including HIV, by suppressing their RNA-dependent polymerase enzymes. It may also inhibit such enzyme (reverse transcriptase) in the human retrotransposons, including human endogenous retroviruses (HERVs). As the activation of retrotransposons can be the major factor of ongoing cancer genome instability and consequently higher aggressiveness of tumors, FNC has a potential to increase the efficacy of multiple anticancer therapies. Furthermore, FNC also showed other aspects of anticancer activity by inhibiting adhesion, migration, invasion, and proliferation of malignant cells. It was also reported to be involved in cell cycle arrest and apoptosis, thereby inhibiting the progression of cancer through different pathways. To the date, the grounds of FNC effects on cancer cells are not fully understood and hence additional studies are needed for better understanding molecular mechanisms of its anticancer activities to support its medical use in oncology.

A Variety of Nucleic Acid Species Are Sensed by cGAS, Implications for Its Diverse Functions

Cyclic GMP-AMP synthase (cGAS) recognizes double-stranded DNA (dsDNA) derived from invading pathogens and induces an interferon response via activation of the key downstream adaptor protein stimulator of interferon genes (STING). This is the most classic biological function of the cGAS-STING signaling pathway and is critical for preventing pathogenic microorganism invasion. In addition, cGAS can interact with various types of nucleic acids, including cDNA, DNA : RNA hybrids, and circular RNA, to contribute to a diverse set of biological functions. An increasing number of studies have revealed an important relationship between the cGAS-STING signaling pathway and autophagy, cellular senescence, antitumor immunity, inflammation, and autoimmune diseases. This review details the mechanism of action of cGAS as it interacts with different types of nucleic acids, its rich biological functions, and the potential for targeting this pathway to treat various diseases.

Deletion and replacement of long genomic sequences using prime editing

Genomic insertions, duplications and insertion/deletions (indels), which account for ~14% of human pathogenic mutations, cannot be accurately or efficiently corrected by current gene-editing methods, especially those that involve larger alterations (>100 base pairs (bp)). Here, we optimize prime editing (PE) tools for creating precise genomic deletions and direct the replacement of a genomic fragment ranging from ~1 kilobases (kb) to ~10 kb with a desired sequence (up to 60 bp) in the absence of an exogenous DNA template. By conjugating Cas9 nuclease to reverse transcriptase (PE-Cas9) and combining it with two PE guide RNAs (pegRNAs) targeting complementary DNA strands, we achieve precise and specific deletion and repair of target sequences via using this PE-Cas9-based deletion and repair (PEDAR) method. PEDAR outperformed other genome-editing methods in a reporter system and at endogenous loci, efficiently creating large and precise genomic alterations. In a mouse model of tyrosinemia, PEDAR removed a 1.38-kb pathogenic insertion within the Fah gene and precisely repaired the deletion junction to restore FAH expression in liver.

Detecting Bacterial-Human Lateral Gene Transfer in Chronic Lymphocytic Leukemia

Chronic lymphocytic leukemia (CLL) is a very common and mostly incurable B-cell malignancy. Recent studies revealed high interpatient mutational heterogeneity and worsened therapy response and survival of patients with complex genomic aberrations. In line with this, a better understanding of the underlying mechanisms of specific genetic aberrations would reveal new prognostic factors and possible therapeutic targets. It is known that chromosomal rearrangements including DNA insertions often play a role during carcinogenesis. Recently it was reported that bacteria (microbiome)–human lateral gene transfer occurs in somatic cells and is enriched in cancer samples. To further investigate this mechanism in CLL, we analyzed paired-end RNA sequencing data of 45 CLL patients and 9 healthy donors, in which we particularly searched for bacterial DNA integrations into the human somatic genome. Applying the Burrows–Wheeler aligner (BWA) first on a human genome and then on bacterial genome references, we differentiated between sequencing reads mapping to the human genome, to the microbiome or to bacterial integrations into the human genome. Our results indicate that CLL samples featured bacterial DNA integrations more frequently (approx. two-fold) compared to normal samples, which corroborates the latest findings in other cancer entities. Moreover, we determined common integration sites and recurrent integrated bacterial transcripts. Finally, we investigated the contribution of bacterial integrations to oncogenesis and disease progression.

The non-LTR retrotransposons of Entamoeba histolytica: genomic organization and biology

Genome sequence analysis of Entamoeba species revealed various classes of transposable elements. While E. histolytica and E. dispar are rich in non-long terminal repeat (LTR) retrotransposons, E. invadens contains predominantly DNA transposons. Non-LTR retrotransposons of E. histolytica constitute three families of long interspersed nuclear elements (LINEs), and their short, nonautonomous partners, SINEs. They occupy ~ 11% of the genome. The EhLINE1/EhSINE1 family is the most abundant and best studied. EhLINE1 is 4.8 kb, with two ORFs that encode functions needed for retrotransposition. ORF1 codes for the nucleic acid-binding protein, and ORF2 has domains for reverse transcriptase (RT) and endonuclease (EN). Most copies of EhLINEs lack complete ORFs. ORF1p is expressed constitutively, but ORF2p is not detected. Retrotransposition could be demonstrated upon ectopic over expression of ORF2p, showing that retrotransposition machinery is functional. The newly retrotransposed sequences showed a high degree of recombination. In transcriptomic analysis, RNA-Seq reads were mapped to individual EhLINE1 copies. Although full-length copies were transcribed, no full-length 4.8 kb transcripts were seen. Rather, sense transcripts mapped to ORF1, RT and EN domains. Intriguingly, there was strong antisense transcription almost exclusively from the RT domain. These unique features of EhLINE1 could serve to attenuate retrotransposition in E. histolytica.

Emerging roles and potential application of PIWI-interacting RNA in urological tumors

The piRNA (PIWI-interacting RNA) is P-Element induced wimpy testis (PIWI)-interacting RNA which is a small molecule, non-coding RNA with a length of 24-32nt. It was originally found in germ cells and is considered a regulator of germ cell function. It can interact with PIWI protein, a member of the Argonaute family, and play a role in the regulation of gene transcription and epigenetic silencing of transposable factors in the nucleus. More and more studies have shown that piRNAs are abnormally expressed in a variety of cancer tissues and patient fluids, and may become diagnostic tools, therapeutic targets, staging markers, and prognostic evaluation tools for cancer. This article reviews the recent research on piRNA and summarizes the structural characteristics, production mechanism, applications, and its role in urological tumors, to provide a reference value for piRNA to regulate urological tumors.

Comparative analysis of transposable elements provides insights into genome evolution in the genus Camelus

BACKGROUND

Transposable elements (TEs) are common features in eukaryotic genomes that are known to affect genome evolution critically and to play roles in gene regulation. Vertebrate genomes are dominated by TEs, which can reach copy numbers in the hundreds of thousands. To date, details regarding the presence and characteristics of TEs in camelid genomes have not been made available.

RESULTS

We conducted a genome-wide comparative analysis of camelid TEs, focusing on the identification of TEs and elucidation of transposition histories in four species: Camelus dromedarius, C. bactrianus, C. ferus, and Vicugna pacos. Our TE library was created using both de novo structure-based and homology-based searching strategies ( https://github.com/kacst-bioinfo-lab/TE_ideintification_pipeline ). Annotation results indicated a similar proportion of each genomes comprising TEs (35-36%). Class I LTR retrotransposons comprised 16-20% of genomes, and mostly consisted of the endogenous retroviruses (ERVs) groups ERVL, ERVL-MaLR, ERV_classI, and ERV_classII. Non-LTR elements comprised about 12% of genomes and consisted of SINEs (MIRs) and the LINE superfamilies LINE1, LINE2, L3/CR1, and RTE clades. Least represented were the Class II DNA transposons (2%), consisting of hAT-Charlie, TcMar-Tigger, and Helitron elements and comprising about 1-2% of each genome.

CONCLUSIONS

The findings of the present study revealed that the distribution of transposable elements across camelid genomes is approximately similar. This investigation presents a characterization of TE content in four camelid to contribute to developing a better understanding of camelid genome architecture and evolution.

Systems Biology to Understand and Regulate Human Retroviral Proinflammatory Response

The majority of human genome are non-coding genes. Recent research have revealed that about half of these genome sequences make up of transposable elements (TEs). A branch of these belong to the endogenous retroviruses (ERVs), which are germline viral infection that occurred over millions of years ago. They are generally harmless as evolutionary mutations have made them unable to produce viral agents and are mostly epigenetically silenced. Nevertheless, ERVs are able to express by still unknown mechanisms and recent evidences have shown links between ERVs and major proinflammatory diseases and cancers. The major challenge is to elucidate a detailed mechanistic understanding between them, so that novel therapeutic approaches can be explored. Here, we provide a brief overview of TEs, human ERVs and their links to microbiome, innate immune response, proinflammatory diseases and cancer. Finally, we recommend the employment of systems biology approaches for future HERV research.

Estimating DNA methylation potential energy landscapes from nanopore sequencing data

High-throughput third-generation nanopore sequencing devices have enormous potential for simultaneously observing epigenetic modifications in human cells over large regions of the genome. However, signals generated by these devices are subject to considerable noise that can lead to unsatisfactory detection performance and hamper downstream analysis. Here we develop a statistical method, CpelNano, for the quantification and analysis of 5mC methylation landscapes using nanopore data. CpelNano takes into account nanopore noise by means of a hidden Markov model (HMM) in which the true but unknown (“hidden”) methylation state is modeled through an Ising probability distribution that is consistent with methylation means and pairwise correlations, whereas nanopore current signals constitute the observed state. It then estimates the associated methylation potential energy function by employing the expectation-maximization (EM) algorithm and performs differential methylation analysis via permutation-based hypothesis testing. Using simulations and analysis of published data obtained from three human cell lines (GM12878, MCF-10A, and MDA-MB-231), we show that CpelNano can faithfully estimate DNA methylation potential energy landscapes, substantially improving current methods and leading to a powerful tool for the modeling and analysis of epigenetic landscapes using nanopore sequencing data.

Transposable elements that have recently been mobile in the human genome

Background

Transposable elements (TE) comprise nearly half of the human genome and their insertions have profound effects to human genetic diversification and as well as disease. Despite their abovementioned significance, there is no consensus on the TE subfamilies that remain active in the human genome. In this study, we therefore developed a novel statistical test for recently mobile subfamilies (RMSs), based on patterns of overlap with > 100,000 polymorphic indels.

Results

Our analysis produced a catalogue of 20 high-confidence RMSs, which excludes many false positives in public databases. Intriguingly though, it includes HERV-K, an LTR subfamily previously thought to be extinct. The RMS catalogue is strongly enriched for contributions to germline genetic disorders (P = 1.1e-10), and thus constitutes a valuable resource for diagnosing disorders of unknown aetiology using targeted TE-insertion screens. Remarkably, RMSs are also highly enriched for somatic insertions in diverse cancers (P = 2.8e-17), thus indicating strong correlations between germline and somatic TE mobility. Using CRISPR/Cas9 deletion, we show that an RMS-derived polymorphic TE insertion increased the expression of RPL17, a gene associated with lower survival in liver cancer. More broadly, polymorphic TE insertions from RMSs were enriched near genes with allele-specific expression, suggesting widespread effects on gene regulation.

Conclusions

By using a novel statistical test we have defined a catalogue of 20 recently mobile transposable element subfamilies. We illustrate the gene regulatory potential of RMS-derived polymorphic TE insertions, using CRISPR/Cas9 deletion in vitro on a specific candidate, as well as by genome wide analysis of allele-specific expression. Our study presents novel insights into TE mobility and regulatory potential and provides a key resource for human disease genetics and population history studies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08085-0.

BST2 suppresses LINE-1 retrotransposition by reducing the promoter activity of LINE-1 5'-UTR

Endogenous retrotransposons are considered the “molecular fossils” of ancient retroviral insertions. Several studies have indicated that host factors restrict both retroviruses and retrotransposons through different mechanisms. Type 1 long interspersed elements (LINE-1 or L1) are the only active retroelements that can replicate autonomously in the human genome. A recent study reported that LINE-1 retrotransposition is potently suppressed by BST2, a host restriction factor that prevents viral release mainly by physically tethering enveloped virions (such as HIV) to the surface of producer cells. However, no endoplasmic membrane structure has been associated with LINE-1 replication, suggesting that BST2 may utilize a distinct mechanism to suppress LINE-1. In this study, we showed that BST2 is a potent LINE-1 suppressor. Further investigations suggested that BST2 reduces the promoter activity of LINE-1 5′ untranslated region (UTR) and lowers the levels of LINE-1 RNA, proteins, and events during LINE-1 retrotransposition. Surprisingly, although BST2 apparently uses different mechanisms against HIV and LINE-1, two membrane-associated domains that are essential for BST2-mediated HIV tethering also proved important for BST2-induced inhibition of LINE-1 5′ UTR. Additionally, by suppressing LINE-1, BST2 prevented LINE-1-induced genomic DNA damage and innate immune activation. Taken together, our data uncovered the mechanism of BST2-mediated LINE-1 suppression and revealed new roles of BST2 as a promoter regulator, genome stabilizer, and innate immune suppressor.

IMPORTANCE BST2 is a potent antiviral protein that suppresses the release of several enveloped viruses, mainly by tethering the envelope of newly synthesized virions and restraining them on the surface of producer cells. In mammalian cells, there are numerous DNA elements replicating through reverse transcription, among which LINE-1 is the only retroelement that can replicate autonomously. Although LINE-1 retrotransposition does not involve the participation of a membrane structure, BST2 has been reported as an efficient LINE-1 suppressor, suggesting a different mechanism for BST2-mediated LINE-1 inhibition and a new function for BST2 itself. We found that BST2 specifically represses the promoter activity of LINE-1 5′ UTR, resulting in decreased levels of LINE-1 transcription, translation, and subsequent retrotransposition. Additionally, by suppressing LINE-1 activity, BST2 maintains genome stability and regulates innate immune activation. These findings expand our understanding of BST2 and its biological significance.

Differential expression of transposable elements in the medaka melanoma model

Malignant melanoma incidence is rising worldwide. Its treatment in an advanced state is difficult, and the prognosis of this severe disease is still very poor. One major source of these difficulties is the high rate of metastasis and increased genomic instability leading to a high mutation rate and the development of resistance against therapeutic approaches. Here we investigate as one source of genomic instability the contribution of activation of transposable elements (TEs) within the tumor. We used the well-established medaka melanoma model and RNA-sequencing to investigate the differential expression of TEs in wildtype and transgenic fish carrying melanoma. We constructed a medaka-specific TE sequence library and identified TE sequences that were specifically upregulated in tumors. Validation by qRT- PCR confirmed a specific upregulation of a LINE and an LTR element in malignant melanomas of transgenic fish.

Alu RNA and their roles in human disease states

Alu RNA are implicated in the poor prognosis of several human disease states. These RNA are transcription products of primate specific transposable elements called Alu elements. These elements are extremely abundant, comprising over 10% of the human genome, and 100 to 1000 cytoplasmic copies of Alu RNA per cell. Alu RNA do not have a single universal functional role aside from selfish self-propagation. Despite this, Alu RNA have been found to operate in a diverse set of translational and transcriptional mechanisms. This review will focus on the current knowledge of Alu RNA involved in human disease states and known mechanisms of action. Examples of Alu RNA that are transcribed in a variety of contexts such as introns, mature mRNA, and non-coding transcripts will be discussed. Past and present challenges in studying Alu RNA, and the future directions of Alu RNA in basic and clinical research will also be examined.

Transposable elements as new players in neurodegenerative diseases

Neurodegenerative diseases (NDs), including the most prevalent Alzheimer's disease and Parkinson disease, share common pathological features. Despite decades of gene-centric approaches, the molecular mechanisms underlying these diseases remain widely elusive. In recent years, transposable elements (TEs), long considered 'junk' DNA, have gained growing interest as pathogenic players in NDs. Age is the major risk factor for most NDs, and several repressive mechanisms of TEs, such as heterochromatinization, fail with age. Indeed, heterochromatin relaxation leading to TE derepression has been reported in various models of neurodegeneration and NDs. There is also evidence that certain pathogenic proteins involved in NDs (e.g., tau, TDP-43) may control the expression of TEs. The deleterious consequences of TE activation are not well known but they could include DNA damage and genomic instability, altered host gene expression, and/or neuroinflammation, which are common hallmarks of neurodegeneration and aging. TEs might thus represent an overlooked pathogenic culprit for both brain aging and neurodegeneration. Certain pathological effects of TEs might be prevented by inhibiting their activity, pointing to TEs as novel targets for neuroprotection.

Open Reading Frame 1 Protein of the Human Long Interspersed Nuclear Element 1 Retrotransposon Binds Multiple Equivalents of Lead

The human long interspersed nuclear element 1 (LINE1) has been implicated in numerous diseases and has been suggested to play a significant role in genetic evolution. Open reading frame 1 protein (ORF1p) is one of the two proteins encoded in this self-replicating mobile genetic element, both of which are essential for retrotransposition. The structure of the three-stranded coiled-coil domain of ORF1p was recently solved and showed the presence of tris-cysteine layers in the interior of the coiled-coil that could function as metal binding sites. Here, we demonstrate that ORF1p binds Pb(II). We designed a model peptide, GRCSL16CL23C, to mimic two of the ORF1p Cys3 layers and crystallized the peptide both as the apo-form and in the presence of Pb(II). Structural comparison of the ORF1p with apo-(GRCSL16CL23C)3 shows very similar Cys3 layers, preorganized for Pb(II) binding. We propose that exposure to heavy metals, such as lead, could influence directly the structural parameters of ORF1p and thus impact the overall LINE1 retrotransposition frequency, directly relating heavy metal exposure to genetic modification.

Characterization of full-length LINE-1 insertions in 154 genomes

Long interspersed nuclear elements (LINEs) are retrotransposons that contribute to genetic variation in the human genome. LINE-1 elements in larger-scale studies are challenging to identify using sequencing technologies due to cost and scalability. We developed an approach using optical mapping for detection of full-length LINE-1 insertions and 10× sequencing for confirmation. We found 51 true positive full-length LINE-1 insertions, of which 4 are novel insertions, in NA12878. Repeating our analysis on a larger sample set representing 26 populations, we identified 329 full-length LINE-1 elements, of which 123 are novel. 24.8% of these 329 LINE-1 insertions were shared amongst all 5 superpopulations (AFR, AMR, EUR, EAS, SAS). The African superpopulation has a higher percentage of population-specific LINE-1 insertions than any other superpopulation. These data indicate that our approach can provide high-speed, cost-effective, and increased accuracy for LINE-1 detection. These data also provide an insight into variations of LINE-1 elements between different populations.

RNA m6A modification orchestrates a LINE-1-host interaction that facilitates retrotransposition and contributes to long gene vulnerability

The molecular basis underlying the interaction between retrotransposable elements (RTEs) and the human genome remains poorly understood. Here, we profiled N⁶-methyladenosine (m⁶A) deposition on nascent RNAs in human cells by developing a new method MINT-Seq, which revealed that many classes of RTE RNAs, particularly intronic LINE-1s (L1s), are strongly methylated. These m⁶A-marked intronic L1s (MILs) are evolutionarily young, sense-oriented to hosting genes, and are bound by a dozen RNA binding proteins (RBPs) that are putative novel readers of m⁶A-modified RNAs, including a nuclear matrix protein SAFB. Notably, m⁶A positively controls the expression of both autonomous L1s and co-transcribed L1 relics, promoting L1 retrotransposition. We showed that MILs preferentially reside in long genes with critical roles in DNA damage repair and sometimes in L1 suppression per se, where they act as transcriptional "roadblocks" to impede the hosting gene expression, revealing a novel host-weakening strategy by the L1s. In counteraction, the host uses the SAFB reader complex to bind m⁶A-L1s to reduce their levels, and to safeguard hosting gene transcription. Remarkably, our analysis identified thousands of MILs in multiple human fetal tissues, enlisting them as a novel category of cell-type-specific regulatory elements that often compromise transcription of long genes and confer their vulnerability in neurodevelopmental disorders. We propose that this m⁶A-orchestrated L1-host interaction plays widespread roles in gene regulation, genome integrity, human development and diseases.

Distinct Retrotransposon Evolution Profile in the Genome of Rabbit (Oryctolagus cuniculus)

Although the rabbit genome has already been annotated, it is mobilome remains largely unknown. Here, multiple pipelines were used to de novo mine and annotate the mobilome in rabbit. Four families and 19 subfamilies of LINE1s, two families and nine subfamilies of SINEs, and 12 ERV families were defined in rabbit based on sequence identity, structural organization, and phylogenetic tree. The analysis of insertion age and polymerase chain reaction suggests that a number of families are very young and may remain active, such as L1B, L1D, OcuSINEA, and OcuERV1. RepeatMasker annotation revealed a distinct transposable element landscape within the genome, with approximately two million copies of SINEs, representing the greatest proportion of the genome (19.61%), followed by LINEs (15.44%), and LTRs (4.11%), respectively, considerably different from most other mammal mobilomes except hedgehog and tree shrew, in which LINEs have the highest proportion. Furthermore, a very high rate of insertion polymorphisms (>85%) for the youngest subfamily (OcuSINEA1) was identified by polymerase chain reaction. The majority of retrotransposon insertions overlapped with protein-coding regions (>80%) and lncRNA (90%) genes. Genomic distribution bias was observed for retrotransposons, with those immediately upstream (-1 kb) and downstream (1 kb) of genes significantly depleted. Local GC content in 50-kb widows had significantly negative correlations with LINE (rs=-0.996) and LTR (rs=-0.829) insertions. The current study revealed a distinct mobilome landscape in rabbit, which will assist in the elucidation of the evolution of the genome of lagomorphs, and even other mammals.

The RNA editing enzyme ADAR2 restricts L1 mobility

Adenosine deaminases acting on RNA (ADARs) are enzymes that convert adenosines to inosines in double-stranded RNAs (RNA editing A-to-I). ADAR1 and ADAR2 were previously reported as HIV-1 proviral factors. The aim of this study was to investigate the composition of the ADAR2 ribonucleoprotein complex during HIV-1 expression. By using a dual-tag affinity purification procedure in cells expressing HIV-1 followed by mass spectrometry analysis, we identified 10 non-ribosomal ADAR2-interacting factors. A significant fraction of these proteins was previously found associated to the Long INterspersed Element 1 (LINE1 or L1) ribonucleoparticles and to regulate the life cycle of L1 retrotransposons. Considering that we previously demonstrated that ADAR1 is an inhibitor of LINE-1 retrotransposon activity, we investigated whether also ADAR2 played a similar function. To reach this goal, we performed specific cell culture retrotransposition assays in cells overexpressing or ablated for ADAR2. These experiments unveil a novel function of ADAR2 as suppressor of L1 retrotransposition. Furthermore, we showed that ADAR2 binds the basal L1 RNP complex.

Overall, these data support the role of ADAR2 as regulator of L1 life cycle.

SINE jumping contributes to large-scale polymorphisms in the pig genomes

BACKGROUND

Molecular markers based on retrotransposon insertion polymorphisms (RIPs) have been developed and are widely used in plants and animals. Short interspersed nuclear elements (SINEs) exert wide impacts on gene activity and even on phenotypes. However, SINE RIP profiles in livestock remain largely unknown, and not be revealed in pigs.

RESULTS

Our data revealed that SINEA1 displayed the most polymorphic insertions (22.5 % intragenic and 26.5 % intergenic), followed by SINEA2 (10.5 % intragenic and 9 % intergenic) and SINEA3 (12.5 % intragenic and 5.0 % intergenic). We developed a genome-wide SINE RIP mining protocol and obtained a large number of SINE RIPs (36,284), with over 80 % accuracy and an even distribution in chromosomes (14.5/Mb), and 74.34 % of SINE RIPs generated by SINEA1 element. Over 65 % of pig SINE RIPs overlap with genes, most of them (> 95 %) are in introns. Overall, about one forth (23.09 %) of the total genes contain SINE RIPs. Significant biases of SINE RIPs in the transcripts of protein coding genes were observed. Nearly half of the RIPs are common in these pig breeds. Sixteen SINE RIPs were applied for population genetic analysis in 23 pig breeds, the phylogeny tree and cluster analysis were generally consistent with the geographical distributions of native pig breeds in China.

CONCLUSIONS

Our analysis revealed that SINEA1-3 elements, particularly SINEA1, are high polymorphic across different pig breeds, and generate large-scale structural variations in the pig genomes. And over 35,000 SINE RIP markers were obtained. These data indicate that young SINE elements play important roles in creating new genetic variations and shaping the evolution of pig genome, and also provide strong evidences to support the great potential of SINE RIPs as genetic markers, which can be used for population genetic analysis and quantitative trait locus (QTL) mapping in pig.

Expression Quantitative Trait Loci (eQTLs) Associated with Retrotransposons Demonstrate their Modulatory Effect on the Transcriptome

Transposable elements (TEs) are repetitive elements that belong to a variety of functional classes and have an important role in shaping genome evolution. Around 50% of the human genome contains TEs, and they have been termed the "dark matter" of the genome because relatively little is known about their function. While TEs have been shown to participate in aberrant gene regulation and the pathogenesis of diseases, only a few studies have explored the systemic effect of TEs on gene expression. In the present study, we analysed whole genome sequences and blood whole transcriptome data from 570 individuals within the Parkinson's Progressive Markers Initiative (PPMI) cohort to identify expression quantitative trait loci (eQTL) regulating genome-wide gene expression associated with TEs. We identified 2132 reference TEs that were polymorphic for their presence or absence in our study cohort. The presence or absence of the TE element could change the expression of the gene or gene clusters from zero to tens of thousands of copies of RNA. The main finding is that many TEs possess very strong regulatory effects, and they have the potential to modulate large genetic networks with hundreds of target genes over the genome. We illustrate the plethora of regulatory mechanisms using examples of their action at the HLA gene cluster and data showing different TEs' convergence to modulate WFS1 gene expression. In conclusion, the presence or absence of polymorphisms of TEs has an eminent genome-wide regulatory function with large effect size at the level of the whole transcriptome. The role of TEs in explaining, in part, the missing heritability for complex traits is convincing and should be considered.

A porcine brain-wide RNA editing landscape

Adenosine-to-inosine (A-to-I) RNA editing, catalyzed by ADAR enzymes, is an essential post-transcriptional modification. Although hundreds of thousands of RNA editing sites have been reported in mammals, brain-wide analysis of the RNA editing in the mammalian brain remains rare. Here, a genome-wide RNA-editing investigation is performed in 119 samples, representing 30 anatomically defined subregions in the pig brain. We identify a total of 682,037 A-to-I RNA editing sites of which 97% are not identified before. Within the pig brain, cerebellum and olfactory bulb are regions with most edited transcripts. The editing level of sites residing in protein-coding regions are similar across brain regions, whereas region-distinct editing is observed in repetitive sequences. Highly edited conserved recoding events in pig and human brain are found in neurotransmitter receptors, demonstrating the evolutionary importance of RNA editing in neurotransmission functions. Although potential data biases caused by age, sex or health status are not considered, this study provides a rich resource to better understand the evolutionary importance of post-transcriptional RNA editing.

Huang et al performed a genome-wide RNA editing investigation in the porcine brain in which they found over 680,000 A-to-I RNA editing sites. They identified conserved recoding events between pig and human brains thus providing an extensive resource to aid our understanding of the evolutionary importance of post-transcriptional RNA editing.

The methylome of a human polar body reflects that of its sibling oocyte and its aberrance may indicate poor embryo development

STUDY QUESTION

Is it possible to evaluate the methylome of individual oocytes to investigate the DNA methylome alterations in metaphase II (MII) oocytes with reduced embryo developmental potential?

SUMMARY ANSWER

The DNA methylome of each human first polar body (PB1) closely mirrored that of its sibling MII oocyte; hypermethylated long interspersed nuclear element (LINE) and long terminal repeats (LTRs) and methylation aberrations in PB1 promoter regions may indicate poor embryo development.

WHAT IS KNOWN ALREADY

The developmental potential of an embryo is determined by the oocyte's developmental competence, and the PB1 is a good substitute to examine the chromosomal status of the corresponding oocyte. However, DNA methylation, a key epigenetic modification, also regulates gene expression and embryo development.

STUDY DESIGN, SIZE, DURATION

Twelve pairs of PB1s and sibling MII oocytes were biopsied and sequenced to compare their methylomes. To further investigate the methylome of PB1s and the potential epigenetic factors that may affect oocyte quality, MII oocytes (n = 74) were fertilized through ICSI, while PB1s were biopsied and profiled to measure DNA methylation. The corresponding embryos were further cultured to track their development potential. The oocytes and sperm samples used in this study were donated by healthy volunteers with signed informed consent.

PARTICIPANTS/MATERIALS, SETTING, METHODS

Single-cell methylome sequencing was applied to obtain the DNA methylation profiles of PB1s and oocytes. The DNA methylome of PB1s was compared between the respective group of oocytes that progressed to blastocysts and the group of oocytes that failed to develop. DNA methylation levels of corresponding regions and differentially methylated regions were calculated using customized Perl and R scripts. RNA-seq data were downloaded from a previously published paper and reanalysed.

MAIN RESULTS AND THE ROLE OF CHANCE

The results from PB1-MII oocyte pair validated that PB1 contains nearly the same methylome (average Pearson correlation is 0.92) with sibling MII oocyte. LINE and LTR expression increased markedly after fertilization. Moreover, the DNA methylation levels in LINE (including LINE1 and LINE2) and LTR were significantly higher in the PB1s of embryos that could not reach the blastocyst stage (Wilcoxon-Matt-Whitney test, P < 0.05). DNA methylation in PB1 promoters correlated negatively with gene expression of MII oocyte. Regarding the methylation status of the promoter regions, 66 genes were hypermethylated in the developmental arrested group, with their related functions (significantly enriched in several Gene Ontology terms) including transcription, positive regulation of adenylate cyclase activity, mitogen-activated protein kinase (MAPK) cascade and intracellular oestrogen receptor signalling pathway.

LARGE SCALE DATA

N/A.

LIMITATIONS, REASONS FOR CAUTION

Data analysis performed in this study focused on the competence of human oocytes and compared them with maternal genetic and epigenetic profiles. Therefore, data regarding the potential regulatory roles of paternal genomes in embryo development are lacking.

WIDER IMPLICATIONS OF THE FINDINGS

The results from PB1-oocyte pairs demonstrated that PB1s shared similar methylomes with their sibling oocytes. The selection of the good embryos for transfer should not only rely on morphology but also consider the DNA methylation of the corresponding PB1 and therefore MII oocyte. The application of early-stage analysis of PB1 offers an option for high-quality oocyte and embryo selection, which provides an additional tool for elective single embryo transfer in assisted reproduction.

STUDY FUNDING/COMPETING INTEREST(S)

This study was supported by the National Key Research and Development Program of China (2018YFC1004003, 2017YFA0103801), the National Natural Science Foundation of China (81730038, 3187144, 81521002) and the Strategic Priority Research Program of the Chinese Academy of Sciences (XDA16020703). The authors have no conflicts of interest to declare.

Transposable element insertion as a mechanism of SMARCB1 inactivation in atypical teratoid/rhabdoid tumor

Atypical teratoid/rhabdoid tumor (AT/RT) is a malignant brain tumor predominantly occurring in infants. Biallelic SMARCB1 mutations causing loss of nuclear SMARCB1/INI1 protein expression represent the characteristic genetic lesion. Pathogenic SMARCB1 mutations comprise single nucleotide variants, small insertions/deletions, large deletions, which may be also present in the germline (rhabdoid tumor predisposition syndrome 1), as well as somatic copy-number neutral loss of heterozygosity (LOH). In some SMARCB1-deficient AT/RT underlying biallelic mutations cannot be identified. Here we report the case of a 24-months-old girl diagnosed with a large brain tumor. The malignant rhabdoid tumor showed loss of nuclear SMARCB1/INI1 protein expression and the diagnosis of AT/RT was confirmed by DNA methylation profiling. While FISH, MLPA, Sanger sequencing and DNA methylation data-based imbalance analysis did not disclose alterations affecting SMARCB1, OncoScan array analysis revealed a 28.29 Mb sized region of copy-number neutral LOH on chromosome 22q involving the SMARCB1 locus. Targeted next-generation sequencing did also not detect a single nucleotide variant but instead revealed insertion of an AluY element into exon 2 of SMARCB1. Specific PCR-based Sanger sequencing verified the Alu insertion (SMARCB1 c.199_200 Alu ins) resulting in a frame-shift truncation not present in the patient's germline. In conclusion, transposable element insertion represents a hitherto not widely recognized mechanism of SMARCB1 disruption in AT/RT, which might not be detected by several widely applied conventional diagnostics assays. This finding has particular clinical implications, if rhabdoid predisposition syndrome 1 is suspected, but germline SMARCB1 alterations cannot be identified.

Haplotype-resolved diverse human genomes and integrated analysis of structural variation

Long-read and strand-specific sequencing technologies together facilitate the de novo assembly of high-quality haplotype-resolved human genomes without parent-child trio data. We present 64 assembled haplotypes from 32 diverse human genomes. These highly contiguous haplotype assemblies (average minimum contig length needed to cover 50% of the genome: 26 million base pairs) integrate all forms of genetic variation, even across complex loci. We identified 107,590 structural variants (SVs), of which 68% were not discovered with short-read sequencing, and 278 SV hotspots (spanning megabases of gene-rich sequence). We characterized 130 of the most active mobile element source elements and found that 63% of all SVs arise through homology-mediated mechanisms. This resource enables reliable graph-based genotyping from short reads of up to 50,340 SVs, resulting in the identification of 1526 expression quantitative trait loci as well as SV candidates for adaptive selection within the human population.

Phase separation of the LINE-1 ORF1 protein is mediated by the N-terminus and coiled-coil domain

Long interspersed nuclear element-1 (L1) is a retrotransposable element that autonomously replicates in the human genome, resulting in DNA damage and genomic instability. Activation of L1 in senescent cells triggers a type I interferon response and age-associated inflammation. Two open reading frames encode an ORF1 protein functioning as messenger RNA chaperone and an ORF2 protein providing catalytic activities necessary for retrotransposition. No function has been identified for the conserved, disordered N-terminal region of ORF1. Using microscopy and NMR spectroscopy, we demonstrate that ORF1 forms liquid droplets in vitro in a salt-dependent manner and that interactions between its N-terminal region and coiled-coil domain are necessary for phase separation. Mutations disrupting blocks of charged residues within the N-terminus impair phase separation, whereas some mutations within the coiled-coil domain enhance phase separation. Demixing of the L1 particle from the cytosol may provide a mechanism to protect the L1 transcript from degradation.

Viruses and Endogenous Retroviruses as Roots for Neuroinflammation and Neurodegenerative Diseases

Many neurodegenerative diseases are associated with chronic inflammation in the brain and periphery giving rise to a continuous imbalance of immune processes. Next to inflammation markers, activation of transposable elements, including long intrespersed nuclear elements (LINE) elements and endogenous retroviruses (ERVs), has been identified during neurodegenerative disease progression and even correlated with the clinical severity of the disease. ERVs are remnants of viral infections in the human genome acquired during evolution. Upon activation, they produce transcripts and the phylogenetically youngest ones are still able to produce viral-like particles. In addition, ERVs can bind transcription factors and modulate immune response. Being between own and foreign, ERVs are reviewed in the context of viral infections of the central nervous system, in aging and neurodegenerative diseases. Moreover, this review tests the hypothesis that viral infection may be a trigger at the onset of neuroinflammation and that ERVs sustain the inflammatory imbalance by summarizing existing data of neurodegenerative diseases associated with viruses and/or ERVs.

Involvement of transposable elements in neurogenesis

The article is about the role of transposons in the regulation of functioning of neuronal stem cells and mature neurons of the human brain. Starting from the first division of the zygote, embryonic development is governed by regular activations of transposable elements, which are necessary for the sequential regulation of the expression of genes specific for each cell type. These processes include differentiation of neuronal stem cells, which requires the finest tuning of expression of neuron genes in various regions of the brain. Therefore, in the hippocampus, the center of human neurogenesis, the highest transposon activity has been identified, which causes somatic mosaicism of cells during the formation of specific brain structures. Similar data were obtained in studies on experimental animals. Mobile genetic elements are the most important sources of long non-coding RNAs that are coexpressed with important brain protein-coding genes. Significant activity of long non-coding RNA was detected in the hippocampus, which confirms the role of transposons in the regulation of brain function. MicroRNAs, many of which arise from transposon transcripts, also play an important role in regulating the differentiation of neuronal stem cells. Therefore, transposons, through their own processed transcripts, take an active part in the epigenetic regulation of differentiation of neurons. The global regulatory role of transposons in the human brain is due to the emergence of protein-coding genes in evolution by their exonization, duplication and domestication. These genes are involved in an epigenetic regulatory network with the participation of transposons, since they contain nucleotide sequences complementary to miRNA and long non-coding RNA formed from transposons. In the memory formation, the role of the exchange of virus-like mRNA with the help of the Arc protein of endogenous retroviruses HERV between neurons has been revealed. A possible mechanism for the implementation of this mechanism may be reverse transcription of mRNA and site-specific insertion into the genome with a regulatory effect on the genes involved in the memory.

HIV suppressors against LINE-1: one functions as two

Exogenous retroviruses are RNA viruses that require reverse transcription for their replication. Among these viruses, human immunodeficiency virus (HIV) is infectious to humans and causes the development of acquired immune deficiency syndrome (AIDS). There are also endogenous retroelements that require reverse transcription for their retrotransposition, among which the type 1 long interspersed element (LINE-1) is the only type of retroelement that can replicate autonomously. It was once believed that retroviruses like HIV and retroelements like LINE-1 share similarities in processes such as reverse transcription and integration. Accordingly, many HIV suppressors are also potent LINE-1 inhibitors. However, in many cases, one suppressor uses two or more distinct mechanisms to repress HIV and LINE-1. In this review, we discuss some of these suppressors, focusing on their alternative mechanisms opposing the replication of HIV and LINE-1. Based on the differences in HIV and LINE-1 activity, the subcellular localization of these suppressors, and the impact of LINE-1 retrotransposition on human cells, we propose possible reasons for the inhibition of HIV and LINE-1 through different pathways by these suppressors, with the hope of accelerating future studies in associated research fields.

Transposable Element Landscape in Drosophila Populations Selected for Longevity

Transposable elements (TEs) inflict numerous negative effects on health and fitness as they replicate by integrating into new regions of the host genome. Even though organisms employ powerful mechanisms to demobilize TEs, transposons gradually lose repression during aging. The rising TE activity causes genomic instability and was implicated in age-dependent neurodegenerative diseases, inflammation, and the determination of lifespan. It is therefore conceivable that long-lived individuals have improved TE silencing mechanisms resulting in reduced TE expression relative to their shorter-lived counterparts and fewer genomic insertions. Here, we test this hypothesis by performing the first genome-wide analysis of TE insertions and expression in populations of Drosophila melanogaster selected for longevity through late-life reproduction for 50–170 generations from four independent studies. Contrary to our expectation, TE families were generally more abundant in long-lived populations compared with nonselected controls. Although simulations showed that this was not expected under neutrality, we found little evidence for selection driving TE abundance differences. Additional RNA-seq analysis revealed a tendency for reducing TE expression in selected populations, which might be more important for lifespan than regulating genomic insertions. We further find limited evidence of parallel selection on genes related to TE regulation and transposition. However, telomeric TEs were genomically and transcriptionally more abundant in long-lived flies, suggesting improved telomere maintenance as a promising TE-mediated mechanism for prolonging lifespan. Our results provide a novel viewpoint indicating that reproduction at old age increases the opportunity of TEs to be passed on to the next generation with little impact on longevity.

Transposable Elements and Stress in Vertebrates: An Overview

Since their identification as genomic regulatory elements, Transposable Elements (TEs) were considered, at first, molecular parasites and later as an important source of genetic diversity and regulatory innovations. In vertebrates in particular, TEs have been recognized as playing an important role in major evolutionary transitions and biodiversity. Moreover, in the last decade, a significant number of papers has been published highlighting a correlation between TE activity and exposition to environmental stresses and dietary factors. In this review we present an overview of the impact of TEs in vertebrate genomes, report the silencing mechanisms adopted by host genomes to regulate TE activity, and finally we explore the effects of environmental and dietary factor exposures on TE activity in mammals, which is the most studied group among vertebrates. The studies here reported evidence that several factors can induce changes in the epigenetic status of TEs and silencing mechanisms leading to their activation with consequent effects on the host genome. The study of TE can represent a future challenge for research for developing effective markers able to detect precocious epigenetic changes and prevent human diseases.

Light and shadow on the mechanisms of integration site selection in yeast Ty retrotransposon families

Transposable elements are ubiquitous in genomes. Their successful expansion depends in part on their sites of integration in their host genome. In Saccharomyces cerevisiae, evolution has selected various strategies to target the five Ty LTR-retrotransposon families into gene-poor regions in a genome, where coding sequences occupy 70% of the DNA. The integration of Ty1/Ty2/Ty4 and Ty3 occurs upstream and at the transcription start site of the genes transcribed by RNA polymerase III, respectively. Ty5 has completely different integration site preferences, targeting heterochromatin regions. Here, we review the history that led to the identification of the cellular tethering factors that play a major role in anchoring Ty retrotransposons to their preferred sites. We also question the involvement of additional factors in the fine-tuning of the integration site selection, with several studies converging towards an importance of the structure and organization of the chromatin.

Invasive and Noninvasive Nonfunctioning Gonadotroph Pituitary Tumors Differ in DNA Methylation Level of LINE-1 Repetitive Elements

PURPOSE

Epigenetic dysregulation plays a role in pituitary tumor pathogenesis. Some differences in DNA methylation were observed between invasive and noninvasive nonfunctioning gonadotroph tumors. This study sought to determine the role of DNA methylation changes in repetitive LINE-1 elements in nonfunctioning gonadotroph pituitary tumors.

METHODS

We investigated LINE-1 methylation levels in 80 tumors and normal pituitary glands with bisulfite-pyrosequencing. Expression of two LINE-1 open reading frames (L1-ORF1 and L1-ORF2) was analyzed with qRT-PCR in tumor samples and mouse gonadotroph pituitary cells treated with DNA methyltransferase inhibitor. Immunohistochemical staining against L1-ORF1p was also performed in normal pituitary glands and tumors.

RESULTS

Hypomethylation of LINE-1 was observed in pituitary tumors. Tumors characterized by invasive growth revealed lower LINE-1 methylation level than noninvasive ones. LINE-1 methylation correlated with overall DNA methylation assessed with HM450K arrays and negatively correlated with L1-ORF1 and L1-ORF2 expression. Treatment of αT3-1 gonadotroph cells with 5-Azacytidine clearly increased the level of L1-ORF1 and L1-ORF2 mRNA; however, its effect on LβT2 cells was less pronounced. Immunoreactivity against L1-ORF1p was higher in tumors than normal tissue. No difference in L1-ORF1p expression was observed in invasive and noninvasive tumors.

CONCLUSION

Hypomethylation of LINE-1 is related to invasive growth and influences transcriptional activity of transposable elements.

Transposable elements in brain health and disease

Transposable elements (TEs) occupy a large fraction of the human genome but only a small proportion of these elements are still active today. Recent works have suggested that TEs are expressed and active in the brain, challenging the dogma that neuronal genomes are static and revealing that they are susceptible to somatic genomic alterations. These new findings have major implications for understanding the neuroplasticity of the brain, which could hypothetically have a role in behavior and cognition, and contribute to vulnerability to disease. As active TEs could induce genetic diversity and mutagenesis, their influences on human brain development and diseases are of great interest. In this review, we will focus on the active TEs in the human genome and discuss in detail their impacts on human brain development. Furthermore, the association between TEs and brain-related diseases is discussed.

A BBS1 SVA F retrotransposon insertion is a frequent cause of Bardet-Biedl syndrome

Bardet-Biedl syndrome (BBS) is a ciliopathy characterized by retinitis pigmentosa, obesity, polydactyly, cognitive impairment and renal failure. Pathogenic variants in 24 genes account for the molecular basis of >80% of cases. Toward saturated discovery of the mutational basis of the disorder, we carefully explored our cohorts and identified a hominid-specific SINE-R/VNTR/Alu type F (SVA-F) insertion in exon 13 of BBS1 in eight families. In six families, the repeat insertion was found in trans with c.1169 T > G, p.Met390Arg and in two families the insertion was found in addition to other recessive BBS loci. Whole genome sequencing, de novo assembly and SNP array analysis were performed to characterize the genomic event. This insertion is extremely rare in the general population (found in 8 alleles of 8 BBS cases but not in >10 800 control individuals from gnomAD-SV) and due to a founder effect. Its 2435 bp sequence contains hallmarks of LINE1 mediated retrotransposition. Functional studies with patient-derived cell lines confirmed that the BBS1 SVA-F is deleterious as evidenced by a significant depletion of both mRNA and protein levels. Such findings highlight the importance of dedicated bioinformatics pipelines to identify all types of variation.

How Retroviruses and Retrotransposons in Our Genome May Contribute to Autoimmunity in Rheumatological Conditions

More than 200 human disorders include various manifestations of autoimmunity. The molecular events that lead to these diseases are still incompletely understood and their causes remain largely unknown. Numerous potential triggers of autoimmunity have been proposed over the years, but very few of them have been conclusively confirmed or firmly refuted. Viruses have topped the lists of suspects for decades, and it seems that many viruses, including those of the Herpesviridae family, indeed can influence disease initiation and/or promote exacerbations by a number of mechanisms that include prolonged anti-viral immunity, immune subverting factors, and mechanisms, and perhaps “molecular mimicry”. However, no specific virus has yet been established as being truly causative. Here, we discuss a different, but perhaps mechanistically related possibility, namely that retrotransposons or retroviruses that infected us in the past and left a lasting copy of themselves in our genome still can provoke an escalating immune response that leads to autoimmune disease. Many of these loci still encode for retroviral proteins that have retained some, or all, of their original functions. Importantly, these endogenous proviruses cannot be eliminated by the immune system the way it can eliminate exogenous viruses. Hence, if not properly controlled, they may drive a frustrated and escalating chronic, or episodic, immune response to the point of a frank autoimmune disorder. Here, we discuss the evidence and the proposed mechanisms, and assess the therapeutic options that emerge from the current understanding of this field.

p53 directly represses human LINE1 transposons

p53 is a potent tumor suppressor and commonly mutated in human cancers. Recently, we demonstrated that p53 genes act to restrict retrotransposons in germline tissues of flies and fish but whether this activity is conserved in somatic human cells is not known. Here we show that p53 constitutively restrains human LINE1s by cooperatively engaging sites in the 5'UTR and stimulating local deposition of repressive histone marks at these transposons. Consistent with this, the elimination of p53 or the removal of corresponding binding sites in LINE1s, prompted these retroelements to become hyperactive. Concurrently, p53 loss instigated chromosomal rearrangements linked to LINE sequences and also provoked inflammatory programs that were dependent on reverse transcriptase produced from LINE1s. Taken together, our observations establish that p53 continuously operates at the LINE1 promoter to restrict autonomous copies of these mobile elements in human cells. Our results further suggest that constitutive restriction of these retroelements may help to explain tumor suppression encoded by p53, since erupting LINE1s produced acute oncogenic threats when p53 was absent.

OneStopRNAseq: A Web Application for Comprehensive and Efficient Analyses of RNA-Seq Data

Over the past decade, a large amount of RNA sequencing (RNA-seq) data were deposited in public repositories, and more are being produced at an unprecedented rate. However, there are few open source tools with point-and-click interfaces that are versatile and offer streamlined comprehensive analysis of RNA-seq datasets. To maximize the capitalization of these vast public resources and facilitate the analysis of RNA-seq data by biologists, we developed a web application called OneStopRNAseq for the one-stop analysis of RNA-seq data. OneStopRNAseq has user-friendly interfaces and offers workflows for common types of RNA-seq data analyses, such as comprehensive data-quality control, differential analysis of gene expression, exon usage, alternative splicing, transposable element expression, allele-specific gene expression quantification, and gene set enrichment analysis. Users only need to select the desired analyses and genome build, and provide a Gene Expression Omnibus (GEO) accession number or Dropbox links to sequence files, alignment files, gene-expression-count tables, or rank files with the corresponding metadata. Our pipeline facilitates the comprehensive and efficient analysis of private and public RNA-seq data.

RNA-cDNA hybrids mediate transposition via different mechanisms

Retrotransposons can represent half of eukaryotic genomes. Retrotransposon dysregulation destabilizes genomes and has been linked to various human diseases. Emerging regulators of retromobility include RNA–DNA hybrid-containing structures known as R-loops. Accumulation of these structures at the transposons of yeast 1 (Ty1) elements has been shown to increase Ty1 retromobility through an unknown mechanism. Here, via a targeted genetic screen, we identified the rnh1Δ rad27Δ yeast mutant, which lacked both the Ty1 inhibitor Rad27 and the RNA–DNA hybrid suppressor Rnh1. The mutant exhibited elevated levels of Ty1 cDNA-associated RNA–DNA hybrids that promoted Ty1 mobility. Moreover, in this rnh1Δ rad27Δ mutant, but not in the double RNase H mutant rnh1Δ rnh201Δ, RNA–DNA hybrids preferentially existed as duplex nucleic acid structures and increased Ty1 mobility in a Rad52-dependent manner. The data indicate that in cells lacking RNA–DNA hybrid and Ty1 repressors, elevated levels of RNA-cDNA hybrids, which are associated with duplex nucleic acid structures, boost Ty1 mobility via a Rad52-dependent mechanism. In contrast, in cells lacking RNA–DNA hybrid repressors alone, elevated levels of RNA-cDNA hybrids, which are associated with triplex nucleic acid structures, boost Ty1 mobility via a Rad52-independent process. We propose that duplex and triplex RNA–DNA hybrids promote transposon mobility via Rad52-dependent or -independent mechanisms.

TIP_finder: An HPC Software to Detect Transposable Element Insertion Polymorphisms in Large Genomic Datasets

Transposable elements (TEs) are non-static genomic units capable of moving indistinctly from one chromosomal location to another. Their insertion polymorphisms may cause beneficial mutations, such as the creation of new gene function, or deleterious in eukaryotes, e.g., different types of cancer in humans. A particular type of TE called LTR-retrotransposons comprises almost 8% of the human genome. Among LTR retrotransposons, human endogenous retroviruses (HERVs) bear structural and functional similarities to retroviruses. Several tools allow the detection of transposon insertion polymorphisms (TIPs) but fail to efficiently analyze large genomes or large datasets. Here, we developed a computational tool, named TIP_finder, able to detect mobile element insertions in very large genomes, through high-performance computing (HPC) and parallel programming, using the inference of discordant read pair analysis. TIP_finder inputs are (i) short pair reads such as those obtained by Illumina, (ii) a chromosome-level reference genome sequence, and (iii) a database of consensus TE sequences. The HPC strategy we propose adds scalability and provides a useful tool to analyze huge genomic datasets in a decent running time. TIP_finder accelerates the detection of transposon insertion polymorphisms (TIPs) by up to 55 times in breast cancer datasets and 46 times in cancer-free datasets compared to the fastest available algorithms. TIP_finder applies a validated strategy to find TIPs, accelerates the process through HPC, and addresses the issues of runtime for large-scale analyses in the post-genomic era.

Promoter methylation, transcription, and retrotransposition of LINE-1 in colorectal adenomas and adenocarcinomas

BACKGROUND

The methylation of the CpG islands of the LINE-1 promoter is a tight control mechanism on the function of mobile elements. However, simultaneous quantification of promoter methylation and transcription of LINE-1 has not been performed in progressive stages of colorectal cancer. In addition, the insertion of mobile elements in the genome of advanced adenoma stage, a precancerous stage before colorectal carcinoma has not been emphasized. In this study, we quantify promoter methylation and transcripts of LINE-1 in three stages of colorectal non-advanced adenoma, advanced adenoma, and adenocarcinoma. In addition, we analyze the insertion of LINE-1, Alu, and SVA elements in the genome of patient tumors with colorectal advanced adenomas.

METHODS

LINE-1 hypomethylation status was evaluated by absolute quantitative analysis of methylated alleles (AQAMA) assay. To quantify the level of transcripts for LINE-1, quantitative RT-PCR was performed. To find mobile element insertions, the advanced adenoma tissue samples were subjected to whole genome sequencing and MELT analysis.

RESULTS

We found that the LINE-1 promoter methylation in advanced adenoma and adenocarcinoma was significantly lower than that in non-advanced adenomas. Accordingly, the copy number of LINE-1 transcripts in advanced adenoma was significantly higher than that in non-advanced adenomas, and in adenocarcinomas was significantly higher than that in the advanced adenomas. Whole-genome sequencing analysis of colorectal advanced adenomas revealed that at this stage polymorphic insertions of LINE-1, Alu, and SVA comprise approximately 16%, 51%, and 74% of total insertions, respectively.

CONCLUSIONS

Our correlative analysis showing a decreased methylation of LINE-1 promoter accompanied by the higher level of LINE-1 transcription, and polymorphic genomic insertions in advanced adenoma, suggests that the early and advanced polyp stages may host very important pathogenic processes concluding to cancer.

Early-onset breast cancer in a woman with a germline mobile element insertion resulting in BRCA2 disruption: a case report

Mobile element insertions (MEIs) contribute to genomic diversity, but they can be responsible for human disease in some cases. Initial clinical testing (BRCA1, BRCA2 and PALB2) in a 40-year-old female with unilateral breast cancer did not detect any pathogenic variants. Subsequent reanalysis for MEIs detected a novel likely pathogenic insertion of the retrotransposon element (RE) c.7894_7895insSVA in BRCA2. This case highlights the importance of bioinformatic pipeline optimization for the detection of MEIs in genes associated with hereditary cancer, as early detection can significantly impact clinical management.

LINE-1 ORF1p does not determine substrate preference for human/orangutan SVA and gibbon LAVA

Background

Non-autonomous VNTR (Variable Number of Tandem Repeats) composite retrotransposons – SVA (SINE-R-VNTR-Alu) and LAVA (L1-Alu-VNTR-Alu) – are specific to hominoid primates. SVA expanded in great apes, LAVA in gibbon. Both SVA and LAVA have been shown to be mobilized by the autonomous LINE-1 (L1)-encoded protein machinery in a cell-based assay in trans. The efficiency of human SVA retrotransposition in vitro has, however, been considerably lower than would be expected based on recent pedigree-based in vivo estimates. The VNTR composite elements across hominoids – gibbon LAVA, orangutan SVA_A descendants and hominine SVA_D descendants – display characteristic structures of the 5′ Alu-like domain and the VNTR. Different partner L1 subfamilies are currently active in each of the lineages. The possibility that the lineage-specific types of VNTR composites evolved in response to evolutionary changes in their autonomous partners, particularly in the nucleic acid binding L1 ORF1-encoded protein, has not been addressed.

Results

Here I report the identification and functional characterization of a highly active human SVA element using an improved mneo retrotransposition reporter cassette. The modified cassette (mneoM) minimizes splicing between the VNTR of human SVAs and the neomycin phosphotransferase stop codon. SVA deletion analysis provides evidence that key elements determining its mobilization efficiency reside in the VNTR and 5′ hexameric repeats. Simultaneous removal of the 5′ hexameric repeats and part of the VNTR has an additive negative effect on mobilization rates. Taking advantage of the modified reporter cassette that facilitates robust cross-species comparison of SVA/LAVA retrotransposition, I show that the ORF1-encoded proteins of the L1 subfamilies currently active in gibbon, orangutan and human do not display substrate preference for gibbon LAVA versus orangutan SVA versus human SVA. Finally, I demonstrate that an orangutan-derived ORF1p supports only limited retrotransposition of SVA/LAVA in trans, despite being fully functional in L1 mobilization in cis.

Conclusions

Overall, the analysis confirms SVA as a highly active human retrotransposon and preferred substrate of the L1-encoded protein machinery. Based on the results obtained in human cells coevolution of L1 ORF1p and VNTR composites does not appear very likely. The changes in orangutan L1 ORF1p that markedly reduce its mobilization capacity in trans might explain the different SVA insertion rates in the orangutan and hominine lineages, respectively.

Human TRIM5α senses and restricts LINE-1 elements

Mobile genetic elements have significantly shaped our genomic landscape. LINE-1 retroelements are the only autonomously active elements left in the human genome. Since new insertions can have detrimental consequences, cells need to efficiently control LINE-1 retrotransposition. Here, we demonstrate that the intrinsic immune factor TRIM5α senses and restricts LINE-1 retroelements. Previously, rhesus TRIM5α has been shown to efficiently block HIV-1 replication, while human TRIM5α was found to be less active. Surprisingly, we found that both human and rhesus TRIM5α efficiently repress human LINE-1 retrotransposition. TRIM5α interacts with LINE-1 ribonucleoprotein complexes in the cytoplasm, which is essential for restriction. In line with its postulated role as pattern recognition receptor, we show that TRIM5α also induces innate immune signaling upon interaction with LINE-1 ribonucleoprotein complexes. The signaling events activate the transcription factors AP-1 and NF-κB, leading to the down-regulation of LINE-1 promoter activity. Together, our findings identify LINE-1 as important target of human TRIM5α, which restricts and senses LINE-1 via two distinct mechanisms. Our results corroborate TRIM5α as pattern recognition receptor and shed light on its previously undescribed activity against mobile genetic elements, such as LINE-1, to protect the integrity of our genome.

The dynamic methylome of islets in health and disease

Background

Epigenetic processes control timing and level of gene expression throughout life, during development, differentiation, and aging, and are the link to adapting gene expression profiles to environmental cues. To qualify for the definition of ‘epigenetic’, a change to a gene's activity must be inherited through at least one mitotic division. Epigenetic mechanisms link changes in the environment to adaptions of the genome that do not rely on changes in the DNA sequence. In the past two decades, multiple studies have aimed to identify epigenetic mechanisms, and to define their role in development, differentiation and disease.

Scope of review

In this review, we will focus on the current knowledge of the epigenetic control of pancreatic beta cell maturation and dysfunction and its relationship to the development of islet cell failure in diabetes. Most of the data currently available have been obtained in mice, but we will summarize studies of human data as well. We will focus here on DNA methylation, as this is the most stable epigenetic mark, and least impacted by the variables inherent in islet procurement, isolation, and culture.

Major conclusions

DNA methylation patterns of beta cell are dynamic during maturation and during the diabetic process. In both cases, the changes occur at cell specific regulatory regions such as enhancers, where the methylation profile is cell type specific. Frequently, the differentially methylated regulatory elements are associated with key function genes such as PDX1, NKX6-1 and TCF7L2. During maturation, enhancers tend to become demethylated in association with increased activation of beta cell function genes and increased functionality, as indicated by glucose stimulated insulin secretion. Likewise, the changes to the DNA methylome that are present in pancreatic islets from diabetic donors are enriched in regulatory regions as well.

Retrocopying expands the functional repertoire of APOBEC3 antiviral proteins in primates

Host-virus arms races are inherently asymmetric; viruses evolve much more rapidly than host genomes. Thus, there is high interest in discovering mechanisms by which host genomes keep pace with rapidly evolving viruses. One family of restriction factors, the APOBEC3 (A3) cytidine deaminases, has undergone positive selection and expansion via segmental gene duplication and recombination. Here, we show that new copies of A3 genes have also been created in primates by reverse transcriptase-encoding elements like LINE-1 or endogenous retroviruses via a process termed retrocopying. First, we discovered that all simian primate genomes retain the remnants of an ancient A3 retrocopy: A3I. Furthermore, we found that some New World monkeys encode up to ten additional APOBEC3G (A3G) retrocopies. Some of these A3G retrocopies are transcribed in a variety of tissues and able to restrict retroviruses. Our findings suggest that host genomes co-opt retroelement activity in the germline to create new host restriction factors as another means to keep pace with the rapid evolution of viruses. (163)