Large Deletions, Cleavage of the Telomeric Repeat Sequence, and Reverse Transcriptase-Mediated DNA Damage Response Associated with Long Interspersed Element-1 ORF2p Enzymatic Activities.

L1 elements can cause DNA damage and genomic variation via retrotransposition and the generation of endonuclease-dependent DNA breaks. These processes require L1 ORF2p protein that contains an endonuclease domain, which cuts genomic DNA, and a reverse transcriptase domain, which synthesizes cDNA. The complete impact of L1 enzymatic activities on genome stability and cellular function remains understudied, and the spectrum of L1-induced mutations, other than L1 insertions, is mostly unknown. Using an inducible system, we demonstrate that an ORF2p containing functional reverse transcriptase is sufficient to elicit DNA damage response even in the absence of the functional endonuclease. Using a TK/Neo reporter system that captures misrepaired DNA breaks, we demonstrate that L1 expression results in large genomic deletions that lack any signatures of L1 involvement. Using an in vitro cleavage assay, we demonstrate that L1 endonuclease efficiently cuts telomeric repeat sequences. These findings support that L1 could be an unrecognized source of disease-promoting genomic deletions, telomere dysfunction, and an underappreciated source of chronic RT-mediated DNA damage response in mammalian cells. Our findings expand the spectrum of biological processes that can be triggered by functional and nonfunctional L1s, which have impactful evolutionary- and health-relevant consequences.

A novel role of TRIM28 B box domain in L1 retrotransposition and ORF2p-mediated cDNA synthesis.

The long interspersed element 1 (LINE-1 or L1) integration is affected by many cellular factors through various mechanisms. Some of these factors are required for L1 amplification, while others either suppress or enhance specific steps during L1 propagation. Previously, TRIM28 has been identified to suppress transposable elements, including L1 expression via its canonical role in chromatin remodeling. Here, we report that TRIM28 through its B box domain increases L1 retrotransposition and facilitates shorter cDNA and L1 insert generation in cultured cells. Consistent with the latter, we observe that tumor specific L1 inserts are shorter in endometrial, ovarian, and prostate tumors with higher TRIM28 mRNA expression than in those with lower TRIM28 expression. We determine that three amino acids in the B box domain that are involved in TRIM28 multimerization are critical for its effect on both L1 retrotransposition and cDNA synthesis. We provide evidence that B boxes from the other two members in the Class VI TRIM proteins, TRIM24 and TRIM33, also increase L1 retrotransposition. Our findings could lead to a better understanding of the host/L1 evolutionary arms race in the germline and their interplay during tumorigenesis.

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.

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.

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

BACKGROUND: Retrotransposons are one of the oldest evolutionary forces shaping mammalian genomes, with the ability to mobilize from one genomic location to another. This mobilization is also a significant factor in human disease. The only autonomous human retroelement, L1, has propagated to make up 17% of the human genome, accumulating over 500,000 copies. The majority of these loci are truncated or defective with only a few reported to remain capable of retrotransposition. We have previously published a strand-specific RNA-Seq bioinformatics approach to stringently identify at the locus-specific level the few expressed full-length L1s using cytoplasmic RNA. With growing repositories of RNA-Seq data, there is potential to mine these datasets to identify and study expressed L1s at single-locus resolution, although many datasets are not strand-specific or not generated from cytoplasmic RNA. RESULTS: We developed whole-cell, cytoplasmic and nuclear RNA-Seq datasets from 22Rv1 prostate cancer cells to test the influence of different preparations on the quality and effort needed to measure L1 expression. We found that there was minimal data loss in the identification of full-length expressed L1 s using whole cell, strand-specific RNA-Seq data compared to cytoplasmic, strand-specific RNA-Seq data. However, this was only possible with an increased amount of manual curation of the bioinformatics output to eliminate increased background. About half of the data was lost when the sequenced datasets were non-strand specific. CONCLUSIONS: The results of these studies demonstrate that with rigorous manual curation the utilization of stranded RNA-Seq datasets allow identification of expressed L1 loci from either cytoplasmic or whole-cell RNA-Seq datasets.

Influence of Daytime LED Light Exposure on Circadian Regulatory Dynamics of Metabolism and Physiology in Mice.

Light is a potent biologic force that profoundly influences circadian, neuroendocrine, and neurobehavioral regulation in animals. Previously we examined the effects of light-phase exposure of rats to white light-emitting diodes (LED), which emit more light in the blue-appearing portion of the visible spectrum (465 to 485 nm) than do broad-spectrum cool white fluorescent (CWF) light, on the nighttime melatonin amplitude and circadian regulation of metabolism and physiology. In the current studies, we tested the hypothesis that exposure to blue-enriched LED light at day (bLAD), compared with CWF, promotes the circadian regulation of neuroendocrine, metabolic, and physiologic parameters that are associated with optimizing homeostatic regulation of health and wellbeing in 3 mouse strains commonly used in biomedical research (C3H [melatonin-producing], C57BL/6, and BALB/c [melatonin-non-producing]). Compared with male and female mice housed for 12 wk under 12:12-h light:dark (LD) cycles in CWF light, C3H mice in bLAD evinced 6-fold higher peak plasma melatonin levels at the middark phase; in addition, high melatonin levels were prolonged 2 to 3 h into the light phase. C57BL/6 and BALB/c strains did not produce nighttime pineal melatonin. Body growth rates; dietary and water intakes; circadian rhythms of arterial blood corticosterone, insulin, leptin, glucose, and lactic acid; pO2 and pCO2; fatty acids; and metabolic indicators (cAMP, DNA, tissue DNA 3H-thymidine incorporation, fat content) in major organ systems were significantly lower and activation of major metabolic signaling pathways (mTOR, GSK3ß, and SIRT1) in skeletal muscle and liver were higher only in C3H mice in bLAD compared with CWF. These data show that exposure of C3H mice to bLAD compared with CWF has a marked positive effect on the circadian regulation of neuroendocrine, metabolic, and physiologic parameters associated with the promotion of animal health and wellbeing that may influence scientific outcomes. The absence of enhancement in amelatonic strains suggests hyperproduction of nighttime melatonin may be a key component of the physiology.

RNA Next-Generation Sequencing and a Bioinformatics Pipeline to Identify Expressed LINE-1s at the Locus-Specific Level.

Long INterspersed Elements-1 (LINEs/L1s) are repetitive elements that can copy and randomly insert in the genome resulting in genomic instability and mutagenesis. Understanding the expression patterns of L1 loci at the individual level will lend to the understanding of the biology of this mutagenic element. This autonomous element makes up a significant portion of the human genome with over 500,000 copies, though 99% are truncated and defective. However, their abundance and dominant number of defective copies make it challenging to identify authentically expressed L1s from L1-related sequences expressed as part of other genes. It is also challenging to identify which specific L1 locus is expressed due to the repetitive nature of the elements. Overcoming these challenges, we present an RNA-Seq bioinformatic approach to identify L1 expression at the locus specific level. In summary, we collect cytoplasmic RNA, select for polyadenylated transcripts, and utilize strand-specific RNA-Seq analyses to uniquely map reads to L1 loci in the human reference genome. We visually curate each L1 locus with uniquely mapped reads to confirm transcription from its own promoter and adjust mapped transcript reads to account for mappability of each individual L1 locus. This approach was applied to a prostate tumor cell line, DU145, to demonstrate the ability of this protocol to detect expression from a small number of the full-length L1 elements.

Mouse germ line mutations due to retrotransposon insertions.

Transposable element (TE) insertions are responsible for a significant fraction of spontaneous germ line mutations reported in inbred mouse strains. This major contribution of TEs to the mutational landscape in mouse contrasts with the situation in human, where their relative contribution as germ line insertional mutagens is much lower. In this focussed review, we provide comprehensive lists of TE-induced mouse mutations, discuss the different TE types involved in these insertional mutations and elaborate on particularly interesting cases. We also discuss differences and similarities between the mutational role of TEs in mice and humans.

Effect of Daytime Blue-enriched LED Light on the Nighttime Circadian Melatonin Inhibition of Hepatoma 7288CTC Warburg Effect and Progression.

Liver cancer is the second leading cause of cancer death worldwide. Metabolic pathways within the liver and liver cancers are highly regulated by the central circadian clock in the suprachiasmatic nuclei (SCN). Daily light and dark cycles regulate the SCN-driven pineal production of the circadian anticancer hormone melatonin and temporally coordinate circadian rhythms of metabolism and physiology in mammals. In previous studies, we demonstrated that melatonin suppresses linoleic acid metabolism and the Warburg effect (aerobic glycolysis)in human breast cancer xenografts and that blue-enriched light (465-485 nm) from light-emitting diode lighting at daytime (bLAD) amplifies nighttime circadian melatonin levels in rats by 7-fold over cool white fluorescent (CWF) lighting. Here we tested the hypothesis that daytime exposure of tissue-isolated Morris hepatoma 7288CTC-bearing male rats to bLAD amplifies the nighttime melatonin signal to enhance the inhibition of tumor growth. Compared with rats housed under a 12:12-h light:dark cycle in CWF light, rats in bLAD light evinced a 7-fold higher peak plasma melatonin level at the mid-dark phase; in addition, high melatonin levels were prolonged until 4 h into the light phase. After implantation of tissue-isolated hepatoma 7288CTC xenografts, tumor growth rates were markedly delayed, and tumor cAMP levels, LA metabolism, the Warburg effect, and growth signaling activities were decreased in rats in bLAD compared with CWF daytime lighting. These data show that the increased nighttime circadian melatonin levels due to bLAD exposure decreases hepatoma metabolic, signaling, and proliferative activities beyond what occurs after normal melatonin signaling under CWF light.

Evaluating different DNA binding domains to modulate L1 ORF2p-driven site-specific retrotransposition events in human cells.

DNA binding domains (DBDs) have been used with great success to impart targeting capabilities to a variety of proteins creating highly useful genomic tools. We evaluated the ability of five types of DBDs and strategies (AAV Rep proteins, Cre, TAL effectors, zinc finger proteins, and Cas9/gRNA system) to target the L1 ORF2 protein to drive retrotransposition of Alu inserts to specific sequences in the human genome. First, we find that the L1 ORF2 protein tolerates the addition of protein domains both at the amino- and carboxy-terminus. Although in some instances retrotransposition efficiencies slightly diminished, all fusion proteins containing an intact ORF2 were capable of driving retrotransposition. Second, the stability of individual ORF2 fusion proteins varies and difficult to predict. Third, DBDs that require the formation of multimers for target recognition are unlikely to modify targeting of ORF2p-driven insertions. Fourth, the more components needed to assemble into a complex to drive targeted retrotransposition, the less likely the strategy will increase targeted insertions. Fifth, abundance of target sequences present in the genome will likely dictate the effectiveness and efficiency of targeted insertions. Lastly, the cleavage capabilities of Cas9 (or a Cas9 nickase variant) are unable to substitute for the L1 ORF2 endonuclease domain functions, suggestive that the endonuclease domain has alternate functions needed for retrotransposition. From these studies, we conclude that the most critical component for the modification of the human L1 ORF2 protein to drive targeted insertions is the selection of the DBD due to the varying functional requirements and impacts on protein stability.

Truncated ORF1 proteins can suppress LINE-1 retrotransposition in trans.

Long interspersed element 1 (L1) is an autonomous non-LTR retroelement that is active in mammalian genomes. Although retrotranspositionally incompetent and functional L1 loci are present in the same genomes, it remains unknown whether non-functional L1s have any trans effect on mobilization of active elements. Using bioinformatic analysis, we identified over a thousand of human L1 loci containing at least one stop codon in their ORF1 sequence. RNAseq analysis confirmed that many of these loci are expressed. We demonstrate that introduction of equivalent stop codons in the full-length human L1 sequence leads to the expression of truncated ORF1 proteins. When supplied in trans some truncated human ORF1 proteins suppress human L1 retrotransposition. This effect requires the N-terminus and coiled-coil domain (C-C) as mutations within the ORF1p C-C domain abolish the suppressive effect of truncated proteins on L1 retrotransposition. We demonstrate that the expression levels and length of truncated ORF1 proteins influence their ability to suppress L1 retrotransposition. Taken together these findings suggest that L1 retrotransposition may be influenced by coexpression of defective L1 loci and that these L1 loci may reduce accumulation of de novo L1 integration events.

Involvement of Conserved Amino Acids in the C-Terminal Region of LINE-1 ORF2p in Retrotransposition.

Long interspersed element 1 (L1) is the only currently active autonomous retroelement in the human genome. Along with the parasitic SVA and short interspersed element Alu, L1 is the source of DNA damage induced by retrotransposition: a copy-and-paste process that has the potential to disrupt gene function and cause human disease. The retrotransposition process is dependent upon the ORF2 protein (ORF2p). However, it is unknown whether most of the protein is important for retrotransposition. In particular, other than the Cys motif, the C terminus of the protein has not been intensely examined in the context of retrotransposition. Using evolutionary analysis and the Alu retrotransposition assay, we sought to identify additional amino acids in the C terminus important for retrotransposition. Here, we demonstrate that Gal4-tagged and untagged C-terminally truncated ORF2p fragments possess residual potential to drive Alu retrotransposition. Using sight-directed mutagenesis we identify that while the Y1180 amino acid is important for ORF2p- and L1-driven Alu retrotransposition, a mutation at this position improves L1 retrotransposition. Even though the mechanism of the contribution of Y1180 to Alu and L1 mobilization remains unknown, experimental evidence rules out its direct involvement in the ability of the ORF2p reverse transcriptase to generate complementary DNA. Additionally, our data support that ORF2p amino acids 1180 and 1250-1262 may be involved in the reported ORF1p-mediated increase in ORF2p-driven Alu retrotransposition.

The Nucleotide Excision Repair Pathway Limits L1 Retrotransposition.

Long interspersed elements 1 (L1) are active mobile elements that constitute almost 17% of the human genome. They amplify through a "copy-and-paste" mechanism termed retrotransposition, and de novo insertions related to these elements have been reported to cause 0.2% of genetic diseases. Our previous data demonstrated that the endonuclease complex ERCC1-XPF, which cleaves a 3' DNA flap structure, limits L1 retrotransposition. Although the ERCC1-XPF endonuclease participates in several different DNA repair pathways, such as single-strand annealing, or in telomere maintenance, its recruitment to DNA lesions is best characterized in the nucleotide excision repair (NER) pathway. To determine if the NER pathway prevents the insertion of retroelements in the genome, we monitored the retrotransposition efficiencies of engineered L1 elements in NER-deficient cells and in their complemented versions. Core proteins of the NER pathway, XPD and XPA, and the lesion binding protein, XPC, are involved in limiting L1 retrotransposition. In addition, sequence analysis of recovered de novo L1 inserts and their genomic locations in NER-deficient cells demonstrated the presence of abnormally large duplications at the site of insertion, suggesting that NER proteins may also play a role in the normal L1 insertion process. Here, we propose new functions for the NER pathway in the maintenance of genome integrity: limitation of insertional mutations caused by retrotransposons and the prevention of potentially mutagenic large genomic duplications at the site of retrotransposon insertion events.

A comprehensive approach to expression of L1 loci.

L1 elements represent the only currently active, autonomous retrotransposon in the human genome, and they make major contributions to human genetic instability. The vast majority of the 500 000 L1 elements in the genome are defective, and only a relatively few can contribute to the retrotransposition process. However, there is currently no comprehensive approach to identify the specific loci that are actively transcribed separate from the excess of L1-related sequences that are co-transcribed within genes. We have developed RNA-Seq procedures, as well as a 1200 bp 5΄ RACE product coupled with PACBio sequencing that can identify the specific L1 loci that contribute most of the L1-related RNA reads. At least 99% of L1-related sequences found in RNA do not arise from the L1 promoter, instead representing pieces of L1 incorporated in other cellular RNAs. In any given cell type a relatively few active L1 loci contribute to the 'authentic' L1 transcripts that arise from the L1 promoter, with significantly different loci seen expressed in different tissues.

Effects of Daytime Exposure to Light from Blue-Enriched Light-Emitting Diodes on the Nighttime Melatonin Amplitude and Circadian Regulation of Rodent Metabolism and Physiology.

Regular cycles of exposure to light and dark control pineal melatonin production and temporally coordinate circadian rhythms of metabolism and physiology in mammals. Previously we demonstrated that the peak circadian amplitude of nocturnal blood melatonin levels of rats were more than 6-fold higher after exposure to cool white fluorescent (CWF) light through blue-tinted (compared with clear) rodent cages. Here, we evaluated the effects of light-phase exposure of rats to white light-emitting diodes (LED), which emit light rich in the blue-appearing portion of the visible spectrum (465-485 nm), compared with standard broadspectrum CWF light, on melatonin levels during the subsequent dark phase and on plasma measures of metabolism and physiology. Compared with those in male rats under a 12:12-h light:dark cycle in CWF light, peak plasma melatonin levels at the middark phase (time, 2400) in rats under daytime LED light were over 7-fold higher, whereas midlight phase levels (1200) were low in both groups. Food and water intakes, body growth rate, and total fatty acid content of major metabolic tissues were markedly lower, whereas protein content was higher, in the LED group compared with CWF group. Circadian rhythms of arterial plasma levels of total fatty acids, glucose, lactic acid, pO2, pCO2, insulin, leptin, and corticosterone were generally lower in LED-exposed rats. Therefore, daytime exposure of rats to LED light with high blue emissions has a marked positive effect on the circadian regulation of neuroendocrine, metabolic, and physiologic parameters associated with the promotion of animal health and wellbeing and thus may influence scientific outcomes.

The importance of L1 ORF2p cryptic sequence to ORF2p fragment-mediated cytotoxicity.

The Long Interspersed Element 1 (LINE1 or L1) ORF2 protein (ORF2p) can cause DNA damage through the activity of its endonuclease domain (EN). The DNA double-strand breaks (DSB) introduced by the ORF2p EN have the potential to be mutagenic. Previously, our lab has shown that ORF2p fragments containing the EN domain could be expressed in mammalian cells and have variable cytotoxicity. Inclusion of the ORF2p sequence C-terminal to the EN domain in these fragments both reduced the cytotoxicity of these fragments and increased their presence in the nucleus as detected by Western blot analysis. Here, we identify the amino acids (aa 270-274) in the newly-identified ORF2p Cryptic region (Cry) that may be important to the subcellular localization and cytotoxic potential of these EN-containing ORF2p fragments.

The endonuclease domain of the LINE-1 ORF2 protein can tolerate multiple mutations.

BACKGROUND: Approximately 17 % of the human genome is comprised of the Long INterspersed Element-1 (LINE-1 or L1) retrotransposon, the only currently active autonomous family of retroelements. Though L1 elements have helped to shape mammalian genome evolution over millions of years, L1 activity can also be mutagenic and result in human disease. L1 expression has the potential to contribute to genomic instability via retrotransposition and DNA double-strand breaks (DSBs). Additionally, L1 is responsible for structural genomic variations induced by other transposable elements such as Alu and SVA, which rely on the L1 ORF2 protein for their propagation. Most of the genomic damage associated with L1 activity originates with the endonuclease domain of the ORF2 protein, which nicks the DNA in preparation for target-primed reverse transcription. RESULTS: Bioinformatic analysis of full-length L1 loci residing in the human genome identified numerous mutations in the amino acid sequence of the ORF2 endonuclease domain. Some of these mutations were found in residues which were predicted to be phosphorylation sites for cellular kinases. We mutated several of these putative phosphorylation sites in the ORF2 endonuclease domain and investigated the effect of these mutations on the function of the full-length ORF2 protein and the endonuclease domain (ENp) alone. Most of the single and multiple point mutations that were tested did not significantly impact expression of the full-length ORF2p, or alter its ability to drive Alu retrotransposition. Similarly, most of those same mutations did not significantly alter expression of ENp, or impair its ability to induce DNA damage and cause toxicity. CONCLUSIONS: Overall, our data demonstrate that the full-length ORF2p or the ENp alone can tolerate several specific single and multiple point mutations in the endonuclease domain without significant impairment of their ability to support Alu mobilization or induce DNA damage, respectively.