L1Base 2: more retrotransposition-active LINE-1s, more mammalian genomes

Subcellular location of L1 retrotransposon-encoded ORF1p, reverse transcription products, and DNA sensors in lupus granulocytes

BACKGROUND

Systemic lupus erythematosus (SLE) is a chronic autoimmune disease with an unpredictable course of recurrent exacerbations alternating with more stable disease. SLE is characterized by broad immune activation and autoantibodies against double-stranded DNA and numerous proteins that exist in cells as aggregates with nucleic acids, such as Ro60, MOV10, and the L1 retrotransposon-encoded ORF1p.

RESULTS

Here we report that these 3 proteins are co-expressed and co-localized in a subset of SLE granulocytes and are concentrated in cytosolic dots that also contain DNA: RNA heteroduplexes and the DNA sensor ZBP1, but not cGAS. The DNA: RNA heteroduplexes vanished from the neutrophils when they were treated with a selective inhibitor of the L1 reverse transcriptase. We also report that ORF1p granules escape neutrophils during the extrusion of neutrophil extracellular traps (NETs) and, to a lesser degree, from neutrophils dying by pyroptosis, but not apoptosis.

CONCLUSIONS

These results bring new insights into the composition of ORF1p granules in SLE neutrophils and may explain, in part, why proteins in these granules become targeted by autoantibodies in this disease.

The transcription factor PAX5 activates human LINE1 retrotransposons to induce cellular senescence

As a hallmark of senescent cells, the derepression of Long Interspersed Elements 1 (LINE1) transcription results in accumulated LINE1 cDNA, which triggers the secretion of the senescence-associated secretory phenotype (SASP) and paracrine senescence in a cGAS-STING pathway-dependent manner. However, transcription factors that govern senescence-associated LINE1 reactivation remain ill-defined. Here, we predict several transcription factors that bind to human LINE1 elements to regulate their transcription by analyzing the conserved binding motifs in the 5'-untranslated regions (UTR) of the commonly upregulated LINE1 elements in different types of senescent cells. Further analysis reveals that PAX5 directly binds to LINE1 5'-UTR and the binding is enhanced in senescent cells. The enrichment of PAX5 at the 5'-UTR promotes cellular senescence and SASP by activating LINE1. We also demonstrate that the longevity gene SIRT6 suppresses PAX5 transcription by directly binding to the PAX5 promoter, and overexpressing PAX5 abrogates the suppressive effect of SIRT6 on stress-dependent cellular senescence. Our work suggests that PAX5 could serve as a potential target for drug development aiming to suppress LINE1 activation and treat senescence-associated diseases.

Transcription-coupled DNA-protein crosslink repair by CSB and CRL4CSA-mediated degradation

DNA-protein crosslinks (DPCs) arise from enzymatic intermediates, metabolism or chemicals like chemotherapeutics. DPCs are highly cytotoxic as they impede DNA-based processes such as replication, which is counteracted through proteolysis-mediated DPC removal by spartan (SPRTN) or the proteasome. However, whether DPCs affect transcription and how transcription-blocking DPCs are repaired remains largely unknown. Here we show that DPCs severely impede RNA polymerase II-mediated transcription and are preferentially repaired in active genes by transcription-coupled DPC (TC-DPC) repair. TC-DPC repair is initiated by recruiting the transcription-coupled nucleotide excision repair (TC-NER) factors CSB and CSA to DPC-stalled RNA polymerase II. CSA and CSB are indispensable for TC-DPC repair; however, the downstream TC-NER factors UVSSA and XPA are not, a result indicative of a non-canonical TC-NER mechanism. TC-DPC repair functions independently of SPRTN but is mediated by the ubiquitin ligase CRL4CSA and the proteasome. Thus, DPCs in genes are preferentially repaired in a transcription-coupled manner to facilitate unperturbed transcription.

LINE-1 mRNA 3' end dynamics shape its biology and retrotransposition potential

LINE-1 (L1) retrotransposons are mobile genetic elements that create new genomic insertions by a copy-paste mechanism involving L1 RNA/RNP intermediates. L1 encodes two ORFs, of which L1-ORF2p nicks genomic DNA and reverse transcribes L1 mRNA using the nicked DNA as a primer which base-pairs with poly(A) tail of L1 mRNA. To better understand the importance of non-templated L1 3' ends' dynamics and the interplay between L1 3' and 5' ends, we investigated the effects of genomic knock-outs and temporal knock-downs of XRN1, DCP2, and other factors. We hypothesized that in the absence of XRN1, the major 5'→3' exoribonuclease, there would be more L1 mRNA and retrotransposition. Conversely, we observed that loss of XRN1 decreased L1 retrotransposition. This occurred despite slight stabilization of L1 mRNA, but with decreased L1 RNP formation. Similarly, loss of DCP2, the catalytic subunit of the decapping complex, lowered retrotransposition despite increased steady-state levels of L1 proteins. In both XRN1 and DCP2 depletions we observed shortening of L1 3' poly(A) tails and their increased uridylation by TUT4/7. We explain the observed reduction of L1 retrotransposition by the changed qualities of non-templated L1 mRNA 3' ends demonstrating the important role of L1 3' end dynamics in L1 biology.

Improved LINE-1 Detection through Pattern Matching by Increasing Probe Length

Long Interspersed Element-1 (LINE-1 or L1) is an autonomous transposable element that accounts for 17% of the human genome. Strong correlations between abnormal L1 expression and diseases, particularly cancer, have been documented by numerous studies. L1PD (LINE-1 Pattern Detection) had been previously created to detect L1s by using a fixed pre-determined set of 50-mer probes and a pattern-matching algorithm. L1PD uses a novel seed-and-pattern-match strategy as opposed to the well-known seed-and-extend strategy employed by other tools. This study discusses an improved version of L1PD that shows how increasing the size of the k-mer probes from 50 to 75 or to 100 yields better results, as evidenced by experiments showing higher precision and recall when compared to the 50-mers. The probe-generation process was updated and the corresponding software is now shared so that users may generate probes for other reference genomes (with certain limitations). Additionally, L1PD was applied to other non-human genomes, such as dogs, horses, and cows, to further validate the pattern-matching strategy. The improved version of L1PD proves to be an efficient and promising approach for L1 detection.

Navigating the brain and aging: exploring the impact of transposable elements from health to disease

Transposable elements (TEs) are mobile genetic elements that constitute on average 45% of mammalian genomes. Their presence and activity in genomes represent a major source of genetic variability. While this is an important driver of genome evolution, TEs can also have deleterious effects on their hosts. A growing number of studies have focused on the role of TEs in the brain, both in physiological and pathological contexts. In the brain, their activity is believed to be important for neuronal plasticity. In neurological and age-related disorders, aberrant activity of TEs may contribute to disease etiology, although this remains unclear. After providing a comprehensive overview of transposable elements and their interactions with the host, this review summarizes the current understanding of TE activity within the brain, during the aging process, and in the context of neurological and age-related conditions.

Locus-specific LINE-1 expression in clinical ovarian cancer specimens at the single-cell level

Long interspersed nuclear elements (LINE-1s/L1s) are a group of retrotransposons that can copy themselves within a genome. In humans, it is the most successful transposon in nucleotide content. L1 expression is generally mild in normal human tissues, but the activity has been shown to increase significantly in many cancers. Few studies have examined L1 expression at single-cell resolution, thus it is undetermined whether L1 reactivation occurs solely in malignant cells within tumors. One of the cancer types with frequent L1 activity is high-grade serous ovarian carcinoma (HGSOC). Here, we identified locus-specific L1 expression with 3' single-cell RNA sequencing in pre- and post-chemotherapy HGSOC sample pairs from 11 patients, and in fallopian tube samples from five healthy women. Although L1 expression quantification with the chosen technique was challenging due to the repetitive nature of the element, we found evidence of L1 expression primarily in cancer cells, but also in other cell types, e.g. cancer-associated fibroblasts. The expression levels were similar in samples taken before and after neoadjuvant chemotherapy, indicating that L1 transcriptional activity was unaffected by clinical platinum-taxane treatment. Furthermore, L1 activity was negatively associated with the expression of MYC target genes, a finding that supports earlier literature of MYC being an L1 suppressor.

Locus-level L1 DNA methylation profiling reveals the epigenetic and transcriptional interplay between L1s and their integration sites

Long interspersed element 1 (L1) retrotransposons are implicated in human disease and evolution. Their global activity is repressed by DNA methylation, but deciphering the regulation of individual copies has been challenging. Here, we combine short- and long-read sequencing to unveil L1 methylation heterogeneity across cell types, families, and individual loci and elucidate key principles involved. We find that the youngest primate L1 families are specifically hypomethylated in pluripotent stem cells and the placenta but not in most tumors. Locally, intronic L1 methylation is intimately associated with gene transcription. Conversely, the L1 methylation state can propagate to the proximal region up to 300 bp. This phenomenon is accompanied by the binding of specific transcription factors, which drive the expression of L1 and chimeric transcripts. Finally, L1 hypomethylation alone is typically insufficient to trigger L1 expression due to redundant silencing pathways. Our results illuminate the epigenetic and transcriptional interplay between retrotransposons and their host genome.

A highly contiguous genome assembly for the pocket mouse Perognathus longimembris longimembris

The little pocket mouse, Perognathus longimembris, and its nine congeners are small heteromyid rodents found in arid and seasonally arid regions of Western North America. The genus is characterized by behavioral and physiological adaptations to dry and often harsh environments, including nocturnality, seasonal torpor, food caching, enhanced osmoregulation, and a well-developed sense of hearing. Here we present a genome assembly of Perognathus longimembris longimembris generated from PacBio HiFi long read and Omni-C chromatin-proximity sequencing as part of the California Conservation Genomics Project. The assembly has a length of 2.35 Gb, contig N50 of 11.6 Mb, scaffold N50 of 73.2 Mb, and includes 93.8% of the BUSCO Glires genes. Interspersed repetitive elements constitute 41.2% of the genome. A comparison with the highly endangered Pacific pocket mouse, P. l. pacificus, reveals broad synteny. These new resources will enable studies of local adaptation, genetic diversity, and conservation of threatened taxa.

Long-run real-time PCR analysis of repetitive nuclear elements as a novel tool for DNA damage quantification in single cells: an approach validated on mouse oocytes and fibroblasts

Since DNA damage is of great importance in various biological processes, its rate is frequently assessed both in research studies and in medical diagnostics. The most precise methods of quantifying DNA damage are based on real-time PCR. However, in the conventional version, they require a large amount of genetic material and therefore their usefulness is limited to multicellular samples. Here, we present a novel approach to long-run real-time PCR-based DNA-damage quantification (L1-LORD-Q), which consists in amplification of long interspersed nuclear elements (L1) and allows for analysis of single-cell genomes. The L1-LORD-Q was compared with alternative methods of measuring DNA breaks (Bioanalyzer system, γ-H2AX foci staining), which confirmed its accuracy. Furthermore, it was demonstrated that the L1-LORD-Q is sensitive enough to distinguish between different levels of UV-induced DNA damage. The method was validated on mouse oocytes and fibroblasts, but the general idea is universal and can be applied to various types of cells and species.

Aging differentially alters the transcriptome and landscape of chromatin accessibility in the male and female mouse hippocampus

Aging-related memory impairment and pathological memory disorders such as Alzheimer's disease differ between males and females, and yet little is known about how aging-related changes in the transcriptome and chromatin environment differ between sexes in the hippocampus. To investigate this question, we compared the chromatin accessibility landscape and gene expression/alternative splicing pattern of young adult and aged mouse hippocampus in both males and females using ATAC-seq and RNA-seq. We detected significant aging-dependent changes in the expression of genes involved in immune response and synaptic function and aging-dependent changes in the alternative splicing of myelin sheath genes. We found significant sex-bias in the expression and alternative splicing of hundreds of genes, including aging-dependent female-biased expression of myelin sheath genes and aging-dependent male-biased expression of genes involved in synaptic function. Aging was associated with increased chromatin accessibility in both male and female hippocampus, especially in repetitive elements, and with an increase in LINE-1 transcription. We detected significant sex-bias in chromatin accessibility in both autosomes and the X chromosome, with male-biased accessibility enriched at promoters and CpG-rich regions. Sex differences in gene expression and chromatin accessibility were amplified with aging, findings that may shed light on sex differences in aging-related and pathological memory loss.

LINE-1 hypomethylation, increased retrotransposition and tumor-specific insertion in upper gastrointestinal cancer

The long interspersed nuclear element-1 (LINE-1) retrotransposons are a major family of mobile genetic elements, comprising approximately 17% of the human genome. The methylation state of LINE-1 is often used as an indicator of global DNA methylation levels and it regulates the retrotransposition and somatic insertion of the genetic element. We have previously reported the significant relationship between LINE-1 hypomethylation and poor prognosis in upper gastrointestinal (GI) cancers. However, the causal relationships between LINE-1 hypomethylation, retrotransposition, and tumor-specific insertion in upper GI cancers remain unknown. We used bisulfite-pyrosequencing and quantitative real-time PCR to verify LINE-1 methylation and copy number in tissue samples of 101 patients with esophageal and 103 patients with gastric cancer. Furthermore, we analyzed the LINE-1 retrotransposition profile with an originally developed L1Hs-seq. In tumor samples, LINE-1 methylation levels were significantly lower than non-tumor controls, while LINE-1 copy numbers were markedly increased. As such, there was a significant inverse correlation between the LINE-1 methylation level and copy number in tumor tissues, with lower LINE-1 methylation levels corresponding to higher LINE-1 copy numbers. Of particular importance is that somatic LINE-1 insertions were more numerous in tumor than normal tissues. Furthermore, we observed that LINE-1 was inserted evenly across all chromosomes, and most often within genomic regions associated with tumor-suppressive genes. LINE-1 hypomethylation in upper GI cancers is related to increased LINE-1 retrotransposition and tumor-specific insertion events, which may collectively contribute to the acquisition of aggressive tumor features through the inactivation of tumor-suppressive genes.

Asymmetric distribution of parental H3K9me3 in S phase silences L1 elements

In eukaryotes, repetitive DNA sequences are transcriptionally silenced through histone H3 lysine 9 trimethylation (H3K9me3). Loss of silencing of the repeat elements leads to genome instability and human diseases, including cancer and ageing1-3. Although the role of H3K9me3 in the establishment and maintenance of heterochromatin silencing has been extensively studied4-6, the pattern and mechanism that underlie the partitioning of parental H3K9me3 at replicating DNA strands are unknown. Here we report that H3K9me3 is preferentially transferred onto the leading strands of replication forks, which occurs predominantly at long interspersed nuclear element (LINE) retrotransposons (also known as LINE-1s or L1s) that are theoretically transcribed in the head-on direction with replication fork movement. Mechanistically, the human silencing hub (HUSH) complex interacts with the leading-strand DNA polymerase Pol ε and contributes to the asymmetric segregation of H3K9me3. Cells deficient in Pol ε subunits (POLE3 and POLE4) or the HUSH complex (MPP8 and TASOR) show compromised H3K9me3 asymmetry and increased LINE expression. Similar results were obtained in cells expressing a MPP8 mutant defective in H3K9me3 binding and in TASOR mutants with reduced interactions with Pol ε. These results reveal an unexpected mechanism whereby the HUSH complex functions with Pol ε to promote asymmetric H3K9me3 distribution at head-on LINEs to suppress their expression in S phase.

Assessing the Expression of Long INterspersed Elements (LINEs) via Long-Read Sequencing in Diverse Human Tissues and Cell Lines

Transposable elements, such as Long INterspersed Elements (LINEs), are DNA sequences that can replicate within genomes. LINEs replicate using an RNA intermediate followed by reverse transcription and are typically a few kilobases in length. LINE activity creates genomic structural variants in human populations and leads to somatic alterations in cancer genomes. Long-read RNA sequencing technologies, including Oxford Nanopore and PacBio, can directly sequence relatively long transcripts, thus providing the opportunity to examine full-length LINE transcripts. This study focuses on the development of a new bioinformatics pipeline for the identification and quantification of active, full-length LINE transcripts in diverse human tissues and cell lines. In our pipeline, we utilized RepeatMasker to identify LINE-1 (L1) transcripts from long-read transcriptome data and incorporated several criteria, such as transcript start position, divergence, and length, to remove likely false positives. Comparisons between cancerous and normal cell lines, as well as human tissue samples, revealed elevated expression levels of young LINEs in cancer, particularly at intact L1 loci. By employing bioinformatics methodologies on long-read transcriptome data, this study demonstrates the landscape of L1 expression in tissues and cell lines.

Locus-resolution analysis of L1 regulation and retrotransposition potential in mouse embryonic development

Mice harbor ∼2800 intact copies of the retrotransposon Long Interspersed Element 1 (L1). The in vivo retrotransposition capacity of an L1 copy is defined by both its sequence integrity and epigenetic status, including DNA methylation of the monomeric units constituting young mouse L1 promoters. Locus-specific L1 methylation dynamics during development may therefore elucidate and explain spatiotemporal niches of endogenous retrotransposition but remain unresolved. Here, we interrogate the retrotransposition efficiency and epigenetic fate of source (donor) L1s, identified as mobile in vivo. We show that promoter monomer loss consistently attenuates the relative retrotransposition potential of their offspring (daughter) L1 insertions. We also observe that most donor/daughter L1 pairs are efficiently methylated upon differentiation in vivo and in vitro. We use Oxford Nanopore Technologies (ONT) long-read sequencing to resolve L1 methylation genome-wide and at individual L1 loci, revealing a distinctive "smile" pattern in methylation levels across the L1 promoter region. Using Pacific Biosciences (PacBio) SMRT sequencing of L1 5' RACE products, we then examine DNA methylation dynamics at the mouse L1 promoter in parallel with transcription start site (TSS) distribution at locus-specific resolution. Together, our results offer a novel perspective on the interplay between epigenetic repression, L1 evolution, and genome stability.

Reference LINE-1 insertion polymorphisms correlate with Parkinson's disease progression and differential transcript expression in the PPMI cohort

Long interspersed nuclear element-1 (LINE-1/L1) retrotransposons make up 17% of the human genome. They represent one class of transposable elements with the capacity to both mobilize autonomously and in trans via the mobilization of other elements, primarily Alu and SVA elements. Reference LINE-1 elements are, by definition, found in the reference genome, however, due to the polymorphic nature of these elements, variation for presence or absence is present within the population. We used a combination of clinical and transcriptomic data from the Parkinson's Progression Markers Initiative (PPMI) and applied matrix expression quantitative trait loci analysis and linear mixed-effects models involving 114 clinical, biochemical and imaging data from the PPMI cohort to elucidate the role of reference LINE-1 insertion polymorphism on both gene expression genome-wide and progression of Parkinson's disease (PD). We demonstrate that most LINE-1 insertion polymorphisms are capable of regulating gene expression, preferentially in trans, including previously identified PD risk loci. In addition, we show that 70 LINE-1 elements were associated with longitudinal changes of at least one PD progression marker, including ipsilateral count density ratio and UPDRS scores which are indicators of degeneration and severity. In conclusion, this study highlights the effect of the polymorphic nature of LINE-1 retrotransposons on gene regulation and progression of PD which underlines the importance of analyzing transposable elements within complex diseases.

Research progress of LINE-1 in the diagnosis, prognosis, and treatment of gynecologic tumors

The retrotransposon known as long interspersed nuclear element-1 (LINE-1), which is currently the sole autonomously mobile transposon in the human genome, can result in insertional mutations, chromosomal rearrangements, and genomic instability. In recent years, numerous studies have shown that LINE-1 is involved in the development of various diseases and also plays an important role in the immune regulation of the organism. The expression of LINE-1 in gynecologic tumors suggests that it is expected to be an independent indicator for early diagnosis and prognosis, and also, as a therapeutic target, LINE-1 is closely associated with gynecologic tumor prognosis. This article discusses the function of LINE-1 in the diagnosis, treatment, and prognosis of ovarian, cervical, and endometrial malignancies, as well as other gynecologic malignancies. It offers fresh perspectives on the early detection of tumors and the creation of novel anti-tumor medications.

Exploratory analysis of L1 retrotransposons expression in autism

BACKGROUND

Autism spectrum disorder (ASD) is a set of highly heterogeneous neurodevelopmental diseases whose genetic etiology is not completely understood. Several investigations have relied on transcriptome analysis from peripheral tissues to dissect ASD into homogenous molecular phenotypes. Recently, analysis of changes in gene expression from postmortem brain tissues has identified sets of genes that are involved in pathways previously associated with ASD etiology. In addition to protein-coding transcripts, the human transcriptome is composed by a large set of non-coding RNAs and transposable elements (TEs). Advancements in sequencing technologies have proven that TEs can be transcribed in a regulated fashion, and their dysregulation might have a role in brain diseases.

METHODS

We exploited published datasets comprising RNA-seq data from (1) postmortem brain of ASD subjects, (2) in vitro cell cultures where ten different ASD-relevant genes were knocked out and (3) blood of discordant siblings. We measured the expression levels of evolutionarily young full-length transposable L1 elements and characterized the genomic location of deregulated L1s assessing their potential impact on the transcription of ASD-relevant genes. We analyzed every sample independently, avoiding to pool together the disease subjects to unmask the heterogeneity of the molecular phenotypes.

RESULTS

We detected a strong upregulation of intronic full-length L1s in a subset of postmortem brain samples and in in vitro differentiated neurons from iPSC knocked out for ATRX. L1 upregulation correlated with an high number of deregulated genes and retained introns. In the anterior cingulate cortex of one subject, a small number of significantly upregulated L1s overlapped with ASD-relevant genes that were significantly downregulated, suggesting the possible existence of a negative effect of L1 transcription on host transcripts.

LIMITATIONS

Our analyses must be considered exploratory and will need to be validated in bigger cohorts. The main limitation is given by the small sample size and by the lack of replicates for postmortem brain samples. Measuring the transcription of locus-specific TEs is complicated by the repetitive nature of their sequence, which reduces the accuracy in mapping sequencing reads to the correct genomic locus.

CONCLUSIONS

L1 upregulation in ASD appears to be limited to a subset of subjects that are also characterized by a general deregulation of the expression of canonical genes and an increase in intron retention. In some samples from the anterior cingulate cortex, L1s upregulation seems to directly impair the expression of some ASD-relevant genes by a still unknown mechanism. L1s upregulation may therefore identify a group of ASD subjects with common molecular features and helps stratifying individuals for novel strategies of therapeutic intervention.

Genome-Wide Young L1 Methylation Profiling by bs-ATLAS-seq

By silencing L1 retrotransposons, DNA methylation protects mammalian genomes from potent endogenous mutagens. However, some loci can escape this repressive mechanism and become active, particularly in carcinomas. Alterations of L1 DNA methylation can also locally influence gene expression. Comprehensive measurement of L1 DNA methylation at the locus level remains challenging. Here, we present bs-ATLAS-seq, a genome-wide approach to locate full-length L1 elements in the human genome, and assess their methylation levels at single-base and single-locus resolutions. This strategy targets the youngest, and only retrotransposition-competent family, L1HS, but also detects a significant fraction of older elements (L1PA2 to L1PA8). Bs-ATLAS-seq evaluates methylation at the first 15 CpGs of L1 5' UTR, which corresponds to the first half of the sense promoter. It relies on random fragmentation of the genomic DNA, adapter ligation, bisulfite treatment and suppression PCR, and ends by asymmetrical paired-end sequencing. A dedicated pipeline provides the location of L1 elements and their methylation status, including for non-reference loci, as well as their single-molecule DNA profiles.

A Method for Detection of Somatic LINE-1 Insertions at the Single-Cell Level from Postmortem Human Brain

Long Interspersed Element-1 (LINE-1, L1) is a retrotransposon that has the ability to amplify its copy in the genome autonomously. L1Hs is a human-specific active subtype of L1 reported to amplify its copy in neural progenitor cells causing genomic mosaicism. This chapter describes a new method named NECO-seq (Novel Elements Concentrated-sequence) to identify the genomic locus of L1Hs insertions at the single-cell level. This protocol contains the steps of (1) preparation of neuronal cell nuclei from a postmortem human brain, (2) whole genome amplification from single neural nuclei (snWGA), (3) single nucleotide polymorphisms (SNPs) genotyping for quality control of snWGA products, (4) library preparation for next-generation sequencing to enrich the genomic locus of L1Hs insertions, and (5) bioinformatic analysis to detect novel somatic L1Hs insertions. This method can detect approximately 97% of L1Hs originally existing in reference human genome and approximately 10-20 newly inserted L1Hs copies in a neuronal cell of a postmortem human brain.

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.

Evolution of the ancestral mammalian karyotype and syntenic regions

Computational reconstruction of ancestral mammalian karyotypes revealed a comprehensive picture of the chromosome rearrangements that occurred over the evolutionary history of mammals. Ancient gene order, in some cases extending to full chromosomes, was found conserved for more than 300 My, demonstrating strong evolutionary constraint against rearrangements in some regions. Conserved segments of chromosomes are enriched for genes that control developmental processes. Therefore, Darwinian selection likely maintains ancient gene combinations while allowing for genomic innovations within or near chromosomal sites that break and rearrange over evolutionary time. The revealed relationship between the three-dimensional structure of chromosomes and the evolutionary stability of chromosome segments provides additional insights into the mechanisms of chromosome evolution and diseases associated with genome rearrangements.

Decrypting the rearrangements that drive mammalian chromosome evolution is critical to understanding the molecular bases of speciation, adaptation, and disease susceptibility. Using 8 scaffolded and 26 chromosome-scale genome assemblies representing 23/26 mammal orders, we computationally reconstructed ancestral karyotypes and syntenic relationships at 16 nodes along the mammalian phylogeny. Three different reference genomes (human, sloth, and cattle) representing phylogenetically distinct mammalian superorders were used to assess reference bias in the reconstructed ancestral karyotypes and to expand the number of clades with reconstructed genomes. The mammalian ancestor likely had 19 pairs of autosomes, with nine of the smallest chromosomes shared with the common ancestor of all amniotes (three still conserved in extant mammals), demonstrating a striking conservation of synteny for ∼320 My of vertebrate evolution. The numbers and types of chromosome rearrangements were classified for transitions between the ancestral mammalian karyotype, descendent ancestors, and extant species. For example, 94 inversions, 16 fissions, and 14 fusions that occurred over 53 My differentiated the therian from the descendent eutherian ancestor. The highest breakpoint rate was observed between the mammalian and therian ancestors (3.9 breakpoints/My). Reconstructed mammalian ancestor chromosomes were found to have distinct evolutionary histories reflected in their rates and types of rearrangements. The distributions of genes, repetitive elements, topologically associating domains, and actively transcribed regions in multispecies homologous synteny blocks and evolutionary breakpoint regions indicate that purifying selection acted over millions of years of vertebrate evolution to maintain syntenic relationships of developmentally important genes and regulatory landscapes of gene-dense chromosomes.

Pattern matching for high precision detection of LINE-1s in human genomes

Background

Long interspersed element 1 (LINE-1 or L1) retrotransposons are mobile elements that constitute 17–20% of the human genome. Strong correlations between abnormal L1 expression and several human diseases have been reported. This has motivated increasing interest in accurate quantification of the number of L1 copies present in any given biologic specimen. A main obstacle toward this aim is that L1s are relatively long DNA segments with regions of high variability, or largely present in the human genome as truncated fragments. These particularities render traditional alignment strategies, such as seed-and-extend inefficient, as the number of segments that are similar to L1s explodes exponentially. This study uses the pattern matching methodology for more accurate identification of L1s. We validate experimentally the superiority of pattern matching for L1 detection over alternative methods and discuss some of its potential applications.

Results

Pattern matching detected full-length L1 copies with high precision, reasonable computational time, and no prior input information. It also detected truncated and significantly altered copies of L1 with relatively high precision. The method was effectively used to annotate L1s in a target genome and to calculate copy number variation with respect to a reference genome. Crucial to the success of implementation was the selection of a small set of k-mer probes from a set of sequences presenting a stable pattern of distribution in the genome. As in seed-and-extend methods, the pattern matching algorithm sowed these k-mer probes, but instead of using heuristic extensions around the seeds, the analysis was based on distribution patterns within the genome. The desired level of precision could be adjusted, with some loss of recall.

Conclusion

Pattern matching is more efficient than seed-and-extend methods for the detection of L1 segments whose characterization depends on a finite set of sequences with common areas of low variability. We propose that pattern matching may help establish correlations between L1 copy number and disease states associated with L1 mobilization and evolution.

Recent, full-length gene retrocopies are common in canids

Gene retrocopies arise from the reverse transcription and insertion into the genome of processed mRNA transcripts. Although many retrocopies have acquired mutations that render them functionally inactive, most mammals retain active LINE-1 sequences capable of producing new retrocopies. New retrocopies, referred to as retro copy number variants (retroCNVs), may not be identified by standard variant calling techniques in high-throughput sequencing data. Although multiple functional FGF4 retroCNVs have been associated with skeletal dysplasias in dogs, the full landscape of canid retroCNVs has not been characterized. Here, retroCNV discovery was performed on a whole-genome sequencing data set of 293 canids from 76 breeds. We identified retroCNV parent genes via the presence of mRNA-specific 30-mers, and then identified retroCNV insertion sites through discordant read analysis. In total, we resolved insertion sites for 1911 retroCNVs from 1179 parent genes, 1236 of which appeared identical to their parent genes. Dogs had on average 54.1 total retroCNVs and 1.4 private retroCNVs. We found evidence of expression in testes for 12% (14/113) of the retroCNVs identified in six Golden Retrievers, including four chimeric transcripts, and 97 retroCNVs also had significantly elevated F ST across dog breeds, possibly indicating selection. We applied our approach to a subset of human genomes and detected an average of 4.2 retroCNVs per sample, highlighting a 13-fold relative increase of retroCNV frequency in dogs. Particularly in canids, retroCNVs are a largely unexplored source of genetic variation that can contribute to genome plasticity and that should be considered when investigating traits and diseases.

Melatonin: Regulation of Viral Phase Separation and Epitranscriptomics in Post-Acute Sequelae of COVID-19

The relentless, protracted evolution of the SARS-CoV-2 virus imposes tremendous pressure on herd immunity and demands versatile adaptations by the human host genome to counter transcriptomic and epitranscriptomic alterations associated with a wide range of short- and long-term manifestations during acute infection and post-acute recovery, respectively. To promote viral replication during active infection and viral persistence, the SARS-CoV-2 envelope protein regulates host cell microenvironment including pH and ion concentrations to maintain a high oxidative environment that supports template switching, causing extensive mitochondrial damage and activation of pro-inflammatory cytokine signaling cascades. Oxidative stress and mitochondrial distress induce dynamic changes to both the host and viral RNA m6A methylome, and can trigger the derepression of long interspersed nuclear element 1 (LINE1), resulting in global hypomethylation, epigenetic changes, and genomic instability. The timely application of melatonin during early infection enhances host innate antiviral immune responses by preventing the formation of "viral factories" by nucleocapsid liquid-liquid phase separation that effectively blockades viral genome transcription and packaging, the disassembly of stress granules, and the sequestration of DEAD-box RNA helicases, including DDX3X, vital to immune signaling. Melatonin prevents membrane depolarization and protects cristae morphology to suppress glycolysis via antioxidant-dependent and -independent mechanisms. By restraining the derepression of LINE1 via multifaceted strategies, and maintaining the balance in m6A RNA modifications, melatonin could be the quintessential ancient molecule that significantly influences the outcome of the constant struggle between virus and host to gain transcriptomic and epitranscriptomic dominance over the host genome during acute infection and PASC.

Dynamic reprogramming of H3K9me3 at hominoid-specific retrotransposons during human preimplantation development

Reprogramming of H3K9me3-dependent heterochromatin is required for early development. How H3K9me3 is involved in early human development remains, however, largely unclear. Here, we resolve the temporal landscape of H3K9me3 during human preimplantation development and its regulation for diverse hominoid-specific retrotransposons. At the 8-cell stage, H3K9me3 reprogramming at hominoid-specific retrotransposons termed SINE-VNTR-Alu (SVA) facilitates interaction between certain promoters and SVA-derived enhancers, promoting the zygotic genome activation. In trophectoderm, de novo H3K9me3 domains prevent pluripotent transcription factors from binding to hominoid-specific retrotransposons-derived regulatory elements for inner cell mass (ICM)-specific genes. H3K9me3 re-establishment at SVA elements in the ICM is associated with higher transcription of DNA repair genes, when compared with naive human pluripotent stem cells. Our data demonstrate that species-specific reorganization of H3K9me3-dependent heterochromatin at hominoid-specific retrotransposons plays important roles during early human development, shedding light on how the epigenetic regulation for early development has evolved in mammals.

ADAR1 masks the cancer immunotherapeutic promise of ZBP1-driven necroptosis

Only a small proportion of patients with cancer show lasting responses to immune checkpoint blockade (ICB)-based monotherapies. The RNA-editing enzyme ADAR1 is an emerging determinant of resistance to ICB therapy and prevents ICB responsiveness by repressing immunogenic double-stranded RNAs (dsRNAs), such as those arising from the dysregulated expression of endogenous retroviral elements (EREs)1-4. These dsRNAs trigger an interferon-dependent antitumour response by activating A-form dsRNA (A-RNA)-sensing proteins such as MDA-5 and PKR5. Here we show that ADAR1 also prevents the accrual of endogenous Z-form dsRNA elements (Z-RNAs), which were enriched in the 3' untranslated regions of interferon-stimulated mRNAs. Depletion or mutation of ADAR1 resulted in Z-RNA accumulation and activation of the Z-RNA sensor ZBP1, which culminated in RIPK3-mediated necroptosis. As no clinically viable ADAR1 inhibitors currently exist, we searched for a compound that can override the requirement for ADAR1 inhibition and directly activate ZBP1. We identified a small molecule, the curaxin CBL0137, which potently activates ZBP1 by triggering Z-DNA formation in cells. CBL0137 induced ZBP1-dependent necroptosis in cancer-associated fibroblasts and reversed ICB unresponsiveness in mouse models of melanoma. Collectively, these results demonstrate that ADAR1 represses endogenous Z-RNAs and identifies ZBP1-mediated necroptosis as a new determinant of tumour immunogenicity masked by ADAR1. Therapeutic activation of ZBP1-induced necroptosis provides a readily translatable avenue for rekindling the immune responsiveness of ICB-resistant human cancers.

SARS-CoV-2 and human retroelements: a case for molecular mimicry?

Background

The factors driving the late phase of COVID-19 are still poorly understood. However, autoimmunity is an evolving theme in COVID-19’s pathogenesis. Additionally, deregulation of human retroelements (RE) is found in many viral infections, and has also been reported in COVID-19.

Results

Unexpectedly, coronaviruses (CoV) – including SARS-CoV-2 – harbour many RE-identical sequences (up to 35 base pairs), and some of these sequences are part of SARS-CoV-2 epitopes associated to COVID-19 severity. Furthermore, RE are expressed in healthy controls and human cells and become deregulated after SARS-CoV-2 infection, showing mainly changes in long interspersed nuclear element (LINE1) expression, but also in endogenous retroviruses.

Conclusion

CoV and human RE share coding sequences, which are targeted by antibodies in COVID-19 and thus could induce an autoimmune loop by molecular mimicry.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-022-01040-2.

Locus specific reduction of L1 expression in the cortices of individuals with amyotrophic lateral sclerosis

The activation and dysregulation of retrotransposons has been identified in the CNS of individuals with the fatal neurodegenerative disorder Amyotrophic lateral sclerosis (ALS). This includes elements from multiple different families and subfamilies of retrotransposons, however there is limited knowledge of the specific loci from which this expression occurs in ALS. The long interspersed element-1 (L1) is the only autonomous retrotransposon in the human genome and members of this family of elements maintain the ability to mobilise. Despite L1s contributing to 17% of the human genome only 80-100 L1s encode the required proteins for mobilisation and are retrotransposition competent. Identifying the specific loci from which L1 expression occurs will inform on the potential functional consequences of their expression, such as the potential for somatic retrotransposition or DNA damage caused by the endonuclease activity of the ORF2 protein of the L1. Here we characterised L1 loci expression using the L1EM tool ( https://github.com/FenyoLab/L1EM ) in RNA sequencing data from 518 samples across four tissues (motor cortex, frontal cortex, cerebellum and cervical spinal cord) in the Target ALS cohort obtained from the New York Genome Center. There was a significant reduction in total intact L1 expression (those that encode functional proteins) in two brain regions of individuals with ALS compared to controls and clustering of the ALS brain regions occurred based on their intact L1 expression profile. Although overall the levels of L1 expression were reduced in ALS/ALS with other neurological disorder (ND) there were individuals in which L1s were expressed at much higher levels than the rest of the ALS/ALSND cohort. Expressed L1 loci were more frequently located in introns compared to those not expressed and the level of L1 expression positively correlated with the expression of the gene in which it was located. Significant differences were observed in the expression profiles of L1s in ALS and specific features of these elements, such as location in the genome and whether or not they are intact, were significantly associated with those that were expressed in the cohort.

ZNF146/OZF and ZNF507 target LINE-1 sequences

Repetitive sequences including transposable elements and transposon-derived fragments account for nearly half of the human genome. While transposition-competent transposable elements must be repressed to maintain genomic stability, mutated and fragmented transposable elements comprising the bulk of repetitive sequences can also contribute to regulation of host gene expression and broader genome organization. Here, we analyzed published ChIP-seq data sets to identify proteins broadly enriched on transposable elements in the human genome. We show 2 of the proteins identified, C2H2 zinc finger-containing proteins ZNF146 (also known as OZF) and ZNF507, are targeted to distinct sites within LINE-1 ORF2 at thousands of locations in the genome. ZNF146 binding sites are found at old and young LINE-1 elements. In contrast, ZNF507 preferentially binds at young LINE-1 sequences correlated to sequence changes in LINE-1 elements at ZNF507's binding site. To gain further insight into ZNF146 and ZNF507 function, we disrupt their expression in HEK293 cells using CRISPR/Cas9 and perform RNA sequencing, finding modest gene expression changes in cells where ZNF507 has been disrupted. We further identify a physical interaction between ZNF507 and PRMT5, suggesting ZNF507 may target arginine methylation activity to LINE-1 sequences.

The Absence of Retroelement Activity Is Characteristic for Childhood Acute Leukemias and Adult Acute Lymphoblastic Leukemia

Retroelements (RE) have been proposed as important players in cancerogenesis. Different cancer types are characterized by a different level of tumor-specific RE insertions. In previous studies, small cohorts of hematological malignancies, such as acute myeloid leukemia, multiple myeloma, and chronic lymphocytic leukemia have been characterized by a low level of RE insertional activity. Acute lymphoblastic leukemia (ALL) in adults and childhood acute leukemias have not been studied in this context. We performed a search for new RE insertions (Alu and L1) in 44 childhood ALL, 14 childhood acute myeloid leukemia, and 14 adult ALL samples using a highly sensitive NGS-based approach. First, we evaluated the method sensitivity revealing the 1% detection threshold for the proportion of cells with specific RE insertion. Following this result, we did not identify new tumor-specific RE insertions in the tested cohort of acute leukemia samples at the established level of sensitivity. Additionally, we analyzed the transcription levels of active L1 copies and found them increased. Thus, the increased transcription of active L1 copies is not sufficient for overt elevation of L1 retrotranspositional activity in leukemia.

Analysis of epigenetic features characteristic of L1 loci expressed in human cells

Only a select few L1 loci in the human genome are expressed in any given cell line or organ, likely to minimize damage done to the genome. The epigenetic features and requirements of expressed L1 loci are currently unknown. Using human cells and comprehensive epigenetic analysis of individual expressed and unexpressed L1 loci, we determined that endogenous L1 transcription depends on a combination of epigenetic factors, including open chromatin, activating histone modifications, and hypomethylation at the L1 promoter. We demonstrate that the L1 promoter seems to require interaction with enhancer elements for optimal function. We utilize epigenetic context to predict the expression status of L1Hs loci that are poorly mappable with RNA-Seq. Our analysis identified a population of ‘transitional’ L1 loci that likely have greater potential to be activated during the epigenetic dysregulation seen in tumors and during aging because they are the most responsive to targeted CRISPR-mediated delivery of trans-activating domains. We demonstrate that an engineered increase in endogenous L1 mRNA expression increases Alu mobilization. Overall, our findings present the first global and comprehensive analysis of epigenetic status of individual L1 loci based on their expression status and demonstrate the importance of epigenetic context for L1 expression heterogeneity.

Retrotransposons as a Source of DNA Damage in Neurodegeneration

The etiology of aging-associated neurodegenerative diseases (NDs), such as Parkinson’s disease (PD) and Alzheimer’s disease (AD), still remains elusive and no curative treatment is available. Age is the major risk factor for PD and AD, but the molecular link between aging and neurodegeneration is not fully understood. Aging is defined by several hallmarks, some of which partially overlap with pathways implicated in NDs. Recent evidence suggests that aging-associated epigenetic alterations can lead to the derepression of the LINE-1 (Long Interspersed Element-1) family of transposable elements (TEs) and that this derepression might have important implications in the pathogenesis of NDs. Almost half of the human DNA is composed of repetitive sequences derived from TEs and TE mobility participated in shaping the mammalian genomes during evolution. Although most TEs are mutated and no longer mobile, more than 100 LINE-1 elements have retained their full coding potential in humans and are thus retrotransposition competent. Uncontrolled activation of TEs has now been reported in various models of neurodegeneration and in diseased human brain tissues. We will discuss in this review the potential contribution of LINE-1 elements in inducing DNA damage and genomic instability, which are emerging pathological features in NDs. TEs might represent an important molecular link between aging and neurodegeneration, and a potential target for urgently needed novel therapeutic disease-modifying interventions.

Genome surveillance by HUSH-mediated silencing of intronless mobile elements

All life forms defend their genome against DNA invasion. Eukaryotic cells recognize incoming DNA and limit its transcription through repressive chromatin modifications. The human silencing hub (HUSH) complex transcriptionally represses long interspersed element-1 retrotransposons (L1s) and retroviruses through histone H3 lysine 9 trimethylation (H3K9me3)1-3. How HUSH recognizes and initiates silencing of these invading genetic elements is unknown. Here we show that HUSH is able to recognize and transcriptionally repress a broad range of long, intronless transgenes. Intron insertion into HUSH-repressed transgenes counteracts repression, even in the absence of intron splicing. HUSH binds transcripts from the target locus, prior to and independent of H3K9me3 deposition, and target transcription is essential for both initiation and propagation of HUSH-mediated H3K9me3. Genomic data reveal how HUSH binds and represses a subset of endogenous intronless genes generated through retrotransposition of cellular mRNAs. Thus intronless cDNA-the hallmark of reverse transcription-provides a versatile way to distinguish invading retroelements from host genes and enables HUSH to protect the genome from 'non-self' DNA, despite there being no previous exposure to the invading element. Our findings reveal the existence of a transcription-dependent genome-surveillance system and explain how it provides immediate protection against newly acquired elements while avoiding inappropriate repression of host genes.

Variant calling: Considerations, practices, and developments

The success of many clinical, association, or population genetics studies critically relies on properly performed variant calling step. The variety of modern genomics protocols, techniques, and platforms makes our choices of methods and algorithms difficult and there is no "one size fits all" solution for study design and data analysis. In this review, we discuss considerations that need to be taken into account while designing the study and preparing for the experiments. We outline the variety of variant types that can be detected using sequencing approaches and highlight some specific requirements and basic principles of their detection. Finally, we cover interesting developments that enable variant calling for a broad range of applications in the genomics field. We conclude by discussing technological and algorithmic advances that have the potential to change the ways of calling DNA variants in the nearest future.

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.

Factors Regulating the Activity of LINE1 Retrotransposons

LINE-1 (L1) is a class of autonomous mobile genetic elements that form somatic mosaicisms in various tissues of the organism. The activity of L1 retrotransposons is strictly controlled by many factors in somatic and germ cells at all stages of ontogenesis. Alteration of L1 activity was noted in a number of diseases: in neuropsychiatric and autoimmune diseases, as well as in various forms of cancer. Altered activity of L1 retrotransposons for some pathologies is associated with epigenetic changes and defects in the genes involved in their repression. This review discusses the molecular genetic mechanisms of the retrotransposition and regulation of the activity of L1 elements. The contribution of various factors controlling the expression and distribution of L1 elements in the genome occurs at all stages of the retrotransposition. The regulation of L1 elements at the transcriptional, post-transcriptional and integration into the genome stages is described in detail. Finally, this review also focuses on the evolutionary aspects of L1 accumulation and their interplay with the host regulation system.

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.

A Model-Driven Quantitative Analysis of Retrotransposon Distributions in the Human Genome

Retrotransposons, DNA sequences capable of creating copies of themselves, compose about half of the human genome and played a central role in the evolution of mammals. Their current position in the host genome is the result of the retrotranscription process and of the following host genome evolution. We apply a model from statistical physics to show that the genomic distribution of the two most populated classes of retrotransposons in human deviates from random placement, and that this deviation increases with time. The time dependence suggests a major role of the host genome dynamics in shaping the current retrotransposon distributions. Focusing on a neutral scenario, we show that a simple model based on random placement followed by genome expansion and sequence duplications can reproduce the empirical retrotransposon distributions, even though more complex and possibly selective mechanisms can have contributed. Besides the inherent interest in understanding the origin of current retrotransposon distributions, this work sets a general analytical framework to analyze quantitatively the effects of genome evolutionary dynamics on the distribution of genomic elements.

Formation of spermatogonia and fertile oocytes in golden hamsters requires piRNAs

PIWI-interacting RNAs (piRNAs) support the germline by suppressing retrotransposons. Studies of the pathway in mice have strongly shaped the view that mammalian piRNAs are essential for male but not for female fertility. Here, we report that the role of the piRNA pathway substantially differs in golden hamsters (Mesocricetus auratus), the piRNA pathway setup of which more closely resembles that of other mammals, including humans. The loss of the Mov10l1 RNA helicase—an essential piRNA biogenesis factor—leads to striking phenotypes in both sexes. In contrast to mice, female Mov10l1^–/– hamsters are sterile because their oocytes do not sustain zygotic development. Furthermore, Mov10l1^–/– male hamsters have impaired establishment of spermatogonia accompanied by transcriptome dysregulation and an expression surge of a young retrotransposon subfamily. Our results show that the mammalian piRNA pathway has essential roles in both sexes and its adaptive nature allows it to manage emerging genomic threats and acquire new critical roles in the germline.

A set of three papers reports that the piRNA pathway is essential for mammalian female fertility based on genetic perturbation experiments performed in golden hamsters.

DeepGRP: engineering a software tool for predicting genomic repetitive elements using Recurrent Neural Networks with attention

BACKGROUND

Repetitive elements contribute a large part of eukaryotic genomes. For example, about 40 to 50% of human, mouse and rat genomes are repetitive. So identifying and classifying repeats is an important step in genome annotation. This annotation step is traditionally performed using alignment based methods, either in a de novo approach or by aligning the genome sequence to a species specific set of repetitive sequences. Recently, Li (Bioinformatics 35:4408-4410, 2019) developed a novel software tool dna-brnn to annotate repetitive sequences using a recurrent neural network trained on sample annotations of repetitive elements.

RESULTS

We have developed the methods of dna-brnn further and engineered a new software tool DeepGRP. This combines the basic concepts of Li (Bioinformatics 35:4408-4410, 2019) with current techniques developed for neural machine translation, the attention mechanism, for the task of nucleotide-level annotation of repetitive elements. An evaluation on the human genome shows a 20% improvement of the Matthews correlation coefficient for the predictions delivered by DeepGRP, when compared to dna-brnn. DeepGRP predicts two additional classes of repeats (compared to dna-brnn) and is able to transfer repeat annotations, using RepeatMasker-based training data to a different species (mouse). Additionally, we could show that DeepGRP predicts repeats annotated in the Dfam database, but not annotated by RepeatMasker. DeepGRP is highly scalable due to its implementation in the TensorFlow framework. For example, the GPU-accelerated version of DeepGRP is approx. 1.8 times faster than dna-brnn, approx. 8.6 times faster than RepeatMasker and over 100 times faster than HMMER searching for models of the Dfam database.

CONCLUSIONS

By incorporating methods from neural machine translation, DeepGRP achieves a consistent improvement of the quality of the predictions compared to dna-brnn. Improved running times are obtained by employing TensorFlow as implementation framework and the use of GPUs. By incorporating two additional classes of repeats, DeepGRP provides more complete annotations, which were evaluated against three state-of-the-art tools for repeat annotation.

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.

Perturbed DNA methylation by Gadd45b induces chromatin disorganization, DNA strand breaks and dopaminergic neuron death

Age is a major risk factor for neurodegenerative diseases like Parkinson's disease, but few studies have explored the contribution of key hallmarks of aging, namely DNA methylation changes and heterochromatin destructuration, in the neurodegenerative process. Here, we investigated the consequences of viral overexpression of Gadd45b, a multifactorial protein involved in DNA demethylation, in the mouse midbrain. Gadd45b overexpression induced global and stable changes in DNA methylation, particularly in introns of genes related to neuronal functions, as well as on LINE-1 transposable elements. This was paralleled by disorganized heterochromatin, increased DNA damage, and vulnerability to oxidative stress. LINE-1 de-repression, a potential source of DNA damage, preceded Gadd45b-induced neurodegeneration, whereas prolonged Gadd45b expression deregulated expression of genes related to heterochromatin maintenance, DNA methylation, or Parkinson's disease. Our data indicates that aging-related alterations contribute to dopaminergic neuron degeneration with potential implications for Parkinson's disease.

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.

Endogenous retroviruses promote homeostatic and inflammatory responses to the microbiota

The microbiota plays a fundamental role in regulating host immunity. However, the processes involved in the initiation and regulation of immunity to the microbiota remain largely unknown. Here, we show that the skin microbiota promotes the discrete expression of defined endogenous retroviruses (ERVs). Keratinocyte-intrinsic responses to ERVs depended on cyclic GMP-AMP synthase (cGAS)/stimulator of interferon genes protein (STING) signaling and promoted the induction of commensal-specific T cells. Inhibition of ERV reverse transcription significantly impacted these responses, resulting in impaired immunity to the microbiota and its associated tissue repair function. Conversely, a lipid-enriched diet primed the skin for heightened ERV- expression in response to commensal colonization, leading to increased immune responses and tissue inflammation. Together, our results support the idea that the host may have co-opted its endogenous virome as a means to communicate with the exogenous microbiota, resulting in a multi-kingdom dialog that controls both tissue homeostasis and inflammation.

Organ-, sex- and age-dependent patterns of endogenous L1 mRNA expression at a single locus resolution

Expression of L1 mRNA, the first step in the L1 copy-and-paste amplification cycle, is a prerequisite for L1-associated genomic instability. We used a reported stringent bioinformatics method to parse L1 mRNA transcripts and measure the level of L1 mRNA expressed in mouse and rat organs at a locus-specific resolution. This analysis determined that mRNA expression of L1 loci in rodents exhibits striking organ specificity with less than 0.8% of loci shared between organs of the same organism. This organ specificity in L1 mRNA expression is preserved in male and female mice and across age groups. We discovered notable differences in L1 mRNA expression between sexes with only 5% of expressed L1 loci shared between male and female mice. Moreover, we report that the levels of total L1 mRNA expression and the number and spectrum of expressed L1 loci fluctuate with age as independent variables, demonstrating different patterns in different organs and sexes. Overall, our comparisons between organs and sexes and across ages ranging from 2 to 22 months establish previously unforeseen dynamic changes in L1 mRNA expression in vivo. These findings establish the beginning of an atlas of endogenous L1 mRNA expression across a broad range of biological variables that will guide future studies.

Elongation factor ELOF1 drives transcription-coupled repair and prevents genome instability

Correct transcription is crucial for life. However, DNA damage severely impedes elongating RNA polymerase II, causing transcription inhibition and transcription-replication conflicts. Cells are equipped with intricate mechanisms to counteract the severe consequence of these transcription-blocking lesions. However, the exact mechanism and factors involved remain largely unknown. Here, using a genome-wide CRISPR-Cas9 screen, we identified the elongation factor ELOF1 as an important factor in the transcription stress response following DNA damage. We show that ELOF1 has an evolutionarily conserved role in transcription-coupled nucleotide excision repair (TC-NER), where it promotes recruitment of the TC-NER factors UVSSA and TFIIH to efficiently repair transcription-blocking lesions and resume transcription. Additionally, ELOF1 modulates transcription to protect cells against transcription-mediated replication stress, thereby preserving genome stability. Thus, ELOF1 protects the transcription machinery from DNA damage via two distinct mechanisms.

Human L1 Transposition Dynamics Unraveled with Functional Data Analysis

Long INterspersed Elements-1 (L1s) constitute >17% of the human genome and still actively transpose in it. Characterizing L1 transposition across the genome is critical for understanding genome evolution and somatic mutations. However, to date, L1 insertion and fixation patterns have not been studied comprehensively. To fill this gap, we investigated three genome-wide data sets of L1s that integrated at different evolutionary times: 17,037 de novo L1s (from an L1 insertion cell-line experiment conducted in-house), and 1,212 polymorphic and 1,205 human-specific L1s (from public databases). We characterized 49 genomic features-proxying chromatin accessibility, transcriptional activity, replication, recombination, etc.-in the ±50 kb flanks of these elements. These features were contrasted between the three L1 data sets and L1-free regions using state-of-the-art Functional Data Analysis statistical methods, which treat high-resolution data as mathematical functions. Our results indicate that de novo, polymorphic, and human-specific L1s are surrounded by different genomic features acting at specific locations and scales. This led to an integrative model of L1 transposition, according to which L1s preferentially integrate into open-chromatin regions enriched in non-B DNA motifs, whereas they are fixed in regions largely free of purifying selection-depleted of genes and noncoding most conserved elements. Intriguingly, our results suggest that L1 insertions modify local genomic landscape by extending CpG methylation and increasing mononucleotide microsatellite density. Altogether, our findings substantially facilitate understanding of L1 integration and fixation preferences, pave the way for uncovering their role in aging and cancer, and inform their use as mutagenesis tools in genetic studies.

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.

A genome-wide strategy to identify causes and consequences of retrotransposon expression finds activation by BRCA1 in ovarian cancer

It is challenging to identify the causes and consequences of retrotransposon expression in human disease due to the hundreds of active genomic copies and their poor conservation across species. We profiled genomic insertions of retrotransposons in ovarian cancer. In addition, in ovarian and breast cancer we analyzed RNAs exhibiting Bayesian correlation with retrotransposon RNA to identify causes and consequences of retrotransposon expression. This strategy finds divergent inflammatory responses associated with retrotransposon expression in ovarian and breast cancer and identifies new factors inducing expression of endogenous retrotransposons including anti-viral responses and the common tumor suppressor BRCA1. In cell lines, mouse ovarian epithelial cells and patient-derived tumor spheroids, BRCA1 promotes accumulation of retrotransposon RNA. BRCA1 promotes transcription of active families of retrotransposons and their insertion into the genome. Intriguingly, elevated retrotransposon expression predicts survival in ovarian cancer patients. Retrotransposons are part of a complex regulatory network in ovarian cancer including BRCA1 that contributes to patient survival. The described strategy can be used to identify the regulators and impacts of retrotransposons in various contexts of biology and disease in humans.