Human Retrotransposon Insertion Polymorphisms Are Associated with Health and Disease via Gene Regulatory Phenotypes.

The human genome hosts several active families of transposable elements (TEs), including the Alu, LINE-1, and SVA retrotransposons that are mobilized via reverse transcription of RNA intermediates. We evaluated how insertion polymorphisms generated by human retrotransposon activity may be related to common health and disease phenotypes that have been previously interrogated through genome-wide association studies (GWAS). To address this question, we performed a genome-wide screen for retrotransposon polymorphism disease associations that are linked to TE induced gene regulatory changes. Our screen first identified polymorphic retrotransposon insertions found in linkage disequilibrium (LD) with single nucleotide polymorphisms that were previously associated with common complex diseases by GWAS. We further narrowed this set of candidate disease associated retrotransposon polymorphisms by identifying insertions that are located within tissue-specific enhancer elements. We then performed expression quantitative trait loci analysis on the remaining set of candidates in order to identify polymorphic retrotransposon insertions that are associated with gene expression changes in B-cells of the human immune system. This progressive and stringent screen yielded a list of six retrotransposon insertions as the strongest candidates for TE polymorphisms that lead to disease via enhancer-mediated changes in gene regulation. For example, we found an SVA insertion within a cell-type specific enhancer located in the second intron of the B4GALT1 gene. B4GALT1 encodes a glycosyltransferase that functions in the glycosylation of the Immunoglobulin G (IgG) antibody in such a way as to convert its activity from pro- to anti-inflammatory. The disruption of the B4GALT1 enhancer by the SVA insertion is associated with down-regulation of the gene in B-cells, which would serve to keep the IgG molecule in a pro-inflammatory state. Consistent with this idea, the B4GALT1 enhancer SVA insertion is linked to a genomic region implicated by GWAS in both inflammatory conditions and autoimmune diseases, such as systemic lupus erythematosus and Crohn's disease. We explore this example and the other cases uncovered by our genome-wide screen in an effort to illuminate how retrotransposon insertion polymorphisms can impact human health and disease by causing changes in gene expression.

Exaptation of protein coding sequences from transposable elements.

The activity of transposable elements (TEs) has had a profound impact on the evolution of eukaryotic genomes. Once thought to be purely selfish genomic entities, TEs are now recognized to occupy a continuum of relationships, ranging from parasitic to mutualistic, with their host genomes. One of the many ways that TEs contribute to the function and evolution of the genomes in which they reside is through the donation of host protein coding sequences (CDSs). In this chapter, we will describe several notable examples of eukaryotic host CDSs that are derived from TEs. Despite the existence of a number of such well-established cases, the overall extent and significance of this phenomenon remains a matter of controversy. Genome-scale computational analyses have yielded vastly different estimates for the fraction of host CDSs that are derived from TEs. We explain how these seemingly contradictory findings are the result of specific ascertainment biases introduced by the different methods used to detect TE-related sequences. In light of this problem, we propose a comprehensive and systematic framework for definitively characterizing the contribution of TEs to eukaryotic CDSs.

Transposable elements donate lineage-specific regulatory sequences to host genomes.

The evolutionary implications of transposable element (TE) influences on gene regulation are explored here. An historical perspective is presented to underscore the importance of TE influences on gene regulation with respect to both the discovery of TEs and the early conceptualization of their potential impact on host genome evolution. Evidence that points to a role for TEs in host gene regulation is reviewed, and comparisons between genome sequences are used to demonstrate the fact that TEs are particularly lineage-specific components of their host genomes. Consistent with these two properties of TEs, regulatory effects and evolutionary specificity, human-mouse genome wide sequence comparisons reveal that the regulatory sequences that are contributed by TEs are exceptionally lineage specific. This suggests a particular mechanism by which TEs may drive the diversification of gene regulation between evolutionary lineages.

Independent evolution of heavy metal-associated domains in copper chaperones and copper-transporting atpases.

Copper chaperones are small cytoplasmic proteins that bind intracellular copper (Cu) and deliver it to Cu-dependent enzymes such as cytochrome oxidase, superoxide dismutase, and amine oxidase. Copper chaperones are similar in sequence and structure to the Cu-binding heavy metal-associated (HMA) domains of Cu-transporting ATPases (Cu-ATPases), and the genes for copper chaperones and Cu-ATPases are often located in the same operon. Phylogenetic analysis shows that Cu chaperones and HMA domains of Cu-ATPases represent ancient and distinct lineages that have evolved largely independently since their initial separation. Copper chaperone-Cu-ATPase operons appear to have evolved independently in different prokaryotic lineages, probably due to a strong selective pressure for coexpression of these genes.

Sequence and structural aspects of functional diversification in class I alpha-mannosidase evolution.

MOTIVATION: Class I alpha-mannosidases comprise a homologous and functionally diverse family of glycoside hydrolases. Phylogenetic analysis based on an amino acid sequence alignment of the catalytic domain of class I alpha-mannosidases reveals four well-supported phylogenetic groups within this family. These groups include a number of paralogous members generated by gene duplications that occurred as far back as the initial divergence of the crown-group of eukaryotes. Three of the four phylogenetic groups consist of enzymes that have group-specific biochemical specificity and/or sites of activity. An attempt has been made to uncover the role that natural selection played in the sequence and structural divergence between the phylogenetically and functionally distinct Endoplasmic Reticulum (ER) and Golgi apparatus groups. RESULTS: Comparison of site-specific amino acid variability profiles for the ER and Golgi groups revealed statistically significant evidence for functional diversification at the sequence level and indicated a number of residues that are most likely to have played a role in the functional divergence between the two groups. The majority of these sites appear to contain residues that have been fixed within one organelle-specific group by positive selection. Somewhat surprisingly these selected residues map to the periphery of the alpha-mannosidase catalytic domain tertiary structure. Changes in these peripherally located residues would not seem to have a gross effect on protein function. Thus diversifying selection between the two groups may have acted in a gradual manner consistent with the Darwinian model of natural selection. CONTACT: bishogr@millsaps.edu.

Lineage-specific gene expansions in bacterial and archaeal genomes.

Gene duplication is an important mechanistic antecedent to the evolution of new genes and novel biochemical functions. In an attempt to assess the contribution of gene duplication to genome evolution in archaea and bacteria, clusters of related genes that appear to have expanded subsequent to the diversification of the major prokaryotic lineages (lineage-specific expansions) were analyzed. Analysis of 21 completely sequenced prokaryotic genomes shows that lineage-specific expansions comprise a substantial fraction (approximately 5%-33%) of their coding capacities. A positive correlation exists between the fraction of the genes taken up by lineage-specific expansions and the total number of genes in a genome. Consistent with the notion that lineage-specific expansions are made up of relatively recently duplicated genes, >90% of the detected clusters consists of only two to four genes. The more common smaller clusters tend to include genes with higher pairwise similarity (as reflected by average score density) than larger clusters. Regardless of size, cluster members tend to be located more closely on bacterial chromosomes than expected by chance, which could reflect a history of tandem gene duplication. In addition to the small clusters, almost all genomes also contain rare large clusters of size > or =20. Several examples of the potential adaptive significance of these large clusters are explored. The presence or absence of clusters and their related genes was used as the basis for the construction of a similarity graph for completely sequenced prokaryotic genomes. The topology of the resulting graph seems to reflect a combined effect of common ancestry, horizontal transfer, and lineage-specific gene loss.

Gene conversions in genes encoding outer-membrane proteins in H. pylori and C. pneumoniae.

Helicobacter pylori and Chlamydia pneumoniae are both pathogenic to humans. Their genomes have recently been completed, allowing detailed study of their evolution and organization. Here we describe an evolutionary analysis of the H. pylori and C. pneumoniae genes that encode their outer-membrane proteins. By comparing complete genome sequences of two H. pylori strains and two C. pneumoniae strains, we identify multiple independent conversions among these genes. Such recombination events might provide a selective advantage for these bacterial pathogens.

Constant relative rate of protein evolution and detection of functional diversification among bacterial, archaeal and eukaryotic proteins.

BACKGROUND: Detection of changes in a protein's evolutionary rate may reveal cases of change in that protein's function. We developed and implemented a simple relative rates test in an attempt to assess the rate constancy of protein evolution and to detect cases of functional diversification between orthologous proteins. The test was performed on clusters of orthologous protein sequences from complete bacterial genomes (Chlamydia trachomatis, C. muridarum and Chlamydophila pneumoniae), complete archaeal genomes (Pyrococcus horikoshii, P. abyssi and P. furiosus) and partially sequenced mammalian genomes (human, mouse and rat). RESULTS: Amino-acid sequence evolution rates are significantly correlated on different branches of phylogenetic trees representing the great majority of analyzed orthologous protein sets from all three domains of life. However, approximately 1% of the proteins from each group of species deviates from this pattern and instead shows variation that is consistent with an acceleration of the rate of amino-acid substitution, which may be due to functional diversification. Most of the putative functionally diversified proteins from all three species groups are predicted to function at the periphery of the cells and mediate their interaction with the environment. CONCLUSIONS: Relative rates of protein evolution are remarkably constant for the three species groups analyzed here. Deviations from this rate constancy are probably due to changes in selective constraints associated with diversification between orthologs. Functional diversification between orthologs is thought to be a relatively rare event. However, the resolution afforded by the test designed specifically for genomic-scale datasets allowed us to identify numerous cases of possible functional diversification between orthologous proteins.

The alpha-mannosidases: phylogeny and adaptive diversification.

alpha-Mannosidase enzymes comprise a class of gylcoside hydrolases involved in the maturation and degradaton of glycoprotein-linked oligosaccharides. Various alpha-mannosidase enzymatic activities are encoded by an ancient and ubiquitous gene superfamily. A comparative sequence analysis was employed here to characterize the evolutionary relationships and dynamics of the alpha-mannosidase superfamily. A series of lineage-specific BLAST searches recovered the first ever recognized archaean and eubacterial alpha-mannosidase sequences, in addition to numerous eukaryotic sequences. Motif-based alignment and subsequent phylogenetic analysis of the entire superfamily revealed the presence of three well-supported monophyletic clades that represent discrete alpha-mannosidase families. The comparative method was used to evaluate the phylogenetic distribution of alpha-mannosidase functional variants within families. Results of this analysis demonstrate a pattern of functional diversification of alpha-mannosidase paralogs followed by conservation of function among orthologs. Nucleotide polymorphism among the most closely related pair of duplicated genes was analyzed to evaluate the role of natural selection in the functional diversification of alpha-mannosidase paralogs. Ratios of nonsynonymous and synonymous variation show an increase in the rate of nonsynonymous change after duplication and a relative excess of fixed nonsynonymous changes between the two groups of paralogs. These data point to a possible role for positive Darwinian selection in the evolution of alpha-mannosidase functional diversification following gene duplication.

Molecular evolution of the Paramyxoviridae and Rhabdoviridae multiple-protein-encoding P gene.

Presented here is an analysis of the molecular evolutionary dynamics of the P gene among 76 representative sequences of the Paramyxoviridae and Rhabdoviridae RNA virus families. In a number of Paramyxoviridae taxa, as well as in vesicular stomatitis viruses of the Rhabdoviridae, the P gene encodes multiple proteins from a single genomic RNA sequence. These products include the phosphoprotein (P), as well as the C and V proteins. The complexity of the P gene makes it an intriguing locus to study from an evolutionary perspective. Amino acid sequence alignments of the proteins encoded at the P and N loci were used in independent phylogenetic reconstructions of the Paramyxoviridae and Rhabdoviridae families. P-gene-coding capacities were mapped onto the Paramyxoviridae phylogeny, and the most parsimonious path of multiple-coding-capacity evolution was determined. Levels of amino acid variation for Paramyxoviridae and Rhabdoviridae P-gene-encoded products were also analyzed. Proteins encoded in overlapping reading frames from the same nucleotides have different levels of amino acid variation. The nucleotide architecture that underlies the amino acid variation was determined in order to evaluate the role of selection in the evolution of the P gene overlapping reading frames. In every case, the evolution of one of the proteins encoded in the overlapping reading frames has been constrained by negative selection while the other has evolved more rapidly. The integrity of the overlapping reading frame that represents a derived state is generally maintained at the expense of the ancestral reading frame encoded by the same nucleotides. The evolution of such multicoding sequences is likely a response by RNA viruses to selective pressure to maximize genomic information content while maintaining small genome size. The ability to evolve such a complex genomic strategy is intimately related to the dynamics of the viral quasispecies, which allow enhanced exploration of the adaptive landscape.

Evidence for the recent horizontal transfer of long terminal repeat retrotransposon.

The evolutionary dynamics existing between transposable elements (TEs) and their host genomes have been likened to an "arms race." The selfish drive of TEs to replicate, in turn, elicits the evolution of host-mediated regulatory mechanisms aimed at repressing transpositional activity. It has been postulated that horizontal (cross-species) transfer may be one effective strategy by which TEs and other selfish genes can escape host-mediated silencing mechanisms over evolutionary time; however, to date, the most definitive evidence that TEs horizontally transfer between species has been limited to class II or DNA-type elements. Evidence that the more numerous and widely distributed retroelements may also be horizontally transferred between species has been more ambiguous. In this paper, we report definitive evidence for a recent horizontal transfer of the copia long terminal repeat retrotransposon between Drosophila melanogaster and Drosophila willistoni.

Genomic demography: a life-history analysis of transposable element evolution.

Retrotransposons are ubiquitous mobile genetic elements that have played a significant role in shaping eukaryotic genome evolution. The genome of the yeast Saccharomyces cerevisiae harbours five families of retrotransposons, Ty1-Ty5. With the publication of the S. cerevisiae genome sequence, for the first time a full genomic complement of retrotransposon sequences is available. Analysis of these sequences promises to yield insight into the nature of host--transposon coevolution. Evolutionary change in Ty elements depends on their replication and excision rates, which have been determined in the laboratory. Rates measured in the laboratory may differ from those that have operated over evolutionary time. Based on an analysis of sequence data for the Ty1, Ty2 and hybrid Ty1/2 families, we develop a novel 'genomic demography' model to estimate long-term transposition and excision rates and to estimate how long ago these elements entered the yeast genome. We find that rates of excision and transposition have averaged 7.2-8.7 x 10(-8) per generation over evolutionary time. Two separate models provide upper- and lower-bound estimates for the age of the system, suggesting that the first elements entered the genome between approximately 50 million and 250 million generations ago.

The role of interelement selection in Saccharomyces cerevisiae Ty element evolution.

Retrotransposons are mobile genetic elements that are ubiquitous components of eukaryotic genomes. The evolutionary success of retrotransposons is explained by their ability to replicate faster than the host genomes in which they reside. Elements with higher rates of genomic replication possess a selective advantage over less active elements. Retrotransposon populations, therefore, are shaped largely by selective forces acting at the genomic level between elements. To evaluate rigorously the effects of selective forces acting on retrotransposons, detailed information on the patterns of molecular variation within and between retrotransposon families is needed. The sequencing of the Saccharomyces cerevisiae genome, which includes the entire genomic complement of yeast retrotransposons, provides an unprecedented opportunity to access and analyze such data. In this study, we analyzed in detail the patterns of nucleotide variation within the open reading frames of two parental (Ty1 and Ty2) and one hybrid (Ty1/2) family of yeast retrotransposons. The pattern and distribution of nucleotide changes on the phylogenetic reconstructions of the three families of Ty elements reveal evidence of negative selection on both internal and external branches of the Ty phylogenies. These results indicate that most, if not all, Ty elements examined represent active or recently active retrotransposon lineages. We discuss the relevance of these findings with respect to the coevolutionary dynamic operating between genomic element populations and the host organisms in which they reside.

Tempo and mode of Ty element evolution in Saccharomyces cerevisiae.

The Saccharomyces cerevisiae genome contains five families of long terminal repeat (LTR) retrotransposons, Ty1-Ty5. The sequencing of the S. cerevisiae genome provides an unprecedented opportunity to examine the patterns of molecular variation existing among the entire genomic complement of Ty retrotransposons. We report the results of an analysis of the nucleotide and amino acid sequence variation within and between the five Ty element families of the S. cerevisiae genome. Our results indicate that individual Ty element families tend to be highly homogenous in both sequence and size variation. Comparisons of within-element 5' and 3' LTR sequences indicate that the vast majority of Ty elements have recently transposed. Furthermore, intrafamily Ty sequence comparisons reveal the action of negative selection on Ty element coding sequences. These results taken together suggest that there is a high level of genomic turnover of S. cerevisiae Ty elements, which is presumably in response to selective pressure to escape host-mediated repression and elimination mechanisms.

Comparative genomics and evolutionary dynamics of Saccharomyces cerevisiae Ty elements.

The availability of the complete genome sequence of Saccharomyces cerevisiae provides the unique opportunity to study an entire genomic complement of retrotransposons from an evolutionary perspective. There are five families of yeast retrotransposons, Ty1-Ty5. We have conducted a series of comparative sequence analyses within and among S. cerevisiae Ty families in an effort to document the evolutionary forces that have shaped element variation. Our results indicate that within families Ty elements vary little in terms of both size and sequence. Furthermore, intra-element 5'-3' long terminal repeat (LTR) sequence comparisons indicate that almost all Ty elements in the genome have recently transposed. For each family, solo LTR sequences generated by intra-element recombination far outnumber full length insertions. Taken together, these results suggest a rapid genomic turnover of S. cerevisiae Ty elements. The closely related Ty1 and Ty2 are the most numerous elements in the genome. Phylogenetic analysis of full length insertions reveals that reverse transcriptase mediated recombination between Ty1 and Ty2 elements has generated a number of hybrid Ty1/2 elements. These hybrid Ty1/2 elements have similar genomic structures with chimeric LTRs and chimeric TYB (pol) genes. Analysis of the levels of nonsynonymous (Ka) and synonymous (Ks) nucleotide variation indicates that Ty1 and Ty2 coding regions have been subject to strong negative (purifying) selection. Distribution of Ka and Ks on Ty1, Ty2 and Ty1/2 phylogenies reveals evidence of negative selection on both internal and external branches. This pattern of variation suggests that the majority of full length Ty1, Ty2 and Ty1/2 insertions represent active or recently active element lineages and is consistent with a high level of genomic turnover. The evolutionary dynamics of S. cerevisae Ty elements uncovered by our analyses are discussed with respect to selection among elements and the interaction between the elements and their host genome.

Interelement selection in the regulatory region of the copia retrotransposon.

We report the results of an analysis of naturally occurring cis-regulatory variation within and between two families of the copia Drosophila long terminal repeat (LTR) retrotransposon. The copia 5' LTR and adjacent untranslated leader region (ULR) consists of a number of well-characterized sequence motifs which play a role in regulating expression of the element. In order to understand the evolutionary forces which may be responsible for generating and maintaining copia regulatory sequence variation, we have quantified levels of naturally occurring copia LTR-ULR nucleotide variation and subjected the data to a series of tests of neutrality. Our analysis indicates that the copia LTR-ULR has been subject to negative purifying selection within families and positive adaptive selection between families. We discuss these findings with respect to the regulatory evolution of retrotransposons and the phenomenon of interelement selection.

Evolution of the copia retrotransposon in the Drosophila melanogaster species subgroup.

We report the results of a phylogenetic survey of the retrotransposon copia in the melanogaster subgroup of the Drosophila genus. The polymerase chain reaction was used to amplify the copia 5' long terminal repeat and the adjacent untranslated leader region from representative melanogaster subgroup species. Restriction and sequence analyses of this region reveal discrete classes of copia size variants within the melanogaster subgroup. Phylogenetic comparisons of copia sequence data indicate that the size variants represent different copia subfamilies which diverged prior to their distribution in the melanogaster subgroup. Our results also suggest that copia elements have been subject to horizontal and vertical transmission during their evolution.

Evidence for the role of recombination in the regulatory evolution of Saccharomyces cerevisiae Ty elements.

The recent completion of the sequencing of the Saccharomyces cerevisiae genome provides a unique opportunity to analyze the evolutionary relationships existing among the entire complement of retrotransposons residing within a single genome. In this article we report the results of such an analysis of two closely related families of yeast long terminal repeat (LTR) retrotransposons, Ty1 and Ty2. In our study, we analyzed the molecular variation existing among the 32 Ty1 and 13 Ty2 elements present within the S. cerevisiae genome recently sequenced within the context of the yeast genome project. Our results indicate that while the Ty1 family is most likely ancestral to Ty2 elements, both families of elements are relatively recent components of the S. cerevisiae genome. Our results also indicate that both families of elements have been subject to purifying selection within their protein coding regions. Finally, and perhaps most interestingly, our results indicate that a relatively recent recombination event has occurred between Ty2 and a subclass of Ty1 elements involving the LTR regulatory region. We discuss the possible biological significance of these findings and, in particular, how they contribute to a better overall understanding of LTR retrotransposon evolution.