Comparison of Pax1/9 locus reveals 500-Myr-old syntenic block and evolutionary conserved noncoding regions

Identification of conserved genomic regions within and between different genomes is crucial when studying genome evolution. Here, we described regions of strong synteny conservation between vertebrate deuterostomes (tetrapods and teleosts) and invertebrate deuterostomes (amphioxus and sea urchin). The shared gene contents across phylogenetically distant species demonstrate that the conservation of the regions stemmed from an ancestral segment instead of a series of independent convergent events. Comparison of the syntenic regions allows us to postulate the primitive gene organization in the last common ancestor of deuterostomes and the evolutionary events that occurred to the 3 distinct lineages of sea urchin, amphioxus, and vertebrates after their separation. In addition, alignment of the syntenic regions led to the identification of 8 noncoding evolutionarily conserved regions shared between amphioxus and vertebrates. To our knowledge, this is the first report of conserved noncoding sequences shared by vertebrates and nonvertebrates. These noncoding sequences have high possibility of being elements that regulate neighboring genes. They are likely to be a factor in the maintenance of conserved synteny over long phylogenetic distance in different deuterostome lineages.     

Evolution of aspartyl proteases by gene duplication: the mouse renin gene is organized in two homologous clusters of four exons.

Overlapping recombinant clones that appear to encompass the entire renin gene, named Ren 1, have been isolated from a library of BALB/c mouse genomic DNA fragments. Based on restriction endonuclease mapping and DNA sequence analysis, Ren 1 spans 9.6 kb and contains nine exons interrupted by eight intervening sequences of highly variable size. The first exon, encoding the signal peptide of preprorenin, is separated from the eight following exons by a 3-kb intron. These eight exons are organized into two clusters of four separated by a 2-kb intron. DNA stretches encoding the aspartyl residues, which are part of the active site of renin, are located at homologous positions in both clusters. Our results show that aspartyl protease genes have arisen by duplication and fusion of an ancestral gene containing five exons. The estimated date of the duplication event of the mouse renin genes Ren 1 and Ren 2 is discussed.     

The primary structure and gene organization of human substance P and neuromedin K receptors.

The gene organization and amino acid sequences of human substance P and neuromedin K receptors (SPR and NKR, respectively) are reported on the basis of molecular cloning and sequence determination of genomic DNA containing the respective receptor gene. The human SPR and NKR genes, unlike many other genes for G-protein-coupled receptors, (G protein, guanyl-nucleotide-binding-regulatory protein), contain introns which interrupt the protein-coding regions into 5 exons. The human SPR and NKR genes extend over 60 kb and 45 kb, respectively and are considerably larger than the human substance K receptor (SKR) gene consisting of 12 kb. All 4 introns, however, are located at equivalent positions of the three tachykinin receptor genes, suggesting that they evolved from a common ancestral gene. Human SPR and NKR consist of 407 and 465 amino acid residues, respectively, each possessing structural features characteristic of the members of G-protein-coupled receptors. The human and rat receptors show a common tendency of distinctly segmented sequence conservation and divergence among the three receptors, and the observed sequence conservation and divergence would contribute to the emergence of similar but distinct properties of the three receptors. Furthermore, the amino acid sequences and the gene sizes are more closely related between SPR and NKR than between SKR and NKR, suggesting that the SPR gene evolved from the primordial NKR gene after a gene duplication to form the NKR and SKR genes.     

Molecular evolution of the mammalian ribosomal protein gene, RPS14

Ribosomal protein S14 genes (RPS14) in eukaryotic species from protozoa to primates exhibit dramatically different intron-exon structures yet share homologous polypeptide-coding sequences. To recognize common features of RPS14 gene architectures in closely related mammalian species and to evaluate similarities in their noncoding DNA sequences, we isolated the intron-containing S14 locus from Chinese hamster ovary (CHO) cell DNA by using a PCR strategy and compared it with human RPS14. We found that rodent and primate S14 genes are composed of identical protein-coding exons interrupted by introns at four conserved DNA sites. However, the structures of corresponding CHO and human RPS14 introns differ significantly. Nonetheless, individual intron splice donor, splice acceptor, and upstream flanking motifs have been conserved within mammalian S14 homologues as well as within RPS14 gene fragments PCR amplified from other vertebrate genera (birds and bony fish). Our data indicate that noncoding, intronic DNA sequences within highly conserved, single-copy ribosomal protein genes are useful molecular landmarks for phylogenetic analysis of closely related vertebrate species.      

Alpha-tropomyosin gene organization. Alternative splicing of duplicated isotype-specific exons accounts for the production of smooth and striated muscle isoforms

We have previously isolated and characterized cloned complementary DNAs (cDNAs) for striated and smooth muscle alpha-tropomyosin. The sequences of these cDNA clones suggested that these two isoforms were encoded by the same gene. Here, we have determined the complete structure of the alpha-tropomyosin (alpha-TM) gene, establishing that a single gene, with a sequence complexity of 28 kilobase pairs, is split into 12 exons and produces the smooth and striated muscle alpha-TM mRNA isoforms by alternative splicing of a minimum of five exchangeable isotype-specific exons. The elucidation of the intron/exon organization of alpha-TM suggests that this gene evolved from an ancestral gene encoding a 21-aa protein that might represent the primordial actin binding domain. Sequence comparison between the pairs of exons coding for the "isotype switch regions" and among the corresponding regions of tropomyosin genes in a variety of species ranging from insects to mammals, suggests that the alternatively spliced exons are very old and might have arisen before the radiation of the arthropods, more than 600 million years ago. Additionally, the examination of the intronic sequences has uncovered potential alternative intramolecular secondary structures (hairpin-loop structures) which might be involved in the tissue-specific expression of the duplicated and mutually exclusive alpha-TM isotype-specific exons.     

Structure, organization, and chromosomal location of the gene encoding a form of rice soluble starch synthase.

A rice (Oryza sativa L.) genomic clone encoding the gene for a form of soluble starch synthase (SSS1) and its 5'- and 3'-flanking regions has been isolated and sequenced. The SSS1 gene contained 15 exons interrupted by 14 introns. The exon/intron organization of the SSS1 gene was divergent from that of the rice Waxy gene coding for granule-bound starch synthase, thus suggesting that the SSS1 and granule-bound starch synthase genes have evolved from an ancestral gene in a different way or that the two genes are products of different ancestral genes that have converged during evolution. However, these two genes were closely located to each other on rice chromosome 6 at an approximate map distance of 5 centimorgans. The nucleotide sequence of the 5'-end region of the gene is unique because of the presence of some repetitive sequences.     

The gene map of the Norway rat (Rattus norvegicus) and comparative mapping with mouse and man

The current status of the rat gene map is presented. Mapping information is now available for a total of 214 loci and the number of mapped genes is increasing steadily. The corresponding number of loci quoted at HGM10 was 128. Genes have been assigned to 20 of the 22 chromosomes in the rat. Some aspects of comparative mapping with mouse and man are also discussed. It was found that there is a good correlation between the morphological homologies detectable in rat and mouse chromosomes, on the one hand, and homology at the gene level on the other. For 10 rat synteny groups all the genes so far mapped are syntenic also in the mouse. For the remaining rat synteny groups it appears that the majority of the genes will be syntenic on specific (homologous) mouse chromosomes, with only a few genes dispersed to other members of the mouse karyotype. Furthermore, the data indicate that mouse chromosome 1 genetically corresponds to two rat chromosomes, viz., 9 and 13, equalizing the difference in chromosome number between the two species. Further mappings will show whether the genetic homology will prove to be as extensive as these preliminary results indicate. As might be expected from evolutionary considerations, rat synteny groups are much more dispersed in the human genome. It is clear, however, that many groups of genes have remained syntenic during the period since man and rat shared a common ancestor. One further point was noted. In two cases groups of genes were syntenic in the mouse but dispersed to two chromosomes in rat and man, whereas in a third case a group of genes was syntenic in the rat but dispersed to two chromosomes in mouse and man. This finding argues in favor of the notion that the original gene groups were on separate ancestral chromosomes, which have fused in one rodent species but remained separate in the other and in man.     

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     

Comparative analysis of the structural organization of two closely related vitellogenin genes in X. laevis

The structural organization of the two closely related vitellogenin genes A1 and A2 has been determined and compared by electron microscopy. In both genes the mRNA-coding sequence of 6 kb is interrupted 33 times, leading to a total gene length of 21 kb for gene A1 and 16 kb for gene A2. Thus both genes have a mean exon length of 0.175 kb, while the mean intron length is 0.45 kb in gene A1 and 0.31 kb in gene A2. Because the introns interrupt the structural sequence at homologous positions in genes A1 and A2, we suggest that these two genes are the products of a duplication of an ancestral gene which had an intron-exon arrangement similar to that of the extant genes. Since the duplication event, the sequence and length of the analogous introns have changed rapidly, whereas homologous exons have diverged to an extent of only 5% of their sequences. The results suggest different mechanisms of evolution for exons and introns. While the exons evolved primarily by point mutations, such mutations, as well as deletion, insertion and duplication events, were important in the evolution of the introns.     

The ovalbumin gene family: Structure of the X gene and evolution of duplicated split genes

The X, Y and ovalbumin genes, which are found within a 40 kb region of the chicken genome, are all expressed in oviduct under steroid hormone control, and share some sequence homologies. We have now cloned the complete X gene and have analyzed its structure. It codes for two RNA species, X and X′; both are coded by eight exons and appear to differ only by the size of their 3′ untranslated region, X′ RNA being 1400 nucleotides longer than X RNA. The striking similarity in the number and length of the exons which constitute the X, Y or ovalbumin genes establishes that they have evolved from a common ancestor gene by duplication events. Comparison of selected regions of the X and ovalbumin genes indicates that the exon sequences coding for protein and the location of the splice junctions have been well-conserved. The introns and the 3′ untranslated exonic sequences have diverged much more rapidly. Four regions of apparently unrelated repetitive sequences are found both outside the X gene and within it (in two introns and in the sequence coding for the 3′ untranslated part of X′RNA). The intragenic repetitive sequences have no counterpart in the ovalbumin and Y genes.     

Frequent appearance of novel protein-coding sequences by frameshift translation

Genomic duplication, followed by divergence, contributes to organismal evolution. Several mechanisms, such as exon shuffling and alternative splicing, are responsible for novel gene functions, but they generate homologous domains and do not usually lead to drastic innovation. Major novelties can potentially be introduced by frameshift mutations and this idea can explain the creation of novel proteins. Here, we employ a strategy using simulated protein sequences and identify 470 human and 108 mouse frameshift events that originate new gene segments. No obvious interspecies overlap was observed, suggesting high rates of acquisition of evolutionary events. This inference is supported by a deficiency of TpA dinucleotides in the protein-coding sequences, which decreases the occurrence of translational termination, even on the complementary strand. Increased usage of the TGA codon as the termination signal in newer genes also supports our inference. This suggests that tolerated frameshift changes are a prevalent mechanism for the rapid emergence of new genes and that protein-coding sequences can be derived from existing or ancestral exons rather than from events that result in noncoding sequences becoming exons.     

Recognition of unknown conserved alternatively spliced exons

The split structure of most mammalian protein-coding genes allows for the potential to produce multiple different mRNA and protein isoforms from a single gene locus through the process of alternative splicing (AS). We propose a computational approach called UNCOVER based on a pair hidden Markov model to discover conserved coding exonic sequences subject to AS that have so far gone undetected. Applying UNCOVER to orthologous introns of known human and mouse genes predicts skipped exons or retained introns present in both species, while discriminating them from conserved noncoding sequences. The accuracy of the model is evaluated on a curated set of genes with known conserved AS events. The prediction of skipped exons in the approximately 1% of the human genome represented by the ENCODE regions leads to more than 50 new exon candidates. Five novel predicted AS exons were validated by RT-PCR and sequencing analysis of 15 introns with strong UNCOVER predictions and lacking EST evidence. These results imply that a considerable number of conserved exonic sequences and associated isoforms are still completely missing from the current annotation of known genes. UNCOVER also identifies a small number of candidates for conserved intron retention.     

Structure of human immunoglobulin gamma genes: implications for evolution of a gene family

We have cloned five human immunoglobulin gamma genes from a fetal liver gene library. Four of them encode the known human immunoglobulin gamma chains gamma 1, gamma 2, gamma 3 and gamma 4. A fifth gamma gene seems to be a pseudogene. Nucleotide sequence determination demonstrates that the gamma 3 gene contains four separate hinge exons. Comparison of these hinge exons with those of the other gamma genes indicates that the first hinge exon is homologous to that of the pseudogene, and that the other three hinge exons are homologous to that of the gamma 1 gene, suggesting that the gamma 3 gene ancestor is a hybrid gene created by unequal crossing-over between the ancestral gamma 1 and psi gamma genes. Amplification of the gamma 1-type hinge exon probably followed to complete the gamma 3 gene. This hypothesis inevitably postulates the gene order 5'-gamma 1-gamma 3-psi gamma-3'. Cloning of overlapping chromosomal segments demonstrates that the gamma 2 gene is located 19 kb 5' to the gamma 4 gene. These analyses indicate that the human gamma-gene family has evolved by several types of DNA rearrangemet, including duplication of a complete gene; duplication of a hinge exon; and reassortment of exons by unequal cross-over between two adjacent genes.     

Patterns of Y and X chromosome DNA sequence divergence during the Felidae radiation.

The 37 species of modern cats have evolved from approximately eight phylogenetic lineages within the past 10 to 15 million years. The Felidae family has been described with multiple measures of morphologic and molecular evolutionary methods that serve as a framework for tracking gene divergence during brief evolutionary periods. In this report, we compare the mode and tempo of evolution of noncoding sequences of a large intron within Zfy (783 bp) and Zfx (854 bp), homologous genes located on the felid Y and X chromosomes, respectively. Zfy sequence variation evolves at about twice the rate of Zfx, and both gene intron sequences track feline hierarchical topologies accurately. As homoplasies are infrequent in patterns of nucleotide substitution, the Y chromosome sequence displays a remarkable degree of phylogenetic consistency among cat species and provides a highly informative glimpse of divergence of sex chromosome sequences in Felidae.     

A comparative study on the genes for three porins of the Escherichia coli outer membrane. DNA sequence of the osmoregulated ompC gene

The DNA sequence of the ompC gene which encodes one of the outer membrane porins has been determined. The gene appears to encode a secretory precursor of OmpC protein consisting of a total of 367 amino acid residues with a signal peptide of 21 amino acid residues at its NH2-terminal end. The 5' end noncoding region including the promoter of the ompC gene is extremely [A-T]-rich, and the codon usage in the ompC gene is unusual as are those in genes for other abundant outer membrane proteins. The promoter sequence of the ompC gene was compared with that of the ompF gene, both of which are controlled by the osmoregulatory operon, ompB. The deduced amino acid sequence of the OmpC protein showed extensive homology with that of the other porins (OmpF and PhoE proteins). The homology in the primary amino acid sequences, as well as the coding DNA sequences among the porins, indicates that the structural genes for the three porins evolved from a common ancestral gene. Comparison of the amino acid sequences among the OmpC, OmpF, and PhoE porins will be discussed with regard to structure and function.     

Rapid evolutionary dynamics in a 2.8‐Mb chromosomal region containing multiple prolamin and resistance gene families in Aegilops tauschii

Prolamin and resistance gene families are important in wheat food use and in defense against pathogen attacks, respectively. To better understand the evolution of these multi‐gene families, the DNA sequence of a 2.8‐Mb genomic region, representing an 8.8 cM genetic interval and harboring multiple prolamin and resistance‐like gene families, was analyzed in the diploid grass Aegilops tauschii, the D‐genome donor of bread wheat. Comparison with orthologous regions from rice, Brachypodium, and sorghum showed that the Ae. tauschii region has undergone dramatic changes; it has acquired more than 80 non‐syntenic genes and only 13 ancestral genes are shared among these grass species. These non‐syntenic genes, including prolamin and resistance‐like genes, originated from various genomic regions and likely moved to their present locations via sequence evolution processes involving gene duplication and translocation. Local duplication of non‐syntenic genes contributed significantly to the expansion of gene families. Our analysis indicates that the insertion of prolamin‐related genes occurred prior to the separation of the Brachypodieae and Triticeae lineages. Unlike in Brachypodium, inserted prolamin genes have rapidly evolved and expanded to encode different classes of major seed storage proteins in Triticeae species. Phylogenetic analyses also showed that the multiple insertions of resistance‐like genes and subsequent differential expansion of each R gene family. The high frequency of non‐syntenic genes and rapid local gene evolution correlate with the high recombination rate in the 2.8‐Mb region with nine‐fold higher than the genome‐wide average. Our results demonstrate complex evolutionary dynamics in this agronomically important region of Triticeae species.     

Exon-intron organization of the human gp130 gene

The exon-intron organization and sequences of the exon-intron boundaries of the human gp130 transmembrane receptor gene have been determined using genomic DNAs as samples. The gp130 gene comprises 17 exons and 16 introns. The positions of the exon-intron boundaries show good correlation to the functional/homology regions of gp130. Exons 3-17 code for the gp130 protein, and each subdomain of the receptor is encoded by a set of exons. The coding potential of exons and the intron phasing of the human gp130 gene conform to the patterns observed previously for other cytokine receptor genes. This supports the notions that the gp130 gene evolved from the same ancestral gene that gave rise to other members of the cytokine receptor family.     

Linkage arrangement of human placental lactogen and growth hormone genes

Human placental lactogen (hPL) and growth hormone (hGH) are two hormones thought to have evolved from a common ancestral gene (along with prolactin), yet they have quite different functions and specificities. The nucleic acid sequences of the respective cDNAs of the two genes share considerable homology, as well as the existence of multiple forms of each gene within the genome. In this study we report on the linkage arrangement of several genes from this group. Two hPL-like genes as well as an hGH gene are shown to be linked within a 38-kilobase pair region of DNA. Linkage between a variant hGH gene and an hPL gene is also shown. The orientation and structural organization of these genes was previously established using 5'- and 3'-specific probes from a placental lactogen cDNA clone and detailed restriction endonuclease mapping. Restriction fragments from the overlapping clones were verified by comparison to digests of high molecular weight genomic DNA. In addition, the location of a specific class of repetitive DNA sequences, the Alu family, was mapped on these clones using the recombinant clone BLUR 8. All members of this multigene family have Alu repeat sequences either immediately flanking their 3' or 5' untranslated regions or within their intervening sequences.