Alu element mutation spectra: molecular clocks and the effect of DNA methylation

Evolution of retrocopies in the context of HUSH silencing

Retrotransposition is one of the main factors responsible for gene duplication and thus genome evolution. However, the sequences that undergo this process are not only an excellent source of biological diversity, but in certain cases also pose a threat to the integrity of the DNA. One of the mechanisms that protects against the incorporation of mobile elements is the HUSH complex, which is responsible for silencing long, intronless, transcriptionally active transposed sequences that are rich in adenine on the sense strand. In this study, broad sets of human and porcine retrocopies were analysed with respect to the above factors, taking into account evolution of these molecules. Analysis of expression pattern, genomic structure, transcript length, and nucleotide substitution frequency showed the strong relationship between the expression level and exon length as well as the protective nature of introns. The results of the studies also showed that there is no direct correlation between the expression level and adenine content. However, protein-coding retrocopies, which have a lower adenine content, have a significantly higher expression level than the adenine-rich non-coding but expressed retrocopies. Therefore, although the mechanism of HUSH silencing may be an important part of the regulation of retrocopy expression, it is one component of a more complex molecular network that remains to be elucidated.

Alu RNA fold links splicing with signal recognition particle proteins

Transcriptomic diversity in primates was considerably expanded by exonizations of intronic Alu elements. To better understand their cellular mechanisms we have used structure-based mutagenesis coupled with functional and proteomic assays to study the impact of successive primate mutations and their combinations on inclusion of a sense-oriented AluJ exon in the human F8 gene. We show that the splicing outcome was better predicted by consecutive RNA conformation changes than by computationally derived splicing regulatory motifs. We also demonstrate an involvement of SRP9/14 (signal recognition particle) heterodimer in splicing regulation of Alu-derived exons. Nucleotide substitutions that accumulated during primate evolution relaxed the conserved left-arm AluJ structure including helix H1 and reduced the capacity of SRP9/14 to stabilize the closed Alu conformation. RNA secondary structure-constrained mutations that promoted open Y-shaped conformations of the Alu made the Alu exon inclusion reliant on DHX9. Finally, we identified additional SRP9/14 sensitive Alu exons and predicted their functional roles in the cell. Together, these results provide unique insights into architectural elements required for sense Alu exonization, identify conserved pre-mRNA structures involved in exon selection and point to a possible chaperone activity of SRP9/14 outside the mammalian signal recognition particle.

Alu-minating the Mechanisms Underlying Primate Cortex Evolution

The higher-order cognitive functions observed in primates correlate with the evolutionary enhancement of cortical volume and folding, which in turn are driven by the primate-specific expansion of cellular diversity in the developing cortex. Underlying these changes is the diversification of molecular features including the creation of human and/or primate-specific genes, the activation of specific molecular pathways, and the interplay of diverse layers of gene regulation. We review and discuss evidence for connections between Alu elements and primate brain evolution, the evolutionary milestones of which are known to coincide along primate lineages. Alus are repetitive elements that contribute extensively to the acquisition of novel genes and the expansion of diverse gene regulatory layers, including enhancers, alternative splicing, RNA editing, and microRNA pathways. By reviewing the impact of Alus on molecular features linked to cortical expansions or gyrification or implications in cognitive deficits, we suggest that future research focusing on the role of Alu-derived molecular events in the context of brain development may greatly advance our understanding of higher-order cognitive functions and neurologic disorders.

Can-SINE dynamics in the giant panda and three other Caniformia genomes

Background

Although repeat sequences constitute about 37% of carnivore genomes, the characteristics and distribution of repeat sequences among carnivore genomes have not been fully investigated. Based on the updated Repbase library, we re-annotated transposable elements (TEs) in four Caniformia genomes (giant panda, polar bear, domestic dog, and domestic ferret) and performed a systematic, genome-wide comparison focusing on the Carnivora-specific SINE family, Can-SINEs.

Results

We found the majority of young recently integrated transposable elements are LINEs and SINEs in carnivore genomes. In particular, SINEC1_AMe, SINEC1B_AMe and SINEC_C1 are the top three most abundant Can-SINE subfamilies in the panda and polar bear genomes. Transposition in transposition analysis indicates that SINEC1_AMe and SINEC1B_AMe are the most active subfamilies in the panda and the polar bear genomes. SINEC2A1_CF and SINEC1A_CF subfamilies show a higher retrotransposition activity in the dog genome, and MVB2 subfamily is the most active Can-SINE in the ferret genome. As the giant panda is an endangered icon species, we then focused on the identification of panda specific Can-SINEs. With the panda-associated two-way genome alignments, we identified 250 putative panda-specific (PPS) elements (139 SINEC1_AMes and 111 SINEC1B_AMes) that inserted in the panda genome but were absent at the orthologous regions of the other three genomes. Further investigation of these PPS elements allowed us to identify a new Can-SINE subfamily, the SINEC1_AMe2, which was distinguishable from the current SINEC1_AMe consensus by four non-CpG sites. SINEC1_AMe2 has a high copy number (> 100,000) in the panda and polar bear genomes and the vast majority (> 96%) of the SINEC1_AMe2 elements have divergence rates less than 10% in both genomes.

Conclusions

Our results suggest that Can-SINEs show lineage-specific retransposition activity in the four genomes and have an important impact on the genomic landscape of different Caniformia lineages. Combining these observations with results from the COSEG, Network, and target site duplication analysis, we suggest that SINEC1_AMe2 is a young mobile element subfamily and currently active in both the panda and polar bear genomes.

Intronic antisense Alu elements have a negative splicing effect on the inclusion of adjacent downstream exons

Alu elements occupy 10% of the human genome. However, although they contribute to genomic and transcriptomic diversity, their function is still not fully understood. We hypothesized that intronic Alu elements may contribute to alternative splicing. We therefore examined their effect on splicing using minigene constructs including exon 9-exon 11 inclusive of ACAT1 with truncated introns 9 and 10. These constructs contained a suboptimal splice acceptor site for intron 9. Insertion of AluY-partial AluSz6-AluSx, originally located in ACAT1 intron 5, in an antisense direction within intron 9 had a negative effect on exon 10 inclusion. This effect was additive with that of an exonic splicing enhancer mutation in exon 10, and was canceled by the substitution of G for C at the first nucleotide of exon 10 which optimized the splice acceptor site of intron 9. A sense AluY-partial AluSz6-AluSx insertion had no effect on exon 10 inclusion, and one antisense AluSx insertion had a similar effect to antisense AluY-partial AluSz6-AluSx insertion. The shorter the distance between the antisense Alu element and exon 10, the greater the negative effect on exon 10 inclusion. This distance effect was more evident for suboptimal than optimal splice sites. Based on our data, we propose that intronic antisense Alu elements contribute to alternative splicing and transcriptomic diversity in some genes, especially when splice acceptor sites are suboptimal.

Integrated genome-wide Alu methylation and transcriptome profiling analyses reveal novel epigenetic regulatory networks associated with autism spectrum disorder

Background

Alu elements are a group of repetitive elements that can influence gene expression through CpG residues and transcription factor binding. Altered gene expression and methylation profiles have been reported in various tissues and cell lines from individuals with autism spectrum disorder (ASD). However, the role of Alu elements in ASD remains unclear. We thus investigated whether Alu elements are associated with altered gene expression profiles in ASD.

Methods

We obtained five blood-based gene expression profiles from the Gene Expression Omnibus database and human Alu-inserted gene lists from the TranspoGene database. Differentially expressed genes (DEGs) in ASD were identified from each study and overlapped with the human Alu-inserted genes. The biological functions and networks of Alu-inserted DEGs were then predicted by Ingenuity Pathway Analysis (IPA). A combined bisulfite restriction analysis of lymphoblastoid cell lines (LCLs) derived from 36 ASD and 20 sex- and age-matched unaffected individuals was performed to assess the global DNA methylation levels within Alu elements, and the Alu expression levels were determined by quantitative RT-PCR.

Results

In ASD blood or blood-derived cells, 320 Alu-inserted genes were reproducibly differentially expressed. Biological function and pathway analysis showed that these genes were significantly associated with neurodevelopmental disorders and neurological functions involved in ASD etiology. Interestingly, estrogen receptor and androgen signaling pathways implicated in the sex bias of ASD, as well as IL-6 signaling and neuroinflammation signaling pathways, were also highlighted. Alu methylation was not significantly different between the ASD and sex- and age-matched control groups. However, significantly altered Alu methylation patterns were observed in ASD cases sub-grouped based on Autism Diagnostic Interview-Revised scores compared with matched controls. Quantitative RT-PCR analysis of Alu expression also showed significant differences between ASD subgroups. Interestingly, Alu expression was correlated with methylation status in one phenotypic ASD subgroup.

Conclusion

Alu methylation and expression were altered in LCLs from ASD subgroups. Our findings highlight the association of Alu elements with gene dysregulation in ASD blood samples and warrant further investigation. Moreover, the classification of ASD individuals into subgroups based on phenotypes may be beneficial and could provide insights into the still unknown etiology and the underlying mechanisms of ASD.

Electronic supplementary material

The online version of this article (10.1186/s13229-018-0213-9) contains supplementary material, which is available to authorized users.

Alu exaptation enriches the human transcriptome by introducing new gene ends

In mammals, transposable elements are largely silenced, but under fortuitous circumstances may be co-opted to play a functional role. Here, we show that when Alu elements are inserted within or nearby genes in sense orientation, they may contribute to the transcriptome diversity by forming new cleavage and polyadenylation sites. We mapped these new gene ends in human onto the Alu sequence and identified three hotspots of cleavage and polyadenylation site formation. Interestingly, the native Alu sequence does not contain any canonical polyadenylation signal. We therefore studied what evolutionary processes might explain the formation of these specific hotspots of novel gene ends. We show that two of the three hotspots might have emerged from mutational processes that turned sequences that resemble polyadenylation signals into full-blown canonical signals, whereas one hotspot is tightly linked to the process of Alu insertion into the genome. Overall, Alu elements may lie behind the formation of 302 new gene end variants, affecting a total of 243 genes. Intergenic Alu elements may elongate genes by creating a downstream cleavage site, intronic Alu elements may lead to gene variants which code for truncated proteins, and 3'UTR Alu elements may result in gene variants with alternative 3'UTR.

3'UTR-located ALU Elements: Donors of Potetial miRNA Target Sites and Mediators of Network miRNA-based Regulatory Interactions

Recent research data reveal complex, network-based interactions between mobile elements and regulatory systems of eukaryotic cells. In this article, we focus on regulatory interactions between Alu elements and micro RNAs (miRNAs). Our results show that the majority of the Alu sequences inserted in 3'UTRs of analyzed human genes carry strong potential target sites for at least 53 different miRNAs. Thus, 3'UTR-located Alu elements may play the role of mobile regulatory modules that supply binding sites for miRNA regulation. Their abundance and ability to distribute a set of certain miRNA target sites may have an important role in establishment, extension, network organization, and, as we suppose – in the regulation and environment-dependent activation/inactivation of some elements of the miRNA regulatory system, as well as for a larger scale RNA-based regulatory interactions. The Alu-miRNA connection may be crucial especially for the primate/human evolution.

Cross platform analysis of transcriptomic data identifies ageing has distinct and opposite effects on tendon in males and females

The development of tendinopathy is influenced by a variety of factors including age, gender, sex hormones and diabetes status. Cross platform comparative analysis of transcriptomic data elucidated the connections between these entities in the context of ageing. Tissue-engineered tendons differentiated from bone marrow derived mesenchymal stem cells from young (20-24 years) and old (54-70 years) donors were assayed using ribonucleic acid sequencing (RNA-seq). Extension of the experiment to microarray and RNA-seq data from tendon identified gender specific gene expression changes highlighting disparity with existing literature and published pathways. Separation of RNA-seq data by sex revealed underlying negative binomial distributions which increased statistical power. Sex specific de novo transcriptome assemblies generated fewer larger transcripts that contained miRNAs, lincRNAs and snoRNAs. The results identify that in old males decreased expression of CRABP2 leads to cell proliferation, whereas in old females it leads to cellular senescence. In conjunction with existing literature the results explain gender disparity in the development and types of degenerative diseases as well as highlighting a wide range of considerations for the analysis of transcriptomic data. Wider implications are that degenerative diseases may need to be treated differently in males and females because alternative mechanisms may be involved.

Romain CP, Gonzalez-Quiroga P, Denham MW, Mierl JR, Konkel MK, Batzer MA. Evolution of Alu Subfamily Structure in the Saimiri Lineage of New World Monkeys

Squirrel monkeys, Saimiri, are commonly found in zoological parks and used in biomedical research. S. boliviensis is the most common species for research; however, there is little information about genome evolution within this primate lineage. Here, we reconstruct the Alu element sequence amplification and evolution in the genus Saimiri at the time of divergence within the family Cebidae lineage. Alu elements are the most successful SINE (Short Interspersed Element) in primates. Here, we report 46 Saimiri lineage specific Alu subfamilies. Retrotransposition activity involved subfamilies related to AluS, AluTa10, and AluTa15. Many subfamilies are simultaneously active within the Saimiri lineage, a finding which supports the stealth model of Alu amplification. We also report a high resolution analysis of Alu subfamilies within the S. boliviensis genome [saiBol1].

Oxidative stress and repetitive element methylation changes in artisanal gold miners occupationally exposed to mercury

Mercury (Hg) exposure is a public health concern due to its persistence in the environment and its high toxicity. Such toxicity has been associated with the generation of oxidative stress in occupationally exposed subjects, such as artisanal gold miners. In this study, we characterize occupational exposure to Hg by measuring blood, urine and hair levels, and investigate oxidative stress and DNA methylation associated with gold mining. To do this, samples from 53 miners and 36 controls were assessed. We show higher levels of oxidative stress marker 8-OHdG in the miners. Differences in LINE1 and Alu(Yb8) DNA methylation between gold miners and control group are present in peripheral blood leukocytes. LINE1 methylation is positively correlated with 8-OHdG levels, while XRCC1 and LINE1 methylation are positively correlated with Hg levels. These results suggest an effect of Hg on oxidative stress and DNA methylation in gold miners that may have an impact on miners' health.

Discovery of rare, diagnostic AluYb8/9 elements in diverse human populations

BACKGROUND

Polymorphic human Alu elements are excellent tools for assessing population structure, and new retrotransposition events can contribute to disease. Next-generation sequencing has greatly increased the potential to discover Alu elements in human populations, and various sequencing and bioinformatics methods have been designed to tackle the problem of detecting these highly repetitive elements. However, current techniques for Alu discovery may miss rare, polymorphic Alu elements. Combining multiple discovery approaches may provide a better profile of the polymorphic Alu mobilome. AluYb8/9 elements have been a focus of our recent studies as they are young subfamilies (~2.3 million years old) that contribute ~30% of recent polymorphic Alu retrotransposition events. Here, we update our ME-Scan methods for detecting Alu elements and apply these methods to discover new insertions in a large set of individuals with diverse ancestral backgrounds.

RESULTS

We identified 5,288 putative Alu insertion events, including several hundred novel AluYb8/9 elements from 213 individuals from 18 diverse human populations. Hundreds of these loci were specific to continental populations, and 23 non-reference population-specific loci were validated by PCR. We provide high-quality sequence information for 68 rare AluYb8/9 elements, of which 11 have hallmarks of an active source element. Our subfamily distribution of rare AluYb8/9 elements is consistent with previous datasets, and may be representative of rare loci. We also find that while ME-Scan and low-coverage, whole-genome sequencing (WGS) detect different Alu elements in 41 1000 Genomes individuals, the two methods yield similar population structure results.

CONCLUSION

Current in-silico methods for Alu discovery may miss rare, polymorphic Alu elements. Therefore, using multiple techniques can provide a more accurate profile of Alu elements in individuals and populations. We improved our false-negative rate as an indicator of sample quality for future ME-Scan experiments. In conclusion, we demonstrate that ME-Scan is a good supplement for next-generation sequencing methods and is well-suited for population-level analyses.

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.

Neurotoxic Doses of Chronic Methamphetamine Trigger Retrotransposition of the Identifier Element in Rat Dorsal Dentate Gyrus

Short interspersed elements (SINEs) are typically silenced by DNA hypermethylation in somatic cells, but can retrotranspose in proliferating cells during adult neurogenesis. Hypomethylation caused by disease pathology or genotoxic stress leads to genomic instability of SINEs. The goal of the present investigation was to determine whether neurotoxic doses of binge or chronic methamphetamine (METH) trigger retrotransposition of the identifier (ID) element, a member of the rat SINE family, in the dentate gyrus genomic DNA. Adult male Sprague-Dawley rats were treated with saline or high doses of binge or chronic METH and sacrificed at three different time points thereafter. DNA methylation analysis, immunohistochemistry and next-generation sequencing (NGS) were performed on the dorsal dentate gyrus samples. Binge METH triggered hypomethylation, while chronic METH triggered hypermethylation of the CpG-2 site. Both METH regimens were associated with increased intensities in poly(A)-binding protein 1 (PABP1, a SINE regulatory protein)-like immunohistochemical staining in the dentate gyrus. The amplification of several ID element sequences was significantly higher in the chronic METH group than in the control group a week after METH, and they mapped to genes coding for proteins regulating cell growth and proliferation, transcription, protein function as well as for a variety of transporters. The results suggest that chronic METH induces ID element retrotransposition in the dorsal dentate gyrus and may affect hippocampal neurogenesis.

Evolutionary dynamics of selfish DNA explains the abundance distribution of genomic subsequences

Since the sequencing of large genomes, many statistical features of their sequences have been found. One intriguing feature is that certain subsequences are much more abundant than others. In fact, abundances of subsequences of a given length are distributed with a scale-free power-law tail, resembling properties of human texts, such as Zipf's law. Despite recent efforts, the understanding of this phenomenon is still lacking. Here we find that selfish DNA elements, such as those belonging to the Alu family of repeats, dominate the power-law tail. Interestingly, for the Alu elements the power-law exponent increases with the length of the considered subsequences. Motivated by these observations, we develop a model of selfish DNA expansion. The predictions of this model qualitatively and quantitatively agree with the empirical observations. This allows us to estimate parameters for the process of selfish DNA spreading in a genome during its evolution. The obtained results shed light on how evolution of selfish DNA elements shapes non-trivial statistical properties of genomes.

Pinpointing the vesper bat transposon revolution using the Miniopterus natalensis genome

BACKGROUND

Around 40 million years ago DNA transposons began accumulating in an ancestor of bats in the family Vespertilionidae. Since that time, Class II transposons have been continuously reinvading and accumulating in vespertilionid genomes at a rate that is unprecedented in mammals. Miniopterus (Miniopteridae), a genus of long-fingered bats that was recently elevated from Vespertilionidae, is the sister taxon to the vespertilionids and is often used as an outgroup when studying transposable elements in vesper bats. Previous wet-lab techniques failed to identify Helitrons, TcMariners, or hAT transposons in Miniopterus. Limitations of those methods and ambiguous results regarding the distribution of piggyBac transposons left some questions as to the distribution of Class II elements in this group. The recent release of the Miniopterus natalensis genome allows for transposable element discovery with a higher degree of precision.

RESULTS

Here we analyze the transposable element content of M. natalensis to pinpoint with greater accuracy the taxonomic distribution of Class II transposable elements in bats. These efforts demonstrate that, compared to the vespertilionids, Class II TEs are highly mutated and comprise only a small portion of the M. natalensis genome. Despite the limited Class II content, M. natalensis possesses a limited number of lineage-specific, low copy number piggyBacs and shares several TcMariner families with vespertilionid bats. Multiple efforts to identify Helitrons, one of the major TE components of vesper bat genomes, using de novo repeat identification and structural based searches failed.

CONCLUSIONS

These observations combined with previous results inform our understanding of the events leading to the unique Class II element acquisition that characterizes vespertilionids. While it appears that a small number of TcMariner and piggyBac elements were deposited in the ancestral Miniopterus + vespertilionid genome, these elements are not present in M. natalensis genome at high copy number. Instead, this work indicates that the vesper bats alone experienced the expansion of TEs ranging from Helitrons to piggyBacs to hATs.

Detection of differential DNA methylation in repetitive DNA of mice and humans perinatally exposed to bisphenol A

Developmental exposure to bisphenol A (BPA) has been shown to induce changes in DNA methylation in both mouse and human genic regions; however, the response in repetitive elements and transposons has not been explored. Here we present novel methodology to combine genomic DNA enrichment with RepeatMasker analysis on next-generation sequencing data to determine the effect of perinatal BPA exposure on repetitive DNA at the class, family, subfamily, and individual insertion level in both mouse and human samples. Mice were treated during gestation and lactation to BPA in chow at 0, 50, or 50,000 ng/g levels and total BPA was measured in stratified human fetal liver tissue samples as low (non-detect to 0.83 ng/g), medium (3.5 to 5.79 ng/g), or high (35.44 to 96.76 ng/g). Transposon methylation changes were evident in human classes, families, and subfamilies, with the medium group exhibiting hypomethylation compared to both high and low BPA groups. Mouse repeat classes, families, and subfamilies did not respond to BPA with significantly detectable differential DNA methylation. In human samples, 1251 individual transposon loci were detected as differentially methylated by BPA exposure, but only 19 were detected in mice. Of note, this approach recapitulated the discovery of a previously known mouse environmentally labile metastable epiallele, Cabp^IAP. Thus, by querying repetitive DNA in both mouse and humans, we report the first known transposons in humans that respond to perinatal BPA exposure.

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.

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.

Drawing a fine line on endogenous retroelement activity

Endogenous retroelements (EREs) are essential motors of evolution yet require careful control to prevent genomic catastrophes, notably during the vulnerable phases of epigenetic reprogramming that occur immediately after fertilization and in germ cells. Accordingly, a variety of mechanisms restrict these mobile genetic units. Previous studies have revealed the importance of KRAB-containing zinc finger proteins (KRAB-ZFPs) and their cofactor, KAP1, in the early embryonic silencing of endogenous retroviruses and so-called SVAs, but the implication of this transcriptional repression system in the control of LINE-1, the only known active autonomous retrotransposon in the human genome, was thought to be marginal. Two recent studies straighten the record by revealing that the KRAB/KAP system is key to the control of L1 in embryonic stem (ES) cells, and go further in demonstrating that DNA methylation and KRAB/KAP1-induced repression contribute to this process in an evolutionally dynamic fashion. These results shed light on the delicate equilibrium between higher vertebrates and endogenous retroelements, which are not just genetic invaders calling for strict control but rather a constantly renewed and nicely exploitable source of evolutionary potential.

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.

Hominoid composite non-LTR retrotransposons-variety, assembly, evolution, and structural determinants of mobilization

SVA (SINE-R-VNTR-Alu) elements constitute the youngest family of composite non-LTR retrotransposons in hominoid primates. The sequence of their assembly, however, remains unclear. Recently, a second family of VNTR-containing composites, LAVA (L1-Alu-VNTR-Alu), has been identified in gibbons. We now report the existence of two additional VNTR composite families, PVA (PTGR2-VNTR-Alu) and FVA (FRAM-VNTR-Alu), in the genome of Nomascus leucogenys. Like LAVA, they share the 5'-Alu-like region and VNTR with SVA, but differ at their 3'-ends. The 3'-end of PVA comprises part of the PTGR2 gene, whereas FVA is characterized by the presence of a partial FRAM element in its 3'-domain. Splicing could be identified as the mechanism of acquisition of the variant 3'-ends in all four families of VNTR composites. SVAs have been shown to be mobilized by the L1 protein machinery in trans. A critical role in this process has been ascribed to their 5'-hexameric repeat/ Alu-like region. The Alu-like region displays specific features in each of the VNTR composite families/subfamilies with characteristic deletions found in the evolutionary younger subfamilies. Using reciprocal exchanges between SVA_E and PVA/FVA elements, we demonstrate that the structure, not the presence of the (CCCTCT)n/ Alu-like region determines mobilization capacity. Combination of LAVA and SVA_E domains does not yield any active elements-suggesting the use of different combinations of host factors for the two major groups of VNTR composites. Finally, we demonstrate that the LAVA 3'-L1ME5 fragment attenuates mobilization capacity.

Evolutionary age of repetitive element subfamilies and sensitivity of DNA methylation to airborne pollutants

Background

Repetitive elements take up >40% of the human genome and can change distribution through transposition, thus generating subfamilies. Repetitive element DNA methylation has associated with several diseases and environmental exposures, including exposure to airborne pollutants. No systematic analysis has yet been conducted to examine the effects of exposures across different repetitive element subfamilies. The purpose of the study is to evaluate sensitivity of DNA methylation in differentially‒evolved LINE, Alu, and HERV subfamilies to different types of airborne pollutants.

Methods

We sampled a total of 120 male participants from three studies (20 high-, 20 low-exposure in each study) of steel workers exposed to metal-rich particulate matter (measured as PM₁₀) (Study 1); gas-station attendants exposed to air benzene (Study 2); and truck drivers exposed to traffic-derived elemental carbon (Study 3). We measured methylation by bisulfite-PCR-pyrosequencing in 10 differentially‒evolved repetitive element subfamilies.

Results

High-exposure groups exhibited subfamily-specific methylation differences compared to low-exposure groups: L1PA2 showed lower DNA methylation in steel workers (P=0.04) and gas station attendants (P=0.03); L1Ta showed lower DNA methylation in steel workers (P=0.02); AluYb8 showed higher DNA methylation in truck drivers (P=0.05). Within each study, dose–response analyses showed subfamily-specific correlations of methylation with exposure levels. Interaction models showed that the effects of the exposures on DNA methylation were dependent on the subfamily evolutionary age, with stronger effects on older LINEs from PM₁₀ (p‒interaction=0.003) and benzene (p‒interaction=0.04), and on younger Alus from PM₁₀ (p-interaction=0.02).

Conclusions

The evolutionary age of repetitive element subfamilies determines differential susceptibility of DNA methylation to airborne pollutants.

Transposable Elements: No More 'Junk DNA'

Since the advent of whole-genome sequencing, transposable elements (TEs), just thought to be 'junk' DNA, have been noticed because of their numerous copies in various eukaryotic genomes. Many studies about TEs have been conducted to discover their functions in their host genomes. Based on the results of those studies, it has been generally accepted that they have a function to cause genomic and genetic variations. However, their infinite functions are not fully elucidated. Through various mechanisms, including de novo TE insertions, TE insertion-mediated deletions, and recombination events, they manipulate their host genomes. In this review, we focus on Alu, L1, human endogenous retrovirus, and short interspersed element/variable number of tandem repeats/Alu (SVA) elements and discuss how they have affected primate genomes, especially the human and chimpanzee genomes, since their divergence.

Orangutan Alu quiescence reveals possible source element: support for ancient backseat drivers

Background

Sequence analysis of the orangutan genome revealed that recent proliferative activity of Alu elements has been uncharacteristically quiescent in the Pongo (orangutan) lineage, compared with all previously studied primate genomes. With relatively few young polymorphic insertions, the genomic landscape of the orangutan seemed like the ideal place to search for a driver, or source element, of Alu retrotransposition.

Results

Here we report the identification of a nearly pristine insertion possessing all the known putative hallmarks of a retrotranspositionally competent Alu element. It is located in an intronic sequence of the DGKB gene on chromosome 7 and is highly conserved in Hominidae (the great apes), but absent from Hylobatidae (gibbon and siamang). We provide evidence for the evolution of a lineage-specific subfamily of this shared Alu insertion in orangutans and possibly the lineage leading to humans. In the orangutan genome, this insertion contains three orangutan-specific diagnostic mutations which are characteristic of the youngest polymorphic Alu subfamily, AluYe5b5_Pongo. In the Homininae lineage (human, chimpanzee and gorilla), this insertion has acquired three different mutations which are also found in a single human-specific Alu insertion.

Conclusions

This seemingly stealth-like amplification, ongoing at a very low rate over millions of years of evolution, suggests that this shared insertion may represent an ancient backseat driver of Alu element expansion.

How does selfing affect the dynamics of selfish transposable elements?

Background

Many theoretical models predicting the dynamics of transposable elements (TEs) in genomes, populations, and species have already been proposed. However, most of them only focus on populations of sexual diploid individuals, and TE dynamics in populations partly composed by autogamous individuals remains poorly investigated. To estimate the impact of selfing on TE dynamics, the short- and long-term evolution of TEs was simulated in outcrossing populations with various proportions of selfing individuals.

Results

Selfing has a deep impact on TE dynamics: the higher the selfing rate, the lower the probability of invasion. Already known non-equilibrium dynamics (complete loss, domestication, cyclical invasion of TEs) can all be described whatever the mating system. However, their pattern and their respective frequencies greatly depend on the selfing rate. For instance, in cyclical dynamics resulting from interactions between autonomous and non-autonomous copies, cycles are faster when the selfing rate increases. Interestingly, an abrupt change in the mating system from sexuality to complete asexuality leads to the loss of all the elements over a few hundred generations. In general, for intermediate selfing rates, the transposition activity remains maintained.

Conclusions

Our theoretical results evidence that a clear and systematic contrast in TE content according to the mating system is expected, with a smooth transition for intermediate selfing rates. Several parameters impact the TE copy number, and all dynamics described in allogamous populations can be also observed in partly autogamous species. This study thus provides new insights to understand the complex signal from empirical comparison of closely related species with different mating systems.

Distribution of a marker of germline methylation differs between major families of transposon-derived repeats in the human genome

A potential relationship between transposon-derived repeats (TDR) and human germline methylation is of biological importance since many genes are flanked by TDR and methylation could affect the expression of nearby genes. Furthermore, DNA methylation has been suggested as a global defense mechanism against genome instability threatened by TDR. We studied the correlation between the density of HapMap methyl-associated SNPs (mSNPs), a marker of germline methylation, and proportion of TDR. After correcting for confounding variables, we found a negative correlation between proportion of Alu repeats and mSNP density for 125-1000 kb windows. Similar results were found for the most active subgroup of repeats. In contrast, a negative correlation between proportion of L1 repeats and mSNP density was found only in the larger 1000 kb windows. Using methylation data on germ cells (sperm) from the Human Epigenome Project, we found a lower proportion of Alu repeats adjacent (3-15 kb) to hypermethylated amplicons. On the contrary, there was a higher proportion of L1 repeats in the 3-5 kb of sequence flanking hypermethylated amplicons but not in the 10-15 kb flanks. Our data indicate a differential response to the major repeat families and that DNA methylation is unlikely to be a uniform global defense system against all TDR. It appears to play a role for the L1 subgroup, with sequences adjacent to L1 repeats methylated in response to their proximity. In contrast, sequences adjacent to Alu repeats appear to be hypomethylated, arguing against a role of methylation in germline defense against those elements.

Short- and long-term evolutionary dynamics of bacterial insertion sequences: insights from Wolbachia endosymbionts

Transposable elements (TE) are one of the major driving forces of genome evolution, raising the question of the long-term dynamics underlying their evolutionary success. Long-term TE evolution can readily be reconstructed in eukaryotes, thanks to many degraded copies constituting genomic fossil records of past TE proliferations. By contrast, bacterial genomes usually experience high sequence turnover and short TE retention times, thereby obscuring ancient TE evolutionary patterns. We found that Wolbachia bacterial genomes contain 52–171 insertion sequence (IS) TEs. IS account for 11% of Wolbachia wRi, which is one of the highest IS genomic coverage reported in prokaryotes to date. We show that many IS groups are currently expanding in various Wolbachia genomes and that IS horizontal transfers are frequent among strains, which can explain the apparent synchronicity of these IS proliferations. Remarkably, >70% of Wolbachia IS are nonfunctional. They constitute an unusual bacterial IS genomic fossil record providing direct empirical evidence for a long-term IS evolutionary dynamics following successive periods of intense transpositional activity. Our results show that comprehensive IS annotations have the potential to provide new insights into prokaryote TE evolution and, more generally, prokaryote genome evolution. Indeed, the identification of an important IS genomic fossil record in Wolbachia demonstrates that IS elements are not always of recent origin, contrary to the conventional view of TE evolution in prokaryote genomes. Our results also raise the question whether the abundance of IS fossils is specific to Wolbachia or it may be a general, albeit overlooked, feature of prokaryote genomes.

Substitution rate variation at human CpG sites correlates with non-CpG divergence, methylation level and GC content

BACKGROUND

A major goal in the study of molecular evolution is to unravel the mechanisms that induce variation in the germ line mutation rate and in the genome-wide mutation profile. The rate of germ line mutation is considerably higher for cytosines at CpG sites than for any other nucleotide in the human genome, an increase commonly attributed to cytosine methylation at CpG sites. The CpG mutation rate, however, is not uniform across the genome and, as methylation levels have recently been shown to vary throughout the genome, it has been hypothesized that methylation status may govern variation in the rate of CpG mutation.

RESULTS

Here, we use genome-wide methylation data from human sperm cells to investigate the impact of DNA methylation on the CpG substitution rate in introns of human genes. We find that there is a significant correlation between the extent of methylation and the substitution rate at CpG sites. Further, we show that the CpG substitution rate is positively correlated with non-CpG divergence, suggesting susceptibility to factors responsible for the general mutation rate in the genome, and negatively correlated with GC content. We only observe a minor contribution of gene expression level, while recombination rate appears to have no significant effect.

CONCLUSIONS

Our study provides the first direct empirical support for the hypothesis that variation in the level of germ line methylation contributes to substitution rate variation at CpG sites. Moreover, we show that other genomic features also impact on CpG substitution rate variation.

Impact of Alu repeats on the evolution of human p53 binding sites

BACKGROUND

The p53 tumor suppressor protein is involved in a complicated regulatory network, mediating expression of ~1000 human genes. Recent studies have shown that many p53 in vivo binding sites (BSs) reside in transposable repeats. The relationship between these BSs and functional p53 response elements (REs) remains unknown, however. We sought to understand whether the p53 REs also reside in transposable elements and particularly in the most-abundant Alu repeats.

RESULTS

We have analyzed ~160 functional p53 REs identified so far and found that 24 of them occur in repeats. More than half of these repeat-associated REs reside in Alu elements. In addition, using a position weight matrix approach, we found ~400,000 potential p53 BSs in Alu elements genome-wide. Importantly, these putative BSs are located in the same regions of Alu repeats as the functional p53 REs - namely, in the vicinity of Boxes A/A' and B of the internal RNA polymerase III promoter. Earlier nucleosome-mapping experiments showed that the Boxes A/A' and B have a different chromatin environment, which is critical for the binding of p53 to DNA. Here, we compare the Alu-residing p53 sites with the corresponding Alu consensus sequences and conclude that the p53 sites likely evolved through two different mechanisms - the sites overlapping with the Boxes A/A' were generated by CG → TG mutations; the other sites apparently pre-existed in the progenitors of several Alu subfamilies, such as AluJo and AluSq. The binding affinity of p53 to the Alu-residing sites generally correlates with the age of Alu subfamilies, so that the strongest sites are embedded in the 'relatively young' Alu repeats.

CONCLUSIONS

The primate-specific Alu repeats play an important role in shaping the p53 regulatory network in the context of chromatin. One of the selective factors responsible for the frequent occurrence of Alu repeats in introns may be related to the p53-mediated regulation of Alu transcription, which, in turn, influences expression of the host genes.

Polymorphic major histocompatibility complex class II Alu insertions at five loci and their association with HLA-DRB1 and -DQB1 in Japanese and Caucasians

We investigated polymorphic Alu insertion (POALIN) frequencies at five loci in the major histocompatibility complex (MHC) class II genomic region to determine their allele and haplotype frequencies and associations with the human leukocyte antigen (HLA)-DRB1 and -DQB1 genes for 100 Japanese, 174 Australian Caucasians and 67 HLA reference cell lines obtained from different ethnic groups. The POALINs varied in frequency between 11% and 57% with significant differences between the Japanese and Caucasians at three loci. One POALIN locus deviated significantly from Hardy-Weinberg equilibrium (HWE) and four POALIN loci were in significant linkage disequilibrium and had a high percentage association with a variety of HLA-DRB1 or -DQB1 two-digit alleles. Inferred haplotype analysis among two-locus, five-locus and seven-locus haplotype structures showed maximum differences between the Japanese and Caucasians with the seven-locus haplotypes. The most common multilocus haplotype in Caucasians was DRB1*1501/DQB1*0602/AluDQ1/AluDRB1/AluORF10/AluDPB2 (6.7%), whereas the second most common allele HLA-DRB1*15 (17.5%) in Japanese was associated with three or four Alu insertions. The HLA class II POALINs also differentiated within and between HLA-DRB1 super-haplotypes DR1, DR8, DR51, DR52 and DR53. This is the first comparative population study of multilocus POALINs in the HLA class II region, which shows that POALINs whether investigated alone or together with the HLA class II alleles are informative genetic markers for the identification of allele and haplotype lineages and variations within the same and/or different populations.

Computational methods for the analysis of primate mobile elements

Transposable elements (TE), defined as discrete pieces of DNA that can move from one site to another site in genomes, represent significant components of eukaryotic genomes, including primates. Comparative genome-wide analyses have revealed the considerable structural and functional impact of TE families on primate genomes. Insights into these questions have come in part from the development of computational methods that allow detailed and reliable identification, annotation, and evolutionary analyses of the many TE families that populate primate genomes. Here, we present an overview of these computational methods and describe efficient data mining strategies for providing a comprehensive picture of TE biology in newly available genome sequences.

The impact of retrotransposons on human genome evolution

Their ability to move within genomes gives transposable elements an intrinsic propensity to affect genome evolution. Non-long terminal repeat (LTR) retrotransposons--including LINE-1, Alu and SVA elements--have proliferated over the past 80 million years of primate evolution and now account for approximately one-third of the human genome. In this Review, we focus on this major class of elements and discuss the many ways that they affect the human genome: from generating insertion mutations and genomic instability to altering gene expression and contributing to genetic innovation. Increasingly detailed analyses of human and other primate genomes are revealing the scale and complexity of the past and current contributions of non-LTR retrotransposons to genomic change in the human lineage.

Loss of epigenetic silencing in tumors preferentially affects primate-specific retroelements

Close to 50% of the human genome harbors repetitive sequences originally derived from mobile DNA elements, and in normal cells, this sequence compartment is tightly regulated by epigenetic silencing mechanisms involving chromatin-mediated repression. In cancer cells, repetitive DNA elements suffer abnormal demethylation, with potential loss of silencing. We used a genome-wide microarray approach to measure DNA methylation changes in cancers of the head and neck and to compare these changes to alterations found in adjacent non-tumor tissues. We observed specific alterations at thousands of small clusters of CpG dinucleotides associated with DNA repeats. Among the 257,599 repetitive elements probed, 5% to 8% showed disease-related DNA methylation alterations. In dysplasia, a large number of local events of loss of methylation appear in apparently stochastic fashion. Loss of DNA methylation is most pronounced for certain members of the SVA, HERV, LINE-1P, AluY, and MaLR families. The methylation levels of retrotransposons are discretely stratified, with younger elements being highly methylated in healthy tissues, while in tumors, these young elements suffer the most dramatic loss of methylation. Wilcoxon test statistics reveals that a subset of primate LINE-1 elements is demethylated preferentially in tumors, as compared to non-tumoral adjacent tissue. Sequence analysis of these strongly demethylated elements reveals genomic loci harboring full length, as opposed to truncated elements, while possible enrichment for functional LINE-1 ORFs is weaker. Our analysis suggests that, in non-tumor adjacent tissues, there is generalized and highly variable disruption of epigenetic control across the repetitive DNA compartment, while in tumor cells, a specific subset of LINE-1 retrotransposons that arose during primate evolution suffers the most dramatic DNA methylation alterations.

Alu-mediated acquisition of unstable ATTCT pentanucleotide repeats in the human ATXN10 gene

Spinocerebellar ataxia type 10 is caused by ATTCT repeat expansion in the ATXN10 gene in humans. We studied the evolutionary history of the human genome to determine the time and mechanism of the acquisition of unstable ATTCT repeats in the genome. We found that long interspersed element-1 (LINE-1) was inserted into ATXN10 intron 9; Alu was then inserted in the middle of LINE-1; and endogenous retrovilcus K was lastly retrotransposed in the middle of Alu. The ATTCT repeat was located on the boundary between the 3'-end of the Alu element and the direct repeat arising from LINE-1. We determined nucleotide sequences of the orthologous region of 50 individuals representing 33 primate species and compared them with the human sequence. The analysis revealed that the ATTCT repeat is present only in human and apes. Old World monkeys also possess pentanucleotide repeats, but their motifs are TGTCT and GGTCT. New World monkeys and prosimians are not informative because they lack the corresponding region in ATXN10 intron 9. Our studies dictate two parsimonious scenarios of evolution. First, a TTTCT motif arose from a TTTTT motif at the junction of Alu and LINE-1, which was followed by introduction of A to make an ATTCT motif in hominoids. Second, an ATTCT motif was directly generated from an ancestral ATTTT motif in the common ancestor of catarrhines. We also demonstrate that orangutan uniquely introduced G to make a GTTCT motif and later C to make a GTTCC motif, where newly introduced nucleotides are underlined. Our studies reveal that nucleotide substitutions in a poly(A) tail of the Alu element and the following amplification of pentanucleotides occurred in the lineages of Old World monkeys and hominoids and that unstable ATTCT pentanucleotide repeats originated in the common ancestor of hominoids. These findings also highlight a new aspect of the role of retrotransposons in human disease and evolution, which might be useful in investigating the mystery of human uniqueness.

Evolutionary breakpoints in the gibbon suggest association between cytosine methylation and karyotype evolution

Gibbon species have accumulated an unusually high number of chromosomal changes since diverging from the common hominoid ancestor 15-18 million years ago. The cause of this increased rate of chromosomal rearrangements is not known, nor is it known if genome architecture has a role. To address this question, we analyzed sequences spanning 57 breaks of synteny between northern white-cheeked gibbons (Nomascus l. leucogenys) and humans. We find that the breakpoint regions are enriched in segmental duplications and repeats, with Alu elements being the most abundant. Alus located near the gibbon breakpoints (<150 bp) have a higher CpG content than other Alus. Bisulphite allelic sequencing reveals that these gibbon Alus have a lower average density of methylated cytosine that their human orthologues. The finding of higher CpG content and lower average CpG methylation suggests that the gibbon Alu elements are epigenetically distinct from their human orthologues. The association between undermethylation and chromosomal rearrangement in gibbons suggests a correlation between epigenetic state and structural genome variation in evolution.

High-throughput sequence-based epigenomic analysis of Alu repeats in human cerebellum

DNA methylation, the only known covalent modification of mammalian DNA, occurs primarily in CpG dinucleotides. 51% of CpGs in the human genome reside within repeats, and 25% within Alu elements. Despite that, no method has been reported for large-scale ascertainment of CpG methylation in repeats. Here we describe a sequencing-based strategy for parallel determination of the CpG-methylation status of thousands of Alu repeats, and a computation algorithm to design primers that enable their specific amplification from bisulfite converted genomic DNA. Using a single primer pair, we generated amplicons of high sequence complexity, and derived CpG-methylation data from 31 178 Alu elements and their 5′ flanking sequences, altogether representing over 4 Mb of a human cerebellum epigenome. The analysis of the Alu methylome revealed that the methylation level of Alu elements is high in the intronic and intergenic regions, but low in the regions close to transcription start sites. Several hypomethylated Alu elements were identified and their hypomethylated status verified by pyrosequencing. Interestingly, some Alu elements exhibited a strikingly tissue-specific pattern of methylation. We anticipate the amplicons herein described to prove invaluable as epigenome representations, to monitor epigenomic alterations during normal development, in aging and in diseases such as cancer.

Source gene composition and gene conversion of the AluYh and AluYi lineages of retrotransposons

Background

Alu elements are a family of SINE retrotransposons in primates. They are classified into subfamilies according to specific diagnostic mutations from the general Alu consensus. It is now believed that there may be several retrotranspositionally-competent source genes within an Alu subfamily. In this study, subfamilies falling on the AluYi and AluYh lineages, and the AluYg6 subfamily, are assessed for the presence of secondary source genes, and the influence of gene conversion on the AluYh and AluYi lineages is also described.

Results

The AluYh7 and AluYi6 subfamilies appear to contain multiple source genes. The novel subfamilies AluYh3a1 and AluYh3a3 are described, for which there is no convincing evidence to suggest the presence of secondary sources. The mutational substructure of AluYh3a3 can be explained completely by inference of single master gene. A complete backwards gene conversion event appears to have inactivated the AluYh3a3 master gene in humans. Polymorphism data suggest a larger number of secondary source elements may be active in the AluYg6 family than previously thought.

Conclusion

It is clear that there is considerable variation in the number of source genes present in each of the young Alu subfamilies. This can range from a single master source gene, as for AluYh3a3, to as many as 14 source elements in AluYi6.

Estimating the age of retrotransposon subfamilies using maximum likelihood

We present a maximum likelihood model to estimate the age of retrotransposon subfamilies. This method is designed around a master gene model which assumes a constant retrotransposition rate. The statistical properties of this model and an ad hoc estimation procedure are compared using two simulated data sets. We also test whether each estimation procedure is robust to violation of the master gene model. According to our results, both estimation procedures are accurate under the master gene model. While both methods tend to overestimate ages under the intermediate model, the maximum likelihood estimate is significantly less inflated than the ad hoc estimate. We estimate the ages of two subfamilies of human-specific LINE-I insertions using both estimation procedures. By calculating confidence intervals around the maximum likelihood estimate, our model can both provide an estimate of retrotransposon subfamily age and describe the range of subfamily ages consistent with the data.

Analysis of newly identified low copy AluYj subfamily

Human specific AluY elements were investigated by comparative analysis between human chromosome 21 and chimpanzee chromosome 22. Human specific AluY element was identified on human chromosome 21q22 (accession no. AL163282), and then that was a new member of AluYj subfamily. From the bioinformatic analysis, AluYj subfamily was investigated in human whole genome using AluYj4 consensus sequence (accession no. AL163282). Thirteen members of the AluYj4 elements (4 diagnostic mutations) and eight members of the AluYj3 elements (3 diagnostic mutations) were identified with distinct diagnostic mutation from AluY consensus sequence. The results of the molecular clock calculation of non-CpG region substitution indicated that, AluYj4 elements (2.1 million years old) may be proliferated more recent time than AluYj3 elements (14.1 million years old). For the verification of recent insertion time, four of AluYj4 elements (ch2-AC017101, ch10-AC044786, ch12-AC007656 and ch21-AL163282) from human chromosomes 2, 10, 12, 21 were analyzed by PCR amplification using various human and primate DNA samples. Though, no polymorphism was detected in human population, we identified the new AluYj4 subfamily as the human specific elements.

Methylation and deamination of CpGs generate p53-binding sites on a genomic scale

The formation of transcription-factor-binding sites is an important evolutionary process. Here, we show that methylation and deamination of CpG dinucleotides generate in vivo p53-binding sites in numerous Alu elements and in non-repetitive DNA in a species-specific manner. In light of this, we propose that the deamination of methylated CpGs constitutes a universal mechanism for de novo generation of various transcription-factor-binding sites in Alus.

Using Alu elements as polyadenylation sites: A case of retroposon exaptation

Of the 1.1 million Alu retroposons in the human genome, about 10,000 are inserted in the 3' untranslated regions (UTR) of protein-coding genes and 1% of these (107 events) are active as polyadenylation sites (PASs). Strikingly, although Alu's in 3' UTR are indifferently inserted in the forward or reverse direction, 99% of polyadenylation-active Alu sequences are forward oriented. Consensus Alu+ sequences contain sites that can give rise to polyadenylation signals and enhancers through a few point mutations. We found that the strand bias of polyadenylation-active Alu's reflects a radical difference in the fitness of sense and antisense Alu's toward cleavage/polyadenylation activity. In contrast to previous beliefs, Alu inserts do not necessarily represent weak or cryptic PASs; instead, they often constitute the major or the unique PAS in a gene, adding to the growing list of Alu exaptations. Finally, some Alu-borne PASs are intronic and produce truncated transcripts that may impact gene function and/or contribute to gene remodeling.

Regulation of mammalian gene expression by retroelements and non-coding tandem repeats

Genomes of higher eukaryotes contain abundant non-coding repeated sequences whose overall biological impact is unclear. They comprise two categories. The first consists of retrotransposon-derived elements. These are three major families of retroelements (LINEs, SINEs and LTRs). SINEs are clustered in gene-rich regions and are found in promoters of genes while LINEs are concentrated in gene-poor regions and are depleted from promoters. The second class consists of non-coding tandem repeats (satellite DNAs and TTAGGG arrays), which are associated with mammalian centromeres, heterochromatin and telomeres. Terminal TTAGGG arrays are involved in telomere capping and satellite DNAs are located in heterochromatin, which is implicated in transcription silencing by gene repositioning (relocalization). It is unknown whether interstitial TTAGGG sequences, which are present in many vertebrates, have a function. Here, evidence will be presented that retroelements and TTAGGG arrays are involved in regulation of gene expression. Retroelements can provide binding sites for transcription factors and protect promoter CpG islands from repressive chromatin modifications, and may be also involved in nuclear compartmentalization of transcriptionally active and inactive domains. Interstitial telomere-like sequences can form dynamically maintained three-dimensional nuclear networks of transcriptionally inactive domains, which may be involved in transcription silencing like classic heterochromatin.

What do unicellular organisms teach us about DNA methylation?

DNA methylation is an epigenetic hallmark that has been studied intensively in mammals and plants. However, knowledge of this phenomenon in unicellular organisms is scanty. Examining epigenetic regulation, and more specifically DNA methylation, in these organisms represents a unique opportunity to better understand their biology. The determination of their methylation status is often complicated by the presence of several differentiation stages in their life cycle. This article focuses on some recent advances that have revealed the unexpected nature of the epigenetic determinants present in protozoa. The role of the enigmatic DNA methyltransferase Dnmt2 in unicellular organisms is discussed.

Multiple waves of recent DNA transposon activity in the bat, Myotis lucifugus

DNA transposons, or class 2 transposable elements, have successfully propagated in a wide variety of genomes. However, it is widely believed that DNA transposon activity has ceased in mammalian genomes for at least the last 40 million years. We recently reported evidence for the relatively recent activity of hAT and Helitron elements, two distinct groups of DNA transposons, in the lineage of the vespertilionid bat Myotis lucifugus. Here, we describe seven additional families that have also been recently active in the bat lineage. Early vespertilionid genome evolution was dominated by the activity of Helitrons, mariner-like and Tc2-like elements. This was followed by the colonization of Tc1-like elements, and by a more recent explosion of hAT-like elements. Finally, and most recently, piggyBac-like elements have amplified within the Myotis genome and our results indicate that one of these families is probably still expanding in natural populations. Together, these data suggest that there has been tremendous recent activity of various DNA transposons in the bat lineage that far exceeds those previously reported for any mammalian lineage. The diverse and recent populations of DNA transposons in genus Myotis will provide an unprecedented opportunity to study the impact of this class of elements on mammalian genome evolution and to better understand what makes some species more susceptible to invasion by genomic parasites than others.

Genome-wide tracking of unmethylated DNA Alu repeats in normal and cancer cells

Methylation of the cytosine is the most frequent epigenetic modification of DNA in mammalian cells. In humans, most of the methylated cytosines are found in CpG-rich sequences within tandem and interspersed repeats that make up to 45% of the human genome, being Alu repeats the most common family. Demethylation of Alu elements occurs in aging and cancer processes and has been associated with gene reactivation and genomic instability. By targeting the unmethylated SmaI site within the Alu sequence as a surrogate marker, we have quantified and identified unmethylated Alu elements on the genomic scale. Normal colon epithelial cells contain in average 25 486 ± 10 157 unmethylated Alu's per haploid genome, while in tumor cells this figure is 41 995 ± 17 187 (P = 0.004). There is an inverse relationship in Alu families with respect to their age and methylation status: the youngest elements exhibit the highest prevalence of the SmaI site (AluY: 42%; AluS: 18%, AluJ: 5%) but the lower rates of unmethylation (AluY: 1.65%; AluS: 3.1%, AluJ: 12%). Data are consistent with a stronger silencing pressure on the youngest repetitive elements, which are closer to genes. Further insights into the functional implications of atypical unmethylation states in Alu elements will surely contribute to decipher genomic organization and gene regulation in complex organisms.