Lineage specific evolution of the VNTR composite retrotransposon central domain and its role in retrotransposition of gibbon LAVA elements

Background

VNTR (Variable Number of Tandem Repeats) composite retrotransposons - SVA (SINE-R-VNTR-Alu), LAVA (LINE-1-Alu-VNTR-Alu), PVA (PTGR2-VNTR-Alu) and FVA (FRAM-VNTR-Alu) - are specific to hominoid primates. Their assembly, the evolution of their 5’ and 3’ domains, and the functional significance of the shared 5’ Alu-like region are well understood. The central VNTR domain, by contrast, has long been assumed to represent a more or less random collection of 30-50 bp GC-rich repeats. It is only recently that it attracted attention in the context of regulation of SVA expression.

Results

Here we provide evidence that the organization of the VNTR is non-random, with conserved repeat unit (RU) arrays at both the 5’ and 3’ ends of the VNTRs of human, chimpanzee and orangutan SVA and gibbon LAVA. The younger SVA subfamilies harbour highly organized internal RU arrays. The composition of these arrays is specific to the human/chimpanzee and orangutan lineages, respectively. Tracing the development of the VNTR through evolution we show for the first time how tandem repeats evolve within the constraints set by a functional, non-autonomous non-LTR retrotransposon in two different families - LAVA and SVA - in different hominoid lineages. Our analysis revealed that a microhomology-driven mechanism mediates expansion/contraction of the VNTR domain at the DNA level.Elements of all four VNTR composite families have been shown to be mobilized by the autonomous LINE1 retrotransposon in trans. In case of SVA, key determinants of mobilization are found in the 5’ hexameric repeat/Alu-like region. We now demonstrate that in LAVA, by contrast, the VNTR domain determines mobilization efficiency in the context of domain swaps between active and inactive elements.

Conclusions

The central domain of VNTR composites evolves in a lineage-specific manner which gives rise to distinct structures in gibbon LAVA, orangutan SVA, and human/chimpanzee SVA. The differences observed between the families and lineages are likely to have an influence on the expression and mobilization of the elements.

Electronic supplementary material

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

Phylogenetic Analysis of a Retroposon Family in African Great Apes

The SINE-R retroposon family has been identified by its relationship with the long terminal repeats (LTRs) of human endogenous retrovirus class K (HERV-K) as a mobile element that has evolved recently in the human genome. Here we examined the recent evolutionary history of this class of elements by a PCR approach to genomic DNA from the African great apes and by phylogenetic analysis including comparison with the HERV K10 parent sequence. With primers derived from a cDNA sequence from human brain, we identified 27 sequences from the chimpanzee and 16 from the gorilla. Phylogenetic comparisons with previously recognized sequences from the human and from the orangutan and gibbon revealed wide overlap of elements across species, suggesting multiple origins in the course of hominoid evolution. Two human elements SINE-R.C2 and HS307 were the furthest removed from the HERV-K10 sequence but these two elements were closely related to three elements from the chimpanzee and four elements from the gorilla. This group of elements (our clusters 14 and 15) appears to have transposed late in hominoid evolution. One element (Ch-M16) showed 99.1% sequence identity with the SINE-R.C2 element, which is human-specific. Thus the SINE-R family appears to have continued to be active in transposition throughout the course of primate evolution. Received: 12 March 1999 / Accepted: 25 May 1999     

SVA elements: a hominid-specific retroposon family

SVA is a composite repetitive element named after its main components, SINE, VNTR and Alu. We have identified 2762 SVA elements from the human genome draft sequence. Genomic distribution analysis indicates that the SVA elements are enriched in G+C-rich regions but have no preferences for inter- or intragenic regions. A phylogenetic analysis of the elements resulted in the recovery of six subfamilies that were named SVA_A to SVA_F. The composition, age and genomic distribution of the subfamilies have been examined. Subfamily age estimates based upon nucleotide divergence indicate that the expansion of four SVA subfamilies (SVA_A, SVA_B, SVA_C and SVA_D) began before the divergence of human, chimpanzee and gorilla, while subfamilies SVA_E and SVA_F are restricted to the human lineage. A survey of human genomic diversity associated with SVA_E and SVA_F subfamily members showed insertion polymorphism frequencies of 37.5% and 27.6%, respectively. In addition, we examined the amplification dynamics of SVA elements throughout the primate order and traced their origin back to the beginnings of hominid primate evolution, approximately 18 to 25 million years ago. This makes SVA elements the youngest family of retroposons in the primate order.     

The evolution of mobile DNAs: when will transposons create phylogenies that look as if there is a master gene?

Some families of mammalian interspersed repetitive DNA, such as the Alu SINE sequence, appear to have evolved by the serial replacement of one active sequence with another, consistent with there being a single source of transposition: the "master gene." Alternative models, in which multiple source sequences are simultaneously active, have been called "transposon models." Transposon models differ in the proportion of elements that are active and in whether inactivation occurs at the moment of transposition or later. Here we examine the predictions of various types of transposon model regarding the patterns of sequence variation expected at an equilibrium between transposition, inactivation, and deletion. Under the master gene model, all bifurcations in the true tree of elements occur in a single lineage. We show that this property will also hold approximately for transposon models in which most elements are inactive and where at least some of the inactivation events occur after transposition. Such tree shapes are therefore not conclusive evidence for a single source of transposition.     

Evolutionary Breakpoints in the Gibbon Suggest Association between Cytosine Methylation and Karyotype Evolution

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

Recently mobilized transposons in the human and chimpanzee genomes

Transposable genetic elements are abundant in the genomes of most organisms, including humans. These endogenous mutagens can alter genes, promote genomic rearrangements, and may help to drive the speciation of organisms. In this study, we identified almost 11,000 transposon copies that are differentially present in the human and chimpanzee genomes. Most of these transposon copies were mobilized after the existence of a common ancestor of humans and chimpanzees, approximately 6 million years ago. Alu, L1, and SVA insertions accounted for >95% of the insertions in both species. Our data indicate that humans have supported higher levels of transposition than have chimpanzees during the past several million years and have amplified different transposon subfamilies. In both species, approximately 34% of the insertions were located within known genes. These insertions represent a form of species-specific genetic variation that may have contributed to the differential evolution of humans and chimpanzees. In addition to providing an initial overview of recently mobilized elements, our collections will be useful for assessing the impact of these insertions on their hosts and for studying the transposition mechanisms of these elements.     

Construction, characterization, and chromosomal mapping of a fosmid library of the white-cheeked gibbon (Nomascus leucogenys)

Gibbons have experienced extensive karyotype rearrangements during evolution and represent an ideal model for studying the underlying molecular mechanism of evolutionary chromosomal rearrangements. It is anticipated that the cloning and sequence characterization of evolutionary chromosomal breakpoints will provide vital insights into the molecular force that has driven such a radical karyotype reshuffle in gibbons. We constructed and characterized a high-quality fosmid library of the white-cheeked gibbon (Nomascus leucogenys) containing 192,000 non- redundant clones with an average insert size of 38 kb and 2.5-fold genome coverage. By end sequencing of 100 randomly selected fosmid clones, we generated 196 sequence tags for the library. These end-sequenced fosmid clones were then mapped onto the chromosomes of the white-cheeked gibbon by fluorescence in situ hybridization, and no spurious chimeric clone was detected. BLAST search against the human genome showed a good correlation between the number of hit clones and the number of chromosomes, an indication of unbiased chromosomal distribution of the fosmid library. The chromosomal distribution of the mapped clones is also consistent with the BLAST search result against human and white-cheeked gibbon genomes. The fosmid library and the mapped clones will serve as a valuable resource for further studying gibbons' chromosomal rearrangements and the underlying molecular mechanism as well as for comparative genomic study in the lesser apes.     

Construction, Characterization, and Chromosomal Mapping of a Fosmid Library of the White-Cheeked Gibbon （Nomascus leucogenys）

Gibbons have experienced extensive karyotype rearrangements during evolution and represent an ideal model for studying the underlying molecular mechanism of evolutionary chromosomal rearrangements. It is anticipated that the cloning and sequence characterization of evolutionary chromosomal breakpoints will provide vital insights into the molecular force that has driven such a radical karyotype reshuffle in gibbons. We constructed and characterized a high-quality fosmid li- brary of the white-cheeked gibbon (Nomascus leucogenys) containing 192,000 non- redundant clones with an average insert size of 38 kb and 2.5-fold genome coverage. By end sequencing of 100 randomly selected fosmid clones, we generated 196 se- quence tags for the library. These end-sequenced fosmid clones were then mapped onto the chromosomes of the white-cheeked gibbon by fluorescence in situ hy- bridization, and no spurious chimeric clone was detected. BLAST search against the human genome showed a good correlation between the number of hit clones and the number of chromosomes, an indication of unbiased chromosomal distribu- tion of the fosmid library. The chromosomal distribution of the mapped clones is also consistent with the BLAST search result against human and white-cheeked gibbon genomes. The fosmid library and the mapped clones will serve as a valu- able resource for further studying gibbons' chromosomal rearrangements and the underlying molecular mechanism as well as for comparative genomic study in the lesser apes.     

The transposable elements resident on the plasmids of Pseudomonas putida strain H, Tn 5501 and Tn 5502 , are cryptic transposons of the Tn 3 family

Genes for (methyl)phenol degradation in Pseudomonas putida strain H (phl genes) are located on the plasmid pPGH1. Adjacent to the phl catabolic operon we identified a cryptic transposon, Tn5501, of the Tn3 family (class II transposons). The genes encoding the resolvase and the transposase are transcribed in the same direction, as is common for the Tn501 subfamily. The enzymes encoded by Tn5501, however, show only the overall homology characteristic for resolvases/integrases and transposases of Tn3-type transposons. Therefore it is likely that Tn5501 is not a member of one of the previously defined subfamilies. Inactivation of the conditional lethal sacB gene was used to detect transposition of Tn5501. While screening for transposition events we found another transposon integrated into sacB in one of the sucrose-resistant survivors. This element, Tn5502, is a composite transposon consisting of Tn5501 and an additional DNA fragment. It is flanked by inverted repeats identical to those of Tn5501 and the additional fragment is separated from the Tn5501 portion by an internal repeat (identical to the left terminal repeat). Transposition of phenol degradation genes could not be detected. Analysis of sequence data revealed that the phl genes are not located on a Tn5501-like transposon. Received: 21 July 1997 / Accepted: 7 July 1998     

Two independent mutational events in the loss of urate oxidase during hominoid evolution

Summary Urate oxidase was lost in hominoids during primate evolution. The mechanism and biological reason for this loss remain unknown. In an attempt to address these questions, we analyzed the sequence of urate oxidase genes from four species of hominoids: human (Homo sapiens), chimpanzee (Pan troglodytes), orangutan (Pongo pygmaeus), and gibbon (Hylobates). Two nonsense mutations at codon positions 33 and 187 and an aberrant splice site were found in the human gene. These three deleterious mutations were also identified in the chimpanzee. The nonsense mutation at codon 33 was observed in the orangutan urate oxidase gene. None of the three mutations was present in the gibbon; in contrast, a 13-bp deletion was identified that disrupted the gibbon urate oxidase reading frame. These results suggest that the loss of urate oxidase during the evolution of hominoids could be caused by two independent events after the divergence of the gibbon lineage; the nonsense mutation at codon position 33 resulted in the loss of urate oxidase activity in the human, chimpanzee, and orangutan, whereas the 13-bp deletion was responsible for the urate oxidase deficiency in the gibbon. Because the disruption of a functional gene by independent events in two different evolutionary lineages is unlikely to occur on a chance basis, our data favor the hypothesis that the loss of urate oxidase may have evolutionary advantages. Offprint requests to: C.T. Caskey     

Phylogeny of SINE-R retroposons in Asian apes.

The SINE-R retroposon family was derived from the long terminal repeats (LTRs) of human endogenous retrovirus K (HERV-K) that had been active during the hominoid evolution. The retroposons and HERV-K LTR elements have potential relevance to structural change and genetic variation of the hominoid genome. In our previous study, we found that the SINE-R retroposons were hominoid specific. Here we identified seventeen new SINE-R retroposons (14 from orangutan and 3 from gibbon) from Asian apes and phylogenetically analysed them in comparison with those of the humans and African great apes. None of the retroposons from Asian apes were closely related to SINE-R.C2 that is human specific, and originally identified in the gene for the C2 component of complement, whereas some retroposons (Ch-M10, Ch-M16, Gor-M, Gor-F1, Gor-M6, and Gor-F9) from African great apes showed very close relationship with that of the SINE-R.C2 retroposon. The phylogenetic tree based on the SINE-R retroposons revealed wide overlap of the retroposons across species, suggesting that the SINE-R retroposons have been evolved parallel pattern in the course of hominoid evolution.     

Evolution of an Alu DNA element of type Sx in the lineage of primates and the origin of an associated tetranucleotide microsatellite.

A 394-bp DNA fragment, which in human is on chromosome 6 near the MOG (myelin oligodendrocyte glycoprotein) gene and encompasses an Alu element and an associated tetranucleotide microsatellite, was sequenced from a large range of primate species to follow its evolutionary divergence and to understand the origin of the microsatellite. This Alu element is found at the same orthologous position in all primates sequenced, but the tetranucleotide repeat is present only in Catarrhini between the 3'-oligo(dA) of the Alu element and the 3' flanking direct repeat. Little intraspecific variation was found. Sequence identity values for this orthologous primate Alu averaged 90% (82-99%) with transitions comprising between 70% and 100% of the observed nucleotide substitutions. Although the insertion of the Alu element predates the separation of these species, the original sequence of the site of integration can still be identified. This identification of the direct repeats suggests an active role of the oligo(dA) of the Alu element in the origin of the tetranucleotide repeats. The microsatellite probably appeared after the insertion of the Alu element, early in the lineage leading to the common ancestor of the hominoids and the Old World monkeys.     

Natural genetic variation caused by transposable elements in humans

Transposons and transposon-like repetitive elements collectively occupy 44% of the human genome sequence. In an effort to measure the levels of genetic variation that are caused by human transposons, we have developed a new method to broadly detect transposon insertion polymorphisms of all kinds in humans. We began by identifying 606,093 insertion and deletion (indel) polymorphisms in the genomes of diverse humans. We then screened these polymorphisms to detect indels that were caused by de novo transposon insertions. Our method was highly efficient and led to the identification of 605 nonredundant transposon insertion polymorphisms in 36 diverse humans. We estimate that this represents 25-35% of approximately 2075 common transposon polymorphisms in human populations. Because we identified all transposon insertion polymorphisms with a single method, we could evaluate the relative levels of variation that were caused by each transposon class. The average human in our study was estimated to harbor 1283 Alu insertion polymorphisms, 180 L1 polymorphisms, 56 SVA polymorphisms, and 17 polymorphisms related to other forms of mobilized DNA. Overall, our study provides significant steps toward (i) measuring the genetic variation that is caused by transposon insertions in humans and (ii) identifying the transposon copies that produce this variation.     

The involucrin gene of the gibbon: the middle region shared by the hominoids

The evolution of the anthropoid involucrin gene has resulted largely from a process of vectorial addition of short tandem repeats. The coding region of the involucrin gene of the gibbon (Hylobates lar), including the segment of repeats, has been cloned and sequenced, and its repeat structure can now be compared with that of the other hominoids. In the gibbon, as in the others, repeat additions in the past can be assigned to early, middle, and late regions of the present-day segment of repeats. All 10 repeats of the gibbon early region were completed in a common anthropoid ancestor. All 17 repeats of the gibbon middle region were completed in a common hominoid ancestor. After divergence of the gibbon lineage, eight repeats were added to the middle region of the great ape-human lineages. Seven of these are shared by two to four species, according to the order of their divergences from each other. After its divergence, the gibbon lineage added a short species-specific late region. The gibbon also possesses an incomplete repeat just 3' of the early region, the only addition in this region in any hominoid. Comparison of the number of repeats added with the number of nucleotides substituted shows an inconstant relation between the two.     

Characterization of a Chlamydomonas Transposon, Gulliver, Resembling Those in Higher Plants

While pursuing a chromosomal walk through the mt+ locus of linkage group VI of Chlamydomonas reinhardtii, I encountered a 12-kb sequence that was found to be present in approximately 12 copies in the nuclear genome. Comparison of various C. reinhardtii laboratory strains provided evidence that the sequence was mobile and therefore a transposon. One of two separate natural isolates interfertile with C. reinhardtii, C. smithii (CC-1373), contained the transposon, but at completely different locations in its nuclear genome than C. reinhardtii; and a second, CC-1952 (S1-C5), lacked the transposon altogether. Genetic analysis indicated that the transposon was found at dispersed sites throughout the genome, but had a conserved structure at each location. Sequence homology between the termini was limited to an imperfect 15-bp inverted repeat. An 8-bp target site duplication was created by insertion; transposon sequences were completely removed upon excision leaving behind both copies of the target site duplication, with minor base changes. The transposon contained an internal region of unique repetitive sequence responsible for restriction fragment length heterogeneity among the various copies of the transposon. In several cases it was possible to identify which of the dozen transposons in a given strain served as the donor when a transposition event occurred. The transposon often moved into a site genetically linked to the donor, and transposition appeared to be nonreplicative. Thus the mechanism of transposition and excision of the transposon, which I have named Gulliver, resembles that of certain higher plant transposons, like the Ac transposon of maize.     

Structure analysis of a class II transposon encoding the mercury resistance of the Gram-positive Bacterium bacillus megaterium MB1, a strain isolated from minamata bay, Japan.

A unique transposon was found in the chromosome of Bacillus megaterium MB1, a Gram-positive bacterium isolated from mercury-polluted sediments of Minamata Bay, Japan. The transposon region of a 14.5kb DNA fragment was amplified by PCR using a single PCR primer designed from the nucleotide sequence of an inverted repeat of class II transposons. The molecular analysis revealed that the PCR-amplified DNA fragment encodes a transposition module similar to that of Tn21. The transposon also encodes a broad-spectrum mercury resistance region having a restriction endonuclease map identical to that of Bacillus cereus RC607, a strain isolated from Boston Harbor, USA. The result of a phylogenetic analysis of the amino acid sequence of putative resolvase of the transposon showed that the transposon is phylogenetically closer to the transposons of Gram-positive bacteria than those of Gram-negative bacteria. Besides the transposition module and mer operon, the transposon encodes a mobile genetic element of bacterial group II introns between the resolvase gene and mer operon. The intron, however, does not intervene in any exon gene. The discovery of this newly found combination of the complex mobile elements may offer a clue to understanding the horizontal dissemination of broad-spectrum mercury resistance among microbes.     

The Chinese hamster Alu-equivalent sequence: a conserved highly repetitious, interspersed deoxyribonucleic acid sequence in mammals has a structure suggestive of a transposable element.

A consensus sequence has been determined for a major interspersed deoxyribonucleic acid repeat in the genome of Chinese hamster ovary cells (CHO cells). This sequence is extensively homologous to (i) the human Alu sequence (P. L. Deininger et al., J. Mol. Biol., in press), (ii) the mouse B1 interspersed repetitious sequence (Krayev et al., Nucleic Acids Res. 8:1201-1215, 1980) (iii) an interspersed repetitious sequence from African green monkey deoxyribonucleic acid (Dhruva et al., Proc. Natl. Acad. Sci. U.S.A. 77:4514-4518, 1980) and (iv) the CHO and mouse 4.5S ribonucleic acid (this report; F. Harada and N. Kato, Nucleic Acids Res. 8:1273-1285, 1980). Because the CHO consensus sequence shows significant homology to the human Alu sequence it is termed the CHO Alu-equivalent sequence. A conserved structure surrounding CHO Alu-equivalent family members can be recognized. It is similar to that surrounding the human Alu and the mouse B1 sequences, and is represented as follows: direct repeat-CHO-Alu-A-rich sequence-direct repeat. A composite interspersed repetitious sequence has been identified. Its structure is represented as follows: direct repeat-residue 47 to 107 of CHO-Alu-non-Alu repetitious sequence-A-rich sequence-direct repeat. Because the Alu flanking sequences resemble those that flank known transposable elements, we think it likely that the Alu sequence dispersed throughout the mammalian genome by transposition.