LINE-1 endonuclease-dependent retrotranspositional events causing human genetic disease: mutation detection bias and multiple mechanisms of target gene disruption

SINEs as Potential Expression Cassettes: Impact of Deletions and Insertions on Polyadenylation and Lifetime of B2 and Ves SINE Transcripts Generated by RNA Polymerase III

Short Interspersed Elements (SINEs) are common in the genomes of most multicellular organisms. They are transcribed by RNA polymerase III from an internal promoter comprising boxes A and B. As transcripts of certain SINEs from mammalian genomes can be polyadenylated, such transcripts should contain the AATAAA sequence as well as those called β- and τ-signals. One of the goals of this work was to evaluate how autonomous and independent other SINE parts are β- and τ-signals. Extended regions outside of β- and τ-signals were deleted from SINEs B2 and Ves and the derived constructs were used to transfect HeLa cells in order to evaluate the relative levels of their transcripts as well as their polyadenylation efficiency. If the deleted regions affected boxes A and B, the 5'-flanking region of the U6 RNA gene with the external promoter was inserted upstream. Such substitution of the internal promoter in B2 completely restored its transcription. Almost all tested deletions/substitutions did not reduce the polyadenylation capacity of the transcripts, indicating a weak dependence of the function of β- and τ-signals on the neighboring sequences. A similar analysis of B2 and Ves constructs containing a 55-bp foreign sequence inserted between β- and τ-signals showed an equal polyadenylation efficiency of their transcripts compared to those of constructs without the insertion. The acquired poly(A)-tails significantly increased the lifetime and thus the cellular level of such transcripts. The data obtained highlight the potential of B2 and Ves SINEs as cassettes for the expression of relatively short sequences for various applications.

Tail Wags Dog’s SINE: Retropositional Mechanisms of Can SINE Depend on Its A-Tail Structure

Simple Summary

The genomes of higher organisms including humans are invaded by millions of repetitive elements (transposons), which can sometimes be deleterious or beneficial for hosts. Many aspects of the mechanisms underlying the expansion of transposons in the genomes remain unclear. Short retrotransposons (SINEs) are one of the most abundant classes of genomic repeats. Their amplification relies on two major processes: transcription and reverse transcription. Here, short retrotransposons of dogs and other canids called Can SINE were analyzed. Their amplification was extraordinarily active in the wolf and, particularly, dog breeds relative to other canids. We also studied a variation of their transcription mechanism involving the polyadenylation of transcripts. An analysis of specific signals involved in this process allowed us to conclude that Can SINEs could alternate amplification with and without polyadenylation in their evolution. Understanding the mechanisms of transposon replication can shed light on the mechanisms of genome function.

Abstract

SINEs, non-autonomous short retrotransposons, are widespread in mammalian genomes. Their transcripts are generated by RNA polymerase III (pol III). Transcripts of certain SINEs can be polyadenylated, which requires polyadenylation and pol III termination signals in their sequences. Our sequence analysis divided Can SINEs in canids into four subfamilies, older a1 and a2 and younger b1 and b2. Can_b2 and to a lesser extent Can_b1 remained retrotranspositionally active, while the amplification of Can_a1 and Can_a2 ceased long ago. An extraordinarily high Can amplification was revealed in different dog breeds. Functional polyadenylation signals were analyzed in Can subfamilies, particularly in fractions of recently amplified, i.e., active copies. The transcription of various Can constructs transfected into HeLa cells proposed AATAAA and (TC)_n as functional polyadenylation signals. Our analysis indicates that older Can subfamilies (a1, a2, and b1) with an active transcription terminator were amplified by the T⁺ mechanism (with polyadenylation of pol III transcripts). In the currently active Can_b2 subfamily, the amplification mechanisms with (T⁺) and without the polyadenylation of pol III transcripts (T⁻) irregularly alternate. The active transcription terminator tends to shorten, which renders it nonfunctional and favors a switch to the T⁻ retrotransposition. The activity of a truncated terminator is occasionally restored by its elongation, which rehabilitates the T⁺ retrotransposition for a particular SINE copy.

Transcriptomic profiles of muscular dystrophy with myositis (mdm) in extensor digitorum longus, psoas, and soleus muscles from mice

Background

Titinopathies are inherited muscular diseases triggered by genetic mutations in the titin gene. Muscular dystrophy with myositis (mdm) is one such disease caused by a LINE repeat insertion, leading to exon skipping and an 83-amino acid residue deletion in the N2A-PEVK region of mouse titin. This region has been implicated in a number of titin—titin ligand interactions, hence are important for myocyte signaling and health. Mice with this mdm mutation develop a severe and progressive muscle degeneration. The range of phenotypic differences observed in mdm mice shows that the deletion of this region induces a cascade of transcriptional changes extending to numerous signaling pathways affected by the titin filament. Previous research has focused on correlating phenotypic differences with muscle function in mdm mice. These studies have provided understanding of the downstream physiological effects resulting from the mdm mutation but only provide insights on processes that can be physiologically observed and measured. We used differential gene expression (DGE) to compare the transcriptomes of extensor digitorum longus (EDL), psoas and soleus muscles from wild-type and mdm mice to develop a deeper understand of these tissue-specific responses.

Results

The overall expression pattern observed shows a well-differentiated transcriptional signature in mdm muscles compared to wild type. Muscle-specific clusters observed within the mdm transcriptome highlight the level of variability of each muscle to the deletion. Differential gene expression and weighted gene co-expression network analysis showed a strong directional response in oxidative respiration-associated mitochondrial genes, which aligns with the poor shivering and non-shivering thermogenesis previously observed. Sln, which is a marker associated with shivering and non-shivering thermogenesis, showed the strongest expression change in fast-fibered muscles. No drastic changes in MYH expression levels were reported, which indicated an absence of major fiber-type switching events. Overall expression shifts in MYH isoforms, MARPs, and extracellular matrix associated genes demonstrated the transcriptional complexity associated with mdm mutation. The expression alterations in mitochondrial respiration and metabolism related genes in the mdm muscle dominated over other transcriptomic changes, and likely account for the late stage cellular responses in the mdm muscles.

Conclusions

We were able to demonstrate that the complex nature of mdm mutation extends beyond a simple rearrangement in titin gene. EDL, psoas and soleus exemplify unique response modes observed in skeletal muscles with mdm mutation. Our data also raises the possibility that failure to maintain proper energy homeostasis in mdm muscles may contribute to the pathogenesis of the degenerative phenotype in mdm mice. Understanding the full disease-causing molecular cascade is difficult using bulk RNA sequencing techniques due to intricate nature of the disease. The development of the mdm phenotype is temporally and spatially regulated, hence future studies should focus on single fiber level investigations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08873-2.

Intragenic L1 Insertion: One Possibility of Brain Disorder

Long interspersed nuclear element 1 (LINE1, L1) is a retrotransposon comprising ~17% of the human genome. A subset of L1s maintains the potential to mobilize and alter the genomic landscape, consequently contributing to the change in genome integrity and gene expression. L1 retrotransposition occurs in the human brain regardless of disease status. However, in the brain of patients with various brain diseases, the expression level and copy number of L1 are significantly increased. In this review, we briefly introduce the methodologies applied to measure L1 mobility and identify genomic loci where new insertion of L1 occurs in the brain. Then, we present a list of genes disrupted by L1 transposition in the genome of patients with brain disorders. Finally, we discuss the association between genes disrupted by L1 and relative brain disorders.

Insertion of LINE-1 Retrotransposon Inducing Exon Inversion Causes a Rotor Syndrome Phenotype

Rotor syndrome, a rare autosomal-recessive genetic disorder characterized by conjugated hyperbilirubinemia, is caused by biallelic pathogenic variants in both SLCO1B1 and SLCO1B3 genes. Long interspersed nuclear elements (LINEs) make up about 17% of the human genome and insertion of LINE-1 in genes can result in genetic diseases. In the current study, we examined SLCO1B1 and SLCO1B3 genes in two Chinese patients diagnosed with Rotor syndrome based on laboratory tests. In one patient, a novel exon 4 inversion variant was identified. This variant may have been induced by LINE-1 retrotransposon insertion into SLCO1B3 intron 3, and was identified using genome walking. Splicing assay results indicated that the exon inversion, resulting in SLCO1B3 exon 4 (122 bp) exclusion in the mature mRNA, might generate a premature termination codon. Here, we describe an exon inversion contributing to the molecular etiology of Rotor syndrome. Our results may inform future diagnoses and guide drug prescriptions and genetic counseling.

DNA methylation of individual repetitive elements in hepatitis C virus infection-induced hepatocellular carcinoma

Background

The two most common repetitive elements (REs) in humans, long interspersed nuclear element-1 (LINE-1) and Alu element (Alu), have been linked to various cancers. Hepatitis C virus (HCV) may cause hepatocellular carcinoma (HCC) by suppressing host defenses, through DNA methylation that controls the mobilization of REs. We aimed to investigate the role of RE methylation in HCV-induced HCC (HCV-HCC).

Results

We studied methylation of over 30,000 locus-specific REs across the genome in HCC, cirrhotic, and healthy liver tissues obtained by surgical resection. Relative to normal liver tissue, we observed the largest number of differentially methylated REs in HCV-HCC followed by alcohol-induced HCC (EtOH-HCC). After excluding EtOH-HCC-associated RE methylation (FDR < 0.001) and those unable to be validated in The Cancer Genome Atlas (TCGA), we identified 13 hypomethylated REs (11 LINE-1 and 2 Alu) and 2 hypermethylated REs (1 LINE-1 and 1 Alu) in HCV-HCC (FDR < 0.001). A majority of these REs were located in non-coding regions, preferentially enriched with chromatin repressive marks H3K27me3, and positively associated with gene expression (median correlation r = 0.32 across REs). We further constructed an HCV-HCC RE methylation score that distinguished HCV-HCC (lowest score), HCV-cirrhosis, and normal liver (highest score) in a dose-responsive manner (p for trend < 0.001). HCV-cirrhosis had a lower score than EtOH-cirrhosis (p = 0.038) and HCV-HCC had a lower score than EtOH-HCC in TCGA (p = 0.024).

Conclusions

Our findings indicate that HCV infection is associated with loss of DNA methylation in specific REs, which could implicate molecular mechanisms in liver cancer development. If our findings are validated in larger sample sizes, methylation of these REs may be useful as an early detection biomarker for HCV-HCC and/or a target for prevention of HCC in HCV-positive individuals.

Electronic supplementary material

The online version of this article (10.1186/s13148-019-0733-y) contains supplementary material, which is available to authorized users.

Structural Characterization of ABCB1, the Gene Underlying the d2 Dwarf Phenotype in Pearl Millet, Cenchrus Americanus (L.) Morrone

Pearl millet is an important food crop in arid and semi-arid regions of South Asia and sub-Saharan Africa and is grown in Australia and the United States as a summer fodder crop. The d2 dwarf germplasm has been widely used in the last half-century to develop high-performing pearl millet hybrids. We previously mapped the d2 phenotype to a 1.6 cM region in linkage group (LG) 4 and identified the ABCB1 gene as a candidate underlying the trait. Here, we report the sequence, structure and expression of ABCB1 in tall (D2D2) and d2 dwarf (d2d2) germplasm. The ABCB1 allele in d2 dwarfs differs from that in tall inbreds by the presence of two different high copy transposable elements, one in the coding region and the second located 664 bp upstream of the ATG start codon. These transposons were present in all d2 dwarfs tested that were reported to be of independent origin and absent in the analyzed wild-type tall germplasm. We also compared the expression profile of this gene in different organs of multiple tall and d2 dwarf inbreds, including the near-isogenic inbreds at the d2 locus, Tift 23B (D2D2) and Tift 23DB (d2d2). Heterologous transformation of the tall (Ca_ABCB1) and the d2 dwarf (Ca_abcb1) pearl millet alleles in the Arabidopsis double mutant abcb1abcb19 showed that the pearl millet D2 but not the d2 allele complements the Arabidopsis abcb1 mutation. Our studies also show the importance of the COOH-terminal 22 amino acids of the ABCB1 protein in either protein function or stability.

A study of meiomitosis and novel pathways of genomic instability in cutaneous T-cell lymphomas (CTCL)

Genomic instability is a hallmark of cancer and an enabling factor for genetic alterations that drive cancer development and progression. The clashing of mitosis and aberrantly expressed meiosis machineries, which may contribute to genomic instability, has been coined cancer “meiomitosis”. LINE-1 retrotransposition, a process active in germ cells, acts outside of the meiotic machinery to create DNA double strand breaks (DNA DSBs) and has played an important role in the evolution of the human genome. We have previously demonstrated that in CTCL several cancer testis/meiotic genes are expressed. Furthermore, this cancer exhibits extensive and ongoing chromosomal/microsatellite instability. In this study we analyzed immortalized patient-derived cells and primary CTCL patient samples using RT-PCR, western blotting and confocal microscopy and found that proteins critically involved in meiosis and LINE-1 retrotransposition are expressed and are associated with chromosomal instability and DNA DSB formation. Using cell cycle synchronization, we show G1/S phase-transition-specific expression of meiosis proteins. Using the Alu retrotransposition assay, we demonstrate the functional activity of LINE-1 retrotransposon in CTCL. Histone acetyltransferase inhibition results in downregulation of the ectopic germ cell programs and concomitant decrease in DNA DSBs foci formation. Notably, LINE-1 and meiosis genes were expressed across a panel of other solid tumor cell lines. Taken together, our results indicate that malignant cells in culture undergo “cancer meiomitosis” rather than the classic mitosis division. The ectopic expression of meiosis genes and reactivation of LINE-1 may be contributing to genomic instability and represent novel targets for immunotherapy in this and other cancers.

Endogenous LINE-1 (Long Interspersed Nuclear Element-1) Reverse Transcriptase Activity in Platelets Controls Translational Events Through RNA-DNA Hybrids

OBJECTIVE

One source of endogenous reverse transcriptase (eRT) activity in nucleated cells is the LINE-1/L1 (long interspersed nuclear element-1), a non-LTR retrotransposon that is implicated in the regulation of gene expression. Nevertheless, the presence and function of eRT activity and LINE-1 in human platelets, an anucleate cell, has not previously been determined.

APPROACH AND RESULTS

We demonstrate that human and murine platelets possess robust eRT activity and identify the source as being LINE-1 ribonucleoprotein particles. Inhibition of eRT in vitro in isolated platelets from healthy individuals or in people with HIV treated with RT inhibitors enhanced global protein synthesis and platelet activation. If HIV patients were treated with reverse transcriptase inhibitor, we found that platelets from these patients had increased basal activation. We next discovered that eRT activity in platelets controlled the generation of RNA-DNA hybrids, which serve as translational repressors. Inhibition of platelet eRT lifted this RNA-DNA hybrid-induced translational block and was sufficient to increase protein expression of target RNAs identified by RNA-DNA hybrid immunoprecipitation.

CONCLUSIONS

Thus, we provide the first evidence that platelets possess L1-encoded eRT activity. We also demonstrate that platelet eRT activity regulates platelet hyperreactivity and thrombosis and controls RNA-DNA hybrid formation and identify that RNA-DNA hybrids function as a novel translational control mechanism in human platelets.

Pivotal Impacts of Retrotransposon Based Invasive RNAs on Evolution

RNAs have long been described as the mediators of gene expression; they play a vital role in the structure and function of cellular complexes. Although the role of RNAs in the prokaryotes is mainly confined to these basic functions, the effects of these molecules in regulating the gene expression and enzymatic activities have been discovered in eukaryotes. Recently, a high-resolution analysis of the DNA obtained from different organisms has revealed a fundamental impact of the RNAs in shaping the genomes, heterochromatin formation, and gene creation. Deep sequencing of the human genome revealed that about half of our DNA is comprised of repetitive sequences (remnants of transposable element movements) expanded mostly through RNA-mediated processes. ORF2 encoded by L1 retrotransposons is a cellular reverse transcriptase which is mainly responsible for RNA invasion of various transposable elements (L1s, Alus, and SVAs) and cellular mRNAs in to the genomic DNA. In addition to increasing retroelements copy number; genomic expansion in association with centromere, telomere, and heterochromatin formation as well as pseudogene creation are the evolutionary consequences of this RNA-based activity. Threatening DNA integrity by disrupting the genes and forming excessive double strand breaks is another effect of this invasion. Therefore, repressive mechanisms have been evolved to control the activities of these invasive intracellular RNAs. All these mechanisms now have essential roles in the complex cellular functions. Therefore, it can be concluded that without direct action of RNA networks in shaping the genome and in the development of different cellular mechanisms, the evolution of higher eukaryotes would not be possible.

Identification of pathogenic retrotransposon insertions in cancer predisposition genes

Cancer risks have been previously reported for some retrotransposon element (RE) insertions; however, detection of these insertions is technically challenging and very few oncogenic RE insertions have been reported. Here we evaluate RE insertions identified during hereditary cancer genetic testing using a comprehensive testing strategy. Individuals who had single-syndrome or pan-cancer hereditary cancer genetic testing from February 2004 to March 2017 were included. RE insertions were identified using Sanger sequencing, Next Generation Sequencing, or multiplex quantitative PCR, and further characterized using targeted PCR and sequencing analysis. Personal cancer history, ancestry, and haplotype were evaluated. A total of 37 unique RE insertions were identified in 10 genes, affecting 211 individuals. BRCA2 accounted for 45.9% (17/37) of all unique RE insertions. Several RE insertions were detected with high frequency in populations of conserved ancestry wherein up to 100% of carriers shared a high degree of haplotype conservation, suggesting founder effects. Our comprehensive testing strategy resulted in a substantial increase in the number of reported oncogenic RE insertions, several of which may have possible founder effects. Collectively, these data show that the detection of RE insertions is an important component of hereditary cancer genetic testing and may be more prevalent than previously reported.

Recessive inheritance of population-specific intronic LINE-1 insertion causes a rotor syndrome phenotype

Sequences of long-interspersed elements (LINE-1, L1) make up ∼17% of the human genome. De novo insertions of retrotransposition-active L1s can result in genetic diseases. It has been recently shown that the homozygous inactivation of two adjacent genes SLCO1B1 and SLCO1B3 encoding organic anion transporting polypeptides OATP1B1 and OATP1B3 causes a benign recessive disease presenting with conjugated hyperbilirubinemia, Rotor syndrome. Here, we examined SLCO1B1 and SLCO1B3 genes in six Japanese diagnosed with Rotor syndrome on the basis of laboratory data and laparoscopy. All six Japanese patients were homozygous for the c.1738C>T nonsense mutation in SLCO1B1 and homozygous for the insertion of a ∼6.1-kbp L1 retrotransposon in intron 5 of SLCO1B3, which altogether make up a Japanese-specific haplotype. RNA analysis revealed that the L1 insertion induced deleterious splicing resulting in SLCO1B3 transcripts lacking exon 5 or exons 5-7 and containing premature stop codons. The expression of OATP1B1 and OATP1B3 proteins was not detected in liver tissues. This is the first documented case of a population-specific polymorphic intronic L1 transposon insertion contributing to molecular etiology of recessive genetic disease. Since L1 activity in human genomes is currently seen as a major source of individual genetic variation, further investigations are warranted to determine whether this phenomenon results in other autosomal-recessive diseases.

The impact of transposable elements in genome evolution and genetic instability and their implications in various diseases

Approximately 45% of the human genome is comprised of transposable elements (TEs). Results from the Human Genome Project have emphasized the biological importance of TEs. Many studies have revealed that TEs are not simply "junk" DNA, but rather, they play various roles in processes, including genome evolution, gene expression regulation, genetic instability, and cancer disposition. The effects of TE insertion in the genome varies from negligible to disease conditions. For the past two decades, many studies have shown that TEs are the causative factors of various genetic disorders and cancer. TEs are a subject of interest worldwide, not only in terms of their clinical aspects but also in basic research, such as evolutionary tracking. Although active TEs contribute to genetic instability and disease states, non-long terminal repeat transposons are well studied, and their roles in these processes have been confirmed. In this review, we will give an overview of the importance of TEs in studying genome evolution and genetic instability, and we suggest that further in-depth studies on the mechanisms related to these phenomena will be useful for both evolutionary tracking and clinical diagnostics.

Reviving the dead: history and reactivation of an extinct l1

Although L1 sequences are present in the genomes of all placental mammals and marsupials examined to date, their activity was lost in the megabat family, Pteropodidae, ∼24 million years ago. To examine the characteristics of L1s prior to their extinction, we analyzed the evolutionary history of L1s in the genome of a megabat, Pteropus vampyrus, and found a pattern of periodic L1 expansion and quiescence. In contrast to the well-characterized L1s in human and mouse, megabat genomes have accommodated two or more simultaneously active L1 families throughout their evolutionary history, and major peaks of L1 deposition into the genome always involved multiple families. We compared the consensus sequences of the two major megabat L1 families at the time of their extinction to consensus L1s of a variety of mammalian species. Megabat L1s are comparable to the other mammalian L1s in terms of adenosine content and conserved amino acids in the open reading frames (ORFs). However, the intergenic region (IGR) of the reconstructed element from the more active family is dramatically longer than the IGR of well-characterized human and mouse L1s. We synthesized the reconstructed element from this L1 family and tested the ability of its components to support retrotransposition in a tissue culture assay. Both ORFs are capable of supporting retrotransposition, while the IGR is inhibitory to retrotransposition, especially when combined with either of the reconstructed ORFs. We dissected the inhibitory effect of the IGR by testing truncated and shuffled versions and found that length is a key factor, but not the only one affecting inhibition of retrotransposition. Although the IGR is inhibitory to retrotransposition, this inhibition does not account for the extinction of L1s in megabats. Overall, the evolution of the L1 sequence or the quiescence of L1 is unlikely the reason of L1 extinction.

Primary immunodeficiency caused by an exonized retroposed gene copy inserted in the CYBB gene

Retrotransposon-mediated insertion of a long interspersed nuclear element (LINE)-1 or an Alu element into a human gene is a well-known pathogenic mechanism. We report a novel LINE-1-mediated insertion of a transcript from the TMF1 gene on chromosome 3 into the CYBB gene on the X-chromosome. In a Dutch male patient with chronic granulomatous disease, a 5.8-kb, incomplete and partly exonized TMF1 transcript was identified in intron 1 of CYBB, in opposite orientation to the host gene. The sequence of the insertion showed the hallmarks of a retrotransposition event, with an antisense poly(A) tail, target site duplication, and a consensus LINE-1 endonuclease cleavage site. This insertion induced aberrant CYBB mRNA splicing, with inclusion of an extra 117-bp exon between exons 1 and 2 of CYBB. This extra exon contained a premature stop codon. The retrotransposition took place in an early stage of fetal development in the mother of the patient, because she showed a somatic mosaicism for the mutation that was not present in the DNA of her parents. However, the mutated allele was not expressed in the patient's mother because the insertion was found only in the methylated fraction of her DNA.

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.

How do mammalian transposons induce genetic variation? A conceptual framework: the age, structure, allele frequency, and genome context of transposable elements may define their wide-ranging biological impacts

In this essay, we discuss new insights into the wide-ranging impacts of mammalian transposable elements (TE) on gene expression and function. Nearly half of each mammalian genome is comprised of these mobile, repetitive elements. While most TEs are ancient relics, certain classes can move from one chromosomal location to another even now. Indeed, striking recent data show that extensive transposition occurs not only in the germline over evolutionary time, but also in developing somatic tissues and particular human cancers. While occasional germline TE insertions may contribute to genetic variation, many other, similar TEs appear to have little or no impact on neighboring genes. However, the effects of somatic insertions on gene expression and function remain almost completely unknown. We present a conceptual framework to understand how the ages, allele frequencies, molecular structures, and especially the genomic context of mammalian TEs each can influence their various possible functional consequences.

Impact of retrotransposons in pluripotent stem cells

Retrotransposons, which constitute approximately 40% of the human genome, have the capacity to 'jump' across the genome. Their mobility contributes to oncogenesis, evolution, and genomic plasticity of the host genome. Induced pluripotent stem cells as well as embryonic stem cells are more susceptible than differentiated cells to genomic aberrations including insertion, deletion and duplication. Recent studies have revealed specific behaviors of retrotransposons in pluripotent cells. Here, we review recent progress in understanding retrotransposons and provide a perspective on the relationship between retrotransposons and genomic variation in pluripotent stem cells.

Alu insertion profiling: array-based methods to detect Alu insertions in the human genome

The analysis of the genetic variability associated to Alu sequences was hampered by the absence of genome-wide methodologies able to efficiently detect new polymorphisms/mutations among these repetitive elements. Here we describe two Alu insertion profiling (AIP) methods based on the hybridization of Alu-flanking genomic fragments on tiling microarrays. Protocols are designed to preferentially detect active Alu subfamilies. We tested AIP methods by analyzing chromosomes 1 and 6 in two genomic samples. In genomic regions covered by array-features, with a sensitivity of 2% (AIP1) -4% (AIP2) and 5% (AIP1) -8% (AIP2) for the old J and S Alu lineages respectively, we obtained a sensitivity of 67% (AIP1) -90% (AIP2) for the young Ya subfamily. Among the loci showing sample-to-sample differences, 5 (AIP1) -8 (AIP2) were associated to known Alu polymorphisms. Moreover, we were able to confirm by PCR and DNA sequencing 4 new intragenic Alu elements, polymorphic in 10 additional individuals.

Mitochondrial pseudogenes in the nuclear genomes of Drosophila

Mitochondrial pseudogenes in nuclear chromosomes (numts) have been detected in the genomes of a diverse range of eukaryotic species. However, the numt content of different genomes and their properties is not uniform, and study of these differences provides insight into the mechanisms and dynamics of genome evolution in different organisms. In the genus Drosophila, numts have previously only been identified on a genome-wide scale in the melanogaster subgroup. The present study extends the identification to 11 species of the Drosophila genus. We identify a total of 302 numts and show that the numt complement is highly variable in Drosophilids, ranging from just 4 in D. melanogaster to 67 in D. willistoni, broadly correlating with genome size. Many numts have undergone large-scale rearrangements in the nucleus, including interruptions, inversions, deletions and duplications of sequence of variable size. Estimating the age of the numts in the nucleus by phylogenetic tree reconstruction reveals the vast majority of numts to be recent gains, 90% having arisen on terminal branches of the species tree. By identifying paralogs and counting duplications among the extant numts we estimate that 23% of extant numts arose through post-insertion duplications. We estimate genus average rates of insertion of 0.75 per million years, and a duplication rate of 0.010 duplications per numt per million years.

The NF1 gene contains hotspots for L1 endonuclease-dependent de novo insertion

Long interspersed (L1) and Alu elements are actively amplified in the human genome through retrotransposition of their RNA intermediates by the -100 still retrotranspositionally fully competent L1 elements. Retrotransposition can cause inherited disease if such an element is inserted near or within a functional gene. Using direct cDNA sequencing as the primary assay for comprehensive NF1 mutation analysis, we uncovered in 18 unrelated index patients splicing alterations not readily explained at the genomic level by an underlying point-mutation or deletion. Improved PCR protocols avoiding allelic drop-out of the mutant alleles uncovered insertions of fourteen Alu elements, three L1 elements, and one poly(T) stretch to cause these splicing defects. Taken together, the 18 pathogenic L1 endonuclease-mediated de novo insertions represent the largest number of this type of mutations characterized in a single human gene. Our findings show that retrotransposon insertions account for as many as -0.4% of all NF1 mutations. Since altered splicing was the main effect of the inserted elements, the current finding was facilitated by the use of RNA-based mutation analysis protocols, resulting in improved detection compared to gDNA-based approaches. Six different insertions clustered in a relatively small 1.5-kb region (NF1 exons 21(16)-23(18)) within the 280-kb NF1 gene. Furthermore, three different specific integration sites, one of them located in this cluster region, were each used twice, i.e. NM_000267.3(NF1):c.1642-1_1642 in intron 14(10c), NM_000267.3(NF1):c.2835_2836 in exon 21(16), and NM_000267.3(NF1):c.4319_4320 in exon 33(25). Identification of three loci that each served twice as integration site for independent retrotransposition events as well as 1.5-kb cluster region harboring six independent insertions supports the notion of non-random insertion of retrotransposons in the human genome. Currently, little is known about which features make sites particularly vulnerable to L1 EN-mediated insertions. The here identified integration sites may serve to elucidate these features in future studies.

A mobile threat to genome stability: The impact of non-LTR retrotransposons upon the human genome

It is now commonly agreed that the human genome is not the stable entity originally presumed. Deletions, duplications, inversions, and insertions are common, and contribute significantly to genomic structural variations (SVs). Their collective impact generates much of the inter-individual genomic diversity observed among humans. Not only do these variations change the structure of the genome; they may also have functional implications, e.g. altered gene expression. Some SVs have been identified as the cause of genetic disorders, including cancer predisposition. Cancer cells are notorious for their genomic instability, and often show genomic rearrangements at the microscopic and submicroscopic level to which transposable elements (TEs) contribute. Here, we review the role of TEs in genome instability, with particular focus on non-LTR retrotransposons. Currently, three non-LTR retrotransposon families - long interspersed element 1 (L1), SVA (short interspersed element (SINE-R), variable number of tandem repeats (VNTR), and Alu), and Alu (a SINE) elements - mobilize in the human genome, and cause genomic instability through both insertion- and post-insertion-based mutagenesis. Due to the abundance and high sequence identity of TEs, they frequently mislead the homologous recombination repair pathway into non-allelic homologous recombination, causing deletions, duplications, and inversions. While less comprehensively studied, non-LTR retrotransposon insertions and TE-mediated rearrangements are probably more common in cancer cells than in healthy tissue. This may be at least partially attributed to the commonly seen global hypomethylation as well as general epigenetic dysfunction of cancer cells. Where possible, we provide examples that impact cancer predisposition and/or development.

Mobile elements create structural variation: analysis of a complete human genome

Structural variants (SVs) are common in the human genome. Because approximately half of the human genome consists of repetitive, transposable DNA sequences, it is plausible that these elements play an important role in generating SVs in humans. Sequencing of the diploid genome of one individual human (HuRef) affords us the opportunity to assess, for the first time, the impact of mobile elements on SVs in an individual in a thorough and unbiased fashion. In this study, we systematically evaluated more than 8000 SVs to identify mobile element-associated SVs as small as 100 bp and specific to the HuRef genome. Combining computational and experimental analyses, we identified and validated 706 mobile element insertion events (including Alu, L1, SVA elements, and nonclassical insertions), which added more than 305 kb of new DNA sequence to the HuRef genome compared with the Human Genome Project (HGP) reference sequence (hg18). We also identified 140 mobile element-associated deletions, which removed approximately 126 kb of sequence from the HuRef genome. Overall, approximately 10% of the HuRef-specific indels larger than 100 bp are caused by mobile element-associated events. More than one-third of the insertion/deletion events occurred in genic regions, and new Alu insertions occurred in exons of three human genes. Based on the number of insertions and the estimated time to the most recent common ancestor of HuRef and the HGP reference genome, we estimated the Alu, L1, and SVA retrotransposition rates to be one in 21 births, 212 births, and 916 births, respectively. This study presents the first comprehensive analysis of mobile element-related structural variants in the complete DNA sequence of an individual and demonstrates that mobile elements play an important role in generating inter-individual structural variation.

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.

An Alu insertion in compound heterozygosity with a microduplication in GNPTAB gene underlies Mucolipidosis II

Mucolipidosis type II (ML II) is a fatal, autosomal recessive, lysosomal storage disorder characterized by severe clinical and radiologic features. ML II results from mutations in alpha and beta subunits, encoded by the GlcNAc-1-phosphotransferase gene (GNPTAB). Most of the 40 different GNPTAB mutations reported so far are insertions and deletions predicting diverse types of aberrant proteins. Alu mobile elements have however never been involved in these events up to now. The Italian ML II patient of our study showed an Alu retrotrasposition in GNPTAB exon 5. The Alu insertion mutation (NM_024312.3:c.555_556insHSU14569) generated a transcript with a skipping of the target exon 5 and a frameshift p.S122fs, causing a premature translation termination codon at position 123. This insertion mutation was found in compound heterozygosity with the frameshift p.S887KfsX33, resulting from a new mono-nucleotide duplication (c.2659dupA) that occurred in GNPTAB exon 13. A possible involvement of cis-splicing elements having an exonic location, such as exon enhancers (ESEs), is discussed as mechanism that led to the production of the aberrant mRNA splicing.

APOBEC3 proteins: major players in intracellular defence against LINE-1-mediated retrotransposition

Mammalian genomes are littered with enormous numbers of transposable elements interspersed within and between single-copy endogenous genes. The only presently spreading class of human transposable elements comprises non-LTR (long terminal repeat) retrotransposons, which cover approx. 34% of the human genome. Non-LTR retrotransposons include the widespread autonomous LINEs (long interspersed nuclear elements) and non-autonomous elements such as processed pseudogenes, SVAs [named after SINE (short interspersed nuclear element), VNTR (variable number of tandem repeats) and Alu] and SINEs. Mobilization of these elements affects the host genome, can be deleterious to the host cell, and cause genetic disorders and cancer. In order to limit negative effects of retrotransposition, host genomes have adopted several strategies to curb the proliferation of transposable elements. Recent studies have demonstrated that members of the human APOBEC3 (apolipoprotein B mRNA editing enzyme catalytic polypeptide 3) protein family inhibit the mobilization of the non-LTR retrotransposons LINE-1 and Alu significantly and participate in the intracellular defence against retrotransposition by mechanisms unknown to date. The striking coincidence between the expansion of the APOBEC3 gene cluster and the abrupt decline in retrotransposon activity in primates raises the possibility that these genes may have been expanded to prevent genomic instability caused by endogenous retroelements.

Progress in understanding the biology of the human mutagen LINE-1

Long interspersed nucleotide element (LINE)-1 retrotransposon (L1) has emerged as the largest contributor to mammalian genome mass, responsible for over 35% of the human genome. Differences in the number and activity levels of L1s contribute to interindividual variation in humans, both by affecting an individual's likelihood of acquiring new L1-mediated mutations, as well as by differentially modifying gene expression. Here, we report on recent progress in understanding L1 biology, with a focus on mechanisms of L1-mediated disease. We discuss known details of L1 life cycle, including L1 structure, transcriptional regulation, and the mechanisms of translation and retrotransposition. Current views on cell type specificity, timing, and control of retrotransposition are put forth. Finally, we discuss the role of L1 as a mutagen, using the latest findings in L1 biology to illuminate molecular mechanisms of L1-mediated gene disruption.

Detection of two Alu insertions in the CFTR gene

BACKGROUND

LINE-1 (long interspersed element-1) or L1-mediated retrotransposition is a potent force in human genome evolution and an occasional cause of human genetic disease. Since the first report of two de novo L1 insertions in the F8 gene causing hemophilia A, more than 50 L1-mediated retrotranspositional events have been identified as causing human genetic disease. However, a significant bias has generally militated against the detection of these pathological events at autosomal loci. Based upon this and other observations, we surmised that some previously unresolved cystic fibrosis chromosomes might carry hitherto undetected L1-mediated retrotranspositional insertions at the CFTR locus. This study represents an attempt to identify such mutational events.

METHODS

100 previously unresolved cystic fibrosis chromosomes were carefully reanalyzed using quantitative high-performance liquid chromatography (QHPLC).

RESULTS

Two simple Alu insertions were identified in the CFTR gene, within exons 16 and 17b respectively.

CONCLUSIONS

Our findings have not only revealed a previously unknown mutational mechanism responsible for cystic fibrosis but also represent an important addition to the already diverse spectrum of known CFTR gene mutations. Experience with the CFTR gene suggests that pathological L1-mediated retrotranspositional events may also have been overlooked at other gene loci and should always be considered in cases that appear to be refractory to analysis.

Mechanism of Alu integration into the human genome

LINE-1 or L1 has driven the generation of at least 10% of the human genome by mobilising Alu sequences. Although there is no doubt that Alu insertion is initiated by L1-dependent target site-primed reverse transcription, the mechanism by which the newly synthesised 3' end of a given Alu cDNA attaches to the target genomic DNA is less well understood. Intrigued by observations made on 28 pathological simple Alu insertions, we have sought to ascertain whether microhomologies could have played a role in the integration of shorter Alu sequences into the human genome. A meta-analysis of the 1624 Alu insertion polymorphisms deposited in the Database of Retrotransposon Insertion Polymorphisms in Humans (dbRIP), when considered together with a re-evaluation of the mechanism underlying how the three previously annotated large deletion-associated short pathological Alu inserts were generated, enabled us to present a unifying model for Alu insertion into the human genome. Since Alu elements are comparatively short, L1 RT is usually able to complete nascent Alu cDNA strand synthesis leading to the generation of full-length Alu inserts. However, the synthesis of the nascent Alu cDNA strand may be terminated prematurely if its 3' end anneals to the 3' terminal of the top strand's 5' overhang by means of microhomology-mediated mispairing, an event which would often lead to the formation of significantly truncated Alu inserts. Furthermore, the nascent Alu cDNA strand may be 'hijacked' to patch existing double strand breaks located in the top-strand's upstream regions, leading to the generation of large genomic deletions.

A large genomic deletion in the PDHX gene caused by the retrotranspositional insertion of a full-length LINE-1 element

The long interspersed element-1 (LINE-1 or L1) retrotransposition has altered the human genome in many ways. In particular, recent in vitro studies have demonstrated that the retrotranspositional insertion of L1 elements has resulted in significant genomic deletions. Here we provide evidence for its operation in the human genome by identifying a approximately 46-kb pathological genomic deletion in the PDHX gene directly linked to the insertion of a full-length L1 element, in a patient with pyruvate dehydrogenase complex (PDHc) deficiency. Both the deduced bottom and top strand cleavage sites in the PDHX gene coincide with the consensus L1 endonuclease (EN) target sequence 5'-TTTT/A-3', while the full-length L1 element is followed by a 67-bp poly(A) tail. Interestingly, two hairpin structures, potentially formed by the inverted repeats present immediately 5' to the top strand nick site and 3' to the bottom strand nick site, may have facilitated the accessibility of L1 EN to the target sequences and also brought the two otherwise distantly located sequences into close proximity. Since the L1 element inserted in the PDHX gene is full-length, we favor the model of the template jumping as opposed to that of the microhomology-mediated end-joining for linking the 5' end of the nascent L1 copy to its genomic target. Our finding not only serves as an important complement to the in vitro approaches to studying L1 retrotransposition, but also reveals a novel mechanism causing human genetic disease.

A novel insertion of a rearranged L1 element in exon 44 of the dystrophin gene: further evidence for possible bias in retroposon integration

L1 elements are mammalian retrotransposons contributing to genome evolution and causing rare mutations in human. We describe a de novo insertion of an L1 element into the dystrophin gene resulting in skipping of exon 44 and causing Duchenne muscular dystrophy in a boy. The L1 element was rearranged due to the twin-priming mechanism, but contrary to all described L1 rearrangements the 5' region of the inverted L1 sequence ended within the poly(A) tail of the element. Furthermore, the target site for the insertion was located only 87 bp from the insertion site in another patient described previously. These findings can contribute to the understanding of the mechanisms of L1 element rearrangement, and may support the notion that some subregions of the human genome could be preferred targets for retroelements using the L1 enzymatic machinery.

A systematic analysis of disease-associated variants in the 3' regulatory regions of human protein-coding genes I: general principles and overview

The 3' regulatory regions (3' RRs) of human genes play an important role in regulating mRNA 3' end formation, stability/degradation, nuclear export, subcellular localization and translation and are consequently rich in regulatory elements. Although 3' RRs contain only approximately 0.2% of known disease-associated mutations, this is likely to represent a rather conservative estimate of their actual prevalence. In an attempt to catalogue 3' RR-mediated disease and also to gain a greater understanding of the functional role of regulatory elements within 3' RRs, we have performed a systematic analysis of disease-associated 3' RR variants; 121 3' RR variants in 94 human genes were collated. These included 17 mutations in the upstream core polyadenylation signal sequence (UCPAS), 81 in the upstream sequence (USS) between the translational termination codon and the UCPAS, 6 in the left arm of the 'spacer' sequence (LAS) between the UCPAS and the pre-mRNA cleavage site (CS), 3 in the right arm of the 'spacer' sequence (RAS) or downstream core polyadenylation signal sequence (DCPAS) and 7 in the downstream sequence (DSS) of the 3'-flanking region, with 7 further mutations being treated as isolated examples. All the UCPAS mutations and the rather unusual cases of the DMPK, SCA8, FCMD and GLA mutations exert a significant effect on the mRNA phenotype and are usually associated with monogenic disease. By contrast, most of the remaining variants are polymorphisms that exert a comparatively minor influence on mRNA expression, but which may nevertheless predispose to or otherwise modify complex clinical phenotypes. Considerable efforts have been made to validate/elucidate the mechanisms through which the 3' untranslated region (3' UTR) variants affect gene expression. It is hoped that the integrative approach employed here in the study of naturally occurring variants of actual or potential pathological significance will serve to complement ongoing efforts to identify all functional regulatory elements in the human genome.