Standardized nomenclature for Alu repeats

Association of angiotensin I converting enzyme gene I/D polymorphism and coronary artery disease in the Pakistani population

ObjectiveThe aim of this study was to investigate the association between angiotensin I converting enzyme gene polymorphism and coronary artery disease in the Pakistani population, given the early onset and aggressive nature of coronary artery disease in this region.MethodsA case-control study was conducted involving 540 Pakistani patients with established coronary artery disease and 224 healthy controls. DNA samples were amplified using polymerase chain reaction with primers targeting the insertion (I)/deletion (D) sites of intron 16 of the angiotensin I converting enzyme gene. The polymerase chain reaction products were analyzed for the presence of 490-bp (II), 190-bp (DD), or both (ID) fragments.ResultsThe frequency of the homozygous insertion (II) genotype was 9% in the control group and 25% in coronary artery disease patients, while the homozygous deletion (DD) genotype was 26% in controls and 24% in patients. A significant association was found between the angiotensin I converting enzyme II genotype and coronary artery disease (odds ratio = 7.963, p < 0.0001).ConclusionThe angiotensin I converting enzyme II genotype is significantly associated with an increased risk of coronary artery disease in the Pakistani population, suggesting a potential genetic predisposition to early and aggressive coronary artery disease in this group.

Inverted Alu repeats in loop-out exon skipping across hominoid evolution

Background

Changes in RNA splicing over the course of evolution have profoundly diversified the functional landscape of the human genome. While DNA sequences proximal to intron-exon junctions are known to be critical for RNA splicing, the impact of distal intronic sequences remains underexplored. Emerging evidence suggests that inverted pairs of intronic Alu elements can promote exon skipping by forming RNA stem-loop structures. However, their prevalence and influence throughout evolution remain unknown.

Results

Here, we present a systematic analysis of inverted Alu pairs across the human genome to assess their impact on exon skipping through predicted RNA stem-loop formation and their relevance to hominoid evolution. We found that inverted Alu pairs, particularly pairs of AluY-AluSx1 and AluSz-AluSx, are enriched in the flanking regions of skippable exons genome-wide and are predicted to form stable stem-loop structures. Exons defined by weak 3' acceptor and strong 5' donor splice sites appear especially prone to this skipping mechanism. Through comparative genome analysis across nine primate species, we identified 67,126 hominoid-specific Alu insertions, primarily from AluY and AluS subfamilies, which form inverted pairs enriched across skippable exons in genes of ubiquitination-related pathways. Experimental validation of exon skipping among several hominoid-specific inverted Alu pairs further reinforced their potential evolutionary significance.

Conclusion

This work extends our current knowledge of the roles of RNA secondary structure formed by inverted Alu pairs and details a newly emerging mechanism through which transposable elements have contributed to genomic innovation across hominoid evolution at the transcriptomic level.

Comparative Genomics Reveals LINE-1 Recombination with Diverse RNAs

Long interspersed element-1 (LINE-1, L1) retrotransposons are the most abundant protein-coding transposable elements (TE) in mammalian genomes, and have shaped genome content over 170 million years of evolution. LINE-1 is self-propagating and mobilizes other sequences, including Alu elements. Occasionally, LINE-1 forms chimeric insertions with non-coding RNAs and mRNAs. U6 spliceosomal small nuclear RNA/LINE-1 chimeras are best known, though there are no comprehensive catalogs of LINE-1 chimeras. To address this, we developed TiMEstamp, a computational pipeline that leverages multiple sequence alignments (MSA) to estimate the age of LINE-1 insertions and identify candidate chimeric insertions where an adjacent sequence arrives contemporaneously. Candidates were refined by detecting hallmark features of L1 retrotransposition, such as target site duplication (TSD). Applying this pipeline to the human genome, we recovered all known species of LINE-1 chimeras and discovered new chimeric insertions involving small RNAs, Alu elements, and mRNA fragments. Some insertions are compatible with known mechanisms, such as RNA ligation. Other structures nominate novel mechanisms, such as trans-splicing. We also see evidence that LINE-1 loci with defunct promoters can acquire regulatory elements from nearby genes to restore retrotransposition activity. These discoveries highlight the recombinatory potential of LINE-1 RNA with implications for genome evolution and TE domestication.

Cancer cells restrict immunogenicity of retrotransposon expression via distinct mechanisms

To thrive, cancer cells must navigate acute inflammatory signaling accompanying oncogenic transformation, such as via overexpression of repeat elements. We examined the relationship between immunostimulatory repeat expression, tumor evolution, and the tumor-immune microenvironment. Integration of multimodal data from a cohort of pancreatic ductal adenocarcinoma (PDAC) patients revealed expression of specific Alu repeats predicted to form double-stranded RNAs (dsRNAs) and trigger retinoic-acid-inducible gene I (RIG-I)-like-receptor (RLR)-associated type-I interferon (IFN) signaling. Such Alu-derived dsRNAs also anti-correlated with pro-tumorigenic macrophage infiltration in late stage tumors. We defined two complementary pathways whereby PDAC may adapt to such anti-tumorigenic signaling. In mutant TP53 tumors, ORF1p from long interspersed nuclear element (LINE)-1 preferentially binds Alus and decreases their expression, whereas adenosine deaminases acting on RNA 1 (ADAR1) editing primarily reduces dsRNA formation in wild-type TP53 tumors. Depletion of either LINE-1 ORF1p or ADAR1 reduced tumor growth in vitro. The fact that tumors utilize multiple pathways to mitigate immunostimulatory repeats implies the stress from their expression is a fundamental phenomenon to which PDAC, and likely other tumors, adapt.

DSB profiles in human spermatozoa highlight the role of TMEJ in the male germline

The male mammalian germline is characterized by substantial chromatin remodeling associated with the transition from histones to protamines during spermatogenesis, followed by the reversal to nucleohistones in the male pronucleus preceding the zygotic genome activation. Both transitions are associated with the extensive formation of DNA double-strand breaks (DSBs), requiring an estimated 5 to 10 million transient DSBs per spermatozoa. Additionally, the high transcription rate in early stages of spermatogenesis leads to transcription-coupled damage preceding meiotic homologous recombination, potentially further contributing to the DSB landscape in mature spermatozoa. Once meiosis is completed, spermatozoa remain haploid and therefore cannot rely on error-free homologous recombination, but instead depend on error-prone classical non-homologous end joining (cNHEJ). This DNA damage/repair-scenario is proposed to be one of the main causes of the observed paternal mutation propensity in human evolution. Recent studies have shown that DSBs in the male pronucleus are repaired by maternally provided Polθ in Caenorhabditis elegans through Polθ-mediated end joining (TMEJ). Additionally, population genetic datasets have revealed a preponderance of TMEJ signatures associated with human variation. Since these signatures are the result of the combined effect of TMEJ and DSB formation in spermatozoa and male pronuclei, we used a BLISS-based protocol to analyze recurrent DSBs in mature human sperm heads as a proxy of the male pronucleus before zygotic chromatin remodeling. The DSBs were found to be enriched in (YR)n short tandem repeats and in evolutionarily young SINEs, reminiscent to patterns observed in murine spermatids, indicating evolutionary hotspots of recurrent DSB formation in mammalian spermatozoa. Additionally, we detected a similar DSB pattern in diploid human IMR90 cells when cNHEJ was selectively inhibited, indicating the significant impact of absent cNHEJ on the sperm DSB landscape. Strikingly, regions associated with most retained histones, and therefore less condensed chromatin, were not strongly enriched with recurrent DSBs. In contrast, the fraction of retained H3K27me3 in the mature spermatozoa displayed a strong association with recurrent DSBs. DSBs in H3K27me3 are associated with a preference for TMEJ over cNHEJ during repair. We hypothesize that the retained H3K27me3 may trigger transgenerational DNA repair by priming maternal Polθ to these regions.

Helenus and Ajax, Two Groups of Non-Autonomous LTR Retrotransposons, Represent a New Type of Small RNA Gene-Derived Mobile Elements

Terminal repeat retrotransposons in miniature (TRIMs) are short non-autonomous long terminal repeat (LTR) retrotransposons found from various eukaryotes. Cassandra is a unique TRIM lineage which contains a 5S rRNA-derived sequence in its LTRs. Here, two new groups of TRIMs, designated Helenus and Ajax, are reported based on bioinformatics analysis and the usage of Repbase. Helenus is found from fungi, animals, and plants, and its LTRs contain a tRNA-like sequence. It includes two LTRs and between them, a primer-binding site (PBS) and polypurine tract (PPT) exist. Fungal and plant Helenus generate 5 bp target site duplications (TSDs) upon integration, while animal Helenus generates 4 bp TSDs. Ajax includes a 5S rRNA-derived sequence in its LTR and is found from two nemertean genomes. Ajax generates 5 bp TSDs upon integration. These results suggest that despite their unique promoters, Helenus and Ajax are TRIMs whose transposition is dependent on autonomous LTR retrotransposon. These TRIMs can originate through an insertion of SINE in an LTR of TRIM. The discovery of Helenus and Ajax suggests the presence of TRIMs with a promoter for RNA polymerase III derived from a small RNA gene, which is here collectively termed TRIMp3.

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.

Cancer cells co-evolve with retrotransposons to mitigate viral mimicry

Overexpression of repetitive elements is an emerging hallmark of human cancers ¹ . Diverse repeats can mimic viruses by replicating within the cancer genome through retrotransposition, or presenting pathogen-associated molecular patterns (PAMPs) to the pattern recognition receptors (PRRs) of the innate immune system ^2-5 . Yet, how specific repeats affect tumor evolution and shape the tumor immune microenvironment (TME) in a pro- or anti-tumorigenic manner remains poorly defined. Here, we integrate whole genome and total transcriptome data from a unique autopsy cohort of multiregional samples collected in pancreatic ductal adenocarcinoma (PDAC) patients, into a comprehensive evolutionary analysis. We find that more recently evolved S hort I nterspersed N uclear E lements (SINE), a family of retrotransposable repeats, are more likely to form immunostimulatory double-strand RNAs (dsRNAs). Consequently, younger SINEs are strongly co-regulated with RIG-I like receptor associated type-I interferon genes but anti-correlated with pro-tumorigenic macrophage infiltration. We discover that immunostimulatory SINE expression in tumors is regulated by either L ong I nterspersed N uclear E lements 1 (LINE1/L1) mobility or ADAR1 activity in a TP53 mutation dependent manner. Moreover, L1 retrotransposition activity tracks with tumor evolution and is associated with TP53 mutation status. Altogether, our results suggest pancreatic tumors actively evolve to modulate immunogenic SINE stress and induce pro-tumorigenic inflammation. Our integrative, evolutionary analysis therefore illustrates, for the first time, how dark matter genomic repeats enable tumors to co-evolve with the TME by actively regulating viral mimicry to their selective advantage.

Framework of the Alu Subfamily Evolution in the Platyrrhine Three-Family Clade of Cebidae, Callithrichidae, and Aotidae

The history of Alu retroposons has been choreographed by the systematic accumulation of inherited diagnostic nucleotide substitutions to form discrete subfamilies, each having a distinct nucleotide consensus sequence. The oldest subfamily, AluJ, gave rise to AluS after the split between Strepsirrhini and what would become Catarrhini and Platyrrhini. The AluS lineage gave rise to AluY in catarrhines and to AluTa in platyrrhines. Platyrrhine Alu subfamilies Ta7, Ta10, and Ta15 were assigned names based on a standardized nomenclature. However, with the subsequent intensification of whole genome sequencing (WGS), large scale analyses to characterize Alu subfamilies using the program COSEG identified entire lineages of subfamilies simultaneously. The first platyrrhine genome with WGS, the common marmoset (Callithrix jacchus; [caljac3]), resulted in Alu subfamily names sf0 to sf94 in an arbitrary order. Although easily resolved by alignment of the consensus sequences, this naming convention can become increasingly confusing as more genomes are independently analyzed. In this study, we reported Alu subfamily characterization for the platyrrhine three-family clade of Cebidae, Callithrichidae, and Aotidae. We investigated one species/genome from each recognized family of Callithrichidae and Aotidae and of both subfamilies (Cebinae and Saimiriinae) of the family Cebidae. Furthermore, we constructed a comprehensive network of Alu subfamily evolution within the three-family clade of platyrrhines to provide a working framework for future research. Alu expansion in the three-family clade has been dominated by AluTa15 and its derivatives.

RNA editing of microtubule-associated protein tau circular RNAs promotes their translation and tau tangle formation

Aggregation of the microtubule-associated protein tau characterizes tauopathies, including Alzheimer's disease and frontotemporal lobar degeneration (FTLD-Tau). Gene expression regulation of tau is complex and incompletely understood. Here we report that the human tau gene (MAPT) generates two circular RNAs (circRNAs) through backsplicing of exon 12 to either exon 7 (12→7 circRNA) or exon 10 (12→10 circRNA). Both circRNAs lack stop codons. The 12→7 circRNA contains one start codon and is translated in a rolling circle, generating a protein consisting of multimers of the microtubule-binding repeats R1-R4. For the 12→10 circRNA, a start codon can be introduced by two FTLD-Tau mutations, generating a protein consisting of multimers of the microtubule-binding repeats R2-R4, suggesting that mutations causing FTLD may act in part through tau circRNAs. Adenosine to inosine RNA editing dramatically increases translation of circRNAs and, in the 12→10 circRNA, RNA editing generates a translational start codon by changing AUA to AUI. Circular tau proteins self-aggregate and promote aggregation of linear tau proteins. Our data indicate that adenosine to inosine RNA editing initiates translation of human circular tau RNAs, which may contribute to tauopathies.

A 69 kb Deletion in chr19q13.42 including PRPF31 Gene in a Chinese Family Affected with Autosomal Dominant Retinitis Pigmentosa

We aimed to identify the genetic cause of autosomal dominant retinitis pigmentosa (adRP) and characterize the underlying molecular mechanisms of incomplete penetrance in a Chinese family affected with adRP. All enrolled family members underwent ophthalmic examinations. Whole-genome sequencing (WGS), multiplex ligation-dependent probe amplification (MLPA), linkage analysis and haplotype construction were performed in all participants. RNA-seq was performed to analyze the regulating mechanism of incomplete penetrance among affected patients, mutation carriers and healthy controls. In the studied family, 14 individuals carried a novel heterozygous large deletion of 69 kilobase (kb) in 19q13.42 encompassing exon 1 of the PRPF31 gene and five upstream genes: TFPT, OSCAR, NDUFA3, TARM1, and VSTM1. Three family members were sequenced and diagnosed as non-penetrant carriers (NPCs). RNA-seq showed significant differential expression of genes in deletion between mutation carriers and healthy control. The RP11 pedigree in this study was the largest pedigree compared to other reported RP11 pedigrees with large deletions. Early onset in all affected members in this pedigree was considered to be a special phenotype and was firstly reported in a RP11 family for the first time. Differential expression of PRPF31 between affected and unaffected subjects indicates a haploinsufficiency to cause the disease in the family. The other genes with significant differential expression might play a cooperative effect on the penetrance of RP11.

Owl Monkey Alu Insertion Polymorphisms and Aotus Phylogenetics

Owl monkeys (genus Aotus), or "night monkeys" are platyrrhine primates in the Aotidae family. Early taxonomy only recognized one species, Aotus trivirgatus, until 1983, when Hershkovitz proposed nine unique species designations, classified into red-necked and gray-necked species groups based predominately on pelage coloration. Recent studies questioned this conventional separation of the genus and proposed designations based on the geographical location of wild populations. Alu retrotransposons are a class of mobile element insertion (MEI) widely used to study primate phylogenetics. A scaffold-level genome assembly for one Aotus species, Aotus nancymaae [Anan_2.0], facilitated large-scale ascertainment of nearly 2000 young lineage-specific Alu insertions. This study provides candidate oligonucleotides for locus-specific PCR assays for over 1350 of these elements. For 314 Alu elements across four taxa with multiple specimens, PCR analyses identified 159 insertion polymorphisms, including 21 grouping A. nancymaae and Aotus azarae (red-necked species) as sister taxa, with Aotus vociferans and A. trivirgatus (gray-necked) being more basal. DNA sequencing identified five novel Alu elements from three different taxa. The Alu datasets reported in this study will assist in species identification and provide a valuable resource for Aotus phylogenetics, population genetics and conservation strategies when applied to wild populations.

Comprehensive In Silico Analysis of Retrotransposon Insertions within the Survival Motor Neuron Genes Involved in Spinal Muscular Atrophy

Transposable elements (TEs) are interspersed repetitive and mobile DNA sequences within the genome. Better tools for evaluating TE-derived sequences have provided insights into the contribution of TEs to human development and disease. Spinal muscular atrophy (SMA) is an autosomal recessive motor neuron disease that is caused by deletions or mutations in the Survival Motor Neuron 1 (SMN1) gene but retention of its nearly perfect orthologue SMN2. Both genes are highly enriched in TEs. To establish a link between TEs and SMA, we conducted a comprehensive, in silico analysis of TE insertions within the SMN1/2 loci of SMA, carrier and healthy genomes. We found an Alu insertion in the promoter region and one L1 element in the 3'UTR that may play an important role in alternative promoter as well as in alternative transcriptional termination. Additionally, several intronic Alu repeats may influence alternative splicing via RNA circularization and causes the presence of new alternative exons. These Alu repeats present throughout the genes are also prone to recombination events that could lead to SMN1 exons deletions and, ultimately, SMA. TE characterization of the SMA genomic region could provide for a better understanding of the implications of TEs on human disease and genomic evolution.

An antisense Alu transposon insertion/deletion polymorphism of ALDH1A1 may functionally associate with Parkinson's disease

BACKGROUND

Aldehyde dehydrogenase 1 (encoded by ALDH1A1) has been shown to protect against Parkinson's disease (PD) by reducing toxic metabolites of dopamine. We herein revealed an antisense Alu element insertion/deletion polymorphism in intron 4 of ALDH1A1, and hypothesized that it might play a role in PD.  METHODS: A Han Chinese cohort comprising 488 PD patients and 515 controls was recruited to validate the Alu insertion/deletion polymorphism following a previous study of tag-single nucleotide polymorphisms, where rs7043217 was shown to be significantly associated with PD. Functional analyses of the Alu element insertion were performed.

RESULTS

The Alu element of ALDH1A1 was identified to be a variant of Yb8 subfamily and termed as Yb8c4. The antisense Yb8c4 insertion/deletion polymorphism (named asYb8c4^ins and asYb8c4^del, respectively) appeared to be in a complete linkage disequilibrium with rs7043217 and was validated to be significantly associated with PD susceptibility with asYb8c4^ins serving as a risk allele (P = 0.030, OR = 1.224, 95% CI = 1.020-1.470). Multiple functional analyses including ALDH1A1 mRNA expression in blood cells of carriers, and reporters of EGFP and luciferase showed that the asYb8c4^ins had a suppressive activity on gene transcription. Mechanistic explorations suggested that the asYb8c4^ins induced no changes in CpG methylation and mRNA splicing of ALDH1A1 and appeared no binding of transcription factors.

CONCLUSIONS

Our results consolidate an involvement of ALDH1 in PD pathogenesis. The asYb8c4 polymorphism may be a functional output of its linkage disequilibrium-linked single nucleotide polymorphisms.

RNA transcription and degradation of Alu retrotransposons depends on sequence features and evolutionary history

Alu elements are one of the most successful groups of RNA retrotransposons and make up 11% of the human genome with over 1 million individual loci. They are linked to genetic defects, increases in sequence diversity, and influence transcriptional activity. Still, their RNA metabolism is poorly understood yet. It is even unclear whether Alu elements are mostly transcribed by RNA Polymerase II or III. We have conducted a transcription shutoff experiment by α-amanitin and metabolic RNA labeling by 4-thiouridine combined with RNA fragmentation (TT-seq) and RNA-seq to shed further light on the origin and life cycle of Alu transcripts. We find that Alu RNAs are more stable than previously thought and seem to originate in part from RNA Polymerase II activity, as previous reports suggest. Their expression however seems to be independent of the transcriptional activity of adjacent genes. Furthermore, we have developed a novel statistical test for detecting the expression of quantitative trait loci in Alu elements that relies on the de Bruijn graph representation of all Alu sequences. It controls for both statistical significance and biological relevance using a tuned k-mer representation, discovering influential sequence features missed by regular motif search. In addition, we discover several point mutations using a generalized linear model, and motifs of interest, which also match transcription factor-binding motifs.

A high-quality, long-read genome assembly of the endangered ring-tailed lemur (Lemur catta)

BACKGROUND

The ring-tailed lemur (Lemur catta) is a charismatic strepsirrhine primate endemic to Madagascar. These lemurs are of particular interest, given their status as a flagship species and widespread publicity in the popular media. Unfortunately, a recent population decline has resulted in the census population decreasing to <2,500 individuals in the wild, and the species's classification as an endangered species by the IUCN. As is the case for most strepsirrhine primates, only a limited amount of genomic research has been conducted on L. catta, in part owing to the lack of genomic resources.

RESULTS

We generated a new high-quality reference genome assembly for L. catta (mLemCat1) that conforms to the standards of the Vertebrate Genomes Project. This new long-read assembly is composed of Pacific Biosciences continuous long reads (CLR data), Optical Mapping Bionano reads, Arima HiC data, and 10X linked reads. The contiguity and completeness of the assembly are extremely high, with scaffold and contig N50 values of 90.982 and 10.570 Mb, respectively. Additionally, when compared to other high-quality primate assemblies, L. catta has the lowest reported number of Alu elements, which results predominantly from a lack of AluS and AluY elements.

CONCLUSIONS

mLemCat1 is an excellent genomic resource not only for the ring-tailed lemur community, but also for other members of the Lemuridae family, and is the first very long read assembly for a strepsirrhine.

Blood Transcriptome Analysis of Septic Patients Reveals a Long Non-Coding Alu-RNA in the Complement C5a Receptor 1 Gene

Many severe inflammation conditions are complement-dependent with the complement component C5a-C5aR1 axis as an important driver. At the RNA level, the blood transcriptome undergoes programmed expression of coding and long non-coding RNAs to combat invading microorganisms. Understanding the expression of long non-coding RNAs containing Alu elements in inflammation is important for reconstructing cell fate trajectories leading to severe disease. We have assembled a pipeline for computation mining of new Alu-containing long non-coding RNAs by intersecting immune genes with known Alu coordinates in the human genome. By applying the pipeline to patient bulk RNA-seq data with sepsis, we found immune genes containing 48 Alu insertion as robust candidates for further study. Interestingly, 1 of the 48 candidates was located within the complement system receptor gene C5aR1 and holds promise as a target for RNA therapeutics.

Recently Integrated Alu Elements in Capuchin Monkeys: A Resource for Cebus/Sapajus Genomics

Capuchins are platyrrhines (monkeys found in the Americas) within the Cebidae family. For most of their taxonomic history, the two main morphological types of capuchins, gracile (untufted) and robust (tufted), were assigned to a single genus, Cebus. Further, all tufted capuchins were assigned to a single species, Cebus apella, despite broad geographic ranges spanning Central and northern South America. In 2012, tufted capuchins were assigned to their genus, Sapajus, with eight currently recognized species and five Cebus species, although these numbers are still under debate. Alu retrotransposons are a class of mobile element insertion (MEI) widely used to study primate phylogenetics. However, Alu elements have rarely been used to study capuchins. Recent genome-level assemblies for capuchins (Cebus imitator; [Cebus_imitator_1.0] and Sapajus apella [GSC_monkey_1.0]) facilitated large scale ascertainment of young lineage-specific Alu insertions. Reported here are 1607 capuchin specific and 678 Sapajus specific Alu insertions along with candidate oligonucleotides for locus-specific PCR assays for many elements. PCR analyses identified 104 genus level and 51 species level Alu insertion polymorphisms. The Alu datasets reported in this study provide a valuable resource that will assist in the classification of archival samples lacking phenotypic data and for the study of capuchin phylogenetic relationships.

The L1-dependant and Pol III transcribed Alu retrotransposon, from its discovery to innate immunity

The 300 bp dimeric repeats digestible by AluI were discovered in 1979. Since then, Alu were involved in the most fundamental epigenetic mechanisms, namely reprogramming, pluripotency, imprinting and mosaicism. These Alu encode a family of retrotransposons transcribed by the RNA Pol III machinery, notably when the cytosines that constitute their sequences are de-methylated. Then, Alu hijack the functions of ORF2 encoded by another transposons named L1 during reverse transcription and integration into new sites. That mechanism functions as a complex genetic parasite able to copy-paste Alu sequences. Doing that, Alu have modified even the size of the human genome, as well as of other primate genomes, during 65 million years of co-evolution. Actually, one germline retro-transposition still occurs each 20 births. Thus, Alu continue to modify our human genome nowadays and were implicated in de novo mutation causing diseases including deletions, duplications and rearrangements. Most recently, retrotransposons were found to trigger neuronal diversity by inducing mosaicism in the brain. Finally, boosted during viral infections, Alu clearly interact with the innate immune system. The purpose of that review is to give a condensed overview of all these major findings that concern the fascinating physiology of Alu from their discovery up to the current knowledge.

Cytoplasmic synthesis of endogenous Alu complementary DNA via reverse transcription and implications in age-related macular degeneration

Alu retroelements propagate via retrotransposition by hijacking long interspersed nuclear element-1 (L1) reverse transcriptase (RT) and endonuclease activities. Reverse transcription of Alu RNA into complementary DNA (cDNA) is presumed to occur exclusively in the nucleus at the genomic integration site. Whether Alu cDNA is synthesized independently of genomic integration is unknown. Alu RNA promotes retinal pigmented epithelium (RPE) death in geographic atrophy, an untreatable type of age-related macular degeneration. We report that Alu RNA-induced RPE degeneration is mediated via cytoplasmic L1-reverse-transcribed Alu cDNA independently of retrotransposition. Alu RNA did not induce cDNA production or RPE degeneration in L1-inhibited animals or human cells. Alu reverse transcription can be initiated in the cytoplasm via self-priming of Alu RNA. In four health insurance databases, use of nucleoside RT inhibitors was associated with reduced risk of developing atrophic macular degeneration (pooled adjusted hazard ratio, 0.616; 95% confidence interval, 0.493-0.770), thus identifying inhibitors of this Alu replication cycle shunt as potential therapies for a major cause of blindness.

Sequence diversity analyses of an improved rhesus macaque genome enhance its biomedical utility

The rhesus macaque (Macaca mulatta) is the most widely studied nonhuman primate (NHP) in biomedical research. We present an updated reference genome assembly (Mmul_10, contig N50 = 46 Mbp) that increases the sequence contiguity 120-fold and annotate it using 6.5 million full-length transcripts, thus improving our understanding of gene content, isoform diversity, and repeat organization. With the improved assembly of segmental duplications, we discovered new lineage-specific genes and expanded gene families that are potentially informative in studies of evolution and disease susceptibility. Whole-genome sequencing (WGS) data from 853 rhesus macaques identified 85.7 million single-nucleotide variants (SNVs) and 10.5 million indel variants, including potentially damaging variants in genes associated with human autism and developmental delay, providing a framework for developing noninvasive NHP models of human disease.

Retroelement-derived RNA and its role in the brain

Comprising ~40% of the human genome, retroelements are mobile genetic elements which are transcribed into RNA, then reverse-transcribed into DNA and inserted into a new site in the genome. Retroelements are referred to as "genetic parasites", residing among host genes and relying on host machinery for transcription and evolutionary propagation. The healthy brain has the highest expression of retroelement-derived sequences compared to other somatic tissue, which leads to the question: how does retroelement-derived RNA influence human traits and cellular states? While the functional importance of upregulating retroelement expression in the brain is an active area of research, RNA species derived from retroelements influence both self- and host gene expression by contributing to chromatin remodeling, alternative splicing, somatic mosaicism and translational repression. Here, we review the emerging evidence that the functional importance of RNA derived from retroelements is multifaceted. Retroelements can influence organismal states through the seeding of epigenetic states in chromatin, the production of structured RNA and even catalytically active ribozymes, the generation of cytoplasmic ssDNA and RNA/DNA hybrids, the production of viral-like proteins, and the generation of somatic mutations. Comparative sequencing suggests that retroelements can contribute to intraspecies variation through these mechanisms to alter transcript identity and abundance. In humans, an increasing number of neurodevelopmental and neurodegenerative conditions are associated with dysregulated retroelements, including Aicardi-Goutieres syndrome (AGS), Rett syndrome (RTT), Amyotrophic Lateral Sclerosis (ALS), Alzheimer's disease (AD), multiple sclerosis (MS), schizophrenia (SZ), and aging. Taken together, these concepts suggest a larger functional role for RNA derived from retroelements. This review aims to define retroelement-derived RNA, discuss how it impacts the mammalian genome, as well as summarize data supporting phenotypic consequences of this unique RNA subset in the brain.

Pathogenicity reclassification of two BRCA1/BRCA2 exonic duplications after identification of genomic breakpoints and tandem orientation

The genomic consequence and clinical interpretation of large duplications are difficult to infer without determining the location and orientation of the duplicated sequence. We aimed to characterize two intragenic duplications detected in two hereditary breast and ovarian cancer syndrome (HBOC) families, namely BRCA1 exon 4 to 6 and BRCA2 exon 17 to 18, previously detected by multiplex ligation probe amplification and initially classified as variants of unknown significance. Using long range PCR, with duplication-specific primers, we were able to ascertain the genomic breakpoints and observed that the two rearrangements occurred in tandem and in direct orientation. The BRCA1 c.134+440_441+870dup and BRCA2 c.7806-2083_8332-1512dup duplications here identified are predicted to cause frameshifts that create a premature stop codon and were reclassified as pathogenic. Furthermore, both families present phenotypic traits typical of HBOC syndrome. We also observed that the genomic breakpoints of these two duplications occurred within highly homologous Alu elements. Concluding, we characterized two in tandem BRCA1 and BRCA2 duplications that likely occurred by Alu-mediated homologous recombination, allowing identification of the underlying cause of the HBOC syndrome in these families.

Alu master copies serve as the drivers of differential SINE transposition in recent primate genomes

Alu elements, averaging ~300bp in length, are a family of primate-specific short intersperse nuclear elements (SINEs) with more than one million copies and contributing to ~11% of primate genomes. Despite mostly being shared among primates, our recent study revealed highly differential recent Alu transposition among the genomes of primates from Hominidae and Cercopithecidae families. To understand the underlying mechanism, we analyzed six primate genomes and revealed species- and lineage-specific Alu profile exclusively defined by AluY composition. Among all Alus from the 6 genomes, we identified 5401 Alu master copies with 99% being from the AluY subfamily. The numbers of Alu master copies are positively correlated to the number of AluY elements in the genomes with the baboon genome having the largest number of most recent Alu master copies at high activities, while the crab-eating macaque genome having a low number of Alu master copies with low activity. Furthermore, the expression level of Alu master copies is positively correlated with their transposition activity. Our results support the concept that Alu transposition in primate genomes is driven by a small number of master copies, the number and relative activity of which contribute to the differential Alu transposition in recent primate genomes.

Highly specific, quantitative polymerase chain reaction probe for the quantification of human cells in cynomolgus monkeys

Quantification of human cells may be performed using quantitative polymerase chain reaction (qPCR). In preclinical studies, the human Alu sequence is widely used as biomarker for human DNA. However, because the Alu gene is shared by primates, its use is limited to non-primate studies. The biodistribution of human cells in primates is also necessary for translational studies. Therefore, we aimed to design a novel, human-specific primer/probe that enables the quantification of human cells in primates and other animal models. A novel primer/probe set was successfully designed based on highly repetitive LINE1 sequences. qPCR efficiency (94.95-99.21%) and linearity of calibration curves (r² = 0.996-0.999) were confirmed in tissue homogenates of cynomolgus monkey. The lower limit of detection was 10 cells per 15-mg tissue sample, a sensitivity that is equivalent to existing Alu primers/probes. The set was also effective in other animal models such as mice, rabbits, pigs, and common marmosets. To our knowledge, this is the first study describing the successful design of a human-specific qPCR primer/probe for human cell quantification in various animals, including non-human primates, using LINE1 sequence. The excellent selectivity, sensitivity, and versatility of the LINE1 primers/probes make it a promising quantification tool in preclinical biodistribution studies.

A comprehensive analysis of chimpanzee (Pan Troglodytes)-specific AluYb8 element

BACKGROUND

Alu elements are most abundant retrotransposons with > 1.2 million copies in the primate genome. AluYb8 subfamily was diverged from AluY lineage, and has accumulated eight diagnostic mutations and 7-bp duplication during primate evolution. A total of 1851 AluYb copies are present in the human genome, and most of them are human-specific. On the other hand, only a few AluYb8 copies were identified in the chimpanzee genome by previous studies on AluYb8. The significantly different number of species-specific AluYb8 elements between human and chimpanzee might result from the incompletion of chimpanzee reference genome sequences at the time of the previous study.

OBJECTIVE

AluYb8 elements could generate genomic structural variations in the chimpanzee genome. This study aimed to identify and characterize chimpanzee-specific AluYb elements using the most updated chimpanzee reference genome sequences (Jan. 2018, panTro6).

METHODS

To identify chimpanzee-specific AluYb8, we carried out genomic comparison with non-chimpanzee primate genome using the UCSC table browser. In addition, chimpanzee-specific AluYb8 candidates were manually inspected and experimentally verified using PCR and Sanger sequencing.

RESULTS

Among a total of 231 chimpanzee-specific AluYb8 candidates, 11 of the candidates are chimpanzee-specific AluYb8, and 29 elements are shared between the chimpanzee and non-chimpanzee primate genomes. Through the sequence analysis of AluYb8 and other Alu subfamilies, we were able to observe various diagnostic mutations and variable length duplications in 7-bp duplication region of AluYb8 element. In addition, we further validated two of the chimpanzee-specific AluYb8 elements (CS8 and CS20) that were not previously discovered by display PCR and Sanger sequencing. Interestingly, we identified a AluYb8 insertion-mediated deletion (CS8 locus) in the chimpanzee genome.

CONCLUSION

Our study found that AluYb8 elements are much more abundant in the human genome than chimpanzee genome, and that it could be due to the absence of hyperactive "master" AluYb8 elements in the chimpanzee genome.

Amplification Dynamics of Platy-1 Retrotransposons in the Cebidae Platyrrhine Lineage

Platy-1 elements are Platyrrhine-specific, short interspersed elements originally discovered in the Callithrix jacchus (common marmoset) genome. To date, only the marmoset genome has been analyzed for Platy-1 repeat content. Here, we report full-length Platy-1 insertions in other New World monkey (NWM) genomes (Saimiri boliviensis, squirrel monkey; Cebus imitator, capuchin monkey; and Aotus nancymaae, owl monkey) and analyze the amplification dynamics of lineage-specific Platy-1 insertions. A relatively small number of full-length and lineage-specific Platy-1 elements were found in the squirrel, capuchin, and owl monkey genomes compared with the marmoset genome. In addition, only a few older Platy-1 subfamilies were recovered in this study, with no Platy-1 subfamilies younger than Platy-1-6. By contrast, 62 Platy-1 subfamilies were discovered in the marmoset genome. All of the lineage-specific insertions found in the squirrel and capuchin monkeys were fixed present. However, ∼15% of the lineage-specific Platy-1 loci in Aotus were polymorphic for insertion presence/absence. In addition, two new Platy-1 subfamilies were identified in the owl monkey genome with low nucleotide divergences compared with their respective consensus sequences, suggesting minimal ongoing retrotransposition in the Aotus genus and no current activity in the Saimiri, Cebus, and Sapajus genera. These comparative analyses highlight the finding that the high number of Platy-1 elements discovered in the marmoset genome is an exception among NWM analyzed thus far, rather than the rule. Future studies are needed to expand upon our knowledge of Platy-1 amplification in other NWM genomes.

ALUminating the Path of Atherosclerosis Progression: Chaos Theory Suggests a Role for Alu Repeats in the Development of Atherosclerotic Vascular Disease

Atherosclerosis (ATH) and coronary artery disease (CAD) are chronic inflammatory diseases with an important genetic background; they derive from the cumulative effect of multiple common risk alleles, most of which are located in genomic noncoding regions. These complex diseases behave as nonlinear dynamical systems that show a high dependence on their initial conditions; thus, long-term predictions of disease progression are unreliable. One likely possibility is that the nonlinear nature of ATH could be dependent on nonlinear correlations in the structure of the human genome. In this review, we show how chaos theory analysis has highlighted genomic regions that have shared specific structural constraints, which could have a role in ATH progression. These regions were shown to be enriched with repetitive sequences of the Alu family, genomic parasites that have colonized the human genome, which show a particular secondary structure and are involved in the regulation of gene expression. Here, we show the impact of Alu elements on the mechanisms that regulate gene expression, especially highlighting the molecular mechanisms via which the Alu elements alter the inflammatory response. We devote special attention to their relationship with the long noncoding RNA (lncRNA); antisense noncoding RNA in the INK4 locus (ANRIL), a risk factor for ATH; their role as microRNA (miRNA) sponges; and their ability to interfere with the regulatory circuitry of the (nuclear factor kappa B) NF-κB response. We aim to characterize ATH as a nonlinear dynamic system, in which small initial alterations in the expression of a number of repetitive elements are somehow amplified to reach phenotypic significance.

The angiotensin-I-converting enzyme insertion/deletion in polymorphic element codes for an AluYa5 RNA that downregulates gene expression

Angiotensin-I-converting enzyme (ACE) is involved in the synthesis and degradation of important bioactive peptides. The ACE gene has a 287-bp insertion/deletion polymorphism that controls ACE expression through a mechanism that remains elusive. In this study, we found that the 287-bp polymorphic element of the ACE gene, a member of the AluYa5 sub-family of Alu elements, codes for an RNA molecule that controls the levels of ACE mRNA. Transient transfection of a plasmid containing a CMV promoter upstream of the ACE polymorphic element resulted in significant expression of an AluYa5 RNA and reduced ACE mRNA expression as well as ACE enzymatic activity in AD 293 cells. The AluYa5 element also independently reduced the expression of other genes, regardless of whether these genes harbored Alu elements within their genomic context. Interestingly, the CMV promoter was not required for the expression of the AluYa5 element in AD 293 cells. The 287-bp sequence was sufficient to produce AluYa5 RNA and led to a significant reduction in ACE gene expression. Moreover, the removal of an 11-bp fragment of the 3' end of the ACE polymorphic sequence, which is specific to this particular AluYa5 element, did not prevent this element from being expressed but did affect its ability to target ACE expression. Thus, the expression of the AluYa5 polymorphic element within the ACE gene could explain why patients carrying the ACE insertion polymorphism have reduced risk of developing several chronic diseases.

Evidence for DNA Sequence Encoding of an Accessible Nucleosomal Array across Vertebrates

Nucleosome-depleted regions around which nucleosomes order following the "statistical" positioning scenario were recently shown to be encoded in the DNA sequence in human. This intrinsic nucleosomal ordering strongly correlates with oscillations in the local GC content as well as with the interspecies and intraspecies mutation profiles, revealing the existence of both positive and negative selection. In this letter, we show that these predicted nucleosome inhibitory energy barriers (NIEBs) with compacted neighboring nucleosomes are indeed ubiquitous to all vertebrates tested. These 1 kb-sized chromatin patterns are widely distributed along vertebrate chromosomes, overall covering more than a third of the genome. We have previously observed in human deviations from neutral evolution at these genome-wide distributed regions, which we interpreted as a possible indication of the selection of an open, accessible, and dynamic nucleosomal array to constitutively facilitate the epigenetic regulation of nuclear functions in a cell-type-specific manner. As a first, very appealing observation supporting this hypothesis, we report evidence of a strong association between NIEB borders and the poly(A) tails of Alu sequences in human. These results suggest that NIEBs provide adequate chromatin patterns favorable to the integration of Alu retrotransposons and, more generally to various transposable elements in the genomes of primates and other vertebrates.

Analysis of lineage-specific Alu subfamilies in the genome of the olive baboon, Papio anubis

Background

Alu elements are primate-specific retroposons that mobilize using the enzymatic machinery of L1 s. The recently completed baboon genome project found that the mobilization rate of Alu elements is higher than in the genome of any other primate studied thus far. However, the Alu subfamily structure present in and specific to baboons had not been examined yet.

Results

Here we report 129 Alu subfamilies that are propagating in the genome of the olive baboon, with 127 of these subfamilies being new and specific to the baboon lineage. We analyzed 233 Alu insertions in the genome of the olive baboon using locus specific polymerase chain reaction assays, covering 113 of the 129 subfamilies. The allele frequency data from these insertions show that none of the nine groups of subfamilies are nearing fixation in the lineage.

Conclusions

Many subfamilies of Alu elements are actively mobilizing throughout the baboon lineage, with most being specific to the baboon lineage.

Analysis of chromatin accessibility in decidualizing human endometrial stromal cells

Spontaneous decidualization of the endometrium in response to progesterone signaling is confined to menstruating species, including humans and other higher primates. During this process, endometrial stromal cells (EnSCs) differentiate into specialized decidual cells that control embryo implantation. We subjected undifferentiated and decidualizing human EnSCs to an assay for transposase accessible chromatin with sequencing (ATAC-seq) to map the underlying chromatin changes. A total of 185,084 open DNA loci were mapped accurately in EnSCs. Altered chromatin accessibility upon decidualization was strongly associated with differential gene expression. Analysis of 1533 opening and closing chromatin regions revealed over-representation of DNA binding motifs for known decidual transcription factors (TFs) and identified putative new regulators. ATAC-seq footprint analysis provided evidence of TF binding at specific motifs. One of the largest footprints involved the most enriched motif-basic leucine zipper-as part of a triple motif that also comprised the estrogen receptor and Pax domain binding sites. Without exception, triple motifs were located within Alu elements, which suggests a role for this primate-specific transposable element (TE) in the evolution of decidual genes. Although other TEs were generally under-represented in open chromatin of undifferentiated EnSCs, several classes contributed to the regulatory DNA landscape that underpins decidual gene expression.-Vrljicak, P., Lucas, E. S., Lansdowne, L., Lucciola, R., Muter, J., Dyer, N. P., Brosens, J. J., Ott, S. Analysis of chromatin accessibility in decidualizing human endometrial stromal cells.

Human transposable elements in Repbase: genomic footprints from fish to humans

Repbase is a comprehensive database of eukaryotic transposable elements (TEs) and repeat sequences, containing over 1300 human repeat sequences. Recent analyses of these repeat sequences have accumulated evidences for their contribution to human evolution through becoming functional elements, such as protein-coding regions or binding sites of transcriptional regulators. However, resolving the origins of repeat sequences is a challenge, due to their age, divergence, and degradation. Ancient repeats have been continuously classified as TEs by finding similar TEs from other organisms. Here, the most comprehensive picture of human repeat sequences is presented. The human genome contains traces of 10 clades (L1, CR1, L2, Crack, RTE, RTEX, R4, Vingi, Tx1 and Penelope) of non-long terminal repeat (non-LTR) retrotransposons (long interspersed elements, LINEs), 3 types (SINE1/7SL, SINE2/tRNA, and SINE3/5S) of short interspersed elements (SINEs), 1 composite retrotransposon (SVA) family, 5 classes (ERV1, ERV2, ERV3, Gypsy and DIRS) of LTR retrotransposons, and 12 superfamilies (Crypton, Ginger1, Harbinger, hAT, Helitron, Kolobok, Mariner, Merlin, MuDR, P, piggyBac and Transib) of DNA transposons. These TE footprints demonstrate an evolutionary continuum of the human genome.

Detection of de novo single nucleotide variants in offspring of atomic-bomb survivors close to the hypocenter by whole-genome sequencing

Ionizing radiation released by the atomic bombs at Hiroshima and Nagasaki, Japan, in 1945 caused many long-term illnesses, including increased risks of malignancies such as leukemia and solid tumours. Radiation has demonstrated genetic effects in animal models, leading to concerns over the potential hereditary effects of atomic bomb-related radiation. However, no direct analyses of whole DNA have yet been reported. We therefore investigated de novo variants in offspring of atomic-bomb survivors by whole-genome sequencing (WGS). We collected peripheral blood from three trios, each comprising a father (atomic-bomb survivor with acute radiation symptoms), a non-exposed mother, and their child, none of whom had any past history of haematological disorders. One trio of non-exposed individuals was included as a control. DNA was extracted and the numbers of de novo single nucleotide variants in the children were counted by WGS with sequencing confirmation. Gross structural variants were also analysed. Written informed consent was obtained from all participants prior to the study. There were 62, 81, and 42 de novo single nucleotide variants in the children of atomic-bomb survivors, compared with 48 in the control trio. There were no gross structural variants in any trio. These findings are in accord with previously published results that also showed no significant genetic effects of atomic-bomb radiation on second-generation survivors.

The Mobile Element Locator Tool (MELT): population-scale mobile element discovery and biology

Mobile element insertions (MEIs) represent ∼25% of all structural variants in human genomes. Moreover, when they disrupt genes, MEIs can influence human traits and diseases. Therefore, MEIs should be fully discovered along with other forms of genetic variation in whole genome sequencing (WGS) projects involving population genetics, human diseases, and clinical genomics. Here, we describe the Mobile Element Locator Tool (MELT), which was developed as part of the 1000 Genomes Project to perform MEI discovery on a population scale. Using both Illumina WGS data and simulations, we demonstrate that MELT outperforms existing MEI discovery tools in terms of speed, scalability, specificity, and sensitivity, while also detecting a broader spectrum of MEI-associated features. Several run modes were developed to perform MEI discovery on local and cloud systems. In addition to using MELT to discover MEIs in modern humans as part of the 1000 Genomes Project, we also used it to discover MEIs in chimpanzees and ancient (Neanderthal and Denisovan) hominids. We detected diverse patterns of MEI stratification across these populations that likely were caused by (1) diverse rates of MEI production from source elements, (2) diverse patterns of MEI inheritance, and (3) the introgression of ancient MEIs into modern human genomes. Overall, our study provides the most comprehensive map of MEIs to date spanning chimpanzees, ancient hominids, and modern humans and reveals new aspects of MEI biology in these lineages. We also demonstrate that MELT is a robust platform for MEI discovery and analysis in a variety of experimental settings.

The Influence of LINE-1 and SINE Retrotransposons on Mammalian Genomes

Transposable elements have had a profound impact on the structure and function of mammalian genomes. The retrotransposon Long INterspersed Element-1 (LINE-1 or L1), by virtue of its replicative mobilization mechanism, comprises ∼17% of the human genome. Although the vast majority of human LINE-1 sequences are inactive molecular fossils, an estimated 80-100 copies per individual retain the ability to mobilize by a process termed retrotransposition. Indeed, LINE-1 is the only active, autonomous retrotransposon in humans and its retrotransposition continues to generate both intra-individual and inter-individual genetic diversity. Here, we briefly review the types of transposable elements that reside in mammalian genomes. We will focus our discussion on LINE-1 retrotransposons and the non-autonomous Short INterspersed Elements (SINEs) that rely on the proteins encoded by LINE-1 for their mobilization. We review cases where LINE-1-mediated retrotransposition events have resulted in genetic disease and discuss how the characterization of these mutagenic insertions led to the identification of retrotransposition-competent LINE-1s in the human and mouse genomes. We then discuss how the integration of molecular genetic, biochemical, and modern genomic technologies have yielded insight into the mechanism of LINE-1 retrotransposition, the impact of LINE-1-mediated retrotransposition events on mammalian genomes, and the host cellular mechanisms that protect the genome from unabated LINE-1-mediated retrotransposition events. Throughout this review, we highlight unanswered questions in LINE-1 biology that provide exciting opportunities for future research. Clearly, much has been learned about LINE-1 and SINE biology since the publication of Mobile DNA II thirteen years ago. Future studies should continue to yield exciting discoveries about how these retrotransposons contribute to genetic diversity in mammalian genomes.

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.

Identification of a recently active mammalian SINE derived from ribosomal RNA

Complex eukaryotic genomes are riddled with repeated sequences whose derivation does not coincide with phylogenetic history and thus is often unknown. Among such sequences, the capacity for transcriptional activity coupled with the adaptive use of reverse transcription can lead to a diverse group of genomic elements across taxa, otherwise known as selfish elements or mobile elements. Short interspersed nuclear elements (SINEs) are nonautonomous mobile elements found in eukaryotic genomes, typically derived from cellular RNAs such as tRNAs, 7SL or 5S rRNA. Here, we identify and characterize a previously unknown SINE derived from the 3'-end of the large ribosomal subunit (LSU or 28S rDNA) and transcribed via RNA polymerase III. This new element, SINE28, is represented in low-copy numbers in the human reference genome assembly, wherein we have identified 27 discrete loci. Phylogenetic analysis indicates these elements have been transpositionally active within primate lineages as recently as 6 MYA while modern humans still carry transcriptionally active copies. Moreover, we have identified SINE28s in all currently available assembled mammalian genome sequences. Phylogenetic comparisons indicate that these elements are frequently rederived from the highly conserved LSU rRNA sequences in a lineage-specific manner. We propose that this element has not been previously recognized as a SINE given its high identity to the canonical LSU, and that SINE28 likely represents one of possibly many unidentified, active transposable elements within mammalian genomes.

Identification of novel exonic mobile element insertions in epithelial ovarian cancers

Mobile elements comprise about half of the human genome. Three active mobile element families (L1, Alu, and SVA) possibly cause diseases such as cancer. We conducted mobile element insertion (MEI) profiling of 44 epithelial ovarian cancers using exome-sequencing data. We identified a total of 106 MEIs using the Mobster program, 8 of which were novel exonic MEIs.

Dynamic Alu methylation during normal development, aging, and tumorigenesis

DNA methylation primarily occurs on CpG dinucleotides and plays an important role in transcriptional regulations during tissue development and cell differentiation. Over 25% of CpG dinucleotides in the human genome reside within Alu elements, the most abundant human repeats. The methylation of Alu elements is an important mechanism to suppress Alu transcription and subsequent retrotransposition. Decades of studies revealed that Alu methylation is highly dynamic during early development and aging. Recently, many environmental factors were shown to have a great impact on Alu methylation. In addition, aberrant Alu methylation has been documented to be an early event in many tumors and Alu methylation levels have been associated with tumor aggressiveness. The assessment of the Alu methylation has become an important approach for early diagnosis and/or prognosis of cancer. This review focuses on the dynamic Alu methylation during development, aging, and tumor genesis. The cause and consequence of Alu methylation changes will be discussed.

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.

Identification of three new Alu Yb subfamilies by source tracking of recently integrated Alu Yb elements

BACKGROUND

Alu elements are the most abundant mobile elements in the human genome, with over 1 million copies and constituting more than 10% of the genome. The majority of these Alu elements were inserted into the primate genome 35 to 60 million years ago, but certain subfamilies of Alu elements are relatively very new and suspected to be still evolving. We attempted to trace the source/master copies of all human-specific members of the Alu Yb lineage using a computational approach by clustering similar Yb elements and constructing an evolutionary relation among the members of a cluster.

RESULTS

We discovered that one copy of Yb8 at 10p14 is the source of several active Yb8 copies, which retrotransposed to generate 712 copies or 54% of all human-specific Yb8 elements. We detected eight other Yb8 elements that had generated ten or more copies, potentially acting as 'stealth drivers'. One Yb8 element at 14q32.31 seemed to act as the source copy for all Yb9 elements tested, having producing 13 active Yb9 elements, and subsequently generated a total of 131 full-length copies. We identified and characterized three new subclasses of Yb elements: Yb8a1, Yb10 and Yb11. Their copy numbers in the reference genome are 75, 8 and 16. We analysed personal genome data from the 1000 Genome Project and detected an additional 6 Yb8a1, 3 Yb10 and 15 Yb11 copies outside the reference genome. Our analysis indicates that the Yb8a1 subfamily has a similar age to Yb9 (1.93 million years and 2.15 million years, respectively), while Yb10 and Yb11 evolved only 1.4 and 0.71 million years ago, suggesting a linear evolutionary path from Yb8a1 to Yb10 and then to Yb11. Our preliminary data indicate that members in Yb10 and Yb11 are mostly polymorphic, indicating their young age.

CONCLUSIONS

Our findings suggest that the Yb lineage is still evolving with new subfamilies being formed. Due to their very young age and the high rate of being polymorphic, insertions from these young subfamilies are very useful genetic markers for studying human population genetics and migration patterns, and the trend for mobile element insertions in the human genome.

Detecting Alu insertions from high-throughput sequencing data

High-throughput sequencing technologies have allowed for the cataloguing of variation in personal human genomes. In this manuscript, we present alu-detect, a tool that combines read-pair and split-read information to detect novel Alus and their precise breakpoints directly from either whole-genome or whole-exome sequencing data while also identifying insertions directly in the vicinity of existing Alus. To set the parameters of our method, we use simulation of a faux reference, which allows us to compute the precision and recall of various parameter settings using real sequencing data. Applying our method to 100 bp paired Illumina data from seven individuals, including two trios, we detected on average 1519 novel Alus per sample. Based on the faux-reference simulation, we estimate that our method has 97% precision and 85% recall. We identify 808 novel Alus not previously described in other studies. We also demonstrate the use of alu-detect to study the local sequence and global location preferences for novel Alu insertions.

Detection of serum Alu element hypomethylation for the diagnosis and prognosis of glioma

Global genomic hypomethylation is a hallmark of cancer in humans. In the present study, the feasibility of measuring hypomethylation of Alu elements (Alu) in serum and its clinical utility were investigated. Tumor tissues and matched serum specimens from 65 glioma patients and serum samples from 30 healthy controls were examined for Alu hypomethylation by bisulfite sequencing. The median serum Alu methylation level was 47.30 % in patients (interquartile range (IQR), 35.40-54.25 %) and 57.90 % in the controls (IQR, 55.25-61.45 %). The median Alu methylation level in tumor samples was 40.30 % (IQR, 36.80-54.20 %), which shows the correlation of Alu hypomethylation between tumor and serum samples (r = 0.882) in the study group. The methylation level was higher in the low-grade glioma group than in the high-grade group both in tumor and serum samples. A correlation between high methylation level and longer survival time was detected in tumor and serum samples. Receiver operating characteristic curve analysis showed that the area under the curve for diagnosis was 0.861 (95 % confidence interval, 0.789-0.933), suggesting that Alu hypomethylation in serum may be of diagnostic value. Our results indicate that the detection of Alu hypomethylation in serum may be clinically useful for the diagnosis and prognosis of glioma.

Alu mobile elements: from junk DNA to genomic gems

Alus, the short interspersed repeated sequences (SINEs), are retrotransposons that litter the human genomes and have long been considered junk DNA. However, recent findings that these mobile elements are transcribed, both as distinct RNA polymerase III transcripts and as a part of RNA polymerase II transcripts, suggest biological functions and refute the notion that Alus are biologically unimportant. Indeed, Alu RNAs have been shown to control mRNA processing at several levels, to have complex regulatory functions such as transcriptional repression and modulating alternative splicing and to cause a host of human genetic diseases. Alu RNAs embedded in Pol II transcripts can promote evolution and proteome diversity, which further indicates that these mobile retroelements are in fact genomic gems rather than genomic junks.

Genomic characterization of two large Alu-mediated rearrangements of the BRCA1 gene

To determine whether a large genomic rearrangement is actually novel and to gain insight about the mutational mechanism responsible for its occurrence, molecular characterization with breakpoint identification is mandatory. We here report the characterization of two large deletions involving the BRCA1 gene. The first rearrangement harbored a 89,664-bp deletion comprising exon 7 of the BRCA1 gene to exon 11 of the NBR1 gene (c.441+1724_oNBR1:c.1073+480del). Two highly homologous Alu elements were found in the genomic sequences flanking the deletion breakpoints. Furthermore, a 20-bp overlapping sequence at the breakpoint junction was observed, suggesting that the most likely mechanism for the occurrence of this rearrangement was nonallelic homologous recombination. The second rearrangement fully characterized at the nucleotide level was a BRCA1 exons 11-15 deletion (c.671-319_4677-578delinsAlu). The case harbored a 23,363-bp deletion with an Alu element inserted at the breakpoints of the deleted region. As the Alu element inserted belongs to a still active AluY family, the observed rearrangement could be due to an insertion-mediated deletion mechanism caused by Alu retrotransposition. To conclude, we describe the breakpoints of two novel large deletions involving the BRCA1 gene and analysis of their genomic context allowed us to gain insight about the respective mutational mechanism.