AsMamDB: an alternative splice database of mammals

The objective of database AsMamDB is to facilitate the systematic study of alternatively spliced genes of mammals. Version 1.0 of AsMamDB contains 1563 alternatively spliced genes of human, mouse and rat, each associated with a cluster of nucleotide sequences. The main information provided by AsMamDB includes gene alternative splicing patterns, gene structures, locations in chromosomes, products of genes and tissues where they express. Alternative splicing patterns are represented by multiple alignments of various gene transcripts and by graphs of their topological structures. Gene structures are illustrated by exon, intron and various regulatory elements distributions. There are 4204 DNAs, 3977 mRNAs, 8989 CDSs and 126 931 ESTs in the current database. More than 130 000 GenBank entries are covered and 4443 MEDLINE records are linked. DNA, mRNA, exon, intron and relevant regulatory element sequences are provided in FASTA format. More information can be obtained by using the web-based multiple alignment tool Asalign and various category lists. AsMamDB can be accessed at http://166.111.30.6 5/ASMAM DB.html.

SpliceDB: database of canonical and non-canonical mammalian splice sites

A database (SpliceDB) of known mammalian splice site sequences has been developed. We extracted 43 337 splice pairs from mammalian divisions of the gene-centered Infogene database, including sites from incomplete or alternatively spliced genes. Known EST sequences supported 22 815 of them. After discarding sequences with putative errors and ambiguous location of splice junctions the verified dataset includes 22 489 entries. Of these, 98.71% contain canonical GT–AG junctions (22 199 entries) and 0.56% have non-canonical GC–AG splice site pairs. The remainder (0.73%) occurs in a lot of small groups (with a maximum size of 0.05%). We especially studied non-canonical splice sites, which comprise 3.73% of GenBank annotated splice pairs. EST alignments allowed us to verify only the exonic part of splice sites. To check the conservative dinucleotides we compared sequences of human non-canonical splice sites with sequences from the high throughput genome sequencing project (HTG). Out of 171 human non-canonical and EST-supported splice pairs, 156 (91.23%) had a clear match in the human HTG. They can be classified after sequence analysis as: 79 GC–AG pairs (of which one was an error that corrected to GC–AG), 61 errors corrected to GT–AG canonical pairs, six AT–AC pairs (of which two were errors corrected to AT–AC), one case was produced from a non-existent intron, seven cases were found in HTG that were deposited to GenBank and finally there were only two other cases left of supported non-canonical splice pairs. The information about verified splice site sequences for canonical and non-canonical sites is presented in SpliceDB with the supporting evidence. We also built weight matrices for the major splice groups, which can be incorporated into gene prediction programs. SpliceDB is available at the computational genomic Web server of the Sanger Centre: http://genomic.sanger.ac.uk/spldb/SpliceDB.html and at http://www.softberry.com/spldb/SpliceDB.html.

Point mutant mice with hypersensitive α4 nicotinic receptors show dopaminergic deficits and increased anxiety

Knock-in mice were generated that harbored a leucine-to-serine mutation in the α4 nicotinic receptor near the gate in the channel pore. Mice with intact expression of this hypersensitive receptor display dominant neonatal lethality. These mice have a severe deficit of dopaminergic neurons in the substantia nigra, possibly because the hypersensitive receptors are continuously activated by normal extracellular choline concentrations. A strain that retains the neo selection cassette in an intron has reduced expression of the hypersensitive receptor and is viable and fertile. The viable mice display increased anxiety, poor motor learning, excessive ambulation that is eliminated by very low levels of nicotine, and a reduction of nigrostriatal dopaminergic function upon aging. These knock-in mice provide useful insights into the pathophysiology of sustained nicotinic receptor activation and may provide a model for Parkinson's disease.

A method for the generation of conditional gene repair mutations in mice

Conditional gene repair mutations in the mouse can assist in cell lineage analyses and provide a valuable complement to conditional gene inactivation strategies. We present a method for the generation of conditional gene repair mutations that employs a loxP-flanked (floxed) selectable marker and transcriptional/translational stop cassette (neostop) located within the first intron of a target gene. In the absence of Cre recombinase, expression of the targeted allele is suppressed generating a null allele, while in the presence of Cre, excision of neostop restores expression to wild-type levels. To test this strategy, we have generated a conditional gene repair allele of the mouse Huntington’s disease gene homolog (Hdh). Insertion of neostop within the Hdh intron 1 generated a null allele and mice homozygous for this allele resembled nullizygous Hdh mutants and died after embryonic day 8.5. In the presence of a cre transgene expressed ubiquitously early in development, excision of neostop restored Hdh expression and rescued the early embryonic lethality. A simple modification of this strategy that permits the generation of conventional gene knockout, conditional gene knockout and conditional gene repair alleles using one targeting construct is discussed.

Formation of a GNRA tetraloop in P5abc can disrupt an interdomain interaction in the Tetrahymena group I ribozyme

The secondary structure of a truncated P5abc subdomain (tP5abc, a 56-nucleotide RNA) of the Tetrahymena thermophila group I intron ribozyme changes when its tertiary structure forms. We have now used heteronuclear NMR spectroscopy to determine its conformation in solution. The tP5abc RNA that contains only secondary structure is extended compared with the tertiary folded form; both forms coexist in slow chemical exchange (the interconversion rate constant is slower than 1 s−1) in the presence of magnesium. Kinetic experiments have shown that tertiary folding of the P5abc subdomain is one of the earliest folding transitions in the group I intron ribozyme, and that it leads to a metastable misfolded intermediate. Previous mutagenesis studies suggest that formation of the extended P5abc structure described here destabilize a misfolded intermediate. This study shows that the P5abc RNA subdomain containing a GNRA tetraloop in P5c (in contrast to the five-nucleotide loop P5c in the tertiary folded ribozyme) can disrupt the base-paired interdomain (P14) interaction between P5c and P2.

A human protein containing multiple types of protease-inhibitory modules

By using sensitive homology-search and gene-finding programs, we have found that a genomic region from the tip of the short arm of human chromosome 16 (16p13.3) encodes a putative secreted protein consisting of a domain related to the whey acidic protein (WAP) domain, a domain homologous with follistatin modules of the Kazal-domain family (FS module), an immunoglobulin-related domain (Ig domain), two tandem domains related to Kunitz-type protease inhibitor modules (KU domains), and a domain belonging to the recently defined NTR-module family (NTR domain). The gene encoding these WAP, FS, Ig, KU, and NTR modules (hereafter referred to as the WFIKKN gene) is intron-depleted—its single 1,157-bp intron splits the WAP module. The validity of our gene prediction was confirmed by sequencing a WFIKKN cDNA cloned from a lung cDNA library. Studies on the tissue-expression pattern of the WFIKKN gene have shown that the gene is expressed primarily in pancreas, kidney, liver, placenta, and lung. As to the function of the WFIKKN protein, it is noteworthy that it contains FS, WAP, and KU modules, i.e., three different module types homologous with domains frequently involved in inhibition of serine proteases. The protein also contains an NTR module, a domain type implicated in inhibition of zinc metalloproteinases of the metzincin family. On the basis of its intriguing homologies, we suggest that the WFIKKN protein is a multivalent protease inhibitor that may control the action of multiple types of serine proteases as well as metalloproteinase(s).

The human XPG gene: gene architecture, alternative splicing and single nucleotide polymorphisms

Defects in the XPG DNA repair endonuclease gene can result in the cancer-prone disorders xeroderma pigmentosum (XP) or the XP–Cockayne syndrome complex. While the XPG cDNA sequence was known, determination of the genomic sequence was required to understand its different functions. In cells from normal donors, we found that the genomic sequence of the human XPG gene spans 30 kb, contains 15 exons that range from 61 to 1074 bp and 14 introns that range from 250 to 5763 bp. Analysis of the splice donor and acceptor sites using an information theory-based approach revealed three splice sites with low information content, which are components of the minor (U12) spliceosome. We identified six alternatively spliced XPG mRNA isoforms in cells from normal donors and from XPG patients: partial deletion of exon 8, partial retention of intron 8, two with alternative exons (in introns 1 and 6) and two that retained complete introns (introns 3 and 9). The amount of alternatively spliced XPG mRNA isoforms varied in different tissues. Most alternative splice donor and acceptor sites had a relatively high information content, but one has the U12 spliceosome sequence. A single nucleotide polymorphism has allele frequencies of 0.74 for 3507G and 0.26 for 3507C in 91 donors. The human XPG gene contains multiple splice sites with low information content in association with multiple alternatively spliced isoforms of XPG mRNA.

Structures of two molluscan hemocyanin genes: Significance for gene evolution

We present here the description of genes coding for molluscan hemocyanins. Two distantly related mollusks, Haliotis tuberculata and Octopus dofleini, were studied. The typical architecture of a molluscan hemocyanin subunit, which is a string of seven or eight globular functional units (FUs, designated a to h, about 50 kDa each), is reflected by the gene organization: a series of eight structurally related coding regions in Haliotis, corresponding to FU-a to FU-h, with seven highly variable linker introns of 174 to 3,198 bp length (all in phase 1). In Octopus seven coding regions (FU-a to FU-g) are found, separated by phase 1 introns varying in length from 100 bp to 910 bp. Both genes exhibit typical signal (export) sequences, and in both cases these are interrupted by an additional intron. Each gene also contains an intron between signal peptide and FU-a and in the 3′ untranslated region. Of special relevance for evolutionary considerations are introns interrupting those regions that encode a discrete functional unit. We found that five of the eight FUs in Haliotis each are encoded by a single exon, whereas FU-f, FU-g, and FU-a are encoded by two, three and four exons, respectively. Similarly, in Octopus four of the FUs each correspond to an uninterrupted exon, whereas FU-b, FU-e, and FU-f each contain a single intron. Although the positioning of the introns between FUs is highly conserved in the two mollusks, the introns within FUs show no relationship either in location nor phase. It is proposed that the introns between FUs were generated as the eight-unit polypeptide evolved from a monomeric precursor, and that the internal introns have been added later. A hypothesis for evolution of the ring-like quaternary structure of molluscan hemocyanins is presented.

Dynamics of mobile element activity in chalcone synthase loci in the common morning glory (Ipomoea purpurea)

Mobile element dynamics in seven alleles of the chalcone synthase D locus (CHS-D) of the common morning glory (Ipomoea purpurea) are analyzed in the context of synonymous nucleotide sequence distances for CHS-D exons. By using a nucleotide sequence of CHS-D from the sister species Ipomoea nil (Japanese morning glory [Johzuka-Hisatomi, Y., Hoshino, A., Mori, T., Habu, Y. & Iida, S. (1999) Genes Genet. Syst. 74, 141–147], it is also possible to determine the relative frequency of insertion and loss of elements within the CHS-D locus between these two species. At least four different types of transposable elements exist upstream of the coding region, or within the single intron of the CHS-D locus in I. purpurea. There are three distinct families of miniature inverted-repeat transposable elements (MITES), and some recent transpositions of Activator/Dissociation (Ac/Ds)-like elements (Tip100), of some short interspersed repetitive elements (SINEs), and of an insertion sequence (InsIpCHSD) found in the neighborhood of this locus. The data provide no compelling evidence of the transposition of the mites since the separation of I. nil and I. purpurea roughly 8 million years ago. Finally, it is shown that the number and frequency of mobile elements are highly heterogeneous among different duplicate CHS loci, suggesting that the dynamics observed at CHS-D are locus-specific.

Origin of Genes

We discuss two tests of the hypothesis that the first genes were assembled from exons. The hypothesis of exon shuffling in the progenote predicts that intron phases will be correlated so that exons will be an integer number of codons and predicts that the exons will be correlated with compact regions of polypeptide chain. These predictions have been tested on ancient conserved proteins (proteins without introns in prokaryotes but with introns in eukaryotes) and hold with high statistical significance. We conclude that introns are correlated with compact features of proteins 15-, 22-, or 30-amino acid residues long, as was predicted by “The Exon Theory of Genes.”

Dynamic evolution of plant mitochondrial genomes: Mobile genes and introns and highly variable mutation rates

We summarize our recent studies showing that angiosperm mitochondrial (mt) genomes have experienced remarkably high rates of gene loss and concomitant transfer to the nucleus and of intron acquisition by horizontal transfer. Moreover, we find substantial lineage-specific variation in rates of these structural mutations and also point mutations. These findings mostly arise from a Southern blot survey of gene and intron distribution in 281 diverse angiosperms. These blots reveal numerous losses of mt ribosomal protein genes but, with one exception, only rare loss of respiratory genes. Some lineages of angiosperms have kept all of their mt ribosomal protein genes whereas others have lost most of them. These many losses appear to reflect remarkably high (and variable) rates of functional transfer of mt ribosomal protein genes to the nucleus in angiosperms. The recent transfer of cox2 to the nucleus in legumes provides both an example of interorganellar gene transfer in action and a starting point for discussion of the roles of mechanistic and selective forces in determining the distribution of genetic labor between organellar and nuclear genomes. Plant mt genomes also acquire sequences by horizontal transfer. A striking example of this is a homing group I intron in the mt cox1 gene. This extraordinarily invasive mobile element has probably been acquired over 1,000 times separately during angiosperm evolution via a recent wave of cross-species horizontal transfers. Finally, whereas all previously examined angiosperm mtDNAs have low rates of synonymous substitutions, mtDNAs of two distantly related angiosperms have highly accelerated substitution rates.

Ullrich scleroatonic muscular dystrophy is caused by recessive mutations in collagen type VI

Ullrich syndrome is a recessive congenital muscular dystrophy affecting connective tissue and muscle. The molecular basis is unknown. Reverse transcription–PCR amplification performed on RNA extracted from fibroblasts or muscle of three Ullrich patients followed by heteroduplex analysis displayed heteroduplexes in one of the three genes coding for collagen type VI (COL6). In patient A, we detected a homozygous insertion of a C leading to a premature termination codon in the triple-helical domain of COL6A2 mRNA. Both healthy consanguineous parents were carriers. In patient B, we found a deletion of 28 nucleotides because of an A → G substitution at nucleotide −2 of intron 17 causing the activation of a cryptic acceptor site inside exon 18. The second mutation was an exon skipping because of a G → A substitution at nucleotide −1 of intron 23. Both mutations are present in an affected brother. The first mutation is also present in the healthy mother, whereas the second mutation is carried by their healthy father. In patient C, we found only one mutation so far—the same deletion of 28 nucleotides found in patient B. In this case, it was a de novo mutation, as it is absent in her parents. mRNA and protein analysis of patient B showed very low amounts of COL6A2 mRNA and of COL6. A near total absence of COL6 was demonstrated by immunofluorescence in fibroblasts and muscle. Our results demonstrate that Ullrich syndrome is caused by recessive mutations leading to a severe reduction of COL6.

The Two Genes Encoding Starch-Branching Enzymes IIa and IIb Are Differentially Expressed in Barley1

The sbeIIa and sbeIIb genes, encoding starch-branching enzyme (SBE) IIa and SBEIIb in barley (Hordeum vulgare L.), have been isolated. The 5′ portions of the two genes are strongly divergent, primarily due to the 2064-nucleotide-long intron 2 in sbeIIb. The sequence of this intron shows that it contains a retro-transposon-like element. Expression of sbeIIb but not sbeIIa was found to be endosperm specific. The temporal expression patterns for sbeIIa and sbeIIb were similar and peaked around 12 d after pollination. DNA gel-blot analysis demonstrated that sbeIIa and sbeIIb are both single-copy genes in the barley genome. By fluorescence in situ hybridization, the sbeIIa and sbeIIb genes were mapped to chromosomes 2 and 5, respectively. The cDNA clones for SBEIIa and SBEIIb were isolated and sequenced. The amino acid sequences of SBEIIa and SBEIIb were almost 80% identical. The major structural difference between the two enzymes was the presence of a 94-amino acid N-terminal extension in the SBEIIb precursor. The (β/α)8-barrel topology of the α-amylase superfamily and the catalytic residues implicated in branching enzymes are conserved in both barley enzymes.

Allele-Dependent Barley Grain β-Amylase Activity1

The wild ancestor of cultivated barley, Hordeum vulgare subsp. spontaneum (K. Koch) A. & Gr. (H. spontaneum), is a source of wide genetic diversity, including traits that are important for malting quality. A high β-amylase trait was previously identified in H. spontaneum strains from Israel, and transferred into the backcross progeny of a cross with the domesticated barley cv Adorra. We have used Southern-blot analysis and β-amy1 gene characterization to demonstrate that the high β-amylase trait in the backcross line is co-inherited with the β-amy1 gene from the H. spontaneum parent. We have analyzed the β-amy1 gene organization in various domesticated and wild-type barley strains and identified three distinct β-amy1 alleles. Two of these β-amy1 alleles were present in modern barley, one of which was specifically found in good malting barley cultivars. The third allele, linked with high grain β-amylase activity, was found only in a H. spontaneum strain from the Judean foothills in Israel. The sequences of three isolated β-amy1 alleles are compared. The involvement of specific intron III sequences, in particular a 126-bp palindromic insertion, in the allele-dependent expression of β-amylase activity in barley grain is proposed.

Three Drought-Responsive Members of the Nonspecific Lipid-Transfer Protein Gene Family in Lycopersicon pennellii Show Different Developmental Patterns of Expression1

Genomic clones of two nonspecific lipid-transfer protein genes from a drought-tolerant wild species of tomato (Lycopersicon pennellii Corr.) were isolated using as a probe a drought- and abscisic acid (ABA)-induced cDNA clone (pLE16) from cultivated tomato (Lycopersicon esculentum Mill.). Both genes (LpLtp1 and LpLtp2) were sequenced and their corresponding mRNAs were characterized; they are both interrupted by a single intron at identical positions and predict basic proteins of 114 amino acid residues. Genomic Southern data indicated that these genes are members of a small gene family in Lycopersicon spp. The 3′-untranslated regions from LpLtp1 and LpLtp2, as well as a polymerase chain reaction-amplified 3′-untranslated region from pLE16 (cross-hybridizing to a third gene in L. pennellii, namely LpLtp3), were used as gene-specific probes to describe expression in L. pennellii through northern-blot analyses. All LpLtp genes were exclusively expressed in the aerial tissues of the plant and all were drought and ABA inducible. Each gene had a different pattern of expression in fruit, and LpLtp1 and LpLtp2, unlike LpLtp3, were both primarily developmentally regulated in leaf tissue. Putative ABA-responsive elements were found in the proximal promoter regions of LpLtp1 and LpLtp2.