Comparative analysis on the expression of L1 loci using various RNA-Seq preparations

Repetitive element transcript accumulation is associated with inflammaging in humans

Chronic, low-grade inflammation increases with aging, contributing to functional declines and diseases that reduce healthspan. Growing evidence suggests that transcripts from repetitive elements (RE) in the genome contribute to this "inflammaging" by stimulating innate immune activation, but evidence of RE-associated inflammation with aging in humans is limited. Here, we present transcriptomic and clinical data showing that RE transcript levels are positively related to gene expression of innate immune sensors, and to serum interleukin 6 (a marker of systemic inflammation), in a large group of middle-aged and older adults. We also: (1) use transcriptomics and whole-genome bisulfite (methylation) sequencing to show that many RE may be hypomethylated with aging, and that aerobic exercise, a healthspan-extending intervention, reduces RE transcript levels and increases RE methylation in older adults; and (2) extend our findings in a secondary dataset demonstrating age-related changes in RE chromatin accessibility. Collectively, our data support the idea that age-related RE transcript accumulation may play a role in inflammaging in humans, and that RE dysregulation with aging may be due in part to upstream epigenetic changes.

TE-Seq: A Transposable Element Annotation and RNA-Seq Pipeline

Motivation

The recognition that transposable elements (TEs) play important roles in many biological processes has elicited growing interest in analyzing sequencing data derived from this 'dark genome'. This is however complicated by the highly repetitive nature of these sequences in genomes, requiring the deployment of several problem-specific tools as well as the curation of appropriate genome annotations. This pipeline aims to make the analysis of TE sequences and their expression more generally accessible.

Results

The TE-Seq pipeline conducts an end-to-end analysis of RNA sequencing data, examining both genes and TEs. It implements the most current computational methods tailor-made for TEs, and produces a comprehensive analysis of TE expression at both the level of the individual element and at the TE clade level. Furthermore, if supplied with long-read DNA sequencing data, it is able to assess TE expression from non-reference (polymorphic) loci. As a demonstration, we analyzed proliferating, early senescent, and late senescent lung fibroblast RNA-Seq data, and created a custom reference genome and annotations for this cell strain using Nanopore sequencing data. We found that several retrotransposable element (RTE) clades were upregulated in senescence, which included non-reference, intact, and potentially active elements.

Availability and implementation

TE-Seq is made available as a Snakemake pipeline which can be obtained at https://github.com/maxfieldk/TE-Seq. All software dependencies besides Snakemake and Docker/Singularity are packaged into a container which is automatically built and deployed by the pipeline at runtime.

LINE-1 transcription activates long-range gene expression

Long interspersed nuclear element-1 (LINE-1 or L1) is a retrotransposon group that constitutes 17% of the human genome and shows variable expression across cell types. However, the control of L1 expression and its function in gene regulation are incompletely understood. Here we show that L1 transcription activates long-range gene expression. Genome-wide CRISPR-Cas9 screening using a reporter driven by the L1 5' UTR in human cells identifies functionally diverse genes affecting L1 expression. Unexpectedly, altering L1 expression by knockout of regulatory genes impacts distant gene expression. L1s can physically contact their distal target genes, with these interactions becoming stronger upon L1 activation and weaker when L1 is silenced. Remarkably, L1s contact and activate genes essential for zygotic genome activation (ZGA), and L1 knockdown impairs ZGA, leading to developmental arrest in mouse embryos. These results characterize the regulation and function of L1 in long-range gene activation and reveal its importance in mammalian ZGA.

Cancer Relevance of Circulating Antibodies Against LINE-1 Antigens in Humans

Long interspersed nuclear element-1 (LINE-1 or L1), the most abundant family of autonomous retrotransposons occupying over 17% of human DNA, is epigenetically silenced in normal tissues by the mechanisms involving p53 but is frequently derepressed in cancer, suggesting that L1-encoded proteins may act as tumor-associated antigens recognized by the immune system. In this study, we established an immunoassay to detect circulating autoantibodies against L1 proteins in human blood. Using this assay in >2,800 individuals with or without cancer, we observed significantly higher IgG titers against L1-encoded ORF1p and ORF2p in patients with lung, pancreatic, ovarian, esophageal, and liver cancers than in healthy individuals. Remarkably, elevated levels of anti-ORF1p-reactive IgG were observed in patients with cancer with disease stages 1 and 2, indicating that the immune response to L1 antigens can occur in the early phases of carcinogenesis. We concluded that the antibody response against L1 antigens could contribute to the diagnosis and determination of immunoreactivity of tumors among cancer types that frequently escape early detection.

SIGNIFICANCE

The discovery of autoantibodies against antigens encoded by L1 retrotransposons in patients with five poorly curable cancer types has potential implications for the detection of an ongoing carcinogenic process and tumor immunoreactivity.

Nuclei on the Rise: When Nuclei-Based Methods Meet Next-Generation Sequencing

In the last decade, we have witnessed an upsurge in nuclei-based studies, particularly coupled with next-generation sequencing. Such studies aim at understanding the molecular states that exist in heterogeneous cell populations by applying increasingly more affordable sequencing approaches, in addition to optimized methodologies developed to isolate and select nuclei. Although these powerful new methods promise unprecedented insights, it is important to understand and critically consider the associated challenges. Here, we provide a comprehensive overview of the rise of nuclei-based studies and elaborate on their advantages and disadvantages, with a specific focus on their utility for transcriptomic sequencing analyses. Improved designs and appropriate use of the various experimental strategies will result in acquiring biologically accurate and meaningful information.

SCIFER: approach for analysis of LINE-1 mRNA expression in single cells at a single locus resolution

BACKGROUND

Endogenous expression of L1 mRNA is the first step in an L1-initiated mutagenesis event. However, the contribution of individual cell types to patterns of organ-specific L1 mRNA expression remains poorly understood, especially at single-locus resolution. We introduce a method to quantify expression of mobile elements at the single-locus resolution in scRNA-Seq datasets called Single Cell Implementation to Find Expressed Retrotransposons (SCIFER). SCIFER aligns scRNA-Seq reads uniquely to the genome and extracts alignments from single cells by cell-specific barcodes. In contrast to the alignment performed using default parameters, this alignment strategy increases accuracy of L1 locus identification by retaining only reads that are uniquely mapped to individual L1 loci. L1 loci expressed in single cells are unambiguously identified using a list of L1 loci manually validated to be expressed in bulk RNA-Seq datasets generated from the same cell line or organ.

RESULTS

Validation of SCIFER using MCF7 cells determined technical parameters needed for optimal detection of L1 expression in single cells. We show that unsupervised analysis of L1 expression in single cells exponentially inflates both the levels of L1 expression and the number of expressed L1 loci. Application of SCIFER to analysis of scRNA-Seq datasets generated from mouse and human testes identified that mouse Round Spermatids and human Spermatogonia, Spermatocytes, and Round Spermatids express the highest levels of L1 mRNA. Our analysis also determined that similar to mice, human testes from unrelated individuals share as much as 80% of expressed L1 loci. Additionally, SCIFER determined that individual mouse cells co-express different L1 sub-families and different families of transposable elements, experimentally validating their co-existence in the same cell.

CONCLUSIONS

SCIFER detects mRNA expression of individual L1 loci in single cells. It is compatible with scRNA-Seq datasets prepared using traditional sequencing methods. Validated using a human cancer cell line, SCIFER analysis of mouse and human testes identified key cell types supporting L1 expression in these species. This will further our understanding of differences and similarities in endogenous L1 mRNA expression patterns in mice and humans.

Urinary extracellular vesicles contain mature transcriptome enriched in circular and long noncoding RNAs with functional significance in prostate cancer

Long noncoding (lnc)RNAs modulate gene expression alongside presenting unexpected source of neoantigens. Despite their immense interest, their ability to be transferred and control adjacent cells is unknown. Extracellular Vesicles (EVs) offer a protective environment for nucleic acids, with pro and antitumourigenic functions by controlling the immune response. In contrast to extracellular nonvesicular RNA, few studies have addressed the full RNA content within human fluids' EVs and have compared them with their tissue of origin. Here, we performed Total RNA-Sequencing on six Formalin-Fixed-Paraffin-Embedded (FFPE) prostate cancer (PCa) tumour tissues and their paired urinary (u)EVs to provide the first whole transcriptome comparison from the same patients. UEVs contain simplified transcriptome with intron-free cytoplasmic transcripts and enriched lnc/circular (circ)RNAs, strikingly common to an independent 20 patients' urinary cohort. Our full cellular and EVs transcriptome comparison within three PCa cell lines identified a set of overlapping 14 uEV-circRNAs characterized as essential for prostate cell proliferation in vitro and 28 uEV-lncRNAs belonging to the cancer-related lncRNA census (CLC2). In addition, we found 15 uEV-lncRNAs, predicted to encode 768 high-affinity neoantigens, and for which three of the encoded-ORF produced detectable unmodified peptides by mass spectrometry. Our dual analysis of EVs-lnc/circRNAs both in urines' and in vitro's EVs provides a fundamental resource for future uEV-lnc/circRNAs phenotypic characterization involved in PCa.

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.

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.

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.

Transcriptomic analysis of Entamoeba histolytica reveals domain-specific sense strand expression of LINE-encoded ORFs with massive antisense expression of RT domain

LINEs are retrotransposable elements found in diverse organisms. Their activity is kept in check by several mechanisms, including transcriptional silencing. Here we have analyzed the transcription status of LINE1 copies in the early-branching parasitic protist Entamoeba histolytica. Full-length EhLINE1 encodes ORF1, and ORF2 with reverse transcriptase (RT) and endonuclease (EN) domains. RNA-Seq analysis of EhLINE1 copies (both truncated and full-length) showed unique features. Firstly, although 20/41 transcribed copies were full-length, we failed to detect any full-length transcripts. Rather, sense-strand transcripts mapped to the functional domains- ORF1, RT and EN. Secondly, there was strong antisense transcription specifically from RT domain. No antisense transcripts were seen from ORF1. Antisense RT transcripts did not encode known functional peptides. They could possibly be involved in attenuating translation of RT domain, as we failed to detect ORF2p, whereas ORF1p was detectable. Lack of full-length transcripts and strong antisense RT expression may serve to limit EhLINE1 retrotransposition.

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.

Repetitive elements as a transcriptomic marker of aging: Evidence in multiple datasets and models

Transcriptomic markers of aging can be useful for studying age‐related processes and diseases. However, noncoding repetitive element (RE) transcripts, which may play an important role in aging, are commonly overlooked in transcriptome studies—and their potential as a transcriptomic marker of aging has not been evaluated. Here, we used multiple RNA‐seq datasets generated from human samples and Caenorhabditis elegans and found that most RE transcripts (a) accumulate progressively with aging; (b) can be used to accurately predict age; and (c) may be a good marker of biological age. The strong RE/aging correlations we observed are consistent with growing evidence that RE transcripts contribute directly to aging and disease.

Measuring and interpreting transposable element expression

Transposable elements (TEs) are insertional mutagens that contribute greatly to the plasticity of eukaryotic genomes, influencing the evolution and adaptation of species as well as physiology or disease in individuals. Measuring TE expression helps to understand not only when and where TE mobilization can occur but also how this process alters gene expression, chromatin accessibility or cellular signalling pathways. Although genome-wide gene expression assays such as RNA sequencing include transposon-derived transcripts, most computational analytical tools discard or misinterpret TE-derived reads. Emerging approaches are improving the identification of expressed TE loci and helping to discriminate TE transcripts that permit TE mobilization from chimeric gene-TE transcripts or pervasive transcription. Here we review the main challenges associated with the detection of TE expression, including mappability, insertional and internal sequence polymorphisms, and the diversity of the TE transcriptional landscape, as well as the different experimental and computational strategies to solve them.