Translocation events in the evolution of aminoacyl-tRNA synthetases.

We have characterized hisS, the gene encoding the histidyl-tRNA synthetase (HisRS) from the tetraodontoid fish Fugu rubripes. The hisS gene is about 3.5 kbp long and contains 13 exons and 12 introns of 172 bp, on average. The Fugu hisS gene encodes a putative protein of 519 amino acids with the three motifs identified as signatures of class 2 aminoacyl-tRNA synthetases. A model for the shifting of intron 8 between Fugu and hamster is proposed based on the successive appearance of a cryptic splicing site followed by an insertion mutation that created a new acceptor site. In addition, sequence comparisons suggest that the hisS gene has undergone a translocation through the first intron. As a result, the Fugu HisRS has an N-terminal sequence markedly different from that in the human and hamster enzymes. We propose that similar events have been responsible for variations at the N-terminal end of other aminoacyl-tRNA synthetases. Our analysis suggests that this involves exchanges through introns of two exons encoding an ancestral 32-amino acid motif.

Identification of interdomain sequences promoting the intronless evolution of a bacterial protein family.

In the evolution of eukaryotic genes, introns are believed to have played a major role in increasing the probability of favorable duplication events, chance recombinations, and exon shuffling resulting in functional hybrid proteins. As a rule, prokaryotic genes lack introns, and the examples of prokaryotic introns described do not seem to have contributed to gene evolution by exon shuffling. Still, certain protein families in modern bacteria evolve rapidly by recombination of genes, duplication of functional domains, and as shown for protein PAB of the anaerobic bacterial species Peptostreptococcus magnus, by the shuffling of an albumin-binding protein module from group C and G streptococci. Characterization of a protein PAB-related gene in a P. magnus strain with less albumin-binding activity revealed that the shuffled module was missing. Based on this fact and observations made when comparing gene sequences of this family of bacterial surface proteins interacting with albumin and/or immunoglobulin, a model is presented that can explain how this rapid intronless evolution takes place. A new kind of genetic element is introduced: the recer sequence promoting interdomain, in frame recombination and acting as a structure-less flexibility-promoting spacer in the corresponding protein. The data presented also suggest that antibiotics could represent the selective pressure behind the shuffling of protein modules in P. magnus, a member of the indigenous bacterial flora.

Complete structure of the human alpha-albumin gene, a new member of the serum albumin multigene family.

The nucleotide sequence of the human alpha-albumin gene, including 887 bp of the 5'-flanking region and 1311 bp of the 3-flanking region (24,454 in total), was determined from three overlapping lambda phage clones. The sequence spans 22,256 bp from the cap site to the polyadenylylation site, revealing a gene structure of 15 exons separated by 14 introns. The methionine initiation codon ATG is within exon 1; the termination codon TGA is within exon 14. Exon 15 is entirely untranslated and contains the polyadenylylation signal AATAAA. The deduced polypeptide chain is composed of a 21-amino-acid leader peptide, followed by 578 amino acids of the mature protein. There are seven repetitive DNA elements (Alu and Kpn) in the introns and 3-flanking region. The sizes of the 15 alpha-albumin exons match closely those of the albumin, alpha-fetoprotein, and vitamin D-binding protein genes. The exons are symmetrically placed within the three domains of the individual proteins, and they share a characteristic codon splitting pattern that is conserved among members of the gene family. The results provide strong evidence that alpha-albumin belongs to, and most likely completes with, the serum albumin gene family. Based on structural similarity, alpha-albumin appears to be most closely related to alpha-fetoprotein. The complete structure of this family of four tandemly linked genes provides a well-characterized approximately 200 kb locus in the 4q subcentromeric region of the human genome.

Structure of murine and human renal type II Na+-phosphate cotransporter genes (Npt2 and NPT2).

Na+-phosphate (Pi) cotransport across the renal brush border membrane is the rate limiting step in the overall reabsorption of filtered Pi. Murine and human renal-specific cDNAs (NaPi-7 and NaPi-3, respectively) related to this cotransporter activity (type II Na+-Pi cotransporter) have been cloned. We now report the cloning and characterization of the corresponding mouse (Npt2) and human (NPT2) genes. The genes were cloned by screening mouse genomic and human chromosome 5-specific libraries, respectively. Both genes are approximately 16 kb and are comprised of 13 exons and 12 introns, the junctions of which conform to donor and acceptor site consensus sequences. Putative CAAT and TATA boxes are located, respectively, at positions -147 and -40 of the Npt2 gene and -143 and -51 of the NPT2 gene, relative to nucleotide 1 of the corresponding cDNAs. The translation initiation site is within exon 2 of both genes. The first 220 bp of the mouse and human promoter regions exhibit 72% identity. Two transcription start sites (at positions -9 and - 10 with respect to nucleotide 1 of NaPi-7 cDNA) and two polyadenylylation signals were identified in the Npt2 gene by primer extension, 5' and 3' rapid amplification of cDNA ends (RACE). A 484-bp 5' flanking region of the Npt2 gene, comprising the CAAT box, TATA box, and exon 1, was cloned upstream of a luciferase reporter gene; this construct significantly stimulated luciferase gene expression, relative to controls, when transiently transfected into OK cells, a renal cell line expressing type II Na+ -Pi cotransporter activity. The present data provide a basis for detailed analysis of cis and trans elements involved in the regulation of Npt2/NPT2 gene transcription and facilitate screening for mutations in the NPT2 gene in patients with autosomally inherited disorders of renal Pi reabsorption.

RIG-E, a human homolog of the murine Ly-6 family, is induced by retinoic acid during the differentiation of acute promyelocytic leukemia cell.

In vivo all-trans-retinoic acid (ATRA), a differentiation inducer, is capable of causing clinical remission in about 90% of patients with acute promyelocytic leukemia (APL). The molecular basis for the differentiation of APL cells after treatment with ATRA remains obscure and may involve genes other than the known retinoid nuclear transcription factors. We report here the ATRA-induced gene expression in a cell line (NB4) derived from a patient with APL. By differential display-PCR, we isolated and characterized a novel gene (RIG-E) whose expression is up-regulated by ATRA. The gene is 4.0 kb long, consisting of four exons and three introns, and is localized on human chromosome region 8q24. The deduced amino acid sequence predicts a cell surface protein containing 20 amino acids at the N-terminal end corresponding to a signal peptide and an extracellular sequence containing 111 amino acids. The RIG-E coded protein shares some homology with CD59 and with a number of growth factor receptors. It shares high sequence homology with the murine LY-6 multigene family, whose members are small cysteine-rich proteins differentially expressed in several hematopoietic cell lines and appear to function in signal transduction. It seems that so far RIG-E is the closest human homolog of the LY-6 family. Expression of RIG-E is not restricted to myeloid differentiation, because it is also present in thymocytes and in a number of other tissues at different levels.

The molecular basis of Sanfilippo syndrome type B.

The Sanfilippo syndrome type B is a lysosomal storage disorder caused by deficiency of alpha-N-acetylglucosaminidase; it is characterized by profound mental deterioration in childhood and death in the second decade. For understanding the molecular genetics of the disease and for future development of DNA-based therapy, we have cloned the cDNA and gene encoding alpha-N-acetylglucosaminidase. Cloning started with purification of the bovine enzyme and use of a conserved oligonucleotide sequence to probe a human cDNA library. The cDNA sequence was found to encode a protein of 743 amino acids, with a 20- to 23-aa signal peptide immediately preceding the amino terminus of the tissue enzyme and with six potential N-glycosylation sites. The 8.5-kb gene (NAGLU), interrupted by 5 introns, was localized to the 5'-flanking sequence of a known gene, EDH17B, on chromosome 17q21. Five mutations were identified in cells of patients with Sanfilippo syndrome type B: 503del10, R297X, R626X, R643H, and R674H. The occurrence of a frameshift and a nonsense mutation in homozygous form confirms the identity of the NAGLU gene.

Genomic cloning of methylthioadenosine phosphorylase: a purine metabolic enzyme deficient in multiple different cancers.

5'-Deoxy-5'-methylthioadenosine phosphorylase (methylthioadeno-sine: ortho-phosphate methylthioribosyltransferase, EC 24.2.28; MTAP) plays a role in purine and polyamine metabolism and in the regulation of transmethylation reactions. MTAP is abundant in normal cells but is deficient in many cancers. Recently, the genes for the cyclin-dependent kinase inhibitors p16 and p15 have been localized to the short arm of human chromosome 9 at band p21, where MTAP and interferon alpha genes (IFNA) also map. Homozygous deletions of p16 and p15 are frequent malignant cell lines. However, the order of the MTAP, p16, p15, and IFNA genes on chromosome 9p is uncertain, and the molecular basis for MTAP deficiency in cancer is unknown. We have cloned the MTAP gene, and have constructed a topologic map of the 9p21 region using yeast artificial chromosome clones, pulse-field gel electrophoresis, and sequence-tagged-site PCR. The MTAP gene consists of eight exons and seven introns. Of 23 malignant cell lines deficient in MTAP protein, all but one had complete or partial deletions. Partial or total deletions of the MTAP gene were found in primary T-cell acute lymphoblastic leukemias (T-ALL). A deletion breakpoint of partial deletions found in cell lines and primary T-ALL was in intron 4. Starting from the centromeric end, the gene order on chromosome 9p2l is p15, p16, MTAP, IFNA, and interferon beta gene (IFNB). These results indicate that MTAP deficiency in cancer is primarily due to codeletion of the MTAP and p16 genes.

Human plectin: organization of the gene, sequence analysis, and chromosome localization (8q24).

Plectin, a 500-kDa intermediate filament binding protein, has been proposed to provide mechanical strength to cells and tissues by acting as a cross-linking element of the cytoskeleton. To set the basis for future studies on gene regulation, tissue-specific expression, and pathological conditions involving this protein, we have cloned the human plectin gene, determined its coding sequence, and established its genomic organization. The coding sequence contains 32 exons that extend over 32 kb of the human genome. Most of the introns reside within a region encoding the globular N-terminal domain of the molecule, whereas the entire central rod domain and the entire C-terminal globular domain were found to be encoded by single exons of remarkable length, >3 kb and >6 kb, respectively. Overall, the organization of the human plectin gene was strikingly similar to that of human bullous pemphigoid antigen 1 (BPAG1), confirming that both proteins belong to the same gene family. Comparison of the deduced protein sequences for human and rat plectin revealed that they were 93% identical. By using fluorescence in situ hybridization, we have mapped the plectin gene to the long arm of chromosome 8 within the telomeric region. This gene locus (8q24) has previously been implicated in the human blistering skin disease epidermolysis bullosa simplex Ogna. Detailed knowledge of the structure of the plectin gene and its chromosome localization will aid in the elucidation of whether this or any other pathological conditions are linked to alterations in the plectin gene.

Homing events in the gyrA gene of some mycobacteria.

The A subunit of DNA gyrase in Mycobacterium leprae, unlike its counterpart in Mycobacterium tuberculosis, is produced by protein splicing as its gene, gyrA, harbors a 1260-bp in-frame insertion encoding an intein, a putative homing endonuclease. Analysis of the gyrA locus from different mycobacterial species revealed the presence of inteins in Mycobacterium flavescens, Mycobacterium gordonae and Mycobacterium kansasii but not in 10 other pathogenic or saprophytic mycobacteria. In all four cases where intein coding sequences were found, they were localized in the same position in gyrA, immediately downstream of the codon for the key active-site residue Tyr-130. The intein products were similar, but not identical, in sequence and the splice junctions displayed all the features found in other polypeptides known to be produced by protein splicing from a precursor protein. Paired motifs, found in homing endonucleases encoded by some group I RNA introns, and inteins showing endonuclease activity, were present in the gyrA inteins as were other intein-specific signatures. Some strains of M. flavescens, M. gordonae, and M. kansasii were shown by PCR analysis to have inteinless gyrA genes, in contrast to the situation in M. leprae where all the isolates possessed insertions in gyrA. Sequencing of the corresponding regions revealed that, although the GyrA protein sequence was conserved, the nucleotide sequences differed in gyrA genes with and without inteins, suggesting that the homing endonuclease displays sequence specificity.

Pre-mRNA splicing in plants: characterization of Ser/Arg splicing factors.

The fact that animal introns are not spliced out in plants suggests that recognition of pre-mRNA splice sites differs between the two kingdoms. In plants, little is known about proteins required for splicing, as no plant in vitro splicing system is available. Several essential splicing factors from animals, such as SF2/ASF and SC-35, belong to a family of highly conserved proteins consisting of one or two RNA binding domain(s) (RRM) and a C-terminal Ser/Arg-rich (SR or RS) domain. These animal SR proteins are required for splice site recognition and spliceosome assembly. We have screened for similar proteins in plants by using monoclonal antibodies specific for a phosphoserine epitope of the SR proteins (mAb1O4) or for SF2/ASF. These experiments demonstrate that plants do possess SR proteins, including SF2/ASF-like proteins. Similar to the animal SR proteins, this group of proteins can be isolated by two salt precipitations. However, compared to the animal SR proteins, which are highly conserved in size and number, SR proteins from Arabidopsis, carrot, and tobacco exhibit a complex pattern of intra- and interspecific variants. These plant SR proteins are able to complement inactive HeLa cell cytoplasmic S1OO extracts that are deficient in SR proteins, yielding functional splicing extracts. In addition, plant SR proteins were active in a heterologous alternative splicing assay. Thus, these plant SR proteins are authentic plant splicing factors.

Structure of a gene encoding a cytosolic acetyl-CoA carboxylase of hexaploid wheat.

An entire gene encoding wheat (var. Hard Red Winter Tam 107) acetyl-CoA carboxylase [ACCase; acetyl-CoA:carbon-dioxide ligase (ADP-forming), EC 6.4.1.2] has been cloned and sequenced. Comparison of the 12-kb genomic sequence with the 7.4-kb cDNA sequence reported previously revealed 29 introns. Within the coding region, the exon sequence is 98% identical to the known wheat cDNA sequence. A second ACCase gene was identified by sequencing fragments of genomic clones that include the first two exons and the first intron. Additional transcripts were detected by 5' and 3' RACE analysis (rapid amplification of cDNA ends). One set of transcripts had a 5' end sequence identical to the cDNA found previously and another set was identical to the gene reported here. The 3' RACE clones fall into four distinguishable sequence sets, bringing the number of ACCase sequences to six. None of these cDNA or genomic clones encodes a chloroplast targeting signal. Identification of six different sequences suggests that either the cytosolic ACCase genes are duplicated in the three chromosome sets in hexaploid wheat or that each of the six alleles of the cytosolic ACCase gene has a readily distinguishable DNA sequence.

The chicken beta 2-microglobulin gene is located on a non-major histocompatibility complex microchromosome: a small, G+C-rich gene with X and Y boxes in the promoter.

beta 2-Microglobulin is an essential subunit of major histocompatibility complex (Mhc) class I molecules, which present antigenic peptides to T lymphocytes. We sequenced a number of cDNAs and two genomic clones corresponding to chicken beta 2-microglobulin. The chicken beta 2-microglobulin gene has a similar genomic organization but smaller introns and higher G+C content than mammalian beta 2-microglobulin genes. The promoter region is particularly G+C-rich and contains, in addition to interferon regulatory elements, potential S/W, X, and Y boxes that were originally described for mammalian class II but not class I alpha or beta 2-microglobulin genes. There is a single chicken beta 2-microglobulin gene that has little polymorphism in the coding region. Restriction fragment length polymorphisms from Mhc homozygous lines, Mhc congenic lines, and backcross families, as well as in situ hybridization, show that the beta 2-microglobulin gene is located on a microchromosome different from the one that contains the chicken Mhc. We propose that the structural similarities between the beta 2-microglobulin and Mhc genes in the chicken are due to their presence on microchromosomes and suggest that these features and the microchromosomes appeared by deletion of DNA in the lineage leading to the birds.

Intron phase correlations and the evolution of the intron/exon structure of genes.

Two issues in the evolution of the intron/exon structure of genes are the role of exon shuffling and the origin of introns. Using a large data base of eukaryotic intron-containing genes, we have found that there are correlations between intron phases leading to an excess of symmetric exons and symmetric exon sets. We interpret these excesses as manifestations of exon shuffling and make a conservative estimate that at least 19% of the exons in the data base were involved in exon shuffling, suggesting an important role for exon shuffling in evolution. Furthermore, these excesses of symmetric exons appear also in those regions of eukaryotic genes that are homologous to prokaryotic genes: the ancient conserved regions. This last fact cannot be explained in terms of the insertional theory of introns but rather supports the concept that some of the introns were ancient, the exon theory of genes.

Maintenance of pre-mRNA secondary structure by epistatic selection.

Linkage disequilibrium between polymorphisms in a natural population may result from various evolutionary forces, including random genetic drift due to sampling of gametes during reproduction, restricted migration between subpopulations in a subdivided population, or epistatic selection. In this report, we present evidence that the majority of significant linkage disequilibria observed in introns of the alcohol dehydrogenase locus (Adh) of Drosophila pseudoobscura are due to epistatic selection maintaining secondary structure of precursor mRNA (pre-mRNA). Based on phylogenetic-comparative analysis and a likelihood approach, we propose secondary structure models of Adh pre-mRNA for the regions of the adult intron and intron 2 where clustering of linkage disequilibria has been observed. Furthermore, we applied the likelihood ratio test to the phylogenetically predicted secondary structure in intron 1. In contrast to the other two structures, polymorphisms associated with the more conserved stem-loop structure of intron 1 are in low frequency, and linkage disequilibria have not been observed. These findings are qualitatively consistent with a model of compensatory fitness interactions. This model assumes that mutations disrupting pairing in a secondary structural element are individually deleterious if they destabilize a functionally important structure; a second "compensatory" mutation, however, may restabilize the structure and restore fitness.

A nuclear gene of eubacterial origin in Euglena gracilis reflects cryptic endosymbioses during protist evolution.

Genes for glycolytic and Calvin-cycle glyceraldehyde-3-phosphate dehydrogenase (GAPDH) of higher eukaryotes derive from ancient gene duplications which occurred in eubacterial genomes; both were transferred to the nucleus during the course of endosymbiosis. We have cloned cDNAs encoding chloroplast and cytosolic GAPDH from the early-branching photosynthetic protist Euglena gracilis and have determined the structure of its nuclear gene for cytosolic GAPDH. The gene contains four introns which possess unusual secondary structures, do not obey the GT-AG rule, and are flanked by 2- to 3-bp direct repeats. A gene phylogeny for these sequences in the context of eubacterial homologues indicates that euglenozoa, like higher eukaryotes, have obtained their GAPDH genes from eubacteria via endosymbiotic (organelle-to-nucleus) gene transfer. The data further suggest that the early-branching protists Giardia lamblia and Entamoeba histolytica--which lack mitochondria--and portions of the trypanosome lineage have acquired GAPDH genes from eubacterial donors which did not ultimately give rise to contemporary membrane-bound organelles. Evidence that "cryptic" (possibly ephemeral) endosymbioses during evolution may have entailed successful gene transfer is preserved in protist nuclear gene sequences.