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.

Altered DNA repair creates novel Alu/Alu repeat-mediated deletions.

Alu elements are the most abundant source of nonallelic homology that influences genetic instability in the human genome. When there is a DNA double-stranded break, the Alu element's high copy number, moderate length and distance and mismatch between elements uniquely influence recombination processes. We utilize a reporter-gene assay to show the complex influence of Alu mismatches on Alu-related repeat-mediated deletions (RMDs). The Alu/Alu heteroduplex intermediate can result in a nonallelic homologous recombination (HR). Alternatively, the heteroduplex can result in various DNA breaks around the Alu elements caused by competing nucleases. These breaks can undergo Alt-nonhomologous end joining to cause deletions focused around the Alu elements. Formation of these heteroduplex intermediates is largely RAD52 dependent. Cells with low ERCC1 levels utilize more of these alternatives resolutions, while cells with MSH2 defects tend to have more RMDs with a specific increase in the HR events. Therefore, Alu elements are expected to create different forms of deletions in various cancers depending on a number of these DNA repair defects.

Long-Distance Relationships: Suppression of Repeat-Mediated Deletions.

The high proportion of repetitive DNA sequences in the human genome provides tremendous opportunities for DNA rearrangements between non-allelic repetitive elements. The genome must use multiple competing and collaborating repair mechanisms to minimize these types of DNA rearrangements, some of which fail in cancer cells where DNA repair pathways are suppressed.

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.

The contribution of alu elements to mutagenic DNA double-strand break repair.

Alu elements make up the largest family of human mobile elements, numbering 1.1 million copies and comprising 11% of the human genome. As a consequence of evolution and genetic drift, Alu elements of various sequence divergence exist throughout the human genome. Alu/Alu recombination has been shown to cause approximately 0.5% of new human genetic diseases and contribute to extensive genomic structural variation. To begin understanding the molecular mechanisms leading to these rearrangements in mammalian cells, we constructed Alu/Alu recombination reporter cell lines containing Alu elements ranging in sequence divergence from 0%-30% that allow detection of both Alu/Alu recombination and large non-homologous end joining (NHEJ) deletions that range from 1.0 to 1.9 kb in size. Introduction of as little as 0.7% sequence divergence between Alu elements resulted in a significant reduction in recombination, which indicates even small degrees of sequence divergence reduce the efficiency of homology-directed DNA double-strand break (DSB) repair. Further reduction in recombination was observed in a sequence divergence-dependent manner for diverged Alu/Alu recombination constructs with up to 10% sequence divergence. With greater levels of sequence divergence (15%-30%), we observed a significant increase in DSB repair due to a shift from Alu/Alu recombination to variable-length NHEJ which removes sequence between the two Alu elements. This increase in NHEJ deletions depends on the presence of Alu sequence homeology (similar but not identical sequences). Analysis of recombination products revealed that Alu/Alu recombination junctions occur more frequently in the first 100 bp of the Alu element within our reporter assay, just as they do in genomic Alu/Alu recombination events. This is the first extensive study characterizing the influence of Alu element sequence divergence on DNA repair, which will inform predictions regarding the effect of Alu element sequence divergence on both the rate and nature of DNA repair events.

Sequencing, identification and mapping of primed L1 elements (SIMPLE) reveals significant variation in full length L1 elements between individuals.

BACKGROUND: There are over a half a million copies of L1 retroelements in the human genome which are responsible for as much as 0.5% of new human genetic diseases. Most new L1 inserts arise from young source elements that are polymorphic in the human genome. Highly active polymorphic "hot" L1 source elements have been shown to be capable of extremely high levels of mobilization and result in numerous instances of disease. Additionally, hot polymorphic L1s have been described to be highly active within numerous cancer genomes. These hot L1s result in mutagenesis by insertion of new L1 copies elsewhere in the genome, but also have been shown to generate additional full length L1 insertions which are also hot and able to further retrotranspose. Through this mechanism, hot L1s may amplify within a tumor and result in a continued cycle of mutagenesis. RESULTS AND CONCLUSIONS: We have developed a method to detect full-length, polymorphic L1 elements using a targeted next generation sequencing approach, Sequencing Identification and Mapping of Primed L1 Elements (SIMPLE). SIMPLE has 94% sensitivity and detects nearly all full-length L1 elements in a genome. SIMPLE will allow researchers to identify hot mutagenic full-length L1s as potential drivers of genome instability. Using SIMPLE we find that the typical individual has approximately 100 non-reference, polymorphic L1 elements in their genome. These elements are at relatively low population frequencies relative to previously identified polymorphic L1 elements and demonstrate the tremendous diversity in potentially active L1 elements in the human population.

Inferring the expression variability of human transposable element-derived exons by linear model analysis of deep RNA sequencing data.

BACKGROUND: The exonization of transposable elements (TEs) has proven to be a significant mechanism for the creation of novel exons. Existing knowledge of the retention patterns of TE exons in mRNAs were mainly established by the analysis of Expressed Sequence Tag (EST) data and microarray data. RESULTS: This study seeks to validate and extend previous studies on the expression of TE exons by an integrative statistical analysis of high throughput RNA sequencing data. We collected 26 RNA-seq datasets spanning multiple tissues and cancer types. The exon-level digital expressions (indicating retention rates in mRNAs) were quantified by a double normalized measure, called the rescaled RPKM (Reads Per Kilobase of exon model per Million mapped reads). We analyzed the distribution profiles and the variability (across samples and between tissue/disease groups) of TE exon expressions, and compared them with those of other constitutive or cassette exons. We inferred the effects of four genomic factors, including the location, length, cognate TE family and TE nucleotide proportion (RTE, see Methods section) of a TE exon, on the exons' expression level and expression variability. We also investigated the biological implications of an assembly of highly-expressed TE exons. CONCLUSION: Our analysis confirmed prior studies from the following four aspects. First, with relatively high expression variability, most TE exons in mRNAs, especially those without exact counterparts in the UCSC RefSeq (Reference Sequence) gene tables, demonstrate low but still detectable expression levels in most tissue samples. Second, the TE exons in coding DNA sequences (CDSs) are less highly expressed than those in 3' (5') untranslated regions (UTRs). Third, the exons derived from chronologically ancient repeat elements, such as MIRs, tend to be highly expressed in comparison with those derived from younger TEs. Fourth, the previously observed negative relationship between the lengths of exons and the inclusion levels in transcripts is also true for exonized TEs. Furthermore, our study resulted in several novel findings. They include: (1) for the TE exons with non-zero expression and as shown in most of the studied biological samples, a high TE nucleotide proportion leads to their lower retention rates in mRNAs; (2) the considered genomic features (i.e. a continuous variable such as the exon length or a category indicator such as 3'UTR) influence the expression level and the expression variability (CV) of TE exons in an inverse manner; (3) not only the exons derived from Alu elements but also the exons from the TEs of other families were preferentially established in zinc finger (ZNF) genes.

RNA truncation by premature polyadenylation attenuates human mobile element activity.

Long interspersed elements (LINE-1s, also called L1s) are the only active members of the autonomous, non-long terminal repeat (LTR) retrotransposon family, which reshapes mammalian genomes in many different ways. LINE-1 expression is low in most differentiated cells but high in some cancer cells, in testis and during embryonic development. To minimize the negative impact on their hosts' genomes, many mobile elements strategically limit their amplification potential, particularly in somatic cells. Here we show that the A-rich coding strand of the human LINE-1 contains multiple functional canonical and noncanonical polyadenylation (poly(A)) signals, resulting in truncation of full-length transcripts by premature polyadenylation. This attenuation lowers the rate of retrotransposition in assays using HeLa cells. It probably also increases the negative effects of LINE-1 insertions into genes.

Dominant and recessive deafness caused by mutations of a novel gene, TMC1, required for cochlear hair-cell function.

Positional cloning of hereditary deafness genes is a direct approach to identify molecules and mechanisms underlying auditory function. Here we report a locus for dominant deafness, DFNA36, which maps to human chromosome 9q13-21 in a region overlapping the DFNB7/B11 locus for recessive deafness. We identified eight mutations in a new gene, transmembrane cochlear-expressed gene 1 (TMC1), in a DFNA36 family and eleven DFNB7/B11 families. We detected a 1.6-kb genomic deletion encompassing exon 14 of Tmc1 in the recessive deafness (dn) mouse mutant, which lacks auditory responses and has hair-cell degeneration. TMC1 and TMC2 on chromosome 20p13 are members of a gene family predicted to encode transmembrane proteins. Tmc1 mRNA is expressed in hair cells of the postnatal mouse cochlea and vestibular end organs and is required for normal function of cochlear hair cells.

Somatic expression of LINE-1 elements in human tissues.

LINE-1 expression damages host DNA via insertions and endonuclease-dependent DNA double-strand breaks (DSBs) that are highly toxic and mutagenic. The predominant tissue of LINE-1 expression has been considered to be the germ line. We show that both full-length and processed L1 transcripts are widespread in human somatic tissues and transformed cells, with significant variation in both L1 expression and L1 mRNA processing. This is the first demonstration that RNA processing is a major regulator of L1 activity. Many tissues also produce translatable spliced transcript (SpORF2). An Alu retrotransposition assay, COMET assays and 53BP1 foci staining show that the SpORF2 product can support functional ORF2 protein expression and can induce DNA damage in normal cells. Tests of the senescence-associated beta-galactosidase expression suggest that expression of exogenous full-length L1, or the SpORF2 mRNA alone in human fibroblasts and adult stem cells triggers a senescence-like phenotype, which is one of the reported responses to DNA damage. In contrast to previous assumptions that L1 expression is germ line specific, the increased spectrum of tissues exposed to L1-associated damage suggests a role for L1 as an endogenous mutagen in somatic tissues. These findings have potential consequences for the whole organism in the form of cancer and mammalian aging.

All y'all need to know 'bout retroelements in cancer.

Genetic instability is one of the principal hallmarks and causative factors in cancer. Human transposable elements (TE) have been reported to cause human diseases, including several types of cancer through insertional mutagenesis of genes critical for preventing or driving malignant transformation. In addition to retrotransposition-associated mutagenesis, TEs have been found to contribute even more genomic rearrangements through non-allelic homologous recombination. TEs also have the potential to generate a wide range of mutations derivation of which is difficult to directly trace to mobile elements, including double strand breaks that may trigger mutagenic genomic rearrangements. Genome-wide hypomethylation of TE promoters and significantly elevated TE expression in almost all human cancers often accompanied by the loss of critical DNA sensing and repair pathways suggests that the negative impact of mobile elements on genome stability should increase as human tumors evolve. The biological consequences of elevated retroelement expression, such as the rate of their amplification, in human cancers remain obscure, particularly, how this increase translates into disease-relevant mutations. This review is focused on the cellular mechanisms that control human TE-associated mutagenesis in cancer and summarizes the current understanding of TE contribution to genetic instability in human malignancies.

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.

Alu distribution and mutation types of cancer genes.

BACKGROUND: Alu elements are the most abundant retrotransposable elements comprising ~11% of the human genome. Many studies have highlighted the role that Alu elements have in genetic instability and how their contribution to the assortment of mutagenic events can lead to cancer. As of yet, little has been done to quantitatively assess the association between Alu distribution and genes that are causally implicated in oncogenesis. RESULTS: We have investigated the effect of various Alu densities on the mutation type based classifications of cancer genes. In order to establish the direct relationship between Alus and the cancer genes of interest, genome wide Alu-related densities were measured using genes rather than the sliding windows of fixed length as the units. Several novel genomic features, such as the density of the adjacent Alu pairs and the number of Alu-Exon-Alu triplets, were developed in order to extend the investigation via the multivariate statistical analysis toward more advanced biological insight. In addition, we characterized the genome-wide intron Alu distribution with a mixture model that distinguished genes containing Alu elements from those with no Alus, and evaluated the gene-level effect of the 5'-TTAAAA motif associated with Alu insertion sites using a two-step regression analysis method. CONCLUSIONS: The study resulted in several novel findings worthy of further investigation. They include: (1) Recessive cancer genes (tumor suppressor genes) are enriched with Alu elements (p < 0.01) compared to dominant cancer genes (oncogenes) and the entire set of genes in the human genome; (2) Alu-related genomic features can be used to cluster cancer genes into biological meaningful groups; (3) The retention of exon Alus has been restricted in the human genome development, and an upper limit to the chromosome-level exon Alu densities is suggested by the distribution profile; (4) For the genes with at least one intron Alu repeat in individual chromosomes, the intron Alu densities can be well fitted by a Gamma distribution; (5) The effect of the 5'-TTAAAA motif on Alu densities varies across different chromosomes.

Near-IR single fluorophore quenching system based on phthalocyanine (Pc) aggregation and its application for monitoring inhibitor/activator action on a therapeutic target: L1-EN.

Controlled H-aggregation of single Pc-labeled oligonucleotides is utilized as a fluorescence quenching system to discern changes in enzyme activity for the discovery of inhibitors for Long Interspersed Element 1 endonuclease (L1-EN), which is involved in genome instability and implicated in many different diseases.

Cross-talk-free dual-color fluorescence cross-correlation spectroscopy for the study of enzyme activity.

We have developed an instrument for spectral cross-talk-free dual-color fluorescence cross-correlation spectroscopy (FCCS), which provides a readout modality for the study of enzyme activity in application areas such as high-throughput screening. Two spectrally distinct (approximately 250 nm) fluorophores, Cy3 and IRD800, were excited simultaneously using two different excitation sources: one poised at 532 nm and the other at 780 nm. The fluorescence information was processed on two different color channels monitored with single-photon avalanche diodes (SPADs) that could transduce events at the single-molecule level. The system provided no color cross-talk (cross-excitation and/or cross-emission) and/or fluorescence resonance energy transfer (FRET), significantly improving data quality. To provide evidence of cross-talk-free operation, the system was evaluated using bright microspheres (lambda(abs) = 532 nm, lambda(em) = 560 nm) and quantum dots (lambda(abs) = 532 nm, lambda(em) = 810 nm). Experimental results indicated that no color leakage from the microspheres or quantum dots into inappropriate color channels was observed. To demonstrate the utility of the system, the enzymatic activity of APE1, which is responsible for nicking the phosphodiester backbone in DNA on the 5' side of an apurinic/apyrimidinic site, was monitored by FCCS using a double-stranded DNA substrate dual labeled with Cy3 and IRD800. Activity of APE1 was also monitored in the presence of an inhibitor (7-nitroindole-2-carboxylic acid) of the enzyme using this cross-talk-free FCCS platform. In all cases, no spectral leakage from single-molecule events into inappropriate color channels was observed.

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.

Rare mutations of FGFR2 causing apert syndrome: identification of the first partial gene deletion, and an Alu element insertion from a new subfamily.

Apert syndrome (AS) is a severe disorder, characterized by craniosynostosis and complex syndactyly of the hands and feet. Two heterozygous gain-of-function substitutions (Ser252Trp and Pro253Arg) in exon IIIa of fibroblast growth factor receptor 2 (FGFR2) are responsible for >98% of cases. Here we describe two novel mutations in FGFR2 in the two patients in whom a mutation had not previously been found in our cohort of 227 AS cases. The first is a 1.93-kb deletion, removing exon IIIc and substantial portions of the flanking introns. This is the first large FGFR2 deletion described in any individual with craniosynostosis. The other mutation is a 5' truncated Alu insertion into exon IIIc. This is the third Alu insertion identified in AS; all have occurred within an interval of only 104 bp, representing an enrichment of over a million-fold compared to the background genomic rate. We show that the inserted Alu element belongs to a small subfamily, not previously known to be mobile, which we term Alu Yk13. Both the deletion and insertion are likely to act by a similar gain-of-function mechanism in which disruption of exon IIIc leads to illegitimate mesenchymal expression of an FGFR2 spliceform containing the alternatively spliced exon IIIb. All the AS-associated Alu insertions have arisen in the paternal germline; we propose that their enrichment in FGFR2 is driven by positive selection of the mutant spermatogonial progenitors, a mechanism analogous to that explaining why the canonical AS nucleotide substitutions also reach exceptionally high levels in sperm.

ERCC1/XPF limits L1 retrotransposition.

Retrotransposons are currently active in the human and mouse genomes contributing to novel disease mutations and genomic variation via de novo insertions. However, little is known about the interactions of non-long terminal repeat (non-LTR) retrotransposons with the host DNA repair machinery. Based on the model of retrotransposition for the human and mouse LINE-1 element, one likely intermediate is an extension of cDNA that is heterologous to the genomic target, a flap intermediate. To determine whether a human flap endonuclease could recognize and process this potential intermediate, the genetic requirement for the ERCC1/XPF heterodimer during LINE-1 retrotransposition was characterized. Reduction of XPF in human cells increased retrotransposition whereas complementation of ERCC1-deficiency in hamster cells reduced retrotransposition. These results demonstrate for the first time that DNA repair enzymes act to limit non-LTR retrotransposition and may provide insight into the genetic instability phenotypes of ercc1 and xpf individuals.

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.