Impact of transposable elements on the evolution of mammalian gene regulation

Transposable elements (TEs) are present in all organisms and nearly half of the human and mouse genome is derived from ancient transpositions. This fact alone suggests that TEs have played a major role in genome organization and evolution. Studies undertaken over the last two decades or so clearly show that TEs of various kinds have played an important role in organism evolution. Here we review the impact TEs have on the evolution of gene regulation and gene function with an emphasis on humans. Understanding the mechanisms resulting in genomic change is central to our understanding of gene regulation, genetic disease and genome evolution. Full comprehension of these biological processes is not possible without an in depth knowledge of how TEs impact upon the genome.     

Vertebrate protamine gene evolution I. Sequence alignments and gene structure

Summary The availability of the amino acid sequence for nine different mammalian P1 family protamines and the revised amino acid sequence of the chicken protamine galline (Oliva and Dixon 1989) reveals a much close relationship between mammalian and avian protamines than was previously thought (Nakano et al. 1976). Dot matrix analysis of all protamine genes for which genomic DNA or cDNA sequence is available reveals both marked sequence similarities in the mammalian protamine gene family and internal repeated sequences in the chicken protamine gene. The detailed alignments of the cis-acting regulatory DNA sequences shows several consensus sequence patterns, particularly the conservation of a cAMP response element (CRE) in all the protamine genes and of the regions flanking the TATA box, CAP site, N-terminal coding region, and polyadenylation signal. In addition we have found a high frequency of the CA dinucleotide immediately adjacent to the CRE element of both the protamine genes and the testis transition proteins, a feature not present in other genes, which suggests the existence of an extended CRE motif involved in the coordinate expression of protamine and transition protein genes during spermatogenesis. Overall these findings suggest the existence of an avian-mammalian P1 protamine gene line and are discussed in the context of different hypotheses for protamine gene evolution and regulation.     

The globin gene family of the cephalochordate amphioxus: implications for chordate globin evolution

Background  

The lancelet amphioxus (Cephalochordata) is a close relative of vertebrates and thus may enhance our understanding of vertebrate gene and genome evolution. In this context, the globins are one of the best studied models for gene family evolution. Previous biochemical studies have demonstrated the presence of an intracellular globin in notochord tissue and myotome of amphioxus, but the corresponding gene has not yet been identified. Genomic resources of Branchiostoma floridae now facilitate the identification, experimental confirmation and molecular evolutionary analysis of its globin gene repertoire.     

The sequence of the gorilla fetal globin genes: evidence for multiple gene conversions in human evolution

Two fetal globin genes (G gamma and A gamma) from one chromosome of a lowland gorilla (Gorilla gorilla gorilla) have been sequenced and compared to three human loci (a G gamma-gene and two A gamma-alleles). A comparison of regions of local homology among these five sequences indicates that long after the duplication that produced the two nonallelic gamma-globin loci of catarrhine primates, about 35 million years (Myr) ago, at least one gene conversion event occurred between these loci. This conversion occurred not long before the ancestral divergence (about 6 Myr ago) of Homo and Gorilla. After this ancestral divergence, a minimum of three more gene conversion events occurred in the human lineage. Each human A gamma-allele shares specific sequence features with the gorilla A gamma-gene; one such distinctive allelic feature involves the simple repeated sequence in IVS 2. This suggests that early in the human lineage the A gamma-genes may have undergone a crossing-over event mediated by this simple repeated sequence. The DNA sequences from coding regions of both G gamma- and A gamma-loci, a comparison of 292 codons in the corresponding gorilla and human genes, show an unusually low evolutionary rate, with only two nonsilent differences and, surprisingly, not even one silent substitution. The two nonsynonymous substitutions observed predict a glycine at codon 73 and an arginine at codon 104 in the gorilla A gamma-sequence rather than aspartic acid and lysine, respectively, in human A gamma. Because only arginine has been found at position 104 in gamma-chains of Old World monkeys, it may represent the ancestral residue lost in gorilla and human G gamma-chains and in the human A gamma-chain. Possibly the arginine codon (AGG) was replaced by the lysine codon (AAG) in the G gamma-gene of a common ancestor of Homo and Gorilla and then was transferred to the A gamma-gene by subsequent conversions in the human lineage. DNA sequence conversions, similar to that attributed to the fetal gamma-globin genes, appear to be relatively frequent phenomena and, if widespread throughout the genome, may have profound evolutionary consequences.      

Retrotransposons and the evolution of mammalian gene expression

Transposable elements, and retroviral-like elements in particular, are a rich potential source of genetic variation within a host's genome. Many mutations of endogenous genes in phylogenetically diverse organisms are due to insertion of elements that affect gene expression by altering the normal pattern of regulation. While few such associations are known to have been maintained over time, two recently elucidated examples suggest transposable elements may have a significant impact in evolution of gene expression. The first example, concerning the mouse sex-limited protein (Slp), clearly establishes that ancient retroviral enhancer sequences now confer hormonal dependence on the adjacent gene. The second example shows that within the human amylase gene family, salivary specific expression has arisen due to inserted sequences, deriving perhaps from a conjunction of two retrotransposable elements.     

The genesis and evolution of homeobox gene clusters

Once called the 'Rosetta stone' of developmental biology, the homeobox continues to fascinate both evolutionary and developmental biologists. The birth of the homeotic, or Hox, gene cluster, and its subsequent evolution, has been crucial in mediating the major transitions in metazoan body plan. Comparative genomics studies indicate that the more recently discovered ParaHox and NK clusters were linked to the Hox cluster early in evolution, and that together they constituted a 'megacluster' of homeobox genes that conspicuously contributed to body-plan evolution.     

Nucleotide sequence analysis of the lemur beta-globin gene family: evidence for major rate fluctuations in globin polypeptide evolution

Lemur beta-related globin genes have been isolated and sequenced. Orthology of prosimian and human epsilon-, gamma-, and beta-related globin genes was established by dot-matrix analysis. All of these lemur globin genes potentially encode functional beta-related globin polypeptides, though precisely when the gamma-globin gene is expressed remains unknown. The organization of the 18-kb brown lemur beta-globin gene cluster (5' epsilon-gamma-[psi eta-delta]-beta 3') is consistent with its evolution by contraction via unequal crossing-over from the putative ancestral mammalian beta-globin gene cluster (5' epsilon-gamma- eta-delta-beta 3'). The dwarf lemur nonadult globin genes are arranged as in the brown lemur. Similar levels of synonymous (silent) nucleotide substitutions and noncoding DNA sequence differences have accumulated between species in all of these genes, suggesting a uniform rate of noncoding DNA divergence throughout primate beta-globin gene clusters. These differences are comparable with those observed in the nonfunctional psi eta pseudogene and have therefore accumulated at the presumably maximal neutral rate. In contrast, nonsynonymous (replacement) nucleotide substitutions show a significant heterogeneity in distribution for both the same gene in different lineages and different genes in the same lineage. These major fluctuations in replacement but not silent substitution rates cannot be attributed to changes in mutation rate, suggesting that changes in the rate of globin polypeptide evolution in primates is not governed solely by variable mutation rates.      

Interacting gene clusters and the evolution of the vertebrate immune system

Unraveling the "code" of genome structure is an important goal of genomics research. Colocalization of genes in eukaryotic genomes may facilitate preservation of favorable allele combinations between epistasic loci or coregulation of functionally related genes. However, the presence of interacting gene clusters in the human genome has remained unclear. We systematically searched the human genome for evidence of closely linked genes whose protein products interact. We find 83 pairs of interacting genes that are located within 1 Mbp in the human genome or 37 if we exclude hub proteins. This number of interacting gene clusters is significantly more than expected by chance and is not the result of tandem duplications. Furthermore, we find that these clusters are significantly more conserved across vertebrate (but not chordate) genomes than other pairs of genes located within 1 Mbp in the human genome. In many cases, the genes are both present but not clustered in older vertebrate lineages. These results suggest gene cluster creation along the human lineage. These clusters are not enriched for housekeeping genes, but we find a significant contribution from genes involved in "response to stimulus." Many of these genes are involved in the immune response, including, but not limited to, known clusters such as the major histocompatibility complex. That these clusters were formed contemporaneously with the origin of adaptive immunity within the vertebrate lineage suggests that novel evolutionary and regulatory constraints were associated with the operation of the immune system.     

Buffering and the evolution of chromosome-wide gene regulation

Copy number variation (CNV) in terms of aneuploidies of both entire chromosomes and chromosomal segments is an important evolutionary driving force, but it is inevitably accompanied by potentially problematic variations in gene doses and genomic instability. Thus, a delicate balance must be maintained between mechanisms that compensate for variations in gene doses (and thus allow such genomic variability) and selection against destabilizing CNVs. In Drosophila, three known compensatory mechanisms have evolved: a general segmental aneuploidy-buffering system and two chromosome-specific systems. The two chromosome-specific systems are the male-specific lethal complex, which is important for dosage compensation of the male X chromosome, and Painting of fourth, which stimulates expression of the fourth chromosome. In this review, we discuss the origin and function of buffering and compensation using Drosophila as a model.     

Two kinetic methods to study the regulation of mammalian hexokinases

Two new kinetic analyses for mammalian hexokinases are presented, which permit one to study the regulation of these enzymes by product inhibition. One method uses the pyruvate kinase-coupled assay and the other the glucose-6-phosphate dehydrogenase-coupled assay. Both methods give simple linear plots, which indicate that the magnesium-ATP complex overcomes the glucose 6-phosphate inhibition competitively, but by atypical kinetics. A new regulation coefficient (Kr) was defined and it was shown that, with both assay methods, the reciprocals of the slopes of the simple linear plots are proportional to Kr.[Mg.ATP].NADP, but not NAD, was found to be a powerful inhibitor of pig heart hexokinase.     

A novel approach to local reliability of sequence alignments

MOTIVATION: The pairwise alignment of biological sequences obtained from an algorithm will in general contain both correct and incorrect parts. Hence, to allow for a valid interpretation of the alignment, the local trustworthiness of the alignment has to be quantified. RESULTS: We present a novel approach that attributes a reliability index to every pair of residues, including gapped regions, in the optimal alignment of two protein sequences. The method is based on a fuzzy recast of the dynamic programming algorithm for sequence alignment in terms of mean field annealing. An extensive evaluation with structural reference alignments not only shows that the probability for a pair of residues to be correctly aligned grows consistently with increasing reliability index, but moreover demonstrates that the value of the reliability index can directly be translated into an estimate of the probability for a correct alignment.     

Concerted evolution led to high expression of a prosimian primate delta globin gene locus

The delta globin gene in simian primates is either weakly expressed (in hominoids and New World monkeys) or silent (in Old World monkeys). In prosimian primates, however, an unequal homologous crossover between the psi eta and delta loci of lemurs produced a hybrid psi eta delta pseudogene locus, whereas in tarsier the delta locus encodes a beta-type chain found in 18% of adult tarsier hemoglobin molecules. In the present study, the nucleotide and amino acid sequences of the galago delta and beta globin genes and their encoded peptides were determined, and evidence is provided showing that the galago delta locus encodes a beta-type chain (beta 2) found in 40% of the galago fetal and postnatal hemoglobin molecules, whereas the beta locus encodes the remaining 60% of the beta-type chain (beta 1). Galago beta 1 and beta 2 chains differ from each other by only one amino acid residue. The homology between the galago delta and beta loci extends from 800 base pairs 5' of the proximal CCAAT element to near the end of exon 3 as a result of a recombination event in which beta sequence replaced delta sequence. After this initial recombination event, concerted evolution between the loci continued over their conserved coding, intron 1, and promoter regions but failed to occur between the two loci in their intron 2 and distal 5'-flanking sequences where the two loci have now diverged by 20%. Calculations based on this divergence value and on a rate of noncoding sequence evolution of 4.2 x 10(-9) to 5.5 x 10(-9) substitutions/site/year for the lorisiform lineage to galago yielded a date of 18-24 million years ago for the initial recombination event. The fact that the promoter sequences of the galago delta locus are the same as that of the galago beta locus may account for the high level of expression of the galago delta gene.     

Organization and evolution of mitochondrial gene clusters in human

Currently, the spatial patterns of mitochondrial genes and how the genomic localization of (pseudo)genes originated from mitochondrial DNA remain largely unexplained. The aim of this study was to elucidate the organization of mitochondrial (pseudo)genes given their evolutionary origin. We used a keyword finding method and a bootstrapping method to estimate parameter values that represent the distribution pattern of mitochondrial genes in the nuclear genome. Almost half of mitochondrial genes showing physical clusters were located in the pericentromeric and subtelomeric regions of the chromosome. Most interestingly, the size of these clusters ranged from 0.085 to 3.2 Mb (average ± SD 1.3 ± 0.73 Mb), which coincides with the size of the evolutionary pocket, or the average size of evolutionary breakpoint regions. Our findings imply that the localization of mitochondrial genes in the human genome is determined independent of adaptation.