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

Human Retrotransposons and Effective Computational Detection Methods for Next-Generation Sequencing Data

Transposable elements (TEs) are classified into two classes according to their mobilization mechanism. Compared to DNA transposons that move by the "cut and paste" mechanism, retrotransposons mobilize via the "copy and paste" method. They have been an essential research topic because some of the active elements, such as Long interspersed element 1 (LINE-1), Alu, and SVA elements, have contributed to the genetic diversity of primates beyond humans. In addition, they can cause genetic disorders by altering gene expression and generating structural variations (SVs). The development and rapid technological advances in next-generation sequencing (NGS) have led to new perspectives on detecting retrotransposon-mediated SVs, especially insertions. Moreover, various computational methods have been developed based on NGS data to precisely detect the insertions and deletions in the human genome. Therefore, this review discusses details about the recently studied and utilized NGS technologies and the effective computational approaches for discovering retrotransposons through it. The final part covers a diverse range of computational methods for detecting retrotransposon insertions with human NGS data. This review will give researchers insights into understanding the TEs and how to investigate them and find connections with research interests.

A comprehensive analysis of gorilla-specific LINE-1 retrotransposons

BACKGROUND

Long interspersed element-1 (LINE-1 or L1) is the most abundant retrotransposons in the primate genome. They have approximately 520,000 copies and make up ~ 17% of the primate genome. Full-length L1s can mobilize to a new genomic location using their enzymatic machinery. Gorilla is the second closest species to humans after the chimpanzee, and human-gorilla split 7-12 million years ago. The gorilla genome provides an opportunity to explore primate origins and evolution.

OBJECTIVE

L1s have contributed to genome diversity and variations during primate evolution. This study aimed to identify gorilla-specific L1s using a more recent version of the gorilla reference genome (Mar. 2016 GSMRT3/gorGor5).

METHODS

We collected gorilla-specific L1 candidates through computational analysis and manual inspection. L1Xplorer was used to identify whether full-length gorilla-specific L1s were intact. In addition, to determine the level of sequence conservation between intact fulllength gorilla-specific L1s, two ORFs of intact L1s were aligned with the L1PA2 consensus sequence.

RESULTS

2002 gorilla-specific L1 candidates were identified through computational analysis. Among them, we manually inspected 1,883 gorilla-specific L1s, among which most of them belong to the L1PA2 subfamily and 12 were intact L1s that could influence genomic variations in the gorilla genome. Interestingly, the 12 intact full-length gorilla-specific L1s have 14 highly conserved nonsynonymous mutations, including 6 mutations and 8 mutations in ORF1 and ORF2, respectively. In comparison to the intact full-length chimpanzee-specific L1s and human-specific hot-L1s, two of these in ORF1 (L256F and E293G) were shown as gorilla-specific nonsynonymous mutations.

CONCLUSION

The gorilla-specific L1s may have had significantly affected the gorilla genome to compose a genome different form that of other primates during primate evolution.

Higher Rates of Processed Pseudogene Acquisition in Humans and Three Great Apes Revealed by Long-Read Assemblies

LINE-1-mediated retrotransposition of protein-coding mRNAs is an active process in modern humans for both germline and somatic genomes. Prior works that surveyed human data mostly relied on detecting discordant mappings of paired-end short reads, or exon junctions contained in short reads. Moreover, there have been few genome-wide comparisons between gene retrocopies in great apes and humans. In this study, we introduced a more sensitive and accurate method to identify processed pseudogenes. Our method utilizes long-read assemblies, and more importantly, is able to provide full-length retrocopy sequences as well as flanking regions which are missed by short-read based methods. From 22 human individuals, we pinpointed 40 processed pseudogenes that are not present in the human reference genome GRCh38 and identified 17 pseudogenes that are in GRCh38 but absent from some input individuals. This represents a significantly higher discovery rate than previous reports (39 pseudogenes not in the reference genome out of 939 individuals). We also provided an overview of lineage-specific retrocopies in chimpanzee, gorilla, and orangutan genomes.

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.

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.

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.

Recently integrated Alu insertions in the squirrel monkey (Saimiri) lineage and application for population analyses

BACKGROUND

The evolution of Alu elements has been ongoing in primate lineages and Alu insertion polymorphisms are widely used in phylogenetic and population genetics studies. Alu subfamilies in the squirrel monkey (Saimiri), a New World Monkey (NWM), were recently reported. Squirrel monkeys are commonly used in biomedical research and often require species identification. The purpose of this study was two-fold: 1) Perform locus-specific PCR analyses on recently integrated Alu insertions in Saimiri to determine their amplification dynamics, and 2) Identify a subset of Alu insertion polymorphisms with species informative allele frequency distributions between the Saimiri sciureus and Saimiri boliviensis groups.

RESULTS

PCR analyses were performed on a DNA panel of 32 squirrel monkey individuals for 382 Alu insertion events ≤2% diverged from 46 different Alu subfamily consensus sequences, 25 Saimiri specific and 21 NWM specific Alu subfamilies. Of the 382 loci, 110 were polymorphic for presence / absence among squirrel monkey individuals, 35 elements from 14 different Saimiri specific Alu subfamilies and 75 elements from 19 different NWM specific Alu subfamilies (13 of 46 subfamilies analyzed did not contain polymorphic insertions). Of the 110 Alu insertion polymorphisms, 51 had species informative allele frequency distributions between Saimiri sciureus and Saimiri boliviensis groups.

CONCLUSIONS

This study confirms the evolution of Alu subfamilies in Saimiri and provides evidence for an ongoing and prolific expansion of these elements in Saimiri with many active subfamilies concurrently propagating. The subset of polymorphic Alu insertions with species informative allele frequency distribution between Saimiri sciureus and Saimiri boliviensis will be instructive for specimen identification and conservation biology.

Romain CP, Bennett EC, Robichaux A, Clement BN, Raveendran M, Worley KC, Phillips-Conroy J, Jolly CJ, Rogers J, Konkel MK, Batzer MA. Alu Insertion Polymorphisms as Evidence for Population Structure in Baboons

Male dispersal from the natal group at or near maturity is a feature of most baboon (Papio) species. It potentially has profound effects upon population structure and evolutionary processes, but dispersal, especially for unusually long distances, is not readily documented by direct field observation. In this pilot study, we investigate the possibility of retrieving baboon population structure in yellow (Papio cynocephalus) and kinda (Papio kindae) baboons from the distribution of variation in a genome-wide set of 494 Alu insertion polymorphisms, made available via the recently completed Baboon Genome Analysis Consortium. Alu insertion variation in a mixed population derived from yellow and olive (Papio anubis) baboons identified each individual's proportion of heritage from either parental species. In an unmixed yellow baboon population, our analysis showed greater similarity between neighboring than between more distantly situated groups, suggesting structuring of the population by male dispersal distance. Finally (and very provisionally), an unexpectedly sharp difference in Alu insertion frequencies between members of neighboring social groups of kinda baboons suggests that intergroup migration may be more rare than predicted in this little known species.