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.     

Overlapping genes as rare genomic markers: the phylogeny of gamma-Proteobacteria as a case study

Phylogenies can be constructed in many ways, including using shared complex characters known as rare genomic changes (RGCs), such as insertions and deletions (indels), retroposon integrations and intron positions. Here, we demonstrate that distance-based phylogenies, which were determined by shared overlapping genes from 13 completely sequenced gamma-Proteobacteria genomes, are consistent with phylogenies based on 16S rRNAs and other robust markers. These findings suggest that overlapping genes could provide interesting additional insights into the phylogenomics of completely sequenced microbial genomes.     

Discovery of a large number of previously unrecognized mitochondrial pseudogenes in fish genomes

Nuclear inserted copies of mitochondrial origin (numts) vary widely among eukaryotes, with human and plant genomes harboring the largest repertoires. Numts were previously thought to be absent from fish species, but the recent release of three fish nuclear genome sequences provides the resource to obtain a more comprehensive insight into the extent of mtDNA transfer in fishes. From the sequence analyses of the genomes of Fugu rubripes, Tetraodon nigroviridis, and Danio rerio, we have identified 2, 5, and 10 recent numt integrations, respectively, which integrated into those genomes less than 0.6 million years (Myr) ago. Such results contradict the hypothesis of absence or rarity of numts in fishes, as (i) the ratio of numts to the total size of the nuclear genome in T. nigroviridis was superior to the ratio observed in several higher vertebrate species (e.g., chicken, mouse, and rat), and only surpassed by humans, and (ii) the mtDNA coverage transferred to the nuclear genome of D. rerio is exceeded only by human and mouse, within the whole range of eukaryotic genomes surveyed for numts. Additionally, 335, 336, and 471 old numts (>12.5 Myr) were detected in F. rubripes, T. nigroviridis, and D. rerio, respectively. Surprisingly, old numts are inserted preferentially into known or predicted genes, as inferred for recent numts in human. However, because in fish genomes such integrations are old, they are likely to represent evolutionary successes and they may be considered a potential important evolutionary mechanism for the enhancement of genomic coding regions.     

Divergent patterns of recent retroviral integrations in the human and chimpanzee genomes: probable transmissions between other primates and chimpanzees

The human genome is littered by endogenous retrovirus sequences (HERVs), which constitute up to 8% of the total genomic sequence. The sequencing of the human (Homo sapiens) and chimpanzee (Pan troglodytes) genomes has facilitated the evolutionary study of ERVs and related sequences. We screened both the human genome (version hg16) and the chimpanzee genome (version PanTro1) for ERVs and conducted a phylogenetic analysis of recent integrations. We found a number of recent integrations within both genomes. They segregated into four groups. Two larger gammaretrovirus-like groups (PtG1 and PtG2) occurred in chimpanzees but not in humans. The PtG sequences were most similar to two baboon ERVs and a macaque sequence but neither to other chimpanzee ERVs nor to any human gammaretrovirus-like ERVs. The pattern was consistent with cross-species transfer via predation. This appears to be an example of horizontal transfer of retroviruses with occasional fixation in the germ line.     

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.     

Unexpected Inheritance: Multiple Integrations of Ancient Bornavirus and Ebolavirus/Marburgvirus Sequences in Vertebrate Genomes

Vertebrate genomes contain numerous copies of retroviral sequences, acquired over the course of evolution. Until recently they were thought to be the only type of RNA viruses to be so represented, because integration of a DNA copy of their genome is required for their replication. In this study, an extensive sequence comparison was conducted in which 5,666 viral genes from all known non-retroviral families with single-stranded RNA genomes were matched against the germline genomes of 48 vertebrate species, to determine if such viruses could also contribute to the vertebrate genetic heritage. In 19 of the tested vertebrate species, we discovered as many as 80 high-confidence examples of genomic DNA sequences that appear to be derived, as long ago as 40 million years, from ancestral members of 4 currently circulating virus families with single strand RNA genomes. Surprisingly, almost all of the sequences are related to only two families in the Order Mononegavirales: the Bornaviruses and the Filoviruses, which cause lethal neurological disease and hemorrhagic fevers, respectively. Based on signature landmarks some, and perhaps all, of the endogenous virus-like DNA sequences appear to be LINE element-facilitated integrations derived from viral mRNAs. The integrations represent genes that encode viral nucleocapsid, RNA-dependent-RNA-polymerase, matrix and, possibly, glycoproteins. Integrations are generally limited to one or very few copies of a related viral gene per species, suggesting that once the initial germline integration was obtained (or selected), later integrations failed or provided little advantage to the host. The conservation of relatively long open reading frames for several of the endogenous sequences, the virus-like protein regions represented, and a potential correlation between their presence and a species'' resistance to the diseases caused by these pathogens, are consistent with the notion that their products provide some important biological advantage to the species. In addition, the viruses could also benefit, as some resistant species (e.g. bats) may serve as natural reservoirs for their persistence and transmission. Given the stringent limitations imposed in this informatics search, the examples described here should be considered a low estimate of the number of such integration events that have persisted over evolutionary time scales. Clearly, the sources of genetic information in vertebrate genomes are much more diverse than previously suspected.     

Unique profile of ordered arrangements of repetitive elements in the C57BL/6J mouse genome implicating their functional roles

The entirety of all protein coding sequences is reported to represent a small fraction (~2%) of the mouse and human genomes; the vast majority of the rest of the genome is presumed to be repetitive elements (REs). In this study, the C57BL/6J mouse reference genome was subjected to an unbiased RE mining to establish a whole-genome profile of RE occurrence and arrangement. The C57BL/6J mouse genome was fragmented into an initial set of 5,321 units of 0.5 Mb, and surveyed for REs using unbiased self-alignment and dot-matrix protocols. The survey revealed that individual chromosomes had unique profiles of RE arrangement structures, named RE arrays. The RE populations in certain genomic regions were arranged into various forms of complexly organized structures using combinations of direct and/or inverse repeats. Some of these RE arrays spanned stretches of over 2 Mb, which may contribute to the structural configuration of the respective genomic regions. There were substantial differences in RE density among the 21 chromosomes, with chromosome Y being the most densely populated. In addition, the RE array population in the mouse chromosomes X and Y was substantially different from those of the reference human chromosomes. Conversion of the dot-matrix data pertaining to a tandem 13-repeat structure within the Ch7.032 genome unit into a line map of known REs revealed a repeat unit of ~11.3 Kb as a mosaic of six different RE types. The data obtained from this study allowed for a comprehensive RE profiling, including the establishment of a library of RE arrays, of the reference mouse genome. Some of these RE arrays may participate in a spectrum of normal and disease biology that are specific for mice.     

Retroviruses and primate evolution

Human endogenous retroviruses (HERVs), probably representing footprints of ancient germ-cell retroviral infections, occupy about 1% of the human genome. HERVs can influence genome regulation through expression of retroviral genes, either via genomic rearrangements following HERV integrations or through the involvement of HERV LTRs in the regulation of gene expression. Some HERVs emerged in the genome over 30 MYr ago, while others have appeared rather recently, at about the time of hominid and ape lineages divergence. HERVs might have conferred antiviral resistance on early human ancestors, thus helping them to survive. Furthermore, newly integrated HERVs could have changed the pattern of gene expression and therefore played a significant role in the evolution and divergence of Hominoidea superfamily. Comparative analysis of HERVs, HERV LTRs, neighboring genes, and their regulatory interplay in the human and ape genomes will help us to understand the possible impact of HERVs on evolution and genome regulation in the primates. BioEssays 22:161-171, 2000.     

A rare event of insertion polymorphism of a HERV-K LTR in the human genome

Human endogenous retroviruses (HERVs), which constitute a significant part of the human genome, might have a serious impact on primate evolution. Over a hundred insertions of HERV-K(HML-2) family members distinguish the human genome from other primate genomes. However, only three cases of insertion polymorphisms have been reported so far, all for endogenous HERV-K proviruses. This suggests that some retroviral integrations occurred rather recently in human genome evolution. In this report, we describe a very rare case of true insertion polymorphism of a solitary HERV-K LTR in the human genome. Distribution of the LTR-containing allele was tested in 5 Africans and 83 individuals from three Russian populations. The allele frequency appeared to be relatively high in populations of both European and Asian origin. The detected polymorphic LTR could be a useful molecular genetic marker of the corresponding genomic region.     

REMiner-II: A tool for rapid identification and configuration of repetitive element arrays from large mammalian chromosomes as a single query

Genes occupy ~3% of the human and mouse genomes whereas repetitive elements (REs), whose biologic functions are largely uncharacterized, constitute greater than 50%. A heterogeneous population of RE arrays (arrangement structures) is formed by combinations of various REs in mammalian genomes. In this study, REMiner-II was refined from the original REMiner for a more efficient identification and configuration of RE arrays from large queries (e.g., human chromosomes) using an unbiased self-alignment protocol. Chromosome-wide RE array profiles for the entire sets of human and mouse chromosomes were obtained using REMiner-II on a personal computer. REMiner-II provides 10 adjustable parameters and three data output modes to accommodate different experimental settings and/or goals. Examination of the human and mouse chromosome data using the REMiner-II viewer revealed species-specific libraries of complexly organized RE arrays. In conclusion, REMiner-II is an efficient tool for chromosome-wide identification and characterization of RE arrays from mammalian genomes.     

The MitoDrome database annotates and compares the OXPHOS nuclear genes of Drosophila melanogaster, Drosophila pseudoobscura and Anopheles gambiae

The oxidative phosphorylation (OXPHOS) is the primary energy-producing process of all aerobic organisms and the only cellular function under the dual control of both the mitochondrial and the nuclear genomes. Functional characterization and evolutionary study of the OXPHOS system is of great importance for the understanding of many as yet unclear aspects of nucleus-mitochondrion genomic co-evolution and co-regulation gene networks. The MitoDrome database is a web-based database which provides genomic annotations about nuclear genes of Drosophila melanogaster encoding for mitochondrial proteins. Recently, MitoDrome has included a new section annotating genomic information about OXPHOS genes in Drosophila pseudoobscura and Anopheles gambiae and their comparative analysis with their Drosophila melanogaster and human counterparts. The introduction of this new comparative annotation section into MitoDrome is expected to be a useful resource for both functional and structural genomics related to the OXPHOS system.     

The Drosophila gene disruption project: progress using transposons with distinctive site specificities

The Drosophila Gene Disruption Project (GDP) has created a public collection of mutant strains containing single transposon insertions associated with different genes. These strains often disrupt gene function directly, allow production of new alleles, and have many other applications for analyzing gene function. Here we describe the addition of ～7600 new strains, which were selected from >140,000 additional P or piggyBac element integrations and 12,500 newly generated insertions of the Minos transposon. These additions nearly double the size of the collection and increase the number of tagged genes to at least 9440, approximately two-thirds of all annotated protein-coding genes. We also compare the site specificity of the three major transposons used in the project. All three elements insert only rarely within many Polycomb-regulated regions, a property that may contribute to the origin of "transposon-free regions" (TFRs) in metazoan genomes. Within other genomic regions, Minos transposes essentially at random, whereas P or piggyBac elements display distinctive hotspots and coldspots. P elements, as previously shown, have a strong preference for promoters. In contrast, piggyBac site selectivity suggests that it has evolved to reduce deleterious and increase adaptive changes in host gene expression. The propensity of Minos to integrate broadly makes possible a hybrid finishing strategy for the project that will bring >95% of Drosophila genes under experimental control within their native genomic contexts.     

Conservation of the MHC-like region throughout evolution

Identification of conserved regions between the genomes of distant species is a crucial step in the reconstruction of the genomic organization of their last common ancestor. Here we confirm for the first time with robust evidence, the existence of a region of conserved synteny between the human genome and the Drosophila genome. This evolutionarily conserved synteny involves the human MHC and paralogous regions, and we identified 19 conserved genes between these two species in a Drosophila genomic region of less than 2 Mb. The statistical analysis of the distribution of these 19 genes between the Drosophila and human genomes shows that it cannot be explained by chance. Our study constitutes a first step towards the reconstruction of the genome of Urbilateria (the ancestor of all bilaterian) and allows for a better understanding of the evolutionary history of our genome as well as other metazoan genomes.     

Genome-wide comparison of differences in the integration sites of interspersed repeats between closely related genomes

A technique for genome-wide detection of differences in the integration site positions of interspersed repeats in related genomes (DiffIR) is described. The technique is based on a whole- genome selective PCR amplification of the repeats’ flanking regions followed by a differential hybridization screening of the arrayed library of the selected amplicons. The technique was successfully applied to the comparison of the integration sites in the human and chimpanzee genomes, allowing us to discover 11 new human-specific integrations of human endogenous retrovirus, K family (HML-2) long terminal repeats.     

Alu recombination-mediated structural deletions in the chimpanzee genome

With more than 1.2 million copies, Alu elements are one of the most important sources of structural variation in primate genomes. Here, we compare the chimpanzee and human genomes to determine the extent of Alu recombination-mediated deletion (ARMD) in the chimpanzee genome since the divergence of the chimpanzee and human lineages (~6 million y ago). Combining computational data analysis and experimental verification, we have identified 663 chimpanzee lineage-specific deletions (involving a total of ~771 kb of genomic sequence) attributable to this process. The ARMD events essentially counteract the genomic expansion caused by chimpanzee-specific Alu inserts. The RefSeq databases indicate that 13 exons in six genes, annotated as either demonstrably or putatively functional in the human genome, and 299 intronic regions have been deleted through ARMDs in the chimpanzee lineage. Therefore, our data suggest that this process may contribute to the genomic and phenotypic diversity between chimpanzees and humans. In addition, we found four independent ARMD events at orthologous loci in the gorilla or orangutan genomes. This suggests that human orthologs of loci at which ARMD events have already occurred in other nonhuman primate genomes may be “at-risk” motifs for future deletions, which may subsequently contribute to human lineage-specific genetic rearrangements and disorders.