Isolation and characterization of the beta and epsilon subunit genes of mouse muscle acetylcholine receptor

The genes coding for the beta and epsilon subunits of the mouse muscle nicotinic acetylcholine receptor (nAChR) were mapped by Southern blot analysis, and the entire loci for both genes cloned. The results indicate that they are single-copy genes. Both were sequenced to determine their size and structural organization. The beta subunit gene spans approximately 8 kilobases and is organized into 11 exons. A region containing cysteines, which are thought to form a disulfide bond and which are highly conserved, is encoded by one exon in all muscle acetylcholine receptor genes with the exception of the beta subunit gene, where it is split into two exons. The epsilon subunit gene spans 4.3 kilobases and contains 12 exons; it has the same structure as the gamma and delta nAChR genes. The intron-exon boundaries and exonic organization of the five known nAChR genes were compared. The analysis showed that the first 4 exons and the last exon of all muscle and brain nAChR subunit genes have the same boundaries, with the exception of a nAChR-related gene in Drosophila.     

Using ESTs to improve the accuracy of de novo gene prediction

Background  

ESTs are a tremendous resource for determining the exon-intron structures of genes, but even extensive EST sequencing tends to leave many exons and genes untouched. Gene prediction systems based exclusively on EST alignments miss these exons and genes, leading to poor sensitivity. De novo gene prediction systems, which ignore ESTs in favor of genomic sequence, can predict such "untouched" exons, but they are less accurate when predicting exons to which ESTs align. TWINSCAN is the most accurate de novo gene finder available for nematodes and N-SCAN is the most accurate for mammals, as measured by exact CDS gene prediction and exact exon prediction.     

Correlation of gene and protein structure of rat and human lipocortin I

Lipocortins (annexins) are a family of calcium-dependent phospholipid-binding proteins with phospholipase A2 inhibitory activity. The characteristic primary structure of members of this family consists of a core structure of four or eight repeated domains, which have been implicated in calcium-dependent phospholipid binding. In two lipocortins (I and II) a short amino-terminal sequence distinct from the core structure has potential regulatory functions which are dependent on its phosphorylation state. We have isolated the rat and the human lipocortin I genes and found that they both consist of 13 exons with a striking conservation of their exon-intron structure and their promoter and amino acid sequences. Both lipocortin I genes are at least 19 kbp in length with exons ranging from 57 to 123 bp interrupted by introns as large as 5 kbp. Each of the four repeat units of lipocortin I are encoded by two consecutive exons while individual exons code for the highly conserved putative calcium-binding domains. The promoter sequences in the rat and in human genes are highly conserved and contain nucleotide sequences characterized as enhancer sequences in other genes. The structure of the lipocortin I gene lends support to the hypothesis that the lipocortin genes arose by a duplication of a single domain.     

Structure and evolution of the bovine prothrombin gene

The cloned bovine prothrombin gene has been characterized by partial DNA sequence analysis, including the 5' and 3' flanking sequences and all the intron-exon junctions. The gene is approximately 15.4 x 10(3) base-pairs in length and comprises 14 exons interrupted by 13 introns. The exons coding for the prepro-leader peptide and the gamma-carboxyglutamic acid-containing region are similar in organization to the corresponding exons in the factor IX and protein C genes. This region has probably evolved as a result of recent gene duplication and exon shuffling events. The exons coding for the kringles and the serine protease region of the prothrombin gene are different in organization from the homologous regions in other genes, suggesting that introns have been inserted into these regions after the initial gene duplication events.     

Gene Conversion and Evolution of Daphniid Hemoglobins (Crustacea, Cladocera)

The extracellular hemoglobins of cladocerans derive from the aggregation of 12 two-domain globin subunits that are apparently encoded by four genes. This study establishes that at least some of these genes occur as a tandem array in both Daphnia magna and Daphnia exilis. The genes share a uniform structure; a bridge intron separates two globin domains which each include three exons and two introns. Introns are small, averaging just 77 bp, but a longer sequence (2.2–3.2 kb) separates adjacent globin genes. A survey of structural diversity in globin genes from other daphniids revealed three independent cases of intron loss, but exon lengths were identical, excepting a 3-bp insertion in exon 5 of Simocephalus. Heterogeneity in the extent of nucleotide divergence was marked among exons, largely as a result of the pronounced diversification of the terminal exon. This variation reflected, in part, varying exposure to concerted evolution. Conversion events were frequent in exons 1–4 but were absent from exons 5 and 6. Because of this difference, the results of phylogenetic analyses were strongly affected by the sequences employed in this construction. Phylogenies based on total nucleotide divergence in exons 1–4 revealed affinities among all genes isolated from a single species, reflecting the impact of gene conversion events. In contrast, phylogenies based on total nucleotide divergence in exons 5 and 6 revealed affinities among orthologous genes from different taxa. Received: 8 March 1999 / Accepted: 14 July 1999     

Origin of structural domains of the serum-albumin gene family and a predicted structure of the gene for vitamin D-binding protein

We have recently determined complete DNA sequences for the human albumin and alpha-fetoprotein [AFP] genes and thus have identified their detailed structures. Each is composed of three domains of four exons, three of which are internal and one of which is a domain-linking exon. Equivalent exons in each domain show sufficient sequence and structural similarity to be considered homologous; additional unique exons at each end of the gene show no similarity to the internal triplicated structures. Since earlier, conflicting evolutionary models were based on analysis of single gene structures, we derived from five genes a series of consensus sequences representing the three internal exons as well as the domain-linking exon. The five genes were human and rat albumin and human, mouse, and rat AFP genes. Structurally equivalent exons of the different domains are shown to have arisen from a single exon in a one-domain precursor. Exons that bridge the domains arose from an unequal crossover that fused two exons of the precursor. Our model suggests that part of the coding sequence of the one-domain precursor may have been derived from an intron, by way of loss of a splice site. The consensus sequences were used to propose an intron-exon structure for the related gene encoding the serum vitamin D-binding protein (DBP). DBP is truncated relative to albumin and AFP, and we submit that this results from deletion of two internal exons in the third domain of the gene rather than from premature termination of the coding sequence.     

Characterization of bovine keratin genes: similarities of exon patterns in genes coding for different keratins.

Four different genomic clones which contain the genes coding for epidermal keratins Ia (mol. wt. approximately 68 000), Ib (68 000), III (60 000) and VIb (54 500) have been selected using cDNA probes and identified by hybrid-selection translation. The genes vary considerably in length, primarily due to differences in intron sizes: keratin Ia, 9.3 kb (approximately 2.55 kb total exons); keratin Ib, 6.0 kb (2.25 kb exons); keratin III, 6.0 kb (2.2 kb exons); keratin VIb, 4.4 kb (1.85 kb exons). The genes for all three representatives of the basic (type II) cytokeratin subfamily, i.e., keratins Ia, Ib and III, contain eight introns of variable sizes (0.1-1.8 kb) and their exon patterns are very similar. The gene coding for keratin VIb, a representative of the acidic (type I) subfamily, contains seven introns, and the size pattern of its five innermost exons closely resembles that of the genes of the type II keratins. Most of the introns are located in regions coding for the alpha-helical cores of these proteins. Mapping of the intron positions by the S1 nuclease technique and sequencing of some exon-intron boundaries has revealed that some of the introns of all four keratin genes have similar positions to each other and to those of the hamster vimentin gene.(ABSTRACT TRUNCATED AT 250 WORDS)     

The universe of exons revisited

We study the distribution of exons in eukaryotic genes to determine whether one can detect the reuse of exon sequences and to use the frequency of such reuse to estimate how many ancestral exon sequences there might have been. We use two databases of exons. One contained 56,276 internal exons from putatively unrelated genes (less than 20% sequence identity) and the second contained 8917 internal exons from regions of these genes that are homologous and colinear with prokaryotic genes; these are ancient conserved regions (ACRs). At the 95% significance level we find 3500 exon-sequence matches in the large database and 500 matches in the ACR database. These matches correspond to groups of similar sequences. The size–rank relationship for these groups follows a power law, the size falling off as the inverse square root of the rank. This form of the power law distribution leads us to make an estimate for the size of a possible universe of ancestral exons. Using the data corresponding to the ACR regions, that universe is estimated to be about 15,000–30,000 in size.     

Gene conversion and natural selection in the evolution of X-linked color vision genes in higher primates

During higher primate evolution, gene conversion seems to have occurred often between the red and green photo-pigment genes, which are tandemly linked on the X chromosome. To understand this phenomenon better, intron 4 sequences of the red and green pigment genes of a male human (an Asian Indian), a male chimpanzee, and a male baboon were amplified by PCR and sequenced. The data show that the intron 4 sequences between the two genes have been strongly or completely homogenized in the three species studied. Apparently recent gene conversion events have occurred in introns 4 of the red and green pigment genes in humans and chimpanzees. Two or more conversion events may have occurred at different times in introns 4 of the two pigment genes in baboons. The divergence between the two genes is significantly lower in intron 4 than in exons 4 and 5 in each species, contrary to the usual situation that introns evolve faster than exons. It is most likely that strong natural selection for maintaining the distinct functions of exons 4 and 5 of the red and green pigment genes has acted against sequence homogenization of these exons.      

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 prediction of exons through an analysis of spliceable open reading frames.

We have developed a computer program which predicts internal exons from naive genomic sequence data and which will run on any IBM-compatible 80286 (or higher) computer. The algorithm searches a sequence for 'spliceable open reading frames' (SORFs), which are open reading frames bracketed by suitable splice-recognition sequences, and then analyzes the region for codon usage. Potential exons are stratified according to the reliability of their prediction, from confidence levels 1 to 5. The program is designed to predict internal exons of length greater than 60 nucleotides. In an analysis of 116 genes of a training set, 384 out of 441 such exons (87.1%) are identified, with 280 (63.5%) of predictions matching the true exon exactly (at both 5' and 3' splice junctions and in the correct reading frame), and with 104 (23.6%) exons matching partially. In a similar analysis of 14 genes in a test set unrelated to the genes used to generate the parameters of the program, 70 out of 80 internal exons greater than 60 bp in length are identified (87.5%), with 47 completely and 23 partially matched. SORFs that partially match true internal exons share at least one splice junction with the exon, or share both splice junctions but are interpreted in an incorrect reading frame. Specificity (the percentage of SORFs that correspond to true exons) varies from 91% at confidence level 1 to 16% at confidence level 5, with an overall specificity of 35-40%. The output displays nucleotide position, confidence level, reading frame phase at the 5' and 3' ends, acceptor and donor sequences and scoring statistics and also gives an amino acid translation of the potential exon. SORFIND compares favourably with other programs currently used to predict protein-coding regions.     

Structure and expression of two porcine genomic clones encoding class I MHC antigens

Two nonallelic porcine class I MHC (SLA) genes have been isolated and characterized. Both genes are expressed in mouse L cells, directing the synthesis of class I SLA molecules that carry common monomorphic determinants but are serologically distinct. The corresponding DNA sequences have been determined. The organization of both of these genes is similar to that of other class I genes: a leader exon, three exons encoding extracellular domains, a transmembrane exon, and three intracytoplasmic exons. The two genes are highly homologous in both exon and intron segments, with average homologies of 88% and 80%, respectively. Nucleotide changes in exon 2 are clustered, whereas those in the other exons are dispersed throughout. Comparison of the swine DNA sequences with class I genes from other species reveals a generally high conservation of exons 2, 3, 4, and 6 with lower homology in the remaining protein-encoding domains. Introns are markedly less well conserved, although moderate homology is found between swine and human class I MHC genes in both introns and 3' flanking regions. Taken together with comparisons of the deduced protein sequences, these data indicate an order of swine greater than human greater than rabbit greater than mouse in the relationship of class I genes.     

The human CYP2C locus: a prototype for intergenic and exon repetition splicing events

In human there are four known CYP2C genes that have been mapped to chromosome 10q24 with the order Cen-2C18-2C19-2C9-2C8-Tel. Previously we have shown that splicing events joining exons from the neighboring 2C18 and 2C19 genes occur in human liver and epidermis. Here evidence is presented that the terminal genes of this cluster, 2C18 and 2C8, are also involved in intergenic splicing. Most interestingly, several of these 2C18/2C8 RNAs were composed of all nine exons, thus conceivably having the potential for coding functional proteins. Moreover, chimeric RNA species consisting of exons originating not only from the CYP2C8 and CYP2C18 genes, but also from the CYP2C19 gene were detected. In all cases the exons from the different CYP2C genes were joined at the correct canonical splice sites. However, the closely linked RBP4 gene is not participating in intergenic splicing with the CYP2C genes. In addition, CYP2C8 gene expression was found to generate a variety of scrambled RNA molecules including species that contained repetitions of certain exons.     

Recognizing translation initiation sites of eukaryotic genes based on the cooperatively scanning model

MOTIVATION: Translation initiation sites (TISs) of genes are the key points of protein synthesis. Exact recognition of TISs in eukaryotic genes is one of the most important tasks in gene-finding algorithms. However, the task has not been satisfactorily fulfilled up to the present. Here, we propose a cooperatively scanning model for recognizing TISs and the first exons of eukaryotic genes on the basis of the structural characteristics of multi-exon genes. RESULTS: The model was employed to cooperatively scan the TISs and 3' splicing sites in eukaryotic genes, and the TISs and the first exons of 132 mammalian gene sequences are identified to evaluate the model. Accuracy of exactly recognizing the TISs and the first exons has been found to amount respectively to 64.4 and 51.5%. We believe that the model will be a useful tool for genome annotation and that it can be easily incorporated into other algorithms to achieve higher accuracy in recognizing TISs and the first exons. AVAILABILITY: The program is available upon request.