A tissue-specific decrease in the pre-mRNA level of L1- and alu-containing alleles of human genes

LINE1 and Alu retroelements occupy approximately 17 and 13% of the human genome, respectively. They include the evolutionarily youngest element groups Ta-L1, AluYa5, and AluYb8, many inserts of which are polymorphous in the Homo sapiens population. Despite the data on the ability of L1 and Alu elements to cause various modifications of the genome, the effects of these retroelements on gene expression have yet not been studied. Using the RT PCR method, we analyzed the pre-mRNA (heterogeneous nuclear RNA) content of allele pairs of four genes in five human cell lines, heterozygous with respect to intronic inserts of L1 and Alu elements. We showed for the first time a tissue-specific decrease in the pre-mRNA content of the gene allele bearing L1 or Alu inserts relative to the other allele of the same gene lacking the retroelement.     

Retroposons in modern human genome evolution

The ascertainment of the rates and driving forces of human genome evolution along with the genetic diversity of populations or separate population groups remains a topical problem of fundamental and applied genomics. According to the results of comparative analysis, the most numerous human genome structure peculiarities are connected with the distribution of mobile genetic retroelements—LTR, LINE1, SVA, and Alu repeats. Due to the wide distribution in different genome loci, conversed retropositional activity, and the retroelements’ regulatory potential, let us regard them as one of the significant evolutionary driving forces and the source of human genome variability. In the current review, we summarize published data and recent results of our research aimed at the analysis of the evolutionary impact of the young retroelements group on the function and variability of the human genome. We examine modern approaches of the polygenomic identification of polymorphic retroelements inserts. Using an original Internet resource, we analyze special features of the genomic polymorphic inserts of AluY repeats. We thoroughly characterize the strategy of large-scale functional analysis of polymorphic retroelement inserts. The presented results confirm the hypothesis of the roles of retroelements as active cis regulatory elements that are able to modulate surrounding genes.     

Most recent AluY insertions in human gene introns reduce the content of the primary transcripts in a cell type specific manner

Being the most effectively transposed primate-specific SINEs, Alu elements are present in more than one million copies in the human genome and include most recently transposed subsets of AluY elements that are polymorphic in humans. Although Alu elements are commonly thought to play an essential role in shaping and functioning of primate genomes, the understanding of the impact of recent Alu insertions on human gene expression is far from being comprehensive. Here we compared hnRNA contents for allele pairs of genes heterozygous for AluY insertions in their introns in human cell lines of various origins. We demonstrated that some AluY insertions correlated with decreased content of the corresponding hnRNAs. The effect observed does not depend on sequences of Alu elements and their orientation but is likely to be cell type specific.     

Analysis of the human Alu Ya-lineage

The Alu Ya-lineage is a group of related, short interspersed elements (SINEs) found in primates. This lineage includes subfamilies Ya1-Ya5, Ya5a2 and others. Some of these subfamilies are still actively mobilizing in the human genome. We have analyzed 2482 elements that reside in the human genome draft sequence and focused our analyses on the 2318 human autosomal Ya Alu elements. A total of 1470 autosomal loci were subjected to polymerase chain reaction (PCR)-based assays that allow analysis of individual Ya-lineage Alu elements. About 22% (313/1452) of the Ya-lineage Alu elements were polymorphic for the insertion presence on human autosomes. Less than 0.01% (5/1452) of the Ya-lineage loci analyzed displayed insertions in orthologous loci in non-human primate genomes. DNA sequence analysis of the orthologous inserts showed that the orthologous loci contained older pre-existing Y, Sc or Sq Alu subfamily elements that were the result of parallel forward insertions or involved in gene conversion events in the human lineage. This study is the largest analysis of a group of "young", evolutionarily related human subfamilies. The size, evolutionary age and variable allele insertion frequencies of several of these subfamilies makes members of the Ya-lineage useful tools for human population studies and primate phylogenetics.     

Selection against LINE-1 retrotransposons results principally from their ability to mediate ectopic recombination

LINE-1 (L1) retrotransposons constitute the most successful family of autonomous retroelements in mammals and they represent at least 17% of the size of the human genome. L1 insertions have occasionally been recruited to perform a beneficial function but the vast majority of L1 inserts are either neutral or deleterious. The basis for the deleterious effect of L1 remains a matter of debate and three possible mechanisms have been suggested: the direct effect of L1 inserts on gene activity, genetic rearrangements caused by L1-mediated ectopic recombination, or the retrotransposition process per se. We performed a genome-wide analysis of the distribution of L1 retrotransposons relative to the local recombination rate and the age and length of the elements. The proportion of L1 elements that are longer than 1.2 Kb is higher in low-recombining regions of the genome than in regions with a high recombination rate, but the genomic distributions of full-length elements (i.e. elements capable of retrotransposition) and long truncated elements were indistinguishable. We also found that the intensity of selection against long elements is proportional to the replicative success of L1 families. This suggests that the deleterious effect of L1 elements results principally from their ability to mediate ectopic recombination.     

The Distribution of L1 and Alu Retroelements in Relation to GC Content on Human Sex Chromosomes Is Consistent with the Ectopic Recombination Model

The distribution of Alu and L1 retroelements in the human genome changes with their age. Active retroelements target AT-rich regions, but their frequency increases in GC- and gene-rich regions of the genome with increasing age of the insertions. Currently there is no consensus on the mechanism generating this pattern. In this paper we test the hypothesis that selection against deleterious deletions caused by ectopic recombination between repeats is the main cause of the inhomogeneous distribution of L1s and Alus, by means of a detailed analysis of the GC distribution of the repeats on the sex chromosomes. We show that (1) unlike on the autosomes and X chromosome, L1s do not accumulate on the Y chromosome in GC-rich regions, whereas Alus accumulate there to a minor extent; (2) on the Y chromosome Alu and L1 densities are positively correlated, unlike the negative correlation on other chromosomes; and (3) in gene-poor regions of chromosome 4 and X, the distribution of Alus and L1s does not shift toward GC-rich regions. In addition, we show that although local GC content of long L1 insertions is lower than average, their selective loss from recombining chromosomes is not the main cause of the enrichment of ancient L1s in GC-rich regions. The results support the hypothesis that ectopic recombination causes the shift of Alu and L1 distributions toward the gene-rich regions of the genome. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. Reviewing Editor: Dr. Deborah Charlesworth     

The length of CpG islands is associated with the distribution of Alu and L1 retroelements

Alu and L1 retroelements have been suggested to initiate the spread of CpG methylation. In this study, the spread of CpG methylation was estimated based on the distance between the CpG islands and the nearest retroelements. All human genes (23,116) were examined and the correlations between the length of the CpG islands and the distance and density of the confronting retroelements were examined using nonoverlapping 5-kb windows. There was a linear relationship between the length of the CpG islands and the density of the Alu elements and an inverse relationship between the CpG islands and the L1 elements located more distantly, suggesting a suppressive effect of the Alu's on the spread of L1 methylation. Methylation analysis of the transitional CpG sites between the CpG islands and the nearest retroelements upstream of 16 genes was then carried out using DNA preparations from 11 different human tissues. Methylation-variable transitional CpGs were observed for the selected genes and the different tissues.     

Characterization of pre-insertion loci of de novo L1 insertions

The human Long Interspersed Element-1 (LINE-1) and the Short Interspersed Element (SINE) Alu comprise 28% of the human genome. They share the same L1-encoded endonuclease for insertion, which recognizes an A+T-rich sequence. Under a simple model of insertion distribution, this nucleotide preference would lead to the prediction that the populations of both elements would be biased towards A+T-rich regions. Genomic L1 elements do show an A+T-rich bias. In contrast, Alu is biased towards G+C-rich regions when compared to the genome average. Several analyses have demonstrated that relatively recent insertions of both elements show less G+C content bias relative to older elements. We have analyzed the repetitive element and G+C composition of more than 100 pre-insertion loci derived from de novo L1 insertions in cultured human cancer cells, which should represent an evolutionarily unbiased set of insertions. An A+T-rich bias is observed in the 50 bp flanking the endonuclease target site, consistent with the known target site for the L1 endonuclease. The L1, Alu, and G+C content of 20 kb of the de novo pre-insertion loci shows a different set of biases than that observed for fixed L1s in the human genome. In contrast to the insertion sites of genomic L1s, the de novo L1 pre-insertion loci are relatively L1-poor, Alu-rich and G+C neutral. Finally, a statistically significant cluster of de novo L1 insertions was localized in the vicinity of the c-myc gene. These results suggest that the initial insertion preference of L1, while A+T-rich in the initial vicinity of the break site, can be influenced by the broader content of the flanking genomic region and have implications for understanding the dynamics of L1 and Alu distributions in the human genome.     

Characterization of pre-insertion loci of de novo L1 insertions

The human Long Interspersed Element-1 (LINE-1) and the Short Interspersed Element (SINE) Alu comprise 28% of the human genome. They share the same L1-encoded endonuclease for insertion, which recognizes an A+T-rich sequence. Under a simple model of insertion distribution, this nucleotide preference would lead to the prediction that the populations of both elements would be biased towards A+T-rich regions. Genomic L1 elements do show an A+T-rich bias. In contrast, Alu is biased towards G+C-rich regions when compared to the genome average. Several analyses have demonstrated that relatively recent insertions of both elements show less G+C content bias relative to older elements. We have analyzed the repetitive element and G+C composition of more than 100 pre-insertion loci derived from de novo L1 insertions in cultured human cancer cells, which should represent an evolutionarily unbiased set of insertions. An A+T-rich bias is observed in the 50 bp flanking the endonuclease target site, consistent with the known target site for the L1 endonuclease. The L1, Alu, and G+C content of 20 kb of the de novo pre-insertion loci shows a different set of biases than that observed for fixed L1s in the human genome. In contrast to the insertion sites of genomic L1s, the de novo L1 pre-insertion loci are relatively L1-poor, Alu-rich and G+C neutral. Finally, a statistically significant cluster of de novo L1 insertions was localized in the vicinity of the c-myc gene. These results suggest that the initial insertion preference of L1, while A+T-rich in the initial vicinity of the break site, can be influenced by the broader content of the flanking genomic region and have implications for understanding the dynamics of L1 and Alu distributions in the human genome.     

Spreading of Alu methylation to the promoter of the MLH1 gene in gastrointestinal cancer

The highly repetitive Alu retroelements are regarded as methylation centres in the genome. Methylation in the gene promoters could be spreading from them. Promoter methylation of MLH1 is frequently detected in cancers, but the underlying mechanism is unclear. The aim of this study is to understand whether the methylation in the Alu elements is associated with promoter methylation in the MLH1 gene. Bisulfite genomic sequencing was used to analyse the CpG sites of the 5' end (promoter, exon 1 and Alu-containing intron 1) of the MLH1 gene in colorectal cancer cells and tissues, and gastric cancer tissues. Hypomethylation in the Alu elements and hypermethylation in the promoters and the regions between the promoters and the Alu elements were detected in two cancer cell lines and seven cancer tissues. However, demethylation or hypomethylation of the MLH1 promoter and regions between promoter and the Alu elements, and hypermethylation in the Alu elements, were identified in the normal tissues. MLH1 promoter methylation may spread from Alu elements that are located in intron 1 of the MLH1 gene. The trans-acting elements binding to the mutation sites could play a role in the methylation spreading.     

Active Alu elements are passed primarily through paternal germlines

Repetitive elements are distributed non-randomly in the human genome but, as reviewed in this paper, biological processes underlying the observed patterns appear to be complex and remain relatively obscure. Recent findings indicate that chromosomal distribution of Alu retroelements deposited in the past is different from the distribution of Alu elements that continue to be inserted in human population. These active elements from AluY sub(sub)families are the major focus of this paper. In particular, we analyzed chromosomal proportions of 19 AluY subfamilies, of which nine are reported for the first time in this paper. These 19 subfamilies contain over 80% of Alu elements that are polymorphic in the human genome. The chromosomal density of these most recent Alu insertions is around three times higher on chromosome Y than on chromosome X and over two times higher than the average density for all human autosomes. Based on this observation and other data we propose that active Alu elements are passed through paternal germlines. There is also some evidence that a small fraction of active Alu elements from less abundant subfamilies can be retroposed in female germlines or in the early embryos. Finally, we propose that the origin of Alu subfamilies in human populations may be related to evolution of chromosome Y.     

Footprint of APOBEC3 on the Genome of Human Retroelements

Almost half of the human genome is composed of transposable elements. The genomic structures and life cycles of some of these elements suggest they are a result of waves of retroviral infection and transposition over millions of years. The reduction of retrotransposition activity in primates compared to that in nonprimates, such as mice, has been attributed to the positive selection of several antiretroviral factors, such as apolipoprotein B mRNA editing enzymes. Among these, APOBEC3G is known to mutate G to A within the context of GG in the genome of endogenous as well as several exogenous retroelements (the underlining marks the G that is mutated). On the other hand, APOBEC3F and to a lesser extent other APOBEC3 members induce G-to-A changes within the nucleotide GA. It is known that these enzymes can induce deleterious mutations in the genome of retroviral sequences, but the evolution and/or inactivation of retroelements as a result of mutation by these proteins is not clear. Here, we analyze the mutation signatures of these proteins on large populations of long interspersed nuclear element (LINE), short interspersed nuclear element (SINE), and endogenous retrovirus (ERV) families in the human genome to infer possible evolutionary pressure and/or hypermutation events. Sequence context dependency of mutation by APOBEC3 allows investigation of the changes in the genome of retroelements by inspecting the depletion of G and enrichment of A within the APOBEC3 target and product motifs, respectively. Analysis of approximately 22,000 LINE-1 (L1), 24,000 SINE Alu, and 3,000 ERV sequences showed a footprint of GG→AG mutation by APOBEC3G and GA→AA mutation by other members of the APOBEC3 family (e.g., APOBEC3F) on the genome of ERV-K and ERV-1 elements but not on those of ERV-L, LINE, or SINE.     

Potential for genomic instability associated with retrotranspositionally-incompetent L1 loci

Expression of the L1 retrotransposon can damage the genome through insertional mutagenesis and the generation of DNA double-strand breaks (DSBs). The majority of L1 loci in the human genome are 5′-truncated and therefore incapable of retrotransposition. While thousands of full-length L1 loci remain, most are retrotranspositionally-incompetent due to inactivating mutations. However, mutations leading to premature stop codons within the L1 ORF2 sequence may yield truncated proteins that retain a functional endonuclease domain. We demonstrate that some truncated ORF2 proteins cause varying levels of toxicity and DNA damage when chronically overexpressed in mammalian cells. Furthermore, transfection of some ORF2 constructs containing premature stop codons supported low levels of Alu retrotransposition, demonstrating the potential for select retrotranspositionally-incompetent L1 loci to generate genomic instability. This result suggests yet another plausible explanation for the relative success of Alu elements in populating the human genome. Our data suggest that a subset of retrotranspositionally-incompetent L1s, previously considered to be harmless to genomic integrity, may have the potential to cause chronic DNA damage by introducing DSBs and mobilizing Alu. These results imply that the number of known L1 loci in the human genome that potentially threaten its stability may not be limited to the retrotranspositionally active loci.     

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.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1374-y) contains supplementary material, which is available to authorized users.