Modeling the amplification dynamics of human Alu retrotransposons

Alu insertion polymorphisms shared by Papio baboons and Theropithecus gelada reveal an intertwined common ancestry

Background

Baboons (genus Papio) and geladas (Theropithecus gelada) are now generally recognized as close phylogenetic relatives, though morphologically quite distinct and generally classified in separate genera. Primate specific Alu retrotransposons are well-established genomic markers for the study of phylogenetic and population genetic relationships. We previously reported a computational reconstruction of Papio phylogeny using large-scale whole genome sequence (WGS) analysis of Alu insertion polymorphisms. Recently, high coverage WGS was generated for Theropithecus gelada. The objective of this study was to apply the high-throughput "poly-Detect" method to computationally determine the number of Alu insertion polymorphisms shared by T. gelada and Papio, and vice versa, by each individual Papio species and T. gelada. Secondly, we performed locus-specific polymerase chain reaction (PCR) assays on a diverse DNA panel to complement the computational data.

Results

We identified 27,700 Alu insertions from T. gelada WGS that were also present among six Papio species, with nearly half (12,956) remaining unfixed among 12 Papio individuals. Similarly, each of the six Papio species had species-indicative Alu insertions that were also present in T. gelada. In general, P. kindae shared more insertion polymorphisms with T. gelada than did any of the other five Papio species. PCR-based genotype data provided additional support for the computational findings.

Conclusions

Our discovery that several thousand Alu insertion polymorphisms are shared by T. gelada and Papio baboons suggests a much more permeable reproductive barrier between the two genera then previously suspected. Their intertwined evolution likely involves a long history of admixture, gene flow and incomplete lineage sorting.

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

Background

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

Results

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

Conclusion

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

Estimating the age of retrotransposon subfamilies using maximum likelihood

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

Mobile DNA elements in primate and human evolution

Roughly 50% of the primate genome consists of mobile, repetitive DNA sequences such as Alu and LINE1 elements. The causes and evolutionary consequences of mobile element insertion, which have received considerable attention during the past decade, are reviewed in this article. Because of their unique mutational mechanisms, these elements are highly useful for answering phylogenetic questions. We demonstrate how they have been used to help resolve a number of questions in primate phylogeny, including the human-chimpanzee-gorilla trichotomy and New World primate phylogeny. Alu and LINE1 element insertion polymorphisms have also been analyzed in human populations to test hypotheses about human evolution and population affinities and to address forensic issues. Finally, these elements have had impacts on the genome itself. We review how they have influenced fundamental ongoing processes like nonhomologous recombination, genomic deletion, and X chromosome inactivation.

The heterochromatic copies of the LTR retrotransposons as a record of the genomic events that have shaped the Drosophila melanogaster genome

Transposable elements, which are major components of most genomes, are known to accumulate in heterochromatic regions in which they have progressively diverged in sequence by mutations and internal deletions and insertions (indels) during the course of evolution. They therefore provide a record of the genomic events that have shaped the genomes, some of which could correspond to speciation events. Using the sequence divergence between the long terminal repeats (LTRs), we estimated the date of the insertion events of the LTR retrotransposon copies embedded within the heterochromatin regions of the Drosophila melanogaster genome. We did not detect traces of any specific waves of mobilization of retrotransposons within heterochromatin, apart from a very recent wave, which corresponds to the numerous LTR retrotransposon copies found in euchromatin.

Analysis of the features and source gene composition of the AluYg6 subfamily of human retrotransposons

Background

Alu elements are a family of SINE retrotransposons in primates. They are classified into subfamilies according to specific diagnostic mutations from the general Alu consensus. It is now believed that there may be several retrotranspositionally-competent source genes within an Alu subfamily. To investigate the evolution of young Alu elements it is critical to have access to complete subfamilies, which, following the release of the final human genome assembly, can now be obtained using in silico methods.

Results

380 elements belonging to the young AluYg6 subfamily were identified in the human genome, a number significantly exceeding prior expectations. An AluYg6 element was also identified in the chimpanzee genome, indicating that the subfamily is older than previously estimated, and appears to have undergone a period of dormancy before its expansion. The relative contributions of back mutation and gene conversion to variation at the six diagnostic positions are examined, and cases of complete forward gene conversion events are reported. Two small subfamilies derived from AluYg6 have been identified, named AluYg6a2 and AluYg5b3, which contain 40 and 27 members, respectively. These small subfamilies are used to illustrate the ambiguity regarding Alu subfamily definition, and to assess the contribution of secondary source genes to the AluYg6 subfamily.

Conclusion

The number of elements in the AluYg6 subfamily greatly exceeds prior expectations, indicating that the current knowledge of young Alu subfamilies is incomplete, and that prior analyses that have been carried out using these data may have generated inaccurate results. A definition of primary and secondary source genes has been provided, and it has been shown that several source genes have contributed to the proliferation of the AluYg6 subfamily. Access to the sequence data for the complete AluYg6 subfamily will be invaluable in future computational analyses investigating the evolution of young Alu subfamilies.

piggyBac transposon mediated transgenesis of the human blood fluke, Schistosoma mansoni

The transposon piggyBac from the genome of the cabbage looper moth Trichoplusia ni has been observed in the laboratory to jump into the genomes of key model and pathogenic eukaryote organisms including mosquitoes, planarians, human and other mammalian cells, and the malaria parasite Plasmodium falciparum. Introduction of exogenous transposons into schistosomes has not been reported but transposon-mediated transgenesis of schistosomes might supersede current methods for functional genomics of this important human pathogen. In the present study we examined whether the piggyBac transposon could deliver reporter transgenes into the genome of Schistosoma mansoni parasites. A piggyBac donor plasmid modified to encode firefly luciferase under control of schistosome gene promoters was introduced along with 7-methylguanosine capped RNAs encoding piggyBac transposase into cultured schistosomula by square wave electroporation. The activity of the helper transposase mRNA was confirmed by Southern hybridization analysis of genomic DNA from the transformed schistosomes, and hybridization signals indicated that the piggyBac transposon had integrated into numerous sites within the parasite chromosomes. piggyBac integrations were recovered by retrotransposon-anchored PCR, revealing characteristic piggyBac TTAA footprints in the vicinity of the endogenous schistosome retrotransposons Boudicca, SR1, and SR2. This is the first report of chromosomal integration of a transgene and somatic transgenesis of this important human pathogen, in this instance accomplished by mobilization of the piggyBac transposon.

Human population genetic structure and diversity inferred from polymorphic L1(LINE-1) and Alu insertions

BACKGROUND/AIMS

The L1 retrotransposable element family is the most successful self-replicating genomic parasite of the human genome. L1 elements drive replication of Alu elements, and both have had far-reaching impacts on the human genome. We use L1 and Alu insertion polymorphisms to analyze human population structure.

METHODS

We genotyped 75 recent, polymorphic L1 insertions in 317 individuals from 21 populations in sub-Saharan Africa, East Asia, Europe and the Indian subcontinent. This is the first sample of L1 loci large enough to support detailed population genetic inference. We analyzed these data in parallel with a set of 100 polymorphic Alu insertion loci previously genotyped in the same individuals.

RESULTS AND CONCLUSION

The data sets yield congruent results that support the recent African origin model of human ancestry. A genetic clustering algorithm detects clusters of individuals corresponding to continental regions. The number of loci sampled is critical: with fewer than 50 typical loci, structure cannot be reliably discerned in these populations. The inclusion of geographically intermediate populations (from India) reduces the distinctness of clustering. Our results indicate that human genetic variation is neither perfectly correlated with geographic distance (purely clinal) nor independent of distance (purely clustered), but a combination of both: stepped clinal.

Different evolutionary fates of recently integrated human and chimpanzee LINE-1 retrotransposons

The long interspersed element-1 (LINE-1 or L1) is a highly successful retrotransposon in mammals. L1 elements have continued to actively propagate subsequent to the human-chimpanzee divergence, approximately 6 million years ago, resulting in species-specific inserts. Here, we report a detailed characterization of chimpanzee-specific L1 subfamily diversity and a comparison with their human-specific counterparts. Our results indicate that L1 elements have experienced different evolutionary fates in humans and chimpanzees within the past approximately 6 million years. Although the species-specific L1 copy numbers are on the same order in both species (1200-2000 copies), the number of retrotransposition-competent elements appears to be much higher in the human genome than in the chimpanzee genome. Also, while human L1 subfamilies belong to the same lineage, we identified two lineages of recently integrated L1 subfamilies in the chimpanzee genome. The two lineages seem to have coexisted for several million years, but only one shows evidence of expansion within the past three million years. These differential evolutionary paths may be the result of random variation, or the product of competition between L1 subfamily lineages. Our results suggest that the coexistence of several L1 subfamily lineages within a species may be resolved in a very short evolutionary period of time, perhaps in just a few million years. Therefore, the chimpanzee genome constitutes an excellent model in which to analyze the evolutionary dynamics of L1 retrotransposons.

Estimating the retrotransposition rate of human Alu elements

Mobile elements such as Alu repeats have substantially altered the architecture of the human genome, and de novo mobile element insertions sometimes cause genetic disorders. Previous estimates for the retrotransposition rate (RR) of Alu elements in humans of one new insertion every approximately 100-125 births were developed prior to the sequencing of the human and chimpanzee genomes. Here, we used two independent methods (based on the new genomic data and on disease-causing de novo Alu insertions) to generate refined Alu RR estimates in humans. Both methods consistently yielded RR on the order of one new Alu insertion every approximately 20 births, despite the fact that the evolutionary-based method represents an average RR over the past approximately 6 million years while the mutation-based method better reflects the current-day RR. These results suggest that Alu elements retrotranspose at a faster rate in humans than previously thought, and support the potential of Alu elements as mutagenic factors in the human genome.