Assignment of orthologous relationships among mammalian alpha-globin genes by examining flanking regions reveals a rapid rate of evolution

In order to study the relationships among mammalian alpha-globin genes, we have determined the sequence of the 3' flanking region of the human alpha 1 globin gene and have made pairwise comparisons between sequenced alpha-globin genes. The flanking regions were examined in detail because sequence matches in these regions could be interpreted with the least complication from the gene duplications and conversions that have occurred frequently in mammalian alpha-like globin gene clusters. We found good matches between the flanking regions of human alpha 1 and rabbit alpha 1, human psi alpha 1 and goat I alpha, human alpha 2 and goat II alpha, and horse alpha 1 and goat II alpha. These matches were used to align the alpha-globin genes in gene clusters from different mammals. This alignment shows that genes at equivalent positions in the gene clusters of different mammals can be functional or nonfunctional, depending on whether they corrected against a functional alpha-globin gene in recent evolutionary history. The number of alpha-globin genes (including pseudogenes) appears to differ among species, although highly divergent pseudogenes may not have been detected in all species examined. Although matching sequences could be found in interspecies comparisons of the flanking regions of alpha- globin genes, these matches are not as extensive as those found in the flanking regions of mammalian beta-like globin genes. This observation suggests that the noncoding sequences in the mammalian alpha-globin gene clusters are evolving at a faster rate than those in the beta-like globin gene clusters. The proposed faster rate of evolution fits with the poor conservation of the genetic linkage map around alpha-globin gene clusters when compared to that of the beta-like globin gene clusters. Analysis of the 3' flanking regions of alpha-globin genes has revealed a conserved sequence approximately 100-150 bp 3' to the polyadenylation site; this sequence may be involved in the expression or regulation of alpha-globin genes.      

Evolution of genes for allelic and isotypic forms of immunoglobulin kappa chains and of the genes for T-cell receptor beta chains in rabbits

New insights into the evolution of the families of genes encoding immunoglobulins and T-cell receptors of rabbits (Oryctolagus cuniculus) have come from molecular genetic studies. In contrast to human and mouse, rabbits were shown to have two genes for the constant region of immunoglobulin light chains (C kappa 1 and C kappa 2 isotypes) and complex allelic variants of K1 (allotypes). Although K1 allotype protein sequences differed at up to 41% of the amino acid positions, 3' untranslated, 5', and 3' flanking regions were conserved, and in the coding regions 78-80% of the codons with differences had replacement changes. Proportions of silent changes and changes in noncoding regions were comparable. Thus, in spite of their markedly different protein sequences, the K1b4, b5, and b9 allotypes appeared to be products of allelic genes. Molecular genetic analyses suggested that they may have undergone rapid divergence after an ancestral K2-like gene duplicated. Some rabbits were found to have two similar T-cell receptor C beta genes as do humans and many strains of mice, but others appeared to have three different C beta. In addition, we found allotypic forms of C beta. Some of the C beta allotypic differences occurred at positions where analogous C kappa allotypic differences were found. We also found V beta in mouse and human that were more similar to rabbit V beta than closely linked rabbit genes were to each other. This contrasts with rabbit immunoglobulin VH gene sequences that reflect concerted evolution. The data suggested that T-cell receptor V beta genes duplicated prior to mammalian radiation.     

Purifying and directional selection in overlapping prokaryotic genes

In overlapping genes, the same DNA sequence codes for two proteins using different reading frames. Analysis of overlapping genes can help in understanding the mode of evolution of a coding region from noncoding DNA. We identified 71 pairs of convergent genes, with overlapping 3' ends longer than 15 nucleotides, that are conserved in at least two prokaryotic genomes. Among the overlap regions, we observed a statistically significant bias towards the 123:132 phase (i.e. the second codon base in one gene facing the degenerate third position in the second gene). This phase ensures the least mutual constraint on nonconservative amino acid replacements in both overlapping coding sequences. The excess of this phase is compatible with directional (positive) selection acting on the overlapping coding regions. This could be a general evolutionary mode for genes emerging from noncoding sequences, in which the protein sequence has not been subject to selection.     

The structure and evolution of the spider monkey delta-globin gene

We have isolated the delta-globin gene of the New-World spider monkey, Ateles geoffroyi, and compared its nucleotide sequence with those of other primate delta- and beta-globin genes. Among primate delta-globin genes, the rate of nonsynonymous substitutions is much less than the rate of synonymous substitutions. This suggests that primate delta- globin genes may remain under evolutionary conservation, perhaps because hemoglobin A2 has an as yet unknown physiological importance.      

The evolution of a plant globin gene family

Summary We have analyzed the sequences of soybean leghemoglobin genes as an initial step toward understanding their mode of evolution. Alignment of the sequences of plant globin genes with those of animals reveals that (i) based on the proportion of nucleotide substitutions that have occurred at the first, second, and third codon positions, the time of divergence of plant and animal globin gene families appears to be extremely remote (between 900 million and 1.4 billion years ago, if one assumes constancy of evolutionary rate in both the plant and animal lineages) and (ii) in addition to the normal regulatory sequences on the 5 end, an approximately 30-base-pair sequence, specific to globin genes, that surrounds the cap site is conserved between the plant and animal globin genes. Comparison of the leghemoglobin sequences with one another shows that (i) the relative amount of sequence divergence in various coding and noncoding regions is roughly similar to that found for animal globin genes and (ii) as in animal globin genes, the positions of insertions and deletions in the intervening sequences often coincide with the locations of direct repeats. Thus, the mode of evolution of the plant globin genes appears to resemble, in many ways, that of their animal counterparts. We contrast the overall intergenic organization of the plant globin genes with that of animal genes, and discuss the possibility of the concerted evolution of the leghemoglobin genes.     

Sequencing and comparative analysis of a conserved syntenic segment in the Solanaceae

Comparative genomics is a powerful tool for gaining insight into genomic function and evolution. However, in plants, sequence data that would enable detailed comparisons of both coding and noncoding regions have been limited in availability. Here we report the generation and analysis of sequences for an unduplicated conserved syntenic segment (CSS) in the genomes of five members of the agriculturally important plant family Solanaceae. This CSS includes a 105-kb region of tomato chromosome 2 and orthologous regions of the potato, eggplant, pepper, and petunia genomes. With a total neutral divergence of 0.73-0.78 substitutions/site, these sequences are similar enough that most noncoding regions can be aligned, yet divergent enough to be informative about evolutionary dynamics and selective pressures. The CSS contains 17 distinct genes with generally conserved order and orientation, but with numerous small-scale differences between species. Our analysis indicates that the last common ancestor of these species lived approximately 27-36 million years ago, that more than one-third of short genomic segments (5-15 bp) are under selection, and that more than two-thirds of selected bases fall in noncoding regions. In addition, we identify genes under positive selection and analyze hundreds of conserved noncoding elements. This analysis provides a window into 30 million years of plant evolution in the absence of polyploidization.     

Concerted evolution of the primate immunoglobulin alpha-gene through gene conversion.

We determined four nucleotide sequences of the hominoid immunoglobulin alpha (C alpha) genes (chimpanzee C alpha 2, gorilla C alpha 2, and gibbon C alpha 1 and C alpha 2 genes), which made possible the examination of gene conversions in all hominoid C alpha genes. The following three methods were used to detect gene conversions: 1) phenetic tree construction; 2) detection of a DNA segment with extremely low variability between duplicated C alpha genes; and 3) a site by site search of shared nucleotide changes between duplicated C alpha genes. Results obtained from method 1 indicated a concerted evolution of the duplicated C alpha genes in the human, chimpanzee, gorilla, and gibbon lineages, while results obtained from method 2 suggested gene conversions in the human, gorilla, and gibbon C alpha genes. With method 3 we identified clusters of shared nucleotide changes between duplicated C alpha genes in human, chimpanzee, gorilla, and gibbon lineages, and in their hypothetical ancestors. In the present study converted regions were identified over the entire C alpha gene region excluding a few sites in the coding region which have escaped from gene conversion. This indicates that gene conversion is a general phenomenon in evolution, that can be clearly observed in non-functional regions.     

Gene conversions in the growth hormone gene family of primates: stronger homogenizing effects in the Hominidae lineage

In humans, the growth hormone/chorionic somatomammotropin gene family is composed of five highly similar genes. We characterized the gene conversions that occurred between the growth hormone genes of 11 primate species. We detected 48 conversions using GENECONV and others were only detected using phylogenetic analyses. Gene conversions were detected in all species analyzed, their average size (± standard deviation) is 197.8 ± 230.4 nucleotides, the size of the conversions is correlated with sequence similarity and converted regions are significantly more GC-rich than non-converted regions. Gene conversions have a stronger homogenizing effect in Hominidae genes than in other primate species. They are also less frequent in conserved gene regions and towards functionally important genes. This suggests that the high degree of sequence similarity observed between the growth hormone genes of primate species is a consequence of frequent gene conversions in gene regions which are under little selective constraints.     

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.     

Molecular evolution of the 5'-flanking regions of the duplicated Amy genes in Drosophila melanogaster species subgroup

The nucleotide sequences of the 5'-flanking regions of the duplicated Amy genes in eight sibling species belonging to the melanogaster species subgroup are analyzed. In Drosophila melanogaster, a region of about 450 bp immediately upstream of the translation initiation site of the two paralogous genes (the proximal and distal genes) has sequence similarities. However, we could not detect any significant sequence similarity in the region more upstream than -450. This result indicates that the coding regions of the ancestral Amy gene were duplicated together with 450 bp of the 5'-flanking region as one unit. Multiple alignment of these 450-bp sequences in the proximal and distal genes of all eight species revealed a mosaic pattern of highly conserved and divergent regions. The conserved regions included almost all the putative regulatory elements identified in previous analyses of the sequences. A phylogenetic analysis of the aligned sequences shows that these 450-bp sequences are clustered into the proximal and the distal groups. As a whole, the divergence between groups in this region is very large in contrast to that in the coding regions. Based on the divergence between groups, the 450-bp region is divided into two subregions. We found that the ratios of the divergence between groups to that within groups differ in the two subregions. From these observations, we discuss a possibility of positive selection acting on the subregion immediately upstream of the Amy coding region to cause divergence of regulatory elements of the paralogous genes.      

Evidence of Gene Conversion Events Between Paralogous Sequences Produced by Tetraploidization in Salmoninae Fish

We investigated the occurrence of gene conversions between paralogous sequences of Salmoninae derived from ancestral tetraploidization and their effect on the evolutionary history of DNA sequences. A microsatellite with long flanking regions (750 bp) including both coding and noncoding sequences was analyzed. Microsatellite size polymorphism was used to detect the alleles of both paralogous counterparts and infer linkage arrangement between loci. DNA sequencing of seven Salmoninae species revealed that paralogous sequences were highly differentiated within species, especially for noncoding regions. Ten gene conversion events between paralogous sequences were inferred. While these events appears to have homogenized regions of otherwise highly differential paralogous sequences, they amplified the differentiation among orthologous sequences. Their effects were larger on coding than on noncoding regions. As a consequence, noncoding sequences grouped by orthologous lineages in phylogenetic trees, whereas coding regions grouped by taxa. Based upon these results, we present a model showing how gene conversion events may also result in the PCR amplification of nonorthologous sequences in different taxa, with obvious complications for phylogenetic inferences, comparative mapping, and population genetic studies. Received: 11 October 2000 / Accepted: 18 September 2001     

Regulation, linkage, and sequence of mouse metallothionein I and II genes.

The mouse metallothionein II (MT-II) gene is located approximately 6 kilobases upstream of the MT-I gene. A comparison of the sequences of mouse MT-I and MT-II genes (as well as those of other mammals) reveals that the coding regions are highly conserved even at "silent" positions but that the noncoding regions and introns are extremely divergent between primates and rodents. There are four blocks of conserved sequences in the promoters of mouse MT-I, mouse MT-II, and human MT-IIA genes; one includes the TATAAA sequence, and another has been implicated in regulation by heavy metals. Mouse MT-I and MT-II mRNAs are induced to approximately the same extent in vivo in response to cadmium, dexamethasone, or lipopolysaccharide. Mouse MT-I and MT-II genes are regulated by metals but not by glucocorticoids after transfection into HeLa cells.     

A novel gene family NBPF: intricate structure generated by gene duplications during primate evolution

Partial and complete genome duplications occurred during evolution and resulted in the creation of new genes and gene families. We identified a novel and intricate human gene family located primarily in regions of segmental duplications on human chromosome 1. We named it NBPF, for neuroblastoma breakpoint family, because one of its members is disrupted by a chromosomal translocation in a neuroblastoma patient. The NBPF genes have a repetitive structure with high intragenic and intergenic sequence similarity in both coding and noncoding regions. These similarities might expose these genomic regions to illegitimate recombination, resulting in structural variation in the NBPF genes. The encoded proteins contain a highly conserved domain of unknown function, which we have named the NBPF repeat. In silico analysis combined with the isolation of multiple full-length cDNA clones showed that several members of this gene family are abundantly expressed in a large variety of tissues and cell lines. Strikingly, no discernable orthologues could be identified in the completed genomes of fruit fly, nematode, mouse, or rat, but sequences with low homology could be isolated from the draft canine and bovine genomes. Interestingly, this gene family shows primate-specific duplications that result in species-specific arrays of NBPF homologous sequences. Overall, this novel NBPF family reflects the continuous evolution of primate genomes that resulted in large physiological differences, and its potential role in this process is discussed.     

Evidence for Widespread Positive and Negative Selection in Coding and Conserved Noncoding Regions of Capsella grandiflora

The extent that both positive and negative selection vary across different portions of plant genomes remains poorly understood. Here, we sequence whole genomes of 13 Capsella grandiflora individuals and quantify the amount of selection across the genome. Using an estimate of the distribution of fitness effects, we show that selection is strong in coding regions, but weak in most noncoding regions, with the exception of 5′ and 3′ untranslated regions (UTRs). However, estimates of selection on noncoding regions conserved across the Brassicaceae family show strong signals of selection. Additionally, we see reductions in neutral diversity around functional substitutions in both coding and conserved noncoding regions, indicating recent selective sweeps at these sites. Finally, using expression data from leaf tissue we show that genes that are more highly expressed experience stronger negative selection but comparable levels of positive selection to lowly expressed genes. Overall, we observe widespread positive and negative selection in coding and regulatory regions, but our results also suggest that both positive and negative selection on plant noncoding sequence are considerably rarer than in animal genomes.     

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.      

Identification and characterization of new long conserved noncoding sequences in vertebrates

Comparative sequence analyses have identified highly conserved genomic DNA sequences, including noncoding sequences, between humans and other species. By performing whole-genome comparisons of human and mouse, we have identified 611 conserved noncoding sequences longer than 500 bp, with more than 95% identity between the species. These long conserved noncoding sequences (LCNS) include 473 new sequences that do not overlap with previously reported ultraconserved elements (UCE), which are defined as aligned sequences longer than 200 bp with 100% identity in human, mouse, and rat. The LCNS were distributed throughout the genome except for the Y chromosome and often occurred in clusters within regions with a low density of coding genes. Many of the LCNS were also highly conserved in other mammals, chickens, frogs, and fish; however, we were unable to find orthologous sequences in the genomes of invertebrate species. In order to examine whether these conserved sequences are functionally important or merely mutational cold spots, we directly measured the frequencies of ENU-induced germline mutations in the LCNS of the mouse. By screening about 40.7 Mb, we found 35 mutations, including mutations at nucleotides that were conserved between human and fish. The mutation frequencies were equivalent to those found in other genomic regions, including coding sequences and introns, suggesting that the LCNS are not mutational cold spots at all. Taken together, these results suggest that mutations occur with equal frequency in LCNS but are eliminated by natural selection during the course of evolution.