Genome-wide targeted search for human specific and polymorphic L1 integrations

Retroelements (REs) occupy up to 40% of the human genome. Newly integrated REs can change the pattern of expression of pre-existing host genes and therefore might play a significant role in evolution. In particular, human- and primate-specific REs could affect the divergence of the Hominoidea superfamily. A comparative genome-wide analysis of RE sites of integration, neighboring genes, and their regulatory interplay in human and ape genomes would be of help in understanding the impact of REs on evolution and genome regulation. We have developed a technique for the genome-wide comparison of the integrations of transposable elements in genomic DNAs of closely related species. The technique called targeted genome differences analysis (TGDA) is also useful for the detection of deletion/insertion polymorphisms of REs. The technique is based on an enhanced version of subtractive hybridization and does not require preliminary knowledge of the genome sequences under comparison. In this report, we describe its application to the detection and analysis of human specific L1 integrations and their polymorphisms. We obtained a library highly enriched in human-specific L1 insertions and identified 24 such new insertions. Many of these insertions are polymorphic in human populations. The total number of human-specific L1 inserts was estimated to be approximately 4000. The results suggest that TGDA is a universal method that can be successfully used for the detection of evolutionary and polymorphic markers in any closely related genomes.     

A novel approach to identification of somatic retroelements’ insertions in human genome

The activity of retroelements is one of the factors leading to genetic variability of the modern humans. Insertions of retroelements may result in alteration of gene expression and functional diversity between cells. In recent years an increasing amount of data indicating an elevated level of retroelements’ mobilisation in some human and animal tissues has been reported. Therefore, the development of a system for the detection of somatic retroposition events is required. Here we describe a novel approach to the whole-genome identification of somatic retroelement insertions in human genome. The developed approach was applied for the comparisons of somatic mosaicism levels in two tissues of the investigated individual. A total of 3410 insertions of retroelements belonging to AluYa5 subfamily were identified.     

A new database on polymorphic retroelements in human genome (PRED)

Comparison of primate genomes sequences has confirmed the evidence that substantial part of intra- and interspecies differences is provided by retroelements. Human genome contains thousands of polymorphic retroelement copies considered to be perspective molecular genetic markers of new generation. However utilization of polymorphic retroelements as molecular genetic markers is limited due to lack of systematic data on their number, genomic context and distribution among human populations. We have created first bilingual (Russian/English) internet-resource devoted to known polymorphic retroelements discovered in human genome by our group as well as by other researchers worldwide. The database contains information about each retroelement copy location, position relative to known and predicted genes, frequency of alleles in human populations and others. Our internet portal allows to perform a search in database using multiple search conditions and available on http://labcfg.ibch.ru/home.html. The database provides an opportunity to investigate distribution of polymorphic retroelements in human genome and to design new genetic markers for various population and medical studies.     

Retroposons in modern human genome evolution

The ascertainment of the rates and driving forces of human genome evolution along with the genetic diversity of populations or separate population groups remains a topical problem of fundamental and applied genomics. According to the results of comparative analysis, the most numerous human genome structure peculiarities are connected with the distribution of mobile genetic retroelements—LTR, LINE1, SVA, and Alu repeats. Due to the wide distribution in different genome loci, conversed retropositional activity, and the retroelements’ regulatory potential, let us regard them as one of the significant evolutionary driving forces and the source of human genome variability. In the current review, we summarize published data and recent results of our research aimed at the analysis of the evolutionary impact of the young retroelements group on the function and variability of the human genome. We examine modern approaches of the polygenomic identification of polymorphic retroelements inserts. Using an original Internet resource, we analyze special features of the genomic polymorphic inserts of AluY repeats. We thoroughly characterize the strategy of large-scale functional analysis of polymorphic retroelement inserts. The presented results confirm the hypothesis of the roles of retroelements as active cis regulatory elements that are able to modulate surrounding genes.     

LTR retrotransposon dynamics in the evolution of the olive (Olea europaea) genome

Improved knowledge of genome composition, especially of its repetitive component, generates important information for both theoretical and applied research. The olive repetitive component is made up of two main classes of sequences: tandem repeats and retrotransposons (REs). In this study, we provide characterization of a sample of 254 unique full-length long terminal repeat (LTR) REs. In the sample, Ty1-Copia elements were more numerous than Ty3-Gypsy elements. Mapping a large set of Illumina whole-genome shotgun reads onto the identified retroelement set revealed that Gypsy elements are more redundant than Copia elements. The insertion time of intact retroelements was estimated based on sister LTR’s divergence. Although some elements inserted relatively recently, the mean insertion age of the isolated retroelements is around 18 million yrs. Gypsy and Copia retroelements showed different waves of transposition, with Gypsy elements especially active between 10 and 25 million yrs ago and nearly inactive in the last 7 million yrs. The occurrence of numerous solo-LTRs related to isolated full-length retroelements was ascertained for two Gypsy elements and one Copia element. Overall, the results reported in this study show that RE activity (both retrotransposition and DNA loss) has impacted the olive genome structure in more ancient times than in other angiosperms.     

A new polymorphic retroelement database (PRED) for the human genome

Comparison of primate genomic sequences has demonstrated that the intra-and interspecific genetic variation is provided by retroelements (REs). The human genome contains many thousands of polymorphic RE copies, which are regarded as a promising source of new generation molecular genetic markers. However, the absence of systematized data on the RE number, distribution, genomic context, and abundance in various human populations limits the use of RE insertion polymorphism. We designed the first bilingual (Russian/English) web resource on the known polymorphic REs discovered both by our team and other researchers. The database contains the information about the genomic location of each RE, its position relative to known and predicted genes, abundance in human populations, and other data. Our web portal () allows a search of the database with user-specified parameters. The database makes it possible to most comprehensively analyze the RE distribution in the human genome and to design molecular genetic markers for studies of human genome diversity and biomedical applications.     

Large-scale determination of the methylation status of retrotransposons in different tissues using a methylation tags approach

A technique for simultaneous determination of the methylation status of numerous loci containing retroelements (REs) is reported. It is based on the observation that methylated and unmethylated areas in the genome are usually extended, and therefore the methylation of particular methyl-sensitive restriction endonuclease recognition sites might reflect the methylation status of DNA regions around them. The method includes dot-blot hybridization of repeat flanking sequences arrayed on a solid support with specifically amplified flanking regions of presumably unmethylated repeats. A multitude of flanking regions of REs adjacent to unmethylated restriction sites are amplified simultaneously, providing a complex hybridization probe. The technique thus allows the determination of the methylation status of restriction sites, which serve as tags of the methylation status of the surrounding regions. The validity of the technique was confirmed by various means, including bisulfite sequencing. The technique was successfully applied to the identification of methylation patterns of the regions surrounding 38 human-specific HERV-K(HML-2) long terminal repeats in cerebellum- and lymph node-derived genomic DNAs. The described technique can be readily adapted to the use of DNA microarray technology.     

A tissue-specific decrease in the pre-mRNA level of L1-and Alu-containing alleles of human genes

LINE1 and Alu retroelements occupy approximately 17 and 13% of the human genome, respectively. They include the evolutionarily youngest element groups Ta-L1, AluYa5, and AluYb8, many inserts of which are polymorphic in the Homo sapiens population. Despite the data on the ability of L1 and Alu elements to cause various modifications of the genome, the effects of these retroelements on gene expression have yet not been studied. Using the RT PCR method, we analyzed the pre-mRNA (heterogeneous nuclear RNA) content of allele pairs of four genes in five human cell lines, heterozygous with respect to intronic inserts of L1 and Alu elements. We showed for the first time a tissue-specific decrease in the pre-mRNA content of the gene allele bearing L1 or Alu inserts relative to the other allele of the same gene lacking the retroelement.     

Sequencing,identification and mapping of primed L1 elements (SIMPLE) reveals significant variation in full length L1 elements between individuals

Background

There are over a half a million copies of L1 retroelements in the human genome which are responsible for as much as 0.5% of new human genetic diseases. Most new L1 inserts arise from young source elements that are polymorphic in the human genome. Highly active polymorphic “hot” L1 source elements have been shown to be capable of extremely high levels of mobilization and result in numerous instances of disease. Additionally, hot polymorphic L1s have been described to be highly active within numerous cancer genomes. These hot L1s result in mutagenesis by insertion of new L1 copies elsewhere in the genome, but also have been shown to generate additional full length L1 insertions which are also hot and able to further retrotranspose. Through this mechanism, hot L1s may amplify within a tumor and result in a continued cycle of mutagenesis.

Results and conclusions

We have developed a method to detect full-length, polymorphic L1 elements using a targeted next generation sequencing approach, Sequencing Identification and Mapping of Primed L1 Elements (SIMPLE). SIMPLE has 94% sensitivity and detects nearly all full-length L1 elements in a genome. SIMPLE will allow researchers to identify hot mutagenic full-length L1s as potential drivers of genome instability. Using SIMPLE we find that the typical individual has approximately 100 non-reference, polymorphic L1 elements in their genome. These elements are at relatively low population frequencies relative to previously identified polymorphic L1 elements and demonstrate the tremendous diversity in potentially active L1 elements in the human population.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1374-y) contains supplementary material, which is available to authorized users.     

Methylation perturbations in retroelements within the genome of a Mus interspecific hybrid correlate with double minute chromosome formation

A reduction in the DNA modification of cytosine methylation has been linked directly to chromosome rearrangements concomitant with retroelement amplification in several marsupial hybrid genomes. While phenotypes observed for interspecific eutherian hybrids are suggestive of methylation perturbations and retroelement instability, no link between retroelements, DNA methylation, and chromosome instability has yet been identified. Previous studies in eutherian hybrids, however, have been limited to a gross examination of methylation using methylation-sensitive restriction enzyme analysis or focused on single-copy genes and/or have avoided examination of repetitive DNA. Methylation changes and retroelements are proposed as mechanisms for double minute chromosome formation and oncogene amplification, both present in the genome of a Mus hybrid model, thus making it an ideal system to evaluate methylation status more closely. We have used the PCR-based methodologies methylation-sensitive amplicon subtraction (MS-AS) and methylation-sensitive representational difference analysis (MS-RDA) to detect differentially methylated sequences between three complex genomes and to isolate methylation perturbations in a Mus musculusxMus caroli hybrid. This novel application of MS-AS and MS-RDA resulted in the isolation of differentially methylated retroelements surrounding the locus on Chromosome 10 responsible for double minute chromosome formation within this interspecific eutherian hybrid.     

The genomic landscape of polymorphic human nuclear mitochondrial insertions

The transfer of mitochondrial genetic material into the nuclear genomes of eukaryotes is a well-established phenomenon that has been previously limited to the study of static reference genomes. The recent advancement of high throughput sequencing has enabled an expanded exploration into the diversity of polymorphic nuclear mitochondrial insertions (NumtS) within human populations. We have developed an approach to discover and genotype novel Numt insertions using whole genome, paired-end sequencing data. We have applied this method to a thousand individuals in 20 populations from the 1000 Genomes Project and other datasets and identified 141 new sites of Numt insertions, extending our current knowledge of existing NumtS by almost 20%. We find that recent Numt insertions are derived from throughout the mitochondrial genome, including the D-loop, and have integration biases that differ in some respects from previous studies on older, fixed NumtS in the reference genome. We determined the complete inserted sequence for a subset of these events and have identified a number of nearly full-length mitochondrial genome insertions into nuclear chromosomes. We further define their age and origin of insertion and present an analysis of their potential impact to ongoing studies of mitochondrial heteroplasmy and disease.     

REMiner: a tool for unbiased mining and analysis of repetitive elements and their arrangement structures of large chromosomes

Repetitive elements (REs) constitute a substantial portion of the genomes of human and other species; however, the RE profiles (type, density, and arrangement) within the individual genomes have not been fully characterized. In this study, we developed an RE analysis tool, called REMiner, for a chromosome-wide investigation into the occurrence of individual REs and arrangement of clusters of REs, and REMiner's functional features were examined using the human chromosome Y. The algorithm implemented by REMiner focused on unbiased mining of REs in large chromosomes and data interface within a viewer. The data from the chromosome demonstrated that REMiner is an efficient tool in regard to its capacity for a large query size and the availability of a high-resolution viewer, featuring instant retrieval of alignment data and control of magnification and identity ratio. The chromosome-wide survey identified a diverse population of ordered RE arrangements, which may participate in the genome biology.     

How repeated retroelements format genome function

Genomes operate as sophisticated information storage systems. Generic repeated signals in the DNA format expression of coding sequence files and organize additional functions essential for genome replication and accurate transmission to progeny cells. Retroelements comprise a major fraction of many genomes and contain a surprising diversity of functional signals. In this article, we summarize some features of the taxonomic distribution of retroelements, especially mammalian SINEs, tabulate functional roles documented for different classes of retroelements, and discuss their potential roles as genome organizers. In particular, the fact that certain retroelements serve as boundaries for heterochromatin domains and provide a significant fraction of scaffolding/matrix attachment regions (S/MARs) suggests that the reversed transcribed component of the genome plays a major architectonic role in higher order physical structuring. Employing an information science model, the "functionalist" perspective on repetitive DNA leads to new ways of thinking about the systemic organization of cellular genomes and provides several novel possibilities involving retroelements in evolutionarily significant genome reorganization.     

Retrofitting the genome: L1 extinction follows endogenous retroviral expansion in a group of muroid rodents

Long interspersed nuclear element 1 (LINE-1; L1) retrotransposons are the most common retroelements in mammalian genomes. Unlike individual families of endogenous retroviruses (ERVs), they have remained active throughout the mammalian radiation and are responsible for most of the retroelement movement and much genome rearrangement within mammals. They can be viewed as occupying a substantial niche within mammalian genomes. Our previous demonstration that L1s and B1 short interspersed nuclear elements (SINEs) are inactive in a group of South American rodents led us to ask if other elements have amplified to fill the empty niche. We identified a novel and highly active family of ERVs (mysTR). To determine whether loss of L1 activity was correlated with expansion of mysTR, we examined mysTR activity in four South American rodent species that have lost L1 and B1 activity and four sister species with active L1s. The copy number of recent mysTR insertions was extremely high, with an average of 4,200 copies per genome. High copy numbers exist in both L1-active and L1-extinct species, so the mysTR expansion appears to have preceded the loss of both SINE and L1 activity rather than to have filled an empty niche created by their loss. It may be coincidental that two unusual genomic events--loss of L1 activity and massive expansion of an ERV family--occur in the same group of mammals. Alternatively, it is possible that this large ERV expansion set the stage for L1 extinction.     

Rapid genome evolution revealed by comparative sequence analysis of orthologous regions from four triticeae genomes

Bread wheat (Triticum aestivum) is an allohexaploid species, consisting of three subgenomes (A, B, and D). To study the molecular evolution of these closely related genomes, we compared the sequence of a 307-kb physical contig covering the high molecular weight (HMW)-glutenin locus from the A genome of durum wheat (Triticum turgidum, AABB) with the orthologous regions from the B genome of the same wheat and the D genome of the diploid wheat Aegilops tauschii (Anderson et al., 2003; Kong et al., 2004). Although gene colinearity appears to be retained, four out of six genes including the two paralogous HMW-glutenin genes are disrupted in the orthologous region of the A genome. Mechanisms involved in gene disruption in the A genome include retroelement insertions, sequence deletions, and mutations causing in-frame stop codons in the coding sequences. Comparative sequence analysis also revealed that sequences in the colinear intergenic regions of these different genomes were generally not conserved. The rapid genome evolution in these regions is attributable mainly to the large number of retrotransposon insertions that occurred after the divergence of the three wheat genomes. Our comparative studies indicate that the B genome diverged prior to the separation of the A and D genomes. Furthermore, sequence comparison of two distinct types of allelic variations at the HMW-glutenin loci in the A genomes of different hexaploid wheat cultivars with the A genome locus of durum wheat indicates that hexaploid wheat may have more than one tetraploid ancestor.     

Sequence analysis of p53 response-elements suggests multiple binding modes of the p53 tetramer to DNA targets

The p53 tetramer recognizes specifically a 20-bp DNA element. Here, we examined symmetries encoded in p53 response elements (p53REs). We analyzed base inversion correlations within the half-site, as well as in the full-site palindrome. We found that p53REs are not only direct repeats of half-sites; rather, two p53 half-sites couple to form a higher order 20bp palindrome. The palindrome couplings between the half-sites are stronger for the human than for the mouse genome. The full-site palindrome and half-site palindrome are controlled by insertions between the two half-sites. The most notable feature is that the full-site palindrome with coupling between quarter-sites one and four (H14 coupling) dominates the p53REs without insertions. The most frequently observed insertion in human p53REs of 3bp enhances the half-site palindrome. The statistical frequencies of the coupling between the half-sites in the human genome correlate with grouped experimental p53 affinities with p53REs. Examination of known p53REs indicates the H14 couplings are stronger for positive regulation than for negatively regulated p53REs, with repressors having the lowest H14 couplings. We propose that the palindromic sequence couplings may encode such potential preferred multiple binding modes of the p53 tetramer to DNA.