Active Alu element "A-tails": size does matter

Variable patterns of retrotransposition in different HeLa strains provide mechanistic insights into SINE RNA mobilization processes

Alu elements are non-autonomous Short INterspersed Elements (SINEs) derived from the 7SL RNA gene that are present at over one million copies in human genomic DNA. Alu mobilizes by a mechanism known as retrotransposition, which requires the Long INterspersed Element-1 (LINE-1) ORF2-encoded protein (ORF2p). Here, we demonstrate that HeLa strains differ in their capacity to support Alu retrotransposition. Human Alu elements retrotranspose efficiently in HeLa-HA and HeLa-CCL2 (Alu-permissive) strains, but not in HeLa-JVM or HeLa-H1 (Alu-nonpermissive) strains. A similar pattern of retrotransposition was observed for other 7SL RNA-derived SINEs and tRNA-derived SINEs. In contrast, mammalian LINE-1s, a zebrafish LINE, a human SINE-VNTR-Alu (SVA) element, and an L1 ORF1-containing mRNA can retrotranspose in all four HeLa strains. Using an in vitro reverse transcriptase-based assay, we show that Alu RNAs associate with ORF2p and are converted into cDNAs in both Alu-permissive and Alu-nonpermissive HeLa strains, suggesting that 7SL- and tRNA-derived SINEs use strategies to 'hijack' L1 ORF2p that are distinct from those used by SVA elements and ORF1-containing mRNAs. These data further suggest ORF2p associates with the Alu RNA poly(A) tract in both Alu-permissive and Alu-nonpermissive HeLa strains, but that Alu retrotransposition is blocked after this critical step in Alu-nonpermissive HeLa strains.

Long-read sequencing improves the genetic diagnosis of retinitis pigmentosa by identifying an Alu retrotransposon insertion in the EYS gene

BACKGROUND

Biallelic variants in EYS are the major cause of autosomal recessive retinitis pigmentosa (arRP) in certain populations, a clinically and genetically heterogeneous disease that may lead to legal blindness. EYS is one of the largest genes (~ 2 Mb) expressed in the retina, in which structural variants (SVs) represent a common cause of disease. However, their identification using short-read sequencing (SRS) is not always feasible. Here, we conducted targeted long-read sequencing (T-LRS) using adaptive sampling of EYS on the MinION sequencing platform (Oxford Nanopore Technologies) to definitively diagnose an arRP family, whose affected individuals (n = 3) carried the heterozygous pathogenic deletion of exons 32-33 in the EYS gene. As this was a recurrent variant identified in three additional families in our cohort, we also aimed to characterize the known deletion at the nucleotide level to assess a possible founder effect.

RESULTS

T-LRS in family A unveiled a heterozygous AluYa5 insertion in the coding exon 43 of EYS (chr6(GRCh37):g.64430524_64430525ins352), which segregated with the disease in compound heterozygosity with the previously identified deletion. Visual inspection of previous SRS alignments using IGV revealed several reads containing soft-clipped bases, accompanied by a slight drop in coverage at the Alu insertion site. This prompted us to develop a simplified program using grep command to investigate the recurrence of this variant in our cohort from SRS data. Moreover, LRS also allowed the characterization of the CNV as a ~ 56.4kb deletion spanning exons 32-33 of EYS (chr6(GRCh37):g.64764235_64820592del). The results of further characterization by Sanger sequencing and linkage analysis in the four families were consistent with a founder variant.

CONCLUSIONS

To our knowledge, this is the first report of a mobile element insertion into the coding sequence of EYS, as a likely cause of arRP in a family. Our study highlights the value of LRS technology in characterizing and identifying hidden pathogenic SVs, such as retrotransposon insertions, whose contribution to the etiopathogenesis of rare diseases may be underestimated.

Variable patterns of retrotransposition in different HeLa strains provide mechanistic insights into SINE RNA mobilization processes

Alu elements are non-autonomous Short INterspersed Elements (SINEs) derived from the 7SL RNA gene that are present at over one million copies in human genomic DNA. Alu mobilizes by a mechanism known as retrotransposition, which requires the Long INterspersed Element-1 (LINE-1 or L1) ORF2 -encoded protein (ORF2p). Here, we demonstrate that HeLa strains differ in their capacity to support Alu retrotransposition. Human Alu elements retrotranspose efficiently in HeLa-HA and HeLa-CCL2 ( Alu -permissive) strains, but not in HeLa-JVM or HeLa-H1 ( Alu -nonpermissive) strains. A similar pattern of retrotransposition was observed for other 7SL RNA -derived SINEs and tRNA -derived SINEs. In contrast, mammalian LINE-1s, a zebrafish LINE, a human SINE-VNTR - Alu ( SVA ) element, and an L1 ORF1 -containing messenger RNA can retrotranspose in all four HeLa strains. Using an in vitro reverse transcriptase-based assay, we show that Alu RNAs associate with ORF2p and are converted into cDNAs in both Alu -permissive and Alu -nonpermissive HeLa strains, suggesting that 7SL - and tRNA -derived SINE RNAs use strategies to 'hijack' L1 ORF2p that are distinct from those used by SVA elements and ORF1 -containing mRNAs. These data further suggest ORF2p associates with the Alu RNA poly(A) tract in both Alu -permissive and Alu -nonpermissive HeLa strains, but that Alu retrotransposition is blocked after this critical step in Alu -nonpermissive HeLa strains.

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.

Extensive Independent Amplification of Platy-1 Retroposons in Tamarins, Genus Saguinus

Platy-1 retroposons are short interspersed elements (SINEs) unique to platyrrhine primates. Discovered in the common marmoset (Callithrix jacchus) genome in 2016, these 100 bp mobile element insertions (MEIs) appeared to be novel drivers of platyrrhine evolution, with over 2200 full-length members across 62 different subfamilies, and strong evidence of ongoing proliferation in C. jacchus. Subsequent characterization of Platy-1 elements in Aotus, Saimiri and Cebus genera, suggested that the widespread mobilization detected in marmoset (family Callithrichidae) was perhaps an anomaly. Two additional Callithrichidae genomes are now available, a scaffold level genome assembly for Saguinus imperator (tamarin; SagImp_v1) and a chromosome-level assembly for Saguinus midas (Midas tamarin; ASM2_v1). Here, we report that each tamarin genome contains over 11,000 full-length Platy-1 insertions, about 1150 are shared by both Saguinus tamarins, 7511 are unique to S. imperator, and another 8187 are unique to S. midas. Roughly 325 are shared among the three callithrichids. We identified six new Platy-1 subfamilies derived from Platy-1-8, with the youngest new subfamily, Platy-1-8c_Saguinus, being the primary source of the Saguinus amplification burst. This constitutes the largest expansion of Platy-1 MEIs reported to date and the most extensive independent SINE amplification between two closely related species.

Disruption of FBN1 by an Alu element insertion: A novel genetic cause of Marfan syndrome

Alu elements are retrotransposons with ubiquitous presence in the human genome that have contributed to human genomic diversity and health. These approximately 300-bp sequences can cause or mediate disease by disrupting coding/splicing regions in the germline, by insertional mutagenesis in somatic cells, and in promoting formation of copy-number variants. Alu elements may also disrupt epigenetic regulation by affecting non-coding regulatory regions. There are increasing reports of apparently sporadic and inherited genetic disorders caused by Alu-related gene disruption, but Marfan syndrome resulting from Alu element insertion has not been previously described. We report a family with classic features of Marfan syndrome whose previous FBN1 genetic testing was inconclusive. Using contemporary next-generation sequencing and bioinformatics analysis, a pathogenic/disruptive Alu insertion occurring in the coding region of the FBN1 gene was identified (c.6564_6565insAlu; p. Glu2189fs) and was confirmed and specified further with Sanger sequencing. This identified the molecular basis of disease in the family that was missed using previous genetic testing technologies and highlights a novel pathogenic mechanism for Marfan syndrome. This case adds to the growing literature of Mendelian diseases caused by Alu retrotransposition, and it also shows the growing capability of genomic technologies for detecting atypical mutation events.

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.

The genetic and molecular features of the intronic pentanucleotide repeat expansion in spinocerebellar ataxia type 10

Spinocerebellar ataxia type 10 (SCA10) is characterized by progressive cerebellar neurodegeneration and, in many patients, epilepsy. This disease mainly occurs in individuals with Indigenous American or East Asian ancestry, with strong evidence supporting a founder effect. The mutation causing SCA10 is a large expansion in an ATTCT pentanucleotide repeat in intron 9 of the ATXN10 gene. The ATTCT repeat is highly unstable, expanding to 280–4,500 repeats in affected patients compared with the 9–32 repeats in normal individuals, one of the largest repeat expansions causing neurological disorders identified to date. However, the underlying molecular basis of how this huge repeat expansion evolves and contributes to the SCA10 phenotype remains largely unknown. Recent progress in next-generation DNA sequencing technologies has established that the SCA10 repeat sequence has a highly heterogeneous structure. Here we summarize what is known about the structure and origin of SCA10 repeats, discuss the potential contribution of variant repeats to the SCA10 disease phenotype, and explore how this information can be exploited for therapeutic benefit.

Alu element in the RNA binding motif protein, X-linked 2 (RBMX2) gene found to be linked to bipolar disorder

Objective

We have used long-read single molecule, real-time (SMRT) sequencing to fully characterize a ~12Mb genomic region on chromosome Xq24-q27, significantly linked to bipolar disorder (BD) in an extended family from a genetic sub-isolate. This family segregates BD in at least four generations with 24 affected individuals.

Methods

We selected 16 family members for targeted sequencing. The selected individuals either carried the disease haplotype, were non-carriers of the disease haplotype, or served as married-in controls. We designed hybrid capture probes enriching for 5-9Kb fragments spanning the entire 12Mb region that were then sequenced to screen for candidate structural variants (SVs) that could explain the increased risk for BD in this extended family.

Results

Altogether, 201 variants were detected in the critically linked region. Although most of these represented common variants, three variants emerged that showed near-perfect segregation among all BD type I affected individuals. Two of the SVs were identified in or near genes belonging to the RNA Binding Motif Protein, X-Linked (RBMX) gene family—a 330bp Alu (subfamily AluYa5) deletion in intron 3 of the RBMX2 gene and an intergenic 27bp tandem repeat deletion between the RBMX and G protein-coupled receptor 101 (GPR101) genes. The third SV was a 50bp tandem repeat insertion in intron 1 of the Coagulation Factor IX (F9) gene.

Conclusions

Among the three genetically linked SVs, additional evidence supported the Alu element deletion in RBMX2 as the leading candidate for contributing directly to the disease development of BD type I in this extended family.

The nuclear and cytoplasmic activities of RNA polymerase III, and an evolving transcriptome for surveillance

A 1969 report that described biochemical and activity properties of the three eukaryotic RNA polymerases revealed Pol III as highly distinguishable, even before its transcripts were identified. Now known to be the most complex, Pol III contains several stably-associated subunits referred to as built-in transcription factors (BITFs) that enable highly efficient RNA synthesis by a unique termination-associated recycling process. In vertebrates, subunit RPC7(α/β) can be of two forms, encoded by POLR3G or POLR3GL, with differential activity. Here we review promoter-dependent transcription by Pol III as an evolutionary perspective of eukaryotic tRNA expression. Pol III also provides nonconventional functions reportedly by promoter-independent transcription, one of which is RNA synthesis from DNA 3'-ends during repair. Another is synthesis of 5'ppp-RNA signaling molecules from cytoplasmic viral DNA in a pathway of interferon activation that is dysfunctional in immunocompromised patients with mutations in Pol III subunits. These unconventional functions are also reviewed, including evidence that link them to the BITF subunits. We also review data on a fraction of the human Pol III transcriptome that evolved to include vault RNAs and snaRs with activities related to differentiation, and in innate immune and tumor surveillance. The Pol III of higher eukaryotes does considerably more than housekeeping.

Post-transcriptional regulation in spermatogenesis: all RNA pathways lead to healthy sperm

Multiple RNA pathways are required to produce functional sperm. Here, we review RNA post-transcriptional regulation during spermatogenesis with particular emphasis on the role of 3' end modifications. From early studies in the 1970s, it became clear that spermiogenesis transcripts could be stored for days only to be translated at advanced stages of spermatid differentiation. The transition between the translationally repressed and active states was observed to correlate with the shortening of the transcripts' poly(A) tail, establishing a link between RNA 3' end metabolism and male germ cell differentiation. Since then, numerous RNA metabolic pathways have been implicated not only in the progression through spermatogenesis, but also in the maintenance of genomic integrity. Recent studies have characterized the elusive 3' biogenesis of Piwi-interacting RNAs (piRNAs), identified a critical role for messenger RNA (mRNA) 3' uridylation in meiotic progression, established the mechanisms that destabilize transcripts with long 3' untranslated regions (3'UTRs) in post-mitotic cells, and defined the physiological relevance of RNA exonucleases and deadenylases in male germ cells. In this review, we discuss RNA processing in the male germline in the light of the most recent findings. A brief recollection of different RNA-processing events will aid future studies exploring post-transcriptional regulation in spermatogenesis.

Analysis of SINE Families B2, Dip, and Ves with Special Reference to Polyadenylation Signals and Transcription Terminators

Short Interspersed Elements (SINEs) are eukaryotic non-autonomous retrotransposons transcribed by RNA polymerase III (pol III). The 3′-terminus of many mammalian SINEs has a polyadenylation signal (AATAAA), pol III transcription terminator, and A-rich tail. The RNAs of such SINEs can be polyadenylated, which is unique for pol III transcripts. Here, B2 (mice and related rodents), Dip (jerboas), and Ves (vespertilionid bats) SINE families were thoroughly studied. They were divided into subfamilies reliably distinguished by relatively long indels. The age of SINE subfamilies can be estimated, which allows us to reconstruct their evolution. The youngest and most active variants of SINE subfamilies were given special attention. The shortest pol III transcription terminators are TCTTT (B2), TATTT (Ves and Dip), and the rarer TTTT. The last nucleotide of the terminator is often not transcribed; accordingly, the truncated terminator of its descendant becomes nonfunctional. The incidence of complete transcription of the TCTTT terminator is twice higher compared to TTTT and thus functional terminators are more likely preserved in daughter SINE copies. Young copies have long poly(A) tails; however, they gradually shorten in host generations. Unexpectedly, the tail shortening below A₁₀ increases the incidence of terminator elongation by Ts thus restoring its efficiency. This process can be critical for the maintenance of SINE activity in the genome.

Haplotype-resolved diverse human genomes and integrated analysis of structural variation

Long-read and strand-specific sequencing technologies together facilitate the de novo assembly of high-quality haplotype-resolved human genomes without parent-child trio data. We present 64 assembled haplotypes from 32 diverse human genomes. These highly contiguous haplotype assemblies (average minimum contig length needed to cover 50% of the genome: 26 million base pairs) integrate all forms of genetic variation, even across complex loci. We identified 107,590 structural variants (SVs), of which 68% were not discovered with short-read sequencing, and 278 SV hotspots (spanning megabases of gene-rich sequence). We characterized 130 of the most active mobile element source elements and found that 63% of all SVs arise through homology-mediated mechanisms. This resource enables reliable graph-based genotyping from short reads of up to 50,340 SVs, resulting in the identification of 1526 expression quantitative trait loci as well as SV candidates for adaptive selection within the human population.

Diversity of short interspersed nuclear elements (SINEs) in lepidopteran insects and evidence of horizontal SINE transfer between baculovirus and lepidopteran hosts

Background

Short interspersed nuclear elements (SINEs) belong to non-long terminal repeat (non-LTR) retrotransposons, which can mobilize dependent on the help of counterpart long interspersed nuclear elements (LINEs). Although 234 SINEs have been identified so far, only 23 are from insect species (SINEbase: http://sines.eimb.ru/).

Results

Here, five SINEs were identified from the genome of Plutella xylostella, among which PxSE1, PxSE2 and PxSE3 were tRNA-derived SINEs, PxSE4 and PxSE5 were 5S RNA-derived SINEs. A total of 18 related SINEs were further identified in 13 lepidopteran insects and a baculovirus. The 3′-tail of PxSE5 shares highly identity with that of LINE retrotransposon, PxLINE1. The analysis of relative age distribution profiles revealed that PxSE1 is a relatively young retrotransposon in the genome of P. xylostella and was generated by recent explosive amplification. Integration pattern analysis showed that SINEs in P. xylostella prefer to insert into or accumulate in introns and regions 5 kb downstream of genes. In particular, the PxSE1-like element, SlNPVSE1, in Spodoptera litura nucleopolyhedrovirus II genome is highly identical to SfSE1 in Spodoptera frugiperda, SlittSE1 in Spodoptera littoralis, and SlituSE1 in Spodoptera litura, suggesting the occurrence of horizontal transfer.

Conclusions

Lepidopteran insect genomes harbor a diversity of SINEs. The retrotransposition activity and copy number of these SINEs varies considerably between host lineages and SINE lineages. Host-parasite interactions facilitate the horizontal transfer of SINE between baculovirus and its lepidopteran hosts.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07543-z.

Analysis of HLA-G long-read genomic sequences in mother-offspring pairs with preeclampsia

Preeclampsia is a pregnancy-induced disorder that is characterized by hypertension and is a leading cause of perinatal and maternal–fetal morbidity and mortality. HLA-G is thought to play important roles in maternal–fetal immune tolerance, and the associations between HLA-G gene polymorphisms and the onset of pregnancy-related diseases have been explored extensively. Because contiguous genomic sequencing is difficult, the association between the HLA-G genotype and preeclampsia onset is controversial. In this study, genomic sequences of the HLA-G region (5.2 kb) from 31 pairs of mother–offspring genomic DNA samples (18 pairs from normal pregnancies/births and 13 from preeclampsia births) were obtained by single-molecule real-time sequencing using the PacBio RS II platform. The HLA-G alleles identified in our cohort matched seven known HLA-G alleles, but we also identified two new HLA-G alleles at the fourth-field resolution and compared them with nucleotide sequences from a public database that consisted of coding sequences that cover the 3.1-kb HLA-G gene span. Intriguingly, a potential association between preeclampsia onset and the poly T stretch within the downstream region of the HLA-G*01:01:01:01 allele was found. Our study suggests that long-read sequencing of HLA-G will provide clues for characterizing HLA-G variants that are involved in the pathophysiology of preeclampsia.

Divergence of 3' ends as a driver of short interspersed nuclear element (SINE) evolution in the Salicaceae

Short interspersed nuclear elements (SINEs) are small, non-autonomous and heterogeneous retrotransposons that are widespread in plants. To explore the amplification dynamics and evolutionary history of SINE populations in representative deciduous tree species, we analyzed the genomes of the six following Salicaceae species: Populus deltoides, Populus euphratica, Populus tremula, Populus tremuloides, Populus trichocarpa, and Salix purpurea. We identified 11 Salicaceae SINE families (SaliS-I to SaliS-XI), comprising 27 077 full-length copies. Most of these families harbor segmental similarities, providing evidence for SINE emergence by reshuffling or heterodimerization. We observed two SINE groups, differing in phylogenetic distribution pattern, similarity and 3' end structure. These groups probably emerged during the 'salicoid duplication' (~65 million years ago) in the Salix-Populus progenitor and during the separation of the genus Salix (45-65 million years ago), respectively. In contrast to conserved 5' start motifs across species and SINE families, the 3' ends are highly variable in sequence and length. This extraordinary 3'-end variability results from mutations in the poly(A) tail, which were fixed by subsequent amplificational bursts. We show that the dissemination of newly evolved 3' ends is accomplished by a displacement of older motifs, leading to various 3'-end subpopulations within the SaliS families.

Insertion of Telomeric Repeats in the Human and Horse Genomes: An Evolutionary Perspective

Interstitial telomeric sequences (ITSs) are short stretches of telomeric-like repeats (TTAGGG)n at nonterminal chromosomal sites. We previously demonstrated that, in the genomes of primates and rodents, ITSs were inserted during the repair of DNA double-strand breaks. These conclusions were derived from sequence comparisons of ITS-containing loci and ITS-less orthologous loci in different species. To our knowledge, insertion polymorphism of ITSs, i.e., the presence of an ITS-containing allele and an ITS-less allele in the same species, has not been described. In this work, we carried out a genome-wide analysis of 2504 human genomic sequences retrieved from the 1000 Genomes Project and a PCR-based analysis of 209 human DNA samples. In spite of the large number of individual genomes analyzed we did not find any evidence of insertion polymorphism in the human population. On the contrary, the analysis of ITS loci in the genome of a single horse individual, the reference genome, allowed us to identify five heterozygous ITS loci, suggesting that insertion polymorphism of ITSs is an important source of genetic variability in this species. Finally, following a comparative sequence analysis of horse ITSs and of their orthologous empty loci in other Perissodactyla, we propose models for the mechanism of ITS insertion during the evolution of this order.

Finding and extending ancient simple sequence repeat-derived regions in the human genome

Background

Previously, 3% of the human genome has been annotated as simple sequence repeats (SSRs), similar to the proportion annotated as protein coding. The origin of much of the genome is not well annotated, however, and some of the unidentified regions are likely to be ancient SSR-derived regions not identified by current methods. The identification of these regions is complicated because SSRs appear to evolve through complex cycles of expansion and contraction, often interrupted by mutations that alter both the repeated motif and mutation rate. We applied an empirical, kmer-based, approach to identify genome regions that are likely derived from SSRs.

Results

The sequences flanking annotated SSRs are enriched for similar sequences and for SSRs with similar motifs, suggesting that the evolutionary remains of SSR activity abound in regions near obvious SSRs. Using our previously described P-clouds approach, we identified ‘SSR-clouds’, groups of similar kmers (or ‘oligos’) that are enriched near a training set of unbroken SSR loci, and then used the SSR-clouds to detect likely SSR-derived regions throughout the genome.

Conclusions

Our analysis indicates that the amount of likely SSR-derived sequence in the human genome is 6.77%, over twice as much as previous estimates, including millions of newly identified ancient SSR-derived loci. SSR-clouds identified poly-A sequences adjacent to transposable element termini in over 74% of the oldest class of Alu (roughly, AluJ), validating the sensitivity of the approach. Poly-A’s annotated by SSR-clouds also had a length distribution that was more consistent with their poly-A origins, with mean about 35 bp even in older Alus. This work demonstrates that the high sensitivity provided by SSR-Clouds improves the detection of SSR-derived regions and will enable deeper analysis of how decaying repeats contribute to genome structure.

Genome-wide analysis of short interspersed nuclear elements provides insight into gene and genome evolution in citrus

Short interspersed nuclear elements (SINEs) are non-autonomous retrotransposons that are highly abundant, but not well annotated, in plant genomes. In this study, we identified 41,573 copies of SINEs in seven citrus genomes, including 11,275 full-length copies. The citrus SINEs were distributed among 12 families, with an average full-length rate of 0.27, and were dispersed throughout the chromosomes, preferentially in AT-rich areas. Approximately 18.4% of citrus SINEs were found in close proximity (≤1 kb upstream) to genes, indicating a significant enrichment of SINEs in promoter regions. Citrus SINEs promote gene and genome evolution by offering exons as well as splice sites and start and stop codons, creating novel genes and forming tandem and dispersed repeat structures. Comparative analysis of unique homologous SINE-containing loci (HSCLs) revealed chromosome rearrangements in sweet orange, pummelo, and mandarin, suggesting that unique HSCLs might be valuable for understanding chromosomal abnormalities. This study of SINEs provides us with new perspectives and new avenues by which to understand the evolution of citrus genes and genomes.

Simultaneous TE Analysis of 19 Heliconiine Butterflies Yields Novel Insights into Rapid TE-Based Genome Diversification and Multiple SINE Births and Deaths

Transposable elements (TEs) play major roles in the evolution of genome structure and function. However, because of their repetitive nature, they are difficult to annotate and discovering the specific roles they may play in a lineage can be a daunting task. Heliconiine butterflies are models for the study of multiple evolutionary processes including phenotype evolution and hybridization. We attempted to determine how TEs may play a role in the diversification of genomes within this clade by performing a detailed examination of TE content and accumulation in 19 species whose genomes were recently sequenced. We found that TE content has diverged substantially and rapidly in the time since several subclades shared a common ancestor with each lineage harboring a unique TE repertoire. Several novel SINE lineages have been established that are restricted to a subset of species. Furthermore, the previously described SINE, Metulj, appears to have gone extinct in two subclades while expanding to significant numbers in others. This diversity in TE content and activity has the potential to impact how heliconiine butterflies continue to evolve and diverge.

Selective elimination of long INterspersed element-1 expressing tumour cells by targeted expression of the HSV-TK suicide gene

In gene therapy, effective and selective suicide gene expression is crucial. We exploited the endogenous Long INterspersed Element-1 (L1) machinery often reactivated in human cancers to integrate the Herpes Simplex Virus Thymidine Kinase (HSV-TK) suicide gene selectively into the genome of cancer cells. We developed a plasmid-based system directing HSV-TK expression only when reverse transcribed and integrated in the host genome via the endogenous L1 ORF1/2 proteins and an Alu element. Delivery of these new constructs into cells followed by Ganciclovir (GCV) treatment selectively induced mortality of L1 ORF1/2 protein expressing cancer cells, but had no effect on primary cells that do not express L1 ORF1/2. This novel strategy for selective targeting of tumour cells provides high tolerability as the HSV-TK gene cannot be expressed without reverse transcription and integration, and high selectivity as these processes take place only in cancer cells expressing high levels of functional L1 ORF1/2.

Transposable elements: genome innovation, chromosome diversity, and centromere conflict

Although it was nearly 70 years ago when transposable elements (TEs) were first discovered “jumping” from one genomic location to another, TEs are now recognized as contributors to genomic innovations as well as genome instability across a wide variety of species. In this review, we illustrate the ways in which active TEs, specifically retroelements, can create novel chromosome rearrangements and impact gene expression, leading to disease in some cases and species-specific diversity in others. We explore the ways in which eukaryotic genomes have evolved defense mechanisms to temper TE activity and the ways in which TEs continue to influence genome structure despite being rendered transpositionally inactive. Finally, we focus on the role of TEs in the establishment, maintenance, and stabilization of critical, yet rapidly evolving, chromosome features: eukaryotic centromeres. Across centromeres, specific types of TEs participate in genomic conflict, a balancing act wherein they are actively inserting into centromeric domains yet are harnessed for the recruitment of centromeric histones and potentially new centromere formation.

Insertion and deletion polymorphisms of the ancient AluS family in the human genome

Background

Polymorphic Alu elements account for 17% of structural variants in the human genome. The majority of these belong to the youngest AluY subfamilies, and most structural variant discovery efforts have focused on identifying Alu polymorphisms from these currently retrotranspositionally active subfamilies. In this report we analyze polymorphisms from the evolutionarily older AluS subfamily, whose peak activity was tens of millions of years ago. We annotate the AluS polymorphisms, assess their likely mechanism of origin, and evaluate their contribution to structural variation in the human genome.

Results

Of 52 previously reported polymorphic AluS elements ascertained for this study, 48 were confirmed to belong to the AluS subfamily using high stringency subfamily classification criteria. Of these, the majority (77%, 37/48) appear to be deletion polymorphisms. Two polymorphic AluS elements (4%) have features of non-classical Alu insertions and one polymorphic AluS element (2%) likely inserted by a mechanism involving internal priming. Seven AluS polymorphisms (15%) appear to have arisen by the classical target-primed reverse transcription (TPRT) retrotransposition mechanism. These seven TPRT products are 3′ intact with 3′ poly-A tails, and are flanked by target site duplications; L1 ORF2p endonuclease cleavage sites were also observed, providing additional evidence that these are L1 ORF2p endonuclease-mediated TPRT insertions. Further sequence analysis showed strong conservation of both the RNA polymerase III promoter and SRP9/14 binding sites, important for mediating transcription and interaction with retrotransposition machinery, respectively. This conservation of functional features implies that some of these are fairly recent insertions since they have not diverged significantly from their respective retrotranspositionally competent source elements.

Conclusions

Of the polymorphic AluS elements evaluated in this report, 15% (7/48) have features consistent with TPRT-mediated insertion, thus suggesting that some AluS elements have been more active recently than previously thought, or that fixation of AluS insertion alleles remains incomplete. These data expand the potential significance of polymorphic AluS elements in contributing to structural variation in the human genome. Future discovery efforts focusing on polymorphic AluS elements are likely to identify more such polymorphisms, and approaches tailored to identify deletion alleles may be warranted.

Electronic supplementary material

The online version of this article (doi:10.1186/s13100-017-0089-9) contains supplementary material, which is available to authorized users.

Conserved 3' UTR stem-loop structure in L1 and Alu transposons in human genome: possible role in retrotransposition

BACKGROUND

In the process of retrotransposition LINEs use their own machinery for copying and inserting themselves into new genomic locations, while SINEs are parasitic and require the machinery of LINEs. The exact mechanism of how a LINE-encoded reverse transcriptase (RT) recognizes its own and SINE RNA remains unclear. However it was shown for the stringent-type LINEs that recognition of a stem-loop at the 3'UTR by RT is essential for retrotransposition. For the relaxed-type LINEs it is believed that the poly-A tail is a common recognition element between LINE and SINE RNA. However polyadenylation is a property of any messenger RNA, and how the LINE RT recognizes transposon and non-transposon RNAs remains an open question. It is likely that RNA secondary structures play an important role in RNA recognition by LINE encoded proteins.

RESULTS

Here we selected a set of L1 and Alu elements from the human genome and investigated their sequences for the presence of position-specific stem-loop structures. We found highly conserved stem-loop positions at the 3'UTR. Comparative structural analyses of a human L1 3'UTR stem-loop showed a similarity to 3'UTR stem-loops of the stringent-type LINEs, which were experimentally shown to be recognized by LINE RT. The consensus stem-loop structure consists of 5-7 bp loop, 8-10 bp stem with a bulge at a distance of 4-6 bp from the loop. The results show that a stem loop with a bulge exists at the 3'-end of Alu. We also found conserved stem-loop positions at 5'UTR and at the end of ORF2 and discuss their possible role.

CONCLUSIONS

Here we presented an evidence for the presence of a highly conserved 3'UTR stem-loop structure in L1 and Alu retrotransposons in the human genome. Both stem-loops show structural similarity to the stem-loops of the stringent-type LINEs experimentally confirmed as essential for retrotransposition. Here we hypothesize that both L1 and Alu RNA are recognized by L1 RT via the 3'-end RNA stem-loop structure. Other conserved stem-loop positions in L1 suggest their possible functions in protein-RNA interactions but to date no experimental evidence has been reported.

Short interspersed nuclear elements (SINEs) are abundant in Solanaceae and have a family-specific impact on gene structure and genome organization

Short interspersed nuclear elements (SINEs) are highly abundant non-autonomous retrotransposons that are widespread in plants. They are short in size, non-coding, show high sequence diversity, and are therefore mostly not or not correctly annotated in plant genome sequences. Hence, comparative studies on genomic SINE populations are rare. To explore the structural organization and impact of SINEs, we comparatively investigated the genome sequences of the Solanaceae species potato (Solanum tuberosum), tomato (Solanum lycopersicum), wild tomato (Solanum pennellii), and two pepper cultivars (Capsicum annuum). Based on 8.5 Gbp sequence data, we annotated 82 983 SINE copies belonging to 10 families and subfamilies on a base pair level. Solanaceae SINEs are dispersed over all chromosomes with enrichments in distal regions. Depending on the genome assemblies and gene predictions, 30% of all SINE copies are associated with genes, particularly frequent in introns and untranslated regions (UTRs). The close association with genes is family specific. More than 10% of all genes annotated in the Solanaceae species investigated contain at least one SINE insertion, and we found genes harbouring up to 16 SINE copies. We demonstrate the involvement of SINEs in gene and genome evolution including the donation of splice sites, start and stop codons and exons to genes, enlargement of introns and UTRs, generation of tandem-like duplications and transduction of adjacent sequence regions.

Get in LINE: Competition for Newly Minted Retrotransposon Proteins at the Ribosome

In this issue, Ahl et al. (2015) and Doucet et al. (2015) illuminate structural and functional features of substrates that promote integration of RNA molecules into the human genome by LINE retrotransposons, contributing to the ∼ 50% of the human genome that has been colonized by mobile genetic elements.

Discovery of a new repeat family in the Callithrix jacchus genome

We identified a novel repeat family, termed Platy-1, in the Callithrix jacchus (common marmoset) genome that arose around the time of the divergence of platyrrhines and catarrhines and established itself as a repeat family in New World monkeys (NWMs). A full-length Platy-1 element is ∼100 bp in length, making it the shortest known short interspersed element (SINE) in primates, and harbors features characteristic of non-LTR retrotransposons. We identified 2268 full-length Platy-1 elements across 62 subfamilies in the common marmoset genome. Our subfamily reconstruction and phylogenetic analyses support Platy-1 propagation throughout the evolution of NWMs in the lineage leading to C. jacchus Platy-1 appears to have reached its amplification peak in the common ancestor of current day marmosets and has since moderately declined. However, identification of more than 200 Platy-1 elements identical to their respective consensus sequence, and the presence of polymorphic elements within common marmoset populations, suggests ongoing retrotransposition activity. Platy-1, a SINE, appears to have originated from an Alu element, and hence is likely derived from 7SL RNA. Our analyses illustrate the birth of a new repeat family and its propagation dynamics in the lineage leading to the common marmoset over the last 40 million years.

DNA methylation variation of human-specific Alu repeats

DNA methylation is the major repression mechanism for human retrotransposons, such as the Alu family. Here, we have determined the methylation levels associated with 5238 loci belonging to 2 Alu subfamilies, AluYa5 and AluYb8, using high-throughput targeted repeat element bisulfite sequencing (HT-TREBS). The results indicate that ∼90% of loci are repressed by high methylation levels. Of the remaining loci, many of the hypomethylated elements are found near gene promoters and show high levels of DNA methylation variation. We have characterized this variation in the context of tumorigenesis and interindividual differences. Comparison of a primary breast tumor and its matched normal tissue revealed early DNA methylation changes in ∼1% of AluYb8 elements in response to tumorigenesis. Simultaneously, AluYa5/Yb8 elements proximal to promoters also showed differences in methylation of up to one order of magnitude, even between normal individuals. Overall, the current study demonstrates that early loss of methylation occurs during tumorigenesis in a subset of young Alu elements, suggesting their potential clinical relevance. However, approaches such as deep-bisulfite-sequencing of individual loci using HT-TREBS are required to distinguish clinically relevant loci from the background observed for AluYa5/Yb8 elements in general with regard to high levels of interindividual variation in DNA methylation.

Insertion of Retrotransposons at Chromosome Ends: Adaptive Response to Chromosome Maintenance

The telomerase complex is a specialized reverse transcriptase (RT) that inserts tandem DNA arrays at the linear chromosome ends and contributes to the protection of the genetic information in eukaryotic genomes. Telomerases are phylogenetically related to retrotransposons, encoding also the RT activity required for the amplification of their sequences throughout the genome. Intriguingly the telomerase gene is lost from the Drosophila genome and tandem retrotransposons replace telomeric sequences at the chromosome extremities. This observation suggests the versatility of RT activity in counteracting the chromosome shortening associated with genome replication and that retrotransposons can provide this activity in case of a dysfunctional telomerase. In this review paper, we describe the major classes of retroelements present in eukaryotic genomes in order to point out the differences and similarities with the telomerase complex. In a second part, we discuss the insertion of retroelements at the ends of chromosomes as an adaptive response for dysfunctional telomeres.

Diversification, evolution and methylation of short interspersed nuclear element families in sugar beet and related Amaranthaceae species

Short interspersed nuclear elements (SINEs) are non-autonomous non-long terminal repeat retrotransposons which are widely distributed in eukaryotic organisms. While SINEs have been intensively studied in animals, only limited information is available about plant SINEs. We analysed 22 SINE families from seven genomes of the Amaranthaceae family and identified 34 806 SINEs, including 19 549 full-length copies. With the focus on sugar beet (Beta vulgaris), we performed a comparative analysis of the diversity, genomic and chromosomal organization and the methylation of SINEs to provide a detailed insight into the evolution and age of Amaranthaceae SINEs. The lengths of consensus sequences of SINEs range from 113 nucleotides (nt) up to 224 nt. The SINEs show dispersed distribution on all chromosomes but were found with higher incidence in subterminal euchromatic chromosome regions. The methylation of SINEs is increased compared with their flanking regions, and the strongest effect is visible for cytosines in the CHH context, indicating an involvement of asymmetric methylation in the silencing of SINEs.

Transposable element detection from whole genome sequence data

The number of software tools available for detecting transposable element insertions from whole genome sequence data has been increasing steadily throughout the last ~5 years. Some of these methods have unique features suiting them for particular use cases, but in general they follow one or more of a common set of approaches. Here, detection and filtering approaches are reviewed in the light of transposable element biology and the current state of whole genome sequencing. We demonstrate that the current state-of-the-art methods still do not produce highly concordant results and provide resources to assist future development in transposable element detection methods.

Polyadenylation of RNA transcribed from mammalian SINEs by RNA polymerase III: Complex requirements for nucleotide sequences

It is generally accepted that only transcripts synthesized by RNA polymerase II (e.g., mRNA) were subject to AAUAAA-dependent polyadenylation. However, we previously showed that RNA transcribed by RNA polymerase III (pol III) from mouse B2 SINE could be polyadenylated in an AAUAAA-dependent manner. Many species of mammalian SINEs end with the pol III transcriptional terminator (TTTTT) and contain hexamers AATAAA in their A-rich tail. Such SINEs were united into Class T(+), whereas SINEs lacking the terminator and AATAAA sequences were classified as T(-). Here we studied the structural features of SINE pol III transcripts that are necessary for their polyadenylation. Eight and six SINE families from classes T(+) and T(-), respectively, were analyzed. The replacement of AATAAA with AACAAA in T(+) SINEs abolished the RNA polyadenylation. Interestingly, insertion of the polyadenylation signal (AATAAA) and pol III transcription terminator in T(-) SINEs did not result in polyadenylation. The detailed analysis of three T(+) SINEs (B2, DIP, and VES) revealed areas important for the polyadenylation of their pol III transcripts: the polyadenylation signal and terminator in A-rich tail, β region positioned immediately downstream of the box B of pol III promoter, and τ region located upstream of the tail. In DIP and VES (but not in B2), the τ region is a polypyrimidine motif which is also characteristic of many other T(+) SINEs. Most likely, SINEs of different mammals acquired these structural features independently as a result of parallel evolution.

Sequence Analysis and Characterization of Active Human Alu Subfamilies Based on the 1000 Genomes Pilot Project

The goal of the 1000 Genomes Consortium is to characterize human genome structural variation (SV), including forms of copy number variations such as deletions, duplications, and insertions. Mobile element insertions, particularly Alu elements, are major contributors to genomic SV among humans. During the pilot phase of the project we experimentally validated 645 (611 intergenic and 34 exon targeted) polymorphic "young" Alu insertion events, absent from the human reference genome. Here, we report high resolution sequencing of 343 (322 unique) recent Alu insertion events, along with their respective target site duplications, precise genomic breakpoint coordinates, subfamily assignment, percent divergence, and estimated A-rich tail lengths. All the sequenced Alu loci were derived from the AluY lineage with no evidence of retrotransposition activity involving older Alu families (e.g., AluJ and AluS). AluYa5 is currently the most active Alu subfamily in the human lineage, followed by AluYb8, and many others including three newly identified subfamilies we have termed AluYb7a3, AluYb8b1, and AluYa4a1. This report provides the structural details of 322 unique Alu variants from individual human genomes collectively adding about 100 kb of genomic variation. Many Alu subfamilies are currently active in human populations, including a surprising level of AluY retrotransposition. Human Alu subfamilies exhibit continuous evolution with potential drivers sprouting new Alu lineages.

A genome survey sequencing of the Java mouse deer (Tragulus javanicus) adds new aspects to the evolution of lineage specific retrotransposons in Ruminantia (Cetartiodactyla)

Ruminantia, the ruminating, hoofed mammals (cow, deer, giraffe and allies) are an unranked artiodactylan clade. Around 50-60 million years ago the BovB retrotransposon entered the ancestral ruminantian genome through horizontal gene transfer. A survey genome screen using 454-pyrosequencing of the Java mouse deer (Tragulus javanicus) and the lesser kudu (Tragelaphus imberbis) was done to investigate and to compare the landscape of transposable elements within Ruminantia. The family Tragulidae (mouse deer) is the only representative of Tragulina and phylogenetically important, because it represents the earliest divergence in Ruminantia. The data analyses show that, relative to other ruminantian species, the lesser kudu genome has seen an expansion of BovB Long INterspersed Elements (LINEs) and BovB related Short INterspersed Elements (SINEs) like BOVA2. In comparison the genome of Java mouse deer has fewer BovB elements than other ruminants, especially Bovinae, and has in addition a novel CHR-3 SINE most likely propagated by LINE-1. By contrast the other ruminants have low amounts of CHR SINEs but high numbers of actively propagating BovB-derived and BovB-propagated SINEs. The survey sequencing data suggest that the transposable element landscape in mouse deer (Tragulina) is unique among Ruminantia, suggesting a lineage specific evolutionary trajectory that does not involve BovB mediated retrotransposition. This shows that the genomic landscape of mobile genetic elements can rapidly change in any lineage.

U6 snRNA Pseudogenes: Markers of Retrotransposition Dynamics in Mammals

Transposable elements comprise more than 45% of the human genome and long interspersed nuclear element 1 (LINE-1 or L1) is the only autonomous mobile element remaining active. Since its identification, it has been proposed that L1 contributes to the mobilization and amplification of other cellular RNAs and more recently, experimental demonstrations of this function has been described for many transcripts such as Alu, a nonautonomous mobile element, cellular mRNAs, or small noncoding RNAs. Detailed examination of the mobilization of various cellular RNAs revealed distinct pathways by which they could be recruited during retrotransposition; template choice or template switching. Here, by analyzing genomic structures and retrotransposition signatures associated with small nuclear RNA (snRNA) sequences, we identified distinct recruiting steps during the L1 retrotransposition cycle for the formation of snRNA-processed pseudogenes. Interestingly, some of the identified recruiting steps take place in the nucleus. Moreover, after comparison to other vertebrate genomes, we established that snRNA amplification by template switching is common to many LINE families from several LINE clades. Finally, we suggest that U6 snRNA copies can serve as markers of L1 retrotransposition dynamics in mammalian genomes.

Ancient traces of tailless retropseudogenes in therian genomes

Transposable elements, once described by Barbara McClintock as controlling genetic units, not only occupy the largest part of our genome but are also a prominent moving force of genomic plasticity and innovation. They usually replicate and reintegrate into genomes silently, sometimes causing malfunctions or misregulations, but occasionally millions of years later, a few may evolve into new functional units. Retrotransposons make their way into the genome following reverse transcription of RNA molecules and chromosomal insertion. In therian mammals, long interspersed elements 1 (LINE1s) self-propagate but also coretropose many RNAs, including mRNAs and small RNAs that usually exhibit an oligo(A) tail. The revitalization of specific LINE1 elements in the mammalian lineage about 150 Ma parallels the rise of many other nonautonomous mobilized genomic elements. We previously identified and described hundreds of tRNA-derived retropseudogenes missing characteristic oligo(A) tails consequently termed tailless retropseudogenes. Additional analyses now revealed hundreds of thousands of tailless retropseudogenes derived from nearly all types of RNAs. We extracted 2,402 perfect tailless sequences (with discernible flanking target site duplications) originating from tRNAs, spliceosomal RNAs, 5S rRNAs, 7SK RNAs, mRNAs, and others. Interestingly, all are truncated at one or more defined positions that coincide with internal single-stranded regions. 5S ribosomal and U2 spliceosomal RNAs were analyzed in the context of mammalian phylogeny to discern the origin of the therian LINE1 retropositional system that evolved in our 150-Myr-old ancestor.

Identification of human-specific AluS elements through comparative genomics

Mobile elements are responsible for ~45% of the human genome. Among them is the Alu element, accounting for 10% of the human genome (>1.1million copies). Several studies of Alu elements have reported that they are frequently involved in human genetic diseases and genomic rearrangements. In this study, we investigated the AluS subfamily, which is a relatively old Alu subfamily and has the highest copy number in primate genomes. Previously, a set of 263 human-specific AluS insertions was identified in the human genome. To validate these, we compared each of the human-specific AluS loci with its pre-insertion site in other primate genomes, including chimpanzee, gorilla, and orangutan. We obtained 24 putative human-specific AluS candidates via the in silico analysis and manual inspection, and then tried to verify them using PCR amplification and DNA sequencing. Through the PCR product sequencing, we were able to detect two instances of near-parallel Alu insertions in nearby sites that led to computational false negatives. Finally, we computationally and experimentally verified 23 human-specific AluS elements. We reported three alternative Alu insertion events, which are accompanied by filler DNA and/or Alu retrotransposition mediated-deletion. Bisulfite sequencing was carried out to examine DNA methylation levels of human-specific AluS elements. The results showed that fixed AluS elements are hypermethylated compared with polymorphic elements, indicating a possible relation between DNA methylation and Alu fixation in the human genome.

An Intronic SINE insertion in FAM161A that causes exon-skipping is associated with progressive retinal atrophy in Tibetan Spaniels and Tibetan Terriers

Progressive retinal atrophy (PRA) in dogs is characterised by the degeneration of the photoreceptor cells of the retina, resulting in vision loss and eventually complete blindness. The condition affects more than 100 dog breeds and is known to be genetically heterogeneous between breeds. Around 19 mutations have now been identified that are associated with PRA in around 49 breeds, but for the majority of breeds the mutation(s) responsible have yet to be identified. Using genome-wide association with 22 Tibetan Spaniel PRA cases and 10 controls, we identified a novel PRA locus, PRA3, on CFA10 (p_raw = 2.01×10⁻⁵, p_genome = 0.014), where a 3.8 Mb region was homozygous within 12 cases. Using targeted next generation sequencing, a short interspersed nuclear element insertion was identified near a splice acceptor site in an intron of a provocative gene, FAM161A. Analysis of mRNA from an affected dog revealed that the SINE causes exon skipping, resulting in a frame shift, leading to a downstream premature termination codon and possibly a truncated protein product. This mutation segregates with the disease in 22 out of 35 cases tested (63%). Of the PRA controls, none are homozygous for the mutation, 15% carry the mutation and 85% are homozygous wildtype. This mutation was also identified in Tibetan Terriers, although our results indicate that PRA is genetically heterogeneous in both Tibetan Spaniels and Tibetan Terriers.

Characterisation of genetic variation in ST8SIA2 and its interaction region in NCAM1 in patients with bipolar disorder

Alpha-2,8-sialyltransferase 2 (ST8SIA2) is an enzyme responsible for the transfer of polysialic acid (PSA) to glycoproteins, principally the neuronal cell adhesion molecule (NCAM1), and is involved in neuronal plasticity. Variants within ST8SIA2 have previously shown association with bipolar disorder, schizophrenia and autism. In addition, altered PSA-NCAM expression in brains of patients with schizophrenia or bipolar disorder indicates a functional dysregulation of glycosylation in mental illness. To explore the role of sequence variation affecting PSA-NCAM formation, we conducted a targeted re-sequencing study of a ∼100 kb region – including the entire ST8SIA2 gene and its region of interaction with NCAM1 – in 48 Caucasian cases with bipolar disorder using the Roche 454 platform. We identified over 400 DNA variants, including 47 putative novel variants not described in dbSNP. Validation of a subset of variants via Sequenom showed high reliability of Roche 454 genotype calls (97% genotype concordance, with 80% of novel variants independently verified). We did not observe major loss-of-function mutations that would affect PSA-NCAM formation, either by ablating ST8SIA2 function or by affecting the ability of NCAM1 to be glycosylated. However, we identified 13 SNPs in the UTRs of ST8SIA2, a synonymous coding SNP in exon 5 (rs2305561, P207P) and many additional non-coding variants that may influence splicing or regulation of ST8SIA2 expression. We calculated nucleotide diversity within ST8SIA2 on specific haplotypes, finding that the diversity on the specific “risk” and “protective” haplotypes was lower than other non-disease-associated haplotypes, suggesting that putative functional variation may have arisen on a spectrum of haplotypes. We have identified common and novel variants (rs11074064, rs722645, 15∶92961050) that exist on a spectrum of haplotypes, yet are plausible candidates for conferring the effect of risk and protective haplotypes via multiple enhancer elements. A Galaxy workflow/pipeline for sequence analysis used herein is available at: https://main.g2.bx.psu.edu/u/a-shaw-neura/p/next-generation-resources.

Alu elements: an intrinsic source of human genome instability

Alu elements are ∼300bp sequences that have amplified via an RNA intermediate leading to the accumulation of over 1 million copies in the human genome. Although a few of the copies are active, Alu germline activity is the highest of all human retrotransposons and does significantly contribute to genetic disease and population diversity. There are two basic mechanisms by which Alu elements contribute to disease: through insertional mutagenesis and as a large source of repetitive sequences that contribute to nonallelic homologous recombination (NAHR) that cause genetic deletions and duplications.

Analysis of western lowland gorilla (Gorilla gorilla gorilla) specific Alu repeats

Background

Research into great ape genomes has revealed widely divergent activity levels over time for Alu elements. However, the diversity of this mobile element family in the genome of the western lowland gorilla has previously been uncharacterized. Alu elements are primate-specific short interspersed elements that have been used as phylogenetic and population genetic markers for more than two decades. Alu elements are present at high copy number in the genomes of all primates surveyed thus far. The AluY subfamily and its derivatives have been recognized as the evolutionarily youngest Alu subfamily in the Old World primate lineage.

Results

Here we use a combination of computational and wet-bench laboratory methods to assess and catalog AluY subfamily activity level and composition in the western lowland gorilla genome (gorGor3.1). A total of 1,075 independent AluY insertions were identified and computationally divided into 10 subfamilies, with the largest number of gorilla-specific elements assigned to the canonical AluY subfamily.

Conclusions

The retrotransposition activity level appears to be significantly lower than that seen in the human and chimpanzee lineages, while higher than that seen in orangutan genomes, indicative of differential Alu amplification in the western lowland gorilla lineage as compared to other Homininae.

RNA-Mediated Gene Duplication and Retroposons: Retrogenes, LINEs, SINEs, and Sequence Specificity

A substantial number of “retrogenes” that are derived from the mRNA of various intron-containing genes have been reported. A class of mammalian retroposons, long interspersed element-1 (LINE1, L1), has been shown to be involved in the reverse transcription of retrogenes (or processed pseudogenes) and non-autonomous short interspersed elements (SINEs). The 3′-end sequences of various SINEs originated from a corresponding LINE. As the 3′-untranslated regions of several LINEs are essential for retroposition, these LINEs presumably require “stringent” recognition of the 3′-end sequence of the RNA template. However, the 3′-ends of mammalian L1s do not exhibit any similarity to SINEs, except for the presence of 3′-poly(A) repeats. Since the 3′-poly(A) repeats of L1 and Alu SINE are critical for their retroposition, L1 probably recognizes the poly(A) repeats, thereby mobilizing not only Alu SINE but also cytosolic mRNA. Many flowering plants only harbor L1-clade LINEs and a significant number of SINEs with poly(A) repeats, but no homology to the LINEs. Moreover, processed pseudogenes have also been found in flowering plants. I propose that the ancestral L1-clade LINE in the common ancestor of green plants may have recognized a specific RNA template, with stringent recognition then becoming relaxed during the course of plant evolution.

Mobile DNA distributions refine the phylogeny of "matsutake" mushrooms, Tricholoma sect. Caligata

"Matsutake" mushrooms are formed by several species of Tricholoma sect. Caligata distributed across the northern hemisphere. A phylogenetic analysis of matsutake based on virtually neutral mutations in DNA sequences resolved robust relationships among Tricholoma anatolicum, Tricholoma bakamatsutake, Tricholoma magnivelare, Tricholoma matsutake, and Tricholoma sp. from Mexico (=Tricholoma sp. Mex). However, relationships among these matsutake and other species, such as Tricholoma caligatum and Tricholoma fulvocastaneum, were ambiguous. We, therefore, analyzed genomic copy numbers of σ marY1 , marY1, and marY2N retrotransposons by comparing them with the single-copy mobile DNA megB1 using real-time polymerase chain reaction (PCR) to clarify matsutake phylogeny. We also examined types of megB1-associated domains, composed of a number of poly (A) and poly (T) reminiscent of RNA-derived DNA elements among these species. Both datasets resolved two distinct groups, one composed of T. bakamatsutake, T. fulvocastaneum, and T. caligatum that could have diverged earlier and the other comprising T. magnivelare, Tricholoma sp. Mex, T. anatolicum, and T. matsutake that could have evolved later. In the first group, T. caligatum was the closest to the second group, followed by T. fulvocastaneum and T. bakamatsutake. Within the second group, T. magnivelare was clearly differentiated from the other species. The data suggest that matsutake underwent substantial evolution between the first group, mostly composed of Fagaceae symbionts, and the second group, comprised only of Pinaceae symbionts, but diverged little within each groups. Mobile DNA markers could be useful in resolving difficult phylogenies due to, for example, closely spaced speciation events.

Non-LTR retrotransposons and microsatellites: Partners in genomic variation

The human genome is laden with both non-LTR (long-terminal repeat) retrotransposons and microsatellite repeats. Both types of sequences are able to, either actively or passively, mutagenize the genomes of human individuals and are therefore poised to dynamically alter the human genomic landscape across generations. Non-LTR retrotransposons, such as L1 and Alu, are a major source of new microsatellites, which are born both concurrently and subsequently to L1 and Alu integration into the genome. Likewise, the mutation dynamics of microsatellite repeats have a direct impact on the fitness of their non-LTR retrotransposon parent owing to microsatellite expansion and contraction. This review explores the interactions and dynamics between non-LTR retrotransposons and microsatellites in the context of genomic variation and evolution.

LINE-1-derived poly(A) microsatellites undergo rapid shortening and create somatic and germline mosaicism in mice

Interspersed and tandem repeat sequences comprise the bulk of mammalian genomes. Interspersed repeats result from successive replication by transposable elements, such as Alu and long interspersed element type 1 (L1). Microsatellites are tandem repeats of 1-6 base pairs, among which poly(A) microsatellites are the most abundant in the human genome. The rise and fall of a microsatellite has been depicted as a life cycle. Previous studies have demonstrated that Alu and L1 insertions are a major source of A-rich microsatellites owing to the concurrent formation of a poly(A) DNA tract at the 3'-end of each insertion. The fate of such poly(A) tracts has been studied by surveying the length distribution of genomic resident Alu and L1 insertions. However, these cross-sectional studies provide no information about the tempo of mutation immediately after birth. In this study, de novo L1 insertions were created using a transgenic L1 mouse model and traced through generations to investigate the early life of poly(A) microsatellites. High frequencies of intra-individual and intergenerational shortening were observed for long poly(A) tracts, creating somatic and germline mosaicism at the insertion site, whereas little variation was observed for short poly(A) alleles. As poly(A) microsatellites are the major intrinsic signal for nucleosome positioning, their remarkable abundance and variability make them a significant source of epigenetic variation. Thus, the birth of poly(A) microsatellites from retrotransposons and the subsequent rapid and variable shortening represent a new way with which retrotransposons can modify the genetic and epigenetic architecture of our genome.

Direct mapping of symbolic DNA sequence into frequency domain in global repeat map algorithm

The main feature of global repeat map (GRM) algorithm (www.hazu.hr/grm/software/win/grm2012.exe) is its ability to identify a broad variety of repeats of unbounded length that can be arbitrarily distant in sequences as large as human chromosomes. The efficacy is due to the use of complete set of a K-string ensemble which enables a new method of direct mapping of symbolic DNA sequence into frequency domain, with straightforward identification of repeats as peaks in GRM diagram. In this way, we obtain very fast, efficient and highly automatized repeat finding tool. The method is robust to substitutions and insertions/deletions, as well as to various complexities of the sequence pattern. We present several case studies of GRM use, in order to illustrate its capabilities: identification of α-satellite tandem repeats and higher order repeats (HORs), identification of Alu dispersed repeats and of Alu tandems, identification of Period 3 pattern in exons, implementation of ‘magnifying glass’ effect, identification of complex HOR pattern, identification of inter-tandem transitional dispersed repeat sequences and identification of long segmental duplications. GRM algorithm is convenient for use, in particular, in cases of large repeat units, of highly mutated and/or complex repeats, and of global repeat maps for large genomic sequences (chromosomes and genomes).

Poly(A) binding protein C1 is essential for efficient L1 retrotransposition and affects L1 RNP formation

Poly(A) binding proteins (PABPs) specifically bind the polyadenosine tail of mRNA and have been shown to be important for RNA polyadenylation, translation initiation, and mRNA stability. Using a modified L1 retrotransposition vector, we examined the effects of two PABPs (encoded by PABPN1 and PABPC1) on the retrotransposition activity of the L1 non-long-terminal-repeat (non-LTR) retrotransposon in both HeLa and HEK293T cells. We demonstrated that knockdown of these two genes by RNA interference (RNAi) effectively reduced L1 retrotransposition by 70 to 80% without significantly changing L1 transcription or translation or the status of the poly(A) tail. We identified that both poly(A) binding proteins were associated with the L1 ribonucleoprotein complex, presumably through L1 mRNA. Depletion of PABPC1 caused a defect in L1 RNP formation. Knockdown of the PABPC1 inhibitor PAIP2 increased L1 retrotransposition up to 2-fold. Low levels of exogenous overexpression of PABPN1 and PABPC1 increased L1 retrotransposition, whereas unregulated overexpression of these two proteins caused pleiotropic effects, such as hypersensitivity to puromycin and decreased L1 activity. Our data suggest that PABPC1 is essential for the formation of L1 RNA-protein complexes and may play a role in L1 RNP translocation in the host cell.

Rescuing Alu: recovery of new inserts shows LINE-1 preserves Alu activity through A-tail expansion

Alu elements are trans-mobilized by the autonomous non-LTR retroelement, LINE-1 (L1). Alu-induced insertion mutagenesis contributes to about 0.1% human genetic disease and is responsible for the majority of the documented instances of human retroelement insertion-induced disease. Here we introduce a SINE recovery method that provides a complementary approach for comprehensive analysis of the impact and biological mechanisms of Alu retrotransposition. Using this approach, we recovered 226 de novo tagged Alu inserts in HeLa cells. Our analysis reveals that in human cells marked Alu inserts driven by either exogenously supplied full length L1 or ORF2 protein are indistinguishable. Four percent of de novo Alu inserts were associated with genomic deletions and rearrangements and lacked the hallmarks of retrotransposition. In contrast to L1 inserts, 5′ truncations of Alu inserts are rare, as most of the recovered inserts (96.5%) are full length. De novo Alus show a random pattern of insertion across chromosomes, but further characterization revealed an Alu insertion bias exists favoring insertion near other SINEs, highly conserved elements, with almost 60% landing within genes. De novo Alu inserts show no evidence of RNA editing. Priming for reverse transcription rarely occurred within the first 20 bp (most 5′) of the A-tail. The A-tails of recovered inserts show significant expansion, with many at least doubling in length. Sequence manipulation of the construct led to the demonstration that the A-tail expansion likely occurs during insertion due to slippage by the L1 ORF2 protein. We postulate that the A-tail expansion directly impacts Alu evolution by reintroducing new active source elements to counteract the natural loss of active Alus and minimizing Alu extinction.

SINEs are mobile elements that are found ubiquitously throughout a large diversity of genomes from plants to mammals. The human SINE, Alu, is among the most successful mobile elements, with more than one million copies in the genome. Due to its high activity and ability to insert throughout the genome, Alu retrotransposition is responsible for the majority of diseases reported to be caused by mobile element activity. To further evaluate the genomic impact of SINEs, we recovered and characterized over 200 de novo Alu inserts under controlled conditions. Our data reinforce observations on the mutagenic potential of Alu, with newly retrotransposed Alu elements favoring insertion into genic and highly conserved elements. Alu-mediated deletions and rearrangements are infrequent and lack the typical hallmarks of TPRT retrotransposition, suggesting the use of an alternate method for resolving retrotransposition intermediates or an atypical insertion mechanism. Our data also provide novel insights into SINE retrotransposition biology. We found that slippage of L1 ORF2 protein during reverse transcription expands the A-tails of de novo insertions. We propose that the L1 ORF2 protein plays a major role in minimizing Alu extinction by reintroducing active Alu elements to counter the natural loss of Alu source elements.

Alu expression in human cell lines and their retrotranspositional potential

BACKGROUND

The vast majority of the 1.1 million Alu elements are retrotranspositionally inactive, where only a few loci referred to as 'source elements' can generate new Alu insertions. The first step in identifying the active Alu sources is to determine the loci transcribed by RNA polymerase III (pol III). Previous genome-wide analyses from normal and transformed cell lines identified multiple Alu loci occupied by pol III factors, making them candidate source elements.

FINDINGS

Analysis of the data from these genome-wide studies determined that the majority of pol III-bound Alus belonged to the older subfamilies Alu S and Alu J, which varied between cell lines from 62.5% to 98.7% of the identified loci. The pol III-bound Alus were further scored for estimated retrotransposition potential (ERP) based on the absence or presence of selected sequence features associated with Alu retrotransposition capability. Our analyses indicate that most of the pol III-bound Alu loci candidates identified lack the sequence characteristics important for retrotransposition.

CONCLUSIONS

These data suggest that Alu expression likely varies by cell type, growth conditions and transformation state. This variation could extend to where the same cell lines in different laboratories present different Alu expression patterns. The vast majority of Alu loci potentially transcribed by RNA pol III lack important sequence features for retrotransposition and the majority of potentially active Alu loci in the genome (scored high ERP) belong to young Alu subfamilies. Our observations suggest that in an in vivo scenario, the contribution of Alu activity on somatic genetic damage may significantly vary between individuals and tissues.