Multiple genome comparison within a bacterial species reveals a unit of evolution spanning two adjacent genes in a tandem paralog cluster

It has been assumed that an open reading frame (ORF) represents a unit of gene evolution as well as a unit of gene expression and function. In the present work, we report a case in which a unit comprising the 3' region of an ORF linked to a downstream intergenic region that is in turn linked to the 5' region of a downstream ORF has been conserved, and has served as the unit of gene evolution. The genes are tandem paralogous genes from the bacterium Staphylococcus aureus, for which more than ten entire genomes have been sequenced. We compared these multiple genome sequences at a locus for the lpl (lipoprotein-like) cluster (encoding lipoprotein homologs presumably related to their host interaction) in the genomic island termed nuSaalpha. A highly conserved nucleotide sequence found within every lpl ORF is likely to provide a site for homologous recombination. Comparison of phylogenies of the 5'-variable region and the 3'-variable region within the same ORF revealed significant incongruence. In contrast, pairs of the 3'-variable region of an ORF and the 5'-variable region of the next downstream ORF gave more congruent phylogenies, with distinct groups of conserved pairs. The intergenic region seemed to have coevolved with the flanking variable regions. Multiple recombination events at the central conserved region appear to have caused various types of rearrangements among strains, shuffling the two variable regions in one ORF, but maintaining a conserved unit comprising the 3'-variable region, the intergenic region, and the 5'-variable region spanning adjacent ORFs. This result has strong impact on our understanding of gene evolution because most gene lineages underwent tandem duplication and then diversified. This work also illustrates the use of multiple genome sequences for high-resolution evolutionary analysis within the same species.     

Expansion and divergence of the GH locus between spider monkey and chimpanzee

Growth hormone (GH) has been previously described as showing distinct evolutionary stories between primates and other mammals. A burst of changes and successive amplification events took place in the primate lineage giving rise to a multigene family in the three Anthropoidea lineages. Polymerase chain reaction (PCR) was used to obtain the genes and the intergenic regions comprising the GH loci of the spider monkey (Ateles geoffroyi), a New-World primate, and of the chimpanzee (Pan troglodytes), an ape. The intergenic sequences of both species were screened by hybridization to detect copies of the Alu family, which have been implicated in the formation of the human GH locus. The GH locus of the spider monkey contains at least six GH-related genes, four of them were cloned. Likewise, five short intergenic sequences of approximately 3 kb were amplified and cloned. On the other hand, in the chimpanzee four new placental lactogen (PL) genes as well as four intergenic regions were amplified. Consequently, in this ape, six genes (two GHs, previously obtained, and four PLs) are clustered, separated by intergenic sequences of different lengths (two short ones of about 5 kb, and at least two long ones between 9 and 13 kb). The presence of Alu sequences within the intergenic regions of both GH loci corroborates the current hypothesis that they acted as a driving force for the locus expansion. GH sequence comparisons reveal that several gene-conversion events might have occurred during the formation of this genome region, which has undergone independent evolution in the three Anthropoidea branches. To establish the GH's evolutionary history may prove to be a difficult task due to these gene-conversion events.     

Structural organization of the barley D-hordein locus in comparison with its orthologous regions of wheat genomes.

D hordein, a prolamin storage protein of barley endosperms, is highly homologous to the high molecular weight (HWM) glutenin subunits, which are the major determinants of bread-making quality in wheat flour. In hexaploid wheat (AABBDD), each genome contains two paralogous copies of HMW-glutenin genes that encode the x- and y-type HMW-glutenin subunits. Previously, we reported the sequence analysis of a 102-kb genomic region that contains the HMW-glutenin locus of the D genome from Aegilops tauschii, the donor of the D genome of hexaploid wheat. Here, we present the sequence analysis of a 120-kb D-hordein region of the barley genome, a more distantly related member of the Triticeae grass tribe. Comparative sequence analysis revealed that gene content and order are generally conserved. Genes included in both of these orthologous regions are arranged in the following order: a Xa21-like receptor kinase, an endosperm globulin, an HMW prolamin, and a serine (threonine) protein kinase. However, in the wheat D genome, a region containing both the globulin and HMW-glutenin gene was duplicated, indicating that this duplication event occurred after the separation of the wheat and barley genomes. The intergenic regions are divergent with regard to the sequence and structural organization. It was found that different types of retroelements are responsible for the intergenic structure divergence in the wheat and barley genomes. In the barley region, we identified 16 long terminal repeat (LTR) retrotransposons in three distinct nested clusters. These retroelements account for 63% of the contig sequence. In addition, barley D hordein was compared with wheat HMW glutenins in terms of cysteine residue conservation and repeat domain organization.     

Evolution of a functionally related lactate dehydrogenase and pyruvate decarboxylase pseudogene complex in maize.

A large proportion of the maize genome is repetitive DNA (60-80%) with retrotransposons contributing significantly to the repetitive DNA component. The majority of retrotransposon DNA is located in intergenic regions and is organized in a nested fashion. Analysis of an 8.2-kb segment of maize genomic DNA demonstrated the presence of three retrotransposons of different reiteration classes in addition to lactate dehydrogenase and pyruvate decarboxylase pseudogenes. Both of the pseudogenes were located within a defective retrotransposon element (LP-like element) which possessed identical long terminal repeats (LTRs) with inverted repeats at each end, a primer binding site, a polypurine tract, and generated a 5-bp target site duplication. A model describing the events leading to the formation of the LP-like element is proposed.     

Rapid genome evolution revealed by comparative sequence analysis of orthologous regions from four triticeae genomes

Bread wheat (Triticum aestivum) is an allohexaploid species, consisting of three subgenomes (A, B, and D). To study the molecular evolution of these closely related genomes, we compared the sequence of a 307-kb physical contig covering the high molecular weight (HMW)-glutenin locus from the A genome of durum wheat (Triticum turgidum, AABB) with the orthologous regions from the B genome of the same wheat and the D genome of the diploid wheat Aegilops tauschii (Anderson et al., 2003; Kong et al., 2004). Although gene colinearity appears to be retained, four out of six genes including the two paralogous HMW-glutenin genes are disrupted in the orthologous region of the A genome. Mechanisms involved in gene disruption in the A genome include retroelement insertions, sequence deletions, and mutations causing in-frame stop codons in the coding sequences. Comparative sequence analysis also revealed that sequences in the colinear intergenic regions of these different genomes were generally not conserved. The rapid genome evolution in these regions is attributable mainly to the large number of retrotransposon insertions that occurred after the divergence of the three wheat genomes. Our comparative studies indicate that the B genome diverged prior to the separation of the A and D genomes. Furthermore, sequence comparison of two distinct types of allelic variations at the HMW-glutenin loci in the A genomes of different hexaploid wheat cultivars with the A genome locus of durum wheat indicates that hexaploid wheat may have more than one tetraploid ancestor.     

Human U1 small nuclear RNA genes: extensive conservation of flanking sequences suggests cycles of gene amplification and transposition.

The DNA immediately flanking the 164-base-pair U1 RNA coding region is highly conserved among the approximately 30 human U1 genes. The U1 multigene family also contains many U1 pseudogenes (designated class I) with striking although imperfect flanking homology to the true U1 genes. Using cosmid vectors, we now have cloned, characterized, and partially sequenced three 35-kilobase (kb) regions of the human genome spanning U1 homologies. Two clones contain one true U1 gene each, and the third bears two class I pseudogenes 9 kb apart in the opposite orientation. We show by genomic blotting and by direct DNA sequence determination that the conserved sequences surrounding U1 genes are much more extensive than previously estimated: nearly perfect sequence homology between many true U1 genes extends for at least 24 kb upstream and at least 20 kb downstream from the U1 coding region. In addition, the sequences of the two new pseudogenes provide evidence that class I U1 pseudogenes are more closely related to each other than to true genes. Finally, it is demonstrated elsewhere (Lindgren et al., Mol. Cell. Biol. 5:2190-2196, 1985) that both true U1 genes and class I U1 pseudogenes map to chromosome 1, but in separate clusters located far apart on opposite sides of the centromere. Taken together, these results suggest a model for the evolution of the U1 multigene family. We speculate that the contemporary family of true U1 genes was derived from a more ancient family of U1 genes (now class I U1 pseudogenes) by gene amplification and transposition. Gene amplification provides the simplest explanation for the clustering of both U1 genes and class I pseudogenes and for the conservation of at least 44 kb of DNA flanking the U1 coding region in a large fraction of the 30 true U1 genes.     

Further evidence of microcolinearity between barley and rice genomes at two orthologous regions

Two genetic markers, BCD135 and RZ567 were used to select clones from genomic BAC libraries of barley and rice for sequencing and subsequent sequence comparisons. A set of two orthologous BACs each from barley and rice was selected by hybridization with BCD135 and RZ567 cDNA probes. A total of 556-kb stretch including two barley BACs (773K135 and 745C13) and two orthologous rice BACs (24K23 and 49D11) was completely sequenced. Comparative sequence analysis between orthologous BACs from the two species revealed presence of two conserved genes at BCD135 region and only one gene at the RZ567 regions. The two conserved genes were in the same order and orientation in both the species however, separated by significantly larger distance in barley. The larger distance between two barley genes was mainly due to presence of different retrotransposable elements and their derivatives (78.9% of the intergenic region) that expanded the barley BCD135 region at the rate of 9.1X. An additional gene of unknown function was also inserted along with several retrotransposable elements between two conserved genes at barley BCD135 region. More genome expansion rate (10X) around barley RZ567 locus was estimated by extremely high proportion (> 70%) of retrotransposons. Among different retrotransposons, the Sabrina elements rather than BARE were more prevalent in both the regions. Contrary to it, the BCD135 region of rice was composed of only 17.1% retrotransposable elements and no significant retrotransposons except 14 miniature inverted transposable elements (MITEs) were observed in its RZ567 region. The sequence comparison between orthologous regions of rice and barley genomes was useful for gene identification and determination of individual gene structure indicating the possibility of effective utilization of rice genome sequences in understanding the large genome of barley. (The sequence data described in this paper have been submitted to the GenBank data library under the accession no. AF474072 (773K14), AF474071 (745C13), AF480497 (24K23) and AF480496 (49D11)).     

Genomic and phylogenic comparisons of the alpha-globin and beta-globin intergenic sequences between zebra fish (Danio rerio) and six closely related Cyprinindae species

The alpha-globin and beta-globin genes of the zebrafish are tightly linked on the same chromosome in a 3'-5' and 5'-3' configuration, respectively. Although the location of the controlling sequences has been mapped to the intergenic region, analysis determining the uniqueness of this unusual arrangement to zebrafish has not been undertaken. To explore this, we isolated, sequenced, and phylogenetically analyzed the intergenic region between globin gene families of seven Cyprinindae species including zebrafish. These species were grouped into an in group (immediate relatives, not so distant relatives), and an out group (distant relative). Cellulose acetate electrophoresis of hemoglobin (Hb) detected multiple variants in each species, but a band with electrophoretic mobility (EM) of 6.7 x 10(-5) cm(2).volt(-1).sec(-1) was shared between species. Polymerase chain reaction (PCR) amplification of the intergenic globin gene region also detected a 1.0-kb fragment that was repeated in the in group and a 1.2-kb fragment in the out group. Sequence comparison confirmed that the genetic orientation and controlling sequences location were conserved throughout this region in all seven species. This phylogenic footprinting indicated that the configuration was not exclusive to zebrafish. To confirm sequence alignment, maximum parsimony phylogenic analysis, was performed. Only one member of that group the giant danio, was not closely clustered, being located almost equidistance between the immediate relative and the species of the other clusters. This may represent an ancestral configuration prior to transposition of the alpha globin and beta-globin genes families to nonsynteny.     

Structural characterization of the dihydroflavonol 4-reductase B (DFR-B) gene in the sweet potato.

A genomic fragment containing the dihydroflavonol 4-reductase B (DFR-B) gene was cloned from the sweet potato (Ipomoea batatas) and its nucleotide sequence was analyzed. The exons and flanking regions were highly homologous to those of previously reported DFR-B genes of the Japanese morning glory, whereas the introns and the intergenic region were less conserved. In addition to the sequences of three miniature inverted-repeat transposable elements (MITEs) and one direct repeat previously reported in the DFR-B gene of Japanese morning glory, two mobile element-like sequences were newly identified in the sweet potato DFR-B gene. At least four allelic sequences were found to exist by amplification of the DFR-B gene from various sweet potato cultivars. One of these allelic sequences had a 2-kb deletion in the intergenic region and was observed in the cultivars with high anthocyanin content in their storage roots.     

Comparison of the maxicircle (mitochondrial) genomes of Leishmania tarentolae and Trypanosoma brucei at the level of nucleotide sequence

The entire 16.7-kilobase (kb) transcribed region of the Leishmania tarentolae maxicircle was compared to the entire 15-kb transcribed region of the Trypanosoma brucei maxicircle at the nucleotide sequence level by dot matrix analysis and by alignments of individual genes. The L. tarentolae NADH dehydrogenase subunit 1 (ND1) gene was identified in a newly obtained 2.9-kb sequence. All but two regions which flank the cytochrome b gene are highly conserved in both species. One 3.1-kb region in L. tarentolae that contains the cytochrome oxidase subunit III (COIII) gene and several open reading frames corresponds to a 2-kb sequence in T. brucei with limited sequence homology that lacks the COIII gene. Another 0.6-kb region that comprises an unidentified open reading frame (open reading frame 12) in L. tarentolae is substituted by a nonhomologous 0.4-kb open reading frame in T. brucei. A short intergenic region between the ND1 gene and the maxicircle unidentified reading frame 1 gene shows limited sequence homology, and the regions between the ND4 and ND5 genes and between the COI and ND4 genes are not conserved. All of the intergenic regions share G + C richness and a similar pattern of G versus C strand bias. 1.8 kb of the L. tarentolae divergent region (DV) and around 3 kb of the T. brucei DV were also obtained. The T. brucei DV sequences were not homologous to the L. tarentolae DV sequence but were organized in a similar fashion with tandem repeats of varying complexity.     

PseudoPipe: an automated pseudogene identification pipeline

MOTIVATION: Mammalian genomes contain many 'genomic fossils' i.e. pseudogenes. These are disabled copies of functional genes that have been retained in the genome by gene duplication or retrotransposition events. Pseudogenes are important resources in understanding the evolutionary history of genes and genomes. RESULTS: We have developed a homology-based computational pipeline ('PseudoPipe') that can search a mammalian genome and identify pseudogene sequences in a comprehensive and consistent manner. The key steps in the pipeline involve using BLAST to rapidly cross-reference potential "parent" proteins against the intergenic regions of the genome and then processing the resulting "raw hits" -- i.e. eliminating redundant ones, clustering together neighbors, and associating and aligning clusters with a unique parent. Finally, pseudogenes are classified based on a combination of criteria including homology, intron-exon structure, and existence of stop codons and frameshifts.     

Examination of four HLA class I pseudogenes. Common events in the evolution of HLA genes and pseudogenes.

The HLA class I gene family in lymphoblastoid cell line 721 has been studied in detail and a number of sequences in addition to the classical genes have been identified. The cloning, characterization, and nucleotide sequences of four sequences, all full length HLA class I pseudogenes, are described in this report. These pseudogenes, contained within 5.4-, 5.9-, 7.0-, and 9.2-kb HindIII fragments, each have the class I exon-intron structure as well as class I homology in their 5' and 3' flanking regions. However, all four sequences have one or more substitutions that perturb the coding region, leaving little doubt that they are in fact pseudogenes. Comparisons among these sequences and the HLA class I genes revealed that their homology with the class I genes is patchwork. Thus, although some regions have diverged, other contiguous intron-exon sequences are highly conserved. Comparisons in the 5' regions indicate that the pseudogene promoters more closely resemble the classical HLA promoters than the nonclassical promoters as none of the unique structural features found in the HLA-E, -F, or -G regulatory regions are present in any of the pseudogene promoters. Further comparisons revealed that at least two putative gene conversion events, similar to those hypothesized to have occurred in the evolution of some HLA genes, may have occurred in the evolution of some of the pseudogenes. These and other hypothetical events in the evolution of the class I gene family are discussed.     

Gene-containing regions of wheat and the other grass genomes

Deletion line-based high-density physical maps revealed that the wheat (Triticum aestivum) genome is partitioned into gene-rich and -poor compartments. Available deletion lines have bracketed the gene-containing regions to about 10% of the genome. Emerging sequence data suggest that these may further be partitioned into "mini" gene-rich and gene-poor regions. An average of about 10% of each gene-rich region seem to contain genes. Sequence analyses in various species suggest that uneven distribution of genes may be a characteristic of all grasses and perhaps all higher organisms. Comparison of the physical maps with genetic linkage maps showed that recombination in wheat and barley (Hordeum vulgare) is confined to the gene-containing regions. Number of genes, gene density, and the extent of recombination vary greatly among the gene-rich regions. The gene order, relative region size, and recombination are highly conserved within the tribe Triticeae and moderately conserved within the family. Gene-poor regions are composed of retrotransposon-like non-transcribing repeats and pseudogenes. Direct comparisons of orthologous regions indicated that gene density in wheat is about one-half compared with rice (Oryza sativa). Genome size difference between wheat and rice is, therefore, mainly because of amplification of the gene-poor regions. Presence of species-, genera-, and family-specific repeats reveal a repeated invasion of the genomes by different retrotransposons over time. Preferential transposition to adjacent locations and presence of vital genes flanking a gene-rich region may have restricted retrotransposon amplification to gene-poor regions, resulting into tandem blocks of non-transcribing repeats. Insertional inactivation by adjoining retro-elements and selection seem to have played a major role in stabilizing genomes.