The structure of the human liver-type phosphofructokinase gene

We have isolated the gene for the human liver-type phosphofructokinase, from upstream to the 5' mRNA terminus to beyond the polyadenylation site. The gene is at least 28 kb long and is divided into 22 exons; it contains conventional splice-junction sequences and one polyadenylation signal. Exons and introns are quite rich in G and C residues; some 60% of all nucleotides are either G or C. Five possible sites of polymorphism have been found. The gene structure reveals no signs of internal similarities despite protein sequence evidence which suggests that the PFK molecule is divided into two similar halves. The structure and organization of the human liver-type PFK gene are shown to be extremely similar to those of the rabbit muscle-type PFK.     

Isolation and sequence analysis of a β-tubulin gene from arbuscular mycorrhizal fungi

A full-length β-tubulin gene has been cloned and sequenced from Gigaspora gigantea and Glomus clarum, two arbuscular mycorrhizal fungi (AMF) species in the phylum Glomeromyota. The gene in both species is organized into five exons and four introns. Both genes are 94.9% similar and encode a 447 amino acid protein. In comparison with other fungal groups, the amino acid sequence is most similar to that of fungi in the Chytridiomycota. The codon usage of the gene in both AMF species is broad and biased in favor of an A or a T in the third position. The four introns varied in length from 87 to 168 bp for G. gigantea and from 90 to 136 bp for G. clarum. Of all fungi in which full-length sequences have been published, only AMF do not have an intron before codon 174. The introns positioned at codons 174 and 257 in AMF match the position of different introns in β-tubulin genes of some Zygomycete, Basidiomycete, and Ascomycete fungi. The 5′ and 3′ splice site consensus sequences are similar to those found in introns of most fungi. Sequence analysis from single-strand conformation polymorphism analysis confirmed the presence of two β-tubulin gene copies in G. clarum, but only one copy was evident in G. gigantea based on Southern hybridization analysis.     

Structure of the human hemopexin gene and evidence for intron-mediated evolution

Summary The human hemopexin gene was isolated and its structure determined. The gene spans approximately 12 kb and is interrupted by nine introns. When the intron/exon pattern was examined with respect to the polypeptide segments they encode, a direct correspondence between exons and the 10 repeating units in the protein was observed. The introns are not randomly placed; they fall in the middle of the region of amino acid sequence homology in strikingly similar locations in 6 of the 10 units and in a symmetrical position in the two halves of the coding sequence. These features strongly support the hypothesis that the gene evolved through intron-mediated duplications of a primordial sequence to a five-exon cluster. A more recent gene duplication led to the present-day gene organization.     

Isolation, sequence analysis, and intron-exon arrangement of the gene encoding bovine rhodopsin

We have isolated cDNA clones generated from the mRNA encoding the opsin apoprotein of bovine rhodopsin and used these cDNAs to isolate genomic DNA clones containing the complete opsin gene. Nucleotide sequence analysis of the cloned DNAs has yielded a complete amino acid sequence for bovine rhodopsin and provided an intron-exon map of its gene. The mRNA homologous sequences in the 6.4 kb gene consist of a 96 bp 5' untranslated region, a 1044 bp coding region, and a surprisingly long approximately 1400 bp 3' untranslated region, and are divided into five exons by four introns that interrupt the coding region. Secondary structure analysis predicts that the bovine rhodopsin chain, like that of bacteriorhodopsin, contains seven transmembrane segments. Interestingly, three of the four introns are immediately distal to the codons for three of these segments, and one of these introns marks the boundary between the C-terminal domain and a transmembrane domain.     

Immunological identification of rat tissue kallikrein cDNA and characterization of the kallikrein gene family

A tissue kallikrein cDNA was identified by direct immunological screening with affinity-purified anti-rat tissue kallikrein antibody from a rat submandibular cDNA library constructed with the expression vector pUC8. Sequence analysis of the kallikrein cDNA revealed an encoded protein 97% homologous to the partial amino acid sequence of rat submandibular kallikrein. This cDNA was used to hybrid-select kallikrein-specific RNA from submandibular gland. Translation of the hybrid-selected RNA in a cell-free assay system resulted in the production of a 37 kDa peptide representing the preproenzyme. In addition, hybrid-selection of RNA under less stringent conditions showed cross-hybridization with other submandibular gland mRNA species. In correlation with these results, analysis of rat genomic DNA showed extensive hybridization, suggesting a family of closely related kallikrein-like genes. Consequently, a Charon 4A rat genomic library was screened for kallikrein genes by hybridization with rat tissue kallikrein cDNA. Thirty-four clones were isolated and found to be highly homologous by hybridization and restriction enzymes analyses. Fourteen unique clones were identified by restriction enzyme site polymorphisms within DNA segments which hybridized to the kallikrein cDNA probe and it was estimated that at least 17 different kallikrein-like genes are present in the rat. Sequence and structural analysis of one of the genomic clones revealed a gene structure similar to that of other serine proteinases. Comparison of the partially sequenced exon regions of the gene with the sequence of rat tissue kallikrein cDNA reveals 89% identity when aligned for the greatest homology. However, the genomic sequence predicts termination codons in all three translational reading frames, implying that this gene is nonfunctional, i.e., a pseudogene. Comparison of the rat genomic sequence to a kallikrein-like gene from the mouse reveals extensive preservation of exons, less identity within introns and no significant homology between extragenic regions.     

Isolation and characterization of the chicken cardiac myosin light chain (L-2A) gene. Evidence for two additional N-terminal amino acids

The contractile proteins of striated muscle are encoded by multigene families and constitute an excellent system to investigate differentiation and developmental control of gene expression. Different forms of myosin light chains are expressed in skeletal muscle as well as in the myocard. To study the gene structure and molecular mechanisms underlying differential gene expression, the structural cardiac myosin light chain 2 (MLC-2A) gene was isolated from a chicken genomic DNA library. Restriction enzyme mapping, electron microscopic analysis, and partial sequencing revealed that the gene coding for the MLC mRNA of 700 nucleotides in length extends over 4.2 kilobases of DNA and is interrupted by 5 introns. Sequence analysis led to the detection of two codons for additional amino acids at the N terminus which were not reported to be present in the mature protein and are presumably removed post-translationally. These two amino acids, methionine and alanine, are coded on two separate exons split by the largest intron of the entire gene. Southern blot analysis of genomic chicken DNA indicates the presence of one MLC-2A gene per haploid chicken genome.     

Genomic organization of the human multidrug resistance (MDR1) gene and origin of P-glycoproteins

The MDR1 gene, responsible for multidrug resistance in human cells, encodes a broad specificity efflux pump (P-glycoprotein). P-glycoprotein consists of two similar halves, each half including a hydrophobic transmembrane region and a nucleotide-binding domain. On the basis of sequence homology between the N-terminal and C-terminal halves of P-glycoprotein, we have previously suggested that this gene arose by duplication of a primordial gene. We have now determined the complete intron/exon structure of the MDR1 gene by direct sequencing of cosmid clones and enzymatic amplification of genomic DNA segments. The MDR1 gene includes 28 introns, 26 of which interrupt the protein-coding sequence. Although both halves of the protein-coding sequence are composed of approximately the same number of exons, only two intron pairs, both within the nucleotide-binding domains, are located at conserved positions in the two halves of the protein. The other introns occur at different locations in the two halves of the protein and in most cases interrupt the coding sequence at different positions relative to the open reading frame. These results suggest that the P-glycoprotein arose by fusion of genes for two related but independently evolved proteins rather than by internal duplication.     

Partial structure of the gene for chicken cartilage proteoglycan core protein

A genomic DNA fragment (gCORE-1), encoding a portion of the cartilage proteoglycan core protein, has been isolated from a phage library using cDNA as a probe. The genomic insert is about 17 kilobase pairs; two BamHI fragments of the insert (1.3 and 4.8 kilobase pairs) contain most of the hybridizable sequences found in the cDNA. Sequence analysis of these fragments shows that they contain a total of five exons that encompass 216 amino acid residues, all of which are identical to those of the corresponding cDNA sequence. Three of the exons, which are adjacent to one another, are very similar to the corresponding exons in the gene of a rat hepatic lectin as well as to an exon in the gene of human pulmonary surfactant-associated protein. There is a strong degree of conservation of amino acid sequences encoded in the three genes, although there is no similarity between their introns. The sizes of the five exons in gCORE-1, except for one (which is indeterminate because only a partial cDNA sequence is available), are less than 184 base pairs, whereas the sizes of the introns range from 218 to greater than 2629 base pairs. Four of the introns interrupt an exon codon at either their donor or acceptor sites, between the first and second nucleotides. Only one intron does not split a codon. Intron and exon boundary sites are in agreement with known consensus sequences for introns. The dispersed distribution and relatively small size of the exons, if representative of the entire gene, suggest that the complete gene which codes for the core protein may be quite sizable.     

Module-intron correlation and intron sliding in family F/10 xylanase genes

Xylanases are classified into two families, numbered F/10 and G/11 according to the similarity of amino acid sequences of their catalytic domain (Henrissat, B., Bairoch, A., 1993. New families in the classification of glycosyl hydrolases based on amino acid sequence similarities. Biochem. J. 293, 781-788). Three-dimensional structure of the catalytic domain of the family F/10 xylanase was reported (White, A., Withers, S.G., Gilkes, N.R., Rose, D.R., 1994. Crystal structure of the catalytic domain of the beta-1,4-glycanase Cex from Cellulomonas fimi. Biochemistry 33, 12546-12552). The domain was decomposed into 22 modules by centripetal profiles (Go, M., Nosaka, M., 1987. Protein architecture and the origin of introns. Cold Spring Harbor Symp. Quant. Biol. 52, 915-924; Noguti, T., Sakakibara, H., Go, M., 1993. Localization of hydrogen-bonds within modules in barnase. Proteins 16, 357-363). A module is a contiguous polypeptide segment of amino acid residues having a compact conformation within a globular domain. Collected 31 intron sites of the family F/10 xylanase genes from fungus were found to be correlated to module boundaries with considerable statistical force (p values <0.001). The relationship between the intron locations and protein structures provides supporting evidence for the ancient origin of introns, because such a relationship cannot be expected by random insertion of introns into eukaryotic genes, but it rather suggests pre-existence of introns in the ancestral genes of prokaryotes and eukaryotes. A phylogenetic tree of the fungal and bacterial xylanase sequences made two clusters; one includes both the bacterial and fungal genes, but the other consists of only fungal genes. The mixed cluster of bacterial genes without introns and the fungal genes with introns further supports the ancient origin of introns. Comparison of the conserved base sequences of introns indicates that sliding of a splice site occurred in Aspergillus kawachii gene by one base from the ancestral position. Substrate-binding sites of xylanase are localized on eight modules, and introns are found at both termini of six out of these functional modules. This result suggests that introns might play a functional role in shuffling the exons encoding the substrate-binding modules.     

Statistical analysis of the exon-intron structure of higher and lower eukaryote genes

Statistics of the exon-intron structure and splicing sites of several diverse eukaryotes was studied. The yeast exon-intron structures have a number of unique features. A yeast gene usually have at most one intron. The branch site is strongly conserved, whereas the polypirimidine tract is short. Long yeast introns tend to have stronger acceptor sites. In other species the branch site is less conserved and often cannot be determined. In non-yeast samples there is an almost universal correlation between lengths of neighboring exons (all samples excluding protists) and correlation between lengths of neighboring introns (human, drosophila, protists). On the average first introns are longer, and anomalously long introns are usually first introns in a gene. There is a universal preference for exons and exon pairs with the (total) length divisible by 3. Introns positioned between codons are preferred, whereas those positioned between the first and second positions in codon are avoided. The choice of A or G at the third position of intron (the donor splice sites generally prefer purines at this position) is correlated with the overall GC-composition of the gene. In all samples dinucleotide AG is avoided in the region preceding the acceptor site.     

Isolation, characterization, cDNA cloning and gene expression of an avian transthyretin. Implications for the evolution of structure and function of transthyretin in vertebrates.

A chicken liver cDNA library was constructed in bacteriophage lambda gt10. A full-length transthyretin cDNA clone was identified by screening with rat transthyretin cDNA and was sequenced. A three-dimensional model of chicken transthyretin was obtained by computer-graphics-based prediction from the derived amino acid sequence for chicken transthyretin and from the structure of human transthyretin determined by X-ray diffraction analysis [Blake, C.C.F., Geisow, M.J., Oatley, S.J., Rérat, B. & Rérat, C. (1978) J. Mol. Biol. 121, 339-356]. The similarity of the amino acid sequences of chicken and human transthyretins was 75% overall and 100% for the central channel containing the thyroxine-binding site. Also, the organization of the transthyretin gene into exons and introns and the tissue specificity of expression of the transthyretin gene were similar in chicken and mammals, despite an evolutionary distance of about 3 x 10(8) years from their common ancestor, the Cotylosaurus. By far the highest levels of transthyretin mRNA were found in choroid plexus. The data suggest a fundamental role for the cerebral expression of transthyretin in all vertebrates. It has been proposed that this role is the transport of thyroxine from the bloodstream to the brain [Schreiber, G., Aldred, A.R., Jaworowski, A., Nilsson, C., Achen, M.G. & Segal, M.B. (1990) Am. J. Physiol. 258, R338-R345].     

Sequence of the chicken delta 2 crystallin gene and its intergenic spacer. Extreme homology with the delta 1 crystallin gene

The chicken delta-crystallin locus consists of 2 nonallelic, tandemly arranged genes (5'-delta 1-delta 2-3'). Only the delta 1 gene is known to be expressed. The nucleotide sequence for the delta 1 gene has been reported recently (Nickerson, J.M., Wawrousek, E.F., Hawkins, J.W., Wakil, A.S., Wistow, G.J., Thomas, G., Norman, B.L., and Piatigorsky, J. (1985) J. Biol. Chem. 260, 9100-9105; Ohno, M., Sakamoto, H., Yasuda, K., Okada, T.S., and Shimura, Y. (1985) Nucleic Acids Res. 13, 1593-1606). We now report the sequence for the delta 2 gene and the 4-kilobase intergenic spacer between the two delta-crystallin genes. The delta 2 gene, like the delta 1 gene, has 17 exons and 16 introns. The homologous exons are remarkably similar: exons 3-17 are identical in size between delta 1 and delta 2, and the sequence homology ranges from 70% (exon 2) to 100% (exons 7, 12, and 15), with the remaining exons having 89-98% identity between the delta 1 and delta 2 genes. Consequently, the encoded delta 2 polypeptide is 91% identical to the delta 1 polypeptide. Considerable similarity also exists between homologous introns of delta 1 and delta 2, with most of the differences accounted for by insertions and/or deletions. The presence of a TATA box, consensus splice junctions (almost identical to the delta 1 gene), lariat branch sequences, and a polyadenylation signal strengthen the possibility that delta 2 is a functional gene.     

A simple program to calculate codon bias index

A computer program (PCBI) was developed to quickly calculate codon bias index (CBI). PCBI can analyze a gene containing introns. The 22 preferred codons defined fromSaccharomyces cerevisiae were used in PCBI as the standard to measure the CBI values. However, users can modify the preferred codons to suit each organism. The data PCBI provides include DNA sequence of open reading frame without introns, amino acid sequence of gene product, a table of amino acid composition, a table of codon usage and (G+C) content, parameters for calculating CBI, and the value of CBI. PCBI runs on DOS or Windows environment, but results can be saved in ASCII text format.     

Nucleotide sequence and intron structure of the apocytochrome b gene of Neurospora crassa mitochondria

The sequence of the apocytochrome b (cob) gene of Neurospora crassa has been determined. The structural gene is interrupted by two intervening sequences of approximately 1260 bp each. The polypeptide encoded by the exons shows extensive homology with the cob proteins of Aspergillus nidulans and Saccharomyces cerevisiae (79% and 60%, respectively). The two introns are, however, located at sites different from those of introns in the cob genes of A. nidulans and S. cerevisiae (which contain highly homologous introns at the same site within the gene). The introns share several short regions of sequence homology (10-12 bp long) with each other and with other fungal mitochondrial introns. Moreover, the second intron contains a 50 nucleotide long sequence that is highly homologous with sequences within every ribosomal intron of fungal mitochondria sequenced to date. The conserved sequences may allow the formation of a core secondary structure, which is nearly identical in many mitochondrial introns. The conserved secondary structure may be required for intron splicing. The second intron contains an open reading frame, continuous with the preceding exon, of approximately 290 codons. Two stretches of 10 amino acid residues, conserved in many introns, are present in the open reading frame.